Designing microbial communities to maximize the thermodynamic driving force for the production of chemicals

DyMMM-LEAPS: An ML-based framework for modulating evenness and stability in synthetic microbial communities

There have been a growing number of computational strategies to aid in the design of synthetic microbial consortia. A framework to identify regions in parametric space to maximize two essential properties, evenness and stability, is critical. In this study, we introduce DyMMM-LEAPS (dynamic multispecies metabolic modeling-locating evenness and stability in large parametric space), an extension of the DyMMM framework. Our method explores the large parametric space of genetic circuits in synthetic microbial communities to identify regions of evenness and stability. Due to the high computational costs of exhaustive sampling, we utilize adaptive sampling and surrogate modeling to reduce the number of simulations required to map the vast space. Our framework predicts engineering targets and computes their operating ranges to maximize the probability of the engineered community to have high evenness and stability. We demonstrate our approach by simulating five cocultures and one three-strain culture with different social interactions (cooperation, competition, and predation) employing quorum-sensing-based genetic circuits. In addition to guiding circuit tuning, our pipeline gives an opportunity for a detailed analysis of pockets of evenness and stability for the circuit under investigation, which can further help dissect the relationship between the two properties. DyMMM-LEAPS is easily customizable and can be expanded to a larger community with more complex interactions.

PyCoMo: a python package for community metabolic model creation and analysis

SUMMARY

PyCoMo is a python package for quick and easy generation of genome-scale compartmentalized community metabolic models that are compliant with current openCOBRA file formats. The resulting models can be used to predict (i) the maximum growth rate at a given abundance profile, (ii) the feasible community compositions at a given growth rate, and (iii) all exchange metabolites and cross-feeding interactions in a community metabolic model independent of the abundance profile; we demonstrate PyCoMo's capability by analysing methane production in a previously published simplified biogas community metabolic model.

AVAILABILITY AND IMPLEMENTATION

PyCoMo is freely available under an MIT licence at http://github.com/univieCUBE/PyCoMo, the Python Package Index, and Zenodo.

Creating new-to-nature carbon fixation: A guide

Synthetic biology aims at designing new biological functions from first principles. These new designs allow to expand the natural solution space and overcome the limitations of naturally evolved systems. One example is synthetic CO2-fixation pathways that promise to provide more efficient ways for the capture and conversion of CO2 than natural pathways, such as the Calvin Benson Bassham (CBB) cycle of photosynthesis. In this review, we provide a practical guideline for the design and realization of such new-to-nature CO2-fixation pathways. We introduce the concept of "synthetic CO2-fixation", and give a general overview over the enzymology and topology of synthetic pathways, before we derive general principles for their design from their eight naturally evolved analogs. We provide a comprehensive summary of synthetic carbon-assimilation pathways and derive a step-by-step, practical guide from the theoretical design to their practical implementation, before ending with an outlook on new developments in the field.

Model-driven insights into the effects of temperature on metabolism

Temperature affects cellular processes at different spatiotemporal scales, and identifying the genetic and molecular mechanisms underlying temperature responses paves the way to develop approaches for mitigating the effects of future climate scenarios. A systems view of the effects of temperature on cellular physiology can be obtained by focusing on metabolism since: (i) its functions depend on transcription and translation and (ii) its outcomes support organisms' development, growth, and reproduction. Here we provide a systematic review of modelling efforts directed at investigating temperature effects on properties of single biochemical reactions, system-level traits, metabolic subsystems, and whole-cell metabolism across different prokaryotes and eukaryotes. We compare and contrast computational approaches and theories that facilitate modelling of temperature effects on key properties of enzymes and their consideration in constraint-based as well as kinetic models of metabolism. In addition, we provide a summary of insights from computational approaches, facilitating integration of omics data from temperature-modulated experiments with models of metabolic networks, and review the resulting biotechnological applications. Lastly, we provide a perspective on how different types of metabolic modelling can profit from developments in machine learning and models of different cellular layers to improve model-driven insights into the effects of temperature relevant for biotechnological applications.

