A model for aryl hydrocarbon receptor-activated gene expression shows potency and efficacy changes and predicts squelching due to competition for transcription co-activators

Connecting Suborganismal Data to Bioenergetic Processes: Killifish Embryos Exposed to a Dioxin-Like Compound

A core challenge for ecological risk assessment is to integrate molecular responses into a chain of causality to organismal or population-level outcomes. Bioenergetic theory may be a useful approach for integrating suborganismal responses to predict organismal responses that influence population dynamics. We describe a novel application of dynamic energy budget (DEB) theory in the context of a toxicity framework (adverse outcome pathways [AOPs]) to make quantitative predictions of chemical exposures to individuals, starting from suborganismal data. We use early-life stage exposure of Fundulus heteroclitus to dioxin-like chemicals (DLCs) and connect AOP key events to DEB processes through "damage" that is produced at a rate proportional to the internal toxicant concentration. We use transcriptomic data of fish embryos exposed to DLCs to translate molecular indicators of damage into changes in DEB parameters (damage increases somatic maintenance costs) and DEB models to predict sublethal and lethal effects on young fish. By changing a small subset of model parameters, we predict the evolved tolerance to DLCs in some wild F. heteroclitus populations, a data set not used in model parameterization. The differences in model parameters point to reduced sensitivity and altered damage repair dynamics as contributing to this evolved resistance. Our methodology has potential extrapolation to untested chemicals of ecological concern. Environ Toxicol Chem 2023;42:2040-2053. © 2023 Oak Ridge National Laboratory and The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Development, scrutiny, and modulation of transient reporter gene assays of the xenobiotic metabolism pathway in zebrafish hepatocytes

The "toxicology in the twenty-first century" paradigm shift demands the development of alternative in vitro test systems. Especially in the field of ecotoxicology, coverage of aquatic species-specific assays is relatively scarce. Transient reporter gene assays could be a quick, economical, and reliable bridging technology. However, the user should be aware of potential pitfalls that are influenced by reporter vector geometry. Here, we report the development of an AhR-responsive transient reporter-gene assay in the permanent zebrafish hepatocytes cell line (ZFL). Additionally, we disclose how viral, constitutive promoters within reporter-gene assay cassettes induce squelching of the primary signal. To counter this, we designed a novel normalization vector, bearing an endogenous zebrafish-derived genomic promoter (zfEF1aPro), which rescues the squelching-delimited system, thus, giving new insights into the modulation of transient reporter systems under xenobiotic stress. Finally, we uncovered how the ubiquitously used ligand BNF promiscuously activates multiple toxicity pathways of the xenobiotic metabolism and cellular stress response in an orchestral manner, presumably leading to a concentration-related inhibition of the AhR/ARNT/XRE-toxicity pathway and non-monotonous concentration-response curves. We named such a multi-level inhibitory mechanism that might mask effects as "maisonette squelching." A transient reporter gene assay in zebrafish cell lines utilizing endogenous regulatory gene elements shows increased in vitro toxicity testing performance. Synthetic and constitutive promotors interfere with signal transduction ("squelching") and might increase cellular stress (cytotoxicity). The squelching phenomenon might occur on multiple levels (toxicity pathway crosstalk and normalization vector), leading to a complete silencing of the reporter signal.

CYP1A1, VEGFA and Adipokine Responses of Human Adipocytes Co-exposed to PCB126 and Hypoxia

