Evaluation and Optimization of Pharmacokinetic Models for in Vitro to in Vivo Extrapolation of Estrogenic Activity for Environmental Chemicals

Informing Hazard Identification and Risk Characterization of Environmental Chemicals by Combining Transcriptomic and Functional Data from Human-Induced Pluripotent Stem-Cell-Derived Cardiomyocytes

Environmental chemicals may contribute to the global burden of cardiovascular disease, but experimental data are lacking to determine which substances pose the greatest risk. Human-induced pluripotent stem cell (iPSC)-derived cardiomyocytes are a high-throughput cardiotoxicity model that is widely used to test drugs and chemicals; however, most studies focus on exploring electro-physiological readouts. Gene expression data may provide additional molecular insights to be used for both mechanistic interpretation and dose-response analyses. Therefore, we hypothesized that both transcriptomic and functional data in human iPSC-derived cardiomyocytes may be used as a comprehensive screening tool to identify potential cardiotoxicity hazards and risks of the chemicals. To test this hypothesis, we performed concentration-response analysis of 464 chemicals from 12 classes, including both pharmaceuticals and nonpharmaceutical substances. Functional effects (beat frequency, QT prolongation, and asystole), cytotoxicity, and whole transcriptome response were evaluated. Points of departure were derived from phenotypic and transcriptomic data, and risk characterization was performed. Overall, 244 (53%) substances were active in at least one phenotype; as expected, pharmaceuticals with known cardiac liabilities were the most active. Positive chronotropy was the functional phenotype activated by the largest number of tested chemicals. No chemical class was particularly prone to pose a potential hazard to cardiomyocytes; a varying proportion (10-44%) of substances in each class had effects on cardiomyocytes. Transcriptomic data showed that 69 (15%) substances elicited significant gene expression changes; most perturbed pathways were highly relevant to known key characteristics of human cardiotoxicants. The bioactivity-to-exposure ratios showed that phenotypic- and transcriptomic-based POD led to similar results for risk characterization. Overall, our findings demonstrate how the integrative use of in vitro transcriptomic and phenotypic data from iPSC-derived cardiomyocytes not only offers a complementary approach for hazard and risk prioritization, but also enables mechanistic interpretation of the in vitro test results to increase confidence in decision-making.

Challenging the status quo: a framework for mechanistic and human-relevant cardiovascular safety screening

Traditional approaches to preclinical drug safety assessment have generally protected human patients from unintended adverse effects. However, these assessments typically occur too late to make changes in the formulation or in phase 1 and beyond, are highly dependent on animal studies and have the potential to lead to the termination of useful drugs due to liabilities in animals that are not applicable in patients. Collectively, these elements come at great detriment to both patients and the drug development sector. This phenomenon is particularly problematic in the area of cardiovascular safety assessment where preclinical attrition is high. We believe that a more efficient and translational approach can be defined. A multi-tiered assessment that leverages our understanding of human cardiovascular biology, applies human cell-based in vitro characterizations of cardiovascular responses to insult, and incorporates computational models of pharmacokinetic relationships would enable earlier and more translational identification of human-relevant liabilities. While this will take time to develop, the ultimate goal would be to implement such assays both in the lead selection phase as well as through regulatory phases.

Interpretation of in vitro concentration-response data for risk assessment and regulatory decision-making: Report from the 2022 IWGT quantitative analysis expert working group meeting

Quantitative risk assessments of chemicals are routinely performed using in vivo data from rodents; however, there is growing recognition that non-animal approaches can be human-relevant alternatives. There is an urgent need to build confidence in non-animal alternatives given the international support to reduce the use of animals in toxicity testing where possible. In order for scientists and risk assessors to prepare for this paradigm shift in toxicity assessment, standardization and consensus on in vitro testing strategies and data interpretation will need to be established. To address this issue, an Expert Working Group (EWG) of the 8th International Workshop on Genotoxicity Testing (IWGT) evaluated the utility of quantitative in vitro genotoxicity concentration-response data for risk assessment. The EWG first evaluated available in vitro methodologies and then examined the variability and maximal response of in vitro tests to estimate biologically relevant values for the critical effect sizes considered adverse or unacceptable. Next, the EWG reviewed the approaches and computational models employed to provide human-relevant dose context to in vitro data. Lastly, the EWG evaluated risk assessment applications for which in vitro data are ready for use and applications where further work is required. The EWG concluded that in vitro genotoxicity concentration-response data can be interpreted in a risk assessment context. However, prior to routine use in regulatory settings, further research will be required to address the remaining uncertainties and limitations.

