Home | About | Journals | Submit | Contact Us | Français |

**|**PLoS Comput Biol**|**v.6(8); 2010 August**|**PMC2924242

Formats

Article sections

- Abstract
- Author Summary
- Introduction
- Results
- Discussion
- Materials and Methods
- Supporting Information
- References

Authors

Related links

PLoS Comput Biol. 2010 August; 6(8): e1000890.

Published online 2010 August 19. doi: 10.1371/journal.pcbi.1000890

PMCID: PMC2924242

Jonathan Saragosti,^{#}^{
1
,}^{
*
} Vincent Calvez,^{#}^{
2
,}^{
3
,}^{
*
} Nikolaos Bournaveas,^{
4
} Axel Buguin,^{
1
} Pascal Silberzan,^{
1
} and Benoît Perthame^{
5
,}^{
6
}

Christopher V. Rao, Editor^{}

University of Illinois at Urbana-Champaign, United States of America

Conceived and designed the experiments: JS AB PS. Performed the experiments: JS. Analyzed the data: JS VC NB AB PS BP. Contributed reagents/materials/analysis tools: VC NB BP. Wrote the paper: JS VC NB AB PS BP.

Received 2010 January 10; Accepted 2010 July 15.

Copyright Saragosti et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.

This article has been cited by other articles in PMC.

The Keller-Segel system has been widely proposed as a model for bacterial waves driven by chemotactic processes. Current experiments on *Escherichia coli* have shown the precise structure of traveling pulses. We present here an alternative mathematical description of traveling pulses at the macroscopic scale. This modeling task is complemented with numerical simulations in accordance with the experimental observations. Our model is derived from an accurate kinetic description of the mesoscopic run-and-tumble process performed by bacteria. This can account for recent experimental observations with *E. coli*. Qualitative agreements include the asymmetry of the pulse and transition in the collective behaviour (clustered motion versus dispersion). In addition, we can capture quantitatively the traveling speed of the pulse as well as its characteristic length. This work opens several experimental and theoretical perspectives since coefficients at the macroscopic level are derived from considerations at the cellular scale. For instance, the particular response of a single cell to chemical cues turns out to have a strong effect on collective motion. Furthermore, the bottom-up scaling allows us to perform preliminary mathematical analysis and write efficient numerical schemes. This model is intended as a predictive tool for the investigation of bacterial collective motion.

Modeling chemotaxis has raised a lot of interest in the applied mathematics community in past decades. The precise description of bacterial pulses traveling in a narrow channel is a challenging issue in the self-organization of cells. Indeed, our biological knowledge of signal integration in *E. coli* has grown in parallel with the development of more involved mathematical models. There exists a hierarchy of models for the analysis of bacteria *E. coli* motion depending on the scale under consideration. In this work, we derive macroscopic equations from the mesoscopic scale. This allows us to perform qualitative and quantitative analysis based on numerical simulations. We compare our predictions with current experiments performed with *E. coli*. The results can be interpreted at the cellular scale due to the bottom-up integration. This approach reveals better agreement with current experiments than the widely used Keller-Segel model. We conclude that the mesoscopic run-and-tumble description is compatible with the propagation of a pulse at the macroscopic scale.

Since Adler's seminal paper [1], several groups have reported the formation and the propagation of concentration waves in bacteria suspensions [2], [3]. Typically, a suspension of swimming bacteria such as *E. coli* self-concentrates in regions where the environment is slightly different such as the entry ports of the chamber (more exposed to oxygen) or regions of different temperatures. After their formation, these high concentration regions propagate along the channel, within the suspension. It is commonly admitted that chemotaxis (motion of cells directed by a chemical signal) is one of the key ingredients triggering the formation of these pulses. We refer to [4] for all biological aspects of *E. coli*.

Our goal is to derive a macroscopic model for these chemotactic pulses based on a mesoscopic underlying description. This approach relies on kinetic theory adapted to the specific run-and-tumble process that bacteria undergo [5], [6]. We base our modeling task on recent experimental evidence for traveling pulses obtained in our group (Fig. 1). These traveling pulses possess the following features which we are able to recover analytically: constant speed, constant amount of cells and strong asymmetry in the profile.

Many other micro-organisms exhibit collective behaviors. For instance, *Dictyostelium discoideum* cells collectively switch their cAMP-mediated signaling activity from stochastic to oscillatory when a concentration threshold is reached [7]. These oscillations, then synchronized at the scale of the population, give rise to non-dissipating waves of cAMP that guide the cells during fruiting body formation. Another example is given by *Myxococcus xanthus* that can grab and pull a neighbor cell by the mean of their pili, resulting in a long range alignment of the population and the formation of aggregates [8]. In the first case the pulsatile aspect is crucial for population scale communication and the speed of the cAMP waves is one order of magnitude larger than the velocities of the individual cells. In the second case, the physical contact between cells is critical for the aggregation.

