Genome-scale model for Clostridium acetobutylicum: Part I. Metabolic network resolution and analysis

Modeling control and transduction of electrochemical gradients in acid-stressed bacteria

Transmembrane electrochemical gradients drive solute uptake and constitute a substantial fraction of the cellular energy pool in bacteria. These gradients act not only as "homeostatic contributors," but also play a dynamic and keystone role in several bacterial functions, including sensing, stress response, and metabolism. At the system level, multiple gradients interact with ion transporters and bacterial behavior in a complex, rapid, and emergent manner; consequently, experiments alone cannot untangle their interdependencies. Electrochemical gradient modeling provides a general framework to understand these interactions and their underlying mechanisms. We quantify the generation, maintenance, and interactions of electrical, proton, and potassium potential gradients under lactic acid-stress and lactic acid fermentation. Further, we elucidate a gradient-mediated mechanism for intracellular pH sensing and stress response. We demonstrate that this gradient model can yield insights on the energetic limitations of membrane transport, and can predict bacterial behavior across changing environments.

A Holistic Approach from Systems Biology Reveals the Direct Influence of the Quorum-Sensing Phenomenon on Pseudomonas aeruginosa Metabolism to Pyoverdine Biosynthesis

Computational modeling and simulation of biological systems have become valuable tools for understanding and predicting cellular performance and phenotype generation. This work aimed to construct, model, and dynamically simulate the virulence factor pyoverdine (PVD) biosynthesis in Pseudomonas aeruginosa through a systemic approach, considering that the metabolic pathway of PVD synthesis is regulated by the quorum-sensing (QS) phenomenon. The methodology comprised three main stages: (i) Construction, modeling, and validation of the QS gene regulatory network that controls PVD synthesis in P. aeruginosa strain PAO1; (ii) construction, curating, and modeling of the metabolic network of P. aeruginosa using the flux balance analysis (FBA) approach; (iii) integration and modeling of these two networks into an integrative model using the dynamic flux balance analysis (DFBA) approximation, followed, finally, by an in vitro validation of the integrated model for PVD synthesis in P. aeruginosa as a function of QS signaling. The QS gene network, constructed using the standard System Biology Markup Language, comprised 114 chemical species and 103 reactions and was modeled as a deterministic system following the kinetic based on mass action law. This model showed that the higher the bacterial growth, the higher the extracellular concentration of QS signal molecules, thus emulating the natural behavior of P. aeruginosa PAO1. The P. aeruginosa metabolic network model was constructed based on the iMO1056 model, the P. aeruginosa PAO1 strain genomic annotation, and the metabolic pathway of PVD synthesis. The metabolic network model included the PVD synthesis, transport, exchange reactions, and the QS signal molecules. This metabolic network model was curated and then modeled under the FBA approximation, using biomass maximization as the objective function (optimization problem, a term borrowed from the engineering field). Next, chemical reactions shared by both network models were chosen to combine them into an integrative model. To this end, the fluxes of these reactions, obtained from the QS network model, were fixed in the metabolic network model as constraints of the optimization problem using the DFBA approximation. Finally, simulations of the integrative model (CCBM1146, comprising 1123 reactions and 880 metabolites) were run using the DFBA approximation to get (i) the flux profile for each reaction, (ii) the bacterial growth profile, (iii) the biomass profile, and (iv) the concentration profiles of metabolites of interest such as glucose, PVD, and QS signal molecules. The CCBM1146 model showed that the QS phenomenon directly influences the P. aeruginosa metabolism to PVD biosynthesis as a function of the change in QS signal intensity. The CCBM1146 model made it possible to characterize and explain the complex and emergent behavior generated by the interactions between the two networks, which would have been impossible to do by studying each system's individual components or scales separately. This work is the first in silico report of an integrative model comprising the QS gene regulatory network and the metabolic network of P. aeruginosa.

Transcriptomic studies of solventogenic clostridia, Clostridium acetobutylicum and Clostridium beijerinckii

Solventogenic clostridia are not a strictly defined group within the genus Clostridium but its representatives share some common features, i.e. they are anaerobic, non-pathogenic, non-toxinogenic and endospore forming bacteria. Their main metabolite is typically 1-butanol but depending on species and culture conditions, they can form other metabolites such as acetone, isopropanol, ethanol, butyric, lactic and acetic acids, and hydrogen. Although these organisms were previously used for the industrial production of solvents, they later fell into disuse, being replaced by more efficient chemical production. A return to a more biological production of solvents therefore requires a thorough understanding of clostridial metabolism. Transcriptome analysis, which reflects the involvement of individual genes in all cellular processes within a population, at any given (sampling) moment, is a valuable tool for gaining a deeper insight into clostridial life. In this review, we describe techniques to study transcription, summarize the evolution of these techniques and compare methods for data processing and visualization of solventogenic clostridia, particularly the species Clostridium acetobutylicum and Clostridium beijerinckii. Individual approaches for evaluating transcriptomic data are compared and their contributions to advancements in the field are assessed. Moreover, utilization of transcriptomic data for reconstruction of computational clostridial metabolic models is considered and particular models are described. Transcriptional changes in glucose transport, central carbon metabolism, the sporulation cycle, butanol and butyrate stress responses, the influence of lignocellulose-derived inhibitors on growth and solvent production, and other respective topics, are addressed and common trends are highlighted.

Novel Drivers of Virulence in Clostridioides difficile Identified via Context-Specific Metabolic Network Analysis

The pathogen Clostridioides difficile causes toxin-mediated diarrhea and is the leading cause of hospital-acquired infection in the United States. Due to growing antibiotic resistance and recurrent infection, targeting C. difficile metabolism presents a new approach to combat this infection. Genome-scale metabolic network reconstructions (GENREs) have been used to identify therapeutic targets and uncover properties that determine cellular behaviors. Thus, we constructed C. difficile GENREs for a hypervirulent isolate (strain [str.] R20291) and a historic strain (str. 630), validating both with in vitro and in vivo data sets. Growth simulations revealed significant correlations with measured carbon source usage (positive predictive value [PPV] ≥ 92.7%), and single-gene deletion analysis showed >89.0% accuracy. Next, we utilized each GENRE to identify metabolic drivers of both sporulation and biofilm formation. Through contextualization of each model using transcriptomes generated from in vitro and infection conditions, we discovered reliance on the pentose phosphate pathway as well as increased usage of cytidine and N-acetylneuraminate when virulence expression is reduced, which was subsequently supported experimentally. Our results highlight the ability of GENREs to identify novel metabolite signals in higher-order phenotypes like bacterial pathogenesis. IMPORTANCE Clostridioides difficile has become the leading single cause of hospital-acquired infections. Numerous studies have demonstrated the importance of specific metabolic pathways in aspects of C. difficile pathophysiology, from initial colonization to regulation of virulence factors. In the past, genome-scale metabolic network reconstruction (GENRE) analysis of bacteria has enabled systematic investigation of the genetic and metabolic properties that contribute to downstream virulence phenotypes. With this in mind, we generated and extensively curated C. difficile GENREs for both a well-studied laboratory strain (str. 630) and a more recently characterized hypervirulent isolate (str. R20291). In silico validation of both GENREs revealed high degrees of agreement with experimental gene essentiality and carbon source utilization data sets. Subsequent exploration of context-specific metabolism during both in vitro growth and infection revealed consistent patterns of metabolism which corresponded with experimentally measured increases in virulence factor expression. Our results support that differential C. difficile virulence is associated with distinct metabolic programs related to use of carbon sources and provide a platform for identification of novel therapeutic targets.

Physiological limitations and opportunities in microbial metabolic engineering

Metabolic engineering can have a pivotal role in increasing the environmental sustainability of the transportation and chemical manufacturing sectors. The field has already developed engineered microorganisms that are currently being used in industrial-scale processes. However, it is often challenging to achieve the titres, yields and productivities required for commercial viability. The efficiency of microbial chemical production is usually dependent on the physiological traits of the host organism, which may either impose limitations on engineered biosynthetic pathways or, conversely, boost their performance. In this Review, we discuss different aspects of microbial physiology that often create obstacles for metabolic engineering, and present solutions to overcome them. We also describe various instances in which natural or engineered physiological traits in host organisms have been harnessed to benefit engineered metabolic pathways for chemical production.

Modeling Growth Kinetics, Interspecies Cell Fusion, and Metabolism of a Clostridium acetobutylicum/Clostridium ljungdahlii Syntrophic Coculture

Clostridium acetobutylicum and Clostridium ljungdahlii grown in a syntrophic culture were recently shown to fuse membranes and exchange cytosolic contents, yielding hybrid cells with significant shifts in gene expression and growth phenotypes. Here, we introduce a dynamic genome-scale metabolic modeling framework to explore how cell fusion alters the growth phenotype and panel of metabolites produced by this binary community. Computational results indicate C. ljungdahlii persists in the coculture through proteome exchange during fusing events, which endow C. ljungdahlii cells with expanded substrate utilization, and access to additional reducing equivalents from C. acetobutylicum-evolved H₂ and through acquisition of C. acetobutylicum-native cofactor-reducing enzymes. Simulations predict maximum theoretical ethanol and isopropanol yields that are increased by 0.64 mmol and 0.39 mmol per mmol hexose sugar consumed, respectively, during exponential growth when cell fusion is active. This modeling effort provides a mechanistic explanation for the metabolic outcome of cellular fusion and altered homeostasis achieved in this syntrophic clostridial community.IMPORTANCE Widespread cell fusion and protein exchange between microbial organisms as observed in synthetic C. acetobutylicum/C. ljungdahlii culture is a novel observation that has not been explored in silico The mechanisms responsible for the observed cell fusion events in this culture are still unknown. In this work, we develop a modeling framework that captures the observed culture composition and metabolic phenotype, use it to offer a mechanistic explanation for how the culture achieves homeostasis, and identify C. ljungdahlii as primary beneficiary of fusion events. The implications for the events described in this study are far reaching, with potential to reshape our understanding of microbial community behavior synthetically and in nature.

Electron availability in CO2 , CO and H2 mixtures constrains flux distribution, energy management and product formation in Clostridium ljungdahlii

Acetogens such as Clostridium ljungdahlii can play a crucial role reducing the human CO2 footprint by converting industrial emissions containing CO2 , CO and H2 into valuable products such as organic acids or alcohols. The quantitative understanding of cellular metabolism is a prerequisite to exploit the bacterial endowments and to fine-tune the cells by applying metabolic engineering tools. Studying the three gas mixtures CO2  + H2 , CO and CO + CO2  + H2 (syngas) by continuously gassed batch cultivation experiments and applying flux balance analysis, we identified CO as the preferred carbon and electron source for growth and producing alcohols. However, the total yield of moles of carbon (mol-C) per electrons consumed was almost identical in all setups which underlines electron availability as the main factor influencing product formation. The Wood-Ljungdahl pathway (WLP) showed high flexibility by serving as the key NAD+ provider for CO2  + H2, whereas this function was strongly compensated by the transhydrogenase-like Nfn complex when CO was metabolized. Availability of reduced ferredoxin (Fdred ) can be considered as a key determinant of metabolic control. Oxidation of CO via carbon monoxide dehydrogenase (CODH) is the main route of Fdred formation when CO is used as substrate, whereas Fdred is mainly regenerated via the methyl branch of WLP and the Nfn complex utilizing CO2  + H2 . Consequently, doubled growth rates, highest ATP formation rates and highest amounts of reduced products (ethanol, 2,3-butanediol) were observed when CO was the sole carbon and electron source.

Towards continuous industrial bioprocessing with solventogenic and acetogenic clostridia: challenges, progress and perspectives

The sustainable production of solvents from above ground carbon is highly desired. Several clostridia naturally produce solvents and use a variety of renewable and waste-derived substrates such as lignocellulosic biomass and gas mixtures containing H2/CO2 or CO. To enable economically viable production of solvents and biofuels such as ethanol and butanol, the high productivity of continuous bioprocesses is needed. While the first industrial-scale gas fermentation facility operates continuously, the acetone-butanol-ethanol (ABE) fermentation is traditionally operated in batch mode. This review highlights the benefits of continuous bioprocessing for solvent production and underlines the progress made towards its establishment. Based on metabolic capabilities of solvent producing clostridia, we discuss recent advances in systems-level understanding and genome engineering. On the process side, we focus on innovative fermentation methods and integrated product recovery to overcome the limitations of the classical one-stage chemostat and give an overview of the current industrial bioproduction of solvents.

Development of a Genome-Scale Metabolic Model of Clostridium thermocellum and Its Applications for Integration of Multi-Omics Datasets and Computational Strain Design

