Skip to main content
Advertisement
  • Loading metrics

Cargo surface fluidity can reduce inter-motor mechanical interference, promote load-sharing and enhance processivity in teams of molecular motors

  • Niranjan Sarpangala,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Physics, and Center for Cellular and Biomolecular Machines, University of California, Merced, California, United States of America

  • Ajay Gopinathan

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    agopinathan@ucmerced.edu

    Affiliation Department of Physics, and Center for Cellular and Biomolecular Machines, University of California, Merced, California, United States of America

Abstract

In cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of cytoskeletal filaments. How motors that typically mechanically interfere with each other, work together as teams is unclear. Here we explored the possibility that purely physical mechanisms, such as cargo surface fluidity, may potentially enhance teamwork, both at the single motor and cargo level. To explore these mechanisms, we developed a three dimensional simulation of cargo transport along microtubules by teams of kinesin-1 motors. We accounted for cargo membrane fluidity by explicitly simulating the Brownian dynamics of motors on the cargo surface and considered both the load and ATP dependence of single motor functioning. Our simulations show that surface fluidity could lead to the reduction of negative mechanical interference between kinesins and enhanced load sharing thereby increasing the average duration of single motors on the filament. This, along with a cooperative increase in on-rates as more motors bind leads to enhanced collective processivity. At the cargo level, surface fluidity makes more motors available for binding, which can act synergistically with the above effects to further increase transport distances though this effect is significant only at low ATP or high motor density. Additionally, the fluid surface allows for the clustering of motors at a well defined location on the surface relative to the microtubule and the fluid-coupled motors can exert more collective force per motor against loads. Our work on understanding how teamwork arises in cargo-coupled motors allows us to connect single motor properties to overall transport, sheds new light on cellular processes, reconciles existing observations, encourages new experimental validation efforts and can also suggest new ways of improving the transport of artificial cargo powered by motor teams.

Author summary

In cells, multiple molecular motors work together as teams to carry cargoes such as vesicles and organelles over long distances to their destinations by stepping along a network of protein filaments. How do motors that interfere with each other when coupled rigidly in laboratory experiments, function well as teams within cells? In this paper, we show, using computer simulations, that the fluid surfaces of cellular cargo reduce the mechanical interference between motors allowing better load sharing and increased duration of motors on the filament, thereby increasing the distance over which they can carry cargo. We additionally found that motors pull cargo closer to the filament increasing the attachment rates of other unbound motors. These effects act synergistically with an increased availability of motors due to fluidity to further increase travel distances. Fluid surfaces also allow clustering of the motors and increased collective force against load. Our work on understanding how teamwork arises in mechanically coupled motors sheds new light on cellular processes, reconciles existing observations, encourages new experimental validation efforts and can also suggest new ways of improving transport of artificial cargo powered by motor teams.

Introduction

In eukaryotic cells, the intracellular transport of material between various organelles and the cellular and nuclear membranes is critical for cellular function and is an active process facilitated by the consumption of energy [1, 2]. A variety of motor proteins including myosins, kinesins and dyneins convert chemical energy from ATP hydrolysis into directed stepping along cytoskeletal protein filaments such as actin filaments and microtubules [1, 3]. These motors move in different directions along the filaments typically carrying lipid bilayer vesicles packed with proteins and signaling molecules or even membrane bound organelles such as Golgi and mitochondria [1, 2]. Defects in motor function can impair normal cell functioning and lead to a variety of pathologies including neurodegenerative diseases like Alzheimer’s disease and ALS [49]. Given its importance in vivo, there has been a lot of work done over the years in understanding motor function in the transport context ranging from single molecule studies that elucidate the detailed mechanistic workings of motors [1017] to the functioning of multiple motors and their co-ordination [1827]. Several theoretical and in vitro experimental studies [2832] indicate that the collective behavior of motors is critically influenced by their coupling to each other resulting in observable effects in transport speeds and run lengths. For example, it is known that non-cooperative kinesin motors coupled together by a rigid cargo interfere with each other’s functioning leading to enhanced detachments and lowered run lengths [2830]. How then are in vivo cargoes such as membrane-bound vesicles and organelles typically carried over long distances by teams of kinesin motors [33, 34]? It is possible that a fluid membrane simply makes more motors available for binding at the filament either by dynamic clustering of motors [35, 36] or membrane induced clustering of motors into micro-domains [37, 38] thereby enhancing processivity. Alternatively or in addition, there could be other coupling-dependent physical mechanisms that directly affect single motor functioning and thereby promote teamwork. Studies on transport in vivo could address these questions directly, but the environmental complexity makes it difficult to disentangle fundamental phenomena of physical origin from the effects of other regulatory mechanisms [39]. In vitro systems, which offer cleaner insights, include microtubule gliding assays on flat bilayers with embedded motors [35, 4042], motor driven nanotube extraction from vesicles [43] and more recently even membrane covered beads carried by kinesin motors [44]. Studies of these systems and related models have generated a variety of interesting and sometimes conflicting results on cargo run lengths and velocities in the presence of membrane fluidity [36, 40, 44, 45]. These results are also obtained under different conditions of cargo geometry, motor density, cargo surface fluidity and environmental factors like ATP concentration. In order to reconcile them, we need an understanding of the relative importance and interplay of these factors and their effects on the number of engaged motors and the loads they experience, which in turn influences their collective speed and processivity. More importantly, to uncover any physical mechanisms present in teams of motors that directly affect single motor functioning and thereby enhance teamwork, we need to be able to monitor both single motor dynamics and collective cargo transport simultaneously.

Here, we used Brownian dynamics simulations to accomplish this goal. Briefly, we model the transport of a spherical cargo with a fluid surface carried by multiple kinesin-1 motors along a microtubule (MT) (see Fig 1 and details in Materials and methods). We explicitly accounted for the movement of the attachment points of the motors on the cargo surface (Ai in Fig 1c) by allowing free diffusion for unbound motors and diffusion biased by exerted forces (spring forces Fi in Fig 1c) for bound motors, with cargo surface fluidity determining the diffusion constant. We also accounted for the force and ATP dependence of the bound motor’s off-rate and stepping rate. The spherical cargo’s position was then updated by Brownian dynamics depending on the net force exerted by all bound motors. We also incorporated the rotational diffusion of the cargo due to the torque from motor forces and thermal fluctuations. Our model is overall quite similar to the Brownian dynamics model of Bovyn et. al. [36] which also implemented overdamped Langevin dynamics of motors on an undeformable, spherical lipid cargo surface. While there are minor differences between our models such as in the details of the load dependence of the motor unbinding rates and the incorporation of ATP, we do not expect significant quantitative differences. Indeed, measurements of the same quantities such as run lengths and single motor binding times are consistent between our studies. The major difference with our work lies in the focus of their study which was the effect of surface fluidity on the availability of motors for both clustered and dispersed motors and how this influenced the trade-off between binding rates from solution and run length. While we also consider the availability of motors, a major focus in this paper is on the effect of fluidity at the single motor level. This includes effects on the mechanical interference between motors which affects their load-sharing and thereby single motor processivity as well as collective force generation. Another difference is that we specifically consider the effects of varying ATP concentration. Finally, we study the relative importance of these effects and how they act synergistically at the level of the collective to influence cargo transport.

thumbnail
Fig 1. Schematics of cargoes considered in the study.

(a) Rigid cargo (membrane-free cargo). Molecular motors are permanently attached to random locations on the cargo surface. (b) Fluid/Lipid cargo (membrane-enclosed cargo). Molecular motors diffuse on the cargo surface. Inset of Fig 1b The lipid bilayer that forms the fluid cargo surface. (c) Schematic that shows how we compute the forces generated by motors. X is the center of mass of the cargo, Ai and Hi represent motor anchor and motor head positions. Motors exert spring like forces only when their lengths Li exceed their rest length Lmot. Here, motor 1 experiences a hindering load (f1 is negative), while motor 2 experiences an assistive force (f2 is positive).

https://doi.org/10.1371/journal.pcbi.1010217.g001

Specifically, our model enabled us to focus on the dynamics of the motors as they diffused on the surface, their binding and unbinding from the filament and the forces experienced by individual motors, as a function of motor diffusion constant, motor density and ATP concentration. We connected these characteristics at the individual motor level to the effect of surface fluidity (inverse viscosity proportional to the diffusion constant) on the collective behavior of motors at the cargo level by analyzing transport properties such as the average number of bound motors, collective force generated against load, speed and distance traveled (run length). We uncovered two salient features, the reduction of interference and the cooperative increase of on-rates. Our simulations showed explicitly that surface fluidity leads to the reduction of negative mechanical interference between kinesins, characterized by lower forces on individual motors and a significant drop in antagonistic forces between motors. This allows teams of fluid-coupled motors to more fully exploit load sharing without inter-motor interference. This decreases single motor off-rates and increases processivity. Our simulations also showed that the fluid surface allows for the clustering of motors at a well defined location on the surface relative to the microtubule and that the fluid-coupled motors can exert more collective force per motor against loads. Interestingly, increasing numbers of bound motors pull the cargo closer to the microtubule, increasing the on-rates of unbound motors, resulting in a cooperative increase in bound motor numbers that depends on 3D cargo geometry. Surface fluidity also makes more motors available for binding, as expected, and as indicated by previous studies [35, 36, 45]. However, this effect, by itself, is significant only at lower ATP concentrations (and/or very high motor numbers) when the effective timescale for diffusion and binding is less than the unbinding time. In fact, we estimate that the reduction in interference is relatively the most important effect under physiological conditions of high ATP and low motor numbers. Taken altogether these effects can cooperatively result in increased processivity and collective mechanical force production against load with an increase in fluidity for teams of molecular motors, with both in vivo and technological implications.

Results and discussion

Surface fluidity reduces negative interference between bound kinesin motors

We first addressed the question of whether the fluidity of the cargo surface has any influence on the functioning of a bound kinesin motor. We expected that the motor attachment point’s freedom to move on the lipid cargo surface in response to forces could suppress strain and reduce the forces experienced by the motor. We also expected that the surface fluidity would allow bound motors to alter their spatial distribution relative to the cargo and each other thereby also changing the forces they experienced. We started by looking at the positions of motor heads on the MT and anchor points on the cargo during a typical cargo run for a rigid and a lipid cargo (Fig 2a; see S1 Fig for a longer time window). From the motor head position subplots, it can be seen that motors on a rigid cargo are more spread apart and can be located opposite to each other relative to the cargo center of mass while on lipid cargos, bound motors tend to be closer and on the same side, in front of the cargo. This supports our expectation that the lipid membrane can alter the relative positioning of the motors which could lead to different strain forces on the motors.

thumbnail
Fig 2. Surface fluidity reduces the negative interference between bound kinesin motors which leads to reduction in motor off-rate.

