Use of Physiologically Based Kinetic Modeling-Based Reverse Dosimetry to Predict in Vivo Toxicity from in Vitro Data

Advancing PFAS risk assessment: Integrative approaches using agent-based modelling and physiologically-based kinetic for environmental and health safety

Per- and polyfluoroalkyl substances (PFAS), ubiquitous in a myriad of consumer and industrial products, and depending on the doses of exposure represent a hazard to both environmental and public health, owing to their persistent, mobile, and bio accumulative properties. These substances exhibit long half-lives in humans and can induce potential immunotoxic effects at low exposure levels, sparking growing concerns. While the European Food Safety Authority (EFSA) has assessed the risk to human health related to the presence of PFAS in food, in which a reduced antibody response to vaccination in infants was considered as the most critical human health effect, a comprehensive grasp of the molecular mechanisms spearheading PFAS-induced immunotoxicity is yet to be attained. Leveraging modern computational tools, including the Agent-Based Model (ABM) Universal Immune System Simulator (UISS) and Physiologically Based Kinetic (PBK) models, a deeper insight into the complex mechanisms of PFAS was sought. The adapted UISS serves as a vital tool in chemical risk assessments, simulating the host immune system's reactions to diverse stimuli and monitoring biological entities within specific adverse health contexts. In tandem, PBK models unravelling PFAS' biokinetics within the body i.e. absorption, distribution, metabolism, and elimination, facilitating the development of time-concentration profiles from birth to 75 years at varied dosage levels, thereby enhancing UISS-TOX's predictive abilities. The integrated use of these computational frameworks shows promises in leveraging new scientific evidence to support risk assessments of PFAS. This innovative approach not only allowed to bridge existing data gaps but also unveiled complex mechanisms and the identification of unanticipated dynamics, potentially guiding more informed risk assessments, regulatory decisions, and associated risk mitigations measures for the future.

Advancing ecotoxicological studies: Utilizing new approach methodologies to enable cross-species extrapolation and reduce avian testing

Ecological risk assessments of agrochemicals have traditionally depended on in vivo guideline tests using northern bobwhite and mallard to provide relevant endpoints for avian species. However, these studies have limitations, including animal welfare concerns, the time and cost involved, limited potential for extrapolation to more realistic exposure conditions, and the lack of mechanistic understanding. The proof-of-concept work presented a case study for thiamethoxam in three avian species, demonstrating the potential of physiologically based kinetic (PBK) modeling to enable dosimetry extrapolations that inform hazard characterization in risk assessment, and reduce the use of avian testing. The model structure for northern bobwhite and mallard contained ten compartments, while an additional ovulation model was included for chicken in the physiological state of egg-laying. The model was first parameterized and evaluated for chicken and northern bobwhite using in vitro kinetic measurements and in vivo toxicokinetic (TK) data. The chicken model was then extrapolated to mallard based on allometric scaling. The models were then used to map the TK profiles across species by simulating internal dose metrics in different avian toxicology studies. These metrics, including peak blood concentrations (Cmax) and area under the curve (AUC) for blood concentration, were determined for acute, subacute, or chronic toxicity endpoints for mallard and northern bobwhite, enabling a quantitative cross-species and cross-route comparison of dosimetry. The results suggested that the chronic toxicological response of birds exposed to thiamethoxam is highly dependent on internal exposure, while mallard appeared to be more dynamically sensitive to thiamethoxam on an acute oral exposure basis. The case study increases the confidence in using new approach methodologies (NAMs) for interpreting avian toxicity studies and facilitating in vitro-in silico-based ecological risk assessments of agrochemicals.

New approach methodologies to confirm developmental toxicity of pharmaceuticals based on weight of evidence

The aim of embryo-fetal developmental toxicity assessments for pharmaceuticals is to inform potential risk of adverse pregnancy outcome, which has traditionally relied on studies in pregnant animals. Recent updates to international safety guidelines (ICH S5R3) have incorporated information on how to use weight of evidence and alternative assays to reduce animal use while still informing risk of fetal harm. Uptake of these alternative approaches has been slow due to limitations in understanding how alternative assays translate to in vivo effects and then relevance to human exposure. To understand the predictivity of new approach methodologies for developmental toxicity (DevTox NAMs), we used two pharmaceutical examples (glasdegib and lorlatinib) to illustrate the value of DevTox NAMs to complement weight of evidence (WoE) assessments while considering the relationship of concentration-effect levels in NAMs to in vivo studies. The in vitro results generated in a battery of assays (mEST, rWEC, zebrafish, and human based stem cells) confirmed the WoE based on literature and further confirmed by preliminary embryo-fetal development data. The data generated for these two compounds supports integrating DevTox NAMs into the developmental toxicity assessment for advanced cancer indications.

Perspectives on Advancing Multimodal Learning in Environmental Science and Engineering Studies

The environment faces increasing anthropogenic impacts, resulting in a rapid increase in environmental issues that undermine the natural capital essential for human wellbeing. These issues are complex and often influenced by various factors represented by data with different modalities. While machine learning (ML) provides data-driven tools for addressing the environmental issues, the current ML models in environmental science and engineering (ES&E) often neglect the utilization of multimodal data. With the advancement in deep learning, multimodal learning (MML) holds promise for comprehensive descriptions of the environmental issues by harnessing data from diverse modalities. This advancement has the potential to significantly elevate the accuracy and robustness of prediction models in ES&E studies, providing enhanced solutions for various environmental modeling tasks. This perspective summarizes MML methodologies and proposes potential applications of MML models in ES&E studies, including environmental quality assessment, prediction of chemical hazards, and optimization of pollution control techniques. Additionally, we discuss the challenges associated with implementing MML in ES&E and propose future research directions in this domain.

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.

Advancing exposure assessment approaches to improve wildlife risk assessment

