Specific proteins are concentrated within primary cilia, whereas others remain excluded. To understand the mechanistic basis of entry into cilia, we developed an in vitro assay using cells in which the plasma membrane was permeabilized, but the ciliary membrane was left intact. Using a diffusion-to-capture system and quantitative analysis, we find that proteins >9 nm in diameter (∼100 kD) are restricted from entering cilia, and we confirm these findings in vivo. Interference with the nuclear pore complex (NPC) or the actin cytoskeleton in permeabilized cells demonstrated that the ciliary diffusion barrier is mechanistically distinct from those of the NPC or the axon initial segment. Moreover, applying a mass transport model to this system revealed diffusion coefficients for soluble and membrane proteins within cilia that are compatible with rapid exploration of the ciliary space in the absence of active transport. Our results indicate that large proteins require active transport for entry into cilia but not necessarily for movement inside cilia.
Introduction
The primary cilium organizes a variety of signaling cascades (G protein–coupled receptors [GPCRs], Hedgehog, and others) by dynamically concentrating signaling components (Goetz and Anderson, 2010; Nachury et al., 2010). A major unanswered question is how cilia concentrate signaling proteins despite topological continuities between plasma membrane and ciliary membrane and between cytosol and ciliary lumen. In the case of membrane proteins, lateral exchange between plasma and ciliary membranes is prevented by a septin-based diffusion barrier at the transition zone, a region at the base of cilia (Hu et al., 2010; Chih et al., 2012; Reiter et al., 2012), and by tethering of some plasma membrane proteins to the actin cytoskeleton (Francis et al., 2011). Conversely, it remains unclear whether a ciliary diffusion barrier exists for soluble proteins and, were it to exist, how it might operate.
On one hand, Kee et al. (2012) proposed that a size-dependent diffusion barrier restricts entry of cytosolic proteins into cilia. Specifically, after microinjection into the cytosol, fluorescent probes larger than ∼40 kD were not detectable in cilia. This study also suggested that nucleoporins (Nups) localize near the base of the cilium to restrict protein entry. On the other hand, Calvert et al. (2010; Najafi et al., 2012) have found no evidence for a diffusion barrier at the connecting cilium of rod photoreceptors, a structure analogous to the transition zone of primary cilia. First, the kinetics and energy independence of arrestin (47 kD) and transducin (27 kD) translocation through the connecting cilium are fully accounted for by free diffusion (Nair et al., 2005; Rosenzweig et al., 2007). Second, proteins 27–81 kD in size were found to cross the connecting cilium at the same rate (Najafi et al., 2012). Yet, at steady state, the larger proteins do not distribute evenly between inner and outer segments (equivalents of the cell body and the distal part of the cilium, respectively). Instead, the tight packing of disk membranes in the outer segment limits the volume accessible to large proteins, and these steric effects result in an apparent decrease in protein concentration in the outer segment. The absence of flux measurements by Kee et al. (2012) and their inability to resolve the base of cilia raise the possibility that steric effects may account for the observed size-dependent distribution of probes in primary cilia versus cytoplasm (Najafi and Calvert, 2012).
Here, we establish and validate a permeabilized cell assay to directly and quantitatively test whether soluble protein entry into mammalian primary cilia is gated by a diffusion barrier. Using this system, we find that primary cilia possess a size-dependent diffusion barrier that is mechanistically distinct from those found at the axon initial segment and the nuclear pore complex (NPC). We anticipate that our assay will be a powerful tool for mechanistic studies of trafficking to cilia and provide a basis for understanding how cilia regulate signal transduction.
Results
A permeabilized cell system for ciliary trafficking
In the course of observations on the lipid composition of primary cilia, we found that very low amounts of the cholesterol-dependent detergent digitonin selectively permeabilize the plasma membrane while leaving the ciliary membrane intact. This selective permeabilization is illustrated by the failure of antibodies against ciliary markers (e.g., acetylated α-tubulin and Arl13b) to stain primary cilia in digitonin-permeabilized cells despite strongly labeling cilia after permeabilization with 0.1% Triton X-100 (Fig. 1 A). In contrast, cytoplasmically exposed Nups at the nuclear envelope are readily detected by the mAb414 antibody after permeabilization with digitonin. Our results therefore suggest that after digitonin permeabilization, (a) the ciliary membrane remains intact, and (b) antibodies are not able to enter cilia from the cytosol.
To verify and extend these results, we stably expressed the somatostatin receptor 3 (Sstr3), a GPCR that is targeted to primary cilia, in IMCD3–Flp-In cells. For visualization, Sstr3 was tagged on its intracellular tail with GFP and the S tag. These cells were permeabilized with digitonin, incubated on ice with antibodies against GFP or the S tag, and then fixed and stained. Using this assay, we found that both anti-GFP and anti–S tag antibodies fail to access ciliary Sstr3-S-GFP after digitonin permeabilization but give strong ciliary staining in a conventional, Triton X-100–based immunofluorescence protocol (Figs. 1, B and C; and S1 B). For both antibodies, the plasma membrane pool of Sstr3-S-GFP is accessible in digitonin-permeabilized cells, thus confirming plasma membrane permeabilization and successful immunological detection of Sstr3-S-GFP under our experimental conditions (Fig. S1, A and B). Similar results were obtained when cells were permeabilized with the cholesterol-dependent pore-forming toxin perfringolysin O (PFO; Figs. 1 C and S1 A).
We next tested whether the size of the GFP binder affects its ability to enter primary cilia of digitonin-permeabilized cells. Here, we took advantage of a miniaturized GFP-binding protein (GBP; Kirchhofer et al., 2010). In contrast to the anti-GFP antibody, the 13.5-kD GBP readily gained access to ciliary Sstr3-GFP in digitonin or PFO-permeabilized cells (Fig. 1 D). A similar size-dependent entry of GBPs is observed when GFP is displayed inside cilia via fusion to the ciliary targeting signal of Pkhd1 (Pkhd1ICD; Fig. S1 C; Follit et al., 2010). Notably, a gradient of GBP signal from the base to the tip of the cilium was often seen in digitonin-permeabilized cells, consistent with GBP entering at the base of the cilium from the cytoplasmic compartment (Fig. 1, B, D, and E). This gradient of signal is not seen in GBP-stained cells after Triton X-100 permeabilization (Fig. 1, D and E), indicating that GBP entry at the base of cilia is specific to digitonin-permeabilized cells.
Permeabilization with PFO consistently yields the same results as digitonin, even though PFO is a pore-forming oligomeric cytolysin, whereas digitonin is a chemical detergent (Heuck et al., 2010). PFO-induced pores have a stereotypical diameter of ∼30 nm (Olofsson et al., 1993; Czajkowsky et al., 2004), whereas the diameter of an antibody is 10 nm, thus ruling out the possibility that PFO pore size accounts for size-dependent entry into cilia. Finally, two independent kinetic experiments (see Figs. 3 C and 8 C) confirm that GBPs entering cilia first appear in proximal regions near the base of the cilium before gaining access to more distal receptors. Collectively, the aforementioned observations confirm that digitonin and PFO selectively permeabilize the plasma membrane while leaving the ciliary membrane intact.
We further validated our permeabilized cell system by assaying two key features of primary cilia. First, FRAP of ciliary Sstr3-GFP in permeabilized cells confirms that the membrane diffusion barrier seen in live cells (Hu et al., 2010) remains functional, as unbleached receptors in the plasma membrane cannot move into the cilium within 5 min of bleaching (Fig. 2 A). As a control, we note that Sstr3 moves freely within the ciliary membrane, as assessed by half-cilium FRAP, and within the plasma membrane, as evidenced by rapid fluorescence recovery in bleached regions of the plasma membrane near the cilium (Fig. 2 A; Ye et al., 2013).
Second, we investigated intraflagellar transport (IFT), the motor-driven movement of clusters of IFT complexes along the axoneme (Pedersen and Rosenbaum, 2008). Live imaging of GFP-tagged IFT88 (an IFT-B complex subunit) revealed that IFT stops in permeabilized cells (Fig. 2 B and Ye et al., 2013). Because ATP is lost from permeabilized cells (Ye et al., 2013) and cytoplasmic factors may similarly be washed away, we reasoned that it should be possible to reactivate IFT by adding cellular extracts and energy as long as ciliary structures have not been profoundly altered by digitonin permeabilization. Using extracts from bovine retina, a rich source of ciliary transport factors (Jin et al., 2010), or from Xenopus laevis eggs, we observed a resumption of GFP-IFT88 movement inside nearly all tested cilia (Fig. 2 B). These results strongly suggest that cilia from permeabilized cells are functionally intact.
Several mechanisms may act in concert to account for the selective permeabilization of the plasma membrane. First, both PFO and digitonin interact with cholesterol and require cholesterol for permeabilization (Schulz, 1990; Flanagan et al., 2009). Thus, a lower amount of cholesterol in the ciliary membrane or a reduced accessibility of cholesterol (e.g., caused by association with lipid raft domains) may prevent permeabilization of the ciliary membrane. Similarly, low cholesterol levels in the nuclear envelope enabled the use of digitonin for nuclear transport assays (Adam et al., 1990). Using filipin to assess cholesterol levels in primary cilia of IMCD3 cells, we found no enrichment of cholesterol in ciliary versus plasma membrane (Fig. S2 A). Second, primary cilia may avoid permeabilization by low amounts of digitonin because of reduced accessibility to added buffer. The majority of cilia were found on the ventral surface of cells under our experimental conditions, as indicated by total internal reflection fluorescence (TIRF) microscopy (Fig. S2 B). Third, it is possible that cilia remain intact because limiting amounts of digitonin are more likely to permeabilize the much larger portion of the cell surface corresponding to the plasma membrane than the comparatively small fraction (∼1%) represented by a single cilium. Indeed, given that a small fraction of cells typically remain unpermeabilized after digitonin treatment, the probability of disrupting the ciliary membrane may be quite low.
Quantitative analysis of the ciliary diffusion barrier
To determine the size cutoff for passive entry into primary cilia, we generated a panel of six GBP fusions with proteins of increasing sizes and measured their hydrodynamic diameters by size exclusion chromatography (Fig. 3 B). Incubation of digitonin-permeabilized cells with dye-labeled GBP fusions revealed a progressive decrease in ciliary entry for proteins above 4.8 nm in diameter and no detectable ciliary signal for proteins above ∼8.6 nm in diameter (Fig. 3, A and B). The larger GBP fusions have the ability to bind Sstr3-GFP because all GBPs stain plasma membrane–localized Sstr3-GFP in digitonin-permeabilized cells (Fig. 3 A) and detect ciliary Sstr3-GFP in cells permeabilized with Triton X-100 (Fig. S3 A). To further test whether size is the primary factor governing ciliary entry, we generated a GBP fusion containing two copies of the IgG-binding Z domain derived from Staphylococcus aureus protein A (Nilsson et al., 1987). The ZZ-GBP fusion alone enters cilia of digitonin-permeabilized cells, but preincubation with rabbit IgG renders the enlarged ZZ-GBP:IgG unable to enter cilia (Fig. S3 B).
We next determined how GBP size affects the rate of entry using live confocal microscopy. As expected, the GBP signal inside cilia appears first at the base, before increasing in intensity and spreading distally toward the tip (Fig. 3 C and Video 1). To determine entry rates, we measured the time-dependent increase of the GBP signal inside cilia; only the signal at the most proximal segment of the cilium was analyzed to isolate entry from diffusional movement of GBPs within the cilium. The resulting data fit a first-order exponential equation, and the fitted entry rate decreases sharply with increasing GBP size (Fig. 3 D). These results echo those obtained in a previous study on rates of diffusion through the NPC (Mohr et al., 2009) and firmly establish that a size-dependent permeability barrier controls protein entry into cilia.
In vivo validation of the permeability barrier
To confirm that the basic properties of the permeability barrier are preserved in semipermeabilized cells, we sought to better characterize the permeability barrier in intact cells. Here, we leveraged a well-characterized system for chemical induction of heterodimerization that consists of the protein moieties FK506-binding protein (FKBP) and FKBP12-rapamycin binding domain (FRB) and the small molecule rapamycin (Komatsu et al., 2010). FRB was fused to Sstr3 together with the fluorescent protein TagRFP-T, and GFP was fused to FKBP (Fig. 4 A). In untreated cells, Sstr3-RFP-FRB is concentrated within the ciliary membrane, whereas GFP-FKBP is diffusely cytoplasmic (Fig. 4 B). We note a slight enrichment of GFP-FKBP near the basal body, in agreement with our observation that GFP alone is slightly enriched at the centrosome (Fig. S4 B). Upon addition of rapamycin, GFP-FKBP accumulated in cilia (Fig. 4 B), confirming that a small protein (molecular mass [MM] = 42 kD) rapidly diffuses into cilia. The accumulation of GFP-FKBP into cilia proceeded with very rapid kinetics, with the intensity of GFP-FKBP inside cilia plateauing within 45 s of rapamycin addition (Fig. 4 C). To assess whether a ciliary permeability barrier exists in intact cells, we fused GFP-FKBP to proteins of increasing size and found that the tetrameric fusions Tetra[FKBP–maltose-binding protein (MBP)-GFP] (MM = 351 kD) and LacZ-GFP-FKBP (MM = 634 kD) are excluded from cilia within our experimental time frame (Fig. 4 D). Importantly, these constructs are functional, as they can be rapidly recruited to the plasma membrane by Lyn-FRB (Fig. S3 D). We limited our analysis to <10 min, as FRAP analysis shows that the ciliary levels of Sstr3 start to increase 10 min after photobleaching, most likely because of the trafficking of newly synthesized Sstr3 to cilia. Meanwhile, the fusions 2×GFP-FKBP (69 kD), MBP-GFP-FKBP (84 kD), and luciferase (Luc)-GFP-FKBP (103 kD) were able to enter cilia within 6 min of rapamycin addition.
Although the molecular weight limit we observed in vivo is slightly greater than the value we found in vitro, the actual differences in hydrodynamic radii between Luc-GFP-FKBP (Stokes radius = 4.4 nm; Lin et al., 2013) and Luc-GBP (Stokes radius = 4.3 nm) are minimal. Furthermore, we surmised that transient transfection of GFP-FKBP fusions generates considerably higher protein concentrations than used in our in vitro system (55 nM) and that such a high concentration of the diffusion probe sufficiently increases the flux into cilia to generate a detectable signal. Consistent with this interpretation, low signals of Luc-GBP can be detected within cilia of semipermeabilized cells when the Luc-GBP concentration is increased sixfold and incubation time is increased fourfold (Fig. S3 C). Together with the IFT reactivation results, our in vivo system shows that the basic properties of primary cilia are conserved between live cells and permeabilized cells.
Unique properties of the ciliary diffusion barrier
Leveraging the ease of manipulation afforded by an in vitro system, we sought to characterize the molecular properties of the ciliary permeability barrier. At the axon initial segment, a barrier of densely cross-linked filaments separates the cell body from the axonal compartment (Winckler et al., 1999; Nakada et al., 2003; Song et al., 2009). Although the axonal permeability is disrupted by the actin poison Cytochalasin D (Song et al., 2009), treatment of IMCD3 cells with Cytochalasin D for 30 min before digitonin permeabilization did not permit ciliary entry of any GBPs larger than the size cutoff defined in Fig. 3 (Fig. 5 A).
Next, we examined the relationship between the ciliary diffusion barrier and the NPC barrier, which is formed by a hydrophobic hydrogel of FG repeats. Elements of the nuclear transport machinery have been implicated in ciliary trafficking, and it was recently proposed that Nups prevent passive protein entry into cilia (Dishinger et al., 2010; Fan et al., 2011; Kee et al., 2012). We tested the possibility that Nups restrict ciliary entry by applying reagents that modulate the NPC hydrogel to our assay. First, we examined 1,2-trans-cyclohexanediol, a “minidetergent” that disrupts hydrophobic interactions within the NPC hydrogel and thereby accelerates passive entry into HeLa cell nuclei (Fig. 5 B; Ribbeck and Görlich, 2002). However, cyclohexanediol treatment has no effect on the rate of MBP-GBP entry into IMCD3 cilia (Fig. 5, C and E). A similar lack of effect on ciliary entry is also seen with two potent inhibitors of entry into the nucleus: a truncated form of the nuclear transport factor Importin-β (Impβ(45–462); Kutay et al., 1997) and the lectin WGA, which binds and cross-links N-acetylglucosamine–modified Nups (Finlay et al., 1987). As previously reported (Mohr et al., 2009), both Impβ(45–462) and WGA block GFP entry into HeLa nuclei (Fig. 5 D, right), but neither reagent prevents ciliary entry of MBP-GBP (Fig. 5 D, left). Although we did observe a marked change in the distribution of MBP-GBP along the cilium in WGA-treated cells (Fig. 5 D and see Figs. 7 and 8), kinetic analysis confirmed that the rate of MBP-GBP arrival in the proximal cilium is unaffected by WGA (Figs. 5 E and S4 A).
As another means to assess the role of Nups in the ciliary diffusion barrier, we used the mAb414 antibody, which recognizes several FG-containing Nups and inhibits nuclear import (Davis and Blobel, 1986, 1987; Michaud and Goldfarb, 1992; Marshallsay and Lührmann, 1994). As seen for Impβ(45–462) and WGA, mAb414 did not inhibit MBP-GBP entry into cilia, despite prominently labeling the nuclear envelope of digitonin-permeabilized cells (Fig. 5 F).
Finally, we did not detect any mAb414-reactive Nups at or near the base of IMCD3 cell primary cilia by conventional immunofluorescence staining (Fig. 6 A). Furthermore, expression of three different Nups fused to GFP revealed no localization to the base of cilia, whereas targeting to the NPC was consistently robust (Fig. 6 B). These Nups belong to distinct subcomplexes of the central NPC scaffold (Hoelz et al., 2011) that had been proposed to localize to the base of the cilium (Kee et al., 2012). Although this difference in localization may be a result of differences in cell types, we note that GFP alone has a weak tendency to accumulate around centrosomes (Fig. S4 B and Fig. 4 B). Collectively, our data indicate that the ciliary diffusion barrier is distinct from the barriers present at the axon initial segment or the NPC.
WGA binds and immobilizes ciliary Sstr3-GFP
In digitonin-permeabilized cells incubated with WGA, we found that GBPs are restricted to the proximal region of the cilium (Fig. 5 D). In evaluating the mechanistic basis for this effect, we asked whether WGA might bind directly to Sstr3-GFP. Many membrane proteins contain N-acetylglucosamine and sialic acid moieties recognized by WGA, and Sstr3 is highly glycosylated and captured by WGA-derivatized resin (Händel et al., 1999; Nehring et al., 2000). Furthermore, because it is a multivalent lectin (Nagata and Burger, 1974), WGA may cross-link and immobilize Sstr3, as reported previously for red blood cell glycoproteins (Golan et al., 1986). Indeed, half-cilium FRAP of Sstr3-GFP confirms that WGA prevents the lateral movement of ciliary Sstr3 in permeabilized cells (Fig. 7, A and B), as it does in live cells (Ye et al., 2013).
To rule out the possibility that WGA alters the GBP progression within cilia through a Nup-related mechanism, we incubated cells with WGA and washed out excess lectin before digitonin addition. Under these conditions, WGA cannot bind to intracellular Nups, but a proximal bias of GBP staining is still observed (Fig. S5 A). We also found a proximal-shifted GBP distribution in cells treated with Maackia amurensis lectin, which recognizes sialic acid–containing glycans (Knibbs et al., 1991) and does not bind Nups (Fig. S5 B; Emig et al., 1995). Lastly, we found no effect of WGA on the GBP distribution when GFP is displayed in cilia via the Pkhd1ICD, which is not glycosylated (Fig. 7 C). Collectively, these data show that WGA reduces the mobility of Sstr3-GFP, which in turn influences the distribution of Sstr3-GFP:GBP complexes within cilia.
Properties of GBP diffusion to capture in the presence of WGA
How does immobilization of Sstr3-GFP modulate GBP capture inside cilia? The rapid progression of GBPs toward the cilium tip in the absence of WGA suggests that Sstr3 mobility accelerates GBP appearance in the distal portion of the cilium. In particular, diffusion within the ciliary membrane may enable Sstr3-GFP:GBP complexes that form in the proximal region of the cilium (near the site of GBP entry) to move toward the tip (Fig. 8 A, left). However, when Sstr3-GFP movement is blocked by WGA, GBPs entering cilia will be captured and immobilized at binding sites near the base, and only after these proximal GBP binding sites are saturated will GBPs gain access to more distal receptors (Fig. 8 A, right).
Consistent with this model, GBP staining in WGA-treated cells extends to increasingly distal regions of the cilium with longer staining times (Fig. 8 B). Furthermore, sequential incubation of WGA-treated cells first with Alexa Fluor 568–GBP and then with Alexa Fluor 647–GBP gives rise to two nonoverlapping domains of staining in which Alexa Fluor 568–GBP is found proximally and Alexa Fluor 647–GBP is found immediately distal to the first segment (Fig. 8 C, right). In cells not treated with WGA, free lateral movement of Sstr3-GFP:GBP complexes allows both GBPs to become evenly distributed along the entire cilium (Fig. 8 C, left).
To obtain a quantitative understanding of how Sstr3-GFP mobility impacts GBP entry and movement within cilia, we assessed the kinetics, concentration dependence, and spatial distribution of GBP staining in WGA-treated cells. For each condition, line intensity plots from base to tip were used to measure the length of the cilium stained by GBP (Fig. 9, A and B). From this analysis, we found that the square of the stained length increases linearly with time (Fig. 9 C). Furthermore, the rate of increase (i.e., slope) is proportional to the concentration of GBP added (Fig. 9 C). Similar results were obtained for thioredoxin (Trx)-GBP and MBP-GBP but with reduced rates (Fig. 9 F). This relationship between the length stained and time suggests that diffusion within cilia has a prominent role in our diffusion-to-capture assay, a possibility we examine by mathematical modeling.
Mathematical modeling of diffusion to capture within the cilium
We developed a mass transport model for the entry and movement of soluble GBPs within the cilium. Our model comprises concentration-dependent GBP entry into cilia, one-dimensional diffusion within cilia, and instantaneous capture by mobile or immobile Sstr3-GFP (Fig. 9 D). The motion of GBP is given by the standard diffusion equation
in which B is the concentration of free GBP, DB is its diffusion constant, and x is the position within the cilium. The flux of GBP entering through the barrier at the base of the cilium is modeled as a first-order kinetic process, occurring with a net rate proportional to the difference in GBP concentrations inside and outside the cilium, with rate constant k. At the base (x = 0), the flux of entering GBPs, J, is equal to the diffusional flux of GBPs, as derived from Fick’s law:
Here, B0 is the cytoplasmic concentration of GBP, assumed to be constant over time.
Our model assumes that the binding of GBP to Sstr3-GFP is irreversible and diffusion limited; indeed, GBP binds GFP with a Kd of ∼0.5 nM, and the half-life for complex dissociation is ∼80 min (Kirchhofer et al., 2010). Consequently, free Sstr3-GFP can never coexist with unbound GBP at the same location, and the cilium can be divided into two regions: a proximal region containing free GBP and a distal region containing free Sstr3-GFP. The boundary between these two regions serves as the site of Sstr3-GFP:GBP complex formation and moves from the base to the tip as more GBP molecules enter the cilium. The boundary position, L(t), is defined such that the concentration of free GBP at L(t) is 0:
When Sstr3-GFP is immobilized by WGA, free GBP arriving at the boundary binds to available receptors, gradually depleting free Sstr3-GFP and causing the boundary to move toward the cilium tip. In a given time interval, the number of GBP molecules arriving at the boundary equals the number of newly occupied Sstr3-GFPs, which can be calculated from the product of two measureable properties: the distance the boundary is displaced distally and the Sstr3-GFP concentration per unit length (S0):
We solve Eqs. 1, 2, 3, and 4 using standard techniques developed for dealing with analogous moving boundary diffusion problems (see Materials and methods; Crank, 2001). The solution exhibits two distinct behaviors of the boundary motion in different time regimes:
For short times, the motion is limited by entry of free GBP through the barrier at the base of the cilium. At long times, the boundary movement is dominated by diffusion of free GBP past the proximal region of occupied receptors to reach free Sstr3-GFPs. In this regimen, the square of the boundary position scales linearly with time and with the cytoplasmic GBP concentration L2(t) = 2[B0/S0]DBt. Note that this expression closely resembles canonical one-dimensional diffusion, <L2> = 2Dt, but with a scaling factor that accounts for the relative concentrations of GBP and Sstr3-GFP. The transition between the two regimes occurs at time t*, corresponding to the time needed for the concentration of free GBP at the cilium base to equilibrate with the cytoplasmic concentration.
To evaluate whether this model accounts for the observed rate and distance dependence of GBP staining, we used a recent tomographic study (Gilliam et al., 2012) to estimate the geometry of the cilium, and we determined the only other parameters not directly known: S0, the concentration of Sstr3-GFP in the ciliary membrane; k, the rate constant for GBP entry; and DB, the diffusion coefficient for GBP within the cilium. By comparing Sstr3-GFP fluorescence to a viral particle reference standard possessing 120 GFP molecules (Charpilienne et al., 2001; Joglekar et al., 2008), we found that S0 = 6,540 ± 2,340 Sstr3-GFP molecules per micrometer length of the cilium (Fig. S5 C). We used this density together with the kinetics of GBP entry measured in Fig. 3 D to determine the initial rate of GPB entry in molecules per second and thereby calculate k (see Materials and methods).
With these values known, we obtain DB by fitting the observed L2 versus t data (Fig. 9 C) to the model described by Eq. 5. The resulting diffusion coefficient for GBP within the cilium equals 7.3 ± 2.8 µm2/s. Interestingly, this value is comparable to cytoplasmic diffusion coefficients that have been measured for proteins of similar size (Elowitz et al., 1999; Luby-Phelps, 2000; Dix and Verkman, 2008), implying that the ciliary environment does not significantly slow down the motion of a small protein. It should also be noted that the transition time t* is ∼33 s, confirming that our experimental measurements fall in the long-time regime.
To assess the accuracy of our model, we plotted the time course of GBP progression into cilia of WGA-treated cells using Eq. 5 and the derived values for B0, S0, k, and DB and compared the results to those obtained experimentally. Importantly, we find a close agreement between the modeled and observed results (P = 0.984 using an F test of the null hypothesis that the model fits the data; Fig. 9, C and E). Thus, a basic mass transport system can accurately describe soluble protein entry and movement within cilia of permeabilized cells.
For the case with mobile Sstr3-GFP, we developed an analogous reaction diffusion model that includes diffusion of Sstr3 through the cilium (see Materials and methods). This model reveals that the approximate GBP boundary position is given by in which DS is the Sstr3 diffusion coefficient. Notably, when Sstr3-GFP can diffuse freely, the boundary motion is dictated by diffusion of Sstr3 rather than the motion of free GBP. Furthermore, comparison with Eq. 5 indicates that GBP progression along the length of the cilium is expected to be significantly faster in the case of mobile Sstr3-GFP even though the diffusion constant is smaller for the transmembrane protein (DS < DB). This effect, observed in Fig. 8 B, arises because the high concentration of Sstr3-GFP relative to GBP acts as a GBP sink that slows boundary movement in the WGA-treated case. Collectively, we find that a mathematical model accurately describes our experimental results and makes it possible to determine key physical properties of ciliary proteins.
Discussion
We have developed a permeabilized cell system that enables the quantitative analysis of protein entry into primary cilia. Several lines of evidence confirm that digitonin selectively permeabilizes the plasma membrane but not the ciliary membrane. First, antibodies cannot access antigens inside cilia of digitonin- or PFO-permeabilized cells but can recognize cytosolically exposed epitopes. Second, large (but not small) GBPs fail to access ciliary Sstr3-GFP after digitonin permeabilization despite readily recognizing plasma membrane Sstr3-GFP. Lastly, kinetic experiments confirm that small GBPs appear first near the base of cilia before gaining access to more distal regions.
By measuring flux into cilia of semipermeabilized cells, we resolve two conflicting studies (Kee et al., 2012; Najafi et al., 2012) and demonstrate that mammalian primary cilia possess a permeability barrier that restricts the entry of soluble proteins in a size-dependent manner. This barrier is evident from the reduction in entry rate as protein size increases, with entry not detectable for proteins with >9 nm hydrodynamic diameter (MM ∼100 kD). The agreement between our permeabilized cell system and in vivo studies by us (Fig. 4) and others (Lin et al., 2013) provides further support for the validity of our in vitro system. For soluble proteins found within cilia that are close to or above the passive diffusion limit (e.g., tubulin dimers, molecular motors, and Gli transcription factors), active transport mechanisms must be required for efficient ciliary targeting. Given that the transport of soluble proteins needs to overcome different constraints than those imposed on membranous cargoes, the identification of the targeting signals present on these ciliary proteins and of the transport machinery that recognizes these signals will be key goals for future studies. Furthermore, it will be critical to determine how trans-acting factors carry large cargoes across the ciliary permeability barrier. To gain mechanistic insights into these fundamental questions, a reconstituted transport system in semipermeabilized cells may provide a unique tool, enabling a range of experimental manipulations not possible in live cells. In particular, the reactivation of IFT observed upon addition of cellular extracts and ATP (Fig. 2 B) offers an appealing system for the mechanistic dissection of IFT-mediated trafficking.
Although the molecular basis of the ciliary permeability barrier awaits further characterization, our data demonstrate that this barrier is distinct from those present at the axon initial segment or the NPC. This latter finding is in contrast to a recent study indicating that Nups at the base of cilia restrict protein entry (Kee et al., 2012). Because Nup-regulated ciliary entry was only functionally supported for kinesin Kif17 (MM = 116 kD), it is possible that our divergent conclusions reflect mechanistic differences between passive diffusion and active transport. Nonetheless, the absence of detectable Nups at the transition zone in IMCD3 cells (Fig. 6) makes this possibility unlikely.
In addition to outstanding questions regarding the molecular basis of the ciliary diffusion barrier, the ultrastructural organization of the barrier is also unknown. The two most distinctive structural elements found at the base of cilia are the transition zone Y links and the transition fibers (Nachury et al., 2010; Reiter et al., 2012). However, neither of these elements is known to constrict to <40 nm. It is therefore likely that the ciliary diffusion barrier is not readily seen by electron microscopy and that new approaches will be needed to characterize its structural features.
Our permeabilized cell system also enabled the measure of key physical properties for ciliary proteins. Specifically, FRAP analysis and a quantitative model for GBP entry made it possible to determine diffusion coefficients for both transmembrane and soluble proteins inside mammalian primary cilia. Notably, because IFT movement stops in permeabilized cells, our measurements reflect diffusion alone and are not impacted by active transport. We find that diffusion within primary cilia is not significantly reduced relative to the cytoplasm, a result in agreement with studies on diffusion of macromolecules and metabolites in sperm flagella and frog photoreceptors (Takao and Kamimura, 2008; Calvert et al., 2010). Specifically, our observed ciliary diffusion coefficients (0.2 µm2/s for Sstr3-GFP and 7 µm2/s for GBP) are similar to reported values for protein diffusion in the plasma membrane (∼0.1–0.4 µm2/s) and the cytoplasm (∼1–20 µm2/s; Elowitz et al., 1999; Luby-Phelps, 2000; Calvert et al., 2001; Suzuki et al., 2005; Dix and Verkman, 2008; Chung et al., 2010; Jaqaman et al., 2011). Interestingly, these diffusion coefficients indicate that soluble and membrane proteins can travel the roughly 4-µm distance from ciliary base to tip in <5 and ∼30 s, respectively, time scales comparable to IFT-mediated transport (∼0.5–1.5 µm/s; unpublished data; Follit et al., 2006; Tran et al., 2008; Besschetnova et al., 2009).
It is therefore possible that IFT complexes are not strictly required for rapid protein movement within cilia. For example, consider a scenario in which tubulin incorporation at the tip of the axoneme is limited solely by the rate of diffusional arrival of tubulin from a pool near the base (analogous to GBP capture by immobilized Sstr3-GFP). In this case, at steady-state, there would be a linear concentration profile of free tubulin from base to tip, and the rate of arrival at the cilium tip would be proportional to the tubulin concentration at the base. Additionally, the arrival rate would be inversely proportional to cilium length, a relationship also predicted by the balance point model for ciliary length regulation proposed by Marshall et al. (2005). Furthermore, a tubulin concentration near the base of ∼0.5 µM (with Dtubulin = 1 µm2/s) would allow diffusional flux to sustain a cilium of length of ∼5 µm and offset continuous turnover of axonemal tubulin (Marshall et al., 2005). Such a model is also consistent with the recent finding that soluble tubulin levels regulate cilium length (Sharma et al., 2011).
Although the regulation of axoneme growth is likely to be considerably more complex in vivo and may use IFT (Hao et al., 2011), our observations nonetheless underscore the need for a better understanding of which ciliary processes are dependent on active transport and which are supported by simple diffusion. The importance of diffusion is highlighted by the recent finding that, in specific genetic contexts, seemingly null mutations in dynein 1b permit normal Hedgehog signaling and cilium assembly (Ocbina et al., 2011) and by the observation that normal cilium morphology is maintained upon acute inactivation of dynein 1b (Engel et al., 2012). In the presence of such dynein mutations, protein movement from tip to base of cilia—to recycle damaged proteins or transmit Hedgehog signals to the cell body—is likely to rely on passive diffusion. Meanwhile, anterograde IFT is indispensable for cilium assembly, likely because of the large amounts of materials that must be delivered to the tip of cilia during axoneme elongation (Ishikawa and Marshall, 2011). IFT may also carry out alternative functions such as promoting the active entry of large proteins into cilia (Pedersen and Rosenbaum, 2008) or participating in scaffolding signal transduction cascades, as suggested by the work of the Snell group on IFT-dependent signaling in flagella of Chlamydomonas reinhardtii (Wang et al., 2006). In future studies, the permeabilized cell system presented here will likely provide a valuable tool for dissecting active transport processes and their contributions to ciliary function.
Materials and methods
Antibodies and reagents
The following antibodies were used for immunofluorescence staining: goat anti-GFP (#600-101-215; Rockland Immunochemicals), mouse anti–S tag (#71549; EMD), mouse antimyc (#sc-40; Santa Cruz Biotechnology, Inc.), rabbit anti-Arl13b (gift from T. Caspary; Emory University, Atlanta, GA), mouse anti–acetylated α-tubulin (#T6793; Sigma-Aldrich), mouse anti–γ-tubulin (#T6557; Sigma-Aldrich), mouse anti-Nup mAb414 (#MMS-120R; Covance), and rabbit antininein (gift from M. Bornens, Institut Curie, UMR144 du Centre National de la Recherche Scientifique, Paris, France). Alexa Fluor 546–labeled phalloidin was purchased from Invitrogen (#A22283). WGA (#L9640 and #L5142), Cytochalasin D (#C8273), cyclohexanediol (#141712), creatine kinase (C3755), phosphocreatine (P7936), and ATP (A7699) were purchased from Sigma-Aldrich. High-purity digitonin (#300410) and rapamycin (#553211) were purchased from EMD Millipore.
Cell culture
IMCD3 cells were cultured in DMEM/F12 medium (Life Technologies) containing 10% FBS (HyClone) at 37°C in 5% CO2. To induce ciliogenesis, cells were grown in serum starvation medium (0.2% FBS) for 24 h. Clonal stable cell lines expressing Sstr3-S-GFP, GFP-S-IFT88, GFP-Pkhd1ICD, and Sstr3-TagRFP-T-FRB under control of the EF1-α promoter were generated using the Flp-In system (Invitrogen) as described previously (Jin et al., 2010). Plasmids encoding TagRFP-T, Sstr3, and the FRB/FKBP domains were provided by R. Tsien (University of California, San Diego, La Jolla, CA), K. Mykytyn (The Ohio State University, Columbus, OH), and T. Meyer (Stanford University, Stanford, CA), respectively. Where indicated, cells were transfected using Lipofectamine 2000 (Life Technologies) with tagged Nups (gift from M. Hetzer, Salk Institute for Biological Studies, La Jolla, CA), the pericentrin-AKAP450 centrosomal targeting domain from pericentrin (gift from S. Munro, Medical Research Council Laboratory of Molecular Biology, Cambridge, England, UK; Gillingham and Munro, 2000) or with FKBP fusions, including the tetramerizing coiled coil from Gcn4 (gift from K.C. Garcia, Stanford University, Stanford, CA).
Cell extracts
Cytosolic extract from bovine retina was prepared as previously described (Jin et al., 2010). In brief, frozen retinas (W L Lawson Co.) were thawed in breaking buffer (25 mM Hepes, pH 7.4, 150 mM KOAc, 210 mM sucrose, 2 mM MgCl2, 2.5 mM CaCl2, and 5 mM EGTA) supplemented with protease inhibitors. After dounce homogenization and removal of cell debris by centrifugation at 1,000 g, the lysate was centrifuged for 1 h in a rotor (TLA100.3; Beckman Coulter) at 100,000 g. The resulting supernatant was supplemented before use with an ATP-regenerating system containing 1 mM ATP, 5 mM phosphocreatine, and 40 U/ml rabbit muscle creatine kinase. Cytosolic extract from Xenopus eggs (low-speed supernatant) was prepared as described previously (Moree et al., 2011). In brief, Xenopus eggs were washed in MMR buffer (5 mM Hepes, 0.1 mM EDTA, 100 mM NaCl, 2 mM KCl, 1 mM MgCl2, and 2 mM CaCl2, pH 7.8) and dejellied in MMR with 2% l-cysteine. Eggs were then washed in CSF-XB buffer (100 mM KCl, 50 mM sucrose, 2 mM MgCl2, 0.1 mM CaCl2, 10 mM Hepes, and 5 mM EGTA, pH 7.7) and washed again in CSF-XB with protease inhibitors. Eggs were placed in a 13 × 51–mm ultraclear tube (Beckman Coulter) and packed by centrifugation. Excess buffer was removed before centrifugation for 15 min at 10,000 rpm in a rotor (SW55 Ti; Beckman Coulter). The soluble cytoplasmic material was then removed and supplemented with 10 mg/ml Cytochalasin D and an energy regenerating system.
Recombinant protein expression, purification, and labeling
PFO with the C459A mutation (gift from A. Johnson, Texas A&M University, College Station, TX) was expressed in Rosetta2(DE3)-pLysS cells (EMD Millipore) and purified as previously described (Shepard et al., 1998). His-tagged EGFP and His-tagged Impβ(45–462) (gift from D. Görlich, Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany) were expressed from pET vectors (EMD Millipore) in Rosetta2(DE3)-pLysS cells. Protein expression was induced with 200 µM IPTG, and cells were grown overnight at 19°C in 2YT medium. EGFP-expressing cells were lysed by sonication in 4XT buffer (80 mM Tris and 800 mM NaCl, pH 7.4) before dilution 1:1 with water. Cells expressing Impβ(45–462) were lysed in 30 mM sodium phosphate and 300 mM NaCl, pH 7.2. Proteins were purified on Ni–nitrilotriacetic acid agarose (QIAGEN) and eluted from resin with 300 mM imidazole.
DNA encoding the camelid GBP (Kirchhofer et al., 2010) was prepared by gene synthesis and cloned into a pET expression vector with a C-terminal 6×His tag. GBP fusions (except for MBP-GBP) were generated by inserting genes of interest upstream of GBP in this vector. MBP-GBP was generated by cloning GBP-6×His into the pMal-C2x MBP-tagging plasmid (New England Biolabs, Inc.). All GBP fusions were expressed in Rosetta2(DE3)-pLysS cells and grown in ZYM-5052 autoinduction medium for 18 h at 22°C. Cells were harvested by centrifugation and lysed by sonication in buffer A (20 mM sodium phosphate, pH 7.2, 500 mM NaCl, 10 mM imidazole, 0.5 mM DTT, 1% NP-40, and 10 µg/ml each of leupeptin, pepstatin A, bestatin, aprotinin, AEBSF, and E-60). Lysates were clarified by centrifugation and bound to Ni–nitrilotriacetic acid agarose. After extensive washing, GBPs were eluted in buffer A with 300 mM imidazole and loaded onto a Superdex 200 column for further purification and Stokes radius measurement. Size exclusion chromatography was performed in buffer B (20 mM Hepes, pH 7.7, 150 mM NaCl, and 0.5 mM DTT) on a purifier (Äkta; GE Healthcare). Calibration standards for Stokes radius determination were purchased from Sigma-Aldrich (blue dextran, thyroglobulin, β-amylase, BSA, and cytochrome c). Peak fractions were pooled and labeled with Alexa Fluor 647 or Alexa Fluor 568 succinimidyl esters according to the manufacturer’s protocol (Life Technologies), and excess dye was separated from labeled protein with columns (NAP-5; GE Healthcare).
Semipermeabilization with digitonin and PFO
Coverslips with serum-starved IMCD3 cells were first placed on an ice-chilled metal block and washed twice with cold assay buffer (20 mM Hepes, pH 7.4, 115 mM KOAc, 1 mM MgCl2, and 1 mM EGTA). To induce permeabilization, coverslips were incubated for 7 min in cold assay buffer supplemented with 30 µg/ml digitonin and protease inhibitors (10 µg/ml each of leupeptin, pepstatin A, bestatin, aprotinin, AEBSF, and E-60). Alternatively, coverslips were incubated on ice for 10 min in assay buffer supplemented with 50 µg/ml PFO. Unbound PFO was removed by washing with cold assay buffer, and pore formation was then induced by incubating coverslips for 15 min at 37°C in assay buffer with protease inhibitors.
After permeabilization, coverslips were washed twice with assay buffer and incubated with GBPs and/or antibodies diluted in assay buffer supplemented with 0.1% BSA. Unless otherwise indicated, incubations with GFP binders were performed for 10 min at room temperature. All GBPs were used at 55 nM for endpoint assays and at 110 nM for time-lapse imaging; antibodies were used at 2–4 µg/ml. Where indicated, rabbit IgG at 50 µg/ml was preincubated with ZZ-GBP for 10 min. For WGA treatment, cells were incubated with WGA at 75 µg/ml for 7 min immediately before GBP staining; in these cases, WGA was also included in assay buffer during staining. After staining, coverslips were washed twice with assay buffer, fixed, and processed for immunofluorescence microscopy.
Immunofluorescence microscopy
Cells grown on acid-washed 12-mm #1.5 coverslips were fixed in assay buffer containing 4% paraformaldehyde for 10 min at room temperature. For γ-tubulin staining, cells were fixed for 5 min with paraformaldehyde, extracted in −20°C methanol for 5 min, and rehydrated in PBS. After cell fixation, cells were permeabilized in PBS with 0.1% Triton X-100 for 10 min, washed twice with PBS, and blocked for 30 min in PBS with 3% BSA and 5% normal donkey serum. Cells were then incubated with primary antibodies diluted in PBS + 3% BSA for 1 h at room temperature, washed five times with PBS + 3% BSA, and incubated with Cy3- and/or Alexa Fluor 647–labeled secondary antibodies (GBPs labeled with Alexa Fluor 647 or Alexa Fluor 568) diluted in PBS + 3% BSA for 30 min. Finally, coverslips were further washed, stained with Hoechst DNA dye, and mounted on slides in 80% glycerol and 10 mM Tris, pH 7.4.
Fixed cells were imaged at room temperature at 350-nm z intervals on a microscope (Axio Imager.M1; Carl Zeiss). Epifluorescent illumination was provided by a light source (Lambda XL; Sutter Instrument) through a 63×, 1.4 NA Plan Apochromat objective. Images were captured with a camera (CoolSNAP HQ2; Photometrics) using SlideBook software (Intelligent Imaging Innovations).
Live-cell imaging and photobleaching
For live-imaging microscopy, cells were cultured on 25-mm circular acid-washed #1.5 glass coverslips and imaged at 37°C. Permeabilized cells and viral particles were imaged in assay buffer; live cells were imaged in live-imaging medium (phenol red-free DMEM/F12 medium supplemented with 25 mM Hepes, pH 7.4).
Time-lapse imaging of GBP entry and imaging of GFP-tagged viral particles were performed through a 63×, 1.4 NA Plan Apochromat objective on an inverted microscope (TE2000; Nikon). Illumination at 488 and 647 nm was provided by a laser (Innova 70C-Spectrum; Coherent), and images were collected through a spinning-disc confocal scanner (CSU10; Yokogawa Corporation of America) and acquired with a charge-coupled device camera (Cascade 512B; Photometrics) using MetaMorph software (Molecular Devices). To minimize background signal from GBPs binding to Sstr3-GFP in the plasma membrane, permeabilized cells were preincubated with unlabeled LacZ-GBP. Images in the GFP and Alexa Fluor 647–GBP channels were acquired in five z planes with 450-nm spacing at the indicated time intervals. Cilia of permeabilized cells were identified from their GFP fluorescence, and time-lapse videos were started upon addition of GBP at 110 nM.
For analysis of ciliary entry in intact cells, GFP-FKBP fusion proteins were expressed from an EF1-α promoter via transient transfection of IMCD3 cells stably expressing Sstr3-TagRFP-T-FRB. Time-lapse imaging was performed after addition of 1 µM rapamycin; alternatively, cells were fixed after 10 min and processed for immunofluorescent staining as described previously. Time-lapse imaging of rapamycin-dependent ciliary entry and of GFP-IFT88 movement was performed on an inverted microscope (DeltaVision Elite; Applied Precision) equipped with a solid state light source (InsightSSI) and a laser module (X4) controlling a 50-mW 488-nm laser (Applied Precision). Images were captured with CoolSNAP HQ2 or sCMOS cameras (Applied Precision) using softWoRx software (Applied Precision). TIRF images of GFP-FKBP fluorescence were acquired through a 60×, 1.49 NA Apochromat normal TIRF objective (Olympus) at 6–9-s intervals for 6 min after rapamycin addition. Additional epifluorescence and TIRF images of Sstr3-RFP-FRB and GFP-FKBP were acquired immediately before and after time-lapse imaging. For analysis of GFP-IFT88 movement, a single z plane was imaged through a 60×, 1.4 NA objective (Olympus) at 1–2 frames/s for 30 s.
FRAP and TIRF microscopy of IMCD3 Sstr3-GFP cells were performed on a microscope (Axiovert 200M; Carl Zeiss). Epifluorescence excitation through a 63×, 1.4 NA oil objective (Olympus) was provided by a Xenon lamp (DG4 300W; Sutter Instrument). TIRF excitation was provided by a 473-nm solid-state laser (CrystaLaser) and a TIRF slider (Carl Zeiss) through a 100× Plan Fluor 1.45 NA objective (Olympus). Photobleaching was performed with a FRAP laser system (MicroPoint; Photonic Instruments). Images were acquired with a camera (CoolSNAP HQ), and the system was controlled with SlideBook software. In a typical bleaching experiment, cells were imaged in a single focal plane at one frame/s for 3 s before bleaching and 60 s after bleaching. For WGA-treated cells, the interval was 0.2 frames/s.
Image analysis
All data analysis described in this section was performed using Matlab (MathWorks) and ImageJ (National Institutes of Health) software. For measurement of GBP entry rates, masks corresponding to cilia were first defined based on Sstr3-GFP signal relative to background. Minor drift in the focal plane during imaging was effectively corrected by selecting at each time point the z section with the greatest ciliary GFP signal. To measure the GBP signal at the base, the GBP intensity was computed for all 2 × 2–pixel areas near the ciliary base, and the mean of the four brightest 2 × 2 areas was used for subsequent data analysis. This GBP intensity, after subtracting the background measured from pixels neighboring the cilium, was determined at each time point and used for fitting to a single exponential equation with rate k. Fitted rates were determined separately for each cilium and averaged to determine the mean and standard deviation of the entry rate for each GBP. We note that similar results were obtained with or without photobleaching correction.
For measurement of GBP staining length in WGA-treated cells, line segments along cilia were defined in ImageJ. Intensity line profiles for each channel along the length of cilia were measured in ImageJ and aligned according to the peak of the basal body staining. The points at which the observed intensity decreased to <25% of the maximal observed intensity were used to define the start and end points of GFP and GBP signals. The length from start to end point along the cilium was measured for each cilium and averaged to obtain the mean length; normalized intensity profiles were averaged to produce a representative line profile for each condition.
For quantification of FRAP data, cilia were first identified from the prebleach Sstr3-GFP signal. For each ciliary pixel, the relative intensity before and after bleaching was determined and used to define the boundary between the bleached and unbleached region (with the boundary typically corresponding to pixels with ∼30% of prebleach intensity). Regions roughly 1 µm proximal to and distal to the boundary were used for measurement of unbleached and bleached pixel intensities. The mean signal in the bleached and unbleached regions after background subtraction was determined at each time point. The mean pixel intensity throughout the cilium was also calculated to determine the expected final intensity in the bleached region (at infinite time and with 100% receptor mobile fraction). To generate mean recovery curves, measured intensities were normalized for each cilium such that the intensity of the bleached region equals 1.0 before bleaching and equals 0.0 immediately after bleaching. The receptor mobile fraction was calculated from the observed magnitude of fluorescence recovery relative to that expected if Sstr3-GFP fluorescence were spread evenly along the cilium. The fraction of molecules within the bleached region that were successfully photobleached was determined from the relative fluorescence in the bleached versus unbleached regions immediately after bleaching. To obtain diffusion coefficients, recovery curves were first normalized such that the mean postbleach intensity in the whole cilium equals 1.0 and the background fluorescence equals 0.0. The cumulative density of bleached molecules between positions 0 and x at time t is given by
in which the region from 0 to a is bleached, L is the length of the cilium, and D is the diffusion coefficient. We fit the measured fluorescence intensity in a narrow segment of the cilium as a function of time after photobleaching to a model of one-dimensional diffusion in a bounded region. If fb is the fraction of molecules bleached within the bleached region, and fm is the fraction of molecules that are mobile, the fractional intensity measured between positions a1 and a2 at time t is given by
This functional form is used in a least-squares fit for the diffusion coefficient D. Fitted values for D were determined for each cilium analyzed and averaged to obtain the mean diffusion coefficient for Sstr3-GFP.
To determine the concentration of Sstr3-GFP molecules in the ciliary membrane, GFP-tagged viral particles were imaged under the same imaging conditions as used for measuring rates of GBP entry in permeabilized cells. Each particle’s intensity was measured in a summed z-stack projection. After background subtraction, the mean particle intensity was used to determine the fluorescence intensity of a single GFP molecule under the imaging conditions used. Similarly, the summed z-projection intensities for cilia were measured, and the Sstr3-GFP intensity per unit length was used to determine the mean number of Sstr3-GFP molecules per micrometer length of the ciliary membrane.
Mathematical modeling of diffusion to capture
Immobile Sstr3-GFP.
The transport of GBP along the cilium where the Sstr3-GFP receptors have been immobilized with WGA is described by the moving-boundary diffusion problem encompassed by Eqs. 1, 2, 3, and 4, together with the initial condition L(0) = 0. We note that our model assumes rapid equilibration of GBP concentration across the cilium cross section, so that only transport along the length of the cilium is considered. The solution to these differential equations is derived analogously to the moving boundary problems described by Crank (2001) and can be expressed as
in which Erf is the error function, Erfc is the complementary error function, ξ satisfies
and t* satisfies
At short times (t << t*), the front position at which the concentration of GBP hits zero moves as L(t) = t/2t*. At long times (t >> t*), this front moves as
In the limiting case in which the Sstr3-GFP receptors are much more concentrated than the GBP (S0 >> B0), this solution reduces to the form given in Eq. 5.
Mobile Sstr3-GFP.
If the Sstr3-GFP receptors are allowed to diffuse over the ciliary membrane, there will be an initial period during which all entering GBP is absorbed by receptors at the base, and there is no free GBP within the cilium. This period ends at time when all the free Sstr3-GFP at the base is depleted. As explained in this section, we find that our measurements fall exclusively within the early time regimen (t << ), and thus, we investigated GBP capture under these conditions. The concentration profile S(x,t) of free Sstr3-GFP in this case must satisfy the diffusion equation and initial condition ∂S/∂t = DS(∂2S/∂x2) and S(x,t = 0) = S0. Furthermore, at the entry point at the base of the cilium, the flux of disappearing free Sstr3-GFP should equal the rate of incoming GBP molecules DS(∂S/∂x) = kB0. The solution to these three equations is given by
The transition time is then found by setting yielding
In the early time regimen, all bound GFP–GBP complexes are produced at the origin, with constant rate kB0. The bound complexes spread diffusively from the origin, with a diffusion constant that is assumed to be equal to that of free Sstr3-GFP (DS). The density profile of the complexes is then given by
We define the approximate distance along the cilium over which the bound complexes have spread as the position such that 75% of the bound complex density lies between the base and This position can then be expressed as in which a ≈ 1.08 satisfies the equation
Evaluation of model parameters
The density of Sstr3-GFP is known from the surface density of Sstr3 receptors determined in Fig. S5 C, yielding S0 = 6.54 receptors/nm (expressed subsequently as nm−1) and DS = 0.21 µm2/s, as measured in Fig. 2 A. To determine the rate constant for GBP entry, k, we used the relative rate of GBP entry measured in Fig. 3 D. This relative rate, k* = 0.031 s−1, can be used together with Nbase, the estimated final number of GBP molecules near the base of the cilium, to determine V0, the initial rate of GBP entry in molecules per second. V0 is in turn related to k as follows: V0 = k* × Nbase = k × B0, yielding k = 28,000 nm/s when Nbase is calculated for a 500-nm segment at the cilium base containing 500 × S0 = 3,270 GBP binding sites.
We note that in this one-dimensional model, the variables B0 and S0 correspond to the linear densities of GBP and Sstr3-GFP, respectively. To calculate these linear densities from the bulk concentration requires an estimate of the cilium geometry (Gilliam et al., 2012). We treat the ciliary lumen as an annular region with an outer radius, R = 151 nm, and an inner radius, r = 78 nm (corresponding to the region excluded by the microtubule core). The density of GBP is then given by B0 (in nm−1) = 0.6 × π(R2 − r2)[B]0, in which [B]0 is the bulk solution concentration expressed in molars. We note that treating the ciliary lumen as a hollow tube with R = 151 nm that contains axonemal microtubules (modeled as nine 22-nm-diameter A tubules and nine 24-nm-diameter B tubules) yields similar results (Dannulus = 7.3 µm2/s; Dtube = 5.9 µm2/s).
To extract an estimate for the diffusion constant of free GBP within the cilium, we fit all data in Fig. 9 C to the full expression for the boundary position over time, obtained by setting Eq. 6 equal to 0. The resulting diffusion constant DB = 7.3 µm2/s was used to generate the theoretical curves in Fig. 9 E. An F test was used to assess goodness of fit of the modeled curves to the observed results. Note that the standard deviation for DB of 2.8 µm2/s also includes error in S0 and B0 (the latter caused by uncertainty in reported values for cilium geometry; Gilliam et al., 2012). From DB = 7.3 µm2/s, we find that t* = 33 s for [B]0 = 27.5 nM; t* = 8 s for [B]0 = 110 nM; = 3 h for [B]0 = 27.5 nM; and = 11 min for [B]0 = 110 nM.
Online supplemental material
Fig. S1 shows additional validation of the semipermeabilized cell assay. Fig. S2 shows filipin staining and TIRF imaging of IMCD3 cell cilia. Fig. S3 shows size-dependent cilia entry of GBP fusion proteins and size-independent access of these proteins as well as large FKBP fusions to the plasma membrane. Fig. S4 shows kinetics of MBP-GBP cilia entry in WGA-treated cells relative to untreated cells as well as weak accumulation of GFP alone near the cilium base. Fig. S5 provides further evidence that WGA cross-links and immobilizes Sstr3-GFP and shows absolute quantification of Sstr3-GFP concentration in cilia by comparison to viral particles bearing exactly 120 GFP molecules. Video 1 shows kinetics of GBP, Trx-GBP, and MBP-GBP capture in cilia of digitonin-permeabilized cells.
Acknowledgments
We thank A. Straight and W.J. Nelson for generously sharing microscopes, J. Liang for help with video editing and development of the in vivo assay for diffusion into cilia, and A. Amodeo for preparation of the Xenopus egg extract. We thank M. Bornens, T. Caspary, K.C. Garcia, D. Görlich, S. Munro, M. Hetzer, A. Johnson, T. Meyer, K. Mykytyn, and R. Tsien for sharing antibodies and plasmids, P. Miller for help with initial characterization of GBPs, H. Jin for generating the GFP-IFT88 cell line, C. Fuller and F. Ye for assistance with live imaging, and members of the Nachury laboratory for helpful discussions.
This work was supported by grants to M.V. Nachury from the National Institutes of Health (GM089933) and the March of Dimes (1-FY11-517) and by the Connie and Bob Lurie Fellowship of the Damon Runyon Cancer Research Foundation to D.K. Breslow (DRG 2087-11).
References
- FRB
FKBP12-rapamycin binding domain
- FKBP
FK506-binding protein
- GBP
GFP-binding protein
- GPCR
G protein–coupled receptor
- IFT
intraflagellar transport
- Luc
luciferase
- MBP
maltose-binding protein
- MM
molecular mass
- NPC
nuclear pore complex
- Nup
nucleoporin
- PFO
perfringolysin O
- Sstr3
somatostatin receptor 3
- TIRF
total internal reflection fluorescence
- Trx
thioredoxin