Visualization of barriers and obstacles to molecular diffusion in live cells by spatial pair-cross-correlation in two dimensions

: Despite recent advances in optical super-resolution, we lack a method that can visualize the path followed by diffusing molecules in the cytoplasm or in the nucleus of cells. Fluorescence correlation spectroscopy (FCS) provides molecular dynamics at the single molecule level by averaging the behavior of many molecules over time at a single spot, thus achieving very good statistics but at only one point in the cell. Earlier image-based methods including raster-scan and spatiotemporal image correlation need spatial averaging over relatively large areas, thus compromising spatial resolution. Here, we use spatial pair-cross-correlation in two dimensions (2D-pCF) to obtain relatively high resolution images of molecular diffusion dynamics and transport in live cells. The 2D-pCF method measures the time for a particle to go from one location to another by cross-correlating the intensity fluctuations at specific points in an image. Hence, a visual map of the average path followed by molecules is created.


Introduction
Signal transduction at a sub-cellular level is a key element in the regulation of cell functions, cell fate and cell migration. In the compartmentalized environment of cells these transduction processes are regulated by a tight spatio-temporal control of molecular flow [1,2]. The regulation of this molecular flow requires a complex interplay between the structure of the surrounding cellular components [3][4][5] and a spatio-temporal control of molecularinteractions [6,7]. The consequence is a molecular transport that depends strictly on the molecular species, on the aggregations state, on the active state of the molecule and on the cell-cycle phase [8][9][10][11][12]. Despite the pivotal role of molecular flow coordination in the living cell, we lack systematic approaches to visualize the complex connectivity of the cell interior with high resolution and, at the same time, on the proper temporal scale. This lack of methods to visualize molecular connectivity is a significant limit for advancing our knowledge in signal transmission. Here, we developed algorithms based on the pair correlation approach to display 2-D molecular connectivity maps of the cell interior. Molecular diffusion and transport processes are fundamental in physical, chemical and biological systems. Most of the methods proposed to measure diffusion in cells are based on fluorescence correlation spectroscopy (FCS) and provide only local information about the diffusivity of molecules in cells and tissues. Fluctuation-based techniques, FCS, STICS (Spatio-Temporal Image Correlation Spectroscopy) and RICS-(Raster-Scan Image Correlation Spectroscopy) are able to measure the diffusion coefficient at single pixels (FCS) or in a small region of a cell (STICS and RICS). Eventually, by measuring many points in sequence, FCS can produce a diffusion map of the cell interior; however, these approaches cannot measure the spatial communication between regions of the cell at the scale of an entire cell. Although local or average values are needed for estimating protein concentration and mobility, we cannot determine from these measurements the path followed by molecules after activation or upon interaction with molecular partners. Here we propose a method to detect molecular transport by measuring the cross-correlation of fluctuations for all the pairs of points in the sample. To this end, we expand the pCF (pair Correlation Function) concept [13] to 2-dimensional images (2D-pCF). The pCF function is defined by among cell compartments [5, 8-10, 14, 15]. Using the pCF method along a line, we were able for the first time to describe the effect of chromatin on the flow of GFP in live cell nuclei [5,8,9]. We also measured fluxes between the cytoplasm and the nucleus and the delayed communication caused by the passage through the nuclear pore complex [7,16]. Other successful applications of the pCF approach have been reported by other labs [11,[17][18][19][20]. However, all previous work was aimed at detecting the existence of correlated fluctuations between two distant points without constructing an image of all correlated points. If there is an obstacle in between the two points molecules will go around and it will take them longer to appear at the other location. Thus, 2D-pCF (Eq. (1)) detects barriers to diffusion and heterogeneity because the time of the correlation maximum is delayed in the presence of diffusion barriers. In a nutshell, this non-invasive, sensitive technique can follow the same molecule over a large area producing a map of molecular diffusion. It does not require isolated molecules; thereby, many molecules can be labeled with the same fluorophore and be within the same volume of excitation (the PSF). The 2D-pCF approach is a change of paradigm from both a technological and a conceptual point of view. From the technological perspective we need very high resolution spatial sampling coupled with a large field of view, and very high time resolution coupled with high sensitivity detection in order to measure the fluctuations due to fast moving and dim molecules in a whole cell. Currently, this can be achieved using fast and highly sensitive cameras. Then, from a conceptual point of view, we need an efficient algorithm able to visualize the connectivity between ~10 5 -10 6 points in one slice of the cell interior, i.e., crosscorrelate the fluctuations between all pairs of points of large images collected for a certain period of time. Although many investigations have focused on the study of the dynamical properties of proteins in live cells, most of the results are local, purely temporal or spatially averaged and in all these cases spatial communication is not visualized [11,[21][22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37]. Also, many of the current approaches require an interpretative model and a consequent analytical relationship between the measured properties (usually a change in fluorescence intensity) and the physical properties of the particle motion (e.g., diffusion coefficient). In order to unveil something unknown, such as the spatial communication in the cell, a model-free approach is required.
In this article we show that average molecular trajectories can be measured by the 2D-pCF approach without specific single molecule localization algorithms as used in single particle tracking. The experimental proof of these concepts shows that the 2D-pCF can be applied in live cells using available fluorescent chimeras of functional proteins and thus utilizes common techniques for protein labeling and expands our ability to quantify key players in cellular flow. Furthermore, the practical demonstration of connectivity maps may result in a breakthrough in the study of biochemical communication in live cells, by unveiling the role of the cell spatial organization in signal spreading, as well as describing the route followed by the signal molecules to spread from the source to the target.

Identifying barriers to diffusion by the 2D-pCF approach
In a spatially heterogeneous environment, the local diffusion coefficient is not sufficient to accurately describe the transport of molecules to neighboring regions. Penetrable or impenetrable barriers can shape the trajectories of molecules, guiding and confining molecular diffusion. Fig. 1. Pictorial representation of the principle of the 2D-pCF to detect molecular transport in cells guided by obstacles and molecular interactions. A) Camera image time stack of the same focal plane (5,000 to 50,000 images) giving access to a grid of points where fluctuations due to the underlying molecular flow can be measured. B) The 2D-pCF is computed between every pair of points of the image and at a given distance. In an isotropic environment. Starting from a point, all the correlation functions should be equal and independent of the direction. The pCF correlation functions can have a maximum. C) If the environment is anisotropic, molecules could take longer to reach a given distance depending on the direction. D) The polar image resulting from plotting all the correlation functions as a function of the angle (2D-pCF) is analyzed using the first and second central moment (Eqs. (2)- (7)). E) The 2D-pCF is plotted using only the long axis of the ellipses shown in D. The definition of the axes lengths is given in Eq. (5). F) By joining the segments corresponding to the long axis we obtain the connectivity map that delineates the presence of barriers shown in red. G) If an underlying structure is present, shown schematically by a maze, the spatial connectivity can be obtained using different radii for the calculation of the 2D-pCF. H) The 2D-pCF function at a distance of 3 pixels for this maze reveals connected points at a small scale. I) The 2D-pCF calculated at a larger distance (4 pixels) show paths connected at a larger scale in the maze, and where there is an obstacle (center of the maze). The distance for the calculation of the 2D-pCF is used to reveal different barriers.
As shown in Fig. 1(A), we start with a stack of images of the same focal plane in the cell. From this stack we isolate a point and then we produce a set of correlation functions 2D-pCF(δr) at a given distance, δr, at different angles around a given point. This set for one particular distance δr can be plotted in a polar plot as a heat map where the hotter colors indicate larger values of the 2D-pCF function ( Fig. 1(B)). We used a log scale for the delay time variable (Fig. 1(B)), as this is common practice in FCS. Note that the 2D-pCF is calculated at a given distance which can be varied producing visualizations of the 2D-pCF that show the path of molecules that travel long distances. For example, if there is an obstacle at a given distance from the original point r 0 , for shorter distance we could observe changes of the correlation function at given angles that indicate the presence and local direction of the obstacle. But if we look at a longer distance we could observe that a particle has reached a distant position at much later times going around the obstacle ( Fig. 1(G)-1(I)). Figure 2 shows 2D-pCFs calculated at specific points in the image and the algorithm for the calculation of parameters sufficient to describe the lack of symmetry of the 2D-pCF ( Fig.  2(J)). (1) using 32 angles at a distance of 4 pixels. B) Intensity of a single image. C) Average intensity of all images D) 2D-pCF at one pixel from the right border of the image at point 1. E) 2D-pCF in a region far from the borders at point 2. F) 2D-pCF at 2 pixels from the left border at point 3. G) 2D-pCF at 2 pixels below the channel, point 4. H) 2D-pCF in the channel at point 5. I) 2D-pCF at 1 pixel above the channel, point 6. The presence of the discontinuity due to the channel can be seen in panels G, H and I although the walls of the channel are invisible since they are not labeled. J) Computed analysis of the shape of the 2D-pCF according to Eqs.
(2)-(7) using the pattern of panel H. For this simulation the PSF was 300 nm, the pixel size was 114 nm and the total field of view was 7.2 µm. The image has 64x64 pixels and 32,768 images were simulated; the particle density for this simulation was 10 molecules/μm 2 , which is a low value compared to experimental situations. Note that the 2D-pCF gives different patterns by moving the center of the calculation by one pixel as shown in panels D and F. Panels D-J show an image representing the 2D-pCF in a polar plot. The axes are in log-delay time. The maximum delay time is equal to 164s.
In Fig. 2 we show pCF polar maps for a simulated situation. We use image processing methods to obtain parameters of the angular distribution of the pCF. For this article, we use the first and second central moments of the spatial distribution; higher moments are ignored.
For each point and for each distance we describe the angular distribution using few parameters ( Fig. 2(J) and Eqs. (2)- (7)). 11 20 02 2 1 arc tan 2 20 02 20 02 I(x,y) in Eq. (2) is the value of the 2D-pCF correlation function at pixels shifts x and y. M ij is the definition of the moments of the distribution of the 2D-pCF. The shift of the center of mass of the distribution (the μ parameters at the first order) can be used to obtain the net velocity of the particle if the evaluation of the pCF is done at multiple distances. From Eq.
(3), the displacement of the center of mass is expressed in the unit of pixels of the image. In Eq. (5), the values of the parameters λ correspond to the long (λ 1 ) and short (λ 2 ) axis of the ellipses shown schematically in Fig. 2(J). In Eq. (5) the units of the length of the ellipses axes (second order moments) are in pixels square and they indicate how much the pCF distribution expands. As previously discussed [5], the pCF is larger and the peak of the pCF will appear at shorter time for faster diffusion. Using this principle, the difference between the long axis and the short axis of the ellipses indicate that in the direction of the long axis the pCF distribution is expanding faster. In this work, we use the long and short axis of the ellipses only as an indication of the extent of the anisotropy of the diffusion and of the direction of the preferential anisotropy. For this reason we use the ratio (or normalized difference) of the long and short axis of the ellipses to construct a map indicating where the pCF distribution is expanding anisotropically. To summarize, the second order moments are analyzed in terms of the anisotropy or eccentricity of the angular distribution (Eqs. (6) and (7)) which are related to the size of the major and minor axis of the ellipses describing the second moment distribution and the angle of the major axis of the ellipses (Eq. (4)) that indicates the local direction of the diffusional motion. The analysis software provides six maps, the intensity at each point, the values of the long and short axis of the ellipses and the angle of the ellipses, the values of the displacement of the center of mass and the average value of the pCF. The long and short axes of the ellipses are used to calculate the eccentricity (or anisotropy) of the diffusion, the two components of the position of the center of mass are used to calculate the direction and the modulus of the displacement of the center of mass of the pCF distribution.
Along each line as a function of the angle, if we change the distance, the time of the maximum of the pCF correlation function could be related to the mean square displacement (MSD) since the pCF represents the time it takes for a particle to reach a certain distance (Fig.  4). From the MSD we could obtain local diffusion coefficients and other parameters such as local confinement or the coefficient of anomalous diffusion. This calculation is not used in this article since the derivation of the diffusion coefficient form the slope of the MSD plot as a function of time delay assumes that the motion is isotropic, which is not the case for the data on cells used in this work. The calculation of the diffusion tensor has recently been presented and we are not repeating here the major findings of that paper [38]. Instead, we will discuss here how to obtain high resolution connectivity maps of the cell interior. The analysis using the moments of the distribution of the radial 2D-pCF function shown in Fig. 2(J) can be used as the set of rules for the classification of the molecular movement in different parts of an image to produce a map of obstacles and a map of the regions where the space is connected. Figure 3 shows the analysis of the channel simulation used in Fig. 2 in which molecules are free to diffuse in a plane in which there is a 300-nm wide channel that simulates the presence of guided diffusion. We applied the 2D-pCF algorithms described in the discussion of Fig. 2 and in Appendix at each pixel of the simulation shown in Fig. 3. The intensity is relatively uniform since the particle density was about the same in the region of simulation (64x64 pixels) as shown in Fig. 3(A). The calculation of the eccentricity (Eq. (6)) is done by first calculating the central moments of the radial 2D-pCF function. For this figure the calculation was done at four pixels distance. The eccentricity is plotted as a heat map on a scale of 0 to 1, where 1 corresponds to the hottest color. For the connectivity map the threshold for the eccentricity was set at 0.45 ( Fig. 3(D)). Only points close to the channel or inside the channel have larger eccentricity than the threshold value. For these points a segment of 10 pixels length is drawn in the direction of the ellipses as described in the Appendix. Figure 3(B) shows that only pixels near the channel wall or inside the channel display anisotropic diffusion. Note that the motion of particles is highly anisotropic in the channel and in the proximity of the channel. Only points with an eccentricity larger than 0.3 are plotted but Fig.  3(B) shows that there is a region outside the cannel where the motion has different degrees of eccentricity. This is because in Eq.    the average eccentricity was 0.31 ± 0.08 (mean ± SD). B) From the same data, the anisotropy was calculated according to the equations Eqs. (2)(3)(4)(5)(6)(7) in the manuscript. The histogram is shown in the middle panel (0.06 ± 0.03). This shows that, even in the absence of barriers, obstacles or directed motion, we get an eccentricity value >0. This is because in Eq. (5) the values of λ 1 is always greater than λ 2 . This value of 0.3 obtained in the case of random motion in the absence of barriers was used as lower threshold for regions of random motion for data visualization. C) To demonstrate that the fluorescence fluctuations were caused by free diffusion of EGFP in solution, the iMSD of the same data was calculated as shown on the right hand side, the resulting diffusion coefficient was 114 µm 2 /s, which is in within the range of values found for EGFP in solution. This example shows that the data used for the 2D-pCF could also be used far from obstacles to measure the diffusion coefficient.

