A model-driven quantitative metabolomics analysis of aerobic and anaerobic metabolism in E. coli K-12 MG1655 that is biochemically and thermodynamically consistent

Model-driven experimental design workflow expands understanding of regulatory role of Nac in Escherichia coli

The establishment of experimental conditions for transcriptional regulator network (TRN) reconstruction in bacteria continues to be impeded by the limited knowledge of activating conditions for transcription factors (TFs). Here, we present a novel genome-scale model-driven workflow for designing experimental conditions, which optimally activate specific TFs. Our model-driven workflow was applied to elucidate transcriptional regulation under nitrogen limitation by Nac and NtrC, in Escherichia coli. We comprehensively predict alternative nitrogen sources, including cytosine and cytidine, which trigger differential activation of Nac using a model-driven workflow. In accordance with the prediction, genome-wide measurements with ChIP-exo and RNA-seq were performed. Integrative data analysis reveals that the Nac and NtrC regulons consist of 97 and 43 genes under alternative nitrogen conditions, respectively. Functional analysis of Nac at the transcriptional level showed that Nac directly down-regulates amino acid biosynthesis and restores expression of tricarboxylic acid (TCA) cycle genes to alleviate nitrogen-limiting stress. We also demonstrate that both TFs coherently modulate α-ketoglutarate accumulation stress due to nitrogen limitation by co-activating amino acid and diamine degradation pathways. A systems-biology approach provided a detailed and quantitative understanding of both TF's roles and how nitrogen and carbon metabolic networks respond complementarily to nitrogen-limiting stress.

Flux balance analysis of metabolic networks for efficient engineering of microbial cell factories

Metabolic engineering principles have long been applied to explore the metabolic networks of complex microbial cell factories under a variety of environmental constraints for effective deployment of the microorganisms in the optimal production of biochemicals like biofuels, polymers, amino acids, recombinant proteins. One of the methodologies used for analyzing microbial metabolic networks is the Flux Balance Analysis (FBA), which employs applications of optimization techniques for forecasting biomass growth and metabolic flux distribution of industrially important products under specified environmental conditions. The in silico flux simulations are instrumental for designing the production-specific microbial cell factories. However, FBA has some inherent limitations. The present review emphasizes how the incorporation of additional kinetic, thermodynamic, expression and regulatory constraints and integration of omics data into the classical FBA platform improve the prediction accuracy of FBA. A programmed comparison of the simulated data with the experimental observations is presented for supporting the claim. The review further accounts for the successful implementation of classical FBA in biotechnological applications and identifies areas in which classical FBA fails to make correct predictions. The analysis of the predictive capabilities of the different FBA strategies presented here is expected to help researchers in finding new avenues in engineering highly efficient microbial metabolic pathways and identify the key metabolic bottlenecks during the process. Based on the appropriate metabolic network design, fermentation engineers will be able to effectively design the bioreactors and optimize large-scale biochemical production through suitable pathway modifications.

Constraint-based metabolic control analysis for rational strain engineering

The advancements in genome editing techniques over the past years have rekindled interest in rational metabolic engineering strategies. While Metabolic Control Analysis (MCA) is a well-established method for quantifying the effects of metabolic engineering interventions on flows in metabolic networks and metabolite concentrations, it does not consider the physiological limitations of the cellular environment and metabolic engineering design constraints. We report here a constraint-based framework, Network Response Analysis (NRA), for rational genetic strain design. NRA is cast as a Mixed-Integer Linear Programming problem that integrates MCA, Thermodynamically-based Flux Analysis (TFA), biologically relevant constraints, as well as genome editing restrictions into a comprehensive platform for identifying metabolic engineering targets. We show that the NRA formulation and its core constraints are equivalent to the ones of Flux Balance Analysis (FBA) and TFA, which allows it to be used for a wide range of optimization criteria and with various physiological constraints. We also show how the parametrization and introduction of biological constraints enhance the NRA formulation compared to the classical MCA approach, and we demonstrate its features and its ability to generate multiple alternative optimal strategies given several user-defined boundaries and objectives. In summary, NRA is a sophisticated alternative to classical MCA for rational metabolic engineering that accommodates the incorporation of physiological data at metabolic flux, metabolite concentration, and enzyme expression levels.

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.

The effects of model complexity and size on metabolic flux distribution and control: case study in Escherichia coli

Background

Significant efforts have been made in building large-scale kinetic models of cellular metabolism in the past two decades. However, most kinetic models published to date, remain focused around central carbon pathways or are built around ad hoc reduced models without clear justification on their derivation and usage. Systematic algorithms exist for reducing genome-scale metabolic reconstructions to build thermodynamically feasible and consistently reduced stoichiometric models. However, it is important to study how network complexity affects conclusions derived from large-scale kinetic models built around consistently reduced models before we can apply them to study biological systems.

