![[Feedback]](/icons/banner/feedbackACT.gif){kind=icon}
![[For Subscribers]](/icons/banner/subscriberACT.gif){kind=icon}
![[Archive]](/icons/banner/archiveACT.gif){kind=icon}
![[Search]](/icons/banner/searchACT.gif){kind=icon}
![[Contents]](/icons/banner/tocACT.gif){kind=icon}
|
|
|
|
|
||
Vol. 15, Issue 10, 4749-4760, October 2004
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Institute of Molecular Medicine, Faculty of Medicine, University of Lisbon, 1649-028 Lisbon, Portugal
Submitted June 17, 2004;
Accepted July 23, 2004
Monitoring Editor: Joseph Gall
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
30% in aqueous solution but has little effect on mobility in the nucleoplasm. This suggests that spatial constraints to diffusion of dextrans inside the nucleus are insensitive to temperature. |
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
) as a tool to study the molecular mobility in several media, including aqueous solutions, gels, and living cells. In cell biology, FRAP was originally used to study membrane diffusion of lipids and proteins coupled to fluorophores (Axelrod et al., 1976; Edidin et al., 1976). More recently, noninvasive fluorescent tagging with the green fluorescent protein (GFP) or its variants has stimulated the use of FRAP to measure protein dynamics inside the living cell (Phair and Misteli, 2000; Reits and Neefjes, 2001).
In a typical FRAP experiment a region of the cell is briefly illuminated with a high-intensity laser beam with a wavelength near the excitation peak of a fluorophore. Consequently, most of the fluorophores inside that region loose their fluorescence irreversibly, a phenomenon known as photobleaching. As fluorophores are attached to molecules that move in and out of the bleached region, fluorescence inside that zone increases, and eventually equilibrium is reached. A higher mobility of the molecules results in a shorter time of recovery. From the recovery curve, it is possible to obtain estimates of the diffusion coefficient and immobile fraction (i.e., the percentage of proteins that are unable to move) (for reviews, see White and Stelzer, 1999; Klonis et al., 2002; Carrero et al., 2003).
Early applications of FRAP involved the use of nonscanning microscopes (Axelrod et al., 1976; Soumpasis, 1983; Tsay and Jacobson, 1991; Berk et al., 1993). In most cases, bleaching was performed with a stationary laser beam (either considered to be uniform or Gaussian) focused toward thin samples where diffusion could be considered to be mostly two-dimensional. Fluorescence recovery was then recorded by the fluorescence microscope. To obtain estimates for the diffusion coefficients and immobile fractions, appropriate theoretical models were developed. One of the most widely used models in cell biology was developed by Axelrod et al., 1976 (Calapez et al., 2002; Wei et al., 2003; Shimi et al., 2004). The model assumed an infinite homogenous medium with cylindrical symmetry. Bleaching was considered to be a first-order linear process, taking a negligible amount of time compared with diffusion. It also was assumed that there was no significant delay between bleaching and scanning. The laser was considered to have a Gaussian intensity profile, for both bleaching and observation of recovery. Using these assumptions, the two-dimensional diffusion equation was solved and a formula for fluorescence evolution was found.
In conjunction with GFP tagging of intracellular proteins, the development of confocal laser scanning microscopes (CLSMs) provided cell biologists with an excellent standard tool to perform FRAP experiments. Currently, most commercially available CLSMs are equipped with acoustic-optic tunable filters (AOTF) that allow bleaching arbitrary regions in the sample. During bleaching, the microscope directs a moving excitation laser beam to the sample. The bleach region is thus scanned point by point and line by line, increasing the incident laser intensity inside the pattern defined (for a thorough theoretical explanation of the bleaching process in a CLSM, see Braeckmans et al., 2003).
In the context of CLSM imaging, some assumptions of the models developed for nonscanning microscopes do not hold, because bleaching and monitoring of fluorescence recovery are done in a fundamentally different way.
Despite the widespread use of photobleaching microscopy, quantitative analysis of the data has been limited by lack of appropriate FRAP models for the CLSM. Although two-dimensional (2D) and three-dimensional (3D) models have been developed based on numerical approaches (Wedekind et al., 1996; Kubitscheck et al., 1994, 1998; Cutts et al., 1995; Peters and Kubitscheck, 1999), these have had little practical application in cell biology. A more practical 3D model was devised for use with objective lenses of low numerical aperture (Blonk et al., 1993). However, this method relies on a stationary laser beam for bleaching and a line-scanning beam for recording the fluorescence recovery, a sequence that is not possible in standard commercial CLSMs (Braeckmans et al., 2003). Recently, a new 3D FRAP model that can be easily applied on almost any modern CLSM was reported (Braeckmans et al., 2003). This method was shown to be valid only when a number of assumptions are met by the experimental conditions. Namely, that an objective lens of low numerical aperture is used, and that the bleaching phase is sufficiently short to avoid recovery during bleach.
However, in most cell biological applications of FRAP microscopy, a region of at least 1 µm2 is bleached for hundreds of milliseconds and acquisition of an image takes at least tens of milliseconds (i.e., the bleached region of interest [ROI], which is typically near the center of the image, is not imaged immediately after the end of bleaching). Thus, it is expected that for molecules diffusing at rates as low as 0.2 µm2 s-1, fluorescence recovery will occur before the bleached zone is effectively imaged for the first time. Here, we demonstrate experimentally that this is indeed the case. We therefore devised a novel 2D and 3D approximation that takes into account diffusion during the bleach period. The new FRAP model can be easily programmed in a standard fitting routine and readily applied on most commercial CLSMs, even when high numerical aperture objectives lenses are used.
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
Enhanced green fluorescent protein plasmid (pEGFP) (BD Biosciences Clontech, Palo Alto, CA) for transfection assays was purified using plasmid DNA midi-prep kit (QIAGEN, Hilden, Germany). HeLa subconfluent cells were transiently transfected with 1 µg of purified pEGFP plasmid by using FuGENE6 reagent (Roche Biochemicals, Indianapolis, IN) according to the manufacturer's protocol.
For imaging and microinjection, the medium was changed to Dulbecco's-MEM/F-12 without phenol red supplemented with 15 mM HEPES buffer (Invitrogen).
Dextran Purification
Fluorescein isothiocyanate (FITC)-labeled dextrans with average molecular sizes of 20, 40, and 70 kDa (Sigma-Aldrich, St. Louis, MO) were dissolved to 10 mg/ml in water. One milliliter of each sample was fractionated on a Superdex 200 column (Amersham Biosciences, Piscataway, NJ). For each dextran, the column fractions with higher absorbency were pooled and lyophilized. The 500-kDa FITC-labeled dextran was used without further purification. For imaging, all samples were diluted to 200 µg/ml in water.
Confocal Microscopy and Image Analysis
Live cells and FITC-dextrans in aqueous solution were imaged at either 37 or 22°C maintained by a heating/cooling frame (LaCon, Staig, Germany), in conjunction with an objective heater (PeCon, Erbach, Germany). Images were acquired on an LSM 510 (Carl Zeiss, Jena, Germany) with the Planapochromat 63x/1.4 objective. FITC fluorescence was detected using the 488 laser line of an Ar laser (25% of 25-mW nominal output) in conjunction with a LP 505 filter. For imaging the 3D bleached volume, 500-kDa FITC-dextrans were immobilized in a 15% polyacrylamide gel, and a z-stack series was acquired. In total, 39 images were taken and the distance between each image in the stack was 0.4 µm. Image size was 256 x 256 pixels, and the pixel width was 114 nm. Background was estimated from a 15% polyacrylamide gel devoid of fluorescent molecules.
Each FRAP experiment started with three image scans, followed by a bleach pulse of 242 ms on a spot with a diameter of 21 pixel (1.19-µm radius). A series of 97 single section images were then collected at 29.82-ms intervals, and the first image started to be acquired 2 ms after the end of bleaching. Image size was 256 x 30 pixels and the pixel width was 114 nm. For imaging the laser power was attenuated to 2% of the bleach intensity.
For FRAP analysis of intracellular mobility, FITC-dextrans (at 200 µg/ml in water) were microinjected into the nucleus of HeLa cells, as described previously (Almeida et al., 1998). In each time series, the background and nuclear regions were identified using an implementation of the ICM segmentation algorithm (Calapez et al., 2002) in Matlab software (Mathworks, Natick, MA). The average fluorescence in the nucleus T(t) and the average fluorescence in the bleached region I(t) were calculated for each background subtracted image at time t after bleaching. FRAP recovery curves were normalized according to Phair and Misteli (2000))
 | (1) |
