Single-molecule tracking of tau reveals fast kiss-and-hop interaction with microtubules in living neurons.

The microtubule-associated phosphoprotein tau regulates microtubule dynamics and is involved in neurodegenerative diseases collectively called tauopathies. It is generally believed that the vast majority of tau molecules decorate axonal microtubules, thereby stabilizing them. However, it is an open question how tau can regulate microtubule dynamics without impeding microtubule-dependent transport and how tau is also available for interactions other than those with microtubules. Here we address this apparent paradox by fast single-molecule tracking of tau in living neurons and Monte Carlo simulations of tau dynamics. We find that tau dwells on a single microtubule for an unexpectedly short time of ∼40 ms before it hops to the next. This dwell time is 100-fold shorter than previously reported by ensemble measurements. Furthermore, we observed by quantitative imaging using fluorescence decay after photoactivation recordings of photoactivatable GFP-tagged tubulin that, despite this rapid dynamics, tau is capable of regulating the tubulin-microtubule balance. This indicates that tau's dwell time on microtubules is sufficiently long to influence the lifetime of a tubulin subunit in a GTP cap. Our data imply a novel kiss-and-hop mechanism by which tau promotes neuronal microtubule assembly. The rapid kiss-and-hop interaction explains why tau, although binding to microtubules, does not interfere with axonal transport.

INTRODUCTION

The cytoskeleton and its dynamic regulation are central for practically all aspects of cellular functions. In neurons, the microtubule-associated protein (MAP) tau regulates microtubule (MT) dynamics and plays a major role in neurodegenerative diseases such as Alzheimer's disease (Ittner and Götz, 2010; Irwin ; Spillantini and Goedert, 2013). The general belief is that the vast majority of tau molecules decorate axonal MTs, thereby stabilizing them. During disease, hyperphosphorylation or aggregation of tau supposedly causes the detachment of tau from MTs, followed by their breakdown (Selkoe, 2013).

In addition to its role in regulating microtubule dynamics, it has also become evident that tau is a multifaceted protein that interacts with various proteins and fulfills diverse functions in different neuronal compartments (Morris ). However, it is unknown how tau exerts its non–MT-related functions in the cell body or dendrites when being attached to MTs (Ittner ). It is also surprising that neither elevated nor reduced tau levels affect axonal transport processes despite tau being a potential obstacle for motor proteins (Yuan ).

To address the apparent paradox that tau regulates MT dynamics but does not interfere with MT-dependent axonal transport, we performed fast single-molecule tracking (SMT) of tau in living neuronally differentiated cells and primary neurons. The spatial resolution of SMT (∼20 nm) was high enough to discriminate between tau being attached to one particular MT or to its neighbor. Step-size distribution (SSD) analysis of pseudotrajectories revealed hops of tau between neighboring MTs, indicating that tau can have multiple binding partners within a few tens of milliseconds. The results were validated by Monte Carlo simulation of tau hopping in a cellular process filled with MTs. Tau-induced changes in MT dynamics were monitored by quantitative imaging using fluorescence decay after photoactivation (FDAP) recordings of photoactivatable green fluorescent protein (PAGFP)–tagged α-tubulin.

This is the first study in which the interaction of tau with MTs in living neurons has been evaluated by direct single-molecule observations. The approach provides a novel method to follow physiological or pathological changes in the tau–MT interaction in living neurons with unprecedented precision.

RESULTS

Single tau molecules show undirected fast movement in processes of neural cells

To assess the behavior of tau in neuronal processes, single, fluorescently labeled tau molecules were tracked with a temporal resolution of 5 ms. HaloTag was fused to the amino terminus of human tau, where it does not interfere with the MT interaction that is mediated by tau's carboxy-terminal half (Gauthier-Kemper ; Figure 1A, top). The focal plane of the total internal reflection fluorescence (TIRF) microscope was adjusted to capture MTs (Figure 1A).

FIGURE 1:

Single-molecule tracking reveals undirected fast movement of tau in processes of neural cells. (A) Schematic representation of the Halo-tagged tau fusion construct (Halo-tau) and the visualization of a thin cellular layer by TIRF microscopy. The microtubule-binding repeats (R1–R4) are indicated by the gray boxes, and the evanescent field close to the cellular attachment site is shown in yellow. Bottom, colocalization of Halo-tau (left) and monomeric enhanced GFP–tagged α-tubulin (EGFP-tub; right) as visualized by TIRF microscopy of the cell body of living PC12 cells. Sum projections of 1000 consecutive frames were recorded. Scale bar, 10 μm. MTOC, microtubule-organizing center. (B) Time series of individual Halo-tau molecules moving in a process of a neuronally differentiated PC12 cell. Motion during the first 85 ms. Dashed lines indicate the border of the process. A trajectory generated from the complete time series (2.2 s = 440 frames) is shown below, indicating undirected fast movement both longitudinally and transversally. The starting point is indicated by a circle and the end by a square. Scale bar, 0.5 μm.