Results

We reduced the iJO1366 Escherichia Coli genome-scale metabolic reconstruction systematically to build three stoichiometric models of different size. Since the reduced models are expansions around the core subsystems for which the reduction was performed, the models are nested. We present a method for scaling up the flux profile and the concentration vector reference steady-states from the smallest model to the larger ones, whilst preserving maximum equivalency. Populations of kinetic models, preserving similarity in kinetic parameters, were built around the reference steady-states and their metabolic sensitivity coefficients (MSCs) were computed. The MSCs were sensitive to the model complexity. We proposed a metric for measuring the sensitivity of MSCs to these structural changes.

Conclusions

We proposed for the first time a workflow for scaling up the size of kinetic models while preserving equivalency between the kinetic models. Using this workflow, we demonstrate that model complexity in terms of networks size has significant impact on sensitivity characteristics of kinetic models. Therefore, it is essential to account for the effects of network complexity when constructing kinetic models. The presented metric for measuring MSC sensitivity to structural changes can guide modelers and experimentalists in improving model quality and guide synthetic biology and metabolic engineering. Our proposed workflow enables the testing of the suitability of a kinetic model for answering certain study-specific questions. We argue that the model-based metabolic design targets that are common across models of different size are of higher confidence, while those that are different could be the objective of investigations for model improvement.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04066-y.

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.

Principles and practice of designing microbial biocatalysts for fuel and chemical production

The finite nature of fossil fuels and the environmental impact of its use have raised interest in alternate renewable energy sources. Specifically, non-food carbohydrates, such as lignocellulosic biomass, can be used to produce next generation biofuels, including cellulosic ethanol and other non-ethanol fuels like butanol. However, currently there is no native microorganism that can ferment all lignocellulosic sugars to fuel molecules. Thus, research is focused on engineering improved microbial biocatalysts for production of liquid fuels at high productivity, titer and yield. A clear understanding and application of the basic principles of microbial physiology and biochemistry are crucial to achieve this goal. In this review, we present and discuss the construction of microbial biocatalysts that integrate these principles with ethanol-producing Escherichia coli as an example of metabolic engineering. These principles also apply to fermentation of lignocellulosic sugars to other chemicals that are currently produced from petroleum.

Clavulanic Acid Production by Streptomyces clavuligerus: Insights from Systems Biology, Strain Engineering, and Downstream Processing

Clavulanic acid (CA) is an irreversible β-lactamase enzyme inhibitor with a weak antibacterial activity produced by Streptomyces clavuligerus (S. clavuligerus). CA is typically co-formulated with broad-spectrum β‑lactam antibiotics such as amoxicillin, conferring them high potential to treat diseases caused by bacteria that possess β‑lactam resistance. The clinical importance of CA and the complexity of the production process motivate improvements from an interdisciplinary standpoint by integrating metabolic engineering strategies and knowledge on metabolic and regulatory events through systems biology and multi-omics approaches. In the large-scale bioprocessing, optimization of culture conditions, bioreactor design, agitation regime, as well as advances in CA separation and purification are required to improve the cost structure associated to CA production. This review presents the recent insights in CA production by S. clavuligerus, emphasizing on systems biology approaches, strain engineering, and downstream processing.

Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis

Metabolic networks are regulated to ensure the dynamic adaptation of biochemical reaction fluxes to maintain cell homeostasis and optimal metabolic fitness in response to endogenous and exogenous perturbations. To this end, metabolism is tightly controlled by dynamic and intricate regulatory mechanisms involving allostery, enzyme abundance and post-translational modifications. The study of the molecular entities involved in these complex mechanisms has been boosted by the advent of high-throughput technologies. The so-called omics enable the quantification of the different molecular entities at different system layers, connecting the genotype with the phenotype. Therefore, the study of the overall behavior of a metabolic network and the omics data integration and analysis must be approached from a holistic perspective. Due to the close relationship between metabolism and cellular phenotype, metabolic modelling has emerged as a valuable tool to decipher the underlying mechanisms governing cell phenotype. Constraint-based modelling and kinetic modelling are among the most widely used methods to study cell metabolism at different scales, ranging from cells to tissues and organisms. These approaches enable integrating metabolomic data, among others, to enhance model predictive capabilities. In this review, we describe the current state of the art in metabolic modelling and discuss future perspectives and current challenges in the field.

The ETFL formulation allows multi-omics integration in thermodynamics-compliant metabolism and expression models

Systems biology has long been interested in models capturing both metabolism and expression in a cell. We propose here an implementation of the metabolism and expression model formalism (ME-models), which we call ETFL, for Expression and Thermodynamics Flux models. ETFL is a hierarchical model formulation, from metabolism to RNA synthesis, that allows simulating thermodynamics-compliant intracellular fluxes as well as enzyme and mRNA concentration levels. ETFL formulates a mixed-integer linear problem (MILP) that enables both relative and absolute metabolite, protein, and mRNA concentration integration. ETFL is compatible with standard MILP solvers and does not require a non-linear solver, unlike the previous state of the art. It also accounts for growth-dependent parameters, such as relative protein or mRNA content. We present ETFL along with its validation using results obtained from a well-characterized E. coli model. We show that ETFL is able to reproduce proteome-limited growth. We also subject it to several analyses, including the prediction of feasible mRNA and enzyme concentrations and gene essentiality.

Accounting for the effects of genetic expression in genome-scale metabolic models is challenging. Here, the authors introduce a model formulation that efficiently simulates thermodynamic-compliant fluxes, enzyme and mRNA concentration levels, allowing omics integration and broad analysis of in silico cellular physiology.

Energy metabolism controls phenotypes by protein efficiency and allocation

Cells require energy for growth and maintenance and have evolved to have multiple pathways to produce energy in response to varying conditions. A basic question in this context is how cells organize energy metabolism, which is, however, challenging to elucidate due to its complexity, i.e., the energy-producing pathways overlap with each other and even intertwine with biomass formation pathways. Here, we propose a modeling concept that decomposes energy metabolism into biomass formation and ATP-producing pathways. The latter can be further decomposed into a high-yield and a low-yield pathway. This enables independent estimation of protein efficiency for each pathway. With this concept, we modeled energy metabolism for Escherichia coli and Saccharomyces cerevisiae and found that the high-yield pathway shows lower protein efficiency than the low-yield pathway. Taken together with a fixed protein constraint, we predict overflow metabolism in E. coli and the Crabtree effect in S. cerevisiae, meaning that energy metabolism is sufficient to explain the metabolic switches. The static protein constraint is supported by the findings that protein mass of energy metabolism is conserved across conditions based on absolute proteomics data. This also suggests that enzymes may have decreased saturation or activity at low glucose uptake rates. Finally, our analyses point out three ways to improve growth, i.e., increasing protein allocation to energy metabolism, decreasing ATP demand, or increasing activity for key enzymes.

A White-Box Machine Learning Approach for Revealing Antibiotic Mechanisms of Action

Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop an integrated "white-box" biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genome-scale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.

Kinetic models of metabolism that consider alternative steady-state solutions of intracellular fluxes and concentrations

Large-scale kinetic models are used for designing, predicting, and understanding the metabolic responses of living cells. Kinetic models are particularly attractive for the biosynthesis of target molecules in cells as they are typically better than other types of models at capturing the complex cellular biochemistry. Using simpler stoichiometric models as scaffolds, kinetic models are built around a steady-state flux profile and a metabolite concentration vector that are typically determined via optimization. However, as the underlying optimization problem is underdetermined, even after incorporating available experimental omics data, one cannot uniquely determine the operational configuration in terms of metabolic fluxes and metabolite concentrations. As a result, some reactions can operate in either the forward or reverse direction while still agreeing with the observed physiology. Here, we analyze how the underlying uncertainty in intracellular fluxes and concentrations affects predictions of constructed kinetic models and their design in metabolic engineering and systems biology studies. To this end, we integrated the omics data of optimally grown Escherichia coli into a stoichiometric model and constructed populations of non-linear large-scale kinetic models of alternative steady-state solutions consistent with the physiology of the E. coli aerobic metabolism. We performed metabolic control analysis (MCA) on these models, highlighting that MCA-based metabolic engineering decisions are strongly affected by the selected steady state and appear to be more sensitive to concentration values rather than flux values. To incorporate this into future studies, we propose a workflow for moving towards more reliable and robust predictions that are consistent with all alternative steady-state solutions. This workflow can be applied to all kinetic models to improve the consistency and accuracy of their predictions. Additionally, we show that, irrespective of the alternative steady-state solution, increased activity of phosphofructokinase and decreased ATP maintenance requirements would improve cellular growth of optimally grown E. coli.

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.

Multiple Optimal Phenotypes Overcome Redox and Glycolytic Intermediate Metabolite Imbalances in Escherichia coli pgi Knockout Evolutions

