1 Quantitative Determinants of Aerobic Glycolysis Identify Flux Through the Enzyme 1 GAPDH As a Limiting Step 2 3 Alexander a Sh

Total Page:16

File Type:pdf, Size:1020Kb

1 Quantitative Determinants of Aerobic Glycolysis Identify Flux Through the Enzyme 1 GAPDH As a Limiting Step 2 3 Alexander a Sh 1 Quantitative determinants of aerobic glycolysis identify flux through the enzyme 2 GAPDH as a limiting step 3 4 Alexander A Shestov1*, Xiaojing Liu1*, Zheng Ser1, Ahmad Cluntun2, Yin P Hung3, Lei 5 Huang4, Dongsung Kim2, Anne Le5, Gary Yellen3, John G Albeck6, Jason W 6 Locasale1,2,4** 7 8 1Division of Nutritional Sciences, Cornell University, Ithaca NY 14850 9 2Field of Biochemistry and Molecular Cell Biology, Department of Molecular Biology 10 and Genetics Cornell University, Ithaca NY 14850 11 3Department of Neurobiology, Harvard Medical School, Boston MA 02115 12 4Field of Computational Biology, Department of Biological Statistics and Computational 13 Biology, Cornell University, Ithaca NY 14850 14 5Department of Pathology and Oncology, Johns Hopkins University School of Medicine, 15 Baltimore, MD 21231 16 6Department of Molecular Cell Biology, University of California Davis, Davis CA 95616 17 * Equal contribution 18 ** Correspondence: [email protected] 19 20 Abstract 21 Aerobic glycolysis or the Warburg Effect (WE) is characterized by the increased 22 metabolism of glucose to lactate. It remains unknown what quantitative changes to the 23 activity of metabolism are necessary and sufficient for this phenotype. We developed a 24 computational model of glycolysis and an integrated analysis using metabolic control 25 analysis (MCA), metabolomics data, and statistical simulations. We identified and 26 confirmed a novel mode of regulation specific to aerobic glycolysis where flux through 27 GAPDH, the enzyme separating lower and upper glycolysis, is the rate-limiting step in 28 the pathway and the levels of fructose (1,6) bisphosphate (FBP), are predictive of the rate 29 and control points in glycolysis. Strikingly, negative flux control was found and 30 confirmed for several steps thought to be rate-limiting in glycolysis. Together these 31 findings enumerate the biochemical determinants of the WE, and suggest strategies for 32 identifying the contexts in which agents that target glycolysis might be most effective. 33 34 Impact Statement 35 A computational model of the Warburg Effect is developed that reveals new mechanisms 36 for the regulation of glycolysis. These findings could lead to the definition of biomarkers 37 that determine the responses to targeting glycolysis in malignancy. 38 39 40 41 42 43 44 45 46 1 47 48 49 50 Introduction 51 52 Proliferating cells increase their glucose consumption and secrete lactate as 53 opposed to completely oxidizing the glucose in the mitochondria(Warburg, Wind et al. 54 1927) and is known as aerobic glycolysis or the Warburg Effect. Currently, this altered 55 metabolism is exploited for diagnostics and is subject to multiple drug development 56 efforts(Koppenol, Bounds et al. 2011; Vander Heiden 2011; Hamanaka and Chandel 57 2012). Numerous studies have identified genes such as KRAS, PIK3CA, and cMYC and 58 microenvironments such as hypoxia and hypoglycemia that promote aerobic glycolysis 59 but a complete understanding of the necessary and sufficient biochemical alterations 60 associated with this phenotype is unknown. Furthermore successful translation for 61 biomedical applications is limited by understanding the contexts in which therapies that 62 target glycolysis might be effective. 63 Computational modeling has a successful history in the study of 64 metabolism(Rapoport, Heinrich et al. 1976; Fell 1992; Schilling, Schuster et al. 1999; 65 Cascante, Boros et al. 2002). Genome-scale stoichiometric models of metabolism have 66 been developed to study the effects of drug targets in human metabolism and have had 67 success in predicting the WE(Molenaar, van Berlo et al. 2009; Vazquez, Liu et al. 2010; 68 Folger, Jerby et al. 2011; Shlomi, Benyamini et al. 2011). However a comprehensive 69 quantitative understanding of the WE requires knowledge of enzyme activities and 70 metabolic control. 71 Therefore, we collected and integrated multiple forms of data into a modeling 72 framework involving flux balances of glycolysis, detailed chemical kinetics based on 2 73 reaction mechanisms and parameters measured, physico-chemical constraints from 74 thermodynamics and mass conservation, metabolic control analysis, and Monte Carlo 75 sampling of parameter space. We next use mass spectrometry and isotope tracing to 76 probe concentrations and fluxes through the pathway and their responses to several 77 perturbations. Together we elaborate the determinants of aerobic glycolysis and identify 78 and confirm novel points of regulation in glycolysis that have remained unidentified for 79 over 50 years since the discovery of the pathway. 80 81 Results 82 Biochemical kinetic model of aerobic glycolysis 83 We investigated the kinetics of the glycolytic pathway from glucose uptake to 84 oxidation of pyruvate in the mitochondria or export of lactate out of the cytosol. We 85 modeled each step of the pathway according to enzymatic mechanism and known modes 86 of allosteric control resulting in a set of differential equations (Fig 1A, Materials and 87 Methods, Supplementary File 1). While it is not possible to model every possible 88 interaction explicitly, the aim is to capture enough of the pathway so that a large range of 89 experimentally realized measurements can be obtained and relationships between 90 variables can be observed. 91 Since glycolysis is the most extensively studied biochemical pathway, there is a 92 wealth of information on the kinetic parameters and enzyme expression that govern the 93 equations. Nevertheless, it is also not impossible to capture cellular physiology in any 94 biochemical model with single values of kinetic parameters(Daniels, Chen et al. 2008). 95 This difficulty arises from tremendous amount of heterogeneity within cells at multiple 3 96 levels. The origins of this heterogeneity vary from genetic variation observed across 97 cancer types, tumor types, differences in signaling mechanisms that affect post- 98 translational modifications in each cell, and the differences in microenvironmental 99 pressures (e.g. the oxygen availability) that each cell within a given tumor experiences, as 100 well as the inherent cell to cell variation common to all cells(Marusyk, Almendro et al. 101 2012). Therefore we developed an integrated algorithm to evaluate the statistics of the 102 kinetics of glycolysis that accounts for the possible variation within metabolism (Fig 1B, 103 Supplementary Information). 104 First the model is constrained using mass conservation constraints that conserve 105 the balance of glucose, redox state, and energy status. Next thermodynamic constraints 106 are used to constrain fluxes according to the free energy of the reactions determined by 107 Haldane relationships. These physical constraints are combined with the kinetic 108 mechanisms that define each step of glycolysis and the chemical reactions involving the 109 redox-associated metabolites NAD+, and NADH and energy-associated metabolites ATP, 110 ADP, and AMP. Next, the model is constrained to expression data so that protein 111 concentrations are subsumed in the Vmax values and chosen from typical concentrations 112 in cancer cells, measured kinetic parameters, and measured concentrations of nutrients 113 such as glucose, oxygen, and total intracellular adeno-nucleotide concentration. 114 At this stage, the model is subjected to a thorough statistical analysis. A Monte- 115 Carlo simulation is conducted for which the parameters within the model are randomized 116 and resulting differential equations solved. The distributions for each parameter are 117 chosen to capture observed ranges of variation. After each simulation, numerical stability 118 and the thermodynamics are assessed and all simulations that are unstable or appear 4 119 thermodynamically infeasible (i.e. a positive net flux through glycolysis is required) are 120 rejected. In each simulation, concentrations, fluxes, metabolic control coefficients, and 121 thermodynamic quantities are computed and recorded. This statistical analysis explores 122 the space of glycolysis in the context of the Warburg Effect. By assessing the statistics, 123 inferences can be made on the determinants of the Warburg Effects and its context 124 dependence to pharmacological intervention and nutrient environment. We first utilized 125 a recently developed a NADH/NAD+ fluorescent reporter(Hung, Albeck et al. 2011) and 126 measured the ratio of NADH/NAD+ across a population of cells where some experience 127 hypoxia and others glucose deprivation (Fig 1C). This variation in nutrient availability 128 across individual cells occurs in epithelial cell cultures due to differences in diffusion and 129 available cell surface area(Sheta, Trout et al. 2001). The distribution is plotted for three 130 media concentrations ranging from extreme hypoglycemia to the hyperglycemic 131 conditions typically used in cell culture (0.5 mM, 5.5mM, and 25mM). 132 The distribution of resulting fluxes (Fig 1D) was first found consistent with 133 known measurements. An analysis of the simulation revealed a distribution of NAD+ 134 levels and NADH/NAD+ redox potential consistent with those observed experimentally 135 (Figs 1E,F). ATP levels peaked around 3mM with little variation (Fig 1G) and the 136 ATP/ADP energy state is observed to be bimodal (Fig 1H). Concentrations along the 137 glycolytic pathway varied with means similar to those measured in normal tissues with 138 fructose-1,6-bisphosphate (FBP) being most variable (Fig 1I). Together, these findings 139 indicate
Recommended publications
  • In Vivo Metabolism of [1,6-13C2]Glucose Reveals Distinct
    H OH metabolites OH Article 13 In Vivo Metabolism of [1,6- C2]Glucose Reveals Distinct Neuroenergetic Functionality between Mouse Hippocampus and Hypothalamus Antoine Cherix 1,†, Rajesh Sonti 1,‡, Bernard Lanz 1 and Hongxia Lei 2,3,*,§ 1 Laboratory of Functional and Metabolic Imaging (LIFMET), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; [email protected] (A.C.); [email protected] (R.S.); bernard.lanz@epfl.ch (B.L.) 2 Animal Imaging and Technology (AIT), Center for Biomedical Imaging (CIBM), Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland 3 Faculty of Medicine, University of Geneva, CH-1206 Geneva, Switzerland * Correspondence: [email protected] † Current address: Wellcome Centre for Integrative Neuroimaging (WIN), Oxford Centre for Functional MRI of the Brain (FMRIB), Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford OX3 9DU, UK. ‡ Current address: Department of Pharmaceutical analysis, National Institute of Pharmaceutical Education and Research (NIPER) Hyderabad- NH 65, Hyderabad 500037, India. § Current address: Wuhan United Imaging Life Science Instruments Ltd., Wuhan 430206, China. Abstract: Glucose is a major energy fuel for the brain, however, less is known about specificities of its metabolism in distinct cerebral areas. Here we examined the regional differences in glucose utilization between the hypothalamus and hippocampus using in vivo indirect 13C magnetic res- 1 13 13 onance spectroscopy ( H-[ C]-MRS) upon infusion of [1,6- C2]glucose. Using a metabolic flux analysis with a 1-compartment mathematical model of brain metabolism, we report that compared to hippocampus, hypothalamus shows higher levels of aerobic glycolysis associated with a marked Citation: Cherix, A.; Sonti, R.; Lanz, gamma-aminobutyric acid-ergic (GABAergic) and astrocytic metabolic dependence.
    [Show full text]
  • Quantitative Analysis of Amino Acid Metabolism in Liver Cancer Links Glutamate Excretion to Nucleotide Synthesis
    Quantitative analysis of amino acid metabolism in liver cancer links glutamate excretion to nucleotide synthesis Avlant Nilssona,1, Jurgen R. Haanstrab,1, Martin Engqvista, Albert Gerdingc,d, Barbara M. Bakkerb,c, Ursula Klingmüllere, Bas Teusinkb, and Jens Nielsena,f,2 aDepartment of Biology and Biological Engineering, Chalmers University of Technology, SE41296 Gothenburg, Sweden; bSystems Biology Lab, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Vrije Universiteit Amsterdam, NL1081HZ Amsterdam, The Netherlands; cLaboratory of Pediatrics, Systems Medicine of Metabolism and Signaling, University of Groningen, University Medical Center Groningen, NL-9713AV Groningen, The Netherlands dDepartment of Laboratory Medicine, University of Groningen, University Medical Center Groningen, NL-9713AV Groningen, The Netherlands; eDivision of Systems Biology and Signal Transduction, German Cancer Research Center, D-69120 Heidelberg, Germany; and fNovo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Kongens Lyngby, DK2800, Denmark Contributed by Jens Nielsen, March 9, 2020 (sent for review November 4, 2019; reviewed by Eytan Ruppin and Matthew G. Vander Heiden) Many cancer cells consume glutamine at high rates; counterintu- the excretion of lactate, glutamate, alanine, and glycine (5). How- itively, they simultaneously excrete glutamate, the first interme- ever, a full genome-wide modeling approach is required to expand diate in glutamine metabolism. Glutamine consumption has been our knowledge of how metabolism functions beyond the canonical linked to replenishment of tricarboxylic acid cycle (TCA) interme- pathways and, in particular, to understand the intricate effects of diates and synthesis of adenosine triphosphate (ATP), but the metabolic compartmentalization and the interplay between ex- reason for glutamate excretion is unclear. Here, we dynamically change fluxes and synthesis of biomass.
    [Show full text]
  • Nicotinamide Adenine Dinucleotide Is Transported Into Mammalian
    RESEARCH ARTICLE Nicotinamide adenine dinucleotide is transported into mammalian mitochondria Antonio Davila1,2†, Ling Liu3†, Karthikeyani Chellappa1, Philip Redpath4, Eiko Nakamaru-Ogiso5, Lauren M Paolella1, Zhigang Zhang6, Marie E Migaud4,7, Joshua D Rabinowitz3, Joseph A Baur1* 1Department of Physiology, Institute for Diabetes, Obesity, and Metabolism, Perelman School of Medicine, University of Pennsylvania, Philadelphia, United States; 2PARC, Perelman School of Medicine, University of Pennsylvania, Philadelphia, United States; 3Lewis-Sigler Institute for Integrative Genomics, Department of Chemistry, Princeton University, Princeton, United States; 4School of Pharmacy, Queen’s University Belfast, Belfast, United Kingdom; 5Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, United States; 6College of Veterinary Medicine, Northeast Agricultural University, Harbin, China; 7Mitchell Cancer Institute, University of South Alabama, Mobile, United States Abstract Mitochondrial NAD levels influence fuel selection, circadian rhythms, and cell survival under stress. It has alternately been argued that NAD in mammalian mitochondria arises from import of cytosolic nicotinamide (NAM), nicotinamide mononucleotide (NMN), or NAD itself. We provide evidence that murine and human mitochondria take up intact NAD. Isolated mitochondria preparations cannot make NAD from NAM, and while NAD is synthesized from NMN, it does not localize to the mitochondrial matrix or effectively support oxidative phosphorylation. Treating cells *For correspondence: with nicotinamide riboside that is isotopically labeled on the nicotinamide and ribose moieties [email protected] results in the appearance of doubly labeled NAD within mitochondria. Analogous experiments with †These authors contributed doubly labeled nicotinic acid riboside (labeling cytosolic NAD without labeling NMN) demonstrate equally to this work that NAD(H) is the imported species.
    [Show full text]
  • Metabolic Flux Analysis—Linking Isotope Labeling and Metabolic Fluxes
    H OH metabolites OH Review Metabolic Flux Analysis—Linking Isotope Labeling and Metabolic Fluxes Yujue Wang 1,2, Fredric E. Wondisford 1, Chi Song 3, Teng Zhang 4 and Xiaoyang Su 1,2,* 1 Department of Medicine, Rutgers-Robert Wood Johnson Medical School, New Brunswick, NJ 08901, USA; [email protected] (Y.W.); [email protected] (F.E.W.) 2 Metabolomics Shared Resource, Rutgers Cancer Institute of New Jersey, New Brunswick, NJ 08903, USA 3 Division of Biostatistics, College of Public Health, The Ohio State University, Columbus, OH 43210, USA; [email protected] 4 Department of Mathematics, University of Central Florida, Orlando, FL 32816, USA; [email protected] * Correspondence: [email protected]; Tel.: +1-732-235-5447 Received: 14 September 2020; Accepted: 4 November 2020; Published: 6 November 2020 Abstract: Metabolic flux analysis (MFA) is an increasingly important tool to study metabolism quantitatively. Unlike the concentrations of metabolites, the fluxes, which are the rates at which intracellular metabolites interconvert, are not directly measurable. MFA uses stable isotope labeled tracers to reveal information related to the fluxes. The conceptual idea of MFA is that in tracer experiments the isotope labeling patterns of intracellular metabolites are determined by the fluxes, therefore by measuring the labeling patterns we can infer the fluxes in the network. In this review, we will discuss the basic concept of MFA using a simplified upper glycolysis network as an example. We will show how the fluxes are reflected in the isotope labeling patterns. The central idea we wish to deliver is that under metabolic and isotopic steady-state the labeling pattern of a metabolite is the flux-weighted average of the substrates’ labeling patterns.
    [Show full text]
  • Prediction of Metabolic Fluxes by Incorporating Genomic Context And
    Prediction of metabolic fluxes by incorporating genomic context and flux-converging pattern analyses Jong Myoung Parka,b, Tae Yong Kima,b, and Sang Yup Leea,b,c,1 aMetabolic and Biomolecular Engineering National Research Laboratory, Department of Chemical and Biomolecular Engineering (BK21 Program), KAIST, Daejeon 305-701, Republic of Korea; bCenter for Systems and Synthetic Biotechnology, Institute for the BioCentury, KAIST, Daejeon 305-701, Republic of Korea; and cDepartment of Bio and Brain Engineering, BioProcess Engineering Research Center and Bioinformatics Research Center, KAIST, Daejeon 305-701, Republic of Korea Edited* by Bernhard Ø. Palsson, University of California at San Diego, La Jolla, CA, and accepted by the Editorial Board July 6, 2010 (received for review March 26, 2010) Flux balance analysis (FBA) of a genome-scale metabolic model However, this approach still works with a flux solution space that allows calculation of intracellular fluxes by optimizing an objective is too large compared with the biologically feasible flux solution function, such as maximization of cell growth, under given con- space of a real organism (15–17). Thus, additional information straints, and has found numerous applications in the field of systems and procedures, such as providing more constraints, are needed to biology and biotechnology. Due to the underdetermined nature of narrow down the set of candidate solutions. the system, however, it has limitations such as inaccurate prediction One example of a suitable constraint that can be applied is the of fluxes and existence of multiple solutions for an optimal objective experimentally measured metabolic flux data obtained from 13C- value. Here, we report a strategy for accurate prediction of meta- based flux analysis and fermentation data (18–21).
    [Show full text]
  • Glucose Flux in Relation to Energy Expenditure in Malnourished Patients with and Without Cancer During Periods of Fasting and Feeding1
    [CANCER RESEARCH 44, 1718-1724, April 1984] Glucose Flux in Relation to Energy Expenditure in Malnourished Patients with and without Cancer during Periods of Fasting and Feeding1 Elisabeth Edén,Staffar) Edström,Klas Bennegárd, Tore Scherstén,and Kent Lundholm2 Surgical Metabolic Research Laboratory, Department of Surgery I [T. S., K. L], Anaesthasiology l-ll [E. E., K. B.] and Otolaryngology [S. £.],Sahlgrenska Hospital, University of Gothenburg, Gothenburg, Sweden ABSTRACT fuel economy in progressive neoplasia (12). Our recent evidence suggests that elevated lactate turnover in cancer patients may Glucose dynamics, energy metabolism, and nitrogen balance be counteracted by depressed recycling of carbon in the glucose- were studied in eight malnourished cancer patients and seven alanine cycle (2), supported by the findings that a decreased malnourished patients without cancer. Glucose flux was mea alanine release occurred in combination with a diminished glu sured by single injection of [6-3H]glucose and [U-14C]glucose. cose uptake across the leg in cancer patients after an overnight Energy expenditure was measured by indirect calorimetry. Each fast (2), and that intrahepatic cycling of glucose carbon was patient was studied after an overnight fast and during constant decreased in sarcoma-bearing mice (21). Thus, the evidence gastric infusion of a formula diet. points to the possibility that net recycling of glucose in cancer Cancer patients had elevated glucose flux when fasting, cor cachexia may not influence the overall fuel economy, even if the responding to 42% of their spontaneous daily intake of glucose. Cori cycle activity is highly elevated (16). Instead, it is rather the At least one-half of the elevated flux in cancer patients compared sum of whole-body gluconeogenesis and glycogenolysis, mea with controls was due to increased recycling of glucose carbon sured as glucose turnover, that may be of quantitative impor after an overnight fast.
    [Show full text]
  • 13C-Metabolic Flux Analysis: an Accurate Approach to Demystify Microbial Metabolism for Biochemical Production
    Review 13C-Metabolic Flux Analysis: An Accurate Approach to Demystify Microbial Metabolism for Biochemical Production Weihua Guo, Jiayuan Sheng and Xueyang Feng * Received: 1 October 2015; Accepted: 18 December 2015; Published: 25 December 2015 Academic Editors: Mark Blenner and Michael D. Lynch Department of Biological Systems Engineering, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA; [email protected] (W.G.); [email protected] (J.S.) * Correspondence: [email protected]; Tel.: +1-540-231-2974 Abstract: Metabolic engineering of various industrial microorganisms to produce chemicals, fuels, and drugs has raised interest since it is environmentally friendly, sustainable, and independent of nonrenewable resources. However, microbial metabolism is so complex that only a few metabolic engineering efforts have been able to achieve a satisfactory yield, titer or productivity of the target chemicals for industrial commercialization. In order to overcome this challenge, 13C Metabolic Flux Analysis (13C-MFA) has been continuously developed and widely applied to rigorously investigate cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, many 13C-MFA studies have been performed in academic labs and biotechnology industries to pinpoint key issues related to microbe-based chemical production. Insightful information about the metabolic rewiring has been provided to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this review, we will introduce the basics of 13C-MFA and illustrate how 13C-MFA has been applied via integration with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production for various host microorganisms Keywords: Bottleneck; isotope; cofactor imbalance; cell metabolism; synthetic biology; biofuels 1.
    [Show full text]
  • Dynamic Flux Balance Analysis to Evaluate the Strain Production Performance on Shikimic Acid Production in Escherichia Coli
    H OH metabolites OH Article Dynamic Flux Balance Analysis to Evaluate the Strain Production Performance on Shikimic Acid Production in Escherichia coli Yuki Kuriya 1 and Michihiro Araki 1,2,3,* 1 Graduate School of Medicine, Kyoto University, 54 ShogoinKawahara-cho, Sakyo-ku, Kyoto 606-8507, Japan; [email protected] 2 National Institutes of Biomedical Innovation, Health and Nutrition, 1-23-1 Toyama, Shinjuku-ku, Tokyo 162-8636, Japan 3 Graduate School of Science, Technology and Innovation, Kobe University, 1-1 Rokkodai, Nada, Kobe, Hyogo 657-8501, Japan * Correspondence: [email protected]; Tel.: +81-3-3203-5721 Received: 17 April 2020; Accepted: 12 May 2020; Published: 15 May 2020 Abstract: Flux balance analysis (FBA) is used to improve the microbial production of useful compounds. However, a large gap often exists between the FBA solution and the experimental yield, because of growth and byproducts. FBA has been extended to dynamic FBA (dFBA), which is applicable to time-varying processes, such as batch or fed-batch cultures, and has significantly contributed to metabolic and cultural engineering applications. On the other hand, the performance of the experimental strains has not been fully evaluated. In this study, we applied dFBA to the production of shikimic acid from glucose in Escherichia coli, to evaluate the production performance of the strain as a case study. The experimental data of glucose consumption and cell growth were used as FBA constraints. Bi-level FBA optimization with maximized growth and shikimic acid production were the objective functions. Results suggest that the shikimic acid concentration in the high-shikimic-acid-producing strain constructed in the experiment reached up to 84% of the maximum value by simulation.
    [Show full text]
  • A Genetically Encoded Tool for Manipulation of NADP+/NADPH in Living Cells
    HHS Public Access Author manuscript Author ManuscriptAuthor Manuscript Author Nat Chem Manuscript Author Biol. Author Manuscript Author manuscript; available in PMC 2018 February 07. Published in final edited form as: Nat Chem Biol. 2017 October ; 13(10): 1088–1095. doi:10.1038/nchembio.2454. A genetically encoded tool for manipulation of NADP+/NADPH in living cells Valentin Cracan1,2,3,†, Denis V. Titov1,2,3,†, Hongying Shen1,2,3, Zenon Grabarek1,*, and Vamsi K. Mootha1,2,3,* 1Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA 2Department of Systems Biology, Harvard Medical School, Boston, MA, USA 3Broad Institute, Cambridge, MA, USA Abstract NADH and NADPH are redox coenzymes broadly required for energy metabolism, biosynthesis and detoxification. Despite detailed knowledge of specific enzymes and pathways that utilize these coenzymes, a holistic understanding of the regulation and compartmentalization of NADH and NADPH-dependent pathways is lacking, in part because of a lack of tools with which to investigate them in living cells. We previously reported the use of the naturally occurring Lactobacillus brevis H2O-forming NADH oxidase (LbNOX) as a genetic tool for manipulation of the NAD+/NADH ratio in human cells. Here we present TPNOX (triphosphopyridine nucleotide oxidase), a rationally designed and engineered mutant of LbNOX that is strictly specific towards NADPH. We characterize the effects of TPNOX expression on cellular metabolism and use it in combination with LbNOX to show how the redox states of mitochondrial NADPH and NADH pools are connected. Graphical Abstract Users may view, print, copy, and download text and data-mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use: http://www.nature.com/authors/editorial_policies/license.html#terms *Corresponding authors: Vamsi K.
    [Show full text]
  • Capacity for Instantaneous Catabolism of Preferred and Non-Preferred
    www.nature.com/scientificreports OPEN Capacity for instantaneous catabolism of preferred and non-preferred carbon sources in Received: 6 April 2018 Accepted: 26 July 2018 Escherichia coli and Bacillus subtilis Published: xx xx xxxx Marieke F. Bufng1,2, Hannes Link3, Dimitris Christodoulou 1,2 & Uwe Sauer1 Making the right choice for nutrient consumption in an ever-changing environment is a key factor for evolutionary success of bacteria. Here we investigate the regulatory mechanisms that enable dynamic adaptation between non-preferred and preferred carbon sources for the model Gram- negative and -positive species Escherichia coli and Bacillus subtilis, respectively. We focus on the ability for instantaneous catabolism of a gluconeogenic carbon source upon growth on a glycolytic carbon source and vice versa. By following isotopic tracer dynamics on a 1–2 minute scale, we show that fux reversal from the preferred glucose to non-preferred pyruvate as the sole carbon source is primarily transcriptionally regulated. In the opposite direction, however, E. coli can reverse its fux instantaneously by means of allosteric regulation, whereas in B. subtilis this fux reversal is transcriptionally regulated. Upon removal of transcriptional regulation, B. subtilis assumes the ability of instantaneous glucose catabolism. Using an approach that combines quantitative metabolomics and kinetic modelling, we then identify the additionally necessary key metabolite-enzyme interactions that implement the instantaneous fux reversal in the transcriptionally deregulated B. subtilis, and validate the most relevant allosteric interactions. Evolutionary adaptation enables bacterial growth and survival in a large range of habitats, as is evident from the ability to exploit diverse and frequently changing energy and nutrient sources1.
    [Show full text]
  • Cofactor Engineering for Advancing Chemical Biotechnology. Yipeng
    Cofactor engineering for advancing chemical biotechnology. Yipeng Wanga, Ka-Yiu Sanb,c and George N. Bennetta a. Department of Biochemistry and Cell Biology, Rice University, Houston, TX USA 77005 b. Department of Bioengineering, Rice University, Houston, TX USA 77005 c. Department of Chemical and Biomolecular Engineering, Rice University, Houston, USA TX 77251 Corresponding author George N. Bennett Telephone:713 348 4920 e-mail: [email protected] fax: 713 348 5154 Abstract Cofactors provide redox carriers for biosynthetic reactions, catabolic reactions and act as important agents in transfer of energy for the cell. Recent advances in manipulating cofactors include culture conditions or additive alterations, genetic modification of host pathways for increased availability of desired cofactor, changes in enzyme cofactor specificity, and introduction of novel redox partners to form effective circuits for biochemical processes and biocatalysts. Genetic strategies to employ ferredoxin, NADH and NADPH most effectively in natural or novel pathways have improved yield and efficiency of large-scale processes for fuels and chemicals and have been demonstrated with a variety of microbial organisms. Introduction Cofactors provide redox carriers for biosynthetic reactions, catabolic reactions and act as important agents in transfer of energy for the cell. Some examples illustrating the requirement for cofactor balance and availability include: the conversion of biomass feedstocks containing xylose to ethanol where the formation of xylitol is a problem [1-7]; as a driving force for more effective production of reduced compounds such as biofuels [8]; in using cytochrome P450s in specific oxidation reactions where the recycling of active enzyme is required [9-11]; and the production of chiral pharmaceutical intermediates where specific reductions require a certain cofactor [12,13].
    [Show full text]
  • Flux Balance Analysis Flux Balance Analysis (FBA) Is a Computational Method for Studying Fluxes Through Metabolic Networks at Steady State
    Flux balance analysis Flux balance analysis (FBA) is a computational method for studying fluxes through metabolic networks at steady state. Here, steady state means that metabolic fluxes (reaction rates) and ​ ​ intracellular metabolite concentration are both constant over time. Therefore, the total rates of production and consumption of any intracellular metabolite must be equal (otherwise the ​ ​ metabolite concentration would change). We will see how this “balancing” of fluxes help us analyze metabolic networks. For each reaction acting on metabolite in a network, we assign a negative stoichiometry coefficient s < 0 if the metabolite is a substrate, and a positive s > 0 if the metabolic is a product ​ ​ ​ ​ of the reaction. For example, for the adenylate kinase reaction AMP + ATP ←→ 2 ADP the stoichiometry coefficient is ­1 for AMP, ­1 for ATP, and +2 for ADP. To deal with metabolic ​ ​ networks consisting of many reactions and metabolites, we collect all these stoichiometry coefficients in a stoichiometry matrix S where metabolites are rows and reactions are columns, ​ ​ ​ ​ and write Smr for the coefficient of metabolite m, reaction r. Below is an example network and ​ ​ ​ ​ ​ ​ ​ corresponding stoichiometry matrix for the urea cycle. For example, the Ass1 has coefficients S33 = ­1 for Citrulline (substrate) and S43 = +1 for ArgSucc (product). Note that uptake and ​ ​ ​ ​ ​ release fluxes (arrows in/out of the “cell”) are not considered as “reactions” here and are not part of the stoichiometry matrix. Now let vr denote the flux through reaction r, in the direction specified by the metabolic network. ​ ​ ​ ​ ​ For example, v1 = 100 in the urea cycle network would denote a flux of 100 (some kind of molar ​ ​ ​ amounts per unit time) through the Cps1 reaction in the forward direction, while negative values mean flux in the reverse direction.
    [Show full text]