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.

MetaBiome: A Multiscale Model Integrating Agent-Based Modeling and Metabolic Networks Reveals Spatial Regulation in Mucosal Microbial Communities

Metabolic interactions shape emergent biofilm structures in a conceptual model of gut mucosal bacterial communities

Metabolic Modeling of Microbial Community Interactions for Health, Environmental and Biotechnological Applications

Proteomics and Metaproteomics Add Functional, Taxonomic and Biomass Dimensions to Modeling the Ecosystem at the Mucosal-luminal Interface

Emergence of microbial diversity due to cross-feeding interactions in a spatial model of gut microbial metabolism

GEM-based computational modeling for exploring metabolic interactions in a microbial community

Virtual Colon: Spatiotemporal modelling of metabolic interactions in a computational colonic environment

Community‐scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

A Top-Down Systems Biology View of Microbiome-Mammalian Metabolic Interactions in A Mouse Model

Metabolic model-based ecological modeling for probiotic design

Agent Based Modeling of Human Gut Microbiome Interactions and Perturbations

Modeling microbial cross-feeding at intermediate scale portrays community dynamics and species coexistence

Metabolic Network Modeling of Microbial Interactions in Natural and Engineered Environmental Systems

Spatio-temporal modeling of the crowding conditions and metabolic variability in microbial communities

Microbial community-scale metabolic modeling predicts personalized short chain fatty acid production profiles in the human gut

Probing interspecies metabolic interactions within a synthetic binary microbiome using genome-scale modeling

Metabolic modeling of species interaction in the human microbiome elucidates community-level assembly rules

Metabolic Modeling of Human Gut Microbiota on a Genome Scale: An Overview

Microbiome modeling: a beginner's guide

Inferring microbial interactions with their environment from genomic and metagenomic data

Spatial transcriptome-guided multi-scale framework connects P . aeruginosa metabolic states to oxidative stress biofilm microenvironment