Limits...
Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration.

Fajer M, Hamelberg D, McCammon JA - J Chem Theory Comput (2008)

Bottom Line: Accelerated molecular dynamics (AMD) is an efficient strategy for accelerating the sampling of molecular dynamics simulations, and observable quantities such as free energies derived on the biased AMD potential can be reweighted to yield results consistent with the original, unmodified potential.We propose a replica exchange of various degrees of acceleration (REXAMD) to retain good statistics while achieving enhanced sampling.The REXAMD method is validated and benchmarked on two simple gas-phase model systems, and two different strategies for computing reweighted averages over a simulation are compared.

View Article: PubMed Central - PubMed

Affiliation: Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California 92093-0365, Center for Theoretical Biological Physics, University of California at San Diego, La Jolla, California 92039-0365, Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-4098, Department of Pharmacology, University of California at San Diego, La Jolla, California 92093-0365, and Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California 92093-0365.

ABSTRACT
Accelerated molecular dynamics (AMD) is an efficient strategy for accelerating the sampling of molecular dynamics simulations, and observable quantities such as free energies derived on the biased AMD potential can be reweighted to yield results consistent with the original, unmodified potential. In conventional AMD the reweighting procedure has an inherent statistical problem in systems with large acceleration, where the points with the largest biases will dominate the reweighted result and reduce the effective number of data points. We propose a replica exchange of various degrees of acceleration (REXAMD) to retain good statistics while achieving enhanced sampling. The REXAMD method is validated and benchmarked on two simple gas-phase model systems, and two different strategies for computing reweighted averages over a simulation are compared.

No MeSH data available.


Block average of the MSA thermodynamic integration results when using the reweighted periods strategy. The symbols show the average value of each simulation type, and the shaded region shows the standard error from the ten duplicate runs.
© Copyright Policy - open-access - ccc-price
Related In: Results  -  Collection

License
getmorefigures.php?uid=PMC2651661&req=5

fig2: Block average of the MSA thermodynamic integration results when using the reweighted periods strategy. The symbols show the average value of each simulation type, and the shaded region shows the standard error from the ten duplicate runs.

Mentions: In MSA both the REXREG and REXAMDt simulations were able to efficiently and exhaustively explore the conformational space (data not shown), and the replica mixing was quite efficient (Table 1) within the 8 ns runs. The exhaustive sampling resulted in converged ΔG values within the first ns of the REXREG and REXAMDt simulations (Figure 2). The ΔG results from the entire 8 ns are summarized in Table 2. Recall that “MSA REXAMDt” refers to results taken from the nonaccelerated state and “MSA REXAMDt-rw” refers to the reweighted results of the accelerated state.


Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration.

Fajer M, Hamelberg D, McCammon JA - J Chem Theory Comput (2008)

Block average of the MSA thermodynamic integration results when using the reweighted periods strategy. The symbols show the average value of each simulation type, and the shaded region shows the standard error from the ten duplicate runs.
© Copyright Policy - open-access - ccc-price
Related In: Results  -  Collection

License
Show All Figures
getmorefigures.php?uid=PMC2651661&req=5

fig2: Block average of the MSA thermodynamic integration results when using the reweighted periods strategy. The symbols show the average value of each simulation type, and the shaded region shows the standard error from the ten duplicate runs.
Mentions: In MSA both the REXREG and REXAMDt simulations were able to efficiently and exhaustively explore the conformational space (data not shown), and the replica mixing was quite efficient (Table 1) within the 8 ns runs. The exhaustive sampling resulted in converged ΔG values within the first ns of the REXREG and REXAMDt simulations (Figure 2). The ΔG results from the entire 8 ns are summarized in Table 2. Recall that “MSA REXAMDt” refers to results taken from the nonaccelerated state and “MSA REXAMDt-rw” refers to the reweighted results of the accelerated state.

Bottom Line: Accelerated molecular dynamics (AMD) is an efficient strategy for accelerating the sampling of molecular dynamics simulations, and observable quantities such as free energies derived on the biased AMD potential can be reweighted to yield results consistent with the original, unmodified potential.We propose a replica exchange of various degrees of acceleration (REXAMD) to retain good statistics while achieving enhanced sampling.The REXAMD method is validated and benchmarked on two simple gas-phase model systems, and two different strategies for computing reweighted averages over a simulation are compared.

View Article: PubMed Central - PubMed

Affiliation: Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California 92093-0365, Center for Theoretical Biological Physics, University of California at San Diego, La Jolla, California 92039-0365, Department of Chemistry, Georgia State University, Atlanta, Georgia 30302-4098, Department of Pharmacology, University of California at San Diego, La Jolla, California 92093-0365, and Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California 92093-0365.

ABSTRACT
Accelerated molecular dynamics (AMD) is an efficient strategy for accelerating the sampling of molecular dynamics simulations, and observable quantities such as free energies derived on the biased AMD potential can be reweighted to yield results consistent with the original, unmodified potential. In conventional AMD the reweighting procedure has an inherent statistical problem in systems with large acceleration, where the points with the largest biases will dominate the reweighted result and reduce the effective number of data points. We propose a replica exchange of various degrees of acceleration (REXAMD) to retain good statistics while achieving enhanced sampling. The REXAMD method is validated and benchmarked on two simple gas-phase model systems, and two different strategies for computing reweighted averages over a simulation are compared.

No MeSH data available.