Community standards to facilitate development and address challenges in metabolic modeling

Genome-scale metabolic model of Staphylococcus epidermidis ATCC 12228 matches in vitro conditions

Staphylococcus epidermidis, a commensal bacterium inhabiting collagen-rich areas like human skin, has gained significance due to its probiotic potential in the nasal microbiome and as a leading cause of nosocomial infections. While infrequently leading to severe illnesses, S. epidermidis exerts a significant influence, particularly in its close association with implant-related infections and its role as a classic opportunistic biofilm former. Understanding its opportunistic nature is crucial for developing novel therapeutic strategies, addressing both its beneficial and pathogenic aspects, and alleviating the burdens it imposes on patients and healthcare systems. Here, we employ genome-scale metabolic modeling as a powerful tool to elucidate the metabolic capabilities of S. epidermidis. We created a comprehensive computational resource for understanding the organism's growth conditions within diverse habitats by reconstructing and analyzing a manually curated and experimentally validated metabolic model. The final network, iSep23, incorporates 1,415 reactions, 1,051 metabolites, and 705 genes, adhering to established community standards and modeling guidelines. Benchmarking with the Metabolic Model Testing suite yields a high score, indicating the model's remarkable semantic quality. Following the findable, accessible, interoperable, and reusable (FAIR) data principles, iSep23 becomes a valuable and publicly accessible asset for subsequent studies. Growth simulations and carbon source utilization predictions align with experimental results, showcasing the model's predictive power. Ultimately, this work provides a robust foundation for future research aimed at both exploiting the probiotic potential and mitigating the pathogenic risks posed by S. epidermidis.

IMPORTANCE

Staphylococcus epidermidis, a bacterium commonly found on human skin, has shown probiotic effects in the nasal microbiome and is a notable causative agent of hospital-acquired infections. While these infections are typically non-life-threatening, their economic impact is considerable, with annual costs reaching billions of dollars in the United States. To better understand its opportunistic nature, we employed genome-scale metabolic modeling to construct a detailed network of S. epidermidis's metabolic capabilities. This model, comprising over a thousand reactions, metabolites, and genes, adheres to established standards and demonstrates solid benchmarking performance. Following the findable, accessible, interoperable, and reusable (FAIR) data principles, the model provides a valuable resource for future research. Growth simulations and predictions closely match experimental data, underscoring the model's predictive accuracy. Overall, this work lays a solid foundation for future studies aimed at leveraging the beneficial properties of S. epidermidis while mitigating its pathogenic potential.

From FAIR to CURE: Guidelines for Computational Models of Biological Systems

Guidelines for managing scientific data have been established under the FAIR principles requiring that data be Findable, Accessible, Interoperable, and Reusable. In many scientific disciplines, especially computational biology, both data and models are key to progress. For this reason, and recognizing that such models are a very special type of "data", we argue that computational models, especially mechanistic models prevalent in medicine, physiology and systems biology, deserve a complementary set of guidelines. We propose the CURE principles, emphasizing that models should be Credible, Understandable, Reproducible, and Extensible. We delve into each principle, discussing verification, validation, and uncertainty quantification for model credibility; the clarity of model descriptions and annotations for understandability; adherence to standards and open science practices for reproducibility; and the use of open standards and modular code for extensibility and reuse. We outline recommended and baseline requirements for each aspect of CURE, aiming to enhance the impact and trustworthiness of computational models, particularly in biomedical applications where credibility is paramount. Our perspective underscores the need for a more disciplined approach to modeling, aligning with emerging trends such as Digital Twins and emphasizing the importance of data and modeling standards for interoperability and reuse. Finally, we emphasize that given the non-trivial effort required to implement the guidelines, the community moves to automate as many of the guidelines as possible.

Exploring the metabolic profile of A. baumannii for antimicrobial development using genome-scale modeling

