The term “cytokine storm” describes an acute pathophysiologic state of the immune system characterized by a burst of cytokine release, systemic inflammatory response, and multiple organ failure, which are crucial determinants of many disease outcomes. In light of the complexity of cytokine storms, specific strategies are needed to prevent and alleviate their occurrence and deterioration. Nuclear factor erythroid 2-related factor 2 (NRF2) is a CNC-basic region-leucine zipper protein that serves as a master transcription factor in maintaining cellular redox homeostasis by orchestrating the expression of many antioxidant and phase II detoxification enzymes. Given that inflammatory response is intertwined with oxidative stress, it is reasonable to assume that NRF2 activation limits inflammation and thus cytokine storms. As NRF2 can mitigate inflammation at many levels, it has emerged as a potential target to prevent and treat cytokine storms. In this review, we summarized the cytokine storms caused by different etiologies and the rationale of interventions, focusing mainly on NRF2 as a potential therapeutic target.

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

Vitamin D has been linked to various physiological functions in pregnant women and their fetuses. Previous studies have suggested that some per- and polyfluoroalkyl substances (PFAS) may alter serum vitamin D concentrations. However, no study has investigated the relationship between PFAS and vitamin D in pregnant women. This study aims to evaluate the associations of serum PFAS with serum total and free 25-hydroxyvitamin D (25(OH)D) during pregnancy in a cohort of African American women in Atlanta, GA. Blood samples from 442 participants were collected in early pregnancy (8–14 weeks of gestation) for PFAS and 25(OH)D measurements, and additional samples were collected in late pregnancy (24–30 weeks) for the second 25(OH)D measurements. We fit multivariable linear regressions and weighted quantile sum (WQS) regressions to estimate the associations of individual PFAS and their mixtures with 25(OH)D concentrations. We found mostly positive associations of total 25(OH)D with PFHxS (perfluorohexane sulfonic acid), PFOS (perfluorooctane sulfonic acid), PFDA (perfluorodecanoic acid), and NMeFOSAA (N-methyl perfluorooctane sulfonamido acetic acid), and negative associations with PFPeA (perfluoropentanoic acid). For free 25(OH)D, positive associations were observed with PFHxS, PFOS, PFOA (perfluorooctanoic acid), and PFDA, and a negative association with PFPeA among the women with male fetuses in the models using 25(OH)D measured in late pregnancy. In mixture models, a quartile increase in WQS index was associated with 2.88 ng/mL (95%CI 1.14–4.59) and 5.68 ng/mL (95%CI 3.31–8.04) increases in total 25(OH)D measured in the early and late pregnancy, respectively. NMeFOSAA, PFDA, and PFOS contributed the most to the overall effects among the eight PFAS. No association was found between free 25(OH)D and the PFAS mixture. These results suggest that PFAS may affect vitamin D biomarker concentrations in pregnant African American women, and some of the associations were modified by fetal sex.

Thyroid hormones (THs), including T4 and T3, are produced and released by the thyroid gland under the stimulation of thyroid-stimulating hormone (TSH). The homeostasis of THs is regulated via the coordination of the hypothalamic-pituitary-thyroid axis, plasma binding proteins, and local metabolism in tissues. TH synthesis and secretion in the thyrocytes-containing thyroid follicles are exquisitely regulated by an elaborate molecular network comprising enzymes, transporters, signal transduction machineries, and transcription factors. In this article, we synthesized the relevant literature, organized and dissected the complex intrathyroidal regulatory network into structures amenable to functional interpretation and systems-level modeling. Multiple intertwined feedforward and feedback motifs were identified and described, centering around the transcriptional and posttranslational regulations involved in TH synthesis and secretion, including those underpinning the Wolff-Chaikoff and Plummer effects and thyroglobulin-mediated feedback regulation. A more thorough characterization of the intrathyroidal network from a systems biology perspective, including its topology, constituent network motifs, and nonlinear quantitative properties, can help us to better understand and predict the thyroidal dynamics in response to physiological signals, therapeutic interventions, and environmental disruptions.

Signal amplification in biomolecular networks converts a linear input to a steeply sigmoid output and is central to a number of cellular functions including proliferation, differentiation, homeostasis, adaptation, and biological rhythms. One canonical signal amplifying motif is zero-order ultrasensitivity that is mediated through the posttranslational modification (PTM) cycle of signaling proteins. The functionality of this signaling motif has been examined conventionally by supposing that the total amount of the protein substrates remains constant, as by the classical Koshland– Goldbeter model. However, covalent modification of signaling proteins often results in changes in their stability, which affects the abundance of the protein substrates. Here, we use mathematical models to explore the signal amplification properties in such scenarios and report some novel aspects. Our analyses indicate that PTM-induced protein stabilization brings the enzymes closer to saturation. As a result, ultrasensitivity may emerge or is greatly enhanced, with a steeper sigmoidal response, higher magnitude, and generally longer response time. In cases where PTM destabilizes the protein, ultrasensitivity can be regained through changes in the activities of the involved enzymes or from increased protein synthesis. Importantly, ultrasensitivity is not limited to modified or unmodified protein substrates—when protein turnover is considered, the total free protein substrate can also exhibit ultrasensitivity under several conditions. When full enzymatic reactions are used instead of Michaelis–Menten kinetics for the modeling, the total free protein substrate can even exhibit nonmonotonic dose–response patterns. It is conceivable that cells use inducible protein stabilization as a strategy in the signaling network to boost signal amplification while saving energy by keeping the protein substrate levels low at basal conditions.

Fine-tuning of osteoclast differentiation (OD) and bone remodeling is crucial for bone homeostasis. Dissecting the mechanisms regulating osteoclastogenesis is fundamental to understanding the pathogenesis of various bone disorders including osteoporosis and arthritis. Nuclear factor erythroid 2-related factor 1 (NFE2L1, also known as NRF1), which belongs to the CNC-bZIP family of transcription factors, orchestrates a variety of physiological processes and stress responses. While Nfe2l1 gene may be transcribed into multiple alternatively spliced isoforms, the biological function of the different isoforms of NFE2L1 in bone metabolism, osteoclastogenesis in particular, has not been reported. Here we demonstrate that knockout of all isoforms of Nfe2l1 transcripts specifically in the myeloid lineage in mice [Nfe2l1(M)-KO] results in increased activity of osteoclasts, decreased bone mass and worsening of osteoporosis induced by ovariectomy and aging. In comparison, LysM-Cre-mediated Nfe2l1 deletion has no significant effect on the osteoblast and osteocytes. Mechanistic investigations using bone marrow cells and RAW 264.7 cells revealed that deficiency of Nfe2l1 leads to accelerated and elevated OD, which is attributed, at least in part, to enhanced accumulation of ROS in the early stage of OD and expression of nuclear factor of activated T cells, cytoplasmic, calcineurin dependent 1α (Nfatc1/α). In addition, NFE2L1 regulates the transcription of multiple antioxidant genes and Nfatc1/α and OD in an isoform-specific manner. While long isoforms of NFE2L1 function as accelerators of induction of Nfatc1/α and antioxidant genes and OD, the short isoform NFE2L1-453 serves as a brake that keeps the long isoforms’ accelerator effects in check. These findings provide a novel insight into the regulatory roles of NFE2L1 in osteoclastogenesis and highlight that NFE2L1 is essential in regulating bone remodeling and thus may be a valuable therapeutic target for bone disorders.

Increasing evidence reveals that a broad spectrum of environmental chemicals and pharmaceutical compounds cause female ovarian toxicity (ovotoxicity). The current gold standard of ovotoxicity testing largely relies on whole laboratory animals, but in vivo models are time consuming, costly, and present animal welfare concerns. We previously demonstrated that the 3D encapsulated in vitro follicle growth (eIVFG) is a robust in vitro model for ovotoxicity testing. However, the follicle preparation process is complex and highly dependent on technical skills. Here, we aimed to use vitrification method to cryopreserve murine immature follicles for a high-content eIVFG, chemical exposure, and ovotoxicity screening. Results indicated that a closed vitrification system combined with optimized vitrification protocols preserved mouse follicle viability and functionality and vitrified follicles exhibited comparable follicle and oocyte reproductive outcomes to freshly harvested follicles during eIVFG, including follicle survival and development, ovarian steroidogenesis, and oocyte maturation and ovulation. Moreover, vitrified follicles consistently responded to ovotoxic chemical, doxorubicin (DOX). We further used vitrified follicles to test the response of microcystins (MCs), an emerging category of environmental contaminants produced by cyanobacteria associated with harmful algal blooms (HABs), and found that different congeners of MCs exhibited differential ovotoxicities. In summary, our study demonstrates that vitrification enables a long-term-storage and ready-to-use ovarian follicle bank for high-throughput ovotoxicity screening, which identifies endocrine disrupting effects of MCs.

Under oxidative and electrophilic stresses, cells launch an NRF2-mediated transcriptional antioxidant program. The activation of NRF2 depends on a redox sensor, KEAP1, which promotes the ubiquitination and degradation of NRF2. While a great deal has been learned about this duo, its quantitative signaling properties are largely unexplored. Here we examined these properties, including half-life, maximal activation, and response steepness (ultrasensitivity) of NRF2, through mathematical modeling. The models describe the binding of KEAP1 and NRF2 via ETGE and DLG motifs, NRF2 production, KEAP1-dependent and independent NRF2 degradation, and perturbations by different classes of NRF2 activators. Simulations revealed at the basal condition, NRF2 is sequestered by KEAP1 and the KEAP1-NRF2 complex is distributed comparably in an ETGE-bound (open) state and an ETGE and DLG dual-bound (closed) state. When two-step ETGE binding is considered, class I–V, electrophilic NRF2 activators shift the balance to a closed state incompetent to degrade NRF2, while the open and closed KEAP1-NRF2 complexes transition from operating in cycle mode to equilibrium mode. Ultrasensitive NRF2 activation (a steep rise of free NRF2) can occur when NRF2 nearly saturates KEAP1. The ultrasensitivity results from zero-order degradation through DLG binding and protein sequestration through ETGE binding. Optimal abundances of cytosolic and nuclear KEAP1 exist to maximize ultrasensitivity. These response characteristics do not require disruption of DLG binding as suggested by the hinge-latch hypothesis. In comparison, class VI NRF2 activators cause a shift to the open KEAP1-NRF2 complex and ultimately its complete dissociation, resulting in a fast release of NRF2 followed by stabilization. However, ultrasensitivity is lost due to decreasing free KEAP1 abundance. In summary, by simulating the dual role of KEAP1, i.e., sequestering and promoting degradation of NRF2, our modeling provides novel quantitative insights into NRF2 activation, which may help design novel NRF2 modulators and understand the oxidative actions of environmental stressors.

The KEAP1-NRF2-ARE signaling pathway plays a central role in mediating the adaptive cellular stress response to oxidative and electrophilic chemicals. This canonical pathway has been extensively studied and reviewed in the past two decades, but rarely was it looked at from a quantitative signaling perspective. Signal amplification, i.e., ultrasensitivity, is crucially important for robust induction of antioxidant genes to appropriate levels that can adequately counteract the stresses. In this review article, we examined a number of well-known molecular events in the KEAP1-NRF2-ARE pathway from a quantitative perspective with a focus on how signal amplification can be achieved. We illustrated, by using a series of mathematical models, that redox-regulated protein sequestration, stabilization, translation, nuclear trafficking, DNA promoter binding, and transcriptional induction – which are embedded in the molecular network comprising KEAP1, NRF2, sMaf, p62, and BACH1 – may generate highly ultrasensitive NRF2 activation and antioxidant gene induction. The emergence and degree of ultrasensitivity depend on the strengths of protein-protein and protein-DNA interaction and protein abundances. A unique, quantitative understanding of signal amplification in the KEAP1-NRF2-ARE pathway will help to identify sensitive targets for the prevention and therapeutics of oxidative stress-related diseases and develop quantitative adverse outcome pathway models to facilitate the health risk assessment of oxidative chemicals.

Background:Resistance to chemotherapy is a major obstacle in the treatment of human hepatocellular carcinoma (HCC). Despite playing an important role in chemoprevention, nuclear factor erythroid 2-related factor 2 (NRF2) also contributes to chemo- and radio-resistance. The current study focusses on camptothecin as a novel NRF2 inhibitor to sensitise HCC to chemotherapy.Methods:The expression and transcriptional activity of NRF2 in human HCC biopsies and camptothecin-treated culture cells were determined using immunostaining, western blot, reverse-transcription quantitative real-time PCR (RT-qPCR) and luciferase reporter assay. The effect of camptothecin on chemosensitivity of cancer cells was assessed in vitro and in xenografts.Results:The expression and transcriptional activity of NRF2 were substantially elevated in HCC biopsies compared with corresponding adjacent tissues, and positively correlated with serum α-fetoprotein, a clinical indicator of pathological progression. In searching chemicals targeting NRF2 for chemotherapy, we discovered that camptothecin is a potent NRF2 inhibitor. Camptothecin markedly suppressed NRF2 expression and transcriptional activity in different types of cancer cells including HepG2, SMMC-7721 and A549. As a result, camptothecin sensitised these cells to chemotherapeutic drugs in vitro and in xenografts.Conclusions:Camptothecin is a novel NRF2 inhibitor that may be repurposed in combination with other chemotherapeutics to enhance their efficacy in treating high NRF2-expressing cancers.