Simulations of molecules diffusing in a plane with a 300 nm channel
In Fig. 5 we show a simulation of particles moving in and around a relatively small square box. The side of the box is 16 pixels which in the simulation (a pixel is 50 nm) correspond to a size of 0.8 μm, comparable to structures found in the cell nucleus or large vesicles. Figure 5 shows that the connectivity pattern caused by a small box has lines that depart from the surface. This is due to the high curvature of the structure simulated. This simulation shows that the simple rules used to explain the eccentricity patterns illustrated in Fig. 2 are valid for borders of low curvature but for highly curved structures additional rules must be considered. The pattern is fully reproducible in the simulation. This is not a random speckle but it is what is expected for the diffusion anisotropy close to corners and at points of high curvature. In principle, we can simulate any structure in both 2D and 3D that will help to recognize the observed patterns. The pattern simulated in Fig. 5 will be discussed later in conjunction with patterns of connectivity found in live cells. Fig. 5. Simulation of molecules moving in a plane where there is an impenetrable box. The size of the box is 0.8 μm x 0.8 μm, comparable to structures found in the cell nucleus and large vesicles. The particle density for this simulation is 8.7 molecules/μm 2 . Local diffusion coefficient is 10 μm 2 /s both inside and outside the box. The box is impenetrable. The 2D-pCF was calculated for a distance of 4 pixels. In the simulation the waist of the PSF was 200 nm. The camera frame rate was set to 100 frames/s and a total of 16,384 frames were analyzed. A) The eccentricity map obtained using the algorithms in the manuscript. The noise is due to the low particle density. B) Connectivity map using a threshold of 0.45 for the eccentricity values to be included in the connectivity map. Note that in this structure and at regions of high curvature there are connectivity lines perpendicular to the box side in addition to connectivity lines parallel to the box sides. C) Eccentricity histogram color coded using the same colors in the eccentricity map.

EGFP molecules moving in artificial grooves
In order to measure the fluctuations on a grid of points, the first requirement is fast frame acquisition of at least one optical section of the 3D intracellular space. By the classical rasterscan confocal approach the acquisition time increases linearly with the number of pixel measured and thus the sampling time required to measure the transport of small molecules (in the ms range) can only be achieved using a very fast frame rate scanner or in one dimension (a line scan). If the sample is inherently flat (<1 µm), there is no need for optical sectioning and a standard epifluorescence illumination with fast, camera-based image acquisition is adequate. As a first experimental verification of the 2D-pCF method for 2D diffusion we prepared a surface with periodic grooves 450 nm deep and 450 nm wide by stamping this pattern into a thin layer of PMMA on top of a glass cover slip. An electron microscopy image of the resulting pattern is shown in Fig. 6(A). To close the channels, a second cover slip was pressed on top of the grooves. Subsequently, the channels were filled with a 30 nM solution of monomeric EGFP and the sample was subjected to epifluorescence imaging with a camera frame rate of 100 Hz. Clearly, the channels exhibit a higher fluorescence intensity indicated by hotter colors in Fig. 6(B). Application of the 2D-pCF algorithm to the image time series yield the eccentricity map shown in Fig. 6(C). All channels show the same directionality indicating guided diffusion as we expected from this pattern. The direction of the molecular flow along the channels was visualized plotting the connectivity map ( Fig. 6(D)) as described in the Appendix. Fig. 6. Measuring the connectivity maps in testing samples. A 30 nM solution of EGFP is spread on top of a slide with grooves 450 nm deep and 450 nm wide. A cover glass is on top of the sample to maintain the solution in the grooves. The image was acquired in a wide field microscope equipped with an EMCCD camera at 100 frames/s; a total of 10,000 frames was collected for a total time of 100 s. The 2D-pCF algorithms were used to calculate the eccentricity and the connectivity maps. A) EM images of the grooves. B) The grooves containing a 30 nM EGFP solution. C) Eccentricity map. D) Connectivity map. EGFP molecules only diffuse in the groove as shown by the connectivity map. E) Intensity image of a GUV membrane labeled with DiOC16 in the sideSPIM microscope equipped with a sCMOS camera. F) Eccentricity map. G) Connectivity map. For the connectivity maps the eccentricity threshold was 0.45.

