Genome-Scale Architecture of Small Molecule Regulatory Networks and the Fundamental Trade-Off between Regulation and Enzymatic Activity

The Calvin Benson cycle in bacteria: New insights from systems biology

The Calvin Benson cycle in phototrophic and chemolithoautotrophic bacteria has ecological and biotechnological importance, which has motivated study of its regulation. I review recent advances in our understanding of how the Calvin Benson cycle is regulated in bacteria and the technologies used to elucidate regulation and modify it, and highlight differences between and photoautotrophic and chemolithoautotrophic models. Systems biology studies have shown that in oxygenic phototrophic bacteria, Calvin Benson cycle enzymes are extensively regulated at post-transcriptional and post-translational levels, with multiple enzyme activities connected to cellular redox status through thioredoxin. In chemolithoautotrophic bacteria, regulation is primarily at the transcriptional level, with effector metabolites transducing cell status, though new methods should now allow facile, proteome-wide exploration of biochemical regulation in these models. A biotechnological objective is to enhance CO2 fixation in the cycle and partition that carbon to a product of interest. Flux control of CO2 fixation is distributed over multiple enzymes, and attempts to modulate gene Calvin cycle gene expression show a robust homeostatic regulation of growth rate, though the synthesis rates of products can be significantly increased. Therefore, de-regulation of cycle enzymes through protein engineering may be necessary to increase fluxes. Non-canonical Calvin Benson cycles, if implemented with synthetic biology, could have reduced energy demand and enzyme loading, thus increasing the attractiveness of these bacteria for industrial applications.

Prediction and integration of metabolite-protein interactions with genome-scale metabolic models

Metabolites, as small molecules, can act not only as substrates to enzymes, but also as effectors of activity of proteins with different functions, thereby affecting various cellular processes. While several experimental techniques have started to catalogue the metabolite-protein interactions (MPIs) present in different cellular contexts, characterizing the functional relevance of MPIs remains a challenging problem. Computational approaches from the constrained-based modeling framework allow for predicting MPIs and integrating their effects in the in silico analysis of metabolic and physiological phenotypes, like cell growth. Here, we provide a classification of all existing constraint-based approaches that predict and integrate MPIs using genome-scale metabolic networks as input. In addition, we benchmark the performance of the approaches to predict MPIs in a comparative study using different features extracted from the model structure and predicted metabolic phenotypes with the state-of-the-art metabolic networks of Escherichia coli and Saccharomyces cerevisiae. Lastly, we provide an outlook for future, feasible directions to expand the consideration of MPIs in constraint-based modeling approaches with wide biotechnological applications.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell's entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

In silico prediction of mutation sites for anthranilate synthase from Serratia marcesens to deregulate tryptophan feedback inhibition

Allosteric feedback inhibition of the committed step in amino acid biosynthetic pathways is a major concern for production of amino acids at industrial scale. Anthranilate synthase (AS) catalyzes the first reaction of tryptophan biosynthetic pathway found in microorganisms and is feedback inhibited by its own product i.e. tryptophan. Here, we identified new mutant sites in AS using computational mutagenesis approach. MD simulations (20 ns) followed by MMPBSA and per residue decomposition energy analysis identified seven amino acid residues with best binding affinity for tryptophan. All 19 mutant structures were generated for each identified amino acid residue followed by simulation to evaluate effect of mutation on protein stability. Later, molecular docking studies were employed to generate mutant-tryptophan complex and structures with binding energies (kcal/mol) much higher than wild-type AS were selected. Finally, two mutants i.e., S37W and S37H were identified on the basis of positive binding scores and loss of tryptophan binding inside pocket. Further, MD simulations run for 200 ns were performed over these mutant-tryptophan complexes followed by RMSD, RMSF, radius of gyration , solvent accessible surface area , intra-protein hydrogen bond numbers, principal component analysis, free energy landscape (FEL) and secondary structure analysis to rationale effect of mutations on stability of protein. Cross correlation analysis of mutant site amino acids (S37W) with key residues of catalytic site (G325, T326, H395 and G482) was done to evaluate the effect of mutations on catalytic site conformation. Current computational mutagenesis approach predicted two mutants S37W and S37H with proposed deregulated feedback inhibition by tryptophan and retained catalytic activity.Communicated by Ramaswamy H. Sarma.

Transcriptome-guided metabolic network analysis reveals rearrangements of carbon flux distribution in Neisseria gonorrhoeae during neutrophil co-culture

The ability of bacterial pathogens to metabolically adapt to the environmental conditions of their hosts is critical to both colonization and invasive disease. Infection with Neisseria gonorrhoeae (the gonococcus, Gc) is characterized by the influx of neutrophils [polymorphonuclear leukocytes (PMNs)], which fail to clear the bacteria and make antimicrobial products that can exacerbate tissue damage. The inability of the human host to clear Gc infection is particularly concerning in light of the emergence of strains that are resistant to all clinically recommended antibiotics. Bacterial metabolism represents a promising target for the development of new therapeutics against Gc. Here, we generated a curated genome-scale metabolic network reconstruction (GENRE) of Gc strain FA1090. This GENRE links genetic information to metabolic phenotypes and predicts Gc biomass synthesis and energy consumption. We validated this model with published data and in new results reported here. Contextualization of this model using the transcriptional profile of Gc exposed to PMNs revealed substantial rearrangements of Gc central metabolism and induction of Gc nutrient acquisition strategies for alternate carbon source use. These features enhanced the growth of Gc in the presence of neutrophils. From these results, we conclude that the metabolic interplay between Gc and PMNs helps define infection outcomes. The use of transcriptional profiling and metabolic modeling to reveal new mechanisms by which Gc persists in the presence of PMNs uncovers unique aspects of metabolism in this fastidious bacterium, which could be targeted to block infection and thereby reduce the burden of gonorrhea in the human population. IMPORTANCE The World Health Organization designated Gc as a high-priority pathogen for research and development of new antimicrobials. Bacterial metabolism is a promising target for new antimicrobials, as metabolic enzymes are widely conserved among bacterial strains and are critical for nutrient acquisition and survival within the human host. Here we used genome-scale metabolic modeling to characterize the core metabolic pathways of this fastidious bacterium and to uncover the pathways used by Gc during culture with primary human immune cells. These analyses revealed that Gc relies on different metabolic pathways during co-culture with human neutrophils than in rich media. Conditionally essential genes emerging from these analyses were validated experimentally. These results show that metabolic adaptation in the context of innate immunity is important to Gc pathogenesis. Identifying the metabolic pathways used by Gc during infection can highlight new therapeutic targets for drug-resistant gonorrhea.