The exposure assessment component of a Wildlife Ecological Risk Assessment aims to estimate the magnitude, frequency, and duration of exposure to a chemical or environmental contaminant, along with characteristics of the exposed population. This can be challenging in wildlife as there is often high uncertainty and error caused by broad-based, interspecific extrapolation and assumptions often because of a lack of data. Both the US Environmental Protection Agency (USEPA) and European Food Safety Authority (EFSA) have broadly directed exposure assessments to include estimates of the quantity (dose or concentration), frequency, and duration of exposure to a contaminant of interest while considering "all relevant factors." This ambiguity in the inclusion or exclusion of specific factors (e.g., individual and species-specific biology, diet, or proportion time in treated or contaminated area) can significantly influence the overall risk characterization. In this review, we identify four discrete categories of complexity that should be considered in an exposure assessment-chemical, environmental, organismal, and ecological. These may require more data, but a degree of inclusion at all stages of the risk assessment is critical to moving beyond screening-level methods that have a high degree of uncertainty and suffer from conservatism and a lack of realism. We demonstrate that there are many existing and emerging scientific tools and cross-cutting solutions for tackling exposure complexity. To foster greater application of these methods in wildlife exposure assessments, we present a new framework for risk assessors to construct an "exposure matrix." Using three case studies, we illustrate how the matrix can better inform, integrate, and more transparently communicate the important elements of complexity and realism in exposure assessments for wildlife. Modernizing wildlife exposure assessments is long overdue and will require improved collaboration, data sharing, application of standardized exposure scenarios, better communication of assumptions and uncertainty, and postregulatory tracking. Integr Environ Assess Manag 2024;20:674-698. © 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.

Unraveling the micro- and nanoplastic predicament: A human-centric insight

Micro- and nanoplastics are vast anthropogenic pollutants in our direct surroundings with a robust environmental stability and a potential for a long-lasting and increasing global circulation. This has raised concerns among the public and policy makers for human health upon exposure to these particles. The micro- and nanoplastic burden on humans is currently under debate, along with criticism on the experimental approaches used in hazard assessment. The present review presents an overview of the human-relevant aspects associated with the current micro-and nanoplastic burden. We focus on environmental circulation and the estimation of exposure quantities to humans, along with a state-of-the-art overview of particle accumulation in over 15 human organs and other specimen. Additionally, data regarding particle characteristics used in toxicity testing was extracted from 91 studies and discussed considering their environmental and human relevance.

Acetylcholinesterase Inhibition in Rats and Humans Following Acute Fenitrothion Exposure Predicted by Physiologically Based Kinetic Modeling-Facilitated Quantitative In Vitro to In Vivo Extrapolation

Worldwide use of organophosphate pesticides as agricultural chemicals aims to maintain a stable food supply, while their toxicity remains a major public health concern. A common mechanism of acute neurotoxicity following organophosphate pesticide exposure is the inhibition of acetylcholinesterase (AChE). To support Next Generation Risk Assessment for public health upon acute neurotoxicity induced by organophosphate pesticides, physiologically based kinetic (PBK) modeling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE) approach was employed in this study, with fenitrothion (FNT) as an exemplary organophosphate pesticide. Rat and human PBK models were parametrized with data derived from in silico predictions and in vitro incubations. Then, PBK model-based QIVIVE was performed to convert species-specific concentration-dependent AChE inhibition obtained from in vitro blood assays to corresponding in vivo dose-response curves, from which points of departure (PODs) were derived. The obtained values for rats and humans were comparable with reported no-observed-adverse-effect levels (NOAELs). Humans were found to be more susceptible than rats toward erythrocyte AChE inhibition induced by acute FNT exposure due to interspecies differences in toxicokinetics and toxicodynamics. The described approach adequately predicts toxicokinetics and acute toxicity of FNT, providing a proof-of-principle for applying this approach in a 3R-based chemical risk assessment paradigm.

Application of high throughput in vitro metabolomics for hepatotoxicity mode of action characterization and mechanistic-anchored point of departure derivation: a case study with nitrofurantoin

Omics techniques have been increasingly recognized as promising tools for Next Generation Risk Assessment. Targeted metabolomics offer the advantage of providing readily interpretable mechanistic information about perturbed biological pathways. In this study, a high-throughput LC-MS/MS-based broad targeted metabolomics system was applied to study nitrofurantoin metabolic dynamics over time and concentration and to provide a mechanistic-anchored approach for point of departure (PoD) derivation. Upon nitrofurantoin exposure at five concentrations (7.5 µM, 15 µM, 20 µM, 30 µM and 120 µM) and four time points (3, 6, 24 and 48 h), the intracellular metabolome of HepG2 cells was evaluated. In total, 256 uniquely identified metabolites were measured, annotated, and allocated in 13 different metabolite classes. Principal component analysis (PCA) and univariate statistical analysis showed clear metabolome-based time and concentration effects. Mechanistic information evidenced the differential activation of cellular pathways indicative of early adaptive and hepatotoxic response. At low concentrations, effects were seen mainly in the energy and lipid metabolism, in the mid concentration range, the activation of the antioxidant cellular response was evidenced by increased levels of glutathione (GSH) and metabolites from the de novo GSH synthesis pathway. At the highest concentrations, the depletion of GSH, together with alternations reflective of mitochondrial impairments, were indicative of a hepatotoxic response. Finally, a metabolomics-based PoD was derived by multivariate PCA using the whole set of measured metabolites. This approach allows using the entire dataset and derive PoD that can be mechanistically anchored to established key events. Our results show the suitability of high throughput targeted metabolomics to investigate mechanisms of hepatoxicity and derive point of departures that can be linked to existing adverse outcome pathways and contribute to the development of new ones.

The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment

Toxicokinetics plays a crucial role in the health risk assessments of xenobiotics. Classical compartmental models are limited in their ability to determine chemical concentrations in specific organs or tissues, particularly target organs or tissues, and their limited interspecific and exposure route extrapolation hinders satisfactory health risk assessment. In contrast, physiologically based toxicokinetic (PBTK) models quantitatively describe the absorption, distribution, metabolism, and excretion of chemicals across various exposure routes and doses in organisms, establishing correlations with toxic effects. Consequently, PBTK models serve as potent tools for extrapolation and provide a theoretical foundation for health risk assessment and management. This review outlines the construction and application of PBTK models in health risk assessment while analyzing their limitations and future perspectives.

Integrating In Vitro Data and Physiologically Based Kinetic Modeling to Predict and Compare Acute Neurotoxic Doses of Saxitoxin in Rats, Mice, and Humans

