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.     

High-resolution solution structure of reduced French bean plastocyanin and comparison with the crystal structure of poplar plastocyanin

The three-dimensional solution structure of reduced (CuI) plastocyanin from French bean leaves has been determined by distance geometry and restrained molecular dynamics methods using constraints obtained from 1H n.m.r. (nuclear magnetic resonance) spectroscopy. A total of 1244 experimental constraints were used, including 1120 distance constraints, 103 dihedral angle constraints and 21 hydrogen bond constraints. Stereospecific assignments were made for 26 methylene groups and the methyls of 11 valines. Additional constraints on copper co-ordination were included in the restrained dynamics calculations. The structures are well defined with average atomic root-mean-square deviations from the mean of 0.45 A for all backbone heavy atoms and 1.08 A for side-chain heavy atoms. French bean plastocyanin adopts a beta-sandwich structure in solution that is similar to the X-ray structure of reduced poplar plastocyanin; the average atomic root-mean-square difference between 16 n.m.r. structures and the X-ray structure is 0.76 A for all backbone heavy atoms. The conformations of the side-chains that constitute the hydrophobic core of French bean plastocyanin are very well defined. Of 47 conserved residues that populate a single chi 1 angle in solution, 43 have the same rotamer in the X-ray structure. Many surface side-chains adopt highly preferred conformations in solution, although the 3J alpha beta coupling constants often indicate some degree of conformational averaging. Some surface side-chains are disordered in both the solution and crystal structures of plastocyanin. There is a striking correlation between measures of side-chain disorder in solution and side-chain temperature factors in the X-ray structure. Side-chains that form a distinctive acidic surface region, believed to be important in binding other electron transfer proteins, appear to be disordered. Fifty backbone amide protons form hydrogen bonds to carbonyls in more than 60% of the n.m.r. structures; 45 of these amide protons exchange slowly with solvent deuterons. Ten hydrogen bonds are formed between side-chain and backbone atoms, eight of which are correlated with decreased proton exchange. Of the 60 hydrogen bonds formed in French bean plastocyanin, 56 occur in the X-ray structure of the poplar protein; two of the missing hydrogen bonds are absent as a result of mutations. It appears that molecular dynamics refinement of highly constrained n.m.r. structures allows accurate prediction of the pattern of hydrogen bonding.     

Application of molecular dynamics with interproton distance restraints to three-dimensional protein structure determination. A model study of crambin

The applicability of restrained molecular dynamics for the determination of three-dimensional protein structures on the basis of short interproton distances (less than 4 A) that can be realistically determined from nuclear magnetic resonance measurements in solution is assessed. The model system used is the 1.2 A resolution crystal structure of the 46 residue protein crambin, from which a set of 240 approximate distance restraints, divided into three ranges (2.5 +/- 0.5, 3.0+0.5(-1.0) and 4 +/- 1 A), is derived. This interproton distance set comprises 159 short-range ([i-j] less than or equal to 5) and 56 ([i-j] greater than 5) long-range inter-residue distances and 25 intra-residue distances. Restrained molecular dynamics are carried out using a number of different protocols starting from two initial structures: a completely extended beta-strand; and an extended structure with two alpha-helices in the same positions as in the crystal structure (residues 7 to 19, and 23 to 30) and all other residues in the form of extended beta-strands. The root-mean-square (r.m.s.) atomic differences between these two initial structures and the crystal structure are 43 A and 23 A, respectively. It is shown that, provided protocols are used that permit the secondary structure elements to form at least partially prior to folding into a tertiary structure, convergence to the correct final structure, both globally and locally, is achieved. The r.m.s. atomic differences between the converged restrained dynamics structures and the crystal structure range from 1.5 to 2.2 A for the backbone atoms and from 2.0 to 2.8 A for all atoms. The r.m.s. atomic difference between the X-ray structure and the structure obtained by first averaging the co-ordinates of the converged restrained dynamics structures is even smaller: 1.0 A for the backbone atoms and 1.6 A for all atoms. These results provide a measure with which to judge future experimental results on proteins whose crystal structures are unknown. In addition, from an examination of the dynamics trajectories, it is shown that the convergence pathways followed by the various simulations are different.     

