ABSTRACT
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.
Subject(s)
Models, Chemical , Organic Chemicals/chemistry , Phospholipids/chemistry , Water/chemistry , Ions , Membrane PotentialsABSTRACT
Results obtained with two computational approaches for the simulation of ion motion at elevated pressure are compared with experimentally derived ion current data. The computational approaches used are charged particle tracings with the software package SIMION ver. 8 and finite element based calculations using the software package Comsol Multiphysics ver. 4.0/4.0a. The experimental setup consisted of a tubular corona discharge ion source coupled to a cylindrical measurement chamber held at atmospheric pressure. Generated ions are flown into the chamber at essentially subsonic laminar isothermal conditions. In the simulations, strictly stationary conditions were assumed. The results show very good agreement between the SIMION/SDS model and experimental data. For the Comsol model, only qualitative agreement is observed.