Metabolic capabilities of Escherichia coli: I. synthesis of biosynthetic precursors and cofactors

Review: Towards the next-generation models of the rumen microbiome for enhancing predictive power and guiding sustainable production strategies

The rumen ecosystem harbours a galaxy of microbes working in syntrophy to carry out a metabolic cascade of hydrolytic and fermentative reactions. This fermentation process allows ruminants to harvest nutrients from a wide range of feedstuff otherwise inaccessible to the host. The interconnection between the ruminant and its rumen microbiota shapes key animal phenotypes such as feed efficiency and methane emissions and suggests the potential of reducing methane emissions and enhancing feed conversion into animal products by manipulating the rumen microbiota. Whilst significant technological progress in omics techniques has increased our knowledge of the rumen microbiota and its genome (microbiome), translating omics knowledge into effective microbial manipulation strategies remains a great challenge. This challenge can be addressed by modelling approaches integrating causality principles and thus going beyond current correlation-based approaches applied to analyse rumen microbial genomic data. However, existing rumen models are not yet adapted to capitalise on microbial genomic information. This gap between the rumen microbiota available omics data and the way microbial metabolism is represented in the existing rumen models needs to be filled to enhance rumen understanding and produce better predictive models with capabilities for guiding nutritional strategies. To fill this gap, the integration of computational biology tools and mathematical modelling frameworks is needed to translate the information of the metabolic potential of the rumen microbes (inferred from their genomes) into a mathematical object. In this paper, we aim to discuss the potential use of two modelling approaches for the integration of microbial genomic information into dynamic models. The first modelling approach explores the theory of state observers to integrate microbial time series data into rumen fermentation models. The second approach is based on the genome-scale network reconstructions of rumen microbes. For a given microorganism, the network reconstruction produces a stoichiometry matrix of the metabolism. This matrix is the core of the so-called genome-scale metabolic models which can be exploited by a plethora of methods comprised within the constraint-based reconstruction and analysis approaches. We will discuss how these methods can be used to produce the next-generation models of the rumen microbiome.

Trade-offs between the instantaneous growth rate and long-term fitness: Consequences for microbial physiology and predictive computational models

Microbial systems biology has made enormous advances in relating microbial physiology to the underlying biochemistry and molecular biology. By meticulously studying model microorganisms, in particular Escherichia coli and Saccharomyces cerevisiae, increasingly comprehensive computational models predict metabolic fluxes, protein expression, and growth. The modeling rationale is that cells are constrained by a limited pool of resources that they allocate optimally to maximize fitness. As a consequence, the expression of particular proteins is at the expense of others, causing trade-offs between cellular objectives such as instantaneous growth, stress tolerance, and capacity to adapt to new environments. While current computational models are remarkably predictive for E. coli and S. cerevisiae when grown in laboratory environments, this may not hold for other growth conditions and other microorganisms. In this contribution, we therefore discuss the relationship between the instantaneous growth rate, limited resources, and long-term fitness. We discuss uses and limitations of current computational models, in particular for rapidly changing and adverse environments, and propose to classify microbial growth strategies based on Grimes's CSR framework.

Simulated Design-Build-Test-Learn Cycles for Consistent Comparison of Machine Learning Methods in Metabolic Engineering

Combinatorial pathway optimization is an important tool in metabolic flux optimization. Simultaneous optimization of a large number of pathway genes often leads to combinatorial explosions. Strain optimization is therefore often performed using iterative design-build-test-learn (DBTL) cycles. The aim of these cycles is to develop a product strain iteratively, every time incorporating learning from the previous cycle. Machine learning methods provide a potentially powerful tool to learn from data and propose new designs for the next DBTL cycle. However, due to the lack of a framework for consistently testing the performance of machine learning methods over multiple DBTL cycles, evaluating the effectiveness of these methods remains a challenge. In this work, we propose a mechanistic kinetic model-based framework to test and optimize machine learning for iterative combinatorial pathway optimization. Using this framework, we show that gradient boosting and random forest models outperform the other tested methods in the low-data regime. We demonstrate that these methods are robust for training set biases and experimental noise. Finally, we introduce an algorithm for recommending new designs using machine learning model predictions. We show that when the number of strains to be built is limited, starting with a large initial DBTL cycle is favorable over building the same number of strains for every cycle.

A neural-mechanistic hybrid approach improving the predictive power of genome-scale metabolic models

Constraint-based metabolic models have been used for decades to predict the phenotype of microorganisms in different environments. However, quantitative predictions are limited unless labor-intensive measurements of media uptake fluxes are performed. We show how hybrid neural-mechanistic models can serve as an architecture for machine learning providing a way to improve phenotype predictions. We illustrate our hybrid models with growth rate predictions of Escherichia coli and Pseudomonas putida grown in different media and with phenotype predictions of gene knocked-out Escherichia coli mutants. Our neural-mechanistic models systematically outperform constraint-based models and require training set sizes orders of magnitude smaller than classical machine learning methods. Our hybrid approach opens a doorway to enhancing constraint-based modeling: instead of constraining mechanistic models with additional experimental measurements, our hybrid models grasp the power of machine learning while fulfilling mechanistic constrains, thus saving time and resources in typical systems biology or biological engineering projects.

Reconstructing Kinetic Models for Dynamical Studies of Metabolism using Generative Adversarial Networks

Kinetic models of metabolism relate metabolic fluxes, metabolite concentrations and enzyme levels through mechanistic relations, rendering them essential for understanding, predicting and optimizing the behaviour of living organisms. However, due to the lack of kinetic data, traditional kinetic modelling often yields only a few or no kinetic models with desirable dynamical properties, making the analysis unreliable and computationally inefficient. We present REKINDLE (Reconstruction of Kinetic Models using Deep Learning), a deep-learning-based framework for efficiently generating kinetic models with dynamic properties matching the ones observed in cells. We showcase REKINDLE's capabilities to navigate through the physiological states of metabolism using small numbers of data with significantly lower computational requirements. The results show that data-driven neural networks assimilate implicit kinetic knowledge and structure of metabolic networks and generate kinetic models with tailored properties and statistical diversity. We anticipate that our framework will advance our understanding of metabolism and accelerate future research in biotechnology and health.

Quantitative prediction of conditional vulnerabilities in regulatory and metabolic networks using PRIME

The ability of Mycobacterium tuberculosis (Mtb) to adopt heterogeneous physiological states underlies its success in evading the immune system and tolerating antibiotic killing. Drug tolerant phenotypes are a major reason why the tuberculosis (TB) mortality rate is so high, with over 1.8 million deaths annually. To develop new TB therapeutics that better treat the infection (faster and more completely), a systems-level approach is needed to reveal the complexity of network-based adaptations of Mtb. Here, we report a new predictive model called PRIME (Phenotype of Regulatory influences Integrated with Metabolism and Environment) to uncover environment-specific vulnerabilities within the regulatory and metabolic networks of Mtb. Through extensive performance evaluations using genome-wide fitness screens, we demonstrate that PRIME makes mechanistically accurate predictions of context-specific vulnerabilities within the integrated regulatory and metabolic networks of Mtb, accurately rank-ordering targets for potentiating treatment with frontline drugs.

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

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

Searching for principles of microbial physiology

Why do evolutionarily distinct microorganisms display similar physiological behaviours? Why are transitions from high-ATP yield to low(er)-ATP yield metabolisms so widespread across species? Why is fast growth generally accompanied with low stress tolerance? Do these regularities occur because most microbial species are subject to the same selective pressures and physicochemical constraints? If so, a broadly-applicable theory might be developed that predicts common microbiological behaviours. Microbial systems biologists have been working out the contours of this theory for the last two decades, guided by experimental data. At its foundations lie basic principles from evolutionary biology, enzyme biochemistry, metabolism, cell composition and steady-state growth. The theory makes predictions about fitness costs and benefits of protein expression, physicochemical constraints on cell growth and characteristics of optimal metabolisms that maximise growth rate. Comparisons of the theory with experimental data indicates that microorganisms often aim for maximisation of growth rate, also in the presence of stresses; they often express optimal metabolisms and metabolic proteins at optimal concentrations. This review explains the current status of the theory for microbiologists; its roots, predictions, experimental evidence and future directions.

Co-feeding glucose with either gluconate or galacturonate during clostridial fermentations provides metabolic fine-tuning capabilities

Clostridium acetobutylicum ATCC 824 effectively utilizes a wide range of substrates to produce commodity chemicals. When grown on substrates of different oxidation states, the organism exhibits different recycling needs of reduced intracellular electron carrying co-factors. Ratios of substrates with different oxidation states were used to modulate the need to balance electron carriers and demonstrate fine-tuned control of metabolic output. Three different oxidized substrates were first fed singularly, then in different ratios to three different strains of Clostridium sp. Growth was most robust when fed glucose in exclusive fermentations. However, the use of the other two more oxidized substrates was strain-dependent in exclusive feeds. In glucose-galacturonate mixed fermentation, the main products (acetate and butyrate) were dependant on the ratios of the substrates. Exclusive fermentation on galacturonate was nearly homoacetic. Co-utilization of galacturonate and glucose was observed from the onset of fermentation in growth conditions using both substrates combined, with the proportion of galacturonate present dictating the amount of acetate produced. For all three strains, increasing galacturonate content (%) in a mixture of galacturonate and glucose from 0 to 50, and 100, resulted in a corresponding increase in the amount of acetate produced. For example, C. acetobutylicum increased from ~ 10 mM to ~ 17 mM, and then ~ 23 mM. No co-utilization was observed when galacturonate was replaced with gluconate in the two substrate co-feed.

Curation and Analysis of a Saccharomyces cerevisiae Genome-Scale Metabolic Model for Predicting Production of Sensory Impact Molecules under Enological Conditions

One approach for elucidating strain-to-strain metabolic differences is the use of genome-scale metabolic models (GSMMs). To date GSMMs have not focused on the industrially important area of flavor production and, as such; do not cover all the pathways relevant to flavor formation in yeast. Moreover, current models for Saccharomyces cerevisiae generally focus on carbon-limited and/or aerobic systems, which is not pertinent to enological conditions. Here, we curate a GSMM (iWS902) to expand on the existing Ehrlich pathway and ester formation pathways central to aroma formation in industrial winemaking, in addition to the existing sulfur metabolism and medium-chain fatty acid (MCFA) pathways that also contribute to production of sensory impact molecules. After validating the model using experimental data, we predict key differences in metabolism for a strain (EC 1118) in two distinct growth conditions, including differences for aroma impact molecules such as acetic acid, tryptophol, and hydrogen sulfide. Additionally, we propose novel targets for metabolic engineering for aroma profile modifications employing flux variability analysis with the expanded GSMM. The model provides mechanistic insights into the key metabolic pathways underlying aroma formation during alcoholic fermentation and provides a potential framework to contribute to new strategies to optimize the aroma of wines.

Providing new insights on the biphasic lifestyle of the predatory bacterium Bdellovibrio bacteriovorus through genome-scale metabolic modeling

In this study we analyze the growth-phase dependent metabolic states of Bdellovibrio bacteriovorus by constructing a fully compartmented, mass and charge-balanced genome-scale metabolic model of this predatory bacterium (iCH457). Considering the differences between life cycle phases driving the growth of this predator, growth-phase condition-specific models have been generated allowing the systematic study of its metabolic capabilities. Using these computational tools, we have been able to analyze, from a system level, the dynamic metabolism of the predatory bacteria as the life cycle progresses. We provide computational evidences supporting potential axenic growth of B. bacteriovorus's in a rich medium based on its encoded metabolic capabilities. Our systems-level analysis confirms the presence of "energy-saving" mechanisms in this predator as well as an abrupt metabolic shift between the attack and intraperiplasmic growth phases. Our results strongly suggest that predatory bacteria's metabolic networks have low robustness, likely hampering their ability to tackle drastic environmental fluctuations, thus being confined to stable and predictable habitats. Overall, we present here a valuable computational testbed based on predatory bacteria activity for rational design of novel and controlled biocatalysts in biotechnological/clinical applications.

From bag-of-genes to bag-of-genomes: metabolic modelling of communities in the era of metagenome-assembled genomes

Metagenomic sequencing of complete microbial communities has greatly enhanced our understanding of the taxonomic composition of microbiotas. This has led to breakthrough developments in bioinformatic disciplines such as assembly, gene clustering, metagenomic binning of species genomes and the discovery of an incredible, so far undiscovered, taxonomic diversity. However, functional annotations and estimating metabolic processes from single species - or communities - is still challenging. Earlier approaches relied mostly on inferring the presence of key enzymes for metabolic pathways in the whole metagenome, ignoring the genomic context of such enzymes, resulting in the 'bag-of-genes' approach to estimate functional capacities of microbiotas. Here, we review recent developments in metagenomic bioinformatics, with a special focus on emerging technologies to simulate and estimate metabolic information, that can be derived from metagenomic assembled genomes. Genome-scale metabolic models can be used to model the emergent properties of microbial consortia and whole communities, and the progress in this area is reviewed. While this subfield of metagenomics is still in its infancy, it is becoming evident that there is a dire need for further bioinformatic tools to address the complex combinatorial problems in modelling the metabolism of large communities as a 'bag-of-genomes'.

Genome-Scale Metabolic Model Driven Design of a Defined Medium for Campylobacter jejuni M1cam

Campylobacter jejuni, the most frequent cause of food-borne bacterial gastroenteritis, is a fastidious organism when grown in the laboratory. Oxygen is required for growth, despite the presence of the metabolic mechanism for anaerobic respiration. Amino acid auxotrophies are variably reported and energy metabolism can occur through several electron donor/acceptor combinations. Overall, the picture is one of a flexible, but vulnerable metabolism. To understand Campylobacter metabolism, we have constructed a fully curated, metabolic model for the reference organism M1 (our variant is M1cam) and validated it through laboratory experiments. Our results show that M1cam is auxotrophic for methionine, niacinamide, and pantothenate. There are complete biosynthesis pathways for all amino acids except methionine and it can produce energy, but not biomass, in the absence of oxygen. M1cam will grow in DMEM/F-12 defined media but not in the previously published Campylobacter specific defined media tested. Using the model, we identified potential auxotrophies and substrates that may improve growth. With this information, we designed simple defined media containing inorganic salts, the auxotrophic substrates, L-methionine, niacinamide, and pantothenate, pyruvate and additional amino acids L-cysteine, L-serine, and L-glutamine for growth enhancement. Our defined media supports a 1.75-fold higher growth rate than Brucella broth after 48 h at 37°C and sustains the growth of other Campylobacter jejuni strains. This media can be used to design reproducible assays that can help in better understanding the adaptation, stress resistance, and the virulence mechanisms of this pathogen. We have shown that with a well-curated metabolic model it is possible to design a media to grow this fastidious organism. This has implications for the investigation of new Campylobacter species defined through metagenomics, such as C. infans.

LK-DFBA: a linear programming-based modeling strategy for capturing dynamics and metabolite-dependent regulation in metabolism

