The PLOD2/succinate axis regulates the epithelial-mesenchymal plasticity and cancer cell stemness.

Aberrant accumulation of succinate has been detected in many cancers. However, the cellular function and regulation of succinate in cancer progression is not completely understood. Using stable isotope-resolved metabolomics analysis, we showed that the epithelial mesenchymal transition (EMT) was associated with profound changes in metabolites, including elevation of cytoplasmic succinate levels. The treatment with cell-permeable succinate induced mesenchymal phenotypes in mammary epithelial cells and enhanced cancer cell stemness. Chromatin immunoprecipitation and sequence analysis showed that elevated cytoplasmic succinate levels were sufficient to reduce global 5-hydroxymethylcytosinene (5hmC) accumulation and induce transcriptional repression of EMT-related genes. We showed that expression of procollagen-lysine,2-oxoglutarate 5-dioxygenase 2 (PLOD2) was associated with elevation of cytoplasmic succinate during the EMT process. Silencing of PLOD2 expression in breast cancer cells reduced succinate levels and inhibited cancer cell mesenchymal phenotypes and stemness, which was accompanied by elevated 5hmC levels in chromatin. Importantly, exogenous succinate rescued cancer cell stemness and 5hmC levels in PLOD2-silenced cells, suggesting that PLOD2 promotes cancer progression at least partially through succinate. These results reveal the previously unidentified function of succinate in enhancing cancer cell plasticity and stemness.

Metabolic reprogramming driven by EZH2 inhibition depends on cell-matrix interactions.

EZH2 (Enhancer of Zeste Homolog 2), a subunit of Polycomb Repressive Complex 2 (PRC2), catalyzes the trimethylation of histone H3 at lysine 27 (H3K27me3), which represses expression of genes. It also has PRC2-independent functions, including transcriptional coactivation of oncogenes, and is frequently overexpressed in lung cancers. Clinically, EZH2 inhibition can be achieved with the FDA-approved drug EPZ-6438 (tazemetostat). To realize the full potential of EZH2 blockade, it is critical to understand how cell-cell/cell-matrix interactions present in 3D tissue and cell culture systems influences this blockade in terms of growth-related metabolic functions. Here, we show that EZH2 suppression reduced growth of human lung adenocarcinoma A549 cells in 2D cultures but stimulated growth in 3D cultures. To understand the metabolic underpinnings, we employed [13C6]-glucose stable isotope-resolved metabolomics to determine the effect of EZH2 suppression on metabolic networks in 2D versus 3D A549 cultures. The Krebs cycle, neoribogenesis, γ-aminobutyrate metabolism, and salvage synthesis of purine nucleotides were activated by EZH2 suppression in 3D spheroids but not in 2D cells, consistent with the growth effect. Using simultaneous 2H7-glucose + 13C5,15N2-Gln tracers and EPZ-6438 inhibition of H3 trimethylation, we delineated the effects on the Krebs cycle, γ-aminobutyrate metabolism, gluconeogenesis, and purine salvage to be PRC2-dependent. Furthermore, the growth/metabolic effects differed for mouse Matrigel versus self-produced A549 extracellular matrix. Thus, our findings highlight the importance of the presence and nature of extracellular matrix in studying the function of EZH2 and its inhibitors in cancer cells for modeling the in vivo outcomes.

Cyclin D1 extensively reprograms metabolism to support biosynthetic pathways in hepatocytes.

