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

Discretised Flux Balance Analysis for Reaction-Diffusion Simulation of Single-Cell Metabolism

Metabolites have to diffuse within the sub-cellular compartments they occupy to specific locations where enzymes are, so reactions could occur. Conventional flux balance analysis (FBA), a method based on linear programming that is commonly used to model metabolism, implicitly assumes that all enzymatic reactions are not diffusion-limited though that may not always be the case. In this work, we have developed a spatial method that implements FBA on a grid-based system, to enable the exploration of diffusion effects on metabolism. Specifically, the method discretises a living cell into a two-dimensional grid, represents the metabolic reactions in each grid element as well as the diffusion of metabolites to and from neighbouring elements, and simulates the system as a single linear programming problem. We varied the number of rows and columns in the grid to simulate different cell shapes, and the method was able to capture diffusion effects at different shapes. We then used the method to simulate heterogeneous enzyme distribution, which suggested a theoretical effect on variability at the population level. We propose the use of this method, and its future extensions, to explore how spatiotemporal organisation of sub-cellular compartments and the molecules within could affect cell behaviour.

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

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

Revealing the Metabolic Alterations during Biofilm Development of Burkholderia cenocepacia Based on Genome-Scale Metabolic Modeling

Burkholderia cenocepacia is among the important pathogens isolated from cystic fibrosis (CF) patients. It has attracted considerable attention because of its capacity to evade host immune defenses during chronic infection. Advances in systems biology methodologies have led to the emergence of methods that integrate experimental transcriptomics data and genome-scale metabolic models (GEMs). Here, we integrated transcriptomics data of bacterial cells grown on exponential and biofilm conditions into a manually curated GEM of B. cenocepacia. We observed substantial differences in pathway response to different growth conditions and alternative pathway susceptibility to extracellular nutrient availability. For instance, we found that blockage of the reactions was vital through the lipid biosynthesis pathways in the exponential phase and the absence of microenvironmental lysine and tryptophan are essential for survival. During biofilm development, bacteria mostly had conserved lipid metabolism but altered pathway activities associated with several amino acids and pentose phosphate pathways. Furthermore, conversion of serine to pyruvate and 2,5-dioxopentanoate synthesis are also identified as potential targets for metabolic remodeling during biofilm development. Altogether, our integrative systems biology analysis revealed the interactions between the bacteria and its microenvironment and enabled the discovery of antimicrobial targets for biofilm-related diseases.

Spatiotemporal metabolic modeling of bacterial life in complex habitats

The combination of genome-scale metabolic networks with spatially explicit representation of microbial habitats (spatiotemporal metabolic network modeling) paves the way to predict complex metabolic landscapes to a hitherto unparalleled detail, thus providing new insights into trophic interactions occurring at different scales. Placing detailed bacterial metabolism in realistic physical environment highlights the roles of physical barriers and diffusional bottlenecks on bacterial community interactions, structure and stability. We review recent advances in spatiotemporal metabolic network modeling using a few illustrative examples that highlight the immense potential of these novel approaches to interpret and design metabolic mediated interactions in structures (natural and engineered) environments.

Automated engineering of synthetic metabolic pathways for efficient biomanufacturing

Metabolic engineering involves the engineering and optimization of processes from single-cell to fermentation in order to increase production of valuable chemicals for health, food, energy, materials and others. A systems approach to metabolic engineering has gained traction in recent years thanks to advances in strain engineering, leading to an accelerated scaling from rapid prototyping to industrial production. Metabolic engineering is nowadays on track towards a truly manufacturing technology, with reduced times from conception to production enabled by automated protocols for DNA assembly of metabolic pathways in engineered producer strains. In this review, we discuss how the success of the metabolic engineering pipeline often relies on retrobiosynthetic protocols able to identify promising production routes and dynamic regulation strategies through automated biodesign algorithms, which are subsequently assembled as embedded integrated genetic circuits in the host strain. Those approaches are orchestrated by an experimental design strategy that provides optimal scheduling planning of the DNA assembly, rapid prototyping and, ultimately, brings forward an accelerated Design-Build-Test-Learn cycle and the overall optimization of the biomanufacturing process. Achieving such a vision will address the increasingly compelling demand in our society for delivering valuable biomolecules in an affordable, inclusive and sustainable bioeconomy.

Artificial consortium demonstrates emergent properties of enhanced cellulosic-sugar degradation and biofuel synthesis