BACKGROUND

The systems-scale analysis of cellular metabolites, "metabolomics," provides data ideal for applications in metabolic engineering. However, many of the computational tools for strain design are built around Flux Balance Analysis (FBA), which makes assumptions that preclude direct integration of metabolomics data into the underlying models. Finding a way to retain the advantages of FBA's linear structure while relaxing some of its assumptions could allow us to account for metabolite levels and metabolite-dependent regulation in strain design tools built from FBA, improving the accuracy of predictions made by these tools. We designed, implemented, and characterized a modeling strategy based on Dynamic FBA (DFBA), called Linear Kinetics-Dynamic Flux Balance Analysis (LK-DFBA), to satisfy these specifications. Our strategy adds constraints describing the dynamics and regulation of metabolism that are strictly linear. We evaluated LK-DFBA against alternative modeling frameworks using simulated noisy data from a small in silico model and a larger model of central carbon metabolism in E. coli, and compared each framework's ability to recapitulate the original system.

RESULTS

In the smaller model, we found that we could use regression from a dynamic flux estimation (DFE) with an optional non-linear parameter optimization to reproduce metabolite concentration dynamic trends more effectively than an ordinary differential equation model with generalized mass action rate laws when tested under realistic data sampling frequency and noise levels. We observed detrimental effects across all tested modeling approaches when metabolite time course data were missing, but found these effects to be smaller for LK-DFBA in most cases. With the E. coli model, we produced qualitatively reasonable results with similar properties to the smaller model and explored two different parameterization structures that yield trade-offs in computation time and accuracy.

CONCLUSIONS

LK-DFBA allows for calculation of metabolite concentrations and considers metabolite-dependent regulation while still retaining many computational advantages of FBA. This provides the proof-of-principle for a new metabolic modeling framework with the potential to create genome-scale dynamic models and the potential to be applied in strain engineering tools that currently use FBA.

An integrated computational and experimental study to investigate Staphylococcus aureus metabolism

Staphylococcus aureus is a metabolically versatile pathogen that colonizes nearly all organs of the human body. A detailed and comprehensive knowledge of staphylococcal metabolism is essential to understand its pathogenesis. To this end, we have reconstructed and experimentally validated an updated and enhanced genome-scale metabolic model of S. aureus USA300_FPR3757. The model combined genome annotation data, reaction stoichiometry, and regulation information from biochemical databases and previous strain-specific models. Reactions in the model were checked and fixed to ensure chemical balance and thermodynamic consistency. To further refine the model, growth assessment of 1920 nonessential mutants from the Nebraska Transposon Mutant Library was performed, and metabolite excretion profiles of important mutants in carbon and nitrogen metabolism were determined. The growth and no-growth inconsistencies between the model predictions and in vivo essentiality data were resolved using extensive manual curation based on optimization-based reconciliation algorithms. Upon intensive curation and refinements, the model contains 863 metabolic genes, 1379 metabolites (including 1159 unique metabolites), and 1545 reactions including transport and exchange reactions. To improve the accuracy and predictability of the model to environmental changes, condition-specific regulation information curated from the existing knowledgebase was incorporated. These critical additions improved the model performance significantly in capturing gene essentiality, substrate utilization, and metabolite production capabilities and increased the ability to generate model-based discoveries of therapeutic significance. Use of this highly curated model will enhance the functional utility of omics data, and therefore, serve as a resource to support future investigations of S. aureus and to augment staphylococcal research worldwide.

Integration of in vivo experiment with a newly developed model of Staphylococcus aureus metabolism helps explore its metabolic versatility. A multidisciplinary team led by Rajib Saha at the University of Nebraska developed a new genome-scale metabolic model of the multi-drug resistant pathogen S. aureus by combining genome annotation data, reaction stoichiometry, and condition- and mutant-specific regulations from biochemical databases and previous strain-specific models. Extensive manual curation and incorporation of newly generated experimental data on growth and metabolite production improved the accuracy and predictability of the model and increased its ability to generate model-based discoveries of therapeutic significance. Use of this highly curated model will enhance the functional utility of omics data and, therefore, serve as a resource to support future investigations of S. aureus and to augment staphylococcal research worldwide.

d-Mannose Treatment neither Affects Uropathogenic Escherichia coli Properties nor Induces Stable FimH Modifications

Urinary tract infections (UTIs) are mainly caused by uropathogenic Escherichia coli (UPEC). Acute and recurrent UTIs are commonly treated with antibiotics, the efficacy of which is limited by the emergence of antibiotic resistant strains. The natural sugar d-mannose is considered as an alternative to antibiotics due to its ability to mask the bacterial adhesin FimH, thereby preventing its binding to urothelial cells. Despite its extensive use, the possibility that d-mannose exerts "antibiotic-like" activity by altering bacterial growth and metabolism or selecting FimH variants has not been investigated yet. To this aim, main bacterial features of the prototype UPEC strain CFT073 treated with d-mannose were analyzed by standard microbiological methods. FimH functionality was analyzed by yeast agglutination and human bladder cell adhesion assays. Our results indicate that high d-mannose concentrations have no effect on bacterial growth and do not interfere with the activity of different antibiotics. d-mannose ranked as the least preferred carbon source to support bacterial metabolism and growth, in comparison with d-glucose, d-fructose, and l-arabinose. Since small glucose amounts are physiologically detectable in urine, we can conclude that the presence of d-mannose is irrelevant for bacterial metabolism. Moreover, d-mannose removal after long-term exposure did not alter FimH's capacity to bind to mannosylated proteins. Overall, our data indicate that d-mannose is a good alternative in the prevention and treatment of UPEC-related UTIs.

Metabolic Modeling Elucidates the Transactions in the Rumen Microbiome and the Shifts Upon Virome Interactions

The complex microbial ecosystem within the bovine rumen plays a crucial role in host nutrition, health, and environmental impact. However, little is known about the interactions between the functional entities within the system, which dictates the community structure and functional dynamics and host physiology. With the advancements in high-throughput sequencing and mathematical modeling, in silico genome-scale metabolic analysis promises to expand our understanding of the metabolic interplay in the community. In an attempt to understand the interactions between microbial species and the phages inside rumen, a genome-scale metabolic modeling approach was utilized by using key members in the rumen microbiome (a bacteroidete, a firmicute, and an archaeon) and the viral phages associated with them. Individual microbial host models were integrated into a community model using multi-level mathematical frameworks. An elaborate and heuristics-based computational procedure was employed to predict previously unknown interactions involving the transfer of fatty acids, vitamins, coenzymes, amino acids, and sugars among the community members. While some of these interactions could be inferred by the available multi-omic datasets, our proposed method provides a systemic understanding of why the interactions occur and how these affect the dynamics in a complex microbial ecosystem. To elucidate the functional role of the virome on the microbiome, local alignment search was used to identify the metabolic functions of the viruses associated with the hosts. The incorporation of these functions demonstrated the role of viral auxiliary metabolic genes in relaxing the metabolic bottlenecks in the microbial hosts and complementing the inter-species interactions. Finally, a comparative statistical analysis of different biologically significant community fitness criteria identified the variation in flux space and robustness of metabolic capacities of the community members. Our elucidation of metabolite exchange among the three members of the rumen microbiome shows how their genomic differences and interactions with the viral strains shape up a highly sophisticated metabolic interplay and explains how such interactions across kingdoms can cause metabolic and compositional shifts in the community and affect the health, nutrition, and pathophysiology of the ruminant animal.

Effects of MCHM on yeast metabolism

