Developing a Physiologically-Based Pharmacokinetic Model Knowledgebase in Support of Provisional Model Construction

Physiologically based kinetic (PBK) modeling of propiconazole using a machine learning-enhanced read-across approach for interspecies extrapolation

A significant challenge in the traditional human health risk assessment of agrochemicals is the uncertainty in quantifying the interspecies differences between animal models and humans. To work toward a more accurate and animal-free risk determination, new approaches such as physiologically based kinetic (PBK) modeling have been used to perform dosimetry extrapolation from animals to humans. However, the regulatory use and acceptance of PBK modeling is limited for chemicals that lack in vivo animal pharmacokinetic (PK) data, given the inability to evaluate models. To address these challenges, this study developed PBK models in the absence of in vivo PK data for the fungicide propiconazole, an activator of constitutive androstane receptor (CAR)/pregnane X receptor (PXR). A fit-for-purpose read-across approach was integrated with hierarchical clustering - an unsupervised machine learning algorithm, to bridge the knowledge gap. The integration allowed the incorporation of a broad spectrum of attributes for analog consideration, and enabled the analog selection in a simple, reproducible, and objective manner. The applicability was evaluated and demonstrated using penconazole (source) and three pseudo-unknown target chemicals (epoxiconazole, tebuconazole and triadimefon). Applying this machine learning-enhanced read-across approach, difenoconazole was selected as the most appropriate analog for propiconazole. A mouse PBK model was developed and evaluated for difenoconazole (source), with the mode of action of CAR/PXR activation incorporated to simulate the in vivo autoinduction of metabolism. The difenoconazole mouse model then served as a template for constructing the propiconazole mouse model. A parallelogram approach was subsequently applied to develop the propiconazole rat and human models, enabling a quantitative assessment of interspecies differences in dosimetry. This integrated approach represents a substantial advancement toward refining risk assessment of propiconazole within the framework of animal alternative safety assessment strategies.

A physiologically based pharmacokinetic/pharmacodynamic model to determine dosage regimens and withdrawal intervals of aditoprim against Streptococcus suis

Introduction: Streptococcus suis (S. suis) is a zoonotic pathogen threatening public health. Aditoprim (ADP), a novel veterinary medicine, exhibits an antibacterial effect against S. suis. In this study, a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model was used to determine the dosage regimens of ADP against S. suis and withdrawal intervals. Methods: The PBPK model of ADP injection can predict drug concentrations in plasma, liver, kidney, muscle, and fat. A semi-mechanistic pharmacodynamic (PD) model, including susceptible subpopulation and resistant subpopulation, is successfully developed by a nonlinear mixed-effect model to evaluate antibacterial effects. An integrated PBPK/PD model is conducted to predict the time-course of bacterial count change and resistance development under different ADP dosages. Results: ADP injection, administrated at 20 mg/kg with 12 intervals for 3 consecutive days, can exert an excellent antibacterial effect while avoiding resistance emergence. The withdrawal interval at the recommended dosage regimen is determined as 18 days to ensure food safety. Discussion: This study suggests that the PBPK/PD model can be applied as an effective tool for the antibacterial effect and safety evaluation of novel veterinary drugs.

Case study on the impact of the source of metabolism parameters in next generation physiologically based pharmacokinetic models: Implications for occupational exposures to trimethylbenzenes

Physiologically based pharmacokinetic (PBPK) models are a means of making important linkages between exposure assessment and in vitro toxicity. A key constraint on rapid application of PBPK models in risk assessment is traditional reliance on substance-specific in vivo toxicokinetic data to evaluate model quality. Bounding conditions, in silico, in vitro, and chemical read-across approaches have been proposed as alternative sources for metabolic clearance estimates. A case study to test consistency of predictive ability across these approaches was conducted using trimethylbenzenes (TMB) as prototype chemicals. Substantial concordance was found among TMB isomers with respect to accuracy (or inaccuracy) of approaches to estimating metabolism; for example, the bounding conditions never reproduced the human in vivo toxicokinetic data within two-fold. Using only approaches that gave acceptable prediction of in vivo toxicokinetics for the source compound (1,2,4-TMB) substantially narrowed the range of plausible internal doses for a given external dose for occupational, emergency response, and environmental/community health risk assessment scenarios for TMB isomers. Thus, risk assessments developed using the target compound models with a constrained subset of metabolism estimates (determined for source chemical models) can be used with greater confidence that internal dosimetry will be estimated with accuracy sufficient for the purpose at hand.

IVIVE: Facilitating the Use of In Vitro Toxicity Data in Risk Assessment and Decision Making