(a) Typical trajectories for a rigid (top row) and a lipid cargo (second row). Motor head position subplots (middle column) include the x-positions of the center of mass of cargo (black lines) and head positions of motors on the microtubule (colored lines). The corresponding x-component of motor anchor positions (right column) on the cargo are also shown. (b) Distributions of the absolute magnitude of the force experienced by motors in rigid and lipid cargoes (including without cargo rotation) at a fixed number of bound motors (n = 3). Distribution for singly bound motor carrying a rigid cargo is also shown for comparison. We first collect S = 10000 force values randomly drawn from 200 cargo runs. We then bootstrap this (bootstrap sample size = 1,000 and the number of bootstrap samples = 10) to get the mean distribution and standard error of the mean. (c) Distribution of magnitude of motor forces when they specifically experience hindering loads (f < 0) and assistive loads (f > 0). Cartoon inside shows a possible scenario where the colored motor experiences hindering and assistive load. (d) The mean variance of forces among bound motors in rigid (D = 0) and lipid cargoes (D = 1 μm2 s−1) with (triangles) and without cargo rotation (circles). Mean variance was calculated as where the summation is over all the sample time points when the number of bound motors is n. is the variance in the magnitude of force experienced by the n bound motors at ith sample. Sample size S = 10000 drawn from 200 cargo runs. (e) Mean value of the correlation between the x-components of motor forces, , averaged over motor pairs and time points with (triangles) and without cargo rotation (circles). Sample size S = 10000 drawn from 200 cargo runs. (f) Average value of motor off-rate as a function of fluidity of cargo surface with (triangles) and without cargo rotation (circles). (S = 580 from 200 cargo runs). Data for all the plots were obtained from the simulation of transport of cargoes with N = 16 motors at high ATP concentration of 2 mM. Data sampling rate was 100 s−1.

https://doi.org/10.1371/journal.pcbi.1010217.g002

To explore the forces experienced by motors more quantitatively, we computed the distribution of forces experienced by a single motor when it is part of a team of n motors that are simultaneously bound to the microtubule. This distribution, shown for n = 3 bound motors in Fig 2b, is broader, reflecting larger forces, for a motor on a rigid cargo compared to any of the lipid cargoes. The lipid cargo motor force distributions is, in fact, comparable to that for a single bound motor carrying a rigid cargo. A rigid coupling of motors, therefore, results in individual bound motors experiencing higher forces than when they are singly bound to cargo, indicating negative interference. The introduction of surface fluidity appears to reduce this interference by effectively decoupling the motors, consistent with our expectations. This reduction is also apparent in the distribution of forces experienced both against (hindering) and in (assistive) the direction of processive motion (Fig 2c).

The reduction in forces due to fluidity is also evident in the fraction of non-negligible forces (S2 Fig), which we take to be forces with a magnitude ≥ 1% of Fs, the stall force of the kinesin motor (S3 and S4 Figs show similar fractions in hindering and assistive directions). For rigid cargo, this fraction increases as a function of the number of bound motors, n, indicating negative interference between motors. In the presence of surface fluidity, however, no such increase in the fraction with n is apparent, indicating a reduction of negative interference. In fact, the fraction decreases with n at high surface fluidity (D = 1 μm2 s−1), which is likely due to the decrease in the average force experienced by the motors in the presence of multiple bound motors (S5 Fig). While the decrease in average force occurs for rigid cargo as well, the increase with n due to negative interference dominates.

The variance between the individual forces experienced by simultaneously bound motors is also a measure of negative interference as it quantifies the deviation from perfect load-sharing where the variance should be zero. The mean variance over all sampled instances with n bound motors shows a dramatic decrease at higher fluidity (D = 1 μm2 s−1) compared to the rigid cargo case (Fig 2d, see S5 Appendix for detailed description of this and other force related metrics used in this section). In fact, even the presence of a small amount of fluidity accounts for most of the significant drop in variance (S6 Fig) indicating that fairly good load-sharing can be achieved with modest fluidity. We also quantified negative interference more directly by studying whether individual motors were acting antagonistically, i.e., exerting forces in opposite directions. To do this, we computed the correlation between the x-components of bound motor forces averaged over different motor pairs (j, k) and time (Fig 2e). We see that this correlation is highly negative for rigid cargo indicating opposing motor forces while it is negligible in the presence of fluidity indicating motors on the same side of the cargo pulling in the same direction. Taken together, our results show that interference dominates over load sharing for rigid cargo while increasing surface fluidity reduces this interference and promotes load-sharing.

Next, we addressed the implications of the reduction in negative interference for the functioning of motors in the context of transport. The off-rate of a kinesin motor increases with the magnitude of the force that it experiences and also depends on its direction [15, 23, 4648]. The observation that surface fluidity reduces negative interference leads us to expect that this will also lead to lower off-rates for individual motors. We measured the mean off-rate of motors as the inverse of the mean time spent by individual motors between a binding and subsequent unbinding event. We found that bound motors on fluid cargoes do indeed have lower off-rates than those on rigid cargoes (Fig 2f). Consistent with our findings for the forces and interference, the introduction of a small amount of fluidity accounts for most of the significant drop in off-rates with no significant variation with increasing diffusivity.

The forces that motors feel and the off-rate of motor is a function of the spring constant of the kinesin coiled-coil region. In fact, previous studies [28, 49] have found that off-rate of ensemble of motors and the transport properties of cargo are sensitive to the stiffness of the motor. Our simulations (see S22 Fig) indicate that average single motor off-rate can increase by almost 20% for a change in the motor stiffness from 0.1 pN/nm to 0.5 pN/nm for rigid cargo, while for lipid cargo the change is not significant, which is again consistent with the picture of the lipid membrane suppressing strain in the motors.

We also ran simulations at different cargo sizes to determine how the single motor metrics change. The distribution of forces (S7(a) and S7(b) Fig) shows no clear trend except perhaps for a marginal increase in larger forces with decreasing radii for lipid cargo, which could be due to the steric constraints becoming more important. Interference between motors, quantified by the force correlations (S7(c) Fig), show negligible correlations between motors for lipid cargo at any size, while for rigid cargo, the interference appears to increase with size, consistent with potentially increased spacing between motors. Motor off-rate (S7d Fig) decreases as a function of cargo size for both rigid and lipid cargo but the off-rate of motor in lipid cargo is lower than rigid cargo for all cargo sizes, which is again consistent with picture of reduced interference and greater load-sharing with a lipid cargo.

We also took the opportunity to explore the impact of cargo rotations on motor forces. While cargoes do rotate considerably during the course of a typical cargo run and especially at intersections or obstacles [50, 51], it is not clear whether this rotation impacts transport metrics like motor off-rate, force distribution and cargo runlength. Answering this question is also interesting from a modeling perspective and helpful for improved analytical approximations. Disallowing rotational motion of the cargo in our simulations leads to a slight decrease in the magnitude of motor forces, as is apparent in n = 3 cases shown in Fig 2b and 2c. In our simulations and in in vitro experiments in aqueous media, the load due to viscous forces is low and the origin of higher forces is thermal fluctuations [52] and negative interference between motors. The higher forces experienced by motors in the presence of rotation is therefore likely due to rotational diffusion adding additional noise into the cargo dynamics. It is to be noted that the reduction in forces due to removing cargo rotations is slightly more pronounced for rigid cargo but in both cases the reduction is small compared to the reduction due to fluidity. The variance of motor forces (Fig 2d) also decreases when we remove rotational diffusion because of the decrease in the noise though the effect is not significant for lipid cargo. Similarly, the magnitude of force correlations (Fig 2e) increased when we removed rotation for rigid cargoes while the correlations are essentially the same and close to zero for lipid cargo. These effects arise because the removal of rotational diffusion removes a source of uncorrelated noise. Finally, we compared the change in motor off-rate with and without incorporating rotation in our model (Fig 2f). We notice that rotation may increase the off-rate slightly (by about 5%) for rigid cargo, while the off-rate of motors on lipid cargo is unaffected by cargo rotation. Thus, it seems that ignoring rotation of cargo leads to no qualitative changes and is a perfectly good quantitative approximation for lipid cargo and a reasonable one that slightly underestimates forces and off-rates for motors on a rigid cargo. It is to be noted that rotations can play a much more important role in enhancing binding rates from solution as seen in [36], but for lipid cargo, the diffusion of motors dominates over rotations of the cargo, and therefore rotations play a negligible role. In what follows, results presented include rotation of the cargo and are essentially very similar to the results without rotation unless explicitly noted.

Overall, our modeling indicates that increasing surface fluidity reduces inter-motor interference and promotes load-sharing leading to reduced off-rates and potentially longer processive runs for individual motors.

On-rate of a motor increases with the number of bound motors for fluid cargo

The run length of a cargo carried by a team of motors depends on the interplay between the off-rates and on-rates of individual motors [23]. The on-rate of a motor, typically considered to be independent of the number of other bound motors, depends on its intrinsic binding rate when it is within reach of the MT. The region of the cargo surface (or access area) from where an anchored kinesin can reach the MT is defined, in our model, by the set of points on the surface at a distance less than or equal to the rest length of the motor (Lmot = 57 nm). For a cargo bound to the MT with one or more motors, this access area is a function of the cargo geometry and the proximity of the cargo to the MT [35, 38, 53]. For rigid cargo, motors can only bind if they happen to be anchored in the access area while for fluid cargo, motors can diffuse in and out of the access area. For a given motor type with a fixed intrinsic binding rate and a given cargo geometry (a sphere here), the access area is then simply a function of the distance of the cargo from the MT. Interestingly, our measurements of the average distance of the center of mass of the cargo from the MT indicated that the cargo comes closer to the microtubule as the number of bound motors increases (Fig 3a). This decrease in distance arises from the increase in the net vertical component of the average force due to an increased number of motors which serves to counteract cargo fluctuations away from the MT. We also note that cargoes with low surface fluidity are closer to the MT than the cargoes with high surface fluidity, especially for high numbers of bound motors. This is because motors on fluid cargoes can relax their tension by sliding on the cargo surface allowing for larger fluctuations (S8 and S9 Figs) and therefore a higher average distance from the MT. We also note that rotations of the cargo seem to allow the cargo to get slightly closer to the MT, presumably because the extra forces due to rotational diffusion serve to increase the net vertical component of the average force as well.

thumbnail
Fig 3. On-rate of a motor increases with the number of bound motors for fluid cargo.

(a) Average measured distance (H) of the center of mass of cargo from MT as a function of the number of bound motors and fluidity. Horizontal line indicates the height below which the cargo experiences steric interaction from the MT. Data were obtained from the simulation of transport of cargoes with N = 16 motors at high ATP concentration of 2 mM. Data sampling rate was 100 s−1. (b) Computational measurements of single motor on-rate (green diamonds) and analytical approximation (black dashed line with diamond), (D = 1 μm2 s−1, S = 10000) using the cargo height for a given n from the data in (a). Error bars for all plots represent the standard error of the mean. Rotation of cargo was not considered in the on-rate measurement simulations for (b).

https://doi.org/10.1371/journal.pcbi.1010217.g003

Thus, for both rigid and fluid cargo, as the number of bound motors increases, the cargo approaches the MT leading to an increase in the access area. We expect that this should lead to an increase in the on-rate of an individual unbound motor because its likelihood of being in the access area is increased. We verified this by measuring the on-rate of a motor diffusing on the surface of a cargo held at specific distances from the MT (Fig 3b) that correspond to different numbers of bound motors (from data in Fig 3a for D = 1 μm2 s−1). Here we chose a high diffusion constant for simplicity and measured the mean time required for a randomly placed motor to diffuse and bind to the MT. Indeed, we see that the effective on-rate of the motor, calculated as the inverse of this mean time, increases with decreasing distance corresponding to increasing numbers of bound motors (Fig 3b). To show that the quantitative increase in on-rate was explained by the decreasing distance, we also analytically estimated the change in on-rate. Since the average time spent by a motor in the access area before binding to the microtubule is approximately proportional to the ratio of access area to total area in the high diffusion constant limit, we estimated that the effective on-rate is πad = π0(Sa/ST), where π0 = 5 s−1 is the intrinsic binding rate, Sa is the access area and ST = 4πR2 is the total surface area of the cargo. Such an assumption is reasonable for cargo surfaces with high diffusion constants where the time scale for motor diffusion across the cargo surface is small as seen in previous studies [35, 38, 54]. By computing the access area numerically for different cargo distances, we were able to estimate πad corresponding to different numbers of bound motors. The rather good agreement with the measured values from simulations (Fig 3b) indicates that the increase in on-rate is captured by the effects of increased access area. It is to be noted that the rate at which the cargo transitions from being bound by n to n + 1 motors can be calculated from this on-rate, the off-rate and the number of unbound motors available for binding (details in S2 Appendix and data in S10 Fig).

Our results reveal a positive feedback effect at play. As the number of bound motors increases, the cargo gets closer to the MT and the on-rates of each of the remaining motors increase which further increases the number of bound motors. Interestingly, the apparently small change in cargo height can have a big effect on the overall processivity of the team as discussed in the subsequent section on the contributions of different effects to overall run length.

However, it should be noted that this positive cooperativity that we observe is conditioned on our assumption that an unbound motor can bind to any position on the microtubule at a distance less than the rest length of the motor (Lmot). For situations where a large number of protofilament tracks are not available, such as myosin motors on actin filaments, steric occlusion could be important. For instance, studies [54] on myosin binding to actin also found that the cargo comes closer to the filament with an increase in the number of bound motors, but this leads to a decrease in the on-rate of motors because of increased steric occlusion from the cargo itself.

Surface fluidity increases the availability of motors for binding

A cargo cannot, however, benefit from the positive feedback noted in the last section, unless it makes a sufficient number of motors available for binding. We expect that the diffusion of motors on the cargo surface should lead to a greater availability of motors in the access region and hence a higher number of bound motors compared to rigid cargo with the same overall number of motors. We also expect that this will be true only if τbind—the typical time for any one of the unbound motors to bind to the microtubule is less than τoff = 1/ϵ0—the typical time for a given bound motor to detach. The relative values of these two timescales are set by parameters such as the number of motors (N), the radius of cargo (R), and the ATP concentration. For instance, raising the motor density (say by raising N) will result in a higher effective on-rate and a decreased τbind, while decreasing the ATP concentration, lowers the off-rate and raises τoff.

We first looked at a parameter set (N = 16 [ATP] = 2 mM) satisfying the condition τbind < τoff (see S2 Appendix for details of estimates). Fig 4a shows that, while the average number of bound motors, , increases with time for both rigid and fluid cargoes, the fluid cargo indeed accumulates more motors. The spatial distribution of motors in Fig 4b further highlights the difference between rigid and fluid cargoes with the probability of finding motors near the microtubule on a fluid cargo being larger at late times than for rigid cargo. It is to be noted that the slight increase at late times for rigid cargoes is simply due to the higher contribution to the average of cargoes, that happen to have motors clustered together, surviving for longer times.

thumbnail
Fig 4. Surface fluidity increases the availability of motors for binding.

(a) Ensemble average of the number of bound motors, as a function of time (s). Error bars represent the standard error of the mean. (b) Spatial distribution of the kinesin motors on cargo, computed as the probability of locating motors as a function of the polar angle θ and time, t. θ is the polar angle of the motor, measured with respect to z axis perpendicular to the microtubule and passing through the center of the cargo. Data was obtained from the transport of cargoes with a total of N = 16 motors at [ATP] = 2 mM.

https://doi.org/10.1371/journal.pcbi.1010217.g004

To further test the dependence of motor accumulation on the relative values of the two time scales, we considered four other sets of physiologically relevant parameters of N and [ATP] that had different relative magnitudes of τoff and τbind (shown in Table A in S2 Appendix) As predicted, there is no statistically significant difference in the average number of bound motors as a function of diffusivity (see Fig 5a) for cases when τoffτbind (N = 4, [ATP] = 2 mM and 100 μM). The number of bound motors, however, increases as a function of diffusion constant (D) for cases where τbind < τoff (N = 16 [ATP] = 2 mM and 100 μM, N = 4 [ATP] = 4.9 μM) as expected. This supports our expectation that fluid cargoes can indeed accumulate higher number of bound motors than rigid cargoes, but only under conditions that ensure that τbind < τoff, i.e. when the number of motors (N) is high or the ATP concentration is low or both.

thumbnail
Fig 5. Fluid cargo have longer run lengths than rigid cargo.

(a) Average number bound motors, (b) Runlength as a function of diffusivity (D) for two different numbers of motors on the cargo (N) and three different ATP concentrations. 200 cargo runs were considered for each parameter set. Error bars represent the standard error of the mean (SEM). The average number of bound motors in (a) was calculated as the mean of the number of bound motors (n) in all time samples (with data sampling rate = 100 s−1) of these 200 cargo runs. Rotational diffusion of cargo was not considered in these simulations. The sample size was large and hence the error bars obtained as the standard error of the mean is very small. (c) Comparison between run lengths from our simulations (N = 4, [ATP] = 4.9 μM, yellow circles) and analytical estimates (maroon triangle, solid, dashed and dash-dotted lines) as described in the text. Left cartoon: A lipid cargo, with access area (Sa, green) and influx area (SI, yellow) shown. Influx area is defined as the region within from the access area.

https://doi.org/10.1371/journal.pcbi.1010217.g005

Fluid cargo have longer run lengths than rigid cargo

So far we have shown that increasing surface fluidity reduces inter-motor interference and promotes load-sharing leading to reduced off-rates. Fluidity also makes more motors available for binding in a cooperative fashion with increased binding leading to higher on-rates. We would expect that having more motors bind that survive longer should lead to increased run lengths for fluid cargo. The run lengths of cargoes measured as a function of diffusion constants (Fig 5b) suggest that this is indeed possible given the right ATP concentration and the number of motors on the cargo (see S13 Fig for additional diffusion constants). Although the run length is unaffected by fluidity when τbind is of the order of τoff (low N and medium/high [ATP]), the run length increases with the fluidity of the cargo surface when τoff > τbind (for high N or low [ATP] or both). The effect is much more pronounced when the disparity between the timescales is higher.

We saw earlier that the average number of bound motors shares similar trends (Fig 5a) as a function of cargo surface fluidity. It is to be noted, however, that cargoes can have the same average number of bound motors but different run lengths. For example, although high fluidity (D = 1 μm2 s−1) cargoes with N = 4 at [ATP] = 4.9 μM and N = 16 at [ATP] = 2 mM have the same average number of bound motors, the lower ATP case with fewer total motors shows a substantially larger run length (Fig 5a and 5b). This is a consequence of the higher tendency of a cargo to survive a 1 motor state at low ATP due to a lower motor detachment rate (more details in S4 Appendix).

For saturating ATP conditions, we note that the observation that a high number of motors (16 here) achieves a runlength of only about 2 μm is consistent with values predicted in other in vitro experiments and models [35]. This is because even though there are 16 motors on the cargo, the average number of motors that can access the microtubule is only about 1.37 (Na = 1 + (Sa/ST)). Our runlength values are also comparable to simulation results by Bovyn et. al. [36] for rigid and lipid cargoes under saturating ATP conditions. Finally, the range of runlength values for different numbers of motors for a cargo of radius 250 nm agrees with the range of values from experimental studies [44], lending further validity to our model.

While the results in Fig 5 are from simulations without considering cargo rotation, we also ran simulations with rotation and compared the runlength and average number of motors between these two cases S20 Fig. There is no noticeable change in these metrics with and without cargo rotation.

Overall, our results indicate that surface fluidity increases run lengths when the number of motors (N) is high or the ATP concentration is low or both. While previous studies have shown that the runlength of multi-kinesin rigid cargoes increases with decreasing ATP concentration [55], our simulations show explicitly that the increase in runlength with decreasing ATP concentration is significantly larger for lipid cargoes than for rigid cargoes.

Contributions of different effects to overall run length

We now consider the relative contribution to the overall run length of the different effects we showed are at play; the reduction in negative interference, the cooperative increase in on-rates and the increased availability of motors.

Analytical expression for run length.

To estimate the magnitude of these contributions, we first consider an analytical expression for the run length [23] (1) where vo, Na, πad, ϵm are the motor speed, the available number of motors for transport, the single motor binding rate and the mean motor off-rate respectively. Based on our results, we now generalize the expression by incorporating the dependence of these parameters on the cargo surface diffusivity D and the number of bound motors.

The mean off-rate, ϵm, is sensitive to the presence of inter-motor interference and the introduction of fluidity almost completely eliminates interference (Fig 2b–2f). This means we can treat ϵm as having two values; a higher one for rigid cargo due to interference and a lower one for fluid cargo (Fig 2f).

The number of available motors, Na, is usually taken to be the number of motors in the access area, Sa, for rigid cargo. For D > 0, however, additional motors from an influx area SI (see Fig 5c) can reach the access area before the cargo unbinds from the microtubule and contribute to transport. Here we estimate SI as defined by a region within a distance of the access area, representing the distance diffused by a motor over the mean motor lifetime . We, therefore, take Na to be the average number of motors in the access and influx areas combined, Na = 1 + (N − 1)(SI + Sa)/ST. Thus the number of available motors starts from a minimum value Na = 1 + (N − 1)(Sa/ST) at D = 0 and increases monotonically as function of D before saturating at Na = N.

Finally, motors bind with a rate πad = π0 = 5 s−1 when they are within the access area and do not bind otherwise. This means we can take the mean on-rate of the available motors to be πad = π0Sa/(SI + Sa). We note that the access area, Sa, depends on cargo geometry and in particular on the height of the cargo above the microtubule. As more motors bind, the height decreases, Sa increases and the effective on-rate increases (Fig 3). Finally, we note that the cargo velocity did not change by more than 5% over the range of parameters we tested (S15 Fig) and so we took vo to be constant.

Approximate bounds for run length.

Using these analytic estimates, we can compute approximate bounds for the relative contributions of the different effects to overall run length. As an example, we do this for the case of low motor density and low ATP (N = 4 [ATP] = 4.9 μM, Fig 5c, other cases shown in S17 Fig). We start with the limiting lower bound case where there is no additional availability of motors due to diffusion (Na = 1 + (N−1)(Sa/ST)), no reduction in interference (ϵm = 0.12 s−1, S16 Fig) and the cargo height is at its maximum (corresponding to having 1 motor bound (Fig 3a)). This yields the horizontal dashed line at the bottom in Fig 5c. A similar upper bound that assumes all motors are available (Na = N), interference is absent (ϵm = 0.1 s−1, see S16 Fig) and the cargo height corresponds to that at n = 2 (which is close to the average number of bound motors for this configuration), yields the horizontal dashed line at the top. While the bounds bracket the run lengths measured from the simulations, we notice that the lower limit is well below the actual values, especially at higher diffusion constants. It is to be noted that the bounds do not really extrapolate well to D = 0 because, for rigid cargo, the average run length is not set by the average number of motors in access area. A better approximation is the weighted average of run lengths for different numbers of motors in the access area weighted by the probability of having that many motors in the access area (maroon triangle in Fig 5c; details in S3 Appendix).

Comparison of contributions from different effects.

We first consider the contribution due to diffusion increasing the availability of motors alone. Allowing Na to increase with D but with no reduction in interference (ϵm = 0.12 s−1) and the maximum cargo height, corresponding to n = 1, yields the lowest dash-dotted curve in Fig 5c. While there is a substantial increase in run length with increasing D, this increase (of order Δa at high D) yields an estimate that is still significantly below the measured values at high D. This indicates that the increased availability of motors alone cannot account for the entire increase in run length. We next incorporate the cooperative increase in the accessible area due to more motors binding and the cargo moving closer to the microtubule by setting the cargo height to that corresponding to n = 2 bound motors. This yields an increase in the run length (middle dash-dotted curve in Fig 5c) but the gain, of order Δh at high D, is also too small to account for the observed increase. We finally incorporate the effect of reduced interference between the motors due to surface fluidity by taking ϵm = 0.1 s−1. This shifts the run lengths upwards considerably, by Δi at high D, and brings the estimates to within the range of the observed values (solid curve in Fig 5c; the upper dash-dotted curve is for a cargo height corresponding to one bound motor). Intriguingly, the surface fluidity based reduction in negative interference has an effect that is comparable or even stronger (see S17 Fig also) than the effect due to the increase in the number of available motors. We note that a reduction in the interference and the resultant reduction in mean off-rate leads to an increase in the mean motor lifetime, τm, and hence an increase in the influx area SI which is set by . Thus, reduction in interference contributes to an increased run length not only due to an increased lifetime of each motor but also due to the resultant increase in the number of motors available for transport. While we used the single motor lifetime to estimate SI, this represents a lower bound since new motors may reach the access area over the entire cargo runtime τc. Using the saturating value of τc (for high D) in place of in the expression for SI therefore gives an upper-bound estimate for the run length (upper dotted curve in Fig 5c). We note that the top three curves for estimates that include the effects of reduced interference and increased availability of motors reflect the measured values of run length fairly well as a function of increasing surface diffusion constant across different conditions of ATP and motor number (see S17 Fig).

Motor teams under external load

So far, we have considered the only load on the motor teams to come from viscous drag on the cargo which is very low for the aqueous medium we have been considering and usually dominated by thermal fluctuations. We now consider the behavior of the motor teams under external loads, specifically looking at motor clustering and collective force generation.

Motors dynamically cluster at specific locations on the lipid cargoes.

Past studies of force generation by motors coupled to lipid membranes have shown that the motor performance (as measured by the velocity of microtubules) is reduced due to “slippage” of motors on the lipid membrane [40]. Also, studies have shown that rigid-like ordered regions on the membrane are beneficial for transport of cargo because the motors can anchor on these ordered regions and use this as a rigid substrate to effectively generate forces [37, 38]. Hence, although motors on lipid cargo can interfere less and survive longer, a question arises as to whether the motors can efficiently generate forces against load while hauling cargo.

An important thing to note is that the geometry of the lipid membrane is a spherical surface, or at least a closed, bounded surface, for cargo in vivo such as vesicles and organelles. For such cargo, there is a limit to how far the motors can slip. Motors can slip only as long as the tangential force is non-zero. Once the tangential component reduces to zero and all the force is radial, a motor on a lipid cargo should be as efficient in generating force as on any other substrate.

To understand whether motors are likely to be in this radial force regime, we looked at the probability distribution of the location of motors on the cargo surface. To obtain this distribution, we ran simulations of cargo transport (N = 16, High [ATP] = 2mM) under different hindering loads on the cargo and recorded (sampling rate of 100 s−1) the anchor positions of MT bound motors relative to the body coordinate axis of the cargo (see schematic in Fig 6). Heat maps in Fig 6 show the probability distributions of the anchor positions of MT bound motors on the cargo surface. We can clearly see that there is a preferred position for motors on the cargo surface at different hindering loads (Fig 6). At zero hindering load, motors experience only the viscous load which is small enough that this distribution is broad and only a diffuse enhancement is visible near the MT (bottom of the cargo) However, as we increase the cargo load, motors increasingly appear at a preferred location on the cargo surface to balance the load. We calculated the location on the cargo where the tangential force reduces to zero and the horizontal component of force balances the hindering load and found that the height (z-coordinate) of this location above the MT increases with increasing hindering load (see S5 Appendix). The range of locations we obtained from this argument were in reasonable agreement with the values we measured from simulations (see S5 Appendix for calculation and values) validating the mechanism for preferential localization of motors. While we see motor density enrichment for both rigid and lipid cargoes with increasing hindering load, we note that the lipid membrane can allow motors to slide to this preferred location while this is not possible in rigid cargo, explaining the higher density enhancement in the lipid case. Motors can therefore dynamically cluster at a preferred location on the lipid cargo which minimizes the tangential component of motor forces and leads to better load-sharing. Our results are conisitent with prior observations of dynein clustering in axonemal endosomes and related modeling [56]. This sort of dynamic clustering can exist in addition to any intrinsic clustering mechanisms that may exist for in vivo lipid membranes [3638] which can enhance processivity in motor teams.

thumbnail
Fig 6. Motors dynamically cluster at specific locations on the lipid cargoes.

Average distribution of bound motors on the cargo surface at different hindering loads on the cargo (fh). Only multiple bound motor cases were considered. Data is from 200 cargo runs with N = 16, high [ATP] = 2 mM with anchor positions of bound motors recorded at a sampling rate of 100 s−1. Anchor positions are measured relative to the body coordinate axis of the cargo (schematic at top-right). Diffusion constant of motors on lipid cargoes was 1 μm2 s−1.

https://doi.org/10.1371/journal.pcbi.1010217.g006

This dynamic clustering of motors may also lead to the deformation of the vesicle. Collective forces due to kinesin motors have been observed to lead to the extraction of membrane tubes out of vesicles [43, 57]. There is a threshold of point force on membrane below which such a membrane extraction is not possible (see S1 Appendix). When the hindering load is low, motors are not localized, the force is low and not concentrated and the threshold is not reached, which is the reason we can ignore deformation of vesicles under normal cargo transport conditions (see S1 Appendix). If the hindering load exceeds the threshold, the clustering of motors may lead to membrane deformation and tube extraction. It should be noted that in our model, motors can come arbitrarily close to each other, while in reality, volume exclusion could broaden the localization region slightly, thereby affecting the threshold for tube extraction.

Motor teams in fluid cargoes generate higher collective forces.

We have predicted that lipid cargoes may have a higher number of MT bound motors which exhibit less negative interference, better load sharing and increased processivity. We have also predicted that motors can cluster together on lipid cargo under constant increased hindering load. Do all these factors imply that lipid membrane coupled motors can also generate more collective force?

Collective forces generated by teams on cargo have been measured by optical trap experiments [58] for motor-teams in vitro as well as in vivo. In vitro experiments have indicated that the kinesin-motor teams generate sub-additive forces [30, 59, 60] whereas in vivo experiments [60, 61] have shown that motors exert additive forces. While this difference in cooperativity might arise from differences in vitro and in vivo environments, we wanted to check if it could arise from the difference in the cargo surface properties in in vitro (rigid cargoes) and in vivo experiments (lipid cargoes).

To this end, we modeled cargo placement in an optical trap of force constant (0.06 pN/nm) [61] (see cartoon in Fig 7a) and ran simulations of cargo transport starting from a single microtubule bound motor for a maximum of 10 s (or until all motors unbind) for a total of 200 cargoes. We chose the parameters N = 4, [ATP] = 2 mM, which we believe is accessible to experiments [44]. We measured the values of the average total x-component of force due to motors on the cargo. We then rescaled the collective force [60] by N × Fs to get the rescaled average force per motor. The rescaled average force tells us how well the individual forces add up and whether or not the motors detach before reaching their stall force, serving as a good measure of the cooperativity between motors.

thumbnail
Fig 7. Motor teams in fluid cargoes generate higher collective forces.

(a) Schematic of the optical trap set up that we used in simulation. (b) Rescaled average collective force as a function of unloaded motor velocity (v0) and stall force of motor (Fs) in rigid (left) and lipid (right) cargo. (c) Ratio of between lipid and rigid cargo. 200 cargo runs were used for each parameter set. For each cargo run, data of motor forces for each microtubule bound motors were recorded at a sampling rate of 100/s for a maximum time of 10 s for each cargo run. Average collective force is calculated as the average of the total x-component of motor forces over all time samples and cargo runs and then rescaled by NFs. N = 4 and [ATP] = 2 mM.

https://doi.org/10.1371/journal.pcbi.1010217.g007

We scanned through the parameter space of single motor unloaded velocity (vo) and stall force (Fs) (see Fig 7b) to get a clear picture of motor cooperativity in different parameter regimes and to allow for comparisons with previous studies [60]. For both rigid and lipid cargoes, we observe that the rescaled average force is higher for lower stall forces and higher velocity in agreement with previous coarse-grained computation for multi-motor teams [60]. However, we find that motors on lipid cargo seem to generate higher rescaled average forces compared to those on rigid cargo virtually over the entire phase space explored (see Fig 7c). Interestingly, the rescaled average force that we measure in the lipid cargo case is comparable to previous estimates that assumes perfect load sharing [60]. This is again an indicator of improved load sharing properties in lipid cargoes.

Thus the observations of sub-additive forces in in vitro experiments [30, 59] and additive forces in in vivo [61] experiments could be due to differences in the fluidity of the cargo surface. More controlled experiments on membrane-bound and membrane-free cargoes [44] could verify our prediction.

Conclusion

Our results show that motors coupled through a rigid cargo experience a broader range of forces when multiple motors are bound, reflecting the geometric constraints on the motor-cargo attachment points leading to antagonistic forces and therefore negative mechanical interference between motors. This interference is consistent with literature findings [2830] and counteracts the expected load-sharing with increasing numbers of motors. We predict that the increased fluidity of the lipid membrane which allows the attachment points to move, suppresses strain, decreases the negative interference between kinesins and permits load-sharing between motors to be more effective. The reduction of negative interference not only contributes almost equally to run length at low ATP (and/or high motor numbers) as increased motor availability but is, in fact, the dominant effect at saturating ATP and low motor numbers (see S17(a) Fig) which is typically the normal physiological (and in vitro experimental) state. Direct measurements of this enhanced load-sharing may be possible using optical trap experiments with an applied load on rigid and fluid cargo. Tuning the load and motor number to a regime where enhanced load-sharing can decrease motor off-rates may be able to resolve the effects in comparisons of rigid versus fluid cargo.

The effects of inter-motor interference in rigid cargo are also expected to be seen in the measurements of speed. However, while our measurements of speeds from the simulations do show a decrease with increasing numbers of bound motors (S15 Fig), it is not as dramatic as reported in [44]. Consistent with [44], we find this decrease to be negated in the presence of fluidity but the effect is not statistically significant. One possibility is that the force-velocity relation we used (Eq 17 in Materials and Methods) is not sensitive enough at low loads. Another more intriguing possibility is that the experimentally observed speed-up of lipid cargo [44, 45] may be due to a dispersion of speeds among motors leading to a bias for faster motors in the lead, and slower motors trailing and unbinding more often. On the other hand, microtubule-assay experiments [40] have indicated that the gliding velcoity of microtubules is decreased when motors were anchored on a membrane because the motors “slip” on the membrane. This seems to contradict velocity measurements from both the membrane-bound cargo transport experiments [44, 45] and our simulation. However, this disagreement can be attributed to the difference in the geometry of lipid bilayer between assay set-up [40] and our cargo transport set-up. In a spherical geometry, like our cargo, motors slip only until the tangential component of force goes to zero, while, on a flat membrane, a motor slips every time it makes a step.

Since the reduction in interference manifests as a decrease in the off-rate of motors, it leads to increasing numbers of bound motors at higher membrane diffusion constants. Interestingly, we found that the on-rate was not constant, but increased with an increasing number of bound motors that effectively pulled the cargo surface closer to the microtubule. Monitoring the binding times of successive motors in an optical trap geometry could potentially allow for the experimental verification of this effect. This interesting cooperative effect on the on-rate that we uncovered has a smaller effect, in general, on the run length than the reduced interference or increased availability of motors but is still significant (≈10%), especially at low ATP and/or high motor numbers (see S17 Fig). While the cooperative increase in on-rates works for both rigid and fluid cargo, the effect is enhanced when more motors are available due to fluidity. Finally, we showed that increased fluidity results in an expected increased recruitment of motors to the microtubule, thereby increasing the number of bound motors, with the effect only being significant at very high motor density or low ATP. We found that the confluence of decreased interference, decreased off-rate, increased on-rate and increased motor recruitment can lead to positive feedback resulting in dramatic effects on the run length, under the right conditions. We predict, for example, that run length is insensitive to cargo surface fluidity for moderate numbers of motors at saturating ATP but can increase by several fold for fluid membranes at low ATP conditions. Taken together, our work reconciles varying experimental results including the observed insensitivity of run length to cargo fluidity at moderate motor densities and high ATP concentrations [44, 45]. Our prediction of a significant enhancement of run length at low ATP is potentially physiologically important as an adaptive mechanism, in vivo, under stress or starvation conditions [62].

We developed a generalized version of the analytical expression for run length that accounts for cargo surface fluidity and expect it to be useful to explore a wide range of parameter spaces for different motor types and cargo geometries going beyond our simulations, which were mostly done with a fixed cargo radius of 250 nm. Cargo geometry enters into the effective on-rate for motors which scales as πad = π0Sa/(SI + Sa). The increase in access area, Sa, as the cargo is pulled closer to the MT is determined by the curvature of the cargo surface and sets the magnitude of the cooperative increase in motor numbers. Even considering a fixed cargo height and high D, the on-rate still depends on the ratio of access area to the total surface area which decreases with increasing radius of curvature of the cargo. Thus, from our generalized version of Eq 1, we expect an increased on-rate and higher run lengths for smaller cargo with a fixed number of motors (see S18 Fig). This is consistent with our simulation results (S21 Fig) and with other numerical studies [36]. The smallest cargos, with sizes comparable to synaptic vesicles of about 100 nm [63], in fact, show an enormous increase in the run length over a range of fluidity that is consistent with the ranges for in vivo membranes in different contexts [40, 64]. For synaptic vesicles with fluid membranes therefore, transport over very large axonal distances, is a natural outcome of our model. If the surface density of motors is fixed, on the other hand, the quadratic increase in the number of available motors compensates for the reduction in on-rate resulting in larger run lengths for larger cargo (S19 Fig). In general, the changes in run length over physiological ranges of fluidity imply that surface fluidity could be used as a control parameter to regulate transport.

We also explored the effects of load on motor team performance and showed that motors on lipid cargo can undergo dynamic clustering and enhanced collective force generation. In vivo, increased load may be due to trapping or steric occlusion of the cargo, in which case motors can potentially cluster together exerting increased directed force against the load, perhaps even freeing the cargo. Such an adaptive response involving motor clustering has been observed in axonemal endosomes [56] and increased collective force generation in response to increased load could potentially be important in maintaining transport in vivo where crowding is high.

Our results on the transport properties of fluid cargo over a wide range of surface fluidity, ATP concentration, motor number and cargo size offer quantitative guidelines for the design of artificial cargo driven by kinesin teams. While all our results are specifically using parameters for kinesin motors, they can be readily generalized to other systems with different geometries of cargos with different motors coupled by a fluid surface. Thus our results could be of general value in designing artificial cargo utilizing teams of engineered versions of non-cooperative motors with fluid coupling to enhance load-sharing and processivity. More generally, coupling intrinsically non-cooperative mechanical force-generating elements via a fluid surface could be a generic mechanism to promote load-sharing that is exploited within the cell.

Materials and methods

We developed a Brownian dynamics model to address the question of how cargo surface fluidity influences transport by teams of kinesin motors. Early models of multi-motor cargo transport used a mean-field approach where all motors were considered to share the load equally [65]. Later stochastic models that assume unequal load sharing [30, 43, 46, 52, 6668] were found to describe experimental observations better [46, 69, 70]. The next improvement to the model was to explicitly consider that motors are bound to a three dimensional spherical cargo surface, that was typically considered rigid [51, 71]. Only a few recent models have attempted to consider the lipid membrane on the cargo surface [36, 45, 54, 56]. For example, Lombardo et. al. [54], in the context of transport of vesicles by teams of myosins, modeled an ideally fluid cargo surface where motors instantaneously relax the tangential component of force. Our model is closer to that of Bovyn et. al. [36] where we consider unequal load sharing with explicit implementation of motor diffusion on the cargo surface.

In our model, we consider a spherical cargo of radius, R (R = 250 nm in most cases) which is decorated with a given number of molecular motors (N) on the surface. N doesn’t change during the cargo run, meaning there is no binding and unbinding of motors between the cargo surface and the solution. Each of these N motors is initially assigned a random, uniformly distributed anchor position on the cargo surface. Molecular motors on rigid cargo are fixed on the cargo surface. So a given initial configuration of motors on the rigid cargo surface persists throughout the cargo run whereas motors on lipid cargo diffuse on the surface.

In this study we considered that all N motors are kinesin motors which have a rest length of Lmot = 57 nm [53]. We assume that an unbound kinesin motor binds to the microtubule with constant rate, π0 = 5 s−1 [43], if some part of the microtubule is within Lmot distance from the anchor position of that unbound motor. In other words, unbound motors bind with a specific rate to the microtubule if they can access the microtubule. The number of motors on the cargo, (N) is assumed to be constant, similar to several other modeling studies [36, 72]. In other words, we assume the timescale for motor detachment from the cargo is much larger than the lifetime of the cargo on the microtubule.

Each cargo run is initiated with at least one motor bound to microtubule and stopped when all the motors detach from the microtubule. This procedure is similar to Bovyn. et. al. [36] allowing a direct comparision between the runlength results between our studies. However, in other models [72] the cargo is tracked even during bulk diffusion when all motors are unbound till some motor binds back to the microtubule.

A microtubule bound motor is characterized by two position vectors, anchor point () on the cargo surface and head () on the microtubule. The head position is defined just by its continuous position along the x-axis, i.e. we do not model the lattice structure of the microtubule unlike in a recent Brownian dynamics model [73]. We justify this based on the fact that the microtubule has multiple protofilaments and when two kinesins are on different protofilaments they can have same x-position on the microtubule. There are also models that have explicitly incorporated the fact that the motor heads cannot step on each other [54, 73], which is likely of more importance when the filament has only a few tracks like actin, which has only two protofilaments in a helix, compared to 10–15 for microtubules.

A microtubule bound motor is assumed to exert a spring-like force when the length of motor exceeds the motor’s rest length (Lmot) with a force constant, kmot = 0.32 pN/nm [74, 75]. Let and . The force exerted by the motor on the cargo is then given by, (2) (3)

The translation velocity of the center of mass of cargo is given by the (overdamped) Langevin equation (4)

Here, is the force exerted on cargo by the jth motor. is the spring-like steric force on the cargo from the microtubule, represented with a high force constant 10kmot. We note that if the vesicle is deformable the steric spring constant could be significantly smaller. This force is present only if the cargo-microtubule distance is less than the sum of the cargo and microtubule radii. is the random force on the cargo due to collisions with the intracellular medium. We assume a normally distributed noise with zero mean, and the fluctuation-dissipation relation, 〈ζμ(t)ζν(t′)〉 = 2γckBμν δ(tt′). γc is the friction co-efficient for the cargo given by γc = 18πηv R, where ηv is the co-efficient of viscosity of cytoplasm experienced by cargo. We approximated ηv to be equal to the co-efficient of viscosity of water in room temperature (20°C), ηv = 10−3 Pa.s. We integrate Eq 4 using the Euler-Maruyama scheme (5) where and ξ1, ξ2, ξ3 are drawn from normal distribution with zero mean and unit variance.

In addition to translational motion, the cargo also has rotational dynamics due to thermal fluctuations and torques from motor forces. Consider a motor at position exerting a force on the cargo. The torque on the cargo due to this motor is . The total torque on the cargo due to all the motor forces is (6)

The angular displacement in time Δt taking into account this torque and thermal fluctuations is (7) Where, γR is the friction co-efficient, given by 8πηv R3. DR is the rotational diffusion constant of the cargo. α is a calibration constant to match the experimentally measured rotational mean square displacement. We set the rotational diffusion constant (DR) for lipid cargo to be equal to that of a free spherical bead in solution, which is kBT/γR = kBT/8πηv R3. The rotational diffusion constant for rigid cargo bound by one motor was measured [76] to be 7 × 10−2 rad2 s−1. In our simulations we used this value to calibrate the rotational diffusion of rigid cargo (see S1 Appendix). in Eq 7 is given by which is a random unit vector in 3-dimensions. To get this, we draw numbers a and b from uniform distribution in the interval [0, 1]. Calculate angles θt = cos−1(2a − 1) and ϕt = 2πb. Then .

Let Δθ be the magnitude and be the direction of this angular displacement vector . The Rodrigues’ rotation matrix corresponding to this rotation is [77] (8) where is the 3 × 3 identity matrix and is given by (9)

We update each anchor point position using this matrix (10)

At each time step we also update the anchor positions of each motor on the lipid cargo surface using a similar Brownian dynamics formalism given by (11) where Δlθ and Δlϕ are the small displacements along and directions in the plane tangential to cargo surface at . Δlθ and Δlϕ are given by (12) where ξa and ξb are random variables obtained from normal distribution with zero mean and unit variance. D is the diffusion constant for motor diffusion on cargo surface and γs is the friction coefficient given by γs = kBT/D. (Fθ, Fϕ) are the components of motor forces along and respectively which can be obtained from the motor force in Cartesian co-ordinates, as follows (13)

At every time step, each bound kinesin motor hydrolyses an ATP molecule with certain probability and attempts to move forward on the microtubule. This stepping probability is a function of the motor force and also the ATP concentration. We have adopted the following relation for the stepping probability [47] (14) where δ is the step size, the distance moved by motor after hydrolyzing one ATP molecule (δ = 8 nm [13, 78, 79]). v is the velocity of kinesin motor which is a function of the ATP concentration and motor force and is taken to be (15) where is the velocity of motor under no-load condition at a given ATP concentration. The ATP dependence of v0 is described by the Michaelis-Menten equation (16)

We considered the no-load velocity at saturated ATP concentration to be vmax = 800 nm/s [47, 48] and Km = 44 μM [80].

gives the force dependence of the velocity. In the hindering direction, we assume [46, 47] (17)

F is the magnitude of motor force, . Fs is the stall force, the value of force beyond which kinesin motor stops walking. We considered Fs = 7 pN [48, 81, 82] and w = 2 [69]. In the assistive direction, velocity is assumed to be independent of force magnitude, [47, 48].

Experimentally it is found that a kinesin motor is more likely to detach from the microtubule when one head is detached from the microtubule while trying to take a step than when both the heads are bound to the microtubule [80, 83, 84]. In our model we assume that a motor can detach only when it tries to take a step. At every time step we first check whether a bound motor tries to make a step with probability using Eq 14. If it tries to take a step, we check whether it detaches from the microtubule before completing the step using a microscopic off-rate whose value is calibrated based on the experimentally observed off-rate as a function of force, F, at saturating ATP concentrations. (18) is the probability to step forward in time step Δt at high ATP concentration of 2 mM.

For hindering forces, we used the following relationship between observed off-rate and magnitude of motor force F developed in [15, 85, 86] based on Kramer’s theory [87] and used in several studies [23, 46, 47] (19) ϵ0 is the off-rate under no load condition. We used ϵ0 = 0.79 s−1 [47, 48]. Fd is the detachment force. We approximated Fd to be equal to the stall force Fs. For assistive forces, the relationship between observed off-rate and magnitude of force is taken to be [47, 48] (20)

Our model could easily be extended to incorporate more detailed chemo-mechanical models of motors [72, 73].

The code for our model and data are available at https://github.com/nsarpangala/lipid-cargo-transport.

Supporting information

S1 Fig. Motor and cargo position and corresponding forces.

Variables and the simulation parameters same as Fig 2a.

https://doi.org/10.1371/journal.pcbi.1010217.s001

(EPS)

S2 Fig. The fraction of the force distribution with a magnitude of force f greater than or equal to 0.01Fs, as a function of the diffusion constant for the different number of bound motors n.

Error bars represent the standard error of the mean obtained by the bootstrap method. Fs is the stall force of kinesin. The fraction increases as a function of n for rigid cargo but decreases for lipid cargo. Required force distributions as a function of n and D were obtained using the same procedure as explained in Fig 2b (simulations without rotational diffusion, N = 16, [ATP] = 2 mM).

https://doi.org/10.1371/journal.pcbi.1010217.s002

(EPS)

S3 Fig. Fraction of force distribution with force magnitude greater than 0.01Fs in the hindering direction.

Required force distributions as a function of n and D were obtained using the same procedure as explained in Fig 2b (simulations without rotational diffusion, N = 16, [ATP] = 2 mM).

https://doi.org/10.1371/journal.pcbi.1010217.s003

(EPS)

S4 Fig. Fraction of force distribution with force magnitude greater than 0.01Fs in the assistive direction.

Required force distributions as a function of n and D were obtained using the same procedure as explained in Fig 2b (simulations without rotational diffusion, N = 16, [ATP] = 2 mM).

https://doi.org/10.1371/journal.pcbi.1010217.s004

(EPS)

S5 Fig. Mean value motor force (in Newtons) experienced by a bound motor as a function of the number of bound motors n and diffusion constant D.

Mean force magnitude is negative which means on an average the motor is in the hindering direction, as one would expect because of active motion of motor on MT. The magnitude of mean force decreases with an increase in the number of bound motors due to load sharing. (simulations without rotational diffusion, N = 16, [ATP] = 2mM).

https://doi.org/10.1371/journal.pcbi.1010217.s005

(EPS)

S6 Fig. The mean variance of forces (in units of N2) among bound motors.

Calculated as where the summation is over all the sample time points where the number of bound motors is n. is the variance in the magnitude of force experienced by the n bound motors at ith data sample. S = 10000 except for (i) D = 0 n = 3, S = 7360 (ii) D = 0 n = 4, S = 1083 (iii) D = 0.001 n = 4, S = 3754 (iv) D = 0.01 n = 4, S = 6396 (v) D = 0.01 n = 4, S = 3508 (vi) D = 0.1 n = 4, S = 4883. (simulations without rotational diffusion, N = 16 [ATP] = 2 mM).

https://doi.org/10.1371/journal.pcbi.1010217.s006

(EPS)

S7 Fig. Simulation results for different cargo radii.

(a) Distribution of the absolute magnitude of motor forces in rigid and lipid cargoes for two different cargo radii, without considering rotational diffusion of cargo. (b) Distribution of the absolute magnitude of motor forces in rigid and lipid cargoes for two different cargo radii, including rotational diffusion of cargo. (c) Mean force correlation between x-components of motor forces (d) Mean single motor off-rate. N = 16, [ATP] = 2 mM.

https://doi.org/10.1371/journal.pcbi.1010217.s007

(EPS)

S8 Fig. Distribution of cargo displacements (position fluctuations) measured in 0.01 s time intervals along (a) x-axis (b) y-axis and the (c) z-axis.

In general, fluctuations increase with the increase in the fluidity of the cargo surface. Vertical lines indicate the mean values of fluctuations, interpreted as veff δt where veff is the effective velocity of cargo and δt is the time interval used for measuring fluctuations (0.01 s). Fluctuations along the x-axis have a positive mean (veff δt ≈ 8nm) because of the active motion due to molecular motors. Fluctuations along the y-axis and z-axis have a 0 mean. Fluctuations along the z-axis are much narrower than fluctuations along the y-axis and x-axis because of the steric force due to the microtubule. Data were obtained from the transport of cargoes (without rotational diffusion) with a total of N = 16 motors at [ATP] = 2 mM. 200 cargo runs were considered for each diffusion constant.

https://doi.org/10.1371/journal.pcbi.1010217.s008

(EPS)

S9 Fig. Standard deviation (nm) of the position fluctuations as a function of diffusivity and the number of bound motors.

Data were obtained from the simulation of the transport of cargoes (without rotational diffusion) with a total of N = 16 motors at [ATP] = 2 mM. 200 cargo runs were considered for each diffusion constant.

https://doi.org/10.1371/journal.pcbi.1010217.s009

(EPS)

S10 Fig.

(a) Conditional on-rate measured from the simulation trajectories. We simulated 200 cargo runs (without rotational diffusion) for N = 16 [ATP] = 2 mM case and recorded the value of the number of bound motors at a sampling rate of 100 s−1. With this data, we identified all the n bound motors to n + 1 bound motor transitions, calculated the mean time for such transitions, and then the mean rate (1/meantime). We then divided this rate by the number of unbound motors Nn to obtain the conditional on-rate per motor. We then repeated the analysis process for a different number of bound motors n and diffusion constants D. (b) Comparing the measured conditional on-rate with analytical expressions. Black triangles with dashed lines provide conditional on-rate where we consider the increase in on-rate of a single motor with an increase in the number of bound motors. Grey circles with a dashed line provide conditional on rate assuming a constant single motor on-rate. It can be inferred that if we measure this quantity in experiments one can expect to see only a slight increase as the cargo comes closer to the microtubule. This is because as the number of bound motors increases, the effective detachment rate of the n-bound state increases, hence the gating time decreases. Our analysis shows that the decrease in gating time contributes more to the increase in the conditional on-rate than the increase in the single motor on-rate. See S2 Appendix for more details.

https://doi.org/10.1371/journal.pcbi.1010217.s010

(EPS)

S11 Fig. Ensemble average (± SEM) of the number of bound motors for three different cargo-motor systems.

Vertical lines indicate the estimated time-scales, mean time for a new motor to bind—τbind (dashed lines) and mean unbinding time of a kinesin motor—τoff (solid black line). N is the total number of motors on cargo, High [ATP] = 2 mM, Low [ATP] = 4.9 μM. The averaging was performed over 200 cargo runs (without rotational diffusion) each. Please refer to Table 1 in S2 Appendix for more information on τbind and τoff.

https://doi.org/10.1371/journal.pcbi.1010217.s011

(EPS)

S12 Fig. The distribution of the angular distance (ΔΘ) between the anchor positions of microtubule bound motors along the great circles connecting them (ΔΘ).

200 cargo runs were considered for each parameter set with data sampling rate = 100 s−1).

https://doi.org/10.1371/journal.pcbi.1010217.s012

(EPS)

S13 Fig. (a) Average number of motors (± SEM) (b) Runlength (± SEM) for more values of diffusion constants.

Averaging was performed over 200 cargo runs in each case.

https://doi.org/10.1371/journal.pcbi.1010217.s013

(EPS)

S14 Fig. Probability distribution of the number of bound motors in three different cargo-motor systems. High [ATP] = 2 mM, Low [ATP] = 4.9 μM.

Data were obtained from 200 cargo runs for each case.

https://doi.org/10.1371/journal.pcbi.1010217.s014

(EPS)

S15 Fig. Velocity of cargo along x-axis, as a function of the number of bound motors and diffusion constant.

Velocity is measured as the ratio of mean displacement along x-axis in a given time window (we took δt = 0.1s) to δt. Cargo position data was obtained from the simulations of the transport of cargoes with a total of N = 16 motors at [ATP] = 2 mM. 200 cargo runs were considered each for each diffusion constant. We recoreded data at a sampling rate of 100 s−1. Data in (b) is for N = 16, [ATP] = 2 mM without rotational diffusion. Error bars represent the standard error of the mean. For figure (a) we considered 10000 random time windows over 200 cargo runs to get the mean displacement for each diffusion constant. For figure (b) we considered all the time window samples with given n over all the cargo runs.

https://doi.org/10.1371/journal.pcbi.1010217.s015

(EPS)

S16 Fig. Average value of motor off-rate as a function of fluidity of cargo surface for N = 4 [ATP] = 4.9 μM.

(Sample size = 530 from 200 cargo runs).

https://doi.org/10.1371/journal.pcbi.1010217.s016

(EPS)

S17 Fig. Comparison between run lengths from our simulations (circles) and analytical estimates (maroon triangle, solid, dashed and dash-dotted lines) as described in the main text.

https://doi.org/10.1371/journal.pcbi.1010217.s017

(EPS)

S18 Fig. Analytically estimated runlengths for different cargo radii for a fixed number of motors on cargo.

The general method for calculating runlength analytically is described in the main text. We have numerically computed the access area, Sa, considering the typical distance between cargo surface and MT for 1 motor bound case. for computing influx area, , was taken to be 1 s which is the typical motor unbinding time at saturating ATP concentration. Yellow vertical bands in (a) and (b) correspond to the range of physiologically relevant diffusion constants of motors.

https://doi.org/10.1371/journal.pcbi.1010217.s018

