Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2010 Apr 8;6(4):e1000735.
doi: 10.1371/journal.pcbi.1000735.

Quantitative modeling of Escherichia coli chemotactic motion in environments varying in space and time

Affiliations

Quantitative modeling of Escherichia coli chemotactic motion in environments varying in space and time

Lili Jiang et al. PLoS Comput Biol. .

Abstract

Escherichia coli chemotactic motion in spatiotemporally varying environments is studied by using a computational model based on a coarse-grained description of the intracellular signaling pathway dynamics. We find that the cell's chemotaxis drift velocity v(d) is a constant in an exponential attractant concentration gradient [L] proportional, variantexp(Gx). v(d) depends linearly on the exponential gradient G before it saturates when G is larger than a critical value G(C). We find that G(C) is determined by the intracellular adaptation rate k(R) with a simple scaling law: G(C) infinity k(1/2)(R). The linear dependence of v(d) on G = d(ln[L])/dx directly demonstrates E. coli's ability in sensing the derivative of the logarithmic attractant concentration. The existence of the limiting gradient G(C) and its scaling with k(R) are explained by the underlying intracellular adaptation dynamics and the flagellar motor response characteristics. For individual cells, we find that the overall average run length in an exponential gradient is longer than that in a homogeneous environment, which is caused by the constant kinase activity shift (decrease). The forward runs (up the gradient) are longer than the backward runs, as expected; and depending on the exact gradient, the (shorter) backward runs can be comparable to runs in a spatially homogeneous environment, consistent with previous experiments. In (spatial) ligand gradients that also vary in time, the chemotaxis motion is damped as the frequency omega of the time-varying spatial gradient becomes faster than a critical value omega(c), which is controlled by the cell's chemotaxis adaptation rate k(R). Finally, our model, with no adjustable parameters, agrees quantitatively with the classical capillary assay experiments where the attractant concentration changes both in space and time. Our model can thus be used to study E. coli chemotaxis behavior in arbitrary spatiotemporally varying environments. Further experiments are suggested to test some of the model predictions.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Figure 1
Figure 1. Illustrations of the E. coli chemotaxis pathway and the SPECS (Signaling Pathway-based E. coli Chemotaxis Simulator) model.
(A) The E. coli chemotaxis signal transduction pathway. The MCP complex receptor-CheW-CheA is the sensor and can be active (green) or inactive (dark brown). Binding of attractant molecules (yellow) decreases the probability of the receptor to be active. Once activated, the histidine kinase CheA quickly autophosphorylates and then transfers the phosphate group to CheY. CheY-p, the response regulator, diffuses to the flagellar motor and modulates its switching probability between CW and CCW rotations. CheZ, the CheY-p phosphatase, serves as the sink of the signal. Adaptation of the system is carried out by methylation and demethylation of the receptor, which are facilitated by the enzymes, CheR and CheB-p, respectively. (B) Flow chart of the SPECS model (reproduced from Figure 3 in [21]). The input for the signaling pathway (inside the box) is the instantaneous ligand concentration formula image. The internal signaling pathway dynamics is described at the coarse-grained level by the interactions between the average receptor methylation level formula image and the kinase activity formula image, which eventually determines the switching probability of the flagellar motor formula image. The switching probability is then used to determine the cell motion (run or tumble), and direction of motion during run fluctuates due to rotational diffusion. The motion of the cell leads it to a new location with a new ligand concentration for the cell and the whole simulation process continues.
Figure 2
Figure 2. Comparison of chemotactic motions in exponential and linear ligand profiles.
(A) Cell motion and intracellular dynamics in exponential ligand concentration profiles: formula image. (B) Cell motion and intracellular dynamics in linear ligand profiles: formula image. In both (A) & (B), the dynamics of the (population) averaged position (formula image); the average receptor methylation level (formula image) and the average kinase activity formula image are shown for different decay lengths formula image. formula image. The population-averaged position increases linearly with time until the methylation level reaches saturation in exponential profiles; while it slows down continuously in the linear profiles. After a transient decrease, the kinase activity stays roughly constant in exponential profiles, while it varies continuously recovering to its pre-stimulus level formula image in linear profiles. (C) Direct comparison of instantaneous velocities between exponential (solid lines) and linear (dotted lines) profiles for formula image (black); 1.0 mm (purple). Within the chemosensitivity range formula image, the instantaneous chemotaxis drift velocity is constant in exponential profiles, while it decreases continuously with [L] in the linear ligand profiles.
Figure 3
Figure 3. Dependence of chemotaxis motion on the adaptation rate.
(A) The chemotaxis drift velocity formula image for different exponential gradient formula image. Different symbols represent different adaptation rates formula image. Note that formula image first increases linearly (dashed line) with formula image before reaching a saturation velocity at a critical gradient formula image,. We can fit formula image with: formula image, in which formula image is the chemotaxis motility constant given by the linear fitting coefficient and the saturation drift velocity is formula image. The dependences of formula image and formula image on formula image are shown in (B) and (C) respectively. For the range of formula image we studied, we found that formula image is roughly independent of formula image and formula image depends on formula image with a simple scaling relation: formula image.
Figure 4
Figure 4. Single cell behaviors in exponential attractant gradients.
(A) Trajectory of a cell for 10 min for formula image with formula image. The random walk motion is biased towards the gradient (arrow). The forward runs up the gradient are in red, and the backward runs down the gradient are in black. (B) Run length distribution for forward runs (red), backward runs (black) in the presence of an exponential gradient and without a gradient (purple). The backward run length distribution is close to the run length distribution in the absence of a gradient, similar to the experimental results from Berg and Brown as reproduced in the inset. (C) Distributions of cell kinase activity for forward (red) and backward (black) runs. A time series of the activity formula image of a single cell is shown: forward is in red and backward is in black. The average activity for backward (black dashed line) is closer to the adapted activity (purple dashed line, 0.5) compared with the average activity for forward (red dashed line). (D) Methylation level of different individual cells at different times and in different exponential gradients (represented by color symbols as in Figure 2) all increase with logarithmic ligand concentration along a universal line, despite large temporal fluctuations in methylation levels and positions for (two) individual cells as shown in the inset.
Figure 5
Figure 5. Responses to oscillating linear gradient.
(A) Time dependence of the average positions of cells for three oscillatory gradients, all with the same amplitude but different frequencies (ω). The responses have the same frequencies as their driving signals, but the response amplitude decreases with the driving frequency. (B) The amplitude of the response decreases with frequency ω. The cross-over at low frequency (formula image) is caused by boundary effects. (C) Upon decreasing adaptation rate, a transition to a steeper decay of the amplitude appears at frequencies higher than a transition frequency formula image within the range of frequencies studied. Three cases with smaller values of formula image are shown, and the dependence of formula image on the adaptation rate formula image is shown in the inset of (C).
Figure 6
Figure 6. Quantitative simulation of the classical capillary assay and comparison with experiments.
(A) Time-dependent ligand concentrations at three different positions in the suspension pool (see inset) from directly solving the ligand diffusion equation. The ligand concentration at a given position peaks at a given time, depending on its location. formula image. (B) The exact ligand profile (solid line) along the center line of the capillary at different times, in comparison with the asymptotic solutions (dashed lines) by Furtelle and Berg . (C) Cell density in the rectangular coordinate is shown together with the contours of the logarithmic ligand concentration (in formula image) at different times. Three individual cell trajectories (starting from circles and ending at squares) are shown. Only the black cell ends in the capillary. (D) Probability distribution of the original positions of cells that end in the capillary. For a cell originally located at position formula image, the probability of it ending in the capillary at a later time (45 min), formula image, is shown. formula image, formula image. (E) Concentration-response curve for the capillary assay. The average number of bacteria in the capillary after 45–50 min subtracted by the number of bacteria in the capillary in the absence of attractant is defined as the response (ordinate). The results from our model with the exact ligand profile are labeled by solid symbols (fitted by a solid line). They agree well with the experimental measurements (hollow squares) of Mesibov et al ,. The results from using the asymptotic ligand profile by Furtelle and Berg are shown by the dashed lines. (F) Response curve for capillary assay with formula image. The solid symbols (fitted with a solid line) represent the model, and the hollow symbols represent the experimental results (both collected at 60 min).

Similar articles

Cited by

References

    1. Sourjik V. Receptor clustering and signal processing in E.coli chemotaxis. Trends Microbiol. 2004;12:569–576. - PubMed
    1. Stock J, Re SD. In: Chemotaxis. Lederberg J, Alexander M, Bloom B, Hopwood D, Hull R, editors. San Diego, CA: Academic Press; 2000. pp. 772–780.
    1. Bren A, Eisenbach M. How signals are heard during bacterial chemotaxis: protein-protein interactions in sensory propagation. J Bacteriol. 2000;182:6865–6873. - PMC - PubMed
    1. Berg HC, Tedesco PM. Transient response to chemotactic stimuli in Escherichia coli. Proc Natl Acad Sci USA. 1975;72:3235–3239. - PMC - PubMed
    1. Segall JE, Block SM, Berg HC. Temporal comparisons in bacterial chemotaxis. Proc Natl Acad Sci USA. 1986;83:8987–8991. - PMC - PubMed

Publication types