Quantitativein vitrotoin vivoextrapolation of cell-based toxicity assay results

Developing a predictive model for blood-brain-barrier permeability to explore relevance of in vitro neurotoxicity data for in vivo risk assessment

Introduction

Despite recent rapid advancements in in vitro toxicology, its application to whole-body health outcomes remains limited. Incorporating factors like internal exposure, such as permeability across biomembranes, could improve its relevance. Notably, there is a lack of data and predictive models for blood-brain barrier (BBB) permeability, a proxy for the exposure of target organs to neurotoxicity. We developed a predictive model for BBB permeability to investigate whether it can strengthen the correlation between in vitro and in vivo neurotoxicity data.

Methods

We collected permeability data from parallel artificial membrane permeability assays for brain membranes (PAMPA-BBB) for 106 compounds with varied physicochemical properties. This was utilized to develop an empirical model to expand the potential coverage of chemicals. A list of 23 chemicals with available in vivo and in vitro neurotoxicity data from EPA IRIS and ToxCast was curated to analyze the correlation in toxicity rankings with the Spearman correlation coefficient, with and without the consideration of permeability from our predictive model.

Results

The PAMPA-BBB predictive model showed promising results, with an R2 of 0.71 (measured vs predicted permeabilities). Considering permeability did not improve the correlation between in vitro and in vivo neurotoxicity (0.01 vs -0.11).

Discussion

This weak correlation may stem from model uncertainty and the exclusion of other toxicokinetic processes, along with interspecies toxicodynamics differences. Our results indicate more detailed information on how neurotoxic substances behave inside the body is essential to better utilize the in vitro neurotoxicity data for predicting in vivo toxicity and assessing the risk to the central nervous system.

In vitro to in vivo extrapolation from 3D hiPSC-derived cardiac microtissues and physiologically based pharmacokinetic modeling to inform next-generation arrhythmia risk assessment

Proarrhythmic cardiotoxicity remains a substantial barrier to drug development as well as a major global health challenge. In vitro human pluripotent stem cell-based new approach methodologies have been increasingly proposed and employed as alternatives to existing in vitro and in vivo models that do not accurately recapitulate human cardiac electrophysiology or cardiotoxicity risk. In this study, we expanded the capacity of our previously established 3D human cardiac microtissue model to perform quantitative risk assessment by combining it with a physiologically based pharmacokinetic model, allowing a direct comparison of potentially harmful concentrations predicted in vitro to in vivo therapeutic levels. This approach enabled the measurement of concentration responses and margins of exposure for 2 physiologically relevant metrics of proarrhythmic risk (i.e. action potential duration and triangulation assessed by optical mapping) across concentrations spanning 3 orders of magnitude. The combination of both metrics enabled accurate proarrhythmic risk assessment of 4 compounds with a range of known proarrhythmic risk profiles (i.e. quinidine, cisapride, ranolazine, and verapamil) and demonstrated close agreement with their known clinical effects. Action potential triangulation was found to be a more sensitive metric for predicting proarrhythmic risk associated with the primary mechanism of concern for pharmaceutical-induced fatal ventricular arrhythmias, delayed cardiac repolarization due to inhibition of the rapid delayed rectifier potassium channel, or hERG channel. This study advances human-induced pluripotent stem cell-based 3D cardiac tissue models as new approach methodologies that enable in vitro proarrhythmic risk assessment with high precision of quantitative metrics for understanding clinically relevant cardiotoxicity.

Systematic evaluation of high-throughput PBK modelling strategies for the prediction of intravenous and oral pharmacokinetics in humans

Physiologically based kinetic (PBK) modelling offers a mechanistic basis for predicting the pharmaco-/toxicokinetics of compounds and thereby provides critical information for integrating toxicity and exposure data to replace animal testing with in vitro or in silico methods. However, traditional PBK modelling depends on animal and human data, which limits its usefulness for non-animal methods. To address this limitation, high-throughput PBK modelling aims to rely exclusively on in vitro and in silico data for model generation. Here, we evaluate a variety of in silico tools and different strategies to parameterise PBK models with input values from various sources in a high-throughput manner. We gather 2000 + publicly available human in vivo concentration-time profiles of 200 + compounds (IV and oral administration), as well as in silico, in vitro and in vivo determined compound-specific parameters required for the PBK modelling of these compounds. Then, we systematically evaluate all possible PBK model parametrisation strategies in PK-Sim and quantify their prediction accuracy against the collected in vivo concentration-time profiles. Our results show that even simple, generic high-throughput PBK modelling can provide accurate predictions of the pharmacokinetics of most compounds (87% of Cmax and 84% of AUC within tenfold). Nevertheless, we also observe major differences in prediction accuracies between the different parameterisation strategies, as well as between different compounds. Finally, we outline a strategy for high-throughput PBK modelling that relies exclusively on freely available tools. Our findings contribute to a more robust understanding of the reliability of high-throughput PBK modelling, which is essential to establish the confidence necessary for its utilisation in Next-Generation Risk Assessment.

A UK framework for the assessment and integration of different scientific evidence streams in chemical risk assessment

BACKGROUND

Few methods are available for transparently combining different evidence streams for chemical risk assessment to reach an integrated conclusion on the probability of causation. Hence, the UK Committees on Toxicity (COT) and on Carcinogenicity (COC) have reviewed current practice and developed guidance on how to achieve this in a transparent manner, using graphical visualisation.

METHODS/APPROACH

All lines of evidence, including toxicological, epidemiological, new approach methodologies, and mode of action should be considered, taking account of their strengths/weaknesses in their relative weighting towards a conclusion on the probability of causation. A qualitative estimate of the probability of causation is plotted for each line of evidence and a combined estimate provided.

DISCUSSION/CONCLUSIONS

Guidance is provided on integration of multiple lines of evidence for causation, based on current best practice. Qualitative estimates of probability for each line of evidence are plotted graphically. This ensures a deliberative, consensus conclusion on likelihood of causation is reached. It also ensures clear communication of the influence of the different lines of evidence on the overall conclusion on causality. Issues on which advice from the respective Committees is sought varies considerably, hence the guidance is designed to be sufficiently flexible to meet this need.

Discovery Phase Agrochemical Predictive Safety Assessment Using High Content In Vitro Data to Estimate an In Vivo Toxicity Point of Departure

Utilization of in vitro (cellular) techniques, like Cell Painting and transcriptomics, could provide powerful tools for agrochemical candidate sorting and selection in the discovery process. However, using these models generates challenges translating in vitro concentrations to the corresponding in vivo exposures. Physiologically based pharmacokinetic (PBPK) modeling provides a framework for quantitative in vitro to in vivo extrapolation (IVIVE). We tested whether in vivo (rat liver) transcriptomic and apical points of departure (PODs) could be accurately predicted from in vitro (rat hepatocyte or human HepaRG) transcriptomic PODs or HepaRG Cell Painting PODs using PBPK modeling. We compared two PBPK models, the ADMET predictor and the httk R package, and found httk to predict the in vivo PODs more accurately. Our findings suggest that a rat liver apical and transcriptomic POD can be estimated utilizing a combination of in vitro transcriptome-based PODs coupled with PBPK modeling for IVIVE. Thus, high content in vitro data can be translated with modest accuracy to in vivo models of ultimate regulatory importance to help select agrochemical analogs in early stage discovery program.

The Potential of Cold Atmospheric Pressure Plasmas for the Direct Degradation of Organic Pollutants Derived from the Food Production Industry

Specialized chemicals are used for intensifying food production, including boosting meat and crop yields. Among the applied formulations, antibiotics and pesticides pose a severe threat to the natural balance of the ecosystem, as they either contribute to the development of multidrug resistance among pathogens or exhibit ecotoxic and mutagenic actions of a persistent character. Recently, cold atmospheric pressure plasmas (CAPPs) have emerged as promising technologies for degradation of these organic pollutants. CAPP-based technologies show eco-friendliness and potency for the removal of organic pollutants of diverse chemical formulas and different modes of action. For this reason, various types of CAPP-based systems are presented in this review and assessed in terms of their constructions, types of discharges, operating parameters, and efficiencies in the degradation of antibiotics and persistent organic pollutants. Additionally, the key role of reactive oxygen and nitrogen species (RONS) is highlighted. Moreover, optimization of the CAPP operating parameters seems crucial to effectively remove contaminants. Finally, the CAPP-related paths and technologies are further considered in terms of biological and environmental effects associated with the treatments, including changes in antibacterial properties and toxicity of the exposed solutions, as well as the potential of the CAPP-based strategies for limiting the spread of multidrug resistance.

Wildlife ecological risk assessment in the 21st century: Promising technologies to assess toxicological effects

Despite advances in toxicity testing and the development of new approach methodologies (NAMs) for hazard assessment, the ecological risk assessment (ERA) framework for terrestrial wildlife (i.e., air-breathing amphibians, reptiles, birds, and mammals) has remained unchanged for decades. While survival, growth, and reproductive endpoints derived from whole-animal toxicity tests are central to hazard assessment, nonstandard measures of biological effects at multiple levels of biological organization (e.g., molecular, cellular, tissue, organ, organism, population, community, ecosystem) have the potential to enhance the relevance of prospective and retrospective wildlife ERAs. Other factors (e.g., indirect effects of contaminants on food supplies and infectious disease processes) are influenced by toxicants at individual, population, and community levels, and need to be factored into chemically based risk assessments to enhance the "eco" component of ERAs. Regulatory and logistical challenges often relegate such nonstandard endpoints and indirect effects to postregistration evaluations of pesticides and industrial chemicals and contaminated site evaluations. While NAMs are being developed, to date, their applications in ERAs focused on wildlife have been limited. No single magic tool or model will address all uncertainties in hazard assessment. Modernizing wildlife ERAs will likely entail combinations of laboratory- and field-derived data at multiple levels of biological organization, knowledge collection solutions (e.g., systematic review, adverse outcome pathway frameworks), and inferential methods that facilitate integrations and risk estimations focused on species, populations, interspecific extrapolations, and ecosystem services modeling, with less dependence on whole-animal data and simple hazard ratios. Integr Environ Assess Manag 2024;20:725-748. © 2023 His Majesty the King in Right of Canada and The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). Reproduced with the permission of the Minister of Environment and Climate Change Canada. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Interspecies scaling of toxicity reference values in human health versus ecological risk assessments: A critical review

Risk assessments that focus on anthropogenic chemicals in environmental media-whether considering human health or ecological effects-often rely on toxicity data from experimentally studied species to estimate safe exposures for species that lack similar data. Current default extrapolation approaches used in both human health risk assessments and ecological risk assessments (ERAs) account for differences in body weight between the test organisms and the species of interest, but the two default approaches differ in important ways. Human health risk assessments currently employ a default based on body weight raised to the three-quarters power. Ecological risk assessments for wildlife (i.e., mammals and birds) are typically based directly on body weight, as measured in the test organism and receptor species. This review describes differences in the experimental data underlying these default practices and discusses the many factors that affect interspecies variability in chemical exposures. The interplay of these different factors can lead to substantial departures from default expectations. Alternative methodologies for conducting more accurate interspecies extrapolations in ERAs for wildlife are discussed, including tissue-based toxicity reference values, physiologically based toxicokinetic and/or toxicodynamic modeling, chemical read-across, and a system of categorical defaults based on route of exposure and toxic mode of action. Integr Environ Assess Manag 2024;20:749-764. © 2023 SETAC.

A novel method to derive a human safety limit for PFOA by gene expression profiling and modelling

Perfluorooctanoic acid (PFOA) is a persistent environmental contaminant that can accumulate in the human body due to its long half-life. This substance has been associated with liver, pancreatic, testicular and breast cancers, liver steatosis and endocrine disruption. PFOA is a member of a large group of substances also known as "forever chemicals" and the vast majority of substances of this group lack toxicological data that would enable their effective risk assessment in terms of human health hazards. This study aimed to derive a health-based guidance value for PFOA intake (ng/kg BW/day) from in vitro transcriptomics data. To this end, we developed an in silico workflow comprising five components: (i) sourcing in vitro hepatic transcriptomics concentration-response data; (ii) deriving molecular points of departure using BMDExpress3 and performing pathway analysis using gene set enrichment analysis (GSEA) to identify the most sensitive molecular pathways to PFOA exposure; (iii) estimating freely-dissolved PFOA concentrations in vitro using a mass balance model; (iv) estimating in vivo doses by reverse dosimetry using a PBK model for PFOA as part of a quantitative in vitro to in vivo extrapolation (QIVIVE) algorithm; and (v) calculating a tolerable daily intake (TDI) for PFOA. Fourteen percent of interrogated genes exhibited in vitro concentration-response relationships. GSEA pathway enrichment analysis revealed that "fatty acid metabolism" was the most sensitive pathway to PFOA exposure. In vitro free PFOA concentrations were calculated to be 2.9% of the nominal applied concentrations, and these free concentrations were input into the QIVIVE workflow. Exposure doses for a virtual population of 3,000 individuals were estimated, from which a TDI of 0.15 ng/kg BW/day for PFOA was calculated using the benchmark dose modelling software, PROAST. This TDI is comparable to previously published values of 1.16, 0.69, and 0.86 ng/kg BW/day by the European Food Safety Authority. In conclusion, this study demonstrates the combined utility of an "omics"-derived molecular point of departure and in silico QIVIVE workflow for setting health-based guidance values in anticipation of the acceptance of in vitro concentration-response molecular measurements in chemical risk assessment.

The effect of chronic lithium treatment on hippocampal progenitor cells: Transcriptomic analysis and systems pharmacology

OBJECTIVE

To identify the genomics underpinning the increased volume of the hippocampus after long-term administration of lithium (Li) in bipolar disorder patients, hypothesizing the possible contribution of cell growth and differentiation pathways to this complication.

METHODS

RNA-seq profiles of four samples of hippocampal progenitor cells chronically treated with a high dose of Li and three samples chronically treated with the therapeutic dose were retrieved from NCBI-GEO. The raw data underwent filtration, quality control, expression fold change, adjusted significance, functional enrichment, and pharmacogenomic analyses.

RESULTS

CCND1, LOXL2, and PRNP were identified as the genes involved in the drug response and the chronic effects of Li in the hippocampal cells. GSK-3β was also a hub in the pharmacogenomic network of Li. In addition, ZMPSTE24 and DHX35 were identified as the important genes in lithium therapy.

CONCLUSIONS

As shown by gene ontology results, these findings conclude that lithium may increase the size of the hippocampus in bipolar patients by stimulating the generation of new neurons and promoting their differentiation into neuroblasts, neurons, or microglia.

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.

Linking Transthyretin-Binding Chemicals and Free Thyroid Hormones: In Vitro to In Vivo Extrapolation Based on a Competitive Binding Model

Large numbers of pollutants competitively inhibit the binding between thyroid hormones and transthyretin (TTR) in vitro. However, the impact of this unintended binding on free thyroid hormones in vivo has not yet been characterized. Herein, we established a quantitative in vitro to in vivo extrapolation (QIVIVE) method based on a competitive binding model to quantify the effect of TTR-binding chemicals on free thyroid hormones in human blood. Twenty-five TTR-binding chemicals including 6 hydroxyl polybromodiphenyl ethers (OH-PDBEs), 6 hydroxyl polychlorobiphenyls (OH-PCBs), 4 halogenphenols, 5 per- and polyfluorinated substances (PFASs), and 4 phenols were selected for investigation. Incorporating the in vitro binding parameters and human exposure data, the QIVIVE model could well predict the in vivo effect on free thyroid hormones. Co-exposure to twenty-five typical TTR-binding chemicals resulted in median increases of 0.080 and 0.060% in circulating levels of free thyroxine (FT4) and free triiodothyronine (FT3) in the general population. Individuals with occupational exposure to TTR-binding chemicals suffered 1.88-32.2% increases in free thyroid hormone levels. This study provides a quantitative tool to evaluate the thyroid-disrupting risks of TTR-binding chemicals and proposes a new framework for assessing the in vivo effects of chemical exposures on endogenous molecules.

