Motivation: The number and size of computational models in biology have drastically increased over the past years and continue to grow. Modeled networks are becoming more complex, and reconstructing them from the beginning in an exchangeable and reproducible manner is challenging. Using precisely defined ontologies enables the encoding of field-specific knowledge and the association of disparate data types. In computational modeling, the medium for representing domain knowledge is the set of orthogonal structured controlled vocabularies named Systems Biology Ontology (SBO). The SBO terms enable modelers to explicitly define and describe model entities, including their roles and characteristics. Results: Here, we present the first standalone tool that automatically assigns SBO terms to multiple entities of a given SBML model, named the SBOannotator. The main focus lies on the reactions, as the correct assignment of precise SBO annotations requires their extensive classification. Our implementation does not consider only top-level terms but examines the functionality of the underlying enzymes to allocate precise and highly specific ontology terms to biochemical reactions. Transport reactions are examined separately and are classified based on the mechanism of molecule transport. Pseudo-reactions that serve modeling purposes are given reasonable terms to distinguish between biomass production and the import or export of metabolites. Finally, other model entities, such as metabolites and genes, are annotated with appropriate terms. Including SBO annotations in the models will enhance the reproducibility, usability, and analysis of biochemical networks. Availability and implementation: SBOannotator is freely available from https://github.com/draeger-lab/SBOannotator/.

Genome-scale metabolic models (GEMs) and other constraint-based models (CBMs) play a pivotal role in understanding biological phenotypes and advancing research in areas like metabolic engineering, human disease modelling, drug discovery, and personalized medicine. Despite their growing application, a significant challenge remains in ensuring the reproducibility of GEMs, primarily due to inconsistent reporting and inadequate model documentation of model results. Addressing this gap, we introduce FROG analysis, a community-driven initiative aimed at standardizing reproducibility assessments of CBMs and GEMs. The FROG framework encompasses four key analyses—Flux variability, Reaction deletion, Objective function, and Gene deletion—to produce standardized, numerically reproducible FROG reports. These reports serve as reference datasets, enabling model evaluators, curators, and independent researchers to verify the reproducibility of GEMs systematically. BioModels, a leading repository of systems biology models, has integrated FROG analysis into its curation workflow, enhancing the reproducibility and reusability of submitted GEMs. In our study evaluating 65 GEM submissions from the community, approximately 40% reproduced without intervention, 28% requiring minor adjustments, and 32% needing input from authors. The standardization introduced by FROG analysis facilitated the detection and resolution of issues, ultimately leading to the successful reproduction of all models. By establishing a standardized and comprehensive approach to evaluating GEM reproducibility, FROG analysis significantly contributes to making CBMs and GEMs more transparent, reusable, and reliable for the broader scientific community.

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, i RM23NL. 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, i RM23NL 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. IMPORTANCE: Cystic 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.

COVID-19 is one of the deadliest respiratory diseases, and its emergence caught the pharmaceutical industry off guard. While vaccines have been rapidly developed, treatment options for infected people remain scarce, and COVID-19 poses a substantial global threat. This study presents a novel workflow to predict robust druggable targets against emerging RNA viruses using metabolic networks and information of the viral structure and its genome sequence. For this purpose, we implemented pymCADRE and PREDICATE to create tissue-specific metabolic models, construct viral biomass functions and predict host-based antiviral targets from more than one genome. We observed that pymCADRE reduces the computational time of flux variability analysis for internal optimizations. We applied these tools to create a new metabolic network of primary bronchial epithelial cells infected with SARS-CoV-2 and identified enzymatic reactions with inhibitory effects. The most promising reported targets were from the purine metabolism, while targeting the pyrimidine and carbohydrate metabolisms seemed to be promising approaches to enhance viral inhibition. Finally, we computationally tested the robustness of our targets in all known variants of concern, verifying our targets’ inhibitory effects. Since laboratory tests are time-consuming and involve complex readouts to track processes, our workflow focuses on metabolic fluxes within infected cells and is applicable for rapid hypothesis-driven identification of potentially exploitable antivirals concerning various viruses and host cell types. The recently emerged human coronavirus SARS-CoV-2 spread worldwide, causing severe challenges in health care, the economy, and society. Developing new vaccines and therapies is essential to prevent the next pandemic efficiently. However, vaccines have the disadvantage of decreased immunity over time, while they lose their efficacy against subsequent mutations and variants. Hence, effective pandemic preparedness requires discovering broadly acting antivirals with high resistance barriers. Here, we accelerate the detection of antiviral drug candidates against RNA viruses by analyzing metabolic changes in infected cells. Viruses rely on the host to acquire all macromolecules needed for their replication, re-programming the cellular metabolism according to their needs. For this reason, host-directed approaches are of great importance. We develop software to reconstruct constraint-based models, simulate infections of a cell, and identify host-based metabolic pathways that can be inhibited to suppress viral replication. We identify promising targets with inhibitory effects across multiple variants facilitating further in vitro and in vivo experiments. Our workflow can be applied to any RNA virus and aims to rapidly identify antiviral targets to better prepare for the next pandemic.

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.

The SARS-CoV-2 pandemic has reemphasized the urgent need for broad-spectrum antiviral therapies. We developed a computational workflow using scRNA-Seq data to assess cellular metabolism during viral infection. With this workflow we predicted the capacity of cells to sustain SARS-CoV-2 virion production in patients and found a tissue-wide induction of metabolic pathways that support viral replication. Expanding our analysis to influenza A and dengue viruses, we identified metabolic targets and inhibitors for potential broad-spectrum antiviral treatment. These targets were highly enriched for known interaction partners of all analyzed viruses. Indeed, phenformin, an NADH:ubiquinone oxidoreductase inhibitor, suppressed SARS-CoV-2 and dengue virus replication. Atpenin A5, blocking succinate dehydrogenase, inhibited SARS-CoV-2, dengue virus, respiratory syncytial virus, and influenza A virus with high selectivity indices. In vivo, phenformin showed antiviral activity against SARS-CoV-2 in a Syrian hamster model. Our work establishes host metabolism as druggable for broad-spectrum antiviral strategies, providing invaluable tools for pandemic preparedness.