Singlemolecule tracking in live cells reveals distinct targetsearch strategies of transcription factors in the nucleus
Abstract
Gene regulation relies on transcription factors (TFs) exploring the nucleus searching their targets. So far, most studies have focused on how fast TFs diffuse, underestimating the role of nuclear architecture. We implemented a singlemolecule tracking assay to determine TFs dynamics. We found that cMyc is a global explorer of the nucleus. In contrast, the positive transcription elongation factor PTEFb is a local explorer that oversamples its environment. Consequently, each cMyc molecule is equally available for all nuclear sites while PTEFb reaches its targets in a positiondependent manner. Our observations are consistent with a model in which the exploration geometry of TFs is restrained by their interactions with nuclear structures and not by exclusion. The geometrycontrolled kinetics of TFs targetsearch illustrates the influence of nuclear architecture on gene regulation, and has strong implications on how proteins react in the nucleus and how their function can be regulated in space and time.
https://doi.org/10.7554/eLife.02230.001eLife digest
Transcription factors are proteins that control the expression of genes in the nucleus, and they do this by binding to other proteins or DNA. First, however, these regulatory proteins need to overcome the challenge of finding their targets in the nucleus, which is crowded with other proteins and DNA.
Much research to date has focused on measuring how fast proteins can diffuse and spread out throughout the nucleus. However these measurements only make sense if these proteins have access to the same space within the nucleus.
Now, Izeddin, Récamier et al. have developed a new technique to track single protein molecules in the nucleus of mammalian cells. A transcription factor called cMyc and another protein called PTEFb were tracked and while they diffused at similar rates, they ‘explored’ the space inside the nucleus in very different ways.
Izeddin, Récamier et al. found that cMyc explores the nucleus in a socalled ‘noncompact’ manner: this means that it can move almost everywhere inside the nucleus, and has an equal chance of reaching any target regardless of its position in this space. PTEFb, on the other hand, searches the nucleus in a ‘compact’ way. This means that it is constrained to follow a specific path through the nucleus and is therefore guided to its potential targets.
Izeddin, Récamier et al. explain that the different ‘search strategies’ used by these two proteins influence how long it takes them to find their targets and how far they can travel in a given time. These findings, together with information about where and when different proteins interact in the nucleus, will be essential to understand how the organization of the genome within the nucleus can control the expression of genes. The next challenge will now be to uncover what determines a protein's search strategy in the nucleus, as well as the potential ways that this strategy might be regulated.
https://doi.org/10.7554/eLife.02230.002Introduction
The nucleus is a complex environment where biochemical reactions are spatially organized in an interaction network devoted to transcription, replication, or repair of the genome (Misteli, 2001). Molecular interactions relevant to gene regulation involve transcription factors (TFs) that bind to specific DNA regulatory sequences or other components of the transcriptional machinery. In order to find their targets, TFs diffuse within the seemingly noncompartmentalized yet highly organized nuclear volume. Since the kinetics of a reaction can be largely determined by the mobility characteristics of the reactants (Rice, 1985; Shlesinger and Zaslavsky, 1993), the targetsearch strategy of TFs is a key element to understand the dynamics of transcriptional activity and regulation.
Over the past decade, the nuclear dynamics of TFs has become an important topic of research and has been investigated with a variety of imaging and biochemical approaches. Overall, these studies have emphasized the high mobility of nuclear factors, which results from a combination of diffusive motion and transient specific and nonspecific interactions with chromatin (Darzacq et al., 2009; Mueller et al., 2010; Normanno et al., 2012). These transient interactions are essential to ensure a fine regulation of binding site occupancy—by competition or by altering the TF concentration—but must also be persistent enough to enable the assembly of multicomponent complexes (Dundr, 2002; Darzacq and Singer, 2008; Gorski et al., 2008; Cisse et al., 2013).
In parallel to the experimental evidence of the fast diffusive motion of nuclear factors, our understanding of the intranuclear space has evolved from a homogeneous environment to an organelle where spatial arrangement among genes and regulatory sequences play an important role in transcriptional control (Heard and Bickmore, 2007). The nucleus of eukaryotes displays a hierarchy of organized structures (Gibcus and Dekker, 2013) and is often referred to as a crowded environment.
How crowding influences transport properties of macromolecules and organelles in the cell is a fundamental question in quantitative molecular biology. While a restriction of the available space for diffusion can slow down transport processes, it can also channel molecules towards their targets increasing their chance to meet interacting partners. A widespread observation in quantitative cell biology is that the diffusion of molecules is anomalous, often attributed to crowding in the nucleoplasm, cytoplasm, or in the membranes of the cell (Höfling and Franosch, 2013). An open debate remains on how to determine whether diffusion is anomalous or normal (Malchus and Weiss, 2009; Saxton, 2012), and the mechanisms behind anomalous diffusion (Saxton, 2007). The answer to these questions bears important consequences for the understanding of the biochemical reactions of the cell.
The problem of diffusing molecules in nonhomogenous media has been investigated in different fields. Following the seminal work of de Gennes (1982a), (1982b) in polymer physics, the study of diffusivity of particles and their reactivity has been generalized to random or disordered media (Kopelman, 1986; Lindenberg et al., 1991). These works have set a framework to interpret the mobility of macromolecular complexes in the cell, and recently in terms of kinetics of biochemical reactions (Condamin et al., 2007). Experimental evidence has also been found, showing the influence of the glasslike properties of the bacterial cytoplasm in the molecular dynamics of intracellular processes (Parry et al., 2014). These studies demonstrate that the geometry of the medium in which diffusion takes place has important repercussions for the search kinetics of molecules. The notion of compact and noncompact exploration was introduced by de Gennes (1982a) in the context of dense polymers and describes two fundamental types of diffusive behavior. While a noncompact explorer leaves a significant number of available sites unvisited, a compact explorer performs a redundant exploration of the space. In chemistry, the influence of compactness is well established to describe dimensional effects on reaction rates (Kopelman, 1986).
In this study, we aim to elucidate the existence of different types of mobility of TFs in the eukaryotic nucleus, as well as the principles governing nuclear exploration of factors relevant to transcriptional control. To this end, we used singlemolecule (SM) imaging to address the relationship between the nuclear geometry and the search dynamics of two nuclear factors having distinct functional roles: the protooncogene cMyc and the positive transcription elongation factor (PTEFb). cMyc is a basic helixloophelix DNAbinding transcription factor that binds to EBoxes; 18,000 Eboxes are found in the genome, and cMyc affects the transcription of numerous genes (Gallant and Steiger, 2009). Recently, cMyc has been demonstrated to be a general transcriptional activator upregulating transcription of nearly all genes (Lin et al., 2012; Nie et al., 2012). PTEFb is an essential actor in the transcription regulation driven by RNA Polymerase II. PTEFb is a cyclindependent kinase, comprising a CDK9 and a Cyclin T subunit. It phosphorylates the elongation control factors SPT5 and NELF to allow productive elongation of class II gene transcription (Wada et al., 1998). The carboxyterminal domain (CTD) of the catalytic subunit RPB1 of polymerase II is also a major target of PTEFb (Zhou et al., 2012). cMyc and PTEFb are therefore two good examples of transcriptional regulators binding to numerous sites in the nucleus; the latter binds to the transcription machinery itself and the former directly to DNA.
Single particle tracking (SPT) constitutes a powerful method to probe the mobility of molecules in living cells (Lord et al., 2010). In the nucleus, SPT has been first employed to investigate the dynamics of mRNAs (Fusco et al., 2003; ShavTal et al., 2004) or for rheological measurements of the nucleoplasm using inert probes (Bancaud et al., 2009). Recently, the tracking of single nuclear factors has been facilitated by the advent of efficient in situ tagging methods such as Halo tags (Mazza et al., 2012). An alternative approach takes advantage of photoconvertible tags (LippincottSchwartz and Patterson, 2009) and photoactivated localization microscopy (PALM) (Betzig et al., 2006; Hess et al., 2006). Single particle tracking PALM (sptPALM) was first used to achieve highdensity diffusion maps of membrane proteins (Manley et al., 2008). However, sptPALM experiments have typically been limited to proteins with slow mobility (Manley et al., 2008) or those that undergo restricted motions (Frost et al., 2010; English et al., 2011). Recently, by inclusion of lightsheet illumination, it has been used to determine the binding characteristics of TFs to DNA (Gebhardt et al., 2013).
In this study, we developed a new sptPALM procedure adapted for the recording of individual proteins rapidly diffusing in the nucleus of mammalian cells. We used the photoconvertible fluorophore Dendra2 (Gurskaya et al., 2006) and took advantage of tilted illumination (Tokunaga et al., 2008). A careful control of the photoconversion rate minimized the background signal due to outoffocus activated molecules, and we could thus follow the motion of individual proteins freely diffusing within the nuclear volume. With this sptPALM technique, we recorded large data sets (on the order of 10^{4} single translocations in a single imaging session), which were essential for a proper statistical analysis of the search dynamics.
We applied our technique to several nuclear proteins and found that diffusing factors do not sense a unique nucleoplasmic architecture: cMyc and PTEFb adopt different nuclear spaceexploration strategies, which drastically change the way they reach their specific targets. The differences observed between the two factors were not due to their diffusive kinetic parameters but to the geometry of their exploration path. cMyc and our control protein, ‘free’ Dendra2, showed free diffusion in a threedimensional nuclear space. In contrast, PTEFb explored the nuclear volume by sampling a space of reduced dimensionality, displaying characteristics of exploration constrained in fractal structures. The role of the spacesampling mode in the search strategy has long been discussed from a theoretical point of view (de Gennes, 1982a; Kopelman, 1986; Lindenberg et al., 1991). Our experimental results support the notion that it could indeed be a key parameter for diffusionlimited chemical reactions in the closed environment of the nucleus (Bénichou et al., 2010). We discuss the implications of our observations in terms of gene expression control, and its relation to the spatial organization of genes within the nucleus.
Results
Intracellular singlemolecule tracking with photoconvertible fluorescent proteins
We developed a simple and versatile approach based on photoconvertible protein tags that extends the use of sptPALM to any protein expressed in mammalian cells. Proteins of interest were fused to the photoconvertible protein Dendra2 (Gurskaya et al., 2006; Figure 1A). A standard widefield configuration of the microscope allowed fast and sensitive acquisition with an EMCCD camera (‘Materials and methods—Singlemolecule imaging and Detection and tracking of single molecules’, Figure 1—figure supplement 1). We used low activation intensity and tilted illumination (Figure 1B) in order to reach the regime of SM detection, characterized by singlestep activation and photobleaching (Figure 1C). Due to activation of outoffocus fluorophores, a decreasing density of detected particles was correlated with an increasing average signaltonoise ratio (SNR) (Figure 1D). We found that activation intensity around 0.01 kW/cm^{2} offered the best tradeoff between the number of detected particles (∼1) and SNR.
Compared to membrane proteins or other proteins with constrained mobility, diffusion dynamics of intracellular molecules is much higher and can exceed 10 μm^{2}/s. Images recorded for such fast moving objects depart from the welldefined point spread function (PSF) of the microscope and exhibit a motion blur that cannot be characterized with standard Gaussian localization algorithms (Thompson et al., 2002). Therefore, we developed new localization and tracking algorithms (‘Materials and methods—Detection and tracking of single molecules’ and Figure 1—figure supplements 1 and 2) and validated them with simulations (‘Materials and methods—Numerical simulations’ and Figure 1—figure supplement 3). We could thus obtain single trajectories formed by individual translocations recorded every 10 ms. 50% of the traces were reconstructed with more than four time points, and some of them were as long as 60 consecutive translocations. The step size of single translocations ranged between tens of nanometers (limited by our localization accuracy of ∼70 nm) and ∼2 μm (Figure 1E and Videos 1–5). Hence, it became possible to track molecules with diffusion coefficients exceeding 10 μm^{2}/s.
System validation using ‘free’ Dendra2 and histone H2B fused to Dendra2
We first investigated two limit cases relevant to protein dynamics in the nucleoplasm: Dendra2 and DNAassociated histone H2B. Dendra2 is the fluorescent label that we fused to all other proteins used in our analysis. Green fluorescent protein (GFP) has no detectable interacting partners in mammalian cells (TrinkleMulcahy et al., 2008), and we therefore considered ‘free’ Dendra2 as a model for freely diffusing particles due to its structural similarity with GFP. In contrast, Dendra2 fused to histone H2B (Dendra2H2B) was expected to insert into chromatin and thus to display restricted motion.
Indeed, from a visual inspection, ‘free’ Dendra2 and Dendra2H2B trajectories (Figure 2A,B, respectively) exhibited obvious differences. Notably, translocation histograms for ‘free’ Dendra2 and for Dendra2H2B were not consistent with a single diffusing species (Figure 2—figure supplement 1, ‘Materials and methods–Cumulative histogram analysis and mean square displacement’), thus suggesting that displacements of these molecules were more complex than anticipated. Three distinct populations were needed to fit the translocation histograms at all time intervals (Figure 2—figure supplement 1).
To complement our analysis of the translocation histograms, we plotted the mean square displacement (MSD) of the molecules as a function of time (‘Materials and methods—Cumulative histogram analysis and mean square displacement’). For DendraH2B, the MSD reached a plateau after ∼20 ms at ∼ 0.5 μm^{2} (Figure 2C), consistent with a confined motion of individual histone molecules inserted into chromatin. The MSD of ‘free’ Dendra2 increased regularly with time. However, it slightly deviated from the linear behavior expected for molecules undergoing normal diffusion. This was attributed to a ‘population exclusion effect’ due to the different defocusing rates of the various diffusive subpopulations of Dendra2.
Because of their threedimensional motion in the nucleus, slow moving particles remained within the focal depth of observation (∼0.5–1 μm) for a longer time than fast moving ones. As a result, fast diffusing molecules contributed comparatively less than the slow ones to the MSD at longer time lags. Note that this effect is inevitable for any singlemolecule experiment involving more than one diffusive population and in which the threedimensional movement of particles is recorded in two dimensions (‘Materials and methods—Numerical simulations’ and Figure 2—figure supplement 2). We therefore adjusted the rates of the different diffusive populations for each molecule, and have used the corrected values through the text and for our analysis. The deviation from linearity of the MSD curve produced by such an exclusion effect clearly illustrates the need to complement the analysis of molecular mobility with other observables, ideally independent of population heterogeneity.
Finally, to carefully establish the range of application of our experimental and analytical methods, we performed numerical simulations (‘Materials and methods—Numerical simulations’). On the one hand, the particle localization precision sets the lower bound to a reliable estimation of the diffusion parameters, that is ∼ 0.04 μm^{2}/s for a pointing accuracy of ∼70 nm. On the other hand, fast moving particles can be tracked with a mobility up to ∼20 μm^{2}/s, beyond the experimental values determined for ‘free’ Dendra2. Altogether, our experimental and numerical results provide a benchmark for studying nuclear factors with a mobility ranging between that of chromatinbound H2B molecules and of ‘free’ proteins such as Dendra2.
cMyc and PTEFb differ in the nature of their diffusion
We next probed the mobility of transcription factors. Dendra2 was fused to the protooncogene cMyc and to the Cyclin T1 subunit of PTEFb. It has recently been shown that, rather than activating new sets of genes in the cell, the role of cMyc is that of an amplifier of transcription of already active genes (Lin et al., 2012; Nie et al., 2012). We thus tested the functionality of cMycDendra2 by performing RTqPCR on a set of active genes in our U2OS cell line. When comparing the wildtype cells and those expressing cMycDendra2, we measured an increase of RNA expression levels in 10 out of 12 tested genes (‘Materials and methods—mRNA expression and cMyc expression amplification analysis’).
Translocation histograms for cMyc were well fit with three diffusive populations (Figure 3—figure supplement 1). The most abundant corresponded to rapidly diffusing particles (13.5 μm^{2}/s, 70% of the molecules) (Figure 3A, black trajectories). In addition, a significant fraction of cMyc was immobile (9.5%) (Figure 3A, green trajectory) or displayed slow diffusion (D_{2} = 0.5 μm^{2}/s, 20.5%) (Figure 3A, blue trajectories). For PTEFb, the typical translocation length and the translocation histograms were comparable to those obtained for cMyc (Figure 3—figure supplement 2).
When plotting the MSD as a function of time for cMyc and PTEFb, we observed a deviation from linearity for both factors (Figure 3C). Such deviation could be due to the ‘population exclusion effect’ described above (‘Materials and methods—Numerical simulations’, Figure 2—figure supplement 2), but, alternatively, it could also be the signature of an anomalous diffusion process. When a particle undergoes anomalous diffusion, the MSD vs time scales as a power law t^{α}, where α < 1 is characteristic of a subdiffusion process (Saxton, 2007). However, neither the ‘free’ Dendra2 nor the cMyc MSD data could be properly fit by such a law (Figure 3D). Similarly to ‘free’ Dendra2, cMyc molecules were distributed between populations of very distinct diffusion coefficients. In contrast, for PTEFb, the MSD variations were remarkably fit by a t^{α} power law with the anomalous coefficient α = 0.6 (Figure 3D). The subdiffusion of PTEFb was also apparent when we plotted the cumulative histograms of the square displacement for multiples of the time interval (Δt) between two frames and rescaled them by the factor t^{α}, with α determined from the fit in Figure 3C. All the rescaled histograms curves collapsed remarkably well for PTEFb but not for cMyc or ‘free’ Dendra2 (Figure 3—figure supplement 3). We therefore concluded that the characteristics of single PTEFb trajectories are consistent with an anomalous diffusive behavior whereas the deviation from linearity of the cMyc MSD curve reflects the heterogeneity of its diffusion dynamics.
Asymmetric distribution of angles between consecutive translocations
Subdiffusion in cells is commonly attributed to one of the following two microscopic processes: a broad distribution of trapping times or an obstructed movement resulting from a reduction of the accessible space (Condamin et al., 2008) (for a discussion about subdiffusion causes, see ‘Materials and methods—Numerical simulations of anomalous diffusion models’). In other words, the subdiffusive behavior, evidenced by the sublinear MSD, is due to either temporal or spatial restrictions. In order to probe the spatial characteristics of the exploration independently of temporal considerations, we analyzed the distribution of angles Θ between two consecutive translocations, an observable that is predominantly sensitive to the geometry of the exploration space (Liao et al., 2012) and able to elucidate complex dynamics of molecules (Burov et al., 2013).
For ‘free’ Dendra2 and cMyc, we found a quasiuniform angular distribution (Figure 4A), as expected for Brownian diffusion. In a threedimensional space, there is no privileged direction and all angles Θ are equiprobable. In contrast, the angular distribution for PTEFb was significantly biased toward 180°, reflecting an anticorrelation between two successive displacements. Such anisotropic angular distribution is consistent with diffusion in a space of reduced dimensionality such as a fractal network (benAvraham and Havlin, 2005). A particle that diffuses in such a structure encounters dead ends, in which case it cannot but return back to previously visited locations (Θ = 180°). Noteworthy, the diffusing subpopulation of H2B molecules also showed a nonuniform angular distribution (Figure 4—figure supplement 1).
Temporal evolution and spatial dependence of the angular distribution
An alternative scenario in which such an asymmetric angular distribution of SM traces may arise is that of confined diffusion. If a diffusing particle is confined in a volume of size comparable to the translocation steps, its repetitive bouncing against the trap walls will produce a relative increase of angles larger than 90°. In such a case, the temporal evolution of the antipersistence reflects the length ratio between the displacement steps and the size of the confining volume. On the other hand, a defining property of fractal structures is their scale invariance (benAvraham and Havlin, 2005), namely the repetition of structural motifs at different length scales. For a particle diffusing in such a fractal structure, we expected the scale invariance to be apparent in the characteristics of the movement of the particle.
We therefore examined the temporal and spatial dependences of the angular distribution in order to further investigate the origin of the antipersistence of the trajectories and the underlying geometry of the space available for exploration. We defined the asymmetry coefficient (AC) as the logarithm to the base 2 of the ratio between the frequency of forward angles (between 0° and 30°) and backward angles (150°–180°) (Figure 4B). The AC is thus negative for angular distributions with a dominant number of backward angles, and it measures the deviation from a homogenous distribution. We calculated the AC for the angles formed at increasing lag times (Figure 4C, Figure 4—figure supplement 1) as well as a function of the average length of the consecutive translocations forming the angle θ (Figure 4D, Figure 4—figure supplement 1). It is important to note that with this analysis, the experimental localization accuracy is reflected in the first data point of the spatial dependence of the AC, and not for the data above the 0–150 nm bin. Also, fewer particles are contributing to the AC at larger times, as can be observed in the angular distribution histograms in Figure 4—figure supplement 1.
We found out that the angular distribution of cMyc deviates from homogeneity at increasing lag times with increasing negative AC (Figure 4C), potentially reflecting a hindrance to the free diffusion of cMyc and its confinement to domains significantly smaller than the nucleus. However, the angular distribution became isotropic (AC = 0) at translocations larger than 300 nm (Figure 4D). This transition can be interpreted as an indication of the upper limit size of the confining volume. Hence, the temporal and spatial evolution of the AC suggest two subpopulations of cMyc, one confined into regions smaller than ∼300 nm and a nonconfined fraction of cMyc molecules. PTEFb, on the other hand, displayed a remarkable constant value of AC for both, time and space (Figure 4C,D), possibly reflecting a lengthinvariant property of the medium in which diffusion takes place.
Numerical simulations: intermittent diffusion and particle diffusion in media of increasing complexity
In order to gain insight about the different scenarios giving rise to the observed angular distributions, we performed numerical simulations of models with increased levels of complexity (see ‘Materials and methods—Numerical simulations’ for details about the numerical simulations). In line with the observation of different diffusing populations even for free Dendra2, we first considered an intermittent diffusion model. Here, particles had a probability to switch from a fast to a slow diffusion coefficient and vice versa. We also considered an intermittent trap model, where diffusing particles with fixed diffusion coefficient have a probability to be confined in a spherical trap. We adjusted the parameters of the models in order to obtain similar translocation histograms to those of cMyc and PTEFb (Figure 5—figure supplement 1, ‘Materials and methods—Numerical simulations’). However, none of these simple intermittent models reproduced the antipersistent characteristics of the experimentally measured trajectories (Figure 5—figure supplement 2, ‘Materials and methods—Numerical simulations’).
We then considered a model that results from a combination of intermittent diffusion and intermittent trap. We performed simulations of fast diffusing particles (diffusion coefficient D_{1}) with a probability Kon to engage into a slower diffusion (D_{2}) confined in a trap of radius R (Figure 5A). Here, the AC decreased with increasing lag times (Figure 5B), reproducing the trend observed in cMyc. Likewise, the AC displayed the same behavior as cMyc, tending to zero for larger values of the translocation steps (Figure 5C). Following this model, cMyc performs thus a free exploration of the nuclear space, combined with slower yet still normal diffusion of confined domains, reflecting its interactions with a multiplicity of partners.
Finally, in order to reproduce the invariant properties of the angular asymmetry observed for PTEFb, we needed to invoke a hierarchical organization of the space. We considered the intermittent trap model, this time with a distribution of trap sizes governed by a Pareto power law (exponent 0.1). With this model, we obtained an antipersistent angular distribution (Figure 5D) and a closer reproduction of the AC behavior observed for PTEFb. Although the AC was not strictly constant with time, it did not show a tendency towards zero (Figure 5E). Moreover, the spatial dependence of the angular asymmetry formed a plateau for translocations larger than 600 nm (Figure 5F).
Such a hierarchy of confining sizes led us to consider a fractal network as an underlying structure on which to simulate the diffusion of particles, also motivated by recent works on the geometry of the nuclear space (Bancaud et al., 2012). We considered a 3D percolation cluster as well as a 2D Sierpinski carpet. The Sierpinski carpet is an exact fractal lattice with multiscale selfsimilarities. Random walks on a Sierpinski lattice are anomalous because its structure induces spatial correlations between successive displacements. The percolation cluster at the critical percolation threshold possesses the property of statistical internal selfsimilarity. As a consequence, the percolation cluster exhibits fractal properties without a defined geometric shape (benAvraham and Havlin, 2005). For both fractal structures, the angular anisotropy was constant with time (Figure 5E), illustrating the scaleinvariant features of fractal structures, as observed in the experimental data of PTEFb. Surprisingly, the AC decreased for larger translocations in the case of the percolation cluster, while the Sierpinski carpet yielded an invariant asymmetry in space. This was interesting because it indicates that the underlying network needs to reproduce a certain degree of geometrical selfsimilarity, as it is the case of the Sierpinski carpet. The percolation cluster, on the other hand, does not conserve its geometry at different scales but rather other features like the local density obey a power law.
Targetsearch and sampling: compact vs noncompact space exploration cMyc and PTEFb adopt opposed search strategies
We have determined that while cMyc undergoes normal diffusion (with a subpopulation seemingly confined in domains smaller than the nucleus), the dynamics of PTEFb is well described by a subdiffusive behavior. In the case of PTEFb, our simulations support the notion that anomalous diffusion is compatible with an obstructed mobility of the proteins, as obtained on a fractal structure (we have ruled out other models of subdiffusion, see Figure 5—figure supplement 3 and ‘Materials and methods—Numerical simulations of anomalous diffusion models’ for a more detailed discussion). As previously described, the exponent α = 0.6 of anomalous diffusion obtained for PTEFb (Figure 3D) is a direct measure of the dimension of the walk D_{w} = 2/α = 3.3. Since the fractal dimension D_{f} has an upper limit at D_{f} = 3, we can therefore conclude that D_{w} > D_{f}, and thus that PTEFb is engaged in a compact exploration of the nucleoplasm. In contrast, the isotropic sampling of space of cMyc excludes a compact mode of exploration; it undergoes normal 3D diffusion irrespective of its confinement, and hence the dimension of the walk is D_{w} = 2, sampling the nucleoplasm in a noncompact manner. These results imply that different factors sense a proteindependent nuclear environment, which can be determinant for their exploration strategy.
The distancedependence of the mean first passage time differs between cMyc and PTEFb
The distinctive properties of compact and noncompact trajectories have potentially important functional consequences on the ability of searchers to find and react with molecular partners. As noted above, a striking difference is the distancedependence of the mean first passage time (MFPT) of the searcher to the target site. The MFPT of noncompact explorers is essentially constant, depending solely on the total volume and not on the distance r to the target. Conversely, in the compact case, the MFPT still scales with the volume but also increases with the distance as ${r}^{\left({D}_{w}{D}_{f}\right)}$.
As an illustration, we computed the MFPT as a function of the distance (see analytical expressions of MFPT in Condamin et al., 2005; Bénichou et al., 2010), using the experimental data for cMyc and PTEFb, two examples of noncompact and compact explorers. For cMyc, which behaves as an ordinary Brownian walker, the fractal dimension is D_{f} = 3, and the dimension of the walk is D_{w} = 2. We used a diffusion coefficient D = 9.8 μm^{2}/s, the value obtained by a weighted average of the diffusion coefficients of the three subpopulations. (It is important to note that the value used for the diffusion coefficient does not affect the dependence of the MFPT on the initial distance to the target.) To calculate the MFPT, we used a nuclear volume of 600 μm^{3} and considered a target in its center. For PTEFb, we did not have direct access to the value of D_{f} and used several values previously reported as estimations in the nucleoplasm (Bancaud et al., 2012). In Figure 6, we used D_{f} = 2.6 and the results were qualitatively similar for values of D_{f} = 2.2, and D_{f} = 3 (Figure 6—figure supplement 1). For both proteins, we also varied the size a of the target between 1 nm (i.e., corresponding to a couple of base pairs), 10 nm (the size of a protein complex), and 100 nm (the size of a large multimolecular complex).
For cMyc, the MFPT was constant, irrespective of the distance r (Figure 6A). However, it was inversely proportional to the size of the target, similar to what is predicted from the diffusionlimited rate of bimolecular reactions (Nelson et al., 2008). In contrast, the MFPT of PTEFb increased with the distance r but did not depend on the target size. The lack of size dependence can be simply viewed as a consequence of the redundant exploration of compact explorers, and reflects the fact that the limiting step to find a target is the time taken to reach its vicinity. We stress that the differences of MFPT can be very significant. For instance, the time needed to find a 10 nm target located at a distance of 250 nm is 68 times longer for cMyc compared to PTEFb (506.1 s for cMyc and 7.4 s for PTEFb). If the target is located at 5 μm of the TF, the difference in the search time is reduced to a factor of 8 (525.3 s for cMyc and 64.6 for PTEFb).
Here, we considered that cMyc has a full access to the nuclear volume. It is interesting to note that if, as suggested by the temporal variance of the angular distribution, cMyc is confined to a smaller domain, the MFPT would scale linearly with this volume.
We also considered the case of a factor susceptible to bind to two different targets T1 and T2 (Figure 6B). To do so, we computed the splitting probability P, that is the probability to reach T1 before T2 as a function of the initial distance to T1. For cMyc, the probability was equal to 0.5 as soon as the initial distance was larger than a few tens of nanometers, in stark contrast with the case of PTEFb, for which P varied almost linearly with the distance.
Overall, our analysis of SM experiments of cMyc and PTEFb reveals two characteristics of TFs diffusion relevant to the understanding of transcription regulation kinetics. First, the exploration geometry of the nucleus by TFs is determined by the function and interactions of the nuclear factor. Rather than being subjected to a universal sampling geometry imposed by the nuclear architecture, cMyc and PTEFb adopt different modes of exploration leading to normal and anomalous diffusion, respectively. Second, despite apparently similar diffusion coefficients, the different exploration strategies of cMyc and PTEFb (noncompact and compact, respectively) can lead to opposite dependence of the search kinetics on the distance to the target and on the target size. The distancedependence of the MFPT has direct implications on the probability of interaction of cMyc and the PTEFb with their respective partners, which in turn may affect transcriptional kinetics and regulation.
Discussion
Proteinspecific sensing of nuclear organization
With the PALM imaging assay adapted for SM detection of intracellular proteins in eukaryotic cells, we probed the spatial dynamics of different proteins in the nucleus of live human cells: ‘free’ Dendra2, histone H2B, the protooncogene cMyc, and the elongation factor PTEFb. The analysis of individual trajectories, supported by numerical simulations of diffusive tracers on free, confined, and fractal structures, and switching between different regimes, shows that these nuclear proteins fundamentally differ in their exploration of the nucleoplasm. Our results on ‘free’ Dendra2 are along the lines of those obtained with microinjected fluorescent streptavidin, which explores all nuclear compartments with three subpopulations having different diffusion characteristics (0.15, 0.8, and 5 μm^{2}/s) (Grünwald et al., 2008), possibly reflecting differences in viscosity and/or crowding in the nucleus. In contrast, FCS experiments using ‘free’ GFPrepeats or SPT tracking of QD aggregates suggested anomalous diffusion (Bancaud et al., 2009).
We determined that ‘free’ Dendra2 and the protooncogene cMyc undergo normal diffusion in 3D, whereas the displacement of PTEFb was accounted for by a subdiffusive movement. This finding was further supported by measurements of the distribution of angles between consecutive translocations. Importantly, this distribution was initially isotropic for Dendra2 and cMyc and an asymmetry towards the return angles increased over time, as expected for confined Brownian motion. Conversely, PTEFb showed a pronounced and timeinvariant anisotropy consistent with the motion on a fractal structure. Thus, the nuclear geometry, or equivalently, the architecture of the space sampled by diffusing factors, is not unique but constitutes a proteinspecific parameter. Furthermore, taking into consideration the diffusion parameters derived from the analysis of the MSD, together with the geometrical aspects of the exploration of cMyc and PTEFb, we determined the mode of exploration of these factors to be noncompact and compact, respectively.
We stress that the distinction between compact and noncompact exploration, rather than the one between anomalous and normal diffusion, is the proper criterion to analyze the search dynamics of transcription factors. The notion of compactness is intimately linked to the geometry and the dimensionality of the sampled space. In this regard, there is a specificity of random motions in a threedimensional medium with respect to the one and bidimensional cases, for which the exploration is always compact since the fractal dimension D_{f} (less or equal to 1 and 2, respectively) is necessarily smaller than D_{w}. Only in the case of 3D search, can both compact and noncompact behaviors be observed. Our data demonstrate the relevance of the notion of compactness for the description of nuclear factor dynamics.
Possible mechanisms controlling the geometry of nuclear explorations
One microscopic mechanism leading to a compact exploration of the nucleus could be a compartmentalization of the nucleoplasm into interconnected domains forming a fractal labyrinth in which molecules diffuse. In our view, such a model assuming that molecules encounter physical barriers is poorly compatible with the dynamic nature of nuclear organization and with the lack of correlation between protein size and mobility in the nucleus (Sprague et al., 2004; Mueller et al., 2008).
Another interpretation is that of a fractal structure restricting the mobility of proteins at its surface. Chromatin has been described as a fractal globule (Grosberg et al., 2007; LiebermanAiden et al., 2009) and transient, nonspecific interactions to a continuum of binding sites would account for the diffusing factors not escaping from their interaction with chromatin. In this scenario, the number of binding sites with which cMyc interacts is not sufficient to restrict its motion to chromatin (36,000 Eboxes in a diploid genome, representing less than 50 sites per μm^{3}).
PTEFb interacts with the CTD of the catalytic subunit of RNA Polymerase II, which contains 52 repetitions of a heptapeptide motif (Taube et al., 2002). The RNA Polymerase CTD is not folded and can occupy the space very efficiently, potentially forming a mesh offering a nuclear continuum of binding sites for PTEFb. Such CTD matrix could have an intrinsic existence or be linked to the chromatin globular organization. The existence of a nuclear protein scaffold or matrix has been speculated for more than half a century (Pederson, 2000) and both our works offer an observation of a functional role for such a structure. Several other studies support this hypothesis, showing that nuclear proteins are in constant interaction with their environment and their motion is governed by specific and nonspecific bindings (Phair et al., 2004; Sprague et al., 2004; Hager et al., 2009; Speil et al., 2011), therefore opening the door for mechanisms where factors are guided on networks of binding sites (Bénichou et al., 2011).
The effect of the exploration strategy on gene regulation by transcription factors
From a general standpoint, the distancedependence of the search kinetics could have strong implications for gene regulation. For example, it has been recently shown that, in Escherichia coli, the spatial distribution of TFs is determined by the local state of DNA (Kuhlman and Cox, 2012). Let us consider the case of TFs coregulating multiple loci; the relative localization of these loci is an important parameter that will play different roles depending on the compact or noncompact exploration of the TFs. Noncompact TFs have a very similar probability to bind to all loci. In other words, all loci will have the same probability to be occupied, regardless of their spatial position. In contrast, compact factors will be preferentially shared between proximal loci, and therefore the probability of a locus to be occupied by a compact explorer is a function of the occupation history of its neighboring sites: it is distance and time dependent. Importantly, this indicates that two loci, such as two regulatory sites located a few tens of kbp away from each other, can transfer information and influence one another without direct physical contact. This spatial relation could underlie the process of sequestration of factors away from their targets (Yao et al., 2011), which would occur only with compact explorers. Such geometrically controlled longdistance interactions are not detectable using conventional chromatin capture assays, which predominantly rely on the chemical crosslinking between contacting sites.
Compact transcription factors and the stability of molecular complexes
A remarkable feature of compact searchers is their propensity to visit their neighboring sites multiple times. As a result, they have a probability equal to one to return to a site that they previously occupied, a property designated as the recurrence of compact trajectories. From a biochemical viewpoint, this property might affect our understanding of the kinetic stability of molecular complexes. Certainly, molecular machines controlling the nuclear functions such as transcription, splicing, and replication are composed of large numbers of molecules. Some of these molecules are stable constituents while others can be rapidly exchanged in order to control the specificity and modulate the activity of a particular complex (Fong et al., 2012). It is therefore important to understand how these molecular machines can assemble from their principal components. For instance, we cannot yet reconcile the need for strong and stable interactions, believed to be required for the viability of such complexes, and the requisite of weak and transient interactions required for molecules to compete for the same target regulating their composition. The observation of compact modes suggests that strong binding, associated to small dissociation rates, is not required to ensure high occupancy.
Compact transcription factors favor transcriptional bursting
Recently, the role and importance of transcriptional fluctuations within a single cell have been extensively studied (Raj et al., 2008; Zenklusen et al., 2008; Larson et al., 2009; English et al., 2011; Itzkovitz and van Oudenaarden, 2011). Using a simple model in which the activation of a gene is controlled by the binding of a single TF to a locus, Meyer et al. (2012) have modeled how the search dynamics of these TFs affects the transcriptional response. In this model, for compact TFs, the recurrence of the trajectories and the facilitated reassociation to the locus would result in transcriptional bursting. In contrast, for the noncompact case, the gene activation rate is determined by the total TF concentration in the nucleus, and the transcriptional activity is uncorrelated in time. This further illustrates how the translocation properties of nuclear factors might underlie the kinetics of functional cellular events.
Conclusion
In this study, we have experimentally demonstrated that different nuclear proteins with different functions sample the nucleoplasm with different search strategies: the exploration geometry of the nucleus is proteindependent. We have also determined that two different universality classes of search modes, namely compact and noncompact explorations, coexist in the nucleoplasm. Our current view of the nucleoplasm and chromatin is that of a structure whose condensation influences its accessibility to transacting factors. Here, we have established that the same target in the nucleoplasm can be visited with different probability and kinetics by different factors depending on how they sample space. In addition to chromatin condensation, the compactness of the exploration itself needs to be taken into account to understand how gene regulation operates.
While the spacesampling mode of a random exploration is either compact or noncompact, the question to be answered in the future is whether one molecule manifests both types of search, exhibiting transitions between them, and whether different dynamics may still arise within each search mode. If that is the case, it will be of paramount importance to understand the level of regulation of such transitions, as well as the implications for the kinetics of the transcription process. The inverse first passage time is a measure of the reaction rate constant. Therefore, the different interaction kinetics that results from compact or noncompact explorations has profound implications in the understanding of the dynamic interactions and reactivity rates between TFs and corresponding regulated genes. For a noncompact explorer like cMyc, the interaction rate is that of a homogenous solution, and thus it will bind with equal probability to any target in the nucleoplasmic volume. On the other hand, for a compact exploration such as the one of PTEFb, the recurrent search of the local environment and the distance dependence of the search time translate into spatial and temporal correlations between binding events. Such spatial correlation can be seen as a mechanism that adds a level of control to the rapid assembly of molecular complexes, reconciling weak and transient interactions with functional stability. This last notion suggests the idea of a regulated level of compactness of TFs both in time and space.
Materials and methods
Cell culture and transfection
Request a detailed protocolU2OS (Human Osteosarcoma) cells were grown in DMEM (Life Technologies, Carlsbad, CA) with 1 g/l glucose and glutamax supplemented with 10% FBS (Fetal Bovine Serum, Life Technologies) and 1% Penicillin/Streptomycin (Life Technologies) at 37°C with 5% CO_{2}. 48 hr prior to the imaging, cells were seeded at 30–40% confluence on a plasmacleaned (2 min with air with Femto model, Diener Electronic, Ebhausen, Germany) and collagencoated (Collagen I from Rat tail, Life Technologies) coverslips (N°1 25 mm, Marienfeld, LaudaKönigshofen, Germany).
The C terminal of cMyc and H2B were fused to Dendra2 and expressed under the control of the CMV promoter. Prior to experiments, U2OS cells were transfected 24 hr before imaging with the plasmid of interest (100 ng/25 mm coverslip) using Fugene 6 (Roche Applied Science, Penzberg, Germany) according to manufacturer's instructions. Clones with very low overexpression of exogenous protein, as judged by low fluorescence intensity of preconverted Dendra2, were used.
Experiments with PTEFb (Cyclin T1 fused to Dendra2 on N terminal) and Dendra2 (alone) were performed on U2OS cell line stably transfected and selected with geneticin (Life Technologies). Clones with very low expression of fluorescent protein (CyclinT1Dendra2 or Dendra2), as judged by low fluorescence intensity of preconverted Dendra2, were used. Transient transfections of Cyclin T1 Dendra2 were also performed and gave the same results.
Singlemolecule imaging
Request a detailed protocolSinglemolecule imaging was performed on an inverted microscope Nikon Ti Eclipse (Nikon Instruments, Tokyo, Japan), with a high numerical aperture objective (1.49 NA) and 100X magnification; extra magnification of 1.5X was used in the tube lens of the microscope, resulting in a total magnification of 150X. We also used perfect focus system (Nikon) designed to avoid drift on the Zaxis (focus) of the objective, relative to the coverslip. The excitation (561 nm) and activation (405 nm) laser beams were injected into a fiber and focused in the back focal plane of the objective, using an appropriate dichroic (Di01R56125x36) (Figure 1—figure supplement 1A). A motorized mirror allowed us to choose between widefield or inclined excitation configurations; a small angle, between 0 and 30°, was typically used to avoid straylight reflections and reduce background from cell autofluorescence. Experiments were acquired under continuous excitation (561 nm laser, 5 kW/cm^{2} on the sample) and pulsed activation (405 nm laser, 1 pulse of 10 ms per second, 0.01 kW/cm^{2} during the pulse on the sample). Fluorescence emission from individual Dendra2 molecules was filtered with a single band emission filter centered at 617 nm and a bandpass of 73 nm and recorded on an EMCCD camera (iXon 897 Andor Technology, Belfast, Ireland). The pixel size of the EMCCD was 16 μm, and we imaged a small region of interest (ROI) of about 100 pixels × 100 pixels. This ROI was sufficient for imaging a large crosssection within the nucleus of single cells, and allowed acquisition rates as fast as 100 Hz (10 ms per frame). Images of the preconverted (green) form of the ensemble fluorescence of Dendra2 were taken using a mercury lamp for illumination (excitation: 485 nm, emission FF01525/30).
Cells were imaged in Leibovitz's L15 medium (Life Technologies) containing 10% FBS (Fetal Bovine Serum, Life Technologies). The sample was placed on the microscope, on a stage heated at 37°C on the microscope. Once an ROI was selected from the preconverted (Dendra2 greenform) fluorescence imaging of the live cells, activation pulses were fired every 100 frames, and videos of several thousands of frames were acquired under continuous 561 nm illumination (typically 2000 to 10,000 frames per cell). Each coverslip was used for a maximum of 45 min after placing them on the scope.
Preconverted Dendra2 imaging
Request a detailed protocolThe same conditions that were used for SM imaging were used to obtain the images of the preconverted ensemble fluorescence of Dendra2, but exchanging the light source for a mercury Lamp (Intensilight, Nikon) and appropriate excitation and emission filters (485/20 nm and 525/30 nm, respectively). In order to compensate the very weak expression levels, images were reconstituted averaging 100 images of a temporal sequence therefore minimizing the noise.
mRNA expression and cmyc expression amplification analysis
Request a detailed protocolBased on RNA Pol II chIPSEQ data available in the laboratory, we selected genes that are expressed in U2OS cells. Those genes were: SPG21, LMF1, BEX1, IGF2R, GAPDH, HMGB2, SOD1, RPL30, ORC3, CUL1, TRAF5, STX11. Using RTqPCR, we compared the mRNA expression of these genes in two conditions: wild type untransfected U2OS and cMycDendra2 transfected U2OS. In order to precisely compare the amount of RNA, we counted and fluorescenceactivated cell sorted (FACS) the same number of untransfected and cMYCDendra2expressing cells. We then performed quantitative PCR experiments and compared the expression levels of the analyzed RNA in the two different conditions.
RNA was purified using TRIzol Reagent (15596018; Invitrogen, Life Technologies) according to the manufacturer's instructions. Total RNA was quantified on a NanoDrop 2000c Spectrophotometer (Thermo Fisher Scientific, Waltham, MA) and their quality was evaluated on RNA Nano Chips (50671511; Agilent 2100 bioanalyzer, Agilent, Santa Clara, CA). Reverse transcription from total RNA to cDNA was done with oligodT (18418020; Invitrogen, Life technologies) using SuperScript III RT (18080085; Invitrogen, Life technologies) and RNAse OUT (10777019; Invitrogen, Life technologies).
Quantitative realtime PCR (qPCR) was done using 5 μl of 1:20 diluted cDNA on a LightCycler480 system (Roche, Basel, Switzerland) using Maxima SYBR Green qPCR Master Mix (K0252; Fermentas, Thermo Fisher Scientific). A final concentration of 500 nM of primer pairs (Eurofins, MWG, Huntsville, Al, designed according to DugastDarzacq and Grange, 2009) was used for each qPCR reaction. The cycling conditions were as follows: 95°C for 10 min, 45 cycles (95°C, 15 s; 58°C, 30 s; 72°C, 20 s) and melting curve analysis. LightCycler 480 SW 1.5 was used to evaluate and to analyze the data.
Analytical methods
Detection and tracking of single molecules
Request a detailed protocolDetection
Request a detailed protocolThe diffusion of the molecules imaged in the nucleus of eukaryotes can be as fast as ∼10 μm^{2}/s. This implies that the detected molecules can travel a distance larger than the diffraction limit of light (∼250 nm) during the characteristic acquisition time (10 ms). The intensity profile is therefore a convolution between the Airy pattern of the point spread function (PSF) and the trajectory of the particle during the 10 ms acquisition time (Figure 1—figure supplement 1B). While such motion blur contains potentially useful information (Elf et al., 2007), it has some detrimental consequences: a decrease of the SNR and the ineffectiveness of traditional Gaussian fit localization algorithms. An approach, demonstrated in bacteria, to minimize the blurring effect consists in illuminating the sample with brief (1 ms or less) and intense (up to 100 kW/cm^{2}) laser pulses (Elf et al., 2007). In bacteria, this stroboscopic method is all the more necessary since the extension of the motion blur is often comparable to the size of the cell itself (section of 1 μm^{2}). However, given the larger dimension (section ∼200 μm^{2}) of a mammalian nucleus, this method requires high laser power, not practical with standard microscopes or live cell microscopy due to phototoxicity effects. Hence, we favored an approach using lower intensity (∼4 kW/cm^{2}) and longer illumination time (∼10 ms), limited by the readout rate of our camera. At this time scale, trajectories are not affected by the nuclear confinement.
In this case, the emission of fast diffusing single fluorophores cannot be detected with traditional twodimensional Gaussian fit algorithms (Cheezum et al., 2001; Abraham et al., 2009). We developed an alternative, comprehensive algorithm capable of detecting fast diffusing molecules that typically have low signaltonoise ratio (SNR) as well as immobile particles with higher SNR.
For each frame, background intensity was estimated at each pixel as the median intensity of the pixel over the entire video. This background was subsequently subtracted from the raw image. Fluorescence signal from individual molecules may still appear as an aggregate of disconnected pixels, therefore a smoothing step was applied using a Gaussian mask with standard deviation of σ = 121 nm. Those pixels with an intensity corresponding to the highest 20% of the nonsmoothed (but background corrected) image were selected (Figure 1—figure supplement 1C). At such threshold level random noise fluctuations were still included in the pixel selection, we therefore disregarded any spot that spanned less than 0.2 µm² (∼20 pixels, or about half the theoretical optical response of the system). Individual pixel aggregates were then selected for each frame, with one additional constraint to account for molecules diffusing outside and back inside the focal plane during the acquisition time; we considered detected spots closer than 500 nm as originating from the same molecule. The position of each spot was calculated as the center of mass of the pixel aggregate, which is a good estimator of the particle position suggested by deconvolution approaches (Michalet, 2010).
We tested our detection algorithm with simulated videos consisting of white noise (without single particles signals) with pixel intensity values and standard deviation comparable to the background noise of our live cell data, resulting on a detection rate of 10^{−6} detections per frame per μm², three orders of magnitude lower than the typical detection rates obtained with the experimental data.
The localization accuracy of the detection algorithm could also be estimated. We calculated the standard deviation of the position coordinates of a H2B molecule, detected in 290 consecutive frames (3 s tracking) with no apparent diffusion. As shown in Figure 1—figure supplement 1D, we obtained a localization accuracy of ∼70 nm.
Tracking
Request a detailed protocolIn order to connect consecutive detections of one given molecule, we defined the maximum distance R allowed for the translocation of a single step of the particle. For each single particle detection, the radius R defined an area centered on the particle position at time T on the consecutive frame at T + Δt. When a detection at T + Δt was found within the area defined by R, the two detections were linked in a trajectory. When two or more particles were detected within this area, the trajectory was truncated and the positions considered as the first detection of new trajectories. When the detection of one particle could be included in two different trajectories, both trajectories were also truncated, and the detection was disregarded. We defined such a restrictive policy of tracking in order to reduce the number of misconnections, or falsepositive tracking connections. A misconnection occurs when two consecutive detections from two different molecules are included in the same tracking sequence. Therefore, when there is any ambiguity between two spatially closed detections, the algorithm truncates the trajectories.
Such restricting tracking policy reduced the number of misconnection but also reduced the total number of traces suitable for analysis. Therefore, in order to set an appropriate maximum radius R, we computed the probability of detecting two different molecules in consecutive frames within a distance lower than R.
Maximum tracking radius R and misconnection probability
Request a detailed protocolConsidering the detection of a given molecule in consecutive frames, we could estimate the probability of tracking error by determining the local density of molecules different than the molecule of interest. In order to do so, we determined the local particle density at a time point where the probability of detecting the same particle is close to zero. We estimated the fluorescence photobleaching characteristic time under our experimental conditions by measuring the fluorescence lifetime of an ensemble of proteins in the nucleus, after a high intensity activation pulse, and under usual imaging conditions (Figure 1—figure supplement 2A). We measured a fluorescence halflife of ∼600 ms, suggesting that the probability of a molecule photobleaching between two consecutive frames is 0.02. After 5 s (476 frames) of the initial detection, it is highly improbable (0.0001 probability) that a detection originates from the same molecule. For each molecule detected, we could therefore calculate the number of detections around the same spatial coordinates but at a time separation of 5 s or more and thus estimate the average local density within a radius R of the molecule.
Considering the set of all the detections $\left({x}_{i},{y}_{i},{t}_{i}\right)$, where (x_{i}, y_{i}) are the spatial coordinates and t_{i} the time, for each detection i, we defined N_{i}, the total number of frames recorded more than 5 s after each given detection i. The detections made during this period are estimators of the local density around detection i. We therefore defined W_{i}(R) as the total number of detections during these timeshifted frames within a distance smaller than R (i.e., the total number of detections within R, after 5 s). Being M the total number of detections, we could calculate the total number of expected misconnections within a radius R as follows:
By comparing the value obtained from this expression to the number of connections we measured, we could estimate the probability to misconnect two detections. The total number of connections C(R) could be then calculated integrating both, the misconnections and the positive translocations. For every detection _{i} we calculated C_{i}(R), the number of detections at the consecutive frame at a distance smaller than R. The sum of all detections was therefore:
In Figure 1—figure supplement 2B, we plotted the measured W(R) and C(R) as well as their difference, for the free fluorophore Dendra2 as well as for all the proteins under study. We observed that for R bigger than ∼2 µm the total number of tracking assignments was dominated by misconnections; we therefore set the maximum allowed radius R for tracking under our imaging conditions to be 2 µm.
A tracking misconnection could occur at the first or last translocations of a trace, or in the middle of the trajectory. If a false connection occurred in the middle of the trace, its origin was the erroneous link of two traces from different molecules, being the first one detected until frame i and the second one starting at frame i + 1, appearing in the vicinity of the first molecule. The probability of such event to happen is very low because the number of single frame detections outnumbered by two orders of magnitude the number of trajectories with at least two consecutive detections. We could therefore consider that the tracking error misconnections occurred mainly at the beginning or end of the trajectories. Taking into account only the detections that are at the extremities of a trajectory, after tracking with a maximum radius of 2 µm, we could therefore estimate the probability of false connection from the fraction of misconnections as:
where W′ and C′ stand for the first order derivative of W and C, which were estimated for mathematical convenience every 0.05 µm and then linearly extrapolated.
To extend the notion of onestep translocation error probability to several steps displacement, the probability for a trajectory to be false was set to be the mean of all the onestep translocations that composed the trajectory.
Cumulative histogram analysis and mean square displacement
Request a detailed protocolThe data of proteins diffusing in the nuclear volume is the 2D projection of a 3D motion. Provided that the nucleus is isotropic along the three spatial axes X, Y, and Z, the XY projection data fully reflect the 3D behavior of the molecules.
The analysis of the cumulative translocations histogram allows for the determination of individual components from a mixed set of translocations, that is, translocation steps that cannot be governed by a single diffusion coefficient. For the 1Δt time step (10.5 ms), the cumulative distribution function (CDF) is a function F_{1}(x) that represents the probability that a random translocation may be found at a distance smaller than x.
The cumulative function weighted with the probability of misconnection is:
where the sum is computed for all the recorded translocations d_{i}, H represents the Heaviside step function, which is 1 for x − d_{i} ≥ 0 and 0 for x − d_{i} < 0 and P^{w} is the misconnection probability described in the previous section.
The probability distribution of 2D translocations for a single diffusion coefficient D and an interframe lag time T is given by $\frac{1}{\sqrt{4DT}}{e}^{\frac{{x}^{2}}{4DT}}$. This imposes a CDF of translocations for a set of translocations that can be described by single population of diffusion to be:
When the measured translocations reflects a pool of molecules with different diffusion kinetics, this single exponential function fails to describe the empirical CDF. In the case of k diffusing species, the empirical CDF is best described by
where a_{i} represents the fraction of translocations in the probability distribution imposed by a diffusion coefficient D_{i}. The normalization condition a_{1} + a_{2} +…+ a_{k} = 1 has to be satisfied, and k, the number of different diffusing populations, is as small as possible.
The evaluation of the CDF for different time lags (Δt = 1, 2, … 10) allowed us to estimate the individual diffusion coefficients (Schütz et al., 1997) (Figure 2—figure supplement 1). This analysis of the cumulative distribution function fits the experimental data with a model of Brownian diffusion of different populations. Further analysis of the mean square displacement of translocations and the step correlation was necessary to determine the nature of diffusion and validate or refuse the simple Brownian model independently for each protein.
Mean square displacement (MSD)
Request a detailed protocolWe first computed the mean square displacement of translocations for each individual trace j (MSD_{j}) of length n, weighted with the probability of misconnection previously described. The MSD_{j} is therefore:
where d_{i,i + t} is the translocation distance between the frames i and i + t. The MSD_{j} for individual traces was then computed for increasing lag times up to 10Δt (t = 1Δt, 2Δt, … 10Δt) where Δt is the experimental interframe time interval of 10.5 ms. We then calculated the average mean square displacement MSD(t) for t = 1Δt, 2Δt, … 10Δt as the mean of all the individual traces MSD_{j}(t) for all the trajectories that had a length of at least equal to t. Finally, error bars for each data point of the average MSD(t) were calculated as the 95% interval of confidence computed by bootstrap resampling of the population.
Numerical simulations
Request a detailed protocolIn order to validate our detection and tracking algorithms and analysis, we performed a series of numerical simulations. These simulations consisted in videos of particles with the optical response of our optical system, diffusing in 3D Brownian motion with a given diffusion coefficient. The signal was then corrupted with additional noise composed by a mixture of Gaussian and shot noise that mimicked our raw experimental data.
Parameters of the simulation
Request a detailed protocolThe PSF of the single particle signal was obtained with the PSF Lab software (Nasse and Woehl, 2010). The parameters used to retrieve the PSF were: emission wavelength 600 nm, objective NA 1.49, coverslip thickness 150 μm, oil refractive index 1.51, coverslip refractive index 1.52, and sample refractive index 1.3. The PSF was computed for a total height of 4 μm in layers of 100 nm and radius of 2 μm, on a pixelated image with pixel size of 107 nm. Intermediate values were estimated by linear interpolation of the eight surrounding points.
In order to estimate the noise, we analyzed the distribution of the pixel intensity values of experimental videos after removing the values of those pixels included in any detection. We then fitted this distribution to a combination of Gaussian and Poisson distributions, with a result of 95% Poisson distribution (λ = 20) multiplied by a factor determined by the camera gain, and 5% white noise.
The movement of the particles was simulated to take place in the interior of a closed box with similar dimensions to those of the eukaryotes nucleus: 10 μm × 10 μm × 6 μm. The particle density inside the box was set to be constant (i.e., the photobleaching rate and the photoactivation rate were the same), and therefore the ratio of disappearance and appearance of a new particle at a random position was set accordingly to the measured photobleaching halflife.
The video images were finally obtained as a convolution of the PSF with the displacement of the particle in the pixelated matrix during the acquisition time. We computed this by estimating the position (x, y, z) every 1 ms (10 estimations per frame) and by adding the convolved PSF at (x, y, z) to the final image. To take into account the displacement during the EMCCD transfer time between two consecutives images, an additional unrecorded movement of 0.5 ms was added to the simulation.
Reconstruction of diffusion
Request a detailed protocolWe then run the simulated videos of one single population of 3D Brownian diffusing particles through our detection and tracking algorithms. The histogram of translocations retrieved from the analysis of simulated films was in very good agreement with the theoretical values for diffusion coefficients between 0.1 µm²/s and 20 µm²/s (Figure 1—figure supplement 3A and 3B).
Diffusion coefficient boundaries
Request a detailed protocolIn our experiments, the minimum interframe displacement was limited by the experimental single molecule localization accuracy. The pointing error can be defined as the distance between the real centroid of the particle and the coordinates of the detection (Figure 1—figure supplement 3C). Using our simulations, we could determine the localization accuracy as the mean value of the pointing error, as a function of the diffusion coefficient. It was estimated to be ∼70 nm with a dramatic increase for particles with diffusion coefficient higher than 10 μm^{2}/s. This is in agreement with the experimental estimation of the localization accuracy retrieved from the consecutive detections of an immobile H2B molecule (Figure 1—figure supplement 1D). The lower bound of a detectable diffusion coefficient was thus ∼0.04 μm^{2}/s.
The analysis of simulated videos also allowed us to determine the percentage of detected particles. We could determine that the percentage of detections followed a Gaussianlike distribution along the optical axis, centered at the focal plane (Figure 1—figure supplement 3D). The width of such distribution determined the focal depth and it is in good agreement with the axial width of the PSF in our experimental conditions (∼600 nm). Moreover, the amplitude of the detection distribution was dependent on the diffusion coefficient of the particles. There is a higher rate of detection for slow particles than for those with higher diffusion coefficient. This effect could also be observed in the in vivo data by plotting the averaged 1Δt displacement as a function of the duration of the trajectory (Figure 2—figure supplement 2A, for Dendra2). Less mobile particles were detected for longer periods of time.
Population exclusion
Request a detailed protocolThe dependency of the percentage of detected particles with the diffusion coefficient of the particles needs to be taken into account when analyzing the mobility of a heterogeneous mixture of molecules with different diffusion coefficients. We performed simulations of an extreme case with 50% of the molecules following Brownian diffusion at 1 µm²/s and 50% at 10 µm²/s. As expected, fast particles had a higher probability of escaping the focal depth of observation between two consecutive frames, and therefore slow particles were overrepresented in the measurement (Figure 1—figure supplement 3C). Such exclusion of the fast diffusing particles population significantly affects population analysis as well as the average MSD. The population analysis of the one step translocation histogram gave a proportion rate of 80% of particles with D = 1 µm²/s and 20% D = 10 µm²/s. Despite the bias on the population, the values of diffusion coefficients were not affected by the population exclusion (Figure 2—figure supplement 2C). Similarly, the resulting MSD analysis of the simulated data showed a deviation from linearity, suggesting an apparent subdiffusive behavior of the ensemble of molecules (Figure 2—figure supplement 2D).
In order to take this bias into account in our analysis, we measured the number of translocations detected on single particle simulation videos with diffusion coefficients ranging from 0.001 µm²/s to 20 µm²/s (Figure 2—figure supplement 2B). We then used this information as the reference curve to correct the proportions of populations retrieved from the analysis of our experimental data.
In order to perform such correction, we considered an arbitrary fit of the one step translocation histogram with three populations: (a_{1}, D_{1}), (a_{2}, D_{2}) and (a_{3}, D_{3}), with a_{1}, a_{2}, and a_{3} representing the fractions of populations and D_{1}, D_{2}, and D_{3} their diffusion coefficients. We found, interpolating the reference curve for each diffusion coefficient, the proportion of the population that was integrated in our study p(D_{1}), p(D_{2}), and p(D_{3}). For instance, when D = 1 µm²/s, we detected p(D) = 97% of the translocations. We then computed the corrected values for a_{1}, a_{2}, and a_{3} as $\frac{{a}_{1}}{\text{p}\left({D}_{1}\right)}$, $\frac{{a}_{2}}{\text{p}\left({D}_{2}\right)}$, and $\frac{{a}_{3}}{\text{p}\left({D}_{3}\right)}$ to obtain the relative rate of diffusive populations.
It is important to note that this correction assumes Brownian diffusion of the molecules, and therefore has to be understood as a first order correction of the population rates in all our experimental data. However, due to this selection bias, fast diffusing molecules, for which the MSD slope is the highest, contribute less to the average MSD at longer time lags. As a result, the average MSD observed for Dendra2 is consistent with a normal diffusive behavior for three species.
Simulations of models with intermittent regimes
Request a detailed protocolWith our experimental system well characterized, we simulated increasingly sophisticated models of diffusion with intermittent regimes. We first considered a model in which particles transition from fast to slow Brownian diffusion. We also considered a model in which particles switch from free diffusion to confined diffusion within a spherical trap. In both cases, we considered an infinite volume in the plane x–y, and to 1 µm in the axial direction. We then recorded the positions of simulated traces until they exited the volume. For each model, we simulated 100,000 trajectories whose recorded position was corrupted by a Gaussian curve with 70 nm standard deviation, simulating the experimental localization accuracy.
We first simulated regime switching between a fast (D1) and a slow (D2) diffusing coefficient. For each translocation, there was a probability Kon to transition from D1 to D2, and a probability Koff to transition from D2 to D1. We adjusted D1 and D2 to the diffusion coefficients obtained from the twopopulation fit of the experimental cumulative histograms of steps, for both cMyc and PTEFb. These values were D1 = 14 μm^{2}/s and D2 = 1 μm^{2}/s for cMyc, and D1 = 15 μm^{2}/s and D2 = 1 μm^{2}/s for PTEFb. Likewise, we determined the ratio Kon/Koff from the retrieved populations of the fit (Kon/Koff = 0.68/0.32 for cMyc, Kon/Koff = 0.54/0.46 for PTEFb). With these constrains, we found the association and dissociation rates that resulted in a translocation of histograms close to those obtained in the experiments (Figure 5—figure supplement 1).
We also investigated switching between a freely diffusing mode and a confined diffusion. The potential of confinement was set to ‘hardwall’ type, with random repositioning of the molecule in the volume in the case of two successive bouncing against the boundaries. A new confinement was created anytime the particle switched between a freely diffusing and confined mode.
With the intermittent diffusion model, we retrieved mild negative values of the AC. These were prominently at small lag times and small translocation steps. The tendency for larger times and steps was always toward AC = 0, that is a symmetric angular distribution (Figure 5—figure supplement 2). We concluded that such scenario could not reproduce the antipersistent characteristics of our data. Conversely, we obtained stronger angular asymmetries with the intermittent trap model. The AC curves also evolved with time and space and were strongly dependent of the size of the trap. Expectedly, the simulated particles undergo confined diffusion only when they are inside the trap. In conclusion, although we reproduced certain characteristics of the angular distribution of our data, these simple simulations failed to reproduce the AC temporal and spatial dependences of cMyc, as well as the scale independent behavior of PTEFb.
Numerical simulations of anomalous diffusion models
Request a detailed protocolModels for subdiffusion
Request a detailed protocolSubdiffusion motion has been frequently reported in SPT experiments (Saxton, 2007). In cells, it is commonly attributed to one of the following two processes: a broad distribution of trapping times or an obstructed movement due to crowding effects. In our experiments on PTEFb, we could rule out the former, often referred to as the continuous time random walk model (CTRW) (Metzler and Klafter, 2000). We addressed this model by simulation to compute its angle distribution. In Figure 5—figure supplement 3, we show the result of a Monte Carlo simulation of a continuous time random walks performed on an infinite cubic lattice. The position was recorded every 1000 steps, and the waiting times were uncorrelated following a discretized heavytailed probability distribution,
with α set as 0.6. The MSD shown in Figure 5—figure supplement 3A is an ensemble MSD, averaged on 10,000 trajectories and rescales as a power law 〈r²(t)〉 ∼ t^{α}. It is noteworthy to point out that the timeaverage MSD of a CTRW realization does not result in a sublinear relationship with time (benAvraham and Havlin, 2005). The experimental MSD curves shown in Figure 2 and Figure 3 were averaged over time and also over the ensemble of all the trajectories, which was an additional indication against the CTRW model for our data.
Random walks on a fractal medium induce spatial correlation between successive displacements imposed by the self similarity of the geometry. For that reason, the fractal model has been applied to the comprehension of random walks on disordered media (Szymanski and Weiss, 2009). Fractal object involves the repetition of the same features of an object at different scales. If the whole object is repeated then the fractal is exact. One example of exact fractal is the Sierpinski gasket, a 2Dembedded fractal lattice of dimension df = log(8)/log(3) ≈ 1.89. An example of a nonexact fractal is the maximum site percolation cluster at percolation threshold. The geometry is not conserved in the cluster at different scales, but rather features such as the local density obey a power law.
We computed a fractal network as the maximum cluster of a cubic lattice at critical site percolation probability (benAvraham and Havlin, 2005). The initial cubic lattice dimensions were 2000 × 2000 × 500 sites. We then removed sites from the lattice according to the critical probability P_{c} = 0.311604 (benAvraham and Havlin, 2005). The size of the maximum cluster was 5,967,870 sites. We then performed random walks on such fractal structure by recording a position every 2000 steps on the lattice.
Angular distribution evaluation by Monte Carlo simulations
Request a detailed protocolSimulations of trajectories on cubic and fractal lattices were performed in order to obtain the angular distribution of consecutive steps, and their temporal evolution. The angular distribution was obtained by Monte Carlo simulations of 10,000 realizations of trajectories of 500 steps with a randomly distributed start. Since the initial mesh was a cube, there were privileged directions with higher number of sites and thus a higher number of possible successive positions. In computations of the angular distributions on such simulated trajectories, this results in an over representation of the directions imposed by the lattice geometry such a 90°angle that vanished with an increasing time lag. This bias would be negligible if the percolation cluster had a large number of sites, such that allowed us to record the position of the trajectory for at time lag significantly larger than 2000 steps. The limitation to perform such a simulation was the random access memory of the computer. We therefore applied a simple correction to the angular distribution of the simulated trajectories. We recorded all the possible translocations from our simulations for time lags 1Δt to 10Δt. These translocations were then shuffled to compute the distribution of the angles that was inherent to the network itself and not to the successive displacement correlation. Such ‘structural distribution’ reflected therefore the anisotropy of the structure due to the finite scale of the computations, and we used it as a normalization distribution. We verified that this ‘structural’ angular distribution flattened at increasing time lag. The angular distributions shown in Figure 5 were therefore rescaled by dividing each bin proportion by the corresponding one in the structural distribution.
References

Quantitative study of single molecule location estimation techniquesOptics Express 17:23352–23373.https://doi.org/10.1364/OE.17.023352

A fractal model for nuclear organization: current evidence and biological implicationsNucleic Acids Research 40:8783–8792.https://doi.org/10.1093/nar/gks586

BookDiffusion and reactions in fractals and disordered systemsCambridge University Press.

Facilitated diffusion of proteins on chromatinPhysical Review Letters 106:038102.https://doi.org/10.1103/PhysRevLett.106.038102

Distribution of directional change as a signature of complex dynamicsProceedings of the National Academy of Sciences of the United States of America 110:19689–19694.https://doi.org/10.1073/pnas.1319473110

Quantitative Comparison of algorithms for tracking single fluorescent particlesBiophysical Journal 81:2378–2388.https://doi.org/10.1016/S00063495(01)758845

Firstpassage times for random walks in bounded domainsPhysical Review Letters 95:260601.https://doi.org/10.1103/PhysRevLett.95.260601

Probing microscopic origins of confined subdiffusion by firstpassage observablesProceedings of the National Academy of Sciences of the United States of America 105:5675–5680.https://doi.org/10.1073/pnas.0712158105

The dynamic range of transcriptionMolecular Cell 30:545–546.https://doi.org/10.1016/j.molcel.2008.05.009

Imaging transcription in living cellsAnnual Review of Biophysics 38:173–196.https://doi.org/10.1146/annurev.biophys.050708.133728

Kinetics of diffusioncontrolled processes in dense polymer systems. I. Nonentangled regimesThe Journal of Chemical Physics 76:3316.https://doi.org/10.1063/1.443328

Kinetics of diffusioncontrolled processes in dense polymer systems. II. Effects of entanglementsThe Journal of Chemical Physics 76:3322.https://doi.org/10.1063/1.443329

MethylQuant: a realtime PCRbased method to quantify DNA methylation at single specific cytosinesMethods in Molecular Biology 507:281–303.https://doi.org/10.1007/9781597455220_21

Singlemolecule investigations of the stringent response machinery in living bacterial cellsProceedings of the National Academy of Sciences of the United States of America 108:E365–E373.https://doi.org/10.1073/pnas.1102255108

Transcriptional regulation by coactivators in embryonic stem cellsTrends in Cell Biology 22:292–298.https://doi.org/10.1016/j.tcb.2012.04.002

The hierarchy of the 3D genomeMolecular Cell 49:773–782.https://doi.org/10.1016/j.molcel.2013.02.011

Crumpled globule model of the threedimensional structure of DNAEurophysics Letters 23:373–378.https://doi.org/10.1209/02955075/23/5/012

Probing intranuclear environments at the singlemolecule levelBiophysical Journal 94:2847–2858.https://doi.org/10.1529/biophysj.107.115014

The ins and outs of gene regulation and chromosome territory organisationCurrent Opinion in Cell Biology 19:311–316.https://doi.org/10.1016/j.ceb.2007.04.016

Ultrahigh resolution imaging by fluorescence photoactivation localization microscopyBiophysical Journal 91:4258–4272.https://doi.org/10.1529/biophysj.106.091116

Anomalous transport in the crowded world of biological cellsReports on Progress in Physics 76:046602.https://doi.org/10.1088/00344885/76/4/046602

Validating transcripts with probes and imaging technologyNature Methods 8:S12–S19.https://doi.org/10.1038/nmeth.1573

Rateprocesses on fractals  Theory, simulations, and experimentsJournal of Statistical Physics 42:185–200.https://doi.org/10.1007/BF01010846

Gene location and DNA density determine transcription factor distributions in Escherichia coliMolecular Systems Biology 8:1–13.https://doi.org/10.1038/msb.2012.42

A single molecule view of gene expressionTrends in Cell Biology 19:630–637.https://doi.org/10.1016/j.tcb.2009.08.008

Scaling properties of diffusionlimited reactions on fractal and euclidean geometriesJournal of Statistical Physics 65:1269–1283.https://doi.org/10.1007/BF01049612

