Modeling microbial chemotaxis in a diffusion gradient chamber

The diffusion gradient chamber (DGC) has proven to be a useful experimental tool for studying population-level microbial growth and chemotaxis. A mathematical model capable of reproducing the population-level patterns formed as a result of cellular growth and chemotaxis in the DGC has been developed. The model consists of coupled partial differential balance equations for cells, chemoattractants, and a nutrient, which are solved simultaneously by the alternating direction implicit method. Modeling simulation results were compared with population-level migration patterns of Escherichia coli growing on glycerol and responding to a gradient of the chemoattractant aspartate for two different initial conditions. To accurately reproduce the experimental results, a second chemoattractant equation was necessary. The second chemoattractant has been identified as oxygen by directly measuring oxygen gradients in the DGC. Important trends observed experimentally and reproduced by the model include the formation of a chemotactic wave, a reduction in the wave velocity as it encounters higher chemoattractant concentrations, and chemotaxis in response to two different chemoattractants simultaneously. The model was also used to study the relative magnitude of cell fluxes due to random motility and chemotaxis, and the suppression of chemotaxis due to receptor saturation.