Next Article in Journal
Experimental Evidence for Diiodohydroxyquinoline-Induced Neurotoxicity: Characterization of Age and Gender as Predisposing Factors
Previous Article in Journal
Utilization of Polymeric Micelles as a Lucrative Platform for Efficient Brain Deposition of Olanzapine as an Antischizophrenic Drug via Intranasal Delivery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Vitro–In Silico Modeling of Caffeine and Diclofenac Permeation in Static and Fluidic Systems with a 16HBE Lung Cell Barrier

1
Department of Clinical Pharmacy, Saarland University, 66123 Saarbrucken, Germany
2
Fraunhofer Institute for Biomedical Engineering IBMT, 66280 Sulzbach, Germany
3
Department of Internal Medicine V, Saarland University, 66421 Homburg, Germany
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2022, 15(2), 250; https://doi.org/10.3390/ph15020250
Submission received: 17 January 2022 / Revised: 11 February 2022 / Accepted: 16 February 2022 / Published: 18 February 2022
(This article belongs to the Section Pharmacology)

Abstract

:
Static in vitro permeation experiments are commonly used to gain insights into the permeation properties of drug substances but exhibit limitations due to missing physiologic cell stimuli. Thus, fluidic systems integrating stimuli, such as physicochemical fluxes, have been developed. However, as fluidic in vitro studies display higher complexity compared to static systems, analysis of experimental readouts is challenging. Here, the integration of in silico tools holds the potential to evaluate fluidic experiments and to investigate specific simulation scenarios. This study aimed to develop in silico models that describe and predict the permeation and disposition of two model substances in a static and fluidic in vitro system. For this, in vitro permeation studies with a 16HBE cellular barrier under both static and fluidic conditions were performed over 72 h. In silico models were implemented and employed to describe and predict concentration–time profiles of caffeine and diclofenac in various experimental setups. For both substances, in silico modeling identified reduced apparent permeabilities in the fluidic compared to the static cellular setting. The developed in vitro–in silico modeling framework can be expanded further, integrating additional cell tissues in the fluidic system, and can be employed in future studies to model pharmacokinetic and pharmacodynamic drug behavior.

1. Introduction

With the introduction of the European Directive 2010/63/EU on the protection of animals used for scientific purposes [1], the 3R principles, which aim to replace, reduce, and refine preclinical animal studies [2], received legal recognition for the first time. For this, further investigations into how in vitro and in silico approaches can support overcoming the limitations of current preclinical methods whilst reducing the need for animal studies are required [3].
In vitro permeation experiments can be used to characterize the transport of drug molecules across biological barriers, such as the pulmonary [4,5,6] or intestinal epithelia [7,8,9] barrier. The information gained in permeation studies on the properties of drug substances or environmental pollutants and their interaction with permeation barriers can be employed to estimate transport through physiologic endothelial and epithelial barriers in vivo [7,10,11,12]. A common approach for investigating the respective apparent permeability through a cell monolayer in vitro is the application of a single-barrier model [9,11,13]. Here, the apparent permeability is a composite parameter that not only comprises diffusion through the apical and basolateral membrane but that also considers additional processes such as diffusion through unstirred aqueous layers or through porous filter membranes and that has been successfully used for the prediction of in vivo absorption processes [7,10,11].
However, while static in vitro cell culture experiments are still widely used to investigate research questions, static monoculture approaches show several limitations regarding missing cell stimuli that are present in vivo [14,15]. One major group of missing stimuli is that of physicochemical fluxes [14,15]. To integrate this stimulus and to generate more predictive in vitro models, the development and application of fluidic and microfluidic systems have been proposed as early stage approaches for the investigation of drug compounds and complex biological processes [14,16,17,18,19,20].
Fluidic systems show distinct advantages over their static counterparts, as they hold the potential to mimic physiologic conditions more closely, replicating the physiological flows of tissue fluids, introducing shear stress, enabling the spatiotemporal control of medium supply, and ensuring waste removal, among other aspects [17,21,22]. Additionally, by integrating multiple tissues in multi-organ fluidic systems, tissue interactions and physiological processes, such as the bioactivation of compounds, can be simulated [21,22,23,24,25]. However, the interpretation and translation of the corresponding in vitro readouts are challenging [19].
Here, in silico modeling techniques that are already widely applied in drug discovery and development [26,27,28] can be helpful tools to support study design, data analysis, and extrapolation of in vitro studies [19,29,30]. The integration of in vitro and in silico modeling approaches offers the opportunity to quantitatively evaluate the results of fluidic in vitro systems and to use simulation techniques to predict yet untested scenarios [7,29,31]. Moreover, in silico models hold the potential to integrate prior knowledge into simulation experiments to help inform decisions related to study design (e.g., time of sampling) and to gain information on a drug’s pharmacokinetic and pharmacodynamic behavior, allowing extrapolation of drug doses and decreasing the need for in vivo studies [7,25,29,31].
In this study, the modular MultiCompartmental Bioreactor (MCmB) array Quasi-Vivo® was employed [16,20,32,33,34], and an in silico model mimicking the in vitro fluidic experimental system was implemented. The in silico model was used to simulate concentration–time profiles of the investigated probe substances caffeine and diclofenac. For model predictions in the fluidic system, the apparent permeability of the substances was estimated based on static in vitro experiments that were conducted to obtain prior knowledge regarding the experimental setup and the model substances.
Thus, the presented work aimed to develop in silico models to describe and predict the permeation and disposition of the two model substances, caffeine and diclofenac, in a static and fluidic in vitro system. For this, a compartmental modeling approach was employed to simulate drug concentration–time profiles over a period of 72 h after single-dose administration. The in vitro experiments were performed both for a control scenario (transfer through a porous membrane without cells) and for permeation through a 16HBE epithelial cell monolayer cultured on a porous membrane. Subsequently, the experimental setups were assessed with the implemented in silico models.

2. Results

In vitro permeation studies with 16HBE cells (cell-based experiments) and without cells (control) were performed under static and fluidic conditions. At eight different time points (15 min–72 h), the apical and basolateral concentrations of the two probe substances (caffeine and diclofenac) were quantified via high-performance liquid chromatography (HPLC). Based on these data, in silico models mimicking the in vitro setup were developed to stimulate concentration–time profiles of the model substances.
Before performing the permeation studies, the 16HBE cell model was established to control the integrity of the cell layer and the viability of the cells. The integrity of the 16HBE cell layer was assessed using transepithelial electrical resistance (TEER) measurements, showing TEER values of 671.46 Ω*cm2 ± 46.21 Ω*cm2. The viability of the lung cells was verified by an alamarBlue™ assay. Compared to the cell viability in the static system, which was defined as 100% with a standard deviation of 15.2% (n = 3), the viability under fluidic conditions was 118.9% ± 7.1% (n = 3), while the difference was not statistically significant (p < 0.05).

2.1. Static Permeation Studies

