Modeling the Contribution of Allosteric Regulation for Flux Control in the Central Carbon Metabolism of E. coli

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

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

Omics data for sampling thermodynamically feasible kinetic models

Kinetic models are key to understanding and predicting the dynamic behaviour of metabolic systems. Traditional models require kinetic parameters which are not always available and are often estimated in vitro. Ensemble models overcome this challenge by sampling thermodynamically feasible models around a measured reference point. However, it is unclear if the convenient distributions used to generate the ensemble produce a natural distribution of model parameters and hence if the model predictions are reasonable. In this paper, we produced a detailed kinetic model for the central carbon metabolism of Escherichia coli. The model consists of 82 reactions (including 13 reactions with allosteric regulation) and 79 metabolites. To sample the model, we used metabolomic and fluxomic data from a single steady-state time point for E. coli K-12 MG1655 growing on glucose minimal M9 medium (average sampling time for 1000 models: 11.21 ± 0.14 min). Afterwards, in order to examine whether our sampled models are biologically sound, we calculated Km, Vmax and kcat for the reactions and compared them to previously published values. Finally, we used metabolic control analysis to identify enzymes with high control over the fluxes in the central carbon metabolism. Our analyses demonstrate that our platform samples thermodynamically feasible kinetic models, which are in agreement with previously published experimental results and can be used to investigate metabolic control patterns within cells. This renders it a valuable tool for the study of cellular metabolism and the design of metabolic pathways.

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

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

Thermodynamic limitations of PHB production from formate and fructose in Cupriavidus necator

The chemolithotroph Cupriavidus necator H16 is known as a natural producer of the bioplastic-polymer PHB, as well as for its metabolic versatility to utilize different substrates, including formate as the sole carbon and energy source. Depending on the entry point of the substrate, this versatility requires adjustment of the thermodynamic landscape to maintain sufficiently high driving forces for biological processes. Here we employed a model of the core metabolism of C. necator H16 to analyze the thermodynamic driving forces and PHB yields from formate for different metabolic engineering strategies. For this, we enumerated elementary flux modes (EFMs) of the network and evaluated their PHB yields as well as thermodynamics via Max-min driving force (MDF) analysis and random sampling of driving forces. A heterologous ATP:citrate lyase reaction was predicted to increase driving force for producing acetyl-CoA. A heterologous phosphoketolase reaction was predicted to increase maximal PHB yields as well as driving forces. These enzymes were then verified experimentally to enhance PHB titers between 60 and 300% in select conditions. The EFM analysis also revealed that PHB production from formate may be limited by low driving forces through citrate lyase and aconitase, as well as cofactor balancing, and identified additional reactions associated with low and high PHB yield. Proteomics analysis of the engineered strains confirmed an increased abundance of aconitase and cofactor balancing. The findings of this study aid in understanding metabolic adaptation. Furthermore, the outlined approach will be useful in designing metabolic engineering strategies in other non-model bacteria.

Metabolic pathway design for growth-associated phenylalanine production using synthetically designed mutualism

Combination of growth-associated pathway engineering based on flux balance analysis (FBA) and adaptive laboratory evolution (ALE) is a powerful approach to enhance the production of useful compounds. However, the feasibility of such growth-associated pathway designs depends on the type of target compound. In the present study, FBA predicted a set of gene deletions (pykA, pykF, ppc, zwf, and adhE) that leads to growth-associated phenylalanine production in Escherichia coli. The knockout strain is theoretically enforced to produce phenylalanine only at high growth yields, and could not be applied to the ALE experiment because of a severe growth defect. To overcome this challenge, we propose a novel approach for ALE based on mutualistic co-culture for coupling growth and production, regardless of the growth rate. We designed a synthetic mutualism of a phenylalanine-producing leucine-auxotrophic strain (KF strain) and a leucine-producing phenylalanine-auxotrophic strain (KL strain) and performed an ALE experiment for approximately 160 generations. The evolved KF strain (KF-E strain) grew in a synthetic medium (with glucose as the main carbon source) supplemented with leucine, while severe growth defects were observed in the parental KF strain. The phenylalanine yield of the KF-E strain was 2.3 times higher than that of the KF strain.

Next-Generation Genome-Scale Metabolic Modeling through Integration of Regulatory Mechanisms

Genome-scale metabolic models (GEMs) are powerful tools for understanding metabolism from a systems-level perspective. However, GEMs in their most basic form fail to account for cellular regulation. A diverse set of mechanisms regulate cellular metabolism, enabling organisms to respond to a wide range of conditions. This limitation of GEMs has prompted the development of new methods to integrate regulatory mechanisms, thereby enhancing the predictive capabilities and broadening the scope of GEMs. Here, we cover integrative models encompassing six types of regulatory mechanisms: transcriptional regulatory networks (TRNs), post-translational modifications (PTMs), epigenetics, protein-protein interactions and protein stability (PPIs/PS), allostery, and signaling networks. We discuss 22 integrative GEM modeling methods and how these have been used to simulate metabolic regulation during normal and pathological conditions. While these advances have been remarkable, there remains a need for comprehensive and widespread integration of regulatory constraints into GEMs. We conclude by discussing challenges in constructing GEMs with regulation and highlight areas that need to be addressed for the successful modeling of metabolic regulation. Next-generation integrative GEMs that incorporate multiple regulatory mechanisms and their crosstalk will be invaluable for discovering cell-type and disease-specific metabolic control mechanisms.

SCOUR: a stepwise machine learning framework for predicting metabolite-dependent regulatory interactions

BACKGROUND

The topology of metabolic networks is both well-studied and remarkably well-conserved across many species. The regulation of these networks, however, is much more poorly characterized, though it is known to be divergent across organisms-two characteristics that make it difficult to model metabolic networks accurately. While many computational methods have been built to unravel transcriptional regulation, there have been few approaches developed for systems-scale analysis and study of metabolic regulation. Here, we present a stepwise machine learning framework that applies established algorithms to identify regulatory interactions in metabolic systems based on metabolic data: stepwise classification of unknown regulation, or SCOUR.

RESULTS

We evaluated our framework on both noiseless and noisy data, using several models of varying sizes and topologies to show that our approach is generalizable. We found that, when testing on data under the most realistic conditions (low sampling frequency and high noise), SCOUR could identify reaction fluxes controlled only by the concentration of a single metabolite (its primary substrate) with high accuracy. The positive predictive value (PPV) for identifying reactions controlled by the concentration of two metabolites ranged from 32 to 88% for noiseless data, 9.2 to 49% for either low sampling frequency/low noise or high sampling frequency/high noise data, and 6.6-27% for low sampling frequency/high noise data, with results typically sufficiently high for lab validation to be a practical endeavor. While the PPVs for reactions controlled by three metabolites were lower, they were still in most cases significantly better than random classification.

CONCLUSIONS

SCOUR uses a novel approach to synthetically generate the training data needed to identify regulators of reaction fluxes in a given metabolic system, enabling metabolomics and fluxomics data to be leveraged for regulatory structure inference. By identifying and triaging the most likely candidate regulatory interactions, SCOUR can drastically reduce the amount of time needed to identify and experimentally validate metabolic regulatory interactions. As high-throughput experimental methods for testing these interactions are further developed, SCOUR will provide critical impact in the development of predictive metabolic models in new organisms and pathways.

Global mapping of protein-metabolite interactions in Saccharomyces cerevisiae reveals that Ser-Leu dipeptide regulates phosphoglycerate kinase activity

