Design principles governing the motility of myosin V

Modeling the motion of disease-associated KIF1A heterodimers

KIF1A is a member of the kinesin-3 motor protein family that transports synaptic vesicle precursors in axons. Mutations in the Kif1a gene cause neuronal diseases. Most patients are heterozygous and have both mutated and intact KIF1A alleles, suggesting that heterodimers composed of wild-type KIF1A and mutant KIF1A are likely involved in pathogenesis. In this study, we propose mathematical models to describe the motility of KIF1A heterodimers composed of wild-type KIF1A and mutant KIF1A. Our models precisely describe run length, run time, and velocity of KIF1A heterodimers using a few parameters obtained from two homodimers. The first model is a simple hand-over-hand model in which stepping and detachment rates from a microtubule of each head are identical to those in the respective homodimers. Although the velocities of heterodimers expected from this model were in good agreement with the experimental results, this model underestimated the run lengths and run times of some heterodimeric motors. To address this discrepancy, we propose the tethered-head affinity model, in which we hypothesize a tethered head, in addition to a microtubule-binding head, contributes to microtubule binding in a vulnerable one-head-bound state. The run lengths and run times of the KIF1A heterodimers predicted by the tethered-head affinity model matched well with experimental results, suggesting a possibility that the tethered head affects the microtubule binding of KIF1A. Our models provide insights into how each head contributes to the processive movement of KIF1A and can be used to estimate motile parameters of KIF1A heterodimers.

Kinetic Asymmetry versus Dissipation in the Evolution of Chemical Systems as Exemplified by Single Enzyme Chemotaxis

Single enzyme chemotaxis is a phenomenon by which a nonequilibrium spatial distribution of an enzyme is created and maintained by concentration gradients of the substrate and product of the catalyzed reaction. These gradients can arise either naturally through metabolism or experimentally, e.g., by flow of materials through microfluidic channels or by use of diffusion chambers with semipermeable membranes. Numerous hypotheses regarding the mechanism of this phenomenon have been proposed. Here, we discuss a mechanism based solely on diffusion and chemical reaction and show that kinetic asymmetry, a difference in the transition state energies for dissociation/association of substrate and product, and diffusion asymmetry, a difference in the diffusivities of the bound and free forms of the enzyme, are the determinates of the direction of chemotaxis and can result in either positive or negative chemotaxis, both of which have been demonstrated experimentally. Exploration of these fundamental symmetries that govern nonequilibrium behavior helps to distinguish between possible mechanisms for the evolution of a chemical system from initial to the steady state and whether the principle that determines the direction a system shifts when exposed to an external energy source is based on thermodynamics or on kinetics with the latter being supported by the results of the present paper. Our results show that, while dissipation ineluctably accompanies nonequilibrium phenomena, including chemotaxis, systems do not evolve to maximize or minimize dissipation but rather to attain greater kinetic stability and accumulate in regions where their effective diffusion coefficient is as small as possible. The chemotactic response to the chemical gradients formed by other enzymes participating in a catalytic cascade provides a mechanism for forming loose associations known as metabolons. Significantly, the direction of the effective force due to these gradients depends on the kinetic asymmetry of the enzyme and so can be nonreciprocal, where one enzyme is attracted to another enzyme, but the other enzyme is repelled by the one, in seeming contradiction to Newtons third law. This nonreciprocity is an important ingredient in the behavior of active matter.

Theory and simulations of condensin mediated loop extrusion in DNA

Condensation of hundreds of mega-base-pair-long human chromosomes in a small nuclear volume is a spectacular biological phenomenon. This process is driven by the formation of chromosome loops. The ATP consuming motor, condensin, interacts with chromatin segments to actively extrude loops. Motivated by real-time imaging of loop extrusion (LE), we created an analytically solvable model, predicting the LE velocity and step size distribution as a function of external load. The theory fits the available experimental data quantitatively, and suggests that condensin must undergo a large conformational change, induced by ATP binding, bringing distant parts of the motor to proximity. Simulations using a simple model confirm that the motor transitions between an open and a closed state in order to extrude loops by a scrunching mechanism, similar to that proposed in DNA bubble formation during bacterial transcription. Changes in the orientation of the motor domains are transmitted over ~50 nm, connecting the motor head and the hinge, thus providing an allosteric basis for LE.

