Stall force measurement of the kinesin-3 motor KIF1A using a programmable DNA origami nanospring

  1. Nobumichi Takamatsu
  2. Hiroko Furumoto
  3. Takayuki Ariga
  4. Mitsuhiro Iwaki  Is a corresponding author
  5. Kumiko Hayashi  Is a corresponding author
  1. The Institute for Solid State Physics, The University of Tokyo, Japan
  2. Department of Complexity Science and Engineering, Graduate School of Frontier Sciences, The University of Tokyo, Japan
  3. Systems Biochemistry in Pathology and Regeneration, Yamaguchi University Graduate School of Medicine, Japan
  4. Graduate School of Frontier Bioscience, The University of Osaka, Japan
  5. Advanced ICT Research Institute, National Institute of Information and Communications Technology, Japan
  6. Immunology Frontier Research Center (IFReC), The University of Osaka, Japan

eLife Assessment

Optical tweezers have been instrumental to the determination of mechanical parameters of molecular motors. This study by Takamatsu et al. reports key mechanical parameters of kinesin KIF1A using fluorescence microscopy, wherein the motor is tethered to a DNA nanospring, without the use of an optical trapping apparatus, which represents an exciting development. The approach and the findings reported change current thinking about KIF1A‑mediated transport, with potential implications for understanding human disease. The findings are important and the strength of the evidence is compelling.

https://doi.org/10.7554/eLife.108477.3.sa0

Abstract

DNA origami technology is a method for designing and constructing nanoscale structures using DNA, and it is being applied across various fields. This technology was advanced by developing the nanospring (NS), a fluorescently visible molecular spring that quantifies forces through its extension and has been used to measure myosin-generated forces. This study aims to measure the force exerted by the kinesin-3 motor protein KIF1A, mutations of which cause KIF1A-associated neurological disorder (KAND) and are associated with reduced force and motility. Unlike kinesin-1, KIF1A detaches easily under perpendicular loads, which can occur in optical tweezers experiments. By applying force parallel to the microtubule using the NS, we were able to precisely measure the stall force even for KAND mutants, for which such measurements are typically challenging. This result highlights the potential of the NS as a new tool for force spectroscopy in biophysics.

Introduction

The kinesin superfamily of microtubule-based motor proteins in humans consists of approximately 45 types and represents a large group of motor proteins with diverse functions (Hirokawa et al., 2010; Vale, 2003). Among them, KIF1A, a member of the kinesin-3 subfamily, plays a critical role in the long-distance transport of synaptic vesicle precursors and other cargos in neurons (Gabrych et al., 2019). This transport is essential for maintaining synaptic activity, supporting neuronal development, and ensuring proper neuronal function. Structurally, KIF1A contains a motor domain that binds to microtubules and hydrolyzes ATP to generate mechanical force, driving movement toward the plus end of microtubules. In recent years, KIF1A has gained significant attention due to the identification of over 100 mutations linked to KIF1A-associated neurological disorder (KAND), a spectrum of neurodevelopmental conditions characterized by spastic paraplegia, intellectual disability, and progressive neurodegeneration (Kaur et al., 2020; Boyle et al., 2021).

Single-molecule measurements of KIF1A mutants associated with KAND have developed (Budaitis et al., 2021; Lam et al., 2021; Anazawa et al., 2022; Chiba et al., 2019). Precise measurements of the velocity of KIF1A single molecules have been achieved using TIRF microscopy, while optical tweezers have enabled force measurements. These techniques have been used to compare the biophysical properties of KIF1A single molecules with those of its mutants that cause KAND. The motility of the single molecules has been reported to correlate with disease severity, highlighting the importance of single-molecule biophysical measurements into precision medicine (Rao et al., 2025). The velocity of KIF1A has been extensively studied; however, there is still room for improvement in force measurements using optical tweezers. To date, optical tweezers have been used to measure the stall force, defined as the maximum force that a motor can generate when its movement is halted, of various microtubule-based motor proteins (Schnitzer et al., 2000; Gennerich et al., 2007; Elshenawy et al., 2019; Mallik et al., 2004; Tomishige et al., 2002; Okada et al., 2003; Nishiyama et al., 2002). Since a submicron-sized bead is attached to the kinesin (~10 nm) and force is applied via the bead to the kinesin in this method, the applied force tends to pull it upward, which can sometimes lead to its detachment from the microtubule (Pyrpassopoulos et al., 2020; Pyrpassopoulos et al., 2023). Since KIF1A is easier to detach from microtubules compared to conventional kinesin-1, it is more challenging to successfully stall KIF1A in optical tweezer experiments without detachment. Indeed, most of the stalls lasted less than 1 s, and the force trajectories exhibited sawtooth-like patterns in previous studies (Budaitis et al., 2021; Lam et al., 2021).

To overcome the challenges associated with force measurements in KAND mutants, we employed an alternative approach based on DNA nanotechnology. In a previous study, Iwaki et al. developed a distinct method for force measurement by using programmable DNA origami to construct nanoscale, coil-like springs known as nanosprings (NSs) (Iwaki et al., 2016; Matsubara et al., 2023). These NSs are fluorescently labeled, allowing their extension to be visualized with conventional fluorescence microscopy. Since the force–extension relationship of the NS is pre-calibrated using acoustic force spectroscopy (AFM; Matsubara et al., 2023), the force exerted on the spring can be quantitatively estimated from its measured extension. In this study, we applied the NS method to measure the force generated by single KIF1A molecules, including both wild-type and KAND mutants. Using the NS system, we were able to apply a horizontal load parallel to the microtubule axis to single KIF1A molecules and clearly observed sustained stalling for several tens of seconds with reduced detachment. Notably, stalling events were also stably observed in KAND mutants such as P305L and V8M, which disrupt critical force-generating regions of the protein (e.g., the K-loop and neck linker) (Budaitis et al., 2021; Lam et al., 2021). We also included the A255V mutant, previously unmeasured in force assays, which is known to affect ATP hydrolysis (Budaitis et al., 2021). Furthermore, the NS allowed for precise stall force measurement, revealing force differences between KAND mutant homodimers and heterodimers composed of a wild-type and the mutant subunit. These results highlight the advantage of our approach in minimizing the force component perpendicular to the microtubule’s long axis.

In previous single-molecule experiments in kinesin, DNA origami was used to control the number of cooperatively moving kinesin motors, and their velocities have been analyzed (Furuta et al., 2013; Derr et al., 2012). As a distinct application of DNA-origami-based motor protein systems, in this study, we successfully measured force using only fluorescence imaging. This study is significant in that it extends fluorescence imaging from merely capturing velocity to also allowing the quantification of force. In the field of motor proteins research, measurements of physical quantities related to distance, such as velocity and stepping motion, have rapidly advanced with the development of cutting-edge microscopy techniques like a super-resolution microscopy technology MINFLUX (Balzarotti et al., 2017), enabling their investigation not only in glass chambers but also inside living cells (Deguchi et al., 2023). However, since both force and velocity are crucial physical parameters to understand the functions of motor proteins, and force–velocity relationships are indeed related to their ATP hydrolysis mechanisms (Schnitzer et al., 2000; Sasaki et al., 2018), we believe that the development of a new type of force measurement is of great significance in this field. The development of NS-based force measurement, which enables force quantification at the nanoscale, is expected to have significant implications not only in the field of motor protein research but also in the broader study of protein biophysics.

Results

The NS–KIF1A–inert KIF5B complex as a force measurement system

To measure the force generated by single KIF1A molecule, we constructed a complex comprising a KIF1A molecule, an inert KIF5B mutant, and an NS. For KIF1A, the motor domain (amino acid residues 1–393) was used, excluding the C-terminal and cargo-binding domains (Figure 1a). Since KIF1A does not form stable dimers in the absence of cargo-binding domains, a leucine zipper domain was incorporated to stabilize KIF1A dimers, mimicked the activated state of full-length KIF1A, as in the previous study (Budaitis et al., 2021). As for the inert KIF5B attached to the one end of the NS, it was engineered to remain immobilized on microtubules by introducing the G234A mutation (Methods). This mutant lacks ATP hydrolysis activity and, as a result, does not exhibit motility, allowing it to remain stably attached to microtubules (Rice et al., 1999). Both kinesins were equipped with SNAP-tags at their N-termini, which enabled chemical coupling to the NS (Figure 1b; Methods). The inert KIF5B mutant served to anchor the KIF1A–NS complex to a microtubule. This anchoring of the KIF5B to the microtubule allowed a KIF1A to stay near the microtubule, increasing the probability of its binding to the microtubule.