It is increasingly recognized that hypoxia may develop in adipose tissue as its mass expands. Adipose tissue is also the main reservoir of lipophilic pollutants, including polychlorinated biphenyls (PCBs). Both hypoxia and PCBs have been shown to alter adipose tissue functions. The signaling pathways induced by hypoxia and pollutants may crosstalk, as they share a common transcription factor: aryl hydrocarbon receptor nuclear translocator (ARNT). Whether hypoxia and PCBs crosstalk and affect adipokine secretion in human adipocytes remains to be explored. Using primary human adipocytes acutely co-exposed to different levels of hypoxia (24 h) and PCB126 (48 h), we observed that hypoxia significantly inhibits the PCB126 induction of cytochrome P450 (CYP1A1) transcription in a dose-response manner, and that Acriflavine (ACF)—an HIF1α inhibitor—partially restores the PCB126 induction of CYP1A1 under hypoxia. On the other hand, exposure to PCB126 did not affect the transcription of the vascular endothelial growth factor-A (VEGFA) under hypoxia. Exposure to hypoxia increased leptin and interleukin-6 (IL-6), and decreased adiponectin levels dose-dependently, while PCB126 increased IL-6 and IL-8 secretion in a dose-dependent manner. Co-exposure to PCB126 and hypoxia did not alter the adipokine secretion pattern observed under hypoxia and PCB126 exposure alone. In conclusion, our results indicate that (1) hypoxia inhibits PCB126-induced CYP1A1 expression at least partly through ARNT-dependent means, suggesting that hypoxia could affect PCB metabolism and toxicity in adipose tissue, and (2) hypoxia and PCB126 affect leptin, adiponectin, IL-6 and IL-8 secretion differently, with no apparent crosstalk between the two factors.

Discrete Biochemical Systems Theory

Almost every biomedical systems analysis requires early decisions regarding the choice of the most suitable representations to be used. De facto the most prevalent choice is a system of ordinary differential equations (ODEs). This framework is very popular because it is flexible and fairly easy to use. It is also supported by an enormous array of stand-alone programs for analysis, including many distinct numerical solvers that are implemented in the main programming languages. Having selected ODEs, the modeler must then choose a mathematical format for the equations. This selection is not trivial as nearly unlimited options exist and there is seldom objective guidance. The typical choices include ad hoc representations, default models like mass-action or Lotka-Volterra equations, and generic approximations. Within the realm of approximations, linear models are typically successful for analyses of engineered systems, but they are not as appropriate for biomedical phenomena, which often display nonlinear features such as saturation, threshold effects or limit cycle oscillations, and possibly even chaos. Power-law approximations are simple but overcome these limitations. They are the key ingredient of Biochemical Systems Theory (BST), which uses ODEs exclusively containing power-law representations for all processes within a model. BST models cover a vast repertoire of nonlinear responses and, at the same time, have structural properties that are advantageous for a wide range of analyses. Nonetheless, as all ODE models, the BST approach has limitations. In particular, it is not always straightforward to account for genuine discreteness, time delays, and stochastic processes. As a new option, we therefore propose here an alternative to BST in the form of discrete Biochemical Systems Theory (dBST). dBST models have the same generality and practicality as their BST-ODE counterparts, but they are readily implemented even in situations where ODEs struggle. As a case study, we illustrate dBST applied to the dynamics of the aryl hydrocarbon receptor (AhR), a signal transduction system that simultaneously involves time delays and stochasticity.

Androgen Receptor-Mediated Transcription in Prostate Cancer

Androgen receptor (AR)-mediated transcription is critical in almost all stages of prostate cancer (PCa) growth and differentiation. This process involves a complex interplay of coregulatory proteins, chromatin remodeling complexes, and other transcription factors that work with AR at cis-regulatory enhancer regions to induce the spatiotemporal transcription of target genes. This enhancer-driven mechanism is remarkably dynamic and undergoes significant alterations during PCa progression. In this review, we discuss the AR mechanism of action in PCa with a focus on how cis-regulatory elements modulate gene expression. We explore emerging evidence of genetic variants that can impact AR regulatory regions and alter gene transcription in PCa. Finally, we highlight several outstanding questions and discuss potential mechanisms of this critical transcription factor.

DNA methylation differentiates smoking from vaping and non-combustible tobacco use