Five-Site Model for Brownian Dynamics Simulations of a Molecular Walker in Three Dimensions

Motor proteins play an important role in many biological processes and have inspired the development of synthetic analogues. Molecular walkers, such as kinesin, dynein, and myosin V, fulfill a diverse set of functions including transporting cargo along tracks, pulling molecules through membranes, and deforming fibers. The complexity of molecular motors and their environment makes it difficult to model the detailed dynamics of molecular walkers over long time scales. In this work, we present a simple, three-dimensional model for a molecular walker on a bead-spring substrate. The walker is represented by five spherically symmetric particles that interact through common intermolecular potentials and can be simulated efficiently in Brownian dynamics simulations. The movement of motor protein walkers entails energy conversion through ATP hydrolysis while artificial motors typically rely on a local conversion of energy supplied through external fields. We model energy conversion through rate equations for mechanochemical states that couple positional and chemical degrees of freedom and determine the walker conformation through interaction potential parameters. We perform Brownian dynamics simulations for two scenarios: In the first, the model walker transports cargo by walking on a substrate whose ends are fixed. In the second, a tethered motor pulls a mobile substrate chain against a variable force. We measure relative displacements and determine the effects of cargo size and retarding force on the efficiency of the walker. We find that, while the efficiency of our model walker is less than for the biological system, our simulations reproduce trends observed in single-molecule experiments on kinesin. In addition, the model and simulation method presented here can be readily adapted to biological and synthetic systems with multiple walkers.

Mechanochemical enzymes and protein machines as hydrodynamic force dipoles: the active dimer model

Mechanochemically active enzymes change their shapes within every turnover cycle. Therefore, they induce circulating flows in the solvent around them and behave as oscillating hydrodynamic force dipoles. Because of non-equilibrium fluctuating flows collectively generated by the enzymes, mixing in the solution and diffusion of passive particles within it are expected to get enhanced. Here, we investigate the intensity and statistical properties of such force dipoles in the minimal active dimer model of a mechanochemical enzyme. In the framework of this model, novel estimates for hydrodynamic collective effects in solution and in lipid bilayers under rapid rotational diffusion are derived, and available experimental and computational data is examined.

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.

Processivity and Velocity for Motors Stepping on Periodic Tracks

Processive molecular motors enable cargo transportation by assembling into dimers capable of taking several consecutive steps along a cytoskeletal filament. In the well-accepted hand-over-hand stepping mechanism, the trailing motor detaches from the track and binds the filament again in the leading position. This requires fuel consumption in the form of ATP hydrolysis and coordination of the catalytic cycles between the leading and the trailing heads. Alternate stepping pathways also exist, including inchworm-like movements, backward steps, and foot stomps. Whether all the pathways are coupled to ATP hydrolysis remains to be determined. Here, to establish the principles governing the dynamics of processive movement, we present a theoretical framework that includes all of the alternative stepping mechanisms. Our theory bridges the gap between the elemental rates describing the biochemical and structural transitions in each head and the experimentally measurable quantities such as velocity, processivity, and probability of backward stepping. Our results, obtained under the assumption that the track is periodic and infinite, provide expressions that hold regardless of the topology of the network connecting the intermediate states, and are therefore capable of describing the function of any molecular motor. We apply the theory to myosin VI, a motor that takes frequent backward steps and moves forward with a combination of hand-over-hand and inchworm-like steps. Our model quantitatively reproduces various observables of myosin VI motility reported by four experimental groups. The theory is used to predict the gating mechanism, the pathway for backward stepping, and the energy consumption as a function of ATP concentration.

Myosin V executes steps of variable length via structurally constrained diffusion

The molecular motor myosin V transports cargo by stepping on actin filaments, executing a random diffusive search for actin binding sites at each step. A recent experiment suggests that the joint between the myosin lever arms may not rotate freely, as assumed in earlier studies, but instead has a preferred angle giving rise to structurally constrained diffusion. We address this controversy through comprehensive analytical and numerical modeling of myosin V diffusion and stepping. When the joint is constrained, our model reproduces the experimentally observed diffusion, allowing us to estimate bounds on the constraint energy. We also test the consistency between the constrained diffusion model and previous measurements of step size distributions and the load dependence of various observable quantities. The theory lets us address the biological significance of the constrained joint and provides testable predictions of new myosin behaviors, including the stomp distribution and the run length under off-axis force.

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 non-tight chemomechanical coupling model for force-dependence of movement dynamics of molecular motors

