Linking Mechanistic Effects of Pharmaceuticals and Personal Care Products to Ecologically Relevant Outcomes: A Decade of Progress

There are insufficient toxicity data to assess the ecological risks of many pharmaceuticals and personal care products (PPCPs). While data limitations are not uncommon for contaminants of environmental concern, PPCPs are somewhat unique in that an a priori understanding of their biological activities in conjunction with measurements of molecular, biochemical, or histological responses could provide a foundation for understanding mode(s) of action and predicting potential adverse apical effects. Over the past decade significant progress has been made in the development of new approach methodologies (NAMs) to efficiently quantify these types of endpoints using computational models and pathway-based in vitro and in vivo assays. The availability of open-access knowledgebases to curate biological response (including NAM) data and sophisticated bioinformatics tools to help interpret the information also has significantly increased. Finally, advances in the development and implementation of the adverse outcome pathway framework provide the critical conceptual underpinnings needed to translate NAM data into predictions of the ecologically relevant outcomes required by risk assessors and managers. The evolution and convergence of these various data streams, tools, and concepts provides the basis for a fundamental change in how ecological risks of PPCPs can be pragmatically assessed. Environ Toxicol Chem 2024;43:537-548. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Comparison of in silico, in vitro, and in vivo toxicity benchmarks suggests a role for ToxCast data in ecological hazard assessment

Large repositories of in vitro bioactivity data such as US EPA's Toxicity Forecaster (ToxCast) provide a wealth of publicly accessible toxicity information for thousands of chemicals. These data can be used to calculate point-of-departure (POD) estimates via concentration-response modeling that may serve as lower bound, protective estimates of in vivo effects. However, the data are predominantly based on mammalian models and discussions to date about their utility have largely focused on potential integration into human hazard assessment, rather than application to ecological risk assessment. The goal of the present study was to compare PODs based on (1) quantitative structure-activity relationships (QSARs), (2) the 5th centile of the activity concentration at cutoff (ACC), and (3) lower-bound cytotoxic burst (LCB) from ToxCast, with the distribution of in vivo PODs compiled in the Ecotoxicology Knowledgebase (ECOTOX). While overall correlation between ToxCast ACC5 and ECOTOX PODs for 649 chemicals was weak, there were significant associations among PODs based on LCB and ECOTOX, LCB and QSARs, and ECOTOX and QSARs. Certain classes of compounds showed moderate correlation across datasets (eg, antimicrobials/disinfectants), while others, such as organophosphate insecticides, did not. Unsurprisingly, more precise classifications of the data based on ECOTOX effect and endpoint type (eg, apical vs biochemical; acute vs chronic) had a significant effect on overall relationships. Results of this research help to define appropriate roles for data from new approach methodologies in chemical prioritization and screening of ecological hazards.

Verification of In Vivo Estrogenic Activity for Four Per- and Polyfluoroalkyl Substances (PFAS) Identified as Estrogen Receptor Agonists via New Approach Methodologies

Given concerns about potential toxicological hazards of the thousands of data-poor per- and polyfluorinated alkyl substances (PFAS) currently in commerce and detected in the environment, tiered testing strategies that employ high-throughput in vitro screening as an initial testing tier have been implemented. The present study evaluated the effectiveness of previous in vitro screening for identifying PFAS capable, or incapable, of inducing estrogenic responses in fish exposed in vivo. Fathead minnows (Pimephales promelas) were exposed for 96 h to five PFAS (perfluorooctanoic acid [PFOA]; 1H,1H,8H,8H-perfluorooctane-1,8-diol [FC8-diol]; 1H,1H,10H,10H-perfluorodecane-1,10-diol [FC10-diol]; 1H,1H,8H,8H-perfluoro-3,6-dioxaoctane-1,8-diol [FC8-DOD]; and perfluoro-2-methyl-3-oxahexanoic acid [HFPO-DA]) that showed varying levels of in vitro estrogenic potency. In agreement with in vitro screening results, exposure to FC8-diol, FC10-diol, and FC8-DOD caused concentration-dependent increases in the expression of transcript coding for vitellogenin and estrogen receptor alpha and reduced expression of insulin-like growth factor and apolipoprotein eb. Once differences in bioconcentration were accounted for, the rank order of potency in vivo matched that determined in vitro. These results provide a screening level benchmark for worst-case estimates of potential estrogenic hazards of PFAS and a basis for identifying structurally similar PFAS to scrutinize for putative estrogenic activity.