Image processing routines also outputted the normalized fluorescence profile of the first postbleach image. From the fit with Eq. 5, parameters wM and KM were obtained.
FRAP recovery curves were then fitted with Eq. 10, yielding parameters 
D, 
, and K. All fitting procedures were performed with NonLinearRegress function of Mathematica 4.0 (Wolfram Research, Champaign, IL).
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
)
 | (2) |
In most commercially available CLSMs, the bleaching laser beam is not stationary (Braeckmans et al., 2003). Rather, these microscopes are equipped with a scanning beam. The user defines the bleach region of interest (ROI), and the laser scans this region line by line. Laser intensity is precisely controlled by an AOTF and is set to a high value only when the laser beam enters the region defined (Figure 1).
|
To directly observe the fluorescence profiles generated by bleaching a region of interest with a scanning laser beam, we immobilized FITC-dextrans in a 15% polyacrylamide gel and imaged the bleached volume along the z-direction. For bleaching, we used a small circular region with radius wB (1 µm).
As depicted in Figure 2, the observed fluorescence profiles can be approximated by the following equation:
 | (3) |
|
Noteworthy, Eq. 3 has the same functional form as Eq. 2, indicating that bleaching profiles generated in a standard CLSM approach those produced by an idealized stationary Gaussian laser beam.
Fitting with Eq. 3 was then tested for different radial values of the bleach region (0.14, 0.26, 0.71, 1.14, and 2.28 µm). Fitting with Eq. 3 was valid for the bleaching spots smaller or equal to 1.14 µm. The values of Kim were contained in a typical range from 0 to 6. For the 2.28-µm bleach region, the fit yielded unexpectedly high values for Kim (60). In this region, the postbleach profile is constant for r <2.0 µm, which is consistent with the uniform disk model (Braeckmans et al., 2003). The maximum size of the bleaching region radius, for which fitting Eq. 3 is valid, will depend essentially on the effective size of the bleaching beam spot, which is related with the numerical aperture of the objective used, and should be experimentally determined for each objective. Because the size of the bleaching laser beam is roughly proportional to 1/numerical aperture (NA), lower NA objective lenses are expected to allow the use of larger bleaching regions.
According to Eq. 3, normalized fluorescence inside a region of radius ws then yields (see Appendix A)
 | (4) |
