scholarly journals Thermodynamic Decomposition of Solvation Free Energies with Particle Mesh Ewald and Long-Range Lennard-Jones Interactions in Grid Inhomogeneous Solvation Theory

Author(s):  
Lieyang Chen ◽  
Anthony Cruz ◽  
Daniel R. Roe ◽  
Andrew C. Simmonett ◽  
Lauren Wickstrom ◽  
...  
2020 ◽  
Author(s):  
Lieyang Chen ◽  
Anthony Cruz ◽  
Daniel R. Roe ◽  
Andrew Simmonett ◽  
Lauren Wickstrom ◽  
...  

Grid Inhomogeneous Solvation Theory (GIST) maps out solvation thermodynamic properties on a fine meshed grid and provides a statistical mechanical formalism for thermodynamic end-state calculations. However, differences in how long-range non-bonded interactions are calculated in molecular dynamics engines and in the current implementation of GIST have prevented precise comparisons between free energies estimated using GIST and those from other free energy methods such as thermodynamic integration (TI). Here, we address this by presenting PME-GIST, a formalism by which particle mesh Ewald (PME) based electrostatic energies and long-range Lennard-Jones (LJ) energies are decomposed and assigned to individual atoms and the corresponding voxels they occupy in a manner consistent with the GIST approach. PME-GIST yields potential energy calculations that are precisely consistent with modern simulation engines and performs these calculations at a dramatically faster speed than prior implementations. Here, we apply PME-GIST end-states analyses to 32 small molecules whose solvation free energies are close to evenly distributed from 2 kcal/mol to -17 kcal/mol and obtain solvation energies consistent with TI calculations (R2 = 0.99, mean unsigned difference 0.8 kcal/mol). We also estimate the entropy contribution from the 2nd and higher order entropy terms that are truncated in GIST by the differences between entropies calculated in TI and GIST. With a simple correction for the high order entropy terms, PME-GIST obtains solvation free energies that are highly consistent with TI calculations (R2 = 0.99, mean unsigned difference = 0.4 kcal/mol) and experimental results (R2 = 0.88, mean unsigned difference = 1.4 kcal/mol). The precision of PME-GIST also enables us to show that the solvation free energy of small hydrophobic and hydrophilic molecules can be largely understood based on perturbations of the solvent in a region extending a few solvation shells from the solute. We have integrated PME-GIST into the open-source molecular dynamics analysis software CPPTRAJ.


2020 ◽  
Author(s):  
Lieyang Chen ◽  
Anthony Cruz ◽  
Daniel R. Roe ◽  
Andrew Simmonett ◽  
Lauren Wickstrom ◽  
...  

Grid Inhomogeneous Solvation Theory (GIST) maps out solvation thermodynamic properties on a fine meshed grid and provides a statistical mechanical formalism for thermodynamic end-state calculations. However, differences in how long-range non-bonded interactions are calculated in molecular dynamics engines and in the current implementation of GIST have prevented precise comparisons between free energies estimated using GIST and those from other free energy methods such as thermodynamic integration (TI). Here, we address this by presenting PME-GIST, a formalism by which particle mesh Ewald (PME) based electrostatic energies and long-range Lennard-Jones (LJ) energies are decomposed and assigned to individual atoms and the corresponding voxels they occupy in a manner consistent with the GIST approach. PME-GIST yields potential energy calculations that are precisely consistent with modern simulation engines and performs these calculations at a dramatically faster speed than prior implementations. Here, we apply PME-GIST end-states analyses to 32 small molecules whose solvation free energies are close to evenly distributed from 2 kcal/mol to -17 kcal/mol and obtain solvation energies consistent with TI calculations (R2 = 0.99, mean unsigned difference 0.8 kcal/mol). We also estimate the entropy contribution from the 2nd and higher order entropy terms that are truncated in GIST by the differences between entropies calculated in TI and GIST. With a simple correction for the high order entropy terms, PME-GIST obtains solvation free energies that are highly consistent with TI calculations (R2 = 0.99, mean unsigned difference = 0.4 kcal/mol) and experimental results (R2 = 0.88, mean unsigned difference = 1.4 kcal/mol). The precision of PME-GIST also enables us to show that the solvation free energy of small hydrophobic and hydrophilic molecules can be largely understood based on perturbations of the solvent in a region extending a few solvation shells from the solute. We have integrated PME-GIST into the open-source molecular dynamics analysis software CPPTRAJ.


