Next Article in Journal
Artificial Intelligence/Machine Learning-Driven Small Molecule Repurposing via Off-Target Prediction and Transcriptomics
Previous Article in Journal
Effectiveness of Tap Water in Reducing the Generation of Ultrafine Wear Particles from the Wheel-Rail Contact by Eliminating the Water Vapor Effect
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment

1
Key Laboratory of Ministry of Education for Gastrointestinal Cancer, School of Basic Medical Sciences, Fujian Medical University, Fuzhou 350108, China
2
Department of Toxicology, School of Public Health, Peking University, Beijing 100191, China
3
Department of Pathology, School of Basic Medical Sciences, Ningxia Medical University, Yinchuan 750004, China
4
Department of Gynecology, Fujian Maternity and Child Health Hospital (Fujian Obstetrics and Gynecology Hospital), Fuzhou 350001, China
5
Key Laboratory of State Administration of Traditional Chinese Medicine for Compatibility Toxicology, Beijing 100191, China
6
Beijing Key Laboratory of Toxicological Research and Risk Assessment for Food Safety, Beijing 100191, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Toxics 2023, 11(10), 874; https://doi.org/10.3390/toxics11100874
Submission received: 13 September 2023 / Revised: 17 October 2023 / Accepted: 19 October 2023 / Published: 21 October 2023
(This article belongs to the Section Exposome Analysis and Risk Assessment)

Abstract

:
Toxicokinetics plays a crucial role in the health risk assessments of xenobiotics. Classical compartmental models are limited in their ability to determine chemical concentrations in specific organs or tissues, particularly target organs or tissues, and their limited interspecific and exposure route extrapolation hinders satisfactory health risk assessment. In contrast, physiologically based toxicokinetic (PBTK) models quantitatively describe the absorption, distribution, metabolism, and excretion of chemicals across various exposure routes and doses in organisms, establishing correlations with toxic effects. Consequently, PBTK models serve as potent tools for extrapolation and provide a theoretical foundation for health risk assessment and management. This review outlines the construction and application of PBTK models in health risk assessment while analyzing their limitations and future perspectives.

1. Introduction

Toxicokinetics (TK) is an interdisciplinary field that integrates the principles and methodologies of pharmacokinetics and toxicology. It plays a crucial role in understanding the absorption, distribution, metabolism, excretion, and toxicity (ADMET) characteristics of xenobiotics in both human and animal models, considering various exposure routes, doses, and frequencies. TK data enables the calculation of kinetic parameters and facilitates adjustments to experimental designs. By establishing the relationship between dose-exposure and exposure-toxicity, TK allows for the assessment of toxicity effects, mechanisms in target organs and tissues, and the associated risks to human health.
Traditional TK studies often involve doses higher than those used in clinical pharmacokinetic research, resulting in saturation of dissolution, absorption, and metabolism processes and leading to nonlinear kinetics models [1,2]. Therefore, the bioavailability, half-life, apparent volume of distribution, clearance, and other parameters calculated from low-dose pharmacokinetic studies may not be applicable to high-dose toxicokinetic investigations. To enable reliable interspecific extrapolation, dose extrapolation, and in vivo-in vitro correlation, it is essential to systematically consider dose-response and dose-effect relationships in TK studies. After a satisfactory evaluation, the PBTK model is utilized to extrapolate and compute internal doses, thereby enhancing the dose-response relationship in toxicity testing and risk assessment [3]. The European Human Biomonitoring Initiative effectively employed the PBTK model to evaluate the association between chemical exposure and observed HBM data [4]. In the drug review process conducted by the United States Food and Drug Administration (FDA), submitted PBTK models from drug developers were assessed for adequacy at different stages of drug development to achieve desired objectives. Previous studies had solely focused on developing rat-based PBTK models to assess risks associated with oral exposure to diethyl phthalate (DEP), and Hu et al. [5] developed a human dermal pathway-specific model for assessing health risks when exposed to DEP, which suggested that dermal exposure to DEP posed a greater risk to human health compared with oral exposure.
In general, four conventional methods are employed for the calculation of pharmacokinetic or toxicokinetic parameters, including the compartmental model, statistical moment analysis of the noncompartmental model, population pharmacokinetics or toxicokinetic analysis, and the physiologically based pharmacokinetics or toxicokinetic (PBPK/PBTK) model. However, the first three methods do not consider physiological characteristics and, thus, can only reflect concentration changes over time in specific tissues or organs, lacking the ability to analyze the distribution and metabolism of xenobiotics simultaneously. As a result, these methods fail to extrapolate prediction results across different species, exposure routes, and doses [1,6,7]. In contrast, PBTK models represent a more complex form of compartmental models, where tissue and organ entities act as interconnected compartments through blood circulation. The parameters in PBTK models possess physiological significance, such as blood flow volume and tissue volume, enabling the simulation of xenobiotic absorption, distribution, metabolism, and excretion in various tissues or organs, thereby facilitating extrapolation. This review primarily focuses on constructing and applying the PBTK model in health risk assessments while analyzing its limitations and future perspectives.

2. The Construction of the PBTK Model

The construction of a PBTK model typically involves five steps, as illustrated in Figure 1: (1) model characterization, (2) definition of the model parameters, (3) model simulation, (4) model evaluation, and (5) model optimization [1,7,8,9].

2.1. Model Characterization

2.1.1. Determination of the Model Structure

The structure of the PBTK model is based on physiological and anatomical considerations; it can be either a whole-body PBTK model or a semi-PBTK model, depending on the complexity required. In the whole-body PBTK model (Figure 2A), each organ or tissue is considered as an independent compartment, connected through blood flow to others [10]. In the semi-PBTK model (Figure 2B), major organs like blood, liver (the main metabolic organ), and kidneys (the main excretory organ) are separate compartments. Other organs not relevant to the study are grouped into one compartment, simplifying the structure and reducing the number of parameters for improved parameter fitting accuracy. Figure 2B illustrates a common human semi-PBTK model with a respiratory exposure route. In this model, organs or tissues not directly relevant to the study are categorized as richly perfused tissues (e.g., brain, heart, spleen, intestines, kidneys, adrenal glands, thyroid, lungs, and bone marrow) or poorly perfused tissues (e.g., muscle and skin) [11].
During the modeling process, consideration must be given to which organs or tissues should be included. Firstly, the PBTK model includes compartments for blood, metabolic organs (such as the liver), and excretory organs (such as the kidneys), as they are closely associated with toxicokinetic studies. Secondly, it is necessary to determine whether other organs or tissues should be included in the chemical ADMET processes. For example, in the case of oral administration, compartments for the stomach and gut should be added. For lipophilic chemicals with high octanol-water partition coefficients (log Kow > 3), an independent compartment for adipose tissue should be included. Thirdly, adherence to the parsimony principle is crucial [9]. As model complexity increases, more data are required for fitting and validation purposes. In the field of chemical toxicology, greater attention has been focused on target organs or tissues such as the liver and kidney, while fewer studies have explored other organs or tissues comprehensively. Due to insufficient data availability that can support a complex PBTK model effectively, it is essential to simplify the model as much as possible in order to enhance its predictive reliability.

2.1.2. Building Mathematical Equations

After determining the model structure, mathematical equations were established to describe the ADME processes of xenobiotics in each organ or tissue. The establishment of the model should follow the principle of mass equation [7,9,12]:
d A t d t = A t , i n A t , o u t A t , e A t , m
where dAt/dt is the amount if change of xenobiotics in a compartment per unit time; At,in is the amount of entry into the compartment; At,out is the quantity leaving the compartment; At,e is the amount excreted by the compartment; and At,m is the amount of compartment metabolism.
In the PBTK model, the absorption process of xenobiotics by organs or tissues obeys the simple diffusion principle in Fick’s law [1,11]. Based on the distribution characteristics of xenobiotics in organs or tissues, PBTK models can be categorized into three types:
(1)
Perfusion-limited model: in this model, xenobiotic concentrations in organs, tissues, and blood reach instant equilibrium without concentration differences. Organs, tissues, and blood are considered as a single compartment with a homogeneous distribution of xenobiotics. The only limiting factor for xenobiotic distribution is blood flow velocity [12,13]. This model is suitable for small lipid-soluble molecules that can readily cross membrane barriers and when organs or tissues are small-sized with high blood flow [7,14,15];
(2)
Permeability-limited model: in this model, xenobiotics penetrate organs or tissues through membrane barriers, resulting in a concentration gradient between organs or tissues and blood. Depending on the number of membrane barriers, the permeability-limited model can have two or three sub-compartments. This model is suitable for molecules with polarity or large molecular weight [14,15];
(3)
Dispersion model: in this model, xenobiotics are distributed in organs or tissues with a gradient, and the degree of dispersion is evaluated using the dispersion coefficient (DN). A higher DN indicates a greater dispersion of xenobiotics in organs or tissues. When DN approaches infinity, the dispersion model is similar to the perfusion-limited model [16,17,18,19,20,21]. This model is suitable for xenobiotics with high hepatic clearance [8].
The metabolism of xenobiotics can be described as first-order kinetic, second-order kinetic, and saturation process [11,15]. Table 1 provides commonly used mathematical equations to describe the ADME processes in PBTK models [15].

2.2. Definition of Model Parameters

As shown in Table 2, the PBTK model mainly contains physiological, physicochemical, and biochemical parameters. The physiological parameters include animal body weight, cardiac output, blood flow through an organ or tissue, the volume of the organ or tissue, and so on [22,23,24,25]. The physiological parameters are obtained from the literature, depending on the specific study objectives, or determined through experiments [26,27,28,29]. Physiological parameters can be influenced by factors such as age, sex, and diseases.
Physicochemical parameters refer to the blood-air distribution coefficient and tissue-blood distribution coefficient, which represent the concentration ratio of xenobiotics in different carriers at steady state and vary with the properties of the xenobiotics. For xenobiotics with insufficient documented physicochemical parameters, data can be obtained through three approaches. Firstly, in vivo experiments can provide distribution coefficients based on steady-state toxicokinetic data from repeated administration or dynamic data from single intravenous administration [11,30,31,32,33,34]. Secondly, in vitro experiments such as vial equilibration [35,36,37,38,39,40,41,42,43,44], equilibrium dialysis [45,46,47], ultrafiltration [46], and others have been reported for determining distribution coefficients. Thirdly, in silico approaches involve computation and simulation methods that determine distribution coefficients based on structural characteristics and physicochemical properties of xenobiotics [48,49,50,51,52,53,54,55,56,57,58].
Biochemical parameters include the absorption rate constant, maximum velocity, Michaelis constant, etc. These parameters can be calculated using plasma or tissue xenobiotic concentration versus time profiles and excretion data from in vivo [59] or in vitro [60,61,62,63,64,65] experiments. They can also be obtained through computation and simulation using quantitative structure–activity relationship models based on the molecular structure of xenobiotics [66,67,68,69]. It is important to note that metabolic parameters obtained from in vitro experiments cannot be directly used in an in vivo PBTK model. Appropriate data transformation should be performed, considering the metabolic differences between in vitro and in vivo studies [11,70,71]. In order to identify the most sensitive parameters within existing mammalian models, Schneckener et al. [72] extrapolated a validated rabbit PBTK model to six other mammalian species using established species-specific models. By calculating sensitivity across various mammalian species, different compounds, and diverse administration routes, only 5.4% (n = 33) parameters had sensitivity values over 0.4.