For the in silico modeling of caffeine and diclofenac permeation in fluidic systems, in vitro studies under static cultivation conditions were performed to obtain prior information on model substance behavior in the experimental setup. In the static permeation experiments, the apical concentrations of both caffeine and diclofenac substantially decreased within the first two hours (first three measurements). Concurrently, the basolateral concentrations increased and reached equilibrium with the apical concentrations after ~22 h (Figure 1 and Figure 2).
For the static permeation experiments without lung cells (control), model parameter estimation resulted in apparent permeabilities (P12) of 3.27 × 10−3 cm/min for caffeine and 2.51× 10−3 cm/min for diclofenac (Table 1). With these input parameters, the observed rapid decline in apical concentrations and the observed rapid incline in basolateral concentrations were well described by the static in silico model (Figure 1). Since no cellular barrier was involved in the static control experiments, an equilibrium was reached between the apical and basolateral compartments that yielded similar concentrations in both chambers at the end of the experiment. Of note, several quantified basolateral diclofenac concentrations (t ≥ 22 h) were slightly elevated compared to the corresponding model-simulated concentrations (Figure 1b). Further, the model simulations suggest that the observed equilibrium might have already been reached after ~10 h for caffeine and after ~14 h for diclofenac (Figure 1).
To model the static in vitro cell permeation studies, the experimental data supported the assumption of equal permeabilities from the apical to the basolateral compartment and vice versa since an equilibrium with similar apical and basolateral concentrations of ~40 µM was reached for both of the modeled compounds after ~24 h (non-polarized permeation, Figure 2). The application of the estimated permeabilities (Table 1) of 2.62 × 10−3 cm/min for caffeine and 1.38 × 10−3 cm/min for diclofenac led to a close concordance between the experimentally quantified (in vitro) and model-simulated (in silico) apical and basolateral concentrations. Similar to the static control permeation experiment, several observed basolateral diclofenac concentrations (t ≥ 22 h) were slightly elevated compared to the model-simulated concentrations (Figure 2b). Additionally, the model simulations suggest that for caffeine, equilibrium might have already been reached after ~12 h (Figure 2a).
The estimated apparent permeabilities and residual variabilities for the static permeation experiments (both control and cell experiments) are shown in Table 1.
Of note, the proportional residual error term in the in silico model for the control experiments was estimated to be equal to 0 and was hence fixed to this value to ensure proper calculation of the covariance matrix and the resulting residual errors.

2.2. Fluidic Permeation Studies

The in silico model for the fluidic system was applied to predict caffeine and diclofenac concentrations in the fluidic system by using the estimated apparent permeabilities derived from the static permeation experiments. All other model input parameters (i.e., compartment volumes and flow rates) were fixed to experimentally determined values. Hence, the depicted model simulations in Figure 3 (control) and Figure 4a,b (cell-based experiments) represent pure predictions of apical and basolateral concentrations based on the static in vitro–in silico modeling platform.
The resulting model simulations for the fluidic study without a lung cell barrier led to precise predictions of the measured caffeine concentrations in both the apical and basolateral circuits (Figure 3a). While the model predictions of the basolateral diclofenac concentrations were in close agreement with the observed data, the measured apical diclofenac concentrations barely decreased during the experiment and were elevated compared to model-simulated concentrations (Figure 3b).
Due to repeated sampling without solvent replacement (eight times) at both the apical and basolateral three-way cock (hereafter referred to as the sampling and dosing port, respectively), the total volume of the apical and basolateral circuits decreased from 15 mL to 11 mL. This resulting reduction in volume (and amount of compound) led to a step-like decrease in simulated apical concentrations at the sampling time points that could also be detected in the observed concentrations (Figure 3).
Using the estimated apparent permeability parameters from the static cell experiments, substance transfer from the apical to the basolateral circuit was overpredicted in the fluidic system for both caffeine and diclofenac (Figure 4a,b). After adapting the apparent permeabilities from the static system (factor (F) of 0.34 for caffeine and 0.17 for diclofenac, respectively), the observed concentration–time courses in the fluidic system could be described well with the in silico model (Figure 4c,d). The resulting permeabilities (P16; apparent permeability between the apical and basolateral chambers in the fluidic system) were 8.92 × 10−4 cm/min for caffeine and 2.35 × 10−4 cm/min for diclofenac.
Model simulations for the fluidic system, which are depicted in Figure 5, show the impact of repeated sampling in the in silico model on apical and basolateral concentrations. While considering the decreasing volume and, consequently, decreasing amount of drug substance led to precise predictions of the measured caffeine in vitro concentrations (Figure 5a, grey profiles), disregarding the sampling would have resulted in an overprediction of the apical concentrations. In contrast, the effect on basolateral concentrations was negligible. Of note, for diclofenac, both implemented approaches led to an underprediction of the apical concentrations (Figure 5b).
In the fluidic experimental setup, the apical concentrations measured at time t = 15 min were lower compared to the concentrations measured at t = 75 min and t = 120 min (Figure 6). This was in agreement with model simulations where apical concentrations at the sampling port increased during the first ~45 min until equilibrium was reached in the apical circuit (depicted in Figure 6, zooming into the first 2.5 h of the fluidic concentration-time profiles).

3. Discussion

