A systematic assessment of current genome-scale metabolic reconstruction tools

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

Systems-Level Modeling for CRISPR-Based Metabolic Engineering

The CRISPR-Cas system has enabled the development of sophisticated, multigene metabolic engineering programs through the use of guide RNA-directed activation or repression of target genes. To optimize biosynthetic pathways in microbial systems, we need improved models to inform design and implementation of transcriptional programs. Recent progress has resulted in new modeling approaches for identifying gene targets and predicting the efficacy of guide RNA targeting. Genome-scale and flux balance models have successfully been applied to identify targets for improving biosynthetic production yields using combinatorial CRISPR-interference (CRISPRi) programs. The advent of new approaches for tunable and dynamic CRISPR activation (CRISPRa) promises to further advance these engineering capabilities. Once appropriate targets are identified, guide RNA prediction models can lead to increased efficacy in gene targeting. Developing improved models and incorporating approaches from machine learning may be able to overcome current limitations and greatly expand the capabilities of CRISPR-Cas9 tools for metabolic engineering.

Introducing 'identification probability' for automated and transferable assessment of metabolite identification confidence in metabolomics and related studies

Methods for assessing compound identification confidence in metabolomics and related studies have been debated and actively researched for the past two decades. The earliest effort in 2007 focused primarily on mass spectrometry and nuclear magnetic resonance spectroscopy and resulted in four recommended levels of metabolite identification confidence - the Metabolite Standards Initiative (MSI) Levels. In 2014, the original MSI Levels were expanded to five levels (including two sublevels) to facilitate communication of compound identification confidence in high resolution mass spectrometry studies. Further refinement in identification levels have occurred, for example to accommodate use of ion mobility spectrometry in metabolomics workflows, and alternate approaches to communicate compound identification confidence also have been developed based on identification points schema. However, neither qualitative levels of identification confidence nor quantitative scoring systems address the degree of ambiguity in compound identifications in context of the chemical space being considered, are easily automated, or are transferable between analytical platforms. In this perspective, we propose that the metabolomics and related communities consider identification probability as an approach for automated and transferable assessment of compound identification and ambiguity in metabolomics and related studies. Identification probability is defined simply as 1/N, where N is the number of compounds in a reference library or chemical space that match to an experimentally measured molecule within user-defined measurement precision(s), for example mass measurement or retention time accuracy, etc. We demonstrate the utility of identification probability in an in silico analysis of multi-property reference libraries constructed from the Human Metabolome Database and computational property predictions, provide guidance to the community in transparent implementation of the concept, and invite the community to further evaluate this concept in parallel with their current preferred methods for assessing metabolite identification confidence.

Multi-omic characterization of antibody-producing CHO cell lines elucidates metabolic reprogramming and nutrient uptake bottlenecks

Characterizing the phenotypic diversity and metabolic capabilities of industrially relevant manufacturing cell lines is critical to bioprocess optimization and cell line development. Metabolic capabilities of production hosts limit nutrient and resource channeling into desired cellular processes and can have a profound impact on productivity. These limitations cannot be directly inferred from measured data such as spent media concentrations or transcriptomics. Here, we present an integrated multi-omic analysis pipeline combining exo-metabolomics, transcriptomics, and genome-scale metabolic network analysis and apply it to three antibody-producing Chinese Hamster Ovary cell lines to identify reprogramming features associated with high-producer clones and metabolic bottlenecks limiting product formation in an industrial bioprocess. Analysis of individual datatypes revealed a decreased nitrogenous byproduct secretion in high-producing clones and the topological changes in peripheral metabolic pathway expression associated with phase shifts. An integrated omics analysis in the context of the genome-scale metabolic model elucidated the differences in central metabolism and identified amino acid utilization bottlenecks limiting cell growth and antibody production that were not evident from exo-metabolomics or transcriptomics alone. Thus, we demonstrate the utility of a multi-omics characterization in providing an in-depth understanding of cellular metabolism, which is critical to efforts in cell engineering and bioprocess optimization.

Panera: An innovative framework for surmounting uncertainty in microbial community modeling using pan-genera metabolic models

Utilization of 16S rRNA data in constraint-based modeling to characterize microbial communities confronts a major hurdle of lack of species-level resolution, impeding the construction of community models. We introduce "Panera," an innovative framework designed to model communities under this uncertainty and yet perform metabolic inferences using pan-genus metabolic models (PGMMs). We demonstrated PGMMs' utility for comprehending the metabolic capabilities of a genus and in characterizing community models using amplicon data. The unique, adaptable nature of PGMMs unlocks their potential in building hybrid communities, combining genome-scale metabolic models (GSMMs) and PGMMs. Notably, these models provide predictions comparable to the standard GSMM-based community models, while achieving a nearly 46% reduction in error compared to the genus model-based communities. In essence, "Panera" presents a potent and effective approach to aid in metabolic modeling by enabling robust predictions of community metabolic potential when dealing with amplicon data, and offers insights into genus-level metabolic landscapes.

PCMD: A multilevel comparison database of intra- and cross-species metabolic profiling in 530 plant species