2.3. Model Simulation

PBTK model simulation involves generating predictive curves of toxicokinetics that fit with measured data by solving ordinary differential equations through mathematical calculations and software.

2.3.1. Algorithm

Common algorithms used to solve ordinary differential equations are Euler, Gear, Runge-Kutta, and predictor-corrector. These algorithms follow the principle:
N e w   v a l u e = o l d   v a l u e + s l o p e × d t
where dt is the calculated step value, and slope is the first derivative of the equation of a curve at a given point. For a particular compartment in the PBTK model,
A t , 1 = A t , 0 + ( d A t / d t ) × d t
where A is the total amount of xenobiotics, and the subscripts t,0 and t,1 represent the start and end times of each calculation, respectively.
The Euler is the most frequently used algorithm to solve ordinary differential equations in the PBTK model [73,74]. The error generated by the Euler algorithm is proportional to the square of the step value (dt2), so a smaller dt2 is needed to reduce this error, although it requires more calculation time. In addition, the Euler algorithm is not applicable to stiff system solutions, for which the Gear algorithm is recommended due to its greater stability [11].

2.3.2. Software

Existing software for PBTK modeling and analysis can be divided into two categories; the first is general mathematical and engineering modeling software, such as MATLAB, Berkeley Madonna, R, and acslX. These software packages require users to build a model framework and write corresponding codes, necessitating high modeling and programming skills. The advantage of these software packages is their high flexibility [75], allowing users to construct customized PBTK models based on their needs and perform complex analyses such as uncertainty analysis, Monte Carlo analysis, and Markov chain Monte Carlo (MCMC) analysis. The second is the professional PBTK and modeling software, such as Simcyp, PK-Sim, GastroPlus, Cloe PK, PKQuest, etc. This software is user-friendly and requires minimal modeling and programming skills. Users can either combine various components of the model or directly use built-in models provided by the software. Some software packages also integrate functions for predicting the physicochemical properties of xenobiotics. These software options are suitable for researchers who are not interested in computation and programming or lack modeling experience [7,8,9,76,77]. Mahdi et al. [78] investigated a suitable solvent for subcutaneous delivery of rifampicin through in silico predictions by GastroPlus, which reduced the cost and time required for formulation development by utilizing in vivo data to simulate in vitro findings. Hanke et al. [79] developed a rosuvastatin model using the open-source PBTK software PK-Sim, incorporating plasma, urine, and feces data of rosuvastatin, positron emission tomography measurements of tissue concentrations, and seven different rosuvastatin drug–drug interaction studies. The PK-Sim expression database provided the relative expression in different organs by analyzing drug transporters and metabolizing enzymes. However, there are certain limitations when using software to predict drug metabolism. These software have a high threshold and need users to receive sufficient professional knowledge on the data and parameters that are to be used in the implementation. The establishment of a PBTK model needs numerous drug parameters, so further validation and optimization of the constructed model are essential to improve its accuracy.

2.4. Model Evaluation

A constructed PBTK model should undergo evaluation under various scenarios, particularly its ability to predict experimental data not used in the modeling process. Additionally, uncertainty, variability, and sensitivity analyses are performed [7,11].

2.4.1. Verification of Measured Data

The predictive performance of the PBTK model is evaluated through visual inspection, statistical tests, and discrepancy indices by comparing predicted values to measured values [1,11,80,81]. Based on visual inspection, the most commonly used method, Figure 3A exhibits superior predictive ability compared to Figure 3B due to the logarithmic scale typically used for the concentration versus time curve.
In the context of PBTK model validation, statistical hypothesis tests and discrepancy tests are not applicable because the null hypothesis assumes that the model is identical to the actual biological system, which is not realistic [80]. Classical statistical methods such as Mann-Whitney, two-sample χ2, and two-sample Kolmogorov are also not suitable for assessing consistency between predicted and measured values due to the self-correlation of concentration values at different times in pharmacokinetics [11]. To assess the predictive ability of the PBTK model, linear regression analysis can be performed to compare the predicted and measured values. A high prediction ability is indicated when the intercept of the regression equation is close to 0 and the slope and correlation coefficient are close to 1 [10,77]. Another method for assessing the predictive ability of the PBTK model is to calculate the chi-square χ2 as Equation (16). This approach is applicable when there are multiple measured values at the same time point. Peters used this method in an overall PBTK model to evaluate the model fitting [82]. An χ2 value close to 1 indicates better model fitting and prediction performance. The Akaike information criterion (AIC) and Bayesian information criterion (BIC) are commonly used for assessing predictive ability [83,84]. AIC and BIC are calculated using Equations (17) and (18), respectively. Smaller values of AIC and BIC indicate a better model fit. When there is a significant difference between the fitted model and the real model, it is primarily reflected in the likelihood function term −2In(L). When the difference in the likelihood function is not significant, the penalty term 2k, accounting for the number of model parameters, becomes influential. Models with better fit and fewer parameters are preferred as the increase in the number of parameters leads to increased 2k and AIC values. The penalty term for BIC is larger than that for AIC, taking into account the number of observations. This prevents excessive model complexity when there is a large number of observations. The last method is to calculate the fold error (FE), as Equation (19). Firstly, the noncompartmental analysis was performed to obtain the peak concentration, peak time, area under the curve, and other parameters of the observations. Then, the fitting values or predicting values of the above parameters were obtained from corresponding models. Due to the different laboratory conditions, detection methods, and other confounding factors, the measured value and the predicted value cannot be consistent completely, and the model is acceptable when the FE ≤ 2 [80].
χ 2 = 1 n i = 1 n ( Δ i 2 σ i 2 )  
In Equation (16), n is the number of observations, Δ is the residual, i.e., the difference between the measured and predicted values, and σ is the standard deviation of the measured values at the same time point.
AIC = 2 k 2 In L  
BIC = In n · k 2 In L  
In Equations (17) and (18), k is the number of model variables, L is the likelihood function, and n is the number of observations.
FE = 10 l o g S i m u l a t e d O b s e r v e d  

2.4.2. Uncertainty Analysis

The PBTK model is constructed according to the physiological and anatomical features artificially, and these parameters are obtained from in vitro and in vivo experiments or in silico predictions. An illogical model structure or inaccurate parameters will result in an inaccurate prediction [7,11]; therefore, the uncertainty analysis is performed to assess the impact of parameter errors on model prediction results. Three methods are used for uncertainty analysis; the first is the Monte Carlo simulation, which calculates an output parameter value by randomly selecting input parameter values fitting specific probability distribution. Cox [85] employed Monte Carlo uncertain analysis within a PBTK model to re-evaluate the risk assessments of benzene metabolism in animals and humans, and the probability distribution of the entire dose–response function for benzene-induced tumors in male B6C3F1 mice was quantified. The second method is probability-bound analysis; in this approach, input parameters are systematically varied from minimum to slightly above minimum, then set at the normal value, followed by slightly below the maximum, and finally at the maximum of the distribution. Concurrently, the corresponding output parameters are observed, representing extreme values of parameter variation. The third is the Fuzzy simulation, in which the parameters are fuzzed after input, and the output is obtained by fuzzy rules. If the distributions of uncertain parameters are known, the Monte Carlo simulation method can be used [86,87,88,89,90]. Otherwise, the fuzzy simulation is suggested [91,92,93]. If only the maximum and minimum values of the parameters are known, the probability-bound analysis is suggested [11,94]. Through the uncertainty analysis, the impact of input parameters on the model is elucidated and improves the prediction ability of the model.

2.4.3. Variation Analysis

When constructing the PBTK model for population studies, the physiological and biochemical parameters are variable due to individual differences. The variation analysis is used to analyze the parameter variability in the PBTK model during the assessment of predicting ability in population toxicokinetics. The methods for uncertainty analysis are also applicable to variation analysis, and the MCMC is the classic method for variability analysis [95,96,97,98]. In their study on the maximum detection time (MDT) of 1,1-difluoroethane in the blood after inhalation abuse among the adult male population in the United States, Huet et al. [99] introduced variability in huffing pattern and body mass index using the Monte Carlo simulation within the PBTK model, and the results indicated that the MDT of 1,1-difluoroethane in blood after abuse ranged from 7.8 to 15.8 h.

2.4.4. Sensitivity Analysis

The sensitivity analysis reflects the impact of parameter changes on model prediction results and is expressed by sensitivity ratio as Equation (16) [76]:
S = O I × I O
In Equation (20), S is the sensitivity value representing the percentage change of the corresponding output value when a certain input parameter changes by 1% and other input parameters remain unchanged. O is the output prediction result of the model, I is the input parameter of the model. S 0.5 means a highly sensitive parameter; 0.2 S < 0.5 means a medium sensitivity parameter; 0.1 S < 0.2 means a low sensitivity parameter; S < 0.1 means an insensitive parameter, and the parameter sensitivity is not significant.
To analyze the interaction effects of parameters on the prediction results, the above sensitivity analysis method is not applicable, and the global sensitivity analysis is suggested. Campolongo and Saltelli proposed the two-step method for global sensitivity analysis. Firstly, the Morris method was utilized for screening the parameters potentially involved in interactions; then, the effects of these interactions were quantified using the extended Fourier amplitude sensitivity test [100], a reliable and efficient method suitable for small sample sizes [101,102]. A high sensitivity value of a certain parameter means that a slight change of the parameter will dramatically affect the prediction results. The parameters with high sensitivity are considered the key factors that should be accurately controlled to improve the prediction accuracy of the models. Jeong et al. [103] used a PBTK approach, coupled with global sensitivity analysis, to predict and evaluate the pharmacokinetic profiles of nafamostat in a virtual healthy population under various dosing regimens, and the results suggested that the hepatic distribution and metabolism had a significant impact on systemic exposure and clearance of nafamostat.

2.4.5. Model Optimization

When the predicted value differs greatly from the measured value, the model structure or parameters may be inappropriate, and model optimization should be performed. In vivo experiments are suggested to clarify the ADME process for a reasonable model structure. For the prediction deviation caused by inappropriate parameters, the uncertainty, variability, and sensitivity analysis can be performed to determine the parameters that greatly impact the prediction value. Then, the parameters are optimized by in vivo, in vitro, and in silico methods to improve the prediction accuracy of the models.

