Noise characteristics of the Escherichia coli rotary motor.  
Jump to Full Text  
MedLine Citation:

PMID: 21951560 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

BACKGROUND: The chemotaxis pathway in the bacterium Escherichia coli allows cells to detect changes in external ligand concentration (e.g. nutrients). The pathway regulates the flagellated rotary motors and hence the cells' swimming behaviour, steering them towards more favourable environments. While the molecular components are well characterised, the motor behaviour measured by tethered cell experiments has been difficult to interpret. RESULTS: We study the effects of sensing and signalling noise on the motor behaviour. Specifically, we consider fluctuations stemming from ligand concentration, receptor switching between their signalling states, adaptation, modification of proteins by phosphorylation, and motor switching between its two rotational states. We develop a model which includes all signalling steps in the pathway, and discuss a simplified version, which captures the essential features of the full model. We find that the noise characteristics of the motor contain signatures from all these processes, albeit with varying magnitudes. CONCLUSIONS: Our analysis allows us to address how celltocell variation affects motor behaviour and the question of optimal pathway design. A similar comprehensive analysis can be applied to other twocomponent signalling pathways. 
Authors:

Diana Clausznitzer; Robert G Endres 
Publication Detail:

Type: Journal Article; Research Support, NonU.S. Gov't Date: 20110927 
Journal Detail:

Title: BMC systems biology Volume: 5 ISSN: 17520509 ISO Abbreviation: BMC Syst Biol Publication Date: 2011 
Date Detail:

Created Date: 20111031 Completed Date: 20120518 Revised Date: 20130627 
Medline Journal Info:

Nlm Unique ID: 101301827 Medline TA: BMC Syst Biol Country: England 
Other Details:

Languages: eng Pagination: 151 Citation Subset: IM 
Affiliation:

Division of Molecular Biosciences, Imperial College London, SW7 2AZ London, UK. 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Chemotaxis Escherichia coli / physiology* Flagella / physiology* Ligands Models, Biological* Signal Transduction Systems Biology 
Grant Support  
ID/Acronym/Agency:

BB/G000131/1//Biotechnology and Biological Sciences Research Council 
Chemical  
Reg. No./Substance:

0/Ligands 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): BMC Syst Biol ISSN: 17520509 Publisher: BioMed Central 
Article Information Download PDF Copyright ©2011 Clausznitzer and Endres; licensee BioMed Central Ltd. openaccess: Received Day: 24 Month: 3 Year: 2011 Accepted Day: 27 Month: 9 Year: 2011 collection publication date: Year: 2011 Electronic publication date: Day: 27 Month: 9 Year: 2011 Volume: 5First Page: 151 Last Page: 151 ID: 3224245 Publisher Id: 175205095151 PubMed Id: 21951560 DOI: 10.1186/175205095151 
Noise characteristics of the Escherichia coli rotary motor  
Diana Clausznitzer123  Email: diana.clausznitzer@bioquant.uniheidelberg.de 
Robert G Endres12  Email: r.endres@imperial.ac.uk 
1Division of Molecular Biosciences, Imperial College London, SW7 2AZ London, UK 

2Centre for Integrative Systems Biology and Bioinformatics, Imperial College London, UK 

3BioQuant, Universität Heidelberg, 69120 Heidelberg, Germany 
Biological systems sense stimuli from their environment using cellsurface receptors, and process this information to make reliable decisions, e.g. where to move, or whether to divide or to express new enzymes. Typically, intracellular signalling molecules are activated by modification, e.g. phosphorylation or methylation, and interact in complicated biochemical reaction networks. The biochemical reactions in such networks rely on probabilistic collisions of a limited number of molecules. Hence, the number of signalling molecules fluctuates with time, making signal processing noisy. The abundance of noise sources in a cell is in stark contrast to the remarkable accuracy with which cells are known to respond to minute amounts of chemical concentration, including growing axons and immune cells [^{1},^{2}].
The high biological relevance of noise has widely been recognised and studied extensively in gene expression [^{3}^{}^{7}]. In contrast, noise in signal transduction is less well characterised, despite its importance for accurate sensing and celldecision making. Examples of eukaryotic systems, in which signalling noise has been considered include the ultrasensitive thresholding cascades [^{8},^{9}], pheromone sensing in yeast [^{10},^{11}], signal transduction in photoreceptors [^{12}], feedback loops for noise suppression [^{13},^{14}], and eukaryotic chemotaxis [^{15}]. Furthermore, signalling noise has been considered in parts of bacterial pathways [^{16}^{}^{19}]. However, such analyses have either been not comprehensive, or signal and noise transmission have not been compared in detail with experimental measurements along the pathway. An important class of signalling pathways are the bacterial twocomponent systems, including hundreds of pathways responsible for wide ranging functions such as sensing of and responding to nutrients, osmolarity, antibiotics, as well as quorum signals [^{20}]. A particularly well characterised example is the chemotaxis pathway in E. coli (Figure 1), allowing cells to swim towards nutrients and away from toxins with high sensitivity over a wide range of ambient concentrations [^{21}^{}^{25}]. Specifically, the kinase CheA autophosphorylates when receptors are active and passes on phosphoryl groups to the response regulators CheY and CheB. Phosphorylated CheY (CheYP) modulates the probability of counterclockwise (CCW) or clockwise (CW) rotation of the motor. The rotational directions of motors correspond to the two swimming modes of the bacterium, namely smooth swimming and tumbling, respectively. Adaptation, i.e. the reversal of the effect of changes in the ligand concentration, is mediated by reversible receptor methylation and demethylation, catalysed by enzymes CheR and phosphorylated CheB (CheBP), respectively. Using the E. coli chemotaxis pathway as an example, we are interested in the behaviour of the rotary motor, i.e. the cell's final output, and how its rotation is affected by signalling and noise.
To characterise signal propagation along the pathway, we use the socalled linear response function. This function characterises the response to a small timevarying input signal, such as impulse, sinusoidal and step stimuli. However, knowing the linear response function allows the calculation of system response to arbitrary input signals, provided they are sufficiently small (see Figure 2 and Methods). Typically, any system is subject to various sources of noise, i.e. random fluctuations in the input, as well as from signal processing. This is true in particular for biological systems, which rely on biochemical reactions and conformational changes of signalling molecules, which are intrinsically probabilistic. Here, we use the noise power spectrum to characterise fluctuations. The power spectrum captures the correlations between fluctuations in motor behaviour at different time points (see Figure 2 and Methods). In order to make the analysis easier, we typically consider the Fourier transforms of the dynamical equations of our model. The Fourier transform is an equivalent description of timedependent data, and represents its decomposition into its frequency components. Using this analysis, there are several fundamental questions we would like to address:
Firstly, what types of signals are transmitted and what types are attenuated by the pathway? Early work showed that the system responds to the timederivative of the input signal [^{26}]. A number of research groups have measured the averaged response of cells to chemotactic signals [^{26}^{}^{28}], and found that slowly, as well as rapidly changing input signals are not transmitted by the pathway. The response to slowly changing signals is attenuated by adaptation, which reverses the activation by ligand binding [^{28}^{}^{30}]. Rapidly changing signals were conjectured to be attenuated by a thirdorder filter [^{26},^{27}]. While the phosphorylation dynamics of CheYP has been shown to contribute a firstorder filter [^{30}], the exact filtering dynamics of the full pathway has not been addressed.
Secondly, how is noise generated, amplified or filtered in the signalling pathway, and how do different sources of noise affect the motor behaviour? The power spectrum of the motor rotation was measured for wildtype cells and mutant cells lacking the chemotaxis signalling pathway [^{31}]. The spectrum was found to have a large lowfrequency component in the wildtype cells, indicating that there is a dominant noise source in the signalling pathway with long correlations. In simulations of the signalling pathway, the adaptation dynamics can play an important role in generating long correlations [^{16},^{31}]. However, only signalling up to CheYP was considered. Other studies include stochastic simulations of the noisy biochemical reactions of the pathway [^{32}], and addressed the mechanism of motor rotation [^{33}^{}^{37}], including the thermodynamics of motor switching [^{38}^{}^{40}]. However, noise generation, filtering and amplification has not been addressed systematically for the various levels of the signalling pathway from chemoreceptors to motors.
Finally, how reliably are concentration input signals transmitted through the pathway in the presence of noise? An important task for the cell is to generate an appropriate motor response to input signals in the presence of fluctuations in the input, as well as due to noise in the biochemical signalling pathway. Furthermore, celltocell variation in protein expression influences signal transmission and noise filtering. Comparing these two aspects of the pathway dynamics, namely chemotactic signal transmission and noise filtering, in close comparison with experimental data along the pathway is a novel perspective in our present study.
In the following, we present a mathematical model for the chemotaxis signalling pathway. A simplified pathway is discussed in the text to provide intuitive understanding. However, results are shown for an extended model for the full pathway provided in the Methods section. We discuss the average (deterministic) response of the signalling pathway to concentration signals. We analyse the noise sources in the signalling pathway and their effects. Finally, we vary pathway parameters and study how they affect signal and noise transmission. We also discuss briefly how our approach can be applied to other twocomponent systems and signalling pathways.
The signal propagation in the chemotaxis pathway has been characterised by the response to small concentration signals (linear response function; see Methods). Specifically, the response has been measured at the level of CheYP using fluorescence resonance energy transfer (FRET) by Shimizu et al. [^{28}]. In that study the system was stimulated by a periodic variation of the concentration of attractant αmethylDLaspartate (MeAsp). Using a series of frequencies of the stimulation, the magnitude (modulus) and phase, i.e. the lag between signal and response, of the response was determined. In celltether experiments of motor rotation, the response to short impulses of attractants was measured at the level of the motor by Block et al. [^{26}] and Segall et al. [^{27}]. Such data determines the linear response function up to a constant factor. Experimental results are shown in Figure 3. Noise propagation in wildtype cells has only been studied at the level of the motor by Korobkova et al. [^{31}]. We use the experimental response functions to calibrate our model, and subsequently study the noise power spectrum and signaltonoise ratio.
Here we consider a simplified pathway to gain intuition of the key processes involved. The simplified pathway consists of chemoreceptor signalling in response to ligand binding and receptor methylation, as well as the rotary motor. Specifically, we use stochastic differential equations in a Langevin approach [^{41}] to describe the dynamics of each type of signalling protein. We assume throughout that fluctuations in concentration are small, allowing us to describe the average behaviour of a signalling molecule by a deterministic dynamics and the fluctuations around the mean by additive noise.
We assume N receptors form cooperative signalling complexes, which can switch between an active (on) and an inactive (off) state. Their activity A is described by the MonodWymanChangeux (MWC) model [^{42}^{}^{47}]. The activity depends on the external ligand concentration c at the receptor complex, as well as the methylation level M of the complex as detailed in Methods.
We consider N_{C }receptor complexes in a cell, and assume that each complex signals independently of the others. The total activity A_{c }of all receptors in a cell is determined by the sum over all signalling complexes j. The dynamics of the total activity is
(1)
dAcdt= ∑j=1NC∂A∂MdMjdt+∂A∂cdcjdt+ηAj(t), 
i.e. the dynamics of the complex activity is affected by changes in the receptor complex methylation level (first term), changes in ligand concentration (second term), as well as fluctuations due to the switching of the complex between its states (last term). All noise terms η(t) introduced in this section are discussed in Methods.
Changes in the concentration originate from timevarying input signals 〈c(t)〉, as well as fluctuations due to ligand diffusion. The dynamics of the concentration at the jth receptor complex is given by
(2)
dcjdt=d⟨c(t)⟩dt+ηcj(t), 
where the first term captures average concentration changes (indicated by angular brackets 〈...〉), affecting all receptors, and the second term describes concentration fluctuations at each receptor complex, assumed to be uncorrelated between different receptor complexes.
Adaptation is provided by reversible receptor methylation and demethylation, whose dynamics is described by the following equation [^{47}]:
(3)
dMjdt=γR(NAj)γBAj3+ηMj(t). 
The total methylation level M_{j }of a receptor complex j is changed by methylation of receptors in the inactive state by CheR (first term) and demethylation (second term). This latter rate is assumed to be strongly dependent on the receptor complex activity as only active receptors are demethylated by phosphorylated demethylation enzymes CheBP. These may act cooperatively, in agreement with time course data of wildtype and a CheB adaptation mutant [^{47}]. The last term represents fluctuations due to the noisy processivity of the methylation and demethylation enzymes.
The motor is described as a twostate system with CW and CCW rotating states, corresponding to running and tumbling modes, respectively. The dynamics of the motor bias X is derived from the Master equation of an ensemble of twostate motors, and is given by
(4)
dXdt=k+(Ac)(1X)k(Ac)X+ηX(t), 
with X the fraction of motors in the CW rotational state and applying the singlemotor limit. In Eq. 4, the first term represents the switching from CCW to CW with the transition rate k_{+}, the second term represents switching from CW to CCW with transition rate k_{}, and the third term describes temporal fluctuations in switching rates due to the stochastic nature of motor switching. Note that the Langevin equation for the twostate motor is exact [^{48}], see also Additional File 1.
Here, transition rates are modulated by the receptor signalling activity A_{c}, whereas in the full pathway model CheYP modulates motor switching. These rates have been experimentally measured using signalling mutants expressing varying amounts of constitutively active signalling molecule CheY [^{39}]. The switching rates, including a fit of the model we used [^{39}] (cf. Methods) to the data, are shown in Figure 3C.
We consider the response to input concentration signals at various levels in the signalling pathway to study how signals are transmitted to the rotary motor (see Methods). Briefly, an input signal Δc(t) is a concentration change relative to a constant background concentration c_{0}, affecting all receptors equally and representing a "meaningful" input to the chemotaxis signalling pathway. Hence, the concentration is given by 〈c(t)〉 = c_{0 }+ Δc(t). Furthermore, cells are assumed to be adapted to the prestimulus concentration c_{0 }with the various levels R of the signalling pathway adapted to their steadystate values R*.
We can analytically calculate the Fourier transformed linear response function from the dynamical equations Eq. 14 without noise (see Methods). We can analyse the filtering of the signal at each level of the pathway. The Fourier transformed linear response function for the total activity of all receptors in a cell is
(5)
χ^Ac(ω)=iωNC∂A∂cωMiω. 
The receptor activity is a highpass filter: The magnitude of the response function is small for frequencies ω below ω_{M }= (γ_{R }+ 3γ_{B}A*^{2})∂A/∂M, which is the characteristic frequency due to adaptation. For frequencies above ω_{M }the response function is a constant, given by the number of receptor complexes N_{C }participating in the response, and their sensitivity ∂A/∂c to ligand, evaluated at steadystate. The sensitivity is proportional to the receptor complex size N, i.e. it describes the amplification of the response of a single receptor.
Similarly, the Fourier transformed response of the motor is given by
(6)
χ^X(ω)=ω2ωXiωχ^Ac(ω). 
The motor is a lowpass filter, i.e. its dynamics introduces a frequencydependent response, which is constant below the characteristic frequency ωX=k+*+k* of the motor due to the steadystate switching rates k+* and k*. The parameter ω_{2 }describes the sensitivity of motor switching with respect to changes in receptor activity (Methods). At frequencies above ω_{X }the response is reduced. From Eq. 6 it is obvious that receptors and motor are in a cascade: The motor response introduces a new filter proportional to (ω_{i }iω)^{1 }which simply multiplies the response function of the response of the receptor activity. The response functions of the full pathway including the phosphorylation reactions are shown in Additional File 1.
For further analysis, we can write the Fourier transformed linear response function as
(7)
χ^R(ω)=χ^R(ω)eiϕR(ω), 
where χ^R is the magnitude and ϕ_{R }is the phase of the response function, which characterise the amplitude and lag of the response behind the input signal, respectively.
Figure 3 shows experimental data for the response function, as well as the fits of our full pathway model. Block et al. [^{26}] and Segall et al. [^{27}] measured the response of the motor using impulses of attractant. For our fit we adjusted adaptation and motor switching rates. Compared to the data by Shimizu et al. [^{28}] at the same temperature, adaptation rates are one order of magnitude larger, i.e. adaptation is faster in these experiments. The parameter ω_{X }of the motor switching is 2.1/s, consistent with switching rates of about 1 Hz [^{29}]. It is not clear from where the difference in adaptation rates between the two sets of experiments originates. However, different strains and media can lead to large variations in receptor expression level [^{49}]. Besides different experimental conditions, Shimizu et al. [^{28}] used populations of cells, whereas measurements by Segall et al. [^{27}] were done on single cells. For the fit of our model to the data by Shimizu et al. [^{28}], we adjusted only the adaptation rates, as measurements were restricted to low frequencies. The fit at 32°C yields the same adaptation parameters as obtained from fitting doseresponse curves of adapting cells [^{47}] (Figure 3B, left). The adaptation rates for room temperature are one order of magnitude smaller. Importantly, fitting to the magnitude of the Fourier transformed response yields a good fit for the phase of the response as well (Figure 3B, right). Fitted parameters are given in the Methods.
Figure 4 shows simulated time courses of the chemotactic response to an concentration impulse and the Fourier transforms of corresponding linear response functions, as well as our analytical results. As can be seen in the figure, linear response functions for the numerically solved nonlinear model indeed match the analytically calculated functions, confirming that our calculation results are valid for the chosen input signal. We observe how the input signal is transmitted through the pathway, with the effective pulse durations becoming progressively longer along the pathway (Figure 4, left), including total receptor activity in a cell (A_{c}), phosphorylated kinase CheA, phosphorylated response regulator CheY, and finally the motor (X). In Figure 4, middle we show the corresponding linear response functions.
The receptor activity acts as a highpass filter, i.e. it transmits highfrequency signals, but not lowfrequency signals. As can be seen from our simple model (cf. Eq. 5), this property is due to adaptation, which introduces the timederivative of the signal Δc(t) up to the characteristic frequency ω_{M}, eliminating the response to slowly changing attractant concentrations. The activity of chemoreceptors is the input to further levels in the pathway. The response of CheAP is fast, and shows no qualitative difference to the response of receptors in the frequency range shown. In contrast, due to the fast but finite rates of phosphorylation and dephosphorylation, preventing the CheYP concentration to respond to rapidly changing input signals, the response at the level of CheY is reduced at high frequencies. Similarly, the motor introduces another highfrequency filter due to slow switching between its two states. This additional filter can be deduced from Eq. 6, where the motor response function takes the response of chemoreceptors as input, and additionally introduces a characteristic cutoff frequency ω_{X }due to slow motor switching rates. Hence, the chemotaxis pathway acts as a bandpass filter [^{26}], which only transmits input signals within a selected frequency range, which is of the order of 1 to 10 s. This time scale corresponds to the average time between two tumbles, allowing sensing of concentration changes during periods of running. As shown in Figure 4, middle the phase tends towards π/2, i.e. a quarter period, at low frequencies. This has been analysed only for the receptor complex activity [^{28}]. This phase difference is due to adaptation and represents the fact that the system takes the time derivative of the stimulus below the characteristic frequency ω_{M }of adaptation. The phase shift of the receptor activity increases to π at high frequencies, indicating that the activity simply follows the output (a negative sign is due to the negative response of the activity to attractant concentration [^{28}]). The phase at high frequencies for the response of CheA follows the phase of the receptor activity, except for a small increase of the phase shift. In contrast, the phase of CheY and the motor increase significantly beyond π indicating that slow rates of modification and motor switching introduce a lag of the response behind the stimulus.
To understand the noise characteristics of the motor, we consider the noise sources and their transmission in the pathway. Each step in the signalling pathway is essentially probabilistic, hence, noisy: ligand diffusion and binding, receptor switching between its functional on and off states, as well as receptor methylation and demethylation, phosphorylation and dephosphorylation of signalling proteins CheA, CheY and CheB, and switching of the rotary motor between its two states, CW and CCW rotation. To characterise fluctuations of the phosphorylated signalling protein δR(t) around its mean value 〈R(t)〉, we use the power spectrum S_{R}(ω) and the variance 〈δR^{2}〉 = 〈R^{2}(t)〉  〈R(t)〉^{2 }(cf. Methods).
Considering Eq. 14 with noise, we can analytically calculate power spectra (see Methods). The power spectrum of activity fluctuations is given by
(8)
SAc(ω)=NCω2Sa(ω)+∂A∂c2Sc(ω)+∂A∂M2QMωM2+ω2. 
In this equation we considered fluctuations from receptor switching (first term in numerator), ligand diffusion (second term), as well as the receptor methylation dynamics (third term) at each of the N_{C }receptor complexes per cell. We have assumed that fluctuations at different receptor complexes are independent. Therefore, we obtain the sum of N_{C }identical spectra for all complexes. The individual terms S_{a}(ω), S_{c}(ω) and Q_{M }are given by Eq. 30, 34 and 36 in Methods. The frequency dependence of the ligand noise, as well as noise from receptor complex switching, indicates filtering of slowly varying fluctuations with frequencies below the characteristic frequency ω_{M }due to adaptation. In contrast, only highfrequency fluctuations from the receptor methylation dynamics are filtered by the adaptation dynamics. This is due to finite rates of methylation and demethylation fluctuations introducing correlations in the receptor methylation level.
The power spectrum of fluctuations in the motor bias is obtained as
(9)
SX(ω)=ω22SAc(ω)+QXω2+ωX2. 
The first term represents transmitted noise from receptor complexes, including the noise power spectrum of the receptor activity and the sensitivity ω22 of motor switching rates to changes in activity. The second term is motor switching noise. Both noises are filtered by the motor, as its finite rates of switching introduce correlations with characteristic frequency ω_{X}. The noise spectra of the full pathway including the phosphorylation reactions are shown in Additional File 1.
In Figure 4, right we show the power spectrum of fluctuations at the various levels of the signalling pathway, i.e. total receptor activity, CheAP, CheYP and the motor. We also plot the individual contributions from processes generating noise, namely ligand diffusion, receptor switching, methylation and demethylation of receptors, and phosphorylation and dephosphorylation of proteins, as well as motor switching. This allows us to follow how noise is generated and transmitted at the various levels of the pathway. The noise spectrum of the receptor activity has its largest contribution at low frequencies, which originates in the receptor methylation and phosphorylation dynamics. Most of the fluctuations from phosphorylation stem from CheB (the separate contributions to the phosphorylation noise are not shown in Figure 4, right). At highfrequencies, the activity noise spectrum is at. This is due to ligand and receptor switching noise, which is removed at low frequencies by adaptation, but not at highfrequencies. The general behaviour of the noise spectrum corresponds to the simplified model (cf. Eq. 8).
The noise spectrum of CheAP has generally the same shape as the activity spectrum with a large lowfrequency component, mainly due to receptor methylation and CheB phosphorylation dynamics. This spectrum also has an almost flat highfrequency behaviour in the frequency range shown. Apart from ligand and receptor switching noise, the flat part of the spectrum is largely determined by fluctuations from CheA autophosphorylation, which has roughly the same shape as activity noise at high frequencies because autophosphorylation depends on the receptor activity.
The noise spectrum of CheYP is also largest at low frequencies. However, at high frequencies the spectrum falls off as noise is filtered due to the finite rates of CheY phosphorylation and dephosphorylation, which introduce correlations in the fluctuations.
The motor introduces another layer of filtering of transmitted noise with the characteristic motor switching frequency ω_{X }(cf. Eq. 9). Hence, transmitted noise is reduced by two filters in the frequency range shown, namely due to the CheYP and motor dynamics. However, the main contribution to the spectrum is due to the motor switching itself, which is reduced only by a firstorder filter with characteristic frequency ω_{X}.
How are the signal response, fluctuations and the signaltonoise ratio (SNR) affected by changing parameters of the pathway such as size of receptor complexes, protein concentrations and reaction rate constants? In this section, we discuss the effect of celltocell variation on the power spectrum of the motor. In the next section, we discuss the SNR and its contributions, and how they depend on receptor complex size and adaptation rates.
According to our model parameters obtained from fitting the Fourier transformed linear response to data, the main contribution to the power spectrum comes from the steadystate switching of the motor between CCW and CW state. However, celltocell variation in protein content and motor switching rates can lead to modifications of the largely Lorentzianshaped spectrum. These modifications are caused by the transmitted noise from receptor methylation and phosphorylation dynamics (green and blue lines in Figure 4, right). Specifically, Figure 5A shows the motor power spectrum for increased motor switching rates as well as reduced adaptation rates and number of chemoreceptors in a cell. In all cases the lowfrequency component of the transmitted noise becomes more prominent.
An increased lowfrequency component has been observed in the motor power spectrum for cells with low motor bias [^{31}]. Both, wildtype cells and mutants lacking the signalling pathway were measured. Hence, the mutant's spectrum represents the component to the power spectrum from steadystate motor switching only. Wildtype cells showed a large lowfrequency component compared to the mutants. Figure 5B shows that our model can reproduce these experimental data (shown in the Inset), provided we assume a low CCW to CW switching rate leading to small motor bias (see Methods for the details of parameters). Specifically, the lowfrequency component of our spectrum originates from noise in the methylation and phosphorylation dynamics in cells with low motor bias. Furthermore, experiments show that the lowfrequency component in the motor power spectrum is reduced by increasing the expression level of CheR [^{31}]. In Figure 5C we show that our model reproduces this experimental finding (shown in the Inset). The lowfrequency component due to noise from receptor methylation and demethylation is effectively reduced by increasing the methylation rate constant.
To characterise how signals are transmitted in the presence of noise, we define the SNR at the level of the motor as
(10)
SNR=ΔX2⟨δX2⟩ 
with ΔX^{2 }and 〈δX^{2}〉 defined in Methods. For optimal signalling this ratio should be maximised. For simplicity, we only discuss the receptor activity in the text, while in the figures we additionally show the contribution from phosphorylation processes as transmitted to the motor in the full pathway model.
Receptor complexes amplify small signals proportionally to their size N. However, also concentration fluctuations are expected to be amplified. Hence, we hypothesise that the receptor complex size could be optimised to yield a balance of advantageous amplification of signals and detrimental amplification of input noise.
In Figure 6A we show the integrated motor response ΔX^{2 }(see Methods) to a step stimulus for varying background concentration and receptor complex size. We assume that the step stimulus size is a constant fraction of 10 percent of the background concentration. The integrated response has a characteristic variation with background concentration with the maximum in the sensitivity range of Tar receptors (indicated by their dissociation constants). Furthermore, the response increases with receptor complex size N. We calculated the integrated signal response of the receptor activity (see Methods). This quantity scales linearly with receptor complex size, ΔAc2∝N, due to coherent addition of the signalling responses of different receptor complexes, amplification of concentration changes by receptor complexes, as well as filtering by adaptation.
In Figure 6B and 6C, we show the variance (i.e., the integrated noise power spectrum, see Methods) of the transmitted noise of the pathway at the level of the motor. Only the contribution to the variance from ligand diffusion depends on the background concentration. Compared to the signal response, the maximum of the variance is shifted to a slightly lower concentration. The contribution to the variance from switching of receptor complexes is relatively small compared to the other contributions and roughly constant with receptor complex size, whereas those from ligand diffusion, receptor methylation and phosphorylation dynamics increase with receptor complex size.
To understand these behaviours of the variance more intuitively, we analysed the receptor activity analytically (for details of the calculation, see Methods). We find, the contribution to the variance of the receptor activity from receptor switching is indeed constant, independent of N. The contribution from ligand diffusion scales steeply as N^{2}, the difference between ligand noise and ligand signal amplification being due to (i) noise from different complexes is added up incoherently, and (ii) the main contribution to the variance coming from highfrequency ligand noise, which is not filtered by adaptation. The contribution from receptor methylation grows approximately linearly with receptor complex size as a result of the incoherent addition of fluctuations at different receptor complexes and the sensitivity of the receptor complex activity with respect to changes in methylation level increasing proportionally with N. The contribution to the variance from phosphorylation processes grows with receptor complex size similar to the contribution from the methylation dynamics. Overall, the total variance of transmitted noise at the level of the motor has contributions from receptor switching, the dynamics of receptor methylation, and phosphorylation. The latter is approximately constant or grows slower than the amplified signal response, whereas the component from ligand diffusion increases steeper than the signal response with growing receptor complex size.
The resulting SNR, i.e. the ratio of integrated signal response and variance of the noise, is shown in Figure 6D and 6E. The SNR is largest at background concentrations in the sensitivity range of the Tar receptor. Furthermore, due to the different dependencies of the signal and the noise on the receptor complex size, the SNR has a maximum at a particular receptor complex size (Figure 6E). The SNR grows below that complex size due to signal amplification, while the amplified ligand noise from ligand diffusion is still below the internal noise level from receptor switching and receptor methylation and phosphorylation dynamics. Above the optimal receptor complex size, the SNR decreases because the ligand noise is amplified more than the signal.
As shown above, adaptation filters slow input signals, with the adaptation speed determining what input frequencies are transmitted by the pathway. Furthermore, the adaptation dynamics filters input noise. Hence, adaptation rates may be expected to be optimised for signal and noise propagation.
Figure 7A shows the integrated signal response at the level of the motor for varying rates of receptor methylation (γ_{R}) and demethylation (γ_{B}). Varying these parameters describes changing the concentrations of receptor modification enzymes CheR and CheB. The integrated signal response is found to be maximal when both rates of methylation and demethylation become small simultaneously. Then, the adapted receptor activity is in the steep region of doseresponse curves, and the initial response to attractant is maximal. At the same time adaptation becomes slow and, therefore, the response lasts long. Interestingly, varying the two parameters independently has different effects on the signalling response: the integrated signal response increases for vanishing γ_{R}, whereas it decreases for vanishing γ_{B}. There are two effects that contribute to this behaviour, illustrated in Figure 7B: Firstly, if the concentration of one of the receptor modification enzymes is reduced, the receptors becomes modified predominantly by the opposing enzyme, hence driving the receptor activity towards saturation (A* = 0 or A* = 1). This effect would tend to quench the response by receptors. Secondly, as the enzyme concentration is reduced, adaptation times increase. Hence, this effect increases the integrated signal response as the time the receptor activity deviates from the adapted state increases. According to calculations shown in Methods for the integrated response of receptors, the first effect dominates in the case of reduced γ_{B}: Due to the strong activity dependence of the demethylation rate, reducing the demethylation rate constant effects the adapted activity of receptors strongly. Hence, receptors are quickly driven into saturation for vanishing γ_{B}. In contrast, in the case of reduced γ_{R }the second effect dominates and the increased adaptation time leads to an increased integrated signal response. At large methylation and demethylation rates, adaptation times are reduced leading to a decreasing integrated signal response.
The variance of fluctuations is shown in Figure 7C and 7D. The individual contributions from transmitted noise at the level of the motor look qualitatively similar. For γ_{R }and γ_{B }such that the adapted receptor activity is in the steep region of doseresponse curves, the variance of transmitted noise is largest. In contrast, for either vanishing γ_{R }or γ_{B }all contributions decrease, consistent with calculations for the variance of the receptor activity in Methods. In these cases, the adapted receptor activity becomes saturated, hence, quenching fluctuations transmitted by receptors.
The resulting SNR is shown in Figure 7E. The SNR increases for vanishing γ_{R }and decreases for vanishing γ_{B}. According to Figure 7E, a large SNR is obtained for small γ_{R }and large γ_{B}, corresponding to the parameters of our model.
Park et al. [^{50}] presented the idea that the signalling response to concentration signals and fluctuations in the chemotaxis pathway are not independent of each other, because they are produced by the same molecular interactions. Specifically, based on measurements at the level of the motor these authors proposed a fluctuationresponse theorem, namely an approximate linear relationship between the adaptation time to step stimuli and the variance of fluctuations in CheYP concentration.
Using our model, we tested this hypothesis and varied the adaptation rates, as well as the total CheY concentration in a cell, resulting in a shifted adapted CheYP concentration at steady state. We find that the variance of CheYP (normalised by the squared adapted value) decreases as the adapted CheYP value increases except for very small adapted CheYP concentrations (Figure 8A), indicating that the relative strength of fluctuations decreases as expected. In Figure 8B we show the adaptation time, approximated by the inverse of the characteristic frequency due to adaptation, plotted against the variance of CheYP. We find that at low adaptation times (thick line styles), the adaptation time increases with the variance of CheYP, indicating that cells with large fluctuations also respond longer to concentration signals. In contrast at long adaptation times, the adaptation time decreases with increasing variance of the pathway (grey parts of the curves). This behaviour can be directly traced back to the nonmonotonic variance shown in Figure 7. It is maximal when the adapted CheYP concentration is about 5 μM, i.e. when typically half of CheY is phosphorylated. The exact relationship depends on what parameters varied, exemplified by the different curves in Figure 8B. For each parameter and small adaptation times, we find an approximate linear relationship in line with Park et al. [^{50}], see Inset.
Biological signalling pathways employ biochemical reaction networks and molecular state transitions to sense and process signals from the environment. Fluctuations inherent in these processes determine the signals which can reliably be transmitted. Here, we studied the signal and noise propagation in the E. coli chemotaxis signalling pathway, which controls the bacterial swimming behaviour in chemical gradients. Specifically, we considered the dynamics of ligand diffusion, receptor methylation and demethylation, receptor complex switching between on and off, phosphorylation and dephosphorylation of the kinase CheA, and response regulators CheY and CheB, as well as from rotary motor switching between CW and CCW direction. We assume cooperative chemoreceptor signalling complexes, whose activities depend on ligand concentrations and receptor methylation level, described by the MWC model [^{42}^{}^{47}].
We formulated a model which includes all processes in the signalling pathway. Not included is the dynamics of gene expression and enzyme localisation to the chemoreceptor cluster, which is assumed to be much slower than the dynamics of signalling processes (discussed below). To make results intuitive we also presented a simplified version of the model, which only includes the dynamics of the activity of chemoreceptors, ligand concentration and receptor methylation level, as well as the motor dynamics. To calibrate the model, we first collected experimental data sets for the signalling pathway and rotary motor [^{26}^{}^{28}], and the motor switching behaviour [^{31},^{51}]. Using the Fourier transformed linear response function, we subsequently fitted our model parameters. We found a range of parameters fitting different data sets, revealing a striking experimental variation, which may require further characterisation in the future.
Despite the fitting, there is a discrepancy of our response function and the data at large frequencies. Block et al. [^{26}] and Segall et al. [^{27}] conjectured that the pathway is a thirdorder lowpass filter. In contrast, we find that the only relevant filters in that frequency range are due to CheYP and motor dynamics, leading to only a secondorder filter. One explanation for the missing filter is that experimental concentration pulses were not short enough, leaving a signature from the input signal at large frequencies. Alternatively, additional processes such as a slow release of CheYP from the chemosensory complexes as discussed in Ref. [^{52}] could lead to an additional filter. However, CheYP/CheZ complex formation and potential oligomerisation of CheYP/CheZ complexes [^{53}^{}^{55}] are not expected to contribute to highfrequency filtering (Additional File 1).
The motor behaviour is the final cell output, which contains characteristic noise signatures of all upstream signalling components, including the receptors. We found that motor switching is the dominant contribution to the spectrum of the fluctuations in motor bias (Figure 5). This may be not surprising as motor switching enables E. coli to tumble and change its swimming direction, and is therefore crucial for its search strategy. However, lowfrequency contributions from signalling and adaptation may become more significant for some cells or under certain conditions (Figure 5). Specifically, we predict that due to celltocell variation of protein contents or fast motor switching, these lowfrequency components become comparable to motor switching noise. For instance, Korobkova et al. [^{51}] measured power spectra in cells with low motor bias and found that the lowfrequency component is increased compared to mutants lacking the signalling pathway, and that increasing the expression level of CheR decreases the lowfrequency component. Our model is able to reproduce these spectra, demonstrating the general applicability of our approach. Long correlations in motor bias may enable subpopulations of cells to swim further without tumbling or to tumble more frequently.
To test our predictions, adaptation rates can be varied using cells expressing different amounts of CheR and CheB from an inducible plasmid. Alternatively, the natural variability in protein expression between cells can be exploited. Numerous wildtype cells could be monitored. By extracting the adaptation times for chemotactic stimuli, the adaptation rate constants for individual cells can be inferred. Subsequently, the same cells can be used to measure long time traces of motor switching and noise spectra can be calculated. While we included a number of noise sources in our model, others have not been included. For instance the binding and unbinding of the adaptational enzymes to the Ctermini of the receptors [^{56}] may introduce lowfrequency fluctuations in the motor noise power spectrum. However, as shown for enzyme CheR in Figure 9A, this contribution is significantly smaller than the methylation and demethylation noise. In the same plot we also show the CheR geneexpression noise, which only contributes at very low frequencies (observation times larger than 10^{3 }s). Furthermore, the ligand noise may be underestimated in our model since fluctuations of the ligand concentration at individual receptor complexes were treated as uncorrelated. If the fluctuations are treated as strongly correlated (on the scale of the polar chemoreceptor cluster), the ligand noise becomes substantially larger as complexes coherently transmit the ligand noise (see Figure 9B). Taken together, we believe that we included the main noise sources in our model although due to subtleties, the magnitude of noise contributions are much harder to pinpoint exactly. Our fullpathway model allows us to connect to a large variety of data and literature. For instance, we also investigated how the response to concentration signals is related to fluctuations in the chemotaxis pathway (Figure 8), similar to ideas presented by Park et al. [^{50}]. We do not find one unifying fluctuationresponse theorem, but a nonlinear trend for the relationship between adaptation times and variance in CheYP. For small adaptation times we find an approximate linear relationship in line with experiments [^{50}].
Although chemotaxis is one of many capabilities a cell has and may not be optimised in isolation without the rest of the cell, we speculate the cell aims to maximise the SNR for most efficient signalling and chemotaxis. We found that the SNR is maximised at particular receptor complex sizes similar to values of receptor cooperativity extracted from FRET doseresponse curves [^{46}]. In line with the data, the "optimal" complex size increases with external ligand concentration, and hence with receptor methylation level. While our complex sizes appear overestimated, noise from ligand molecules rebinding to the same receptor complex [^{57}] has not been considered here. This may well increase the noise level from external sources and hence decrease the predicted optimal receptor complex size. Using our model, we also analysed the effect of varying the methylation and demethylation rate constants. We found that a large SNR is obtained for small methylation and large demethylation rate constant, corresponding to our fitted model parameters from FRET doseresponse curves [^{47}].
To describe the noise, the Langevin approximation is expected to work for the phosphorylation and dephosphorylation of the abundant protein CheY. However, its applicability is less clear for receptor signalling due to both extrinsic ligand noise and intrinsic noise from receptor methylation. Furthermore, the switching of the binary motor may constitute relatively large noise. As shown in Additional File 1 the Master equation and Langevin approximation yield the same results for receptor signalling [^{58}]. As for the motor, the Langevin equation is exact [^{48}]. However, we explicitly tested that the statistical properties of the time series obtained for twostate switching and Langevin equation are the same. For constant rates, as well as for noisy rates due to fluctuations in CheYP concentration, the power spectra obtained for the two processes are the same (Additional File 1). In Additional File 1 we further show that noise terms are indeed sufficiently small that linearisation of the pathway equations is justified.
We presented a comprehensive analysis of signal and noise propagation in the chemotaxis pathway of E. coli. Our fullpathway model allowed us to successfully explain a large variety of experimental data from the literature. However, our work goes beyond E. coli chemotaxis as similar methods can be applied to many other pathways.
The E. coli chemotaxis pathway is a member of the large class of twocomponent systems, containing hundreds of closely related pathways involved in stress response, virulence and intercell communication [^{20},^{59}]. In these pathways, activation of a sensor histidine kinase results in its autophosphorylation, and subsequently in phosphorylation of a response regulator, which typically binds to DNA and regulates gene expression. The final output, i.e. activation of gene expression, is again binary and hence similar to the bacterial chemotaxis pathway (Figure 10). There are particularly wellstudied examples of twocomponent systems: (i) the VanS (kinase)/VanR (response regulator) system conferring vancomycin resistance in Gram positive bacteria [^{60}], (ii) quorum sensing in Vibrio harveyi, where the three kinases LuxN, LuxQ and CqsS respond to different autoinducers and first phosphorylate the phosphotransferase LuxO (which has no equivalent in the chemotaxis pathway), which then phosphorylates the response regulator LuxU [^{61}], and (iii) the phosphorelay controlling sporulation in Bacillus subtilis. The relay contains at least four kinases KinAKinB and the phosphotransferase Spo0F, which phosphorylates the response regulator Spo0A [^{62}].
Our work may also be applicable to other bacterial chemotaxis pathways, which are often considerably more complex than E. coli's pathway and still need to be understood in detail. For instance, the photosynthetic bacterium Rhodobacter sphaeroides has several homologues of each of the chemotaxis proteins in E. coli [^{63}]. Interestingly, Rhodobacter has two chemotaxis receptor clusters, one polar cluster similar to E. coli and one cytoplasmic cluster, which is thought to sense the metabolic state of the cell. Both clusters need to be present for chemotaxis [^{63}]. The soil bacterium Bacillus subtilis has three adaptation systems [^{64}]: one based on methylation and demethylation of receptors similar to E. coli and two independent of receptor methylation, the CheC/CheD system and the CheV system. Furthermore, in B. subtilis, sensory adaptation is not determined by the level of receptor methylation but the location of methylation sites on the receptors. Hence, our analysis presented here may also help elucidate the design of these pathways and clarify the computational problems these bacteria try to solve.
The MWC model describes signalling by receptor complexes, which can switch between their on and an off states. The average activity of a complex is given by
(11)
A=N1+eF(c,M), 
ranging from zero to N. The freeenergy difference F(c, M) between the on and off state is
(12)
F(c,M)=N12M+Nνa ln1+c∕Kaoff1+c∕Kaoff+νs ln1+c∕Ksoff1+c∕Ksoff, 
which is a function of the concentration c present at the receptor complex site and the methylation level M of the receptor complex. The methylation level of a complex is the sum of methylation levels of all receptors in a complex. Here, we consider two receptor types, Tar (indicated by index a) with fraction υ_{a }of receptors in the complex, and Tsr (indicated by index s) with fraction υ_{s }of receptors. Receptors are sensitive to attractant MeAsp with dissociation constants K^{on }and K^{off }in the on and off state, respectively. We use the following parameters for the MWC model for receptor complexes: Kaoff=0.02mM, Kaon=0.5mM, Ksoff=100mM and Kson=106mM[^{44},^{47}].
Equations 14 represent a simplified model of the chemotaxis pathway to illustrate signalling and noise transmission. In the following, our model used for the full signalling pathway is summarised. Equations 1 and 2 equally apply to the full pathway model. However, instead of Eq. 3, the dynamics of the methylation level of complex j in the full pathway model is described by
(13)
dMjdt=γR(NAj)γBAjBp2+ηMj(t)=γR(NAj)γBVcell2AjNBp2+ηMj(t) 
Here, we explicitly include the number of CheBP (B_{p}) molecules NBp in the demethylation term, with V_{cell }the cell volume.
In addition, we take into account phosphorylation and dephosphorylation of CheA (A_{p}), CheY (Y_{p}) and CheB (B_{p}), which are described by the following equations:
(14)
dNApdt=AckANCN(NA,totNAp)+kyVcell(NY,totNYp)NAp+kbVcell(NB,totNBp)NAp++ηA,p(t)+ηA,Yp(t)+ηA,Bp(t) 
(15)
dNYpdt=kyVcell(NY,totNYp)NApkyNYp+ηA,Yp(t)+ηYp(t) 
(16)
dNBpdt=kbVcell(NB,totNBp)NApkbNBp+ηA,Bp(t)+ηBp(t) 
with N_{i }the number of molecules of species i in a cell volume. Note, that we neglected the binding of CheYP to its phosphatase CheZ for simplicity, and describe dephosphorylation of CheYP by the effective dephosphorylation rate kyNYp. As shown in Additional File 1 this simplification has no qualitative effect on the response function. The term η_{A,p }describes the noise associated with CheA autophosphorylation. The terms ηA,Bp(t) and ηA,Yp(t) represent the noise generated in phosphorylation of CheB and CheY by CheA, respectively. The terms ηBp(t) and ηYp(t) describe the noise associated with dephosphorylation.
Note that some noise terms appear in two equations. This is due to the fact that we assign noise terms to a specific process, e.g. phosphorylation of CheY by CheA. Hence, the corresponding noise term ηA,Yp(t) appears in the dynamics of CheAP and of CheYP. As a positive fluctuation in the dynamics of CheAP due to phosphorylation of CheY corresponds to a negative fluctuation in the dynamics of CheYP, these noise terms appear with opposite signs in the two equations. The noise intensities and parameter values of the model are summarised below. Finally, instead of Eq. 4 the dynamics of the motor is now described in terms of
(17)
dXdt=k+(NYp)(1X)k(NYp)X+ηX(t), 
where we use experimentally derived switching rates k_{+ }and k_{}as a function of CheYP concentration (cf. Figure 3C and below).
Turner et al. presented a model for motor switching to explain the observed motor switching rates [^{39}]. The model for motor switching is an MWC model, where 26 subunits of the motor assume one of two states corresponding to CW and CCW rotation. While these subunits bind the molecule CheYP independently of each other, the switching of states occurs cooperatively. The authors derive the overall rates of switching (averaging over all possible CheYP occupancy states) as
(18)
k+(Yp)=k+(0)⋅1+μYpKCCW1+YpKCCWmcoop 
(19)
k(Yp)=k(0)⋅1+μYpKCCW1+YpKCWmcoop, 
where m_{coop }= 26 is the number of motor subunits, k_{B}T ln(μ) is the freeenergy difference of switching per molecule of CheYP, and K_{CCW }and K_{CW }are the dissociation constants for binding CheYP in the CCW and CW state, respectively.
Motor switching rates k_{+ }and k_{}have been derived experimentally as a function of the concentration of a signalling mutant CheY**, which is constitutively active [^{39}], as shown in Figure 3C. The cellular concentration of Y** varied as a result of induced expression in cells in which wildtype CheY and other cytoplasmic chemotaxis proteins had been deleted. To obtain the switching rates in terms of CheYP, rather than signalling mutant CheY**, we rescaled the dissociation constants of CheY binding to the motor such that the switching rates are equal, i.e. CW bias about 1/2, at CheYP concentration 3.2 μM [^{65}]. We fitted the above model to the experimental data in Figure 3C and the CW bias at 33° [^{39}] and used the rate constants k_{+ }and k_{}in our full pathway model (Eq. 17).
The mean, or deterministic response of a system to an arbitrary, sufficiently small signal is described in terms of a characteristic function of the system, the linear response function χ_{R}. The timedependent response ΔR(t), i.e. the deviation from the adapted state due to a small input signal Δc(t), is linear and determined by
(20)
ΔR(t)= ∫ ∞tχR(tτ)Δc(τ)dτ. 
Hence, the time course of the response is determined by the convolution of the linear response function and the input signal. The linear response function describes the dynamics of the pathway and the convolution with the input signal represents the fact that the current state of the system is determined by the history of the input signal [^{66}]. Importantly, for a nonlinear system to respond linearly to a stimulus, the input signal has to be small. For the receptor complex activity, a signal is small if the change in ligand concentration is small compared to the apparent dissociation constants for the on and off states. These dissociation constants depend on the background concentration as the system adapts to the background concentration c_{0}, resulting in a particular freeenergy difference F*. Thus, the concentration change has to be small compared to c0+Kaon and c0+Kaoff. In this case, the change in freeenergy difference due to the concentration change is much smaller than k_{B}T. In Figure 4, the background concentration c_{0 }= 0 and the peak concentration of the impulse is 10^{3 }mM, compared to Kaon=0.5mM and Kaoff=0.02mM. Hence, the input signal is indeed small.
The Fourier transform of Eq. 20 reads more simply
(21)
ΔR^(ω)=χ^R(ω)Δĉ(ω). 
χ^R(ω) is also called the frequencydependent gain [^{12}]. The modulus χ^R(ω) describes what frequencies of the input signal are transmitted well, and which ones are attenuated.
Typically, finite activation rates of the system limit the response to rapidly changing input signals, i.e. highfrequency signals. In this case, the Fourier transformed linear response functions falls off at highfrequencies, and the system is called a lowpass filter. If lowfrequency components of the input signal are filtered out rather than high frequencies, the system is called a highpass filter. The chemotaxis pathway is a bandpass filter (see Figure 2), filtering out low and highfrequency components.
To obtain a succinct measure for the signalling response due to an input concentration change Δc(t), we define ΔR^{2 }the integral of the response over frequency
(22)
ΔR2=12π∫ ∞∞dωχ^R(ω)Δĉ(ω)2. 
Fluctuations δR(t) around the mean value 〈R(t)〉 can be characterised by their correlations at different time points. The autocorrelation function K is defined as
(23)
K(τ)=⟨δR(t)δR(t+τ)⟩, 
i.e. the average value of the product of fluctuations at two time points. It only depends on the interval between time points if the dynamics of R(t) is stationary, i.e. the mean value 〈R(t)〉 and variance 〈δR^{2}(t)〉 do not vary with time t. Averaging over different measurements of R(t) is indicated by angular brackets. Typically, correlations decay with the interval length τ between time points. Here, we use the power spectrum to characterise fluctuations. According to the WienerKinchin theorem the power spectrum is the Fourier transform of the autocorrelation function [^{67}],
(24)
SR(ω)= ∫ ∞∞K(t)eiωtdt. 
For exponentially decaying correlations as in Figure 2, the power spectrum is Lorentzian, i.e. has the frequency dependency
(25)
SR(ω)∝1ω2+ωα2. 
The power spectrum can be calculated from the absolute square of the Fourier transform of time series δR(t) measured or simulated over a time interval T,
(26)
SR(ω)= limT→∞⟨δR^(ω)δR^*(ω)⟩T, 
where the Fourier transform is defined on the finite measurement interval T and the average 〈·〉 is over multiple time series.
The variance of a stationary process can be calculated as the integral of the power spectrum over frequency,
(27)
⟨δR2⟩=12π∫ ∞∞dωSR(ω). 
Here, we use the Langevin approximation to describe fluctuations (noise) around average molecule numbers. For instance, consider the following general stochastic differential equation:
(28)
dRdt=r1r2+η(t). 
The rates r_{1 }and r_{2 }typically depend on the concentrations of proteins in the signalling network. The noise term η(t) is composed of two terms η_{1}(t) and η_{2}(t), which are associated with the rates r_{1 }and r_{2}, respectively. We assume η_{1 }and η_{2 }to be independent, i.e. 〈η_{1}(t) η_{2}(t')〉 = 0. In general, this is justified as different reactions are catalysed by different proteins. Using 〈η_{j}(t)η_{j}(t')〉 = Q_{j}δ(t t'), the noise intensities can be calculated if we make the assumption that fluctuations are due to socalled birth and death processes, i.e. creation and destruction of the molecules with average rates r1* and r2*. Then the associated noise intensities are Q1=r1* and Q2=r2*[^{8}]. The intensity of the total noise η(t) is the sum Q = Q_{1 }+ Q_{2 }due to the independence of the two noises. As forward and backward rate are equal at steady state, Q is twice the reaction rate in one direction at steadystate.
The switching noise η_{A}(t) in Eq. 1 is due to the switching of each receptor complex between on and off states. We assume the switching to be a fast process, which can be described by the following dynamics of the receptorcomplex onbias a:
(29)
dadt=k1(Na)k2a+ηa(t). 
The noise term η_{a}(t) is a Gaussian white noise with zero mean and noise intensity Q_{a }= 2k_{2}A*, where we used that the receptor complex activity A = 〈a〉 which is equal to the (quasi) steadystate activity of a, and A = A* when adapted. The power spectrum of a due to switching between on and off states is
(30)
Sa(ω)=Qaω2+(k1+k2)2, 
where k_{1 }+ k_{2 }is the characteristic frequency of switching. Hence, the highfrequency component of fluctuations δa(t) is reduced due to averaging by the finite rates of switching. Hence, the power spectrum of activity fluctuations η_{A}(t), which are due to fluctuations in the rate of receptor complex switching da/dt, is
(31)
SηA(ω)=ω2Sa(ω). 
This result is due to the Fourier transform of the rate of receptor complex switching da/dt yielding iωa^. Hence, the power spectra of a and da/dt differ by a factor ω^{2 }as the power spectrum is proportional to the magnitude squared of the Fourier transform according to Eq. 26.
The number of ligand molecules in the vicinity of a receptor complex fluctuates due to binding/unbinding, and potential rebinding of previously bound ligand molecules at this complex, as well as diffusion [^{57},^{68}]. Here, we use a simplified description of diffusion to calculate the spectrum of noise in the ligand dynamics η_{c}(t) in Eq. 2. Consider a volume whose dimensions are given by the diameter of a receptor complex s=NsR, where s_{R }= 1 nm is the size of a receptor dimer [^{69}]. The change of ligandmolecule number L in this volume is determined by the exchange rate k_{D }≈ D/(2s^{2}) due to diffusion [^{70}]:
(32)
dLdt=kD(c0s3L)+ηL(t) 
where k_{D}L is the rate of molecules moving out of the volume by diffusion, and k_{D }times the mean concentration c_{0 }in solution serves as a proxy of the rate of ligand molecules moving into the volume. The noise term η_{L}(t) is assumed to be Gaussian and white, with zero mean and noise intensity Q_{L }= Dsc_{0}.
The power spectra of the number L and concentration c of molecules at receptor complex j are respectively
(33)
SL(ω)=Dsc0ω2+kD2, 
(34)
Sc(ω)=SL(ω)s6, 
where s^{6 }is the squared volume given by the dimension of the receptor complex. The zerofrequency limit of the power spectrum of the ligand concentration S_{c}(0) = c_{0}/(Ds), which corresponds to calculations by Berg and Purcell [^{71}] and Bialek and Setayeshgar [^{57}] for the uncertainty in sensing ligand concentration.
The noise η_{c}(t) in Eq. 2 is related to rate of change of the ligand concentration, similar to the considerations of the switching noise above. Hence, the power spectrum of the ligand fluctuations η_{c}(t) is
(35)
Sηc(ω)=ω2Sc(ω). 
The size of fluctuations in the rate of methylation of a receptor complex j in Eq. 3 is estimated from the average rates of methylation and demethylation at the adapted state, respectively. The noise η_{M }(t) is assumed to be Gaussian and white, with zero mean, noise intensity Q_{M }= 2γ_{R}(N A*) and power spectrum
(36)
SηM(ω)=QM. 
The noise in motor switching rate in Eq. 4 is assumed to be a Gaussian white noise term with zero mean, noise intensity QX=2k+(Ac*)(1PCW*) and power spectrum
(37)
SηX(ω)=QX. 
After linearising around the steady state and inserting the Fourier transforms we obtain for the simplified model
(38)
iωΔA^=iω∂A∂MΔM^iω∂A∂cΔĉ 
(39)
iωΔM^=ω1ΔA^ 
(40)
iωΔP^CW=ω2ΔA^cωXΔP^CW, 
where
(41)
ω1=γR+3γBA*2=γR(32Ar*)∕Ar* 
with A*=N⋅Ar*≈N∕3[^{72}] the adapted activity of a receptor complex, Ar* denoting the adapted activity of individual receptors. In the second equality we have used that at the adapted state γR(NA*)=γBA*3. The parameter ω2=(1PCW*)∂k+∂AcPCW*∂k∂Ac is the derivative of the motor switching rates with respect to activity, and ω_{X }= k_{+}* + k_{}* is a characteristic frequency due to motor switching at steady state. ΔA is the response of every receptor signalling complex, and ΔA_{c }= N_{C}ΔA is the activity response of all receptor complexes in a cell. Solving for ΔA^c and ΔP^CW, and division by the stimulus Δĉ yields the response functions in Eq. 5 and 6.
To calculate spectra, we linearise the deterministic parts of Eq. 14 similar to the calculation of the response functions, and formally Fourier transform the equations. We obtain
(42)
iωδA^c=iω∂A∂M∑jδM^j+∂A∂c∑jη^cj++ ∑jη^Aj 
(43)
iωδM^j=ω1δA^j+η^Mj 
(44)
iωδP^CW=ω2δA^cωXδP^CW+η^X. 
We solve for the Fourier transformed activity fluctuations δA^c and obtain
(45)
δA^c=∂A∂M∑jη^Mj+∂A∂c∑jηcj+ ∑jη^AjωMiω, 
which yields the power spectrum in Eq. 8. The parameter ω_{M }= ω_{1}∂A/∂M, and we used Eq. 31 and 35. From Eq. 44 we obtain for the Fourier transformed fluctuations in the motor bias δP^CW
(46)
δP^CW=ω2δA^c+η^XωXiω, 
and their power spectrum is given by Eq. 9.
The integrated response of the receptor activity to a step stimulus is
(47)
ΔAc2= ∫ ∞∞dωχ^Ac(ω)Δĉ(ω)2=πNC2∂A∂c2(αc)2ωM, 
where we inserted Eq. 5. Hence, the activity response scales as ΔAc2∝Ntot∕N2N22∕N∝N, where we used that N_{C }= N_{tot}/N with N_{tot }the total number of receptors in a cell.
The variance of the receptor activity is given by the integral of the power spectrum of activity fluctuations Eq. 8
(48)
⟨δAc2⟩=NC2π∫ τ1τ1dωω2Sa(ω)+∂A∂c2Sc(ω)ω2+ωM2+NC2π∫ τ1τ1dω∂A∂M2QMω2+ωM2, 
where we consider the frequency range relevant for motor switching indicated by τ^{1 }≈ 0.1 ... 1 Hz. The contribution from receptor switching is
(49)
⟨δAc2⟩a=NC2π∫ τ1τ1dωω2Sa(ω)ω2+ωM2≈2k2Ar*Ntotπτ(k1+k2)2 
where we used Q_{a }and inserted Eq. 30 for the power spectrum of receptor switching noise and used that it is almost constant and equal to its zerofrequency value over the integration range. Furthermore, the factor ω2∕(ω2+ωM2)≈1 and Ar*=A*∕N is the adapted activity of an individual receptor. Hence, according to this simple calculation the contribution to the variance from receptor switching is roughly constant with receptor complex size.
The contribution from ligand diffusion is
(50)
⟨δAc2⟩c=NC2π∫ τ1τ1dωω2∂A∂c2Sc(ω)ω2+ωM2≈NC∂A∂c2⟨δc2⟩, 
where 〈δc^{2}〉 = c_{0}/(Dsτ) is the variance of the ligand concentration measured during the time interval τ. We used Eq. 34 and the same argument as for the switching noise to calculate the integral. Hence, the contribution to the variance from the ligand diffusion grows as ⟨δAc2⟩c∝N5∕2 as a result of incoherent addition of noise from different receptor complexes, the sensitivity ∂A/∂c increasing as N^{2}, and the size dependence of receptor complexes s ∝ N^{1/2}.
The contribution to the variance from receptor methylation is
(51)
⟨δAc2⟩M=NC2π∂A∂M2 ∫ τ1τ1dωQMω2+ω1∂A∂M2≈2NtotγR(1Ar*)Ar*ω1∂A∂M 
where we defined ω1=γR+3γBN2(Ar*)2, inserted QM=2γRN(1Ar*) and ω_{M }= ω_{1}(∂A/∂M). Hence, ⟨δAc2⟩M grows approximately linearly with receptor complex size.
The SNR grows linearly with N for small complex sizes, and decreases as N^{2 }for larger complex sizes, resulting in an optimal medium receptor complex size, in qualitative agreement with Figure 6.
The integrated signal response of the receptor activity Eq. 47
(52)
ΔAc2=πNC2∂A∂c2(αc)2ω1∂A∂M, 
where the numerator expresses the initial response of receptors of concentration changes and the denominator the filtering by adaptation. The sensitivity ∂A∕∂c=NAr*(1Ar*)h(c), where h(c) = ∂F/∂c, ω1=γR+3γBN2(Ar*)2, and ∂A∕∂M=NAr*(1Ar*)∕2. The adapted activity can be obtained analytically for our simplified model from the steady state of the methylation dynamics Eq. 3,
(53)
Ar*=12β+β24+β3273β312β+β24+β3273, 
and is only a function of the ratio β = γ_{R}/γ_{B}. Expanding the adapted activity around Ar*=0 (for γ_{R }→ 0) yields Ar*∝γR1∕3, and around Ar*=1 (for γ_{B }→ 0) yields Ar*∝γB. Similarly, ω1∝γR2∕3(const.+γB4∕3). Hence, ∂A∕∂c∝γR1∕3(γB) and ω1∂A∕∂M∝γR4∕3(γB).
The initial response to concentration changes decreases slower than adaptation times, resulting in an increased signal response for vanishing γ_{R}. For vanishing γ_{B}, the initial response to concentration changes decreases faster than adaptation speed, hence yielding a vanishing signal response. The overall dependence of the integrated signal response is ΔAc2∝γR1∕3(γB) for γ_{R }→ 0(γ_{B }→ 0). For the contributions to the variance of the receptor activity from receptor switching, ligand diffusion and receptor methylation dynamics we obtain ⟨δAc2⟩a∝γR1∕3(γB), ⟨δAc2⟩c∝γR2∕3(γB2)and ⟨δAc2⟩M∝γR(γB2∕3), respectively. Hence, according to our simplified model the SNR of the receptor activity goes as SNR∝γR2∕3(γB4∕3), in qualitative agreement with Figure 7.
Only the full pathway model without noise was solved numerically (cf. Figure 4). We used an Euler method to integrate the differential equations 1317 without noise. As concentration input signal we used a triangular impulse of duration 0.02 s, c_{0 }= 0 and peak concentration 10^{3 }mM. The integration time step was Δt = 10^{3}, maximum integration time was 100 s. In each integration step, first the receptor complex activity was calculated according to Eq. 11 based on the current ligand concentration and receptor methylation level from the previous iteration. All other variables x were updated in each time step i according to x_{i+1 }= x_{i }+ Δt · rate, where rate is given by the r.h.s. of the differential equations for the respective variable.
Rate constants and total cell concentrations of proteins for the full pathway model are given in Table 1. The noise terms ηAj, ηcj, ηMj and η_{X }are the same as in Eq. 24 and 2729 and their power spectra are given there. The noise associated with phosphorylation and dephosphorylation ηAp, ηA,Bp, ηA,Yp, ηBp and ηYp are assumed to be Gaussian white noise terms with zero mean and autocorrelations 〈η_{i}(t)η_{i}(t')〉 = Q_{i}δ(t t') with noise intensities Q_{i }given in Table 2. Fitting parameters of the Fourier transformed response function Figure 3 are listed in Table 3. Parameters for Figure 5 are listed in Tables 4 and 5, and those for Figure 8 are listed in Table 6.
DC and RGE conceived and designed the study, performed analytical calculations, analysed the data, and wrote the paper. DC performed computer simulations. All authors read and approved the final manuscript.
Supplementary information. Supplementary information in .pdf format, including the extended model for the chemotaxis pathway, additional mathematical derivations and validity checks.
Click here for additional data file (175205095151S1.PDF)
We thank Richard Berry, Martin Buck, Tom Duke and William Ryu for helpful discussions, and Suhail Islam for computational support. RGE was supported by Biotechnological and Biological Sciences Research Council grant BB/G000131/1 and the Centre for Integrative Systems Biology and Bioinformatics (CISBIO) at Imperial College.
References
Sykulev Y,Joo M,Vturina I,Tsomides TJ,Eisen HN,Evidence that a Single Peptide MHC Complex on a Target Cell Can Elicit a Cytolytic T Cell ResponseImmunityYear: 1996456557110.1016/S10747613(00)8048358673703  
Mortimer D,Pujic Z,Vaughan T,Thompson AW,Feldner J,Vetter I,Goodhill GJ,Axon guidance by growthrate modulationProc Natl Acad Sci USAYear: 20101075202520710.1073/pnas.090925410720194766  
Elowitz MB,Levine AJ,Siggia ED,Swain PS,Stochastic Gene Expression in a Single CellScienceYear: 20022971183118610.1126/science.107091912183631  
Pedraza JM,van Oudenaarden A,Noise propagation in gene networksScienceYear: 20053071965196910.1126/science.110909015790857  
Paulsson J,Models of stochastic gene expressionPhys Life RevYear: 2005215717510.1016/j.plrev.2005.03.003  
Acar M,Pando BF,Arnold FH,Elowitz MB,van Oudenaarden A,A General Mechanism for NetworkDosage Compensation in Gene CircuitsScienceYear: 20103291656166010.1126/science.119054420929850  
Eldar A,Elowitz MB,Functional roles for noise in genetic circuitsNatureYear: 201046716717310.1038/nature0932620829787  
Thattai M,van Oudenaarden A,Attenuation of noise in ultrasensitive signaling cascadesBiophys JYear: 2002822943295010.1016/S00063495(02)75635X12023217  
Shibata T,Fujimoto K,Noisy signal amplification in ultrasensitive signal transductionProc Natl Acad Sci USAYear: 200510233133610.1073/pnas.040335010215625116  
ColmanLerner A,Gordon A,Serra E,Chin T,Resnekov O,Endy D,Pesce CG,Brent R,Regulated celltocell variation in a cellfate decision systemNatureYear: 200543769970610.1038/nature0399816170311  
Taylor RJ,Falconnet D,Niemistö A,Ramsey SA,Prinz S,Shmulevich I,Galitski T,Hansen CL,Dynamic analysis of MAPK signaling using a highthroughput microfluidic singlecell imaging platformProc Natl Acad Sci USAYear: 20091063758376310.1073/pnas.081341610619223588  
Detwiler PB,Ramanathan S,Sengupta A,Shraiman BI,Engineering aspects of enzymatic signal transduction: photoreceptors in the retinaBiophys JYear: 2000792801281710.1016/S00063495(00)76519211106590  
Hornung G,Barkai N,Noise Propagation and Signaling Sensitivity in Biological Networks: A Role for Positive FeedbackPLoS Comput BiolYear: 20084e810.1371/journal.pcbi.004000818179281  
Lestas I,Vinnicombe G,Paulsson J,Fundamental limits on the suppression of molecular fluctuationsNatureYear: 201046717417810.1038/nature0933320829788  
Ueda M,Shibata T,Stochastic Signal Processing and Transduction in Chemotactic Response of Eukaryotic CellsBiophys JYear: 200793112010.1529/biophysj.106.10026317416630  
Emonet T,Cluzel P,Relationship between cellular response and behavioral variability in bacterial chemotaxisProc Natl Acad Sci USAYear: 20081053304330910.1073/pnas.070546310518299569  
Mehta P,Goyal S,Wingreen NS,A quantitative comparison of sRNAbased and proteinbased gene regulationMol Syst BiolYear: 2008422118854820  
Tostevin F,ten Wolde PR,Mutual Information between Input and Output Trajectories of Biochemical NetworksPhys Rev LettYear: 200910221810119519137  
Mora T,Wingreen NS,Limits of Sensing Temporal Concentration Changes by Single CellsPhys Rev LettYear: 201010424810120867338  
Laub MT,Goulian M,Specificity in TwoComponent Signal Transduction PathwaysAnn Rev GenetYear: 20074112114510.1146/annurev.genet.41.042007.17054818076326  
Berg HC,Motile behavior of bacteriaPhys TodayYear: 2000532429  
Falke JJ,Hazelbauer GL,Transmembrane signaling in bacterial chemoreceptorsTrends Biochem SciYear: 20012625726510.1016/S09680004(00)01770911295559  
Sourjik V,Receptor clustering and signal processing in E. coli chemotaxisTrends MicrobiolYear: 20041256957610.1016/j.tim.2004.10.00315539117  
Wadhams GH,Armitage JP,Making sense of it all: bacterial chemotaxisNat Rev Mol Cell BiolYear: 200451024103710.1038/nrm152415573139  
Baker MD,Wolanin PM,Stock JB,Systems biology of bacterial chemotaxisCurr Opin MicrobiolYear: 2006918719210.1016/j.mib.2006.02.00716529985  
Block SM,Segall JE,Berg HC,Impulse responses in bacterial chemotaxisCellYear: 19823121522610.1016/00928674(82)9042146760985  
Segall JE,Block SM,Berg HC,Temporal comparisons in bacterial chemotaxisProc Natl Acad Sci USAYear: 1986838987899110.1073/pnas.83.23.89873024160  
Shimizu TS,Tu Y,Berg HC,A modular gradientsensing network for chemotaxis in Escherichia coli revealed by responses to timevarying stimuliMol Syst BiolYear: 2010638220571531  
Block SM,Segall JE,Berg HC,Adaptation kinetics in bacterial chemotaxisJ BacteriolYear: 19831543123236339475  
Tu Y,Shimizu TS,Berg HC,Modeling the chemotactic response of Escherichia coli to timevarying stimuliProc Natl Acad Sci USAYear: 2008105148551486010.1073/pnas.080756910518812513  
Korobkova E,Emonet T,Vilar J,Shimizu T,Cluzel P,From molecular noise to behavioral variability in a single bacteriumNatureYear: 200442857457810.1038/nature0240415058306  
MortonFirth CJ,Bray D,Predicting temporal fluctuations in an intracellular signalling pathwayJ Theor BiolYear: 199819211712810.1006/jtbi.1997.06519628844  
Xing J,Bai F,Berry R,Oster G,Torquespeed relationship of the bacterial flagellar motorProc Natl Acad Sci USAYear: 20061031260126510.1073/pnas.050795910316432218  
Meacci G,Tu Y,Dynamics of the bacterial flagellar motor with multiple statorsProc Natl Acad Sci USAYear: 20091063746375110.1073/pnas.080992910619234112  
Mora T,Yu H,Sowa Y,Wingreen NS,Steps in the bacterial flagellar motorPLoS Comput BiolYear: 20095e100054010.1371/journal.pcbi.100054019851449  
Mora T,Yu H,Wingreen NS,Modeling torque versus speed, shot noise, and rotational diffusion of the bacterial flagellar motorPhys Rev LettYear: 200910324810220366231  
van Albada SB,TănaseNicola S,ten Wolde PR,The switching dynamics of the bacterial flagellar motorMol Syst BiolYear: 2009531619888211  
Scharf BE,Fahrner KA,Turner L,Berg HC,Control of direction of flagellar rotation in bacterial chemotaxisProc Natl Acad Sci USAYear: 19989520120610.1073/pnas.95.1.2019419353  
Turner L,Samuel AD,Stern AS,Berg HC,Temperature dependence of switching of the bacterial flagellar motor by the protein CheY(13DK106YW)Biophys JYear: 19997759760310.1016/S00063495(99)76916X10388784  
Tu Y,Grinstein G,How White Noise Generates PowerLaw Switching in Bacterial Flagellar MotorsPhys Rev LettYear: 20059420810116090291  
van Kampen NG,Stochastic processes in physics and chemistry. North HollandYear: 2007  
Sourjik V,Berg HC,Functional interactions between receptors in bacterial chemotaxisNatureYear: 200442843744110.1038/nature0240615042093  
Mello BA,Tu Y,An allosteric model for heterogeneous receptor complexes: understanding bacterial chemotaxis responses to multiple stimuliProc Natl Acad Sci USAYear: 2005102173541735910.1073/pnas.050696110216293695  
Keymer JE,Endres RG,Skoge M,Meir Y,Wingreen NS,Chemosensing in Escherichia coli: two regimes of twostate receptorsProc Natl Acad Sci USAYear: 20061031786179110.1073/pnas.050743810316446460  
Endres RG,Wingreen NS,Precise adaptation in bacterial chemotaxis through "assistance neighborhoods"Proc Natl Acad Sci USAYear: 2006103130401304410.1073/pnas.060310110316924119  
Endres RG,Oleksiuk O,Hansen CH,Meir Y,Sourjik V,Wingreen NS,Variable sizes of Escherichia coli chemoreceptor signaling teamsMol Syst BiolYear: 2008421118682701  
Clausznitzer D,Oleksiuk O,Sourjik V,Endres RG,Chemotactic response and adaptation dynamics in Escherichia coliPLoS Comput BiolYear: 20106e100078410.1371/journal.pcbi.100078420502674  
Stock G,Ghosh K,Dill KA,Maximum Caliber: A variational approach applied to twostate dynamicsJ Chem PhysYear: 200812819410210.1063/1.291834518500851  
Li M,Hazelbauer GL,Cellular Stoichiometry of the Components of the Chemotaxis Signaling ComplexJ BacteriolYear: 20041863687369410.1128/JB.186.12.36873694.200415175281  
Park H,Pontius W,Guet CC,Marko JF,Emonet T,Cluzel P,Interdependence of behavioural variability and response to small stimuli in bacteriaNatureYear: 201046881982310.1038/nature0955121076396  
Korobkova E,Emonet T,Park H,Cluzel P,Hidden stochastic nature of a single bacterial motorPhys Rev LettersYear: 20069658105  
Blat Y,Gillespie B,Bren A,Dahlquist FW,Eisenbach M,Regulation of phosphatase activity in bacterial chemotaxisJ Mol BiolYear: 19982841191119910.1006/jmbi.1998.22249837737  
Blat Y,Eisenbach M,Oligomerization of the Phosphatase CheZ Upon Interaction with the Phosphorylated Form of CheYJ Biol ChemYear: 19962711226123110.1074/jbc.271.2.12268557654  
Blat Y,Eisenbach M,Mutants with Defective Phosphatase Activity Show No Phosphorylationdependent Oligomerization of CheZJ Biol ChemYear: 19962711232123610.1074/jbc.271.2.12328557655  
Eisenbach M,Eisenbach MBacterial chemotaxisChemotaxisYear: 2004Imperial College Press53215  
Hansen CH,Endres RG,Wingreen NS,Chemotaxis in Escherichia coli: a molecular model for robust precise adaptationPLoS Comput BiolYear: 20084e110.1371/journal.pcbi.004000118179279  
Bialek W,Setayeshgar S,Physical limits to biochemical signalingProc Natl Acad Sci USAYear: 2005102100401004510.1073/pnas.050432110216006514  
Aquino G,Clausznitzer D,Tollis S,Endres RG,Optimal receptorcluster size determined by intrinsic and extrinsic noisePhys Rev EYear: 201183021914  
Sourjik V,Armitage JP,Spatial organization in bacterial chemotaxisEMBO JYear: 2010292724273310.1038/emboj.2010.17820717142  
Hutchings MI,Hong HJ,Buttner MJ,The vancomycin resistance VanRS twocomponent signal transduction system of Streptomyces coelicolorMol MicrobiolYear: 2006599233510.1111/j.13652958.2005.04953.x16420361  
Henke JM,Bassler BL,Three parallel quorumsensing systems regulate gene expression in Vibrio harveyiJ BacteriolYear: 20041866902691410.1128/JB.186.20.69026914.200415466044  
Jiang M,Shao W,Perego M,Hoch JA,Multiple histidine kinases regulate entry into stationary phase and sporulation in Bacillus subtilisMol MicrobiolYear: 20003853534210.1046/j.13652958.2000.02148.x11069677  
Porter SL,Wadhams GH,Armitage JP,Rhodobacter sphaeroides: complexity in chemotactic signallingTrends in MicrobiologyYear: 20081625126010.1016/j.tim.2008.02.00618440816  
Glekas GD,Cates JR,Cohen TM,Rao C,Ordal GW,Sitespecific methylation in Bacillus subtilis chemotaxis: effect of covalent modifications to the chemotaxis receptor McpBMicrobiologyYear: 2011157566510.1099/mic.0.044685020864474  
Cluzel P,Surette M,Leibler S,An ultrasensitive bacterial motor revealed by monitoring signaling proteins in single cellsScienceYear: 20002871652165510.1126/science.287.5458.165210698740  
Kubo R,Statistical Mechanical Theory of Irreversible Processes IJournal of the Physical Society of JapanYear: 19571257058610.1143/JPSJ.12.570  
Stratonovich R,Topics in the Theory of Random Noise, Volume IYear: 1963New York: Gordon and Breach  
Endres RG,Wingreen NS,Accuracy of direct gradient sensing by cellsurface receptorsProg Biophys Mol BiolYear: 2009100333910.1016/j.pbiomolbio.2009.06.00219523978  
Hazelbauer GL,Bacterial chemoreceptorsCurr Opin Struct BiolYear: 1992250551010.1016/0959440X(92)90079M  
Berg HC,Random Walks in BiologyPrinceton University Press, Princeton 1993 chap. Diffusion: Macroscopic Theory  
Berg HC,Purcell EM,Physics of chemoreceptionBiophys JYear: 19772019321910.1016/S00063495(77)855446911982  
Sourjik V,Berg HC,Receptor sensitivity in bacterial chemotaxisProc Natl Acad Sci USAYear: 20029912312710.1073/pnas.01158999811742065  
Sourjik V,Berg HC,Binding of the Escherichia coli response regulator CheY to its target measured in vivo by fluorescence resonance energy transferProc Natl Acad Sci USAYear: 200299126691267410.1073/pnas.19246319912232047  
Shapovalov G,Lester HA,Gating Transitions in Bacterial Ion Channels Measured at 3 us ResolutionJ Gen PhysYear: 200412415116110.1085/jgp.200409087  
Wolanin PM,Baker MD,Francis NR,Thomas DR,DeRosier DJ,Stock JB,Selfassembly of receptor/signaling complexes in bacterial chemotaxisProc Natl Acad Sci USAYear: 2006103143131431810.1073/pnas.060635010316973743  
Stewart RC,Jahreis K,Parkinson JS,Rapid phosphotransfer to CheY from a CheA protein lacking the CheYbinding domainBiochemistryYear: 200039131571316510.1021/bi001100k11052668  
Bray D,Bourret RB,Computer analysis of the binding reactions leading to a transmembrane receptorlinked multiprotein complex involved in bacterial chemotaxisMol Biol CellYear: 19956136713808573792  
Stewart RC,Activating and inhibitory mutations in the regulatory domains of the methylesterase in bacterial chemotaxisJ Biol ChemYear: 1993268192119308420965 
Figures
[Figure ID: F1] 
Figure 1
Chemotaxis pathways of E. coli. The pathway consists of transmembrane chemoreceptors, which constantly undergo molecular transitions between their on (black) and off (white) conformations. Receptors signal to CheA, which subsequently autophosphorylates. The average receptor activity is modulated by ligand binding and unbinding, as well receptors methylation and demethylation. CheA phosphorylates the response regulator CheY, which diffuses through the cell and binds to the rotary motors. Upon CheYP binding, motors switch from their default state of CCW rotation (i.e. running mode) to CW rotation (i.e. tumbling mode). In addition, CheYP is dephosphorylated by its phosphatase CheZ. Receptor methylation is catalysed by CheR, which preferentially modifies inactive receptors. Receptor demethylation is catalysed by CheB, which is activated by phosphorylation, and modifies preferentially active receptors. 
[Figure ID: F2] 
Figure 2
Characterisation of signal and noise propagation. Signal response and fluctuations can be analysed in the time domain or frequency domain, the latter allowing for analytical treatment. Analysis of signal propagation: A small stimulus Δc(t) (Input) is applied, which results in a measurable response ΔR(t) (Output). The response ΔR(t) of the system to an impulse input represents the linear response function χ_{R}(t) (up to a constant factor). In the frequency domain, this stimulus is a constant. The Fourier transformed linear response function ΔR^∝χ^R(ω) can be analysed for its frequencyresolved transmission behaviour. Noise propagation: Fluctuations are characterised by their correlations over the time interval τ. The autocorrelation function K(τ) (Inset) typically decreases as a function of interval length. In the frequency domain, the noise power spectrum S_{R}(ω), which is the Fourier transform of the autocorrelation function, characterises the frequency components of the noise. 
[Figure ID: F3] 
Figure 3
Calibration of the model. (A and B) Fourier transformed linear response function. (A) Magnitude of the response function of the rotary motor measured by Block et al. [^{26}] (circles) and Segall et al. [^{27}] (plus symbols). The fit of our model is shown by the solid line. The dashed line represents a 3rdorder filter for comparison. (B) Magnitude (left) and phase (right) of the response function at the level of the response regulator CheY measured by Shimizu et al. [^{28}]. The symbols are measurements at 22°C (circles) and 32°C (plus symbols). Grey lines represent the fit of our model to the magnitude of the response. (C) Switching rates of the motor from CCW to CW rotation k_{+ }(squares) and from CW to CCW rotation k_{}(circles) as a function of the concentration of signallingactive unphosphorylated doublemutant Y**. The cellular concentration of Y** varied as a result of induced expression in cells in which wildtype CheY and other cytoplasmic chemotaxis proteins had been deleted. A fit using the model of Turner et al. [^{39}] is shown as well (solid and dashed lines; cf. Methods). 
[Figure ID: F4] 
Figure 4
Chemotaxis pathway response and noise spectra in the full pathway model. (Left and Middle) Response upon impulse stimulation with attractant MeAsp. (Left) Time courses for MeAsp concentration c, total activity of receptors A_{c}, CheAP and CheYP, and motor bias X obtained from simulation of the full pathway model using the Euler method for integration of the differential equations without noise. (Middle) Fourier transformed ligand signal, as well as response functions. Symbols correspond to the numerical simulation in the left panel, and solid lines to analytically calculated functions. (Right) Analytically calculated noise spectra of ligand and for the total activity of receptors A_{c}, CheAP and CheYP, and motor bias X (thick black lines). Also shown are the contributions to the spectrum from ligand binding (L; thin solid black lines), receptor switching (RS; thin red lines and symbols), receptor methylation and demethylation (M; green lines), as well as phosphorylation and dephosphorylation of CheA, CheY and CheB (P; blue lines) and motor switching (MS; dashed grey line). 
[Figure ID: F5] 
Figure 5
Effect of celltocell variation on power spectrum of the motor in the full pathway model. (A) Predictions about how different cell parameters affect the motor power spectrum, in particular its lowfrequency component. The black line (wildtype WT1) is the same as the total motor spectrum in Fig. 4, right. The motor spectra for increased motor switching rates (MS; red line), as well as reduced rates of receptor methylation and demethylation (M; green line) and the total number of receptors (blue line) are shown as well. Arrows indicate the features in the spectra that are affected. Dashed box is the area shown in panel (B). (B) Motor spectrum of cells with low motor bias (black line), as well as the component from steadystate motor switching only (grey line). (Inset) Measured power spectra for wildtype cells (WT2) with low motor bias (black) and mutant lacking the signalling pathway (grey). Axes are the same as in the main panel. Error bars indicate the measurement uncertainty. Spectra were traced from data presented by Korobkova et al. [^{31}]. Model parameters are listed in the Methods. (C) Motor spectrum of cells with varying CheR expression level. Shown are WT2 (black line; same as in panel B), as well as 2 (red), 4(green) and 10fold (blue) CheRexpression level compared to wildtype. (Inset) Measured power spectra of cells with corresponding CheRexpression level (see main panel). Axes are the same as in the main panel. Spectra were traced from data presented by Korobkova et al. [^{31}]. Model parameters are listed in the Methods. 
[Figure ID: F6] 
Figure 6
Integrated response, variance due to transmitted noise and SNR in the full pathway model when varying ambient concentration and receptor complex size. (A) Integrated response of the motor. (B) Variance of the motor including only contributions from receptor switching, ligand diffusion, methylation and phosphorylation. (C) Individual contributions to the variance of the motor. (D) SNR based on the signal response and variance shown in (A) and (B). (E) SNR as a function of receptor complex size at ambient concentration 0.02 (solid), 0.03 (dashed) and 0.05 mM (dotted line). Optimal complex size is indicated by arrows. 
[Figure ID: F7] 
Figure 7
Integrated response, variance due to transmitted noise and SNR in the full pathway model when varying receptor methylation and demethylation rate constants γ_{R }and γ_{B}, respectively. (A) Integrated response of the motor. (B) Illustration of the effects of vanishing γ_{R }and γ_{B }on adapted activity (indicated by dot and circle along doseresponse curve; left), as well as on time courses (right) for three cases, γ_{R }→ 0 (top right), γ_{R}/γ_{B }= const (middle right) and γ_{B }→ 0 (bottom right). For further explanation see text. (C) Variance of the motor including only from receptor switching, ligand diffusion, methylation and phosphorylation. (D) Individual contributions to the variance of the motor. (E) SNR based on the signal response and variance shown in (A) and (C). 
[Figure ID: F8] 
Figure 8
Fluctuationresponse relationships. (A) Variance of CheYP fluctuations (normalised by the squared adapted CheYP value) as a function of adapted CheYP value for varying adaptation rates γ_{R }(solid line) and γ_{B }(dashed line) and total CheY concentration (dotted line). The filled circle marks the point, where parameters for the three parameter variations coincide. (B) Adaptation time, calculated as (λ_{1}∂A/∂M)^{1}, plotted against the variance of CheYP fluctuations. Adapted CheYP concentration of 5μM is indicated by the open circles. Legend is the same as in (A). Our parameters are given in Additional File 1. (Inset) Variance of CheYP concentrations plotted against the adaptation time as extracted from measurements of motor rotation by Park et al. [^{50}]. The line is a linear fit through the data. 
[Figure ID: F9] 
Figure 9
Effects of CheR localisation and gene expression noise, as well as correlated ligand noise on noise power spectrum. (A) Total power spectrum of the motor excluding CheR localisation and gene expression noise (black line), and contributions to the motor power spectrum from receptor methylation and demethylation (green solid line), CheR localisation (dashed line) and gene expression noise (dotted line). (B) Total power spectrum of the motor excluding correlated ligand noise (solid black line), and contributions to the motor power spectrum from uncorrelated ligand noise at receptor complexes (solid grey line; original model) and correlated ligand noise across the whole receptor cluster (dashed grey line). For details, see Additional File 1. 
[Figure ID: F10] 
Figure 10
Comparison of E. coli chemotaxis and other twocomponent systems. In the chemotaxis pathway (left panel) and other twocomponent systems (right panel), a sensor kinase is activated by a cellexternal signal, upon which it autophosphorylates and passes on a phosphoryl group to its response regulator, which typically induces a transcriptional response. The time course of the final output, i.e. gene expression, can be directly mapped onto the binary output of the chemotaxis pathway. 
Tables
Parameters of the full pathway model, including references to literature.
Parameter  Value  Reference 

[A]_{tot}  5 μM  [^{73}] 
[B]_{tot}  0.28 μM  [^{49}] 
[R]_{tot}  0.16 μM  [^{49}] 
[Y]_{tot}  9.7 μM  [^{49}] 
V_{cell}  1.4 fl  [^{73}] 
N_{A,tot}  4215  calculated from above 
N_{B,tot}  236  calculated from above 
N_{R,tot}  135  calculated from above 
N_{Y,tot}  8177  calculated from above 
N_{A,tot }= NN_{C}  7027  [^{49},^{73}] 
k_{2}  10^{3 }s^{1}  [^{74}] 
k_{A}  10^{3 }s^{1}  [^{75}] 
k_{y}  100 μM^{1 }s^{1}  [^{76}] 
k_{B}  15 μM^{1 }s^{1}  [^{76}] 
k_{Y}  5 s^{1}  adjusted to yield steadystate value 
k_{B}  1.35 s^{1}  (0.35 s^{1}) [^{77},^{78}] 
γ_{R}  0.0061 s^{1}  [^{47}] 
γ_{B}  3.14 μM^{2 }s^{1}  [^{47}] 
The literature values are given in parentheses where different from our parameter values. k_{}_{Y }was determined by the condition that at steadystate with AR*=1∕3, the concentration [Y_{p}]* = [Y]_{tot}/3[^{73}].
Intensities of Gaussian white noise terms in the full pathway model.
process  index i  noise intensity Q_{i} 

receptor switching  a  2k_{2}A* 
ligand diffusion  L  2Dsc_{0} 
receptor de/methylation  M  2γ_{R}(N A*) 
CheA autophosphorylation  A_{p}  Ac*kANCN(NA,totNAp*) 
CheY phosphorylation  A, Y_{p}  kyVcell(NY,totNYp*)NAp* 
CheB phosphorylation  A, B_{p}  kbVcell(NB,totNBp*)NAp* 
CheY dephosphorylation  Y_{p}  kyNYp* 
CheB dephosphorylation  B_{p}  kbNBp* 
motor switching  X  2k+*k*k+*+k* 
Index i represents noise term η_{i}.
Fitting parameters for response function of the full pathway model for Fig. 3.
Parameter  Block et al., Segall et al. [^{26},^{27}][s^{1}]  Shimizu et al. [^{28}]  

32°C [s^{1}]  22°C [s^{1}]  
adaptation:  
λ_{1}(∂A=∂M )  0.178  0.018  0.0039 
λ_{9}  0.0263  0.0027  5.6 10^{4} 


motor switching:  
λ_{7}  4.4 10^{4}     
λ_{8}  2.111     
Motor switching rates where not adjusted when fitting to the data by Shimizu et al. [^{28}] as the highfrequency response was not measured in these experiments.
Parameters for celltocell variation in Fig. 5A.
Parameter  WT1 (black line)  red line  green line  blue line 

k+*[s1]  1.05  52.4  1.05  1.05 
k*[s1]  1.06  53.0  1.06  1.06 
γ_{R }[s^{1}]  0.0061  0.0061  6.1 10^{5}  0.0061 
γ_{B }[μM^{2 }s^{1}]  3.14  3.14  3.14 10^{2}  3.14 
N_{tot}  7000  7000  7000  70 
For additional parameters see Table 1.
Parameters for cells with low motor bias in Fig. 5B and C.
Parameter  value 

N_{tot}  4200 
Y_{tot }[μM]  2 
k+*[s1]  0.015 
k*[s1]  4.75 
∂k_{+}/∂Y_{p }[s^{1 }μM^{1}]  4.75 
∂k_{}/∂Y_{p }[s^{1 }μM^{1}]  5.5 
For additional parameters see Table 1. Similar spectra with increased lowfrequency component have been found for cells with other motor biases [^{50}].
Parameters for the fluctuationresponse theorem in Fig. 8.
Parameter  varying γ_{R }(solid)  varying γ_{B }(dashed)  varying [Y]_{tot }(dotted) 

N_{tot}  700  700  700 
γ_{R }[s^{1}]  varied  0.0061  0.0061 
γ_{B }[s^{1}]  3.14  varied  3.14 
[Y]_{tot }[μM]  9.7  9.7  Varied 
For additional parameters see Table 1.
Article Categories:

Previous Document: Ankles back in randomized controlled trial (ABrCt): braces versus neuromuscular exercises for the se...
Next Document: Kinship of longterm associates in the highly social sperm whale.