A simple kinetic model describes the processivity of myosin-v

Dynamics of ATP-dependent and ATP-independent steppings of myosin-V on actin: catch-bond characteristics

An analytical theory is presented for the dynamics of myosin-V molecular motor, where both ATP-dependent and ATP-independent steppings are taken into account. Specifically, the dependences of velocity, run length and unbinding rate upon both forward and backward loads and ATP concentration are studied, explaining quantitatively the diverse available single-molecule data and providing predicted results. The results show that the unbinding rate increases with the increase of ATP concentration and levels off at both low and high ATP concentrations. More interestingly, at an ATP concentration that is not very low, the unbinding rate exhibits characteristics of a catch-slip bond under backward load, with the unbinding rate decreasing rapidly with the increase of the backward load in the range smaller than about 2.5 pN and then increasing slowly with the further increase of the backward load. By contrast, under forward load the unbinding rate exhibits a slip-bond characteristic.

How kinesin waits for ATP affects the nucleotide and load dependence of the stepping kinetics

Conventional kinesin, responsible for directional transport of cellular vesicles, takes multiple nearly uniform 8.2-nm steps by consuming one ATP molecule per step as it walks toward the plus end of the microtubule (MT). Despite decades of intensive experimental and theoretical studies, there are gaps in the elucidation of key steps in the catalytic cycle of kinesin. How the motor waits for ATP to bind to the leading head is controversial. Two experiments using a similar protocol have arrived at different conclusions. One asserts that kinesin waits for ATP in a state with both the heads bound to the MT, whereas the other shows that ATP binds to the leading head after the trailing head detaches. To discriminate between the 2 scenarios, we developed a minimal model, which analytically predicts the outcomes of a number of experimental observable quantities (the distribution of run length, the distribution of velocity [[Formula: see text]], and the randomness parameter) as a function of an external resistive force (F) and ATP concentration ([T]). The differences in the predicted bimodality in [Formula: see text] as a function of F between the 2 models may be amenable to experimental testing. Most importantly, we predict that the F and [T] dependence of the randomness parameters differ qualitatively depending on the waiting states. The randomness parameters as a function of F and [T] can be quantitatively measured from stepping trajectories with very little prejudice in data analysis. Therefore, an accurate measurement of the randomness parameter and the velocity distribution as a function of load and nucleotide concentration could resolve the apparent controversy.

Characteristics of Myosin V Processivity

Myosin V is a processive doubled-headed biomolecular motor involved in many intracellular organelle and vesicle transport. The unidirectional movement is coupled with the adenosine triphosphate (ATP) hydrolysis and product release cycle. With the progress of experimental techniques and the enhancement of measuring directness, detailed knowledge of the motility of myosin V has been obtained. Following the ATPase cycle, the 4-state mechanochemical model of the myosin V's processive movement is used. The transitions between various states take place in a stochastic manner. We can use the master equation to analyze and calculate quantitatively. Meanwhile, the effect of the reverse reaction is taken fully into account. We fit the mean velocity, the mean dwell time, the mean run length, and the ratio of forward/backward steps as a functionof ATP, ADP, and Pi concertration. The theoretical curves are generally in line with the experimental data. This work provides a new insight for the characteristic of myosin V.

A model for the chemomechanical coupling of myosin-V molecular motors

Herein, a model for the chemomechanical coupling of dimeric myosin-V motors is presented. Based on this model and the proposal that the rate constants of the ATPase activity of the two heads are independent of an external force in a range smaller than the stall force, we analytically studied the dynamics of the motor, such as the stepping ratio, dwell time between two mechanical steps, and velocity, under varying force and ATP concentrations. The theoretical results well reproduce the diverse available single-molecule experimental data. In particular, the experimental data showing that at a low ATP concentration, the dwell time and velocity have less force dependency than at a high ATP concentration is explained quantitatively. Moreover, the dependency of the chemomechanical coupling ratio on the force and ATP concentration was studied.

The paper presents a model of chemomechanical coupling of myosin-V motor, explaining the dynamics under varying force and ATP concentrations.

A Unified Walking Model for Dimeric Motor Proteins

Dimeric motor proteins, kinesin-1, cytoplasmic dynein-1, and myosin-V, move stepwise along microtubules and actin filaments with a regular step size. The motors take backward as well as forward steps. The step ratio r and dwell time τ, which are the ratio of the number of backward steps to the number of forward steps and the time between consecutive steps, respectively, were observed to change with the load. To understand the movement of motor proteins, we constructed a unified and simple mathematical model to explain the load dependencies of r and of τ measured for the above three types of motors quantitatively. Our model consists of three states, and the forward and backward steps are represented by the cycles of transitions visiting different pairs of states among the three, implying that a backward step is not the reversal of a forward step. Each of r and τ is given by a simple expression containing two exponential functions. The experimental data for r and τ for dynein available in the literature are not sufficient for a quantitative analysis, which is in contrast to those for kinesin and myosin-V. We reanalyze the data to obtain r and τ of native dynein to make up the insufficient data to fit them to the model. Our model successfully describes the behavior of r and τ for all of the motors in a wide range of loads from large assisting loads to superstall loads.

Energetic Costs, Precision, and Transport Efficiency of Molecular Motors

An efficient molecular motor would deliver cargo to the target site at a high speed and in a punctual manner while consuming a minimal amount of energy. According to a recently formulated thermodynamic principle, referred to as the thermodynamic uncertainty relation, the travel distance of a motor and its variance are, however, constrained by the free energy being consumed. Here we use the principle underlying the uncertainty relation to quantify the transport efficiency of molecular motors for varying ATP concentration ([ATP]) and applied load (f). Our analyses of experimental data find that transport efficiencies of the motors studied here are semioptimized under the cellular condition. The efficiency is significantly deteriorated for a kinesin-1 mutant that has a longer neck-linker, which underscores the importance of molecular structure. It is remarkable to recognize that, among many possible directions for optimization, biological motors have evolved to optimize the transport efficiency in particular.

Allocating and Splitting Free Energy to Maximize Molecular Machine Flux

Biomolecular machines transduce between different forms of energy. These machines make directed progress and increase their speed by consuming free energy, typically in the form of nonequilibrium chemical concentrations. Machine dynamics are often modeled by transitions between a set of discrete metastable conformational states. In general, the free-energy change associated with each transition can increase the forward rate constant, decrease the reverse rate constant, or both. In contrast to previous optimizations, we find that in general flux is maximized neither by devoting all free-energy changes to increasing forward rate constants nor by solely decreasing reverse rate constants. Instead, the optimal free-energy splitting depends on the detailed dynamics. Extending our analysis to machines with vulnerable states (from which they can break down), in the strong driving corresponding to in vivo cellular conditions, processivity is maximized by reducing the occupation of the vulnerable state.

Allocating dissipation across a molecular machine cycle to maximize flux

Biomolecular machines consume free energy to break symmetry and make directed progress. Nonequilibrium ATP concentrations are the typical free energy source, with one cycle of a molecular machine consuming a certain number of ATP, providing a fixed free energy budget. Since evolution is expected to favor rapid-turnover machines that operate efficiently, we investigate how this free energy budget can be allocated to maximize flux. Unconstrained optimization eliminates intermediate metastable states, indicating that flux is enhanced in molecular machines with fewer states. When maintaining a set number of states, we show that-in contrast to previous findings-the flux-maximizing allocation of dissipation is not even. This result is consistent with the coexistence of both "irreversible" and reversible transitions in molecular machine models that successfully describe experimental data, which suggests that, in evolved machines, different transitions differ significantly in their dissipation.

Cargo Transport by Two Coupled Myosin Va Motors on Actin Filaments and Bundles

Myosin Va (myoVa) is a processive, actin-based molecular motor essential for intracellular cargo transport. When a cargo is transported by an ensemble of myoVa motors, each motor faces significant physical barriers and directional challenges created by the complex actin cytoskeleton, a network of actin filaments and actin bundles. The principles that govern the interaction of multiple motors attached to the same cargo are still poorly understood. To understand the mechanical interactions between multiple motors, we developed a simple in vitro model in which two individual myoVa motors labeled with different-colored Qdots are linked via a third Qdot that acts as a cargo. The velocity of this two-motor complex was reduced by 27% as compared to a single motor, whereas run length was increased by only 37%, much less than expected from multimotor transport models. Therefore, at low ATP, which allowed us to identify individual motor steps, we investigated the intermotor dynamics within the two-motor complex. The randomness of stepping leads to a buildup of tension in the linkage between motors-which in turn slows down the leading motor-and increases the frequency of backward steps and the detachment rate. We establish a direct relationship between the velocity reduction and the distribution of intermotor distances. The analysis of run lengths and dwell times for the two-motor complex, which has only one motor engaged with the actin track, reveals that half of the runs are terminated by almost simultaneous detachment of both motors. This finding challenges the assumptions of conventional multimotor models based on consecutive motor detachment. Similar, but even more drastic, results were observed with two-motor complexes on actin bundles, which showed a run length that was even shorter than that of a single motor.

Tight-binding approach to overdamped Brownian motion on a bichromatic periodic potential

We present a theoretical treatment of overdamped Brownian motion on a time-independent bichromatic periodic potential with spatially fast- and slow-changing components. In our approach, we generalize the Wannier basis commonly used in the analysis of periodic systems to define a basis of S states that are localized at local minima of the potential. We demonstrate that the S states are orthonormal and complete on the length scale of the periodicity of the fast-changing potential, and we use the S-state basis to transform the continuous Smoluchowski equation for the system to a discrete master equation describing hopping between local minima. We identify the parameter regime where the master equation description is valid and show that the interwell hopping rates are well approximated by Kramers' escape rate in the limit of deep potential minima. Finally, we use the master equation to explore the system dynamics and determine the drift and diffusion for the system.

Skewness and Kurtosis in Statistical Kinetics

We obtain lower and upper bounds on the skewness and kurtosis associated with the cycle completion time of unicyclic enzymatic reaction schemes. Analogous to a well-known lower bound on the randomness parameter, the lower bounds on skewness and kurtosis are related to the number of intermediate states in the underlying chemical reaction network. Our results demonstrate that evaluating these higher order moments with single molecule data can lead to information about the enzymatic scheme that is not contained in the randomness parameter.

Detachment, futile cycling, and nucleotide pocket collapse in myosin-V stepping

