More Than a Metabolic Enzyme: MTHFD2 as a Novel Target for Anticancer Therapy?

MTHFD2-mediated redox homeostasis promotes gastric cancer progression under hypoxic conditions

OBJECTIVES

Cancer cells undergo metabolic reprogramming to adapt to high oxidative stress, but little is known about how metabolic remodeling enables gastric cancer cells to survive stress associated with aberrant reactive oxygen species (ROS) production. Here, we aimed to identify the key metabolic enzymes that protect gastric cancer (GC) cells from oxidative stress.

METHODS

ROS level was detected by DCFH-DA probes. Multiple cell biological studies were performed to identify the underlying mechanisms. Furthermore, cell-based xenograft and patient-derived xenograft (PDX) model were performed to evaluate the role of MTHFD2 in vivo.

RESULTS

We found that overexpression of MTHFD2, but not MTHFD1, is associated with reduced overall and disease-free survival in gastric cancer. In addition, MTHFD2 knockdown reduces the cellular NADPH/NADP+ ratio, colony formation and mitochondrial function, increases cellular ROS and cleaved PARP levels and induces in cell death under hypoxia, a hallmark of solid cancers and a common inducer of oxidative stress. Moreover, genetic or pharmacological inhibition of MTHFD2 reduces tumor burden in both tumor cell lines and patient-derived xenograft-based models.

DISCUSSION

our study highlights the crucial role of MTHFD2 in redox regulation and tumor progression, demonstrating the therapeutic potential of targeting MTHFD2.

Loss of the histone chaperone UNC-85/ASF1 inhibits the epigenome-mediated longevity and modulates the activity of one-carbon metabolism

Histone H3/H4 chaperone anti-silencing function 1 (ASF1) is a conserved factor mediating nucleosomal assembly and disassembly, playing crucial roles in processes such as replication, transcription, and DNA repair. Nevertheless, its involvement in aging has remained unclear. Here, we utilized the model organism Caenorhabditis elegans to demonstrate that the loss of UNC-85, the homolog of ASF1, leads to a shortened lifespan in a multicellular organism. Furthermore, we show that UNC-85 is required for epigenome-mediated longevity, as knockdown of the histone H3 lysine K4 methyltransferase ash-2 does not extend the lifespan of unc-85 mutants. In this context, we found that the longevity-promoting ash-2 RNA interference enhances UNC-85 activity by increasing its nuclear localization. Finally, our data indicate that the loss of UNC-85 increases the activity of one-carbon metabolism, and that downregulation of the one-carbon metabolism component dao-3/MTHFD2 partially rescues the short lifespan of unc-85 mutants. Together, these findings reveal UNC-85/ASF1 as a modulator of the central metabolic pathway and a factor regulating a pro-longevity response, thus shedding light on a mechanism of how nucleosomal maintenance associates with aging.

Immunometabolism in cancer: basic mechanisms and new targeting strategy

Maturing immunometabolic research empowers immune regulation novel approaches. Progressive metabolic adaptation of tumor cells permits a thriving tumor microenvironment (TME) in which immune cells always lose the initial killing capacity, which remains an unsolved dilemma even with the development of immune checkpoint therapies. In recent years, many studies on tumor immunometabolism have been reported. The development of immunometabolism may facilitate anti-tumor immunotherapy from the recurrent crosstalk between metabolism and immunity. Here, we discuss clinical studies of the core signaling pathways of immunometabolism and their inhibitors or agonists, as well as the specific functions of these pathways in regulating immunity and metabolism, and discuss some of the identified immunometabolic checkpoints. Understanding the comprehensive advances in immunometabolism helps to revise the status quo of cancer treatment. An overview of the new landscape of immunometabolism. The PI3K pathway promotes anabolism and inhibits catabolism. The LKB1 pathway inhibits anabolism and promotes catabolism. Overactivation of PI3K/AKT/mTOR pathway and IDO, IL4I1, ACAT, Sirt2, and MTHFD2 promote immunosuppression of TME formation, as evidenced by increased Treg and decreased T-cell proliferation. The LKBI-AMPK pathway promotes the differentiation of naive T cells to effector T cells and memory T cells and promotes anti-tumor immunity in DCs.

Radiological and Molecular Analysis of Radioiodinated Anastrozole and Epirubicin as Innovative Radiopharmaceuticals Targeting Methylenetetrahydrofolate Dehydrogenase 2 in Solid Tumors

In the dynamic field of radiopharmaceuticals, innovating targeted agents for cancer diagnosis and therapy is crucial. Our study enriches this evolving landscape by evaluating the potential of radioiodinated anastrozole ([125I]anastrozole) and radioiodinated epirubicin ([125I]epirubicin) as targeting agents against MTHFD2-driven tumors. MTHFD2, which is pivotal in one-carbon metabolism, is notably upregulated in various cancers, presenting a novel target for radiopharmaceutical application. Through molecular docking and 200 ns molecular dynamics (MD) simulations, we assess the binding efficiency and stability of [125I]anastrozole and [125I]epirubicin with MTHFD2. Molecular docking illustrates that [125I]epirubicin has a superior binding free energy (∆Gbind) of -41.25 kJ/mol compared to -39.07 kJ/mol for [125I]anastrozole and -38.53 kJ/mol for the control ligand, suggesting that it has a higher affinity for MTHFD2. MD simulations reinforce this, showing stable binding, as evidenced by root mean square deviation (RMSD) values within a narrow range, underscoring the structural integrity of the enzyme-ligand complexes. The root mean square fluctuation (RMSF) analysis indicates consistent dynamic behavior of the MTHFD2 complex upon binding with [125I]anastrozole and [125I]epirubicin akin to the control. The radius of gyration (RG) measurements of 16.90 Å for MTHFD2-[125I]anastrozole and 16.84 Å for MTHFD2-[125I]epirubicin confirm minimal structural disruption upon binding. The hydrogen bond analysis reveals averages of two and three stable hydrogen bonds for [125I]anastrozole and [125I]epirubicin complexes, respectively, highlighting crucial stabilizing interactions. The MM-PBSA calculations further endorse the thermodynamic favorability of these interactions, with binding free energies of -48.49 ± 0.11 kJ/mol for [125I]anastrozole and -43.8 kJ/mol for MTHFD2-. The significant contribution of Van der Waals and electrostatic interactions to the binding affinities of [125I]anastrozole and [125I]epirubicin, respectively, underscores their potential efficacy for targeted tumor imaging and therapy. These computational findings lay the groundwork for the future experimental validation of [125I]anastrozole and [125I]epirubicin as MTHFD2 inhibitors, heralding a notable advancement in precision oncology tools. The data necessitate subsequent in vitro and in vivo assays to corroborate these results.

Revolutionizing breast cancer treatment: Harnessing the related mechanisms and drugs for regulated cell death (Review)

Breast cancer arises from the malignant transformation of mammary epithelial cells under the influence of various carcinogenic factors, leading to a gradual increase in its prevalence. This disease has become the leading cause of mortality among female malignancies, posing a significant threat to the health of women. The timely identification of breast cancer remains challenging, often resulting in diagnosis at the advanced stages of the disease. Conventional therapeutic approaches, such as surgical excision, chemotherapy and radiotherapy, exhibit limited efficacy in controlling the progression and metastasis of the disease. Regulated cell death (RCD), a process essential for physiological tissue cell renewal, occurs within the body independently of external influences. In the context of cancer, research on RCD primarily focuses on cuproptosis, ferroptosis and pyroptosis. Mounting evidence suggests a marked association between these specific forms of RCD, and the onset and progression of breast cancer. For example, a cuproptosis vector can effectively bind copper ions to induce cuproptosis in breast cancer cells, thereby hindering their proliferation. Additionally, the expression of ferroptosis‑related genes can enhance the sensitivity of breast cancer cells to chemotherapy. Likewise, pyroptosis‑related proteins not only participate in pyroptosis, but also regulate the tumor microenvironment, ultimately leading to the death of breast cancer cells. The present review discusses the unique regulatory mechanisms of cuproptosis, ferroptosis and pyroptosis in breast cancer, and the mechanisms through which they are affected by conventional cancer drugs. Furthermore, it provides a comprehensive overview of the significance of these forms of RCD in modulating the efficacy of chemotherapy and highlights their shared characteristics. This knowledge may provide novel avenues for both clinical interventions and fundamental research in the context of breast cancer.

Construction of an oxidative phosphorylation-related gene signature for predicting prognosis and identifying immune infiltration in osteosarcoma

BACKGROUND

