From random to directed motion: Understanding chemotaxis in E. Coli within a simplified model

The bacterium E.Coli swims in a zig-zag manner, in a series of straight runs and tumbles occurring alternately, with the run-durations dependent on the local spatial gradient of chemo-attractants/repellants. This enables the organism to move towards nutrient sources and move away from toxins. The signal transduction network of E.Coli has been well-characterized, and theoretical modeling has been used, with some success, in understanding its many remarkable features, including the near-perfect adaptation to spatially uniform stimulus. We study a reduced form of this network, with 3 methylation states for the receptor instead of 5. We derive an analytical form of the response function of the tumbling rate and use it to compute the drift velocity of the bacterium in the presence of a weak spatial attractant gradient.

💡 Research Summary

The paper investigates the chemotactic behavior of Escherichia coli by constructing a simplified version of the well‑studied signal‑transduction network that controls the run‑tumble motility pattern. In the canonical description, the chemoreceptor can exist in five methylation states, each influencing the activity of the kinase CheA and consequently the probability of a tumble. Experimental observations, however, suggest that only a subset of these states is functionally relevant under typical laboratory conditions. Motivated by this, the authors reduce the receptor model to three methylation levels (0, 1, 2), thereby decreasing the dimensionality of the system while preserving its essential adaptive features.

The authors first formulate the kinetic equations governing the occupancy of each methylation state, the ligand‑binding dynamics, and the transition rates between states (methylation rate k_m and demethylation rate k_d). They then linearize the system around a homogeneous steady state, assuming that the external attractant concentration varies only weakly in space (i.e., a shallow gradient). By applying Laplace transforms, they derive an analytical expression for the response function R(s) that relates a small perturbation in attractant concentration c(s) to the modulation of the tumble rate α(s):

  α(s) = R(s) c(s).

A crucial property of R(s) is that R(0)=0, guaranteeing “near‑perfect adaptation”: the tumble rate returns to its baseline value when the stimulus becomes spatially uniform. The time‑domain form of R(t) consists of an exponentially decaying component combined with a short‑lived overshoot, reflecting the rapid initial response followed by a slower adaptation phase.

Using the linear response framework, the authors compute the average drift velocity v_d of a bacterium moving in a weak attractant gradient G = ∇c. The drift emerges because the run duration is lengthened when the cell moves up the gradient and shortened when it moves down. The analytical result takes the familiar form

  v_d = χ G,

where the chemotactic sensitivity coefficient χ is expressed in terms of the model parameters: the basal swimming speed v₀, the mean run time τ_run, the methylation/demethylation rates, and the derivative of the response function at zero frequency R′(0). Explicitly,

  χ = (v₀ τ_run² R′(0)) /