; Braeckmans et al., 2003).). It also is assumed that all image points are taken at the same time, but in a CLSM pixels are acquired sequentially, so the bleached region is imaged a few milliseconds after the start of image acquisition. For simplicity, the time taken from the start of the bleaching phase until the first effective scan of the ROI will be hereafter referred as total bleach time (Figure 1C). To test whether if diffusion during this period can be neglected in a typical cell biological application of FRAP, we imaged the radial fluorescence profiles of FITC-dextrans immediately after bleaching (Figure 3). Dextrans of 40 and 500 kDa in aqueous solution were used. Assuming that no significant diffusion occurs during total bleach time, fluorescence profiles should be similar for both dextrans. As shown in Figure 3, this is clearly not the case. The 40-kDa dextran molecules diffuse faster than those of 500 kDa and by the time the postbleach profile is acquired by the CLSM, a larger fraction of the 40-kDa bleached molecules have moved away from the original bleach region. Consequently, the fluorescence reduction is less pronounced and the radius of the bleached area is larger.
|
For mobile molecules, we found that with the following approximation good fitting results were obtained for the postbleach radial profile:
 | (5) |
Assuming isotropic diffusion, we expect that the best 3D fitting results are obtained with the following approximation:
 | (6) |
0 we have wM wim, zM zim and KM Kim, and thus this approximation approaches Eq. 3.
Derivation of the Model: 3D Recovery after Bleaching
As in most previous models (Tsay and Jacobson, 1991; Blonk et al., 1993; Braeckmans et al., 2003), we consider diffusion as the only relevant transport process and that the medium is homogenous. Thus, after bleaching, the system behaves according to the classical diffusion equation
 | (7) |
is the concentration of the unbleached fluorophores at position 
at time t.
The partial differential equation above may be written in terms of observed fluorescence and with Eq. 6 as initial value, it has a known analytical solution (Blonk etal., 1993). Normalized fluorescence recovery inside a region of radius wS (which in general can be different from wB) is thus given by (see Appendix B):
 | (8) |
 |
The square root term in Eq. 8 takes the finite axial size of the bleached volume into account. If the sample thickness is smaller than that, the system will behave like in a 2D situation. The 2D counterpart of Eq. 8 is found by zim+
 | (9) |
So far, we assumed that all the molecules are diffusing. But if a fraction of molecules are immobile then the normalized recovery curve will be
 | (10) |
im is the normalized fluorescence of the immobile molecules as determined in the previous section (Eq. 4), and M is the normalized fluorescence of the mobile fraction (from either Eqs. 8 or 9). The variant using the 3D correction (Eq. 8) will be referred as 3D variant, whereas its 2D counterpart (Eq. 9) will be referred as 2D variant. This equation computes the fluorescence recovery curve for molecules in which some of them are diffusing while the others are immobile. The first term of Eq. 10 is the contribution of the fluorescence of the mobile fraction to the normalized fluorescence recovery, whereas the second term is a theoretical estimate of the postbleach total fluorescence of the immobile fraction inside the ROI (which is assumed to remain constant during a FRAP experiment).
It should be emphasized that wM refers to the profile width of the mobile population. However, in the presence of an unknown immobile fraction, it is not possible to directly estimate wM because the mobile and immobile populations are indistinguishable in an image. The postbleach concentration profile width of a mixture is always narrower than the purely mobile one (i.e., with = 0). Thus, an underestimation of the diffusion coefficient is expected especially in the presence of a high immobile fraction.
To overcome this problem we propose to make the estimation in two steps. In the first step, the recovery curve is fitted with Eq. 10, by using the fluorescence profile width, as shown by the first postbleach image, as a first estimate of wM. Numerical simulations show that immobile fraction estimates are relatively independent of the real value of wM (our unpublished data). Thus the immobile fraction can be estimated with a reasonable degree of accuracy. This immobile fraction value would be used to make an estimation of the real mobile fraction concentration, according to the relation
 | (11) |
Validation of the Model by using Numerical Simulations
Considering diffusion as the only relevant transport process and assuming homogeneity of the medium, the evolution of the concentration of unbleached fluorophore is described by the following partial differential equation accounting for diffusion during bleach:
 | (12) |
is the 3D laser intensity distribution at time t of the bleaching beam, and µ the fluorophore bleach rate constant. To check the influence of the intensity distribution on parameter estimation several configurations were used: stationary Gaussian beam, scanning Gaussian beam, modified Gaussian beam (Kubitscheck et al., 1998) and a "stationary" intensity distribution derived from experimental data obtained with immobile molecules. The modified Gaussian beam approximates the true illumination profile by
 | (13) |
 | (13a) |
being the wavelength of the laser and w0 its radial width at z = 0.
The characteristic recovery time is
 | (13b) |
). In simulations employing the "stationary" intensity distribution, the total bleaching time chosen was 240 ms. For slowly diffusing molecules (<0.1 µm2 s-1), this time is insignificant in comparison with the characteristic recovery time. In this situation, the stationary intensity distribution derived from the experimental data provides a more realistic bleaching profile.
After bleaching, the system behaves according to the classical diffusion equation
 | (14) |
 | (15) |
