The short-term toxicity and metabolome of dicyclopentadiene

Dicyclopentadiene (DCPD) was investigated in a 14-day oral rat toxicity study based on the OECD 407 guideline in combination with plasma metabolomics. Wistar rats received the compound daily via gavage at dose levels of 0, 50 and 150 mg/kg bw. The high dose induced transient clinical signs of toxicity and in males only reduced body weight gain. High dose liver changes were characterized by altered clinical chemistry parameters in both sexes and pathological changes in females. In high dose males an accumulation of alpha-2 u-globulin in the kidney was noted. Comparing the DCPD metabolome with previously established specific metabolome patterns in the MetaMap® Tox data base suggested that the high dose would result in liver enzyme induction leading to increased breakdown of thyroid hormones for males and females. An indication for liver toxicity in males was also noted. Metabolomics also suggested an effect on the functionality of the adrenals in high dose males, which together with published data, is suggestive of a stress related effect in this organ. The results of the present 14-day combined toxicity and metabolome investigations were qualitatively in line with literature data from subchronic oral studies in rats with DCPD. Importantly no other types of organ toxicity, or hormone dysregulation beyond the ones associated with liver enzyme induction and stress were indicated, again in line with results of published 90-day studies. It is therefore suggested that short term "smart" studies, combining classical toxicity with 'omics technologies, could be a 2 R (refine and reduce) new approach method allowing for the reduction of in vivo toxicity testing.

Analysis of health concerns not addressed by REACH for low tonnage chemicals and opportunities for new approach methodology

In Registration, Evaluation, Authorisation and Restriction of Chemicals (REACH) the criterion for deciding the studies that must be performed is the annual tonnage of the chemical manufactured or imported into the EU. The annual tonnage may be considered as a surrogate for levels of human exposure but this does not take into account the physico-chemical properties and use patterns that determine exposure. Chemicals are classified using data from REACH under areas of health concern covering effects on the skin and eye; sensitisation; acute, repeated and prolonged systemic exposure; effects on genetic material; carcinogenicity; and reproduction and development. We analysed the mandated study lists under REACH for each annual tonnage band in terms of the information they provide on each of the areas of health concern. Using the European Chemicals Agency (ECHA) REACH Registration data base of over 20,000 registered substances, we found that only 19% of registered substances have datasets on all areas of health concern. Information limited to acute exposure, sensitisation and genotoxicity was found for 62%. The analysis highlighted the shortfall of information mandated for substances in the lower tonnage bands. Deploying New Approach Methodologies (NAMs) at this lower tonnage band to assess health concerns which are currently not covered by REACH, such as repeat and extended exposure and carcinogenicity, would provide additional information and would be a way for registrants and regulators to gain experience in the use of NAMs. There are currently projects in Europe aiming to develop NAM-based assessment frameworks and they could find their first use in assessing low tonnage chemicals once confidence has been gained by their evaluation with data rich chemicals.

A high-throughput metabolomics in vitro platform for the characterization of hepatotoxicity

Cell-based metabolomics provides multiparametric physiologically relevant readouts that can be highly advantageous for improved, biologically based decision making in early stages of compound development. Here, we present the development of a 96-well plate LC-MS/MS-based targeted metabolomics screening platform for the classification of liver toxicity modes of action (MoAs) in HepG2 cells. Different parameters of the workflow (cell seeding density, passage number, cytotoxicity testing, sample preparation, metabolite extraction, analytical method, and data processing) were optimized and standardized to increase the efficiency of the testing platform. The applicability of the system was tested with seven substances known to be representative of three different liver toxicity MoAs (peroxisome proliferation, liver enzyme induction, and liver enzyme inhibition). Five concentrations per substance, aimed at covering the complete dose-response curve, were analyzed and 221 uniquely identified metabolites were measured, annotated, and allocated in 12 different metabolite classes such as amino acids, carbohydrates, energy metabolism, nucleobases, vitamins and cofactors, and diverse lipid classes. Multivariate and univariate analyses showed a dose response of the metabolic effects, a clear differentiation between liver toxicity MoAs and resulted in the identification of metabolite patterns specific for each MoA. Key metabolites indicative of both general and mechanistic specific hepatotoxicity were identified. The method presented here offers a multiparametric, mechanistic-based, and cost-effective hepatotoxicity screening that provides MoA classification and sheds light into the pathways involved in the toxicological mechanism. This assay can be implemented as a reliable compound screening platform for improved safety assessment in early compound development pipelines.

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.

Predicting acute paraquat toxicity using physiologically based kinetic modelling incorporating in vitro active renal excretion via the OCT2 transporter

Including active renal excretion in physiologically based kinetic (PBK) models can improve their use in quantitative in vitro- in vivo extrapolation (QIVIVE) as a new approach methodology (NAM) for predicting the acute toxicity of organic cation transporter 2 (OCT2) substrates like paraquat (PQ). To realise this NAM, kinetic parameters Vmax and Km for in vitro OCT2 transport of PQ were obtained from the literature. Appropriate scaling factors were applied to translate the in vitro Vmax to an in vivo Vmax. in vitro cytotoxicity data were defined in the rat RLE-6TN and L2 cell lines and the human A549 cell line. The developed PQ PBK model was used to apply reverse dosimetry for QIVIVE translating the in vitro cytotoxicity concentration-response curves to predicted in vivo toxicity dose-response curves after which the lower and upper bound benchmark dose (BMD) for 50% lethality (BMDL50 and BMDU50) were derived by applying BMD analysis. Comparing the predictions to the in vivo reported LD50 values resulted in a conservative prediction for rat and a comparable prediction for human showing proof of principle on the inclusion of active renal excretion and prediction of PQ acute toxicity for the developed NAM.

