Secondary structure of acylphosphatase from rabbit skeletal muscle. A nuclear magnetic resonance study

Sequence-specific assignment of 1H nuclear magnetic resonance spectra of acylphosphatase (EC 3.6.1.7) isolated from rabbit skeletal muscle have made it possible to identify short distance constraints from nuclear Overhauser enhancement spectra, to evaluate spin-spin coupling constants of many backbone amide hydrogens and to assess their slow exchange with deuterons in 2H2O solution. Analysis of these data show that the major regular secondary structure of the enzyme consists of five extended beta-strands, four of which are arranged in an antiparallel beta-sheet, while the fifth is attached parallel. A helix consisting of 11 residues has also been identified. Consideration of additional distance constraints between sequentially remote residues has allowed us to give an outline of the overall fold of the protein.     

A new method for building protein conformations from sequence alignments with homologues of known structure

We describe a largely automatic procedure for building protein structures from sequence alignments with homologues of known structure. This procedure uses simple rules by which multiple sequence alignments can be translated into distance and chirality constraints, which are then used as input for distance geometry calculations. By this means one obtains an ensemble of conformations for the unknown structure that are compatible with the rules employed, and the differences among these conformations provide an indication of the reliability of the structure prediction. The overall approach is demonstrated here by applying it to several Kazal-type trypsin inhibitors, for which experimentally determined structures are available. On the basis of our experience with these test problems, we have further predicted the conformation of the human pancreatic secretory trypsin inhibitor, for which no experimentally determined structure is presently available.     

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.     

Solution structure of [d(GCGTATACGC)]2.

The solution structure of the alternating pyrimidine-purine DNA duplex [d(GCGTATACGC)]2 has been determined using two-dimensional nuclear magnetic resonance techniques and distance geometry methods. Backbone distance constraints derived from experimental nuclear Overhauser enhancement and J-coupling torsion angle constraints were required to adequately define the conformation of the inter-residue backbone linkages and to avoid underwinding of the duplex. The distance geometry structures were further refined by back-calculation of the two-dimensional nuclear Overhauser enhancement spectra to correct spin-diffusion distance errors. Fifteen final structures for [d(GCGTATACGC)]2 were generated from the refined experimental distance bounds. These structures all exhibit fully wound B-form geometry with small penalty values (< 1.5 A) against the distance bounds and small pair-wise root-mean-square deviation values (typically 0.6 A to 1.5 A). The final structures exhibit positive base-pair inclination with respect to the helix axis, a marked alternation in rise and twist, and are shorter and wider than classical fiber B-form DNA. The purines were found to adopt a sugar pucker close to the C-2'-endo conformation while pyrimidine sugars exhibited significantly lower pseudorotation phase angles in the C-1'-exo to C-2'-endo range. The minor groove cross-strand steric clashes at pyrimidine-purine steps that would exist in pure B-DNA are attenuated by an increased rise at these steps (and an increased roll angle at TpA steps). Concomitantly the backbone torsion angles of the pyrimidine moieties have larger gamma values, larger epsilon values, and smaller zeta values than the purines. The structures generated by distance geometry methods were also compared with those obtained from restrained molecular dynamics with empirical force-field potentials. The results indicate that the nuclear magnetic resonance/distance geometry approach alone is capable of elucidating most of the salient structural features of double-stranded helical nucleic acids in solution without resorting to empirical energy potentials and without using any structural assumptions from crystallographic data.     

Solution structure of alpha t alpha, a helical hairpin peptide of de novo design.

alpha t alpha is a 38-residue peptide designed to adopt a helical hairpin conformation in solution (Fezoui Y, Weaver DL Osterhout JJ, 1995, Protein Sci 4:286-295). A previous study of the carboxylate form of alpha t alpha by CD and two-dimensional NMR indicated that the peptide was highly helical and that the helices associated in approximately the intended orientation (Fezoui Y, Weaver DL, Osterhout JJ, 1994, Proc Natl Acad Sci USA 91:3675-3679). Here, the solution structure of alpha t alpha as determined by two-dimensional NMR is reported. A total of 266 experimentally derived distance restraints and 20 dihedral angle restraints derived from J-couplings were used. One-hundred initial structures were generated by distance geometry and refined by dynamical simulated annealing. Twenty-three of the lowest-energy structures consistent with the experimental restraints were analyzed. The results presented here show that alpha t alpha is comprised of two associating helices connected by a turn region.     

Mobility of secondary structure units of horse-muscle acylphosphatase. Relation to antigenicity

The antigenic properties of acylphosphatase are compared with its various sequential characteristics (hydrophobicity, chemical shift of the main-chain 1H-NMR resonances, numbers and intensities of the nuclear Overhauser enhancements, hydrogen-deuterium exchange and sequential arrangement of the secondary structure units). The discussion is based on the complete sequential assignment of the 1H-NMR spectrum and the knowledge of the three-dimensional fold of the protein obtained by NMR spectroscopy from distance geometry calculations. Regions with very different degrees of mobility can be distinguished. It is found that all major antigenic sites are located in the most mobile surface loops.     

Three-dimensional solution structure of apo-neocarzinostatin from Streptomyces carzinostaticus determined by NMR spectroscopy.

The three-dimensional solution structure of apo-neocarzinostatin has been resolved from nuclear magnetic resonance spectroscopy data. Up to 1034 constraints were used to generate an initial set of 45 structures using a distance geometry algorithm (DSPACE). From this set, ten structures were subjected to refinement by restrained energy minimization and molecular dynamics. The average atomic root mean square deviations between the final ten structures and the mean structure obtained by averaging their coordinates run from 0.085 nm for the best defined beta-sheet regions of the protein to 0.227 nm for the side chains of the most flexible loops. The solution structure of apo-neocarzinostatin is closely similar to that of the related proteins, macromomycin and actinoxanthin. It contains a seven-stranded antiparallel beta-barrel which forms, together with two external loops, a deep cavity that is the chromophore binding site. It is noteworthy that aromatic side chains extend into the binding cleft. They may be responsible for the stabilization of the holo-protein complex and for the chromophore specificity within the antitumoral family.     

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

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

Three-dimensional structure in solution of barwin, a protein from barley seed.

The solution structure of a 125-residue basic protein, barwin, has been determined using 1H nuclear magnetic resonance spectroscopy. This protein is closely related to domains in proteins encoded by wound-induced genes in plants. Analysis of the 1H nuclear Overhauser spectrum revealed the assignment of more than 1400 nuclear Overhauser effects. Twenty structures were calculated based on 676 nontrivial distance restraints, 152 torsion angle restraints (92 phi, 56 chi 1, and 4 omega for proline), and stereospecific assignments of 38 chiral centers, using distance geometry, simulated annealing, and restrained energy minimization. None of the distance restraints was violated by more than 0.5 A in any of the 20 structures, and none of the torsion angle restraints was violated by more than 1 degree in any of the structures. The RMS difference between the calculated and target interproton distance restraints is 0.033 A, and the average atomic RMS differences between the 20 structures and their geometric average are 1.23 A for backbone atoms and 1.73 A for all heavy atoms. The dominating structural feature of the protein is a well-defined four-stranded antiparallel beta-sheet, two parallel beta-sheets packed antiparallel to each other and four short alpha-helices. The binding site of barwin to the tetramer N-acetylglucosamine has been qualitatively investigated, and the dissociation constant of the complex has been determined using one-dimensional 1H nuclear magnetic resonance spectroscopy.     

Three-dimensional solution structure of Ca(2+)-loaded porcine calbindin D9k determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of native, intact porcine calbindin D9k has been determined by distance geometry and restrained molecular dynamics calculations using distance and dihedral angle constraints obtained from 1H NMR spectroscopy. The protein has a well-defined global fold consisting of four helices oriented in a pairwise antiparallel manner such that two pairs of helix-loop-helix motifs (EF-hands) are joined by a linker segment. The two EF-hands are further coupled through a short beta-type interaction between the two Ca(2+)-binding loops. Overall, the structure is very similar to that of the highly homologous native, minor A form of bovine calbindin D9k determined by X-ray crystallography [Szebenyi, D. M. E., & Moffat, K. (1986) J. Biol. Chem. 261, 8761-8776]. A model structure built from the bovine calbindin D9k crystal structure shows several deviations larger than 2 A from the experimental distance constraints for the porcine protein. These structural differences are efficiently removed by subjecting the model structure to the experimental distance and dihedral angle constraints in a restrained molecular dynamics protocol, thereby generating a model that is very similar to the refined distance geometry derived structures. The N-terminal residues of the intact protein that are absent in the minor A form appear to be highly flexible and do not influence the structure of other regions of the protein. This result is important because it validates the conclusions drawn from the wide range of studies that have been carried out on minor A forms rather than the intact calbindin D9k.     

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 the EcoRI DNA sequence: refinement of NMR-derived distance geometry structures by NOESY spectrum back-calculations