1 µm) is considered to be much smaller than the size of the compartment under study (i.e., 6 µm of radius in a HeLa cell nucleus). First, we used simulations to assess how the scanning laser beam shape influenced the postbleach profiles (for simulation details, see Appendix C). Results indicate that the bleach profile generated by a scanning Gaussian laser is still properly approximated by Eq. 6 if the size of the bleach ROI is smaller or comparable to the size of the Gaussian bleaching laser (Figure 4, A and B). A scanning beam deviating strongly from a Gaussian will generate accordingly a non-Gaussian bleach pattern (Figure 4, C and D). Together with the experimental measurements of the bleached volume, simulations further support the idea that, for our purposes, the bleaching laser beam can be considered to be approximately Gaussian, even when using a high NA lens.
|
Figure 5 summarizes simulation results, comparing the performance of the two 3D variants of the new method for several configurations. For the stationary Gaussian configuration, in A, for the diffusion coefficients range tested (between 0.1 and 24 µm2 s-1), the maximum error in the estimates is 12% (for D = 10 µm2 s-1), decreasing for higher diffusion coefficients. For high diffusion coefficients (D > 5 µm2 s-1), the 3D variant has a better performance than the 3D variant with IF. In B, the influence of the immobile fraction is tested. It is shown that the estimates error grows continuously for the 3D variant, whereas for the 3D variant with IF the error is <10% until = 70%. In C, the influence of the 3D shape of the illuminating beam in scanning configurations is tested. If the beam deviates strongly from the Gaussian approximation, the error on the estimated D is much larger than in the Gaussian case. The error for the Gaussian beam is relatively small, as the bleach pattern generated by this beam could be effectively approximated by an exponential of a Gaussian function. This situation is similar to what was obtained experimentally, so we expect that when the new method is applied to a real situation its estimation performance is also similar.
|
We also tested the method with the experimentally obtained bleach profile, for diffusion coefficients ranging from 0.01 µm2 s-1 to 0.1 µm2 s-1. Both 3D methods yielded similar estimates and diffusion coefficients were slightly overestimated by 5–9%.
Simulations were performed with a bleach ROI radius much larger (5x) than the radial extent of the bleaching laser. For the example depicted in Figure 4, G and H, the error on the diffusion coefficients estimate was 6%. One reason for this low error is that the fluorescence profiles, due to diffusion, rapidly approach a curve similar to Eq. 6, approaching the 3D method assumptions. For a bleach time small compared with characteristic diffusion time (i.e., for D <0.50 µm2 s-1, in our case), the bleaching profile obtained (Figure 4, E and F) is consistent with the predictions from the uniform disk model (Braeckmans et al., 2003). This should be the method of choice in these conditions.
Our previous methodology based on a model proposed by Axelrod et al., 1976, underestimated significantly diffusion coefficients, especially when simulating highly mobile molecules. This method yielded the best results for simulations performed under its assumptions (our unpublished data). This also shows that both numerical simulations and fitting routines were working properly.
Experimental Validation of the method
Assuming that dextrans in aqueous solution have a dynamical behavior similar to random walking particles, we estimate their diffusion coefficient (D) using the Stokes-Einstein equation:
 | (16) |