Principles of metabolome conservation in animals

Metabolite levels shape cellular physiology and disease susceptibility, yet the general principles governing metabolome evolution are largely unknown. Here, we introduce a measure of conservation of individual metabolite levels among related species. By analyzing multispecies tissue metabolome datasets in phylogenetically diverse mammals and fruit flies, we show that conservation varies extensively across metabolites. Three major functional properties, metabolite abundance, essentiality, and association with human diseases predict conservation, highlighting a striking parallel between the evolutionary forces driving metabolome and protein sequence conservation. Metabolic network simulations recapitulated these general patterns and revealed that abundant metabolites are highly conserved due to their strong coupling to key metabolic fluxes in the network. Finally, we show that biomarkers of metabolic diseases can be distinguished from other metabolites simply based on evolutionary conservation, without requiring any prior clinical knowledge. Overall, this study uncovers simple rules that govern metabolic evolution in animals and implies that most tissue metabolome differences between species are permitted, rather than favored by natural selection. More broadly, our work paves the way toward using evolutionary information to identify biomarkers, as well as to detect pathogenic metabolome alterations in individual patients.

Insight into de-regulation of amino acid feedback inhibition: a focus on structure analysis method

Regulation of amino acid's biosynthetic pathway is of significant importance to maintain homeostasis and cell functions. Amino acids regulate their biosynthetic pathway by end-product feedback inhibition of enzymes catalyzing committed steps of a pathway. Discovery of new feedback resistant enzyme variants to enhance industrial production of amino acids is a key objective in industrial biotechnology. Deregulation of feedback inhibition has been achieved for various enzymes using in vitro and in silico mutagenesis techniques. As enzyme's function, its substrate binding capacity, catalysis activity, regulation and stability are dependent on its structural characteristics, here, we provide detailed structural analysis of all feedback sensitive enzyme targets in amino acid biosynthetic pathways. Current review summarizes information regarding structural characteristics of various enzyme targets and effect of mutations on their structures and functions especially in terms of deregulation of feedback inhibition. Furthermore, applicability of various experimental as well as computational mutagenesis techniques to accomplish feedback resistance has also been discussed in detail to have an insight into various aspects of research work reported in this particular field of study.

Prediction of metabolite-protein interactions based on integration of machine learning and constraint-based modeling

Motivation

Metabolite-protein interactions play an important role in regulating protein functions and metabolism. Yet, predictions of metabolite-protein interactions using genome-scale metabolic networks are lacking. Here, we fill this gap by presenting a computational framework, termed SARTRE, that employs features corresponding to shadow prices determined in the context of flux variability analysis to predict metabolite-protein interactions using supervised machine learning.

Results

By using gold standards for metabolite-protein interactomes and well-curated genome-scale metabolic models of Escherichia coli and Saccharomyces cerevisiae, we found that the implementation of SARTRE with random forest classifiers accurately predicts metabolite-protein interactions, supported by an average area under the receiver operating curve of 0.86 and 0.85, respectively. Ranking of features based on their importance for classification demonstrated the key role of shadow prices in predicting metabolite-protein interactions. The quality of predictions is further supported by the excellent agreement of the organism-specific classifiers on unseen interactions shared between the two model organisms. Further, predictions from SARTRE are highly competitive against those obtained from a recent deep-learning approach relying on a variety of protein and metabolite features. Together, these findings show that features extracted from constraint-based analyses of metabolic networks pave the way for understanding the functional roles of the interactions between proteins and small molecules.

Availability and implementation

https://github.com/fayazsoleymani/SARTRE.

A Reduction of Transcriptional Regulation in Aquatic Oligotrophic Microorganisms Enhances Fitness in Nutrient-Poor Environments

In this review, we consider the regulatory strategies of aquatic oligotrophs, microbial cells that are adapted to thrive under low-nutrient concentrations in oceans, lakes, and other aquatic ecosystems. Many reports have concluded that oligotrophs use less transcriptional regulation than copiotrophic cells, which are adapted to high nutrient concentrations and are far more common subjects for laboratory investigations of regulation. It is theorized that oligotrophs have retained alternate mechanisms of regulation, such as riboswitches, that provide shorter response times and smaller amplitude responses and require fewer cellular resources. We examine the accumulated evidence for distinctive regulatory strategies in oligotrophs. We explore differences in the selective pressures copiotrophs and oligotrophs encounter and ask why, although evolutionary history gives copiotrophs and oligotrophs access to the same regulatory mechanisms, they might exhibit distinctly different patterns in how these mechanisms are used. We discuss the implications of these findings for understanding broad patterns in the evolution of microbial regulatory networks and their relationships to environmental niche and life history strategy. We ask whether these observations, which have emerged from a decade of increased investigation of the cell biology of oligotrophs, might be relevant to recent discoveries of many microbial cell lineages in nature that share with oligotrophs the property of reduced genome size.

Toward mechanistic modeling and rational engineering of plant respiration

Plant respiration not only provides energy to support all cellular processes, including biomass production, but also plays a major role in the global carbon cycle. Therefore, modulation of plant respiration can be used to both increase the plant yield and mitigate the effects of global climate change. Mechanistic modeling of plant respiration at sufficient biochemical detail can provide key insights for rational engineering of this process. Yet, despite its importance, plant respiration has attracted considerably less modeling effort in comparison to photosynthesis. In this update review, we highlight the advances made in modeling of plant respiration, emphasizing the gradual but important change from phenomenological to models based on first principles. We also provide a detailed account of the existing resources that can contribute to resolving the challenges in modeling plant respiration. These resources point at tangible improvements in the representation of cellular processes that contribute to CO2 evolution and consideration of kinetic properties of underlying enzymes to facilitate mechanistic modeling. The update review emphasizes the need to couple biochemical models of respiration with models of acclimation and adaptation of respiration for their effective usage in guiding breeding efforts and improving terrestrial biosphere models tailored to future climate scenarios.