During the past few decades, the science of toxicology has been undergoing a transformation from observational to predictive science. New approach methodologies (NAMs), including in vitro assays, in silico models, read-across, and in vitro to in vivo extrapolation (IVIVE), are being developed to reduce, refine, or replace whole animal testing, encouraging the judicious use of time and resources. Some of these methods have advanced past the exploratory research stage and are beginning to gain acceptance for the risk assessment of chemicals. A review of the recent literature reveals a burst of IVIVE publications over the past decade. In this review, we propose operational definitions for IVIVE, present literature examples for several common toxicity endpoints, and highlight their implications in decision-making processes across various federal agencies, as well as international organizations, including those in the European Union (EU). The current challenges and future needs are also summarized for IVIVE. In addition to refining and reducing the number of animals in traditional toxicity testing protocols and being used for prioritizing chemical testing, the goal to use IVIVE to facilitate the replacement of animal models can be achieved through their continued evolution and development, including a strategic plan to qualify IVIVE methods for regulatory acceptance.

In Silico Prediction of Pharmacokinetic Profile for Human Oral Drug Candidates Which Lack Clinical Pharmacokinetic Experiment Data

BACKGROUNDS AND OBJECTIVES

In silico methods which can generate high-quality physiologically based pharmacokinetic (PBPK) models for arbitrary drug candidates are greatly needed to select developable drug candidates that escape drug attrition because of the poor pharmacokinetic profile. The purpose of this study is to develop a novel protocol to preliminarily predict the concentration profile of a target drug based on the PBPK model of a structurally similar template drug by combining two software platforms for PBPK modeling, the SimCYP simulator and ADMET Predictor.

METHODS

The method was evaluated by utilizing 13 drug pairs from 18 drugs in the built-in database of the SimCYP software. All drug pairs have Tanimoto scores (TS) no less than 0.5. As each drug in a drug pair can serve as both target and template, 26 sets were studied in this work. Three versions (V1, V2 and V3) of models for the target drug were constructed by replacing the corresponding parameters of the template drug step by step with those predicted by ADMET Predictor for the target drug. V1 represents the replacement of molecular weight (MW), V2 includes the replacement of parameter MW, fraction unbound in plasma (f_u), blood-to-plasma partition ratio (B/P), logarithm of the octanol-buffer partition coefficient (log P_o:w) and acid dissociation constant (pK_a). In V3, all above-mentioned parameters as well as human jejunum effective permeability (P_eff), V_d and cytochrome P450 (CYP) metabolism parameters (K_m, V_max or CL_int) are modified. Normalized root mean square error (NRMSE) was used for the evaluation of the model performance.

RESULTS

We found that the performance of the three versions of the models depends on structural similarity of the drug pairs. For Group I drug pairs (TS ≤ 0.7), V2 and V3 performed better than V1 in terms of NRMSE; for Group II drug pairs (0.7 < TS ≤ 0.9), 8 out of 10 V3 models had NRMSE < 0.2, the cutoff we applied to judge whether the simulated concentration-time (C-T) curve was satisfactory or not. V3 outperformed the V1 and V2 versions. For the two drug pairs belonging to Group III (TS > 0.9), V2 outperformed V1 and V3, suggesting more unnecessary replacement can lower the performance of PBPK models. We also investigated how the prediction accuracy of ADMET Predictor as well as its collaboration with SimCYP influences the quality of PBPK models constructed using SimCYP.

CONCLUSION

In conclusion, we generated practical guidance on applying two mainstream software packages, ADMET Predictor and SimCYP, to construct PBPK models for drugs or drug candidates that lack ADME parameters in model construction.

A Systematic Review of Published Physiologically-based Kinetic Models and an Assessment of their Chemical Space Coverage

Across multiple sectors, including food, cosmetics and pharmaceutical industries, there is a need to predict the potential effects of xenobiotics. These effects are determined by the intrinsic ability of the substance, or its derivatives, to interact with the biological system, and its concentration-time profile at the target site. Physiologically-based kinetic (PBK) models can predict organ-level concentration-time profiles, however, the models are time and resource intensive to generate de novo. Read-across is an approach used to reduce or replace animal testing, wherein information from a data-rich chemical is used to make predictions for a data-poor chemical. The recent increase in published PBK models presents the opportunity to use a read-across approach for PBK modelling, that is, to use PBK model information from one chemical to inform the development or evaluation of a PBK model for a similar chemical. Essential to this process, is identifying the chemicals for which a PBK model already exists. Herein, the results of a systematic review of existing PBK models, compliant with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) format, are presented. Model information, including species, sex, life-stage, route of administration, software platform used and the availability of model equations, was captured for 7541 PBK models. Chemical information (identifiers and physico-chemical properties) has also been recorded for 1150 unique chemicals associated with these models. This PBK model data set has been made readily accessible, as a Microsoft Excel^® spreadsheet, providing a valuable resource for those developing, using or evaluating PBK models in industry, academia and the regulatory sectors.