the viscosity of the medium, and R the hydrodynamic radius of the molecule. At 37°C (310 K), water has a viscosity of 0.6915 mPa.s and at 22°C (295 K) its viscosity is 0.955 mPa.s. Because dextran molecules are not perfect spheres, the theoretical values must be considered as an approximation of the expected diffusion rates.
First, we used a method based on Axelrod formula (Axelrod et al., 1976) and that has been widely used to perform quantitative FRAP in cell biology (Phair and Misteli, 2000, Calapez et al., 2002; Wei et al., 2003, Shimi et al., 2004), to determine the diffusion coefficient of a FITC-labeled 20-kDa dextran in aqueous solution. The estimated value was 21 µm2 s-1, contrasting with an expected value of 97 µm2 s-1 (as calculated from the Stokes-Einstein equation). Clearly, this method is underestimating the diffusion coefficient of fast-moving molecules. Next, we estimated the diffusion coefficient values for each dextran by using Eq. 10 (3D variant), as sample thickness is of the order of 100 µm, larger than the axial extent of the bleached volume (9 µm). Nearly complete fluorescence recovery for all the FITC-dextrans was observed (>95% of the initial fluorescence), which is consistent with previous reports (Seksek et al., 1997). Table 1 shows, for dextrans of different sizes dissolved in water, the theoretically expected values of diffusion coefficients and the experimentally obtained values. As indicated by the ratio between expected and estimated values, the diffusion coefficients measured by the new FRAP model are always very close to the theoretical predictions (Figure 6).
|
|
Having experimentally validated the model, the different dextrans were then microinjected into the nucleus of HeLa cells. As the average thickness of a HeLa cell in the z-direction (5 µm) is almost 2 times smaller than the axial extent of the bleached volume, the 2D variant of the model was selected (Eq. 9). The diffusion coefficient values and IF values estimated in cells maintained at either 37 or 22°C are shown in Table 2 (also see Figure 7). As can be seen from the ratio between values estimated in water and in cell nucleus at 22 or 37°C, the mobility impairment inside the nucleus is higher for the higher molecular weight dextrans.
|
|
Finally, using the 2D variant of the new method we estimated a diffusion coefficient of 33.3 ± 3.6 µm2 s-1 for GFP in the HeLa cell nucleus, with an IF of 1.0 ± 0.1%. This is consistent with recently reported values of 22 ± 2 µm2 s-1 at 22°C (with an of IF = 4.0 ± 0.3%) (Potma et al., 2001) and 32 ± 4 µm2 s-1 at 37°C (Beaudouin, 2003).
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
The GFP is a 28-kDa protein with a Stokes' radius of 2.35 nm (Ribbeck and Görlich, 2001). This corresponds to an expected diffusion coefficient in water at 37°C of 139.2 µm2 s-1. Confocal spot photobleaching recovery measurements performed with a specifically modified microscopic system yielded GFP diffusion coefficients of 87 ± 3 µm2 s-1 in aqueous solution and 24 ± 2 µm2 s-1 in the cytoplasm of Dictyostelium discoideum cells kept at room temperature (Potma et al., 2001). In our previous work, we used a standard CLSM to perform quantitative FRAP, and we estimated an apparent diffusion coefficient of 5 µm2 s-1 for GFP in the nucleus of HeLa cells (Calapez et al., 2002). We believe that this value was significantly underestimated because the mathematical model applied considered that fluorescence recovery during the bleach phase is negligible, an assumption that is not met by the experimental conditions in the case of highly mobile molecules (Figure 3). Using the same CLSM instrument and the new FRAP model, we estimate a diffusion coefficient of 33.3 ± 3.6 µm2 s-1 for GFP in the nucleus of HeLa cells maintained at 37°C.
To develop a more accurate practical approach that can be readily applied by cell biologists interested in performing quantitative FRAP analysis with a standard CLSM, we first analyzed experimentally the fluorescence profile generated by bleaching a region of interest with a scanning laser beam. FITC-labeled dextrans were immobilized in a polyacrylamide gel, and the bleached volume was imaged along the z-direction. The results show that, for bleaching spots up to 1 µm in radius, the observed fluorescence profiles can be effectively approximated by the theoretical models that consider bleaching generated by a stationary Gaussian laser. Our data thus validate the application of the Gaussian approximation to postbleach radial fluorescence profiles generated with the CLSM, as assumed previously (Calapez et al., 2002). However, our results further indicate that in order to apply the Gaussian approximation the bleach region should not exceed 1 µm in radius (the maximum radius for which fitting is valid depends on the numerical aperture of the objective lens). For much larger bleach regions, a better approximation is the uniform disk profile, as predicted from the model proposed by Braeckmans et al., 2003.
To date, most FRAP models assume that both bleaching and image acquisition are sufficiently fast to avoid diffusion during those periods (Braeckmans et al., 2003). As a rule, for diffusion to be neglected, the total bleaching time should be at least 15 times smaller than the recovery time (Meyvis et al., 1999). If we consider a molecule that diffuses at 0.5 µm2 s-1 and a bleaching region with a radius of 1 µm, then the characteristic recovery time will be 0.5 s. The bleaching time should therefore be <33 ms. In practice, most cell biological FRAP applications using a similar bleach ROI size use bleach periods of 100 ms or longer (Phair and Misteli, 2000; Shimi et al., 2004). One may therefore predict that in most FRAP experiments the assumption that negligible fluorescence recovery occurs during bleaching is not correct.
A direct demonstration that before the first postbleach image is acquired significant diffusion takes place (namely, for a 40-kDa dextran) is depicted in Figure 3. Diffusion during bleach is obviously less important for slower moving molecules (i.e., a 500-kDa dextran).
Using the approximations mentioned above. it was possible to derive analytical formulas for the normalized fluorescence recovery curve. The three-dimensional correction is important for measurements in bulk solutions. In cellular measurements, it will contribute significantly only if the cell thickness is larger than the axial extent of the bleach volume. In a thin sample, the 2D variant should be used instead.
For simultaneously high values of the immobile fraction and of the diffusion coefficient the capability of the 3D (or 2D) method to correctly estimate these parameters is reduced. This happens because the mobile fraction concentration profile is not directly identified, being wM estimated from the fluorescence profile of the mixture of the mobile and immobile fractions. Thus, a higher level of immobilized molecules leads to an underestimation of the mobile fraction profile radius, and consequently, to an underestimation of the diffusion coefficient. A way to circumvent this problem was devised, by making a fitting procedure (with Eq. 10) in two steps. The first run would be used to estimate the immobile fraction only, computing then the mobile fraction fluorescence profile (Eq. 11), and in the second run corrected profile width values would be used. It is important to highlight that the approach describe here requires a new calibration every time there is a change in bleaching disk radius, laser beam power, or fluorophore.
To be sure that these expressions were adequate for analysis of real FRAP experiments, we generated recovery curves from the simulations and fitted those curves with the formula. The estimated diffusion coefficients and immobile fractions were always close to the parameters used in the simulations (Figure 5).
The effect of noise was tested by generating several curves with the same level of Gaussian noise added to a simulated recovery curve. This procedure led to the conclusion that noise has a relevant role in the quality of the estimates, introducing some variability in the estimation. This was especially noticeable when simulating molecules with high diffusion coefficients. The best way to improve the estimation was to average the highest number possible of experiments.
As a first biological application, the new FRAP method was used to compare the diffusion rates of different size macromolecules in aqueous solution and in the nucleus of living HeLa cells. FITC-dextrans were either directly imaged in solution or microinjected into the nucleus. A rather homogeneous fluorophore distribution was observed, suggesting that the dextrans spread freely throughout both the aqueous sample and the nucleoplasm. In aqueous solution, the ratios between the diffusion coefficient values theoretically expected from the Stokes-Einstein equation and those estimated by FRAP both at 37 and 22°C were close to 1. As expected, in the nucleoplasm the diffusion coefficient values decreased relative to water. Noteworthy, higher molecular weight dextrans were proportionally more retarded in the nucleoplasm, at 37 and 22°C, than smaller molecular weight molecules (Table 2). This means that nucleoplasm deviates from an ideal liquid behavior, with the effective viscosity increasing with the size of the molecules. A similar observation was reported for FITC-dextrans injected into the eye vitreous, but not for dextrans diffusing in cystic fibrosis sputum (Braeckmans et al., 2003), where the decrease in diffusion coefficient values of the FITC-dextrans seemed rather independent of their size. In contrast with the aqueous sputum, the vitreous is composed of a meshwork formed by polymers of hyaluronic acid. Most probably, this meshwork causes a sterical hindrance that is stronger for larger molecules (Braeckmans et al., 2003). A similar situation is likely to occur in the nucleoplasm, where sterical hindrance is caused by the meshwork composed of chromatin and nonchromatin nuclear components. Most interestingly, the sterical hindrance effect inside the nucleus is similar at 37 and at 22°C (1.7± 0.3 and 1.5 ± 0.3 µm2 s-1 for 500-kDa dextrans). Yet, decreasing the temperature from 37 to 22°C reduces the expected mobility rate of the dextrans in aqueous solution by 30% (from 32.8 ± 1.9 to 23.2 ± 1.1 µm2 s-1 for 500-kDa dextrans). We have previously observed that when cells are depleted of ATP or incubated at 22°C, messenger ribonucleoprotein (mRNP) particles show significantly reduced mobility rates in the nucleus, whereas large-molecular-weight dextrans are not much affected (Calapez et al., 2002). The results reported here reinforce the view that spatial constraints to diffusion of dextrans inside the nucleus are insensitive to temperature and therefore energy independent. In contrast, energy-dependent processes are possibly involved in facilitating the diffusion of mRNP complexes in transit to the cytoplasm (Calapez et al., 2002; Carmo-Fonseca et al., 2002).
Inside the cell, binding events are expected to slow down the diffusion dynamics of macromolecules. In the case that the binding reaction is faster than the typical times involved in the diffusive process, an effective diffusion coefficient can be defined (Crank, 1975). Such effective diffusion coefficient integrates the absolute diffusion coefficient and the rates of association and dissociation of the macromolecule. However, further work is needed to address the complex interplay between macromolecular diffusion and binding events that takes place in the living cell (Phair et al., 2004; Sprague et al., 2004).
|
Appendix |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
; Braeckmans et al., 2003). Accordingly, we averaged our fluorescence signal over a circular region of radius wS from the images taken at the focal plane (z = 0). wS in general can be different from wB. The normalized fluorescence is thus
 | (A.17) |
