Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2874687

Formats

Article sections

Authors

Related links

Phys Rev Lett. Author manuscript; available in PMC 2010 May 22.

Published in final edited form as:

Published online 2009 December 8.

PMCID: PMC2874687

NIHMSID: NIHMS202243

The publisher's final edited version of this article is available at Phys Rev Lett

See other articles in PMC that cite the published article.

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 *p*H 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* [1]. Motor rotation is typically driven by H^{+} ions that generate torque by passing into the cell via the motor, down an electrochemical 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.

(color). Schematic model of the bacterial flagellar motor. Left: The passage of a proton, or possibly a group of protons, through a torque-generating unit (a MotA/B stator—only one stator of about 10 is shown) causes a protein spring to stretch **...**

The flagellar motor operates with near-perfect efficiency at low speeds [2]. 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 [3–6] and modeling progress [6–11], the origin of the knee is still poorly understood. In [11,12], 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 in which the knee arises from the crossover between two time scales, one governing mechanical relaxation and the other proton translocation. In our model, proton translocation, which is assumed to be the rate-limiting step, is modeled by a barrier crossing event as in [11], with the difference that the barrier height depends on the mechanical tension between stator and rotor. 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 2 [13]. (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) [10,14]. This assumption enforces the tight coupling between proton current and rotation speed [2]. (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 nonlinear dependence of the TSR on proton concentrations [Fig. 2]. (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}_{\text{ext}}^{n}/({K}_{\text{ext}}^{n}+{H}_{\text{ext}}^{n})$, 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(

$${J}_{\text{in}}=n{J}_{0}\frac{{H}_{\text{ext}}^{n}}{{K}_{\text{ext}}^{n}+{H}_{\text{ext}}^{n}}\frac{{K}_{\text{int}}^{n}}{{K}_{\text{int}}^{n}+{H}_{\text{int}}^{n}}\text{exp}\left(-n\frac{{U}_{\text{in}}}{{k}_{B}T}\right),$$

(1)

*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

$${J}_{\text{in}}-{J}_{\text{out}}=n{J}_{0}{e}^{-n{U}_{\text{in}}^{0}/{k}_{B}T}\frac{{H}_{\text{ext}}^{n}}{{K}_{\text{ext}}^{n}+{H}_{\text{ext}}^{n}}\frac{{K}_{\text{int}}^{n}}{{K}_{\text{int}}^{n}+{H}_{\text{int}}^{n}}\phantom{\rule{thinmathspace}{0ex}}\times \phantom{\rule{thinmathspace}{0ex}}{e}^{-n(\alpha \mathrm{\Delta}\psi +\beta \tau \delta \theta )/{k}_{B}T}\left[1-\text{exp}\left(n\frac{e\mathrm{\Delta}p+\tau \delta \theta}{{k}_{B}T}\right)\right].$$

(2)

Δ*p* := Δψ + (*k _{B}T/e*) log(

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

$$\frac{d\theta}{\mathit{\text{dt}}}=\frac{1}{\nu}{\displaystyle \sum _{i=1}^{N}{\tau}_{i}+\xi (t),}$$

(3)

$$\frac{d{\tau}_{i}}{\mathit{\text{dt}}}=k({\tau}_{i})\phantom{\rule{thinmathspace}{0ex}}\left[-\frac{d\theta}{\mathit{\text{dt}}}+\mathrm{\Omega}({\tau}_{i})+{\xi}_{i}^{\tau}(t)\right].$$

(4)

ν 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 (backreaction of the rotor onto stators, −*d*θ/*dt*), but gets restretched by proton translocations (Ω(τ_{i}) := (*J*_{in} − *J*_{out})δθ). ξ(*t*) is the thermal noise on the load and satisfies the Einstein relation: ξ(*t*)ξ(*t*′) = 2(*k _{B}T*/ν)δ(

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°, *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.

Figure 2 shows fits of TSRs measured under various *p*H conditions [5]. The electric potential Δψ was not measured and so was used as a fitting parameter for each set of *p*H conditions. Our fit indicates that |Δψ| increases with *p*H (Fig. S1 [17]), consistently with previous measurements [18,19].

Electrorotation experiments [3,7] 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 [7]. When the motor is driven to speeds larger then the maximum operating speed (lower-right quadrant of the TSR), 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 [3]. 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. 2.

Within our model, the proton flux is limited by the loading of the external and internal gates and by the barrier crossing. This limitation on flux accounts for the knee of the TSR. In Fig. 2, the position of the knee strongly depends on *p*H 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 torque.

Our model also accounts quantitatively for measurements of the TSR at different temperatures (Fig. S2 [17]) [4], as well as for measurements with different numbers of stators [Inset of Fig. 2] [6]. We compared the model to earlier measurements at different temperatures [3] (data reported in [4]), with no additional fitting parameters, and found excellent agreement (Fig. S3 [17]).

Experiments by Gabel and Berg (Fig. 3) [20] 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 [20]. 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. 3. These two loads correspond to the two dashed lines in the upper inset of Fig. 3. As cells were deenergized 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 deenergization the electric and chemical parts of the PMF regressed in fixed ratio to each other. The model fits are consistent with the data, both at 24 °C (Fig. S6 [17]) and at 16.2 °C (Fig. 3). However, a further reduction of the low load would be predicted to lead to a strong deviation from proportionality (dotted line in Fig. 3).

(color). Relationship between the high-speed (low-load) and low-speed (high-load) regimes of the *E. coli* flagellar motor, as the PMF varies from −150 mV to 0, at temperature *T* = 16.2 °C. Symbols: Experimental data from individual cells **...**

Adding cooperativity (*n*>1) to the model still captures the general shape of the TSR, as well as its temperature dependence (Fig. S4 [17]), but agrees poorly with the observed *p*H dependence (Fig. S5 [17]), as it implies a stronger dependence of proton current on ion concentrations. The model with *n*>1 is consistent with the Gabel-Berg data for most cells (Fig. S7 [17]), 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 ${\xi}_{i}^{\tau}(t)$ in the stretching of the protein spring at each stator *i*. We can approximate the shot noise as Gaussian white noise: ${\xi}_{i}^{\tau}(t){\xi}_{i}^{\tau}(t\prime )=2{D}_{\text{shot}}\delta (t-t\prime )$, with *D*_{shot} = (1/2)(*J*_{in} + *J*_{out})*n*δθ^{2} (this approximation is valid for times much larger than ${J}_{\text{in}}^{-1}\phantom{\rule{thinmathspace}{0ex}}\text{and}\phantom{\rule{thinmathspace}{0ex}}{J}_{\text{out}}^{-1}$). Solving Eqs. (3) and (4), including shot noise ${\xi}_{i}^{\tau}(t)$ and thermal noise ξ(*t*), we find an exact expression for the effective diffusion coefficient of the rotor angle:

$$\begin{array}{cc}{D}_{\text{eff}}\hfill & :=\underset{t\to \infty}{\text{lim}}\frac{1}{2t}(\theta {(t)}^{2}-\theta {(t)2}^{})\hfill & =\frac{{k}_{B}T}{\nu}{\left(1-\frac{N\mu}{\nu +N\mu}\right)}^{2}+\frac{{D}_{\text{shot}}}{N}{\left(\frac{N\mu}{\nu +N\mu}\right)}^{2},\hfill \hfill \end{array}$$

(5)

where *N* is the number of stators and µ(τ) = −(*d*Ω/*d*τ)^{−1} is minus the local slope of the TSR. Figure 4 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.

(color). Effective diffusion as a function of rotation speed for a single stator (*N* = 1), with different ion translocation cooperativities *n* = 1, … , 4 (bottom to top). Inset: schematic of cooperative ion translocation.

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:

Equation (5) could be used to infer *n* experimentally. In the low-load limit, we have *J*_{in} *J*_{out}, and therefore *n* ≈ 2*D*_{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 [21,22], 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 Avidgor Eldar, Yigal Meir, and Anirvan Sengupta for helpful suggestions. T. M. was supported by the Human Frontier Science Program and N. S.W. by National Institutes of Health Grant No. R01 GM082938.

1. Sowa Y, Berry RM. Q. Rev. Biophys. 2008;41:103. [PubMed]

2. Meister M, Berg HC. Biophys. J. 1987;52:413. [PubMed]

3. Berg HC, Turner L. Biophys. J. 1993;65:2201. [PubMed]

4. Chen X, Berg HC. Biophys. J. 2000;78:1036. [PubMed]

5. Nakamura S, Kami-ike N, Yokota JP, Kudo S, Minamino T, Namba K. J. Mol. Biol. 2009;386:332. [PubMed]

6. Ryu WS, Berry RM, Berg HC. Nature. Vol. 403. London: 2000. p. 444. [PubMed]

7. Berry RM, Berg HC. Biophys. J. 1999;76:580. [PubMed]

8. Iwazawa J, Imae Y, Kobayasi S. Biophys. J. 1993;64:925. [PubMed]

9. Meacci G, Tu Y. Proc. Natl. Acad. Sci. U.S.A. 2009;106:3746. [PubMed]

10. Meister M, Caplan SR, Berg HC. Biophys. J. 1989;55:905. [PubMed]

11. Xing J, Bai F, Berry R, Oster G. Proc. Natl. Acad. Sci. U.S.A. 2006;103:1260. [PubMed]

12. Bai F, Lo C-J, Berry RM, Xing J. Biophys. J. 2009;96:3154. [PubMed]

13. Yuan J, Berg HC. Proc. Natl. Acad. Sci. U.S.A. 2008;105:1182. [PubMed]

14. Mora T, Yu H, Sowa Y, Wingreen NS. PLoS Comp. Biol. 2009;5:e1000540. [PMC free article] [PubMed]

15. Fung DC, Berg HC. Nature. Vol. 375. London: 1995. p. 809. [PubMed]

16. Meister M, Lowe G, Berg HC. Cell. 1987;49:643. [PubMed]

17.
See EPAPS Document No. E-PRLTAO-103-010950 for supplementary figures. For more information on EPAPS, see http://www.aip.org/pubservs/epaps.html

18. Lo C-J, Leake MC, Pilizota T, Berry RM. Biophys. J. 2007;93:294. [PubMed]

19. Minamino T, Imae Y, Oosawa F, Kobayashi Y, Oosawa K. J. Bacteriol. 2003;185:1190. [PMC free article] [PubMed]

20. Gabel CV, Berg HC. Proc. Natl. Acad. Sci. U.S.A. 2003;100:8748. [PubMed]

21. Samuel AD, Berg HC. Proc. Natl. Acad. Sci. U.S.A. 1995;92:3502. [PubMed]

22. Samuel AD, Berg HC. Biophys. J. 1996;71:918. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |