Prediction of Phospholipid-Water Partition Coefficients of Ionic Organic Chemicals Using the Mechanistic Model COSMOmic.

The partition coefficient of chemicals from water to phospholipid membrane, K(lipw), is of central importance for various fields. For neutral organic molecules, log K(lipw) correlates with the log of bulk solvent-water partition coefficients such as the octanol-water partition coefficient. However, this is not the case for charged compounds, for which a mechanistic modeling approach is highly necessary. In this work, we extend the model COSMOmic, which adapts the COSMO-RS theory for anisotropic phases and has been shown to reliably predict K(lipw) for neutral compounds, to the use of ionic compounds. To make the COSMOmic model applicable for ionic solutes, we implemented the internal membrane dipole potential in COSMOmic. We empirically optimized the potential with experimental K(lipw) data of 161 neutral and 75 ionic compounds, yielding potential shapes that agree well with experimentally determined potentials from the literature. This model refinement has no negative effect on the prediction accuracy of neutral compounds (root-mean-square error, RMSE = 0.62 log units), while it highly improves the prediction of ions (RMSE = 0.70 log units). The refined COSMOmic is, to our knowledge, the first mechanistic model that predicts K(lipw) of both ionic and neutral species with accuracies better than 1 log unit.