Ensuring mixing efficiency of replica-exchange molecular dynamics simulations

被引:63
作者
Abraham, Mark J. [1 ]
Gready, Jill E. [1 ]
机构
[1] Australian Natl Univ, John Curtin Sch Med Res, Computat Prote Grp, Canberra, ACT 2601, Australia
关键词
D O I
10.1021/ct800016r
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
We address the question of constructing a protocol for replica-exchange molecular dynamics (REMD) simulations that make efficient use of the replica space, assess whether published applications are achieving such "mixing" efficiency, and provide a how-to guide to assist users to plan efficient REMD simulations. To address our first question, we introduce and discuss three metrics for assessing the number of replica-exchange attempts required to justify the use of a replica scheme and define a "transit number" as the lower bound for the length of an efficient simulation. Our literature survey of applications of REMD simulations of peptides in explicit solvent indicated that authors are not routinely reporting sufficient details of their simulation protocols to allow readers to make independent assessments of the impact of the method on their results, particularly whether mixing efficiency has been achieved. Necessary details include the expected or observed replica-exchange probability, together with the total number of exchange attempts, the exchange period, and estimates of the autocorrelation time of the potential energy. Our analysis of cases where the necessary information was reported suggests that in many of these simulations there are insufficient exchanges. attempted or an insufficiently long period between them to provide confidence that the simulation length justifies the size of the replica scheme. We suggest guidelines for designing REMD simulation protocols to ensure mixing efficiency. Two key recommendations are that the exchange period should in general be larger than 1 ps and the number of exchange attempts should be chosen to significantly exceed the transit number for the replica scheme.
引用
收藏
页码:1119 / 1128
页数:10
相关论文
共 67 条
[1]   A novel Hamiltonian replica exchange MD protocol to enhance protein conformational space sampling [J].
Affentranger, R ;
Tavernelli, I ;
Di Iorio, EE .
JOURNAL OF CHEMICAL THEORY AND COMPUTATION, 2006, 2 (02) :217-228
[2]   Amyloid β-protein monomer structure:: A computational and experimental study [J].
Baumketner, A ;
Bernstein, SL ;
Wyttenbach, T ;
Bitan, G ;
Teplow, DB ;
Bowers, MT ;
Shea, JE .
PROTEIN SCIENCE, 2006, 15 (03) :420-428
[3]   Structure of the 21-30 fragment of amyloid β-protein [J].
Baumketner, Andrij ;
Bernstein, Summer L. ;
Wyttenbach, Thomas ;
Lazo, Noel D. ;
Teplow, David B. ;
Bowers, Michael T. ;
Shea, Joan-Emma .
PROTEIN SCIENCE, 2006, 15 (06) :1239-1247
[4]   GROMACS - A MESSAGE-PASSING PARALLEL MOLECULAR-DYNAMICS IMPLEMENTATION [J].
BERENDSEN, HJC ;
VANDERSPOEL, D ;
VANDRUNEN, R .
COMPUTER PHYSICS COMMUNICATIONS, 1995, 91 (1-3) :43-56
[5]   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
[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]   Modified replica exchange simulation methods for local structure refinement [J].
Cheng, XL ;
Cui, GL ;
Hornak, V ;
Sinnnerling, C .
JOURNAL OF PHYSICAL CHEMISTRY B, 2005, 109 (16) :8220-8230
[8]   Multigraining: An algorithm for simultaneous fine-grained and coarse-grained simulation of molecular systems [J].
Christen, M ;
van Gunsteren, WF .
JOURNAL OF CHEMICAL PHYSICS, 2006, 124 (15)
[9]   A SMOOTH PARTICLE MESH EWALD METHOD [J].
ESSMANN, U ;
PERERA, L ;
BERKOWITZ, ML ;
DARDEN, T ;
LEE, H ;
PEDERSEN, LG .
JOURNAL OF CHEMICAL PHYSICS, 1995, 103 (19) :8577-8593
[10]   On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction [J].
Fukunishi, H ;
Watanabe, O ;
Takada, S .
JOURNAL OF CHEMICAL PHYSICS, 2002, 116 (20) :9058-9067