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

**|**Scientific Reports**|**PMC4913325

Formats

Article sections

- Abstract
- State of the Art
- Asymmetric Torsional Energy of the Coiled Coil
- Two-state Stochastic Brownian Ratchet Motor
- F1 Motor Operating Against External Torque
- Discussion
- Additional Information
- References

Authors

Related links

Sci Rep. 2016; 6: 28180.

Published online 2016 June 20. doi: 10.1038/srep28180

PMCID: PMC4913325

Received 2016 February 2; Accepted 2016 May 31.

Copyright © 2016, Macmillan Publishers Limited

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

This article has been cited by other articles in PMC.

F_{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. 400rad/s anticlockwise at 4mM 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 H^{+} concentration gradient on the two sides of the membrane). The pmf arising from this concentration gradient is able to deliver the sufficient torque on the *c*-ring, which is then passed via the *γ*-shaft into the cavity of the *α*_{3}*β*_{3} complex and forces the 3-stage sequential structural isomerization (opening) of the *β*-subunits to facilitate the capture of ADP and phosphate ion, the enzymatic reaction ADP+Pi, and the release of ATP.

The action mechanism of the F_{1} motor is largely unknown, in spite of many attempts to explain it and the corresponding large literature. It is known to be powered by the ATP hydrolysis in the *β*-subunits of the *α*_{3}*β*_{3} complex (i.e. the reverse of the ADP phosphorilation sequence) and produces the opposite (anticlockwise) rotation of the *c*-ring, which is therefore forced to rotate in spite of the energy barrier of non-protonated Asp61, and thus deliver H^{+} ions into the high-concentration channel, while taking them away from the low-concentration volume. In this way ATPase can pump protons and build up the transmembrane pH gradient using the energy stored in ATP. A fully analogous V1 ATPase system also actis to transport protons against the chemical potential gradient across a membrane^{12}^{,13}^{,14}, while other variants of F-type ATPase are used to pump Na^{+} ions^{15}. There are many interesting and detailed papers in the literature, attempting to discern how this anticlockwise torque can be generated in the junction between the *γ*-shaft and *α*_{3}*β*_{3} complex, but none are satisfactory. They often start from the original mechanistic “binding change” ideas of P. Boyer^{16} and develop towards the power-stroke mechanism of Wang and Oster^{17}^{,18}. The latter models in this vein are very complex, with up to 64 states and many coupled stochastic variables^{17}^{,19}. With many free parameters they do succeed at matching the phenomenology of the F_{1} motor^{20}. There nonetheless remains a need to produce a concept model which explains the physical origin of the unidirectional motion and corresponding useful work generated by F_{1}, and many analogous ATP-driven rotary motors, and the direction of this rotation in terms of the microscopic configuration and the physics of the system. There is also an often-quoted “100% efficiency” of the F_{1}-ATPase motor^{3}^{,21}^{,22}^{,23}, a phenomenon that would contradict basic thermodynamics. The origin of such a confusion is that people incorrectly claim that the “useful work” the motor produces is the work against viscous drag – and thus estimate the motor efficiency from the ratio of this to the rate of ATP hydrolysis. One will obviously get 100% every time, if the experiment is done correctly, because all of the ATP-released energy eventually dissipates into heat. We will show below that the energy influx from the ATP hydrolysis only serves to switch between the states of the system, while the driving power arises from the Brownian motion itself. We will also illustrate that the ideas of “power stroke” (i.e. a short period of time when a unidirectional torque is exerted by the motor) do not correspond to observations – instead, the stochastic Brownian-ratchet motor generates the average motion (but not torque).

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 *ξ**ω*=40pN.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 40pN.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 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 2mM^{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 and33 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 and55 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 ~610^{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}. No doubt motivated by Michaelis-Menten, people frequently use a fitting function *ω*=(1/*v*_{noload}+2*πξ*/*M*)^{−1} where *M* is the “assumed torque” generated by the motor. At a high [ATP]=2mM and low friction, the value of *v*_{noload} was found to be ~130 rps^{33}. (3) The value of the “assumed torque” is frequently quoted as 40 pN.nm^{22}^{,40}, although the earlier work of Junge *et al*. produced ~56pN.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–610^{7}M^{−1}s^{−1} by a number of different groups^{33}^{,39}^{,40}^{,41}. These values are our initial ‘target’ in this paper.

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.2nm; its persistence length has been extensively studied, providing a reasonably accurate value for the bending modulus *B*≈570 pNnm^{2} (equivalent to a persistence length of ~100nm 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 ~150nm. 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 *C*=*B*/(1+*σ*), 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.77*B*, 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.46nm and a pitch *p* in the range between 11nm^{49} and 14nm^{50}. For our calculations we shall take *p*=11nm. 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:

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}, the constraint or 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.6nm and 3.6 residues per turn, making the step along its axis of *h*=0.6/3.6≈0.167nm 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.5nm), these relations determine the equilibrium pitch (*p*=11nm) 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 *L*≈6.72nm, 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.5nm (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 ) and adding the stretching energy *E*_{seam}(*θ*) to the twisting energy of two *α*-helices 2*E*_{α}(*θ*), 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 published experimental results constrain them very accurately. Junge *et al*. have measured the linear torsional modulus of the *γ*-shaft by observing free thermal fluctuations of its end^{7}^{,35}^{,38}. A wide range of values is reported, depending on which section of the filament one is tracking in elaborate and delicate experiments, but the order of magnitude remains well defined. A bigger problem is that the torsional modulus has to be presented in the units of [energy/angle^{2}] while all the reported values for the ‘modulus’ are given in the dimensions of energy alone [pN.nm]. The MD simulation study of *γ*-shaft elasticity^{36} also has a typo, presenting the modulus in [pN.nm/rad]. In spite of all that, the core experimental result of Junge *et al*. is clear and unambiguous: they have mapped the probability 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.4pN.nm/deg^{2}=1300pN.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ε*≈250*k*_{B}*T* (or 147kcal/mol, i.e. *ε*=9.8kcal/mol=69pN.nm per bonded pair of the heptad), and *H*_{0}≈6.33nm. So in equilibrium, the seam line is slightly stretched (to 6.5nm) 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}≈136*k*_{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.

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}/4*Dt*]. 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*≈210^{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 *ξ*=210^{−4} pN.nm.s, giving the diffusion constant *D*≈2.110^{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.1rad^{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 (Δ*t*=1/*k*_{on}[*ATP*]+*τ*_{2}), giving *ω*_{0}=(*a*+*b*)(*p*_{+}−*p*_{−})/Δ*t*, 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:

where the subscript in *ω*_{0} indicates the absence of an external torque: a true “no-load” condition of the motor. This rotation velocity is nominally measured in [rad/s]; to convert it into [rps] one has to divide it by 2*π*. Here the rate of ATP binding *k*_{on}[*ATP*]=*κ*_{0}exp[Δ*μ*_{ATP}/*k*_{B}*T*] can be expressed via the chemical potential Δ*μ*=*k*_{B}*T*ln[*ATP*]. Taking the specific rate of ATP binding *k*_{on}=410^{7}M^{−1}s^{−1} allows converting the ATP concentration of 2mM to the actual rate of binding given by the product *k*_{on}[*ATP*]≈810^{4}s^{−1}, or the time interval between the ATP binding events of ~13*μ*s.

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 40nm gold bead marker, or the wild type motor, *D*=210^{4}rad^{2}/s, and several realistic values of [ATP] given the *k*_{on}=410^{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 2*a*^{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 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.

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 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}/4*Dt*+*Mθ*/2*k*_{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 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 and the external torque not exceeding pN.nm (which is equal to 12*k*_{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 *ξ*=210^{−4} pN.nm.s, corresponding to the rotational diffusion *D*=210^{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}=410^{7}M^{−1}s^{−1}), and also the excited state life-time *τ*_{2}.

The zoomed-in plot in Fig. 10(b) shows the details of what happens at no or very low external torque. The ‘no-load’ values of the average angular velocity *ω*_{0} depend on the [ATP] as well as on the friction constant as the results in Figs 7 and and88 have shown. Since the ‘stall’ condition occurs at quite a low value of external torque, a simple expansion of Eq. (8) is justified, and it gives an approximate expression: =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–1000s^{−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}(*θ*).

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: Δ*G*_{1}=Δ*G*_{0}+*k*_{B}*T*ln([ADP][Pi]/[ATP]), with Δ*G*_{0}=−50pN.nm (or equivalently: 12*k*_{B}*T*) is the free-energy change per molecule of ATP hydrolysis at pH=7. Taking intracellular concentrations [ATP]≈[Pi]1mM^{22}, we obtain the chemical energy released per step: Δ*G*_{1}=−97pN.nm () at [ADP]=10*μ*M, Δ*G*_{1}=−88pN.nm (=21.2*k*_{B}*T*) at [ADP]=0.1mM, and Δ*G*_{1}=−78pN.nm (=18.9*k*_{B}*T*) at [ADP]=1mM, i.e. not changing very dramatically. Per unit time, the chemical energy input into F_{1} operation is therefore equal to: [ATP]Δ*G*_{1}; for [ATP]=1mM this gives an estimate 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 , where Eq. (8) has to be employed and the velocity in [radian/s] units must be used. The motor efficiency in this mode of operation is . The useful work (and the efficiency *η*) is zero when *M*=0 (in ‘no-load’ conditions), and equally zero at the ‘stall’ point where *ω*(*M*)=0. The maximum of the useful work lies between these two limits, and Fig. 12 shows its value and position for the wild-type level of friction ( pN.nm.s), the life-time of the excited state of *β*-subunit isomerization *τ*_{2}=50*μ*s, and several values of [ATP]. For instance, when [ATP]=[ADP]=1mM, the maximum rate of useful work rate is pN.nm/s.

The resulting calculated ‘efficiency’ of the F_{1} motor is very low, contrary to many statements in the literature– and 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 4*e*ΔpH per step, where *e* is the unit charge. Equivalently, the rate of work 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 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: pH/Δ*G*_{1}. We do not know at which ΔpH this occurs, but to test the values let us take an example value ΔpH=100mV: this gives 4*e*ΔpH=16*k*_{B}*T* and 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 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.

**How to cite this article**: Kulish, O. *et al*. F_{1} rotary motor of ATP synthase is driven by the torsionally-asymmetric drive shaft. *Sci. Rep.*
**6**, 28180; doi: 10.1038/srep28180 (2016).

The authors have benefited from extensive discussions with J. R. Blundell, C. Prior, and G. Fraser, as well as the conceptual input from J. E. Walker (who has originally suggested that the torsional energy of the *γ*–shaft might be asymmetric). This work has been funded by the {100+100+100} program by the Ukrainian Government, and the EPSRC Critical Mass Grant for Cambridge Theoretical Condensed Matter EP/J017639.

**Author Contributions** E.M.T. developed the concept, O.K., A.D.W. and E.M.T. carried out the analysis. All authors contributed to writing the paper.

- Stock D., Gibbons C., Arechaga I., Leslie A. G. W. & Walker J. E. The rotary mechanism of ATP synthase. Curr. Opin. Struct. Biol. 10, 672–679 (2000). [PubMed]
- Yoshida M., Muneyuki E. & Hisabori T. ATP synthase - a marvellous rotary engine of the cell. Nature Rev. Mol. Cell Biol. 2, 669–677 (2001). [PubMed]
- Junge W., Sielaff H. & Engelbrecht S.
Torque generation and elastic power transmission in the rotary F
_{0}F_{1}-ATPase. Nature 459, 364–370 (2009). [PubMed] - Noji H., Yasuda R., Yoshida M. & Kinosita K.
Direct observation of the rotation of F
_{1}-ATPase. Nature 386, 299–302 (1997). [PubMed] - Itoh H.
et al. . Mechanically driven ATP synthesis by F
_{1}-ATPase. Nature 427, 465–468 (2004). [PubMed] - Hausrath A. C., Grüber G., Matthews B. W. & Capaldi R. A.
Structural features of the
*γ*subunit of the escherichia coli F_{1}ATPase revealed by a 4.4-Å resolution map obtained by x-ray crystallography. Proc. Natl. Acad. Sci. USA 96, 3697–3702 (1999). [PubMed] - Sielaff H., Rennekamp H., Engelbrecht S. & Junge W.
Functional halt positions of rotary F
_{0}F_{1}-ATPase correlated with crystal structures. Biophys. J. 95, 4979–4987 (2008). [PubMed] - Sun S., Chandler D., Dinner A. & Oster G.
Elastic energy storage in
*β*-sheets with application to F_{1}-ATPase. Eur. Biophys. J. 32, 676–683 (2003). [PubMed] - Stock D., Leslie A. G. W. & Walker J. E. Molecular architecture of the rotary motor in ATP synthase. Science 286, 1700–1705 (1999). [PubMed]
- Abrahams J. P., Leslie A. G. W., Lutter R. & Walker J.
Structure at 2.8-Å resolution of F
_{1}-ATPase from bovine heart - mitochondria. Nature 370, 621–628 (1994). [PubMed] - Elston T., Wang H. & Oster G. Energy transduction in ATP synthase. Nature 391, 510–513 (1998). [PubMed]
- Imamura H.
et al. . Evidence for rotation of V
_{1}-ATPase. Proc. Natl. Acad. Sci USA 100, 2312–2315 (2003). [PubMed] - Beyenbach K. W. & Wieczorek H. The V-type H+ ATPase: molecular structure and function, physiological roles and regulation. J. Exp. Biol. 209, 577–589 (2006). [PubMed]
- Minagawa Y.
et al. . Basic properties of rotary dynamics of the molecular motor enterococcus hirae V
_{1}-ATPase. J. Biol. Chem. 288, 32700–32707 (2013). [PMC free article] [PubMed] - Thomas Meier1, K. D., W. W. P. D. Patrick Polzer2. Structure of the rotor ring of F-type Na+ − ATPase from ilyobacter tartaricus. Science 308, 659–662 (2005). [PubMed]
- Boyer P. D. The ATP synthase - a splendid molecular machine. Annu. Rev. Biochem. 66, 717–749 (1997). [PubMed]
- Wang H. & Oster G.
Energy transduction in the F
_{1}motor of ATP synthase. Nature 396, 279–282 (1998). [PubMed] - Wang H. & Oster G. Ratchets, power strokes, and molecular motors. Appl. Phys. A 75, 315–323 (2002).
- Gaspard P. & Gerritsma E.
The stochastic chemomechanics of the F
_{1}-ATPase molecular motor. J. Theor. Biol. 247, 672–686 (2007). [PubMed] - Kawaguchi K., Sasa S. & Sagawa T.
Nonequilibrium dissipation-free transport in F
_{1}ATPase and the thermodynamic role of asymmetric allosterism. Biophys. J. 106, 2450–2457 (2014). [PubMed] - Yasuda R., Noji H., Kinosita K. & Yoshida M.
F
_{1}-ATPase is a highly efficient molecular motor that rotates with discrete 120° steps. Cell 93, 1117–1124 (1998). [PubMed] - Kinosita K., Yasuda R., Noji H. & Adachi K. A rotary molecular motor that can work at near 100% efficiency. Philos. Trans. R. Soc. Lond. B 355, 473–489 (2000). [PMC free article] [PubMed]
- Toyabe S., Watanabe-Nakayama T., Okamoto T., Kudo S. & Muneyuki E.
Thermodynamic efficiency and mechanochemical coupling of F
_{1}-ATPase. Proc. Natl. Acad. Sci. USA 108, 17951–17956 (2011). [PubMed] - Okazaki K.-I. & Hummer G.
Phosphate release coupled to rotary motion of F
_{1}-ATPase. Proc. Natl. Acad. Sci. USA 110, 16468–16473 (2013). [PubMed] - Mekherjee S. & Warshel A.
Dissecting the role of the
*γ*-subunit in the rotary-chemical coupling and torque generation of F_{1}-ATPase. Proc. Natl. Acad. Sci. USA 112, 2746–2751 (2015). [PubMed] - Reimann P. Brownian motors: noisy transport far from equilibrium. Phys. Rep. 361, 57–265 (2002).
- Parmeggiani A., Jülicher F., Ajdari A. & Prost J. Energy transduction of isothermal ratchets. Phys. Rev. E 60, 2127–2140 (1999). [PubMed]
- Toyabe S. & Muneyuki E.
Single molecule thermodynamics of ATP synthesis by F
_{1}-ATPase. New J. Phys. 17, 015008 (2015). - Hayashi K., Ueno H., Iino R. & Noji H.
Fluctuation theorem applied to F
_{1}-ATPase. Phys. Rev. Lett. 104, 218103 (2010). [PubMed] - Watanabe-Nakayama T.
et al. . Effect of external torque on the ATP-driven rotation of F
_{1}-ATPase. Biochem. Biophys. Res. Comm. 366, 951–957 (2008). [PubMed] - Furuike S.
et al. . Axle-less F
_{1}-ATPase rotates in the correct direction. Science 319, 955–958 (2008). [PubMed] - Uchihashi T., Iino R., Ando T. & Noji H.
High-speed atomic force microscopy reveals rotary catalysis of rotorless F
_{1}-ATPase. Science 333, 755–758 (2011). [PubMed] - Yasuda R., Noji H., Yoshida M., Kinosita K. & Itoh H.
Resolution of distinct rotational substeps by submillisecond kinetic analysis of F
_{1}-ATPase. Nature 410, 898–904 (2001). [PubMed] - Jülicher F., Ajdari A. & Prost J. Modeling molecular motors. Rev. Mod. Phys. 69, 1269–1281 (1997).
- Sielaff H.
et al. . Domain compliance and elastic power transmission in rotary F
_{0}F_{1}-ATPase. Proc. Natl. Acad. Sci. USA 105, 17760–17765 (2008). [PubMed] - Czub J. & Grubmüller H.
Torsional elasticity and energetics of F
_{1}-ATPase). Proc. Natl. Acad. Sci. USA 108, 7408–7413 (2011). [PubMed] - Sun S. X., Wang H. & Oster G.
Asymmetry in the F
_{1}-ATPase and its implications for the rotational cycle. Biophys. J. 86, 1373–1384 (2004). [PubMed] - Pänke O., Cherepanov D. A., Gumbiowski K., Engelbrecht S. & Junge W. Viscoelastic dynamics of actin filaments coupled to rotary f-ATPase: Torque profile of the enzyme. Biophys. J. 81, 1220–1233 (2001). [PubMed]
- Furuike S. et al. . Resolving stepping rotation in thermus thermophilus ATP synthase with an essentially drag-free probe. Nature Comm. 6, 233(1–9) (2011). [PMC free article] [PubMed]
- Sakaki N.
et al. . One rotary mechanism for F
_{1}-ATPase over ATP concentrations from millimolar down to nanomolar. Biophys. J. 88, 2047–2056 (2005). [PubMed] - Bilyard T.
et al. . High-resolution single-molecule characterization of the enzymatic states in e-coli F
_{1}-ATPase. Philos. Trans. R. Soc. Lond. B 368, 20120023 (2013). [PMC free article] [PubMed] - Neukirch S., Goriely A. & Hausrath A. C. Chirality of coiled coils: elasticity matters. Phys. Rev. Lett. 100, 038105 (2008). [PubMed]
- Neukirch S., Goriely A. & Hausrath A. C. Elastic coiled-coils act as energy buffers in the ATP synthase. Intl. J. Non-Linear Mech. 43, 1064–1073 (2008).
- Suezaki Y. & Go N. Fluctuations and mechanical strength of α-helices of polyglycine and poly (L-alanine). Biopolymers 15, 2137–2153 (1976). [PubMed]
- Choe S. & Sun S. X. The elasticity of alpha-helices. J. Chem. Phys. 122, 244912 (2005). [PubMed]
- Landau L. D. & Lifshitz E. M.
Theory of Elasticity, vol. 7 of
*A Course of Theoretical Physics*(Pergamon Press, Oxford, 1970). - Evans K. E., Nkansah M. A., Hutchinson I. J. & Rogers S. C. Molecular network design. Nature 353, 124 (1991).
- Lakes R. S. Foam structures with a negative Poisson’s ratio. Science 235, 1038 (1987). [PubMed]
- Alberts B., Bray D. & Lewis J. Essential Cell Biology (Garland Science, New York, 2013), 4 edn.
- Gaspari Z. & Nyitray L. Coiled coils as possible models of protein structure evolution. Bio. Mol. Concepts 2, 199–210 (2011). [PubMed]
- Crick F. H. C.
The packing of
*α*-helices: Simple coiled coils. Acta Crystallogr. 6, 689–697 (1953). - Cohen C. & Parry D. A. D.
*α*-helical coiled coils and bundles: how to design an*α*-helical protein. Proteins 7, 1–15 (1990). [PubMed] - Wolgemuth C. W. & Sun S. X.
Elasticity of
*α*-helical coiled coils. Phys. Rev. Lett. 97, 248101 (2006). [PubMed] - Yogurtcua O. N., Wolgemuth C. W. & Sun S. X.
Mechanical response and conformational amplification in
*α*-helical coiled coils. Biophys. J. 99, 3895–3904 (2010). [PubMed] - Yamakawa H. Helical wormlike chains in polymer solutions (Springer, Heidelberg, 1997).
- Hawkins R. J. & McLeish T. C. B.
Dynamic allostery of protein
*α*-helical coiled-coils. J. R. Soc. Interface 3, 125–138 (2006). [PMC free article] [PubMed] - Liu J. et al. . A seven-helix coiled coil. Proc. Natl. Acad. Sci. USA 103, 15457–15462 (2006). [PubMed]
- Blundell J. R. & Terentjev E. M. Stretching semiflexible filaments and their networks. Macromolecules 42, 5388–5394 (2009).
- Nelson P. Biological Physics (W. H. Freeman, New York, 2007).
- Feynman R. P., Leighton R. B. & Sands M. The Feynman Lectures on Physics, vol. 1 (Addison-Wesley, 1963).
- Beyenbach K. W. & Wieczorek H. The V-type H+ ATPase: molecular structure and function, physiological roles and regulation. J. Exp. Biol. 209, 577–589 (2006). [PubMed]
- Murata T., Yamato I., Kakinuma Y., Leslie A. G. W. & Walker J. E. Structure of the rotor of the V-type Na+ ATPase from enterococcus hirae. Science 308, 654–659 (2005). [PubMed]

Articles from Scientific Reports are provided here courtesy of **Nature Publishing Group**

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. |