Systems-level analysis of mechanisms regulating yeast metabolic flux

Dissecting N-Glycosylation Dynamics in Chinese Hamster Ovary Cells Fed-batch Cultures using Time Course Omics Analyses

N-linked glycosylation affects the potency, safety, immunogenicity, and pharmacokinetic clearance of several therapeutic proteins including monoclonal antibodies. A robust control strategy is needed to dial in appropriate glycosylation profile during the course of cell culture processes accurately. However, N-glycosylation dynamics remains insufficiently understood owing to the lack of integrative analyses of factors that influence the dynamics, including sugar nucleotide donors, glycosyltransferases, and glycosidases. Here, an integrative approach involving multi-dimensional omics analyses was employed to dissect the temporal dynamics of glycoforms produced during fed-batch cultures of CHO cells. Several pathways including glycolysis, tricarboxylic citric acid cycle, and nucleotide biosynthesis exhibited temporal dynamics over the cell culture period. The steps involving galactose and sialic acid addition were determined as temporal bottlenecks. Our results show that galactose, and not manganese, is able to mitigate the temporal bottleneck, despite both being known effectors of galactosylation. Furthermore, sialylation is limited by the galactosylated precursors and autoregulation of cytidine monophosphate-sialic acid biosynthesis.

•

Major glycosylated species exhibit temporal dynamics during fed-batch processes

•

Key metabolic pathways linked to N-glycosylation exhibit significant temporal dynamics

•

Dynamics in nucleotide sugar donors (NSDs) directly influences glycoform heterogeneity

•

Glycoform heterogeneity can be mitigated by supplementing NSD biosynthetic precursors

Using the Human Genome-Scale Metabolic Model Recon 2 for Steady-State Flux Analysis of Cancer Cell Metabolism

Flux analysis is performed to infer intracellular metabolic activity using measured rates. By applying the highly curated human metabolic reconstruction Recon 2 as the reference model, the investigation of cancer cell metabolic fluxes can encompass the full metabolic potential of a human cell. But in its full form, Recon 2 is unsuitable for conventional metabolic flux analysis due to a large number of redundant elements. Here, we describe a procedure to reduce Recon 2 to an appropriate scale for cancer cell flux analysis such that calculated flux intervals are still informative, without compromising the opportunity to explore alternative pathways encoded in Recon 2 that may reveal novel metabolic features.

Genome-Scale Metabolic Modeling from Yeast to Human Cell Models of Complex Diseases: Latest Advances and Challenges

Genome-scale metabolic models (GEMs) are mathematical models that enable systematic analysis of metabolism. This modeling concept has been applied to study the metabolism of many organisms including the eukaryal model organism, the yeast Saccharomyces cerevisiae, that also serves as an important cell factory for production of fuels and chemicals. With the application of yeast GEMs, our knowledge of metabolism is increasing. Therefore, GEMs have also been used for modeling human cells to study metabolic diseases. Here we introduce the concept of GEMs and provide a protocol for reconstructing GEMs. Besides, we show the historic development of yeast GEMs and their applications. Also, we review human GEMs as well as their uses in the studies of complex diseases.

Metabolic programming a lean phenotype by deregulation of RNA polymerase III

As a master negative regulator of RNA polymerase (Pol) III, Maf1 modulates transcription in response to nutrients and stress to balance the production of highly abundant tRNAs, 5S rRNA, and other small noncoding RNAs with cell growth and maintenance. This regulation of Pol III transcription is important for energetic economy as mice lacking Maf1 are lean and resist weight gain on normal and high fat diets. The lean phenotype of Maf1 knockout (KO) mice is attributed in part to metabolic inefficiencies which increase the demand for cellular energy and elevate catabolic processes, including autophagy/lipophagy and lipolysis. A futile RNA cycle involving increased synthesis and turnover of Pol III transcripts has been proposed as an important driver of these changes. Here, using targeted metabolomics, we find changes in the liver of fed and fasted Maf1 KO mice consistent with the function of mammalian Maf1 as a chronic Pol III repressor. Differences in long-chain acylcarnitine levels suggest that energy demand is higher in the fed state of Maf1 KO mice versus the fasted state. Quantitative metabolite profiling supports increased activity in the TCA cycle, the pentose phosphate pathway, and the urea cycle and reveals changes in nucleotide levels and the creatine system. Metabolite profiling also confirms key predictions of the futile RNA cycle hypothesis by identifying changes in many metabolites involved in nucleotide synthesis and turnover. Thus, constitutively high levels of Pol III transcription in Maf1 KO mice reprogram central metabolic pathways and waste metabolic energy through a futile RNA cycle.

