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.
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–3]. Despite the importance of receptor clustering in accounting for these signaling properties [4–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–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–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]. 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 [8,34–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 indirectly study the cooperative behavior of the specific receptors inside the cells 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.
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 K D on and K D off for a specific ligand. As the attractant affinity of inactive receptors is higher than for active receptors (K D on ≫ K D off), the state ratio tips towards inactive receptors with increasing ligand concentration c. In contrast, receptor modification m favors the active state in the absence of ligands represented by an energy offset Δϵ (m). The resulting single-dimer free energies in the active and inactive states are given by f on = Δ ϵ (m) - ln (1 + c K D on) + μ f off = - ln (1 + c K D off) + μ, (1) with μ the chemical potential of the receptors in the membrane. All energies are expressed in units of the thermal energy, kB T. In our approach, we allow for an ensemble of different complexes with varying complex size x (i. e. number of connected TDs) and partially developed linkers as rest groups R (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 μW and μA2 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 μ W = ln ( (K D W · K D A) 1 / 2 [ W ]) μ A 2 = ln (K D A [ A ]), (2) where [W] and [A] indicate monomer concentrations and K D W and K D A are dissociation constants for CheW–receptor and CheW–CheA2 binding, respectively. In particular [W] and [A] are expressed as fractional changes i and j of wild-type expression levels [W]0 and [A]0, respectively: [ W ] (i) = i · [ W ] 0 [ A ] (j) = j · [ A ] 0. (3) The TD is assumed to be the smallest receptor unit [9,42], 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 J 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 x and rest group R are given by (cf. Fig 1D) F on (x, R) = 3 x f on + (x - 1) (μ A 2 + 2 μ W + J) + R (μ A 2, μ W) F off (x, R) = 3 x f off + (x - 1) (μ A 2 + 2 μ W) + R (μ A 2, μ W), (4) with 3x receptor dimers per complex of size x and x − 1 linkers. Such a complex has x + 2 rest groups with R (μA2, μW) given by R 1 = 0 R 2 (μ W) = (x + 2) μ W R 3 (μ A 2, μ W) = (x + 2) (μ W + μ A 2), (5) for Eq (1) no rest group, Eq (2) a CheW and Eq (3) a CheW and a CheA dimer, respectively. The probability PS for a certain complex type S (x, R) and its probability P S on of being active follow from standard combinatorial reasoning and the partition function Z Z ≡ 1 + ∑ S (e - F on (S) + e - F off (S) ) (6) P S = e - F on (S) + e - F off (S) Z (7) P S on = (1 + e F on (S) - F off (S) ) - 1, (8) where the number 1 in the partition function Z reflects the possibility of an empty membrane site. Assuming the FRET signal to report the number nA2 (S) of CheA2 dimers within an active complex, we define the receptor activity as A = ∑ S P (S, on) · n A 2 (S) = ∑ S P S · P S on · n A 2 (S). (9) In contrast, the classical MWC model for coupled receptors describes the response of a single complex of N TDs to a change in ligand concentration. Without incorporating the receptor coupling explicitly, the corresponding activity A ˇ reads [23] A ˇ = (1 + exp [ N { Δ ϵ (m) + log (1 + c / K D off 1 + c / K D on) } ]) - 1. (10) In the past, the Hill coefficient nH and complex size N have broadly been treated as equivalent to quantify the cooperative behavior of receptor complexes, and in [23], an increase in N 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 nH = N requires c ⪡ K D on [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 N for all modification levels. This treatment results in a similar quality of fit when relating the reduced number of parameters to the new χ2 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 N in the classical MWC model is naturally generalized by a constant average complex size 〈N〉 with respect to ligand concentration and receptor-modification state. The average complex size, which we term receptor density ρ, is given by ρ = ∑ S 3 · x · P S = 3 ⟨ x ⟩ ≡ constant, (11) with x being the number of dimers of a given complex type S. 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.
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). Fig 2B displays the fitted receptor density ρ next to the Hill coefficients nH of the experimental curves (see Materials and Methods) and the complex size N 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. 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 μ ≡ ∂F/∂N reflects the amount of energy required for adding a particle to a system with free energy F. 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 c < cH, with cH being the half-maximum concentration obtained from Hill fits, the cost necessary to maintain a constant density increases rapidly for ligand concentrations beyond cH. In this second regime, the curves for all modification levels m are of the form f (c) = f0 + ln c, which is the functional description of an ideal chemical potential. Although the slope in the second regime is the same for all values of m, the different offsets f0 (m) reflect the modification-dependent energy Δϵ (m). 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. 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. Fig 3A shows modeled dose-response curves for different ρ values of 1. 5 ⋅ ρ0, ρ0 and 0. 5 ⋅ ρ0 with ρ0 = 7. 5 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. Fig 3B reflects the associated trend in cooperativity by comparing density ρ and Hill coefficient nH. 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]. 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 [W]0 and [A]0, 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 (Δϵ, ρ, K D, Tsr on, K D, Tsr off, J, μ W 0 and μ A 2 0). Multiplication of the calculated activities with scaling parameters sA and sW, 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]. Fig 4A and 4B 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 and 4D). Although the nH values from the model change significantly with expression level [W] at a 95% confidence level, which is in qualitative agreement with the experimental data, especially with respect to the highest CheW expression level, the change in nH 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. 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, 5B and 5C). 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. The opposing trends in nH 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 〈ρ〉 = 3〈x〉 = 9 (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. 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 μW2; we keep the previously introduced rest groups for simplicity. Fig 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. 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 μ W 0, μ A 2 0 and μ W 2 0, 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].
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 μW and μA2 respectively for each CheW and CheA2 molecule within a complex as part of a fully or partially developed linker. A coupling energy J < 0 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 c. 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 N and the curves’ Hill coefficients nH within the classical MWC model. An increase in N is not directly associated with an increase in nH. In our statistical-mechanics model, the approximately constant nH 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 J ^ (m) without modeling the actual protein connectivity. Limited conformational spread and hence a finite complex size is achieved by using a long-range repulsion energy U between all trimers within a complex. In contrast, our model is simpler while providing valuable insights. Neither μW and μA2 nor J 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 μ (c) provides insights into the energetic cost of insertion of receptors into the membrane and its dependence on ligand concentration c, albeit based on an equilibrium mechanism. For constant J 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]). 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] (Figs 8C and 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 [29–31], 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].
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 A (c) Eq (12) with amplitude A0, half-maximum concentration cH and Hill coefficient nH are fitted to the model evaluated at 50 logarithmically spaced concentrations between c = 0. 001mM and c = 1mM. The Hill coefficients in the comparative plot Fig 2B result from direct fitting to the experimental data. A (c) = A 0 1 + (c c H) n H (12) 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 χ2 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.
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.
Protein Connectivity in Chemotaxis Receptor Complexes