Dynamics of cooperative transport by multiple kinesin motors and diffusing microtubule-associated proteins

Yao Wang,Yu-Ru Liu,Peng-Ye Wangand Ping Xie,*

1 Key Laboratory of Soft Matter Physics,Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China

2 School of Physical Sciences,University of Chinese Academy of Sciences,Beijing 100049,China

Abstract In eukaryote cells,cargos are often transported cooperatively by kinesin motors and nonmotor microtubule-associated proteins(MAPs).The prior in vitro experimental data showed that the velocity of the cargo transported by kinesin motors and Ndc80(a member of MAP)proteins of truncated coiled-coil stalks decreases sensitively with the increase of the ratio of Ndc80 to motor number.However,the underlying mechanism of Ndc80 affecting sensitively the cooperative cargo transport by kinesin motors is unclear.To understand the mechanism,here we study numerically the cooperative cargo transport by kinesin motors and Ndc80 proteins.Our results showed that for the case of the motors and Ndc80 proteins with truncated short stalks,as used in the experiments,the calculated results reproduce quantitatively the prior experimental data.The mechanism of the cargo velocity decreasing sensitively with the ratio of Ndc80 to motor number is revealed.By contrast,for the case of the motors and Ndc80 proteins with full-length long stalks,the velocity of the cargo decreases slowly with the increase in the ratio of Ndc80 to kinesin number.Our results thus give an explanation of why the kinesin motors working in the cell have long stalks.

Supplementary material for this article is available online

Keywords:intracellular transport,molecular motor,kinesin,microtubule-associated protein,Ndc80

1.Introduction

Efficient cargo transports in cells are performed by motor proteins such as kinesin motors[1–6].Nonmotor microtubule-associated proteins(MAPs)are also involved in intracellular transports[6–14].Kinesin motors can move directionally along microtubule(MT)filaments towards the plus end by the hydrolysis of ATP[3–5].Many MAPs without ATPase activity can diffuse unbiasedly along the MT filaments[13–15].Ndc80 protein is a member of diffusing MAPs[16,17].It is the key MT-binding component of the kinetochore,participating in the transport of kinetochores[17–21].During prometaphase in cell division,Ndc80 proteins cooperate with kinesin-7 CENP-E motors to transport the kinetochore towards the MT plus end[22–24].Recently,Chakraborty et al[17]studied in detail the cooperative cargo transport by multiple CENP-E motors and Ndc80 proteins,determining quantitatively the effect of the Ndc80 proteins on the dynamics of the cargo transport.For comparison,they also made studies with CENP-E motors being replaced by kinesin-1 motors[17].Puzzlingly,they found experimentally that the presence of Ndc80 proteins has a sensitive effect on the reduction of the velocity of the cargo transport by multiple kinesin motors[17].For example,in the case of kinesin-1 motors,the velocity decreases linearly with the increase of the ratio of the number of Ndc80 proteins to motors,and at the ratio of 0.8,the velocity decreases by about 2.7-fold relative to that in the absence of the Ndc80[17].For the case of CENP-E motors,as the ratio increases the velocity decreases much quicker and at the ratio of only 1,the velocity becomes nearly 0[17].By contrast,the experimental data showed that the velocity of the cooperative transports by multiple kinesin-1 motors is kept nearly unchanged with the motor number[25]and the velocity of the cooperative transports by multiple CENP-E motors at high motor concentrations decreases by about 1.4-fold relative to the velocity of the single motor[26].

While the dynamics of the cooperative cargo transport by multiple molecular motors have been studied intensively and extensively[27–41],the theoretical and numerical study of the dynamics of the cooperative cargo transport by multiple molecular motors and diffusing MAPs has not been paid much attention to.The underlying mechanism of the effect of the diffusing MAPs on the cooperative cargo transport by kinesin motors is unclear.In particular,the puzzling experimental results showing that the presence of Ndc80 proteins has a sensitive effect on the reduction of the cargo velocity have not been explained up to now[17].In addition,the difference in the feature for the motor-number dependence of the velocity of the cooperative transports by multiple kinesin-1 motors from that by multiple CENP-E motors[25,26]has not been explained.

