Logic of a mammalian metabolic cycle: an oscillated NAD+/NADH redox signaling regulates coordinated histone expression and S-phase progression

STAT5 Is Necessary for the Metabolic Switch Induced by IL-2 in Cervical Cancer Cell Line SiHa

The tumor cells reprogram their metabolism to cover their high bioenergetic demands for maintaining uncontrolled growth. This response can be mediated by cytokines such as IL-2, which binds to its receptor and activates the JAK/STAT pathway. Some reports show a correlation between the JAK/STAT pathway and cellular metabolism, since the constitutive activation of STAT proteins promotes glycolysis through the transcriptional activation of genes related to energetic metabolism. However, the role of STAT proteins in the metabolic switch induced by cytokines in cervical cancer remains poorly understood. In this study, we analyzed the effect of IL-2 on the metabolic switch and the role of STAT5 in this response. Our results show that IL-2 induces cervical cancer cell proliferation and the tyrosine phosphorylation of STAT5. Also, it induces an increase in lactate secretion and the ratio of NAD+/NADH, which suggest a metabolic reprogramming of their metabolism. When STAT5 was silenced, the lactate secretion and the NAD+/NADH ratio decreased. Also, the expression of HIF1α and GLUT1 decreased. These results indicate that STAT5 regulates IL-2-induced cell proliferation and the metabolic shift to aerobic glycolysis by regulating genes related to energy metabolism. Our results suggest that STAT proteins modulate the metabolic switch in cervical cancer cells to attend to their high demand of energy required for cell growth and proliferation.

Pyruvate kinase, a metabolic sensor powering glycolysis, drives the metabolic control of DNA replication

Background

In all living organisms, DNA replication is exquisitely regulated in a wide range of growth conditions to achieve timely and accurate genome duplication prior to cell division. Failures in this regulation cause DNA damage with potentially disastrous consequences for cell viability and human health, including cancer. To cope with these threats, cells tightly control replication initiation using well-known mechanisms. They also couple DNA synthesis to nutrient richness and growth rate through a poorly understood process thought to involve central carbon metabolism. One such process may involve the cross-species conserved pyruvate kinase (PykA) which catalyzes the last reaction of glycolysis. Here we have investigated the role of PykA in regulating DNA replication in the model system Bacillus subtilis.

Results

On analysing mutants of the catalytic (Cat) and C-terminal (PEPut) domains of B. subtilis PykA we found replication phenotypes in conditions where PykA is dispensable for growth. These phenotypes are independent from the effect of mutations on PykA catalytic activity and are not associated with significant changes in the metabolome. PEPut operates as a nutrient-dependent inhibitor of initiation while Cat acts as a stimulator of replication fork speed. Disruption of either PEPut or Cat replication function dramatically impacted the cell cycle and replication timing even in cells fully proficient in known replication control functions. In vitro, PykA modulates activities of enzymes essential for replication initiation and elongation via functional interactions. Additional experiments showed that PEPut regulates PykA activity and that Cat and PEPut determinants important for PykA catalytic activity regulation are also important for PykA-driven replication functions.

Conclusions

We infer from our findings that PykA typifies a new family of cross-species replication control regulators that drive the metabolic control of replication through a mechanism involving regulatory determinants of PykA catalytic activity. As disruption of PykA replication functions causes dramatic replication defects, we suggest that dysfunctions in this new family of universal replication regulators may pave the path to genetic instability and carcinogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01278-3.

p53 is regulated by aerobic glycolysis in cancer cells by the CtBP family of NADH-dependent transcriptional regulators

High rates of glycolysis in cancer cells are a well-established characteristic of many human tumors, providing rapidly proliferating cancer cells with metabolites that can be used as precursors for anabolic pathways. Maintenance of high glycolytic rates depends on the lactate dehydrogenase-catalyzed regeneration of NAD⁺ from GAPDH-generated NADH because an increased NADH:NAD⁺ ratio inhibits GAPDH. Here, using human breast cancer cell models, we identified a pathway in which changes in the extramitochondrial-free NADH:NAD⁺ ratio signaled through the CtBP family of NADH-sensitive transcriptional regulators to control the abundance and activity of p53. NADH-free forms of CtBPs cooperated with the p53-binding partner HDM2 to suppress p53 function, and loss of these forms in highly glycolytic cells resulted in p53 accumulation. We propose that this pathway represents a "glycolytic stress response" in which the initiation of a protective p53 response by an increased NADH:NAD⁺ ratio enables cells to avoid cellular damage caused by mismatches between metabolic supply and demand.

Oxidative stress-induced structural changes in the microtubule-associated flavoenzyme Irc15p from Saccharomyces cerevisiae

The genome of the yeast Saccharomyces cerevisiae encodes a canonical lipoamide dehydrogenase (Lpd1p) as part of the pyruvate dehydrogenase complex and a highly similar protein termed Irc15p (increased recombination centers 15). In contrast to Lpd1p, Irc15p lacks a pair of redox active cysteine residues required for the reduction of lipoamide and thus it is very unlikely that Irc15p performs a similar dithiol-disulfide exchange reaction as reported for lipoamide dehydrogenases. We expressed IRC15 in Escherichia coli and purified the produced protein to conduct a detailed biochemical characterization. Here, we show that Irc15p is a dimeric protein with one FAD per protomer. Photoreduction of the protein generates the fully reduced hydroquinone without the occurrence of a flavin semiquinone radical. Similarly, reduction with NADH or NADPH yields the flavin hydroquinone without the occurrence of intermediates as observed for lipoamide dehydrogenase. The redox potential of Irc15p was -313 ± 1 mV and is thus similar to lipoamide dehydrogenase. Reduced Irc15p is oxidized by several artificial electron acceptors such as potassium ferricyanide, 2,6-dichlorophenol-indophenol, 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide, and menadione. However, disulfides such as cystine, glutathione, and lipoamide were unable to react with reduced Irc15p. Limited proteolysis and SAXS-measurements revealed that the NADH-dependent formation of hydrogen peroxide caused a substantial structural change in the dimeric protein. Therefore, we hypothesize that Irc15p undergoes a conformational change in the presence of elevated levels of hydrogen peroxide, which is a putative biomarker of oxidative stress. This conformational change may in turn modulate the interaction of Irc15p with other key players involved in regulating microtubule dynamics.

Keap1-MCM3 interaction is a potential coordinator of molecular machineries of antioxidant response and genomic DNA replication in metazoa

Coordination of DNA replication and cellular redox homeostasis mechanisms is essential for the sustained genome stability due to the sensitivity of replicating DNA to oxidation. However, substantial gaps remain in our knowledge of underlying molecular pathways. In this study, we characterise the interaction of Keap1, a central antioxidant response regulator in Metazoa, with the replicative helicase subunit protein MCM3. Our analysis suggests that structural determinants of the interaction of Keap1 with its critical downstream target - Nrf2 master transactivator of oxidative stress response genes – may have evolved in evolution to mimic the conserved helix-2-insert motif of MCM3. We show that this has led to a competition between MCM3 and Nrf2 proteins for Keap1 binding, and likely recruited MCM3 for the competitive binding dependent modulation of Keap1 controlled Nrf2 activities. We hypothesise that such mechanism could help to adjust the Keap1-Nrf2 antioxidant response pathway according to the proliferative and replicative status of the cell, with possible reciprocal implications also for the regulation of cellular functions of MCM3. Altogether this suggests about important role of Keap1-MCM3 interaction in the cross-talk between replisome and redox homeostasis machineries in metazoan cells.

Spatial perspectives in the redox code-Mass spectrometric proteomics studies of moonlighting proteins

The Redox Code involves specific, reversible oxidative changes in proteins that modulate protein tertiary structure, interactions, trafficking, and activity, and hence couple the proteome to the metabolic/oxidative state of cells. It is currently a major focus of study in cell biology. Recent studies of dynamic cellular spatial reorganization with MS-based subcellular-spatial-razor proteomics reveal that protein constituents of many subcellular structures, including mitochondria, the endoplasmic reticulum, the plasma membrane, and the extracellular matrix, undergo changes in their subcellular abundance/distribution in response to oxidative stress. These proteins are components of a diverse variety of functional processes spatially distributed across cells. Many of the same proteins are involved in response to suppression of DNA replication indicate that oxidative stress is strongly intertwined with DNA replication/proliferation. Both are replete with networks of moonlighting proteins that show coordinated changes in subcellular location and that include primary protein actuators of the redox code involved in the processing of NAD⁺ /NADH, NADP⁺ /NADPH, Cys/CySS, and GSH/GSSG redox couples. Small groups of key proteins such as {KPNA2, KPNB1, PCNA, PTMA, SET} constitute "spatial switches" that modulate many nuclear processes. Much of the functional response involves subcellular protein trafficking, including nuclear import/export processes, vesicle-mediated trafficking, the endoplasmic reticulum/Golgi pathway, chaperone-assisted processes, and other transport systems. This is not visible to measurements of total protein abundance by transcriptomics or proteomics. Comprehensive pictures of cellular function will require collection of data on the subcellular transport and local functions of many moonlighting proteins, especially of those with critical roles in spatial coordination across cells. The proteome-wide analysis of coordinated changes in abundance and trafficking of proteins offered by MS-based proteomics has a unique, crucial role to play in deciphering the complex adaptive systems that underlie cellular function. © 2016 Wiley Periodicals, Inc. Mass Spec Rev.