Advances in metabolic flux analysis toward genome-scale profiling of higher organisms

Methodological and technological advances have recently paved the way for metabolic flux profiling in higher organisms, like plants. However, in comparison with omics technologies, flux profiling has yet to provide comprehensive differential flux maps at a genome-scale and in different cell types, tissues, and organs. Here we highlight the recent advances in technologies to gather metabolic labeling patterns and flux profiling approaches. We provide an opinion of how recent local flux profiling approaches can be used in conjunction with the constraint-based modeling framework to arrive at genome-scale flux maps. In addition, we point at approaches which use metabolomics data without introduction of label to predict either non-steady state fluxes in a time-series experiment or flux changes in different experimental scenarios. The combination of these developments allows an experimentally feasible approach for flux-based large-scale systems biology studies.

Evolution of gene knockout strains of E. coli reveal regulatory architectures governed by metabolism

Biological regulatory network architectures are multi-scale in their function and can adaptively acquire new functions. Gene knockout (KO) experiments provide an established experimental approach not just for studying gene function, but also for unraveling regulatory networks in which a gene and its gene product are involved. Here we study the regulatory architecture of Escherichia coli K-12 MG1655 by applying adaptive laboratory evolution (ALE) to metabolic gene KO strains. Multi-omic analysis reveal a common overall schema describing the process of adaptation whereby perturbations in metabolite concentrations lead regulatory networks to produce suboptimal states, whose function is subsequently altered and re-optimized through acquisition of mutations during ALE. These results indicate that metabolite levels, through metabolite-transcription factor interactions, have a dominant role in determining the function of a multi-scale regulatory architecture that has been molded by evolution.

The function of metabolic genes in the context of regulatory networks is not well understood. Here, the authors investigate the adaptive responses of E. coli after knockout of metabolic genes and highlight the influence of metabolite levels in the evolution of regulatory function.

Computational prediction of changes in brain metabolic fluxes during Parkinson's disease from mRNA expression

Background

Parkinson’s disease is a widespread neurodegenerative disorder which affects brain metabolism. Although changes in gene expression during disease are often measured, it is difficult to predict metabolic fluxes from gene expression data. Here we explore the hypothesis that changes in gene expression for enzymes tend to parallel flux changes in biochemical reaction pathways in the brain metabolic network. This hypothesis is the basis of a computational method to predict metabolic flux changes from post-mortem gene expression measurements in Parkinson’s disease (PD) brain.

Results

We use a network model of central metabolism and optimize the correspondence between relative changes in fluxes and in gene expression. To this end we apply the Least-squares with Equalities and Inequalities algorithm integrated with Flux Balance Analysis (Lsei-FBA). We predict for PD (1) decreases in glycolytic rate and oxygen consumption and an increase in lactate production in brain cortex that correspond with measurements (2) relative flux decreases in ATP synthesis, in the malate-aspartate shuttle and midway in the TCA cycle that are substantially larger than relative changes in glucose uptake in the substantia nigra, dopaminergic neurons and most other brain regions (3) shifts in redox shuttles between cytosol and mitochondria (4) in contrast to Alzheimer’s disease: little activation of the gamma-aminobutyric acid shunt pathway in compensation for decreased alpha-ketoglutarate dehydrogenase activity (5) in the globus pallidus internus, metabolic fluxes are increased, reflecting increased functional activity.

Conclusion

Our method predicts metabolic changes from gene expression data that correspond in direction and order of magnitude with presently available experimental observations during Parkinson’s disease, indicating that the hypothesis may be useful for some biochemical pathways. Lsei-FBA generates predictions of flux distributions in neurons and small brain regions for which accurate metabolic flux measurements are not yet possible.

Trans-omic Analysis Reveals Selective Responses to Induced and Basal Insulin across Signaling, Transcriptional, and Metabolic Networks

The concentrations of insulin selectively regulate multiple cellular functions. To understand how insulin concentrations are interpreted by cells, we constructed a trans-omic network of insulin action in FAO hepatoma cells using transcriptomic data, western blotting analysis of signaling proteins, and metabolomic data. By integrating sensitivity into the trans-omic network, we identified the selective trans-omic networks stimulated by high and low doses of insulin, denoted as induced and basal insulin signals, respectively. The induced insulin signal was selectively transmitted through the pathway involving Erk to an increase in the expression of immediate-early and upregulated genes, whereas the basal insulin signal was selectively transmitted through a pathway involving Akt and an increase of Foxo phosphorylation and a reduction of downregulated gene expression. We validated the selective trans-omic network in vivo by analysis of the insulin-clamped rat liver. This integrated analysis enabled molecular insight into how liver cells interpret physiological insulin signals to regulate cellular functions.

Machine Learning Predicts the Yeast Metabolome from the Quantitative Proteome of Kinase Knockouts

A challenge in solving the genotype-to-phenotype relationship is to predict a cell’s metabolome, believed to correlate poorly with gene expression. Using comparative quantitative proteomics, we found that differential protein expression in 97 Saccharomyces cerevisiae kinase deletion strains is non-redundant and dominated by abundance changes in metabolic enzymes. Associating differential enzyme expression landscapes to corresponding metabolomes using network models provided reasoning for poor proteome-metabolome correlations; differential protein expression redistributes flux control between many enzymes acting in concert, a mechanism not captured by one-to-one correlation statistics. Mapping these regulatory patterns using machine learning enabled the prediction of metabolite concentrations, as well as identification of candidate genes important for the regulation of metabolism. Overall, our study reveals that a large part of metabolism regulation is explained through coordinated enzyme expression changes. Our quantitative data indicate that this mechanism explains more than half of metabolism regulation and underlies the interdependency between enzyme levels and metabolism, which renders the metabolome a predictable phenotype.

•

The proteome of kinase knockouts is dominated by enzyme abundance changes

•

The enzyme expression profiles of kinase knockouts are non-redundant

•

Metabolism is regulated by many expression changes acting in concert

•

Machine learning accurately predicts the metabolome from enzyme abundance

Four Key Steps Control Glycolytic Flux in Mammalian Cells

Altered glycolysis is a hallmark of diseases including diabetes and cancer. Despite intensive study of the contributions of individual glycolytic enzymes, systems-level analyses of flux control through glycolysis remain limited. Here, we overexpress in two mammalian cell lines the individual enzymes catalyzing each of the 12 steps linking extracellular glucose to excreted lactate, and find substantial flux control at four steps: glucose import, hexokinase, phosphofructokinase, and lactate export (and not at any steps of lower glycolysis). The four flux-controlling steps are specifically upregulated by the Ras oncogene: optogenetic Ras activation rapidly induces the transcription of isozymes catalyzing these four steps and enhances glycolysis. At least one isozyme catalyzing each of these four steps is consistently elevated in human tumors. Thus, in the studied contexts, flux control in glycolysis is concentrated in four key enzymatic steps. Upregulation of these steps in tumors likely underlies the Warburg effect.

Dynamic Metabolomics Reveals that Insulin Primes the Adipocyte for Glucose Metabolism

Insulin triggers an extensive signaling cascade to coordinate adipocyte glucose metabolism. It is considered that the major role of insulin is to provide anabolic substrates by activating GLUT4-dependent glucose uptake. However, insulin stimulates phosphorylation of many metabolic proteins. To examine the implications of this on glucose metabolism, we performed dynamic tracer metabolomics in cultured adipocytes treated with insulin. Temporal analysis of metabolite concentrations and tracer labeling revealed rapid and distinct changes in glucose metabolism, favoring specific glycolytic branch points and pyruvate anaplerosis. Integrating dynamic metabolomics and phosphoproteomics data revealed that insulin-dependent phosphorylation of anabolic enzymes occurred prior to substrate accumulation. Indeed, glycogen synthesis was activated independently of glucose supply. We refer to this phenomenon as metabolic priming, whereby insulin signaling creates a demand-driven system to "pull" glucose into specific anabolic pathways. This complements the supply-driven regulation of anabolism by substrate accumulation and highlights an additional role for insulin action in adipocyte glucose metabolism.

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.

Transcriptomic dissection of Bradyrhizobium sp. strain ORS285 in symbiosis with Aeschynomene spp. inducing different bacteroid morphotypes with contrasted symbiotic efficiency

To circumvent the paucity of nitrogen sources in the soil legume plants establish a symbiotic interaction with nitrogen-fixing soil bacteria called rhizobia. During symbiosis, the plants form root organs called nodules, where bacteria are housed intracellularly and become active nitrogen fixers known as bacteroids. Depending on their host plant, bacteroids can adopt different morphotypes, being either unmodified (U), elongated (E) or spherical (S). E- and S-type bacteroids undergo a terminal differentiation leading to irreversible morphological changes and DNA endoreduplication. Previous studies suggest that differentiated bacteroids display an increased symbiotic efficiency (E > U and S > U). In this study, we used a combination of Aeschynomene species inducing E- or S-type bacteroids in symbiosis with Bradyrhizobium sp. ORS285 to show that S-type bacteroids present a better symbiotic efficiency than E-type bacteroids. We performed a transcriptomic analysis on E- and S-type bacteroids formed by Aeschynomene afraspera and Aeschynomene indica nodules and identified the bacterial functions activated in bacteroids and specific to each bacteroid type. Extending the expression analysis in E- and S-type bacteroids in other Aeschynomene species by qRT-PCR on selected genes from the transcriptome analysis narrowed down the set of bacteroid morphotype-specific genes. Functional analysis of a selected subset of 31 bacteroid-induced or morphotype-specific genes revealed no symbiotic phenotypes in the mutants. This highlights the robustness of the symbiotic program but could also indicate that the bacterial response to the plant environment is partially anticipatory or even maladaptive. Our analysis confirms the correlation between differentiation and efficiency of the bacteroids and provides a framework for the identification of bacterial functions that affect the efficiency of bacteroids.© 2018 Society for Applied Microbiology and John Wiley & Sons Ltd.

Mass Spectrometry-Based Method to Study Inhibitor-Induced Metabolic Redirection in the Central Metabolism of Cancer Cells

Cancer cells often respond to chemotherapeutic inhibitors by redirecting carbon flow in the central metabolism. To understand the metabolic redirections of inhibitor treatment on cancer cells, this study established a 13C-metabolic flux analysis (13C-MFA)-based method to evaluate metabolic redirection in MCF-7 breast cancer cells using mass spectrometry. A metabolic stationary state necessary for accurate 13C-MFA was confirmed during an 8-24 h window using low-dose treatments of various metabolic inhibitors. Further 13C-labeling experiments using [1-13C]glucose and [U-13C]glutamine, combined with gas chromatography-mass spectrometry (GC-MS) analysis of mass isotopomer distributions (MIDs), confirmed that an isotopic stationary state of intracellular metabolites was reached 24 h after treatment with paclitaxel (Taxol), an inhibitor of mitosis used for cancer treatment. Based on these metabolic and isotopic stationary states, metabolic flux distribution in the central metabolism of paclitaxel-treated MCF-7 cells was determined by 13C-MFA. Finally, estimations of the 95% confidence intervals showed that tricarboxylic acid cycle metabolic flux increased after paclitaxel treatment. Conversely, anaerobic glycolysis metabolic flux decreased, revealing metabolic redirections by paclitaxel inhibition. The gap between total regeneration and consumption of ATP in paclitaxel-treated cells was also found to be 1.2 times greater than controls, suggesting ATP demand was increased by paclitaxel treatment, likely due to increased microtubule polymerization. These data confirm that 13C-MFA can be used to investigate inhibitor-induced metabolic redirection in cancer cells. This will contribute to future pharmaceutical developments and understanding variable patient response to treatment.

13C-metabolic flux analysis of ethanol-assimilating Saccharomyces cerevisiae for S-adenosyl-L-methionine production

BACKGROUND

Saccharomyces cerevisiae is a host for the industrial production of S-adenosyl-L-methionine (SAM), which has been widely used in pharmaceutical and nutritional supplement industries. It has been reported that the intracellular SAM content in S. cerevisiae can be improved by the addition of ethanol during cultivation. However, the metabolic state in ethanol-assimilating S. cerevisiae remains unclear. In this study, ¹³C-metabolic flux analysis (¹³C-MFA) was conducted to investigate the metabolic regulation responsible for the high SAM production from ethanol.

RESULTS

The comparison between the metabolic flux distributions of central carbon metabolism showed that the metabolic flux levels of the tricarboxylic acid cycle and glyoxylate shunt in the ethanol culture were significantly higher than that of glucose. Estimates of the ATP balance from the ¹³C-MFA data suggested that larger amounts of excess ATP was produced from ethanol via increased oxidative phosphorylation. The finding was confirmed by the intracellular ATP level under ethanol-assimilating condition being similarly higher than glucose.

CONCLUSIONS

These results suggest that the enhanced ATP regeneration due to ethanol assimilation was critical for the high SAM accumulation.

A machine learning approach to predict metabolic pathway dynamics from time-series multiomics data

New synthetic biology capabilities hold the promise of dramatically improving our ability to engineer biological systems. However, a fundamental hurdle in realizing this potential is our inability to accurately predict biological behavior after modifying the corresponding genotype. Kinetic models have traditionally been used to predict pathway dynamics in bioengineered systems, but they take significant time to develop, and rely heavily on domain expertise. Here, we show that the combination of machine learning and abundant multiomics data (proteomics and metabolomics) can be used to effectively predict pathway dynamics in an automated fashion. The new method outperforms a classical kinetic model, and produces qualitative and quantitative predictions that can be used to productively guide bioengineering efforts. This method systematically leverages arbitrary amounts of new data to improve predictions, and does not assume any particular interactions, but rather implicitly chooses the most predictive ones.

1H-NMR as implemented in several origin of life studies artificially implies the absence of metabolism-like non-enzymatic reactions by being signal-suppressed

Background. Life depends on small subsets of chemically possible reactions. A chemical process can hence be prebiotically plausible, yet be unrelated to the origins of life. An example is the synthesis of nucleotides from hydrogen cyanide, considered prebiotically plausible, but incompatible with metabolic evolution. In contrast, only few metabolism-compatible prebiotic reactions were known until recently. Here, we question whether technical limitations may have contributed to the situation. Methods: Enzymes evolved to accelerate and control biochemical reactions. This situation dictates that compared to modern metabolic pathways, precursors to enzymatic reactions have been slower and less efficient, yielding lower metabolite quantities. This situation demands for the application of highly sensitive analytical techniques for studying ‘proto-metabolism’. We noticed that a set of proto-metabolism studies derive conclusions from the absence of metabolism-like signals, yet do not report detection limits. We here benchmark the respective 1H-NMR implementation for the ability to detect Krebs cycle intermediates, considered examples of plausible metabolic precursors. Results: Compared to metabolomics ‘gold-standard’ methods, 1H-NMR as implemented is i) at least one hundred- to thousand-fold less sensitive, ii) prone to selective metabolite loss, and iii) subject to signal suppression by Fe(II) concentrations as extrapolated from Archean sediment. In sum these restrictions mount to huge sensitivity deficits, so that even highly concentrated Krebs cycle intermediates are rendered undetectable unless the method is altered to boost sensitivity. Conclusions 1H-NMR as implemented in several origin of life studies does not achieve the sensitivity to detect cellular metabolite concentrations, let alone evolutionary precursors at even lower concentration. These studies can hence not serve as proof-of-absence for metabolism-like reactions. Origin of life theories that essentially depend on this assumption, i.e. those that consider proto-metabolism to consist of non-metabolism-like reactions derived from non-metabolic precursors like hydrogen cyanide, may have been derived from a misinterpretation of negative analytical results.

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

Metabolic flux is in part regulated by endogenous small molecules that modulate the catalytic activity of an enzyme, e.g., allosteric inhibition. In contrast to transcriptional regulation of enzymes, technical limitations have hindered the production of a genome-scale atlas of small molecule-enzyme regulatory interactions. Here, we develop a framework leveraging the vast, but fragmented, biochemical literature to reconstruct and analyze the small molecule regulatory network (SMRN) of the model organism Escherichia coli, including the primary metabolite regulators and enzyme targets. Using metabolic control analysis, we prove a fundamental trade-off between regulation and enzymatic activity, and we combine it with metabolomic measurements and the SMRN to make inferences on the sensitivity of enzymes to their regulators. Generalizing the analysis to other organisms, we identify highly conserved regulatory interactions across evolutionarily divergent species, further emphasizing a critical role for small molecule interactions in the maintenance of metabolic homeostasis.

Metabolic engineering of Pichia pastoris

Besides its use for efficient production of recombinant proteins the methylotrophic yeast Pichia pastoris (syn. Komagataella spp.) has been increasingly employed as a platform to produce metabolites of varying origin. We summarize here the impressive methodological developments of the last years to model and analyze the metabolism of P. pastoris, and to engineer its genome and metabolic pathways. Efficient methods to insert, modify or delete genes via homologous recombination and CRISPR/Cas9, supported by modular cloning techniques, have been reported. An outstanding early example of metabolic engineering in P. pastoris was the humanization of protein glycosylation. More recently the cell metabolism was engineered also to enhance the productivity of heterologous proteins. The last few years have seen an increased number of metabolic pathway design and engineering in P. pastoris, mainly towards the production of complex (secondary) metabolites. In this review, we discuss the potential role of P. pastoris as a platform for metabolic engineering, its strengths, and major requirements for future developments of chassis strains based on synthetic biology principles.

RapidRIP quantifies the intracellular metabolome of 7 industrial strains of E. coli

Fast metabolite quantification methods are required for high throughput screening of microbial strains obtained by combinatorial or evolutionary engineering approaches. In this study, a rapid RIP-LC-MS/MS (RapidRIP) method for high-throughput quantitative metabolomics was developed and validated that was capable of quantifying 102 metabolites from central, amino acid, energy, nucleotide, and cofactor metabolism in less than 5 minutes. The method was shown to have comparable sensitivity and resolving capability as compared to a full length RIP-LC-MS/MS method (FullRIP). The RapidRIP method was used to quantify the metabolome of seven industrial strains of E. coli revealing significant differences in glycolytic, pentose phosphate, TCA cycle, amino acid, and energy and cofactor metabolites were found. These differences translated to statistically and biologically significant differences in thermodynamics of biochemical reactions between strains that could have implications when choosing a host for bioprocessing.