Solving environmental and social challenges such as climate change requires a shift from our current non-renewable manufacturing model to a sustainable bioeconomy. To lower carbon emissions in the production of fuels and chemicals, plant biomass feedstocks can replace petroleum using microorganisms as biocatalysts. The anaerobic thermophile Clostridium thermocellum is a promising bacterium for bioconversion due to its capability to efficiently degrade lignocellulosic biomass. However, the complex metabolism of C. thermocellum is not fully understood, hindering metabolic engineering to achieve high titers, rates, and yields of targeted molecules. In this study, we developed an updated genome-scale metabolic model of C. thermocellum that accounts for recent metabolic findings, has improved prediction accuracy, and is standard-conformant to ensure easy reproducibility. We illustrated two applications of the developed model. We first formulated a multi-omics integration protocol and used it to understand redox metabolism and potential bottlenecks in biofuel (e.g., ethanol) production in C. thermocellum. Second, we used the metabolic model to design modular cells for efficient production of alcohols and esters with broad applications as flavors, fragrances, solvents, and fuels. The proposed designs not only feature intuitive push-and-pull metabolic engineering strategies, but also present novel manipulations around important central metabolic branch-points. We anticipate the developed genome-scale metabolic model will provide a useful tool for system analysis of C. thermocellum metabolism to fundamentally understand its physiology and guide metabolic engineering strategies to rapidly generate modular production strains for effective biosynthesis of biofuels and biochemicals from lignocellulosic biomass.

Pathway dissection, regulation, engineering and application: lessons learned from biobutanol production by solventogenic clostridia

The global energy crisis and limited supply of petroleum fuels have rekindled the interest in utilizing a sustainable biomass to produce biofuel. Butanol, an advanced biofuel, is a superior renewable resource as it has a high energy content and is less hygroscopic than other candidates. At present, the biobutanol route, employing acetone-butanol-ethanol (ABE) fermentation in Clostridium species, is not economically competitive due to the high cost of feedstocks, low butanol titer, and product inhibition. Based on an analysis of the physiological characteristics of solventogenic clostridia, current advances that enhance ABE fermentation from strain improvement to product separation were systematically reviewed, focusing on: (1) elucidating the metabolic pathway and regulation mechanism of butanol synthesis; (2) enhancing cellular performance and robustness through metabolic engineering, and (3) optimizing the process of ABE fermentation. Finally, perspectives on engineering and exploiting clostridia as cell factories to efficiently produce various chemicals and materials are also discussed.

Gsmodutils: a python based framework for test-driven genome scale metabolic model development

MOTIVATION

Genome scale metabolic models (GSMMs) are increasingly important for systems biology and metabolic engineering research as they are capable of simulating complex steady-state behaviour. Constraints based models of this form can include thousands of reactions and metabolites, with many crucial pathways that only become activated in specific simulation settings. However, despite their widespread use, power and the availability of tools to aid with the construction and analysis of large scale models, little methodology is suggested for their continued management. For example, when genome annotations are updated or new understanding regarding behaviour is discovered, models often need to be altered to reflect this. This is quickly becoming an issue for industrial systems and synthetic biotechnology applications, which require good quality reusable models integral to the design, build, test and learn cycle.

RESULTS

As part of an ongoing effort to improve genome scale metabolic analysis, we have developed a test-driven development methodology for the continuous integration of validation data from different sources. Contributing to the open source technology based around COBRApy, we have developed the gsmodutils modelling framework placing an emphasis on test-driven design of models through defined test cases. Crucially, different conditions are configurable allowing users to examine how different designs or curation impact a wide range of system behaviours, minimizing error between model versions.

AVAILABILITY AND IMPLEMENTATION

The software framework described within this paper is open source and freely available from http://github.com/SBRCNottingham/gsmodutils.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Kinetic modeling of Stickland reactions-coupled methanogenesis for a methanogenic culture

Studying amino acid catabolism-coupled methanogenesis is the important standpoints to decipher the metabolic behavior of a methanogenic culture. l-Glycine and l-alanine are acted as sole carbon and nitrogen sources for acidogenic bacteria. One amino acid is oxidized and another one is reduced for acetate production via pyruvate by oxidative deamination process in the Stickland reactions. Herein, we have developed a kinetic model for the Stickland reactions-coupled methanogenesis (SRCM) and simulated objectively to maximize the rate of methane production. We collected the metabolic information from enzyme kinetic parameters for amino acid catabolism of Clostridium acetobutylicum ATCC 824 and methanogenesis of Methanosarcina acetivorans C2A. The SRCM model of this study consisted of 18 reactions and 61 metabolites with enzyme kinetic parameters derived experimental data. The internal or external metabolic flux rate of this system found to control the acidogenesis and methanogenesis in a methanogenic culture. Using the SRCM model, flux distributions were calculated for each reaction and metabolite in order to maximize the methane production rate from the glycine–alanine pair. Results of this study, we demonstrated the metabolic behavior, metabolite pairing while mutually interact, and advantages of syntrophic metabolism of amino acid-directed methane production in a methanogenic starter culture.

Metabolism the Difficile Way: The Key to the Success of the Pathogen Clostridioides difficile

Strains of Clostridioides difficile cause detrimental diarrheas with thousands of deaths worldwide. The infection process by the Gram-positive, strictly anaerobic gut bacterium is directly related to its unique metabolism, using multiple Stickland-type amino acid fermentation reactions coupled to Rnf complex-mediated sodium/proton gradient formation for ATP generation. Major pathways utilize phenylalanine, leucine, glycine and proline with the formation of 3-phenylproprionate, isocaproate, butyrate, 5-methylcaproate, valerate and 5-aminovalerate. In parallel a versatile sugar catabolism including pyruvate formate-lyase as a central enzyme and an incomplete tricarboxylic acid cycle to prevent unnecessary NADH formation completes the picture. However, a complex gene regulatory network that carefully mediates the continuous adaptation of this metabolism to changing environmental conditions is only partially elucidated. It involves the pleiotropic regulators CodY and SigH, the known carbon metabolism regulator CcpA, the proline regulator PrdR, the iron regulator Fur, the small regulatory RNA CsrA and potentially the NADH-responsive regulator Rex. Here, we describe the current knowledge of the metabolic principles of energy generation by C. difficile and the underlying gene regulatory scenarios.

Gender Differences in User Satisfaction of Mobile Touch Screen Interfaces: University Students’ Service Sites

The rapid development of mobile user interfaces for students’ websites and the constant utilization of such interfaces by students have witnessed a significant upsurge in growth. However, mobile service providers may lack valuable feedback on user satisfaction, particularly for Arabic users, because the sites are designed and implemented without students’ participation. This paper empirically investigates the user satisfaction of a mobile banner system for the University of Ha’il in Saudi Arabia. Users’ satisfaction was evaluated across six scales: overall reactions, screens, terminology and system information, learning, system capabilities, and technical manuals and online help. A quantitative research method was utilized, involving a questionnaire survey of 235 students. We found that female students have significant concerns about user satisfaction. The paper proposes theoretical and practical implications for future work.