To understand the above-mentioned mechanism and to explain the puzzling experimental data,in this work we analyze the cooperative cargo transport by multiple kinesin motors and Ndc80 proteins.As done in the experiments of Chakraborty et al[17],we also consider two families of kinesin motors—kinesin-1 and kinesin-7 CENP-E.Since the two heads of the motor are connected to the common coiledcoil stalk via the flexible neck linkers(NLs),it is argued that the relatively rigid stalk should rotate freely in a large range of the angle relative to the heads.Thus,we first use all-atom molecular dynamics(MD)simulations to validate the argument.The large free rotation of the stalk indicates that in the cooperative cargo transport by multiple motors a motor can move freely in a range of the distance on MT without experiencing any external force that arises from the stretching of its stalk,as argued before[30–36].With the consideration of the presence of such a range,we show that the velocity of the cargo cooperatively transported by multiple kinesin motors versus the motor number can reproduce well the available experimental data.On the basis of the above studies,we then study the cooperative cargo transport by multiple kinesin motors and Ndc80 proteins with truncated stalks,as done in the experiments of Chakraborty et al[17],explaining quantitatively the experimental data.Furthermore,we study the cooperative cargo transport by multiple kinesin motors and Ndc80 proteins with full-length stalks,showing that the cargo velocity decreases only slowly with the ratio of Ndc80 to motor number.

2.Models for dynamics of single kinesin motors and dynamics of single Ndc80 proteins

2.1.Kinesin motor

Consider an external force,F,on the stalk of the motor,with F resisting the forward movement of the motor being defined to have a positive value.As it is known,during the processive stepping of the motor each step consists of a long dwell period,during which the position of the motor is kept unchanged,and a very short moving period,during which the position of the motor varies.In the single-molecule optical trapping experiments,the motor experiences a constant external force F during the dwell period.During the moving period,here we consider two models for the external force experienced by the motor to determine the stepping rate.

The first model is that during the moving period of one step the motor experiences the constant external force F that is equal to that during the dwell period before the occurrence of the moving period,as generally studied in the literature[27,28,34–41].For the convenience of writing,this model is called the constant-force model here.Based on the model,the forward and backward stepping of the motor as the function of F can be written as[42–44]

where d=8.2 nm is the step size[45,46].

The second model is that during the moving period of one step the motor experiences the varying external force,as studied recently[32,47].This is because even in the single molecule optical trapping experiments with the feedback to ensure a ‘constant’ load on the motor the moving time in one step is much shorter than the response time of the feedback,as discussed before[47–51].For the cooperative cargo transport by multiple motors and MAPs,during the moving period of one motor in one step the other motors and MAPs are unmoved and thus the moving motor would experience a varying force arising from the stretching of the linkers connecting the common cargo and motors and MAPs.This model is called the varying-force model here.Based on the model,the forward and backward stepping rates can be written as[47]

Figure 1.Dynamics of the single kinesin motors in the optical trapping experiments.(a)Velocity versus external force for kinesin-1.The thick black line represents theoretical results obtained with the varying-force model.The thin green line represents theoretical results obtained with the constant-force model.Circles are experimental data at saturating ATP(1 mM)measured by Nishiyama et al[49].(b)Velocity versus external force for kinesin-7 CENP-E.The thick black line represents theoretical results obtained with the varying-force model.The thin green line represents theoretical results obtained with the constant-force model.Circles are experimental data for truncated CENP-E at saturating ATP(2 mM)measured by Gudimchuk et al[46].(c)Detachment time versus external force for kinesin-7 CENP-E.The line represents theoretical results.Circles are experimental data for truncated CENP-E at saturating ATP(2 mM)measured by Gudimchuk et al[46].

Table 1.Parameter values for kinesin motors.

2.2.Ndc80 protein

Under no external force,Ndc80 makes an unbiased onedimensional diffusion on MT[16,17].We still denote byε0the potential energy of Ndc80 in the initial position,byε1the potential energy after Ndc80 takes the forward step and byε-1the potential energy after Ndc80 takes the backward step.As done in the literature[15,56],the forward and backward stepping rates can be written as

3.Results

3.1.The kinesin stalk can rotate freely in a large range of the angle relative to its head

Due to the flexibility of the NL,it is argued here that the relatively rigid rod-like stalk that connects the two heads via their NLs can rotate freely in a large range of angle relative to the heads.To validate this argument,we resort to all-atom MD simulations(see Methods section).The system for our MD simulations is the x-ray structure of the single kinesin head(pdb:2KIN)[57].The structure of 2KIN is composed of a motor domain with the docked NL and an α7 helix that can form the coiled-coil stalk with that of the partner head.

We perform MD simulations of the system for a time of 250 ns(see movie S1 in supplementary information).In figures 2(a)–(d)we show some snapshots of the structure at different simulation times,and in figure 2(e)we show the temporal evolution of the angles of the helix relative to the three planes,y’o’z’,x’o’y’and x’o’z’,that are fixed relative to the head.From figure 2 and movie S1 it is seen that indeed the α7 helix can rotate freely in a large range of the angle relative to the head.The rotation of the α7 helix relative to the head is equivalent to that of the stalk relative to the two heads of the motor in the dimeric form.This large rotation of the stalk can have a large effect on the cargo transport by multiple kinesin motors and MAPs.

3.2.Velocity of cargo driven by multiple kinesin motors of truncated stalks

Figure 2.All-atom MD simulation results for the kinesin head with a segment of α7 helix.(a)–(d)Structures at different simulation times,with t=0 corresponding to the initial simulation time.The head is drawn in cyan and the α7 helix is drawn in red.(e)Temporal evolution of the angle of the α7 helix relative to the three planes,y’o’z’,x’o’y’ and x’o’z’.

As shown in the above section,the stalk of the kinesin motor can rotate freely in a large range of angle relative to the two heads.This implies that for a cargo driven by multiple kinesin motors moving on an MT filament,as schematically shown in figure 3,when the angle of the stalk of a motor(e.g.motor 2)relative to the y axis in xoy plane is in a range of-θ0＜θ＜θ0,within which the stalk can rotate freely,the motor can move freely on the MT filament in a range of-l0＜x＜l0,within which its stalk is not stretched.In other words,during the processive stepping in the range of–l0＜x＜l0the motor experiences no external force.By contrast,if the angle of the stalk of a motor(e.g.motor 1 or motor 3)relative to the y axis in xoy plane is smaller than-θ0or larger thanθ0,namely with∣x∣＞l0,the motor experiences the external force of the magnitude ofK1(∣x∣-l0)arising from the stretching of the stalk,where K1is the spring coefficient of the stalk.Under this consideration,we study numerically the cooperative transport of the cargo by multiple kinesin motors in this section.We use both the constant-force model and the varying-force model.

The simulation methods are described in the Methods section.The parameter values for the kinesin motors are listed in table 1.After detachment from MT,the rebinding rate of the kinesin-1 motor is taken as μ = 5 s-1(table 1),which is consistent with the available experimental data[58],and the rebinding rate of the CENP-E motor is taken as μ = 0.4 s-1(table 1),as taken before[17].The spring coefficient of the stalk,K1,is chosen as follows.In the full-length form,kinesin-1 has a stalk length of about 110 nm[34].The available experimental data indicated that the spring coefficient of the stalk of the full-length kinesin-1 is about 0.3 pN nm-1[52].In the experiments of Chakraborty et al[17],the truncated kinesin-1 motors,CENP-1 motors and Ndc80 proteins of the stalk length of about 30 nm were used.Here,we also consider the truncated kinesin motors having a stalk length of 30 nm.The spring coefficient of the truncated kinesin-1 and CENP-E is then estimated to be about 1.1 pN nm-1.Thus,we take K1=1 pN nm-1in the calculation.

First,we use the constant-force model,namely,we use equation(1)to calculate the forward stepping rate of one motor.In figure 4(a)(left panel)we show the calculated results of the velocity of the cargo driven by multiple kinesin-1 motors versus the number,N,of the motors connected to the cargo for different values of l0.In some prior theoretical and numerical studies[27–29],it was considered that no range on the MT filament is present,in which the motors can step freely without experiencing any external force.The results with l0=0 shown in figure 4(a)correspond to the case of those prior studies.For comparison,the available experimental data on the velocity versus the number of kinesin-1 motors are shown in the right panel of figure 4(a)[25].The corresponding calculated results for CENP-E are shown in figure 4(b)(left panel).For comparison,the available experimental data on the velocity versus the concentration of CENP-E motors are shown in the right panel of figure 4(b)[26].From figure 4(a)(left panel)it is seen that for any value of l0,the calculated velocity of the cooperative transport by multiple kinesin-1 motors increases slightly with the increase of the motor number N.From figure 4(b)(left panel)it is seen that for any value of l0,the calculated velocity of the cooperative transport by multiple CENP-E motors is insensitive to the motor number N.By contrast,the experimental data(right panel of figure 4(b))showed that the velocity of the cooperative transport by multiple CENP-E motors decreases evidently with the increase of the motor concentration and becomes leveled off at high motor concentrations[26](noting that the velocity becoming leveled off at high motor concentrations implies that the decrease of the velocity does not arise from crowding of the motors).The inconsistency between the numerical and experimental results of CENP-E motors indicates that the constant-force model is not reasonable.

Figure 3.Schematic diagram of the cargo driven by three kinesin motors.When the angle of the stalk of a motor(e.g.motor 2)relative to the y axis in xoy plane is in the range of- θ0＜θ＜θ 0,within which the stalk of the motor is not stretched,the motor can move freely on the MT filament in the range of–l0＜x＜l0 without experiencing the external force arising from the stretching of the stalk.Here,motor 2 is drawn in two different positions,with one at x=-l 0 and the other at x=l0.When the angle of the stalk of a motor(e.g.motor 1 or motor 3)relative to the y axis in xoy plane is smaller than- θ0 or larger thanθ 0 ,the motor experiences the external force arising from the stretching of the stalk.

Second,we use the varying-force model,namely,we use equation(4)to calculate the forward stepping rate of the moving motor.The free energy changeε1-ε0in equation(4)corresponds to the change in the total elastic energy of stretching of all stalks of the motors when the moving motor takes one step because the stepping of the motor is required to stretch all stalks of the motors involved.In figure 5(a)(left panel)we show the calculated results of the velocity of the cargo by multiple kinesin-1 motors versus the motor number N for different values of l0.For comparison,the available experimental data on the velocity versus the number of kinesin-1 motors are shown in the right panel of figure 5(a)[25].From figure 5(a)(left panel)it is seen that with l0=0 the velocity decreases evidently with the increase of the motor number.Interestingly,with l0≥d,the velocity of the cargo is insensitive to the motor number.This implies that for the cooperative transport by multiple kinesin-1 motors,the presence of a range in which the motors can step freely without experiencing any external force has the advantage over the absence of such a range.The corresponding calculated results for CENP-E are shown in figure 5(b)(left panel).For comparison,the available experimental data on the velocity versus the concentration of the CENP-E motors are shown in the right panel of figure 5(b)[26].From figure 5(b)(left panel)it is seen that with l0=0 the velocity decreases largely with the increase of the motor number.With l0≥d/2 the velocity at large N has a much smaller decrease relative to that at N=1 than with l0=0.The larger the l0is,the smaller the velocity decreases.This implies that for the cooperative transport by multiple CENP-E motors,the presence of a range in which the motors can step freely without experiencing any external force also has the advantage over the absence of such a range.More interestingly,from figure 5 it is seen that for both the cooperative transport by kinesin-1 motors and that by CENPE motors,the calculated results on the curve form of velocity versus N with l0=d(red dots in the left panels)are in good agreement with the experimental results.Speaking concretely,for the kinesin-1 motors,both the calculated velocity with l0=d and the experimental velocity are kept nearly unchanged with the change of the motor number or concentration.Note that by tuning parameter k(+)(with k(+)=85 s-1)while keeping other parameters unchanged for the kinesin-1 motors we can make the calculated velocity(with l0=d)be in quantitative agreement with the experimental data(see figure S2 in supplementary information).For the CENP-E motors,both the calculated velocity with l0=d at N=1 and that at large N are in quantitative agreement with the experimental velocity at N=1 or at the very low motor concentration and that at the high motor concentration,respectively.These results support strongly the varying-force model.

Taken together,the studies in this section show that the varying-force model is more reasonable than the constantforce model to describe the cooperative transport by multiple kinesin motors.Thus,in the following studies,we will only use the varying-force model.

Figure 4.Results for velocity of cooperative transport by multiple kinesin motors of the truncated stalks studied with the constant-force model.(a)The results for the case of kinesin-1 motors.The left panel shows the calculated velocity versus the number of motors connected to the cargo for different values of l0.The right panel shows the experimental data taken from Derr et al[25].(b)The results for the case of CENP-E motors.The left panel shows the calculated velocity versus the number of motors connected to the cargo for different values of l0.The right panel shows the experimental data taken from Chakraborty et al[26].

3.3.Velocity of cargo driven by multiple kinesin motors and Ndc80 proteins of truncated stalks

In this section,we study the cooperative cargo transport by multiple kinesin motors and Ndc80 proteins with truncated stalks of the length of 30 nm,as done in the experiments of Chakraborty et al[17]Thus,we take the spring coefficient K2=1 pN nm-1for Ndc80 proteins having the same value as K1=1 pN nm-1for motors,as done in the above section.For the motor,we consider that there is a range of–l0＜x＜l0on the MT filament,in which the motor can step freely without experiencing any external force.Different from the structure of the kinesin motor,showing that the connection between the stalk and head is flexible,the structure of Ndc80 protein indicated that the connection between the stalk and head is not flexible[59–61].Thus,for the Ndc80 protein,there is no range on the MT filament,in which the protein can diffuse freely without experiencing any external force.We use the varying-force model here.As discussed in the above section,the energy changesε1-ε0andε-1-ε0in equations(4),(10)and(11)correspond to the change in the total elastic energy of stretching of all stalks of the motors and Ndc80 proteins when the moving motor or Ndc80 takes one step.The simulation methods are described in the Methods section.The parameter values for the kinesin motors are given in table 1.The parameter values for Ndc80 are1264 s-1andλ=0.5,as mentioned above.For simplicity,the detachment of Ndc80 is not considered here.In the system,the number of motors connected to the cargo,N,is kept fixed while the number of Ndc80 proteins connected to the cargo is varied,and thus the ratio of Ndc80 to kinesin-1 number is varied,as done in the experiments of Chakraborty et al[17].As the available experimental data indicated about 8–30 Ndc80 complexes per microtubule at kinetochores[20,62,63],we take N=9,15,24,30 in our calculations.

Figure 5.Results for velocity of cooperative transport by multiple kinesin motors of the truncated stalks studied with the varying-force model.(a)The results for the case of kinesin-1 motors.The left panel shows the calculated velocity versus the number of motors connected to the cargo for different values of l0.The right panel shows the experimental data taken from Derr et al[25].(b)The results for the case of CENP-E motors.The left panel shows the calculated velocity versus the number of motors connected to the cargo for different values of l0.The right panel shows the experimental data taken from Chakraborty et al[26].

In figures 6(a)–(c),we show the calculated results of the velocity of the cargo driven by multiple kinesin-1 motors and Ndc80 proteins versus the ratio of Ndc80 to kinesin-1 number for different values of l0and different kinesin-1 number N,where the velocity is normalized by the velocity in the absence of Ndc80 with the same N.The corresponding results for cargo transport by multiple CENP-E motors and Ndc80 proteins are shown in figures 6(d)–(f).For comparison,in figure 6 we also show the corresponding experimental data measured by Chakraborty et al[17].From figure 6 it is seen interestingly that the results for the normalized velocity are insensitive to l0.This can be understood as follows.Although at a given ratio of Ndc80 to motor number the cargo velocity by multiple motors and Ndc80 proteins increases with the increase of l0,the cargo velocity by multiple motors in the absence of Ndc80 also increases with the increase of l0(see figure 5(b)).This thus gives the normalized velocity at a given ratio of Ndc80 to motor number to be insensitive to l0.More interestingly,from figures 6(a)–(c)it is seen that for the case of cargo transport by multiple kinesin-1 motors and Ndc80 proteins,the results are in good agreement with the experimental data for any N(≥9).When the ratio of Ndc80 to kinesin-1 number becomes equal to one,implying that both the motors and Ndc80 proteins have the same number,the cargo velocity is reduced by about 5-fold relative to that when the ratio of Ndc80 to kinesin-1 number is equal to zero,implying that only the motors are present.From figures 6(d)–(f)it is seen that for the case of cargo transport by multiple CENP-E motors and Ndc80 proteins,the results are in good agreement with the experimental data when N≥24.When the ratio of Ndc80 to CENP-E number becomes equal to one,namely both the motors and Ndc80 proteins having the same number,the cargo velocity becomes nearly equal to zero.The underlying mechanism of the sensitive decrease in the velocity to the increase in the ratio of Ndc80 to kinesin number is as follows.Since the short truncated stalks for both motors and Ndc80 proteins have the large spring coefficients K1and K2,the energy changeε1-ε0caused by the forward stepping of the moving motor is large,thus resulting in a large decrease in the forward stepping rate of the motor,as noted from equation(4).

In addition,it is noted that at a given ratio of Ndc80 to CENP-E number the normalized velocity for l0=d=8.2 nm(figure 6(d))is very similar to that for l0=2d=16.4 nm(figure 6(f)).Moreover,it is noted that at a given ratio of Ndc80 to CENP-E number the difference between the normalized velocity for l0=d=8.2 nm(figure 6(d))and that for l0=3d/2=12.3 nm(figure 6(e))is larger than the difference between the normalized velocity for l0=d=8.2 nm(figure 6(d))and that for l0=2d=16.4 nm(figure 6(f)).These results can be understood as follows.Consider that one motor is initially at the equilibrium position x=0 withθ=0(see figure 3).For the case of l0=nd(n is an integer),during the movement from x=0 to x=l0=nd the motor experiences no force arising from the stretching of the stalks of the motor and Ndc80,while throughout the next whole step from x=l0=nd to x=(n+1)d the movement of the motor is affected by the energy change arising from the stretching of the stalks of the motor and Ndc80.Thus,it is approximately expected that for any n the normalized velocity of the cargo is very similar.By contrast,for the case of l0=(n+1/2)d,during the movement from x=0 to x=l0=(n+1/2)d the motor experiences no force arising from the stretching of the stalks of the motor and Ndc80,while throughout the next half step from x=(n+1/2)d to x=(n+1)d the movement of the motor is affected by the energy change arising from the stretching of the stalks of the motor and Ndc80.Thus,it is approximately expected that the difference between the normalized velocity for l0=d(figure 6(d))and that for l0=3d/2(figure 6(e))is larger than the difference between the normalized velocity for l0=d(figure 6(d))and that for l0=2d(figure 6(f)).

It is noted that for simplicity we have not considered the detachment of the Ndc80 protein in figure 6.With the consideration of the detachment of the Ndc80 protein,the calculated results are also in good agreement with the experimental data(see section S1 in supplementary information).

Figure 6.Results for velocity of cooperative transport by multiple kinesin motors and Ndc80 proteins of the truncated stalks studied with the varying-force model.The experimental data are taken from Chakraborty et al[17],with the errors of the experimental data being calculated with where v and v0 denote the velocities with the presence and absence of Ndc80 proteins,respectively.(a)–(c)The normalized velocity of the cooperative transport by multiple kinesin-1 motors and Ndc80 proteins versus the ratio of Ndc80 to kinesin-1 number for different values of l0 and different kinesin-1 number N,where the velocity is normalized by the velocity in the absence of Ndc80 with the same N.(d)–(f)The normalized velocity of the cooperative transport by multiple CENP-E motors and Ndc80 proteins versus the ratio of Ndc80 to CENP-E number for different values of l0 and different CENP-E number N,where the velocity is normalized by the velocity in the absence of Ndc80 with the same N.

3.4.Velocity of cargo driven by multiple kinesin motors and Ndc80 proteins of full-length stalks

In this section,we make similar studies to those presented in the above section except that we use the full-length kinesin motors and Ndc80 proteins.Thus,we take K=0.3 pN nm-1for the full-length kinesin-1 motors,as mentioned above.Since the full-length CENP-E has a length of about 200–230 nm[64],we take K1=0.15 pN nm-1for CENP-E motors.Since the full-length Ndc80 has a length of about 60 nm[57,65],we take K2=0.55 pN nm-1for Ndc80 proteins.For simplicity,we take the truncated and full-length kinesin motors having the same values for other parameters,as used in the above section.

In figures 7(a)–(c),we show the calculated results of the velocity of the cargo driven by multiple full-length kinesin-1 motors and Ndc80 proteins versus the ratio of Ndc80 to kinesin-1 number for different values of l0and different kinesin-1 number N,where the velocity is normalized by the velocity in the absence of Ndc80 with the same N.The corresponding results for cargo transport by multiple fulllength CENP-E motors and Ndc80 proteins are shown in figures 7(d)–(f).From figure 7 it is seen interestingly that the velocity decreases only slowly with the increase in the ratio of Ndc80 to kinesin number,which is contrary to the case of the kinesin motors of the truncated stalks(figure 6).The small decrease of the velocity for the case of full-length stalks is due to the small elastic energy changeε1-ε0when one motor takes a forward step,which is in turn due to the small spring coefficients K1and K2of the stalks for both motors and Ndc80 proteins.Thus,the kinesin motors in the full-length form,together with MAPs,can transport cargo much more efficiently than the truncated ones with the shorter stalk length.This gives an explanation of why the full-length kinesin motor,and in particular the CENP-E motor,has a long stalk.

4.Discussion

In some prior theoretical and numerical studies of cargo transport by multiple kinesin motors[27–29],it was assumed that there is an equilibrium position of one motor on MT,at which the stalk of the motor is not stretched and thus the motor experiences no external force.However,no range for the equilibrium position of the motor is present,within which the stalk of the motor is not stretched and thus the motor can step freely without experiencing any external force.On the contrary,the all-atom MD simulations here indicate that the stalk of the kinesin motor can rotate freely in a large range of the angle relative to the head(figure 2 and Movie S1).This implies that for the cargo transport by multiple kinesin motors and MAPs,a range for the equilibrium position of a motor is present,within which the stalk of the motor is not stretched and thus the motor can step freely without experiencing any external force,as argued in other prior studies[30–36].The studies here show that the presence of the range of the equilibrium position has an advantage over the absence of such the range in the cooperative transport of the cargo by multiple motors,with the velocity of the cargo for the former being larger than that for the latter(figure 5).

In prior theoretical and numerical studies of cargo transport by multiple kinesin motors[27,28,34–41],the force dependence of the velocity for the single motor was generally described by the constant-force model,in which it was assumed that during the moving period in one step the motor experiences the constant external force F that is equal to that during the dwell period before the moving period.The studies here show that with the constant-force model for the single motor,the calculated results of the velocity of the cargo transport by multiple CENP-E motors versus the motor number or concentration are inconsistent with the available experimental data(figure 4).Recently,the varying-force model was used to study the dynamics of the single motor in the single molecule optical trapping experiments,in which it was proposed that during the moving period in one step the motor experiences the varying external force although during the dwell period the motor experiences the constant external force[32,47].Here,with the varying-force model for the motor,both the calculated results on the motor-number dependence of the velocity of the cargo transport by multiple kinesin-1 motors and those by multiple CENP-E motors are in good agreement with the available experimental data(figure 5).

It is noted that for the cargo transport by multiple kinesin motors and Ndc80 proteins with the truncated stalks,the cargo velocity decreases significantly with the increase of the ratio of Ndc80 to kinesin number.For example,in the case of kinesin-1,when the ratio is 0.8 the cargo velocity normalized by that with the ratio equal to zero is only about 0.37(figure 6).For the case of CENP-E,when the ratio is 1 the cargo velocity becomes nearly zero(figure 6).Here,we give a quantitative explanation of these results.These results indicate that when the ratio of Ndc80 to kinesin number is high the cargo transport is very inefficient.In particular,when the ratio of Ndc80 number to CENP-E number is 1,the cargo is unable to be transported.Interestingly,we show that the velocity of the cargo transport by multiple kinesin motors and Ndc80 proteins with the full-length stalks decreases only slowly with the increase of the ratio of Ndc80 to kinesin number(figure 7).In particular,even when the ratio of Ndc80 to CENP-E number is 1,the cargo velocity is decreased by only about 1.2-fold,implying a very efficient cargo transport.This thus gives an explanation of why CENP-E has a long stalk.

Figure 7.Results for velocity of cooperative transport by multiple kinesin motors and Ndc80 proteins of the full-length stalks studied with the varying-force model.For comparison,the experimental data that are the same as those shown in figure 6 are reshown here.(a)–(c)The normalized velocity of the cooperative transport by multiple kinesin-1 motors and Ndc80 proteins versus the ratio of Ndc80 to kinesin-1 number for different values of l0 and different kinesin-1 number N,where the velocity is normalized by the velocity in the absence of Ndc80 with the same N.(d)–(f)The normalized velocity of the cooperative transport by multiple CENP-E motors and Ndc80 proteins versus the ratio of Ndc80 to CENP-E number for different values of l0 and different CENP-E number N,where the velocity is normalized by thevelocity in the absence of Ndc80 with the same N.

In conclusion,using the varying-force model and with the argument of the presence of the range of the equilibrium position of the motor on MT,which is verified by our allatom MD simulations,our calculated results for the velocity in the cooperative transport by multiple kinesin motors versus the motor number reproduce well the available experimental data.Moreover,we study numerically the cooperative cargo transport by multiple kinesin motors and Ndc80 proteins with the truncated stalks of the large spring coefficient,explaining quantitatively the available experimental results on the velocity versus the ratio of Ndc80 to kinesin number.The sensitive decrease of the velocity to the ratio of Ndc80 to motor number is due to the large decrease of the forward stepping rate of the moving motor,which arises from the large change in the elastic potential energy of stretching the stalks of the motors and Ndc80 proteins when the moving motor makes a forward stepping.The cooperative cargo transport by multiple kinesin motors and Ndc80 proteins with the full-length stalks of the small spring coefficient is also studied numerically.The results show that the velocity of the cargo decreases only slowly with the increase in the ratio of Ndc80 to kinesin number.This is due to the small change in the elastic potential energy of stretching the stalks of the motors and Ndc80 proteins when the moving motor makes a forward stepping.Our results thus give an explanation of why the full-length kinesin motor,and in particular the CENP-E motor,has a long stalk.The predicted results with the full-length stalks can be easily tested by future experiments.

5.Methods

5.1.All-atom MD simulations

We use the structure of the single rat kinesin head(pdb:2KIN)[57]as the starting model in our all-atom MD simulations.The MD simulation procedure is described in detail in our previous works[66,67],which is re-described briefly here.The missing atoms in the structure are added by using the software Swiss-PdbViewer.Counter-ions are added to neutralize the system.The kinesin head is solvated in a water box(using the TIP3P water model[68]).In our simulations,the distance between the protein complex and the boundary of the box is about 4 nm to avoid the edge effect.The MD simulations are run in Gromacs4.6[69]using the AMBER99SB force field[70]with a time step of 2 fs.All chemical bonds are constrained using the LINCS algorithm[71].The cutoff for van der Waals interaction and short-range electrostatics interaction is set as 1.4 nm.The long-range electrostatic interaction is calculated using the particle mesh Ewald(PME)algorithm[72].Velocity-rescaling temperature coupling[73]and Berendsen pressure coupling[74]are used.We first perform the energy minimizations twice.The energy minimizations are performed for 5×104steps using the steepest descent method.Then the systems are equilibrated for 100 ps at 300 K and 1 bar pressure in the NVT ensemble and NPT ensemble,respectively.In the MD simulations,we restrict residues 223–227 in the head.These restricted residues are far away from the neck so that the restriction of them would have no influence on the motion of the neck domain.The simulated results are analyzed by VMD1.9.2[75].

5.2.Monte-Carlo simulations

With the above-mentioned Monte-Carlo simulation procedure,the method for simulations of the cargo transport by multiple kinesin motors and Ndc80 proteins is described in the supplementary information(section S2).With the varying-force model for kinesin motors,the method to calculate the stepping rate of one motor in the cargo transport by multiple kinesin motors and Ndc80 proteins is described in the supplementary information(section S3).The method to calculate the stepping rate of one Ndc80 protein in the cargo transport by multiple kinesin motors and Ndc80 proteins is also described in the supplementary information(section S3).

ORCID iDs