Mathematical models for chemotaxis highlight a positive feedback which counteracts dispersion of individuals and may eventually lead to aggregation. There is a large amount of literature dealing with this subtle mathematical phenomenon (cf. [9], [10] and the references therein, see also [11] for alternative models which are closer to our approach). Self-induced chemotaxis following the Keller-Segel model has been shown successful for modeling self-organization of various cell populations undergoing aggregation [12]–[15].

In particular the Keller-Segel model has been proposed as a basis for modeling the propagation of traveling waves [16]–[19]. We refer to [20] for a complete review of contributions to this modeling issue. It has been postulated that a single chemotactic signal, namely the nutrient, could be responsible for the motion of the wave. However it is required that the chemosensitivity function is singular when the nutrient concentration vanishes. Our approach is more robust as we give a large class of fluxes for which traveling pulses do propagate. Furthermore these fluxes are derived from an accurate mesoscopic description of bacterial interactions.

In addition to chemotaxis, the contribution of cell division has been considered by many authors (cf. [21]–[23] and the references therein). Following the theory of reaction-diffusion equations, these authors have demonstrated the existence of traveling waves under general assumptions. However taking into account population growth seems unreasonable in view of the time scale of the experimental setting we aim at describing.

An extension of the classical Keller-Segel model was also proposed in seminal paper by Brenner *et al.*
[24] for the self-organization of *E. coli*. Production of the chemoattractant by the bacteria triggers consumption of an external field (namely the succinate). Their objective is to accurately describe aggregation of bacteria along rings or spots, as observed in earlier experiments by Budrene and Berg that were performed over the surface of gels [2]. However the experimental setting we are based on is quite different from Budrene and Berg's experiments: for the experiments discussed in the present paper, the bacteria swim in a liquid medium and not on agar plates. Therefore we will not follow [24]. On the other hand Salman et al. [25] consider an experimental setting very similar to ours. However the model they introduce to account for their observations is not expected to exhibit pulse waves (although the mathematical analysis would be more complex in its entire form than in [18]).

A new class of models for the collective motion of cells (*e.g.* swimming bacteria, the slime mold *D. discoideum*) has emerged recently. It differs significantly from the classical Keller-Segel model. Rather than following intuitive rules (or first order approximations), the chemotactic fluxes (, say: being the concentration of a chemotactic cue) are derived analytically from a mesoscopic description of the run-and-tumble dynamics at the individual level and possibly involving internal molecular pathways, [9], [26]–[34]. The upscaling limit which links the macroscopic flux to the kinetic description is now well understood since the pioneering works [5], [6], [26]. Here we propose to follow the analysis in [11], [31]. We write accordingly the macroscopic chemotactic flux in full generality as:

(1)

where denotes the concentration of chemoattractant. We shall derive an explicit formulation for the macroscopic quantity . Indeed it contains the microscopic features that stem from the precise response of a single bacterium to a change in the concentration of the chemoattractant in its surrounding environment. The upscaling limit is based on the following experimental fact: the (collective) pulse speed and the (individual) speed of bacteria differ by one order of magnitude. To the best of our knowledge, this is the first work where this powerful approach has been applied to the propagation of bands in populations of *E. coli*
[35].

When confined in micro environment, motile populations of *Escherichia coli* exhibit robust collective behaviours in the form of propagation of concentration waves. If this phenomenon is relatively easy to observe, its quantitative study requires a reproducible preparation of the system. To do so, we perform the following experiment. Fluorescent bacteria are grown in a nutritive medium until they reach a sufficient density and a good motility. We then fill PDMS/glass micro channels directly with this suspension, or after resuspension in a different medium. The channels are then sealed with epoxy resin thus confining the homogeneous suspension of motile bacteria. The centrifugation of this system reproducibly accumulate bacteria at one end of the channel while preserving the motility. When the centrifugation is stopped, a sharp pulse forms and propagates along the channel. Fluorescence video microscopy allows the measurement of the speed and shape of the traveling pulse. Precise experimental details are given in material and methods.