3. Applications of Physiologically Based Toxicokinetics (PBTK) in Health Risk Assessment

3.1. Exposure Assessment

Forward dosimetry is used to predict the internal dose (tissular chemical concentration) by external dose (exposure chemical concentration) and is called reverse dosimetry and vice versa [104]. The PBTK model establishes a reliable correlation between the external dose and internal dose so that reverse dosimetry is particularly suitable for exposure assessment [105]. It is suitable for practical applications in exposure assessment, such as population-based studies on herbicide fluorine [106]. If the external dose is linear with the internal dose, it can be calculated by the following equation:
Exposure   dose = C X  
where C is the measured concentration of the chemical or metabolite in tissues, and [X] is the theoretical value in a biological sample at a simulated exposure dose of one unit. [X] represents the maximum concentration, area under the curve, or other toxicokinetic parameters.

3.2. Extrapolation

Compared to classical compartment models and population pharmacokinetic models, the PBTK model, which is based on physiological and anatomical features, is more suitable for extrapolation. Many parameters in the PBTK model cannot be directly obtained through in vivo experiments; therefore, in vitro and in silico assays are necessary. In the PBTK model, various extrapolations such as in vitro-in vivo extrapolation (IVIVE), dose extrapolation, exposure route extrapolation, and interspecific extrapolation are performed to predict toxicokinetic processes across different species and exposure conditions. This approach significantly reduces the reliance on animal experiments while conserving human and material resources, aligning with the “3R” principle of animal ethics.

3.2.1. IVIVE

The in vitro cellular assays used for predicting the in vivo ADME processes of chemicals in IVIVE do not account for active transport, biotransformation, and other factors, leading to uncertainty in this extrapolation method. Poulin et al. [107] highlighted the need for adjustments when applying intracellular free concentrations obtained from in vitro assays to in vivo studies. While novel cell culture methods such as 3D cell culture and coculture provide better simulation of the in vivo state compared to traditional monolayer cell culture, their complexity and cost restrict their utilization for high-throughput screening. In this context, the PBTK model used in IVIVE offers a more accurate and convenient prediction of chemical ADME processes in vivo [108].
The reverse dosimetry approach based on PBTK modeling enables establishing a relationship between in vitro concentrations and human exposure doses, facilitating a dose-response assessment using IVIVE. In their study, Fabian et al. [109] employed the PBTK model to extrapolate the lowest observed effect concentrations of ten environmental endocrine disruptors determined in vitro to the equivalent oral doses in rats and compared them with the lowest observed effect levels (LOEL) values derived from animal studies. Results indicated that predicted LOEL values for six out of ten chemicals were consistent with those obtained from animal experiments within an acceptable range, while the remaining four differed by more than 10-fold. This investigation highlights the potential utility of reverse dosimetry based on the PBTK model for conducting IVIVE analysis, although further improvements in accuracy are needed.
Dejongh et al. [110] developed a PBTK model using data from in vitro assays, such as tissue-blood partition coefficients and liver metabolism parameters. The in vivo neurotoxicity of eight neurotoxicants (benzene, toluene, lindane, acrylamide, parathion/paraoxon, caffeine, diazepam, and phenytoin) was predicted and compared to literature-based in vivo values. The predicted values showed a correlation with the experimental values, with greater accuracy observed for low neurotoxic chemicals (approximately 2-fold deviation) compared to high neurotoxic chemicals (approximately 10-fold deviation). This study suggested the potential application of PBTK models in IVIVE studies. Brinkmann et al. [111] focused on two classes of chemicals, using 7-ethoxyresorufin-O-deethylase and vitellogenin as exposure biomarkers in Oncorhynchus mykiss. The half-effect concentrations (EC50) of these chemicals were determined through both in vitro and in vivo assays. The in vitro EC50 was used as a dose parameter in the PBTK model to predict chemical concentrations in the blood and liver of Oncorhynchus mykiss. Linear correlations were observed between in vitro EC50 and the concentrations of chemicals in the blood and liver, indicating that combining the PBTK model with in vitro assays can accurately predict results from in vivo experiments.

3.2.2. Dose Extrapolation

In toxicology experiments, a high dose exposure is commonly employed; however, actual exposure to xenobiotics typically occurs at low doses. Consequently, dose extrapolation becomes crucial in risk assessment, particularly when transitioning from high to low doses. Simple linear extrapolation often leads to significant deviations due to potential saturation of absorption, metabolism, and excretion during high dose exposure. To address this issue more accurately, the PBTK model incorporates nonlinear equations such as the Michaelis-Menten equation to describe the nonlinear processes of chemicals in vivo. For example, Wang et al. [112] employed the PBTK model to simulate the toxicokinetic of vinyl chloride in the human body after exposure to 100 ppm for 6 h. Subsequently, this model was utilized to predict the exhaled breath concentration of vinyl chloride following exposures of 59 ppm and 261 ppm for 7.5 h. The PBTK model’s predicted values aligned closely with the measured values, thereby demonstrating its precise predictive capability for in vivo concentrations of inhaled vinyl chloride.

3.2.3. Exposure Route Extrapolation

Exposure route extrapolation involves developing a model based on intravenous administration and incorporating absorption processes to extrapolate to non-intravenous routes such as oral, subcutaneous, or inhalation exposure. Traditional extrapolation methods often assume complete absorption in different routes, disregarding factors like the first-pass effect in oral administration. The PBTK model can simulate the absorption process in various exposure routes, offering a more scientifically sound and accurate extrapolation approach.
Gajewska et al. [113] employed a PBTK model to extrapolate the oral no observed adverse effect level (NOAEL) of three cosmetic ingredients (coumarin, hydroquinone, and caffeine) in animal studies to estimate the corresponding oral NOAEL in humans. Based on this, they predicted the concentration in the human body following percutaneous absorption. The findings revealed that for low exposure doses, the oral route did not consistently result in higher in vivo concentrations compared to the percutaneous route.

3.2.4. Interspecific Extrapolation

In toxicological experiments, ethical considerations often require testing on experimental animals, and the findings need to be extrapolated to humans for risk assessment purposes. Traditional extrapolation methods rely on scaling factors such as body weight or body surface area, which fail to account for metabolic differences between species and consequently introduce significant uncertainty. Interspecific extrapolation using the PBTK model allows the accurate prediction of toxicokinetic processes across different species by adjusting species-specific parameters such as physiological and metabolic rate parameters.
The human PBTK model for acrylamide and its metabolite epoxyproacrylamide was developed by Sweeney et al. [114], based on the rat PBTK model, with adjustments made to certain parameters using human data. The model successfully predicted the toxicokinetic processes of these compounds in humans, yielding results consistent with literature measurements. Furthermore, this model was employed to investigate the safety dose of acrylamide in humans.

3.3. PBTK of the Special Population

The primary distinction between the special population and the general population lies in the alterations observed in physiological structure and parameters. Consequently, a PBTK model designed for a specific population can be developed. This involves modifying the general population model framework, replacing standard physiological parameters with those specific to the particular population, and integrating distinctive mechanisms related to specific chemicals. Commonly employed PBTK models for special populations encompass three categories: pregnant women and fetuses, children, as well as patients afflicted with hepatic or renal impairment.
For pregnant women and fetuses, they are often modeled together rather than separately in PBTK models. Abduljalil et al. [115] investigated the anatomical, physiological, and metabolic changes occurring in pregnant women throughout gestational weeks. Studies have demonstrated significant alterations in crucial parameters such as cardiac output, protein binding, and metabolic enzyme activity during pregnancy. Moreover, many physiological parameters change non-linearly over time and require description through distinct algorithms. Xia et al. [116] presented a relatively simple structure of the PBTK model for pregnancy, incorporating essential physiological parameters like renal function and CYP450 enzyme system activity. The model was validated using four reference drugs. Wu et al. [117] expanded the PBTK model constructed by Xia et al. by adding physiological parameters such as fetal placenta, fat, and plasma volume for the application in pregnant women. Dallmann et al. [118] provided a comprehensive review of existing PBTK models during pregnancy, characterizing the model structures for pregnant women and fetuses. In this model, the fetus was connected to the mother through the umbilical artery and vein of the placenta. Corley et al. [119] assessed the applicability of PBTK models for pregnant and lactating women in child risk assessment. The PBTK model developed by Lu et al. [120] integrated variations in enzyme activity and accurately predicted the therapeutic doses of three drugs, namely caffeine, metoprolol, and midazolam, which were metabolized by CYP1A2, CYP2D6, and CYP3A4 respectively, in pregnant women. The findings revealed a two-fold increase in caffeine exposure among pregnant women compared to their non-pregnant counterparts.
For children, physiological parameters vary significantly across the different stages of development from infancy to 18-year-old adolescence. Therefore, it is more appropriate to express these parameters as age-related. Edgington et al. [121] summarized changes in relevant physiological parameters in children, such as body weight, height, organ weight, and organ-specific blood flow. They modified adult PBTK models by incorporating age-related adjustments and proposed age-corrected models to predict chemical concentrations in plasma for children at corresponding ages. Kovar et al. [122] developed an adult PBTK model for buprenorphine and norbuprenorphine, and this model was also suitable for children and preterm neonates since it considered age-related changes. Khalil et al. [123] provided a summary of the methodology, application, and limitations of PBTK models for pediatric medication. Generally, the PBTK model proves valuable in simulating and predicting pharmacokinetics after medication administration in children by providing guidance on optimal timing, dosage, and frequency while avoiding ethical concerns. However, there are limitations due to the incompleteness of physiological process descriptions and variability and uncertainty in the population, limiting the clinical application of these models.
In the case of hepatic and renal impairment, the primary impact is on the clearance of chemicals from the body, leading to a slower decrease in plasma concentration. To predict the temporal profile of chemical concentrations in patients with hepatic and renal impairment, it is necessary to adjust physiological parameters affected by liver and kidney damage without modifying the model structure [124]. Compared to a healthy state, parameters such as hepatic and renal blood flow, liver and kidney function, plasma protein binding rate, enzyme activity, and elimination processes are primarily affected. Furthermore, enzyme activity, CL, and biliary excretion are affected during oral exposure in cases of hepatic impairment [125,126]. Additionally, chemicals primarily metabolized by the liver are significantly affected by renal impairment. This condition decreases glomerular filtration rate (GFR), renal tubular secretion, and protein binding and has repercussions on intestinal and hepatic [127]. Strougo et al. [124] proposed a semi-PBTK model for patients with both hepatic and renal insufficiencies to accommodate situations where data availability is limited since full PBTK models necessitate an extensive. Li et al. [128] used the PBTK model to predict the intravenous and oral pharmacokinetics of bisoprolol across multiple dose levels in normal adult human populations and patients with renal impairment. During the model construction process, parameters such as GFR, creatinine clearance rate, age, and tissue blood flow were adjusted to represent patients with impaired renal function.