Key challenges for in vitro testing of tobacco products for regulatory applications: Recommendations for dosimetry

The Institute for In Vitro Sciences (IIVS) is sponsoring a series of workshops to develop recommendations for optimal scientific and technical approaches for conducting in vitro assays to assess potential toxicity within and across tobacco and various next-generation products (NGPs) including heated tobacco products (HTPs) and electronic nicotine delivery systems (ENDSs). This publication was developed by a working group of the workshop members in conjunction with the sixth workshop in that series entitled "Dosimetry for conducting in vitro evaluations" and focuses on aerosol dosimetry for aerosol exposure to combustible cigarettes, HTP, and ENDS aerosolized tobacco products and summarizes the key challenges as well as documenting areas for future research.

The phytoestrogenic potential of flavonoid glycosides from Selaginella moellendorffii via ERα-dependent signaling pathway

ETHNOPHARMACOLOGICAL RELEVANCE

Selaginella moellendorffii (SM) has been applied as an ethnic drug to treat conditions such as osteoporosis, idiopathic thrombocytopenic purpura, and chronic inflammation. It is known to be rich in flavonoids, including apigenin glycosides and unique elements of bioflavonoids.

AIM OF THE STUDY

To investigate estrogen-like constituents of SM and the possible mechanism.

MATERIALS AND METHODS

We identified the main components in liquid chromatography and liquid chromatography-mass spectrometry. The estrogenic effects were examined using a recombinant yeast screening assay, an E-screen cell proliferation assay, and an in vivo uterotrophic assay.

RESULTS

Flavonoid glycosides extract, some flavonoid glycosides, and apigenin showed estrogen agonistic activity in the yeast screening assay. They also induced cell proliferation in estrogen receptor-positive (ER+) cells but not in estrogen receptor-negative (ER-) cells. Consistently, the protein expression of ERα, phosphorylation protein kinase B (p-AKT), phosphatidylinositol 3 kinase (PI3K), phosphorylation mammalian target of rapamycin (p-mTOR), phosphorylation 38,000-Da protein (p-P38), and phosphorylation extracellular-regulated kinase 1/2 (p-ERK1/2) elevated following treatment with flavonoid glycoside extract (P < 0.01 or P < 0.05). These effects could be blocked by ER antagonist or ERα antagonist but not be blocked by ERβ antagonist. In vivo assay, flavonoid glycoside extract could significantly increase body weight, serum estradiol level, uterine wet weight, alter uterine morphology, and promote ERα protein expression (P < 0.01 or P < 0.05).

CONCLUSIONS

ERα induction via mitogen-activated protein kinases (MAPK) and PI3K/Akt/mTOR pathways might be the possible mechanism underlying the phytoestrogen effect of SM, and the flavonoid glycosides might be the critical estrogenic constituents.

Human health risk assessment of 6:2 Cl-PFESA through quantitative in vitro to in vivo extrapolation by integrating cell-based assays, an epigenetic key event, and physiologically based pharmacokinetic modeling

Human health risk assessment of chemicals is essential but often relies on time-consuming and animal and labor-extensive procedures. Here, we develop a population-based, quantitative in vitro to in vivo extrapolation (QIVIVE) approach which depended on cellular effects monitored by in vitro assays, considered chemical internal concentration determined by LC-MS/MS, extrapolated into in vivo target tissue concentration through physiologically based pharmacokinetic (PBPK) modelling, and assessed populational health risk using in silico modelling. By applying this QIVIVE approach to 6:2 chlorinated polyfluorinated ether sulfonate (6:2 Cl-PFESA), as a representative of the emerging pollutants, we find that 6:2 Cl-PFESA disturbed lipid homeostasis in HepG2 cells through enhancement of lipid accumulation and fatty acid β-oxidation, during which miR-93-5p served as a key event towards toxicity and thus, could serve as an efficient toxicity marker for risk assessment; further, the disruption potency of lipid homeostasis of 6:2 Cl-PFESA for the most of studied populations in China might be of moderate concern. Together, our approach improved the reliability of QIVIVE during human health risk assessment, which can readily be used for other chemicals.