Development and intercomparison of single and multicompartment physiologically-based toxicokinetic models: Implications for model selection and tiered modeling frameworks

This study describes the development and intercomparison of generic physiologically-based toxicokinetic (PBTK) models for humans comprised of internally consistent one-compartment (1Co-) and multi-compartment (MCo-) implementations (G-PBTK). The G-PBTK models were parameterized for an adult male (70 kg) using common physiological parameters and in vitro biotransformation rate estimates and subsequently evaluated using independent concentration versus time data (n = 6) and total elimination half-lives (n = 15) for diverse organic chemicals. The model performance is acceptable considering the inherent uncertainty in the biotransformation rate data and the absence of model calibration. The G-PBTK model was then applied using hypothetical neutral organics, acidic ionizable organics and basic ionizable organics (IOCs) to identify combinations of partitioning properties and biotransformation rates leading to substantial discrepancies between 1Co- and MCo-PBTK calculations for whole body concentrations and half-lives. The 1Co- and MCo-PBTK model calculations for key toxicokinetic parameters are broadly consistent unless biotransformation is rapid (e.g., half-life less than five days). When half-lives are relatively short, discrepancies are greatest for the neutral organics and least for the acidic IOCs which follows from the estimated volumes of distribution (e.g., VD_SS = 9.6-15.4 L/kg vs 0.3-1.6 L/kg for the neutral and acidic compounds respectively) and the related approach to internal chemical equilibrium. The model intercomparisons demonstrate that 1Co-PBTK models can be applied with confidence to many exposure scenarios, particularly those focused on chronic or repeat exposures and for prioritization and screening-level decision contexts. However, MCo-PBTK models may be necessary in certain contexts, particularly for intermittent, short-term and highly variable exposures. A key recommendation to guide model selection and the development of tiered PBTK modeling frameworks that emerges from this study is the need to harmonize models with respect to parameterization and process descriptions to the greatest extent possible when proceeding from the application of simpler to more complex modeling tools as part of chemical assessment activities.

DES-Tcell is a knowledgebase for exploring immunology-related literature

T-cells are a subtype of white blood cells circulating throughout the body, searching for infected and abnormal cells. They have multifaceted functions that include scanning for and directly killing cells infected with intracellular pathogens, eradicating abnormal cells, orchestrating immune response by activating and helping other immune cells, memorizing encountered pathogens, and providing long-lasting protection upon recurrent infections. However, T-cells are also involved in immune responses that result in organ transplant rejection, autoimmune diseases, and some allergic diseases. To support T-cell research, we developed the DES-Tcell knowledgebase (KB). This KB incorporates text- and data-mined information that can expedite retrieval and exploration of T-cell relevant information from the large volume of published T-cell-related research. This KB enables exploration of data through concepts from 15 topic-specific dictionaries, including immunology-related genes, mutations, pathogens, and pathways. We developed three case studies using DES-Tcell, one of which validates effective retrieval of known associations by DES-Tcell. The second and third case studies focuses on concepts that are common to Grave’s disease (GD) and Hashimoto’s thyroiditis (HT). Several reports have shown that up to 20% of GD patients treated with antithyroid medication develop HT, thus suggesting a possible conversion or shift from GD to HT disease. DES-Tcell found miR-4442 links to both GD and HT, and that miR-4442 possibly targets the autoimmune disease risk factor CD6, which provides potential new knowledge derived through the use of DES-Tcell. According to our understanding, DES-Tcell is the first KB dedicated to exploring T-cell-relevant information via literature-mining, data-mining, and topic-specific dictionaries.

High-throughput PBTK models for in vitro to in vivo extrapolation

INTRODUCTION

Toxicity data are unavailable for many thousands of chemicals in commerce and the environment. Therefore, risk assessors need to rapidly screen these chemicals for potential risk to public health. High-throughput screening (HTS) for in vitro bioactivity, when used with high-throughput toxicokinetic (HTTK) data and models, allows characterization of these thousands of chemicals.

AREAS COVERED

This review covers generic physiologically based toxicokinetic (PBTK) models and high-throughput PBTK modeling for in vitro-in vivo extrapolation (IVIVE) of HTS data. We focus on 'httk', a public, open-source set of computational modeling tools and in vitro toxicokinetic (TK) data.

EXPERT OPINION

HTTK benefits chemical risk assessors with its ability to support rapid chemical screening/prioritization, perform IVIVE, and provide provisional TK modeling for large numbers of chemicals using only limited chemical-specific data. Although generic TK model design can increase prediction uncertainty, these models provide offsetting benefits by increasing model implementation accuracy. Also, public distribution of the models and data enhances reproducibility. For the httk package, the modular and open-source design can enable the tool to be used and continuously improved by a broad user community in support of the critical need for high-throughput chemical prioritization and rapid dose estimation to facilitate rapid hazard assessments.

Integration of data across toxicity endpoints for improved safety assessment of chemicals: the example of carcinogenicity assessment

In view of the need to enhance the assessment of consumer products called for in the EU Chemicals Strategy for Sustainability, we developed a methodology for evaluating hazard by combining information across different systemic toxicity endpoints and integrating the information with new approach methodologies. This integrates mechanistic information with a view to avoiding redundant in vivo studies, minimising reliance on apical endpoint tests and ultimately devising efficient testing strategies. Here, we present the application of our methodology to carcinogenicity assessment, mapping the available information from toxicity test methods across endpoints to the key characteristics of carcinogens. Test methods are deconstructed to allow the information they provide to be organised in a systematic way, enabling the description of the toxicity mechanisms leading to the adverse outcome. This integrated approach provides a flexible and resource-efficient means of fully exploiting test methods for which test guidelines are available to fulfil regulatory requirements for systemic toxicity assessment as well as identifying where new methods can be integrated.

Assessment of the predictive capacity of a physiologically based kinetic model using a read-across approach

•

Potential regulatory application of PBK modelling information to assist read-across.

•

Presents workflow to read across PBK model information from data-rich to data-poor chemicals.

•

Describes appropriate analogue selection based on a set of specific criteria.

•

Uses estragole and safrole as source chemicals for a target chemical - methyleugenol.

•

Example of PBK model validation where in vivo kinetic data are lacking.

With current progress in science, there is growing interest in developing and applying Physiologically Based Kinetic (PBK) models in chemical risk assessment, as knowledge of internal exposure to chemicals is critical to understanding potential effects in vivo. In particular, a new generation of PBK models is being developed in which the model parameters are derived from in silico and in vitro methods. To increase the acceptance and use of these “Next Generation PBK models”, there is a need to demonstrate their validity. However, this is challenging in the case of data-poor chemicals that are lacking in kinetic data and for which predictive capacity cannot, therefore, be assessed. The aim of this work is to lay down the fundamental steps in using a read across framework to inform modellers and risk assessors on how to develop, or evaluate, PBK models for chemicals without in vivo kinetic data. The application of a PBK model that takes into account the absorption, distribution, metabolism and excretion characteristics of the chemical reduces the uncertainties in the biokinetics and biotransformation of the chemical of interest. A strategic flow-charting application, proposed herein, allows users to identify the minimum information to perform a read-across from a data-rich chemical to its data-poor analogue(s). The workflow analysis is illustrated by means of a real case study using the alkenylbenzene class of chemicals, showing the reliability and potential of this approach. It was demonstrated that a consistent quantitative relationship between model simulations could be achieved using models for estragole and safrole (source chemicals) when applied to methyleugenol (target chemical). When the PBK model code for the source chemicals was adapted to utilise input values relevant to the target chemical, simulation was consistent between the models. The resulting PBK model for methyleugenol was further evaluated by comparing the results to an existing, published model for methyleugenol, providing further evidence that the approach was successful. This can be considered as a “read-across” approach, enabling a valid PBK model to be derived to aid the assessment of a data poor chemical.

Modeling Pharmacokinetic Natural Product-Drug Interactions for Decision-Making: A NaPDI Center Recommended Approach

The popularity of botanical and other purported medicinal natural products (NPs) continues to grow, especially among patients with chronic illnesses and patients managed on complex prescription drug regimens. With few exceptions, the risk of a given NP to precipitate a clinically significant pharmacokinetic NP-drug interaction (NPDI) remains understudied or unknown. Application of static or dynamic mathematical models to predict and/or simulate NPDIs can provide critical information about the potential clinical significance of these complex interactions. However, methods used to conduct such predictions or simulations are highly variable. Additionally, published reports using mathematical models to interrogate NPDIs are not always sufficiently detailed to ensure reproducibility. Consequently, guidelines are needed to inform the conduct and reporting of these modeling efforts. This recommended approach from the Center of Excellence for Natural Product Drug Interaction Research describes a systematic method for using mathematical models to interpret the interaction risk of NPs as precipitants of potential clinically significant pharmacokinetic NPDIs. A framework for developing and applying pharmacokinetic NPDI models is presented with the aim of promoting accuracy, reproducibility, and generalizability in the literature. SIGNIFICANCE STATEMENT: Many natural products (NPs) contain phytoconstituents that can increase or decrease systemic or tissue exposure to, and potentially the efficacy of, a pharmaceutical drug; however, no regulatory agency guidelines exist to assist in predicting the risk of these complex interactions. This recommended approach from a multi-institutional consortium designated by National Institutes of Health as the Center of Excellence for Natural Product Drug Interaction Research provides a framework for modeling pharmacokinetic NP-drug interactions.

