In silico prediction of drug solubility: 1. Free energy of hydration

As a first step in the computational prediction of drug solubility the free energy of hydration, DeltaG*(vw) in TIP4P water has been computed for a data set of 48 drug molecules using the free energy of perturbation method and the optimized potential for liquid simulations all-atom force field. The simulations were performed in two steps, where first the Coulomb and then the Lennard-Jones interactions between the solute and the water molecules were scaled down from full to zero strength to provide physical understanding and simpler predictive models. The results have been interpreted using a theory assuming DeltaG*(vw) = A(MS)gamma + E(LJ) + E(C)/2 where A(MS) is the molecular surface area, gamma is the water-vapor surface tension, and E(LJ) and E(C) are the solute-water Lennard-Jones and Coulomb interaction energies, respectively. It was found that by a proper definition of the molecular surface area our results as well as several results from the literature were found to be in quantitative agreement using the macroscopic surface tension of TIP4P water. This is in contrast to the surface tension for water around a spherical cavity that previously has been shown to be dependent on the size of the cavity up to a radius of approximately 1 nm. The step of scaling down the electrostatic interaction can be represented by linear response theory.