Leveraging ToxCast data and protein sequence conservation to complement aquatic life criteria derivation

The USEPA's 1985 guidelines for the derivation of aquatic life criteria (ALC) are robust but data-intensive. For many chemicals, the extensive in vivo data sets required for ALC derivation are not available. Thus, alternative analyses and processes that can provide provisional values to guide states, tribes, and other stakeholders while data accumulate and more rigorous criteria are derived would be beneficial. The overarching purpose of this study was to assess the feasibility of using data from new approach methodologies (NAMs) like ToxCast to derive first-pass, provisional values to guide chemical prioritization and resource management as a complement to traditional ALC derivation. To address this goal, the study objectives were to (1) estimate chemical potency using data from NAMs for nine compounds with available aquatic benchmarks, (2) evaluate the utility of using NAM data to elucidate potential mechanisms of toxicity to guide problem formulation, and (3) determine the species relevance of toxicity pathways for compounds with clearly defined mechanisms of action as a means to evaluate whether minimum data requirements could potentially be waived when deriving a more formal ALC. Points of departure were derived from ToxCast data based on the fifth percentile of the distribution of activity concentration above cutoff values falling below the cytotoxic burst. Mechanistic inferences were made based on active target hits in ToxCast and, where applicable, assessed for taxonomic conservation using SeqAPASS. ToxCast-based point-of-departure aligned relatively closely (six of nine test chemicals within a factor of 10; eight of nine within a factor of 100) with aquatic benchmarks from the USEPA and US Department of Energy (DOE). Moreover, pathways of toxicity gleaned from NAM data were reflective of in vivo-based findings from the literature. These results, while preliminary, and based on a limited number of substances, support the potential application of NAM data to complement traditional ALC derivation approaches and prioritization. Integr Environ Assess Manag 2023;19:224-238. © 2022 Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Persistence of improved glucose homeostasis in Gclm null mice with age and cadmium treatment

Antioxidant signaling/communication is among the most important cellular defense and survival pathways, and the importance of redox signaling and homeostasis in aging has been well-documented. Intracellular levels of glutathione (GSH), a very important endogenous antioxidant, both govern and are governed by the Nrf2 pathway through expression of genes involved in its biosynthesis, including the subunits of the rate-limiting enzyme (glutamate cysteine ligase, GCL) in GSH production, GCLC and GCLM. Mice homozygous null for the Gclm gene are severely deficient in GSH compared to wild-type controls, expressing approximately 10% of normal GSH levels. To compensate for GSH deficiency, Gclm null mice have upregulated redox-regulated genes, and, surprisingly, are less susceptible to certain types of oxidative damage. Furthermore, young Gclm null mice display an interesting lean phenotype, resistance to high fat diet-induced diabetes and obesity, improved insulin and glucose tolerance, and decreased expression of genes involved in lipogenesis. However, the persistence of this phenotype has not been investigated into old age, which is important in light of studies which suggest aging attenuates antioxidant signaling, particularly in response to exogenous stimuli. In this work, we addressed whether aging compromises the favorable phenotype of increased antioxidant activity and improved glucose homeostasis observed in younger Gclm null mice. We present data showing that under basal conditions and in response to cadmium exposure (2 mg/kg, dosed once via intraperitoneal injection), the phenotype previously described in young (<6 months) Gclm null mice persists into old age (24+ months). We also provide evidence that transcriptional activation of the Nrf2, AMPK, and PPARγ pathways underlie the favorable metabolic phenotype observed previously in young Gclm null mice.

Carbonyl Reductase 1 Plays a Significant Role in Converting Doxorubicin to Cardiotoxic Doxorubicinol in Mouse Liver, but the Majority of the Doxorubicinol-Forming Activity Remains Unidentified