With the emergence of multidrug-resistant bacteria, the World Health Organization published a catalog of microorganisms urgently needing new antibiotics, with the carbapenem-resistant Acinetobacter baumannii designated as "critical". Such isolates, frequently detected in healthcare settings, pose a global pandemic threat. One way to facilitate a systemic view of bacterial metabolism and allow the development of new therapeutics is to apply constraint-based modeling. Here, we developed a versatile workflow to build high-quality and simulation-ready genome-scale metabolic models. We applied our workflow to create a metabolic model for A. baumannii and validated its predictive capabilities using experimental nutrient utilization and gene essentiality data. Our analysis showed that our model iACB23LX could recapitulate cellular metabolic phenotypes observed during in vitro experiments, while positive biomass production rates were observed and experimentally validated in various growth media. We further defined a minimal set of compounds that increase A. baumannii's cellular biomass and identified putative essential genes with no human counterparts, offering new candidates for future antimicrobial development. Finally, we assembled and curated the first collection of metabolic reconstructions for distinct A. baumannii strains and analyzed their growth characteristics. The presented models are in a standardized and well-curated format, enhancing their usability for multi-strain network reconstruction.

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.

Genome-scale model of Rothia mucilaginosa predicts gene essentialities and reveals metabolic capabilities

Cystic fibrosis (CF), an inherited genetic disorder caused by mutations in the cystic fibrosis transmembrane conductance regulator gene, results in sticky and thick mucosal fluids. This environment facilitates the colonization of various microorganisms, some of which can cause acute and chronic lung infections, while others may positively impact the disease. Rothia mucilaginosa, an oral commensal, is relatively abundant in the lungs of CF patients. Recent studies have unveiled its anti-inflammatory properties using in vitro three-dimensional lung epithelial cell cultures and in vivo mouse models relevant to chronic lung diseases. Apart from this, R. mucilaginosa has been associated with severe infections. However, its metabolic capabilities and genotype-phenotype relationships remain largely unknown. To gain insights into its cellular metabolism and genetic content, we developed the first manually curated genome-scale metabolic model, iRM23NL. Through growth kinetics and high-throughput phenotypic microarray testings, we defined its complete catabolic phenome. Subsequently, we assessed the model's effectiveness in accurately predicting growth behaviors and utilizing multiple substrates. We used constraint-based modeling techniques to formulate novel hypotheses that could expedite the development of antimicrobial strategies. More specifically, we detected putative essential genes and assessed their effect on metabolism under varying nutritional conditions. These predictions could offer novel potential antimicrobial targets without laborious large-scale screening of knockouts and mutant transposon libraries. Overall, iRM23NL demonstrates a solid capability to predict cellular phenotypes and holds immense potential as a valuable resource for accurate predictions in advancing antimicrobial therapies. Moreover, it can guide metabolic engineering to tailor R. mucilaginosa's metabolism for desired performance.IMPORTANCECystic fibrosis (CF) is a genetic disorder characterized by thick mucosal secretions, leading to chronic lung infections. Rothia mucilaginosa is a common bacterium found in various parts of the human body, acting as a normal part of the flora. In people with weakened immune systems, it can become an opportunistic pathogen, while it is prevalent and active in CF airways. Recent studies have highlighted its anti-inflammatory properties in the lower pulmonary system, indicating the intricate relationship between microbes and human health. Herein, we have developed the first manually curated metabolic model of R. mucilaginosa. Our study examined the previously unknown relationships between the bacterium's genotype and phenotype and identified essential genes that impact the metabolism under various conditions. With this, we opt for paving the way for developing new strategies in antimicrobial therapy and metabolic engineering, leading to enhanced therapeutic outcomes in cystic fibrosis and related conditions.

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.

The transition from genomics to phenomics in personalized population health

Modern health care faces several serious challenges, including an ageing population and its inherent burden of chronic diseases, rising costs and marginal quality metrics. By assessing and optimizing the health trajectory of each individual using a data-driven personalized approach that reflects their genetics, behaviour and environment, we can start to address these challenges. This assessment includes longitudinal phenome measures, such as the blood proteome and metabolome, gut microbiome composition and function, and lifestyle and behaviour through wearables and questionnaires. Here, we review ongoing large-scale genomics and longitudinal phenomics efforts and the powerful insights they provide into wellness. We describe our vision for the transformation of the current health care from disease-oriented to data-driven, wellness-oriented and personalized population health.

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.