Network-wide thermodynamic constraints shape NAD(P)H cofactor specificity of biochemical reactions

The ubiquitous coexistence of the redox cofactors NADH and NADPH is widely considered to facilitate an efficient operation of cellular redox metabolism. However, it remains unclear what shapes the NAD(P)H specificity of specific redox reactions. Here, we present a computational framework to analyze the effect of redox cofactor swaps on the maximal thermodynamic potential of a metabolic network and use it to investigate key aspects of redox cofactor redundancy in Escherichia coli. As one major result, our analysis suggests that evolved NAD(P)H specificities are largely shaped by metabolic network structure and associated thermodynamic constraints enabling thermodynamic driving forces that are close or even identical to the theoretical optimum and significantly higher compared to random specificities. Furthermore, while redundancy of NAD(P)H is clearly beneficial for thermodynamic driving forces, a third redox cofactor would require a low standard redox potential to be advantageous. Our approach also predicts trends of redox-cofactor concentration ratios and could facilitate the design of optimal redox cofactor specificities.

A structured evaluation of genome-scale constraint-based modeling tools for microbial consortia

Harnessing the power of microbial consortia is integral to a diverse range of sectors, from healthcare to biotechnology to environmental remediation. To fully realize this potential, it is critical to understand the mechanisms behind the interactions that structure microbial consortia and determine their functions. Constraint-based reconstruction and analysis (COBRA) approaches, employing genome-scale metabolic models (GEMs), have emerged as the state-of-the-art tool to simulate the behavior of microbial communities from their constituent genomes. In the last decade, many tools have been developed that use COBRA approaches to simulate multi-species consortia, under either steady-state, dynamic, or spatiotemporally varying scenarios. Yet, these tools have not been systematically evaluated regarding their software quality, most suitable application, and predictive power. Hence, it is uncertain which tools users should apply to their system and what are the most urgent directions that developers should take in the future to improve existing capacities. This study conducted a systematic evaluation of COBRA-based tools for microbial communities using datasets from two-member communities as test cases. First, we performed a qualitative assessment in which we evaluated 24 published tools based on a list of FAIR (Findability, Accessibility, Interoperability, and Reusability) features essential for software quality. Next, we quantitatively tested the predictions in a subset of 14 of these tools against experimental data from three different case studies: a) syngas fermentation by C. autoethanogenum and C. kluyveri for the static tools, b) glucose/xylose fermentation with engineered E. coli and S. cerevisiae for the dynamic tools, and c) a Petri dish of E. coli and S. enterica for tools incorporating spatiotemporal variation. Our results show varying performance levels of the best qualitatively assessed tools when examining the different categories of tools. The differences in the mathematical formulation of the approaches and their relation to the results were also discussed. Ultimately, we provide recommendations for refining future GEM microbial modeling tools.

Optimal evaluation of energy yield and driving force in microbial metabolic pathway variants

This work presents a methodology to evaluate the bioenergetic feasibility of alternative metabolic pathways for a given microbial conversion, optimising their energy yield and driving forces as a function of the concentration of metabolic intermediates. The tool, based on thermodynamic principles and multi-objective optimisation, accounts for pathway variants in terms of different electron carriers, as well as energy conservation (proton translocating) reactions within the pathway. The method also accommodates other constraints, some of them non-linear, such as the balance of conserved moieties. The approach involves the transformation of the maximum energy yield problem into a multi-objective mixed-integer linear optimisation problem which is then subsequently solved using the epsilon-constraint method, highlighting the trade-off between yield and rate in metabolic reactions. The methodology is applied to analyse several pathway alternatives occurring during propionate oxidation in anaerobic fermentation processes, as well as to the reverse TCA cycle pathway occurring during autotrophic microbial CO2 fixation. The results obtained using the developed methodology match previously reported literature and bring about insights into the studied pathways.

Genome-scale community modeling for deciphering the inter-microbial metabolic interactions in fungus-farming termite gut microbiome