3.4. Metabolic Characteristics and Mechanisms

Complex models can be built based on the ADME processes of chemicals to simulate the effects of transporters, metabolic enzymes, and receptors, thereby predicting the internal dose of specific organs or tissues. Through the integration of TK and toxicodynamic (TD) models, a concentration–effect relationship can be established to aid in investigating the mechanism underlying xenobiotic toxicity. For example, glycyrrhetinic acid (GA) inhibits renal 11β-hydroxysteroid dehydrogenase 2 (11β-HSD 2) and enhances mineralocorticoid receptor activity, leading to pseudo-hyperaldosteronism. Xu et al. [129] developed a human PBTK model of GA and a renal TK/TD model that can describe the effects of 11β-HSD 2 and mineralocorticoid receptors. The biomarkers such as urine cortisol and cortisone concentrations, as well as endpoints including angiotensin II, aldosterone, potassium ions, and sodium ions of GA were predicted using this model. This integrated model can be utilized to determine the safe dosage of GA.

3.5. Mixture Risk Assessment

In daily life, exposure to multiple chemicals simultaneously is a common occurrence. Synergistic and antagonistic interactions between these chemicals can either enhance or weaken their toxic effects. Therefore, accurately assessing the risk of chemical mixtures requires considering these interactions rather than simply calculating the total toxic effects through addition or subtraction. The PBTK model can account for the pharmacokinetic characteristics of chemicals and accurately evaluate the toxic effects of mixtures by incorporating mixture interactions. In studies involving multiple mixtures, it is common practice to establish a binary mixture PBTK model and subsequently build upon it to develop a multivariate mixture PBTK model. For instance, Haddad et al. [130] constructed a PBTK model to simulate the toxicokinetics of mixtures composed of benzene (B), toluene (T), ethylbenzene (E), and m-xylene (X) in rats. Initially, pairwise chemical combinations were simulated with optimized parameters obtained from experimental data for B-T, B-E, B-X, T-E, T-X, and E-X. The mechanism of interactions between binary mixtures was determined by evaluating the metabolic inhibition constant Ki for competitive or noncompetitive inhibition scenarios. Subsequently, using this information as a foundation, they developed a comprehensive PBTK model for all four chemicals in the mixture, which effectively predicted their respective toxicokinetic processes.

4. Prospects

The PBTK model has undergone gradual development since its initial proposal by Teorell in the 1930s [131], particularly from the 1960s onwards. In the 1980s, it found application in toxicology and became known as the PBTK model. With advancements in computer science and technology, the PBPK/PBTK model has experienced rapid growth and widespread utilization. Its applications encompass chemical risk assessment, new drug development, and generic drug research. Compared to traditional compartmental models, the PBTK model offers enhanced accuracy in dose determination, exposure route analysis, and interspecific extrapolation, rendering it more suitable for health risk assessment purposes. A robust PBTK model should possess four essential characteristics [11]: (1) a physiologically-based structure with appropriate assumptions, (2) accurate mathematical equations representing the model, (3) precise parameter values within the model, and (4) comprehensive evaluation and verification of the model’s performance. However, meeting these characteristics consistently throughout the construction processes of a PBTK model poses inherent challenges:
(1)
The PBTK model is based on the understanding of the in vivo ADME processes of chemicals. However, it is challenging to construct a corresponding PBTK model for certain drugs, such as Class 3 and 4 drugs in the Biopharmaceutics Classification System (BCS), traditional Chinese medicine ingredients, heavy metals, etc., due to limited knowledge about their ADME processes.
(2)
The parameters required by the PBTK model are typically obtained through in vivo, in vitro, and in silico assays. However, uncertainties may exist regarding parameters measured using in vitro and in silico assays that require validation with further in vivo experimental data. Moreover, special population groups, such as different races, children, pregnant women, obese individuals, and those who are ill, exhibit distinct physiological and ADME characteristics. Therefore, constructing PBTK models for these special groups necessitates extensive experimental research support.
(3)
Validating exposure characteristics among different subjects under various exposure conditions is crucial for establishing a reliable PBTK model. Unfortunately, in toxicological studies, the lack of relevant experimental data poses challenges for verifying the accuracy of PBTK models.
(4)
Constructing PBTK models can be relatively complex as it requires researchers to possess fundamental knowledge of toxicokinetics, toxicology, physiology, mathematics, and modeling. This complexity limits accessibility to these models.
(5)
For mixtures composed of different types of chemicals, it becomes difficult to construct mixed PBTK models due to variations observed during their respective ADME processes within an organism and also because interaction processes between them can become intricate.
To support scientific researchers and risk assessment practitioners, the US Environmental Protection Agency (EPA) [132] and the World Health Organization (WHO) [80] have issued guidance documents on the PBPK/PBTK model in 2006 and 2010, respectively. These documents provide modeling methods, applications in risk assessment, and standardized model description paradigms to facilitate understanding and qualification by regulatory authorities. Since 2014, the FDA organized multiple PBTK workshops, and modeling reports have been submitted to health authorities to optimize clinical trials or seek clinical study waivers [133]. Recognizing the potential of the PBPK/PBTK model in pharmaceutical development, the FDA and the European Medicines Agency (EMA) released guidance documents in 2016. FDA focused on the standard format and structure of the report [105], and EMA issued the submission of supporting documentation to demonstrate model reliability [106]. The Organization for Economic Cooperation and Development (OECD) has issued a guidance document on PBTK models for regulatory purposes concerning chemicals [134]. The OECD not only provides the process and self-check form for constructing PBTK models but also offers 7541 existing PBTK models for 1150 chemicals and 13 typical cases as references.
Future research in PBTK modeling aims to develop more efficient algorithms that reduce simulation time and user-friendly software with simpler interfaces and lower costs. Integration with quantitative pharmacology, systems biology, cellular pharmacokinetics, and other disciplines will expand the capabilities of the PBTK model. Recognition and acceptance by government and regulatory authorities will play a vital role in promoting its application and development. The advancement and utilization of the PBTK model require interdisciplinary collaboration among physiology, biochemistry, pharmacology, toxicology, mathematics, computer science, and software engineering, as well as involvement from government agencies, academia, research institutions, and enterprises. This collective effort will drive the development and application of the PBTK model, including a more complex mixture of PBTK models, PBTK/TD models, and the transition from animal-based testing methods to non-animal alternatives in health risk assessment. However, when utilizing PBTK models, it is crucial to consider various biological variabilities, differences in exposure scenarios, and the accuracy of model parameters to ensure the reliability and precision of the predictions. It is noteworthy that PBTK models predominantly focus on vertebrates, raising concerns about their applicability to other acquisition of corresponding model parameters when constructing PBTK models for non-vertebrate species within the ecosystem present critical considerations. These issues necessitate attention in PBTK models targeting ecological systems and require resolution in future research.

5. Conclusions

This review systematically delineates the characteristics of PBTK models, providing a comprehensive overview covering aspects of modeling, application, and prospects. The application of PBTK models has not only deepened our understanding of biological processes within the body but also provided substantial support for fields such as drug exploitation and environmental toxicology. The continuous development and refinement of PBTK models will offer more accurate and reliable tools for health risk assessment.

Author Contributions