Genome-scale metabolic models consistently predict in vitro characteristics of Corynebacterium striatum

Introduction: Genome-scale metabolic models (GEMs) are organism-specific knowledge bases which can be used to unravel pathogenicity or improve production of specific metabolites in biotechnology applications. However, the validity of predictions for bacterial proliferation in in vitro settings is hardly investigated. Methods: The present work combines in silico and in vitro approaches to create and curate strain-specific genome-scale metabolic models of Corynebacterium striatum. Results: We introduce five newly created strain-specific genome-scale metabolic models (GEMs) of high quality, satisfying all contemporary standards and requirements. All these models have been benchmarked using the community standard test suite Metabolic Model Testing (MEMOTE) and were validated by laboratory experiments. For the curation of those models, the software infrastructure refineGEMs was developed to work on these models in parallel and to comply with the quality standards for GEMs. The model predictions were confirmed by experimental data and a new comparison metric based on the doubling time was developed to quantify bacterial growth. Discussion: Future modeling projects can rely on the proposed software, which is independent of specific environmental conditions. The validation approach based on the growth rate calculation is now accessible and closely aligned with biological questions. The curated models are freely available via BioModels and a GitHub repository and can be used. The open-source software refineGEMs is available from https://github.com/draeger-lab/refinegems.

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.

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.

Standards, dissemination, and best practices in systems biology

The reproducibility of scientific research is crucial to the success of the scientific method. Here, we review the current best practices when publishing mechanistic models in systems biology. We recommend, where possible, to use software engineering strategies such as testing, verification, validation, documentation, versioning, iterative development, and continuous integration. In addition, adhering to the Findable, Accessible, Interoperable, and Reusable modeling principles allows other scientists to collaborate and build off of each other's work. Existing standards such as Systems Biology Markup Language, CellML, or Simulation Experiment Description Markup Language can greatly improve the likelihood that a published model is reproducible, especially if such models are deposited in well-established model repositories. Where models are published in executable programming languages, the source code and their data should be published as open-source in public code repositories together with any documentation and testing code. For complex models, we recommend container-based solutions where any software dependencies and the run-time context can be easily replicated.

Mycobacterial metabolic model development for drug target identification

Antibiotic resistance is increasing at an alarming rate, and three related mycobacteria are sources of widespread infections in humans. According to the World Health Organization, Mycobacterium leprae, which causes leprosy, is still endemic in tropical countries; Mycobacterium tuberculosis is the second leading infectious killer worldwide after COVID-19; and Mycobacteroides abscessus, a group of non-tuberculous mycobacteria, causes lung infections and other healthcare-associated infections in humans. Due to the rise in resistance to common antibacterial drugs, it is critical that we develop alternatives to traditional treatment procedures. Furthermore, an understanding of the biochemical mechanisms underlying pathogenic evolution is important for the treatment and management of these diseases. In this study, metabolic models have been developed for two bacterial pathogens, M. leprae and My. abscessus, and a new computational tool has been used to identify potential drug targets, which are referred to as bottleneck reactions. The genes, reactions, and pathways in each of these organisms have been highlighted; the potential drug targets can be further explored as broad-spectrum antibacterials and the unique drug targets for each pathogen are significant for precision medicine initiatives. The models and associated datasets described in this paper are available in GigaDB, Biomodels, and PatMeDB repositories.

Towards the human nasal microbiome: Simulating D. pigrum and S. aureus