Osteosarcoma is a prevalent malignant tumor that originates from mesenchymal tissue. It typically affects children and adolescents. Although it is known that the growth of osteosarcoma relies on oxidative phosphorylation for energy production, limited attention has been paid to exploring the potential of oxidative phosphorylation-related genes in predicting the prognosis of individuals suffering from osteosarcoma.

METHODS

All the data were retrieved from the UCSC Xena and GEO (GENE EXPRESSION OMNIBUS). Identification of the oxidative phosphorylation genes linked to the prognosis of individuals with osteosarcoma was done by means of univariate COX and LASSO regression analyses. Following that, patients were categorized into a high-risk group and a low-risk group as per the risk score determined by the identified oxidative phosphorylation genes. Furthermore, a comparison was made in terms of the survival and immune infiltration between both groups, and the prognostic model was established.

RESULTS

Five oxidative phosphorylation genes (ATP6V0D1, LHPP, COX6A2, MTHFD2, NDUFB9) associated with the prognosis of individuals with osteosarcoma were identified and the risk prognostic models were constructed. In the current research, the analysis of the ROC curves indicated a superior predictive accuracy exhibited by the risk model. The prognosis was adversely affected by immune infiltration in the high-risk group in comparison with the low-risk group. The function of the oxidative phosphorylation-related prognostic gene set was verified by GO and KEGG analysis. Furthermore, the link between oxidative phosphorylation-related genes and osteosarcoma immune infiltration was examined by GSEA analysis.

CONCLUSIONS

In this study, a prognostic model that demonstrated good predictive performance was constructed. Additionally, this study highlighted a correlation between oxidative phosphorylation-related genes and immune infiltration.

Metabolic subtypes and immune landscapes in esophageal squamous cell carcinoma: prognostic implications and potential for personalized therapies

BACKGROUND

This study aimed to identify metabolic subtypes in ESCA, explore their relationship with immune landscapes, and establish a metabolic index for accurate prognosis assessment.

METHODS

Clinical, SNP, and RNA-seq data were collected from 80 ESCA patients from the TCGA database and RNA-seq data from the GSE19417 dataset. Metabolic genes associated with overall survival (OS) and progression-free survival (PFS) were selected, and k-means clustering was performed. Immune-related pathways, immune infiltration, and response to immunotherapy were predicted using bioinformatic algorithms. Weighted gene co-expression network analysis (WGCNA) was conducted to identify metabolic genes associated with co-expression modules. Lastly, cell culture and functional analysis were performed using patient tissue samples and ESCA cell lines to verify the identified genes and their roles.

RESULTS

Molecular subtypes were identified based on the expression profiles of metabolic genes, and univariate survival analysis revealed 163 metabolic genes associated with ESCA prognosis. Consensus clustering analysis classified ESCA samples into three distinct subtypes, with MC1 showing the poorest prognosis and MC3 having the best prognosis. The subtypes also exhibited significant differences in immune cell infiltration, with MC3 showing the highest scores. Additionally, the MC3 subtype demonstrated the poorest response to immunotherapy, while the MC1 subtype was the most sensitive. WGCNA analysis identified gene modules associated with the metabolic index, with SLC5A1, NT5DC4, and MTHFD2 emerging as prognostic markers. Gene and protein expression analysis validated the upregulation of MTHFD2 in ESCA. MTHFD2 promotes the progression of ESCA and may be a potential therapeutic target for ESCA.

CONCLUSION

The established metabolic index and identified metabolic genes offer potential for prognostic assessment and personalized therapeutic interventions for ESCA, underscoring the importance of targeting metabolism-immune interactions in ESCA. MTHFD2 promotes the progression of ESCA and may be a potential therapeutic target for ESCA.

A Randomized Phase III Study of Arfolitixorin versus Leucovorin with 5-Fluorouracil, Oxaliplatin, and Bevacizumab for First-Line Treatment of Metastatic Colorectal Cancer: The AGENT Trial

PURPOSE

Suboptimal treatment outcomes with 5-fluorouracil (5-FU)/folate, the standard of care for metastatic colorectal cancer (mCRC), have generated interest in optimizing the folate. Arfolitixorin ([6R]-5,10-methylene-tetrahydrofolate) is an immediately active folate and may improve outcomes over the existing standard of care (leucovorin).

EXPERIMENTAL DESIGN

AGENT was a randomized, phase III study (NCT03750786). Patients with mCRC were randomized to arfolitixorin (120 mg/m2 given as two intravenous bolus doses of 60 mg/m2) or leucovorin (400 mg/m2 given as a single intravenous infusion) plus 5-FU, oxaliplatin, and bevacizumab. Assessments were performed every 8 weeks. The primary endpoint was the superiority of arfolitixorin for overall response rate (ORR).

RESULTS

Between February 2019 and April 2021, 490 patients were randomized (245 to each arm). After a median follow-up of 266 days, the primary endpoint of superiority for ORR was not achieved (48.2% for arfolitixorin vs. 49.4% for leucovorin, Psuperiority = 0.57). Outcomes were not achieved for median progression-free survival (PFS; 12.8 and 11.6 months, P = 0.38), median duration of response (12.2 and 12.9 months, P = 0.40), and median overall survival (23.8 and 28.0 months, P = 0.78). The proportion of patients with an adverse event of grade ≥3 severity was similar between arms (68.7% and 67.2%, respectively), as was quality of life. BRAF mutations and MTHFD2 expression were both associated with a lower PFS with arfolitixorin.

CONCLUSIONS

The study failed to demonstrate clinical benefit of arfolitixorin (120 mg/m2) over leucovorin. However, it provides some useful insights from the first-line treatment setting, including the effect of gene expression on outcomes.

SIGNIFICANCE

This phase III study compared arfolitixorin, a direct-acting folate, with leucovorin in FOLFOX plus bevacizumab in mCRC. Arfolitixorin (120 mg/m2) did not improve the ORR, potentially indicating a suboptimal dose.

Supramolecular self-assembled peptide-engineered nanofibers: A propitious proposition for cancer therapy

Cancer is a devastating disease that causes a substantial number of deaths worldwide. Current therapeutic interventions for cancer include chemotherapy, radiation therapy, or surgery. These conventional therapeutic approaches are associated with disadvantages such as multidrug resistance, destruction of healthy tissues, and tissue toxicity. Therefore, there is a paradigm shift in cancer management wherein nanomedicine-based novel therapeutic interventions are being explored to overcome the aforementioned disadvantages. Supramolecular self-assembled peptide nanofibers are emerging drug delivery vehicles that have gained much attention in cancer management owing to their biocompatibility, biodegradability, biomimetic property, stimuli-responsiveness, transformability, and inherent therapeutic property. Supramolecules form well-organized structures via non-covalent linkages, the intricate molecular arrangement helps to improve tissue permeation, pharmacokinetic profile and chemical stability of therapeutic agents while enabling targeted delivery and allowing efficient tumor imaging. In this review, we present fundamental aspects of peptide-based self-assembled nanofiber fabrication their applications in monotherapy/combinatorial chemo- and/or immuno-therapy to overcome multi-drug resistance. The role of self-assembled structures in targeted/stimuli-responsive (pH, enzyme and photo-responsive) drug delivery has been discussed along with the case studies. Further, recent advancements in peptide nanofibers in cancer diagnosis, imaging, gene therapy, and immune therapy along with regulatory obstacles towards clinical translation have been deliberated.

Overexpression of MTHFD2 represents an inflamed tumor microenvironment and precisely predicts the molecular subtype and immunotherapy response of bladder cancer

Introduction

Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2), whose aberrant expression is common in cancers, has recently been identified as a potential regulator of immune response. However, its immune-related role in bladder cancer (BLCA) and its association with immunotherapy efficacy remain unclear.

Methods

RNA sequencing data from The Cancer Genome Atlas (TCGA) was applied to analyze the immunological roles and prognostic value of MTHFD2 in pan-cancers. The association of MTHFD2 with several immunological features of tumor microenvironment (TME), including cancer-immunity cycle, immune cells infiltration, immune checkpoints expression, and T cell inflamed score was analyzed in TCGA-BLCA cohort. The predictors of cancer treatments effectiveness, including the expression and mutation of certain genes, molecular subtypes, and several signatures were evaluated as well. These results were validated by another independent cohort (GSE48075). Finally, the predictive value of MTHFD2 for TME and immunotherapy efficacy were validated using immunohistochemistry assay and RNA sequencing data from IMvigor210 cohort, respectively.

Results

MTHFD2 was found to be positively associated with several immunological features of an inflamed tumor microenvironment (TME) in various cancers and could predict BLCA patients' prognosis. In BLCA, high expression of MTHFD2 was observed to be positively related with the cancer-immunity cycle, the infiltration of several immune cells, and the expression of immunoregulators and T-cell inflamed scores, indicating a positive correlation with the inflamed TME. Moreover, patients with high MTHFD2 expression were more likely to be basal-like subtypes and respond to BLCA treatments, including immunotherapy, chemotherapy, and target therapy. The clinical data of the IMvigor210 cohort confirmed the higher response rates and better survival benefits of immunotherapy in high-MTHFD2-expression patients.

Conclusion

Collectively, high MTHFD2 predicts an inflamed TME, a basal-like subtype, and a better response to various therapeutic strategies, especially the ICB therapy, in bladder cancer.

Mitochondrial subtype MB-G3 contains potential novel biomarkers and therapeutic targets associated with prognosis of medulloblastoma

BACKGROUND

Medulloblastoma is the most common malignant brain tumor in children. There are four groups, each with different causal mutations, affected pathways and prognosis. Here, we investigated the role of mitochondria in medulloblastoma and whether there are differences between the different groups.

METHODS

We compared the gene expression levels in the four different medulloblastoma groups (MB-WNT, MB-SHH, MB-G3 and MB-G4), with the focus on genes associated with mitochondria. We used several tools including Salmon, Tximeta, DESeq2, BiomaRt, STRING, Ggplot2, EnhancedVolcano, Venny 2.1 and Metscape.

RESULTS

A total of 668 genes were differentially expressed and the most abundant genes were associated with cell division pathway followed by modulation of chemical synaptic transmission. We also identified several genes (ABAT, SOX9, ALDH5A, FOXM1, ABL1, NHLH1, NEUROD1 and NEUROD2) known to play vital role in medulloblastoma. Comparative expression analysis revealed OXPHOS complex-associated proteins of mitochondria. The most significantly expressed genes in the MB-SHH and MB-G4 groups were AHCYL1 and SFXN5 while PAICS was significantly upregulated in MB-WNT group. Notably, MB-G3 contained the most downregulated genes from the OXPHOS complexes, except COX6B2 which was strongly upregulated.

CONCLUSIONS

We show the importance of mitochondria and compare their role in the four different medulloblastoma groups.

Mitohormesis

Perturbation of mitochondrial function can trigger a host of cellular responses that seek to restore cellular metabolism, cytosolic proteostasis, and redox homeostasis. In some cases, these responses persist even after the stress is relieved, leaving the cell or tissue in a less vulnerable state. This process-termed mitohormesis-is increasingly viewed as an important aspect of normal physiology and a critical modulator of various disease processes. Here, we review aspects of mitochondrial stress signaling that, among other things, can rewire the cell's metabolism, activate the integrated stress response, and alter cytosolic quality-control pathways. We also discuss how these pathways are implicated in various disease states from pathogen challenge to chemotherapeutic resistance and how their therapeutic manipulation can lead to new strategies for a host of chronic conditions including aging itself.

MTHFD2 ablation in T cells protects against heart transplant rejection by perturbing IRF4/PD-1 pathway through the metabolic-epigenetic nexus

BACKGROUND

One-carbon metabolism supports the activation, proliferation, and function of multiple immune cells. However, researchers have not clearly determined whether and how one-carbon metabolic enzymes contribute to heart transplant rejection.

METHODS

We investigated the dynamic metabolic adaptation in grafts during heart transplant rejection by conducting transcriptomics, metabolomics and single-cell RNA sequencing studies of cardiac tissue from human and mouse heart transplant recipients. We also assessed the expression of the one-carbon metabolic enzyme methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) in cardiac grafts by immunofluorescence and flow cytometry assays. Then we constructed a murine heart transplant model with T cell-specific Mthfd2 knockout mice, analyzed T cells function by flow cytometry assays and enzyme-linked immunospot assays, and studied the mechanism by Cleavage Under Targets and Tagmentation assays. Finally, we studied the effect of a pharmacological inhibitor of MTHFD2 in humanized skin transplant model.

RESULTS

We revealed that the one-carbon metabolism enzyme MTHFD2 was a hallmark of alloreactive T cells and was linked to T cell proliferation and function after exposure to alloantigen. And, Mthfd2 ablation prevented murine heart transplant rejection. Mechanistically, we found Mthfd2 ablation affected the interferon regulatory factor 4/programmed death-1 pathway through a metabolic-epigenetic mechanism involving H3K4me3. Furthermore, we found that inhibiting MTHFD2 attenuated human allograft rejection in a humanized skin transplant model.

CONCLUSIONS

These data show that the one-carbon metabolic enzyme MTHFD2 serves as a metabolic checkpoint of alloreactive T cells and suggest that it may be a potential therapeutic target for heart transplant rejection.

Targeting allosteric binding site in methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) to identify natural product inhibitors via structure-based computational approach

Cancer has been viewed as one of the deadliest diseases worldwide. Among various types of cancer, breast cancer is the most common type of cancer in women. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is a promising druggable target and is overexpressed in cancerous cells, like, breast cancer. We conducted structure-based modeling on the allosteric site of the enzyme. Targeting the allosteric site avoids the problem of drug resistance. Pharmacophore modeling, molecular docking, HYDE assessment, drug-likeness, ADMET predictions, simulations, and free-energy calculations were performed. The RMSD, RMSF, RoG, SASA, and Hydrogen-bonding studies showed that seven candidates displayed stable behaviour. As per the literature, average superimposed simulated structures revealed a similar protein conformational change in the αE'-βf' loop, causing its displacement away from the allosteric site. The MM-PBSA showed tight binding of six compounds with the allosteric pocket. The effect of inhibitors interacting in the allosteric site causes a decrease in the binding energy of J49 (active-site inhibitor), suggesting the effect of allosteric binding. The PCA and FEL analysis revealed the significance of the docked compounds in the stable behaviour of the complexes. The outcome can contribute to the development of potential natural products with drug-like properties that can inhibit the MTHFD2 enzyme.

Constructing a prognostic model for head and neck squamous cell carcinoma based on glucose metabolism related genes

Background

Glucose metabolism (GM) plays a crucial role in cancer cell proliferation, tumor growth, and survival. However, the identification of glucose metabolism-related genes (GMRGs) for effective prediction of prognosis in head and neck squamous cell carcinoma (HNSC) is still lacking.

Methods

We conducted differential analysis between HNSC and Normal groups to identify differentially expressed genes (DEGs). Key module genes were obtained using weighted gene co-expression network analysis (WGCNA). Intersection analysis of DEGs, GMRGs, and key module genes identified GMRG-DEGs. Univariate and multivariate Cox regression analyses were performed to screen prognostic-associated genes. Independent prognostic analysis of clinical traits and risk scores was implemented using Cox regression. Gene set enrichment analysis (GSEA) was used to explore functional pathways and genes between high- and low-risk groups. Immune infiltration analysis compared immune cells between the two groups in HNSC samples. Drug prediction was performed using the Genomics of Drug Sensitivity in Cancer (GDSC) database. Quantitative real-time fluorescence PCR (qRT-PCR) validated the expression levels of prognosis-related genes in HNSC patients.

Results

We identified 4973 DEGs between HNSC and Normal samples. Key gene modules, represented by black and brown module genes, were identified. Intersection analysis revealed 76 GMRG-DEGs. Five prognosis-related genes (MTHFD2, CDKN2A, TPM2, MPZ, and DNMT1) were identified. A nomogram incorporating age, lymph node status (N), and risk score was constructed for survival prediction in HNSC patients. Immune infiltration analysis showed significant differences in five immune cell types (Macrophages M0, memory B cells, Monocytes, Macrophages M2, and Dendritic resting cells) between the high- and low-risk groups. GDSC database analysis identified 53 drugs with remarkable differences between the groups, including A.443654 and AG.014699. DNMT1 and MTHFD2 were up-regulated, while MPZ was down-regulated in HNSC.

Conclusion

Our study highlights the significant association of five prognosis-related genes (MTHFD2, CDKN2A, TPM2, MPZ, and DNMT1) with HNSC. These findings provide further evidence of the crucial role of GMRGs in HNSC.

MTHFD2 promotes PD-L1 expression via activation of the JAK/STAT signalling pathway in bladder cancer