Cell proliferation requires metabolic reprogramming to accommodate biosynthesis of new cell components, and similar alterations occur in cancer cells. However, the mechanisms linking the cell cycle machinery to metabolism are not well defined. Cyclin D1, along with its main partner cyclin-dependent kinase 4 (Cdk4), is a pivotal cell cycle regulator and driver oncogene that is overexpressed in many cancers. Here, we examine hepatocyte proliferation to define novel effects of cyclin D1 on biosynthetic metabolism. Metabolomic studies reveal that cyclin D1 broadly promotes biosynthetic pathways including glycolysis, the pentose phosphate pathway, and the purine and pyrimidine nucleotide synthesis in hepatocytes. Proteomic analyses demonstrate that overexpressed cyclin D1 binds to numerous metabolic enzymes including those involved in glycolysis and pyrimidine synthesis. In the glycolysis pathway, cyclin D1 activates aldolase and GAPDH, and these proteins are phosphorylated by cyclin D1/Cdk4 in vitro. De novo pyrimidine synthesis is particularly dependent on cyclin D1. Cyclin D1/Cdk4 phosphorylates the initial enzyme of this pathway, carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase (CAD), and metabolomic analysis indicates that cyclin D1 depletion markedly reduces the activity of this enzyme. Pharmacologic inhibition of Cdk4 along with the downstream pyrimidine synthesis enzyme dihydroorotate dehydrogenase synergistically inhibits proliferation and survival of hepatocellular carcinoma cells. These studies demonstrate that cyclin D1 promotes a broad network of biosynthetic pathways in hepatocytes, and this model may provide insights into potential metabolic vulnerabilities in cancer cells.

Polarization and ß-Glucan Reprogram Immunomodulatory Metabolism in Human Macrophages and Ex Vivo in Human Lung Cancer Tissues.

Immunomodulatory (IM) metabolic reprogramming in macrophages (Mϕs) is fundamental to immune function. However, limited information is available for human Mϕs, particularly in response plasticity, which is critical to understanding the variable efficacy of immunotherapies in cancer patients. We carried out an in-depth analysis by combining multiplex stable isotope-resolved metabolomics with reversed phase protein array to map the dynamic changes of the IM metabolic network and key protein regulators in four human donors' Mϕs in response to differential polarization and M1 repolarizer ß-glucan (whole glucan particles [WGPs]). These responses were compared with those of WGP-treated ex vivo organotypic tissue cultures (OTCs) of human non-small cell lung cancer. We found consistently enhanced tryptophan catabolism with blocked NAD+ and UTP synthesis in M1-type Mϕs (M1-Mϕs), which was associated with immune activation evidenced by increased release of IL-1ß/CXCL10/IFN-γ/TNF-α and reduced phagocytosis. In M2a-Mϕs, WGP treatment of M2a-Mϕs robustly increased glucose utilization via the glycolysis/oxidative branch of the pentose phosphate pathway while enhancing UDP-N-acetyl-glucosamine turnover and glutamine-fueled gluconeogenesis, which was accompanied by the release of proinflammatory IL-1ß/TNF-α to above M1-Mϕ's levels, anti-inflammatory IL-10 to above M2a-Mϕ's levels, and attenuated phagocytosis. These IM metabolic responses could underlie the opposing effects of WGP, i.e., reverting M2- to M1-type immune functions but also boosting anti-inflammation. Variable reprogrammed Krebs cycle and glutamine-fueled synthesis of UTP in WGP-treated OTCs of human non-small cell lung cancer were observed, reflecting variable M1 repolarization of tumor-associated Mϕs. This was supported by correlation with IL-1ß/TNF-α release and compromised tumor status, making patient-derived OTCs unique models for studying variable immunotherapeutic efficacy in cancer patients.

Bayesian kinetic modeling for tracer-based metabolomic data.

