Leakiness and size exclusion of paracellular channels in cultured epithelial cell monolayers-interlaboratory comparison

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.

Solubility vs Dissolution in Physiological Bicarbonate Buffer

BACKGROUND

Phosphate buffer is often used as a replacement for the physiological bicarbonate buffer in pharmaceutical dissolution testing, although there are some discrepancies in their properties making it complicated to extrapolate dissolution results in phosphate to the in vivo situation. This study aims to characterize these discrepancies regarding solubility and dissolution behavior of ionizable compounds.

METHODS

The dissolution of an ibuprofen powder with a known particle size distribution was simulated in silico and verified experimentally in vitro at two different doses and in two different buffers (5 mM pH 6.8 bicarbonate and phosphate).

RESULTS

The results showed that there is a solubility vs. dissolution mismatch in the two buffers. This was accurately predicted by the in-house simulations based on the reversible non-equilibrium (RNE) and the Mooney models.

CONCLUSIONS

The results can be explained by the existence of a relatively large gap between the initial surface pH of the drug and the bulk pH at saturation in bicarbonate but not in phosphate, which is caused by not all the interfacial reactions reaching equilibrium in bicarbonate prior to bulk saturation. This means that slurry pH measurements, while providing surface pH estimates for buffers like phosphate, are poor indicators of surface pH in the intestinal bicarbonate buffer. In addition, it showcases the importance of accounting for the H2CO3-CO2 interconversion kinetics to achieve good predictions of intestinal drug dissolution.

A Refined Thin-Film Model for Drug Dissolution Considering Radial Diffusion - Simulating Powder Dissolution

PURPOSE

We aim to present a refined thin-film model describing the drug particle dissolution considering radial diffusion in spherical boundary layer, and to demonstrate the ability of the model to describe the dissolution behavior of bulk drug powders.

METHODS

The dissolution model introduced in this study was refined from a radial diffusion-based model previously published by our laboratory (So et al. in Pharm Res. 39:907-17, 2022). The refined model was created to simulate the dissolution of bulk powders, and to account for the evolution of particle size and diffusion layer thickness during dissolution. In vitro dissolution testing, using fractionated hydrochlorothiazide powders, was employed to assess the performance of the model.

RESULTS

Overall, there was a good agreement between the experimental dissolution data and the predicted dissolution profiles using the proposed model across all size fractions of hydrochlorothiazide. The model over-predicted the dissolution rate when the particles became smaller. Notably, the classic Nernst-Brunner formalism led to an under-estimation of the dissolution rate. Additionally, calculation based on the equivalent particle size derived from the specific surface area substantially over-predicted the dissolution rate.

CONCLUSION

The study demonstrated the potential of the radial diffusion-based model to describe dissolution of drug powders. In contrast, the classic Nernst-Brunner equation could under-estimate drug dissolution rate, largely due to the underlying assumption of translational diffusion. Moreover, the study indicated that not all surfaces on a drug particle contribute to dissolution. Therefore, relying on the experimentally-determined specific surface area for predicting drug dissolution is not advisable.

Predicting the intrinsic membrane permeability of Caco-2/MDCK cells 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.

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.

Revisiting the pKa-Flux method for determining intrinsic membrane permeability

Intrinsic membrane permeability is one of several factors that critically determine the intestinal absorption of a chemical. The intrinsic membrane permeability of a chemical is usually extracted from transwell experiments with Caco-2 or MDCK cells, preferably by the pKa-Flux method, which is considered the method of choice when aqueous boundary layer effects need to be excluded. The pKa-Flux method has two variants, the iso-pH method, where apical and basolateral pH are equal, and the gradient-pH method, where apical and basolateral pH are different. The most commonly used method is the gradient-pH method, as it is intended to reflect the pH-conditions in the gastrointestinal tract. However, concentration-shift effects caused by the applied pH-difference between apical and basolateral compartment in the gradient-pH method have not been considered in the evaluation of the experimental data in the past. Consequently, incorrect intrinsic membrane permeabilities have been determined. In this work, we present a revised method for extracting the intrinsic membrane permeability from gradient-pH data that considers concentration-shift effects in the basolateral aqueous boundary layer and filter as well as in the cytosol. Furthermore, we propose the use of the iso-pH method, where only concentration-shift effects in the cytosol need to be considered, as an alternative to the gradient-pH method. We use the five lipophilic bases amantadine, chloroquine, propranolol, venlafaxine and verapamil as examples to compare gradient-pH method and iso-pH method with regard to the extractability of the intrinsic membrane permeability. For lipophilic bases, the iso-pH method proves to be advantageous. All intrinsic membrane permeabilities determined in this work were substantially higher than the intrinsic membrane permeabilities reported in literature.

Predictions from First-Principles of Membrane Permeability to Small Molecules: How Useful Are They in Practice?

Predicting from first-principles the rate of passive permeation of small molecules across the biological membrane represents a promising strategy for screening lead compounds upstream in the drug-discovery and development pipeline. One popular avenue for the estimation of permeation rates rests on computer simulations in conjunction with the inhomogeneous solubility-diffusion model, which requires the determination of the free-energy change and position-dependent diffusivity of the substrate along the translocation pathway through the lipid bilayer. In this Perspective, we will clarify the physical meaning of the membrane permeability inferred from such computer simulations, and how theoretical predictions actually relate to what is commonly measured experimentally. We will also examine why these calculations remain both technically challenging and overly computationally expensive, which has hitherto precluded their routine use in nonacademic settings. We finally synopsize possible research directions to meet these challenges, increase the predictive power of physics-based rates of passive permeation, and, by ricochet, improve their practical usefulness.

Prediction of Liquid-Liquid Phase Separation at the Dissolving Drug Salt Particle Surface

