A stochastic version of the Barkai-Leibler model of chemotaxis receptors in {it E. coli} is studied here to elucidate the effects of intrinsic network noise in their conformational dynamics. It was originally proposed to explain the robust and near-perfect adaptation of {it E. coli} observed across a wide range of spatially uniform attractant/repellent (ligand) concentrations. A receptor is either active or inactive and can stochastically switch between the two states. Enzyme CheR methylates inactive receptors while CheB demethylates active ones and the probability for it to be active depends on its level of methylation and ligandation. A simple version of the model with two methylation sites per receptor (M=2) shows zero-order ultrasensitivity (ZOU) akin to the classical 2-state model of covalent modification studied by Goldbeter and Koshland (GK). For extremely small and large ligand concentrations, the system reduces to two 2-state GK modules. A quantitative measure of the spontaneous fluctuations in activity (variance) estimated mathematically under linear noise approximation (LNA) is found to peak near the ZOU transition. The variance is a weak, non-monotonic and decreasing functions of ligand and receptor concentrations. Gillespie simulations for M=2 show excellent agreement with analytical results obtained under LNA. Numerical results for M=2, 3 and 4 show ZOU in mean activity; the variance is found to be smaller for larger M. The magnitude of receptor noise deduced from available experimental data is consistent with our predictions. A simple analysis of the downstream signaling pathway shows that this noise is large enough to have a beneficial effect on the motility of the organism. The response of mean receptor activity to small time-dependent changes in the external ligand concentration, computed within linear response theory, is found to have a bilobe form.
Download