Conceptualization, T.Z., A.Z. and Q.W.; funding acquisition, A.Z. and Q.W.; writing—original draft, M.C., R.D., C.L., W.B., F.X., S.H., Q.Y. and L.C.; writing—review and editing, A.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (82104520, 82174068), the Natural Science Foundation of Fujian Province (2021J05045), the Fujian Medical University High-level Talent Research Startup Funding Project (XRCZX2020014), and the Open Research Fund of Key Laboratory of Gastrointestinal Cancer (Fujian Medical University), the Ministry of Education (FMUGIC-202101), and the Startup Fund for Scientific Research of Fujian Medical University (2021QH1004).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are included in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dixit, R.; Riviere, J.; Krishnan, K.; Andersen, M.E. Toxicokinetics and physiologically based toxicokinetics in toxicology and risk assessment. J. Toxicol. Environ. Health B Crit. Rev. 2003, 6, 1–40. [Google Scholar] [CrossRef] [PubMed]
  2. Yu, M.; Li, Z. Status and prospects in drug toxicokinetics. Chin. J. New Drugs 2017, 26, 742–748. [Google Scholar]
  3. Lipscomb, J.C.; Haddad, S.; Poet, T.; Krishnan, K. Physiologically-based pharmacokinetic (PBPK) models in toxicity testing and risk assessment. Adv. Exp. Med. Biol. 2012, 745, 76–95. [Google Scholar] [PubMed]
  4. Sarigiannis, D.A.; Karakitsios, S.; Dominguez-Romero, E.; Papadaki, K.; Brochot, C.; Kumar, V.; Schuhmacher, M.; Sy, M.; Mielke, H.; Greiner, M.; et al. Physiology-based toxicokinetic modelling in the frame of the European Human Biomonitoring Initiative. Environ. Res. 2019, 172, 216–230. [Google Scholar] [CrossRef]
  5. Hu, M.; Zhang, Y.; Zhan, M.; He, G.; Qu, W.; Zhou, Y. Physiologically-based toxicokinetic modeling of human dermal exposure to diethyl phthalate: Application to health risk assessment. Chemosphere 2022, 307, 135931. [Google Scholar] [CrossRef]
  6. Yao, M.; Jiang, X.; Lu, Y.; Qiao, Y. PBPK Model and Its Application in TCM Studies. Mod. Tradi. Chin. Med. Mater. Medica-World Sci.Tech. 2007, 3, 55–59. [Google Scholar]
  7. Zhang, S.; Wang, Z.; Chen, J. Application of physiologically based toxicokinetics models in risk assessment of chemicals. Chin. Sci. Bull. 2017, 62, 4139–4150. [Google Scholar] [CrossRef]
  8. Jin, Y.; Ma, Y. Progress in methodology of establishing physiologically based pharmacokinetic models. Acta Pharm. Sinica 2014, 49, 16–22. [Google Scholar]
  9. Niu, Z.; Zang, X. Research advances of the physioloaically based pharmacokinetic(PBPK) models and their applications in health risk assessment. J. Health Saf. Environ. 2015, 15, 30–34. [Google Scholar]
  10. Peters, S.A. Evaluation of a generic physiologically based pharmacokinetic model for lineshape analysis. Clin. Pharmacokinet 2008, 47, 261–275. [Google Scholar] [CrossRef]
  11. Krishnan, K.; Johanson, G. Physiologically-based pharmacokinetic and toxicokinetic models in cancer risk assessment. J. Environ. Sci. Health C Environ. Carcinog. Ecotoxicol. Rev. 2005, 23, 31–53. [Google Scholar] [CrossRef] [PubMed]
  12. Dong, Y.; Zhao, L.; Wu, P.; Wang, J. Characteristics and Research Progress of Physiologically Based Pharmacokinetic Model. Chin. J. ETMF 2012, 18, 247–250. [Google Scholar]
  13. Nestorov, I. Whole body pharmacokinetic models. Clin. Pharmacokinet 2003, 42, 883–908. [Google Scholar] [CrossRef] [PubMed]
  14. Jones, H.M.; Gardner, I.B.; Watson, K.J. Modelling and PBPK simulation in drug discovery. AAPS J. 2009, 11, 155–166. [Google Scholar] [CrossRef] [PubMed]
  15. Krishnan, K.; Loizou, G.D.; Spendiff, M.; Lipscomb, J.C.; Andersen, M.E. Quantitative Modeling in Toxicology, 3rd ed.; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2010; pp. 24–26. [Google Scholar]
  16. Barboriak, D.P.; MacFall, J.R.; Viglianti, B.L.; Dewhirst Dvm, M.W. Comparison of three physiologically-based pharmacokinetic models for the prediction of contrast agent distribution measured by dynamic MR imaging. J. Magn. Reson. Imaging 2008, 27, 1388–1398. [Google Scholar] [CrossRef] [PubMed]
  17. Oliver, R.E.; Jones, A.F.; Rowland, M. A whole-body physiologically based pharmacokinetic model incorporating dispersion concepts: Short and long time characteristics. J. Pharmacokinet. Pharmacodyn. 2001, 28, 27–55. [Google Scholar] [CrossRef]
  18. Pang, K.S.; Durk, M.R. Physiologically-based pharmacokinetic modeling for absorption, transport, metabolism and excretion. J. Pharmacokinet. Pharmacodyn. 2010, 37, 591–615. [Google Scholar] [CrossRef] [PubMed]
  19. Roberts, M.S.; Rowland, M. A dispersion model of hepatic elimination: 1. Formulation of the model and bolus considerations. J. Pharmacokinet. Biopharm. 1986, 14, 227–260. [Google Scholar] [CrossRef]
  20. Roberts, M.S.; Rowland, M. A dispersion model of hepatic elimination: 2. Steady-state considerations--influence of hepatic blood flow, binding within blood, and hepatocellular enzyme activity. J. Pharmacokinet. Biopharm. 1986, 14, 261–288. [Google Scholar] [CrossRef]
  21. Roberts, M.S.; Rowland, M. A dispersion model of hepatic elimination: 3. Application to metabolite formation and elimination kinetics. J. Pharmacokinet. Biopharm. 1986, 14, 289–308. [Google Scholar] [CrossRef]
  22. Brown, R.P.; Delp, M.D.; Lindstedt, S.L.; Rhomberg, L.R.; Beliles, R.P. Physiological parameter values for physiologically based pharmacokinetic models. Toxicol. Ind. Health 1997, 13, 407–484. [Google Scholar] [CrossRef] [PubMed]
  23. Davies, B.; Morris, T. Physiological parameters in laboratory animals and humans. Pharm. Res. 1993, 10, 1093–1095. [Google Scholar] [CrossRef] [PubMed]
  24. Johns, D.O.; Owens, E.O.; Thompson, C.M.; Sonawane, B.; Hattis, D.; Krishnan, K. Physiological Parameters and Databases for PBPK Modeling. In Quantitative Modeling in Toxicology, 3rd ed.; Krishnan, K., Loizou, G.D., Spendiff, M., Lipscomb, J.C., Andersen, M.E., Eds.; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2010; pp. 107–134. [Google Scholar]
  25. Thompson, C.M.; Johns, D.O.; Sonawane, B.; Barton, H.A.; Hattis, D.; Tardif, R.; Krishnan, K. Database for physiologically based pharmacokinetic (PBPK) modeling: Physiological data for healthy and health-impaired elderly. J. Toxicol. Environ. Health B. Crit. Rev. 2009, 12, 1–24. [Google Scholar] [CrossRef] [PubMed]
  26. Caster, W.O.; Poncelet, J.; Simon, A.B.; Armstrong, W.D. Tissue weights of the rat. I. Normal values determined by dissection and chemical methods. Proc. Soc. Exp. Biol. Med. 1956, 91, 122–126. [Google Scholar] [CrossRef] [PubMed]
  27. Domenech, R.J.; Hoffman, J.I.; Noble, M.I.; Saunders, K.B.; Henson, J.R.; Subijanto, S. Total and regional coronary blood flow measured by radioactive microspheres in conscious and anesthetized dogs. Circ. Res. 1969, 25, 581–596. [Google Scholar] [CrossRef] [PubMed]
  28. Mauderly, J.L. Measurement of Respiration and Respiratory Responses During Inhalation Exposures. J. Am. Coll. Toxicol. 1990, 9, 397–405. [Google Scholar] [CrossRef]
  29. Schoeffner, D.J.; Warren, D.A.; Muralidara, S.; Bruckner, J.V.; Simmons, J.E. Organ weights and fat volume in rats as a function of strain and age. J. Toxicol. Environ. Health A 1999, 56, 449–462. [Google Scholar] [CrossRef] [PubMed]
  30. Chen, H.S.; Gross, J.F. Estimation of tissue-to-plasma partition coefficients used in physiological pharmacokinetic models. J. Pharmacokinet. Biopharm. 1979, 7, 117–125. [Google Scholar] [CrossRef]
  31. Dedrick, R.L.; Zaharko, D.S.; Lutz, R.J. Transport and binding of methotrexate in vivo. J. Pharm. Sci. 1973, 62, 882–890. [Google Scholar] [CrossRef]
  32. Gabrielsson, J.L.; Paalzow, L.K.; Nordström, L. A physiologically based pharmacokinetic model for theophylline disposition in the pregnant and nonpregnant rat. J. Pharmacokinet. Biopharm. 1984, 12, 149–165. [Google Scholar] [CrossRef]
  33. Gallo, J.M.; Lam, F.C.; Perrier, D.G. Area method for the estimation of partition coefficients for physiological pharmacokinetic models. J. Pharmacokinet. Biopharm. 1987, 15, 271–280. [Google Scholar] [CrossRef] [PubMed]
  34. Lam, G.; Chen, M.L.; Chiou, W.L. Determination of tissue to blood partition coefficients in physiologically-based pharmacokinetic studies. J. Pharm. Sci. 1982, 71, 454–456. [Google Scholar] [CrossRef] [PubMed]
  35. Bertelsen, S.L.; Hoffman, A.D.; Gallinat, C.A.; Elonen, C.M.; Nichols, J.W. Evaluation of log KOW and tissue lipid content as predictors of chemical partitioning to fish tissues. Environ. Toxicol. Chem. 1998, 17, 1447–1455. [Google Scholar] [CrossRef]
  36. Endo, S.; Bauerfeind, J.; Goss, K.U. Partitioning of neutral organic compounds to structural proteins. Environ. Sci. Technol. 2012, 46, 12697–12703. [Google Scholar] [CrossRef] [PubMed]
  37. Endo, S.; Escher, B.I.; Goss, K.U. Capacities of membrane lipids to accumulate neutral organic chemicals. Environ. Sci. Technol. 2011, 45, 5912–5921. [Google Scholar] [CrossRef] [PubMed]
  38. Endo, S.; Goss, K.U. Serum albumin binding of structurally diverse neutral organic compounds: Data and models. Chem. Res. Toxicol. 2011, 24, 2293–2301. [Google Scholar] [CrossRef] [PubMed]
  39. Fiserova-Bergerova, V.; Diaz, M.L. Determination and prediction of tissue-gas partition coefficients. Int. Arch. Occup. Environ. Health 1986, 58, 75–87. [Google Scholar] [CrossRef] [PubMed]
  40. Gargas, M.L.; Burgess, R.J.; Voisard, D.E.; Cason, G.H.; Andersen, M.E. Partition coefficients of low-molecular-weight volatile chemicals in various liquids and tissues. Toxicol. Appl. Pharmacol. 1989, 98, 87–99. [Google Scholar] [CrossRef]
  41. Gargas, M.L.; Seybold, P.G.; Andersen, M.E. Modeling the tissue solubilities and metabolic rate constant (Vmax) of halogenated methanes, ethanes, and ethylenes. Toxicol. Lett. 1988, 43, 235–256. [Google Scholar] [CrossRef]
  42. Geisler, A.; Endo, S.; Goss, K.U. Partitioning of organic chemicals to storage lipids: Elucidating the dependence on fatty acid composition and temperature. Environ. Sci. Technol. 2012, 46, 9519–9524. [Google Scholar] [CrossRef]
  43. Johanson, G.; Dynésius, B. Liquid/air partition coefficients of six commonly used glycol ethers. Br. J. Ind. Med. 1988, 45, 561–564. [Google Scholar] [CrossRef] [PubMed]
  44. Sato, A.; Nakajima, T. Partition coefficients of some aromatic hydrocarbons and ketones in water, blood and oil. Br. J. Ind. Med. 1979, 36, 231–234. [Google Scholar] [CrossRef] [PubMed]
  45. Igari, Y.; Sugiyama, Y.; Sawada, Y.; Iga, T.; Hanano, M. Prediction of diazepam disposition in the rat and man by a physiologically based pharmacokinetic model. J. Pharmacokinet. Biopharm. 1983, 11, 577–593. [Google Scholar] [CrossRef] [PubMed]
  46. Lin, J.H.; Sugiyama, Y.; Awazu, S.; Hanano, M. In vitro and in vivo evaluation of the tissue-to-blood partition coefficient for physiological pharmacokinetic models. J. Pharmacokinet. Biopharm. 1982, 10, 637–647. [Google Scholar] [CrossRef] [PubMed]
  47. Sultatos, L.G.; Kim, B.; Woods, L. Evaluation of estimations in vitro of tissue/blood distribution coefficients for organothiophosphate insecticides. Toxicol. Appl. Pharmacol. 1990, 103, 52–55. [Google Scholar] [CrossRef] [PubMed]
  48. Assmus, F.; Houston, J.B.; Galetin, A. Incorporation of lysosomal sequestration in the mechanistic model for prediction of tissue distribution of basic drugs. Eur. J. Pharm. Sci. 2017, 109, 419–430. [Google Scholar] [CrossRef] [PubMed]
  49. Gombar, V.K.; Hall, S.D. Quantitative structure-activity relationship models of clinical pharmacokinetics: Clearance and volume of distribution. J. Chem. Inf. Model 2013, 53, 948–957. [Google Scholar] [CrossRef] [PubMed]
  50. Jansson, R.; Bredberg, U.; Ashton, M. Prediction of drug tissue to plasma concentration ratios using a measured volume of distribution in combination with lipophilicity. J. Pharm. Sci. 2008, 97, 2324–2339. [Google Scholar] [CrossRef]
  51. Poulin, P. Drug Distribution to Human Tissues: Prediction and Examination of the Basic Assumption in In Vivo Pharmacokinetics-Pharmacodynamics (PK/PD) Research. J. Pharm. Sci. 2015, 104, 2110–2118. [Google Scholar] [CrossRef]
  52. Poulin, P.; Haddad, S. Advancing prediction of tissue distribution and volume of distribution of highly lipophilic compounds from a simplified tissue-composition-based model as a mechanistic animal alternative method. J. Pharm. Sci. 2012, 101, 2250–2261. [Google Scholar] [CrossRef]
  53. Poulin, P.; Krishnan, K. An algorithm for predicting tissue: Blood partition coefficients of organic chemicals from n-octanol: Water partition coefficient data. J. Toxicol. Environ. Health 1995, 46, 117–129. [Google Scholar] [CrossRef] [PubMed]
  54. Poulin, P.; Krishnan, K. Molecular Structure-Based Prediction of the Partition Coefficients of Organic Chemicals for Physiological Pharmacokinetic Models. Toxicol. Mech. Methods 1996, 6, 117–137. [Google Scholar] [CrossRef]
  55. Poulin, P.; Theil, F.P. Prediction of pharmacokinetics prior to in vivo studies. 1. Mechanism-based prediction of volume of distribution. J. Pharm. Sci. 2002, 91, 129–156. [Google Scholar] [CrossRef] [PubMed]
  56. Poulin, P.; Theil, F.P. A priori prediction of tissue:plasma partition coefficients of drugs to facilitate the use of physiologically-based pharmacokinetic models in drug discovery. J. Pharm. Sci. 2000, 89, 16–35. [Google Scholar] [CrossRef]
  57. 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, 1259–1276. [Google Scholar] [CrossRef] [PubMed]
  58. Rodgers, T.; Rowland, M. Mechanistic approaches to volume of distribution predictions: Understanding the processes. Pharm. Res. 2007, 24, 918–933. [Google Scholar] [CrossRef] [PubMed]
  59. DeJongh, J.; Blaauboer, B.J. Simulation of toluene kinetics in the rat by a physiologically based pharmacokinetic model with application of biotransformation parameters derived independently in vitro and in vivo. Fundam. Appl. Toxicol. 1996, 32, 260–268. [Google Scholar] [CrossRef]
  60. Bauman, J.N.; Kelly, J.M.; Tripathy, S.; Zhao, S.X.; Lam, W.W.; Kalgutkar, A.S.; Obach, R.S. Can in vitro metabolism-dependent covalent binding data distinguish hepatotoxic from nonhepatotoxic drugs? An analysis using human hepatocytes and liver S-9 fraction. Chem. Res. Toxicol. 2009, 22, 332–340. [Google Scholar] [CrossRef]
  61. Iwatsubo, T.; Suzuki, H.; Sugiyama, Y. Prediction of species differences (rats, dogs, humans) in the in vivo metabolic clearance of YM796 by the liver from in vitro data. J. Pharmacol. Exp. Ther. 1997, 283, 462–469. [Google Scholar]
  62. Kedderis, G.L.; Held, S.D. Prediction of furan pharmacokinetics from hepatocyte studies: Comparison of bioactivation and hepatic dosimetry in rats, mice, and humans. Toxicol. Appl. Pharmacol. 1996, 140, 124–130. [Google Scholar] [CrossRef]
  63. Lipscomb, J.C.; Fisher, J.W.; Confer, P.D.; Byczkowski, J.Z. In vitro to in vivo extrapolation for trichloroethylene metabolism in humans. Toxicol. Appl. Pharmacol. 1998, 152, 376–387. [Google Scholar] [CrossRef] [PubMed]
  64. Mortensen, B.; Nilsen, O.G. Allometric species comparison of toluene and n-hexane metabolism: Prediction of hepatic clearance in man from experiments with rodent liver S9 in a head space vial equilibration system. Pharmacol. Toxicol. 1998, 82, 183–188. [Google Scholar] [CrossRef] [PubMed]
  65. Obach, R.S.; Kalgutkar, A.S.; Soglia, J.R.; Zhao, S.X. Can in vitro metabolism-dependent covalent binding data in liver microsomes distinguish hepatotoxic from nonhepatotoxic drugs? An analysis of 18 drugs with consideration of intrinsic clearance and daily dose. Chem. Res. Toxicol. 2008, 21, 1814–1822. [Google Scholar] [CrossRef] [PubMed]
  66. Doddareddy, M.R.; Cho, Y.S.; Koh, H.Y.; Kim, D.H.; Pae, A.N. In silico renal clearance model using classical Volsurf approach. J. Chem. Inf. Model 2006, 46, 1312–1320. [Google Scholar] [CrossRef] [PubMed]
  67. Manga, N.n.; Duffy, J.C.; Rowe, P.H.; Cronin, M.T.D. A Hierarchical QSAR Model for Urinary Excretion of Drugs in Humans as a Predictive Tool for Biotransformation. Mol. Inform. 2003, 22, 263–273. [Google Scholar] [CrossRef]
  68. Paine, S.W.; Barton, P.; Bird, J.; Denton, R.; Menochet, K.; Smith, A.; Tomkinson, N.P.; Chohan, K.K. A rapid computational filter for predicting the rate of human renal clearance. J. Mol. Graph. Model. 2010, 29, 529–537. [Google Scholar] [CrossRef] [PubMed]
  69. Price, K.; Krishnan, K. An integrated QSAR-PBPK modelling approach for predicting the inhalation toxicokinetics of mixtures of volatile organic chemicals in the rat. SAR QSAR Environ. Res. 2011, 22, 107–128. [Google Scholar] [CrossRef] [PubMed]
  70. Nichols, J.W.; Fitzsimmons, P.N.; Burkhard, L.P. In vitro-in vivo extrapolation of quantitative hepatic biotransformation data for fish. II. Modeled effects on chemical bioaccumulation. Environ. Toxicol. Chem. 2007, 26, 1304–1319. [Google Scholar] [CrossRef]
  71. Nichols, J.W.; Schultz, I.R.; Fitzsimmons, P.N. In vitro-in vivo extrapolation of quantitative hepatic biotransformation data for fish. I. A review of methods, and strategies for incorporating intrinsic clearance estimates into chemical kinetic models. Aquat. Toxicol. 2006, 78, 74–90. [Google Scholar] [CrossRef]
  72. Schneckener, S.; Preuss, T.G.; Kuepfer, L.; Witt, J. A workflow to build PBTK models for novel species. Arch. Toxicol. 2020, 94, 3847–3860. [Google Scholar] [CrossRef]
  73. Haddad, S.; Pelekis, M.; Krishnan, K. A methodology for solving physiologically based pharmacokinetic models without the use of simulation softwares. Toxicol. Lett. 1996, 85, 113–126. [Google Scholar] [CrossRef]
  74. Johanson, G.; Näslund, P.H. Spreadsheet programming—A new approach in physiologically based modeling of solvent toxicokinetics. Toxicol. Lett. 1988, 41, 115–127. [Google Scholar] [CrossRef] [PubMed]
  75. Lin, Z.; Jaberi-Douraki, M.; He, C.; Jin, S.; Yang, R.S.H.; Fisher, J.W.; Riviere, J.E. Performance Assessment and Translation of Physiologically Based Pharmacokinetic Models From acslX to Berkeley Madonna, MATLAB, and R Language: Oxytetracycline and Gold Nanoparticles As Case Examples. Toxicol. Sci. 2017, 158, 23–35. [Google Scholar] [CrossRef] [PubMed]
  76. Peters, S.A. Physiologically-based pharmacokinetic (PBPK) modeling and simulations: Principles, methods, and applications in the pharmaceutical industry. In Quantitative Modeling in Toxicology, 3rd ed.; Krishnan, K., Loizou, G.D., Spendiff, M., Lipscomb, J.C., Andersen, M.E., Eds.; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2012; pp. 158–159. [Google Scholar]
  77. Wang, D.; Yang, F.; Sun, L.; Bai, J.; LI, Y. Physiologically based pharmacokinetic modeling based on different softwares. Chin. J. Vet. Med. 2015, 51, 58–60. [Google Scholar]
  78. Mahdi, W.A.; Hussain, A.; Altamimi, M.A.; Alshehri, S.; Bukhari, S.I.; Ahsan, M.N. Experimental Solubility, Thermodynamic/Computational Validations, and GastroPlus-Based In Silico Prediction for Subcutaneous Delivery of Rifampicin. AAPS Pharm. 2021, 22, 116. [Google Scholar] [CrossRef] [PubMed]
  79. Hanke, N.; Gómez-Mantilla, J.D.; Ishiguro, N.; Stopfer, P.; Nock, V. Physiologically Based Pharmacokinetic Modeling of Rosuvastatin to Predict Transporter-Mediated Drug-Drug Interactions. Pharm. Res. 2021, 38, 1645–1661. [Google Scholar] [CrossRef] [PubMed]
  80. WHO. Characterization and Application of Physiologically Based Pharmacokinetic Models in Risk Assessment; WHO: Geneva, Switzerland, 2010; pp. 21–37. [Google Scholar]
  81. Zhao, H.; Ge, P.; Chang, X.; Ren, X.; Chen, L.; Gao, J.; Hu, J.; Li, X. Physiologically based toxicokinetic model for nickel. J. Hyg. Res. 2017, 46, 797–801. [Google Scholar]
  82. Vrieze, S.I. Model selection and psychological theory: A discussion of the differences between the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). Psychol. Methods 2012, 17, 228–243. [Google Scholar] [CrossRef] [PubMed]
  83. Lin, Z.; Li, M.; Gehring, R.; Riviere, J.E. Development and application of a multiroute physiologically based pharmacokinetic model for oxytetracycline in dogs and humans. J. Pharm. Sci. 2015, 104, 233–243. [Google Scholar] [CrossRef]
  84. Portet, S. A primer on model selection using the Akaike Information Criterion. Infect. Dis. Model. 2020, 5, 111–128. [Google Scholar] [CrossRef]
  85. Cox, L.A., Jr. Reassessing benzene risks using internal doses and Monte-Carlo uncertainty analysis. Environ. Health Perspect. 1996, 104, 1413–1429. [Google Scholar] [PubMed]
  86. Farrar, D.; Allen, B.; Crump, K.; Shipp, A. Evaluation of uncertainty in input parameters to pharmacokinetic models and the resulting uncertainty in output. Toxicol. Lett. 1989, 49, 371–385. [Google Scholar] [CrossRef] [PubMed]
  87. Iman, R.L.; Helton, J.C. An Investigation of Uncertainty and Sensitivity Analysis Techniques for Computer Models. Risk Anals. 1988, 8, 71–90. [Google Scholar] [CrossRef]
  88. Khazaee, M.; Ng, C.A. Evaluating parameter availability for physiologically based pharmacokinetic (PBPK) modeling of perfluorooctanoic acid (PFOA) in zebrafish. Environ. Sci. Process Impacts 2018, 20, 105–119. [Google Scholar] [CrossRef] [PubMed]
  89. Krewski, D.; Wang, Y.; Bartlett, S.; Krishnan, K. Uncertainty, variability, and sensitivity analysis in physiological pharmacokinetic models. J. Biopharm. Stat. 1995, 5, 245–271. [Google Scholar] [CrossRef] [PubMed]
  90. Radomyski, A.; Giubilato, E.; Ciffroy, P.; Critto, A.; Brochot, C.; Marcomini, A. Modelling ecological and human exposure to POPs in Venice lagoon—Part II: Quantitative uncertainty and sensitivity analysis in coupled exposure models. Sci. Total Environ. 2016, 569–570, 1635–1649. [Google Scholar] [CrossRef] [PubMed]
  91. Gueorguieva, I.I.; Nestorov, I.A.; Rowland, M. Fuzzy simulation of pharmacokinetic models: Case study of whole body physiologically based model of diazepam. J. Pharmacokinet. Pharmacodyn. 2004, 31, 185–213. [Google Scholar] [CrossRef]
  92. Nestorov, I. Modelling and simulation of variability and uncertainty in toxicokinetics and pharmacokinetics. Toxicol. Lett. 2001, 120, 411–420. [Google Scholar] [CrossRef]
  93. Seng, K.Y.; Nestorov, I.; Vicini, P. Simulating pharmacokinetic and pharmacodynamic fuzzy-parameterized models: A comparison of numerical methods. J. Pharmacokinet. Pharmacodyn. 2007, 34, 595–621. [Google Scholar] [CrossRef]
  94. Nong, A.; Krishnan, K. Estimation of interindividual pharmacokinetic variability factor for inhaled volatile organic chemicals using a probability-bounds approach. Regul. Toxicol. Pharmacol. 2007, 48, 93–101. [Google Scholar] [CrossRef]
  95. Bois, F.Y.; Jamei, M.; Clewell, H.J. PBPK modelling of inter-individual variability in the pharmacokinetics of environmental chemicals. Toxicology 2010, 278, 256–267. [Google Scholar] [CrossRef] [PubMed]
  96. Krauss, M.; Burghaus, R.; Lippert, J.; Niemi, M.; Neuvonen, P.; Schuppert, A.; Willmann, S.; Kuepfer, L.; Görlitz, L. Using Bayesian-PBPK modeling for assessment of inter-individual variability and subgroup stratification. Silico. Pharmacol. 2013, 1, 6. [Google Scholar] [CrossRef] [PubMed]
  97. Krauss, M.; Tappe, K.; Schuppert, A.; Kuepfer, L.; Goerlitz, L. Bayesian Population Physiologically-Based Pharmacokinetic (PBPK) Approach for a Physiologically Realistic Characterization of Interindividual Variability in Clinically Relevant Populations. PLoS ONE 2015, 10, e0139423. [Google Scholar] [CrossRef] [PubMed]
  98. Weijs, L.; Yang, R.S.; Das, K.; Covaci, A.; Blust, R. Application of Bayesian population physiologically based pharmacokinetic (PBPK) modeling and Markov chain Monte Carlo simulations to pesticide kinetics studies in protected marine mammals: DDT, DDE, and DDD in harbor porpoises. Environ. Sci. Technol. 2013, 47, 4365–4374. [Google Scholar] [CrossRef] [PubMed]
  99. Huet, R.; Johanson, G. 1,1-Difluoroethane Detection Time in Blood after Inhalation Abuse Estimated by Monte Carlo PBPK Modeling. Pharmaceutics 2020, 12, 997. [Google Scholar] [CrossRef] [PubMed]
  100. Saltelli, A.; Ratto, M.; Tarantola, S.; Campolongo, F. Update 1 of: Sensitivity analysis for chemical models. Chem. Rev. 2012, 112, Pr1–Pr21. [Google Scholar] [CrossRef]
  101. Brochot, C.; Smith, T.J.; Bois, F.Y. Development of a physiologically based toxicokinetic model for butadiene and four major metabolites in humans: Global sensitivity analysis for experimental design issues. Chem. Biol. Interact. 2007, 167, 168–183. [Google Scholar] [CrossRef] [PubMed]
  102. Saltelli, A.; Tarantola, S.; Chan, K.P.S. A Quantitative Model-Independent Method for Global Sensitivity Analysis of Model Output. Technometrics 1999, 41, 39–56. [Google Scholar] [CrossRef]
  103. Jeong, H.C.; Chae, Y.J.; Shin, K.H. Predicting the systemic exposure and lung concentration of nafamostat using physiologically-based pharmacokinetic modeling. Transl. Clin. Pharmacol. 2022, 30, 201–211. [Google Scholar] [CrossRef]
  104. Yu, L.; Guo, J.; Guo, J.; Yin j Cui, L.; Yang, Y.; Zhao, J.; Peng, S. Prediction of drug toxicity by reverse dosimetry based on PBPK modeling. In Proceedings of the 3rd International Conference on Toxicity Testing Alternative and Translational Toxicology, Nanjing, China, 9–12 July 2017. [Google Scholar]
  105. Louisse, J.; Beekmann, K.; Rietjens, I.M. Use of Physiologically Based Kinetic Modeling-Based Reverse Dosimetry to Predict in Vivo Toxicity from in Vitro Data. Chem. Res. Toxicol. 2017, 30, 114–125. [Google Scholar] [CrossRef]
  106. Cooper, A.B.; Aggarwal, M.; Bartels, M.J.; Morriss, A.; Terry, C.; Lord, G.A.; Gant, T.W. PBTK model for assessment of operator exposure to haloxyfop using human biomonitoring and toxicokinetic data. Regul. Toxicol. Pharmacol. 2019, 102, 1–12. [Google Scholar] [CrossRef] [PubMed]
  107. Poulin, P.; Burczynski, F.J.; Haddad, S. The Role of Extracellular Binding Proteins in the Cellular Uptake of Drugs: Impact on Quantitative In Vitro-to-In Vivo Extrapolations of Toxicity and Efficacy in Physiologically Based Pharmacokinetic-Pharmacodynamic Research. J. Pharm. Sci. 2016, 105, 497–508. [Google Scholar] [CrossRef] [PubMed]
  108. Jaroch, K.; Jaroch, A.; Bojko, B. Cell cultures in drug discovery and development: The need of reliable in vitro-in vivo extrapolation for pharmacodynamics and pharmacokinetics assessment. J. Pharm. Biomed. Anal. 2018, 147, 297–312. [Google Scholar] [CrossRef] [PubMed]
  109. Fabian, E.; Gomes, C.; Birk, B.; Williford, T.; Hernandez, T.R.; Haase, C.; Zbranek, R.; van Ravenzwaay, B.; Landsiedel, R. In vitro-to-in vivo extrapolation (IVIVE) by PBTK modeling for animal-free risk assessment approaches of potential endocrine-disrupting compounds. Arch. Toxicol. 2019, 93, 401–416. [Google Scholar] [CrossRef] [PubMed]
  110. Dejongh, J.; Forsby, A.; Houston, J.B.; Beckman, M.; Combes, R.; Blaauboer, B.J. An Integrated Approach to the Prediction of Systemic Toxicity using Computer-based Biokinetic Models and Biological In vitro Test Methods: Overview of a Prevalidation Study Based on the ECITTS Project. Toxicol. Vitr. 1999, 13, 549–554. [Google Scholar] [CrossRef] [PubMed]
  111. Brinkmann, M.; Eichbaum, K.; Buchinger, S.; Reifferscheid, G.; Bui, T.; Schäffer, A.; Hollert, H.; Preuss, T.G. Understanding receptor-mediated effects in rainbow trout: In vitro-in vivo extrapolation using physiologically based toxicokinetic models. Environ. Sci. Technol. 2014, 48, 3303–3309. [Google Scholar] [CrossRef] [PubMed]
  112. Wang, Y.; Liu, M. Solution of lnternal Doses for lnhaled Vinyl Chloride by Physiologically Based Toxicokinetic (PBTK) Model. Ind. Health Occup. Dis. 2009, 35, 280–284. [Google Scholar]
  113. Gajewska, M.; Worth, A.; Urani, C.; Briesen, H.; Schramm, K.W. Application of physiologically-based toxicokinetic modelling in oral-to-dermal extrapolation of threshold doses of cosmetic ingredients. Toxicol. Lett. 2014, 227, 189–202. [Google Scholar] [CrossRef]
  114. Sweeney, L.M.; Kirman, C.R.; Gargas, M.L.; Carson, M.L.; Tardiff, R.G. Development of a physiologically-based toxicokinetic model of acrylamide and glycidamide in rats and humans. Food Chem. Toxicol. 2010, 48, 668–685. [Google Scholar] [CrossRef]
  115. Abduljalil, K.; Furness, P.; Johnson, T.N.; Rostami-Hodjegan, A.; Soltani, H. Anatomical, physiological and metabolic changes with gestational age during normal pregnancy: A database for parameters required in physiologically based pharmacokinetic modelling. Clin. Pharmacokinet 2012, 51, 365–396. [Google Scholar] [CrossRef]
  116. Xia, B.; Heimbach, T.; Gollen, R.; Nanavati, C.; He, H. A simplified PBPK modeling approach for prediction of pharmacokinetics of four primarily renally excreted and CYP3A metabolized compounds during pregnancy. AAPS J. 2013, 15, 1012–1024. [Google Scholar] [CrossRef] [PubMed]
  117. Wu, Q.; Shi, A.X. Current application of physiologically based pharmacokinetic model and its applications in pregnant women’s drug research. Chin. J. Clin. Pharmacol. 2017, 33, 2209–2211. [Google Scholar]
  118. Dallmann, A.; Pfister, M.; van den Anker, J.; Eissing, T. Physiologically Based Pharmacokinetic Modeling in Pregnancy: A Systematic Review of Published Models. Clin. Pharmacol. Ther. 2018, 104, 1110–1124. [Google Scholar] [CrossRef] [PubMed]
  119. Corley, R.A.; Mast, T.J.; Carney, E.W.; Rogers, J.M.; Daston, G.P. Evaluation of physiologically based models of pregnancy and lactation for their application in children’s health risk assessments. Crit. Rev. Toxicol 2003, 33, 137–211. [Google Scholar] [CrossRef] [PubMed]
  120. Gaohua, L.; Abduljalil, K.; Jamei, M.; Johnson, T.N.; Rostami-Hodjegan, A. A pregnancy physiologically based pharmacokinetic (p-PBPK) model for disposition of drugs metabolized by CYP1A2, CYP2D6 and CYP3A4. Br. J. Clin. Pharmacol. 2012, 74, 873–885. [Google Scholar] [CrossRef] [PubMed]
  121. Edginton, A.N.; Schmitt, W.; Willmann, S. Development and evaluation of a generic physiologically based pharmacokinetic model for children. Clin. Pharmacokinet 2006, 45, 1013–1034. [Google Scholar] [CrossRef] [PubMed]
  122. Kovar, L.; Schräpel, C.; Selzer, D.; Kohl, Y.; Bals, R.; Schwab, M.; Lehr, T. Physiologically-Based Pharmacokinetic (PBPK) Modeling of Buprenorphine in Adults, Children and Preterm Neonates. Pharmaceutics 2020, 12, 578. [Google Scholar] [CrossRef]
  123. Khalil, F.; Läer, S. Physiologically based pharmacokinetic modeling: Methodology, applications, and limitations with a focus on its role in pediatric drug development. J. Biomed. Biotechnol. 2011, 2011, 907461. [Google Scholar] [CrossRef]
  124. Strougo, A.; Yassen, A.; Krauwinkel, W.; Danhof, M.; Freijer, J. A semiphysiological population model for prediction of the pharmacokinetics of drugs under liver and renal disease conditions. Drug Metab. Dispos. 2011, 39, 1278–1287. [Google Scholar] [CrossRef]
  125. Edginton, A.N.; Willmann, S. Physiology-based simulations of a pathological condition: Prediction of pharmacokinetics in patients with liver cirrhosis. Clin. Pharmacokinet 2008, 47, 743–752. [Google Scholar] [CrossRef]
  126. Schuppan, D.; Afdhal, N.H. Liver cirrhosis. Lancet 2008, 371, 838–851. [Google Scholar] [CrossRef] [PubMed]
  127. Nolin, T.D.; Naud, J.; Leblond, F.A.; Pichette, V. Emerging evidence of the impact of kidney disease on drug metabolism and transport. Clin. Pharmacol. Ther. 2008, 83, 898–903. [Google Scholar] [CrossRef] [PubMed]
  128. Li, G.F.; Wang, K.; Chen, R.; Zhao, H.R.; Yang, J.; Zheng, Q.S. Simulation of the pharmacokinetics of bisoprolol in healthy adults and patients with impaired renal function using whole-body physiologically based pharmacokinetic modeling. Acta Pharmacol. Sin. 2012, 33, 1359–1371. [Google Scholar] [CrossRef] [PubMed]
  129. Xu, R.; Liu, X.; Yang, J. A semi-physiologically based pharmacokinetic pharmacodynamic model for glycyrrhizin-induced pseudoaldosteronism and prediction of the dose limit causing hypokalemia in a virtual elderly population. PLoS ONE 2014, 9, e114049. [Google Scholar] [CrossRef] [PubMed]
  130. Haddad, S.; Tardif, R.; Charest-Tardif, G.; Krishnan, K. Physiological modeling of the toxicokinetic interactions in a quaternary mixture of aromatic hydrocarbons. Toxicol. Appl. Pharmacol. 1999, 161, 249–257. [Google Scholar] [CrossRef] [PubMed]
  131. Plowchalk, D.R.; Andersen, M.E.; Bogdanffy, M.S. Physiologically based modeling of vinyl acetate uptake, metabolism, and intracellular pH changes in the rat nasal cavity. Toxicol. Appl. Pharmacol. 1997, 142, 386–400. [Google Scholar] [CrossRef] [PubMed]
  132. McNally, K.; Cotton, R.; Loizou, G.D. A Workflow for Global Sensitivity Analysis of PBPK Models. Front. Pharmacol. 2011, 2, 31. [Google Scholar] [CrossRef] [PubMed]
  133. Lin, W.; Chen, Y.; Unadkat, J.D.; Zhang, X.; Wu, D.; Heimbach, T. Applications, Challenges, and Outlook for PBPK Modeling and Simulation: A Regulatory, Industrial and Academic Perspective. Pharm. Res. 2022, 39, 1701–1731. [Google Scholar] [CrossRef]
  134. OECD. Guidance Document on the Characterization, Validation and Reporting of Physiologically Based Kinetic (PBK) Models for Regulatory Purposes; OECD Series on Testing and Assessment, No. 331, Environment, Health and Safety, Environment Directorate; OECD: Paris, France, 2021. [Google Scholar]
Figure 1. The diagram for the construction of the PBTK model.
Figure 1. The diagram for the construction of the PBTK model.
Toxics 11 00874 g001
Figure 2. The whole PBTK model (A) and semi-PBTK model (B) of humans. IV, intravenous injection; IA, arterial injection; Q, blood flow; M, metabolites; U, urine; CLint, intrinsic clearance; Ke, excretion rate. The subscripts LU, HT, BR, MU, AD, SK, SP, PA, HA, ST, GU, BO, KI, TH, RP, PP, and LI refer to lung, heart, brain, muscle, fat, skin, spleen, pancreas, hepatic artery, stomach, gut, bone, kidney, thymus, richly perfused tissue, poorly perfused tissue, hepatic vein, respectively.
Figure 2. The whole PBTK model (A) and semi-PBTK model (B) of humans. IV, intravenous injection; IA, arterial injection; Q, blood flow; M, metabolites; U, urine; CLint, intrinsic clearance; Ke, excretion rate. The subscripts LU, HT, BR, MU, AD, SK, SP, PA, HA, ST, GU, BO, KI, TH, RP, PP, and LI refer to lung, heart, brain, muscle, fat, skin, spleen, pancreas, hepatic artery, stomach, gut, bone, kidney, thymus, richly perfused tissue, poorly perfused tissue, hepatic vein, respectively.
Toxics 11 00874 g002
Figure 3. In PBTK models, the left panel (A) is superior to the right panel (B) on the consistency between the predicted value (solid line) and the measured value (hollow circle) by visual inspection.
Figure 3. In PBTK models, the left panel (A) is superior to the right panel (B) on the consistency between the predicted value (solid line) and the measured value (hollow circle) by visual inspection.
Toxics 11 00874 g003
Table 1. The mathematical equations in the PBTK model [15].
Table 1. The mathematical equations in the PBTK model [15].
Toxicokinetic ProcessEquation
Absorption
Respiratory tract C a = Q p · C i n h + Q c · C v Q c + Q p / P b (2)
Percutaneous d A s k d t = K p · S C a i r C s k P s : a + Q s k · C a C s k P s : b (3)
Oral d A o d t = K o · A A o (4)
Intravenous C v = K z + t n Q t · C v t Q c (5)
Distribution
Protein binding C b = n · β · K d · C f 1 + K d · C f (6)
Irrigation rate-limiting structure d A t d t = Q t · C a C v t (7)
Membrane rate-limiting tissue d A t d t = P A t · C v t C t P t (8)
Metabolism
First-order kinetic d A m e t d t = K f · C v t · V t = C L · C v t (9)
Second-order kinetic d A m e t d t = K s · C v t · V t · C c f (10)
Saturation process d A m e t d t = V m a x · C v t K m + C v t (11)
Excretion
Kidney d A r c d t = G F R T m K t + C p · C p (12)
Lung C x = 0.7 · C a P b + 0.3 · C i n h (13)
Equation (2): Qc, cardiac output; Qp, alveolar ventilation; Pb, blood-air distribution coefficient; Ca, xenobiotics concentration in arterial blood; Cinh, xenobiotics concentration in inhaled gas; Cv, xenobiotics concentration in mixed venous blood. Equation (3): Ask, total amount of skin xenobiotics exposed; t, time; Kp, skin permeability coefficient; S, skin exposure area; Cair, concentration of xenobiotics in the air; Csk, concentration of exposed skin xenobiotics; Ps:a, skin-air distribution coefficient; Qsk, skin blood flow; Ps:b, skin-blood partition coefficient. Equation (4): Ao, total amount of xenobiotics absorbed; Ko, oral absorption rate constant; A, oral exposure dose. Equation (5): Kz, intravenous administration rate; Qt, blood flow in “t” chamber; Cvt, xenobiotics concentration in venous blood of outflow chamber “t”. Equation (6): Cb, binding xenobiotics concentration; n∙β, maximum binding rate; Kd, dissociation constant; Cf, free xenobiotics concentration. Equation (7): At, total amount of xenobiotics in “t” chamber. Equation (8): PAt, mass transfer coefficient; Pt, organ/tissue–blood allocation coefficient. Equation (9): Amet, total amount of xenobiotics metabolism; Kf, first-order metabolic rate constant; Vt, “t” chamber volume; CL, clearance. Equation (10): Ks, second-order me tabolic rate constant; Ccf, cofactor concentration. Equation (11): Vmax, maximum velocity of enzyme-catalysis; Km, Michaelis constant. Equation (12): Arc, total amount of xenobiotics in kidney; GFR, glomerular filtration rate; Tm, apparent maximum transport rate of the carrier system; Kt, apparent Mieman constant; Cp, xenobiotics concentration in plasma. Equation (13): Cx, xenobiotics concentration in exhaled gas.
Table 2. The physiological, physicochemical, and biochemical parameters in the PBTK model.
Table 2. The physiological, physicochemical, and biochemical parameters in the PBTK model.
ParameterData Source
Physiological parameterBody weight
Cardiac output
Blood flow to organ or tissue
Volume of an organ or tissue
Alveolar ventilation
Literature
In vivo experiment.
Physicochemical parameterBlood-air distribution coefficient
Tissue-blood distribution coefficient
Literature
In vivo experiment
In vitro experiment
In silico prediction
Biochemical parameterMaximum velocity
Michaelis constant
First-order rate constant
Second-order rate constant
Literature
In vivo experiment
In vitro experiment
In silico prediction
For the parameters unable to be obtained, the accessible data can be fixed in the PBTK model, and added data can be measured, thus simulating and obtaining uncertain data.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, M.; Du, R.; Zhang, T.; Li, C.; Bao, W.; Xin, F.; Hou, S.; Yang, Q.; Chen, L.; Wang, Q.; et al. The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment. Toxics 2023, 11, 874. https://doi.org/10.3390/toxics11100874

AMA Style

Chen M, Du R, Zhang T, Li C, Bao W, Xin F, Hou S, Yang Q, Chen L, Wang Q, et al. The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment. Toxics. 2023; 11(10):874. https://doi.org/10.3390/toxics11100874

Chicago/Turabian Style

Chen, Mengting, Ruihu Du, Tao Zhang, Chutao Li, Wenqiang Bao, Fan Xin, Shaozhang Hou, Qiaomei Yang, Li Chen, Qi Wang, and et al. 2023. "The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment" Toxics 11, no. 10: 874. https://doi.org/10.3390/toxics11100874

APA Style

Chen, M., Du, R., Zhang, T., Li, C., Bao, W., Xin, F., Hou, S., Yang, Q., Chen, L., Wang, Q., & Zhu, A. (2023). The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment. Toxics, 11(10), 874. https://doi.org/10.3390/toxics11100874

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