Spatial transcriptome-guided multi-scale framework connects P . aeruginosa metabolic states to oxidative stress biofilm microenvironment
Tracy J. Kuper,Mohammad Mazharul Islam,Shayn M. Peirce-Cottler,Jason A. Papin,Roseanne M Ford
DOI: https://doi.org/10.1371/journal.pcbi.1012031
2024-04-27
PLoS Computational Biology
Abstract:With the generation of spatially resolved transcriptomics of microbial biofilms, computational tools can be used to integrate this data to elucidate the multi-scale mechanisms controlling heterogeneous biofilm metabolism. This work presents a Multi-scale model of Metabolism In Cellular Systems (MiMICS) which is a computational framework that couples a genome-scale metabolic network reconstruction (GENRE) with Hybrid Automata Library (HAL), an existing agent-based model and reaction-diffusion model platform. A key feature of MiMICS is the ability to incorporate multiple -omics-guided metabolic models, which can represent unique metabolic states that yield different metabolic parameter values passed to the extracellular models. We used MiMICS to simulate Pseudomonas aeruginosa regulation of denitrification and oxidative stress metabolism in hypoxic and nitric oxide (NO) biofilm microenvironments. Integration of P . aeruginosa PA14 biofilm spatial transcriptomic data into a P . aeruginosa PA14 GENRE generated four PA14 metabolic model states that were input into MiMICS. Characteristic of aerobic, denitrification, and oxidative stress metabolism, the four metabolic model states predicted different oxygen, nitrate, and NO exchange fluxes that were passed as inputs to update the agent's local metabolite concentrations in the extracellular reaction-diffusion model. Individual bacterial agents chose a PA14 metabolic model state based on a combination of stochastic rules, and agents sensing local oxygen and NO. Transcriptome-guided MiMICS predictions suggested microscale denitrification and oxidative stress metabolic heterogeneity emerged due to local variability in the NO biofilm microenvironment. MiMICS accurately predicted the biofilm's spatial relationships between denitrification, oxidative stress, and central carbon metabolism. As simulated cells responded to extracellular NO, MiMICS revealed dynamics of cell populations heterogeneously upregulating reactions in the denitrification pathway, which may function to maintain NO levels within non-toxic ranges. We demonstrated that MiMICS is a valuable computational tool to incorporate multiple -omics-guided metabolic models to mechanistically map heterogeneous microbial metabolic states to the biofilm microenvironment. Microbes secrete and respond to environmental metabolite signals, resulting in the spatial organization of heterogeneous physiological metabolic states within infectious microbial communities. Despite experimental advances, it is difficult to measure the connected, dynamic processes that control microbial community organization across multiple spatial and time scales. Thus, we developed an extendable multi-scale computational framework to simulate metabolic processes spanning intracellular to extracellular scales. We used this framework to simulate a Pseudomonas aeruginosa , a bacterium that causes lung infections, 3D biofilm community containing tens of thousands of single-cells. This is the first framework to integrate spatially resolved gene expression data. This data integration was advantageous to simulate cells regulating anaerobic, toxic metabolite byproduct secretion, and antioxidant metabolic states in response to a heterogeneous oxygen and toxic byproduct biofilm microenvironment. These cellular metabolic states were predicted to co-exist in biofilm spatial niches due to microscale variations in extracellular oxygen and toxic metabolite signals. Difficult to measure in experiments, this framework revealed multi-scale mechanisms underlying the emergent dynamics of each metabolic state, which may help to inform P . aeruginosa biofilm treatment strategies. We believe this multi-scale framework is valuable to integrate microbial data to uncover the multi-scale mechanisms regulating heterogenous metabolic processes in microbial communities.
biochemical research methods,mathematical & computational biology