The Knowns and Unknowns in Protein-Metabolite Interactions

Increasing attention has been focused on the study of protein-metabolite interactions (PMI), which play a key role in regulating protein functions and directing an orchestra of cellular processes. The investigation of PMIs is complicated by the fact that many such interactions are extremely short-lived, which requires very high resolution in order to detect them. As in the case of protein-protein interactions, protein-metabolite interactions are still not clearly defined. Existing assays for detecting protein-metabolite interactions have an additional limitation in the form of a limited capacity to identify interacting metabolites. Thus, although recent advances in mass spectrometry allow the routine identification and quantification of thousands of proteins and metabolites today, they still need to be improved to provide a complete inventory of biological molecules, as well as all interactions between them. Multiomic studies aimed at deciphering the implementation of genetic information often end with the analysis of changes in metabolic pathways, as they constitute one of the most informative phenotypic layers. In this approach, the quantity and quality of knowledge about PMIs become vital to establishing the full scope of crosstalk between the proteome and the metabolome in a biological object of interest. In this review, we analyze the current state of investigation into the detection and annotation of protein-metabolite interactions, describe the recent progress in developing associated research methods, and attempt to deconstruct the very term "interaction" to advance the field of interactomics further.

A Systems Biology Approach To Disentangle the Direct and Indirect Effects of Global Transcription Factors on Gene Expression in Escherichia coli

Delineating the pleiotropic effects of global transcriptional factors (TFs) is critical for understanding the system-wide regulatory response in a particular environment. Currently, with the availability of genome-wide TF binding and gene expression data for Escherichia coli, several gene targets can be assigned to the global TFs, albeit inconsistently. Here, using a systematic integrated approach with emphasis on metabolism, we characterized and quantified the direct effects as well as the growth rate-mediated indirect effects of global TFs using deletion mutants of FNR, ArcA, and IHF regulators (focal TFs) under glucose fermentative conditions. This categorization enabled us to disentangle the dense connections seen within the transcriptional regulatory network (TRN) and determine the exact nature of focal TF-driven epistatic interactions with other global and pathway-specific local regulators (iTFs). We extended our analysis to combinatorial deletions of these focal TFs to determine their cross talk effects as well as conserved patterns of regulatory interactions. Moreover, we predicted with high confidence several novel metabolite-iTF interactions using inferred iTF activity changes arising from the allosteric effects of the intracellular metabolites perturbed as a result of the absence of focal TFs. Further, using compendium level computational analyses, we revealed not only the coexpressed genes regulated by these focal TFs but also the coordination of the direct and indirect target expression in the context of the economy of intracellular metabolites. Overall, this study leverages the fundamentals of TF-driven regulation, which could serve as a better template for deciphering mechanisms underlying complex phenotypes. IMPORTANCE Understanding the pleiotropic effects of global TFs on gene expression and their relevance underlying a specific response in a particular environment has been challenging. Here, we distinguish the TF-driven direct effects and growth rate-mediated indirect effects on gene expression using single- and double-deletion mutants of FNR, ArcA, and IHF regulators under anaerobic glucose fermentation. Such dissection assists us in unraveling the precise nature of interactions existing between the focal TF(s) and several other TFs, including those altered by allosteric effects of intracellular metabolites. We were able to recapitulate the previously known metabolite-TF interactions and predict novel interactions with high confidence. Furthermore, we determined that the direct and indirect gene expression have a strong connection with each other when analyzed using the coexpressed- or coregulated-gene approach. Deciphering such regulatory patterns explicitly from the metabolism point of view would be valuable in understanding other unpredicted complex regulation existing in nature.

Protein-Metabolite Interactions Shape Cellular Metabolism and Physiology

Protein-metabolite interactions regulate many important cellular processes but still remain understudied. Recent technological advancements are gradually uncovering the complexity of the protein-metabolite interactome. Here, we highlight some classic and recent examples of how protein metabolite interactions regulate metabolism, both locally and globally, and how this contributes to cellular physiology. We also discuss the importance of these interactions in diseases such as cancer.

Characterization of an Entner-Doudoroff pathway-activated Escherichia coli

BACKGROUND

Escherichia coli have both the Embden-Meyerhof-Parnas pathway (EMPP) and Entner-Doudoroff pathway (EDP) for glucose breakdown, while the EDP primarily remains inactive for glucose metabolism. However, EDP is a more favorable route than EMPP for the production of certain products.

RESULTS

EDP was activated by deleting the pfkAB genes in conjunction with subsequent adaptive laboratory evolution (ALE). The evolved strains acquired mutations in transcriptional regulatory genes for glycolytic process (crp, galR, and gntR) and in glycolysis-related genes (gnd, ptsG, and talB). The genotypic, transcriptomic and phenotypic analyses of those mutations deepen our understanding of their beneficial effects on cellulosic biomass bio-conversion. On top of these scientific understandings, we further engineered the strain to produce higher level of lycopene and 3-hydroxypropionic acid.

CONCLUSIONS

These results indicate that the E. coli strain has innate capability to use EDP in lieu of EMPP for glucose metabolism, and this versatility can be harnessed to further engineer E. coli for specific biotechnological applications.

Quantitative modeling of pentose phosphate pathway response to oxidative stress reveals a cooperative regulatory strategy

Living cells use signaling and regulatory mechanisms to adapt to environmental stresses. Adaptation to oxidative stress involves the regulation of many enzymes in both glycolysis and pentose phosphate pathways (PPP), so as to support PPP-driven NADPH recycling for antioxidant defense. The underlying regulatory logic is investigated by developing a kinetic modeling approach fueled with metabolomics and 13C-fluxomics datasets from human fibroblast cells. Bayesian parameter estimation and phenotypic analysis of models highlight complementary roles for several metabolite-enzyme regulations. Specifically, carbon flux rerouting into PPP involves a tight coordination between the upregulation of G6PD activity concomitant to a decreased NADPH/NADP+ ratio and the differential control of downward and upward glycolytic fluxes through the joint inhibition of PGI and GAPD enzymes. Such functional interplay between distinct regulatory feedbacks promotes efficient detoxification and homeostasis response over a broad range of stress level, but can also explain paradoxical pertubation phenotypes for instance reported for 6PGD modulation in mammalian cells.