Figure 1 with 2 supplements see all
Experimental design for the stall force measurement of KIF1A using an NS.

(a) Schematic of the domain structure of full-length KIF1A and the recombinant construct used in the experiments. To stabilize KIF1A dimers, which do not form stably without cargo-binding domains, a leucine zipper was incorporated, and SNAP-tags were added at the N-termini to enable chemical coupling to the NS. (b) Schematic of the chemical modifications required for coupling the NS to kinesin (Methods). (c) An inert KIF5B is anchored to the microtubule, and the NS extends as a KIF1A moves toward the plus end. The NS is uniformly labeled with Cy3 fluorophores, allowing force to be calculated from its extension. The micrographs depict the NS in the retracted and extended states. Here, the microtubule axis is defined as the x-direction, and the direction perpendicular to the microtubule is defined as the z-direction. (d) Force–extension relationship of the NS, showing nonlinear elastic behavior, calibrated by acoustic force spectroscopy (AFS) (Matsubara et al., 2023) and fitted with an exponential function (Methods).

Figure 1—source data 1

Excel file containing force–extension measurements of the nanospring.

https://cdn.elifesciences.org/articles/108477/elife-108477-fig1-data1-v1.xlsx

At equilibrium, in the absence of tension applied by KIF1A, the NS, uniformly labeled with 124 Cy3 fluorophores, exhibited a resting length of approximately 300 nm and appeared as a diffraction-limited circular spot in a fluorescence image (Figure 1c). In the presence of ATP, the NS extended as the KIF1A moved against the load, and its fluorescence image became an elongated shape (Figure 1c). To calibrate the pulling force acting on the NS from its extension, the AFS was used as in the previous study (Matsubara et al., 2023). Figure 1d presents the relationship between NS extension and force measured by the AFS, revealing its nonlinear elastic behavior. Its nonlinear force–extension relationship was fitted with an exponential function (Methods).

Simulation-based evaluation of fluorescence image analysis methods

Before analyzing the fluorescence images of an NS, we generated artificial images that mimic experimental conditions in order to develop a method for estimating the length of a rod-shaped object from its fluorescence image. The artificial images were generated by placing 116 two-dimensional isotropic Gaussian functions aligned along a straight line (Figure 2a), to mimic a DNA calibration rod labeled with 116 fluorescent molecules, as introduced in the following section. The value of the variance in the Gaussian function was determined based on the approximate shape of the point spread function (PSF) of the Cy3 fluorophore (Figure 2b). While the PSF is, in principle, influenced by the optical system, here it was used as a reference value for generating the artificial images. White noise was then added to the images. When the 116 Gaussian functions were superimposed, they appeared as a single bright spot with an elliptical shape (Figure 2c, top). For this artificial image, we present a comparison of the results obtained using two methods: the Gaussian fitting method and the chain fitting method.

Length estimation using DNA calibration rods.

(a) The fluorescence intensity of the DNA calibration rod was theoretically modeled as a superposition of Gaussian functions with variance σ2 aligned along a straight line, and an artificial image was generated accordingly. (b) The σ value of the Gaussian function used in the model was fitted to match the point spread function of Cy3 fluorescent dye (black line), resulting in a value of 90.65 nm. (c) Artificial image generated by placing 116 two-dimensional Gaussian functions (top). When 116 Gaussian functions were superimposed, they appeared as a single bright spot with an elliptical shape. The image (top) was fitted by the Gaussian fitting method (Equation 1) (middle) and the chain fitting method (Equation 2) (bottom). (d) Estimated length plotted against the true length of the model (L) using the Gaussian fitting method (Equation 1) (dark colors) and the chain fitting method (Equation 2) (bright colors), respectively. The dotted line represents the linear equation y=x. The chain fitting model provides estimates of the true value of L in the simulation. Note that we generated 30 artificial videos, each consisting of 300 frames. (e) Fluorescence micrographs of the DNA calibration rods (Matsubara et al., 2023) obtained in real experiments, with lengths of 398, 501, 599, and 658 nm. The scale bars indicate 2 μm. (f) Estimated length plotted against the true length of the rods using the Gaussian fitting method (Equation 1) (dark colors) and the chain fitting method (Equation 2) (bright colors), respectively. The dotted line represents the linear equation y=x. The chain fitting method provides estimates closer to the true value of L. Note that approximately 30 videos were recorded, each containing about 300 frames for each DNA calibration rod.

Figure 2—source data 1

Excel file containing estimated lengths of fluorescence images of DNA rods for simulations and experiments.

https://cdn.elifesciences.org/articles/108477/elife-108477-fig2-data1-v1.xlsx

In the Gaussian fitting method, the whole fluorescence spot is approximated using a single two-dimensional Gaussian function Ig(x,y), and the estimated standard deviation along the long axis (σlong) (Figure 2c, middle).

(1) Ig(x,y)=Cg exp[(xx0)22σlong+(yy0)22σshort].

Here, Cg is a constant, and x and y represent the major and minor axis directions, respectively. Figure 2d (gray symbols) shows the relationship between the true length L of the rod-shaped object and σlong. If the linear relationship between L and σlong is determined in advance, the value of L can be estimated from the measured σlong (Matsubara et al., 2023).

While the Gaussian fitting method provided reasonable results, we considered an alternative approach, called the chain fitting method, to directly estimate the length L of the model. In the Gaussian fitting method, proposed in the previous study (Matsubara et al., 2023), the true length L of the rod is not included as a parameter in the fitting function (Equation 1). The chain fitting method is designed based on the structure in which fluorescent molecules are arranged in a line, and because the fitting parameters directly include the length L, it is expected to provide a more accurate estimation than the Gaussian fitting method. The fluorescence intensity of this model can be described by the following equation:

(2) Ic(x,y)=CcL2L2exp((xx)2+y22σ2)dx=π2Cc σexp(y22σ2)[Φ(x+L22σ)Φ(xL22σ)],

where Cc is a constant, and Φ is the error function. Figure 2c, bottom, shows the resulting image obtained by the chain fitting method. The simulation results indeed demonstrated that the estimated length by the chain fitting method corresponded with the true length L (Figure 2d, color symbols).

Experimental validation of the chain fitting method using DNA calibration rods

We experimentally validated our approach of the chain fitting method using DNA calibration rods with known lengths (Matsubara et al., 2023). Figure 2e shows the fluorescence images of biotinylated DNA calibration rods immobilized on a glass surface via biotin–BSA and streptavidin (Methods), recorded under the same conditions used for observation of the NSs. As illustrated in the schematic, a total of 116 Cy3 molecules were conjugated to the DNA calibration rod at equal intervals (Matsubara et al., 2023). The length was varied by adjusting the interval between fluorescent dyes. The end-to-end distances between fluorescent molecules were designed to be 398, 501, 599, and 658 nm, respectively. Since the NSs were labeled with 124 Cy3 fluorescent molecules, we considered the fluorescence images obtained from these rods to be similar to those of the NSs.

Figure 2f shows the σlong and length L estimated from experimental fluorescence images of DNA calibration rods using the Gaussian fitting method (Equation 1) and chain fitting method (Equation 2), respectively. There is a proportional relationship between σlong and the true length. On the other hand, L estimated by Equation 2 provides a direct estimate close to the true length of the DNA calibration rod. However, the fact that L was approximately 10% shorter than the true length is considered to be likely due to slight bending of the DNA-origami structures caused by thermal fluctuations. Indeed, in the case of longer rods shown in Figure 2e, a slight bending was observed. Accordingly, the length was corrected by a factor of 1.07 (Methods). In the following, the extension of the NS was primarily estimated using the chain fitting method. Note that results obtained using the conventional Gaussian fitting method are also presented in Figure 3—figure supplement 2. The results from both methods were consistent.

Observation of stall behavior in wild-type KIF1A

Video 1 shows a representative example of the movement of an NS with KIF1A and KIF5B. At first, an inert KIF5B is firmly fixed to the microtubule. Then, the diffusive behavior of a KIF1A—attached to the NS, which is anchored to the microtubule via the inert KIF5B—as it searches for a binding site on the microtubule, is observed. Once the KIF1A binds to the microtubule, it starts moving along the microtubule and stretches the spring. When the tension in the NS balances with the maximum force generated by the KIF1A, the extension temporarily halts during the time interval (attachment duration). Eventually, the KIF1A detaches from the microtubule, the tension is released, and the NS returns to its retracted state. This cycle event was repeatedly observed.

Video 1
Wild-type KIF1A pulling a nanospring.

The video is shown at 10× speed.

The extension of the NS was estimated by the chain fitting method (Equation 2; Figure 3a, bottom). For the event during a time interval (2.8st21.2 s), the histogram of the spring length was calculated (Figure 3a, left). The histogram exhibits a bimodal Gaussian distribution, with each peak corresponding to the extended and retracted states of the NS shown in the micrographs (Figure 3a left). Lstall was calculated by averaging the NS extension during the force-plateau periods within each attachment event. The attachment of KIF1A to the microtubule was identified by the marked suppression of angular fluctuations of the NS relative to the microtubule (Methods and Figure 3a, top), which clearly distinguished the bound state from the diffusive search behavior. Within this attachment duration, the force plateau (stall duration t) was detected by analyzing the rate of relative increase in NS’s length (Methods). The mean NS extension over these plateau segments was defined as Lstall.

Figure 3 with 3 supplements see all
Stall force measurement of wild-type KIF1A homodimers using NSs.

(a) Time course of NS extension (L(t)) in the case of wild-type KIF1A. The black line (trace) represents the average over 10 frames. As illustrated in the schematic in Figure 1c, the NS is stretched (micrograph, top) as a KIF1A moves toward the plus end of the microtubule. When the load reaches the maximum force that the KIF1A can generate, a stall is observed, followed by the detachment of the KIF1A from the microtubule. The NS then returns to its original retracted state (micrograph, bottom). The stall duration t (violet region) was defined based on the angular fluctuation and the rate of relative increase in NS’s length (Methods), where the red regions represent the attachment durations decided based on the angular fluctuations. For each stall event, the histogram of NS extension exhibits a bimodal Gaussian distribution, with the higher peak corresponding to the stall length Lstall (left panel). (b) Histogram of Lstall values calculated from 86 stall events (left). The right panel shows the distribution of stall forces for each KIF1A molecule in which six or more stall events were observed. (c) Lstall as a function of stall duration t with (n = 105) and without (n = 86) PEG. The correlation coefficient (R) between Lstall and t is shown in the figure (left). The right panel presents the magnified view of the blue rectangle in the left panel and clearly indicates that R is small. (d) Comparison of Lstall and t with and without PEG (Mann–Whitney U test, p = 0.4718 for Lstall, p = 0.0616 for t). n.s., not significant (p ≥ 0.05). The green bars indicate the median values along with the first and third quartiles.

Figure 3—source code 1

Time-course analysis of single-molecule videos.

https://cdn.elifesciences.org/articles/108477/elife-108477-fig3-code1-v1.zip
Figure 3—source data 1

Excel file containing the time course, Lstall and t from the stall force experiment of wild-type KIF1A.

https://cdn.elifesciences.org/articles/108477/elife-108477-fig3-data1-v1.xlsx

In our experiments, we analyzed a total of 86 events from 24 KIF1A molecules (Figure 3b, left). From the mean value of Lstall (562±48nm), the stall force of the wild-type KIF1A was estimated to be 4.7 pN (Table 1) by using the force–extension calibration of the NS (Figure 1d). This value is similar to the stall force values of the wild-type KIF1A(1–393) measured previously (Pyrpassopoulos et al., 2023), and larger than the detachment force values (Budaitis et al., 2021; Lam et al., 2021). The right panel of Figure 3b shows the distribution of stall forces for each KIF1A molecule in which six or more stall events were observed. The broad distribution of stall forces arises not only from differences between molecules, but also from the variability of stall forces within a single molecule. The broad distribution is caused by the presence of step-back events (e.g., a green arrow in Figure 3a); however, the chemical state underlying these step backs in KIF1A remains unknown.

Table 1
Stall force of KIF1A(1–393) estimated from the average value of Lstall.

The stall force values were calculated from the mean values of Lstall by using the force–extension relationship of the NS (Figure 1d). The error of Lstall value represents the standard deviation (SD).

Stall force (pN)Detachment force (pN)Termination force (pN)
WT/WT4.7 (n = 86)
(Lstall=562±48 nm)
2.7*
2.2
4
6 §
P305L/P305L0.2 (n = 58)
(Lstall=375±49 nm)
0.7-
P305L/WT0.3 (n = 41)
(Lstall=402±54 nm)
--
V8M/V8M2.2 (n = 95)
(Lstall=516±57 nm)
1.9*-
V8M/WT1.0 (n = 61)
(Lstall=469±27 nm)
--
A255V/A255V3.0 (n = 83)
(Lstall=535±46 nm)
--
A255V/WT2.5 (n = 40)
(Lstall=524±26 nm)
--
  1. Detachment force and termination force measured by using optical tweezers.

  2. *

    Reported in reference (Budaitis et al., 2021).

  3. Reported in reference (Lam et al., 2021).

  4. Reported in reference (Pyrpassopoulos et al., 2023) for the single-bead assay.

  5. §

    Reported in reference (Pyrpassopoulos et al., 2023) for the three-bead assay.

In the results presented here, the extension of the NS was estimated using the chain fitting method (Equation 2). For comparison, we also present the results obtained using the conventional Gaussian fitting method (Equation 1) in Figure 3—figure supplement 2. The two methods provided consistent values, suggesting that future analyses can be performed without the need for prior calibration using DNA rods, which is required by the Gaussian fitting method. In addition, while direct velocity measurements were difficult in the NS-based experiments, we analyzed the behavior under low-load conditions by overlaying the time courses of 10 stall events (Figure 4—figure supplement 1).

Stall duration Δt

In previous studies (Budaitis et al., 2021; Lam et al., 2021), most of the stalls lasted less than 1 s, and the force trajectories exhibited sawtooth-like patterns, which was a characteristic feature of force generation by KIF1A unlike kinesin-1, which exhibits clear stalling behavior (Kojima et al., 1997). However, in our experiments, we observed clear stall states lasting for several tens of seconds, similar to the stalling behavior observed with kinesin-1 (Figure 3a). To quantify this observation, we measured the stall durations (t) and analyzed the relationship between Lstall and t (Figure 3c).

In a single-bead optical tweezers experiment, it was reported that t was correlated with the detachment force when the detachment of a KIF1A from the microtubule was likely accelerated by forces applied perpendicular to the microtubule’s long axis (z-axis depicted in the schematics of Figure 1c), due to the contact between a single bead and an underlying microtubule (Pyrpassopoulos et al., 2020). By using the three-bead optical tweezers system (Pyrpassopoulos et al., 2020; Pyrpassopoulos et al., 2023), which was a modified version of the conventional optical trapping assay, to minimize the z-component of the force, it was found that t did not depend on the detachment force (schematic in Figure 3c; Pyrpassopoulos et al., 2020). Our experimental results also indicated that t was independent of Lstall, as the correlation coefficient was 0.2 (Figure 3c), suggesting that the z-component of the force was effectively minimized.

Effect of electrostatic repulsion on force measurement

Because the DNA-based NS is negatively charged, we investigated the effect of its electrostatic repulsion with negatively charged microtubules. This repulsion could potentially promote the detachment of KIF1A from the microtubule, raising a concern for the measurement. To minimize repulsive interactions between the NS and the microtubule, we coated the NS with oligolysine conjugated to polyethylene glycol (PEG) to render the spring electrically neutral (Methods). Based on a statistical test, there was little difference in Lstall and stall duration (t) between the conditions with and without PEG (Figure 3c, d; Mann–Whitney U test, p = 0.4718 for Lstall, p = 0.0616 for t). These results suggest that the negative charge of the NS had little influence on the stall force measurements.

