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.     

Determination of a high-quality nuclear magnetic resonance solution structure of the bovine pancreatic trypsin inhibitor and comparison with three crystal structures.

A high-quality three-dimensional structure of the bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution was determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and compared to the three available high-resolution X-ray crystal structures. A newly collected input of 642 distance constraints derived from nuclear Overhauser effects and 115 dihedral angle constraints was used for the structure calculations with the program DIANA, followed by restrained energy minimization with the program AMBER. The BPTI solution structure is represented by a group of 20 conformers with an average root-mean-square deviation (RMSD) relative to the mean solution structure of 0.43 A for backbone atoms and 0.92 A for all heavy atoms of residues 2 to 56. The pairwise RMSD values of the three crystal structures relative to the mean solution structure are 0.76 to 0.85 A for the backbone atoms and 1.24 to 1.33 A for all heavy atoms of residues 2 to 56. Small local differences in backbone atom positions between the solution structure and the X-ray structures near residues 9, 25 to 27, 46 to 48 and 52 to 58, and conformational differences for individual amino acid side-chains were analyzed for possible correlations with intermolecular protein-protein contacts in the crystal lattices, using the pairwise RMSD values among the three crystal structures as a reference.     

Structural analysis of zinc-substituted cytochrome c

Zinc-substituted cytochrome c has been widely used in studies of protein-protein interactions and photo-induced electron transfer reactions between proteins. However, the coordination geometry of zinc in zinc-substituted cyt c has not yet been determined; two different opinions about the coordination have been reached. Here the solution structures of zinc-substituted cytochrome c that might be five-coordinated and six-coordinated have been refined separately by using (1)H NMR spectroscopy, and the zinc coordination geometry was determined just by NOE distance constraints. Structural analysis of the energy-minimized average solution structures of both the pentacoordinated and hexacoordinated geometries indicate that that zinc in zinc-substituted cyt c should be bound to both His18 and Met80, which means that the zinc is six-coordinated. RMSD values of the family of 25 six-coordinated structures from the average structure are 0.66+/-0.13 A and 1.09+/-0.16 A for the backbone and all heavy atoms, respectively. A statistical analysis of the structure indicates its satisfactory quality. Comparison of the solution structure of the six-coordinated energy-minimized average structure of zinc-substituted cytochrome c with the solution structure of reduced cytochrome c reveals that for the overall folding the secondary structure elements are very close. The availability of the structure provides for a better understanding of the protein-protein complex and for electron transfer processes between Zn cyt c and other metalloproteins.     

NMR structure of lung surfactant peptide SP-B(11-25)

Surfactant protein B (SP-B) is a 79-residue essential component of lung surfactant, the film of lipid and protein lining the alveoli, and is the subject of great interest for its role in lung surfactant replacement therapies. Here we report circular dichroism results and the solution NMR structure of SP-B(11-25) (CRALIKRIQAMIPKG) dissolved in CD(3)OH at 5 degrees C. This is the first report of NMR data related to the protein SP-B, whose structure promises to help elucidate the mechanism of its function. Sequence-specific resonance assignments were made for all observable (1)H NMR signals on the basis of standard 2D NMR methods. Structures were determined by the simulated annealing method using restraints derived from 2D NOESY data. The calculations yielded 17 energy-minimized structures, three of which were subjected to 0.95 ns of restrained dynamics to assess the relevance of the static structures to more realistic dynamic behavior. Our CD and NMR data confirm that this segment is an amphiphilic alpha helix from approximately residue L14 through M21. The backbone heavy-atom RMSD for residues L14 through M21 is 0.09 +/- 0.12 A, and the backbone heavy-atom RMSD for the whole peptide is 0.96 +/- 2.45 A, the difference reflecting fraying at the termini. Aside from the disordered termini, the minimized structures represent dynamic structures well. Structural similarity to the homologous regions of related saposin-like proteins and the importance of the distribution of polar residues about the helix axis are discussed.     

Flexibility of BIV TAR-Tat: models of peptide binding

A new approach in determining local residue flexibility from base-amino acid contact frequencies is applied to the twelve million lattice chains modeling BIV Tat peptide binding to TAR RNA fragment. Many of the resulting key features in flexibility correspond to RMSD calculations derived from a set of five NMR derived structures (X. Ye, R. A. Kumar, and D. J. Patel, Protein Data Bank: Database of three-dimensional structures determined from NMR (1996)) and binding studies of mutants (L. Chen and A. D. Frankel, Proc. Natl. Acad. Sci. USA 92, 5077-5081 (1995)). The lattice and RMSD calculations facilitate the identification of peptide hinge regions that can best utilize the introduction of Gly or other flexible residues. This approach for identifying potential sites amenable to substitution of more flexible residues to enhance peptide binding to RNA targets could be a useful design tool.     

The solution structure of eglin c based on measurements of many NOEs and coupling constants and its comparison with X-ray structures.

A high-precision solution structure of the elastase inhibitor eglin c was determined by NMR and distance geometry calculations. A large set of 947 nuclear Overhauser (NOE) distance constraints was identified, 417 of which were quantified from two-dimensional NOE spectra at short mixing times. In addition, a large number of homonuclear 1H-1H and heteronuclear 1H-15N vicinal coupling constants were used, and constraints on 42 chi 1 and 38 phi angles were obtained. Structure calculations were carried out using the distance geometry program DG-II. These calculations had a high convergence rate, in that 66 out of 75 calculations converged with maximum residual NOE violations ranging from 0.17 A to 0.47 A. The spread of the structures was characterized with average root mean square deviations () between the structures and a mean structure. To calculate the unbiased toward any single structure, a new procedure was used for structure alignment. A canonical structure was calculated from the mean distances, and all structures were aligned relative to that. Furthermore, an angular order parameter S was defined and used to characterize the spread of structures in torsion angle space. To obtain an accurate estimate of the precision of the structure, the number of calculations was increased until the and the angular order parameters stabilized. This was achieved after approximately 40 calculations. The structure consists of a well-defined core whose backbone deviates from the canonical structure ca. 0.4 A, a disordered N-terminal heptapeptide whose backbone deviates by 0.8-12 A, and a proteinase-binding loop whose backbone deviates up to 3.0 A. Analysis of the angular order parameters and inspection of the structures indicates that a hinge-bending motion of the binding loop may occur in solution. Secondary structures were analyzed by comparison of dihedral angle patterns. The high precision of the structure allows one to identify subtle differences with four crystal structures of eglin c determined in complexes with proteinases.     

A large data set comparison of protein structures determined by crystallography and NMR: statistical test for structural differences and the effect of crystal packing

The existence of a large number of proteins for which both nuclear magnetic resonance (NMR) and X-ray crystallographic coordinates have been deposited into the Protein Data Bank (PDB) makes the statistical comparison of the corresponding crystal and NMR structural models over a large data set possible, and facilitates the study of the effect of the crystal environment and other factors on structure. We present an approach for detecting statistically significant structural differences between crystal and NMR structural models which is based on structural superposition and the analysis of the distributions of atomic positions relative to a mean structure. We apply this to a set of 148 protein structure pairs (crystal vs NMR), and analyze the results in terms of methodological and physical sources of structural difference. For every one of the 148 structure pairs, the backbone root-mean-square distance (RMSD) over core atoms of the crystal structure to the mean NMR structure is larger than the average RMSD of the members of the NMR ensemble to the mean, with 76% of the structure pairs having an RMSD of the crystal structure to the mean more than a factor of two larger than the average RMSD of the NMR ensemble. On average, the backbone RMSD over core atoms of crystal structure to the mean NMR is approximately 1 A. If non-core atoms are included, this increases to 1.4 A due to the presence of variability in loops and similar regions of the protein. The observed structural differences are only weakly correlated with the age and quality of the structural model and differences in conditions under which the models were determined. We examine steric clashes when a putative crystalline lattice is constructed using a representative NMR structure, and find that repulsive crystal packing plays a minor role in the observed differences between crystal and NMR structures. The observed structural differences likely have a combination of physical and methodological causes. Stabilizing attractive interactions arising from intermolecular crystal contacts which shift the equilibrium of the crystal structure relative to the NMR structure is a likely physical source which can account for some of the observed differences. Methodological sources of apparent structural difference include insufficient sampling or other issues which could give rise to errors in the estimates of the precision and/or accuracy.     

