Run-and-tumble motion

Biological examples

Run-and-tumble motion is found in many peritrichous bacteria, including E. coli, Salmonella typhimurium, and Bacillus subtilis. It has also been observed in the alga Chlamydomonas reinhardtii and the cyanobacterium Synechocystis.

Directed motility (taxis)

Genetically diverse groups of microorganisms rely upon directed motility (taxis), such as chemotaxis or phototaxis, to optimally navigate through complex environments or colonise host tissues. In the model organisms Escherichia coli and Salmonella, bacteria swim in a random pattern produced by alternating counterclockwise (CCW) and clockwise (CW) flagellar rotation. Chemoreceptors detect attractants or repellents and stimulate responses through a signalling cascade that controls the direction of the flagellar motor. This can result in a chemotaxis, where attractant gradients extend the length of time flagellar motors rotate CCW, resulting in more smooth swimming in a favourable direction, while repellents cause an increase of CW rotations, resulting in more tumbling and changes in direction. The cyanobacterium Synechocystis uses run-and-tumbling in a manner which can result in phototaxis.

Escherichia coli

An archetype of bacterial swimming is represented by the well-studied model organism Escherichia coli. With its peritrichous flagellation, E. coli performs a run-and-tumble swimming pattern, as shown in the diagrams below. Counterclockwise rotation of the flagellar motors leads to flagellar bundle formation that pushes the cell in a forward run, parallel to the long axis of the cell. Clockwise rotation disassembles the bundle and the cell rotates randomly (tumbling). After the tumbling event, straight swimming is recovered in a new direction. That is, counterclockwise rotation results in steady motion and clockwise rotation in tumbling; counterclockwise rotation in a given direction is maintained longer in the presence of molecules of interest (like sugars or aminoacids).

In a uniform medium, run-and-tumble trajectories appear as a sequence of nearly straight segments interspersed by erratic reorientation events, during which the bacterium remains stationary. The straight segments correspond to the runs, and the reorientation events correspond to the tumbles. Because they exist at low Reynolds number, bacteria starting at rest quickly reach a fixed terminal velocity, so the runs can be approximated as constant velocity motion. The deviation of real-world runs from straight lines is usually attributed to rotational diffusion, which causes small fluctuations in the orientation over the course of a run.

In contrast with the more gradual effect of rotational diffusion, the change in orientation (turn angle) during a tumble is large; for an isolated E. Coli in a uniform aqueous medium, the mean turn angle is about 70 degrees, with a relatively broad distribution. In more complex environments, the tumbling distribution and run duration may depend on the agent's local environment, which allows for goal-oriented navigation (taxis). For example, a tumbling distribution that depends on a chemical gradient can guide bacteria toward a food source or away from a repellant, a behavior referred to as chemotaxis. Tumbles are typically faster than runs: tumbling events of E. Coli last about 0.1 seconds, compared to ~ 1 second for a run.

Synechocystis

Another example is Synechocystis, a genus of cyanobacterium. Cyanobacterium do not have flagella. Nonetheless, Synechocystis species can move in cell suspensions and on moist surfaces and by using retractile type IV pili, displaying an intermittent two phase motion; a high-motility run and a low-motility tumble (see diagram). The two phases can be modified under various external stressors. Increasing the light intensity, uniformly over the space, increases the probability of Synechocystis being in the run state randomly in all directions. This feature, however, vanishes after a typical characteristic time of about one hour, when the initial probability is recovered. These results were well described by a mathematical model based on a linear response theory proposed by Vourc’h et al.

Synechocystis cells can also undergo biased motility under directional illumination. Under directional light flux, Synehcocystis cells perform phototactic motility and head toward the light source (in positive phototaxis). Vourc’h et al. (2020) showed that this biased motility stems from the averaged displacements during run periods, which is no longer random (as it was in the uniform illumination). They showed the bias is the result of the number of runs, which is greater toward the light source, and not of longer runs in this direction. Brought together, these results suggest distinct pathways for the recognition of light intensity and light direction in this prokaryotic microorganism. This effect can be used in the active control of bacterial flows.

