+ Site Statistics
+ Search Articles
+ PDF Full Text Service
How our service works
Request PDF Full Text
+ Follow Us
Follow on Facebook
Follow on Twitter
Follow on LinkedIn
+ Subscribe to Site Feeds
Most Shared
PDF Full Text
+ Translate
+ Recently Requested

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation



Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation



Journal of Chemical Theory and Computation 10(7): 2677-2689



Accelerated molecular dynamics (aMD) simulations greatly improve the efficiency of conventional molecular dynamics (cMD) for sampling biomolecular conformations, but they require proper reweighting for free energy calculation. In this work, we systematically compare the accuracy of different reweighting algorithms including the exponential average, Maclaurin series, and cumulant expansion on three model systems: alanine dipeptide, chignolin, and Trp-cage. Exponential average reweighting can recover the original free energy profiles easily only when the distribution of the boost potential is narrow (e.g., the range ≤20kBT) as found in dihedral-boost aMD simulation of alanine dipeptide. In dual-boost aMD simulations of the studied systems, exponential average generally leads to high energetic fluctuations, largely due to the fact that the Boltzmann reweighting factors are dominated by a very few high boost potential frames. In comparison, reweighting based on Maclaurin series expansion (equivalent to cumulant expansion on the first order) greatly suppresses the energetic noise but often gives incorrect energy minimum positions and significant errors at the energy barriers (∼2-3kBT). Finally, reweighting using cumulant expansion to the second order is able to recover the most accurate free energy profiles within statistical errors of ∼kBT, particularly when the distribution of the boost potential exhibits low anharmonicity (i.e., near-Gaussian distribution), and should be of wide applicability. A toolkit of Python scripts for aMD reweighting "PyReweighting" is distributed free of charge at http://mccammon.ucsd.edu/computing/amdReweighting/.

Please choose payment method:






(PDF emailed within 0-6 h: $19.90)

Accession: 053742757

Download citation: RISBibTeXText

PMID: 25061441

DOI: 10.1021/ct500090q


Related references

Free energy calculations from adaptive molecular dynamics simulations with adiabatic reweighting. Journal of Chemical Physics 140(10): 104108, 2014

Replica Exchange Gaussian Accelerated Molecular Dynamics: Improved Enhanced Sampling and Free Energy Calculation. Journal of Chemical Theory and Computation 14(4): 1853-1864, 2018

Calculation of the entropy and free energy of peptides by molecular dynamics simulations using the hypothetical scanning molecular dynamics method. Journal of Chemical Physics 125(2): 24905, 2006

Free energy calculation of permeant-membrane interactions using molecular dynamics simulations. Methods in Molecular Biology 926: 189-202, 2012

Gaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling and Free Energy Calculation. Journal of Chemical Theory and Computation 11(8): 3584-3595, 2015

Gaussian Accelerated Molecular Dynamics: Unconstrained Enhanced Sampling and Free Energy Calculation. Journal of Chemical Theory and Computation 11(8): 3584-3595, 2015

Sampling enhancement for the quantum mechanical potential based molecular dynamics simulations: a general algorithm and its extension for free energy calculation on rugged energy surface. Journal of Chemical Physics 126(11): 114104, 2007

Insight to the molecular mechanisms of the osmolyte effects on Mycobacterium tuberculosis pyrazinamidase stability using experimental studies, molecular dynamics simulations, and free energy calculation. International Journal of Mycobacteriology 7(3): 268-274, 2018

Molecular mechanism of R-bicalutamide switching from androgen receptor antagonist to agonist induced by amino acid mutations using molecular dynamics simulations and free energy calculation. Journal of Computer-Aided Molecular Design 30(12): 1189-1200, 2016

On the improvement of free-energy calculation from steered molecular dynamics simulations using adaptive stochastic perturbation protocols. Plos one 9(9): E101810, 2014

On the improvement of free-energy calculation from steered molecular dynamics simulations using adaptive stochastic perturbation protocols. Plos one 9(9): E101810, 2014

Multiscale Free Energy Simulations: An Efficient Method for Connecting Classical MD Simulations to QM or QM/MM Free Energies Using Non-Boltzmann Bennett Reweighting Schemes. Journal of Chemical Theory and Computation 10(4): 1406-1419, 2014

An efficient protocol for obtaining accurate hydration free energies using quantum chemistry and reweighting from molecular dynamics simulations. Bioorganic and Medicinal Chemistry 24(20): 4988-4997, 2016

Computational study on the interaction between CCR5 and HIV-1 entry inhibitor maraviroc: insight from accelerated molecular dynamics simulation and free energy calculation. Physical Chemistry Chemical Physics 16(44): 24332-8, 2014

Calculation of the standard binding free energy of sparsomycin to the ribosomal peptidyl-transferase P-site using molecular dynamics simulations with restraining potentials. Journal of Molecular Recognition 23(2): 128-141, 2010