TY - JOUR
T1 - FF12MC
T2 - A revised AMBER forcefield and new protein simulation protocol
AU - Pang, Yuan Ping
N1 - Funding Information:
Yuan-Ping Pang acknowledges the support of this work from the US Defense Advanced Research Projects Agency (DAAD19-01–1-0322), the US Army Medical Research Material Command (W81XWH-04–2-0001), the US Army Research Office (DAAD19-03–1-0318, W911NF-09–1-0095, and W911NF-16–1-0264), the US Department of Defense High Performance Computing Modernization Office, and the Mayo Foundation for Medical Education and Research. The author remains in debt to the late Professor Peter A. Kollman for teaching him the minimalist approach to forcefield development during a one-year sabbatical in the Kollman group 1994–1995. The author is also in debt to the late Professor Shneior Lifson for a stimulating discussion on forcefield development during his visit to the Weizmann Institute of Science, Rehovot, Israel in 1996. The author thanks four anonymous reviewers for their comments and suggestions. The contents of this article are the sole responsibility of the author and do not necessarily represent the official views of the funders.
Publisher Copyright:
© 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.
PY - 2016/10/1
Y1 - 2016/10/1
N2 - Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened CH bonds, (ii) removal of torsions involving a nonperipheral sp3 atom, and (iii) reduced 1–4 interaction scaling factors of torsions ϕ and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric–isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left- and right-handed configurations for C14–C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp-cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general-purpose AMBER forcefield FF14SB locks the C14–C38 bond to the right-handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric–isothermal MD simulations 2–4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490–1516.
AB - Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened CH bonds, (ii) removal of torsions involving a nonperipheral sp3 atom, and (iii) reduced 1–4 interaction scaling factors of torsions ϕ and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric–isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left- and right-handed configurations for C14–C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp-cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general-purpose AMBER forcefield FF14SB locks the C14–C38 bond to the right-handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric–isothermal MD simulations 2–4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490–1516.
KW - BPTI
KW - CLN025
KW - Trp-cage
KW - chignolin
KW - force field
KW - molecular dynamics simulation
KW - protein dynamics
KW - protein folding
KW - protein simulation
KW - protein structure refinement
UR - http://www.scopus.com/inward/record.url?scp=84978937989&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84978937989&partnerID=8YFLogxK
U2 - 10.1002/prot.25094
DO - 10.1002/prot.25094
M3 - Article
C2 - 27348292
AN - SCOPUS:84978937989
SN - 0887-3585
VL - 84
SP - 1490
EP - 1516
JO - Proteins: Structure, Function and Bioinformatics
JF - Proteins: Structure, Function and Bioinformatics
IS - 10
ER -