Extracting Microscopic Energies from Oil-Phase Solvation Experiments

Abstract

There is a large literature that reports oil−water partition coefficients of small molecule compounds. A goal of much of this work is to obtain hydrophobicity parameters for models of folding, docking, binding, conformational changes, partitioning, and flow, mainly of biomolecules and drugs. The quantity obtained from experiments is the difference between a solute's oil- and water-phase chemical potentials Δμ ⋮ μwat − μoil, but the quantity needed for models is a contact free energy per unit area. Usually, solvation modelers find this quantity in two steps: (i) subtract off some assumed “irrelevant” solute entropies from the experimental Δμ's and then (ii) divide the result by the solute area. Here we show that such a procedure can lead to large errors. We use a two -dimensional lattice model of hydrocarbon liquids in NVT Monte Carlo simulations, where the underlying contact energy is rigorously known in advance, to test strategies for converting partitioning data to contact energies. We find that rather...