+ 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

Footprinting molecular electrostatic potential surfaces for calculation of solvation energies

Footprinting molecular electrostatic potential surfaces for calculation of solvation energies

Physical Chemistry Chemical Physics 15(41): 18262-18273

A liquid is composed of an ensemble of molecules that populate a large number of different states, so calculation of the solvation energy of a molecule in solution requires a method for summing the interactions with the environment over all of these states. The surface site interaction model for the properties of liquids at equilibrium (SSIMPLE) simplifies the surface of a molecule to a discrete number of specific interaction sites (SSIPs). The thermodynamic properties of these interaction sites can be characterised experimentally, for example, through measurement of association constants for the formation of simple complexes that feature a single H-bonding interaction. Correlation of experimentally determined solution phase H-bond parameters with gas phase ab initio calculations of maxima and minima on molecular electrostatic potential surfaces (MEPS) provides a method for converting gas phase calculations on isolated molecules to parameters that can be used to estimate solution phase interaction free energies. This approach has been generalised using a footprinting technique that converts an MEPS into a discrete set of SSIPs (each described by a polar interaction parameter, εi). These SSIPs represent the molecular recognition properties of the entire surface of the molecule. For example, water is described by four SSIPs, two H-bond donor sites and two H-bond acceptor sites. A liquid mixture is described as an ensemble of SSIPs that represent the components of the mixture at appropriate concentrations. Individual SSIPs are assumed to be independent, so speciation of SSIP contacts can be calculated based on properties of the individual SSIP interactions, which are given by the sum of a polar (εiεj) and a non-polar (E(vdW)) interaction term. Results are presented for calculation the free energies of transfer of a range of organic molecules from the pure liquid into water, from the pure liquid into n-hexadecane, from n-hexadecane into water, from n-octanol into water, and for the transfer of water from pure water into a range of organic liquids. The agreement with experiment is accurate to within 1.6-3.9 kJ mol(-1) root mean square difference, which suggests that the SSIMPLE approach is a promising method for estimation of solvation energies in more complex systems.

Please choose payment method:

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

Accession: 053279278

Download citation: RISBibTeXText

PMID: 24064723

DOI: 10.1039/c3cp53158a

Related references

Calculation of the total electrostatic energy of a macromolecular system: solvation energies, binding energies, and conformational analysis. Proteins 4(1): 7, 1988

Computation of methodology-independent ionic solvation free energies from molecular simulations. I. The electrostatic potential in molecular liquids. Journal of Chemical Physics 124(12): 124106, 2006

Performance comparison of generalized born and Poisson methods in the calculation of electrostatic solvation energies for protein structures. Journal of Computational Chemistry 25(2): 265-284, 2004

Calculation of solvation interaction energies for protein adsorption on polymer surfaces. Journal of Biomaterials Science. Polymer Edition 3(2): 127-147, 1991

Bounding the electrostatic free energies associated with linear continuum models of molecular solvation. Journal of Chemical Physics 130(10): 104108, 2009

Prediction of Aqueous Solvation Free Energies from Properties of Solute Molecular Surface Electrostatic Potentials. Journal of Physical Chemistry A 103(12): 1853-1856, 1999

Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions. Journal of Chemical Physics 147(16): 161716, 2017

Solvation model based on order parameters and a fast sampling method for the calculation of the solvation free energies of peptides. Journal of Physical Chemistry. B 110(3): 1476-1484, 2006

Solvation free energies for periodic surfaces: comparison of implicit and explicit solvation models. Physical Chemistry Chemical Physics 18(46): 31850-31861, 2016

Accurate first principles calculation of molecular charge distributions and solvation energies from Ab initio quantum mechanics and continuum dielectric theory. Journal of the American Chemical Society 116(26): 11875-11882, 1994

Free Energies of Solvation with Surface, Volume, and Local Electrostatic Effects and Atomic Surface Tensions to Represent the First Solvation Shell. Journal of Chemical Theory and Computation 6(4): 1109-1117, 2010

Molecular Potential Energies from Experimental Electric Field and Electrostatic Potential at Nuclei. Journal of Physical Chemistry. a 123(10): 2018-2024, 2019

Progress in ab initio QM/MM free-energy simulations of electrostatic energies in proteins: accelerated QM/MM studies of pKa, redox reactions and solvation free energies. Journal of Physical Chemistry. B 113(5): 1253-1272, 2009

Quantitative analysis of molecular surfaces: areas, volumes, electrostatic potentials and average local ionization energies. Journal of Molecular Modeling 16(11): 1679-1691, 2010

Calculation of electrostatic interaction energies in molecular dimers from atomic multipole moments obtained by different methods of electron density partitioning. Journal of Computational Chemistry 25(7): 921-934, 2004