# Migration and accumulation of bacteria with chemotaxis and chemokinesis

Theresa Jakuszeit, James Lindsey-Jones, François J Peaudecerf, Ottavio A Croze

► To cite this version:

Theresa Jakuszeit, James Lindsey-Jones, François J Peaudecerf, Ottavio A Croze. Migration and accumulation of bacteria with chemotaxis and chemokinesis. European Physical Journal E: Soft matter and biological physics, 2021, 44 (3), pp.32. 10.1140/epje/s10189-021-00009-w. hal-03897600

HAL Id: hal-03897600

https://hal.science/hal-03897600

Submitted on 14 Dec 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Migration and accumulation of bacteria with chemotaxis and chemokinesis

Theresa Jakuszeit<fnref n="1" />,ª,*, James Lindsey-Jones<fnref n="1" />, François J. Peaudecerf<fnref n="2" />, and Ottavio A. Croze<fnref n="1" />,<fnref n="3" />,ᵇ

Received 17 September 2020 / Accepted 5 January 2021

© The Author(s) 2021

**Abstract** Bacteria can chemotactically migrate up attractant gradients by controlling run-and-tumble motility patterns. In addition to this well-known chemotactic behaviour, several soil and marine bacterial species perform chemokinesis; they adjust their swimming speed according to the local concentration of chemoeffector, with higher speed at higher concentration. A field of attractant then induces a spatially varying swimming speed, which results in a drift towards lower attractant concentrations—contrary to the drift created by chemotaxis. Here, to explore the biological benefits of chemokinesis and investigate its impact on the chemotactic response, we extend a Keller-Segel-type model to include chemokinesis. We apply the model to predict the dynamics of bacterial populations capable of chemokinesis and chemotaxis in chemoeffector fields inspired by microfluidic and agar plate migration assays. We find that chemokinesis combined with chemotaxis not only may enhance the population response with respect to pure chemotaxis, but also modifies it qualitatively. We conclude presenting predictions for bacteria around dynamic finite-size nutrient sources, simulating, e.g. a marine particle or a root. We show that chemokinesis can reduce the measuring bias that is created by a decaying attractant gradient.

## 1 Introduction

Many bacteria are able to swim by rotating helical flagella distributed on their cell body, and control their swimming pattern by modulating the speed or direction of rotation of their flagellar motors. For example, in the model run-and-tumble motion of *Escherichia coli* [1], a bacterium swims approximately straight in a ‘run’ by rotating its flagella in a bundle. When some flagella change their rotation direction, the bundle comes apart and the cell randomly changes direction in a ‘tumble’. In the absence of external bias, this microscopic pattern resembles a random walk, and leads to macroscopic diffusion of a bacterial population. In the presence of a chemical gradient, the random walk is biased, a response known as chemotaxis. As shown in Fig. 1a, a bacterium achieves the biased motion up a gradient of attractant by varying the frequency of tumbles in its random walk; if the bacterium moves up the gradient, the tumbling rate $\alpha$ decreases and, thus the run length increases, while the swimming speed remains constant.

The run-and-tumble model was originally introduced for enteric bacteria such as *E. coli* and *Salmonella*

*typhimurium* [2,3], which commonly live in nutrient-rich environments, such as the gut. Marine and soil bacteria, however, often experience heterogeneous and nutrient-scarce environments, and have been found to display different motility patterns. For example, several species living in these harsher environments respond to higher concentrations of attractant by increasing their speed [4–6]. This response, known as ‘chemokinesis’, modifies the swimming speed in response to the local chemical concentration without affecting the tumbling rate, as shown in Fig. 1b. A positive chemokinetic response leads to a higher swimming speed at higher attractant concentrations, whereas a negative response lowers the speed at those concentrations. The strength of positive chemokinetic response can be defined as relative increase in swimming speed over the speed in the absence of chemokinetic effector. A wide range of the chemokinetic response strength has been reported, even for a single species. The responses have been found to vary for the symbiotic soil bacteria *Sinorhizobium meliloti* and *Azospirillum basilense* from 7.5 to 35% [7–9] and 40 to 77% [10], respectively; 7.5 to 73% for the soil and freshwater purple bacterium *Rhodobacter sphaeroides* [11,12]; 26 to 53% for the enterobacterium *E. coli* [13]; 48% [5] or 6 to 64% [14] for the marine pathogen *Vibrio coralliilyticus*. The marine bacterium *V.alginolyticus* showed an increase of up to 80% upon

<footnote n="1">¹ Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, UK</footnote>

<footnote n="2">² Department of Civil, Environmental and Geomatic Engineering, Institute of Environmental Engineering, ETH Zürich, 8093 Zürich, Switzerland</footnote>

<footnote n="3">³ Present address: School of Mathematics, Statistics and Physics, Newcastle University, Newcastle upon Tyne NE1 7RU, UK</footnote>

<footnote>ªe-mail: theresa.jakuszeit@curie.fr (corresponding author)</footnote>

<footnote>ᵇe-mail: o.croze@physics.org (corresponding author)</footnote>
stimulation with glucose in [15]. However, to the best of our knowledge, the pure chemokinetic speed increase as a function of attractant concentration has not been systematically measured for any of these species.

While the role of chemokinesis has been studied extensively in *Paramecium* spp. and other protozoa [16–19], the biological significance of the chemokinetic response of marine and soil bacteria has yet to be fully elucidated. Based on the environment that chemokinetic bacteria have been found in, (positive) chemokinesis might be beneficial in heterogeneous environments with scarce sources of nutrients (attractants). For example, alga-sized microbeads coated with various amino acids were used to study the response of marine bacteria to point-like sources of attractants [4]. All marine bacteria studied were observed to accumulate in bands around the point-like sources while displaying a chemokinetic response. Furthermore, chemokinesis could allow marine bacteria to track algae, helping to foster symbioses with these microorganisms, as well as permitting to respond quickly to short bursts of nutrients, such as those generated from lysing algae [20]. Another example of a chemokinetic marine bacterium is the coral pathogen *V. coralliolyticus*. Microfluidic experiments on this pathogenic bacterium in combination with mathematical modelling have suggested that the maximum accumulation in response to chemical cues produced by heat-stressed coral hosts is larger and is reached faster than in the absence of chemokinesis [5, 14]. As heat-stressed corals are more susceptible to pathogens, chemokinesis could be a crucial evolutionary advantage in oceans heating up due to climate change. In fact, the chemokinetic response was shown to be even stronger at elevated temperatures increasing from 6% at 20°C to 64% at 30°C [14].

Recent interest in chemokinesis has also been sparked by synthetic microswimmers, such as Janus particles. Janus particles are synthetic colloids in a bath of fuel (e.g. $H_2O_2$) that propel due to an asymmetric chemical reaction on their surface [21]. These particles show a positive chemokinetic response since their swimming speed increases with increasing fuel concentration [22], and therefore accumulate in areas of lower fuel concentration.

To date, theoretical work on the combination of positive chemokinesis and chemotaxis has focussed on single-cell level using agent-based models [5, 15, 23]. The chemokinetic response of the marine pathogen *V. coralliolyticus* has been modelled as a step increase in swimming speed beyond a threshold attractant concentration [5]. This model was used to analyse the chemotactic response to a transient attractant gradient in a microfluidic device after an initial release of attractant. The chemotactic index (i.e. the enhancement in cell concentration over a control region) suggests that chemokinesis enables a stronger and faster response. This model has been further adapted to include speed dependent changes in the probability of flicking and reorientation frequency of *V. alginolyticus* [15, 24]. In this particular case, the speed induced changes in motility pattern are responsible for a significant part of the

<figure><img src="image_2.png" /><figcaption>**Fig. 1** Chemotaxis versus chemokinesis: a Chemotaxis is the biased movement of bacteria up a chemical gradient of attractant by reducing the tumbling rate $\alpha$ and, thereby increasing the length of runs in a favourable direction. b Positive chemokinesis leads to an increase in swimming speed $v$ in response to an increase in the local attractant concentration</figcaption></figure>

chemotaxis improvement as shown by the agent-based model in [15].

In this work, we use a continuum model to study the spatio-temporal dynamics of bacterial populations with chemokinetic and chemotactic responses. We incorporate chemokinesis into the standard Keller–Segel model for chemotaxis by deriving the model from microscopic run and tumble dynamics. The model is then used to obtain analytical conditions for chemokinetic drift, and solved numerically for three different example attractant distributions that are inspired by existing experimental systems.

