Simulated and NMR-Derived Backbone Dynamics of a Protein with Significant Flexibility: A Comparison of Spectral Densities for the βARK1 PH Domain
- 9 March 2001
- journal article
- Published by American Chemical Society (ACS) in Journal of the American Chemical Society
- Vol. 123 (13) , 3021-3036
- https://doi.org/10.1021/ja0031117
Abstract
A 7.6 ns molecular dynamics trajectory of the betaARK1 PH domain in explicit water with appropriate ions was calculated at 300 K. Spectral densities at omega = 0, omega(N), and 0.87omega(H) and the model-free parameters were evaluated from the experimental as well as the simulated data, taking the anisotropic overall motion of the protein into account. Experimental and simulated spectral densities are in reasonable general agreement for NH bond vectors, where the corresponding motions have converged within the simulation time. A sufficient sampling of the motions for NH bonds within flexible parts of the protein requires a longer simulation time. The simulated spectral densities J(0) and J(omega(N)) are, on average, 4.5% and 16% lower than the experimental data; the corresponding numbers for the core residues are about 6%; the high-frequency spectral densities J(0.87omega(H)) are lower by, on average, 16% (21% for the core). The simulated order parameters, S(2), are also lower, although the overall disagreement between the simulation and experiment is less pronounced: 1% for all residues and 6% for the core. The observed systematic decrease of simulated spectral density and the order parameters compared to the experimental data can be partially attributed to the ultrafast librational motion of the NH bonds with respect to their peptide plane, which was analyzed in detail. This systematic difference is most pronounced for J(0.87omega(H)), which appears to be most sensitive to the slow, subnanosecond time scale of internal motion, whereas J(0) and J(omega(N)) are dominated by the overall rotational tumbling of the protein. Similar discrepancies are observed between the experimentally measured (15)N relaxation parameters (R(1), R(2), NOE) and their values calculated from the simulated spectral densities. The analysis of spectral densities provides additional information regarding the comparison of the simulated and experimental data, not available from the model-free analysis.Keywords
This publication has 45 references indexed in Scilit:
- Non-native interactions in protein folding intermediates: molecular dynamics simulations of hen lysozymeJournal of Molecular Biology, 1998
- The Solution Structure and Dynamics of the Pleckstrin Homology Domain of G Protein-coupled Receptor Kinase 2 (β-Adrenergic Receptor Kinase 1)Journal of Biological Chemistry, 1998
- The main-chain dynamics of the dynamin pleckstrin homology (PH) domain in solution: analysis of 15N relaxation with monomer/dimer equilibrationJournal of Molecular Biology, 1997
- AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of moleculesComputer Physics Communications, 1995
- Long-Range Motional Restrictions in a Multidomain Zinc-Finger Protein from Anisotropic TumblingScience, 1995
- Backbone Dynamics of Ribonuclease HI: Correlations with Structure and Function in an Active EnzymeJournal of Molecular Biology, 1995
- Characterization of biomolecular structure and dynamics by NMR cross relaxationProgress in Nuclear Magnetic Resonance Spectroscopy, 1994
- A 500 ps molecular dynamics simulation study of interleukin-1β in waterJournal of Molecular Biology, 1992
- Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics, 1983
- Discussion des modes normaux des groupements amides cis et trans à partir des champs de force du σ-valérolactame et du N methylacetamideSpectrochimica Acta Part A: Molecular Spectroscopy, 1973