Stall force measurement of the KIF1A disease-associated variant P305L

P305 of the motor domain of KIF1A is located adjacent to the positively charged K-loop insertion region of KIF1A (Figure 4a). The positively charged K-loop interacts with the negatively charged C-terminal loop of tubulin in microtubules, promoting a high binding affinity to microtubules (Tomishige et al., 2002; Soppina and Verhey, 2014). Therefore, a mutation at P305 directly affects the interaction between KIF1A and microtubules. In other words, the P305L mutation could impair the function of the K-loop, leading to a reduced binding affinity to microtubules (Lam et al., 2021; Anazawa et al., 2022). It was reported that the velocity and run length were reduced by half, and the force decreased fourfold, indicating that the P305L mutation impaired the motility of the motor (Lam et al., 2021). The decrease in force generation in the P305L mutant is a serious defect.

Figure 4 with 2 supplements see all
Stall force measurement of KAND mutants KIF1A homodimers and heterodimers using NSs.

(a) Schematic of KIF1A domain structure showing the functions affected by the P305L, V8M, and A255V mutations (Budaitis et al., 2021). Representative traces of NS extension are shown for homodimers and WT-mutant heterodimers of P305L, V8M, and A255V (b–g). The black lines (traces) represent the average over 10 frames. The lighter-colored regions in the graph represent the attachment durations, while the darker-colored regions indicate the stall durations. The identification of these durations is described in the Methods section.

Figure 4—source data 1

Excel file containing the time course data from the stall force experiment of KIF1A mutants.

https://cdn.elifesciences.org/articles/108477/elife-108477-fig4-data1-v1.xlsx

The wild-type KIF1A depicted in Figure 1b, c was replaced with the P305L mutant, and force measurements were performed using an NS for the mutant (Figure 4b), as shown in Video 2. Although it has been reported that the landing rate of P305L on microtubules was reduced by 97% (Lam et al., 2021; Anazawa et al., 2022), in this experiment, an inert KIF5B firmly anchored to the microtubule at the end of the NS, which kept the P305L mutant in close proximity to the microtubule, thereby effectively increasing its binding probability. In other words, the anchoring effect enabled clear observation of microtubule binding by the P305L mutant. This high binding rate is one of the advantages of the NS–KIF1A–inert KIF5B complex, when used as a force measurement system. Based on the value of Lstall, the force generated by the mutant was reduced to be 0.2 pN (Table 1). Note that the histogram of Lstall is shown in Figure 4—figure supplement 2.

Video 2
The P305L mutant pulling a nanospring.

The video is shown at 1×speed.

KIF1A is considered to primarily form dimers to anterogradely transport synaptic cargos in neurons. In recent years, the motility of heterodimers composed of wild-type and KAND-associated mutant subunits has begun to be investigated experimentally (Budaitis et al., 2021; Anazawa et al., 2022; Kita et al., 2023). In heterozygous patients, WT/WT, P305L/WT, and P305L/P305L dimers coexist, highlighting the importance of studying the heterodimers. As in the previous study (Anazawa et al., 2022), we constructed heterodimers composed of wild-type and P305L subunits (Methods) and performed stall force experiments using these heterodimers (Figure 4c, Table 1).

Stall force measurement of the KIF1A disease-associated variants V8M and A255V

Like kinesin-1, KIF1A also utilizes conformational changes in the neck linker to generate force (Budaitis et al., 2021). The V8M mutation may hinder the neck linker from accessing its docking pocket, which in turn disrupts docking and impairs force generation (Figure 4a). Using the NS-based stall force measurement, we also observed that the stall force of the V8M homodimer was reduced by approximately half (Figure 4d, Table 1), consistent with the previous study (Budaitis et al., 2021).

The A255 residue is located in loop L11 of the switch II cluster and may affect KIF1A function by altering the structure of L11 within switch II or the back door (Lee et al., 2015; Figure 4a). In the single-molecule experiments, although its force has not yet been measured, a decrease in velocity has been reported (Chiba et al., 2019; Guedes-Dias et al., 2019). Loss of function in the switch II region or the back door may impair the conversion of ATP hydrolysis-driven conformational changes into mechanical force. Indeed, our NS-based stall force measurements revealed that the stall force was reduced by approximately half compared to the wild-type (Figure 4f, Table 1). We also performed stall force measurements for heterodimers composed of the KAND mutants V8M (Figure 4e, Table 1) or A255V (Figure 4g, Table 1) and the wild-type subunit. Note that the histogram of Lstall for the mutants in this section is shown in Figure 4—figure supplement 2.

Lstall values for different KIF1A variants

All Lstall values for the KIF1A mutants are summarized in Figure 5a. The values increased in the order of P305L, V8M, and A255V. We then converted the Lstall values into forces and compared them with the clinical severity of KAND (Figure 5—figure supplement 1). Specifically, we compared them with REVEL, CADD v1.4, and the evolutionary scale model (ESM) score reported in the reference (Boyle et al., 2021). Among these metrics, the ensemble prediction score REVEL, which integrates multiple pathogenicity prediction algorithms, showed the correlation. It is important to expand the number of variants analyzed to further investigate the relationship between stall force and the severity of KAND.

Figure 5 with 1 supplement see all
Comparison of Lstall and t among KAND mutants.

Lstall (a) and Δt (b) for homodimers and heterodimers of P305L, V8M, and A255V, compared with WT.

Figure 5—source data 1

Excel file containing Lstall and Δt from the stall force experiment of KIF1A mutants.

https://cdn.elifesciences.org/articles/108477/elife-108477-fig5-data1-v1.xlsx

t values for different KIF1A variants

As shown in Figure 5b, stall duration (t) was compared among KAND mutants (P305L, V8M, and A255V). Because stall behavior has rarely been observed for KIF1A in single-bead optical tweezers experiments, comparison of t has been limited in previous studies (Budaitis et al., 2021; Lam et al., 2021). t did not show as clear differences among the mutants as Lstall did. For WT, which provided a sufficient number of stall events, we investigated the cumulative distribution of t (Figure 3—figure supplement 3). Understanding KIF1A’s behavior associated with these time constants will be important in future studies.

Discussion

We measured the force generated by KIF1A using the NSs developed by Iwaki et al., 2016; Matsubara et al., 2023; Fujita et al., 2019; Figure 1. Because KIF1A more readily detaches from microtubules compared to kinesin-1, and the z-component (refer to Figure 1c for the z-direction) of the applied forces often causes detachment in optical tweezers experiments, making it difficult to observe its stall behavior. In contrast, by using the NS, we were able to apply load to a single KIF1A molecule solely in the horizontal direction (x-axis) along microtubules, enabling clear observation of the stall events (Figure 3). Note that, while performing these measurements, we also improved the image analysis method (Figure 2). In a previous study (Matsubara et al., 2023), the length of the NS was estimated using the Gaussian fitting method (Equation 1). In contrast, in the present study, we adopted the chain fitting method (Equation 2), which allowed for a more direct estimation of the length L from the images. Below, we discuss the advantages and future prospects of NS-based force measurements.

In KAND mutants, ESM scores have been compared with clinical phenotype scores such as the Vineland Adaptive Behavior Composite (VABS-ABC) scores (Sudnawa et al., 2024). ESM refers to large-scale protein language models based on Transformer architectures, trained on millions of sequences to predict structure, function, and mutation effects from primary amino acid sequences (Rives et al., 2021). However, due to the currently low correlation between ESM scores and clinical severity, incorporating biophysical parameters such as velocity and force may provide additional predictive insight, particularly for de novo mutations—a direction already explored in part by a previous study (Rao et al., 2025). In the present study, we propose a high-throughput method for accurately measuring the force generated by KIF1A mutants using the NSs. Anchoring a KIF1A–NS complex to the microtubule via an inert KIF5B facilitated the observation of stall events even for KAND-related KIF1A mutants (Figure 4), which were otherwise prone to detachment due to the dysfunctions caused by the mutations. Indeed, we successfully measured the stall force of several KAND-associated mutants, including P305L, V8M, and A255V. P305L, V8M, and A255V are representative mutants that primarily impair microtubule binding, force generation near the neck linker, and ATP hydrolysis near the Switch II region, respectively (Figure 4a). Notably, stall force measurements for the homodimeric A255V and heterodimeric combinations P305L/WT, V8M/WT, and A255V/WT had not previously been reported using optical tweezers. Building a database of such biophysical measurements will be essential for linking molecular dysfunction to clinical outcomes. The NS-based force assays offer a powerful platform for this purpose. In Figure 5—figure supplement 1, we additionally compared these three mutants with the REVEL, CADD v1.4, and ESM scores (Boyle et al., 2021). It is important to further investigate the relationship between stall force and the severity of KAND.

