NMR-pseudoenergy approach to the solution structure of acyl carrier protein

A method for protein structure determination from two-dimensional NMR cross-relaxation data is presented and explored by using short amino acid segments from acyl carrier protein as a test case. The method is based on a molecular mechanics program and incorporates NMR distance constraints in the form of a pseudoenergy term that accurately reflects the distance-dependent precision of NMR cross-relaxation data. When it is used in an indiscriminant fashion, the method has a tendency to produce structures representing local energy minima near starting structures, rather than structures representing a global energy minimum. However, stepwise inclusion of energy terms, beginning with a function heavily weighted by backbone distance constraints, appears to simplify the potential energy surface to a point where convergence to a common backbone structure from a variety of starting structures is possible. In the case of the segment from residues 3 to 15 in acyl carrier protein, a nearly perfect alpha-helix is produced starting with a linear chain, an alpha-helical chain, or a chain having residues with alternating linear and alpha-helical backbone torsional angles. In the case of the segment from residues 26 to 36 a structure having a right-handed loop is produced.     

Solution structure of murine epidermal growth factor determined by NMR spectroscopy and refined by energy minimization with restraints.

The solution structure of murine epidermal growth factor (mEGF) at pH 3.1 and a temperature of 28 degrees C has been determined from NMR data, using distance geometry calculations and restrained energy minimization. The structure determination is based on 730 conformational constraints derived from NMR data, including 644 NOE-derived upper bound distance constraints, constraints on the ranges of 32 dihedral angles based on measurements of vicinal coupling constants, and 54 upper and lower bound constraints associated with nine hydrogen bonds and the three disulfide bonds. The distance geometry interpretation of the NMR data is based on previously published sequence-specific 1H resonance assignments [Montelione et al. (1988) Biochemistry 27, 2235-2243], supplemented here with individual assignments for some side-chain amide, methylene, and isopropyl methyl protons. The molecular architecture of mEGF is the same as that described previously [Montelione et al. (1987) Proc. Natl. Acad. Sci. U.S.A. 84, 5226-5230], but the structure is overall more precisely determined by a more extensive set of NMR constraints. Analysis of proton NMR line widths, amide proton exchange rates, and side-chain 3J(H alpha-H beta) coupling constants provides evidence for internal motion in several regions of the mEGF molecule. Because mEGF is one member of a large family of homologous growth factors and protein domains for which X-ray crystal structures are not yet available, the atomic coordinates resulting from the present structure refinement (which we have deposited in the Brookhaven Protein Data Bank) are important data for understanding the structures of EGF-like proteins and for further detailed comparisons of these structures with mEGF.     

The war of tools: how can NMR spectroscopists detect errors in their structures?

Protein structure determination by NMR methods has started in the mid-eighties and has been growing steadily since then. Ca. 14% of the protein structures deposited in the PDB have been solved by NMR. The evaluation of the quality of NMR structures however is still lacking a well-established practice. In this work, we examined various tools for the assessment of structural quality to ascertain the extent to which these tools could be applied to detect flaws in NMR structures. In particular, we investigated the variation in the scores assigned by these programs as a function of the deviation of the structures induced by errors in assignments or in the upper distance limits used. These perturbations did not distort radically the protein fold, but resulted in backbone RMS deviations up to 3 Å, which is in line with errors highlighted in the available literature. We found that it is quite difficult to discriminate the structures perturbed because of misassignments from the original ones, also because the spread in score over the conformers of the original bundle is relatively large. ?–ψ distributions and normality scores related to the backbone conformation and to the distribution of side-chain dihedral angles are the most sensitive indicators of flaws.     

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.     

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.     

Automated protein fold determination using a minimal NMR constraint strategy

