On relation between the free-energy perturbation and Bennett’s acceptance ratio methods: Tracing the influence of the energy gap

2010 ◽  
Vol 132 (19) ◽  
pp. 194104 ◽  
Author(s):  
V. B. Luzhkov
2021 ◽  
Vol 7 (2) ◽  
pp. 69-75
Author(s):  
S. P. Khanal ◽  
B. Poudel ◽  
R. P. Koirala ◽  
N. P. Adhikari

In the present work, we have used an alchemical approach for calculating solvation free energy of protonated lysine in water from molecular dynamics simulations. These approaches use a non-physical pathway between two end states in order to compute free energy difference from the set of simulations. The solute is modeled using bonded and non-bonded interactions described by OPLS-AA potential, while four different water models: TIP3P, SPC, SPC/E and TIP4P are used. The free energy of solvation of protonated lysine in water has been estimated using thermodynamic integration, free energy perturbation, and Bennett acceptance ratio methods at 310 K temperature. The contributions to the free energy due to van der Waals and electrostatics parameters are also separately computed. The estimated values of free energy of solvation using different methods are in well agreement with previously reported experimental value within 14 %.


Molecules ◽  
2019 ◽  
Vol 24 (4) ◽  
pp. 681 ◽  
Author(s):  
Fiona Kearns ◽  
Luke Warrensford ◽  
Stefan Boresch ◽  
H. Woodcock

Indirect (S)QM/MM free energy simulations (FES) are vital to efficiently incorporating sufficient sampling and accurate (QM) energetic evaluations when estimating free energies of practical/experimental interest. Connecting between levels of theory, i.e., calculating Δ A l o w → h i g h , remains to be the most challenging step within an indirect FES protocol. To improve calculations of Δ A l o w → h i g h , we must: (1) compare the performance of all FES methods currently available; and (2) compile and maintain datasets of Δ A l o w → h i g h calculated for a wide-variety of molecules so that future practitioners may replicate or improve upon the current state-of-the-art. Towards these two aims, we introduce a new dataset, “HiPen”, which tabulates Δ A g a s M M → 3 o b (the free energy associated with switching from an M M to an S C C − D F T B molecular description using the 3ob parameter set in gas phase), calculated for 22 drug-like small molecules. We compare the calculation of this value using free energy perturbation, Bennett’s acceptance ratio, Jarzynski’s equation, and Crooks’ equation. We also predict the reliability of each calculated Δ A g a s M M → 3 o b by evaluating several convergence criteria including sample size hysteresis, overlap statistics, and bias metric ( Π ). Within the total dataset, three distinct categories of molecules emerge: the “good” molecules, for which we can obtain converged Δ A g a s M M → 3 o b using Jarzynski’s equation; “bad” molecules which require Crooks’ equation to obtain a converged Δ A g a s M M → 3 o b ; and “ugly” molecules for which we cannot obtain reliably converged Δ A g a s M M → 3 o b with either Jarzynski’s or Crooks’ equations. We discuss, in depth, results from several example molecules in each of these categories and describe how dihedral discrepancies between levels of theory cause convergence failures even for these gas phase free energy simulations.


BIBECHANA ◽  
2018 ◽  
Vol 16 ◽  
pp. 92-105
Author(s):  
Sunil Pokharel ◽  
Shyam Prakash Khanal ◽  
N P Adhikari

Computer simulations of molecular models are powerful technique that have improved the under- standing of many biochemical phenomena. The method is frequently applied to study the motions of biological macromolecules such as protein and nucleic acids, which can be useful for interpreting the results of certain biophysical experiments. In this work, we have estimated the solvation free energy for light alkane (methane, ethane, propane and n-butane) dissolved in water and methanol respectively over a broad range of temperatures, from 275 K to 375 K, using molecular dynamics simulations. The alkane (methane, ethane, propane and n-butane), and methanol molecules are described by the OPLS-AA (Optimized Potentials for Liquid Simulations-All Atom) potential, while water is modeled by TIP3P (Transferable Intermolecular Potential with 3-Points) model. We have used the free energy perturbation method (Bennett Acceptance Ratio (BAR) method) for the calculation of free energy of solvation. The estimated values of solvation free energy of alkane in the corresponding solvents agree well with the available experimental data.BIBECHANA 16 (2019) 91-104


1991 ◽  
Vol 266 (18) ◽  
pp. 11801-11809
Author(s):  
N. Mizushima ◽  
D. Spellmeyer ◽  
S. Hirono ◽  
D. Pearlman ◽  
P. Kollman

2021 ◽  
Vol 7 (3) ◽  
pp. 467-475 ◽  
Author(s):  
Chun-Hui Zhang ◽  
Elizabeth A. Stone ◽  
Maya Deshmukh ◽  
Joseph A. Ippolito ◽  
Mohammad M. Ghahremanpour ◽  
...  

2002 ◽  
Vol 45 (23) ◽  
pp. 4995-5004 ◽  
Author(s):  
Cristiano Ruch Werneck Guimarães ◽  
Ricardo Bicca de Alencastro

Sign in / Sign up

Export Citation Format

Share Document