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.

A Domain Decomposition Method for Solution of a PDE-Constrained Generalized Nash Equilibrium Model of Biofilm Community Metabolism

Inverse problems for a model of biofilm growth

Stability analysis of hyperbolic-parabolic free boundary problems modelling biofilms

Spatiotemporal modeling of microbial metabolism

Game-theoretic Flux Balance Analysis Model for Predicting Stable Community Composition

Mathematical Modeling of Heavy Metal Biosorption in Multispecies Biofilms

Analysis of a finite-volume scheme for a single-species biofilm model

A chemostat model with variable dilution rate due to biofilm growth

Stochastic Modeling of Biofilm Formation with Bacterial Quorum Sensing

Continuum and discrete approach in modeling biofilm development and structure: a review

Designing Metabolic Division of Labor in Microbial Communities

Invasion moving boundary problem for a biofilm reactor model

Biofilm.jl: a fast solver for one-dimensional biofilm chemistry and ecology

Immersed boundary approach to biofilm spread on surfaces

Free boundary problem for the role of planktonic cells in biofilm formation and development

A Mathematical Model of Diel Activity and Long Time Survival in Phototrophic Mixed-Species Subaerial Biofilms

Metabolic Co-Dependence Gives Rise to Collective Oscillations Within Biofilms

Mathematical analysis of two microbial species competing for two complementary resources with internal storages and different removal , dilution rates

Mathematical treatment of PDE model describing chemotactic E. coli colonies

Mathematical Analysis Of Two Microbial Species Competing For Two Complementary Resources With Internal Storage And Different Removal Rates

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