## 2 Model

We derive a model for chemotaxis in combination with positive chemokinesis by considering a one-dimensional system in which a cell can move either to the right or left with speed $v$. In this system, a right (left) moving particle changes direction with rate $\alpha_R$ ($\alpha_L$). Following previous approaches [25, 26], the one-particle probability density for bacteria evolves according to $\partial b/\partial t = -\partial J/\partial x$, and the bacterial flux $J$ can be derived as

$$J(x) = -D_b(x) \frac{\partial b}{\partial x} + V(x)b \quad (1)$$

$$\text{with} \quad D_b(x) = \frac{v^2}{2\alpha}, \quad V(x) = V_k + V_\chi, \quad (2)$$

where $\alpha_R + \alpha_L = 2\alpha$, and where we have defined the chemokinetic

$$V_k = -\frac{v}{\alpha} \frac{\partial v}{\partial x} \quad (3)$$

and chemotactic
$$V_{\chi} = \frac{v}{2\alpha}(\alpha_{L} - \alpha_{R}), \qquad (4)$$

drift speeds, respectively.

We will now connect the bacterial flux (1) to the commonly used Keller–Segel model of chemotaxis, adapted to the description of chemokinetic populations in dynamic environments.

## 2.1 Chemotaxis

For chemotactic populations, the drift-diffusion flux (1) is coupled to a chemoattractant density field c via the chemotactic drift speed, given by Eq. (4). In the standard Keller–Segel model, it is phenomenologically asserted that this chemotactic drift speed is proportional to the change in the chemical attractant in space, $V_{\chi} = \chi \nabla f_{\chi}$, where $\chi$ is the chemotactic sensitivity parameter and $f_{\chi}$ is a function of c that ensures that the chemotactic drift is biased towards higher attractant concentrations [27]. This definition of the chemotactic drift speed assumes that the chemical attractant profile changes in space but not in time. However, as mentioned, known chemokinetic bacteria inhabit dynamic environments, such as in the ocean or soil. As pointed out by Hein et al [28], the effective gradient perceived by a bacterium changes in a temporally varying attractant profile depending on the direction of its run. Consider a source of attractant at one point of space, and the associated negative gradients of concentration as one gets away from it. If these gradients are steady in time, for example as in a microfluidic setting, a bacterium exploiting the attractant landscape will correctly detect the gradient and swim towards the source, as in standard chemotaxis. However, if the source corresponds to a single point-like release of attractant, and the concentration at the source position also decays due to diffusion, a bacterium travelling towards the source perceives a smaller increase (or even a decrease) in concentration compared to the steady case. On the other hand, a bacterium moving away from the source perceives a decrease, which is reduced in magnitude compared to the steady case. We show in “Appendix A” how this influences the mean run length of a bacterium, and, in light of this, modify the chemotactic drift speed to

$$V_{\chi} = \chi \left( \nabla f_{\chi} + \frac{1}{v} \partial_{t} f_{\chi} \right), \qquad (5)$$

where we see that the presence of a temporally increasing (decreasing) gradient increases (decreases) the chemotactic drift speed. Further, we note that the perturbation to the chemotactic drift speed is smaller the larger the value of the swimming speed. For the chemotactic function $f_{\chi}$, we may choose, for example [29],

$$f_{\chi} = \frac{c(x, t)}{c(x, t) + k_{\chi}}. \qquad (6)$$

Following [30], using relationship (4) for the chemotactic speed, the microscopic swimming parameters v

and $\alpha$ can be related to the macroscopic parameter $\chi$ (see “Appendix A” for details) as

$$\chi = \frac{v^2}{\alpha} \beta, \qquad (7)$$

where $\beta$ is constant dependent on the chemotactic response (again, see “Appendix A” for details).

## 2.2 Chemokinesis

In this section, we consider how chemokinesis can be modelled. We can modify Eq. (7) to include chemokinesis (a spatially varying swimming speed) to give

$$\chi(x) = \chi_0 \frac{v(x)^2}{v_0^2}, \qquad (8)$$

in which the subscript refers to the parameters in the absence of chemokinesis. (The interested reader can find a detailed derivation in “Appendix A”.) Here we assumed that the swimming speed is constant during a run.

In addition to modifying the chemotactic drift $\chi$, the spatially varying swimming speed of chemokinetic bacteria causes an additional chemokinetic drift to arise with a speed given by Eq. (3). When acting alone, this drift drives the cells towards regions of lower speed [25].

Next we consider how to quantify the chemokinetic coupling between speed and local concentration of attractant. We need to make assumptions about the relationship between speed and attractant concentration, as experimental studies have not been carried out to provide this. Firstly, we assume that cells swim at a base-level speed, $v_0$. Secondly, we reasonably posit that chemokinesis monotonically increases the swimming speed up to a maximum speed denoted as $v_0 + v_c$. The dynamics between these limits are given by an unknown function characterizing the chemokinetic response. In the following, we choose a Hill-type equation

$$v(x) = v_0 + v_c \frac{c(x)^n}{c(x)^n + k_c^n} \qquad (9)$$

to approximate the chemokinetic response. The Hill parameter n allows us to introduce an inflection point and change the gradient $\partial v / \partial c$. As shown in Fig. 2, the speed increases monotonically with increasing attractant concentrations c for all n and $v_c > 0$, where the half-maximum speed is reached at the attractant concentration $k_c$. For any n, $v = v_0 + k_c/2$ at c = $k_c$ as all functions covered by Eq. (9) have the same half-saturation constant. Note that upon setting n = 1, one recovers a Michaelis-Menten-type response; for n → ∞, Eq. (9) approaches a step function, which has been used previously to approximate the chemokinetic response [5,15].

Finally, in situations where bacterial growth cannot be ignored (e.g. migration across agar plates considered
<figure><img src="image_3.png" /><figcaption>Fig. 2 Chemokinetic response function. The swimming speed increases from a reference speed, $v_0$, by at most $v_c$, depending on the local concentration of the attractant according to a Hill-type function. The parameter $n$ determines the strength of the gradient at the half-saturation concentration, $k_c$</figcaption></figure>

in Sect. 3.2), bacteria are also assumed to undergo logistic growth, which comprises growth and death terms [29]. We note, however, that the chemoeffectors eliciting chemotaxis and chemokinesis need not in general be metabolizable nutrients that induce growth. The full model equations for the chemical attractant field and the chemotactic and chemokinetic bacterial population field, therefore, are

$$ \frac{\partial c}{\partial t} = D_c \frac{\partial^2 c}{\partial x^2} - k_g g(c) \frac{b}{Y} \quad (10a) $$

$$ \frac{\partial b}{\partial t} = -\frac{\partial J}{\partial x} + k_g g(c)b \left(1 - \frac{b}{k_b}\right) \quad (10b) $$

$$ J(x) = -\frac{v^2}{\alpha} \frac{\partial b}{\partial x} - \frac{v}{\alpha} \frac{\partial v}{\partial x} b + \chi(x) \left( \nabla f_{\chi} + \frac{1}{v} \partial_t f_{\chi} \right) b, \quad (10c) $$

where $k_g$ is the maximum growth rate, $Y$ is the bacterial yield, and $k_b$ is the carrying capacity. Furthermore, $g(c)$ is chosen as a Monod-type growth function $g(c) = c/(c+k_s)$ with the half-saturation constant $k_s$. Note that $\chi(x)$ is given by Eq. (8), $f_\chi$ by Eq. (6) and $v(x)$ by Eq. (9).

Let us summarize the effects of chemokinesis in our model. A spatially varying speed affects all three terms of the bacterial flux $J$ in Eq. (10c) as: (i) there are regions with a higher diffusivity since $D_b \sim v^2$ (first term); (ii) it introduces a drift where $\partial v/\partial x \neq 0$ (second term); and (iii) there are regions with a larger chemotactic drift as $\chi \sim v^2$ (third term).

## 2.3 Non-dimensionalization

The system of partial differential equations (10) is non-dimensionalized using the characteristic time and

length scales $t_0 = k_g^{-1}$ and $x_0 = \sqrt{t_0 D_b^0}$, where the bacterial diffusivity is $D_b^0 = v_0^2 \alpha^{-1}$. We rescale the attractive and bacterial densities by their respective initial densities, $c_0$ and $b_0$. The system of PDEs in dimensionless form thus reads