Determination of precise and accurate protein structures by NMR generally requires weeks or even months to acquire and interpret all the necessary NMR data. However, even medium-accuracy fold information can often provide key clues about protein evolution and biochemical function(s). In this article we describe a largely automatic strategy for rapid determination of medium-accuracy protein backbone structures. Our strategy derives from ideas originally introduced by other groups for determining medium-accuracy NMR structures of large proteins using deuterated, (13)C-, (15)N-enriched protein samples with selective protonation of side-chain methyl groups ((13)CH(3)). Data collection includes acquiring NMR spectra for automatically determining assignments of backbone and side-chain (15)N, H(N) resonances, and side-chain (13)CH(3) methyl resonances. These assignments are determined automatically by the program AutoAssign using backbone triple resonance NMR data, together with Spin System Type Assignment Constraints (STACs) derived from side-chain triple-resonance experiments. The program AutoStructure then derives conformational constraints using these chemical shifts, amide (1)H/(2)H exchange, nuclear Overhauser effect spectroscopy (NOESY), and residual dipolar coupling data. The total time required for collecting such NMR data can potentially be as short as a few days. Here we demonstrate an integrated set of NMR software which can process these NMR spectra, carry out resonance assignments, interpret NOESY data, and generate medium-accuracy structures within a few days. The feasibility of this combined data collection and analysis strategy starting from raw NMR time domain data was illustrated by automatic analysis of a medium accuracy structure of the Z domain of Staphylococcal protein A.     

Simulated annealing with restrained molecular dynamics using a flexible restraint potential: theory and evaluation with simulated NMR constraints.

A new functional representation of NMR-derived distance constraints, the flexible restraint potential, has been implemented in the program CONGEN (Bruccoleri RE, Karplus M, 1987, Biopolymers 26:137-168) for molecular structure generation. In addition, flat-bottomed restraint potentials for representing dihedral angle and vicinal scalar coupling constraints have been introduced into CONGEN. An effective simulated annealing (SA) protocol that combines both weight annealing and temperature annealing is described. Calculations have been performed using ideal simulated NMR constraints, in order to evaluate the use of restrained molecular dynamics (MD) with these target functions as implemented in CONGEN. In this benchmark study, internuclear distance, dihedral angle, and vicinal coupling constant constraints were calculated from the energy-minimized X-ray crystal structure of the 46-amino acid polypeptide crambin (ICRN). Three-dimensional structures of crambin that satisfy these simulated NMR constraints were generated using restrained MD and SA. Polypeptide structures with extended backbone and side-chain conformations were used as starting conformations. Dynamical annealing calculations using extended starting conformations and assignments of initial velocities taken randomly from a Maxwellian distribution were found to adequately sample the conformational space consistent with the constraints. These calculations also show that loosened internuclear constraints can allow molecules to overcome local minima in the search for a global minimum with respect to both the NMR-derived constraints and conformational energy. This protocol and the modified version of the CONGEN program described here are shown to be reliable and robust, and are applicable generally for protein structure determination by dynamical simulated annealing using NMR data.     

Solution NMR studies of the integral membrane proteins OmpX and OmpA from Escherichia coli

Membrane proteins are usually solubilized in polar solvents by incorporation into micelles. Even for small membrane proteins these mixed micelles have rather large molecular masses, typically beyond 50000 Da. The NMR technique TROSY (transverse relaxation-optimized spectroscopy) has been developed for studies of structures of this size in solution. In this paper, strategies for the use of TROSY-based NMR experiments with membrane proteins are discussed and illustrated with results obtained with the Escherichia coli integral membrane proteins OmpX and OmpA in mixed micelles with the detergent dihexanoylphosphatidylcholine (DHPC). For OmpX, complete sequence-specific NMR assignments have been obtained for the polypeptide backbone. The 13C chemical shifts and nuclear Overhauser effect data then resulted in the identification of the regular secondary structure elements of OmpX/DHPC in solution, and in the collection of an input of conformational constraints for the computation of the global fold of the protein. For OmpA, the NMR assignments are so far limited to about 80% of the polypeptide chain, indicating different dynamic properties of the reconstituted OmpA beta-barrel from those of OmpX. Overall, the present data demonstrate that relaxation-optimized NMR techniques open novel avenues for studies of structure, function and dynamics of integral membrane proteins.     

Solution conformation of thymosin beta 4: a nuclear magnetic resonance and simulated annealing study