Multi-omics-based label-free metabolic flux inference reveals obesity-associated dysregulatory mechanisms in liver glucose metabolism

Glucose homeostasis is maintained by modulation of metabolic flux. Enzymes and metabolites regulate the involved metabolic pathways. Dysregulation of glucose homeostasis is a pathological event in obesity. Analyzing metabolic pathways and the mechanisms contributing to obesity-associated dysregulation in vivo is challenging. Here, we introduce OMELET: Omics-Based Metabolic Flux Estimation without Labeling for Extended Trans-omic Analysis. OMELET uses metabolomic, proteomic, and transcriptomic data to identify relative changes in metabolic flux, and to calculate contributions of metabolites, enzymes, and transcripts to the changes in metabolic flux. By evaluating the livers of fasting ob/ob mice, we found that increased metabolic flux through gluconeogenesis resulted primarily from increased transcripts, whereas that through the pyruvate cycle resulted from both increased transcripts and changes in substrates of metabolic enzymes. With OMELET, we identified mechanisms underlying the obesity-associated dysregulation of metabolic flux in the liver.

•

We developed OMELET to infer metabolic flux from label-free multi-omic data

•

Contributions of metabolites, enzymes, and transcripts for flux were inferred

•

Gluconeogenic flux increased in fasting ob/ob mice by increased transcripts

•

Increased pyruvate cycle fluxes were led by increased transcripts and substrates

Global pleiotropic effects in adaptively evolved Escherichia coli lacking CRP reveal molecular mechanisms that define the growth physiology

Evolution facilitates emergence of fitter phenotypes by efficient allocation of cellular resources in conjunction with beneficial mutations. However, system-wide pleiotropic effects that redress the perturbations to the apex node of the transcriptional regulatory networks remain unclear. Here, we elucidate that absence of global transcriptional regulator CRP in Escherichia coli results in alterations in key metabolic pathways under glucose respiratory conditions, favouring stress- or hedging-related functions over growth-enhancing functions. Further, we disentangle the growth-mediated effects from the CRP regulation-specific effects on these metabolic pathways. We quantitatively illustrate that the loss of CRP perturbs proteome efficiency, as evident from metabolic as well as ribosomal proteome fractions, that corroborated with intracellular metabolite profiles. To address how E. coli copes with such systemic defect, we evolved Δcrp mutant in the presence of glucose. Besides acquiring mutations in the promoter of glucose transporter ptsG, the evolved populations recovered the metabolic pathways to their pre-perturbed state coupled with metabolite re-adjustments, which altogether enabled increased growth. By contrast to Δcrp mutant, the evolved strains remodelled their proteome efficiency towards biomass synthesis, albeit at the expense of carbon efficiency. Overall, we comprehensively illustrate the genetic and metabolic basis of pleiotropic effects, fundamental for understanding the growth physiology.

Protein allocation and utilization in the versatile chemolithoautotroph Cupriavidus necator

Bacteria must balance the different needs for substrate assimilation, growth functions, and resilience in order to thrive in their environment. Of all cellular macromolecules, the bacterial proteome is by far the most important resource and its size is limited. Here, we investigated how the highly versatile 'knallgas' bacterium Cupriavidus necator reallocates protein resources when grown on different limiting substrates and with different growth rates. We determined protein quantity by mass spectrometry and estimated enzyme utilization by resource balance analysis modeling. We found that C. necator invests a large fraction of its proteome in functions that are hardly utilized. Of the enzymes that are utilized, many are present in excess abundance. One prominent example is the strong expression of CBB cycle genes such as Rubisco during growth on fructose. Modeling and mutant competition experiments suggest that CO₂-reassimilation through Rubisco does not provide a fitness benefit for heterotrophic growth, but is rather an investment in readiness for autotrophy.

Conservation of metabolic regulation by phosphorylation and non-covalent small-molecule interactions

Here, we review extant observations of protein phosphorylation and small-molecule interactions in metabolism and ask which of their specific regulatory functions are conserved in Escherichia coli and Homo sapiens. While the number of phosphosites is dramatically higher in humans, the number of metabolite-protein interactions remains largely constant. Moreover, we found the regulatory logic of metabolite-protein interactions, and in many cases also the effector molecules, to be conserved. Post-translational regulation through phosphorylation does not appear to replace this regulation in human but rather seems to add additional opportunities for fine-tuning and more complex responses. The abundance of metabolite-protein interactions in metabolism, their conserved cross-species abundance, and the apparent conservation of regulatory logic across enormous phylogenetic distance demonstrate their relevance for maintaining cellular homeostasis in these ancient biological processes.

Combination of network and molecule structure accurately predicts competitive inhibitory interactions

Mining of metabolite-protein interaction networks facilitates the identification of design principles underlying the regulation of different cellular processes. However, identification and characterization of the regulatory role that metabolites play in interactions with proteins on a genome-scale level remains a pressing task. Based on availability of high-quality metabolite-protein interaction networks and genome-scale metabolic networks, here we propose a supervised machine learning approach, called CIRI that determines whether or not a metabolite is involved in a competitive inhibitory regulatory interaction with an enzyme. First, we show that CIRI outperforms the naive approach based on a structural similarity threshold for a putative competitive inhibitor and the substrates of a metabolic reaction. We also validate the performance of CIRI on several unseen data sets and databases of metabolite-protein interactions not used in the training, and demonstrate that the classifier can be effectively used to predict competitive inhibitory interactions. Finally, we show that CIRI can be employed to refine predictions about metabolite-protein interactions from a recently proposed PROMIS approach that employs metabolomics and proteomics profiles from size exclusion chromatography in E. coli to predict metabolite-protein interactions. Altogether, CIRI fills a gap in cataloguing metabolite-protein interactions and can be used in directing future machine learning efforts to categorize the regulatory type of these interactions.

Global Transcriptional Regulators Fine-Tune the Translational and Metabolic Efficiency for Optimal Growth of Escherichia coli

