Hybrid Distance Geometry-dynamical Simulated Annealing Research Articles

Solution conformation of cobrotoxin: A nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing study

The solution conformation of cobrotoxin has been determined by using proton nuclear magnetic resonance spectroscopy. With the combination of various two-dimensional NMR techniques, the 1H-NMR spectrum of cobrotoxin was completely assigned (Yu et al., 1990). A set of 435 approximate interproton distance restraints was derived from nuclear Overhauser enhancement (NOE) measurements. These NOE constraints, in addition to the 29 dihedral angle constraints (from coupling constant measurements) and 26 hydrogen bonding restraints (from the pattern of short-range NOEs), form the basis of 3-D structure determination by the hybrid distance geometry-dynamical simulated annealing method. The 23 structures that were obtained satisfy the experimental restraints, display small deviation from idealized covalent geometry, and possess good nonbonded contacts. Analysis of converged structures indicated that there are two antiparallel beta sheets (double and triple stranded), duly confirming our earlier observations. These are well defined in terms of both atomic root mean square (RMS) differences and backbone torsional angles. The average backbone RMS deviation between the calculated structures and the mean structure, for the beta-sheet regions, is 0.92 A. The mean solution structure was compared with the X-ray crystal structure of erabutoxin b, the homologous protein. This yielded information that both structures resemble each other except at the exposed loop/surface regions, where the solution structure seems to possess more flexibility.

Solution structure of a polypeptide dimer comprising the fourth calcium-binding site of troponin C by nuclear magnetic resonance spectroscopy

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

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

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

High-resolution three-dimensional structure of a single zinc finger from a human enhancer binding protein in solution

The three-dimensional structure of a 30-residue synthetic peptide containing the carboxy-terminal "zinc finger" motif of a human enhancer binding protein has been determined by two-dimensional nuclear magnetic resonance (2D NMR) spectroscopy and hybrid distance geometry-dynamical simulated annealing calculations. The structure determination is based on 487 approximate interproton distance and 63 torsion angle (phi, psi, and chi 1) restraints. A total of 40 simulated annealing structures were calculated, and the atomic rms distribution about the mean coordinate positions (excluding residues 29 and 30 which are ill-defined) is 0.4 A for the backbone atoms, 0.8 A for all atoms, and 0.41 A for all atoms excluding the lysine and arginine side chains, which are disordered. The solution structure of the zinc finger consists of two irregular antiparallel beta-strands connected by an atypical turn (residues 3-12) and a classical alpha-helix (residues 14-24). The zinc is tetrahedrally coordinated to the sulfur atoms of two cysteines (Cys-5 and Cys-8) and to the N epsilon 2 atoms of two histidines (His-21 and His-27). The two cysteine residues are located in the turn connecting the two beta-strands (residues 5-8); one of the histidine ligands (His-21) is in the alpha-helix, while the second histidine (His-27) is at the end of a looplike structure (formed by the end of the alpha-helix and a turn). The general architecture is qualitatively similar to two previously determined low-resolution Cys2-His2 zinc finger structures, although distinct differences can be observed in the beta-strands and turn and in the region around the two histidines coordinated to zinc. Comparison of the overall polypeptide fold of the enhancer binding protein zinc finger with known structures in the crystallographic data base reveals a striking similarity to one region (residues 23-44) of the X-ray structure of proteinase inhibitor domain III of Japanese quail ovomucoid [Papamokos, E., Weber, E., Bode, W., Huber, R., Empie, M. W., Kato, I., & Laskowski, M. (1982) J. Mol. Biol. 158, 515-537], which could be superimposed with a backbone atomic rms difference of 0.95 A on residues 3-25 (excluding residue 6) of the zinc finger from the enhancer binding protein. The presence of structural homology between two proteins of very different function may indicate that the so-called zinc finger motif is not unique for a class of DNA binding proteins but may represent a general folding motif found in a variety of proteins irrespective of their function.

Low resolution structure of interleukin-1β in solution derived from 1H 15N heteronuclear three-dimensional nuclear magnetic resonance spectroscopy

A low resolution solution structure of the cytokine interleukin-1β, a 153 residue protein of molecular weight 17,400, has been determined on the basis of 446 nuclear Overhauser effect (NOE) derived approximate interproton distance restraints involving solely NH, C αH and C βH protons, supplemented by 90 distance restraints for 45 hydrogen bonds, and 79 φ torsion angle restraints. With the exception of 27 C αHC αH NOEs, all the NOEs were assigned from a three-dimensional 1H 1H NOE 15N 1H heteronuclear multiple quantum coherence (HMQC) spectrum. The torsion angle restraints were obtained from accurate 3 J HNα coupling constants measured from a HMQC-J spectrum, while the hydrogen bonds were derived from a qualitative analysis of the NOE, coupling constant and amide exchange data. A total of 20 simulated annealing (SA) structures was computed using the hybrid distance geometry-dynamical simulated annealing method. The solution structure of IL-1β comprises 12 β-strands arranged in three pseudo-symmetrical topological units (each consisting of 5 anti-parallel β-strands), joined by turns, short loops and long loops. The core of the structure, which is made up of the 12 β-strands, together with the turns joining strands I and II, strands VIII and IX and strands X and XI, is well determined with a backbone atomic root-mean-square (r.m.s.) distribution about the mean co-ordinate positions of 1.2(± 0.1) Å. The loop conformations, on the other hand, are poorly determined by the current data. A comparison of the core of the low resolution solution structure of IL-1β with that of the X-ray structure indicates that they are similar, with a backbone atomic r.m.s. difference of only 1.5 Å between the co-ordinates of the restrained minimized mean of the SA structures and the X-ray structure.

Three-dimensional structure of interleukin 8 in solution

The solution structure of the interleukin 8 (IL-8) dimer has been solved by nuclear magnetic resonance (NMR) spectroscopy and hybrid distance geometry-dynamical simulated annealing calculations. The structure determination is based on a total of 1880 experimental distance restraints (of which 82 are intersubunit) and 362 torsion angle restraints (comprising phi, psi, and chi 1 torsion angles). A total of 30 simulated annealing structures were calculated, and the atomic rms distribution about the mean coordinate positions (excluding residues 1-5 of each subunit) is 0.41 +/- 0.08 A for the backbone atoms and 0.90 +/- 0.08 A for all atoms. The three-dimensional solution structure of the IL-8 dimer reveals a structural motif in which two symmetry-related antiparallel alpha-helices, approximately 24 A long and separated by about 14 A, lie on top of a six-stranded antiparallel beta-sheet platform derived from two three-stranded Greek keys, one from each monomer unit. The general architecture is similar to that of the alpha 1/alpha 2 domains of the human class I histocompatibility antigen HLA-A2. It is suggested that the two alpha-helices form the binding site for the cellular receptor and that the specificity of IL-8, as well as that of a number of related proteins involved in cell-specific chemotaxis, mediation of cell growth, and the inflammatory response, is achieved by the distinct distribution of charged and polar residues at the surface of the helices.

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 φ 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) Å for the backbone atoms and 0.89(±0.17) Å for all heavy atoms. The precision of the structure determination is discussed and correlated to the experimental input data.

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)

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

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

Solution structure of recombinant hirudin and the Lys-47 .fwdarw. Glu mutants: a nuclear magnetic resonance and hybrid distance geometry-dynamical simulated annealing study

The solution structure of recombinant wild-type hirudin and of the putative active site mutant Lys-47----Glu has been investigated by nuclear magnetic resonance (NMR) spectroscopy at 600 MHz. The 1H NMR spectra of the two hirudin variants are assigned in a sequential manner with a combination of two-dimensional NMR techniques. Some assignments made in our previous paper [Sukumaran, D. K., Clore, G. M., Preuss, A., Zarbock, J., & Gronenborn, A. M. (1987) Biochemistry 26, 333-338] were found to be incorrect and are now corrected. Analysis of the NOE data indicates that hirudin consists of an N-terminal compact domain (residues 1-49) held together by three disulfide linkages and a disordered C-terminal tail (residues 50-65) which does not fold back on the rest of the protein. This last observation corrects conclusions drawn by us previously on hirudin extracted from its natural source, the leech Hirudo medicinalis. The improved sensitivity of the 600-MHz spectrometer relative to that of our old 500-MHz spectrometer, the availability of two variants with slightly different chemical shifts, and the additional information arising from stereospecific assignments of methylene beta-protons and methyl protons of valine have permitted the determination of the solution structure of hirudin with much greater precision than before. Structure calculations on the N-terminal domain using the hybrid distance geometry-dynamical simulated annealing method were based on 685 and 661 approximate interproton distance restraints derived from nuclear Overhauser enhancement (NOE) data for the wild-type and mutant hirudin, respectively, together with 16 distance restraints for 8 backbone hydrogen bonds identified on the basis of NOE and amide NH exchange data and 26 phi backbone and 18 chi 1 side-chain torsion angle restraints derived from NOE and three-bond coupling constant data. A total of 32 structures were computed for both the wild-type and mutant hirudin. The structure of residues 2-30 and 37-48 which form the core of the N-terminal domain is well determined in both cases with an average atomic rms difference between the individual structures and the respective mean structures of approximately 0.7 A for the backbone atoms and approximately 1 A for all atoms. As found previously, the orientation of the exposed finger of antiparallel beta-sheet (residues 31-36) with respect to the core could not be determined on the basis of the present data due to the absence of any long-range NOEs between the exposed finger and the core.(ABSTRACT TRUNCATED AT 250 WORDS)

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

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

Determination of three-dimensional structures of proteins from interproton distance data by hybrid distance geometry-dynamical simulated annealing calculations

A new hybrid distance space-real space method for determining three-dimensional structures of proteins on the basis of interproton distance restraints is presented. It involves the following steps: (i) the approximate polypeptide fold is obtained by generating a set of substructures comprising only a small subset of atoms by projection from multi-dimensional distance space into three-dimensional cartesian coordinate space using a procedure known as ‘embedding’; (ii) all remaining atoms are then added by best fitting extended amino acids one residue at a time to the substructures; (iii) the resulting structures are used as the starting point for real space dynamical simulated annealing calculations. The latter involve heating the system to a high temperature followed by slow cooling in order to overcome potential barriers along the pathway towards the global minimum region. This is carried out by solving Newton's equations of motion. Unlike conventional restrained molecular dynamics, however, the non-bonded interactions are represented by a simple van der Waals repulsion term. The method is illustrated by calculations on crambin (46 residues) and the globular domain of histone H5 (79 residues). It is shown that the hybrid method is more efficient computationally and samples a larger region of conformational space consistent with the experimental data than full metric matrix distance geometry calculations alone, particularly for large systems.