A diffusing tau molecule is so fast that it produces spatially blurred fluorescence that submerges in the background. In fact, a freely diffusing tau molecule would cover a distance of ∼400 nm during 5 ms. Thus the signal intensity clearly surmounts the background only if tau is bound to a MT filament, which allows tracking subsequent MT-binding events of single tau molecules. Individual tau molecules bound to MTs could be localized with a precision of ∼20 nm by determining the centroid of its diffraction-limited fluorescence signal. Unexpectedly, we observed tau molecules rapidly moving within some tens of milliseconds between binding sites that were several hundred nanometers apart (Figure 1B, top, and Supplemental Movie S1). A typical trajectory of 440 frames revealed undirected, fast, longitudinal and transversal displacements of fluorescent tau in the cellular process (Figure 1B, bottom). This appears inconsistent with a previous study (Konzack ). By using fluorescence recovery after photobleaching (FRAP), the authors recorded the recovery of fluorescence in a large ensemble of tau molecules and, by fitting their data by some reaction-diffusion model, they estimated a rather long dwell time of tau on MTs, namely ∼4 s. It was also speculated that bound tau might slowly slide on the surface of MTs (Hinrichs ). Previous treatments of FDAP/FRAP data referred to the whole MT bundle. Here, for the first time, the binding to single MT filaments is resolved in living cells. Deficiencies of the previous models and a rigorous reaction-diffusion model accounting for the bundle-structure of MT filaments, which leads to a 100-fold correction to the long apparent dwell-time, will be detailed elsewhere (Igaev ).

Tau interacts with microtubules via a kiss-and-hop mechanism in processes of neural cells

The spatial resolution of single-molecule tracking (∼20 nm) was high enough to discriminate between tau being attached to one particular MT or to its neighbor in processes of PC12 cells (Figure 2A). Note that the spacing between MTs in PC12 neurites is ∼70 nm (Jacobs and Stevens, 1986), a much higher value than the typical MT–MT distance in axons (Hirokawa and Takemura, 2005) due to the fact that PC12 cells do not develop axonal compartmentalization. Figure 2B (see also Supplemental Movie S2) documents the “pseudotrajectory” of binding events of tau over time (for a real trajectory, see Figure 2C). It is also evident that tau binds to several neighboring MTs in the transversal direction during the recording time (Figure 2B).

FIGURE 2:

Tau interacts with multiple microtubules in processes of PC12 cells. (A) Schematic representation showing the geometry and MT distribution in a process of a neuronally differentiated PC12 cell. (B) Localization of a single Halo-tau molecule in a PC12 cell process over time (for a high-resolution version see Supplemental Figure S1). For comparison, thickness and density of MTs are schematically indicated by green bars. Top, enlargement of the indicated segment. (C) Schematic representation of a real trajectory of a molecule moving between microtubules in a bundle is shown in black. Yellow and black circles indicate binding events of tau with microtubules. The recorded pseudotrajectory is indicated in red. Only the binding events within the time frame (Δt) are detected and indicated by the yellow circles. (D) Time series showing the movement of a syp-mCherry–tagged vesicle in a PC12 cell process. In contrast to tau, localization of syp-mCherry over time indicates directional movement along a single MT. Scale bar, 1 μm.

As a control, we analyzed the behavior of synaptophysin (syp), a marker for fast axonal transport (Supplemental Movie S3). If tau hops between microtubules on a stationary microtubule array, we would expect that syp would not reveal transversal displacement between adjacent MTs. Indeed, a representative two-dimensional (2D) projection of binding positions of syp over time (Figure 2D, top, and Supplemental Movie S4) showed directional lateral movement on a single MT filament with small transversal spread. This behavior clearly differed from what was observed for the tau protein.

We observed fast longitudinal and transversal displacements of tau in the cellular process (Figure 1B). To characterize quantitatively the displacements of tau, we calculated SSDs. Here the probability of a particle moving over a certain distance within a fixed time interval ∆t was evaluated. Figure 3A represents an SSD histogram plotted using >1000 2D pseudotrajectories for the first time step (see Materials and Methods). The SSD analysis of the pseudotrajectories revealed the presence of three peaks, at ∼95, ∼190, and ∼285 nm. Remarkably, these values correspond to the distance between the centers of neighboring MTs, reflecting hops 1 → 2 and multiples thereof (hops 1 → 3, 1 → 4; Figure 3A), indicating that tau can have multiple binding partners within a few tens of milliseconds. Note that peaks at intermediate locations, which might be expected due to different z-positions of neighboring MTs, are merged due to the 2D projection and contribute to the widening of the peaks.

FIGURE 3:

Step-size distribution analysis and Monte Carlo simulation indicate jumps of tau between microtubules. (A) SSD analysis of >1000 pseudotrajectories indicating subpopulations of jumps between MTs. The peaks in the histogram correspond to 1 → 2, 1 → 3, and 1 → 4 hops between MTs. (B) SSD analysis of simulated pseudotrajectories at different pseudoequilibrium constants (kon*/koff). Note that the histogram at kon*/koff = 102 closely resembles the data shown in A. The Monte Carlo simulation revealed an additional peak (indicated in red), which was not resolved by imaging.

To validate this result, we simulated by Monte Carlo algorithm the hopping of tau in a cellular process filled with microtubules. Two input parameters were fixed, namely, the diffusion constant of free tau (14.4 μm2/s; see Materials and Methods) and the average spacing of MTs of 70 nm. The rate constants for association and dissociation of tau were varied (Figure 3B). The simulated SSD resembles the data shown in Figure 3A for kon*/koff = 102. This value implies 99% of tau being bound, which is compatible with data on tau–MT binding from cell-free experiments at cellular tubulin concentration (Brandt and Lee, 1993b). In vitro, tau binds to microtubules with relatively weak affinity (Kd = 1.1 × 10−7 M; Butner and Kirschner, 1991) and the saturating molar ratio measured for tau protein binding to polymerized tubulin has been determined as 1:5 (Cleveland ; Hirokawa ). On the basis of these numbers, one calculates kon*/koff as ∼50, which is also close to the ratio where the simula­tion revealed peaks. Note that the simulation revealed an additional peak at ∼20 nm. This peak represents hops between MTs whose projections on the 2D plane closely neighbor but can be more distant in three dimensions. This peak remains unresolved in the experiment because it was close to the localization precision of ∼20 nm.

Single-molecule experiments allow estimating the dwell time of tau on a single microtubule. The respective frequency histograms of residence times of tau on MTs could be best fitted to a single exponential, indicating that tau has only one major interaction partner in neuritic shafts, namely, MTs. Time constants of the exponential fits were referred to as dwell times (Figure 4). We found that tau dwells on a single MT for ∼40 ms, which is shorter by two orders of magnitude than previously reported (Konzack ).

FIGURE 4:

Representative histogram of residence times of tau on microtubules in PC12 cells. The monoexponential fit is shown by a red line, and the time constant (dwell time) is indicated on top.

Tau promotes microtubule assembly despite its rapid kiss-and-hop dynamics

The data show that tau only briefly resides on one particular MT until it hops to a neighboring one. We asked whether this short-lived interaction is compatible with an effect of tau on microtubule polymerization in living cells. To address this question, we coexpressed mCherry-tagged tau or, as a control, 3×mCherry together with PAGFP-tagged α-tubulin in neuronally differentiated PC12 cells. Tau-induced changes in MT dynamics were then monitored by FDAP measurements from cells in which the fluorescence of PAGFP had been focally activated in the middle of a process (Gauthier-Kemper ; Figure 5, A and B). The normalized intensity of the control (3×PAGFP) showed a considerably faster decay than that of α-tubulin (Figure 5C), indicating reaction-dependent diffusion of tubulin versus free diffusion of 3×PAGFP. The presence of mCherry-tau significantly slowed the decay of α-tubulin, whereas 3×mCherry did not have any effect (Figure 5C, left vs. right). To determine the fractions of free and bound tubulin, we fitted the FDAP curves to a reaction-diffusion model (see Materials and Methods). In the absence of exogenous tau, about two-thirds of PAGFP-tubulin was polymerized in the cell processes. The exogenous mCherry-tau increased polymerization by ∼10% (Figure 5C, bottom). This was consistent with previous results from in vitro experiments, in which tau was shown to induce microtubule polymerization in a concentration-dependent manner (Brandt and Lee, 1993a). Taken together, the data indicate that tau is capable of promoting MT assembly in neurites despite its rapid kiss-and-hop dynamics. However, it should be taken into account that the experiment uses an overexpression system in which tubulin and tau are expressed on top of the respective endogenous proteins. It remains to be shown whether it is also the physiological function of tau to promote microtubule assembly in neurites.

FIGURE 5:

Tau promotes microtubule assembly in neurites of living cells. (A) Live-cell imaging of a PC12 cell coexpressing mCherry-tagged tau (left, inset) and PAGFP-tagged α-tubulin (PAGFP-tub) before and after fluorescence photoactivation. The position of activation is indicated by a violet box. Fluorescence decay of PAGFP-tub in the activated region is indicated by the arrows. Scale bar, 10 μm. (B) Determination of the axes of symmetry (left) and the plot of normalized intensity distributions (right) as they were used to define the ROI. (C) FDAP curves of PAGFP-tub in the ROI. The presence of mCherry-tau decreases fluorescence decay of PAGFP-tub (left), whereas a control protein (3×mCherry) has no effect (right). For comparison, FDAP of a soluble protein (3×PAGFP) is shown. Values are expressed as mean ± SEM (n = 29–36). Fractions of polymerized PAGFP-tub were estimated based on fitting of the respective FDAP curves to a reaction-diffusion model (bottom).

Tau exhibits kiss-and-hop behavior also in axons of primary neurons

PC12 cells are a well-characterized neuronal model with a rather homogeneous cell population and well-defined microtubule distribution in processes. They do not, however, develop axonal-somatodendritic polarity. To examine the behavior of tau in axons, we prepared primary mouse cortical cultures and introduced Halo-tagged tau by lentiviral gene transfer (Bakota ; Figure 6A, left). To avoid potential interference of endogenous mouse tau, we prepared cultures from TAU− animals. A Western blot confirmed 1) the absence of endogenous mouse tau in the TAU− cultures (arrowheads, Figure 6A, right) and 2) the expression of Halo-tau at a level comparable to that of the endogenous one (arrow, Figure 6A, right).

FIGURE 6:

Tau exhibits kiss-and-hop behavior also in axons of primary neurons. (A) Triple fluorescence micrographs showing expression of Halo-tau (red) in lentivirally infected wild-type mouse cortical neurons. Staining against tubulin (DM1A, green) and nuclei (4′,6-diamidino-2-phenylindole [DAPI], blue) is shown. Right, Western blot indicating expression of exogenous Halo-tau (arrow) in cultures from wild-type (wt) or TAU− mice (KO). The presence of endogenous tau in wt mice is indicated by arrowheads. Bottom, staining for actin as a loading control. Scale bar, 20 μm. (B) Triple fluorescence micrographs after a fixation-extraction protocol, which reveals cytoskeletal association. Cortical neurons from TAU− mice were lentivirally infected to express Halo-tau (red). Top, overview; bottom, individual growth cone region. Note the enrichment of Halo-tau at the transition between the axon shaft and the growth cone as indicated by arrows. Staining against tubulin (DM1A, green) and nuclei (DAPI, blue) is shown. Scale bar, 20 μm (top), 10 μm (bottom). (C) Localization of a single Halo-tau molecule in an axon of a TAU− neuron over time. For comparison, thickness and density of axonal MTs are schematically indicated by the green bars. Top, enlargement of the indicated segment. Note that the MT–MT distance in axons of cortical neurons is much smaller than in neurites of PC12 cells. (D) Dwell times of Halo-tau in neurites of PC12 cells (39 ± 4 ms; n = 7) and axons of primary cultures (36 ± 5 ms; n = 4). Data represent mean ± SEM.

Schematic representation illustrating features and consequences of the kiss-and-hop behavior of tau as identified in this study. P, phosphorylation; D, dephosphorylation.

Halo-tau was ubiquitously present in the neurons, similar to what has been observed after expression of exogenous tau in cultured neurons or transgenic animals (Konzack ; Hundelt ). This might be due to the fact that the expression constructs lack the 3′ untranslated region of the tau mRNA, which targets the mRNA to the proximal axon (Litman ; Aronov ). To prove that Halo-tau binds to axonal microtubules, we used a previously described extraction-fixation protocol to visualize tau's microtubule interaction. We observed that Halo-tau clearly associated with axonal MTs (Figure 6B). Furthermore, binding was highest at the transition between the axon shaft and the growth cones (Figure 6B, arrows), similar to what was previously observed for endogenous tau in cultured neurons (Kempf ). Thus the experimental approach appears to be appropriate to analyze the mode of the interaction of tau with axonal microtubules. Note that primary neurons have a much higher density of MTs in their axons than is present in the processes of PC12 cells. The typical axonal MT–MT distance was reported to be ∼20 nm (Hirokawa and Takemura, 2005). However, despite this difference, a representative 2D projection of binding positions of a single tau molecule revealed very similar kiss-and-hop behavior, indicating that it is independent of MT density and the presence of endogenous tau (Figure 6C and Supplemental Movie S5). Moreover, we found a dwell time of ∼40 ms on the MT surface (Figure 6D; see Materials and Methods). Therefore there was no significant difference between the dwell times of tau molecules in processes of PC12 cells and axons of primary cultures.

DISCUSSION

The foregoing data represent the first recording on a single-molecule level of the interaction of a structural MAP with the MT network in living neurons. The major findings are the following: 1) the dwell time of bound tau on each MT filament (∼40 ms) is shorter by two orders of magnitude than previously believed (Konzack ), 2) tau hops between MTs in an undirected manner, and 3) despite its high mobility, tau is capable of shifting the equilibrium between free tubulin and the MT network toward polymerization.

Our findings have several consequences for the physiological role of tau (schematized in Figure 7). In vivo the velocity of axonal transport is virtually unaffected by tau (Yuan ), which is fully compatible with its kiss-and-hop kinetics. Indeed, knowing the concentration of tau molecules in neurons and the ratio of bound molecules, one calculates the mean lateral distance between the two closest tau mole­cules bound to the same MT simultaneously to be ∼45 nm (see Materials and Methods). Kinesins move on the MT surface with an average velocity of ∼1 μm/s (Kaether ). During the mean dwell time of tau (∼40 ms), they progress by ∼40 nm. Therefore kinesins can readily proceed on the MT without being impeded by tau standing in their way. If the dwell time was much longer, tau, which shares binding sites with kinesins on MTs (Seeger and Rice, 2010), would impair the axonal transport. It is noteworthy that the rapid dynamic behavior of tau is also compatible with its non–MT-related functions and its multiple interactions in the cell (Morris ).

The view that tau binding to MTs is regulated by phosphorylation and dephosphorylation cycles (Ballatore ) obviously conflicts with tau's rapid dynamic behavior. The catalytic rate of a typical protein kinase is ∼10 reactions/s (Newton, 1995). Given the average number of MT interactions that a single tau molecule experiences per second (∼2000 reactions/s), phosphorylation cannot directly modulate the rapid tau–MT binding kinetics within one phosphorylation cycle.

We provided evidence that tau is capable of promoting MT assembly in neurites despite its rapid kiss-and-hop dynamics. According to a study by Gardner , the MT self-assembly kinetics is strongly dependent on the concentration of free tubulin subunits and is an order of magnitude more rapid than previously believed (the lifetime of a tubulin subunit in a GTP cap is ∼1 ms at the physiological concentration of tubulin subunits). It appears unlikely that such a dynamic system as the MT network is regulated by a MAP whose dwell time on MTs is in the range of several seconds. This rapid MT kinetics also implies that even a weak suppression of the tubulin off-rate strongly affects catastrophe and rescue events in the axon. Our data suggest a new mechanism of axonal MT regulation by tau, namely, the short interaction of tau with MTs (∼40 ms) suppresses the natural rapid dissociation of tubulin dimers at the MT tip, leading to an increase of the MT GTP cap size. It is noteworthy that such a function of tau is fully compatible with its distribution in the axon, since we observed that the association of tau with the cytoskeleton is highest at the transition between the axon shaft and the growth cone (Figure 6B), where microtubules are known to be most dynamic.

Thus the revealed dwell time of tau (∼40 ms) appears optimal for the following two reasons. On the one hand, it is short enough not to inhibit axonal transport, and on the other hand, it is long enough to regulate MT dynamics. Pathological shifts in tau's MT association or dwell times (e.g., induced by long-term abnormal phosphorylation of tau) may provoke changes in the dynamic equilibrium of tau with MTs and hence result in a disturbed tubulin–MT equilibrium.

MATERIALS AND METHODS

Materials and antibodies

Chemicals were obtained from Sigma-Aldrich (Deisenhofen, Germany), cell culture media and supplements were obtained from Sigma-Aldrich and Invitrogen, and culture flasks, plates, and dishes were obtained from Thermo Fisher Scientific (Schwerte, Germany), unless stated otherwise. The following antibodies were used: anti-tubulin (mouse, DM1A), Tau5 (mouse; BD Biosciences, Heidelberg, Germany) and anti-actin (mouse; Calbiochem/Merck, Darmstadt, Germany). As secondary antibodies, peroxidase- and Alexa Fluor 488–conjugated anti-mouse antibodies (Jackson ImmunoResearch Laboratories, West Grove, PA) were used.

Construction of expression vectors, transfections, and lentiviral infections

For bright and highly photostable signals from single tau molecules during live-cell imaging, we used the enzymatic HaloTag system and the HaloTag TMR ligand (Los ). The eukaryotic expression plasmid for adult (441 amino acids) human tau with amino-terminally fused HaloTag tag were constructed from pRc/cytomegalovirus (CMV)-htau441wt expression vectors (Gauthier-Kemper ) by addition of restriction sites NotI and XhoI flanking to tau. The resulting PCR product was sequentially cut with the respective enzymes (Fermentas, St. Leon-Rot, Germany) and cloned into pSems-Halo-1-26m (Wilmes ). For generating pSems-mEGFP-α-tubulin the respective α-tubulin sequence was amplified by PCR from pIREShyg-PAGFP-α-tubulin (Tulu ; plasmid 12296; Addgene) and inserted into pSems-mEGFP restricted by BamHI and NotI (Wilmes ). Synaptophysin-mCherry was synthesized with the restriction sites BamHI and XbaI from Life Technologies (Darmstadt, Germany) and cloned into pRc/CMV (Gauthier-Kemper ). The pCMV-3×PAGFP plasmid was described previously (Weissmann ).

To obtain the L22-3×mCherry construct, named pJJH1392, in vivo recombination was used employing triple mammalian/yeast/bacterial shuttle vectors as described previously (Bakota ). Specifically, the coding sequence for 3×mCherry was PCR amplified from pJJH1255, a derivative of pCM79, kindly provided by M. Knop (University of Heidelberg, Germany; Maeder ), using the oligonucleotide pair 11.357/11.358 (AATCTTTCACAAATTTTGTAATCCAGAGGTTGATTATCGATAAGAAGTGGCGCGCTCACTCGAG, TAGTTCTGGGGG­CAGCGGGGATCCCCGGGTACCGGTCGCCACCGAAGCTTCGTACGCTGCAGG). The PCR product was cotransformed into the diploid yeast strain DHD5 with pJJH1308 linearized with BamHI and selected for growth on medium lacking uracil (Bakota ). Isolation of plasmid DNA, amplification in Escherichia coli, and sequencing confirmed that the recombination occurred as expected. For construction of the L22-mCherry-htau441wt viral vector, named pJJH1311, a similar strategy was followed based on in vivo recombination with pJJH1308 linearized with BamHI. As a template to obtain the mCherry-htau441wt coding sequence, we used pJJH1309 with the oligonucleotide pair 11.045/11.046 (GGGCAGGCGAGTGGCCCCTAGTTCTGGGGG­CAGCGGGGATATGGTGAGCAAGGGCGAGGAGG, AGGTTGATTATCGATAAGCTTGATATCGAATTCTAGAaccggtgcgcggatc­CGACCTATCACAAACCCTGC, with the lowercase letters indicating sequences inserted for introduction of restriction sites, which appear neither in the template nor in the recipient vector). pJJH1309 was also obtained by in vivo recombination with linearized pJJH900 (containing a Flag-htau441wt sequence cloned under the control of the yeast GAL1/10 promoter as a BamHI/SphI fragment into pJJH447; Lorberg ). The mCherry sequence was introduced as a PCR product using pAJ001 (Jendretzki ) as a template with the oligonucleotides 11.043/11.044 (ATCCAAAAAAAAAGTAAGAATTTTTGAAAATTCAATATAAATGGTGAGCAAGGGCGAGGAGG; CGAACTCCTGGCGGGGCTCAG­CCATCGATTTATCGTCATCCTTGTACAGCTCGTCCATGCCGCC) and substituted the Flag tag in pJJH900. pJJH1311 was then linearized with Sbf1, and the mCherry sequence was substituted for a Halo-tag by in vivo recombination in yeast with a PCR product obtained from pSems-Halo-tau using the oligonucleotides 11.037/11.519 (cgcgCTCGAGTTTCAGGCCAGCGTCCGTGTCAC, binding with the uppercase letter sequences within the tau sequence, and AGCACGGGCAGGCGAGTGGCCCCTAGTTCTGGGGGCAGCGGGGAATGGGATCCGAAATCGGTACTGG, which provides the homology to the CamKII promoter in the vector sequence). The resulting mCherry-htau441 construct was named pJJH1440 and used to generate viral vectors.

