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

Logic programming-based Minimal Cut Sets reveal consortium-level therapeutic targets for chronic wound infections

Minimal Cut Sets (MCSs) identify sets of reactions which, when removed from a metabolic network, disable certain cellular functions. The traditional search for MCSs within genome-scale metabolic models (GSMMs) targets cellular growth, identifies reaction sets resulting in a lethal phenotype if disrupted, and retrieves a list of corresponding gene, mRNA, or enzyme targets. Using the dual link between MCSs and Elementary Flux Modes (EFMs), our logic programming-based tool aspefm was able to compute MCSs of any size from GSMMs in acceptable run times. The tool demonstrated better performance when computing large-sized MCSs than the mixed-integer linear programming methods. We applied the new MCSs methodology to a medically-relevant consortium model of two cross-feeding bacteria, Staphylococcus aureus and Pseudomonas aeruginosa. aspefm constraints were used to bias the computation of MCSs toward exchanged metabolites that could complement lethal phenotypes in individual species. We found that interspecies metabolite exchanges could play an essential role in rescuing single-species growth, for instance inosine could complement lethal reaction knock-outs in the purine synthesis, glycolysis, and pentose phosphate pathways of both bacteria. Finally, MCSs were used to derive a list of promising enzyme targets for consortium-level therapeutic applications that cannot be circumvented via interspecies metabolite exchange.

Integration of graph neural networks and genome-scale metabolic models for predicting gene essentiality

Genome-scale metabolic models are powerful tools for understanding cellular physiology. Flux balance analysis (FBA), in particular, is an optimization-based approach widely employed for predicting metabolic phenotypes. In model microbes such as Escherichia coli, FBA has been successful at predicting essential genes, i.e. those genes that impair survival when deleted. A central assumption in this approach is that both wild type and deletion strains optimize the same fitness objective. Although the optimality assumption may hold for the wild type metabolic network, deletion strains are not subject to the same evolutionary pressures and knock-out mutants may steer their metabolism to meet other objectives for survival. Here, we present FlowGAT, a hybrid FBA-machine learning strategy for predicting essentiality directly from wild type metabolic phenotypes. The approach is based on graph-structured representation of metabolic fluxes predicted by FBA, where nodes correspond to enzymatic reactions and edges quantify the propagation of metabolite mass flow between a reaction and its neighbours. We integrate this information into a graph neural network that can be trained on knock-out fitness assay data. Comparisons across different model architectures reveal that FlowGAT predictions for E. coli are close to those of FBA for several growth conditions. This suggests that essentiality of enzymatic genes can be predicted by exploiting the inherent network structure of metabolism. Our approach demonstrates the benefits of combining the mechanistic insights afforded by genome-scale models with the ability of deep learning to infer patterns from complex datasets.

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.

Napthyridine-derived compounds as promising inhibitors for Staphylococcus aureus CrtM: a primer for the discovery of potential anti-Staphylococcus aureus agents

The disease-free existence of humans is constantly under attack by a variety of infections caused by a variety of organisms including bacteria. Notable among the bacteria is Staphylococcus aureus which is an etiological organism for infections including impetigo, folliculitis, and furuncles. The response of the human immune system against this disease is often neutralized by the production of a pigment called Staphyloxanthin (STX) via a series of reactions mediated by several enzymes. Among these enzymes, dehydrosqualene synthase, also known as CrtM, has emerged as a viable drug target due to its role in mediating the first step of the pathway. Consequently, this study employs molecular modeling approaches including molecular docking, quantum mechanical calculations, and molecular dynamics (MD) simulations among others to investigate the potential of napthyridine derivatives to serve as inhibitors of the CrtM. The results of the study revealed the high binding affinities of the compounds for the target as demonstrated by their docking scores, while further subjection to screening pipeline aimed at determining their fitness for development into drugs revealed just one compound namely 6-[[1-[(2-fluorophenyl) methyl]triazol-4-yl]methoxy]-4-oxo-1H-1,5-naphthyridine-3-carboxylic acid as the compound with good drug-like, pharmacokinetics, and toxicity properties profiles. A 100 ns-long MD simulation of the complexes formed after molecular docking revealed the stable interaction of the compound with the target. Ultimately, this study can be a promising outlet to discover a weapon to fight against clinically resistant bacteria, however, further experimental studies are suggested to carry out in the wet lab, pre-clinical, and clinical levels.

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.

Comparison of Different Label-Free Raman Spectroscopy Approaches for the Discrimination of Clinical MRSA and MSSA Isolates

