Myosin V executes steps of variable length via structurally constrained diffusion
Abstract
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 offaxis force.
Introduction
Molecular motors are cellular machines that function by converting chemical energy into mechanical work (Schliwa and Woehlke, 2003). Motors play key roles in many intracellular biological processes, including signaling and the transport of cargo (Sun and Goldman, 2011). Members of the myosin superfamily, one class of molecular motors, perform these functions by binding to actin filaments and generating energy through ATP hydrolysis. Myosin V, a dimeric transport motor, is composed of two stiff polymer chains joined at a pivot with an actinbinding head at the end of each chain (ReckPeterson et al., 2000). The motor walks forward along the actin, stepping handoverhand, by alternating head detachment, with the free head performing a diffusive search for actin binding sites during each step (Shiroguchi and Kinosita, 2007). Such unidirectional motility requires coordination between the two heads with preferential detachment of the rear head, the so called ‘gating’ mechanism, which is regulated by the strain within the lever arms while the heads are bound to actin (Veigel et al., 2002; Veigel et al., 2005; Purcell et al., 2005; Sakamoto et al., 2008). Myosin V propels itself toward the plus (barbed) end of the actin using two changes in the lever arm orientation. The power stroke, executed by an actinbound head, swings the lever arm forward, while the recovery stroke, executed during diffusion, returns the lever to its original orientation, which favors binding to forward actin sites (Shiroguchi et al., 2011). Most frequently, myosin V takes ≈ 74 nm steps, roughly equal in length to the actin helical pitch, but shorter and longer steps also occur (Yildiz et al., 2003). The near correspondence of the step size with the pitch, as well as the narrow step distribution, allows myosin V to approximately maintain its azimuthal orientation with respect to the actin over multiple steps (Sun and Goldman, 2011). Mutant myosins with altered lever arms show a linear relation between mean step size and arm length (Sakamoto et al., 2005; Oke et al., 2010). Under moderate backward force myosin V remains highly processive, up to a stall force ≈ 1.9  3 pN, (Mehta et al., 1999; Veigel et al., 2002; Kad et al., 2008; Uemura et al., 2004; Cappello et al., 2007; Gebhardt et al., 2006), at which the mean velocity of the motor goes to zero.
While myosin V is one of the most extensively studied motors, new functional features continue to be discovered as the spatiotemporal resolution of experimental imaging improves. Most recently Andrecka et al. (2015) achieved simultaneous millisecond temporal and nanometer spatial resolution with interferometric scattering microscopy in which the position of the gold nanoparticle attached to one of the motor heads was used to track the diffusion of the free head during its step. This measurement and recent electron micrographs of freely floating myosin taken by Takagi et al. (2014) indicate that the joint between the myosin V lever arms does not rotate freely, but instead has structural constraints giving rise to a preferred interarm joint angle. The presence of a joint constraint seems to contradict previous diffusion measurements by Dunn and Spudich (2007) and other experiments (Shiroguchi and Kinosita, 2007; Fujita et al., 2012; Beausang et al., 2013) which indicated a freely rotating joint. Further, a number of theoretical and numerical studies based on free diffusion models have been remarkably successful in quantitatively describing various aspects of myosin V motility (Hinczewski et al., 2013; Craig and Linke, 2009; Mukherjee et al., 2017).
To address these apparent conflicts, we extend a minimal model of myosin V previously introduced by Hinczewski et al. (2013), incorporating a constraint on the relative orientation of the two lever arms. The model combines a coarsegrained polymeric description of the diffusive search (Thirumalai and Ha, 1998) with the reaction network of discrete states taken by the motor heads during the mechanochemical stepping cycle (Vilfan, 2005a; Bierbaum and Lipowsky, 2011; Sumi, 2017). The large persistence length of the myosin V lever arms allows us to derive an approximate but accurate semianalytical expression for the equilibrium distribution of positions occupied by the free head during the diffusive search. The kinetic network accounts for not only forward steps, but also foot stomps (the head reattaching near the site of detachment) and backward steps which have been observed experimentally (Kodera et al., 2010) and become more prominent as the resistive force increases. In contrast to the simplified model of Hinczewski et al. (2013), here we include the full set of available binding sites on the doublehelical actin filaments, enabling a description of the distributions of steps and stomps taken by myosin V. The backward force applied by the load induces conformational changes in the lever arms, altering the diffusive search for binding sites and the associated binding probabilities. The effects of the magnitude of the resistive force and direction are easily incorporated into the theory. We supplement the analytical theory with Brownian dynamics (BD) simulations of myosin stepping dynamics, the results of which largely concur with analytical predictions. In addition to addressing the constrained diffusion hypothesis, our model provides the most comprehensive accounting to date of the full range of substall experimental data, including stepsize distributions and the load dependence of several physical observables.
The polymer model gives direct insights into the connection between structural features of myosin V, including the interarm joint constraint, and both the diffusive search and kinetics of the motor. With a joint constraint, our model predicts diffusion profiles similar to that observed by Andrecka et al. (2015). By computing the changes in diffusion as the constraint strength is varied, we estimate upper and lower bounds on the constraint energy. Fitting the model to experimental measurements of myosin V step distributions (Yildiz et al., 2003; Sakamoto et al., 2005; Oke et al., 2010) and the force dependence of the backwardtoforward step ratio (Kad et al., 2008) and mean run length/velocity (Mehta et al., 1999; Uemura et al., 2004; Clemen et al., 2005; Gebhardt et al., 2006; Kad et al., 2008), we confirm the consistency of the constrained diffusion picture with previous experimental data on myosin V. Interestingly, while the joint constraint considerably alters the diffusive search space, it has relatively small influence on the myosin V step size and force response. A free diffusion model, for instance, produces similar kinetic behaviors. Our model allows us to address questions related to the biological significance of the joint constraint. We find, for example, that the constraint does not necessarily speed up the binding time, despite narrowing the space of the diffusive search. However it does narrow the width of the forward step distributions on actin, allowing the head to more consistently target the actin binding sites at halfhelical length intervals. Finally, the model provides testable predictions of new quantities yet to be probed by experiments, including the stomp distribution near the stall regime and the robust run length under offaxis forces.
Results
Theoretical model for myosin V dynamics
In the following sections we describe the main features of our polymer structural model for myosin V, the actin filament geometry, and how the model allows us to predict diffusion and stepping behavior, including the probabilities of various kinetic pathways. The full mathematical details of the analytical theory can be found in Appendix 1. The details of the BD simulations, which we used to validate the analytical theory results, are described in Appendix 2.
Polymer model and actin filament geometry
Following Hinczewski et al. (2013), we model the actinbinding heads and lever arm domains of myosin V as semiflexible polymer chains with length $L$ and persistence length ${l}_{p}$. We will denote each polymer chain as a ‘leg’ of the motor (consisting of the combined head and lever arm domains), and refer to trailing or leading legs depending on whether a given leg is further away from or closer to the barbed (plus) end of actin, respectively. Several experiments measuring myosin step distributions and other properties have been carried out comparing wildtype myosin V to mutants where the lever arm length is altered through addition or deletion of IQ motifs (Yildiz et al., 2003; Sakamoto et al., 2005; Oke et al., 2010). We assume the actinbinding head and IQ motifs are each approximately 5 nm in length, so that wildtype myosin V (with 6IQ motifs) has $L=35\mathrm{nm}$ nm and the 4IQ and 8IQ mutants have $L=25\mathrm{nm}$ nm and $L=45\mathrm{nm}$ nm respectively. The persistence length has estimated values ranging from ${l}_{p}\approx 100\mathrm{nm}$ nm (Howard and Spudich, 1996) to ${l}_{p}\approx 375\mathrm{nm}$ nm (Vilfan, 2005a). While fitting our model to experimental data we allow the persistence length to vary within this range, though the model predictions are qualitatively similar for any ${l}_{p}\gg L$ in the stiff leg regime.
The two polymer legs are connected at a joint forming the myosin V dimer (see Figure 1 for an illustration of the geometry). Recent experimental evidence suggests this this joint does not not rotate freely, but instead has a preferred joint angle ${\theta}_{p}$ giving rise to constrained diffusion (Takagi et al., 2014; Andrecka et al., 2015). In our model, this preferred angle is enforced by a potential ${\mathscr{H}}_{J}={\mu}_{c}{k}_{B}T[1\mathrm{cos}({\theta}_{J}{\theta}_{p})]$, which is minimum when the interleg angle at the joint ${\theta}_{J}$ is equal to the preferred angle. The parameter ${\mu}_{c}$ is the constraint strength, designating the energy cost to deviations from ${\theta}_{p}$, while ${k}_{B}$ is Boltzmann’s constant and T is temperature. Note that for small angle differences ${\theta}_{J}{\theta}_{p}$, the constraint potential ${\mathscr{H}}_{J}\approx \frac{1}{2}{\mu}_{c}{k}_{B}T{({\theta}_{J}{\theta}_{p})}^{2}$ is approximately harmonic. In the limit ${\mu}_{c}\to 0$ the joint becomes freely rotating with no preferred angle, the case considered previously Hinczewski et al. (2013). This cosine potential is used throughout the paper, but below we briefly discuss how the form of the interleg potential affects the diffusion of the free head. The tail domain, which attaches to cargo, transmits a load force $\mathbf{\mathbf{F}}$ to the joint. The direction of the force is parameterized by ${\theta}_{\text{F}}$, measured clockwise from the $\widehat{\mathbf{\mathbf{z}}}$ axis, and ${\varphi}_{\text{F}}$, measured counterclockwise from the $\widehat{\mathbf{\mathbf{x}}}$ axis. Our main focus is on the behavior of myosin V under zero force and backward force (${\theta}_{\text{F}}={\varphi}_{\text{F}}=0$), but we also consider offaxis forces at the end of the Results section.
The actinbinding heads, located at the ends of the polymer legs in our model, can bind to various sites along the doublehelical filamentous actin structure. Actin is composed of two filaments each containing 13 actin subunits per helical rotation, with one binding site per subunit. The filaments run parallel to the $\widehat{\mathbf{\mathbf{z}}}$ axis, leading to a geometry in which the binding sites ${\mathbf{\mathbf{r}}}_{n}$, $n=0,\pm 1,\pm 2,\mathrm{\dots}$ have positions
where $R=5.5\mathrm{nm}$ nm is the radius of the helix, $\mathrm{\Delta}z=72/13\approx 5.5\mathrm{nm}$ nm is the size of each actin subunit, and ${\varphi}_{n}=12\pi n/13$ is the angle between adjacent subunits (Lan and Sun, 2006). Equation 1 and all the other key analytical quantities in our theory are summarized in Table 1 for ease of reference. Even and odd $n$ respectively correspond to subunits on the 1^{st} and 2^{nd} filaments of the double helix. Wildtype myosin V steps most frequently to the halfhelical sites $n=\pm 13$ located at $z=\pm 36\mathrm{nm}$ nm, while mutants favor other sites depending on their lever arm length (Sakamoto et al., 2005; Oke et al., 2010).
When myosin V is bound to an actin filament, the lever arm can be in two orientations. After the leg has bound, but before the power stroke has been executed, the lever arm points toward the pointed (minus) end of the actin. After the power stroke, the leg rotates toward the barbed (plus) end of the actin and is held at an angle ${\theta}_{c}$ above the actin. Similar to the interleg constraint described above, in our model the preferred forward tilting angle is enacted through a harmonic potential, ${\mathscr{H}}_{c}=\frac{1}{2}{\nu}_{c}{k}_{B}T{({\widehat{\mathbf{\mathbf{u}}}}_{0}{\widehat{\mathbf{\mathbf{u}}}}_{c})}^{2}$. The constraint has strength ${\nu}_{c}$, the unit vector ${\widehat{\mathbf{\mathbf{u}}}}_{0}$ is the tangent to the bound leg at the binding point, and the unit vector ${\widehat{\mathbf{\mathbf{u}}}}_{c}$ defines the preferred direction, which lies in the $\widehat{\mathbf{\mathbf{x}}}\widehat{\mathbf{\mathbf{z}}}$ plane at an angle ${\theta}_{c}$ from the $\widehat{\mathbf{\mathbf{z}}}$ axis. As we will see below, steps and stomps occur only when the bound leg is in the postpower stroke orientation. For the purposes of modeling these aspects of the dynamics we do not need to consider a separate constraint potential for diffusion while the bound leg has the prepower stroke orientation.
Volume exclusion effects introduce another constraint on the orientation of the myosin legs during diffusion. For instance, the myosin heads are unlikely to be found in close proximity due to steric repulsion. Such exclusion interactions apply not only to the myosin heads, but also to the legs, which have to be brought close together in order to accommodate small head separations. Though these interactions can’t be explicitly included in the coarsegrained polymer model, we capture the effective repulsion between the myosin heads using the potential ${\mathscr{H}}_{V}={k}_{B}T{({d}_{\text{V}}/r)}^{6}$, where $r$ is the distance between the bound and free myosin heads and ${d}_{\text{V}}$ is the effective length scale of the repulsion. The magnitude of ${d}_{\text{V}}$ depends on the details of the interacting legs between the heads, in particular their length $L$L. Using BD simulations described below (which explicitly include all volume exclusion interactions between the myosin legs), we estimate ${d}_{\text{V}}\approx$ 20, 27.5, and 35 nm for 4IQ, 6IQ, and 8IQ myosin respectively.
Brownian Dynamics simulations
We also study the stepping dynamics of myosin V using BD simulations. In the simulations, myosin V is treated as two connected and interacting polymer chains. The Hamiltonian that prescribes the interactions and the resulting structural features like the persistence length, is outlined in Appendix 2. A BD trajectory begins with the trailing lever arm unbinding from actin (mimicking the effect of ATP binding). A power stroke moves the unbound head quickly forward, after which it executes a diffusive search until it finds an energetically favorable actin binding site. As in the analytical model, the actin binding sites are given by Equation 1. We trace the diffusive search and record the binding locations. The reported data is averaged over 200 randomized trajectories. The simulations allow us to include the finite sizes of the lever arms, the effect of the volume excluded by actin, and a glass coverslip that is used in certain experimental setups. Additionally we can test different functional forms for the binding criterion and power stroke execution.
Kinetic pathways
Myosin V can exhibit a wide range of behaviors, shown schematically in Figure 2, including steps and stomps involving each polymer leg. Before following a particular kinetic pathway, Myosin V starts in a waiting state where both legs are strongly bound to the actin and each head has an associated ADP molecule. Each leg is in the postpower stroke orientation, pointed toward the plus end of the actin, and the leading leg is bent backward under tension. The motor then goes through one of many kinetic pathways, which can be broken into the five categories described below.
Forward steps
The myosin sits in the waiting state until the ADP molecule unbinds from the trailing head, allowing an ATP molecule to take its place. With the bound ATP, the trailing leg is only weakly associated to the actin and therefore quickly detaches. We assume (except for modeling one run velocity experiment described below) that the ATP concentration is near saturation levels, so that ATP binding occurs rapidly. Then the detachment process occurs roughly at the rate of ADP dissociation ${t}_{\text{d1}}^{1}=12{\mathrm{s}}^{1}$ s^{–1}, which has been measured experimentally (De La Cruz et al., 1999). Following detachment of the trailing leg, the leading leg, previously under tension, undergoes rapid diffusive relaxation: since the leading leg had been bent backwards, away from its preferred postpower stroke orientation, the leg relaxes toward the preferred orientation, thus leading to the entire system swinging forward toward the plus end of actin (Hinczewski et al., 2013; Dunn and Spudich, 2007). The timescale ${t}_{r}$ of this relaxation depends on the load, but is generally $\lesssim$ 5 μs, based on both theoretical and numerical considerations Hinczewski et al. (2013). Over much longer timescales, the free leg diffusively searches until it reaches a binding site on the actin. If we let the bound leg be attached at the origin $\mathbf{\mathbf{r}}=0$ then the available binding sites are located at ${\mathbf{\mathbf{r}}}_{n}$ given by Equation 1, and forward steps occur when the free leg reaches a site with $n>0$, assuming the following additional binding criteria are also fulfilled.
Before the detached head can bind to actin, the ATP must hydrolyze, ATP → ADP + ${\text{P}}_{\text{i}}$, which has two primary functions. First, during hydrolysis the free head undergoes the recovery stroke, rotating the head into the prepower stroke orientation, so that binding to forward sites along the actin is conformationally favorable. Second, the hydrolysis produces an ADP molecule, which is required for the head to strongly associate with actin and successfully bind. The ATP hydrolysis occurs at rate $t}_{\text{h}}^{1$ = 750 s^{—1} (De La Cruz et al., 1999) and we assume the reverse reaction rate is negligible.
Once ATP hydrolysis is completed, the free head must diffuse close enough to the actin to appreciably interact and bind. We assume that binding occurs when the head is within a distance $a$ from an actin binding site. An upper bound for the capture radius $a$ is the Debye screening length, which under physiological conditions (KCL concentrations of 25 – 400 mM) is ${\lambda}_{\text{D}}\approx 1.90.5\mathrm{nm}$ nm (Craig and Linke, 2009). On these length scales, the detached head can interact with the binding site and we expect near certain binding at slightly shorter distances. Below we find that for $a=0.4\mathrm{nm}$ nm the model produces quantitative agreement with experimental data. Finally, in addition to requiring close proximity of the free head to an actin binding site, the detached myosin head must also be in the correct orientation for binding to occur (the full details of this binding criterion are described in the section ‘First Passage Times and the Binding Acceptance Region’ below.)
Trailing leg stomps
This kinetic pathway starts identically to the forward steps described above, namely the trailing leg dissociation is followed by ATP binding and hydrolysis. In this case, however, the diffusive search finds a site ${\mathbf{\mathbf{r}}}_{n}$ with $n<0$. Instead of executing a handoverhand step, the myosin stomps and the heads remain in the same relative order on the actin, with a comparably small change in the center of mass position of the myosin. As before, we assume the binding occurs once the head is within a distance $a$ from a binding site with the free leg in a sufficiently acceptable orientation.
The recovery stroke is executed during ATP hydrolysis and before binding, so that the prepower stroke orientation of the head favors binding to forward sites. After a forward step, the lever arm points toward the minus end of the actin in a relaxed prepower stroke state. To bind to sites behind the bound leg, however, the free lever arm must bend and point toward the plus end of the actin despite its prepower stroke orientation with respect to the head. This unnatural configuration puts the free leg under additional strain, so that there is an energy barrier to binding at backward actin sites. We model this effective barrier by reducing the probability of binding when the free head reaches a radius $a$ of actin sites with $n<0$. Instead of binding with certainty, the myosin head binds with probability $b<1$, which we call the binding penalty. The reduced binding probability, a consequence of the recovery stroke, increases the relative probability that the myosin will step forward, contributing to the biased motion along the actin filament.
Backward steps
The other possible kinetic pathways start with the leading head detaching from the actin, in contrast to the trailing head detachment which initiates forward steps and trailing leg stomps. After the power stroke, the leading lever arm is under considerable backward tension. This tension dramatically suppresses the release of ADP from the leading head by a factor of 50–70 (Kodera et al., 2010; Rosenfeld and Sweeney, 2004). Since these events are very rare, an alternative detachment pathway is dominant: the leading head dissociates from the actin but retains the associated ADP molecule (Kodera et al., 2010; Purcell et al., 2005). This detachment pathway occurs at a slower rate than trailing head detachment, ${t}_{\text{d2}}^{1}={(g{t}_{\text{d1}})}^{1}$, where $g>1$ is the gating ratio. Under backward force of ∼2 pN, singleheaded myosin V has been observed to detach from actin at a rate of 1.5 s^{—1 }(Purcell et al., 2005), while previous theoretical calculations using the above described polymer model estimate the backward force on the leading leg to be 2.7 pN (Hinczewski et al., 2013) when myosin V sits in the waiting state with both legs bound. We therefore use ${t}_{\text{d2}}^{1}=1.5{\mathrm{s}}^{1}$ s^{—1}, which corresponds to the gating ratio $g=8$.
Once the leading head unbinds from the actin, it executes a diffusive search until it reaches an actin binding site. Note that the bound leg is in the postpower stroke state, so the diffusion is statistically identical to that which occurs during forward steps and trailing leg stomps. In this case, however, the free head has an associated ADP molecule immediately after detachment, so the free leg can rebind without undergoing ATP hydrolysis. Furthermore, because ATP hydrolysis is not executed, the free head remains in the postpower stroke state orientation which favors binding to backward sites. If we let the bound leg be at the origin $\mathbf{\mathbf{r}}=0$, a backward step occurs when the diffusive search finds an actin site ${\mathbf{\mathbf{r}}}_{n}$ with $n<0$. Hence, the free head orientation favors binding to these sites and binding occurs with probability 1 when the head diffuses within a distance $a$ of the target site and the free leg is within the acceptance region.
Leading leg stomps
Similar to the trailing leg stomps described above, the myosin V can also execute a leading leg stomp (Kodera et al., 2010; Andrecka et al., 2015). This kinetic pathway begins with leading leg detachment, identical to the backward step. To perform a stomp, the free head diffuses to an actin binding site ${\mathbf{\mathbf{r}}}_{n}$ with $n>0$. Since the free head is in the prepower stroke state, binding to these forward sites requires an unnatural deformation of the free leg, which as in the case of trailing leg stomps, introduces an effective energy barrier to these events. Therefore, when the free head diffuses within a radius $a$ of a target site with the free leg in the acceptance region, binding occurs with probability $b<1$, given by the binding penalty.
Termination
The final kinetic pathway occurs when the bound myosin head detaches before the free head finds an open actin binding site, terminating the processive run of myosin V. We assume that the bound leg has a constant detachment rate equal to ${t}_{\text{d1}}^{1}$, independent of the backward load force exerted on the motor.
These five pathways complete the kinetic description of myosin V procession. The first four pathways, forward steps, leading and trailing leg stomps, and backward steps, each return the myosin to its waiting state with both legs bound, while termination is a complete dissociation from the actin. Therefore the myosin will continue to execute steps and stomps until termination which ends the run. For notational convenience below, we will denote the binding penalty of a given actin site ${\mathbf{\mathbf{r}}}_{n}$ following trailing head detachment as ${b}_{n}$. Then ${b}_{n}=1$ for $n>0$, there is no penalty for forward steps, and ${b}_{n}=b<1$ if $n<0$. The binding penalty following leading head detachment is ${b}_{n}$. If we sum the probabilities over all binding sites possible in a given pathway, we get the overall probability of each type of step and stomp. We define ${\mathcal{P}}_{\text{f}}$, ${\mathcal{P}}_{\text{Ls}}$, ${\mathcal{P}}_{\text{Ts}}$, ${\mathcal{P}}_{\text{b}}$, and ${\mathcal{P}}_{\text{t}}$ to be the overall probability of myosin executing a forward step, leading/trailing leg stomp, backward step, and termination respectively, starting from the waiting state. The above pathways are dominant when the backward load forces are below or slightly above the stall force, for which experimental estimates range from $\approx 1.93\mathrm{pN}$ pN (Mehta et al., 1999; Veigel et al., 2002; Kad et al., 2008; Uemura et al., 2004; Cappello et al., 2007; Gebhardt et al., 2006). For large superstall forces ($\gtrsim 4$ pN), the myosin can experience power stroke reversal, swinging the leading lever arm back toward the minus end of the actin (Sellers and Veigel, 2010). At these high forces the diffusive search can therefore occur with the bound leg in the prepower stroke orientation, so additional pathways must be added to the above model to accurately describe the superstall behavior of myosin V. We restrict our focus to the substall and stall regimes.
First passage times and the binding acceptance region
The key quantities necessary for our theoretical description of myosin V processivity are the mean first passage times, ${t}_{\text{fp}}^{n}$, of the free head to each of the actin binding sites ${\mathbf{\mathbf{r}}}_{n}$. These determine the relative binding probabilities for each of the many possible target actin sites. Combined with the above described step and stomp pathways, this gives us a complete kinetic description of myosin V informed by the polymer nature of the lever arms. Below we use the first passage times to compute various experimentally observable quantities including step distributions, the forward/backward step ratio, and the mean run length and velocity.
Our polymer model admits an accurate approximate analytical expression for the mean first passage time. In past theoretical work by Hinczewski et al. (2013), it was shown that the mean first passage time is approximately
where $a$ is the capture radius, $\mathcal{P}(\mathbf{\mathbf{r}})$ is the probability density of the free head being at position $\mathbf{\mathbf{r}}$, and ${D}_{\text{h}}=5.7\times {10}^{7}{\mathrm{cm}}^{2}/\mathrm{s}$ cm^{2}/s is the diffusion constant of the myosin V head which was estimated using the program HYDROPRO (Ortega et al., 2011) with the Protein Data Bank structure 1W8J (Coureux et al., 2004). Here $\mathcal{P}(\mathbf{\mathbf{r}})$ is the equilibrium distribution achieved at timescales greater than ${t}_{r}$, the relaxation timescale of the polymer. This result relies on the separation of timescales between the polymer relaxation of the lever arms and the first passage of the diffusive search. Previous BD simulations (Hinczewski et al., 2013) found that the lever arm relaxation timescale is ${t}_{\text{r}}\lesssim 5$ μs, which is two orders of magnitude smaller than the smallest mean first passage times, ${t}_{\text{fp}}^{\text{min}}\sim \mathcal{O}(0.1\text{}\mathrm{ms})$. Another requirement is that the time to diffuse the distance $a$, ${t}_{a}={a}^{2}/{D}_{\text{h}}\approx 2.8\mathrm{ns}$ ns be much smaller than the relaxation time ${t}_{\text{r}}$. This condition, which is clearly satisfied since ${t}_{a}/{t}_{r}\approx 5\times {10}^{4}\ll 1$, guarantees that after the free head reaches the capture radius it can undergo fast microscopic rearrangements required to bind without significant conformational changes in the rest of the lever arm.
In our model there is an additional condition necessary for binding: the free leg must be within the angular acceptance region with respect to the actin subunit. Therefore, we are actually interested in the mean first passage time to finding a binding site and simultaneously having the correct orientation. Electron and atomic force microscopy imaging of myosin V indicate that while bound the myosin head attaches approximately perpendicular to the outer surface of an actin subunit (Oke et al., 2010; Kodera et al., 2010). This binding, which involves the interaction of a specific region of the head with a corresponding region on the actin subunit, is mimicked in our coarsegrained model through an angular criterion: we require that the angle between the free leg and the outward pointing normal (${\widehat{\mathbf{\mathbf{r}}}}_{n}\times \widehat{\mathbf{\mathbf{z}}}$) at the target actin site be smaller than $\delta {\varphi}_{\text{ac}}$, which defines a conical acceptance region in which binding is allowed. Based on fits to experiments described below, we set $\delta {\varphi}_{\text{ac}}={55.6}^{\circ}$. We assume binding occurs with probability 1 when the free head is within a distance $a$ from a binding site and the free leg is inside the acceptance region.
Though this angular criterion is straightforward to implement in our BD simulations, for the analytical theory it significantly complicates the calculation. Hence we simplify the angular criterion in the analytical model, in a way that preserves most of the physical effect while making the derivation of the mean first passage time tractable. Instead of a conical acceptance region around the outward pointing normal to the target actin site, we only require an azimuthal angle similar to the normal. This amounts to replacing the probability density $\mathcal{P}({\mathbf{\mathbf{r}}}_{n})$ in Equation 2 with the joint density of finding the free head at position ${\mathbf{\mathbf{r}}}_{n}$ and the free leg simultaneously having azimuthal orientation similar to that of the binding site, $\mathcal{P}({\mathbf{\mathbf{r}}}_{n},\delta {\varphi}_{1}>({\varphi}_{f}+\pi ){\varphi}_{n}>\delta {\varphi}_{2})$. Here ${\varphi}_{f}$ is the azimuthal angle of the endtoend vector for the free myosin leg and ${\varphi}_{n}=12\pi n/13$ is the azimuthal angle of an actin binding site (note the factor of $\pi $ enters because the free leg points toward the binding site). The angles are shown in Figure 1B. The acceptance region defined by $(\delta {\varphi}_{1},\delta {\varphi}_{2})$ will for simplicity be taken as symmetric, $\delta {\varphi}_{1}=\delta {\varphi}_{2}=\delta {\varphi}_{\text{ac}}$, but we do not rule out the possibility that asymmetry in the myosin head or actin binding pocket favors binding from one direction. Similar effects have been previously observed for myosin V, for instance the suppression of ADP dissociation and head detachment depends asymmetrically on the direction of offaxis forces (Oguchi et al., 2010). The acceptance region requirement increases the first passage time, because in some cases the free head will diffuse to a binding site but have the wrong orientation and fail to bind. Therefore, the required separation of timescales still holds and, with a symmetric acceptance region, the mean first passage time is,
The use of the simplified angular criterion in the analytical theory gives results similar to the BD simulations, where the full conical acceptance region was used.
The properties of the polymer legs as well as the influence of the binding acceptance conditions are implemented through the joint distribution $\mathcal{P}({\mathbf{\mathbf{r}}}_{n},\delta {\varphi}_{\text{ac}}>({\varphi}_{f}+\pi ){\varphi}_{n})$. Adapting a polymer mean field theory (Hinczewski et al., 2013; Thirumalai and Ha, 1998) and exploiting the fact that the myosin lever arms are in the stiff regime ${l}_{p}\gg L$, we derive an approximate analytical expression for this distribution. We also derive an expression for $\mathcal{P}(\mathbf{\mathbf{r}})$ at any point $\mathbf{\mathbf{r}}$ in space, which describes the equilibrium distribution of the free leg during diffusion and captures the effects of the interleg structural constraint, the bound leg constraint, and the load force (Appendix 1). Assuming the binding events along with trailing/leading head detachment and hydrolysis are each Poisson processes with rates ${({t}_{\text{fp}}^{n})}^{1}$, ${t}_{\text{d1}}^{1}$, ${t}_{\text{d2}}^{1}$, and ${t}_{\text{h}}^{1}$ respectively, we derive the probabilities of passage through each kinetic pathway and analytical expressions for observable quantities including step distributions as well as mean run length and velocity.
Estimating the joint constraint strength from the free head spatial distribution during stepping
Until recently, the joint between the myosin V lever arms was believed to be freely rotating, allowing the detached head to perform a full threedimensional diffusive random walk while searching for actin binding sites. This hypothesis was supported by the work of Dunn and Spudich (2007), who used darkfield microscopy to image (from above) the diffusive path of 40 nm gold nanoparticles attached to the leg of myosin V. Their measured diffusion contours are radially symmetric about the interleg joint implying a freely rotating hinge. Other experiments also suggest a freely rotating joint (Shiroguchi and Kinosita, 2007; Fujita et al., 2012; Beausang et al., 2013), while theoretical and numerical models based on the free diffusion hypothesis largely agree with experimental data (Hinczewski et al., 2013; Craig and Linke, 2009; Mukherjee et al., 2017). Despite this body of evidence, recent experiments with improved spatiotemporal imaging resolution indicate that the myosin V joint might not freely rotate, instead preferring a particular interleg angle due to structural constraints. In particular, electron micrographs of Takagi et al. (2014) show that in the absence of actin, myosin V has a weakly preferred interleg angle of approximately 110°. Further, by tracking 20 nm gold particles attached to the Nterminus of the myosin head with interferometric scattering microscopy, Andrecka et al. (2015) observed a multipeaked diffusion contour, contradicting the presence of a freely rotating joint.
Our polymer model accounts for structural constraints on diffusion encoded through the interleg angle preference ${\theta}_{p}$ and constraint strength ${\mu}_{c}$. We let ${\theta}_{p}={83.0}^{\circ}$, as determined by fitting to step distribution and force dependence experiments discussed below. While this value is somewhat smaller that the actinfree Takagi et al. (2014) observations, we do not expect an exact correspondence. In particular, the observed angle will be larger than ${\theta}_{p}$ due to volume exclusion interactions and the effective preferred angle in the presence of actin may be further altered through the influence of the bound leg on the overall conformation of the protein. The presence of a joint constraint changes the diffusion of the free leg so that it swings through a nearly onedimensional arc while searching for binding sites rather than exploring the full threedimensional space. In order to visualize the diffusion, we use the probability density of the location of the free head $\mathcal{P}(\mathbf{\mathbf{r}})$ projected onto the two dimensional $zx$, $zy$, and $yx$ planes as well as the cylindrical plane $z\rho $, where $\rho ={x}^{2}+{y}^{2}$. The projection is performed by integrating over the remaining degree of freedom, for instance $\mathcal{P}(x,z)=\int \mathit{d}y\mathcal{P}(\mathbf{\mathbf{r}})$.
The diffusion contour measured by Andrecka et al. (2015) and shown in Figure 3I, corresponds to the $zx$ projection since they imaged through a glass surface parallel to both the myosin bound leg and direction of motion. We evaluated this contour using our analytical theory for several values of the interleg constraint strength, ${\mu}_{c}=0$, 3, 5, and 12 (Figure 3A–D). These calculations agree well with diffusion contours measured from Brownian dynamics simulations for the same values of ${\mu}_{c}$ (Figure 3E–H). Our results show that a nonzero interleg constraint ${\mu}_{c}>0$ is required to produce a multipeaked contour similar to that measured by Andrecka et al. (2015). In particular, when ${\mu}_{c}=5$ the predicted diffusion agrees well with the experimental data, while ${\mu}_{c}=3$ and ${\mu}_{c}=12$ respectively produce contours with lower and higher peaks than the measured distribution. Therefore, we estimate the energy of the interleg constraint lies within this range $3{K}_{B}T\lesssim {\mu}_{c}{K}_{B}K\lesssim 12{K}_{B}T$. Incidentally, the upper bound 12 k_{B}T is also approximately the energy at which our model can no longer accurately predict the multipeaked step distribution of 8IQ mutants discussed in the next section.
For the cases where the constraint is strong enough to induce a bimodal distribution, the two peaks are roughly equal in height both in the analytical theory and in simulations. In constrast, the experimental distribution is asymmetric, with the peak closer to the actin about 1.5 times higher than the one further away from the actin. While we cannot rule out that some aspect of this asymmetry may be due to additional structural features not present in the model, it is also likely to be at least in part an artifact of the finite time resolution of the experiment and possibly consequences of attaching a gold nanoparticle to infer the motility of the detached motor head. The peak region further away from actin is the one that is initially visited after detachment of the trailing head. If one cannot resolve the precise moment of detachment, one might not be able to fully capture this portion of the distribution in the experimental data. In fact, Andrecka et al. (2015) measured the distribution for a range of ATP concentrations, and as the ATP was progressively increased from 1µM to 1 mM, the height of the far peak decreased relative to the near peak, even though the distribution should in principle be independent of ATP. Detachment is much more rapid at high ATP, however, and the resulting practical difficulties in collecting trajectories were mentioned in the paper.
The other projections complete the picture of the free head diffusion. These are shown for the constrained diffusion model (${\mu}_{c}=5$) in Figure 3—figure supplement 1 and for the free diffusion model (${\mu}_{c}=0$) in Figure 3—figure supplement 2. When the interleg constraint is present, the $yx$ (front view) and $zy$ (top view) projections together show that during its diffusive search the free leg swings out away from the actin axis at a particular height, with motion similar to that of a drawing compass. We also computed and simulated diffusion contours for myosin V under 2 pN backward force (Figure 3—figure supplement 3) as well as for the 4IQ ((Figure 3—figure supplement 4) and 8IQ mutants (Figure 3—figure supplement 5). Under force, the free head distribution rotates toward the minus end of the actin changing which target binding sites are most likely to be found during a step. The diffusion of mutants is almost identical to wildtype myosin V, but is scaled up or down based on lever arm length. This suggests that qualitative differences in stepping between mutants (discussed below) are due entirely to actin geometry.
Our calculation of the free head spatial distribution allows us to test how different forms of the interleg potential affect the diffusive step of myosin V. Figure 3—figure supplement 6 shows the KullbackLeibler (KL) divergence ${D}_{KL}({\mathcal{P}}_{\mathrm{cos}}\mathcal{P})=\int {d}^{3}\mathbf{\mathbf{r}}{\mathcal{P}}_{\mathrm{cos}}(\mathbf{\mathbf{r}})\mathrm{log}\left[{\mathcal{P}}_{\mathrm{cos}}(\mathbf{\mathbf{r}})/\mathcal{P}(\mathbf{\mathbf{r}})\right]$ between the free head distribution ${\mathcal{P}}_{\mathrm{cos}}(\mathbf{\mathbf{r}})$ arising from the cosine joint potential ${\mathcal{\mathscr{H}}}_{J}={\mu}_{c}{k}_{B}T[1\mathrm{cos}({\theta}_{J}{\theta}_{p})]$ introduced above, and the distribution $\mathcal{P}(\mathbf{\mathbf{r}})$ arising from a general quartic interleg potential ${\mathcal{\mathscr{H}}}_{J}={\mu}_{c}{k}_{B}T[(\mathrm{\Delta}\theta {)}^{2}/2+{h}_{3}(\mathrm{\Delta}\theta {)}^{3}/3!+{h}_{4}(\mathrm{\Delta}\theta {)}^{4}/4!]$, where $\mathrm{\Delta}\theta ={\theta}_{J}{\theta}_{p}$. Also shown in this figure are the $xz$ diffusion contours for a few representative alternative potentials. This calculation provides a quantitative picture of how much the diffusion changes over a range of joint potentials. Specifically, our results are insensitive to the form of the potential when the anharmonic terms (${h}_{3}$, ${h}_{4}$, etc.) are small enough to not introduce a new energy minimum in the physically relevant range of interleg joint angles ($\mathrm{\Delta}\theta \lesssim \pi /2$). In this case, the resulting diffusion contour is bimodal and qualitatively similar to that shown in Figure 3.
The imaging techniques used by Andrecka et al. (2015) require a glass coverslip which excludes half the space available for the myosin diffusive search. They report no evidence of interactions between the myosin and the surface. However, entropic forces due to volume exclusion do influence the diffusion. While it is unlikely these effects would induce a multipeaked diffusion contour, we explicitly check this by adding a potential barrier excluding half of space in Brownian dynamics simulations (see Figure 3—figure supplement 7). Our results rule out the appearance of a multipeaked distribution due solely to volume exclusion and confirm the coverslip only slightly changes the shape of the diffusion contours.
The discrepancy between the singlepeaked distribution of the Dunn and Spudich (2007) experiment and the the multipeaked distribution of Andrecka et al. (2015) remains a controversy, but as noted in the latter work it can be partially resolved by adding localization noise and rebinning their data, accounting for the differences in gold nanoparticle size and the associated measurement precision between the two experiments. The fact that Dunn and Spudich (2007) attached their gold nanoparticle labels on the lever arm closer to the joint may also have contributed to concealing the constrained diffusion: the closer the label is to the joint, the closer together the two peaks appear in the distribution, and the harder it becomes to distinguish one from two peaks given finite experimental spatial resolution.
In the following sections, we use our model to test the constrained diffusion hypothesis against the large body of existing experimental data on myosin V. In the Discussion, we compare the constrained and free diffusion models and suggest further experiments to more conclusively discern between these competing myosin V diffusion hypotheses.
Constrained diffusion model predicts zero force step distributions
The spatial fluctuations of motor proteins, among them distributions of step sizes, have long been an intriguing aspect of their behavior (Das and Kolomeisky, 2008). Let us first consider the step size distribution under zero backward force, one of the most commonly measured aspects of myosin V dynamics. The classic experiment by Yildiz et al. (2003) measured one of the first myosin step distributions using fluorescence imaging with one nanometer accuracy (FIONA). Later experiments using FIONA (Sakamoto et al., 2005) and electron microscopy (Oke et al., 2010) determined step distributions for both wildtype myosin as well as 4IQ and 8IQ mutants. The FIONA experiments measure the distance travelled by the trailing head while executing a forward step, while the electron micrographs were used to directly determine the number of actin sites separating legs while bound. These experiments qualitatively agree, with the wildtype distribution peaked at the halfhelical site ($z$ = 36 nm) and an approximately linear relation between step size and lever arm length, though Sakamoto et al. (2005) measured slightly shorter 4IQ steps than Oke et al. (2010). Furthermore, Sakamoto et al. (2005) found that, due to the longer leg length, 8IQ mutants can reach the full helical actin sites ($z$ = 72 nm), giving rise to multipeaked step distributions. This effect was not observed by Oke et al. (2010), which they note is likely because their image processing technique did not resolve large steps for which the myosin is stretched out very close to the actin.
Our polymer model gives the step distributions in terms of the structural parameters and kinetic rates described in the preceding sections. As shown in Appendix 1, after trailing leg detachment the probability of a step/stomp to site n is
where ${r}_{\text{T}}={\sum}_{n}{b}_{n}{({t}_{\text{fp}}^{n})}^{1}$. The step/stomp probabilities for the leading leg ${\mathcal{P}}_{\text{L}}^{n}$ are given by the same expression with the substitutions, ${b}_{n}\to {b}_{n}$, ${t}_{\text{h}}\to 0$, and ${r}_{\text{T}}\to {r}_{\text{L}}={\sum}_{n}{b}_{n}{({t}_{\text{fp}}^{n})}^{1}$ to account for the lack of ATP hydrolysis and recovery stroke in these pathways. The trailing leg detachment occurs with probability $g{(1+g)}^{1}$ while the leading leg detachment occurs with probability ${(1+g)}^{1}$. Thus, the probability of observing a bound myosin with $n$ subunits between the leading and trailing heads is
This distribution (normalized to remove termination pathways), describes the probability of finding $n$ actin sites between bound heads, which was measured in the Oke et al. (2010) experiment. The overall step sizes, that is the distances traveled by the trailing leg in executing a forward step, have a distribution given by the convolution ${\mathcal{P}}_{\text{T}}({z}_{n})=({\mathcal{P}}_{\text{dist}}*{\mathcal{P}}_{\text{T}})[n]$, where ${z}_{n}=n\mathrm{\Delta}z/2$. The distribution ${\mathcal{P}}_{\text{T}}({z}_{n})$ includes both steps and stomps. For comparison to forward step distributions measured in experiments, we exclude the stomps by disallowing binding behind the bound leg and setting ${\mathcal{P}}_{\text{T}}^{n}=0$ for $n<0$ in the convolution. To account for the measurement resolution in FIONA experiments, we convolve this distribution with 1 nm Gaussian noise and bin the data identically to Sakamoto et al. (2005).
To fit step distributions, we fix the binding penalty $b$, capture radius $a$, and power stroke effectiveness defined as $\mathcal{T}=1+20{\nu}_{v}/(20+7\kappa {\nu}_{c})$, where $\kappa =L/{l}_{p}$, determined via fitting to force dependence data (described below). The variable $\mathcal{T}$ measures the energy loaded by the power stroke in the effective spring formed by the lever arms (Hinczewski et al., 2013). We vary the bound leg constraint angle ${\theta}_{c}$, the acceptance region size $\delta {\varphi}_{\text{ac}}$, and interleg preferred angle ${\theta}_{p}$ to minimize the KullbackLeibler (KL) divergence between experimental step distributions and theoretical predictions for 4, 6, and 8 IQ myosin V. We also let the persistence length ${l}_{p}$ and bound leg constraint strength ${\nu}_{c}$ vary along curves of constant $\mathcal{T}$, but the model predictions are robust to such parameter variation. Informed by these fits, we choose a set of parameters for which agreement between theory and experiment is excellent for all distributions. We further refine the parameter choices by alternating between fitting step distributions and force dependence behavior, feeding the parameters determined by one fit into the other. We find that ${\theta}_{c}={65.0}^{\circ}$, $\delta {\varphi}_{\text{ac}}={55.6}^{\circ}$ and ${\theta}_{p}={83.0}^{\circ}$ optimizes the agreement between our theory and the full set of experimental step distributions. These and all the other parameter values in the model are summarized in Table 2. This constraint angle is similar to that indicated by imaging experiments (Walker et al., 2000; Kodera et al., 2010; Lewis et al., 2012), while the preferred interleg angle is slightly smaller than that observed for unbound myosin V (Takagi et al., 2014), which, as noted above, could be due to volume exclusion effects and the influence of binding to actin. The azimuthal constraint $\delta {\varphi}_{\text{ac}}$ corresponds to an acceptance region of ∼111°, which plays a key role in restricting binding too far from the halfhelical and fullhelical actin sites. Step distributions and direct measurements of the azimuthal trajectory of myosin indicate steps with large azimuthal rotations are rare (Lewis et al., 2012), likely due to orientational binding constraints.
The theoretical head separation distributions are shown in Figure 4, row 1 alongside data from Oke et al. (2010) while in Figure 4, row 2 we plot the full convolved step distributions. Since Oke et al. (2010) were not able to resolve large head separations for 8IQ mutants, we compute an alternative ‘small steps’ distribution in which large head separation (> 18 actin subunits) is not considered. The theoretical distribution peaks and width agree well with Oke et al. (2010), especially for the 4IQ mutant. For 6IQ and 8IQ myosin V, the distributions are slightly broader, more similar to Sakamoto et al. (2005). In particular, we capture the multipeaked 8IQ distribution, with peaks at about 78 nm, 110 nm, and 135 nm (indicated by arrows in panel F). The first peak comes from steps in which the free leg detaches from the halfhelical site behind the bound leg and steps to the halfhelical site ahead of the bound leg. Similarly, the second peak corresponds to steps with one headtohead distance (during detachment or binding) being the full helical length and the other being the halfhelical length. Finally, the third, smallest peak corresponds to the free leg starting and ending at the full helical distance. One might expect the rich stepping behavior of the 8IQ mutant to be inconsistent with the constrained diffusion hypothesis. If the joint does not freely rotate, can the myosin reach both the half and fullhelical actin subunits? Our model shows that a 5 K_{B}T joint constraint energy is sufficiently weak to allow multipeaked step distributions, but strong enough to considerably alter the diffusion (see previous section). We also tested large constraint energies and found that for ${\mu}_{c}{k}_{B}T\approx 1012{k}_{B}T$ we could no longer simultaneously fit the multipeaked 8IQ step distribution and other step distributions. This upper bound on the joint constraint energy is similar to that estimated from diffusion contours.
The BD simulation step distributions are consistent with the results of the analytical model, as seen in Figure 4. Additionally, the simulations allow us to directly visualize examples of individual stepping trajectories. Figure 4—video 1 through Figure 4—video 3 show steps of three different sizes for the 6IQ wildtype without the interleg constraints (${\mu}_{c}=0$). Figure 4—video 4 through Figure 4—video 6 are analogous, but with the constraint present (${\mu}_{c}=5$).
Myosin V exhibits robust step distributions under load
As a backward load force is applied to myosin V, forward steps decrease slightly in size while trailing leg stomps and backward steps become more likely. To elucidate this transition, we derive the combined leading and trailing leg step distribution (including stomps). To begin, in addition to the trailing leg steps considered above, we must also consider the distributions of steps and stomps originating from leading leg detachment. In this case, the bound leg position is behind the free leg starting position. So the initial distance between the heads is distributed as ${\mathcal{P}}_{\text{dist}}^{n}={\mathcal{P}}_{\text{dist}}^{n}$. The leading leg step distribution is then given by the convolution ${\mathcal{P}}_{\text{L}}({z}_{n})=({\mathcal{P}}_{\text{dist}}*{\mathcal{P}}_{\text{L}})[n]$, with ${z}_{n}=n\mathrm{\Delta}z/2$, and the combined leading/trailing leg step distribution is,
where ${\mathcal{P}}_{\text{T}}({z}_{n})$ is as defined in the preceding section. The force dependence of this distribution comes from the influence of the force on the myosin lever arm, which bends under the load giving rise to a new effective constraint direction ${\widehat{\mathbf{\mathbf{u}}}}_{c}^{\prime}$ and new power stroke effectiveness ${\mathcal{T}}^{\prime}$,
where $F$ is the magnitude of the load force with direction $\widehat{\mathbf{\mathbf{F}}}$, $L$ is the myosin leg length, and $\beta =1/({k}_{B}T)$. For backward forces, $\widehat{\mathbf{\mathbf{F}}}=\widehat{\mathbf{\mathbf{z}}}$, the lever arm bends toward the minus end of the actin increasing the effective constraint angle, but decreasing the power stroke effectiveness.
In Figure 5A, we show the combined distribution Equation 6 for substall backward forces, $F$ = 0.0 and 1.0 pN, at stall force, $F$ = 1.9 pN, and above stall force, $F$ = 2.5 pN. In this plot, peaks near +72 nm are forward steps, peaks near 0 nm are trailing or leading stomps, and peaks near –72 nm are backward steps. At zero force forward steps are dominant, exhibiting the experimentally observed step distribution, while forward leg stomps also occur less frequently. The myosin remains resilient as the force increases to 1 pN, with the step distribution shifting backward by approximately a single actin binding site. This result is consistent with the step distribution under 1 pN backward force measured by Clemen et al. (2005), which was nearly identical to zero force step distributions.
Further increasing the force to stall at 1.9 pN, we see that the probability of forward steps dramatically decreases. Interestingly, even though forward steps are no longer kinetically dominant, the location and shape of their distribution remains robust, as shown in Figure 5B (BD simulations show the same qualitative behavior, see Figure 5—figure supplement 1). The stall force is defined to be the force at which the expected run length of the myosin is zero. Our calculation shows that stall is due to the emergence of a backward step distribution identical to the that for forward steps. Thus, at stall force the myosin takes forward and backward steps with equal probability, making no progress along the actin. In this regime, however, the dominant kinetic pathway is trailing leg stomps which give rise to the large peak around 0 in the step/stomp distribution. Finally, increasing the force further eliminates all forward steps, increasing the prevalence of trailing stomps while shifting the backward step distribution by about one actin binding site.
An interesting prediction made by the model is the distribution of stomps, which is particularly apparent above stall force and has comparable width to the zero force step distribution. With ever improving imaging technology, it may soon be possible to reliably detect such variation in stomping experimentally. Because stomps are the dominant kinetic pathway in this force regime, many stomping events could be observed if sufficient imaging precision is achieved. Such measurements would provide further insights into how myosin V executes foot stomps. If the measured distribution significantly differs from our theoretical prediction, perhaps alternative mechanisms are involved in foot stomping. Since the leading head retains its ADP molecule after unbinding, it can in principle rebind quickly. Therefore, it is possible that leading foot stomps are dependent on nonequilibrium diffusion that occurs before the myosin legs relax to equilibrium. Further experiments are required to complete our understanding of myosin foot stomping.
The predictions made in Figure 5 will require optical trap experiments in which a load force is applied to a cargo. If the contours of the distribution of the free head are to be simultaneously measured then it becomes necessary to also attach a gold nanoparticle (Andrecka et al., 2015). The analyses of the results of such experiments, which could test our predictions, might be complicated due to hydrodynamics effects of both the cargo and the gold nanoparticle.
Probing the biological function of the joint constraint: effects on timing and consistency of stepping
In addition to making behavioral predictions, our model also provides insights into the potential biological function of the structural joint constraint. As discussed above, our computed diffusion contours and the measurements of Andrecka et al. (2015) indicate that the joint constraint reduces the diffusion search space from fully threedimensional to nearly onedimensional, with the free head locations concentrated along a curve. How does this effect influence the stepping behavior of myosin V? One might expect that restricting the search space would affect both mean binding times and the step distributions.
Let us first consider the subject of binding times. As discussed above, binding in the model is more complex than just waiting for the head to diffuse within the capture radius $a$ of a binding site on actin. For the trailing leg, successful binding cannot occur until the head has hydrolyzed ATP and undergone the recovery stroke. Thus the mean timescale for ATP hydrolysis, ${t}_{h}$, sets the lower bound on the mean binding time. After the recovery stroke, if the trailing leg diffuses near one of the backward ($n<0$) binding sites, the probability of binding is reduced by a factor $b<1$ because the postrecoverystroke head is not in an orientation that favors backward binding. Since polymer relaxation is fast between capture attempts, in effect this approximates partially absorptive reaction kinetics (Šolc and Stockmayer, 1971): the postrecoverystroke head must make on average $1/b$ diffusive excursions near an $n<0$ site before it is captured. The net result of these two constraints, as derived in Appendix 1, is that the mean binding time for the trailing leg is:
Here, ${b}_{n}=b$ for $n<0$, ${b}_{n}=1$ for $n>0$. In contrast, if one is interested in the mean time for the head to diffuse within a distance $a$ of any of the binding sites, this is given by ${t}_{\text{diff}}={\left[{\sum}_{n}{({t}_{\text{fp}}^{n})}^{1}\right]}^{1}$ (irrespective of whether it is a trailing or leading head). Figure 6 shows ${t}_{T}$ versus ${t}_{\text{diff}}$ as a function of applied force F, and we always see ${t}_{T}>{t}_{\text{diff}}$, as expected from the above constraints. For small $F$ the trailing leg kinetics is dominated by forward steps. ${t}_{\text{diff}}<{t}_{h}$ in this regime, but even though the head can diffuse rapidly to forward binding sites, it does not bind until hydrolysis occurs. As the force is increased toward the stall regime, trailing leg stomps become the dominant pathway. The mean binding time ${t}_{T}$ increases by an order of magnitude, from about 1.6 ms at $F=0$ to a peak of 25 ms at $F=1.8$ pN. This is because the force biases the system toward backward binding sites (trailing leg stomps), but binding to those sites requires multiple diffusive attempts even after hydrolysis has occurred, because of the head orientation.
This increase may be reflected in an interesting experimental observation. Veigel et al. (2002) used an actin filament attached to beads in an optical trap, and studied the interactions between individual myosin V motors and the filament. By sinusoidally oscillating one bead (and hence the filament) they could use the amplitude of the resulting fluctuations to determine the stiffness of the motoractin complex over time. Intervals of reduced stiffness were interpreted as times when only one head was bound, and hence the unbound head is undergoing a diffusive search. At low loads ($F<1$ pN) these intervals averaged around ∼18 ms, while at nearstall loads the intervals lasted for hundreds of milliseconds. This order of magnitude relative increase is consistent with our prediction for ${t}_{T}$ as a function of $F$, with the interpretation that what they were observing at high loads was primarily trailing head stomping. The absolute timescales in the experiment are larger than in the theory (for example the predicted $t}_{T$ = 1.6 ms at $F=0$ versus the experimental value of ∼18 ms), but this is likely due to the fact that the sinusoidal forcing (at 75 Hz with peaktopeak amplitude 250 nm) is a substantial perturbation to the system that makes it harder for the head to bind (hence increasing mean binding times). However the relative increase of the binding times with force is captured by our model. Other experiments, where a gold nanoparticle is attached to one of the motor legs (Dunn and Spudich, 2007; Andrecka et al., 2015), have found binding timescales on the order of tens of milliseconds in the low load regime. Here there is another experimental artifact in play: as discussed in Hinczewski et al. (2013), the drag from the gold nanoparticle can substantially increase diffusion times, with ${t}_{T}$ becoming much larger than ${t}_{h}$ at $F=0$ because of slow diffusion.
To complete the description of the binding times, we note that the situation after leading leg detachment has several differences. Because the ADP molecule is retained, ATP hydrolysis is not necessary for rebinding, and the post power stroke orientation of the head favors backward rather than forward sites (with the latter penalized by a factor of $b$). The resulting leading leg mean binding time ${t}_{L}$ is:
At small loads the primary pathway for the leading leg is stomping. Since the head orientation necessitates multiple diffusive attempts to rebind for leading stomps, ${t}_{L}$ is much larger than ${t}_{\text{diff}}$ for small $F$, as seen in Figure 6. At larger loads backward stepping becomes the dominant pathway after leading leg detachment, and the combination of the backwards force and the favorable head orientation makes backward stepping nearly diffusionlimited, with ${t}_{L}$ approaching ${t}_{\text{diff}}$.
Based on the above discussion of binding times, we can guess that the effect of the structural constraint ${\mu}_{c}$ might be relatively limited: while ${\mu}_{c}$ does indeed affect the diffusive search space, for forward stepping at low loads the binding timescale is dominated by the ATP hydrolysis time ${t}_{h}$. Despite the fact that increasing ${\mu}_{c}$ restricts the search space (and hence potentially makes diffusion faster), we find that ${t}_{\text{T}}$ at $F=0$ does not substantially decrease and can increase depending on the preferred interleg angle (Figure 7, dashed lines). While adding the constraint makes it easier to find a few target binding sites, it simultaneously makes it harder to find other sites, so that the overall change in binding time is small. Furthermore, we emphasize that head detachment (not the binding time) is the rate limiting action in the stepping cycle at low loads. Thus, the mean binding time only influences the motor dynamics indirectly, through determination of the mean number of steps taken before complete dissociation. To have any impact on motor performance, the mean binding time would have to decrease more substantially than we have found here.
So if the timing of the steps is not significantly affected, what about the step distributions? We computed these for several values of ${\mu}_{c}$ and found that increasing the joint constraint energy narrowed the distributions. Figure 7 (solid lines) shows the standard deviation of the forward step distribution ${\sigma}_{\text{step}}$ as a function of ${\mu}_{c}$ for preferred interleg angle ${\theta}_{p}={83}^{\circ}$ (the value used in above fits), 65°, and 95°. The standard deviation decreases, independent of angle preference, by approximately 0.1 nm per ${k}_{B}T$ of constraint energy. The constraint therefore plays a role in improving the consistency with which steps are made to a particular few target actin binding sites near the half helical length.
Finally, we conjecture that the joint constraint may also play a role in obstacle avoidance inside the crowded confines of a cell. Our diffusion contours (Figure 3, Figure 3—figure supplement 1, Figure 3—figure supplement 2) indicate that in the constrained diffusion model, the free head is significantly closer to the actin filaments throughout the diffusion. We believe that crowding could further restrict the conformational space explored by the free head, thus resulting in an even narrower step size distribution in vivo. Myosin V procession in the presence of obstacles would be an interesting topic for future study.
Load dependence: step ratio, run length, and run velocity
For our final test of the constrained diffusion model we compare theoretical predictions to experiments on the load dependence of the forwardbackward step ratio, the mean run length, and the mean run velocity (Figure 8). To fit the load dependence data, we require that the stall force and zero force run length agree with experimentally determined values, as further described below. This is achieved by varying the power stroke effectiveness $\mathcal{T}$, the binding penalty $b$, and capture radius $a$. All other parameters were held fixed as determined by experimental data or fitting the step distributions. As noted previously, we alternate between fitting the step distributions and force dependence to achieve optimal agreement between the theory and experiments.
Step ratio
The backwardtoforward step ratio is defined as ${\mathcal{P}}_{b}/{\mathcal{P}}_{f}$, where ${\mathcal{P}}_{b}=1/(1+g){\sum}_{n<0}{\mathcal{P}}_{L}^{n}$ is the probability of taking a backward step and ${\mathcal{P}}_{f}=g/(1+g){\sum}_{n>0}{\mathcal{P}}_{T}^{n}$ is the probability of taking a forward step. The step ratio measured by Kad et al. (2008) exponentially increases for large load force as the myosin transitions from forward to backward stepping. The step ratio ${\mathcal{P}}_{b}/{\mathcal{P}}_{f}=1$ at $F=1.9\mathrm{pN}$ pN, when forward and backward steps are equally likely.
For the full range of physically reasonable parameters, we find that the force at which the step ratio is one is nearly identical to the stall force (defined as the force for which the run length ${z}_{\text{run}}=0$). Therefore, as part of our fitting process, we require the stall force be F_{stall} = 1.9 pN. This value is within the range measured in other experiments. The agreement between the stall force and the force at which the step ratio is 1, is a consequence of the symmetry between the forward and backward step distributions shown in Figure 5A.
Interestingly, for intermediate forces beginning around F ≈ 1 pN, there is a small increase in the step ratio before the exponential divergence near the stall force. We find that this effect is due to the binding penalty and only occurs for $b\lesssim 0.1$. Taking $b=0.045$, the step ratio predicted by the model agrees well with the experimental data as shown in Figure 8A.
Run length
The run length is the distance travelled by the myosin V along the actin in a given run. By averaging over the step distribution (Equation 6), we compute the mean run length,
where ${\mathcal{P}}_{t}$ is the termination probability and ${z}_{n}=n\mathrm{\Delta}z/2$ as above. In this expression the numerator is the mean step size while the denominator accounts for the mean number of steps occurring in a run and normalization of the distribution $\mathcal{P}({z}_{n})$. In Appendix 1 we show that the mean number of steps is $\u27e8{N}_{\text{run}}\u27e9=1/{\mathcal{P}}_{t}$. Finally, the factor of 1/2 accounts for the fact that the center of mass of the myosin moves half the distance of the head domain in a given step.
Experimental estimates, mostly performed with zero load force, report mean run lengths in the range 0.7 – 2.4 μm , with the large variation likely due to measurement conditions (Sakamoto et al., 2000; Baker et al., 2004; Pierobon et al., 2009). We choose a representative value ${z}_{\mathrm{r}\mathrm{u}\mathrm{n}}=1.3$ μm for fitting the model. With this choice, the computed run length under force also is in reasonable agreement with measurements by Clemen et al. (2005) as shown in Figure 8B.
Run velocity
The mean run velocity is ${v}_{\text{run}}={z}_{\text{run}}/{t}_{\text{run}}$, where ${t}_{\text{run}}$ is the mean time elapsed during a processive run. In Appendix 1, we find
where the mean trailing/leading binding times ${t}_{T}$ and ${t}_{L}$ are given by Equations 8 and 9. The other timescale in this expression ${t}_{\text{wait}}=g{t}_{d1}/(1+g)$ is the mean waiting time between detachment events. Each time is weighted by the respective probability of occurrence and the mean number of steps taken.
The theoretical mean run velocity agrees well with the experimental data shown in Figure 8C. In addition to fitting the drop to zero run length at the stall force, our model also captures the gradual decrease in run velocity at low forces ($\lesssim 1\mathrm{p}\mathrm{N}$), which is due to small backward shifts in the step distribution under load. This effect is not seen in the previous simplified model in which only halfhelical steps were considered. We also compare to the large force measurements of Gebhardt et al. (2006), which observed myosin V walking backward with velocity ≈ 90 nm/s under 3 pN backward force. This experiment was performed at 1 μM ATP concentration, which lowers the ATP binding rate and associated trailing leg dissociation rate. Therefore, for this comparison we set t_{d1} = 2.2 s^{—1} as estimated from experimental kinetics (Hinczewski et al., 2013; Gebhardt et al., 2006). The low ATP run velocity is shown as the dashed curve in Figure 8C. The backward velocity under 3 pN force is comparable, though slightly smaller than the measured value. This small discrepancy, as well as disagreement with high force data (at 5 and 10 pN force), is likely due to additional kinetic pathways in the superstall regime, such as the power stroke reversal discussed above.
Robust procession under off axis loads
While we have focused on verifying the consistency of structurally constrained diffusion with the full range of experiments, our model also provides a fast analytical framework for studying complex behaviors of myosin V. As an example, we consider procession under offaxis load forces, which are not parallel to the axis of the actin filaments. Offaxis loads are likely especially relevant in vivo where the myosin navigates a crowded environment.
Experiments by Oguchi et al. (2010) measured the distance myosin V walked against offaxis loads applied by a bead in an optical trap. To apply the offaxis force, the actin was attached to a glass stage which was shifted perpendicularly to the actin axis after the myosin had taken two steps. For forces that were approximately ±25° offaxis on average, the myosin walked slightly further (∼10 nm on average) than under backward force. The force at termination as well as the onaxis component of the termination force were also larger for offaxis loads. We emulated these experiments using Monte Carlo simulations, adjusting the load force after each step or stomp based on the distance from the starting location and shifting the force off axis after two steps. The average simulated run length and termination force each increased under offaxis loads (not shown), agreeing qualitatively with experiments.
We also directly compute the mean run length under a constant 1 pN load using Equation 10 for a large range of offaxis forces parameterized by the spherical angles ${\theta}_{\text{F}}$ and ${\varphi}_{\text{F}}$. The percent difference between the run length with offaxis force $({\theta}_{\text{F}},{\varphi}_{\text{F}})$ and the run length under backward force is shown in Figure 9. Our model estimates that at worst, changing the direction of a constant load force can decrease the run length by ∼15%, while many offaxis directions lead to considerable increases in run length. The shortest run length for 1 pN force occurs at $({\varphi}_{\text{F}}=0$, ${\theta}_{\text{F}}\approx {20}^{\circ})$. This particular ${\theta}_{\text{F}}$ is slightly larger than that (${\theta}_{\text{F}}\approx {6}^{\circ}$) which maximizes the effective constraint angle ${\theta}_{c}^{\prime}=\mathrm{arctan}(\widehat{\mathbf{\mathbf{x}}}\cdot {\widehat{\mathbf{\mathbf{u}}}}_{c}^{\prime}/\widehat{\mathbf{\mathbf{z}}}\cdot {\widehat{\mathbf{\mathbf{u}}}}_{c}^{\prime})$, but has a larger ${\mathcal{T}}^{\prime}$, which strengthens the preference for locations along the constraint direction so that diffusion to forward actin sites becomes more difficult. The competition between these two effects leads to the minimized run length. For fully offaxis force (${\varphi}_{\text{F}}=\pm {90}^{\circ}$, ${\theta}_{\text{F}}>0$), all run lengths are larger than those under backward force. These results further corroborate the robust processivity of myosin V under offaxis loads observed experimentally by Oguchi et al. (2010).
Discussion
Constrained versus free diffusion: simplified effective theories
Our theoretical analysis presented above largely supports the constrained diffusion hypothesis. The $zx$ diffusion contour produced by the constrained diffusion model closely resembles that measured by Andrecka et al. (2015) and we find strong quantitative agreement between the model predictions and experimental data for zero force step distributions and force dependence of the step ratio, run length, and run velocity. While this evidence is convincing, we emphasize that with relatively minor parameters changes the free diffusion model also makes similarly accurate predictions, with the exception of the diffusion contours. As is particularly apparent from comparison of the $zx$ and $yx$ projections (Figure 3—figure supplement 2), free diffusion doesn’t produce multipeaked distributions, implying that the freely rotating myosin explores the entire three dimensional space, contrary to the results of Andrecka et al. (2015). With respect to other experiments, however, the free diffusion model performs favorably, agreeing with step distribution and force dependence data nearly as well as the constrained diffusion model (see Figure 4—figure supplement 1 and Figure 8—figure supplement 1).
The similarity between the predictions of the constrained and free diffusion models can be explained through further analysis of the diffusion contours close to the actin. While the global diffusion patterns are very different from one another, close to the actin they look quite similar. Therefore, both the onedimensional path under constrained diffusion and the full threedimensional exploration under free diffusion favor binding to similar actin sites, namely those near the halfhelical length of actin. This means that the free diffusion model accurately describes both step distributions and force dependence behavior. Since the diffusion time scale is much smaller than the head detachment time, the run velocity as a function of force is also quite similar between models.
The similarity between the free and constrained diffusion model predictions for onactin behavior indicates that experimental measurements of such quantities (step distributions, run length, etc.) cannot effectively discern the true structure of the myosin V joint. The diffusive search, however, is dramatically influenced even by a small joint structural constraint. Therefore, to further probe the presence and consequences of structural constraints, experiments should directly measure the free leg diffusion. In particular, the $yx$ projection of the diffusion shows particularly stark contrast between the constrained and free diffusion models. This contour corresponds to imaging the myosin from the front along the axis of the actin.
While our analysis and the most recent experiments suggest that myosin V has structurally constrained diffusion, the free diffusion model is still useful as a simplified effective theory for stepping behavior on the actin filaments. The precise values of physical parameters likely cannot be accurately determined with the free diffusion model, but extrapolation from fits to existing experimental data allows for valuable behavioral predictions. Both the free and constrained diffusion models can contribute to our understanding of the function of myosin V in complicated environments, for instance with offaxis load forces (see above). In fact, previous work (Hinczewski et al., 2013), indicates that the force dependence of average quantities, such as the forward step probability, step ratio, or run length, are well described by a further simplified model that only allows steps to the halfhelical actin sites. Such a model cannot describe step distributions, but does contribute to our understanding of myosin V’s resilience under backward force and robust motility under perturbations to structural parameters. Since our computed and simulated diffusion contours indicate a similarity between the constrained diffusion model and simplified free diffusion models near the actin, the latter remain quite useful for studying certain behaviors of myosin V.
Conclusion
We have developed a comprehensive lowforce model of myosin V, incorporating the polymeric nature of the lever arms, the joint angle preference which gives rise to a structurally constrained diffusive search, and the full set of kinetic pathways involving all actin binding sites. The analytical model allows us to compute bounds on the joint constraint energy and captures experimental results for step distributions and the force dependence of the step ratio, run length, and run velocity. While our results are largely in support of the constrained diffusion hypothesis, the theory also provides insight about how simplified models (eg. free diffusion) can provide a useful analytical description of some myosin V behaviors. Finally, using the model we can make predictions about experimentally measurable quantities including stomp distributions and robust run length under offaxis forces. To conclude we discuss limitations and potentially interesting extensions to the model for future studies.
Throughout the present study, we have assumed that myosin V walks along a single static actin filament. In reality, the many actin filaments within a cell can come together forming crossing and branching network structures (Pollard et al., 2000). Furthermore, bending and rotational fluctuations of the actin also occur and have been characterized experimentally (Egelman et al., 1982). A recent study observed myosin V walking on actin rafts (Bao et al., 2013), while fluctuations have been considered in previous theoretical work on step distributions (Vilfan, 2005b). New actin geometries and fluctuations could be easily incorporated in our model by convolving the distribution of actin binding site locations with the distribution for the position of the free myosin head.
One regime in which our model fails is under very large superstall loads. In particular, the large backward run velocity measured by Gebhardt et al. (2006) at 5 pN and 10 pN backward force are not captured by the model. This alternative behavior is likely due to the power stroke reversal noted above and observed in experiments (Sellers and Veigel, 2010), which would promote more frequent large backward steps. It is also likely that under such extreme forces the detachment rates of the myosin heads are altered. By adding new power stroke reversal kinetic pathways, we expect our model to better capture the large force behavior of myosin V.
Finally, our methods may be more broadly applied to the large class of processive molecular motors that operate through the combination of chemical reactions and a diffusive Brownian search. Of particular interest are myosins VI and X, which operate through a similar set of kinetic pathways, but have heterogenous lever arms with both stiff and flexible components (Sun and Goldman, 2011). Perhaps applying our approach can resolve the controversy over the conformation of the myosin VI lever arm (Spink et al., 2008; Mukherjea et al., 2009; Thirumalai and Zhang, 2010). Beyond the myosin superfamily, we anticipate our approach may also prove useful for understanding kinesin and dynein motors as well as other biomolecules that can be suitably modeled by fluctuating polymer chains.
Materials and methods
The theoretical methods used above are fully described in Appendix 1, while the details of the Brownian dynamics simulations are in Appendix 2.
Appendix 1
In this appendix we develop the analytical theory of myosin V. In the first two sections, we define the actin geometry and solve the first passage problem for binding to sites along the actin double helix. From this we derive step distributions, the backwardforward step ratio, the mean run length, and the mean run velocity. This analysis is quite general: it applies to any actinbased motor with a similar mechanochemical cycle, assuming one can compute the equilibrium probability density for the positions of the free head during diffusion. In later sections, we compute this probability density for myosin V, incorporating effects of an interleg constraint potential, volume exclusion, orientation requirements for binding, and load force. The combined result gives a complete diffusioninformed kinetic description of myosin V.
Actin geometry
The actin double helix is composed of two filaments each containing 13 actin monomers per helical rotation. As it walks, the myosin heads can bind to any monomer along the actin strand. As discussed in the main text, the geometry of the actin helix is given by
where $R$ = 5.5 nm is the radius of the helix, $\mathrm{\Delta}z=72/13\approx 5.5\phantom{\rule{thinmathspace}{0ex}}\mathrm{n}\mathrm{m}$ is the size of each actin monomer, and ${\varphi}_{n}=12\pi n/13$ is the angle between adjacent monomers (Lan and Sun, 2006). Even and odd $n$ respectively designate the positions of monomers on the 1^{st} and 2^{nd} filament of the double helix. In the following analysis, we assume the bound leg is attached to the $n=0$ at position (0, 0, 0) on the first filament. Our work is a generalization of that by Hinczewski et al. (2013), in which only the half helical binding sites ${\mathbf{\mathbf{r}}}_{\pm 13}$ ($z=\pm 36\phantom{\rule{thinmathspace}{0ex}}\mathrm{n}\mathrm{m}$) were considered.
First passage time analysis
We begin by applying first passage time analysis to derive the stepping and stomping probabilities at each binding site along the actin. The first passage time of the free myosin leg to a given binding site along the actin ${\mathbf{\mathbf{r}}}_{n}$ is approximately,
where $a=0.4\mathrm{nm}$ nm is the capture radius around the actin binding site, $D}_{h$ = 5.7 × 10^{–7} cm^{2}/s is the diffusion constant of the myosin V head, and $\mathcal{P}({\mathbf{\mathbf{r}}}_{n})$ is the equilibrium probability distribution for the free myosin head evaluated at the actin binding site ${\widehat{\mathbf{\mathbf{r}}}}_{n}$. For a detailed calculation of this result using the renewal approach (Van Kampen, 2007), see Hinczewski et al. (2013). The accuracy of this approximation relies on the assumption that $t}_{a}={a}^{2}/{D}_{h}\approx 0.18\phantom{\rule{thinmathspace}{0ex}}\mathrm{n}\mathrm{s$, the time scale describing the diffusion of the myosin head a distance $a$, is much less than the relaxation time of the myosin polymer legs ${t}_{r}\sim \mathcal{O}(1\text{}\mu \mathrm{s})$ which was calculated by Hinczewski et al. (2013) using Brownian molecular dynamics simulations (Ermak and McCammon, 1978).
The actin binding sites are located on the outer surface of each actin subunit. As discussed in the main text, we therefore require that the orientation of the free leg, defined in terms of its azimuthal orientation ${\varphi}_{f}$, must be similar to that of the target actin binding site. Thus, for the description of binding events, we use the joint density $\mathcal{P}({\mathbf{\mathbf{r}}}_{n},\delta {\varphi}_{1}>({\varphi}_{f}+\pi ){\varphi}_{n}>\delta {\varphi}_{2})$. The two angles $(\delta {\varphi}_{1},\delta {\varphi}_{2})$ define a (possibly asymmetric) acceptance region, that is the deviations from the actin orientation ${\varphi}_{n}=12\pi n/13$ for which the head can still bind. The mean first passage time of the free head to a given binding site with the correct orientation is therefore,
This expression is used below in our derivations of physical observables for myosin V, connecting these quantities to mechanical properties of the polymer myosin lever arm.
While in the waiting state, both myosin heads are bound to the actin with ADP. The myosin can initiate a diffusive search through two kinetic pathways: (I) ADP is released from the trailing head (TH) and replaced by ATP, releasing the head from the actin, (II) the leading head (LH) detaches from the actin without ADP release. In the second scenario, dissociation of the LH occurs at a rate ${t}_{d2}^{1}$ much less than ${t}_{d1}^{1}$, the rate of ADP release/ATP binding in the first scenario. The probability of (I) and (II) occurring are then $g{(1+g)}^{1}$ and ${(1+g)}^{1}$ respectively, in terms of the gating parameter $g={t}_{d2}/{t}_{d1}\gg 1$.
After the TH is released, it must hydrolyze ATP to execute the recovery stroke, which occurs at the rate ${t}_{h}^{1}$. The free leg then diffusively searches for a binding site, reaching the capture radius of the n^{th} site with rate ${({t}_{\text{fp}}^{n})}^{1}$. Because the recovery stroke orients the myosin head to favor forward binding, we implement a binding penalty $b<1$ for sites behind the leading leg ($n<0$) designating the probability of binding once the site is reached. The binding penalty at site $n$ is denoted by,
With this definition, the effective diffusion/binding rate to each site is ${b}_{n}{({t}_{\text{fp}}^{n})}^{1}$. Given the first passage times to each of the binding sites, we can calculate the distribution of free leg binding times to the ${\mathbf{\mathbf{r}}}_{n}$ site,
where ${r}_{T}={\sum}_{n}{b}_{n}{({t}_{\text{fp}}^{n})}^{1}$. Here the sum is over all binding sites which the myosin head can reach, that is those sites with ${\mathbf{\mathbf{r}}}_{n}<2L$, where $L$ is the length of the myosin arms. The integral in Equation 16 is the convolution of the probability of hydrolysis occurring by some time ${t}^{\prime}$ and the probability of binding to the actin site over the interval time interval $t{t}^{\prime}$. These individual probabilities are modeled as exponential random variables with rate constants given by ${({t}_{h})}^{1}$ and ${b}_{n}{({t}_{\text{fp}}^{n})}^{1}$ respectively. The overall probability of binding to site $n$ given trailing head dissociation is:
where the exponential term implements the constraint that the free head must find a binding site before the bound head dissociates. From Equation 16 we can also see that the average time to bind ${t}_{T}$ is independent of target site:
The diffusive search of the LH is analogous with the following changes: (1) The LH retains its ADP molecule, so ATP hydrolysis is not required, (2) the LH is in the post power stroke orientation and favors binding to backward sites, behind the bound trailing leg. Thus, we obtain the LH binding probability and mean binding time by substituting ${t}_{h}\to 0$, ${b}_{n}\to {b}_{n}$, ${r}_{T}\to {r}_{L}={\sum}_{n}{b}_{n}{({t}_{\text{fp}}^{n})}^{1}$ in the above equations,
From the two distributions, ${\mathcal{P}}_{T}^{n}$ and ${\mathcal{P}}_{L}^{n}$, we can calculate the probabilities of each kinetic pathway. The probability of forward steps, trailing stomps, leading stomps and backward steps are respectively
where the factors $g{(1+g)}^{1}$ and ${(g+1)}^{1}$ account for the probability of trailing versus leading head detachement. The final pathway is termination, in which the free leg cannot complete its diffusive search before the bound leg detaches, leading to complete dissociation of the myosin from actin. The termination probability is simply, ${\mathcal{P}}_{t}=1{\mathcal{P}}_{f}{\mathcal{P}}_{Ts}{\mathcal{P}}_{Ls}{\mathcal{P}}_{b}$.
For comparisons to electron microscopy experiments (Oke et al., 2010), we require the distribution of head separation distances. The probability of observing bound myosin with $n$ actin subunits between the head is just the probability of taking a step/stomp of size $n$. This probability has contributions from each of the four primary kinetic pathways,
The full step size distribution, measured by florescence imaging experiments (Yildiz et al., 2003; Sakamoto et al., 2005), is then the convolution of the binding distribution ${\mathcal{P}}_{T}^{n}$ with ${\mathcal{P}}_{\text{dist}}^{n}$. This gives,
where ${z}_{n}=n\mathrm{\Delta}z/2$ and $\mathrm{\Delta}z$ is the size of the actin monomers. Similarly, before the release of the LH, the distances are distributed according to Equation 21, but now the distance to the bound TH is negative, so substituting ${\mathcal{P}}_{\text{dist}}^{m}\to {\mathcal{P}}_{\text{dist}}^{m}$ and summing over negative $m$, the full step size distribution for the leading head is given by,
For comparison to experiments, we only consider the trailing leg portion of step distribution, ${\mathcal{P}}_{T}({z}_{n})$, which we normalize appropriately. The combined leading and trailing leg step distribution, however, allows us to study the intricate changes in forward and backward steps/stomps in response to an applied load force. This behavior is highlighted in Figure 5 in the main text. The required distribution is
where again trailing leg and leading leg events are weighted by the gating ratio. This distribution will also allow us to accurately calculate run length and velocity, accounting for the full distribution of step sizes and possibility of forward or backward motion due to stomping.
Given the termination probability ${\mathcal{P}}_{t}$ of each diffusive search, the average number of steps taken is then ${\sum}_{n=1}^{\mathrm{\infty}}n{(1{\mathcal{P}}_{t})}^{n1}{\mathcal{P}}_{t}=1/{\mathcal{P}}_{t}$. Each step of size $n\mathrm{\Delta}z/2$ occurs with probability $\mathcal{P}({z}_{n})/(1{\mathcal{P}}_{t})$. The total run length is defined as the change in position of the center of mass of the myosin, which is half the change in position of the free leg. Summing over all step sizes, we obtain
The average run velocity is defined as ${v}_{\text{run}}={z}_{\text{run}}/{t}_{\text{run}}$, where ${t}_{\text{run}}$ is the average runtime. Before release of the TH or LH, the average waiting time bound to the actin is ${t}_{d1}{t}_{d2}/({t}_{d1}+{t}_{d2})=g{t}_{d1}/(1+g)$. After release, the mean binding times of the TH and LH are ${t}_{T}$ and ${t}_{L}$ respectively. The contributions of these times to the average run time is weighted by the sum of the TH/LH binding probabilities,
Equilibrium probability of the free myosin head
The mean first passage time, ${t}_{\text{fp}}^{n}$ used in the preceding analysis relied on knowledge of the free leg equilibrium probability distribution, $\mathcal{P}({\mathbf{\mathbf{r}}}_{n})$ evaluated at actin site $n$. In this section we derive the general equilibrium probability $\mathcal{P}(\mathbf{\mathbf{r}})$ for an arbitrary position $\mathbf{\mathbf{r}}$. Projection of this distribution onto, for instance, the $xz$ plane, gives a diffusion contour, like that measured by Andrecka et al. (2015). Data from this experiment indicate that there may be an interleg joint potential ${\mathscr{H}}_{J}$ constraining the diffusion of myosin V. The precise form of the potential is discussed later, but we assume it is only a function of the interleg angle, or equivalently ${\widehat{\mathbf{\mathbf{r}}}}_{f}\cdot {\widehat{\mathbf{\mathbf{r}}}}_{b}$, where ${\mathbf{\mathbf{r}}}_{f}$ and ${\mathbf{\mathbf{r}}}_{b}$ are the endtoend vectors for the free and bound legs. We also include a volume exclusion potential to capture the effects of steric repulsion between the myosin legs. Volume exclusion effects give rise to an effective repulsion between the bound and free heads given by the potential ${\mathscr{H}}_{V}={k}_{B}T{({d}_{V}/r)}^{6}$, where $r$ is the distance between the bound and free heads and ${d}_{V}$ is the effective length scale of the repulsion (see main text for more details). The probability can be obtained from the convolution of the endtoend probability of the bound leg, ${\mathcal{P}}_{b}({\mathbf{\mathbf{r}}}_{b})$ with that of the free leg ${\mathcal{P}}_{f}({\mathbf{\mathbf{r}}}_{f})$, while accounting for the joint potential and volume exclusion through the appropriate Boltzmann factors,
where the constraint $\mathbf{\mathbf{r}}={\mathbf{\mathbf{r}}}_{b}+{\mathbf{\mathbf{r}}}_{f}$ is enforced with the $\delta $function and $\beta =1/({k}_{B}T)$. Due to the volume exclusion term and because the joint constraint is a nonuniform weight in the convolution, we also require an overall normalization constant $A$. Each leg is modeled as a inextensible polymer with length L and persistence length ${l}_{p}$, and as noted above, the bound leg is assumed to be fixed at the origin $\mathbf{\mathbf{r}}=0$.
The free leg endtoend equilibrium probability may be accurately approximated using polymer mean field theory (Thirumalai and Ha, 1998),
where $\kappa =L/{l}_{p}$, ${\xi}_{f}=1{r}_{f}^{2}/{L}^{2}$, and ${A}_{f}$ is a normalization constant. Because the polymer is inextensible the dimensionless parameter, ${\xi}_{f}$ is constrained to be between 0 and 1. The endtoend vector ${\mathbf{\mathbf{r}}}_{f}$ and be fully specified using ${\xi}_{f}$ as well as polar and azimuthal angles, ${\theta}_{f}$ and ${\varphi}_{f}$. The normalization of the distribution then requires
which gives the normalization constant
For stiff polymers ($\kappa \approx 0$) relevant to myosin V, the distribution approximates a delta function at ${r}_{f}=L$. In the flexible regime ($\kappa \gg 1$), which may accurately describe motors with flexible components, the distribution becomes a spherically symmetric Gaussian centered at $\mathbf{\mathbf{r}}=0$.
The bound leg is assumed to be constrained in some direction ${\widehat{\mathbf{\mathbf{u}}}}_{c}$ by the power stroke. The tangent vector of this leg at the origin ${\widehat{\mathbf{\mathbf{u}}}}_{0}$ is therefore subject to a harmonic constraint with energy ${\mathscr{H}}_{c}=\frac{1}{2}{k}_{B}T{\nu}_{c}{({\widehat{\mathbf{\mathbf{u}}}}_{0}{\widehat{\mathbf{\mathbf{u}}}}_{c})}^{2}$, where there parameter ${\nu}_{c}$ is the strength of the constraint. This constraint biases the probability in the direction of ${\widehat{\mathbf{\mathbf{u}}}}_{c}$. We assume this effect is captured by the following ansatz:
where ${\xi}_{b}=1{r}_{b}^{2}/{L}^{2}$, ${A}_{b}$ is a normalization constant, and $\mathcal{T}$ is a function of ${\nu}_{c}$ satisfying $\mathcal{T}=0$ when ${\nu}_{c}=0$ so that in the limit of no constraint, the distribution reduces to that of a free leg. With this addition, the normalization constant becomes
The function $\mathcal{T}$ was determined in Hinczewski et al. (2013) by matching the moments of this distribution to the exact moments for a semiflexible polymer under harmonic constraint. In the limit that $\kappa \ll 1$ and ${\nu}_{c}\gg 1$ (large stiffness and strong constraint), the result is
With the two probability distributions ${\mathcal{P}}_{f}({\mathbf{\mathbf{r}}}_{f})$ and ${\mathcal{P}}_{b}({\mathbf{\mathbf{r}}}_{b})$, we can now evaluate the convolution in Equation 27,
In the second line we integrate over the bound leg vector ${\mathbf{\mathbf{r}}}_{b}$ with the delta function constraining ${\mathbf{\mathbf{r}}}_{b}=\mathbf{\mathbf{r}}{\mathbf{\mathbf{r}}}_{f}$ so that
where ${\theta}_{f}$ is the angle between the free leg end to end vector ${\mathbf{\mathbf{r}}}_{f}$ and the target position $\mathbf{\mathbf{r}}$.
The bound leg is constrained by the harmonic potential to the constraint direction, ${\widehat{\mathbf{\mathbf{u}}}}_{c}$, which has angular representation $({\theta}_{c},{\varphi}_{c}=0)$ with respect to the zaxis. To proceed, however, we consider the axis along the direction $\widehat{\mathbf{\mathbf{r}}}$. With respect to this axis, we can represent the unit vectors ${\widehat{\mathbf{\mathbf{r}}}}_{f}$ and ${\widehat{\mathbf{\mathbf{u}}}}_{c}$ by their polar and azimuthal angles $({\theta}_{f},{\varphi}_{f})$ and $({\stackrel{~}{\theta}}_{c},{\stackrel{~}{\varphi}}_{c})$. Then we may write the dot products as,
and
The integration element is then $d{\mathbf{\mathbf{r}}}_{f}={r}_{f}^{2}d{r}_{f}d\mathrm{cos}{\theta}_{f}d{\varphi}_{f}$ and after integrating over ${\varphi}_{f}$, Equation 34 becomes
where ${\mathcal{T}}_{z}=\mathcal{T}\mathrm{cos}{\stackrel{~}{\theta}}_{c}$ and ${\mathcal{T}}_{x}=\mathcal{T}\mathrm{sin}{\stackrel{~}{\theta}}_{c}$. Performing the change of variables,
we see that the main contributions to the integral in the stiff limit $\kappa \ll 1$ occur when ${\xi}_{f}\ll 1$ and ${\xi}_{b}\ll 1$. In this limit, the change of variables becomes ${r}_{f}\to L$, and $\mathrm{cos}{\theta}_{f}\to r/(2L)$ with Jacobian determinant $\mathrm{d}\mathrm{e}\mathrm{t}\phantom{\rule{thinmathspace}{0ex}}J={L}^{2}/(4r)$ and upper integration limits $u=2r/L{r}^{2}/{L}^{2}$. Then Equation 38 becomes
After transforming back to the zaxis frame of reference, using $\mathrm{cos}{\stackrel{~}{\theta}}_{c}=\widehat{\mathbf{\mathbf{r}}}\cdot {\widehat{\mathbf{\mathbf{u}}}}_{c}$, the integral Equation 40 evaluates to,
This expression gives an analytical description of the free head diffusion as the myosin takes a step. Projecting the distribution onto a given plane produces a 2D diffusion contour, like that measured by Andrecka et al. (2015). Depending on the joint constraint Hamiltonian and volume exclusion term, the normalization $A$ must generally be determined numerically. If ${\mathscr{H}}_{J}=0$ and ${d}_{V}=0$, then $A=1$. In this case and for $n=\pm 13$, we have ${\widehat{\mathbf{\mathbf{r}}}}_{\pm 13}\cdot {\widehat{\mathbf{\mathbf{u}}}}_{c}=\pm \mathrm{cos}{\theta}_{c}$ and $u\approx 1$, so that the result of Equation 41 reduces to that previously obtained for the half helical binding sites (Hinczewski et al., 2013).
Joint constraint Hamiltonian
In light of experiments by Andrecka et al. (2015), we use a joint potential ${\mathscr{H}}_{J}$ with an energy minimum at a certain preferred angle ${\theta}_{p}$. Instead of a freely rotating joint, this potential has an energy cost for deviations from ${\theta}_{p}$. Specifically we use the potential,
where ${\mu}_{c}$ is the dimensionless constraint strength and ${\theta}_{J}$ is the interleg angle at the joint. As discussed in the main text, this choice reduces to a harmonic potential when ${\theta}_{J}\approx {\theta}_{p}$. In large persistence length limit, $\kappa \ll 1$, we have $\mathrm{cos}{\theta}_{J}\approx {\widehat{\mathbf{\mathbf{r}}}}_{f}\cdot {\widehat{\mathbf{\mathbf{r}}}}_{b}$. Above we found that in this limit, the potential becomes ${\mathscr{H}}_{J}[{\widehat{\mathbf{\mathbf{r}}}}_{f}\cdot {\widehat{\mathbf{\mathbf{r}}}}_{b}]\approx {\mathscr{H}}_{J}[{r}^{2}/(2{L}^{2})1]$. Expanding the cosine in our potential we arrive at,
Orientation constraint for binding
The above model gives a description of the myosin diffusion: positions in the three dimensional space are visited by the free myosin head with probability $\mathcal{P}(r\widehat{\mathbf{\mathbf{r}}})$. To bind to actin, the free head must come sufficiently close to a binding size ${\mathbf{\mathbf{r}}}_{n}$ and have the correct orientation, pointing toward the outer surface of actin subunit.
As discussed above, the density relevant for the first passage time to reach a binding site with the correct orientation is ${\mathcal{\wp}}_{\text{bind}}({\mathbf{r}}_{n})=\mathcal{\wp}({\mathbf{r}}_{n},\delta {\varphi}_{1}>({\varphi}_{f}+\pi ){\varphi}_{n}>\delta {\varphi}_{2})$. To evaluate this quantity we repeat the derivation from the previous section, but integrate ${\varphi}_{f}$ only over the acceptance region $({\varphi}_{n}\delta {\varphi}_{2},{\varphi}_{n}+\delta {\varphi}_{1})$. Note we can choose our axes so that ${\stackrel{~}{\varphi}}_{c}=0$. In principle we should use ${\stackrel{~}{\varphi}}_{n}$, the azimuthal orientation of actin subunit $n$ with respect to the $\hat{\mathbf{r}}}_{n$ axis. The sites with nonnegligible probability of binding, however, have $z\gg x,y$, so that ${\stackrel{~}{\varphi}}_{n}\approx {\varphi}_{n}$. Carrying out the derivation, we arrive at
To evaluate this probability the constant $A$ must be calculated numerically and the integral over ${\varphi}_{f}$ is computed numerically for each target binding site. This quantity gives us the first passage times in Equation 14 and allows for computation of step distributions and other physical observables derived above.
Equilibrium probability of the free myosin head under load
If we apply a load force $\mathbf{\mathbf{F}}$ exerted at the end of the bound the myosin leg, the distribution Equation 31 is multiplied by a factor $\mathrm{exp}(\beta F{r}_{b}\widehat{\mathbf{\mathbf{F}}}\cdot {\widehat{\mathbf{\mathbf{r}}}}_{b})$. This term is simply the Boltzmann factor describing the energy cost required for the myosin leg to move a distance ${r}_{b}$ in the direction ${\widehat{\mathbf{\mathbf{r}}}}_{b}$ under the load force. Again assuming the stiff limit, (${r}_{b}\approx L$), this expression becomes $\mathrm{exp}(\beta FL\widehat{\mathbf{\mathbf{F}}}\cdot {\widehat{\mathbf{\mathbf{r}}}}_{b})$. Thus, the load force simply changes the effective constraint direction and strength, ${\widehat{\mathbf{\mathbf{u}}}}_{c}\to {\widehat{\mathbf{\mathbf{u}}}}_{c}^{\prime}$ and $\mathcal{T}\to {\mathcal{T}}^{\prime}$ where
Making these substitutions in Equation 44, gives the equilibrium freeleg probability of being at a binding site with the orientation required to bind for myosin V under load force $\mathbf{\mathbf{F}}$.
Appendix 2
Brownian dynamics simulations
In order to explore the dynamics of a myosin motor and explicitly incorporate excluded volume and other constraints, we used Brownian Dynamics (BD) simulations. In the simulations, myosin is modeled as a chain of beads of diameter d_{M} = 2 nm. The wildtype myosin V dimer has two lever arms that include 6 IQ domains each. We model the dimer as a polymer of N = 35 beads. To study how the lever arm length affects the step size distributions, we created myosin V models with 4 IQ domains and 8 IQ domains. The models consisted of 25 and 45 beads, respectively.
The Hamiltonian that governs the interactions between the constituent beads of myosin V and myosinactin interactions in the BD simulations is given as:
The five contributions to the Hamiltonian are listed below. First, the bond potential between myosin V beads is given as:
where ${r}_{i,i+1}$ is the centertocenter distance between the beads $i$ and $i+1$. The bond strength in our simulations is k_{s} = 700 kcal/ mol/ nm^{2}. In practice, ${k}_{s}$ must be large enough to prevent the contour length of the polymer from changing significantly. Next,
where ${r}_{i,j}$ is the distance between two beads. This term ensures that there are no unphysical overlaps between any beads. In the equation above, ${d}_{M}$ = 2 nm which corresponds to the diameter of a bead. The potential strength $\u03f5$ = 1 kcal/mol.
A similar term, ${H}_{AM,repel}$, is included to ensure that a myosin lever arm is sterically repulsed by the actin filament. For ${H}_{AM,repel}$, the potential has the a similar form as Equation 48 with the parameter, $d}_{AM$ = 5.5 nm (the sum of the radii of a bead and the actin filament) and the analog to ${r}_{i,j}$ is calculated as the perpendicular distance from the center of a lever arm bead to the actin filament.
Actinmyosin binding is given by a screened Coulomb potential,
where the potential strength, $V}_{0$ = 3.3 kcal/mol, the screening length $1/\alpha$ = 0.5 nm, and $\rho $ is the distance between the binding region on the stepping head of myosin V and any actin binding site. The actin binding region of the myosin motor is located on the surface of the last bead. We consider myosin bound and a trajectory complete when the binding energy becomes larger than 2.0 kcal/mol. We ran 200 stepping trajectories for each different setup we studied. Each trajectory ended when the stepping head was bound at an actin binding site.
The bending potential is given by:
where the bending constant, ${k}_{bend,i}$ is bond specific. Incorporating a heterogeneous bending constant allows us to create a flexible joint in the middle of the polymer and introduce varying amount of flexibility throughout the polymer structure. The bending constant is related to the persistence length ${L}_{p}$ of a polymer as: ${k}_{bend}={L}_{p}{k}_{B}T/{d}_{M}^{3}$, where ${k}_{B}T$ determines the thermal energy scale. The values of ${k}_{bend}$ used in the simulations are given in the table below.
Bead sequence number  

6 IQ model  4 IQ model  8 IQ model  ${k}_{bend}$ in kcal/mol 
1  1  1  150 
2  2  2  34.9 
3–17  3–12  3–22  23.7 
18  13  23  0 
19–33  14–23  24–43  23.7 
34  24  44  34.9 
35  25  45  N/A 
Power stroke
The first term in the bending potential, ${k}_{bend,1}({d}_{M}^{2}{\overrightarrow{r}}_{1,2}\cdot {\overrightarrow{r}}_{0,1})$, is used to create a power stroke by including a vector, ${\overrightarrow{r}}_{1,0}$ that makes an angle ${\theta}_{c}={65}^{\circ}$ with the actin filament. Specifically, the first term in the bending potential is zero when the angle between the vector created by the first two beads and the imposed direction is zero; the energy penalty increases as the attached lever arm direction deviates from the direction imposed by constraint angle, ${\theta}_{c}$. In practice, a large ${k}_{bend,1}$ keeps the bound head in the postpowerstroke orientation.
External force
To analyze the effect of external pulling force, we add a constant force in the $\hat{\mathbf{z}}$ direction to the myosin lever arm junction (i. e. the middle bead). The other forces and parameters remain unchanged.
Joint constraint
To study a model for myosin V that steps in a compasslike motion, we introduce an additional harmonic constraint between the two beads adjacent to the junction. We report data where the spring constant of the constraint is 5 k_{B}T/nm^{2} and the equilibrium length of the spring corresponds to a 90° angle between the two lever arms.
Glass cover
To investigate the effect of a cover slip (or any other excluded space) on the dynamics of myosin V dimer, we excluded half of the search space using a softcore repulsive potential similar to Equation 48): ${H}_{glass}=\u03f5{\sum}_{i=1}^{N}{\left(\frac{{d}_{C}}{{r}_{i}}\right)}^{6}$. The data reported in this paper was obtained using $d}_{C$ = 1 nm, $\u03f5$ = 5 kcal/mol, and ${r}_{i}$ is measured as a perpendicular distance from a bead to the glass plate that is directly under the actin filament in the $xz$plane.
References
 1
 2
 3
 4

5
Chemomechanical coupling and motor cycles of myosin VBiophysical Journal 100:1747–1755.https://doi.org/10.1016/j.bpj.2011.02.012
 6

7
Forcedependent stepping kinetics of myosinVBiophysical Journal 88:4402–4410.https://doi.org/10.1529/biophysj.104.053504
 8
 9

10
Spatial fluctuations affect the dynamics of motor proteinsThe Journal of Physical Chemistry B 112:11112–11121.https://doi.org/10.1021/jp800982b
 11

12
Dynamics of the unbound head during myosin V processive translocationNature Structural & Molecular Biology 14:246–248.https://doi.org/10.1038/nsmb1206
 13

14
Brownian dynamics with hydrodynamic interactionsThe Journal of Chemical Physics 69:1352–1360.https://doi.org/10.1063/1.436761
 15
 16
 17
 18

19
Load and pi control flux through the branched kinetic cycle of myosin VThe Journal of Biological Chemistry 283:17477–17484.https://doi.org/10.1074/jbc.M800539200
 20
 21

22
The azimuthal path of myosin V and its dependence on leverarm lengthThe Journal of General Physiology 139:101–120.https://doi.org/10.1085/jgp.201110715
 23
 24
 25

26
Robust processivity of myosin V under offaxis loadsNature Chemical Biology 6:300–305.https://doi.org/10.1038/nchembio.322
 27
 28

29
Velocity, Processivity, and individual steps of single myosin V molecules in live cellsBiophysical Journal 96:4268–4275.https://doi.org/10.1016/j.bpj.2009.02.045

30
Molecular mechanisms controlling actin filament dynamics in nonmuscle cellsAnnual Review of Biophysics and Biomolecular Structure 29:545–576.https://doi.org/10.1146/annurev.biophys.29.1.545
 31

32
Class V myosinsBiochimica Et Biophysica Acta (BBA)  Molecular Cell Research 1496:36–51.https://doi.org/10.1016/S01674889(00)000070

33
A model of myosin V processivityJournal of Biological Chemistry 279:40100–40111.https://doi.org/10.1074/jbc.M402583200

34
Direct observation of processive movement by individual myosin V moleculesBiochemical and Biophysical Research Communications 272:586–590.https://doi.org/10.1006/bbrc.2000.2819

35
Stepsize is determined by neck length in myosin VBiochemistry 44:16203–16210.https://doi.org/10.1021/bi0512086
 36
 37

38
Direct observation of the myosinVa power stroke and its reversalNature Structural & Molecular Biology 17:590–595.https://doi.org/10.1038/nsmb.1820
 39
 40

41
Kinetics of diffusion‐controlled reaction between chemically asymmetric molecules. I. general theoryThe Journal of Chemical Physics 54:2981–2988.https://doi.org/10.1063/1.1675283

42
Long single αhelical tail domains bridge the gap between structure and function of myosin VINature Structural & Molecular Biology 15:591–597.https://doi.org/10.1038/nsmb.1429
 43

44
Leverarm mechanics of processive myosinsBiophysical Journal 101:1–11.https://doi.org/10.1016/j.bpj.2011.05.026
 45
 46

47
Myosin VI: how do charged tails exert control?Structure 18:1393–1394.https://doi.org/10.1016/j.str.2010.10.004

48
Mechanochemical coupling of two substeps in a single myosin V motorNature Structural & Molecular Biology 11:877–883.https://doi.org/10.1038/nsmb806
 49

50
The gated gait of the processive molecular motor, myosin VNature Cell Biology 4:59–65.https://doi.org/10.1038/ncb732

51
Loaddependent kinetics of myosinV can explain its high processivityNature Cell Biology 7:861–869.https://doi.org/10.1038/ncb1287

52
Elastic leverarm model for myosin VBiophysical Journal 88:3792–3805.https://doi.org/10.1529/biophysj.104.046763

53
Influence of fluctuations in actin structure on myosin V step sizeJournal of Chemical Information and Modeling 45:1672–1675.https://doi.org/10.1021/ci050182m
 54
 55
Decision letter

Aleksandra M WalczakSenior Editor; École Normale Supérieure, France

Taekjip HaReviewing Editor; Johns Hopkins University, United States

Anatoly KolomeiskyReviewer; Rice University, United States

James R SellersReviewer
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
This manuscript revisits the modeling of the stepping mechanics of the processive actindependent motor, myosin V, by introducing a structural constraint wherein the angle between the lever arms connecting the two heads is relatively fixed as suggested by a more recent experimental paper (Andrecka et al., 2015). That paper suggested that, rather than function as a freely jointed swivel the junction between the two lever arms was relatively fixed and that the free head of myosin V did not explore random space while searching for a new, forward actin binding site. Rather, it occupied an offactin dwell space from which it periodically explored the actin landscape. The new model introduces this structural constraint and examines the consequences in terms of free head localization, selection of forward actin monomer binding site, step size and step directionality. In addition, the authors examine the effect of modeled load on several parameters, including the effect of off axis load. The new model predicts the experimental free head localization data obtained in Andrecka et al., 2015, as well as the kinetic and mechanical data found in many other myosin V papers.
Decision letter after peer review:
Thank you for submitting your article "Myosin V executes steps of variable length via structurally constrained diffusion" for consideration by eLife. Your article has been reviewed by two peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Aleksandra Walczak as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Anatoly Kolomeisky (Reviewer #1); James R. Sellers (Reviewer #2).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
This manuscript revisits the modeling of the stepping mechanics of the processive actindependent motor, myosin V, by introducing a structural constraint wherein the angle between the lever arms connecting the two heads is relatively fixed as suggested by a more recent experimental paper (Andrecka et al., 2015). That paper suggested that, rather than function as a freely jointed swivel the junction between the two lever arms was relatively fixed and that the free head of myosin V did not explore random space while searching for a new, forward actin binding site. Rather, it occupied an offactin dwell space from which it periodically explored the actin landscape. The new model introduces this structural constraint and examines the consequences in terms of free head localization, selection of forward actin monomer binding site, step size and step directionality. In addition, they examine the effect of modeled load on several parameters, including the effect of off axis load. The conclusion is that their new model predicts the experimental free head localization data obtained in Andrecka et al., 2015, as well as the kinetic and mechanical data found in many other myosin V papers. However, they add the comment that a freely diffusional model also fits most of the data in the literature with the exception of the discrete location for the dwell state.
This is a comprehensive theoretical work on uncovering mechanism of the motion of myosinV motor proteins. The subject is very important, and this work makes a huge step in our understanding of this difficult problem. The paper is very nicely written and all the details are clearly explained. The authors are noted for their selfcritical approach and for critically analyzing many experimental and theoretical results. For example, the discussion of asymmetry in diffusion contours is very honest and very objective. They combine computer simulations and analytical theory, the latter of which providing a better molecular understanding. This work should stimulate a lot of new experiments.
Essential revisions:
1) The use of mean firstpassage times as estimates for the binding reaction processes assumes that the reaction is taking place after a single collision (diffusion limit). While this might be the case, it is not proven in the case of motor proteins. I would add a caution at this point about the application of mean firstpassage times. By the way, this might also be the reason why the theory does not work at large forces. Some discussion might be useful here. For example, authors should note that Dunn et al., 2007, Andrecka et al., 2015, and Veigel et al., 2002, all estimated the actual transit time from detachment of the trail head to its reattachment at a new forward position to be much greater than the theoretic first passage time.
2) Figure 5. The Veigel et al., 2002 paper showed evidence of stiffness changes at stall forces that may represent the proposal of trail head stomping. While the duration of these events were not quantified in this paper, lifetimes of 200—500 ms are shown in Figure 5 of that paper at 100 μm ATP. Could the authors comment on whether their model would predict such long lifetimes of singleheaded attachment?
https://doi.org/10.7554/eLife.51569.sa1Author response
Essential revisions:
1) The use of mean firstpassage times as estimates for the binding reaction processes assumes that the reaction is taking place after a single collision (diffusion limit). While this might be the case, it is not proven in the case of motor proteins. I would add a caution at this point about the application of mean firstpassage times. By the way, this might also be the reason why the theory does not work at large forces. Some discussion might be useful here. For example, authors should note that Dunn et al., 2007, Andrecka et al., 2015, and Veigel et al., 2002, all estimated the actual transit time from detachment of the trail head to its reattachment at a new forward position to be much greater than the theoretic first passage time.
The reviewers raise an important point that was not sufficiently clear in the original manuscript. Binding in our model is not simply diffusionlimited, but this was not discussed explicitly in the main text, and we did not show any binding time results. For that reason we have added a new section to the main text, “Probing the biological function of the joint constraint: effects on timing and consistency of stepping”. Part of this section is devoted to a detailed description of binding timescales, along with a new figure (Figure 6) that shows a comparison of diffusion, hydrolysis, and binding timescales as a function of load force. Quoting from that new section:
“As discussed above, binding in the model is more complex than just waiting for the head to diffuse within the capture radius a of a binding site on actin. […] Figure 6 shows t_{T} versus t_{diff} as a function of applied force F, and we always see t_{T} > t_{diff}, as expected from the above constraints.”
The discussion of the experimental estimates for transit times is later in the same section, and we quote it in response to the next reviewer comment.
2) Figure 5. The Veigel et al., 2002 paper showed evidence of stiffness changes at stall forces that may represent the proposal of trail head stomping. While the duration of these events were not quantified in this paper, lifetimes of 200—500 ms are shown in Figure 5 of that paper at 100 μm ATP. Could the authors comment on whether their model would predict such long lifetimes of singleheaded attachment?
This is indeed an interesting experimental observation which we believe is linked to trailing leg stomping. We now discuss this experiment in our new section on timing that we have added to the main text. Quoting from that section:
“For small F the trailing leg kinetics is dominated by forward steps. t_{diff} < t_{h} in this regime, but even though the head can diffuse rapidly to forward binding sites, it does not bind until hydrolysis occurs. […] Here there is another experimental artifact in play: as discussed in Hinczewski et al.,2013, the drag from the gold nanoparticle can substantially increase diffusion times, with t_{T} becoming much larger than t_{h} at F = 0 because of slow diffusion.”
https://doi.org/10.7554/eLife.51569.sa2Article and author information
Author details
Funding
National Science Foundation (DGE1650441)
 David Hathcock
National Science Foundation (CHE 1900093)
 Dave Thirumalai
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
DH was supported by the Presidential Life Science Fellowship from Cornell and an NSF Graduate Research Fellowship (DGE1650441). DT acknowledges support from the National Science Foundation (CHE 19–00033) and the CollieWelch Chair (F0019).
Senior Editor
 Aleksandra M Walczak, École Normale Supérieure, France
Reviewing Editor
 Taekjip Ha, Johns Hopkins University, United States
Reviewers
 Anatoly Kolomeisky, Rice University, United States
 James R Sellers
Publication history
 Received: September 3, 2019
 Accepted: January 14, 2020
 Accepted Manuscript published: January 15, 2020 (version 1)
 Version of Record published: March 3, 2020 (version 2)
Copyright
© 2020, Hathcock et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 859
 Page views

 133
 Downloads

 2
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Download citations (links to download the citations from this article in formats compatible with various reference manager tools)
Open citations (links to open the citations from this article in various online reference manager services)
Further reading

 Physics of Living Systems
Advances in DNA sequencing have revolutionized our ability to read genomes. However, even in the most wellstudied of organisms, the bacterium Escherichia coli, for ≈ 65% of promoters we remain ignorant of their regulation. Until we crack this regulatory Rosetta Stone, efforts to read and write genomes will remain haphazard. We introduce a new method, RegSeq, that links massivelyparallel reporter assays with mass spectrometry to produce a base pair resolution dissection of more than 100 E. coli promoters in 12 growth conditions. We demonstrate that the method recapitulates known regulatory information. Then, we examine regulatory architectures for more than 80 promoters which previously had no known regulatory information. In many cases, we also identify which transcription factors mediate their regulation. This method clears a path for highly multiplexed investigations of the regulatory genome of model organisms, with the potential of moving to an array of microbes of ecological and medical relevance.

 Developmental Biology
 Physics of Living Systems
The coordination of cell proliferation and migration in growing tissues is crucial in development and regeneration but remains poorly understood. Here, we find that, while expanding with an edge speed independent of initial conditions, millimeterscale epithelial monolayers exhibit internal patterns of proliferation and migration that depend not on the current but on the initial tissue size, indicating memory effects. Specifically, the core of large tissues becomes very dense, almost quiescent, and ceases cellcycle progression. In contrast, initiallysmaller tissues develop a local minimum of cell density and a tissuespanning vortex. To explain vortex formation, we propose an active polar fluid model with a feedback between cell polarization and tissue flow. Taken together, our findings suggest that expanding epithelia decouple their internal and edge regions, which enables robust expansion dynamics despite the presence of size and historydependent patterns in the tissue interior.