During the dissolution of drug salt particles, liquid-liquid phase separation (LLPS) of a free form can occur within the unstirred water layer (UWL) of the particles (UWL-LLPS). Theoretically, UWL-LLPS occurs when the free form concentration at the salt particle surface (C₀) exceeds the intrinsic LLPS concentration (S₀^LLPS) of the free form. In the present study, we attempted to predict UWL-LLPS based on the intrinsic physicochemical properties of drugs. Cyproheptadine hydrochloride (CPH-HCl), diclofenac sodium (DCF-Na), papaverine hydrochloride (PAP-HCl), and propafenone hydrochloride (PRF-HCl) were selected as model drug salts. The pH₀ and C₀ values at pHs 4.0-9.5 (citric acid, phosphoric acid, and boric acid, buffer capacity = ca. 4 mM/ΔpH) were calculated using the pK_a, solubility product (K_sp), and diffusion coefficient (D) of a drug. S₀^LLPS was measured using the pH-shift method. UWL-LLPS was predicted to occur when C₀ ≥ S₀^LLPS. The prediction result was then compared with UWL-LLPS observed at each pH by polarized light microscopy (PLM). The pH-LLPS concentration (S_pH^LLPS) profile of each drug was also measured. UWL-LLPS was approximately correctly predicted for CPH-HCl, DCF-Na, and PRF-HCl. However, UWL-LLPS was not observable when C₀ was close to S₀^LLPS. Furthermore, UWL-LLPS was not accurately predicted in the case of PAP-HCl. The pH-S_pH^LLPS profile of PAP did not follow the Henderson-Hasselbalch equation, probably because of the formation of cationic aggregates. In conclusion, UWL-LLPS was approximately predictable for drug salts using their intrinsic physicochemical properties (K_sp, pK_a, D, and S₀^LLPS), except for PAP-HCl.

In silico model-based exploration of the applicability of parallel artificial membrane permeability assay (PAMPA) to screen chemicals of environmental concern

Parallel Artificial Membrane Permeability Assay (PAMPA) is an in vitro laboratory method for screening the transmembrane permeability of chemicals. Stemming from medicinal chemistry, PAMPA has the potential for use in the cost-effective high-throughput evaluation of chemicals of environmental concern. However, many chemicals of environmental concern differ substantially from pharmaceuticals in hydrophobicity and volatility. Here, we develop an in silico mass balance model to explore the impacts of chemical properties on chemical mass transfer in PAMPA and PAMPA's applicability to hydrophobic or volatile chemicals of environmental concern. The model's performance is evaluated by agreement between predicted and measured permeabilities of 1383 chemicals. The model predicts that the PAMPA measured permeability can be highly uncertain for hydrophobic chemicals because of considerable retention by the artificial membrane and for volatile chemicals because of substantial volatilization to the headspace. Notably, the permeabilities of hydrophobic chemicals are remarkably sensitive to specific experimental conditions, for example, the frequency of stirring and incubation time, challenging the comparison between measurements under different conditions. For hydrophobic chemicals, the PAMPA measured permeability may largely indicate the permeability of the unstirred water layer over the membrane, instead of the "intrinsic" permeability of the membrane, and therefore, may not be of interest for environmental exposure and risk assessments. The model also predicts that the time for mass transfer of highly hydrophobic chemicals to reach the steady state likely exceeds the incubation time, which violates the steady-state assumption used in calculating permeability from measured concentrations. Overall, our theoretical analysis underscores the importance to consider chemical properties when applying the current design of PAMPA to chemicals of environmental concern.

Impact of cholesterol and sphingomyelin on intrinsic membrane permeability

Transwell experiments with Caco-2 or MDCK cells are the gold standard for determining the intestinal permeability of chemicals. The intrinsic membrane permeability (P₀), that can be extracted from these experiments, might be comparable to P₀ measured in black lipid membrane (BLM) experiments and P₀ predicted by the solubility-diffusion model. Unfortunately, the overlap between experimental P_{0,Caco-2/MDCK} and P_0,BLM data is very small. So far, differences between both approaches have been attributed to the cholesterol and sphingomyelin content of cell membranes, but the database is too sparse to thoroughly test this theory. To create a diverse dataset, we measured P_0,BLM of ten chemicals in BLM experiments using DPhPC and DPhPC/cholesterol/sphingomyelin membranes. The results were compared to predicted BLM data and experimental Caco-2/MDCK data obtained from literature. While P_0,BLM of all chemicals was well predicted by the solubility-diffusion model, P_{0,Caco-2/MDCK} was only predictable for rather hydrophilic compounds with logarithmic hexadecane/water partition coefficients below -0.5. The effect of cholesterol and sphingomyelin on P_0,BLM was negligibly small.

Non-Effective Improvement of Absorption for Some Nanoparticle Formulations Explained by Permeability under Non-Sink Conditions

We evaluated the in vitro permeability of nanoparticle formulations of high and low lipophilic compounds under non-sink conditions, wherein compounds are not completely dissolved. The permeability of the highly lipophilic compound, griseofulvin, was improved by about 30% due to nanonization under non-sink conditions. Moreover, this permeability was about 50% higher than that under sink conditions. On the other hand, for the low lipophilic compound, hydrocortisone, there was no difference in permeability between micro-and nano-sized compounds under non-sink conditions. The nanonization of highly lipophilic compounds improves the permeability of the unstirred water layer (UWL), which in turn improves overall permeability. On the other hand, because the rate-limiting step in permeation for the low lipophilic compounds is the diffusion of the compounds in the membrane, the improvement of UWL permeability by nanonization does not improve the overall permeability. Based on this mechanism, nanoparticle formulations are not effective for low lipophilic compounds. To accurately predict the absorption of nanoparticle formulations, it is necessary to consider their permeability under non-sink conditions which reflect in vivo conditions.

Predicting Pharmacokinetics of Multisource Acyclovir Oral Products Through Physiologically Based Biopharmaceutics Modeling