Myosin-V is a highly processive dimeric protein that walks with 36-nm steps along actin tracks, powered by coordinated adenosine triphosphate (ATP) hydrolysis reactions in the two myosin heads. No previous theoretical models of the myosin-V walk reproduce all the observed trends of velocity and run length with adenosine diphosphate (ADP), ATP and external forcing. In particular, a result that has eluded all theoretical studies based upon rigorous physical chemistry is that run length decreases with both increasing [ADP] and [ATP]. We systematically analyze which mechanisms in existing models reproduce which experimental trends and use this information to guide the development of models that can reproduce them all. We formulate models as reaction networks between distinct mechanochemical states with energetically determined transition rates. For each network architecture, we compare predictions for velocity and run length to a subset of experimentally measured values, and fit unknown parameters using a bespoke Monte Carlo simulated annealing optimization routine. Finally we determine which experimental trends are replicated by the best-fit model for each architecture. Only two models capture them all: one involving [ADP]-dependent mechanical detachment, and another including [ADP]-dependent futile cycling and nucleotide pocket collapse. Comparing model-predicted and experimentally observed kinetic transition rates favors the latter.

Thermal fluctuation statistics in a molecular motor described by a multidimensional master equation

We present a theoretical investigation of thermal fluctuation statistics in a molecular motor. Energy transfer in the motor is described using a multidimensional discrete master equation with nearest-neighbor hopping. In this theory, energy transfer leads to statistical correlations between thermal fluctuations in different degrees of freedom. For long times, the energy transfer is a multivariate diffusion process with constant drift and diffusion. The fluctuations and drift align in the strong-coupling limit enabling a one-dimensional description along the coupled coordinate. We derive formal expressions for the probability distribution and simulate single trajectories of the system in the near- and far-from-equilibrium limits both for strong and weak coupling. Our results show that the hopping statistics provide an opportunity to distinguish different operating regimes.

Motor proteins and molecular motors: how to operate machines at the nanoscale

Several classes of biological molecules that transform chemical energy into mechanical work are known as motor proteins or molecular motors. These nanometer-sized machines operate in noisy stochastic isothermal environments, strongly supporting fundamental cellular processes such as the transfer of genetic information, transport, organization and functioning. In the past two decades motor proteins have become a subject of intense research efforts, aimed at uncovering the fundamental principles and mechanisms of molecular motor dynamics. In this review, we critically discuss recent progress in experimental and theoretical studies on motor proteins. Our focus is on analyzing fundamental concepts and ideas that have been utilized to explain the non-equilibrium nature and mechanisms of molecular motors.

Extracting signal from noise: kinetic mechanisms from a Michaelis-Menten-like expression for enzymatic fluctuations

Enzyme-catalyzed reactions are naturally stochastic, and precision measurements of these fluctuations, made possible by single-molecule methods, promise to provide fundamentally new constraints on the possible mechanisms underlying these reactions. We review some aspects of statistical kinetics: a new field with the goal of extracting mechanistic information from statistical measures of fluctuations in chemical reactions. We focus on a widespread and important statistical measure known as the randomness parameter. This parameter is remarkably simple in that it is the squared coefficient of variation of the cycle completion times, although it places significant limits on the minimal complexity of possible enzymatic mechanisms. Recently, a general expression has been introduced for the substrate dependence of the randomness parameter that is for rate fluctuations what the Michaelis-Menten expression is for the mean rate of product generation. We discuss the information provided by the new kinetic parameters introduced by this expression and demonstrate that this expression can simplify the vast majority of published models.

Design principles governing the motility of myosin V

The molecular motor myosin V (MyoV) exhibits a wide repertoire of pathways during the stepping process, which is intimately connected to its biological function. The best understood of these is the hand-over-hand stepping by a swinging lever arm movement toward the plus end of actin filaments. Single-molecule experiments have also shown that the motor "foot stomps," with one hand detaching and rebinding to the same site, and back-steps under sufficient load. The complete taxonomy of MyoV's load-dependent stepping pathways, and the extent to which these are constrained by motor structure and mechanochemistry, are not understood. Using a polymer model, we develop an analytical theory to describe the minimal physical properties that govern motor dynamics. We solve the first-passage problem of the head reaching the target-binding site, investigating the competing effects of backward load, strain in the leading head biasing the diffusion in the direction of the target, and the possibility of preferential binding to the forward site due to the recovery stroke. The theory reproduces a variety of experimental data, including the power stroke and slow diffusive search regimes in the mean trajectory of the detached head, and the force dependence of the forward-to-backward step ratio, run length, and velocity. We derive a stall force formula, determined by lever arm compliance and chemical cycle rates. By exploring the MyoV design space, we predict that it is a robust motor whose dynamical behavior is not compromised by reasonable perturbations to the reaction cycle and changes in the architecture of the lever arm.

The load dependence of the physical properties of a molecular motor

The physical properties of a molecular motor with load changing in a wide range will be discussed in this study, in particular the mean velocity, output power and energy efficiency. The main difficulty of this study is that both the states of the molecular motor and the energy barriers between them change with the loading force. Moreover, with the change of load, the number of motor states may also change, so different models should be used to calculate the corresponding physical quantities in different ranges of load. The results show that, in contrast to the usual intuition, the mean velocity and output power of the molecular motor do not change continuously with load.

Bounds and phase diagram of efficiency at maximum power for tight-coupling molecular motors

The efficiency at maximum power (EMP) for tight-coupling molecular motors is investigated within the framework of irreversible thermodynamics. It is found that the EMP depends merely on the constitutive relation between the thermodynamic current and force. The motors are classified into four generic types (linear, superlinear, sublinear, and mixed types) according to the characteristics of the constitutive relation, and then the corresponding ranges of the EMP for these four types of molecular motors are obtained. The exact bounds of the EMP are derived and expressed as the explicit functions of the free energy released by the fuel in each motor step. A phase diagram is constructed which clearly shows how the region where the parameters (the load distribution factor and the free energy released by the fuel in each motor step) are located can determine whether the value of the EMP is larger or smaller than 1/2. This phase diagram reveals that motors using ATP as fuel under physiological conditions can work at maximum power with higher efficiency (> 1/2) for a small load distribution factor (< 0.1).

A branched kinetic scheme describes the mechanochemical coupling of Myosin Va processivity in response to substrate

Myosin Va is a double-headed cargo-carrying molecular motor that moves processively along cellular actin filaments. Long processive runs are achieved through mechanical coordination between the two heads of myosin Va, which keeps their ATPase cycles out of phase, preventing both heads detaching from actin simultaneously. The biochemical kinetics underlying processivity are still uncertain. Here we attempt to define the biochemical pathways populated by myosin Va by examining the velocity, processive run-length, and individual steps of a Qdot-labeled myosin Va in various substrate conditions (i.e., changes in ATP, ADP, and P(i)) under zero load in the single-molecule total internal reflection fluorescence microscopy assay. These data were used to globally constrain a branched kinetic scheme that was necessary to fit the dependences of velocity and run-length on substrate conditions. Based on this model, myosin Va can be biased along a given pathway by changes in substrate concentrations. This has uncovered states not normally sampled by the motor, and suggests that every transition involving substrate binding and release may be strain-dependent.

Phenomenological analysis of ATP dependence of motor proteins

In this study, through phenomenological comparison of the velocity-force data of processive motor proteins, including conventional kinesin, cytoplasmic dynein and myosin V, I found that, the ratio between motor velocities of two different ATP concentrations is almost invariant for any substall, superstall or negative external loads. Therefore, the velocity of motors can be well approximated by a Michaelis-Menten like formula , with the step size, and the external load dependent rate of one mechanochemical cycle of motor motion in saturated ATP solution. The difference of Michaelis-Menten constant for substall, superstall and negative external load indicates, the configurations at which ATP molecule can bind to motor heads for these three cases might be different, though the expression of as a function of might be unchanged for any external load . Verifications of this Michaelis-Menten like formula has also been done by fitting to the recent experimental data.

Kinesins with extended neck linkers: a chemomechanical model for variable-length stepping

We develop a stochastic model for variable-length stepping of kinesins engineered with extended neck linkers. This requires that we consider the separation in microtubule binding sites between the heads of the motor at the beginning of a step. We show that this separation is stationary and can be included in the calculation of standard experimental quantities. We also develop a corresponding matrix computational framework for conducting computer experiments. Our matrix approach is more efficient computationally than large-scale Monte Carlo simulation. This efficiency greatly eases sensitivity analysis, an important feature when there is considerable uncertainty in the physical parameters of the system. We demonstrate the application and effectiveness of our approach by showing that the worm-like chain model for the neck linker can explain recently published experimental data. While we have focused on a particular scenario for kinesins, these methods could also be applied to myosin and other processive motors.

Chemomechanical coupling and motor cycles of myosin V

The molecular motor myosin V has been studied extensively both in bulk and single molecule experiments. Based on the chemical states of the motor, we construct a systematic network theory that includes experimental observations about the stepping behavior of myosin V. We utilize constraints arising from nonequilibrium thermodynamics to determine motor parameters and demonstrate that the motor behavior is governed by three chemomechanical motor cycles. The competition between these cycles can be understood via the influence of external load forces onto the chemical transition rates for the binding of adenosine triphosphate and adenosine diphosphate. In addition, we also investigate the functional dependence of the mechanical stepping rates on these forces. For substall forces, the dominant pathway of the motor network is profoundly different from the one for superstall forces, which leads to a stepping behavior that is in agreement with the experimental observations. Our theory provides a unified description of the experimental data as obtained for myosin V in single motor experiments.

Coarse-grained modeling of conformational transitions underlying the processive stepping of myosin V dimer along filamentous actin

To explore the structural basis of processive stepping of myosin V along filamentous actin, we have performed comprehensive modeling of its key conformational states and transitions with an unprecedented residue level of details. We have built structural models for a myosin V monomer complexed with filamentous actin at four biochemical states [adenosine diphosphate (ATP)-, adenosine diphosphate (ADP)-phosphate-, ADP-bound or nucleotide-free]. Then we have modeled a myosin V dimer (consisting of lead and rear head) at various two-head-bound states with nearly straight lever arms rotated by intramolecular strain. Next, we have performed transition pathway modeling to determine the most favorable sequence of transitions (namely, phosphate release at the lead head followed by ADP release at the rear head, while ADP release at the lead head is inhibited), which underlie the kinetic coordination between the two heads. Finally, we have used transition pathway modeling to reveal the order of structural changes during three key biochemical transitions (phosphate release at the lead head, ADP release and ATP binding at the rear head), which shed lights on the strain-dependence of the allosterically coupled motions at various stages of myosin V's work cycle. Our modeling results are in agreement with and offer structural insights to many results of kinetic, single-molecule and structural studies of myosin V.

