Multiscale Multiobjective Systems Analysis (MiMoSA): an advanced metabolic modeling framework for complex systems

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.

Investigating the Unique Ability of Trichodesmium To Fix Carbon and Nitrogen Simultaneously Using MiMoSA

The open ocean is an extremely competitive environment, partially due to the dearth of nutrients. Trichodesmium erythraeum, a marine diazotrophic cyanobacterium, is a keystone species in the ocean due to its ability to fix nitrogen and leak 30 to 50% into the surrounding environment, providing a valuable source of a necessary macronutrient to other species. While there are other diazotrophic cyanobacteria that play an important role in the marine nitrogen cycle, Trichodesmium is unique in its ability to fix both carbon and nitrogen simultaneously during the day without the use of specialized cells called heterocysts to protect nitrogenase from oxygen. Here, we use the advanced modeling framework called multiscale multiobjective systems analysis (MiMoSA) to investigate how Trichodesmium erythraeum can reduce dimolecular nitrogen to ammonium in the presence of oxygen. Our simulations indicate that nitrogenase inhibition is best modeled as Michealis-Menten competitive inhibition and that cells along the filament maintain microaerobia using high flux through Mehler reactions in order to protect nitrogenase from oxygen. We also examined the effect of location on metabolic flux and found that cells at the end of filaments operate in distinctly different metabolic modes than internal cells despite both operating in a photoautotrophic mode. These results give us important insight into how this species is able to operate photosynthesis and nitrogen fixation simultaneously, giving it a distinct advantage over other diazotrophic cyanobacteria because they can harvest light directly to fuel the energy demand of nitrogen fixation. IMPORTANCE Trichodesmium erythraeum is a marine cyanobacterium responsible for approximately half of all biologically fixed nitrogen, making it an integral part of the global nitrogen cycle. Interestingly, unlike other nitrogen-fixing cyanobacteria, Trichodesmium does not use temporal or spatial separation to protect nitrogenase from oxygen poisoning; instead, it operates photosynthesis and nitrogen fixation reactions simultaneously during the day. Unfortunately, the exact mechanism the cells utilize to operate carbon and nitrogen fixation simultaneously is unknown. Here, we use an advanced metabolic modeling framework to investigate and identify the most likely mechanisms Trichodesmium uses to protect nitrogenase from oxygen. The model predicts that cells operate in a microaerobic mode, using both respiratory and Mehler reactions to dramatically reduce intracellular oxygen concentrations.

Modeling approaches for probing cross-feeding interactions in the human gut microbiome

Microbial communities perform emergent activities that are essentially different from those carried by their individual members. The gut microbiome and its metabolites have a significant impact on the host, contributing to homeostasis or disease. Food molecules shape this community, being fermented through cross-feeding interactions of metabolites such as lactate, acetate, and amino acids, or products derived from macromolecule degradation. Mathematical and experimental approaches have been applied to understand and predict the interactions between microorganisms in complex communities such as the gut microbiota. Rational and mechanistic understanding of microbial interactions is essential to exploit their metabolic activities and identify keystone taxa and metabolites. The latter could be used in turn to modulate or replicate the metabolic behavior of the community in different contexts. This review aims to highlight recent experimental and modeling approaches for studying cross-feeding interactions within the gut microbiome. We focus on short-chain fatty acid production and fiber fermentation, which are fundamental processes in human health and disease. Special attention is paid to modeling approaches, particularly kinetic and genome-scale stoichiometric models of metabolism, to integrate experimental data under different diet and health conditions. Finally, we discuss limitations and challenges for the broad application of these modeling approaches and their experimental verification for improving our understanding of the mechanisms of microbial interactions.

Metagenome-Scale Metabolic Network Suggests Folate Produced by Bifidobacterium longum Might Contribute to High-Fiber-Diet-Induced Weight Loss in a Prader-Willi Syndrome Child

Gut-microbiota-targeted nutrition intervention has achieved success in the management of obesity, but its underlying mechanism still needs extended exploration. An obese Prader-Willi syndrome boy lost 25.8 kg after receiving a high-fiber dietary intervention for 105 days. The fecal microbiome sequencing data taken from the boy on intervention days 0, 15, 30, 45, 60, 75, and 105, along with clinical indexes, were used to construct a metagenome-scale metabolic network. Firstly, the abundances of the microbial strains were obtained by mapping the sequencing reads onto the assembly of gut organisms through use of reconstruction and analysis (AGORA) genomes. The nutritional components of the diet were obtained through the Virtual Metabolic Human database. Then, a community model was simulated using the Microbiome Modeling Toolbox. Finally, the significant Spearman correlations among the metabolites and the clinical indexes were screened and the strains that were producing these metabolites were identified. The high-fiber diet reduced the overall amount of metabolite secretions, but the secretions of folic acid derivatives by Bifidobacterium longum strains were increased and were significantly relevant to the observed weight loss. Reduced metabolites might also have directly contributed to the weight loss or indirectly contribute by enhancing leptin and decreasing adiponectin. Metagenome-scale metabolic network technology provides a cost-efficient solution for screening the functional microbial strains and metabolic pathways that are responding to nutrition therapy.

Application of attached algae flow-ways for coupling biomass production with the utilization of dilute non-point source nutrients in the Upper Laguna Madre, TX

The purpose of this study is to determine the potential for an attached algae flow-way system to efficiently produce algal biomass in estuarine surface waters by utilizing dilute non-point source nutrients from local urban, industrial, and agricultural discharges into the Upper Laguna Madre, Corpus Christi, Texas. The study was conducted over the course of two years to establish seasonal base-line biomass productivity and composition for bioproducts applications, and to identify key environmental factors and flow-way cohorts impacting biomass production. For the entire cultivation period, continuous ash-free biomass production at 4 to 10 g/m²/day (corresponding to nutrient recovery at 300 to 500 mg of nitrogen/m²/day and 15 to 30 mg of phosphorus/m²/day) was successfully achieved without system restart. Upon start-up, a latency period was observed which indicates roles for species succession from relatively low productivity, high ash content pioneer periphytic culture composed primarily of benthic diatoms from the source waters to higher productivity, reduced ash content, and more resilient culture mainly composed of filamentous chlorophyta, Ulva lactuca. Principal Component Analysis (PCA) was used to identify environmental factors driving biomass production, and machine learning (ML) models were constructed to assess the predictive capability of the data set for system performance using the local multi-season environmental variations. Environmental datasets were segregated for ML training, validation, and testing using three methods: regression tree, ensemble regression, and Gaussian process regression (GPR). The predicted ash-free biomass productivity using ML models resulted in root-squared-mean-errors (RSME) from 1.78 to 1.86 g/m²/day, and R² values from 0.67 to 0.75 using different methods. The greatest contributor to net productivity was total solar irradiation, followed by air temperature, salinity, and pH. The results of the study should be useful as a decision-making tool to application of attached algae flow-ways for biomass production while preventing algal blooms in the environment.

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.

Understanding the host-microbe interactions using metabolic modeling

The human gut harbors an enormous number of symbiotic microbes, which is vital for human health. However, interactions within the complex microbiota community and between the microbiota and its host are challenging to elucidate, limiting development in the treatment for a variety of diseases associated with microbiota dysbiosis. Using in silico simulation methods based on flux balance analysis, those interactions can be better investigated. Flux balance analysis uses an annotated genome-scale reconstruction of a metabolic network to determine the distribution of metabolic fluxes that represent the complete metabolism of a bacterium in a certain metabolic environment such as the gut. Simulation of a set of bacterial species in a shared metabolic environment can enable the study of the effect of numerous perturbations, such as dietary changes or addition of a probiotic species in a personalized manner. This review aims to introduce to experimental biologists the possible applications of flux balance analysis in the host-microbiota interaction field and discusses its potential use to improve human health. Video abstract.

Inferring composition and function of the human gut microbiome in time and space: A review of genome-scale metabolic modelling tools

The human gut hosts a complex community of microorganisms that directly influences gastrointestinal physiology, playing a central role in human health. Because of its importance, the metabolic interplay between the gut microbiome and host metabolism has gained special interest. While there has been great progress in the field driven by metagenomics and experimental studies, the mechanisms underpinning microbial composition and interactions in the microbiome remain poorly understood. Genome-scale metabolic models are mathematical structures capable of describing the metabolic potential of microbial cells. They are thus suitable tools for probing the metabolic properties of microbial communities. In this review, we discuss the most recent and relevant genome-scale metabolic modelling tools for inferring the composition, interactions, and ultimately, biological function of the constituent species of a microbial community with special emphasis in the gut microbiota. Particular attention is given to constraint-based metabolic modelling methods as well as hybrid agent-based methods for capturing the interactions and behavior of the community in time and space. Finally, we discuss the challenges hindering comprehensive modelling of complex microbial communities and its application for the in-silico design of microbial consortia with therapeutic functions.