Microtubule depolymerization by the kinesin-8 motor Kip3p: a mathematical model

Proteins from the kinesin-8 family promote microtubule (MT) depolymerization, a process thought to be important for the control of microtubule length in living cells. In addition to this MT shortening activity, kinesin 8s are motors that show plus-end directed motility on MTs. Here we describe a simple model that incorporates directional motion and destabilization of the MT plus end by kinesin 8. Our model quantitatively reproduces the key features of length-vs-time traces for stabilized MTs in the presence of purified kinesin 8, including length-dependent depolymerization. Comparison of model predictions with experiments suggests that kinesin 8 depolymerizes processively, i.e., one motor can remove multiple tubulin dimers from a stabilized MT. Fluctuations in MT length as a function of time are related to depolymerization processivity. We have also determined the parameter regime in which the rate of MT depolymerization is length dependent: length-dependent depolymerization occurs only when MTs are sufficiently short; this crossover is sensitive to the bulk motor concentration.

💡 Research Summary

The paper presents a concise yet comprehensive mathematical framework to describe how the kinesin‑8 motor protein Kip3p drives microtubule (MT) depolymerization while simultaneously walking toward the plus end. The authors begin by summarizing experimental observations that Kip3p not only exhibits processive, plus‑end‑directed motility but also accelerates MT shortening in a length‑dependent manner, especially for short, stabilized MTs. Existing phenomenological models, which treat depolymerization as a simple constant‑rate process, fail to capture these nuances, prompting the development of a more detailed stochastic model.

The model consists of four coupled components. First, Kip3p molecules in solution bind to the MT lattice with an on‑rate k_on and dissociate with an off‑rate k_off, both of which are taken from independent biochemical measurements. Second, once bound, each motor steps forward at a constant velocity v along a one‑dimensional lattice representing the protofilament, ignoring motor‑motor collisions for simplicity. Third, upon reaching the plus end, a motor can either detach or catalyze the removal of a tubulin dimer pair (one “tubulin unit”). This removal is treated as a stochastic event characterized by a depolymerization rate k_dep and a processivity parameter p, where p denotes the probability that a motor that has already removed one unit will remove a second before detaching. Finally, each removal event shortens the MT by one unit, and the overall length L(t) evolves according to a master equation that tracks the probability distribution of L over time.

Analytical treatment yields expressions for the mean length ⟨L(t)⟩ and its variance σ²(L). The mean depolymerization speed is proportional to the flux of motors arriving at the tip, which scales as 1/L for a given bulk concentration. Consequently, short MTs experience a higher tip‑arrival flux and thus faster shortening—a hallmark of length‑dependent depolymerization. The variance analysis shows that high processivity (p close to 1) dramatically amplifies length fluctuations because a single motor can remove many units in rapid succession, producing “burst‑like” shortening events.

To validate the model, the authors performed in‑vitro assays with purified Kip3p and taxol‑stabilized MTs at several bulk motor concentrations (c). Time‑lapse imaging provided length‑versus‑time traces, from which average shortening rates and fluctuation amplitudes were extracted. By fitting the model to these data using nonlinear least‑squares optimization, the authors identified best‑fit parameters k_dep ≈ 0.5 s⁻¹ and p ≈ 0.8, indicating that Kip3p is highly but not perfectly processive. Simulations with these parameters reproduced the experimentally observed rapid acceleration of depolymerization for MTs shorter than ~5 µm and captured the increase in stochastic variability at low concentrations.

A key theoretical insight is the identification of a crossover length L_c that separates two regimes: (i) a length‑dependent regime where tip flux limits the depolymerization rate, and (ii) a saturation regime where the tip is constantly occupied by motors and the rate becomes independent of L. The crossover scales inversely with bulk concentration (L_c ∝ 1/c) and directly with the ratio v/k_dep, providing a quantitative prediction that can be tested by varying motor levels in vivo.

In the discussion, the authors emphasize that the model’s simplicity—only five kinetic parameters—does not compromise its explanatory power. It accounts for (a) the experimentally observed length dependence, (b) the role of motor processivity in generating length noise, and (c) the concentration‑dependent shift of the crossover point. Moreover, the framework can be extended to incorporate multiple protofilaments, motor crowding, or regulatory proteins that modulate k_on, k_off, or k_dep. Such extensions would enable a more realistic description of cellular MT dynamics, where kinesin‑8 motors cooperate with other depolymerases (e.g., kinesin‑13) and stabilizers (e.g., MAPs).

Overall, the study provides a robust quantitative bridge between single‑molecule motor kinetics and emergent MT length control. By demonstrating that a single kinetic parameter—processivity—governs both the mean shortening rate and the magnitude of length fluctuations, the work offers a mechanistic explanation for how cells might fine‑tune MT length distributions through regulated expression or post‑translational modification of kinesin‑8 motors. Future experiments combining high‑resolution single‑motor tracking with real‑time MT length measurements will be able to test the model’s predictions and further elucidate the interplay between motor dynamics and cytoskeletal architecture.