An agent-based model of metabolic signaling oscillations in Bacillus subtilis biofilms

Microbes of nearly every species can form biofilms, communities of cells bound together by a self-produced matrix. It is not understood how variation at the cellular level impacts putatively beneficial, colony-level behaviors, such as cell-to-cell signaling. Here we investigate this problem with an agent-based computational model of metabolically driven electrochemical signaling in Bacillus subtilis biofilms. In this process, glutamate-starved interior cells release potassium, triggering a depolarizing wave that spreads to exterior cells and limits their glutamate uptake. More nutrients diffuse to the interior, temporarily reducing glutamate stress and leading to oscillations. In our model, each cell has a membrane potential coupled to metabolism. As a simulated biofilm grows, collective membrane potential oscillations arise spontaneously as cells deplete nutrients and trigger potassium release, reproducing experimental observations. We further validate our model by comparing spatial signaling patterns and cellular signaling rates with those observed experimentally. By oscillating external glutamate and potassium, we find that biofilms synchronize to external potassium, providing a potential mechanism for previously observed oscillation synchronization. By tracking cellular glutamate concentrations, we find that oscillations evenly distribute nutrients in space: non-oscillating biofilms have an external layer of well-fed cells surrounding a starved core, whereas oscillating biofilms exhibit a relatively uniform distribution of glutamate. Our work shows the potential of agent-based models to connect cellular properties to collective phenomena and facilitates studies of how inheritance of cellular traits can affect the evolution of group behaviors.