The conformation of the polypeptide thymosin beta 4 in solutions of 60% (v/v) trifluoroethanol-d3 and 50% (v/v) hexafluoroisopropyl-d2 alcohol in water is investigated by nuclear magnetic resonance (NMR) spectroscopy. Under these conditions thymosin beta 4 adopts an ordered structure. By use of a combination of two-dimensional NMR techniques, the 1H NMR spectrum of thymosin beta 4 is assigned. A set of 180 approximate interproton distance constraints is derived from nuclear Overhauser enhancement (NOE) measurements. These, together with 33 phi constraints obtained for JNH alpha coupling data and the 23 psi dihedral angles identified on the basis of the pattern of short-range NOEs, form the basis of a three-dimensional structure determination by dynamical simulated annealing. The calculations are carried out starting from three initial structures, an alpha-helix, an extended beta-strand, and a mixed alpha/beta structure. Ten independent structures are computed from each starting structure by using different random number seeds for the assignments of the initial velocities. All 30 calculated structures satisfy the experimental constraints, display very small deviations from idealized covalent geometry, and possess good nonbonded contacts. Analysis of the 30 converged structures indicates that there are two helical regions extending from residues 4-16 and from residues 30-40, which are well defined both in terms of 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 2 A. The two helices exhibit typical amino acid preferences for specific locations at the ends of helices.(ABSTRACT TRUNCATED AT 250 WORDS)     

Determining the three-dimensional fold of a protein from approximate constraints

We propose a new approach for calculating the three-dimensional (3D) structure of a protein from distance and dihedral angle constraints derived from experimental data. We suggest that such constraints can be obtained from experiments such as tritium planigraphy, chemical or enzymatic cleavage of the polypeptide chain, paramagnetic perturbation of nuclear magnetic resonance (NMR) spectra, measurement of hydrogen-exchange rates, mutational studies, mass spectrometry, and electron paramagnetic resonance. These can be supplemented with constraints from theoretical prediction of secondary structures and of buried/exposed residues. We report here distance geometry calculations to generate the structures of a test protein Staphylococcal nuclease (STN), and the HIV-1 rev protein (REV) of unknown structure. From the available 3D atomic coordinates of STN, we set up simulated data sets consisting of varying number and quality of constraints, and used our group's Self Correcting Distance Geometry (SECODG) program DIAMOD to generate structures. We could generate the correct tertiary fold from qualitative (approximate) as well as precise distance constraints. The root mean square deviations of backbone atoms from the native structure were in the range of 2.0 A to 8.3 A, depending on the number of constraints used. We could also generate the correct fold starting from a subset of atoms that are on the surface and those that are buried. When we used data sets containing a small fraction of incorrect distance constraints, the SECODG technique was able to detect and correct them. In the case of REV, we used a combination of constraints obtained from mutagenic data and structure predictions. DIAMOD generated helix-loop-helix models, which, after four self-correcting cycles, populated one family exclusively. The features of the energy-minimized model are consistent with the available data on REV-RNA interaction. Our method could thus be an attractive alternative for calculating protein 3D structures, especially in cases where the traditional methods of X-ray crystallography and multidimensional NMR spectroscopy have been unsuccessful.     

Determination of the complete three-dimensional structure of the trypsin inhibitor from squash seeds in aqueous solution by nuclear magnetic resonance and a combination of distance geometry and dynamical simulated annealing

The complete three-dimensional structure of the trypsin inhibitor from seeds of the squash Cucurbita maxima in aqueous solution was determined on the basis of 324 interproton distance constraints, 80 non-nuclear Overhauser effect distances, and 22 hydrogen-bonding constraints, supplemented by 27 phi backbone angle constraints derived from nuclear magnetic resonance measurements. The nuclear magnetic resonance input data were converted to the distance constraints in a semiquantitative manner after a sequence specific assignment of 1H spectra was obtained using two-dimensional nuclear magnetic resonance techniques. Stereospecific assignments were obtained for 17 of the 48 prochiral centers of the squash trypsin inhibitor using the floating chirality assignment introduced at the dynamical simulated annealing stage of the calculations. A total of 34 structures calculated by a hybrid distance geometry-dynamical simulated annealing method exhibit well-defined positions for both backbone and side-chain atoms. The average atomic root-mean-square difference between the individual structures and the minimized mean structure is 0.35(+/- 0.08) A for the backbone atoms and 0.89(+/- 0.17) A for all heavy atoms. The precision of the structure determination is discussed and correlated to the experimental input data.     

