Physiologically based pharmacokinetic (PBPK) modeling is a
mathematical modeling
A mathematical model is a description of a system using mathematical concepts and language. The process of developing a mathematical model is termed mathematical modeling. Mathematical models are used in the natural sciences (such as physics, b ...
technique for predicting the absorption, distribution, metabolism and excretion (ADME) of synthetic or natural chemical substances in humans and other animal species. PBPK modeling is used in pharmaceutical research and drug development, and in health
risk assessment
Broadly speaking, a risk assessment is the combined effort of:
# identifying and analyzing potential (future) events that may negatively impact individuals, assets, and/or the environment (i.e. hazard analysis); and
# making judgments "on the to ...
for cosmetics or general chemicals.
PBPK models strive to be mechanistic by mathematically transcribing anatomical, physiological, physical, and chemical descriptions of the phenomena involved in the complex ADME processes. A large degree of residual simplification and empiricism is still present in those models, but they have an extended domain of applicability compared to that of classical, empirical function based,
pharmacokinetic
Pharmacokinetics (from Ancient Greek ''pharmakon'' "drug" and ''kinetikos'' "moving, putting in motion"; see chemical kinetics), sometimes abbreviated as PK, is a branch of pharmacology dedicated to determining the fate of substances administered ...
models. PBPK models may have purely predictive uses, but other uses, such as statistical inference, have been made possible by the development of
Bayesian
Thomas Bayes (/beɪz/; c. 1701 – 1761) was an English statistician, philosopher, and Presbyterian minister.
Bayesian () refers either to a range of concepts and approaches that relate to statistical methods based on Bayes' theorem, or a followe ...
statistical tools able to deal with complex models. That is true for both toxicity risk assessment and therapeutic drug development.
PBPK models try to rely ''a priori'' on the anatomical and physiological structure of the body, and to a certain extent, on biochemistry. They are usually
multi-compartment model
A multi-compartment model is a type of mathematical model used for describing the way materials or energies are transmitted among the ''compartments'' of a system. Sometimes, the physical system that we try to model in equations is too complex, so ...
s, with compartments corresponding to predefined organs or tissues, with interconnections corresponding to blood or lymph flows (more rarely to diffusions). A system of
differential equations
In mathematics, a differential equation is an equation that relates one or more unknown functions and their derivatives. In applications, the functions generally represent physical quantities, the derivatives represent their rates of change, an ...
for concentration or quantity of substance on each compartment can be written, and its parameters represent blood flows, pulmonary ventilation rate, organ volumes etc., for which information is available in scientific publications. Indeed, the description they make of the body is simplified and a balance needs to be struck between complexity and simplicity. Besides the advantage of allowing the recruitment of ''a priori'' information about parameter values, these models also facilitate inter-species transpositions or extrapolation from one mode of administration to another (''e.g.'', inhalation to oral). An example of a 7-compartment PBPK model, suitable to describe the fate of many solvents in the mammalian body, is given in the Figure on the right.
History
The first pharmacokinetic model described in the scientific literature
was in fact a PBPK model. It led, however, to computations intractable at that time. The focus shifted then to simpler models
,
for which analytical solutions could be obtained (such solutions were sums of exponential terms, which led to further simplifications.) The availability of computers and numerical integration algorithms marked a renewed interest in
physiological models in the early 1970s.
For substances with complex kinetics, or when inter-species extrapolations were required, simple models were insufficient and research continued on physiological models
.
By 2010, hundreds of scientific publications have described and used PBPK models, and at least two private companies are basing their business on their expertise in this area.
Building a PBPK model
The model equations follow the principles of mass transport, fluid dynamics, and biochemistry in order to simulate the fate of a substance in the body
.
Compartments are usually defined by grouping organs or tissues with similar blood perfusion rate and lipid content (''i.e.'' organs for which chemicals' concentration ''vs.'' time profiles will be similar). Ports of entry (lung, skin, intestinal tract...), ports of exit (kidney, liver...) and target organs for therapeutic effect or toxicity are often left separate. Bone can be excluded from the model if the substance of interest does not distribute to it. Connections between compartment follow physiology (''e.g.'', blood flow in exit of the gut goes to liver, ''etc.'')
Basic transport equations
Drug distribution into a tissue can be rate-limited by either perfusion or permeability. Perfusion-rate-limited kinetics apply when the tissue membranes present no barrier to diffusion. Blood flow, assuming that the drug is transported mainly by blood, as is often the case, is then the limiting factor to distribution in the various cells of the body. That is usually true for small lipophilic drugs. Under perfusion limitation, the instantaneous rate of entry for the quantity of drug in a compartment is simply equal to (blood) volumetric flow rate through the organ times the incoming blood concentration. In that case; for a generic compartment ''i'', the differential equation for the quantity ''Q
i'' of substance, which defines the rate of change in this quantity, is:
where ''F
i'' is blood flow (noted ''Q'' in the Figure above), ''C
art'' incoming arterial blood concentration, ''P
i'' the tissue over blood partition coefficient and ''V
i'' the volume of compartment ''i''.
A complete set of differential equations for the 7-compartment model shown above could therefore be given by the following table:
The above equations include only transport terms and do not account for inputs or outputs.
Those can be modeled with specific terms, as in the following.
Modeling inputs
Modeling inputs is necessary to come up with a meaningful description of a chemical's pharmacokinetics. The following examples show how to write the corresponding equations.
Ingestion
When dealing with an oral
bolus dose (''e.g.'' ingestion of a tablet), first order absorption is a very common assumption. In that case the gut equation is augmented with an input term, with an absorption rate constant ''K
a'':
That requires defining an equation for the quantity ingested and present in the gut lumen:
In the absence of a gut compartment, input can be made directly in the liver. However, in that case local metabolism in the gut may not be correctly described. The case of approximately continuous absorption (''e.g. via'' drinking water) can be modeled by a zero-order absorption rate (here ''R
ing'' in units of mass over time):
More sophisticated gut absorption model can be used. In those models, additional compartments describe the various sections of the gut lumen and tissue. Intestinal pH, transit times and presence of active transporters can be taken into account
.
Skin depot
The absorption of a chemical deposited on skin can also be modeled using first order terms. It is best in that case to separate the skin from the other tissues, to further differentiate exposed skin and non-exposed skin, and differentiate viable skin (dermis and epidermis) from the
stratum corneum
The stratum corneum (Latin for 'horny layer') is the outermost layer of the epidermis. The human stratum corneum comprises several levels of flattened corneocytes that are divided into two layers: the ''stratum disjunctum'' and ''stratum compact ...
(the actual skin upper layer exposed). This is the approach taken in [Bois F., Diaz Ochoa J.G. Gajewska M., Kovarich S., Mauch K., Paini A., Péry A., Sala Benito J.V., Teng S., Worth A., in press, Multiscale modelling approaches for assessing cosmetic ingredients safety, Toxicology. doi: 10.1016/j.tox.2016.05.026]
Unexposed
stratum corneum
The stratum corneum (Latin for 'horny layer') is the outermost layer of the epidermis. The human stratum corneum comprises several levels of flattened corneocytes that are divided into two layers: the ''stratum disjunctum'' and ''stratum compact ...
simply exchanges with the underlying viable skin by diffusion:
where
is the partition coefficient,
is the total skin surface area,
the fraction of skin surface area exposed, ...
For the viable skin unexposed:
For the skin
stratum corneum
The stratum corneum (Latin for 'horny layer') is the outermost layer of the epidermis. The human stratum corneum comprises several levels of flattened corneocytes that are divided into two layers: the ''stratum disjunctum'' and ''stratum compact ...
exposed:
for the viable skin exposed:
dt(QSkin_u) and dt(QSkin_e) feed from arterial blood and back to venous blood.
More complex diffusion models have been published
eference to add
Intra-venous injection
Intravenous injection is a common clinical route of administration. (to be completed)
Inhalation
Inhalation occurs through the lung and is hardly dissociable from exhalation (to be completed)
Modelling metabolism
There are several ways metabolism can be modeled. For some models, a linear excretion rate is preferred. This can be accomplished with a simple differential equation. Otherwise a
Michaelis-Menten equation, as follows, is generally appropriate for a more accurate result.
:
.
Uses of PBPK modeling
PBPK models are compartmental models like many others, but they have a few advantages over so-called "classical" pharmacokinetic models, which are less grounded in physiology. PBPK models can first be used to abstract and eventually reconcile disparate data (from physicochemical or biochemical experiments, ''
in vitro
''In vitro'' (meaning in glass, or ''in the glass'') studies are performed with microorganisms, cells, or biological molecules outside their normal biological context. Colloquially called "test-tube experiments", these studies in biology an ...
'' or ''
in vivo
Studies that are ''in vivo'' (Latin for "within the living"; often not italicized in English) are those in which the effects of various biological entities are tested on whole, living organisms or cells, usually animals, including humans, and ...
'' pharmacological or toxicological experiments, ''etc.'') They give also access to internal body concentrations of chemicals or their metabolites, and in particular at the site of their effects, be it therapeutic or toxic. Finally they also help interpolation and extrapolation of knowledge between:
*Doses: e.g., from the high concentrations typically used in laboratory experiments to those found in the environment
*Exposure duration: ''e.g.'', from continuous to discontinuous, or single to multiple exposures
*Routes of administration: ''e.g.'', from inhalation exposures to ingestion
*Species: ''e.g.'', transpositions from rodents to human, prior to giving a drug for the first time to subjects of a clinical trial, or when experiments on humans are deemed unethical, such as when the compound is toxic without therapeutic benefit
*Individuals: ''e.g.'', from males to females, from adults to children, from non-pregnant women to pregnant
*
From ''in vitro'' to ''in vivo''.
Some of these extrapolations are "parametric" : only changes in input or parameter values are needed to achieve the extrapolation (this is usually the case for dose and time extrapolations). Others are "nonparametric" in the sense that a change in the model structure itself is needed (''e.g.'', when extrapolating to a pregnant female, equations for the foetus should be added).
Owing to the mechanistic basis of PBPK models, another potential use of PBPK modeling is hypothesis testing. For example, if a drug compound showed lower-than-expected oral bioavailability, various model structures (i.e., hypotheses) and parameter values can be evaluated to determine which models and/or parameters provide the best fit to the observed data. If the hypothesis that metabolism in the intestines was responsibility for the low bioavailability yielded the best fit, then the PBPK modeling results support this hypothesis over the other hypotheses evaluated.
As such, PBPK modeling can be used, ''inter alia'', to evaluate the involvement of carrier-mediated transport, clearance saturation, enterohepatic recirculation of the parent compound, extra-hepatic/extra-gut elimination; higher ''in vivo'' solubility than predicted ''in vitro''; drug-induced gastric emptying delays; gut loss and regional variation in gut absorption.
Limits and extensions of PBPK modeling
Each type of modeling technique has its strengths and limitations. PBPK modeling is no exception. One limitation is the potential for a large number of parameters, some of which may be correlated. This can lead to the issues of parameter identifiability and redundancy. However, it is possible (and commonly done) to model explicitly the correlations between parameters (for example, the non-linear relationships between age, body-mass, organ volumes and blood flows).
After numerical values are assigned to each PBPK model parameter, specialized or general computer software is typically used to numerically integrate a set of ordinary differential equations like those described above, in order to calculate the numerical value of each compartment at specified values of time (see Software). However, if such equations involve only linear functions of each compartmental value, or under limiting conditions (e.g., when input values remain very small) that guarantee such linearity is closely approximated, such equations may be solved analytically to yield explicit equations (or, under those limiting conditions, very accurate approximations) for the time-weighted average (TWA) value of each compartment as a function of the TWA value of each specified input (see, e.g.,).
PBPK models can rely on chemical property prediction models (
QSAR models or predictive chemistry models) on one hand. For example, QSAR models can be used to estimate partition coefficients. They also extend into, but are not destined to supplant,
systems biology
Systems biology is the computational modeling, computational and mathematical analysis and modeling of complex biological systems. It is a biology-based interdisciplinary field of study that focuses on complex interactions within biological syst ...
models of metabolic pathways. They are also parallel to
physiome The physiome of an individual's or species' physiological state is the description of its functional behavior. The physiome describes the physiological dynamics of the normal intact organism and is built upon information and
structure (genome, prote ...
models, but do not aim at modelling physiological functions beyond fluid circulation in detail. In fact the above four types of models can reinforce each other when integrated.
References
Further references:
*
*
*
*
*Reddy M. ''et al.'' (2005) ''Physiologically Based Pharmacokinetic Modeling : Science and Applications'', Wiley-Interscience.
*Peters S.A (2012) ''Physiologically-Based Pharmacokinetic (PBPK) Modeling and Simulations'', Wiley.
{{refend
Forums
Ecotoxmodelsis a website on mathematical models in ecotoxicology.
Software
Dedicated software:
BioDMET*
GastroPlus Maxsim2 PK-Sim PKQuestPSE: gCOAS*
Simcyp Simulator
Simcyp Limited is a research-based company which provides modelling and simulation software to the pharmaceutical industry for use during drug development. Simcyp is based in Sheffield, UK.
Simcyp’s Simulators allow '' in silico'' prediction of ...
ADME Workbench
General software:
ADAPT 5*
Berkeley Madonna
Berkeley Madonna is a mathematical modelling software package, developed at the University of California at Berkeley by Robert Macey and George Oster. It numerically solves ordinary differential equations and difference equations, originally deve ...
COPASI: Biochemical System Simulator*
Ecolego
Ecolego is a simulation software tool that is used for creating dynamic models and performing deterministic and probabilistic simulations. It is also used for conducting risk assessments of complex dynamic systems evolving over time.
Ecolego c ...
Free simulation software GNU MCSIM GNU Octave MatlabPottersWheel
PottersWheel is a MATLAB toolbox for mathematical modeling of time-dependent dynamical systems that can be expressed as chemical reaction networks or ordinary differential equations (ODEs). It allows the automatic calibration of model parameters ...
ModelMaker PhysioLab R deSolve package SAAM II Phoenix WinNonlin/NLME/IVIVC/Trial Simulator
Toxicology
Pharmacokinetics
Pharmaceutics