Increasing use of non-combusted forms of nicotine such as e-cigarettes poses important public health questions regarding their specific risks relative to combusted tobacco products such as cigarettes. To fully delineate these risks, improved biomarkers that can distinguish between these forms of nicotine use are needed. Prior work has suggested that methylation status at cg05575921 may serve as a specific biomarker of combusted tobacco smoke exposure. We hypothesized combining this epigenetic biomarker with conventional metabolite assays could classify the type of nicotine product consumption. Therefore, we determined DNA methylation and serum cotinine values in samples from 112 smokers, 35 e-cigarette users, 19 smokeless tobacco users, and 269 controls, and performed mass spectroscopy analyses of urine samples from all nicotine users and 22 verified controls to determine urinary levels of putatively nicotine product-specific substances; propylene glycol, 2-cyanoethylmercapturic acid (CEMA), and anabasine. 1) Cigarette smoking was associated with a dose dependent demethylation of cg05575921 and increased urinary CEMA and anabasine levels, 2) e-cigarette use did not demethylate cg05575921, 3) smokeless tobacco use also did not demethylate cg05575921 but was positively associated with anabasine levels 4) CEMA and cg05575921 levels were highly correlated and 5) propylene glycol levels did not reliably distinguish use groups. Cg05575921 assessments distinguish exposure to tobacco smoke from smokeless sources of nicotine including e-cigarettes and smokeless tobacco, neither of which are associated with cg05575921 demethylation. A combination of methylomic and metabolite profiling may allow for accurate classification use status of a variety of nicotine containing products.

Potentials and pitfalls of transient in vitro reporter bioassays: interference by vector geometry and cytotoxicity in recombinant zebrafish cell lines

The water framework directive re-evaluation proposes the integration of effect-based tools, increasing the need for alternative methods. Especially within aquatic toxicology, coverage of specific toxicity pathways is scarce, and most applications are based on mammalian or bacterial models, not reflecting realistic exposure scenarios. The use of transient reporter gene assays in cells from organisms of interest could be a quick and inexpensive solution. However, interference with cellular homeostasis may impact the system beyond the function of the manipulated gene and thus lead to non-specific results. We describe how varying vector geometry and different regulatory gene elements on plasmids used for transfection in zebrafish hepatocytes and embryonic fibroblasts may lead up to a tenfold difference in potency. Cells were transiently co-transfected with an Nrf2-responsive Firefly luciferase reporter plasmid and eight different Renilla luciferase normalization plasmids. Transfected cells were exposed to two different regimes (0.1–100 µM and 7.8–250 µM) of the oxidative stress-inducing compounds, sulforaphane, tertbutylhydroquinone, and metazachlor. Nrf2 activity was measured in dual-luciferase assays. In parallel, cytotoxicity was assessed for different endpoints (energy metabolism, protein amount, membrane stability, and cell proliferation) in non-transfected cells and cells co-transfected with constructs of increasing size, to be used for normalization. Transfected cells were more susceptible to cytotoxicity in a vector size-dependent manner. Conclusively, we report that vector geometries (size, backbones, gene-regulatory units), cell line (tissue origin), applied transfection methods, and signal normalization may alter the sensitivity of reporter bioassays in a synergistic manner. Further, we propose that thorough bioassay design is needed to ensure reliability and regulatory acceptance.

Aryl Hydrocarbon Receptor Diet and Breast Cancer Risk

Breast cancer is the most common type of cancer and leading cause of cancer mortality among women worldwide. However, the majority of breast malignancies are of sporadic etiology. Therefore, identifying risk-mitigating factors may significantly decrease the burden of breast cancer. Diet can have both a predisposing and protective role in breast tumorigenesis. However, establishing efficacy of dietary constituents for cancer prevention has been limited by suboptimal dietary assessment. There is a need to acquire new experimental evidence that can be used to discriminate beneficial from harmful dietary constituents. The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that is recognized as the mediator of halogenated and polycyclic aromatic hydrocarbon toxicities. Importantly, evidence points to a breast tumor-promoting role for the AhR. Preclinical and clinical studies suggest that the AhR is overexpressed in advanced and triple negative breast cancers. Several dietary constituents, namely flavonoid compounds, have demonstrated inhibitory effects on AhR activation. Given this background, in this paper we elaborate on the working hypothesis that a diet rich in AhR food agonists favors breast tumor development, whereas a diet rich in AhR food antagonists is protective. As an initial approach to developing an AhR diet hypothesis, we conducted a review of published studies reporting on the association between intake of AhR inhibitory foods and risk of breast cancer. To assist the reader with interpretation of the concepts leading to the AhR diet hypothesis, we have preceded this review with an overview of AhR biology and its role in breast cancer development.