The human nose harbors various microbes that decisively influence the wellbeing and health of their host. Among the most threatening pathogens in this habitat is Staphylococcus aureus. Multiple epidemiological studies identify Dolosigranulum pigrum as a likely beneficial bacterium based on its positive association with health, including negative associations with S. aureus. Carefully curated GEMs are available for both bacterial species that reliably simulate their growth behavior in isolation. To unravel the mutual effects among bacteria, building community models for simulating co-culture growth is necessary. However, modeling microbial communities remains challenging. This article illustrates how applying the NCMW fosters our understanding of two microbes' joint growth conditions in the nasal habitat and their intricate interplay from a metabolic modeling perspective. The resulting community model combines the latest available curated GEMs of D. pigrum and S. aureus. This uses case illustrates how to incorporate genuine GEM of participating microorganisms and creates a basic community model mimicking the human nasal environment. Our analysis supports the role of negative microbe-microbe interactions involving D. pigrum examined experimentally in the lab. By this, we identify and characterize metabolic exchange factors involved in a specific interaction between D. pigrum and S. aureus as an in silico candidate factor for a deep insight into the associated species. This method may serve as a blueprint for developing more complex microbial interaction models. Its direct application suggests new ways to prevent disease-causing infections by inhibiting the growth of pathogens such as S. aureus through microbe-microbe interactions.

Dissecting carbon metabolism of Yarrowia lipolytica type strain W29 using genome-scale metabolic modelling

Yarrowia lipolytica is a widely-used chassis cell in biotechnological applications. It has recently gained extensive research interest owing to its extraordinary ability of producing industrially valuable biochemicals from a variety of carbon sources. Genome-scale metabolic models (GSMMs) enable analyses of cellular metabolism for engineering various industrial hosts. In the present study, we developed a high-quality GSMM iYli21 for Y. lipolytica type strain W29 by extensive manual curation with Biolog experimental data. The model showed a high accuracy of 85.7% in predicting nutrient utilization. Transcriptomics data were integrated to delineate cellular metabolism of utilizing six individual metabolites as sole carbon sources. Comparisons showed that 302 reactions were commonly used, including those from TCA cycle, oxidative phosphorylation, and purine metabolism for energy and material supply. Whereas glycolytic reactions were employed only when glucose and glycerol used as sole carbon sources, gluconeogenesis and fatty acid oxidation reactions were specifically employed when fatty acid, alkane and glycerolipid were the sole carbon sources. Further test of 46 substrates for generating 5 products showed that hexanoate outcompeted other compounds in terms of maximum theoretical yield owing to the lowest carbon loss for energy supply. This newly generated model iYli21 will be a valuable tool in dissecting metabolic mechanism and guiding metabolic engineering of this important industrial cell factory.

A Genome-Scale Metabolic Model of Methanoperedens nitroreducens: Assessing Bioenergetics and Thermodynamic Feasibility

Methane is an abundant low-carbon fuel that provides a valuable energy resource, but it is also a potent greenhouse gas. Therefore, anaerobic oxidation of methane (AOM) is an essential process with central features in controlling the carbon cycle. Candidatus 'Methanoperedens nitroreducens' (M. nitroreducens) is a recently discovered methanotrophic archaeon capable of performing AOM via a reverse methanogenesis pathway utilizing nitrate as the terminal electron acceptor. Recently, reverse methanogenic pathways and energy metabolism among anaerobic methane-oxidizing archaea (ANME) have gained significant interest. However, the energetics and the mechanism for electron transport in nitrate-dependent AOM performed by M. nitroreducens is unclear. This paper presents a genome-scale metabolic model of M. nitroreducens, iMN22HE, which contains 813 reactions and 684 metabolites. The model describes its cellular metabolism and can quantitatively predict its growth phenotypes. The essentiality of the cytoplasmic heterodisulfide reductase HdrABC in the reverse methanogenesis pathway is examined by modeling the electron transfer direction and the specific energy-coupling mechanism. Furthermore, based on better understanding electron transport by modeling, a new energy transfer mechanism is suggested. The new mechanism involves reactions capable of driving the endergonic reactions in nitrate-dependent AOM, including the step reactions in reverse canonical methanogenesis and the novel electron-confurcating reaction HdrABC. The genome metabolic model not only provides an in silico tool for understanding the fundamental metabolism of ANME but also helps to better understand the reverse methanogenesis energetics and its thermodynamic feasibility.