The presented work aimed to develop in silico models to describe and predict the permeation and disposition of model substances in a static and fluidic in vitro system. For this, a compartmental modeling approach was used to simulate drug concentration–time profiles over a period of 72 h. In vitro experiments were performed with 16HBE epithelial cells, which were qualified for the permeation studies based on characteristic TEER values [4,35,36] and the non-significant difference in cell viability between static and fluidic cultivation conditions.
Based on the implemented in vitro–in silico modeling framework, concentration–time profiles of the two model substances in both static and fluidic conditions could be described and predicted. Caffeine, a rather hydrophilic (logP~−0.2) and uncharged molecule [37,38], and diclofenac, a highly lipophilic (logP~4) and mainly charged molecule [38,39], were selected as model substances to evaluate the applicability of the in vitro–in silico system.
The experimental data that were collected in the static cellular study suggested that neither caffeine nor diclofenac were subject to active transport regarding the human bronchial epithelial cell line 16HBE since an equilibrium with similar apical and basolateral concentrations was reached at the end of the in vitro experiments. Hence, the permeabilities from the apical to the basolateral compartment and vice versa were assumed to be equal. The bidirectional single-barrier model was suitable to describe the concentration–time profiles in the static system for both the control and the cell-based studies and successfully led to the estimation of apparent permeabilities regarding the transfer of both model substances.
Using the estimated apparent permeabilities from the static in vitro experiments, the fluidic control experiments could be well predicted for the apical and basolateral caffeine concentrations as well as for the basolateral diclofenac concentrations.
In contrast, when employing the estimated apparent permeabilities from the static in vitro study for the fluidic 16HBE cell experiments, caffeine and diclofenac transfer to the basolateral circuit was overpredicted. This could be due to the fact that cellular behavior can differ in fluidic compared to static systems, as shown for cell proliferation, biotransformation, and/or up-regulation of proteins and metabolic enzymes, among others [21,22,40,41,42]. As a result, transfer through the lung cell layer might have been altered in the fluidic compared to the static system due to changes in, for example, protein expression and cell morphology. After empirically adapting the permeabilities from the static in vitro studies, the in silico model could also precisely describe the experimental results from the lung cell barrier study. Here, the estimated adaption factors were 0.34 for caffeine and 0.17 for diclofenac, implying a highly reduced apparent permeability compared to the static in vitro experiments. This case study suggests that estimated apparent permeabilities from static in vitro systems could require refinement or adaptation before application for predictions in fluidic systems.
The estimated apparent permeabilities of 2.62 × 10−3 cm/min (static system) and 8.92 × 10−4 cm/min (fluidic system) for caffeine in the cell-based experiments are within the same order of magnitude with a published value from a study by Mathias et al. (1.31 × 10−3 cm/min) that used the human bronchial epithelial cell line Calu-3 [5]. The apparent permeabilities for diclofenac were estimated to be 1.38 × 10−3 (static system) and 2.35 × 10−4 cm/min (fluidic system) in the cell-based experiments, and, consequently, were about 2–3.5 times lower compared to the permeabilities for caffeine.
Fluidic permeation studies were performed using the MCmB QuasiVivo® system, that has also been employed in other studies [32,33,34]. Interestingly, in all of the fluidic experiments, the apical concentration measured at time t = 15 min was lower compared to the concentrations measured at t = 75 min and t = 120 min. This observation was assumed to be due to the fluidic experimental setup. After dose administration and additional rinsing with 0.5 mL, drug molecules were pumped from the medium reservoir in the direction of the apical part of the lung chamber and back to the apical dosing/sampling port. As a result, the concentrations at the dosing/sampling port were expected to increase in the first few minutes of the experiments until equilibrium was reached in the apical circuit. This assumption could be supported by in silico model simulations, where the concentrations in C2 rapidly increased after dose application, reaching its maximum concentration (Cmax) of ~99.5 µM after ~40 min. Thereafter, the concentrations in the apical circuit continuously decreased due to the transfer of drug molecules into the basolateral compartment.
As described by Sun et al. [43], the implementation of medium volume and drug mass removal might be important in in silico modeling approaches. In the fluidic control experiments, the pairwise observed concentrations at 22 h/24 h and 46 h/48 h showed a distinct stepwise decrease in concentrations even though the concentrations were measured only two hours apart. This stepwise reduction could be due to the decrease in volume during sampling. For example, the 0.5 mL sample at 48 h reduced the total apical volume by 4.2%. Hence, a correction of the volume and mass due to sampling was integrated into the in silico modeling process. Of note, the stepwise decrease was also observed in the fluidic cell experiment for caffeine at 46 h/48 h and for diclofenac at 22 h/24 h, while a stepwise increase was observed for caffeine at 22 h/24 h and for diclofenac at 46 h/48 h.
Some discrepancies between the measured and model-simulated concentrations could be observed. Simulated “steady-state” (t > 24 h) apical and basolateral diclofenac concentrations in the static experiments were slightly lower than the observed concentrations. This may be ascribed to analytical misspecifications, as the average sum of the measured amount of diclofenac in the apical and basolateral compartments at these time points was ~14% (control) and ~10% (cell experiment) higher than the actual amount of substance applied at the beginning of the experiment. Additionally, in the fluidic control experiments, the basolateral concentrations that were observed within the first two hours were slightly underpredicted, while the basolateral concentrations at later time points were predicted precisely.
The residual variability estimated from the static permeation experiments was carried over to the fluidic model without acknowledging potential additional system-dependent variability from the fluidic setup. Comparing fluidic model simulations and the observed concentrations, the predicted variability in the in silico model also covered most of the observed variability in the fluidic scenario, with some exceptions for the apical concentrations and caffeine basolateral concentrations in the control experiments, with some exceptions.
As previously described, the static permeation study was implemented with a bidirectional single-barrier model [11]. While more detailed models that incorporate additional compartments, representing unstirred water layers or cellular compartments, for example, have been used in recent studies [11,13,44], the employed two-compartment approach was a reasonable choice to describe the measured in vitro concentration–time profiles. Here, the single-barrier in silico model was selected, aiming for minimal sufficient model complexity to describe substance permeation. As a result of the employed model, the estimated apparent permeabilities represent composites of all permeation barriers involved, including the apical and basolateral unstirred water layers, the porous membrane, and the lung cell barrier if applicable [11].
In this study, a single-barrier approach was used to model substance permeation. Following this, in a second step, the intracellular pulmonary concentrations could be estimated and integrated in the in silico model. These would be helpful for the mechanistic implementation of active processes for substances that are subject to active transport, intracellular sequestration, or metabolism [10]. Here, potential metabolism of caffeine and diclofenac in the pulmonary epithelial cells seemed to be negligible, as the mean mass balance of the static cellular experiments showed a recovery that was close to the applied dose (39.4 µmol caffeine, 44.2 µmol diclofenac at steady state; applied dose: 40 µmol).
While physicochemical fluxes were integrated via the employed fluidic system, two other major stimuli acting on cells in vivo—biochemical stimuli from other cells (i.e., cytokines, hormones, etc.) and structural stimuli from the three-dimensional (3D) microenvironment [14,15]—could not be considered in this work. Studies have demonstrated more physiologically relevant behavior regarding morphology, cell–substrate interactions, and proliferation rates, among others, in 3D cell cultures [45,46]. In lung cells, differences have been observed regarding the maturation and functionality of the epithelial barrier, with lower TEER and reduced expression of tight junction proteins under 2D conditions [47]. Here, the 16HBE epithelial cell line was used as a proxy to mimic the pulmonary barrier in the lung chamber. Additionally, the study was performed under submerged conditions (apical and basolateral compartment embedded in a medium-filled fluid circuit) and not as an air–liquid interface model. The in vitro system could be adapted in future approaches, integrating, for example, cocultures of multiple cell types in a 3D system to increase physiologic relevance [14,15,48].
Furthermore, the integration of primary cells and additional tissues in multi-organ models could enhance physiologic characteristics in future studies [14,15,49]. Here, a second chamber (C8) was already integrated into the in vitro–in silico model of this study, providing the option to integrate additional tissues in the fluidic system. When experimental data from a fluidic system integrating lung cells and additional cell lines (e.g., hepatic cells) become available, the in vitro–in silico framework can be further evaluated for the prediction of these fluidic scenarios. With that, combination of in vitro–in silico models holds the potential to investigate tissue–tissue interactions in multi-organ models and to use modeling and simulation techniques to predict yet untested scenarios, such as drug−drug interactions or pharmacodynamic effects after the hepatic bioactivation of a prodrug molecule [7,25,29,31].

4. Materials and Methods

4.1. Standard Cell Cultivation

The bronchial cell line 16HBE was cultivated in Dulbecco’s modified Eagle medium (DMEM/F12 with L-glutamine and HEPES; Gibco 11039-021, Thermo Fisher Scientific, Waltham, MA, USA), supplemented with 10% fetal bovine serum (FBS, S0615-500mL, Sigma Aldrich, St. Louis, MO, USA), 100 U/mL penicillin, and 100 µg/mL streptomycin (Gibco 15140-122, Thermo Fisher Scientific). The cells were cultivated in a humidified atmosphere at 37 °C and 5% CO2. The bronchial cell line 16HBE was obtained from Dr. Dieter C. Gruenert (Cardiovascular Research Institute, University of California) [50].

4.2. Transepithelial Electric Resistance (TEER) Measurements

TEER measurements were performed with a cellZscope+ (Nanoanalytics GmbH, Münster, Germany). The TEER of the in vitro model represents a parameter for the integrity of the cell barrier and the functionality of the tight junctions [51]. Measurements were performed at 37 °C every hour for five days.

4.3. Cell Viability Assay

To verify that the fluidic conditions did not negatively affect cell viability, the viability of the 16HBE cells was determined by the alamarBlue assay (10099022, Invitrogen/Thermo Fisher Scientific) according to the manufacturer’s operating instructions three days after cultivating the cells under fluidic conditions (as described below). The fluorescence intensity was determined at ex/em 560 nm/610 nm. The viability of the cells that had been cultured under static conditions was set to 100% as a reference.

4.4. Sample Preparation for In Vitro Permeation Studies

For the in vitro permeation studies, a stock solution of caffeine (Sigma Aldrich #27602-250G) and diclofenac sodium salt (Sigma Aldrich #93484-100MG) was prepared in water (caffeine) and methanol (J.T. Baker, Phillipsburg, NJ, USA) (diclofenac), respectively, with a concentration of 1 mM. Solubilization of caffeine was enhanced using an ultrasonic bath (Elma S15 Elmasonic, Singen, Germany) for 10 min. For each in vitro study, the concentration of the applied solution was measured by HPLC to verify the solubility of the chemicals.

4.5. Permeation Studies under Static Conditions

Static permeation studies were performed in cell culture inserts with a pore size of 0.4 μm and a growth area of 0.6 cm2 (PIHP 012 50, 12 mm Diameter Millicell Culture Inserts, Merck Millipore, Burlington, MA, USA) (Figure 7a). For the in vitro cell permeation studies, 4.3 × 104 16HBE cells were seeded in the cell culture insert. 400 µL cell culture medium (CCM) was added on the apical side, and 600 µL was added on the basolateral side. The cells were cultivated in a 24-well plate (662160, Greiner bio-one, Frickenhausen, Germany) at 37 °C for five days under a physiological pH in the range of 7.0–7.4 with one CCM exchange after two days. Thereafter, CCM was removed, the stock solution was diluted 1:10 with CCM, and the cells were incubated with 100 µM caffeine or diclofenac on the apical side of the insert. The basolateral compartment was filled with fresh CCM. After model substance exposure, total apical, and basolateral CCM were collected each in a 1.5 mL tube (Kisker Biotech GmbH & Co. KG, Steinfurt, Germany) at the following eight time points: 15 min, 1 h 15 min, 2 h, 22 h, 24 h, 46 h, 48 h, and 72 h. For each time point, a separate cell culture insert was used. Samples were then stored at 4 °C until further analysis via HPLC.
In addition to the setup with 16HBE cells, control experiments (without cells) were performed with caffeine and diclofenac (porous membrane as permeation barrier) by following the same procedure described above for the in vitro cell-based study.

4.6. Permeation Studies under Fluidic Conditions

Fluidic permeation studies were performed in the MCmB QuasiVivo® QV600 (Kirkstall), which was assembled according to the manufacturer’s instructions. An apical and basolateral flow with a flow rate of 200 µL/min was used (Figure 7b). Continuous flow was achieved with a peristaltic pump (ISM597, Ismatec, Cole-Parmer GmbH, Wertheim, Germany). The apical circuit consisted of the apical reservoir bottle, the apical part of the pulmonary chamber, the pump, a three-way cock (B.Braun 4095111; dosing/sampling port; in order of medium flow), and connecting tubes. The basolateral circuit comprised the basolateral reservoir bottle, the basolateral part of the pulmonary chamber, a second chamber allowing the integration of additional cell tissues, the pump, the basolateral sampling port (in order of medium flow), and connecting tubes (Figure 8). Different tubing (Kirkstall QV TS-2 Tubing set; Cole-Parmer GZ06422-01; VWR 228-1066, Radnor, PA, USA) was used to establish the optimal connection between the reservoir bottles, cell chambers, pump, and sampling ports.
For the in vitro cell permeation studies, 4.3 × 104 16HBE cells were seeded in a cell culture insert with a pore size of 0.4 μm and a growth area of 0.6 cm2 (PIHP 012 50, Merck Millipore) and cultivated at 37 °C for five days. Thereafter, the cell culture medium was removed completely, and the insert was positioned in the fluidic QuasiVivo® QV600 (Kirkstall) setup. Both the apical and basolateral circuits were filled with CCM. To adapt the cells in the QV600, the cells were cultured for 30 min with a flow rate of 200 µL/min. For the application of the model substances (caffeine or diclofenac), the flow was stopped, and the stock solution of the substance was injected via the sampling port with a dilution resulting in 100 µM caffeine/diclofenac in the apical circuit. Dosing resulted in both the apical and basolateral circuits being filled with 15 mL CCM. Subsequently, the cells were cultivated for three days under fluidic conditions at a constant flow rate of 200 µL/min. At each of the eight time points (15 min–72 h), samples were collected by opening the sampling port, so that the media could drip from the circuit into a 1.5 mL tube (Kisker). After collecting 500 µL CCM from the apical circuit, the sampling port was closed, and the procedure was repeated for the basolateral circuit.
In addition to the setup with 16HBE cells, control experiments (without cells) were performed with caffeine and diclofenac (porous membrane as permeation barrier), following the same procedure as the one described above for the in vitro cell-based study.

4.7. High-Performance Liquid Chromatography (HPLC) Analysis

After the static and fluidic permeation studies, the model substances in the apical and basolateral CCM were quantified via HPLC measurements as described in Elberskirch et al. [52]. The amounts of caffeine and diclofenac were determined by an Agilent 1260 HPLC (Agilent Technologies, Santa Clara, CA, USA). Separation was performed using a Poroshell 120 EC-C18 column (2.1 × 100 mm, 2.7 µm; 695775-902, Agilent Technologies) connected to a guard column (Poroshell 120 EC-C18, 2.1 × 5 mm, 2.7 μm, Agilent Technologies). The mobile phase consisted of three running agents: water (solvent A; sterile filtered 0.2 µm pore), acetonitrile (Solvent B; J.T. Baker), and 0.1% trifluoroacetic acid (TFA) in water (Solvent C; Fisher Scientific). The gradient was programed as presented in Table 2. The flow rate was set to 0.4 mL/min with 5 µL of the sample injected. The column temperature was set to 60 °C (caffeine) and 30 °C (diclofenac), respectively. For detection, a diode array detector (DAD) was used at a wavelength of 275 nm.

4.8. Software

For in silico modeling, NONMEM® version 7.4 (Icon Development Solutions, Ellicott City, MD, USA) and the programming language R (version 3.6.3, R Foundation for Statistical Computing, Vienna, Austria [53]) with the mrgsolve package version 0.11.1 [54] were used. Further, R was employed for data processing and plot generation.

4.9. In Silico Model for the Static System

To model the static in vitro permeation experiments (control and cell experiments), a bidirectional single-barrier model was implemented as a two-compartment ODE system. Since an equilibrium with similar apical and basolateral concentrations was reached for both of the modeled compounds after ~24 h, permeabilities from the apical to the basolateral compartment and vice versa were assumed to be equal. The ODEs describing the changes in concentrations of the model substances in the apical and basolateral compartments over time are depicted in Equations (1) and (2)
d c 1 ( t ) d t = P 12 ,   i × A × ( c 2 c 1 ) V 1  
d c 2 ( t ) d t = P 12 ,   i × A × ( c 1 c 2 ) V 2
with the initial conditions c 1 ( 0 ) = D s V 1 and c 2 ( 0 ) = 0 .
Here, P12,i represents the permeability of compound i (i.e., caffeine or diclofenac) in the control (only porous membrane) and 16HBE cell experiment, respectively. Thereby, P12,i reflects the compound permeability through all present permeation barriers in the respective experimental setup, such as unstirred water layers, the porous filter membrane, and cell monolayer (apparent permeability) [11]. A represents the area of the filter membrane in the lung chamber (here, 0.6 cm2), c 1 represents the drug substance concentration in the apical compartment, and c 2 represents the concentration in the basolateral compartment. V1 (0.4 mL) and V2 (0.6 mL) represent the corresponding volumes of the compartments, while Ds is the dose applied to the apical compartment at time t = 0 (here, 40 nmol).

4.10. In Silico Model for the Fluidic System

Figure 8 depicts a schematic overview of the developed in silico model for the fluidic system. Here, C1 represents the apical part of the lung chamber and C6 the basolateral part. C2 comprises the volume of the sampling port and the tubes connecting the lung chamber and the apical sampling port, including the tubes that are clamped in the apical pump. C3 represents the tubes connecting the apical sampling port and the apical medium reservoir (C4), while C5 represents the tubes connecting the apical medium reservoir and the apical part of the lung chamber. The basolateral circuit comprises the basolateral part of the lung chamber (C6), a second chamber (C8) that provides the option to integrate additional cell tissue cultures in the fluidic system, the basolateral sampling port, including the tubes clamped in the basolateral pump (C9), a basolateral medium reservoir (C11), as well as connecting tubes (C7, C10, and C12).
In the in silico model for the fluidic system, the compartment volumes were set to the experimentally determined volumes of medium present in the respective compartment at t = 0: 1.5 mL (V1), 0.53 mL (V2), 0.53 mL (V3), 11.62 mL (V4), 0.82 mL (V5), 2 mL (V6), 0.44 mL (V7), 2.5 mL (V8), 0.49 mL (V9), 0.53 mL (V10), 8.49 mL (V11), and 0.55 mL (V12), with Vj representing the volume of compartment Cj. The compartments in the apical circuit and the basolateral circuit, respectively, were connected via medium flows (cf. Figure 8), and the flow rate was fixed to the experimentally determined value of 200 µL/min. The respective ODEs describing the changes in concentration of the model substances in the different compartments over time are depicted in Equations (3)–(5):
d c 1 ( t ) d t = P 16 , i × A × ( c 6 c 1 ) + Q × ( c 5 c 1 ) V 1
d c 6 ( t ) d t = P 16 , i × A × ( c 1 c 6 ) + Q × ( c 12 c 6 ) V 6
d c j ( t ) d t = Q × ( c j 1 c j ) V j   j { 2..5 ,   7..12 }
with the initial conditions c 4 ( 0 ) = D s V 4 and c j ( 0 ) = 0     j { 1..3 ,   5..12 } .
Here, P16,i is the apparent permeability of compound i (i.e., caffeine or diclofenac). A represents the area of the filter membrane in the lung chamber, which equals 0.6 cm2. cj is the concentration of the drug substance in compartment j, and Vj is the corresponding volume of the compartment. Moreover, Q represents the flow rate in both the apical and basolateral circuit (200 µL/min), and Ds is the dose of 1.5 µmol applied to the apical circuit at the beginning of the experiments. The dose was applied via a three-way cock followed by 0.5 mL CCM for rinsing. As a result, the applied dose was assumed to be located in the apical reservoir (C4) at the beginning of the fluidic experiments.
During sampling, the medium flow was redirected from C2 (apical) and C9 (basolateral) to a sampling container. Due to continuous flow, the volume of the sampling compartments refilled continuously, allowing the removal of the sample volume. At each sampling time point, 500 µL of medium was extracted from both the apical and basolateral circuits without replacement, reducing their volumes over time. This stepwise decrease in volume in the two reservoirs (V4 and V11) at each of the sampling time points was incorporated in the in silico model as a zero-order process over the sampling period of 2.5 min.

4.11. Model Parameter Estimation and Simulation

During parameter estimation in the static permeation study, permeability parameters (P12,i) were estimated by fitting the model-simulated concentrations to the observed experimental data via first-order conditional estimation method with interaction (FOCE-I) in NONMEM®, using a combined proportional and additive residual error model (Equation (6)).
Y i j = c i j × ( 1 + ε p i j ) + ε a i j
Here, Yij represents the jth observed concentration for compound i (i.e., caffeine or diclofenac) in the apical or basolateral compartment. c i j is the corresponding model estimated concentration of the drug substance, and ε p i j and ε a i j are the residual errors for the proportional and the additive components of the model with means of zero and variances of σ a 2 and σ p 2 , respectively (ε~N [0, σ 2 ]).
For the stochastic model simulations, the ODEs for the static and fluidic system were implemented in R within the mrgsolve framework and used to estimate the 68% prediction intervals (n = 500 replicates). Here, the estimated values for the proportional and additive residual errors from the in silico model for the static experiments were carried over to the in silico model for the fluidic system.
For refinement of model predictions in the fluidic cell experiments, a factor was estimated with the model implemented in mrgsolve and multiplied with the obtained apparent permeability values from the static in vitro experiments. The factor was estimated via maximum likelihood estimation, assuming a combined proportional and additive error model, minimizing Equation (7) with the L-BFGS-B optimizer from the nloptr package [55,56,57].
2 L L = i = 1 n ( ( y ^ y ) 2 ( y ^ 2   σ p 2 + σ a 2 ) + l o g ( y ^ 2 σ p 2 + σ a 2 ) ) + n l o g ( 2 π )
Here, y represents the observed concentration, y ^ is the corresponding predicted concentration, and n is the number of observed concentrations. Further, σ a 2 represents the variance of the additive, and σ p 2 represents the variance of the proportional residual error.

4.12. Statistical Analysis

In vitro measurements are presented as arithmetic mean with standard deviation of three independent experiments (n = 3) unless mentioned otherwise. Cell viability under fluidic conditions was compared to cell viability under static conditions using an unpaired two-sample t-test performed in Excel (Microsoft Office Professional Plus 2016, Microsoft Corporation, Redmond, WA, USA).

5. Conclusions

In summary, this work presents a case example of in silico modeling to describe and predict the concentration–time course of two model substances in a fluidic in vitro permeation system. Here, the estimated apparent permeabilities from static in vitro experiments could be used to predict the behavior of the compounds in the fluidic system under control conditions (with the filter membrane as permeation barrier). In contrast, a decreased permeability in the fluidic setting compared to the static experimental design had to be estimated for a precise description of the observed caffeine and diclofenac concentrations in the fluidic cell experiments with 16HBE cells. The developed models can be further expanded and employed to model the kinetics and the drug effects of pharmacologically active substances in future studies, including experimental scenarios with multiple tissues integrated into the fluidic system. With that, the in vitro–in silico models can contribute to a reduction in preclinical animal studies in line with the 3R principle.

Author Contributions

Conceptualization, L.K., L.W., D.S., Y.K., R.B. and T.L.; funding acquisition, Y.K., R.B. and T.L.; software, L.K. and D.S.; investigation, L.K., L.W., D.S., Y.K. and T.L.; methodology, L.K., L.W., D.S., Y.K. and T.L.; visualization, L.K., L.W. and D.S.; writing—original draft, L.K., L.W., D.S., Y.K. and T.L.; writing—review and editing, L.K., L.W., D.S., Y.K., R.B. and T.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Federal Ministry of Education and Research (BMBF), 031L0153 “Alternativmethoden zum Tierversuch”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Directive 2010/63/EU of the European Parliament and of the Coundcil of 22 September 2010 on the Protection of Animals used for Scientific Purposes. Off. J. Eur. Union 2010, L276, 33–79.
  2. Kirk, R.G.W. Recovering The Principles of Humane Experimental Technique. Sci. Technol. Hum. Values 2018, 43, 622–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Movia, D.; Prina-Mello, A. Preclinical Development of Orally Inhaled Drugs (OIDs)—Are Animal Models Predictive or Shall We Move Towards In Vitro Non-Animal Models? Animals 2020, 10, 1259. [Google Scholar] [CrossRef] [PubMed]
  4. Togami, K.; Chono, S.; Morimoto, K. Transport characteristics of clarithromycin, azithromycin and telithromycin, antibiotics applied for treatment of respiratory infections, in Calu-3 cell monolayers as model lung epithelial cells. Pharmazie 2012, 67, 389–393. [Google Scholar] [CrossRef]
  5. Mathia, N.R.; Timoszyk, J.; Stetsko, P.I.; Megill, J.R.; Smith, R.L.; Wall, D.A. Permeability characteristics of calu-3 human bronchial epithelial cells: In vitro-in vivo correlation to predict lung absorption in rats. J. Drug Target. 2002, 10, 31–40. [Google Scholar] [CrossRef] [PubMed]
  6. Foster, K.A.; Avery, M.L.; Yazdanian, M.; Audus, K.L. Characterization of the Calu-3 cell line as a tool to screen pulmonary drug delivery. Int. J. Pharm. 2000, 208, 1–11. [Google Scholar] [CrossRef]
  7. Tatosian, D.A.; Shuler, M.L. A novel system for evaluation of drug mixtures for potential efficacy in treating multidrug resistant cancers. Biotechnol. Bioeng. 2009, 103, 187–198. [Google Scholar] [CrossRef]
  8. Balimane, P.V.; Chong, S. Cell culture-based models for intestinal permeability: A critique. Drug Discov. Today 2005, 10, 335–343. [Google Scholar] [CrossRef]
  9. Heikkinen, A.T.; Korjamo, T.; Lepikkö, V.; Mönkkönen, J. Effects of experimental setup on the apparent concentration dependency of active efflux transport in in vitro cell permeation experiments. Mol. Pharm. 2010, 7, 605–617. [Google Scholar] [CrossRef]
  10. Heikkinen, A.T.; Mönkkönen, J.; Korjamo, T. Determination of permeation resistance distribution in in vitro cell monolayer permeation experiments. Eur. J. Pharm. Sci. 2010, 40, 132–142. [Google Scholar] [CrossRef]
  11. Heikkinen, A.T.; Korjamo, T.; Mönkkönen, J. Modelling of Drug Disposition Kinetics in In Vitro Intestinal Absorption Cell Models. Basic Clin. Pharmacol. Toxicol. 2010, 106, 180–188. [Google Scholar] [CrossRef] [PubMed]
  12. Lozoya-Agullo, I.; Araújo, F.; González-Álvarez, I.; Merino-Sanjuán, M.; González-Álvarez, M.; Bermejo, M.; Sarmento, B. Usefulness of Caco-2/HT29-MTX and Caco-2/HT29-MTX/Raji B Coculture Models To Predict Intestinal and Colonic Permeability Compared to Caco-2 Monoculture. Mol. Pharm. 2017, 14, 1264–1270. [Google Scholar] [CrossRef] [PubMed]
  13. Heikkinen, A.T.; Mönkkönen, J.; Korjamo, T. Kinetics of cellular retention during Caco-2 permeation experiments: Role of lysosomal sequestration and impact on permeability estimates. J. Pharmacol. Exp. Ther. 2009, 328, 882–892. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Mattei, G.; Giusti, S.; Ahluwalia, A. Design criteria for generating physiologically relevant in vitro models in bioreactors. Processes 2014, 2, 548–569. [Google Scholar] [CrossRef] [Green Version]
  15. Di Nardo, P.; Minieri, M.; Ahluwalia, A. Engineering the Stem Cell Niche and the Differentiative Micro- and Macroenvironment: Technologies and Tools for Applying Biochemical, Physical and Structural Stimuli and Their Effects on Stem Cells. In Stem Cell Engineering; Artmann, G.M., Minger, S., Hescheler, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; pp. 41–59. ISBN 978-3-642-11864-7. [Google Scholar]
  16. Sbrana, T.; Ahluwalia, A. Engineering quasi-vivo® in vitro organ models. Adv. Exp. Med. Biol. 2012, 745, 138–153. [Google Scholar] [CrossRef]
  17. Chan, C.Y.; Huang, P.-H.; Guo, F.; Ding, X.; Kapur, V.; Mai, J.D.; Yuen, P.K.; Huang, T.J. Accelerating drug discovery via organs-on-chips. Lab Chip 2013, 13, 4697–4710. [Google Scholar] [CrossRef] [Green Version]
  18. Alborzinia, H.; Can, S.; Holenya, P.; Scholl, C.; Lederer, E.; Kitanovic, I.; Wölfl, S. Real-time monitoring of cisplatin-induced cell death. PLoS ONE 2011, 6, e19714. [Google Scholar] [CrossRef]
  19. Esch, E.W.; Bahinski, A.; Huh, D. Organs-on-chips at the frontiers of drug discovery. Nat. Rev. Drug Discov. 2015, 14, 248–260. [Google Scholar] [CrossRef] [Green Version]
  20. Mazzei, D.; Guzzardi, M.A.; Giusti, S.; Ahluwalia, A. A low shear stress modular bioreactor for connected cell culture under high flow rates. Biotechnol. Bioeng. 2010, 106, 127–137. [Google Scholar] [CrossRef]
  21. Theobald, J.; Abu El Maaty, M.A.; Kusterer, N.; Wetterauer, B.; Wink, M.; Cheng, X.; Wölfl, S. In vitro metabolic activation of vitamin D3 by using a multi-compartment microfluidic liver-kidney organ on chip platform. Sci. Rep. 2019, 9, 4616. [Google Scholar] [CrossRef]
  22. Theobald, J.; Ghanem, A.; Wallisch, P.; Banaeiyan, A.A.; Andrade-Navarro, M.A.; Taškova, K.; Haltmeier, M.; Kurtz, A.; Becker, H.; Reuter, S.; et al. Liver-Kidney-on-Chip To Study Toxicity of Drug Metabolites. ACS Biomater. Sci. Eng. 2018, 4, 78–89. [Google Scholar] [CrossRef] [PubMed]
  23. Schimek, K.; Frentzel, S.; Luettich, K.; Bovard, D.; Rütschle, I.; Boden, L.; Rambo, F.; Erfurth, H.; Dehne, E.-M.; Winter, A.; et al. Human multi-organ chip co-culture of bronchial lung culture and liver spheroids for substance exposure studies. Sci. Rep. 2020, 10, 7865. [Google Scholar] [CrossRef] [PubMed]
  24. Sung, J.H.; Shuler, M.L. A micro cell culture analog (microCCA) with 3-D hydrogel culture of multiple cell lines to assess metabolism-dependent cytotoxicity of anti-cancer drugs. Lab Chip 2009, 9, 1385–1394. [Google Scholar] [CrossRef] [PubMed]
  25. Low, L.A.; Mummery, C.; Berridge, B.R.; Austin, C.P.; Tagle, D.A. Organs-on-chips: Into the next decade. Nat. Rev. Drug Discov. 2021, 20, 345–361. [Google Scholar] [CrossRef]
  26. Marshall, S.; Madabushi, R.; Manolis, E.; Krudys, K.; Staab, A.; Dykstra, K.; Visser, S.A.G. Model-Informed Drug Discovery and Development: Current Industry Good Practice and Regulatory Expectations and Future Perspectives. CPT Pharmacomet. Syst. Pharmacol. 2019, 8, 87–96. [Google Scholar] [CrossRef]
  27. EFPIA MID3 Workgroup; Marshall, S.F.; Burghaus, R.; Cosson, V.; Cheung, S.Y.A.; Chenel, M.; DellaPasqua, O.; Frey, N.; Hamrén, B.; Harnisch, L.; et al. Good Practices in Model-Informed Drug Discovery and Development: Practice, Application, and Documentation. CPT Pharmacomet. Syst. Pharmacol. 2016, 5, 93–122. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, Y.; Zhu, H.; Madabushi, R.; Liu, Q.; Huang, S.; Zineh, I. Model-Informed Drug Development: Current US Regulatory Practice and Future Considerations. Clin. Pharmacol. Ther. 2019, 105, 899–911. [Google Scholar] [CrossRef]
  29. Sung, J.H.; Esch, M.B.; Shuler, M.L. Integration of in silico and in vitro platforms for pharmacokinetic–pharmacodynamic modeling. Expert Opin. Drug Metab. Toxicol. 2010, 6, 1063–1081. [Google Scholar] [CrossRef]
  30. Sung, J.H.; Srinivasan, B.; Esch, M.B.; McLamb, W.T.; Bernabini, C.; Shuler, M.L.; Hickman, J.J. Using physiologically-based pharmacokinetic-guided “body-on-a-chip” systems to predict mammalian response to drug and chemical exposure. Exp. Biol. Med. 2014, 239, 1225–1239. [Google Scholar] [CrossRef] [Green Version]
  31. Edington, C.D.; Chen, W.L.K.; Geishecker, E.; Kassis, T.; Soenksen, L.R.; Bhushan, B.M.; Freake, D.; Kirschner, J.; Maass, C.; Tsamandouras, N.; et al. Interconnected Microphysiological Systems for Quantitative Biology and Pharmacology Studies. Sci. Rep. 2018, 8, 4530. [Google Scholar] [CrossRef] [Green Version]
  32. Giusti, S.; Sbrana, T.; La Marca, M.; Di Patria, V.; Martinucci, V.; Tirella, A.; Domenici, C.; Ahluwalia, A. A novel dual-flow bioreactor simulates increased fluorescein permeability in epithelial tissue barriers. Biotechnol. J. 2014, 9, 1175–1184. [Google Scholar] [CrossRef] [PubMed]
  33. Chandorkar, P.; Posch, W.; Zaderer, V.; Blatzer, M.; Steger, M.; Ammann, C.G.; Binder, U.; Hermann, M.; Hörtnagl, P.; Lass-Flörl, C.; et al. Fast-track development of an in vitro 3D lung/immune cell model to study Aspergillus infections. Sci. Rep. 2017, 7, 11644. [Google Scholar] [CrossRef] [Green Version]
  34. Martin, K.C.; Yuan, X.; Stimac, G.; Bannerman, K.; Anderson, J.; Roy, C.; Glykofrydis, F.; Yin, H.; Davies, J.A. Symmetry-breaking in branching epithelia: Cells on micro-patterns under flow challenge the hypothesis of positive feedback by a secreted autocrine inhibitor of motility. J. Anat. 2017, 230, 766–774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Callaghan, P.J.; Ferrick, B.; Rybakovsky, E.; Thomas, S.; Mullin, J.M. Epithelial barrier function properties of the 16HBE14o- human bronchial epithelial cell culture model. Biosci. Rep. 2020, 40, BSR20201532. [Google Scholar] [CrossRef] [PubMed]
  36. González-Alvarez, I.; Fernández-Teruel, C.; Garrigues, T.M.; Casabo, V.G.; Ruiz-García, A.; Bermejo, M. Kinetic modelling of passive transport and active efflux of a fluoroquinolone across Caco-2 cells using a compartmental approach in NONMEM. Xenobiotica 2005, 35, 1067–1088. [Google Scholar] [CrossRef]
  37. DrugBank, Entry for Caffeine. Available online: https://go.drugbank.com/drugs/DB00201 (accessed on 7 October 2021).
  38. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef]
  39. DrugBank, Entry for Diclofenac. Available online: https://go.drugbank.com/drugs/DB00586 (accessed on 7 October 2021).
  40. Prot, J.M.; Aninat, C.; Griscom, L.; Razan, F.; Brochot, C.; Guillouzo, C.G.; Legallais, C.; Corlu, A.; Leclerc, E. Improvement of HepG2/C3a cell functions in a microfluidic biochip. Biotechnol. Bioeng. 2011, 108, 1704–1715. [Google Scholar] [CrossRef]
  41. Baudoin, R.; Griscom, L.; Prot, J.M.; Legallais, C.; Leclerc, E. Behavior of HepG2/C3A cell cultures in a microfluidic bioreactor. Biochem. Eng. J. 2011, 53, 172–181. [Google Scholar] [CrossRef]
  42. Snouber, L.C.; Letourneur, F.; Chafey, P.; Broussard, C.; Monge, M.; Legallais, C.; Leclerc, E. Analysis of transcriptomic and proteomic profiles demonstrates improved Madin-Darby canine kidney cell function in a renal microfluidic biochip. Biotechnol. Prog. 2012, 28, 474–484. [Google Scholar] [CrossRef]
  43. Sun, H.; Zhang, L.; Chow, E.C.Y.; Lin, G.; Zuo, Z.; Pang, K.S. A catenary model to study transport and conjugation of baicalein, a bioactive flavonoid, in the Caco-2 cell monolayer: Demonstration of substrate inhibition. J. Pharmacol. Exp. Ther. 2008, 326, 117–126. [Google Scholar] [CrossRef] [Green Version]
  44. Kalvass, J.C.; Pollack, G.M. Kinetic considerations for the quantitative assessment of efflux activity and inhibition: Implications for understanding and predicting the effects of efflux inhibition. Pharm. Res. 2007, 24, 265–276. [Google Scholar] [CrossRef] [PubMed]
  45. Jensen, C.; Teng, Y. Is It Time to Start Transitioning From 2D to 3D Cell Culture? Front. Mol. Biosci. 2020, 7, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Mirbagheri, M.; Adibnia, V.; Hughes, B.R.; Waldman, S.D.; Banquy, X.; Hwang, D.K. Advanced cell culture platforms: A growing quest for emulating natural tissues. Mater. Horizons 2019, 6, 45–71. [Google Scholar] [CrossRef]
  47. Bhowmick, R.; Gappa-Fahlenkamp, H. Cells and Culture Systems Used to Model the Small Airway Epithelium. Lung 2016, 194, 419–428. [Google Scholar] [CrossRef]
  48. van Midwoud, P.M.; Verpoorte, E.; Groothuis, G.M.M. Microfluidic devices for in vitro studies on liver drug metabolism and toxicity. Integr. Biol. 2011, 3, 509–521. [Google Scholar] [CrossRef]
  49. Pedersen, J.M.; Shim, Y.-S.; Hans, V.; Phillips, M.B.; Macdonald, J.M.; Walker, G.; Andersen, M.E.; Clewell, H.J.; Yoon, M. Fluid Dynamic Modeling to Support the Development of Flow-Based Hepatocyte Culture Systems for Metabolism Studies. Front. Bioeng. Biotechnol. 2016, 4, 72. [Google Scholar] [CrossRef] [Green Version]
  50. Cozens, A.L.; Yezzi, M.J.; Kunzelmann, K.; Ohrui, T.; Chin, L.; Eng, K.; Finkbeiner, W.E.; Widdicombe, J.H.; Gruenert, D.C. CFTR expression and chloride secretion in polarized immortal human bronchial epithelial cells. Am. J. Respir. Cell Mol. Biol. 1994, 10, 38–47. [Google Scholar] [CrossRef]
  51. Srinivasan, B.; Kolli, A.R.; Esch, M.B.; Abaci, H.E.; Shuler, M.L.; Hickman, J.J. TEER Measurement Techniques for In Vitro Barrier Model Systems. J. Lab. Autom. 2015, 20, 107–126. [Google Scholar] [CrossRef] [Green Version]
  52. Elberskirch, L.; Knoll, T.; Moosmann, A.; Wilhelm, N.; von Briesen, H.; Wagner, S. A novel microfluidic mucus-chip for studying the permeation of compounds over the mucus barrier. J. Drug Deliv. Sci. Technol. 2019, 54, 101248. [Google Scholar] [CrossRef]
  53. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  54. Elmokadem, A.; Riggs, M.M.; Baron, K.T. Quantitative Systems Pharmacology and Physiologically-Based Pharmacokinetic Modeling With mrgsolve: A Hands-On Tutorial. CPT Pharmacomet. Syst. Pharmacol. 2019, 8, 883–893. [Google Scholar] [CrossRef] [Green Version]
  55. Johnson, S.G. The NLopt Nonlinear-Optimization Package. Available online: https://github.com/stevengj/nlopt (accessed on 14 October 2021).
  56. Nocedal, J. Updating Quasi-Newton Matrices with Limited Storage. Math. Comput. 1980, 35, 773. [Google Scholar] [CrossRef]
  57. Liu, D.C.; Nocedal, J. On the limited memory BFGS method for large scale optimization. Math. Program. 1989, 45, 503–528. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the static permeation study without lung cells (control). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 5 for caffeine, n = 3 for diclofenac). Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Figure 1. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the static permeation study without lung cells (control). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 5 for caffeine, n = 3 for diclofenac). Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Pharmaceuticals 15 00250 g001