Aryl hydrocarbon receptor signaling modulates antiviral immune responses: ligand metabolism rather than chemical source is the stronger predictor of outcome

The aryl hydrocarbon receptor (AHR) offers a compelling target to modulate the immune system. AHR agonists alter adaptive immune responses, but the consequences differ across studies. We report here the comparison of four agents representing different sources of AHR ligands in mice infected with influenza A virus (IAV): TCDD, prototype exogenous AHR agonist; PCB126, pollutant with documented human exposure; ITE, novel pharmaceutical; and FICZ, degradation product of tryptophan. All four compounds diminished virus-specific IgM levels and increased the proportion of regulatory T cells. TCDD, PCB126 and ITE, but not FICZ, reduced virus-specific IgG levels and CD8⁺ T cell responses. Similarly, ITE, PCB126, and TCDD reduced Th1 and Tfh cells, whereas FICZ increased their frequency. In Cyp1a1-deficient mice, all compounds, including FICZ, reduced the response to IAV. Conditional Ahr knockout mice revealed that all four compounds require AHR within hematopoietic cells. Thus, differences in the immune response to IAV likely reflect variances in quality, magnitude, and duration of AHR signaling. This indicates that binding affinity and metabolism may be stronger predictors of immune effects than a compound’s source of origin, and that harnessing AHR will require finding a balance between dampening immune-mediated pathologies and maintaining sufficient host defenses against infection.

Integration of novel approaches demonstrates simultaneous metabolic inactivation and CAR-mediated hepatocarcinogenesis of a nitrification inhibitor

Nitrapyrin, a nitrification inhibitor, produces liver tumors in mice at high doses. Several experiments were performed to investigate molecular, cellular, and apical endpoints to define the key events leading to the tumor formation. These data support a mode-of-action (MoA) characterized by constitutive androstane receptor (CAR) nuclear receptor activation, increased hepatocellular proliferation leading to hepatocellular foci and tumor formation. Specifically, nitrapyrin induced a dose-related increase in the Cyp2b10/CAR-associated transcript and protein. Interestingly, the corresponding enzyme activity (7-pentoxyresorufin-O-dealkylase (PROD) was not enhanced due to nitrapyrin-mediated suicide inhibition of PROD activity. Nitrapyrin exposure elicited a clear dose-responsive increase in hepatocellular proliferation in wild-type mice, but not in CAR knock-out mice, informing that CAR activation is an obligatory key event in this test material-induced hepatocarcinogenesis. Furthermore, nitrapyrin exposure induced a clear, concentration-responsive increase in cell proliferation in mouse, but not human, hepatocytes in vitro. Evaluation of the data from repeat dose and MoA studies by the Bradford Hill criteria and a Human Relevance Framework (HRF) suggested that nitrapyrin-induced mouse liver tumors are not relevant to human health risk assessment because of qualitative differences between these two species.

How well can carcinogenicity be predicted by high throughput "characteristics of carcinogens" mechanistic data?

IARC has begun using ToxCast/Tox21 data in efforts to represent key characteristics of carcinogens to organize and weigh mechanistic evidence in cancer hazard determinations and this implicit inference approach also is being considered by USEPA. To determine how well ToxCast/Tox21 data can explicitly predict cancer hazard, this approach was evaluated with statistical analyses and machine learning prediction algorithms. Substances USEPA previously classified as having cancer hazard potential were designated as positives and substances not posing a carcinogenic hazard were designated as negatives. Then ToxCast/Tox21 data were analyzed both with and without adjusting for the cytotoxicity burst effect commonly observed in such assays. Using the same assignments as IARC of ToxCast/Tox21 assays to the seven key characteristics of carcinogens, the ability to predict cancer hazard for each key characteristic, alone or in combination, was found to be no better than chance. Hence, we have little scientific confidence in IARC's inference models derived from current ToxCast/Tox21 assays for key characteristics to predict cancer. This finding supports the need for a more rigorous mode-of-action pathway-based framework to organize, evaluate, and integrate mechanistic evidence with animal toxicity, epidemiological investigations, and knowledge of exposure and dosimetry to evaluate potential carcinogenic hazards and risks to humans.