Modeling the pharmacodynamics of passive membrane permeability

Small molecule permeability through cellular membranes is critical to a better understanding of pharmacodynamics and the drug discovery endeavor. Such permeability may be estimated as a function of the free energy change of barrier crossing by invoking the barrier domain model, which posits that permeation is limited by passage through a single “barrier domain” and assumes diffusivity differences among compounds of similar structure are negligible. Inspired by the work of Rezai and co-workers (JACS 128:14073–14080, 2006), we estimate this free energy change as the difference in implicit solvation free energies in chloroform and water, but extend their model to include solute conformational affects. Using a set of eleven structurally diverse FDA approved compounds and a set of thirteen congeneric molecules, we show that the solvation free energies are dominated by the global minima, which allows solute conformational distributions to be effectively neglected. For the set of tested compounds, the best correlation with experiment is obtained when the implicit chloroform global minimum is used to evaluate the solvation free energy difference. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s10822-011-9480-7) contains supplementary material, which is available to authorized users.