1H 2D NMR and distance geometry study of the folding of Ecballium elaterium trypsin inhibitor, a member of the squash inhibitors family

The solution conformation of synthetic Ecballium elaterium trypsin inhibitor II, a 28-residue peptide with 3 disulfide bridges, has been studied by 1H 2D NMR measurements. Secondary structure elements were determined: a miniantiparallel beta-sheet Met 7-Cys 9 and Gly 25-Cys 27, a beta-hairpin 20-28 with beta-turn 22-25, and two tight turns Asp 12-Cys 15 and Leu 16-Cys 19. A set of interproton distance restraints deduced from two-dimensional nuclear Overhauser enhancement spectra and 13 phi backbone torsion angles restraints were used as the basis of three-dimensional structure computations including disulfide bridges arrangement by using distance geometry calculations. Computations for the 15 possible S-S linkage combinations lead to the proposal of the array 2-19, 9-21, 15-27 as the most probable structure for EETI II.     

Effects of glycine substitutions on the structure and function of gramicidin a channels

Tryptophan residues often are found at the lipid-aqueous interface region of membrane-spanning proteins, including ion channels, where they are thought to be important determinants of protein structure and function. To better understand how Trp residues modulate the function of membrane-spanning channels, we have examined the effects of Trp replacements on the structure and function of gramicidin A channels. Analogues of gramicidin A in which the Trp residues at positions 9, 11, 13, and 15 were sequentially replaced with Gly were synthesized, and the three-dimensional structure of each analogue was determined using a combination of two-dimensional NMR techniques and distance geometry-simulated annealing structure calculations. Though Trp --> Gly substitutions destabilize the beta6.3-helical gA channel structure, it is possible to determine the structure of analogues with Trp --> Gly substitutions at positions 11, 13, and 15, but not for the analogue with the Trp --> Gly substitution at position 9. The Gly11-, Gly13-, and Gly15-gA analogues form channels that adopt a backbone fold identical to that of native gramicidin A, with only small changes in the side chain conformations of the unsubstituted residues. Single-channel current measurements show that the channel function and lifetime of the analogues are significantly affected by the Trp --> Gly replacements. The conductance variations appear to be caused by sequential removal of the Trp dipoles, which alter the ion-dipole interactions that modulate ion movement. The lifetime variations did not appear to follow a clear pattern.     

Conformation of recombinant desulfatohirudin in aqueous solution determined by nuclear magnetic resonance

The three-dimensional structure of recombinant desulfatohirudin in aqueous solution was determined by 1H nuclear magnetic resonance at 600 MHz and distance geometry calculations with the program DISMAN. The input for the structure calculations was prepared on the basis of complete sequence-specific resonance assignments at pH 4.5 and 22 degrees C and consisted of 425 distance constraints from nuclear Overhauser enhancements and 159 supplementary constraints from spin-spin coupling constants and from the identification of intramolecular hydrogen bonds. Residues 3-30 and 37-48 form a molecular core with two antiparallel beta-sheets and several well-defined turns. The three disulfide bonds 6-14, 16-28, and 22-39 were identified by NMR. In contrast to this well-defined molecular core, with an average root mean square distance for the polypeptide backbone of 0.8 A for a group of nine DISMAN solutions, no preferred conformation was found for the C-terminal segment 49-65, and a loop consisting of residues 31-36 is not uniquely constrained by the NMR data either. These structural properties of recombinant desulfatohirudin coincide closely with the previously described solution conformation of natural hirudin, but the presence of localized differences is indicated by chemical shift differences for residues Asp 5, Ser 9, Leu 15, Asp 53, Gly 54, and Asp 55.     

Solution conformations of human growth hormone releasing factor: comparison of the restrained molecular dynamics and distance geometry methods for a system without long-range distance data