We describe the population of bacteria by its density (at time and position ). We consider here short timescales, hence cell division is assumed to be negligible. The cell density follows a drift-diffusion equation, combining brownian diffusion together with directed fluxes being the chemotactic contributions. This is coupled to reaction-diffusion equations driving the external chemical concentrations. In this paper we consider the influence of two chemical species, namely the chemoattractant signal , and the nutrient . Although this is a very general framework, it has been shown in close but different conditions that glycine can play the role of the chemoattractant [25]. Similarly, glucose is presumed to be the nutrient. The exact nature of the chemical species has very little influence on our modeling process. In fact there is no need to know precisely the mechanisms of signal integration at this stage. The model reads as follows:

(2)

The chemoattractant is assumed to be secreted by the bacteria (at a constant rate ), and is naturally degraded at rate , whereas the nutrient is consumed at rate . Both chemical species diffuse with possibly different molecular diffusion coefficients. We assume a linear integration of the signal at the microscopic scale, resulting in a summation of two independent contributions for the directed part of the motion expressed by the fluxes and . We expect that the flux will contribute to gather the cell density and create a pulse. The flux will be responsible for the motion of this pulse towards higher nutrient levels.

The fluxes and are built from the kinetic description of motion at the mesoscopic scale (see Materials and Methods). To summarize we assume that bacteria follow a run-and-tumble process mediated by the chemical micro-environment. The tumbling rate is dependent upon the material derivatives and (see [11], [31] for related works), where , and denotes the cell velocity. Namely, we assume that the tumbling rate writes as follows: . Here is the basal rate of tumbling in the absence of chemoattractant and is a decreasing function: tumble is more likely to occur if the chemoattractant concentration decreases along the trajectory [36], [37]. The (small) parameter accounts for the small variations of tumbling rates which have been measured experimentally (results not shown). The synthesis of these phenonena yields a macroscopic equation for the cell density (2), where the chemical drift is given by

(3)

where denotes the set of possible velocities. The same holds for . The dependency upon the time derivative disappears due to time/space scaling. We could keep this dependency at first order, but we omit it for the sake of clarity.

Several systems such as (2) have been proposed and the upmost classical is the so-called Keller-Segel equation [12], [16]. In the latter, the fluxes are proportional to the gradient of the chemical: , *resp.*
. Such a coupling is known to possibly drive the system into aggregated configurations for which the density of cells can become unbounded [9]. Notice that the two possible choices coincide in the linear regime, *i.e.* for small amplitudes of . They strongly differ however far from the linear regime. Especially the flux given by (3) is bounded by the individual speed of bacteria, whereas the chemotactic flux in the Keller-Segel model generally becomes unbounded when aggregative instability occurs, which is a strong obstacle to the existence of traveling pulses.

We restrict our attention to the one-dimensional case due to the specific geometry of the channels. It is usually impossible to compute explicitely traveling pulse solutions for general systems such as (2). To obtain qualitative properties is also a difficult problem: we refer to [17], [18], [23] for examples of rigorous results in this direction.

Here, we are able to handle analytical computations in the limiting case where the signal response function is indeed a step function. This owes to the assumption of high sensitivity of bacteria or large gradients of chemical. Then the fluxes (1) are given by the expression (3) which reduces to

(4)

We seek traveling pulses, in other words particular solutions of the form , , where denotes the speed of the wave. This reduces (2) to a new system with a single variable ,

(5)

We prescribe the following conditions at infinity

(6)

We impose without loss of generality. This means that the fresh nutrient is located on the right side, and thus we look for an increasing nutrient concentration . We expect that the chemoattractant profile exhibits a maximum coinciding with the cell density peak (say at ), and we look for a solution where changes sign only once at . Then, the fluxes (4) express under the traveling wave ansatz as

Integrating once the cell density equation in (5) we obtain

The flux takes two values (with a jump at ), whereas the flux is constant. Therefore the cell density is a combination of two exponential distributions

(7)

This combination of two exponentials matches with the numerical simulations (Fig. 2) and the experimental observations (Fig. 2).

To close the analysis it remains to recover the two unknowns: the maximum cell density and the speed , given the mass and the constraint that vanishes at (because reaches a maximum at this location). We have the following implicit formula for the speed of the pulse (see Text S1 for details):

(8)

We deduce from monotonicity arguments that there is a unique positive traveling speed .

On the other hand, the asymmetry factor is given by

(9)

This is a key macroscopic quantity as it enables to retrieve some parameters from experimental measurements. Interestingly enough, the speed and the asymmetry factor do not depend on the number of bacteria when the response function is stiff.

Mittal *et al.* have presented remarkable experiments where bacteria *E. coli* self-organize in coherent aggregated structures due to chemotaxis [38]. The cluster diameters are shown essentially not to depend on the quantity of cells being trapped. This experimental observation can be recovered from direct numerical simulations of random walks [39].

