Enhanced sampling molecular dynamics simulation captures experimentally suggested intermediate and unfolded states in the folding pathway of Trp-cage miniprotein.

The ability of molecular dynamics simulation to capturing the transient states within the folding pathway of protein is important to the understanding of protein folding mechanism. In the present study, the integrated-tempering-sampling molecular dynamics (ITS-MD) simulation was performed to investigate the transient states including intermediate and unfolded ones in the folding pathway of a miniprotein, Trp-cage. Three force fields (FF03, FF99SB, and FF96) were tested, and both intermediate and unfolded states with their characteristics in good agreement with experiments were observed during the simulations, which supports the hypothesis that observable intermediates might present in the folding pathway of small polypeptides. In addition, it was demonstrated that FF03 force field as combined with ITS-MD is in overall a more proper force field than the others in reproducing experimentally recorded properties in UVRS, ECD, and NMR, Photo-CIDNP NMR, and IR T-jump experiments, and the folding∕unfolding thermodynamics parameters, such as ΔG(U), ΔC(p), and ΔH(U) (T(m)). In summary, the present study showed that using suitable force field and energy sampling method, molecular dynamics simulation could capture the transient states within the folding pathway of protein which are consistent with the experimental measurements, and thus provide information of protein folding mechanism and thermodynamics.

[1]  Seokmin Shin,et al.  Folding simulations with novel conformational search method. , 2007, The Journal of chemical physics.

[2]  R. Berry,et al.  Investigations into sequence and conformational dependence of backbone entropy, inter-basin dynamics and the Flory isolated-pair hypothesis for peptides. , 2003, Journal of molecular biology.

[3]  J. H. Zhang,et al.  Simulation of the thermodynamics of folding and unfolding of the Trp-cage mini-protein TC5b using different combinations of force fields and solvation models , 2010 .

[4]  Qiang Shao,et al.  Molecular dynamics and ion mobility spectrometry study of model β-hairpin peptide, trpzip1. , 2011, The journal of physical chemistry. A.

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

[6]  H. Roder,et al.  Kinetic role of early intermediates in protein folding. , 1997, Current opinion in structural biology.

[7]  Qiang Shao,et al.  Comparison between integrated and parallel tempering methods in enhanced sampling simulations. , 2009, The Journal of chemical physics.

[8]  Lijiang Yang,et al.  Thermodynamics and kinetics simulations of multi-time-scale processes for complex systems , 2008 .

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

[10]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[11]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[12]  George I Makhatadze,et al.  Unfolding thermodynamics of Trp-cage, a 20 residue miniprotein, studied by differential scanning calorimetry and circular dichroism spectroscopy. , 2007, Biochemistry.

[13]  S. Buldyrev,et al.  Folding Trp-cage to NMR resolution native structure using a coarse-grained protein model. , 2004, Biophysical journal.

[14]  Petra Rovó,et al.  Cooperativity network of Trp‐cage miniproteins: probing salt‐bridges , 2011, Journal of peptide science : an official publication of the European Peptide Society.

[15]  Qiang Shao,et al.  The relative helix and hydrogen bond stability in the B domain of protein A as revealed by integrated tempering sampling molecular dynamics simulation. , 2011, The Journal of chemical physics.

[16]  A. Fersht,et al.  Protein Folding and Unfolding at Atomic Resolution , 2002, Cell.

[17]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[18]  M. Zacharias,et al.  Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations , 2009, International journal of molecular sciences.

[19]  W. Eaton,et al.  Characterizing the unfolded states of proteins using single-molecule FRET spectroscopy and molecular simulations , 2007, Proceedings of the National Academy of Sciences.

[20]  Qiang Shao,et al.  Structure change of β-hairpin induced by turn optimization: an enhanced sampling molecular dynamics simulation study. , 2011, The Journal of chemical physics.

[21]  B. Berne,et al.  Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). , 2011, The journal of physical chemistry. B.

[22]  Asim Okur,et al.  Improving the description of salt bridge strength and geometry in a Generalized Born model. , 2011, Journal of molecular graphics & modelling.

[23]  R. Zhou Trp-cage: Folding free energy landscape in explicit water , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Qiang Shao,et al.  Effects of turn stability and side-chain hydrophobicity on the folding of β-structures. , 2010, Journal of molecular biology.

[25]  K. H. Mok,et al.  A pre-existing hydrophobic collapse in the unfolded state of an ultrafast folding protein , 2007, Nature.

[26]  Klaus Schulten,et al.  Challenges in protein-folding simulations , 2010 .

[27]  Jeffery G Saven,et al.  Ultrafast folding of a computationally designed Trp-cage mutant: Trp2-cage. , 2006, The journal of physical chemistry. B.

[28]  William Swope,et al.  Understanding folding and design: Replica-exchange simulations of ``Trp-cage'' miniproteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[29]  A. Garcia,et al.  Microsecond simulations of the folding/unfolding thermodynamics of the Trp‐cage miniprotein , 2010, Proteins.

[30]  A. Garcia,et al.  Computing the stability diagram of the Trp-cage miniprotein , 2008, Proceedings of the National Academy of Sciences.

[31]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

[32]  Asim Okur,et al.  Evaluation of Salt Bridge Structure and Energetics in Peptides Using Explicit, Implicit, and Hybrid Solvation Models. , 2008, Journal of chemical theory and computation.

[33]  F. Gai,et al.  Achieving secondary structural resolution in kinetic measurements of protein folding: a case study of the folding mechanism of Trp-cage. , 2011, Angewandte Chemie.

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

[35]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[36]  Haruki Nakamura,et al.  Energy landscape of a peptide consisting of α‐helix, 310‐helix, β‐turn, β‐hairpin, and other disordered conformations , 2001 .

[37]  Alessandro Laio,et al.  A Kinetic Model of Trp-Cage Folding from Multiple Biased Molecular Dynamics Simulations , 2009, PLoS Comput. Biol..

[38]  S. Asher,et al.  UV-resonance raman thermal unfolding study of Trp-cage shows that it is not a simple two-state miniprotein. , 2005, Journal of the American Chemical Society.

[39]  Kevin W Plaxco,et al.  Protein folding: Defining a “standard” set of experimental conditions and a preliminary kinetic data set of two‐state proteins , 2005, Protein science : a publication of the Protein Society.

[40]  Y. Duan,et al.  Characterizing the rate-limiting step of Trp-cage folding by all-atom molecular dynamics simulations , 2004 .

[41]  Lauren Wickstrom,et al.  Secondary structure bias in generalized Born solvent models: comparison of conformational ensembles and free energy of solvent polarization from explicit and implicit solvation. , 2007, The journal of physical chemistry. B.

[42]  A. Roitberg,et al.  Smaller and faster: the 20-residue Trp-cage protein folds in 4 micros. , 2002, Journal of the American Chemical Society.

[43]  P. Bolhuis,et al.  Sampling the multiple folding mechanisms of Trp-cage in explicit solvent , 2006, Proceedings of the National Academy of Sciences.

[44]  Xu Zhang,et al.  Dynamics and cooperativity of Trp-cage folding. , 2008, Archives of biochemistry and biophysics.

[45]  Sören Doose,et al.  A microscopic view of miniprotein folding: enhanced folding efficiency through formation of an intermediate. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[46]  Haruki Nakamura,et al.  Peptide free‐energy profile is strongly dependent on the force field: Comparison of C96 and AMBER95 , 2000 .

[47]  Y. Gao An integrate-over-temperature approach for enhanced sampling. , 2008, The Journal of chemical physics.

[48]  Qiang Shao,et al.  A test of implicit solvent models on the folding simulation of the GB1 peptide. , 2009, The Journal of chemical physics.

[49]  Sheena E Radford,et al.  Intermediates: ubiquitous species on folding energy landscapes? , 2007, Current opinion in structural biology.

[50]  Haruki Nakamura,et al.  Conformational transition states of a β‐hairpin peptide between the ordered and disordered conformations in explicit water , 2002, Protein science : a publication of the Protein Society.

[51]  Valerie Daggett,et al.  The complete folding pathway of a protein from nanoseconds to microseconds , 2003, Nature.

[52]  Temperature Dependence of Hydrogen-Bond Stability in β-Hairpin Structures , 2010 .