Protein-metabolite interactions are of crucial importance for all cellular processes but remain understudied. Here, we applied a biochemical approach named PROMIS, to address the complexity of the protein-small molecule interactome in the model yeast Saccharomyces cerevisiae. By doing so, we provide a unique dataset, which can be queried for interactions between 74 small molecules and 3982 proteins using a user-friendly interface available at https://promis.mpimp-golm.mpg.de/yeastpmi/ . By interpolating PROMIS with the list of predicted protein-metabolite interactions, we provided experimental validation for 225 binding events. Remarkably, of the 74 small molecules co-eluting with proteins, 36 were proteogenic dipeptides. Targeted analysis of a representative dipeptide, Ser-Leu, revealed numerous protein interactors comprising chaperones, proteasomal subunits, and metabolic enzymes. We could further demonstrate that Ser-Leu binding increases activity of a glycolytic enzyme phosphoglycerate kinase (Pgk1). Consistent with the binding analysis, Ser-Leu supplementation leads to the acute metabolic changes and delays timing of a diauxic shift. Supported by the dipeptide accumulation analysis our work attests to the role of Ser-Leu as a metabolic regulator at the interface of protein degradation and central metabolism.

A review of methods for the reconstruction and analysis of integrated genome-scale models of metabolism and regulation

The current survey aims to describe the main methodologies for extending the reconstruction and analysis of genome-scale metabolic models and phenotype simulation with Flux Balance Analysis mathematical frameworks, via the integration of Transcriptional Regulatory Networks and/or gene expression data. Although the surveyed methods are aimed at improving phenotype simulations obtained from these models, the perspective of reconstructing integrated genome-scale models of metabolism and gene expression for diverse prokaryotes is still an open challenge.

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.

OptRAM: In-silico strain design via integrative regulatory-metabolic network modeling

The ultimate goal of metabolic engineering is to produce desired compounds on an industrial scale in a cost effective manner. To address challenges in metabolic engineering, computational strain optimization algorithms based on genome-scale metabolic models have increasingly been used to aid in overproducing products of interest. However, most of these strain optimization algorithms utilize a metabolic network alone, with few approaches providing strategies that also include transcriptional regulation. Moreover previous integrated approaches generally require a pre-existing regulatory network. In this study, we developed a novel strain design algorithm, named OptRAM (Optimization of Regulatory And Metabolic Networks), which can identify combinatorial optimization strategies including overexpression, knockdown or knockout of both metabolic genes and transcription factors. OptRAM is based on our previous IDREAM integrated network framework, which makes it able to deduce a regulatory network from data. OptRAM uses simulated annealing with a novel objective function, which can ensure a favorable coupling between desired chemical and cell growth. The other advance we propose is a systematic evaluation metric of multiple solutions, by considering the essential genes, flux variation, and engineering manipulation cost. We applied OptRAM to generate strain designs for succinate, 2,3-butanediol, and ethanol overproduction in yeast, which predicted high minimum predicted target production rate compared with other methods and previous literature values. Moreover, most of the genes and TFs proposed to be altered by OptRAM in these scenarios have been validated by modification of the exact genes or the target genes regulated by the TFs, for overproduction of these desired compounds by in vivo experiments cataloged in the LASER database. Particularly, we successfully validated the predicted strain optimization strategy for ethanol production by fermentation experiment. In conclusion, OptRAM can provide a useful approach that leverages an integrated transcriptional regulatory network and metabolic network to guide metabolic engineering applications.

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

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

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.

Small molecule signaling, regulation, and potential applications in cellular therapeutics

Small molecules have many important roles across the tree of life: they regulate processes from metabolism to transcription, they enable signaling within and between species, and they serve as the biochemical building blocks for cells. They also represent valuable phenotypic endpoints that are promising for use as biomarkers of disease states. In the context of engineering cell-based therapeutics, they hold particularly great promise for enabling finer control over the therapeutic cells and allowing them to be responsive to extracellular cues. The natural signaling and regulatory functions of small molecules can be harnessed and rewired to control cell activity and delivery of therapeutic payloads, potentially increasing efficacy while decreasing toxicity. To that end, this review considers small molecule-mediated regulation and signaling in bacteria. We first discuss some of the most prominent applications and aspirations for responsive cell-based therapeutics. We then describe the transport, signaling, and regulation associated with three classes of molecules that may be exploited in the engineering of therapeutic bacteria: amino acids, fatty acids, and quorum-sensing signaling molecules. We also present examples of existing engineering efforts to generate cells that sense and respond to levels of different small molecules. Finally, we discuss future directions for how small molecule-mediated regulation could be harnessed for therapeutic applications, as well as some critical considerations for the ultimate success of such endeavors. WIREs Syst Biol Med 2018, 10:e1405. doi: 10.1002/wsbm.1405 This article is categorized under: Biological Mechanisms > Cell Signaling Biological Mechanisms > Metabolism Translational, Genomic, and Systems Medicine > Therapeutic Methods.

Modeling and simulation of the redox regulation of the metabolism in Escherichia coli at different oxygen concentrations

BACKGROUND

Microbial production of biofuels and biochemicals from renewable feedstocks has received considerable recent attention from environmental protection and energy production perspectives. Many biofuels and biochemicals are produced by fermentation under oxygen-limited conditions following initiation of aerobic cultivation to enhance the cell growth rate. Thus, it is of significant interest to investigate the effect of dissolved oxygen concentration on redox regulation in Escherichia coli, a particularly popular cellular factory due to its high growth rate and well-characterized physiology. For this, the systems biology approach such as modeling is powerful for the analysis of the metabolism and for the design of microbial cellular factories.

RESULTS

Here, we developed a kinetic model that describes the dynamics of fermentation by taking into account transcription factors such as ArcA/B and Fnr, respiratory chain reactions and fermentative pathways, and catabolite regulation. The hallmark of the kinetic model is its ability to predict the dynamics of metabolism at different dissolved oxygen levels and facilitate the rational design of cultivation methods. The kinetic model was verified based on the experimental data for a wild-type E. coli strain. The model reasonably predicted the metabolic characteristics and molecular mechanisms of fnr and arcA gene-knockout mutants. Moreover, an aerobic-microaerobic dual-phase cultivation method for lactate production in a pfl-knockout mutant exhibited promising yield and productivity.

CONCLUSIONS

It is quite important to understand metabolic regulation mechanisms from both scientific and engineering points of view. In particular, redox regulation in response to oxygen limitation is critically important in the practical production of biofuel and biochemical compounds. The developed model can thus be used as a platform for designing microbial factories to produce a variety of biofuels and biochemicals.

The self-inhibitory nature of metabolic networks and its alleviation through compartmentalization

Metabolites can inhibit the enzymes that generate them. To explore the general nature of metabolic self-inhibition, we surveyed enzymological data accrued from a century of experimentation and generated a genome-scale enzyme-inhibition network. Enzyme inhibition is often driven by essential metabolites, affects the majority of biochemical processes, and is executed by a structured network whose topological organization is reflecting chemical similarities that exist between metabolites. Most inhibitory interactions are competitive, emerge in the close neighbourhood of the inhibited enzymes, and result from structural similarities between substrate and inhibitors. Structural constraints also explain one-third of allosteric inhibitors, a finding rationalized by crystallographic analysis of allosterically inhibited L-lactate dehydrogenase. Our findings suggest that the primary cause of metabolic enzyme inhibition is not the evolution of regulatory metabolite–enzyme interactions, but a finite structural diversity prevalent within the metabolome. In eukaryotes, compartmentalization minimizes inevitable enzyme inhibition and alleviates constraints that self-inhibition places on metabolism.

Metabolites act as enzyme inhibitors, but their global impact on metabolism has scarcely been considered. Here, the authors generate a human genome-wide metabolite-enzyme inhibition network, and find that inhibition occurs largely due to limited structural diversity of metabolites, leading to a global constraint on metabolism which subcellular compartmentalization minimizes.

Metabolic Complementation in Bacterial Communities: Necessary Conditions and Optimality

Bacterial communities may display metabolic complementation, in which different members of the association partially contribute to the same biosynthetic pathway. In this way, the end product of the pathway is synthesized by the community as a whole. However, the emergence and the benefits of such complementation are poorly understood. Herein, we present a simple model to analyze the metabolic interactions among bacteria, including the host in the case of endosymbiotic bacteria. The model considers two cell populations, with both cell types encoding for the same linear biosynthetic pathway. We have found that, for metabolic complementation to emerge as an optimal strategy, both product inhibition and large permeabilities are needed. In the light of these results, we then consider the patterns found in the case of tryptophan biosynthesis in the endosymbiont consortium hosted by the aphid Cinara cedri. Using in-silico computed physicochemical properties of metabolites of this and other biosynthetic pathways, we verified that the splitting point of the pathway corresponds to the most permeable intermediate.

Development of an accurate kinetic model for the central carbon metabolism of Escherichia coli

Background

A kinetic model provides insights into the dynamic response of biological systems and predicts how their complex metabolic and gene regulatory networks generate particular functions. Of many biological systems, Escherichia coli metabolic pathways have been modeled extensively at the enzymatic and genetic levels, but existing models cannot accurately reproduce experimental behaviors in a batch culture, due to the inadequate estimation of a specific cell growth rate and a large number of unmeasured parameters.

Results

In this study, we developed a detailed kinetic model for the central carbon metabolism of E. coli in a batch culture, which includes the glycolytic pathway, tricarboxylic acid cycle, pentose phosphate pathway, Entner-Doudoroff pathway, anaplerotic pathway, glyoxylate shunt, oxidative phosphorylation, phosphotransferase system (Pts), non-Pts and metabolic gene regulations by four protein transcription factors: cAMP receptor, catabolite repressor/activator, pyruvate dehydrogenase complex repressor and isocitrate lyase regulator. The kinetic parameters were estimated by a constrained optimization method on a supercomputer. The model estimated a specific growth rate based on reaction kinetics and accurately reproduced the dynamics of wild-type E. coli and multiple genetic mutants in a batch culture.

Conclusions

This model overcame the intrinsic limitations of existing kinetic models in a batch culture, predicted the effects of multilayer regulations (allosteric effectors and gene expression) on central carbon metabolism and proposed rationally designed fast-growing cells based on understandings of molecular processes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0511-x) contains supplementary material, which is available to authorized users.

The Landscape of Evolution: Reconciling Structural and Dynamic Properties of Metabolic Networks in Adaptive Diversifications

The network of the interactions among genes, proteins, and metabolites delineates a range of potential phenotypic diversifications in a lineage, and realized phenotypic changes are the result of differences in the dynamics of the expression of the elements and interactions in this deterministic network. Regulatory mechanisms, such as hormones, mediate the relationship between the structural and dynamic properties of networks by determining how and when the elements are expressed and form a functional unit or state. Changes in regulatory mechanisms lead to variable expression of functional states of a network within and among generations. Functional properties of network elements, and the magnitude and direction of evolutionary change they determine, depend on their location within a network. Here, we examine the relationship between network structure and the dynamic mechanisms that regulate flux through a metabolic network. We review the mechanisms that control metabolic flux in enzymatic reactions and examine structural properties of the network locations that are targets of flux control. We aim to establish a predictive framework to test the contributions of structural and dynamic properties of deterministic networks to evolutionary diversifications.

Metabolic modeling of clostridia: current developments and applications

Anaerobic Clostridium spp. is an important bioproduction microbial genus that can produce solvents and utilize a broad spectrum of substrates including cellulose and syngas. Genome-scale metabolic (GSM) models are increasingly being put forth for various clostridial strains to explore their respective metabolic capabilities and suitability for various bioconversions. In this study, we have selected representative GSM models for six different clostridia (Clostridium acetobutylicum, C. beijerinckii, C. butyricum, C. cellulolyticum, C. ljungdahlii and C. thermocellum) and performed a detailed model comparison contrasting their metabolic repertoire. We also discuss various applications of these GSM models to guide metabolic engineering interventions as well as assessing cellular physiology.