Genome-scale metabolic model of Staphylococcus epidermidis ATCC 12228 matches in vitro conditions
1Computational Systems Biology of Infections and Antimicrobial-Resistant Pathogens, Institute for Bioinformatics and Medical Informatics (IBMI), Eberhard Karl University of Tübingen, 72076 Tübingen, DE
2Department of Computer Science, Eberhard Karl University of Tübingen, 72076 Tübingen, DE
3Cluster of Excellence ‘Controlling Microbes to Fight Infections’, Eberhard Karl University of Tübingen, DE
4German Center for Infection Research (DZIF), partner site Tübingen, DE
5Institute for Pharmaceutical Microbiology, University of Bonn, University Hospital Bonn, Bonn, DE
6German Center for Infection Research (DZIF), partner site Bonn-Cologne, Bonn, DE
7Structural and Computational Biology Unit, European Molecular Biology Laboratory (EMBL), 69117 Heidelberg, DE
*Correspondence: nantia.leonidou@uni-tuebingen.deABSTRACT
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 lifestyle and 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 MEMOTE test suite yields a high score, highlighting the model’s high semantic quality. Following the 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. This metabolic model advances our understanding of S. epidermidis as a commensal and potential probiotic and enhances insights into its opportunistic pathogenicity against other microorganisms.
Article notes
Competing Interest Statement
The authors have declared no competing interest.
Introduction
A prevalent constituent of the human skin flora is the coagulase-negative commensal Staphylococcus epidermidis1, 2. This Gram-positive coccus predominantly inhabits the skin and mucosal membranes in areas such as the axillae, head, legs, arms, and nares. S. epidermidis plays a crucial role in maintaining a balanced microbiome within the human nasal cavity, where harmful pathogens like Staphylococcus aureus commonly establish colonization. There is ongoing discourse regarding whether S. epidermidis, through competition in nutritionally scarce environments like the human nose, may exhibit probiotic effects against formidable pathogens such as S. aureus2, 3. Nevertheless, S. epidermidis is recognized as a significant causative agent of nosocomial infections under specific conditions4. Notably, S. epidermidis stands out as the primary source of infections associated with indwelling medical devices, including intravascular catheters and implants such as prosthetic joints1, 5, 6. The high occurrence of these nosocomial infections is attributed to S. epidermidis’s ubiquitous presence on the human skin, increasing the likelihood of contamination during the insertion of medical devices7. Upon infection, S. epidermidis strains are capable of forming biofilms that shield them from antibiotics and host defense mechanisms, rendering S. epidermidis infections resistant and challenging to eliminate1, 7. Often, removing the foreign material becomes necessary to combat the infection effectively. While S. epidermidis infections seldom lead to life-threatening conditions, their impact on patients and the public health system is substantial. In the United States alone, the annual economic burden of S. epidermidis vascular catheter-related bloodstream infections is estimated to be around $2 billion1, 6. Besides biofilm formation, also other specific molecular determinants contribute to the pathogenicity of this particular pathogen, enabling immune evasion. Therefore, there is an urgent need for a more comprehensive understanding of S. epidermidis and its opportunistic characteristics to identify novel therapeutic strategies1, 7.
One way to better understand an organism’s lifestyle and capabilities is the reconstruction and analysis of genome-scale metabolic models (GEMs). These models rely on the annotated genome sequence of the organism in question. Specifically, genes encoding proteins with metabolic significance are allocated to their respective reactions through gene-protein-reaction associations (GPRs). Within the resulting network, biochemical reactions establish connections between metabolites, with enzymatic activities guided by genes associated with these reactions. Such models enable the comprehension of an organism’s metabolism at a systemic level. Díaz Calvo et al. reconstructed the metabolic network of RP62A, a slime-producing and methicillin-resistant biofilm isolate8. However, the resulted model is available only upon request. Figure 1 summarizes the computational and experimental approach of this article. This work introduces iSep23, the first manually curated and experimentally validated GEM of S. epidermidis ATCC 12228. The model comprises 1,415 reactions, 1,051 metabolites, and 705 genes and is freely available from BioModels Database9 with the accession identifier MODEL2012220002. Moreover, it aligns with current community standards10, 11, 12 and modeling guidelines13, 14. Semantic benchmarking was conducted utilizing the MEMOTE genome-scale metabolic model test suite15. Consequently, iSep23 upholds the Findable, Accessible, Interoperable, and Reusable (FAIR) data principles16, rendering it a valuable resource for subsequent research17, 18. To assess the predictive capacity of the model, growth simulations in various media were compared against laboratory experiments. The model’s predictions regarding the utilization of diverse carbon sources were cross-referenced with experimental findings. Altogether, our model establishes a foundation for improved comprehension of the organism’s phenotypes and behavior under different nutritional conditions.
Results
Properties of the constructed GEM
The initial CarveMe draft comprised 1,295 reactions, 933 metabolites, and 722 genes, yielding a Metabolic Model Testing (MEMOTE) 15 score of 36 %. Subsequent manual refinement involved the addition of 120 reactions, 118 metabolites, and 63 genes, as illustrated in Figure 2, resulting in an overall MEMOTE score of 88 %. The 63 mass- and charge-imbalanced reactions were reduced to one mass-imbalanced and nine charge-imbalanced reactions, resulting in a MEMOTE mass balance score of 99.7 % and a charge balance score of 99.3 %. Based on literature evidence, we corrected the directionality of 34 enzymatic reactions in the model to ensure proper constraints during model simulations. Moreover, the final metabolic network does not include infeasible energy generating cycle (EGC) that could inflate the simulation results (see Materials and Methods). We annotated the model instances with cross-references to various databases and additional information to increase the model’s interoperability and re-usability. The reaction annotations are divided into three different biological qualifier types:
- The cross-references to the nine databases are stored under the biological qualifier type
BQB_IS. - The ECO terms are stored under
BQB_IS_DESCRIED_BY. - Pathways associated with a reaction are saved with the biological qualifier type
BQB_OCCURS_IN.
The metabolites and genes were annotated with twelve and three external databases, respectively, using the biological qualifier type BQB_IS (Table 1). The inclusion of ECO terms ensures a comprehensive understanding of evidence and assertion methodologies36, thereby facilitating robust quality control measures and evidence queries. The ECO term with the lowest evidence level is ECO:0000001, coding for inference from background scientific knowledge (Figure 2). This term was ascribed to 30.2 % of the biochemical reactions within the network. Notably, this percentage encompasses pseudo-reactions, such as exchanges, sinks, demands, and the biomass function. Within the group of 431 reactions associated with this ECO term, 170 pertained to pseudo reactions. The ECO term ECO:0000251 denotes similarity evidence used in automatic assertion and was assigned to 28.5 % of all reactions. Moreover, the terms ECO:0(computational inference used in automatic assertion) and ECO:0000044 (sequence similarity evidence) annotated 9.3 % and 31.9 % of all reactions, respectively. A minimal fraction (0.1 %) of reactions exhibits protein assay evidence, identified by the ECO:0000039 term. Additionally, the SBOannotator was utilized to annotate the model with precise and descriptive SBO terms19 (Figure 2). Totally, 25 terms were incorporated describing classes of bio-chemical reactions and further model elements.
The final curated metabolic model was stored as a Systems Biology Markup Language (SBML) Level 337 file. This format version supports the integration of various plugins, such as the fbc package38 and the groups package39, which are both enabled in iSep23. The groups package facilitates the incorporation of additional information without impacting the mathematical interpretation of the model. We defined all pathways and subsystems identified from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database24 as an individual group and added corresponding reactions as members. Overall, we added 99 distinct groups to the model that facilitate pathway-related analysis.
Discussion
Here, we present a manually curated GEM of S. epidermidis ATCC 12228, iSep23. Literature-based corrections and meticulous manual curation ensured accurate representation of enzymatic reaction directions, essential for precise constraints during simulations. Overall, our model aligns with experimental data and offers a comprehensive platform for exploring S. epidermidis’s metabolic capabilities and behavior under diverse conditions. The inconsitency between the in silco and in vitro results reagrding the AAM-in the presence of glucose could be attributed to factors beyond the metabolic scope. For instance, non-metabolic factors could be regulatory mechanisms and Post-translational modifications. The observed discrepancy suggests a need for a more detailed understanding of the regulatory and metabolic factors influencing S. epidermidis growth in AAM-. Further experimental validation and exploration of regulatory mechanisms are crucial for resolving the observed differences between in silico predictions and experimental outcomes.
All in all, the refined network serves as a powerful tool for exploring S. epidermidis’s metabolic capabilities and behavior under diverse conditions. Future perspectives involve leveraging the model for targeted studies, such as investigating metabolic pathways, assessing the impact of genetic modifications, and exploring potential drug targets. The model’s compatibility with the fbc and groups packages in the SBML Level 3 Version 112 format enhances its flexibility, enabling the integration of additional plugins for more intricate analyses. Including 99 distinct groups representing pathways and subsystems from the KEGG database provides a foundation for comprehensive pathway-related analyses. Altogether, iSep23 aligns with experimental data and lays the groundwork for future investigations into the bacterium’s metabolism. Its accuracy, comprehensibility, and flexibility make it a valuable resource for advancing our understanding of microbial physiology and metabolic engineering applications.
Materials and Methods
Reconstructing the draft model of S. epidermidis
The reconstruction of the GEM is based on protocols described in previous studies45, 46. The fast and automated reconstruction tool CarveMe47 curates genome-scale metabolic models of microbial species and communities47. During the initial curation phase, a universal model was systematically compared to the annotated genome sequence of the species of interest, facilitating the construction of individual single-species metabolic models. In this study, we utilized CarveMe version 1.2.2 and the annotated genome sequence of S. epidermidis ATCC 12228 with the RefSeq48 accession ID NC_004461.1 that covers the bacterial chromosome. Throughout the drafting process and subsequent model iterations, rigorous monitoring and benchmarking were conducted using MEMOTE 15. MEMOTE performs standardized metabolic tests across four key domains: annotation, basic tests, biomass reaction, and stoichiometry. The results are stored in a comprehensive report that includes the model’s overall performance assessed by a metric called MEMOTE score (denoted as a percentage with 100 %). A higher MEMOTE score correlates with enhanced annotation quality, greater consistency, and formal correctness of the model in SBML49 format. To refine the initial model automatically, the ModelPolisher50 was employed in a preliminary step. Leveraging the Biochemical, Genetical, and Genomical (BiGG) Models database20 identifiers of the model instances, the ModelPolisher systematically accessed the BiGG Models database, assimilating all available information for these instances into the network as annotations.
Evaluation and validation of growth capabilities
Different growth media
The growth behavior of S. epidermidis was assessed in three distinct synthetic minimal media initially formulated for investigating the metabolic requirements of S. aureus. These are the: (i) SMM40, (ii) AAM41, and (iii) AAM-42; a modified version of the AAM medium. The concentrations of the various components served as lower bounds for the corresponding exchange reactions of metabolites, as detailed in Table 2. In addition to the already provided salts and ions, we added minimal traces of zinc (EX_zn2_e), cobalt (EX_cobalt2_e), and copper (EX_cu2_e) to the simulated medium to enable growth. The lower bound of these reactions was set to −0.0001 mmol/(gDW · h). Oxygen availability was defined by setting the lower bound of the exchange reaction to −20 mmol/(gDW · h). The initial formulation of the three media involved the use of nicotinic acid. However, as nicotinic acid was substituted with nicotinamide in laboratory experiments, our simulated media also incorporated nicotinamide. In addition to the three minimal media, we tested S. epidermidis’s growth on the LB47. The lower bounds of the compounds’ exchange reactions listed in the LB were set to −10 mmol/(gDW · h). All in silico simulations were evaluated with and without d-glucose as a carbon source.
Different carbon sources
Twelve different sugars were tested for their potential role as a carbon source: d-glucose, d-arabinose, maltose, lactose, raffinose, d-sucrose, trehalose, d-xylose, d-cellobiose, fructose, mannose, and d-ribose. For the growth simulations in different carbon sources, we used the SMM with nicotinamide instead of nicotinic acid as a basis (see Table 3). The concentrations reported in the medium were established as lower bounds for the simulation. The concentrations of the listed carbon sources were calculated to be equivalent in carbon content to the initial 5 g/L of glucose used in the defined SMM.
Laboratory validation
Media preparation
The minimal media AAM, AAM-, and SMM were prepared as carbon-source free base media following the methods provided by Machado et al. after omitting glucose as the default carbon source40. The carbohydrates to replace glucose as alternative carbon sources were dissolved in their respective base medium, and the resulting media were sterile filtered. Carbohydrates were obtained from Carl Roth (d-arabinose, d-glucose, trehalose, lactose, sucrose, raffinose), EMD-Millipore (fructose), Fluka (maltose, d-cellobiose), and Sigma Aldrich (mannose, d-ribose, d-xylose) in purity grades of ≥ 98 %. LB was prepared following the standard formulation of 10 g/L tryptone (MP Biomedicals), 10 g/L sodium chloride (Carl Roth), 5 g/L yeast extract (Carl Roth), and 5 g/L glucose when required.
Growth experiments
Cultures of S. epidermidis ATCC 12228 were initiated by inoculating overnight precultures in LB at 37 °C. Subsequently, primary cultures in LB were established from them and allowed to grow to an optical density (OD) at 600 nm (OD600 nm) of 0.5. Cell harvesting was achieved through centrifugation and two washes with the carbon-source-free medium. The cells were then resuspended to an OD600 nm of 0.05 in media containing the respective carbon source. Growth was assessed by determination of the OD after a 24 h-incubation at 37 °C. Growth experiments were performed in at least three biological replicates in a 96-well plate format. OD measurements were performed with a Tecan Spark microplate reader.
Data availability
Acknowledgments
This work was funded by the Deutsche Forschungsgemein-schaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2124 – 390838134 and supported by the Cluster of Excellence ‘Controlling Microbes to Fight Infections’ (CMFI). F.G. and A.D. is supported by the German Center for Infection Research (DZIF, doi: 10.13039/100009139) within the Deutsche Zentren der Gesundheitsforschung (BMBF-DZG, German Centers for Health Research of the Federal Ministry of Education and Research (BMBF)), grants № 8020708703 and № 8016708710. The authors acknowledge the support by the Open Access Publishing Fund of the University of Tübingen (https://uni-tuebingen.de/en/216529).
Competing interests
The authors declare no conflict of interest.
List of Abbreviations
- API
- Application Programming transfer Interface
- BiGG
- Biochemical, Genetical, and Genomical
- BMBF
- Federal Ministry of Education and Research (Bundesministerium für Bildung und Forschung)
- BMBF-DZG
- Deutsche Zentren der Gesundheitsforschung
- CMFI
- Controlling Microbes to Fight Infections
- CV
- controlled vocabulary
- DFG
- Deutsche Forschungsgemeinschaft
- DZIF
- German Center for Infection Research
- ECO
- Evidence and Conclusion Ontology
- EGC
- energy generating cycle
- EMBL
- European Molecular Biology Laboratory
- FAIR
- Findable, Accessible, Interoperable, and Reusable
- GEM
- genome-scale metabolic model
- GPR
- gene-protein-reaction association
- KEGG
- Kyoto Encyclopedia of Genes and Genomes
- LB
- lysogeny broth
- MEMOTE
- Metabolic Model Testing
- NCBI
- National Centre for Biotechnology Information
- OD
- optical density
- REST
- Representational State Transfer
- SBML
- Systems Biology Markup Language
- SBO
- Systems Biology Ontology
- SMM
- synthetic minimal medium