Current climate trends are likely to expand the geographic distribution of the toxigenic microalgae and concomitant phycotoxins, making intoxications by such toxins a global phenomenon. Among various phycotoxins, saxitoxin (STX) acts as a neurotoxin that might cause severe neurological symptoms in mammals following consumptions of contaminated seafood. To derive a point of departure (POD) for human health risk assessment upon acute neurotoxicity induced by oral STX exposure, a physiologically based kinetic (PBK) modeling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE) approach was employed. The PBK models for rats, mice, and humans were built using parameters from the literature, in vitro experiments, and in silico predictions. Available in vitro toxicity data for STX were converted to in vivo dose-response curves via the PBK models established for these three species, and POD values were derived from the predicted curves and compared to reported in vivo toxicity data. Interspecies differences in acute STX toxicity between rodents and humans were found, and they appeared to be mainly due to differences in toxicokinetics. The described approach resulted in adequate predictions for acute oral STX exposure, indicating that new approach methodologies, when appropriately integrated, can be used in a 3R-based chemical risk assessment paradigm.

Explore the Dosimetric Relationship between the Intake of Chemical Contaminants and Their Occurrence in Blood and Urine

The dosimetric relationship between the human intake dose of a chemical contaminant (an "external dose") and its concentrations in bodily fluids such as blood and urine (related to an "internal dose"), often characterized by a dose-to-concentration ratio, has critical applications in exposure science, toxicology, and risk assessment, especially in the "new approach methods" era. However, there is a lack of a mechanistic, systematic understanding of how such a dosimetric relationship depends on fundamental chemical properties, such as partition coefficients and biotransformation half-lives. Here, we investigate this issue using a well-evaluated toxicokinetic model, which links external and internal doses by quantifying the absorption and elimination of chemicals. Results are visualized in a series of chemical partitioning space plots, whereby a chemical's dose-to-concentration ratio can be approximately predicted based on its partitioning between air, water, and octanol phases. Our results indicate that when taken in equal doses, chemicals with low volatility and moderate to high hydrophobicity exhibit the highest concentrations in the blood, and chemicals undergoing significant biotransformation tend to exhibit lower concentrations in comparison to their counterparts undergoing negligible biotransformation but possessing similar partitioning properties. Chemicals with high hydrophilicity have the highest concentrations in urine. Such revealed property dependence is similar for both adults and children and for individuals with normal body weights and with obesity. Overall, insights gained from this study are important in predicting blood and urinary concentrations from exposure information and in determining the exposure rate that produces the blood or urinary concentrations observed in biomonitoring studies.

Towards a science-based testing strategy to identify maternal thyroid hormone imbalance and neurodevelopmental effects in the progeny - Part IV: the ECETOC and CLE Proposal for a Thyroid Function-Related Neurodevelopmental Toxicity Testing and Assessment Scheme (Thyroid-NDT-TAS)

Following the European Commission Endocrine Disruptor Criteria, substances shall be considered as having endocrine disrupting properties if they (a) elicit adverse effects, (b) have endocrine activity, and (c) the two are linked by an endocrine mode-of-action (MoA) unless the MoA is not relevant for humans. A comprehensive, structured approach to assess whether substances meet the Endocrine Disruptor Criteria for the thyroid modality (EDC-T) is currently unavailable. Here, the European Centre for Ecotoxicology and Toxicology of Chemicals Thyroxine Task Force and CropLife Europe propose a Thyroid Function-Related Neurodevelopmental Toxicity Testing and Assessment Scheme (Thyroid-NDT-TAS). In Tier 0, before entering the Thyroid-NDT-TAS, all available in vivo, in vitro and in silico data are submitted to weight-of-evidence (WoE) evaluations to determine whether the substance of interest poses a concern for thyroid disruption. If so, Tier 1 of the Thyroid-NDT-TAS includes an initial MoA and human relevance assessment (structured by the key events of possibly relevant adverse outcome pathways) and the generation of supportive in vitro/in silico data, if relevant. Only if Tier 1 is inconclusive, Tier 2 involves higher-tier testing to generate further thyroid- and/or neurodevelopment-related data. Tier 3 includes the final MoA and human relevance assessment and an overarching WoE evaluation to draw a conclusion on whether, or not, the substance meets the EDC-T. The Thyroid-NDT-TAS is based on the state-of-the-science, and it has been developed to minimise animal testing. To make human safety assessments more accurate, it is recommended to apply the Thyroid-NDT-TAS during future regulatory assessments.

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.

New approach methodologies: A quantitative in vitro to in vivo extrapolation case study with PFASs

PER: and polyfluoroalkyl substances (PFASs) have been associated with increased blood lipids in humans. Perfluorooctanoic acid (PFOA) has been also linked with elevated alanine transferase (ALT) serum levels in humans, and in rodents the liver is a main target organ for many PFASs. With the focus on New Approach Methodologies, the chronic oral equivalent effect doses were calculated for PFOA, PFNA (perfluorononanoic acid), PFHxS (perfluorohexanesulfonic acid) and PFOS (perfluorooctane sulfonic acid) based on in vitro effects measured in the HepaRG cell line. Selected in vitro readouts were considered biomarkers for lipid disturbances and hepatotoxicity. Concentration-response data obtained from HepaRG cells on triglyceride (TG) accumulation and expression changes of 12 selected genes (some involved in cholesterol homeostasis) were converted into corresponding human dose-response data, using physiologically based kinetic (PBK) model-facilitated reverse dosimetry. Next to this, the biokinetics of the chemicals were studied in the cell system. The current European dietary PFASs exposure overlaps with the calculated oral equivalent effect doses, indicating that the latter may lead to interference with hepatic gene expression and lipid metabolism. These findings illustrate an in vitro-in silico methodology, which can be applied for more PFASs, to select those that should be prioritized for further hazard characterization.

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.

Cytotoxicity and Thermal Characterization Assessment of Excipients for the Development of Innovative Lyophilized Formulations for Oncological Applications

