{'fraction_non_alphanumeric': 0.04484132082744443, 'fraction_numerical': 0.013490941796222537, 'mean_word_length': 4.320232399179767, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Modeling torque versus speed, shot noise, and rotational diffusion of the bacterial flagellar motor

We present a minimal physical model for the flagellar motor that enables bacteria to swim. Our model explains the experimentally measured torque-speed relationship of the proton-driven E. coli motor at various pH and temperature conditions. In particular, the dramatic drop of torque at high rotation speeds (the"knee") is shown to arise from saturation of the proton flux. Moreover, we show that shot noise in the proton current dominates the diffusion of motor rotation at low loads. This suggests a new way to probe the discreteness of the energy source, analogous to measurements of charge quantization in superconducting tunnel junctions.

The bacterial flagellar motor is a molecular machine that rotates a helical filament and thereby powers the swimming of bacteria like E. coli [19]. Motor rotation is typically driven by H + ions that generate torque by passing into the cell via the motor, down an electro-chemical gradient called the proton motive force (PMF). Although much work has been devoted to understanding proton translocation and its coupling to torque generation, biochemical details are lacking and many questions remain unanswered. An important one is whether ion translocation is cooperative, i.e. whether protons translocate individually or in groups. Here, we present a minimal physical model for torque generation (Fig. 1) that not only explains a variety of previous experimental observations, but also suggests a way to measure the cooperativity of proton translocation. Specifically, the model predicts that at low loads, motor diffusion is dominated by proton shot noise with a strong (quadratic) dependence on proton cooperativity.
The flagellar motor operates with near-perfect efficiency at low speeds [10]. As the speed is increased, e.g. by reducing the load, the torque and efficiency initially remain high-the "plateau" of the torque-speed relationship (TSR)-and then drop abruptly at a "knee" (cf. Fig. 2). This knee occurs at higher speeds as temperature is increased. Despite much experimental [2,4,15,16] and modeling progress [3,7,9,11,16,20], the origin of the knee is still poorly understood. In [1,20] the cause of the knee was argued to be the gating of proton translocation by the relative position between stator and rotor. In [9], a detailed model of motor kinetics was proposed to explain the observation that motor speed is independent of the number of stators at low loads [21]. In our model, proton translocation, which is assumed to be the rate-limiting step, is modeled by a barrier crossing event. The knee in the TSR then arises from the kinetically limited rate of proton translocation. Importantly, our model fully incorporates proton thermodynamics and yields the separate dependence of the TSR on the electrical and chemical parts of the PMF.
Three ingredients underlie our model: (i) Each torque- generating unit (MotA/B stator) contributes independently and additively to the total torque, in agreement with experimental observations (Fig. 2b) [21]. (ii) The torque from each MotA/B stator is applied to the rotor by a protein spring. Proton translocation into the cell causes the stretching of a protein spring to its next attachment site ( Fig. 1) [11,14]. This assumption enforces the tight coupling between proton current and rotation speed [10]. (iii) Assuming a cooperativity index of n, translocation occurs through three reversible steps: first n protons load into an external gate, then all n cross an energy barrier to an internal gate, and finally all n are released into the cell. The barrier crossing event is the rate-limiting step. The external and internal gates are necessary to explain the non-linear dependence of the arXiv:0910.0835v1 [q-bio.SC] 5 Oct 2009 TSR on proton concentrations (Fig. 2a). (We define the distance between two attachment sites as nδθ, so that the average displacement per proton is δθ.) The external and internal gates are assumed to be in fast equilibrium with the external and internal proton concentrations, respectively. Their dissociation constants are denoted by K ext and K int , so that the occupancy of the external gate is H n ext /(K n ext + H n ext ), where H ext is the proton concentration outside the cell, and similarly for the internal gate. The energy difference between the internal and external gates is n × [e∆ψ + τ δθ + k B T log(K int /K ext )]. ∆ψ is the transmembrane electric potential (ψ int − ψ ext ) and τ δθ is the work necessary to stretch the protein spring. When ∆ψ = 0 and τ = 0, the energy barrier per proton to inward translocation is U 0 in , and the barrier to outward translocation is U 0 out , with U 0 in − U 0 out = k B T log(K int /K ext ). In general, some fraction α of the electric potential and some fraction β of the work contribute to the inward barrier, so that the barriers per proton to inward and outward translocations are, respectively, For simplicity we have neglected any dependence of K int,ext on ∆ψ and τ . Then the rate of inward proton translocations is: J 0 is a kinetic constant (in Hz), and the other prefactors represent the occupancies of the external and internal gates. The outward rate J out is given by a similar expression, so that the net inward proton flux is: ∆p := ∆ψ + (k B T /e) log(H int /H ext ) is the PMF composed of the electrical and chemical potential differences. Its value is approximately −150 mV in normal conditions. To account for the data, we assume that the height of the barrier may depend on temperature, and we expand the prefactor to linear order in temperature: Fig. 2b, is chosen as a reference temperature.
Rotation is then described by coupled stochastic equations for the angular position of the rotor θ and the stretching of the protein springs i = 1, . . . , N , where N is the number of stators, each exerting a torque τ i on the rotor: ν is the frictional drag coefficient of the load, and k(τ ) is the spring constant of the (possibly non-Hookean) protein springs. The spring constant need not be specified as none of the observables computed below depend on it. In the second equation, each spring relaxes as the rotor moves (back-reaction of the rotor onto stators, −dθ/dt), but gets restreched by proton translocations is the thermal noise on the load and satisfies the Einstein relation: . ξ τ i is the shot noise at each stator i, due to the randomness of proton translocation events. To obtain the average speed ω and torque per stator τ = τ i , we solve Eqs. (3,4) at steady state (dτ i /dt = 0) and in the absence of noise, yielding: ω = dθ/dt = N τ /ν and τ i = τ , with Ω(τ ) = ω. Along with Eq. (2), this gives a closed set of equations from which we obtain the TSR. Note that the resulting expression depends separately on the electrical and chemical potential for protons. At stall (ν → ∞) the system is in equilibrium. The energy necessary for a protein spring to move to its next attachment site, which is proportional to the torque τ it exerts on the rotor, is matched by the PMF, τ δθ + e∆p = 0. Consequently the total torque grows linearly with the PMF, N τ = −N e∆p/δθ, and the efficiency near stall is ≈ 100%, in agreement with experiments [5,12].
Our model with no cooperativity (n = 1) can fit all existing measured TSR of the E. coli motor. Some our model's parameters are fixed properties of the motor and thus are fit by single values: (δθ = 4.6 o , K int = 1.2·10 −8 , K ext = 2 · 10 −7 ,J 0 = 670 Hz, α = 0.2, β = 0.078, η = 0.11 K −1 ) while others depend on conditions (T , H int,ext , ∆ψ, number of stators N ), and may or may not have been measured in the experiments. Fig. 2a shows fits of TSRs measured under various pH conditions [15]. The electric potential ∆ψ was not measured and so was used as a fitting parameter for each set of pH conditions. Our fit indicates that |∆ψ| increases with pH (Fig. 3), consistently with previous measurements [8,13].
Electrorotation experiments [2,3] have been used to apply an external torque on the load via an oscillating field. When the motor is driven backwards (upper-left quadrant of the TSR), the internal torque is approximately equal to its stall value up to speeds of −100 Hz [3]. When the motor is driven to speeds larger then the maximum operating speed (lower-right quadrant of theTSR), the motor resists rotation, resulting in a negative internal torque. In this regime, the slope of the TSR remains approximately the same as for positive torques beyond the knee [2]. Our model agrees with measurements in both regimes. The absence of a barrier to backward rotation follows from the reversibility of proton translocation. For negative torques, the model predicts an inflection of the TSR, as seen for the red curve in Fig. 2a.
Within our model, the proton flux is limited by the loading of the external and internal gates and by the  [15]; solid curves: model fits. b. Rotation speed measured with 0.25-0.36 µm diameter beads attached to flagellar stubs for a wide range of viscosities, at three different temperatures [4]; circles: experimental data, solid curves: model fits. Inset: Total torque vs. speed from [16], normalized by the number of stators. Data collapse indicates that stators contribute independently and additively to the total torque. Solid curve: model TSR for a single stator using the same parameters as in the main figure, with the temperature fit as 21 o C and the stall torque fit as 300 pN.nm. barrier crossing. This limitation on flux accounts for the knee of the TSR. In Fig. 2a, the position of the knee strongly depends on pH values: as the internal proton concentration increases, the internal gate gets saturated, preventing protons from translocating inwards, and thus limiting the stretching of stator springs and the applied Electric potential difference versus external pH, as obtained from the fit to data from [15] (Fig. 2a) using ∆ψ = ∆p − (kBT /e) log(Hint/Hext), and assuming that the maximum torque (2500 pN.nm) is reached at ∆p = −150 mV. The apparent linear dependence of ∆ψ on external pH is consistent with previous data [8,13].

torque.
Our model also accounts quantitatively for measurements of the TSR at different temperatures (Fig. 2b) [4], as well as for measurements with different numbers of stators (Inset of Fig. 2b) [16]. As a test, we compared the model to earlier measurements at different temperatures [2] (data reported in [4]), with no additional fitting parameters, and found excellent agreement (Fig. 4).
Experiments by Gabel and Berg (Fig. 5) [6] have been interpreted to imply that the rotation speed is proportional to PMF, even at high speeds beyond the knee of the TSR. Our model predicts that speed is proportional to PMF at low speeds, in the plateau region of the TSR. However at high speeds, the torque is limited by the proton flux, and therefore both torque and speed grow sublinearly with PMF. Nevertheless, our model (with n = 1) is fully consistent with the measurements reported in [6]. Experimentally, speeds were simultaneously recorded for two motors of the same cell, one rotating the cell itself (high load) and the other rotating a small polystyrene bead (low load), as shown schematically in the lower inset of Fig. 5. These two loads correspond to the two dashed lines in the upper inset of Fig. 5. As cells were de-energized by the introduction of a respiratory poison, the PMF ∆p regressed from −150 mV to 0, and the motors slowed, with the two speeds approximately proportional to each other, even at low temperature where the low-load, high-speed motor was in the kinetically-limited regime. We fitted the data for each cell, using N/ν cell and N/ν bead as free parameters, and assuming that during de-energization the electric and chemical parts of the PMF regressed in fixed ratio to each other. The model fits are consistent with the data in the considered parameter regime, both at 24 o C (Fig. 8) and at 16.2 o C (Fig. 5). However, a further reduction of the low load would be predicted to lead to a strong deviation from proportionality (dotted line in Fig. 5).
Adding cooperativity (n > 1) to the model still captures the general shape of the TSR, as well as its temperature dependence (Fig. 6), but agrees poorly with the observed pH dependence (Fig. 7), as it implies a stronger dependence of proton current on ion concentrations. The model with n > 1 is consistent with Gabel and Berg's data for most cells (Fig. 9), although it breaks down for the two cells with the largest high speeds.
Although n = 1 appears to best explain the data, measurements of average rotation speeds do not allow us to discriminate with certainty the proton cooperativity. In contrast, our model predicts that diffusion of the rotor angle at low loads should depend very strongly on proton cooperativity. The discrete nature of proton translocations implies the existence of shot noise ξ τ i (t) in the stretching of the protein spring at each stator i. We can approximate the shot noise as Gaussian white noise: ξ τ i (t)ξ τ i (t ) = 2D shot δ(t − t ), with D shot = (1/2)(J in + J out ) nδθ 2 (this approximation is valid for times much larger than J −1 in and J −1 out ). Solving Eqs. (3,4) including shot noise ξ τ i (t) and thermal noise ξ(t), we find an exact expression for the effective diffusion coefficient of the rotor angle: where N is the number of stators and µ(τ ) = −(dΩ/dτ ) −1 is minus the local slope of the TSR. Fig. 10 shows the effective diffusion coefficient as a function of motor speed for different proton cooperativities. Parameters were chosen so that the speed at zero torque is 200 Hz.
At high loads (ν µ), diffusion is entirely due to thermal noise: D eff ≈ k B T /ν. However, at low loads (ν µ), diffusion is dominated by shot noise: D eff ≈ D shot /N . In fact, the thermal noise is completely suppressed in the low-load limit: e.g., a small thermally induced backward jump in rotor angle causes the stretching of all springs, which then rapidly pull the rotor forward, thus canceling the jump. Notice that in the low-load limit, the shot-  [4]. Right panels: comparison between model prediction and electrorotation data from [2] with the same parameter values as in the left panels. noise contribution is inversely proportional to the number of stators. Intuitively, a small jump in the angular stretch ∆θ i of one protein spring ultimately only causes the rotor to move ∆θ i /N because the rotor is equally coupled to all N stators. The variance per jump is therefore (∆θ i /N ) 2 and with N independent stators, the resulting  diffusion scales as 1/N .
Eq. (5) could be used to infer n experimentally. In the low-load limit, we have J in J out , and therefore n ≈ 2D eff /(ωδθ). δθ can be determined by measuring the torque per stator and the PMF at stall: δθ = −e∆p/τ . Our analysis of rotor diffusion suggests a novel experimental test to investigate the cooperativity of proton translocation. Some rotational diffusion measurements have already been made [17,18], but not in the regime of very low load, where shot noise is expected to dominate. Although we have derived the expression for diffusion in the specific framework of our minimal spring model, the same approach is generalizable to more detailed models of the bacterial flagellar motor.
We thank Avigdor Eldar and Anirvan Sengupta for helpful suggestions. T.M. was supported by a Human Frontier Science Program fellowship, and N.S.W. by National Institutes of Health Grant R01 GM082938.