Mechanistic constraints from the substrate concentration dependence of enzymatic fluctuations

The time it takes an enzyme to complete its reaction is a stochastic quantity governed by thermal fluctuations. With the advent of high-resolution methods of single-molecule manipulation and detection, it is now possible to observe directly this natural variation in the enzymatic cycle completion time and extract kinetic information from the statistics of its fluctuations. To this end, the inverse of the squared coefficient of variation, which we term n(min), is a useful measure of fluctuations because it places a strict lower limit on the number of kinetic states in the enzymatic mechanism. Here we show that there is a single general expression for the substrate dependence of n(min) for a wide range of kinetic models. This expression is governed by three kinetic parameters, which we term N(L), N(S), and alpha. These parameters have simple geometric interpretations and provide clear constraints on possible kinetic mechanisms. As a demonstration of this analysis, we fit the fluctuations in the dwell times of the packaging motor of the bacteriophage varphi29, revealing additional features of the nucleotide loading process in this motor. Because a diverse set of kinetic models display the same substrate dependence for their fluctuations, the expression for this general dependence may prove of use in the characterization and study of the dynamics of a wide range of enzymes.

Comprehensive physical mechanism of two-headed biomotor myosin V

Two-headed biomotor myosin V autonomously coordinates its two identical heads in fuel consumption and mechanical stepping, so that the dimerized motor as a whole gains the capability of processive, unidirectional movement along cytoskeletal filament. How the dimer-level functions like sustained direction rectification and autonomous coordination emerge out of physical principles poses an outstanding question pertinent to motor protein biology as well as the nascent field of bioinspired nanomotors. Here the comprehensive physical mechanism for myosin V motor is identified by a dimer-level free-energy analysis that is methodologically calibrated against experimental data. A hallmark of the identified mechanism is a mechanically mediated symmetry breaking that occurs at the dimer level and prevails against ubiquitous thermal fluctuations. Another character is the onset of substantial free-energy gaps between major dimer-track binding configurations. The symmetry breaking is the basis for myosin V's directional rectification, and the energy gaps facilitate autonomous head-head coordination. The mechanism explains the experimental finding that myosin V makes ATP-independent consecutive steps under high opposing loads but not under pushing loads. Interestingly, myosin V and another major biomotor kinesin 1 are found to share essentially the same core mechanism but for distinctly different working regimes.

Extensibility of the extended tail domain of processive and nonprocessive myosin V molecules

Myosin V is a single-molecule motor that moves organelles along actin. When myosin V pulls loads inside the cell in a highly viscous environment, the force on the motor is unlikely to be constant. We propose that the tether between the single-molecule motor and the cargo (i.e., the extended tail domain of the molecule) must be able to absorb the sudden mechanical motions of the motor and allow smooth relaxation of the motion of the cargo to a new position. To test this hypothesis, we compared the elastic properties of the extended tail domains of processive (mouse myosin Va) and nonprocessive (Drosophila myosin V) molecular motors. The extended tail domain of these myosins consists of mechanically strong coiled-coil regions interspersed with flexible loops. In this work we explored the mechanical properties of coiled-coil regions using atomic force microscopy. We found that the processive and nonprocessive coiled-coil fragments display different unfolding patterns. The unfolding of coiled-coil structures occurs much later during the atomic force microscopy stretch cycle for processive myosin Va than for nonprocessive Drosophila myosin V, suggesting that this elastic tether between the cargo and motor may play an important role in sustaining the processive motions of this single-molecule motor.

Biological proton pumping in an oscillating electric field

Time-dependent external perturbations provide powerful probes of the function of molecular machines. Here we study biological proton pumping in an oscillating electric field. The protein cytochrome c oxidase is the main energy transducer in aerobic life, converting chemical energy into an electric potential by pumping protons across a membrane. With the help of master-equation descriptions that recover the key thermodynamic and kinetic properties of this biological "fuel cell," we show that the proton pumping efficiency and the electronic currents in steady state depend significantly on the frequency and amplitude of the applied field, allowing us to distinguish between different microscopic mechanisms of the machine. A spectral analysis reveals dominant reaction steps consistent with an electron-gated pumping mechanism.

Kinetic gating of the proton pump in cytochrome c oxidase

Cytochrome c oxidase (CcO), the terminal enzyme of the respiratory chain, reduces oxygen to water and uses the released energy to pump protons across a membrane. Here, we use kinetic master equations to explore the energetic and kinetic control of proton pumping in CcO. We construct models consistent with thermodynamic principles, the structure of CcO, experimentally known proton affinities, and equilibrium constants of intermediate reactions. The resulting models are found to capture key properties of CcO, including the midpoint redox potentials of the metal centers and the electron transfer rates. We find that coarse-grained models with two proton sites and one electron site can pump one proton per electron against membrane potentials exceeding 100 mV. The high pumping efficiency of these models requires strong electrostatic couplings between the proton loading (pump) site and the electron site (heme a), and kinetic gating of the internal proton transfer. Gating is achieved by enhancing the rate of proton transfer from the conserved Glu-242 to the pump site on reduction of heme a, consistent with the predictions of the water-gated model of proton pumping. The model also accounts for the phenotype of D-channel mutations associated with loss of pumping but retained turnover. The fundamental mechanism identified here for the efficient conversion of chemical energy into an electrochemical potential should prove relevant also for other molecular machines and novel fuel-cell designs.

Enhancement of cargo processivity by cooperating molecular motors

Cellular cargo can be bound to cytoskeletal filaments by one or more active or passive molecular motors. Recent experiments have shown that the presence of auxiliary, nondriving motors results in an enhanced processivity of the cargo, compared to the case of a single active driving motor alone. We model the observed cooperative transport process using a stochastic model that describes the dynamics of two molecular motors, an active one that moves cargo unidirectionally along a filament track, and a passive one that acts as a tether. Analytical expressions obtained from our analysis are fit to experimental data to estimate the microscopic kinetic parameters of our model. Our analysis reveals two qualitatively distinct processivity-enhancing mechanisms: the passive tether can decrease the typical detachment rate of the active motor from the filament track or it can increase the corresponding reattachment rate. Comparing analytical results with experimental data, we can show unambiguously that in the case of kinesin transport on microtubules, a higher average run length arises mainly from the ability of the passive motor to keep the cargo close to the filament, enhancing the reattachment rate of recently detached active kinesin motors. On the other hand, in the case of myosin-driven transport along actin, the passive motor tightly tethers the cargo to the filament, suppressing the detachment rate of the active myosin.

The energetics of allosteric regulation of ADP release from myosin heads

Myosin molecules are involved in a wide range of transport and contractile activities in cells. A single myosin head functions through its ATPase reaction as a force generator and as a mechanosensor, and when two or more myosin heads work together in moving along an actin filament, the interplay between these mechanisms contributes to collective myosin behaviors. For example, the interplay between force-generating and force-sensing mechanisms coordinates the two heads of a myosin V molecule in its hand-over-hand processive stepping along an actin filament. In muscle, it contributes to the Fenn effect and smooth muscle latch. In both examples, a key force-sensing mechanism is the regulation of ADP release via interhead forces that are generated upon actin-myosin binding. Here we present a model describing the mechanism of allosteric regulation of ADP release from myosin heads as a change, DeltaDeltaG(-D), in the standard free energy for ADP release that results from the work, Deltamicro(mech), performed by that myosin head upon ADP release, or DeltaDeltaG(-D) = Deltamicro(mech). We show that this model is consistent with previous measurements for strain-dependent kinetics of ADP release in both myosin V and muscle myosin II. The model makes explicit the energetic cost of accelerating ADP release, showing that acceleration of ADP release during myosin V processivity requires approximately 4 kT of energy whereas the energetic cost for accelerating ADP release in a myosin II-based actin motility assay is only approximately 0.4 kT. The model also predicts that the acceleration of ADP release involves a dissipation of interhead forces. To test this prediction, we use an in vitro motility assay to show that the acceleration of ADP release from both smooth and skeletal muscle myosin II correlates with a decrease in interhead force. Our analyses provide clear energetic constraints for models of the allosteric regulation of ADP release and provide novel, testable insights into muscle and myosin V function.

Models of protein linear molecular motors for dynamic nanodevices

Protein molecular motors are natural nano-machines that convert the chemical energy from the hydrolysis of adenosine triphosphate into mechanical work. These efficient machines are central to many biological processes, including cellular motion, muscle contraction and cell division. The remarkable energetic efficiency of the protein molecular motors coupled with their nano-scale has prompted an increasing number of studies focusing on their integration in hybrid micro- and nanodevices, in particular using linear molecular motors. The translation of these tentative devices into technologically and economically feasible ones requires an engineering, design-orientated approach based on a structured formalism, preferably mathematical. This contribution reviews the present state of the art in the modelling of protein linear molecular motors, as relevant to the future design-orientated development of hybrid dynamic nanodevices.

Regular gaits and optimal velocities for motor proteins

It has been observed in numerical experiments that adding a cargo to a motor protein can regularize its gait. Here we explain these results via asymptotic analysis on a general stochastic motor protein model. This analysis permits a computation of various observables (e.g., the mean velocity) of the motor protein and shows that the presence of the cargo also makes the velocity of the motor nonmonotone in certain control parameters (e.g., ATP concentration). As an example, we consider the case of a single myosin-V protein transporting a cargo and show that, at realistic concentrations of ATP, myosin-V operates in the regime which maximizes motor velocity. Our analysis also suggests an experimental regimen which can test the efficacy of any specific motor protein model to a greater degree than was heretofore possible.

Interaction between motor domains can explain the complex dynamics of heterodimeric kinesins

Motor proteins are active enzyme molecules that play a crucial role in many biological processes. They transform chemical energy into mechanical work and move unidirectionally along rigid cytoskeleton filaments. Single-molecule experiments indicate that motor proteins, consisting of two motor domains, move in a hand-over-hand mechanism where each subunit changes between trailing and leading positions in alternating steps, and it is assumed that these subunits do not interact with each other. However, recent experiments on heterodimeric kinesins suggest that the motion of motor domains is not independent, but rather strongly coupled and coordinated, although the mechanism of these interactions is not known. We propose a simple discrete stochastic model to describe the dynamics of homodimeric and heterodimeric two-headed motor proteins. It is argued that interactions between motor domains modify original free energy landscapes for each motor subunit, while motor proteins still move via the hand-over-hand mechanism but with different transition rates specified by the new free energy profiles. Our calculations of biophysical properties agree with experimental observations. Several ways to test the theoretical model are proposed.

Myosin V passing over Arp2/3 junctions: branching ratio calculated from the elastic lever arm model

Myosin V is a two-headed processive motor protein that walks in a hand-over-hand fashion along actin filaments. When it encounters a filament branch, formed by the Arp2/3 complex, it can either stay on the straight mother filament, or switch to the daughter filament. We study both probabilities using the elastic lever arm model for myosin V. We calculate the shapes and bending energies of all relevant configurations in which the trail head is bound to the actin filament before Arp2/3 and the lead head is bound either to the mother or to the daughter filament. Based on the assumption that the probability for a head to bind to a certain actin subunit is proportional to the Boltzmann factor obtained from the elastic energy, we calculate the mother/daughter filament branching ratio. Our model predicts a value of 27% for the daughter and 73% for the mother filament. This result is in good agreement with recent experimental data.

Mechanoenzymes under superstall and large assisting loads reveal structural features

Single-molecule experiments on the motor protein kinesin have observed runs of backsteps and thus a negative, that is, reverse mean velocity, V, under superstall loads, F; but, counterintuitively, beyond stall, V(F) displays a shallow minimum and then decreases in magnitude. Conversely, under assisting loads V(F) rises to a maximum before decreasing monotonically. By contrast, while the velocity of myosin V also saturates under assisting loads, the motor moves backward increasingly rapidly under superstall loads. For both kinesin and myosin V this behavior is implied remarkably well by simple two-state kinetic models when extrapolated to large loads. To understand the origins of such results in general mechanoenzymes, biochemical kinetic descriptions are discussed on the basis of a free-energy landscape picture. It transpires that the large-load performance is determined by the geometrical placement of the intermediate mechanochemical states of the enzymatic cycles relative to the associated transition states. Explicit criteria are presented for N-state sequential kinetics, including side-reaction chains, etc., and for parallel-pathway models. Physical colocalization of biochemically distinct states generally implies large-load velocity saturation.

Regularity and synchrony in motor proteins

We investigate the origin of the regularity and synchrony which have been observed in numerical experiments of two realistic models of molecular motors, namely Fisher-Kolomeisky's model of myosin V for vesicle transport in cells and Duke's model of myosin II for sarcomere shortening in muscle contraction. We show that there is a generic organizing principle behind the emergence of regular gait for a motor pulling a large cargo and synchrony of action of many motors pulling a single cargo. These results are surprising in that the models used are inherently stochastic, and yet they display regular behaviors in the parameter range relevant to experiments. Our results also show that these behaviors are not tied to the particular models used in these experiments, but rather are generic to a wide class of motor protein models.

A kinetic model of coordinated myosin V

We present a kinetic model for the walking of myosin V on actin under conditions of zero external force. The model includes three pathways and the termination of the processivity. Experimentally measured kinetic parameters are used in the model to obtain quantitative results. Using the model and associated parameters, we compute the proportion of the pathway containing an intermediate state, as well as the walking velocities and run lengths at various concentrations of ATP and ADP. The resulting trends agree with experimental data. The model explains the surprising experimental finding that myosin walks at a faster speed but for a shorter distance as the ATP concentration increases in the absence of ADP. It also suggests that under physiological condition ([ADP] approximately 12-50 microM), myosin walks with a higher speed and for longer distances when ATP is more abundant.

Molecular Motors: A Theorist's Perspective

Individual molecular motors, or motor proteins, are enzymatic molecules that convert chemical energy, typically obtained from the hydrolysis of ATP (adenosine triphosphate), into mechanical work and motion. Processive motor proteins, such as kinesin, dynein, and certain myosins, step unidirectionally along linear tracks, specifically microtubules and actin filaments, and play a crucial role in cellular transport processes, organization, and function. In this review some theoretical aspects of motor-protein dynamics are presented in the light of current experimental methods that enable the measurement of the biochemical and biomechanical properties on a single-molecule basis. After a brief discussion of continuum ratchet concepts, we focus on discrete kinetic and stochastic models that yield predictions for the mean velocity, V(F, [ATP], ...), and other observables as a function of an imposed load force F, the ATP concentration, and other variables. The combination of appropriate theory with single-molecule observations should help uncover the mechanisms underlying motor-protein function.

Dwell time symmetry in random walks and molecular motors