Transfections of PC12 cells were performed with Lipofectamine 2000 (Invitrogen) essentially as described previously (Fath ). For production of lentivirus, human embryonic kidney cells 293FT (Invitrogen) were transfected with the expression vector, and two helper plasmids and viral particles from the supernatant were concentrated by ultracentrifugation as described previously (Bakota ).

Cell culture and live-cell labeling

PC12 cells were cultured in serum-DME as described previously (Fath ). Undifferentiated cells were plated on 35-mm poly-l-lysine– and collagen-coated glass-bottom culture dishes (MatTek, Ashland, MA) at 103 cells/cm2 and cultured in DME with 1% (vol/vol) serum. For TIRF imaging, cells were plated on glass coverslips (24 mm, No. 1; VWR, Radnor, PA) that were coated with a poly-l-lysine-graft-(polyethylene glycol) copolymer functionalized with RGD as described previously (VandeVondele ). Cells were flattened by addition of 100 ng/ml 7S mouse nerve growth factor (NGF; Alomone Laboratories, Munich, Germany) for 1–2 d. For neuronal differentiation, NGF treatment was extended to 4–6 d, and the medium was changed every 2–3 d. Before imaging, the medium was exchanged against the same medium containing DME without phenol red. Primary cortical cultures were prepared from cerebral cortices of mouse embryos (days 14–16 of gestation) and cultured as described previously (Bakota ). The cultures were obtained from C57BL/6 and TAU− mice (Dawson ). Cells were plated at 2.5 × 104 cells/cm2 on polylysine- and laminin-coated coverslips.

Labeling of PC12 cells and primary cultures was performed by incubation with serum-DME containing 1–5 nM TMR-HTL for 15 min at 37°C. This corresponded to a labeling efficiency of <0.02% of the exogenous Halo-tau. Subsequently, cells were washed with serum-DME without phenol red and then left to recover for 1/2 h in the incubator, followed by an additional medium exchange. All animals were maintained and killed according to National Institutes of Health guidelines and German animal care regulations.

Immunocytochemistry

Immunocytochemistry after fixation with 4% paraformaldehyde and a combined NP-40 detergent extraction-fixation procedure to reveal cytoskeletal association were performed as described previously (Gauthier-Kemper ). Fluorescence microscopy used a dry 40× (numerical aperture [NA] 0.75) objective lens on a fluorescence microscope (Eclipse TE2000-U; Nikon, Tokyo, Japan) equipped with a digital camera (COOL-1300; Vosskühler, Osnabrück, Germany). Image analysis used the program Lucia G (Nikon).

Photoactivation and live imaging

Live imaging was performed on a laser scanning microscope (Eclipse TE2000-U inverted; Nikon) equipped with argon (488 nm), helium/neon (543 nm), and blue diode (405 nm) lasers. PAGFP-tubulin expressing cells were visualized with a Fluor 60× (NA 1.4) ultraviolet-corrected (VC) objective lens. The microscope was enclosed in an incubation chamber maintained at 37°C and 5% CO2 (Solent Scientific, Fareham, UK). Automated image acquisition after photoactivation was essentially performed as described previously (Weissmann ). Photoactivation was performed with the blue diode, producing an activation spot 5 μm in diameter. Frames were obtained at a frequency of 1 frame/s, and 112 frames were collected per experiment. Standard series were collected at a resolution of 256 × 256 pixels.

Single-molecule microscopy, localization, and tracking

For recording single molecules, TIRF microscopy was performed using an Olympus excellence cell^TIRFM microscope equipped with 405-nm (200 mW), 488-nm (200 mW), 561-nm (200 mW), and 647-nm (140 mW) lasers (Olympus, Tokyo, Japan), as well as a backilluminated electron multiplied charge-coupled device camera (C9100-13; Hamamatsu, Hamamatsu City, Japan). Fast single-molecule tracking used a digital complementary metal–oxide–semiconductor camera (ORCA-Flash4.0 V2 C11440-22CU; Hamamatsu). A 150× magnification objective with NA 1.45 (UAPON 150×/1.45, Olympus) was used for TIR illumination. The excitation beam was reflected into the objective by a quad-line dichroic beam splitter for reflection at 405, 488, 568, and 647 nm (Di01 R405/488/561/647; Semrock, Rochester, NY). The emitted light from the sample was filtered using a quad-band bandpass filter (FF01 446/523/600/677; Semrock). The microscope was enclosed in an incubation chamber maintained at 37°C and 5% CO2 (Olympus).

Localization of single molecules and single-molecule tracking were carried out as previously described (Hess ; Sergé ; Wilmes ). The theoretical localization precision of Halo-tau and syp-mCherry was 16.7 ± 0.2 and 14.2 ± 0.2 nm, respectively. For further analysis, a mask was placed over the processes, and localizations/trajectories outside of this mask were removed. In addition, all trajectories were aligned with the horizontal axis.

SSD analysis

SSDs were calculated as described before (Weimann ). To exclude potential false connections produced by the tracking algorithm, we restricted the duration of the trajectories to the first three steps. The obtained shortened trajectories were used to plot the SSDs for the first time step, t = ∆t.

Dwell-time determination

The dwell time of tau using single-molecule localizations was determined by counting how many consecutive frames a single localized signal stayed within a radius of 50 nm. The respective frequency histograms were fitted to a single exponential to obtain the dwell time.

FDAP data analysis

To plot normalized decay curves, we first extracted the first frame from the image sequence and chose the region of interest (ROI) containing a part of the cellular process with the desired intensity distribution I(x, y). The image was thresholded to suppress the background noise, and the axes of symmetry of I(x, y) were found. The longitudinal (along the process) and transversal (across the process) axes were then used to plot the intensity profile distributions. The distributions were fitted to slightly modified Gaussian profiles:

where f and f are offset parameters and A > 0 and A > 0 are the maximal amplitudes of the intensity profiles. The parameters σ > 0 and σ> 0 represent the spatial half-widths of the longitudinal and transversal distributions, respectively. Then we used σ and σ to define the recording area—a rectangle with the sides 2σ in length and 2σ in width. The average intensity was then acquired from the recording area. More than 20 independent FDAPs for PAGFP-tub coexpressed with either 3×mCherry or mCherry-tau were performed to form the final single FDAP curves. In all cases, the preactivation frame was subtracted from the rest of the frames in order to exclude the preactivation fluorescence from the desired signal. The obtained raw curves were then averaged and normalized on the maximal intensity.

Monte Carlo simulations of tau dynamics

We simulated tau's reaction-diffusion dynamics in a process of a PC12 cell. The simulation space was defined as a tube of length L = 100 μm and radius R = 500 nm. The tube was filled with 60 parallel MTs of RMT = 12.5 nm, with an intermicrotubule distance of 70 nm. The geometry chosen was in good correspondence with the actual one revealed by Jacobs and Stevens (1986) for cellular processes of PC12 cells after 4 d of NGF treatment. At the initial moment, 1000 particles were randomly positioned within the simulation space. Each particle had a probability kon*/(kon* + koff) of being bound to a MT and a probability koff/(kon* + koff) of being placed anywhere in the cytosol.

Once the geometry and the populations had been set, the time loop was started. The time step was equal to ∆t = 1 μs. At each time step, each bound particle had a probability poff of unbinding. In case it was decided for unbinding, the detachment of the particle occurred in the direction locally perpendicular to the MT surface to which this particle was attached. On the contrary, each free particle was first checked for the act of hitting a MT surface. When hit, the particle had a probability pon of being attached. In case of not attaching, the particle again came off in the direction locally perpendicular to the surface. If the free particle did not hit a MT surface or failed attaching to it and came off, a random displacement (∆x, ∆y, ∆z) for this particle was generated. The displacements ∆x, ∆y, and ∆z were calculated as normally distributed random numbers with the SD equal to , where D is the diffusion constant. Reflection boundary conditions at the outer boundary of the simulation space were imposed. The procedure was iterated for 106–107 time steps.

Diffusion constant estimation for PAGFP-htau and PAGFP-tub

The Stokes–Einstein equation, D = kBT/6πηRH, for diffusion of spherical particles in liquid with low Reynolds numbers was used to estimate diffusion constants. In this formula, kB is the Boltzmann constant, T is absolute temperature of the surrounding medium, η is the viscosity of the medium, and RH is the hydrodynamic radius of the molecule of interest. We assumed the medium to have a viscosity of 3.4 cP at T = 310 K, as previously measured for the cytosol of neuronal processes (Popov and Poo, 1992). The hydrodynamic radius, in turn, depends on the structural properties of the molecule of interest. First, we estimated the diffusion constant for GFP and its derivatives (enhanced GFP, PAGFP). According to Hink (2000), the hydrodynamic radius of GFP is ∼2.3 nm. This is in good correspondence with an empirical formula for globular proteins, RH = 0.595M0.427 (mass in kilodaltons), which yields 2.4 nm for 27-kDa mass of GFP. This radius corresponds to D = 28 μm2/s. Several research groups have reported various experimental estimates for the diffusion constant of GFP in different cellular compartments (Swaminathan ; Hink, 2000; Sprague ; Michelman-Ribeiro ). These values varied from 15 to 27 μm2/s. Our FDAP experiments yielded a consensus value of ∼20 μm2/s for PAGFP based on the measurements for 3×PAGFP (D = 13.9 ± 1.0 μm2/s) in processes of neuronal PC12 cells and subsequent recalculation for a smaller mass, since D ∼ M−1/3. We then estimated the diffusion coefficient for htau alone, which has a molecular weight of 45.85 kDa. Previous small-angle x-ray scattering experiments revealed a radius of gyration of 6.5 nm for 441-residue htau (Mylonas ). The ratio of the radius of gyration to the hydrodynamic radius lies in the range 1.2–1.5 for intrinsically disordered proteins (Mukrasch ), thus yielding RH = 5.4 nm for htau. This value is in accordance with the value of 5.76 nm given by an empirical formula for random coils (Mylonas ), Rg = 0.1927n0.588, where n is the number of residues, although htau is not completely unstructured. The experimental estimate for the hydrodynamic radius of htau hence gives a diffusion constant of ∼12 μm2/s. Finally, we estimated the diffusion constant for PAGFP-htau. To our knowledge, the hydrodynamic radius for this fusion has not been measured experimentally. However, we estimated maximal and minimal values RH can take. The maximal value follows from the assumption that the fusion PAGFP-htau is completely unstructured. Using the empirical formula for random coils yields RH,max =7.4 nm. The minimal value RH,min = 3.7 nm follows from the empirical formula for globular proteins. Thus the actual diffusion constant for PAGFP-htau should be between 9 and 18 μm2/s. On the basis of D = 13.9 ± 1.0 μm2/s for 3×PAGFP, we recalculated the diffusion constant for PAGFP-htau (), which was equal to ∼14.4 μm2/s.