Global transcriptional regulators coordinate complex genetic interactions that bestow better adaptability for an organism against external and internal perturbations. These transcriptional regulators are known to control an enormous array of genes with diverse functionalities. However, regulator-driven molecular mechanisms that underpin precisely tuned translational and metabolic processes conducive for rapid exponential growth remain obscure. Here, we comprehensively reveal the fundamental role of global transcriptional regulators FNR, ArcA, and IHF in sustaining translational and metabolic efficiency under glucose fermentative conditions in Escherichia coli By integrating high-throughput gene expression profiles and absolute intracellular metabolite concentrations, we illustrate that these regulators are crucial in maintaining nitrogen homeostasis, govern expression of otherwise unnecessary or hedging genes, and exert tight control on metabolic bottleneck steps. Furthermore, we characterize changes in expression and activity profiles of other coregulators associated with these dysregulated metabolic pathways, determining the regulatory interactions within the transcriptional regulatory network. Such systematic findings emphasize their importance in optimizing the proteome allocation toward metabolic enzymes as well as ribosomes, facilitating condition-specific phenotypic outcomes. Consequentially, we reveal that disruption of this inherent trade-off between ribosome and metabolic proteome economy due to the loss of regulators resulted in lowered growth rates. Moreover, our findings reinforce that the accumulations of intracellular metabolites in the event of proteome repartitions negatively affects the glucose uptake. Overall, by extending the three-partition proteome allocation theory concordant with multi-omics measurements, we elucidate the physiological consequences of loss of global regulators on central carbon metabolism restraining the organism to attain maximal biomass synthesis.IMPORTANCE Cellular proteome allocation in response to environmental or internal perturbations governs their eventual phenotypic outcome. This entails striking an effective balance between amino acid biosynthesis by metabolic proteins and its consumption by ribosomes. However, the global transcriptional regulator-driven molecular mechanisms that underpin their coordination remains unexplored. Here, we emphasize that global transcriptional regulators, known to control expression of a myriad of genes, are fundamental for precisely tuning the translational and metabolic efficiencies that define the growth optimality. Towards this, we systematically characterized the single deletion effect of FNR, ArcA, and IHF regulators of Escherichia coli on exponential growth under anaerobic glucose fermentative conditions. Their absence disrupts the stringency of proteome allocation, which manifests as impairment in key metabolic processes and an accumulation of intracellular metabolites. Furthermore, by incorporating an extension to the empirical growth laws, we quantitatively demonstrate the general design principles underlying the existence of these regulators in E. coli.

Integration of relative metabolomics and transcriptomics time-course data in a metabolic model pinpoints effects of ribosome biogenesis defects on Arabidopsis thaliana metabolism

Ribosome biogenesis is tightly associated to plant metabolism due to the usage of ribosomes in the synthesis of proteins necessary to drive metabolic pathways. Given the central role of ribosome biogenesis in cell physiology, it is important to characterize the impact of different components involved in this process on plant metabolism. Double mutants of the Arabidopsis thaliana cytosolic 60S maturation factors REIL1 and REIL2 do not resume growth after shift to moderate 10 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }\hbox {C}$$\end{document}∘C chilling conditions. To gain mechanistic insights into the metabolic effects of this ribosome biogenesis defect on metabolism, we developed TC-iReMet2, a constraint-based modelling approach that integrates relative metabolomics and transcriptomics time-course data to predict differential fluxes on a genome-scale level. We employed TC-iReMet2 with metabolomics and transcriptomics data from the Arabidopsis Columbia 0 wild type and the reil1-1 reil2-1 double mutant before and after cold shift. We identified reactions and pathways that are highly altered in a mutant relative to the wild type. These pathways include the Calvin–Benson cycle, photorespiration, gluconeogenesis, and glycolysis. Our findings also indicated differential NAD(P)/NAD(P)H ratios after cold shift. TC-iReMet2 allows for mechanistic hypothesis generation and interpretation of system biology experiments related to metabolic fluxes on a genome-scale level.

Enzyme activities predicted by metabolite concentrations and solvent capacity in the cell

Experimental measurements or computational model predictions of the post-translational regulation of enzymes needed in a metabolic pathway is a difficult problem. Consequently, regulation is mostly known only for well-studied reactions of central metabolism in various model organisms. In this study, we use two approaches to predict enzyme regulation policies and investigate the hypothesis that regulation is driven by the need to maintain the solvent capacity in the cell. The first predictive method uses a statistical thermodynamics and metabolic control theory framework while the second method is performed using a hybrid optimization-reinforcement learning approach. Efficient regulation schemes were learned from experimental data that either agree with theoretical calculations or result in a higher cell fitness using maximum useful work as a metric. As previously hypothesized, regulation is herein shown to control the concentrations of both immediate and downstream product concentrations at physiological levels. Model predictions provide the following two novel general principles: (1) the regulation itself causes the reactions to be much further from equilibrium instead of the common assumption that highly non-equilibrium reactions are the targets for regulation; and (2) the minimal regulation needed to maintain metabolite levels at physiological concentrations maximizes the free energy dissipation rate instead of preserving a specific energy charge. The resulting energy dissipation rate is an emergent property of regulation which may be represented by a high value of the adenylate energy charge. In addition, the predictions demonstrate that the amount of regulation needed can be minimized if it is applied at the beginning or branch point of a pathway, in agreement with common notions. The approach is demonstrated for three pathways in the central metabolism of E. coli (gluconeogenesis, glycolysis-tricarboxylic acid (TCA) and pentose phosphate-TCA) that each require different regulation schemes. It is shown quantitatively that hexokinase, glucose 6-phosphate dehydrogenase and glyceraldehyde phosphate dehydrogenase, all branch points of pathways, play the largest roles in regulating central metabolism.

Regulatory interventions improve the biosynthesis of limiting amino acids from methanol carbon to improve synthetic methylotrophy in Escherichia coli