2020 ◽  
Author(s):  
Lieyang Chen ◽  
Anthony Cruz ◽  
Daniel R. Roe ◽  
Andrew Simmonett ◽  
Lauren Wickstrom ◽  
...  

Grid Inhomogeneous Solvation Theory (GIST) maps out solvation thermodynamic properties on a fine meshed grid and provides a statistical mechanical formalism for thermodynamic end-state calculations. However, differences in how long-range non-bonded interactions are calculated in molecular dynamics engines and in the current implementation of GIST have prevented precise comparisons between free energies estimated using GIST and those from other free energy methods such as thermodynamic integration (TI). Here, we address this by presenting PME-GIST, a formalism by which particle mesh Ewald (PME) based electrostatic energies and long-range Lennard-Jones (LJ) energies are decomposed and assigned to individual atoms and the corresponding voxels they occupy in a manner consistent with the GIST approach. PME-GIST yields potential energy calculations that are precisely consistent with modern simulation engines and performs these calculations at a dramatically faster speed than prior implementations. Here, we apply PME-GIST end-states analyses to 32 small molecules whose solvation free energies are close to evenly distributed from 2 kcal/mol to -17 kcal/mol and obtain solvation energies consistent with TI calculations (R2 = 0.99, mean unsigned difference 0.8 kcal/mol). We also estimate the entropy contribution from the 2nd and higher order entropy terms that are truncated in GIST by the differences between entropies calculated in TI and GIST. With a simple correction for the high order entropy terms, PME-GIST obtains solvation free energies that are highly consistent with TI calculations (R2 = 0.99, mean unsigned difference = 0.4 kcal/mol) and experimental results (R2 = 0.88, mean unsigned difference = 1.4 kcal/mol). The precision of PME-GIST also enables us to show that the solvation free energy of small hydrophobic and hydrophilic molecules can be largely understood based on perturbations of the solvent in a region extending a few solvation shells from the solute. We have integrated PME-GIST into the open-source molecular dynamics analysis software CPPTRAJ.


2005 ◽  
Vol 70 (11) ◽  
pp. 1769-1786 ◽  
Author(s):  
Luc A. Vannier ◽  
Chunxiang Yao ◽  
František Tureček

A computational study at correlated levels of theory is reported to address the structures and energetics of transient radicals produced by hydrogen atom abstraction from C-1, C-2, C-3, C-4, C-5, O-1, O-3, and O-5 positions in 2-deoxyribofuranose in the gas phase and in aqueous solution. In general, the carbon-centered radicals are found to be thermodynamically and kinetically more stable than the oxygen-centered ones. The most stable gas-phase radical, 2-deoxyribofuranos-5-yl (5), is produced by H-atom abstraction from C-5 and stabilized by an intramolecular hydrogen bond between the O-5 hydroxy group and O-1. The order of radical stabilities is altered in aqueous solution due to different solvation free energies. These prefer conformers that lack intramolecular hydrogen bonds and expose O-H bonds to the solvent. Carbon-centered deoxyribose radicals can undergo competitive dissociations by loss of H atoms, OH radical, or by ring cleavages that all require threshold dissociation or transition state energies >100 kJ mol-1. This points to largely non-specific dissociations of 2-deoxyribose radicals when produced by exothermic hydrogen atom abstraction from the saccharide molecule. Oxygen-centered 2-deoxyribose radicals show only marginal thermodynamic and kinetic stability and are expected to readily fragment upon formation.


2014 ◽  
Vol 28 (8) ◽  
pp. 825-829 ◽  
Author(s):  
Sreeja Parameswaran ◽  
David L. Mobley

2006 ◽  
Vol 419 (1-3) ◽  
pp. 240-244 ◽  
Author(s):  
Takumi Hori ◽  
Hideaki Takahashi ◽  
Masayoshi Nakano ◽  
Tomoshige Nitta ◽  
Weitao Yang

2014 ◽  
Vol 28 (3) ◽  
pp. 135-150 ◽  
Author(s):  
David L. Mobley ◽  
Karisa L. Wymer ◽  
Nathan M. Lim ◽  
J. Peter Guthrie

Sign in / Sign up

Export Citation Format

Share Document