Genome-scale modeling using flux ratio constraints to enable metabolic engineering of clostridial metabolism in silico

Simulating compatible solute biosynthesis using a metabolic flux model of the biomining acidophile, Acidithiobacillus ferrooxidans ATCC 23270

Halotolerant, acidophilic, bioleaching microorganisms are crucial to biomining operations that utilize saline water. Compatible solutes play an important role in the adaptation of these microorganisms to saline environments. Acidithiobacillus ferrooxidans ATCC 23270, an iron- and sulfur-oxidizing acidophilic bacterium, synthesizes trehalose as its native compatible solute but is still sensitive to salinity. Recently, halotolerant bioleaching bacteria were found to use ectoine as their key compatible solute. Previously, bioleaching bacteria were recalcitrant to genetic manipulation; however, recent advancements in genetic tools and techniques allow successful genetic modification of A. ferrooxidans ATCC 23270. Therefore, this study aimed to test, in silico, the effect of native and synthetic compatible solute biosynthesis by A. ferrooxidans ATCC 23270 on its growth and metabolism. Metabolic network flux modelling was used to provide a computational framework for the prediction of metabolic fluxes during production of native and synthetic compatible solutes by A. ferrooxidans ATCC 23270, in silico. Complete pathways for trehalose biosynthesis by the bacterium are proposed and captured in the updated metabolic model including a newly discovered UDP-dependent trehalose synthesis pathway. Finally, the effect of nitrogen sources on compatible solute production was simulated and showed that using nitrogen gas as the sole nitrogen source enables the ectoine-producing 'engineered' microbe to oxidize up to 20% more ferrous iron in comparison to the native microbe that only produces trehalose. Therefore, the predictive outcomes of the model have the potential to guide the design and optimization of a halotolerant strain of A. ferrooxidans ATCC 23270 for saline bioleaching operations.

FBA-PRCC. Partial Rank Correlation Coefficient (PRCC) Global Sensitivity Analysis (GSA) in Application to Constraint-Based Models

Background: Whole-genome models (GEMs) have become a versatile tool for systems biology, biotechnology, and medicine. GEMs created by automatic and semi-automatic approaches contain a lot of redundant reactions. At the same time, the nonlinearity of the model makes it difficult to evaluate the significance of the reaction for cell growth or metabolite production. Methods: We propose a new way to apply the global sensitivity analysis (GSA) to GEMs in a straightforward parallelizable fashion. Results: We have shown that Partial Rank Correlation Coefficient (PRCC) captures key steps in the metabolic network despite the network distance from the product synthesis reaction. Conclusions: FBA-PRCC is a fast, interpretable, and reliable metric to identify the sign and magnitude of the reaction contribution to various cellular functions.

CONTRABASS: exploiting flux constraints in genome-scale models for the detection of vulnerabilities

MOTIVATION

Despite the fact that antimicrobial resistance is an increasing health concern, the pace of production of new drugs is slow due to the high cost and uncertain success of the process. The development of high-throughput technologies has allowed the integration of biological data into detailed genome-scale models of multiple organisms. Such models can be exploited by means of computational methods to identify system vulnerabilities such as chokepoint reactions and essential reactions. These vulnerabilities are appealing drug targets that can lead to novel drug developments. However, the current approach to compute these vulnerabilities is only based on topological data and ignores the dynamic information of the model. This can lead to misidentified drug targets.

RESULTS

This work computes flux constraints that are consistent with a certain growth rate of the modelled organism, and integrates the computed flux constraints into the model to improve the detection of vulnerabilities. By exploiting these flux constraints, we are able to obtain a directionality of the reactions of metabolism consistent with a given growth rate of the model, and consequently, a more realistic detection of vulnerabilities can be performed. Several sets of reactions that are system vulnerabilities are defined and the relationships among them are studied. The approach for the detection of these vulnerabilities has been implemented in the Python tool CONTRABASS. Such tool, for which an online web server has also been implemented, computes flux constraints and generates a report with the detected vulnerabilities.

AVAILABILITY AND IMPLEMENTATION

CONTRABASS is available as an open source Python package at https://github.com/openCONTRABASS/CONTRABASS under GPL-3.0 License. An online web server is available at http://contrabass.unizar.es.