It has also been observed that very strong local illumination inactivates the motility apparatus. Increasing the light intensity of more than ~475 μmol m⁻² s⁻¹ reverses the direction of Synechocystis cells to move away from the high levels of radiation source. Moreover, Synechocystis cells show a negative phototaxis behavior under ultraviolet radiation as an effective escape mechanism to avoid damage to DNA and other cellular components of Synechocystis. Contrary to the run phase that can extend from a fraction of a second to several minutes, the tumble lasts only a fraction of a second. The tumbling phase is a clockwise rotation that allows the cell to change the motility direction of the next run.

Chemotaxis is another scheme that allows an organism to move toward or away from gradients of nutrients or other chemical stimuli. Detecting by transmembrane chemoreceptors the microorganism performs a three-dimensional random walk is observed in a homogenous environment, and the direction of each run is identified after a tumble.

Mathematical modeling

Theoretically and computationally, run-and-tumble motion can be modeled as a stochastic process. One of the simplest models is based on the following assumptions:

• Runs are straight and performed at constant velocity v₀ (initial speed-up is instantaneous)
• Tumbling events are uncorrelated and occur at average rate α, i.e., the number of tumbling events in a given time interval has a Poisson distribution. This implies that the run durations are exponentially distributed with mean α⁻¹.
• Tumble duration is negligible
• Interactions with other agents are negligible (dilute limit)

With a few other simplifying assumptions, an integro-differential equation can be derived for the probability density function f (r, ŝ, t), where r is the particle position and ŝ is the unit vector in the direction of its orientation. In d-dimensions, this equation is

${\displaystyle \frac{\mathrm{\partial }f(\mathbf{r},\hat{\mathbf{s}},t)}{\mathrm{\partial }t}+v_0\hat{\mathbf{s}}\cdot \mathrm{\nabla }f(\mathbf{r},\hat{\mathbf{s}},t)={\xi }^{-1}\mathrm{\nabla }\cdot (\mathrm{\nabla }V(\mathbf{r})f(\mathbf{r},\hat{\mathbf{s}},t))-\alpha f(\mathbf{r},\hat{\mathbf{s}},t)+\frac{\alpha }{{\mathrm{\Omega }}_d}\int g(\hat{\mathbf{s}}-{\hat{\mathbf{s}}}^{\prime })f(\mathbf{r},{\hat{\mathbf{s}}}^{\prime },t)d{\hat{\mathbf{s}}}^{\prime }}$

where Ω_d = 2π^d/2/Γ(d/2) is the d-dimensional solid angle, V(r) is an external potential, ξ is the friction, and the function g (ŝ - ŝ') is a scattering cross section describing transitions from orientation ŝ' to ŝ. For complete reorientation, g = 1. The integral is taken over all possible unit vectors, i.e., the d-dimensional unit sphere.

In free space (far from boundaries), the mean squared displacement ⟨r(t)²⟩ generically scales as ⟨r(t)²⟩ ~ t² for small t and ⟨r(t)²⟩ ~ t for large t. In two dimensions, the mean squared displacement corresponding to initial condition f (r, ŝ, 0) = δ(r)/(2π) is

${\displaystyle ⟨{\mathbf{r}}^2⟩=\frac{2v_0^2}{{\alpha }^2(1-{\sigma }_1)}\left[\alpha t+\frac{e^{-\alpha (1-{\sigma }_1)t}-1}{1-{\sigma }_1}\right]}$

where

${\displaystyle {\sigma }_1={\int }_0^{2\pi }g(\hat{\mathbf{s}})(\cos \theta )d\theta }$

with ŝ parametrized as ŝ = (cos θ, sin θ).

In real-world systems, more complex models may be required. In such cases, specialized analysis methods have been developed to infer model parameters from experimental trajectory data.

The mathematical abstraction of run-and-tumble motion also appears outside of biology—for example, in idealized models of radiative transfer and neutron transport.