F1 rotary motor of ATP synthase is driven by the torsionally-asymmetric drive shaft

F1F0 ATP synthase (ATPase) either facilitates the synthesis of ATP in a process driven by the proton moving force (pmf), or uses the energy from ATP hydrolysis to pump protons against the concentration gradient across the membrane. ATPase is composed of two rotary motors, F0 and F1, which compete for control of their shared γ -shaft. We present a self-consistent physical model of F1 motor as a simplified two-state Brownian ratchet using the asymmetry of torsional elastic energy of the coiled-coil γ -shaft. This stochastic model unifies the physical concepts of linear and rotary motors, and explains the stepped unidirectional rotary motion. Substituting the model parameters, all independently known from recent experiments, our model quantitatively reproduces the ATPase operation, e.g. the ‘no-load’ angular velocity is ca. 400 rad/s anticlockwise at 4 mM ATP. Increasing the pmf torque exerted by F0 can slow, stop and overcome the torque generated by F1, switching from ATP hydrolysis to synthesis at a very low value of ‘stall torque’. We discuss the motor efficiency, which is very low if calculated from the useful mechanical work it produces - but is quite high when the ‘useful outcome’ is measured in the number of H+ pushed against the chemical gradient.

1 F 0 ATP synthase (ATPase) either facilitates the synthesis of ATP in a process driven by the proton moving force (pmf), or uses the energy from ATP hydrolysis to pump protons against the concentration gradient across the membrane. ATPase is composed of two rotary motors, F 0 and F 1 , which compete for control of their shared γ -shaft. We present a self-consistent physical model of F 1 motor as a simplified two-state Brownian ratchet using the asymmetry of torsional elastic energy of the coiled-coil γ -shaft. This stochastic model unifies the physical concepts of linear and rotary motors, and explains the stepped unidirectional rotary motion. Substituting the model parameters, all independently known from recent experiments, our model quantitatively reproduces the ATPase operation, e.g. the 'no-load' angular velocity is ca. 400 rad/s anticlockwise at 4 mM ATP. Increasing the pmf torque exerted by F 0 can slow, stop and overcome the torque generated by F 1 , switching from ATP hydrolysis to synthesis at a very low value of 'stall torque'. We discuss the motor efficiency, which is very low if calculated from the useful mechanical work it produces -but is quite high when the 'useful outcome' is measured in the number of H + pushed against the chemical gradient.
Adenosine triphosphate (ATP), the universal fuel of the cell, is synthesized in the membranes of mitochondria, chloroplasts and bacteria by the ATP synthase complex. Already in its abbreviation, F 0 F 1 -ATPase, one can see its dual functionality 1 , which arises from the two competing rotary motors (F 0 and F 1 ) that sit on the opposite ends of the shared drive shaft 2,3 . In one of its operation modes, this multi-subunit enzyme uses free energy stored in a transmembrane electrochemical gradient (pmf) to drive the synthesis of ATP at the site of the F 0 motor. In the other mode, it reverses the reaction, hydrolyzing ATP and utilizing the chemical energy to drive the F 1 motor and pump protons against their concentration gradient 4 . Which of the two competing motors "wins", and consequently in which mode the ATPase complex will operate, is determined by the ratio of the ATP/ADP available to the F 1 part, and the torque exerted by the pmf-driven F 0 . It has been demonstrated that when a clockwise rotation is induced artificially by magnetic tweezers, instead of the F 0 natural action, the ATP synthesis is similarly promoted on the F 1 side of the complex 5 .
In engineering, if two rotary motors share the same drive shaft and are not attached to any support, the outcome depends on their relative mass (inertia moments); this would be impossible in a heavily overdamped molecular system. Thankfully, in its natural setting, ATPase is anchored in the membrane by the a-subunit of its stator; this constraint is then passed through the b-and δ-subunits to immobilize the α 3 β 3 hexamer complex of the F 1 motor. So the only relevant moving parts in the complex are the central γ-shaft, which is made of a coiled-coil protein α-helix, and the disc-like c-subunit, which is rigidly attached to it and which can rotate relative to the a-subunit anchored in the middle of membrane bilayer plane, see Fig. 1(a). The other end of the γ-shaft is inserted into the central cavity of the α 3 β 3 complex, with three identical equilibrium positions separated by 120° representing the catalytic sites on the β-subunits. The "cam" end of the γ-shaft: the off-axis eccentric bulge region of the γ-subunit [6][7][8] , is able to dock into this cavity in three equivalent orientations, when it fits in one of the indentations, Fig. 1(b). These are the distilled structural features of ATPase that we shall need to develop our physical model of the Brownian ratchet F 1 motor; these structural details are extensively characterized in many famous papers 1,3,9,10 .
The mechanism of the F 0 motor is relatively well understood, in both its physical principles and quantitative predictions, following the original work of G. Oster et al. 11 , where it is shown how an asymmetry of H + ion channels in the a-subunit can selectively protonate the Asp61 site on the c-rotor and thus constrain its rotational Brownian motion preferentially in the clockwise direction as seen from the membrane side (assuming the natural

State of the Art
Of the large literature on the ATP-driven rotary motors such as F 1 -ATPase, the most recent and most useful overview is found in papers by Okazaki and Hummer 24 and by Mukherjee and Warshel 25 . In the first, the detailed review of key experimental observations about the motor is assembled and the atomic-level resolution of the α 3 β 3 complex -γ-shaft interaction is presented. The second paper 25 gives a careful comparison of various models of torque generation. The paper itself extends the idea of chemical-rotational free energy landscape that should generate the torque on one 120° leg of the rotation cycle.
In relation to the ideas of torque generation, we would like to point out that in the case of linear motors it is now well established that it is not the chemical-mechanical coupling that produces what people often call the "power stroke". Instead, the whole process is diffusive (not dynamic: instead driven by the Brownian motion), and the switching between conformational states maintains the out-of-equilibrium conditions that allow unidirectional transport 26 . This last point needs to be emphasized since it is a widely spread issue in the literature. Normally, in a dynamical system, if there is a friction drag opposing the motion (rotational friction in our case, expressed via the coefficient of friction and angular velocity: − ξω), then there has to be an external force applied (torque or moment of force in our case: M) in order to sustain the constant-speed motion in what is called the 'terminal velocity' regime: M = ξω. Then the power loss against dissipative friction is the product P = Mω, and the source of this energy is coming from the agent that delivers the torque. The situation is different on the molecular scale, when the physics is dominated by the Brownian motion, i.e. the thermal stochastic force acting on all elements of the system. Later in the text we will derive the expression for the average angular velocity of rotation of the motor, which is strongly affected by the rotational friction via the diffusion constant D = k B T/ξ. In the overdamped regime there is no momentum transfer! (For the same reason, there can be no 'turbine' action as well, even though it is attractive to relate a linear flow of ions (e.g. protons as in F 0 or flagellar motor) to a rotary shaft motion by a turbine association (see: Boyer and Kohlbrenner 1981, and much subsequent literature). In a turbine the rotary motion is induced by the lateral force (leading to torque on the axle) arising from the momentum transfer of the flow through the inclined blade: a process impossible on the molecular scale in a dense fluid environment).
It is incorrect to assume that the motor generates a torque merely because a constant average velocity is observed. The motion occurs due to the thermal noise, while the ATP hydrolysis is needed to take the system out of the detailed-balance equilibrium conditions. In order to generate useful work there has to be a counter-torque against which the stochastic motor has to push, as it is carefully studied by Prost et al. 27 . Later in the text we will estimate the efficiency, which depends on conditions (e.g. becomes zero at the stall torque) but is certainly much less than 100% in all situations. The efficiency is also zero when there is no external torque (free rotation of the motor against friction) since there is no work done (we already pointed out that it is incorrect to interpret the dissipative loss as 'useful work'). Recently, there were several remarkable experiments when an external mechanical torque was applied to the shaft of F 1 motor, generated by magnetic or electric field on a probe particle 5,23,28 , and the balance of energy was calculated. However, again the claim of 100% efficiency was made. We shall discuss these issues and calculation in greater detail at the end of the paper. Figure 2(a) illustrates some details of the above argument. Here two typical traces of the rotor angle are sketched, one at higher the other at lower [ATP], which are a sequence of consecutive 120° steps. This allows to determine the average motor velocity 〈 ω〉 as the average slope of such a trace. Each step is assumed to be initiated by the ATP binding, so the time interval between steps is the inverse reaction rate (discused at length below in the main text). In the past, many papers related this average velocity to the friction drag: M = ξ〈 ω〉 21 , which is clearly in error. The simplest way to see this is by asking why would 〈 ω〉 itself depend on the friction coefficient (as prominently demonstrated in many classical papers 21,33 ), when it is evidently the ratio of 120° and the time interval between the steps (approximately, for the sake of simplicity ignoring the short duration of the step itself for the moment). Curiously, these papers, and many that followed, have reported the torque ξ〈 ω〉 = 40 pN.nm (asserting that this is also what the motor exerts). More recently, a very careful and detailed paper by the group of H. Noji have studied the 'power stroke' step (inset in the figure) and concluded that it is here that the torque M is applied, working against friction to produce a step velocity ω ps 29 . Curiously, this torque was also reported as 40 pN.nm, and discussed broadly in subsequent experimental and theoretical papers 25,30 . If this was such a balance of torque vs. friction, then the (very short) duration of this 'power stroke' would strongly vary with the changing friction coefficient (given that the magnitude of the step remains constant): this is not what observed in experiment.
At this stage it is appropriate to make another comment. A number of papers cited above, both experimental and theoretical, work on the assumption of a "power stroke" produced by the asymmetric structural features of the α 3 β 3 cavity and its interaction with the base of the γ-shaft. But there are many detailed experimental traces of The essential features of experimentally measured rotation of the ATPase rotor, at higher and lower ATP concentration. The sequence of steps of 120° are apparent and the average rotation velocity 〈 ω〉 is determined by the average slope. The inset shows a zoomed-in region of the step, where the apparent slope defines another angular velocity, ω ps . (b) The essence of the stochastic motor mechanism described in this paper. Most of the ATP-binding events end up back in the same rotational step; the occasional jumps forward are determined by the asymmetry of integral probabilities (p + − p − ) of the process. In both cases, the 'active period' duration is τ 2 , outside of which the γ-shaft is confined and its angular fluctuations are bound. rotational angle vs. time (such as the ones sketched in Fig. 2 for illustration) that show an occasional backward step at low [ATP]. The presence of such steps is an evidence that the base of the γ-shaftis actually able to rotate both ways, and this cannot happen if there was a power stroke based on definite structural features.
The other panel of Fig. 2, plot (b), illustrates a scheme of the stochastic motor described in this paper. There are two characteristic times: the inverse rate of ATP hydrolysis and the life-time of the isomerized state τ 2 (which is the duration of the γ-shaft release from its confinement). Outside of this 'active period' the base of the γ-shaft is tightly locked in the α 3 β 3 cavity and the rotational fluctuation of c-ring is controlled by the torsional modulus (measured by Junge et al., see the next section for detail). Although most of the hydrolysis 'events' end up with the rotor in the same position, there is a statistical imbalance towards anticlockwise steps, determined by the torsional asymmetry of the coiled coil γ-shaft and the rotational diffusion constant D = k B T/ξ: that is how our average velocity 〈 ω〉 becomes a natural function of friction ξ. In practice, the experiment is unable to resolve the fast thermal fluctuations, giving only the coarse-grained trace with correct average features.
Among the large amount of experimental work, recently there have been two very important studies of the F 1 motor, examining the role of the γ-shaft, both extremely delicate and sophisticated in their techniques: the group of K. Kinosita 31 have sequentially chopped off the length of the γ-shaft and monitored the change in the motor action, in each case attaching a marker to the shorter and shorter shaft. The group of H. Noji 32 have managed to remove the γ-shaft altogether and monitored the unidirectional rotation of the hydrolysis site on the α 3 β 3 complex by imaging in a fast AFM. The latter paper claims disagreement with the former (and many others), pointing at the apparent unidirectional motor action without any rotor -a remarkable observation in itself. Yet, on careful analysis, we conclude that there is no discrepancy and the results are in fact consistent. It was found that by making the γ-shaft sequentially shorter, the average speed of the motor has become consistently lower (in spite of the friction resistance of the rotating marker being kept the same): changing from well over 100 rps (revolutions per second) for the wild type down to ~0.8 rps for the shortest mutant at a high ATP concentration of 2 mM 31 . The experiment with γ-shaft rotor removed did not have any friction, and recorded an average unidirectional sequence of binding on β subunits with a speed ~0.4 rps at 2 μM ATP. This needs to be compared with the speed of 2-5 rps in the wild-type F 1 motor observed at the same 2 μM ATP with an actin filament creating high friction drag 21 or ~20 rps with a low-drag gold bead 33 . It is clear from both experiments on γ-shaft reduction that there is a complicated (and not known) machinery in the α 3 β 3 complex that maintains the transfer of the hydrolysis site, slowly but consistently in counterclockwise direction. But it is equally clear that the γ-shaft plays a crucial role in achieving the required speeds and the motor efficiency, and its reduction dramatically reduces both. We also note that a coiled-coil filament is a definitive common feature of all other rotary motors driven by ATP hydrolysis (in contrast to pmf-driven motors like F 0 or flagellar).
In this paper we develop a concept of such a rotary motor radically different from what was considered previously: a model that is based on the very generic asymmetric torsional elasticity of the coiled coil γ-shaft. We demonstrate that the motor mechanism is robust and generic, not much influenced by the details of molecular structure (except the 120° periodicity and the parameters of the γ-shaft, which are both crucial). Our idea of a two-state Brownian ratchet motor is not original (only its rotary aspect is): a long time ago it was developed by Jülicher, Ajdari and Prost, and is now an accepted standard for linear motors such as myosin on actin, or kinesin on microtubules 27,34 . For such a motor, at the most basic level one needs just two elements: a periodic profile of the ground-state potential energy with left-right asymmetry (which is the source of symmetry breaking for unidirectional motion), and the external energy input to disturb the equilibrium detailed balance and bring the mechanical system into the upper short-lived 'excited' state (which is what the ATP hydrolysis provides in all cases). For F 1 ATPase the required asymmetric potential energy is the torsional elastic energy of twisting of the coiled-coil γ-shaft about its axis, Figs 1 and 3 below. The linear torsional elasticity of the γ-subunit has been studied in great detail by monitoring the thermal fluctuations of the angle of the c-ring freely rotating with respect to the α 3 β 3 complex immobilized on a substrate. However, the fact that it costs more energy to under-twist the coiled coil, than to over-twist it, is not a priori obvious. Here we model the α-helical coiled coil and produce the distinctly asymmetric torsional potential energy E(θ). The 3-fold, 120° periodicity is inherent in the α 3 β 3 structure, and we plot this potential energy against the angle θ with respect to the neutral equilibrium position of its cam docked in a given slot, finding an asymmetry as illustrated in Figs 4(b) and 5 below. For the model of torsional elasticity of the coiled coil to match the geometry and the measured torsional modulus of the γ-shaft 18,35,36 , there is very little freedom in its choice of parameters. Note that this physical reason for the rotation symmetry breaking of the F 1 -ATPase motor is different from the asymmetry suggested by Oster et al. 37 , where the 'culprit' was proposed to be the shape of the β-subunit confining the bulge-end of the γ-shaft. We do not argue against such an asymmetry of the 'lock' shape (in fact, experiments of Noji et al. 32 on the unidirectional binding with no rotor surely prove its validity), merely point out that it cannot produce a power stroke in overdamped conditions where no momentum transfer can occur.
The next section gives a brief description of the asymmetric torsion of the coiled coil. We then show how the rotary motion is induced in such an asymmetric periodic potential, when a second 'excited' state (corresponding to the undocked cam and the unrestricted rotational diffusion of the c-ring) can be reached by an ATP energy influx. The final section demonstrates how such a rotary motor would operate against a counter-torque arising from the pmf-driven F 0 motor, and find the stall conditions. At the end we compare theoretical predictions with experimental measurement data on ATP synthase F 1 motor action 33,38,39 , and find a remarkable level of agreement given that we have practically no fitting parameters (meaning that most of the model parameters are in fact determined from independent experiments). The key observations that we wish to match are: (1) At low [ATP] there is an almost linear dependence of average rotation velocity with ATP concentration, practically independent on the degree of frictional resistance, with a slope ~6 ⋅ 10 6 rps/M 33,40 . (2) At high [ATP] the rotation rate saturates at a value determined by the friction resistance, in many different experiments this rotational friction coefficient has spanned the range 10 −4 -10 2 pN.nm.s 33  torque" is frequently quoted as 40 pN.nm 22,40 , although the earlier work of Junge et al. produced ~56 pN.nm 3 . We already explained the errors in interpreting a product ξ〈 ω〉 of a stochastic motor as an actual torque generated, nevertheless, the actual experimental values of ξ and 〈 ω〉 need to be matched. All of this happens on ATP hydrolysis and the binding rates are quoted in the range 2-6 ⋅ 10 7 M −1 s −1 by a number of different groups 33,[39][40][41] . These values are our initial 'target' in this paper.

Asymmetric Torsional Energy of the Coiled Coil
The central γ-shaft is a left-handed coiled coil of two α-helical proteins. Here we calculate how the elastic energy increases as a function of the angle of twist of such a coiled coil filament. Note that there are other, sometimes more sophisticated (and always more complex) elastic models of helical filaments 42,43 , but we choose to stay with very simple, qualitatively clear description.
We can regard the α-helix as an elastic filament with known characteristics. Its outer diameter is 1.2 nm; its persistence length has been extensively studied, providing a reasonably accurate value for the bending modulus B ≈ 570 pN ⋅ nm 2 (equivalent to a persistence length of ~100 nm at room temperature 44,45 ). Surprisingly for a cylindrical elastic filament, the torsional elasticity of an α-helix is greater: C ≈ (3/2)B with the corresponding torsional persistence length of ~150 nm. The 'surprise' is because in a homogeneous cylinder made of an isotropic elastic material, the torsional modulus is related to the bending modulus by an equation where σ is the Poisson ratio of the material 46 . So for an incompressible material one expects C = (2/3)B and for a typical Poisson ratio of crystalline solids (σ = 0.3) we would have C = 0.77B, while the observed factor of 3/2 implies a negative Poisson ratio σ = −1/3. This is a very unusual case in elasticity, often referred to as an "auxetic material" 47,48 . On the other hand, it may not be so surprising if we consider the main elastic elements in the α-helix are the hydrogen bonds in the outer shell of the cylinder, which are aligned with the helix axis: such a configuration has all the makings of a locally auxetic material. These facts and a discussion are well presented in a paper by Sun et al. 45 .
When a pair of α-helices is wound into a two-strand coiled coil tertiary structure, the axis of each of the α-helices makes a helical curve of its own, with a radius R = 0.46 nm and a pitch p in the range between 11 nm 49 and 14 nm 50 . For our calculations we shall take p = 11 nm. There is an extensive knowledge of the coiled coil geometry, going back to the work of Pauling and Crick in the 1950s 51,52 ; this geometry and its parameters are well summarized in recent papers by Sun et al. 53,54 . We might take a view that the α-helix is an elastic filament that is straight in equilibrium, but is forced to make a left-handed helical curve in space when it is twisting in the coiled coil configuration. In calculating the elastic energy we use the approach of Yamakawa 55 , by expressing the geometry of a helical curve via the two characteristic parameters: curvature and torsion (κ and τ), which are directly related to the radius R and pitch p of the coiled coil, Fig. 3: In equilibrium each α-helical filament has κ = 0 and τ = 0, but when they are twisted into the coiled coil, the two space curves r 1 and r 2 represent the centerlines of each filament: 1 2 where s is the arc length along the curves of α-helix centerline, ω = 2π/p is the rate of helical winding and q = pω/2π is the rate of advance of the curve along z-axis. Note that since the α-helix is essentially inextensible 45 2 is in place. The winding angle at top of the coiled coil is what has been acquired when s = L, that is, θ = ωL. The Frenet-Serret equations allow calculating the curvature and torsion for a given helix configuration: which shows that for such overconstrained space curves the curvature and torsion are not independent: The elastic energy of each such filament (as it is forced to wind into a coiled coil) as a function of deviation from its natural equilibrium is measured by κ and τ: The reason the centerline of an α-helix forms a helical curve in the tertiary structure of coiled coil is because a mechanical frustration is imposed on it by the hydrophobic bonding of residues with the parallel second α-helix 51,53,56 . A typical α-helix has a helical pitch of 0.6 nm and 3.6 residues per turn, making the step along its axis of h = 0.6/3.6 ≈ 0.167 nm per residue. Hence the length of the α-helical filament with N residues has a fixed value of L = Nh, and we already stated that it is essentially inextensible 45 . However, when a second α-helix is aligned parallel to it, the hydrophobic pairing of apolar residues of the so-called heptad repeat (7:2 pairs) 52,57 forms a contact line or "seam" between two α-helices, which has 3.5 residues per turn, bonding the two residues (a,d) of each heptad (the study of coiled coil elasticity 42 calls this line the "interface curve"). Since the length along this contact line of paired residues on the side of the cylindrical α-helical filament is shorter than the natural length of the centerline of this elastic filament, the rest of the filament coils into a helix -a phenomenon familiar in telephone cords, plant tendrils and curly hair. The key to the resulting geometry is that the height of the coiled coil (defined as H in Fig. 3: the length along the seam line) is shorter than the contour length and their observed ratio is H/L = 3.5/3.6 = 0.97. Since the height H of the γ-shaft is known from experiment 6 (H ≈ 6.5 nm), these relations determine the equilibrium pitch (p = 11 nm) or the equilibrium twist (θ 0 = 2πH/p = 3.71 rad, or 210°) of the coiled coil pair in the γ-shaft. The elastically active length of α-helices is Scientific RepoRts | 6:28180 | DOI: 10.1038/srep28180 L ≈ 6.72 nm, which makes approximately m = 2(6.72/0.167)/7 ≈ 15 hydrophobic bonds along the seam line of the γ-shaft, assuming that two residues (a,d) of each heptad are paired.
However, we cannot assume the hydrophobically-bonded pairs of matching (a,d) heptad members experience a harmonic potential with respect to stretching/compressing the seam line. This line is made of a sequence of pairs of hydrogen bonds in each α-helix, linked by the matching (a,d) heptad members and nearly parallel to the line itself: it is much harder to stretch the sequence of bonded pairs than to compress it. In stretching, or elongating the distance between residues along the axis of the α-helix, the existing hydrogen bonds that are aligned along the axis of the α-helix need to be stretched, which is a hard proposition. In contrast, the shortening of this line can be much more easily accommodated by tilting of the α-helical bonds away from being parallel to the helical centerline. One can see an analogy with the classical Euler problem of strut elasticity: it is hard to stretch an elastic rod, but on compression it can easily buckle and thus respond with a much lower force. The expression for such a stretching/compressing energy, as function of the length of interface curve, x, given its equilibrium length H 0 , and the overall length of the filament L, should take the form analogous to the general energy of a stiff filament 58 : This stretching energy E seam is plotted in Fig. 4(a), illustrating a typical response of a stiff filament, which is hard to stretch but easy to buckle. The experimentally observed coiled coil height H = 6.5 nm (or the corresponding twisting angle θ 0 = 3.71 rad) are obtained by converting from the seam line length x to the coiled coil twisting angle θ (via θ = − x L R ( ) 2 2 ) and adding the stretching energy E seam (θ) to the twisting energy of two α-helices 2E α (θ), which is plotted as a total in Fig. 4(b). The minimum of this energy defines the equilibrium shape of the coiled coil: a pair of naturally straight elastic filaments (α-helices) frustrated (contorted) by the added elastic energy of the seam line.
It may seem that there are many unknown (free) parameters in this simple elastic model, but in fact the pub-   distribution of angle θ of the freely rotationalkly fluctuating γ-shaft and obtained a Gaussian shape, which allows to determine the modulus directly, giving a value ~0.4 pN.nm/deg 2 = 1300 pN.nm/rad 2 (in fact, a range of values within this order of magnitude, depending on which segment of γ-shaft is examined). A combination of this modulus (i.e. the curvature of the torsional energy near the minimum θ 0 , and the position of this minimum at θ 0 = 3.71 rad, are enough to determine the parameters mε ≈ 250k B T (or 147 kcal/mol, i.e. ε = 9.8 kcal/ mol = 69 pN.nm per bonded pair of the heptad), and H 0 ≈ 6.33 nm. So in equilibrium, the seam line is slightly stretched (to 6.5 nm) at the expense of coiling the two α-helix filaments. Figure 5 zooms into the region near the minimum of this torsional energy, which we label as E 1 (θ) as the first level of the two-state model for the motor, and we can clearly see the asymmetry: the energy rises steeper for unwinding the coiled coil (i.e. attempting to stretch the seam line). The torsional fluctuations of the coiled coil about this minimum are markedly asymmetric, which is the key ingredient of our model of rotary motor. It is remarkable in the hindsight that the experimental results of Junge et al. are in fact showing the slight asymmetry of their probability distributions! The analysis in this section provides the values of asymmetry about the minimum of torsional elastic energy. As we have discussed earlier and illustrated in Fig. 1(b), there is a 120°-periodicity of the cam end of the γ-shaft confinement inside the α 3 β 3 cavity. Taking this 'window' of 120° = 2π/3 for each cycle, the asymmetry of the torsional potential makes the angular intervals: a ≈ 1.2 and b ≈ 0.9 [rad], as illustrated in Fig. 5. Of course, the sum a + b = 2π/3, but the crucial parameter is the difference: a − b ≈ 0.3 [rad], or 17°. The height of energy barrier at the boundaries of each angular cycle is Δ E 0 ≈ 136k B T within this model. This is a high energy barrier, so one can be assured that no spontaneous 'slippage' of the of the γ-shaft can occur while its cam is confined inside the α 3 β 3 cavity.

Two-state Stochastic Brownian Ratchet Motor
In this section we re-write the original model of Jülicher, Ajdari and Prost 27,34 for the case of rotary motion, in the simplest possible form. In thermodynamic equilibrium, without external torque, the asymmetry of potential energy E 1 (θ) is irrelevant: the detailed balance ensures that the diffusion flux in both directions must be equal. So in spite of the rotational symmetry broken by the slanted potential energy in each period, the thermal motion alone cannot produce unidirectional rotation. In this 120°-periodic low-energy state E 1 (θ) the bulge in the γ-shaft (the cam tooth) is locked in one of the three positions around the circle -and the c-ring only fluctuates confined near the minimum of this torsional elastic energy. These are the fluctuations measured by Junge et al. 7,35,38 and simulated by Czub and Grubmüller 36 , which we discussed in the previous section.
When an ATP molecule is hydrolyzed in the α 3 β 3 complex (with the binding rate k on ), the confinement of the γ-shaft cam changes: here we assume it is released and becomes free to rotate inside the α 3 β 3 cavity. This means that the system enters the "excited state" E 2 with no rotational confinement -this assumption is an obvious limitation made to simplify the model and make it easily tractable, and we shall discuss later in the text how relevant it is (it turns out that it is not). Free rotational diffusion of the c-ring and γ-shaft occurs then, with the Gaussian distribution of probability p(Δ θ) ~ exp[− (θ − θ 0 ) 2 /4Dt]. In the wild type ATPase embedded in its proper membrane, the rotational diffusion constant of the c-ring rotor was measured by several authors and reported, e.g., by Elston et al. 11 : D ≈ 2 ⋅ 10 4 rad 2 /s. In many famous experiments observing the F 1 rotation in vitro, people were attaching either an actin filament or a gold bead in place of the c-ring -and observing the rotation of such a marker against viscous friction in solution 21,33,40 . The diffusion constant D = k B T/ξ has been varied over nearly 6 orders of magnitude, depending on the geometry of the marker 33 . For instance, at room temperature the 40-nm gold bead of Yasuda et al. 33 has ξ = 2 ⋅ 10 −4 pN.nm.s, giving the diffusion constant D ≈ 2.1 ⋅ 10 4 rad 2 /s: and almost exact match with the membrane friction of the c-ring. Using a 1-μm long actin filament as a marker gives ξ ≈ 2 pN.nm.s 21 , producing D ≈ 2.1 rad 2 /s. The authors in the past have incorrectly used the idea of "no-load" velocity at very low friction, but we see that all the increasing friction does is to slow down the free rotational diffusion in the upper unconstrained state R 2 . We shall examine the effect of changing friction on the motor velocity later in the text.
After a characteristic life-time τ 2 of this isomerized state where the γ-shaft could freely diffuse, the α 3 β 3 complex returns to its natural state and the cam returns back to a locked position; this may be the same position as when it was released -but due to the free rotational diffusion in the E 2 state, it will more frequently drop into the neighboring position to the left (the closest to the original, as b < a). We can estimate the average angular velocity of the resulting rotation as the ratio of angular step (a + b = 120°) and the total time of this cycle where the probabilities of dropping into the right-and left-side pockets of E 1 are determined by the corresponding error functions from the limited integrals of p(Δ θ). In the limit when the torsional asymmetry (a − b) is small, this velocity takes a more simple form:  Figure 6 shows the pictorial scheme of this process, while the plots in Fig. 7(a,b) show how the average anticlockwise "no-load" velocity 〈 ω 0 〉 depends on the key parameters. First of all, there is the life-time in the isomerized free-rotation state τ 2 . This time is not generally known, and the only close measurement reported by Furuike et al. 39 is that of a dwell time of several milliseconds (also see 25 ). We assume that the actual life time of the isomerized β-subunit is much shorter. The plot 7(a) is made for the low-friction rotation corresponding to a 40 nm gold bead marker, or the wild type motor, D = 2 ⋅ 10 4 rad 2 /s, and several realistic values of [ATP] given the k on = 4 ⋅ 10 7 M −1 s −1 . First of all, we see that our concept model, in spite of all its simplifications, is quantitatively producing the values of average rotation velocity in the observed range (spinning at over a 100 rps is clearly realistic). On the other hand, the dependence on the poorly known life time τ 2 is very sharp and if one 'misses' the optimal value, the rotation velocity drops rapidly. This dependence is easily understood: at very short life-time the cam cannot diffuse too far from its original position in the E 1 periodic potential, and it's most likely to re-engage in the same position giving no net rotation -at very long τ 2 the cam has time to diffuse very far and it becomes relatively equal-probability to re-engage to the left and to the right. It is clear that this is a crucial parameter for evolution to adjust for any given motor and its specific use in the cell.
Plot 7(b) shows the dependence on [ATP] concentration, which enters in our model as the rate of E 1 → E 2 transition. As expected, and in good agreement with 'canonical' experiments 33,40 , at low ATP concentration the average velocity varies linearly with concentration -while at high concentration it saturates at a maximal value ω max which is controlled by the rate of the reverse E 2 → E 1 transition, 1/τ 2 . The change between the two regimes occurs at k on [ATP] ≈ 1/τ 2 . It is common in the literature to fit this type of curve to a Michelis-Menten equation, which is indeed what Eq. (5) gives, with (still in the nominal units of [rad/s]) maintaining the simplified expanded form for the small torsional asymmetry (a − b) ≪ 1. Note that the ratio 2a 2 /D is approximately the time for a non-engaged γ-shaft to diffuse a single 120° period, the ratio of which to τ 2 is what effectively controls the magnitude of these expressions. Let us now examine how this motor operates in different friction regimes. As explained above, for the stochastic motor like this the (rotational) friction does not provide a 'load' , but simply acts through the changing rate of diffusion, D = k B T/ξ. Wishing to compare with the famous experiments of Yasuda et al. 21,33 which studied a wide range of markers providing different rotational friction constant, we plot the dependence of average velocity 〈 ω 0 〉 on this friction coefficient (assuming room temperature of 23 °C), rather than the diffusion constant D, see Fig. 8. The log-linear version of this plot (a) shows the exponential decay of 'no-load' 〈 ω 0 〉 with increasing friction, which is in fact obvious from examining Eqs (5) and (6). The log-log plot (b) enhances the small change in 〈 ω 0 〉 at very low friction (perhaps at unreasonably low value used merely for completeness of exposure). This may seem not in perfect agreement with the results of Yasuda et al. 33 that appear to plateau at ξ → 0. However, we would point that our model was ultimately simplified for clarity of concept, and it certainly underestimates the friction in the 'free-diffusion' excited state (when we took that the γ-shaft has no constraint whatsoever). Equally, there are difficulties to actually achieve vanishing friction in experiment, so this small decrease might not have been noticed (as well as the reliability of accurately determine the friction constant limited).
Yasuda et al. 21,33 fit the curve of 〈 ω 0 〉 vs. friction constant ξ by an interpolated formula: 〈 ω 0 〉 = (1/V noload + 2πξ/M) −1 , where their 'assumed torque' was 40 pN.nm. We already explained that a stochastic molecular motor does not exert a torque in the unloaded state -it just makes more random steps in one direction than in the other. In the high-friction limit when ξ ≫ k B Tτ 2 , our simplified Eq. (5) predicts an exponential drop in the average rotation velocity 〈 ω 0 〉 , while the data of Yasuda et al. 21,33 is probably going broadly over the crossover region.

F 1 Motor Operating Against External Torque
In real circumstances, the F 0 motor driven by the pmf, originally described by Elston et al. 11 and reviewed in several topical reviews on the subject, exerts a torque on the c-ring, which is passed through the γ-shaft to the Scientific RepoRts | 6:28180 | DOI: 10.1038/srep28180 confining region inside the α 3 β 3 complex. This torque, or the turning moment M, is in the direction opposite to the natural anticlockwise rotation of the F 1 motor, and so the two have to compete. Figure 9 shows the pictorial effect of such an external torque on the F 1 operation, slanting both energy states E 1 and E 2 and adding two additional contributions to the average angular velocity of the shaft: 〈 ω(M)〉 . One of these is a drift in the clockwise direction in the free-diffusion excited state E 2 , the other is diffusion to the right in the ground-state periodic potential E 1 (θ) because the energy barriers for the forward and backward motion are now different. The first effect (of the free diffusion envelope drift) is accounted for by the shift in the Gaussian exponent of the probability p(Δ θ) = exp[− (θ − θ 0 ) 2 /4Dt + Mθ/2k B T] and is naturally extending the Brownian ratchet approach that has led to Eq. (5). The second effect of the diffusion in a tilted periodic potential is a classical problem carefully studied in the textbook by Nelson 59 , improving on the original Feynman's treatment of the ratchet and pawl problem 60 . The simple Feynman formula has the exponentially diverging velocity at large torque, while the Nelson approach correctly accounts for the viscous drag and results with a linear 'force-velocity' relationship in the overdamped Figure 6. A scheme of two-state Brownian ratchet motor. In the low-energy state the bulge in the γ-shaft is locked in one of the three positions 120° around the circle -and the c-ring only fluctuates confined near the minimum of this torsional elastic energy. In thermal equilibrium, there can be no unidirectional motion in spite of the chiral anisotropy. An energy input from ATP hydrolysis causes α 3 β 3 isomerization and can disengage the confinement, allowing the c-ring to fluctuate freely while in this unconstrained excited state. The rate of E 1 → E 2 transition is k on [ATP], which is an activation function of ATP chemical potential Δ μ, while the rate of the reverse E 2 → E 1 transition is controlled by the life-time of the isomerized state of β-subunit (τ 2 ). On return to its cam-engaged state, the coiled coil has a higher probability (p − ) to end up at − 120° from the original state, than in the + 120° state (p + ), which results in the average anticlockwise rotation. Scientific RepoRts | 6:28180 | DOI: 10.1038/srep28180 molecular system. The formula for the S-ratchet derived in 59 has to be modified because there are two slopes of the potential from its minimum (towards a and b), so the final expression becomes a little bit more involved, although the underlying principle of diffusive motion across the effective periodic landscape remains the same. It is also instructive to compare the two versions of this part of this clockwise angular velocity at small torque M, i.e. in the linear response regime: However, for our practical purposes in this paper, it turns out that with the energy barrier ∆  E k T 136 B 0 and the external torque not exceeding  M 50 pN.nm (which is equal to 12k B T) this contribution to the motor action is completely negligible.
The drift of the free-diffusing distribution can be directly incorporated into the expression for the average angular velocity due to the Prost mechanism 34 , shifting the mid-point of the spreading free-diffusion distribution p(θ, t) in the excited state E 2 , giving the expression, which at M = 0 reduces to Eq. (5): The plots in Fig. 10 illustrate the predictions of Eq. (8) when the F 1 motor is subjected to a counter-torque M (clockwise, arising from the F 0 motor). Both plots are computed for the 'near-optimal' value of excited state life-time τ 2 = 50 μs, see Fig. 7(a) and the friction constant ξ = 2 ⋅ 10 −4 pN.nm.s, corresponding to the rotational diffusion D = 2 ⋅ 10 4 rad 2 /s (as in the wild type, or low-friction 40-nm gold bead of Yasuda et al. 33 ). The plot 10(a) shows a 'big picture' for a wide range of torques, but not permitting to see what happens around M = 0. We see that the external torque essentially drives the shaft -in the opposite clockwise directiion showing as the negative 〈 ω(M)〉 in the plot, when the torque is against the natural anticlockwise direction of the F 1 spin, or in the positive anticlockwise direction when the external torque works in the same direction as F 1 . At high enough torque the velocity saturates at a constant plateau value, which is determined by the rate of ATP binding (we continue using the 'standard' value for the specific binding rate k on = 4 ⋅ 10 7 M −1 s −1 ), and also the excited state life-time τ 2 .
The zoomed-in plot in Fig. 10 = 0.7 pN.nm (this value does not depend on κ on [ATP], or the energy barrier Δ E 0 , because it is the drift in the assumed 'free' state E 2 that M affects). This is a very low torque needed to stop the F 1 , hence a fully-driven F 0 pmf motor will always 'win' and rotate the γ-shaft clockwise; only a very weak pmf and high [ATP] would allow the F 1 motor to drive the c-ring anticlockwise and increase the pH imbalance by pumping H + ions across the membrane against the concentration gradient. We will discuss the F 1 'efficiency' in the next section.
An important observation we can make from Fig. 11(a) is that the whole ATPase must have some amount of ATP hydrolysing, in our model -generating 'excitation' out of the confinement in the periodic potential E 1 (θ). When there is no or little ATP-activation, the motor will not turn even under a significant torque. This is because the elastic energy barrier Δ E 0 is so high that the moderate external torque cannot break through these barriers and force the backward motion. However, once the rate of ATP binding increases above 100-1000 s −1 , the frequent periods of un-constrained γ-shaft allow a rapid backward drift in the excited state E 2 . Another unexpected result of our model is illustrated in Fig. 11(b), pointing that at very high friction the counter-torque M has very little effect: the clockwise drift in the un-constrained state is too slow, while the forward (anticlockwise) bias is determined by the fixed elastic asymmetry of E 1 (θ).

Discussion
Let us now discuss the energy balance of the F 1 motor operation and the motor efficiency. It is a traditional question to ask, when dealing with an engine, and in the case of F 1 it has several misconceptions in the literature that we would like to address.
The chemical energy obtained from the hydrolysis of one ATP molecule under intracellular conditions is given by the expression: = 1 mM this gives an estimate ⋅ Q 3 10 6 pN.nm/s. The useful work produced by the F 1 motor in its active operation regime is determined by the average angular speed 〈 ω〉 against the counter-torque M, while it rotates anticlockwise (i.e. able to drive against M). In contrast, passive regimes of this motor are those where the average rotation spin is along the direction of external torque, in which case the external mechanical work performed on F 1 is dissipated into heat (or leads to the ADP + Pi → ATP synthesis and chemical energy storage -in that case the useful work of the F 0 motor is of interest, see 11 for detail). These regimes are clear in Fig. 10, where the active regime is in the sector of positive M and 〈 ω〉 in our notation. In this case the useful work per unit time is   The resulting calculated 'efficiency' of the F 1 motor is very low, contrary to many statements in the literatureand contrary to a generic expectation for a biological machine to be very efficient. However, we must ask: what is the F 1 motor for? It counterpart F 0 has a clear purpose: to 'recycle' ADP, i.e. to drive the γ-shaft clockwise and 'forcefully' induce the ADP + Pi → ATP synthesis in the β-subunits. F 1 does not have a purpose to be 'strong' enough to out-perform the working F 0 and induce the anticlockwise rotation: in natural conditions it should only 'win' when the F 0 is dormant, i.e. there is an insufficient pmf to drive it. This is consistent with a very low vale of 'stall' torque M stall calculated earlier. Therefore it is misleading to evaluate its efficiency by counting the useful mechanical work against the counter-torque M. The purpose of F 1 operation is to rebuild the H + gradient across the membrane, and the relevant efficiency has to be evaluated with this 'useful outcome' in mind (not with the mechanical work against the counter-torque M as we had to do in the estimate above). This is discussed in detail by G. Oster et al. 11 in their theory of F 0 motor. Given that 4 protons pass through the rotor per ATP step against the membrane potential Δ pH, the useful work is 4eΔ pH per step, where e is the unit charge. Equivalently, the rate of work π ω = ⋅∆  W e (12 /2 ) p pH, where the average velocity is in [radian/s] units. When there is no H + gradient at all, this work vanishes and the efficiency is very low. When Δ pH reaches the value to generate  M 1 stall pN.nm from the F 0 , the rotation stops and the useful work rate vanishes again. The maximum efficiency of proton transport is, again, in between these two limits: η = = ∆   W Q e / 4 p p pH/Δ G 1 . We do not know at which Δ pH this occurs, but to test the values let us take an example value Δ pH = 100 mV: this gives 4eΔ pH = 16k B T and η  75% given the values of chemical energy Δ G 1 discussed above. Obviously in most cases the membrane potential will be much lower than this, because it is to precisely rebuild this proton gradient that F 1 will work. This work has presented a self-consistent analytical model of the ATP-driven F 1 motor of the ATP synthase complex: something that has been missing in understanding this remarkable molecular machine (in contrast, the pmf-driven F 0 motor has been understood well for many years). It is important to emphasize that our physical model does not use any detailed molecular structure of the machine -it only needs several key elements that are present and not disputed: the two-strand coiled-coil filament providing the torsional elastic bias, the 120° periodicity of the γ-shaft confinement, and the notion that on ATP binding this confinement is released. These elements, with the added Brownian motion, are enough to produce a working physical model of the rotary motor.
The most important and new contribution is the finding that the torsional elastic energy of the coiled coil filament (such as the γ-shaft here) is asymmetric. The very small fluctuations about the minimum of the this energy  Fig. 7(b). As the opposite torque increases the spinning direction reverses into the negative (clockwise) velocity: the F 0 motor takes over. (b) The effect of changing friction, again plotting the average velocity for several values of M, τ 2 = 50 μs, and [ATP] = 1 mM (corresponding to the binding rate k on [ATP] = 4 ⋅ 10 4 s −1 ). The M = 0 curve here is the same as the 50 μs curve in Fig. 8(b). Increasing counter-torque brings forward the 'stall' point of the motor. are described by a proper linear response, with the torsional modulus well-studied in many experiments and simulations. However, as soon as the twisting of the coiled coil becomes more significant, the asymmetry between clockwise and anticlockwise motion becomes relevant. In fact, on careful examination of the experimental probability distributions of the free-fluctuating filament 36,38 one can actually see this asymmetry at larger angles, so we believe this is a genuine effect -probably with applications other than this particular motor.
After over 20 years of detailed experimental studies of ATPase we were able to find the data for all elements of the model in these experiments and so there are really no free parameters -with a possible exception of τ 2 : the life time of the 'excited' second state when the cam-end of the γ-shaft is not tightly confined inside the α 3 β 3 cavity. We were tempted to take the value of τ 2 ≈ 10 μs that optimizes the motor velocity almost universally at any ATP concentration, but we are aware that the only related measurement we know (by Furuike et al.) -that of a dwell time -is much longer, so we used τ 2 ≈ 50 μs in all subsequent graphs. Still, we find it remarkable that our model (with all its simplifications and streamlining) is predicting the rotation velocities and torque, and their dependence on ATP, with quantitative accuracy.
The greatest 'weakness' of our model is the assumption that the γ-shaft is unconstrained in the ATP-binded state and is completely free to rotate (only against the viscous friction provided by the c-ring in the membrane, or any artificially added construct in-vitro). We realize that in reality the torsional energy of this 'excited' state E 2 is likely to be θ-dependent as well, which would affect the effective diffusion constant we used. One might say that we always underestimate the friction in the model, although this is a less straightforward non-linear effect. This would explain the small but noticeable differences in our results when plotted against friction constant ξ in Fig. 8 and the experimental results of Yasuda et al. (which were discussed in the text). Our aim was the maximally streamlined, simplified model that would highlight the key principles of ATP-driven rotary motor based on the torsionally-asymmetric filament axle, and further work should certainly improve on its various aspects using the details of molecular structure in that state.
Although in this paper we discussed the F-type ATPase, there are of course very similar rotary motors related to it -in particular the vacuolar V-ATPases 61,62 . Without going into details, the key structural elements required for our motor mechanism are recognizably the same in the V-ATPase: the stator with a coiled coil drive filament inserted and constrained in its cavity, attached to a ring rotor embedded in the membrane. It would be interesting to re-examine the action of these other rotary motors driven by ATP in view of the findings here. The H + or the Na + pumping resulting from driving this rotor may have specifics in different biological situations, but the fundamental principle of the rotary motor driven by Brownian motion of the torsionally asymmetric filament switching between two states -confined and free-should be universal.