Optimal search in E. coli chemotaxis

Fluctuating landscapes and heavy tails in animal behavior

Animal behavior is shaped by a myriad of mechanisms acting on a wide range of scales. This immense variability hampers quantitative reasoning and renders the identification of universal principles elusive. Through data analysis and theory, we here show that slow non-ergodic drives generally give rise to heavy-tailed statistics in behaving animals. We leverage high-resolution recordings of C. elegans locomotion to extract a self-consistent reduced order model for an inferred reaction coordinate, bridging from sub-second chaotic dynamics to long-lived stochastic transitions among metastable states. The slow mode dynamics exhibits heavy-tailed first passage time distributions and correlation functions, and we show that such heavy tails can be explained by dynamics on a time-dependent potential landscape. Inspired by these results, we introduce a generic model in which we separate faster mixing modes that evolve on a quasi-stationary potential, from slower non-ergodic modes that drive the potential landscape, and reflect slowly varying internal states. We show that, even for simple potential landscapes, heavy tails emerge when barrier heights fluctuate slowly and strongly enough. In particular, the distribution of first passage times and the correlation function can asymptote to a power law, with related exponents that depend on the strength and nature of the fluctuations. We support our theoretical findings through direct numerical simulations.

Effect of receptor cooperativity on methylation dynamics in bacterial chemotaxis with weak and strong gradient

We study methylation dynamics of the chemoreceptors as an Escherichia coli cell moves around in a spatially varying chemoattractant environment. We consider attractant concentration with strong and weak spatial gradient. During the uphill and downhill motion of the cell along the gradient, we measure the temporal variation of average methylation level of the receptor clusters. Our numerical simulations we show that the methylation dynamics depends sensitively on the size of the receptor clusters and also on the strength of the gradient. At short times after the beginning of a run, the methylation dynamics is mainly controlled by short runs which are generally associated with high receptor activity. This results in demethylation at short times. But for intermediate or large times, long runs play an important role and depending on receptor cooperativity or gradient strength, the qualitative variation of methylation can be completely different in this time regime. For weak gradient, both for uphill and downhill runs, after the initial demethylation, we find methylation level increases steadily with time for all cluster sizes. Similar qualitative behavior is observed for strong gradient during uphill runs as well. However, the methylation dynamics for downhill runs in strong gradient show highly nontrivial dependence on the receptor cluster size. We explain this behavior as a result of interplay between the sensing and adaptation modules of the signaling network.

Effect of receptor clustering on chemotactic performance of E. coli: Sensing versus adaptation

We show how the competition between sensing and adaptation can result in a performance peak in Escherichia coli chemotaxis using extensive numerical simulations in a detailed theoretical model. Receptor clustering amplifies the input signal coming from ligand binding which enhances chemotactic efficiency. But large clusters also induce large fluctuations in total activity since the number of clusters goes down. The activity and hence the run-tumble motility now gets controlled by methylation levels which are part of adaptation module rather than ligand binding. This reduces chemotactic efficiency.

First-passage-time distribution in a moving parabolic potential with spatial roughness

In this paper, we investigate the first-passage-time distribution (FPTD) within a time-dependent parabolic potential in the presence of roughness with two methods: the Kramers theory and a nonsingular integral equation. By spatially averaging, the rough potential is equivalent to the combination of an effective smooth potential and an effective diffusion coefficient. Based on the Kramers theory, we first obtain Kramers approximations (KAs) of FPTD for both smooth and rough potentials. As expected, KA is valid only for high barriers and small external forces, and generally applicable for high barriers in rough potentials. To overcome the shortcoming of KA, a probability asymptotic approximation (PAA) based on an integral equation is proposed, which uses the transient probability density function (PDF) of the natural boundary conditions instead of the absorbing boundary conditions. We find that PAA fits very well even for large external forces. This enables us to analytically solve the FPTD for large external forces and low barriers as a strong extension to KA. In addition, we show that in the presence of a rough potential, the PAA of FPTD is in good agreement with numerical simulations for low barrier potentials. The PAA makes it possible to investigate the first-passage problem with ultrafast varying potentials and short exiting time. Thus, KA and PAA are complementary in determining the FPTD both for various barriers and external forces. Finally, the mean first-passage time (MFPT) is studied, which illustrates that the PAA of MFPT is effective almost in the whole range of external forces, while the KA of MFPT is valid only for small external forces.

A hydrodynamic-stochastic model of chemotactic ciliated microorganisms

