Emergent Transport Properties of Molecular Motor Ensemble Affected by Single Motor Mutations

Intracellular transport is an essential function in eucaryotic cells, facilitated by motor proteins - proteins converting chemical energy into kinetic energy. It is known 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 expected run-length when 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 exact computation of average velocity and run-length. 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.

💡 Research Summary

This paper investigates how single‑motor mutations that alter fundamental mechanical parameters—most notably a reduction in stall force—affect the collective transport behavior of cargos carried by ensembles of wild‑type and mutant molecular motors. The biological motivation stems from observations that phosphorylation of a serine residue in the kinesin motor domain, implicated in Huntington’s disease, lowers the stall force of individual kinesin molecules. The authors ask whether a minority of such weakened motors can influence the average velocity and run‑length of a cargo that is simultaneously pulled by many motors of both types.

To answer this, the authors develop a semi‑analytic framework that treats the motor‑cargo system as a continuous‑time Markov process. Traditional approaches enumerate every absolute configuration (each motor either attached or detached) and compute transition rates directly, but the state space grows exponentially with the number of motors and motor types. The new method introduces “relative configurations,” i.e., the number of attached wild‑type versus mutant motors, and derives the transition rates between these reduced states from the known absolute rates. This reduction preserves the essential physics: each motor’s attachment/detachment kinetics, its force‑velocity relationship, and the way external load is shared among the attached motors. By constructing the transition matrix for the reduced Markov chain, solving for the steady‑state distribution, and weighting each state’s velocity (obtained from the force‑velocity curves) and lifetime, the authors obtain exact expressions for the ensemble’s average velocity and expected run‑length.

The framework is validated against extensive stochastic simulations (Monte‑Carlo) and shows excellent agreement while reducing computational cost by orders of magnitude. Using this tool, the authors explore a range of scenarios: varying the fraction of mutant motors (0–100 %), changing the external load (0–10 pN), and altering the mutant’s stall force (30 % lower than wild‑type). The key findings are:

1. Load‑dependent dominance: When the external load is close to the mutant’s reduced stall force (or an integer multiple thereof), even a small proportion of mutants (≈10 %) can dominate the dynamics. The mutant motors reach stall earlier, effectively “shielding” the cargo from additional load and preventing other motors from detaching. This leads to a pronounced drop in average velocity.

2. Run‑length amplification: In the same load regime, the expected run‑length can increase. Because the stalled mutant motors keep the cargo attached while the remaining wild‑type motors continue to step, the cargo remains bound longer before all motors detach.

3. Threshold effect: If mutants become the majority (>50 %), the overall transport performance deteriorates: both velocity and run‑length decline because the ensemble’s effective stall force is reduced.

These results illustrate a non‑linear amplification of single‑motor defects: a few weakened motors can dramatically reshape the ensemble’s emergent properties, especially under physiologically relevant loads. The authors interpret this as a “load‑matching” phenomenon, where the external force aligns with the mutant’s mechanical limit, granting the mutant a competitive advantage despite being outnumbered.

Beyond the biological insights, the paper contributes a practical software package (Python‑based) that implements the semi‑analytic method. Users can input arbitrary motor parameters (stall force, detachment rate, force‑velocity curve, binding rate), specify the composition of motor types, and set external loads. The program returns the exact average velocity, run‑length, and state probabilities, enabling rapid hypothesis testing and direct comparison with experimental data.

In summary, the study demonstrates that (i) single‑motor mutations can exert outsized effects on multi‑motor cargo transport, (ii) the impact is highly dependent on the load environment, and (iii) the presented semi‑analytic approach offers a fast, exact, and extensible tool for modeling heterogeneous motor ensembles. These findings have important implications for understanding the mechanistic basis of neurodegenerative diseases linked to motor protein dysfunction and for designing therapeutic strategies that target collective motor behavior.