Application of sparse NMR restraints to large-scale protein structure prediction

The protein structure prediction algorithm TOUCHSTONEX that uses sparse distance restraints derived from NMR nuclear Overhauser enhancement (NOE) data to predict protein structures at low-to-medium resolution was evaluated as follows: First, a representative benchmark set of the Protein Data Bank library consisting of 1365 proteins up to 200 residues was employed. Using N/8 simulated long-range restraints, where N is the number of residues, 1023 (75%) proteins were folded to a C(alpha) root-mean-square deviation (RMSD) from native <6.5 A in one of the top five models. The average RMSD of the models for all 1365 proteins is 5.0 A. Using N/4 simulated restraints, 1206 (88%) proteins were folded to a RMSD <6.5 A and the average RMSD improved to 4.1 A. Then, 69 proteins with experimental NMR data were used. Using long-range NOE-derived restraints, 47 proteins were folded to a RMSD <6.5 A with N/8 restraints and 61 proteins were folded to a RMSD <6.5 A with N/4 restraints. Thus, TOUCHSTONEX can be a tool for NMR-based rapid structure determination, as well as used in other experimental methods that can provide tertiary restraint information.     

Molecular dynamics re-refinement of two different small RNA loop structures using the original NMR data suggest a common structure

Restrained molecular dynamics simulations are a robust, though perhaps underused, tool for the end-stage refinement of biomolecular structures. We demonstrate their utility-using modern simulation protocols, optimized force fields, and inclusion of explicit solvent and mobile counterions-by re-investigating the solution structures of two RNA hairpins that had previously been refined using conventional techniques. The structures, both domain 5 group II intron ribozymes from yeast ai5γ and Pylaiella littoralis, share a nearly identical primary sequence yet the published 3D structures appear quite different. Relatively long restrained MD simulations using the original NMR restraint data identified the presence of a small set of violated distance restraints in one structure and a possibly incorrect trapped bulge nucleotide conformation in the other structure. The removal of problematic distance restraints and the addition of a heating step yielded representative ensembles with very similar 3D structures and much lower pairwise RMSD values. Analysis of ion density during the restrained simulations helped to explain chemical shift perturbation data published previously. These results suggest that restrained MD simulations, with proper caution, can be used to "update" older structures or aid in the refinement of new structures that lack sufficient experimental data to produce a high quality result. Notable cautions include the need for sufficient sampling, awareness of potential force field bias (such as small angle deviations with the current AMBER force fields), and a proper balance between the various restraint weights.     

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.     

Structural consequences of metallothionein dimerization: solution structure of the isolated Cd4-alpha-domain and comparison with the holoprotein dimer

The NMR determination of the structure of Cd(7)-metallothionein was done previously using a relatively large protein concentration that favors dimer formation. The reactivity of the protein is also affected under this condition. To examine the influence of protein concentration on metallothionein conformation, the isolated Cd(4)-alpha-domain was prepared from rabbit metallothionein-2 (MT 2), and its three-dimensional structure was determined by heteronuclear, (1)H-(111)Cd, and homonuclear, (1)H-(1)H NMR, correlation experiments. The three-dimensional structure was refined using distance and angle constraints derived from these two-dimensional NMR data sets and a distance geometry/simulated annealing protocol. The backbone superposition of the alpha-domain from rabbit holoprotein Cd(7)-MT 2 and the isolated rabbit Cd(4)-alpha was measured at a RMSD of 2.0 A. Nevertheless, the conformations of the two Cd-thiolate clusters were distinctly different at two of the cadmium centers. In addition, solvent access to the sulfhydryl ligands of the isolated Cd(4)-alpha cluster was 130% larger due to this small change in cluster geometry. To probe whether these differences were an artifact of the structure calculation, the Cd(4)-alpha-domain structure in rabbit Cd(7)-MT 2 was redetermined, using the previously defined set of NOEs and the present calculation protocol. All calculations employed the same ionic radius for Cd(2+) and same cadmium-thiolate bond distance. The newly calculated structure matched the original with an RMSD of 1.24 A. It is hypothesized that differences in the two alpha-domain structures result from a perturbation of the holoprotein structure because of head-to-tail dimerization under the conditions of the NMR experiments.     

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

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

