Dynamics of the bacterial flagellar motor: the effects of stator compliance, back steps, temperature, and rotational asymmetry

The rotation of a bacterial flagellar motor (BFM) is driven by multiple stators tethered to the cell wall. Here, we extend a recently proposed power-stroke model to study the BFM dynamics under different biophysical conditions. Our model explains several key experimental observations and reveals their underlying mechanisms. 1), The observed independence of the speed at low load on the number of stators is explained by a force-dependent stepping mechanism that is independent of the strength of the stator tethering spring. Conversely, without force-dependent stepping, an unrealistically weak stator spring is required. 2), Our model with back-stepping naturally explains the observed absence of a barrier to backward rotation. Using the same set of parameters, it also explains BFM behaviors in the high-speed negative-torque regime. 3), From the measured temperature dependence of the maximum speed, our model shows that stator-stepping is a thermally activated process with an energy barrier. 4), The recently observed asymmetry in the torque-speed curve between counterclockwise- and clockwise-rotating BFMs can be quantitatively explained by the asymmetry in the stator-rotor interaction potentials, i.e., a quasilinear form for the counterclockwise motor and a quadratic form for the clockwise motor.     

How Cytoplasmic Dynein Couples ATP Hydrolysis Cycle to Diverse Stepping Motions: Kinetic Modeling

Cytoplasmic dynein is a two-headed molecular motor that moves to the minus end of a microtubule by ATP hydrolysis free energy. By employing its two heads (motor domains), cytoplasmic dynein exhibits various bipedal stepping motions: inchworm and hand-over-hand motions, as well as nonalternating steps of one head. However, the molecular basis to achieve such diverse stepping manners remains unclear because of the lack of an experimental method to observe stepping and the ATPase reaction of dynein simultaneously. Here, we propose a kinetic model for bipedal motions of cytoplasmic dynein and perform Gillespie Monte Carlo simulations that qualitatively reproduce most experimental data obtained to date. The model represents the status of each motor domain as five states according to conformation and nucleotide- and microtubule-binding conditions of the domain. In addition, the relative positions of the two domains were approximated by three discrete states. Accompanied by ATP hydrolysis cycles, the model dynein stochastically and processively moved forward in multiple steps via diverse pathways, including inchworm and hand-over-hand motions, similarly to experimental data. The model reproduced key experimental motility-related properties, including velocity and run length, as functions of the ATP concentration and external force, therefore providing a plausible explanation of how dynein achieves various stepping manners with explicit characterization of nucleotide states. Our model highlights the uniqueness of dynein in the coupling of ATPase with its movement during both inchworm and hand-over-hand stepping.     

The role of thermal activation in motion and force generation by molecular motors

The currently accepted mechanism for ATP-driven motion of kinesin is called the hand-over-hand model, where some chemical transition during the ATP hydrolysis cycle stretches a spring, and motion and force production result from the subsequent relaxation. It is essential in this mechanism for the moving head of kinesin to dissociate, while the other head remains firmly attached to the microtubule. Here we propose an alternative Brownian motor model where the action of ATP modulates the interaction potential between kinesin and the microtubule rather than a spring internal to the kinesin molecule alone. In this model neither head need dissociate (which predicts that under some circumstances a single-headed kinesin can display processive motion) and the transitions by which the motor moves are best described as thermally activated steps. This model is consistent with a wide range of experimental data on the force-velocity curves, the one ATP to one-step stoichiometry observed at small load, and the stochastic properties of the stepping.     

Time-integrated position error accounts for sensorimotor behavior in time-constrained tasks