Planktonic cultures, of a rationally designed consortium, demonstrated emergent properties that exceeded the sums of monoculture properties, including a >200% increase in cellobiose catabolism, a >100% increase in glycerol catabolism, a >800% increase in ethanol production, and a >120% increase in biomass productivity. The consortium was designed to have a primary and secondary-resource specialist that used crossfeeding with a positive feedback mechanism, division of labor, and nutrient and energy transfer via necromass catabolism. The primary resource specialist was Clostridium phytofermentans (a.k.a. Lachnoclostridium phytofermentans), a cellulolytic, obligate anaerobe. The secondary-resource specialist was Escherichia coli, a versatile, facultative anaerobe, which can ferment glycerol and byproducts of cellobiose catabolism. The consortium also demonstrated emergent properties of enhanced biomass accumulation when grown as biofilms, which created high cell density communities with gradients of species along the vertical axis. Consortium biofilms were robust to oxic perturbations with E. coli consuming O₂, creating an anoxic environment for C. phytofermentans. Anoxic/oxic cycling further enhanced biomass productivity of the biofilm consortium, increasing biomass accumulation ~250% over the sum of the monoculture biofilms. Consortium emergent properties were credited to several synergistic mechanisms. E. coli consumed inhibitory byproducts from cellobiose catabolism, driving higher C. phytofermentans growth and higher cellulolytic enzyme production, which in turn provided more substrate for E. coli. E. coli necromass enhanced C. phytofermentans growth while C. phytofermentans necromass aided E. coli growth via the release of peptides and amino acids, respectively. In aggregate, temporal cycling of necromass constituents increased flux of cellulose-derived resources through the consortium. The study establishes a consortia-based, bioprocessing strategy built on naturally occurring interactions for improved conversion of cellulose-derived sugars into bioproducts.

Genome-scale reconstruction of Paenarthrobacter aurescens TC1 metabolic model towards the study of atrazine bioremediation

Atrazine is an herbicide and a pollutant of great environmental concern that is naturally biodegraded by microbial communities. Paenarthrobacter aurescens TC1 is one of the most studied degraders of this herbicide. Here, we developed a genome scale metabolic model for P. aurescens TC1, iRZ1179, to study the atrazine degradation process at organism level. Constraint based flux balance analysis and time dependent simulations were used to explore the organism’s phenotypic landscape. Simulations aimed at designing media optimized for supporting growth and enhancing degradation, by passing the need in strain design via genetic modifications. Growth and degradation simulations were carried with more than 100 compounds consumed by P. aurescens TC1. In vitro validation confirmed the predicted classification of different compounds as efficient, moderate or poor stimulators of growth. Simulations successfully captured previous reports on the use of glucose and phosphate as bio-stimulators of atrazine degradation, supported by in vitro validation. Model predictions can go beyond supplementing the medium with a single compound and can predict the growth outcomes for higher complexity combinations. Hence, the analysis demonstrates that the exhaustive power of the genome scale metabolic reconstruction allows capturing complexities that are beyond common biochemical expertise and knowledge and further support the importance of computational platforms for the educated design of complex media. The model presented here can potentially serve as a predictive tool towards achieving optimal biodegradation efficiencies and for the development of ecologically friendly solutions for pollutant degradation.

In Silico Metabolic Design of Two-Strain Biofilm Systems Predicts Enhanced Biomass Production and Biochemical Synthesis

Engineered biofilm consortia have the potential to solve important biotechnological problems that have proved difficult for monoculture biofilms and planktonic consortia, such as conversion of lignocellulosic material to useful biochemicals. While considerable experimental progress has been reported for engineering and characterizing biofilm consortia, the field still lacks in silico tools for simulation, design, and optimization of stable, robust, and productive designed consortia. We developed biofilm consortia metabolic models for two coculture systems centered around the ecological design motif of a primary cell type that utilizes a supplied electron donor and secretes acetate as a byproduct and a secondary cell type that consumes the acetate, relieving byproduct inhibition on the primary cell type and enhancing overall system biomass. The models presented in this paper predict that distinct metabolic niches for the two cell types could be established by supplying electron donors and acceptors at opposite ends of the biofilm and that acetate consumption by the secondary cell type could increase total biomass accumulation and the synthesis of valuable biochemicals, such as isobutanol, by the primary cell type. System tunability is enhanced when each cell type is supplied with a unique terminal electron acceptor at opposite ends of the biofilm rather than competing for a common electron acceptor. Our model provides good qualitative agreement with data for a synthetic Escherichia coli coculture system, suggesting that the proposed design rules may have wide applicability to engineered biofilm consortia.