Based on the available experimental evidence, we present a simple and general model to describe the movement dynamics of molecular motors that can move processively on their linear tracks by using the chemical energy derived from ATP hydrolysis. An important aspect of the model is the non-tight coupling between the ATP hydrolysis and mechanical stepping, in contrast to the prevailing models presented in the literature that assume the tight chemomechanical coupling. With kinesin as an example, based on the current model, we study in detail its movement dynamics under a backward load, reproducing well the diverse available single-molecule experimental data such as the forward to backward step ratio, velocity, dwell time, randomness, run length, etc., versus the load. Moreover, predicted results are provided on the force-dependence of the mean number of ATP molecules consumed per mechanical step. Additionally, the theoretical data for the dynamics of myosin-V obtained based on the model are also in good agreement with the available experimental data.

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.

Dynamics of Allosteric Transitions in Dynein

Cytoplasmic dynein, whose motor domain belongs to the AAA+ family, walks on microtubules toward the minus end. Using the available structures in different nucleotide states, we performed simulations of a coarse-grained model to elucidate the dynamics of allosteric transitions. Binding of ATP closes the cleft between the AAA1 and AAA2 domains, triggering conformational changes in the rest of the motor domain, thus forming the pre-power stroke state. Interactions with the microtubule, modeled implicitly, enhance ADP release rate, and the formation of the post-power stroke state. The dynamics of the linker (LN), which reversibly changes from a straight to a bent state, is heterogeneous. Persistent interactions between the LN and the insert loops in the AAA2 domain prevent the formation of pre-power stroke state when ATP is bound to AAA3, thus locking dynein in a repressed non-functional state. Application of mechanical force to the LN restores motility in the repressed state.

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.

Parsing the roles of neck-linker docking and tethered head diffusion in the stepping dynamics of kinesin

Kinesin walks processively on microtubules (MTs) in an asymmetric hand-over-hand manner consuming one ATP molecule per 16-nm step. The individual contributions due to docking of the approximately 13-residue neck linker to the leading head (deemed to be the power stroke) and diffusion of the trailing head (TH) that contributes in propelling the motor by 16 nm have not been quantified. We use molecular simulations by creating a coarse-grained model of the MT-kinesin complex, which reproduces the measured stall force as well as the force required to dislodge the motor head from the MT, to show that nearly three-quarters of the step occurs by bidirectional stochastic motion of the TH. However, docking of the neck linker to the leading head constrains the extent of diffusion and minimizes the probability that kinesin takes side steps, implying that both the events are necessary in the motility of kinesin and for the maintenance of processivity. Surprisingly, we find that during a single step, the TH stochastically hops multiple times between the geometrically accessible neighboring sites on the MT before forming a stable interaction with the target binding site with correct orientation between the motor head and the [Formula: see text] tubulin dimer.

Robust mechanobiological behavior emerges in heterogeneous myosin systems

Biological complexity presents challenges for understanding natural phenomenon and engineering new technologies, particularly in systems with molecular heterogeneity. Such complexity is present in myosin motor protein systems, and computational modeling is essential for determining how collective myosin interactions produce emergent system behavior. We develop a computational approach for altering myosin isoform parameters and their collective organization, and support predictions with in vitro experiments of motility assays with α-actinins as molecular force sensors. The computational approach models variations in single myosin molecular structure, system organization, and force stimuli to predict system behavior for filament velocity, energy consumption, and robustness. Robustness is the range of forces where a filament is expected to have continuous velocity and depends on used myosin system energy. Myosin systems are shown to have highly nonlinear behavior across force conditions that may be exploited at a systems level by combining slow and fast myosin isoforms heterogeneously. Results suggest some heterogeneous systems have lower energy use near stall conditions and greater energy consumption when unloaded, therefore promoting robustness. These heterogeneous system capabilities are unique in comparison with homogenous systems and potentially advantageous for high performance bionanotechnologies. Findings open doors at the intersections of mechanics and biology, particularly for understanding and treating myosin-related diseases and developing approaches for motor molecule-based technologies.