Highly variable disposition after oral ingestion of acyclovir has been reported, although little is known regarding the underlying mechanisms. Different studies using the same reference product (Zovirax ®) showed that C_max and AUC were respectively 44 and 35% lower in Saudi Arabians than Europeans, consistent with higher frequencies of reduced-activity polymorphs of the organic cation transporter (OCT1) in Europeans. In this study, the contribution of physiology (i.e., OCT1 activity) to the oral disposition of acyclovir immediate release (IR) tablets was hypothesized to be greater than dissolution. The potential role of OCT1 was studied in a validated physiologically-based biopharmaceutics model (PBBM), while dissolution of two Chilean generics (with demonstrated bioequivalence) and the reference product was assessed in vitro. The PBBM suggested that OCT1 activity could partially explain population-related pharmacokinetic differences. Further, dissolution of generics was slower than the regulatory criterion for BCS III IR products. Remarkably, virtual bioequivalence (incorporating in vitro dissolution into the PBBM) correctly and robustly predicted the bioequivalence of these products, showcasing its value in support of failed BCS biowaivers. These findings suggest that very-rapid dissolution for acyclovir IR products may not be critical for BCS biowaiver. They also endorse the relevance of cross-over designs in bioequivalence trials.

Is equilibrium slurry pH a good surrogate for solid surface pH during drug dissolution?

The purpose of the present study was to investigate the suitability of equilibrium slurry pH (pH_eq) as a surrogate of solid surface pH during drug dissolution (pH₀). A comprehensive calculation scheme for pH_eq and pH₀ was formalized based on the principle of charge neutrality (equilibrium charge neutrality for pH_eq and charge flux neutrality for pH₀). The formalized scheme was then used to investigate the validity of pH₀ ≈ pH_eq approximation. The approximation of pH₀ ≈ pH_eq was suggested to be accurate for small molecules (ca. MW = 150) in high concentration buffer media (ca. buffer capacity = 30 mM/ΔpH). In addition, it is valid provided no precipitation of its free form for salts (vice versa for free forms) in both the slurry pH measurement and at the dissolving drug surface. The formalized calculation scheme is simple and applicable to free and salt form drugs in unbuffered and buffered media including bicarbonate buffer. The computational expense is very small so that it is applicable to various computer simulations such as biopharmaceutics modeling and simulation.

Global Analysis of Models for Predicting Human Absorption: QSAR, In Vitro, and Preclinical Models

Models intended to predict intestinal absorption are an essential part of the drug development process. Although many models exist for capturing intestinal absorption, many questions still exist around the applicability of these models to drug types like "beyond rule of 5" (bRo5) and low absorption compounds. This presents a challenge as current models have not been rigorously tested to understand intestinal absorption. Here, we assembled a large, structurally diverse dataset of ∼1000 compounds with known in vitro, preclinical, and human permeability and/or absorption data. In silico (quantitative structure-activity relationship), in vitro (Caco-2), and in vivo (rat) models were statistically evaluated for predictive performance against this human intestinal absorption dataset. We expect this evaluation to serve as a resource for DMPK scientists and medicinal/computational chemists to increase their understanding of permeability and absorption model utility and applications for academia and industry.

In Vitro Inhibition of Renal OCT2 and MATE1 Secretion by Antiemetic Drugs

The organic cation transporter 2 (OCT2) and multidrug and toxin extrusion protein 1 (MATE1) mediate the renal secretion of drugs. Recent studies suggest that ondansetron, a 5-HT₃ antagonist drug used to prevent nausea and vomiting, can inhibit OCT2- and MATE1-mediated transport. The purpose of this study was to test the ability of five 5-HT₃ antagonist drugs to inhibit the OCT2 and MATE1 transporters. The transport of the OCT2/MATE1 probe substrate ASP⁺ was assessed using two models: (1) HEK293 kidney cells overexpressing human OCT2 or MATE1, and (2) MDCK cells transfected with human OCT2 and MATE1. In HEK293 cells, the inhibition of ASP⁺ uptake by OCT2 listed in order of potency was palonosetron (IC₅₀: 2.6 μM) > ondansetron > granisetron > tropisetron > dolasetron (IC₅₀: 85.4 μM) and the inhibition of ASP⁺ uptake by MATE1 in order of potency was ondansetron (IC₅₀: 0.1 μM) > palonosetron = tropisetron > granisetron > dolasetron (IC₅₀: 27.4 μM). Ondansetron (0.5–20 μM) inhibited the basolateral-to-apical transcellular transport of ASP⁺ up to 64%. Higher concentrations (10 and 20 μM) of palonosetron, tropisetron, and dolasetron similarly reduced the transcellular transport of ASP⁺. In double-transfected OCT2-MATE1 MDCK cells, ondansetron at concentrations of 0.5 and 2.5 μM caused significant intracellular accumulation of ASP⁺. Taken together, these data suggest that 5-HT₃ antagonist drugs may inhibit the renal secretion of cationic drugs by interfering with OCT2 and/or MATE1 function.

Dose-Dependent Solubility-Permeability Interplay for Poorly Soluble Drugs under Non-Sink Conditions

We investigated the solubility-permeability interplay using a solubilizer additive under non-sink conditions. Sodium lauryl sulfate (SLS) was used as a solubilizer additive. The solubility and permeability of two poorly soluble drugs at various doses, with or without SLS, were evaluated by flux measurements. The total permeated amount of griseofulvin, which has high permeability, increased by the addition of SLS. On the other hand, triamcinolone, which has low permeability, showed an almost constant rate of permeation regardless of the SLS addition. The total permeated amount of griseofulvin increased by about 20-30% when the dose amount exceeded its solubility, whereas its concentration in the donor chamber remained almost constant. However, the total permeated amount of triamcinolone was almost constant regardless of dose amount. These results suggest that the permeability of the unstirred water layer (UWL) may be affected by SLS and solid drugs for high-permeable drugs. The effect of solid drugs could be explained by a reduction in the apparent UWL thickness. For the appropriate evaluation of absorption, it would be essential to consider these effects.