Rapid equilibrium dialysis, ultrafiltration or ultracentrifugation? Evaluation of methods to quantify the unbound fraction of substances in plasma

In pharmacokinetics plasma protein binding (PPB) is a well-established parameter impacting drug disposition. The unbound fraction (f_u) is arguably regarded the effective concentration at the target site. Pharmacology and toxicology, increasingly use in vitro models. The translation of in vitro concentrations to in vivo doses can be supported by toxicokinetic modelling, e.g. physiologically based toxicokinetic models (PBTK). PPB of a test substance is an input parameter for PBTK. We compared three methods to quantify f_u: rapid equilibrium dialysis (RED), ultrafiltration (UF) and ultracentrifugation (UC) using twelve substances covering a wide range of Log P_ow (-0.1 to 6.8) and molecular weights (151 and 531 g/mol): Acetaminophen, Bisphenol A, Caffeine, Colchicine, Fenarimol, Flutamide, Genistein, Ketoconazole, α-Methyltestosterone, Tamoxifen, Trenbolone and Warfarin. After RED and UF separation, three polar substances (Log P_ow < 2) were largely unbound (f_u > 70%), while more lipophilic substances were largely bound (f_u < 33%). Compared to RED or UF, UC resulted in a generally higher f_u of lipophilic substances. f_u obtained after RED and UF were more consistent with published data. For half of the substances, UC resulted in f_u higher than the reference data. UF, RED and both UF and UC, resulted in lower f_u of Flutamide, Ketoconazole and Colchicine, respectively. For f_u quantifications, the separation method should be selected according to the test substance's properties. Based on our data, RED is suitable for a broader range of substances while UC and UF are suitable for polar substances.

Closing the gap between in vivo and in vitro omics: using QA/QC to strengthen ex vivo NMR metabolomics

Metabolomics aims to achieve a global quantitation of the pool of metabolites within a biological system. Importantly, metabolite concentrations serve as a sensitive marker of both genomic and phenotypic changes in response to both internal and external stimuli. NMR spectroscopy greatly aids in the understanding of both in vitro and in vivo physiological systems and in the identification of diagnostic and therapeutic biomarkers. Accordingly, NMR is widely utilized in metabolomics and fluxomics studies due to its limited requirements for sample preparation and chromatography, its non-destructive and quantitative nature, its utility in the structural elucidation of unknown compounds, and, importantly, its versatility in the analysis of in vitro, in vivo, and ex vivo samples. This review provides an overview of the strengths and limitations of in vitro and in vivo experiments for translational research and discusses how ex vivo studies may overcome these weaknesses to facilitate the extrapolation of in vitro insights to an in vivo system. The application of NMR-based metabolomics to ex vivo samples, tissues, and biofluids can provide essential information that is close to a living system (in vivo) with sensitivity and resolution comparable to those of in vitro studies. The success of this extrapolation process is critically dependent on high-quality and reproducible data. Thus, the incorporation of robust quality assurance and quality control checks into the experimental design and execution of NMR-based metabolomics experiments will ensure the successful extrapolation of ex vivo studies to benefit translational medicine.

Prediction of in vivo prenatal chlorpyrifos exposure leading to developmental neurotoxicity in humans based on in vitro toxicity data by quantitative in vitro-in vivo extrapolation

Introduction: Epidemiological studies in children suggested that in utero exposure to chlorpyrifos (CPF), an organophosphate insecticide, may cause developmental neurotoxicity (DNT). We applied quantitative in vitro-in vivo extrapolation (QIVIVE) based on in vitro concentration and non-choline esterase-dependent effects data combined with Benchmark dose (BMD) modelling to predict oral maternal CPF exposure during pregnancy leading to fetal brain effect concentration. By comparing the results with data from epidemiological studies, we evaluated the contribution of the in vitro endpoints to the mode of action (MoA) for CPF-induced DNT. Methods: A maternal-fetal PBK model built in PK-Sim^® was used to perform QIVIVE predicting CPF concentrations in a pregnant women population at 15 weeks of gestation from cell lysate concentrations obtained in human induced pluripotent stem cell-derived neural stem cells undergoing differentiation towards neurons and glia exposed to CPF for 14 days. The in vitro concentration and effect data were used to perform BMD modelling. Results: The upper BMD was converted into maternal doses which ranged from 3.21 to 271 mg/kg bw/day. Maternal CPF blood levels from epidemiological studies reporting DNT findings in their children were used to estimate oral CPF exposure during pregnancy using the PBK model. It ranged from 0.11 to 140 μg/kg bw/day. Discussion: The effective daily intake doses predicted from the in vitro model were several orders of magnitude higher than exposures estimated from epidemiological studies to induce developmental non-cholinergic neurotoxic responses, which were captured by the analyzed in vitro test battery. These were also higher than the in vivo LOEC for cholinergic effects. Therefore, the quantitative predictive value of the investigated non-choline esterase-dependent effects, although possibly relevant for other chemicals, may not adequately represent potential key events in the MoA for CPF-associated DNT.

A Physiologically Based Pharmacokinetic (PBPK) Modeling Framework for Mixtures of Dioxin-like Compounds

Humans are exposed to persistent organic pollutants, such as dioxin-like compounds (DLCs), as mixtures. Understanding and predicting the toxicokinetics and thus internal burden of major constituents of a DLC mixture is important for assessing their contributions to health risks. PBPK models, including dioxin models, traditionally focus on one or a small number of compounds; developing new or extending existing models for mixtures often requires tedious, error-prone coding work. This lack of efficiency to scale up for multi-compound exposures is a major technical barrier toward large-scale mixture PBPK simulations. Congeners in the DLC family, including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), share similar albeit quantitatively different toxicokinetic and toxicodynamic properties. Taking advantage of these similarities, here we reported the development of a human PBPK modeling framework for DLC mixtures that can flexibly accommodate an arbitrary number of congeners. Adapted from existing TCDD models, our mixture model contains the blood and three diffusion-limited compartments-liver, fat, and rest of the body. Depending on the number of congeners in a mixture, varying-length vectors of ordinary differential equations (ODEs) are automatically generated to track the tissue concentrations of the congeners. Shared ODEs are used to account for common variables, including the aryl hydrocarbon receptor (AHR) and CYP1A2, to which the congeners compete for binding. Binary and multi-congener mixture simulations showed that the AHR-mediated cross-induction of CYP1A2 accelerates the sequestration and metabolism of DLC congeners, resulting in consistently lower tissue burdens than in single exposure, except for the liver. Using dietary intake data to simulate lifetime exposures to DLC mixtures, the model demonstrated that the relative contributions of individual congeners to blood or tissue toxic equivalency (TEQ) values are markedly different than those to intake TEQ. In summary, we developed a mixture PBPK modeling framework for DLCs that may be utilized upon further improvement as a quantitative tool to estimate tissue dosimetry and health risks of DLC mixtures.

The potential toxicity of polystyrene nanoplastics to human trophoblasts in vitro