Reexamining the origin of the directionality of myosin V

The nature of the conversion of chemical energy to directional motion in myosin V is examined by careful simulations that include two complementary methods: direct Langevin Dynamics (LD) simulations with a scaled-down potential that provided a detailed time-resolved mechanism, and kinetic equations solution for the ensemble long-time propagation (based on information collected for segments of the landscape using LD simulations and experimental information). It is found that the directionality is due to the rate-limiting ADP release step rather than the potential energy of the lever arm angle. We show that the energy of the power stroke and the barriers involved in it are of minor consequence to the selectivity of forward over backward steps and instead suggest that the selective release of ADP from a postrigor myosin motor head promotes highly selective and processive myosin V. Our model is supported by different computational methods-LD simulations, Monte Carlo simulations, and kinetic equations solution-as well as by structure-based binding energy calculations.

Simulating the dynamics of the mechanochemical cycle of myosin-V

The detailed dynamics of the cycle of myosin-V are explored by simulation approaches, examining the nature of the energy-driven motion. Our study started with Langevin dynamics (LD) simulations on a very coarse landscape with a single rate-limiting barrier and reproduced the stall force and the hand-over-hand dynamics. We then considered a more realistic landscape and used time-dependent Monte Carlo (MC) simulations that allowed trajectories long enough to reproduce the force/velocity characteristic sigmoidal correlation, while also reproducing the hand-over-hand motion. Overall, our study indicated that the notion of a downhill lever-up to lever-down process (popularly known as the powerstroke mechanism) is the result of the energetics of the complete myosin-V cycle and is not the source of directional motion or force generation on its own. The present work further emphasizes the need to use well-defined energy landscapes in studying molecular motors in general and myosin in particular.

Quantifying the Heat Dissipation from a Molecular Motor's Transport Properties in Nonequilibrium Steady States

Theoretical analysis, which maps single-molecule time trajectories of a molecular motor onto unicyclic Markov processes, allows us to evaluate the heat dissipated from the motor and to elucidate its dependence on the mean velocity and diffusivity. Unlike passive Brownian particles in equilibrium, the velocity and diffusion constant of molecular motors are closely inter-related. In particular, our study makes it clear that the increase of diffusivity with the heat production is a natural outcome of active particles, which is reminiscent of the recent experimental premise that the diffusion of an exothermic enzyme is enhanced by the heat released from its own catalytic turnover. Compared with freely diffusing exothermic enzymes, kinesin-1, whose dynamics is confined on one-dimensional tracks, is highly efficient in transforming conformational fluctuations into a locally directed motion, thus displaying a significantly higher enhancement in diffusivity with its turnover rate. Putting molecular motors and freely diffusing enzymes on an equal footing, our study offers a thermodynamic basis to understand the heat-enhanced self-diffusion of exothermic enzymes.

Strain Mediated Adaptation Is Key for Myosin Mechanochemistry: Discovering General Rules for Motor Activity

A structure-based model of myosin motor is built in the same spirit of our early work for kinesin-1 and Ncd towards physical understanding of its mechanochemical cycle. We find a structural adaptation of the motor head domain in post-powerstroke state that signals faster ADP release from it compared to the same from the motor head in the pre-powerstroke state. For dimeric myosin, an additional forward strain on the trailing head, originating from the postponed powerstroke state of the leading head in the waiting state of myosin, further increases the rate of ADP release. This coordination between the two heads is the essence of the processivity of the cycle. Our model provides a structural description of the powerstroke step of the cycle as an allosteric transition of the converter domain in response to the P_i release. Additionally, the variation in structural elements peripheral to catalytic motor domain is the deciding factor behind diverse directionalities of myosin motors (myosin V & VI). Finally, we observe that there are general rules for functional molecular motors across the different families. Allosteric structural adaptation of the catalytic motor head in different nucleotide states is crucial for mechanochemistry. Strain-mediated coordination between motor heads is essential for processivity and the variation of peripheral structural elements is essential for their diverse functionalities.

Molecular motors are perhaps the most important proteins present in the cell. The importance specifically lies with the fact that these proteins use the chemical energy source (such as ATP) of the cell to generate mechanical work and perform a wide range of functionalities. In this article, we generalize the idea of using structure-based models to explore the mechanochemistry of myosin molecular motors in structural terms. We find that a structural adaptation of the motor head domain in post-powerstroke state signals faster ADP release from the trailing head to maintain its processivity while directionality arises from a careful design of peripheral structural elements. These results along with our earlier results on other motors provide a general rule for motor activity.

Actomyosin contraction, aggregation and traveling waves in a treadmilling actin array

We use perturbation theory to derive a continuum model for the dynamic actomyosin bundle/ring in the regime of very strong crosslinking. Actin treadmilling is essential for contraction. Linear stability analysis and numerical solutions of the model equations reveal that when the actin treadmilling is very slow, actin and myosin aggregate into equidistantly spaced peaks. When treadmilling is significant, actin filament of one polarity are distributed evenly, while filaments of the opposite polarity develop a shock wave moving with the treadmilling velocity. Myosin aggregates into a sharp peak surfing the crest of the actin wave. Any actomyosin aggregation diminishes contractile stress. The easiest way to maintain higher contraction is to upregulate the actomyosin turnover which destabilizes nontrivial patterns and stabilizes the homogeneous actomyosin distributions. We discuss the model's implications for the experiment.

Protein folding guides disulfide bond formation

The Anfinsen principle that the protein sequence uniquely determines its structure is based on experiments on oxidative refolding of a protein with disulfide bonds. The problem of how protein folding drives disulfide bond formation is poorly understood. Here, we have solved this long-standing problem by creating a general method for implementing the chemistry of disulfide bond formation and rupture in coarse-grained molecular simulations. As a case study, we investigate the oxidative folding of bovine pancreatic trypsin inhibitor (BPTI). After confirming the experimental findings that the multiple routes to the folded state contain a network of states dominated by native disulfides, we show that the entropically unfavorable native single disulfide [14-38] between Cys14 and Cys38 forms only after polypeptide chain collapse and complete structuring of the central core of the protein containing an antiparallel β-sheet. Subsequent assembly, resulting in native two-disulfide bonds and the folded state, involves substantial unfolding of the protein and transient population of nonnative structures. The rate of [14-38] formation increases as the β-sheet stability increases. The flux to the native state, through a network of kinetically connected native-like intermediates, changes dramatically by altering the redox conditions. Disulfide bond formation between Cys residues not present in the native state are relevant only on the time scale of collapse of BPTI. The finding that formation of specific collapsed native-like structures guides efficient folding is applicable to a broad class of single-domain proteins, including enzyme-catalyzed disulfide proteins.

Structural dynamics of myosin 5 during processive motion revealed by interferometric scattering microscopy

Myosin 5a is a dual-headed molecular motor that transports cargo along actin filaments. By following the motion of individual heads with interferometric scattering microscopy at nm spatial and ms temporal precision we found that the detached head occupies a loosely fixed position to one side of actin from which it rebinds in a controlled manner while executing a step. Improving the spatial precision to the sub-nm regime provided evidence for an ångstrom-level structural transition in the motor domain associated with the power stroke. Simultaneous tracking of both heads revealed that consecutive steps follow identical paths to the same side of actin in a compass-like spinning motion demonstrating a symmetrical walking pattern. These results visualize many of the critical unknown aspects of the stepping mechanism of myosin 5 including head-head coordination, the origin of lever-arm motion and the spatiotemporal dynamics of the translocating head during individual steps.

Structural dynamics of myosin 5 during processive motion revealed by interferometric scattering microscopy

Myosin 5a is a dual-headed molecular motor that transports cargo along actin filaments. By following the motion of individual heads with interferometric scattering microscopy at nm spatial and ms temporal precision we found that the detached head occupies a loosely fixed position to one side of actin from which it rebinds in a controlled manner while executing a step. Improving the spatial precision to the sub-nm regime provided evidence for an ångstrom-level structural transition in the motor domain associated with the power stroke. Simultaneous tracking of both heads revealed that consecutive steps follow identical paths to the same side of actin in a compass-like spinning motion demonstrating a symmetrical walking pattern. These results visualize many of the critical unknown aspects of the stepping mechanism of myosin 5 including head-head coordination, the origin of lever-arm motion and the spatiotemporal dynamics of the translocating head during individual steps.

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.

