Theoretical results for chemotactic response and drift of E. coli in a weak attractant gradient

The bacterium Escherichia coli (E. coli) moves in its natural environment in a series of straight runs, interrupted by tumbles which cause change of direction. It performs chemotaxis towards chemo-attractants by extending the duration of runs in the direction of the source. When there is a spatial gradient in the attractant concentration, this bias produces a drift velocity directed towards its source, whereas in a uniform concentration, E.coli adapts, almost perfectly in case of methyl aspartate. Recently, microfluidic experiments have measured the drift velocity of E.coli in precisely controlled attractant gradients, but no general theoretical expression for the same exists. With this motivation, we study an analytically soluble model here, based on the Barkai-Leibler model, originally introduced to explain the perfect adaptation. Rigorous mathematical expressions are obtained for the chemotactic response function and the drift velocity in the limit of weak gradients and under the assumption of completely random tumbles. The theoretical predictions compare favorably with experimental results, especially at high concentrations. We further show that the signal transduction network weakens the dependence of the drift on concentration, thus enhancing the range of sensitivity.

💡 Research Summary

The paper addresses a long‑standing gap in the quantitative theory of bacterial chemotaxis: a general expression for the drift velocity of Escherichia coli moving in a weak attractant gradient. Building on the well‑known Barkai‑Leibler (BL) model, which successfully explains perfect adaptation to uniform concentrations of methyl‑aspartate, the authors develop an analytically tractable framework that couples the BL signaling network to the physical run‑and‑tumble motility pattern. The key assumptions are (i) a spatially linear, weak gradient (∇c → 0) so that linear response theory applies, and (ii) completely random tumbles, meaning that after each tumble the swimming direction is uniformly redistributed.

Starting from the BL differential equations for receptor activity, the authors perform a Laplace‑transform analysis to obtain the chemotactic response function χ(t). Remarkably, χ(t) emerges as a sum of two exponential terms, reflecting the fast ligand‑binding dynamics and the slower methylation/demethylation adaptation processes. This response function is then convolved with the instantaneous swimming velocity to calculate the modulation of run duration as a function of the angle between the swimming direction and the gradient. By averaging over all possible directions, they derive a closed‑form expression for the mean drift velocity:

 v_d = κ(c₀, λ, k_m, …) · ∇c

where κ is a composite coefficient that depends on the background attractant concentration c₀, the tumble frequency λ, the methylation rate k_m, and other biochemical parameters of the BL network. Importantly, κ decreases at high c₀, indicating a saturation effect that stems from the network’s near‑perfect adaptation; this attenuation broadens the concentration range over which the cell remains responsive.

The theoretical predictions are compared with recent microfluidic experiments that measured E. coli drift speeds under precisely controlled gradients of methyl‑aspartate. In the low‑gradient regime the model slightly underestimates the drift, likely because real tumbles are not perfectly random. However, at higher concentrations the predicted velocities match the experimental data almost exactly, confirming that the BL‑based linear response captures the essential physics of chemotactic drift. The analysis also quantifies how variations in λ or k_m would alter κ, offering a direct route for experimental manipulation or synthetic biology applications.

The authors discuss several limitations. The assumption of completely random tumbles neglects rotational diffusion and inertia that have been observed in high‑resolution tracking studies. The linear‑response framework cannot address steep or time‑varying gradients, and the model does not incorporate cell‑cell interactions, viscosity changes, or temperature effects. Despite these constraints, the work provides the first general analytical formula linking the biochemical adaptation network to macroscopic drift, bridging a gap between molecular signaling and population‑level behavior.

In conclusion, the study demonstrates that a BL‑based linear response theory yields a robust, experimentally validated expression for E. coli chemotactic drift in weak attractant gradients. The signal transduction network’s adaptation dynamics attenuate the concentration dependence of drift, thereby extending the sensory range. Future extensions to non‑linear gradients, partially biased tumbles, and collective dynamics are suggested as promising directions to deepen our understanding of bacterial navigation in complex environments.