Visualizing the connectivity map in biological samples
In order to apply the 2D-pCF method to biological samples, we will take advantage of the widefield and SPIM microscope methods available in many laboratories today. For example, the SPIM microscope uses a low NA objective to produce a single plane of excitation light in the sample and a second (high NA) objective perpendicular to the first to collect the fluorescence light produced in the sample (see Figs. 8(A), 8(B)). By measuring the emission with a camera, the frame acquisition time in this microscope configuration can be pushed to the millisecond range. In the classical SPIM setup, the light sheet was obtained by placing a cylindrical lens in the excitation light path [39]; in the SPIM we built in our lab instead, the light sheet is obtained by rapid scanning of a focalized cone of light. This detail in the design of the microscope holds several advantages. For the purpose of this paper it is important that we can control the width of the light sheet by limiting the scan range. Moreover, the real exposure time of the sample can be controlled independently from the camera readout, avoiding the blurring effect caused both by particle motion during the exposure time and charge transfer during camera readout in rolling shutter mode. In fact, both EMCCD and sCMOS based detector arrays apply readout protocols that are "shutter free" in order to improve the read out speed. In other words, during shifting of charges produced by the incoming light in the sensor, the sensor is still exposed to light. Thus, a continuous exposure will introduce a trail in the image. Moreover, particularly in the case of an EMCCD camera, in order to decrease the acquisition time below 1 ms in the readout protocol, the photoelectrons are dynamically stored in a region of the chip that is outside the region of interest (ROI) but exposed to the incoming light. In order to avoid image overlapping, the chip outside the ROI has to be covered. Using a scanned beam, we can simply limit the size of the light sheet. The overall performance of our setup for detecting fast molecular dynamics in live cells has been already demonstrated in a recent publication [40].

Diffusion of DiOC16 in a GUV measured with the sideSPIM microscope
So far, we have shown connectivity maps for measurements done for simulated molecules moving in a channel and EGFP molecules diffusing in physical grooves where we had a clear expectation of the path followed by the molecules. In both cases the measured connectivity map was consistent with the expected pattern. Another model system where we can easily predict the direction of the diffusion are giant unilamellar vesicles (GUVs). We prepared GUVs labeled with dye molecules that are diffusing in the membrane so that the connectivity should only occur in the membrane itself (Figs. 6(E)-(G)). Figure 6(G) shows that, indeed, there is connectivity only along the membrane bilayer. The data for the GUVs experiments were collected with the sideSPIM microscope [41].

U2OS cells labeled with DiOC18 measured on a widefield microscope using epiillumination
In this measurement, the dye is labeling all membranes of the cell as shown in Fig. 7(A). The eccentricity map (Fig. 7(B)) shows strong anisotropic motion associated with membrane features and the connectivity map (Figs. 7(C), 7(D)) shows that the diffusive pattern is mainly along the membrane features. Note that each membrane seems to show two regions of anisotropic diffusion similarly to what was obtained in the simulations of diffusing molecules in the presence of a channel ( Fig. 7(C)).  Fig. 5 for a simulation of molecules moving around or inside a box that show that connectivity lines departing from a surface can be observed if the surface has high curvature. J) At pixel resolution the image shows the diffusion of the dye along the membrane at the junction between the two cells. K) The same junction explored at a longer distance (pCF at 8 pixels). L) Schematic representation of the flow expected at the two cell junction and plasma membrane where red is the barriers and the white lines indicate the connectivity. M) Envelope of histograms of eccentricity for pCF(4) and pCF (8). The vertical line indicates the threshold used for the connectivity maps shown in C and D. For this experiment 10,000 images were collected within a total time of 100 s. The background colors in the zoomed images G,H,J and K are the values of the anisotropy.
The difference in the connectivity maps calculated at a distance of 4 and 8 pixels is notable and it shows that the motion in the membranes is very local. The white arrow in Fig.  7(A) indicates a possible obstacle the dye molecules have to move around. Figures 7(G), 7(H) show that the 2D-pCF can reveal the occurrence of barriers and obstacles. At 4 pixels pCF distance, the two bright spots look connected and at 8 pixels distance the connectivity map strongly supports this hypothesis. Moreover, the comparison between pCF(4) and pCF (8) shows the occurrence of an impenetrable barrier for the dye. These results can be interpreted as the occurrence of a non-lipidic barrier around of a piled up membrane structure. The analysis of Figs. 7(J), 7(K) yield a clear result, supporting diffusion of the dye along the membrane at the cells junction. It is possible to identify that the dye cannot cross this section and this result is more evident comparing the results in pCF(4) vs pCF(8) pixels distances.

Elucidating obstacles and barrier in vivo of a cell expressing EGFP
Next, we measured EGFP diffusing in cells obtained with the sideSPIM microscope ( Fig.  8(A) 8(B)). Figure 8 shows the eccentricity of the diffusion and the connectivity map. In this case we do not have specific expectations about paths the EGFP molecules will follow in their diffusive motions detected by the algorithm for the pCF calculation other than that the pCF should be sensitive to barriers and obstacles for diffusion. By calculating the pCF at different distances we should be able to detect different paths for diffusion. In this case the system is fully three dimensional and we only measured an optical section of the cell. We note that in the intensity image ( Fig. 8(C)) there are no evident objects or barriers. Nevertheless, from the connectivity map ( Fig. 8(D)) we can see obvious obstacles for diffusion due to the external cell membrane and due to internal structures that shows up as molecules are going around an object as indicated by the white arrow ( Fig. 8(C)). By using longer distances for the computation of the 2D-pCF the occurrence of barriers that restrict the diffusion of the EGFP was evident and at the distance we used for the analysis these barriers were impenetrable (Figs. 8(D), 8(E)). A model for the obstacles highlighted by the white arrow is proposed in Fig. 8(F). This obstacle is located inside the nucleus, which suggests the presence of a nucleolus. Also, in Fig. 9, we show the nucleus contour using a white dotted line obtained from the CMS of the cell. The superposition of the connectivity map on the CMS image allowed us to highlight the obstacles and barriers inside the cell (Figs. 9(A), 9(B)). In Fig. 9(C) we show schematically that the presence of a barrier to diffusion results in a change of the center of mass of the angular distribution of the pCF function because the distribution is essentially one sided, so that the center of mass is shifted away with respect to the barrier.  (4) in a MB231 cell expressing EGFP. C) Schematic representation of the CMS concept. The pink stick represents an impenetrable barrier and the hemi circle represents the 2D-pCF. The CMS is defined as the displacement of the first order moment of the pCF distribution. The displacement of the center of mass of the 2D-pCF distribution will occur in the proximity of a barrier since the cross-correlation between on point in one side of the barrier and another point in the other part of the barrier will be zero.

Conclusions
In this paper we describe a method based on the 2D-pCF to produce high resolution maps of molecular movements in cells. We proceeded following three logical steps. First, we presented the 2D-pCF algorithms and discussed the algorithm for construction of the connectivity map using simulations starting from the 2D-pCF. Second, we tested the algorithm in physical systems where EGFP molecules were constrained to move in artificially produced grooves and in the membrane of a GUV. Third, we used cells with a membrane dye and cells expressing the enhanced green fluorescent protein (EGFP) which is known to have minimal interactions with cellular scaffolds and membranes. For a proof of principle, we first collected data of a thin portion of two cells touching each other using a wide field microscope with epi-fluorescence lamp illumination. Finally, we used SPIM microscopy where the spatio-temporal fluorescence fluctuation can be measured at any position in a thick sample due to the optical sectioning capability of this instrument with the required spatial and temporal resolution. In each case we obtained connectivity maps using the algorithm described in the method section showing barriers for diffusion and obstacles that force molecules to go around. The algorithm to measure the 2D-pCF in an image is the basis for visualizing the spatial connectivity map in the living cell.
For the 2D-pCF approach, we present preliminary simulations to demonstrate the algorithms that we are developing to visualize connectivity maps. There are two notable consequences of the definition we are using for the 2D-pCF: i) it contains the information on the diffusion tensor of the molecule and ii) it is sensitive to displacements smaller than the optical resolution. This is shown in Fig. 2(D) and 2(E) that differ by one pixel in the distance from the barrier. In the simulation the PSF was 300nm. A relevant result of this work is the demonstration that large frame sCMOS cameras for the SPIM microscope can be used as a sensitive detector for molecular fluctuations. While the EMCCD camera has a gain that is uniform and also the noise is uniform across the camera field, the sCMOS camera has a gain (sensitivity) that is different for each pixel. In addition, the noise is not Poissonian and pixels are correlated at a distance (mainly along rows and columns). Therefore, the use of these cameras for spatio-temporal fluctuations is not trivial. We have partially solved the problem by constructing a matrix of gain for each pixel and by avoiding bad pixel. Another way to deal with the noise of the sCMOS camera is to randomize the intensity of two adjacent pixels, which is the method used in this paper. The correction routine is part of the software to calculate the 2D-pair correlation function. For our sCMOS camera we have analyzed the dark noise and we have found that some pixels have a large background. We marked these pixels as bad so that they are not used in the calculations. In addition, the sCMOS camera has correlations between lines, likely due to the process of reading the camera. Before calculation of the 2D-pCF correlation function we perform a local randomization of the points along a line with the adjacent points (simply we exchange the order, but randomly). If this procedure is not applied, some extra correlations (along the lines) are visible only if the intensity is very low. Since the routine is applied if selected, the user can choose not to apply the routine and observe the differences. The ability to use sCMOS cameras for fluctuation spectroscopy is a crucial development for the entire community since these cameras are fast and relatively inexpensive.
One important issue for the application of the 2D-pCF algorithm is that the system must be stationary for the duration of the measurement. For the measurement reported in this manuscript the duration was about 100s. If motion of the cell or part of the cell occurs during the observation time, we could have measurement artifacts. There are general procedures to detect motion artifacts, such as fiduciary markers or motion correction algorithms. Although these methods were not applied to the data shown in this work, we have implemented a "moving average" method in our software based on dividing the total duration in shorter segments and comparing the data for each of the shorter segments. This allows us to detect motion artifact and possibly excluding the segments where the artifacts are present. Based on simulated data, the effect of motion is not detectable if the total drift is less than 1-2 pixels (typically a pixel is about 100 nm) during a data segment. However, if the motion is about 10 pixels or more then the motion appears as an additional eccentricity in the direction of the motion. In this case we disregard the segment or the entire measurement.
A final comment is about comparison with other methods that have been used to map the molecular connectivity in a cell, notable methods based on FLIP-photobleaching [43]. A fundamental difference between fluctuation and photobleaching methods is that fluctuations are measured at equilibrium and on minimally perturbed systems. Fluctuations methods can be used continuously and over an extend period of time in the SPIM microscope, allowing to follow stimulation of events until the stimulus has finished. Fluctuation methods are intrinsically single molecule approaches and in this respect can be compared with single particle tracking. Photobleaching is an ensemble approach and, although very powerful and used by many, it provides average information. This average information could be all you need but we are pointing out here the difference between the 2D-pCF and the bleaching based approaches.
MA-USA) was used. A 10 ms exposure time was used for the 2D-pCF data acquisition. A detailed description can be obtained from [41].
Measurement on widefield instrumentation were done in a commercial Olympus TIRF microscope (Olympus, USA), using mercury lamp illumination and two filter set, for EGFP a HQ500/20 and 515 nm long pass filter split by a Q515lp dichroic mirror (Chroma Technology Corporation., Bellows Falls, VT-USA), for DiOC16 a Texas red D560/40 nm and 610 nm long pass filter split by a 600 nm long pass dichroic mirror (Chroma Technology Corporation., Bellows Falls, VT-USA). The signal was collected using an Olympus TIRF UIS2 PlanApo N 60x (NA = 1.45) oil-immersion objective and fluorescence fluctuation was recorded by an EMCCD camera 512B Cascade (Photometrics, Tucson, AZ, USA). The frame rate was 10 ms and the image size 128x128 pixels.

Optimization and testing of the 2D-pCF algorithm using a simulated experiment
The algorithm to calculate the pCF(δr) (Eq. (1)) as a function of the distance, δr, for each point of an image has been parallelized in order to reduce the computational time and the memory requirements. The algorithm constructs what we call a "sprite' at each point of the image. A "sprite' is a small image of 32x32 points. The first index is the angle and the second index is the delay time in log scale. The sprite image is what is shown in the polar images in Fig. 2 (panels D to I). Given a point, the algorithm calculates a sprite for each distance δr. So for an image of 256x256 pixels there are about 2 16 sprites, and the calculation is repeated for 4 to 8 distances δr, for example, for each δr we produce a sprite file of about 2 26 or 256 MB (4 bytes per entry). The sprites, calculated only once, are compressed and stored in the relatively small sprite file. The visualization of the results of the 2D-pCF calculation is done by computing the central moments of each sprite polar image. There are a total of 5 parameters from the central moment's calculation (Center-x, Center-y, ellipses long axis, ellipse short axis and ellipses direction) for each sprite (Eqs (2)- (7)). These 5 primitive maps are used to construct all the derived maps including diffusion eccentricity and anisotropy, time delay long, time delay short, angle, velocity modulus and angle of velocity, and the connectivity map. The calculation of the maps from the sprite file is virtually instantaneous because there is no fit, but only calculation of the second order moments for each sprite map. Also the 5 primitive maps are compressed and stored with a typical size of less than 1 MB total. The highly parallelized calculation of the sprites for a 256x256 image stack of 32,768 frames (about 4 GB of data) takes less than 60 s in a laptop and produces a reduced representation of the diffusion eccentricity and connectivity in files of less than 1 MB. The visualization of connectivity uses the eccentricity and the time delay maps and it does not require recalculating the sprite file again, since it is a simple visualization of a combination of the maps. For the connectivity map we first establish a threshold value for the eccentricity, for example 0.45. Then at each pixel we plot a segment that has the orientation of the diffusion anisotropy (the ellipses direction) at that pixel and a length proportional to the long axis of the ellipses as shown in Fig. 2. The algorithms for the calculation of the sprites can be found in the program SimFCS available at www.lfd.uci.edu.