For immobile molecules, if bleaching is performed on a small bleaching region (i.e., when fitting with Eq. 3 is valid) we have
 | (A.17a) |
B. Mobile Fraction: Fluorescence Recovery Inside a Disk of Radius wS
In any fluorescence microscope, an image created by an object is the convolution of the fluorophore concentration and the point spread function (PSF) of the system (Berk et al., 1993). If the system behaves according to the classical diffusion equation, we have from the convolution properties
 | (A.17b) |
 | (A.18) |
 | (A.19) |
 | (A.20) |
), yielding
 | (A.21) |
 | (A.22) |
 | (A.22a) |
After bleaching, by the time the bleached region is imaged the expected theoretical profile width should be approximately (similarly to Jain et al., 1990):
 | (A.23) |
T is the difference between total bleach time and the bleach phase (Fig. 1). As expected, this formula states that the speed by which profiles spread is higher for faster diffusive rates. To correct for slight mismatches between the real width and the estimated width in the axial direction, we multiply Eq.A.7 by the ratio between real (wM) and estimated values in the radial direction. So zM becomes:
 | (A.24) |
Substituting this in Eq.A.22 yields Eq. 8 in the text.
C. Numerical Simulations
To obtain simulated recovery curves, Eqs. 12 and 14 were solved numerically in the three spatial dimensions (x,y,z) plus time, by a finite-differencing scheme implemented in Mathematica 4.0 (Wolfram Research). The user controls fluorophore and bleaching beam parameters (the beam width in the axial and radial direction, the laser intensity, the size of the bleach ROI, the number of lines, the time taken per line, the start and the end points of the scanning lines, and the bleach rate constant), the time of bleach phase, the dynamical properties of the molecule (D, ), acquisition parameters (time between images, the number of bleach iterations, the start of imaging phase and the total duration of the simulated FRAP experiment) and numerical parameters (the grid points number in each spatiotemporal direction, and the location of the boundaries). The bleach region was located at the center of the volume, and the boundaries were set to a distance at least two times larger than the bleaching laser width in that direction. The grid points number was selected to obtain stable numerical solutions.
Simulations parameters chosen are close to those of a real FRAP experiment (see Materials and Methods). For details, see below.
The time between images was set to 30 ms (an image takes 28 ms and the next image starts to be scanned after 2 ms), and there is a time delay between the end of bleach and the start of imaging of 2 ms. The bleached zone is, on average, in the center of the imaging area, so the first postbleach image of the bleached area starts 28/2 + 2 ms = 16 ms after the end of bleaching.
In the scanning laser configurations, the time taken for each line is 1 ms. The bleach region is a circle of 1.19-µm radius (containing 20 scanning lines). This region is scanned for two iterations, for a total bleaching phase time of 40 ms. Although in practice we use longer bleach times, this is still longer than the characteristic recovery time of fast diffusing molecules. The choice of smaller times in these simulations is due to computer memory constraints.
Laser parameters were chosen to get a bleach immobile profile similar (in terms of spatial extension) to the experimentally obtained for dextrans immobilized in a polyacrylamide gel.
Usually, the set of pixels from which fluorescence recovery is computed is the same as the bleach ROI defined. According to the new method its size is independent of bleach ROI defined, changing its size had no significant influence on the results (our unpublished data). If the radius of the monitoring zone is increased the recovery curve becomes sigmoidal, which was also predicted from the analytical expression.
Gaussian noise was added to the simulated recovery curves and fluorescence profiles. The fluorescence profiles were fitted with either Eq. 3 (for immobile molecules) or Eq. 6 (for mobile molecules). The recovery curves were fitted with Eq. 10 by using NonlinearRegress, yielding estimates for KM, D, and .
|
ACKNOWLEDGMENTS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
|
Footnotes |
|---|
* Corresponding author. E-mail address: josebraga{at}fm.ul.pt.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND METHODSRESULTSDISCUSSIONAppendixACKNOWLEDGMENTSREFERENCES |
|---|
Axelrod, D., Koppel, D.E., Schlessinger, J., Elson, E., and Webb, W.W. ((1976). ). Mobility measurement by analysis of fluorescence photobleaching recovery kinetics. Biophys J. 16, , 1055-1069.[Medline]
Beaudouin, J. ((2003). ). Structural and molecular dynamics of nuclear proteins revealed by fluorescence microscopy. Ph.D. dissertation. Available at http://library.embl-heidelberg.de/pdf/beaudouin.pdf.
Berk, D.A., Yuan, F., Leunig, M., and Jain, R.K. ((1993). ). Fluorescence photobleaching with spatial Fourier analysis: measurement of diffusion of lightscattering media. Biophys J. 65, , 2428-2436.[Medline]
Blonk, J.C.G., Don, A., Van Aalst, H., and Birmingham, J.J. ((1993). ). Fluorescence photobleaching recovery in the confocal scanning light microscope. J. Microsc. 169, , 363-374.
Braeckmans, K., Peeters, L., Sanders, N.N., De Smedt, S.C., and Demeester, J. ((2003). ). Three-dimensional fluorescence recovery after photobleaching with the confocal scanning laser microscope. Biophys J. 85, , 2240-2252.[Medline]
Calapez, A., Pereira, H.M., Calado, A., Braga, J., Rino, J., Carvalho, C., Tavanez, J.P., Wahle, E., Rosa, A.C., and Carmo-Fonseca, M. ((2002). ). The intranuclear mobility of messenger RNA binding proteins is ATP dependent and temperature sensitive. J. Cell Biol. 159, , 795-805.
Carmo-Fonseca, M., Platani, M., and Swedlow, J.R. ((2002). ). Macromolecular mobility inside the cell nucleus. Trends Cell Biol. 12, , 491-495.[CrossRef][Medline]
Carrero, G., McDonald, D., Crawford, E., de Vries, G., and Hendze, M.J. ((2003). ). Using FRAP and mathematical modeling to determine the in vivo kinetics of nuclear proteins. Methods 29, , 14-28.[CrossRef][Medline]
Crank, J. ((1975). ). The Mathematics of Diffusion, 2nd ed. Clarendon Press, Oxford, United Kingdom.
Cutts, L.S., Roberts, P.A., Adler, J., Davies, M.C., and Melia, C.D. ((1995). ). Determination of localized diffusion coefficients in gels using confocal scanning laser microscopy. J Microsc. 180, , 131-139.
Edidin, M., Zagyansky, Y., and Lardner, T.J. ((1976). ). Measurement of membrane protein lateral diffusion in single cells. Science 191, , 466-468.
Jain, R.K., Stock, R.J., Chary, S.R., and Rueter, M. ((1990). ). Convection and diffusion measurements using fluorescence recovery after photobleaching and video image analysis: in vitro calibration and assessment. Microvasc. Res. 39, , 77-93.[CrossRef][Medline]
Klonis, N., Rug, M., Harper, I., Wickham, M., Cowman, A., and Tilley, L. ((2002). ). Fluorescence photobleaching analysis for the study of cellular dynamics. Eur. Biophys. J. 31, , 36-51.[CrossRef][Medline]
Kubitscheck, U., Wedekind, P., and Peters, R. ((1994). ). Lateral diffusion measurement at high spatial resolution by scanning microphotolysis in a confocal microscope. Biophys J. 67, , 948-956.[Medline]
Kubitscheck, U., Wedekind, P., and Peters, R. ((1998). ). Three-dimensional diffusion measurements by scanning microphotolysis. J. Microsc. 192, , 126-138.[CrossRef]
Meyvis, T.K.L., De Smedt, S.C., Van Oostveldt, P., and Demeester, J. ((1999). ). Fluorescent recovery after photobleaching: a versatile tool for mobility and interaction measurements in pharmaceutical research. Pharm. Res. 16, , 1153-1162.[CrossRef][Medline]
Peters, R., and Kubitscheck, U. ((1999). ). Scanning microphotolysis: three-dimensional diffusion measurement and optical single-transporter recording. Methods 18, , 508-517.[CrossRef][Medline]
Peters, R., Peters, J., Tews, K.H., and Bahr, W. ((1974). ). A microfluorimetric study of translational diffusion in erythrocyte membranes. Biochem. Biophys. Acta 367, , 282-294.[Medline]
Phair, R.D., and Misteli, T. ((2000). ). High mobility of proteins in the mammalian cell nucleus. Nature 404, , 604-609.[CrossRef][Medline]
Phair, R.D., Scaffidi, P., Elbi, C., Vecerová, J., Dey, A., Ozato, K., Brown, D.T., Hager, G., Bustin, M., and Misteli, T. ((2004). ). Global nature of dynamic protein-chromatin interactions in vivo: three-dimensional genome scanning and dynamic interaction networks of chromatin proteins. Mol. Cell. Biol. 24, , 6393-6402.
Potma, E.O., de Boeij, W.P., Bosgraaf, L., Roelofs, J., van Haastert, P.J.M., and Wiersma, D.A. ((2001). ). Reduced protein diffusion rate by cytoskeleton in vegetative and polarized Dyctyostelium cells. Biophys. J. 81, , 2010-2019.[Medline]
Reits, E.A., and Neefjes, J.J. ((2001). ). From fixed to FRAP: measuring protein mobility and activity in living cells. Nat. Cell Biol. 3, , E145-147.[CrossRef][Medline]
Ribbeck, K., and Görlich, D. ((2001). ). Kinetic analysis of translocation through nuclear pore complexes. EMBO J. 20, , 1320-1330.[CrossRef][Medline]
Seksek, O., Biwersi, J., and Verkman, A.S. ((1997). ). Translational diffusion of macromolecule-sized solutes in cytoplasm and nucleus. J. Cell Biol. 138, , 131-142.
Shimi, T., Koujin, T., Segura-Totten, M., Wilson, K.L., Haraguchi, T., and Hiraoka, Y. ((2004). ). Dynamic interaction between BAF and emerin revealed by FRAP, FLIP, and FRET analyses in living HeLa cells. J. Struct. Biol. 147, , 31-41.[CrossRef][Medline]
Sheppard, C.J.R., and Torok, P. ((1997). ). Effects of specimen refractive index on confocal imaging. J. Microsc. 185, , 366-374.[CrossRef]
Soumpasis, D.M. ((1983). ). Theoretical analysis of fluorescence photobleaching recovery experiments. Biophys. J. 41, , 95-97.[Medline]
Sprague, B.L., Pego, R.L., Stavreva, D.A., and McNally, J.G. ((2004). ). Analysis of binding reactions by fluorescence recovery after photobleaching. Biophys. J. 86, , 3473-3495.[CrossRef][Medline]
Tsay, T.T., and Jacobson, K.A. ((1991). ). Spatial Fourier analysis of video photobleaching measurements–Principles and optimization. Biophys. J. 60, , 360-368.[Medline]
Wedekind, P., Kubitscheck, U., Heinrich, O., and Peters, R. ((1996). ). Line-scanning microphotolysis for diffraction-limited measurements of lateral diffusion. Biophys. J. 71, , 1621-1632.[Medline]
Wei, X., Henke, V.G., Strübing, C., Brown, E.B., and Clapham, D.E. ((2003). ). Real-time imaging of nuclear permeation by EGFP in single intact cells. Biophys. J. 84, , 1317-1327.[Medline]
White, J., and Stelzer, E. ((1999). ). Photobleaching GFP reveals protein dynamics inside live cells. Trends Cell Biol. 9, , 61-65. 3[CrossRef][Medline]
This article has been cited by other articles:
|
|
 |