Automated assignment and 3D structure calculations using combinations of 2D homonuclear and 3D heteronuclear NOESY spectra

The NOAH/DIAMOD suite uses feedback filtering and self-correcting distance geometry to generate 3D structures from unassigned NOESY spectra. In this study we determined the minimum set of experiments needed to generate a high quality structure bundle. Different combinations of 3D ¹⁵N-edited, ¹³C-edited HSQC-NOESY and 2D homonuclear ¹H-¹H NOESY spectra of the 77 amino acid protein, myeloid progenitor inhibitory factor-1 (MPIF-1) were used as input for NOAH/DIAMOD calculations. The quality of the assignments of NOESY cross peaks and the accuracy of the automatically generated 3D structures were compared to those obtained with a conventional manual procedure. Combining data from two types of experiments synergistically increased the number of peaks assigned unambiguously in both individual spectra. As a general trend for the accuracy of the structures we observed structural variations in the backbone fold of the final structures of about 2 Å for single spectral data, of 1 Å to 1.5 Å for double spectral data, and of 0.6 Å for triple spectral data sets. The quality of the assignments and 3D structures from the optimal data using all three spectra were similar to those obtained from traditional assignment methods with structural variations within the bundle of 0.6 Å and 1.3 Å for backbone and heavy atoms, respectively. Almost all constraints (97%) of the automatic NOESY cross peak assignments were cross compatible with the structures from the conventional manual assignment procedure, and an even larger proportion (99%) of the manually derived constraints were compatible with the automatically determined 3D structures. The two mean structures determined by both methods differed only by 1.3 Å rmsd for the backbone atoms in the well-defined regions of the protein. Thus NOAD/DIAMOD analysis of spectra from labeled proteins provides a reliable method for high throughput analysis of genomic targets.     

A computational method for NMR-constrained protein threading.

Protein threading provides an effective method for fold recognition and backbone structure prediction. But its application is currently limited due to its level of prediction accuracy and scope of applicability. One way to significantly improve its usefulness is through the incorporation of underconstrained (or partial) NMR data. It is well known that the NMR method for protein structure determination applies only to small proteins and that its effectiveness decreases rapidly as the protein mass increases beyond about 30 kD. We present, in this paper, a computational framework for applying underconstrained NMR data (that alone are insufficient for structure determination) as constraints in protein threading and also in all-atom model construction. In this study, we consider both secondary structure assignments from chemical shifts and NOE distance restraints. Our results have shown that both secondary structure assignments and a small number of long-range NOEs can significantly improve the threading quality in both fold recognition and threading-alignment accuracy, and can possibly extend threading's scope of applicability from homologs to analogs. An accurate backbone structure generated by NMR-constrained threading can then provide a great amount of structural information, equivalent to that provided by many NMR data; and hence can help reduce the number of NMR data typically required for an accurate structure determination. This new technique can potentially accelerate current NMR structure determination processes and possibly expand NMR's capability to larger proteins.     

1H NMR assignments of apo calcyclin and comparative structural analysis with calbindin D9k and S100 beta.