Development of a physiologically-based pharmacokinetic pediatric brain model for prediction of cerebrospinal fluid drug concentrations and the influence of meningitis

Different pediatric physiologically-based pharmacokinetic (PBPK) models have been described incorporating developmental changes that influence plasma drug concentrations. Drug disposition into cerebrospinal fluid (CSF) is also subject to age-related variation and can be further influenced by brain diseases affecting blood-brain barrier integrity, like meningitis. Here, we developed a generic pediatric brain PBPK model to predict CSF concentrations of drugs that undergo passive transfer, including age-appropriate parameters. The model was validated for the analgesics paracetamol, ibuprofen, flurbiprofen and naproxen, and for a pediatric meningitis population by empirical optimization of the blood-brain barrier penetration of the antibiotic meropenem. Plasma and CSF drug concentrations derived from the literature were used to perform visual predictive checks and to calculate ratios between simulated and observed area under the concentration curves (AUCs) in order to evaluate model performance. Model-simulated concentrations were comparable to observed data over a broad age range (3 months–15 years postnatal age) for all drugs investigated. The ratios between observed and simulated AUCs (AUCo/AUCp) were within 2-fold difference both in plasma (range 0.92–1.09) and in CSF (range 0.64–1.23) indicating acceptable model performance. The model was also able to describe disease-mediated changes in neonates and young children (<3m postnatal age) related to meningitis and sepsis (range AUCo/AUCp plasma: 1.64–1.66, range AUCo/AUCp CSF: 1.43–1.73). Our model provides a new computational tool to predict CSF drug concentrations in children with and without meningitis and can be used as a template model for other compounds that passively enter the CNS.

Developmental processes in children affect pharmacokinetics and should ideally be taken into account when establishing drug dosing regimens. One way to incorporate developmental differences is by making use of physiologically-based pharmacokinetic (PBPK) models in which kinetic equations are used to describe drug disposition processes and developmental biology. With these equations the absorption of drugs into the model, the flow of drugs between different compartments (representing major organs/tissues), and excretion from the model are predicted. PBPK models can also be used to describe drug concentrations in different target tissues, which often correlate better with the clinical effects. Here, we developed a generic pediatric PBPK model of drug disposition in the cerebrospinal fluid (CSF), that was able to describe clinically measured drug concentrations of several drugs in neonates and children. The model could be useful in predicting CSF concentrations of other drugs in pediatric populations where clinical data is often sparse or absent and by this means guide first-in-child dose recommendations.

Role of Physiologically Based Kinetic modelling in addressing environmental chemical mixtures - A review

•

The availability and applicability of Physiologically Based Kinetic (PBK) models for mixtures is reviewed.

•

PBK models can support risk assessment of mixtures by incorporating the toxicokinetic processes.

•

Quantitative structure-activity relationship (QSAR) models can be used to fill data gaps in PBK modelling.

•

PBK models for mixtures can be improved by including various types of interactions.

The role of Physiologically Based Kinetic (PBK) modelling in assessing mixture toxicology has been growing for the last three decades. It has been widely used to investigate and address interactions in mixtures. This review describes the current state-of-the-art of PBK models for chemical mixtures and to evaluate the applications of PBK modelling for mixtures with emphasis on their role in chemical risk assessment. A total of 35 mixture PBK models were included after searching web resources (Scopus, PubMed, Web of Science, and Google Scholar), screening for duplicates, and excluding articles based on eligibility criteria. Binary mixtures and volatile organic compounds accounted for two-thirds of the chemical mixtures identified. The most common exposure route and modelled system were found to be inhalation and rats respectively. Twenty two (22) models were for binary mixtures, 5 for ternary mixtures, 3 for quaternary mixtures, and 5 for complex mixtures. Both bottom-up and top-down PBK modelling approaches are described. Whereas bottom-up approaches are based on a series of binary interactions, top-down approaches are based on the lumping of mixture components. Competitive inhibition is the most common type of interaction among the various types of mixtures, and usually becomes a concern at concentrations higher than environmental exposure levels. It leads to reduced biotransformation that either means a decrease in the amount of toxic metabolite formation or an increase in toxic parent chemical accumulation. The consequence is either lower or higher toxicity compared to that estimated for the mixture based on the additivity principle. Therefore, PBK modelling can play a central role in predicting interactions in chemical mixture risk assessment.

