Human risk assessment of di-isobutyl phthalate through the application of a developed physiologically based pharmacokinetic model of di-isobutyl phthalate and its major metabolite mono-isobutyl phthalate

Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics

Physiologically based pharmacokinetic/toxicokinetic (PBPK/PBTK) models are designed to elucidate the mechanism of chemical compound action in organisms based on the physiological, biochemical, anatomical, and thermodynamic properties of organisms. After nearly a century of research and practice, good results have been achieved in the fields of medicine, environmental science, and ecology. However, there is currently a lack of a more systematic review of progress in the main research directions of PBPK models, especially a more comprehensive understanding of the application in aquatic environmental research. In this review, a total of 3974 articles related to PBPK models from 1996 to 24 March 2024 were collected. Then, the main research areas of the PBPK model were categorized based on the keyword co-occurrence maps and cluster maps obtained by CiteSpace. The results showed that research related to medicine is the main application area of PBPK. Four major research directions included in the medical field were "drug assessment", "cross-species prediction", "drug-drug interactions", and "pediatrics and pregnancy drug development", in which "drug assessment" accounted for 55% of the total publication volume. In addition, bibliometric analyses indicated a rapid growth trend in the application in the field of environmental research, especially in predicting the residual levels in organisms and revealing the relationship between internal and external exposure. Despite facing the limitation of insufficient species-specific parameters, the PBPK model is still an effective tool for improving the understanding of chemical-biological effectiveness and will provide a theoretical basis for accurately assessing potential risks to ecosystems and human health. The combination with the quantitative structure-activity relationship model, Bayesian method, and machine learning technology are potential solutions to the previous research gaps.

Sex differences in 4-tert-octylphenol toxicokinetics: Exploration of sex as an effective covariate through an in vivo modeling approach

4-tert-octylphenol (4-tert-OP) is a potentially harmful substance, which is found widely in the environment. Nevertheless, information on the in vivo toxicokinetics of 4-tert-OP is lacking, and quantitative risk assessment studies are urgently needed. Therefore, we aimed to quantitatively identify differences in the toxicokinetics of 4-tert-OP and its distribution among tissues between sexes. To this end, following exposure of male and female rats to 10 or 50 mg/kg 4-tert-OP orally and 4 or 8 mg/kg 4-tert-OP intravenously, we conducted a quantitative analysis of samples using ultra-high performance liquid chromatography-tandem mass spectrometry. The results revealed that the 4-tert-OP plasma concentration profiles differed between sexes; however, systemic absorption of 4-tert-OP through the gastrointestinal tract occurred within 0.5 h of exposure in both sexes. Although small, the excretion percentage of 4-tert-OP in urine and feces was lower in males than females (0.06-0.08% vs. 0.82-1.11% of exposure). Significant sex differences were also confirmed in the tissue distribution patterns of 4-tert-OP, and overall, the average tissue distribution in males was lower than that in females. The distribution of 4-tert-OP to liver, adipose, spleen, kidney, brain, and lung in both sexes was predominant. A covariate exploration modeling approach revealed that sex explained the differences in 4-tert-OP toxicokinetics between sexes. These significant differences in the toxicokinetics and tissue distribution of 4-tert-OP between sexes will be important for the scientific precision human risk assessment of 4-tert-OP.

Inter-individual exposure variability interpretation through reflection of biological age algorithm in physiologically based toxicokinetic model: Application to human risk assessment of di-isobutyl-phthalate

Age-related changes and interindividual variability in the degree of exposure to hazardous substances in the environment are pertinent factors to be considered in human risk assessment. Existing risk assessments remain in a one-size-fits-all approach, often without due consideration of inter-individual toxicokinetic variability factors, such as age. The purpose of this study was to advance from the existing risk assessment of hazardous substances based on toxicokinetics to a precise human risk assessment by additionally considering the effects of physiologic and metabolic fluctuations and interindividual variability in age. Qualitative age-associated physiologic and metabolic changes in humans, obtained through a meta-analysis, were quantitatively modeled to produce the final biological age algorithm (BAA). The developed BAAs (for males) were extended and applied to the reported testicular reproductive toxicity-focused di-isobutyl-phthalate (DiBP)-mono-isobutyl-phthalate (MiBP) physiologically based toxicokinetic (PBTK) model in males. The advanced PBTK model combined with the BAA was applied to the human risk assessment based on MiBP biomonitoring data. As a result, the specialized DiBP external exposure values for each age could be estimated. Additionally, by applying the Monte Carlo simulation, the distribution of internal exposure diversity among individuals according to the same external exposure dose could be estimated. The contributions of physiologic and metabolic factors to the age-dependent toxicokinetic changes were approximately 93.41-99.99 and 0.01-6.59%, respectively. In addition, the relative contribution of metabolic factors was major in infants and continued to decrease as age increased (up to about age 30 years). This study provides a step-by-step platform that can be widely applied to overcome the limitations of existing toxicokinetic models that still require interindividual pharmacokinetic variability explanations. This will be important for the rationalization and explanation of inter-individual variability in the pharmacokinetics of many substances.

Torsemide Pharmacometrics in Healthy Adult Populations Including CYP2C9 Genetic Polymorphisms and Various Patient Groups through Physiologically Based Pharmacokinetic-Pharmacodynamic Modeling

Torsemide is a widely used diuretic in clinical practice. In this study, pharmacokinetic (PK) and pharmacodynamic (PD) simulations of torsemide for various population groups and exposure scenarios were performed through human-scale physiologically-based PK-PD (PBPK-PD) modeling of torsemide. For PBPK-PD modeling of torsemide, invitro and clinical data of torsemide reported previously were used. After exposure to clinical doses of torsemide, observed plasma (or serum) concentration and urine torsemide excretion profiles were used as PK-data, and observed urinary sodium excretion rate was used as PD-data. The model was then extended to take into account physiological and biochemical factors according to different CYP2C9 phenotypes or patient populations. The established model captured various torsemide clinical results well. Differences in torsemide PKs and PDs between patient groups or CYP2C9 genetic polymorphisms were modelologically identified. It was confirmed that degrees of differences in torsemide PKs and PDs by disease groups were greater than those according to different CYP2C9 phenotypes. According to torsemide administration frequency or dose change, it was confirmed that although the difference in plasma PKs between groups (healthy adult and patient groups) could increase to 14.80 times, the difference in PDs was reduced to 1.01 times. Results of this study suggested that it is very important to consider disease groups in the setting of torsemide clinical therapy and that it is difficult to predict PD proportionally with only differences in PKs of torsemide between population groups. The PBPK-PD model established in this study is expected to be utilized for various clinical cases involving torsemide application in the future, enabling optimal drug therapy.

Physiologically-based toxicokinetic modeling of human dermal exposure to diethyl phthalate: Application to health risk assessment

Diethyl phthalate (DEP) has been most frequently detected in personal care products (PCPs) as a solvent followed by indoor air as one of the semi-volatile organic compounds (SVOCs). Human exposure to DEP predominantly occurs via dermal uptake. However, the available physiologically based toxicokinetics (PBTK) models are developed in rats for risk assessment of DEP exposure resulting from the oral than dermal pathway. To address this issue, DEP in simulated PCPs was dermally administrated to five adult volunteers at real population levels. Following the construction of a dermal absorption model for DEP, the dermal PBTK modeling of DEP involving PCPs and air-to-skin exposure routes in humans was developed for the first time. The data of monoethyl phthalate (MEP) in serum or urine obtained from published human studies and this study were applied to calibrate and validate the developed dermal PBTK model. Monte Carlo simulation was used to evaluate model uncertainty. The dermal absorption fraction of DEP was obtained to be 56.2% for PCPs exposure and 100% for air-to-skin exposure, respectively. Approximate 24.9% of DEP in exposed skin became absorbed into systemic circulation. Model predictions were generally within 2-fold of the observed MEP levels in human serum or urine. Uncertainty analysis showed 90% of the predicted variability (P₉₅/P₅) fell within less than one order of magnitude. Assuming human intake of 5 mg/kg bw per day, the predicted serum area under the curve at steady state of DEP from the dermal route was 1.7 (PCPs) and 2.4 (air) times of those from the peroral route, respectively. It suggested that dermal exposure to DEP would pose greater risk to human health compared with oral exposure. The application of the developed dermal PBTK model provides a valuable insight into health risk assessment of DEP in humans.

Development of physiologically-based toxicokinetic-toxicodynamic (PBTK-TD) model for 4-nonylphenol (4-NP) reflecting physiological changes according to age in males: Application as a new risk assessment tool with a focus on toxicodynamics

Environmental exposure to 4-nonylphenol (4-NP) is extensive, and studies related to human risk assessment must continue. Especially, prediction of toxicodynamics (TDs) related to reproductive toxicity in males is very important in risk-level assessment and management of 4-NP. This study aimed to develop a physiologically-based-toxicokinetic-toxicodynamic (PBTK-TD) model that added a TD prostate model to the previously reported 4-n-nonylphenol (4-n-NP) physiologically-based-pharmacokinetic (PBPK) model. Modeling was performed under the assumption of similar TKs between 4-n-NP and 4-NP because TK experiments on 4-NP, a random-mixture, are practically difficult. This study was very important to quantitatively predict the TKs and TDs of 4-NP by age at exposure using an advanced PBTK-TD model that reflected physiological-changes according to age. TD-modeling was performed based on the reported toxic effects of 4-NP on RWPE-1 cells, a human-prostate-epithelial-cell-line. Through a meta-analysis of reported human physiological data, body weight, tissue volume, and blood flow rate patterns according to age were mathematically modeled. These relationships were reflected in the PBTK-TD model for 4-NP so that the 4-NP TK and TD changes according to age and their differences could be confirmed. Differences in TK and TD parameters of 4-NP at various ages were not large, within 3.61-fold. Point-of-departure (POD) and reference-doses for each age estimated using the model varied as 426.37-795.24 and 42.64-79.52 μg/kg/day, but the differences (in POD or reference doses between ages) were not large, at less than 1.87-times. The PBTK-TD model simulation predicted that even a 100-fold 4-NP POD_man dose would not have large toxicity to the prostate. With a focus on TDs, the predicted maximum possible exposure of 4-NP was as high as 6.06-23.60 mg/kg/day. Several toxicity-related values estimated by the dose-response curve were higher than those calculated, depending upon the PK or TK, which would be useful as a new exposure limit for prostate toxicity of 4-NP.

Human risk assessment of 4-n-nonylphenol (4-n-NP) using physiologically based pharmacokinetic (PBPK) modeling: analysis of gender exposure differences and application to exposure analysis related to large exposure variability in population

As a toxic substance, 4-n-nonylphenol (4-n-NP) or 4-nonylphenol (4-NP) is widely present in the environment. 4-n-NP is a single substance with a linear-alkyl side chain, but 4-NP usually refers to a random mixture containing various branched types. Unfortunately, human risk assessment and/or exposure level analysis for 4-n-NP (or 4-NP) were almost nonexistent, and related research was urgently needed. This study aimed to analyze the various exposures of 4-n-NP (or 4-NP) through development of a physiologically based-pharmacokinetic (PBPK) model considering gender difference in pharmacokinetics of 4-n-NP and its application to human risk assessment studies. A PBPK model was newly developed considering gender differences in 4-n-NP pharmacokinetics and applied to a human risk assessment for each gender. Exposure analysis was performed using a PBPK model that considered gender differences in 4-n-NP (or 4-NP) exposure and high variabilities in several countries. Furthermore, an extended application was attempted as a human risk assessment for random mixture 4-NP, which is difficult to accurately evaluate in reality. External-exposure and margin-of-safety estimated with the same internal exposure amount differed between genders, meaning the need for a differentiated risk assessment considering gender. Exposure analysis based on biomonitoring data confirmed large variability in exposure to 4-n-NP (or 4-NP) by country, group, and period. External-exposures estimated using PBPK model varied widely, ranging from 0.039 to 63.875 mg/kg/day (for 4-n-NP or 4-NP). By country, 4-n-NP (or 4-NP) exposure was higher in females than in males and the margin-of-safety tended to be low. Overall, exposure to 4-n-NP (or 4-NP) in populations was largely not safe, suggesting need for ongoing management and monitoring. Considering low in vivo accumulation confirmed by PBPK model, risk reduction of 4-n-NP is possible by reducing its use.

Toxicokinetics of di-isodecyl phthalate and its major metabolites in rats through the application of a developed and validated UHPLC-ESI-MS/MS method

Di-isodecyl phthalate (DiDP) is a high-molecular-weight phthalate that is mainly used as a plasticizer for plastics. Therefore, exposure to DiDP in the environment has become common with the increasing use of plastics around the world. Environmental regulations and scientific risk management for DiDP, which can be associated with endocrine disruption and various metabolic diseases, are urgently needed. The purpose of this study was to provide useful reference material for future human DiDP risk assessments by conducting toxicokinetic studies on DiDP. Rats were given 100 mg/kg of DiDP orally or intravenously, and plasma, urine, feces, and various tissues were sampled at preset times. DiDP and its major metabolites mono-isodecyl-phthalate (MiDP), mono-hydroxy-isodecyl-phthalate (MHiDP), mono-carboxy-isononyl-phthalate (MCiNP), and mono-oxo-isodecyl-phthalate (MOiDP) were simultaneously quantified from collected biological samples through the application of a newly developed and verified ultrahigh-performance liquid chromatography-electrospray ionization-tandem mass spectrometer (UHPLC-ESI-MS/MS) method. Based on the quantitative results for each analyte, toxicokinetic analyses were performed. DiDP was rapidly and extensively metabolized to MiDP, MHiDP, MCiNP, and MOiDP. The major metabolite excreted in the urine was MCiNP, suggesting that it could be a useful biomarker. The conjugated forms of DiDP and its metabolites have been significantly quantified in the plasma, urine, and feces. DiDP and its major metabolites were also distributed in various tissues in significant quantities. The toxicokinetic properties of DiDP, which have not been clearly reported previously, were identified through this study. This report will serve as a useful reference for future DiDP environmental regulation and scientific human risk assessment studies.