Toward Developing Discriminating Dissolution Methods for Formulations Containing Nanoparticulates in Solution: The Impact of Particle Drift and Drug Activity in Solution

Enabling formulations are an attractive approach to increase the dissolution rate, solubility, and oral bioavailability of poorly soluble compounds. With the growing prevalence of poorly soluble drug compounds in the pharmaceutical pipeline, supersaturating drug delivery systems (SDDS), a subset of enabling formulations, have grown in popularity due to their properties allowing for drug concentrations greater than the corresponding crystalline solubility. However, the extent of supersaturation generated as the enabling formulation traverses the gastrointestinal (GI) tract is dynamic and poorly understood. The dynamic nature of supersaturation is a result of several competing kinetic processes such as dissolution, solubilization by formulation and endogenous surfactants, crystallization, and absorption. Ultimately, the free drug concentration, which is equivalent to the drug's inherent thermodynamic activity amid these kinetic processes, defines the true driving force for drug absorption. However, in cases where solubilizing agents are present (i.e., surfactants and bile salts), drug molecules may associate with colloidal nanoscale species, complicating drug activity determination. These nanoscale species can drift into the aqueous boundary layer (ABL), increasing the local API activity at the membrane surface, resulting in increased bioavailability. Herein, a novel approach was developed to accurately measure thermodynamic drug activity in complex media containing drug distributed in nanoparticulate species. This approach captures the influence of the ABL on the observed flux and, ultimately, the predicted unbound drug concentration. The results demonstrate that this approach can help to (1) measure the true extent of local supersaturation in complex systems containing solubilizing excipients and (2) elucidate the mechanisms by which colloidal aggregates can modulate the drug activity in solution and potentially enhance the flux observed across a membrane. The utilization of these techniques may provide development scientists with a strategy to evaluate formulation sensitivity to nanospeciation and allow formulators to maximize the driving force for absorption in a complex environment, perhaps enabling the development of dissolution methods with greater discrimination and correlation to pre-clinical and clinical data sets.

Prediction of negative food effect induced by bile micelle binding on oral absorption of hydrophilic cationic drugs

The purpose of the present study was to quantitatively predict the negative food effect induced by bile micelle binding on the oral absorption of hydrophilic cationic drugs. The intrinsic membrane permeability and bile micelle unbound fraction of 12 model drugs (7 tertiary amines, 3 quaternary ammoniums, and 2 neutral drugs) were calculated from the experimental Caco-2 permeability data (P_app) under fasted and fed conditions. From these input data, the fraction of a dose absorbed (Fa) was predicted using the gastrointestinal unified theoretical framework, a mechanism-based oral absorption model. The predicted Fa ratio (fed/fasted) was then compared with the in vivo fed/fasted area under the plasma concentration-time curve ratio (AUCr). The AUCr values of tertiary amines and neutral drugs were appropriately predicted (absolute average fold error (AAFE) = 1.19), whereas those of quaternary ammoniums were markedly underestimated (AAFE = 4.70). The P_app ratio (fed/fasted) predicted AUCr less quantitatively (AAFE = 1.30 for tertiary amines and neutral drugs). The results of the present study would lead to a better understanding of negative food effect on oral drug absorption.

Prediction Characteristics of Oral Absorption Simulation Software Evaluated Using Structurally Diverse Low-Solubility Drugs

The purpose of the present study was to characterize current biopharmaceutics modeling and simulation software regarding the prediction of the fraction of a dose absorbed (Fa) in humans. As commercial software products, GastroPlus™ and Simcyp® were used. In addition, the gastrointestinal unified theoretical framework, a simple and publicly accessible model, was used as a benchmark. The Fa prediction characteristics for a total of 96 clinical Fa data of 27 model drugs were systematically evaluated using the default settings of each software product. The molecular weight, dissociation constant, octanol-water partition coefficient, solubility in biorelevant media, dose, and particle size of model drugs were used as input data. Although the same input parameters were used, GastroPlus™, Simcyp®, and the gastrointestinal unified theoretical framework showed different Fa prediction characteristics depending on the rate-limiting steps of oral drug absorption. The results of the present study would be of great help for the overall progression of physiologically based absorption models.

Structural dynamics of tight junctions modulate the properties of the epithelial barrier

Tight junctions are dynamic structures that are crucial in establishing the diffusion and electrical barrier of epithelial monolayers. Dysfunctions in the tight junctions can impede this barrier function and lead to many pathological conditions. Unfortunately, detailed understanding of the non-specific permeation pathway through the tight junctions, the so-called leak pathway, is lacking. We created computational models of the leak pathway to describe the two main barrier measures, molecular permeability and transepithelial electric resistance while using common structural dynamics. Our results showed that the proposed alternatives for the leak pathway, the bicellular strand opening dynamics and the tricellular pores, contribute together with distinct degrees, depending on the epithelium. The models can also capture changes in the tight junction barrier caused by changes in tight junction protein composition. In addition, we observed that the molecular permeability was markedly more sensitive to changes in the tight junction structure and strand dynamics compared with transepithelial electric resistance. The results highlight that our model creates a good methodological framework to integrate knowledge on the tight junction structure as well as to provide insights and tools to advance tight junction research.

Prediction Accuracy of Mechanism-Based Oral Absorption Model for Dogs

The purpose of the present study was to evaluate the prediction accuracy of a mechanism-based oral absorption model for the fraction of a dose absorbed (Fa) in dogs, focusing on poorly soluble drugs. As an open mechanism-based model, the gastrointestinal unified theoretical framework was used in this study. The prediction accuracy of the gastrointestinal unified theoretical framework was evaluated using Fa data in dogs (63 data sets for marketed drugs and proprietary compounds). For neutral compounds, Fa was accurately predicted, suggesting that the physiological parameters of dogs were appropriate except for gastrointestinal pH. An extensive literature survey on the small intestinal pH of dogs was then conducted. The result suggested that the pH value ranged between 6.5 and 7.5, with the midst value of 7.0, but there was a great variation among the literature. To confirm the appropriateness of this pH value, the Fa of free acid compounds was predicted by setting the small intestinal pH to 6.5, 7.0, and 7.5. The proportions of compounds with <2-fold error were 57%, 90%, and 76%, respectively. The results of the present study would enable an appropriate use of a mechanism-based model for drug discovery and development.