A mechanistic understanding of how new phenotypes develop to overcome the loss of a gene product provides valuable insight on both the metabolic and regulatory functions of the lost gene. The pgi gene, whose product catalyzes the second step in glycolysis, was deleted in a growth-optimized Escherichia coli K-12 MG1655 strain. The initial knockout (KO) strain exhibited an 80% drop in growth rate that was largely recovered in eight replicate, but phenotypically distinct, cultures after undergoing adaptive laboratory evolution (ALE). Multi-omic data sets showed that the loss of pgi substantially shifted pathway usage, leading to a redox and sugar phosphate stress response. These stress responses were overcome by unique combinations of innovative mutations selected for by ALE. Thus, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.IMPORTANCE A mechanistic understanding of how microbes are able to overcome the loss of a gene through regulatory and metabolic changes is not well understood. Eight independent adaptive laboratory evolution (ALE) experiments with pgi knockout strains resulted in eight phenotypically distinct endpoints that were able to overcome the gene loss. Utilizing multi-omics analysis, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.

Metabolomics analysis of Pseudomonas chlororaphis JK12 algicidal activity under aerobic and micro-aerobic culture condition

Utilization of algicidal bacteria as a biological agent have been receiving significant interest for controlling harmful algal blooms. While various algicidal bacterial strains have been identified, limited studies have explored the influence of bacterial culture conditions on its algicidal activity. Here, the effect of oxygen on the algicidal activity of a novel bacterium JK12, against a model diatom, Phaeodactylum tricornutum (P. tricornutum) was studied. Strain JK12 showed high algicidal activity against P. tricornutum and was identified as Pseudomonas chlororaphis (P. chlororaphis) by 16S ribosomal RNA gene analysis. JK12 culture supernatant exhibited strong algicidal activity while washed JK12 cells showed no obvious activity, indicating that JK12 indirectly attacks algae by secreting extracellular algicidal metabolites. Micro-aerobic culture condition dramatically enhanced the algicidal activity of JK12 by 50%, compared to that cultured under aerobic condition in 24 h. Extracellular metabolomic profiling of JK12 using gas chromatography–mass spectrometry and liquid chromatography–mass spectrometry analysis revealed significantly higher amounts of allantoic acid, urocanic acid, cytidine 2′,3′-cyclic phosphate, uridine 2′,3′-cyclic phosphate, and chlorinated tryptophan in the micro-aerobic culture. This is the first report to demonstrate the important role of oxygen on the algicidal activity of a non-pathogenic strain P. chlororaphis. In addition, the metabolomics analysis provided insights into the algicidal mechanism of P. chlororaphis.

Metabolome analysis-based design and engineering of a metabolic pathway in Corynebacterium glutamicum to match rates of simultaneous utilization of D-glucose and L-arabinose

Background

l-Arabinose is the second most abundant component of hemicellulose in lignocellulosic biomass, next to d-xylose. However, few microorganisms are capable of utilizing pentoses, and catabolic genes and operons enabling bacterial utilization of pentoses are typically subject to carbon catabolite repression by more-preferred carbon sources, such as d-glucose, leading to a preferential utilization of d-glucose over pentoses. In order to simultaneously utilize both d-glucose and l-arabinose at the same rate, a modified metabolic pathway was rationally designed based on metabolome analysis.

Results

Corynebacterium glutamicum ATCC 31831 utilized d-glucose and l-arabinose simultaneously at a low concentration (3.6 g/L each) but preferentially utilized d-glucose over l-arabinose at a high concentration (15 g/L each), although l-arabinose and d-glucose were consumed at comparable rates in the absence of the second carbon source. Metabolome analysis revealed that phosphofructokinase and pyruvate kinase were major bottlenecks for d-glucose and l-arabinose metabolism, respectively. Based on the results of metabolome analysis, a metabolic pathway was engineered by overexpressing pyruvate kinase in combination with deletion of araR, which encodes a repressor of l-arabinose uptake and catabolism. The recombinant strain utilized high concentrations of d-glucose and l-arabinose (15 g/L each) at the same consumption rate. During simultaneous utilization of both carbon sources at high concentrations, intracellular levels of phosphoenolpyruvate declined and acetyl-CoA levels increased significantly as compared with the wild-type strain that preferentially utilized d-glucose. These results suggest that overexpression of pyruvate kinase in the araR deletion strain increased the specific consumption rate of l-arabinose and that citrate synthase activity becomes a new bottleneck in the engineered pathway during the simultaneous utilization of d-glucose and l-arabinose.

Conclusions