Clustering algorithms for identifying core atom sets and for assessing the precision of protein structure ensembles

An important open question in the field of NMR-based biomolecular structure determination is how best to characterize the precision of the resulting ensemble of structures. Typically, the RMSD, as minimized in superimposing the ensemble of structures, is the preferred measure of precision. However, the presence of poorly determined atomic coordinates and multiple "RMSD-stable domains"--locally well-defined regions that are not aligned in global superimpositions--complicate RMSD calculations. In this paper, we present a method, based on a novel, structurally defined order parameter, for identifying a set of core atoms to use in determining superimpositions for RMSD calculations. In addition we present a method for deciding whether to partition that core atom set into "RMSD-stable domains" and, if so, how to determine partitioning of the core atom set. We demonstrate our algorithm and its application in calculating statistically sound RMSD values by applying it to a set of NMR-derived structural ensembles, superimposing each RMSD-stable domain (or the entire core atom set, where appropriate) found in each protein structure under consideration. A parameter calculated by our algorithm using a novel, kurtosis-based criterion, the epsilon-value, is a measure of precision of the superimposition that complements the RMSD. In addition, we compare our algorithm with previously described algorithms for determining core atom sets. The methods presented in this paper for biomolecular structure superimposition are quite general, and have application in many areas of structural bioinformatics and structural biology.     

NMR solution structure of an oxaliplatin 1,2-d(GG) intrastrand cross-link in a DNA dodecamer duplex

We have determined, at high resolution, the NMR solution structure of an oxaliplatin-GG DNA dodecamer in the AGGC sequence context by 2D NMR studies. Homonuclear assignment strategies resulted in unambiguous assignment of 203 out of 249 protons, which corresponds to assignment of approximately 81% of the protons. Assignments of H5' and H5" protons were tentative due to resonance overlap. The structure of the oxaliplatin duplex was calculated using the program CNS with a simulated annealing protocol. A total of 510 experimental restraints were employed in the structure calculation. Of 20 calculated structures, the 15 with the lowest energy were accepted as a family. The RMSD of the 15 lowest energy structures was 0.68 A, indicating good structural convergence. The theoretical NOESY spectrum obtained by back-calculation from the final average structure showed excellent agreement with the experimental data, indicating that the final structure was in good agreement with the experimental NMR data. Significant conformational differences were observed between the oxaliplatin-GG 12-mer DNA we studied and all previous solution structures of cisplatin-GG DNA duplexes. For example, the oxaliplatin-GG adduct shows much less distortion at the AG base-pair step than the cisplatin-GG adducts. In addition, the oxaliplatin-GG structure also has a narrow minor groove and an overall axis bend of about 31 degrees, both of which are very different from the recent NMR structures for the cisplatin-GG adducts. These structural differences may explain some of the biological differences between oxaliplatin- and cisplatin-GG adducts.     

High resolution solution structure of apo calcyclin and structural variations in the S100 family of calcium-binding proteins

The three-dimensional solution structure of apo rabbit lung calcyclin has been refined to high resolution through the use of heteronuclear NMR spectroscopy and 13C,15N- enriched protein. Upon completing the assignment of virtually all of the 15N, 13C and 1H NMR resonances, the solution structure was determined from a combination of 2814 NOE- derived distance constraints, and 272 torsion angle constraints derived from scalar couplings. A large number of critical inter- subunit NOEs (386) were identified from 13C- select,13C-filtered NOESY experiments, providing a highly accurate dimer interface. The combination of distance geometry and restrained molecular dynamics calculations yielded structures with excellent agreement with the experimental data and high precision (rmsd from the mean for the backbone atoms in the eight helices: 0.33 Å). Calcyclin exhibits a symmetric dimeric fold of two identical 90 amino acid subunits, characteristic of the S100 subfamily of EF-hand Ca2+-binding proteins. The structure reveals a readily identified pair of putative sites for binding of Zn2+. In order to accurately determine the structural features that differentiate the various S100 proteins, distance difference matrices and contact maps were calculated for the NMR structural ensembles of apo calcyclin and rat and bovine S100B. These data show that the most significant variations among the structures are in the positioning of helix III and in loops, the regions with least sequence similarity. Inter-helical angles and distance differences for the proteins show that the positioning of helix III of calcyclin is most similar to that of bovine S100B, but that the helix interfaces are more closely packed in calcyclin than in either S100B structure. Surprisingly large differences were found in the positioning of helix III in the two S100B structures, despite there being only four non-identical residues, suggesting that one or both of the S100B structures requires further refinement.     