We can recover this feature in our analytical context using a model similar to (2) derived from a kinetic description. We compute the solutions of (5) in the absence of nutrient (assuming again a stiff response function). Observe that stationary solutions correspond here to zero-speed traveling pulses, that is

(10)

We assume again that . This simply leads to,

This is compatible with the postulate that changes sign only once, at (the source being even). The typical size of the clusters is of the order , which does not depend on the total number of cells. This is in good quantitative agreement with experiments exhibited in [38]. The fact that we can recover them from numerical simulations indicates that these stationary states are expected to be stable.

Cluster formation provides a good framework for investigating pattern formation when we relax the stiffness assumption on the response function . We introduce the stiffness parameter through its derivative at the transition between unfavourable and favourable regimes: . The case corresponds to a step response function.

We get from the dispersion relation (see Text S1) that the constant stationary state is linearly stable if and only if the following condition is fulfilled:

(11)

where the constant depends on the other parameters (including the mean square velocity , see Text S1). Here denotes the size of the channel and is the range of action of the chemical signal (namely ). The picture is not complete as we have not investigated the stability of the non-trivial steady-state. However this indicates that the stiffness parameter plays an important role regarding cluster formation. We show below that stiffness plays an important role for coherent motion of a pulse too.

We complete the theoretical analysis with some numerical simulations of the full model (2)–(3) exhibiting propagation of pulses (or not) in regimes where analytical solutions are not available (Fig. 3). The set of parameters is given in Table 1. The two parameters subject to variation are the stiffness parameter and the initial level of nutrient .

We can draw the following conclusions from our numerical simulations. The first remarkable fact is that we do observe traveling pulses (Fig. 3). Dispersion effects are counterbalanced by self-attraction due to the signal . These traveling pulses possess the correct asymmetry in the profile, and the speed matches experimental observations.

When the stiffness assumption for the internal response function is relaxed, so that dispersion effects become too strong, no pulse propagation is observed numerically (Fig. 4). This is in agreement with analytical results obtained for the zero-speed solution in the absence of nutriment. Indeed the cluster becomes unstable as gets large (11).

When the initial level of nutrient is low (or equivalently the consumption rate is high), and conditions for a pulse to travel are fulfilled, then only part of the bacterial population leaves the initial bump (Fig. 5). The solution appears to be the superposition of a traveling pulse and a stationary state (admissible in the absence of nutrient). Solitary modes with smaller amplitudes may appear at the back of the leading one (not shown). To predict which fraction of mass starts traveling turns out to be a difficult question.

We present in this article a simple mathematical description for the collective motion of bacterial pulses with constant speed and asymmetric profile in a channel. The nature of this model significantly differs from the classical Keller-Segel system although it belongs to the same class of drift-diffusion equations. Our model is formally derived from a mesoscopic description of the bacterial density, which allows for a more accurate expression of the cell flux at the macroscopic level.

The main conclusion of our work is the compatibility of the description of individual cell motion at the mesoscopic scale with the macroscopic observations of collective cell movements. The run-and-tumble process is qualitatively and quantitatively consistent with the propagation of a pulse at constant speed.

We do not discuss the conditions ensuring the existence of a traveling pulse solution to the system (5). This has been performed in the case of a stiff response function when computations are tractable analytically. We conjecture that existence holds under rather general assumptions, but such mathematical developments are beyond the scope of this work.

We point out the theoretical connection between the present work and the observations of Mittal *et al.*
[38]. The latter corresponds somehow to zero-speed traveling pulses, namely stationary clusters of bacteria. Our approach can be summarized as follows: a nutrient is added to pull chemotactic clusters of cells. This creates an imbalance in the fluxes which induces the asymmetry of the traveling profile.

We are able to compute the quantitative features of the traveling pulse in the case of a stiff response function. According to (8) the theoretical pulse speed does not depend upon the total number of cells. This can be related to experimental evidence by Mittal *et al.*
[38] where bacteria self-organize into size-independent clusters. In the case of a smooth tumbling kernel in (3), our model would predict a dependency of the speed upon the quantity of cells. But this analysis suggests that the number of cells is presumably not a sensitive biophysical parameter.

The speed also does not depend on the effective diffusion coefficient of bacteria when the response function is stiff. Therefore we expect to get the same formula if we follow the hyperbolic approach of [31] in order to derive a macroscopic model. Indeed the main difference is the diffusion coefficient which is very small in the hyperbolic scaling. Nevertheless, the density distribution would be very different, being much more confined when described by the hyperbolic system. Furthermore, scaling back the system to its original variables, we would obtain a pulse speed being comparable to the individual speed of bacteria (Materials and Methods). This is clearly not the case.