SUPPLEMENTARY INFORMATION

A glossary of terms are available at Bioinformatics online.

A Computational Toolbox to Investigate the Metabolic Potential and Resource Allocation in Fission Yeast

The fission yeast, Schizosaccharomyces pombe, is a popular eukaryal model organism for cell division and cell cycle studies. With this extensive knowledge of its cell and molecular biology, S. pombe also holds promise for use in metabolism research and industrial applications. However, unlike the baker's yeast, Saccharomyces cerevisiae, a major workhorse in these areas, cell physiology and metabolism of S. pombe remain less explored. One way to advance understanding of organism-specific metabolism is construction of computational models and their use for hypothesis testing. To this end, we leverage existing knowledge of S. cerevisiae to generate a manually curated high-quality reconstruction of S. pombe's metabolic network, including a proteome-constrained version of the model. Using these models, we gain insights into the energy demands for growth, as well as ribosome kinetics in S. pombe. Furthermore, we predict proteome composition and identify growth-limiting constraints that determine optimal metabolic strategies under different glucose availability regimes and reproduce experimentally determined metabolic profiles. Notably, we find similarities in metabolic and proteome predictions of S. pombe with S. cerevisiae, which indicate that similar cellular resource constraints operate to dictate metabolic organization. With these cases, we show, on the one hand, how these models provide an efficient means to transfer metabolic knowledge from a well-studied to a lesser-studied organism, and on the other, how they can successfully be used to explore the metabolic behavior and the role of resource allocation in driving different strategies in fission yeast. IMPORTANCE Our understanding of microbial metabolism relies mostly on the knowledge we have obtained from a limited number of model organisms, and the diversity of metabolism beyond the handful of model species thus remains largely unexplored in mechanistic terms. Computational modeling of metabolic networks offers an attractive platform to bridge the knowledge gap and gain new insights into physiology of lesser-studied organisms. Here we showcase an example of successful knowledge transfer from the budding yeast Saccharomyces cerevisiae to a popular model organism in molecular and cell biology, fission yeast Schizosaccharomyces pombe, using computational models.

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.

Kinetic model of Clostridium beijerinckii's Acetone-Butanol-Ethanol fermentation considering metabolically diverse cell types

Clostridium beijerinckii population branches into metabolically diverse cell types in batch cultures. Here, we present a new kinetic model of C. beijerinckii's Acetone-Butanol-Ethanol fermentation that considers three cell types: producers of acids (acidogenic), consumer of acids and producers of solvents (solventogenic), and spores cells. The model accurately recapitulates batch culture data. Also, the model estimates cell type-specific kinetic parameters, which can be helpful to improve the operation of the ABE fermentation and give a framework to study acidogenic and solventogenic metabolic pathways. To exemplify the latter, we used a constraint-based model to study how the ABE pathways are used among acidogenic and solventogenic cell types. We found that among both cell types, glycolytic production of ATP and consumption of NAD⁺ varies widely during the fermentation, with their maximum production/consumption rates happening when acidogenic and solventogenic growth rates were at their highest. However, acidogenic cells use the ABE pathway to contribute with an extra 12.5% of the total production of ATP, whereas solventogenic cell types use the ABE pathway to supply more than 75% of the demand for NAD⁺, alternating between the production of lactate and butyrate, being both coupled to the production of NAD⁺.

Towards a Systems Biology Approach to Understanding the Lichen Symbiosis: Opportunities and Challenges of Implementing Network Modelling

Lichen associations, a classic model for successful and sustainable interactions between micro-organisms, have been studied for many years. However, there are significant gaps in our understanding about how the lichen symbiosis operates at the molecular level. This review addresses opportunities for expanding current knowledge on signalling and metabolic interplays in the lichen symbiosis using the tools and approaches of systems biology, particularly network modelling. The largely unexplored nature of symbiont recognition and metabolic interdependency in lichens could benefit from applying a holistic approach to understand underlying molecular mechanisms and processes. Together with ‘omics’ approaches, the application of signalling and metabolic network modelling could provide predictive means to gain insights into lichen signalling and metabolic pathways. First, we review the major signalling and recognition modalities in the lichen symbioses studied to date, and then describe how modelling signalling networks could enhance our understanding of symbiont recognition, particularly leveraging omics techniques. Next, we highlight the current state of knowledge on lichen metabolism. We also discuss metabolic network modelling as a tool to simulate flux distribution in lichen metabolic pathways and to analyse the co-dependence between symbionts. This is especially important given the growing number of lichen genomes now available and improved computational tools for reconstructing such models. We highlight the benefits and possible bottlenecks for implementing different types of network models as applied to the study of lichens.

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.

Reconstructing organisms in silico: genome-scale models and their emerging applications

Escherichia coli is considered to be the best-known microorganism given the large number of published studies detailing its genes, its genome and the biochemical functions of its molecular components. This vast literature has been systematically assembled into a reconstruction of the biochemical reaction networks that underlie E. coli's functions, a process which is now being applied to an increasing number of microorganisms. Genome-scale reconstructed networks are organized and systematized knowledge bases that have multiple uses, including conversion into computational models that interpret and predict phenotypic states and the consequences of environmental and genetic perturbations. These genome-scale models (GEMs) now enable us to develop pan-genome analyses that provide mechanistic insights, detail the selection pressures on proteome allocation and address stress phenotypes. In this Review, we first discuss the overall development of GEMs and their applications. Next, we review the evolution of the most complete GEM that has been developed to date: the E. coli GEM. Finally, we explore three emerging areas in genome-scale modelling of microbial phenotypes: collections of strain-specific models, metabolic and macromolecular expression models, and simulation of stress responses.

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.

Distributed flux balance analysis simulations of serial biomass fermentation by two organisms

Intelligent biorefinery design that addresses both the composition of the biomass feedstock as well as fermentation microorganisms could benefit from dedicated tools for computational simulation and computer-assisted optimization. Here we present the BioLego Vn2.0 framework, based on Microsoft Azure Cloud, which supports large-scale simulations of biomass serial fermentation processes by two different organisms. BioLego enables the simultaneous analysis of multiple fermentation scenarios and the comparison of fermentation potential of multiple feedstock compositions. Thanks to the effective use of cloud computing it further allows resource intensive analysis and exploration of media and organism modifications. We use BioLego to obtain biological and validation results, including (1) exploratory search for the optimal utilization of corn biomasses-corn cobs, corn fiber and corn stover-in fermentation biorefineries; (2) analysis of the possible effects of changes in the composition of K. alvarezi biomass on the ethanol production yield in an anaerobic two-step process (S. cerevisiae followed by E. coli); (3) analysis of the impact, on the estimated ethanol production yield, of knocking out single organism reactions either in one or in both organisms in an anaerobic two-step fermentation process of Ulva sp. into ethanol (S. cerevisiae followed by E. coli); and (4) comparison of several experimentally measured ethanol fermentation rates with the predictions of BioLego.

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.

Seeing the wood for the trees: a forest of methods for optimization and omic-network integration in metabolic modelling

Metabolic modelling has entered a mature phase with dozens of methods and software implementations available to the practitioner and the theoretician. It is not easy for a modeller to be able to see the wood (or the forest) for the trees. Driven by this analogy, we here present a 'forest' of principal methods used for constraint-based modelling in systems biology. This provides a tree-based view of methods available to prospective modellers, also available in interactive version at http://modellingmetabolism.net, where it will be kept updated with new methods after the publication of the present manuscript. Our updated classification of existing methods and tools highlights the most promising in the different branches, with the aim to develop a vision of how existing methods could hybridize and become more complex. We then provide the first hands-on tutorial for multi-objective optimization of metabolic models in R. We finally discuss the implementation of multi-view machine learning approaches in poly-omic integration. Throughout this work, we demonstrate the optimization of trade-offs between multiple metabolic objectives, with a focus on omic data integration through machine learning. We anticipate that the combination of a survey, a perspective on multi-view machine learning and a step-by-step R tutorial should be of interest for both the beginner and the advanced user.

CiliateGEM: an open-project and a tool for predictions of ciliate metabolic variations and experimental condition design

BACKGROUND

