Ndel1-derived peptides modulate bidirectional transport of injected beads in the squid giant axon

Long-distance intracellular transport is typically conducted by microtubule-associated molecular motors. Most members of the kinesin superfamily move toward the plus ends of microtubules, while the dynein complex is the main motor moving toward the minus ends of microtubules. Numerous vesicular and molecular cargoes are attached to these motors. Neurons, and in particular axons, which contain long, unipolar parallel microtubule arrays primarily orienting with the plus end towards the synapse, are very sensitive to impaired motor activity. Aberrant intracellular transport has been suggested as a common underlying theme for multiple neurodegenerative diseases such as Alzheimer's disease, ALS (amyotrophic lateral sclerosis), Parkinson disease and Huntington disease (Chevalier-Larsen and Holzbaur, 2006; Colin et al., 2008; Dompierre et al., 2007; Duncan and Goldstein, 2006; Holzbaur, 2004; Levy and Holzbaur, 2006; Perlson et al., 2009; Stokin and Goldstein, 2006). One of the major challenges is to characterize normal transport. However, this task is not simple. There is a growing realization that many cargoes bind multiple motors in vivo and show complex transport patterns involving movements in both anterograde and retrograde directions (Ally et al., 2009; Encalada et al., 2011; Guerin et al., 2010; Holzbaur and Goldman, 2010; Ma and Chisholm, 2002; Rogers et al., 1997; Scott et al., 2011; Shubeita et al., 2008).

Our understanding of the way in which bidirectional transport is successfully achieved in vivo is still rudimentary, although several models and experimental designs have been tested (Gross, 2004; Muller et al., 2008; Muller et al., 2010). How transport directionality is regulated temporally is even less understood. One direction looks to biochemical regulation, perhaps involving factors that affect the activities of kinesin and dynein motors, such as tau (Dixit et al., 2008; Mandelkow et al., 2004; Vershinin et al., 2007; Vershinin et al., 2008). Another direction suggests mechanical competition between sets of anterograde and retrograde directed motors (Guerin et al., 2010; Muller et al., 2008). One molecular complex that regulates different motors is dynactin; it contributes to dynein-based motility (Schroer, 2004; Vallee et al., 2004) but also regulates kinesins 2 and 5 (Berezuk and Schroer, 2007; Deacon et al., 2003; Gross, 2003; Haghnia et al., 2007; Martin et al., 1999). The LIS1-containing complex is also known to regulate dynein activity. LIS1 and its interacting partners Nde1 and Ndel1 directly affect the activity of cytoplasmic dynein (McKenney et al., 2010; Mesngon et al., 2006; Shmueli et al., 2010; Yamada et al., 2008). Nde1 interacts with LIS1, several centrosomal proteins, and dynein light and intermediate chains (Feng et al., 2000; Feng and Walsh, 2004; Hirohashi et al., 2006a; Hirohashi et al., 2006b; Stehman et al., 2007). Ndel1 is found in complex with LIS1, and with dynein heavy and intermediate chains (Niethammer et al., 2000; Sasaki et al., 2000). Thus, it is possible that Ndel1 is one of the biochemical factors that influence competition between forces generated by distinct molecular motor types.

We have tested the role of mammalian Ndel1 and specific peptides derived from it in the squid giant axon, which is an excellent model system in which to study neuronal transport (Allen et al., 1982; Bearer et al., 2000; Terada et al., 2000; Terasaki et al., 1995; Vale et al., 1985b). The axon contains plumes of organized microtubules (Bearer and Reese, 1999) and expresses members of both the kinesin superfamily (DeGiorgis et al., 2008; Vale et al., 1985a) and cytoplasmic dynein (Schnapp and Reese, 1989) molecular motors. Squid axoplasm has been used for in vitro assays (Brady et al., 1982; Brady et al., 1985; DeGiorgis et al., 2008; Kuznetsov et al., 1992; Vale et al., 1985a; Vale et al., 1985c), while the intact squid axon has been used as an in vivo system to follow transport of fluorescent proteins or negatively charged beads (Bearer et al., 2000; Terada et al., 2000; Terasaki et al., 1995). The motility of the population of negatively charged beads within injected axons has been shown to fit slow motility (Terada et al., 2000; Terasaki et al., 1995). A recent study has demonstrated that many cytosolic proteins travel in mammalian axons using a similar slow motor dependent transport (Scott et al., 2011). We injected negatively charged fluorescent beads into the squid giant axon together with mammalian Ndel1, Ndel1-derived peptides, or control proteins, and monitored their transport using time-lapse confocal microscopy as previous described (Bearer et al., 2000; Satpute-Krishnan et al., 2003; Satpute-Krishnan et al., 2006). Our analysis involved the use of automated tracking combined with large-scale, unbiased statistical analysis. The results confirm slow but processive motility of negatively charged beads and reveal distinct effects of the individual peptides on bead movements.