Freeze-drying, also known as lyophilization, significantly improves the storage, stability, shelf life, and clinical translation of biopharmaceuticals. On the downside, this process faces complex challenges, i.e., the presence of freezing and drying stresses for the active compounds, the uniformity and consistency of the final products, and the efficiency and safety of the reconstituted lyophilized formulations. All these requirements can be addressed by adding specific excipients that can protect and stabilize the active ingredient during lyophilization, assisting in the formation of solid structures without interfering with the biological and/or pharmaceutical action of the reconstituted products. However, these excipients, generally considered safe and inert, could play an active role in the formulation interacting with the biological cellular machinery and promoting toxicity. Any side effects should be carefully identified and characterized to better tune any treatments in terms of concentrations and administration times. In this work, various concentrations in the range of 1 to 100 mg/mL of cellobiose, lactose, sucrose, trehalose, isoleucine, glycine, methionine, dextran, mannitol, and (2-hydroxypropyl)-β-cyclodextrin were evaluated in terms of their ability to create uniform and solid lyophilized structures. The freeze-dried products were then reconstituted in the appropriate cell culture media to assess their in vitro cytotoxicity on both a healthy cell line (B-lymphocytes) and their tumoral lymphoid counterpart (Daudi). Results showed that at 10 mg/mL, all the excipients demonstrated suitable lyophilized solid structures and high tolerability by both cell lines, while dextran was the only excipient well-tolerated also up to 100 mg/mL. An interesting result was shown for methionine, which even at 10 mg/mL, selectively affected the viability of the cancerous cell line only, opening future perspectives for antitumoral applications.

Use of Physiologically Based Kinetic Modeling-Based Reverse Dosimetry to Predict In Vivo Nrf2 Activation by EGCG and Its Colonic Metabolites in Humans

(-)-Epigallocatechin gallate (EGCG) is prone to microbial metabolism when reaching the colon. This study aimed to develop a human physiologically based kinetic (PBK) model for EGCG, with sub-models for its colonic metabolites gallic acid and pyrogallol. Results show that the developed PBK model could adequately predict in vivo time-dependent blood concentrations of EGCG after either the single or repeated oral administration of EGCG under both fasting and non-fasting conditions. The predicted in vivo blood C_max of EGCG indicates that the Nrf2 activation is limited, while concentrations of its metabolites in the intestinal tract may reach levels that are higher than that of EGCG and also high enough to activate Nrf2 gene transcription. Taken together, combining in vitro data with a human PBK model allowed the prediction of a dose-response curve for EGCG-induced Nrf2-mediated gene expression in humans and provided insights into the contribution of gut microbial metabolites to this effect.

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.

Population pharmacokinetic model to generate mechanistic insights in bile acid homeostasis and drug-induced cholestasis

Bile acids (BA) fulfill a wide range of physiological functions, but are also involved in pathologies, such as cholestasis. Cholestasis is characterized by an intrahepatic accumulation of BAs and subsequent spillage to the systemic circulation. The aim of the present study was to develop physiologically based kinetic (PBK) models that would provide a tool to predict dose-dependent BA accumulation in humans upon treatment with a Bile Salt Export Pump (BSEP) inhibitor. We developed a PBK model describing the BA homeostasis using glycochenodeoxycholic acid as an exemplary BA. Population wide distributions of BSEP abundances were incorporated in the PBK model using Markov Chain Monte Carlo simulations, and alternatively the total amount of BAs was scaled empirically to describe interindividual differences in plasma BA levels. Next, the effects of the BSEP inhibitor bosentan on the BA levels were simulated. The PBK model developed adequately predicted the in vivo BA dynamics. Both the Markov Chain Monte Carlo simulations based on a distribution of BSEP abundances and empirical scaling of the total BA pool readily described the variations within and between data in human volunteers. Bosentan treatment disproportionally increased the maximum BA concentration in individuals with a large total BA pool or low BSEP abundance. Especially individuals having a large total BA pool size and a low BSEP abundance were predicted to be at risk for rapid saturation of BSEP and subsequent intrahepatic BA accumulation. This model provides a first estimate of personalized safe therapeutic external dose levels of compounds with BSEP-inhibitory properties.

Pluripotent stem cell assays: Modalities and applications for predictive developmental toxicity

This manuscript provides a review focused on embryonic stem cell-based models and their place within the landscape of alternative developmental toxicity assays. Against the background of the principles of developmental toxicology, the wide diversity of alternative methods using pluripotent stem cells developed in this area over the past half century is reviewed. In order to provide an overview of available models, a systematic scoping review was conducted following a published protocol with inclusion criteria, which were applied to select the assays. Critical aspects including biological domain, readout endpoint, availability of standardized protocols, chemical domain, reproducibility and predictive power of each assay are described in detail, in order to review the applicability and limitations of the platform in general and progress moving forward to implementation. The horizon of innovative routes of promoting regulatory implementation of alternative methods is scanned, and recommendations for further work are given.

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.

Cytochrome P450 isoforms contribution, plasma protein binding, toxicokinetics of enniatin A in rats and in vivo clearance prediction in humans

Emerging mycotoxins, such as enniatin A (ENNA), are becoming a worldwide concern owing to their presence in different types of food and feed. However, comprehensive toxicokinetic data that links intake, exposure and toxicological effects of ENNA has not been elucidated yet. Therefore, the present study investigated the in vitro (rat and human) and in vivo (rat) toxicokinetic properties of ENNA. Towards this, an easily applicable and sensitive bioanalytical method was developed and validated for the estimation of ENNA in rat plasma. ENNA exhibited high plasma protein binding (99%), high hepatic clearance and mainly underwent metabolism via CYP3A4 (74%). The in-house predicted hepatic clearance (54 mL/min/kg) and observed in vivo rat clearance (55 mL/min/kg) were comparable. The predicted in vivo human hepatic clearance was 18 mL/min/kg. ENNA underwent slow absorption (T_max = 4 h) and rapid elimination following oral administration to rats. The absolute oral bioavailability was 47%. The toxicokinetic findings for ENNA from this study will help in designing and interpreting toxicological studies in rats. Besides, these findings could be used in physiologically based toxicokinetic (PBTK) model development for exposure predictions and risk assessment for ENNA in humans.

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).

The Role of Kinetics as Key Determinant in Toxicity of Pyrrolizidine Alkaloids and Their N-Oxides

Pyrrolizidine alkaloids (PAs) are a large group of plant constituents of which especially the 1,2- unsaturated PAs raise a concern because of their liver toxicity and potential genotoxic carcinogenicity. This toxicity of PAs depends on their kinetics. Differences in absorption, distribution, metabolism, and excretion (ADME) characteristics of PAs may substantially alter the relative toxicity of PAs. As a result, kinetics will also affect relative potency (REP) values. The present review summarizes the current state-of-the art on PA kinetics and resulting consequences for toxicity and illustrates how physiologically-based kinetic (PBK) modelling can be applied to take kinetics into account when defining the relative differences in toxicity between PAs in the in vivo situation. We conclude that toxicokinetics play an important role in the overall toxicity of pyrrolizidine alkaloids. and that kinetics should therefore be considered when defining REP values for combined risk assessment. New approach methodologies (NAMs) can be of use to quantify these kinetic differences between PAs and their N-oxides, thus contributing to the 3Rs (Replacement, Reduction and Refinement) in animal studies.

Novel Insights into Pyrrolizidine Alkaloid Toxicity and Implications for Risk Assessment: Occurrence, Genotoxicity, Toxicokinetics, Risk Assessment-A Workshop Report

This paper reports on the major contributions and results of the 2nd International Workshop of Pyrrolizidine Alkaloids held in September 2020 in Kaiserslautern, Germany. Pyrrolizidine alkaloids are among the most relevant plant toxins contaminating food, feed, and medicinal products of plant origin. Hundreds of PA congeners with widespread occurrence are known, and thousands of plants are assumed to contain PAs. Due to certain PAs' pronounced liver toxicity and carcinogenicity, their occurrence in food, feed, and phytomedicines has raised serious human health concerns. This is particularly true for herbal teas, certain food supplements, honey, and certain phytomedicinal drugs. Due to the limited availability of animal data, broader use of in vitro data appears warranted to improve the risk assessment of a large number of relevant, 1,2-unsaturated PAs. This is true, for example, for the derivation of both toxicokinetic and toxicodynamic data. These efforts aim to understand better the modes of action, uptake, metabolism, elimination, toxicity, and genotoxicity of PAs to enable a detailed dose-response analysis and ultimately quantify differing toxic potencies between relevant PAs. Accordingly, risk-limiting measures comprising production, marketing, and regulation of food, feed, and medicinal products are discussed.

Integrating in vitro chemical transplacental passage into a generic PBK model: A QIVIVE approach

With the increasing application of cell culture models as primary tools for predicting chemical safety, the quantitative extrapolation of the effective dose from in vitro to in vivo (QIVIVE) is of increasing importance. For developmental toxicity this requires scaling the in vitro observed dose-response characteristics to in vivo fetal exposure, while integrating maternal in vivo kinetics during pregnancy, in particular transplacental transfer. Here the transfer of substances across the placental barrier, has been studied using the in vitro BeWo cell assay and six embryotoxic compounds of different kinetic complexity. The BeWo assay results were incorporated in an existing generic Physiologically Based Kinetic (PBK) model which for this purpose was extended with rat pregnancy. Finally, as a "proof of principle", the BeWo PBK model was used to perform a QIVIVE based on developmental toxicity as observed in various different in vitro toxicity assays. The BeWo results illustrated different transport profiles of the chemicals across the BeWo monolayer, allocating the substances into two distinct groups: the 'quickly-transported' and the 'slowly-transported'. BeWo PBK exposure simulations during gestation were compared to experimentally measured maternal blood and fetal concentrations and a reverse dosimetry approach was applied to translate in vitro observed embryotoxicity into equivalent in vivo dose-response curves. This approach allowed for a direct comparison of the in vitro dose-response characteristics as observed in the Whole Embryo Culture (WEC), and the Embryonic Stem Cell test (cardiac:ESTc and neural:ESTn) with in vivo rat developmental toxicity data. Overall, the in vitro to in vivo comparisons suggest a promising future for the application of such QIVIVE methodologies for screening and prioritization purposes of developmental toxicants. Nevertheless, the clear need for further improvements is acknowledged for a wider application of the approach in chemical safety assessment.

EDC-induced mechanisms of immunotoxicity: a systematic review

Endocrine-disrupting chemicals (EDCs) refer to a group of chemicals that cause adverse effects in human health, impairing hormone production and regulation, resulting in alteration of homeostasis, reproductive, and developmental, and immune system impairments. The immunotoxicity of EDCs involves many mechanisms altering gene expression that depend on the activation of nuclear receptors such as the aryl hydrocarbon receptor (AHR), the estrogen receptor (ER), and the peroxisome proliferator-activated receptor (PPAR), which also results in skin and intestinal disorders, microbiota alterations and inflammatory diseases. This systematic review aims to review different mechanisms of immunotoxicity and immunomodulation of T cells, focusing on T regulatory (Treg) and Th17 subsets, B cells, and dendritic cells (DCs) caused by specific EDCs such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), bisphenols (BPs) and polyfluoroalkyl substances (PFASs). To achieve this objective, a systematic study was conducted searching various databases including PubMed and Scopus to find in-vitro, in-vivo, and biomonitoring studies that examine EDC-dependent mechanisms of immunotoxicity. While doing the systematic review, we found species- and cell-specific outcomes and a translational gap between in-vitro and in-vivo experiments. Finally, an adverse outcome pathway (AOP) framework is proposed, which explains mechanistically toxicity endpoints emerging from different EDCs having similar key events and can help to improve our understanding of EDCs mechanisms of immunotoxicity. In conclusion, this review provides insights into the mechanisms of immunotoxicity mediated by EDCs and will help to improve human health risk assessment.

Integrating toxicokinetics into toxicology studies and the human health risk assessment process for chemicals: Reduced uncertainty, better health protection

The database of practical examples where toxicokinetic (TK) data has benefitted all stages of the human health risk assessment process are increasingly being published and accepted. This review aimed to highlight and summarise notable examples and to describe the "state of the art" in this field. The overall recommendation is that for any in vivo animal study conducted, measurements of TK should be very carefully considered for inclusion as the numerous benefits this brings continues to grow, particularly during the current march towards animal free toxicology testing and ambitions to eventually conduct human health risk assessments entirely based upon non-animal methods.

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.

Derivation of a Human In Vivo Benchmark Dose for Perfluorooctanoic Acid 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) modeling, global sensitivity analysis (GSA), approximate Bayesian computation (ABC), and Markov Chain Monte Carlo (MCMC) simulation was developed to facilitate quantitative in vitro to in vivo extrapolation (QIVIVE). The workflow accounts for parameter and model uncertainty within a computationally efficient framework. The workflow was tested using a human PBK model for perfluorooctanoic acid (PFOA) and high throughput screening (HTS) in vitro concentration–response data, determined in a human liver cell line, from the ToxCast/Tox21 database. In vivo benchmark doses (BMDs) for PFOA intake (ng/kg BW/day) and drinking water exposure concentrations (µg/L) were calculated from the in vivo dose responses and compared to intake values derived by the European Food Safety Authority (EFSA). The intake benchmark dose lower confidence limit (BMDL₅) of 0.82 was similar to 0.86 ng/kg BW/day for altered serum cholesterol levels derived by EFSA, whereas the intake BMDL₅ of 6.88 was six-fold higher than the value of 1.14 ng/kg BW/day for altered antibody titer also derived by the EFSA. Application of a chemical-specific adjustment factor (CSAF) of 1.4, allowing for inter-individual variability in kinetics, based on biological half-life, gave an intake BMDL₅ of 0.59 for serum cholesterol and 4.91 (ng/kg BW/day), for decreased antibody titer, which were 0.69 and 4.31 the EFSA-derived values, respectively. The corresponding BMDL₅ for drinking water concentrations, for estrogen receptor binding activation associated with breast cancer, pregnane X receptor binding associated with altered serum cholesterol levels, thyroid hormone receptor α binding leading to thyroid disease, and decreased antibody titer (pro-inflammation from cytokines) were 0.883, 0.139, 0.086, and 0.295 ng/ml, respectively, with application of no uncertainty factors. These concentrations are 5.7-, 36-, 58.5-, and 16.9-fold lower than the median measured drinking water level for the general US population which is approximately, 5 ng/ml.

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.

Physiologically based kinetic modelling based prediction of in vivo rat and human acetylcholinesterase (AChE) inhibition upon exposure to diazinon

The present study predicts in vivo human and rat red blood cell (RBC) acetylcholinesterase (AChE) inhibition upon diazinon (DZN) exposure using physiological based kinetic (PBK) modelling-facilitated reverse dosimetry. Due to the fact that both DZN and its oxon metabolite diazoxon (DZO) can inhibit AChE, a toxic equivalency factor (TEF) was included in the PBK model to combine the effect of DZN and DZO when predicting in vivo AChE inhibition. The PBK models were defined based on kinetic constants derived from in vitro incubations with liver fractions or plasma of rat and human, and were used to translate in vitro concentration-response curves for AChE inhibition obtained in the current study to predicted in vivo dose-response curves. The predicted dose-response curves for rat matched available in vivo data on AChE inhibition, and the benchmark dose lower confidence limits for 10% inhibition (BMDL10 values) were in line with the reported BMDL10 values. Humans were predicted to be 6-fold more sensitive than rats in terms of AChE inhibition, mainly because of inter-species differences in toxicokinetics. It is concluded that the TEF-coded DZN PBK model combined with quantitative in vitro to in vivo extrapolation (QIVIVE) provides an adequate approach to predict RBC AChE inhibition upon acute oral DZN exposure, and can provide an alternative testing strategy for derivation of a point of departure (POD) in risk assessment.

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.

Prediction of dose-dependent in vivo acetylcholinesterase inhibition by profenofos in rats and humans using physiologically based kinetic (PBK) modeling-facilitated reverse dosimetry

Organophosphate pesticides (OPs) are known to inhibit acetylcholine esterase (AChE), a critical effect used to establish health-based guidance values. This study developed a combined in vitro–in silico approach to predict AChE inhibition by the OP profenofos in rats and humans. A physiologically based kinetic (PBK) model was developed for both species. Parameter values for profenofos conversion to 4-bromo-2-chlorophenol (BCP) were derived from in vitro incubations with liver microsomes, liver cytosol, and plasma from rats (catalytic efficiencies of 1.1, 2.8, and 0.19 ml/min/mg protein, respectively) and humans (catalytic efficiencies of 0.17, 0.79, and 0.063 ml/min/mg protein, respectively), whereas other chemical-related parameter values were derived using in silico calculations. The rat PBK model was evaluated against literature data on urinary excretion of conjugated BCP. Concentration-dependent inhibition of rat and human AChE was determined in vitro and these data were translated with the PBK models to predicted dose-dependent AChE inhibition in rats and humans in vivo. Comparing predicted dose-dependent AChE inhibition in rats to literature data on profenofos-induced AChE inhibition revealed an accurate prediction of in vivo effect levels. Comparison of rat predictions (BMDL10 of predicted dose–response data of 0.45 mg/kg bw) and human predictions (BMDL10 of predicted dose–response data of 0.01 mg/kg bw) suggests that humans are more sensitive than rats, being mainly due to differences in kinetics. Altogether, the results demonstrate that in vivo AChE inhibition upon acute exposure to profenofos was closely predicted in rats, indicating the potential of this novel approach method in chemical hazard assessment.

A versatile, compartmentalised gut-on-a-chip system for pharmacological and toxicological analyses

A novel, integrated, in vitro gastrointestinal (GI) system is presented to study oral bioavailability parameters of small molecules. Three compartments were combined into one hyphenated, flow-through set-up. In the first compartment, a compound was exposed dynamically to enzymatic digestion in three consecutive microreactors, mimicking the processes of the mouth, stomach, and intestine. The resulting solution (chyme) continued to the second compartment, a flow-through barrier model of the intestinal epithelium allowing absorption of the compound and metabolites thereof. The composition of the effluents from the barrier model were analysed either offline by electrospray-ionisation-mass spectrometry (ESI-MS), or online in the final compartment using chip-based ESI-MS. Two model drugs, omeprazole and verapamil, were used to test the integrated model. Omeprazole was shown to be broken down upon treatment with gastric acid, but reached the cell barrier unharmed when introduced to the system in a manner emulating an enteric-coated formulation. In contrast, verapamil was unaffected by digestion. Finally, a reduced uptake of verapamil was observed when verapamil was introduced to the system dissolved in apple juice, a simple food matrix. It is envisaged that this integrated, compartmentalised GI system has potential for enabling future research in the fields of pharmacology, toxicology, and nutrition.

Species Differences in in vitro and Estimated in vivo Kinetics for Intestinal Microbiota Mediated Metabolism of Acetyl-deoxynivalenols

SCOPE

Deoxynivalenol (DON) and its acetylated derivatives 3-acetyl-DON (3-Ac-DON) and 15-acetyl-DON (15-Ac-DON) are important mycotoxins of concern in the modern food chain.

METHODS AND RESULTS

The present study reveals that the rate of de-acetylation in in vitro anaerobic fecal incubations decreased in the order rat > mouse > human > pig for 3-Ac-DON, and mouse > human > rat > pig for 15-Ac-DON. The ratio between the de-acetylation rate of 3-Ac-DON and 15-Ac-DON varies with the species. Scaling of the kinetic parameters to the in vivo situation results in catalytic efficiencies decreasing in the order human > rat > pig > mouse for 3-Ac-DON and human > pig > rat > mouse for 15-Ac-DON. The results obtained indicate that in mice, 3-Ac-DON can be fully deconjugated while 15-Ac-DON cannot. In rats, pigs, and humans, both 3-Ac-DON and 15-Ac-DON can be totally transformed by gut fecal microbiota during the estimated intestinal residence time. A correlation analysis between the deacetylation rate and the relative abundance of the microbiome suggests Lachnospiraceae may be involved in the deacetylation process.

CONCLUSION

It is concluded that interspecies differences in deacetylation of acetylated DONs exist but that in risk assessment assumption of complete intestinal deconjugation provides an adequate approach.

Incorporating renal excretion via the OCT2 transporter in physiologically based kinetic modelling to predict in vivo kinetics of mepiquat in rat

The present study aimed at incorporating active renal excretion via the organic cation transporter 2 (OCT2) into a generic rat physiologically based kinetic (PBK) model using an in vitro human renal proximal tubular epithelial cell line (SA7K) and mepiquat chloride (MQ) as the model compound. The Vmax (10.5 pmol/min/mg protein) and Km (20.6 μM) of OCT2 transport of MQ were determined by concentration-dependent uptake in SA7K cells using doxepin as inhibitor. PBK model predictions incorporating these values in the PBK model were 6.7-8.4-fold different from the reported in vivo data on the blood concentration of MQ in rat. Applying an overall scaling factor that also corrects for potential differences in OCT2 activity in the SA7K cells and in vivo kidney cortex and species differences resulted in adequate predictions for in vivo kinetics of MQ in rat (2.3-3.2-fold). The results indicate that using SA7K cells to define PBK parameters for active renal OCT2 mediated excretion with adequate scaling enables incorporation of renal excretion via the OCT2 transporter in PBK modelling to predict in vivo kinetics of mepiquat in rat. This study demonstrates a proof-of-principle on how to include active renal excretion into generic PBK models.

Prediction of the dose range for adverse neurological effects of amiodarone in patients from an in vitro toxicity test by in vitro-in vivo extrapolation

Amiodarone is an antiarrhythmic agent inducing adverse effects on the nervous system, among others. We applied physiologically based pharmacokinetic (PBPK) modeling combined with benchmark dose modeling to predict, based on published in vitro data, the in vivo dose of amiodarone which may lead to adverse neurological effects in patients. We performed in vitro–in vivo extrapolation (IVIVE) from concentrations measured in the cell lysate of a rat brain 3D cell model using a validated human PBPK model. Among the observed in vitro effects, inhibition of choline acetyl transferase (ChAT) was selected as a marker for neurotoxicity. By reverse dosimetry, we transformed the in vitro concentration–effect relationship into in vivo effective human doses, using the calculated in vitro area under the curve (AUC) of amiodarone as the pharmacokinetic metric. The upper benchmark dose (BMDU) was calculated and compared with clinical doses eliciting neurological adverse effects in patients. The AUCs in the in vitro brain cell culture after 14-day repeated dosing of nominal concentration equal to 1.25 and 2.5 µM amiodarone were 1.00 and 1.99 µg*h/mL, respectively. The BMDU was 385.4 mg for intravenous converted to 593 mg for oral application using the bioavailability factor of 0.65 as reported in the literature. The predicted dose compares well with neurotoxic doses in patients supporting the hypothesis that impaired ChAT activity may be related to the molecular/cellular mechanisms of amiodarone neurotoxicity. Our study shows that predicting effects from in vitro data together with IVIVE can be used at the initial stage for the evaluation of potential adverse drug reactions and safety assessment in humans.

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.

Novel testing strategy for prediction of rat biliary excretion of intravenously administered estradiol-17β glucuronide

The aim of the present study was to develop a generic rat physiologically based kinetic (PBK) model that includes a novel testing strategy where active biliary excretion is incorporated using estradiol-17β glucuronide (E217βG) as the model substance. A major challenge was the definition of the scaling factor for the in vitro to in vivo conversion of the PBK-model parameter Vmax. In vitro values for the Vmax and Km for transport of E217βG were found in the literature in four different studies based on experiments with primary rat hepatocytes. The required scaling factor was defined based on fitting the PBK model-based predicted values to reported experimental data on E217βG blood levels and cumulative biliary E217βG excretion. This resulted in a scaling factor of 129 mg protein/g liver. With this scaling factor the PBK model predicted the in vivo data for blood and cumulative biliary E217βG levels with on average of less than 1.8-fold deviation. The study provides a proof of principle on how biliary excretion can be included in a generic PBK model using primary hepatocytes to define the kinetic parameters that describe the biliary excretion.

An in vitromodel to quantify interspecies differences in kinetics for intestinal microbial bioactivation and detoxification of zearalenone

