Description of hydration free energy density as a function of molecular physical properties

Abstract

A method to calculate the solvation free energy density (SFED) at any point in the cavity surface or solvent volume surrounding a solute is proposed. In the special case in which the solvent is water, the SFED is referred to as the hydration free energy density (HFED). The HFED is described as a function of some physical properties of the molecules. These properties are represented by simple basis functions. The hydration free energy of a solute was obtained by integrating the HFED over the solvent volume surrounding the solute, using a grid model. Of 34 basis functions that were introduced to describe the HFED, only six contribute significantly to the HFED. These functions are representations of the surface area and volume of the solute, of the polarization and dispersion of the solute, and of two types of electrostatic interactions between the solute and its environment. The HFED is described as a linear combination of these basis functions, evaluated by summing the interaction energy between each atom of the solute with a grid point in the solvent, where each grid point is a representation of a finite volume of the solvent. The linear combination coefficients were determined by minimizing the error between the calculated and experimental hydration free energies of 81 neutral organic molecules that have a variety of functional groups. The calculated hydration free energies agree well with the experimental results. The hydration free energy of any other solute molecule can then be calculated by summing the product of the linear combination coefficients and the basis functions for the solute.