Photoactivatable fluorescent proteins for diffractionlimited and superresolution imagingTrends in Cell Biology 19:555–565.https://doi.org/10.1016/j.tcb.2009.09.003

Singlemolecule Spectroscopy and imaging of Biomolecules in living cellsAnalytical Chemistry 82:2192–2203.https://doi.org/10.1021/ac9024889

A benchmark for chromatin binding measurements in live cellsNucleic Acids Research 40:e119.https://doi.org/10.1093/nar/gks701

The Random walk's guide to anomalous diffusionThe Random walk's guide to anomalous diffusion.

GeometryInduced Bursting dynamics in gene expressionBiophysical Journal 102:2186–2191.https://doi.org/10.1016/j.bpj.2012.03.060

FRAP and kinetic modeling in the analysis of nuclear protein dynamics: what do we really know?Current Opinion in Cell Biology 22:403–411.https://doi.org/10.1016/j.ceb.2010.03.002

Realistic modeling of the illumination point spread function in confocal scanning optical microscopyJournal of the Optical Society of America A, Optics, Image Science, and Vision 27:295–302.https://doi.org/10.1364/JOSAA.27.000295

Intranuclear mobility and target search mechanisms of transcription factors: a singlemolecule perspective on gene expressionBBA–Gene Regulatory Mechanisms 1819:482–493.https://doi.org/10.1016/j.bbagrm.2012.02.001

Half a century of “the nuclear matrix”Molecular Biology of the Cell 11:799–805.https://doi.org/10.1091/mbc.11.3.799

A biological interpretation of transient anomalous subdiffusion. I. Qualitative modelBiophysical Journal 92:1178–1191.https://doi.org/10.1529/biophysj.106.092619

Wanted: a positive control for anomalous subdiffusionBiophysical Journal 103:2411–2422.https://doi.org/10.1016/j.bpj.2012.10.038

Singlemolecule microscopy on model membranes reveals anomalous diffusionBiophysical Journal 73:1073–1080.https://doi.org/10.1016/S00063495(97)781396

Analysis of binding reactions by fluorescence recovery after photobleachingBiophysical Journal 86:3473–3495.https://doi.org/10.1529/biophysj.103.026765

Elucidating the origin of anomalous diffusion in crowded fluidsPhysical Review Letters 103:038102.https://doi.org/10.1103/PhysRevLett.103.038102

Precise nanometer localization analysis for individual fluorescent probesBiophysical Journal 82:2775–2783.https://doi.org/10.1016/S00063495(02)75618X

Identifying specific protein interaction partners using quantitative mass spectrometry and bead proteomesThe Journal of Cell Biology 183:223–239.https://doi.org/10.1083/jcb.200805092

