Structure determination of the Antp (C39----S) homeodomain from nuclear magnetic resonance data in solution using a novel strategy for the structure calculation with the programs DIANA, CALIBA, HABAS and GLOMSA.

The structure of a mutant Antennapedia homeodomain, Antp(C39----S), from Drosophila melanogaster was determined using a set of new programs introduced in the accompanying paper. An input dataset of 957 distance constraints and 171 dihedral angle constraints was collected using two-dimensional n.m.r. experiments with the 15N-labeled protein. The resulting high quality structure for Antp(C39----S), with an average root-mean-square deviation of 0.53 A between the backbone atoms of residues 7 to 59 in 20 energy-refined distance geometry structures and the mean structure, is nearly identical to the previously reported structure of the wild-type Antp homeodomain. The only significant difference is in the connection between helices III and IV, which was found to be less kinked than was indicated by the structure determination for Antp. The main emphasis of the presentation in this paper is on a detailed account of the practical use of a novel strategy for the computation of nuclear magnetic resonance structures of proteins with the combined use of the programs DIANA, CALIBA, HABAS and GLOMSA.     

Refined solution structure and ligand-binding properties of PDC-109 domain b. A collagen-binding type II domain.

We have determined, via 1H-n.m.r., the solution conformation of the collagen-binding b-domain of the bovine seminal fluid protein PDC-109 (PDC-109/b). The structure determination is based on 341 interproton distance estimates and 42 dihedral angle estimates: a set of 24 initial structures were computed; 12 using the variable target function program DIANA, and 12 using the metric matrix program DISGEO. These structures were optimized by restrained energy minimization and dynamic simulated annealing using the CHARMM and X-PLOR programs. The average pairwise root-mean-square difference (r.m.s.d) between the optimized DIANA (DISGEO) structures is 0.71 A (0.82 A) for the backbone atoms, and 1.73 A (2.03 A) for all atoms. Both sets of structures exhibit the same global fold, secondary structure and placement of most non-polar side-chains. Two central antiparallel beta-sheets, which lie roughly perpendicular to each other, and two irregular loops support a large, partially exposed, hydrophobic surface that defines a putative binding site. A test of a hybrid relaxation matrix-based distance refinement protocol (MIDGE program) was performed using a normalized 250 millisecond NOESY spectrum. The resulting distances were input to the molecular mechanics/dynamics procedures mentioned above in order to optimize the DIANA structures. Our results indicate that relaxation matrix refinement of distances is most useful when used conservatively for identifying underestimated distance constraints. 1H-n.m.r. monitored ligand titration experiments revealed definite, albeit weak, binding interactions for phenethylamine and leucine analogs (Ka less than or equal to 25 M-1). Residues perturbed by ligand binding include Tyr7, Trp26, Tyr33, Asp34 and Trp39. These results suggest that PDC-109/b may recognize specific leucine and/or isoleucine-containing sequences within collagen.     

Improved efficiency of protein structure calculations from NMR data using the program DIANA with redundant dihedral angle constraints

Summary A new strategy for NMR structure calculations of proteins with the variable target function method (Braun, W. and Go, N. (1985)J. Mol. Biol.,186, 611) is described, which makes use of redundant dihedral angle constraints (REDAC) derived from preliminary calculations of the complete structure. The REDAC approach reduces the computation time for obtaining a group of acceptable conformers with the program DIANA 5-100-fold, depending on the complexity of the protein structure, and retains good sampling of conformation space.Dedicated to the memory of Professor V.F. Bystrov     

Three-dimensional structure of the neurotoxin ATX Ia from Anemonia sulcata in aqueous solution determined by nuclear magnetic resonance spectroscopy