The solution structure of the self-complementary DNA duplex [d(CGCGAATTCGCG)]2, which contains the EcoRI restriction site sequence GAATTC at the center, has been studied by two-dimensional nuclear magnetic resonance spectroscopy. Time-dependent nuclear Overhauser effect spectra were used to obtain the initial cross-relaxation rates between 155 pairs of protons. These initial cross-relaxation rates were converted into interproton distances and entered into a distance (bounds) matrix. A distance geometry algorithm (DSPACE) was used to create embedded starting structures and to refine these structures until they showed good agreement with the distance matrix; symmetry constraints were included in the refinement procedure, making the two strands in the refined distance geometry structures virtually identical and significantly improving the agreement with the distance matrix. The NOESY spectrum for one of these distance geometry structures was then calculated from the explicit coordinates by numerically integrating all the z-magnetization transfer pathways among neighboring protons within a specified radius. Distances in this distance geometry structure that did not agree with the experimental NOESY time course were then adjusted accordingly. This process was iterated until a good agreement between calculated and experimental NOESY spectra was reached. The final structure, which generates good agreement with the experimental NOESY spectrum, displays kinks at the C3-G4 base step and at the A6-T7 base step that appear to be similar to those reported for the EcoRI restriction site DNA bound to its endonuclease. The solution structure is not the same as the crystal structure of this DNA duplex.     

Alpha-cobratoxin: proton NMR assignments and solution structure.

The solution structure of alpha-cobratoxin, a neurotoxin purified from the venom of the snake Naja naja siamensis, at pH 3.2 is reported. Sequence-specific assignments of the NMR resonances was attained by a combination of a generalized main-chain-directed strategy and of the sequential method. The NMR data show the presence of a triple-stranded beta-sheet (residues 19-25, 36-41, and 52-57), a short helix, and turns. An extensive number of NOE cross peaks were identified in the NOESY NMR maps. These were applied as distance constraints in a molecular modeling protocol which includes distance geometry and dynamical simulated annealing calculations. A single family of structures is observed which fold in such a way that three major loops emerge from a globular head. The solution and crystal structures of alpha-cobratoxin are very similar. This is in clear contrast to results reported for alpha-bungarotoxin where significant differences exist.     

Design, synthesis and solution structure of a renin inhibitor. Structural constraints from NOE, and homonuclear and heteronuclear coupling constants combined with distance geometry calculations.

A macrocyclic renin inhibitor was designed using molecular modeling and a model of human renin. The synthesized molecular displayed poor binding affinity. To investigate the reasons for the observed inactivity, the structure of the compound has been studied by NMR spectroscopy and distance geometry. Structural constraints for distance geometry calculations were derived from nuclear Overhauser effects and homonuclear and heteronuclear three bond coupling constants. Homonuclear coupling constants were measured directly from the resolution-enhanced proton spectra and heteronuclear coupling constants were measured from the natural abundance 15N- and 13C-edited TOCSY experiments. One phi angle was determined uniquely by this method and two were reduced to two possible values each. By using a statistical analysis of 400 structures generated with distance geometry, two families of structures were found to be consistent with the NMR data. The solution structures so derived were different from the originally designed structure, including an internal hydrogen bond. This provides a possible explanation for the lack of effectiveness of this compound.     

Solution structure of neurotoxin I from the sea anemone Stichodactyla helianthus. A nuclear magnetic resonance, distance geometry, and restrained molecular dynamics study