Development of a Chlamydomonas reinhardtii metabolic network dynamic model to describe distinct phenotypes occurring at different CO2 levels

The increase in atmospheric CO2 due to anthropogenic activities is generating climate change, which has resulted in a subsequent rise in global temperatures with severe environmental impacts. Biological mitigation has been considered as an alternative for environmental remediation and reduction of greenhouse gases in the atmosphere. In fact, the use of easily adapted photosynthetic organisms able to fix CO2 with low-cost operation is revealing its high potential for industry. Among those organism, the algae Chlamydomonas reinhardtii have gain special attention as a model organism for studying CO2 fixation, biomass accumulation and bioenergy production upon exposure to several environmental conditions. In the present study, we studied the Chlamydomonas response to different CO2 levels by comparing metabolomics and transcriptomics data with the predicted results from our new-improved genomic-scale metabolic model. For this, we used in silico methods at steady dynamic state varying the levels of CO2. Our main goal was to improve our capacity for predicting metabolic routes involved in biomass accumulation. The improved genomic-scale metabolic model presented in this study was shown to be phenotypically accurate, predictive, and a significant improvement over previously reported models. Our model consists of 3726 reactions and 2436 metabolites, and lacks any thermodynamically infeasible cycles. It was shown to be highly sensitive to environmental changes under both steady-state and dynamic conditions. As additional constraints, our dynamic model involved kinetic parameters associated with substrate consumption at different growth conditions (i.e., low CO2-heterotrophic and high CO2-mixotrophic). Our results suggest that cells growing at high CO2 (i.e., photoautotrophic and mixotrophic conditions) have an increased capability for biomass production. In addition, we have observed that ATP production also seems to be an important limiting factor for growth under the conditions tested. Our experimental data (metabolomics and transcriptomics) and the results predicted by our model clearly suggest a differential behavior between low CO2-heterotrophic and high CO2-mixotrophic growth conditions. The data presented in the current study contributes to better dissect the biological response of C. reinhardtii, as a dynamic entity, to environmental and genetic changes. These findings are of great interest given the biotechnological potential of this microalga for CO2 fixation, biomass accumulation, and bioenergy production.

Understanding regulation in substrate dependent modulation of growth and production of alcohols in Clostridium sporogenes NCIM 2918 through metabolic network reconstruction and flux balance analysis

Flux Balance Analysis was performed for Clostridium sporogenes NCIM 2918 grown on sole glucose and glycerol or glucose-glycerol combinations at varied concentrations. During acidogenesis, glucose and glucose-glycerol combinations favored improved growth and butyric acid production. Glycerol fermentation was however marked by reduced growth and predominant ethanol synthesis. Further, with increase of glycerol fraction in glucose-glycerol blend, flux towards ethanol synthesis linearly increased with simultaneous decrease in butanol flux. Elevated ATP demand due to improved growth was satisfied by upregulated carbon flux towards butyric acid synthesis during both glucose and dual substrate fermentations. Possible repression of pyruvate carboxylase by glycerol resulting in downturn of carbon uptake flux towards TCA cycle through anaplerotic reaction may be responsible for reduced growth in glycerol fermentation. Ammonium acetate mediated induction of acetic acid utilization, during acidogenesis, led to excess acetyl-CoA generation and its subsequent metabolism to lesser reduced products, butyric acid or ethanol.

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

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

Amino acid catabolism-directed biofuel production in Clostridium sticklandii: An insight into model-driven systems engineering

•

Model-driven systems engineering has been more fascinating process for microbial biofuel production.

•

Clostridium sticklandii is a potential strain for the solventogenesis and acidogenesis.

•

The present review provides an insight for the protein catabolism-directed biofuel production.

Model-driven systems engineering has been more fascinating process for the microbial production of biofuel and bio-refineries in chemical and pharmaceutical industries. Genome-scale modeling and simulations have been guided for metabolic engineering of Clostridium species for the production of organic solvents and organic acids. Among them, Clostridium sticklandii is one of the potential organisms to be exploited as a microbial cell factory for biofuel production. It is a hyper-ammonia producing bacterium and is able to catabolize amino acids as important carbon and energy sources via Stickland reactions and the development of the specific pathways. Current genomic and metabolic aspects of this bacterium are comprehensively reviewed herein, which provided information for learning about protein catabolism-directed biofuel production. It has a metabolic potential to drive energy and direct solventogenesis as well as acidogenesis from protein catabolism. It produces by-products such as ethanol, acetate, n-butanol, n-butyrate and hydrogen from amino acid catabolism. Model-driven systems engineering of this organism would improve the performance of the industrial sectors and enhance the industrial economy by using protein-based waste in environment-friendly ways.

Clostridium butyricum maximizes growth while minimizing enzyme usage and ATP production: metabolic flux distribution of a strain cultured in glycerol

BACKGROUND

The increase in glycerol obtained as a byproduct of biodiesel has encouraged the production of new industrial products, such as 1,3-propanediol (PDO), using biotechnological transformation via bacteria like Clostridium butyricum. However, despite the increasing role of Clostridium butyricum as a bio-production platform, its metabolism remains poorly modeled.

RESULTS

We reconstructed iCbu641, the first genome-scale metabolic (GSM) model of a PDO producer Clostridium strain, which included 641 genes, 365 enzymes, 891 reactions, and 701 metabolites. We found an enzyme expression prediction of nearly 84% after comparison of proteomic data with flux distribution estimation using flux balance analysis (FBA). The remaining 16% corresponded to enzymes directionally coupled to growth, according to flux coupling findings (FCF). The fermentation data validation also revealed different phenotype states that depended on culture media conditions; for example, Clostridium maximizes its biomass yield per enzyme usage under glycerol limitation. By contrast, under glycerol excess conditions, Clostridium grows sub-optimally, maximizing biomass yield while minimizing both enzyme usage and ATP production. We further evaluated perturbations in the GSM model through enzyme deletions and variations in biomass composition. The GSM predictions showed no significant increase in PDO production, suggesting a robustness to perturbations in the GSM model. We used the experimental results to predict that co-fermentation was a better alternative than iCbu641 perturbations for improving PDO yields.

CONCLUSIONS