Nanoplastics (NPs), the emerging contaminants in recent years, widely distributed in the environment and are bioaccumulated and biomagnified in organisms through food chain. A growing number of studies have detected plastic particulates in human placenta and blood. However, few studies have focused on their effects during human pregnancy. Herein, human trophoblast HTR-8/Svneo cells were used to evaluate the effects and the possible mechanism of 100-nm polystyrene NPs on placental trophoblasts at the maternal-fetal interface. The results showed that NPs entered the trophoblastic cytoplasm, decreased cell viability, caused cell cycle arrest, reduced the cell migration and invasion abilities, increased level of intracellular reactive oxygen species and the production of proinflammatory cytokines (TNF-α and IFN-γ) in a dose-dependent manner. Furthermore, global transcriptome sequencing (RNA-Seq) was performed on HTR-8/Svneo cells with or without 100 μg/mL PS-NP exposure for 24 h. A total of 344 differentially expressed genes were detected. The gene functions for regulation of leukocyte differentiation, response to stimulus, cell cycle, apoptotic process, and cell adhesion were enriched. Thyroid hormone, Hippo, TGF-β and FoxO signaling pathways were activated. Collectively, our data provided evidences for the adverse consequences of NPs on the biological functions of trophoblasts, which provided new insights into the potential trophoblast toxicity of NPs in mammals.

NAM-based Prediction of Point-of-contact Toxicity in the Lung: A Case Example With 1,3-dichloropropene

Time, cost, ethical, and regulatory considerations surrounding in vivo testing methods render them insufficient to meet existing and future chemical safety testing demands. There is a need for the development of in vitro and in silico alternatives to replace traditional in vivo methods for inhalation toxicity assessment. Exposures of differentiated airway epithelial cultures to gases or aerosols at the air-liquid interface (ALI) can assess tissue responses and in vitro to in vivo extrapolation can align in vitro exposure levels with in-life exposures expected to give similar tissue exposures. Because the airway epithelium varies along its length, with various regions composed of different cell types, we have introduced a known toxic vapor to five human-derived, differentiated, in vitro airway epithelial cell culture models-MucilAir of nasal, tracheal, or bronchial origin, SmallAir, and EpiAlveolar-representing five regions of the airway epithelium-nasal, tracheal, bronchial, bronchiolar, and alveolar. We have monitored toxicity in these cultures 24hours after acute exposure using an assay for transepithelial conductance (for epithelial barrier integrity) and the lactate dehydrogenase (LDH) release assay (for cytotoxicity). Our vapor of choice in these experiments was 1,3-dichloropropene (1,3-DCP). Finally, we have developed an airway dosimetry model for 1,3-DCP vapor to predict in vivo external exposure scenarios that would produce toxic local tissue concentrations as determined by in vitro experiments. Measured in vitro points of departure (PoDs) for all tested cell culture models were similar. Calculated rat equivalent inhaled concentrations varied by model according to position of the modeled tissue within the airway, with nasal respiratory tissue being the most proximal and most sensitive tissue, and alveolar epithelium being the most distal and least sensitive tissue. These predictions are qualitatively in accordance with empirically determined in vivo PoDs. The predicted PoD concentrations were close to, but slightly higher than, PoDs determined by in vivo subchronic studies.

Towards best use and regulatory acceptance of generic physiologically based kinetic (PBK) models for in vitro-to-in vivo extrapolation (IVIVE) in chemical risk assessment

With an increasing need to incorporate new approach methodologies (NAMs) in chemical risk assessment and the concomitant need to phase out animal testing, the interpretation of in vitro assay readouts for quantitative hazard characterisation becomes more important. Physiologically based kinetic (PBK) models, which simulate the fate of chemicals in tissues of the body, play an essential role in extrapolating in vitro effect concentrations to in vivo bioequivalent exposures. As PBK-based testing approaches evolve, it will become essential to standardise PBK modelling approaches towards a consensus approach that can be used in quantitative in vitro-to-in vivo extrapolation (QIVIVE) studies for regulatory chemical risk assessment based on in vitro assays. Based on results of an ECETOC expert workshop, steps are recommended that can improve regulatory adoption: (1) define context and implementation, taking into consideration model complexity for building fit-for-purpose PBK models, (2) harmonise physiological input parameters and their distribution and define criteria for quality chemical-specific parameters, especially in the absence of in vivo data, (3) apply Good Modelling Practices (GMP) to achieve transparency and design a stepwise approach for PBK model development for risk assessors, (4) evaluate model predictions using alternatives to in vivo PK data including read-across approaches, (5) use case studies to facilitate discussions between modellers and regulators of chemical risk assessment. Proof-of-concepts of generic PBK modelling approaches are published in the scientific literature at an increasing rate. Working on the previously proposed steps is, therefore, needed to gain confidence in PBK modelling approaches for regulatory use.

Enhanced antitumor activity of combined methotrexate and histone deacetylase inhibitor valproic acid on mammary cancer in vitro and in vivo

Histone deacetylase inhibitors (HDACIs) act as antiproliferative agents by promoting differentiation and inducing apoptosis. Valproic acid (VPA) is an HDACI that shows promising chemotherapeutic effect in several tumor cells. The present study aimed to investigate the inhibitory effect of VPA on the viability of mammary cancer cells and its enhancing effect with methotrexate (MTX) in vitro and in vivo. Treatment with VPA or MTX alone induced concentration-dependent cytotoxic effects in two breast cancer cell lines. VPA significantly increased the cytotoxicity of MTX 3 times against MCF7. VPA addition to MTX, however, did not produce any significant changes on MTX cytotoxicity against MDA-MB231. VPA (150 and 200 mg/kg) significantly inhibited the growth of IP and SC Ehrlich ascites carcinoma tumor mouse models and improved results were achieved for tumor inhibition when VPA was combined with MTX (1 and 2 mg/kg) in vivo. The antitumor activity was not associated with a significant increase in toxicity or mice mortality rate. All these findings suggest that the combination of MTX and VPA may have clinical and/or adjuvant therapeutic application in the treatment of mammary cancer.

Cigarette Smoke Extract, but Not Electronic Cigarette Aerosol Extract, Inhibits Monoamine Oxidase in vitro and Produces Greater Acute Aversive/Anhedonic Effects Than Nicotine Alone on Intracranial Self-Stimulation in Rats

Conventional tobacco cigarettes appear to have greater abuse liability than non-combusted products such as electronic cigarettes (ECs) and nicotine replacement therapy (NRT). This may be due to the higher levels of behaviorally active non-nicotine constituents [e.g., monoamine oxidase (MAO) inhibitors such as β-carbolines] in cigarette smoke (CS) compared to non-combusted products. To evaluate this hypothesis, the current studies compared the relative abuse liability of CS and EC aerosol extracts containing nicotine and a range of non-nicotine constituents to that of nicotine alone (NRT analog) using intracranial self-stimulation (ICSS) in rats. Effects of formulations on brain MAO activity in vitro and ex vivo were also studied to evaluate the potential role of MAO inhibition in the ICSS study. CS extract contained higher levels of several behaviorally active non-nicotine constituents (e.g., the β-carbolines norharmane and harmane) than EC extract. Nicotine alone reduced ICSS thresholds at a moderate nicotine dose, suggesting a reinforcement-enhancing effect that may promote abuse liability, and elevated ICSS thresholds at a high nicotine dose, suggesting an aversive/anhedonic effect that may limit abuse liability. CS extract elevated ICSS thresholds to a greater degree than nicotine alone at high nicotine doses. Effects of EC extract on ICSS did not differ from those of nicotine alone. Finally, CS extract significantly inhibited MAO-A and MAO-B activity in vitro, whereas EC extract and nicotine alone did not. None of the formulations inhibited MAO measured ex vivo. These findings indicate greater acute aversive/anhedonic effects for CS extract compared to nicotine alone, suggesting lower abuse liability. Although confirmation of our findings using other dosing regimens, preclinical addiction models, and tobacco product extracts is needed, these findings suggest that the centrally-mediated effects of MAO inhibitors and other non-nicotine constituents may not account for the greater abuse liability of cigarettes compared to non-combusted products. Nonetheless, identifying the specific constituent(s) mediating the effects of CS extracts in this study could help clarify mechanisms mediating tobacco addiction and inform FDA product standards.

