Zwitterions contain both positively and negatively charged functional groups, resulting in an overall net neutral charge. Nevertheless, the membrane permeability of the zwitterionic form of a compound is assumed to be much lower than the permeability of the uncharged neutral form. Although a significant proportion of pharmaceuticals are zwitterionic, it has not been clear so far whether their permeability is dominated by the permeation of the zwitterionic or the neutral form, since neutral fractions are often quite low as compared to the zwitterionic fraction. This complicates the in silico prediction of the permeability of zwitterionic compounds. In this work, we re-evaluated existing in vitro permeability data from literature measured with Caco-2/MDCK cell assays, using more strict exclusion criteria for effects like diffusion limitation by the aqueous boundary layers, paracellular transport, active transport and retention. Using this re-evaluated data set, we show that extracted intrinsic permeabilities of the neutral fraction are well predicted by the solubility-diffusion model (RMSE = 1.21; n = 18) if the permeability of the zwitterionic species is assumed negligible. Our work thus suggests that only the neutral species is relevant for the membrane permeability of zwitterionic compounds, and that membrane permeability of zwitterionic compounds is indeed predictable by the solubility-diffusion model.

Membrane permeability is one of the main determinants for the absorption, distribution, metabolism and excretion of compounds and is therefore of crucial importance for successful drug development. Experiments with artificial phospholipid membranes have shown that the intrinsic membrane permeability (P0) of compounds is well-predicted by the solubility-diffusion model (SDM). However, using the solubility-diffusion model to predict the P0 of biological Caco-2 and MDCK cell membranes has proven unreliable so far. Recent publications revealed that many published P0 extracted from Caco-2 and MDCK experiments are incorrect. In this work, we therefore used a small self-generated set as well as a large revised set of experimental Caco-2 and MDCK data from literature to compare experimental and predicted P0. The P0 extracted from Caco-2 and MDCK experiments were systematically lower than the P0 predicted by the solubility-diffusion model. However, using the following correlation: log P0,Caco-2/MDCK = 0.84 log P0,SDM - 1.85, P0 of biological Caco-2 and MDCK cell membranes was well-predicted by the solubility-diffusion model.

This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.