Evolutionary Convergence of Pathway-Specific Enzyme Expression Stoichiometry

Coexpression of proteins in response to pathway-inducing signals is the founding paradigm of gene regulation. However, it remains unexplored whether the relative abundance of co-regulated proteins requires precise tuning. Here, we present large-scale analyses of protein stoichiometry and corresponding regulatory strategies for 21 pathways and 67-224 operons in divergent bacteria separated by 0.6-2 billion years. Using end-enriched RNA-sequencing (Rend-seq) with single-nucleotide resolution, we found that many bacterial gene clusters encoding conserved pathways have undergone massive divergence in transcript abundance and architectures via remodeling of internal promoters and terminators. Remarkably, these evolutionary changes are compensated post-transcriptionally to maintain preferred stoichiometry of protein synthesis rates. Even more strikingly, in eukaryotic budding yeast, functionally analogous proteins that arose independently from bacterial counterparts also evolved to convergent in-pathway expression. The broad requirement for exact protein stoichiometries despite regulatory divergence provides an unexpected principle for building biological pathways both in nature and for synthetic activities.

Eight Kinetically Stable but Thermodynamically Activated Molecules that Power Cell Metabolism

Contemporary analyses of cell metabolism have called out three metabolites: ATP, NADH, and acetyl-CoA, as sentinel molecules whose accumulation represent much of the purpose of the catabolic arms of metabolism and then drive many anabolic pathways. Such analyses largely leave out how and why ATP, NADH, and acetyl-CoA (Figure 1 ) at the molecular level play such central roles. Yet, without those insights into why cells accumulate them and how the enabling properties of these key metabolites power much of cell metabolism, the underlying molecular logic remains mysterious. Four other metabolites, S-adenosylmethionine, carbamoyl phosphate, UDP-glucose, and Δ2-isopentenyl-PP play similar roles in using group transfer chemistry to drive otherwise unfavorable biosynthetic equilibria. This review provides the underlying chemical logic to remind how these seven key molecules function as mobile packets of cellular currencies for phosphoryl transfers (ATP), acyl transfers (acetyl-CoA, carbamoyl-P), methyl transfers (SAM), prenyl transfers (IPP), glucosyl transfers (UDP-glucose), and electron and ADP-ribosyl transfers (NAD(P)H/NAD(P)+) to drive metabolic transformations in and across most primary pathways. The eighth key metabolite is molecular oxygen (O2), thermodynamically activated for reduction by one electron path, leaving it kinetically stable to the vast majority of organic cellular metabolites.

Competitive resource allocation to metabolic pathways contributes to overflow metabolisms and emergent properties in cross-feeding microbial consortia

Resource scarcity is a common stress in nature and has a major impact on microbial physiology. This review highlights microbial acclimations to resource scarcity, focusing on resource investment strategies for chemoheterotrophs from the molecular level to the pathway level. Competitive resource allocation strategies often lead to a phenotype known as overflow metabolism; the resulting overflow byproducts can stabilize cooperative interactions in microbial communities and can lead to cross-feeding consortia. These consortia can exhibit emergent properties such as enhanced resource usage and biomass productivity. The literature distilled here draws parallels between in silico and laboratory studies and ties them together with ecological theories to better understand microbial stress responses and mutualistic consortia functioning.

Advances in analytical tools for high throughput strain engineering

The emergence of inexpensive, base-perfect genome editing is revolutionising biology. Modern industrial biotechnology exploits the advances in genome editing in combination with automation, analytics and data integration to build high-throughput automated strain engineering pipelines also known as biofoundries. Biofoundries replace the slow and inconsistent artisanal processes used to build microbial cell factories with an automated design-build-test cycle, considerably reducing the time needed to deliver commercially viable strains. Testing and hence learning remains relatively shallow, but recent advances in analytical chemistry promise to increase the depth of characterization possible. Analytics combined with models of cellular physiology in automated systems biology pipelines should enable deeper learning and hence a steeper pitch of the learning cycle. This review explores the progress, advances and remaining bottlenecks of analytical tools for high throughput strain engineering.

Reverse engineering the cancer metabolic network using flux analysis to understand drivers of human disease