We used the same strategy to estimate the diffusion coefficient for PAGFP-tub, which after forming a dimer with endogenous β-tubulin has a mass of 137 kDa. However, unlike htau, tubulin is a globular protein, which encourages use of the empirical formula for globular proteins to estimate D for PAGFP-tub. Doing so yields ∼14 μm/s. Previous measurements using injection of rhodamine-labeled tubulin into squid giant axons gave D = 8.6 μm2/s (Galbraith ). Using the diffusion constant of 3×PAGFP (D = 13.9 ± 1.0 μm2/s), we get 11.8 μm2/s for PAGFP-tub ().

Reaction-diffusion model for FDAP

We consider a photoactivatable molecule that diffuses in an infinite cylinder with a diffusion constant D and obeys first-order kinetics,

where C represents the free molecule species at concentration c, B represents binding sites at concentration b, and S represents the bound molecular species at concentration s. We make the following assumptions:

Before the photoactivation, the molecular ensemble has reached equilibrium with its binding sites, and the equilibrium concentrations ceq, beq, and seq have been settled.

The photoactivation does not disturb the concentration of binding sites (b = beq).

Diffusion is one dimensional and proceeds along the longitudinal axis Oz. The lateral diffusion is negligible.

The molecule of interest diffuses only when free in cytosol.

The PAGFP tag does not affect binding.

The reaction-diffusion equations then read as follows:

where kon* = konbeq is the pseudo-first-order association rate. The relative equilibrium concentrations are given by

where ctot is the total concentration of the molecules of interest. Here we derive the solutions to the system of Eqs. 3 for two simplified cases. The first case relates to pure diffusion, that is, to a scenario when the majority of molecules are free (ceq > seq). In this case, the system of Eqs. 3 reduces to a simple diffusion equation,

This scenario is supposed to describe a freely diffusing protein like 3×PAGFP. The second case concerns a scenario in which diffusion is very fast compared with the time scale of reaction (τ 1/kon*, where τ is the characteristic time scale of diffusion over the activation region). Thus diffusion and binding are largely separable processes, and the free and bound molecules behave as uncoupled species. The system of Eqs. 3 then turns into

This scenario is expected to describe the behavior of PAGFP-tub. In the FDAP experiment, such behavior was typical (Figure 5A). We then calculate the theoretical expression for the normalized FDAPs defined as (where 2σ is the width of the photoactivation region) for the two simplified scenarios:

where τ. The experimental FDAP curves for 3×PAGFP and PAGFP-tub were fitted by the corresponding theoretical scenario of Eqs. 7 by the method of least squares. Fitting with the pure-diffusion scenario yields the diffusion constant, whereas fitting with the second scenario can additionally yield information about the relative free and bound fractions. The fitting was performed using a self-written Python script (van Rossum and de Boer, 1991; Oliphant, 2007). To calculate 95% confidence intervals for the estimated parameters, we used the formulas for nonlinear data fitting given in Bates and Watts (1988).

Calculation of the minimal distance between two bound tau molecules

We first calculated the average area occupied by a bound tau molecule on a single MT:

where SMT is the total area of the MT surfaces in the cellular process and Ntau is the total number of bound tau molecules in this process. Given the total concentration of tau, ctot, the ratio at which it is bound on average over time, , and the geometry of the cellular process as shown in Figure 2A, one can expand Eq. 8 as follows:

where N is the total number of MTs in the process, RMT is the radius of a single MT, R is the radius of the cellular process, and h is the length of an arbitrary process segment. Because 〈A〉 scales as l2, where l is the linear size of the area , one can estimate the minimal distance between two bound tau molecules by extracting the square root from 〈A〉. Substituting the numerical values for the parameters in Eq. 9, one obtains l ≈ 45 nm.

Other methods

Western blot analysis was performed as described elsewhere (Gauthier-Kemper ).