The homodimeric S100 protein calcyclin has been studied in the apo state by two-dimensional 1H NMR spectroscopy. Using a combination of scalar correlation and NOE experiments, sequence-specific 1H NMR assignments were obtained for all but one backbone and > 90% of the side-chain resonances. To our knowledge, the 2 x 90 residue (20 kDa) calcyclin dimer is the largest protein system for which such complete assignments have been made by purely homonuclear methods. Sequential and medium-range NOEs and slowly exchanging backbone amide protons identified directly the four helices and the short antiparallel beta-type interaction between the two binding loops that comprise each subunit of the dimer. Further analysis of NOEs enabled the unambiguous assignment of 556 intrasubunit distance constraints, 24 intrasubunit hydrogen bonding constraints, and 2 x 26 intersubunit distance constraints. The conformation of the monomer subunit was refined by distance geometry and restrained molecular dynamics calculations using the intrasubunit constraints only. Calculation of the dimer structure starting from this conformational ensemble has been reported elsewhere. The extent of structural homology among the apo calcyclin subunit, the monomer subunit of apo S100 beta, and monomeric apo calbindin D9k has been examined in detail by comparing 1H NMR chemical shifts and secondary structures. This analysis was extended to a comprehensive comparison of the three-dimensional structures of the calcyclin monomer subunit and calbindin D9k, which revealed greater similarity in the packing of their hydrophobic cores than was anticipated previously. Together, these results support the hypothesis that all members of the S100 family have similar core structures and similar modes of dimerization. Analysis of the amphiphilicity of Helix IV is used to explain why calbindin D9k is monomeric, but full-length S100 proteins form homodimers.     

The solution structure of a [3Fe-4S] ferredoxin: oxidised ferredoxin II from Desulfovibrio gigas

The use of standard 2D NMR experiments in combination with 1D NOE experiments allowed the assignment of 51 of the 58 spin systems of oxidised [3Fe-4S] ferredoxin isolated from Desulfovibrio gigas. The NMR solution structure was determined using data from 1D NOE and 2D NOESY spectra, as distance constraints, and information from the X-ray structure for the spin systems not detected by NMR in torsion angle dynamics calculations to produce a family of 15 low target function structures. The quality of the NMR family, as judged by the backbone r.m.s.d. values, was good (0.80?Å), with the majority of φ/ψ angles falling within the allowed region of the Ramachandran plot. A comparison with the X-ray structure indicated that the overall global fold is very similar in solution and in the solid state. The determination of the solution structure of ferredoxin II (FdII) in the oxidised state (FdII_ox) opens the way for the determination of the solution structure of the redox intermediate state of FdII (FdII_int), for which no X-ray structure is available.     

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.     

Determination of the three-dimensional solution structure of the C-terminal domain of cellobiohydrolase I from Trichoderma reesei. A study using nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing

The solution structure of a synthetic 36-residue polypeptide comprising the C-terminal cellulose binding domain of cellobiohydrolase I (CT-CBH I) from Trichoderma reesei was investigated by nuclear magnetic resonance (NMR) spectroscopy. The 1H NMR spectrum was completely assigned in a sequential manner by two-dimensional NMR techniques. A large number of stereospecific assignments for beta-methylene protons, as well as ranges for the phi, psi, and chi 1 torsion angles, were obtained on the basis of sequential and intraresidue nuclear Overhauser enhancement (NOE) and coupling constant data in combination with a conformational data base search. The structure calculations were carried out in an iterative manner by using the hybrid distance geometry-dynamical simulated annealing method. This involved computing a series of initial structures from a subset of the experimental data in order to resolve ambiguities in the assignments of some NOE cross-peaks arising from chemical shift degeneracy. Additionally, this permitted us to extend the stereospecific assignments to the alpha-methylene protons of glycine using information on phi torsion angles derived from the initial structure calculations. The final experimental data set consisted of 554 interproton distance restraints, 24 restraints for 12 hydrogen bonds, and 33 phi, 24 psi, and 25 chi 1 torsion angle restraints. CT-CBH I has two disulfide bridges whose pairing was previously unknown. Analysis of structures calculated with all three possible combinations of disulfide bonds, as well as without disulfide bonds, indicated that the correct disulfide bridge pairing was 8-25 and 19-35. Forty-one structures were computed with the 8-25 and 19-35 disulfide bridges, and the average atomic rms difference between the individual structures and the mean structure obtained by averaging their coordinates was 0.33 +/- 0.04 A for the backbone atoms and 0.52 +/- 0.06 A for all atoms. The protein has a wedgelike shape with an amphiphilic character, one face being predominantly hydrophilic and the other mainly hydrophobic. The principal element of secondary structure is made up of an irregular triple-stranded antiparallel beta-sheet composed of residues 5-9 (beta 1), 24-28 (beta 2), and 33-36 (beta 3) in which strand beta 3 is hydrogen bonded to the other two strands.(ABSTRACT TRUNCATED AT 400 WORDS)     