Although combination chemotherapy is widely used for bladder cancer (BC) treatment, the recurrence and progression rates remain high. Therefore, novel therapeutic targets are required. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) contributes to tumourigenesis and immune evasion in several cancers; however, its biological function in BC remains unknown. This study aimed to investigate the expression, prognostic value and protumoural function of MTHFD2 in BC and elucidate the mechanism of programmed death-ligand 1 (PD-L1) upregulation by MTHFD2. An analysis using publicly available databases revealed that a high MTHFD2 expression was correlated with clinical features and a poor prognosis in BC. Furthermore, MTHFD2 promoted the growth, migration, invasion and tumourigenicity and decreased the apoptosis of BC cells in vivo and in vitro. The results obtained from databases showed that MTHFD2 expression was correlated with immune infiltration levels, PD-L1 expression, and the Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway. The expression of MTHFD2, PD-L1 and JAK/STAT signalling pathway-related proteins increased after interferon gamma treatment and decreased after MTHFD2 knockdown. Moreover, addition of a JAK/STAT pathway activator partially reduced the effect of MTHFD2 knockdown on BC cells. Collectively, our findings suggest that MTHFD2 promotes the expression of PD-L1 through the JAK/STAT signalling pathway in BC.

Induction of Reactive Bone Stromal Fibroblasts in 3D Models of Prostate Cancer Bone Metastases

A dynamic interplay between prostate cancer (PCa) cells and reactive bone stroma modulates the growth of metastases within the bone microenvironment. Of the stromal cells, metastasis-associated fibroblasts (MAFs) are known to contribute but are the least studied cell type in PCa tumour progression. It is the aim of the current study to establish a biologically relevant 3D in vitro model that mimics the cellular and molecular profiles of MAFs found in vivo. Using 3D in vitro cell culture models, the bone-derived fibroblast cell line, HS-5, was treated with conditioned media from metastatic-derived PCa cell lines, PC3 and MDA-PCa 2b, or mouse-derived fibroblasts 3T3. Two corresponding reactive cell lines were propagated: HS5-PC3 and HS5-MDA, and evaluated for alterations in morphology, phenotype, cellular behaviour, plus protein and genomic profiles. HS5-PC3 and HS5-MDA displayed distinct alterations in expression levels of N-Cadherin, non-functional E-Cadherin, alpha-smooth muscle actin (α-SMA), Tenascin C, and vimentin, along with transforming growth factor receptor expression (TGF β R1 and R2), consistent with subpopulations of MAFs reported in vivo. Transcriptomic analysis revealed a reversion of HS5-PC3 towards a metastatic phenotype with an upregulation in pathways known to regulate cancer invasion, proliferation, and angiogenesis. The exploitation of these engineered 3D models could help further unravel the novel biology regulating metastatic growth and the role fibroblasts play in the colonisation process.

Delineating the glioblastoma stemness by genes involved in cytoskeletal rearrangements and metabolic alterations

Literature data on glioblastoma ongoingly underline the link between metabolism and cancer stemness, the latter is one responsible for potentiating the resistance to treatment, inter alia due to increased invasiveness. In recent years, glioblastoma stemness research has bashfully introduced a key aspect of cytoskeletal rearrangements, whereas the impact of the cytoskeleton on invasiveness is well known. Although non-stem glioblastoma cells are less invasive than glioblastoma stem cells (GSCs), these cells also acquire stemness with greater ease if characterized as invasive cells and not tumor core cells. This suggests that glioblastoma stemness should be further investigated for any phenomena related to the cytoskeleton and metabolism, as they may provide new invasion-related insights. Previously, we proved that interplay between metabolism and cytoskeleton existed in glioblastoma. Despite searching for cytoskeleton-related processes in which the investigated genes might have been involved, not only did we stumble across the relation to metabolism but also reported genes that were found to be implicated in stemness. Thus, dedicated research on these genes in GSCs seems justifiable and might reveal novel directions and/or biomarkers that could be utilized in the future. Herein, we review the previously identified cytoskeleton/metabolism-related genes through the prism of glioblastoma stemness.

Regulation of RNA methylation by therapy treatment, promotes tumor survival

Our data previously revealed that chemosurviving cancer cells translate specific genes. Here, we find that the m6A-RNA-methyltransferase, METTL3, increases transiently in chemotherapy-treated breast cancer and leukemic cells in vitro and in vivo. Consistently, m6A increases on RNA from chemo-treated cells, and is needed for chemosurvival. This is regulated by eIF2α phosphorylation and mTOR inhibition upon therapy treatment. METTL3 mRNA purification reveals that eIF3 promotes METTL3 translation that is reduced by mutating a 5'UTR m6A-motif or depleting METTL3. METTL3 increase is transient after therapy treatment, as metabolic enzymes that control methylation and thus m6A levels on METTL3 RNA, are altered over time after therapy. Increased METTL3 reduces proliferation and anti-viral immune response genes, and enhances invasion genes, which promote tumor survival. Consistently, overriding phospho-eIF2α prevents METTL3 elevation, and reduces chemosurvival and immune-cell migration. These data reveal that therapy-induced stress signals transiently upregulate METTL3 translation, to alter gene expression for tumor survival.

The Metabolic Landscape of Breast Cancer and Its Therapeutic Implications

Breast cancer is the most common malignant tumor globally as of 2020 and remains the second leading cause of cancer-related death among female individuals worldwide. Metabolic reprogramming is well recognized as a hallmark of malignancy owing to the rewiring of multiple biological processes, notably, glycolysis, oxidative phosphorylation, pentose phosphate pathway, as well as lipid metabolism, which support the demands for the relentless growth of tumor cells and allows distant metastasis of cancer cells. Breast cancer cells are well documented to reprogram their metabolism via mutations or inactivation of intrinsic factors such as c-Myc, TP53, hypoxia-inducible factor, and the PI3K/AKT/mTOR pathway or crosstalk with the surrounding tumor microenvironments, including hypoxia, extracellular acidification and interaction with immune cells, cancer-associated fibroblasts, and adipocytes. Furthermore, altered metabolism contributes to acquired or inherent therapeutic resistance. Therefore, there is an urgent need to understand the metabolic plasticity underlying breast cancer progression as well as to dictate metabolic reprogramming that accounts for the resistance to standard of care. This review aims to illustrate the altered metabolism in breast cancer and its underlying mechanisms, as well as metabolic interventions in breast cancer treatment, with the intention to provide strategies for developing novel therapeutic treatments for breast cancer.

Identifying multiscale translational safety biomarkers using a network-based systems approach

Animal testing is the current standard for drug and chemicals safety assessment, but hazards translation to human is uncertain. Human in vitro models can address the species translation but might not replicate in vivo complexity. Herein, we propose a network-based method addressing these translational multiscale problems that derives in vivo liver injury biomarkers applicable to in vitro human early safety screening. We applied weighted correlation network analysis (WGCNA) to a large rat liver transcriptomic dataset to obtain co-regulated gene clusters (modules). We identified modules statistically associated with liver pathologies, including a module enriched for ATF4-regulated genes as associated with the occurrence of hepatocellular single-cell necrosis, and as preserved in human liver in vitro models. Within the module, we identified TRIB3 and MTHFD2 as a novel candidate stress biomarkers, and developed and used BAC-eGFPHepG2 reporters in a compound screening, identifying compounds showing ATF4-dependent stress response and potential early safety signals.

The Influence of Mitochondrial Energy and 1C Metabolism on the Efficacy of Anticancer Drugs: Exploring Potential Mechanisms of Resistance

Mitochondria are the main energy factory in living cells. To rapidly proliferate and metastasize, neoplastic cells increase their energy requirements. Thus, mitochondria become one of the most important organelles for them. Indeed, much research shows the interplay between cancer chemoresistance and altered mitochondrial function. In this review, we focus on the differences in energy metabolism between cancer and normal cells to better understand their resistance and how to develop drugs targeting energy metabolism and nucleotide synthesis. One of the differences between cancer and normal cells is the higher nicotinamide adenine dinucleotide (NAD+) level, a cofactor for the tricarboxylic acid cycle (TCA), which enhances their proliferation and helps cancer cells survive under hypoxic conditions. An important change is a metabolic switch called the Warburg effect. This effect is based on the change of energy harvesting from oxygen-dependent transformation to oxidative phosphorylation (OXPHOS), adapting them to the tumor environment. Another mechanism is the high expression of one-carbon (1C) metabolism enzymes. Again, this allows cancer cells to increase proliferation by producing precursors for the synthesis of nucleotides and amino acids. We reviewed drugs in clinical practice and development targeting NAD+, OXPHOS, and 1C metabolism. Combining novel drugs with conventional antineoplastic agents may prove to be a promising new way of anticancer treatment.

Identification of potential microRNAs regulating metabolic plasticity and cellular phenotypes in glioblastoma