Passive Permeability of Planar Lipid Bilayers to Organic Anions

The membrane permeability P of organic ions was reported to be governed by the structure of the permeating molecule. Thus far, it is unclear whether the ion structure alters membrane partition or translocation proper across the membrane. Here, we obtained P values for 24 anionic compounds (18 concrete values, 6 upper limits) measuring the current that they carry through folded planar lipid bilayers. The P values range over more than 10 log units. Our measured permeability values correlate well (r = 0.95; logRMSE 0.74) with the hexadecane/water partition coefficients of the respective chemicals predicted by the COSMO-RS theory. Other attempts to predict P from the partition coefficient of the neutral molecule and from the solvation energy (Born energy) that opposes transfer into the membrane once the molecule is charged were unsuccessful. The uncertainties in assigning an effective radius to nonspherical molecules were much too large. The observation underlines that the actual structure of the molecules needs to be considered to predict partition and thus P by the solubility-diffusion model.

Physicochemical QSAR Analysis of Passive Permeability Across Caco-2 Monolayers

Caco-2 cell line is frequently used as a simplified in vitro model of intestinal absorption. In this study, a database of 1366 Caco-2 permeability coefficients (P_e) for 768 diverse drugs and drug-like compounds was compiled from public sources. The collected data represent permeation rates measured at varying experimental conditions (pH from 4.0 to 8.0, and stirring rates from 0 to >1000 rpm) that presumably account for passive diffusion across mucosal epithelium. These data were subjected to multistep nonlinear regression analysis using a minimal set of physicochemical descriptors (octanol-water log D, pKa, hydrogen bonding potential, and molecular size). The model was constructed in a mechanistic manner incorporating the following components: (i) a hydrodynamic equation of size- and charge-specific along with nonspecific diffusion across the paracellular pathway; (ii) transcellular diffusion represented by thermodynamic membrane/water partitioning ratio; (iii) stirring-dependent limit of maximum achievable permeability due to the presence of unstirred water layer. The obtained model demonstrates good accuracy of log P_e predictions with a residual mean square error <0.5 log units for all training and validation sets. Given its robust performance and straightforward interpretation in terms of simple physicochemical properties, the proposed model may serve as a valuable tool to guide drug discovery efforts toward readily absorbable compounds.

MDCKpred: a web-tool to calculate MDCK permeability coefficient of small molecule using membrane-interaction chemical features

Structure-based models to understand the transport of small molecules through biological membrane can be developed by enumerating intermolecular interactions of the small molecule with a biological membrane, usually a dimyristoylphosphatidylcholine (DMPC) monolayer. This ADME (absorption, distribution, metabolism, and excretion) property based on Madin-Darby Canine Kidney (MDCK) cell line demonstrates intestinal drug absorption of small molecules and correlated to human intestinal absorption which acts as a determining factor to forecast small-molecule prioritization in drug-discovery projects. We present here the development of MDCKpred web-tool which calculates MDCK permeability coefficient of small molecule based on the regression model, developed using membrane-interaction chemical features. The web-tool allows users to calculate the MDCK permeability coefficient (nm/s) instantly by providing simple descriptor inputs. The chemical-interaction features are derived from different parts of the DMPC molecule viz. head, middle, and tail regions and accounts overall intermolecular contacts of the small molecule when passively diffused through the phospholipid-rich biological membrane. The MDCKpred model is both internally (R² = .76; [Formula: see text]= .68; R_train = .87; R_test = .69) and externally (R_ext = .55) validated. Furthermore, we used natural molecules as application examples to demonstrate its utility in lead exploration and optimization projects. The MDCKpred web-tool can be accessed freely at http://www.mdckpred.in . This web-tool is designed to offer an intuitive way of prioritizing small molecules based on calculated MDCK permeabilities.

The impact of food intake on the luminal environment and performance of oral drug products with a view to in vitro and in silico simulations: a PEARRL review

Objectives

Using the type of meal and dosing conditions suggested by regulatory agencies as a basis, this review has two specific objectives: first, to summarize our understanding on the impact of food intake on luminal environment and drug product performance and second, to summarize the usefulness and limitations of available in vitro and in silico methodologies for the evaluation of drug product performance after food intake.

Key findings

Characterization of the luminal environment and studies evaluating product performance in the lumen, under conditions suggested by regulatory agencies for simulating the fed state, are limited. Various in vitro methodologies have been proposed for evaluating drug product performance in the fed state, but systematic validation is lacking. Physiologically based pharmacokinetic (PBPK) modelling approaches require the use of in vitro biorelevant data and, to date, have been used primarily for investigating the mechanisms via which an already observed food effect is mediated.

Summary

Better understanding of the impact of changes induced by the meal administration conditions suggested by regulatory agencies on the luminal fate of the drug product is needed. Relevant information will be useful for optimizing the in vitro test methods and increasing the usefulness of PBPK modelling methodologies.

Ranking Itraconazole Formulations Based on the Flux through Artificial Lipophilic Membrane

PURPOSE

The goal of the study was to evaluate a miniaturized dissolution-permeation apparatus (μFLUX™ apparatus) for its ability to benchmark several itraconazole (ITZ) formulations for which in vivo PK data was available in the literature.

METHOD

Untreated and micronized powders of ITZ and various enabling formulations of ITZ (commercial Sporanox® solid dispersion, a Soluplus®-based solid dispersion and a nanosuspension) were introduced to the donor compartment of μFLUX™ apparatus. Donor and acceptor chambers were divided from each other by a lipophilic membrane. In addition to the flux evaluations, changes in solid state as a function of time were investigated to gain further insight into the flux changes observed over time for the solid dispersion formulations.

