Conformation in solution of porcine brain natriuretic peptide determined by combined use of nuclear magnetic resonance and distance geometry

The conformation in solution of porcine brain natriuretic peptide was determined by combined use of NMR spectroscopy and distance geometry. A set of 157 inter-proton-distance constraints was derived from the two-dimensional NOE spectra, and further a set of three hydrogen bond constraints was obtained from analysis of the temperature dependence of labile protons. The five structures with minimal violations were selected after performing distance-geometry calculations starting from 40 random initial conformations. The distance-geometry structures were further refined by the use of restrained energy minimization and restrained molecular dynamics. This structure shows a compact conformation with the carboxy-terminal region, Asn21-Tyr26, folded back to the disulfide-linked loop region, Cys4-Cys20. The characteristics of the conformation determined are as follows: conformations of the three segments interposed by glycine residues, which are Arg7-Ile12, Ser14-Leu18 and Cys20-Arg25, were well defined and the segments Arg7-Ile12 and Cys20-Arg25 are rather close to each other and nearly parallel. The biological significance of these local conformations is discussed on the basis of comparisons with those of atrial natriuretic peptide reported by Kobayashi et al.     

Three-dimensional structure of acyl carrier protein in solution determined by nuclear magnetic resonance and the combined use of dynamical simulated annealing and distance geometry

The solution conformation of acyl carrier protein from Escherichia coli (77 residues) has been determined on the basis of 423 interproton-distance restraints and 32 hydrogen-bonding restraints derived from NMR measurements. A total of nine structures were computed using a hybrid approach combining metric matrix distance geometry and dynamic simulated annealing. The polypeptide fold is well defined with an average backbone atomic root-mean-square difference of 0.20 +/- 0.03 nm between the final nine converged structures and the mean structure obtained by averaging their coordinates. The principal structural motif is composed of three helices: 1 (residues 3-12), 2 (residues 37-47) and 4 (residues 65-75) which line a hydrophobic cavity. Helices 2 and 4 are approximately parallel to each other and anti-parallel at an angle of approximately equal to 150 degrees to helix 1. The smaller helix 3 (residues 56-63) is at an angle of approximately equal to 100 degrees to helix 4.     

The conformations of hirudin in solution: a study using nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The solution conformations of the protein hirudin have been investigated by the combined use of distance geometry and restained molecular dynamics calculations. The basis for the structure determination comprised 359 approximate inter-proton distance restrains and 10 phi backbone torsion angle restrains derived from n.m.r. measurements. It is shown that hirudin is composed of three domains: a central core made up of residues 3-30, 37-46 and 56-57; a protruding 'finger' (residues 31-36) consisting of the tip of an antiparallel beta sheet, and an exposed loop (residues 47-55). The structure of each individual domain is relatively well defined with average backbone atomic r.m.s. differences of <2 A between the final seven converged restrained dynamic structures and the mean structure obtained by averaging their coordinates. The orientation of the two minor domains relative to the central core, however, could not be determined as no long-range (i-h >5) interdomain proton-proton contacts could be observed in the two-dimensional nuclear Overhauser enhancement spectra. From the restrained molecular dynamics calculations it appears that the two minor domains exhibit large rigid-body motions relative to the central core.     

The three-dimensional structure of alpha1-purothionin in solution: combined use of nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The determination of the three-dimensional solution structure of α1-purothionin using a combination of metric matrix distance geometry and restrained molecular dynamics calculations based on n.m.r. data is presented. The experimental data comprise complete sequence-specific proton resonance assignments, a set of 310 approximate interproton distance restraints derived from nuclear Overhauser effects, 27 Ø backbone torsion angle restraints derived from vicinal coupling constants, 4 distance restraints from hydrogen bonds and 12 distance restraints from disulphide bridges. The average atomic rms difference between the final nine converged structures and the mean structure obtained by averaging their coordinates is 1.5 ± 0.1 å for the backbone atoms and 2.0 ± 0.1 å for all atoms. The overall shape of α1-purothionin is that of the capital letter L, similar to that of crambin, with the longer arm comprising two approximately parallel α-helices and the shorter arm a strand and a mini anti-parallel β sheet.     

Pseudo-structures for the 20 common amino acids for use in studies of protein conformations by measurements of intramolecular proton-proton distance constraints with nuclear magnetic resonance

"Pseudo-structures" of the 20 common amino acid residues are introduced for use in protein spatial structure determinations, which rely on the use of intramolecular proton-proton distance constraints determined by nuclear Overhauser effects as input for distance geometry calculations. The proposed structures satisfy requirements for the initial structural interpretation of the nuclear magnetic resonance data that arise from the absence of stereospecific assignments and/or limited spectral resolution for certain resonance lines. The pseudo-atoms used as reference points for the experimental distance constraints can be used in conjunction with the real amino acid structures representing the van der Waals' constraints on the spatial molecular structure, or with simplified models in order to reduce the computing time for the distance geometry calculations.     

Solution conformation of endothelin determined by nuclear magnetic resonance and distance geometry

The solution conformation of endothelium-derived vasoconstrictor peptide, endothelin, has been determined by two-dimensional 1H-NMR spectroscopy and distance geometry. Conformation in the N-terminal core region (residues 1-15) is well-defined and a characteristic is the helix-like conformation in the segment from Lys9 to Cys15. Contrarily, the C-terminal tail region (residues 16-21) does not assume a defined conformation and there are no specific interactions between the core and the tail regions.     

Solution structure of apamin determined by nuclear magnetic resonance and distance geometry

The solution structure of the bee venom neurotoxin apamin has been determined with a distance geometry program using distance constraints derived from NMR. Twenty embedded structures were generated and refined by using the program DSPACE. After error minimization using both conjugate gradient and dynamics algorithms, six structures had very low residual error. Comparisons of these show that the backbone of the peptide is quite well-defined with the largest rms difference between backbone atoms in these structures of 1.34 A. The side chains have far fewer constraints and show greater variability in their positions. The structure derived here is generally consistent with the qualitative model previously described, with most differences occurring in the loop between the beta-turn (residues 2-5) and the C-terminal alpha-helix (residues 9-17). Comparisons are made with previously derived models from NMR data and other methods.     

Solution structure of tertiapin determined using nuclear magnetic resonance and distance geometry

The solution structure of tertiapin, a 21-residue bee venom peptide, has been characterized by circular dichroism (CD), two-dimensional nuclear magnetic resonance (NMR) spectroscopy, and distance geometry. A total of 21 lowest error structures were obtained from distance geometry calculations. Superimposition of these structures shows that the backbone of tertiapin is very well defined. One type-I reverse turn from residue 4 to 7 and an α-helix from residue 12 to 19 exist in the structure of tertiapin. The α-helical region is best defined from both conformational analysis and structural superimposition. The overall three-dimensional structure of tertiapin is highly compact resulting from side chain interactions. The structural information obtained from CD and NMR are compared for both tertiapin and apamin (ref. 3), another bee venom peptide. Tertiapin and apamin have some similar secondary structure, but display different tertiary structures. © 1993 Wiley-Liss, Inc.     

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.     

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.     

Three-dimensional structure of the wild-type lac Pribnow promoter DNA in solution. Two-dimensional nuclear magnetic resonance studies and distance geometry calculations

The solution structure of a 12 base-pair DNA duplex containing the wt-lac promoter Pribnow sequence TATGTT has been studied by two-dimensional nuclear magnetic resonance spectroscopy. Proton assignments for the 24 sugar and base residues were obtained from two-dimensional correlated nuclear magnetic resonance and two-dimensional nuclear Overhauser effect spectra in both 2H2O and H2O, and by two-dimensional relayed coherence transfer nuclear magnetic resonance spectroscopy experiments. Time-dependent, two-dimensional nuclear Overhauser effect spectra were used to determine the initial cross-relaxation rates between 212 pairs of assigned protons, leading to 212 interproton distances in the double helix (8 to 9 per nucleotide). These distance constraints, and known bond lengths and angles, were entered into a distance matrix. After smoothing the bounds of the distance matrix, 12 trial matrices within the bounds constraints were independently generated and embedded in three-dimensional space using a distance geometry algorithm, to generate 12 trial structures. These trial structures were then refined until they no longer violated the distance matrix. The resulting structures are very similar at the local base-pair and nearest-neighbor base-pair level, but exhibit increasing variation at more distant and global levels. At the nearest-neighbor level, the A to T step and the G to T step within the Pribnow hexamer, as well as the G to T step preceding the hexamer, all exhibit very low screw pitch, i.e. 5(+/- 6) degrees. Conversely, the T to G step in the center of the promoter has a large screw pitch (47(+/- 2) degrees) and the T to G step at the 3' end of the promoter has a very large screw pitch (60(+/- 3) degrees). The limitations of nuclear magnetic resonance spectroscopy distance determination of structure are discussed in terms of resolution and spectral overlap of two-dimensional nuclear Overhauser effect crosspeaks. In the present duplex, the inability to measure several 1'-2' and 1'-2" distances resulted in underdetermination of the precise local sugar conformation for seven of the 24 residues, although the spatial position of all sugars was well defined.     

Application of nuclear magnetic resonance for the determination of the structure of proteins in solution