A Pragmatic Framework for the Application of New Approach Methodologies in One Health Toxicological Risk Assessment

Globally, industries and regulatory authorities are faced with an urgent need to assess the potential adverse effects of chemicals more efficiently by embracing new approach methodologies (NAMs). NAMs include cell and tissue methods (in vitro), structure-based/toxicokinetic models (in silico), methods that assess toxicant interactions with biological macromolecules (in chemico), and alternative models. Increasing knowledge on chemical toxicokinetics (what the body does with chemicals) and toxicodynamics (what the chemicals do with the body) obtained from in silico and in vitro systems continues to provide opportunities for modernizing chemical risk assessments. However, directly leveraging in vitro and in silico data for derivation of human health-based reference values has not received regulatory acceptance due to uncertainties in extrapolating NAM results to human populations, including metabolism, complex biological pathways, multiple exposures, interindividual susceptibility and vulnerable populations. The objective of this article is to provide a standardized pragmatic framework that applies integrated approaches with a focus on quantitative in vitro to in vivo extrapolation (QIVIVE) to extrapolate in vitro cellular exposures to human equivalent doses from which human reference values can be derived. The proposed framework intends to systematically account for the complexities in extrapolation and data interpretation to support sound human health safety decisions in diverse industrial sectors (food systems, cosmetics, industrial chemicals, pharmaceuticals etc.). Case studies of chemical entities, using new and existing data, are presented to demonstrate the utility of the proposed framework while highlighting potential sources of human population bias and uncertainty, and the importance of Good Method and Reporting Practices.

Quantitative in vitro to in vivo extrapolation for developmental toxicity potency of valproic acid analogues

BACKGROUND

The developmental toxicity potential (dTP) concentration from the devTOX quickPredict (devTOX^qP ) assay, a metabolomics-based human induced pluripotent stem cell assay, predicts a chemical's developmental toxicity potency. Here, in vitro to in vivo extrapolation (IVIVE) approaches were applied to address whether the devTOX^qP assay could quantitatively predict in vivo developmental toxicity lowest effect levels (LELs) for the prototypical teratogen valproic acid (VPA) and a group of structural analogues.

METHODS

VPA and a series of structural analogues were tested with the devTOX^qP assay to determine dTP concentration and we estimated the equivalent administered doses (EADs) that would lead to plasma concentrations equivalent to the in vitro dTP concentrations. The EADs were compared to the LELs in rat developmental toxicity studies, human clinical doses, and EADs reported using other in vitro assays. To evaluate the impact of different pharmacokinetic (PK) models on IVIVE outcomes, we compared EADs predicted using various open-source and commercially available PK and physiologically based PK (PBPK) models. To evaluate the effect of in vitro kinetics, an equilibrium distribution model was applied to translate dTP concentrations to free medium concentrations before subsequent IVIVE analyses.

RESULTS

The EAD estimates for the VPA analogues based on different PK/PBPK models were quantitatively similar to in vivo data from both rats and humans, where available, and the derived rank order of the chemicals was consistent with observed in vivo developmental toxicity. Different models were identified that provided accurate predictions for rat prenatal LELs and conservative estimates of human safe exposure. The impact of in vitro kinetics on EAD estimates is chemical-dependent. EADs from this study were within range of predicted doses from other in vitro and model organism data.

CONCLUSIONS

This study highlights the importance of pharmacokinetic considerations when using in vitro assays and demonstrates the utility of the devTOX^qP human stem cell-based platform to quantitatively assess a chemical's developmental toxicity potency.

Dynamic Mass Balance Modeling for Chemical Distribution Over Time in In Vitro Systems With Repeated Dosing