Nuclear magnetic resonance solution and X-ray structures of squash trypsin inhibitor exhibit the same conformation of the proteinase binding loop

A comparison of the solution nuclear magnetic resonance (n.m.r.) structures of squash trypsin inhibitor from seeds of the squash Cucurbita maxima with the X-ray structure of a trypsin complex of the inhibitor shows that the n.m.r. and X-ray structures are similar in terms of the global folding and secondary structure. The average atomic root-mean-square difference between the 36 n.m.r. structures on the one hand and the X-ray structure is 0.96 A for the backbone atoms and 1.95 A for all heavy atoms. The n.m.r. and X-ray structures exhibit extremely similar conformations of the primary proteinase binding loop. Despite the overall similarity, there are small differences between the mean computed structure and the X-ray structure. The n.m.r. structures have slightly different positions of the segments from residues 16 to 18, and 24 and 25. The n.m.r. results show that the inclusion of stereospecific assignments and precise distance constraints results in a significant improvement in the definition of the n.m.r. structure, making possible a detailed analysis of the local conformations in the protein.     

Comparison of the solution and X-ray structures of barley serine proteinase inhibitor 2

A comparison of the solution n.m.r. structures of barley serine protease inhibitor 2 (BSPI-2) with the X-ray structures of both subtilisin complexed and native BSPI-2 is presented. It is shown that the n.m.r. and X-ray structures are very similar in terms of overall shape, size, polypeptide fold and secondary structure. The average atomic rms difference between the 11 restrained dynamics structures on the one hand and the two X-ray structures on the other is 1.9 +/- 0.2 A for the backbone atoms and 3.0 +/- 0.3 A for all atoms. The corresponding values for the restrained energy minimized mean dynamics structure are 1.5 and 2.4 A, respectively.     

Solution of the phase problem in the X-ray diffraction method for proteins with the nuclear magnetic resonance solution structure as initial model. Patterson search and refinement for the alpha-amylase inhibitor tendamistat

Patterson search calculations using the three-dimensional structure of the alpha-amylase inhibitor from Streptomyces tendae obtained from experimental nuclear magnetic resonance (n.m.r.) data were performed to study the possibility of solving the phase problem in the X-ray diffraction method with protein structures determined by n.m.r. Using all heavy atoms (C, N, O, S) of the residues 5 to 73 in the best n.m.r. structure of the alpha-amylase inhibitor (520 out of the 558 heavy atoms in the complete polypeptide chain), the maximum of the rotation function corresponded to the correct solution obtained by the previous independent determination of the crystal structure. However, additional local maxima, which are not significantly lower than the global maximum, also showed up. Performing the Patterson search with a model containing the backbone atoms and the heavy atoms of only the interior side-chains (399 atoms), which are much better defined by the n.m.r. data, the correct maximum was significantly higher than all other maxima. A translation search for the best orientation of the latter model yielded the correct solution. The energy-restrained crystallographic refinement was performed with this model to an R-factor of 26%. This corresponds approximately to the R-factor calculated for the X-ray crystal structure previously determined using the isomorphous replacement technique, if the residues 1 to 4 and 74 and all localized solvent molecules were removed from this structure. During the refinement the root-mean-square deviation between the two structures decreased from 1.03 A to 0.26 A for the polypeptide backbone and from 1.64 A to 0.73 A for all heavy atoms. There are no major local conformational differences between the two structures, with the single exception of the side-chain of Gln52.     

Comparison of the solution nuclear magnetic resonance and crystal structures of interleukin-8. Possible implications for the mechanism of receptor binding

