Predicting Nash equilibria for microbial metabolic interactions

Deciphering and designing microbial communities by genome-scale metabolic modelling

Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.

Leveraging genome-scale metabolic models to understand aerobic methanotrophs

Genome-scale metabolic models (GEMs) are valuable tools serving systems biology and metabolic engineering. However, GEMs are still an underestimated tool in informing microbial ecology. Since their first application for aerobic gammaproteobacterial methane oxidizers less than a decade ago, GEMs have substantially increased our understanding of the metabolism of methanotrophs, a microbial guild of high relevance for the natural and biotechnological mitigation of methane efflux to the atmosphere. Particularly, GEMs helped to elucidate critical metabolic and regulatory pathways of several methanotrophic strains, predicted microbial responses to environmental perturbations, and were used to model metabolic interactions in cocultures. Here, we conducted a systematic review of GEMs exploring aerobic methanotrophy, summarizing recent advances, pointing out weaknesses, and drawing out probable future uses of GEMs to improve our understanding of the ecology of methane oxidizers. We also focus on their potential to unravel causes and consequences when studying interactions of methane-oxidizing bacteria with other methanotrophs or members of microbial communities in general. This review aims to bridge the gap between applied sciences and microbial ecology research on methane oxidizers as model organisms and to provide an outlook for future studies.

Microbial interactions from a new perspective: reinforcement learning reveals new insights into microbiome evolution

MOTIVATION

Microbes are essential part of all ecosystems, influencing material flow and shaping their surroundings. Metabolic modeling has been a useful tool and provided tremendous insights into microbial community metabolism. However, current methods based on flux balance analysis (FBA) usually fail to predict metabolic and regulatory strategies that lead to long-term survival and stability especially in heterogenous communities.

RESULTS

Here, we introduce a novel reinforcement learning algorithm, Self-Playing Microbes in Dynamic FBA, which treats microbial metabolism as a decision-making process, allowing individual microbial agents to evolve by learning and adapting metabolic strategies for enhanced long-term fitness. This algorithm predicts what microbial flux regulation policies will stabilize in the dynamic ecosystem of interest in the presence of other microbes with minimal reliance on predefined strategies. Throughout this article, we present several scenarios wherein our algorithm outperforms existing methods in reproducing outcomes, and we explore the biological significance of these predictions.

AVAILABILITY AND IMPLEMENTATION

The source code for this article is available at: https://github.com/chan-csu/SPAM-DFBA.

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.

Assumptions on decision making and environment can yield multiple steady states in microbial community models

BACKGROUND

Microbial community simulations using genome scale metabolic networks (GSMs) are relevant for many application areas, such as the analysis of the human microbiome. Such simulations rely on assumptions about the culturing environment, affecting if the culture may reach a metabolically stationary state with constant microbial concentrations. They also require assumptions on decision making by the microbes: metabolic strategies can be in the interest of individual community members or of the whole community. However, the impact of such common assumptions on community simulation results has not been investigated systematically.

RESULTS

Here, we investigate four combinations of assumptions, elucidate how they are applied in literature, provide novel mathematical formulations for their simulation, and show how the resulting predictions differ qualitatively. Our results stress that different assumption combinations give qualitatively different predictions on microbial coexistence by differential substrate utilization. This fundamental mechanism is critically under explored in the steady state GSM literature with its strong focus on coexistence states due to crossfeeding (division of labor). Furthermore, investigating a realistic synthetic community, where the two involved strains exhibit no growth in isolation, but grow as a community, we predict multiple modes of cooperation, even without an explicit cooperation mechanism.

CONCLUSIONS

Steady state GSM modelling of microbial communities relies both on assumed decision making principles and environmental assumptions. In principle, dynamic flux balance analysis addresses both. In practice, our methods that address the steady state directly may be preferable, especially if the community is expected to display multiple steady states.

Two-species community design of lactic acid bacteria for optimal production of lactate