Rapid analysis of protein backbone resonance assignments using cryogenic probes,a distributed Linux-based computing architecture,and an integrated set of spectral analysis tools

Rapid data collection, spectral referencing, processing by time domain deconvolution, peak picking and editing, and assignment of NMR spectra are necessary components of any efficient integrated system for protein NMR structure analysis. We have developed a set of software tools designated AutoProc, AutoPeak, and AutoAssign, which function together with the data processing and peak-picking programs NMRPipe and Sparky, to provide an integrated software system for rapid analysis of protein backbone resonance assignments. In this paper we demonstrate that these tools, together with high-sensitivity triple resonance NMR cryoprobes for data collection and a Linux-based computer cluster architecture, can be combined to provide nearly complete backbone resonance assignments and secondary structures (based on chemical shift data) for a 59-residue protein in less than 30 hours of data collection and processing time. In this optimum case of a small protein providing excellent spectra, extensive backbone resonance assignments could also be obtained using less than 6 hours of data collection and processing time. These results demonstrate the feasibility of high throughput triple resonance NMR for determining resonance assignments and secondary structures of small proteins, and the potential for applying NMR in large scale structural proteomics projects.Abbreviations: BPTI – bovine pancreatic trypsin inhibitor; LP – linear prediction; FT – Fourier transform; S/N – signal-to-noise ratio; FID – free induction decay     

Determining stereo-specific 1H nuclear magnetic resonance assignments from distance geometry calculations

Stereo-specific 1H nuclear magnetic resonance assignments can be obtained following distance geometry structure calculations. The key to this method is to allow stereo-related atoms or methyls to float between pro-R and pro-S configurations, the final configuration being determined by the experimental constraints. Resonances from stereo-related pairs are given initial random assignments (either pro-R or pro-S) for identifying nuclear Overhauser effects (NOEs). A list of distance constraints using these assignments is compiled and a series of structures calculated where the chirality of non-C alpha chiral centers is not constrained; no pseudoatom corrections are required. Calculated structures are both locally and globally well-determined since the assignments rely upon the structure determination rather than the structure quality relying upon stereo-specific assignments. The method represents a global approach to determining stereo-specific assignments versus previously reported methods where only intraresidue NOEs and J-coupling information are used.     

Solution structures of proteins from NMR data and modeling: alternative folds for neutrophil peptide 5

The structure of neutrophil peptide 5 in solution has recently been reported (Pardi et al., 1988). The structure determination was accomplished by using a distance geometry algorithm and 107 interproton distance constraints obtained from 2D NMR data. In each of the eight independent solutions to the distance geometry equations, the overall fold of the polypeptide backbone was identical and the root mean square (rms) deviation between backbone atoms of the superimposed structures was small (approximately 2.4 A). In this paper we report additional NP-5 structures obtained by using a new structure generation algorithm: a Monte Carlo search in torsion angle space. These structures have a large rms backbone deviation from the distance geometry structures (approximately 5.0 A). The backbone topologies differ in significant respects from the distance geometry structures and from each other. Structures are found that are pseudo mirror images of part or all of the fold corresponding to that first obtained with the distance geometry procedure. For small proteins, the problem of distinguishing the correct structure among pseudo mirror images is likely to be greater than previously recognized. When a set of test distance constraints constructed from a novel Monte Carlo structure is used as input in the distance geometry algorithm, the fold of the resulting structure does not correspond to that of the target. The results also demonstrate that the previously accepted criteria (the magnitude of the rms deviation between multiple solutions of the distance geometry equations) for defining the accuracy and precision of a peptide structure generated from NMR data are inadequate. An energetic analysis of structures corresponding to the different folding topologies has been carried out. The molecular mechanics energies obtained by minimization and molecular dynamics refinement provide sufficient information to eliminate certain alternative structures. On the basis of a careful comparison of the different trial structures with the experimental data, it is concluded that the NP-5 peptide fold which was originally reported is most consistent with the data. An alternative fold corresponding to structures with low energies and small total distance violations is ruled out because for this fold predicted NOEs are not observed experimentally.