A comparison of the solution structure of the interleukin-8 dimer determined by nuclear magnetic resonance spectroscopy with that of the 2 A resolution X-ray structure, solved by molecular replacement using the solution structure as a starting model, is presented. At the monomer level the atomic root-mean-square difference between the two structures for residues 7 to 72 is approximately 1.1 A for the backbone atoms, approximately 1.6 A for all atoms, and approximately 1 A for all atoms of the internal residues. There are two main regions of difference in the monomer. In the X-ray structure residues 4 to 6 are well ordered and the charged groups of Glu4 of one subunit and Lys23' of the other are in close enough proximity to form an electrostatic interaction. In contrast, these residues are partially disordered in solution and the electrostatic interaction involving Glu4 is replaced by one between Glu29 of one subunit and Lys23' of the other. In the loop comprising residues 31 to 36, His33 accepts a hydrogen bond from the backbone amide group of Gln8 in the solution structure, but donates a hydrogen bond to the backbone carbonyl group of Glu29 in the X-ray structure. There is also a difference in the quaternary structure with regard to the relative orientation of the two subunits produced by a rigid body rotation about the C2 axis that alters the angle between the central beta-strands (formed by residues 23 to 29 of the 2 subunits) at the dimer interface, without breaking the symmetry. In the solution structure this angle has a value of 168 degrees, while in the X-ray structure the central strands are essentially flat, with an angle of 179 degrees. As a result, the separation between the two anti-parallel helices, which lie at an angle of about 60 degrees to the underlying beta-strands, is decreased from 14.8 A in the solution structure to 11.1 A in the X-ray structure. The quaternary structural difference is related to the different conformations of the N terminus and the 31 to 36 loop, both of which display different interactions with respect to the ends of the central beta-strands in the two structures. These findings indicate that interleukin-8 has the potential to undergo conformational transitions that may be of functional significance.     

Solution structure and backbone dynamics of the catalytic domain of matrix metalloproteinase-2 complexed with a hydroxamic acid inhibitor

MMP-2 is a member of the matrix metalloproteinase family that has been implicated in tumor cell metastasis and angiogenesis. Here, we describe the solution structure of a catalytic domain of MMP-2 complexed with a hydroxamic acid inhibitor (SC-74020), determined by three-dimensional heteronuclear NMR spectroscopy. The catalytic domain, designated MMP-2C, has a short peptide linker replacing the internal fibronectin-domain insertion and is enzymatically active. Distance geometry-simulated annealing calculations yielded 14 converged structures with atomic root-mean-square deviations (r.m.s.d.) of 1.02 and 1.62 A from the mean coordinate positions for the backbone and for all heavy atoms, respectively, when 11 residues at the N-terminus are excluded. The structure has the same global fold as observed for other MMP catalytic domains and is similar to previously solved crystal structures of MMP-2. Differences observed between the solution and the crystal structures, near the bottom of the S1' specificity loop, appear to be induced by the large inhibitor present in the solution structure. The MMP-2C solution structure is compared with MMP-8 crystal structure bound to the same inhibitor to highlight the differences especially in the S1' specificity loop. The finding provides a structural explanation for the selectivity between MMP-2 and MMP-8 that is achieved by large inhibitors.     

The determination of the three-dimensional structure of barley serine proteinase inhibitor 2 by nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The solution structure of the 64 residue structured domain (residues 20-83) of barley serine proteinase inhibitor 2 (BSPI-2) is determined on the basis of 403 interproton distance, 34 phi backbone torsion angle and 26 hydrogen bonding restraints derived from n.m.r. measurements. A total of 11 converged structures were computed using a metric matrix distance geometry algorithm and refined by restrained molecular dynamics. The average rms difference between the final 11 structures and the mean structure obtained by averaging their coordinates is 1.4 +/- 0.2 A for the backbone atoms and 2.1 +/- 0.1 A for all atoms. The overall structure, which is almost identical to that found by X-ray crystallography, is disc shaped and consists of a central four component mixed parallel and antiparallel beta-sheet flanked by a 13 residue alpha-helix on one side and the reactive site loop on the other.     

A molecular dynamics simulation of bacteriophage T4 lysozyme.