|
  E. C. Freundt, L. Yu, E. Park, M. J. Lenardo, and X.-N. Xu Molecular Determinants for Subcellular Localization of the Severe Acute Respiratory Syndrome Coronavirus Open Reading Frame 3b Protein J. Virol., July 1, 2009; 83(13): 6631 - 6640. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. Rino, J. M. P. Desterro, T. R. Pacheco, T. W. J. Gadella Jr., and M. Carmo-Fonseca Splicing Factors SF1 and U2AF Associate in Extraspliceosomal Complexes Mol. Cell. Biol., May 1, 2008; 28(9): 3045 - 3057. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 D. Grunwald, B. Spottke, V. Buschmann, and U. Kubitscheck Intranuclear Binding Kinetics and Mobility of Single Native U1 snRNP Particles in Living Cells Mol. Biol. Cell, December 1, 2006; 17(12): 5017 - 5027. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 P. R. Darrah, M. Tlalka, A. Ashford, S. C. Watkinson, and M. D. Fricker The vacuole system is a significant intracellular pathway for longitudinal solute transport in basidiomycete fungi. Eukaryot. Cell, July 1, 2006; 5(7): 1111 - 1125. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 J. Perez-Vilar, J. C. Olsen, M. Chua, and R. C. Boucher pH-dependent Intraluminal Organization of Mucin Granules in Live Human Mucous/Goblet Cells J. Biol. Chem., April 29, 2005; 280(17): 16868 - 16881. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
