Dataset Information

Fluctuations and correlations in crystalline protein dynamics: a simulation analysis of staphylococcal nuclease.

ABSTRACT: Understanding collective motions in protein crystals is likely to furnish insight into functional protein dynamics and will improve models for refinement against diffraction data. Here, four 10 ns molecular dynamics simulations of crystalline Staphylococcal nuclease are reported and analyzed in terms of fluctuations and correlations in atomic motion. The simulation-derived fluctuations strongly correlate with, but are slightly higher than, the values derived from the experimental B-factors. Approximately 70% of the atomic fluctuations are due to internal protein motion. For 65% of the protein atoms the internal fluctuations converge on the nanosecond timescale. Convergence is much slower for the elements of the interatomic displacement correlation matrix--of these, >80% converge within 1 ns for interatomic distances less, approximately <6 A, but only 10% for separations approximately =12 A. Those collective motions that converged on the nanosecond timescale involve mostly correlations within the beta-barrel or between alpha-helices of the protein. The R-factor with the experimental x-ray diffuse scattering for the crystal, which is determined by the displacement variance-covariance matrix, decreases to 8% after 10 ns simulation. Both the number of converged correlation matrix elements and the R-factor depend logarithmically on time, consistent with a model in which the number of energy minima sampled depends exponentially on the maximum energy barrier crossed. The logarithmic dependence is also extrapolated to predict a convergence time for the whole variance-covariance matrix of approximately 1 micros.

SUBMITTER: Meinhold L

PROVIDER: S-EPMC1305352 | biostudies-literature | 2005 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Fluctuations and correlations in crystalline protein dynamics: a simulation analysis of staphylococcal nuclease.

Meinhold Lars L Smith Jeremy C JC

Biophysical journal 20050128 4

Understanding collective motions in protein crystals is likely to furnish insight into functional protein dynamics and will improve models for refinement against diffraction data. Here, four 10 ns molecular dynamics simulations of crystalline Staphylococcal nuclease are reported and analyzed in terms of fluctuations and correlations in atomic motion. The simulation-derived fluctuations strongly correlate with, but are slightly higher than, the values derived from the experimental B-factors. Appr ...[more]

PMID: 15681654

Similar Datasets

Project description:To search for submolecular foldon units, the spontaneous reversible unfolding and refolding of staphylococcal nuclease under native conditions was studied by a kinetic native-state hydrogen exchange (HX) method. As for other proteins, it appears that staphylococcal nuclease is designed as an assembly of well-integrated foldon units that may define steps in its folding pathway and may regulate some other functional properties. The HX results identify 34 amide hydrogens that exchange with solvent hydrogens under native conditions by way of large transient unfolding reactions. The HX data for each hydrogen measure the equilibrium stability (Delta G(HX)) and the kinetic unfolding and refolding rates (k(op) and k(cl)) of the unfolding reaction that exposes it to exchange. These parameters separate the 34 identified residues into three distinct HX groupings. Two correspond to clearly defined structural units in the native protein, termed the blue and red foldons. The remaining HX grouping contains residues, not well separated by their HX parameters alone, that represent two other distinct structural units in the native protein, termed the green and yellow foldons. Among these four sets, a last unfolding foldon (blue) unfolds with a rate constant of 6 x 10(-6) s(-1) and free energy equal to the protein's global stability (10.0 kcal/mol). It represents part of the beta-barrel, including mutually H-bonding residues in the beta 4 and beta 5 strands, a part of the beta 3 strand that H-bonds to beta 5, and residues at the N-terminus of the alpha2 helix that is capped by beta 5. A second foldon (green), which unfolds and refolds more rapidly and at slightly lower free energy, includes residues that define the rest of the native alpha2 helix and its C-terminal cap. A third foldon (yellow) defines the mutually H-bonded beta1-beta2-beta 3 meander, completing the native beta-barrel, plus an adjacent part of the alpha1 helix. A final foldon (red) includes residues on remaining segments that are distant in sequence but nearly adjacent in the native protein. Although the structure of the partially unfolded forms closely mimics the native organization, four residues indicate the presence of some nonnative misfolding interactions. Because the unfolding parameters of many other residues are not determined, it seems likely that the concerted foldon units are more extensive than is shown by the 34 residues actually observed.

Project description:The nucleases A produced by two strains of Staphylococcus aureus, which have different stabilities, differ only in the identity of the single amino acid at residue 124. The nuclease from the Foggi strain of S. aureus (by convention nuclease WT), which contains His124, is 1.9 kcal.mol-1 less stable (at pH 5.5 and 20 degrees C) than the nuclease from the V8 strain (by convention nuclease H124L), which contains Leu124. In addition, the population of the trans conformer at the Lys116-Pro117 peptide bond, as observed by NMR spectroscopy, is different for the two variants: about 15% for nuclease WT and 9% for nuclease H124L. In order to improve our understanding of the origin of these differences, we compared the properties of WT and H124L with those of the H124A and H124I variants. We discovered a correlation between effects of different residues at this position on protein stability and on stabilization of the cis configuration of the Lys116-Pro117 peptide bond. In terms of free energy, approximately 17% of the increase in protein stability manifests itself as stabilization of the cis configuration at Lys116-Pro117. This result implies that the differences in stability arise mainly from structural differences between the cis configurational isomers at Pro117 of the different variants at residue 124. We solved the X-ray structure of the cis form of the most stable variant, H124L, and compared it with the published high-resolution X-ray structure of the cis form of the most stable variant, WT (Hynes TR, Fox RO, 1991, Proteins Struct Funct Genet 10:92-105). The two structures are identical within experimental error, except for the side chain at residue 124, which is exposed in the models of both variants. Thus, the increased stability and changes in the trans/cis equilibrium of the Lys116-Pro117 peptide bond observed in H124L relative to WT are due to subtle structural changes that are not observed by current structure determination technique. Residue 124 is located in a helix. However, the stability changes are too large and follow the wrong order of stability to be explained simply by differences in helical propensity. A second site of conformational heterogeneity in native nuclease is found at the His46-Pro47 peptide bond, which is approximately 80% trans in both WT and H124L. Because proline to glycine substitutions at either residue 47 or 117 remove the structural heterogeneity at that position and increase protein stability, we determined the X-ray structures of H124L + P117G and H124L + P47G + P117G and the kinetic parameters of H124L, H124L + P47G, H124L + P117G, and H124L + P47G + P117G. The individual P117G and P47G mutations cause decreases in nuclease activity, with kcat affected more than Km, and their effects are additive. The P117G mutation in nuclease H124L leads to the same local conformational rearrangement described for the P117G mutant of WT (Hynes TR, Hodel A, Fox RO, 1994, Biochemistry 33:5021-5030). In both P117G mutants, the loop formed by residues 112-117 is located closer to the adjacent loop formed by residues 77-85, and residues 115-118 adopt a type I' beta-turn conformation with the Lys116-Gly117 peptide bond in the trans configuration, as compared with the parent protein in which these residues have a typeVIa beta-turn conformation with the Lys116-Pro117 peptide bond in the cis configuration. Addition of the P47G mutation appears not to cause any additional structural changes. However, the electron density for part of the loop containing this peptide bond was not strong enough to be interpreted.

Dataset Information

Fluctuations and correlations in crystalline protein dynamics: a simulation analysis of staphylococcal nuclease.

Publications

Fluctuations and correlations in crystalline protein dynamics: a simulation analysis of staphylococcal nuclease.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets