Reconstruction and modeling protein translocation and compartmentalization in Escherichia coli at the genome-scale

Precise metabolic modeling in post-omics era: accomplishments and perspectives

Microbes have been extensively utilized for their sustainable and scalable properties in synthesizing desired bio-products. However, insufficient knowledge about intracellular metabolism has impeded further microbial applications. The genome-scale metabolic models (GEMs) play a pivotal role in facilitating a global understanding of cellular metabolic mechanisms. These models enable rational modification by exploring metabolic pathways and predicting potential targets in microorganisms, enabling precise cell regulation without experimental costs. Nonetheless, simplified GEM only considers genome information and network stoichiometry while neglecting other important bio-information, such as enzyme functions, thermodynamic properties, and kinetic parameters. Consequently, uncertainties persist particularly when predicting microbial behaviors in complex and fluctuant systems. The advent of the omics era with its massive quantification of genes, proteins, and metabolites under various conditions has led to the flourishing of multi-constrained models and updated algorithms with improved predicting power and broadened dimension. Meanwhile, machine learning (ML) has demonstrated exceptional analytical and predictive capacities when applied to training sets of biological big data. Incorporating the discriminant strength of ML with GEM facilitates mechanistic modeling efficiency and improves predictive accuracy. This paper provides an overview of research innovations in the GEM, including multi-constrained modeling, analytical approaches, and the latest applications of ML, which may contribute comprehensive knowledge toward genetic refinement, strain development, and yield enhancement for a broad range of biomolecules.

Establishing comprehensive quaternary structural proteomes from genome sequence

A critical body of knowledge has developed through advances in protein microscopy, protein-fold modeling, structural biology software, availability of sequenced bacterial genomes, large-scale mutation databases, and genome-scale models. Based on these recent advances, we develop a computational framework that; i) identifies the oligomeric structural proteome encoded by an organism's genome from available structural resources; ii) maps multi-strain alleleomic variation, resulting in the structural proteome for a species; and iii) calculates the 3D orientation of proteins across subcellular compartments with residue-level precision. Using the platform, we; iv) compute the quaternary E. coli K-12 MG1655 structural proteome; v) use a dataset of 12,000 mutations to build Random Forest classifiers that can predict the severity of mutations; and, in combination with a genome-scale model that computes proteome allocation, vi) obtain the spatial allocation of the E. coli proteome. Thus, in conjunction with relevant datasets and increasingly accurate computational models, we can now annotate quaternary structural proteomes, at genome-scale, to obtain a molecular-level understanding of whole-cell functions.

Significance

Advancements in experimental and computational methods have revealed the shapes of multi-subunit proteins. The absence of a unified platform that maps actionable datatypes onto these increasingly accurate structures creates a barrier to structural analyses, especially at the genome-scale. Here, we describe QSPACE, a computational annotation platform that evaluates existing resources to identify the best-available structure for each protein in a user's query, maps the 3D location of actionable datatypes ( e.g. , active sites, published mutations) onto the selected structures, and uses third-party APIs to determine the subcellular compartment of all amino acids of a protein. As proof-of-concept, we deployed QSPACE to generate the quaternary structural proteome of E. coli MG1655 and demonstrate two use-cases involving large-scale mutant analysis and genome-scale modelling.

Cardinality optimization in constraint-based modelling: application to human metabolism

MOTIVATION

Several applications in constraint-based modelling can be mathematically formulated as cardinality optimization problems involving the minimization or maximization of the number of nonzeros in a vector. These problems include testing for stoichiometric consistency, testing for flux consistency, testing for thermodynamic flux consistency, computing sparse solutions to flux balance analysis problems and computing the minimum number of constraints to relax to render an infeasible flux balance analysis problem feasible. Such cardinality optimization problems are computationally complex, with no known polynomial time algorithms capable of returning an exact and globally optimal solution.

RESULTS

By approximating the zero-norm with nonconvex continuous functions, we reformulate a set of cardinality optimization problems in constraint-based modelling into a difference of convex functions. We implemented and numerically tested novel algorithms that approximately solve the reformulated problems using a sequence of convex programs. We applied these algorithms to various biochemical networks and demonstrate that our algorithms match or outperform existing related approaches. In particular, we illustrate the efficiency and practical utility of our algorithms for cardinality optimization problems that arise when extracting a model ready for thermodynamic flux balance analysis given a human metabolic reconstruction.

AVAILABILITY AND IMPLEMENTATION

Open source scripts to reproduce the results are here https://github.com/opencobra/COBRA.papers/2023_cardOpt with general purpose functions integrated within the COnstraint-Based Reconstruction and Analysis toolbox: https://github.com/opencobra/cobratoolbox.

Construction and Analysis of an Enzyme-Constrained Metabolic Model of Corynebacterium glutamicum

The genome-scale metabolic model (GEM) is a powerful tool for interpreting and predicting cellular phenotypes under various environmental and genetic perturbations. However, GEM only considers stoichiometric constraints, and the simulated growth and product yield values will show a monotonic linear increase with increasing substrate uptake rate, which deviates from the experimentally measured values. Recently, the integration of enzymatic constraints into stoichiometry-based GEMs was proven to be effective in making novel discoveries and predicting new engineering targets. Here, we present the first genome-scale enzyme-constrained model (ecCGL1) for Corynebacterium glutamicum reconstructed by integrating enzyme kinetic data from various sources using a ECMpy workflow based on the high-quality GEM of C. glutamicum (obtained by modifying the iCW773 model). The enzyme-constrained model improved the prediction of phenotypes and simulated overflow metabolism, while also recapitulating the trade-off between biomass yield and enzyme usage efficiency. Finally, we used the ecCGL1 to identify several gene modification targets for l-lysine production, most of which agree with previously reported genes. This study shows that incorporating enzyme kinetic information into the GEM enhances the cellular phenotypes prediction of C. glutamicum, which can help identify key enzymes and thus provide reliable guidance for metabolic engineering.

Characterizing the Interplay of Rubisco and Nitrogenase Enzymes in Anaerobic-Photoheterotrophically Grown Rhodopseudomonas palustris CGA009 through a Genome-Scale Metabolic and Expression Model