The asymmetry factor is another key outcome of the experimental observations. We are able to give a formula for this asymmetry when the response function is stiff. It turns out that asymmetry is favoured when is negligible with respect to the speed of the pulse (9). The former parameter is known as the propagation speed of a reaction-diffusion front [12], [40], except that the sign of is the opposite.

Although the chemotactic equation of (1)–(2) is significantly different from the standard Keller-Segel model, they coincide in the linear regime. It is well known that the Keller-Segel system is subject to a bifurcation phenomenon due to its quadratic nonlinearity [9], [10]. In the context of cluster formation, we learn from (11) that the stiffness parameter plays an important role in the stability of the homogeneous (flat) state. In other words, it is required that the bacteria are sufficiently sensitive in order to form a stable cluster. Clearly the same kind of mechanism acts here (Fig. 3 as opposed to Fig. 4). However there is no mathematical argumentation to sustain those numerical and intuitive evidence yet.

The influence of the stiffness property of the signal integration process is clear from numerical simulations of the full model (1)–(2). When the response function is smooth, dispersion effects are too strong and the population spreads out (Fig. 4). On the other hand, a stiff response function enables the cells to remain packed under the effect of the self-attractive chemical potential . Establishing the exact conditions that guarantee the propagation of a traveling pulse seems to be a challenging task.

Dynamics of the nutrient have no influence when the response function is stiff (only the sign of the gradient is important). However the evolution of the nutrient plays an important role when the response function is not stiff. It may happen than only part of the population starts traveling when the nutrient is initially at a low level (or is consumed with fast rate). A fraction remains trapped on the boundary, in a cluster configuration, while the rest of the population travels independently with constant speed (Fig. 5).