Biological systems like ciliated microorganisms are capable of responding to the external chemical gradients, a process known as chemotaxis. In this process, the internal signaling network of the microorganism is triggered due to binding of the chemoattractant molecules with the receptors on the surface of the body. This can alter the activity at the surface of the microorganism. We study the chemotaxis of ciliated microorganisms using the chiral squirmer model, a spherical body with a surface slip velocity. In the presence of a chemical gradient, the coefficients of the slip velocity get modified resulting in a change in the path followed by the body. We observe that the strength of the gradient is not the only parameter which controls the dynamics of the body but also the adaptation time plays a very significant role in the success of chemotaxis. The trajectory of the body is smooth if we ignore the discreteness in the ligand-receptor binding which is stochastic in nature. In the presence of the latter, the path is not only irregular but the whole dynamics of the body changes. We calculate the mean first passage time, by varying the strength of the chemical gradient and the adaptation time, to determine the success rate of chemotaxis.

Run-and-tumble motion with steplike responses to a stochastic input

We study a simple run-and-tumble random walk whose switching frequencies between run mode and tumble mode depend on a stochastic signal. We consider a particularly sharp, steplike dependence, where the run-to-tumble switching probability jumps from zero to one as the signal crosses a particular value (say y_{1}) from below. Similarly, tumble-to-run switching probability also shows a jump like this as the signal crosses another value (y_{2}<y_{1}) from above. We are interested in characterizing the effect of signaling noise on the long-time behavior of the random walker. We consider two different time-evolutions of the stochastic signal. In one case, the signal dynamics is an independent stochastic process and does not depend on the run-and-tumble motion. In this case we can analytically calculate the mean value and the complete distribution function of the run duration and tumble duration. In the second case, we assume that the signal dynamics is influenced by the spatial location of the random walker. For this system, we numerically measure the steady state position distribution of the random walker. We discuss some similarities and differences between our system and Escherichia coli chemotaxis, which is another well-known run-and-tumble motion encountered in nature.

Optimal methylation noise for best chemotactic performance of E. coli

In response to a concentration gradient of chemoattractant, E. coli bacterium modulates the rotational bias of flagellar motors which control its run-and-tumble motion, to migrate towards regions of high chemoattractant concentration. Presence of stochastic noise in the biochemical pathway of the cell has important consequences on the switching mechanism of motor bias, which in turn affects the runs and tumbles of the cell in a significant way. We model the intracellular reaction network in terms of coupled time evolution of three stochastic variables-kinase activity, methylation level, and CheY-P protein level-and study the effect of methylation noise on the chemotactic performance of the cell. In presence of a spatially varying nutrient concentration profile, a good chemotactic performance allows the cell to climb up the concentration gradient quickly and localize in the nutrient-rich regions in the long time limit. Our simulations show that the best performance is obtained at an optimal noise strength. While it is expected that chemotaxis will be weaker for very large noise, it is counterintuitive that the performance worsens even when noise level falls below a certain value. We explain this striking result by detailed analysis of CheY-P protein level statistics for different noise strengths. We show that when the CheY-P level falls below a certain (noise-dependent) threshold the cell tends to move down the concentration gradient of the nutrient, which has a detrimental effect on its chemotactic response. This threshold value decreases as noise is increased, and this effect is responsible for noise-induced enhancement of chemotactic performance. In a harsh chemical environment, when the nutrient degrades with time, the amount of nutrient intercepted by the cell trajectory is an effective performance criterion. In this case also, depending on the nutrient lifetime, we find an optimum noise strength when the performance is at its best.

Twitching motility of bacteria with type-IV pili: Fractal walks, first passage time, and their consequences on microcolonies

A human pathogen, Neisseria gonorrhoeae (NG), moves on surfaces by attaching and retracting polymeric structures called Type IV pili. The tug-of-war between the pili results in a two-dimensional stochastic motion called twitching motility. In this paper, with the help of real-time NG trajectories, we develop coarse-grained models for their description. The fractal properties of these trajectories are determined and their influence on first passage time and formation of bacterial microcolonies is studied. Our main observations are as follows: (i) NG performs a fast ballistic walk on small time scales and a slow diffusive walk over long time scales with a long crossover region; (ii) there exists a characteristic persistent length l_{p}^{*}, which yields the fastest growth of bacterial aggregates or biofilms. Our simulations reveal that l_{p}^{*}∼L^{0.6}, where L×L is the surface on which the bacteria move; (iii) the morphologies have distinct fractal characteristics as a consequence of the ballistic and diffusive motion of the constituting bacteria.