Abstract
Purpose
Validate and exemplify a discrete, componentized, in silico, transwell device (ISTD) capable of mimicking the in vitro passive transport properties of compounds through cell monolayers. Verify its use for studying drug–drug interactions.
Methods
We used the synthetic modeling method. Specialized software components represented spatial and functional features including cell components, semi-porous tight junctions, and metabolizing enzymes. Mobile components represented drugs. Experiments were conducted and analyzed as done in vitro.
Results
Verification experiments provided data analogous to those in the literature. ISTD parameters were tuned to simulate and match in vitro urea transport data; the objects representing tight junction (effective radius of 6.66 Å) occupied 0.066% of the surface area. That ISTD was then tuned to simulate pH-dependent, in vitro alfentanil transport properties. The resulting ISTD predicted the passive transport properties of 14 additional compounds, individually and all together in one in silico experiment. The function of a two-site enzymatic component was cross-validated with a kinetic model and then experimentally validated against in vitro benzyloxyresorufin metabolism data. Those components were used to exemplify drug–drug interaction studies.
Conclusions
The ISTD is an example of a new class of simulation models capable of realistically representing complex drug transport and drug–drug interaction phenomena.
Similar content being viewed by others
Abbreviations
- CYP:
-
a P450 enzyme
- ISTD:
-
in silico transwell device
- MW:
-
molecular weight
- PCP:
-
physicochemical properties
- PRN:
-
pseudorandom number
- SM:
-
supplementary material
- TD:
-
transwell device
References
Y. Liu and C. A. Hunt. Studies of intestinal drug transport using an in silico epithelio-mimetic device. Biosystems 82:154–167 (2005).
Y. Liu and C. A. Hunt. Mechanistic study of the cellular interplay of transport and metabolism using the synthetic modeling method. Pharm. Res. 23:493–505 (2006).
G. E. Ropella, C. A. Hunt, and D. A. Nag. Using heuristic models to bridge the gap between analytic and experimental models in biology. In L. Yilmaz (ed), Agent-Directed Simulation Symposium, Vol. Simulation Series 37(2) (ADS’05). SCS Press, San Diego, CA, 2005, pp. 182–190.
G. E. P. Ropella, C. A. Hunt, and S. Sheikh-Bahaei. Methodological considerations of heuristic modeling of biological systems. 9th World Multi-Conference on Systemics, Cybernetics and Informatics, Orlando, FL, 2005.
C. A. Hunt, G. E. Ropella, L. Yan, D. Y. Hung, and M. S. Roberts. Physiologically based synthetic models of hepatic disposition. J. Pharmacokinet. Pharmacodyn. 33:737–772 (2006).
S. Sheikh-Bahaei and C. A. Hunt. Prediction of in vitro hepatic biliary excretion using stochastic agent-based modeling and fuzzy clustering. In: L. F. Perrone et al. (eds.), Proc. of the 37th Conference on Winter Simulation, 1617–24 (2006).
S. Railsback et al. Agent-based simulation platforms: review and development recommendations. Simulation 82:609–623 (2006).
G. T. Knipp, N. F. Ho, C. L. Barsuhn, and R. T. Borchardt. Paracellular diffusion in Caco-2 cell monolayers: effect of perturbation on the transport of hydrophilic compounds that vary in charge and size. J. Pharm. Sci. 86:1105–1110 (1997).
R. Hine. Membrane, The Facts on File Dictionary of Biology, Checkmark, New York, 1999, pp. 198.
P. Ball. Introduction to discrete event simulation, 2nd DYCOMANS workshop on "Management and Control: Tools in Action,” Algarve, Portugal, 1996, pp. 367–76.
T. H. Cormen, C. E. Leiserson, and R. L. Rivest. Introduction to algorithms, McGraw-Hill, New York, 1990.
G. Camenisch, G. Folkers, and H. van de Waterbeemd. Review of theoretical passive drug absorption models: historical background, recent developments and limitations. Pharm. Acta Helv. 71:309–327 (1996).
K. Palm, K. Luthman, J. Ros, J. Grasjo, and P. Artursson. Effect of molecular charge on intestinal epithelial drug transport: pH-dependent transport of cationic drugs. J. Pharmacol. Exp. Ther. 291:435–443 (1999).
K. Obata, K. Sugano, R. Saitoh, A. Higashida, Y. Nabuchi, M. Machida, and Y. Aso. Prediction of oral drug absorption in humans by theoretical passive absorption model. Int. J. Pharm. 293:183–192 (2005).
S. Narasimhulu. Differential behavior of the sub-sites of cytochrome 450 active site in binding of substrates, and products (implications for coupling/uncoupling). Biochim. Biophys. Acta 1770:360–375 (2007).
W. M. Atkins. Non-Michaelis–Menten kinetics in cytochrome P450-catalyzed reactions. Annu. Rev. Pharmacol. Toxicol. 45:291–310 (2005).
M. Shou, Y. Lin, P. Lu, C. Tang, Q. Mei, D. Cui, W. Tang, J. S. Ngui, C. C. Lin, R. Singh, B. K. Wong, J. A. Yergey, J. H. Lin, P. G. Pearson, T. A. Baillie, A. D. Rodrigues, and T. H. Rushmore. Enzyme kinetics of cytochrome P450-mediated reactions. Current Drug Metabolism 2:17–36 (2001).
K. R. Korzekwa, N. Krishnamachary, M. Shou, A. Ogai, R. A. Parise, A. E. Rettie, F. J. Gonzalez, and T. S. Tracy. Evaluation of atypical cytochrome P450 kinetics with two-substrate models: evidence that multiple substrates can simultaneously bind to cytochrome P450 active sites. Biochemistry 37:4137–4147 (1998).
G. Camenisch, G. Folkers, and H. van de Waterbeemd. Shapes of membrane permeability-lipophilicity curves: extension of theoretical models with an aqueous pore pathway. Eur. J. Pharm. Sci. 6:325–329 (1998).
H. Lennernas. Does fluid flow across the intestinal mucosa affect quantitative oral drug absorption? Is it time for a reevaluation? Pharm. Res. 12:1573–1582 (1995).
P. Artursson, K. Palm, and K. Luthman. Caco-2 monolayers in experimental and theoretical predictions of drug transport. Adv. Drug Deliv. Rev. 46:27–43 (2001).
Defense Modeling and Simulation Office, United States Department of Defense. Online M&S Glossary. https://www.dmso.mil/public/resources/glossary/, Cited 3/27/07. Part of the Defense Modeling and Simulation Office Web site, https://www.dmso.mil/public/, Cited 3/27/07.
A. Avdeef. Physicochemical profiling (solubility, permeability and charge state). Curr. Top. Med. Chem. 1:277–351 (2001).
A. Barat, H. J. Ruskin, and M. Crane. Probabilistic models for drug dissolution. Part 1. Review of Monte Carlo and stochastic cellular automata approaches. Simul. Model. Pract. Theory 14:843–856 (2006).
V. Grimm, E. Revilla, U. Berger, F. Jeltsch, W. M. Mooij, S. F. Railsback, H. H. Thulke, J. Weiner, T. Wiegand, and D. L. DeAngelis. Pattern-oriented modeling of agent-based complex systems: lessons from ecology. Science 310:987–991 (2005).
B. Faller. Enhancing knowledge through combination of in-silico and in-vitro ADME profiling assays, 1st annual phys. chem. symposium: early ADME and medical chemistry, Obernai, France, (2004). http://www.acdlabs.com/download/publ/2004/eum04/adme.pdf/. Cited 3/27/07.
R. Schulz and D. Winne. Relationship between antipyrine absorption and blood flow rate in rat jejunum, ileum, and colon. Naunyn–Schmiedeberg’s Arch. Pharmacol. 335:97–102 (1987).
N. Mork, H Bundgaard, M Shalmi, S Christensen. Furosemide prodrugs: synthesis, enzymatic hydrolysis and solubility of various furosemide esters. Int. J. Pharm. 60:163–169 (1990).
M. R. Grant, K. E. Mostov, T. D. Tlsty, and C. A. Hunt. Simulating properties of in vitro epithelial cell morphogenesis. PLoS Comput. Biol. 2:e129 (2006).
J. Tang, K. F. Ley, and C. A. Hunt. Dynamics of in silico leukocyte rolling, activation, and adhesion. BMC Syst. Biol. 1:14 (2007).
B. P. Zeigler, H. Praehofer, and T. G. Kim. Theory of modeling and simulation: integrating discrete event and continuous complex dynamic systems, Academic, San Diego, 2000.
Acknowledgment
This research was funded in part by grant R21CDH00101 from the CDH Research Foundation (of which CAH is a trustee). We thank Professors John Verboncoeur and George Sensabaugh for encouragement. We thank the members of Biosystems Group for helpful discussions, with special thanks going Teddy Lam, Sean Kim, Mark Grant, and Li Yan for commentary on the manuscript, to Sunwoo Park for technical support, and Glen Ropella for his support and keen technical and theoretical insights.
Author information
Authors and Affiliations
Corresponding author
Electronic supplementary material
Below is the linked to the electronic supplementary material:
Appendix
Appendix
Given the five expectations in the Introduction, we worked backwards and specified capabilities needed for their achievement, as done by Grant et al. (29), Tang et al. (30), and Hunt et al. (5) in other modeling contexts. It will require the design and validation of new types of in silico devices that can exhibit at least nine capabilities.
-
1.
Turing test: In silico “transport” data are, to a domain expert (in a type of Turing test), experimentally indistinguishable from the referent wet-lab transport data; achieving this capability requires the in silico device to be suitable for experimentation.
-
2.
Phenotype overlap: There is meaningful overlap between in silico and in vitro system phenotypes. Consequently, the events occurring during simulated transport, for example, accurately represent corresponding in vitro events at a level of detail (granularity) sufficient to mimic or predict specific phenotypic attributes. The concepts underlying designing analogues to achieve phenotype overlap (29,30) are similar to those on which Pattern-Oriented Modeling is based (31).
-
3.
Mappings: Observables in silico are designed to be consistent with in vitro observables. Doing so enables clear mappings between in vitro and in silico components and mechanisms.
-
4.
Transparency: Simulation details as they unfold need to be visualizable, measurable, and comparable where feasible to those of the in vitro system.
-
5.
Flexibility: It must be relatively simple to increase or decrease detail in order to simulate an additional phenotypic attribute or change usage and assumptions, without requiring significant reengineering. Achieving this capability will enable cycles of scientific testing and device refinement, as new wet-lab data becomes available. Having flexible methods will lead to improved, more realistic and heuristic analogue systems.
-
6.
Reusability: Device components can be designed to be autonomous and thus can be easily reconfigured to represent different cell types, experimental conditions, and compounds, alone or in combination. It must be easy to simulate and analyze outcomes of several different in silico experiments in a fraction of the time (and at a fraction of the cost) required to complete the wet-lab experiments.
-
7.
Adaptability: In addition to flexibility and reusability, the components must be constructed so that they can be adapted to function as components in physiologically based, intestine and whole organism models.
-
8.
Components articulate: Because the system attributes of interest may result from interactions of autonomous components, to study alternative mechanisms it must be easy to join, disconnect, and replace them without having to reengineer the in silico device or its components.
-
9.
Local mechanisms: as in vitro, the behaviors that emerge during a simulation are the consequences of local mechanisms–local component interactions.
To enable these capabilities, the in silico device must be synthetic, as is the in vitro model (4,5), and its framework must use discrete interactions. Temporally and spatially discrete models offer theoretical and practical advantages because they allow, even facilitate achieving the preceding capabilities. Moreover, in the case where a continuous model is desirable, a discrete model can simulate or approximate it to arbitrary precision (25).
The equations to determine D (including D a and D m) are:
The first equation is the Stokes–Einstein equation to calculate the diffusion coefficient for Brownian spheres. K B is the Boltzmans constant, ϖ is a remedy coefficient (ϖ = 1 for the cell membranes and 1.67 for aqueous spaces). η is the viscosity at 37°C, η = 5e − 3 Pa s for membranes and 0.69e − 3 Pa s for aqueous spaces). r is the radius of the referent molecule assuming it is an ideal sphere and V is the molecular van der Waals volume. V is proportional to MW with λ = 1.6755 according to the data in (12).
The equations to determine D eff in tight junctions are from (8) and (14):
D eff is the effective diffusion coefficient through the tight junctions. ɛ is the proportion of surface area covered by the tight junctions. D is the aqueous diffusion coefficient. F is the Renkin factor to describe the sieving effect of tight junctions on molecules of different sizes. τ is the tortuousity of the tight junction passage, which is estimated as 2.5–3.0 times the height of tight junctions. κ is the electrochemical energy. f 0 is the fraction of unionized molecules and \(f_{\kappa } \) is the fraction of charged molecules. R is the estimated average tight junction pore size. e is the unit charge of an electron. z is the valence of a molecule. Δψ is the electrical potential difference across the membrane. K is the thermal constant, and T is the Kelvin temperature converted from 37°C.
The mathematical equations derived from (8) to compute P total are:
M R is the mass in the acceptor space; M D(0) is the initial mass in the donor space; H R (2.6/4.71 cm) and H D (1.5/4.71 cm) are the thicknesses of the receiver and donor spaces, respectively; T is time.
The analytical solution to the substrate initial reaction rate ν in the presence of inhibitors for the two-site enzyme is based on the kinetic model in Fig. 11, and is (17):
V max is the maximal reaction rate of the substrate without inhibitors; K S is the Michaelis–Menten constant of the substrate; K i is the Michaelis–Menten constant of the inhibitor; α′ is the enhancing factor for the second site binding when one and only one site is bound with a molecule; \( \beta \prime {\left( {0 \leqslant \beta \prime \leqslant 1} \right)} \) is the inhibition factor of generating metabolites. The parameters in Table III used the equation above. Dissociation constants for all species in Fig. 11 are defined by \( K_{S} = {k_{{ - 1}} } \mathord{\left/ {\vphantom {{k_{{ - 1}} } {k_{1} }}} \right. \kern-\nulldelimiterspace} {k_{1} } \), \( K_{i} = {k_{{ - 2}} } \mathord{\left/ {\vphantom {{k_{{ - 2}} } {k_{2} }}} \right. \kern-\nulldelimiterspace} {k_{2} } \), \( \alpha \prime K_{i} = {k_{{ - 3}} } \mathord{\left/ {\vphantom {{k_{{ - 3}} } {k_{3} }}} \right. \kern-\nulldelimiterspace} {k_{3} } \) and \( \alpha \prime K_{\operatorname{S} } = {k_{{ - 4}} } \mathord{\left/ {\vphantom {{k_{{ - 4}} } {k_{4} }}} \right. \kern-\nulldelimiterspace} {k_{4} } \).
Rights and permissions
About this article
Cite this article
Garmire, L.X., Garmire, D.G. & Hunt, C.A. An In Silico Transwell Device for the Study of Drug Transport and Drug–Drug Interactions. Pharm Res 24, 2171–2186 (2007). https://doi.org/10.1007/s11095-007-9391-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11095-007-9391-4