Specialized microbial communities in the fungus-farming termite gut and fungal comb microbiome help maintain host nutrition through interactive biochemical activities of complex carbohydrate degradation. Numerous research studies have been focused on identifying the microbial species in the termite gut and fungal comb microbiota, but the community-wide metabolic interaction patterns remain obscure. The inter-microbial metabolic interactions in the community environment are essential for executing biochemical processes like complex carbohydrate degradation and maintaining the host's physicochemical homeostasis. Recent progress in high-throughput sequencing techniques and mathematical modeling provides suitable platforms for constructing multispecies genome-scale community metabolic models that can render sound knowledge about microbial metabolic interaction patterns. Here, we have implemented the genome-scale metabolic modeling strategy to map the relationship between genes, proteins, and reactions of 12 key bacterial species from fungal cultivating termite gut and fungal comb microbiota. The resulting individual genome-scale metabolic models (GEMs) have been analyzed using flux balance analysis (FBA) to optimize the metabolic flux distribution pattern. Further, these individual GEMs have been integrated into genome-scale community metabolic models where a heuristics-based computational procedure has been employed to track the inter-microbial metabolic interactions. Two separate genome-scale community metabolic models were reconstructed for the O. badius gut and fungal comb microbiome. Analysis of the community models showed up to ∼167% increased flux range in lignocellulose degradation, amino acid biosynthesis, and nucleotide metabolism pathways. The inter-microbial metabolic exchange of amino acids, SCFAs, and small sugars was also upregulated in the multispecies community for maximum biomass formation. The flux variability analysis (FVA) has also been performed to calculate the feasible flux range of metabolic reactions. Furthermore, based on the calculated metabolic flux values, newly defined parameters, i.e., pairwise metabolic assistance (PMA) and community metabolic assistance (CMA) showed that the microbial species are getting up to 15% higher metabolic benefits in the multispecies community compared to pairwise growth. Assessment of the inter-microbial metabolic interaction patterns through pairwise growth support index (PGSI) indicated an increased mutualistic interaction in the termite gut environment compared to the fungal comb. Thus, this genome-scale community modeling study provides a systematic methodology to understand the inter-microbial interaction patterns with several newly defined parameters like PMA, CMA, and PGSI. The microbial metabolic assistance and interaction patterns derived from this computational approach will enhance the understanding of combinatorial microbial activities and may help develop effective synergistic microcosms to utilize complex plant polymers.

Artificial microbial consortia for bioproduction processes

The application of artificial microbial consortia for biotechnological production processes is an emerging field in research as it offers great potential for the improvement of established as well as the development of novel processes. In this review, we summarize recent highlights in the usage of various microbial consortia for the production of, for example, platform chemicals, biofuels, or pharmaceutical compounds. It aims to demonstrate the great potential of co-cultures by employing different organisms and interaction mechanisms and exploiting their respective advantages. Bacteria and yeasts often offer a broad spectrum of possible products, fungi enable the utilization of complex lignocellulosic substrates via enzyme secretion and hydrolysis, and microalgae can feature their abilities to fixate CO2 through photosynthesis for other organisms as well as to form lipids as potential fuelstocks. However, the complexity of interactions between microbes require methods for observing population dynamics within the process and modern approaches such as modeling or automation for process development. After shortly discussing these interaction mechanisms, we aim to present a broad variety of successfully established co-culture processes to display the potential of artificial microbial consortia for the production of biotechnological products.

Toward merging bottom-up and top-down model-based designing of synthetic microbial communities

The increasing interest of microbial communities as promising biocatalyst is leading an intense effort into the development of computational frameworks assisting the analysis and rational engineering of such complex ecosystems. Here, we critically review the recent computational and model-guided advances in the system-level engineering of microbiome, including both the rational bottom-up and the evolutionary top-down approaches. Furthermore, we highlight modeling and computational methods supporting both engineering paradigms. Finally, we discuss the advantages of combining both strategies into a hybrid top-down/bottom-up (middle-out) strategy to engineer synthetic microbial communities with improved performance and scope.