Challenges in working towards an internal threshold of toxicological concern (iTTC) for use in the safety assessment of cosmetics: Discussions from the Cosmetics Europe iTTC Working Group workshop

The Threshold of Toxicological Concern (TTC) is an important risk assessment tool which establishes acceptable low-level exposure values to be applied to chemicals with limited toxicological data. One of the logical next steps in the continued evolution of TTC is to develop this concept further so that it is representative of internal exposures (TTC based on plasma concentration). An internal TTC (iTTC) would provide threshold values that could be utilized in exposure-based safety assessments. As part of a Cosmetics Europe (CosEu) research program, CosEu has initiated a project that is working towards the development of iTTCs that can be used for the human safety assessment. Knowing that the development of an iTTC is an ambitious and broad-spanning topic, CosEu organized a Working Group comprised a balance of multiple stakeholders (cosmetics and chemical industries, the EPA and JRC and academia) with relevant experience and expertise and workshop to critically evaluate the requirements to establish an iTTC. Outcomes from the workshop included an evaluation on the current state of the science for iTTC, the overall iTTC strategy, selection of chemical databases, capture and curation of chemical information, ADME and repeat dose data, expected challenges, as well as next steps and ongoing work.

Challenges Associated With Applying Physiologically Based Pharmacokinetic Modeling for Public Health Decision-Making

The development and application of physiologically based pharmacokinetic (PBPK) models in chemical toxicology have grown steadily since their emergence in the 1980s. However, critical evaluation of PBPK models to support public health decision-making across federal agencies has thus far occurred for only a few environmental chemicals. In order to encourage decision-makers to embrace the critical role of PBPK modeling in risk assessment, several important challenges require immediate attention from the modeling community. The objective of this contemporary review is to highlight 3 of these challenges, including: (1) difficulties in recruiting peer reviewers with appropriate modeling expertise and experience; (2) lack of confidence in PBPK models for which no tissue/plasma concentration data exist for model evaluation; and (3) lack of transferability across modeling platforms. Several recommendations for addressing these 3 issues are provided to initiate dialog among members of the PBPK modeling community, as these issues must be overcome for the field of PBPK modeling to advance and for PBPK models to be more routinely applied in support of public health decision-making.

Next generation physiologically based kinetic (NG-PBK) models in support of regulatory decision making

The fields of toxicology and chemical risk assessment seek to reduce, and eventually replace, the use of animals for the prediction of toxicity in humans. In this context, physiologically based kinetic (PBK) modelling based on in vitro and in silico kinetic data has the potential to a play significant role in reducing animal testing, by providing a methodology capable of incorporating in vitro human data to facilitate the development of in vitro to in vivo extrapolation of hazard information. In the present article, we discuss the challenges in: 1) applying PBK modelling to support regulatory decision making under the toxicology and risk-assessment paradigm shift towards animal replacement; 2) constructing PBK models without in vivo animal kinetic data, while relying solely on in vitro or in silico methods for model parameterization; and 3) assessing the validity and credibility of PBK models built largely using non-animal data. The strengths, uncertainties, and limitations of PBK models developed using in vitro or in silico data are discussed in an effort to establish a higher degree of confidence in the application of such models in a regulatory context. The article summarises the outcome of an expert workshop hosted by the European Commission Joint Research Centre (EC-JRC) - European Union Reference Laboratory for Alternatives to Animal Testing (EURL ECVAM), on "Physiologically-Based Kinetic modelling in risk assessment - reaching a whole new level in regulatory decision-making" held in Ispra, Italy, in November 2016, along with results from an international survey conducted in 2017 and recently reported activities occurring within the PBK modelling field. The discussions presented herein highlight the potential applications of next generation (NG)-PBK modelling, based on new data streams.

A framework for application of quantitative property-property relationships (QPPRs) in physiologically based pharmacokinetic (PBPK) models for high-throughput prediction of internal dose of inhaled organic chemicals

