Parallel Cascade Selection Molecular Dynamics (PaCS-MD) to generate conformational transition pathway.

Parallel Cascade Selection Molecular Dynamics (PaCS-MD) is proposed as a molecular simulation method to generate conformational transition pathway under the condition that a set of "reactant" and "product" structures is known a priori. In PaCS-MD, the cycle of short multiple independent molecular dynamics simulations and selection of the structures close to the product structure for the next cycle are repeated until the simulated structures move sufficiently close to the product. Folding of 10-residue mini-protein chignolin from the extended to native structures and open-close conformational transition of T4 lysozyme were investigated by PaCS-MD. In both cases, tens of cycles of 100-ps MD were sufficient to reach the product structures, indicating the efficient generation of conformational transition pathway in PaCS-MD with a series of conventional MD without additional external biases. Using the snapshots along the pathway as the initial coordinates, free energy landscapes were calculated by the combination with multiple independent umbrella samplings to statistically elucidate the conformational transition pathways.

[1]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

[2]  Jianpeng Ma,et al.  The unfolding action of GroEL on a protein substrate. , 2004, Biophysical journal.

[3]  Amedeo Caflisch,et al.  Calculation of conformational transitions and barriers in solvated systems: Application to the alanine dipeptide in water , 1999 .

[4]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[5]  Shinya Honda,et al.  10 residue folded peptide designed by segment statistics. , 2004, Structure.

[6]  E. Evans,et al.  Sensitive force technique to probe molecular adhesion and structural linkages at biological interfaces. , 1995, Biophysical journal.

[7]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[8]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[9]  Steven M. Block,et al.  Force and velocity measured for single kinesin molecules , 1994, Cell.

[10]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[11]  I. Tinoco,et al.  Equilibrium Information from Nonequilibrium Measurements in an Experimental Test of Jarzynski's Equality , 2002, Science.

[12]  B. Matthews,et al.  Conservation of solvent‐binding sites in 10 crystal forms of T4 lysozyme , 1994, Protein science : a publication of the Protein Society.

[13]  Jürgen Schlitter,et al.  Targeted Molecular Dynamics Simulation of Conformational Change-Application to the T ↔ R Transition in Insulin , 1993 .

[14]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[15]  Kentaro Shimizu,et al.  Folding free‐energy landscape of a 10‐residue mini‐protein, chignolin , 2006, FEBS letters.

[16]  B. Matthews,et al.  A mutant T4 lysozyme displays five different crystal conformations , 1990, Nature.

[17]  Fumio Hirata,et al.  The effects of solvent on the conformation and the collective motions of protein: normal mode analysis and molecular dynamics simulations of melittin in water and in vacuum , 1991 .

[18]  A. Kitao,et al.  The Fast-Folding Mechanism of Villin Headpiece Subdomain Studied by Multiscale Distributed Computing. , 2012, Journal of chemical theory and computation.

[19]  Makoto Taiji,et al.  Folding Dynamics of 10‐Residue β‐Hairpin Peptide Chignolin , 2007 .

[20]  K. Schulten,et al.  Steered molecular dynamics investigations of protein function. , 2001, Journal of molecular graphics & modelling.

[21]  K. Schulten,et al.  Extraction of Lipids from Phospholipid Membranes by Steered Molecular Dynamics , 1997 .

[22]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[23]  M. Karplus,et al.  Locally accessible conformations of proteins: Multiple molecular dynamics simulations of crambin , 1998, Protein science : a publication of the Protein Society.

[24]  P. Tavan,et al.  Ligand Binding: Molecular Mechanics Calculation of the Streptavidin-Biotin Rupture Force , 1996, Science.

[25]  O. Berger,et al.  Adhesion forces of lipids in a phospholipid membrane studied by molecular dynamics simulations. , 1998, Biophysical journal.

[26]  E. Vanden-Eijnden,et al.  On-the-fly string method for minimum free energy paths calculation , 2007 .

[27]  S. Nosé A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[28]  Christoph Dellago,et al.  Efficient transition path sampling: Application to Lennard-Jones cluster rearrangements , 1998 .

[29]  N. Go,et al.  Effect of solvent on collective motions in globular protein. , 1993, Journal of molecular biology.

[30]  Y. Joti,et al.  Detecting coupled collective motions in protein by independent subspace analysis. , 2010, The Journal of chemical physics.

[31]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[32]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[33]  V. Pande,et al.  The Trp cage: folding kinetics and unfolded state topology via molecular dynamics simulations. , 2002, Journal of the American Chemical Society.

[34]  Caitlin M. Davis,et al.  Raising the speed limit for β-hairpin formation. , 2012, Journal of the American Chemical Society.

[35]  M. Parrinello,et al.  The unfolded ensemble and folding mechanism of the C-terminal GB1 beta-hairpin. , 2008, Journal of the American Chemical Society.

[36]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[37]  A. Kitao,et al.  Erratum to 'Multi-Scale Free Energy Landscape calculation method by combination of coarse-grained and all-atom models' [Chem. Phys. Lett. 503 (2011) 145-152] , 2011 .

[38]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[39]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[40]  V. Muñoz,et al.  Folding dynamics and mechanism of β-hairpin formation , 1997, Nature.

[41]  B. Matthews,et al.  Structure of bacteriophage T4 lysozyme refined at 1.7 A resolution. , 1987, Journal of molecular biology.

[42]  Ryuhei Harada,et al.  Exploring the folding free energy landscape of a β-hairpin miniprotein, chignolin, using multiscale free energy landscape calculation method. , 2011, The journal of physical chemistry. B.

[43]  B. Matthews,et al.  Protein flexibility and adaptability seen in 25 crystal forms of T4 lysozyme. , 1995, Journal of molecular biology.

[44]  C. Dellago,et al.  Transition path sampling and the calculation of rate constants , 1998 .

[45]  Teresa Head-Gordon,et al.  Instantaneous normal modes as an unforced reaction coordinate for protein conformational transitions. , 2010, Biophysical journal.

[46]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[47]  C. Jarzynski Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[48]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[49]  M. Seibert,et al.  Reproducible polypeptide folding and structure prediction using molecular dynamics simulations. , 2005, Journal of molecular biology.

[50]  K. Schulten,et al.  Steered molecular dynamics and mechanical functions of proteins. , 2001, Current opinion in structural biology.

[51]  H. Berendsen,et al.  Domain motions in bacteriophage T4 lysozyme: A comparison between molecular dynamics and crystallographic data , 1998, Proteins.

[52]  Jan F. Prins,et al.  SMD: visual steering of molecular dynamics for protein design , 1996 .

[53]  A. Kitao Transform and relax sampling for highly anisotropic systems: application to protein domain motion and folding. , 2011, The Journal of chemical physics.

[54]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[55]  A. Kidera,et al.  Multicanonical Ensemble Generated by Molecular Dynamics Simulation for Enhanced Conformational Sampling of Peptides , 1997 .

[56]  K Schulten,et al.  Stability and dynamics of G-actin: back-door water diffusion and behavior of a subdomain 3/4 loop. , 1997, Biophysical journal.