(EPS)

S19 Fig. Analytically estimated runlengths as a function of cargo radius (R) for a fixed surface motor density (σ).

The general method for calculating runlength analytically is described in the main text. We have numerically computed access area, Sa, considering the typical distance between cargo surface and MT for 1 motor bound case. for computing influx area, , was taken to be 1 s which is the typical motor unbinding time at saturating ATP concentration. Yellow vertical bands in (a) and (b) correspond to the range of physiologically relevant diffusion constants of motors.

https://doi.org/10.1371/journal.pcbi.1010217.s019

(EPS)

S20 Fig. Comparison of the number of bound motors and cargo runlength without and with rotational diffusion in the model.

(a) Average number of bound motors (b) Runlength. Error bars represent the standard error of the mean. [ATP] = 2 mM.

https://doi.org/10.1371/journal.pcbi.1010217.s020

(EPS)

S21 Fig. Simulation results for different cargo radii.

Runlength, and average number of bound motors in rigid and lipid cargies as a function of the cargo radius. N = 16, [ATP] = 2 mM.

https://doi.org/10.1371/journal.pcbi.1010217.s021

(EPS)

S22 Fig. Simulation results for different motor stiffness (kmot).

Runlength, average number of motors and motor off-rate for different motor stiffness. N = 16 [ATP] = 2 mM. Darkgreen is lipid cargo (D = 1μm2 s−1) and grey is rigid cargo (D = 0).

https://doi.org/10.1371/journal.pcbi.1010217.s022

(EPS)

S23 Fig. The flow chart of the simulation.

The cargo is held near the microtubule until at least one motor is bound. Once a motor is bound, we loop over the series of steps shown in the middle until either all motors are unbound or maximum time is reached. At regular time intervals, we record relevant data like the cargo center of mass, anchor positions of motors on cargo surface, head positions of bound motors on the microtubule.

https://doi.org/10.1371/journal.pcbi.1010217.s023

(EPS)

S1 Video. Visualization of our cargo transport model.

The light brown colored sphere represents the cargo. Small red spheres are unbound motors. Unbound motors diffuse on the cargo surface and bind to the microtubule with a rate of 5 s−1 when they come closer than their rest length to the microtubule. Blue sticks represent bound motors that walk along the microtubule with a rate dependent on the force that they experience. Bound motors unbind with a force-dependent rate. The microtubule is represented as the dark-green horizontal line. Data for this movie was obtained from the simulation of a cargo with 16 motors. The data sampling rate was 105 s−1. We used Mayavi [88] for visualization.

https://doi.org/10.1371/journal.pcbi.1010217.s024

(MOV)

S2 Video. Cargo dynamics viewed at a shorter time and length scale.

The data for this movie is the same as S1 Video.

https://doi.org/10.1371/journal.pcbi.1010217.s025

(MOV)

S1 Appendix. Notes on deformation of vesicle due to motor forces and rotational diffusion of cargo.

https://doi.org/10.1371/journal.pcbi.1010217.s026

(PDF)

S2 Appendix. Estimation of time scales τbind and τoff and a note on conditional on-rate.

https://doi.org/10.1371/journal.pcbi.1010217.s027

(PDF)

S3 Appendix. Analytical estimation of run length of a rigid cargo.

https://doi.org/10.1371/journal.pcbi.1010217.s028

(PDF)

S4 Appendix. Discussion on how the runlength of cargo can be different even though the average number of motors are the same.

https://doi.org/10.1371/journal.pcbi.1010217.s029

(PDF)

S5 Appendix. Details of different force related metrics used in Fig 2.

https://doi.org/10.1371/journal.pcbi.1010217.s030

(PDF)

Acknowledgments

We are grateful for helpful discussions with David Quint and Jing Xu and insightful comments and suggestions from referees.