Doxorubicin is a widely used cancer therapeutic, but its effectiveness is limited by cardiotoxic side effects. Evidence suggests cardiotoxicity is due not to doxorubicin, but rather its metabolite, doxorubicinol. Identification of the enzymes responsible for doxorubicinol formation is important in developing strategies to prevent cardiotoxicity. In this study, the contributions of three murine candidate enzymes to doxorubicinol formation were evaluated: carbonyl reductase (Cbr) 1, Cbr3, and thioredoxin reductase 1 (Tr1). Analyses with purified proteins revealed that all three enzymes catalyzed doxorubicin-dependent NADPH oxidation, but only Cbr1 and Cbr3 catalyzed doxorubicinol formation. Doxorubicin-dependent NADPH oxidation by Tr1 was likely due to redox cycling. Subcellular fractionation results showed that doxorubicin-dependent redox cycling activity was primarily microsomal, whereas doxorubicinol-forming activity was exclusively cytosolic, as were all three enzymes. An immunoclearing approach was used to assess the contributions of the three enzymes to doxorubicinol formation in the complex milieu of the cytosol. Immunoclearing Cbr1 eliminated 25% of the total doxorubicinol-forming activity in cytosol, but immunoclearing Cbr3 had no effect, even in Tr1 null livers that overexpressed Cbr3. The immunoclearing results constituted strong evidence that Cbr1 contributed to doxorubicinol formation in mouse liver but that enzymes other than Cbr1 also played a role, a conclusion supported by ammonium sulfate fractionation results, which showed that doxorubicinol-forming activity was found in fractions that contained little Cbr1. In conclusion, the results show that Cbr1 accounts for 25% of the doxorubicinol-forming activity in mouse liver cytosol but that the majority of the doxorubicinol-forming activity remains unidentified. SIGNIFICANCE STATEMENT: Earlier studies suggested carbonyl reductase (Cbr) 1 plays a dominant role in converting chemotherapeutic doxorubicin to cardiotoxic doxorubicinol, but a new immunoclearing approach described herein shows that Cbr1 accounts for only 25% of the doxorubicinol-forming activity in mouse liver cytosol, that two other candidate enzymes-Cbr3 and thioredoxin reductase 1-play no role, and that the majority of the activity remains unidentified. Thus, targeting Cbr1 is necessary but not sufficient to eliminate doxorubicinol-associated cardiotoxicity; identification of the additional doxorubicinol-forming activity is an important next challenge.

Kinetics of Glutathione Depletion and Antioxidant Gene Expression as Indicators of Chemical Modes of Action Assessed in Vitro in Mouse Hepatocytes with Enhanced Glutathione Synthesis

Here we report a vertically integrated in vitro - in silico study that aims to elucidate the molecular initiating events involved in the induction of oxidative stress (OS) by seven diverse chemicals (cumene hydroperoxide, t-butyl hydroperoxide, hydroquinone, t-butyl hydroquinone, bisphenol A, Dinoseb, and perfluorooctanoic acid). To that end, we probe the relationship between chemical properties, cell viability, glutathione (GSH) depletion, and antioxidant gene expression. Concentration-dependent effects on cell viability were assessed by MTT assay in two Hepa-1 derived mouse liver cell lines: a control plasmid vector transfected cell line (Hepa-V), and a cell line with increased glutamate-cysteine ligase (GCL) activity and GSH content (CR17). Changes to intracellular GSH content and mRNA expression levels for the Nrf2-driven antioxidant genes Gclc, Gclm, heme oxygenase-1 ( Hmox1), and NADPH quinone oxidoreductase-1 ( Nqo1) were monitored after sublethal exposure to the chemicals. In silico models of covalent and redox reactivity were used to rationalize differences in activity of quinones and peroxides. Our findings show CR17 cells were generally more resistant to chemical toxicity and showed markedly attenuated induction of OS biomarkers; however, differences in viability effects between the two cell lines were not the same for all chemicals. The results highlight the vital role of GSH in protecting against oxidative stress-inducing chemicals as well as the importance of probing molecular initiating events in order to identify chemicals with lower potential to cause oxidative stress.

Functional kinomics identifies candidate therapeutic targets in head and neck cancer

PURPOSE

To identify novel therapeutic drug targets for p53-mutant head and neck squamous cell carcinoma (HNSCC).

EXPERIMENTAL DESIGN

RNAi kinome viability screens were performed on HNSCC cells, including autologous pairs from primary tumor and recurrent/metastatic lesions, and in parallel on murine squamous cell carcinoma (MSCC) cells derived from tumors of inbred mice bearing germline mutations in Trp53, and p53 regulatory genes: Atm, Prkdc, and p19(Arf). Cross-species analysis of cell lines stratified by p53 mutational status and metastatic phenotype was used to select 38 kinase targets. Both primary and secondary RNAi validation assays were performed on additional HNSCC cell lines to credential these kinase targets using multiple phenotypic endpoints. Kinase targets were also examined via chemical inhibition using a panel of kinase inhibitors. A preclinical study was conducted on the WEE1 kinase inhibitor, MK-1775.

RESULTS

Our functional kinomics approach identified novel survival kinases in HNSCC involved in G2-M cell-cycle checkpoint, SFK, PI3K, and FAK pathways. RNAi-mediated knockdown and chemical inhibition of the WEE1 kinase with a specific inhibitor, MK-1775, had a significant effect on both viability and apoptosis. Sensitivity to the MK-1775 kinase inhibitor is in part determined by p53 mutational status, and due to unscheduled mitotic entry. MK-1775 displays single-agent activity and potentiates the efficacy of cisplatin in a p53-mutant HNSCC xenograft model.

CONCLUSIONS

WEE1 kinase is a potential therapeutic drug target for HNSCC. This study supports the application of a functional kinomics strategy to identify novel therapeutic targets for cancer.

Metabolism of doxorubicin to the cardiotoxic metabolite doxorubicinol is increased in a mouse model of chronic glutathione deficiency: A potential role for carbonyl reductase 3

Doxorubicin is highly effective at inducing DNA double-strand breaks in rapidly dividing cells, which has led to it being a widely used cancer chemotherapeutic. However, clinical administration of doxorubicin is limited by off-target cardiotoxicity, which is thought to be mediated by doxorubicinol, the primary alcohol metabolite of doxorubicin. Carbonyl reductase 1 (CBR1), a well-characterized monomeric enzyme present at high basal levels in the liver, is known to exhibit activity toward doxorubicin. Little is known about a closely related enzyme, carbonyl reductase 3 (CBR3), which is present in the liver at low basal levels but is highly inducible by the transcription factor Nrf2. Genetic polymorphisms in CBR3, but not CBR1, are associated with differential cardiac outcomes in doxorubicin treated pediatric patients. Cbr3 mRNA and CBR3 protein are highly expressed in the livers of Gclm-/- mice (a mouse model of glutathione deficiency) relative to wild type mice. In the present study, we first investigated the ability of CBR3 to metabolize doxorubicin. Incubations of doxorubicin and purified recombinant murine CBR3 (mCBR3) were analyzed for doxorubicinol formation using HPLC, revealing for the first time that doxorubicin is a substrate of mCBR3. Moreover, hepatocytes from Gclm-/- mice produced more doxorubicinol than Gclm+/+ hepatocytes. In addition, differentiated rat myoblasts (C2C12 cells) co-cultured with primary Gclm-/- murine hepatocytes were more sensitive to doxorubicin-induced cytostasis/cytotoxicity than incubations with Gclm+/+ hepatocytes. Our results indicate a potentially important role for CBR3 in doxorubicin-induced cardiotoxicity. Because there is likely to be variability in hepatic CBR3 activity in humans (due to either genetic or epigenetic influences on its expression), these data also suggest that inhibition of CBR3 may provide protection from doxorubicinol cardiotoxicity.

Of mice, rats, and men: could Nrf2 activation protect against aflatoxin heptocarcinogenesis in humans?

In this issue, Johnson and colleagues provide a remarkable demonstration of the potential for "chemoprevention" of cancer from mutagenic chemicals. The authors demonstrated complete protection of rats from a carcinogenic treatment regimen with the potent dietary mutagen and hepatocarcinogen, aflatoxin B1 (AFB) by pretreatment with a synthetic oleanane triterpenoid, 1-[2-cyano-3-,12-dioxooleana-1,9(11)-dien-28-oyl]imidazole (CDDO-Im). This study is notable for two reasons: (i) Activation of the Nrf2/Keap1/ARE "antioxidant response" pathway by CDDO-Im conferred complete protection against AFB-induced hepatocellular carcinomas in the Fisher F344 rat (a strain frequently used in life-time carcinogenicity bioassays), and (ii) extensive AFB-DNA adduct formation was seen in all animals at early time points, including those treated with CDDO-Im, albeit at lower levels (∼30% of the untreated animals), suggesting a strong divergence in the association between early DNA-damaging events, and tumor formation later in life. The authors suggest that this provides compelling experimental support for the concept of carcinogenic "thresholds" for mutagenic chemicals, because the treatment reduced persistent, mutagenic adducts (AFB-FAPyr adducts) only by 70%, but nearly completely eliminated tumors after approximately 2 years and preneoplastic lesions 6 weeks after the last dose of AFB.