$$ \frac{\partial C}{\partial T} = N \nabla^2 C - HB \frac{C}{C + K_S} \quad (11a) $$

$$ \frac{\partial B}{\partial T} = -\nabla J + B \frac{C}{C + K_S} (1 - B) \quad (11b) $$

$$ J = -\mathcal{V}(X)^2 \nabla B + V_k B + V_\chi B \quad (11c) $$

$$ V_k = -\mathcal{V}(X) \frac{n\eta\omega^n C^{n-1}}{(C^n + \omega^n)^2} \nabla C \quad (11d) $$

$$ V_\chi = \mathcal{V}(X)^2 \frac{\delta_0 K_\chi}{(C + K_\chi)^2} \left( \nabla C + \frac{\zeta}{\mathcal{V}(X)} \frac{\partial C}{\partial T} \right) \quad (11e) $$

with non-dimensional parameters $B = b/b_0$, $C = c/c_0$, $N = D_c/D_b^0$, $H = b_0/(Yc_0)$, $K_S = k_s/c_0$, $\eta = v_c/v_0$, $\omega = k_c/c_0$, $\delta_0 = \chi_0/D_b^0$, $\zeta = v_0/(\alpha x_0)$ and $K_\chi = k_\chi/c_0$, and non-dimensional time $T = t/t_0$ and space $X = x/x_0$. The non-dimensional speed function in Eq. (11) is

$$ \mathcal{V}(X) = 1 + \eta \frac{C^n}{C^n + \omega^n}, \quad (12) $$

where $\eta$ is the maximum increase in swimming speed, $\omega$ is the attractant concentration at which the half maximum speed increase is reached, and $n$ is the Hill parameter. Note that $\omega$ corresponds to the inflection point of Eq. (12). The chemokinetic response is positive for $\eta > 0$.

Finally, the model is extended to the 2D axisymmetric case by introducing polar coordinates as

$$ \nabla B = \frac{\partial B}{\partial R} \quad (13a) $$

$$ \nabla^2 B = \frac{\partial^2 B}{\partial R^2} + \frac{1}{R} \frac{\partial B}{\partial R}, \quad (13b) $$

and the equivalent equations for the chemical field, $C$. The details of the numerical solution and simulations of the model are described in “Appendix B”, including a summary of parameters used. The parameters for the results presented in the main text were chosen to best illustrate the chemokinetic effect (see “Appendix B” for a discussion of parameter values and the Supplementary Information for further simulations with different parameter sets, including $n$). In Sect. 3, we will present the simulation results for three different types of attractive gradient.

## 2.4 Condition for dominant chemokinetic drift

Before solving the extended Keller–Segel model numerically, we can use Eq. (11) to analytically derive a con-
dition on the relative importance of chemokinetic and chemotactic contributions to the drift of the bacterial population. The drift due to a spatially varying swimming speed causes cells to accumulate in regions where they have low speeds. For $\eta > 0$, by construction of the velocity function (12), the speed is low at low attractant concentrations. The chemotactic drift, on the other hand, is directed towards higher attractant concentrations by virtue of Eq. (6). Hence, the bacterial density is governed by two competing drifts, as can be seen from the opposing signs in Eqs. (11d) and (11e). If the chemokinetic drift is larger than the chemotactic drift for a large part of the spatial domain, this could lead to accumulation at low attractant concentrations, instead of the biologically desirable accumulation at high concentrations. Assuming a stationary and linear attractant profile (i.e. $\partial C/\partial T = 0$ and $\partial C/\partial X = \text{const}$), we have from Eqs. (11d) and (11e) that the chemokinetic drift is larger than the chemotactic drift if

$$ \frac{n\eta\omega^n}{(C^n + \omega^n)^2} \nu^{-1} C^{n-1} > \frac{\delta_0 K_\chi}{(C + K_\chi)^2} \quad \forall X \in \Omega, \quad (14) $$

where $\Omega$ is the spatial domain. In the case of a linear attractant profile, we know from Eq. (12) (and Fig. 2) that the gradient $\partial V/\partial C$ is maximum close to the half-saturation constant $\omega$. Thus, we evaluate condition (14) with $C = \omega$, which yields the condition

$$ n > \frac{4\delta_0\omega K_\chi}{(\omega + K_\chi)^2} \left( \frac{1}{\eta} + \frac{1}{2} \right) \quad \forall X \in \Omega : C(X) = \omega. \quad (15) $$

If the Hill parameter $n$ exceeds this threshold, the chemokinetic drift is predicted to be larger than the chemotactic drift at the attractant concentration $C^* = \omega$ for $\eta > 0$. Note that, for a step function, condition (15) is always met at the threshold concentration $C^* = \omega$ since $n \to \infty$. Conversely, if $\eta < 0$ (i.e. modelling a negative chemokinetic response), condition (15) will never be met as the chemotactic and chemokinetic drift have the same direction (chemokinesis in this case is stabilizing).

# 3 Results

## 3.1 Steady linear attractant profile

The full model Eq. (11) includes the effect of growth and consumption as well as chemotaxis and chemokinesis. Thus, any change in the spatial distribution of the bacterial population due to chemotaxis and chemokinesis will feed back onto the attractant distribution due to consumption by the bacteria. In order to identify the influence of chemokinesis without such additional complications, we first solve the model for a steady attractant gradient, i.e. $\partial C/\partial T = 0$. We furthermore assume a linear profile such that $\partial C/\partial X = \text{const}$ and ignore consumption, and thus population growth. This situ-

<figure><img src="image_4.png" /><figcaption>Fig. 3 Chemokinesis in steady linear attractant profile. Bacterial response (bottom row) to a fixed linear attractant profile (top row) using pure chemotaxis (blue curve) versus chemotaxis with chemokinesis (orange curve) for three time points. Note the changing range of the y-axis in the bottom row. The initial bacterial profile is indicated by the grey line. The position at which $C = \omega$ is highlighted by a dotted red line. Parameters $H = 0$, $N = 0$, $K_\chi = 0.53$, $\delta_0 = 50$, $\eta = 2$, $\omega = 0.2$, $n = 5$, $T = 1, 9, 17$</figcaption></figure>

ation might be achieved experimentally in microfluidic devices [31], where the gradient may be fixed and bacterial growth can be neglected on experimental timescales short compared to growth timescales $\sim k_g^{-1}$.

Figure 3 compares the response of a purely chemotactic population to the response of a chemotactic-chemokinetic population. Chemokinesis leads to a stronger and faster accumulation than in the purely chemotactic case. At the critical concentration, $C = \omega$, however, the chemokinetic drift holds back a subset of the population because it is directed towards lower attractant concentrations as described in Sect. 2.4. This can be seen in the form of an accumulation of cells at low attractant concentrations. As the chemotactic sensitivity parameter, $\delta_0$, is large in this simulation, the population subset overcomes the drift and accumulates at high attractant concentrations at long times. However, if condition (15) is fulfilled, the chemokinetic drift is larger than the chemotactic drift. Thus, there is a sub-population driven to small attractant concentrations by chemokinesis. The effect of varying the Hill parameter $n$ in Eq. (12) is shown in Fig. 4. For the parameters chosen in Fig. 3, condition (15) is met for $n > 39.78$. In experiments, this would require observing the transient bacterial concentration profiles in addition to the commonly reported steady-state profiles.

As the attractant concentration is fixed, we can determine a steady-state for the bacterial population, i.e. by setting $\partial B/\partial T = \partial J/\partial X = 0$. Due to the homogeneous Neumann boundary conditions of the problem, we have $J = 0$ in Eq. (11), which yields the differential equation
<figure><img src="image_5.png" /><figcaption>Fig. 4 Effect of Hill parameter on chemokinesis in steady linear attractant profile. Chemotactic–chemokinetic bacterial response to a fixed linear attractant profile at $T = 0.5$ (same constant attractant profile as in Fig. 3). The Hill parameter $n$ in Eq. (12) and, thus, the speed gradient is varied. The position at which $C = \omega$ is highlighted by a dotted red line. Parameters: $H = 0$, $N = 0$, $K_{\chi} = 0.53$, $\delta_0 = 50$, $\eta = 2$, $\omega = 0.2$, $n = 1, 5, 10, 40$. For the chosen parameters, $n = 40$ is just above the threshold given by (15)</figcaption></figure>

$$ \frac{\partial B}{\partial X} + B \left[ \frac{1}{2D} \frac{\partial D}{\partial X} - \delta_0 \frac{d}{dX} \left( \frac{C}{C+K_{\chi}} \right) \right] = 0, \quad (16) $$

where we used $D(X) = V(X)^2$. This equation can be integrated to give the steady state

$$ \frac{B}{B^*} = \frac{V^*}{V} \exp \left\{ \delta_0 \left( \frac{C}{C+K_{\chi}} - \frac{C^*}{C^*+K_{\chi}} \right) \right\}, \quad (17) $$

where $B^*, C^*, V^*$ are reference values at a chosen reference point $X^*$. It is thus clear that in addition to the influence on the dynamics, chemokinesis affects the steady-state solution via the term $V^*/V$, where $V$ varies in space due to chemokinesis. In the case of chemokinesis but no chemotaxis (i.e. $\delta_0 = 0$), the steady-state is determined by the inverse of the speed distribution, i.e. the bacteria accumulate at low speed, as expected and shown previously [25,32,33]. For nonzero $\delta_0$, if the speed is uniform in space, $V^*/V = 1$ and Eq. (17) reduces to the chemotactic steady-state solution. The exponential term in Eq. (17) represents the chemotactic contribution to the steady state, which does not depend on the swimming speed. Thus, the increase in chemotactic sensitivity $\chi$ [see Eq. (8)] must be balanced by the increase in diffusivity at steady state in a fixed chemical gradient. However, chemokinesis still affects the steady state via the term $V^*/V$. This chemokinetic effect in the steady-state may only be detectable in experiments with small $\delta_0 = \chi_0/D_b^0$, since the chemotactic exponential term will dominate $V^*/V$ for large $\delta_0$.

## 3.2 Self-generated gradient: agar plate

We now consider an evolving attractant field with consumption and growth of bacteria. The attractant is ini-

<figure><img src="image_6.png" /><figcaption>Fig. 5 Self-generated gradient. Bacterial populations (bottom row) create an attractant gradient (top row) via consumption, which they respond to with chemotaxis at base speed $V = 1$ (blue curve) or chemotaxis–chemokinesis (orange curve). Chemokinesis leads to a faster but also broader, less pronounced bacterial wave. The chemotactic population travelling at constant speed $V = 1 + \eta$ (purple curve) has the fastest travelling pulse. Parameters $H = 3.5$, $K_S = 1$, $N = 0.5$, $K_{\chi} = 0.53$, $\delta_0 = 105$, $\eta = 0.5$, $\omega = 0.5$, $n = 5$, $T = 1, 8.8, 16.4$</figcaption></figure>

tially uniformly distributed in a 2D axisymmetric setting. This set-up is reminiscent of the classical agar plate experiments, in which bacteria are inoculated in the centre of a nutrient agar plate, see e.g. [29,34]. While growing and consuming nutrient, the population creates a gradient of attractant, which it then follows outwards in a chemotactic wave. The attractant profile is a travelling wave itself, and we assume here that the profile relative to the bacterial travelling wave is stationary, i.e. $\zeta = 0$ in the chemotactic drift (11e).

In Fig. 5, we compare two chemotactic populations to a chemotactic–chemokinetic population. The chemotactic populations travel at a constant speed, either $V = 1$ (blue curve) or $V = 1 + \eta$ (purple curve). Both populations develop a sharp travelling wave, with a larger wave speed for the population at speed $V = 1 + \eta$. The chemotactic–chemokinetic bacterial population, on the other hand, develops a broader wave profile. The peak of the wave front is smaller and is followed by a plateau. This effect is more pronounced at late times, as can be seen in the third lower panel of Fig. 5. The reduced pulse also travels slower than the pulse of the chemotactic population at elevated speed $V = 1 + \eta$ because the front speed scales with the number of bacteria in the pulse [35]. This observation might explain why in agar plate experiments testing for chemotaxis, chemokinetic species such as *Sinorhizobium meliloti* lack the sharp bands [8,36], which are known to be a hallmark of chemotaxis for other species, e.g. *E. coli* [29,34].

Chemokinesis confers an additional biological advantage in the form of increased population growth as shown in Fig. 6, which shows the integrated number of cells over time. At any point in time, the chemotactic population is smaller than the chemotactic–
<figure><img src="image_7.png" /><figcaption>**Fig. 6** Population growth. The size of the populations in Fig. 5 is the bacterial density integrated over the simulation domain. Faster travelling waves of the chemokinetic-chemotactic population (orange curve) and the chemotactic population at constant speed $V = 1 + \eta$ (purple curve) also induce a faster population growth due to consumption of nutrients, compared to the chemotactic population traveling at $V = 1$ (blue curve). Parameters as in Fig. 5</figcaption></figure>

chemokinetic population. However, the chemotactic population at elevated speed, $V = 1 + \eta$, (purple curve) shows a stronger population growth than the chemotactic-chemokinetic population due to the faster travelling pulse, which is caused by the increased swimming speed. This increase in swimming speed is associated with metabolic cost [37]. Thus, permanently swimming faster independent of the attractant concentration could be a beneficial strategy, if the metabolically available energy is not constrained by nutrient supply. When nutrient concentrations are low, on the other hand, increasing swimming speed provides no benefit to bacteria and metabolism is a limiting factor. Chemokinesis could provide an advantageous speed enhancement when it is both metabolically affordable and beneficial [4,20]. While the situation considered in this section assumed an abundant supply of chemoeffectors (and in this case nutrients), the next section will consider the response to a transient burst of chemoeffectors.

### 3.3 Transient source

A localized burst of chemoeffector may, e.g., occur in the sea if algae/phytoplankton lyse and release their content, as has been recently studied in the laboratory [38], or when marine particles exude plumes of chemoeffector [39]. In soil, plant roots exude sugars and other potential nutrients, which locally create a high concentration of chemical attractants. In the following, we consider a single strong pulse of chemoeffector originating from a finite-size axisymmetric source that dissipates via diffusion, modelling a potential dynamic environment around roots or marine particles. The attractant profile that develops is $C(R, T) = S(4\pi NT)^{-1} \exp(-R^2/4NT)$, with S representing the amount of chemoeffector contained in the pulse in non-

<figure><img src="image_8.png" /><figcaption>**Fig. 7** Time evolution of the bacterial response to a diffusing attractant from a transient source. Bacterial populations (bottom row) are attracted to source of diffusing attractant (top row). The chemokinetic-chemotactic population (orange curve) shows a faster and stronger accumulation than the purely chemotactic population (blue curve). Parameters $H = 3.5$, $K_S = 1$, $N = 0.5$, $K_\chi = 0.53$, $\delta_0 = 50$, $\eta = 2$, $\omega = 0.2$, $n = 1$, $S = 0.5$, $T = 0.01, 0.05, 0.64$; no bacterial growth</figcaption></figure>

dimensional units. The attractant profile is shown in Fig. 7 (top row), while the bacteria are initially uniformly distributed in the domain at concentration 0.2. To model the response to such a transient attractant profile, we need to include the chemotactic drift velocity, Eq. (11e) with $\zeta > 0$, modified to account for the effective gradient perceived by bacteria as they traverse the temporally varying pulse (see “Appendix A”).

As shown in Fig. 7, bacteria with chemotaxis and chemokinesis display a faster and stronger response to a chemoeffector pulse than those with chemotaxis alone. This strong accumulation occurs in spite of the fact that, because of chemokinesis, diffusivity close to the source is higher for these bacteria. For the transient pulse under consideration, temporal variations in the chemoeffector concentration ($\sim \partial_t f_\chi$) need to be considered when modelling the chemotactic response. Indeed, the run duration for a bacterium travelling up/down the gradient is $\tau^{L,R} \propto \pm(v\nabla f_\chi + \partial_t f_\chi)$ (see “Appendix A”). Thus, when v is higher, as for chemokinetic bacteria, the effect of temporal variation on the chemotactic bias of tumbles is reduced, increasing the accuracy of the chemotactic response. This effect contributes to accounting for the stronger accumulation of chemotactic-chemokinetic bacteria. To illustrate this further, in Fig. 8, we plot the number of bacteria accumulated at the attractant source, $B_S$ (the maximum of the bacterial profiles shown in Fig. 7), as a function of time, for chemotactic and chemotactic-chemokinetic populations. The plot displays bacterial accumulation when the temporal perturbation to the chemotactic response is included in the model [$\zeta > 0$ in (11e)] and when it is not ($\zeta = 0$). In the case of a purely chemotactic population, it can be seen that the predicted amount of bacteria accumulated at the source is lower for a
<figure><img src="image_9.png" /><figcaption>**Fig. 8** Accumulation at a transient source. The bacterial accumulation at the source is reduced due to the reduced drift $V_\chi$ if the transient nature of the attractant profile is taken into account (i.e. $\zeta \neq 0$). Chemokinesis (CK) mitigates for this effect. Parameters $H = 3.5$, $K_S = 1$, $N = 0.5$, $K_\chi = 0.53$, $\delta_0 = 50$, $\zeta = 8.164 \cdot 10^{-3}$, $\eta = 2$ ($\eta = 0$ if no CK), $\omega = 0.2$, $n = 1$, $S = 0.5$; no bacterial growth. If $V_\chi$ includes the effect of transient source, $\zeta = 8.164 \cdot 10^{-3}$, otherwise $\zeta = 0$</figcaption></figure>

model that ignores the temporal perturbation than for one that includes it. Chemokinesis, on the other hand, reduces the relative effect of temporal perturbation so much that there is very little difference between model predictions with $\zeta = 0$ and $\zeta > 0$.

These results suggest that a chemokinetic population might be able to overtake purely chemotactic competitors in response to a sudden nutrient release. While the difference observed in Fig. 8 of at most ~12% may seem small, a corresponding boost to the growth rate can be sufficient to outcompete a purely chemotactic strain within a few generations. In such transient nutrient landscapes, we have further shown that chemokinesis can reduce the adverse effect that a temporal change in attractant profile can have on the chemotactic response.

## 4 Discussion

Chemokinesis is a known response for many environmentally relevant bacteria, yet its consequences for bacterial population dynamics have been little explored. In this work, using a modified Keller–Segel continuum model, we have shown how chemokinesis significantly affects both the dynamics and steady-states of bacterial populations capable of chemotactic and chemokinetic behaviour. Our model incorporates the effects of a concentration dependent speed, including an increase in the chemotactic sensitivity, and a recently suggested modification to the chemotactic response in dynamic environments [28], which was derived adapting the microscopic model first suggested by de Gennes [30].

We have solved our model numerically to explore the effect of chemokinesis on migration and accumulation in experimentally realistic gradients. In a fixed attractant

gradient, our results show that chemokinesis can lead to two subpopulations travelling at different speeds, with the slower one being held back by the chemokinetic drift. In the case of agar plate migration, where bacteria inoculated onto the plate generate their own gradient by consuming nutrients, we find that populations with chemokinesis migrate out from the inoculation point in waves that are faster, but broader than purely chemotactic migrating populations. While the increase in front speed could be explained by a population at a uniformly increased swimming speed, the broadening is only observed in the chemokinetic population. It is a new feature not predicted by previous studies using agent based simulations [15]. It is possible that this broadening might explain why the classic chemotactic Adler bands observed for E. coli [34] are not observed for chemokinetic soil bacteria such as Sinorhizobium meliloti [8,36]. Furthermore, chemokinesis increases the population growth significantly in comparison with purely chemotactic migration.

Our simulations also considered the case of a transient source of nutrients, e.g., a lysed algal cell. In this case, our results show how bacteria with chemokinesis and chemotaxis accumulate faster and more strongly around the source, while concentrations of nutrients are high, with respect to purely chemotactic bacteria. This chemokinetic advantage is both due to the enhanced migration discussed in the previous examples, but also to the fact that chemokinesis mitigates the perturbation to the chemotactic response due to the transient nature of the attractant profile. We note that, while our model includes the effect of transient chemical fields on tumbles, it does not include a recently modelled effect of chemokinesis on the precision of chemosensing [40]. It will be interesting to include this additional effect, which could lead to further enhancements in chemokinetic accumulations, in future elaborations of our model.

The role of metabolism is an important consideration for chemokinesis. For example, positive chemokinesis might be caused purely by an increase in nutrient availability. A resulting increase in the energy level of the cell may lead to more energy being available for flagella rotation, which allows the cell to swim faster. As the chemotactic sensitivity scales with the swimming speed as $\chi \propto v^2$, a faster swimming population will always show a stronger chemotactic drift. However, swimming is associated with a considerable metabolic cost [37, 41, 42]. Indeed, the energetic cost of swimming increases quadratically with the swimming speed [43]. Therefore, it might be beneficial to swim faster (and, thereby, improve chemotaxis) only if favourable nutrient conditions are available. While the energy level of a cell may influence the extent of the chemokinetic response, there are several indications of metabolism-independent chemokinesis for different bacterial species [5, 12, 13]. In a dynamic environment such as considered in Sect. 3.3, chemokinesis can then help to reduce temporal bias and improve the chemotactic response. As pointed out by Hein et al, the threshold for detecting absolute concentration is smaller than for gradi-
ent detection [28]. Thus, chemokinesis can take place
at lower background concentrations to improve chemo-
taxis only when needed.

The predictions of our model include interesting qual-
itative effects, which have not previously been observed
in agent-based models: the slower subpopulation in the
fixed attractant profile, and a broadening of the trav-
elling wave in a self-generated gradient. To test these
predictions experimentally, chemotaxis and chemokine-
sis should be addressed independently. For example,
a recently developed *E.coli* system with a swimming
speed that is controlled via light [33] could be used to
engineer populations with a swimming speed that can
be controlled independently of chemotaxis.

Future theoretical investigations would benefit sig-
nificantly from the experimental measurement of the
chemokinetic response function relating swimming speed
and local attractant concentration. In this work, we
assumed a smooth change from a reference speed to an
increased speed, where the degree of change in swim-
ming speed changed with a single parameter. As we
have shown, a very steep change in swimming speed
(e.g. in form of a step change as assumed previously
in agent-based models [14,15]) could actually inhibit
chemotaxis rather than promote it. Experimental work
so far has been restricted to measurements at very few
different attractant concentrations, which makes it dif-
ficult to deduce a functional relationship between speed
and attractant concentration. Thus, further work is
required to determine the function v(c) for chemoki-
netic bacterial species. Such measurements would also
allow elucidation of the rate of adaptation, i.e. how
quickly the swimming speed adapts to its new value
both for an increase and decrease in attractant con-
centration, and which we have here assumed instanta-
neous. The increase in the population-averaged swim-
ming speed in response to a uniform addition of effec-
tor occurred on the time-scale of 100-200 s for V. algi-
nolyticus [15]. In [12], on the other hand, the chemoki-
netic response of R. sphaeroides was measured within
10s upon uniform addition of chemoeffector, at which
point the swimming speed had already adapted to its
increased level. Furthermore, a desensitization to a sus-
tained higher level of attractant, as observed for chemo-
taxis, does not seem to occur as the swimming speed
remained at elevated levels for hours in R. sphaeroides
and A. brasilense [10,12]. The discrepancy in the order
of magnitude in the response times might be caused by
the experimental set-ups, e.g. in [15], the effector first
needs to diffuse from two sides in a wide microfluidic
channel before a uniform population response can be
measured. To conclude, further experiments could shed
some light on the chemokinetic response function and
adaptation time, which would benefit the further devel-
opment of the model and its assumptions to understand
how bacteria make use of chemokinesis in dynamic envi-
ronments.

**Supplementary information** The online version con-
tains supplementary material available at https://doi.org/
10.1140/epje/s10189-021-00009-w.

**Acknowledgements** We acknowledge financial support from EPSRC EP/L504920/1 and EP/N509620/1 (T.J.) and the Winton Programme for the Physics of Sustainability (T.J., O.A.C.), and funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie Grant Agreement No 798411 (F.J.P.).

Author contribution statement

TJ carried out analytical calculations and simulations.
JL carried out initial calculations and simulations.
FP assisted with initial simulations. TJ, FP and OC
interpreted the simulations. TJ and OC wrote the
manuscript, which was critically revised by FP. All
authors gave final approval for publication and agree
to be held accountable for the work performed therein.