RESULTS

Initial flux values from Sporanox®, the nanosuspension and the micronized ITZ showed ratios of 52/4/1 with a decreasing flux from nanosuspension and both solid dispersions after 2.5-3 h. Although the initial flux from the Soluplus® formulation was 2.2 times lower than the one observed for Sporanox®, the decrease in flux observed was milder and became ~ 2 times higher than Sporanox® after approximately 2.5 h. The total amounts of ITZ in the receiver compartment after 240 min showed the same rank order as the rodent AUCs of these formulations reported in literature.

CONCLUSIONS

It was demonstrated that in vitro flux measurements using lipophilic artificial membranes could correctly reproduce the rank order of PK results for ITZ formulations. The drop in flux over time for solid dispersions could be backed by experimental indications of crystallization.

Closed-Loop Doluisio (Colon, Small Intestine) and Single-Pass Intestinal Perfusion (Colon, Jejunum) in Rat-Biophysical Model and Predictions Based on Caco-2

PURPOSE

The effective rat intestinal permeability (P _eff ) was deconvolved using a biophysical model based on parameterized paracellular, aqueous boundary layer, transcellular permeabilities, and the villus-fold surface area expansion factor.

METHODS

Four types of rat intestinal perfusion data were considered: single-pass intestinal perfusion (SPIP) in the jejunum (n = 40), and colon (n = 15), closed-loop (Doluisio type) in the small intestine (n = 78), and colon (n = 74). Moreover, in vitro Caco-2 permeability values were used to predict rat in vivo values in the rat data studied.

RESULTS

Comparable number of molecules permeate via paracellular water channels as by the lipoidal transcellular route in the SPIP method, although in the closed-loop method, the paracellular route appears dominant in the colon. The aqueous boundary layer thickness in the small intestine is comparable to that found in unstirred in vitro monolayer assays; it is thinner in the colon. The mucosal surface area in anaesthetized rats is 0.96-1.4 times the smooth cylinder calculated value in the colon, and it is 3.1-3.6 times in the small intestine. The paracellular permeability of the intestine appeared to be greater in rat than human, with the colon showing more leakiness (higher P _para ) than the small intestine.

CONCLUSION

Based on log intrinsic permeability values, the correlations between the in vitro and in vivo models ranged from r² 0.82 to 0.92. The SPIP-Doluisio method comparison indicated identical log permeability selectivity trend with negligible bias.

Impact of capillary flow hydrodynamics on carrier-mediated transport of opioid derivatives at the blood-brain barrier, based on pH-dependent Michaelis-Menten and Crone-Renkin analyses

Most studies of blood-brain barrier (BBB) permeability and transport are conducted at a single pH, but more detailed information can be revealed by using multiple pH values. A pH-dependent biophysical model was applied to the mechanistic analysis of published pH-dependent BBB luminal uptake data from three opioid derivatives in rat: pentazocine (Suzuki et al., 2002a, 2002b), naloxone (Suzuki et al., 2010a), and oxycodone (Okura et al., 2008). Two types of data were processed: in situ brain perfusion (ISBP) and brain uptake index (BUI). The published perfusion data were converted to apparent luminal permeability values, P_app, and analyzed by the pCEL-X program (Yusof et al., 2014), using the pH-dependent Crone-Renkin equation (pH-CRE) to determine the impact of cerebrovascular flow on the Michaelis-Menten transport parameters (Avdeef and Sun, 2011). For oxycodone, the ISBP data had been measured at pH7.4 and 8.4. The present analysis indicates a 7-fold lower value of the cerebrovascular flow velocity, F_pf, than that expected in the original study. From the pyrilamine-inhibited data, the flow-corrected passive intrinsic permeability value was determined to be P₀=398×10^-6cm·s^-1. The uptake data indicate that the neutral form of oxycodone is affected by a transporter at pH8.4. The extent of the cation uptake was less certain from the available data. For pentazocine, the brain uptake by the BUI method had been measured at pH5.5, 6.5, and 7.4, in a concentration range 0.1-40mM. Under similar conditions, ISBP data were also available. The pH-CRE determined values of F_pf from both methods were nearly the same, and were smaller than the expected value in the original publication. The transport of the cationic pentazocine was not fully saturated at pH5.5 at 40mM. The transport of the neutral species at pH7.4 appeared to reach saturation at 40mM pentazocine concentration, but not at 12mM. In the case of naloxone, a pH-dependent Michaelis-Menten equation (pH-MME) analysis of the data indicated a smooth sigmoidal transition from a higher capacity uptake process affecting cationic naloxone (pH5.0-7.0) to a lower capacity uptake process affecting the neutral drug (pH8.0-8.5), with cross-over point near pH7.4. Evidently, measurements at multiple pH values can reveal important information about both cerebrovascular flow and BBB transport kinetics.

Application of a Human Intestinal Epithelial Cell Monolayer to the Prediction of Oral Drug Absorption in Humans as a Superior Alternative to the Caco-2 Cell Monolayer

A human small intestinal epithelial cell (HIEC) monolayer was recently established in our laboratories as a novel system to evaluate the Papp (apparent permeability coefficient) of compounds during their absorption in humans. An effusion-based analysis using polyethylene glycol oligomers with molecular weights ranging from 194-898 indicated that HIEC and Caco-2 cell monolayers both had paracellular pores with 2 distinct radiuses (∼ 5 and 9-14 Å), whereas the porosity of large pores was 11-fold higher in the HIEC monolayer (44 × 10(-8)) than in the Caco-2 cells (4 × 10(-8)). A comparison between the fraction-absorbed (Fa) values observed in humans and those predicted from Papp values in both monolayers indicated that the HIEC monolayer had markedly higher precision to predict Fa values with root mean square error of 9.40 than the Caco-2 cells (root mean square error = 16.90) for 10 paracellularly absorbed compounds. Furthermore, the accuracy of the HIEC monolayer to classify the absorption of 23 test drugs with diverse absorption properties, including different pathways in the presence or absence of susceptibility to efflux transporters, was higher than that of the Caco-2 cell monolayer. In conclusion, the HIEC monolayer exhibited advantages over Caco-2 cells in the ranking and prediction of absorption of compounds in humans.