Knowledge of three-dimensional structure is a key factor in protein engineering. It is useful, for example, in predicting and understanding the functional consequences of specific substitution of one or more amino acids of the polypeptide chain. It is also necessary for the design of new effectors or analogs of the substrates of enzymes and receptors. X-ray diffraction by crystals of the biomolecule was for a long time the only method of determining three-dimensional structures. In the last 5 years, it has been joined by a new technique, two-dimensional nuclear magnetic resonance (2D NMR), which can resolve the structure of middle-sized proteins (less than 10 kilodaltons). The technique is applied on solutions whose pH, ionic strength, and temperature can be chosen and changed. The two basic measurements, COSY and NOESY, detect respectively the systems of hydrogen nuclei, or protons, coupled through covalent bonds, and those in which the interproton distances are less than 0.5 nm. A systematic strategy leads from resonance assignments of the two-dimensional spectrum to molecular modeling with constraints and finally to the determination of the molecular structure in the solution. Much sophistication is needed even today for the first task, the assignment of the resonances. Each of the COSY and NOESY spectra is a two-dimensional map, where the diagonal line is the one-dimensional spectrum, and the off-diagonal peaks indicate connectives between protons. Peak assignment to a specific type of amino acid is based on the pattern of scalar couplings observed in the COSY spectrum. Next, the amino acids are positioned in the primary sequence, using the spatial proximities of polypeptide chain protons, as observed in the NOESY spectrum. The principal secondary structures (alpha helix, beta sheets, etc.) are then identified by their specific connectivities. The tertiary structure is detected by NOESY connectivities between protons of different amino acids which are far apart in the primary sequence. The distance constraints from the NOESY connectivities also provide the starting point for modeling the tertiary structure. This is then refined using distance geometry and molecular dynamics algorithms. The resolution of the structures obtained with the help of recent algorithmic developments may be comparable to that provided by X-ray diffraction. The COSY measurement can be completed or substituted by other measurements, useful albeit more complex. For example, the HOHAHA experiment, currently in wide use, gives the correlations through multiple covalent bonds. Multiquanta experiments, which select systems of a given number of coupled spins, provide spectral simplification.(ABSTRACT TRUNCATED AT 400 WORDS)     

Structure determination of biological macromolecules in solution using nuclear magnetic resonance spectroscopy

A detailed understanding of the function of a biological macromolecule requires knowledge of its three-dimensional structure. Most atomic-resolution structures of biological macromolecules have been solved either by X-ray diffraction in single crystals or by nuclear magnetic resonance (NMR) in solution. This review surveys the method of NMR structure determination. First, a brief introduction to NMR and its basic concepts is presented. The main part of the article deals with the individual steps necessary for an NMR structure determination. At the end, the discussion turns to considerations on the influence of the molecular size of the macromolecules on the structure determination by NMR. New techniques are discussed that greatly enhance the possibilities of applying NMR to large molecular systems.     

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.     

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.     

Solution conformation of conotoxin GI determined by 1H nuclear magnetic resonance spectroscopy and distance geometry calculations

Conformational analysis of conotoxin GI, one of the neurotoxic peptides produced by a marine snail, genus Conus, was performed by a combination of nuclear magnetic resonance spectroscopy (NMR) and distance geometry calculations. The resulting conformers on minimization of the target function were classified into two groups. The difference in the structures of the conformers is mainly due to the difference in the orientation of the side chain of the tyrosyl residue. The results show that the solution structure of conotoxin GI satisfies the conformational requirements for the biological activity of an antagonist toward nicotinic cholinergic receptors elucidated in a series of studies on alkaloids. The structure is discussed on the basis of the results of comparison of the atomic arrangements of the active sites of snake venom peptides and molecular models based on the results of secondary structure prediction.     

Quantitative determination of conformations of flexible molecules in solution using lanthanide ions as nuclear magnetic resonance probes: application to adenosine-5'-monophosphate

A procedure is described here whereby the conformation, of a flexible molecule in solution can be found. The method depends on the study of the nuclear magnetic resonance spectrum of the molecule in the presence of perturbations due to specifically bound lanthanide cations. The magnetic perturbations are of two kinds: shifts of nuclear magnetic resonance spectral lines in the presence of cations such as Eu³⁺ and changes in relaxation rates of the nuclear magnetic resonance excitations in the presence of cations such as Gd³⁺. Suitable expressions are given for the relation between the magnitude of the perturbations and the geometry of the lanthanide complex in the absence of through-bond perturbations and for an axially symmetric system. It is proved that the spectral changes described here are not due to through-bond (contact) effects. The circumstances, in which the anisotropy of the magnetic susceptibility tensor, as seen in the nuclear magnetic resonance spectra, is of axial symmetry, are defined. The experimental systems described are of this kind. A computer program has been devised that searches for the conformations of the molecule which fit the nuclear magnetic resonance data.We outline here the principles of the method and how we have used a combination of relaxation and shift probes to obtain the conformation of adenosine-5′-monophosphate at pH 2. It is shown that a small family of closely related conformations fit the nuclear magnetic resonance data. These conformations are very similar to that of the crystal structure of AMP.     

Combined use of proton-proton Overhauser enhancements and a distance geometry algorithm for determination of polypeptide conformations. Application to micelle-bound glucagon

In a new approach for the determination of polypeptide conformation, experimental data on intramolecular distances between pairs of hydrogen atoms obtained from nuclear Overhauser enhancement studies are used as input for a distance geometry algorithm. The algorithm determines the limits of the conformation space occupied by the polypeptide chain. The experimental data are used in such a way that the real conformation should in all cases be within these limits. Two important features of the method are that the results do not depend critically on the accuracy of the distance measurements by nuclear Overhauser enhancement studies and that internal mobility of the polypeptide conformation is explicitly taken into consideration. The use of this new procedure is illustrated with a structural study of the region 19-27 of glucagon bound to perdeuterated dodecylphosphocholine micelles.