A series of three-dimensional structures of the 1-29 fragment of human growth hormone releasing factor in trifluoroethanol have been determined by molecular dynamics and distance geometry methods. The resulting structures satisfy information from nuclear Overhauser effect (NOE) distance data and an empirical potential energy function. Although the polypeptide was found to have an ordered structure in all simulations, the NOE data were not sufficient for global convergence to a unique three-dimensional geometry. Several satisfactory structures have been determined, all of which are extended conformations consisting of a short beta-strand and two alpha-helices (residues 6-13 and residues 16-29) connected by short segments of less well defined secondary structure. Because of the lack of NOE data connecting the helix segments, their relative orientation is not uniquely determined.     

Conformation of [Cd7]-metallothionein-2 from rat liver in aqueous solution determined by nuclear magnetic resonance spectroscopy

The three-dimensional structure of [Cd7]-metallothionein-2 from rat liver was determined in aqueous solution, using nuclear magnetic resonance spectrometry and distance geometry calculations. The experimental data provided proton-proton distance constraints from measurements of nuclear Overhauser effects, constraints on the geometry of the metal-cysteine clusters determined by heteronuclear correlation spectroscopy, and dihedral angle constraints derived from both coupling constants and nuclear Overhauser effects. The structure calculations were performed with the program DISMAN. As in previous studies with rabbit liver metallothionein-2a, the structure calculations were performed separately for the alpha and beta-domains containing the 4 and 3-metal clusters, respectively, since no interdomain constraints were found. For both domains, the global polypeptide fold, the location of polypeptide secondary structure elements, the architecture of the metal-sulfur cluster and the local chirality of the metal co-ordination are very similar to the solution structure of rabbit metallothionein-2a, but show considerable difference relative to the crystal structure of rat metallothionein-2.     

The structure of the human retinoic acid receptor-beta DNA-binding domain determined by NMR.

The solution structure of the DNA-binding domain (DBD) of the human retinoic acid receptor-beta (hRAR-beta) has been determined by nuclear magnetic resonance (NMR) spectroscopy and distance geometry (DG). The assignments of 1H and 15N resonances were carried out by the use of 1H homonuclear and 15N-1H heteronuclear two- and three-dimensional NMR experiments. The structure of RAR DBD has been obtained on the basis of distance constrains derived from NMR experiments. The structure shows that two "zinc-finger" domains of the protein are followed by two perpendicular alpha-helices and a short beta-sheet near the N-terminus. Apolar residues in both helices form a hydrophobic core. Binding models of RAR DBD to its inverted and direct repeat response elements have been constructed based on this three-dimensional structure.     

Refinement of the three-dimensional solution structure of barley serine proteinase inhibitor 2 and comparison with the structures in crystals.

The three-dimensional structure of barley serine proteinase inhibitor, CI-2, has been determined using nuclear magnetic resonance spectroscopy. The present structure determination is a refinement of the structure previously determined by us, using in the present case stereo-specific assignments, and a virtually complete set of assignments of the two-dimensional nuclear Overhauser spectrum. The structure determination is based on the identification of more than 1300 nuclear Overhauser effects, of which 961 were used in the structure calculation as distance restraints, and on 94 dihedral angle restraints, of which 31 are for chi 1 angles in defined chiral centers. These have been used to calculate a series of 20 three-dimensional structures using a combination of distance geometry, simulated annealing and restrained molecular dynamics. Each of the 20 structures was in agreement within less than 0.5 A of each of the distance restraints and with all dihedral angle restraints. When compared to the geometric average structure of the 20 refined structures the root-mean-square differences for the backbone atoms were 0.8 (+/- 0.2) A and for all atoms were 1.6 (+/- 0.2) A. By comparison, the values obtained for the structures determined previously were 1.4 (+/- 0.2) A and 2.1 (+/- 0.1) A, respectively. The structures were also compared to the structure determined in the crystalline state by X-ray diffraction showing root-mean-square differences of 1.6 (+/- 0.2) A and 2.8 (+/- 0.2) A for the backbone and all atoms, respectively. Common features of the solution structure and the two crystal structures are the four-stranded beta-structure, composed of a pair of parallel strands, and three pairs of antiparallel beta-strands flanked on one side by a 12-residue alpha-helix and on the other side by a loop containing the serine proteinase binding site. The new analysis of the structure has revealed an additional pair of antiparallel beta-strands, consisting of residues 65 to 67 and 81 to 83, that was not seen in either of the crystal structures or the previous solution structure. Identification of this was based on nuclear magnetic resonance evidence for the hydrogen bond (67HN to 81CO) not reported previously. Also the presence of a bifurcated hydrogen bond involving Phe69 CO and HN atoms of Ala77 and Gln78 was observed in solution but not in crystals. Minor differences between the two structures were observed in the phi-angles of residues Met59 and Glu60 in the inhibitory site.     

Structural basis for thermostability of beta-glycosidase from the thermophilic eubacterium Thermus nonproteolyticus HG102

The three-dimensional structure of a thermostable beta-glycosidase (Gly(Tn)) from the thermophilic eubacterium Thermus nonproteolyticus HG102 was determined at a resolution of 2.4 A. The core of the structure adopts the (betaalpha)(8) barrel fold. The sequence alignments and the positions of the two Glu residues in the active center indicate that Gly(Tn) belongs to the glycosyl hydrolases of retaining family 1. We have analyzed the structural features of Gly(Tn) related to the thermostability and compared its structure with those of other mesophilic glycosidases from plants, eubacteria, and hyperthermophilic enzymes from archaea. Several possible features contributing to the thermostability of Gly(Tn) were elucidated.     

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

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

Determination of three-dimensional protein structures from nuclear magnetic resonance data using fragments of known structures

A method to build a three-dimensional protein model from nuclear magnetic resonance (NMR) data using fragments from a data base of crystallographically determined protein structures is presented. The interproton distances derived from the nuclear Overhauser effect (NOE) data are compared to the precalculated distances in the known protein structures. An efficient search algorithm is used, which arranges the distances in matrices akin to a C alpha diagonal distance plot, and compares the NOE distance matrices for short sequential zones of the protein to the data base matrices. After cluster analysis of the fragments found in this way, the structure is built by aligning fragments in overlapping zones. The sequentially long-range NOEs cannot be used in the initial fragments search but are vital to discriminate between several possible combinations of different groups of fragments. The method has been tested on one simulated NOE data set derived from a crystal structure and one experimental NMR data set. The method produces models that have good local structure, but may contain larger global errors. These models can be used as the starting point for further refinement, e.g., by restrained molecular dynamics or interactive graphics.     

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.     

Determination of the positions of bound water molecules in the solution structure of reduced human thioredoxin by heteronuclear three-dimensional nuclear magnetic resonance spectroscopy.

The presence of bound water molecules in the solution structure of reduced human thioredoxin has been investigated using three-dimensional 1H rotating frame Overhauser 1H-15N multiple quantum coherence spectroscopy. It is demonstrated that the backbone amide protons of Lys21, Lys39, Lys82, Gly83 and Asn102, as well as the side-chain amide group of Asn102, are in close proximity to bound water molecules. Examination of the high-resolution solution structure of reduced human thioredoxin reveals that these results are best accounted for by four bound water molecules. Subsequent simulated annealing calculations carried out on the basis of interproton distance and hydrogen bonding restraints to the bound water molecules, supplemented by the original set of experimental restraints used in the calculation of the three-dimensional structure of human thioredoxin, permit a more precise localization of the bound water positions. Potential hydrogen bonds to these water molecules are described and a comparison is made to corresponding bound water molecules in the crystal structure of oxidized Escherichia coli thioredoxin.     

Three-dimensional structure of human [113Cd7]metallothionein-2 in solution determined by nuclear magnetic resonance spectroscopy

The three-dimensional structure of human [113Cd7]metallothionein-2 was determined by nuclear magnetic resonance spectroscopy in solution. Sequence-specific 1H resonance assignments were obtained using the sequential assignment method. The input for the structure calculations consisted of the metal-cysteine co-ordinative bonds identified with heteronuclear correlation spectroscopy, 1H-1H distance constraints from nuclear Overhauser enhancement spectroscopy, and spin-spin coupling constants 3JHN alpha and 3J alpha beta. The molecule consists of two domains, the beta-domain including amino acid residues 1 to 30 and three metal ions, and the alpha-domain including residues 31 to 61 and four metal ions. The nuclear magnetic resonance data present no evidence for a preferred relative orientation of the two domains. The polypeptide-to-metal co-ordinative bonds in human metallothionein-2 are identical to those in the previously determined solution structures of rat metallothionein-2 and rabbit metallothionein-2a, and the polypeptide conformations in the three proteins are also closely similar.     

1H-NMR assignment and secondary structure of human insulin-like growth factor-I (IGF-I) in solution.

Human insulin-like growth factor-I (IGF-I) was studied by two-dimensional 1H-NMR spectroscopy. Resonance assignments were obtained for all the backbone protons and almost all of the sidechain protons of the total 70 amino acid residues, using sequence-specific assignment procedures. The secondary structure elements of human IGF-I were identified by investigation of the sequential and medium range NOEs as a preliminary step in determining the three-dimensional structure of this protein by means of distance geometry calculations. The typical NOEs of d alpha beta(i,i + 3) and d alpha N(i,i + 3), as well as the successive strong NOEs of dNN connectivities and slowly exchanging amide protons confirmed the presence of three helical segments corresponding to the sequence regions, Ala8-Cys18, Gly42-Cys48, and Leu54-Cys61, and the existence of a beta-turn in the Gly19-Gly22 region. Our results definitely indicate that the secondary structure of human IGF-I in solution is consistent with that of insulin in the crystalline state.     

Crystal structure of peptidyl-tRNA hydrolase from mycobacterium smegmatis reveals novel features related to enzyme dynamics

Peptidyl-tRNA hydrolase from Mycobacterium smegmatis is a single domain 21 kDa protein involved in the hydrolysis of prematurely produced peptidyl-tRNAs to ensure the viability of cells in bacteria, thus making it a potentially important drug target. In order to aid the development of potent drugs for controlling bacterial infections, the three-dimensional structure of peptidyl-tRNA hydrolase from Mycobacterium smegmatis has been determined. The protein adopts a compact α/β globular fold with a twisted β-sheet surrounded by α-helices. The functionally important C-terminal stretch has been unambiguously modeled for the first time in the unliganded structure of peptidyl-tRNA hydrolase. The segment, Gly138 - Val150 is mobile because it lacks significant interactions with the rest of the protein molecule. This conformational flexibility is reflected through different values of distances between a reference atom Ala147 C^α of the segment Gly138 - Val150 to Gly114 C^α from another segment from opposite side of the substrate binding channel in Mycobacterium smegmatis (7.8 Ǻ), Mycobacterium tuberculosis (9.5 Ǻ) and Escherichia coli (11.8 Ǻ). Similarly, the conformation of loop Gly109 - Gly117 with respect to another loop Asp95 - Asp100 also shows variability of the substrate binding cleft as the distance between Asp98 O^δ2 to Gly113 C^α in Mycobacterium smegmatis is 4.5 Ǻ while the corresponding distances in Mycobacterium tuberculosis and Escherichia coli are 3.1 Ǻ and 6.7 Ǻ respectively. The hydrogen bonded interactions between Asn116, His22 and Asp95 indicate a stereochemically favorable arrangement of these residues for catalytic action.     

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.     

Proline cis-trans isomers in calbindin D9k observed by X-ray crystallography.

