The reliable response to weak biological signals requires that they be amplified with fidelity. In E. coli, the flagellar motors that control swimming can switch direction in response to very small changes in the concentration of the signaling protein CheY-P, but how this works is not well understood. A recently proposed allosteric model based on cooperative conformational spread in a ring of identical protomers seems promising as it is able to qualitatively reproduce switching, locked state behavior and Hill coefficient values measured for the rotary motor. In this paper we undertook a comprehensive simulation study to analyze the behavior of this model in detail and made predictions on three experimentally observable quantities: switch time distribution, locked state interval distribution, Hill coefficient of the switch response. We parameterized the model using experimental measurements, finding excellent agreement with published data on motor behavior. Analysis of the simulated switching dynamics revealed a mechanism for chemotactic ultrasensitivity, in which cooperativity is indispensable for realizing both coherent switching and effective amplification. These results showed how cells can combine elements of analog and digital control to produce switches that are simultaneously sensitive and reliable.
Bacterial chemotaxis enables the cell to move towards favorable environments. This sensing ability relies closely on collective coordination of several operation modules in the signal transduction pathway (reviewed in [1][2]). The first component of this system is responsible for detecting environmental signals and converting them into intracellular signals. At the surface of the cell, detection of attractants and repellents is mediated by a series of chemoreceptors in the cytoplasmic membrane, the methyl-accepting chemotaxis proteins (MCPs). The second component is the intracellular chemotactic pathway, which processes extracellular signal and converts it into one that is used to determine the behavior of the bacterial flagellar motors: the concentration of the soluble cytoplasmic protein CheY. Binding of repellents induces phosphorylation of CheY, whereas binding of attractants results in CheY dephosphorylation. At the end of the chemotactic pathway lies the final component of the system – the motor block – which changes its switching bias in response to changes in CheY-P (phosphorylated CheY) concentration. On a typical E. coli cell surface, there are 4–5 functioning bacterial flagellar motors. When most of the motors on the membrane spin counterclockwise (CCW), flagellar filaments form a bundle and propel the cell steadily forward; if a few motors (can be as few as one) spin clockwise (CW), flagellar filaments fly apart and the cell tumbles. Therefore the cell repeats a ‘run’-‘tumble’-‘run’ pattern to perform a biased random walk for chemotaxis in a low Reynolds number world [3]. The essential feature of the motor that allows effective chemotaxis is its ability to switch direction quickly and reliably in response to small changes in environmental conditions. Previous studies have revealed that the motor switching responds ultrasensitively to changes in intracellular CheY-P concentration: a high concentration of CheY-P in the cytoplasm of the cell stimulates more CW rotation, while a low concentration of CheY-P results in more CCW rotation. In WT E. coli, the cytoplasmic concentration of CheY-P is around 3 mM and the flagellar motors show stochastic reversals of rotation every second or so [4]. A small change in CheY-P concentration up or down disrupts this equilibrium and produces a large shift toward either CW or CCW rotation. The sensitivity coefficient for the change in rotational bias (time spent in CCW vs. CW) as a function of CheY-P concentration (the Hill coefficient) is ∼10 at the most sensitive part of the region of operation [5]. How the flagellar motor accomplishes this switching behavior is not fully understood, partly because structural data are difficult to obtain. It is known that CheY-P molecules interact with a ring-shaped assembly of about 34 identical FliM protein subunits and this unit is believed to be responsible for determining the direction of rotation [6][7]. For several decades, a series of models have attempted to explain the dynamic behavior of the motor switch and identify the underlying kinetic mechanisms that control the steady state behavior of motor switches [8][9]. Tu and Grinstein [10] used a theoretical argument to suggest that in a dynamical two-state (CW and CCW) model, temporal changes in CheY-P concentration drive the switching behavior of the motor at long time scales and produces a power-law distribution for the durations of the CCW states. Bialek et al. [11] used the bacterial motor as a model system to evaluate the noise limitation of intracellular signaling, concluding that the motor switch operates close to the theoretical limit imposed by diffusive counting noise. A key test for any model of motor switching is the ability to explain how small changes in extracellular concentration are converted into large changes in motor output. To explain this ultrasensitivity, the possibility of cooperative binding of CheY-P to the FliM subunits of the motor switch complex has been suggested [12]. However, studies focused on this binding step [13][14] have determined a Hill coefficient of ∼1 for it, which eliminates the possibility that the amplification is driven by cooperative CheY-P binding to the motor and suggests that a separate, post-binding step within the switch complex is responsible. Duke et al. [15] described a stochastic allosteric model that qualitatively reproduces the ultrasensitive switching and locked state behavior of the motors assuming energetic coupling between neighbor units on the FliM ring inspired by the classic Ising phase transition theory. In particular, this model can reproduce the Hill coefficient of the switch, the nonlinear dependence of rotational bias on CheY-P concentration and the equilibrium between the CW and CCW locked states. The model was based on two assumptions: (a) each subunit of the ring can exist in one of two conformations: CCW and CW state and undergoes a conformational change catalyzed by the binding of CheY-P and (b) a coupling between neighboring subunits favors a coherent configuration and this leads to the propagation of conformational changes along the ring. Although this model is able to qualitatively reproduce the equilibrium behavior of the motor switch, further work is needed to test its ability to reproduce the dynamics of the switching behavior. Here we investigated in detail the behavior of the conformational spread switching model and its ability to reproduce measurements of locked state intervals, the Hill coefficient and other measures of motor dynamics. We then performed a parameter space search to identify the parameters required to best match experimental findings and make further predictions.
Our Monte Carlo model is based on the approach of Duke et al. [15] and our previous work [16], which we briefly describe here. In addition, we make it more general by extending the assumption of symmetry in their original model to include asymmetric cases. The centerpiece of the model is a multi-protein complex or oligomer (to simulate the FliM ring), the individual protomers of which are identical to one another and arranged in a closed ring of size 34. The interface between adjoining ring units represents domains at the boundary between proteins in a biological multi-protein complex. Each protomer can at any time be in either an active (here denoted A and shaded dark in Figure 1a) or inactive (here denoted I and left unshaded in Figure 1a) state, leading to CW and CCW rotation state, respectively. Each protomer can also be bound (here denoted B) or not bound (here denoted N) to a single CheY-P molecule. Then each protomer can make transitions between four possible states, AB↔AN↔IN↔IB↔AB. The model assumes that each protomer can flip reversibly between the two mechanical conformations (CW and CCW) and ligand binding/unbinding changes its chemical conformations, all of which together contribute to a free energy diagram shown in Figure 1a. To reproduce high sensitivity, the model further assumed a coupling energy between the mechanical conformations of adjacent protomers, which favors alike conformations between neighbors, but that the rate constant for CheY-P binding to a protomer is affected only by the conformation of the protomer itself. In the original model of Duke et al. [15], it is assumed that the free energy of the active state, relative to that of the inactive state, changes from +EA to −EA when a protomer binds ligand, for simplicity (Figure 1a). In our model, to make it more general, we introduce two separate energy differences between the active and inactive states: EA0 when a protomer is unliganded and EA1 when it is liganded (Figure 1b). Under the assumption of energy symmetry, the free energy change associated with CheY-P binding can be modeled as, where is the CheY-P concentration required for neutral bias. In the asymmetric case, we use the same definition of, however, does not lead to neutral bias since. In later calculations, we use a numerical method to search for c0. 5 (asymmetric) as a function of. Finally, the model includes a cooperative energy term (EJ, here called cooperativity) so that the free energy of a protomer is lowered by EJ for each neighbor that is in the same state (Figure 1c). This interaction is crucial because it leads to the stochastic creation of semi-stable ‘domains’: regions of the ring whose constituents are either all in the active or all in the inactive state. These domains can then either (a) shrink and disappear, returning the ring to its previous coherent state or (b) grow to encompass the entire ring, a state in which it will remain until another stochastically growing domain of the opposite type will lead to another ring switch. Given all these energy combinations, a protomer can make transitions between all possible states at rate constants describing mechanical conformational changes by: is the sum of the free energy changes associated with changes in activity and interaction. The fundamental flipping frequency, ωa, was set as 104 s−1, a typical rate of protein conformational change and consistent with previous modeling of the switch complex [15]. Lacking information about λa, the parameter that specifies the degree to which changes in the free energy affect forwards as opposed to backward rate constants, an intermediate value of λa = 0. 5 was selected. The free energy associated with CheY-P binding depends only on the conformation of the protomer bound, not on adjacent protomers. So the rate constants describing chemical conformational changes are: where c is the concentration of CheY-P, c0. 5 is the concentration of ligand at which protomers of the ring are 50% occupied on average under the symmetry assumption. ωb is the characteristic binding rate and ΔG (N→B) is the free energy associated with CheY-P binding. A value of ωb = 10 s−1 was selected based on the experimentally determined CheY-P binding rate [17], and consistent with previous modeling of the switch complex [15][16], and λb = 0 such that the binding rate is independent of protomer conformation. In the case of asymmetric EA, the CheY-P concentration corresponding to neutral bias can only be solved numerically. We use custom written C++ code to generate a Monte-Carlo simulation of the conformational spread model [16] (including cases of both symmetric and asymmetric energy). At the beginning of each simulation, each protomer on the ring is set to active and with CheY-P bound. Later on, each protomer n of the ring is assigned two transition times, An and Bn, at which it will undergo a conformational change associated with (A) change between CCW and CW state and (B) CheY-P molecule binds on or off. The program progresses iteratively by locating the event in the [A1, A2……A34, B1, B2……B34] array with the earliest execution time, and after change its state accordingly (either mechanical state or chemical state), new transition time An and Bn of that protomer is updated by t−t0 = −ln (rand) /k, where k is the rate constant for the next transition, t0 is the simulation time when the calculation is made and rand is a random number generated in the interval 0 to 1 [18]. If the transition was associated with a change in mechanical conformation of that protomer, then transition times An+1 and An−1 for the two adjacent protomers are also recalculated (for a closed ring of protomers, we defined An+1 for n = 34 to be A1 and An−1 for n = 1 to be A34). The activity of all protomers on the ring is recorded at integer number MΔt, where Δt = 0. 1 ms as the output sampling time interval and M goes up to 50,000,000. The algorithm continues to update protomer activities on the ring until the simulation time exceeds a specified maximum. Following Duke et al. [15], we assume that switching in the bacterial motor is controlled by the C-ring in the motor complex, which contains 34 copies of the protein FliM and therefore set n = 34 in our model unless otherwise stated. In our model we have in total 4 free parameters: EA0, EA1, EJ, c. The CheY-P concentration c is expressed in the unit of c0. 5 and when we change it we see the ring operate at different bias and therefore the response curve can be plotted. In the following sections, when we make predictions about ring switching time, switching interval etc. , we searched the parameter space EA0, EA1 and EJ across the ranges 0. 5 kBT≤EA0≤1. 5 kBT, 0. 5 kBT≤EA1≤1. 5 kBT, and 3. 5 kBT≤EJ≤4. 5 kBT (shown in Table 1), but for each parameter set, we only present results at neutral bias for simplicity.
A typical screenshot of the ring with multiple domains, labeled with a symbol legend, is shown in Figure 2a. We simulated the qualitative behavior of the ring for a few carefully chosen special cases under the symmetry assumption. Typical screenshots of the ring representing different regimes in the parameter space are shown in Figure 2b. If the activation energy is zero (EA = 0, Figure 2b, top row), growing domains can only form at random through cooperativity between neighbors, but are unstable and unable to grow sufficiently quickly to encompass the ring because any one of their constituent protomers has a high probability of flipping. When the cooperativity energy becomes high, a coherent ring conformation starts to emerge as the coupling between neighboring protomers is sufficiently strong to lock the whole ring in one conformation. However, as the activation energy is zero, the switch complex loses its ability to respond to ligand concentration changes and switching between coherent inactive and coherent active states can be very slow. In the absence of cooperativity (EJ = 0, Figure 2b, middle row), the ring displays random salt-and-pepper patterns reflecting the underlying stochasticity of the ligand binding and unbinding process. When the activation energy becomes high, absolute coupling between chemical conformation and mechanical conformation starts to emerge, and the protomers exist in inactive form only when there is no ligand bound and change to active form once ligand binds. In this case, a coherent active conformation of the ring only exists when there are 34 ligands bound to the ring and for a coherent inactive conformation of the ring we find 0 ligand bound. When cooperativity and activation energy are both present at appropriate magnitudes (EA = 1 kBT, EJ = 4 kBT as discussed in reference [15], Figure 2b, bottom row), the ring spends most of its time locked in either a coherent inactive or active conformation, with transitions between the two (switches) accomplished rapidly by means of a spreading domain. In order to achieve both ring stability and coherent switching, cooperativity is needed to ensure the presence and growth of domains, and activation energy is needed to stabilize these domains by ligand binding. In this energy regime of the model parameter space, the conformational spread model best simulate the performance of the flagellar motor switching responding to external signals. Ring activity, locked states and switching. We simulated the behavior of the 34-protomer ring with the method introduced earlier and EA0 = EA1 = 1 kBT, EJ = 4 kBT. A typical result is shown in Figure 3. The top panel is a time series of the number of active ring protomers (34 active protomers correspond to the CW state and 0 to CCW in our convention). The middle panel shows the number of protomers with bound ligand. This graph matches that of the locked states (i. e. the two fit on top of each other if superimposed) along both axes: ligand binding makes the active state more favorable and the active state binds ligand more strongly and the two effects cooperate to produce locked state and switching behavior. The bottom panel shows the number of independent domains (see Figure 2a for an illustration of typical domain formation present on the ring at any time). Domains appear within a locked CCW or CW state because of stochastic flipping events in protomers and their growth is driven by ligand binding and unbinding and protomer-protomer cooperativity. The domains are transient features of the ring' s behavior. They can either (a) disappear or (b) grow (alone or by fusing with nearby domains) to encompass the whole ring, with these latter events corresponding to switches and occurring very rarely: <1% of domains lead to a switch. We find that at the parameter values identified by Duke et al. [15], the number of independent domains almost never exceed 6 (and that such a ring state only exists ∼0. 0279% of the time). The vast majority of the time, the ring either contains one domain (coherent state, 83. 77% of the time) or contains two domains (15. 34% of the time). Four domains are present 0. 8685% of the time. In addition to being able to reproduce the locked coherent state on the ring and fast switching behavior, a separate key test of the conformational spread model is its ability to reproduce the relationship between changes in CheY-P concentration and motor bias. The Hill coefficient (the maximum sensitivity of the switch) is defined in this case using the relation: where Y is the CW bias, h is the Hill coefficient, c is the concentration of the CheY-P and c0. 5 is the concentration required for neutral bias (in the case of asymmetric model, replace c0. 5 to c0. 5 (asymmetric) ). Here we estimated the Hill coefficient by fitting a linear equation to a plot of log[Y/ (1−Y) ] against log (c/c0. 5), with the slope of this line corresponding to h. For each parameter set, its characteristic Hill curve can be generated by long time simulation with varying c, and plot CW bias of the simulated trace as a function of c. The sensitivity of ring activity to changes in ligand concentration depends more strongly on the activation energy and considerably less on the cooperativity. A lower sensitivity can be brought about by a lower activation energy or by a lower cooperativity, with the activation energy having the dominating influence. However, with a 34 protomer ring, the cooperativity can be no less than the critical cooperativity required for coherent switching to occur, i. e. EJ>3. 5 kBT [15]. In table 1, we present the Hill coefficient calculated for parameters EA0, EA1 and EJ across the ranges 0. 5 kBT≤EA0≤1. 5 kBT, 0. 5 kBT≤EA1≤1. 5 kBT, and 3. 5 kBT≤EJ≤4. 5 kBT. We see that the experimentally determined Hill coefficient ∼10 can be reproduced by a large parameter sets. The behavioral features of the ring can be further characterized by the distributions of (a) the times spent in the two locked states and (b) the times required for both CCW→CW and CW→CCW switches. Here we used simulated ring state data to obtain the theoretical length of the locked state intervals predicted by the model. Because we have direct access to the fundamental protomer states, filtering and threshold algorithms are not needed to identify the intervals (and switches, respectively). Because transitions between the two locked states are not instantaneous, we needed an unambiguous way to define CCW and CW intervals, respectively. We defined such an interval as the time (in simulation steps) between when the ring enters a fully locked state (0 or 34 active protomers, respectively) and when it next enters the other fully locked state (i. e. 34 or 0 active protomers, respectively). Distributions of locked state intervals obtained from simulation traces (EA0 = EA1 = 1 kBT and EJ = 4 kBT at neutral bias) equivalent to 30000 seconds of real time are shown in Figure 4 (a). To make a comparison, the log-linear plot of the distributions at low (0. 2) and high (0. 8) CW biases are also shown in Figure 4 (b). Least-squares fitting of exponential curves to the simulation data are shown overlaid. We see that the locked state distribution follows an exponential distribution. In table 1, we presented the mean locked state interval values calculated for parameter EA0, EA1 and EJ across the ranges 0. 5 kBT≤EA0≤1. 5 kBT, 0. 5 kBT≤EA1≤1. 5 kBT, and 3. 5 kBT≤EJ≤4. 5 kBT. Within the parameter range of our simulation, the minimum value of mean locked state time is 0. 13 s and the maximum value is 22. 17 s (shown in Table 1). The mean locked state time increases when the energy of activation or cooperativity is increased, with the activation energy EA having the dominant influence. The essential feature of interest of the model proposed by Duke et al. [15] is that the ring can simultaneously achieve very rapid switches and very stable locked states. This qualitatively matches what is observed in the rotary motors of flagellar bacteria such as E. coli, which can rotate at hundreds of RPM stably for a long period but switch direction quickly (on the order of ms) and stochastically. Distinct from the classic MWC model, which requires coherent switches to happen instantaneously, in our model switches occur by a mechanism of conformational spread. We defined a switch time as the time (in simulation steps) between when the ring leaves a fully locked state (0 or 34 active protomers, respectively) and when it next enters the other fully locked state (i. e. 34 or 0 active protomers, respectively). We simulated the behavior of the ring at the optimal activation energy and cooperativity values identified earlier (EA0 = EA1 = 1 kBT, EJ = 4 kBT) for different values of bias. The empirical distributions thus determined are shown in Figure 5. In contrast to the distributions of locked state intervals, the switch times follow a peaked gamma distribution. At the parameter value chosen, the mean lies between 58–61 ms for low (0. 2), middle (0. 5) and high (0. 8) biases and these are statistically independent of bias and of direction of switch. In table 1, we presented the mean switch time values calculated for parameter EA0, EA1 and EJ across the ranges 0. 5 kBT≤EA0≤1. 5 kBT, 0. 5 kBT≤EA1≤1. 5 kBT, and 3. 5 kBT≤EJ≤4. 5 kBT. The result of our simulation shows that the mean switch time values changes across the ranges from 12. 03 ms to 322. 65 ms. The mean switch time increases when activation or cooperativity energy is increased, with the activation energy having the dominant influence. To confirm that typically one domain of opposite conformation (rather than several) grows to encompass the entire ring, we also computed power spectra for the ring activity traces in order to characterize the spectral properties of the ring switch complex. If switching events are associated with a single nucleation event (a Possion step), we expect the power spectra of the trace to be monotonically decreasing with a ‘knee’, i. e. display a Lorenzian profile. In contrast, if switching events are associated with multiple hidden steps, as for example in a closed biochemical system with hidden reactions, then we expect a non-Lorentzian profile with a peak (a local maximum) at a characteristic frequency related to the number of steps involved [19]. Our simulation results (Figure 6) show the spectra thus obtained are Lorentzian without a local maximum at long times. This behavior is observed at different values of bias. These results offer an internal confirmation of the model results shown in Figure 4, which indicate that the distributions of times spent in the locked states are exponential. Such a system would be expected to display power spectra with Lorentzian profiles. However, because the power spectra and locked state time distributions are computed independently and by different methods, the result that they predict the same behavior is an important internal test of the model. In particular, the power spectra results confirm that the locked state time distributions are not an artifact of our algorithm for detecting the start and end of a locked state. In our recent experimental paper [16], we used a high-resolution optical system to measure the switching time and locked state interval of bacterial flagellar motors. The experimental observations confirmed that the switching time distribution follows a broad gamma distribution with mean switch time 18. 72 ms and the locked state interval follows an exponential distribution with mean interval value 0. 75 s at neutral bias. Hence we use mean switch time (∼18 ms), mean locked state time (∼0. 75 s) and Hill coefficient (∼10) to parameterize our model. Table 1 shows a coarse parameter search of our model with predictions of the switching time, locked state interval, and Hill coefficient. We identify the 0. 55 kBT≤EA0≤0. 95 kBT, 0. 55 kBT≤EA1≤0. 95 kBT, and 4. 05 kBT≤EJ≤4. 25 kBT region for a fine parameter search (Table 2). We see with only a few parameter sets, the conformational spread model is able to reproduce the three experimentally determined quantities. With the results shown in Table 2, we find 3 groups of values that fit the experimental values determined by Bai et al. [16]. They are EA0 = 0. 55kBT EA1 = 0. 75kBT EJ = 4. 15 kBT, EA0 = 0. 75kBT EA1 = 0. 55kBT EJ = 4. 15 kBT, EA0 = EA1 = 0. 65kBT EJ = 4. 15 kBT. We therefore expect a conformational spread model with activation energy ∼0. 65 kBT and coupling energy ∼4. 15 kBT can well reproduce experimental observations. Please see Figure 7 for a visual summary of our computational results. For simplicity, we only showed those values with EA0 = EA1 = EA and the best-fit parameter set has been labeled by a square. Although we have identified a best-fit parameter set that can well reproduce experimental findings, we have to note that these fit values are sensitive to the parameters fixed earlier, especially to the fundamental flipping frequency ωa. Here we investigate how mean locked state time and mean switch time respond to changes of ωa while other parameters remain fixed. We see in Figure 8 that the fundamental flipping frequency is a scaling factor of the system. Both mean locked state time and mean switch time are inversely scaled by the flipping frequency. When flipping frequency is higher, each protomer on the ring makes more attempts to flip to the opposite conformation and therefore the locked state becomes less stable (hence mean locked state time decreases) and transition becomes much faster (hence mean switch time decreases); when flipping frequency is lower, each protomer on the ring makes fewer attempts to flip to the opposite conformation and therefore locked state becomes more stable (hence mean locked state time increases) and transition becomes much shorter (hence mean switch time increases). In the above sections, we have determined the best-fit model parameters using experimental results. It will be interesting to test the ring behavior at different sizes using those values. When Duke et al. [15] first proposed the conformational spread model, they identified that EJ>kBT ln N (N is the size of the ring) is the condition under which a large ring has the characteristic of a coherent switch. In the case of 34 protomers, this condition requires that EJ>3. 5kBT. When this condition is met, in time series of ring activity, we see for the majority of time that, the ring stays in complete active (active protomer = 34) or complete inactive (active protomer = 0) state. This invokes an empirical mathematical definition of ‘coherent switch’: the active number of protomers on the ring has to be in 0 or N for greater than 65% of the total simulation time. We then simulated the ring activity at sizes of 10,60,100 protomers with activation energy 0. 65 kBT and coupling energy 4. 15 kBT at neutral bias and the result is shown in Figure 9. From the requirement of EJ>kBT ln N we expect to see coherent switch behavior for ring sizes at 10,34, and 60, but not at 100. Indeed, in Figure 9, we see with the same parameter set, the smaller the ring is, the easier a switch happens. At ring size of 10,34,60 protomers, we see clear locked states in the trace and the switching events are fast. However, at the ring size of 100, a switch across the ring becomes very difficult and the time spent during a switch is comparable to the time the ring stays in a locked state. In Table 3, we made predictions about the mean locked state intervals, mean switch times and Hill coefficient of the switch response with activation energy 0. 65 kBT and coupling energy 4. 15 kBT at different ring sizes.
In this study, we undertook a comprehensive numerical simulation analysis of a general model of stochastic allostery in a protein ring and evaluated the ability of such a model to explain the switching, sensitivity and locked state behavior of the rotary bacterial motor. We modeled the gearbox of the motor as a ring of 34 identical protomers, a geometry inspired by the FliM structure in the motor complex, believed to be responsible for motor switching. The model is able to qualitatively reproduce the motor behavior, such as locked rotation in CCW or CW state and fast switching between the two. Furthermore, based on a comprehensive parameter space search, the model can also quantitatively account for the experimentally determined switch time, locked state interval and Hill coefficient of the motor. Specifically, we found a unique set of values that fit the experimental value best, activation energy must be around 0. 65kBT and the cooperativity around 4. 15kBT. The bounds around these values are tight. Smaller or larger energies result in rings that either (a) spend too long or too little time in the locked states, (b) do not have the required sensitivity or (c) far away from this parameter regime, fail to switch coherently. With the ring operating in this parameter set, time traces of ring state (measured as the number of active protomers) indicate that the ring spends most of its time in one of the two locked states, with rare (every 0. 5–2 seconds) switches between the two being accomplished very rapidly (on the order of milliseconds). The trace of ligand activity (measured as the number of ring protomers having bound ligand molecules) mirrors the ring state, with the two driving each other: binding of more ligand drives active domain formation, which in turn leads to a preference for having ligand bound and conversely. Rather than being completely locked in one stable state with all protomers being either active or inactive, the ring displays constant activity in the form of nascent domains of the opposite state to the locked state, seen as ‘noise’. For the vast majority of the time, only one such domain exists, and the presence of two (but no more) growing domains is frequently, but not always, associated with a switching event. The model predicts that the time spent by the ring in the locked states corresponding to CW (all protomers active) and CCW (all protomers inactive) is exponentially distributed. The model can also predict the Hill coefficient (∼10) measured for the sigmoidal curve that relates CheY-P concentration to motor bias. Near the optimal parameter point identified, the distributions of the switching times are gamma-like with a peak around 5–8 ms. To be effective, a switch must achieve two globally conflicting properties. It must accomplish sensitivity by amplifying small changes in the effector, but only over a narrow critical range (the switching point). Outside this range, it must accomplish reliability by being unresponsive to changes in the effector. The allosteric switching model explains how the motor simultaneously meets these competing design requirements. Near the critical CheY-P concentration, a highly cooperative mechanism (EJ≫1kBT) is used to amplify small stochastically occurring nascent domains that can rapidly grow to encompass the entire switching complex. The resulting digital switch displays the desired selective ultrasensitivity but switches chaotically. In order to ensure switch reliability, CheY-P binding must moderately stabilize the protomer active state, providing a mechanism for biasing the whole switch complex by continuously varying the CheY concentration over a large range: a strategy typical of analog control. The values and ratio of the strengths of the two mechanisms must be tightly controlled in order for the switch complex to be functional. We hypothesize that this control is accomplished through the biochemical structure of the protomers and ring, which are genetically determined and so robust to intracellular noise during the cell' s life. In light of recent studies of digital cellular signaling [20], we wish to further suggest that the combination of analog and digital control here proposed to explain the behavior of the bacterial switch complex may be a motif typical of biological switch design. In this study we have focused on a ring consisting of 34 protomers because it is believed that the C-ring in the E. coli flagellar motor, which consists of 34 copies of the protein FliM, acts as the motor direction switch. However, numerous examples of protein rings and other interconnected protein complex geometries are known, including DNA polymerase sliding clamps, voltage-gated ion channels, ATP synthase etc. Each of these rings may hypothetically accomplish its function using a conformational spread mechanism, but would consist of different numbers of protomers. In our model this can be simulated by changing N, the number of elements in the ring. In this paper, we have narrowed our study to a closed protein ring. However, we have to point out that the conformational spread model as well as the numerical method we presented here can be easily modified to describe one dimensional allostery regulation in a protein chain or a strand of DNA molecules. The model can also be modified to describe signal transduction and amplification on a two dimensional plane, which will be of great use in studying functions of cellular receptors. Allostery is a widespread mechanism in biology and conformational change is the basis for a large subset of all protein function. Since protein complexes are the workhorses of the cell, we expect models similar to this and the idea of conformational spread in general to be increasingly important in systems biology and biophysics. Investigating the applicability of conformational spread models to other biological systems will be the subject of future work. Bacterial chemotactic exploration depends on the ability of the flagellar motors at the base of the flagella to perform two tasks: (1) remain stable in their current direction of rotation for long periods (seconds) as required and (2) switch quickly between the two directions in response to the environmental changes detected by the chemotaxis pathway. These properties make the bacterial switch an exquisite computational element that combines ultrasensitivity and reliability. In this paper we presented an analysis of a model featuring conformational spread that aims to explain the mechanism of the motor switch. Simulations confirm that this model is able to reproduce the characteristics of the motor observed in experiments. We speculated that stochastic models of conformational spread will be a common theme in protein allostery and signal transduction.
Bacteria swim to find nutrients or to avoid toxins. Their swimming is powered by the rotation of flagella (hair-like structures) that act as propellers. Each flagellum is driven by a rotary molecular engine (the bacterial flagellar motor) that can rotate in either a counterclockwise or clockwise direction and switches between the two directions are frequent and rapid. Although the motor has been studied in detail, we do not understand how it is able to reliably switch direction – a critical function that gives bacteria the ability to steer. In this paper we examined a mathematical model describing how a potential gearbox in the motor might work inside a ring of identical proteins. We compared the output of this model with experimental data on switching speed and other measures of motor function, finding excellent agreement. This is an exciting finding not only because the operation of the motor itself is important, but also because protein complexes play an important and ubiquitous role in cellular signal transduction and therefore, “conformational spread” may be a widespread mechanism for signal propagation in biology.
theoretical biology
biology
computational biology
biophysics
Conformational Spread in the Flagellar Motor Switch: A Model Study