Metabolome analysis identified potential bottlenecks in d-glucose and l-arabinose metabolism and was then applied to the following rational metabolic engineering. Manipulation of only two genes enabled simultaneous utilization of d-glucose and l-arabinose at the same rate in metabolically engineered C. glutamicum. This is the first report of rational metabolic design and engineering for simultaneous hexose and pentose utilization without inactivating the phosphotransferase system.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0927-6) contains supplementary material, which is available to authorized users.

Promises and pitfalls of untargeted metabolomics

Metabolomics is one of the newer omics fields, and has enabled researchers to complement genomic and protein level analysis of disease with both semi-quantitative and quantitative metabolite levels, which are the chemical mediators that constitute a given phenotype. Over more than a decade, methodologies have advanced for both targeted (quantification of specific analytes) as well as untargeted metabolomics (biomarker discovery and global metabolite profiling). Untargeted metabolomics is especially useful when there is no a priori metabolic hypothesis. Liquid chromatography coupled to mass spectrometry (LC-MS) has been the preferred choice for untargeted metabolomics, given the versatility in metabolite coverage and sensitivity of these instruments. Resolving and profiling many hundreds to thousands of metabolites with varying chemical properties in a biological sample presents unique challenges, or pitfalls. In this review, we address the various obstacles and corrective measures available in four major aspects associated with an untargeted metabolomics experiment: (1) experimental design, (2) pre-analytical (sample collection and preparation), (3) analytical (chromatography and detection), and (4) post-analytical (data processing).

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.

Bio-production of gaseous alkenes: ethylene, isoprene, isobutene

To reduce emissions from petrochemical refinement, bio-production has been heralded as a way to create economically valuable compounds with fewer harmful effects. For example, gaseous alkenes are precursor molecules that can be polymerized into a variety of industrially significant compounds and have biological production pathways. Production levels, however, remain low, thus enhancing bio-production of gaseous petrochemicals for chemical precursors is critical. This review covers the metabolic pathways and production levels of the gaseous alkenes ethylene, isoprene, and isobutene. Techniques needed to drive production to higher levels are also discussed.

Antibiotic-Induced Changes to the Host Metabolic Environment Inhibit Drug Efficacy and Alter Immune Function

Bactericidal antibiotics alter microbial metabolism as part of their lethality and can damage mitochondria in mammalian cells. In addition, antibiotic susceptibility is sensitive to extracellular metabolites, but it remains unknown whether metabolites present at an infection site can affect either treatment efficacy or immune function. Here, we quantify local metabolic changes in the host microenvironment following antibiotic treatment for a peritoneal Escherichia coli infection. Antibiotic treatment elicits microbiome-independent changes in local metabolites, but not those distal to the infection site, by acting directly on host cells. The metabolites induced during treatment, such as AMP, reduce antibiotic efficacy and enhance phagocytic killing. Moreover, antibiotic treatment impairs immune function by inhibiting respiratory activity in immune cells. Collectively, these results highlight the immunomodulatory potential of antibiotics and reveal the local metabolic microenvironment to be an important determinant of infection resolution.

Targeted redox and energy cofactor metabolomics in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum

BACKGROUND

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are prominent candidate biocatalysts that, together, can enable the direct biotic conversion of lignocellulosic biomass to ethanol. The imbalance and suboptimal turnover rates of redox cofactors are currently hindering engineering efforts to achieve higher bioproductivity in both organisms. Measuring relevant intracellular cofactor concentrations will help understand redox state of these cofactors and help identify a strategy to overcome these limitations; however, metabolomic determinations of these labile metabolites have historically proved challenging.

RESULTS

Through our validations, we verified the handling and storage stability of these metabolites, and verified extraction matrices and extraction solvent were not suppressing mass spectrometry signals. We recovered adenylate energy charge ratios (a main quality indicator) above 0.82 for all extractions. NADH/NAD+ values of 0.26 and 0.04 for an adhE-deficient strain of C. thermocellum and its parent, respectively, reflect the expected shift to a more reduced redox potential when a species lacks the ability to re-oxidize NADH by synthesizing ethanol. This method failed to yield reliable results with C. bescii and poor-growing strains of T. saccharolyticum.

CONCLUSIONS

Our validated protocols demonstrate and validate the extraction and analysis of selected redox and energy-related metabolites from two candidate consolidated bioprocessing biocatalysts, C. thermocellum and T. saccharolyticum. This development and validation highlights the important, but often neglected, need to optimize and validate metabolomic protocols when adapting them to new cell or tissue types.

Engineering cofactor flexibility enhanced 2,3-butanediol production in Escherichia coli

