Protein connectivity in chemotaxis receptor complexes


[f]encoding = *, family = *

Protein connectivity in chemotaxis receptor complexes

Stephan Eismann1,2, Robert G Endres2,*,

1 Department of Physics and Astronomy, University of Heidelberg, Heidelberg, Germany

2 Department of Life Sciences and Centre for Integrative Systems Biology and Bioinformatics, Imperial College London, London, UK

* r.endres@imperial.ac.uk


The chemotaxis sensory system allows bacteria such as Escherichia coli to swim towards nutrients and away from repellents. The underlying pathway is remarkably sensitive in detecting chemical gradients over a wide range of ambient concentrations. Interactions among receptors, which are predominantly clustered at the cell poles, are crucial to this sensitivity. Although it has been suggested that the kinase CheA and the adapter protein CheW are integral for receptor connectivity, the exact coupling mechanism remains unclear. Here, we present a statistical-mechanics approach to model the receptor linkage mechanism itself, building on nanodisc and electron cryotomography experiments. Specifically, we investigate how the sensing behavior of mixed receptor clusters is affected by variations in the expression levels of CheA and CheW at a constant receptor density in the membrane. Our model compares favorably with dose-response curves from in vivo Förster resonance energy transfer (FRET) measurements, demonstrating that the receptor-methylation level has only minor effects on receptor cooperativity. Importantly, our model provides an explanation for the non-intuitive conclusion that the receptor cooperativity decreases with increasing levels of CheA, a core signaling protein associated with the receptors, whereas the receptor cooperativity increases with increasing levels of CheW, a key adapter protein. Finally, we propose an evolutionary advantage as explanation for the recently suggested CheW-only linker structures.

Author Summary

Receptor clusters of the bacterial chemotaxis sensory system act as antennae to amplify tiny changes in concentrations in the chemical environment of the cell, ultimately steering the cell towards nutrients and away from toxins. Despite bacterial chemotaxis being the most widely studied sensory pathway, the exact architecture of the receptor clusters remains speculative, with understanding suffering from a number of paradoxical observations. To address these issues with respect to the protein arrangement in the linkers connecting receptors, we present a statistical-mechanics model that combines insights from electron cryotomography on the linker architecture with results from fluorescence imaging of signaling in living cells. Although the signaling data for different expression levels of key molecular components in the linkers seems contradictory at first, our model reconciles these predictions with structural and biochemical data. Finally, we provide an evolutionary explanation for the observation that some of the incorporated linkers do not seem to transmit signals from the receptors.


Escherichia coli cells are able to sense changes in the chemical environment, allowing the bacteria to move towards higher concentrations of attractants and lower concentrations of repellents. The chemotaxis system is remarkable for its high sensitivity, wide dynamic range, and precise adaptation while only involving a small number of molecular components [1, 2, 3]. Despite the importance of receptor clustering in accounting for these signaling properties [4, 5, 6, 7], there are still unresolved issues with the clusters, in particular with respect to the nature of the coupling mechanism between receptors [8]. It has been proposed that receptors assemble into larger arrays via the connection of the kinase CheA and the adapter protein CheW [9, 10], with potentially complementary effects of membrane-mediated interactions [11]. Unexpectedly, in vivo Förster resonance energy transfer (FRET) shows that increasing the expression level of CheA of engineered non-adapting receptors decreases the cooperativity among receptors. In contrast, expressing more CheW increases the cooperativity, albeit in different ranges of expression levels [12]. This raises the question of how these different observations can be reconciled.

In E. coli, there are four types of methyl-accepting chemoreceptors: the high-abundance Tar and Tsr receptors that sense serine and aspartate, respectively, and the low-abundance Trg and Tap receptors [13, 14]. In addition, Aer is a chemoreceptor-like sensor of redox potential [15]. The chemoreceptors form homodimers, which assemble into trimers of dimers (TDs) [16, 17]. On a larger scale, these TDs cluster at cell poles [18, 19, 20]. CheW and CheA, which interact with the cytoplasmic domain of the receptors [21], are involved in the stabilization of these clusters [22], which in turn consist of smaller complexes (signaling teams) [6, 23, 24]. Signal transduction is triggered by ligand-receptor binding, which leads to a conformational change in the cytoplasmic domains of the receptors [25, 26, 27]. The removal of attractant (or addition of repellent) activates autophosphorylation of the kinase CheA, which is associated with the receptors via the adapter protein CheW (Fig. 1A). The phosphoryl group is then transferred to the response regulator protein CheY, which diffuses through the cytoplasm. CheY-P binds to the flagellar motors to induce clockwise rotation and tumbling of the cell. In contrast, addition of attractant (or removal of repellent) inhibits autophosphorylation of CheA. CheY-P dephosphorylation by phosphatase CheZ leads to counterclockwise rotation and straight swimming [1].

Figure 1: Schematic of receptor clustering in E. coli. (A) Association and dissociation of adapter protein CheW and kinase CheA2 with the complexes. CheW interacts directly with the receptors, but the interaction of CheA with the receptors is largely mediated by CheW (see below). (B) Electron cryotomography images show that trimers of dimers (TDs) of chemoreceptors cluster at the cell poles in a hexagonal manner [28, 29, 30, 31]. (C) Nanodisc experiments propose that two TDs are connected by a linker consisting of two CheW monomers and one CheA2 dimer [9]. In reality, the P5 domain of CheA may also contact the trimer [32, 31, 21], although this binding may be an order of magnitude weaker than CheW-trimer binding [9, 29]. (D) Top view of an ensemble of different sizes of receptor complexes in the cytoplasmic membrane. Active receptors are shown in red, and inactive receptors are shown in blue. Each linker between active TDs contributes a coupling energy .

To avoid saturation of the sensory system, adaptation is implemented via covalent receptor modification. This is achieved through changing the receptor-methylation level by the activities of the methyltransferase CheR and the methylesterase CheB, which antagonistically add and remove, respectively, methyl groups at four or five, depending on the receptor, specific glutamate residues on each receptor monomer [33], respectively. Methylation by CheR increases the activity of CheA, i.e., its autophosphorylation rate, thus counteracting the effect of attractant binding. In contrast, CheB activation by phosphorylation by CheA-P decreases CheA activity [12]. Through genetic engineering, the glutamate residues (E) can be replaced by one to four glutamine residues (Q) to mimic increasing receptor-methylation levels in the absence of CheR and CheB [2]. The E. coli chemotaxis pathway is exceptionally well characterized and is thus amenable to modeling at a high quantitative level.

To explain the receptor cooperativity, which generates the high sensitivity of the system, the mechanism of receptor-receptor coupling has attracted much interest [34, 8, 35, 36, 37]. Electron cryotomography (EC) images of the TDs in quick-frozen cells led to the idea that TDs form densely packed hexagonal ‘honeycomb’ arrays (Fig. 1B) [28, 32, 38]. These and other in vitro experiments using nanodiscs and nanoscale plugs to imitate cellular membranes suggest that –CheW–CheA2–CheW– is the structural core unit linking two TDs (see Fig. 1C for a simplified depiction) [9]. An approach to study the cooperative behavior of the specific receptors inside the cells indirectly is to monitor the signaling activity of CheY-P/CheZ pairs via FRET, with the FRET signal being proportional to the overall CheA activity [39]. An increase in the concentration of CheW was observed to enhance the cooperativity of the FRET response mechanism, whereas, unexpectedly, an increase in CheA concentration led to the opposite effect [12]. It is well known that multimeric protein complexes can be inhibited by high concentrations of one of their components, similar to the prozone phenomenon in precipitin tests [40]. However, it is unclear how the FRET results relate to other experimental observations, including the proposed linker and lattice structures.

Here, we use statistical-mechanics modeling within the framework of the Monod-Wyman-Changeux (MWC) model [41] for cooperative receptor complexes to unify the assumed linker and lattice structures with the seemingly contradictory FRET results. By implementing the linker structure we initially fit our model of receptor complexes of up to four TDs to FRET data obtained with cells that express only the Tar receptor in different non-adapting modification states. Next, we apply our model to Tar–Tsr–Tap and Tsr–only cells in the non-adapting QEQE modification state, which mimics half-methylated receptors. As a result we recover the experimentally observed decrease in cooperativity of the response to serine with increasing CheA concentration, whereas increasing CheW yields the observed enhanced cooperativity. Note, other higher order effects of protein overexpression, such as membrane invaginations or interference of CheA/ CheW with clustering, are not included. Our results surmise that the observed opposing trends in cooperativity are based on a critical combination of the correct linker architecture and a constant average complex size.


Statistical-mechanics model of chemotaxis receptors

At the heart of our approach lies the MWC model [5, 6, 12]. Chemoreceptors are regarded as two-state systems being either active (on) or inactive (off), with conformation-dependent dissociation constants and for a specific ligand. As the attractant affinity of inactive receptors is higher than for active receptors (), the state ratio tips towards inactive receptors with increasing ligand concentration . In contrast, receptor modification favors the active state in the absence of ligands represented by an energy offset . The resulting single-dimer free energies in the active and inactive states are given by


with the chemical potential of the receptors in the membrane. All energies are expressed in units of the thermal energy, . In our approach, we allow for an ensemble of different complexes with varying complex size (i.e. number of connected TDs) and partially developed linkers as rest groups (Fig. 1D). All receptors within a complex are assumed to share the same conformational state because of tight coupling. For simplicity, we consider the –CheW–CheA2–CheW– linker structure [9], which we incorporate by assigning energies and for each CheW and CheA2 molecule integrated in a specific receptor-complex type (see Discussion section for an alternative linker structure). These energies are of the forms


where and indicate monomer concentrations and and are dissociation constants for CheW–receptor and CheW–CheA2 binding, respectively. In particular and are expressed as fractional changes and of wild-type expression levels and , respectively:


The TD is assumed to be the smallest receptor unit [42, 9], and the maximal number of connected TDs is restricted to four, in line with observed Hill coefficients from FRET [12, 23]. (Including larger complex sizes does not alter the model predictions, but increases the computational complexity significantly; see Materials and Methods.) Each dimer can maximally bind to one molecule of CheW, whereas CheA is assumed to not interact with receptor dimers directly. In order to restrain the combinatorial complexity partially developed linkers are only considered in a symmetric manner, i.e. all rest groups are assumed to be identical in a complex. Furthermore, we attribute an attractant energy to each linker within an active complex, a treatment in line with the previously proposed enhanced coupling among active receptor dimers [24], albeit independent of receptor-modification level.

The resulting free energies for a complex of size and rest group are given by (cf. Fig. 1D)


with receptor dimers per complex of size and linkers. Such a complex has rest groups with given by


for (1) no rest group, (2) a CheW and (3) a CheW and a CheA dimer, respectively. The probability for a certain complex type and its probability of being active follow from standard combinatorial reasoning and the partition function


where the number 1 in the partition function reflects the possibility of an empty membrane site.

Assuming the FRET signal to report the number of CheA2 dimers within an active complex, we define the receptor activity as


In contrast, the classical MWC model for coupled receptors describes the response of a single complex of TDs to a change in ligand concentration. Without incorporating the receptor coupling explicitly, the corresponding activity reads [23]


In the past, the Hill coefficient and complex size have broadly been treated as equivalent to quantify the cooperative behavior of receptor complexes, and in [23], an increase in with receptor-modification level was equated with an increase in receptor cooperativity. However, both quantities are not necessarily the same as approximating Eq. 10 by a Hill function with requires [6]. We found that, in the classical MWC model, the response of differently modified Tar receptors to MeAsp, a non-metabolizable analog of aspartate, can also be described with a fixed for all modification levels. This treatment results in a similar quality of fit when relating the reduced number of parameters to the new goodness-of-fit value (see S1 Fig..). As our model incorporates an ensemble of complexes of varying sizes, the finding of a constant complex size in the classical MWC model is naturally generalized by a constant average complex size with respect to ligand concentration and receptor-modification state. The average complex size, which we term receptor density , is given by


with being the number of dimers of a given complex type . The chemical potential in Eq. 1 is adjusted throughout the simulation to fulfill this condition, reflecting anticipated regulation of the receptor-expression level by the cell. Biologically, a constant receptor density can be achieved by random receptor insertion into a growing membrane at constant rate [35]. Since wild-type cells express and insert receptors in the QEQE modification state [2], we do not expect a modification-dependent insertion rate. Although allowing for a modification-dependent would increase the quality of fit because of an increased number of fitting parameters, our minimal model with constant can describe the data very well.


Receptor-modification level may not determine cooperative behavior of complexes

In order to test our model, we firstly applied it to FRET data of Tar-only receptors in different non-adapting receptor-modification states from Ref. [23] i.e. Tar{QEQE}, Tar{QEQQ} and Tar{QQQQ}. The dose-response curves of the chemoreceptors match closely the statistical-mechanics model with fixed receptor density, and hence fixed average complex size (Fig. 2A). Figure 2B displays the fitted receptor density next to the Hill coefficients of the experimental curves (see Materials and Methods) and the complex size of the classical MWC model, taken from [23]. Although the classical MWC model predicts a rise in complex size with modification level [23], including its implementation based on a dynamic Ising model [24], this is not true for the Hill coefficients (see also S1 Fig..). This finding shows that receptor modification is not the main determinant of receptor cooperativity.

Figure 2: Kinase activity for different Tar-modification levels. (A) Kinase activity for Tar receptors in QEQE (black), QEQQ (blue) and QQQQ (green) modification states as a function of MeAsp concentration. The curves are normalized with respect to QQQQ activity at concentration mM. (B) All model curves share the same receptor density (light green), which reflects the average complex size, depicted next to the corresponding Hill coefficients (dark green). Parameter (yellow) of the classical MWC model (taken from [23]) is shown for comparison. (C) The chemical potential , shown as a function of ligand concentration for the three modification levels, is adjusted throughout the simulation to ensure constant at all concentrations. Color coding is the same as in panel A. The vertical dashed lines indicate half-maximum concentrations from the corresponding Hill fits. While for is approximately constant, the curves follow a logarithmic function in the regime of . The dotted blue line shows for comparison. Model parameters: , , , , , , , , and . Values for , , and are shared with curves shown in Figs. 4 and 5. The superscript 0 indicates wild-type expression levels for CheA/CheW.

In our model, the chemical potential can be regarded as the cost function for the cell to provide a constant complex size in the membrane. By definition, the chemical potential reflects the amount of energy required for adding a particle to a system with free energy . Although the value of the parameter , introduced to ensure constant receptor density , is gained by solving a highly nonlinear equation, its behavior with respect to ligand concentration is very homogeneous and characterized by two regimes, as shown in Fig. 2C. While this cost is approximately constant for , with being the half-maximum concentration obtained from Hill fits, the cost necessary to maintain a constant density increases rapidly for ligand concentrations beyond . In this second regime, the curves for all modification levels are of the form , which is the functional description of an ideal chemical potential. Although the slope in the second regime is the same for all values of , the different offsets reflect the modification-dependent energy . Note, if we were instead to keep constant (and not ), then bumps would appear in the dose-response curves as a result of the receptor density increasing with ligand concentration (see S2 Fig..).

In summary, our model is capable of quantitatively describing dose-response curves from in vivo FRET, in particular the receptor-receptor cooperativity. Although in spirit similar to other recent statistical-mechanics models, most noticeably by Hansen et al. [24] and Lan et al. [43], only our model addresses the protein connectivity in receptor complexes.

Receptor density governs cooperative behavior of complexes

While the receptor density is assumed to be constant on a short time scale, the rate of receptor expression and insertion into the membrane can be regulated by the cell on a longer time scale. Hence, as a further test of our statistical-mechanics model, we investigated how a change in receptor density affects CheA activity at wild-type expression levels for CheA and CheW. Figure 3A shows modeled dose-response curves for different values of , and with the wild-type receptor density and otherwise using the same parameter set as in Fig. 2. An increase in receptor density is directly associated with an enhanced signal amplitude because more CheA molecules are incorporated into the complexes. Figure 3B reflects the associated trend in cooperativity by comparing density and Hill coefficient . In qualitative agreement with experimental observations [12] and in line with previous modeling [6], larger complex sizes lead to higher sensitivities and hence steeper dose-response curves given a certain receptor-modification state. Since the expression level of receptors (and other chemotaxis proteins) is highest under nutrient-poor conditions, the resulting increase in receptor density and cooperativity leads to enhanced sensitivity when it is most crucial for cell survival [44].

Figure 3: Cooperativity increases with receptor density. (A) Model curves based on previously fitted parameters of Tar–only cells in QEQE modification state (Fig. 2) for different receptor densities (black), (blue) and (green) with the wild-type receptor density . (B) Hill coefficients derived by fitting to the model outcome and corresponding receptor densities .

Increasing the CheW level increases receptor cooperativity

To gain insight into the role of CheA and CheW in forming receptor complexes, we varied the expression levels [A] and [W] to study the effect on receptor activity. According to the experimental observations in [12], we set the CheW concentrations to 0.7, 0.1 and 0.01 and the CheA concentrations to 8, 0.3 and 0.25 times the wild-type values and , respectively. This allowed us to make the comparison with experimental dose-response curves from FRET of Tsr–only cells (for varying CheW) and Tar–Tsr–Tap cells (for varying CheA), both in the non-adapting QEQE modification state. To keep the overall number of parameters small, the data for changes in [A] and [W] was fitted with the same parameter set (, , , , , and ). Multiplication of the calculated activities with scaling parameters and , respectively allows for comparison with the FRET signal amplitudes. Subsequently, a Hill function was fitted to the model curves and the model Hill parameters were compared with the experimental values. Note that our minimal model does not account for alternative forms of signaling disruption upon over- or underexpression of CheA/CheW, such as zipper-like invaginations of the cell membrane [45] or interference with trimer formation [16].

Figures 4A,B show the model data next to the experimentally determined Hill curves for variations in [W]. Enhanced CheW expression results in raised activity amplitudes and Hill coefficients (Fig. 4C,D). Although the values from the model change significantly with expression level [W] at a confidence level, which is in qualitative agreement with the experimental data, especially with respect to the highest CheW expression level, the change in is less pronounced for the model than the experimental data. The positive correlation between kinase activity and amount of available CheW becomes evident in the distribution of complex species at half-maximum concentration (Fig. 4E). Whereas low levels of [W] favor independent, single TDs, larger complexes are more likely to form for larger [W]. As the probability for an empty membrane site also increases, the receptor density remains constant.

Figure 4: Cooperativity increases with the expression level of CheW. (A,B) Model fit (A) and FRET data [12] (B) for different expression levels of CheW at 0.7 (dark blue), 0.1 (blue) and 0.01 (azure) times the native level of . (C) Hill curves were fitted to the model outcome to allow for comparison with experimental results. Hill coefficient pairs (model/experiment) in order of increasing [W] are (3.0/3.0), (4.2/5.0) and (4.6/9.6). (D) The corresponding Hill amplitudes in order of increasing [W] are (0.005/0.006), (0.033/0.036) and (0.049/0.044). (E) Distribution of complex types present at half-maximum concentration for 0.7 (dark blue) and 0.01 (azure) times the native concentration . Model parameters: , , , , , , and . Parameters are shared with model for variation in [A] (Fig. 5).

Increasing the CheA level decreases receptor cooperativity

Changing [A] in our model has the opposite effect on the Hill coefficient as changing [W]. This result is in line with experimental data (Fig. 5A,B,C). The activity amplitude reflecting the amount of active CheA molecules benefits from higher CheA levels, as one would expect (Fig. 5D). In contrast, Hill coefficients are higher for smaller [A], recovering the naively unexpected experimental observations (Fig. 5C). Looking at the distribution of complexes at half-maximum ligand concentration (Fig. 5E), we note that although high CheA concentrations favor rest groups including CheA, complex sizes of 3 and 4 TDs are more likely at lower concentrations of CheA.

Figure 5: Cooperativity decreases with the expression level of CheA. (A,B) Model fit (A) and FRET data [12] (B) for different expression levels of CheA at 8 (red), 0.3 (orange) and 0.25 (yellow) times native level . (C) Hill curves were fitted to the model outcome to allow for comparison with experimental results. Hill coefficient pairs (model/experiment) in order of increasing [A] are (5.5/5.0), (5.5/5.0) and (3.4/2.0). (D) The corresponding Hill amplitudes in order of increasing [A] are (0.023/0.022), (0.026/0.026) and (0.030/0.031). (E) Distribution of complex types present at half-maximum ligand concentration for 8 (red) and 0.25 (yellow) times native concentration . Model parameters: , , , , , , and . Parameters are shared with model for variation in [W] (Fig. 4).

The opposing trends in concerning variations in [A] and [W] are a direct result of the linker stoichiometry and fixed average complex size. For complexes with rest groups, the ratio of CheW molecules per TD is independent of the complex size (Fig. 6A). However, for species without rest groups, this ratio increases with the number of coupled TDs. As a result, an enhancement in [W] yields larger complexes that directly incorporate more CheA molecules. Furthermore, empty sites ensure a constant receptor density even when expression levels of CheW and CheA are extremely low. In this case, the receptor density still remains constant as empty sites can be occupied by individual TDs. This requires a dilute membrane, i.e., a receptor density not much larger than (see Fig. 2B).

In contrast, the corresponding ratio of CheA dimers per TD is highest for single TDs with full rest groups and decreases with increasing complex size (Fig. 6B). The CheA molecules within the rest groups contribute to the FRET amplitude but not to the receptor cooperativity. An accompanying rise in the number of occupied membrane sites ensures a constant receptor density.

Figure 6: CheW and CheA counts per TD show different trends. (A) For complexes with rest groups, the number of CheW molecules per TD is independent of complex size. For species with no rest groups, this ratio increases with the number of TDs in the species. Raising the expression level of CheW results in the formation of larger complex sizes. (B) The number of CheA dimers per TD shows two opposing trends with respect to complex size. As the ratio is highest for single TDs with CheA2-including rest groups, raising the expression level of CheA results in smaller complex sizes but an increased number of CheA molecules contributing to the signal amplitude.

Electron cryotomography suggests the existence of CheW-only linkers

Our model qualitatively reproduces the experimental results obtained when the expression levels of CheW and CheA were changed. However, there are quantitative differences, especially with respect to the change in cooperativity as a function of the expression level of CheW. This change is less pronounced in the model than in the experiment. Recent findings from electron cryotomography may shed light on the reasons for these discrepancies. Although both studies stressed the importance of one dimeric CheA and two CheWs as the minimal unit needed for kinase activation, Briegel et al.[30] and Liu et al.[31] proposed additional CheW-only linkers, underlining the role of CheW in the cooperative behavior of TDs. Such structures could explain how increased levels of CheW contribute to the cooperativity of TDs. In order to quantify this effect, we allowed for additional CheW-only linkers in our model (Fig. 7). The dimeric appearance of CheW in the linker is accounted for by a new parameter ; we keep the previously introduced rest groups for simplicity.

Figure 8 shows the results for varying expression levels of CheW and CheA. The dose-response curves of the new model exhibit the same trends in Hill coefficient and amplitude for variation in [W] (Fig. 8A) and [A] (Fig. 8B) as before, in agreement with experimental results (see also S3 Fig..). However, the difference in behavior is manifested in the comparison panels below. The previously obtained minor changes in receptor cooperativity as a function of [W] are now much more pronounced (Fig. 8C), although the modeled Hill coefficients for [A] variation are larger than the experimental ones (Fig. 8D). The excess CheW leads to formation of CheW-only linkers and hence larger complex sizes when the amount of available CheA is held constant.

Figure 7: CheW-only linkers. (A) In agreement with recent electron cryotomography experiments [31, 30], we allow for additional CheW-only linkers (yellow) connecting two TDs (black and blue). (B) Exemplary ensemble of complexes in the cytoplasmic membrane. The two linkers are represented by solid (–CheW–CheA2–CheW–) and dashed (–CheW–CheW2–CheW–) lines. Active and inactive TDs are shown in red and blue, respectively. Each linker between active TDs contributes an additional coupling energy .
Figure 8: CheW-only linkers can explain the large enhancement in cooperativity with increasing [W]. (A,B) Modeled dose-response curves for different levels of expression of CheW (A) and CheA (B) as multiples of wild-type levels and , respectively. (C,D) Experimental results in gray are superimposed with parameters inferred from Hill curves fitted to the model outcome. Parameters for variation in [W] and [A] are shown in panels (C) and (D), respectively. Model parameters: , , , , , , , , and . The data for variations in CheA and CheW levels was fitted with receptor densities and , respectively.

In order to make predictions beyond the data used to fit the model, we created surface plots of amplitudes and Hill coefficients covering several orders of magnitude for expression levels of CheW and CheA (Fig. 9). The receptor activity and hence amplitude increases monotonically with the level of CheA, whereas the increase in amplitude with respect to the level of CheW is only pronounced in a subspace around the experimental data (Fig. 9A). In the case of high CheA levels, CheW-only linkers exclude CheA from signaling. This also occurs at the wild-type CheA level, although the extent of the effect strongly relies on model parameters. The surface plot showing the Hill coefficients as a function of the expression levels of CheW and CheA has a saddle-like form (Fig. 9B). Although the right flank is consistent with the FRET data at high levels of CheA (small Hill coefficients), the Hill coefficient also decreases at very low levels of CheA as the receptor activity diminishes. To test to what extent the model predictions depend on the actual values of parameters , and , we varied these parameters and found that the general shape of the surface plot was preserved. Taken together, these observations suggest the need for regulation of both CheW and CheA by the cell to balance signaling amplitude and sensitivity. Indeed, as CheW and CheA are required in comparable amounts [13], both are expressed from the same operon [46].

Figure 9: Varying CheW and CheA expression levels over a wide range shows trends beyond the FRET data. (A) Surface plot showing the amplitudes of simulated dose-response curves for different expression levels of CheW and CheA (in units of wild-type levels and , respectively). (B) The corresponding surface plot for the Hill coefficient has a saddle-like form. Simulations were performed using the parameter set of Fig. 8 with .