Short folded peptides derived from a protein are likely to reveal functional activities, which may be masked in full-length proteins due to changes in conformation. We chose our peptides based on their structure within the full-length protein. The target of our study, the Ndel1 protein was predicted to be unstructured using the FoldIndex prediction algorithm (http://bip.weizmann.ac.il/fldbin/findex) (Prilusky et al., 2005), with the exception of approximately fifty amino acids surrounding amino acid 250 (Fig. 1A). This domain is predicted to form a beta-sheet (http://www.ics.uci.edu/~baldig/betasheet.html). In addition, it contains the cysteine (C273) that can undergo lipid modification (Shmueli et al., 2010). We have demonstrated that palmitoylation of C273 negatively regulates the interaction between Ndel1 and cytoplasmic dynein (Shmueli et al., 2010). Other programs such as COILS (http://www.ch.embnet.org/software/COILS_form.html), and more importantly direct structural data for residues 6–166, indicate that the N-terminal part of Ndel1 forms a coiled-coil structure divided into three regions (Fig. 1B) (Derewenda et al., 2007). The dimerization domain of Ndel1 is composed of residues 8–99 (Derewenda et al., 2007). Residues 103–153 contain the minimal LIS1-binding domain within the N-terminal part of Ndell; deletion of twenty amino acids (114–133) was sufficient to abrogate the interaction with LIS1 (Yan et al., 2003). In this study we were interested in investigating LIS1-independent activities of Ndel1. Therefore, based on the above information we generated two peptides from Ndel1: DID, consisting of amino acids 4–103 and pep3, consisting of amino acids 238–284 (schematically shown in Fig. 1C). Based on the limited amino acid similarity between Nde1 and Ndel1 within the DID domain (supplementary material Fig. S1), we also generated one corresponding peptide from Nde1 DID amino acids 3–102 (Nde1-DID). All peptides were fused to gluthathione S-transferase (GST) for easy purification and solubility. GST was used as a control peptide in all experiments.

Protein interactions of these peptides were examined by pull-down experiments (Fig. 1D–E). Our results indicate that Ndel1 and the three peptides pulled down a relatively small amount of conventional kinesin in brain extract, nonetheless no interaction was observed with the control GST peptide (Fig. 1D), suggesting specificity. Ndel1, DID (derived from Ndel1) and the Nde1-DID peptide pulled down dynein intermediate chain (DIC), whereas no signal was noticed with either pep3 or the control GST peptide. Similar results were obtained using extracts of HEK293 or COS-7 cells. The interaction between Ndel1 with pep3 in a cell lysate was enhanced when Ndel1 was palmitoylated (Fig. 1E). Ndel1 and pep3 can interact directly as demonstrated by a GST pull down experiment, whereas GST-pep3 pulled down 6×His Ndel1 but the GST protein did not (Fig. 1F). Furthermore, pep3 can self-interact since GST-pep3 pulled down 6-His-pep3, but GST did not (Fig. 1G). The DID peptide pulled down detectable amounts of dynein heavy chain and dynein intermediate chain identified by mass-spectrometry sequencing of the indicated bands (Fig. 1H). Our results indicate that the DID peptide interacts with dynein (therefore named Dynein Interacting Domain) and pep3 is an additional Ndel1 self-association domain, which is sensitive to the palmitoylation status of Ndel1.

The squid genome is not yet available, however, it has been shown to contain multiple molecular motors including cytoplasmic dynein and members of the kinesin superfamily of proteins, which exhibit high degree of similarity to the mammalian ones (Degiorgis et al., 2011). Although the Ndel1 ortholog has not been detected in squid so far, the squid genome is likely to contain one based on Western blot analysis (Fig. 1I). Using polyclonal anti-Ndel1 antibodies we detected a cross-reactive band in a protein lysate from the squid axoplasm (Fig. 1I). This band has the same molecular weight as Ndel1 detected in a mouse brain protein lysate (Fig. 1I). Furthermore, our previous bioinformatic analysis revealed that Ndel1 and Nde1 are highly conserved proteins (Shmueli et al., 2010). Nevertheless, our interpretation of the results will include the possibility that the squid genome does not contain an ortholog of Ndel1.

We further examined whether Ndel1 and its related peptides may interact with conventional kinesin using far-Western analysis (Fig. 1J). The different tested proteins of kinesin included α2β2 that is the recombinant wild-type kinesin consisting of two heavy chains and two light chains, α2 is kinesin heavy chain homodimer and ΔTail is kinesin heavy chain deleted after residue M559. Our results indicate that the tested Ndel1/Nde1-derived peptides interacted with all the tested kinesin proteins whereas the GST control did not. We are not sure what are the functional implications of these interactions since NDEL1 did not activate kinesin in single-molecule assays.

To evaluate the role of Ndel1 on transport in vivo, we injected either GST, GST-Ndel1 (Ndel1), GST-Ndel1-DID (DID), the related peptide from Nde1 (Nde1-DID) and GST-pep3 (Pep3) mixed with 100 nm diameter carboxylated fluorescent beads into the giant axon of the squid (Fig. 2A–B). Transport of the fluorescent beads was monitored using confocal microscopy (supplementary Movies 1–3) (Fig. 2C–D). Movements of the fluorescent beads were analyzed using the automatic spots-tracking module of the Imaris software package (Bitplane, Inc.) to generate a set of position lists. The generated lists consisted of a track index, followed by x and y coordinates and a time stamp. The deduced tracks enabled to resolve the direction and velocity of each bead during each captured time frame. The anterograde direction was given a positive sign while the retrograde direction has a negative sign.

Each experiment was conducted in 4–6 different axons, recorded at 4 sec intervals for 100 frames. The motile beads were tracked resulting in 1194–4046 tracks. The tracks were separated into segments that were defined by continuous movements of a bead in one direction (4765–12324 segments per treatment). The movements from one frame to the next were used to calculate instantaneous movements (20,304–49,640 per treatment). The summary of the data is shown in Table 1. A few representative tracks are shown in Fig. 3A. Beads appear and disappear as they enter and leave the focal plane. Some beads moved in one direction only, along with pauses of zero velocity (the green, turquoise and purple traces, Fig. 3A). Other beads changed directions, either with or without intervening pauses (the red and blue traces, Fig. 3A). Thus several tracks include both anterograde and retrograde segments. We noticed a large variability of the velocity within the individual segments; beads were often accelerating or slowing down even while moving in one direction. To distinguish the effects of the various peptides on motor-induced bead motion a number of statistical tests to the data were applied. First, we compared the individual, instantaneous movements at the measurement interval from the four treatments in comparison to GST (Fig. 3B–E). The histograms have been normalized to the same areas, so that heights of the bars can be interpreted as probabilities. The horizontal axes show velocities in units of microns per second. The GST control appears in turquoise in all panels for comparison. The first and most notable feature is that there are no peaks at the canonical free motor velocities, which we would expect to be approximately 0.5∼1 µm/sec based on in vitro measurements (Howard et al., 1989; King and Schroer, 2000; Toba et al., 2006; Vale et al., 1996). Instead there is a tall central peak at zero, surrounded by a broad shoulder to the anterograde side and a slim shoulder to the retrograde side. These shoulders decay to a maximum velocity of approximately 1 µm/sec in both retrograde and anterograde directions, just where we might have expected to see peaks due to single motor activity. The shoulders are not symmetric, but show a skew toward the anterograde direction. This is consistent with the visual observation that the cloud of injected beads moves and spreads mainly toward the axon tip, away from the cell body. Overall, the instantaneous velocities are much smaller than expected from in vitro motility assays. This might indicate a drag against the motion, or the effect of competing motors. Subtle differences appear in the shoulders of the distributions with the various peptides. Specifically, the full-length Ndel1 protein increased instantaneous velocities slightly in the anterograde direction (Fig. 3B), while the DID and Nde1-DID peptides increased the instantaneous velocities in both directions with respect to the GST control (Fig. 3D,E). Pep3 had a more pronounced effect in comparison with that of full-length Ndel1, namely it increased the velocity of anterograde steps without an effect on the retrograde steps (Fig. 3C). Furthermore, addition of the peptides, but not of the full-length Ndel1, significantly reduced the number of immobile beads (Fig. 3C–D and Table 2). DID had the most pronounced effect in this regard and its addition decreased the percentage of paused beads from 58% to 38%. Addition of Pep3 or Nde1-DID resulted in 45% of paused beads in comparison with 58% with the GST control. As a quick test of significance we split each of the datasets into six random sub-groups and confirmed that the visual differences remained qualitatively unchanged among them. In order to confirm quantitatively that GST and treatment step velocity distributions are distinct, we applied the Kolmogorov-Smirnov test, which yielded extraordinarily small p values (<10⁻²¹).

The recorded tracks were further divided into retrograde and anterograde segments; when the movement of the bead stopped or changed direction a new segment was defined. To demonstrate the differences in the distributions of the segments more effectively we turned to a two-dimensional, color-based presentation (Fig. 4). For each treatment the segment duration (in seconds) is plotted along the x-axis, and the directional run length (in microns) along the y-axis. These plots show the correlation between persistence in length and temporal duration. In particular, each point in 2D shows the total length and duration of the segments it represents and the coloring represents the percentage of the segments with the particular duration and velocity. This presentation contains more information than conventional bar histograms as it shows how the segment lengths depend on their duration. We see clearly that anterograde segments are more numerous than retrograde in the GST control and in all treatments (about 2/3 of the total segments following exclusion of the very short segments, Fig. 4A left panel). In order to emphasize the differences, thin lines were drawn over the 2D plots representing bounds of 85% in the distributions of both length and duration of continuous runs (Fig. 4 left panels). A shift of the vertical line to the left indicates that a persistent motion is shorter on average; while the vertical position of the two horizontal lines indicates represents the segment lengths in either direction. The intersections of the horizontal and vertical lines with the y and the x axes indicate the duration of segments and the displacement of the segments and their confidence interval are shown in Table 3. The duration of segments in the anterograde direction was slightly affected by the treatments. It shortens in the presence of full length Ndel1 and extended a little with pep3 addition. All treatments decreased the duration of the retrograde segments, mainly Nde1-DID and DID. The displacement of the retrograde tracks was increased in all the treatments, where the most evident effects were noted following addition of DID and Nde1-DID.

In addition, we examined the average velocity per segment. The peak follows the lightest-colored ridges in Fig. 4 (right panels). The diagonal bounding lines are drawn to include 95% of the distributions (turquoise) and the median, 50% (grey). The median shows us the typical or most probable value, while the limit shows the average velocity of the fastest segments. All points that lie along a diagonal line starting from the origin have the same ratio of length to duration, so the slope represents a similar average velocity within the segments. The average velocities and their confidence interval were calculated at the 95% and 50% criteria (Table 3), and were normalized in relation to the retrograde or anterograde average velocities of control GST treated segments (Table 3). Overall, the effect of the peptides was more pronounced when the median velocities are inspected in comparison to the 95 percentile. All of the treatments increased the median velocities of the retrograde segments where DID and Nde1-DID exhibited a 4.5 fold increase in comparison with GST whereas the duration decreased (Table 3). In addition, all of the treatments resulted in increased velocities of the anterograde segments. A modest increase was noted when Ndel1 was added (13%), Pep3 and Nde1-DID each increased the anterograde velocities by 20% and DID contributed additional 40%.

How balanced bidirectional transport is generated in vivo is a topic of active investigation and several possible models have been proposed (Gross, 2004; Muller et al., 2008). There may be a mechanical competition that leads to temporary domination of one motor family over the other. In the “tug of war” model, for example, several motors of both polarities are available to engage the microtubule tracks. In balance the two types will cause a stall. If one side begins to detach, however, the fewer opposing motors remaining will require less and less force to detach one by one. This generates a detachment cascade that leads to one motor family dominating the motion at a time. Our data suggest a model in which oppositely directed motors remain engaged with the microtubules simultaneously. One side dominates the motion, while motors in the opposing direction confer a drag. This is consistent with the slow velocities we observe compared to in vitro studies. Our results fit well with previously measured velocities of negatively charged fluorescent beads in the squid axon. Previous publications measuring movements of beads in the squid axon have using manual tracking and noted that the average velocity was 0.02 µm/sec and the maximum velocity was 0.4 µm/sec (Terada et al., 2000; Terasaki et al., 1995). A recent article, which analyzed movements of PrP^c vesicles in mammalian hippocampal axons demonstrated that particles moving in the retrograde direction had shorter mean run lengths than those moving in the anterograde direction (4.8 versus 6.2 µm) (Encalada et al., 2011). Our results also demonstrate a reduced displacement for retrograde segments (3.7–6.3 versus 7.5–11 µm, Table 3). In addition, the kinetics of cargoes transported within an axon as well as the anterograde bias of this population in mammalian axons fitted the slow motility in our observed large scale tracking of individual beads (Scott et al., 2011). In the context of molecular motors, movement involves unit steps between structurally fixed binding sites along the microtubule. Free-running, the ATP turnover rate determines the maximal velocity as a stepping frequency. An opposing force, or drag, cannot change the step size, but rather affects strongly the time interval between steps, resulting in a reduced velocity. In our observations, the motor movements in both directions show large variation in velocity within a single track. Typically, the beads are not moving at constant velocity as would be expected from motor activity in vitro; rather they are accelerating or slowing down. Stall periods may represent a rigor state with both motors attached but unable to move. This is equivalent to a tug of war model in which the detachment cascade is not complete, perhaps a “kick and scream” variant to represent the competition before one side wins, as shown in Fig. 5. We will discuss below how Ndel1 related peptides affect bidirectional transport within this model.

In addition to force-based competition bidirectional transport may involve specific regulation of the motors themselves. Motors may be composed of different isoforms (King et al., 2002; Myers et al., 2007), they may be modified following translation (Song et al., 2007; Vaughan et al., 2001) and the interactions between the different motors and other regulatory molecules may vary (Ligon et al., 2004). The motor modifiers may have a direct effect on ATPase activity and force generation. Alternatively, they may affect the strength of motor interaction with microtubules, or the interactions between the motors and the various cargoes. Post-translational modifications of tubulin or binding of other microtubule associated proteins (MAPs), such as tau, may also affect the binding of motors to the tracks in a differential manner (Dixit et al., 2008; Dompierre et al., 2007; Erck et al., 2005; Ikegami et al., 2007; Mandelkow et al., 2004; Peris et al., 2006; Reed et al., 2006; Seitz et al., 2002). A few molecules are known to interact with different motor families and to affect their activity (Gross et al., 2002; Kardon and Vale, 2009). One such example is dynactin, an integral component of the dynein motor that allows it to bind a variety of cargoes, regulates its activity directly, and enhances its processivity (Schroer, 2004; Vallee et al., 2004). In addition, dynactin interacts with kinesin-2 and kinesin-5, two plus end-directed microtubule motors (Berezuk and Schroer, 2007; Blangy et al., 1997; Deacon et al., 2003). Dynactin coordinates between dynein and kinesin activities in several systems (Deacon et al., 2003; Gross, 2003; Haghnia et al., 2007; Martin et al., 1999).

Ndel1 and Nde1 are best known for their interaction with cytoplasmic dynein and LIS1 (Feng et al., 2000; Niethammer et al., 2000; Sasaki et al., 2000), through which they mediate their effects on dynein and retrograde transport (Feng and Walsh, 2004; Liang et al., 2004; McKenney et al., 2010; Shim et al., 2008; Stehman et al., 2007; Yamada et al., 2008; Yan et al., 2003; Zhang et al., 2009). We identified dynein-binding activity in the N-terminal DID domain and the addition of this peptide enhanced retrograde movements but the duration decreased. This activity may require dimerization properties embedded within the domain that were defined in structural and interaction-based studies (Bradshaw et al., 2009; Derewenda et al., 2007). An independent recent study demonstrated that the first eighty amino acids of Ndel1 (included in DID) binds directly to dynein intermediate chain and is essential for spindle pole organization in Xenopus egg extracts (Wang and Zheng, 2010). An additional recent study also demonstrated dynein binding activity in the first ninety-nine amino acids of Ndel1 and also Nde1 and further demonstrated that the interaction with dynein is regulated by phosphorylation in the C-terminus (Zylkiewicz et al., 2011). Therefore, even if the squid genome does not contain an ortholog of Ndel1 (which we believe is unlikely), the introduced DID peptide is likely to bind to the squid cytoplasmic dynein molecular motor which is similar to the mammalian one (Degiorgis et al., 2011). The second domain studied here, pep3 (amino acids 256–291), exposed a new dimerization domain that is sensitive to protein palmitoylation. If the squid does not contain an ortholog of Ndel1 we must speculate that the observed effects are due to some unknown protein-protein interactions. In addition, our results indicated that all the peptides used in this study are capable of forming a complex with conventional kinesin though the functional results of these interactions are not clear.

Treatment with Ndel1-derived peptides may affect the competition between dynein and kinesin in a number of ways: they may affect one or the other motor directly, they may affect the strength of the motor interaction with microtubules, or they may affect the strength of the motor interaction with the bead cargo. A telling observation lies in the position traces of Fig. 3A. Some tracks show movements in one direction only, along with pauses of zero velocity. Other traces show movements that pass smoothly through zero velocity, i.e., they slow down and reverse direction continuously without an intermediate stall. This observation is hard to reconcile with a model explaining bidirectional transport by strict alternation of kinesin and dynein motor activity. It would be a natural outcome of competition in which both motors remain engaged, at least transiently, during the excursions in either direction. Progressively increasing velocity would be a result of progressive detachment of opposing motors, while attachment of opposing motors would decrease the velocity. It is significant in this regard that those reversals in direction may or may not involve an intermediate stall.

Surprisingly, our results show that Ndel1 related peptides enhance primarily the anterograde movements of the negatively charged beads. The major influence on anterograde movement was unanticipated given that Ndel1 interacts directly with dynein, yet the anterograde enhancement appears in many aspects of the motion. When looking at the distribution of instantaneous movements we see a relative increase in the number of anterograde steps. This increase comes primarily at the expense of stall periods; the length distribution of individual retrograde steps weakly affected if at all. Furthermore, we see an enhancement of anterograde persistence. By bounding the segment distributions at 85% of each population, we see an increase in overall length of the anterograde movements under DID and Pep3 treatments. Interestingly, the full-length Ndel1 protein did not have a strong effect in most of the examined aspects. We assume that this is due to the fact that Ndel1 is overall an unstructured protein and its activities are exposed when it is found in the appropriate protein complexes. Dynamic post-translational modifications such as phosphorylation and palmitoylation are likely to play an important role in determining which protein complexes Ndel1 will participate. The retrograde segments were affected in a different way. All treatments decreased the duration of movements, while the displacement increased most following addition of DID and Nde1. Finally, there was an increase in both median and maximal average velocities per segment following all treatments in the retrograde direction, and in the median in the anterograde direction. This finding points to enhancement of some parameters relating to both motors. The most dramatic effects are observed in the velocities of the segments in the retrograde direction. The average control velocity is very slow and following the treatments increased 3–4.5 fold. This may hint to a distinct mechanism of action of the added peptides on the different motor proteins.

In this case how the Ndel1-derived peptides affect bidirectional transport may be interpreted that an inhibition of dynein attachment or stroke would appear as an enhancement of anterograde movement. We distinguish between the effects of pep3 and DID. Our data suggests that pep3 may affect anterograde and retrograde movements indirectly via the interaction of the pep3 peptide with Ndel1 (Fig. 5G–I). As mentioned above, in case that the squid does not contain an ortholog of Ndel1, the effects will be mediated through other, currently unknown protein interactions. The DID peptide has the capability to interact both with Ndel1 and directly with cytoplasmic dynein (Fig. 5J–L). As such, it may possibly enhance dynein's activity or affect the number or conformation of engaged retrograde motors. Accumulating data suggest that changes in the number of actively engaged motors directly affect the motility and directionality of different cargoes (Ally et al., 2009; Encalada et al., 2011; Hendricks et al., 2010; Shubeita et al., 2008). Nevertheless, the main consequences following reduction of the levels of either dynein or kinesin were decreases in the length of bidirectional runs, and increases in pause frequencies, consistent with down-regulation of opposing motor activity. Surprisingly however, mean segmental velocity, a parameter also influenced by motor activity, was largely unaffected (Encalada et al., 2011). Our study also suggests that active bidirectional motors are moving one cargo where each movement results from the net competition between the motors, which was displayed as a modified version of the “tug of war”, namely, “kick and scream”. Furthermore, our study demonstrated how additional regulatory proteins may affect retrograde and anterograde transport in a different way than directly affecting the number of motors. Addition of the peptides decreased the number of paused beads and increased the average velocities both in the anterograde and retrograde direction. Run length in the anterograde direction was largely unaffected, whereas in the retrograde direction it was shortened. Thus, our study shows an in vivo demonstration how bidirectional transport can be regulated using motor-binding proteins.

The back stellar nerve or giant axons were dissected from squid (Loligo pealleii) in running sea water. The axon was moved to calcium free seawater, surrounding fibers and connective tissues were removed using a stereomicroscope. Injections were assisted with a micromanipulator (Narishige SM20, Model SM-20) using pulled glass pipettes filled with a mercury drop attached to a Micrometer syringe (Gilmont S-1200, GS-1200). FluoSpheres® Fluorescent Microspheres (Molecular probes, Invitrogen, Carlsbad, CA) carboxylate-modified, 0.1 µm, red (580/605) were washed using Microcon YM-30 filter tubes (Millipore, Billerica, MA) with double-distilled water by centrifugation at 13,000 g, for 3 minutes at room temperature. The fluorescent beads were retrieved by inverting the filter followed by centrifugation at 1000 g for 3 minutes at room temperature. Double-distilled water was added to restore the initial volume, the protein was diluted in 1×PBS. Beads were mixed 1:1 with the recombinant proteins up to thirty minutes prior to the injection (3.85×10⁻⁶ % of the carboxylated binding sites on the beads are occupied by the recombinant protein, based on the assumption that one binding site can bind one protein molecule). All recombinant proteins and peptides were tagged with GST and identical molar quantities were taken for injection. Each beads:protein sample was loaded onto the injection pipette between the two oil droplets (Dimethylpolysiloxane, Sigma-Aldrich, St. Louis, MO) which were used for marking the injection point, monitoring the injection stability and preventing the penetration of sea water into the axoplasm. The axon and the injection pipette were placed on a stage of a regular up right light microscope (Axiscope, Carl Zeiss Inc., Thornwood, NY), the tip of the injection pipette was placed in close vicinity to the membrane of the axon. The pipette is then inserted into the axoplasm and the oil and the sample were injected into the middle of the axoplasm. The injection occurred 12–15 minutes after axon retrieval, after which the axon was placed onto a cover slip suitable for confocal microscopy imaging.

Imaging was done on a LSM 710 confocal microscope (Carl Zeiss Inc., Thornwood, NY) using channels; DIC (differential interference contrast) and a fluorescent channel (excitation 594 nm, 3% intensity, 597–677 nm filters). The injection spot was located using a 10× objective. Fluorescent beads were imaged using a 40× objective (LD C-Apochromat 40×/1.1 W korr, pixel size x, y 0.131203 µm), where frames were captured every 4 seconds one hundred times similar to previously described conditions (Satpute-Krishnan et al., 2006) and re-analyzed (Vermot et al., 2008). The first time-lapse series from each axon was analyzed.

Movements of the axon over time were corrected using the AutoAligner software (Bitplane Inc., Zurich, Switzerland) by aligning the oil bubble indicating the injection site. The movements of the fluorescent particles were analyzed using the automatic spots-tracking module of the Imaris software package to generate a set of position lists (Bitplane Inc., Zurich, Switzerland). The diameter of the spots was measured as an average of 0.4 micron and the threshold filter used was the spot quality. Tracks were generated using the autoregressive motion algorithm using 5 micron maximal distance and no gap allowed. Tracks appeared shorter than 3 frames or those exhibiting a high y displacement value and make few steps only in the y axis (track AR1 Y filter) were eliminated. Note that a track as recognized by the software may consist of several segments moving back and forth, until the particle moves out of the plane of focus and disappears. Following automated tracking, artificial tracks were removed manually from the data set.

All statistical analysis, histograms and heat maps were done by Matlab software (R2009a, The MathWorks, Natick, MA). Track traces were done by plotting directed velocities in the y axis and time in seconds on the x axis. Instantaneous velocities were analyzed using velocity data for each bead in each time frame. Positive values represent anterograde movements and negative values represent retrograde movements.

Histogram plots of velocity values for each treatment were generated by binning velocity values into equally spaced intervals of 0.04 µm/sec and normalized to an area of 1.

The track traces were split into the anterograde and the retrograde segments. The end of each segment was defined by either stop or change of direction of the movement. Run length and the duration of each segment were binned into two-dimensional histogram and shown as a heat map: x axis represents segment duration in seconds and y axis represent total run length in microns. The coloring represents the relative abundance of segments with the designated run length and duration, ranging from zero in black to 1% colored white.

85% cutoffs for run length and segment duration were calculated separately after excluding the segments that belong to bins with more than 1% of the data (normally segments with a close to zero run length).

To calculate the median and 95-quantile velocities we excluded the segments with a close to zero run length and calculated the average velocity for each of the remaining segments. p-values for cutoffs were estimated by bootstrapping with re-sampling and with replacement.

GST proteins were isolated from bacteria and purified on glutathione agarose according the manufacturer's instructions (Sigma-Aldrich, St. Louis, MO). Recombinant proteins (5 µg) were incubated with freshly prepared brain lysate, cell lysates or other recombinant proteins for 2 hours using slow rotation at 4°C. Lysate were prepared by homogenization of brain or cells in Tris-Triton buffer (20 mM Tris pH 8, 100 mM NaCl, 1 mM EDTA, 0.1% Triton with protease inhibitor). After two hours washed glutathione agarose beads were added and incubated for another hour at 4°C. The agarose beads were pelleted in a centrifuge at 2400 g for 5 minutes and washed three times. Following the third wash 2× protein sample buffer was added to the pellet and the samples were heated to 95°C for five minutes. Proteins were separated on 10–12% SDS-PAGE acryl amide gels and either transferred to Nitrocellulose membranes for Western blotting or stained directly with Coomassie Brilliant Blue (Sigma-Aldrich, St. Louis, MO) for mass spectrometry identification. Antibodies used for Western blot analysis included anti kinesin heavy chain (Millipore, Billerica, MA) and anti-dynein intermediate chain (Santa Cruz Inc., Santa Cruz, Ca). Loading control of GST proteins is done by ponceau staining of membrane, all other proteins were checked with specific antibodies (anti-Ndel1, Anti-GFP and Anti His) indicated as “input”.

Palmitoylation of Ndel1 was followed using metabolic labeling with 17-Octadecynoic acid (17-ODYA, Cayman Chemical, Ann Arbor, MI) of HEK293 cell line over expressing mCherry-Ndel1, GFP-pep3 with or without HA-DHHC7 (Martin and Cravatt, 2009). Cells were metabolically labeled by addition of 25 µM 17-ODYA to the media eight hours after calcium-phosphate transfection for overnight. Following which, cells were washed and lysed. Ndel1 was immunoprecipitated with rat anti-Ndel1 antibodies (received from Dr Aoki, Tokyo University) as previously described (Shmueli et al., 2010). The immunoprecipiate was subjection to click chemistry which involved the addition 1 mM Tris (2-arboxyethyl) phosphine (TCEP, Sigma-Aldrich, St. Louis, MO) dissolved in water, 100 µM Tris [(1-benzyl-1H-1,2,3-triazol-4-yl) methyl] amine (TBTA, Sigma-Aldrich, St. Louis, MO) dissolved in DMSO/t-butanol (20%/80%), 1 mM CuSO4 and 20 µM Azide modified Alexa Flour® 647 (655 emission/650 excitation Molecular probes, Invitrogen, Carlsbad, CA) in PBS. Reaction was left on the bead-proteins complex for one hour at room temperature and then washed. Proteins were separated from beads to 2× protein sample buffer by heating to 95°C for five minutes followed by SDS-PAGE. Gels were scanned for a fluorescent signal with Typhoon™ 9400 (GE Healthcare Bio-Sciences Corp., Piscataway, NJ) 670 emission/633 excitation. Western blots were conducted using Rat anti Ndel1 and mouse anti GFP (Sigma-Aldrich, St. Louis, MO).

His-tag conventional kinesin proteins previously described (Coy et al., 1999) were isolated from bacteria and purified using HisPur™ Ni-NTA Resin (Thermo scientific). Three constructs were used: full length kinesin heavy chain homodimer (α₂), the recombinant wild-type kinesin consisting of two heavy chains and two light chains (α₂β₂) and delta-tail construct which is kinesin heavy chain deleted after residue M559. Five microgram of each protein were separated on 10% SDS-PAGE and transferred onto a nitrocellulose membrane. The membranes were then incubated with declining concentrations of guanidine-HCl (6M, 3M, 1M and 0.1M) in AC buffer (10% glycerol, 100 mM NaCl, 20 mM Tris pH 7.5, 1 mM EDTA, 0.1% Tween-20) with 2% milk and 1 mM Dithiothreitol for 30 minutes for each concentration. Blocking was done in AC buffer with 2% milk and 1 mM Dithiothreitol for at least one hour. Five microgram per milliliter of each of the recombinant proteins probes (GST, GST-Ndel1, GST-Pep3, GST-DID and GST-Nde1) were incubated with the membranes overnight at 4°C rotating. Afterwards, the membranes were washed 3 times for 20 minutes in AC buffer and incubated with mouse anti-GST antibodies (UC Davis/NINDS/NIMH NeuroMab Facility). Following incubation with HRP-anti-mouse secondary antibodies the blots were developed with ECL and exposed to films. The membranes were then stained with Ponceau S to visualize the bait proteins.

We thank current and previous lab members for their contribution, Dr Elaine Bearer for her generous introduction to the squid system and usage of her instruments, Derek Nobrega for showing us the squid axon related experimental procedure, Joe DeGiorgis, Captain Henry Klimm, Captain Bill Klimm, Ed Enos and Louis Kerr for help in Woods Hole. O.R. is an Incumbent of the Berstein-Mason professorial chair of Neurochemistry.

M.E. is supported by the Gerhardt M. J. Schmidt Minerva Center for Supramolecular Architecture, and by the historic generosity of the Harold Perlman Family. Our research (O.R. and M.S.) has been supported in part by the Gruss Lipper Family Foundation Fellowship, the Israel Science Foundation (grant no. 270/04), the Legacy Heritage Biomedical Program of the Israel Science Foundation (grant no. 1062/08), BSF grant 2007081, Minerva foundation with funding from the Federal German Ministry for Education and Research, ERANET-NEURON (DISCover, IMOS 3-00000-6785), the Benoziyo Center for Neurological diseases, the Helen and Martin Kimmel Stem Cell Research Institute, and the David and Fela Shapell Family Center for Genetic Disorders Research.