Large-scale 13C-flux analysis reveals distinct transcriptional control of respiratory and fermentative metabolism in Escherichia coli

Systems Metabolic Engineering to Elucidate and Enhance Intestinal Metabolic Activities of Escherichia coli Nissle 1917

Escherichia coli Nissle 1917 (EcN) is one of the most widely used probiotics to treat gastrointestinal diseases. Recently, many studies have engineered EcN to release therapeutic proteins to treat specific diseases. However, because EcN exhibits intestinal metabolic activities, it is difficult to predict outcomes after administration. In silico and fermentation profiles revealed mucin metabolism of EcN. Multiomics revealed that fucose metabolism contributes to the intestinal colonization of EcN by enhancing the synthesis of flagella and nutrient uptake. The multiomics results also revealed that excessive intracellular trehalose synthesis in EcN, which is responsible for galactose metabolism, acts as a metabolic bottleneck, adversely affecting growth. To improve the ability of EcN to metabolize galactose, otsAB genes for trehalose synthesis were deleted, resulting in the ΔotsAB strain; the ΔotsAB strain exhibited a 1.47-fold increase in the growth rate and a 1.37-fold increase in the substrate consumption rate relative to wild-type EcN.

Highly sensitive and accurate measurement of underivatized phosphoenolpyruvate in plasma and serum via EDTA-facilitated hydrophilic interaction liquid chromatography-tandem mass spectrometry

Phosphoenolpyruvate (PEP) is an essential intermediate metabolite that is involved in various vital biochemical reactions. However, achieving the direct and accurate quantification of PEP in plasma or serum poses a significant challenge owing to its strong polarity and metal affinity. In this study, a sensitive method for the direct determination of PEP in plasma and serum based on ethylenediaminetetraacetic acid (EDTA)-facilitated hydrophilic interaction liquid chromatography-tandem mass spectrometry was developed. Superior chromatographic retention and peak shapes were achieved using a zwitterionic stationary-phase HILIC column with a metal-inert inner surface. Efficient dechelation of PEP-metal complexes in serum/plasma samples was achieved through the introduction of EDTA, resulting in a significant enhancement of the PEP signal. A PEP isotopically labelled standard was employed as a surrogate analyte for the determination of endogenous PEP, and validation assessments proved the sensitivity, selectivity, and reproducibility of this method. The method was applied to the comparative quantification of PEP in plasma and serum samples from mice and rats, as well as in HepG2 cells, HEK293T cells, and erythrocytes; the results confirmed its applicability in PEP-related biomedical research. The developed method can quantify PEP in diverse biological matrices, providing a feasible opportunity to investigate the role of PEP in relevant biomedical research.

Metabolomics and Microbial Metabolism: Toward a Systematic Understanding

Over the past decades, our understanding of microbial metabolism has increased dramatically. Metabolomics, a family of techniques that are used to measure the quantities of small molecules in biological samples, has been central to these efforts. Advances in analytical chemistry have made it possible to measure the relative and absolute concentrations of more and more compounds with increasing levels of certainty. In this review, we highlight how metabolomics has contributed to understanding microbial metabolism and in what ways it can still be deployed to expand our systematic understanding of metabolism. To that end, we explain how metabolomics was used to (a) characterize network topologies of metabolism and its regulation networks, (b) elucidate the control of metabolic function, and (c) understand the molecular basis of higher-order phenomena. We also discuss areas of inquiry where technological advances should continue to increase the impact of metabolomics, as well as areas where our understanding is bottlenecked by other factors such as the availability of statistical and modeling frameworks that can extract biological meaning from metabolomics data.

Non-consecutive enzyme interactions within TCA cycle supramolecular assembly regulate carbon-nitrogen metabolism

Enzymes of the central metabolism tend to assemble into transient supramolecular complexes. However, the functional significance of the interactions, particularly between enzymes catalyzing non-consecutive reactions, remains unclear. Here, by co-localizing two non-consecutive enzymes of the TCA cycle from Bacillus subtilis, malate dehydrogenase (MDH) and isocitrate dehydrogenase (ICD), in phase separated droplets we show that MDH-ICD interaction leads to enzyme agglomeration with a concomitant enhancement of ICD catalytic rate and an apparent sequestration of its reaction product, 2-oxoglutarate. Theory demonstrates that MDH-mediated clustering of ICD molecules explains the observed phenomena. In vivo analyses reveal that MDH overexpression leads to accumulation of 2-oxoglutarate and reduction of fluxes flowing through both the catabolic and anabolic branches of the carbon-nitrogen intersection occupied by 2-oxoglutarate, resulting in impeded ammonium assimilation and reduced biomass production. Our findings suggest that the MDH-ICD interaction is an important coordinator of carbon-nitrogen metabolism.

Metabolic Tug-of-War between Glycolysis and Translation Revealed by Biochemical Reconstitution

Inside cells, various biological systems work cooperatively for homeostasis and self-replication. These systems do not work independently as they compete for shared elements like ATP and NADH. However, it has been believed that such competition is not a problem in codependent biological systems such as the energy-supplying glycolysis and the energy-consuming translation system. In this study, we biochemically reconstituted the coupling system of glycolysis and translation using purified elements and found that the competition for ATP between glycolysis and protein synthesis interferes with their coupling. Both experiments and simulations revealed that this interference is derived from a metabolic tug-of-war between glycolysis and translation based on their reaction rates, which changes the threshold of the initial substrate concentration for the success coupling. By the metabolic tug-of-war, translation energized by strong glycolysis is facilitated by an exogenous ATPase, which normally inhibits translation. These findings provide chemical insights into the mechanism of competition among biological systems in living cells and provide a framework for the construction of synthetic metabolism in vitro.

Autotrophic growth of Escherichia coli is achieved by a small number of genetic changes

Synthetic autotrophy is a promising avenue to sustainable bioproduction from CO2. Here, we use iterative laboratory evolution to generate several distinct autotrophic strains. Utilising this genetic diversity, we identify that just three mutations are sufficient for Escherichia coli to grow autotrophically, when introduced alongside non-native energy (formate dehydrogenase) and carbon-fixing (RuBisCO, phosphoribulokinase, carbonic anhydrase) modules. The mutated genes are involved in glycolysis (pgi), central-carbon regulation (crp), and RNA transcription (rpoB). The pgi mutation reduces the enzyme's activity, thereby stabilising the carbon-fixing cycle by capping a major branching flux. For the other two mutations, we observe down-regulation of several metabolic pathways and increased expression of native genes associated with the carbon-fixing module (rpiB) and the energy module (fdoGH), as well as an increased ratio of NADH/NAD+ - the cycle's electron-donor. This study demonstrates the malleability of metabolism and its capacity to switch trophic modes using only a small number of genetic changes and could facilitate transforming other heterotrophic organisms into autotrophs.

Proteome efficiency of metabolic pathways in Escherichia coli increases along the nutrient flow

IMPORTANCE