On January 2014 approximately 10,000 gallons of crude 4-Methylcyclohexanemethanol (MCHM) and propylene glycol phenol ether (PPH) were accidentally released into the Elk River, West Virginia, contaminating the tap water of around 300,000 residents. Crude MCHM is an industrial chemical used as flotation reagent to clean coal. At the time of the spill, MCHM's toxicological data were limited, an issue that has been addressed by different studies focused on understanding the immediate and long-term effects of MCHM on human health and the environment. Using S. cerevisiae as a model organism we study the effect of acute exposure to crude MCHM on metabolism. Yeasts were treated with MCHM 550 ppm in YPD for 30 minutes. Polar and lipid metabolites were extracted from cells by a chloroform-methanol-water mixture. The extracts were then analyzed by direct injection ESI-MS and by GC-MS. The metabolomics analysis was complemented with flux balance analysis simulations done with genome-scale metabolic network models (GSMNM) of MCHM treated vs non-treated control. We integrated the effect of MCHM on yeast gene expression from RNA-Seq data within these GSMNM. A total of 215 and 73 metabolites were identified by the ESI-MS and GC-MS procedures, respectively. From these 26 and 23 relevant metabolites were selected from ESI-MS and GC-MS respectively, for 49 unique compounds. MCHM induced amino acid accumulation, via its effects on amino acid metabolism, as well as a potential impairment of ribosome biogenesis. MCHM affects phospholipid biosynthesis, with a potential impact on the biophysical properties of yeast cellular membranes. The FBA simulations were able to reproduce the deleterious effect of MCHM on cellular growth and suggest that the effect of MCHM on ubiquinol:ferricytochrome c reductase reaction, caused by the under-expression of CYT1 gene, could be the driven force behind the observed effect on yeast metabolism and growth.

BOFdat: Generating biomass objective functions for genome-scale metabolic models from experimental data

Genome-scale metabolic models (GEMs) are mathematically structured knowledge bases of metabolism that provide phenotypic predictions from genomic information. GEM-guided predictions of growth phenotypes rely on the accurate definition of a biomass objective function (BOF) that is designed to include key cellular biomass components such as the major macromolecules (DNA, RNA, proteins), lipids, coenzymes, inorganic ions and species-specific components. Despite its importance, no standardized computational platform is currently available to generate species-specific biomass objective functions in a data-driven, unbiased fashion. To fill this gap in the metabolic modeling software ecosystem, we implemented BOFdat, a Python package for the definition of a Biomass Objective Function from experimental data. BOFdat has a modular implementation that divides the BOF definition process into three independent modules defined here as steps: 1) the coefficients for major macromolecules are calculated, 2) coenzymes and inorganic ions are identified and their stoichiometric coefficients estimated, 3) the remaining species-specific metabolic biomass precursors are algorithmically extracted in an unbiased way from experimental data. We used BOFdat to reconstruct the BOF of the Escherichia coli model iML1515, a gold standard in the field. The BOF generated by BOFdat resulted in the most concordant biomass composition, growth rate, and gene essentiality prediction accuracy when compared to other methods. Installation instructions for BOFdat are available in the documentation and the source code is available on GitHub (https://github.com/jclachance/BOFdat).

Stoichiometric analysis of the energetics and metabolic impact of photorespiration in C3 plants

Analysis of the impact of photorespiration on plant metabolism is usually based on manual inspection of small network diagrams. Here we create a structural metabolic model that contains the reactions that participate in photorespiration in the plastid, peroxisome, mitochondrion and cytosol, and the metabolite exchanges between them. This model was subjected to elementary flux modes analysis, a technique that enumerates all the component, minimal pathways of a network. Any feasible photorespiratory metabolism in the plant will be some combination of the elementary flux modes (EFMs) that contain the Rubisco oxygenase reaction. Amongst the EFMs we obtained was the classic photorespiratory cycle, but there were also modes that involve photorespiration coupled with mitochondrial metabolism and ATP production, the glutathione-ascorbate cycle and nitrate reduction to ammonia. The modes analysis demonstrated the underlying basis of the metabolic linkages with photorespiration that have been inferred experimentally. The set of reactions common to all the elementary modes showed good agreement with the gene products of mutants that have been reported to have a defective phenotype in photorespiratory conditions. Finally, the set of modes provided a formal demonstration that photorespiration itself does not impact on the CO₂ :O₂ ratio (assimilation quotient), except in those modes associated with concomitant nitrate reduction.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

If we tried to list every known chemical reaction within an organism–human, plant or even bacteria–we would get quite a long and confusing read. But when this information is represented in so-called genome-scale metabolic networks, we have the means to access computationally each of those reactions and their interconnections. Some parts of the network have alternatives, while others are unique and therefore can be essential for growth. Here, we simulate growth and compare essential reactions and genes for the simplest type of unicellular species–prokaryotes–to understand which parts of their metabolism are universally essential and potentially ancestral. We show that similar patterns of essential reactions echo phylogenetic relationships (this makes sense, as the genome provides the building plan for the enzymes that perform those reactions). Our computational predictions correlate strongly with experimental essentiality data. Finally, we show that a crucial step of protein synthesis (tRNA charging) and the synthesis and transformation of small molecules that enzymes require (cofactors) are the most essential and conserved parts of metabolism in prokaryotes. Our results are a step further in understanding the biology and evolution of prokaryotes but can also be relevant in applied studies including metabolic engineering and antibiotic design.

A linked organ-on-chip model of the human neurovascular unit reveals the metabolic coupling of endothelial and neuronal cells

The neurovascular unit (NVU) regulates metabolic homeostasis as well as drug pharmacokinetics and pharmacodynamics in the central nervous system. Metabolic fluxes and conversions over the NVU rely on interactions between brain microvascular endothelium, perivascular pericytes, astrocytes and neurons, making it difficult to identify the contributions of each cell type. Here we model the human NVU using microfluidic organ chips, allowing analysis of the roles of individual cell types in NVU functions. Three coupled chips model influx across the blood-brain barrier (BBB), the brain parenchymal compartment and efflux across the BBB. We used this linked system to mimic the effect of intravascular administration of the psychoactive drug methamphetamine and to identify previously unknown metabolic coupling between the BBB and neurons. Thus, the NVU system offers an in vitro approach for probing transport, efficacy, mechanism of action and toxicity of neuroactive drugs.

A low-complexity metabolic network model for the respiratory and fermentative metabolism of Escherichia coli

Over the last decades, predictive microbiology has made significant advances in the mathematical description of microbial spoiler and pathogen dynamics in or on food products. Recently, the focus of predictive microbiology has shifted from a (semi-)empirical population-level approach towards mechanistic models including information about the intracellular metabolism in order to increase model accuracy and genericness. However, incorporation of this subpopulation-level information increases model complexity and, consequently, the required run time to simulate microbial cell and population dynamics. In this paper, results of metabolic flux balance analyses (FBA) with a genome-scale model are used to calibrate a low-complexity linear model describing the microbial growth and metabolite secretion rates of Escherichia coli as a function of the nutrient and oxygen uptake rate. Hence, the required information about the cellular metabolism (i.e., biomass growth and secretion of cell products) is selected and included in the linear model without incorporating the complete intracellular reaction network. However, the applied FBAs are only representative for microbial dynamics under specific extracellular conditions, viz., a neutral medium without weak acids at a temperature of 37℃. Deviations from these reference conditions lead to metabolic shifts and adjustments of the cellular nutrient uptake or maintenance requirements. This metabolic dependency on extracellular conditions has been taken into account in our low-complex metabolic model. In this way, a novel approach is developed to take the synergistic effects of temperature, pH, and undissociated acids on the cell metabolism into account. Consequently, the developed model is deployable as a tool to describe, predict and control E. coli dynamics in and on food products under various combinations of environmental conditions. To emphasize this point,three specific scenarios are elaborated: (i) aerobic respiration without production of weak acid extracellular metabolites, (ii) anaerobic fermentation with secretion of mixed acid fermentation products into the food environment, and (iii) respiro-fermentative metabolic regimes in between the behaviors at aerobic and anaerobic conditions.

Integrated constraints based analysis of an engineered violacein pathway in Escherichia coli

Strategies towards optimal violacein biosynthesis, a potential drug molecule, need systems level coordination of enzymatic activities of individual genes in a multigene operon vioABCDE. Constraints-based flux balance analysis of an extended iAF1260 model (iAF1260vio) with a reconstructed violacein module predicted growth and violacein yields in Escherichia coli accurately. Shadow price (SP) analysis identified tryptophan metabolism and NADPH as limiting. Increased tryptophan levels in Δpgi & ΔpheA were validated using in silico gene deletion analysis. Phenotypic phase plane (PhPP) analysis highlighted sensitivity between tryptophan and NADPH for violacein synthesis at molar growth yields. A synthetic VioABCDE operon (SYNO) sequence was designed to maximize Codon Adaptive Index (CAI: 0.9) and tune translation initiation rates (TIR: 2-50 fold higher) in E. coli. All pSYN E. coli transformants produced higher violacein, with a maximum six-fold increase in yields. The rational design E. coli: ΔpheA SYN: gave the highest violacein titers (33.8 mg/l). Such integrated approaches targeting multiple molecular hierarchies in the cell can be extended further to increase violacein yields.

DMPy: a Python package for automated mathematical model construction of large-scale metabolic systems

Background

Models of metabolism are often used in biotechnology and pharmaceutical research to identify drug targets or increase the direct production of valuable compounds. Due to the complexity of large metabolic systems, a number of conclusions have been drawn using mathematical methods with simplifying assumptions. For example, constraint-based models describe changes of internal concentrations that occur much quicker than alterations in cell physiology. Thus, metabolite concentrations and reaction fluxes are fixed to constant values. This greatly reduces the mathematical complexity, while providing a reasonably good description of the system in steady state. However, without a large number of constraints, many different flux sets can describe the optimal model and we obtain no information on how metabolite levels dynamically change. Thus, to accurately determine what is taking place within the cell, finer quality data and more detailed models need to be constructed.

Results

In this paper we present a computational framework, DMPy, that uses a network scheme as input to automatically search for kinetic rates and produce a mathematical model that describes temporal changes of metabolite fluxes. The parameter search utilises several online databases to find measured reaction parameters. From this, we take advantage of previous modelling efforts, such as Parameter Balancing, to produce an initial mathematical model of a metabolic pathway. We analyse the effect of parameter uncertainty on model dynamics and test how recent flux-based model reduction techniques alter system properties. To our knowledge this is the first time such analysis has been performed on large models of metabolism. Our results highlight that good estimates of at least 80% of the reaction rates are required to accurately model metabolic systems. Furthermore, reducing the size of the model by grouping reactions together based on fluxes alters the resulting system dynamics.

Conclusion

The presented pipeline automates the modelling process for large metabolic networks. From this, users can simulate their pathway of interest and obtain a better understanding of how altering conditions influences cellular dynamics. By testing the effects of different parameterisations we are also able to provide suggestions to help construct more accurate models of complete metabolic systems in the future.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0584-8) contains supplementary material, which is available to authorized users.

