Development, Testing, Parameterisation and Calibration of a Human PBPK Model for the Plasticiser, Di-(2-propylheptyl) Phthalate (DPHP) Using in Silico, in vitro and Human Biomonitoring Data

Human biomonitoring and toxicokinetics as key building blocks for next generation risk assessment

Human health risk assessment is historically built upon animal testing, often following Organisation for Economic Co-operation and Development (OECD) test guidelines and exposure assessments. Using combinations of human relevant in vitro models, chemical analysis and computational (in silico) approaches bring advantages compared to animal studies. These include a greater focus on the human species and on molecular mechanisms and kinetics, identification of Adverse Outcome Pathways and downstream Key Events as well as the possibility of addressing susceptible populations and additional endpoints. Much of the advancement and progress made in the Next Generation Risk Assessment (NGRA) have been primarily focused on new approach methodologies (NAMs) and physiologically based kinetic (PBK) modelling without incorporating human biomonitoring (HBM). The integration of toxicokinetics (TK) and PBK modelling is an essential component of NGRA. PBK models are essential for describing in quantitative terms the TK processes with a focus on the effective dose at the expected target site. Furthermore, the need for PBK models is amplified by the increasing scientific and regulatory interest in aggregate and cumulative exposure as well as interactions of chemicals in mixtures. Since incorporating HBM data strengthens approaches and reduces uncertainties in risk assessment, here we elaborate on the integrated use of TK, PBK modelling and HBM in chemical risk assessment highlighting opportunities as well as challenges and limitations. Examples are provided where HBM and TK/PBK modelling can be used in both exposure assessment and hazard characterization shifting from external exposure and animal dose/response assays to animal-free, internal exposure-based NGRA.

Physiologically based pharmacokinetic modeling to characterize enterohepatic recirculation and predict food effect on the pharmacokinetics of hyzetimibe

BACKGROUND AND OBJECTIVE

Hyzetimibe is a cholesterol absorption inhibitor indicated for the treatment of hypercholesterolemia. This study aims to describe the multiple-peak pharmacokinetics (PK) of hyzetimibe and its active metabolite M1 through physiologically-based pharmacokinetic (PBPK) modeling, and to compare the model predictions of a virtual food effect study with the results of a clinical food effect study.

METHODS

The plasma concentration data used for PBPK modeling were obtained from a single-dose, two-period crossover bioequivalence study in the fasted state. Advanced Compartmental Absorption and Transit model was used for absorption. Enterohepatic recirculation process was modeled by changing the gut physiological state from fasted to fed at meal time. Based on the established PBPK models, a virtual food effect study was simulated. A clinical food effect study was used for model external validation.

RESULTS

PK profiles of hyzetimibe and M1 under fasting condition could be well described by the PBPK model, and the errors of Cmax, AUC0-∞, and AUC0-t were within the two-fold range. Simulated geometric mean ratios (GMRs, fed/fasted) showed that a high-fat breakfast slightly affected the PK of hyzetimibe, expressed as increased Cmax of hyzetimibe (130.6%). Simulated GMRs and 90% confidence intervals of AUC were within the preset bioequivalent range. The results of the simulated virtual food effect trial were consistent with those of the clinical food effect trial.

CONCLUSIONS

The established PBPK model could describe the concentration-time profiles of hyzetimibe and M1 well with good prediction performance. A fully mechanistic model of enterohepatic recirculation warrants further investigation.

The Role of "Physiologically Based Pharmacokinetic Model (PBPK)" New Approach Methodology (NAM) in Pharmaceuticals and Environmental Chemical Risk Assessment

Physiologically Based Pharmacokinetic (PBPK) models are mechanistic tools generally employed in the pharmaceutical industry and environmental health risk assessment. These models are recognized by regulatory authorities for predicting organ concentration-time profiles, pharmacokinetics and daily intake dose of xenobiotics. The extension of PBPK models to capture sensitive populations such as pediatric, geriatric, pregnant females, fetus, etc., and diseased populations such as those with renal impairment, liver cirrhosis, etc., is a must. However, the current modelling practices and existing models are not mature enough to confidently predict the risk in these populations. A multidisciplinary collaboration between clinicians, experimental and modeler scientist is vital to improve the physiology and calculation of biochemical parameters for integrating knowledge and refining existing PBPK models. Specific PBPK covering compartments such as cerebrospinal fluid and the hippocampus are required to gain mechanistic understanding about xenobiotic disposition in these sub-parts. The PBPK model assists in building quantitative adverse outcome pathways (qAOPs) for several endpoints such as developmental neurotoxicity (DNT), hepatotoxicity and cardiotoxicity. Machine learning algorithms can predict physicochemical parameters required to develop in silico models where experimental data are unavailable. Integrating machine learning with PBPK carries the potential to revolutionize the field of drug discovery and development and environmental risk. Overall, this review tried to summarize the recent developments in the in-silico models, building of qAOPs and use of machine learning for improving existing models, along with a regulatory perspective. This review can act as a guide for toxicologists who wish to build their careers in kinetic modeling.

Development, testing, parameterisation, and calibration of a human PBK model for the plasticiser, di (2-ethylhexyl) adipate (DEHA) using in silico, in vitro and human biomonitoring data

Introduction: A physiologically based biokinetic model for di (2-ethylhexyl) adipate (DEHA) based on a refined model for di-(2-propylheptyl) phthalate (DPHP) was developed to interpret the metabolism and biokinetics of DEHA following a single oral dosage of 50 mg to two male and two female volunteers. Methods: The model was parameterized using in vitro and in silico methods such as, measured intrinsic hepatic clearance scaled from in vitro to in vivo and algorithmically predicted parameters such as plasma unbound fraction and tissue:blood partition coefficients (PCs). Calibration of the DEHA model was achieved using concentrations of specific downstream metabolites of DEHA excreted in urine. The total fractions of ingested DEHA eliminated as specific metabolites were estimated and were sufficient for interpreting the human biomonitoring data. Results: The specific metabolites of DEHA, mono-2-ethyl-5-hydroxyhexyl adipate (5OH-MEHA), mono-2-ethyl-5-oxohexyl adipate (5oxo-MEHA), mono-5-carboxy-2-ethylpentyl adipate (5cx-MEPA) only accounted for ∼0.45% of the ingested DEHA. Importantly, the measurements of adipic acid, a non-specific metabolite of DEHA, proved to be important in model calibration. Discussion: The very prominent trends in the urinary excretion of the metabolites, 5cx-MEPA and 5OH-MEHA allowed the important absorption mechanisms of DEHA to be modelled. The model should be useful for the study of exposure to DEHA of the general human population.

Development, testing, parameterisation, and calibration of a human PBPK model for the plasticiser, di-(2-ethylhexyl) terephthalate (DEHTP) using in silico, in vitro and human biomonitoring data

A physiologically based pharmacokinetic model for di-(2-ethylhexyl) terephthalate (DEHTP) based on a refined model for di-(2-propylheptyl) phthalate (DPHP) was developed to interpret the metabolism and biokinetics of DEHTP following a single oral dose of 50 mg to three male volunteers. In vitro and in silico methods were used to generate parameters for the model. For example, measured intrinsic hepatic clearance scaled from in vitro to in vivo and plasma unbound fraction and tissue:blood partition coefficients (PCs) were predicted algorithmically. Whereas the development and calibration of the DPHP model was based upon two data streams, blood concentrations of parent chemical and first metabolite and the urinary excretion of metabolites, the model for DEHTP was calibrated against a single data stream, the urinary excretion of metabolites. Despite the model form and structure being identical significant quantitative differences in lymphatic uptake between the models were observed. In contrast to DPHP the fraction of ingested DEHTP entering lymphatic circulation was much greater and of a similar magnitude to that entering the liver with evidence for the dual uptake mechanisms discernible in the urinary excretion data. Further, the absolute amounts absorbed by the study participants, were much higher for DEHTP relative to DPHP. The in silico algorithm for predicting protein binding performed poorly with an error of more than two orders of magnitude. The extent of plasma protein binding has important implications for the persistence of parent chemical in venous blood-inferences on the behaviour of this class of highly lipophilic chemicals, based on calculations of chemical properties, should be made with extreme caution. Attempting read across for this class of highly lipophilic chemicals should be undertaken with caution since basic adjustments to PCs and metabolism parameters would be insufficient, even when the structure of the model itself is appropriate. Therefore, validation of a model parameterized entirely with in vitro and in silico derived parameters would need to be calibrated against several human biomonitoring data streams to constitute a data rich source chemical to afford confidence for future evaluations of other similar chemicals using the read-across approach.

A review of common non-ortho-phthalate plasticizers for use in food contact materials

Several non-ortho phthalate plasticizers, including ATBC, DEHA, DINCH, DOTP, and ESBO, are currently used in flexible PVC applications for food packaging and processing. The aim of this review is to summarize the available toxicity, migration, and human biomonitoring data. Available assessments from US CPSC, EFSA, other governmental and non-governmental organizations, and published toxicology studies were used to show that these plasticizers are generally well-studied and demonstrate low toxicity with a focus on potential carcinogenicity, reproductive, developmental, and endocrine related adverse effects as well as biodegradation, aquatic toxicity, and bioaccumulation. Seven other plasticizers, 2EHESBO, ASE, CMCHA, DBT, DEHCH, PETV, and TOTM, have at least some recent but limited food contact clearances; assessments from CPSC, EFSA, and robust summaries in the REACH dossiers were reviewed for these products. Data gaps were found for some of these products; however, there were no concerns raised by the existing data, and they for now have limited use in food contact applications. Migration of ASE, COMGHA, DINCH, DOTP, DEHCH, and TOTM in simulants for aqueous and low alcohol foods ranged from <0.02 to 0.165 mg/kg, which showed they are below established migration limits and well-suited for these applications. Human biomonitoring data are available for DINCH, DOTP, DEHA, DINA, and TOTM, and are essential for determining exposure from all uses.

Derivation of a Human In Vivo Benchmark Dose for Bisphenol A from ToxCast In Vitro Concentration Response Data Using a Computational Workflow for Probabilistic Quantitative In Vitro to In Vivo Extrapolation

A computational workflow which integrates physiologically based kinetic (PBK) modelling; global sensitivity analysis (GSA), Approximate Bayesian Computation (ABC), Markov Chain Monte Carlo (MCMC) simulation and the Virtual Cell Based Assay (VCBA) for the estimation of the active, free in vitro concentration of chemical in the reaction medium was developed to facilitate quantitative in vitro to in vivo extrapolation (QIVIVE). The workflow was designed to estimate parameter and model uncertainty within a computationally efficient framework. The workflow was tested using a human PBK model for bisphenol A (BPA) and high throughput screening (HTS) in vitro concentration-response data, for estrogen and pregnane X receptor activation determined in human liver and kidney cell lines, from the ToxCast/Tox21 database. In vivo benchmark dose 10% lower confidence limits (BMDL10) for oral uptake of BPA (ng/kg BW/day) were calculated from the in vivo dose-responses and compared to the human equivalent dose (HED) BMDL10 for relative kidney weight change in the mouse derived by European Food Safety Authority (EFSA). Three from four in vivo BMDL10 values calculated in this study were similar to the EFSA values whereas the fourth was much smaller. The derivation of an uncertainty factor (UF) to accommodate the uncertainties associated with measurements using human cell lines in vitro, extrapolated to in vivo, could be useful for the derivation of Health Based Guidance Values (HBGV).