Comparative analyses of parasites with a comprehensive database of genome-scale metabolic models

Protozoan parasites cause diverse diseases with large global impacts. Research on the pathogenesis and biology of these organisms is limited by economic and experimental constraints. Accordingly, studies of one parasite are frequently extrapolated to infer knowledge about another parasite, across and within genera. Model in vitro or in vivo systems are frequently used to enhance experimental manipulability, but these systems generally use species related to, yet distinct from, the clinically relevant causal pathogen. Characterization of functional differences among parasite species is confined to post hoc or single target studies, limiting the utility of this extrapolation approach. To address this challenge and to accelerate parasitology research broadly, we present a functional comparative analysis of 192 genomes, representing every high-quality, publicly-available protozoan parasite genome including Plasmodium, Toxoplasma, Cryptosporidium, Entamoeba, Trypanosoma, Leishmania, Giardia, and other species. We generated an automated metabolic network reconstruction pipeline optimized for eukaryotic organisms. These metabolic network reconstructions serve as biochemical knowledgebases for each parasite, enabling qualitative and quantitative comparisons of metabolic behavior across parasites. We identified putative differences in gene essentiality and pathway utilization to facilitate the comparison of experimental findings and discovered that phylogeny is not the sole predictor of metabolic similarity. This knowledgebase represents the largest collection of genome-scale metabolic models for both pathogens and eukaryotes; with this resource, we can predict species-specific functions, contextualize experimental results, and optimize selection of experimental systems for fastidious species.

Genome-scale modeling of yeast metabolism: retrospectives and perspectives

Yeasts have been widely used for production of bread, beer and wine, as well as for production of bioethanol, but they have also been designed as cell factories to produce various chemicals, advanced biofuels and recombinant proteins. To systematically understand and rationally engineer yeast metabolism, genome-scale metabolic models (GEMs) have been reconstructed for the model yeast Saccharomyces cerevisiae and nonconventional yeasts. Here, we review the historical development of yeast GEMs together with their recent applications, including metabolic flux prediction, cell factory design, culture condition optimization and multi-yeast comparative analysis. Furthermore, we present an emerging effort, namely the integration of proteome constraints into yeast GEMs, resulting in models with improved performance. At last, we discuss challenges and perspectives on the development of yeast GEMs and the integration of proteome constraints.

Elucidating Plant-Microbe-Environment Interactions Through Omics-Enabled Metabolic Modelling Using Synthetic Communities

With a growing world population and increasing frequency of climate disturbance events, we are in dire need of methods to improve plant productivity, resilience, and resistance to both abiotic and biotic stressors, both for agriculture and conservation efforts. Microorganisms play an essential role in supporting plant growth, environmental response, and susceptibility to disease. However, understanding the specific mechanisms by which microbes interact with each other and with plants to influence plant phenotypes is a major challenge due to the complexity of natural communities, simultaneous competition and cooperation effects, signalling interactions, and environmental impacts. Synthetic communities are a major asset in reducing the complexity of these systems by simplifying to dominant components and isolating specific variables for controlled experiments, yet there still remains a large gap in our understanding of plant microbiome interactions. This perspectives article presents a brief review discussing ways in which metabolic modelling can be used in combination with synthetic communities to continue progress toward understanding the complexity of plant-microbe-environment interactions. We highlight the utility of metabolic models as applied to a community setting, identify different applications for both flux balance and elementary flux mode simulation approaches, emphasize the importance of ecological theory in guiding data interpretation, and provide ideas for how the integration of metabolic modelling techniques with big data may bridge the gap between simplified synthetic communities and the complexity of natural plant-microbe systems.

eQuilibrator 3.0: a database solution for thermodynamic constant estimation