Metabolic dysfunction has reemerged as an essential hallmark of tumorigenesis, and metabolic phenotypes are increasingly being integrated into pre-clinical models of disease. The complexity of these metabolic networks requires systems-level interrogation, and metabolic flux analysis (MFA) with stable isotope tracing present a suitable conceptual framework for such systems. Here we review efforts to elucidate mechanisms through which metabolism influences tumor growth and survival, with an emphasis on applications using stable isotope tracing and MFA. Through these approaches researchers can now quantify pathway fluxes in various in vitro and in vivo contexts to provide mechanistic insights at molecular and physiological scales respectively. Knowledge and discoveries in cancer models are paving the way toward applications in other biological contexts and disease models. In turn, MFA approaches will increasingly help to uncover new therapeutic opportunities that enhance human health.

Stoichiometric Correlation Analysis: Principles of Metabolic Functionality from Metabolomics Data

Recent advances in metabolomics technologies have resulted in high-quality (time-resolved) metabolic profiles with an increasing coverage of metabolic pathways. These data profiles represent read-outs from often non-linear dynamics of metabolic networks. Yet, metabolic profiles have largely been explored with regression-based approaches that only capture linear relationships, rendering it difficult to determine the extent to which the data reflect the underlying reaction rates and their couplings. Here we propose an approach termed Stoichiometric Correlation Analysis (SCA) based on correlation between positive linear combinations of log-transformed metabolic profiles. The log-transformation is due to the evidence that metabolic networks can be modeled by mass action law and kinetics derived from it. Unlike the existing approaches which establish a relation between pairs of metabolites, SCA facilitates the discovery of higher-order dependence between more than two metabolites. By using a paradigmatic model of the tricarboxylic acid cycle we show that the higher-order dependence reflects the coupling of concentration of reactant complexes, capturing the subtle difference between the employed enzyme kinetics. Using time-resolved metabolic profiles from Arabidopsis thaliana and Escherichia coli, we show that SCA can be used to quantify the difference in coupling of reactant complexes, and hence, reaction rates, underlying the stringent response in these model organisms. By using SCA with data from natural variation of wild and domesticated wheat and tomato accession, we demonstrate that the domestication is accompanied by loss of such couplings, in these species. Therefore, application of SCA to metabolomics data from natural variation in wild and domesticated populations provides a mechanistic way to understanding domestication and its relation to metabolic networks.

Designing and interpreting 'multi-omic' experiments that may change our understanding of biology

Most biological mechanisms involve more than one type of biomolecule, and hence operate not solely at the level of either genome, transcriptome, proteome, metabolome or ionome. Datasets resulting from single-omic analysis are rapidly increasing in throughput and quality, rendering multi-omic studies feasible. These should offer a comprehensive, structured and interactive overview of a biological mechanism. However, combining single-omic datasets in a meaningful manner has so far proved challenging, and the discovery of new biological information lags behind expectation. One reason is that experiments conducted in different laboratories can typically not to be combined without restriction. Second, the interpretation of multi-omic datasets represents a significant challenge by nature, as the biological datasets are heterogeneous not only for technical, but also for biological, chemical, and physical reasons. Here, multi-layer network theory and methods of artificial intelligence might contribute to solve these problems. For the efficient application of machine learning however, biological datasets need to become more systematic, more precise - and much larger. We conclude our review with basic guidelines for the successful set-up of a multi-omic experiment.

Preserving protein function through reversible aggregation

It is generally accepted that protein function depends on a defined 3D structure, with unfolding and aggregation dealing a final blow to functionality. A study now shows that the regulated exposure of an unstructured region in yeast pyruvate kinase triggers reversible aggregation to preserve protein function under stress.

Formulation, construction and analysis of kinetic models of metabolism: A review of modelling frameworks

Kinetic models are critical to predict the dynamic behaviour of metabolic networks. Mechanistic kinetic models for large networks remain uncommon due to the difficulty of fitting their parameters. Recent modelling frameworks promise new ways to overcome this obstacle while retaining predictive capabilities. In this review, we present an overview of the relevant mathematical frameworks for kinetic formulation, construction and analysis. Starting with kinetic formalisms, we next review statistical methods for parameter inference, as well as recent computational frameworks applied to the construction and analysis of kinetic models. Finally, we discuss opportunities and limitations hindering the development of larger kinetic reconstructions.

An engineered photoswitchable mammalian pyruvate kinase

