Estimation of the concentration boundary layer adjacent to a flat surface using computational fluid dynamics

Dissolution-permeation (D/P) experiments are widely used during preclinical development due to producing results with better predictability than traditional monophasic experiments. However, it is difficult to compare absorption across in vitro setups given the propensity to only report apparent permeability. We therefore developed an approach to predict the concentration boundary layer for any D/P device by using computational fluid dynamics (CFD). The Navier-Stokes and continuity equation in 2D were solved numerically in MATLAB and by finite element methods in COMSOL v6.1 to predict the momentum [Formula: see text] and concentration ηg boundary layer for a flow over a flat plate, i.e. the classical Blasius boundary layer flow. A MATLAB algorithm was developed to calculate the edge of either boundary layer. The methodology to determine the concentration boundary layer based on Blasius's analysis provided an accurate estimate for both [Formula: see text] and ηg, resulting in, [Formula: see text] , at high Schmidt numbers (Sc ∼ 1000) within 14 % of the Blasius solution and 6.6 % of the accepted Schmidt number correlation ( [Formula: see text] ). The methodology based on the Blasius analysis of the concentration boundary layer using velocity and concentration profiles computed using CFD presented herein will enable characterization/analysis of complex D/P apparatuses used in preclinical development, where an analytical solution may not be available.

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.

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.

Desorption kinetics of hydrophobic organic contaminants from marine plastic pellets

This study investigated the desorption behavior of polychlorinated biphenyls (PCBs) from marine plastic pellets. Long-term desorption experiments were conducted using field-collected polyethylene (PE) pellets. The results indicate that the desorption kinetics highly depends on the PE-water partition coefficients of PCB congeners. After 128 d of the experiment, the smallest congener considered (CB 8) had desorbed nearly completely (98%), whereas major fractions (90-99%) of highly chlorinated congeners remained in the pellets. An intraparticle diffusion model mostly failed to reproduce the desorption kinetics, whereas an aqueous boundary layer (ABL) diffusion model well approximated the data. The desorption half-lives are estimated to 14 d to 210 years for CB 8 to CB 209 in an actively stirred solution (ABL thickness: 30 μm). Addition of methanol to water enhanced the desorption to a large extent. A need for further work to explore roles of organic matter in facilitating solute transfer is suggested.