As toxicologists and risk assessors move away from animal testing and more toward using in vitro models and biological modeling, it is necessary to produce tools to quantify the chemical distribution within the in vitro environment prior to extrapolating in vitro concentrations to human equivalent doses. Although models predicting chemical distribution in vitro have been developed, very little has been done for repeated dosing scenarios, which are common in prolonged experiments where the medium needs to be refreshed. Failure to account for repeated dosing may lead to inaccurate estimations of exposure and introduce bias into subsequent in vitro to in vivo extrapolations. Our objectives were to develop a dynamic mass balance model for repeated dosing in in vitro systems; to evaluate model accuracy against experimental data; and to perform illustrative simulations to assess the impact of repeated doses on predicted cellular concentrations. A novel dynamic in vitro partitioning mass balance model (IV-MBM DP v1.0) was created based on the well-established fugacity approach. We parameterized and applied the dynamic mass balance model to single dose and repeat dosing scenarios, and evaluated the predicted medium and cellular concentrations against available empirical data. We also simulated repeated dosing scenarios for organic chemicals with a range of partitioning properties and compared the in vitro distributions over time. In single dose scenarios, for which only medium concentrations were available, simulated concentrations predicted measured concentrations with coefficients of determination (R2) of 0.85–0.89, mean absolute error within a factor of two and model bias of nearly one. Repeat dose scenario simulations displayed model bias &lt;2 within the cell lysate, and ∼1.5-3 in the medium. The concordance between simulated and available experimental data supports the predictive capacity of the IV-MBM DP v1.0 tool, but further evaluation as empirical data becomes available is warranted, especially for cellular concentrations.

Impact of High-Throughput Model Parameterization and Data Uncertainty on Thyroid-Based Toxicological Estimates for Pesticide Chemicals

Chemical-induced alteration of maternal thyroid hormone levels may increase the risk of adverse neurodevelopmental outcomes in offspring. US federal risk assessments rely almost exclusively on apical endpoints in animal models for deriving points of departure (PODs). New approach methodologies (NAMs) such as high-throughput screening (HTS) and mechanistically informative in vitro human cell-based systems, combined with in vitro to in vivo extrapolation (IVIVE), supplement in vivo studies and provide an alternative approach to calculate/determine PODs. We examine how parameterization of IVIVE models impacts the comparison between IVIVE-derived equivalent administered doses (EADs) from thyroid-relevant in vitro assays and the POD values that serve as the basis for risk assessments. Pesticide chemicals with thyroid-based in vitro bioactivity data from the US Tox21 HTS program were included (n = 45). Depending on the model structure used for IVIVE analysis, up to 35 chemicals produced EAD values lower than the POD. A total of 10 chemicals produced EAD values higher than the POD regardless of the model structure. The relationship between IVIVE-derived EAD values and the in vivo-derived POD values is highly dependent on model parameterization. Here, we derive a range of potentially thyroid-relevant doses that incorporate uncertainty in modeling choices and in vitro assay data.

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.

Considerations for Improving Metabolism Predictions for In Vitro to In Vivo Extrapolation

High-throughput (HT) in vitro to in vivo extrapolation (IVIVE) is an integral component in new approach method (NAM)-based risk assessment paradigms, for rapidly translating in vitro toxicity assay results into the context of in vivo exposure. When coupled with rapid exposure predictions, HT-IVIVE supports the use of HT in vitro assays for risk-based chemical prioritization. However, the reliability of prioritization based on HT bioactivity data and HT-IVIVE can be limited as the domain of applicability of current HT-IVIVE is generally restricted to intrinsic clearance measured primarily in pharmaceutical compounds. Further, current approaches only consider parent chemical toxicity. These limitations occur because current state-of-the-art HT prediction tools for clearance and metabolite kinetics do not provide reliable data to support HT-IVIVE. This paper discusses current challenges in implementation of IVIVE for prioritization and risk assessment and recommends a path forward for addressing the most pressing needs and expanding the utility of IVIVE.

A framework for chemical safety assessment incorporating new approach methodologies within REACH

The long-term investment in new approach methodologies (NAMs) within the EU and other parts of the world is beginning to result in an emerging consensus of how to use information from in silico, in vitro and targeted in vivo sources to assess the safety of chemicals. However, this methodology is being adopted very slowly for regulatory purposes. Here, we have developed a framework incorporating in silico, in vitro and in vivo methods designed to meet the requirements of REACH in which both hazard and exposure can be assessed using a tiered approach. The outputs from each tier are classification categories, safe doses, and risk assessments, and progress through the tiers depends on the output from previous tiers. We have exemplified the use of the framework with three examples. The outputs were the same or more conservative than parallel assessments based on conventional studies. The framework allows a transparent and phased introduction of NAMs in chemical safety assessment and enables science-based safety decisions which provide the same level of public health protection using fewer animals, taking less time, and using less financial and expert resource. Furthermore, it would also allow new methods to be incorporated as they develop through continuous selective evolution rather than periodic revolution.

Update and Evaluation of a High-Throughput In Vitro Mass Balance Distribution Model: IV-MBM EQP v2.0

This study demonstrates the utility of an updated mass balance model for predicting the distribution of organic chemicals in in vitro test systems (IV-MBM EQP v2.0) and evaluates its performance with empirical data. The IV-MBM EQP v2.0 tool was parameterized and applied to four independent data sets with measured ratios of bulk medium or freely-dissolved to initial nominal concentrations (e.g., C24/C0 where C24 is the measured concentration after 24 h of exposure and C0 is the initial nominal concentration). Model performance varied depending on the data set, chemical properties (e.g., "volatiles" vs. "non-volatiles", neutral vs. ionizable organics), and model assumptions but overall is deemed acceptable. For example, the r² was greater than 0.8 and the mean absolute error (MAE) in the predictions was less than a factor of two for most neutral organics included. Model performance was not as good for the ionizable organic chemicals included but the r² was still greater than 0.7 and the MAE less than a factor of three. The IV-MBM EQP v2.0 model was subsequently applied to several hundred chemicals on Canada's Domestic Substances List (DSL) with nominal effects data (AC50s) reported for two in vitro assays. We report the frequency of chemicals with AC50s corresponding to predicted cell membrane concentrations in the baseline toxicity range (i.e., >20-60 mM) and tabulate the number of chemicals with "volatility issues" (majority of chemical in headspace) and "solubility issues" (freely-dissolved concentration greater than water solubility after distribution). In addition, the predicted "equivalent EQP blood concentrations" (i.e., blood concentration at equilibrium with predicted cellular concentration) were compared to the AC50s as a function of hydrophobicity (log octanol-water partition or distribution ratio). The predicted equivalent EQP blood concentrations exceed the AC50 by up to a factor of 100 depending on hydrophobicity and assay conditions. The implications of using AC50s as direct surrogates for human blood concentrations when estimating the oral equivalent doses using a toxicokinetic model (i.e., reverse dosimetry) are then briefly discussed.

An integrated exposure and pharmacokinetic modeling framework for assessing population-scale risks of phthalates and their substitutes

To effectively incorporate in vitro-in silico-based methods into the regulation of consumer product safety, a quantitative connection between product phthalate concentrations and in vitro bioactivity data must be established for the general population. We developed, evaluated, and demonstrated a modeling framework that integrates exposure and pharmacokinetic models to convert product phthalate concentrations into population-scale risks for phthalates and their substitutes. A probabilistic exposure model was developed to generate the distribution of multi-route exposures based on product phthalate concentrations, chemical properties, and human activities. Pharmacokinetic models were developed to simulate population toxicokinetics using Bayesian analysis via the Markov chain Monte Carlo method. Both exposure and pharmacokinetic models demonstrated good predictive capability when compared with worldwide studies. The distributions of exposures and pharmacokinetics were integrated to predict the population distributions of internal dosimetry. The predicted distributions showed reasonable agreement with the U.S. biomonitoring surveys of urinary metabolites. The "source-to-outcome" local sensitivity analysis revealed that food contact materials had the greatest impact on body burden for di(2-ethylhexyl) adipate (DEHA), di-2-ethylhexyl phthalate (DEHP), di(isononyl) cyclohexane-1,2-dicarboxylate (DINCH), and di(2-propylheptyl) phthalate (DPHP), whereas the body burden of diethyl phthalate (DEP) was most sensitive to the concentration in personal care products. The upper bounds of predicted plasma concentrations showed no overlap with ToxCast in vitro bioactivity values. Compared with the in vitro-to-in vivo extrapolation (IVIVE) approach, the integrated modeling framework has significant advantages in mapping product phthalate concentrations to multi-route risks, and thus is of great significance for regulatory use with a relatively low input requirement. Further integration with new approach methodologies will facilitate these in vitro-in silico-based risk assessments for a broad range of products containing an equally broad range of chemicals.

Cardiotoxicity Hazard and Risk Characterization of ToxCast Chemicals Using Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes from Multiple Donors