Comparative metabolomics plays a crucial role in investigating gene function, exploring metabolite evolution, and accelerating crop genetic improvement. However, a systematic platform for intra- and cross-species comparison of metabolites is currently lacking. Here, we report the Plant Comparative Metabolome Database (PCMD; http://yanglab.hzau.edu.cn/PCMD), a multilevel comparison database based on predicted metabolic profiles of 530 plant species. The PCMD serves as a platform for comparing metabolite characteristics at various levels, including species, metabolites, pathways, and biological taxonomy. The database also provides a number of user-friendly online tools, such as species comparison, metabolite enrichment, and ID conversion, enabling users to perform comparisons and enrichment analyses of metabolites across different species. In addition, the PCMD establishes a unified system based on existing metabolite-related databases by standardizing metabolite numbering. The PCMD is the most species-rich comparative plant metabolomics database currently available, and a case study demonstrates its ability to provide new insights into plant metabolic diversity.

Semi-Automatic Detection of Errors in Genome-Scale Metabolic Models

Background

Genome-Scale Metabolic Models (GSMMs) are used for numerous tasks requiring computational estimates of metabolic fluxes, from predicting novel drug targets to engineering microbes to produce valuable compounds. A key limiting step in most applications of GSMMs is ensuring their representation of the target organism's metabolism is complete and accurate. Identifying and visualizing errors in GSMMs is complicated by the fact that they contain thousands of densely interconnected reactions. Furthermore, many errors in GSMMs only become apparent when considering pathways of connected reactions collectively, as opposed to examining reactions individually.

Results

We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a collection of algorithms for detecting errors in GSMMs. The relative frequencies of errors we detect in manually curated GSMMs appear to reflect the different approaches used to curate them. Changing the method used to automatically create a GSMM from a particular organism's genome can have a larger impact on the kinds of errors in the resulting GSMM than using the same method with a different organism's genome. Our algorithms are particularly capable of identifying errors that are only apparent at the pathway level, including loops, and nontrivial cases of dead ends.

Conclusions

MACAW is capable of identifying inaccuracies of varying severity in a wide range of GSMMs. Correcting these errors can measurably improve the predictive capacity of a GSMM. The relative prevalence of each type of error we identify in a large collection of GSMMs could help shape future efforts for further automation of error correction and GSMM creation.

Microbiome modeling: a beginner's guide

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

Emerging methods for genome-scale metabolic modeling of microbial communities

Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels.

Comparative analysis of metabolic models of microbial communities reconstructed from automated tools and consensus approaches

Genome-scale metabolic models (GEMs) of microbial communities offer valuable insights into the functional capabilities of their members and facilitate the exploration of microbial interactions. These models are generated using different automated reconstruction tools, each relying on different biochemical databases that may affect the conclusions drawn from the in silico analysis. One way to address this problem is to employ a consensus reconstruction method that combines the outcomes of different reconstruction tools. Here, we conducted a comparative analysis of community models reconstructed from three automated tools, i.e. CarveMe, gapseq, and KBase, alongside a consensus approach, utilizing metagenomics data from two marine bacterial communities. Our analysis revealed that these reconstruction approaches, while based on the same genomes, resulted in GEMs with varying numbers of genes and reactions as well as metabolic functionalities, attributed to the different databases employed. Further, our results indicated that the set of exchanged metabolites was more influenced by the reconstruction approach rather than the specific bacterial community investigated. This observation suggests a potential bias in predicting metabolite interactions using community GEMs. We also showed that consensus models encompassed a larger number of reactions and metabolites while concurrently reducing the presence of dead-end metabolites. Therefore, the usage of consensus models allows making full and unbiased use from aggregating genes from the different reconstructions in assessing the functional potential of microbial communities.

Engineering Gut Symbionts: A Way to Promote Bee Growth?

Bees play a crucial role as pollinators, contributing significantly to ecosystems. However, the honeybee population faces challenges such as global warming, pesticide use, and pathogenic microorganisms. Promoting bee growth using several approaches is therefore crucial for maintaining their roles. To this end, the bacterial microbiota is well-known for its native role in supporting bee growth in several respects. Maximizing the capabilities of these microorganisms holds the theoretical potential to promote the growth of bees. Recent advancements have made it feasible to achieve this enhancement through the application of genetic engineering. In this review, we present the roles of gut symbionts in promoting bee growth and collectively summarize the engineering approaches that would be needed for future applications. Particularly, as the engineering of bee gut symbionts has not been advanced, the dominant gut symbiotic bacteria Snodgrassella alvi and Gilliamella apicola are the main focus of the paper, along with other dominant species. Moreover, we propose engineering strategies that will allow for the improvement in bee growth with listed gene targets for modification to further encourage the use of engineered gut symbionts to promote bee growth.

Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models

The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.

Using metabolic networks to predict cross-feeding and competition interactions between microorganisms

Understanding the interactions between microorganisms and their impact on bacterial behavior at the community level is a key research topic in microbiology. Different methods, relying on experimental or mathematical approaches based on the diverse properties of bacteria, are currently employed to study these interactions. Recently, the use of metabolic networks to understand the interactions between bacterial pairs has increased, highlighting the relevance of this approach in characterizing bacteria. In this study, we leverage the representation of bacteria through their metabolic networks to build a predictive model aimed at reducing the number of experimental assays required for designing bacterial consortia with specific behaviors. Our novel method for predicting cross-feeding or competition interactions between pairs of microorganisms utilizes metabolic network features. Machine learning classifiers are employed to determine the type of interaction from automatically reconstructed metabolic networks. Several algorithms were assessed and selected based on comprehensive testing and careful separation of manually compiled data sets obtained from literature sources. We used different classification algorithms, including K Nearest Neighbors, XGBoost, Support Vector Machine, and Random Forest, tested different parameter values, and implemented several data curation approaches to reduce the biological bias associated with our data set, ultimately achieving an accuracy of over 0.9. Our method holds substantial potential to advance the understanding of community behavior and contribute to the development of more effective approaches for consortia design.IMPORTANCEUnderstanding bacterial interactions at the community level is critical for microbiology, and leveraging metabolic networks presents an efficient and effective approach. The introduction of this novel method for predicting interactions through machine learning classifiers has the potential to advance the field by reducing the number of experimental assays required and contributing to the development of more effective bacterial consortia.

Integrating biological knowledge for mechanistic inference in the host-associated microbiome

Advances in high-throughput technologies have enhanced our ability to describe microbial communities as they relate to human health and disease. Alongside the growth in sequencing data has come an influx of resources that synthesize knowledge surrounding microbial traits, functions, and metabolic potential with knowledge of how they may impact host pathways to influence disease phenotypes. These knowledge bases can enable the development of mechanistic explanations that may underlie correlations detected between microbial communities and disease. In this review, we survey existing resources and methodologies for the computational integration of broad classes of microbial and host knowledge. We evaluate these knowledge bases in their access methods, content, and source characteristics. We discuss challenges of the creation and utilization of knowledge bases including inconsistency of nomenclature assignment of taxa and metabolites across sources, whether the biological entities represented are rooted in ontologies or taxonomies, and how the structure and accessibility limit the diversity of applications and user types. We make this information available in a code and data repository at: https://github.com/lozuponelab/knowledge-source-mappings. Addressing these challenges will allow for the development of more effective tools for drawing from abundant knowledge to find new insights into microbial mechanisms in disease by fostering a systematic and unbiased exploration of existing information.

Genome-scale metabolic modeling of the human milk oligosaccharide utilization by Bifidobacterium longum subsp. infantis

Bifidobacterium longum subsp. infantis is a representative and dominant species in the infant gut and is considered a beneficial microbe. This organism displays multiple adaptations to thrive in the infant gut, regarded as a model for human milk oligosaccharides (HMOs) utilization. These carbohydrates are abundant in breast milk and include different molecules based on lactose. They contain fucose, sialic acid, and N-acetylglucosamine. Bifidobacterium metabolism is complex, and a systems view of relevant metabolic pathways and exchange metabolites during HMO consumption is missing. To address this limitation, a refined genome-scale network reconstruction of this bacterium is presented using a previous reconstruction of B. infantis ATCC 15967 as a template. The latter was expanded based on an extensive revision of genome annotations, current literature, and transcriptomic data integration. The metabolic reconstruction (iLR578) accounted for 578 genes, 1,047 reactions, and 924 metabolites. Starting from this reconstruction, we built context-specific genome-scale metabolic models using RNA-seq data from cultures growing in lactose and three HMOs. The models revealed notable differences in HMO metabolism depending on the functional characteristics of the substrates. Particularly, fucosyl-lactose showed a divergent metabolism due to a fucose moiety. High yields of lactate and acetate were predicted under growth rate maximization in all conditions, whereas formate, ethanol, and 1,2-propanediol were substantially lower. Similar results were also obtained under near-optimal growth on each substrate when varying the empirically observed acetate-to-lactate production ratio. Model predictions displayed reasonable agreement between central carbon metabolism fluxes and expression data across all conditions. Flux coupling analysis revealed additional connections between succinate exchange and arginine and sulfate metabolism and a strong coupling between central carbon reactions and adenine metabolism. More importantly, specific networks of coupled reactions under each carbon source were derived and analyzed. Overall, the presented network reconstruction constitutes a valuable platform for probing the metabolism of this prominent infant gut bifidobacteria.IMPORTANCEThis work presents a detailed reconstruction of the metabolism of Bifidobacterium longum subsp. infantis, a prominent member of the infant gut microbiome, providing a systems view of its metabolism of human milk oligosaccharides.

mergem: merging, comparing, and translating genome-scale metabolic models using universal identifiers

Numerous methods exist to produce and refine genome-scale metabolic models. However, due to the use of incompatible identifier systems for metabolites and reactions, computing and visualizing the metabolic differences and similarities of such models is a current challenge. Furthermore, there is a lack of automated tools that can combine the strengths of multiple reconstruction pipelines into a curated single comprehensive model by merging different drafts, which possibly use incompatible namespaces. Here we present mergem, a novel method to compare, merge, and translate two or more metabolic models. Using a universal metabolic identifier mapping system constructed from multiple metabolic databases, mergem robustly can compare models from different pipelines, merge their common elements, and translate their identifiers to other database systems. mergem is implemented as a command line tool, a Python package, and on the web-application Fluxer, which allows simulating and visually comparing multiple models with different interactive flux graphs. The ability to merge, compare, and translate diverse genome scale metabolic models can facilitate the curation of comprehensive reconstructions and the discovery of unique and common metabolic features among different organisms.

Microbial Pathway Thermodynamics: Stoichiometric Models Unveil Anabolic and Catabolic Processes

The biotechnological exploitation of microorganisms enables the use of metabolism for the production of economically valuable substances, such as drugs or food. It is, thus, unsurprising that the investigation of microbial metabolism and its regulation has been an active research field for many decades. As a result, several theories and techniques were developed that allow for the prediction of metabolic fluxes and yields as biotechnologically relevant output parameters. One important approach is to derive macrochemical equations that describe the overall metabolic conversion of an organism and basically treat microbial metabolism as a black box. The opposite approach is to include all known metabolic reactions of an organism to assemble a genome-scale metabolic model. Interestingly, both approaches are rather successful at characterizing and predicting the expected product yield. Over the years, macrochemical equations especially have been extensively characterized in terms of their thermodynamic properties. However, a common challenge when characterizing microbial metabolism by a single equation is to split this equation into two, describing the two modes of metabolism, anabolism and catabolism. Here, we present strategies to systematically identify separate equations for anabolism and catabolism. Based on metabolic models, we systematically identify all theoretically possible catabolic routes and determine their thermodynamic efficiency. We then show how anabolic routes can be derived, and we use these to approximate biomass yield. Finally, we challenge the view of metabolism as a linear energy converter, in which the free energy gradient of catabolism drives the anabolic reactions.

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.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

State-of the-Art Constraint-Based Modeling of Microbial Metabolism: From Basics to Context-Specific Models with a Focus on Methanotrophs

Methanotrophy is the ability of an organism to capture and utilize the greenhouse gas, methane, as a source of energy-rich carbon. Over the years, significant progress has been made in understanding of mechanisms for methane utilization, mostly in bacterial systems, including the key metabolic pathways, regulation and the impact of various factors (iron, copper, calcium, lanthanum, and tungsten) on cell growth and methane bioconversion. The implementation of -omics approaches provided vast amount of heterogeneous data that require the adaptation or development of computational tools for a system-wide interrogative analysis of methanotrophy. The genome-scale mathematical modeling of its metabolism has been envisioned as one of the most productive strategies for the integration of muti-scale data to better understand methane metabolism and enable its biotechnological implementation. Herein, we provide an overview of various computational strategies implemented for methanotrophic systems. We highlight functional capabilities as well as limitations of the most popular web resources for the reconstruction, modification and optimization of the genome-scale metabolic models for methane-utilizing bacteria.

Evaluating E. coli genome-scale metabolic model accuracy with high-throughput mutant fitness data

The Escherichia coli genome-scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision-recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene-protein-reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high-throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

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.

Improving pathway prediction accuracy of constraints-based metabolic network models by treating enzymes as microcompartments

Metabolic network models have become increasingly precise and accurate as the most widespread and practical digital representations of living cells. The prediction functions were significantly expanded by integrating cellular resources and abiotic constraints in recent years. However, if unreasonable modeling methods were adopted due to a lack of consideration of biological knowledge, the conflicts between stoichiometric and other constraints, such as thermodynamic feasibility and enzyme resource availability, would lead to distorted predictions. In this work, we investigated a prediction anomaly of EcoETM, a constraints-based metabolic network model, and introduced the idea of enzyme compartmentalization into the analysis process. Through rational combination of reactions, we avoid the false prediction of pathway feasibility caused by the unrealistic assumption of free intermediate metabolites. This allowed us to correct the pathway structures of l-serine and l-tryptophan. A specific analysis explains the application method of the EcoETM-like model and demonstrates its potential and value in correcting the prediction results in pathway structure by resolving the conflict between different constraints and incorporating the evolved roles of enzymes as reaction compartments. Notably, this work also reveals the trade-off between product yield and thermodynamic feasibility. Our work is of great value for the structural improvement of constraints-based models.

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Exploring the human microbiome - A step forward for precision medicine in breast cancer

BACKGROUND

The second most frequent cancer in the world and the most common malignancy in women is breast cancer. Breast cancer is a significant health concern in India with a high mortality-to-incidence ratio and presentation at a younger age.

RECENT FINDINGS

Recent studies have identified gut microbiota as a significant factor that can have an influence on the development, treatment, and prognosis of breast cancer. This review article aims to describe the influence of microbial dysbiosis on breast cancer occurrence and the possible interactions between oncobiome and specific breast cancer molecular subtypes. The review further also discusses the role of epigenetics and diet/nutrition in the regulation of the gut and breast microbiome and its association with breast cancer prevention, therapy, and recurrence. Additionally, the recent technological advances in microbiome research, including next-generation sequencing (NGS) technologies, genome sequencing, single-cell sequencing, and microbial metabolomics along with recent advances in artificial intelligence (AI) have also been reviewed. This is an attempt to present a comprehensive status of the microbiome as a key cancer biomarker.

CONCLUSION

We believe that correlating microbiome and carcinogenesis is important as it can provide insights into the mechanisms by which microbial dysbiosis can influence cancer development and progression, leading to the potential use of the microbiome as a tool for prognostication and personalized therapy.

An automated model annotation system (AMAS) for SBML models

MOTIVATION

Annotations of biochemical models provide details of chemical species, documentation of chemical reactions, and other essential information. Unfortunately, the vast majority of biochemical models have few, if any, annotations, or the annotations provide insufficient detail to understand the limitations of the model. The quality and quantity of annotations can be improved by developing tools that recommend annotations. For example, recommender tools have been developed for annotations of genes. Although annotating genes is conceptually similar to annotating biochemical models, there are important technical differences that make it difficult to directly apply this prior work.

RESULTS

We present AMAS, a system that predicts annotations for elements of models represented in the Systems Biology Markup Language (SBML) community standard. We provide a general framework for predicting model annotations for a query element based on a database of annotated reference elements and a match score function that calculates the similarity between the query element and reference elements. The framework is instantiated to specific element types (e.g. species, reactions) by specifying the reference database (e.g. ChEBI for species) and the match score function (e.g. string similarity). We analyze the computational efficiency and prediction quality of AMAS for species and reactions in BiGG and BioModels and find that it has subsecond response times and accuracy between 80% and 95% depending on specifics of what is predicted. We have incorporated AMAS into an open-source, pip-installable Python package that can run as a command-line tool that predicts and adds annotations to species and reactions to an SBML model.

AVAILABILITY AND IMPLEMENTATION

Our project is hosted at https://github.com/sys-bio/AMAS, where we provide examples, documentation, and source code files. Our source code is licensed under the MIT open-source license.

Bioinformatic approaches for studying the microbiome of fermented food

High-throughput DNA sequencing-based approaches continue to revolutionise our understanding of microbial ecosystems, including those associated with fermented foods. Metagenomic and metatranscriptomic approaches are state-of-the-art biological profiling methods and are employed to investigate a wide variety of characteristics of microbial communities, such as taxonomic membership, gene content and the range and level at which these genes are expressed. Individual groups and consortia of researchers are utilising these approaches to produce increasingly large and complex datasets, representing vast populations of microorganisms. There is a corresponding requirement for the development and application of appropriate bioinformatic tools and pipelines to interpret this data. This review critically analyses the tools and pipelines that have been used or that could be applied to the analysis of metagenomic and metatranscriptomic data from fermented foods. In addition, we critically analyse a number of studies of fermented foods in which these tools have previously been applied, to highlight the insights that these approaches can provide.

An integrated systems biology approach reveals differences in formate metabolism in the genus Methanothermobacter

Methanogenesis allows methanogenic archaea to generate cellular energy for their growth while producing methane. Thermophilic hydrogenotrophic species of the genus Methanothermobacter have been recognized as robust biocatalysts for a circular carbon economy and are already applied in power-to-gas technology with biomethanation, which is a platform to store renewable energy and utilize captured carbon dioxide. Here, we generated curated genome-scale metabolic reconstructions for three Methanothermobacter strains and investigated differences in the growth performance of these same strains in chemostat bioreactor experiments with hydrogen and carbon dioxide or formate as substrates. Using an integrated systems biology approach, we identified differences in formate anabolism between the strains and revealed that formate anabolism influences the diversion of carbon between biomass and methane. This finding, together with the omics datasets and the metabolic models we generated, can be implemented for biotechnological applications of Methanothermobacter in power-to-gas technology, and as a perspective, for value-added chemical production.

Bactabolize is a tool for high-throughput generation of bacterial strain-specific metabolic models

Metabolic capacity can vary substantially within a bacterial species, leading to ecological niche separation, as well as differences in virulence and antimicrobial susceptibility. Genome-scale metabolic models are useful tools for studying the metabolic potential of individuals, and with the rapid expansion of genomic sequencing there is a wealth of data that can be leveraged for comparative analysis. However, there exist few tools to construct strain-specific metabolic models at scale. Here, we describe Bactabolize, a reference-based tool which rapidly produces strain-specific metabolic models and growth phenotype predictions. We describe a pan reference model for the priority antimicrobial-resistant pathogen, Klebsiella pneumoniae, and a quality control framework for using draft genome assemblies as input for Bactabolize. The Bactabolize-derived model for K. pneumoniae reference strain KPPR1 performed comparatively or better than currently available automated approaches CarveMe and gapseq across 507 substrate and 2317 knockout mutant growth predictions. Novel draft genomes passing our systematically defined quality control criteria resulted in models with a high degree of completeness (≥99% genes and reactions captured compared to models derived from matched complete genomes) and high accuracy (mean 0.97, n=10). We anticipate the tools and framework described herein will facilitate large-scale metabolic modelling analyses that broaden our understanding of diversity within bacterial species and inform novel control strategies for priority pathogens.

Reconciliation and evolution of Penicillium rubens genome-scale metabolic networks-What about specialised metabolism?

In recent years, genome sequencing of filamentous fungi has revealed a high proportion of specialised metabolites with growing pharmaceutical interest. However, detecting such metabolites through in silico genome analysis does not necessarily guarantee their expression under laboratory conditions. However, one plausible strategy for enabling their production lies in modifying the growth conditions. Devising a comprehensive experimental design testing in different culture environments is time-consuming and expensive. Therefore, using in silico modelling as a preliminary step, such as Genome-Scale Metabolic Network (GSMN), represents a promising approach to predicting and understanding the observed specialised metabolite production in a given organism. To address these questions, we reconstructed a new high-quality GSMN for the Penicillium rubens Wisconsin 54-1255 strain, a commonly used model organism. Our reconstruction, iPrub22, adheres to current convention standards and quality criteria, incorporating updated functional annotations, orthology searches with different GSMN templates, data from previous reconstructions, and manual curation steps targeting primary and specialised metabolites. With a MEMOTE score of 74% and a metabolic coverage of 45%, iPrub22 includes 5,192 unique metabolites interconnected by 5,919 reactions, of which 5,033 are supported by at least one genomic sequence. Of the metabolites present in iPrub22, 13% are categorised as belonging to specialised metabolism. While our high-quality GSMN provides a valuable resource for investigating known phenotypes expressed in P. rubens, our analysis identifies bottlenecks related, in particular, to the definition of what is a specialised metabolite, which requires consensus within the scientific community. It also points out the necessity of accessible, standardised and exhaustive databases of specialised metabolites. These questions must be addressed to fully unlock the potential of natural product production in P. rubens and other filamentous fungi. Our work represents a foundational step towards the objective of rationalising the production of natural products through GSMN modelling.

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.

An Automated Model Annotation System (AMAS) for SBML Models

Motivation

Annotations of biochemical models provide details of chemical species, documentation of chemical reactions, and other essential information. Unfortunately, the vast majority of biochemical models have few, if any, annotations, or the annotations provide insufficient detail to understand the limitations of the model. The quality and quantity of annotations can be improved by developing tools that recommend annotations. For example, recommender tools have been developed for annotations of genes. Although annotating genes is conceptually similar to annotating biochemical models, there are important technical differences that make it difficult to directly apply this prior work.

Results

We present AMAS, a system that predicts annotations for elements of models represented in the Systems Biology Markup Language (SBML) community standard. We provide a general framework for predicting model annotations for a query element based on a database of annotated reference elements and a match score function that calculates the similarity between the query element and reference elements. The framework is instantiated to specific element types (e.g., species, reactions) by specifying the reference database (e.g., ChEBI for species) and the match score function (e.g., string similarity). We analyze the computational efficiency and prediction quality of AMAS for species and reactions in BiGG and BioModels and find that it has sub-second response times and accuracy between 80% and 95% depending on specifics of what is predicted. We have incorporated AMAS into an open-source, pip-installable Python package that can run as a command-line tool that predicts and adds annotations to species and reactions to an SBML model.

Availability

Our project is hosted at https://github.com/sys-bio/AMAS, where we provide examples, documentation, and source code files. Our source code is licensed under the MIT open-source license.

Metabolic interaction models recapitulate leaf microbiota ecology

Resource allocation affects the structure of microbiomes, including those associated with living hosts. Understanding the degree to which this dependency determines interspecies interactions may advance efforts to control host-microbiome relationships. We combined synthetic community experiments with computational models to predict interaction outcomes between plant-associated bacteria. We mapped the metabolic capabilities of 224 leaf isolates from Arabidopsis thaliana by assessing the growth of each strain on 45 environmentally relevant carbon sources in vitro. We used these data to build curated genome-scale metabolic models for all strains, which we combined to simulate >17,500 interactions. The models recapitulated outcomes observed in planta with >89% accuracy, highlighting the role of carbon utilization and the contributions of niche partitioning and cross-feeding in the assembly of leaf microbiomes.

From genotype to phenotype: computational approaches for inferring microbial traits relevant to the food industry

When selecting microbial strains for the production of fermented foods, various microbial phenotypes need to be taken into account to achieve target product characteristics, such as biosafety, flavor, texture, and health-promoting effects. Through continuous advances in sequencing technologies, microbial whole-genome sequences of increasing quality can now be obtained both cheaper and faster, which increases the relevance of genome-based characterization of microbial phenotypes. Prediction of microbial phenotypes from genome sequences makes it possible to quickly screen large strain collections in silico to identify candidates with desirable traits. Several microbial phenotypes relevant to the production of fermented foods can be predicted using knowledge-based approaches, leveraging our existing understanding of the genetic and molecular mechanisms underlying those phenotypes. In the absence of this knowledge, data-driven approaches can be applied to estimate genotype-phenotype relationships based on large experimental datasets. Here, we review computational methods that implement knowledge- and data-driven approaches for phenotype prediction, as well as methods that combine elements from both approaches. Furthermore, we provide examples of how these methods have been applied in industrial biotechnology, with special focus on the fermented food industry.

Reconstructor: a COBRApy compatible tool for automated genome-scale metabolic network reconstruction with parsimonious flux-based gap-filling

MOTIVATION

Genome-scale metabolic network reconstructions (GENREs) are valuable for understanding cellular metabolism in silico. Several tools exist for automatic GENRE generation. However, these tools frequently (i) do not readily integrate with some of the widely-used suites of packaged methods available for network analysis, (ii) lack effective network curation tools, (iii) are not sufficiently user-friendly, and (iv) often produce low-quality draft reconstructions.

RESULTS

Here, we present Reconstructor, a user-friendly, COBRApy-compatible tool that produces high-quality draft reconstructions with reaction and metabolite naming conventions that are consistent with the ModelSEED biochemistry database and includes a gap-filling technique based on the principles of parsimony. Reconstructor can generate SBML GENREs from three input types: annotated protein .fasta sequences (Type 1 input), a BLASTp output (Type 2), or an existing SBML GENRE that can be further gap-filled (Type 3). While Reconstructor can be used to create GENREs of any species, we demonstrate the utility of Reconstructor with bacterial reconstructions. We demonstrate how Reconstructor readily generates high-quality GENRES that capture strain, species, and higher taxonomic differences in functional metabolism of bacteria and are useful for further biological discovery.

AVAILABILITY AND IMPLEMENTATION

The Reconstructor Python package is freely available for download. Complete installation and usage instructions and benchmarking data are available at http://github.com/emmamglass/reconstructor.

Extracting functionally accurate context-specific models of Atlantic salmon metabolism

Constraint-based models (CBMs) are used to study metabolic network structure and function in organisms ranging from microbes to multicellular eukaryotes. Published CBMs are usually generic rather than context-specific, meaning that they do not capture differences in reaction activities, which, in turn, determine metabolic capabilities, between cell types, tissues, environments, or other conditions. Only a subset of a CBM's metabolic reactions and capabilities are likely to be active in any given context, and several methods have therefore been developed to extract context-specific models from generic CBMs through integration of omics data. We tested the ability of six model extraction methods (MEMs) to create functionally accurate context-specific models of Atlantic salmon using a generic CBM (SALARECON) and liver transcriptomics data from contexts differing in water salinity (life stage) and dietary lipids. Three MEMs (iMAT, INIT, and GIMME) outperformed the others in terms of functional accuracy, which we defined as the extracted models' ability to perform context-specific metabolic tasks inferred directly from the data, and one MEM (GIMME) was faster than the others. Context-specific versions of SALARECON consistently outperformed the generic version, showing that context-specific modeling better captures salmon metabolism. Thus, we demonstrate that results from human studies also hold for a non-mammalian animal and major livestock species.

The media composition as a crucial element in high-throughput metabolic network reconstruction

In recent years, metagenome-assembled genomes (MAGs) have provided glimpses into the intra- and interspecies genetic diversity and interactions that form the bases of complex microbial communities. High-throughput reconstruction of genome-scale metabolic networks (GEMs) from MAGs is a promising avenue to disentangle the myriad trophic interactions stabilizing these communities. However, high-throughput reconstruction of GEMs relies on accurate gap filling of metabolic pathways using automated algorithms. Here, we systematically explore how the composition of the media (specification of the available nutrients and metabolites) during gap filling influences the resulting GEMs concerning predicted auxotrophies for fully sequenced model organisms and environmental isolates. We expand this analysis by using 106 MAGs from the same species with differing quality. We find that although the completeness of MAGs influences the fraction of gap-filled reactions, the composition of the media plays the dominant role in the accurate prediction of auxotrophies that form the basis of myriad community interactions. We propose that constraining the media composition for gap filling through both experimental approaches and computational approaches will increase the reliability of high-throughput reconstruction of genome-scale metabolic models from MAGs and paves the way for culture independent prediction of trophic interactions in complex microbial communities.

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.

Multi-dimensional experimental and computational exploration of metabolism pinpoints complex probiotic interactions

Multi-strain probiotics are widely regarded as effective products for improving gut microbiota stability and host health, providing advantages over single-strain probiotics. However, in general, it is unclear to what extent different strains would cooperate or compete for resources, and how the establishment of a common biofilm microenvironment could influence their interactions. In this work, we develop an integrative experimental and computational approach to comprehensively assess the metabolic functionality and interactions of probiotics across growth conditions. Our approach combines co-culture assays with genome-scale modelling of metabolism and multivariate data analysis, thus exploiting complementary data- and knowledge-driven systems biology techniques. To show the advantages of the proposed approach, we apply it to the study of the interactions between two widely used probiotic strains of Lactobacillus reuteri and Saccharomyces boulardii, characterising their production potential for compounds that can be beneficial to human health. Our results show that these strains can establish a mixed cooperative-antagonistic interaction best explained by competition for shared resources, with an increased individual exchange but an often decreased net production of amino acids and short-chain fatty acids. Overall, our work provides a strategy that can be used to explore microbial metabolic fingerprints of biotechnological interest, capable of capturing multifaceted equilibria even in simple microbial consortia.

A Genome-Scale Metabolic Model of Marine Heterotroph Vibrio splendidus Strain 1A01

While Vibrio splendidus is best known as an opportunistic pathogen in oysters, Vibrio splendidus strain 1A01 was first identified as an early colonizer of synthetic chitin particles incubated in seawater. To gain a better understanding of its metabolism, a genome-scale metabolic model (GSMM) of V. splendidus 1A01 was reconstructed. GSMMs enable us to simulate all metabolic reactions in a bacterial cell using flux balance analysis. A draft model was built using an automated pipeline from BioCyc. Manual curation was then performed based on experimental data, in part by gap-filling metabolic pathways and tailoring the model's biomass reaction to V. splendidus 1A01. The challenges of building a metabolic model for a marine microorganism like V. splendidus 1A01 are described. IMPORTANCE A genome-scale metabolic model of V. splendidus 1A01 was reconstructed in this work. We offer solutions to the technical problems associated with model reconstruction for a marine bacterial strain like V. splendidus 1A01, which arise largely from the high salt concentration found in both seawater and culture media that simulate seawater.

Interactions between Culturable Bacteria Are Predicted by Individual Species' Growth

Predicting interspecies interactions is a key challenge in microbial ecology given that interactions shape the composition and functioning of microbial communities. However, predicting microbial interactions is challenging because they can vary considerably depending on species' metabolic capabilities and environmental conditions. Here, we employ machine learning models to predict pairwise interactions between culturable bacteria based on their phylogeny, monoculture growth capabilities, and interactions with other species. We trained our models on one of the largest available pairwise interactions data set containing over 7,500 interactions between 20 species from two taxonomic groups that were cocultured in 40 different carbon environments. Our models accurately predicted both the sign (accuracy of 88%) and the strength of effects (R² of 0.87) species had on each other's growth. Encouragingly, predictions with comparable accuracy could be made even when not relying on information about interactions with other species, which are often hard to measure. However, species' monoculture growth was essential to the model, as predictions based solely on species' phylogeny and inferred metabolic capabilities were significantly less accurate. These results bring us one step closer to a predictive understanding of microbial communities, which is essential for engineering beneficial microbial consortia. IMPORTANCE In order to understand the function and structure of microbial communities, one must know all pairwise interactions that occur between the different species within the community, as these interactions shape the community's structure and functioning. However, measuring all pairwise interactions can be an extremely difficult task especially when dealing with big complex communities. Because of that, predicting interspecies interactions is a key challenge in microbial ecology. Here, we use machine learning models in order to accurately predict the type and strength of interactions. We trained our models on one of the largest available pairwise interactions data set, containing over 7,500 interactions between 20 different species that were cocultured in 40 different environments. Our results show that, in general, accurate predictions can be made, and that the ability of each species to grow on its own in the given environment contributes the most to predictions. Being able to predict microbial interactions would put us one step closer to predicting the functionality of microbial communities and to rationally microbiome engineering.

Systematizing Microbial Bioplastic Production for Developing Sustainable Bioeconomy: Metabolic Nexus Modeling, Economic and Environmental Technologies Assessment

The excessive usage of non-renewable resources to produce plastic commodities has incongruously influenced the environment's health. Especially in the times of COVID-19, the need for plastic-based health products has increased predominantly. Given the rise in global warming and greenhouse gas emissions, the lifecycle of plastic has been established to contribute to it significantly. Bioplastics such as polyhydroxy alkanoates, polylactic acid, etc. derived from renewable energy origin have been a magnificent alternative to conventional plastics and reconnoitered exclusively for combating the environmental footprint of petrochemical plastic. However, the economically reasonable and environmentally friendly procedure of microbial bioplastic production has been a hard nut to crack due to less scouted and inefficient process optimization and downstream processing methodologies. Thereby, meticulous employment of computational tools such as genome-scale metabolic modeling and flux balance analysis has been practiced in recent times to understand the effect of genomic and environmental perturbations on the phenotype of the microorganism. In-silico results not only aid us in determining the biorefinery abilities of the model microorganism but also curb our reliance on equipment, raw materials, and capital investment for optimizing the best conditions. Additionally, to accomplish sustainable large-scale production of microbial bioplastic in a circular bioeconomy, extraction, and refinement of bioplastic needs to be investigated extensively by practicing techno-economic analysis and life cycle assessment. This review put forth state-of-the-art know-how on the proficiency of these computational techniques in laying the foundation of an efficient bioplastic manufacturing blueprint, chiefly focusing on microbial polyhydroxy alkanoates (PHA) production and its efficacy in outplacing fossil based plastic products.

Metatranscriptomics-guided genome-scale metabolic modeling of microbial communities

Multi-omics data integration via mechanistic models of metabolism is a scalable and flexible framework for exploring biological hypotheses in microbial systems. However, although most microorganisms are unculturable, such multi-omics modeling is limited to isolate microbes or simple synthetic communities. Here, we developed an approach for modeling microbial activity and interactions that leverages the reconstruction of metagenome-assembled genomes and associated genome-centric metatranscriptomes. At its core, we designed a method for condition-specific metabolic modeling of microbial communities through the integration of metatranscriptomic data. Using this approach, we explored the behavior of anaerobic digestion consortia driven by hydrogen availability and human gut microbiota dysbiosis associated with Crohn's disease, identifying condition-dependent amino acid requirements in archaeal species and a reduced short-chain fatty acid exchange network associated with disease, respectively. Our approach can be applied to complex microbial communities, allowing a mechanistic contextualization of multi-omics data on a metagenome scale.

Context-Specific Genome-Scale Metabolic Modelling and Its Application to the Analysis of COVID-19 Metabolic Signatures

Genome-scale metabolic models (GEMs) have found numerous applications in different domains, ranging from biotechnology to systems medicine. Herein, we overview the most popular algorithms for the automated reconstruction of context-specific GEMs using high-throughput experimental data. Moreover, we describe different datasets applied in the process, and protocols that can be used to further automate the model reconstruction and validation. Finally, we describe recent COVID-19 applications of context-specific GEMs, focusing on the analysis of metabolic implications, identification of biomarkers and potential drug targets.

Investigation of two metabolic engineering approaches for (R,R)-2,3-butanediol production from glycerol in Bacillus subtilis

BACKGROUND

Flux Balance Analysis (FBA) is a well-known bioinformatics tool for metabolic engineering design. Previously, we have successfully used single-level FBA to design metabolic fluxes in Bacillus subtilis to enhance (R,R)-2,3-butanediol (2,3-BD) production from glycerol. OptKnock is another powerful technique for devising gene deletion strategies to maximize microbial growth coupling with improved biochemical production. It has never been used in B. subtilis. In this study, we aimed to compare the use of single-level FBA and OptKnock for designing enhanced 2,3-BD production from glycerol in B. subtilis.

RESULTS

Single-level FBA and OptKnock were used to design metabolic engineering approaches for B. subtilis to enhance 2,3-BD production from glycerol. Single-level FBA indicated that deletion of ackA, pta, lctE, and mmgA would improve the production of 2,3-BD from glycerol, while OptKnock simulation suggested the deletion of ackA, pta, mmgA, and zwf. Consequently, strains LM01 (single-level FBA-based) and MZ02 (OptKnock-based) were constructed, and their capacity to produce 2,3-BD from glycerol was investigated. The deletion of multiple genes did not negatively affect strain growth and glycerol utilization. The highest 2,3-BD production was detected in strain LM01. Strain MZ02 produced 2,3-BD at a similar level as the wild type, indicating that the OptKnock prediction was erroneous. Two-step FBA was performed to examine the reason for the erroneous OptKnock prediction. Interestingly, we newly found that zwf gene deletion in strain MZ02 improved lactate production, which has never been reported to date. The predictions of single-level FBA for strain MZ02 were in line with experimental findings.

CONCLUSIONS

We showed that single-level FBA is an effective approach for metabolic design and manipulation to enhance 2,3-BD production from glycerol in B. subtilis. Further, while this approach predicted the phenotypes of generated strains with high precision, OptKnock prediction was not accurate. We suggest that OptKnock modelling predictions be evaluated by using single-level FBA to ensure the accuracy of metabolic pathway design. Furthermore, the zwf gene knockout resulted in the change of metabolic fluxes to enhance the lactate productivity.

Genome-Scale Modeling and Systems Metabolic Engineering of Vibrio natriegens for the Production of 1,3-Propanediol

The fastest-growing bacterium Vibrio natriegens is a highly promising next-generation workhorse for molecular biology and industrial biotechnology. In this work, we described the workflows for developing genome-scale metabolic models and genome-editing protocols for engineering Vibrio natriegens. A case study for metabolic engineering of Vibrio natriegens for the production of 1,3-propanediol was also presented.

ChiMera: an easy to use pipeline for bacterial genome based metabolic network reconstruction, evaluation and visualization

BACKGROUND

Genome-scale metabolic reconstruction tools have been developed in the last decades. They have helped to reconstruct eukaryotic and prokaryotic metabolic models, which have contributed to fields, e.g., genetic engineering, drug discovery, prediction of phenotypes, and other model-driven discoveries. However, the use of these programs requires a high level of bioinformatic skills. Moreover, the functionalities required to build models are scattered throughout multiple tools, requiring knowledge and experience for utilizing several tools.

RESULTS

Here we present ChiMera, which combines tools used for model reconstruction, prediction, and visualization. ChiMera uses CarveMe in the reconstruction module, generating a gap-filled draft reconstruction able to produce growth predictions using flux balance analysis for gram-positive and gram-negative bacteria. ChiMera also contains two modules for metabolic network visualization. The first module generates maps for the most important pathways, e.g., glycolysis, nucleotides and amino acids biosynthesis, fatty acid oxidation and biosynthesis and core-metabolism. The second module produces a genome-wide metabolic map, which can be used to retrieve KEGG pathway information for each compound in the model. A module to investigate gene essentiality and knockout is also present.

CONCLUSIONS

Overall, ChiMera uses automation algorithms to combine a variety of tools to automatically perform model creation, gap-filling, flux balance analysis (FBA), and metabolic network visualization. ChiMera models readily provide metabolic insights that can aid genetic engineering projects, prediction of phenotypes, and model-driven discoveries.

Analysis of the Propionate Metabolism in Bacillus subtilis during 3-Indolacetic Production

The genera Bacillus belongs to the group of microorganisms that are known as plant growth-promoting bacteria, their metabolism has evolved to produce molecules that benefit the growth of the plant, and the production of 3-indole acetic acid (IAA) is part of its secondary metabolism. In this work, Bacillus subtilis was cultivated in a bioreactor to produce IAA using propionate and glucose as carbon sources in an M9-modified media; in both cases, tryptophan was added as a co-substrate. The yield of IAA using propionate is 17% higher compared to glucose. After 48 h of cultivation, the final concentration was 310 mg IAA/L using propionate and 230 mg IAA/L using glucose, with a concentration of 500 mg Trp/L. To gain more insight into propionate metabolism and its advantages, the genome-scale metabolic model of B. subtilis (iBSU 1147) and computational analysis were used to calculate flux distribution and evaluate the metabolic capabilities to produce IAA using propionate. The metabolic fluxes demonstrate that propionate uptake favors the production of precursors needed for the synthesis of the hormone, and the sensitivity analysis shows that the control of a specific growth rate has a positive impact on the production of IAA.

Metabolic Modeling and Bidirectional Culturing of Two Gut Microbes Reveal Cross-Feeding Interactions and Protective Effects on Intestinal Cells

The gut microbiota is constituted by thousands of microbial interactions, some of which correspond to the exchange of metabolic by-products or cross-feeding. Inulin and xylan are two major dietary polysaccharides that are fermented by members of the human gut microbiota, resulting in different metabolic profiles. Here, we integrated community modeling and bidirectional culturing assays to study the metabolic interactions between two gut microbes, Phocaeicola dorei and Lachnoclostridium symbiosum, growing in inulin or xylan, and how they provide a protective effect in cultured cells. P. dorei (previously belonging to the Bacteroides genus) was able to consume inulin and xylan, while L. symposium only used certain inulin fractions to produce butyrate as a major end product. Constrained-based flux simulations of refined genome-scale metabolic models of both microbes predicted high lactate and succinate cross-feeding fluxes between P. dorei and L. symbiosum when growing in each fiber. Bidirectional culture assays in both substrates revealed that L. symbiosum growth increased in the presence of P. dorei. Carbohydrate consumption analyses showed a faster carbohydrate consumption in cocultures compared to monocultures. Lactate and succinate concentrations in bidirectional cocultures were lower than in monocultures, pointing to cross-feeding as initially suggested by the model. Butyrate concentrations were similar across all conditions. Finally, supernatants from both bacteria cultured in xylan in bioreactors significantly reduced tumor necrosis factor-α-induced inflammation in HT-29 cells and exerted a protective effect against the TcdB toxin in Caco-2 epithelial cells. Surprisingly, this effect was not observed in inulin cocultures. Overall, these results highlight the predictive value of metabolic models integrated with microbial culture assays for probing microbial interactions in the gut microbiota. They also provide an example of how metabolic exchange could lead to potential beneficial effects in the host. IMPORTANCE Microbial interactions represent the inner connections in the gut microbiome. By integrating mathematical modeling tools and microbial bidirectional culturing, we determined how two gut commensals engage in the exchange of cross-feeding metabolites, lactate and succinate, for increased growth in two fibers. These interactions underpinned butyrate production in cocultures, resulting in a significant reduction in cellular inflammation and protection against microbial toxins when applied to cellular models.