A study on the combined interplay between stochastic fluctuations and the number of flagella in bacterial chemotaxis

The chemotactic pathway allows bacteria to respond and adapt to environmental changes, by tuning the tumbling and running motions that are due to clockwise and counterclockwise rotations of their flagella. The pathway is tightly regulated by feedback mechanisms governed by the phosphorylation and methylation of several proteins. In this paper, we present a detailed mechanistic model for chemotaxis, that considers all of its transmembrane and cytoplasmic components, and their mutual interactions. Stochastic simulations of the dynamics of a pivotal protein, CheYp, are performed by means of tau leaping algorithm. This approach is then used to investigate the interplay between the stochastic fluctuations of CheYp amount and the number of cellular flagella. Our results suggest that the combination of these factors might represent a relevant component for chemotaxis. Moreover, we study the pathway under various conditions, such as different methylation levels and ligand amounts, in order to test its adaptation response. Some issues for future work are finally discussed.

💡 Research Summary

This paper presents a comprehensive mechanistic model of the bacterial chemotaxis signaling pathway, incorporating all known transmembrane receptors (MCPs), cytoplasmic proteins (CheA, CheW, CheY, CheZ, CheR, CheB), and their biochemical interactions such as phosphorylation, dephosphorylation, methylation, and demethylation. The authors implement a stochastic simulation of the pivotal response regulator CheY‑p using the τ‑leaping algorithm, which allows efficient handling of the wide range of reaction time scales present in the network while preserving the intrinsic noise of molecular events.

The simulation framework is used to explore how fluctuations in CheY‑p concentration interact with the number of flagella per cell. By varying the flagellar count from one to six (or more) in silico, the authors demonstrate that cells with few flagella exhibit large CheY‑p variance, leading to frequent and erratic switches between clockwise (CW) and counter‑clockwise (CCW) rotation. Consequently, their swimming trajectories become highly random, reducing chemotactic efficiency. In contrast, cells equipped with multiple flagella display an averaging effect: the torque generated by many independently rotating filaments smooths out CheY‑p noise, stabilizing the CCW “run” state and producing more directed movement.

The study further investigates the role of receptor methylation levels, controlled by the opposing activities of CheR (methyltransferase) and CheB (demethylase). High methylation suppresses CheA kinase activity, lowering the mean CheY‑p level and extending runs, whereas low methylation enhances CheA activity, raising CheY‑p and increasing tumble frequency. The authors systematically vary the CheR/CheB ratio and quantify both the mean and variance of CheY‑p, revealing that methylation not only shifts the average chemotactic bias but also modulates the amplitude of stochastic fluctuations.

Environmental perturbations are modeled by abrupt changes in external ligand concentration. Upon ligand addition, MCPs bind the attractant, quickly inhibiting CheA and causing a sharp drop in CheY‑p. The subsequent adaptation phase is driven by increased CheR activity, which methylates receptors and gradually restores CheA activity. The opposite response occurs when ligand is removed: CheB becomes active, demethylates receptors, and CheA activity recovers. Importantly, the adaptation dynamics differ with flagellar number: multi‑flagellated cells exhibit smoother, faster recovery because the collective torque mitigates transient CheY‑p spikes, whereas single‑flagellated cells experience pronounced overshoots and delayed return to baseline, impairing chemotactic performance.

Beyond these core findings, the paper discusses limitations and future directions. The current model treats flagella as independent torque generators and does not incorporate mechanical coupling, filament bundling, or hydrodynamic interactions that are known to affect swimming behavior. Additionally, cellular energy status (ATP levels), temperature effects, and spatial gradients of signaling components are omitted. The authors propose extending the framework to a multi‑scale model that couples stochastic intracellular signaling with continuum fluid dynamics and flagellar mechanics, and they suggest validation against single‑cell tracking experiments using fluorescence reporters for CheY‑p.

In summary, the work provides a novel quantitative link between stochastic intracellular signaling noise and a structural cellular parameter—flagellar count—highlighting how their combined interplay can shape chemotactic efficiency and adaptation. This integrative approach advances our understanding of bacterial navigation beyond deterministic models, offering a platform for future experimental and theoretical investigations into the robustness of microbial sensory systems.