MicroRNAs (miRNAs) play important role in regulating cellular metabolism, and are currently being explored in cancer. As metabolic reprogramming in cancer is a major mediator of phenotypic plasticity, understanding miRNA-regulated metabolism will provide opportunities to identify miRNA targets that can regulate oncogenic phenotypes by taking control of cellular metabolism. In the present work, we studied the effect of differentially expressed miRNAs on metabolism, and associated oncogenic phenotypes in glioblastoma (GBM) using patient-derived data. Networks of differentially expressed miRNAs and metabolic genes were created and analyzed to identify important miRNAs that regulate major metabolism in GBM. Graph network-based approaches like network diffusion, backbone extraction, and different centrality measures were used to analyze these networks for identification of potential miRNA targets. Important metabolic processes and cellular phenotypes were annotated to trace the functional responses associated with these miRNA-regulated metabolic genes and associated phenotype networks. miRNA-regulated metabolic gene subnetworks of cellular phenotypes were extracted, and important miRNAs regulating these phenotypes were identified. The most important outcome of the study is the target miRNA combinations predicted for five different oncogenic phenotypes that can be tested experimentally for miRNA-based therapeutic design in GBM. Strategies implemented in the study can be used to generate testable hypotheses in other cancer types as well, and design context-specific miRNA-based therapy for individual patient. Their usability can be further extended to other gene regulatory networks in cancer and other genetic diseases.

Increasing Dosage of Leucovorin Results in Pharmacokinetic and Gene Expression Differences When Administered as Two-Hour Infusion or Bolus Injection to Patients with Colon Cancer

The combination of 5-fluorouracil (5-FU) and leucovorin (LV) forms the chemotherapy backbone for patients with colorectal cancer. However, the LV administration is often standardized and not based on robust scientific data. To address these issues, a randomized pharmacokinetics study was performed in patients with colon cancer. Thirty patients were enrolled, receiving 60, 200 or 500 mg/m² LV as a single two-hour infusion. Blood, tumor, mucosa, and resection margin biopsies were collected. Folate concentrations were analyzed with LC-MS/MS and gene expression with qPCR. Data from a previous study where patients received LV as bolus injections were used as comparison. Saturation of methylenetetrahydrofolate (MeTHF) and tetrahydrofolate (THF) levels was seen after two-hour infusion and polyglutamated MeTHF + THF levels in tumors decreased with increasing LV dosage. The decrease was associated with decreased FPGS and increased GGH expression, which was not observed after LV bolus injection. In the bolus group, results indicate activation of a metabolic switch possibly promoting TYMS inhibition in response to 5-FU. Different metabolic mechanisms appear to be induced when LV is administered as infusion and bolus injection. Since maximal inhibition of TYMS by the 5-FU metabolite 5-fluoro-2'-deoxyuridine 5'-monophosphate (FdUMP) requires excess polyglutamated MeTHF, the results point in favor of the bolus regimen.

Folate enzyme MTHFD2 links one-carbon metabolism to unfolded protein response in glioblastoma

The mitochondrial folate enzyme methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) has shown oncogenic roles in various cancers and may have non-metabolic functions. This study investigated the role of MTHFD2 in glioblastoma pathogenesis. We find that MTHFD2 expression is enriched in gliomas by analysing public databases and clinical specimens. RNA interference (RNAi) and inhibitor of MTHFD2 hamper the proliferation of glioblastoma and induce apoptosis in cell lines, glioma stem-like cells (GSCs) and patient-derived xenografts (PDX). Metabolomic analyses show that MTHFD2 depletion suppresses the central carbon metabolic pathways, including glycolysis, the pentose phosphate pathway (PPP), and the tricarboxylic acid (TCA) cycle. GSEA reveals a novel non-metabolic function of MTHFD2 in association with the unfolded protein response (UPR). MTHFD2 depletion activates the PERK/eIF2α axis which contributes to translation inhibition and apoptosis; these effects are attenuated by a PERK inhibitor. Mechanistically, MTHFD2 may be linked to UPR via the post-transcriptionally regulation of chaperone protein GRP78. In conclusion, MTHFD2 could be a promising therapeutic target for glioblastoma. Besides its canonical role, MTHFD2 may contribute to glioblastoma pathogenesis via UPR, highlighting a newly identified functional link between one-carbon metabolism and cell stress response.

Metabolic collateral lethal target identification reveals MTHFD2 paralogue dependency in ovarian cancer

Recurrent loss-of-function deletions cause frequent inactivation of tumour suppressor genes but often also involve the collateral deletion of essential genes in chromosomal proximity, engendering dependence on paralogues that maintain similar function. Although these paralogues are attractive anticancer targets, no methodology exists to uncover such collateral lethal genes. Here we report a framework for collateral lethal gene identification via metabolic fluxes, CLIM, and use it to reveal MTHFD2 as a collateral lethal gene in UQCR11-deleted ovarian tumours. We show that MTHFD2 has a non-canonical oxidative function to provide mitochondrial NAD+, and demonstrate the regulation of systemic metabolic activity by the paralogue metabolic pathway maintaining metabolic flux compensation. This UQCR11-MTHFD2 collateral lethality is confirmed in vivo, with MTHFD2 inhibition leading to complete remission of UQCR11-deleted ovarian tumours. Using CLIM's machine learning and genome-scale metabolic flux analysis, we elucidate the broad efficacy of targeting MTHFD2 despite distinct cancer genetic profiles co-occurring with UQCR11 deletion and irrespective of stromal compositions of tumours.

ATF4/MYC Regulates MTHFD2 to Promote NSCLC Progression by Mediating Redox Homeostasis

Purpose. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) has been reported to be overexpressed in non-small-cell lung cancer (NSCLC) and to correlate with malignant proliferation. However, the mechanism of high MTHFD2 expression in NSCLC has not been clarified. Methods. qPCR, western blot, and immunofluorescence experiments were used to measure the expression of related mRNAs and proteins. Cell apoptosis was measured by flow cytometry and TUNEL assays. The CCK-8 assay was used to determine cell viability. Flow cytometry was used to analyze the cell cycle. ROS, H2O2, MDA, SOD, and NADPH/NADP+ were evaluated by relevant assay kits. Transfection of siRNA or vectors was used to downregulate or upregulate gene expression. Dual-luciferase reporter gene assays were used to evaluate the regulated relationship between MTHFD2 and ATF4 or MYC. Results. MTHFD2 was highly expressed in NSCLC cells. Knockdown of MTHFD2 inhibited proliferation and increased apoptosis. Furthermore, oxidative factors significantly increased, while antioxidant factors significantly decreased in NSCLC cells with MTHFD2 knockdown, indicating that MTHFD2 was involved in NSCLC progression through the redox pathway. Although MTHFD2 was downregulated with ATF4 silencing, the dual-luciferase reporter assay suggested that ATF4 did not directly mediate MTHFD2 transcription. Further studies revealed that MYC had a transcriptional effect on MTHFD2 and was also regulated by ATF4. PCR, and western blotting experiments with ATF4 knockdown and MYC overexpression as well as ATF4 overexpression and MYC knockdown proved that ATF4 stimulated MTHFD2 through MYC mediation. Conclusions. ATF4 promoted high expression of MTHFD2 in NSCLC dependent on MYC.

[Association of maternal MTHFD1 and MTHFD2 gene polymorphisms with congenital heart disease in offspring]

OBJECTIVES

To study the association of maternal methylenetetrahydrofolate dehydrogenase 1 (MTHFD1) and methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) gene polymorphisms with congenital heart disease (CHD) in offspring.

METHODS

A hospital-based case-control study was conducted. The mothers of 683 children with CHD alone who attended Hunan Children's Hospital, from November 2017 to March 2020 were enrolled as the case group, and the mothers of 740 healthy children who attended the same hospital during the same period and did not have any deformity were enrolled as the control group. A questionnaire survey was performed to collect related exposure data, and then venous blood samples (5 mL) were collected from the mothers to detect MTHFD1 and MTHFD2 gene polymorphisms. A multivariate logistic regression analysis was used to evaluate the association of MTHFD1 and MTHFD2 gene polymorphisms with CHD. The four-gamete test in Haploview 4.2 software was used to construct haplotypes and evaluate the association between haplotypes and CHD. The generalized multifactor dimensionality reduction method and logistic regression analysis were used to examine gene-gene interaction and its association with CHD.

RESULTS

The multivariate logistic regression analysis showed that maternal MTHFD1 gene polymorphisms at rs11849530 (GA vs AA: OR=1.49; GG vs AA: OR=2.04) andat rs1256142 (GA vs GG: OR=2.34; AA vs GG: OR=3.25) significantly increased the risk of CHD in offspring (P<0.05), while maternal MTHFD1 gene polymorphisms at rs1950902 (AA vs GG: OR=0.57) and MTHFD2 gene polymorphisms at rs1095966 (CA vs CC: OR=0.68) significantly reduced the risk of CHD in offspring (P<0.05). The haplotypes of G-G-G (OR=1.86) and G-A-G (OR=1.35) in mothers significantly increased the risk of CHD in offspring (P<0.05). The gene-gene interaction analyses showed that the first-order interaction between MTHFD1 rs1950902 and MTHFD1 rs2236222 and the second-order interaction involving MTHFD1 rs1950902, MTHFD1 rs1256142, and MTHFD2 rs1095966 might be associated with risk of CHD (P<0.05).

CONCLUSIONS

Maternal MTHFD1 and MTHFD2 gene polymorphisms and their haplotypes, as well as the interaction between MTHFD1 rs1950902 and MTHFD1 rs2236222 and between MTHFD1 rs1950902, MTHFD1 rs1256142, and MTHFD2 rs1095966, are associated with the risk of CHD in offspring.

Targeting mitochondrial one-carbon enzyme MTHFD2 together with pemetrexed confers therapeutic advantages in lung adenocarcinoma

Metabolic remodeling is the fundamental molecular feature of malignant tumors. Cancer cells require sufficient energy supplies supporting their high proliferative rate. MTHFD2, a mitochondrial one-carbon metabolic enzyme, is dysregulated in several malignancies and may serve as a promising therapeutic candidate in cancer treatment. Here, our data confirmed that MTHFD2 gene and protein was upregulated in the cancerous tissues of LUAD patients and was correlated with a poor survival in LUAD. MTHFD2 was involved in lung cancer cell proliferation, migration, and apoptosis by mediating its downstream molecules, such as DNA helicases (MCM4 and MCM7), as well as ZEB1, Vimentin and SNAI1, which contributed to tumor cell growth and epithelial-to-mesenchymal transition (EMT) process. Moreover, we identified that miRNA-99a-3p appeared to be an upstream mediator directly regulating MTHFD2 and MCM4 expression. Moreover, specific inhibition of MTHFD2 functions by siRNA or a chemical compound, improved anti-tumor sensitivities induced by pemetrexed in LUAD. Taken together, our study revealed the underlying molecular mechanisms of MTHFD2 in regulating cell proliferation and identified a novel therapeutic strategy improving the treatment efficacies in LUAD.

MTHFD2 is a potential oncogene for its strong association with poor prognosis and high level of immune infiltrates in urothelial carcinomas of bladder

BACKGROUND

The bifunctional methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 2, methenyltetrahydrofolate cyclohydrolase (MTHFD2) has been reported to play an oncogenic role in various types of cancers. However, the function of MTHFD2 in urothelial carcinomas of bladder (UCB) and its association with tumor immune infiltration remains unknown. We aim to examine the suitability of MTHFD2 to be a novel biomarker of bladder cancer and whether MTHFD2 is linked to immune infiltration.

METHODS

RNA sequencing data and clinical information (bladder cancer samples: normal samples = 414: 19) were downloaded from The Cancer Genome Atlas official website. Western blot analysis was performed to detect MTHFD2 expression in human bladder cancer (BLCA) cells and normal urothelial cell line SV-HUC-1. Associations between MTHFD2 expression and clinicopathological features were analyzed using Mann Whitney U test or Kruskal-Wallis H test. The "survival" and "survminer" packages were utilized to plot Kaplan-Meier survival curves. Moreover, the gene set enrichment analysis (GSEA) was conducted using a clusterProfiler package. The correlation of MTHFD2 expression with immune infiltration level was estimated using the single sample GSEA (ssGSEA) algorithm. Furthermore, associations between MTHFD2 and immune checkpoint genes were evaluated using the correlation analysis.

RESULTS

Transcriptome analysis manifested that MTHFD2 was highly expressed in UCB tissues than normal bladder tissues, which was further confirmed by western blot analysis in human BLCA cells and SV-HUC-1 cells. Moreover, MTHFD2 high expression was significantly associated with the advanced disease progression. Also, the high expression of MTHFD2 was correlated with poor prognosis, and MTHFD2 was considered as an independent prognostic factor for disease specific survival. Furthermore, a number of cancer-related pathways were enriched in MTHFD2 high group, including NF-κB activation, JAK/STAT, and cancer immunotherapy by PD1 blockade. Several immune checkpoint molecules were also strongly associated with MTHFD2 expression, including PDCD1, CD274, CTLA4, CD276, LAG3, HAVCR2, and TIGIT.

CONCLUSIONS

MTHFD2 expression was remarkably elevated in UCB, suggesting that MTHFD2 could be a promising biomarker for BLCA as well as novel target for anti-cancer immunotherapy since its close association with immune infiltration.

Proteomic profiling of metformin effects in 3T3-L1 adipocytes by SILAC-based quantification

Metformin is a common and generally the first medication prescribed for treatment of type 2 diabetes. Its mechanism involves affecting pathways that regulate glucose and lipid metabolism in metabolic cells such as that of muscle and liver cells. In spite of various studies exploring its effects, the proteome changes in adipocytes in response to metformin remains poorly understood. In this study, we performed SILAC-based quantitative proteomic profiling to study the effects of metformin specifically on 3T3-L1 adipocytes. We define proteins that exhibited altered levels with metformin treatment, 400 of them showing statistically significant changes in our study. Our results suggest that metformin affects not only the PPARγ signaling pathway, as well as glucose and lipid metabolism, but also protein folding, endoplasmic reticulum stress, negative regulation of appetite, and one-carbon folate metabolism in adipocytes. This proteomic investigation provides important insight into effects of metformin in adipocytes. This article is protected by copyright. All rights reserved.

MTHFD2 is a metabolic checkpoint controlling effector and regulatory T cell fate and function

Antigenic stimulation promotes T cell metabolic reprogramming to meet increased biosynthetic, bioenergetic, and signaling demands. We show that the one-carbon (1C) metabolism enzyme methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) regulates de novo purine synthesis and signaling in activated T cells to promote proliferation and inflammatory cytokine production. In pathogenic T helper-17 (Th17) cells, MTHFD2 prevented aberrant upregulation of the transcription factor FoxP3 along with inappropriate gain of suppressive capacity. MTHFD2 deficiency also promoted regulatory T (Treg) cell differentiation. Mechanistically, MTHFD2 inhibition led to depletion of purine pools, accumulation of purine biosynthetic intermediates, and decreased nutrient sensor mTORC1 signaling. MTHFD2 was also critical to regulate DNA and histone methylation in Th17 cells. Importantly, MTHFD2 deficiency reduced disease severity in multiple in vivo inflammatory disease models. MTHFD2 is thus a metabolic checkpoint to integrate purine metabolism with pathogenic effector cell signaling and is a potential therapeutic target within 1C metabolism pathways.

One-Carbon Metabolism: Pulling the Strings behind Aging and Neurodegeneration

One-carbon metabolism (OCM) is a network of biochemical reactions delivering one-carbon units to various biosynthetic pathways. The folate cycle and methionine cycle are the two key modules of this network that regulate purine and thymidine synthesis, amino acid homeostasis, and epigenetic mechanisms. Intersection with the transsulfuration pathway supports glutathione production and regulation of the cellular redox state. Dietary intake of micronutrients, such as folates and amino acids, directly contributes to OCM, thereby adapting the cellular metabolic state to environmental inputs. The contribution of OCM to cellular proliferation during development and in adult proliferative tissues is well established. Nevertheless, accumulating evidence reveals the pivotal role of OCM in cellular homeostasis of non-proliferative tissues and in coordination of signaling cascades that regulate energy homeostasis and longevity. In this review, we summarize the current knowledge on OCM and related pathways and discuss how this metabolic network may impact longevity and neurodegeneration across species.

A novel oral inhibitor for one-carbon metabolism and checkpoint kinase 1 inhibitor as a rational combination treatment for breast cancer