Methicillin-resistant Staphylococcus aureus (MRSA) is classified as one of the priority pathogens that threaten human health. Resistance detection with conventional microbiological methods takes several days, forcing physicians to administer empirical antimicrobial treatment that is not always appropriate. A need exists for a rapid, accurate, and cost-effective method that allows targeted antimicrobial therapy in limited time. In this pilot study, we investigate the efficacy of three different label-free Raman spectroscopic approaches to differentiate methicillin-resistant and -susceptible clinical isolates of S. aureus (MSSA). Single-cell analysis using 532 nm excitation was shown to be the most suitable approach since it captures information on the overall biochemical composition of the bacteria, predicting 87.5% of the strains correctly. UV resonance Raman microspectroscopy provided a balanced accuracy of 62.5% and was not sensitive enough in discriminating MRSA from MSSA. Excitation of 785 nm directly on the petri dish provided a balanced accuracy of 87.5%. However, the difference between the strains was derived from the dominant staphyloxanthin bands in the MRSA, a cell component not associated with the presence of methicillin resistance. This is the first step toward the development of label-free Raman spectroscopy for the discrimination of MRSA and MSSA using single-cell analysis with 532 nm excitation. IMPORTANCE Label-free Raman spectra capture the high chemical complexity of bacterial cells. Many different Raman approaches have been developed using different excitation wavelength and cell analysis methods. This study highlights the major importance of selecting the most suitable Raman approach, capable of providing spectral features that can be associated with the cell mechanism under investigation. It is shown that the approach of choice for differentiating MRSA from MSSA should be single-cell analysis with 532 nm excitation since it captures the difference in the overall biochemical composition. These results should be taken into consideration in future studies aiming for the development of label-free Raman spectroscopy as a clinical analytical tool for antimicrobial resistance determination.

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.

A New Perspective on the Antimicrobial Mechanism of Berberine Hydrochloride Against Staphylococcus aureus Revealed by Untargeted Metabolomic Studies

Berberine hydrochloride (BBR) is a natural product widely used in clinical medicine and animal production. It has a variety of antimicrobial effects, but its complex antimicrobial mechanism has not been clarified. This study aimed to discover the metabolic markers and gain a new perspective on the antibacterial mechanism of BBR. The effects of different inhibitory concentrations of BBR on the survival and growth of standard strain Staphylococcus aureus ATCC 25923 were analyzed by the bacteriostatic activity test. Differences in intracellular metabolites of S. aureus following 19 μg/ml BBR exposure for 1 h were investigated by combining non-targeted metabolomics techniques of gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS). The results showed that the minimum inhibitory concentration of BBR against S. aureus was 51 μg/ml. A total of 368 and 3,454 putative metabolites were identified by GC-MS and LC-MS analyses, respectively. Principal component analysis showed the separation of intracellular metabolite profiles between BBR-exposed samples and non-exposed controls. Pathway activity profiling analysis indicated a global inhibition of metabolisms by BBR exposure, while enhancement was also found in nucleic acid metabolism, amino sugar, and nucleotide sugar metabolism. Several metabolic markers were screened out mainly based on their variable importance of projection values. Two pyridine dicarboxylic acids were significantly downregulated, suggesting the reduction of stress resistance. The oxidized phospholipid (PHOOA-PE) was accumulated, while lipid antioxidant gamma-tocopherol was decreased, and farnesyl PP, the synthetic precursor of another antioxidant (staphyloxanthin), was decreased below the detection threshold. This evidence indicates that BBR reduced the antioxidant capacity of S. aureus. Accumulation of the precursors (UDP-GlcNAc, CDP-ribitol, and CDP-glycerol) and downregulation of the key metabolite D-Ala-D-Ala suggest the inhibition of cell wall synthesis, especially the peptidoglycan synthesis. Metabolites involved in the shikimate pathway (such as 3-dehydroshikimate) and downstream aromatic amino acid synthesis were disturbed. This study provides the first metabolomics information on the antibacterial mechanism of BBR against S. aureus. The key metabolic markers screened in this study suggest that the shikimate pathway, staphyloxanthin synthesis, and peptidoglycan biosynthesis are new directions for further study of BBR antibacterial mechanism in the future.

Nutritional Interactions between Bacterial Species Colonising the Human Nasal Cavity: Current Knowledge and Future Prospects

The human nasal microbiome can be a reservoir for several pathogens, including Staphylococcus aureus. However, certain harmless nasal commensals can interfere with pathogen colonisation, an ability that could be exploited to prevent infection. Although attractive as a prophylactic strategy, manipulation of nasal microbiomes to prevent pathogen colonisation requires a better understanding of the molecular mechanisms of interaction that occur between nasal commensals as well as between commensals and pathogens. Our knowledge concerning the mechanisms of pathogen exclusion and how stable community structures are established is patchy and incomplete. Nutrients are scarce in nasal cavities, which makes competitive or mutualistic traits in nutrient acquisition very likely. In this review, we focus on nutritional interactions that have been shown to or might occur between nasal microbiome members. We summarise concepts of nutrient release from complex host molecules and host cells as well as of intracommunity exchange of energy-rich fermentation products and siderophores. Finally, we discuss the potential of genome-based metabolic models to predict complex nutritional interactions between members of the nasal microbiome.

Rhea, the reaction knowledgebase in 2022

Rhea (https://www.rhea-db.org) is an expert-curated knowledgebase of biochemical reactions based on the chemical ontology ChEBI (Chemical Entities of Biological Interest) (https://www.ebi.ac.uk/chebi). In this paper, we describe a number of key developments in Rhea since our last report in the database issue of Nucleic Acids Research in 2019. These include improved reaction coverage in Rhea, the adoption of Rhea as the reference vocabulary for enzyme annotation in the UniProt knowledgebase UniProtKB (https://www.uniprot.org), the development of a new Rhea website, and the designation of Rhea as an ELIXIR Core Data Resource. We hope that these and other developments will enhance the utility of Rhea as a reference resource to study and engineer enzymes and the metabolic systems in which they function.