Investigating mixotrophic metabolism in the model diatom Phaeodactylum tricornutum

Diatoms are prominent marine microalgae, interesting not only from an ecological point of view, but also for their possible use in biotechnology applications. They can be cultivated in phototrophic conditions, using sunlight as the sole energy source. Some diatoms, however, can also grow in a mixotrophic mode, wherein both light and external reduced carbon contribute to biomass accumulation. In this study, we investigated the consequences of mixotrophy on the growth and metabolism of the pennate diatom Phaeodactylum tricornutum, using glycerol as the source of reduced carbon. Transcriptomics, metabolomics, metabolic modelling and physiological data combine to indicate that glycerol affects the central-carbon, carbon-storage and lipid metabolism of the diatom. In particular, provision of glycerol mimics typical responses of nitrogen limitation on lipid metabolism at the level of triacylglycerol accumulation and fatty acid composition. The presence of glycerol, despite provoking features reminiscent of nutrient limitation, neither diminishes photosynthetic activity nor cell growth, revealing essential aspects of the metabolic flexibility of these microalgae and suggesting possible biotechnological applications of mixotrophy.This article is part of the themed issue 'The peculiar carbon metabolism in diatoms'.

Computational Approaches on Stoichiometric and Kinetic Modeling for Efficient Strain Design

Engineering biological systems that are capable of overproducing products of interest is the ultimate goal of any biotechnology application. To this end, stoichiometric (or steady state) and kinetic models are increasingly becoming available for a variety of organisms including prokaryotes, eukaryotes, and microbial communities. This ever-accelerating pace of such model reconstructions has also spurred the development of optimization-based strain design techniques. This chapter highlights a number of such frameworks developed in recent years in order to generate testable hypotheses (in terms of genetic interventions), thus addressing the challenges in metabolic engineering. In particular, three major methods are covered in detail including two methods for designing strains (i.e., one stoichiometric model-based and the other by integrating kinetic information into a stoichiometric model) and one method for analyzing microbial communities.

Reconstruction and Analysis of Central Metabolism in Microbes

Genome-scale metabolic models (GEMs) generated from automated reconstruction pipelines often lack accuracy due to the need for extensive gapfilling and the inference of periphery metabolic pathways based on lower-confidence annotations. The central carbon pathways and electron transport chains are among the most well-understood regions of microbial metabolism, and these pathways contribute significantly toward defining cellular behavior and growth conditions. Thus, it is often useful to construct a simplified core metabolic model (CMM) that is comprised of only the high-confidence central pathways. In this chapter, we discuss methods for producing core metabolic models (CMM) based on genome annotations. With its reduced scope compared to GEMs, CMM reconstruction focuses on accurate representation of the central metabolic pathways related to energy biosynthesis and accurate energy yield predictions. We demonstrate the reconstruction and analysis of CMMs using the DOE Systems Biology Knowledgebase (KBase). The complete workflow is available at http://kbase.us/core-models/.

lumpGEM: Systematic generation of subnetworks and elementally balanced lumped reactions for the biosynthesis of target metabolites

In the post-genomic era, Genome-scale metabolic networks (GEMs) have emerged as invaluable tools to understand metabolic capabilities of organisms. Different parts of these metabolic networks are defined as subsystems/pathways, which are sets of functional roles to implement a specific biological process or structural complex, such as glycolysis and TCA cycle. Subsystem/pathway definition is also employed to delineate the biosynthetic routes that produce biomass building blocks. In databases, such as MetaCyc and SEED, these representations are composed of linear routes from precursors to target biomass building blocks. However, this approach cannot capture the nested, complex nature of GEMs. Here we implemented an algorithm, lumpGEM, which generates biosynthetic subnetworks composed of reactions that can synthesize a target metabolite from a set of defined core precursor metabolites. lumpGEM captures balanced subnetworks, which account for the fate of all metabolites along the synthesis routes, thus encapsulating reactions from various subsystems/pathways to balance these metabolites in the metabolic network. Moreover, lumpGEM collapses these subnetworks into elementally balanced lumped reactions that specify the cost of all precursor metabolites and cofactors. It also generates alternative subnetworks and lumped reactions for the same metabolite, accounting for the flexibility of organisms. lumpGEM is applicable to any GEM and any target metabolite defined in the network. Lumped reactions generated by lumpGEM can be also used to generate properly balanced reduced core metabolic models.

Stoichiometric models have been used in the area of metabolic engineering and systems biology for many decades. The early examples of these models include simplified ad hoc built metabolic pathways, and biomass compositions. The development of genome scale models (GEMs) brought a standard to metabolic network modeling. However, the vast amount of detailed biochemistry in GEMs makes it necessary to develop methods to manage the complexity in them. In this study, we developed lumpGEM, a tool that can systematically identify subnetworks from metabolic networks that can perform certain tasks, such as biosynthesis of a biomass building block and any other target metabolite. By generating alternative subnetworks, lumpGEM also accounts for the redundancy in metabolic networks. We applied lumpGEM on latest E. coli GEM iJO1366 and identified subnetworks/lumped reactions for every biomass building block defined in its biomass formulation. We also compared the results from lumpGEM with existing knowledge in the literature. The lumped reactions generated by lumpGEM can be used to generate consistently reduced metabolic network models.

Ethanol production improvement driven by genome-scale metabolic modeling and sensitivity analysis in Scheffersomyces stipitis

The yeast Scheffersomyces stipitis naturally produces ethanol from xylose, however reaching high ethanol yields is strongly dependent on aeration conditions. It has been reported that changes in the availability of NAD(H/+) cofactors can improve fermentation in some microorganisms. In this work genome-scale metabolic modeling and phenotypic phase plane analysis were used to characterize metabolic response on a range of uptake rates. Sensitivity analysis was used to assess the effect of ARC on ethanol production indicating that modifying ARC by inhibiting the respiratory chain ethanol production can be improved. It was shown experimentally in batch culture using Rotenone as an inhibitor of the mitochondrial NADH dehydrogenase complex I (CINADH), increasing ethanol yield by 18%. Furthermore, trajectories for uptakes rates, specific productivity and specific growth rate were determined by modeling the batch culture, to calculate ARC associated to the addition of CINADH inhibitor. Results showed that the increment in ethanol production via respiratory inhibition is due to excess in ARC, which generates an increase in ethanol production. Thus ethanol production improvement could be predicted by a change in ARC.

Perspectives and Challenges in Microbial Communities Metabolic Modeling

Bacteria have evolved to efficiently interact each other, forming complex entities known as microbial communities. These "super-organisms" play a central role in maintaining the health of their eukaryotic hosts and in the cycling of elements like carbon and nitrogen. However, despite their crucial importance, the mechanisms that influence the functioning of microbial communities and their relationship with environmental perturbations are obscure. The study of microbial communities was boosted by tremendous advances in sequencing technologies, and in particular by the possibility to determine genomic sequences of bacteria directly from environmental samples. Indeed, with the advent of metagenomics, it has become possible to investigate, on a previously unparalleled scale, the taxonomical composition and the functional genetic elements present in a specific community. Notwithstanding, the metagenomic approach per se suffers some limitations, among which the impossibility of modeling molecular-level (e.g., metabolic) interactions occurring between community members, as well as their effects on the overall stability of the entire system. The family of constraint-based methods, such as flux balance analysis, has been fruitfully used to translate genome sequences in predictive, genome-scale modeling platforms. Although these techniques have been initially developed for analyzing single, well-known model organisms, their recent improvements allowed engaging in multi-organism in silico analyses characterized by a considerable predictive capability. In the face of these advances, here we focus on providing an overview of the possibilities and challenges related to the modeling of metabolic interactions within a bacterial community, discussing the feasibility and the perspectives of this kind of analysis in the (near) future.

A systems-wide understanding of photosynthetic acclimation in algae and higher plants

The ability of phototrophs to colonise different environments relies on robust protection against oxidative stress, a critical requirement for the successful evolutionary transition from water to land. Photosynthetic organisms have developed numerous strategies to adapt their photosynthetic apparatus to changing light conditions in order to optimise their photosynthetic yield, which is crucial for life on Earth to exist. Photosynthetic acclimation is an excellent example of the complexity of biological systems, where highly diverse processes, ranging from electron excitation over protein protonation to enzymatic processes coupling ion gradients with biosynthetic activity, interact on drastically different timescales from picoseconds to hours. Efficient functioning of the photosynthetic apparatus and its protection is paramount for efficient downstream processes, including metabolism and growth. Modern experimental techniques can be successfully integrated with theoretical and mathematical models to promote our understanding of underlying mechanisms and principles. This review aims to provide a retrospective analysis of multidisciplinary photosynthetic acclimation research carried out by members of the Marie Curie Initial Training Project, AccliPhot, placing the results in a wider context. The review also highlights the applicability of photosynthetic organisms for industry, particularly with regards to the cultivation of microalgae. It intends to demonstrate how theoretical concepts can successfully complement experimental studies broadening our knowledge of common principles in acclimation processes in photosynthetic organisms, as well as in the field of applied microalgal biotechnology.

Metabolic regulation is sufficient for global and robust coordination of glucose uptake, catabolism, energy production and growth in Escherichia coli

The metabolism of microorganisms is regulated through two main mechanisms: changes of enzyme capacities as a consequence of gene expression modulation (“hierarchical control”) and changes of enzyme activities through metabolite-enzyme interactions. An increasing body of evidence indicates that hierarchical control is insufficient to explain metabolic behaviors, but the system-wide impact of metabolic regulation remains largely uncharacterized. To clarify its role, we developed and validated a detailed kinetic model of Escherichia coli central metabolism that links growth to environment. Metabolic control analyses confirm that the control is widely distributed across the network and highlight strong interconnections between all the pathways. Exploration of the model solution space reveals that several robust properties emerge from metabolic regulation, from the molecular level (e.g. homeostasis of total metabolite pool) to the overall cellular physiology (e.g. coordination of carbon uptake, catabolism, energy and redox production, and growth), while allowing a large degree of flexibility at most individual metabolic steps. These properties have important physiological implications for E. coli and significantly expand the self-regulating capacities of its metabolism.

Metabolism is a fundamental biochemical process that enables cells to operate and grow by converting nutrients into ‘building blocks’ and energy. Metabolism happens through the work of enzymes, which are encoded by genes. Thus, genes and their regulation are often thought of controlling metabolism, somewhat at the top of a hierarchical control system. However, an increasing body of evidence indicates that metabolism plays an active role in the control of its own operation via a dense network of metabolite-enzyme interactions. The system-wide role of metabolic regulation is hard to dissect and so far remains largely uncharacterized. To better understand its role, we constructed a detailed kinetic model of the carbon and energy metabolism of the bacterium Escherichia coli, a model organism in Systems and Synthetic biology. Model simulations indicate that kinetic considerations of metabolism alone can explain data from hundreds of experiments, without needing to invoke regulation of gene expression. In particular, metabolic regulation is sufficient to coordinate carbon utilization, redox and energy production, and growth, while maintaining local flexibility at individual metabolic steps. These findings indicate that the self-regulating capacities of E. coli metabolism are far more significant than previously expected, and improve our understanding on how cells work.

Metabolic network modeling with model organisms

Flux balance analysis (FBA) with genome-scale metabolic network models (GSMNM) allows systems level predictions of metabolism in a variety of organisms. Different types of predictions with different accuracy levels can be made depending on the applied experimental constraints ranging from measurement of exchange fluxes to the integration of gene expression data. Metabolic network modeling with model organisms has pioneered method development in this field. In addition, model organism GSMNMs are useful for basic understanding of metabolism, and in the case of animal models, for the study of metabolic human diseases. Here, we discuss GSMNMs of most highly used model organisms with the emphasis on recent reconstructions.

Integration of Biomass Formulations of Genome-Scale Metabolic Models with Experimental Data Reveals Universally Essential Cofactors in Prokaryotes

The composition of a cell in terms of macromolecular building blocks and other organic molecules underlies the metabolic needs and capabilities of a species. Although some core biomass components such as nucleic acids and proteins are evident for most species, the essentiality of the pool of other organic molecules, especially cofactors and prosthetic groups, is yet unclear. Here we integrate biomass compositions from 71 manually curated genome-scale models, 33 large-scale gene essentiality datasets, enzyme-cofactor association data and a vast array of publications, revealing universally essential cofactors for prokaryotic metabolism and also others that are specific for phylogenetic branches or metabolic modes. Our results revise predictions of essential genes in Klebsiella pneumoniae and identify missing biosynthetic pathways in models of Mycobacterium tuberculosis. This work provides fundamental insights into the essentiality of organic cofactors and has implications for minimal cell studies as well as for modeling genotype-phenotype relations in prokaryotic metabolic networks.

•

Seventy one biomass equations of manually curated genome-scale metabolic models are compared.

•

Eight classes of universally essential prokaryotic organic cofactors are proposed.

•

Conditionally essential organic cofactors are presented and discussed.

•

Gene essentiality predictions for Klebsiella pneumoniae are revised.

•

A missing essential pathway in models of Mycobacterium tuberculosis is predicted.

Modeling central metabolism and energy biosynthesis across microbial life

Background

Automatically generated bacterial metabolic models, and even some curated models, lack accuracy in predicting energy yields due to poor representation of key pathways in energy biosynthesis and the electron transport chain (ETC). Further compounding the problem, complex interlinking pathways in genome-scale metabolic models, and the need for extensive gapfilling to support complex biomass reactions, often results in predicting unrealistic yields or unrealistic physiological flux profiles.

Results

To overcome this challenge, we developed methods and tools (http://coremodels.mcs.anl.gov) to build high quality core metabolic models (CMM) representing accurate energy biosynthesis based on a well studied, phylogenetically diverse set of model organisms. We compare these models to explore the variability of core pathways across all microbial life, and by analyzing the ability of our core models to synthesize ATP and essential biomass precursors, we evaluate the extent to which the core metabolic pathways and functional ETCs are known for all microbes. 6,600 (80 %) of our models were found to have some type of aerobic ETC, whereas 5,100 (62 %) have an anaerobic ETC, and 1,279 (15 %) do not have any ETC. Using our manually curated ETC and energy biosynthesis pathways with no gapfilling at all, we predict accurate ATP yields for nearly 5586 (70 %) of the models under aerobic and anaerobic growth conditions. This study revealed gaps in our knowledge of the central pathways that result in 2,495 (30 %) CMMs being unable to produce ATP under any of the tested conditions. We then established a methodology for the systematic identification and correction of inconsistent annotations using core metabolic models coupled with phylogenetic analysis.

Conclusions

We predict accurate energy yields based on our improved annotations in energy biosynthesis pathways and the implementation of diverse ETC reactions across the microbial tree of life. We highlighted missing annotations that were essential to energy biosynthesis in our models. We examine the diversity of these pathways across all microbial life and enable the scientific community to explore the analyses generated from this large-scale analysis of over 8000 microbial genomes.

Electronic supplementary material

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

Yeast petites and small colony variants: for everything there is a season

The yeast petite mutant was first found in the yeast Saccharomyces cerevisiae. The colony is small because of a block in the aerobic respiratory chain pathway, which generates ATP. The petite yeasts are thus unable to grow on nonfermentable carbon sources (such as glycerol or ethanol), and form small anaerobic-sized colonies when grown in the presence of fermentable carbon sources (such as glucose). The petite phenotype results from mutations in the mitochondrial genome, loss of mitochondria, or mutations in the host cell genome. The latter mutations affect nuclear-encoded genes involved in oxidative phosphorylation and these mutants are termed neutral petites. They all produce wild-type progeny when crossed with a wild-type strain. The staphylococcal small colony variant (SCV) is a slow-growing mutant that typically exhibits the loss of many phenotypic characteristics and pathogenic traits. SCVs are mostly small, nonpigmented, and nonhaemolytic. Their small size is often due to an inability to synthesize electron transport chain components and so cannot generate ATP by oxidative phosphorylation. Evidence suggests that they are responsible for persistent and/or recurrent infections. This chapter compares the physiological and genetic basis of the petite mutants and SCVs. The review focuses principally on two representatives, the eukaryote S. cerevisiae and the prokaryote Staphylococcus aureus. There is, clearly, commonality in the physiological response. Interestingly, the similarity, based on their physiological states, has not been commented on previously. The finding of an overlapping physiological response that occurs across a taxonomic divide is novel.

Photorespiratory Bypasses Lead to Increased Growth in Arabidopsis thaliana: Are Predictions Consistent with Experimental Evidence?

Arguably, the biggest challenge of modern plant systems biology lies in predicting the performance of plant species, and crops in particular, upon different intracellular and external perturbations. Recently, an increased growth of Arabidopsis thaliana plants was achieved by introducing two different photorespiratory bypasses via metabolic engineering. Here, we investigate the extent to which these findings match the predictions from constraint-based modeling. To determine the effect of the employed metabolic network model on the predictions, we perform a comparative analysis involving three state-of-the-art metabolic reconstructions of A. thaliana. In addition, we investigate three scenarios with respect to experimental findings on the ratios of the carboxylation and oxygenation reactions of Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). We demonstrate that the condition-dependent growth phenotypes of one of the engineered bypasses can be qualitatively reproduced by each reconstruction, particularly upon considering the additional constraints with respect to the ratio of fluxes for the RuBisCO reactions. Moreover, our results lend support for the hypothesis of a reduced photorespiration in the engineered plants, and indicate that specific changes in CO₂ exchange as well as in the proxies for co-factor turnover are associated with the predicted growth increase in the engineered plants. We discuss our findings with respect to the structure of the used models, the modeling approaches taken, and the available experimental evidence. Our study sets the ground for investigating other strategies for increase of plant biomass by insertion of synthetic reactions.

Zooming-in on cancer metabolic rewiring with tissue specific constraint-based models

The metabolic rearrangements occurring in cancer cells can be effectively investigated with a Systems Biology approach supported by metabolic network modeling. We here present tissue-specific constraint-based core models for three different types of tumors (liver, breast and lung) that serve this purpose. The core models were extracted and manually curated from the corresponding genome-scale metabolic models in the Human Metabolic Atlas database with a focus on the pathways that are known to play a key role in cancer growth and proliferation. Along similar lines, we also reconstructed a core model from the original general human metabolic network to be used as a reference model. A comparative Flux Balance Analysis between the reference and the cancer models highlighted both a clear distinction between the two conditions and a heterogeneity within the three different cancer types in terms of metabolic flux distribution. These results emphasize the need for modeling approaches able to keep up with this tumoral heterogeneity in order to identify more suitable drug targets and develop effective treatments. According to this perspective, we identified key points able to reverse the tumoral phenotype toward the reference one or vice-versa.

Simulating Serial-Target Antibacterial Drug Synergies Using Flux Balance Analysis

Flux balance analysis (FBA) is an increasingly useful approach for modeling the behavior of metabolic systems. However, standard FBA modeling of genetic knockouts cannot predict drug combination synergies observed between serial metabolic targets, even though such synergies give rise to some of the most widely used antibiotic treatments. Here we extend FBA modeling to simulate responses to chemical inhibitors at varying concentrations, by diverting enzymatic flux to a waste reaction. This flux diversion yields very similar qualitative predictions to prior methods for single target activity. However, we find very different predictions for combinations, where flux diversion, which mimics the kinetics of competitive metabolic inhibitors, can explain serial target synergies between metabolic enzyme inhibitors that we confirmed in Escherichia coli cultures. FBA flux diversion opens the possibility for more accurate genome-scale predictions of drug synergies, which can be used to suggest treatments for infections and other diseases.

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.

Design Principles as a Guide for Constraint Based and Dynamic Modeling: Towards an Integrative Workflow

During the last 10 years, systems biology has matured from a fuzzy concept combining omics, mathematical modeling and computers into a scientific field on its own right. In spite of its incredible potential, the multilevel complexity of its objects of study makes it very difficult to establish a reliable connection between data and models. The great number of degrees of freedom often results in situations, where many different models can explain/fit all available datasets. This has resulted in a shift of paradigm from the initially dominant, maybe naive, idea of inferring the system out of a number of datasets to the application of different techniques that reduce the degrees of freedom before any data set is analyzed. There is a wide variety of techniques available, each of them can contribute a piece of the puzzle and include different kinds of experimental information. But the challenge that remains is their meaningful integration. Here we show some theoretical results that enable some of the main modeling approaches to be applied sequentially in a complementary manner, and how this workflow can benefit from evolutionary reasoning to keep the complexity of the problem in check. As a proof of concept, we show how the synergies between these modeling techniques can provide insight into some well studied problems: Ammonia assimilation in bacteria and an unbranched linear pathway with end-product inhibition.

Reconstruction and Use of Microbial Metabolic Networks: the Core Escherichia coli Metabolic Model as an Educational Guide

Biochemical network reconstructions have become popular tools in systems biology. Metabolicnetwork reconstructions are biochemically, genetically, and genomically (BiGG) structured databases of biochemical reactions and metabolites. They contain information such as exact reaction stoichiometry, reaction reversibility, and the relationships between genes, proteins, and reactions. Network reconstructions have been used extensively to study the phenotypic behavior of wild-type and mutant stains under a variety of conditions, linking genotypes with phenotypes. Such phenotypic simulations have allowed for the prediction of growth after genetic manipulations, prediction of growth phenotypes after adaptive evolution, and prediction of essential genes. Additionally, because network reconstructions are organism specific, they can be used to understand differences between organisms of species in a functional context.There are different types of reconstructions representing various types of biological networks (metabolic, regulatory, transcription/translation). This chapter serves as an introduction to metabolic and regulatory network reconstructions and models and gives a complete description of the core Escherichia coli metabolic model. This model can be analyzed in any computational format (such as MATLAB or Mathematica) based on the information given in this chapter. The core E. coli model is a small-scale model that can be used for educational purposes. It is meant to be used by senior undergraduate and first-year graduate students learning about constraint-based modeling and systems biology. This model has enough reactions and pathways to enable interesting and insightful calculations, but it is also simple enough that the results of such calculations can be understoodeasily.

Genome-scale modeling for metabolic engineering

We focus on the application of constraint-based methodologies and, more specifically, flux balance analysis in the field of metabolic engineering, and enumerate recent developments and successes of the field. We also review computational frameworks that have been developed with the express purpose of automatically selecting optimal gene deletions for achieving improved production of a chemical of interest. The application of flux balance analysis methods in rational metabolic engineering requires a metabolic network reconstruction and a corresponding in silico metabolic model for the microorganism in question. For this reason, we additionally present a brief overview of automated reconstruction techniques. Finally, we emphasize the importance of integrating metabolic networks with regulatory information-an area which we expect will become increasingly important for metabolic engineering-and present recent developments in the field of metabolic and regulatory integration.

Validation of an FBA model for Pichia pastoris in chemostat cultures

Background

Constraint-based metabolic models and flux balance analysis (FBA) have been extensively used in the last years to investigate the behavior of cells and also as basis for different industrial applications. In this context, this work provides a validation of a small-sized FBA model of the yeast Pichia pastoris. Our main objective is testing how accurate is the hypothesis of maximum growth to predict the behavior of P. pastoris in a range of experimental environments.

Results

A constraint-based model of P. pastoris was previously validated using metabolic flux analysis (MFA). In this paper we have verified the model ability to predict the cells behavior in different conditions without introducing measurements, experimental parameters, or any additional constraint, just by assuming that cells will make the best use of the available resources to maximize its growth. In particular, we have tested FBA model ability to: (a) predict growth yields over single substrates (glucose, glycerol, and methanol); (b) predict growth rate, substrate uptakes, respiration rates, and by-product formation in scenarios where different substrates are available (glucose, glycerol, methanol, or mixes of methanol and glycerol); (c) predict the different behaviors of P. pastoris cultures in aerobic and hypoxic conditions for each single substrate. In every case, experimental data from literature are used as validation.

Conclusions

We conclude that our predictions based on growth maximisation are reasonably accurate, but still far from perfect. The deviations are significant in scenarios where P. pastoris grows on methanol, suggesting that the hypothesis of maximum growth could be not dominating in these situations. However, predictions are much better when glycerol or glucose are used as substrates. In these scenarios, even if our FBA model is small and imposes a strong assumption regarding how cells will regulate their metabolic fluxes, it provides reasonably good predictions in terms of growth, substrate preference, product formation, and respiration rates.

Electronic supplementary material

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

Stoichiometric model and flux balance analysis for a mixed culture of Leptospirillum ferriphilum and Ferroplasma acidiphilum

The oxidation process of sulfide minerals in natural environments is achieved by microbial communities from the Archaea and Bacteria domains. A metabolic reconstruction of two dominant species, Leptospirillum ferriphilum and Ferroplasma acidiphilum, which are always found together as a mixed culture in this natural environments, was made. The metabolic model, composed of 152 internal reactions and 29 transport reactions, describes the main interactions between these species, assuming that both use ferrous iron as energy source, and F. acidiphilum takes advantage of the organic compounds secreted by L. ferriphilum for chemomixotrophic growth. A first metabolic model for a mixed culture used in bacterial leaching is proposed in this article, which pretends to represent the characteristics of the mixed culture in a simplified manner. It was evaluated with experimental data through flux balance analysis (FBA) using as objective function the maximization of biomass. The growth yields on ferrous iron obtained for each microorganism are consistent with experimental data, and the flux distribution obtained allows understanding of the metabolic capabilities of both microorganisms growing together in a bioleaching process. The model was used to simulate the growth of F. acidiphilum on different substrates, to determine in silico which compounds maximize cell growth, and which are essential. Knockout simulations were carried out for L. ferriphilum and F. acidiphilum metabolic models, predicting key enzymes of central metabolism. The results of this analysis are consistent with experimental data from literature, showing a robust behavior of the metabolic model.

Model-guided identification of gene deletion targets for metabolic engineering in Saccharomyces cerevisiae

Identification of metabolic engineering strategies for rerouting intracellular fluxes towards a desired product is often a challenging task owing to the topological and regulatory complexity of metabolic networks. Genome-scale metabolic models help tackling this complexity through systematic consideration of mass balance and reaction directionality constraints over the entire network. Here, we describe how genome-scale metabolic models can be used for identifying gene deletion targets leading to increased production of the desired product. Vanillin production in Saccharomyces cerevisiae is used as a case study throughout this chapter.