Three-dimensional structure of potato carboxypeptidase inhibitor in solution. A study using nuclear magnetic resonance, distance geometry, and restrained molecular dynamics

The solution conformation of potato carboxypeptidase inhibitor (CPI) has been investigated by 1H NMR spectroscopy. The spectrum is assigned in a sequential manner by using two-dimensional NMR techniques to identify through-bond and through-space (less than 5 A) connectivities. A set of 309 approximate interproton distance restraints is derived from the two-dimensional nuclear Overhauser enhancement spectra and used as the basis of a three-dimensional structure determination by a combination of metric matrix distance geometry and restrained molecular dynamics calculations. A total of 11 converged distance geometry structures were computed and refined by using restrained molecular dynamics. The average atomic root mean square (rms) difference between the final 11 structures and the mean structure obtained by averaging their coordinates is 1.4 +/- 0.3 A for residues 2-39 and 0.9 +/- 0.2 A for residues 5-37. The corresponding values for all atoms are 1.9 +/- 0.3 and 1.4 +/- 0.2 A, respectively. The larger values for residues 2-38 relative to those for residues 5-37 arise from the fact that the positions of the N- (residues 1-4) and C- (residues 38-39) terminal tails are rather poorly determined, whereas those of the core of the protein (residues 5-37) are well determined by the experimental interproton distance data. The computed structures are very close to the X-ray structure of CPI in its complex with carboxypeptidase, and the backbone atomic rms difference between the mean of the computed structures and the X-ray structure is only 1.2 A. Nevertheless, there are some real differences present which are evidenced by significant deviations between the experimental upper interproton distance limits and the corresponding interproton distances derived from the X-ray structure. These principally occur in two regions, residues 18-20 and residues 28-30, the latter comprising part of the region of secondary contacts between CPI and carboxypeptidase in the X-ray structure.     

Determination of the three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata: a study using nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing

The three-dimensional solution structure of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata has been determined on the basis of 489 interproton and 24 hydrogen-bonding distance restraints supplemented by 23 phi backbone and 21 chi 1 side-chain torsion angle restraints derived from nuclear magnetic resonance (NMR) measurements. A total of 42 structures is calculated by a hybrid metric matrix distance geometry-dynamical simulated annealing approach. Both the backbone and side-chain atom positions are well defined. The average atomic rms difference between the 42 individual SA structures and the mean structure obtained by averaging their coordinates is 0.67 +/- 0.12 A for the backbone atoms and 0.90 +/- 0.17 A for all atoms. The core of the protein is formed by a triple-stranded antiparallel beta-sheet composed of residues 14-16 (strand 1), 30-34 (strand 2), and 37-41 (strand 3) with an additional mini-antiparallel beta-sheet at the N-terminus (residues 6-9). The first and second strands of the triple-stranded antiparallel beta-sheet are connected by a long exposed loop (residues 17-30). A number of side-chain interactions are discussed in light of the structure.     

The determination of the three-dimensional structure of barley serine proteinase inhibitor 2 by nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The solution structure of the 64 residue structured domain (residues 20-83) of barley serine proteinase inhibitor 2 (BSPI-2) is determined on the basis of 403 interproton distance, 34 phi backbone torsion angle and 26 hydrogen bonding restraints derived from n.m.r. measurements. A total of 11 converged structures were computed using a metric matrix distance geometry algorithm and refined by restrained molecular dynamics. The average rms difference between the final 11 structures and the mean structure obtained by averaging their coordinates is 1.4 +/- 0.2 A for the backbone atoms and 2.1 +/- 0.1 A for all atoms. The overall structure, which is almost identical to that found by X-ray crystallography, is disc shaped and consists of a central four component mixed parallel and antiparallel beta-sheet flanked by a 13 residue alpha-helix on one side and the reactive site loop on the other.     

The influence of stereospecific assignments on the determination of three-dimensional structures of proteins by nuclear magnetic resonance spectroscopy. Application to the sea anemone protein BDS-I

The influence of the stereospecific assignments of beta-methylene protons and the classification of chi 1 torsion angles on the definition of the three-dimensional structures of proteins determined from NMR data is investigated using the sea anemone protein BDS-I (43 residues) as a model system. Two sets of structures are computed. The first set comprises 42 converged structures (denoted STEREO structures) calculated on the basis of the complete list of restraints derived from the NMR data, consisting of 489 interproton and 24 hydrogen bonding distance restraints, supplemented by 23 phi backbone and 21 chi 1 side chain torsion angle restraints. The second set comprises 31 converged structures (denoted NOSTEREO structures) calculated from a reduced data set in which those restraints arising from stereospecific assignments, and the corresponding chi 1 torsion angle restraints, are explicitly omitted. The results show that the inclusion of the stereospecific restraints leads to a significant improvement in the definition of the structure of BDS-I, both with respect to the backbone and the detailed arrangement of the side chains. Average atomic rms differences between the individual structures and the mean structures for the backbone atoms are 0.67 +/- 0.12 A and 0.93 +/- 0.16 A for the STEREO and NOSTEREO structures, respectively; the corresponding values for all atoms are 0.90 +/- 0.17 A and 1.17 +/- 0.17 A, respectively. In addition, while the overall fold remains unchanged, there is a small but significant atomic displacement between the two sets of structures.     

Application of molecular dynamics with interproton distance restraints to three-dimensional protein structure determination. A model study of crambin

The applicability of restrained molecular dynamics for the determination of three-dimensional protein structures on the basis of short interproton distances (less than 4 A) that can be realistically determined from nuclear magnetic resonance measurements in solution is assessed. The model system used is the 1.2 A resolution crystal structure of the 46 residue protein crambin, from which a set of 240 approximate distance restraints, divided into three ranges (2.5 +/- 0.5, 3.0+0.5(-1.0) and 4 +/- 1 A), is derived. This interproton distance set comprises 159 short-range ([i-j] less than or equal to 5) and 56 ([i-j] greater than 5) long-range inter-residue distances and 25 intra-residue distances. Restrained molecular dynamics are carried out using a number of different protocols starting from two initial structures: a completely extended beta-strand; and an extended structure with two alpha-helices in the same positions as in the crystal structure (residues 7 to 19, and 23 to 30) and all other residues in the form of extended beta-strands. The root-mean-square (r.m.s.) atomic differences between these two initial structures and the crystal structure are 43 A and 23 A, respectively. It is shown that, provided protocols are used that permit the secondary structure elements to form at least partially prior to folding into a tertiary structure, convergence to the correct final structure, both globally and locally, is achieved. The r.m.s. atomic differences between the converged restrained dynamics structures and the crystal structure range from 1.5 to 2.2 A for the backbone atoms and from 2.0 to 2.8 A for all atoms. The r.m.s. atomic difference between the X-ray structure and the structure obtained by first averaging the co-ordinates of the converged restrained dynamics structures is even smaller: 1.0 A for the backbone atoms and 1.6 A for all atoms. These results provide a measure with which to judge future experimental results on proteins whose crystal structures are unknown. In addition, from an examination of the dynamics trajectories, it is shown that the convergence pathways followed by the various simulations are different.     

Solution structure of a polypeptide dimer comprising the fourth Ca(2+)-binding site of troponin C by nuclear magnetic resonance spectroscopy

The structure of a 39 amino acid proteolytic fragment of rabbit skeletal troponin C containing the fourth Ca(2+)-binding site has been determined by an approach involving nuclear magnetic resonance (NMR) spectroscopy combined with hybrid distance geometry-dynamical simulated annealing calculations. Hydrodynamic and NMR evidence establishes unambiguously that the fragment forms a stable dimer in solution in the presence of excess Ca2+. The calculation of the dimeric structure is based on a total of 1056 experimental restraints comprising 422 interproton distances, 35 phi, 28 psi, and 28 chi 1 torsion angle restraints within each subunit, 30 intermonomer distance restraints, and 6 Ca2+ restraints per subunit. A total of 48 final structures were calculated having an rms deviation about the mean atomic backbone coordinate positions of 1.0 A for residues Asp128-Glu156. The solution structure consists of a dimer of helix-loop-helix motifs related by a 2-fold axis of symmetry. The overall architecture of the dimer is very similar to the C-terminal domain in the crystal structure of chicken skeletal troponin C.     

Determinants of protein hyperthermostability: purification and amino acid sequence of rubredoxin from the hyperthermophilic archaebacterium Pyrococcus furiosus and secondary structure of the zinc adduct by NMR

The purification, amino acid sequence, and two-dimensional 1H NMR results are reported for the rubredoxin (Rd) from the hyperthermophilic archaebacterium Pyrococcus furiosus, an organism that grows optimally at 100 degrees C. The molecular mass (5397 Da), iron content (1.2 +/- 0.2 g-atom of Fe/mol), UV-vis spectrophotometric properties, and amino acid sequence (60% sequence identity with Clostridium pasteurianum Rd) are found to be typical of this class of redox protein. However, P. furiosus Rd is remarkably thermostable, being unaffected after incubation for 24 h at 95 degrees C. One- and two-dimensional 1H nuclear magnetic resonance spectra of the oxidized [Fe(III)Rd] and reduced [Fe(II)Rd] forms of P. furiosus Rd exhibited substantial paramagnetic line broadening, and this precluded detailed 3D structural studies. The apoprotein was not readily amenable to NMR studies due to apparent protein oxidation involving the free cysteine sulfhydryls. However, high-quality NMR spectra were obtained for the Zn-substituted protein, Zn(Rd), enabling detailed NMR signal assignment for all backbone amide and alpha and most side-chain protons. Secondary structural elements were determined from qualitative analysis of 2D Overhauser effect spectra. Residues A1-K6, Y10-E14, and F48-E51 form a three-strand antiparallel beta-sheet, which comprises ca. 30% of the primary sequence. Residues C5-Y10 and C38-A43 form types I and II amide-sulfur tight turns common to iron-sulfur proteins. These structural elements are similar to those observed by X-ray crystallography for native Rd from the mesophile C. pasteurianum. However, the beta-sheet domain in P. furiosus Rd is larger than that in C. pasteurianum Rd and appears to begin at the N-terminal residue. From analysis of the secondary structure, potentially stabilizing electrostatic interactions involving the charged groups of residues Ala(1), Glu(14), and Glu(52) are proposed. These interactions, which are not present in rubredoxins from mesophilic organisms, may prevent the beta-sheet from "unzipping" at elevated temperatures.     

Determination of three-dimensional structures of proteins by simulated annealing with interproton distance restraints. Application to crambin, potato carboxypeptidase inhibitor and barley serine proteinase inhibitor 2

An automated method, based on the principle of simulated annealing, is presented for determining the three-dimensional structures of proteins on the basis of short (less than 5 A) interproton distance data derived from nuclear Overhauser enhancement (NOE) measurements. The method makes use of Newton's equations of motion to increase temporarily the temperature of the system in order to search for the global minimum region of a target function comprising purely geometric restraints. These consist of interproton distances supplemented by bond lengths, bond angles, planes and soft van der Waals repulsion terms. The latter replace the dihedral, van der Waals, electrostatic and hydrogen-bonding potentials of the empirical energy function used in molecular dynamics simulations. The method presented involves the implementation of a number of innovations over our previous restrained molecular dynamics approach [Clore, G.M., Brünger, A.T., Karplus, M. and Gronenborn, A.M. (1986) J. Mol. Biol., 191, 523-551]. These include the development of a new effective potential for the interproton distance restraints whose functional form is dependent on the magnitude of the difference between calculated and target values, and the design and implementation of robust and fully automatic protocol. The method is tested on three systems: the model system crambin (46 residues) using X-ray structure derived interproton distance restraints, and potato carboxypeptidase inhibitor (CPI; 39 residues) and barley serine proteinase inhibitor 2 (BSPI-2; 64 residues) using experimentally derived interproton distance restraints. Calculations were carried out starting from the extended strands which had atomic r.m.s. differences of 57, 38 and 33 A with respect to the crystal structures of BSPI-2, crambin and CPI respectively. Unbiased sampling of the conformational space consistent with the restraints was achieved by varying the random number seed used to assign the initial velocities. This ensures that the different trajectories diverge during the early stages of the simulations and only converge later as more and more interproton distance restraints are satisfied. The average backbone atomic r.m.s. difference between the converged structures is 2.2 +/- 0.3 A for crambin (nine structures), 2.4 +/- 0.3 A for CPI (eight structures) and 2.5 +/- 0.2 A for BSPI-2 (five structures). The backbone atomic r.m.s. difference between the mean structures derived by averaging the coordinates of the converged structures and the corresponding X-ray structures is 1.2 A for crambin, 1.6 A for CPI and 1.7 A for BSPI-2.     

Determination of the three-dimensional solution structure of barnase using nuclear magnetic resonance spectroscopy

The solution conformation of the ribonuclease barnase has been determined by using 1H nuclear magnetic resonance (NMR) spectroscopy. The 20 structures were calculated by using 853 interproton distance restraints obtained from analyses of two-dimensional nuclear Overhauser spectra, 72 phi and 53 chi 1 torsion angle restraints, and 17 hydrogen-bond distance restraints. The calculated structures contain two alpha-helices (residues 6-18 and 26-34) and a five-stranded antiparallel beta-sheet (residues 50-55, 70-75, 85-91, 94-101, and 105-108). The core of the protein is formed by the packing of one of the alpha-helices (residues 6-18) onto the beta-sheet. The average RMS deviation between the calculated structures and the mean structure is 1.11 A for the backbone atoms and 1.75 A for all atoms. The protein is least well-defined in the N-terminal region and in three large loops. When these regions are excluded, the average RMS deviation between the calculated structures and the mean structure for residues 5-34, 50-56, 71-76, 85-109 is 0.62 A for the backbone atoms and 1.0 A for all atoms. The NMR-derived structure has been compared with the crystal structure of barnase [Mauguen et al. (1982) Nature (London) 297, 162-164].     

The polypeptide fold of the globular domain of histone H5 in solution. A study using nuclear magnetic resonance, distance geometry and restrained molecular dynamics.

The polypeptide fold of the 79-residue globular domain of chicken histone H5 (GH5) in solution has been determined by the combined use of distance geometry and restrained molecular dynamics calculations. The structure determination is based on 307 approximate interproton distance restraints derived from n.m.r. measurements. The structure is composed of a core made up of residues 3-18, 23-34, 37-60 and 71-79, and two loops comprising residues 19-22 and 61-70. The structure of the core is well defined with an average backbone atomic r.m.s. difference of 2.3 +/- 0.3 A between the final eight converged restrained dynamics structures and the mean structure obtained by averaging their coordinates best fitted to the core residues. The two loops are also well defined locally but their orientation with respect to the core could not be determined as no long range ([i-j[ greater than 5) proton-proton contacts could be observed between the loop and core residues in the two-dimensional nuclear Overhauser enhancement spectra. The structure of the core is dominated by three helices and has a similar fold to the C-terminal DNA binding domain of the cAMP receptor protein.     

Analysis of the relative contributions of the nuclear Overhauser interproton distance restraints and the empirical energy function in the calculation of oligonucleotide structures using restrained molecular dynamics

The relative contributions of the interproton distance restraints derived from nuclear Overhauser enhancement measurements and of the empirical energy function in the determination of oligonucleotide structures by restrained molecular dynamics are investigated. The calculations are based on 102 intraresidue and 126 interresidue interproton distance restraints derived from short mixing time two-dimensional nuclear Overhauser enhancement data on the dodecamer 5'd(CGCGPATTCGCG)2 [Clore, G.M., Oschkinat, H., McLaughlin, L.W., Benseler, F., Scalfi Happ, C., Happ, E., & Gronenborn, A.M. (1988) Biochemistry 27, 4185-4197]. Eight interproton distance restraint lists were made up with errors ranging from -0.1/+0.2 to -1.2/+1.3 A for r less than 2.5 A and from -0.2/+0.3 to -1.3/+1.4 A for r greater than or equal to 2.5 A. These restraints were incorporated into the total energy function of the system in the form of square-well potentials with force constants set sufficiently high to ensure that the deviations between calculated distances and experimental restraints were very small (average interproton distance rms deviation of less than 0.06 A). For each data set, six calculations were carried out, three starting from classical A-DNA and three from classical B-DNA. The results show that structural changes occurring during the course of restrained molecular dynamics and the degree of structural convergence are determined by the interproton distance restraints. All the structures display similar small deviations from idealized geometry and have the same values for the nonbonding energy terms comprising van der Waals, electrostatic, and hydrogen-bonding components. Thus, the function of the empirical energy function is to maintain near perfect stereochemistry and nonbonded interactions.(ABSTRACT TRUNCATED AT 250 WORDS)     

Determination of the spatial structure of insectotoxin 15A from Buthus erpeus by (1)H-NMR spectroscopy data]

The solution structure of insectotoxin 15A (35 residues) from scorpion Buthus eupeus was determined on the basis of 386 interproton distance restraints 12 hydrogen-bonding restraints and 113 dihedral angle restraints derived from 1H NMR experiments. A group of 20 structures was calculated with the distance geometry program DIANA followed by the restrained energy minimization with the program CHARMM. The atomic RMS distribution about the mean coordinate position is 0.64 +/- 0.11 A for the backbone atoms and 1.35 +/- 0.20 A for all atoms. The structure contains an alpha-helix (residues 10-20) and a three-stranded antiparallel beta-sheet (residues 2-5, 24-28 and 29-33). A pairing of the eight cysteine residues of insectotoxin 15A was established basing on NMR data. Three disulfide bridges (residues 2-19, 16-31 and 20-33) connect the alpha-helix with the beta-sheet, and the fourth one (5-26) joins beta-strands together. The spatial fold of secondary structure elements (the alpha-helix and the beta-sheet) of the insectotoxin 15A is very similar to those of the other short and long scorpion toxins in spite of a low (about 20%) sequence homology.     

High-resolution three-dimensional structure of reduced recombinant human thioredoxin in solution

The solution structure of recombinant human thioredoxin (105 residues) has been determined by nuclear magnetic resonance (NMR) spectroscopy combined with hybrid distance geometry-dynamical simulated annealing calculations. Approximate interproton distance restraints were derived from nuclear Overhauser effect (NOE) measurements. In addition, a large number of stereospecific assignments for beta-methylene protons and torsion angle restraints for phi, psi, and chi 1 were obtained by using a conformational grid search on the basis of the intraresidue and sequential NOE data in conjunction with 3JHN alpha and 3J alpha beta coupling constants. The structure calculations were based on 1983 approximate interproton distance restraints, 52 hydrogen-bonding restraints for 26 hydrogen bonds, and 98 phi, 71 psi, and 72 chi 1 torsion angle restraints. The 33 final simulated annealing structures obtained had an average atomic rms distribution of the individual structures about the mean coordinate positions of 0.40 +/- 0.06 A for the backbone atoms and 0.78 +/- 0.05 A for all atoms. The solution structure of human thioredoxin consists of a five-stranded beta-sheet surrounded by four alpha-helices, with an active site protrusion containing the two redox-active cysteines. The overall structure is similar to the crystal and NMR structures of oxidized [Katti, S. K., LeMaster, D. M., & Eklund, H. (1990) J. Mol. Biol. 212, 167-184] and reduced [Dyson, J. H., Gippert, G. P., Case, D. A., Holmgren, A., & Wright, P. (1990) Biochemistry 29, 4129-4136] Escherichia coli thioredoxin, respectively, despite the moderate 25% amino acid sequence homology. Several differences, however, can be noted. The human alpha 1 helix is a full turn longer than the corresponding helix in E. coli thioredoxin and is characterized by a more regular helical geometry. The helix labeled alpha 3 in human thioredoxin has its counterpart in the 3(10) helix of the E. coli protein and is also longer in the human protein. In contrast to these structural differences, the conformation of the active site loop in both proteins is very similar, reflecting the perfect sequence identity for a stretch of eight amino acid residues around the redox-active cysteines.     

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa 2. Three-dimensional solution structure.

The three-dimensional solution structure of conotoxin TVIIA, a 30-residue polypeptide from the venom of the piscivorous cone snail Conus tulipa, has been determined using 2D 1H NMR spectroscopy. TVIIA contains six cysteine residues which form a 'four-loop' structural framework common to many peptides from Conus venoms including the omega-, delta-, kappa-, and muO-conotoxins. However, TVIIA does not belong to these well-characterized pharmacological classes of conotoxins, but displays high sequence identity with conotoxin GS, a muscle sodium channel blocker from Conus geographus. Structure calculations were based on 562 interproton distance restraints inferred from NOE data, together with 18 backbone and nine side-chain torsion angle restraints derived from spin-spin coupling constants. The final family of 20 structures had mean pairwise rms differences over residues 2-27 of 0.18+/-0.05 A for the backbone atoms and 1.39+/-0.33 A for all heavy atoms. The structure consists of a triple-stranded, antiparallel beta sheet with +2x, -1 topology (residues 7-9, 16-20 and 23-27) and several beta turns. The core of the molecule is formed by three disulfide bonds which form a cystine knot motif common to many toxic and inhibitory polypeptides. The global fold, molecular shape and distribution of amino-acid sidechains in TVIIA is similar to that previously reported for conotoxin GS, and comparison with other four-loop conotoxin structures provides further indication that TVIIA and GS represent a new and distinct subgroup of this structural family. The structure of TVIIA determined in this study provides the basis for determining a structure-activity relationship for these molecules and their interaction with target receptors.     

X-ray crystal structures of the oxidized and reduced forms of the rubredoxin from the marine hyperthermophilic archaebacterium Pyrococcus furiosus.

The structures of the oxidized and reduced forms of the rubredoxin from the archaebacterium, Pyrococcus furiosus, an organism that grows optimally at 100 degrees C, have been determined by X-ray crystallography to a resolution of 1.8 A. Crystals of this rubredoxin grow in space group P2(1)2(1)2(1) with room temperature cell dimensions a = 34.6 A, b = 35.5 A, and c = 44.4 A. Initial phases were determined by the method of molecular replacement using the oxidized form of the rubredoxin from the mesophilic eubacterium, Clostridium pasteurianum, as a starting model. The oxidized and reduced models of P. furiosus rubredoxin each contain 414 nonhydrogen protein atoms comprising 53 residues. The model of the oxidized form contains 61 solvent H2O oxygen atoms and has been refined with X-PLOR and TNT to a final R = 0.178 with root mean square (rms) deviations from ideality in bond distances and bond angles of 0.014 A and 2.06 degrees, respectively. The model of the reduced form contains 37 solvent H2O oxygen atoms and has been refined to R = 0.193 with rms deviations from ideality in bond lengths of 0.012 A and in bond angles of 1.95 degrees. The overall structure of P. furiosus rubredoxin is similar to the structures of mesophilic rubredoxins, with the exception of a more extensive hydrogen-bonding network in the beta-sheet region and multiple electrostatic interactions (salt bridge, hydrogen bonds) of the Glu 14 side chain with groups on three other residues (the amino-terminal nitrogen of Ala 1; the indole nitrogen of Trp 3; and the amide nitrogen group of Phe 29). The influence of these and other features upon the thermostability of the P. furiosus protein is discussed.     

Solution structure of the PDZ2 domain from human phosphatase hPTP1E and its interactions with C-terminal peptides from the Fas receptor

The solution structure of the second PDZ domain (PDZ2) from human phosphatase hPTP1E has been determined using 2D and 3D heteronuclear NMR experiments. The binding of peptides derived from the C-terminus of the Fas receptor to PDZ2 was studied via changes in backbone peptide and protein resonances. The structure is based on a total of 1387 nonredundant experimental NMR restraints including 1261 interproton distance restraints, 45 backbone hydrogen bonds, and 81 torsion angle restraints. Analysis of 30 lowest-energy structures resulted in rmsd values of 0.41 +/- 0.09 A for backbone atoms (N, Calpha, C') and 1.08 +/- 0.10 A for all heavy atoms, excluding the disordered N- and C-termini. The hPTP1E PDZ2 structure is similar to known PDZ domain structures but contains two unique structural features. In the peptide binding domain, the first glycine of the GLGF motif is replaced by a serine. This serine appears to replace a bound water observed in PDZ crystal structures that hydrogen bonds to the bound peptide's C-terminus. The hPTP1E PDZ2 structure also contains an unusually large loop following strand beta2 and proximal to the peptide binding site. This well-ordered loop folds back against the PDZ domain and contains several residues that undergo large amide chemical shift changes upon peptide binding. Direct observation of peptide resonances demonstrates that as many as six Fas peptide residues interact with the PDZ2 domain.     

The three-dimensional structure of guanine-specific ribonuclease F1 in solution determined by NMR spectroscopy and distance geometry.

Two-dimensional 1H-NMR studies have been performed on ribonuclease F1 (RNase F1), which contains 106 amino acid residues. Sequence-specific resonance assignments were accomplished for the backbone protons of 99 amino acid residues and for most of their side-chain protons. The three-dimensional structures were constructed on the basis of 820 interproton-distance restraints derived from NOE, 64 distance restraints for 32 hydrogen bonds and 33 phi torsion-angle restraints. A total of 40 structures were obtained by distance geometry and simulated-annealing calculations. The average root-mean-square deviation (residues 1-106) between the 40 converged structures and the mean structure obtained by averaging their coordinates was 0.116 +/- 0.018 nm for the backbone atoms and 0.182 +/- 0.015 nm for all atoms including the hydrogen atoms. RNase F1 was determined to be an alpha/beta-type protein. A well-defined structure constitutes the core region, which consists of a small N-terminal beta-sheet (beta 1, beta 2) and a central five-stranded beta-sheet (beta 3-beta 7) packed on a long helix. The structure of RNase F1 has been compared with that of RNase T1, which was determined by X-ray crystallography. Both belong to the same family of microbial ribonucleases. The polypeptide backbone fold of RNase F1 is basically identical to that of RNase T1. The conformation-dependent chemical shifts of the C alpha protons are well conserved between RNase F1 and RNase T1. The residues implicated in catalysis are all located on the central beta-sheet in a geometry similar to that of RNase T1.     

A proton nuclear magnetic resonance study of the conformation of bovine anaphylatoxin C5a in solution

The solution conformation of bovine anaphylatoxin C5a has been investigated by nuclear magnetic resonance (NMR) spectroscopy. The 1H-NMR spectrum is assigned in a sequential manner using a variety of two-dimensional NMR techniques. A qualitative interpretation of the short range nuclear Overhauser enhancement data involving the NH, C alpha H and C beta H protons suggests that C5a has four helices comprising residues 5-11, 15-25, 33-39 and 46-61, and is composed of a globular head (residues 5-61) and a C-terminal tail. The polypeptide fold was determined by hybrid distance geometry-dynamical simulated annealing calculations on the basis of 203 approximate interproton distance restraints, 22 distance restraints for 11 intrahelical hydrogen bonds (identified on the basis of the pattern of short range NOEs and slowly exchanging backbone amide protons) and restraints for the 3 disulfide bridges. The overall polypeptide fold is similar to that of the sequence related human recombinant anaphylatoxin C5a [(1988) Proteins 3, 139-145].     

The solution conformation of the antibacterial peptide cecropin A: a nuclear magnetic resonance and dynamical simulated annealing study

The solution conformation of the antibacterial polypeptide cecropin A from the Cecropia moth is investigated by nuclear magnetic resonance (NMR) spectroscopy under conditions where it adopts a fully ordered structure, as judged by previous circular dichroism studies [Steiner, H. (1982) FEBS Lett. 137, 283-287], namely, 15% (v/v) hexafluoroisopropyl alcohol. By use of a combination of two-dimensional NMR techniques the 1H NMR spectrum of cecropin A is completely assigned. A set of 243 approximate interproton distance restraints is derived from nuclear Overhauser enhancement (NOE) measurements. These, together with 32 distance restraints for the 16 intrahelical hydrogen bonds identified on the basis of the pattern of short-range NOEs, form the basis of a three-dimensional structure determination by dynamical simulated annealing [Nilges, M., Clore, G.M., & Gronenborn, A.M. (1988) FEBS Lett. 229, 317-324]. The calculations are carried out starting from three initial structures, an alpha-helix, an extended beta-strand, and a mixed alpha/beta structure. Seven independent structures are computed from each starting structure by using different random number seeds for the assignments of the initial velocities. All 21 calculated structures satisfy the experimental restraints, display very small deviations from idealized covalent geometry, and possess good nonbonded contacts. Analysis of the 21 converged structure indicates that there are two helical regions extending from residues 5 to 21 and from residues 24 to 37 which are very well defined in terms of both atomic root mean square differences and backbone torsion angles. For the two helical regions individually the average backbone rms difference between all pairs of structures is approximately 1 A. The long axes of the two helices lie in two planes, which are at an angle of 70-100 degrees to each other. The orientation of the helices within these planes, however, cannot be determined due to the paucity of NOEs between the two helices.     

Strategy for supplementing structure calculations using limited data with hydrophobic distance restraints

Introductory biochemistry texts often note that the fold of a protein is completely defined when the dihedral angles phi and psi are known for each amino acid. This assertion was examined with torsion angle dynamics and simulated annealing (TAD/SA) calculations of protein G using only dihedral angle restraints. When all dihedral angles were restrained to within 1 degrees of the values of the X-ray structure, the TAD/SA structures gave a backbone root mean square deviation to the target of 4 A. Factors that contributed to divergence from the correct solution include deviations of peptide bonds from planarity, internal conflicts resulting from the nonuniform energies of different phi, psi combinations, and relaxation to extended conformations in the absence of long-range constraints. Simulations including hydrogen-bond restraints showed that even a few long-range contacts constrain the fold better than a complete set of accurate dihedral restraints. A procedure is described for TAD/SA calculations using hydrogen-bond restraints, idealized dihedral restraints for residues in regular secondary structures, and "hydrophobic distance restraints" derived from the positions of hydrophobic residues in the amino acid sequence. The hydrogen-bond restraints are treated as inviolable, whereas violated hydrophobic restraints are removed following reduction of restraint upper bounds from 2 to 1 times the predicted radius of gyration. The strategy was tested with simulated restraints from X-ray structures of proteins from different fold classes and NMR data for cold shock protein A that included only backbone chemical shifts and hydrogen bonds obtained from a long-range HNCO experiment.