An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. The simulation was carried out with all hydrogen atoms modeled explicitly, the inclusion of all 152 crystallographic waters and at a temperature of 300 K. Temporal analysis of the trajectory versus energy, hydrogen bond stability, r.m.s. deviation from the starting crystal structure and radius of gyration, demonstrates that the simulation was both stable and representative of the average experimental structure. Average structural properties were calculated from the enzyme trajectory and compared with the crystal structure. The mean value of the C alpha displacements of the average simulated structure from the X-ray structure was 1.1 +/- 0.1 A; differences of the backbone phi and psi angles between the average simulated structure and the crystal structure were also examined. Thermal-B factors were calculated from the simulation for heavy and backbone atoms and both were in good agreement with experimental values. Relationships between protein secondary structure elements and internal motions were studied by examining the positional fluctuations of individual helix, sheet and turn structures. The structural integrity in the secondary structure units was preserved throughout the simulation; however, the A helix did show some unusually high atomic fluctuations. The largest backbone atom r.m.s. fluctuations were found in non-secondary structure regions; similar results were observed for r.m.s. fluctuations of non-secondary structure phi and psi angles. In general, the calculated values of r.m.s. fluctuations were quite small for the secondary structure elements. In contrast, surface loops and turns exhibited much larger values, being able to sample larger regions of conformational space. The C alpha difference distance matrix and super-positioning analyses comparing the X-ray structure with the average dynamics structure suggest that a 'hinge-bending' motion occurs between the N- and C-terminal domains.     

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

A low resolution solution structure of the cytokine interleukin-1 beta, 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 alpha H and C beta H protons, supplemented by 90 distance restraints for 45 hydrogen bonds, and 79 phi torsion angle restraints. With the exception of 27 C alpha H-C alpha 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 3JHN alpha 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 beta comprises 12 beta-strands arranged in three pseudo-symmetrical topological units (each consisting of 5 anti-parallel beta-strands), joined by turns, short loops and long loops. The core of the structure, which is made up of the 12 beta-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) A. 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 beta with that of the X-ray structure indicates that they are similar, with a backbone atomic r.m.s. difference of only 1.5 A between the co-ordinates of the restrained minimized mean of the SA structures and the X-ray structure.     

Three-dimensional solution structure of a curaremimetic toxin from Naja nigricollis venom: a proton NMR and molecular modeling study.

The solution conformation of toxin alpha from Naja nigricollis (61 amino acids and four disulfides), a snake toxin which specifically blocks the activity of the nicotinic acetylcholine receptor (AcChoR), has been determined using nuclear magnetic resonance spectroscopy and molecular modeling. The solution structures were calculated using 409 distance and 73 dihedral angle restraints. The average atomic rms deviation between the eight refined structures and the mean structure is approximately 0.5 A for the backbone atoms. The overall folding of toxin alpha consists of three major loops which are stabilized by three disulfide bridges and one short C terminal loop stabilized by a fourth disulfide bridge. All the disulfides are grouped in the same region of the molecule, forming a highly constrained structure from which the loops protrude. As predicted, this structure appears to be very similar to the 1.4-A resolution crystal structure of another snake neurotoxin, namely, erabutoxin b from Laticauda semifasciata. The atomic rms deviation for the backbone atoms between the solution and crystal structures is approximately 1.7 A. The minor differences which are observed between the two structures are partly related to the deletion of one residue from the chain of toxin alpha. It is notable that, although the two toxins differ from each other by 16 amino acid substitutions, their side chains have an essentially similar spatial organization. However, most of the side chains which constitute the presumed AcChoR binding site for the curaremimetic toxins are poorly resolved in toxin alpha.(ABSTRACT TRUNCATED AT 250 WORDS)     

Factors influencing accuracy of computer-built models: a study based on leucine zipper GCN4 structure.

A three-dimensional model of the leucine zipper GCN4 built from its amino acid sequence had been reported previously by us. When the two alternative x-ray structures of the GCN4 dimer became available, the root mean square (r.m.s.) shifts between our model and the structures were determined as approximately 2.7 A on all atoms. These values are similar to the r.m.s. shift of 2.8 A between the two GCN4 structures in the different crystal forms (C2 and P2(1)2(1)2(1)). CONGEN conformational searches were run to better understand the conditions that may determine the preference of different conformers in different environments and to test the sensitivity of our current modeling techniques. With a judicious choice of CONGEN search parameters, the backbone r.m.s. deviation improved to 0.8 A and 2.5 A on all atoms. The side-chain conformations of Val and Leu at the helical interface were well reproduced (1.2 A r.m.s.), and the large side-chain misplacements occurred with only a small number of charged amino acids and a tyrosine. Inclusion of the crystal environment (C2 symmetry), as a passive background, into the side-chain conformational search further improved the accuracy of the model to an r.m.s. deviation of 2.1 A. Conformational searches carried out in the two different crystal environments and employing the AMBER protein/DNA forcefield, as implemented in CONGEN, gave the r.m.s. values of 2.2 A (for the C2 symmetry) and 2.5 A (for the P2(1)2(1)2(1) symmetry). In the C2 symmetry crystal, as much as 40% of the surface of each dimer was involved in crystal contacts with other dimers, and the charged residues on the surface often interacted with immobilized water molecules. Thus, occasional large r.m.s. deviations between the model and the x-ray side chains were due to specific conditions that did not occur in solution.     

