Equation-Based Integration of Flux Balance Analysis with Diffusion for Spatio-Temporal Simulation of Microbial Communities

Abstract

Spatio-temporal interactions shape microbial community dynamics. Metabolism, through competition and cross-feeding, is a foundational mechanism of these interactions. Flux balance analysis enables efficient simulation of steady-state metabolism. Integrating these simulations through time, using dynamic flux balance analysis, provides temporal predictions of growth and metabolism. Incorporating spatial context, through partial differential equations, enables spatio-temporal simulation of microbial communities. In this chapter, we step through this sequential process, moving from steady-state, to temporal, to spatio-temporal simulation of microbial community metabolism. We provide an illustrative example using the modeling software COMETS (Computation of Microbial Ecosystems in Time and Space) to simulate interacting bacterial colonies of Bifidobacterium longum subsp. infantis and Anaerobutyricum hallii (previously Eubacterium hallii). Within this simulation, both competition and cross-feeding influenced the production of butyrate leading to an intermediate optimal interaction distance for metabolite production. We outline each step and provide open-source code such that this simulation can serve as a template for future spatio-temporal simulations of microbial community metabolism.