Stochastic kinetics of a single headed motor protein: dwell time distribution of KIF1A

KIF1A, a processive single headed kinesin superfamily motor, hydrolyzes Adenosine triphosphate (ATP) to move along a filamentous track called microtubule. The stochastic movement of KIF1A on the track is characterized by an alternating sequence of pause and translocation. The sum of the durations of pause and the following translocation defines the dwell time. Using the NOSC model (Nishinari et. al. PRL, {\bf 95}, 118101 (2005)) of individual KIF1A, we systematically derive an analytical expression for the dwell time distribution. More detailed information is contained in the probability densities of the “conditional dwell times” $\tau_{\pm\pm}$ in between two consecutive steps each of which could be forward (+) or backward (-). We calculate the probability densities $\Xi_{\pm\pm}$ of these four conditional dwell times. However, for the convenience of comparison with experimental data, we also present the two distributions $\Xi_{\pm}^{*}$ of the times of dwell before a forward (+) and a backward (-) step. In principle, our theoretical prediction can be tested by carrying out single-molecule experiments with adequate spatio-temporal resolution.

💡 Research Summary

The paper presents a comprehensive stochastic analysis of the dwell‑time distribution of the single‑headed kinesin motor KIF1A as it walks along microtubules. KIF1A moves by hydrolyzing ATP, and its trajectory consists of alternating pauses and directed steps. The authors adopt the NOSC model originally proposed by Nishinari et al. (Phys. Rev. Lett. 95, 118101, 2005), which describes the motor as alternating between two chemical states (strongly and weakly bound to the filament) with four kinetic rates (k₁–k₄) that depend on ATP concentration and filament interaction strength.

Starting from the master equations for the probability of finding the motor in each state‑position combination, the authors perform a Laplace‑transform analysis to obtain exact analytical expressions for the time‑dependent probabilities. By integrating over the two sub‑intervals that compose a full dwell (the pause and the subsequent translocation), they derive the probability density function (PDF) for the total dwell time τ. Crucially, they distinguish four conditional dwell times: τ₊₊ (forward step followed by another forward step), τ₊₋ (forward then backward), τ₋₊ (backward then forward), and τ₋₋ (backward then backward). For each case they obtain a closed‑form PDF, denoted Ξ₊₊(t), Ξ₊₋(t), Ξ₋₊(t), and Ξ₋₋(t). These PDFs are typically expressed as sums of one or two exponential terms whose decay constants are the eigenvalues of the kinetic matrix. The relative weights of the exponentials are functions of the kinetic rates, allowing a direct mapping from experimentally measured dwell‑time histograms to underlying molecular parameters.

To facilitate comparison with single‑molecule data, the authors further combine the conditional PDFs into two experimentally relevant distributions: Ξ₊(t), the dwell‑time distribution preceding a forward step, and Ξ₋(t), the dwell‑time distribution preceding a backward step. The combination is performed by weighting the conditional PDFs with the probabilities that a given conditional sequence actually occurs (e.g., p₊₊ = k₁/(k₁+k₂)·k₃/(k₃+k₄) for a forward‑forward sequence). Consequently, Ξ₊(t) and Ξ₋(t) are also represented by simple exponential mixtures, whose parameters can be extracted by fitting to experimental dwell‑time histograms.

The paper includes a systematic parameter study. Increasing ATP concentration raises k₁, shortening the average forward‑preceding dwell time and increasing the forward‑step probability. Raising the backward transition rate k₄ lengthens the backward‑preceding dwell time and reduces overall processivity. The authors present analytical formulas for the mean dwell time ⟨τ⟩ and its variance for each conditional case, illustrating how the shape of the distribution (e.g., the presence of a long tail) reflects the balance between forward and backward kinetic pathways.

In the discussion, the authors acknowledge the simplifying assumptions of the NOSC framework: a single head, independent transitions, and a homogeneous microtubule lattice. They outline possible extensions to incorporate multiple heads, external loads, lattice defects, and non‑linear thermal noise. Importantly, they propose experimental protocols—high‑speed total internal reflection fluorescence (TIRF) microscopy or optical tweezers with millisecond resolution—to directly measure the four conditional dwell times. By fitting the measured PDFs to the derived analytical forms, one can infer the kinetic rates k₁–k₄, thereby providing a quantitative link between molecular biochemistry and motor mechanics.

Overall, the study delivers an exact, analytically tractable description of KIF1A dwell‑time statistics, bridges theory with experimentally accessible observables, and offers a robust framework for testing motor‑protein models and for designing drugs that modulate kinesin activity.