As noted in the introduction, optical tweezers have been the primary method for measuring forces generated by motor proteins (Schnitzer et al., 2000; Gennerich et al., 2007; Elshenawy et al., 2019; Mallik et al., 2004; Tomishige et al., 2002; Svoboda et al., 1993; Finer et al., 1994; Kaya and Higuchi, 2010; Ariga et al., 2018). An advantage of the NS is that it enables precise investigation of the force exerted by a single motor protein, whereas multiple kinesin molecules may attach to a bead in optical tweezers experiments. Moreover, unlike optical tweezers, which require individual manipulation of each bead and can measure the force of only one motor protein at a time, the NS-based assay allows simultaneous recording and force measurements of multiple NS–KIF1A complexes within the field of view. This results in significantly higher throughput and represents a key advantage over optical tweezers. From the perspective of compatibility with other measurement techniques, the near-infrared laser used in optical tweezers readily causes photobleaching of organic dyes, making single-molecule imaging of dye-labeled proteins challenging. For example, this means that optical tweezers cannot be used simultaneously with MINFLUX (Deguchi et al., 2023). Moreover, the NS can also be combined with high-speed AFM or electron microscopy, enabling structural analysis of proteins under quantitatively controlled mechanical load. In fact, although not NSs, DNA origami-based thick filaments with attached myosin have been investigated using AFM (Fujita et al., 2019). Therefore, the DNA origami-based NSs are also expected to be used with AFM.

Although NSs offer several advantages, they also have disadvantages. While the optical trap can be approximated as a linear spring, the NS behaves as a nonlinear spring (Figure 1d). Consequently, fluctuations in its length can sometimes lead to large errors in force measurements. Therefore, it is crucial to adjust the spring stiffness to match the force range generated by the motor protein. In this study, we used NSs with different physical properties from those used in the previous force measurements on myosin VI (Schnitzer et al., 2000). Wild-type KIF1A generates a larger stall force than myosin, requiring a different dynamic range for the nonlinear behavior of the spring. Owing to DNA origami technology, the mechanical properties of the NS can be tuned to match the appropriate force range for measuring molecular stiffness.

The reason why the stall force exhibits a wide range of values may be not only the non-linearity of the NS, but also the complexity of the underlying kinetics of a single KIF1A, such as step-back events (e.g., the green arrow in Figure 3a). Indeed, analysis of the t distribution (Figure 3—figure supplement 3) suggests that the motion during stall events can be characterized by two (or more) reaction rate constants (Figure 3—figure supplement 3). Determining the detailed chemical reaction model underlying these stall events will be an important subject for future studies.

With the recent development of the super-resolution microscope MINFLUX, live-cell tracking of single fluorescent molecules has become possible, enabling the precise observation of kinesin-1 stepping motion in living cells with nanometer spatial and millisecond temporal resolution (Deguchi et al., 2023). To date, models of chemo-mechanical energy conversion in motor proteins have been proposed based on the load dependence of their stepping motion (Sasaki et al., 2018). The intracellular environment is highly crowded with biomolecules and also contains active fluctuations (Nishizawa et al., 2017), resulting in complex viscoelastic properties and an inhomogeneous refractive index. This leads to challenges in using optical tweezers, such as accurately estimating the trap stiffness of laser-trapped vesicles in cells. On the other hand, the NSs offer a method of measuring force solely through fluorescence imaging, providing a potential solution to this problem. In relation to cellular applications, the NSs have already been used as force sensors to visualize the dynamics of integrins (Matsubara et al., 2023). Expanding such applications, the introduction of the NS–kinesin complex into living cells offers a promising approach to uncover the load-dependent dynamics of kinesin in its native environment.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Recombinant DNA reagentPlasmid pSN672: pET21a(+) human KIF1A(1–393)-LZ-6HisAnazawa et al., 2022RRID:Addgene_177362Template for SNAP-tag insertion
Recombinant DNA reagentPlasmid pSN643: pET28a(+) human KIF1A(1–393)-LZ-mScarlet-I-Strep-tagIIAnazawa et al., 2022Template for SNAP-tag insertion
Recombinant DNA reagentPlasmid pSNAP-tag (T7)-2 vectorNew England Biolabs#N9181STemplate for SNAP-tag cDNA amplification
Recombinant DNA reagentPlasmid pET21a(+) human KIF1A(1–393)-LZ-SNAP-6HisThis studyMotor domain construct for force measurement
Recombinant DNA reagentPlasmid pET28a(+) human KIF1A(1–393)-LZ-SNAP-Strep-tagIIThis studyUsed for heterodimer preparation
Recombinant DNA reagentPlasmid pET21a(+) human KIF1A(1–393, P305L)-LZ-SNAP-6HisThis studyKIF1A mutant construct associated with KAND individuals
Recombinant DNA reagentPlasmid pET21a(+) human KIF1A(1–393, V8M)-LZ-SNAP-6HisThis studyKIF1A mutant construct associated with KAND individuals
Recombinant DNA reagentPlasmid pET21a(+) human KIF1A(1–393, A255V)-LZ-SNAP-6HisThis studyKIF1A mutant construct associated with KAND individuals
Recombinant DNA reagentPlasmid pYS05: KIF5B(1–560,G234A)::RA::SNAP::6HisShimamoto et al., 2015Inactive kinesin-5B mutant; gift from Dr. Yuta Shimamoto
Strain, strain background (Escherichia coli)BL21(DE3)Novagen#69450
Peptide, recombinant proteinTubulin (porcine brain)Tokyo Shibaura OrganPurified in-house
Peptide, recombinant proteinNeutrAvidinThermo Scientific31000DNA calibration rod immobilization
Peptide, recombinant proteinBiotinylated BSAThermo Scientific29130Surface coating
Chemical compound, drugIPTGSigma-AldrichProtein expression inducer
Chemical compound, drugTaxolFUJIFILM Wako163–28163Microtubule stabilization
Chemical compound, drugCaseinSigma-AldrichC5890-500GSurface blocking
Chemical compound, drugBG-GLA-NHSNew England BiolabsS9151SSNAP-tag labeling
Chemical compound, drugd-DesthiobiotinSigma-AldrichD1411Strep-tag elution
Chemical compound, drugK10–PEG5KAlamanda PolymersmPEG5K-b-PLKC10DNA origami surface passivation
Sequence-based reagentp8064
(DNA scaffold)
Tilibit NanosystemsDNA origami scaffold
Sequence-based reagentCore and handle staples
(oligonucleotides)
Integrated DNA TechnologiesSequences in Tables S2–S7
Sequence-based reagentDNA nanospring
(DNA nanostructure)
This studyDesigned with caDNAno
Sequence-based reagentDNA calibration rod
(DNA nanostructure)
Matsubara et al., 2023Used for force calibration
Software, algorithmcaDNAnoDouglas et al., 2009DNA origami design
Software, algorithmHigh Speed RecordingHamamatsu PhotonicsVideo acquisition
Software, algorithmPythonPython Software FoundationRRID:SCR_008394Version 3.11.3
Software, algorithmChain fitting analysisThis studyCustom Python code using scipy.optimize

Plasmid constructs

Request a detailed protocol