Heart disease remains a significant human health burden worldwide with a significant fraction of morbidity attributable to environmental exposures. However, the extent to which the thousands of chemicals in commerce and the environment may contribute to heart disease morbidity is largely unknown, because in contrast to pharmaceuticals, environmental chemicals are seldom tested for potential cardiotoxicity. Human induced pluripotent stem cell (iPSC)-derived cardiomyocytes have become an informative in vitro model for cardiotoxicity testing of drugs with the availability of cells from multiple individuals allowing in vitro testing of population variability. In this study, we hypothesized that a panel of iPSC-derived cardiomyocytes from healthy human donors can be used to screen for the potential cardiotoxicity hazard and risk of environmental chemicals. We conducted concentration-response testing of 1029 chemicals (drugs, pesticides, flame retardants, polycyclic aromatic hydrocarbons (PAHs), plasticizers, industrial chemicals, food/flavor/fragrance agents, etc.) in iPSC-derived cardiomyocytes from 5 donors. We used kinetic calcium flux and high-content imaging to derive quantitative measures as inputs into Bayesian population concentration-response modeling of the effects of each chemical. We found that many environmental chemicals pose a hazard to human cardiomyocytes in vitro with more than half of all chemicals eliciting positive or negative chronotropic or arrhythmogenic effects. However, most of the tested environmental chemicals for which human exposure and high-throughput toxicokinetics data were available had wide margins of exposure and, thus, do not appear to pose a significant human health risk in a general population. Still, relatively narrow margins of exposure (<100) were estimated for some perfuoroalkyl substances and phthalates, raising concerns that cumulative exposures may pose a cardiotoxicity risk. Collectively, this study demonstrated the value of using a population-based human in vitro model for rapid, high-throughput hazard and risk characterization of chemicals for which little to no cardiotoxicity data are available from guideline studies in animals.

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.

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.

Predictive modeling of biological responses in the rat liver using in vitro Tox21 bioactivity: Benefits from high-throughput toxicokinetics

Computational methods are needed to more efficiently leverage data from in vitro cell-based models to predict what occurs within whole body systems after chemical insults. This study set out to test the hypothesis that in vitro high-throughput screening (HTS) data can more effectively predict in vivo biological responses when chemical disposition and toxicokinetic (TK) modeling are employed. In vitro HTS data from the Tox21 consortium were analyzed in concert with chemical disposition modeling to derive nominal, aqueous, and intracellular estimates of concentrations eliciting 50% maximal activity. In vivo biological responses were captured using rat liver transcriptomic data from the DrugMatrix and TG-Gates databases and evaluated for pathway enrichment. In vivo dosing data were translated to equivalent body concentrations using HTTK modeling. Random forest models were then trained and tested to predict in vivo pathway-level activity across 221 chemicals using in vitro bioactivity data and physicochemical properties as predictor variables, incorporating methods to address imbalanced training data resulting from high instances of inactivity. Model performance was quantified using the area under the receiver operator characteristic curve (AUC-ROC) and compared across pathways for different combinations of predictor variables. All models that included toxicokinetics were found to outperform those that excluded toxicokinetics. Biological interpretation of the model features revealed that rather than a direct mapping of in vitro assays to in vivo pathways, unexpected combinations of multiple in vitro assays predicted in vivo pathway-level activities. To demonstrate the utility of these findings, the highest-performing model was leveraged to make new predictions of in vivo biological responses across all biological pathways for remaining chemicals tested in Tox21 with adequate data coverage (n = 6617). These results demonstrate that, when chemical disposition and toxicokinetics are carefully considered, in vitro HT screening data can be used to effectively predict in vivo biological responses to chemicals.

An integrated chemical environment with tools for chemical safety testing

Moving toward species-relevant chemical safety assessments and away from animal testing requires access to reliable data to develop and build confidence in new approaches. The Integrated Chemical Environment (ICE) provides tools and curated data centered around chemical safety assessment. This article describes updates to ICE, including improved accessibility and interpretability of in vitro data via mechanistic target mapping and enhanced interactive tools for in vitro to in vivo extrapolation (IVIVE). Mapping of in vitro assay targets to toxicity endpoints of regulatory importance uses literature-based mode-of-action information and controlled terminology from existing knowledge organization systems to support data interoperability with external resources. The most recent ICE update includes Tox21 high-throughput screening data curated using analytical chemistry data and assay-specific parameters to eliminate potential artifacts or unreliable activity. Also included are physicochemical/ADME parameters for over 800,000 chemicals predicted by quantitative structure-activity relationship models. These parameters are used by the new ICE IVIVE tool in combination with the U.S. Environmental Protection Agency's httk R package to estimate in vivo exposures corresponding to in vitro bioactivity concentrations from stored or user-defined assay data. These new ICE features allow users to explore the applications of an expanded data space and facilitate building confidence in non-animal approaches.