The statistics of steps and dwell times in reversible molecular motors differ from those of cycle completion in enzyme kinetics. The reason is that a step is only one of several transitions in the mechanochemical cycle. As a result, theoretical results for cycle completion in enzyme kinetics do not apply to stepping data. To allow correct parameter estimation, and to guide data analysis and experiment design, a theoretical treatment is needed that takes this observation into account. In this article, we model the distribution of dwell times and number of forward and backward steps using first passage processes, based on the assumption that forward and backward steps correspond to different directions of the same transition. We extend recent results for systems with a single cycle and consider the full dwell time distributions as well as models with multiple pathways, detectable substeps, and detachments. Our main results are a symmetry relation for the dwell time distributions in reversible motors, and a relation between certain relative step frequencies and the free energy per cycle. We demonstrate our results by analyzing recent stepping data for a bacterial flagellar motor, and discuss the implications for the efficiency and reversibility of the force-generating subunits.

Kinetic models of redox-coupled proton pumping

Cytochrome c oxidase, the terminal enzyme of the respiratory chain, pumps protons across the inner mitochondrial membrane against an opposing electrochemical gradient by reducing oxygen to water. To explore the fundamental mechanisms of such redox-coupled proton pumps, we develop kinetic models at the single-molecule level consistent with basic physical principles. We demonstrate that pumping against potentials >150 mV can be achieved purely through electrostatic couplings, given an asymmetric arrangement of charge centers; however, nonlinear gates are essential for highly efficient real enzymes. The fundamental requirements for proton pumping identified here highlight a possible evolutionary origin of cytochrome c oxidase pumping. The general design principles are relevant also for other molecular machines and suggest future applications in biology-inspired fuel cells.

Back-stepping, hidden substeps, and conditional dwell times in molecular motors

Processive molecular motors take more-or-less uniformly sized steps, along spatially periodic tracks, mostly forwards but increasingly backwards under loads. Experimentally, the major steps can be resolved clearly within the noise but one knows biochemically that one or more mechanochemical substeps remain hidden in each enzymatic cycle. In order to properly interpret experimental data for back-to-forward step ratios, mean conditional step-to-step dwell times, etc., a first-passage analysis has been developed that takes account of hidden substeps in N -state sequential models. The explicit, general results differ significantly from previous treatments that identify the observed steps with complete mechanochemical cycles; e.g., the mean dwell times tau(+) and tau(-) prior to forward and back steps, respectively, are normally unequal although the dwell times tau(++) and tau(--) between successive forward and back steps are equal. Illustrative (N=2) -state examples display a wide range of behavior. The formulation extends to the case of two or more detectable transitions in a multistate cycle with hidden substeps.

The dissociation rate affects the directed motion of the molecular motor

Master equation is used in the paper under the condition of dissociation rate to study the directed motion of molecular motor. The drift velocity, the diffusion constant, the dwell time and the mean length of the motor are calculated. We found the dissociation rate affects the directed motion of the molecular motor.

Simple growth models of rigid multifilament biopolymers

The growth dynamics of rigid biopolymers, consisting of N parallel protofilaments, is investigated theoretically using simple approximate models. In our approach, the structure of a polymer's growing end and lateral interactions between protofilaments are explicitly taken into account, and it is argued that only few configurations are important for a biopolymer's growth. As a result, exact analytic expressions for growth velocity and dispersion are obtained for any number of protofilaments and arbitrary geometry of the growing end of the biopolymer. Our theoretical predictions are compared with a full description of biopolymer growth dynamics for the simplest N=2 model. It is found that the results from the approximate theory are approaching the exact ones for large lateral interactions between the protofilaments. Our theory is also applied to analyze the experimental data on the growth of microtubules.

Flexible light-chain and helical structure of F-actin explain the movement and step size of myosin-VI

Myosin-VI is a dimeric isoform of unconventional myosins. Single molecule experiments indicate that myosin-VI and myosin-V are processive molecular motors, but travel toward opposite ends of filamentous actin. Structural studies show several differences between myosin-V and VI, including a significant difference in the light-chain domain connecting the motor domains. Combining the measured kinetics of myosin-VI with the elasticity of the light chains, and the helical structure of F-actin, we compare and contrast the motility of myosin-VI with myosin-V. We show that the elastic properties of the light-chain domain control the stepping behavior of these motors. Simple models incorporating the motor elastic energy can quantitatively capture most of the observed data. Implications of our result for other processive motors are discussed.

An elastically tethered viscous load imposes a regular gait on the motion of myosin-V. Simulation of the effect of transient force relaxation on a stochastic process

Myosin-V is a processive molecular motor that moves membrane vesicles along actin tracks. In the simple model for motor and cargo motion investigated here, an elastic connection between motor and cargo transiently absorbs the abrupt mechanical transitions of the motor, and allows smooth relaxation of the cargo to a new position. We use a stochastic description to model motor stepping, with kinetics that depends on the instantaneous force exerted on the motor through the elastic connection. Tether relaxation is modelled as a continuous process, in which the rate is determined by the viscous drag of the cargo and the stiffness profile of the connection. Quantitative combined stochastic-continuous simulation of the dynamics of this system shows that bulky loads can impose a highly regular gait on the motor. If the characteristics of the elastic connection are similar to those of the myosin-II coiled-coil domain, the myosin-V motor, tether and cargo form a true escapement, in which the motor only escapes from its current position after one or more force thresholds have been crossed. Multiple thresholds limit the variation in tether length to values below that of the total step size.