Microbial communities that metabolise pentose and hexose sugars are useful in producing high-value chemicals, resulting in the effective conversion of raw materials to the product, a reduction in the production cost, and increased yield. Here, we present a computational analysis approach called CAMP (Co-culture/Community Analyses for Metabolite Production) that simulates and identifies appropriate communities to produce a metabolite of interest. To demonstrate this approach, we focus on the optimal production of lactate from various Lactic Acid Bacteria. We used genome-scale metabolic models (GSMMs) belonging to Lactobacillus, Leuconostoc, and Pediococcus species from the Virtual Metabolic Human (VMH; https://vmh.life/) resource and well-curated GSMMs of L. plantarum WCSF1 and L. reuteri JCM 1112. We analysed 1176 two-species communities using a constraint-based modelling method for steady-state flux-balance analysis of communities. Flux variability analysis was used to detect the maximum lactate flux in the communities. Using glucose or xylose as substrates separately or in combination resulted in either parasitism, amensalism, or mutualism being the dominant interaction behaviour in the communities. Interaction behaviour between members of the community was deduced based on variations in the predicted growth rates of monocultures and co-cultures. Acetaldehyde, ethanol, acetate, among other metabolites, were found to be cross-fed between community members. L. plantarum WCSF1 was found to be a member of communities with high lactate yields. In silico community optimisation strategies to predict reaction knock-outs for improving lactate flux were implemented. Reaction knock-outs of acetate kinase, phosphate acetyltransferase, and fumarate reductase in the communities were found to enhance lactate production.

Predicting Microbiome Metabolism and Interactions through Integrating Multidisciplinary Principles

In this Commentary, we will discuss some of the current trends and challenges in modeling microbiome metabolism. A focus will be the state of the art in the integration of metabolic networks, ecological and evolutionary principles, and spatiotemporal considerations, followed by envisioning integrated frameworks incorporating different principles and data to generate predictive models in the future.

Using Cartesian Doubt To Build a Sequencing-Based View of Microbiology

The technological leap of DNA sequencing generated a tension between modern metagenomics and historical microbiology. We are forcibly harmonizing the output of a modern tool with centuries of experimental knowledge derived from culture-based microbiology. As a thought experiment, we borrow the notion of Cartesian doubt from philosopher Rene Descartes, who used doubt to build a philosophical framework from his incorrigible statement that "I think therefore I am." We aim to cast away preconceived notions and conceptualize microorganisms through the lens of metagenomic sequencing alone. Specifically, we propose funding and building analysis and engineering methods that neither search for nor rely on the assumption of independent genomes bound by lipid barriers containing discrete functional roles and taxonomies. We propose that a view of microbial communities based in sequencing will engender novel insights into metagenomic structure and may capture functional biology not reflected within the current paradigm.

Nash equilibria in human sensorimotor interactions explained by Q-learning with intrinsic costs

The Nash equilibrium concept has previously been shown to be an important tool to understand human sensorimotor interactions, where different actors vie for minimizing their respective effort while engaging in a multi-agent motor task. However, it is not clear how such equilibria are reached. Here, we compare different reinforcement learning models to human behavior engaged in sensorimotor interactions with haptic feedback based on three classic games, including the prisoner's dilemma, and the symmetric and asymmetric matching pennies games. We find that a discrete analysis that reduces the continuous sensorimotor interaction to binary choices as in classical matrix games does not allow to distinguish between the different learning algorithms, but that a more detailed continuous analysis with continuous formulations of the learning algorithms and the game-theoretic solutions affords different predictions. In particular, we find that Q-learning with intrinsic costs that disfavor deviations from average behavior explains the observed data best, even though all learning algorithms equally converge to admissible Nash equilibrium solutions. We therefore conclude that it is important to study different learning algorithms for understanding sensorimotor interactions, as such behavior cannot be inferred from a game-theoretic analysis alone, that simply focuses on the Nash equilibrium concept, as different learning algorithms impose preferences on the set of possible equilibrium solutions due to the inherent learning dynamics.

Advances in Genome-Scale Metabolic Modeling toward Microbial Community Analysis of the Human Microbiome

A genome-scale metabolic model (GEM) represents metabolic pathways of an organism in a mathematical form and can be built using biochemistry and genome annotation data. GEMs are invaluable for understanding organisms since they analyze the metabolic capabilities and behaviors quantitatively and can predict phenotypes. The development of high-throughput data collection techniques led to an immense increase in omics data such as metagenomics, which expand our knowledge on the human microbiome, but this also created a need for systematic analysis of these data. In recent years, GEMs have also been reconstructed for microbial species, including human gut microbiota, and methods for the analysis of microbial communities have been developed to examine the interaction between the organisms or the host. The purpose of this review is to provide a comprehensive guide for the applications of GEMs in microbial community analysis. Starting with GEM repositories, automatic GEM reconstruction tools, and quality control of models, this review will give insights into microbe-microbe and microbe-host interaction predictions and optimization of microbial community models. Recent studies that utilize microbial GEMs and personalized models to infer the influence of microbiota on human diseases such as inflammatory bowel diseases (IBD) or Parkinson's disease are exemplified. Being powerful system biology tools for both species-level and community-level analysis of microbes, GEMs integrated with omics data and machine learning techniques will be indispensable for studying the microbiome and their effects on human physiology as well as for deciphering the mechanisms behind human diseases.