Several studies have shown that human motor behavior can be successfully described using optimal control theory, which describes behavior by optimizing the trade-off between the subject's effort and performance. This approach predicts that subjects reach the goal exactly at the final time. However, another strategy might be that subjects try to reach the target position well before the final time to avoid the risk of missing the target. To test this, we have investigated whether minimizing the control effort and maximizing the performance is sufficient to describe human motor behavior in time-constrained motor tasks. In addition to the standard model, we postulate a new model which includes an additional cost criterion which penalizes deviations between the position of the effector and the target throughout the trial, forcing arrival on target before the final time. To investigate which model gives the best fit to the data and to see whether that model is generic, we tested both models in two different tasks where subjects used a joystick to steer a ball on a screen to hit a target (first task) or one of two targets (second task) before a final time. Noise of different amplitudes was superimposed on the ball position to investigate the ability of the models to predict motor behavior for different levels of uncertainty. The results show that a cost function representing only a trade-off between effort and accuracy at the end time is insufficient to describe the observed behavior. The new model correctly predicts that subjects steer the ball to the target position well before the final time is reached, which is in agreement with the observed behavior. This result is consistent for all noise amplitudes and for both tasks.     

Productive Cooperation among Processive Motors Depends Inversely on Their Mechanochemical Efficiency

Subcellular cargos are often transported by teams of processive molecular motors, which raises questions regarding the role of motor cooperation in intracellular transport. Although our ability to characterize the transport behaviors of multiple-motor systems has improved substantially, many aspects of multiple-motor dynamics are poorly understood. This work describes a transition rate model that predicts the load-dependent transport behaviors of multiple-motor complexes from detailed measurements of a single motor's elastic and mechanochemical properties. Transition rates are parameterized via analyses of single-motor stepping behaviors, load-rate-dependent motor-filament detachment kinetics, and strain-induced stiffening of motor-cargo linkages. The model reproduces key signatures found in optical trapping studies of structurally defined complexes composed of two kinesin motors, and predicts that multiple kinesins generally have difficulties in cooperating together. Although such behavior is influenced by the spatiotemporal dependence of the applied load, it appears to be directly linked to the efficiency of kinesin's stepping mechanism, and other types of less efficient and weaker processive motors are predicted to cooperate more productively. Thus, the mechanochemical efficiencies of different motor types may determine how effectively they cooperate together, and hence how motor copy number contributes to the regulation of cargo motion.     

A model for the coordinated stepping of cytoplasmic dynein

Cytoplasmic dynein play an important role in transporting various intracellular cargos by coupling their ATP hydrolysis cycle with their conformational changes. Recent experimental results showed that the cytoplasmic dynein had a highly variable stepping pattern including “hand-over-hand”, “inchworm” and “nonalternating-inchworm”. Here, we developed a model to describe the coordinated stepping patterns of cytoplasmic dynein, based on its working cycle, construction and the interaction between its leading head and tailing head. The kinetic model showed how change in the distance between the two heads influences the rate of cytoplasmic dynein under different stepping patterns. Numerical simulations of the distribution of step size and striding rate are in good quantitative agreement with experimental observations. Hence, our coordinated stepping model for cytoplasmic dynein successfully explained its diverse stepping patterns as a molecular motor. The cooperative mechanism carried out by the two heads of cytoplasmic dynein shed light on the strategies adopted by the cytoplasmic dynein in executing various functions.     

Torque generation by the Fo motor of the sodium ATPase

Based on recent structural and functional findings, we have constructed a mathematical model for the sodium-driven Fo motor of the F1Fo-ATPase from the anaerobic bacterium Propionigenium modestum. The model reveals the mechanochemical principles underlying the Fo motor's operation, and explains all of the existing experimental data on wild-type and mutant Fo motors. In particular, the model predicts a nonmonotonic dependence of the ATP hydrolysis activity on the sodium concentration, a prediction confirmed by new experiments. To explain experimental observations, the positively charged stator residue (R227) must assume different positions in the ATP synthesis and hydrolysis directions. This work also illustrates how to extract a motor mechanism from dynamical experimental observations in the absence of complete structural information.     

Direct observation of speed fluctuations of flagellar motor rotation at extremely low load close to zero

The bacterial flagellar motor accommodates ten stator units around the rotor to produce large torque at high load. But when external load is low, some previous studies showed that a single stator unit can spin the rotor at the maximum speed, suggesting that the maximum speed does not depend on the number of active stator units, whereas others reported that the speed is also dependent on the stator number. To clarify these two controversial observations, much more precise measurements of motor rotation would be required at external load as close to zero as possible. Here, we constructed a Salmonella filament-less mutant that produces a rigid, straight, twice longer hook to efficiently label a 60 nm gold particle and analyzed flagellar motor dynamics at low load close to zero. The maximum motor speed was about 400 Hz. Large speed fluctuations and long pausing events were frequently observed, and they were suppressed by either over-expression of the MotAB stator complex or increase in the external load, suggesting that the number of active stator units in the motor largely fluctuates near zero load. We conclude that the lifetime of the active stator unit becomes much shorter when the motor operates near zero load.     

Motor Dynamics Underlying Cargo Transport by Pairs of Kinesin-1 and Kinesin-3 Motors

Intracellular cargo transport by kinesin family motor proteins is crucial for many cellular processes, particularly vesicle transport in axons and dendrites. In a number of cases, the transport of specific cargo is carried out by two classes of kinesins that move at different speeds and thus compete during transport. Despite advances in single-molecule characterization and modeling approaches, many questions remain regarding the effect of intermotor tension on motor attachment/reattachment rates during cooperative multimotor transport. To understand the motor dynamics underlying multimotor transport, we analyzed the complexes of kinesin-1 and kinesin-3 motors attached through protein scaffolds moving on immobilized microtubules in vitro. To interpret the observed behavior, simulations were carried out using a model that incorporated motor stepping, attachment/detachment rates, and intermotor force generation. In single-molecule experiments, isolated kinesin-3 motors moved twofold faster and had threefold higher landing rates than kinesin-1. When the positively charged loop 12 of kinesin-3 was swapped with that of kinesin-1, the landing rates reversed, indicating that this “K-loop” is a key determinant of the motor reattachment rate. In contrast, swapping loop 12 had negligible effects on motor velocities. Two-motor complexes containing one kinesin-1 and one kinesin-3 moved at different speeds depending on the identity of their loop 12, indicating the importance of the motor reattachment rate on the cotransport speed. Simulations of these loop-swapped motors using experimentally derived motor parameters were able to reproduce the experimental results and identify best fit parameters for the motor reattachment rates for this geometry. Simulation results also supported previous work, suggesting that kinesin-3 microtubule detachment is very sensitive to load. Overall, the simulations demonstrate that the transport behavior of cargo carried by pairs of kinesin-1 and -3 motors are determined by three properties that differ between these two families: the unloaded velocity, the load dependence of detachment, and the motor reattachment rate.     

Bacterial Flagellar Motor Switch in Response to CheY-P Regulation and Motor Structural Alterations

The bacterial flagellar motor (BFM) is a molecular machine that rotates the helical filaments and propels the bacteria swimming toward favorable conditions. In our previous works, we built a stochastic conformational spread model to explain the dynamic and cooperative behavior of BFM switching. Here, we extended this model to test whether it can explain the latest experimental observations regarding CheY-P regulation and motor structural adaptivity. We show that our model predicts a strong correlation between rotational direction and the number of CheY-Ps bound to the switch complex, in agreement with the latest finding from Fukuoka et al. It also predicts that the switching sensitivity of the BFM can be fine-tuned by incorporating additional units into the switch complex, as recently demonstrated by Yuan et al., who showed that stoichiometry of FliM undergoes dynamic change to maintain ultrasensitivity in the motor switching response. In addition, by locking some rotor switching units on the switch complex into the stable clockwise-only conformation, our model has accurately simulated recent experiments expressing clockwise-locked FliG(ΔPAA) into the switch complex and reproduced the increased switching rate of the motor.     

Force-induced bidirectional stepping of cytoplasmic dynein

Cytoplasmic dynein is a minus-end-directed microtubule motor whose mechanism of movement remains poorly understood. Here, we use optical tweezers to examine the force-dependent stepping behavior of yeast cytoplasmic dynein. We find that dynein primarily advances in 8 nm increments but takes other sized steps (4-24 nm) as well. An opposing force induces more frequent backward stepping by dynein, and the motor walks backward toward the microtubule plus end at loads above its stall force of 7 pN. Remarkably, in the absence of ATP, dynein steps processively along microtubules under an external load, with less force required for minus-end- than for plus-end-directed movement. This nucleotide-independent walking reveals that force alone can drive repetitive microtubule detachment-attachment cycles of dynein's motor domains. These results suggest a model for how dynein's two motor domains coordinate their activities during normal processive motility and provide new clues for understanding dynein-based motility in living cells.     

Helix rotation model of the flagellar rotary motor

A new model of the flagellar motor is proposed that is based on established dynamics of the KcsA potassium ion channel and on known genetic, biochemical, and biophysical facts, which accounts for the mechanics of torque generation, force transmission, and reversals of motor rotation. It predicts that proton (or in some species sodium ion) flow generates short, reversible helix rotations of the MotA-MotB channel complex (the stator) that are transmitted by Coulomb forces to the FliG segments at the rotor surface. Channels are arranged as symmetric pairs, S and T, that swing back and forth in synchrony. S and T alternate in attaching to the rotor, so that force transmission proceeds in steps. The sense of motor rotation can be readily reversed by conformationally switching the position of charged groups on the rotor so that they interact with the stator during the reverse rather than forward strokes. An elastic device accounts for the observed smoothness of rotation and a prolonged attachment of the torque generators to the rotor, i.e., a high duty ratio of each torque-generating unit.     

Steps in the Bacterial Flagellar Motor

The bacterial flagellar motor is a highly efficient rotary machine used by many bacteria to propel themselves. It has recently been shown that at low speeds its rotation proceeds in steps. Here we propose a simple physical model, based on the storage of energy in protein springs, that accounts for this stepping behavior as a random walk in a tilted corrugated potential that combines torque and contact forces. We argue that the absolute angular position of the rotor is crucial for understanding step properties and show this hypothesis to be consistent with the available data, in particular the observation that backward steps are smaller on average than forward steps. We also predict a sublinear speed versus torque relationship for fixed load at low torque, and a peak in rotor diffusion as a function of torque. Our model provides a comprehensive framework for understanding and analyzing stepping behavior in the bacterial flagellar motor and proposes novel, testable predictions. More broadly, the storage of energy in protein springs by the flagellar motor may provide useful general insights into the design of highly efficient molecular machines.     

A kinetic and stochastic analysis of crossbridge-type stepping mechanisms in rotary molecular motors

The bacterial flagellar motor is generally supposed to be a stepping mechanism. The main evidence for this is based on a fluctuation analysis of experiments with tethered bacteria in which rotation frequency was varied by applying an external torque: the variance in time taken for a fixed number of revolutions was found to be essentially proportional to the inverse square of the frequency. This behavior was shown to characterize a Poissonian stepper. Here we present a rigorous kinetic and stochastic analysis of elastic crossbridge stepping in tethered bacteria. We demonstrate that Poissonian stepping is a virtually unachievable limit. To the extent that a system may approach Poissonian stepping it cannot be influenced by an externally applied torque; stepping mechanisms capable of being so influenced are necessarily non-Poissonian and exhibit an approximately inverse cubic dependence. This conclusion applies whatever the torsional characteristics of the tether may be, and contrary to claims, no perceptible relaxation of the tether following each step is found. Furthermore, the inverse square dependence is a necessary but not sufficient condition for Poissonian stepping, since a nonstepping mechanism, which closely reproduces most experimental data, also fulfills this condition. Hence the inference that crossbridge-type stepping occurs is not justified.     

Evaluation of the Duty Ratio of the Bacterial Flagellar Motor by Dynamic Load Control

The bacterial flagellar motor is one of the most complex and sophisticated nanomachineries in nature. A duty ratio D is a fraction of time that the stator and the rotor interact and is a fundamental property to characterize the motor but remains to be determined. It is known that the stator units of the motor bind to and dissociate from the motor dynamically to control the motor torque depending on the load on the motor. At low load, at which the kinetics such as proton translocation speed limits the rotation rate, the dependency of the rotation rate on the number of stator units N implies D: the dependency becomes larger for smaller D. Contradicting observations supporting both the small and large D have been reported. A dilemma is that it is difficult to explore a broad range of N at low load because the stator units easily dissociate, and N is limited to one or two at vanishing load. Here, we develop an electrorotation method to dynamically control the load on the flagellar motor of Salmonella with a calibrated magnitude of the torque. By instantly reducing the load for keeping N high, we observed that the speed at low load depends on N, implying a small duty ratio. We recovered the torque-speed curves of individual motors and evaluated the duty ratio to be 0.14 ± 0.04 from the correlation between the torque at high load and the rotation rate at low load.     

Biomechanical influences on balance recovery by stepping.

Stepping represents a common means for balance recovery after a perturbation to upright posture. Yet little is known regarding the biomechanical factors which determine whether a step succeeds in preventing a fall. In the present study, we developed a simple pendulum-spring model of balance recovery by stepping, and used this to assess how step length and step contact time influence the effort (leg contact force) and feasibility of balance recovery by stepping. We then compared model predictions of step characteristics which minimize leg contact force to experimentally observed values over a range of perturbation strengths. At all perturbation levels, experimentally observed step execution times were higher than optimal, and step lengths were smaller than optimal. However, the predicted increase in leg contact force associated with these deviations was substantial only for large perturbations. Furthermore, increases in the strength of the perturbation caused subjects to take larger, quicker steps, which reduced their predicted leg contact force. We interpret these data to reflect young subjects' desire to minimize recovery effort, subject to neuromuscular constraints on step execution time and step length. Finally, our model predicts that successful balance recovery by stepping is governed by a coupling between step length, step execution time, and leg strength, so that the feasibility of balance recovery decreases unless declines in one capacity are offset by enhancements in the others. This suggests that one's risk for falls may be affected more by small but diffuse neuromuscular impairments than by larger impairment in a single motor capacity.     

Collective dynamics of elastically coupled Myosin v motors

Characterization of the collective behaviors of different classes of processive motor proteins has become increasingly important to understand various intracellular trafficking and transport processes. This work examines the dynamics of structurally-defined motor complexes containing two myosin Va (myoVa) motors that are linked together via a molecular scaffold formed from a single duplex of DNA. Dynamic changes in the filament-bound configuration of these complexes due to motor binding, stepping, and detachment were monitored by tracking the positions of different color quantum dots that report the position of one head of each myoVa motor on actin. As in studies of multiple kinesins, the run lengths produced by two myosins are only slightly larger than those of single motor molecules. This suggests that internal strain within the complexes, due to asynchronous motor stepping and the resultant stretching of motor linkages, yields net negative cooperative behaviors. In contrast to multiple kinesins, multiple myosin complexes move with appreciably lower velocities than a single-myosin molecule. Although similar trends are predicted by a discrete state stochastic model of collective motor dynamics, these analyses also suggest that multiple myosin velocities and run lengths depend on both the compliance and the effective size of their cargo. Moreover, it is proposed that this unique collective behavior occurs because the large step size and relatively small stalling force of myoVa leads to a high sensitivity of motor stepping rates to strain.     

The Winch Model Can Explain both Coordinated and Uncoordinated Stepping of Cytoplasmic Dynein

Cytoplasmic dynein moves processively along microtubules, but the mechanism of how its heads use the energy from ATP hydrolysis, coupled to a linker swing, to achieve directed motion, is still unclear. In this article, we present a theoretical model based on the winch mechanism in which the principal direction of the linker stroke is toward the microtubule-binding domain. When mechanically coupling two identical heads (each with postulated elastic properties and a minimal ATPase cycle), the model reproduces stepping with 8-nm steps (even though the motor itself is much larger), interhead coordination, and processivity, as reported for mammalian dyneins. Furthermore, when we loosen the elastic connection between the heads, the model still shows processive directional stepping, but it becomes uncoordinated and the stepping pattern shows a greater variability, which reproduces the properties of yeast dyneins. Their slower chemical kinetics allows processive motility and a high stall force without the need for coordination.