The polypeptide fold of the globular domain of histone H5 in solution. A study using nuclear magnetic resonance, distance geometry and restrained molecular dynamics.

The polypeptide fold of the 79-residue globular domain of chicken histone H5 (GH5) in solution has been determined by the combined use of distance geometry and restrained molecular dynamics calculations. The structure determination is based on 307 approximate interproton distance restraints derived from n.m.r. measurements. The structure is composed of a core made up of residues 3-18, 23-34, 37-60 and 71-79, and two loops comprising residues 19-22 and 61-70. The structure of the core is well defined with an average backbone atomic r.m.s. difference of 2.3 +/- 0.3 A between the final eight converged restrained dynamics structures and the mean structure obtained by averaging their coordinates best fitted to the core residues. The two loops are also well defined locally but their orientation with respect to the core could not be determined as no long range ([i-j[ greater than 5) proton-proton contacts could be observed between the loop and core residues in the two-dimensional nuclear Overhauser enhancement spectra. The structure of the core is dominated by three helices and has a similar fold to the C-terminal DNA binding domain of the cAMP receptor protein.     

High-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a hydroxamic acid inhibitor

The high-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a sulfonamide derivative of a hydroxamic acid compound (WAY-151693) has been determined by multidimensional heteronuclear NMR. A total of 30 structures were calculated for residues 7-164 by means of hybrid distance geometry-simulated annealing using a total of 3280 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures is 0.43(+/-0.05) A for the backbone atoms, 0.80(+/-0.09) A for all atoms, and 0.47(+/-0.04) A for all atoms excluding disordered side-chains. The overall structure of MMP-13 is composed of a beta-sheet consisting of five beta-strands in a mixed parallel and anti-parallel arrangement and three alpha-helices where its overall fold is consistent with previously solved MMP structures. A comparison of the NMR structure of MMP-13 with the published 1.6 A resolution X-ray structure indicates that the major differences between the structures is associated with loop dynamics and crystal-packing interactions. The side-chains of some active-site residues for the NMR and X-ray structures of MMP-13 adopt distinct conformations. This is attributed to the presence of unique inhibitors in the two structures that encounter distinct interactions with MMP-13. The major structural difference observed between the MMP-13 and MMP-1 NMR structures is the relative size and shape of the S1' pocket where this pocket is significantly longer for MMP-13, nearly reaching the surface of the protein. Additionally, MMP-1 and MMP-13 exhibit different dynamic properties for the active-site loop and the structural Zn-binding region. The inhibitor WAY-151693 is well defined in the MMP-13 active-site based on a total of 52 distance restraints. The binding motif of WAY-151693 in the MMP-13 complex is consistent with our previously reported MMP-1:CGS-27023A NMR structure and is similar to the MMP-13: RS-130830 X-ray structure.     

Three-dimensional solution structure of the E3-binding domain of the dihydrolipoamide succinyltransferase core from the 2-oxoglutarate dehydrogenase multienzyme complex of Escherichia coli.

The three-dimensional solution structure of a 51-residue synthetic peptide comprising the dihydrolipoamide dehydrogenase (E3)-binding domain of the dihydrolipoamide succinyltransferase (E2) core of the 2-oxoglutarate dehydrogenase multienzyme complex of Escherichia coli has been determined by nuclear magnetic resonance spectroscopy and hybrid distance geometry-dynamical simulated annealing calculations. The structure is based on 630 approximate interproton distance and 101 torsion angle (phi, psi, chi 1) restraints. A total of 56 simulated annealing structures were calculated, and the atomic rms distribution about the mean coordinate positions for residues 12-48 of the synthetic peptide is 1.24 A for the backbone atoms, 1.68 A for all atoms, and 1.33 A for all atoms excluding the six side chains which are disordered at chi 1 and the seven which are disordered at chi 2; when the irregular partially disordered loop from residues 31 to 39 is excluded, the rms distribution drops to 0.77 A for the backbone atoms, 1.55 A for all atoms, and 0.89 A for ordered side chains. Although proton resonance assignments for the N-terminal 11 residues and the C-terminal 3 residues were obtained, these two segments of the polypeptide are disordered in solution as evidenced by the absence of nonsequential nuclear Overhauser effects. The solution structure of the E3-binding domain consists of two parallel helices (residues 14-23 and 40-48), a short extended strand (24-26), a five-residue helical-like turn, and an irregular (and more disordered) loop (residues 31-39). This report presents the first structure of an E3-binding domain from a 2-oxo acid dehydrogenase complex.(ABSTRACT TRUNCATED AT 250 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.     

Determination of the nuclear magnetic resonance solution structure of the DNA-binding domain (residues 1 to 69) of the 434 repressor and comparison with the X-ray crystal structure.

The DNA-binding domain of the phage 434 repressor consisting of N-terminal residues 1 to 69 (434 repressor(1-69)), was expressed in Escherichia coli with natural isotope abundance, uniform 15N-labeling and biosynthetically directed fractional 13C-labeling in extent of about 10%. With these protein preparations the three-dimensional structure was determined in solution. The techniques used were nuclear magnetic resonance (n.m.r.) spectroscopy for the collection of conformational constraints, calculation of the protein structure from the n.m.r. data with the program DIANA and structure refinements by restrained energy minimization with a modified version of the program AMBER. A group of 20 conformers characterizes a well-defined structure for residues 1 to 63, with an average of 0-6 A for the root-mean-square deviations (RMSD) calculated for the backbone atoms of the individual conformers relative to the mean co-ordinates. The spatial structure of C-terminal residues 64 to 69 is not defined by the n.m.r. data. The molecular architecture of the 434 repressor(1-69) in solution includes five alpha-helices extending from residues 2 to 13, 17 to 24, 28 to 35, 45 to 52 and 56 to 60, which enclose a well-defined hydrophobic core. The n.m.r. structure is closely similar to the reported crystal structure of the 434 repressor(1-69), with an RMSD value of 1.1 A for the backbone atoms of residues 1 to 63. Small differences between the two structures in regions of the first helix and the loop between helices 3 and 4 were analyzed relative to possible correlations with protein-protein contacts in the crystal lattice and the different milieus of pH and ionic strength in the crystals and n.m.r. samples. Further systematic comparisons of local conformational features indicated that there are correlations between amino acid types, local precision of the structure determination by both techniques and local differences between the structures in the crystals and in solution. Overall, hydrophobic residues are most precisely characterized and agree most closely in the two environments.     

Structure and internal mobility of proteins: a molecular dynamics study of hen egg white lysozyme

The structure and internal motions of the protein hen egg white lysozyme are studied by analysis of simulation and experimental data. A molecular dynamics simulation and an energy minimization of the protein in vacuum have been made and the results compared with high-resolution structures and temperature factors of hen egg white lysozyme in two different crystal forms and of the homologous protein human lysozyme. The structures obtained from molecular dynamics and energy minimization have root-mean-square deviations for backbone atoms of 2.3 Å and 1.1–1.3 Å, respectively, relative to the crystal structures; the different crystal structures have root-mean-square deviations of 0.73–0.81 Å for the backbone atoms. In comparing the backbone dihedral angles, the difference between the dynamics and the crystal structure on which it is based is the same as that between any two crystal structures. The internal fluctuations of atomic positions calculated from the molecular dynamics trajectory agree well with the temperature factors from the three structures. Simulation and crystal results both show that there are large motions for residues involved in exposed turns of the backbone chain, relatively smaller motions for residues involved in the middle of helices or β-sheet structures, and relatively small motions of residues near disulfide bridges. Also, both the simulation and crystal data show that side-chain atoms have larger fluctuations than main-chain atoms. Moreover, the regions that have large deviations among the x-ray crystal structures, which indicates flexibility, are found to have large fluctuations in the simulation.     

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.