The agreement between the predicted and experimental values allows the use of the GSM model constructed for the PDO-producing Clostridium strain to propose new scenarios for PDO production, such as dynamic simulations, thereby reducing the time and costs associated with experimentation.

Model-based quantification of metabolic interactions from dynamic microbial-community data

An important challenge in microbial ecology is to infer metabolic-exchange fluxes between growing microbial species from community-level data, concerning species abundances and metabolite concentrations. Here we apply a model-based approach to integrate such experimental data and thereby infer metabolic-exchange fluxes. We designed a synthetic anaerobic co-culture of Clostridium acetobutylicum and Wolinella succinogenes that interact via interspecies hydrogen transfer and applied different environmental conditions for which we expected the metabolic-exchange rates to change. We used stoichiometric models of the metabolism of the two microorganisms that represents our current physiological understanding and found that this understanding - the model - is sufficient to infer the identity and magnitude of the metabolic-exchange fluxes and it suggested unexpected interactions. Where the model could not fit all experimental data, it indicates specific requirement for further physiological studies. We show that the nitrogen source influences the rate of interspecies hydrogen transfer in the co-culture. Additionally, the model can predict the intracellular fluxes and optimal metabolic exchange rates, which can point to engineering strategies. This study therefore offers a realistic illustration of the strengths and weaknesses of model-based integration of heterogenous data that makes inference of metabolic-exchange fluxes possible from community-level experimental data.

Mathematical modelling of clostridial acetone-butanol-ethanol fermentation

Clostridial acetone-butanol-ethanol (ABE) fermentation features a remarkable shift in the cellular metabolic activity from acid formation, acidogenesis, to the production of industrial-relevant solvents, solventogensis. In recent decades, mathematical models have been employed to elucidate the complex interlinked regulation and conditions that determine these two distinct metabolic states and govern the transition between them. In this review, we discuss these models with a focus on the mechanisms controlling intra- and extracellular changes between acidogenesis and solventogenesis. In particular, we critically evaluate underlying model assumptions and predictions in the light of current experimental knowledge. Towards this end, we briefly introduce key ideas and assumptions applied in the discussed modelling approaches, but waive a comprehensive mathematical presentation. We distinguish between structural and dynamical models, which will be discussed in their chronological order to illustrate how new biological information facilitates the ‘evolution’ of mathematical models. Mathematical models and their analysis have significantly contributed to our knowledge of ABE fermentation and the underlying regulatory network which spans all levels of biological organization. However, the ties between the different levels of cellular regulation are not well understood. Furthermore, contradictory experimental and theoretical results challenge our current notion of ABE metabolic network structure. Thus, clostridial ABE fermentation still poses theoretical as well as experimental challenges which are best approached in close collaboration between modellers and experimentalists.

Enhanced butanol fermentation using metabolically engineered Clostridium acetobutylicum with ex situ recovery of butanol

In this study, metabolic target reactions for strain engineering were searched via intracellular coenzyme A (CoA) metabolite analysis. The metabolic reactions catalyzed by thiolase (AtoB) and aldehyde-alcohol dehydrogenase (AdhE1) were considered potential rate-limiting steps. In addition, CoA transferase (CtfAB) was highlighted as being important for the assimilation of organic acids, in order to achieve high butanol production. Based on this quantitative analysis, the BEKW_E1AB-atoB strain was constructed by overexpressing the thl (atoB), adhE1, and ctfAB genes in Clostridium acetobutylicum strain BEKW, which has the phosphotransacetylase (pta) and butyrate kinase (buk) genes knocked out. After 100h of continuous fermentation coupled with adsorptive ex situ butanol recovery, the concentrations found after considering desorption, yield, and productivity for the BEKW_E1AB-atoB strain were 55.7g/L, 0.38g/g, and 2.64g/L/h, respectively. The level of butanol production achieved (2.64g/L/h) represents the highest reported value obtained after adsorptive, long-term fermentation.

Biomass, strain engineering, and fermentation processes for butanol production by solventogenic clostridia

Butanol is considered an attractive biofuel and a commercially important bulk chemical. However, economical production of butanol by solventogenic clostridia, e.g., via fermentative production of acetone-butanol-ethanol (ABE), is hampered by low fermentation performance, mainly as a result of toxicity of butanol to microorganisms and high substrate costs. Recently, sugars from marine macroalgae and syngas were recognized as potent carbon sources in biomass feedstocks that are abundant and do not compete for arable land with edible crops. With the aid of systems metabolic engineering, many researchers have developed clostridial strains with improved performance on fermentation of these substrates. Alternatively, fermentation strategies integrated with butanol recovery processes such as adsorption, gas stripping, liquid-liquid extraction, and pervaporation have been designed to increase the overall titer of butanol and volumetric productivity. Nevertheless, for economically feasible production of butanol, innovative strategies based on recent research should be implemented. This review describes and discusses recent advances in the development of biomass feedstocks, microbial strains, and fermentation processes for butanol production.

Investigating Moorella thermoacetica metabolism with a genome-scale constraint-based metabolic model

Moorella thermoacetica is a strictly anaerobic, endospore-forming, and metabolically versatile acetogenic bacterium capable of conserving energy by both autotrophic (acetogenesis) and heterotrophic (homoacetogenesis) modes of metabolism. Its metabolic diversity and the ability to efficiently convert a wide range of compounds, including syngas (CO + H2) into acetyl-CoA have made this thermophilic bacterium a promising host for industrial biotechnology applications. However, lack of detailed information on M. thermoacetica's metabolism is a major impediment to its use as a microbial cell factory. In order to overcome this issue, a genome-scale constraint-based metabolic model of Moorella thermoacetica, iAI558, has been developed using its genome sequence and physiological data from published literature. The reconstructed metabolic network of M. thermoacetica comprises 558 metabolic genes, 705 biochemical reactions, and 698 metabolites. Of the total 705 model reactions, 680 are gene-associated while the rest are non-gene associated reactions. The model, in addition to simulating both autotrophic and heterotrophic growth of M. thermoacetica, revealed degeneracy in its TCA-cycle, a common characteristic of anaerobic metabolism. Furthermore, the model helped elucidate the poorly understood energy conservation mechanism of M. thermoacetica during autotrophy. Thus, in addition to generating experimentally testable hypotheses regarding its physiology, such a detailed model will facilitate rapid strain designing and metabolic engineering of M. thermoacetica for industrial applications.

System-level modeling of acetone-butanol-ethanol fermentation

Acetone-butanol-ethanol (ABE) fermentation is a metabolic process of clostridia that produces bio-based solvents including butanol. It is enabled by an underlying metabolic reaction network and modulated by cellular gene regulation and environmental cues. Mathematical modeling has served as a valuable strategy to facilitate the understanding, characterization and optimization of this process. In this review, we highlight recent advances in system-level, quantitative modeling of ABE fermentation. We begin with an overview of integrative processes underlying the fermentation. Next we survey modeling efforts including early simple models, models with a systematic metabolic description, and those incorporating metabolism through simple gene regulation. Particular focus is given to a recent system-level model that integrates the metabolic reactions, gene regulation and environmental cues. We conclude by discussing the remaining challenges and future directions towards predictive understanding of ABE fermentation.

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.

In silico analysis of Clostridium acetobutylicum ATCC 824 metabolic response to an external electron supply

The biological production of butanol has become an important research field and thanks to genome sequencing and annotation; genome-scale metabolic reconstructions have been developed for several Clostridium species. This work makes use of the iCAC490 model of Clostridium acetobutylicum ATCC 824 to analyze its metabolic capabilities and response to an external electron supply through a constraint-based approach using the Constraint-Based Reconstruction Analysis Toolbox. Several analyses were conducted, which included sensitivity, production envelope, and phenotypic phase planes. The model showed that the use of an external electron supply, which acts as co-reducing agent along with glucose-derived reducing power (electrofermentation), results in an increase in the butanol-specific productivity. However, a proportional increase in the butyrate uptake flux is required. Besides, the uptake of external butyrate leads to the coupling of butanol production and growth, which coincides with results reported in literature. Phenotypic phase planes showed that the reducing capacity becomes more limiting for growth at high butyrate uptake fluxes. An electron uptake flux allows the metabolism to reach the growth optimality line. Although the maximum butanol flux does not coincide with the growth optimality line, a butyrate uptake combined with an electron uptake flux would result in an increased butanol volumetric productivity, being a potential strategy to optimize the production of butanol by C. acetobutylicum ATCC 824.

Genome-scale modelling of microbial metabolism with temporal and spatial resolution

Most natural microbial systems have evolved to function in environments with temporal and spatial variations. A major limitation to understanding such complex systems is the lack of mathematical modelling frameworks that connect the genomes of individual species and temporal and spatial variations in the environment to system behaviour. The goal of this review is to introduce the emerging field of spatiotemporal metabolic modelling based on genome-scale reconstructions of microbial metabolism. The extension of flux balance analysis (FBA) to account for both temporal and spatial variations in the environment is termed spatiotemporal FBA (SFBA). Following a brief overview of FBA and its established dynamic extension, the SFBA problem is introduced and recent progress is described. Three case studies are reviewed to illustrate the current state-of-the-art and possible future research directions are outlined. The author posits that SFBA is the next frontier for microbial metabolic modelling and a rapid increase in methods development and system applications is anticipated.

A Quantitative System-Scale Characterization of the Metabolism of Clostridium acetobutylicum

Engineering industrial microorganisms for ambitious applications, for example, the production of second-generation biofuels such as butanol, is impeded by a lack of knowledge of primary metabolism and its regulation. A quantitative system-scale analysis was applied to the biofuel-producing bacterium Clostridium acetobutylicum, a microorganism used for the industrial production of solvent. An improved genome-scale model, iCac967, was first developed based on thorough biochemical characterizations of 15 key metabolic enzymes and on extensive literature analysis to acquire accurate fluxomic data. In parallel, quantitative transcriptomic and proteomic analyses were performed to assess the number of mRNA molecules per cell for all genes under acidogenic, solventogenic, and alcohologenic steady-state conditions as well as the number of cytosolic protein molecules per cell for approximately 700 genes under at least one of the three steady-state conditions. A complete fluxomic, transcriptomic, and proteomic analysis applied to different metabolic states allowed us to better understand the regulation of primary metabolism. Moreover, this analysis enabled the functional characterization of numerous enzymes involved in primary metabolism, including (i) the enzymes involved in the two different butanol pathways and their cofactor specificities, (ii) the primary hydrogenase and its redox partner, (iii) the major butyryl coenzyme A (butyryl-CoA) dehydrogenase, and (iv) the major glyceraldehyde-3-phosphate dehydrogenase. This study provides important information for further metabolic engineering of C. acetobutylicum to develop a commercial process for the production of n-butanol.

Currently, there is a resurgence of interest in Clostridium acetobutylicum, the biocatalyst of the historical Weizmann process, to produce n-butanol for use both as a bulk chemical and as a renewable alternative transportation fuel. To develop a commercial process for the production of n-butanol via a metabolic engineering approach, it is necessary to better characterize both the primary metabolism of C. acetobutylicum and its regulation. Here, we apply a quantitative system-scale analysis to acidogenic, solventogenic, and alcohologenic steady-state C. acetobutylicum cells and report for the first time quantitative transcriptomic, proteomic, and fluxomic data. This approach allows for a better understanding of the regulation of primary metabolism and for the functional characterization of numerous enzymes involved in primary metabolism.

Parallel labeling experiments for pathway elucidation and (13)C metabolic flux analysis

Metabolic pathway models provide the foundation for quantitative studies of cellular physiology through the measurement of intracellular metabolic fluxes. For model organisms metabolic models are well established, with many manually curated genome-scale model reconstructions, gene knockout studies and stable-isotope tracing studies. However, for non-model organisms a similar level of knowledge is often lacking. Compartmentation of cellular metabolism in eukaryotic systems also presents significant challenges for quantitative (13)C-metabolic flux analysis ((13)C-MFA). Recently, innovative (13)C-MFA approaches have been developed based on parallel labeling experiments, the use of multiple isotopic tracers and integrated data analysis, that allow more rigorous validation of pathway models and improved quantification of metabolic fluxes. Applications of these approaches open new research directions in metabolic engineering, biotechnology and medicine.

Salmonella typhimurium and Escherichia coli dissimilarity: Closely related bacteria with distinct metabolic profiles

Live attenuated strains of Salmonella typhimurium have been extensively investigated as vaccines for a number of infectious diseases. However, there is still little information available concerning aspects of their metabolism. S. typhimurium and Escherichia coli show a high degree of similarity in terms of their genome contents and metabolic networks. However, this work presents experimental evidence showing that significant differences exist in their abilities to direct carbon fluxes to biomass and energy production. It is important to study the metabolism of Salmonella to elucidate the formation of acetate and other metabolites involved in optimizing the production of biomass, essential for the development of recombinant vaccines. The metabolism of Salmonella under aerobic conditions was assessed using continuous cultures performed at dilution rates ranging from 0.1 to 0.67 h(-1), with glucose as main substrate. Acetate assimilation and glucose metabolism under anaerobic conditions were also investigated using batch cultures. Chemostat cultivations showed deviation of carbon towards acetate formation, starting at dilution rates above 0.1 h(-1). This differed from previous findings for E. coli, where acetate accumulation was only detected at dilution rates exceeding 0.4 h(-1), and was due to the lower rate of acetate assimilation by S. typhimurium under aerobic conditions. Under anaerobic conditions, both microorganisms mainly produced ethanol, acetate, and formate. A genome-scale metabolic model, reconstructed for Salmonella based on an E. coli model, provided a poor description of the mixed fermentation pattern observed during Salmonella cultures, reinforcing the different patterns of carbon utilization exhibited by these closely related bacteria.

Predicting internal cell fluxes at sub-optimal growth

Background

Flux Balance Analysis (FBA) is a widely used tool to model metabolic behavior and cellular function. Applications of FBA span a breadth of research from synthetic engineering of biofuels to understanding evolutionary adaptations. FBA predicts metabolic reaction fluxes that optimize a given objective. This objective is generally defined for unicellular organisms by a theoretical reaction which simulates biomass production. FBA has been extremely successful at predicting in E. coli growth rates under different media and gene essentiality, amongst other things. In order to improve predictions, additional constraints are coupled with optimization of the biomass function. Studies have suggested, however, that unicellular organisms - like multicellular organisms - do not grow at optimal rates. To further improve FBA predictions, particularly of internal cell fluxes, new techniques to explore the sub-optimal solution space need to be developed.

Results

We present an innovative FBA method called corsoFBA based on the optimization of protein cost at sub-optimal objective levels. Our method shows good agreement with experimental data of E. coli grown at different dilution rates. Maintaining the objective function close to its maximum value predicts metabolic states that closely resemble low dilution rates; while higher dilution rates can be mirrored by lowering the biomass production value. By using a modified version of Extreme Pathways, we are also able to quantify the energy production and overall protein cost for all possible pathways in the central carbon metabolism.

Conclusion

Metabolic flux distributions at the optimal objective can be substantially different from the near-optimal distributions. Importantly, the behavior of E. coli central carbon metabolism can be better predicted by exploring the sub-optimal FBA solution space. The corsoFBA method presented here is able to predict the behavior of PEP Carboxylase, the glyoxylate shunt and the Entner-Doudoroff pathway at different glucose levels, a behavior not predicted by the minimization of metabolic steps and FBA alone. This technique can be used to better predict internal cell fluxes under different conditions, and corsoFBA will be of great help for the study of cells from multicellular organisms using Flux Balance Analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0153-3) contains supplementary material, which is available to authorized users.

Flux balance analysis of different carbon source fermentation with hydrogen producing Clostridium butyricum using Cell Net Analyzer

A metabolic network model for Clostridium butyricum was developed using six different carbon sources (sucrose, fructose, galactose, mannose, trehalose and ribose) to study the fermentative H2 production. The model was used for investigation of H2 production and the ability of growth on different substrates to predict the maximum abilities for fermentative H2 production that each substrate can support. NADH fluxes were calculated by the model as an important cofactor affecting on H2 production. Butyrate and acetate production were used as model assumptions and biomass formation was chosen as the objective function for flux analysis calculations. Among the substrates selected, sucrose and trehalose supported the maximum growth and H2 yields. The Cell Net Analyzer metabolic network model developed for H2 estimation revealed good correlation with experimental data and could be further used to study the effect of environmental conditions and substrates concentration on H2 yield.

Low-Carbon Fuel and Chemical Production by Anaerobic Gas Fermentation

World energy demand is expected to increase by up to 40% by 2035. Over this period, the global population is also expected to increase by a billion people. A challenge facing the global community is not only to increase the supply of fuel, but also to minimize fossil carbon emissions to safeguard the environment, at the same time as ensuring that food production and supply is not detrimentally impacted. Gas fermentation is a rapidly maturing technology which allows low carbon fuel and commodity chemical synthesis. Unlike traditional biofuel technologies, gas fermentation avoids the use of sugars, relying instead on gas streams rich in carbon monoxide and/or hydrogen and carbon dioxide as sources of carbon and energy for product synthesis by specialized bacteria collectively known as acetogens. Thus, gas fermentation enables access to a diverse array of novel, large volume, and globally available feedstocks including industrial waste gases and syngas produced, for example, via the gasification of municipal waste and biomass. Through the efforts of academic labs and early stage ventures, process scale-up challenges have been surmounted through the development of specialized bioreactors. Furthermore, tools for the genetic improvement of the acetogenic bacteria have been reported, paving the way for the production of a spectrum of ever-more valuable products via this process. As a result of these developments, interest in gas fermentation among both researchers and legislators has grown significantly in the past 5 years to the point that this approach is now considered amongst the mainstream of emerging technology solutions for near-term low-carbon fuel and chemical synthesis.

A genome-scale metabolic model of Methanococcus maripaludis S2 for CO2 capture and conversion to methane

Methane is a major energy source for heating and electricity. Its production by methanogenic bacteria is widely known in nature. M. maripaludis S2 is a fully sequenced hydrogenotrophic methanogen and an excellent laboratory strain with robust genetic tools. However, a quantitative systems biology model to complement these tools is absent in the literature. To understand and enhance its methanogenesis from CO2, this work presents the first constraint-based genome-scale metabolic model (iMM518). It comprises 570 reactions, 556 distinct metabolites, and 518 genes along with gene-protein-reaction (GPR) associations, and covers 30% of open reading frames (ORFs). The model was validated using biomass growth data and experimental phenotypic studies from the literature. Its comparison with the in silico models of Methanosarcina barkeri, Methanosarcina acetivorans, and Sulfolobus solfataricus P2 shows M. maripaludis S2 to be a better organism for producing methane. Using the model, genes essential for growth were identified, and the efficacies of alternative carbon, hydrogen and nitrogen sources were studied. The model can predict the effects of reengineering M. maripaludis S2 to guide or expedite experimental efforts.

A curated C. difficile strain 630 metabolic network: prediction of essential targets and inhibitors

BACKGROUND

Clostridium difficile is the leading cause of hospital-borne infections occurring when the natural intestinal flora is depleted following antibiotic treatment. Current treatments for Clostridium difficile infections present high relapse rates and new hyper-virulent and multi-resistant strains are emerging, making the study of this nosocomial pathogen necessary to find novel therapeutic targets.

RESULTS

We present iMLTC806cdf, an extensively curated reconstructed metabolic network for the C. difficile pathogenic strain 630. iMLTC806cdf contains 806 genes, 703 metabolites and 769 metabolic, 117 exchange and 145 transport reactions. iMLTC806cdf is the most complete and accurate metabolic reconstruction of a gram-positive anaerobic bacteria to date. We validate the model with simulated growth assays in different media and carbon sources and use it to predict essential genes. We obtain 89.2% accuracy in the prediction of gene essentiality when compared to experimental data for B. subtilis homologs (the closest organism for which such data exists). We predict the existence of 76 essential genes and 39 essential gene pairs, a number of which are unique to C. difficile and have non-existing or predicted non-essential human homologs. For 29 of these potential therapeutic targets, we find 125 inhibitors of homologous proteins including approved drugs with the potential for drug repositioning, that when validated experimentally could serve as starting points in the development of new antibiotics.

CONCLUSIONS

We created a highly curated metabolic network model of C. difficile strain 630 and used it to predict essential genes as potential new therapeutic targets in the fight against Clostridium difficile infections.

Parallel labeling experiments validate Clostridium acetobutylicum metabolic network model for (13)C metabolic flux analysis

In this work, we provide new insights into the metabolism of Clostridium acetobutylicum ATCC 824 obtained using a systematic approach for quantifying fluxes based on parallel labeling experiments and (13)C-metabolic flux analysis ((13)C-MFA). Here, cells were grown in parallel cultures with [1-(13)C]glucose and [U-(13)C]glucose as tracers and (13)C-MFA was used to quantify intracellular metabolic fluxes. Several metabolic network models were compared: an initial model based on current knowledge, and extended network models that included additional reactions that improved the fits of experimental data. While the initial network model did not produce a statistically acceptable fit of (13)C-labeling data, an extended network model with five additional reactions was able to fit all data with 292 redundant measurements. The model was subsequently trimmed to produce a minimal network model of C. acetobutylicum for (13)C-MFA, which could still reproduce all of the experimental data. The flux results provided valuable new insights into the metabolism of C. acetobutylicum. First, we found that TCA cycle was effectively incomplete, as there was no measurable flux between α-ketoglutarate and succinyl-CoA, succinate and fumarate, and malate and oxaloacetate. Second, an active pathway was identified from pyruvate to fumarate via aspartate. Third, we found that isoleucine was produced exclusively through the citramalate synthase pathway in C. acetobutylicum and that CAC3174 was likely responsible for citramalate synthase activity. These model predictions were confirmed in several follow-up tracer experiments. The validated metabolic network model established in this study can be used in future investigations for unbiased (13)C-flux measurements in C. acetobutylicum.

ITEP: an integrated toolkit for exploration of microbial pan-genomes

Background

Comparative genomics is a powerful approach for studying variation in physiological traits as well as the evolution and ecology of microorganisms. Recent technological advances have enabled sequencing large numbers of related genomes in a single project, requiring computational tools for their integrated analysis. In particular, accurate annotations and identification of gene presence and absence are critical for understanding and modeling the cellular physiology of newly sequenced genomes. Although many tools are available to compare the gene contents of related genomes, new tools are necessary to enable close examination and curation of protein families from large numbers of closely related organisms, to integrate curation with the analysis of gain and loss, and to generate metabolic networks linking the annotations to observed phenotypes.

Results

We have developed ITEP, an Integrated Toolkit for Exploration of microbial Pan-genomes, to curate protein families, compute similarities to externally-defined domains, analyze gene gain and loss, and generate draft metabolic networks from one or more curated reference network reconstructions in groups of related microbial species among which the combination of core and variable genes constitute the their "pan-genomes". The ITEP toolkit consists of: (1) a series of modular command-line scripts for identification, comparison, curation, and analysis of protein families and their distribution across many genomes; (2) a set of Python libraries for programmatic access to the same data; and (3) pre-packaged scripts to perform common analysis workflows on a collection of genomes. ITEP’s capabilities include de novo protein family prediction, ortholog detection, analysis of functional domains, identification of core and variable genes and gene regions, sequence alignments and tree generation, annotation curation, and the integration of cross-genome analysis and metabolic networks for study of metabolic network evolution.

Conclusions

ITEP is a powerful, flexible toolkit for generation and curation of protein families. ITEP's modular design allows for straightforward extension as analysis methods and tools evolve. By integrating comparative genomics with the development of draft metabolic networks, ITEP harnesses the power of comparative genomics to build confidence in links between genotype and phenotype and helps disambiguate gene annotations when they are evaluated in both evolutionary and metabolic network contexts.

Capturing the response of Clostridium acetobutylicum to chemical stressors using a regulated genome-scale metabolic model

BACKGROUND

Clostridia are anaerobic Gram-positive Firmicutes containing broad and flexible systems for substrate utilization, which have been used successfully to produce a range of industrial compounds. In particular, Clostridium acetobutylicum has been used to produce butanol on an industrial scale through acetone-butanol-ethanol (ABE) fermentation. A genome-scale metabolic (GSM) model is a powerful tool for understanding the metabolic capacities of an organism and developing metabolic engineering strategies for strain development. The integration of stress-related specific transcriptomics information with the GSM model provides opportunities for elucidating the focal points of regulation.

RESULTS

We describe here the construction and validation of a GSM model for C. acetobutylicum ATCC 824, iCac802. iCac802 spans 802 genes and includes 1,137 metabolites and 1,462 reactions, along with gene-protein-reaction associations. Both (13)C-MFA and gene deletion data in the ABE fermentation pathway were used to test the predicted flux ranges allowed by the model. We also describe the CoreReg method, introduced in this paper, to integrate transcriptomic data and identify core sets of reactions that, when their flux was selectively restricted, reproduced flux and biomass-formation ranges seen under all regulatory constraints. CoreReg was used in response to butanol and butyrate stress to tighten bounds for 50 reactions within the iCac802 model. These bounds affected the flux of tens of reactions in core metabolism. The model, incorporating the regulatory restrictions from CoreReg under chemical stress, exhibited an approximate 70% reduction in biomass yield for most stress conditions.

CONCLUSIONS

The regulation placed on the model for the two stresses using CoreReg identified differences in the respective responses, including distinct core sets and the restriction of biomass production similar to experimental observations. Given the core sets predicted by the CoreReg method, remedial actions can be taken to counteract the effect of stress on metabolism. For less well-known systems, plausible regulatory loops can be suggested around the affected metabolic reactions, and the hypotheses can be tested experimentally.