New generation of toxicological tests and assessment strategies require validated toxicokinetic data or models that are lacking for most chemicals. This study aimed at developing a quantitative property-property relationship (QPPR)-based human physiologically based pharmacokinetic (PBPK) modeling framework for high-throughput predictions of inhalation toxicokinetics of organic chemicals. A PBPK model was parameterized with QPPR-derived values for hepatic clearance (CL_h) and partition coefficients (P) [blood:air (P_ba) and tissue:air (P_ta) and tissue:blood (P_tb)]. The model was initially applied to an evaluation dataset of 40 organic chemicals in the applicability domain, and then to an expanded dataset of 249 organic chemicals from diverse chemical classes. 'Batch' analyses were performed for rapid assessments of hundreds of chemicals. The simulations of inhalation toxicokinetics following an 8-h exposure to 1 ppm of each chemical were successful. The mean ratios of their predicted-to-experimental values were within a factor of 1.36-2.36 for P_tb and 1.18 for CL_h, for 80% of the chemicals in the evaluation dataset. The predicted 24-h area under the venous blood concentration-time curve (AUC₂₄) values were within the predicted envelopes obtained while using experimental values of P_ba and considering either no or maximal hepatic extraction. The reliability analysis (based on combined sensitivity and uncertainty analyses) indicated that AUC₂₄ predictions for 55% of the expanded dataset were moderately to highly reliable, with 46% exhibiting highly reliable values. Overall, the modeling framework suggests that molecular structure and chemical properties can together be effectively used to obtain first-cut estimates of the toxicokinetics of data-poor organic chemicals for screening and prioritization purposes.

A Semi-Physiologically Based Pharmacokinetic Model Describing the Altered Metabolism of Midazolam Due to Inflammation in Mice

PURPOSE

To investigate influence of inflammation on metabolism and pharmacokinetics (PK) of midazolam (MDZ) and construct a semi-physiologically based pharmacokinetic (PBPK) model to predict PK in mice with inflammatory disease.

METHODS

Glucose-6-phosphate isomerase (GPI)-mediated inflammation was used as a preclinical model of arthritis in DBA/1 mice. CYP3A substrate MDZ was selected to study changes in metabolism and PK during the inflammation. The semi-PBPK model was constructed using mouse physiological parameters, liver microsome metabolism, and healthy animal PK data. In addition, serum cytokine, and liver-CYP (cytochrome P450 enzymes) mRNA levels were examined.

RESULTS

The in vitro metabolite formation rate was suppressed in liver microsomes prepared from the GPI-treated mice as compared to the healthy mice. Further, clearance of MDZ was reduced during inflammation as compared to the healthy group. Finally, the semi-PBPK model was used to predict PK of MDZ after GPI-mediated inflammation. IL-6 and TNF-α levels were elevated and liver-cyp3a11 mRNA was reduced after GPI treatment.

CONCLUSION

The semi-PBPK model successfully predicted PK parameters of MDZ in the disease state. The model may be applied to predict PK of other drugs under disease conditions using healthy animal PK and liver microsomal data as inputs.

Modeling liver electrical conductivity during hypertonic injection

Metastases in the liver frequently grow as scattered tumor nodules that neither can be removed by surgical resection nor focally ablated. Previously, we have proposed a novel technique based on irreversible electroporation that may be able to simultaneously treat all nodules in the liver while sparing healthy tissue. The proposed technique requires increasing the electrical conductivity of healthy liver by injecting a hypersaline solution through the portal vein. Aiming to assess the capability of increasing the global conductivity of the liver by means of hypersaline fluids, here, it is presented a mathematical model that estimates the NaCl distribution within the liver and the resulting conductivity change. The model fuses well-established compartmental pharmacokinetic models of the organ with saline injection models used for resuscitation treatments, and it considers changes in sinusoidal blood viscosity because of the hypertonicity of the solution. Here, it is also described a pilot experimental study in pigs in which different volumes of NaCl 20% (from 100 to 200 mL) were injected through the portal vein at different flow rates (from 53 to 171 mL/minute). The in vivo conductivity results fit those obtained by the model, both quantitatively and qualitatively, being able to predict the maximum conductivity with a 14.6% average relative error. The maximum conductivity value was 0.44 second/m, which corresponds to increasing 4 times the mean basal conductivity (0.11 second/m). The results suggest that the presented model is well suited for predicting on liver conductivity changes during hypertonic saline injection.

Fitting tissue chips and microphysiological systems into the grand scheme of medicine, biology, pharmacology, and toxicology