The next step consists in working at the kinetic level. Much has to be done for the design of efficient numerical methods for the description of collective motion of cells subject to chemotactic interactions. It would also be feasible to point out the dependency of the tumbling operator upon some internal variable (*e.g.* the cytoplasmic concentration of the phosphorylated form of the protein CheY, which is responsible for the reversal of motors). This approach carries out the coupling between an internal protein network and the external chemoattractant signals [39], [41]. Kinetic models are also relevant for describing this microscopic mechanism [30], [42] (the network is basically transported along the cells' trajectories). However the increase in complexity forces to reduce the size of the network, and to use rather caricatural systems mimicking high sensitivity to small temporal variations (excitation) and adaptation to constant levels of the chemoattractant.

Assuming independent integration of the chemical signals constitutes a strong hypothesis of our model. There exist two main membranous receptors triggering chemotaxis, namely Tar and Tsr. As the signals which act in the present experiments are not perfectly determined, we have considered the simplest configuration. To further analyse the interaction between the external signals, one should include more in-depth biological description of the competition for a single class of receptor [43].

We use the *E. coli* strain RP437 which is considered wild type for motility and chemotaxis. It is transformed with a pZE1R-gfp maintained by a resistance to ampicillin. The bacteria constitutively express a high level of green fluorescent protein which is necessary for low magnification fluorescence video microscopy.

We grow the bacteria with ampicillin on LB agar petri dishes at 37°C and keep them for a maximum of 5 days at 4°C. The “unlimited nutrient” culture medium is M9 supplemented with 4% D-Glucose, 1% Bacto Casamino Acids and . Before each experiment described here, a single colony is inoculated in of this medium (and ampicillin) and grown at 30°C under agitation until an OD600 of 0.5 is reached. We use falcon tube with two positions caps to make sure that oxygen is not limited during growth.

The fabrication of micro channels are based on usual soft lithography techniques [44]. high patterns are micro fabricated on silicon wafers using SU-8 100 resin from MICROCHEM. The PDMS is molded on the wafer and peeled after curing. A clean glass slide and the micro patterned PDMS are plasma treated for and directly placed in contact thereby forming a PDMS/glass micro channel. The result is a channel (widthheightlength) that is then filled with the homogeneous suspension of motile bacteria and sealed with a fast curing epoxy resin.

The glass silde is gently centrifuged (, from the axis) at room temperature for . The bacteria accumulate at the end of a channel and stay motile.

The channels is then placed in a chamber at constant temperature (30°C) under a Leica MZ16F stereomicroscope equipped for fluorescence. A CCD camera (CoolSnapHQ, Roper Scientific) records one image every of the fluorescence signal in the channel.

The movie is then processed using Matlab. We detect in each frame the position of the pulse by its maximum and extract its speed by fitting the successive positions by a linear regression.

The classical theory of drift-diffusion limit for kinetic modeling of bacterial chemotaxis is a way to compute the macroscopic fluxes , in (2) [27]. Because we assume a linear integration of the different signals for each individual, we restrict the following presentation to the action of a single chemical species .

The population of bacteria can be described at the mesoscopic scale by its local density of cells located at the position and swimming with velocity . The kinetic equation proposed in the pioneering works of Alt, Dunbar and Othmer [5], [6] combines free runs at speed , and tumbling events changing velocity from (anterior) to (posterior), according to the Boltzman type equation:

(12)

where the tumbling rate satisfies . The velocity space is bounded and symmetric, usually or (bacteria having presumably constant speed). As we deal with the idealization of a two-dimensional phenomenon in one dimension of space, we shall perform our computations for , but the results contained in this paper do not depend on this particular choice. Kinetic models of chemotaxis have been studied recently in [42], [45], [46].

The turning kernel describes the frequency of changing trajectories, from to . It expresses the way external chemicals may influence cell trajectories. A single bacterium is able to sense time variations of a chemical along its trajectory (through a time convolution whose kernel is well described since the experiments performed by Segall *et al.*
[37]). For the sake of simplicity we neglect any memory effect, and we assume that a cell is able of sensing the variation of the chemical concentration along its trajectory. Following [31], this is to say that is given by the expression

(13)

The signal integration function is non-negative and decreasing, expressing that cells are less likely to tumble (thus perform longer runs) when the external chemical signal increases (see Fig. 6 for such a tumbling kernel in the context of the present application). It is expected to have a stiff transition at 0, when the directional time derivative of the signal changes sign [37], [39], [41]. Our study in Section ‘Numerical insights’ boils down to the influence of the stiffness, by introducing a one parameter family of functions .

The main parameters of the model are the total number of bacteria which is conserved, the maximum speed of a single bacterium , and the mean turning frequency (where denotes the dimension of space according to our discussion above). The main unknown is the speed of the traveling pulse, denoted by . We rescale the kinetic model (12) into a nondimensional form as follows:

We aim at describing traveling pulses in the regime . Experimental evidence show that the bulk velocity is much lower than the speed of a single bacterium . This motivates to introduce the ratio . According to experimental measurements, we have . The kinetic equation writes:

(14)

where . Following the experimental setting (see Introduction, Fig. 1 and Fig. 2) and the biological knowledge [4], we choose the scales , , and . Hence . Therefore we rewrite this ratio as:

where the nondimensional coefficient is of order 1.

To perform a drift-diffusion limit when (cf. [9], [27], [28], [32], and [29], [31] for other scaling limits, *e.g.* hyperbolic), we shall assume that the variations of around its meanvalue are of amplitude at most. It writes in the nondimensional version as follows: . Hence the chemotactic contribution is a perturbation of order of a unbiased process which is constant in our case because the turning kernel does not depend on the posterior velocity and the first order contribution is required to be symmetric with respect to . This hypothesis is in agreement with early biological measurements. It is also relevant from the mathematical viewpoint as we are looking for a traveling pulse regime where the speed of the expected pulse is much slower than the speed of a single individual. This argues in favour of a parabolic scaling as performed here.

The resulting macroscopic equation writes as follows, with the position along the channel

(15)

Unlike the classical Keller-Segel model (used for instance by Salman et al. [25]), singularities cannot form (excessively populated aggregates) with the chemotactic flux given in (3). This is because the latter remains uniformly bounded (see also Mittal *et al.*
[38] where clusters emerge which are plateaus and thus not as singular as described for KS system in a mathematical sense).

We explain in the Text S1 how to derive the parabolic equation from the nondimensional kinetic equation (14). We arrive to equation (15) where the bacterial diffusion coefficient and the chemotactic flux are explicitely given by

(16)

In the limiting case where the internal response function is bivaluated: , the flux rewrites simply as

For the sake of comparison, we highlight the corresponding expressions which have been obtained by Dolak and Schmeiser. In [31] authors perform a hyperbolic scaling limit leading to the following chemotactic equation for the density of bacteria

where is an anisotropic diffusion tensor and the chemotactic flux is given by

for some renormalizing factor . The two approaches do not differ that much at first glance (in particular when is bivaluated). Notice however that the “small” parameter does not appear at the same location: in front of the diffusion coefficient in the hyperbolic limit and inside the chemotactic flux in the parabolic limit.

The macroscopic observable quantities are: the shape of the profile, namely the decay rates and , and the pulse speed . On the other hand, there are three parameters which we were unable to retrieve from the literature: the chemical degradation rate and the effective chemotaxis speeds and (although [25] indicates without reference). We deduce from the three constitutive relations (7), (8), the following formulas:

We get from experimental measurements the following values for the observable quantities: , and .

System (2) is solved using the MATLAB software. The drift-diffusion equation is discretized on a regular grid following a semi-implicit finite-difference scheme. The initial conditions are as follows: a decreasing exponential function centered on the left side of the channel for the cell density, no chemical signal, and a constant level of nutriment . The length of the computational channel is .

The Text S1 consists in three parts. First, we provide analytical details yielding Formulae (8) and (9) in section “Results”. Second, we perform the linear stability analysis referred to in section “Results”. Last, we perform the drift-diffusion limit which yields to equations (15)–(16) in Materials and Methods.

(0.07 MB PDF)

Click here for additional data file.^{(67K, pdf)}

This work is partly supported by the PIC “Maladies nosocomiales” – Institut Curie. We thank H C Berg, J S Parkinson for providing us with the strains used in this work, and C Beloin, J M Chigo for the plasmid. We are grateful to the anonymous reviewers for valuable comments.

The authors have declared that no competing interests exist.

This work is partly supported by the PIC “Maladies nosocomiales,” Institut Curie (www.curie.fr), and the Research team INRIA “NUMED,” INRIA Rhone-Alpes (www.inrialpes.fr). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

1. Adler J. Chemotaxis in bacteria. Science. 1966;153:708–716. [PubMed]

2. Budrene E, Berg H. Dynamics of formation of symmetrical patterns by chemotactic bacteria. Nature. 1995;376:49–53. [PubMed]

3. Park S, Wolanin P, Yuzbashyan E, Lin H, Darnton N, et al. Influence of topology on bacterial social interaction. Proc Natl Acad Sci USA. 2003;100:13910–15. [PubMed]

4. Berg H. *E. coli* in motion. Springer; 2004.

5. Alt W. Biased random walk models for chemotaxis and related diffusion approximations. J Math Biol. 1980;9:147–177. [PubMed]

6. Othmer H, Dunbar S, Alt W. Models of dispersal in biological systems. J Math Biol. 1988;26:263–298. [PubMed]

7. Gregor T, Fujimoto K, Masaki N, Sawai S. The onset of collective behavior in social amoebae. Science; 2010. (in press) [PMC free article] [PubMed]

8. Holmes A, Kalvala S, Whitworth D. Spatial simulations of myxobacterial development. PLoS Comput Biol. 2010;6 doi: 10.1371/journal.pcbi.1000686. [PMC free article] [PubMed]

9. Perthame B. Transport Equations in Biology. Frontiers in Mathematics. Birkhäuser Verlag; 2007.

10. Hillen T, Painter K. A user's guide to PDE models for chemotaxis. J Math Biol. 2009;58:183–217. [PubMed]

11. Mazzag B, Zhulin I, Mogilner A. Model of bacterial band formation in aerotaxis. Biophys J. 2003;85:3558–74. [PubMed]

12. Murray J. Mathematical biology, volume 2, Springer; 2002. 3rd edition.

13. Keller E, Segel L. Initiation of slime mold aggregation viewed as an instability. J Theor Biol. 1970;26:399–415. [PubMed]

14. Höfer T, Sherratt J, Maini P. Dictyostelium discoideum: Cellular self-organisation in an excitable medium. Proc Roy Soc Lond B. 1995;259:249–257. [PubMed]

15. Khonsari H, Calvez V. The origins of concentric demyelination: self-organization in the human brain. Plos ONE. 2007 doi: 10.1371/journal.pone.0000150. [PMC free article] [PubMed]

16. Keller E, Segel L. Traveling bands of chemotactic bacteria: A theoretical analysis. J Theor Biol. 1971;30:235–248. [PubMed]

17. Nagai T, Ikeda T. Traveling waves in a chemotactic model. J Math Biol. 1991;30:169–184. [PubMed]

18. Horstmann D, Stevens A. A constructive approach to traveling waves in chemotaxis. J Nonlinear Sci. 2004;14:1–25.

19. Li T, Wang Z. Nonlinear stability of traveling waves to a hyperbolic-parabolic system modeling chemotaxis. SIAM J Appl Math. 2009;70:1522–1541.

20. Tindall M, Maini P, Porter S, Armitage J. Overview of mathematical approaches used to model bacterial chemotaxis ii: bacterial populations. Bull Math Biol. 2008;70:1570–1607. [PubMed]

21. Landman KA, Simpson MJ, Slater JL, Newgreen DF. Diffusive and chemotactic cellular migration: smooth and discontinuous traveling wave solutions. SIAM J Appl Math. 2005;65:1420.

22. Funaki M, Mimura M, Tsujikawa T. Traveling front solutions arising in the chemotaxis-growth model. Interfaces Free Bound. 2006;8

23. Nadin G, Perthame B, Ryzhik L. Traveling waves for the Keller-Segel system with Fisher birth terms. Interfaces Free Bound. 2008;10:517–538.

24. Brenner M, Levitov L, Budrene E. Physical mechanisms for chemotactic pattern formation by bacteria. Biophys J. 1998;74:1677–93. [PubMed]

25. Salman H, Zilman A, Loverdo C, Jeffroy M, Libchaber A. Solitary modes of bacterial culture in a temperature gradient. Phys Rev Lett. 2006;97:118101. [PubMed]

26. Rivero MA, Tranquillo RT, Buettner HM, Lauffenburger DA. Transport models for chemotactic cell populations based on individual cell behavior. Chem Eng Sci. 1989;44:2881–2897.

27. Hillen T, Othmer H. The diffusion limit of transport equations derived from velocity-jump processes. SIAM J Appl Math. 2000;61:751–775.

28. Chalub F, Markowich P, Perthame B, Schmeiser C. Kinetic models for chemotaxis and their drift-diffusion limits. Monatsh Math. 2004;142:123–141.

29. Filbet F, Laurençot P, Perthame B. Derivation of hyperbolic models for chemosensitive movement. J Math Biol. 2005;50:189–207. [PubMed]

30. Erban R, Othmer H. From individual to collective behavior in bacterial chemotaxis. SIAM J Appl Math. 2004;65:361–391.

31. Dolak Y, Schmeiser C. Kinetic models for chemotaxis: hydrodynamic limits and spatio-temporal mechanisms. J Math Biol. 2006;51:595–615. [PubMed]

32. Chalub F, Dolak-Struss Y, Markowich P, Oelz D, Schmeiser C, et al. Model hierarchies for cell aggregation by chemotaxis. Math Models Methods Appl Sci. 2006;16:1173–97.

33. Erban R, Othmer H. Taxis equations for amoeboid cells. J Math Biol. 2007;54:847–885. [PubMed]

34. Xue C, Hwang H, Painter K, Erban R. Travelling waves in hyperbolic chemotaxis equations. Preprint 2010 [PubMed]

35. After submission of their manuscript, the authors have been aware of a similar work by Xue et al. they obtain analogous results based on an alternative approach.

36. Brown D, Berg H. Temporal stimulation of chemotaxis in Escherichia coli. Proc Natl Acad Sci USA. 1974;71 [PubMed]

37. Segall J, Block S, Berg H. Temporal comparisons in bacterial chemotaxis. Proc Natl Acad Sci USA. 1986;83:8987–91. [PubMed]

38. Mittal N, Budrene E, Brenner M, Oudenaarden AV. Motility of Escherichia coli cells in clusters formed by chemotactic aggregation. Proc Natl Acad Sci USA. 2003;100:13259–63. [PubMed]

39. Inoue M, Kaneko K. Conditions for self-consistent aggregation by chemotactic particles. Phys Rev E. 2008;77:041916. [PubMed]

40. Fisher R. The advance of advantageous genes. Ann Eugenics. 1937;7:335–369.

41. Spiro P, Parkinson J, Othmer H. A model of excitation and adaptation in bacterial chemotaxis. Proc Natl Acad Sci USA. 1997;94:7263–68. [PubMed]

42. Bournaveas N, Calvez V. Global existence for the kinetic chemotaxis model without pointwise memory effects, and including internal variables. Kinet Relat Models. 2008;1:29–48.

43. Keymer J, Endres R, Skoge M, Meir Y, Wingreen N. Chemosensing in Escherichia coli: two regimes of two-state receptors. Proc Natl Acad Sci USA. 2006;103:1786–91. [PubMed]

44. Xia Y, Whitesides G. Soft lithography. Annu Rev Mat Sci. 1998;28:153–184.

45. Hwang H, Kang K, Stevens A. Global solutions of nonlinear transport equations for chemosensitive movement. SIAM J Math Anal. 2005;36

46. Bournaveas N, Calvez V, Gutiérrez S, Perthame B. Global existence for a kinetic model of chemotaxis via dispersion and Strichartz estimates. Comm Partial Differential Equations. 2008;33:79–95.

47. Berg H. Random walks in biology. Princeton University Press; 1993. 2nd edition.

Articles from PLoS Computational Biology are provided here courtesy of **Public Library of Science**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |