Estimating kinetic rates from accelerated molecular dynamics simulations: Alanine dipeptide in explicit solvent as a case study

被引:41
作者
De Oliveira, Cesar Augusto F. [1 ]
Hamelberg, Donald
McCammon, J. Andrew
机构
[1] Univ Calif San Diego, Howard Hughes Med Inst, Ctr Theoret Biol Phys, Dept Chem & Biochem, La Jolla, CA 92093 USA
[2] Univ Calif San Diego, Dept Pharmacol, La Jolla, CA 92093 USA
基金
美国国家卫生研究院; 美国国家科学基金会;
关键词
D O I
10.1063/1.2794763
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
Molecular dynamics (MD) simulation is the standard computational technique used to obtain information on the time evolution of the conformations of proteins and many other molecular systems. However, for most biological systems of interest, the time scale for slow conformational transitions is still inaccessible to standard MD simulations. Several sampling methods have been proposed to address this issue, including the accelerated molecular dynamics method. In this work, we study the extent of sampling of the phi/psi space of alanine dipeptide in explicit water using accelerated molecular dynamics and present a framework to recover the correct kinetic rate constant for the helix to beta-strand transition. We show that the accelerated MD can drastically enhance the sampling of the phi/psi conformational phase space when compared to normal MD. In addition, the free energy density plots of the phi/psi space show that all minima regions are accurately sampled and the canonical distribution is recovered. Moreover, the kinetic rate constant for the helix to beta-strand transition is accurately estimated from these simulations by relating the diffusion coefficient to the local energetic roughness of the energy landscape. Surprisingly, even for such a low barrier transition, it is difficult to obtain enough transitions to accurately estimate the rate constant when one uses normal MD.
引用
收藏
页数:8
相关论文
共 52 条
[1]   Molecular dynamics: Survey of methods for simulating the activity of proteins [J].
Adcock, Stewart A. ;
McCammon, J. Andrew .
CHEMICAL REVIEWS, 2006, 106 (05) :1589-1615
[2]  
Allen M. P., 2009, Computer Simulation of Liquids
[3]   MOLECULAR-DYNAMICS WITH COUPLING TO AN EXTERNAL BATH [J].
BERENDSEN, HJC ;
POSTMA, JPM ;
VANGUNSTEREN, WF ;
DINOLA, A ;
HAAK, JR .
JOURNAL OF CHEMICAL PHYSICS, 1984, 81 (08) :3684-3690
[4]   CONFORMATIONAL SAMPLING USING HIGH-TEMPERATURE MOLECULAR-DYNAMICS [J].
BRUCCOLERI, RE ;
KARPLUS, M .
BIOPOLYMERS, 1990, 29 (14) :1847-1862
[5]   Equilibrium free energies from nonequilibrium metadynamics [J].
Bussi, G ;
Laio, A ;
Parrinello, M .
PHYSICAL REVIEW LETTERS, 2006, 96 (09)
[6]   Free-energy landscape for β hairpin folding from combined parallel tempering and metadynamics [J].
Bussi, Giovanni ;
Gervasio, Francesco Luigi ;
Laio, Alessandro ;
Parrinello, Michele .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 2006, 128 (41) :13435-13441
[7]  
CASE D, UNPUB
[8]   Long-time conformational transitions of alanine dipeptide in aqueous solution: Continuous and discrete-state kinetic models [J].
Chekmarev, DS ;
Ishida, T ;
Levy, RM .
JOURNAL OF PHYSICAL CHEMISTRY B, 2004, 108 (50) :19487-19495
[9]   On the application of accelerated molecular dynamics to liquid water simulations [J].
de Oliveira, Cesar Augusto F. ;
Hamelberg, Donald ;
McCammon, J. Andrew .
JOURNAL OF PHYSICAL CHEMISTRY B, 2006, 110 (45) :22695-22701
[10]   ATOMIC LEVEL SIMULATIONS ON A MILLION PARTICLES - THE CELL MULTIPOLE METHOD FOR COULOMB AND LONDON NONBOND INTERACTIONS [J].
DING, HQ ;
KARASAWA, N ;
GODDARD, WA .
JOURNAL OF CHEMICAL PHYSICS, 1992, 97 (06) :4309-4315