Changes in allosteric regulation of glycolytic enzymes have been linked to metabolic reprogramming involved in cancer. Remarkably, allosteric mechanisms control enzyme function at significantly shorter time-scales compared to the long-term effects of metabolic reprogramming on cell proliferation. It remains unclear if and how the speed and reversibility afforded by rapid allosteric control of metabolic enzymes is important for cell proliferation. Tools that allow specific, dynamic modulation of enzymatic activities in mammalian cells would help address this question. Towards this goal, we have used molecular dynamics simulations to guide the design of mPKM2 internal light/oxygen/voltage-sensitive domain 2 (LOV2) fusion at position D24 (PiL[D24]), an engineered pyruvate kinase M2 (PKM2) variant that harbours an insertion of the light-sensing LOV2 domain from Avena Sativa within a region implicated in allosteric regulation by fructose 1,6-bisphosphate (FBP). The LOV2 photoreaction is preserved in the PiL[D24] chimera and causes secondary structure changes that are associated with a 30% decrease in the Km of the enzyme for phosphoenolpyruvate resulting in increased pyruvate kinase activity after light exposure. Importantly, this change in activity is reversible upon light withdrawal. Expression of PiL[D24] in cells leads to light-induced increase in labelling of pyruvate from glucose. PiL[D24] therefore could provide a means to modulate cellular glucose metabolism in a remote manner and paves the way for studying the importance of rapid allosteric phenomena in the regulation of metabolism and enzyme control.

Genome-scale modeling of yeast: chronology, applications and critical perspectives

Over the last 15 years, several genome-scale metabolic models (GSMMs) were developed for different yeast species, aiding both the elucidation of new biological processes and the shift toward a bio-based economy, through the design of in silico inspired cell factories. Here, an historical perspective of the GSMMs built over time for several yeast species is presented and the main inheritance patterns among the metabolic reconstructions are highlighted. We additionally provide a critical perspective on the overall genome-scale modeling procedure, underlining incomplete model validation and evaluation approaches and the quest for the integration of regulatory and kinetic information into yeast GSMMs. A summary of experimentally validated model-based metabolic engineering applications of yeast species is further emphasized, while the main challenges and future perspectives for the field are finally addressed.

Primordial Krebs-cycle-like non-enzymatic reactions detected by mass spectrometry and nuclear magnetic resonance

Background: Metabolism is the process of nutrient uptake and conversion, and executed by the metabolic network. Its evolutionary precursors most likely originated in non-enzymatic chemistry. To be exploitable in a Darwinian process that forms a metabolic pathway, non-enzymatic reactions need to form a chemical network that produces advantage-providing metabolites in a single, life compatible condition. In a hypothesis-generating, large-scale experiment, we recently screened iron and sulfur-rich solutions, and report that upon the formation of sulfate radicals, Krebs cycle intermediates establish metabolism-like non-enzymatic reactivity. A challenge to our results claims that the results obtained by liquid chromatography-selective reaction monitoring (LC-SRM) would not be reproducible by nuclear magnetic resonance spectroscopy (1H-NMR). Methods: This study compared the application of the two techniques to the relevant samples. Results: We detect hundred- to thousand-fold differences in the specific limits of detection between LC-SRM and 1H-NMR to detect Krebs cycle intermediates. Further, the use of 1H-NMR was found generally problematic to characterize early metabolic reactions, as Archean-sediment typical iron concentrations cause paramagnetic signal suppression. Consequently, we selected non-enzymatic Krebs cycle reactions that fall within the determined technical limits. We confirm that these proceed unequivocally as evidenced by both LC-SRM and 1H-NMR. Conclusions: These results strengthen our previous conclusions about the existence of unifying reaction conditions that enables a series of co-occurring metabolism-like non-enzymatic Krebs cycle reactions. We further discuss why constraints applying to metabolism disentangle concentration from importance of any reaction intermediates, and why evolutionary precursors to metabolic pathways must have had much lower metabolite concentrations compared to modern metabolic networks. Research into the chemical origins of life will hence miss out on the chemistry relevant for metabolism if its focus is restricted solely to highly abundant and unreactive metabolites, including when it ignores life-compatibility of the reaction conditions as an essential constraint in enzyme evolution.

Maintenance of ATP Homeostasis Triggers Metabolic Shifts in Gas-Fermenting Acetogens

Acetogens are promising cell factories for producing fuels and chemicals from waste feedstocks via gas fermentation, but quantitative characterization of carbon, energy, and redox metabolism is required to guide their rational metabolic engineering. Here, we explore acetogen gas fermentation using physiological, metabolomics, and transcriptomics data for Clostridium autoethanogenum steady-state chemostat cultures grown on syngas at various gas-liquid mass transfer rates. We observe that C. autoethanogenum shifts from acetate to ethanol production to maintain ATP homeostasis at higher biomass concentrations but reaches a limit at a molar acetate/ethanol ratio of ∼1. This regulatory mechanism eventually leads to depletion of the intracellular acetyl-CoA pool and collapse of metabolism. We accurately predict growth phenotypes using a genome-scale metabolic model. Modeling revealed that the methylene-THF reductase reaction was ferredoxin reducing. This work provides a reference dataset to advance the understanding and engineering of arguably the first carbon fixation pathway on Earth.