Synthetic methylotrophy aims to engineer methane and methanol utilization pathways in platform hosts like Escherichia coli for industrial bioprocessing of natural gas and biogas. While recent attempts to engineer synthetic methylotrophs have proved successful, autonomous methylotrophy, that is, the ability to utilize methane or methanol as sole carbon and energy substrates, has not yet been realized. Here, we address an important limitation of autonomous methylotrophy in E. coli: the inability of the organism to synthesize several amino acids when grown on methanol. We targeted global and local amino acid regulatory networks. Those include removal of amino acid allosteric feedback inhibition (argA^H15Y , ilvA^L447F , hisG^E271K , leuA^G462D , proB^D107N , thrA^S345F , trpE^S40F ), knockouts of transcriptional repressors (ihfA, metJ); and overexpression of amino acid biosynthetic operons (hisGDCBHAFI, leuABCD, thrABC, trpEDCBA) and transcriptional regulators (crp, purR). Compared to the parent methylotrophic E. coli strain that was unable to synthesize these amino acids from methanol carbon, these strategies resulted in improved biosynthesis of limiting proteinogenic amino acids (histidine, leucine, lysine, methionine, phenylalanine, threonine, tyrosine) from methanol carbon. In several cases, improved amino acid biosynthesis from methanol carbon led to improvements in methylotrophic growth in methanol minimal medium supplemented with a small amount of yeast extract. This study addresses a key limitation currently preventing autonomous methylotrophy in E. coli and possibly other synthetic methylotrophs and provides insight as to how this limitation can be alleviated via global and local regulatory modifications.

Compartment and hub definitions tune metabolic networks for metabolomic interpretations

BACKGROUND

Metabolic networks represent all chemical reactions that occur between molecular metabolites in an organism's cells. They offer biological context in which to integrate, analyze, and interpret omic measurements, but their large scale and extensive connectivity present unique challenges. While it is practical to simplify these networks by placing constraints on compartments and hubs, it is unclear how these simplifications alter the structure of metabolic networks and the interpretation of metabolomic experiments.

RESULTS

We curated and adapted the latest systemic model of human metabolism and developed customizable tools to define metabolic networks with and without compartmentalization in subcellular organelles and with or without inclusion of prolific metabolite hubs. Compartmentalization made networks larger, less dense, and more modular, whereas hubs made networks larger, more dense, and less modular. When present, these hubs also dominated shortest paths in the network, yet their exclusion exposed the subtler prominence of other metabolites that are typically more relevant to metabolomic experiments. We applied the non-compartmental network without metabolite hubs in a retrospective, exploratory analysis of metabolomic measurements from 5 studies on human tissues. Network clusters identified individual reactions that might experience differential regulation between experimental conditions, several of which were not apparent in the original publications.

CONCLUSIONS

Exclusion of specific metabolite hubs exposes modularity in both compartmental and non-compartmental metabolic networks, improving detection of relevant clusters in omic measurements. Better computational detection of metabolic network clusters in large data sets has potential to identify differential regulation of individual genes, transcripts, and proteins.

Coupled metabolic-hydrodynamic modeling enabling rational scale-up of industrial bioprocesses

Metabolomics aims to address what and how regulatory mechanisms are coordinated to achieve flux optimality, different metabolic objectives as well as appropriate adaptations to dynamic nutrient availability. Recent decades have witnessed that the integration of metabolomics and fluxomics within the goal of synthetic biology has arrived at generating the desired bioproducts with improved bioconversion efficiency. Absolute metabolite quantification by isotope dilution mass spectrometry represents a functional readout of cellular biochemistry and contributes to the establishment of metabolic (structured) models required in systems metabolic engineering. In industrial practices, population heterogeneity arising from fluctuating nutrient availability frequently leads to performance losses, that is reduced commercial metrics (titer, rate, and yield). Hence, the development of more stable producers and more predictable bioprocesses can benefit from a quantitative understanding of spatial and temporal cell-to-cell heterogeneity within industrial bioprocesses. Quantitative metabolomics analysis and metabolic modeling applied in computational fluid dynamics (CFD)-assisted scale-down simulators that mimic industrial heterogeneity such as fluctuations in nutrients, dissolved gases, and other stresses can procure informative clues for coping with issues during bioprocessing scale-up. In previous studies, only limited insights into the hydrodynamic conditions inside the industrial-scale bioreactor have been obtained, which makes case-by-case scale-up far from straightforward. Tracking the flow paths of cells circulating in large-scale bioreactors is a highly valuable tool for evaluating cellular performance in production tanks. The "lifelines" or "trajectories" of cells in industrial-scale bioreactors can be captured using Euler-Lagrange CFD simulation. This novel methodology can be further coupled with metabolic (structured) models to provide not only a statistical analysis of cell lifelines triggered by the environmental fluctuations but also a global assessment of the metabolic response to heterogeneity inside an industrial bioreactor. For the future, the industrial design should be dependent on the computational framework, and this integration work will allow bioprocess scale-up to the industrial scale with an end in mind.

Measuring glycolytic flux in single yeast cells with an orthogonal synthetic biosensor

Metabolic heterogeneity between individual cells of a population harbors significant challenges for fundamental and applied research. Identifying metabolic heterogeneity and investigating its emergence require tools to zoom into metabolism of individual cells. While methods exist to measure metabolite levels in single cells, we lack capability to measure metabolic flux, i.e., the ultimate functional output of metabolic activity, on the single-cell level. Here, combining promoter engineering, computational protein design, biochemical methods, proteomics, and metabolomics, we developed a biosensor to measure glycolytic flux in single yeast cells. Therefore, drawing on the robust cell-intrinsic correlation between glycolytic flux and levels of fructose-1,6-bisphosphate (FBP), we transplanted the B. subtilis FBP-binding transcription factor CggR into yeast. With the developed biosensor, we robustly identified cell subpopulations with different FBP levels in mixed cultures, when subjected to flow cytometry and microscopy. Employing microfluidics, we were also able to assess the temporal FBP/glycolytic flux dynamics during the cell cycle. We anticipate that our biosensor will become a valuable tool to identify and study metabolic heterogeneity in cell populations.

Systematic identification of metabolites controlling gene expression in E. coli