Protein translation is the most expensive cellular process in fast-growing bacteria, and efficient proteome usage should thus be under strong natural selection. However, recent studies show that a considerable part of the proteome is unneeded for instantaneous cell growth in Escherichia coli. We still lack a systematic understanding of how this excess proteome is distributed across different pathways as a function of the growth conditions. We estimated the minimal required proteome across growth conditions in E. coli and compared the predictions with experimental data. We found that the proteome allocated to the most expensive internal pathways, including translation and the synthesis of amino acids and cofactors, is near the minimally required levels. In contrast, transporters and central carbon metabolism show much higher proteome levels than the predicted minimal abundance. Our analyses show that the proteome fraction unneeded for instantaneous cell growth decreases along the nutrient flow in E. coli.

A neural-mechanistic hybrid approach improving the predictive power of genome-scale metabolic models

Constraint-based metabolic models have been used for decades to predict the phenotype of microorganisms in different environments. However, quantitative predictions are limited unless labor-intensive measurements of media uptake fluxes are performed. We show how hybrid neural-mechanistic models can serve as an architecture for machine learning providing a way to improve phenotype predictions. We illustrate our hybrid models with growth rate predictions of Escherichia coli and Pseudomonas putida grown in different media and with phenotype predictions of gene knocked-out Escherichia coli mutants. Our neural-mechanistic models systematically outperform constraint-based models and require training set sizes orders of magnitude smaller than classical machine learning methods. Our hybrid approach opens a doorway to enhancing constraint-based modeling: instead of constraining mechanistic models with additional experimental measurements, our hybrid models grasp the power of machine learning while fulfilling mechanistic constrains, thus saving time and resources in typical systems biology or biological engineering projects.

Mitigating biomass composition uncertainties in flux balance analysis using ensemble representations

The biomass equation is a critical component in genome-scale metabolic models (GEMs): it is used as the de facto objective function in flux balance analysis (FBA). This equation accounts for the quantities of all known biomass precursors that are required for cell growth based on the macromolecular and monomer compositions measured at certain conditions. However, it is often reported that the macromolecular composition of cells could change across different environmental conditions and thus the use of the same single biomass equation in FBA, under multiple conditions, is questionable. Herein, we first investigated the qualitative and quantitative variations of macromolecular compositions of three representative host organisms, Escherichia coli, Saccharomyces cerevisiae and Cricetulus griseus, across different environmental/genetic variations. While macromolecular building blocks such as RNA, protein, and lipid composition vary notably, changes in fundamental biomass monomer units such as nucleotides and amino acids are not appreciable. We also observed that flux predictions through FBA is quite sensitive to macromolecular compositions but not the monomer compositions. Based on these observations, we propose ensemble representations of biomass equation in FBA to account for the natural variation of cellular constituents. Such ensemble representations of biomass better predicted the flux through anabolic reactions as it allows for the flexibility in the biosynthetic demands of the cells. The current study clearly highlights that certain component of the biomass equation indeed vary across different conditions, and the ensemble representation of biomass equation in FBA by accounting for such natural variations could avoid inaccuracies that may arise from in silico simulations.

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.

Complex regulation in a Comamonas platform for diverse aromatic carbon metabolism

Critical to a sustainable energy future are microbial platforms that can process aromatic carbons from the largely untapped reservoir of lignin and plastic feedstocks. Comamonas species present promising bacterial candidates for such platforms because they can use a range of natural and xenobiotic aromatic compounds and often possess innate genetic constraints that avoid competition with sugars. However, the metabolic reactions of these species are underexplored, and the regulatory mechanisms are unknown. Here we identify multilevel regulation in the conversion of lignin-related natural aromatic compounds, 4-hydroxybenzoate and vanillate, and the plastics-related xenobiotic aromatic compound, terephthalate, in Comamonas testosteroni KF-1. Transcription-level regulation controls initial catabolism and cleavage, but metabolite-level thermodynamic regulation governs fluxes in central carbon metabolism. Quantitative 13C mapping of tricarboxylic acid cycle and cataplerotic reactions elucidates key carbon routing not evident from enzyme abundance changes. This scheme of transcriptional activation coupled with metabolic fine-tuning challenges outcome predictions during metabolic manipulations.

Regulation of phosphate starvation-specific responses in Escherichia coli

Toxic agents added into the medium of rapidly growing Escherichia coli induce specific stress responses through the activation of specialized transcription factors. Each transcription factor and downstream regulon (e.g. SoxR) are linked to a unique stress (e.g. superoxide stress). Cells starved of phosphate induce several specific stress regulons during the transition to stationary phase when the growth rate is steadily declining. Whereas the regulatory cascades leading to the expression of specific stress regulons are well known in rapidly growing cells stressed by toxic products, they are poorly understood in cells starved of phosphate. The intent of this review is to both describe the unique mechanisms of activation of specialized transcription factors and discuss signalling cascades leading to the induction of specific stress regulons in phosphate-starved cells. Finally, I discuss unique defence mechanisms that could be induced in cells starved of ammonium and glucose.

The integration host factor regulates multiple virulence pathways in bacterial pathogen Dickeya zeae MS2

Dickeya zeae is an aggressive bacterial phytopathogen that infects a wide range of host plants. It has been reported that integration host factor (IHF), a nucleoid-associated protein consisting of IHFα and IHFβ subunits, regulates gene expression by influencing nucleoid structure and DNA bending. To define the role of IHF in the pathogenesis of D. zeae MS2, we deleted either and both of the IHF subunit encoding genes ihfA and ihfB, which significantly reduced the production of cell wall-degrading enzymes (CWDEs), an unknown novel phytotoxin and the virulence factor-modulating (VFM) quorum-sensing (QS) signal, cell motility, biofilm formation, and thereafter the infection ability towards both potato slices and banana seedlings. To characterize the regulatory pathways of IHF protein associated with virulence, IHF binding sites (consensus sequence 5'-WATCAANNNNTTR-3') were predicted and 272 binding sites were found throughout the genome. The expression of 110 tested genes was affected by IHF. Electrophoretic mobility shift assay (EMSA) showed direct interaction of IhfA protein with the promoters of vfmE, speA, pipR, fis, slyA, prtD, hrpL, hecB, hcp, indA, hdaA, flhD, pilT, gcpJ, arcA, arcB, and lysR. This study clarified the contribution of IHF in the pathogenic process of D. zeae by controlling the production of VFM and putrescine QS signals, phytotoxin, and indigoidine, the luxR-solo system, Fis, SlyA, and FlhD transcriptional regulators, and secretion systems from type I to type VI. Characterization of the regulatory networks of IHF in D. zeae provides a target for prevention and control of plant soft rot disease.

Two key Geobacter species of wastewater-enriched electroactive biofilm respond differently to electric field

Electroactive biofilms have attracted increasing attention due to their unique ability to exchange electrons with electrodes. Geobacter spp. are widely found to be dominant in biofilms in acetate-rich environments when an appropriate voltage is applied, but it is still largely unknown how these bacteria are selectively enriched. Herein, two key Geobacter spp. that have been demonstrated predominant in wastewater-enriched electroactive biofilm after long-term operation, G. sulfurreducens and G. anodireducens, responded to electric field (EF) differently, leading to a higher abundance of EF-sensitive G. anodireducens in the strong EF region after cocultivation with G. sulfurreducens. Transcriptome analysis indicated that two-component systems containing sensor histidine kinases and response regulators were the key for EF sensing in G. anodireducens rather than in G. sulfurreducens, which are closely connected to chemotaxis, c-di-GMP, fatty acid metabolism, pilus, oxidative phosphorylation and transcription, resulting in an increase in extracellular polymeric substance secretion and rapid cell proliferation. Our data reveal the mechanism by which EF select specific Geobacter spp. over time, providing new insights into Geobacter biofilm formation regulated by electricity.

mfapy: An open-source Python package for 13C-based metabolic flux analysis

¹³C-based metabolic flux analysis (¹³C-MFA) is an essential tool for estimating intracellular metabolic flux levels in metabolic engineering and biology. In ¹³C-MFA, a metabolic flux distribution that explains the observed isotope labeling data was computationally estimated using a non-linear optimization method. Herein, we report the development of mfapy, an open-source Python package developed for more flexibility and extensibility for ¹³C-MFA. mfapy compels users to write a customized Python code by describing each step in the data analysis procedures of the isotope labeling experiments. The flexibility and extensibility provided by mfapy can support trial-and-error performance in the routine estimation of metabolic flux distributions, experimental design by computer simulations of ¹³C-MFA experiments, and development of new data analysis techniques for stable isotope labeling experiments. mfapy is available to the public from the Github repository (https://github.com/fumiomatsuda/mfapy).

•

An open-source Python package, mfapy, is developed for ¹³C-MFA.

•

mfapy enables users to write Python codes for data analysis procedures of ¹³C-MFA.

•

mfapy has a flexibility and extensibility to support various data analysis procedures.

•

Computer simulations of ¹³C-MFA experiments is supported for experimental design.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism

To reproduce, prokaryotic viruses must hijack the cellular machinery of their hosts and redirect it toward the production of viral particles. While takeover of the host replication and protein synthesis apparatus has long been considered an essential feature of infection, recent studies indicate that extensive reprogramming of host primary metabolism is a widespread phenomenon among prokaryotic viruses that is required to fulfill the biosynthetic needs of virion production. In this review we provide an overview of the most significant recent findings regarding virus-induced reprogramming of prokaryotic metabolism and suggest how quantitative systems biology approaches may be used to provide a holistic understanding of metabolic remodeling during lytic viral infection. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Acetate overflow metabolism regulates a major metabolic shift after glucose depletion in Escherichia coli

Acetate overflow refers to the metabolism by which a large part of carbon incorporated as glucose into Escherichia coli cells is catabolized and excreted as acetate into the medium. We previously found that mutants for the acetate overflow pathway enzymes phosphoacetyltransferase (Pta) and acetate kinase (AckA) showed significant diauxic growth after glucose depletion in E. coli. Here, we analyzed the underlying mechanism in the pta mutant. Proteomic and other analyses revealed an increase in pyruvate dehydrogenase complex subunits and a decrease in glyoxylate shunt enzymes, which resulted from pyruvate accumulation. Since restoration of these enzyme levels by overexpressing PdhR (pyruvate-sensing transcription factor) or deleting iclR (gene encoding a pyruvate- and glyoxylate-sensing transcription factor) alleviated the growth lag of the pta mutant after glucose depletion, these changes were considered as the reason for the phenotype. Given the evidence for decreased coenzyme A (HS-CoA) levels in the pta mutant, the growth inhibition after glucose depletion was partly explained by limited availability of HS-CoA in the cell. The findings provide insights into the role of acetate overflow in metabolic regulation, which may be useful for biotechnological applications.

Pareto optimality between growth-rate and lag-time couples metabolic noise to phenotypic heterogeneity in Escherichia coli

Despite mounting evidence that in clonal bacterial populations, phenotypic variability originates from stochasticity in gene expression, little is known about noise-shaping evolutionary forces and how expression noise translates to phenotypic differences. Here we developed a high-throughput assay that uses a redox-sensitive dye to couple growth of thousands of bacterial colonies to their respiratory activity and show that in Escherichia coli, noisy regulation of lower glycolysis and citric acid cycle is responsible for large variations in respiratory metabolism. We found that these variations are Pareto optimal to maximization of growth rate and minimization of lag time, two objectives competing between fermentative and respiratory metabolism. Metabolome-based analysis revealed the role of respiratory metabolism in preventing the accumulation of toxic intermediates of branched chain amino acid biosynthesis, thereby supporting early onset of cell growth after carbon starvation. We propose that optimal metabolic tradeoffs play a key role in shaping and preserving phenotypic heterogeneity and adaptation to fluctuating environments.

Exploring the Glucose Fluxotype of the E. coli y-ome Using High-Resolution Fluxomics

We have developed a robust workflow to measure high-resolution fluxotypes (metabolic flux phenotypes) for large strain libraries under fully controlled growth conditions. This was achieved by optimizing and automating the whole high-throughput fluxomics process and integrating all relevant software tools. This workflow allowed us to obtain highly detailed maps of carbon fluxes in the central carbon metabolism in a fully automated manner. It was applied to investigate the glucose fluxotypes of 180 Escherichia coli strains deleted for y-genes. Since the products of these y-genes potentially play a role in a variety of metabolic processes, the experiments were designed to be agnostic as to their potential metabolic impact. The obtained data highlight the robustness of E. coli’s central metabolism to y-gene deletion. For two y-genes, deletion resulted in significant changes in carbon and energy fluxes, demonstrating the involvement of the corresponding y-gene products in metabolic function or regulation. This work also introduces novel metrics to measure the actual scope and quality of high-throughput fluxomics investigations.

Global coordination of metabolic pathways in Escherichia coli by active and passive regulation

Microorganisms adjust metabolic activity to cope with diverse environments. While many studies have provided insights into how individual pathways are regulated, the mechanisms that give rise to coordinated metabolic responses are poorly understood. Here, we identify the regulatory mechanisms that coordinate catabolism and anabolism in Escherichia coli. Integrating protein, metabolite, and flux changes in genetically implemented catabolic or anabolic limitations, we show that combined global and local mechanisms coordinate the response to metabolic limitations. To allocate proteomic resources between catabolism and anabolism, E. coli uses a simple global gene regulatory program. Surprisingly, this program is largely implemented by a single transcription factor, Crp, which directly activates the expression of catabolic enzymes and indirectly reduces the expression of anabolic enzymes by passively sequestering cellular resources needed for their synthesis. However, metabolic fluxes are not controlled by this regulatory program alone; instead, fluxes are adjusted mostly through passive changes in the local metabolite concentrations. These mechanisms constitute a simple but effective global regulatory program that coarsely partitions resources between different parts of metabolism while ensuring robust coordination of individual metabolic reactions.

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.

Proteome Regulation Patterns Determine Escherichia coli Wild-Type and Mutant Phenotypes

It is generally recognized that proteins constitute the key cellular component in shaping microbial phenotypes. Due to limited cellular resources and space, optimal allocation of proteins is crucial for microbes to facilitate maximum proliferation rates while allowing a flexible response to environmental changes. To account for the growth condition-dependent proteome in the constraint-based metabolic modeling of Escherichia coli, we consolidated a coarse-grained protein allocation approach with the explicit consideration of enzymatic constraints on reaction fluxes. Besides representing physiologically relevant wild-type phenotypes and flux distributions, the resulting protein allocation model (PAM) advances the predictability of the metabolic responses to genetic perturbations. A main driver of mutant phenotypes was ascribed to inherited regulation patterns in protein distribution among metabolic enzymes. Moreover, the PAM correctly reflected metabolic responses to an augmented protein burden imposed by the heterologous expression of green fluorescent protein. In summary, we were able to model the effects of important and frequently applied metabolic engineering approaches on microbial metabolism. Therefore, we want to promote the integration of protein allocation constraints into classical constraint-based models to foster their predictive capabilities and application for strain analysis and engineering purposes.

IMPORTANCE Predictive metabolic models are important, e.g., for generating biological knowledge and designing microbes with superior performance for target compound production. Yet today’s whole-cell models either show insufficient predictive capabilities or are computationally too expensive to be applied to metabolic engineering purposes. By linking the inherent genotype-phenotype relationship to a complete representation of the proteome, the PAM advances the accuracy of simulated phenotypes and intracellular flux distributions of E. coli. Being equally computationally lightweight as classical stoichiometric models and allowing for the application of established in silico tools, the PAM and related simulation approaches will foster the use of a model-driven metabolic research. Applications range from the investigation of mechanisms of microbial evolution to the determination of optimal strain design strategies in metabolic engineering, thus supporting basic scientists and engineers alike.

The nucleoid-associated protein IHF acts as a 'transcriptional domainin' protein coordinating the bacterial virulence traits with global transcription

Bacterial pathogenic growth requires a swift coordination of pathogenicity function with various kinds of environmental stress encountered in the course of host infection. Among the factors critical for bacterial adaptation are changes of DNA topology and binding effects of nucleoid-associated proteins transducing the environmental signals to the chromosome and coordinating the global transcriptional response to stress. In this study, we use the model phytopathogen Dickeya dadantii to analyse the organisation of transcription by the nucleoid-associated heterodimeric protein IHF. We inactivated the IHFα subunit of IHF thus precluding the IHFαβ heterodimer formation and determined both phenotypic effects of ihfA mutation on D. dadantii virulence and the transcriptional response under various conditions of growth. We show that ihfA mutation reorganises the genomic expression by modulating the distribution of chromosomal DNA supercoils at different length scales, thus affecting many virulence genes involved in both symptomatic and asymptomatic phases of infection, including those required for pectin catabolism. Altogether, we propose that IHF heterodimer is a 'transcriptional domainin' protein, the lack of which impairs the spatiotemporal organisation of transcriptional stress-response domains harbouring various virulence traits, thus abrogating the pathogenicity of D. dadantii.

Metabolic flux configuration determination using information entropy

Constraint-based models use steady-state mass balances to define a solution space of flux configurations, which can be narrowed down by measuring as many fluxes as possible. Due to loops and redundant pathways, this process typically yields multiple alternative solutions. To address this ambiguity, flux sampling can estimate the probability distribution of each flux, or a flux configuration can be singled out by further minimizing the sum of fluxes according to the assumption that cellular metabolism favors states where enzyme-related costs are economized. However, flux sampling is susceptible to artifacts introduced by thermodynamically infeasible cycles and is it not clear if the economy of fluxes assumption (EFA) is universally valid. Here, we formulated a constraint-based approach, MaxEnt, based on the principle of maximum entropy, which in this context states that if more than one flux configuration is consistent with a set of experimentally measured fluxes, then the one with the minimum amount of unwarranted assumptions corresponds to the best estimation of the non-observed fluxes. We compared MaxEnt predictions to Escherichia coli and Saccharomyces cerevisiae publicly available flux data. We found that the mean square error (MSE) between experimental and predicted fluxes by MaxEnt and EFA-based methods are three orders of magnitude lower than the median of 1,350,000 MSE values obtained using flux sampling. However, only MaxEnt and flux sampling correctly predicted flux through E. coli’s glyoxylate cycle, whereas EFA-based methods, in general, predict no flux cycles. We also tested MaxEnt predictions at increasing levels of overflow metabolism. We found that MaxEnt accuracy is not affected by overflow metabolism levels, whereas the EFA-based methods show a decreasing performance. These results suggest that MaxEnt is less sensitive than flux sampling to artifacts introduced by thermodynamically infeasible cycles and that its predictions are less susceptible to overfitting than EFA-based methods.

A guide to metabolic flux analysis in metabolic engineering: Methods, tools and applications

The field of metabolic engineering is primarily concerned with improving the biological production of value-added chemicals, fuels and pharmaceuticals through the design, construction and optimization of metabolic pathways, redirection of intracellular fluxes, and refinement of cellular properties relevant for industrial bioprocess implementation. Metabolic network models and metabolic fluxes are central concepts in metabolic engineering, as was emphasized in the first paper published in this journal, "Metabolic fluxes and metabolic engineering" (Metabolic Engineering, 1: 1-11, 1999). In the past two decades, a wide range of computational, analytical and experimental approaches have been developed to interrogate the capabilities of biological systems through analysis of metabolic network models using techniques such as flux balance analysis (FBA), and quantify metabolic fluxes using constrained-based modeling approaches such as metabolic flux analysis (MFA) and more advanced experimental techniques based on the use of stable-isotope tracers, i.e. ¹³C-metabolic flux analysis (¹³C-MFA). In this review, we describe the basic principles of metabolic flux analysis, discuss current best practices in flux quantification, highlight potential pitfalls and alternative approaches in the application of these tools, and give a broad overview of pragmatic applications of flux analysis in metabolic engineering practice.

Diauxie and co-utilization of carbon sources can coexist during bacterial growth in nutritionally complex environments

It is commonly thought that when multiple carbon sources are available, bacteria metabolize them either sequentially (diauxic growth) or simultaneously (co-utilization). However, this view is mainly based on analyses in relatively simple laboratory settings. Here we show that a heterotrophic marine bacterium, Pseudoalteromonas haloplanktis, can use both strategies simultaneously when multiple possible nutrients are provided in the same growth experiment. The order of nutrient uptake is partially determined by the biomass yield that can be achieved when the same compounds are provided as single carbon sources. Using transcriptomics and time-resolved intracellular ¹H-¹³C NMR, we reveal specific pathways for utilization of various amino acids. Finally, theoretical modelling indicates that this metabolic phenotype, combining diauxie and co-utilization of substrates, is compatible with a tight regulation that allows the modulation of assimilatory pathways.

It is thought that when multiple carbon sources are available, bacteria metabolize them either sequentially or simultaneously. Here, the authors show that a marine bacterium can use a mixed strategy when multiple possible nutrients are provided, and analyse the metabolic pathways involved.

Engineering microbial chemical factories using metabolic models

Living organisms in analogy with chemical factories use simple molecules such as sugars to produce a variety of compounds which are necessary for sustaining life and some of which are also commercially valuable. The metabolisms of simple (such as bacteria) and higher organisms (such as plants) alike can be exploited to convert low value inputs into high value outputs. Unlike conventional chemical factories, microbial production chassis are not necessarily tuned for a single product overproduction. Despite the same end goal, metabolic and industrial engineers rely on different techniques for achieving productivity goals. Metabolic engineers cannot affect reaction rates by manipulating pressure and temperature, instead they have at their disposal a range of enzymes and transcriptional and translational processes to optimize accordingly. In this review, we first highlight how various analytical approaches used in metabolic engineering and synthetic biology are related to concepts developed in systems and control engineering. Specifically, how algorithmic concepts derived in operations research can help explain the structure and organization of metabolic networks. Finally, we consider the future directions and challenges faced by the field of metabolic network modeling and the possible contributions of concepts drawn from the classical fields of chemical and control engineering. The aim of the review is to offer a current perspective of metabolic engineering and all that it entails without requiring specialized knowledge of bioinformatics or systems biology.

Automated generation of bacterial resource allocation models

Resource Balance Analysis (RBA) is a computational method based on resource allocation, which performs accurate quantitative predictions of whole-cell states (i.e. growth rate, metabolic fluxes, abundances of molecular machines including enzymes) across growth conditions. We present an integrated workflow of RBA together with the Python package RBApy. RBApy builds bacterial RBA models from annotated genome-scale metabolic models by adding descriptions of cellular processes relevant for growth and maintenance. The package includes functions for model simulation and calibration and for interfacing to Escher maps and Proteomaps for visualization. We demonstrate that RBApy faithfully reproduces results obtained by a hand-curated and experimentally validated RBA model for Bacillus subtilis. We also present a calibrated RBA model of Escherichia coli generated from scratch, which obtained excellent fits to measured flux values and enzyme abundances. RBApy makes whole-cell modelling accessible for a wide range of bacterial wild-type and engineered strains, as illustrated with a CO₂-fixing Escherichia coli strain. AVAILABILITY: RBApy is available at /https://github.com/SysBioInra/RBApy, under the licence GNU GPL version 3, and runs on Linux, Mac and Windows distributions.

BOFdat: Generating biomass objective functions for genome-scale metabolic models from experimental data

Genome-scale metabolic models (GEMs) are mathematically structured knowledge bases of metabolism that provide phenotypic predictions from genomic information. GEM-guided predictions of growth phenotypes rely on the accurate definition of a biomass objective function (BOF) that is designed to include key cellular biomass components such as the major macromolecules (DNA, RNA, proteins), lipids, coenzymes, inorganic ions and species-specific components. Despite its importance, no standardized computational platform is currently available to generate species-specific biomass objective functions in a data-driven, unbiased fashion. To fill this gap in the metabolic modeling software ecosystem, we implemented BOFdat, a Python package for the definition of a Biomass Objective Function from experimental data. BOFdat has a modular implementation that divides the BOF definition process into three independent modules defined here as steps: 1) the coefficients for major macromolecules are calculated, 2) coenzymes and inorganic ions are identified and their stoichiometric coefficients estimated, 3) the remaining species-specific metabolic biomass precursors are algorithmically extracted in an unbiased way from experimental data. We used BOFdat to reconstruct the BOF of the Escherichia coli model iML1515, a gold standard in the field. The BOF generated by BOFdat resulted in the most concordant biomass composition, growth rate, and gene essentiality prediction accuracy when compared to other methods. Installation instructions for BOFdat are available in the documentation and the source code is available on GitHub (https://github.com/jclachance/BOFdat).

Mapping Salmonella typhimurium pathways using 13C metabolic flux analysis

In the last years, Salmonella has been extensively studied not only due to its importance as a pathogen, but also as a host to produce pharmaceutical compounds. However, the full exploitation of Salmonella as a platform for bioproduct delivery has been hampered by the lack of information about its metabolism. Genome-scale metabolic models can be valuable tools to delineate metabolic engineering strategies as long as they closely represent the actual metabolism of the target organism. In the present study, a ¹³C-MFA approach was applied to map the fluxes at the central carbon pathways of S. typhimurium LT2 growing at glucose-limited chemostat cultures. The experiments were carried out in a 2L bioreactor, using defined medium enriched with 20% ¹³C-labeled glucose. Metabolic flux distributions in central carbon pathways of S. typhimurium LT2 were estimated using OpenFLUX2 based on the labeling pattern of biomass protein hydrolysates together with biomass composition. The results suggested that pentose phosphate is used to catabolize glucose, with minor fluxes through glycolysis. In silico simulations, using Optflux and pFBA as simulation method, allowed to study the performance of the genome-scale metabolic model. In general, the accuracy of in silico simulations was improved by the superimposition of estimated intracellular fluxes to the existing genome-scale metabolic model, showing a better fitting to the experimental extracellular fluxes, whereas the intracellular fluxes of pentose phosphate and anaplerotic reactions were poorly described.

13C metabolic flux analysis-guided metabolic engineering of Escherichia coli for improved acetol production from glycerol

BACKGROUND

Bioprocessing offers a sustainable and green approach to manufacture various chemicals and materials. Development of bioprocesses requires transforming common producer strains to cell factories. ¹³C metabolic flux analysis (¹³C-MFA) can be applied to identify relevant targets to accomplish the desired phenotype, which has become one of the major tools to support systems metabolic engineering. In this research, we applied ¹³C-MFA to identify bottlenecks in the bioconversion of glycerol into acetol by Escherichia coli. Valorization of glycerol, the main by-product of biodiesel, has contributed to the viability of biofuel economy.

RESULTS

We performed ¹³C-MFA and measured intracellular pyridine nucleotide pools in a first-generation acetol producer strain (HJ06) and a non-producer strain (HJ06C), and identified that engineering the NADPH regeneration is a promising target. Based on this finding, we overexpressed nadK encoding NAD kinase or pntAB encoding membrane-bound transhydrogenase either individually or in combination with HJ06, obtaining HJ06N, HJ06P and HJ06PN. The step-wise approach resulted in increasing the acetol titer from 0.91 g/L (HJ06) to 2.81 g/L (HJ06PN). To systematically characterize and the effect of mutation(s) on the metabolism, we also examined the metabolomics and transcriptional levels of key genes in four strains. The pool sizes of NADPH, NADP⁺ and the NADPH/NADP⁺ ratio were progressively increased from HJ06 to HJ06PN, demonstrating that the sufficient NADPH supply is critical for acetol production. Flux distribution was optimized towards acetol formation from HJ06 to HJ06PN: (1) The carbon partitioning at the DHAP node directed gradually more carbon from the lower glycolytic pathway through the acetol biosynthetic pathway; (2) The transhydrogenation flux was constantly increased. In addition, ¹³C-MFA showed the rigidity of upper glycolytic pathway, PP pathway and the TCA cycle to support growth. The flux patterns were supported by most metabolomics data and gene expression profiles.

CONCLUSIONS

This research demonstrated how ¹³C-MFA can be applied to drive the cycles of design, build, test and learn implementation for strain development. This succeeding engineering strategy can also be applicable for rational design of other microbial cell factories.

Model metabolic strategy for heterotrophic bacteria in the cold ocean based on Colwellia psychrerythraea 34H

Colwellia psychrerythraea 34H is a model psychrophilic bacterium found in the cold ocean-polar sediments, sea ice, and the deep sea. Although the genomes of such psychrophiles have been sequenced, their metabolic strategies at low temperature have not been quantified. We measured the metabolic fluxes and gene expression of 34H at 4 °C (the mean global-ocean temperature and a normal-growth temperature for 34H), making comparative analyses at room temperature (above its upper-growth temperature of 18 °C) and with mesophilic Escherichia coli When grown at 4 °C, 34H utilized multiple carbon substrates without catabolite repression or overflow byproducts; its anaplerotic pathways increased flux network flexibility and enabled CO₂ fixation. In glucose-only medium, the Entner-Doudoroff (ED) pathway was the primary glycolytic route; in lactate-only medium, gluconeogenesis and the glyoxylate shunt became active. In comparison, E. coli, cold stressed at 4 °C, had rapid glycolytic fluxes but no biomass synthesis. At their respective normal-growth temperatures, intracellular concentrations of TCA cycle metabolites (α-ketoglutarate, succinate, malate) were 4-17 times higher in 34H than in E. coli, while levels of energy molecules (ATP, NADH, NADPH) were 10- to 100-fold lower. Experiments with E. coli mutants supported the thermodynamic advantage of the ED pathway at cold temperature. Heat-stressed 34H at room temperature (2 hours) revealed significant down-regulation of genes associated with glycolytic enzymes and flagella, while 24 hours at room temperature caused irreversible cellular damage. We suggest that marine heterotrophic bacteria in general may rely upon simplified metabolic strategies to overcome thermodynamic constraints and thrive in the cold ocean.

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.

Effect of global transcriptional regulators on kinetic behavior of Escherichia coli under anaerobic fermentation conditions

The phenotype of Escherichia coli is governed by global transcriptional regulators under variable environmental conditions. Fnr, ArcA, IhfA-B, Crp, and Fis are amongst the major global transcription regulators that change their activity across the range of aeration, hence forming the core transcriptional network responsible for survival under changing aeration conditions. Effect of deletion of these global transcription factors on the kinetics of cell growth and mixed acid production under anaerobic fermentation conditions has not been characterized. To quantify the kinetic parameters in the absence of global transcription factors, experiments were performed using single deletion mutants of the above-mentioned global transcription regulators. The absence of global transcription regulators resulted in a relatively higher glucose uptake rate than that required for the observed growth rate. This further resulted in a higher yield of mixed acids per unit biomass in mutants as compared to the parent strain (E. coli BW25113). Thus, the increased channeling of carbon towards mixed acid secretion resulted in a lower growth rate in the mutants.

Data Reduction Approaches for Dissecting Transcriptional Effects on Metabolism

The availability of high-throughput data from transcriptomics and metabolomics technologies provides the opportunity to characterize the transcriptional effects on metabolism. Here we propose and evaluate two computational approaches rooted in data reduction techniques to identify and categorize transcriptional effects on metabolism by combining data on gene expression and metabolite levels. The approaches determine the partial correlation between two metabolite data profiles upon control of given principal components extracted from transcriptomics data profiles. Therefore, they allow us to investigate both data types with all features simultaneously without doing preselection of genes. The proposed approaches allow us to categorize the relation between pairs of metabolites as being under transcriptional or post-transcriptional regulation. The resulting classification is compared to existing literature and accumulated evidence about regulatory mechanism of reactions and pathways in the cases of Escherichia coli, Saccharomycies cerevisiae, and Arabidopsis thaliana.

A guide to 13C metabolic flux analysis for the cancer biologist

Cancer metabolism is significantly altered from normal cellular metabolism allowing cancer cells to adapt to changing microenvironments and maintain high rates of proliferation. In the past decade, stable-isotope tracing and network analysis have become powerful tools for uncovering metabolic pathways that are differentially activated in cancer cells. In particular, ¹³C metabolic flux analysis (¹³C-MFA) has emerged as the primary technique for quantifying intracellular fluxes in cancer cells. In this review, we provide a practical guide for investigators interested in getting started with ¹³C-MFA. We describe best practices in ¹³C-MFA, highlight potential pitfalls and alternative approaches, and conclude with new developments that can further enhance our understanding of cancer metabolism.

Tracing tagged molecules can help researchers understand the altered metabolism of cancer cells. The abilities of cancer cells to multiply rapidly and invade new tissues are supported by metabolic alterations, which can be investigated by feeding tagged molecules to cells and tracing how they are metabolized. These techniques, such as ¹³C metabolic flux analysis (¹³C-MFA), have been perceived as difficult to use, but recent advances are making them more accessible. Maciek Antoniewicz, University of Delaware, Newark, USA, has published a practical guide for researchers wanting to use 13C-MFA. The review includes best practices, pitfalls, alternative approaches, and new developments, especially new user-friendly software that allows researchers without extensive training in mathematics, statistics, or coding to perform 13C-MFA. Broadening access to tools for investigating altered metabolic pathways may spur development of new cancer therapies targeting these pathways.

The effect of global transcriptional regulators on the anaerobic fermentative metabolism of Escherichia coli

Global transcription factors are known to regulate the anaerobic growth of Escherichia coli on glucose. These transcription factors help the organism to sense oxygen and accordingly regulate the synthesis of mixed acid producing enzymes. Five global transcription factors, namely ArcA, Fnr, IhfA-B, Crp and Fis, are known to play an important role in the growth phenotype of the organism in the transition from anaerobic to aerobic conditions. The effect of deletion of most of these global transcription factors on the growth phenotype has not been characterized under strict anaerobic fermentation conditions. In order to enumerate the role of global transcription factors in central carbon metabolism, experiments were performed using single deletion mutants of the above mentioned global transcription regulators. The mutants demonstrated lower growth rates, ranging from 3-75% lower growth as compared to the wild-type strain along with varying glucose uptake rates. Global transcription regulators help in lowering formate and acetate synthesis, thereby effectively channeling the carbon towards redox balance (through ethanol formation) and biomass synthesis. Flux analysis of mutant strains indicated that deletion of a single transcription factor alone does not play a significant role in the normalized flux distribution of the central carbon metabolism.

Catabolite repression in Campylobacter jejuni correlates with intracellular succinate levels

Bacteria have evolved different mechanisms to catabolize carbon sources from nutrient mixtures. They first consume their preferred carbon source, before others are used. Regulatory mechanisms adapt the metabolism accordingly to maximize growth and to outcompete other organisms. The human pathogen Campylobacter jejuni is an asaccharolytic Gram-negative bacterium that catabolizes amino acids and organic acids for growth. It prefers serine and aspartate as carbon sources, however it lacks all regulators known to be involved in regulating carbon source utilization in other organisms. In which manner C. jejuni adapts its metabolism towards the presence or absence of preferred carbon sources is unknown. In this study, we show with transcriptomic analysis and enzyme assays how C. jejuni adapts its metabolism in response to its preferred carbon sources. In the presence of serine as well as lactate and pyruvate C. jejuni inhibits the utilization of other carbon sources, by repressing the expression of a number of central metabolic enzymes. The regulatory proteins RacR, Cj1000 and CsrA play a role in the regulation of these metabolic enzymes. This metabolism dependent transcriptional repression correlates with an accumulation of intracellular succinate. Hence, we propose a demand-based catabolite repression mechanism in C. jejuni, depended on intracellular succinate levels.

Metabolic enzyme cost explains variable trade-offs between microbial growth rate and yield

Microbes may maximize the number of daughter cells per time or per amount of nutrients consumed. These two strategies correspond, respectively, to the use of enzyme-efficient or substrate-efficient metabolic pathways. In reality, fast growth is often associated with wasteful, yield-inefficient metabolism, and a general thermodynamic trade-off between growth rate and biomass yield has been proposed to explain this. We studied growth rate/yield trade-offs by using a novel modeling framework, Enzyme-Flux Cost Minimization (EFCM) and by assuming that the growth rate depends directly on the enzyme investment per rate of biomass production. In a comprehensive mathematical model of core metabolism in E. coli, we screened all elementary flux modes leading to cell synthesis, characterized them by the growth rates and yields they provide, and studied the shape of the resulting rate/yield Pareto front. By varying the model parameters, we found that the rate/yield trade-off is not universal, but depends on metabolic kinetics and environmental conditions. A prominent trade-off emerges under oxygen-limited growth, where yield-inefficient pathways support a 2-to-3 times higher growth rate than yield-efficient pathways. EFCM can be widely used to predict optimal metabolic states and growth rates under varying nutrient levels, perturbations of enzyme parameters, and single or multiple gene knockouts.

Pyruvate cycle increases aminoglycoside efficacy and provides respiratory energy in bacteria

The emergence and ongoing spread of multidrug-resistant bacteria puts humans and other species at risk for potentially lethal infections. Thus, novel antibiotics or alternative approaches are needed to target drug-resistant bacteria, and metabolic modulation has been documented to improve antibiotic efficacy, but the relevant metabolic mechanisms require more studies. Here, we show that glutamate potentiates aminoglycoside antibiotics, resulting in improved elimination of antibiotic-resistant pathogens. When exploring the metabolic flux of glutamate, it was found that the enzymes that link the phosphoenolpyruvate (PEP)-pyruvate-AcCoA pathway to the TCA cycle were key players in this increased efficacy. Together, the PEP-pyruvate-AcCoA pathway and TCA cycle can be considered the pyruvate cycle (P cycle). Our results show that inhibition or gene depletion of the enzymes in the P cycle shut down the TCA cycle even in the presence of excess carbon sources, and that the P cycle operates routinely as a general mechanism for energy production and regulation in Escherichia coli and Edwardsiella tarda These findings address metabolic mechanisms of metabolite-induced potentiation and fundamental questions about bacterial biochemistry and energy metabolism.

Evolutionary adaptations to new environments generally reverse plastic phenotypic changes

Organismal adaptation to a new environment may start with plastic phenotypic changes followed by genetic changes, but whether the plastic changes are stepping stones to genetic adaptation is debated. Here we address this question by investigating gene expression and metabolic flux changes in the two-phase adaptation process using transcriptomic data from multiple experimental evolution studies and computational metabolic network analysis, respectively. We discover that genetic changes more frequently reverse than reinforce plastic phenotypic changes in virtually every adaptation. Metabolic network analysis reveals that, even in the presence of plasticity, organismal fitness drops after environmental shifts, but largely recovers through subsequent evolution. Such fitness trajectories explain why plastic phenotypic changes are genetically compensated rather than strengthened. In conclusion, although phenotypic plasticity may serve as an emergency response to a new environment that is necessary for survival, it does not generally facilitate genetic adaptation by bringing the organismal phenotype closer to the new optimum.

Phenotypic plasticity has been suggested to facilitate survival in new environments and subsequent adaptation. Here, the authors reanalyze transcriptomic data from experimental evolution studies in combination with computational metabolic network analysis and show that genetic adaptation tends to reverse plastic changes in order to recover fitness.

Why Do Fast-Growing Bacteria Enter Overflow Metabolism? Testing the Membrane Real Estate Hypothesis

Bacteria and other cells show a puzzling behavior. At high growth rates, E. coli switch from respiration (which is ATP-efficient) to using fermentation for additional ATP (which is inefficient). This overflow metabolism results in a several-fold decrease in ATP produced per glucose molecule provided as food. By integrating diverse types of experimental data into a simple biophysical model, we give evidence that this onset is the result of the membrane real estate hypothesis: Fast growth drives cells to be bigger, reducing their surface-to-volume ratios. This decreases the membrane area available for respiratory proteins despite growing demand, causing increased crowding. Only when respiratory proteins reach their crowding limit does the cell activate fermentation, since fermentation allows faster ATP production per unit membrane area. Surface limitation thus creates a Pareto trade-off between membrane efficiency and ATP yield that links metabolic choice to the size and shape of a bacterial cell. By exploring the predictions that emerge from this trade-off, we show how consideration of molecular structures, energetics, rates, and equilibria can provide important insight into cellular behavior.

redGEM: Systematic reduction and analysis of genome-scale metabolic reconstructions for development of consistent core metabolic models

Genome-scale metabolic reconstructions have proven to be valuable resources in enhancing our understanding of metabolic networks as they encapsulate all known metabolic capabilities of the organisms from genes to proteins to their functions. However the complexity of these large metabolic networks often hinders their utility in various practical applications. Although reduced models are commonly used for modeling and in integrating experimental data, they are often inconsistent across different studies and laboratories due to different criteria and detail, which can compromise transferability of the findings and also integration of experimental data from different groups. In this study, we have developed a systematic semi-automatic approach to reduce genome-scale models into core models in a consistent and logical manner focusing on the central metabolism or subsystems of interest. The method minimizes the loss of information using an approach that combines graph-based search and optimization methods. The resulting core models are shown to be able to capture key properties of the genome-scale models and preserve consistency in terms of biomass and by-product yields, flux and concentration variability and gene essentiality. The development of these “consistently-reduced” models will help to clarify and facilitate integration of different experimental data to draw new understanding that can be directly extendable to genome-scale models.

Reduced models are used commonly to understand the metabolism of organisms and to integrate experimental data for many different studies such as physiology, fluxomics and metabolomics. Without consistent or clear criteria on how these reduced models are actually developed, it is difficult to ensure that they reflect the detailed knowledge that is kept in genome scale metabolic network models (GEMs). The redGEM algorithm presented here allows us to systematically develop consistently reduced metabolic models from their genome-scale counterparts. We applied redGEM for the construction of a core model for E. coli central carbon metabolism. We constructed the core model irJO1366 based on the latest genome-scale E. coli metabolic reconstruction (iJO1366). irJO1366 contains the central carbon pathways and other immediate pathways that must be connected to them for consistency with the iJO1366. irJO1366 can be used to understand metabolism of the organism and also to provide guidance for metabolic engineering purposes. The algorithm is also designed to be modular so that heterologous reactions or pathways can be appended to the core model akin to a “plug-and-play”, synthetic biology approach. The algorithm is applicable to any compartmentalized or non-compartmentalized GEM.

Assessing glycolytic flux alterations resulting from genetic perturbations in E. coli using a biosensor

We describe the development of an optimized glycolytic flux biosensor and its application in detecting altered flux in a production strain and in a mutant library. The glycolytic flux biosensor is based on the Cra-regulated ppsA promoter of E. coli controlling fluorescent protein synthesis. We validated the glycolytic flux dependency of the biosensor in a range of different carbon sources in six different E. coli strains and during mevalonate production. Furthermore, we studied the flux-altering effects of genome-wide single gene knock-outs in E. coli in a multiplex FlowSeq experiment. From a library consisting of 2126 knock-out mutants, we identified 3 mutants with high-flux and 95 mutants with low-flux phenotypes that did not have severe growth defects. This approach can improve our understanding of glycolytic flux regulation improving metabolic models and engineering efforts.

Transcriptional reprogramming in yeast using dCas9 and combinatorial gRNA strategies

BACKGROUND

Transcriptional reprogramming is a fundamental process of living cells in order to adapt to environmental and endogenous cues. In order to allow flexible and timely control over gene expression without the interference of native gene expression machinery, a large number of studies have focused on developing synthetic biology tools for orthogonal control of transcription. Most recently, the nuclease-deficient Cas9 (dCas9) has emerged as a flexible tool for controlling activation and repression of target genes, by the simple RNA-guided positioning of dCas9 in the vicinity of the target gene transcription start site.

RESULTS

In this study we compared two different systems of dCas9-mediated transcriptional reprogramming, and applied them to genes controlling two biosynthetic pathways for biobased production of isoprenoids and triacylglycerols (TAGs) in baker's yeast Saccharomyces cerevisiae. By testing 101 guide-RNA (gRNA) structures on a total of 14 different yeast promoters, we identified the best-performing combinations based on reporter assays. Though a larger number of gRNA-promoter combinations do not perturb gene expression, some gRNAs support expression perturbations up to ~threefold. The best-performing gRNAs were used for single and multiplex reprogramming strategies for redirecting flux related to isoprenoid production and optimization of TAG profiles. From these studies, we identified both constitutive and inducible multiplex reprogramming strategies enabling significant changes in isoprenoid production and increases in TAG.

CONCLUSION

Taken together, we show similar performance for a constitutive and an inducible dCas9 approach, and identify multiplex gRNA designs that can significantly perturb isoprenoid production and TAG profiles in yeast without editing the genomic context of the target genes. We also identify a large number of gRNA positions in 14 native yeast target pomoters that do not affect expression, suggesting the need for further optimization of gRNA design tools and dCas9 engineering.

Metabolic regulation is sufficient for global and robust coordination of glucose uptake, catabolism, energy production and growth in Escherichia coli

The metabolism of microorganisms is regulated through two main mechanisms: changes of enzyme capacities as a consequence of gene expression modulation (“hierarchical control”) and changes of enzyme activities through metabolite-enzyme interactions. An increasing body of evidence indicates that hierarchical control is insufficient to explain metabolic behaviors, but the system-wide impact of metabolic regulation remains largely uncharacterized. To clarify its role, we developed and validated a detailed kinetic model of Escherichia coli central metabolism that links growth to environment. Metabolic control analyses confirm that the control is widely distributed across the network and highlight strong interconnections between all the pathways. Exploration of the model solution space reveals that several robust properties emerge from metabolic regulation, from the molecular level (e.g. homeostasis of total metabolite pool) to the overall cellular physiology (e.g. coordination of carbon uptake, catabolism, energy and redox production, and growth), while allowing a large degree of flexibility at most individual metabolic steps. These properties have important physiological implications for E. coli and significantly expand the self-regulating capacities of its metabolism.

Metabolism is a fundamental biochemical process that enables cells to operate and grow by converting nutrients into ‘building blocks’ and energy. Metabolism happens through the work of enzymes, which are encoded by genes. Thus, genes and their regulation are often thought of controlling metabolism, somewhat at the top of a hierarchical control system. However, an increasing body of evidence indicates that metabolism plays an active role in the control of its own operation via a dense network of metabolite-enzyme interactions. The system-wide role of metabolic regulation is hard to dissect and so far remains largely uncharacterized. To better understand its role, we constructed a detailed kinetic model of the carbon and energy metabolism of the bacterium Escherichia coli, a model organism in Systems and Synthetic biology. Model simulations indicate that kinetic considerations of metabolism alone can explain data from hundreds of experiments, without needing to invoke regulation of gene expression. In particular, metabolic regulation is sufficient to coordinate carbon utilization, redox and energy production, and growth, while maintaining local flexibility at individual metabolic steps. These findings indicate that the self-regulating capacities of E. coli metabolism are far more significant than previously expected, and improve our understanding on how cells work.

Few regulatory metabolites coordinate expression of central metabolic genes in Escherichia coli

Transcription networks consist of hundreds of transcription factors with thousands of often overlapping target genes. While we can reliably measure gene expression changes, we still understand relatively little why expression changes the way it does. How does a coordinated response emerge in such complex networks and how many input signals are necessary to achieve it? Here, we unravel the regulatory program of gene expression in Escherichia coli central carbon metabolism with more than 30 known transcription factors. Using a library of fluorescent transcriptional reporters, we comprehensively quantify the activity of central metabolic promoters in 26 environmental conditions. The expression patterns were dominated by growth rate‐dependent global regulation for most central metabolic promoters in concert with highly condition‐specific activation for only few promoters. Using an approximate mathematical description of promoter activity, we dissect the contribution of global and specific transcriptional regulation. About 70% of the total variance in promoter activity across conditions was explained by global transcriptional regulation. Correlating the remaining specific transcriptional regulation of each promoter with the cell's metabolome response across the same conditions identified potential regulatory metabolites. Remarkably, cyclic AMP, fructose‐1,6‐bisphosphate, and fructose‐1‐phosphate alone explained most of the specific transcriptional regulation through their interaction with the two major transcription factors Crp and Cra. Thus, a surprisingly simple regulatory program that relies on global transcriptional regulation and input from few intracellular metabolites appears to be sufficient to coordinate E. coli central metabolism and explain about 90% of the experimentally observed transcription changes in 100 genes.