The study of cell metabolism is becoming central in several fields such as biotechnology, evolution/adaptation and human disease investigations. Here we present CiliateGEM, the first metabolic network reconstruction draft of the freshwater ciliate Tetrahymena thermophila. We also provide the tools and resources to simulate different growth conditions and to predict metabolic variations. CiliateGEM can be extended to other ciliates in order to set up a meta-model, i.e. a metabolic network reconstruction valid for all ciliates. Ciliates are complex unicellular eukaryotes of presumably monophyletic origin, with a phylogenetic position that is equal from plants and animals. These cells represent a new concept of unicellular system with a high degree of species, population biodiversity and cell complexity. Ciliates perform in a single cell all the functions of a pluricellular organism, including locomotion, feeding, digestion, and sexual processes.

RESULTS

After generating the model, we performed an in-silico simulation with the presence and absence of glucose. The lack of this nutrient caused a 32.1% reduction rate in biomass synthesis. Despite the glucose starvation, the growth did not stop due to the use of alternative carbon sources such as amino acids.

CONCLUSIONS

The future models obtained from CiliateGEM may represent a new approach to describe the metabolism of ciliates. This tool will be a useful resource for the ciliate research community in order to extend these species as model organisms in different research fields. An improved understanding of ciliate metabolism could be relevant to elucidate the basis of biological phenomena like genotype-phenotype relationships, population genetics, and cilia-related disease mechanisms.

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.

Insights into metabolic osmoadaptation of the ectoines-producer bacterium Chromohalobacter salexigens through a high-quality genome scale metabolic model

Background

The halophilic bacterium Chromohalobacter salexigens is a natural producer of ectoines, compatible solutes with current and potential biotechnological applications. As production of ectoines is an osmoregulated process that draws away TCA intermediates, bacterial metabolism needs to be adapted to cope with salinity changes. To explore and use C. salexigens as cell factory for ectoine(s) production, a comprehensive knowledge at the systems level of its metabolism is essential. For this purpose, the construction of a robust and high-quality genome-based metabolic model of C. salexigens was approached.

Results

We generated and validated a high quality genome-based C. salexigens metabolic model (iFP764). This comprised an exhaustive reconstruction process based on experimental information, analysis of genome sequence, manual re-annotation of metabolic genes, and in-depth refinement. The model included three compartments (periplasmic, cytoplasmic and external medium), and two salinity-specific biomass compositions, partially based on experimental results from C. salexigens. Using previous metabolic data as constraints, the metabolic model allowed us to simulate and analyse the metabolic osmoadaptation of C. salexigens under conditions for low and high production of ectoines. The iFP764 model was able to reproduce the major metabolic features of C. salexigens. Flux Balance Analysis (FBA) and Monte Carlo Random sampling analysis showed salinity-specific essential metabolic genes and different distribution of fluxes and variation in the patterns of correlation of reaction sets belonging to central C and N metabolism, in response to salinity. Some of them were related to bioenergetics or production of reducing equivalents, and probably related to demand for ectoines. Ectoines metabolic reactions were distributed according to its correlation in four modules. Interestingly, the four modules were independent both at low and high salinity conditions, as they did not correlate to each other, and they were not correlated with other subsystems.

Conclusions

Our validated model is one of the most complete curated networks of halophilic bacteria. It is a powerful tool to simulate and explore C. salexigens metabolism at low and high salinity conditions, driving to low and high production of ectoines. In addition, it can be useful to optimize the metabolism of other halophilic bacteria for metabolite production.

Electronic supplementary material

The online version of this article (10.1186/s12934-017-0852-0) contains supplementary material, which is available to authorized users.

In silico engineering of Pseudomonas metabolism reveals new biomarkers for increased biosurfactant production

BACKGROUND

Rhamnolipids, biosurfactants with a wide range of biomedical applications, are amphiphilic molecules produced on the surfaces of or excreted extracellularly by bacteria including Pseudomonas aeruginosa. However, Pseudomonas putida is a non-pathogenic model organism with greater metabolic versatility and potential for industrial applications.

METHODS

We investigate in silico the metabolic capabilities of P. putida for rhamnolipids biosynthesis using statistical, metabolic and synthetic engineering approaches after introducing key genes (RhlA and RhlB) from P. aeruginosa into a genome-scale model of P. putida. This pipeline combines machine learning methods with multi-omic modelling, and drives the engineered P. putida model toward an optimal production and export of rhamnolipids out of the membrane.