Transepithelial Transport of PAMAM Dendrimers Across Isolated Human Intestinal Tissue

Poly(amido amine) (PAMAM) dendrimers have shown transepithelial transport across intestinal epithelial barrier in rats and across Caco-2 cell monolayers. Caco-2 models innately lack mucous barriers, and rat isolated intestinal tissue has been shown to overestimate human permeability. This study is the first report of transport of PAMAM dendrimers across isolated human intestinal epithelium. It was observed that FITC labeled G4-NH2 and G3.5-COOH PAMAM dendrimers at 1 mM concentration do not have a statistically higher permeability compared to free FITC controls in isolated human jejunum and colonic tissues. Mannitol permeability was increased at 10 mM concentrations of G3.5-COOH and G4-NH2 dendrimers. Significant histological changes in human colonic and jejunal tissues were observed at G3.5-COOH and G4-NH2 concentrations of 10 mM implying that dose limiting toxicity may occur at similar concentrations in vivo. The permeability through human isolated intestinal tissue in this study was compared to previous rat and Caco-2 permeability data. This study implicates that PAMAM dendrimer oral drug delivery may be feasible, but it may be limited to highly potent drugs.

Atenolol Renal Secretion Is Mediated by Human Organic Cation Transporter 2 and Multidrug and Toxin Extrusion Proteins

Atenolol is a β-blocker widely used in the treatment of hypertension. Atenolol is cleared predominantly by the kidney by both glomerular filtration and active secretion, but the molecular mechanisms involved in its renal secretion are unclear. Using a panel of human embryonic kidney cell lines stably expressing human organic cation transporter (hOCT) 1-3, human organic anion transporter (hOAT) 1, hOAT3, human multidrug and toxin extrusion protein (hMATE) 1, and hMATE2-K, we found that atenolol interacted with both organic cation and anion transporters. However, it is transported by hOCT1, hOCT2, hMATE1, and hMATE2-K, but not by hOCT3, hOAT1, and hOAT3. A detailed kinetic analysis coupled with absolute quantification of membrane transporter proteins by liquid chromatography-tandem mass spectrometry revealed that atenolol is an excellent substrate for the renal transporters hOCT2, hMATE1, and hMATE2-K. The Km values for hOCT2, hMATE1, and hMATE2-K are 280 ± 4, 32 ± 5, and 76 ± 14 μM, respectively, and the calculated turnover numbers are 2.76, 0.41, and 2.20 s(-1), respectively. To demonstrate unidirectional transepithelial transport of atenolol, we developed and functionally validated a hOCT2/hMATE1 double-transfected Madin-Darby canine kidney cell culture model. Transwell studies showed that atenolol transport in the basal (B)-to-apical (A) direction is 27-fold higher than in the A-to-B direction, whereas its B-to-A/A-to-B transport ratio was only 2 in the vector-transfected control cells. The overall permeability of atenolol in the B-to-A direction in hOCT2/hMATE1 cells was 44-fold higher than in control cells. Together, our data support that atenolol tubular secretion is mediated through the hOCT2/hMATEs secretion pathway and suggest a significant role of organic cation transporters in the disposition of an important antihypertensive drug.

In vitro porcine blood-brain barrier model for permeability studies: pCEL-X software pKa(FLUX) method for aqueous boundary layer correction and detailed data analysis

In vitro blood-brain barrier (BBB) models from primary brain endothelial cells can closely resemble the in vivo BBB, offering valuable models to assay BBB functions and to screen potential central nervous system drugs. We have recently developed an in vitro BBB model using primary porcine brain endothelial cells. The model shows expression of tight junction proteins and high transendothelial electrical resistance, evidence for a restrictive paracellular pathway. Validation studies using small drug-like compounds demonstrated functional uptake and efflux transporters, showing the suitability of the model to assay drug permeability. However, one limitation of in vitro model permeability measurement is the presence of the aqueous boundary layer (ABL) resulting from inefficient stirring during the permeability assay. The ABL can be a rate-limiting step in permeation, particularly for lipophilic compounds, causing underestimation of the permeability. If the ABL effect is ignored, the permeability measured in vitro will not reflect the permeability in vivo. To address the issue, we explored the combination of in vitro permeability measurement using our porcine model with the pKa(FLUX) method in pCEL-X software to correct for the ABL effect and allow a detailed analysis of in vitro (transendothelial) permeability data, Papp. Published Papp using porcine models generated by our group and other groups are also analyzed. From the Papp, intrinsic transcellular permeability (P0) is derived by simultaneous refinement using a weighted nonlinear regression, taking into account permeability through the ABL, paracellular permeability and filter restrictions on permeation. The in vitro P0 derived for 22 compounds (35 measurements) showed good correlation with P0 derived from in situ brain perfusion data (r(2)=0.61). The analysis also gave evidence for carrier-mediated uptake of naloxone, propranolol and vinblastine. The combination of the in vitro porcine model and the software analysis provides a useful tool to better predict BBB permeability in vivo and gain better mechanistic information about BBB permeation.

Critical analysis of 3-D organoid in vitro cell culture models for high-throughput drug candidate toxicity assessments

Drug failure due to toxicity indicators remains among the primary reasons for staggering drug attrition rates during clinical studies and post-marketing surveillance. Broader validation and use of next-generation 3-D improved cell culture models are expected to improve predictive power and effectiveness of drug toxicological predictions. However, after decades of promising research significant gaps remain in our collective ability to extract quality human toxicity information from in vitro data using 3-D cell and tissue models. Issues, challenges and future directions for the field to improve drug assay predictive power and reliability of 3-D models are reviewed.

Prediction of passive drug permeability across the blood-retinal barrier

PURPOSE

The purpose of this study is to develop a computational model of the physical barrier function of the outer blood-retinal barrier (BRB), which is vital for normal retinal function. To our best knowledge no comprehensive models of BRB has been reported.

METHODS

The model construction is based on the three-layered structure of the BRB: retinal pigment epithelium (RPE), Bruch's membrane and choriocapillaris endothelium. Their permeabilities were calculated based on the physical theories and experimental material and permeability studies in the literature, which were used to describe diffusional hindrance in specific environments.

RESULTS

Our compartmental BRB model predicts permeabilities with magnitudes similar to the experimental values in the literature. However, due to the small number and varying experimental conditions there is a large variability in the available experimental data, rendering validation of the model difficult. The model suggests that the paracellular pathway of the RPE largely defines the total BRB permeability.

CONCLUSIONS

Our model is the first BRB model of its level and combines the present knowledge of the BRB barrier function. Furthermore, the model forms a platform for the future model development to be used for the design of new drugs and drug administration systems.

Evaluation of the effects of absorption enhancers on Caco-2 cell monolayers by using a pore permeation model involving two different sizes

We applied a parallel pore permeation model based on the Renkin molecular sieving function by using two different-sized pathways to analyze the permeation-enhancing effects of poly-L-arginine (PLA) or a mixed system of spermine (SPM) and sodium taurocholate (STC). Four paracellular markers were simultaneously applied to Caco-2 cell monolayers, and a set of apparent permeability coefficient (P) values was used to obtain membrane parameters. For PLA treatment, the pore occupancy/length ratio (ε/L) of the large pathways increased while the pore radius (R) did not, suggesting that the number of large pathways for the relatively large hydrophilic molecules in the monolayers could be increased by the addition of PLA. In contrast, application of the mixed system comprising SPM and STC significantly increased not only the R of the large pathways but also ε/L of the small pathways. Such changes in membrane parameters could be related to the enhancing mechanism of these compounds. The simulation curves for molecular weight (MW)-P calculated from the membrane parameters could be used to predict the P of drugs with different MWs.

Physiologically based pharmacokinetic modelling of drug penetration across the blood-brain barrier--towards a mechanistic IVIVE-based approach

Predicting the penetration of drugs across the human blood-brain barrier (BBB) is a significant challenge during their development. A variety of in vitro systems representing the BBB have been described, but the optimal use of these data in terms of extrapolation to human unbound brain concentration profiles remains to be fully exploited. Physiologically based pharmacokinetic (PBPK) modelling of drug disposition in the central nervous system (CNS) currently consists of fitting preclinical in vivo data to compartmental models in order to estimate the permeability and efflux of drugs across the BBB. The increasingly popular approach of using in vitro-in vivo extrapolation (IVIVE) to generate PBPK model input parameters could provide a more mechanistic basis for the interspecies translation of preclinical models of the CNS. However, a major hurdle exists in verifying these predictions with observed data, since human brain concentrations can't be directly measured. Therefore a combination of IVIVE-based and empirical modelling approaches based on preclinical data are currently required. In this review, we summarise the existing PBPK models of the CNS in the literature, and we evaluate the current opportunities and limitations of potential IVIVE strategies for PBPK modelling of BBB penetration.

Dynamic contrast enhanced optical imaging of cervix, in vivo: a paradigm for mapping neoplasia-related parameters

We present a novel biophotonic method and imaging modality for estimating and mapping neoplasia-specific functional and structural parameters of the cervical precancerous epithelium. Estimations were based on experimental data obtained from dynamic contrast-enhanced optical imaging of cervix, in vivo. We have developed a pharmacokinetic, in silico, model of the optical tracer's uptake by the epithelium. We have identified that the kinetic parameters of the model correlate well with pathologic alterations in both metabolic and structural characteristics of the tissue, associated with the neoplasia progress. Global sensitivity analysis and global optimization methods were employed for identifying the key determinant set of biological parameters that dictate the model's output. Particularly, the shuffled complex evolution algorithm converged to a set of four parameters that can be estimated with an error of 7%, indicating a good accuracy and precision. These results are unique in the sense that for the first time functional and microstructural parameter maps can be estimated and displayed together, thus maximizing the diagnostic information. The quantity and the quality of this information are unattainable by other invasive and non invasive methods.

In situ artificial membrane permeation assay under hydrodynamic control: correlation between drug in vitro permeability and fraction absorbed in humans

The purpose of this study was to develop an in vitro permeation model that will predict the fraction of drugs absorbed in humans. A rotating-diffusion cell with two aqueous compartments, separated by a lipid-impregnated artificial membrane, was used to determine the permeability of drugs under conditions of controlled hydrodynamics. The measured effective permeability coefficient was modified to include the paracellular transport derived from a previously reported colorectal adenocarcinoma epithelial cell line (Caco-2) permeability study and the effects of unstirred water layer anticipated in vivo. Permeability data were collected for 31 different marketed drugs with known absolute oral bioavailability and human hepatic clearance data. Literature bioavailability values were corrected for the first pass hepatic clearance thus obtaining the fraction absorbed from intestinal lumen (fraction absorbed), F(a), while assuming that the fraction escaping intestinal extraction, F(g), was approximately ~1. Permeability obtained under conditions of controlled hydrodynamics was compared with the permeability measured under unstirred conditions. It is shown that the optimized effective permeability correlates with the fraction absorbed. In contrast, permeability data obtained under unstirred conditions does not show a good correlation. The in vitro permeation model developed in this study predicts the fraction absorbed of the selected drugs in humans within experimental uncertainty. It has been demonstrated that the correlation with the fraction absorbed is greatly improved using the permeability data obtained under controlled hydrodynamics with paracellular transport included in the model.