eQuilibrator (equilibrator.weizmann.ac.il) is a database of biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts around 1,000 distinct monthly users, its design could not accommodate larger compound databases and it lacked a scalable Application Programming Interface (API) for integration into other tools developed by the systems biology community. Here, we report on the recent updates to the database as well as the addition of a new Python-based interface to eQuilibrator that adds many new features such as a 100-fold larger compound database, the ability to add novel compounds, improvements in speed and memory use, and correction for Mg²⁺ ion concentrations. Moreover, the new interface can compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modelling. We foresee that these improvements will make thermodynamic modelling more accessible and facilitate the integration of eQuilibrator into other software platforms.

High-Quality Genome-Scale Reconstruction of Corynebacterium glutamicum ATCC 13032

Corynebacterium glutamicum belongs to the microbes of enormous biotechnological relevance. In particular, its strain ATCC 13032 is a widely used producer of L-amino acids at an industrial scale. Its apparent robustness also turns it into a favorable platform host for a wide range of further compounds, mainly because of emerging bio-based economies. A deep understanding of the biochemical processes in C. glutamicum is essential for a sustainable enhancement of the microbe's productivity. Computational systems biology has the potential to provide a valuable basis for driving metabolic engineering and biotechnological advances, such as increased yields of healthy producer strains based on genome-scale metabolic models (GEMs). Advanced reconstruction pipelines are now available that facilitate the reconstruction of GEMs and support their manual curation. This article presents iCGB21FR, an updated and unified GEM of C. glutamicum ATCC 13032 with high quality regarding comprehensiveness and data standards, built with the latest modeling techniques and advanced reconstruction pipelines. It comprises 1042 metabolites, 1539 reactions, and 805 genes with detailed annotations and database cross-references. The model validation took place using different media and resulted in realistic growth rate predictions under aerobic and anaerobic conditions. The new GEM produces all canonical amino acids, and its phenotypic predictions are consistent with laboratory data. The in silico model proved fruitful in adding knowledge to the metabolism of C. glutamicum: iCGB21FR still produces L-glutamate with the knock-out of the enzyme pyruvate carboxylase, despite the common belief to be relevant for the amino acid's production. We conclude that integrating high standards into the reconstruction of GEMs facilitates replicating validated knowledge, closing knowledge gaps, and making it a useful basis for metabolic engineering. The model is freely available from BioModels Database under identifier MODEL2102050001.

SBMLWebApp: Web-Based Simulation, Steady-State Analysis, and Parameter Estimation of Systems Biology Models

In systems biology, biological phenomena are often modeled by Ordinary Differential Equations (ODEs) and distributed in the de facto standard file format SBML. The primary analyses performed with such models are dynamic simulation, steady-state analysis, and parameter estimation. These methodologies are mathematically formalized, and libraries for such analyses have been published. Several tools exist to create, simulate, or visualize models encoded in SBML. However, setting up and establishing analysis environments is a crucial hurdle for non-modelers. Therefore, easy access to perform fundamental analyses of ODE models is a significant challenge. We developed SBMLWebApp, a web-based service to execute SBML-based simulation, steady-state analysis, and parameter estimation directly in the browser without the need for any setup or prior knowledge to address this issue. SBMLWebApp visualizes the result and numerical table of each analysis and provides a download of the results. SBMLWebApp allows users to select and analyze SBML models directly from the BioModels Database. Taken together, SBMLWebApp provides barrier-free access to an SBML analysis environment for simulation, steady-state analysis, and parameter estimation for SBML models. SBMLWebApp is implemented in Java™ based on an Apache Tomcat® web server using COPASI, the Systems Biology Simulation Core Library (SBSCL), and LibSBMLSim as simulation engines. SBMLWebApp is licensed under MIT with source code freely available. At the end of this article, the Data Availability Statement gives the internet links to the two websites to find the source code and run the program online.

An updated genome-scale metabolic network reconstruction of Pseudomonas aeruginosa PA14 to characterize mucin-driven shifts in bacterial metabolism