In a structure of recombinant bovine calbindin D9k, determined crystallographically to 1.6 A resolution, a proline in mixed, approximately equally populated, cis and trans conformation is observed. Isomers of this kind have not been reported in structure determinations of calbindin D9k to 2.3 A resolution or in any other crystallographically determined protein structure. The cis-trans isomerization occurs at the peptide bond between Gly42 and Pro43, which is in agreement with results from two-dimensional 1H nuclear magnetic resonance spectroscopy experiments on solutions of calbindin D9k. Alternative backbone stretches have been modeled and refined by stereochemical restrained least-squares refinement for the segment Lys41 to Pro43. The final R-value was 0.188. The structural perturbations accompanying the cis-trans isomerization are found to be very localized. The largest positional differences are observed at residue Gly42, in which the alternative positions of the oxygen atom are 3.6 A apart.     

Activated conformations of the ras-gene-encoded p21 protein. 1. An energy-refined structure for the normal p21 protein complexed with GDP.

A complete three-dimensional structure for the ras-gene-encoded p21 protein with Gly 12 and Gln 61, bound to GDP, has been constructed in four stages using the available alpha-carbon coordinates as deposited in the Brookhaven National Laboratories Protein Data Bank. No all-atom structure has been made available despite the fact that the first crystallographic structure for the p21 protein was reported almost four years ago. In the p21 protein, if amino acid substitutions are made at any one of a number of different positions in the amino acid sequence, the protein becomes permanently activated and causes malignant transformation of normal cells or, in some cell lines, differentiation and maturation. For example, all amino acids except Gly and Pro at position 12 result in an oncogenic protein; all amino acids except Gln, Glu and Pro at position 61 likewise cause malignant transformation of cells. We have constructed our all-atom structure of the non-oncogenic protein from the x-ray structure in order to determine how oncogenic amino acid substitutions affect the three-dimensional structure of this protein. In Stage 1 we generated a poly-alanine backbone (except at Gly and Pro residues) through the alpha-carbon structure, requiring the individual Ala, Pro or Gly residues to conform to standard amino acid geometry and to form trans-planar peptide bonds. Since no alpha-carbon coordinates for residues 60-65 have been determined, these residues were modeled by generating them in the extended conformation and then subjecting them to molecular dynamics using the computer application DISCOVER and energy minimization using DISCOVER and the ECEPP (Empirical Conformational Energies for Peptides Program). In Stage 2, the positions of residues that are homologous to corresponding residues of bacterial elongation factor Tu (EF-Tu) to which p21 bears an overall 40% sequence homology, were determined from their corresponding positions in a high-resolution structure of EF-Tu. Non-homologous loops were taken from the structure generated in Stage 1 and were placed between the appropriate homologous segments so as to connect them. In Stage 3, all bad contacts that occurred in this resulting structure were removed, and the coordinates of the alpha-carbon atoms were forced to superimpose as closely as possible on the corresponding atoms of the reference (x-ray) structure. Then the side chain positions of residues of the non-homologous loop regions were modeled using a combination of molecular dynamics and energy minimization using DISCOVER and ECEPP respectively.(ABSTRACT TRUNCATED AT 400 WORDS)     

Three-dimensional structure of soybean trypsin/chymotrypsin Bowman-Birk inhibitor in solution.

The three-dimensional structure of soybean trypsin/chymotrypsin Bowman-Birk inhibitor in solution has been determined by two-dimensional 1H nuclear magnetic resonance spectroscopy and dynamical simulated annealing using the program XPLOR. The structure was defined by 907 NOEs involving intra- and interresidue contacts which served as distance constraints for a protocol of dynamical simulated annealing. In addition, 48 phi angle constraints involving non-proline amino acids, 29 chi angle constraints, six omega angle constraints for the X-Pro peptide bond, and 35 stereoassignments for prochiral centers were incorporated during the course of the calculation. The protein is characterized by two distinct binding domains for serine protease. Each domain is comprised of a beta-hairpin (antiparallel beta-sheet and a cis-proline-containing type VIb reverse turn) with a short segment making a third strand of antiparallel beta-sheet. The structure determination and refinement are described, and the structure is compared to other structures of Bowman-Birk inhibitors as well as other families of serine protease inhibitors.