The plasmids used in this study are listed in Supplementary file 1. A recombinant human KIF1A(1–393)-LZ-SNAP-6His construct containing the motor domain (amino acids 1–361), neck linker, neck coil domain, and a GCN4 leucine zipper (LZ) for stabilized dimerization was used for NS force generated study (Figure 1a). To generate this construct, a human KIF1A(1–393)-LZ-6His construct named pSN672 (Addgene #177362) (Anazawa et al., 2022) was modified as follows: SNAP-tag cDNA was amplified by polymerase chain reaction (PCR) with high-fidelity PCR enzyme (Prime STAR MAX, Takara Bio Inc, Kusatsu, Japan, Cat#R045), using pSNAP-tag(T7)-2 vector (New England Biolabs, N9181S) as a template. Then the PCR product added with XhoI sites at both ends was inserted into the XhoI site of the plasmid pSN672 by using In-Fusion HD Cloning Kit (Takara Bio Inc, Cat#639649). This KIF1A (1–393)-LZ-SNAP-6His construct was then used as a template to generate the three aberrant motors by site-directed mutagenesis using the KOD Plus mutagenesis kit (TOYOBO, Cat#SMK-101). A recombinant human KIF1A(1–393)-LZ-SNAP-Strep-tagII construct was used to develop heterodimers of wild-type and mutant KIF1A. To create this KIF1A motor domain construct, a human KIF1A(1–393)-LZ-mScarlet-I-Strep-tagII construct, namely pSN643, was modified as follows: The mScarlet-I sequence of pSN643 (Anazawa et al., 2022) was replaced by SNAP cDNA sequence. All constructs were confirmed by sequencing. Inactive mutant kinesin 5 B (G234A) expression vector was a gift from Dr. Yuta Shimamoto (National Institute of Genetics), which has been described (pYS05 K560 pET-17b KIF5B (1–560, G234A)::LZ::His-tag) (Rice et al., 1999; Shimamoto et al., 2015).

Recombinant protein expression and purification (homodimers)

Request a detailed protocol

Proteins were expressed in BL21(DE3) and SNAP-tagged KIF1A motor domain homodimers and inactive mutant KIF5B (G234A) homodimer were purified as previously described (Anazawa et al., 2022). Briefly, His-tagged SNAP-tagged KIF1A motor domain plasmids were transformed into BL21(DE3) (Novagen #69450), and the cells were cultured on LB agar supplemented with ampicillin at 37°C overnight. Colonies were picked and cultured in 10 ml LB medium supplemented with ampicillin overnight. Next morning, the culture was transferred to 750 ml of 2.5× YT (20 g/l Tryptone, 12.5 g/l Yeast Extract, 6.5 g/l NaCl) supplemented with 10 mM phosphate buffer (pH 7.4) and 100 μg/ml ampicillin in 2 l flasks and shaken at 37°C. Two flasks were routinely prepared. When OD600 reached 0.6, the flasks were cooled in ice-cold water for 30 min, then isopropyl-β-D-thiogalactoside was added to the cooled culture to final concentration of 0.2 mM to induce expression. The culture was shaken at 18°C overnight (~16 hr), and the bacteria expressing recombinant proteins were harvested by centrifugation at 4000 rpm, 20 min, 4°C (HITACHI himac CR-21). The cell pellet was resuspended in phosphate-buffered saline; PBS(−) and centrifuged again (4000 rpm, 20 min, 4°C). Supernatant was discarded completely, then cell pellets were frozen and stored at –80°C until use. To purify recombinant protein, bacteria pellets were resuspended in protein buffer (50 mM HEPES, pH 8.0, 150 mM KCH3COO, 2 mM MgSO4, 10% glycerol), supplemented with 1 mM ATP, protease inhibitors, such as 4-(2-Aminoethyl) benzenesulfonylfluoride (FUJIFILM), Leupeptin (FUJIFILM, Wako #336-40413), Aprotinin from Bovine Lung (FUJIFILM, Wako #013-28311), and Pepstatin A (FUJIFILM, Wako #330-43973), using 5 ml protein buffer per gram of wet cell paste. After gentle vortex, bacteria resuspension was with Benzonase nuclease (Novagen) and Lysozyme (Merch), followed by sonication using Branson SONIFIER 250. Lysate was obtained by centrifugation (40,000 rpm, 30 min, 4°C, HITACHI himac CP-80α). After centrifugation, the supernatant was loaded on TALON Metal Affinity Resin (Takara Bio Inc, Cat#635502), and incubated for 30 min at 4°C with mild shaking under platform shaker. The resin was washed twice with His-tag wash buffer (50 mM HEPES, pH 8.0, 450 mM KCH3COO, 2 mM MgSO4, 10% glycerol, 10 mM imidazole), containing 100 μM ATP. Protein was eluted with 10 ml His-tag elution buffer (50 mM HEPES, pH 8.0, 150 mM KCH3COO, 2 mM MgSO4, 10% glycerol, 500 mM imidazole) containing 100 μM ATP. Eluted fractions were collected and concentrated to ~300 μl using Amicon Ultra centrifugal filters (Merck, Darmstadt, Germany). The affinity-purified protein was further separated by gel filtration using a Superdex 200 Increase 10/300 GL column (Cytiva) in protein buffer on an NGC Chromatography system (Bio-Rad Laboratories). Peak fractions were collected and aliquoted and snap-frozen in liquid nitrogen. Protein concentration was assessed by Bradford method (Bio-Rad Protein assay dye Reagent Concentrates, Cat#5000006). Each step of purification was analyzed by SDS–polyacrylamide gel electrophoresis (SDS–PAGE) (Figure 1—figure supplement 1). All recombinant DNA experiments conducted in accordance with the Cartagena Protocol on Biosafety and all procedures were approved by the Genetic Modification Safety Committee of Yamaguchi University (J18016). Preparation of all plasmids and recombinant proteins was performed at Yamaguchi University Center for Gene Research.

Purification of heterodimers

Request a detailed protocol

BL21(DE3) cells transformed with KIF1A(1–393)::LZ::SNAP:: Strep-tag II plasmid were cultured in LB supplemented with kanamycin at 37°C. Competent cells were prepared using a Mix&Go kit (Zymogen). The competent cells were further transformed with mutant KIF1A(1–393)::LZ::SNAP::His plasmid and selected on LB agar supplemented with ampicillin and kanamycin. Colonies were picked and cultured in 10 ml LB medium supplemented with ampicillin and kanamycin overnight. Next morning, 10 ml of the medium was transferred to 750 ml 2.5× YT supplemented with carbenicillin and kanamycin in a 2-l flask and shaken at 37°C. Two flasks were routinely prepared. The procedures for protein expression in bacteria and preparation of bacterial lysate were the same as for the purification of homodimers except that the component of purification buffer was 50 mM HEPES, pH 8.0, 150 mM KCH3COO, 2 mM MgSO4, 1 mM EGTA, 10% glycerol. Lysate was loaded on Strep-Tactin resin (IBA, GmbH)(bead volume: 2 ml) and incubated for 1 hr at 4°C with mild shaking under platform shaker. The resin was washed with 40 ml wash buffer. Protein was eluted with 10 ml protein buffer supplemented with 2.5 mM d-Desthibiotin (Merck, Sigma-Aldrich, D1411). The eluted solution was then loaded on Ni-NTA agarose (QIAGEN, #30210) (bead volume: 2 ml) and incubated for 1 hr at 4°C with mild shaking under platform shaker. The resin was washed with 40 ml His-tag wash buffer (50 mM HEPES, pH 8.0, 450 mM KCH3COO, 2 mM MgSO4, 10 mM imidazole, 10% glycerol), containing 100 μM ATP. Protein was eluted with His-tag elution buffer (50 mM HEPES, pH 8.0, 150 mM KCH3COO, 2 mM MgSO4, 10% glycerol, 500 mM imidazole) containing 100 μM ATP. Eluted fractions were collected and concentrated to ~300 μl using Amicon Ultra centrifugal filters (Merck), then subjected to separation on NGC chromatography system (Bio-Rad) equipped with a Superdex 200 Increase 10/300 GL column (Cytiva). Peak fractions were pooled, concentrated in an Amicon filter again, and snap frozen in liquid nitrogen. The concentration and quality of the protein were assessed by Bradford method and SDS–PAGE, respectively (Figure 1—figure supplement 2).

Microtubule preparation

Request a detailed protocol

Tubulin was purified from porcine brain (Tokyo Shibaura Organ Co, Ltd) through two cycles of polymerization and depolymerization using the revised method of reference (Castoldi and Popov, 2003), and stored at −80°C in PEM buffer (100 mM PIPES, 1 mM MgCl2, 1 mM EGTA, pH 6.9, adjusted with KOH). For polymerization, the tubulin solution was mixed with 2 mM GTP, 10 mM MgCl2, and 15% DMSO, and was then incubated at 37°C for 30 min. Prior to polymerization, the tubulin solution was centrifuged at 180,000 × g for 5 min at 2°C, and the supernatant was collected. The polymerized solution was then incubated at room temperature for 20 min in BRB80 buffer (80 mM PIPES, 1 mM MgCl2, 1 mM EGTA, pH 6.8 w/ KOH) supplemented with 10 µM Taxol (163-28163; FUJIFILM) to stabilize the microtubules. The microtubules were pelleted by centrifugation and then resuspended in BRB80 buffer. Microtubules were stored in the dark at room temperature for up to 2 weeks.

Preparation of NSs and DNA calibration rod

Request a detailed protocol

NSs were designed using caDNAno software (Douglas et al., 2009), as described in the method from the reference (Matsubara et al., 2023). Briefly, NSs were designed by introducing a negative superhelical strain into a four-helix bundle arranged on a square lattice. To impose the torsional strain, two additional nucleotides were periodically inserted every 32 bp in the upper pair of helices, while two nucleotides were deleted at the corresponding positions in the lower pair of helices. This design yielded a uniform coil structure with a diameter of ~35 nm. For the kinesin–NS conjugation, 32-base single-stranded DNA (ssDNA) handles (handle 32A and 32B in Supplementary file 2) were extended from both ends of the NS. To construct the NS, 10 nM scaffold (p8064, tilibit nanosystems) was mixed with 100 nM core staples (Supplementary file 3), 100 nM handle staples (Supplementary file 2), and 12.4 mM Cy3-labeled antihandle (Supplementary file 4). The folding reaction was carried out in folding buffer (5 mM Tris pH 8.0, 1 mM EDTA, and 14 mM MgCl2) with rapid heating to 80°C and cooling in single degree steps to 60°C over 2 hr followed by additional cooling in single degree steps to 25°C over another 24 hr. DNA calibration rods were also prepared as previously described (Matsubara et al., 2023). Briefly, 10 nM scaffold (p8064) was mixed with 100 nM core staples (Supplementary file 5), 100 nM handle staples (Supplementary file 6), 12.4 mM Cy3-labeled antihandle and 800 nM biotin-labeled antihandle (Supplementary file 7). The folding reaction was carried out in folding buffer (5 mM Tris pH 8.0, 1 mM EDTA, and 14 mM MgCl2), and the heating and cooling condition was the same as that for the NS. The folded DNA nanostructures were purified by glycerol gradient ultracentrifugation as previously described (Matsubara et al., 2023). Briefly, 15–45% (vol/vol) gradient glycerol solutions in buffer A (1× TE buffer containing 11 mM MgCl2) were made, and the glycerol fractions containing monomeric nanostructures were determined by agarose gel electrophoresis. The concentration of the nanostructures was determined with a Nanodrop spectrophotometer (Thermo Scientific), and the solution was aliquoted and stored at −80°C until use. To minimize electrostatic interference between NSs and microtubules, NSs were coated with oligolysine conjugated to PEG (K10 PEG5K, purchased from Alamanda Polymers) (Castoldi and Popov, 2003) as follows. 10 µl of 4.2 nM NSs were mixed with 0.45 µl of K10–PEG5K at a concentration corresponding to an N:P ratio (ratio of nitrogen in amines to phosphates in DNA) of 0.5:1. The mixture was then incubated at room temperature for 30 min.

Kinesin–NS conjugation

Request a detailed protocol

A 32-base amine-modified DNA oligonucleotide (Hokkaido System Science) (oligo 32A* and 32B* in Supplementary file 2) was conjugated with BG-GLA-NHS (S9151S; New England Biolabs) through the reaction between the ester and amine groups (BG-oligonucleotide). Here, the oligo 32A* and 32B* was complementary to the ssDNA handle sequence (handle 32A and 32B) of the NS. Then, the 35 μM BG-oligonucleotide (2 μl) and 1 μM SNAP-tagged KIF1A (50 μl) or 1 μM SNAP-tagged inert KIF5B mutant (G234A) (50 μl) were mixed and incubated at room temperature for 1 hr. The oligo-labeled kinesin was aliquoted and stored at –80°C. When the oligo-labeled kinesin and NS were conjugated, 5 μl of oligo 32A*-labeled KIF1A or oligo 32B*-labeled inert KIF5B (0.5–1 μM) was mixed with 10 μl of NSs (~5 nM) in SRP90 buffer (90 mM HEPES, 50 mM CH3COOK, 2 mM Mg(CH3COO)2, 1 mM EGTA, pH 7.6, adjusted with KOH), and incubated on ice for 30 min.

Single-molecule experiment on kinesin–NS conjugates

Request a detailed protocol

A flow chamber was constructed from two uncoated coverslips of different sizes: an 18 mm × 18 mm coverslip (2918COVER18-18, Iwaki) was placed on top, and a 24 mm × 24 mm coverslip (C024241, Matsunami) was placed at the bottom, separated by two spacers of ~50 µm thickness, resulting in a volume of approximately 10 µl. Unless otherwise noted, the assay buffer was the SRP90 buffer supplemented with 20 µM Taxol, 0.4 mM 2-mercaptoethanol (131-14572, FUJIFILM), and 0.2 mg/ml casein (C5890-500G, Sigma-Aldrich). Microtubules were diluted in the buffer (without casein) and flowed into the chamber, allowing them to attach to the glass surface of the chamber through nonspecific interaction during a 5-min incubation. Then, the chamber was incubated with the blocking buffer (assay buffer containing 0.5 mg/ml casein) for 5 min to prevent nonspecific binding, during which microtubules that did not adhere to the glass surface were washed away using the blocking buffer. Next, the assay buffer containing kinesin–NS conjugates was flowed into the chamber, followed by a 5-min incubation to allow the kinesin–NS conjugates to attach to the microtubules. Finally, the assay buffer, additionally supplemented with 1 mM ATP, 3 mg/ml glucose, 0.1 mg/ml glucose oxidase, and 0.02 mg/ml catalase, was flowed into the chamber. Note that the latter three regents are an oxygen scavenger system for fluorescence antifade protection. This step allowed for the removal of kinesin–NS conjugates that had not attached to the microtubules. Finally, to prevent the evaporation of the buffer inside the chamber during the experiment, the open ends of the flow chamber were sealed with clear nail polish. The motion of kinesin–NS conjugates was observed using a fluorescence microscope (IX83, EVIDENT) equipped with a 100×oil immersion objective lens (UPlanFL N 100×/1.30 Oil, EVIDENT) and a fluorescence mirror unit (TRITC-B, Semrock), and was recorded at 33 frames per second (fps) with an effective pixel size of 65 nm by a CMOS camera (ORCA-Flash4.0, Hamamatsu Photonics). Experiments were also performed at 100 fps in addition to the recording rate of 33 fps (Figure 3—figure supplement 1). The flow chamber was maintained at 25°C using a heating plate (HP-R-Z002, Live Cell Instrument) during observation.

Observation of DNA calibration rod

Request a detailed protocol

A flow chamber was prepared in the same manner as described for the observation of kinesin–NS conjugates. Biotinylated BSA (0.5 mg/ml, 29130, Thermo Scientific) was incubated for 5 min to coat the glass surface of the chamber. Excess biotinylated BSA was washed by the assay buffer (SRP90 buffer supplemented with 10 mM Mg(CH3COO)2). Subsequently, neutravidin (0.5 mg/ml, 31000, Thermo Scientific) was flowed into the chamber and incubated for 3 min. Unbound neutravidin was removed by the assay buffer. Finally, ~6 pM biotinylated DNA calibration rods (Matsubara et al., 2023) were flowed into the chamber and incubated for 5 min to allow their attachment to the glass surface. Unbound DNA calibration rods were washed out by the assay buffer containing the oxygen scavenger system for fluorescence antifade protection (3 mg/ml glucose, 0.1 mg/ml glucose oxidase, and 0.02 mg/ml catalase).

Conversion from image to force

Request a detailed protocol

Video acquisition was conducted using the High Speed Recording software (Hamamatsu Photonics), and the exported 8-bit video data in uncompressed AVI format were imported into a custom-written Python program for further analysis. The fluorescent spots corresponding to the DNA calibration rods and NSs were analyzed using the chain fit method (Equation 2). To estimate their lengths (L) using Equation 2, parameter fitting was performed using curve fitting via the ‘curve_fit’ function from the ‘scipy.optimize’ module of Python (version 3.11.3). Data acquisition of the NS fluorescence images was performed under the same optical conditions (e.g., excitation laser intensity, dichroic mirror, filter, and camera settings) as those used for the DNA calibration rod. A correction was applied to the NS’s length (L) to account for possible bending of the DNA origami structure due to thermal fluctuations (Figure 2e, f). The modified length LmC is calculated as LmC=1.07×L-2.35 where L is the estimated value by the chain fitting method (Equation 2). For the parameter σlong estimated by the Gaussian fitting method (Equation 1), the length LmG is calculated as LmG=3.67×σlong-224. These empirical linear relationships between LmC and L, and between LmG and L were obtained from the DNA calibration rod experiment presented in Figure 2f. Then, the force value corresponding to the length value (LmC or LmG) was determined from the force–extension curve (Figure 1d of the main text). Here, the force–extension curve was fitted with an exponential function (Fx=4.08×10-4exp(16.6x)), where x is the extension in µm and Fx is the force in pN.

Angular fluctuation

Request a detailed protocol

Let θt denotes the angle between the major axes of a microtubule and an NS images (the major axis of the ellipse fitted to the fluorescence image of the NS) in frame t. Then angle θt is converted to ϕt=2θt[π,π], and mapped to a unit vector vt=(cosϕt,sinϕt). The mean of vt over 10 frames is defined as Rt=1ωj=tω/2t+ω/2+1vt. The angular fluctuation is quantified as 1|Rt|. Here, the relation 1|Rt|0 indicates that an NS is aligned with the direction of a microtubule, while 1|Rt| approaching 1 indicates that the NS is oriented away from the direction of the microtubule. We considered the NS to be attached to the microtubule when 1|Rt|<0.1. If this condition lasted for at least 5 consecutive frames, the corresponding time interval was measured as the attachment duration.

Rate of relative increase in NS’s length

Request a detailed protocol

The attachment of KIF1A to the microtubule was identified by the marked suppression of angular fluctuations of the NS relative to the microtubule, which clearly distinguished the bound state from the diffusive search behavior (e.g., red regions in Figure 3a). Within this attachment duration, the force plateau was detected by analyzing the temporal change in NS length (e.g., violet regions in Figure 3a): the 10-frame moving average of the extension, L¯(τ)(=(t=10(τ1)+110τL(t))/10) where t represents a video frame, was computed, and the regions where the relative increase rate (L¯(τ)L¯(τ1))/L¯(τ1)0.02 were regarded as plateau segments. The mean NS extension over these plateau segments was defined as Lstall.

Statistical analysis

Request a detailed protocol

Statistical analyses were performed using Python (version 3.11.3) with scipy.stats (for Shapiro–Wilk, Levene, Mann–Whitney U, and Kruskal–Wallis tests) and scikit-posthocs (for Dunn’s post hoc test with Bonferroni correction). Statistical comparisons were conducted for Lstall and t measurements. The normality of each dataset was assessed using the Shapiro–Wilk test. All datasets were non-normally distributed, except for Lstall values obtained from the homozygous P305L mutant (P305L/P305L) and heterozygous V8M mutant (V8M/WT). Homogeneity of variances was evaluated using Levene’s test. Homogeneity of variance was confirmed for Lstall between the w/o PEG and w/ PEG conditions, and among WT/WT and A255V genotypes. However, variance differed significantly among WT/WT and P305L genotypes and among WT/WT and V8M genotypes. All t datasets showed unequal variances. Based on these results, non-parametric statistical tests were applied. The Mann–Whitney U test was used for comparing Lstall between w/o PEG and w/ PEG, and for comparing t between homozygous and heterozygous mutant combinations. The Kruskal–Wallis test was applied to compare Lstall across WT/WT, mutant/mutant, and mutant/WT combinations for each KIF1A mutant (P305L, V8M, and A255V), followed by Dunn’s test with Bonferroni correction.

Data availability

All data generated or analyzed during this study are included in the manuscript and supporting files; source data files have been provided. The primary dataset has been deposited in Dryad (https://doi.org/10.5061/dryad.6q573n6d4).

The following data sets were generated
    1. Takamatsu N
    2. Furumoto H
    3. Ariga T
    4. Iwaki M
    5. Hayashi K
    (2026) Dryad Digital Repository
    Data from: Stall force measurement of the kinesin-3 motor KIF1A using a programmable DNA origami nanospring.
    https://doi.org/10.5061/dryad.6q573n6d4

References

Article and author information

Author details

  1. Nobumichi Takamatsu

    1. The Institute for Solid State Physics, The University of Tokyo, Chiba, Japan
    2. Department of Complexity Science and Engineering, Graduate School of Frontier Sciences, The University of Tokyo, Chiba, Japan
    Contribution
    Data curation, Software, Formal analysis, Methodology, Writing – original draft
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0009-0001-3048-3467
  2. Hiroko Furumoto

    Systems Biochemistry in Pathology and Regeneration, Yamaguchi University Graduate School of Medicine, Ube, Japan
    Contribution
    Resources, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0009-0002-9587-8686
  3. Takayuki Ariga

    Graduate School of Frontier Bioscience, The University of Osaka, Osaka, Japan
    Contribution
    Conceptualization, Resources, Funding acquisition, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2045-5761
  4. Mitsuhiro Iwaki

    1. Graduate School of Frontier Bioscience, The University of Osaka, Osaka, Japan
    2. Advanced ICT Research Institute, National Institute of Information and Communications Technology, Kobe, Japan
    3. Immunology Frontier Research Center (IFReC), The University of Osaka, Osaka, Japan
    Contribution
    Conceptualization, Resources, Funding acquisition, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    iwakim@nict.go.jp
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8893-1103
  5. Kumiko Hayashi

    1. The Institute for Solid State Physics, The University of Tokyo, Chiba, Japan
    2. Department of Complexity Science and Engineering, Graduate School of Frontier Sciences, The University of Tokyo, Chiba, Japan
    Contribution
    Conceptualization, Supervision, Funding acquisition, Writing – original draft, Project administration
    For correspondence
    hayashi@issp.u-tokyo.ac.jp
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2895-8945

Funding

Japan Society for the Promotion of Science (23H02442)

  • Kumiko Hayashi

Japan Society for the Promotion of Science (21H01053)

  • Mitsuhiro Iwaki

Japan Science and Technology Agency

https://doi.org/10.52926/jpmjcr2023
  • Mitsuhiro Iwaki

Japan Science and Technology Agency

https://doi.org/10.52926/JPMJPR21E2
  • Takayuki Ariga

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

We acknowledge Dr. Yuta Shimamoto for providing plasmids, and Ms. Kimiko Nagino and Mr. Gai Ohashi for purifying microtubules. This work was supported by JSPS KAKENHI (Grant No. 23H02442), the Precise Measurement Technology Promotion Foundation (PMTP-F), and a KIF1A.org Mini Grant to KH, JSPS KAKENHI (Grant No. 21H01053) and JST, CREST (Grant No. JPMJCR2023) to MI, JST PRESTO (Grant No. JPMJPR21E2) to TA.

Version history

  1. Sent for peer review:
  2. Preprint posted:
  3. Reviewed Preprint version 1:
  4. Reviewed Preprint version 2:
  5. Version of Record published:

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.108477. This DOI represents all versions, and will always resolve to the latest one.

Copyright

© 2025, Takamatsu 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

  • 1,292
    views
  • 62
    downloads
  • 0
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Nobumichi Takamatsu
  2. Hiroko Furumoto
  3. Takayuki Ariga
  4. Mitsuhiro Iwaki
  5. Kumiko Hayashi
(2026)
Stall force measurement of the kinesin-3 motor KIF1A using a programmable DNA origami nanospring
eLife 14:RP108477.
https://doi.org/10.7554/eLife.108477.3

Share this article

https://doi.org/10.7554/eLife.108477