Mucins are present in mucosal membranes throughout the body and play a key role in the microbe clearance and infection prevention. Understanding the metabolic responses of pathogens to mucins will further enable the development of protective approaches against infections. We update the genome-scale metabolic network reconstruction (GENRE) of one such pathogen, Pseudomonas aeruginosa PA14, through metabolic coverage expansion, format update, extensive annotation addition, and literature-based curation to produce iPau21. We then validate iPau21 through MEMOTE, growth rate, carbon source utilization, and gene essentiality testing to demonstrate its improved quality and predictive capabilities. We then integrate the GENRE with transcriptomic data in order to generate context-specific models of P. aeruginosa metabolism. The contextualized models recapitulated known phenotypes of unaltered growth and a differential utilization of fumarate metabolism, while also revealing an increased utilization of propionate metabolism upon MUC5B exposure. This work serves to validate iPau21 and demonstrate its utility for providing biological insights.

Integrated human/SARS-CoV-2 metabolic models present novel treatment strategies against COVID-19

The coronavirus disease 2019 (COVID-19) pandemic caused by the new coronavirus (SARS-CoV-2) is currently responsible for more than 3 million deaths in 219 countries across the world and with more than 140 million cases. The absence of FDA-approved drugs against SARS-CoV-2 has highlighted an urgent need to design new drugs. We developed an integrated model of the human cell and SARS-CoV-2 to provide insight into the virus' pathogenic mechanism and support current therapeutic strategies. We show the biochemical reactions required for the growth and general maintenance of the human cell, first, in its healthy state. We then demonstrate how the entry of SARS-CoV-2 into the human cell causes biochemical and structural changes, leading to a change of cell functions or cell death. A new computational method that predicts 20 unique reactions as drug targets from our models and provides a platform for future studies on viral entry inhibition, immune regulation, and drug optimisation strategies. The model is available in BioModels (https://www.ebi.ac.uk/biomodels/MODEL2007210001) and the software tool, findCPcli, that implements the computational method is available at https://github.com/findCP/findCPcli.

Curating and comparing 114 strain-specific genome-scale metabolic models of Staphylococcus aureus

Staphylococcus aureus is a high-priority pathogen causing severe infections with high morbidity and mortality worldwide. Many S. aureus strains are methicillin-resistant (MRSA) or even multi-drug resistant. It is one of the most successful and prominent modern pathogens. An effective fight against S. aureus infections requires novel targets for antimicrobial and antistaphylococcal therapies. Recent advances in whole-genome sequencing and high-throughput techniques facilitate the generation of genome-scale metabolic models (GEMs). Among the multiple applications of GEMs is drug-targeting in pathogens. Hence, comprehensive and predictive metabolic reconstructions of S. aureus could facilitate the identification of novel targets for antimicrobial therapies. This review aims at giving an overview of all available GEMs of multiple S. aureus strains. We downloaded all 114 available GEMs of S. aureus for further analysis. The scope of each model was evaluated, including the number of reactions, metabolites, and genes. Furthermore, all models were quality-controlled using MEMOTE, an open-source application with standardized metabolic tests. Growth capabilities and model similarities were examined. This review should lead as a guide for choosing the appropriate GEM for a given research question. With the information about the availability, the format, and the strengths and potentials of each model, one can either choose an existing model or combine several models to create models with even higher predictive values. This facilitates model-driven discoveries of novel antimicrobial targets to fight multi-drug resistant S. aureus strains.

Mechanistic models of microbial community metabolism

Microbial communities affect many facets of human health and well-being. Naturally occurring bacteria, whether in nature or the human body, rarely exist in isolation. A deeper understanding of the metabolic functions of these communities is now possible with emerging computational models. In this review, we summarize frameworks for constructing mechanistic models of microbial community metabolism and discuss available algorithms for model analysis. We highlight essential decision points that greatly influence algorithm selection, as well as model analysis. Polymicrobial metabolic models can be utilized to gain insights into host-pathogen interactions, bacterial engineering, and many more translational applications.

Publishing reproducible dynamic kinetic models

Publishing repeatable and reproducible computational models is a crucial aspect of the scientific method in computational biology and one that is often forgotten in the rush to publish. The pressures of academic life and the lack of any reward system at institutions, granting agencies and journals means that publishing reproducible science is often either non-existent or, at best, presented in the form of an incomplete description. In the article, we will focus on repeatability and reproducibility in the systems biology field where a great many published models cannot be reproduced and in many cases even repeated. This review describes the current landscape of software tooling, model repositories, model standards and best practices for publishing repeatable and reproducible kinetic models. The review also discusses possible future remedies including working more closely with journals to help reviewers and editors ensure that published kinetic models are at minimum, repeatable. Contact: hsauro@uw.edu.

First Genome-Scale Metabolic Model of Dolosigranulum pigrum Confirms Multiple Auxotrophies

Dolosigranulum pigrum is a quite recently discovered Gram-positive coccus. It has gained increasing attention due to its negative correlation with Staphylococcus aureus, which is one of the most successful modern pathogens causing severe infections with tremendous morbidity and mortality due to its multiple resistances. As the possible mechanisms behind its inhibition of S. aureus remain unclear, a genome-scale metabolic model (GEM) is of enormous interest and high importance to better study its role in this fight. This article presents the first GEM of D. pigrum, which was curated using automated reconstruction tools and extensive manual curation steps to yield a high-quality GEM. It was evaluated and validated using all currently available experimental data of D. pigrum. With this model, already predicted auxotrophies and biosynthetic pathways could be verified. The model was used to define a minimal medium for further laboratory experiments and to predict various carbon sources' growth capacities. This model will pave the way to better understand D. pigrum's role in the fight against S. aureus.

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

MetaNetX/MNXref: unified namespace for metabolites and biochemical reactions in the context of metabolic models

MetaNetX/MNXref is a reconciliation of metabolites and biochemical reactions providing cross-links between major public biochemistry and Genome-Scale Metabolic Network (GSMN) databases. The new release brings several improvements with respect to the quality of the reconciliation, with particular attention dedicated to preserving the intrinsic properties of GSMN models. The MetaNetX website (https://www.metanetx.org/) provides access to the full database and online services. A major improvement is for mapping of user-provided GSMNs to MXNref, which now provides diagnostic messages about model content. In addition to the website and flat files, the resource can now be accessed through a SPARQL endpoint (https://rdf.metanetx.org).

MASSpy: Building, simulating, and visualizing dynamic biological models in Python using mass action kinetics

Mathematical models of metabolic networks utilize simulation to study system-level mechanisms and functions. Various approaches have been used to model the steady state behavior of metabolic networks using genome-scale reconstructions, but formulating dynamic models from such reconstructions continues to be a key challenge. Here, we present the Mass Action Stoichiometric Simulation Python (MASSpy) package, an open-source computational framework for dynamic modeling of metabolism. MASSpy utilizes mass action kinetics and detailed chemical mechanisms to build dynamic models of complex biological processes. MASSpy adds dynamic modeling tools to the COnstraint-Based Reconstruction and Analysis Python (COBRApy) package to provide an unified framework for constraint-based and kinetic modeling of metabolic networks. MASSpy supports high-performance dynamic simulation through its implementation of libRoadRunner: the Systems Biology Markup Language (SBML) simulation engine. Three examples are provided to demonstrate how to use MASSpy: (1) a validation of the MASSpy modeling tool through dynamic simulation of detailed mechanisms of enzyme regulation; (2) a feature demonstration using a workflow for generating ensemble of kinetic models using Monte Carlo sampling to approximate missing numerical values of parameters and to quantify biological uncertainty, and (3) a case study in which MASSpy is utilized to overcome issues that arise when integrating experimental data with the computation of functional states of detailed biological mechanisms. MASSpy represents a powerful tool to address challenges that arise in dynamic modeling of metabolic networks, both at small and large scales.