Abstract
Computational pharmacokinetic (PK) modeling gives access to drug concentration vs. time profiles in target organs and allows better interpretation of clinical observations of therapeutic or toxic effects. Physiologically-based PK (PBPK) models in particular, based on mechanistic descriptions of the body anatomy and physiology, may also help to extrapolate in vitro or animal data to human.
Once in the systemic circulation, a chemical has access to the microvasculature of every organ or tissue. However, its penetration in the brain, retina, thymus, spinal cord, testis, placenta,… may be limited or even fully prevented by dynamic physiological blood-tissue barriers. Those barriers are both physical (involving tight junctions between adjacent cells) and biochemical (involving metabolizing enzymes and transporters).
On those cases, correct mechanistic characterization of the passage (or not) of molecules through the barrier can be crucial for improved PBPK modeling and prediction.
In parallel, attempts to understand and quantitatively characterize the processes involved in drug penetration of physiological barriers have led to the development of several in vitro experimental models. Data from such assays are very useful to calibrate PBPK models.
We review here those in vitro and computational models, highlighting the challenges and perspectives for in vitro and computational models to better assess drug availability to target tissues.
Similar content being viewed by others
Review
Computational tools permit to go beyond the frontiers of feasible experiments (McLanahan et al. 2012). They allow the generation and validation of mechanistic hypotheses and the number of results from simulations is virtually infinite. Pharmacokinetics (PK) models, and by extension physiologically-based PK (PBPK) models, aim at simulating quickly and at low cost the time course of the absorption, distribution, metabolism, and excretion (ADME) of a given drug in the body (Rostami-Hodjegan 2012). In doing so, they can provide predictions of a drug concentration in any organ or tissue (defined as a compartment) at any point in time. Access to drug concentration at the target cell level is important for understanding and predicting therapeutic or toxic effects; PBPK models are thus both important at the early stage of drug development and when assessing potential toxicity targets (Benjamin et al. 2010).
PBPK model equations and parameters characterize the various ADME processes and their interactions. Some of the absorption and distribution parameters may be estimated on the basis of tissue composition data and of the drug’s physicochemical properties via quantitative structure-properties models (Schmitt 2008). However, for estimating many of the PBPK parameters, in vitro experimental models have been developed and are essential (Cai et al. 2006).
There are many options when designing a PBPK model: the number of compartments is not limited, and many refinements are possible (Lee et al. 2009). Among them, the mechanistic description of passage through biological blood-tissue barriers appears very promising for both drug discovery and toxicity assessment. The challenges there are to characterize and predict permeation across the barriers, to design molecules which cross (or not) those barriers, and to have access to the effective chemical concentration in the target tissue.
In this review, we first recall the physiological basis for chemical distribution in tissues protected by biological barriers. We then describe the in vitro and computational tools to assess and predict barrier permeability. Finally, we provide an overview of challenges and perspectives in this area.
Chemical distribution and physiological barriers
ADME processes
After entering the body, a drug follows an ADME scheme (Willmann et al. 2005; Leahy 2003). Absorption corresponds to the process by which the compound enters the systemic circulation. This process crucially varies according to the administration route, dose, and form. For example, the rate-limiting step for absorption following oral administration may be either the dissolution rate (function of drug physicochemical characteristics and the physiological environment), or the transport rate (permeability) across the intestinal epithelium (Lennernas 2007). Distribution involves mechanisms of drug dispersion and transport throughout the fluids and tissues of the body. Distribution can be limited by either perfusion (when the tissues present no barrier to diffusion), or permeability across vascular/tissue barrier or across cell membranes inside tissues (Geldof et al. 2008). Metabolism deals with the biotransformation of parent drugs into metabolites, by metabolic enzymes such as cytochromes P450, dehydrogenases, transferases… (Emoto et al. 2010; Yengi et al. 2007). Finally, excretion is the removal of the drugs (or their metabolites) from the body (Aimone 2005).
According to these ADME processes, the free concentration of a drug in a specific tissue usually depends on its plasma free concentration, the plasma/tissue barrier permeability, its tissue binding, cellular membrane permeability, and metabolic modifications by cellular enzymes. Physiological barriers may be encountered at all absorption, distribution, and excretion steps (Kitamura et al. 2008), from the skin and the intestinal barriers regulating absorption, to distribution at the level of several target tissues like brain or testis, and to excretion in kidney, intestine,… An exhaustive review of all those biological barriers is out of the scope of this work and we will focus on blood-tissue barriers limiting distribution, for which details are given next.
Blood-tissue barriers
The key role of blood-tissue barriers is to modulate and restrain permeability (Alexis et al. 2008). They are both physical and biochemical. Physical barriers consist in a layer of cells with closely associated membranes between adjacent cells. Membrane occlusion is mediated by protein complexes, like occludins and claudins for the tight junctions, cadherins for adherens junctions, and connexins for gap junctions. Biochemical barriers involve the metabolism of chemicals by metabolizing enzymes like cytochromes P450, the influx or efflux of chemicals by carrier protein like ATP-Binding cassette transporters such as P-glycoprotein (P-gp). A complex and dynamic multi-pathways process is involved in passage of molecules across biological compartments (Figure 1).
The transcellular pathway corresponds to the passage of molecules through cells (Cheung and Brace 2008). Passive diffusion through cell membranes is the preferential route for small lipophilic molecules (< 400 Da) (Levin 1980). Passage can also be carrier-mediation via facilitated diffusion or an active (or secondary active) process (Lockman et al. 2008). Facilitated diffusion corresponds to the transport of a molecule following its concentration gradient. Active processes involve energy to move a molecule against its electrochemical gradient. This energy can be provided 1) directly by the transporter itself which is able to hydrolyze ATP, or 2) secondary following the use of an ion concentration gradient. The compounds using those specific transporter systems are mostly hydrophilic.
The paracellular pathway permits the passage of small hydrophilic molecules through gaps existing between cells, despites cell junctions (Vandenbroucke et al. 2008). This kind of transport is considered to be passive (driven by diffusion). However, charge selectivity often occurs, notably for modulating the passage of molecules via pores between junctional proteins. Tight junctions between cells are most of the time organized in complexes: tight and adherens junctions for brain (Stamatovic et al. 2008); tight, adherens, and gap junctions for testes (Wong et al. 2004). Tight junctions fuse membranes of two adjacent cells. Adherens junctions physically connect the cytosquelettons of neighboring cells. Gap junctions connect two adjacent cells and form pores allowing small molecules to pass between the cytoplasm of neighboring cells (Li et al. 2012).
Another process involved in limitating the passage of molecules is efflux transport, by active transporters such as P-gp, which extrude back molecules after they have entered the cells (Edwards et al. 2005).
The presence of metabolizing enzymes in the layer of cells forming the barrier may also modulate chemical availability at the tissue level (Miksys and Tyndale 2009). For example, chemicals can eventually be metabolized into a more active compound, a fact exploited by pro-drugs (Cucullo et al. 2012).
The presence of barriers between blood and tissues is a challenge for PBPK modeling, because of their relative complexity. Fortunately, in vitro and computational methods intending to predict permeability have been developed and can be put to use.
In vitro methods to evaluate chemical permeability through physiological barriers
In vitro models have been extensively used as alternative methodologies to in vivo models for evaluating the bioavailability at the level of organs protected by physiological barriers. Even if it is impossible to reproduce the complexity of in vivo systems, in vitro models can be designed with sufficient relevant features (Cucullo et al. 2011). An exhaustive review of all the available in vitro methods for assessing chemical permeability is out of the scope of this article (several reviews exist, see for example Sarmento et al. (2012)), and we chose to describe the main kinds of in vitro models, focusing on the endpoints measured and their relevance to PBPK model inputs.
Overview of the most commonly used in vitro models
According to the endpoint investigated, several in vitro models have been developed, ranging from the determination of partition coefficients to more mechanistic experiments (Table 1).
Cell-based assays consist of biological material from different levels of organization (from subcellular fractions to tissue fragments) and characterization (from well-described cellular transport processes to general information as used when determining partition coefficients). Tissue homogenates or slices can be used to estimate a drug’s partition coefficient between tissues and an extracellular medium (Friden et al. 2007). Cell membrane preparations are used to investigate particular pathways, transporters or receptors (Miller et al. 2011).
Some endothelial and epithelial cell cultures, when grown on permeable supports, spontaneously form monolayers and express functional junctions. Different levels of refinement are possible, from simple monolayers (Artursson et al. 2001; Trickler et al. 2010) to co-cultures (Antunes et al. 2013) to dynamic systems (Cucullo et al. 2011). These models allow dynamic permeability measurement through the monolayer. They also permit detailed characterization of biochemical mechanisms, such as receptor binding and uptake, or identification of relevant signaling pathways with mRNA and protein expression data (Seki et al. 2006).
Cell-based assays, however, are expensive and labor-intensive. In order to provide rapid, low-cost and automation-friendly tools to measure passive permeability, methods mimicking biological barriers with mixtures of lipids and organic solvents have been developed (Kv et al. 2008). For example, the immobilized artificial membranes (IAM) with HPLC (high-performance liquid chromatography) columns (Carrara et al. 2007) or the parallel artificial membrane permeability assay (PAMPA) give passive permeation estimates which correlate well with cell-based assays (Masungi et al. 2008).
Most of these tools can be automated for high throughput applications (Garberg et al. 1999). Because it is very difficult to develop a single in vitro system that can simulate the human in vivo setting, various in vitro assays are usually performed to investigate specific mechanisms (Abbott et al. 2008). Here, the key is to have a good knowledge of the barrier organization in vivo, both in terms of presence/absence of pathways and components, and in terms of their relative quantities.
Limitations of in vitro models and challenges for computational models
In vitro models have been optimized, becoming more and more complex in order to mimic as closely as possible biological barriers, and to permit investigation of several permeation mechanisms (Wuest et al. 2012; Hilgendorf et al. 2000).
Yet, in vitro systems remain quite simple and homogeneous, compared to the in vivo reality which involves several cell types, numerous processes, complex spatial structures, and much variability. It is difficult to obtain all the needed information in a given in vitro model, and computational tools allowing some representation of the complexity (i.e., membrane transport systems, metabolism pathways, cell polarity, and extracellular composition) would be very useful.
A universal issue with in vitro systems is the relevance of their results to in vivo settings, and the conditions for validity of their extrapolation to in vivo. Physiological barriers are most of time complexes of junction proteins and involve both physical and biochemical mechanisms. It is rarely possible to assess the effects of all this components on resulting permeability in a unique stand-alone in vitro experiment. Firstly, it has been shown that junctions are tighter in vivo than in vitro (Garberg et al. 2005); the in vitro configuration, although being improved with tridimensional structures, still lacks relevance. Also, barrier protein functionality is under strong dependence of intercellular signals (like in the blood-testis barrier (Siu et al. 2009)); the in vitro conditions may not allow information to be shared between cells, as well as functional proteins to be expressed (Schug et al. 2013). Finally, the link between any in vitro model and in vivo conditions has to be clearly defined and quantified, ideally with well-justified mathematical descriptions. Spatial organization may also be illustrated computationally.
Computational methods to describe and predict drug permeability through physiological barriers
PBPK models describe specific organs, tissues, and subcellular localizations as a set of pre-defined compartments linked together by the vascular system. In all those models, drug transport occurs at least via blood (Pilari and Huisinga 2010). Formally, they correspond to systems of differential equations for the concentrations or quantities of a given drug in each compartment. These compartments can be general (blood, “highly-perfused tissues”, “poorly perfused tissues”…) or very well described (Graf et al. 2012). In flow-limited models, the derivative for the quantity X T of a drug X in tissue T is typically defined as:
where QT is the blood flow rate, Cart the arterial blood concentration, V T the tissue volume, and P T the tissue-to-blood partition coefficient. The partition coefficient is still a meaningful parameter for an organ protected by a passive barrier, even if the model is quite simplistic. The rate of entry in the tissue will be over-estimated, however, thereby over-estimating tissue exposure.
A better approach for describing a permeability-limited transport (either transcellular, paracellular, or both) is to sub-compartmentalize the PBPK model: organs get usually sub-divided into three (vascular, extracellular and intracellular) compartments (Campbell 2009). Exchange rates between those compartments then fall back on linear processes, such as Fick’s law of diffusion (Kramer et al. 2009), or saturable transport if need is. A simpler sub-division into two sub-compartments will suffice to illustrate the model equations and introduce the needed parameters. The derivative for the quantity X TB in the tissue vascular blood compartment of tissue T can be calculated as:
where, in addition to the parameters and variables defined in eq. 1, V TB is the vascular blood volume in the tissue, PS T the apparent permeability-surface area product, f ub the unbound fraction of the compound in blood, f ut the unbound fraction of the compound in tissue, X TI the quantity of X in the interstitial and intracellular compartment and V TI the volume of that compartment (V TI = V T – V TB ).
For X TI , the corresponding differential equation is:
To model active transport, a transport rate term J (nmol/min/mg protein) can be added to equation 3 and subtracted from equation 2 (for influx, the reverse if efflux is considered). J is usually based on conventional Michaelis-Menten kinetics:
where V max is the maximum transport capacity (nmol/min/mg protein), and K m the half-saturation concentration. For efflux, the TB subscripts would be replaced by TI.
This leaves us with specific parameters to estimate. Such estimates can be obtained in vitro. Partition coefficients can be obtained either by model-based predictions according to the tissue composition (where the equation parameters are given by in vitro experiments on drug lipophilicity and plasma protein binding), or by in vitro direct determination of the drug's concentrations ratio in the buffer and tissue at steady-state (Poulin and Theil 2002b2002a; Rodgers et al. 2005; Rodgers and Rowland 2006). In cell monolayers or membranes in vitro, permeability can be evaluated using the equations for chemical flux, based on Fick’s first law for passive diffusion (Kramer et al. 2009). In that case, in vitro experiments allow a direct determination of an apparent permeability coefficient P app (cm/s):
where dQ/dt (mol/s) is the increase in the amount of drug/tracer molecule in the receiver compartment after a small time interval dt since time zero, S (cm2) the exchange surface, and C 0 (mol/cm3) the initial drug/tracer molecule concentration in the donor chamber.
The apparent permeability-surface area product, PS T , is simply (Koda et al. 2007):
Since junctions between adjacent cells restrict paracellular passage, a correction factor can be used when extrapolating from one barrier condition to another, for example when tight junctions are affected (Kondoh et al. 2012). Experimental permeability measurements may also correspond to the sum of transcellular, paracellular, or carrier-mediated passages (Amasheh et al. 2009). A specific permeability rate for paracellular passage can be assessed by using a tracer molecule which does not penetrate in cells (e.g., lucifer yellow) (Inokuchi et al. 2009).
For active transport, determination of V max and K m follows simply the lines of Michaelis-Menten parameters estimation (Weiss and Kang 2002). When both active and passive transport are present, it is possible to fit a mixed model to the in vitro concentration time course data, at early times:
The above parameters can also be estimated by fitting the whole PBPK model to in vivo PK data, but that is obviously not feasible in the absence of such data, as when in vitro to in vivo extrapolation is sought. The third option is to obtain those parameter estimates from computational "sub-models". Two main kinds of sub-models are used for that purpose: statistical models based on quantitative structure property relationship (QSPR), and mathematical models based on mechanistic description of biological processes.
Quantitative structure property relationship (QSPR)-based permeability models
One way to estimate the permeability of drugs through a given barrier is to link their chemical structures to the property of interest: permeability. The principle of QSPR models is that similar chemical structures should lead to similar properties (Sheridan et al. 2009). These are typical empirical statistical models, calibrated on a training dataset of chemicals of known structures and properties. Structures have to be somehow quantified through “descriptors” which enter the model as input variables, and the model gives a value for the property of interest (Neely et al. 2009). For a new drug, the calibrated model can therefore predict permeability, for example, simply on the basis of its chemical structure (provided that the new drug structure is not too far away from the structure of the drugs used in the training set).
The ability of drugs to traverse a tissue barrier is conditioned by tissue blood flow and several physicochemical properties: molecular size, lipophilicity, plasma protein binding, efflux pump affinity, molecular charge (Giaginis et al. 2009). QSPR modeling efforts for barrier-crossing have notably concentrated on passive permeability, for which enough data was available. As a further simplification, they concentrate on the prediction of partition coefficients (Chuman 2008), and mostly use molecular size and lipophilicity as descriptors. The impact of molecular size on paracellular permeability makes sense intuitively: the larger the molecule, the lower its ability to diffuse through the tight junctions of the tissue barrier. Lipophilicity property, often described by the octanol/water partition, Log P, influences the transcellular passage through lipid membranes.
Because the data necessary for predicting active or facilitated transport processes is currently insufficient, only a few QSPR models exist for those transport processes (Friedrichsen et al. 2001). Attempts to develop QSPR models for P-gp efflux have highlighted difficulties due to the broad specificity of this transporter (Gombar et al. 2004). Some QSPR models have also been developed for predicting enzymatic reactions like, for example, metabolic inactivation (Ekuase et al. 2011).
In spite of their usefulness for high-throughput screening of compounds on the basis of their physicochemical properties, QSPR models may lack predictive capacity (Chen et al. 2007). Mechanism-based models, which should lead to better predictions and on a time and space scale which would not be limited to the available information, are useful complements.
Mechanism-based permeability models
Another way to predict permeability is to describe its mechanisms according to biochemical and physical laws. A finer description of mechanisms, taking receptor uptake and transfer through membrane by carriers into account, can bring high value to PBPK models, in particular in the high dose range used in therapeutic applications (Tanaka et al. 1999).
There have been several efforts to develop computational ways to gain mechanistic insight in permeability processes. For example, Garmire et al. (2007) developed an in silico transwell device, incorporating both spatial (tight junctions) and functional (metabolizing enzymes) barrier features, for mimicking the in vitro passive transport properties through cell monolayers. Model predictions of passive transport were validated across in vitro data, and were reasonable approximations.
Another computational development is that of Dolghih and Jacobson (2013). They developed a computational approach of the blood–brain barrier combining two mechanism-based models: for passive permeation and for active efflux by P-gp. This model is a good illustration of the importance to combine permeability mechanisms for obtaining meaningful predictions.
Perspectives for quantitative predictions of bioavailability to barriers-protected organs
Conceptualization
A well thought conceptualization is the key for developing a useful model. We have seen that there are several ways to build a PBPK model, from description of few global compartments to precise compartmentalization at the tissue or even at the cell levels. The same is true for parameters which can be estimated at whole tissue level or at a finer level.
One important thing in PBPK modeling is to have an integrated view of the whole body behavior. To model every detail is not feasible, and abstraction is essential to understand. Hence, an approach coupling general chemical behavior in the whole body and refinement through description at lower levels in tissues of special interest can give very useful insight about relevant mechanisms, while being predictive at a higher level.
The success of QSPR-based models depends both on the accuracy of algorithms and on the quality of input data. Because of the methodological and experimental variability present in published data, QSPR models may be poorly predictive. This highlights the need of a balance between i) the understanding of the reality of complex physiologic events to provide accurate systems information and ii) the simplicity required for computational modeling feasibility.
Different tools detailing each step of ADME processes are available and can be integrated into a coherent iterative approach (Honorio et al. 2012; Bois et al. 2010; Jamei et al. 2009). Yet again, the experimental evidence is the limit in developing mechanistic models. Thus, when designing a model, the balance has to be made between i) the available information and ii) the model complexity to meet the objectives.
Refinement
The increasing understanding of mechanisms intervening in chemical distribution allows the refinement of predictions about a chemical availability at the site of its effect, be it therapeutic or toxic. Another scale of complexity for chemical distribution may be to work at the level of cell organelles (Zheng et al. 2011).
As the area of systems biology is growing, investigating barrier formation mechanisms via signaling pathways can be of interest to predict behavior at a higher scale and earlier (Lee et al. 2010). Furthermore, besides the intrinsic functionality of physiological barriers, systems biology can also describe dynamical/homeostasis barrier mechanisms (Smallwood 2009).
Validation
As usual with computational models, usefulness and relevance have to be carefully assessed. Indeed, in the case of permeability models as in general, models of biological systems are strong on assumptions and weak on validation. A general framework for validation is provided by Sornette et al. (Sornette et al. 2007). They propose a formal iterative process, leading to the model rejection or progressive refinement and validation. In the particular case of computational models of biological systems, robust validation techniques against biological models have to be employed. Smallwood et al. (Smallwood et al. 2004) describe a modeling paradigm for developing a relevant predictive computational model of cellular interaction. In their example, the key is to understand the in vitro behavior. Cellular processes are stochastic, and the generation of distributions using Monte Carlo techniques appears to be the most relevant. Of course, the ideal would be to have a set of in vitro experiments and to apply a comparison metrics.
Conclusions
Time-course of a drug concentration in organs protected by biological barriers can be obtained thanks to i) the time-course of concentration in blood and ii) the characterization of passage through those barriers. The first kind of information can be obtained from whole-body PBPK models, the second may come from several methods, in vitro experiments, QSPR models, or mechanism-based models defining a set of precise and highly specific parameters.
Both understanding of mechanisms and estimating model parameters in the field of barrier permeability can be done thanks to in vitro models. Of course, the complexity of the living systems and the imperfect predictive power of computational models mean that in vitro assays will be still used for a long time.
As long as experimental evidence come and scientific hypotheses are validated, the mathematical models for predicting passage across biological barriers will tend to be increasingly complex and realistic. They should thus lead to better extrapolation, from in vitro to in vivo, or from animal to human.
References
Abbott NJ, Dolman DE, Patabendige AK: Assays to predict drug permeation across the blood–brain barrier, and distribution to brain. Curr Drug Metab 2008,9(9):901–910. 10.2174/138920008786485182
Aimone LD: Overview of pharmacokinetics. Curr Protoc Pharmacol 2005, 7: Unit7.
Alexis F, Pridgen E, Molnar LK, Farokhzad OC: Factors affecting the clearance and biodistribution of polymeric nanoparticles. Mol Pharm 2008,5(4):505–515. 10.1021/mp800051m
Amasheh M, Grotjohann I, Amasheh S, Fromm A, Soderholm JD, Zeitz M, Fromm M, Schulzke JD: Regulation of mucosal structure and barrier function in rat colon exposed to tumor necrosis factor alpha and interferon gamma in vitro: a novel model for studying the pathomechanisms of inflammatory bowel disease cytokines. Scand J Gastroenterol 2009,44(10):1226–1235. 10.1080/00365520903131973
Antunes F, Andrade F, Ferreira D, Nielsen HM, Sarmento B: Models to Predict Intestinal Absorption of Therapeutic Peptides and Proteins. Curr Drug Metab 2013,14(1):4–20. 10.2174/138920013804545160
Artursson P, Palm K, Luthman K: Caco-2 monolayers in experimental and theoretical predictions of drug transport. Adv Drug Deliv Rev 2001,46(1–3):27–43.
Benjamin B, Barman TK, Chaira T, Paliwal JK: Integration of physicochemical and pharmacokinetic parameters in lead optimization: a physiological pharmacokinetic model based approach. Curr Drug Discov Tech 2010,7(3):143–153.
Bois FY, Jamei M, Clewell HJ: PBPK modelling of inter-individual variability in the pharmacokinetics of environmental chemicals. Toxicology 2010,278(3):256–267. 10.1016/j.tox.2010.06.007
Cai H, Stoner C, Reddy A, Freiwald S, Smith D, Winters R, Stankovic C, Surendran N: Evaluation of an integrated in vitro-in silico PBPK (physiologically based pharmacokinetic) model to provide estimates of human bioavailability. Int J Pharm 2006,308(1–2):133–139.
Campbell A: Development of PBPK model of molinate and molinate sulfoxide in rats and humans. Regul Toxicol Pharmacol 2009,53(3):195–204. 10.1016/j.yrtph.2009.01.003
Carrara S, Reali V, Misiano P, Dondio G, Bigogno C: Evaluation of in vitro brain penetration: optimized PAMPA and MDCKII-MDR1 assay comparison. Int J Pharm 2007,345(1–2):125–133.
Chen CY, Ko CW, Lee PI: Toxicity of substituted anilines to Pseudokirchneriella subcapitata and quantitative structure-activity relationship analysis for polar narcotics. Environ Toxicol Chem /SETAC 2007,26(6):1158–1164. 10.1897/06-293R.1
Cheung CY, Brace RA: Unidirectional transport across cultured ovine amniotic epithelial cell monolayer. Reprod Sci (Thousand Oaks, Calif 2008,15(10):1054–1058. 10.1177/1933719108322426
Chuman H: Toward basic understanding of the partition coefficient log P and its application in QSAR. SAR QSAR Environ Res 2008,19(1–2):71–79.
Cucullo L, Marchi N, Hossain M, Janigro D: A dynamic in vitro BBB model for the study of immune cell trafficking into the central nervous system. J Cereb Blood Flow Metab 2011,31(2):767–777. 10.1038/jcbfm.2010.162
Dolghih E, Jacobson MP: Predicting efflux ratios and blood–brain barrier penetration from chemical structure: combining passive permeability with active efflux by p-glycoprotein. ACS Chem Neurosci 2013,4(2):361–367. 10.1021/cn3001922
Edwards JE, Alcorn J, Savolainen J, Anderson BD, McNamara PJ: Role of P-glycoprotein in distribution of nelfinavir across the blood-mammary tissue barrier and blood–brain barrier. Antimicrob Agents Chemother 2005,49(4):1626–1628. 10.1128/AAC.49.4.1626-1628.2005
Ekuase EJ, Liu Y, Lehmler HJ, Robertson LW, Duffel MW: Structure-activity relationships for hydroxylated polychlorinated biphenyls as inhibitors of the sulfation of dehydroepiandrosterone catalyzed by human hydroxysteroid sulfotransferase SULT2A1. Chem Res Toxicol 2011,24(10):1720–1728. 10.1021/tx200260h
Emoto C, Murayama N, Rostami-Hodjegan A, Yamazaki H: Methodologies for investigating drug metabolism at the early drug discovery stage: prediction of hepatic drug clearance and P450 contribution. Curr Drug Metab 2010,11(8):678–685. 10.2174/138920010794233503
Friden M, Gupta A, Antonsson M, Bredberg U, Hammarlund-Udenaes M: In vitro methods for estimating unbound drug concentrations in the brain interstitial and intracellular fluids. Drug Metabol Dispos: the biological fate of chemicals 2007,35(9):1711–1719. 10.1124/dmd.107.015222
Friedrichsen GM, Jakobsen P, Taub M, Begtrup M: Application of enzymatically stable dipeptides for enhancement of intestinal permeability. Synthesis and in vitro evaluation of dipeptide-coupled compounds. Bioorg Med Chem 2001,9(10):2625–2632. 10.1016/S0968-0896(01)00066-9
Garberg P, Eriksson P, Schipper N, Sjostrom B: Automated absorption assessment using Caco-2 cells cultured on both sides of polycarbonate membranes. Pharm Res 1999,16(3):441–445. 10.1023/A:1018838121975
Garberg P, Ball M, Borg N, Cecchelli R, Fenart L, Hurst RD, Lindmark T, Mabondzo A, Nilsson JE, Raub TJ, Stanimirovic D, Terasaki T, Oberg JO, Osterberg T: In vitro models for the blood–brain barrier. Toxicol In Vitro 2005,19(3):299–334. 10.1016/j.tiv.2004.06.011
Garmire LX, Garmire DG, Hunt CA: An in silico transwell device for the study of drug transport and drug-drug interactions. Pharm Res 2007,24(12):2171–2186. 10.1007/s11095-007-9391-4
Geldof M, Freijer J, van Beijsterveldt L, Danhof M: Pharmacokinetic modeling of non-linear brain distribution of fluvoxamine in the rat. Pharm Res 2008,25(4):792–804. 10.1007/s11095-007-9390-5
Giaginis C, Zira A, Theocharis S, Tsantili-Kakoulidou A: Application of quantitative structure-activity relationships for modeling drug and chemical transport across the human placenta barrier: a multivariate data analysis approach. J Appl Toxicol 2009,29(8):724–733. 10.1002/jat.1466
Gombar VK, Polli JW, Humphreys JE, Wring SA, Serabjit-Singh CS: Predicting P-glycoprotein substrates by a quantitative structure-activity relationship model. J Pharm Sci 2004,93(4):957–968. 10.1002/jps.20035
Graf JF, Scholz BJ, Zavodszky MI: BioDMET: a physiologically based pharmacokinetic simulation tool for assessing proposed solutions to complex biological problems. J Pharmacokinet Pharmacodyn 2012,39(1):37–54. 10.1007/s10928-011-9229-x
Hilgendorf C, Spahn-Langguth H, Regardh CG, Lipka E, Amidon GL, Langguth P: Caco-2 versus Caco-2/HT29-MTX co-cultured cell lines: permeabilities via diffusion, inside- and outside-directed carrier-mediated transport. J Pharm Sci 2000,89(1):63–75. 10.1002/(SICI)1520-6017(200001)89:1<63::AID-JPS7>3.0.CO;2-6
Honorio KM, Moda TL, Andricopulo AD: Pharmacokinetic Properties and In Silico Adme Modeling In Drug Discovery. Medicinal chemistry, Shariqah (United Arab Emirates); 2012.
Inokuchi H, Takei T, Aikawa K, Shimizu M: The effect of hyperosmosis on paracellular permeability in Caco-2 cell monolayers. Biosci Biotechnol Biochem 2009,73(2):328–334. 10.1271/bbb.80538
Jamei M, Marciniak S, Feng K, Barnett A, Tucker G, Rostami-Hodjegan A: The Simcyp population-based ADME simulator. Expert Opin Drug Metab Toxicol 2009,5(2):211–223. 10.1517/17425250802691074
Kitamura S, Maeda K, Sugiyama Y: Recent progresses in the experimental methods and evaluation strategies of transporter functions for the prediction of the pharmacokinetics in humans. Naunyn Schmiedebergs Arch Pharmacol 2008,377(4–6):617–628.
Koda Y, Shiotani K, Toth I, Tsuda Y, Okada Y, Blanchfield JT: Comparison of the in vitro apparent permeability and stability of opioid mimetic compounds with that of the native peptide. Bioorg Med Chem Lett 2007,17(7):2043–2046. 10.1016/j.bmcl.2007.01.051
Kondoh M, Takahashi A, Yagi K: Spiral progression in the development of absorption enhancers based on the biology of tight junctions. Adv Drug Deliv Rev 2012,64(6):515–522. 10.1016/j.addr.2011.07.004
Kramer SD, Lombardi D, Primorac A, Thomae AV, Wunderli-Allenspach H: Lipid-bilayer permeation of drug-like compounds. Chem Biodivers 2009,6(11):1900–1916. 10.1002/cbdv.200900122
Kv S, Devi GS, Mathew ST: Liposomal formulations of serratiopeptidase: in vitro studies using PAMPA and Caco-2 models. Mol Pharm 2008,5(1):92–97. 10.1021/mp700090r
Leahy DE: Progress in simulation modelling for pharmacokinetics. Curr Top Med Chem 2003,3(11):1257–1268. 10.2174/1568026033451961
Lee SK, Hamer D, Bedwell CL, Lohitnavy M, Yang RS: Effect of PCBs on the lactational transfer of methyl mercury in mice: PBPK modeling. Environ Toxicol Pharmacol 2009,27(1):75–83. 10.1016/j.etap.2008.08.014
Lee SE, Choi KJ, Menon GK, Kim HJ, Choi EH, Ahn SK, Lee SH: Penetration pathways induced by low-frequency sonophoresis with physical and chemical enhancers: iron oxide nanoparticles versus lanthanum nitrates. J Investig Dermatol 2010,130(4):1063–1072. 10.1038/jid.2009.361
Lennernas H: Modeling gastrointestinal drug absorption requires more in vivo biopharmaceutical data: experience from in vivo dissolution and permeability studies in humans. Curr Drug Metab 2007,8(7):645–657. 10.2174/138920007782109823
Levin VA: Relationship of octanol/water partition coefficient and molecular weight to rat brain capillary permeability. J Med Chem 1980,23(6):682–684. 10.1021/jm00180a022
Li MW, Mruk DD, Cheng CY: Gap junctions and blood-tissue barriers. Adv Exp Med Biol 2012, 763: 260–280.
Lockman PR, Manda VK, Geldenhuys WJ, Mittapalli RK, Thomas F, Albayati ZF, Crooks PA, Dwoskin LP, Allen DD: Carrier-mediated transport of the quaternary ammonium neuronal nicotinic receptor antagonist n, n'-dodecylbispicolinium dibromide at the blood–brain barrier. J Pharmacol Exp Ther 2008,324(1):244–250.
Masungi C, Mensch J, Van Dijck A, Borremans C, Willems B, Mackie C, Noppe M, Brewster ME: Parallel artificial membrane permeability assay (PAMPA) combined with a 10-day multiscreen Caco-2 cell culture as a tool for assessing new drug candidates. Die Pharmazie 2008,63(3):194–199.
McLanahan ED, El-Masri HA, Sweeney LM, Kopylev LY, Clewell HJ, Wambaugh JF, Schlosser PM: Physiologically based pharmacokinetic model use in risk assessment–Why being published is not enough. Toxicol Sci 2012,126(1):5–15. 10.1093/toxsci/kfr295
Miksys S, Tyndale RF: Brain drug-metabolizing cytochrome P450 enzymes are active in vivo, demonstrated by mechanism-based enzyme inhibition. Neuropsychopharmacology 2009,34(3):634–640. 10.1038/npp.2008.110
Miller DW, Hinton M, Chen F: Evaluation of drug efflux transporter liabilities of darifenacin in cell culture models of the blood–brain and blood-ocular barriers. Neurourol Urodyn 2011,30(8):1633–1638. 10.1002/nau.21110
Neely BJ, Madihally SV, Robinson RL Jr, Gasem KA: Nonlinear quantitative structure–property relationship modeling of skin permeation coefficient. J Pharm Sci 2009,98(11):4069–4084. 10.1002/jps.21678
Pilari S, Huisinga W: Lumping of physiologically-based pharmacokinetic models and a mechanistic derivation of classical compartmental models. J Pharmacokinet Pharmacodyn 2010,37(4):365–405. 10.1007/s10928-010-9165-1
Poulin P, Theil FP: Prediction of pharmacokinetics prior to in vivo studies. 1. Mechanism-based prediction of volume of distribution. J Pharm Sci 2002,91(1):129–156. 10.1002/jps.10005
Poulin P, Theil FP: Prediction of pharmacokinetics prior to in vivo studiesII. Generic physiologically based pharmacokinetic models of drug disposition. J Pharmaceut Sci 2002,91(5):1358–1370. 10.1002/jps.10128
Rodgers T, Rowland M: Physiologically based pharmacokinetic modelling 2: predicting the tissue distribution of acids, very weak bases, neutrals and zwitterions. J Pharm Sci 2006,95(6):1238–1257. 10.1002/jps.20502
Rodgers T, Leahy D, Rowland M: Physiologically based pharmacokinetic modeling 1: predicting the tissue distribution of moderate-to-strong bases. J Pharm Sci 2005,94(6):1259–1276. 10.1002/jps.20322
Rostami-Hodjegan A: Physiologically based pharmacokinetics joined with in vitro-in vivo extrapolation of ADME: a marriage under the arch of systems pharmacology. Clin Pharmacol Ther 2012,92(1):50–61. 10.1038/clpt.2012.65
Sarmento B, Andrade F, da Silva SB, Rodrigues F, Das Neves J, Ferreira D: Cell-based in vitro models for predicting drug permeability. Expert Opin Drug Metab Toxicol 2012,8(5):607–621. 10.1517/17425255.2012.673586
Schmitt W: General approach for the calculation of tissue to plasma partition coefficients. Toxicol In Vitro 2008,22(2):457–467. 10.1016/j.tiv.2007.09.010
Schug M, Stober R, Heise T, Mielke H, Gundert-Remy U, Godoy P, Reif R, Blaszkewicz M, Ellinger-Ziegelbauer H, Ahr HJ, Selinski S, Gunther G, Marchan R, Blaszkewicz M, Sachinidis A, Nussler A, Oberemm A, Hengstler JG: Pharmacokinetics explain in vivo/in vitro discrepancies of carcinogen-induced gene expression alterations in rat liver and cultivated hepatocytes. Arch Toxicol 2013,87(2):337–345. 10.1007/s00204-012-0999-8
Seki T, Kanbayashi H, Nagao T, Chono S, Tabata Y, Morimoto K: Effect of cationized gelatins on the paracellular transport of drugs through caco-2 cell monolayers. J Pharm Sci 2006,95(6):1393–1401. 10.1002/jps.20616
Sheridan RP, Nam K, Maiorov VN, McMasters DR, Cornell WD: QSAR models for predicting the similarity in binding profiles for pairs of protein kinases and the variation of models between experimental data sets. J Chem Inf Model 2009,49(8):1974–1985. 10.1021/ci900176y
Siu ER, Wong EW, Mruk DD, Porto CS, Cheng CY: Focal adhesion kinase is a blood-testis barrier regulator. Proc Natl Acad Sci USA 2009,106(23):9298–9303. 10.1073/pnas.0813113106
Smallwood R: Computational modeling of epithelial tissues. Wiley Interdiscipl Rev 2009,1(2):191–201.
Smallwood RH, Holcombe WM, Walker DC: Development and validation of computational models of cellular interaction. J Mol Histol 2004,35(7):659–665. 10.1007/s10735-004-2660-1
Sornette D, Davis AB, Ide K, Vixie KR, Pisarenko V, Kamm JR: Algorithm for model validation: theory and applications. Proc Natl Acad Sci USA 2007,104(16):6562–6567. 10.1073/pnas.0611677104
Stamatovic SM, Keep RF, Andjelkovic AV: Brain endothelial cell-cell junctions: how to "open" the blood brain barrier. Curr Neuropharmacol 2008,6(3):179–192. 10.2174/157015908785777210
Tanaka C, Kawai R, Rowland M: Physiologically based pharmacokinetics of cyclosporine A: reevaluation of dose-nonlinear kinetics in rats. J Pharmacokinet Biopharm 1999,27(6):597–623. 10.1023/A:1020978509566
Trickler WJ, Lantz SM, Murdock RC, Schrand AM, Robinson BL, Newport GD, Schlager JJ, Oldenburg SJ, Paule MG, Slikker W Jr, Hussain SM, Ali SF: Silver nanoparticle induced blood–brain barrier inflammation and increased permeability in primary rat brain microvessel endothelial cells. Toxicol Sci 2010,118(1):160–170. 10.1093/toxsci/kfq244
Vandenbroucke E, Mehta D, Minshall R, Malik AB: Regulation of endothelial junctional permeability. Ann N Y Acad Sci 2008, 1123: 134–145. 10.1196/annals.1420.016
Weiss M, Kang W: P-glycoprotein inhibitors enhance saturable uptake of idarubicin in rat heart: pharmacokinetic/pharmacodynamic modeling. J Pharmacol Exp Ther 2002,300(2):688–694. 10.1124/jpet.300.2.688
Willmann S, Lippert J, Schmitt W: From physicochemistry to absorption and distribution: predictive mechanistic modelling and computational tools. Expert Opin Drug Metab Toxicol 2005,1(1):159–168. 10.1517/17425255.1.1.159
Wong CH, Mruk DD, Lui WY, Cheng CY: Regulation of blood-testis barrier dynamics: an in vivo study. J Cell Sci 2004,117(Pt 5):783–798.
Wuest DM, Wing AM, Lee KH: Membrane configuration optimization for a murine in vitro blood–brain barrier model. J Neurosci Methods 2012,212(2):211–221.
Yang Y, Voak A, Wilkinson SR, Hu L: Design, synthesis, and evaluation of potential prodrugs of DFMO for reductive activation. Bioorg Med Chem Lett 2012,22(21):6583–6586. 10.1016/j.bmcl.2012.09.005
Yengi LG, Leung L, Kao J: The evolving role of drug metabolism in drug discovery and development. Pharm Res 2007,24(5):842–858. 10.1007/s11095-006-9217-9
Zheng N, Tsai HN, Zhang X, Rosania GR: The subcellular distribution of small molecules: from pharmacokinetics to synthetic biology. Mol Pharm 2011,8(5):1619–1628. 10.1021/mp200092v
Acknowledgements
The author is grateful to Frédéric Bois and the reviewers for critical reading of the manuscript.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The author declares that she has no competing interest.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Quignot, N. Modeling bioavailability to organs protected by biological barriers. In Silico Pharmacol. 1, 8 (2013). https://doi.org/10.1186/2193-9616-1-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/2193-9616-1-8