BACKGROUND: Stable Isotope Resolved Metabolomics (SIRM) is a new biological approach that uses stable isotope tracers such as uniformly [Formula: see text]-enriched glucose ([Formula: see text]-Glc) to trace metabolic pathways or networks at the atomic level in complex biological systems. Non-steady-state kinetic modeling based on SIRM data uses sets of simultaneous ordinary differential equations (ODEs) to quantitatively characterize the dynamic behavior of metabolic networks. It has been increasingly used to understand the regulation of normal metabolism and dysregulation in the development of diseases. However, fitting a kinetic model is challenging because there are usually multiple sets of parameter values that fit the data equally well, especially for large-scale kinetic models. In addition, there is a lack of statistically rigorous methods to compare kinetic model parameters between different experimental groups. RESULTS: We propose a new Bayesian statistical framework to enhance parameter estimation and hypothesis testing for non-steady-state kinetic modeling of SIRM data. For estimating kinetic model parameters, we leverage the prior distribution not only to allow incorporation of experts' knowledge but also to provide robust parameter estimation. We also introduce a shrinkage approach for borrowing information across the ensemble of metabolites to stably estimate the variance of an individual isotopomer. In addition, we use a component-wise adaptive Metropolis algorithm with delayed rejection to perform efficient Monte Carlo sampling of the posterior distribution over high-dimensional parameter space. For comparing kinetic model parameters between experimental groups, we propose a new reparameterization method that converts the complex hypothesis testing problem into a more tractable parameter estimation problem. We also propose an inference procedure based on credible interval and credible value. Our method is freely available for academic use at https://github.com/xuzhang0131/MCMCFlux . CONCLUSIONS: Our new Bayesian framework provides robust estimation of kinetic model parameters and enables rigorous comparison of model parameters between experimental groups. Simulation studies and application to a lung cancer study demonstrate that our framework performs well for non-steady-state kinetic modeling of SIRM data.

NMR-based isotope editing, chemoselection and isotopomer distribution analysis in stable isotope resolved metabolomics.

NMR is a very powerful tool for identifying and quantifying compounds within complex mixtures without the need for individual standards or chromatographic separation. Stable Isotope Resolved Metabolomics (or SIRM) is an approach to following the fate of individual atoms from precursors through metabolic transformation, producing an atom-resolved metabolic fate map. However, extracts of cells or tissue give rise to very complex NMR spectra. While multidimensional NMR experiments may partially overcome the spectral overlap problem, additional tools may be needed to determine site-specific isotopomer distributions. NMR is especially powerful by virtue of its isotope editing capabilities using NMR active nuclei such as 13C, 15N, 19F and 31P to select molecules containing just these atoms in a complex mixture, and provide direct information about which atoms are present in identified compounds and their relative abundances. The isotope-editing capability of NMR can also be employed to select for those compounds that have been selectively derivatized with an NMR-active stable isotope at particular functional groups, leading to considerable spectral simplification. Here we review isotope analysis by NMR, and methods of chemoselection both for spectral simplification, and for enhanced isotopomer analysis.

A negative reciprocal regulatory axis between cyclin D1 and HNF4α modulates cell cycle progression and metabolism in the liver.

Hepatocyte nuclear factor 4α (HNF4α) is a master regulator of liver function and a tumor suppressor in hepatocellular carcinoma (HCC). In this study, we explore the reciprocal negative regulation of HNF4α and cyclin D1, a key cell cycle protein in the liver. Transcriptomic analysis of cultured hepatocyte and HCC cells found that cyclin D1 knockdown induced the expression of a large network of HNF4α-regulated genes. Chromatin immunoprecipitation-sequencing (ChIP-seq) demonstrated that cyclin D1 inhibits the binding of HNF4α to thousands of targets in the liver, thereby diminishing the expression of associated genes that regulate diverse metabolic activities. Conversely, acute HNF4α deletion in the liver induces cyclin D1 and hepatocyte cell cycle progression; concurrent cyclin D1 ablation blocked this proliferation, suggesting that HNF4α maintains proliferative quiescence in the liver, at least, in part, via repression of cyclin D1. Acute cyclin D1 deletion in the regenerating liver markedly inhibited hepatocyte proliferation after partial hepatectomy, confirming its pivotal role in cell cycle progression in this in vivo model, and enhanced the expression of HNF4α target proteins. Hepatocyte cyclin D1 gene ablation caused markedly increased postprandial liver glycogen levels (in a HNF4α-dependent fashion), indicating that the cyclin D1-HNF4α axis regulates glucose metabolism in response to feeding. In AML12 hepatocytes, cyclin D1 depletion led to increased glucose uptake, which was negated if HNF4α was depleted simultaneously, and markedly elevated glycogen synthesis. To summarize, mutual repression by cyclin D1 and HNF4α coordinately controls the cell cycle machinery and metabolism in the liver.

Resolving Enantiomers of 2-Hydroxy Acids by Nuclear Magnetic Resonance.

Biologically important 2-hydroxy carboxylates such as lactate, malate, and 2-hydroxyglutarate exist in two enantiomeric forms that cannot be distinguished under achiral conditions. The D and L (or R, S) enantiomers have different biological origins and functions, and therefore, there is a need for a simple method for resolving, identifying, and quantifying these enantiomers. We have adapted and improved a chiral derivatization technique for nuclear magnetic resonance (NMR), which needs no chromatography for enantiomer resolution, with greater than 90% overall recovery. This method was developed for 2-hydroxyglutarate (2HG) to produce diastereomers resolvable by column chromatography. We have applied the method to lactate, malate, and 2HG. The limit of quantification was determined to be about 1 nmol for 2HG with coefficients of variation of less than 5%. We also demonstrated the method on an extract of a renal carcinoma bearing an isocitrate dehydrogenase-2 (IDH2) variant that produces copious quantities of 2HG and showed that it is the D enantiomer that was exclusively produced. We also demonstrated in the same experiment that the lactate produced in the same sample was the L enantiomer.

Single cell metabolism: current and future trends.

Single cell metabolomics is an emerging and rapidly developing field that complements developments in single cell analysis by genomics and proteomics. Major goals include mapping and quantifying the metabolome in sufficient detail to provide useful information about cellular function in highly heterogeneous systems such as tissue, ultimately with spatial resolution at the individual cell level. The chemical diversity and dynamic range of metabolites poses particular challenges for detection, identification and quantification. In this review we discuss both significant technical issues of measurement and interpretation, and progress toward addressing them, with recent examples from diverse biological systems. We provide a framework for further directions aimed at improving workflow and robustness so that such analyses may become commonly applied, especially in combination with metabolic imaging and single cell transcriptomics and proteomics.

Apolipoprotein E genotype-dependent nutrigenetic effects to prebiotic inulin for modulating systemic metabolism and neuroprotection in mice via gut-brain axis.

