+ 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

Efficiency of alchemical free energy simulations. I. A practical comparison of the exponential formula, thermodynamic integration, and Bennett's acceptance ratio method



Efficiency of alchemical free energy simulations. I. A practical comparison of the exponential formula, thermodynamic integration, and Bennett's acceptance ratio method



Journal of Computational Chemistry 32(7): 1303-1319



We investigate the relative efficiency of thermodynamic integration, three variants of the exponential formula, also referred to as thermodynamic perturbation, and Bennett's acceptance ratio method to compute relative and absolute solvation free energy differences. Our primary goal is the development of efficient protocols that are robust in practice. We focus on minimizing the number of unphysical intermediate states (λ-states) required for the computation of accurate and precise free energy differences. Several indicators are presented which help decide when additional λ-states are necessary. In all tests Bennett's acceptance ratio method required the least number of λ-states, closely followed by the "double-wide" variant of the exponential formula. Use of the exponential formula in only strict "forward" or "backward" mode was not found to be competitive. Similarly, the performance of thermodynamic integration in terms of efficiency was rather poor. We show that this is caused by the use of the trapezoidal rule as method of numerical quadrature. A systematic study focusing on the optimization of thermodynamic integration is presented in a companion paper.

Please choose payment method:






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

Accession: 052906387

Download citation: RISBibTeXText

PMID: 21425288

DOI: 10.1002/jcc.21713


Related references

Comparison of efficiency and bias of free energies computed by exponential averaging, the Bennett acceptance ratio, and thermodynamic integration. Journal of Chemical Physics 122(14): 144107, 2005

Efficiency of alchemical free energy simulations. II. Improvements for thermodynamic integration. Journal of Computational Chemistry 32(7): 1320-1333, 2011

Comparison of thermodynamic integration and Bennett acceptance ratio for calculating relative protein-ligand binding free energies. Journal of Computational Chemistry 34(12): 1024-1034, 2013

Ensemble variance in free energy calculations by thermodynamic integration: theory, optimal "Alchemical" path, and practical solutions. Journal of Computational Chemistry 25(7): 985-993, 2004

Alchemical Free Energy Differences in Flexible Molecules from Thermodynamic Integration or Free Energy Perturbation Combined with Driven Adiabatic Dynamics. Journal of Chemical Theory and Computation 8(10): 3504-3512, 2012

Efficient computation of free energy of crystal phases due to external potentials by error-biased Bennett acceptance ratio method. Journal of Chemical Physics 132(8): 084101, 2010

Thermodynamic reversibility and statistical convergence of alchemical free energy simulations. Biophysical Journal 72(2 Part 2): A397, 1997

Bennett acceptance ratio method to calculate the binding free energy of BACE1 inhibitors: Theoretical model and design of new ligands of the enzyme. Chemical Biology and Drug Design 93(6): 1117-1128, 2019

On the use of Bennett's acceptance ratio method in multi-canonical-type simulations. Journal of Chemical Physics 120(7): 3066-3074, 2004

On relation between the free-energy perturbation and Bennett's acceptance ratio methods: Tracing the influence of the energy gap. Journal of Chemical Physics 132(19): 194104, 2010

Alchemical Free-Energy Calculations by Multiple-Replica λ-Dynamics: The Conveyor Belt Thermodynamic Integration Scheme. Journal of Chemical Theory and Computation 15(4): 2392-2419, 2019

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

Chaperoned alchemical free energy simulations: a general method for QM, MM, and QM/MM potentials. Journal of Chemical Physics 120(20): 9450-9453, 2004

Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis. Journal of Physical Chemistry. B 122(41): 9435-9442, 2018

Enhanced Sampling in Free Energy Calculations: Combining SGLD with the Bennett's Acceptance Ratio and Enveloping Distribution Sampling Methods. Journal of Chemical Theory and Computation 8(10): 3650-3662, 2012