Subnuclear segregation of genes and core promoter factors in myogenesisGenes & Development 25:569–580.https://doi.org/10.1101/gad.2021411

SingleRNA counting reveals alternative modes of gene expression in yeastNature Structural & Molecular Biology 15:1263–1271.https://doi.org/10.1038/nsmb.1514

RNA polymerase ii elongation controlAnnual Review of Biochemistry 81:119–143.https://doi.org/10.1146/annurevbiochem052610095910
Article and author information
Author details
Funding
Agence Nationale de la Recherche (L' Agence Nationale de la Recherche) (PCV DynaFT)
 Raphaël Voituriez
 Olivier Bensaude
 Maxime Dahan
 Xavier Darzacq
Nikon France
 Ignacio Izeddin
 Maxime Dahan
 Xavier Darzacq
Netherlands Organisation for Scientific Research (NWO) (Rubicon)
 Ignacio Izeddin
Fondation pour la Recherche Médicale (Foundation for Medical Research in France)
 Vincent Récamier
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Dan Larson, John Lis, Florian Mueller, Yitzhak Rabin, Robert Singer, and Robert Tjian for discussions and comments. Leonid Mirny for an in depth review and his challenge and help with the modeling approach. We are grateful to Sarah Moorehead and Mohamed El Beheiry for critical reading and to Daniel Ciepielewski and Philippe Rideau for discussions on the microscopy. VR acknowledges financial support from FRM, and II from NWO. Work presented here was supported by the PCV DYNAFT 08PCVI0013 and DynamIC ANR12BSV50018 from Agence Nationale pour la Recherche, and a research contract with Nikon France to XD and MD.
Version history
 Received: January 7, 2014
 Accepted: June 11, 2014
 Accepted Manuscript published: June 12, 2014 (version 1)
 Version of Record published: July 15, 2014 (version 2)
Copyright
© 2014, Izeddin et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics

 10,920
 views

 1,434
 downloads

 233
 citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading

 Biochemistry and Chemical Biology
 Structural Biology and Molecular Biophysics
NADPH oxidases (NOX) are transmembrane proteins, widely spread in eukaryotes and prokaryotes, that produce reactive oxygen species (ROS). Eukaryotes use the ROS products for innate immune defense and signaling in critical (patho)physiological processes. Despite the recent structures of human NOX isoforms, the activation of electron transfer remains incompletely understood. SpNOX, a homolog from Streptococcus pneumoniae, can serves as a robust model for exploring electron transfers in the NOX family thanks to its constitutive activity. Crystal structures of SpNOX fulllength and dehydrogenase (DH) domain constructs are revealed here. The isolated DH domain acts as a flavin reductase, and both constructs use either NADPH or NADH as substrate. Our findings suggest that hydride transfer from NAD(P)H to FAD is the ratelimiting step in electron transfer. We identify significance of F397 in nicotinamide access to flavin isoalloxazine and confirm flavin binding contributions from both DH and Transmembrane (TM) domains. Comparison with related enzymes suggests that distal access to heme may influence the final electron acceptor, while the relative position of DH and TM does not necessarily correlate with activity, contrary to previous suggestions. It rather suggests requirement of an internal rearrangement, within the DH domain, to switch from a resting to an active state. Thus, SpNOX appears to be a good model of active NOX2, which allows us to propose an explanation for NOX2’s requirement for activation.

 Cell Biology
 Structural Biology and Molecular Biophysics
Acetylation of αtubulin at the lysine 40 residue (αK40) by αTAT1/MEC17 acetyltransferase modulates microtubule properties and occurs in most eukaryotic cells. Previous literatures suggest that acetylated microtubules are more stable and damage resistant. αK40 acetylation is the only known microtubule luminal posttranslational modification site. The luminal location suggests that the modification tunes the lateral interaction of protofilaments inside the microtubule. In this study, we examined the effect of tubulin acetylation on the doublet microtubule (DMT) in the cilia of Tetrahymena thermophila using a combination of cryoelectron microscopy, molecular dynamics, and mass spectrometry. We found that αK40 acetylation exerts a smallscale effect on the DMT structure and stability by influencing the lateral rotational angle. In addition, comparative mass spectrometry revealed a link between αK40 acetylation and phosphorylation in cilia.