New approach methodologies in human regulatory toxicology - Not if, but how and when!

The predominantly animal-centric approach of chemical safety assessment has increasingly come under pressure. Society is questioning overall performance, sustainability, continued relevance for human health risk assessment and ethics of this system, demanding a change of paradigm. At the same time, the scientific toolbox used for risk assessment is continuously enriched by the development of "New Approach Methodologies" (NAMs). While this term does not define the age or the state of readiness of the innovation, it covers a wide range of methods, including quantitative structure-activity relationship (QSAR) predictions, high-throughput screening (HTS) bioassays, omics applications, cell cultures, organoids, microphysiological systems (MPS), machine learning models and artificial intelligence (AI). In addition to promising faster and more efficient toxicity testing, NAMs have the potential to fundamentally transform today's regulatory work by allowing more human-relevant decision-making in terms of both hazard and exposure assessment. Yet, several obstacles hamper a broader application of NAMs in current regulatory risk assessment. Constraints in addressing repeated-dose toxicity, with particular reference to the chronic toxicity, and hesitance from relevant stakeholders, are major challenges for the implementation of NAMs in a broader context. Moreover, issues regarding predictivity, reproducibility and quantification need to be addressed and regulatory and legislative frameworks need to be adapted to NAMs. The conceptual perspective presented here has its focus on hazard assessment and is grounded on the main findings and conclusions from a symposium and workshop held in Berlin in November 2021. It intends to provide further insights into how NAMs can be gradually integrated into chemical risk assessment aimed at protection of human health, until eventually the current paradigm is replaced by an animal-free "Next Generation Risk Assessment" (NGRA).

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.

In vitro models to measure effects on intestinal deconjugation and transport of mixtures of bile acids

Bile acid metabolism and transport are critical to maintain bile acid homeostasis and host health. In this study, it was investigated if effects on intestinal bile acid deconjugation and transport can be quantified in vitro model systems using mixtures of bile acids instead of studying individual bile acids. To this end deconjugation of mixtures of selected bile acids in anaerobic rat or human fecal incubations and the effect of the antibiotic tobramycin on these reactions was studied. In addition, the effect of tobramycin on the transport of the bile acids in isolation or in a mixture across Caco-2 cell layers was characterized. The results demonstrate that both the reduction of bile acid deconjugation and transport by tobramycin can be adequately detected in vitro systems using a mixture of bile acids, thus eliminating the need to characterize the effects for each bile acid in separate experiments. Subtle differences between the experiments with single or combined bile acids point at mutual competitive interactions and indicate that the use of bile acid mixtures is preferred over use of single bile acid given that also in vivo bile acids occurs in mixtures.

Connecting Gut Microbial Diversity with Plasma Metabolome and Fecal Bile Acid Changes Induced by the Antibiotics Tobramycin and Colistin Sulfate

The diversity of microbial species in the gut has a strong influence on health and development of the host. Further, there are indications that the variation in expression of gut bacterial metabolic enzymes is less diverse than the taxonomic profile, underlying the importance of microbiome functionality, particularly from a toxicological perspective. To address these relationships, the gut bacterial composition of Wistar rats was altered by a 28 day oral treatment with the antibiotics tobramycin or colistin sulfate. On the basis of 16S marker gene sequencing data, tobramycin was found to cause a strong reduction in the diversity and relative abundance of the microbiome, whereas colistin sulfate had only a marginal impact. Associated plasma and fecal metabolomes were characterized by targeted mass spectrometry-based profiling. The fecal metabolome of tobramycin-treated animals had a high number of significant alterations in metabolite levels compared to controls, particularly in amino acids, lipids, bile acids (BAs), carbohydrates, and energy metabolites. The accumulation of primary BAs and significant reduction of secondary BAs in the feces indicated that the microbial alterations induced by tobramycin inhibit bacterial deconjugation reactions. The plasma metabolome showed less, but still many alterations in the same metabolite groups, including reductions in indole derivatives and hippuric acid, and furthermore, despite marginal effects of colistin sulfate treatment, there were nonetheless systemic alterations also in BAs. Aside from these treatment-based differences, we also uncovered interindividual differences particularly centering on the loss of Verrucomicrobiaceae in the microbiome, but with no apparent associated metabolite alterations. Finally, by comparing the data set from this study with metabolome alterations in the MetaMapTox database, key metabolite alterations were identified as plasma biomarkers indicative of altered gut microbiomes resulting from a wide activity spectrum of antibiotics.

Gut Microbiota as Well as Metabolomes of Wistar Rats Recover within Two Weeks after Doripenem Antibiotic Treatment