RESULTS

We identify a substantial increase in synthesis of rhamnolipids by the engineered model compared to the control model. We apply statistical and machine learning techniques on the metabolic reaction rates to identify distinct features on the structure of the variables and individual components driving the variation of growth and rhamnolipids production. We finally provide a computational framework for integrating multi-omics data and identifying latent pathways and genes for the production of rhamnolipids in P. putida.

CONCLUSIONS

We anticipate that our results will provide a versatile methodology for integrating multi-omics data for topological and functional analysis of P. putida toward maximization of biosurfactant production.

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.

Clostridioides difficile 630Δerm in silico and in vivo - quantitative growth and extensive polysaccharide secretion

Antibiotic-associated infections with Clostridioides difficile are a severe and often lethal risk for hospitalized patients, and can also affect populations without these classical risk factors. For a rational design of therapeutical concepts, a better knowledge of the metabolism of the pathogen is crucial. Metabolic modeling can provide a simulation of quantitative growth and usage of metabolic pathways, leading to a deeper understanding of the organism. Here, we present an elaborate genome-scale metabolic model of C. difficile 630Δerm. The model iHD992 includes experimentally determined product and substrate uptake rates and is able to simulate the energy metabolism and quantitative growth of C. difficile. Dynamic flux balance analysis was used for time-resolved simulations of the quantitative growth in two different media. The model predicts oxidative Stickland reactions and glucose degradation as main sources of energy, while the resulting reduction potential is mostly used for acetogenesis via the Wood-Ljungdahl pathway. Initial modeling experiments did not reproduce the observed growth behavior before the production of large quantities of a previously unknown polysaccharide was detected. Combined genome analysis and laboratory experiments indicated that the polysaccharide is an acetylated glucose polymer. Time-resolved simulations showed that polysaccharide secretion was coupled to growth even during unstable glucose uptake in minimal medium. This is accomplished by metabolic shifts between active glycolysis and gluconeogenesis which were also observed in laboratory experiments.

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.

Exploring complex cellular phenotypes and model-guided strain design with a novel genome-scale metabolic model of Clostridium thermocellum DSM 1313 implementing an adjustable cellulosome

BACKGROUND

Clostridium thermocellum is a gram-positive thermophile that can directly convert lignocellulosic material into biofuels. The metabolism of C. thermocellum contains many branches and redundancies which limit biofuel production, and typical genetic techniques are time-consuming. Further, the genome sequence of a genetically tractable strain C. thermocellum DSM 1313 has been recently sequenced and annotated. Therefore, developing a comprehensive, predictive, genome-scale metabolic model of DSM 1313 is desired for elucidating its complex phenotypes and facilitating model-guided metabolic engineering.

RESULTS

We constructed a genome-scale metabolic model iAT601 for DSM 1313 using the KEGG database as a scaffold and an extensive literature review and bioinformatic analysis for model refinement. Next, we used several sets of experimental data to train the model, e.g., estimation of the ATP requirement for growth-associated maintenance (13.5 mmol ATP/g DCW/h) and cellulosome synthesis (57 mmol ATP/g cellulosome/h). Using our tuned model, we investigated the effect of cellodextrin lengths on cell yields, and could predict in silico experimentally observed differences in cell yield based on which cellodextrin species is assimilated. We further employed our tuned model to analyze the experimentally observed differences in fermentation profiles (i.e., the ethanol to acetate ratio) between cellobiose- and cellulose-grown cultures and infer regulatory mechanisms to explain the phenotypic differences. Finally, we used the model to design over 250 genetic modification strategies with the potential to optimize ethanol production, 6155 for hydrogen production, and 28 for isobutanol production.

CONCLUSIONS

Our developed genome-scale model iAT601 is capable of accurately predicting complex cellular phenotypes under a variety of conditions and serves as a high-quality platform for model-guided strain design and metabolic engineering to produce industrial biofuels and chemicals of interest.

Thermochemical hydrolysis of macroalgae Ulva for biorefinery: Taguchi robust design method

Understanding the impact of all process parameters on the efficiency of biomass hydrolysis and on the final yield of products is critical to biorefinery design. Using Taguchi orthogonal arrays experimental design and Partial Least Square Regression, we investigated the impact of change and the comparative significance of thermochemical process temperature, treatment time, %Acid and %Solid load on carbohydrates release from green macroalgae from Ulva genus, a promising biorefinery feedstock. The average density of hydrolysate was determined using a new microelectromechanical optical resonator mass sensor. In addition, using Flux Balance Analysis techniques, we compared the potential fermentation yields of these hydrolysate products using metabolic models of Escherichia coli, Saccharomyces cerevisiae wild type, Saccharomyces cerevisiae RN1016 with xylose isomerase and Clostridium acetobutylicum. We found that %Acid plays the most significant role and treatment time the least significant role in affecting the monosaccharaides released from Ulva biomass. We also found that within the tested range of parameters, hydrolysis with 121 °C, 30 min 2% Acid, 15% Solids could lead to the highest yields of conversion: 54.134–57.500 gr ethanol kg⁻¹Ulva dry weight by S. cerevisiae RN1016 with xylose isomerase. Our results support optimized marine algae utilization process design and will enable smart energy harvesting by thermochemical hydrolysis.

Synthetic biology and regulatory networks: where metabolic systems biology meets control engineering

Metabolic pathways can be engineered to maximize the synthesis of various products of interest. With the advent of computational systems biology, this endeavour is usually carried out through in silico theoretical studies with the aim to guide and complement further in vitro and in vivo experimental efforts. Clearly, what counts is the result in vivo, not only in terms of maximal productivity but also robustness against environmental perturbations. Engineering an organism towards an increased production flux, however, often compromises that robustness. In this contribution, we review and investigate how various analytical approaches used in metabolic engineering and synthetic biology are related to concepts developed by systems and control engineering. While trade-offs between production optimality and cellular robustness have already been studied diagnostically and statically, the dynamics also matter. Integration of the dynamic design aspects of control engineering with the more diagnostic aspects of metabolic, hierarchical control and regulation analysis is leading to the new, conceptual and operational framework required for the design of robust and productive dynamic pathways.

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.

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.

Perturbation Experiments: Approaches for Metabolic Pathway Analysis in Bioreactors

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

CBFA: phenotype prediction integrating metabolic models with constraints derived from experimental data

Background

Flux analysis methods lie at the core of Metabolic Engineering (ME), providing methods for phenotype simulation that allow the determination of flux distributions under different conditions. Although many constraint-based modeling software tools have been developed and published, none provides a free user-friendly application that makes available the full portfolio of flux analysis methods.

Results

This work presents Constraint-based Flux Analysis (CBFA), an open-source software application for flux analysis in metabolic models that implements several methods for phenotype prediction, allowing users to define constraints associated with measured fluxes and/or flux ratios, together with environmental conditions (e.g. media) and reaction/gene knockouts. CBFA identifies the set of applicable methods based on the constraints defined from user inputs, encompassing algebraic and constraint-based simulation methods. The integration of CBFA within the OptFlux framework for ME enables the utilization of different model formats and standards and the integration with complementary methods for phenotype simulation and visualization of results.

Conclusions

A general-purpose and flexible application is proposed that is independent of the origin of the constraints defined for a given simulation. The aim is to provide a simple to use software tool focused on the application of several flux prediction methods.

Electronic supplementary material

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

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.

A review of metabolic and enzymatic engineering strategies for designing and optimizing performance of microbial cell factories

Microbial cell factories (MCFs) are of considerable interest to convert low value renewable substrates to biofuels and high value chemicals. This review highlights the progress of computational models for the rational design of an MCF to produce a target bio-commodity. In particular, the rational design of an MCF involves: (i) product selection, (ii) de novo biosynthetic pathway identification (i.e., rational, heterologous, or artificial), (iii) MCF chassis selection, (iv) enzyme engineering of promiscuity to enable the formation of new products, and (v) metabolic engineering to ensure optimal use of the pathway by the MCF host. Computational tools such as (i) de novo biosynthetic pathway builders, (ii) docking, (iii) molecular dynamics (MD) and steered MD (SMD), and (iv) genome-scale metabolic flux modeling all play critical roles in the rational design of an MCF. Genome-scale metabolic flux models are of considerable use to the design process since they can reveal metabolic capabilities of MCF hosts. These can be used for host selection as well as optimizing precursors and cofactors of artificial de novo biosynthetic pathways. In addition, recent advances in genome-scale modeling have enabled the derivation of metabolic engineering strategies, which can be implemented using the genomic tools reviewed here as well.

Coenzyme A-transferase-independent butyrate re-assimilation in Clostridium acetobutylicum-evidence from a mathematical model

The hetero-dimeric CoA-transferase CtfA/B is believed to be crucial for the metabolic transition from acidogenesis to solventogenesis in Clostridium acetobutylicum as part of the industrial-relevant acetone-butanol-ethanol (ABE) fermentation. Here, the enzyme is assumed to mediate re-assimilation of acetate and butyrate during a pH-induced metabolic shift and to faciliate the first step of acetone formation from acetoacetyl-CoA. However, recent investigations using phosphate-limited continuous cultures have questioned this common dogma. To address the emerging experimental discrepancies, we investigated the mutant strain Cac-ctfA398s::CT using chemostat cultures. As a consequence of this mutation, the cells are unable to express functional ctfA and are thus lacking CoA-transferase activity. A mathematical model of the pH-induced metabolic shift, which was recently developed for the wild type, is used to analyse the observed behaviour of the mutant strain with a focus on re-assimilation activities for the two produced acids. Our theoretical analysis reveals that the ctfA mutant still re-assimilates butyrate, but not acetate. Based upon this finding, we conclude that C. acetobutylicum possesses a CoA-tranferase-independent butyrate uptake mechanism that is activated by decreasing pH levels. Furthermore, we observe that butanol formation is not inhibited under our experimental conditions, as suggested by previous batch culture experiments. In concordance with recent batch experiments, acetone formation is abolished in chemostat cultures using the ctfa mutant.

Computational approaches to metabolic engineering utilizing systems biology and synthetic biology

Metabolic engineering modifies cellular function to address various biochemical applications. Underlying metabolic engineering efforts are a host of tools and knowledge that are integrated to enable successful outcomes. Concurrent development of computational and experimental tools has enabled different approaches to metabolic engineering. One approach is to leverage knowledge and computational tools to prospectively predict designs to achieve the desired outcome. An alternative approach is to utilize combinatorial experimental tools to empirically explore the range of cellular function and to screen for desired traits. This mini-review focuses on computational systems biology and synthetic biology tools that can be used in combination for prospective in silico strain design.

Metabolic changes in Klebsiella oxytoca in response to low oxidoreduction potential, as revealed by comparative proteomic profiling integrated with flux balance analysis

Oxidoreduction potential (ORP) is an important physiological parameter for biochemical production in anaerobic or microaerobic processes. However, the effect of ORP on cellular physiology remains largely unknown, which hampers the design of engineering strategies targeting proteins associated with ORP response. Here we characterized the effect of altering ORP in a 1,3-propanediol producer, Klebsiella oxytoca, by comparative proteomic profiling combined with flux balance analysis. Decreasing the extracellular ORP from -150 to -240 mV retarded cell growth and enhanced 1,3-propanediol production. Comparative proteomic analysis identified 61 differentially expressed proteins, mainly involved in carbohydrate catabolism, cellular constituent biosynthesis, and reductive stress response. A hypothetical oxidoreductase (HOR) that catalyzes 1,3-propanediol production was markedly upregulated, while proteins involved in biomass precursor synthesis were downregulated. As revealed by subsequent flux balance analysis, low ORP induced a metabolic shift from glycerol oxidation to reduction and rebalancing of redox and energy metabolism. From the integrated protein expression profiles and flux distributions, we can construct a rational analytic framework that elucidates how (facultative) anaerobes respond to extracellular ORP changes.

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.

Integration of a constraint-based metabolic model of Brassica napus developing seeds with (13)C-metabolic flux analysis

The use of large-scale or genome-scale metabolic reconstructions for modeling and simulation of plant metabolism and integration of those models with large-scale omics and experimental flux data is becoming increasingly important in plant metabolic research. Here we report an updated version of bna572, a bottom-up reconstruction of oilseed rape (Brassica napus L.; Brassicaceae) developing seeds with emphasis on representation of biomass-component biosynthesis. New features include additional seed-relevant pathways for isoprenoid, sterol, phenylpropanoid, flavonoid, and choline biosynthesis. Being now based on standardized data formats and procedures for model reconstruction, bna572+ is available as a COBRA-compliant Systems Biology Markup Language (SBML) model and conforms to the Minimum Information Requested in the Annotation of Biochemical Models (MIRIAM) standards for annotation of external data resources. Bna572+ contains 966 genes, 671 reactions, and 666 metabolites distributed among 11 subcellular compartments. It is referenced to the Arabidopsis thaliana genome, with gene-protein-reaction (GPR) associations resolving subcellular localization. Detailed mass and charge balancing and confidence scoring were applied to all reactions. Using B. napus seed specific transcriptome data, expression was verified for 78% of bna572+ genes and 97% of reactions. Alongside bna572+ we also present a revised carbon centric model for (13)C-Metabolic Flux Analysis ((13)C-MFA) with all its reactions being referenced to bna572+ based on linear projections. By integration of flux ratio constraints obtained from (13)C-MFA and by elimination of infinite flux bounds around thermodynamically infeasible loops based on COBRA loopless methods, we demonstrate improvements in predictive power of Flux Variability Analysis (FVA). Using this combined approach we characterize the difference in metabolic flux of developing seeds of two B. napus genotypes contrasting in starch and oil content.

Constraint-based genome-scale metabolic modeling of Clostridium acetobutylicum behavior in an immobilized column

In this study a step-wise optimization procedure was developed to predict solvent production using continuous ABE fermentation with immobilized cells. The modeling approach presented here utilizes previously published constraint-based metabolic model for Clostridium acetobutylicum without direct flux constraints. A recently developed flux ratio constraint method was adopted for the model. An experimental data set consisting of 25 experiments using different sugar mixtures as substrates and differing dilution rates was simulated successfully with the modeling approach. Converted to end product concentrations the mean absolute error for acetone was 0.31 g/l, for butanol 0.49 g/l, and for ethanol 0.17 g/l. The modeling approach was validated with another data set from similar experimental setup. The model errors for the validation data set was 0.24 g/l, 0.60 g/l, and 0.17 g/l for acetone, butanol, and ethanol, respectively.

A metabolite-centric view on flux distributions in genome-scale metabolic models

Background

Genome-scale metabolic models are important tools in systems biology. They permit the in-silico prediction of cellular phenotypes via mathematical optimisation procedures, most importantly flux balance analysis. Current studies on metabolic models mostly consider reaction fluxes in isolation. Based on a recently proposed metabolite-centric approach, we here describe a set of methods that enable the analysis and interpretation of flux distributions in an integrated metabolite-centric view. We demonstrate how this framework can be used for the refinement of genome-scale metabolic models.

Results

We applied the metabolite-centric view developed here to the most recent metabolic reconstruction of Escherichia coli. By compiling the balance sheets of a small number of currency metabolites, we were able to fully characterise the energy metabolism as predicted by the model and to identify a possibility for model refinement in NADPH metabolism. Selected branch points were examined in detail in order to demonstrate how a metabolite-centric view allows identifying functional roles of metabolites. Fructose 6-phosphate aldolase and the sedoheptulose bisphosphate bypass were identified as enzymatic reactions that can carry high fluxes in the model but are unlikely to exhibit significant activity in vivo. Performing a metabolite essentiality analysis, unconstrained import and export of iron ions could be identified as potentially problematic for the quality of model predictions.

Conclusions

The system-wide analysis of split ratios and branch points allows a much deeper insight into the metabolic network than reaction-centric analyses. Extending an earlier metabolite-centric approach, the methods introduced here establish an integrated metabolite-centric framework for the interpretation of flux distributions in genome-scale metabolic networks that can complement the classical reaction-centric framework. Analysing fluxes and their metabolic context simultaneously opens the door to systems biological interpretations that are not apparent from isolated reaction fluxes. Particularly powerful demonstrations of this are the analyses of the complete metabolic contexts of energy metabolism and the folate-dependent one-carbon pool presented in this work. Finally, a metabolite-centric view on flux distributions can guide the refinement of metabolic reconstructions for specific growth scenarios.