References

  1. 1. Howard J. Mechanics of motor proteins and the cytoskeleton. Sinauer Associates (OUP); 2001.
  2. 2. O’Connor CM, Adams JU. Essentials of cell biology. Cambridge, MA: NPG Education; 2010. Available from: https://www.nature.com/scitable/ebooks/essentials-of-cell-biology-14749010.
  3. 3. Ross JL, Ali MY, Warshaw DM. Cargo transport: molecular motors navigate a complex cytoskeleton. Curr Opin Cell Biol. 2008;20(1):41–107. pmid:18226515
  4. 4. Sabharwal V, Koushika SP. Crowd Control: Effects of Physical Crowding on Cargo Movement in Healthy and Diseased Neurons. Front Cell Neurosci. 2019;13:470. pmid:31708745
  5. 5. Lakadamyali M. Navigating the cell: how motors overcome roadblocks and traffic jams to efficiently transport cargo. Phys Chem Chem Phys. 2014;16:5907. pmid:24557020
  6. 6. Grundke-Iqbal I, Iqbal K, Tung YC, Quinlan M, Wisniewski HM, Binder LI. Abnormal phosphorylation of the microtubule-associated protein tau (tau) in Alzheimer cytoskeletal pathology. Proc Natl Acad Sci USA. 1986;83(13):4913–4917. pmid:3088567
  7. 7. Iqbal K, Liu F, Gong CX, Grundke-Iqbal I. Tau in Alzheimer Disease and Related Tauopathies. Curr Alzheimer Res. 2010;7(8):656–664. pmid:20678074
  8. 8. Chevalier-Larsen E, Holzbaur ELF. Axonal transport and neurodegenerative disease. Biochim Biophys Acta Mol Basis Dis. 2006;1762(11):1094—1108. pmid:16730956
  9. 9. Millecamps S, Julien JP. Axonal transport deficits and neurodegenerative diseases. Nat Rev Neurosci. 2013;14(3):161–176. pmid:23361386
  10. 10. Langford GM. Actin- and microtubule-dependent organelle motors: interrelationships between the two motility systems. Curr Opin Cell Biol. 1995;7(1):82–88. pmid:7755993
  11. 11. Svoboda K, Block SM. Force and velocity measured for single kinesin molecules. Cell. 1994;77(5):773–784. pmid:8205624
  12. 12. Visscher K, Schnitzer MJ, Block SM. Single kinesin molecules studied with a molecular force clamp. Nature. 1999;400(6740):184–189. pmid:10408448
  13. 13. Yildiz A, Tomishige M, Vale RD, Selvin PR. Kinesin Walks Hand-Over-Hand. Science. 2004;303(5658):676–678.
  14. 14. Andreasson JOL, Milic B, Chen GY, Guydosh NR, Hancock WO, Block SM. Examining kinesin processivity within a general gating framework. eLife. 2015;4. pmid:25902401
  15. 15. Schnitzer MJ, Visscher K, Block SM. Force production by single kinesin motors. Nat Cell Biol. 2000;2(10):718–723. pmid:11025662
  16. 16. Sellers JR, Veigel C. Walking with myosin V. Curr Opin Cell Biol. 2006;18(1):68—73. pmid:16378722
  17. 17. Block SM. Kinesin Motor Mechanics: Binding, Stepping, Tracking, Gating, and Limping. Biophys J. 2007;92(9):2986–2995. pmid:17325011
  18. 18. Ally S, Larson AG, Barlan K, Rice SE, Gelfand VI. Opposite-polarity motors activate one another to trigger cargo transport in live cells. J Cell Biol. 2009;187(7):1071–82. pmid:20038680
  19. 19. Holzbaur EL, Goldman YE. Coordination of molecular motors: from in vitro assays to intracellular dynamics. Curr Opin Cell Biol. 2010;22(1):4–13. pmid:20102789
  20. 20. Shubeita GT, Tran SL, Xu J, Vershinin M, Cermelli S, Cotton SL, et al. Consequences of Motor Copy Number on the Intracellular Transport of Kinesin-1-Driven Lipid Droplets. Cell. 2008;135(6):1098–1107. pmid:19070579
  21. 21. Mallik R, Gross SP. Molecular motors: strategies to get along. Curr Biol. 2004;14(22):R971–R982.
  22. 22. Klumpp S, Nieuwenhuizen TM, Lipowsky R. Movements of molecular motors: ratchets, random walks and traffic phenomena. Physica E. 2005;29(1–2):380–389.
  23. 23. Klumpp S, Lipowsky R. Cooperative Cargo Transport by Several Molecular Motors. Proc Natl Acad Sci USA. 2005;102(48):17284–17289.
  24. 24. Beeg J, Klumpp S, Dimova R, Garcia RS, Unger E. Transport of Beads by Several Kinesin Motors. Biophys J. 2008;94(2):532–541. pmid:17872957
  25. 25. Müller MJI, Klumpp S, Lipowsky R. Tug-of-war as a cooperative mechanism for bidirectional cargo transport by molecular motors. Proc Natl Acad Sci USA. 2008;105(12):4609–4614. pmid:18347340
  26. 26. Ando D, Mattson MK, Xu J, Gopinathan A. Cooperative protofilament switching emerges from inter-motor interference in multiple-motor transport. Sci Rep. 2014;4:7255. pmid:25434968
  27. 27. Huang KC, Vega C, Gopinathan A. Conformational changes, diffusion and collective behavior in monomeric kinesin-based motility. J Phys Condens Matter. 2011;23(37):374106. pmid:21862841
  28. 28. Driver JW, Rogers AR, Jamison DK, Das RK, Kolomeisky AB, Diehl MR. Coupling between motor proteins determines dynamic behaviors of motor protein assemblies. Phys Chem Chem Phys. 2010;12(35):10398. pmid:20582368
  29. 29. Rogers AR, Driver JW, Constantinou PE, Jamison DK, Diehl MR. Negative interference dominates collective transport of kinesin motors in the absence of load. Phys Chem Chem Phys. 2009;11(24):4800–4803.
  30. 30. Jamison DK, Driver JW, Rogers AR, Constantinou PE, Diehl MR. Two kinesins transport cargo primarily via the action of one motor: Implications for intracellular transport. Biophys J. 2010;99(9):2967–2977.
  31. 31. Peker I, Granek R. Multimotor Driven Cargos: From Single Motor under Load to the Role of Motor-Motor Coupling. J Phys Chem B. 2016;120(26):6319–6326.
  32. 32. McLaughlin RT, Diehl MR, Kolomeisky AB. Collective dynamics of processive cytoskeletal motors. Soft Matter. 2016;12(1):14–21. pmid:26444155
  33. 33. Maday S, Twelvetrees AE, Moughamian AJ, Holzbaur ELF. Axonal Transport: Cargo-Specific Mechanisms of Motility and Regulation. Neuron. 2014;84(2):292–309.
  34. 34. Barlan K, Rossow MJ, Gelfand VI. The journey of the organelle: teamwork and regulation in intracellular transport. Curr Opin Cell Biol. 2013;25(4):483–488.
  35. 35. Jiang R, Vandal S, Park S, Majd S, Tüzel E, Hancock WO. Microtubule binding kinetics of membrane-bound kinesin-1 predicts high motor copy numbers on intracellular cargo. Proc Nat Acad Sci USA. 2019;116(52):26564–26570.
  36. 36. Bovyn M, Janakaloti Narayanareddy BR, Gross S, Allard J. Diffusion of kinesin motors on cargo can enhance binding and run lengths during intracellular transport. Mol Biol Cell. 2021;32(9):984–994. pmid:33439674
  37. 37. Klopfenstein DR, Tomishige M, Stuurman N, Vale RD. Role of Phosphatidylinositol(4, 5)bisphosphate Organization in Membrane Transport by the Unc104 Kinesin Motor. Cell. 2002;109(3):347–358. pmid:12015984
  38. 38. Rai A, Pathak D, Thakur S, Singh S, Dubey AK, Mallik R. Dynein Clusters into Lipid Microdomains on Phagosomes to Drive Rapid Transport toward Lysosomes. Cell. 2016;164(4):722–734. pmid:26853472
  39. 39. Efremov AK, Radhakrishnan A, Tsao DS, Bookwalter CS, Trybus KM, Diehl MR. Delineating cooperative responses of processive motors in living cells. Proc Natl Acad Sci USA. 2014;111(3):E334–E343. pmid:24402168
  40. 40. Grover R, Fischer J, Schwarz FW, Walter WJ, Schwille P, Diez S. Transport efficiency of membrane-anchored kinesin-1 motors depends on motor density and diffusivity. Proc Natl Acad Sci USA. 2016;113(46):E7185–E7193. pmid:27803325
  41. 41. Lopes J, Quint DA, Chapman DE, Xu M, Gopinathan A, Hirst LS. Membrane mediated motor kinetics in microtubule gliding assays. Sci Rep. 2019;9(1). pmid:31270348
  42. 42. Memarian FL, Lopes JD, Schwarzendahl FJ, Athani MG, Sarpangala N, Gopinathan A, et al. Active nematic order and dynamic lane formation of microtubules driven by membrane-bound diffusing motors. Proceedings of the National Academy of Sciences. 2021;118(52). pmid:34934005
  43. 43. Leduc C, Campàs O, Zeldovich KB, Roux A, Jolimaitre P, Bourel-Bonnet L, et al. Cooperative extraction of membrane nanotubes by molecular motors. Proc Natl Acad Sci USA. 2004;101(49):17096–17101.
  44. 44. Li Q, Tseng KF, King SJ, Qiu W, Xu J. A fluid membrane enhances the velocity of cargo transport by small teams of kinesin-1. J Chem Phys. 2018;148(12):123318. pmid:29604873
  45. 45. Nelson SR, Trybus KM, Warshaw DM. Motor coupling through lipid membranes enhances transport velocities for ensembles of myosin Va. Proc Natl Acad Sci USA. 2014;111(38):E3986–E3995. pmid:25201964
  46. 46. Kunwar A, Mogilner A. Robust transport by multiple motors with nonlinear force– velocity relations and stochastic load sharing. Phys Biol. 2010;7(1):16012. pmid:20147778
  47. 47. Arpag G, Shastry S, Hancock WO, Tüzel E. Transport by populations of fast and slow kinesins uncovers novel family-dependent motor characteristics important for in vivo function. Biophys J. 2014;107(8):1896–1904. pmid:25418170
  48. 48. Andreasson JOL. Single-Molecule Biophysics of Kinesin Family Motor Proteins. Stanford University, Stanford,CA; 2013.
  49. 49. Berger F, Keller C, Klumpp S, Lipowsky R. Distinct Transport Regimes for Two Elastically Coupled Molecular Motors. Phys Rev Lett. 2012;108:208101. pmid:23003191
  50. 50. Gao Y, Anthony SM, Yu Y, Yi Y, Yu Y. Cargos Rotate at Microtubule Intersections during Intracellular Trafficking. Biophysical Journal. 2018;114(12):2900–2909.
  51. 51. Bergman JP, Bovyn MJ, Doval FF, Sharma A, Gudheti MV, Gross SP, et al. Cargo navigation across 3D microtubule intersections. Proc Natl Acad Sci USA. 2018;115(3):537–542. pmid:29295928
  52. 52. Wilson JO, Quint DA, Gopinathan A, Xu J. Cargo diffusion shortens single-kinesin runs at low viscous drag. Sci Rep. 2019;9(1):1–12.
  53. 53. Li Q, King SJ, Gopinathan A, Xu J. Quantitative Determination of the Probability of Multiple-Motor Transport in Bead-Based Assays. Biophys J. 2016;110(12):2720–2728.
  54. 54. Lombardo AT, Nelson SR, Ali MY, Kennedy GG, Trybus KM, Walcott S, et al. Myosin Va molecular motors manoeuvre liposome cargo through suspended actin filament intersections in vitro. Nature Communications. 2017;8(1).
  55. 55. Takshak A, Mishra N, Kulkarni A, Kunwar A. Effect of fuel concentration on cargo transport by a team of Kinesin motors. Physica A: Statistical Mechanics and its Applications. 2017;467:395–406.
  56. 56. Chowdary PD, Kaplan L, Che DL, Cui B. Dynamic Clustering of Dyneins on Axonal Endosomes: Evidence from High-Speed Darkfield Imaging. Biophysical Journal. 2018;115(2):230–241. pmid:29933888
  57. 57. Coordination of Kinesin Motors Pulling on Fluid Membranes. Biophysical Journal. 2008;94(12):5009–5017. pmid:18310242
  58. 58. Veigel C, Schmidt CF. Moving into the cell: single-molecule studies of molecular motors in complex environments. Nature Reviews Molecular Cell Biology. 2011;12(3):163–176. pmid:21326200
  59. 59. Furuta K, Furuta A, Toyoshima YY, Amino M, Oiwa K, Kojima H. Measuring collective transport by defined numbers of processive and nonprocessive kinesin motors. Proceedings of the National Academy of Sciences. 2013;110(2):501–506. pmid:23267076
  60. 60. Uçar MC, Lipowsky R. Collective Force Generation by Molecular Motors Is Determined by Strain-Induced Unbinding. Nano Letters. 2020;20(1):669–676. pmid:31797672
  61. 61. Leidel C, Longoria RA, Gutierrez FM, Shubeita GT. Measuring Molecular Motor Forces In Vivo: Implications for Tug-of-War Models of Bidirectional Transport. Biophysical Journal. 2012;103(3):492–500. pmid:22947865
  62. 62. Joyner RP, Tang JH, Helenius J, Dultz E, Brune C, Holt LJ, et al. A glucose-starvation response regulates the diffusion of macromolecules. eLife. 2016;5:e09376. pmid:27003290
  63. 63. Chen X, Williams JA. Exocytosis. In: Encyclopedia of Gastroenterology; 2004.
  64. 64. Klaus CJS, Raghunathan K, DiBenedetto E, Kenworthy AK. Analysis of diffusion in curved surfaces and its application to tubular membranes. Mol Biol Cell. 2016;27(24):3937–3946. pmid:27733625
  65. 65. Müller MJI, Klumpp S, Lipowsky R. Bidirectional Transport by Molecular Motors: Enhanced Processivity and Response to External Forces. Biophys J. 2010;98(11):2610–2618. pmid:20513405
  66. 66. Kunwar A, Vershinin M, Xu J, Gross SP. Stepping, Strain Gating, and an Unexpected Force-Velocity Curve for Multiple-Motor-Based Transport. Curr Biol. 2008;18(16):1173–1183. pmid:18701289
  67. 67. Bouzat S, Falo F. The influence of direct motor-motor interaction in models for cargo transport by a single team of motors. Phys Biol. 2010;7(4). pmid:21098940
  68. 68. Zhang Y. Cargo transport by several motors. Phys Rev E. 2011;83(1):011909. pmid:21405715
  69. 69. Kunwar A, Tripathy SK, Xu J, Mattson MK, Anand P, Sigua R, et al. Mechanical stochastic tug-of-war models cannot explain bidirectional lipid-droplet transport. Proc Natl Acad Sci USA. 2011;108(47):18960–18965. pmid:22084076
  70. 70. Bouzat S, Falo F. Tug of war of molecular motors: The effects of uneven load sharing. Phys Biol. 2011;8(6). pmid:22120209
  71. 71. Erickson RP, Jia Z, Gross SP, Yu CC. How molecular motors are arranged on a cargo is important for vesicular transport. PLoS Comput Biol. 2011;7(5):e1002032. pmid:21573204
  72. 72. Chen K, Nam W, Epureanu BI. Effects of neighboring microtubules on the microscopic dynamics of kinesin transport. Phys Rev E. 2018;98:052412.
  73. 73. Chen K, Nam W, Epureanu BI. Collective intracellular cargo transport by multiple kinesins on multiple microtubules. Phys Rev E. 2020;101:052413.
  74. 74. Coppin CM, Pierce DW, Hsu L, Vale RD. The load dependence of kinesin’s mechanical cycle. Proc Natl Acad Sci USA. 1997;94(16):8539–8544. pmid:9238012
  75. 75. Coppin CM, Finer JT, Spudich JA, Vale RD. Detection of sub-8-nm movements of kinesin by high-resolution optical-trap microscopy. Proc Natl Acad Sci USA. 1996;93(5):1913–1917. pmid:8700858
  76. 76. Gutiérrez-Medina B, Fehr AN, Block SM. Direct measurements of kinesin torsional properties reveal flexible domains and occasional stalk reversals during stepping. Proc Natl Acad Sci U S A. 2009;106(40):17007–17012. pmid:19805111
  77. 77. Serge B. Rodrigues’ Rotation Formula;. Available from: https://mathworld.wolfram.com/RodriguesRotationFormula.html.
  78. 78. Block SM, Schnitzer MJ. Kinesin hydrolyses one ATP per 8-nm step. Nature. 1997;388(6640):386–390. pmid:9237757
  79. 79. Coy DL, Wagenbach M, Howard J. Kinesin takes one 8-nm step for each ATP that it hydrolyzes. J Biol Chem. 1999;274(6):3667–3671. pmid:9920916
  80. 80. Xu J, Shu Z, King SJ, Gross SP. Tuning Multiple Motor Travel via Single Motor Velocity. Traffic. 2012;13(9):1198–1205. pmid:22672518
  81. 81. Carter NJ, Cross RA. Mechanics of the kinesin step. Nature. 2005;435(7040):308–312. pmid:15902249
  82. 82. Block SM, Asbury CL, Shaevitz JW, Lang MJ. Probing the kinesin reaction cycle with a 2D optical force clamp. Proc Natl Acad Sci USA. 2003;100(5):2351–2356. pmid:12591957
  83. 83. Yajima J, Alonso MC, Cross RA, Toyoshima YY. Direct Long-Term Observation of Kinesin Processivity at Low Load. Curr Biol. 2002;12(4):301–306.
  84. 84. Seitz A, Surrey T. Processive movement of single kinesins on crowded microtubules visualized using quantum dots. EMBO J. 2006;25(2):267–277. pmid:16407972
  85. 85. Bell GI. Models for the Specific Adhesion of Cells to Cells. Science. 1978;200(4342):618–627. pmid:347575
  86. 86. Evans E, Ritchie K. Dynamic strength of molecular adhesion bonds. Biophysical Journal. 1997;72(4):1541–1555. pmid:9083660
  87. 87. Kramers HA. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica. 1940;7(4):284–304.
  88. 88. Ramachandran P, Varoquaux G. Mayavi: 3D Visualization of Scientific Data. Comput Sci Eng. 2011;13(2):40–51.