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.

Multiscale flux-based modeling of biofilm communities

Multiscale Computational Modeling of Biofilm

Biophysical metabolic modeling of complex bacterial colony morphology

Anomalous transparency of water-air interface for low-frequency sound

Rheumatoid arthritis--a disease of dysregulated glycosylation?

Validity of Commercial Activity Trackers in Children With Congenital Heart Disease.

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

Flux sampling in genome-scale metabolic modeling of microbial communities

Pore‐Scale Coupling of Flow, Biofilm Growth, and Nutrient Transport: A Microcontinuum Approach

Microbial communities and their interactions in biofilm systems: an overview

Modeling multispecies biofilms including new bacterial species invasion

Unification of cell-scale metabolic activity with biofilm behavior by integration of advanced flow and reactive-transport modeling and microfluidic experiments

A multiscale hybrid modelling methodology for cell cultures enabled by enzyme-constrained dynamic metabolic flux analysis under uncertainty

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

Metabolic Flux Modeling in Marine Ecosystems

Coupled metabolic‐hydrodynamic modeling enabling rational scale‐up of industrial bioprocesses

Is it selfish to be filamentous in biofilms? Individual-based modeling links microbial growth strategies with morphology using the new and modular iDynoMiCS 2.0

Understanding photosynthetic biofilm productivity and structure through 2D simulation

Microfluidic and mathematical modeling of aquatic microbial communities

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

Biofilm Growth in Porous Media Well Approximated by Fractal Multirate Mass Transfer With Advective‐Diffusive Solute Exchange