High-content image analysis (HCIA) assay has the highest correlation with direct counting cell suspension compared to the ATP, WST-8 and Alamar blue assays for measurement of cytotoxicity

INTRODUCTION

Various cytotoxicity assays measuring indicators such as enzyme activity, dye uptake, or cellular ATP content are often performed using 96-well microplates. However, recent reports show that cytotoxicity assays such as the ATP assay and MTS assay underestimate cytotoxicity when compounds such as anti-cancer drugs or mutagens induce cell hypertrophy whilst increasing intracellular ATP content. Therefore, we attempted to evaluate the reliability of a high-content image analysis (HCIA) assay to count cell number in a 96-well microplate automatically without using a cell-number indicator.

METHODS

We compared cytotoxicity results of 25 compounds obtained from ATP, WST-8, Alamar blue, and HCIA assays with those directly measured using an automatic cell counter, and repeating individual experiments thrice.

RESULTS

The number of compounds showing low correlation in cell viability measured using cytotoxicity assays compared to automatic cell counting (r²<0.8, at least 2 of 3 experiments) were follows: ATP assay; 7; WST-8 assay, 2; Alamar blue assay, 3; HCIA cytotoxicity assay, 0. Compounds for which correlation was poor in 3 assays, except the HCIA assay, induced an increase in nuclear and cell size. However, correlation between cell viability measured by automatic cell counter and the HCIA assay was strong regardless of nuclear and cell size. Additionally, correlation coefficients between IC₅₀ values obtained from automatic cell counter and from cytotoxicity assays were as follows: ATP assay, 0.80; WST-8 assay, 0.84; Alamar blue assay, 0.84; and HCIA assay, 0.98.

DISCUSSION

From the above, we showed that the HCIA cytotoxicity assay produces similar data to the automatic cell counter and is highly accurate in measuring cytotoxicity.

The inhibition of lactate dehydrogenase A hinders the transcription of histone 2B gene independently from the block of aerobic glycolysis

Most cancer cells use aerobic glycolysis to fuel their growth and many efforts are made to selectively block this metabolic pathway in cancer cells by inhibiting lactate dehydrogenase A (LDHA). However, LDHA is a moonlighting protein which exerts functions also in the nucleus as a factor associated to transcriptional complexes. Here we found that two small molecules which inhibit the enzymatic activity of LDHA hinder the transcription of histone 2B gene independently from the block of aerobic glycolysis. Moreover, we observed that silencing this gene reduces cell replication, hence suggesting that the inhibition of LDHA can also affect the proliferation of normal non-glycolysing dividing cells.

The Redox Status of Cancer Cells Supports Mechanisms behind the Warburg Effect

To better understand the energetic status of proliferating cells, we have measured the intracellular pH (pHi) and concentrations of key metabolites, such as adenosine triphosphate (ATP), nicotinamide adenine dinucleotide (NAD), and nicotinamide adenine dinucleotide phosphate (NADP) in normal and cancer cells, extracted from fresh human colon tissues. Cells were sorted by elutriation and segregated in different phases of the cell cycle (G0/G1/S/G2/M) in order to study their redox (NAD, NADP) and bioenergetic (ATP, pHi) status. Our results show that the average ATP concentration over the cell cycle is higher and the pHi is globally more acidic in normal proliferating cells. The NAD⁺/NADH and NADP⁺/NADPH redox ratios are, respectively, five times and ten times higher in cancer cells compared to the normal cell population. These energetic differences in normal and cancer cells may explain the well-described mechanisms behind the Warburg effect. Oscillations in ATP concentration, pHi, NAD⁺/NADH, and NADP⁺/NADPH ratios over one cell cycle are reported and the hypothesis addressed. We also investigated the mitochondrial membrane potential (MMP) of human and mice normal and cancer cell lines. A drastic decrease of the MMP is reported in cancer cell lines compared to their normal counterparts. Altogether, these results strongly support the high throughput aerobic glycolysis, or Warburg effect, observed in cancer cells.

ABCE1 is essential for S phase progression in human cells

ABCE1 is a highly conserved protein universally present in eukaryotes and archaea, which is crucial for the viability of different organisms. First identified as RNase L inhibitor, ABCE1 is currently recognized as an essential translation factor involved in several stages of eukaryotic translation and ribosome biogenesis. The nature of vital functions of ABCE1, however, remains unexplained. Here, we study the role of ABCE1 in human cell proliferation and its possible connection to translation. We show that ABCE1 depletion by siRNA results in a decreased rate of cell growth due to accumulation of cells in S phase, which is accompanied by inefficient DNA synthesis and reduced histone mRNA and protein levels. We infer that in addition to the role in general translation, ABCE1 is involved in histone biosynthesis and DNA replication and therefore is essential for normal S phase progression. In addition, we analyze whether ABCE1 is implicated in transcript-specific translation via its association with the eIF3 complex subunits known to control the synthesis of cell proliferation-related proteins. The expression levels of a few such targets regulated by eIF3A, however, were not consistently affected by ABCE1 depletion.

Cell cycle progression is regulated by intertwined redox oscillators

The different phases of the eukaryotic cell cycle are exceptionally well-preserved phenomena. DNA decompaction, RNA and protein synthesis (in late G1 phase) followed by DNA replication (in S phase) and lipid synthesis (in G2 phase) occur after resting cells (in G0) are committed to proliferate. The G1 phase of the cell cycle is characterized by an increase in the glycolytic metabolism, sustained by high NAD+/NADH ratio. A transient cytosolic acidification occurs, probably due to lactic acid synthesis or ATP hydrolysis, followed by cytosolic alkalinization. A hyperpolarized transmembrane potential is also observed, as result of sodium/potassium pump (NaK-ATPase) activity. During progression of the cell cycle, the Pentose Phosphate Pathway (PPP) is activated by increased NADP+/NADPH ratio, converting glucose 6-phosphate to nucleotide precursors. Then, nucleic acid synthesis and DNA replication occur in S phase. Along with S phase, unpublished results show a cytosolic acidification, probably the result of glutaminolysis occurring during this phase. In G2 phase there is a decrease in NADPH concentration (used for membrane lipid synthesis) and a cytoplasmic alkalinization occurs. Mitochondria hyperfusion matches the cytosolic acidification at late G1/S transition and then triggers ATP synthesis by oxidative phosphorylation. We hypothesize here that the cytosolic pH may coordinate mitochondrial activity and thus the different redox cycles, which in turn control the cell metabolism.

Glyceraldehyde 3-phosphate dehydrogenase-telomere association correlates with redox status in Trypanosoma cruzi

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is a classical metabolic enzyme involved in energy production and plays a role in additional nuclear functions, including transcriptional control, recognition of misincorporated nucleotides in DNA and maintenance of telomere structure. Here, we show that the recombinant protein T. cruzi GAPDH (rTcGAPDH) binds single-stranded telomeric DNA. We demonstrate that the binding of GAPDH to telomeric DNA correlates with the balance between oxidized and reduced forms of nicotinamide adenine dinucleotides (NAD⁺/NADH). We observed that GAPDH-telomere association and NAD⁺/NADH balance changed throughout the T. cruzi life cycle. For example, in replicative epimastigote forms of T. cruzi, which show similar intracellular concentrations of NAD⁺ and NADH, GAPDH binds to telomeric DNA in vivo and this binding activity is inhibited by exogenous NAD⁺. In contrast, in the T. cruzi non-proliferative trypomastigote forms, which show higher NAD⁺ concentration, GAPDH was absent from telomeres. In addition, NAD⁺ abolishes physical interaction between recombinant GAPDH and synthetic telomere oligonucleotide in a cell free system, mimicking exogenous NAD⁺ that reduces GAPDH-telomere interaction in vivo. We propose that the balance in the NAD⁺/NADH ratio during T. cruzi life cycle homeostatically regulates GAPDH telomere association, suggesting that in trypanosomes redox status locally modulates GAPDH association with telomeric DNA.

Integration of the metabolic/redox state, histone gene switching, DNA replication and S-phase progression by moonlighting metabolic enzymes

The concept of one-protein–multiple-function, i.e. moonlighting proteins, is an ever-expanding paradigm. We obtained compelling evidence that an array of ‘cytoplasmic’ metabolic enzymes can enter the nuclei to carry out moonlighting transcription functions; this phenomenon is conserved from Drosophila to humans. Of particular interest are the classical glycolytic enzymes GAPDH (glyceraldehyde-3-phosphate dehydrogenase) and LDH (lactate dehydrogenase), which utilize NAD(H) as coenzymes and not only moonlight (in their nuclear forms) to regulate the transcription of S-phase-specific histone genes, but also act as metabolic/redox sensors that link histone gene switching to DNA replication and S-phase progression.

CBP/p300 and SIRT1 are involved in transcriptional regulation of S-phase specific histone genes

BACKGROUND

Histones constitute a type of essential nuclear proteins important for chromatin structure and functions. The expression of major histones is strictly confined to the S phase of a cell cycle and tightly coupled to DNA replication.

METHODOLOGY/PRINCIPAL FINDINGS

With RT-qPCR and ChIP assays, we investigated transcriptional regulation of the S-phase specific histone genes and found that the acetylation level of histones on core histone gene promoters fluctuated during cell cycle in a pattern similar to RNA polymerase II association. Further, we showed that CBP/p300 and SIRT1 were recruited to histone gene promoters in an NPAT-dependent manner, knockdown of which affected histone acetylation on histone gene promoters and histone gene transcription.

SIGNIFICANCE

These observations contribute to further understanding of the mechanism by which the expression of canonical histone genes is regulated, and also implicate a link between histone expression and DNA damage repair and cell metabolism.

Conserved catalytic and C-terminal regulatory domains of the C-terminal binding protein corepressor fine-tune the transcriptional response in development

Transcriptional corepressors play complex roles in developmental gene regulation. These proteins control transcription by recruiting diverse chromatin-modifying enzymes, but it is not known whether corepressor activities are finely regulated in different developmental settings or whether their basic activities are identical in most contexts. The evolutionarily conserved C-terminal binding protein (CtBP) is recruited by a variety of transcription factors that play crucial roles in development and disease. CtBP contains a central NAD(H) binding core domain that is homologous to D2 hydroxy acid dehydrogenase enzymes, as well as an unstructured C-terminal domain. NAD(H) binding is important for CtBP function, but the significance of its intrinsic dehydrogenase activity, as well as that of the unstructured C terminus, is poorly understood. To clarify the biological relevance of these features, we established genetic rescue assays to determine how different forms of CtBP function in the context of Drosophila melanogaster development. The mutant phenotypes and specific gene regulatory effects indicate that both the catalytic site of CtBP and the C-terminal extension play important, if nonessential roles in development. Our results indicate that the structural and enzymatic features of CtBP, previously thought to be dispensable for overall transcriptional control, are critical for modulating this protein's activity in diverse developmental settings.

Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes

Neutral nucleotide substitutions occur at varying rates along genomes, and it remains a major issue to unravel the mechanisms that cause these variations and to analyze their evolutionary consequences. Here, we study the role of replication in the neutral substitution pattern. We obtained a high-resolution replication timing profile of the whole human genome by massively parallel sequencing of nascent BrdU-labeled replicating DNA. These data were compared to the neutral substitution rates along the human genome, obtained by aligning human and chimpanzee genomes using macaque and orangutan as outgroups. All substitution rates increase monotonously with replication timing even after controlling for local or regional nucleotide composition, crossover rate, distance to telomeres, and chromatin compaction. The increase in non-CpG substitution rates might result from several mechanisms including the increase in mutation-prone activities or the decrease in efficiency of DNA repair during the S phase. In contrast, the rate of C --> T transitions in CpG dinucleotides increases in later-replicating regions due to increasing DNA methylation level that reflects a negative correlation between timing and gene expression. Similar results are observed in the mouse, which indicates that replication timing is a main factor affecting nucleotide substitution dynamics at non-CpG sites and constitutes a major neutral process driving mammalian genome evolution.

Drosophila octamer elements and Pdm-1 dictate the coordinated transcription of core histone genes

We reveal a set of divergent octamer elements in Drosophila melanogaster (dm) core histone gene promoters. These elements recruit transcription factor POU-domain protein in D. melanogaster 1 (Pdm-1), which along with co-activator dmOct-1 coactivator in S-phase (dmOCA-S), activates transcription from at least the Drosophila histone 2B (dmH2B) and 4 (dmH4) promoters in a fashion similar to the transcription of mammalian histone 2B (H2B) gene activated by octamer binding transcription factor 1 (Oct-1) and Oct-1 coactivator in S-phase (OCA-S). The expression of core histone genes in both kingdoms is coordinated; however, although the expression of mammalian histone genes involves subtype-specific transcription factors and/or co-activator(s), the expression of Drosophila core histone genes is regulated by a common module (Pdm-1/dmOCA-S) in a directly coordinated manner. Finally, dmOCA-S is recruited to the Drosophila histone locus bodies in the S-phase, marking S-phase-specific transcription activation of core histone genes.