Reconstructing the Equilibrium Boltzmann Distribution from Well-Tempered Metadynamics

被引:271
作者
Bonomi, M. [1 ]
Barducci, A. [1 ]
Parrinello, M. [1 ]
机构
[1] Swiss Fed Inst Technol, Dept Chem & Appl Biosci, CH-6900 Lugano, Switzerland
关键词
molecular dynamics; metadynamics; Boltzmann distribution; reweighting; umbrella sampling; MOLECULAR-DYNAMICS SIMULATION; FREE-ENERGY CALCULATIONS; CONFORMATIONAL TRANSITIONS; PHASE-TRANSITIONS; RATE CONSTANTS; BETA-HAIRPIN; MEAN FORCE; PROTEIN; SYSTEMS; ALGORITHM;
D O I
10.1002/jcc.21305
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
Metadynamics is a widely used and successful method for reconstructing the free-energy surface of complex systems as a function of a small number of suitably chosen collective variables. This is achieved by biasing the dynamics of the system. The bias acting on the collective variables distorts the probability distribution of the other variables. Here we present a simple reweighting algorithm for recovering the unbiased probability distribution of any variable from a well-tempered metadynamics simulation. We show the efficiency of the reweighting procedure by reconstructing the distribution of the four backbone dihedral angles of alanine dipeptide from two and even one dimensional metadynamics simulation. (C) 2009 Wiley Periodicals, Inc. J Comput Chem 30: 1615-1621, 2009
引用
收藏
页码:1615 / 1621
页数:7
相关论文
共 62 条
[1]   Sampling rare switching events in biochemical networks [J].
Allen, RJ ;
Warren, PB ;
ten Wolde, PR .
PHYSICAL REVIEW LETTERS, 2005, 94 (01)
[2]   Calculation of conformational transitions and barriers in solvated systems: Application to the alanine dipeptide in water [J].
Apostolakis, J ;
Ferrara, P ;
Caflisch, A .
JOURNAL OF CHEMICAL PHYSICS, 1999, 110 (04) :2099-2108
[3]   Metadynamics simulation of prion protein:: β-structure stability and the early stages of misfolding [J].
Barducci, A ;
Chelli, R ;
Procacci, P ;
Schettino, V ;
Gervasio, FL ;
Parrinello, M .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 2006, 128 (08) :2705-2710
[4]   Well-tempered metadynamics: A smoothly converging and tunable free-energy method [J].
Barducci, Alessandro ;
Bussi, Giovanni ;
Parrinello, Michele .
PHYSICAL REVIEW LETTERS, 2008, 100 (02)
[5]  
Bartels C, 1997, J COMPUT CHEM, V18, P1450, DOI 10.1002/(SICI)1096-987X(199709)18:12<1450::AID-JCC3>3.0.CO
[6]  
2-I
[7]   CALCULATION OF THE RELATIVE CHANGE IN BINDING FREE-ENERGY OF A PROTEIN-INHIBITOR COMPLEX [J].
BASH, PA ;
SINGH, UC ;
BROWN, FK ;
LANGRIDGE, R ;
KOLLMAN, PA .
SCIENCE, 1987, 235 (4788) :574-576
[8]   Metadynamics simulations of the high-pressure phases of silicon employing a high-dimensional neural network potential [J].
Behler, Joerg ;
Martonak, Roman ;
Donadio, Davide ;
Parrinello, Michele .
PHYSICAL REVIEW LETTERS, 2008, 100 (18)
[9]   Protein Conformational Transitions: The Closure Mechanism of a Kinase Explored by Atomistic Simulations [J].
Berteotti, Anna ;
Cavalli, Andrea ;
Branduardi, Davide ;
Gervasio, Francesco Luigi ;
Recanatini, Maurizio ;
Parrinello, Michele .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 2009, 131 (01) :244-250
[10]   Transition path sampling: Throwing ropes over rough mountain passes, in the dark [J].
Bolhuis, PG ;
Chandler, D ;
Dellago, C ;
Geissler, PL .
ANNUAL REVIEW OF PHYSICAL CHEMISTRY, 2002, 53 :291-318