Dataset Information

Picosecond fluctuating protein energy landscape mapped by pressure temperature molecular dynamics simulation.

ABSTRACT: Microscopic statistical pressure fluctuations can, in principle, lead to corresponding fluctuations in the shape of a protein energy landscape. To examine this, nanosecond molecular dynamics simulations of lysozyme are performed covering a range of temperatures and pressures. The well known dynamical transition with temperature is found to be pressure-independent, indicating that the effective energy barriers separating conformational substates are not significantly influenced by pressure. In contrast, vibrations within substates stiffen with pressure, due to increased curvature of the local harmonic potential in which the atoms vibrate. The application of pressure is also shown to selectively increase the damping of the anharmonic, low-frequency collective modes in the protein, leaving the more local modes relatively unaffected. The critical damping frequency, i.e., the frequency at which energy is most efficiently dissipated, increases linearly with pressure. The results suggest that an invariant description of protein energy landscapes should be subsumed by a fluctuating picture and that this may have repercussions in, for example, mechanisms of energy dissipation accompanying functional, structural, and chemical relaxation.

SUBMITTER: Meinhold L

PROVIDER: S-EPMC2077243 | biostudies-literature | 2007 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Picosecond fluctuating protein energy landscape mapped by pressure temperature molecular dynamics simulation.

Meinhold Lars L Smith Jeremy C JC Kitao Akio A Zewail Ahmed H AH

Proceedings of the National Academy of Sciences of the United States of America 20071023 44

Microscopic statistical pressure fluctuations can, in principle, lead to corresponding fluctuations in the shape of a protein energy landscape. To examine this, nanosecond molecular dynamics simulations of lysozyme are performed covering a range of temperatures and pressures. The well known dynamical transition with temperature is found to be pressure-independent, indicating that the effective energy barriers separating conformational substates are not significantly influenced by pressure. In co ...[more]

PMID: 17956984

Similar Datasets

Project description:RNA hairpin loops are one of the most abundant secondary structure elements and participate in RNA folding and protein-RNA recognition. To characterize the free energy surface of RNA hairpin folding at an atomic level, we calculated the potential of mean force (PMF) as a function of the end-to-end distance, by using umbrella sampling simulations in explicit solvent. Two RNA hairpins containing tetraloop cUUCGg and cUUUUg are studied with AMBER ff99 and CHARMM27 force fields. Experimentally, the UUCG hairpin is known to be significantly more stable than UUUU. In this study, the calculations using AMBER force field give a qualitatively correct description for the folding of two RNA hairpins, as the calculated PMF confirms the global stability of the folded structures and the resulting relative folding free energy is in quantitative agreement with the experimental result. The hairpin stabilities are also correctly differentiated by the more rapid molecular mechanics-Poisson Boltzmann-surface area approach, but the relative free energy estimated from this method is overestimated. The free energy profile shows that the native state basin and the unfolded state plateau are separated by a wide shoulder region, which samples a variety of native-like structures with frayed terminal basepair. The calculated PMF lacks major barriers that are expected near the transition regions, and this is attributed to the limitation of the 1-D reaction coordinate. The PMF results are compared with other studies of small RNA hairpins using kinetics method and coarse grained models. The two RNA hairpins described by CHARMM27 are significantly more deformable than those represented by AMBER. Compared with the AMBER results, the CHARMM27 calculated DeltaG(fold) for the UUUU tetraloop is in better agreement with the experimental results. However, the CHARMM27 calculation does not confirm the global stability of the experimental UUCG structure; instead, the extended conformations are predicted to be thermodynamically stable in solution. This finding is further supported by separate unrestrained CHARMM27 simulations, in which the UUCG hairpin unfolds spontaneously within 10 ns. The instability of the UUCG hairpin originates from the loop region, and propagates to the stem. The results of this study provide a molecular picture of RNA hairpin unfolding and reveal problems in the force field descriptions for the conformational energy of certain RNA hairpin.

Dataset Information

Picosecond fluctuating protein energy landscape mapped by pressure temperature molecular dynamics simulation.

Publications

Picosecond fluctuating protein energy landscape mapped by pressure temperature molecular dynamics simulation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets