Biophysical metabolic modeling of complex bacterial colony morphology

Microbial colony growth is shaped by the physics of biomass propagation and nutrient diffusion, and by the metabolic reactions that organisms activate as a function of the surrounding environment. While microbial colonies have been explored using minimal models of growth and motility, full integration of biomass propagation and metabolism is still lacking. Here, building upon our framework for Computation of Microbial Ecosystems in Time and Space (COMETS), we combine dynamic flux balance modeling of metabolism with collective biomass propagation and demographic fluctuations to provide nuanced simulations of E. coli colonies. Simulations produced realistic colony morphology, consistent with our experiments. They characterize the transition between smooth and furcated colonies and the decay of genetic diversity. Furthermore, we demonstrate that under certain conditions, biomass can accumulate along "metabolic rings" that are reminiscent of coffee-stain rings, but have a completely different origin. Our approach is a key step towards predictive microbial ecosystems modeling.