The three-dimensional structure of the sea anemone polypeptide Stichodactyla helianthus neurotoxin I in aqueous solution has been determined using distance geometry and restrained molecular dynamics simulations based on NMR data acquired at 500 MHz. A set of 470 nuclear Overhauser enhancement values was measured, of which 216 were used as distance restraints in the structure determination along with 15 dihedral angles derived from coupling constants. After restrained molecular dynamics refinement, the eight structures that best fit the input data form a closely related family. They describe a structure that consists of a core of twisted, four-stranded, antiparallel beta-sheet encompassing residues 1-3, 19-24, 29-34, and 40-47, joined by three loops, two of which are well defined by the NMR data. The third loop, encompassing residues 7-16, is poorly defined by the data and is assumed to undergo conformational averaging in solution. Pairwise root mean square displacement values for the backbone heavy atoms of the eight best structures are 1.3 +/- 0.2A when the poorly defined loop is excluded and 3.6 +/- 1.0A for all backbone atoms. Refinement using restrained molecular dynamics improved the quality of the structures generated by distance geometry calculations with respect to the number of nuclear Overhauser enhancements violated, the size of the total distance violations and the total potential energies of the structures. The family of structures for S. heliathus neurotoxin I is compared with structures of related sea anemone proteins that also bind to the voltage-gated sodium channel.     

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.     

Distance geometry generates native-like folds for small helical proteins using the consensus distances of predicted protein structures.

For successful ab initio protein structure prediction, a method is needed to identify native-like structures from a set containing both native and non-native protein-like conformations. In this regard, the use of distance geometry has shown promise when accurate inter-residue distances are available. We describe a method by which distance geometry restraints are culled from sets of 500 protein-like conformations for four small helical proteins generated by the method of Simons et al. (1997). A consensus-based approach was applied in which every inter-Calpha distance was measured, and the most frequently occurring distances were used as input restraints for distance geometry. For each protein, a structure with lower coordinate root-mean-square (RMS) error than the mean of the original set was constructed; in three cases the topology of the fold resembled that of the native protein. When the fold sets were filtered for the best scoring conformations with respect to an all-atom knowledge-based scoring function, the remaining subset of 50 structures yielded restraints of higher accuracy. A second round of distance geometry using these restraints resulted in an average coordinate RMS error of 4.38 A.     

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.     

Tertiary structure of human complement component C5a in solution from nuclear magnetic resonance data

The tertiary structure for the region 1-63 of the 74 amino acid human complement protein C5a in solution was calculated from a large number of distance constraints derived from nuclear Overhauser effects with an angular distance geometry algorithm. The protein consists of four helices juxtaposed in an approximately antiparallel topology connected by peptide loops located at the surface of the molecule. The structures obtained for the helices are compatible with alpha-helical hydrogen-bonding patterns, which provides an explanation for the observed slow solvent exchange kinetics of the amide protons in these peptide regions. In contrast to the peptide region 1-63, no defined structure could be assigned to the C-terminal region 64-74, which increasingly acquires dynamic random coil characteristics as the end of the peptide chain is approached. An average root-mean-square deviation of 1.6 A was obtained for the alpha-carbons of the first 63 residues in the calculated ensemble of C5a structures, while the alpha-helices were determined with an average root-mean-square deviation of 0.8 A for the alpha-carbons. A comparison between the solution structure of C5a and the crystal structure of the functionally related C3a protein, as well as inferences for the interaction of C5a with its receptor on polymorphonuclear leukocytes, is discussed.     

Three-dimensional solution structure of the HIV-1 protease complexed with DMP323, a novel cyclic urea-type inhibitor, determined by nuclear magnetic resonance spectroscopy.

The three-dimensional solution structure of the HIV-1 protease homodimer, MW 22.2 kDa, complexed to a potent, cyclic urea-based inhibitor, DMP323, is reported. This is the first solution structure of an HIV protease/inhibitor complex that has been elucidated. Multidimensional heteronuclear NMR spectra were used to assemble more than 4,200 distance and angle constraints. Using the constraints, together with a hybrid distance geometry/simulated annealing protocol, an ensemble of 28 NMR structures was calculated having no distance or angle violations greater than 0.3 A or 5 degrees, respectively. Neglecting residues in disordered loops, the RMS deviation (RMSD) for backbone atoms in the family of structures was 0.60 A relative to the average structure. The individual NMR structures had excellent covalent geometry and stereochemistry, as did the restrained minimized average structure. The latter structure is similar to the 1.8-A X-ray structure of the protease/DMP323 complex (Chang CH et al., 1995, Protein Science, submitted); the pairwise backbone RMSD calculated for the two structures is 1.22 A. As expected, the mismatch between the structures is greatest in the loops that are disordered and/or flexible. The flexibility of residues 37-42 and 50-51 may be important in facilitating substrate binding and product release, because these residues make up the respective hinges and tips of the protease flaps. Flexibility of residues 4-8 may play a role in protease regulation by facilitating autolysis.