Coupling of lever arm swing and biased Brownian motion in actomyosin

An important unresolved problem associated with actomyosin motors is the role of Brownian motion in the process of force generation. On the basis of structural observations of myosins and actins, the widely held lever-arm hypothesis has been proposed, in which proteins are assumed to show sequential structural changes among observed and hypothesized structures to exert mechanical force. An alternative hypothesis, the Brownian motion hypothesis, has been supported by single-molecule experiments and emphasizes more on the roles of fluctuating protein movement. In this study, we address the long-standing controversy between the lever-arm hypothesis and the Brownian motion hypothesis through in silico observations of an actomyosin system. We study a system composed of myosin II and actin filament by calculating free-energy landscapes of actin-myosin interactions using the molecular dynamics method and by simulating transitions among dynamically changing free-energy landscapes using the Monte Carlo method. The results obtained by this combined multi-scale calculation show that myosin with inorganic phosphate (P_i) and ADP weakly binds to actin and that after releasing P_i and ADP, myosin moves along the actin filament toward the strong-binding site by exhibiting the biased Brownian motion, a behavior consistent with the observed single-molecular behavior of myosin. Conformational flexibility of loops at the actin-interface of myosin and the N-terminus of actin subunit is necessary for the distinct bias in the Brownian motion. Both the 5.5–11 nm displacement due to the biased Brownian motion and the 3–5 nm displacement due to lever-arm swing contribute to the net displacement of myosin. The calculated results further suggest that the recovery stroke of the lever arm plays an important role in enhancing the displacement of myosin through multiple cycles of ATP hydrolysis, suggesting a unified movement mechanism for various members of the myosin family.

Myosin II is a molecular motor that is fueled by ATP hydrolysis and generates mechanical force by interacting with actin filament. Comparison among various myosin structures obtained by X-ray and electron microscope analyses has led to the hypothesis that structural change of myosin in ATP hydrolysis cycle is the driving mechanism of force generation. However, single-molecule experiments have suggested an alternative mechanism in which myosin moves stochastically in a biased direction along actin filament. Computer simulation serves as a platform for assessing these hypotheses by revealing the prominent features of the dynamically changing landscape of actin-myosin interaction. The calculated results show that myosin binds to actin at different locations of actin filament in the weak- and strong-binding states and that the free energy has a global gradient from the weak-binding site to the strong-binding site. Myosin relaxing into the strong-binding state therefore necessarily shows the biased Brownian motion toward the strong-binding site. Lever-arm swing is induced during this relaxation process; therefore, lever-arm swing and the biased Brownian motion are coupled to contribute to the net displacement of myosin. This coupling should affect the dynamical behaviors of muscle and cardiac systems.

Electrostatic origin of the unidirectionality of walking myosin V motors

Understanding the basis for the action of myosin motors and related molecular machines requires a quantitative energy-based description of the overall functional cycle. Previous theoretical attempts to do so have provided interesting insights on parts of the cycle but could not generate a structure-based free energy landscape for the complete cycle of myosin. In particular, a nonphenomenological structure/energy-based understanding of the unidirectional motion is still missing. Here we use a coarse-grained model of myosin V and generate a structure-based free energy surface of the largest conformational change, namely the transition from the post- to prepowerstroke movement. We also couple the observed energetics of ligand binding/hydrolysis and product release to that of the conformational surface and reproduce the energetics of the complete mechanochemical cycle. It is found that the release in electrostatic free energy upon changing the conformation of the lever arm and the convertor domain from its post- to prepowerstroke state provides the necessary energy to bias the system towards the unidirectional movement of myosin V on the actin filament. The free energy change of 11 kcal is also in the range of ∼2-3 pN, which is consistent with the experimentally observed stalling force required to stop the motor completely on its track. The conformational-chemical coupling generating a successful powerstroke cycle is believed to be conserved among most members of the myosin family, thus highlighting the importance of the previously unknown role of electrostatics free energy in guiding the functional cycle in other actin-based myosin motors.