Simulating the conformational dynamics of biomolecules is extremely difficult due to the rugged nature of their free energy landscapes and multiple long-lived, or metastable, states. Generalized ensemble (GE) algorithms, which have become popular in recent years, attempt to facilitate crossing between states at low temperatures by inducing a random walk in temperature space. Enthalpic barriers may be crossed more easily at high temperatures; however, entropic barriers will become more significant. This poses a problem because the dominant barriers to conformational change are entropic for many biological systems, such as the short RNA hairpin studied here. We present a new efficient algorithm for conformational sampling, called the adaptive seeding method (ASM), which uses nonequilibrium GE simulations to identify the metastable states, and seeds short simulations at constant temperature from each of them to quantitatively determine their equilibrium populations. Thus, the ASM takes advantage of the broad sampling possible with GE algorithms but generally crosses entropic barriers more efficiently during the seeding simulations at low temperature. We show that only local equilibrium is necessary for ASM, so very short seeding simulations may be used. Moreover, the ASM may be used to recover equilibrium properties from existing datasets that failed to converge, and is well suited to running on modern computer clusters.
%0 Journal Article
%1 Huang2009GEEquilibration
%A Huang, Xuhui
%A Bowman, Gregory R.
%A Bacallado, Sergio
%A Pande, Vijay S.
%D 2009
%I National Academy of Sciences
%J Proceedings of the National Academy of Sciences
%K equilibrium generalized-ensemble molecular-dynamics nonequilibrium-to-equilibrium relaxation replica-exchange
%N 47
%P 19765--19769
%R 10.1073/pnas.0909088106
%T Rapid equilibrium sampling initiated from nonequilibrium data
%U http://www.pnas.org/content/106/47/19765
%V 106
%X Simulating the conformational dynamics of biomolecules is extremely difficult due to the rugged nature of their free energy landscapes and multiple long-lived, or metastable, states. Generalized ensemble (GE) algorithms, which have become popular in recent years, attempt to facilitate crossing between states at low temperatures by inducing a random walk in temperature space. Enthalpic barriers may be crossed more easily at high temperatures; however, entropic barriers will become more significant. This poses a problem because the dominant barriers to conformational change are entropic for many biological systems, such as the short RNA hairpin studied here. We present a new efficient algorithm for conformational sampling, called the adaptive seeding method (ASM), which uses nonequilibrium GE simulations to identify the metastable states, and seeds short simulations at constant temperature from each of them to quantitatively determine their equilibrium populations. Thus, the ASM takes advantage of the broad sampling possible with GE algorithms but generally crosses entropic barriers more efficiently during the seeding simulations at low temperature. We show that only local equilibrium is necessary for ASM, so very short seeding simulations may be used. Moreover, the ASM may be used to recover equilibrium properties from existing datasets that failed to converge, and is well suited to running on modern computer clusters.
@article{Huang2009GEEquilibration,
abstract = {Simulating the conformational dynamics of biomolecules is extremely difficult due to the rugged nature of their free energy landscapes and multiple long-lived, or metastable, states. Generalized ensemble (GE) algorithms, which have become popular in recent years, attempt to facilitate crossing between states at low temperatures by inducing a random walk in temperature space. Enthalpic barriers may be crossed more easily at high temperatures; however, entropic barriers will become more significant. This poses a problem because the dominant barriers to conformational change are entropic for many biological systems, such as the short RNA hairpin studied here. We present a new efficient algorithm for conformational sampling, called the adaptive seeding method (ASM), which uses nonequilibrium GE simulations to identify the metastable states, and seeds short simulations at constant temperature from each of them to quantitatively determine their equilibrium populations. Thus, the ASM takes advantage of the broad sampling possible with GE algorithms but generally crosses entropic barriers more efficiently during the seeding simulations at low temperature. We show that only local equilibrium is necessary for ASM, so very short seeding simulations may be used. Moreover, the ASM may be used to recover equilibrium properties from existing datasets that failed to converge, and is well suited to running on modern computer clusters.},
added-at = {2018-05-11T21:51:17.000+0200},
author = {Huang, Xuhui and Bowman, Gregory R. and Bacallado, Sergio and Pande, Vijay S.},
biburl = {https://www.bibsonomy.org/bibtex/2dca6281931491de5cecb8ebdcad940c5/salotz},
doi = {10.1073/pnas.0909088106},
eprint = {http://www.pnas.org/content/106/47/19765.full.pdf},
interhash = {4870ec5ef80b9b6edae27ce0026a4c35},
intrahash = {dca6281931491de5cecb8ebdcad940c5},
issn = {0027-8424},
journal = {Proceedings of the National Academy of Sciences},
keywords = {equilibrium generalized-ensemble molecular-dynamics nonequilibrium-to-equilibrium relaxation replica-exchange},
number = 47,
pages = {19765--19769},
publisher = {National Academy of Sciences},
timestamp = {2018-05-11T21:51:17.000+0200},
title = {Rapid equilibrium sampling initiated from nonequilibrium data},
url = {http://www.pnas.org/content/106/47/19765},
volume = 106,
year = 2009
}