Peptide folding kinetics from replica exchange molecular dynamics.
We show how accurate kinetic information, such as the rates of protein folding and unfolding, can be extracted from replica-exchange molecular dynamics (REMD) simulations. From the brief and continuous trajectory segments between replica exchanges, we estimate short-time propagators in conformation space and use them to construct a master equation. For a helical peptide in explicit water, we determine the rates of transitions both locally between microscopic conformational states and globally for folding and unfolding. We show that accurate rates in the approximately 1(100 ns) to approximately 1(1 ns) range can be obtained from REMD with exchange times of 5 ps , in excellent agreement with results from long equilibrium molecular dynamics.
Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Building 5, Bethesda, Maryland 20892-0520, USA.
Nicolae-Viorel
NV
Buchete
Gerhard
G
Hummer
D001703:Biophysics / Q000379:methods*; D003198:Computer Simulation; D007700:Kinetics; D008390:Markov Chains; D008956:Models, Chemical; D008958:Models, Molecular; D015233:Models, Statistical; D008968:Molecular Conformation; D010455:Peptides / Q000737:chemistry; D011487:Protein Conformation; D011489:Protein Denaturation; D017510:Protein Folding; D013816:Thermodynamics; D013997:Time Factors
D010455:Peptides