With the aid of 1H nuclear magnetic resonance (NMR) spectroscopy, the three-dimensional structure in aqueous solution was determined for ATX Ia, which is a 46 residue polypeptide neurotoxin of the sea anemone Anemonia sulcata. The input for the structure calculations consisted of 263 distance constraints from nuclear Overhauser effects (NOE) and 76 vicinal coupling constants. For the structure calculation several new or ammended programs were used in a revised strategy consisting of five successive computational steps. First, the program HABAS was used for a complete search of all backbone and chi 1 conformations that are compatible with the intraresidual and sequential NMR constraints. Second, using the program DISMAN, we extended this approach to pentapeptides by extensive sampling of all conformations that are consistent with the local and medium-range NMR constraints. Both steps resulted in the definition of additional dihedral angle constraints and in stereospecific assignments for a number of beta-methylene groups. In the next two steps DISMAN was used to obtain a group of eight conformers that contain no significant residual violations of the NMR constraints or van der Waals contacts. Finally, these structures were subjected to restrained energy refinement with a modified version of the molecular mechanics module of AMBER, which in addition to the energy force field includes potentials for the NOE distance constraints and the dihedral angle constraints. The average of the pairwise minimal RMS distances between the resulting refined conformers calculated for the well defined molecular core, which contains the backbone atoms of 35 residues and 20 interior side chains, is 1.5 +/- 0.3 A. This core is formed by a four-stranded beta-sheet connected by two well-defined loops, and there is an additional flexible loop consisting of the eleven residues 8-18. The core of the protein is stabilized by three disulfide bridges, which are surrounded by hydrophobic residues and shielded on one side by hydrophilic residues.     

The NMR solution structure of the pheromone Er-1 from the ciliated protozoan Euplotes raikovi.

The 3-dimensional structure of the pheromone Er-1 isolated from the ciliated protozoan Euplotes raikovi has been determined in aqueous solution by 1H NMR spectroscopy. The structure of this 40-residue protein was calculated with the distance geometry program DIANA on the basis of 503 upper distance constraints derived from nuclear Overhauser effects and 77 dihedral angle constraints derived from spin-spin coupling constants, and refined by restrained energy minimization with the program OPAL. The Er-1 solution structure is represented by a group of 20 conformers with an average RMS deviation relative to the mean structure of 0.55 A for the backbone atoms N, C alpha, and C', and 0.93 A for all heavy atoms of the complete polypeptide chain, residues 1-40. The molecular architecture is dominated by an up-down-up bundle of 3 alpha-helices formed by residues 2-9, 12-19, and 24-33. Although this core part coincides closely with the previously determined structure of the homologous pheromone Er-10, the C-terminal peptide segment adopts a novel conformation. This is of interest in view of previous suggestions, based on sequence comparisons, that this molecular region may be important for the different specificity of receptor recognition by different pheromones.     

Determination of a high-quality nuclear magnetic resonance solution structure of the bovine pancreatic trypsin inhibitor and comparison with three crystal structures.

A high-quality three-dimensional structure of the bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution was determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and compared to the three available high-resolution X-ray crystal structures. A newly collected input of 642 distance constraints derived from nuclear Overhauser effects and 115 dihedral angle constraints was used for the structure calculations with the program DIANA, followed by restrained energy minimization with the program AMBER. The BPTI solution structure is represented by a group of 20 conformers with an average root-mean-square deviation (RMSD) relative to the mean solution structure of 0.43 A for backbone atoms and 0.92 A for all heavy atoms of residues 2 to 56. The pairwise RMSD values of the three crystal structures relative to the mean solution structure are 0.76 to 0.85 A for the backbone atoms and 1.24 to 1.33 A for all heavy atoms of residues 2 to 56. Small local differences in backbone atom positions between the solution structure and the X-ray structures near residues 9, 25 to 27, 46 to 48 and 52 to 58, and conformational differences for individual amino acid side-chains were analyzed for possible correlations with intermolecular protein-protein contacts in the crystal lattices, using the pairwise RMSD values among the three crystal structures as a reference.     

Solution structure of the phosphocarrier protein HPr from Bacillus subtilis by two-dimensional NMR spectroscopy.

The solution structure of the phosphocarrier protein, HPr, from Bacillus subtilis has been determined by analysis of two-dimensional (2D) NMR spectra acquired for the unphosphorylated form of the protein. Inverse-detected 2D (1H-15N) heteronuclear multiple quantum correlation nuclear Overhauser effect (HMQC NOESY) and homonuclear Hartmann-Hahn (HOHAHA) spectra utilizing 15N assignments (reported here) as well as previously published 1H assignments were used to identify cross-peaks that are not resolved in 2D homonuclear 1H spectra. Distance constraints derived from NOESY cross-peaks, hydrogen-bonding patterns derived from 1H-2H exchange experiments, and dihedral angle constraints derived from analysis of coupling constants were used for structure calculations using the variable target function algorithm, DIANA. The calculated models were refined by dynamical simulated annealing using the program X-PLOR. The resulting family of structures has a mean backbone rmsd of 0.63 A (N, C alpha, C', O atoms), excluding the segments containing residues 45-59 and 84-88. The structure is comprised of a four-stranded antiparallel beta-sheet with two antiparallel alpha-helices on one side of the sheet. The active-site His 15 residue serves as the N-cap of alpha-helix A, with its N delta 1 atom pointed toward the solvent to accept the phosphoryl group during the phosphotransfer reaction with enzyme I. The existence of a hydrogen bond between the side-chain oxygen atom of Tyr 37 and the amide proton of Ala 56 is suggested, which may account for the observed stabilization of the region that includes the beta-turn comprised of residues 37-40. If the beta alpha beta beta alpha beta (alpha) folding topology of HPr is considered with the peptide chain polarity reversed, the protein fold is identical to that described for another group of beta alpha beta beta alpha beta proteins that include acylphosphatase and the RNA-binding domains of the U1 snRNP A and hnRNP C proteins.     

Determination of a high precision structure of a novel protein, Linum usitatissimum trypsin inhibitor (LUTI), using computer-aided assignment of NOESY cross-peaks

The solution structure of a novel 69 residue proteinase inhibitor, Linum usitatissimum trypsin inhibitor (LUTI), was determined using a method based on computer aided assignment of nuclear Overhauser enhancement spectroscopy (NOESY) data. The approach applied uses the program NOAH/DYANA for automatic assignment of NOESY cross-peaks. Calculations were carried out using two unassigned NOESY peak lists and a set of determined dihedral angle restraints. In addition, hydrogen bonds involving amide protons were identified during calculations using geometrical criteria and values of HN temperature coefficients. Stereospecific assignment of beta-methylene protons was carried out using a standard procedure based on nuclear Overhauser enhancement intensities and 3J(alpha)(beta) coupling constants. Further stereospecific assignment of methylene protons and diastereotopic methyl groups were established upon structure-based method available in the program GLOMSA and chemical shift calculations. The applied algorithm allowed us to assign 1968 out of 2164 peaks (91%) derived from NOESY spectra recorded in H2O and 2H2O. The final experimental data input consisted of 1609 interproton distance restraints, 88 restraints for 44 hydrogen bonds, 63 torsion angle restraints and 32 stereospecifically assigned methylene proton pairs and methyl groups. The algorithm allowed the calculation of a high precision protein structure without the laborious manual assignment of NOESY cross-peaks. For the 20 best conformers selected out of 40 refined ones in the program CNS, the calculated average pairwise rmsd values for residues 3 to 69 were 0.38 A (backbone atoms) and 1.02 A (all heavy atoms). The three-dimensional LUTI structure consists of a mixed parallel and antiparallel beta-sheet, a single alpha-helix and shows the fold of the potato 1 family of proteinase inhibitors. Compared to known structures of the family, LUTI contains Arg and Trp residues at positions P6' and P8', respectively, instead of two Arg residues, involved in the proteinase binding loop stabilization. A consequence of the ArgTrp substitution at P8' is a slightly more compact conformation of the loop relative to the protein core.     

The Nuclear Magnetic Resonance Solution Structure of the Mixed Disulfide between Escherichia coli Glutaredoxin(C14S) and Glutathione

The determination of the nuclear magnetic resonance (NMR) solution structure of the mixed disulfide between the mutant Escherichia coli glutaredoxin Grx(C14S) and glutathione (GSH), Grx(C14S)-SG, is described, the binding site for GSH on Grx(C14S) is located, and the non-bonding interactions between -SG and the protein are characterized. Based on nearly complete sequence-specific NMR assignments, 1010 nuclear Overhauser enhancement upper distance constraints and 116 dihedral angle constraints were obtained as the input for the structure calculations, for which the distance geometry program DIANA was used followed by energy minimization in a waterbath with the AMBER force field in the program OPAL. The -SG moiety was found to be localized on the surface of the protein in a cleft bounded by the amino acid residues Y13, T58, V59, Y72, T73 and D74. Hydrogen bonds have been identified between -SG and the residues V59 and T73 of Grx(C14S), and the formation of an additional hydrogen bond with Y72 and electrostatic interactions with the side-chains of D74 and K45 are also compatible with the NMR, conformational constraints. Comparison of the reduced and oxidized forms of Grx with Grx(C14S)-SG shows that the mixed disulfide more closely resembles the oxidized form of the protein. Functional implications of this observation are discussed. Comparisons are also made with the related proteins bacteriophage T4 glutaredoxin and glutathione S-transferase.     

The NMR solution structure of the pheromone Er-11 from the ciliated protozoan Euplotes raikovi.

The NMR solution structure of the pheromone Er-11, a 39-residue protein from the ciliated protozoan Euplotes raikovi, was calculated with the distance geometry program DIANA from 449 NOE upper distance constraints and 97 dihedral angle constraints, and the program OPAL was employed for structure refinement by molecular mechanics energy minimization in a water bath. For a group of 20 conformers used to characterize the solution structure, the average of the pairwise RMS deviations from the mean structure calculated for the backbone heavy atoms N, C alpha, and C' of residues 2-38 was 0.30 A. The molecular architecture is dominated by an up-down-up bundle of three short helices with residues 2-9, 12-19, and 22-32, which is closely similar to the previously determined structures of the homologous pheromones Er-1, Er-2, and Er-10. This finding provides structural evidence for the capability shown by these pheromones to compete with each other in binding reactions to their cell-surface receptors.     

Stereospecific assignments in proteins using exact NOEs

Recently developed methods to measure distances in proteins with high accuracy by “exact” nuclear Overhauser effects (eNOEs) make it possible to determine stereospecific assignments, which are particularly important to fully exploit the accuracy of the eNOE distance measurements. Stereospecific assignments are determined by comparing the eNOE-derived distances to protein structure bundles calculated without stereospecific assignments, or an independently determined crystal structure. The absolute and relative CYANA target function difference upon swapping the stereospecific assignment of a diastereotopic group yields the respective stereospecific assignment. We applied the method to the eNOE data set that has recently been obtained for the third immunoglobulin-binding domain of protein G (GB3). The 884 eNOEs provide relevant data for 47 of the total of 75 diastereotopic groups. Stereospecific assignments could be established for 45 diastereotopic groups (96 %) using the X-ray structure, or for 27 diastereotopic groups (57 %) using structures calculated with the eNOE data set without stereospecific assignments, all of which are in agreement with those determined previously. The latter case is relevant for structure determinations based on eNOEs. The accuracy of the eNOE distance measurements is crucial for making stereospecific assignments because applying the same method to the traditional NOE data set for GB3 with imprecise upper distance bounds yields only 13 correct stereospecific assignments using the X-ray structure or 2 correct stereospecific assignments using NMR structures calculated without stereospecific assignments.     

Complete relaxation matrix refinement of NMR structures of proteins using analytically calculated dihedral angle derivatives of NOE intensities

Summary A new method for refining three-dimensional (3D) NMR structures of proteins is described, which takes account of the complete relaxation pathways. Derivatives of the NOE intensities with respect to the dihedral angles are analytically calculated, and efficiently evaluated with the use of a filter technique for identifying the dominant terms of these derivatives. This new method was implemented in the distance geometry program DIANA. As an initial test, we refined 30 rigid distorted helical structures, using a simulated data set of NOE distance constraints for a rigid standard -helix. The final root-mean-square deviations of the refined structures relative to the standard helix were less than 0.1 Å, and the R-factors dropped from values between 7% and 32% to values of less than 0.5% in all cases, which compares favorably with the results from distance geometry calculations. In particular, because spin diffusion was not explicitly considered in the evaluation of exact¹H–¹H distances corresponding to the simulated NOE intensities, a group of nearly identical distance geometry structures was obtained which had about 0.5 Å root-mean-square deviation from the standard -helix. Further test calculations using an experimental NOE data set recorded for the protein trypsin inhibitor K showed that the complete relaxation matrix refinement procedure in the DIANA program is functional also with systems of practical interest.Abbreviations RMSD root-mean-square deviation - NOE nuclear Overhauser enhancement - NOESY 2-dimensional nuclear Overhauser enhancement spectroscopy - CPU central processing unit     

Sequence-specific 1H n.m.r. assignments and determination of the three-dimensional structure of reduced Escherichia coli glutaredoxin

The determination of the nuclear magnetic resonance structure of reduced E. coli glutaredoxin in aqueous solution is described. Based on nearly complete, sequence-specific resonance assignments, 813 nuclear Overhauser effect distance constraints and 191 dihedral angle constraints were employed as the input for the structure calculations, for which the distance geometry program DIANA was used followed by simulated annealing with the program X-PLOR. The molecular architecture of reduced glutaredoxin is made up of three helices and four-stranded beta-sheet. The first strand of the beta-sheet (residues 2 to 7) runs parallel to the second strand (32 to 37) and antiparallel to the third strand (61 to 64), and the sheet is extended in an antiparallel fashion with a fourth strand (67 to 69). The first helix with residues 13 to 28 and the last helix (71 to 83) run parallel to each other on one side of the beta-sheet, with their direction opposite to that of the two parallel beta-strands, and the helix formed by residues 44 to 53 fills space available due to the twist of the beta-sheet and the reduced length of the last two beta-strands. The active site Cys11-Pro-Tyr-Cys14 is located after the first beta-strand and occupies the latter part of the loop connecting this strand with the first helix.     

Solution structures of the rabbit neutrophil defensin NP-5

Solution structures of the rabbit neutrophil defensin NP-5 have been determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and distance geometry techniques. This 33 amino acid peptide is part of the oxygen-independent mammalian defense system against microbial infection. The structures were generated from 107 n.m.r. derived inter-residue proton-proton distance constraints. A distance geometry algorithm was then used to determine the range of structures consistent with these distance constraints. These distance geometry calculations employed an improved algorithm that allowed the chirality constraints to be relaxed on prochiral centers when it was not possible to make stereo-specific assignments of protons on these centers. This procedure gave superior results compared with standard distance geometry methods and also produced structures that were more consistent with the original n.m.r. data. Analysis of the NP-5 structures shows that the overall folding of the peptide backbone is well defined by the n.m.r. distance information but that the side-chain group conformations are generally less well defined.     

Three-dimensional structure of ectatomin from Ectatomma tuberculatum ant venom

Summary Two-dimensional ¹H NMR techniques were used to determine the spatial structure of ectatomin, a toxin from the venom of the ant Ectatomma tuberculatum. Nearly complete proton resonance assignments for two chains of ectatomin (37 and 34 amino acid residues, respectively) were obtained using 2D TOCSY, DQF-COSY and NOESY experiments. The cross-peak volumes in NOESY spectra were used to define the local structure of the protein and generate accurate proton-proton distance constraints employing the MARDIGRAS program. Disulfide bonds were located by analyzing the global fold of ectatomin, calculated with the distance geometry program DIANA. These data, combined with data on the rate of exchange of amide protons with deuterium, were used to obtain a final set of 20 structures by DIANA. These structures were refined by unrestrained energy minimization using the CHARMm program. The resulting rms deviations over 20 structures (excluding the mobile N- and C-termini of each chain) are 0.75 ? for backbone heavy atoms, and 1.25 ? for all heavy atoms. The conformations of the two chains are similar. Each chain consists of two α-helices and a hinge region of four residues; this forms a hairpin structure which is stabilized by disulfide bridges. The hinge regions of the two chains are connected together by a third disulfide bridge. Thus, ectatomin forms a four-α-helical bundle structure.     

Extensive distance geometry calculations with different NOE calibrations: New criteria for structure selection applied to Sandostatin and BPTI

Summary To generate structures efficiently, a version of the distance geometry program DIANA for a parallel computer was developed, new objective criteria for the selection of NMR solution structures are presented, and the influence of using different calibrations of NOE intensities on the final structures are described. The methods are applied to the structure determination of Sandostatin, a disulfide-bridge octapeptide, and to model calculations of BPTI. On an Alliant FX2800 computer using 10 processors in parallel, the calculations were done 9.2 times faster than with a single processor. Up to 7000 Sandostatin structures were calculated with distance and angular constraints. The procedure for selecting acceptable structures is based on the maximum values of pairwise RMSDs between structures. Suitable target function cut-offs are defined independent of the number of starting structures. The method allowed for an objective comparison of three groups of Sandostatin structures that were calculated from different sets of upper distance constraints which were derived from the same NOE intensity data using three empirical calibration curves. The number of converged structures and the target function values differed significantly among the three groups, but the structures were qualitatively and quantitatively very similar. The conformation is well determined in the cyclic region Cys²–Cys⁷ and adopts a -turn centered at d-Trp⁴–Lys⁵. The criteria for structure selection were further tested with BPTI. Results obtained from sets of structures calculated with and without using the REDAC strategy are consistent and suggest that the structure selection method is objective and generally applicable.     

Determination of the complete three-dimensional structure of the alpha-amylase inhibitor tendamistat in aqueous solution by nuclear magnetic resonance and distance geometry

The complete three-dimensional structure of the alpha-amylase inhibitor Tendamistat in aqueous solution was determined by 1H nuclear magnetic resonance and distance geometry calculations using the program DISMAN. Compared to an earlier, preliminary determination of the polypeptide backbone conformation, stereo-specific assignments were obtained for 41 of the 89 prochiral groups in the protein, and a much more extensive set of experimental constraints was collected, including 842 distance constraints from nuclear Overhauser effects and over 100 supplementary constraints from spin-spin coupling constants and the identification of intramolecular hydrogen bonds. The complete protein molecule, including the amino acid side-chains is characterized by a group of nine structures corresponding to the results of the nine DISMAN calculations with minimal residual error functions. The average of the pairwise minimal root-mean-square distances among these nine structures is 0.85 A for the polypeptide backbone, and 1.52 A for all the heavy atoms. The procedures used for the structure determination are described and a detailed analysis is presented of correlations between the experimental input data and the precision of the structure determination.     

Nuclear magnetic resonance solution structure of the alpha-neurotoxin from the black mamba (Dendroaspis polylepis polylepis).

The three-dimensional structure in solution of the alpha-neurotoxin from the black mamba (Dendroaspis polylepis polylepis) has been determined by nuclear magnetic resonance spectroscopy. A high quality structure for this 60-residue protein was obtained from 656 NOE distance constraints and 143 dihedral angle constraints, using the distance geometry program DIANA for the structure calculation and AMBER for restrained energy minimization. For a group of 20 conformers used to represent the solution structure, the average root-mean-square deviation value calculated for the polypeptide backbone heavy atoms relative to the mean structure was 0.45 A. The protein consists of a core region from which three finger-like loops extend outwards. It includes a short, two-stranded antiparallel beta-sheet of residues 1-5 and 13-17, a three-stranded antiparallel beta-sheet involving residues 23-31, 34-42 and 51-55, and four disulfide bridges in the core region. There is also extensive non-regular hydrogen bonding between the carboxy-terminal tail of the polypeptide chain and the rest of the core region. Comparison with the crystal structure of erabutoxin-b indicates that the structure of alpha-neurotoxin is quite similar to other neurotoxin structures, but that local structural differences are seen in regions thought to be important for binding of neurotoxins to the acetylcholine receptor. For two regions of the alpha-neurotoxin structure there is evidence for an equilibrium between multiple conformations, which might be related to conformational rearrangements upon binding to the receptor. Overall, the alpha-neurotoxin presents itself as a protein with a stable core and flexible surface areas that interact with the acetylcholine receptor in such a way that high affinity binding is achieved by conformational rearrangements of the deformable regions of the neurotoxin structure.     

NMR structure of oxidized Escherichia coli glutaredoxin: comparison with reduced E. coli glutaredoxin and functionally related proteins.

The determination of the NMR structure of oxidized Escherichia coli glutaredoxin in aqueous solution is described, and comparisons of this structure with that of reduced E. coli glutaredoxin and the related proteins E. coli thioredoxin and T4 glutaredoxin are presented. Based on nearly complete sequence-specific 1H-NMR assignments, 804 nuclear Overhauser enhancement distance constraints and 74 dihedral angle constraints were obtained as the input for the structure calculations, for which the distance geometry program DIANA was used followed by simulated annealing with the program X-PLOR. The molecular architecture of oxidized glutaredoxin is made up of three helices and a four-stranded beta-sheet. The three-dimensional structures of oxidized and the recently described reduced glutaredoxin are very similar. Quantitative analysis of the exchange rates of 34 slowly exchanging amide protons from corresponding series of two-dimensional [15N,1H]-correlated spectra of oxidized and reduced glutaredoxin showed close agreement, indicating almost identical hydrogen-bonding patterns. Nonetheless, differences in local dynamics involving residues near the active site and the C-terminal alpha-helix were clearly manifested. Comparison of the structure of E. coli glutaredoxin with those of T4 glutaredoxin and E. coli thioredoxin showed that all three proteins have a similar overall polypeptide fold. An area of the protein surface at the active site containing Arg 8, Cys 11, Pro 12, Tyr 13, Ile 38, Thr 58, Val 59, Pro 60, Gly 71, Tyr 72, and Thr 73 is proposed as a possible site for interaction with other proteins, in particular ribonucleotide reductase. It was found that this area corresponds to previously proposed interaction sites in T4 glutaredoxin and E. coli thioredoxin. The solvent-accessible surface area at the active site of E. coli glutaredoxin showed a general trend to increase upon reduction. Only the sulfhydryl group of Cys 11 is exposed to the solvent, whereas that of Cys 14 is buried and solvent inaccessible.     

Calculation of protein conformations by proton-proton distance constraints. A new efficient algorithm

We have developed a method to determine the three-dimensional structure of a protein molecule from such a set of distance constraints as can be determined by nuclear magnetic resonance studies. The currently popular methods for distance geometry based on the use of the metric matrix are applicable only to small systems. The method developed here is applicable to large molecules, such as proteins, with all atoms treated explicitly. This method works in the space of variable dihedral angles and determines a three-dimensional structure by minimization of a target function. We avoid difficulties hitherto inherent in this type of approach by two new devices: the use of variable target functions; and a method of rapid calculation of the gradient of the target functions. The method is applied to the determination of the structures of a small globular protein, bovine pancreatic trypsin inhibitor, from several artificial sets of distance constraints extracted from the X-ray crystal structure of this molecule. When a good set of constraints was available for both short- and long-range distances, the crystal structure was regenerated nearly exactly. When some ambiguities, such as those expected in experimental information, are allowed, the protein conformation can be determined up to a few local deformations. These ambiguities are mainly associated with the low resolving power of the short-range information.