Engineering Robust Production Microbes for Large-Scale Cultivation

Systems biology and synthetic biology are increasingly used to examine and modulate complex biological systems. As such, many issues arising during scaling-up microbial production processes can be addressed using these approaches. We review differences between laboratory-scale cultures and larger-scale processes to provide a perspective on those strain characteristics that are especially important during scaling. Systems biology has been used to examine a range of microbial systems for their response in bioreactors to fluctuations in nutrients, dissolved gases, and other stresses. Synthetic biology has been used both to assess and modulate strain response, and to engineer strains to improve production. We discuss these approaches and tools in the context of their use in engineering robust microbes for applications in large-scale production.

Competitive resource allocation to metabolic pathways contributes to overflow metabolisms and emergent properties in cross-feeding microbial consortia

Resource scarcity is a common stress in nature and has a major impact on microbial physiology. This review highlights microbial acclimations to resource scarcity, focusing on resource investment strategies for chemoheterotrophs from the molecular level to the pathway level. Competitive resource allocation strategies often lead to a phenotype known as overflow metabolism; the resulting overflow byproducts can stabilize cooperative interactions in microbial communities and can lead to cross-feeding consortia. These consortia can exhibit emergent properties such as enhanced resource usage and biomass productivity. The literature distilled here draws parallels between in silico and laboratory studies and ties them together with ecological theories to better understand microbial stress responses and mutualistic consortia functioning.

The best models of metabolism

Biochemical systems are among of the oldest application areas of mathematical modeling. Spanning a time period of over one hundred years, the repertoire of options for structuring a model and for formulating reactions has been constantly growing, and yet, it is still unclear whether or to what degree some models are better than others and how the modeler is to choose among them. In fact, the variety of options has become overwhelming and difficult to maneuver for novices and experts alike. This review outlines the metabolic model design process and discusses the numerous choices for modeling frameworks and mathematical representations. It tries to be inclusive, even though it cannot be complete, and introduces the various modeling options in a manner that is as unbiased as that is feasible. However, the review does end with personal recommendations for the choices of default models. WIREs Syst Biol Med 2017, 9:e1391. doi: 10.1002/wsbm.1391 For further resources related to this article, please visit the WIREs website.

Cellular trade-offs and optimal resource allocation during cyanobacterial diurnal growth

Cyanobacteria are an integral part of Earth's biogeochemical cycles and a promising resource for the synthesis of renewable bioproducts from atmospheric CO₂ Growth and metabolism of cyanobacteria are inherently tied to the diurnal rhythm of light availability. As yet, however, insight into the stoichiometric and energetic constraints of cyanobacterial diurnal growth is limited. Here, we develop a computational framework to investigate the optimal allocation of cellular resources during diurnal phototrophic growth using a genome-scale metabolic reconstruction of the cyanobacterium Synechococcus elongatus PCC 7942. We formulate phototrophic growth as an autocatalytic process and solve the resulting time-dependent resource allocation problem using constraint-based analysis. Based on a narrow and well-defined set of parameters, our approach results in an ab initio prediction of growth properties over a full diurnal cycle. The computational model allows us to study the optimality of metabolite partitioning during diurnal growth. The cyclic pattern of glycogen accumulation, an emergent property of the model, has timing characteristics that are in qualitative agreement with experimental findings. The approach presented here provides insight into the time-dependent resource allocation problem of phototrophic diurnal growth and may serve as a general framework to assess the optimality of metabolic strategies that evolved in phototrophic organisms under diurnal conditions.

Diet shifts provoke complex and variable changes in the metabolic networks of the ruminal microbiome

BACKGROUND

Grazing mammals rely on their ruminal microbial symbionts to convert plant structural biomass into metabolites they can assimilate. To explore how this complex metabolic system adapts to the host animal's diet, we inferred a microbiome-level metabolic network from shotgun metagenomic data.

RESULTS

Using comparative genomics, we then linked this microbial network to that of the host animal using a set of interface metabolites likely to be transferred to the host. When the host sheep were fed a grain-based diet, the induced microbial metabolic network showed several critical differences from those seen on the evolved forage-based diet. Grain-based (e.g., concentrate) diets tend to be dominated by a smaller set of reactions that employ metabolites that are nearer in network space to the host's metabolism. In addition, these reactions are more central in the network and employ substrates with shorter carbon backbones. Despite this apparent lower complexity, the concentrate-associated metabolic networks are actually more dissimilar from each other than are those of forage-fed animals. Because both groups of animals were initially fed on a forage diet, we propose that the diet switch drove the appearance of a number of different microbial networks, including a degenerate network characterized by an inefficient use of dietary nutrients. We used network simulations to show that such disparate networks are not an unexpected result of a diet shift.

CONCLUSION

We argue that network approaches, particularly those that link the microbial network with that of the host, illuminate aspects of the structure of the microbiome not seen from a strictly taxonomic perspective. In particular, different diets induce predictable and significant differences in the enzymes used by the microbiome. Nonetheless, there are clearly a number of microbiomes of differing structure that show similar functional properties. Changes such as a diet shift uncover more of this type of diversity.

Toward Multiscale Models of Cyanobacterial Growth: A Modular Approach

Oxygenic photosynthesis dominates global primary productivity ever since its evolution more than three billion years ago. While many aspects of phototrophic growth are well understood, it remains a considerable challenge to elucidate the manifold dependencies and interconnections between the diverse cellular processes that together facilitate the synthesis of new cells. Phototrophic growth involves the coordinated action of several layers of cellular functioning, ranging from the photosynthetic light reactions and the electron transport chain, to carbon-concentrating mechanisms and the assimilation of inorganic carbon. It requires the synthesis of new building blocks by cellular metabolism, protection against excessive light, as well as diurnal regulation by a circadian clock and the orchestration of gene expression and cell division. Computational modeling allows us to quantitatively describe these cellular functions and processes relevant for phototrophic growth. As yet, however, computational models are mostly confined to the inner workings of individual cellular processes, rather than describing the manifold interactions between them in the context of a living cell. Using cyanobacteria as model organisms, this contribution seeks to summarize existing computational models that are relevant to describe phototrophic growth and seeks to outline their interactions and dependencies. Our ultimate aim is to understand cellular functioning and growth as the outcome of a coordinated operation of diverse yet interconnected cellular processes.

Metabolic modeling of a chronic wound biofilm consortium predicts spatial partitioning of bacterial species

Background

Chronic wounds are often colonized by consortia comprised of different bacterial species growing as biofilms on a complex mixture of wound exudate. Bacteria growing in biofilms exhibit phenotypes distinct from planktonic growth, often rendering the application of antibacterial compounds ineffective. Computational modeling represents a complementary tool to experimentation for generating fundamental knowledge and developing more effective treatment strategies for chronic wound biofilm consortia.

Results

We developed spatiotemporal models to investigate the multispecies metabolism of a biofilm consortium comprised of two common chronic wound isolates: the aerobe Pseudomonas aeruginosa and the facultative anaerobe Staphylococcus aureus. By combining genome-scale metabolic reconstructions with partial differential equations for metabolite diffusion, the models were able to provide both temporal and spatial predictions with genome-scale resolution. The models were used to analyze the metabolic differences between single species and two species biofilms and to demonstrate the tendency of the two bacteria to spatially partition in the multispecies biofilm as observed experimentally. Nutrient gradients imposed by supplying glucose at the bottom and oxygen at the top of the biofilm induced spatial partitioning of the two species, with S. aureus most concentrated in the anaerobic region and P. aeruginosa present only in the aerobic region. The two species system was predicted to support a maximum biofilm thickness much greater than P. aeruginosa alone but slightly less than S. aureus alone, suggesting an antagonistic metabolic effect of P. aeruginosa on S. aureus. When each species was allowed to enhance its growth through consumption of secreted metabolic byproducts assuming identical uptake kinetics, the competitiveness of P. aeruginosa was further reduced due primarily to the more efficient lactate metabolism of S. aureus. Lysis of S. aureus by a small molecule inhibitor secreted from P. aeruginosa and/or P. aeruginosa aerotaxis were predicted to substantially increase P. aeruginosa competitiveness in the aerobic region, consistent with in vitro experimental studies.

Conclusions

Our biofilm modeling approach allows the prediction of individual species metabolism and interspecies interactions in both time and space with genome-scale resolution. This study yielded new insights into the multispecies metabolism of a chronic wound biofilm, in particular metabolic factors that may lead to spatial partitioning of the two bacterial species. We believe that P. aeruginosa lysis of S. aureus combined with nutrient competition is a particularly relevant scenario for which model predictions could be tested experimentally.

Electronic supplementary material

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