Can membrane permeability of zwitterionic compounds be predicted by the solubility-diffusion model?

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.

Pitfalls in evaluating permeability experiments with Caco-2/MDCK cell monolayers

When studying the transport of molecules across biological membranes, intrinsic membrane permeability (P0) is more informative than apparent permeability (Papp), because it eliminates external (setup-specific) factors, provides consistency across experiments and mechanistic insight. It is thus an important building block for modeling the total permeability in any given scenario. However, extracting P0 is often difficult, if not impossible, when the membrane is not the dominant transport resistance. In this work, we set out to analyze Papp values measured with Caco-2/MDCK cell monolayers of 69 literature references. We checked the Papp values for a total of 318 different compounds for the extractability of P0, considering possible limitations by aqueous boundary layers, paracellular transport, recovery issues, active transport, a possible proton flux limitation, and sink conditions. Overall, we were able to extract 77 reliable P0 values, which corresponds to about one quarter of the total compounds analyzed, while about half were limited by the diffusion through the aqueous layers. Compared to an existing data set of P0 values published by Avdeef, our approach resulted in a much higher exclusion of compounds. This is a consequence of stricter compound- and reference-specific exclusion criteria, but also because we considered possible concentration-shift effects due to different pH values in the aqueous layers, an effect only recently described in literature. We thus provide a consistent and reliable set of P0, e.g. as a basis for future modeling.

Human gut epithelium features recapitulated in MINERVA 2.0 millifluidic organ-on-a-chip device

We developed an innovative millifluidic organ-on-a-chip device, named MINERVA 2.0, that is optically accessible and suitable to serial connection. In the present work, we evaluated MINERVA 2.0 as millifluidic gut epithelium-on-a-chip by using computational modeling and biological assessment. We also tested MINERVA 2.0 in a serially connected configuration prodromal to address the complexity of multiorgan interaction. Once cultured under perfusion in our device, human gut immortalized Caco-2 epithelial cells were able to survive at least up to 7 days and form a three-dimensional layer with detectable tight junctions (occludin and zonulin-1 positive). Functional layer development was supported by measurable trans-epithelial resistance and FITC-dextran permeability regulation, together with mucin-2 expression. The dynamic culturing led to a specific transcriptomic profile, assessed by RNASeq, with a total of 524 dysregulated transcripts (191 upregulated and 333 downregulated) between static and dynamic condition. Overall, the collected results suggest that our gut-on-a-chip millifluidic model displays key gut epithelium features and, thanks to its modular design, may be the basis to build a customizable multiorgan-on-a-chip platform.

Predicting the Permeability of Macrocycles from Conformational Sampling - Limitations of Molecular Flexibility

Macrocycles constitute superior ligands for targets that have flat binding sites but often require long synthetic routes, emphasizing the need for property prediction prior to synthesis. We have investigated the scope and limitations of machine learning classification models and of regression models for predicting the cell permeability of a set of denovo-designed, drug-like macrocycles. 2D-Based classification models, which are fast to calculate, discriminated between macrocycles that had low-medium and high permeability and may be used as virtual filters in early drug discovery projects. Importantly, stereo- and regioisomer were correctly classified. QSPR studies of two small sets of comparator drugs suggested that use of 3D descriptors, calculated from biologically relevant conformations, would allow development of more precise regression models for late phase drug projects. However, a 3D permeability model could only be developed for a rigid series of macrocycles. Comparison of NMR based conformational analysis with in silico conformational sampling indicated that this shortcoming originates from the inability of the molecular mechanics force field to identify the relevant conformations for flexible macrocycles. We speculate that a Kier flexibility index of ≤10 constitutes a current upper limit for reasonably accurate 3D prediction of macrocycle cell permeability.

Permeability prediction for zwitterions via chromatographic indexes and classification into 'certain' and 'uncertain'

Background: The development of zwitterions to a drug is likely to be more challenging than compounds of other charge types. Results: Two chromatographic indexes (log k'80 PLRP-S and log K_W^IAM) can be successfully used as permeability classifiers of ampholytes. Moreover, a pragmatic classification into ordinary ampholytes; zwitterions 'certain' (i.e., the zwitterionic species is dominant in the physiological pH range); and zwitterions 'uncertain' (multiple species are present in the physiological pH range) enables to study the permeability of ampholytic compounds in relation to species distribution. Methodology: Potentiometry (pK_a), reversed-phase (RP)-chromatography, tri-layer parallel artificial membrane permeability assays, quantitative structure-property relationships (QSPR) and block relevance (BR) analysis, receiver operating characteristic (ROC) curves. Conclusion: Structures considered as poorly permeable like zwitterions can be integrated in drug discovery programs by applying ad hoc experimental and computational tools.

Rationalization of the Membrane Permeability Differences in a Series of Analogue Cyclic Decapeptides

Cyclization and selected backbone N-methylations are found to be often necessary but not sufficient conditions for peptidic drugs to have a good bioavailability. Thus, the design of cyclic peptides with good passive membrane permeability and good solubility remains a challenge. The backbone scaffold of a recently published series of cyclic decapeptides with six selected backbone N-methylations was designed to favor the adoption of a closed conformation with β-turns and four transannular hydrogen bonds. Although this conformation was indeed adopted by the peptides as determined by NMR measurements, substantial differences in the membrane permeability were observed. In this work, we aim to rationalize the impact of discrete side chain modifications on membrane permeability for six of these cyclic decapeptides. The thermodynamic and kinetic properties were investigated using molecular dynamics simulations and Markov state modeling in water and chloroform. The study highlights the influence that side-chain modifications can have on the backbone conformation. Peptides with a d-proline in the β-turns were more likely to adopt, even in water, the closed conformation with transannular hydrogen bonds, which facilitates transition through the membrane. The population of the closed conformation in water was found to correlate positively with PAMPA log P_e.

A Fast Chromatographic Method for Estimating Lipophilicity and Ionization in Nonpolar Membrane-Like Environment

This study describes the design and implementation of a new chromatographic descriptor called log k'80 PLRP-S that provides information about the lipophilicity of drug molecules in the nonpolar environment, both in their neutral and ionized form. The log k'80 PLRP-S obtained on a polymeric column with acetonitrile/water mobile phase is shown to closely relate to log Ptoluene (toluene dielectric constant ε ∼ 2). The main intermolecular interactions governing log k'80 PLRP-S were deconvoluted using the Block Relevance (BR) analysis. The information provided by this descriptor was compared to ElogD and calclog Ptol, and the differences are highlighted. The "charge-flush" concept is introduced to describe the sensitivity of log k'80 PLRP-S to the ionization state of compounds in the pH range 2 to 12. The ability of log k'80 PLRP-S to indicate the propensity of neutral molecules and monoanions to form Intramolecular Hydrogen Bonds (IMHBs) is proven through a number of examples.

Heteroaryl chalcone allied triazole conjugated organosilatranes: synthesis, spectral analysis, antimicrobial screening, photophysical and theoretical investigations

A series of heteroaryl tethered triazole conjoined organosilatranes were synthesized and studied for their solvatochromism experimentally and theoretically followed by antimicrobial screening.