Figure 2. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the static cell permeation study. Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 5 for caffeine, n = 3 for diclofenac). Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Figure 2. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the static cell permeation study. Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 5 for caffeine, n = 3 for diclofenac). Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Pharmaceuticals 15 00250 g002
Figure 3. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the fluidic permeation study without lung cells (control). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and predicted. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 3). Model-predicted concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartment. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Figure 3. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the fluidic permeation study without lung cells (control). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and predicted. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 3). Model-predicted concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartment. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Pharmaceuticals 15 00250 g003
Figure 4. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the fluidic cell permeation study without (a,b) and with adjusted permeabilities (c,d). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 4 for caffeine, n = 3 for diclofenac). Model-predicted concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartments. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval. F, estimated factor to adjust the apparent permeabilities in the fluidic compared to the static setting.
Figure 4. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the fluidic cell permeation study without (a,b) and with adjusted permeabilities (c,d). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 4 for caffeine, n = 3 for diclofenac). Model-predicted concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartments. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval. F, estimated factor to adjust the apparent permeabilities in the fluidic compared to the static setting.
Pharmaceuticals 15 00250 g004
Figure 5. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the fluidic control permeation study with (grey) and without (orange) considering volume and substance removal. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 3). Model-predicted concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartments. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Figure 5. Apical and basolateral concentration–time profiles for caffeine (a) and diclofenac (b) in the fluidic control permeation study with (grey) and without (orange) considering volume and substance removal. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 3). Model-predicted concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartments. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval.
Pharmaceuticals 15 00250 g005
Figure 6. Zoom into the first 2.5 h of apical and basolateral concentration–time profiles for caffeine (a,c) and diclofenac (b,d) in the fluidic permeation study without lung cells (control; a,b) and with lung cells and adjusted permeabilities (c,d). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 3 in (a,b,d) and n = 4 in (c)). Model-simulated concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartments. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval. F, estimated factor to adjust the apparent permeabilities in the fluidic compared to the static setting.
Figure 6. Zoom into the first 2.5 h of apical and basolateral concentration–time profiles for caffeine (a,c) and diclofenac (b,d) in the fluidic permeation study without lung cells (control; a,b) and with lung cells and adjusted permeabilities (c,d). Darker colors show apical concentrations, while lighter colors show basolateral concentrations, both observed and simulated. Observed concentrations are depicted as solid circles (apical) and open triangles (basolateral) (arithmetic mean and standard deviation, n = 3 in (a,b,d) and n = 4 in (c)). Model-simulated concentrations represent concentrations in the C2 (apical) and C9 (basolateral) compartments. Model simulations are shown as lines, and shaded areas depict the 68% simulation prediction interval. F, estimated factor to adjust the apparent permeabilities in the fluidic compared to the static setting.
Pharmaceuticals 15 00250 g006
Figure 7. Schematic overview of the static (a) and the fluidic (b) setup for in vitro permeation studies. Under static conditions an insert was placed in a 24-well plate. The test substance was applied to the apical side (1) and permeated to the basolateral side (2) over time. Under fluidic conditions an insert was placed in a fluidic chamber (QuasiVivo® QV600, Kirkstall) and connected to a fluidic circuit. The test substance was applied to the apical circulation and permeated through the barrier to the basolateral side over time. Black arrows represent the apical (1) and basolateral (2) circuits. Both circuits were connected to a peristaltic pump that ensured a continuous flow of 200 µL/min during the experiment. Pink arrows depict permeation through the barrier.
Figure 7. Schematic overview of the static (a) and the fluidic (b) setup for in vitro permeation studies. Under static conditions an insert was placed in a 24-well plate. The test substance was applied to the apical side (1) and permeated to the basolateral side (2) over time. Under fluidic conditions an insert was placed in a fluidic chamber (QuasiVivo® QV600, Kirkstall) and connected to a fluidic circuit. The test substance was applied to the apical circulation and permeated through the barrier to the basolateral side over time. Black arrows represent the apical (1) and basolateral (2) circuits. Both circuits were connected to a peristaltic pump that ensured a continuous flow of 200 µL/min during the experiment. Pink arrows depict permeation through the barrier.
Pharmaceuticals 15 00250 g007
Figure 8. Schematic overview of the in silico model compartments for the fluidic system. Dark blue boxes represent compartments of the apical circuit, and light blue boxes represent compartments of the basolateral circuit. Gray solid lines depict medium flow. The blue dashed line represents the porous membrane without (control) or with the 16HBE cellular barrier. C1 and C6: apical and basolateral compartments of the lung chamber; C8: second chamber; C2 and C9: apical and basolateral sampling ports including preceding tubes; C4 and C11: apical and basolateral medium reservoirs; C3 C5, C7, C10 and C12: connecting tubes.
Figure 8. Schematic overview of the in silico model compartments for the fluidic system. Dark blue boxes represent compartments of the apical circuit, and light blue boxes represent compartments of the basolateral circuit. Gray solid lines depict medium flow. The blue dashed line represents the porous membrane without (control) or with the 16HBE cellular barrier. C1 and C6: apical and basolateral compartments of the lung chamber; C8: second chamber; C2 and C9: apical and basolateral sampling ports including preceding tubes; C4 and C11: apical and basolateral medium reservoirs; C3 C5, C7, C10 and C12: connecting tubes.
Pharmaceuticals 15 00250 g008
Table 1. Parameter estimates for caffeine and diclofenac in the static permeation study.
Table 1. Parameter estimates for caffeine and diclofenac in the static permeation study.
ParameterValueRSE (%)
Control experiments (without lung cell barrier)
P12, caffeine (cm/min)3.27 × 10−30.1
P12, diclofenac (cm/min)2.51 × 10−30
σ a 2 (µM)21.435
σ p 2 0 (fixed)-
Cell experiments (with lung cell barrier)
P12, caffeine (cm/min)2.62 × 10−30.7
P12, diclofenac (cm/min)1.38 × 10−30.9
σ a 2 (µM)2.3830
σ p 2 0.0059623
P12, apparent permeability from the apical to the basolateral compartment and vice versa; RSE, relative standard error for the fitted parameter; σ a 2 , variance of the additive residual error; σ p 2 , variance of the proportional residual error.
Table 2. HPLC gradient conditions according to Elberskirch et al. [51]. Solvent A: water; solvent B: acetonitrile; solvent C: 0.1% trifluoroacetic acid in water.
Table 2. HPLC gradient conditions according to Elberskirch et al. [51]. Solvent A: water; solvent B: acetonitrile; solvent C: 0.1% trifluoroacetic acid in water.
Time (min)Solvent A (%)Solvent B (%)Solvent C (%)
076.58.515
0.576.58.515
2.5423820
4057.542.5
10057.542.5
1476.58.515
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kovar, L.; Wien, L.; Selzer, D.; Kohl, Y.; Bals, R.; Lehr, T. In Vitro–In Silico Modeling of Caffeine and Diclofenac Permeation in Static and Fluidic Systems with a 16HBE Lung Cell Barrier. Pharmaceuticals 2022, 15, 250. https://doi.org/10.3390/ph15020250

AMA Style

Kovar L, Wien L, Selzer D, Kohl Y, Bals R, Lehr T. In Vitro–In Silico Modeling of Caffeine and Diclofenac Permeation in Static and Fluidic Systems with a 16HBE Lung Cell Barrier. Pharmaceuticals. 2022; 15(2):250. https://doi.org/10.3390/ph15020250

Chicago/Turabian Style

Kovar, Lukas, Lena Wien, Dominik Selzer, Yvonne Kohl, Robert Bals, and Thorsten Lehr. 2022. "In Vitro–In Silico Modeling of Caffeine and Diclofenac Permeation in Static and Fluidic Systems with a 16HBE Lung Cell Barrier" Pharmaceuticals 15, no. 2: 250. https://doi.org/10.3390/ph15020250

APA Style

Kovar, L., Wien, L., Selzer, D., Kohl, Y., Bals, R., & Lehr, T. (2022). In Vitro–In Silico Modeling of Caffeine and Diclofenac Permeation in Static and Fluidic Systems with a 16HBE Lung Cell Barrier. Pharmaceuticals, 15(2), 250. https://doi.org/10.3390/ph15020250

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop