Interrogating Emergent Transport Properties for Molecular Motor Ensembles: A Semi-analytical Approach

Intracellular transport is an essential function in eucaryotic cells, facilitated by motor proteins—proteins converting chemical energy into kinetic energy. It is understood that motor proteins work in teams enabling unidirectional and bidirectional transport of intracellular cargo over long distances. Disruptions of the underlying transport mechanisms, often caused by mutations that alter single motor characteristics, are known to cause neurodegenerative diseases. For example, phosphorylation of kinesin motor domain at the serine residue is implicated in Huntington’s disease, with a recent study of phosphorylated and phosphomimetic serine residues indicating lowered single motor stalling forces. In this article we report the effects of mutations of this nature on transport properties of cargo carried by multiple wild-type and mutant motors. Results indicate that mutants with altered stall forces might determine the average velocity and run-length even when they are outnumbered by wild type motors in the ensemble. It is shown that mutants gain a competitive advantage and lead to an increase in the expected run-length when the load on the cargo is in the vicinity of the mutant’s stalling force or a multiple of its stalling force. A separate contribution of this article is the development of a semi-analytic method to analyze transport of cargo by multiple motors of multiple types. The technique determines transition rates between various relative configurations of motors carrying the cargo using the transition rates between various absolute configurations. This enables a computation of biologically relevant quantities like average velocity and run-length without resorting to Monte Carlo simulations. It can also be used to introduce alterations of various single motor parameters to model a mutation and to deduce effects of such alterations on the transport of a common cargo by multiple motors. Our method is easily implementable and we provide a software package for general use.


Introduction
Motor proteins-kinesin, dynein and myosins-are nanoscale machines that are the main effectors of intracellular transport.They play a critical role in the growth and sustenance of healthy cells by enabling a transport of intracellular cargo over networks of microtubules [1].Disruption of the functions performed by these molecular motors is linked to neurodegenerative diseases such as Huntington's, Parkinson's and Alzheimer's Disease [2,3], muscular disorders such as heart disease, uterine complications and high blood pressure.The mechano-chemical behavior of a single motor moving along the microtubule substrate to transport a cargo is relatively well understood [4].There is considerable prior work on understanding how a singlemotor transports a cargo on a microtubule.In [5] the Fokker-Planck equation for a single particle moving in a one-dimensional potential is obtained and used to determine the probability distribution of the position of the particle in a single repeat length of the polymer track via discretization.The underlying approach was extended to a motor connected to a cargo by a linear spring [6,7] while travelling under the influence of a tilted periodic potential.In vivo it is known that multiple motors work in teams to transport a common cargo [16]; how multiple motor-proteins coordinate to transport a common cargo is not well understood [16][17][18].Many studies employ a probabilistic description of the behavior of a single motor protein to construct models that describe how multiple motors transport a common cargo [16,[19][20][21].In [22], Gross and coworkers employed Monte-Carlo simulation studies built on a model of a single kinesin in to explore how multiple identical kinesin motors might interact to transport a cargo against a hindering load force.Their work indicates, counter intuitively, the existence of a form of strain-gating, where the motors of an ensemble share loads unequally enabling cargo transport over longer distances.Xu and coworkers examined the effects of ATP concentration on the transport of cargo carried by single and two motors in [23].At decreased levels of ATP concentration, the velocities of cargoes transported by single and two motor proteins decreases.
Coincidentally at decreased ATP concentration there was an an appreciable increase in the run-length of cargoes transported by two motors, while no such effect was seen in the case of transport by one motor.The authors proposed that the increased run-length observed in the presence of two motors results from the lowered dissociation of each motor from the microtubule at decreased ATP concentrations and the increased probability that the cargo stayed bound to the microtubule.Studies such as [21,29] using probabilistic models of single motors have also predicted in that an ensemble of kinesin motors is a robust system and the robustness increases under high loads [37].
The study of cargo transport by a heterogeneous ensemble of motor proteins composed of both wild type and mutant motors is important to inform our understanding of how mutant motors impact intracellular transport and lead to an onset of diseases.Recent studies have implied that alterations in the kinesin-1 motor domain may have a role in impaired axonal transport.Phosphorylation of a mammal kinesin motor domain by kinase c-Jun N-terminal kinase-3 (JNK3) at a conserved serine residue (Ser-176 in A and C isoforms and Ser-175 in B isoform) is implicated in Huntington's disease [24].However the mechanisms affected by Ser-175 phosphorylation are not well understood.An experimental study by Selvin and coworkers in [25] reported that a negative charge at Ser-175, acquired through mutation or phosphorylation, leads to a lower stall force and decreased velocity under external loads of 1 pN or more, while leaving the ATPase, microtubule-binding affnity and processivity unchanged.
Using a semi-analytical method, we reveal surprising emergent transport behaviors arising when a cargo is transported by multiple motor proteins, some of which are mutated and some are not.In particular we analyze the impact of Ser-175 kinesin mutation such as those reported in [25] on cargo transport in the presence of wild type motors.The detailed investigation made possible by our method leads us to hypothesize that under certain conditions, a 'mooring mechanism' is activated where proteins moor the cargo to the microtubule and prevent it from being lost.While these mooring proteins do not contribute to the motion of the cargo, they enhance the probability of attachment of other cargo-bound motor proteins to the microtubule that subsequently contributes to an increase in average cargo displacement.The activation of mooring mechanism depends on a number of external factors such as load force and ATP concentration.However, it is also determined by intrinsic properties of the motor protein such as the stall force of the individual motors.Remarkably, mutant motors that have stalling forces matched to the external load force can act dominantly and determine emergent transport properties such as longer run-length, even when they are outnumbered in the ensemble by wildtype motors.Such mechanisms could potentially point to how diseased states emerge and progress coincident with the accumulation of the mutant motor species.
A separate contribution of the article is a semi-analytical method for determining the probability distribution of various configurations of a cargo carried by multiple number and types of motor-proteins.A detailed experimental study of the various modalities of transport by multiple motor proteins (homologous or otherwise) requires significant instrumental resolution than what is needed to investigate single motor behavior.As a consequence, observing the transport dynamics of multiple motors is experimentally challenging [26,27].It is further compounded by the combinatorial complexity introduced by the multitude of scenarios possible when many motors and motor types participate in transport.Such challenges motivate the use of analytical and computational tools.The mean-field approach in [19,21] makes use of simplifying assumptions, such as equal load sharing among all motors, to achieve analytical results, thereby sacrificing significant detail for computational benefits.The Monte Carlo approach in [22], provides better fidelity where complex models can be employed; however, the accuracy of results depend on the number of iterations and on the rarity of the events that occur.
Unlike the Monte-Carlo simulations or any currently implemented simulation methods that study multi-motor ensembles, our Master Equation based method [28] analytically solves for the probability distributions of all possible scenarios at any time point.The methodology is uniquely powerful and enables the calculation of various biologically relevant quantities such as average velocity and run-length, for reasonably sized ensembles and with high accuracy while using lesser computational resources and time as compared to Monte Carlo simulations (the methodology proposed in [6,7] also obtains results without resorting to Monte Carlo simulations with a single motor attached to the cargo).The underlying concepts behind our methodology are motivated by earlier work reported in [29].The key enabling concept is that of 'relative configuration' of motor proteins, determining the transition probabilities between relative configurations from the transition probabilities of the absolute configuration space and subsequently determining the biologically relevant quantities from the relative configuration space.The computation engine is implemented using MATLAB and can be used to simulate cargo transport by any two unidirectional species.
Our method provides a general platform to study the transport of cargo by multiple motors of two different types where each type of motor protein can be individually characterized by a probabilistic model describing its stepping, detachment and attachment rates.For this article, the technique has been utilized to introduce alterations of various parameters from the nominal ones to model a mutation of the serine residue and compute the effect of such a mutation on cargo carried by a mixture of wild-type and mutant motor-proteins.In summary, we developed a simulation engine to study the transport of cargo by multiple motor proteins with distinct properties that in concert can exhibit emergent transport behaviors.

Results and Discussion
In this study, we investigate the impact of a previously reported kinesin mutation on the transport of a motor ensemble and its attached cargo.The mutation, located within the motor domain of the kinesin, mimics the phosphorylated state of Ser-175 [25].The motor domain phosphomimetic mutation does not affect the ATPase, microtubule-binding affinity or processivity of the motor, but does reduce the stall force and velocity of kinesin under a load force.The in vitro phosphorylation of Ser-175 for a full length kinesin similarly reduces stall force and velocity of the motor.Both the mutant and in vitro phosphorylated kinesin showed no other aberrant single motor behavior.
Here we use computational modelling to analyze the impact of the Ser-175 mutant kinesin on a heterogenous ensemble of motors and its transported cargo.The wild-type and mutant kinesin motors have a different stalling force; the wild type kinesin has a stall force of F s = 6pN and the mutated kinesin has a reduced stall force of " F s ¼ 5:5pN(see [25] for more detail).In our analysis we considered cargoes transported by the following motor ensembles: cargo with two wild-type (WW) motors, one wild-type and one mutant (WM) motors, two mutant (MM) motors, three wild-type (WWW) motors, two wild-type and one mutant (WWM) motors, one wild-type and two mutant (WMM)motors, and three mutant (MMM) motors.

How External Loading Conditions Affect Behavior of Heterogenous Motor Ensemble
Motor configuration when engaged to the microtubule approaches steady state.For cargoes carried by multiple motor proteins, transport occurs as each motor steps along the microtubule lattice, sequentially binding, translocating and releasing from the microtubule in a well defined mechanochemical cycle [16,21,30,31].Eventually, the condition where no motors are attached to the microtubule will be reached and the cargo will diffuse away from the microtubule.Thus it is evident that the only steady state is the condition when none of the motors of the ensemble carrying the cargo remain attached to the microtubule.Indeed, in all the cases considered where there are three or less motors on the cargo (i.e.M 3) it is seen that for any specified initial probability distribution of motors engaged to the microtubule, the steady-state probability distribution has no motors attached to the microtubule i.e. cargo is eventually lost.Given that the steady-state probability distribution is the trivial case of no motors attached, nothing much of value can be derived regarding the motor behavior from the trivial steadystate probability distribution.However the conditional probability distribution of the various configurations of motor proteins carrying a common cargo, given that at least one motor is attached to the microtubule (i.e. when the motor-cargo complex is not completely detached from the microtubule), is of significant interest (see analogous description in [14] on discretization that separates the effects of attachment and cycling through the kinetic steps of a motor).Here it is apriori not evident whether the conditional probability distribution has a steady state.
Remarkably, the results show that under a constant load force on the cargo the conditional probability distribution of the various motor configurations reaches a steady state.The steady state distribution is independent of the initial probability distribution chosen (see for example, Fig 1(b) for the WM case).This indicates that the cargo-motor ensemble consisting of motors with different stalling forces is a robust system which, starting from any arbitrary initial condition, behaves according to a fixed distribution after some transient period.Furthermore, the existence of such a steady-state distribution, termed here as the steady-state conditional pdf, allows the computation of most of the variables of biological interest.For a homologous ensemble of molecular motors, recent results have shown that the associated Markov chain is ergodic [13,15], implying that the states eventually reach a unique steady state distribution.
Load forces on cargo affect conditional steady state distribution of motor ensemble.The conditional steady-state distribution is used to analyze transport of cargo with two attached motors and assumes that at least one motor is attached (Fig 2).Fig 2 shows the probabilities of the number of wild-type and mutated motor proteins that remain attached to the microtubule under varying load forces.The analysis indicates that when the cargo has two wild type motors attached (WW case), the probability that both motors remain attached to the microtubule peaks when load force is close to 12pN.This equals twice the stalling force of a single wild type motor; indicating a relationship between number of motors attached (obtained using the steady state distribution) and applied load.Thus at these loads there exists a propensity to retain and not loose the cargo.Similarly, in the MM case (where the cargo is attached to two mutant motors) the probability that both motors are attached to the microtubule is maximum near 11pN which is twice the stalling force of the mutated motor.The WM ensemble displays lower peaks for the dual attachment of motors at both load forces of 11pN and 12pN.
It is therefore evident that for ensembles containing two motors, the conditional steady state distributions of the number of motors attached is not significantly affected until the load force on the cargo is close to twice the stalling force values of the wild -type or mutant motor (2F s = 12pN or 2 " F s ¼ 11pN).A similar observation is made for ensembles with three motors and for load force values close to twice and thrice the stalling force of a single motor or 2F s , 2 " F s , 3F s (= 18pN) and 3 " F s ð¼ 16:5pNÞ.We hypothesize a 'mooring mechanism' to explain the observations where the probability of a motor remaining attached to the microtubule is increased when the load force on the motor is close to its stalling force, thereby mooring the motor to the microtubule.The moored motor prevents the motor-cargo complex from being detached from the microtubule.The mooring mechanism is evident in the case of a wild-type motor, where the probability of the wild-type motor remaining attached to the microtubule is high when the load force on the motor is close to F S = 6 pN Eq (16).The WW curve in Fig 2 thus indicates that when load force is close to 12pN, both the motors face a load force less than but close to F S = 6 pN allowing for the cargo to be 'moored' to the microtubule through both the motors.A similar reasoning can explain the MM curve in Fig 2 at load force close to 11pN, where the mooring mechanism for both the mutant motors is activated when the load force on each motor is close to " F S ¼ 5:5 pN Eq (17).It indicates that for load forces close to multiples of single motor stalling forces, equal sharing of loads is preferred by the motors in the ensemble.Thus a cooperative mode between an ensemble of motors is encouraged when the load forces are close to the stallforce or multiples of the stall forces.For example, if a WW ensemble is subjected to load force close to 2F s = 12 pN, unequal load sharing would subject one of the motors to a force greater than its stalling force F s = 6 pN.This will lead to its detachment and eventual transfer of the entire load onto the remaining motor.This, in turn, would gradually lead to the remaining single motor being detached from the microtubule, and the resultant loss of cargo attachment to the microtubule substrate.On the other hand, an equal load sharing will nearly stall both the motors but ensure the retention of the cargo.A similar inference can be made for cargo transport by three motors and load forces approaching the values F s , " F s , 2F s or 2 " F s .Note that, it is possible when the load force is close to twice the stalling force, the cargo detaches from the microtubule with a high probability.In such a case, the previously mentioned conclusions reached about the probabilities by applying the condition that the cargomotor complex does not detach from the microtubule (i.e. at least one motor remains engaged to the microtubule) will not carry much significance.We then utilized the probabilities obtained without the conditioning to reach similar conclusions.Figs 3 and 4 show that the time taken for the cargo to be lost (i.e.time taken for the probability of detachment of the cargo to be close to one) is higher when the load force is close to but less than the appropriate stalling force or twice the appropriate stalling force value.The time taken for detachment is lesser for a load force of an intermediate value (compare with Fig 5).The figures show that if the load force is matched to an appropriate stall force or its multiples (for example, close to 2F s for WW or 2 " F s for MM ensemble), the probability of detachment of the cargo from the microtubule becomes appreciable only after a long time has passed.Thus, the conclusions reached using the conditional distribution are also supported by results inferred from the entire configuration space with no restriction that the cargo remains attached.It implies that the study on the condition that 'the cargo-motor complex is not detached from the microtubule' provides important insights.
In conclusion, load force on the cargo has a significant effect on the steady state distribution of motor ensembles and can provide an insight into how multiple motor ensembles function under varying loads.In particular, when a cargo is subjected to load forces that are multiples of the stalling force for individual motors, a cooperative behaviour in the form of equal load sharing is observed.This averts a scenario where unequal load sharing would subject some of the motors to loads greater than their stall forces, leading to their detachment and the eventual loss of cargo from the microtubule.Such cooperative behavior possibly helps retain the cargo.It is thus more probable that cargoes attached to a higher number of wild-type motors are retained under load forces close to multiples of the stalling force of a wild-type motor, making the presence of more mutants disadvantageous in these load force regimes.

Hindering Load Force on the Cargo Can Tune Multiple Motor Travel
The knowledge of the steady state conditional probability distributions allows one to compute biologically relevant quantities such as average velocity and run-length.The results obtained for ensembles containing two and three motors for varying load forces are reported here.
Average velocity typically reduces with increasing number of mutant motors, but varies as load forces approach multiples of stall forces.The effect of load force on the average velocity of a cargo is described in Figs 6 and 7 for a total of two and three attached motors, respectively.It is evident that for any given ensemble, average cargo velocity mostly reduces with increase in the hindering load force on the cargo.It is further evident from Fig 6 that for a given load force on the cargo, the average velocities are highest for the WW ensemble, followed by WM and the lowest for MM ensemble.For the same load force and number of motors in the ensemble, it follows that the average velocity decreases with increased participation of mutant motors (which have a lower stall force).A similar conclusion is reached when three motors are attached to the cargo (see Fig 7).Thus, for a fixed load force and total number of motors, average cargo velocity reduces with increasing number of mutant motors in the ensemble.Inferences of this nature can be verified in-vitro by using laser-trap experiments to exert appropriate load forces on a cargo being carried by a quantified number of motor proteins and measuring the average cargo velocity for varying proportion of wild-type and mutant motors carrying the cargo.
However, for some motors ensembles the trend of decreasing cargo velocity at elevated load forces is not observed.12.1 pN.In the case of WWW ensemble, as the load forces increase to 12 pN (which is equal to 2F s ) the average velocity increases with increasing load force.Our semi-analytic method allows for a fine analysis of such counter-intuitive observations, where it is possible to extract the precise motor ensemble configurations that contribute to the observed effect.
Upon analysis it is observed that as the load force approaches the appropriate multiples of stalling force (2F s or 2 " F s in Fig 8), detachment of motors from the microtubule becomes less likely and the attached motors tend to cluster close together on the microtubule.This is similar to an observation made in [13,15,22,29] where under high loads, the motors carrying a common cargo tend to remain clustered together while under low loads they tend to spread apart.Since the detachments are less likely, backward transitions of the cargo induced by single motor detachment are less frequent.Moreover, if a detachment were to occur, the magnitude of load induced backward transition of the cargo is less due to the clustering of motors under high loads.On average, this contributes to an increase in average velocity despite increasing load.This explains the increase in average velocities of the three motor ensembles at load forces approaching 12 pN.
Interestingly, the increase in average velocities at load forces approaching 2F s = 12 pN is more prominent in the WWW ensemble than in MMM ensemble, indicating that it is the detachment of the wild-type motors in the ensemble that becomes less unlikely.Thus, as the number of wild type motors increases, detachment events at load forces close to 2F s = 12 pN become less likely.Similarly, at load forces approaching 2 " F s ¼ 11 pN the average cargo velocity increases with increasing load force; with the increase becoming more prominent for ensembles having more numbers of mutant motors (i.e. higher for MMM ensemble than WWW ensemble).
Significantly, this observation suggests that by tuning the load forces on a cargo to values near the multiples of the stall force, one could counter a decreased velocity resulting from an increase in load force or the presence of a higher number of mutant motors.
An important feature of the semi-analytic method to be emphasized is the ease with which detection of rare events such as those described above are made possible.Monte Carlo based methodology, for example, would not only take significant computations to simulate rare events with a high degree of confidence, but would also be unable to determine the incidence causes of rare events.Thus our model allows the user to easily unearth the cause of a rare event.
In conclusion, load forces approaching multiples of the appropriate stall force for a cargo decrease the probability of detachment of associated motors from the microtubule.Since detachment of motors in a multiple motor ensemble can lead to backward motion of the cargo (because kinesin doesn't actively step backward i.e. towards minus end), a decrease in the detachment probability leads to higher forward motion of cargo on average, leading to an increase in average velocity even though the load increases to approach multiples of stalling forces.If load forces are close to multiples of F s ðor " F s Þ detachment of wild-type (or mutant) motors in the ensemble decreases and cargoes with higher wild-type (or mutant) motors demonstrate greater increase in average velocity.
Cargo run-lengths impacted by heterogenity in the motor ensemble.The dependence of run-length on the load force when the cargo has two and three motors attached to it is shown in Figs 9 and 10 respectively.A general trend is that the average run-length reduces with an increase in the antagonistic load force on the cargo.Moreover, for a fixed number of motors in  (16), it is evident that the probability of the wild type motor remaining engaged with the microtubule lattice is high when the load force is close to its stalling force F S = 6 pN.It can thus be inferred that under this condition, the mooring mechanism is activated by the wild-type motor and the cargo remains bound to the microtubule with a higher probability.As a consequence, the other disengaged mutant motor has a greater opportunity to reattach to the microtubule, overcoming the stalled condition of the leading motor and contributing to the increase in run-length.Thus the run-length under these conditions will increase with an increase in the number of wild-type motors, since it is the wild-type motors that will tether the cargo to the microtubule lattice, favoring its continued translocation.The peak in the run-length observed in Fig 10 inset(C) for a load force close to 12 pN can similarly be explained.As the two leading wild-type motors stall and remain attached to the microtubule, the third motor will reattach, overcome the stalled condition and promote an extension of the run-length.
The same load dependency of cargo run-lengths is also observed for mutant motors with reduced stall forces.At load force equal to the reduced stall force of " F S , a mutant motor exhibits an increased probability for remaining attached to the microtubule lattice (from Eq (17)) and activates the mooring mechanism, preventing the cargo from diffusing away from the microtubule (Fig 11, see peak in green trace at 5.5pN).The tethering by the mutant motor increases the probability that the disengaged wild-type motor will reattach, overcome the stalled mutant motor and extend the cargo runlength.The same effect also explains the positive correlation between run-length and number of mutant motors at load forces close to " F s and multiples of " F s , as the probability of having mutant motors that are stalled increases, making the loss of cargo more unlikely.
Significantly, in these load force ranges and for a fixed number of motors attached to the cargo, the cargo with more mutant motors attached will travel further than other cargoes with fewer or no mutant motors attached.Remarkably, the run-length of the cargo is predicted by the single motor characteristics, namely the stalling force, of the mutant motor.One can test these inferences in-vitro by attaching a quantified number of wild-type and mutant motors to a cargo and using an optical tweezer to exert appropriate load forces to measure the effect of ensemble heterogenity on average distance traversed by the cargo.The relevance of studying cargo under various load forces to in-vivo functionality is detailed in [12].Note that, the sharpness of the transitions around stalling forces (F S and " F s ) and its multiples, as seen in Figs 9 and 10, is not being highlighted here.Instead, we observe that, the run-lengths of multi-motor ensembles tend to increase when the load forces are closer to stall forces or their multiples.
ATP concentration can tune multiple motor travel.It is known that kinesin hydrolyzes one ATP per step [32].Thus, the rate at which kinesin steps depends on its rate of ATP hydrolysis.It was shown experimentally in [23] that under no load conditions, average velocity for ensembles having one and two motors decreases with decreasing ATP concentration.It was also shown that run-length was not appreciably affected by changing ATP concentrations when the cargo had only one kinesin attached.However, when the cargo has two kinesin motors attached, run-length was demonstrated to have a strong negative correlation with single motor velocity (and single motor velocity was shown to have a positive correlation with ATP concentration).Our semi-analytic method was used to assess the conclusions reached in [23] and to further understand how the presence of mutant motors (with a different stalling force " F s ) can impact the ATP dependence of run length.The altering of this behavior due to the presence of more than two motors in the ensemble is also studied.
Average run-length for multiple motor ensemble negatively correlates to ATP concentration when restricted to certain regimes.We first considered the affect of ATP concentration on run-length of a cargo containing only a single wild type motor (Fig 12).Our analysis shows that as ATP concentration increases, the distance traveled by the cargo also increases monotonically until ATP concentration saturates.When the cargo carried by a single wild-type motor is subjected to 0.2 pN load force, the effect of ATP concentration on run-length is not substantial for concentrations above 10 μM (611 nm at 10 μM ATP and 786 nm at 2 mM ATP).Thus run-length of a cargo with a single wild-type motor is not substantially affected by ATP concentration, agreeing with the trend observed in [23].
However, the effect of ATP concentration is more pronounced for ensembles containing more than one motor.For ensembles comprised of multiple motors, decreased ATP concentrations under certain conditions leads to increased average run-length of the cargo.For example, in both the WW (Fig , indicating that the increased run-length is not due to increasing velocity but instead reflects the increased probability that the cargo remains attached to the microtubule for a longer time.As is evident from Fig 14 for WW ensemble, not only is the probability of cargo remaining attached via two motors more than via one motor, but the probability that the two motors remain attached to the microtubule also increases as the ATP concentration is decreased.A similar argument is provided for the experimental observations made by Xu and team in [23] for ensemble WW, that reported a negative correlation between single-motor velocities and run-length of cargo transport by ensemble WW.
Furthermore, the results also predict a negative correlation will persist until a certain minimum ATP concentration, below which the average run-length starts decreasing.To summarize, the runlength cannot continue to increase as ATP concentration falls, since at a certain threshold the low ATP concentrations will hinder ATP hydrolysis and the stepping of the motor domains along the microtubule lattice.At this limiting ATP concentration, regardless of how long the cargo remains attached, the motors will no longer be able to translocate along the microtubule lattice and runlength will decrease.
In conclusion, the ATP concentration in certain regimes negatively correlates with cargo run-length, as long as there are more than one motor attached to the cargo.Such behavior observed for a WW ensemble agrees with results reported in [23].
The heterogeneity of motors within an ensemble does not alter the effect of ATP concentration on run-length.To analyze the impact of ATP concentration on the behavior of ensembles with multiple motors, we characterized the average run-length and velocity for ensembles with two motors (Fig 13   Ensembles with mutant motors also exhibited a similar trend of increasing run-length as ATP concentration was reduced under fixed load forces (Fig 15, WM ensemble).When compared with other two-motor ensembles (Fig 13, load force is 0.2pN), it is seen that at constant load force and ATP concentrations, the cargoes carried by ensembles WW, WM and MM exhibit almost the same average run-length.Thus, if the ATP concentration and number of motors within the ensemble is constant and if the load force does not approach the stall force (or multiples of stall force) for either wild type or mutant (i.e.F s or " F s ), then the average runlength will not be affected by the number of mutant motors in the ensemble.Similar conclusions hold for ensembles with three motors.
When load forces are close to, but less than, the stalling force F s (or its multiples) for the same number of motors in the ensemble, then irrespective of the ATP concentration the ensembles with more wild-type motors will dominate cargo transport and exhibit longer runlengths.In contrast, if load forces are close to, but less than " F s (or its multiples), then ensembles with more mutant motors will dominate cargo transport and exhibit longer run-lengths.

Conclusions
Single motor mutation influences the emergent properties of cargo transport.We developed a simulation engine to study the transport properties of a cargo carried by multiple motor proteins of different types.To demonstrate the efficacy of our methodology we capitalized on a previous study by Selvin and coworkers [25] which characterized a kinesin motor that exhibits a reduced stall force.The study reported that placing a negative charge at the Ser-175 position of the kinesin motor domain, through either mutation or phosphorylation, reduces the stalling force of the motor and decreases its average velocity under external loads.We incorporate this into our method to enable an analysis of the behavior of ensembles comprising of several wild-type and mutant kinesin motors, where the mutant motors have lower stalling forces.The analysis provides support for the following conclusions: 1.A mutation that modulates stalling force of the motors impacts the emergent transport properties of a molecular motor ensemble as follows: a.For a cargo bound to a fixed number of motors and under a constant load force away from the mutant motor stalling force " F s (or its multiples), the average velocity and runlength decreases with an increase in the number of mutant motors in the ensemble.A similar trend is observed for a load force close to the wild-type motor's stalling force F s (or its multiples).Such a behavior is expected, since cargoes carried by ensembles having a higher proportion of mutant motors (that have a reduced stalling force) would traverse shorter runlengths on average.Under these conditions, ensembles with a higher mutant population will eventually lag behind similar-sized ensembles with a lower mutant population.b.Remarkably, for a cargo subjected to load forces close to the mutant motor's stalling force " F s (or its multiples), the average run-length increases with an increase in the number of mutant motors in the ensemble.We hypothesize the existence of a mooring mechanism, where mutant motors tether the cargo to the microtubule, providing the other detached motors with a greater opportunity to reattach to the microtubule and move the cargo forward.It consequently leads to a higher average run-length of the cargo.Thus, surprisingly under these conditions, ensembles with a higher population of mutants gain an aggressive edge over similar sized ensembles with a lower mutant population.These results are similar to the observations in [9,10], where transport of microtubules by varying populations of two kinesin motors, that have a large disparity in average velocities and force dependent dissociation rates, is dictated by a small proportion of faster motors.
2. ATP concentration modulates cargo transport by multi-motor ensembles as follows: a. ATP concentration in specific ranges negatively correlates with average run-length of the cargo under a constant load force, if the cargo has more than one motor attached.This finding agrees with the experimental observation made by Xu and coworkers for a WW ensemble in [23].
b.For a fixed number of motors bound to the cargo, the effect of ATP concentration on run-length is independent of the ratio of wild-type and mutant motors.
3. While we used a kinesin phosphorylation mutant [25] to explore the consequences of a multi-motor ensemble containing motors of lower stall forces, similar analysis and the computational methodology can be utilized to obtain insights on the effect of different types of motors modeled by different sets of parameters on the transport of a common cargo.Moreover, mutations that impact other single motor parameters, including on/off rates, directionality, and elasticity, could be evaluated for impacts on ensemble behavior.
It is likely that the aberrant behavior caused due to the presence of mutant motors with lower stalling forces, when occurring over multitudes of ensembles, can contribute to the disruption of cargo traffic, impediment of neuronal function and the emergence of neurodegeneration.
Simulation engine developed for evaluating the impact of single motor parameters on behavior of motor ensembles.The computational methodology proposed by Gross and coworkers in [22] employed Monte Carlo simulations to offer several insights regarding transport properties of molecular motor ensembles.In this article, we have developed a novel, semianalytical approach to study the effect of single motor mutations on the behavior of motor ensembles and cargo transport.This methodology allows us to compute the biologically relevant quantities such as average run-length and velocity for an ensemble of motors using experimental parameters derived from single motor experiments.The single motor parameters chosen in this article for kinesin-1 motors are obtained from experimental data in [22].The model correctly reproduces the force-velocity curve and run lengths as a function of load for single kinesin-1 motor as reported in [34].The results for ensembles WW and WWW (i.e.all wild-type motors) are in exact quantitative agreement with those obtained via Monte Carlo simulations as presented in [22] as well as [29], fully validating the methodology.For a two-motor ensembles, the results obtained using our model agree with the inferences made in [11,12] using the same (coarse) resolution of load forces.Note that the behavior of multiple motor ensembles around stalling forces (for average run-length and velocity) reported in this manuscript is obtained at a much finer resolution of load forces and remains to be experimentally validated.
Furthermore, unlike Monte Carlo approaches, this model is computationally less extensive and its efficiency is independent of the number of iterations.The model also offers ways to investigate rare events and can be extended to any species of motor proteins, given the knowledge of their individual parameters.It is also easily possible to interrogate specific single motor parameters and to determine their effect on ensemble behavior, thus making it a useful tool to study the contribution of other single motor mutations as well as post translational modifications to the transport behavior of cargoes.

Methods
The Master Equation based methodology used to obtain the aforementioned results is described in this section.This method is used to study emergent properties of an ensemble of multiple motors of two types, that can each take a step on, detach from or reattach onto the microtubule.The knowledge of transition probabilities of stepping, detachment and attachment enable the determination of transition rates between various absolute configurations of the motors, allowing for the calculation of transition rates between the corresponding relative configurations.These rates enable the calculation of the probability distribution of the various ensemble configurations, thereby facilitating the computation of several biologically relevant quantities such as average velocity, run-length and number of attached motors.
We begin with the construction of the relative state space, along with calculations necessary to arrive at several biologically relavent quantities.Then, a general methodology to obtain the transition rates between absolute configurations given the knowledge of the probability rates of stepping, detachment and reattachment to the microtubule for a wild-type motor(P S , P D , P A ) and probability rates for a mutant motor ð " P S ; " P D ; " P A Þ is presented.Finally, the model used to determine the probability rates for kinesin motor proteins is detailed.

Construction of Relative Configuration Space
Consider a cargo that is carried by both wild-type and mutant motor proteins on a microtubule.The microtubule is modeled as directed linear lattice formed by equally sized dimers with dimension d.Here the k th dimer is located at location " a k ¼ kd and indexed by the set of integers I = {. .., −2, −1, 0, 1, 2, . ..}.Each motor protein bound to the cargo can attach, take a for-  The relative configuration of an ensemble of motors is represented using a string of three symbols.Given an absolute configuration we first identify the rearguard motor which is the motor that is attached to the microtubule and lags behind all the other motors on the microtubule.Using the location of the reargurad motor as a reference, the relative configuration ϑ is obtained as a string of three symbols 'M h ' , 'M d ' and '|' , where 'M h ' and 'M d ' denote wild-type and mutant motors respectively, with '|' denoting a separator that distinguishes different microtubule locations.The motor located the farthest from the rearguard motor on the microtubule is identified as the vanguard motor.
For example, the relative configuration of the ensemble in Both the mutant as well as the wild-type motor proteins are characterized via their own set of stepping, attachment, and detachment probabilities(for wild-type (P S , P D , P A ) and for mutant ð " P S ; " P D ; " P A Þ).The individual motors for both species are modeled as hookean springs when stretched that offer no resistance when compressed.A single motor is assumed to have a linkage rest length L 0 and spring stiffness constant K e .Motors of both the species are assumed to not step backward and are bound to the cargo particle irreversibly.More complex models of motor-proteins can be easily accommodated.It is further assumed that the there exists a force F stall called the stalling force, where if the force on the motor protein F ! F stall then the motor does not take a forward step and stepping probability is zero [31].The stalling force can be that is also used to estimate how many motors are carrying the cargo.For an ensemble of wild-type and mutant motors carrying a cargo, the following result holds: Result 1: Given an ensemble of M molecular motors attached to a common cargo that is subjected to a load force F load , the distance between the rearguard the vanguard motor is bound by where F s is the minimum load force for which the stepping probability of the wild-type motor protein becomes zero (i.e. the stalling force for the wild-type motor protein), " F s is the minimum load force for which the stepping probability of the mutant motor protein becomes zero (i.e. the stalling force for the mutated motor protein), L 0 is the rest length of the motor linkage, K e is the linkage stiffness and d is the step-size of the motor.
A detailed derivation is provided in the S1 Text.
It is to be noted that the absolute configuration space admits infinitely many representations as there is always a small probability of finding the cargo at any location on the microtubule.However, the above result concludes, that given a stall force for both wild-type and mutant motors the relative configuration space is finite, since there are no motors beyond n units away from the rearguard motor in any relative configuration.
Transition probabilities between relative configurations.Next, we determine the probability of transitioning from a relative configuration to another from the transition probabilities in the absolute configuration space.
It is evident that an absolute configuration can be associated with a unique relative configuration.For example, the relative configuration associated with Fig 16 We denote the projection operator Υ that maps an absolute configuration O to a corresponding relative configuration by Υ(O).The probability that the absolute configuration is O 0 at time t + Δt conditioned on an initial configuration of O at t is represented as P O (O 0 , t + Δt|O, t).It is assumed that the transition probability P O (O 0 , O) between the two absolute configurations O and O 0 is given by ν O (O 0 , O) Δt for a small time interval Δt, where the notation ν O (O 0 , O) represents the probability rate of transition between O and O 0 .The underlying assumption that the rate ν O (O 0 , O) is independent of the time instant t holds true, since the ν O only depends upon the initial configuration and the type of transition (motor stepping, detachment or attachment) from O to O 0 .
The knowledge of the transition rate ν O in the absolute configuration space enables the following result: Result 2: The rate of transition between relative configurations ϑ and ϑ 0 is given by: where O is any absolute configuration that satisfies Υ(O) = ϑ, O 0 is any absolute configuration that satisfies Υ(O 0 ) = ϑ 0 , (τ β O 0 ) is an absolute configuration obtained after linearly shifting all the motors in O 0 by β locations on the microtubule towards the right and d is the dimension of a single microtubule dimer A detailed derivation of Eq (2) is provided in S2 Text.
As n is finite for a finite number of molecular motors attached to the cargo (from Eq (1)), Eq (2) involves only a finite number of computations.The transition rates between absolute configurations are obtainable using the chemical kinetics of the motor protein as it steps, detaches or attaches to the microtubule.The rate ν ϑ (ϑ 0 , ϑ) is obtainable using Result 2 and the knowledge of the rates ν O between corresponding absolute configurations.
Results 1 and 2 together imply that given a finite number of molecular motors carrying a common cargo, the relative configuration space is finite.Thus, by projecting or aggregating the states of the infinite dimensional model (consisting of absolute configuration O 0 s) into 'macrostates' , one can obtain a finite dimensional model (see [8] for related work).Thus, if H is the set of all the possible relative configurations, the knowledge of transition rates from Eq (2) can be used to show that the probability P ϑ (ϑ, t) of the relative configuration being ϑ at time t satisfies the master equation, Evolution of probability distribution of relative configuration.Solution for the master Eq (3) determines the time evolution of probabilities of all the relative configurations.Note that Eq (3) is an exact description of how the probability distribution evolves in the relative configuration space without any approximations.Consider a ordering of all relative configurations given by ϑ 1 , . .., ϑ N (where N is the finite number of relative configurations) where the probability of finding the motors in a relative configuration ν i at time t is denoted by P i (t) and let P(t) = [P 1 (t), . .., P N (t)] T .Using the expression of transition rates n W ðW j ; W 0 i Þ and Eq (3), it can be shown that the dynamics of the model describing the vector P(t) is given by: where G 2 R NÂN is a sparse stochastic matrix which is determined by the transition rates . Given a specified initial probability vector P(t 0 ), it follows that When specified for kinesin-1 motors and realistic values of number of motors (M 5) and system parameters, the dimension of Γ is N × N where N lies between 10 5 and 10 7 , making the evaluation of e Γ(t − t 0 ) manageable using a standard computer(the results in this article are obtained using Intel quad core i5 processor, 3.4 GHz, RAM 8 GB).The software is easily implementable using platforms like MATLAB and is faster and computationally more efficient than Monte-Carlo based approaches.Computer clusters can be used to manage more complex scenarios involving larger number of motors.
Calculating biologically relevant quantities.Note that in obtaining the transition rates between the relative configurations, some amount of information pertaining to the transitions between the absolute configurations is lost.However as alluded to earlier, the solution to Eq (3) provides an exact description of the probability vector of P(t), which can then be used to obtain expressions of several biologically relevant quantities in the following manner.
Average where d avg (ϑ 0 , ϑ) is the expected change in cargo position, when the initial and final relative configurations at t and t + Δt are restricted to being ϑ and ϑ 0 respectively and is given by The expression d(O 0 , O) in Eq (8) is the change in cargo equilibrium position when the absolute configuration changes from O to O 0 .The detailed formulation of Eqs ( 7) and ( 8) is provided in S3 Text.Average run-length is then calculated by summing the average velocity over time Transition Rates between Absolute Configurations In this section we present a general scheme for determining transition rates between absolute configurations.We begin with a structural model for single motor protein that consists of motor head, stalk and cargo binding tail domain.The linkage between the motor-heads and tail for single motor is modeled as a hookean spring when stretched, that has a rest length L 0 .It offers no resistance when compressed [22].The motor heads move along the microtubules exerting a force F on a cargo that is expressed as a function of its length L by, Z eq is the mean position of the cargo that is its equilibrium position determined by the forces exerted by the motors on the cargo through their linkages and the load force F load on the cargo.The mean cargo position for a fixed F load > 0 is a function of the absolute configuration i.e.Z eq = Z eq (O).If the cargo position is assumed to follow a truncated Gaussian distribution Θ(z) with variance σ, its probability density Θ(z) for |z| < 3σ is given by, The effect of thermal noise can be incorporated by determining the steady state variance σ of the cargo position.We assume that when a motor in the ensemble takes a forward step or detaches, the probability distribution of the cargo position reaches a new distribution with negligible transient.Here, the time scale of the cargo dynamics is assumed to be faster than the rate of change of motor configurations.
A transition to another configuration O 0 occurs if either the wild type or mutant motor at a location " a k steps forward to " a kþ1 , detaches from the location " a k or reattaches to the location " a k on the microtubule.By representing O 0 as O + S, S is a sequence that corresponds to the type of transition(step, detach or attach) and the type of motor(wild-type or mutant) that has transitioned.The transition rate from O to O 0 is determined by averaging the associated probability rate over the position of the cargo.
The model of a single motor-protein is specified via the probability P S (F) of the motor taking a step, the detachment probability, P D (F), of the motor detaching from the microtubule, and the probability of attachment P A of an unattached motor-protien to the microtubule, per second.Here F is the force acting on the motor which is considered positive if it is directed opposite to the motor stepping direction (e.g.kinesin forward stepping is towards the mictorubule + end).Here in order to calculate the transition rates between absolute configurations it is assumed that the probability rates of step, detachment and attachment are known; later we illustrate a way to compute these probabilities for kinesin motors.
Transition rate for stepping.The stepping transition of a wild-type motor from the location " a k to " a kþ1 is represented as where O h,k is the number of wild-type motors located at " a k and P S is the probability of a wild type motor taking a step from the location " a k on the microtubule per second.The stepping transition of a mutant motor from the location " a k to " a kþ1 is represented as where the plus sign is for attachment and minus sign is for detachment.where O d,k is the number of mutant motors located at " a k and " P D is the probability of a mutant motor detaching from the location " a k on the microtubule per second.For both the wild-type and mutant motors, a constant probability rate of reattachment to the microtubule, P A and " P A , are assumed.If a motor is linked to the cargo, it is assumed to attach to the microtubule without stretching its linkage.Thereby, the only locations where the motor reattachment can occur are located within a distance L 0 of the cargo.

Probability Rates for Kinesin
The results put forth in this article correspond to an instantiation of our methodology for kinesin motor protein.The probability rates of stepping, detachment and attachment of a single kinesin motor are determined using several available studies [16,21,30,31].
Probability of stepping, per second.Kinesin takes a step on the microtubule by hydrolyzing an ATP molecule [32] in the following manner: where ΔE is the energy released.From [33], the ATP hydrolysis rate predicted using Michaelis-Menten dynamics relates to the probability rate stepping for a wild-type kinesin by assuming that the free motor head binds to the microtubule location with defined probability η.P S is then expressed as where k m ¼ k cat þk off k on : From [22], the force F exerted by the cargo on the motor affects motor dynamics by modifying η as, Furthermore [22] assumes that the force F influences the kinetics of ATP hydrolysis by affecting k off in the following manner: where k 0,off is the backward reaction rate of hydrolysis when F = 0, δ l is an experimentally determinable parameter, T is the temperature and K b is the Boltzmann constant.Thus, the probability rate of stepping for a wild-type motor under a constant force F is given by, In the case of a mutant motor protein with an altered stalling force " F s , the probability rate of stepping under a force F is given by, where the probability " Z that the free mutant motor head binds to the microtubule location is written as Probability of detachment, per second.Kinesin motors are processive species that take a certain number of steps during their ATP driven movement along microtubule before dissociating from it.The processivity (denoted by L) represents how far the motor can move, on average, before its detachment from microtubule.From definition of L in [34], where A, B and δ l are experimentally determinable parameters.The relation between the probability rate of stepping P S and detachment P D is given by, As long as the wild-type motor is not stalled i.e.F < F s , the probability rate of detachment is given by, When the load force on the wild-type motor equals to or exceeds its stalling force i.e. for F ! F s , a constant detachment rate is assumed in [22] as, P D ðFÞ ¼ P back ¼ 2=sec: For a mutant motor that is not stalled i.e.F < " F s , the probability rate of detachment is given by, When the load force on the mutant motor equals to or exceeds its stalling force i.e. for F ! " F s , the constant detachment rate is assumed to be " P D ðFÞ ¼ P back ¼ 2=sec: For F ! F s or F ! " F s , an assumption that the single motor detachment rate increases exponentially with F as opposed to a discontinuous jump to a constant value of 2/sec, does not alter the nature of the results obtained.Furthermore, assuming that the stepping probability of the single motor under load force equal to its stall force does not become zero but a small non-zero value, only reduces the sharpness of the peaks in Figs 9 and 10; it also does not affect the nature of the results reported.The outcomes remained unchanged even after added simulations were performed (not shown here) to test the robustness of our results, where different single motor models were used while retaining the same fundamental structural assumptions as the single motor model used in this manuscript.
Probability of attachment, per second.Probability of attachment is experimentally found to be P A % 5/sec [35,36].In this article it is assumed that the probability of attachment for a mutant motor, " P A % 5=sec.Numerical parameters for Kinesin.The numerical parameters considered for wild-type motors in this article are the same as that used in [22].Specifically k on = 2.10 6 M −1 s −1 , k off = 55 s −1 , k cat = 105 s −1 , F s = 6 pN, d = 8 nm, d l = 1.6 nm, δ l = 1.3 nm, A = 107, B = 0.029 μM, [ATP] = 2 mM, T = 300K and K el = 0.32 pN/nm.All the parameters for mutant motors except for the stalling force, are assumed to not deviate significantly from the parameter set for wildtype motors.The stalling force for mutant motor is taken as " F s ¼ 5:5 pN.
Fig 1(a) shows the probability distribution computed for the WM case, where with increasing time the probability of no motors attached (red curve) increases and eventually reaches one.

Fig 1 .
Fig 1. Probability of motors being bound to the microtubule at time t under load force of 1.2 pN for the WM case.In (a) the probability of one bound motor (blue) and two bound motors (green) approaches 0 as time t ! 1. Probability of no bound motors (red) (i.e. the cargo dissociates from the microtubule) approaches 1 as t ! 1.In (b) the probability of one bound motor (blue) and two bound motors (red) under the condition that cargo is not lost, reaches a steady state.doi:10.1371/journal.pcbi.1005152.g001

Fig 3 .Fig 4 .
Fig 3. Probability of the cargo-motor complex getting detached from the microtubule, for ensembles WW, WM and MM at load forces close to but less than F s and 2F s (F s = 6 pN).doi:10.1371/journal.pcbi.1005152.g003 Fig 8 shows the average velocity for load forces between 10.4 pN and

Fig 5 .
Fig 5. Probability of the cargo-motor complex getting detached from the microtubule, for ensembles WW, WM and MM at load forces 3 pN and 8 pN.doi:10.1371/journal.pcbi.1005152.g005

Fig 8 .Fig 10 .
Fig 8. Fig showing variation of average velocity with load force on the cargo for ensembles containing 3 motors, for load forces between 10.4 pN and 12.1 pN.doi:10.1371/journal.pcbi.1005152.g008

Fig 11 .
Fig 11.Probability of wild-type and mutant motor being attached for ensemble WM. doi:10.1371/journal.pcbi.1005152.g011 13) and WWW (S3 Fig) ensembles, we observe that with reducing ATP concentration, the run-length first increases, peaks at a certain value and then decreases.The average velocity of cargo transport will decrease with decreasing ATP concentration (see S1 Fig for WW and S2 Fig for WWW) and S1 Fig) and three motors (S3 and S2 Figs).

Fig 12 .
Fig 12.Effect of ATP concentration on run-length for ensemble containing one wild-type motor, for various hindering load forces on the cargo.doi:10.1371/journal.pcbi.1005152.g012

Fig 15 .
Fig 15.Effect of ATP concentration on run-length for ensemble WM, for various hindering load forces on the cargo.doi:10.1371/journal.pcbi.1005152.g015 ward step or detach from the microtubule.The absolute configuration O ≔ O h;k arrangement on the microtubule specifies the number, O h,k , of wild-type motor proteins and the number, O d,k , of mutant motor proteins at the k th location on the microtubule.For example, the absolute configuration of the ensemble of motors illustrated in Fig 16 is given by

Fig 16 .
Fig 16.Locations of wild-type and mutant motors on a section of the microtubule lattice.The microtubule is modelled as a directed linear lattice " a k ¼ " a 0 þ kd, where " a k is the position of the k th location.If all the motors carry a common cargo i.e. the wild-type motor at a 0 is the rearguard motor, the string representation is |M h ||M h M d ||M d ||M h | doi:10.1371/journal.pcbi.1005152.g016 Fig 17(a) is the string '|M h M d || M d ||M h |' .The configuration that results after the furthermost mutant motor in Fig 17(a) takes a step is shown in Fig 17(b) which has a relative configuration given by '|M h M d |||M d |M h |' .

Fig 17 .
Fig 17.Locations of wild-type and mutant motors on a section of the microtubule lattice.The string representation for the two configurations is (a) |M h M d ||M d ||M h | and (b) |M h M d |||M d |M h | doi:10.1371/journal.pcbi.1005152.g017

:
As the transition O !ATT d =DET d O AE S ðattÞ d;k occurs if any of the O d,k mutant motors located at " a k detach from the microtubule, the associated transition rate for detachment, n O ðO À S ðattÞ d;k ; OÞ is determined by averaging over the position of the cargo as, n O ðO À S ðattÞ d;k ; OÞ ¼ O d;k Z Z eq ðOÞþ3s Z eq ðOÞÀ 3s " P D ðFðz À " a k ÞÞYðz À Z eq ðOÞÞdz

Number of Engaged Motors (Wild-type/Mutant)
Average number m h (t) of wild-type motors attached to the microtubule is given by (ϑ i ) represents the number of wild-type motors in the relative configuration ϑ i or the number of M h symbols in the representation ϑ i .In a similar manner, the average number of mutant motors m d (t), attached to the microtubule is given by (ϑ i ) represents the number of mutant motors in the configuration ϑ i or the number of M d symbols in the representation ϑ i .
h;k occurs if any of the O h,k wild-type motors located at " S ðstepÞ d;k occurs if any of the O d,k mutant motors located at " a k step forward to the position " a kþ1 , the associated transition rate for stepping, n O ðO þ S ðstepÞ d;k ; OÞ is determined by averaging over the position of the cargo as, As the transition O !ATT h =DET h O AE S ðattÞ h;k occurs if any of the O h,k wild-type motors located at " a k detach from the microtubule, the associated transition rate for detachment, n O ðO À S ðattÞ h;k ; OÞ is calculated by averaging over the position of the cargo as, ÞÞYðz À Z eq ðOÞÞdz where O h,k is the number of wild-type motors located at " a k and P D is the probability of a wild type motor detaching from the location " a k on the microtubule per second.The attachment/detachment transition of mutant motor at location " a k is represented as