An understanding of the changes in gut microbiome composition and its associated metabolic functions is important to assess the potential implications thereof on host health. Thus, to elucidate the connection between the gut microbiome and the fecal and plasma metabolomes, two poorly bioavailable carbapenem antibiotics (doripenem and meropenem), were administered in a 28-day oral study to male and female Wistar rats. Additionally, the recovery of the gut microbiome and metabolomes in doripenem-exposed rats were studied one and two weeks after antibiotic treatment (i.e., doripenem-recovery groups). The 16S bacterial community analysis revealed an altered microbial population in all antibiotic treatments and a recovery of bacterial diversity in the doripenem-recovery groups. A similar pattern was observed in the fecal metabolomes of treated animals. In the recovery group, particularly after one week, an over-compensation was observed in fecal metabolites, as they were significantly changed in the opposite direction compared to previously changed metabolites upon 28 days of antibiotic exposure. Key plasma metabolites known to be diagnostic of antibiotic-induced microbial shifts, including indole derivatives, hippuric acid, and bile acids were also affected by the two carbapenems. Moreover, a unique increase in the levels of indole-3-acetic acid in plasma following meropenem treatment was observed. As was observed for the fecal metabolome, an overcompensation of plasma metabolites was observed in the recovery group. The data from this study provides insights into the connectivity of the microbiome and fecal and plasma metabolomes and demonstrates restoration post-antibiotic treatment not only for the microbiome but also for the metabolomes. The importance of overcompensation reactions for health needs further studies.

In vitro models to detect in vivo bile acid changes induced by antibiotics

Bile acid homeostasis plays an important role in many biological activities through the bile-liver-gut axis. In this study, two in vitro models were applied to further elucidate the mode of action underlying reported in vivo bile acid changes induced by antibiotics (colistin sulfate, tobramycin, meropenem trihydrate, and doripenem hydrate). 16S rRNA analysis of rat fecal samples anaerobically incubated with these antibiotics showed that especially tobramycin induced changes in the gut microbiota. Furthermore, tobramycin was shown to inhibit the microbial deconjugation of taurocholic acid (TCA) and the transport of TCA over an in vitro Caco-2 cell layer used as a model to mimic intestinal bile acid reuptake. The effects induced by the antibiotics in the in vitro model systems provide novel and complementary insight explaining the effects of the antibiotics on microbiota and fecal bile acid levels upon 28-day in vivo treatment of rats. In particular, our results provide insight in the mode(s) of action underlying the increased levels of TCA in the feces upon tobramycin exposure. Altogether, the results of the present study provide a proof-of-principle on how in vitro models can be used to elucidate in vivo effects on bile acid homeostasis, and to obtain insight in the mode(s) of action underlying the effect of an antibiotic, in this case tobramycin, on bile acid homeostasis via effects on intestinal bile acid metabolism and reuptake.

Towards a science-based testing strategy to identify maternal thyroid hormone imbalance and neurodevelopmental effects in the progeny-part III: how is substance-mediated thyroid hormone imbalance in pregnant/lactating rats or their progeny related to neurodevelopmental effects?

This review investigated which patterns of thyroid- and brain-related effects are seen in rats upon gestational/lactational exposure to 14 substances causing thyroid hormone imbalance by four different modes-of-action (inhibition of thyroid peroxidase, sodium-iodide symporter and deiodinase activities, enhancement of thyroid hormone clearance) or to dietary iodine deficiency. Brain-related parameters included motor activity, cognitive function, acoustic startle response, hearing function, periventricular heterotopia, electrophysiology and brain gene expression. Specific modes-of-action were not related to specific patterns of brain-related effects. Based upon the rat data reviewed, maternal serum thyroid hormone levels do not show a causal relationship with statistically significant neurodevelopmental effects. Offspring serum thyroxine together with offspring serum triiodothyronine and thyroid stimulating hormone appear relevant to predict the likelihood for neurodevelopmental effects. Based upon the collated database, thresholds of ≥60%/≥50% offspring serum thyroxine reduction and ≥20% and statistically significant offspring serum triiodothyronine reduction indicate an increased likelihood for statistically significant neurodevelopmental effects; accuracies: 83% and 67% when excluding electrophysiology (and gene expression). Measurements of brain thyroid hormone levels are likely relevant, too. The extent of substance-mediated thyroid hormone imbalance appears more important than substance mode-of-action to predict neurodevelopmental impairment in rats. Pertinent research needs were identified, e.g. to determine whether the phenomenological offspring thyroid hormone thresholds are relevant for regulatory toxicity testing. The insight from this review shall be used to suggest a tiered testing strategy to determine whether gestational/lactational substance exposure may elicit thyroid hormone imbalance and potentially also neurodevelopmental effects.

Gut microbiome and plasma metabolome changes in rats after oral gavage of nanoparticles: sensitive indicators of possible adverse health effects

BACKGROUND

The oral uptake of nanoparticles is an important route of human exposure and requires solid models for hazard assessment. While the systemic availability is generally low, ingestion may not only affect gastrointestinal tissues but also intestinal microbes. The gut microbiota contributes essentially to human health, whereas gut microbial dysbiosis is known to promote several intestinal and extra-intestinal diseases. Gut microbiota-derived metabolites, which are found in the blood stream, serve as key molecular mediators of host metabolism and immunity.

RESULTS

Gut microbiota and the plasma metabolome were analyzed in male Wistar rats receiving either SiO₂ (1000 mg/kg body weight/day) or Ag nanoparticles (100 mg/kg body weight/day) during a 28-day oral gavage study. Comprehensive clinical, histopathological and hematological examinations showed no signs of nanoparticle-induced toxicity. In contrast, the gut microbiota was affected by both nanoparticles, with significant alterations at all analyzed taxonomical levels. Treatments with each of the nanoparticles led to an increased abundance of Prevotellaceae, a family with gut species known to be correlated with intestinal inflammation. Only in Ag nanoparticle-exposed animals, Akkermansia, a genus known for its protective impact on the intestinal barrier was depleted to hardly detectable levels. In SiO₂ nanoparticles-treated animals, several genera were significantly reduced, including probiotics such as Enterococcus. From the analysis of 231 plasma metabolites, we found 18 metabolites to be significantly altered in Ag-or SiO₂ nanoparticles-treated rats. For most of these metabolites, an association with gut microbiota has been reported previously. Strikingly, both nanoparticle-treatments led to a significant reduction of gut microbiota-derived indole-3-acetic acid in plasma. This ligand of the arylhydrocarbon receptor is critical for regulating immunity, stem cell maintenance, cellular differentiation and xenobiotic-metabolizing enzymes.

CONCLUSIONS

The combined profiling of intestinal microbiome and plasma metabolome may serve as an early and sensitive indicator of gut microbiome changes induced by orally administered nanoparticles; this will help to recognize potential adverse effects of these changes to the host.

Predictive Performance of Next Generation Physiologically Based Kinetic (PBK) Model Predictions in Rats Based on In Vitro and In Silico Input Data

The goal of the present study was to assess the predictive performance of a minimal generic rat physiologically based kinetic (PBK) model based on in vitro and in silico input data to predict peak plasma concentrations (Cmax) upon single oral dosing. To this purpose, a dataset was generated of 3960 Cmax predictions for 44 compounds, applying different combinations of in vitro and in silico approaches for chemical parameterization, and comparison of the predictions to reported in vivo data. Best performance was obtained when (1) the hepatic clearance was parameterized based on in vitro measured intrinsic clearance values, (2) the method of Rodgers and Rowland for calculating partition coefficients, and (3) in silico calculated fraction unbound plasma and Papp values (the latter especially for very lipophilic compounds). Based on these input data, the median Cmax of 32 compounds could be predicted within 10-fold of the observed Cmax, with 22 out of these 32 compounds being predicted within 5-fold, and 8 compounds within 2-fold. Overestimations of more than 10-fold were observed for 12 compounds, whereas no underestimations of more than 10-fold occurred. Median Cmax predictions were frequently found to be within 10-fold of the observed Cmax when the scaled unbound hepatic intrinsic clearance (Clint,u) was either higher than 20 l/h or lower than 1 l/h. Similar findings were obtained with a test set of 5 in-house BASF compounds. Overall, this study provides relevant insights in the predictive performance of a minimal PBK model based on in vitro and in silico input data.

Use of physiologically based kinetic modeling-facilitated reverse dosimetry to predict in vivo acute toxicity of tetrodotoxin in rodents

In this study, the ability of a new in vitro/in silico quantitative in vitro–in vivo extrapolation (QIVIVE) methodology was assessed to predict the in vivo neurotoxicity of tetrodotoxin (TTX) in rodents. In vitro concentration–response data of TTX obtained in a multielectrode array assay with primary rat neonatal cortical cells and in an effect study with mouse neuro-2a cells were quantitatively extrapolated into in vivo dose–response data, using newly developed physiologically based kinetic (PBK) models for TTX in rats and mice. Incorporating a kidney compartment accounting for active renal excretion in the PBK models proved to be essential for its performance. To evaluate the predictions, QIVIVE-derived dose–response data were compared with in vivo data on neurotoxicity in rats and mice upon oral and parenteral dosing. The results revealed that for both rats and mice the predicted dose–response data matched the data from available in vivo studies well. It is concluded that PBK modeling-based reserve dosimetry of in vitro TTX effect data can adequately predict the in vivo neurotoxicity of TTX in rodents, providing a novel proof-of-principle for this methodology.

Ontogeny of renal, hepatic, and placental expression of ATP-binding cassette and solute carrier transporters in the rat and the rabbit

Species differences in developmental toxicity can be due to varying expression of xenobiotic transporters. Hence, knowledge on the ontogeny of these transporters, especially in human, rat and rabbit, is pivotal. Two superfamilies of transporters, the ATP-binding cassette (ABC) and the solute carrier (SLC) transporters, are well known for their role in the absorption, distribution and/or elimination of xenobiotics and endogenous substances. The aim of this study was to compare the expression levels of these xenobiotic transporters in liver, kidney and placenta of man, Wistar rat and New Zealand White rabbit during pre- and postnatal development. For this purpose, qPCR experiments were performed for rat and rabbit tissues and the gene expression profiles were compared with literature data from man, rat and rabbit. Data analysis showed large differences in transporter expression in development and between species. These results can be used to better understand developmental toxicity findings in non-clinical species and their relevance for man.

Uterotrophic and in vitro screening for (anti)estrogenic activity of dipyrone

Dipyrone is a commonly used analgesic in many countries and there is limited data on its possible endocrine disrupting effects. We performed a screening for in vivo and in vitro anti(estrogenic) activity of dipyrone. For the in vivo uterotrophic assay, immature female rats (22-days-old) were treated daily by oral gavage for three days with different doses of dipyrone alone (50, 100, 200 mg/kg/day) and associated with three ethynylestradiol (EE) doses (1, 3 and 10 μg/kg/day), which were based on a dose-response curve experiment. The uterine weight was used as a biomarker for estrogenicity. In a parallel in vitro approach, we used a yeast-based transcriptional activation reporter gene assay (Yeast Estrogen Screening - YES) for assessment of estrogenic agonistic and antagonistic effects of dipyrone and its main metabolites 4-methylaminoantipyrine (MAA) and 4-aminoantipyrine (AA). In the uterotrophic assay, animals that received EE at 1, 3 and 10 μg/kg/day showed an increase in relative uterine weight compared with vehicle-only rats (canola oil). Dipyrone did not increase uterine weight at any dose tested (50, 100 and 200 mg/kg/day) in relation to vehicle control, indicating absence of estrogenic activity. Furthermore, co-administration of dipyrone (50 and 200 mg/kg/day) and EE (1, 3 or 10 μg/kg/day) was unable to block EE estrogenic action in comparison to the groups treated with EE alone, indicating absence of antiestrogenic activity. In the YES assay dipyrone and its metabolites did not demonstrate estrogen agonistic or antagonistic properties in the yeast cells. These results suggest that dipyrone and its metabolites do not produce (anti)estrogenic effects in vivo or in vitro.

Use of in vitro metabolomics in NRK cells to help predicting nephrotoxicity and differentiating the MoA of nephrotoxicants

We describe a strategy using an in vitro metabolomics assay with tubular rat NRK-52E cells to investigate the Modes of Action (MoAs) of nephrotoxic compounds. Chemicals were selected according to their MoAs based on literature information: acetaminophen, 4-aminophenol and S-(trichlorovinyl-)L-cysteine (TCVC), (covalent protein binding); gentamycin, vancomycin, polymycin B and CdCl₂ (lysosomal overload) and tenofovir and cidofovir (mitochondrial DNA-interaction). After treatment and harvesting of the cells, intracellular endogenous metabolites were quantified relative to vehicle control. Metabolite patterns were evaluated in a purely data-driven pattern generation process excluding published information. This strategy confirmed the assignment of the chemicals to the respective MoA except for TCVC and CdCl₂. Finally, TCVC was defined as unidentified and CdCl₂ was reclassified to the MoA "covalent protein binding". Hierarchical cluster analysis of 58 distinct metabolites from the patterns enabled a clear visual separation of chemicals in each MoA. The assay reproducibility was very good and metabolic responses were consistent. These results support the use of metabolome analysis in NRK-52E cells as a suitable tool for understanding and investigating the MoA of nephrotoxicants. This assay could enable the early identification of nephrotoxic compounds and finally reduce animal testing.

Towards a science-based testing strategy to identify maternal thyroid hormone imbalance and neurodevelopmental effects in the progeny - part II: how can key events of relevant adverse outcome pathways be addressed in toxicological assessments?

The current understanding of thyroid-related adverse outcome pathways (AOPs) with adverse neurodevelopmental outcomes in mammals has been reviewed. This served to establish if standard rodent toxicity test methods and in vitro assays allow identifying thyroid-related modes-of-action potentially leading to adverse neurodevelopmental outcomes, and the human relevance of effects - in line with the European Commission's Endocrine Disruptor Criteria. The underlying hypothesis is that an understanding of the key events of relevant AOPs provides insight into differences in incidence, magnitude, or species sensitivity of adverse outcomes. The rodent studies include measurements of serum thyroid hormones, thyroid gland pathology and neurodevelopmental assessments, but do not directly inform on specific modes-of-action. Opportunities to address additional non-routine parameters reflecting critical events of AOPs in toxicological assessments are presented. These parameters appear relevant to support the identification of specific thyroid-related modes-of-action, provided that prevailing technical limitations are overcome. Current understanding of quantitative key event relationships is often weak, but would be needed to determine if the triggering of a molecular initiating event will ultimately result in an adverse outcome. Also, significant species differences in all processes related to thyroid hormone signalling are evident, but the biological implications thereof (including human relevance) are often unknown. In conclusion, careful consideration of the measurement (e.g. timing, method) and interpretation of additional non-routine parameters is warranted. These findings will be used in a subsequent paper to propose a testing strategy to identify if a substance may elicit maternal thyroid hormone imbalance and potentially also neurodevelopmental effects in the progeny.

Development of a Generic Physiologically Based Kinetic Model to Predict In Vivo Uterotrophic Responses Induced by Estrogenic Chemicals in Rats Based on In Vitro Bioassays

The present study assessed the potential of a generic physiologically based kinetic (PBK) model to convert in vitro data for estrogenicity to predict the in vivo uterotrophic response in rats for diethylstibestrol (DES), ethinylestradiol (EE2), genistein (GEN), coumestrol (COU), and methoxychlor (MXC). PBK models were developed using a generic approach and in vitro concentration-response data from the MCF-7 proliferation assay and the yeast estrogen screening assay were translated into in vivo dose-response data. Benchmark dose analysis was performed on the predicted data and available in vivo uterotrophic data to evaluate the model predictions. The results reveal that the developed generic PBK model adequate defines the in vivo kinetics of the estrogens. The predicted dose-response data of DES, EE2, GEN, COU, and MXC matched the reported in vivo uterus weight response in a qualitative way, whereas the quantitative comparison was somewhat hampered by the variability in both in vitro and in vivo data. From a safety perspective, the predictions based on the MCF-7 proliferation assay would best guarantee a safe point of departure for further risk assessment although it may be conservative. The current study indicates the feasibility of using a combination of in vitro toxicity data and a generic PBK model to predict the relative in vivo uterotrophic response for estrogenic chemicals.

Succinate dehydrogenase inhibitors: in silico flux analysis and in vivo metabolomics investigations show no severe metabolic consequences for rats and humans

Succinate dehydrogenase complex II inhibitors (SDHIs) are widely used fungicides since the 1960s. Recently, based on published in vitro cell viability data, potential health effects via disruption of the mitochondrial respiratory chain and tricarboxylic acid cycle have been postulated in mammalian species. As primary metabolic impact of SDH inhibition, an increase in succinate, and compensatory ATP production via glycolysis resulting in excess lactate levels was hypothesized. To investigate these hypotheses, genome-scale metabolic models of Rattus norvegicus and Homo sapiens were used for an in silico analysis of mammalian metabolism. Moreover, plasma samples from 28-day studies with the SDHIs boscalid and fluxapyroxad were subjected to metabolome analyses, to assess in vivo metabolite changes induced by SDHIs. The outcome of in silico analyses indicated that mammalian metabolic networks are robust and able to compensate different types of metabolic perturbation, e.g., partial or complete SDH inhibition. Additionally, the in silico comparison of rat and human responses suggested no noticeable differences between both species, evidencing that the rat is an appropriate testing organism for toxicity of SDHIs. Since no succinate or lactate accumulation were found in rats, such an accumulation is also not expected in humans as a result of SDHI exposure.

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.

Toward a science-based testing strategy to identify maternal thyroid hormone imbalance and neurodevelopmental effects in the progeny - part I: which parameters from human studies are most relevant for toxicological assessments?

The 2018 European Food Safety Authority/European Chemicals Agency Guidance on the Identification of Endocrine Disruptors lacks clarity on how the presence or absence of substance-induced maternal thyroid hormone imbalance, or the potential for subsequent deleterious consequences in child neurodevelopment, should be established by toxicological assessments. To address these uncertainties, this narrative review evaluates human evidence on how altered maternal thyroid function may be associated with child neurodevelopmental outcomes; and seeks to identify parameters in human studies that appear most relevant for toxicological assessments. Serum levels of free thyroxine (fT4) and thyroid stimulating hormone (TSH) are most frequently measured when assessing thyroid function in pregnant women, whereas a broad spectrum of neurodevelopmental parameters is used to evaluate child neurodevelopment. The human data confirms an association between altered maternal serum fT4 and/or TSH and increased risk for child neurodevelopmental impairment. Quantitative boundaries of effects indicative of increased risks need to be established. Moreover, it is unknown if altered serum levels of total T4, free or total triiodothyronine, or parameters unrelated to serum thyroid hormones might be more relevant indicators of such effects. None of the human studies established a link between substance-mediated liver enzyme induction and increased serum thyroid hormone clearance, let alone further to child neurodevelopmental impairment. This review identifies research needs to contribute to the development of toxicity testing strategies, to reliably predict whether substances have the potential to impair child neurodevelopment via maternal thyroid hormone imbalance.

A Workflow for the Performance of the Differential Ovarian Follicle Count Using Deep Neuronal Networks

In order to automate the counting of ovarian follicles required in multigeneration reproductive studies performed in the rat according to Organization for Economic Co-operation and Development guidelines 443 and 416, the application of deep neural networks was tested. The manual evaluation of the differential ovarian follicle count is a tedious and time-consuming task that requires highly trained personnel. In this regard, deep learning outputs provide overlay pictures for a more detailed documentation, together with an increased reproducibility of the counts. To facilitate the planned good laboratory practice (GLP) validation a workflow was set up using MLFlow to make all steps from generating of scans, training of the neural network, uploading of study images to the neural network, generation and storage of the results in a compliant manner controllable and reproducible. PyTorch was used as main framework to build the Faster region-based convolutional neural network for the training. We compared the performances of different depths of ResNet models with specific regard to the sensitivity, specificity, accuracy of the models. In this paper, we describe all steps from data labeling, training of networks, and the performance metrics chosen to evaluate different network architectures. We also make recommendation on steps, which should be taken into consideration when GLP validation is aimed for.

Antibiotic-Induced Changes in Microbiome-Related Metabolites and Bile Acids in Rat Plasma

Various environmental factors can alter the gut microbiome’s composition and functionality, and modulate host health. In this study, the effects of oral and parenteral administration of two poorly bioavailable antibiotics (i.e., vancomycin and streptomycin) on male Wistar Crl/Wi(Han) rats for 28 days were compared to distinguish between microbiome-derived or -associated and systemic changes in the plasma metabolome. The resulting changes in the plasma metabolome were compared to the effects of a third reference compound, roxithromycin, which is readily bioavailable. A community analysis revealed that the oral administration of vancomycin and roxithromycin in particular leads to an altered microbial population. Antibiotic-induced changes depending on the administration routes were observed in plasma metabolite levels. Indole-3-acetic acid (IAA) and hippuric acid (HA) were identified as key metabolites of microbiome modulation, with HA being the most sensitive. Even though large variations in the plasma bile acid pool between and within rats were observed, the change in microbiome community was observed to alter the composition of the bile acid pool, especially by an accumulation of taurine-conjugated primary bile acids. In-depth investigation of the relationship between microbiome variability and their functionality, with emphasis on the bile acid pool, will be necessary to better assess the potential adverseness of environmentally induced microbiome changes.

Body and organ weight data in 28-day toxicological studies in two mouse strains

Toxicological studies were performed in an AAALAC (American Association for Laboratory Animal Care)-approved laboratory at BASF SE, Ludwigshafen, Germany, in accordance with the German Animal Welfare Act and the effective European Council Directive 2010/63/EU. Data were recorded in the BASF SE pathology data capture system.

Historical control data (2008–2013) were compiled for

a) Twelve 28-day studies performed according to OECD 407 with mice from Janvier C57BL/j Rj (J) and Charles River CD-1 (CRL), in total 73 control females and 73 control males, 5–8 weeks old at the beginning of the studies. Data collected: terminal body weight, organ weights of adrenal glands, brain, epididymides, heart, kidneys, liver, ovaries, prostate, seminal vesicles (with coagulating glands), spleen, testes, thymus, uterus.

b) Eight immunotoxicity studies (duration of 28 days) performed according to EPA Health Effects Test Guidelines OPPTS 870.7800 with mice from Janvier C57BL/j Rj (J). 48 control females and 16 control males 5–7 weeks old at the beginning of the studies. Data collected: terminal body weight, mean absolute and relative weights of spleen and thymus.

This data helps interpreting effects caused by treatment in toxicology studies in the mouse.

Coefficients of variation were calculated and discussed in the accompanying research paper:

“Variance of body and organ weights in 28-day studies in mice” (Marxfeld et al. 2019).

Use cases, best practice and reporting standards for metabolomics in regulatory toxicology

Metabolomics is a widely used technology in academic research, yet its application to regulatory science has been limited. The most commonly cited barrier to its translation is lack of performance and reporting standards. The MEtabolomics standaRds Initiative in Toxicology (MERIT) project brings together international experts from multiple sectors to address this need. Here, we identify the most relevant applications for metabolomics in regulatory toxicology and develop best practice guidelines, performance and reporting standards for acquiring and analysing untargeted metabolomics and targeted metabolite data. We recommend that these guidelines are evaluated and implemented for several regulatory use cases.

The in vivo developmental toxicity of diethylstilbestrol (DES) in rat evaluated by an alternative testing strategy

In the present study, we evaluated an alternative testing strategy to quantitatively predict the in vivo developmental toxicity of the synthetic hormone diethylstilbestrol (DES). To this end, a physiologically based kinetic (PBK) model was defined that was subsequently used to translate concentration-response data for the in vitro developmental toxicity of DES, obtained in the ES-D3 cell differentiation assay, into predicted in vivo dose-response data for developmental toxicity. The previous studies showed that the PBK model-facilitated reverse dosimetry approach is a useful approach to quantitatively predict the developmental toxicity of several developmental toxins. The results obtained in the present study show that the PBK model adequately predicted DES blood concentrations in rats. Further studies revealed that DES tested positive in the ES-D3 differentiation assay and that DES-induced inhibition of the ES-D3 cell differentiation could be counteracted by the estrogen receptor alpha (ERα) antagonist fulvestrant, indicating that the in vitro ES-D3 cell differentiation assay was able to mimic the role of ERα reported in the mode of action underlying the developmental toxicity of DES in vivo. In spite of this, combining these in vitro data with the PBK model did not adequately predict the in vivo developmental toxicity of DES in a quantitative way. It is concluded that although the EST qualifies DES as a developmental toxin and detects the role of ERα in this process, the ES-D3 cell differentiation assay of the EST apparently does not adequately capture the processes underlying DES-induced developmental toxicity in vivo.

Intrinsic Xenobiotic Metabolizing Enzyme Activities in Early Life Stages of Zebrafish (Danio rerio)

Early life stages of zebrafish (Danio rerio, zf) are gaining attention as an alternative invivo test system for drug discovery, early developmental toxicity screenings and chemical testing in ecotoxicological and toxicological testing strategies. Previous studies have demonstrated transcriptional evidence for xenobiotic metabolizing enzymes (XME) during early zf development. However, elaborate experiments on XME activities during development are incomplete. In this work, the intrinsic activities of representative phase I and II XME were monitored by transformation of putative zf model substrates analyzed using photometry and high pressure liquid chromatography techniques. Six different defined stages of zf development (between 2.5 h postfertilization (hpf) to 120 hpf) were investigated by preparing a subcellular fraction from whole organism homogenates. We demonstrated that zf embryos as early as 2.5 hpf possess intrinsic metabolic activities for esterase, Aldh, Gst, and Cyp1a above the methodological detection limit. The activities of the enzymes Cyp3a and Nat were measurable during later stages in development. Activities represent dynamic patterns during development. The role of XME activities revealed in this work is relevant for the assessing toxicity in this test system and therefore contributes to a valuable characterization of zf embryos as an alternative testing organism in toxicology.

Hexamoll® DINCH: Lack of in vivo evidence for obesogenic properties

Hexamoll^® DINCH is an important alternative to phthalate plasticizers. Although regulatory reviews have not identified any potential hazards even in sensitive populations, an in vitro study by Campioli et al. (2015) suggested Hexamoll^® DINCH might alter fat storage in adipocytes resulting in obesity. To evaluate this hypothesis, data from studies with Hexamoll^® DINCH were reviewed for evidence of deposition in fat, changes in body weight, or changes in serum chemistry reflecting altered metabolic status. Body weights of F1 and F2 pups in a two-generation study did not differ from controls even at 1000 mg Hexamoll^® DINCH/kg body weight. Mean relative liver weights from the 1000 and 300 mg/kg bw groups were increased, but without histopathologic changes. Triglyceride and cholesterol levels in serum were not affected. In addition, subchronic and chronic studies in rats did not give evidence of an obesogenic effect. Radioactivity from 20 or 1000 mg/kg bw ¹⁴C-labelled Hexamoll^® DINCH dosed orally remained 2-3 times longer in adipose tissue than in well-perfused tissues; however, levels were 20-500% below other tissues at 1 and 8 h post dosing. Radioactivity concentrations in organs and tissues excluding the GI tract declined rapidly and continuously, and decreased in parallel to the concentration in plasma during the following 20 h. Both, initial and terminal half-lives of radioactivity concentration do not indicate a potential for accumulation. Furthermore, a metabolomic comparison of Hexamoll^® DINCH with DEHP and other phthalates shows complete separation of the metabolomic profile of these two chemical classes, meaning that their effects on the body and the body's reaction to the substance are different. Hence, comprehensive in vivo data do not show any evidence of Hexamoll^® DINCH altering fat metabolism or having obesogenic properties.

Towards a generic physiologically based kinetic model to predict in vivo uterotrophic responses in rats by reverse dosimetry of in vitro estrogenicity data

Physiologically based kinetic (PBK) modelling-based reverse dosimetry is a promising tool for the prediction of in vivo developmental toxicity using in vitro concentration–response data. In the present study, the potential of this approach to predict the dose-dependent increase of uterus weight in rats upon exposure to estrogenic chemicals was assessed. In vitro concentration–response data of 17β-estradiol (E2) and bisphenol A (BPA) obtained in the MCF-7/BOS proliferation assay, the U2OS ER-CALUX assay and the yeast estrogen screen (YES) assay, were translated into in vivo dose–response data in rat, using a PBK model with a minimum number of in vitro and in silico determined parameter values. To evaluate the predictions made, benchmark dose (BMD) analysis was performed on the predicted dose–response data and the obtained BMDL₁₀ values were compared with BMDL₁₀ values derived from data on the effects of E2 and BPA in the uterotrophic assay reported in the literature. The results show that predicted dose–response data of E2 and BPA matched with the data from in vivo studies when predictions were made based on YES assay data. The YES assay-based predictions of the BMDL₁₀ values differed 3.9-fold (E2) and 4.7- to 13.4-fold (BPA) from the BMDL₁₀ values obtained from the in vivo data. The present study provides the proof-of-principle that PBK modelling-based reverse dosimetry of YES assay data using a minimum PBK model can predict dose-dependent in vivo uterus growth caused by estrogenic chemicals. In future studies, the approach should be extended to include other estrogens.

Investigations on the dose-response relationship of combined exposure to low doses of three anti-androgens in Wistar rats

The current investigation examines whether combined exposure to three anti-androgens (flutamide, prochloraz, vinclozolin) result in interference with endocrine homeostasis when applied at very low dose levels, and whether the results of combined exposure are more pronounced than to the individual compounds. A pre-post-natal in vivo study design was chosen with more parameters than regulatory testing protocols require (additional endpoints addressing hormone levels, morphology and histopathological examinations). Dose levels were chosen to represent the lowest observed adverse effect level (LOAEL), the no observed adverse effect level (NOAEL), and the acceptable daily intake for each individual substance. Anti-androgenic changes were observable at the effect level (LOAEL) but not at lower exposures. Nipple/areola counts appeared to be a sensitive measure of effect, in addition to male sex organ weights at sexual maturation, and finally gross findings. The results indicate the absence of evidence for effects at low or very low dose levels. No (adverse) effects were seen at the NOAEL dose. A non-monotonic dose-response relationship was not evident. Combined exposure at LOAEL level resulted in enhanced responses for anogenital index, number of areolas/nipples, delayed preputial separation and reduced ventral prostate weight in comparison to the individual compounds.