Rhodopseudomonas palustris CGA009 is a Gram-negative purple nonsulfur bacterium that grows phototrophically by fixing carbon dioxide and nitrogen or chemotrophically by fixing or catabolizing a wide array of substrates, including lignin breakdown products for its carbon and fixing nitrogen for its nitrogen requirements. It can grow aerobically or anaerobically and can use light, inorganic, and organic compounds for energy production. Due to its ability to convert different carbon sources into useful products during anaerobic growth, this study reconstructed a metabolic and expression (ME) model of R. palustris to investigate its anaerobic-photoheterotrophic growth. Unlike metabolic (M) models, ME models include transcription and translation reactions along with macromolecules synthesis and couple these reactions with growth rate. This unique feature of the ME model led to nonlinear growth curve predictions, which matched closely with experimental growth rate data. At the theoretical maximum growth rate, the ME model suggested a diminishing rate of carbon fixation and predicted malate dehydrogenase and glycerol-3 phosphate dehydrogenase as alternate electron sinks. Moreover, the ME model also identified ferredoxin as a key regulator in distributing electrons between major redox balancing pathways. Because ME models include the turnover rate for each metabolic reaction, it was used to successfully capture experimentally observed temperature regulation of different nitrogenases. Overall, these unique features of the ME model demonstrated the influence of nitrogenases and rubiscos on R. palustris growth and predicted a key regulator in distributing electrons between major redox balancing pathways, thus establishing a platform for in silico investigation of R. palustris metabolism from a multiomics perspective. IMPORTANCE In this work, we reconstructed the first ME model for a purple nonsulfur bacterium (PNSB). Using the ME model, different aspects of R. palustris metabolism were examined. First, the ME model was used to analyze how reducing power entering the R. palustris cell through organic carbon sources gets partitioned into biomass, carbon dioxide fixation, and nitrogen fixation. Furthermore, the ME model predicted electron flux through ferredoxin as a major bottleneck in distributing electrons to nitrogenase enzymes. Next, the ME model characterized different nitrogenase enzymes and successfully recapitulated experimentally observed temperature regulations of those enzymes. Identifying the bottleneck responsible for transferring an electron to nitrogenase enzymes and recapitulating the temperature regulation of different nitrogenase enzymes can have profound implications in metabolic engineering, such as hydrogen production from R. palustris. Another interesting application of this ME model can be to take advantage of its redox balancing strategy to gain an understanding of the regulatory mechanism of biodegradable plastic production precursors, such as polyhydroxybutyrate (PHB).

Construction of Multiscale Genome-Scale Metabolic Models: Frameworks and Challenges

Genome-scale metabolic models (GEMs) are effective tools for metabolic engineering and have been widely used to guide cell metabolic regulation. However, the single gene–protein-reaction data type in GEMs limits the understanding of biological complexity. As a result, multiscale models that add constraints or integrate omics data based on GEMs have been developed to more accurately predict phenotype from genotype. This review summarized the recent advances in the development of multiscale GEMs, including multiconstraint, multiomic, and whole-cell models, and outlined machine learning applications in GEM construction. This review focused on the frameworks, toolkits, and algorithms for constructing multiscale GEMs. The challenges and perspectives of multiscale GEM development are also discussed.

Genome-Scale Metabolic Modeling Enables In-Depth Understanding of Big Data

Genome-scale metabolic models (GEMs) enable the mathematical simulation of the metabolism of archaea, bacteria, and eukaryotic organisms. GEMs quantitatively define a relationship between genotype and phenotype by contextualizing different types of Big Data (e.g., genomics, metabolomics, and transcriptomics). In this review, we analyze the available Big Data useful for metabolic modeling and compile the available GEM reconstruction tools that integrate Big Data. We also discuss recent applications in industry and research that include predicting phenotypes, elucidating metabolic pathways, producing industry-relevant chemicals, identifying drug targets, and generating knowledge to better understand host-associated diseases. In addition to the up-to-date review of GEMs currently available, we assessed a plethora of tools for developing new GEMs that include macromolecular expression and dynamic resolution. Finally, we provide a perspective in emerging areas, such as annotation, data managing, and machine learning, in which GEMs will play a key role in the further utilization of Big Data.

Constraint-based modeling of yeast mitochondria reveals the dynamics of protein import and iron-sulfur cluster biogenesis

Mitochondria are a hallmark of eukaryal cells and play an important role in cellular metabolism. There is a vast amount of knowledge available on mitochondrial metabolism and essential mitochondrial functions, such as protein import and iron-sulfur cluster biosynthesis, including multiple studies on the mitochondrial proteome. Therefore, there is a need for in silico approaches to facilitate the analysis of these data. Here, we present a detailed model of mitochondrial metabolism Saccharomyces cerevisiae, including protein import, iron-sulfur cluster biosynthesis, and a description of the coupling between charge translocation processes and ATP synthesis. Model analysis implied a dual dependence of absolute levels of proteins in protein import, iron-sulfur cluster biogenesis and cluster abundance on growth rate and respiratory activity. The model is instrumental in studying dynamics and perturbations in these processes and given the high conservation of mitochondrial metabolism in humans, it can provide insight into their role in human disease.

•

Reconstruction of mitochondrial protein import and cofactor metabolism in yeast

•

Quantification of the energy cost of metabolite transport

•

Protein import activity depends on growth rate and respiratory activity

•

Quantification iron-sulfur cluster requirements show growth rate dependence

Genome-Scale Reconstruction of Microbial Dynamic Phenotype: Successes and Challenges

This review is a part of the SI 'Genome-Scale Modeling of Microorganisms in the Real World'. The goal of GEM is the accurate prediction of the phenotype from its respective genotype under specified environmental conditions. This review focuses on the dynamic phenotype; prediction of the real-life behaviors of microorganisms, such as cell proliferation, dormancy, and mortality; balanced and unbalanced growth; steady-state and transient processes; primary and secondary metabolism; stress responses; etc. Constraint-based metabolic reconstructions were successfully started two decades ago as FBA, followed by more advanced models, but this review starts from the earlier nongenomic predecessors to show that some GEMs inherited the outdated biokinetic frameworks compromising their performances. The most essential deficiencies are: (i) an inadequate account of environmental conditions, such as various degrees of nutrients limitation and other factors shaping phenotypes; (ii) a failure to simulate the adaptive changes of MMCC (MacroMolecular Cell Composition) in response to the fluctuating environment; (iii) the misinterpretation of the SGR (Specific Growth Rate) as either a fixed constant parameter of the model or independent factor affecting the conditional expression of macromolecules; (iv) neglecting stress resistance as an important objective function; and (v) inefficient experimental verification of GEM against simple growth (constant MMCC and SGR) data. Finally, we propose several ways to improve GEMs, such as replacing the outdated Monod equation with the SCM (Synthetic Chemostat Model) that establishes the quantitative relationships between primary and secondary metabolism, growth rate and stress resistance, process kinetics, and cell composition.

Combining Kinetic and Constraint-Based Modelling to Better Understand Metabolism Dynamics

To understand the phenotypic capabilities of organisms, it is useful to characterise cellular metabolism through the analysis of its pathways. Dynamic mathematical modelling of metabolic networks is of high interest as it provides the time evolution of the metabolic components. However, it also has limitations, such as the necessary mechanistic details and kinetic parameters are not always available. On the other hand, large metabolic networks exhibit a complex topological structure which can be studied rather efficiently in their stationary regime by constraint-based methods. These methods produce useful predictions on pathway operations. In this review, we present both modelling techniques and we show how they bring complementary views of metabolism. In particular, we show on a simple example how both approaches can be used in conjunction to shed some light on the dynamics of metabolic networks.

Kinetic, metabolic, and statistical analytics: addressing metabolic transport limitations among organelles and microbial communities

Microbial organisms engage in a variety of metabolic interactions. A crucial part of these interactions is the exchange of molecules between different organelles, cells, and the environment. The main forces mediating this metabolic exchange are transporters. This transport can be difficult to measure experimentally because several transport mechanisms remain opaque. However, theoretical calculations about the inputs and outputs of cells via metabolic exchanges have enabled the successful inference of the workings of intra-organismal and inter-organismal systems. Kinetic, metabolic, and statistical modeling approaches in combination with omics data are enhancing our knowledge and understanding about metabolic exchange and mass resource allocation. This model-driven analytics approach can guide effective experimental design and yield new insights into biological function and control.

Genome-scale metabolic modeling reveals key features of a minimal gene set

Mesoplasma florum, a fast-growing near-minimal organism, is a compelling model to explore rational genome designs. Using sequence and structural homology, the set of metabolic functions its genome encodes was identified, allowing the reconstruction of a metabolic network representing ˜ 30% of its protein-coding genes. Growth medium simplification enabled substrate uptake and product secretion rate quantification which, along with experimental biomass composition, were integrated as species-specific constraints to produce the functional iJL208 genome-scale model (GEM) of metabolism. Genome-wide expression and essentiality datasets as well as growth data on various carbohydrates were used to validate and refine iJL208. Discrepancies between model predictions and observations were mechanistically explained using protein structures and network analysis. iJL208 was also used to propose an in silico reduced genome. Comparing this prediction to the minimal cell JCVI-syn3.0 and its parent JCVI-syn1.0 revealed key features of a minimal gene set. iJL208 is a stepping-stone toward model-driven whole-genome engineering.

Integrating thermodynamic and enzymatic constraints into genome-scale metabolic models

Stoichiometric genome-scale metabolic network models (GEMs) have been widely used to predict metabolic phenotypes. In addition to stoichiometric ratios, other constraints such as enzyme availability and thermodynamic feasibility can also limit the phenotype solution space. Extended GEM models considering either enzymatic or thermodynamic constraints have been shown to improve prediction accuracy. In this paper, we propose a novel method that integrates both enzymatic and thermodynamic constraints in a single Pyomo modeling framework (ETGEMs). We applied this method to construct the EcoETM (E. coli metabolic model with enzymatic and thermodynamic constraints). Using this model, we calculated the optimal pathways for cellular growth and the production of 22 metabolites. When comparing the results with those of iML1515 and models with one of the two constraints, we observed that many thermodynamically unfavorable and/or high enzyme cost pathways were excluded from EcoETM. For example, the synthesis pathway of carbamoyl-phosphate (Cbp) from iML1515 is both thermodynamically unfavorable and enzymatically costly. After introducing the new constraints, the production pathways and yields of several Cbp-derived products (e.g. L-arginine, orotate) calculated using EcoETM were more realistic. The results of this study demonstrate the great application potential of metabolic models with multiple constraints for pathway analysis and phenotype prediction.

Computation of condition-dependent proteome allocation reveals variability in the macro and micro nutrient requirements for growth

Sustaining a robust metabolic network requires a balanced and fully functioning proteome. In addition to amino acids, many enzymes require cofactors (coenzymes and engrafted prosthetic groups) to function properly. Extensively validated resource allocation models, such as genome-scale models of metabolism and gene expression (ME-models), have the ability to compute an optimal proteome composition underlying a metabolic phenotype, including the provision of all required cofactors. Here we apply the ME-model for Escherichia coli K-12 MG1655 to computationally examine how environmental conditions change the proteome and its accompanying cofactor usage. We found that: (1) The cofactor requirements computed by the ME-model mostly agree with the standard biomass objective function used in models of metabolism alone (M-models); (2) ME-model computations reveal non-intuitive variability in cofactor use under different growth conditions; (3) An analysis of ME-model predicted protein use in aerobic and anaerobic conditions suggests an enrichment in the use of peroxyl scavenging acids in the proteins used to sustain aerobic growth; (4) The ME-model could describe how limitation in key protein components affect the metabolic state of E. coli. Genome-scale models have thus reached a level of sophistication where they reveal intricate properties of functional proteomes and how they support different E. coli lifestyles.

Escherichia coli is capable of growing in many environments, each of which requires a different collection of enzymes to metabolize the nutrients within that environment. Each individual enzyme requires its own set of amino acids and oftentimes cofactors, which are accessory molecules essential for the enzyme to function. Thus, the composition of the micronutrients (amino acids, cofactors, etc.) within a cell will differ depending on its metabolic needs. The presented work is the first effort to employ metabolic models to probe the connection between E. coli’s diverse growth environments and its biomass composition. We first show how differences in model-predicted enzyme use for aerobic or anaerobic growth results in distinct amino acid and cofactor usage. Alternatively, we show that the metabolic models can predict how modifying the cell’s biomass composition will affect growth. For example, by modeling the exposure of E. coli to trimethoprim or sulfamethoxazole—two antibiotics that target folate (vitamin B9) synthesis—we predicted how E. coli could adapt to grow under folate-limited conditions. This work demonstrates how models can be used to study antibiotic resistance of drugs that target amino acid or cofactor synthesis.

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

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.

Leveraging genome-scale metabolic models for human health applications

Genome-scale metabolic modeling is a scalable and extensible computational method for analyzing and predicting biological function. With the ongoing improvements in computational methods and experimental capabilities, genome-scale metabolic models (GEMs) are demonstrating utility in addressing human health applications. The initial areas of highest impact are likely to be health applications where disease states involve metabolic changes. In this review, we focus on recent application of GEMs to studying cancer and the human microbiome by describing the enabling methodologies and outcomes of these studies. We conclude with proposing some areas of research that are likely to arise as a result of recent methodological advances.

Synthesizing Systems Biology Knowledge from Omics Using Genome-Scale Models

Omic technologies have enabled the complete readout of the molecular state of a cell at different biological scales. In principle, the combination of multiple omic data types can provide an integrated view of the entire biological system. This integration requires appropriate models in a systems biology approach. Here, genome-scale models (GEMs) are focused upon as one computational systems biology approach for interpreting and integrating multi-omic data. GEMs convert the reactions (related to metabolism, transcription, and translation) that occur in an organism to a mathematical formulation that can be modeled using optimization principles. A variety of genome-scale modeling methods used to interpret multiple omic data types, including genomics, transcriptomics, proteomics, metabolomics, and meta-omics are reviewed. The ability to interpret omics in the context of biological systems has yielded important findings for human health, environmental biotechnology, bioenergy, and metabolic engineering. The authors find that concurrent with advancements in omic technologies, genome-scale modeling methods are also expanding to enable better interpretation of omic data. Therefore, continued synthesis of valuable knowledge, through the integration of omic data with GEMs, are expected.

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.

Modeling of nitrogen fixation and polymer production in the heterotrophic diazotroph Azotobacter vinelandii DJ

Nitrogen fixation is an important metabolic process carried out by microorganisms, which converts molecular nitrogen into inorganic nitrogenous compounds such as ammonia (NH₃). These nitrogenous compounds are crucial for biogeochemical cycles and for the synthesis of essential biomolecules, i.e. nucleic acids, amino acids and proteins. Azotobacter vinelandii is a bacterial non-photosynthetic model organism to study aerobic nitrogen fixation (diazotrophy) and hydrogen production. Moreover, the diazotroph can produce biopolymers like alginate and polyhydroxybutyrate (PHB) that have important industrial applications. However, many metabolic processes such as partitioning of carbon and nitrogen metabolism in A. vinelandii remain unknown to date.

Genome-scale metabolic models (M-models) represent reliable tools to unravel and optimize metabolic functions at genome-scale. M-models are mathematical representations that contain information about genes, reactions, metabolites and their associations. M-models can simulate optimal reaction fluxes under a wide variety of conditions using experimentally determined constraints. Here we report on the development of a M-model of the wild type bacterium A. vinelandii DJ (iDT1278) which consists of 2,003 metabolites, 2,469 reactions, and 1,278 genes. We validated the model using high-throughput phenotypic and physiological data, testing 180 carbon sources and 95 nitrogen sources. iDT1278 was able to achieve an accuracy of 89% and 91% for growth with carbon sources and nitrogen source, respectively. This comprehensive M-model will help to comprehend metabolic processes associated with nitrogen fixation, ammonium assimilation, and production of organic nitrogen in an environmentally important microorganism.

•

Genome-scale metabolic model of Azotobacter vinelandii DJ achives over 90% accuracy.

•

iDT1278 is the most comprehensive model to simulate diazotrophy.

•

Determining the most suitable culture conditions to produce polymers A. vinelandii.

•

Constraint-based modeling unravels links among nitrogen fixation and production of organic nitrogen.

Adaptations of Escherichia coli strains to oxidative stress are reflected in properties of their structural proteomes

BACKGROUND

The reconstruction of metabolic networks and the three-dimensional coverage of protein structures have reached the genome-scale in the widely studied Escherichia coli K-12 MG1655 strain. The combination of the two leads to the formation of a structural systems biology framework, which we have used to analyze differences between the reactive oxygen species (ROS) sensitivity of the proteomes of sequenced strains of E. coli. As proteins are one of the main targets of oxidative damage, understanding how the genetic changes of different strains of a species relates to its oxidative environment can reveal hypotheses as to why these variations arise and suggest directions of future experimental work.

RESULTS

Creating a reference structural proteome for E. coli allows us to comprehensively map genetic changes in 1764 different strains to their locations on 4118 3D protein structures. We use metabolic modeling to predict basal ROS production levels (ROStype) for 695 of these strains, finding that strains with both higher and lower basal levels tend to enrich their proteomes with antioxidative properties, and speculate as to why that is. We computationally assess a strain's sensitivity to an oxidative environment, based on known chemical mechanisms of oxidative damage to protein groups, defined by their localization and functionality. Two general groups - metalloproteins and periplasmic proteins - show enrichment of their antioxidative properties between the 695 strains with a predicted ROStype as well as 116 strains with an assigned pathotype. Specifically, proteins that a) utilize a molybdenum ion as a cofactor and b) are involved in the biogenesis of fimbriae show intriguing protective properties to resist oxidative damage. Overall, these findings indicate that a strain's sensitivity to oxidative damage can be elucidated from the structural proteome, though future experimental work is needed to validate our model assumptions and findings.

CONCLUSION

We thus demonstrate that structural systems biology enables a proteome-wide, computational assessment of changes to atomic-level physicochemical properties and of oxidative damage mechanisms for multiple strains in a species. This integrative approach opens new avenues to study adaptation to a particular environment based on physiological properties predicted from sequence alone.

Genome-scale reconstructions of the mammalian secretory pathway predict metabolic costs and limitations of protein secretion

In mammalian cells, >25% of synthesized proteins are exported through the secretory pathway. The pathway complexity, however, obfuscates its impact on the secretion of different proteins. Unraveling its impact on diverse proteins is particularly important for biopharmaceutical production. Here we delineate the core secretory pathway functions and integrate them with genome-scale metabolic reconstructions of human, mouse, and Chinese hamster ovary cells. The resulting reconstructions enable the computation of energetic costs and machinery demands of each secreted protein. By integrating additional omics data, we find that highly secretory cells have adapted to reduce expression and secretion of other expensive host cell proteins. Furthermore, we predict metabolic costs and maximum productivities of biotherapeutic proteins and identify protein features that most significantly impact protein secretion. Finally, the model successfully predicts the increase in secretion of a monoclonal antibody after silencing a highly expressed selection marker. This work represents a knowledgebase of the mammalian secretory pathway that serves as a novel tool for systems biotechnology.

Genome-scale model of metabolism and gene expression provides a multi-scale description of acid stress responses in Escherichia coli

Response to acid stress is critical for Escherichia coli to successfully complete its life-cycle by passing through the stomach to colonize the digestive tract. To develop a fundamental understanding of this response, we established a molecular mechanistic description of acid stress mitigation responses in E. coli and integrated them with a genome-scale model of its metabolism and macromolecular expression (ME-model). We considered three known mechanisms of acid stress mitigation: 1) change in membrane lipid fatty acid composition, 2) change in periplasmic protein stability over external pH and periplasmic chaperone protection mechanisms, and 3) change in the activities of membrane proteins. After integrating these mechanisms into an established ME-model, we could simulate their responses in the context of other cellular processes. We validated these simulations using RNA sequencing data obtained from five E. coli strains grown under external pH ranging from 5.5 to 7.0. We found: i) that for the differentially expressed genes accounted for in the ME-model, 80% of the upregulated genes were correctly predicted by the ME-model, and ii) that these genes are mainly involved in translation processes (45% of genes), membrane proteins and related processes (18% of genes), amino acid metabolism (12% of genes), and cofactor and prosthetic group biosynthesis (8% of genes). We also demonstrated several intervention strategies on acid tolerance that can be simulated by the ME-model. We thus established a quantitative framework that describes, on a genome-scale, the acid stress mitigation response of E. coli that has both scientific and practical uses.

The y-ome defines the 35% of Escherichia coli genes that lack experimental evidence of function

Experimental studies of Escherichia coli K-12 MG1655 often implicate poorly annotated genes in cellular phenotypes. However, we lack a systematic understanding of these genes. How many are there? What information is available for them? And what features do they share that could explain the gap in our understanding? Efforts to build predictive, whole-cell models of E. coli inevitably face this knowledge gap. We approached these questions systematically by assembling annotations from the knowledge bases EcoCyc, EcoGene, UniProt and RegulonDB. We identified the genes that lack experimental evidence of function (the ‘y-ome’) which include 1600 of 4623 unique genes (34.6%), of which 111 have absolutely no evidence of function. An additional 220 genes (4.7%) are pseudogenes or phantom genes. y-ome genes tend to have lower expression levels and are enriched in the termination region of the E. coli chromosome. Where evidence is available for y-ome genes, it most often points to them being membrane proteins and transporters. We resolve the misconception that a gene in E. coli whose primary name starts with ‘y’ is unannotated, and we discuss the value of the y-ome for systematic improvement of E. coli knowledge bases and its extension to other organisms.

Involvement of AmpG in mediating a dynamic relationship between serine beta-lactamase induction and biofilm-forming ability of Escherichia coli

AmpG permease is implicated both in beta-lactamase induction and peptidoglycan recycling in enterobacterial isolates. Here, physiological studies using molecular genetics show that deletion of AmpG permease dramatically increases beta-lactam susceptibility even in the presence of AmpC, TEM-1 and OXA beta-lactamases. Also, there is an appreciable decrease in the biofilm-forming ability of strains lacking this protein. Expression of this permease in excess probably compromises the integrity of the bacterial cells, leading to cell lysis. Based on these results, we propose that AmpG permease may be used as a potential antibiotic target and its suppression could efficiently inhibit both beta-lactamase induction and biofilm formation.

A system-wide network reconstruction of gene regulation and metabolism in Escherichia coli

Genome-scale metabolic models have become a fundamental tool for examining metabolic principles. However, metabolism is not solely characterized by the underlying biochemical reactions and catalyzing enzymes, but also affected by regulatory events. Since the pioneering work of Covert and co-workers as well as Shlomi and co-workers it is debated, how regulation and metabolism synergistically characterize a coherent cellular state. The first approaches started from metabolic models, which were extended by the regulation of the encoding genes of the catalyzing enzymes. By now, bioinformatics databases in principle allow addressing the challenge of integrating regulation and metabolism on a system-wide level. Collecting information from several databases we provide a network representation of the integrated gene regulatory and metabolic system for Escherichia coli, including major cellular processes, from metabolic processes via protein modification to a variety of regulatory events. Besides transcriptional regulation, we also take into account regulation of translation, enzyme activities and reactions. Our network model provides novel topological characterizations of system components based on their positions in the network. We show that network characteristics suggest a representation of the integrated system as three network domains (regulatory, metabolic and interface networks) instead of two. This new three-domain representation reveals the structural centrality of components with known high functional relevance. This integrated network can serve as a platform for understanding coherent cellular states as active subnetworks and to elucidate crossover effects between metabolism and gene regulation.

Predicting proteome allocation, overflow metabolism, and metal requirements in a model acetogen

The unique capability of acetogens to ferment a broad range of substrates renders them ideal candidates for the biotechnological production of commodity chemicals. In particular the ability to grow with H2:CO2 or syngas (a mixture of H2/CO/CO2) makes these microorganisms ideal chassis for sustainable bioproduction. However, advanced design strategies for acetogens are currently hampered by incomplete knowledge about their physiology and our inability to accurately predict phenotypes. Here we describe the reconstruction of a novel genome-scale model of metabolism and macromolecular synthesis (ME-model) to gain new insights into the biology of the model acetogen Clostridium ljungdahlii. The model represents the first ME-model of a Gram-positive bacterium and captures all major central metabolic, amino acid, nucleotide, lipid, major cofactors, and vitamin synthesis pathways as well as pathways to synthesis RNA and protein molecules necessary to catalyze these reactions, thus significantly broadens the scope and predictability. Use of the model revealed how protein allocation and media composition influence metabolic pathways and energy conservation in acetogens and accurately predicted secretion of multiple fermentation products. Predicting overflow metabolism is of particular interest since it enables new design strategies, e.g. the formation of glycerol, a novel product for C. ljungdahlii, thus broadening the metabolic capability for this model microbe. Furthermore, prediction and experimental validation of changing secretion rates based on different metal availability opens the window into fermentation optimization and provides new knowledge about the proteome utilization and carbon flux in acetogens.

DynamicME: dynamic simulation and refinement of integrated models of metabolism and protein expression

BACKGROUND

Genome-scale models of metabolism and macromolecular expression (ME models) enable systems-level computation of proteome allocation coupled to metabolic phenotype.

RESULTS

We develop DynamicME, an algorithm enabling time-course simulation of cell metabolism and protein expression. DynamicME correctly predicted the substrate utilization hierarchy on a mixed carbon substrate medium. We also found good agreement between predicted and measured time-course expression profiles. ME models involve considerably more parameters than metabolic models (M models). We thus generate an ensemble of models (each model having its rate constants perturbed), and then analyze the models by identifying archetypal time-course metabolite concentration profiles. Furthermore, we use a metaheuristic optimization method to calibrate ME model parameters using time-course measurements such as from a (fed-) batch culture. Finally, we show that constraints on protein concentration dynamics ("inertia") alter the metabolic response to environmental fluctuations, including increased substrate-level phosphorylation and lowered oxidative phosphorylation.

CONCLUSIONS

Overall, DynamicME provides a novel method for understanding proteome allocation and metabolism under complex and transient environments, and to utilize time-course cell culture data for model-based interpretation or model refinement.

Identification of growth-coupled production strains considering protein costs and kinetic variability

Conversion of renewable biomass to useful molecules in microbial cell factories can be approached in a rational and systematic manner using constraint-based reconstruction and analysis. Filtering for high confidence in silico designs is critical because in vivo construction and testing of strains is expensive and time consuming. As such, a workflow was devised to analyze the robustness of growth-coupled production when considering the biosynthetic costs of the proteome and variability in enzyme kinetic parameters using a genome-scale model of metabolism and gene expression (ME-model). A collection of 2632 unfiltered knockout designs in Escherichia coli was evaluated by the workflow. A ME-model was used in the workflow to test the designs' growth-coupled production in addition to a less complex genome-scale metabolic model (M-model). The workflow identified 634 M-model growth-coupled designs which met the filtering criteria and 42 robust designs, which met growth-coupled production criteria using both M and ME-models. Knockouts were found to follow a pattern of controlling intermediate metabolite consumption such as pyruvate consumption and high flux subsystems such as glycolysis. Kinetic parameter sampling using the ME-model revealed how enzyme efficiency and pathway tradeoffs can affect growth-coupled production phenotypes.

Modeling the multi-scale mechanisms of macromolecular resource allocation

As microbes face changing environments, they dynamically allocate macromolecular resources to produce a particular phenotypic state. Broad 'omics' data sets have revealed several interesting phenomena regarding how the proteome is allocated under differing conditions, but the functional consequences of these states and how they are achieved remain open questions. Various types of multi-scale mathematical models have been used to elucidate the genetic basis for systems-level adaptations. In this review, we outline several different strategies by which microbes accomplish resource allocation and detail how mathematical models have aided in our understanding of these processes. Ultimately, such modeling efforts have helped elucidate the principles of proteome allocation and hold promise for further discovery.

Distribution of mechanical stress in the Escherichia coli cell envelope

The cell envelope in Gram-negative bacteria comprises two distinct membranes with a cell wall between them. There has been a growing interest in understanding the mechanical adaptation of this cell envelope to the osmotic pressure (or turgor pressure), which is generated by the difference in the concentration of solutes between the cytoplasm and the external environment. However, it remains unexplored how the cell wall, the inner membrane (IM), and the outer membrane (OM) effectively protect the cell from this pressure by bearing the resulting surface tension, thus preventing the formation of inner membrane bulges, abnormal cell morphology, spheroplasts and cell lysis. In this study, we have used molecular dynamics (MD) simulations combined with experiments to resolve how and to what extent models of the IM, OM, and cell wall respond to changes in surface tension. We calculated the area compressibility modulus of all three components in simulations from tension-area isotherms. Experiments on monolayers mimicking individual leaflets of the IM and OM were also used to characterize their compressibility. While the membranes become softer as they expand, the cell wall exhibits significant strain stiffening at moderate to high tensions. We integrate these results into a model of the cell envelope in which the OM and cell wall share the tension at low turgor pressure (0.3 atm) but the tension in the cell wall dominates at high values (>1 atm).

In Silico Prediction of Large-Scale Microbial Production Performance: Constraints for Getting Proper Data-Driven Models

Industrial bioreactors range from 10.000 to 700.000 L and characteristically show different zones of substrate availabilities, dissolved gas concentrations and pH values reflecting physical, technical and economic constraints of scale-up. Microbial producers are fluctuating inside the bioreactors thereby experiencing frequently changing micro-environmental conditions. The external stimuli induce responses on microbial metabolism and on transcriptional regulation programs. Both may deteriorate the expected microbial production performance in large scale compared to expectations deduced from ideal, well-mixed lab-scale conditions. Accordingly, predictive tools are needed to quantify large-scale impacts considering bioreactor heterogeneities. The review shows that the time is right to combine simulations of microbial kinetics with calculations of large-scale environmental conditions to predict the bioreactor performance. Accordingly, basic experimental procedures and computational tools are presented to derive proper microbial models and hydrodynamic conditions, and to link both for bioreactor modeling. Particular emphasis is laid on the identification of gene regulatory networks as the implementation of such models will surely gain momentum in future studies.

COBRAme: A computational framework for genome-scale models of metabolism and gene expression

Genome-scale models of metabolism and macromolecular expression (ME-models) explicitly compute the optimal proteome composition of a growing cell. ME-models expand upon the well-established genome-scale models of metabolism (M-models), and they enable a new fundamental understanding of cellular growth. ME-models have increased predictive capabilities and accuracy due to their inclusion of the biosynthetic costs for the machinery of life, but they come with a significant increase in model size and complexity. This challenge results in models which are both difficult to compute and challenging to understand conceptually. As a result, ME-models exist for only two organisms (Escherichia coli and Thermotoga maritima) and are still used by relatively few researchers. To address these challenges, we have developed a new software framework called COBRAme for building and simulating ME-models. It is coded in Python and built on COBRApy, a popular platform for using M-models. COBRAme streamlines computation and analysis of ME-models. It provides tools to simplify constructing and editing ME-models to enable ME-model reconstructions for new organisms. We used COBRAme to reconstruct a condensed E. coli ME-model called iJL1678b-ME. This reformulated model gives functionally identical solutions to previous E. coli ME-models while using 1/6 the number of free variables and solving in less than 10 minutes, a marked improvement over the 6 hour solve time of previous ME-model formulations. Errors in previous ME-models were also corrected leading to 52 additional genes that must be expressed in iJL1678b-ME to grow aerobically in glucose minimal in silico media. This manuscript outlines the architecture of COBRAme and demonstrates how ME-models can be created, modified, and shared most efficiently using the new software framework.

Modelling the metabolism of protein secretion through the Tat route in Streptomyces lividans

Background

Streptomyces lividans has demonstrated its value as an efficient host for protein production due to its ability to secrete functional proteins directly to the media. Secretory proteins that use the major Sec route need to be properly folded outside the cell, whereas secretory proteins using the Tat route appear outside the cell correctly folded. This feature makes the Tat system very attractive for the production of natural or engineered Tat secretory proteins. S. lividans cells are known to respond differently to overproduction and secretion of Tat versus Sec proteins. Increased understanding of the impact of protein secretion through the Tat route can be obtained by a deeper analysis of the metabolic impact associated with protein production, and its dependence on protein origin, composition, secretion mechanisms, growth phases and nutrients. Flux Balance Analysis of Genome-Scale Metabolic Network models provides a theoretical framework to investigate cell metabolism under different constraints.

Results

We have built new models for various S. lividans strains to better understand the mechanisms associated with overproduction of proteins secreted through the Tat route. We compare models of an S. lividans Tat-dependent agarase overproducing strain with those of the S. lividans wild-type, an S. lividans strain carrying the multi-copy plasmid vector and an α-amylase Sec-dependent overproducing strain. Using updated genomic, transcriptomic and experimental data we could extend existing S. lividans models and produce a new model which produces improved results largely extending the coverage of S. lividans strains, the number of genes and reactions being considered, the predictive behaviour and the dependence on specification of exchange constraints. Comparison of the optimized solutions obtained highlights numerous changes between Tat- and Sec-dependent protein secreting strains affecting the metabolism of carbon, amino acids, nucleotides, lipids and cofactors, and variability analysis predicts a large potential for protein overproduction.

Conclusions

This work provides a detailed look to metabolic changes associated to Tat-dependent protein secretion reproducing experimental observations and identifying changes that are specific to each secretory route, presenting a novel, improved, more accurate and strain-independent model of S. lividans, thus opening the way for enhanced metabolic engineering of protein overproduction in S. lividans.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1199-3) contains supplementary material, which is available to authorized users.

Predicting Dynamic Metabolic Demands in the Photosynthetic Eukaryote Chlorella vulgaris

Phototrophic organisms exhibit a highly dynamic proteome, adapting their biomass composition in response to diurnal light/dark cycles and nutrient availability. Here, we used experimentally determined biomass compositions over the course of growth to determine and constrain the biomass objective function (BOF) in a genome-scale metabolic model of Chlorella vulgaris UTEX 395 over time. Changes in the BOF, which encompasses all metabolites necessary to produce biomass, influence the state of the metabolic network thus directly affecting predictions. Simulations using dynamic BOFs predicted distinct proteome demands during heterotrophic or photoautotrophic growth. Model-driven analysis of extracellular nitrogen concentrations and predicted nitrogen uptake rates revealed an intracellular nitrogen pool, which contains 38% of the total nitrogen provided in the medium for photoautotrophic and 13% for heterotrophic growth. Agreement between flux and gene expression trends was determined by statistical comparison. Accordance between predicted flux trends and gene expression trends was found for 65% of multisubunit enzymes and 75% of allosteric reactions. Reactions with the highest agreement between simulations and experimental data were associated with energy metabolism, terpenoid biosynthesis, fatty acids, nucleotides, and amino acid metabolism. Furthermore, predicted flux distributions at each time point were compared with gene expression data to gain new insights into intracellular compartmentalization, specifically for transporters. A total of 103 genes related to internal transport reactions were identified and added to the updated model of C. vulgaris, iCZ946, thus increasing our knowledgebase by 10% for this model green alga.

Elucidation of complexity and prediction of interactions in microbial communities

Microorganisms engage in complex interactions with other members of the microbial community, higher organisms as well as their environment. However, determining the exact nature of these interactions can be challenging due to the large number of members in these communities and the manifold of interactions they can engage in. Various omic data, such as 16S rRNA gene sequencing, shotgun metagenomics, metatranscriptomics, metaproteomics and metabolomics, have been deployed to unravel the community structure, interactions and resulting community dynamics in situ. Interpretation of these multi-omic data often requires advanced computational methods. Modelling approaches are powerful tools to integrate, contextualize and interpret experimental data, thus shedding light on the underlying processes shaping the microbiome. Here, we review current methods and approaches, both experimental and computational, to elucidate interactions in microbial communities and to predict their responses to perturbations.

Metabolism of the fast-growing bacterium Vibrio natriegens elucidated by 13C metabolic flux analysis

Vibrio natriegens is a fast-growing, non-pathogenic bacterium that is being considered as the next-generation workhorse for the biotechnology industry. However, little is known about the metabolism of this organism which is limiting our ability to apply rational metabolic engineering strategies. To address this critical gap in current knowledge, here we have performed a comprehensive analysis of V. natriegens metabolism. We constructed a detailed model of V. natriegens core metabolism, measured the biomass composition, and performed high-resolution ¹³C metabolic flux analysis (¹³C-MFA) to estimate intracellular fluxes using parallel labeling experiments with the optimal tracers [1,2-¹³C]glucose and [1,6-¹³C]glucose. During exponential growth in glucose minimal medium, V. natriegens had a growth rate of 1.70 1/h (doubling time of 24min) and a glucose uptake rate of 3.90g/g/h, which is more than two 2-fold faster than E. coli, although slower than the fast-growing thermophile Geobacillus LC300. ¹³C-MFA revealed that the core metabolism of V. natriegens is similar to that of E. coli, with the main difference being a 33% lower normalized flux through the oxidative pentose phosphate pathway. Quantitative analysis of co-factor balances provided additional insights into the energy and redox metabolism of V. natriegens. Taken together, the results presented in this study provide valuable new information about the physiology of V. natriegens and establish a solid foundation for future metabolic engineering efforts with this promising microorganism.

Fast growth phenotype of E. coli K-12 from adaptive laboratory evolution does not require intracellular flux rewiring

Adaptive laboratory evolution (ALE) is a widely-used method for improving the fitness of microorganisms in selected environmental conditions. It has been applied previously to Escherichia coli K-12 MG1655 during aerobic exponential growth on glucose minimal media, a frequently used model organism and growth condition, to probe the limits of E. coli growth rate and gain insights into fast growth phenotypes. Previous studies have described up to 1.6-fold increases in growth rate following ALE, and have identified key causal genetic mutations and changes in transcriptional patterns. Here, we report for the first time intracellular metabolic fluxes for six such adaptively evolved strains, as determined by high-resolution ¹³C-metabolic flux analysis. Interestingly, we found that intracellular metabolic pathway usage changed very little following adaptive evolution. Instead, at the level of central carbon metabolism the faster growth was facilitated by proportional increases in glucose uptake and all intracellular rates. Of the six evolved strains studied here, only one strain showed a small degree of flux rewiring, and this was also the strain with unique genetic mutations. A comparison of fluxes with two other wild-type (unevolved) E. coli strains, BW25113 and BL21, showed that inter-strain differences are greater than differences between the parental and evolved strains. Principal component analysis highlighted that nearly all flux differences (95%) between the nine strains were captured by only two principal components. The distance between measured and flux balance analysis predicted fluxes was also investigated. It suggested a relatively wide range of similar stoichiometric optima, which opens new questions about the path-dependency of adaptive evolution.

Improving the phenotype predictions of a yeast genome-scale metabolic model by incorporating enzymatic constraints

Genome-scale metabolic models (GEMs) are widely used to calculate metabolic phenotypes. They rely on defining a set of constraints, the most common of which is that the production of metabolites and/or growth are limited by the carbon source uptake rate. However, enzyme abundances and kinetics, which act as limitations on metabolic fluxes, are not taken into account. Here, we present GECKO, a method that enhances a GEM to account for enzymes as part of reactions, thereby ensuring that each metabolic flux does not exceed its maximum capacity, equal to the product of the enzyme's abundance and turnover number. We applied GECKO to a Saccharomyces cerevisiae GEM and demonstrated that the new model could correctly describe phenotypes that the previous model could not, particularly under high enzymatic pressure conditions, such as yeast growing on different carbon sources in excess, coping with stress, or overexpressing a specific pathway. GECKO also allows to directly integrate quantitative proteomics data; by doing so, we significantly reduced flux variability of the model, in over 60% of metabolic reactions. Additionally, the model gives insight into the distribution of enzyme usage between and within metabolic pathways. The developed method and model are expected to increase the use of model-based design in metabolic engineering.

Progress in Molecular Dynamics Simulations of Gram-Negative Bacterial Cell Envelopes

Bacteria are protected by complex molecular architectures known as the cell envelope. The cell envelope is composed of regions with distinct chemical compositions and physical properties, namely, membranes and a cell wall. To develop novel antibiotics to combat pathogenic bacteria, molecular level knowledge of the structure, dynamics, and interplay between the chemical components of the cell envelope that surrounds bacterial cells is imperative. In addition, conserved molecular patterns associated with the bacterial envelope are recognized by receptors as part of the mammalian defensive response to infection, and an improved understanding of bacteria-host interactions would facilitate the search for novel immunotherapeutics. This Perspective introduces an emerging area of computational biology: multiscale molecular dynamics simulations of chemically complex models of bacterial lipids and membranes. We discuss progress to date, and identify areas for future development that will enable the study of aspects of the membrane components that are as yet unexplored by computational methods.

Application of theoretical methods to increase succinate production in engineered strains

Computational methods have enabled the discovery of non-intuitive strategies to enhance the production of a variety of target molecules. In the case of succinate production, reviews covering the topic have not yet analyzed the impact and future potential that such methods may have. In this work, we review the application of computational methods to the production of succinic acid. We found that while a total of 26 theoretical studies were published between 2002 and 2016, only 10 studies reported the successful experimental implementation of any kind of theoretical knowledge. None of the experimental studies reported an exact application of the computational predictions. However, the combination of computational analysis with complementary strategies, such as directed evolution and comparative genome analysis, serves as a proof of concept and demonstrates that successful metabolic engineering can be guided by rational computational methods.