Receptor coupling plays a key role in the remarkable sensing and signaling properties of bacterial chemotaxis. These networks can explain the high sensitivity, wide dynamic range and precise adaptation. In this work we present a statistical-mechanics model of different complex sizes, modeling for the first time a molecular linker architecture consistent with (i) FRET dose-response curves, (ii) cryotomography data and (iii) nanodisc experiments. The linker –CheW–CheA2–CheW– proposed by Li and Hazelbauer [9] is incorporated by assigning expression level-dependent energies and respectively for each CheW and CheA2 molecule within a complex as part of a fully or partially developed linker. A coupling energy attributed to linkers between active TDs indicates that the coupling between active trimers is stronger than between inactive trimers, in agreement with previous modeling [24]. Although the actual distribution of complex sizes is influenced by expression levels [W] and [A], a readily adapted chemical potential ensures a fixed average complex size with respect to ligand concentration .

Our model was first applied to describe the dose-response of Tar receptors in different modification states to MeAsp, a non-metabolizable analog of aspartate. We mainly considered a constant, modification-independent , a constraint that not only reduces the number of parameters but also calls into question that the complex size increases with receptor-modification level [23]. In our work we discovered the discrepancies between the number of connected TDs and the curves’ Hill coefficients within the classical MWC model. An increase in is not directly associated with an increase in . In our statistical-mechanics model, the approximately constant is explained by a constant average complex size across all receptor-modification levels. Indeed, experiments show that both the level of expression of receptors and the insertion of newly synthesized receptors into the inner membrane by the Sec-machinery are highly regulated [47, 48].

Hansen et al. [24] previously presented a dynamic-signaling-team approach to describe the data obtained with Tar-only cells in which the allosteric coupling among trimers is represented by a modification-dependent trimer-trimer interaction energy without modeling the actual protein connectivity. Limited conformational spread and hence a finite complex size is achieved by using a long-range repulsion energy between all trimers within a complex. In contrast, our model is simpler while providing valuable insights. Neither and nor in our ensemble model depend on the modification state of the receptor, and ensures constant average complex size without introducing a repulsive term. Furthermore, the chemical potential provides insights into the energetic cost of insertion of receptors into the membrane and its dependence on ligand concentration , albeit based on an equilibrium mechanism.

For constant and , we conclude that receptor modification mainly governs the ‘turn off’-ligand concentration, whereas its influence on receptor clustering is limited. This finding is supported by Briegel et al.[49], who found that the receptor array order and the spacing of receptors in different modification states were indistinguishable. This is in stark contrast to Hansen et al. [24], who predict a strong increase in average complex size with increasing receptor-modification level. High-resolution imaging of equilibrated receptors in artificial membranes by electron or total internal reflection fluorescence (TIRF) microscopy may allow direct determination of receptor-complex distributions and their dependence on receptor-modification level and ligand concentration. Using photoactivated localization microscopy (PALM) [35] or quantitative immunoblotting [13], such an investigation could also be performed on intact cells.

Although CheA and CheW have long been known to mediate receptor interactions [12, 22], an increase in the expression level of CheA leads to a reduction in receptor cooperativity [12]. Varying expression levels of CheA and CheW in our model produced results in agreement with experimental data of Sourjik and Berg [12], thereby supporting the linker architecture we employed. The striking observation that increased CheA levels lead to higher kinase activities but lower cooperativity is based on the fact that the number of CheA dimers per TD is highest for single trimers with almost fully developed linker rest groups (Fig. 6B). Hence, overexpression of CheA, a bridging molecule at the center of the linker, promotes smaller complex sizes. CheA molecules within the rest groups do not contribute to TD coupling and curve steepness, but nevertheless add to the activity of the FRET signal.

In contrast to what is observed with CheA, raising the level of CheW leads to larger complex sizes and an increased number of empty membrane sites. Again, this behavior becomes comprehensible when the number of CheW molecules per TD (Fig. 6A) is taken into account. While this ratio is constant for complexes with rest groups, it increases with complex size in the absence of partially developed linkers. Larger complexes directly incorporate more CheA to enhance cooperativity as well as the amplitudes of FRET signals observed both in the model and experimentally. In light of our model the experimental observations are produced by a combination of constant receptor density and (partial) linkers. Although partial linkers play a crucial role in the mechanism of our model, their inclusion might appear arbitrary at first. Interestingly, Briegel et al.[30] recently observed a range of assembly intermediates and partial receptor hexagons forming when [W] and [A] were varied. Our surface plots of amplitudes and Hill coefficients also make testable predictions for wide-ranging CheA and CheW expression levels (Fig. 9). Is there any evidence to suggest that remains constant when CheA and CheW expression levels change? First, CheA and CheW binding to the receptors occurs after insertion of the receptors into the membrane. Second, increasing the expression of a protein, e.g., of CheW, should remove ribosomes from translating receptor mRNA [50, 51]. Although expected to be a minor perturbation, this may lead to a reduced receptor density and hence cooperativity. However, the opposite trend is observed in FRET experiments [12].

Although our assumed linear linker structure –CheW–CheA2–CheW– matches observed stoichiometries [9, 13], electron cryotomography images suggest that reality is more complicated [30, 52]. Modeling of the electron density and spin-labeling studies suggest that CheW and the P5 domain of CheA form alternating CheW/CheA rings connecting the trimers, with P5 occupying positions approximately equivalent to CheW (see Fig. 10). This arrangement is consistent with the strong structural homology between P5 and CheW. However, to describe the FRET data obtained with cells with overexpressed CheA and CheW [12], our model predicts that CheA2 has the role of a bridging molecule and connects trimers via a CheW associated with each trimer. Indeed, an alternative linker with direct receptor-CheA binding and hence symmetric roles of CheA and CheW upon clustering does not match the FRET data (see panel D in S4 Fig..). This view is supported by binding assays, which show that CheW binds much firmer to receptor trimers than CheA to trimers (see Fig. 5A,B in [9] and also discussion in [29]).

Figure 10: Structural insights from electron cryotomography. (A) Our linker –CheW–CheA2–CheW– based on nanodisc experiments (black line) [9] appears to be more complicated in reality, where the P5 domain of CheA (homologous to CheW) may also contact the trimers of dimers directly (contact 1 in red) [21, 32, 31]. Such an alternative linker, defined as =CheW/CheA2/CheW= is explored in S4 Fig.. However, the binding of CheW to the trimers (contact 2 in orange) is presumably much stronger (see Fig. 5A,B in [9] and discussion in [29]), rendering CheA2 effectively a bridging molecule. (B) Hexagonally packed trimer-of-dimers structure in which the inner connecting ring is formed by alternating CheW/CheA (P5) units [28, 29, 30, 31, 52].

Although our model qualitatively reproduces the experimental FRET data, the change in cooperativity with variation in [W] is less pronounced in the simulation than in experiments. Recent findings based on electron cryotomography offer a possible explanation for this shortcoming. Briegel et al.[30] and Liu et al.[31] stress the importance of the implemented core unit stoichiometry, but they propose a second type of linker that only involves CheW, with P5/CheW interactions replaced by CheW/CheW interactions [31]. To investigate the consequences of these findings for signaling behavior, we allowed for an additional –CheW–CheW2–CheW– linker in our model. The simulated dose-response curves show a greatly enhanced change of cooperativity with variation in [W] (Fig. 8C and Fig. 11A). The generally increased Hill coefficients, and hence sensitivity, may reveal an evolutionary advantage that is not apparent in the tomography images but is detected by FRET. However, whereas CheW-only linkers fit the FRET observations, their incorporation into complexes needs to be tightly regulated. Moreover, in addition to excluding CheA from signaling (Fig. 11B), high levels of CheW were also claimed to disrupt receptor clustering [53]. Taken together these observations suggest that an optimal level of CheW is required for cooperative signaling by receptors (Fig. 11C).

In conclusion, our work integrates functional (FRET) and structural (nanodisc and electron cryotomography) data, explains the paradoxes that increased levels of CheA lead to less cooperativity, and provides a functional role for CheW-only linkers. Our proposed linker –CheW–CheA2–CheW– is consistent both with the data from experiments with nanodiscs [9] and with images from electron cryotomography [31, 29, 30], if the P5 domain of CheA binds more weakly to the receptor than does CheW. We predict that the observed tetrameric CheW linker, if incorporated at an optimal level, increases the cooperativity while keeping the receptor activity at a sufficiently high level. An increased understanding of the protein connectivity in receptor clusters may aid not only in describing the fundamental biology of receptor signaling, including the role of cytoplasmic receptor clusters in Rhodobacter sphaeroides and Vibrio cholerae [52], but may also contribute to the design of novel biosensors [54].

Figure 11: Simulations suggest that an optimal level of CheW is required for highly cooperative signaling. (A) Modeled Hill coefficients for different expression levels of CheW as multiples of its native level . CheA is modeled at its native level for all charts. (B) The relative number of CheA dimers per linker is simulated for the different levels of CheW expression. Results are compared to the native level . (C) Although increased levels of CheW lead to larger clusters, formation of CheW-only linkers also excludes CheA from signaling. These findings suggest that an optimal CheW level is required to balance signaling sensitivity and magnitude. Model parameters are as in Fig. 8.

Materials and Methods

Keeping constant requires nonlinear optimization of at every ligand concentration. For performance reasons we therefore chose to implement the model in C# and used a custom-written toolbox to connect to MATLAB 2014a for parameter optimization and plotting. The value for is determined based on Brent’s method for root-finding [55]. Fitting of model parameters employs Global Search from MATLAB Global Optimization Toolbox. Multiple start points are generated using scatter-search options (5000 trial points). For the different start points square deviations from experimental data are minimized using the function fmincon with interior point optimization. Note while the number of molecular species in the model increases linearly with the maximal complex size, the computational time is determined by the root finding. The latter becomes considerably harder with additional exponentials of increasing arguments in Eqs. 6, 7 and 11.

In order to quantify the cooperative behavior of the complexes, Hill functions (Eq. 12) with amplitude , half-maximum concentration and Hill coefficient are fitted to the model evaluated at 50 logarithmically spaced concentrations between mM and mM. The Hill coefficients in the comparative plot Fig. 2B result from direct fitting to the experimental data.


Though parameter confidence intervals can be calculated based on robust regression and the resulting covariance matrix, especially for highly nonlinear models as ours their validity is questionable given the underlying linear theory [56]. We therefore decided against including confidence intervals except for the fitted Hill curves.

We note that for all simulations with variations in expression of CheA and CheW the Hill amplitudes match quantitatively much better their experimental counterparts than do the Hill coefficients. This observation is partly owed to the fitting routine. With logarithmically spaced concentrations, a difference in amplitude between model and experimental curve directly impacts the corresponding goodness-of-fit value. In contrast, a small variation in the Hill coefficient only influences the slope of the curve within a relatively narrow range of ligand concentrations and hence is less reflected in the optimization function value.

Supporting Information

S1 Fig.

Kinase activity for different Tar-modification levels can be described with a constant receptor-complex size. (A,B) Kinase activity for Tar receptors in QEQE (black), QEQQ (blue) and QQQQ (green) modification states fitted in the classical MWC model with (A) constant and (B) variable receptor-complex size . The fitting based on Principal Component Analysis (PCA) follows Ref. [23]. Relating the resulting values to the degrees of freedom, here calculated as the number of included PCA components minus the number of model parameters, results in similar goodness-of-fit values with subscript ‘red’ describing the reduced . However, it should be noted that the actual here is rather a supportive argument to the apparent similarity of both fits, as the number of degrees of freedom is not well defined for nonlinear models [57]. (C) Comparison of receptor-complex size and Hill coefficient for fits with constant (left, panel A) and variable (right, panel B).

S2 Fig.

Receptor density increases with ligand concentration for constant chemical potential. (A) Receptor density as a function of ligand concentration for (solid) and (dashed). For a constant chemical potential , the values of the single dimer energies and (Eq. 1) decrease with increasing ligand concentration c. The decrease in the resulting complex energies and (Eq. 4) is stronger for larger complexes. Hence, larger complexes are favored with increasing ligand concentration (Eq. 7), resulting in an increased receptor density (Eq. 11). The interim decrease in for is the result of an ensemble effect. While the probabilities of all complex sizes increase with c, the increase for larger complexes starts at higher c values. Starting off at a smaller receptor density, this effect is not visible for . Finally both densities asymptotically approach the maximal value of 12. (B) Normalized kinase activity as a function of ligand concentration for (solid) and (dashed). In the case of , the increase in receptor density is not apparent in the dose-response curve as the receptors ‘turn off’ before the density increase comes into effect. For , however, the increase in receptor density yields a ‘bump’ in the dose-response curve. All plots were generated using the same parameters for QEQE as in Fig. 2.

S3 Fig.

Kinase activity for different Tar-modification levels with additional CheW-only linkers. Plot following Fig. 2A showing kinase activity for Tar receptors in QEQE (black), QEQQ (blue) and QQQQ (green) modification states. Here the model includes both linkers (–CheW–CheW2–CheW– and –CheW–CheA2–CheW–). For simplicity parameters are the same as in Fig. 2 with the additional value for in agreement with the value used in Fig. 8.

S4 Fig.

Alternative model with both CheW and CheA binding to trimers does not explain FRET data. (A) Schematics of an alternative linker =CheW/CheA2/CheW= with both CheA and CheW contacting the trimers directly. (B) Exemplary ensemble of complexes in a membrane. The two linkers are represented by solid (=CheW/CheA2/CheW=) and dashed (=CheW/CheW2/CheW=) lines. Active and inactive TDs are shown in red and blue, respectively. Each linker between active TDs contributes a coupling energy . As monomeric CheA binds directly to trimers, all linker molecule energies are indicated for monomers, hence the linker energy contributions in Eqs. 4 become (standard linker) and (CheW-only linker). (C,D) In analogy to Fig. 8, we fitted the alternative model to the experimental data for varied expression levels of CheW and CheA using a global optimization routine (see Materials and Methods). While the alternative model is qualitatively able to describe the effect of changing CheW levels correctly (C), it falls short of reproducing the cooperativity decrease for increasing CheA levels (D) with nearly identical curves as best fitting result. Model parameters: , , , , , , and .


The authors thank Moritz Beutel for designing the custom-written toolbox connecting C# and MATLAB, and Thomas Rogerson and Victor Sourjik for helpful discussions and a critical reading of the manuscript.


  1. Sourjik V, Armitage JP. Spatial organization in bacterial chemotaxis. EMBO J. 2010 Aug;29(16):2724–2733.
  2. Sourjik V, Berg HC. Receptor sensitivity in bacterial chemotaxis. Proc Natl Acad Sci. 2002 Jan;99(1):123–127.
  3. Endres RG. Physical principles in sensing and signaling. 1st ed. Oxford University Press; 2013.
  4. Bray D, Levin MD, Morton-Firth CJ. Receptor clustering as a cellular mechanism to control sensitivity. Nature. 1998 May;393(6680):85–88.
  5. Mello BA, Tu Y. An allosteric model for heterogeneous receptor complexes: Understanding bacterial chemotaxis responses to multiple stimuli. Proc Natl Acad Sci. 2005 Nov;102(48):17354–17359.
  6. Keymer JE, Endres RG, Skoge M, Meir Y, Wingreen NS. Chemosensing in Escherichia coli: two regimes of two-state receptors. Proc Natl Acad Sci. 2006 Feb;103(6):1786–1791.
  7. Shimizu TS, Tu Y, Berg HC. A modular gradient-sensing network for chemotaxis in Escherichia coli revealed by responses to time-varying stimuli. Mol Syst Biol. 2010 Jun;6(1):382.
  8. Endres RG. Polar chemoreceptor clustering by coupled trimers of dimers. Biophys J. 2009 Jan;96(2):453–463.
  9. Li M, Hazelbauer GL. Core unit of chemotaxis signaling complexes. Proc Natl Acad Sci. 2011 Jun;108(23):9390–9395.
  10. Bhatnagar J, Borbat PP, Pollard AM, Bilwes AM, Freed JH, Crane BR. Structure of the ternary complex formed by a chemotaxis receptor signaling domain, the CheA histidine kinase, and the coupling protein CheW as determined by pulsed dipolar ESR spectroscopy. Biochem. 2010 Mar;49(18):3824–3841.
  11. Haselwandter CA, Wingreen NS. The role of membrane-mediated interactions in the assembly and architecture of chemoreceptor lattices. PLoS Comp Biol. 2014 Dec;10(12):e1003932.
  12. Sourjik V, Berg HC. Functional interactions between receptors in bacterial chemotaxis. Nature. 2004 Mar;428(6981):437–441.
  13. Li M, Hazelbauer GL. Cellular stoichiometry of the components of the chemotaxis signaling complex. J Bacteriol. 2004 Jun;186(12):3687–3694.
  14. Neumann S, Hansen CH, Wingreen NS, Sourjik V. Differences in signalling by directly and indirectly binding ligands in bacterial chemotaxis. EMBO J. 2010 Oct;29(20):3484–3495.
  15. Bibikov SI, Biran R, Rudd KE, Parkinson JS. A signal transducer for aerotaxis in Escherichia coli. J Bacteriol. 1997 Jun;179(12):4075–4079.
  16. Studdert CA, Parkinson JS. Crosslinking snapshots of bacterial chemoreceptor squads. Proc Natl Acad Sci. 2004 Feb;101(7):2117–2122.
  17. Boldog T, Grimme S, Li M, Sligar SG, Hazelbauer GL. Nanodiscs separate chemoreceptor oligomeric states and reveal their signaling properties. Proc Natl Acad Sci. 2006 Aug;103(31):11509–11514.
  18. Maddock, Shapiro L. Polar location of the chemoreceptor complex in the Escherichia coli cell. Science. 1993 Mar;259(5102):1717–1723.
  19. Sourjik V, Berg HC. Localization of components of the chemotaxis machinery of Escherichia coli using fluorescent protein fusions. Mol Microbiol. 2000 Aug;37(4):740–751.
  20. Ames P, Studdert CA, Reiser RH, Parkinson JS. Collaborative signaling by mixed chemoreceptor teams in Escherichia coli. Proc Natl Acad Sci. 2002 May;99(10):7060–7065.
  21. Piasta KN, Falke JJ. Increasing and decreasing the ultrastability of bacterial chemotaxis core signaling complexes by modifying protein–protein contacts. Biochem. 2014 Sep;53(35):5592–5600.
  22. Kentner D, Thiem S, Hildenbeutel M, Sourjik V. Determinants of chemoreceptor cluster formation in Escherichia coli. Mol Microbiol. 2006 Jul;61(2):407–417.
  23. Endres RG, Oleksiuk O, Hansen CH, Meir Y, Sourjik V, Wingreen NS. Variable sizes of Escherichia coli chemoreceptor signaling teams. Mol Syst Biol. 2008 Aug;4:211.
  24. Hansen CH, Sourjik V, Wingreen NS. A dynamic-signaling-team model for chemotaxis receptors in Escherichia coli. Proc Natl Acad Sci. 2010 Oct;107(40):17170–17175.
  25. Peach ML, Hazelbauer GL, Lybrand TP. Modeling the transmembrane domain of bacterial chemoreceptors. Protein Sci. 2002 Apr;11(4):912–923.
  26. Vaknin A, Berg HC. Osmotic stress mechanically perturbs chemoreceptors in Escherichia coli. Proc Natl Acad Sci. 2006 Jan;103(3):592–596.
  27. Ottemann KM, Xiao W, Shin YK, Koshland DE. A piston model for transmembrane signaling of the aspartate receptor. Science. 1999 Sep;285(5434):1751–1754.
  28. Briegel A, Ortega DR, Tocheva EI, Wuichet K, Li Z, Chen S, et al. Universal architecture of bacterial chemoreceptor arrays. Proc Natl Acad Sci. 2009 Oct;106(40):17181–17186.
  29. Briegel A, Li X, Bilwes AM, Hughes KT, Jensen GJ, Crane BR. Bacterial chemoreceptor arrays are hexagonally packed trimers of receptor dimers networked by rings of kinase and coupling proteins. Proc Natl Acad Sci. 2012 Mar;109(10):3766–3771.
  30. Briegel A, Wong ML, Hodges HL, Oikonomou CM, Piasta KN, Harris MJ, et al. New insights into bacterial chemoreceptor array structure and assembly from electron cryotomography. Biochem. 2014 Mar;53(10):1575–1585.
  31. Liu J, Hu B, Morado DR, Jani S, Manson MD, Margolin W. Molecular architecture of chemoreceptor arrays revealed by cryoelectron tomography of Escherichia coli minicells. Proc Natl Acad Sci. 2012 Jun;109(23):E1481–E1488.
  32. Wang X, Vu A, Lee K, Dahlquist FW. CheA–receptor interaction sites in bacterial chemotaxis. J Mol Biol. 2012 Sep;422(2):282–290.
  33. Mello BA, Tu Y. Effects of adaptation in maintaining high sensitivity over a wide range of backgrounds for Escherichia coli chemotaxis. Biophys J. 2007 Apr;92(7):2329–2337.
  34. Shimizu TS, Le Novère N, Levin MD, Beavil AJ, Sutton BJ, Bray D. Molecular model of a lattice of signalling proteins involved in bacterial chemotaxis. Nat Cell Biol. 2000 Nov;2(11):792–796.
  35. Greenfield D, McEvoy AL, Shroff H, Crooks GE, Wingreen NS, Betzig E, et al. Self-organization of the Escherichia coli chemotaxis network imaged with super-resolution light microscopy. PLoS Biol. 2009 Jun;7(6):e1000137.
  36. Endres RG, Falke JJ, Wingreen NS. Chemotaxis receptor complexes: from signaling to assembly. PLoS Comp Biol. 2007 Jul;3(7):e150.
  37. Duke TAJ, Bray D. Heightened sensitivity of a lattice of membrane receptors. Proc Natl Acad Sci. 1999 Aug;96(18):10104–10108.
  38. Vu A, Wang X, Zhou H, Dahlquist FW. The receptor–CheW binding interface in bacterial chemotaxis. J Mol Biol. 2012 Jan;415(4):759–767.
  39. Sourjik V, Vaknin A, Shimizu TS, Berg HC. In vivo measurement by FRET of pathway activity in bacterial chemotaxis. Meth Enzymol. 2007 Feb;423:365–391.
  40. Bray D, Lay S. Computer-based analysis of the binding steps in protein complex formation. Proc Natl Acad Sci. 1997 Dec;94(25):13493–13498.
  41. Monod J, Wyman J, Changeux JP. On the nature of allosteric transitions: a plausible model. J Mol Biol. 1965 May;12(1):88–118.
  42. Vaknin A, Berg HC. Physical responses of bacterial chemoreceptors. J Mol Biol. 2007 Mar;366(5):1416–1423.
  43. Lan G, Schulmeister S, Sourjik V, Tu Y. Adapt locally and act globally: Strategy to maintain high chemoreceptor sensitivity in complex environments. Mol Syst Biol. 2011 Mar;7(1):475.
  44. Khursigara CM, Lan G, Neumann S, Wu X, Ravindran S, Borgnia MJ, et al. Lateral density of receptor arrays in the membrane plane influences sensitivity of the E. coli chemotaxis response. EMBO J. 2011 Mar;30(9):1719–1729.
  45. Zhang P, Khursigara CM, Hartnell LM, Subramaniam S. Direct visualization of Escherichia coli chemotaxis receptor arrays using cryo-electron microscopy. PNAS. 2007 Mar;104(10):3777–3781.
  46. Kollmann M, Løvdok L, Bartholomé K, Timmer J, Sourjik V. Design principles of a bacterial signalling network. Nature. 2005 Nov;438(7067):504–507.
  47. Gebert JF, Overhoff B, Manson MD, Boos W. The Tsr chemosensory transducer of Escherichia coli assembles into the cytoplasmic membrane via a SecA-dependent process. J Biol Chem. 1988 Nov;263(32):16652–16660.
  48. Shiomi D, Yoshimoto M, Homma M, Kawagishi I. Helical distribution of the bacterial chemoreceptor via colocalization with the Sec protein translocation machinery. Mol Microbiol. 2006 Apr;60(4):894–906.
  49. Briegel A, Ames P, Gumbart JC, Oikonomou CM, Parkinson JS, Jensen GJ. The mobility of two kinase domains in the Escherichia coli chemoreceptor array varies with signalling state. Mol Microbiol. 2013 Sep;89(5):831–841.
  50. Scott M, Klumpp S, Mateescu EM, Hwa T. Emergence of robust growth laws from optimal regulation of ribosome synthesis. Mol Syst Biol. 2014 Aug;10(8):747.
  51. Weiße AY, Oyarzún DA, Danos V, Swain PS. Mechanistic links between cellular trade-offs, gene expression, and growth. PNAS. 2015 Mar;112(9):E1038–E1047.
  52. Briegel A, Ladinsky MS, Oikonomou C, Jones CW, Harris MJ, Fowler DJ, et al. Structure of bacterial cytoplasmic chemoreceptor arrays and implications for chemotactic signaling. eLIFE. 2014 Mar;3:e02151.
  53. Cardozo MJ, Massazza DA, Parkinson JS, Studdert CA. Disruption of chemoreceptor signalling arrays by high levels of CheW, the receptor–kinase coupling protein. Mol Microbiol. 2010 Mar;75(5):1171–1181.
  54. Donaldson T, Dattelbaum J. Development of whole cell biosensors mediated by bacteria chemoreceptors. FASEB J. 2014 Apr;28:614.5.
  55. Brent RP. Algorithms for minimization without derivatives. Mineola, N.Y: Dover Pubn Inc; 2003.
  56. Pomerantsev AL. Confidence intervals for nonlinear regression extrapolation. Chemometrics and Intelligent Laboratory Systems. 1999 Sep;49(1):41–48.
  57. Andrae R, Schulze-Hartung T, Melchior P. Dos and dont’s of reduced chi-squared. arXiv. 2010 Dec;1012.3754.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description