Transcriptional changes caused by estrogenic endocrine disrupting chemicals in gonad-mesonephros complexes of genetic male Xenopus laevis: Multiple biomarkers for early detection of testis differentiation disruption

Our recent study revealed some early molecular and cellular events in which 17β-estradiol (E2) disrupted testis differentiation and resulted in feminization in Xenopus laevis (the African clawed frog), an ideal species for studying reproductive endocrine disruption by estrogenic endocrine disrupting chemicals (EDCs). On this basis, we aimed to develop multiple biomarkers for early detection of testis differentiation disruption by estrogenic EDCs in X. laevis. Tadpoles at stage 45/46 were exposed to four known estrogenic EDCs with different estrogenic activities, including E2, diethylstilbestrol (DES), mestranol (MES) and 4-n-nonyphenol (NP). At stage 53, gonadal morphological and histological changes as well as altered sex-dimorphic gene expression in gonad-mesonephros complexes (GMCs) showed that these estrogenic EDCs disrupted testis differentiation and caused feminization to different degrees. Then we measured transcriptional changes of 48 candidate genes, which are believed to be associated with E2-induced testis differentiation alterations, in GMCs at stage 50. As a result, 19 genes were found to be transcriptionally altered by all test chemicals and proposed as promising biomarkers for early detection of testis differentiation disruption by estrogenic EDCs. Finally, all biomarker responses were integrated as integrated biomarker response (IBR) index to characterize testis differentiation disruption by these estrogenic EDCs in X. laevis. Compared with the methods used in previous studies, the multiple biomarker test using X. laevis at early developmental stages largely shortens the exposure duration, thereby achieving the goal of rapid detection. Certainly, the biomarker test needs further validations in the future study.

The EU-ToxRisk method documentation, data processing and chemical testing pipeline for the regulatory use of new approach methods

Hazard assessment, based on new approach methods (NAM), requires the use of batteries of assays, where individual tests may be contributed by different laboratories. A unified strategy for such collaborative testing is presented. It details all procedures required to allow test information to be usable for integrated hazard assessment, strategic project decisions and/or for regulatory purposes. The EU-ToxRisk project developed a strategy to provide regulatorily valid data, and exemplified this using a panel of > 20 assays (with > 50 individual endpoints), each exposed to 19 well-known test compounds (e.g. rotenone, colchicine, mercury, paracetamol, rifampicine, paraquat, taxol). Examples of strategy implementation are provided for all aspects required to ensure data validity: (i) documentation of test methods in a publicly accessible database; (ii) deposition of standard operating procedures (SOP) at the European Union DB-ALM repository; (iii) test readiness scoring accoding to defined criteria; (iv) disclosure of the pipeline for data processing; (v) link of uncertainty measures and metadata to the data; (vi) definition of test chemicals, their handling and their behavior in test media; (vii) specification of the test purpose and overall evaluation plans. Moreover, data generation was exemplified by providing results from 25 reporter assays. A complete evaluation of the entire test battery will be described elsewhere. A major learning from the retrospective analysis of this large testing project was the need for thorough definitions of the above strategy aspects, ideally in form of a study pre-registration, to allow adequate interpretation of the data and to ensure overall scientific/toxicological validity.

Comparative analysis of Rapid Equilibrium Dialysis (RED) and solid phase micro-extraction (SPME) methods for In Vitro-In Vivo extrapolation of environmental chemicals

In vitro to in vivo extrapolation (IVIVE) is a critical component of the efforts to prioritize and assess environmental chemicals using high-throughput in vitro assays. The plasma unbound fraction (F_ub) is a key toxicokinetic parameter in IVIVE, and is usually measured via the Rapid Equilibrium Dialysis (RED) assay widely used for pharmaceuticals. However, pharmaceuticals have a narrower range of physicochemical properties than environmental chemicals. Motivated by the observation that high LogK_OW compounds appeared to have disproportionately low F_ub measurements using RED, we added a protein-free control in order to verify equilibration to 100% unbound in the absence of proteins. We found that many high LogK_OW non-pharmaceuticals fail to equilibrate in RED in protein-free controls, and thus had apparent values of F_ub = 0 in plasma. In these cases, Solid Phase Microextraction (SPME) as an alternative method provided an accurate, though more time-consuming, alternative to accurately determine F_ub. We propose an updated IVIVE workflow that adds a protein-free control to the RED protocol, with the use of alternative approaches, such as SPME, in cases where compounds fail to adequately equilibrate. These refinements will provide additional confidence in the use of IVIVE as part of high-throughput screening programs of chemicals.

Using the concordance of in vitro and in vivo data to evaluate extrapolation assumptions

Linking in vitro bioactivity and in vivo toxicity on a dose basis enables the use of high-throughput in vitro assays as an alternative to traditional animal studies. In this study, we evaluated assumptions in the use of a high-throughput, physiologically based toxicokinetic (PBTK) model to relate in vitro bioactivity and rat in vivo toxicity data. The fraction unbound in plasma (f_up) and intrinsic hepatic clearance (Cl_int) were measured for rats (for 67 and 77 chemicals, respectively), combined with f_up and Cl_int literature data for 97 chemicals, and incorporated in the PBTK model. Of these chemicals, 84 had corresponding in vitro ToxCast bioactivity data and in vivo toxicity data. For each possible comparison of in vitro and in vivo endpoint, the concordance between the in vivo and in vitro data was evaluated by a regression analysis. For a base set of assumptions, the PBTK results were more frequently better associated than either the results from a “random” model parameterization or direct comparison of the “untransformed” values of AC₅₀ and dose (performed best in 51%, 28%, and 21% of cases, respectively). We also investigated several assumptions in the application of PBTK for IVIVE, including clearance and internal dose selection. One of the better assumptions sets–restrictive clearance and comparing free in vivo venous plasma concentration with free in vitro concentration–outperformed the random and untransformed results in 71% of the in vitro-in vivo endpoint comparisons. These results demonstrate that applying PBTK improves our ability to observe the association between in vitro bioactivity and in vivo toxicity data in general. This suggests that potency values from in vitro screening should be transformed using in vitro-in vivo extrapolation (IVIVE) to build potentially better machine learning and other statistical models for predicting in vivo toxicity in humans.

Evaluation of androgen assay results using a curated Hershberger database

A set of 39 reference chemicals with reproducible androgen pathway effects in vivo, identified in the companion manuscript [1], were used to interrogate the performance of the ToxCast/Tox 21 androgen receptor (AR) model based on 11 high throughput assays. Cytotoxicity data and specificity confirmation assays were used to distinguish assay loss-of-function from true antagonistic signaling suppression. Overall agreement was 66% (19/29), with ten additional inconclusive chemicals. Most discrepancies were explained using in vitro to in vivo extrapolation to estimate equivalent administered doses. The AR model had 100% positive predictive value for the in vivo response, i.e. there were no false positives, and chemicals with conclusive AR model results (agonist or antagonist) were consistently positive in vivo. Considering the lack of reproducibility of the in vivo Hershberger assay, the in vitro AR model may better predict specific AR interaction and can rapidly and cost-effectively screen thousands of chemicals without using animals.

In vitro to in vivo extrapolation for high throughput prioritization and decision making

In vitro chemical safety testing methods offer the potential for efficient and economical tools to provide relevant assessments of human health risk. To realize this potential, methods are needed to relate in vitro effects to in vivo responses, i.e., in vitro to in vivo extrapolation (IVIVE). Currently available IVIVE approaches need to be refined before they can be utilized for regulatory decision-making. To explore the capabilities and limitations of IVIVE within this context, the U.S. Environmental Protection Agency Office of Research and Development and the National Toxicology Program Interagency Center for the Evaluation of Alternative Toxicological Methods co-organized a workshop and webinar series. Here, we integrate content from the webinars and workshop to discuss activities and resources that would promote inclusion of IVIVE in regulatory decision-making. We discuss properties of models that successfully generate predictions of in vivo doses from effective in vitro concentration, including the experimental systems that provide input parameters for these models, areas of success, and areas for improvement to reduce model uncertainty. Finally, we provide case studies on the uses of IVIVE in safety assessments, which highlight the respective differences, information requirements, and outcomes across various approaches when applied for decision-making.