In vitro to in vivo extrapolation to support the development of the next generation risk assessment (NGRA) strategy for nanomaterials

There is growing interest in developing novel strategies to support assessment of human health risks due to chemicals. Regulatory and decision-making agencies have recommended that non-animal-based alternatives should be applied whenever possible instead of experimentation on living animals. These alternative methods are beneficial because they are ethical, inexpensive, and rapid. Herein, we review recent activities aimed at developing in vitro to in vivo extrapolation (IVIVE) models as a part of the Next Generation Risk Assessment (NGRA) of nanomaterials. In this context, we show the adverse outcome pathway (AOP)-based methodology for the identification of mechanistically relevant events serving as biomarkers for the targeted selection of in vitro assays. Considered events need to be biologically plausible, regulatory relevant, and crucial for the examination of occurrence of adverse outcomes. The promising advantages of using high-throughout-based omics data are highlighted. Furthermore, the application of 3D in vitro models and nano genome atlases to study nanoparticle toxicity is briefly summarized. Additionally, the challenges related to the extrapolation of in vitro doses into in vivo-relevant responses are presented. We also discuss the limitations of models applied thus far to study the fate of chemicals in the human body, which exist due to the lack of available knowledge regarding transformations of nanomaterials occurring in biological systems.

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.

Investigation of the Associations between a Nanomaterial's Microrheology and Toxicology

With the advent of Nanotechnology, the use of nanomaterials in consumer products is increasing on a daily basis, due to which a deep understanding and proper investigation regarding their safety and risk assessment should be a major priority. To date, there is no investigation regarding the microrheological properties of nanomaterials (NMs) in biological media. In our study, we utilized in silico models to select the suitable NMs based on their physicochemical properties such as solubility and lipophilicity. Then, we established a new method based on dynamic light scattering (DLS) microrheology to get the mean square displacement (MSD) and viscoelastic property of two model NMs that are dendrimers and cerium dioxide nanoparticles in Dulbecco's Modified Eagle Medium (DMEM) complete media at three different concentrations for both NMs. Subsequently, we established the cytotoxicological profiling using water-soluble tetrazolium salt-1 (WST-1) and a reactive oxygen species (ROS) assay. To take one step forward, we further looked into the tight junction properties of the cells using immunostaining with Zonula occluden-1 (ZO-1) antibodies and found that the tight junction function or transepithelial resistance (TEER) was affected in response to the microrheology and cytotoxicity. The quantitative polymerase chain reaction (q-PCR) results in the gene expression of ZO-1 after the 24 h treatment with NPs further validates the findings of immunostaining results. This new method that we established will be a reference point for other NM studies which are used in our day-to-day consumer products.

Application of an Accessible Interface for Pharmacokinetic Modeling and In Vitro to In Vivo Extrapolation

