Whole-Colony Modeling of Escherichia coli

Abstract

AbstractBacterial behavior is the outcome of both molecular mechanisms within each cell and interactions between cells in the context of their environment. Whereas whole-cell models simulate a single cell’s behavior using molecular mechanisms, agent-based models simulate many agents independently acting and interacting to generate complex collective phenomena. To synthesize agent-based and whole-cell modeling, we used a novel model integration software, called Vivarium, to construct an agent-based model of E. coli colonies where each agent is represented by a current source code snapshot from the E. coli Whole-Cell Modeling Project and interacts with other cells in a shared spatial environment. The result is the first “whole-colony” computational model that mechanistically links expression of individual proteins to a population-level phenotype. Simulated colonies exhibit heterogeneous effects on their environments, heterogeneous gene expression, and media-dependent growth. Extending the cellular model with mechanisms of antibiotic susceptibility and resistance, our model also suggested that variation in the expression level of the betalactamase AmpC, and not of the multi-drug efflux pump AcrAB-TolC, was the key mechanistic driver of survival in the presence of nitrocefin. We see this as a significant step forward in the creation of more comprehensive multi-scale models, and it broadens the range of phenomena that can be modeled in mechanistic terms.Author summaryThis work combines several models of molecular and physical processes that impact the physiology and behavior of the common microbe Escherichia coli into a multiscale model. Colonies comprised of multiple individual cells are simulated as they grow and divide—each with complex internal mechanisms, and with physical interactions and molecular diffusion in their environments. The integrative modeling methodology supports the addition of new submodels. The flexibility of this methodology is demonstrated by adding models of antibiotic resistance and simulating the colony’s response to antibiotic treatment.