Metabolism controls gene expression through allosteric interactions between metabolites and transcription factors. These interactions are usually measured with in vitro assays, but there are no methods to identify them at a genome-scale in vivo. Here we show that dynamic transcriptome and metabolome data identify metabolites that control transcription factors in E. coli. By switching an E. coli culture between starvation and growth, we induce strong metabolite concentration changes and gene expression changes. Using Network Component Analysis we calculate the activities of 209 transcriptional regulators and correlate them with metabolites. This approach captures, for instance, the in vivo kinetics of CRP regulation by cyclic-AMP. By testing correlations between all pairs of transcription factors and metabolites, we predict putative effectors of 71 transcription factors, and validate five interactions in vitro. These results show that combining transcriptomics and metabolomics generates hypotheses about metabolism-transcription interactions that drive transitions between physiological states.

Interactions between metabolites and transcription factors are known to control gene expression but analyzing these events at genome-scale is challenging. Here, the authors integrate dynamic metabolome and transcriptome data from E.coli to predict regulatory metabolite-transcription factor interactions.

From Escherichia coli mutant 13C labeling data to a core kinetic model: A kinetic model parameterization pipeline

Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13C-Metabolic Flux Analysis (13C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated.

Reserve Flux Capacity in the Pentose Phosphate Pathway by NADPH Binding Is Conserved across Kingdoms

All organisms evolved defense mechanisms to counteract oxidative stress and buildup of reactive oxygen species (ROS). To test whether a potentially conserved mechanism exists for the rapid response, we investigated immediate metabolic dynamics of Escherichia coli, yeast, and human dermal fibroblasts to oxidative stress that we found to be conserved between species. To elucidate the regulatory mechanisms that implement this metabolic response, we developed mechanistic kinetic models for each organism's central metabolism and systematically tested activation and inactivation of each irreversible reaction by each metabolite. This ensemble modeling predicts in vivo relevant metabolite-enzyme interactions based on their ability to quantitatively describe metabolite dynamics. All three species appear to inhibit their oxidative pentose phosphate pathway during normal growth by the redox cofactor NADPH and relieve this inhibition to increase the pathway flux for detoxification of ROS during stress, with the sole exception of yeast when exposed to high levels of stress.

•

Characterization of immediate metabolic response to oxidative stress

•

The metabolic response in glycolysis and PP pathway depends on stress severity

•

Identification of NADPH feedback inhibition on G6PDH as key regulatory interaction

•

The identified oxidative stress regulatory interaction is conserved across kingdoms

Systematic mapping of protein-metabolite interactions in central metabolism of Escherichia coli

Metabolite binding to proteins regulates nearly all cellular processes, but our knowledge of these interactions originates primarily from empirical in vitro studies. Here, we report the first systematic study of interactions between water-soluble proteins and polar metabolites in an entire biological subnetwork. To test the depth of our current knowledge, we chose to investigate the well-characterized Escherichia coli central metabolism. Using ligand-detected NMR, we assayed 29 enzymes towards binding events with 55 intracellular metabolites. Focusing on high-confidence interactions at a false-positive rate of 5%, we detected 98 interactions, among which purine nucleotides accounted for one-third, while 50% of all metabolites did not interact with any enzyme. In contrast, only five enzymes did not exhibit any metabolite binding and some interacted with up to 11 metabolites. About 40% of the interacting metabolites were predicted to be allosteric effectors based on low chemical similarity to their target's reactants. For five of the eight tested interactions, in vitro assays confirmed novel regulatory functions, including ATP and GTP inhibition of the first pentose phosphate pathway enzyme. With 76 new candidate regulatory interactions that have not been reported previously, we essentially doubled the number of known interactions, indicating that the presently available information about protein-metabolite interactions may only be the tip of the iceberg.

Metabolic profiling of cancer cells reveals genome-wide crosstalk between transcriptional regulators and metabolism

Transcriptional reprogramming of cellular metabolism is a hallmark of cancer. However, systematic approaches to study the role of transcriptional regulators (TRs) in mediating cancer metabolic rewiring are missing. Here, we chart a genome-scale map of TR-metabolite associations in human cells using a combined computational-experimental framework for large-scale metabolic profiling of adherent cell lines. By integrating intracellular metabolic profiles of 54 cancer cell lines with transcriptomic and proteomic data, we unraveled a large space of associations between TRs and metabolic pathways. We found a global regulatory signature coordinating glucose- and one-carbon metabolism, suggesting that regulation of carbon metabolism in cancer may be more diverse and flexible than previously appreciated. Here, we demonstrate how this TR-metabolite map can serve as a resource to predict TRs potentially responsible for metabolic transformation in patient-derived tumor samples, opening new opportunities in understanding disease etiology, selecting therapeutic treatments and in designing modulators of cancer-related TRs.

Aberrant gene expression in cancer coincides with drastic changes in metabolism. Here, the authors combined metabolome, transcriptome and proteome data in 54 cancer cell lines to uncover a genome-scale network of associations between transcriptional regulators and metabolites.

Allosteric Feedback Inhibition Enables Robust Amino Acid Biosynthesis in E. coli by Enforcing Enzyme Overabundance

Microbes must ensure robust amino acid metabolism in the face of external and internal perturbations. This robustness is thought to emerge from regulatory interactions in metabolic and genetic networks. Here, we explored the consequences of removing allosteric feedback inhibition in seven amino acid biosynthesis pathways in Escherichia coli (arginine, histidine, tryptophan, leucine, isoleucine, threonine, and proline). Proteome data revealed that enzyme levels decreased in five of the seven dysregulated pathways. Despite that, flux through the dysregulated pathways was not limited, indicating that enzyme levels are higher than absolutely needed in wild-type cells. We showed that such enzyme overabundance renders the arginine, histidine, and tryptophan pathways robust against perturbations of gene expression, using a metabolic model and CRISPR interference experiments. The results suggested a sensitive interaction between allosteric feedback inhibition and enzyme-level regulation that ensures robust yet efficient biosynthesis of histidine, arginine, and tryptophan in E. coli.

•

Amino acid biosynthesis enzymes do not normally operate at maximum capacity

•

Allosteric feedback inhibition ensures that enzymes are overabundant

•

Enzyme overabundance provides robustness against decreases in gene expression

Thermodynamic constraints on the regulation of metabolic fluxes

Nutrition and metabolism are fundamental to cellular function. Metabolic activity (i.e. rates of flow, most commonly referred to as flux) is constrained by thermodynamics and regulated by the activity of enzymes. The general principles that relate biological and physical variables to metabolic control are incompletely understood. Using metabolic control analysis and computer simulations in several models of simplified metabolic pathways, we derive analytical expressions that define relationships between thermodynamics, enzyme activity, and flux control. The relationships are further analyzed in a mathematical model of glycolysis as an example of a complex biochemical pathway. We show that metabolic pathways that are very far from equilibrium are controlled by the activity of upstream enzymes. However, in general, regulation of metabolic fluxes by an enzyme has a more adaptable pattern, which relies more on distribution of free energy among reaction steps in the pathway than on the thermodynamic properties of the given enzyme. These findings show how the control of metabolic pathways is shaped by thermodynamic constraints of the given pathway.

MoVE identifies metabolic valves to switch between phenotypic states

Metabolism is highly regulated, allowing for robust and complex behavior. This behavior can often be achieved by controlling a small number of important metabolic reactions, or metabolic valves. Here, we present a method to identify the location of such valves: the metabolic valve enumerator (MoVE). MoVE uses a metabolic model to identify genetic intervention strategies which decouple two desired phenotypes. We apply this method to identify valves which can decouple growth and production to systematically improve the rate and yield of biochemical production processes. We apply this algorithm to the production of diverse compounds and obtained solutions for over 70% of our targets, identifying a small number of highly represented valves to achieve near maximal growth and production. MoVE offers a systematic approach to identify metabolic valves using metabolic models, providing insight into the architecture of metabolic networks and accelerating the widespread implementation of dynamic flux redirection in diverse systems.

Integrating -omics data into genome-scale metabolic network models: principles and challenges

At genome scale, it is not yet possible to devise detailed kinetic models for metabolism because data on the in vivo biochemistry are too sparse. Predictive large-scale models for metabolism most commonly use the constraint-based framework, in which network structures constrain possible metabolic phenotypes at steady state. However, these models commonly leave many possibilities open, making them less predictive than desired. With increasingly available -omics data, it is appealing to increase the predictive power of constraint-based models (CBMs) through data integration. Many corresponding methods have been developed, but data integration is still a challenge and existing methods perform less well than expected. Here, we review main approaches for the integration of different types of -omics data into CBMs focussing on the methods' assumptions and limitations. We argue that key assumptions - often derived from single-enzyme kinetics - do not generally apply in the context of networks, thereby explaining current limitations. Emerging methods bridging CBMs and biochemical kinetics may allow for -omics data integration in a common framework to provide more accurate predictions.

DMPy: a Python package for automated mathematical model construction of large-scale metabolic systems

Background

Models of metabolism are often used in biotechnology and pharmaceutical research to identify drug targets or increase the direct production of valuable compounds. Due to the complexity of large metabolic systems, a number of conclusions have been drawn using mathematical methods with simplifying assumptions. For example, constraint-based models describe changes of internal concentrations that occur much quicker than alterations in cell physiology. Thus, metabolite concentrations and reaction fluxes are fixed to constant values. This greatly reduces the mathematical complexity, while providing a reasonably good description of the system in steady state. However, without a large number of constraints, many different flux sets can describe the optimal model and we obtain no information on how metabolite levels dynamically change. Thus, to accurately determine what is taking place within the cell, finer quality data and more detailed models need to be constructed.

Results

In this paper we present a computational framework, DMPy, that uses a network scheme as input to automatically search for kinetic rates and produce a mathematical model that describes temporal changes of metabolite fluxes. The parameter search utilises several online databases to find measured reaction parameters. From this, we take advantage of previous modelling efforts, such as Parameter Balancing, to produce an initial mathematical model of a metabolic pathway. We analyse the effect of parameter uncertainty on model dynamics and test how recent flux-based model reduction techniques alter system properties. To our knowledge this is the first time such analysis has been performed on large models of metabolism. Our results highlight that good estimates of at least 80% of the reaction rates are required to accurately model metabolic systems. Furthermore, reducing the size of the model by grouping reactions together based on fluxes alters the resulting system dynamics.

Conclusion

The presented pipeline automates the modelling process for large metabolic networks. From this, users can simulate their pathway of interest and obtain a better understanding of how altering conditions influences cellular dynamics. By testing the effects of different parameterisations we are also able to provide suggestions to help construct more accurate models of complete metabolic systems in the future.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0584-8) contains supplementary material, which is available to authorized users.