Three-dimensional solution structure and stability of thioredoxin m from spinach.

Proton NMR spectral resonances of thioredoxin m from spinach have been assigned, and its solution structure has been determined on the basis of 1156 nuclear Overhauser effect- (NOE-) derived distance constraints by using restrained molecular dynamics calculations. The average pairwise root-mean-square deviation (RMSD) for the 25 best NMR structures for the backbone was 1.0 +/- 0.1, when the structurally well-defined residues were considered. The N- and C-terminal segments (1-13 and 118-119) and residues 41-49, comprising the active site, are highly disordered. At the time of concluding this work, a crystal structure of this protein was reported, in which thioredoxin m was found to crystallize as noncovalent dimers. Although the solution and crystal structures are very similar, no evidence was found about the existence of dimers in solution, thus confirming that dimerization is not needed for the regulatory activity of thioredoxin m. The spinach thioredoxin m does not unfold by heat in the range 25-85 degrees C, as revealed by thermal circular dichroic (CD) measurements. However, its unfolding free energy (9.1 +/- 0.8 kcal mol(-1), at pH 5.3 and 25 degrees C) could be determined by extrapolating the free energy values obtained at different concentrations of guanidinium chloride (GdmCl). The folding-unfolding process is two-state as indicated by the coincidence of the CD denaturation curves obtained at far and near UV. The H/D exchange behavior of backbone amide protons was analyzed. The slowest-exchanging protons, requiring a global-unfolding mechanism in order to exchange, are those from beta2, beta3, and beta4, the central strands of the beta-sheet, which constitute the main element of the core of the protein. The free energies obtained from exchange measurements of protons belonging to the alpha-helices are lower than those derived from GdmCl denaturation studies, indicating that those protons exchange by local-unfolding mechanisms.     

Three-dimensional structure of the neurotoxin ATX Ia from Anemonia sulcata in aqueous solution determined by nuclear magnetic resonance spectroscopy

With the aid of 1H nuclear magnetic resonance (NMR) spectroscopy, the three-dimensional structure in aqueous solution was determined for ATX Ia, which is a 46 residue polypeptide neurotoxin of the sea anemone Anemonia sulcata. The input for the structure calculations consisted of 263 distance constraints from nuclear Overhauser effects (NOE) and 76 vicinal coupling constants. For the structure calculation several new or ammended programs were used in a revised strategy consisting of five successive computational steps. First, the program HABAS was used for a complete search of all backbone and chi 1 conformations that are compatible with the intraresidual and sequential NMR constraints. Second, using the program DISMAN, we extended this approach to pentapeptides by extensive sampling of all conformations that are consistent with the local and medium-range NMR constraints. Both steps resulted in the definition of additional dihedral angle constraints and in stereospecific assignments for a number of beta-methylene groups. In the next two steps DISMAN was used to obtain a group of eight conformers that contain no significant residual violations of the NMR constraints or van der Waals contacts. Finally, these structures were subjected to restrained energy refinement with a modified version of the molecular mechanics module of AMBER, which in addition to the energy force field includes potentials for the NOE distance constraints and the dihedral angle constraints. The average of the pairwise minimal RMS distances between the resulting refined conformers calculated for the well defined molecular core, which contains the backbone atoms of 35 residues and 20 interior side chains, is 1.5 +/- 0.3 A. This core is formed by a four-stranded beta-sheet connected by two well-defined loops, and there is an additional flexible loop consisting of the eleven residues 8-18. The core of the protein is stabilized by three disulfide bridges, which are surrounded by hydrophobic residues and shielded on one side by hydrophilic residues.