Enzymatic reduction of acetoin into 2,3-butanediol (2,3-BD) typically requires the reduced nicotinamide adenine dinucleotide (NADH) or its phosphate form (NADPH) as electron donor. Efficiency of 2,3-BD biosynthesis, therefore, is heavily influenced by the enzyme specificity and the cofactor availability which varies dynamically. This work describes the engineering of cofactor flexibility for 2,3-BD production by simultaneous overexpression of an NADH-dependent 2,3-BD dehydrogenase from Klebsiella pneumoniae (KpBudC) and an NADPH-specific 2,3-BD dehydrogenase from Clostridium beijerinckii (CbAdh). Co-expression of KpBudC and CbAdh not only enabled condition versatility for 2,3-BD synthesis via flexible utilization of cofactors, but also improved production stereo-specificity of 2,3-BD without accumulation of acetoin. With optimization of medium and fermentation condition, the co-expression strain produced 92 g/L of 2,3-BD in 56 h with 90% stereo-purity for (R,R)-isoform and 85% of maximum theoretical yield. Incorporating cofactor flexibility into the design principle should benefit production of bio-based chemical involving redox reactions.

Reexamination of the Physiological Role of PykA in Escherichia coli Revealed that It Negatively Regulates the Intracellular ATP Levels under Anaerobic Conditions

Pyruvate kinase is one of the three rate-limiting glycolytic enzymes that catalyze the last step of glycolysis, conversion of phosphoenolpyruvate (PEP) into pyruvate, which is associated with ATP generation. Two isozymes of pyruvate kinase, PykF and PykA, are identified in Escherichia coli PykF is considered important, whereas PykA has a less-defined role. Prior studies inactivated the pykA gene to increase the level of its substrate, PEP, and thereby increased the yield of end products derived from PEP. We were surprised when we found a pykA::Tn5 mutant in a screen for increased yield of an end product derived from pyruvate (n-butanol), suggesting that the role of PykA needs to be reexamined. We show that the pykA mutant exhibited elevated intracellular ATP levels, biomass concentrations, glucose consumption, and n-butanol production. We also discovered that the pykA mutant expresses higher levels of a presumed pyruvate transporter, YhjX, permitting the mutant to recapture and metabolize excreted pyruvate. Furthermore, we demonstrated that the nucleotide diphosphate kinase activity of PykA leads to negative regulation of the intracellular ATP levels. Taking the data together, we propose that inactivation of pykA can be considered a general strategy to enhance the production of pyruvate-derived metabolites under anaerobic conditions.IMPORTANCE This study showed that knocking out pykA significantly increased the intracellular ATP level and thus significantly increased the levels of glucose consumption, biomass formation, and pyruvate-derived product formation under anaerobic conditions. pykA was considered to be encoding a dispensable pyruvate kinase; here we show that pykA negatively regulates the anaerobic glycolysis rate through regulating the energy distribution. Thus, knocking out pykA can be used as a general strategy to increase the level of pyruvate-derived fermentative products.

Elucidating dynamic metabolic physiology through network integration of quantitative time-course metabolomics

The increasing availability of metabolomics data necessitates novel methods for deeper data analysis and interpretation. We present a flux balance analysis method that allows for the computation of dynamic intracellular metabolic changes at the cellular scale through integration of time-course absolute quantitative metabolomics. This approach, termed "unsteady-state flux balance analysis" (uFBA), is applied to four cellular systems: three dynamic and one steady-state as a negative control. uFBA and FBA predictions are contrasted, and uFBA is found to be more accurate in predicting dynamic metabolic flux states for red blood cells, platelets, and Saccharomyces cerevisiae. Notably, only uFBA predicts that stored red blood cells metabolize TCA intermediates to regenerate important cofactors, such as ATP, NADH, and NADPH. These pathway usage predictions were subsequently validated through ¹³C isotopic labeling and metabolic flux analysis in stored red blood cells. Utilizing time-course metabolomics data, uFBA provides an accurate method to predict metabolic physiology at the cellular scale for dynamic systems.

Quantification of Microbial Phenotypes

Metabolite profiling technologies have improved to generate close to quantitative metabolomics data, which can be employed to quantitatively describe the metabolic phenotype of an organism. Here, we review the current technologies available for quantitative metabolomics, present their advantages and drawbacks, and the current challenges to generate fully quantitative metabolomics data. Metabolomics data can be integrated into metabolic networks using thermodynamic principles to constrain the directionality of reactions. Here we explain how to estimate Gibbs energy under physiological conditions, including examples of the estimations, and the different methods for thermodynamics-based network analysis. The fundamentals of the methods and how to perform the analyses are described. Finally, an example applying quantitative metabolomics to a yeast model by 13C fluxomics and thermodynamics-based network analysis is presented. The example shows that (1) these two methods are complementary to each other; and (2) there is a need to take into account Gibbs energy errors. Better estimations of metabolic phenotypes will be obtained when further constraints are included in the analysis.

Thermodynamics-based Metabolite Sensitivity Analysis in metabolic networks

The increasing availability of large metabolomics datasets enhances the need for computational methodologies that can organize the data in a way that can lead to the inference of meaningful relationships. Knowledge of the metabolic state of a cell and how it responds to various stimuli and extracellular conditions can offer significant insight in the regulatory functions and how to manipulate them. Constraint based methods, such as Flux Balance Analysis (FBA) and Thermodynamics-based flux analysis (TFA), are commonly used to estimate the flow of metabolites through genome-wide metabolic networks, making it possible to identify the ranges of flux values that are consistent with the studied physiological and thermodynamic conditions. However, unless key intracellular fluxes and metabolite concentrations are known, constraint-based models lead to underdetermined problem formulations. This lack of information propagates as uncertainty in the estimation of fluxes and basic reaction properties such as the determination of reaction directionalities. Therefore, knowledge of which metabolites, if measured, would contribute the most to reducing this uncertainty can significantly improve our ability to define the internal state of the cell. In the present work we combine constraint based modeling, Design of Experiments (DoE) and Global Sensitivity Analysis (GSA) into the Thermodynamics-based Metabolite Sensitivity Analysis (TMSA) method. TMSA ranks metabolites comprising a metabolic network based on their ability to constrain the gamut of possible solutions to a limited, thermodynamically consistent set of internal states. TMSA is modular and can be applied to a single reaction, a metabolic pathway or an entire metabolic network. This is, to our knowledge, the first attempt to use metabolic modeling in order to provide a significance ranking of metabolites to guide experimental measurements.

Homeostasis of metabolites in Escherichia coli on transition from anaerobic to aerobic conditions and the transient secretion of pyruvate

We have developed a method for rapid quenching of samples taken from chemostat cultures of Escherichia coli that gives reproducible and reliable measurements of extracellular and intracellular metabolites by ¹H NMR and have applied it to study the major central metabolites during the transition from anaerobic to aerobic growth. Almost all metabolites showed a gradual change after perturbation with air, consistent with immediate inhibition of pyruvate formate-lyase, dilution of overflow metabolites and induction of aerobic enzymes. Surprisingly, although pyruvate showed almost no change in intracellular concentration, the extracellular concentration transiently increased. The absence of intracellular accumulation of pyruvate suggested that one or more glycolytic enzymes might relocate to the cell membrane. To test this hypothesis, chromosomal pyruvate kinase (pykF) was modified to express either PykF-green fluorescent protein or PykF-FLAG fusion proteins. Measurements showed that PykF-FLAG relocates to the cell membrane within 5 min of aeration and then slowly returns to the cytoplasm, suggesting that on aeration, PykF associates with the membrane to facilitate secretion of pyruvate to maintain constant intracellular levels.

Networks of energetic and metabolic interactions define dynamics in microbial communities

Microorganisms form diverse communities that have a profound impact on the environment and human health. Recent technological advances have enabled elucidation of community diversity at high resolution. Investigation of microbial communities has revealed that they often contain multiple members with complementing and seemingly redundant metabolic capabilities. An understanding of the communal impacts of redundant metabolic capabilities is currently lacking; specifically, it is not known whether metabolic redundancy will foster competition or motivate cooperation. By investigating methanogenic populations, we identified the multidimensional interspecies interactions that define composition and dynamics within syntrophic communities that play a key role in the global carbon cycle. Species-specific genomes were extracted from metagenomic data using differential coverage binning. We used metabolic modeling leveraging metatranscriptomic information to reveal and quantify a complex intertwined system of syntrophic relationships. Our results show that amino acid auxotrophies create additional interdependencies that define community composition and control carbon and energy flux through the system while simultaneously contributing to overall community robustness. Strategic use of antimicrobials further reinforces this intricate interspecies network. Collectively, our study reveals the multidimensional interactions in syntrophic communities that promote high species richness and bolster community stability during environmental perturbations.

Heading in the right direction: thermodynamics-based network analysis and pathway engineering

Thermodynamics-based network analysis through the introduction of thermodynamic constraints in metabolic models allows a deeper analysis of metabolism and guides pathway engineering. The number and the areas of applications of thermodynamics-based network analysis methods have been increasing in the last ten years. We review recent applications of these methods and we identify the areas that such analysis can contribute significantly, and the needs for future developments. We find that organisms with multiple compartments and extremophiles present challenges for modeling and thermodynamics-based flux analysis. The evolution of current and new methods must also address the issues of the multiple alternatives in flux directionalities and the uncertainties and partial information from analytical methods.

Escher: A Web Application for Building, Sharing, and Embedding Data-Rich Visualizations of Biological Pathways

Escher is a web application for visualizing data on biological pathways. Three key features make Escher a uniquely effective tool for pathway visualization. First, users can rapidly design new pathway maps. Escher provides pathway suggestions based on user data and genome-scale models, so users can draw pathways in a semi-automated way. Second, users can visualize data related to genes or proteins on the associated reactions and pathways, using rules that define which enzymes catalyze each reaction. Thus, users can identify trends in common genomic data types (e.g. RNA-Seq, proteomics, ChIP)—in conjunction with metabolite- and reaction-oriented data types (e.g. metabolomics, fluxomics). Third, Escher harnesses the strengths of web technologies (SVG, D3, developer tools) so that visualizations can be rapidly adapted, extended, shared, and embedded. This paper provides examples of each of these features and explains how the development approach used for Escher can be used to guide the development of future visualization tools.

We are now in the age of big data. More than ever before, biological discoveries require powerful and flexible tools for managing large datasets, including both visual and statistical tools. Pathway-based visualization is particularly powerful since it enables one to analyze complex datasets within the context of actual biological processes and to elucidate how each change in a cell effects related processes. To facilitate such approaches, we present Escher, a web application that can be used to rapidly build pathway maps. On Escher maps, diverse datasets related to genes, reactions, and metabolites can be quickly contextualized within metabolism and, increasingly, beyond metabolism. Escher is available now for free use (under the MIT license) at https://escher.github.io.

Systems biology of host-microbe metabolomics

The human gut microbiota performs essential functions for host and well‐being, but has also been linked to a variety of disease states, e.g., obesity and type 2 diabetes. The mammalian body fluid and tissue metabolomes are greatly influenced by the microbiota, with many health‐relevant metabolites being considered ‘mammalian–microbial co‐metabolites’. To systematically investigate this complex host–microbial co‐metabolism, a systems biology approach integrating high‐throughput data and computational network models is required. Here, we review established top‐down and bottom‐up systems biology approaches that have successfully elucidated relationships between gut microbiota‐derived metabolites and host health and disease. We focus particularly on the constraint‐based modeling and analysis approach, which enables the prediction of mechanisms behind metabolic host–microbe interactions on the molecular level. We illustrate that constraint‐based models are a useful tool for the contextualization of metabolomic measurements and can further our insight into host–microbe interactions, yielding, e.g., in potential novel drugs and biomarkers. WIREs Syst Biol Med 2015, 7:195–219. doi: 10.1002/wsbm.1301

For further resources related to this article, please visit the WIREs website.

Conflict of interest: The authors have declared no conflicts of interest for this article.

Construction of Escherichia Coli Cell Factories for Production of Organic Acids and Alcohols

Production of bulk chemicals from renewable biomass has been proved to be sustainable and environmentally friendly. Escherichia coli is the most commonly used host strain for constructing cell factories for production of bulk chemicals since it has clear physiological and genetic characteristics, grows fast in minimal salts medium, uses a wide range of substrates, and can be genetically modified easily. With the development of metabolic engineering, systems biology, and synthetic biology, a technology platform has been established to construct E. coli cell factories for bulk chemicals production. In this chapter, we will introduce this technology platform, as well as E. coli cell factories successfully constructed for production of organic acids and alcohols.

Perturbation Experiments: Approaches for Metabolic Pathway Analysis in Bioreactors

In the last decades, targeted metabolic engineering of microbial cells has become one of the major tools in bioprocess design and optimization. For successful application, a detailed knowledge is necessary about the relevant metabolic pathways and their regulation inside the cells. Since in vitro experiments cannot display process conditions and behavior properly, process data about the cells' metabolic state have to be collected in vivo. For this purpose, special techniques and methods are necessary. Therefore, most techniques enabling in vivo characterization of metabolic pathways rely on perturbation experiments, which can be divided into dynamic and steady-state approaches. To avoid any process disturbance, approaches which enable perturbation of cell metabolism in parallel to the continuing production process are reasonable. Furthermore, the fast dynamics of microbial production processes amplifies the need of parallelized data generation. These points motivate the development of a parallelized approach for multiple metabolic perturbation experiments outside the operating production reactor. An appropriate approach for in vivo characterization of metabolic pathways is presented and applied exemplarily to a microbial L-phenylalanine production process on a 15 L-scale.