Reserve Flux Capacity in the Pentose Phosphate Pathway Enables Escherichia coli's Rapid Response to Oxidative Stress

To counteract oxidative stress and reactive oxygen species (ROS), bacteria evolved various mechanisms, primarily reducing ROS through antioxidant systems that utilize cofactor NADPH. Cells must stabilize NADPH levels by increasing flux through replenishing metabolic pathways like pentose phosphate (PP) pathway. Here, we investigate the mechanism enabling the rapid increase in NADPH supply by exposing Escherichia coli to hydrogen peroxide and quantifying the immediate metabolite dynamics. To systematically infer active regulatory interactions governing this response, we evaluated ensembles of kinetic models of glycolysis and PP pathway, each with different regulation mechanisms. Besides the known inactivation of glyceraldehyde 3-phosphate dehydrogenase by ROS, we reveal the important allosteric inhibition of the first PP pathway enzyme by NADPH. This NADPH feedback inhibition maintains a below maximum-capacity PP pathway flux under non-stress conditions. Relieving this inhibition instantly increases PP pathway flux upon oxidative stress. We demonstrate that reducing cells' capacity to rapidly reroute their flux through the PP pathway increases their oxidative stress sensitivity.

New concepts in feedback regulation of glucose metabolism

Glycolysis, the breakdown of glucose, is one of the most conserved and extensively studied biochemical pathways. Designing principles from chemistry and thermodynamics allow for energy production, biosynthesis and cellular communication. However, the kinetics or metabolic flux through the pathway also determines its function. Recently, there have been numerous developments that establish new allosteric interactions of glycolytic enzymes with small molecule metabolites and other mechanisms that may cooperate to allow for addition complex regulation of glycolysis. This review surveys these newfound sources of glycolysis regulation and discusses their possible roles in establishing kinetic design principles of glycolysis.