## Abstract Replica exchange accelerated molecular dynamics (REXAMD) is a method that enhances conformational sampling while retaining at least one replica on the original potential, thus avoiding the statistical problems of exponential reweighting. In this article, we study three methods that can
Efficient use of nonequilibrium measurement to estimate free energy differences for molecular systems
✍ Scribed by F. Marty Ytreberg; Daniel M. Zuckerman
- Publisher
- John Wiley and Sons
- Year
- 2004
- Tongue
- English
- Weight
- 802 KB
- Volume
- 25
- Category
- Article
- ISSN
- 0192-8651
No coin nor oath required. For personal study only.
✦ Synopsis
Abstract
A promising method for calculating free energy differences ΔF is to generate nonequilibrium data via “fast‐growth” simulations or by experiments—and then use Jarzynski's equality. However, a difficulty with using Jarzynski's equality is that ΔF estimates converge very slowly and unreliably due to the nonlinear nature of the calculation—thus requiring large, costly data sets. The purpose of the work presented here is to determine the best estimate for ΔF given a (finite) set of work values previously generated by simulation or experiment. Exploiting statistical properties of Jarzynski's equality, we present two fully automated analyses of nonequilibrium data from a toy model, and various simulated molecular systems. Both schemes remove at least several k~B~T of bias from ΔF estimates, compared to direct application of Jarzynski's equality, for modest sized data sets (100 work values), in all tested systems. Results from one of the new methods suggest that good estimates of ΔF can be obtained using 5–40‐fold less data than was previously possible. Extending previous work, the new results exploit the systematic behavior of bias due to finite sample size. A key innovation is better use of the more statistically reliable information available from the raw data. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1749–1759, 2004
📜 SIMILAR VOLUMES
## Abstract We present results showing the importance of appropriate treatment of atomic masses in molecular dynamics (MD)‐based single topology free‐energy perturbations (FEPs) on small molecule systems. The reversibility of gas phase simulations is significantly improved by scaling the atomic mas
## Abstract It is demonstrated that in a case where neither the proton nor the natural‐abundance ^13^C‐satellite spectra of a partially oriented molecule carry enough structural information, one can determine the entire molecular geometry by the combined use of several liquid crystals as solvents.