OBJECTIVE: The goal of the study was to identify the potential nutrigenetic effects to inulin, a prebiotic fiber, in mice with different human apolipoprotein E (APOE) genetic variants. Specifically, we compared responses to inulin for the potential modulation of the systemic metabolism and neuroprotection via gut-brain axis in mice with human APOE ϵ3 and ϵ4 alleles. METHOD: We performed experiments with young mice expressing the human APOE3 (E3FAD mice and APOE4 gene (E4FAD mice). We fed mice with either inulin or control diet for 16 weeks starting from 3 months of age. We determined gut microbiome diversity and composition using16s rRNA sequencing, systemic metabolism using in vivo MRI and metabolomics, and blood-brain barrier (BBB) tight junction expression using Western blot. RESULTS: In both E3FAD and E4FAD mice, inulin altered the alpha and beta diversity of the gut microbiome, increased beneficial taxa of bacteria and elevated cecal short chain fatty acid and hippocampal scyllo-inositol. E3FAD mice had altered metabolism related to tryptophan and tyrosine, while E4FAD mice had changes in the tricarboxylic acid cycle, pentose phosphate pathway, and bile acids. Differences were found in levels of brain metabolites related to oxidative stress, and levels of Claudin-1 and Claudin-5 BBB tight junction expression. DISCUSSION: We found that inulin had many similar beneficial effects in the gut and brain for both E3FAD and E4FAD mice, which may be protective for brain functions and reduce risk for neurodegeneration. . E3FAD and E4FAD mice also had distinct responses in several metabolic pathways, suggesting an APOE-dependent nutrigenetic effects in modulating systemic metabolism and neuroprotection.

2'-Alkynyl spin-labelling is a minimally perturbing tool for DNA structural analysis.

The determination of distances between specific points in nucleic acids is essential to understanding their behaviour at the molecular level. The ability to measure distances of 2-10 nm is particularly important: deformations arising from protein binding commonly fall within this range, but the reliable measurement of such distances for a conformational ensemble remains a significant challenge. Using several techniques, we show that electron paramagnetic resonance (EPR) spectroscopy of oligonucleotides spin-labelled with triazole-appended nitroxides at the 2' position offers a robust and minimally perturbing tool for obtaining such measurements. For two nitroxides, we present results from EPR spectroscopy, X-ray crystal structures of B-form spin-labelled DNA duplexes, molecular dynamics simulations and nuclear magnetic resonance spectroscopy. These four methods are mutually supportive, and pinpoint the locations of the spin labels on the duplexes. In doing so, this work establishes 2'-alkynyl nitroxide spin-labelling as a minimally perturbing method for probing DNA conformation.

An Ion Chromatography-Ultrahigh-Resolution-MS1/Data-Independent High-Resolution MS2 Method for Stable Isotope-Resolved Metabolomics Reconstruction of Central Metabolic Networks.

The metabolome comprises a complex network of interconnecting enzyme-catalyzed reactions that involve transfers of numerous molecular subunits. Thus, the reconstruction of metabolic networks requires metabolite substructures to be tracked. Subunit tracking can be achieved by tracing stable isotopes through metabolic transformations using NMR and ultrahigh -resolution (UHR)-mass spectrometry (MS). UHR-MS1 readily resolves and counts isotopic labels in metabolites but requires tandem MS to help identify isotopic enrichment in substructures. However, it is challenging to perform chromatography-based UHR-MS1 with its long acquisition time, while acquiring MS2 data on many coeluting labeled isotopologues for each metabolite. We have developed an ion chromatography (IC)-UHR-MS1/data-independent(DI)-HR-MS2 method to trace the fate of 13C atoms from [13C6]-glucose ([13C6]-Glc) in 3D A549 spheroids in response to anticancer selenite and simultaneously 13C/15N atoms from [13C5,15N2]-glutamine ([13C5,15N2]-Gln) in 2D BEAS-2B cells in response to arsenite transformation. This method retains the complete isotopologue distributions of metabolites via UHR-MS1 while simultaneously acquiring substructure label information via DI-MS2. These details in metabolite labeling patterns greatly facilitate rigorous reconstruction of multiple, intersecting metabolic pathways of central metabolism, which are illustrated here for the purine/pyrimidine nucleotide biosynthesis. The pathways reconstructed based on subunit-level isotopologue analysis further reveal specific enzyme-catalyzed reactions that are impacted by selenite or arsenite treatments.

NMR Analysis of Carboxylate Isotopomers of 13C-Metabolites by Chemoselective Derivatization with 15N-Cholamine.

A substantial fraction of common metabolites contains carboxyl functional groups. Their 13C isotopomer analysis by nuclear magnetic resonance (NMR) is hampered by the low sensitivity of the 13C nucleus, the slow longitudinal relaxation for the lack of an attached proton, and the relatively low chemical shift dispersion of carboxylates. Chemoselective (CS) derivatization is a means of tagging compounds in a complex mixture via a specific functional group. 15N1-cholamine has been shown to be a useful CS agent for carboxylates, producing a peptide bond that can be detected via 15N-attached H with high sensitivity in heteronuclear single quantum coherence experiments. Here, we report an improved method of derivatization and show how 13C-enrichment at the carboxylate and/or the adjacent carbon can be determined via one- and two-bond coupling of the carbons adjacent to the cholamine 15N atom in the derivatives. We have applied this method for the determination of 13C isotopomer distribution in the extracts of A549 cell culture and liver tissue from a patient-derived xenograft mouse.

JAK2-mutant hematopoietic cells display metabolic alterations that can be targeted to treat myeloproliferative neoplasms.

Increased energy requirement and metabolic reprogramming are hallmarks of cancer cells. We show that metabolic alterations in hematopoietic cells are fundamental to the pathogenesis of mutant JAK2-driven myeloproliferative neoplasms (MPNs). We found that expression of mutant JAK2 augmented and subverted metabolic activity of MPN cells, resulting in systemic metabolic changes in vivo, including hypoglycemia, adipose tissue atrophy, and early mortality. Hypoglycemia in MPN mouse models correlated with hyperactive erythropoiesis and was due to a combination of elevated glycolysis and increased oxidative phosphorylation. Modulating nutrient supply through high-fat diet improved survival, whereas high-glucose diet augmented the MPN phenotype. Transcriptomic and metabolomic analyses identified numerous metabolic nodes in JAK2-mutant hematopoietic stem and progenitor cells that were altered in comparison with wild-type controls. We studied the consequences of elevated levels of Pfkfb3, a key regulatory enzyme of glycolysis, and found that pharmacological inhibition of Pfkfb3 with the small molecule 3PO reversed hypoglycemia and reduced hematopoietic manifestations of MPNs. These effects were additive with the JAK1/2 inhibitor ruxolitinib in vivo and in vitro. Inhibition of glycolysis by 3PO altered the redox homeostasis, leading to accumulation of reactive oxygen species and augmented apoptosis rate. Our findings reveal the contribution of metabolic alterations to the pathogenesis of MPNs and suggest that metabolic dependencies of mutant cells represent vulnerabilities that can be targeted for treating MPNs.

De novo synthesis of serine and glycine fuels purine nucleotide biosynthesis in human lung cancer tissues.

Nucleotide synthesis is essential to proliferating cells, but the preferred precursors for de novo biosynthesis are not defined in human cancer tissues. We have employed multiplexed stable isotope-resolved metabolomics to track the metabolism of [13C6]glucose, D2-glycine, [13C2]glycine, and D3-serine into purine nucleotides in freshly resected cancerous and matched noncancerous lung tissues from nonsmall cell lung cancer (NSCLC) patients, and we compared the metabolism with established NSCLC PC9 and A549 cell lines in vitro Surprisingly, [13C6]glucose was the best carbon source for purine synthesis in human NSCLC tissues, in contrast to the noncancerous lung tissues from the same patient, which showed lower mitotic indices and MYC expression. We also observed that D3-Ser was preferentially incorporated into purine rings over D2-glycine in both tissues and cell lines. MYC suppression attenuated [13C6]glucose, D3-serine, and [13C2]glycine incorporation into purines and reduced proliferation in PC9 but not in A549 cells. Using detailed kinetic modeling, we showed that the preferred use of glucose as a carbon source for purine ring synthesis in NSCLC tissues involves cytoplasmic activation/compartmentation of the glucose-to-serine pathway and enhanced reversed one-carbon fluxes that attenuate exogenous serine incorporation into purines. Our findings also indicate that the substrate for de novo nucleotide synthesis differs profoundly between cancer cell lines and fresh human lung cancer tissues; the latter preferred glucose to exogenous serine or glycine but not the former. This distinction in substrate utilization in purine synthesis in human cancer tissues should be considered when targeting one-carbon metabolism for cancer therapy.

Applications of Chromatography-Ultra High-Resolution MS for Stable Isotope-Resolved Metabolomics (SIRM) Reconstruction of Metabolic Networks.

Metabolism is a complex network of compartmentalized and coupled chemical reactions, which often involve transfers of substructures of biomolecules, thus requiring metabolite substructures to be tracked. Stable isotope resolved metabolomics (SIRM) enables pathways reconstruction, even among chemically identical metabolites, by tracking the provenance of stable isotope-labeled substructures using NMR and ultrahigh resolution (UHR) MS. The latter can resolve and count isotopic labels in metabolites and can identify isotopic enrichment in substructures when operated in tandem MS mode. However, MS2 is difficult to implement with chromatography-based UHR-MS due to lengthy MS1 acquisition time that is required to obtain the molecular isotopologue count, which is further exacerbated by the numerous isotopologue source ions to fragment. We review here recent developments in tandem MS applications of SIRM to obtain more detailed information about isotopologue distributions in metabolites and their substructures.

Epigenetic Modifications of Cytosine: Biophysical Properties, Regulation, and Function in Mammalian DNA.

To decode the function and molecular recognition of several recently discovered cytosine derivatives in the human genome - 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxylcytosine - a detailed understanding of their effects on the structural, chemical, and biophysical properties of DNA is essential. Here, we review recent literature in this area, with particular emphasis on features that have been proposed to enable the specific recognition of modified cytosine bases by DNA-binding proteins. These include electronic factors, modulation of base-pair stability, flexibility, and radical changes in duplex conformation. We explore these proposals and assess whether or not they are supported by current biophysical data. This analysis is focused primarily on the properties of epigenetically modified DNA itself, which provides a basis for discussion of the mechanisms of recognition by different proteins.

UCP2 Overexpression Redirects Glucose into Anabolic Metabolic Pathways.

Uncoupling protein 2 (UCP2) is often upregulated in cancer cells. The UCP2 upregulation is positively correlated with enhanced proliferation, tumorigenesis, and metabolic alterations, thus suggesting that UCP2 upregulation can play a key role in sensing metabolic changes to promote tumorigenesis. To determine the global metabolic impact of UCP2 upregulation, 13 C6 glucose as a source molecule is used to "trace" the metabolic fate of carbon atoms derived from glucose. UCP2 overexpression in skin epidermal cells enhances the incorporation of 13 C label to pyruvate, tricarboxylic acid cycle intermediates, nucleotides, and amino acids, suggesting that UCP2 upregulation reprograms cellular metabolism toward macromolecule synthesis. To the best of our knowledge, this is the first study to bring to light the overall metabolic differences caused by UCP2 upregulation.

Inhibition of Anaplerotic Glutaminolysis Underlies Selenite Toxicity in Human Lung Cancer.

Large clinical trials and model systems studies suggest that the chemical form of selenium dictates chemopreventive and chemotherapeutic efficacy. Selenite induces excess ROS production, which mediates autophagy and eventual cell death in non-small cell lung cancer adenocarcinoma A549 cells. As the mechanisms underlying these phenotypic effects are unclear, the clinical relevance of selenite for cancer therapy remains to be determined. The authors' previous stable isotope-resolved metabolomics and gene expression analysis showed that selenite disrupts glycolysis, the Krebs cycle, and polyamine metabolism in A549 cells, potentially through perturbed glutaminolysis, a vital anaplerotic process for proliferation of many cancer cells. Herein, the role of the glutaminolytic enzyme glutaminase 1 (GLS1) in selenite's toxicity in A549 cells and in patient-derived lung cancer tissues is investigated. Using [13 C6 ]-glucose and [13 C5 ,15 N2 ]-glutamine tracers, selenite's action on metabolic networks is determined. Selenite inhibits glutaminolysis and glutathione synthesis by suppressing GLS1 expression, and blocks the Krebs cycle, but transiently activates pyruvate carboxylase activity. Glutamate supplementation partially rescues these anti-proliferative and oxidative stress activities. Similar metabolic perturbations and necrosis are observed in selenite-treated human patients' cancerous lung tissues ex vivo. The results support the hypothesis that GLS1 suppression mediates part of the anti-cancer activity of selenite both in vitro and ex vivo.