Microphysiological systems (MPS), which include engineered organoids (EOs), single organ/tissue chips (TCs), and multiple organs interconnected to create miniature in vitro models of human physiological systems, are rapidly becoming effective tools for drug development and the mechanistic understanding of tissue physiology and pathophysiology. The second MPS thematic issue of Experimental Biology and Medicine comprises 15 articles by scientists and engineers from the National Institutes of Health, the IQ Consortium, the Food and Drug Administration, and Environmental Protection Agency, an MPS company, and academia. Topics include the progress, challenges, and future of organs-on-chips, dissemination of TCs into Pharma, children's health protection, liver zonation, liver chips and their coupling to interconnected systems, gastrointestinal MPS, maturation of immature cardiomyocytes in a heart-on-a-chip, coculture of multiple cell types in a human skin construct, use of synthetic hydrogels to create EOs that form neural tissue models, the blood-brain barrier-on-a-chip, MPS models of coupled female reproductive organs, coupling MPS devices to create a body-on-a-chip, and the use of a microformulator to recapitulate endocrine circadian rhythms. While MPS hardware has been relatively stable since the last MPS thematic issue, there have been significant advances in cell sourcing, with increased reliance on human-induced pluripotent stem cells, and in characterization of the genetic and functional cell state in MPS bioreactors. There is growing appreciation of the need to minimize perfusate-to-cell-volume ratios and respect physiological scaling of coupled TCs. Questions asked by drug developers are followed by an analysis of the potential value, costs, and needs of Pharma. Of highest value and lowest switching costs may be the development of MPS disease models to aid in the discovery of disease mechanisms; novel compounds including probes, leads, and clinical candidates; and mechanism of action of drug candidates. Impact statement Microphysiological systems (MPS), which include engineered organoids and both individual and coupled organs-on-chips and tissue chips, are a rapidly growing topic of research that addresses the known limitations of conventional cellular monoculture on flat plastic - a well-perfected set of techniques that produces reliable, statistically significant results that may not adequately represent human biology and disease. As reviewed in this article and the others in this thematic issue, MPS research has made notable progress in the past three years in both cell sourcing and characterization. As the field matures, currently identified challenges are being addressed, and new ones are being recognized. Building upon investments by the Defense Advanced Research Projects Agency, National Institutes of Health, Food and Drug Administration, Defense Threat Reduction Agency, and Environmental Protection Agency of more than $200 million since 2012 and sizable corporate spending, academic and commercial players in the MPS community are demonstrating their ability to meet the translational challenges required to apply MPS technologies to accelerate drug development and advance toxicology.

Investigating the state of physiologically based kinetic modelling practices and challenges associated with gaining regulatory acceptance of model applications

Physiologically based kinetic (PBK) models are used widely throughout a number of working sectors, including academia and industry, to provide insight into the dosimetry related to observed adverse health effects in humans and other species. Use of these models has increased over the last several decades, especially in conjunction with emerging alternative methods to animal testing, such as in vitro studies and data-driven in silico quantitative-structure-activity-relationship (QSAR) predictions. Experimental information derived from these new approach methods can be used as input for model parameters and allows for increased confidence in models for chemicals that did not have in vivo data for model calibration. Despite significant advancements in good modelling practice (GMP) for model development and evaluation, there remains some reluctance among regulatory agencies to use such models during the risk assessment process. Here, the results of a survey disseminated to the modelling community are presented in order to inform the frequency of use and applications of PBK models in science and regulatory submission. Additionally, the survey was designed to identify a network of investigators involved in PBK modelling and knowledgeable of GMP so that they might be contacted in the future for peer review of PBK models, especially in regards to vetting the models to such a degree as to gain a greater acceptance for regulatory purposes.

•

Physiologically Based kinetic (PBK) models are used widely in academia, industry, and government.

•

Good modelling practice (GMP) for model development and evaluation continues to expand.

•

Further guidance for establishing GMP is called for.

•

There remains some reluctance among regulatory agencies to use PBK models.

•

The next generation of PBK models could be developed using only data from in vitro and in silico methods.

Automated workflows for modelling chemical fate, kinetics and toxicity

Automation is universal in today's society, from operating equipment such as machinery, in factory processes, to self-parking automobile systems. While these examples show the efficiency and effectiveness of automated mechanical processes, automated procedures that support the chemical risk assessment process are still in their infancy. Future human safety assessments will rely increasingly on the use of automated models, such as physiologically based kinetic (PBK) and dynamic models and the virtual cell based assay (VCBA). These biologically-based models will be coupled with chemistry-based prediction models that also automate the generation of key input parameters such as physicochemical properties. The development of automated software tools is an important step in harmonising and expediting the chemical safety assessment process. In this study, we illustrate how the KNIME Analytics Platform can be used to provide a user-friendly graphical interface for these biokinetic models, such as PBK models and VCBA, which simulates the fate of chemicals in vivo within the body and in vitro test systems respectively.

•

The VCBA is a mathematical model that simulates in vitro fate of chemicals and the corresponding cellular effect.

•

The VCBA has been implemented in an open access web-based KNIME platform for ease of use.

•

KNIME Analytics Platform can be used to provide a user-friendly graphical interface for biokinetic models.