**Data Availability Statement** This manuscript has asso-
ciated data in a data repository. [Authors' comment:  
Code and data associated with this manuscript are avail-  
able from the Zenodo repository: https://doi.org/10.5281/  
zenodo.4323421.]

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Appendix A: Chemotactic drift and sensitiv-
ity

Here, we derive new relations for the chemotactic drift
speed $V_\chi$ and sensitivity parameter $\chi$ for bacteria undergo-
ing chemotaxis and chemokinesis in dynamic environments,
connecting the macroscopic parameter $\chi$ with the micro-
scopic swimming speed $v$ and tumbling rate $\alpha$. We follow
the approach of de Gennes [29,30,44], considering run-and-
tumble bacteria with a simplified *E. coli*-like chemotactic
response. Our contribution for the purpose of this work is
to include: (i) a swimming speed that is a function of posi-
tion; (ii) a temporally varying contribution to the gradient
perceived by the bacteria. As in [29,30], we consider a one-
dimensional model of run-and-tumble bacteria propagating
in gradient of attractant, to which bacteria are exhibiting a
'small response', that is the tumbling response remains close
to the adapted value [45]. Chemotactic memory is mod-
elled through kernel integral, reflecting the fact that bacteria 'remember' their chemical environment over a characteristic delay time of a few seconds. We neglect directional persistence and rotational diffusion, which can be addressed in the same framework [44], but are here ignored for simplicity. The modified de Gennes model we employ could be applicable to chemokinetic species like rhizobia, which display run and tumble dynamics [6]. The model neglects features of the motility pattern of chemokinetic marine bacteria, such as a run-reverse-flick mechanism or the influence of a change in swimming speed on reorientation frequency, as observed for V. alginolyticus [15], but instead focus on the pure effect of speed change.

As the tumbling frequency follows an exponential distribution, tumbling events can be treated as independent events of a Poisson process with rate $\alpha(t) = \alpha_0 e^{-\Delta(t)}$, where $\alpha_0$ is the tumbling rate in the absence of a gradient and where $\Delta(t)$ is the chemotactic bias given by the memory integral

$$ \Delta(t) = \int_{-\infty}^{t} dt' K(t-t') f_{\chi}(x(t')). \quad (A1) $$

The concentration function was chosen as in the main text as $f_{\chi} = c(x, t)/(c(x, t) + k_{\chi})$. The memory kernel $K(t)$ was first measured for E. coli [46] and more recently also for V. alginolyticus [47]. Based on these experiments we assume that the kernel obeys $\int_0^{\infty} K(t)dt = 0$, i.e. the tumbling rate perfectly adapts. As mentioned, we shall assume a small response ($\Delta(t) \ll 1$), so that the bacterial tumble rate can be linearized to

$$ \alpha(t) \approx \alpha_0(1 - \Delta(t)), \quad (A2) $$

Considering a run that starts at $t = 0$, the probability density of a tumble event in the interval $[t, t + dt]$ is $\alpha(t) \exp(-\int_0^t dt' \alpha(t'))$. The mean run duration is given by

$$ \tau = \left\langle \int_{0}^{\infty} \alpha(t) \exp \left( - \int_{0}^{t} dt' \alpha(t') \right) t dt \right\rangle_{\text{paths}}, \quad (A3) $$

where angled brackets denote averaging over all possible bacterial swimming paths, as the tumble rate is path-dependent. (The 'paths' subscript will be omitted henceforth.) Integrating (A3) by parts and recalling (A2), we can write

$$ \tau = \left\langle \int_{0}^{\infty} e^{-\alpha_0 t} \exp \left( \int_{0}^{t} dt' \Delta(t') \right) dt \right\rangle_{\text{paths}}. \quad (A4) $$

Next, since $\Delta(t) \ll 1$, we can linearize the exponential integral to obtain

$$ \tau = \frac{1}{\alpha_0} + \alpha_0 \int_{0}^{\infty} dt e^{-\alpha_0 t} \left\langle \int_{0}^{t} dt' \Delta(t') \right\rangle, \quad (A5) $$

where we have brought the angled brackets inside the time integral to surround only path-dependent quantities. Substituting the expression for the bias (A1) and performing a change of variables in the memory integral (A2) by defining $u = t - t'$, we then obtain

$$ \tau = \frac{1}{\alpha_0} + \alpha_0 \int_0^\infty dt e^{-\alpha_0 t} \times \left\langle \int_0^t dt' \int_0^\infty du K(u) f_\chi(x(t' - u)) \right\rangle. \quad (A6) $$

Next, as the time interval of interest is small compared to the gradient variations, we can Taylor expand the concentration function about a reference position and time:

$$ f_{\chi}(x(t-u)) \approx \text{const} + x(t-u) \nabla f_{\chi} |_{t} + (t-u) \partial_{t} f_{\chi} |_{x(t)}. \quad (A7) $$

We note that the constant term does not influence the integral in (A6), since the response function $K(u)$ integrates to zero. Further analysis is simplified by considering the special response function

$$ K(u) = A\delta(u - \theta), \quad (A8) $$

where $\theta$ is a single delay time [30]. We note that this function does not represent a physical response and is considered solely for the purpose of facilitating the calculation. Then, substituting Eqs. (A8) and (A7) into (A6) and integrating over the delta function, we obtain

$$ \tau \approx \frac{1}{\alpha_0} + A\alpha_0 \int_0^\infty e^{-\alpha_0 t} [I_1(t) + I_2(t)] dt \quad (A9) $$

where we have defined

$$ \begin{aligned} I_1(t) &= \left\langle \int_0^t \nabla f_\chi x(t' - \theta) dt' \right\rangle \\ &= \nabla f_\chi \int_0^t \langle x(t' - \theta) \rangle dt'. \end{aligned} \quad (A10) $$

The average over paths can be taken inside the time integral since tumbles for $t > 0$ are treated as having no effect on cell motion, so that paths are independent of time on the interval considered. Proceeding similarly, we define

$$ I_2(t) = \int_0^t (t' - \theta) \langle \partial_t f_\chi |_{x(t')} \rangle dt'. \quad (A11) $$

To evaluate (A10), we recall that the bacteria we are considering do not undergo rotational diffusion and do not possess directional persistence. Thus, following [30], we see that for times preceding a run (when a bacterium is tumbling) the position of a bacterium is, when averaging over paths, not correlated to the velocity. On the other hand, during a run, the position correlates with velocity. In this way, in (A10), we have

$$ \langle x(t' - \theta) \rangle = \begin{cases} 0, & t' < \theta \\ \pm v(x)(t' - \theta), & t' > \theta \end{cases} \quad (A12) $$

So that integral (A10) becomes:

$$ I_1(t) = \pm \nabla f_\chi v(x) \int_0^t (t' - \theta) dt' = \pm \nabla f_\chi v(x) \frac{(t' - \theta)^2}{2} \quad (A13) $$

We proceed analogously to evaluate integral (A11). In this case, the integral over paths for times preceding a run does not average to zero, but sums up the contributions from the temporal variation of the gradient during tumbles. During a run, on the other hand, the temporal variation is evaluated for when bacteria are travelling up (down) the gradient, providing positive (negative) weights to the path integration, which provides

$$ \langle \partial_t f_\chi |_{x(t')} \rangle = \begin{cases} F(t', \theta), & t' < \theta \\ \pm \partial_t f_\chi, & t' > \theta \end{cases} \quad (A14) $$
where $F(t', \theta)$ is a function of time and delay. In this way, we can carry out integral (A11)

$$ 
\begin{aligned}
I_2(t) &= \int_0^\theta F(t', \theta)(t' - \theta) dt' \pm \partial_t f_\chi |_x \int_\theta^t (t' - \theta) dt' \\
&= G(\theta, t) \pm \partial_t f_\chi |_x \frac{(t' - \theta)^2}{2}
\end{aligned} 
\quad (A15) 
$$

where we have defined $G(\theta, t) = \int_0^\theta F(t', \theta)(t' - \theta) dt'$. Because it does not depend on direction up/down the gradient, this function (and its integrals) will cancel out in the evaluation of the run time, and will henceforth be omitted from the derivation. Inserting expressions (A13) and (A15) back into (A9), we can thus write the mean run durations up (R) and down (L) the gradient as

$$ 
\begin{aligned}
\tau^{L,R} \approx & \frac{1}{\alpha_0} + A\alpha_0 \int_0^\infty dt e^{-\alpha_0 t} \\
& \times \left[ v(x) \nabla f_\chi + \partial_t f_\chi \right] \frac{(t-\theta)^2}{2},
\end{aligned} 
\quad (A16) 
$$

which can be integrated to obtain

$$ \tau^{L,R} \approx \frac{1}{\alpha_0} \pm \frac{v(x)}{\alpha_0^2} \left[ \nabla f_\chi + \frac{1}{v(x)} \partial_t f_\chi \right] A e^{-\alpha_0 \theta}. \quad (A17) $$

Extending to general response function $K(\theta)$ with a distribution of delays [30], Eq. (A17)

$$ \tau^{L,R} \approx \frac{1}{\alpha_0} \pm \frac{v(x)}{\alpha_0^2} \left[ \nabla f_\chi + \frac{1}{v(x)} \partial_t f_\chi \right] \int_0^\infty d\theta K(\theta) e^{-\alpha_0 \theta}. \quad (A18) $$

Since $\alpha_{L,R} = 1/\tau_{L,R}$ and (see main text)

$$ V_{\chi} = v(x) \frac{\alpha_{L} - \alpha_{R}}{\alpha_{R} + \alpha_{L}}, \quad (A19) $$

we then arrive at an expression for the chemotactic drift

$$ V_{\chi} = \frac{v(x)^2}{\alpha_0} \beta \left[ \nabla f_{\chi} + \frac{1}{v(x)} \partial_t f_{\chi} \right], \quad (A20) $$

where we have defined the constant $\beta = \int_0^\infty d\theta K(\theta) e^{-\alpha_0 \theta}$. The standard definition of the chemotactic sensitivity is obtained from the empirical drift $V_\chi = \chi \nabla f_\chi$, where as previously $\chi$ is the chemotactic sensitivity. We extend here the definition to include temporal gradients, and define

$$ \chi(x) = \frac{v(x)^2}{\alpha_0} \beta, \quad (A21) $$

so that chemotactic drift is given by Eq. (5) in the main text. In the absence of chemokinetic alterations to the swimming speed, the sensitivity is simply

$$ \chi_0 = \frac{v_0^2}{\alpha_0} \beta, \quad (A22) $$

where subscripts denote a constant swimming speed. Dividing (A21) by (A22) provides relationship (8) in the main text.

<table><caption>**Table 1** Parameters used in simulations</caption><thead><tr><th>Parameter</th><th>Value</th><th>Figures</th><th>References</th></tr></thead><tbody><tr><td rowspan="2">H</td><td>0.0</td><td>3, 4</td><td rowspan="2">[29]</td></tr><tr><td>3.5</td><td>5, 7, 8</td></tr><tr><td rowspan="2">N</td><td>0.0</td><td>3, 4</td><td rowspan="2">[29]</td></tr><tr><td>0.5</td><td>5, 7, 8</td></tr><tr><td rowspan="2">δ<sub>0</sub></td><td>50.0</td><td>3, 4, 7</td><td rowspan="2">[29]</td></tr><tr><td>105.0</td><td>5</td></tr><tr><td rowspan="2">ω</td><td>0.2</td><td>3, 4, 7</td><td rowspan="2">[29]</td></tr><tr><td>0.5</td><td>5</td></tr><tr><td rowspan="2">η</td><td>2.0</td><td>3, 4, 7</td><td rowspan="2">[5]</td></tr><tr><td>0.5</td><td>5</td></tr><tr><td rowspan="4">n</td><td>1</td><td>4, 7</td><td rowspan="4">[29]</td></tr><tr><td>5</td><td>3, 4, 5</td></tr><tr><td>10</td><td>4</td></tr><tr><td>40</td><td>4</td></tr><tr><td rowspan="2">ζ</td><td>0.0</td><td>3, 4, 5</td><td rowspan="2">[29]</td></tr><tr><td>8.164 × 10<sup>-3</sup></td><td>7</td></tr><tr><td>K<sub>S</sub></td><td>1.0</td><td>3, 4, 5, 7</td><td>[29]</td></tr><tr><td>K<sub>χ</sub></td><td>0.53</td><td>3, 4, 5, 7</td><td>[29]</td></tr><tr><td>S</td><td>0.5</td><td>7 8</td><td></td></tr></tbody></table>

Please see main text for parameter definitions

## Appendix B: Numerical solution and parameter values

A finite difference scheme was chosen to compute the numerical solution of the system of PDEs (11). The domain $\Omega$ is represented by a vector of $M$ equally spaced grid points, $X_1, \dots, X_M$. A fourth-order scheme in space and a first order scheme in time are used to approximate the derivatives. Homogeneous Neumann boundary conditions are imposed for both substrate and bacterial concentration, i.e. $\frac{\partial C}{\partial X}|_{X=1,M} = 0$ and $\frac{\partial B}{\partial X}|_{X=1,M} = 0$.

In the case of steady attractant gradient (Figs. 3, 4), the initial attractant distribution is $C(X, 0) = 0.01X$, while the bacteria are uniformly distributed $B(X, 0) = 1.0$. For the self-generated gradient (Figs. 5, 6) we assume a Gaussian bacterial inoculum, $B(R, 0) = \exp(-R^2/\sigma^2)$, with a corresponding attractant distribution, $C(R, 0) = 1 - \exp(-R^2/\sigma^2)$ [29]. Finally, for Figs. 7 and 8, the bacteria are initially uniformly distributed $B(R, 0) = 0.2$, while the attractant distribution is $C(R, 0) = S(4\pi NT)^{-1} \exp(-R^2/4TN)$, where $S = 0.5$ and $T = 0.02$.

The parameters used in the non-dimensional model are summarized in Table 1. The parameters are mostly based on literature values for E. coli. Due to the lack of experimental data required for a full estimate of the functional form of the chemokinetic response of, e.g. a single species, the values for the chemokinetic response ($\eta, \omega, n$) are motivated by several studies and chosen to illustrate the response. For example, $\eta = 0.5$ corresponds to a 50% maximum increase in the swimming speed, which is on the order of magnitude that has been reported for several species [5, 12, 13]. The largest increase reported, to the best of our knowledge, is 80% for V. alginolyticus [15]. Thus, $\eta = 2$ is chosen as an extreme value to illustrate the effect of chemokinesis more clearly. The half-saturation constant $\omega$, which is required for the chemokinetic function given by Eq. (12), has not been reported directly. Upon (visual) inspection of results,
half-saturation constants may be about 0.1mM of glucose and 0.1mM of acetate for E.coli [13] and R. sphaeroides [12], respectively. In [5], the half maximum speed seems to be reached at about 30% relative mucus concentration, while the maximum is reached at about 60%. As pointed out in [28], absolute concentration above a threshold can be detected before gradients can be accurately measured. Therefore, we assumed that the half-saturation constant of chemokinetic response (ω) should be below or on the same order as the chemotactic half-saturation constant (Kₓ). As there has not been a functional fit to chemokinesis measurements, an estimate for the Hill factor *n* is difficult. From visual inspection of results in [5, 13], we may assume a factor between 1 and 3. Simulations with larger Hill factors were performed to compare to results based on the assumption of a step-change in swimming speed (i.e. *n* → ∞) as done in agent-based models [5, 15]. The parameters used to produce the figures in the main text were chosen to illustrate the effect of chemokinesis best. As a comparison, other parameter combinations are given in the Supplementary Material.

## References

1. H.C. Berg, *Random walks in biology* (Princeton University Press, 1993)

2. H.C.E. Berg, *coli in motion* (Springer, 2004)

3. M.J. Tindall, P.K. Maini, S.L. Porter, J.P. Armitage, Overview of Mathematical Approaches Used to Model Bacterial Chemotaxis II: Bacterial Populations. B Math Biol. **70**, 1570–1607 (2008)

4. G.M. Barbara, J.G. Mitchell, Marine bacterial organisation around point-like sources of amino acids. FEMS Microbiology Ecology. **43**(1), 99–109 (2003a)

5. M. Garren, K. Son, J.B. Raina, R. Rusconi, F. Menolascina, O.H. Shapiro et al., A bacterial pathogen uses dimethylsulfoniopropionate as a cue to target heat-stressed corals. ISME J. **8**(5), 999–1007 (2014)

6. J.P. Armitage, R. Schmitt, Bacterial chemotaxis: Rhodobacter sphaeroides and Sinorhizobium meliloti - variations on a theme? *Microbiology*. **143**, 3671–3682 (1997)

7. U. Attmannspacher, B. Scharf, R. Schmitt, Control of speed modulation (chemokinesis) in the unidirectional rotary motor of Sinorhizobium meliloti. *Molecular microbiology*. **56**(3), 708–718 (2005)

8. V.M. Meier, P. Muschler, B.E. Scharf, Functional analysis of nine putative chemoreceptor proteins in Sinorhizobium meliloti. J Bacteriol. **189**(5), 1816–1826 (2007)

9. V. Sourjik, R. Schmitt, Different roles of CheY1 and CheY2 in the chemotaxis of Rhizobium meliloti. Molecular microbiology. **22**(3), 427–436 (1996)

10. I.B. Zhulin, J.P. Armitage, Motility, chemokinesis, and methylation-independent chemotaxis in Azospirillum brasilense. J Bacteriol. **175**(4), 952–958 (1993)

11. S. Brown, P.S. Poole, W. Jezierska, J.P. Armitage, Chemokinesis in Rhodobacter sphaeroides is the result of a long term increase in the rate of flagellar rotation. Biochim Biophys Acta. **1141**, 309–312 (1993)

12. H.L. Packer, J.P. Armitage, The chemokinetic and chemotactic behavior of Rhodobacter sphaeroides: Two independent responses. J Bacteriol. **176**(1), 206–212 (1994)

13. D. Deepika, R. Karmakar, M.S. Tirumkudulu, K.V. Venkatesh, Variation in swimming speed of Escherichia coli in response to attractant. Arch Microbiol. **197**(2), 211–222 (2014)

14. M. Garren, K. Son, J. Tout, J.R. Seymour, R. Stocker, Temperature-induced behavioral switches in a bacterial coral pathogen. ISME J. **10**, 1363–1372 (2016)

15. K. Son, F. Menolascina, R. Stocker, Speed-dependent chemotactic precision in marine bacteria. Proc Natl Acad Sci USA **113**(31), 1–6 (2016)

16. J. Van Houten, Two mechanisms of chemotaxis in Paramecium. J Comp Physiol. **127**(3), 167–174 (1978)

17. J. Van Houten, D. Hauser, M. Levandowsky, Chemosensory behavior in protozoa. Biochemistry and physiology of protozoa. **4**, 67–124 (1981)

18. J. Van Houten, J. Van Houten, Computer simulation of Paramecium chemokinesis behavior. Journal of Theoretical Biology. **98**(3), 453–468 (1982)

19. C. Giuffre, P. Hinow, R. Vogel, T. Ahmed, R. Stocker, T.R. Consi et al., The ciliate Paramecium shows higher motility in non-uniform chemical landscapes. PloS one. **6**(4), e15274 (2011)

20. G.M. Barbara, J.G. Mitchell, Bacterial tracking of motile algae. FEMS Microbiology Ecology. **44**(1), 79–87 (2003b)

21. A. Walther, A.H. Müller, Janus particles: synthesis, self-assembly, physical properties, and applications. Chem Rev. **113**(7), 5194–5261 (2013)

22. J.R. Howse, R.A. Jones, A.J. Ryan, T. Gough, R. Vafabakhsh, R. Golestanian, Self-motile colloidal particles: from directed propulsion to random walk. Phys Rev Lett. **99**(4), 048102 (2007)

23. G.A. Jackson, Simulating chemosensory responses of marine microorganisms. Limnol Oceanogr. **32**(6), 1253–1266 (1987)

24. K. Son, J.S. Guasto, R. Stocker, Bacteria can exploit a flagellar buckling instability to change direction. Nature Phys. **9**(8), 494–498 (2013)

25. M.J. Schnitzer, Theory of continuum random walks and application to chemotaxis. Phys Rev E. **48**(4), 2553–2568 (1993)

26. H.G. Othmer, T. Hillen, The diffusion limit of transport equations II: Chemotaxis equations. SIAM Journal on Applied Mathematics. **62**(4), 1222–1250 (2002)

27. M.E. Cates, Diffusive transport without detailed balance in motile bacteria: does microbiology need statistical physics? Rep Prog Phys. **75**(4), 042601 (2012)

28. A.M. Hein, D.R. Brumley, F. Carrara, R. Stocker, S.A. Levin, Physical limits on bacterial navigation in dynamic environments. J Royal Soc Interface. **13**(114), 20150844 (2016)

29. O.A. Croze, G.P. Ferguson, M.E. Cates, W.C.K. Poon. Migration of chemotactic bacteria in soft agar: Role of gel concentration. Biophys J. **101**(3), 525–534 (2011)

30. P.G. De Gennes, Chemotaxis: The role of internal delays. Eur Biophys J. **33**(8), 691–693 (2004)

31. R. Rusconi, M. Garren, R. Stocker, Microfluidics Expanding the Frontiers of Microbial Ecology. Annu Rev Biophys. **43**(1), 65–91 (2014)

32. Othmer HG, Xue C. The mathematical analysis of biological aggregation and dispersal: progress, problems and perspectives. In: Dispersal, individual movement and spatial ecology. Springer; 2013. p. 79–127
33. J. Arlt, V.A. Martinez, A. Dawson, T. Pilizota, W.C. Poon, Painting with light-powered bacteria. Nat Commun. **9**(1), 768 (2018)

34. J. Adler, Chemotaxis in bacteria. Science. **153**, 708–716 (1966)

35. M. Seyrich, A. Palugniok, H. Stark. Traveling concentration pulses of bacteria in a generalized Keller–Segel model. New J Phys. 2019 oct;21(10):103001

36. J. Platzer, W. Sterr, M. Hausmann, R. Schmitt, Three genes of a motility operon and their role in flagellar rotary speed variation in Rhizobium meliloti. J Bacteriol. **179**(20), 6391–6399 (1997)

37. E. Martínez-García, P.I. Nikel, M. Chavarría, V. de Lorenzo, The metabolic cost of flagellar motion in Pseudomonas putida KT 2440. Environmental microbiology. **16**(1), 291–303 (2014)

38. S. Smriga, V.I. Fernandez, J.G. Mitchell, R. Stocker, Chemotaxis toward phytoplankton drives organic matter partitioning among marine bacteria. Proc Natl Acad Sci USA **113**(6), 1576–1581 (2016)

39. R. Stocker, J.R. Seymour, A. Samadani, D.E. Hunt, M.F. Polz, Rapid chemotactic response enables marine bacteria to exploit ephemeral microscale nutrient patches. Proc Natl Acad Sci USA **105**(11), 4209–4214 (2008)

40. D.R. Brumley, F. Carrara, A.M. Hein, Y. Yawata, S.A. Levin, R. Stocker, Bacteria push the limits of chemotactic precision to navigate dynamic chemical gradients. Proc Natl Acad Sci USA **116**(22), 10792–10797 (2019)

41. J.G. Mitchell, The energetics and scaling of search strategies in bacteria. The American Naturalist. **160**(6), 727–740 (2002)

42. B. Ni, R. Colin, H. Link, R.G. Endres, V. Sourjik, Growth-rate dependent resource investment in bacterial motile behavior quantitatively follows potential benefit of chemotaxis. Proceedings of the National Academy of Sciences. **117**(1), 595–601 (2020)

43. J.R. Taylor, R. Stocker, Trade-offs of chemotactic foraging in turbulent water. Science. **338**(6107), 675–679 (2012)

44. J.T. Locsei, Persistence of direction increases the drift velocity of run and tumble chemotaxis. J Math Biol. **55**, 41–60 (2007)

45. C. Xue, Macroscopic equations for bacterial chemotaxis: integration of detailed biochemistry of cell signaling. J Math Biol. **70**(1–2), 1–44 (2015)

46. J.E. Segall, S.M. Block, H.C. Berg, Temporal comparisons in bacterial chemotaxis. Proc Natl Acad Sci USA **83**(23), 8987–8991 (1986)

47. L. Xie, C. Lu, X.L. Wu, Marine bacterial chemoresponse to a stepwise chemoattractant stimulus. Biophysical journal. **108**(3), 766–774 (2015)