Regulatory toxicology testing has traditionally relied on in vivo methods to inform decision-making. However, scientific, practical, and ethical considerations have led to an increased interest in the use of in vitro and in silico methods to fill data gaps. While in vitro experiments have the advantage of rapid application across large chemical sets, interpretation of data coming from these non-animal methods can be challenging due to the mechanistic nature of many assays. In vitro to in vivo extrapolation (IVIVE) has emerged as a computational tool to help facilitate this task. Specifically, IVIVE uses physiologically based pharmacokinetic (PBPK) models to estimate tissue-level chemical concentrations based on various dosing parameters. This approach is used to estimate the administered dose needed to achieve in vitro bioactivity concentrations within the body. IVIVE results can be useful to inform on metrics such as margin of exposure or to prioritize potential chemicals of concern, but the PBPK models used in this approach have extensive data requirements. Thus, access to input parameters, as well as the technical requirements of applying and interpreting models, has limited the use of IVIVE as a routine part of in vitro testing. As interest in using non-animal methods for regulatory and research contexts continues to grow, our perspective is that access to computational support tools for PBPK modeling and IVIVE will be essential for facilitating broader application and acceptance of these techniques, as well as for encouraging the most scientifically sound interpretation of in vitro results. We highlight recent developments in two open-access computational support tools for PBPK modeling and IVIVE accessible via the Integrated Chemical Environment (https://ice.ntp.niehs.nih.gov/), demonstrate the types of insights these tools can provide, and discuss how these analyses may inform in vitro-based decision making.

Metabolomic Approaches for Detection and Identification of Biomarkers and Altered Pathways in Bladder Cancer

Metabolomic analysis has proven to be a useful tool in biomarker discovery and the molecular classification of cancers. In order to find new biomarkers, and to better understand its pathological behavior, bladder cancer also has been studied using a metabolomics approach. In this article, we review the literature on metabolomic studies of bladder cancer, focusing on the different available samples (urine, blood, tissue samples) used to perform the studies and their relative findings. Moreover, the multi-omic approach in bladder cancer research has found novel insights into its metabolic behavior, providing excellent start-points for new diagnostic and therapeutic strategies. Metabolomics data analysis can lead to the discovery of a “signature pathway” associated with the progression of bladder cancer; this aspect could be potentially valuable in predictions of clinical outcomes and the introduction of new treatments. However, further studies are needed to give stronger evidence and to make these tools feasible for use in clinical practice.

PBK modelling of topical application and characterisation of the uncertainty of Cmax estimate: A case study approach

Combined with in vitro bioactivity data, physiologically based kinetic (PBK) models has increasing applications in next generation risk assessment for animal-free safety decision making. A tiered framework of building PBK models for such application has been developed with increasing complexity and refinements, as model parameters determined in silico, in vitro, and with human pharmacokinetic data become progressively available. PBK modelling has been widely applied for oral/intravenous administration, but less so on topically applied chemicals. Therefore, building PBK models for topical applications and characterizing their uncertainties in the tiered approach is critical to safety decision making. The purpose of this study was to assess the confidence of PBK modelling of topically applied chemicals following the tiered framework, using non-animal methods derived parameters. Prediction of maximum plasma concentration (C_max) and area under the curve were compared to observed kinetics from published dermal clinical studies for five chemicals (diclofenac, salicylic acid, coumarin, nicotine, caffeine). A bespoke Bayesian statistical model was developed to describe the distributions of C_max errors between the predicted and observed data. We showed a general trend that confidence in model predictions increases when more quality in vitro data, particularly those on hepatic clearance and dermal absorption, are available as model input. The overall fold error distributions are useful for characterizing model uncertainty. We concluded that by identifying and quantifying the uncertainties in the tiered approach, we can increase the confidence in using PBK modelling to help make safety decisions on topically applied chemicals in the absence of human pharmacokinetic data.

Extension of the Virtual Cell Based Assay from a 2-D to a 3-D Cell Culture Model

Prediction of chemical toxicity is very useful in risk assessment. With the current paradigm shift towards the use of in vitro and in silico systems, we present herein a theoretical mathematical description of a quasi-diffusion process to predict chemical concentrations in 3-D spheroid cell cultures. By extending a 2-D Virtual Cell Based Assay (VCBA) model into a 3-D spheroid cell model, we assume that cells are arranged in a series of concentric layers within the sphere. We formulate the chemical quasi-diffusion process by simplifying the spheroid with respect to the number of cells in each layer. The system was calibrated and tested with acetaminophen (APAP). Simulated predictions of APAP toxicity were compared with empirical data from in vitro measurements by using a 3-D spheroid model. The results of this first attempt to extend the VCBA model are promising - they show that the VCBA model simulates close correlation between the influence of compound concentration and the viability of the HepaRG 3-D cell culture. The 3-D VCBA model provides a complement to current in vitro procedures to refine experimental setups, to fill data gaps and help in the interpretation of in vitro data for the purposes of risk assessment.

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

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

Use of Physiologically-Based Kinetics Modelling to Reliably Predict Internal Concentrations of the UV Filter, Homosalate, After Repeated Oral and Topical Application

Ethical and legal considerations have led to increased use of non-animal methods to evaluate the safety of chemicals for human use. We describe the development and qualification of a physiologically-based kinetics (PBK) model for the cosmetic UV filter ingredient, homosalate, to support its safety without the need of generating further animal data. The intravenous (IV) rat PBK model, using PK-Sim®, was developed and validated using legacy in vivo data generated prior to the 2013 EU animal-testing ban. Input data included literature or predicted physicochemical and pharmacokinetic properties. The refined IV rat PBK model was subject to sensitivity analysis to identify homosalate-specific sensitive parameters impacting the prediction of Cmax (more sensitive than AUC(0-∞)). These were then considered, together with population modeling, to calculate the confidence interval (CI) 95% Cmax and AUC(0-∞). Final model parameters were established by visual inspection of the simulations and biological plausibility. The IV rat model was extrapolated to oral administration, and used to estimate internal exposures to doses tested in an oral repeated dose toxicity study. Next, a human PBK dermal model was developed using measured human in vitro ADME data and a module to represent the dermal route. Model performance was confirmed by comparing predicted and measured values from a US-FDA clinical trial (Identifier: NCT03582215, https://clinicaltrials.gov/). Final exposure estimations were obtained in a virtual population and considering the in vitro and input parameter uncertainty. This model was then used to estimate the Cmax and AUC(0-24 h) of homosalate according to consumer use in a sunscreen. The developed rat and human PBK models had a good biological basis and reproduced in vivo legacy rat and human clinical kinetics data. They also complied with the most recent WHO and OECD recommendations for assessing the confidence level. In conclusion, we have developed a PBK model which predicted reasonably well the internal exposure of homosalate according to different exposure scenarios with a medium to high level of confidence. In the absence of in vivo data, such human PBK models will be the heart of future completely non-animal risk assessments; therefore, valid approaches will be key in gaining their regulatory acceptance. Clinical Trial Registration: https://clinicaltrials.gov/, identifier, NCT03582215.

(Eco)toxicological tests for assessing impacts of chemical stress to aquatic ecosystems: Facts, challenges, and future

Monitoring of chemicals in the aquatic environment by chemical analysis alone cannot completely assess and predict the effects of chemicals on aquatic species and ecosystems. This is primarily because of the increasing number of (unknown) chemical stressors and mixture effects present in the environment. In addition, the ability of ecological indices to identify underlying stressors causing negative ecological effects is limited. Therefore, additional complementary methods are needed that can address the biological effects in a direct manner and provide a link to chemical exposure, i.e. (eco)toxicological tests. (Eco)toxicological tests are defined as test systems that expose biological components (cells, individuals, populations, communities) to (environmental mixtures of) chemicals to register biological effects. These tests measure responses at the sub-organismal (biomarkers and in vitro bioassays), whole-organismal, population, or community level. We performed a literature search to obtain a state-of-the-art overview of ecotoxicological tests available for assessing impacts of chemicals to aquatic biota and to reveal datagaps. In total, we included 509 biomarkers, 207 in vitro bioassays, 422 tests measuring biological effects at the whole-organismal level, and 78 tests at the population- community- and ecosystem-level. Tests at the whole-organismal level and biomarkers were most abundant for invertebrates and fish, whilst in vitro bioassays are mostly based on mammalian cell lines. Tests at the community- and ecosystem-level were almost missing for organisms other than microorganisms and algae. In addition, we provide an overview of the various extrapolation challenges faced in using data from these tests and suggest some forward looking perspectives. Although extrapolating the measured responses to relevant protection goals remains challenging, the combination of ecotoxicological experiments and models is key for a more comprehensive assessment of the effects of chemical stressors to aquatic ecosystems.

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.

Lithium treatment and human hippocampal neurogenesis

Lithium is a first-line treatment for bipolar disorder, where it acts as a mood-stabilizing agent. Although its precise mechanism remains unclear, neuroimaging studies have shown that lithium accumulates in the hippocampus and that chronic use amongst bipolar disorder patients is associated with larger hippocampal volumes. Here, we tested the chronic effects of low (0.75 mM) and high (2.25 mM) doses of lithium on human hippocampal progenitor cells and used immunocytochemistry to investigate the effects of lithium on cell parameters implicated in neurogenesis. Corresponding RNA-sequencing and gene-set enrichment analyses were used to evaluate whether genes affected by lithium in our model overlap with those regulating the volume of specific layers of the dentate gyrus. We observed that high-dose lithium treatment in human hippocampal progenitors increased the generation of neuroblasts (P ≤ 0.01), neurons (P ≤ 0.01), and glia (P ≤ 0.001), alongside the expression of genes, which regulate the volume of the molecular layer of the dentate gyrus. This study provides empirical support that adult hippocampal neurogenesis and gliogenesis are mechanisms that could contribute to the effects of lithium on human hippocampal volume.

Opportunities and challenges related to saturation of toxicokinetic processes: Implications for risk assessment

Top dose selection for repeated dose animal studies has generally focused on identification of apical endpoints, use of the limit dose, or determination of a maximum tolerated dose (MTD). The intent is to optimize the ability of toxicity tests performed in a small number of animals to detect effects for hazard identification. An alternative approach, the kinetically derived maximum dose (KMD), has been proposed as a mechanism to integrate toxicokinetic (TK) data into the dose selection process. The approach refers to the dose above which the systemic exposures depart from being proportional to external doses. This non-linear external-internal dose relationship arises from saturation or limitation of TK process(es), such as absorption or metabolism. The importance of TK information is widely acknowledged when assessing human health risks arising from exposures to environmental chemicals, as TK determines the amount of chemical at potential sites of toxicological responses. However, there have been differing opinions and interpretations within the scientific and regulatory communities related to the validity and application of the KMD concept. A multi-stakeholder working group, led by the Health and Environmental Sciences Institute (HESI), was formed to provide an opportunity for impacted stakeholders to address commonly raised scientific and technical issues related to this topic and, more specifically, a weight of evidence approach is recommended to inform design and dose selection for repeated dose animal studies. Commonly raised challenges related to the use of TK data for dose selection are discussed, recommendations are provided, and illustrative case examples are provided to address these challenges or refute misconceptions.

An integrative approach to define chemical exposure threshold limits for endangered sea turtles

Environmental contaminants pose serious health threats to marine megafauna species, yet methods defining exposure threshold limits are lacking. Here, a three-pillar chemical risk assessment framework is presented based on (1) species- and chemical-specific lifetime bioaccumulation modelling, (2) non-destructive in vitro and in vivo toxicity threshold assessment, and (3) chemical risk quantification. We used the effects of cadmium (Cd) in green sea turtles (Chelonia mydas) as a proof of concept to evaluate the quantitative mechanistic modelling approach. A physiologically-based kinetic (PBK) model simulated Cd tissue concentrations (liver, kidney, muscle, fat, brain, scute, and 'rest of the body') in C.mydas. The validated PBK model then translated species-specific in vitro results to in vivo effects. The results showed that the resilience of C.mydas towards Cd kidney toxicity is age-dependent and differs with changing physiology and feeding ecology. Using the model in reverse mode, a steady-state exposure threshold of 0.1 µg/g dry weight Cd in forage was derived and compared to real-world exposure scenarios. Three out of the four globally distinct C.mydas populations assessed are exposed to Cd levels above this threshold limit. This approach can be adapted to other marine species and chemicals to prioritize measures for managing potentially harmful chemical exposures.

Quantitative in Vitro-to-in Vivo Extrapolation for Mixtures: A Case Study of Superfund Priority List Pesticides

In vitro cell-based toxicity testing methods generate large amounts of data informative for risk-based evaluations. To allow extrapolation of the quantitative outputs from cell-based tests to the equivalent exposure levels in humans, reverse toxicokinetic (RTK) modeling is used to conduct in vitro-to-in vivo extrapolation (IVIVE) from in vitro effective concentrations to in vivo oral dose equivalents. IVIVE modeling approaches for individual chemicals are well-established; however, the potential implications of chemical-to-chemical interactions in mixture settings on IVIVE remains largely unexplored. We hypothesized that chemical co-exposures could modulate both protein binding efficiency and hepatocyte clearance of the chemicals in a mixture, which would in turn affect the quantitative IVIVE toxicokinetic parameters. To test this hypothesis, we used 20 pesticides from the Agency for Toxic Substances and Disease Registry (ATSDR) Substance Priority List, both individually and as equimolar mixtures, and investigated the concentration-dependent effects of chemical interactions on in vitro toxicokinetic parameters. Plasma protein binding efficiency was determined by using ultracentrifugation, and hepatocyte clearance was estimated in suspensions of cryopreserved primary human hepatocytes. We found that for single chemicals, the protein binding efficiencies were similar at different test concentrations. In a mixture, however, both protein binding efficiency and hepatocyte clearance were affected. When IVIVE was conducted using mixture-derived toxicokinetic data, more conservative estimates of Activity-to-Exposure Ratios (AERs) were produced as compared to using data from single chemical experiments. Because humans are exposed to mixtures of chemicals, this study is significant as it demonstrates the importance of incorporating mixture-derived parameters into IVIVE for in vitro bioactivity data in order to accurately prioritize risks and facilitate science-based decision-making.

Current EU regulatory requirements for the assessment of chemicals and cosmetic products: challenges and opportunities for introducing new approach methodologies

The EU Directive 2010/63/EU   on the protection of animals used for scientific purposes and other EU regulations, such as REACH and the Cosmetic Products Regulation advocate for a change in the way toxicity testing is conducted. Whilst the Cosmetic Products Regulation bans animal testing altogether, REACH aims for a progressive shift from in vivo testing towards quantitative in vitro and computational approaches. Several endpoints can already be addressed using non-animal approaches including skin corrosion and irritation, serious eye damage and irritation, skin sensitisation, and mutagenicity and genotoxicity. However, for systemic effects such as acute toxicity, repeated dose toxicity and reproductive and developmental toxicity, evaluation of chemicals under REACH still heavily relies on animal tests. Here we summarise current EU regulatory requirements for the human health assessment of chemicals under REACH and the Cosmetic Products Regulation, considering the more critical endpoints and identifying the main challenges in introducing alternative methods into regulatory testing practice. This supports a recent initiative taken by the International Cooperation on Alternative Test Methods (ICATM) to summarise current regulatory requirements specific for the assessment of chemicals and cosmetic products for several human health-related endpoints, with the aim of comparing different jurisdictions and coordinating the promotion and ultimately the implementation of non-animal approaches worldwide. Recent initiatives undertaken at European level to promote the 3Rs and the use of alternative methods in current regulatory practice are also discussed.

Development and Application of a Life-Stage Physiologically Based Pharmacokinetic (PBPK) Model to the Assessment of Internal Dose of Pyrethroids in Humans

To address concerns around age-related sensitivity to pyrethroids, a life-stage physiologically based pharmacokinetic (PBPK) model, supported by in vitro to in vivo extrapolation (IVIVE) was developed. The model was used to predict age-dependent changes in target tissue exposure of 8 pyrethroids; deltamethrin (DLM), cis-permethrin (CPM), trans-permethrin, esfenvalerate, cyphenothrin, cyhalothrin, cyfluthrin, and bifenthrin. A single model structure was used based on previous work in the rat. Intrinsic clearance (CLint) of each individual cytochrome P450 or carboxylesterase (CES) enzyme that are active for a given pyrethroid were measured in vitro, then biologically scaled to obtain in vivo age-specific total hepatic CLint. These IVIVE results indicate that, except for bifenthrin, CES enzymes are largely responsible for human hepatic metabolism (>50% contribution). Given the high efficiency and rapid maturation of CESs, clearance of the pyrethroids is very efficient across ages, leading to a blood flow-limited metabolism. Together with age-specific physiological parameters, in particular liver blood flow, the efficient metabolic clearance of pyrethroids across ages results in comparable to or even lower internal exposure in the target tissue (brain) in children than that in adults in response to the same level of exposure to a given pyrethroid (Cmax ratio in brain between 1- and 25-year old = 0.69, 0.93, and 0.94 for DLM, bifenthrin, and CPM, respectively). Our study demonstrated that a life-stage PBPK modeling approach, coupled with IVIVE, provides a robust framework for evaluating age-related differences in pharmacokinetics and internal target tissue exposure in humans for the pyrethroid class of chemicals.

Evaluation and comparison of in vitro intrinsic clearance rates measured using cryopreserved hepatocytes from humans, rats, and rainbow trout

In vitro and in silico methods that can reduce the need for animal testing are being used with increasing frequency to assess chemical risks to human health and the environment. The rate of hepatic biotransformation is an important species-specific parameter for determining bioaccumulation potential and extrapolating in vitro bioactivity to in vivo effects. One approach to estimating hepatic biotransformation is to employ in vitro systems derived from liver tissue to measure chemical (substrate) depletion over time which can then be translated to a rate of intrinsic clearance (CL_int). In the present study, cryopreserved hepatocytes from humans, rats, and rainbow trout were used to measure CL_int values for 54 industrial and pesticidal chemicals at starting test concentrations of 0.1 and 1 μM. A data evaluation framework that emphasizes the behavior of Heat-Treated Controls (HTC) was developed to identify datasets suitable for rate reporting. Measured or estimated ("greater than" or "less than") CL_int values were determined for 124 of 226 (55 %) species-chemical-substrate concentration datasets with acceptable analytical chemistry. A large percentage of tested chemicals exhibited low HTC recovery values, indicating a substantial abiotic loss of test chemical over time. An evaluation of K_OW values for individual chemicals suggested that in vitro test performance declined with increasing chemical hydrophobicity, although differences in testing devices for mammals and fish also likely played a role. The current findings emphasize the value of negative controls as part of a rigorous approach to data quality assessment for in vitro substrate depletion studies. Changes in current testing protocols can be expected to result in the collection of higher quality data. However, poorly soluble chemicals are likely to remain a challenge for CL_int determination.

Next generation risk assessment (NGRA): Bridging in vitro points-of-departure to human safety assessment using physiologically-based kinetic (PBK) modelling - A case study of doxorubicin with dose metrics considerations

Using the chemical doxorubicin (DOX), the objective of the present study was to evaluate the impact of dose metrics selection in the new approach method of integrating physiologically-based kinetic (PBK) modelling and relevant human cell-based assays to inform a priori the point of departure for human health risk. We reviewed the literature on the clinical consequences of DOX treatment to identify dosing scenarios with no or mild cardiotoxicity observed. Key concentrations of DOX that induced cardiomyocyte toxicity in vitro were derived from studies of our own and others. A human population-based PBK model of DOX was developed and verified against pharmacokinetic data. The model was then used to predict plasma and extracellular and intracellular heart concentrations of DOX under selected clinical settings and compared with in vitro outcomes, based on several dose metrics: C_max (maximum concentration) or AUC (area under concentration-time curve) in free or total form of DOX. We found when using in vitro assays to predict cardiotoxicity for DOX, AUC is a better indicator. Our study illustrates that when appropriate dose metrics are used, it is possible to combine PBK modelling with in vitro-derived toxicity information to define margins of safety and predict low-risk human exposure levels.

Quantitative In Vitro-to-In Vivo Extrapolation: Nominal versus Freely Dissolved Concentration

Discussions are ongoing on which dose metric should be used for quantitative in vitro-to-in vivo extrapolation (QIVIVE) of in vitro bioassay data. The nominal concentration of the test chemicals is most commonly used and easily accessible, while the concentration freely dissolved in the assay medium is considered to better reflect the bioavailable concentration but is tedious to measure. The aim of this study was to elucidate how much QIVIVE results will differ when using either nominal or freely dissolved concentrations. QIVIVE_nom and QIVIVE_free ratios, that is, the ratios of plasma concentrations divided by in vitro effect concentrations, were calculated for 10 pharmaceuticals using previously published nominal and freely dissolved effect concentrations for the activation of the peroxisome proliferator-activated receptor gamma (PPARγ) and the activation of oxidative stress response. The QIVIVE_nom ratios were higher than QIVIVE_free ratios by up to a factor of 60. The risk of in vivo effects was classified as being high or low for four chemicals using the QIVIVE_nom and for three chemicals using QIVIVE_free ratios. Unambiguous classification was possible for nine chemicals by combining the QIVIVE_nom or QIVIVE_free ratios with the respective specificity ratios (SR_nom or SR_free) of the in vitro effect data, which helps to identify whether the specific effect was influenced by cytotoxicity. QIVIVE_free models should be preferred as they account for differences in bioavailability between in vitro and in vivo, but QIVIVE_nom may still be useful for screening the effects of large numbers of chemicals because it is generally more conservative. The use of SR of the in vitro effect data as a second classification factor is recommended for QIVIVE_nom and QIVIVE_free models because a clearer picture can be obtained with respect to the likelihood that a biological effect will occur and that it is not caused by nonspecific cytotoxicity.

Development of a Web-Based Toolbox to Support Quantitative In-Vitro-to-In-Vivo Extrapolations (QIVIVE) within Nonanimal Testing Strategies

The goal of the present study was to develop an online web-based toolbox that contains generic physiologically based kinetic (PBK) models for rats and humans, including underlying calculation tools to predict plasma protein binding and tissue:plasma distribution, to be used for quantitative in-vitro-to-in-vivo extrapolations (QIVIVE). The PBK models within the toolbox allow first estimations of internal plasma and tissue concentrations of chemicals to be made, based on the logP and pKa of the chemicals and values for intestinal uptake and intrinsic hepatic clearance. As a case study, the toolbox was used to predict oral equivalent doses of in vitro ToxCast bioactivity data for the food additives methylparaben, propyl gallate, octyl gallate, and dodecyl gallate. These oral equivalent doses were subsequently compared with human exposure estimates, as a low tier assessment allowing prioritization for further assessment. The results revealed that daily intake levels of especially propyl gallate can lead to internal plasma concentrations that are close to in vitro biological effect concentrations, particularly with respect to the inhibition of human thyroid peroxidase (TPO). Estrogenic effects were not considered likely to be induced by the food additives, as daily exposure levels of the different compounds remained 2 orders of magnitude below the oral equivalent doses for in vitro estrogen receptor activation. Overall, the results of the study show how the toolbox, which is freely accessible through www.qivivetools.wur.nl, can be used to obtain initial internal dose estimates of chemicals and to prioritize chemicals for further assessment, based on the comparison of oral equivalent doses of in vitro biological activity data with human exposure levels.

Modeling Bioavailable Concentrations in Zebrafish Cell Lines and Embryos Increases the Correlation of Toxicity Potencies across Test Systems

Linking cellular toxicity to low-tier animal toxicity and beyond is crucial within the adverse outcome pathway concept and the 3R framework. This study aimed to determine and compare the bioavailable effect concentrations in zebrafish cell lines and embryos. Acute, short-term toxicity (48 h) of eight veterinary pharmaceuticals was measured in two zebrafish cell lines (hepatocytes, fibroblasts) and zebrafish embryos. Seven endpoints of cytotoxicity were recorded. The fish embryo acute toxicity test was modified by adding sublethal endpoints. Chemical distribution modeling (mass balance) was applied to compute the bioavailable compound concentrations in cells (Cfree) and embryos (Cint;aq) based on nominal effect concentrations (Cnom). Effect concentration ratios were calculated (cell effects/embryo effects). A low correlation was observed between cytotoxicity and embryo toxicity when nominal concentrations were used. Modeled bioavailable effect concentrations strongly increased correlations and placed regression lines close to the line of unity and axis origin. Cytotoxicity endpoints showed differences in sensitivity and predictability. The hepatocyte cell line depicted closer proximity to the embryo data. Conclusively, the high positive correlation between the cell- and embryo-based test systems emphasizes the appropriate modulation of toxicity when linked to bioavailable concentrations. Furthermore, it highlights the potential of fish cell lines to be utilized in integrated testing strategies.

Optimization of a pre-metabolization procedure using rat liver S9 and cell-extracted S9 in the Ames fluctuation test

Many environmental pollutants pose a toxicological hazard only after metabolic activation. In vitro bioassays using cell lines or bacteria have often no or reduced metabolic activity, which impedes their use in the risk assessment. To improve the predictive capability of in vitro assays, external metabolization systems like the liver S9 fraction are frequently combined with in vitro toxicity assays. While it is typical for S9 fractions that samples and testing systems are combined in the same exposure system, we propose to separate the metabolism step and toxicity measurement. This allows for a modular combination of metabolic activation by enzymes isolated from rat liver (S9) or a biotechnological alternative (ewoS9^R) with in vitro bioassays that lack metabolic capacity. Benzo(a)pyrene and 2-aminoanthracene were used as model compounds to optimize the conditions for the S9 metabolic degradation/activation step. The Ames assay with Salmonella typhimurium strains TA98 and TA100 was applied to validate the set-up of decoupling the S9 activation/metabolism from the bioassay system. S9 protein concentration of 0.25 mg_protein/mL, a supplement of 0.13 mM NADPH and a pre-incubation time of 100 min are recommended for activation of samples prior to dosing them to in vitro bioassays using the regular dosing protocols of the respective bioassay. EwoS9^R performed equally well as Moltox S9, which is a step forward in developing true animal-free in vitro bioassays. After pre-incubation with S9 fraction, chemicals induced bacteria revertants in both the TA98 and the TA100 assay as efficiently as the standard Ames assay. The pre-incubation of chemicals with S9 fraction could serve for a wide range of cellular in vitro assays to efficiently combine activation and toxicity measurement, which may greatly facilitate the application of these assays for chemical hazard assessment and monitoring of environmental samples.

Utilization of Physiologically Based Pharmacokinetic Modeling in Clinical Pharmacology and Therapeutics: an Overview

The purpose of this review was to assess the advancement of applications for physiologically based pharmacokinetic (PBPK) modeling in various therapeutic areas. We conducted a PubMed search, and 166 articles published between 2012 and 2018 on FDA-approved drug products were selected for further review. Qualifying publications were summarized according to therapeutic area, medication(s) studied, pharmacokinetic model type utilized, simulator program used, and the applications of that modeling. The results showed a 13-fold increase in the number of papers published from 2012 to 2018, with the largest proportion of articles dedicated to the areas of infectious diseases, oncology, and neurology, and application extensions including prediction of drug-drug interactions due to metabolism and/or transporter-mediated effects and understanding drug kinetics in special populations. In addition, we profiled several high-impact studies whose results were used to guide package insert information and formulate dose recommendations. These results show that while utilization of PBPK modeling has drastically increased over the past several years, regulatory support, lack of easy-to-use systems for clinicians, and challenges with model validation remain major challenges for the widespread adoption of this practice in institutional and ambulatory settings. However, PBPK modeling will continue to be a useful tool in the future to assess therapeutic drug monitoring and the growing field of personalized medicine.