Patients with triple-negative breast cancer have a poor prognosis as only a few efficient targeted therapies are available. Cancer cells are characterized by their unregulated proliferation and require large amounts of nucleotides to replicate their DNA. One-carbon metabolism contributes to purine and pyrimidine nucleotide synthesis by supplying one carbon atom. Although mitochondrial one-carbon metabolism has recently been focused on as an important target for cancer treatment, few specific inhibitors have been reported. In this study, we aimed to examine the effects of DS18561882 (DS18), a novel, orally active, specific inhibitor of methylenetetrahydrofolate dehydrogenase (MTHFD2), a mitochondrial enzyme involved in one-carbon metabolism. Treatment with DS18 led to a marked reduction in cancer-cell proliferation; however, it did not induce cell death. Combinatorial treatment with DS18 and inhibitors of checkpoint kinase 1 (Chk1), an activator of the S phase checkpoint pathway, efficiently induced apoptotic cell death in breast cancer cells and suppressed tumorigenesis in a triple-negative breast cancer patient-derived xenograft model. Mechanistically, MTHFD2 inhibition led to cell cycle arrest and slowed nucleotide synthesis. This finding suggests that DNA replication stress occurs due to nucleotide shortage and that the S-phase checkpoint pathway is activated, leading to cell-cycle arrest. Combinatorial treatment with both inhibitors released cell-cycle arrest, but induced accumulation of DNA double-strand breaks, leading to apoptotic cell death. Collectively, a combination of MTHFD2 and Chk1 inhibitors would be a rational treatment option for patients with triple-negative breast cancer.

Methylene tetrahydrofolate dehydrogenase 2 (MTHFD2) is overexpressed in head and neck squamous cell carcinoma (HNSCC) and correlated with patient’s poor prognosis

Objective

To investigate methylene tetrahydrofolate dehydrogenase 2 (MTHFD2) expression, biological function, and correlation with head and neck squamous cell carcinoma (HNSCC) patient’s prognosis.

Methods

The relative expression levels of MTHFD2 gene mRNA in tumor tissues of HNSCC and adjacent normal tissues were analyzed in the Cancer Genome Atlas and oncomine database. MTHFD2 protein relative expression in tumor tissue of HNSCC patients was analyzed in human proteome database. Protein–protein interaction (PPI) network of MTHFD2 and correlated genes were constructed in STRING database. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway of MTHFD2 and relevant proteins involved in the PPI network was enriched. The Tumor Immune Estimation Resource database was used to analyze the relationship between MTHFD2 expression and immune infiltration. Overall survival (OS) and progression-free survival (PFS) for MTHFD2 high and low expression groups were investigated in the Kaplan–Meier Plotter database.

Results

In HNSCC, MTHFD2 mRNA relative expression level in tumor tissue was significantly higher than the corresponding normal tissue with statistical difference (p < 0.05). In the PPI network, 21 protein coding genes were involved in the network with 124 edges, which indicated that the enrichment was significant (p < 0.05). MTHFD2 and PPI network involved genes were mainly enriched in tetrahydrofolate metabolic process, one-carbon metabolic process biological process. In KEGG pathway, MTHFD2 and PPI network involved genes were mainly enriched in one-carbon pool by folate, metabolic pathways, glyoxylate, and dicarboxylate metabolism, and carbon metabolism. The relative expression level of MTHFD2 gene was correlated with immune infiltration of macrophage (r = 0.712, p < 0.05), neutrophil (r = 0.158, p < 0.05), dendritic cell (r = 0.1825, p < 0.05), and CD4+ T lymph cell (r = 0.1825, p < 0.05). HNSCC patients with high expression MTHFD2 had low OS compared to low expression cases (hazard ratio = 1.53, 95% CI: 1.16–2.02, p < 0.05).

Conclusion

MTHFD2 is overexpressed in HNSCC and correlated with patient’s prognosis. MTHFD2 maybe a potential target for HNSCC target treatment and provides a possible direction for the research and development of related targeted drugs.

Non-metabolic function of MTHFD2 activates CDK2 in bladder cancer

Bladder cancer is a common tumor with a high recurrence rate and high fatality rate, and its mechanism of occurrence and development remains unclear. Many proteins and metabolites reprogram at different stages of tumor development to support tumor cell growth. The moonlighting effect happens when a protein performs multiple functions simultaneously in a cell. In this study, we identified a metabolic protein, MTHFD2, which participates in the cell cycle by binding to CDK2 in bladder cancer. MTHFD2 has been shown to affect bladder cancer cell growth, which is independent of its metabolic function. We found that MTHFD2 was involved in cell cycle regulation and could encourage cell cycle progression by activating CDK2 and sequentially affecting E2F1 activation. In addition, moonlighting MTHFD2 might be regulated by the dynamics of the mitochondria. In conclusion, MTHFD2 localizes in the nucleus to perform a distinct function of catalyzing metabolic reactions. Moreover, the nuclear MTHFD2 activates CDK2 and promotes bladder cancer cell growth by modulating the cell cycle.

D-SCRIPT translates genome to phenome with sequence-based, structure-aware, genome-scale predictions of protein-protein interactions

We combine advances in neural language modeling and structurally motivated design to develop D-SCRIPT, an interpretable and generalizable deep-learning model, which predicts interaction between two proteins using only their sequence and maintains high accuracy with limited training data and across species. We show that a D-SCRIPT model trained on 38,345 human PPIs enables significantly improved functional characterization of fly proteins compared with the state-of-the-art approach. Evaluating the same D-SCRIPT model on protein complexes with known 3D structure, we find that the inter-protein contact map output by D-SCRIPT has significant overlap with the ground truth. We apply D-SCRIPT to screen for PPIs in cow (Bos taurus) at a genome-wide scale and focusing on rumen physiology, identify functional gene modules related to metabolism and immune response. The predicted interactions can then be leveraged for function prediction at scale, addressing the genome-to-phenome challenge, especially in species where little data are available.

Transcriptome Analysis of Effects of Folic Acid Supplement on Gene Expression in Liver of Broiler Chickens

Folic acid is a water-soluble B vitamin, and plays an important role in regulating gene expression and methylation. The liver is the major site of lipid biosynthesis in the chicken. Nevertheless, how gene expression and regulatory networks are affected by folic acid in liver of broilers are poorly understood. This paper conducted the RNA-seq technology on the liver of broilers under folic acid challenge investigation. First, 405 differentially expressed genes (DEGs), including 157 significantly upregulated and 248 downregulated, were detected between the control group (C) and the 5 mg folic acid group (M). Second, 68 upregulated DEGs and 142 downregulated DEGs were determined between C group and 10 mg folic acid group (H). Third, there were 165 upregulated genes and 179 downregulated genes between M and H groups. Of these DEGs, 903 DEGs were successfully annotated in the public databases. The functional classification based on GO and KEEGG showed that “general function prediction only” represented the largest functional classes, “cell cycle” (C vs. M; M vs. H), and “neuroactive ligand-receptor interaction” (C vs. H) were the highest unique sequences among three groups. SNP analysis indicated that numbers of C, M and H groups were 145,450, 146,131, and 123,004, respectively. Total new predicted alternative splicing events in C, M, and H groups were 9,521, 9,328, and 8,929, respectively. A protein-protein interaction (PPI) network was constructed, and the top 10 hub genes were evaluated among three groups. The results of real time PCR indicated that mRNA abundance of PPARγ and FAS in abdominal fat of M and H groups were reduced compared with the C group (P < 0.05). Ultramicroscopy results showed that folic acid could reduce lipid droplets in livers from chickens. Finally, contents of LPL, PPARγ, and FAS in abdominal fat were decreased with the folic acid supplmented diets (P < 0.01). These findings reveal the effects of folic acid supplemention on gene expression in liver of broilers, which can provide information for understanding the molecular mechanisms of folic acid regulating liver lipid metabolism.

Proof of Gene Doping in a Mouse Model with a Human Erythropoietin Gene Transferred Using an Adenoviral Vector

Despite the World Anti-Doping Agency (WADA) ban on gene doping in the context of advancements in gene therapy, the risk of EPO gene-based doping among athletes is still present. To address this and similar risks, gene-doping tests are being developed in doping control laboratories worldwide. In this regard, the present study was performed with two objectives: to develop a robust gene-doping mouse model with the human EPO gene (hEPO) transferred using recombinant adenovirus (rAdV) as a vector and to develop a detection method to identify gene doping by using this model. The rAdV including the hEPO gene was injected intravenously to transfer the gene to the liver. After injection, the mice showed significantly increased whole-blood red blood cell counts and increased expression of hematopoietic marker genes in the spleen, indicating successful development of the gene-doping model. Next, direct and potentially indirect proof of gene doping were evaluated in whole-blood DNA and RNA by using a quantitative PCR assay and RNA sequencing. Proof of doping could be detected in DNA and RNA samples from one drop of whole blood for approximately a month; furthermore, the overall RNA expression profiles showed significant changes, allowing advanced detection of hEPO gene doping.

Insights on Metabolic Reprogramming and Its Therapeutic Potential in Acute Leukemia

Acute leukemias, classified as acute myeloid leukemia and acute lymphoblastic leukemia, represent the most prevalent hematologic tumors in adolescent and young adults. In recent years, new challenges have emerged in order to improve the clinical effectiveness of therapies already in use and reduce their side effects. In particular, in this scenario, metabolic reprogramming plays a key role in tumorigenesis and prognosis, and it contributes to the treatment outcome of acute leukemia. This review summarizes the latest findings regarding the most relevant metabolic pathways contributing to the continuous growth, redox homeostasis, and drug resistance of leukemia cells. We describe the main metabolic deregulations in acute leukemia and evidence vulnerabilities that could be exploited for targeted therapy.

Xanthine Derivatives Reveal an Allosteric Binding Site in Methylenetetrahydrofolate Dehydrogenase 2 (MTHFD2)

Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) plays an important role in one-carbon metabolism. The MTHFD2 gene is upregulated in various cancers but very low or undetectable in normal proliferating cells, and therefore a potential target for cancer treatment. In this study, we present the structure of MTHFD2 in complex with xanthine derivative 15, which allosterically binds to MTHFD2 and coexists with the substrate analogue. A kinetic study demonstrated the uncompetitive inhibition of MTHFD2 by 15. Allosteric inhibitors often provide good selectivity and, indeed, xanthine derivatives are highly selective for MTHFD2. Moreover, several conformational changes were observed upon the binding of 15, which impeded the binding of the cofactor and phosphate to MTHFD2. To the best of our knowledge, this is the first study to identify allosteric inhibitors targeting the MTHFD family and our results would provide insights on the inhibition mechanism of MTHFD proteins and the development of novel inhibitors.

PLEK2, RRM2, GCSH: A Novel WWOX-Dependent Biomarker Triad of Glioblastoma at the Crossroads of Cytoskeleton Reorganization and Metabolism Alterations

Glioblastoma is one of the deadliest human cancers. Its malignancy depends on cytoskeleton reorganization, which is related to, e.g., epithelial-to-mesenchymal transition and metastasis. The malignant phenotype of glioblastoma is also affected by the WWOX gene, which is lost in nearly a quarter of gliomas. Although the role of WWOX in the cytoskeleton rearrangement has been found in neural progenitor cells, its function as a modulator of cytoskeleton in gliomas was not investigated. Therefore, this study aimed to investigate the role of WWOX and its collaborators in cytoskeleton dynamics of glioblastoma. Methodology on RNA-seq data integrated the use of databases, bioinformatics tools, web-based platforms, and machine learning algorithm, and the obtained results were validated through microarray data. PLEK2, RRM2, and GCSH were the most relevant WWOX-dependent genes that could serve as novel biomarkers. Other genes important in the context of cytoskeleton (BMP4, CCL11, CUX2, DUSP7, FAM92B, GRIN2B, HOXA1, HOXA10, KIF20A, NF2, SPOCK1, TTR, UHRF1, and WT1), metabolism (MTHFD2), or correlation with WWOX (COL3A1, KIF20A, RNF141, and RXRG) were also discovered. For the first time, we propose that changes in WWOX expression dictate a myriad of alterations that affect both glioblastoma cytoskeleton and metabolism, rendering new therapeutic possibilities.

Metabolic Reprogramming of Colorectal Cancer Cells and the Microenvironment: Implication for Therapy

Colorectal carcinoma (CRC) is one of the most frequently diagnosed carcinomas and one of the leading causes of cancer-related death worldwide. Metabolic reprogramming, a hallmark of cancer, is closely related to the initiation and progression of carcinomas, including CRC. Accumulating evidence shows that activation of oncogenic pathways and loss of tumor suppressor genes regulate the metabolic reprogramming that is mainly involved in glycolysis, glutaminolysis, one-carbon metabolism and lipid metabolism. The abnormal metabolic program provides tumor cells with abundant energy, nutrients and redox requirements to support their malignant growth and metastasis, which is accompanied by impaired metabolic flexibility in the tumor microenvironment (TME) and dysbiosis of the gut microbiota. The metabolic crosstalk between the tumor cells, the components of the TME and the intestinal microbiota further facilitates CRC cell proliferation, invasion and metastasis and leads to therapy resistance. Hence, to target the dysregulated tumor metabolism, the TME and the gut microbiota, novel preventive and therapeutic applications are required. In this review, the dysregulation of metabolic programs, molecular pathways, the TME and the intestinal microbiota in CRC is addressed. Possible therapeutic strategies, including metabolic inhibition and immune therapy in CRC, as well as modulation of the aberrant intestinal microbiota, are discussed.

A Review of Small-Molecule Inhibitors of One-Carbon Enzymes: SHMT2 and MTHFD2 in the Spotlight

Metabolic reprogramming is a key hallmark of cancer and shifts cellular metabolism to meet the demands of biomass production necessary for abnormal cell reproduction. One-carbon metabolism (1CM) contributes to many biosynthetic pathways that fuel growth and is comprised of a complex network of enzymes. Methotrexate and 5-fluorouracil were pioneering drugs in this field and are still widely used today as anticancer agents as well as for other diseases such as arthritis. Besides dihydrofolate reductase and thymidylate synthase, two other enzymes of the folate cycle arm of 1CM have not been targeted clinically: serine hydroxymethyltransferase (SHMT) and methylenetetrahydrofolate dehydrogenase (MTHFD). An increasing body of literature suggests that the mitochondrial isoforms of these enzymes (SHMT2 and MTHFD2) are clinically relevant in the context of cancer. In this review, we focused on the 1CM pathway as a target for cancer therapy and, in particular, SHMT2 and MTHFD2. The function, regulation, and clinical relevance of SHMT2 and MTHFD2 are all discussed. We expand on previous clinical studies and evaluate the prognostic significance of these critical enzymes by performing a pan-cancer analysis of patient data from the The Cancer Genome Atlas and a transcriptional coexpression network enrichment analysis. We also provide an overview of preclinical and clinical inhibitors targeting the folate pathway, the methionine cycle, and folate-dependent purine biosynthesis enzymes.

One-carbon metabolism in cancer cells: a critical review based on a core model of central metabolism

One-carbon metabolism (1C-metabolism), also called folate metabolism because the carbon group is attached to folate-derived tetrahydrofolate, is crucial in metabolism. It is at the heart of several essential syntheses, particularly those of purine and thymidylate. After a short reminder of the organization of 1C-metabolism, I list its salient features as reported in the literature. Then, using flux balance analysis, a core model of central metabolism and the flux constraints for an 'average cancer cell metabolism', I explore the fundamentals underlying 1C-metabolism and its relationships with the rest of metabolism. Some unreported properties of 1C-metabolism emerge, such as its potential roles in mitochondrial NADH exchange with cytosolic NADPH, participation in NADH recycling, and optimization of cell proliferation.

Context-Dependent Roles of RNA Modifications in Stress Responses and Diseases

RNA modifications are diverse post-transcriptional modifications that regulate RNA metabolism and gene expression. RNA modifications, and the writers, erasers, and readers that catalyze these modifications, serve as important signaling machineries in cellular stress responses and disease pathogenesis. In response to stress, RNA modifications are mobilized to activate or inhibit the signaling pathways that combat stresses, including oxidative stress, hypoxia, therapeutic stress, metabolic stress, heat shock, DNA damage, and ER stress. The role of RNA modifications in response to these cellular stressors is context- and cell-type-dependent. Due to their pervasive roles in cell biology, RNA modifications have been implicated in the pathogenesis of different diseases, including cancer, neurologic and developmental disorders and diseases, and metabolic diseases. In this review, we aim to summarize the roles of RNA modifications in molecular and cellular stress responses and diseases.

A systematic analysis of a potential metabolism-related prognostic signature for breast cancer patients

Background

Metabolic pathways play an essential role in breast cancer. However, the role of metabolism-related genes in the early diagnosis of breast cancer remains unknown.

Methods

In our study, RNA sequencing (RNA-seq) expression data and clinicopathological information from The Cancer Genome Atlas (TCGA) and GSE20685 were obtained. Univariate cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses were performed on the differentially expressed metabolism-related genes. Then, the formula of the metabolism-related risk model was composed, and the risk score of each patient was calculated. The breast cancer patients were divided into high-risk and low-risk groups with a cutoff of the median expression value of the risk score, and the prognostic analysis was also used to analyze the survival time between these two groups. In the end, we also analyzed the expression, interaction, and correlation among genes in the metabolism-related gene risk model.

Results

The results from the prognostic analysis indicated that the survival was significantly poorer in the high-risk group than in the low-risk group in both TCGA and GSE20685 datasets. In addition, after adjusting for different clinicopathological features in multivariate analysis, the metabolism-related risk model remained an independent prognostic indicator in TCGA dataset.

Conclusions

In summary, we systematically developed a potential metabolism-related gene risk model for predicting prognosis in breast cancer patients.