Zearalenone (ZEN) is a mycotoxin known for its estrogenic activities. The metabolism of ZEN plays a role in the interspecies differences in sensitivity to ZEN, and is known to occur in the liver and via the intestinal microbiota, although the relative contribution of these two pathways remains to be characterized. In the present study a fecal in vitro model was optimized and used to quantify the interspecies differences in kinetics of the intestinal microbial metabolism of ZEN in rat, pig and human. V_max, K_m, and catalytic efficiencies (k_cat) were determined, and results obtained reveal that the k_cat values for formation of α-ZEL and β-ZEL amounted to 0.73 and 0.12 mL/h/kg bw for human microbiota, 2.6 and 1.3 mL/h/kg bw for rat microbiota and 9.4 and 6.3 mL/h/kg bw for pig microbiota showing that overall ZEN metabolism increased in the order human < rat < pig microbiota. Expressed per kg bw the k_cat for ZEN metabolism by the liver surpassed that of the intestinal microbiota in all three species. In conclusion, it is estimated that the activity of the intestinal colon microbiome may be up to 36 % of the activity of the liver, and that it can additionally contribute to the species differences in bioactivation and detoxification and thus the toxicity of ZEN in pigs and rats but not in humans. The results highlight the importance of the development of human specific models for the assessment of the metabolism of ZEN.

Predicting the Acute Liver Toxicity of Aflatoxin B1 in Rats and Humans by an In Vitro-In Silico Testing Strategy

SCOPE

High-level exposure to aflatoxin B1 (AFB1) is known to cause acute liver damage and fatality in animals and humans. The intakes actually causing this acute toxicity have so far been estimated based on AFB1 levels in contaminated foods or biomarkers in serum. The aim of the present study is to predict the doses causing acute liver toxicity of AFB1 in rats and humans by an in vitro-in silico testing strategy.

METHODS AND RESULTS

Physiologically based kinetic (PBK) models for AFB1 in rats and humans are developed. The models are used to translate in vitro concentration-response curves for cytotoxicity in primary rat and human hepatocytes to in vivo dose-response curves using reverse dosimetry. From these data, the dose levels at which toxicity would be expected are obtained and compared to toxic dose levels from available rat and human case studies on AFB1 toxicity. The results show that the in vitro-in silico testing strategy can predict dose levels causing acute toxicity of AFB1 in rats and human.

CONCLUSIONS

Quantitative in vitro in vivo extrapolation (QIVIVE) using PBK modeling-based reverse dosimetry can predict AFB1 doses that cause acute liver toxicity in rats and human.

Integrating in vitro data and physiologically based kinetic modeling-facilitated reverse dosimetry to predict human cardiotoxicity of methadone

Development of novel testing strategies to detect adverse human health effects is of interest to replace in vivo-based drug and chemical safety testing. The aim of the present study was to investigate whether physiologically based kinetic (PBK) modeling-facilitated conversion of in vitro toxicity data is an adequate approach to predict in vivo cardiotoxicity in humans. To enable evaluation of predictions made, methadone was selected as the model compound, being a compound for which data on both kinetics and cardiotoxicity in humans are available. A PBK model for methadone in humans was developed and evaluated against available kinetic data presenting an adequate match. Use of the developed PBK model to convert concentration–response curves for the effect of methadone on human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) in the so-called multi electrode array (MEA) assay resulted in predictions for in vivo dose–response curves for methadone-induced cardiotoxicity that matched the available in vivo data. The results also revealed differences in protein plasma binding of methadone to be a potential factor underlying variation between individuals with respect to sensitivity towards the cardiotoxic effects of methadone. The present study provides a proof-of-principle of using PBK modeling-based reverse dosimetry of in vitro data for the prediction of cardiotoxicity in humans, providing a novel testing strategy in cardiac safety studies.

An integrated pathway based on in vitro data for the human hazard assessment of nanomaterials

In line with the 3R concept, nanotoxicology is shifting from a phenomenological to a mechanistic approach based on in vitro and in silico methods, with a consequent reduction in animal testing. Risk Assessment (RA) and Life Cycle Assessment (LCA) methodologies, which traditionally rely on in vivo toxicity studies, will not be able to keep up with the pace of development of new nanomaterials unless they adapt to use this new type of data. While tools and models are already available and show a great potential for future use in RA and LCA, currently none is able alone to quantitatively assess human hazards (i.e. calculate chronic NOAEL or ED₅₀ values). By highlighting which models and approaches can be used in a quantitative way with the available knowledge and data, we propose an integrated pathway for the use of in vitro data in RA and LCA. Starting with the characterization of nanoparticles' properties, the pathway then investigates how to select relevant in vitro human data, and how to bridge in vitro dose-response relationships to in vivo effects. If verified, this approach would allow RA and LCA to stir up the development of nanotoxicology by giving indications about the data and quality requirements needed in risk methodologies.

Use of Physiologically Based Kinetic Modeling to Predict Rat Gut Microbial Metabolism of the Isoflavone Daidzein to S-Equol and Its Consequences for ERα Activation

SCOPE

To predict gut microbial metabolism of xenobiotics and the resulting plasma concentrations of metabolites formed, an in vitro-in silico-based testing strategy is developed using the isoflavone daidzein and its gut microbial metabolite S-equol as model compounds.

METHODS AND RESULTS

Anaerobic rat fecal incubations are optimized and performed to derive the apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for gut microbial conversion of daidzein to dihydrodaidzein, S-equol, and O-desmethylangolensin, which are input as parameters for a physiologically based kinetic (PBK) model. The inclusion of gut microbiota in the PBK model allows prediction of S-equol concentrations and slightly reduced predicted maximal daidzein concentrations from 2.19 to 2.16 µm. The resulting predicted concentrations of daidzein and S-equol are comparable to in vivo concentrations reported.

CONCLUSION

The optimized in vitro approach to quantify kinetics for gut microbial conversions, and the newly developed PBK model for rats that includes gut microbial metabolism, provide a unique tool to predict the in vivo consequences of daidzein microbial metabolism for systemic exposure of the host to daidzein and its metabolite S-equol. The predictions reveal a dominant role for daidzein in ERα-mediated estrogenicity despite the higher estrogenic potency of its microbial metabolite S-equol.