Markov dynamic models for long-timescale protein motion

Chiang, Tsung-Han; Hsu, David; Latombe, Jean-Claude
June 2010
Bioinformatics;Jun2010, Vol. 26 Issue 12, pi269
Academic Journal
Molecular dynamics (MD) simulation is a well-established method for studying protein motion at the atomic scale. However, it is computationally intensive and generates massive amounts of data. One way of addressing the dual challenges of computation efficiency and data analysis is to construct simplified models of long-timescale protein motion from MD simulation data. In this direction, we propose to use Markov models with hidden states, in which the Markovian states represent potentially overlapping probabilistic distributions over protein conformations. We also propose a principled criterion for evaluating the quality of a model by its ability to predict long-timescale protein motions. Our method was tested on 2D synthetic energy landscapes and two extensively studied peptides, alanine dipeptide and the villin headpiece subdomain (HP-35 NleNle). One interesting finding is that although a widely accepted model of alanine dipeptide contains six states, a simpler model with only three states is equally good for predicting long-timescale motions. We also used the constructed Markov models to estimate important kinetic and dynamic quantities for protein folding, in particular, mean first-passage time. The results are consistent with available experimental measurements.


Related Articles

  • Identification of the protein folding transition state from molecular dynamics trajectories. Muff, S.; Caflisch, A. // Journal of Chemical Physics;3/28/2009, Vol. 130 Issue 12, p125104 

    The rate of protein folding is governed by the transition state so that a detailed characterization of its structure is essential for understanding the folding process. In vitro experiments have provided a coarse-grained description of the folding transition state ensemble (TSE) of small...

  • Mimicking coarse-grained simulations without coarse-graining: Enhanced sampling by damping short-range interactions. Dongshan Wei; Feng Wang // Journal of Chemical Physics;8/28/2010, Vol. 133 Issue 8, p084101 

    The damped-short-range-interaction (DSRI) method is proposed to mimic coarse-grained simulations by propagating an atomistic scale system on a smoothed potential energy surface. The DSRI method has the benefit of enhanced sampling provided by a typical coarse-grained simulation without the need...

  • Molecular dynamics of folding of secondary structures in Go-type models of proteins. Hoang, Trinh Xuan; Trinh Xuan Hoang; Cieplak, Marek // Journal of Chemical Physics;4/15/2000, Vol. 112 Issue 15 

    We consider six different secondary structures of proteins and construct two types of Go-type off-lattice models: with the steric constraints and without. The basic amino acid-amino acid potential is Lennard-Jones for the native contacts and a soft repulsion for the non-native contacts. The...

  • Absolute comparison of simulated and experimental protein-folding dynamics. Snow, Christopher D.; Nguyen, Houbi; Pande, Vijay S.; Gruebele, Martin // Nature;11/7/2002, Vol. 420 Issue 6911, p102 

    Protein folding is difficult to simulate with classical molecular dynamics. Secondary structure motifs such as α-helices and βhairpins can form in 0.1-10 µs (ref. 1), whereas small proteins have been shown to fold completely in tens of microseconds. The longest folding simulation to...

  • Note: Network random walk model of two-state protein folding: Test of the theory. Berezhkovskii, Alexander M.; Murphy, Ronan D.; Buchete, Nicolae-Viorel // Journal of Chemical Physics;1/21/2013, Vol. 138 Issue 3, p036101 

    We study two-state protein folding in the framework of a toy model of protein dynamics. This model has an important advantage: it allows for an analytical solution for the sum of folding and unfolding rate constants [A. M. Berezhkovskii, F. Tofoleanu, and N.-V. Buchete, J. Chem. Theory Comput....

  • Conformational Transition Pathway and Free Energy Analyses of Proteins by Parallel Cascade Selection Molecular Dynamics (PaCS-MD). Ryuhei Harada; Yasutaka Nishihara; Nobuhiko Wakai; Akio Kitao // AIP Conference Proceedings;2014, Vol. 1618, p86 

    Parallel Cascade Selection Molecular Dynamics (PaCS-MD) was recently proposed to generate conformational transition pathways of proteins under the condition that a certain target quantity to be reached is introduced (R. Harada and A. Kitao, J. Chem. Phys., 139, 035103 2013). In PaCS-MD, the...

  • An intermediate seeks instant gratification. Gruebele, Martin // Nature Structural Biology;Mar2002, Vol. 9 Issue 3, p154 

    Focuses on a φ-value analysis of the four-helix protein Im7 which reveals that a folding intermediate with non-native structure rapidly forms en route to the native state. Factors influencing the folding dynamics; General kinetic equations connecting three states; Range of φ values...

  • Physics and biophysics of solvent induced forces: hydrophobic interactions and context-dependent hydration. San Biagio, P. L.; Bulone, D.; Martorana, V.; Palma-Vittorelli, M. B.; Palma, M. U. // European Biophysics Journal;1998, Vol. 27 Issue 3, p183 

    Solvent induced forces (SIFs) among solutes derive from solvent structural modification due to solutes, and consequent thermodynamic drive towards minimization of related free energy costs. The role of SIFs in biomolecular conformation and function is appreciated by observing that typical SIF...

  • Theory for the rate of contact formation in a polymer chain with local conformational transitions. Zhou, Huan-Xiang // Journal of Chemical Physics;1/22/2003, Vol. 118 Issue 4, p2010 

    I derive an expression for the rate of contact formation between two residues in a polymer chain when both residues undergo native to non-native conformational transitions. A contact can only form when the two residues are brought into contact by interresidue diffusion and are in the native...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics