Analysis of fluorophore diffusion by continuous distributions of diffusion coefficients: application to photobleaching measurements of multicomponent and anomalous diffusion.

Fluorescence recovery after photobleaching (FRAP) is widely used to measure fluorophore diffusion in artificial solutions and cellular compartments. Two new strategies to analyze FRAP data were investigated theoretically and applied to complex systems with anomalous diffusion or multiple diffusing species: 1) continuous distributions of diffusion coefficients, alpha(D), and 2) time-dependent diffusion coefficients, D(t). A regression procedure utilizing the maximum entropy method was developed to resolve alpha(D) from fluorescence recovery curves, F(t). The recovery of multi-component alpha(D) from simulated F(t) with random noise was demonstrated and limitations of the method were defined. Single narrow Gaussian alpha(D) were recovered for FRAP measurements of thin films of fluorescein and size-fractionated FITC-dextrans and Ficolls, and multi-component alpha(D) were recovered for defined fluorophore mixtures. Single Gaussian alpha(D) were also recovered for solute diffusion in viscous media containing high dextran concentrations. To identify anomalous diffusion from FRAP data, a theory was developed to compute F(t) and alpha(D) for anomalous diffusion models defined by arbitrary nonlinear mean-squared displacement <x2> versus time relations. Several characteristic alpha(D) profiles for anomalous diffusion were found, including broad alpha(D) for subdiffusion, and alpha(D) with negative amplitudes for superdiffusion. A method to deduce apparent D(t) from F(t) was also developed and shown to provide useful complementary information to alpha(D). alpha(D) and D(t) were determined from photobleaching measurements of systems with apparent anomalous subdiffusion (nonuniform solution layer) and superdiffusion (moving fluid layer). The results establish a practical strategy to characterize complex diffusive phenomena from photobleaching recovery measurements.     

Three-dimensional fluorescence recovery after photobleaching with the confocal scanning laser microscope

Confocal scanning laser microscopes (CSLMs) are equipped with the feature to photobleach user-defined regions. This makes them a handy tool to perform fluorescence recovery after photobleaching (FRAP) measurements. To allow quantification of such FRAP experiments, a three-dimensional model has been developed that describes the fluorescence recovery process for a disk-shaped geometry that is photobleached by the scanning beam of a CSLM. First the general mathematical basis is outlined describing the bleaching process for an arbitrary geometry bleached by a scanning laser beam. Next, these general expressions are applied to the bleaching by a CSLM of a disk-shaped geometry and an analytical solution is derived that describes three-dimensional fluorescence recovery in the bleached area as observed by the CSLM. The FRAP model is validated through both the Stokes-Einstein relation and the comparison of the measured diffusion coefficients with their theoretical estimates. Finally, the FRAP model is used to characterize the transport of FITC-dextrans through bulk three-dimensional biological materials: vitreous body isolated from bovine eyes, and lung sputum expectorated by cystic fibrosis patients. The decrease in the diffusion coefficient relative to its value in solution was dependent on the size of the FITC-dextrans in vitreous, whereas it was size-independent in cystic fibrosis sputum.     

Accurate quantification of diffusion and binding kinetics of non-integral membrane proteins by FRAP

Non-integral membrane proteins frequently act as transduction hubs in vital signaling pathways initiated at the plasma membrane (PM). Their biological activity depends on dynamic interactions with the PM, which are governed by their lateral and cytoplasmic diffusion and membrane binding/unbinding kinetics. Accurate quantification of the multiple kinetic parameters characterizing their membrane interaction dynamics has been challenging. Despite a fair number of approximate fitting functions for analyzing fluorescence recovery after photobleaching (FRAP) data, no approach was able to cope with the full diffusion-exchange problem. Here, we present an exact solution and matlab fitting programs for FRAP with a stationary Gaussian laser beam, allowing simultaneous determination of the membrane (un)binding rates and the diffusion coefficients. To reduce the number of fitting parameters, the cytoplasmic diffusion coefficient is determined separately. Notably, our equations include the dependence of the exchange kinetics on the distribution of the measured protein between the PM and the cytoplasm, enabling the derivation of both k(on) and k(off) without prior assumptions. After validating the fitting function by computer simulations, we confirm the applicability of our approach to live-cell data by monitoring the dynamics of GFP-N-Ras mutants under conditions with different contributions of lateral diffusion and exchange to the FRAP kinetics.     

Challenges and artifacts in quantitative photobleaching experiments

Confocal fluorescence recovery after photobleaching (FRAP) is today the prevalent tool when studying the diffusional and kinetic properties of proteins in living cells. Obtaining quantitative data for diffusion coefficients via FRAP, however, is challenged by the fact that both bleaching and scanning take a finite time. Starting from an experimental case, it is shown by means of computer simulations that this intrinsic temporal limitation can lead to a gross underestimation of diffusion coefficients. Determining the binding kinetics of proteins to membranes with FRAP is further shown to be severely hampered by additional diffusional contributions, e.g. diffusion-limited binding. In some cases, the binding kinetics may even be masked entirely by diffusion. As current efforts to approach biological problems with biophysical models have to rely on experimentally determined model parameters, e.g. binding rates and diffusion constants, it is proposed that the accuracy in evaluating FRAP measurements can be improved by means of accompanying computer simulations.     

Mobility measurement by analysis of fluorescence photobleaching recovery kinetics.

Fluorescence photobleaching recovery (FPR) denotes a method for measuring two-dimensional lateral mobility of fluorescent particles, for example, the motion of fluorescently labeled molecules in approximately 10 mum2 regions of a single cell surface. A small spot on the fluorescent surface is photobleached by a brief exposure to an intense focused laser beam, and the subsequent recovery of the fluorescence is monitored by the same, but attenuated, laser beam. Recovery occurs by replenishment of intact fluorophore in the bleached spot by lateral transport from the surrounding surface. We present the theoretical basis and some practical guidelines for simple, rigorous analysis of FPR experiments. Information obtainable from FPR experiments includes: (a) identification of transport process type, i.e. the admixture of random diffusion and uniform directed flow; (b) determination of the absolute mobility coefficient, i.e. the diffusion constant and/or flow velocity; and (c) the fraction of total fluorophore which is mobile. To illustrate the experimental method and to verify the theory for diffusion, we describe some model experiments on aqueous solutions of rhodamine 6G.     

Fluorescence Recovery after Merging a Droplet to Measure the Two-dimensional Diffusion of a Phospholipid Monolayer

We introduce a new method to measure the lateral diffusivity of a surfactant monolayer at the fluid-fluid interface, called fluorescence recovery after merging (FRAM). FRAM adopts the same principles as the fluorescence recovery after photobleaching (FRAP) technique, especially for measuring fluorescence recovery after bleaching a specific area, but FRAM uses a drop coalescence instead of photobleaching dye molecules to induce a chemical potential gradient of dye molecules. Our technique has several advantages over FRAP: it only requires a fluorescence microscope rather than a confocal microscope equipped with high power lasers; it is essentially free from the selection of fluorescence dyes; and it has far more freedom to define the measured diffusion area. Furthermore, FRAM potentially provides a route for studying the mixing or inter-diffusion of two different surfactants, when the monolayers at a surface of droplet and at a flat air/water interface are prepared with different species, independently.     

Lowering the barriers to random walks on the cell surface

We used fluorescence recovery after photobleaching (FRAP) and single particle tracking (SPT) techniques to compare diffusion of class I major histocompatibility complex molecules (MHC) on normal and alpha-spectrin-deficient murine erythroleukemia (MEL) cells. Because the cytoskeleton mesh acts as a barrier to lateral mobility of membrane proteins, we expected that diffusion of membrane proteins in alpha-spectrin-deficient MEL cells would differ greatly from that in normal MEL cells. In the event, diffusion coefficients derived from either FRAP or SPT analysis were similar for alpha-spectrin-deficient and normal MEL cells, differing by a factor of approximately 2, on three different timescales: tens of seconds, 1-10 s, and 100 ms. SPT analysis showed that the diffusion of most class I MHC molecules was confined on both cell types. On the normal MEL cells, the mean diagonal length of the confined area was 330 nm with a mean residency time of 40s. On the alpha-spectrin-deficient MEL cells, the mean diagonal length was 650 nm with a mean residency time of 45s. Thus there are fewer barriers to lateral diffusion on cytoskeleton mutant MEL cells than on normal MEL cells, but this difference does not strongly affect lateral diffusion on the scales measured here.     

Intracellular macromolecular mobility measured by fluorescence recovery after photobleaching with confocal laser scanning microscopes

Fluorescence recovery after photobleaching (FRAP) is a widely used tool for estimating mobility parameters of fluorescently tagged molecules in cells. Despite the widespread use of confocal laser scanning microscopes (CLSMs) to perform photobleaching experiments, quantitative data analysis has been limited by lack of appropriate practical models. Here, we present a new approximate FRAP model for use on any standard CLSM. The main novelty of the method is that it takes into account diffusion of highly mobile molecules during the bleach phase. In fact, we show that by the time the first postbleach image is acquired in a CLSM a significant fluorescence recovery of fast-moving molecules has already taken place. The model was tested by generating simulated FRAP recovery curves for a wide range of diffusion coefficients and immobile fractions. The method was further validated by an experimental determination of the diffusion coefficient of fluorescent dextrans and green fluorescent protein. The new FRAP method was used to compare the mobility rates of fluorescent dextrans of 20, 40, 70, and 500 kDa in aqueous solution and in the nucleus of living HeLa cells. Diffusion coefficients were lower in the nucleoplasm, particularly for higher molecular weight dextrans. This is most likely caused by a sterical hindrance effect imposed by nuclear components. Decreasing the temperature from 37 to 22 degrees C reduces the dextran diffusion rates by approximately 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.     

Micrometer-scale domains in fibroblast plasma membranes

We have used the technique of fluorescence photobleaching recovery to measure the lateral diffusion coefficients and the mobile fractions of a fluorescent lipid probe, 1-acyl-2-(12-[(7-nitro-2-1, 3-benzoxadiazol-4-yl)aminododecanoyl]) phosphatidylcholine (NBD-PC), and of labeled membrane proteins of human fibroblasts. Values for mobile fractions decrease monotonically with increasing size of the laser spot used for the measurements, over a range of 0.35-5.0 microns. Values for NBD-PC diffusion coefficients increase in part of this range to reach a plateau at larger laser spots. This variation is not an artifact of the measuring system, since the effects are not seen if diffusion of the probe is measured in liposomes. We also find that the distribution of diffusion coefficients measured with small laser spots is heterogeneous indicating that these small spots can sample different regions of the membrane. These regions appear to differ in protein concentration. Our data strongly indicate that fibroblast surface membranes consist of protein-rich domains approximately 1 micron in diameter, embedded in a relatively protein-poor lipid continuum. These features appear in photographs of labeled cell surfaces illuminated by the expanded laser beam.     

The effect of ligand affinity on integrins’ lateral diffusion in cultured cells

The role of ligand affinity in altering αPS2CβPS integrins’ lateral mobility was studied using single particle tracking (SPT) with ligand-functionalized quantum dots (QDs) and fluorescence recovery after photobleaching (FRAP) with fluorescent protein tagged integrins. Integrins are ubiquitous transmembrane proteins that are vital for numerous cellular functions, including bidirectional signaling and cell anchorage. Wild-type and high ligand affinity mutant (αPS2CβPS-V409D) integrins were studied in S2 cells. As measured by SPT, the integrin mobile fraction decreased by 22 % and had a 4× slower diffusion coefficient for αPS2CβPS-V409D compared to wild-type integrins. These differences are partially the result of αPS2CβPS-V409D integrins’ increased clustering. For the wild-type integrins, the average of all diffusion coefficients measured by SPT was statistically similar to the ensemble FRAP results. A 75 % slower average diffusion coefficient was measured by SPT compared to FRAP for αPS2CβPS-V409D integrins, and this may be the result of SPT measuring only ligand-bound integrins, in contrast all ligand-bound and ligand-unbound integrins are averaged in FRAP measurements. Specific binding of the ligand-functionalized QDs was 99 % for integrin expressing cells. The results prove that the ligand binding affinity affects the lateral dynamics of a subset of integrins based on the complementary SPT and FRAP data.     

Effects of organelle shape on fluorescence recovery after photobleaching

The determination of diffusion coefficients from fluorescence recovery data is often complicated by geometric constraints imposed by the complex shapes of intracellular compartments. To address this issue, diffusion of proteins in the lumen of the endoplasmic reticulum (ER) is studied using cell biological and computational methods. Fluorescence recovery after photobleaching (FRAP) experiments are performed in tissue culture cells expressing GFP-KDEL, a soluble, fluorescent protein, in the ER lumen. The three-dimensional (3D) shape of the ER is determined by confocal microscopy and computationally reconstructed. Within these ER geometries diffusion of solutes is simulated using the method of particle strength exchange. The simulations are compared to experimental FRAP curves of GFP-KDEL in the same ER region. Comparisons of simulations in the 3D ER shapes to simulations in open 3D space show that the constraints imposed by the spatial confinement result in two- to fourfold underestimation of the molecular diffusion constant in the ER if the geometry is not taken into account. Using the same molecular diffusion constant in different simulations, the observed speed of fluorescence recovery varies by a factor of 2.5, depending on the particular ER geometry and the location of the bleached area. Organelle shape considerably influences diffusive transport and must be taken into account when relating experimental photobleaching data to molecular diffusion coefficients. This novel methodology combines experimental FRAP curves with high accuracy computer simulations of diffusion in the same ER geometry to determine the molecular diffusion constant of the solute in the particular ER lumen.     

Protein diffusion in mammalian cell cytoplasm

We introduce a new method for mesoscopic modeling of protein diffusion in an entire cell. This method is based on the construction of a three-dimensional digital model cell from confocal microscopy data. The model cell is segmented into the cytoplasm, nucleus, plasma membrane, and nuclear envelope, in which environment protein motion is modeled by fully numerical mesoscopic methods. Finer cellular structures that cannot be resolved with the imaging technique, which significantly affect protein motion, are accounted for in this method by assigning an effective, position-dependent porosity to the cell. This porosity can also be determined by confocal microscopy using the equilibrium distribution of a non-binding fluorescent protein. Distinction can now be made within this method between diffusion in the liquid phase of the cell (cytosol/nucleosol) and the cytoplasm/nucleoplasm. Here we applied the method to analyze fluorescence recovery after photobleach (FRAP) experiments in which the diffusion coefficient of a freely-diffusing model protein was determined for two different cell lines, and to explain the clear difference typically observed between conventional FRAP results and those of fluorescence correlation spectroscopy (FCS). A large difference was found in the FRAP experiments between diffusion in the cytoplasm/nucleoplasm and in the cytosol/nucleosol, for all of which the diffusion coefficients were determined. The cytosol results were found to be in very good agreement with those by FCS.     

Mobility in the mitochondrial electron transport chain

The role of lateral diffusion in mitochondrial electron transport has been investigated by measuring the diffusion coefficients for lipid, cytochrome c, and cytochrome oxidase in membranes of giant mitoplasts from cuprizone-fed mice using the technique of fluorescence redistribution after photobleaching (FRAP). The diffusion coefficient of the phospholipid analogue N-(7-nitro-2,1,3-benzoxadiazol-4-yl)phosphatidylethanolamine is dependent on the technique used to remove the outer mitochondrial membrane. A sonication technique yields mitoplasts with monophasic recovery of the lipid probe (D = 6 X 10(-9) cm2/s), while digitonin-treated mitochondria show biphasic recoveries (D1 = 5 X 10(-9) cm2/s; D2 = 1 X 10(-9) cm2/s). Digitonin appears to incorporate into mitoplasts, giving rise to decreased lipid mobility concomitant with increased rates of electron transfer from succinate to oxygen, in a manner reminiscent of the effects of cholesterol incorporation [Schneider, H., Lemasters, J. J., Hochli, M., & Hackenbrock, C. R. (1980) J. Biol. Chem. 255, 3748-3756]. FRAP measurements on tetramethylrhodamine cytochrome c modified at lysine-39 and on a mixture of active morpholinorhodamine derivatives of cytochrome c gave diffusion coefficients of (3.5-7) X 10(-10) cm2/s depending on the assay medium. With morpholinorhodamine-labeled antibodies purified on a cytochrome oxidase affinity column, the diffusion coefficient for cytochrome oxidase was determined to be 1.5 X 10(-10) cm2/s. The results are discussed in terms of a dynamic aggregate model in which an equilibrium exists between freely diffusing and associated electron-transfer components.     

Measurement of the lateral mobility of cell surface components in single living cells by fluorescence recovery after photobleaching

The use of fluorescence recovery after photobleaching (FRAP) techniques to monitor the lateral mobility of plant lectin-receptor complexes on the surface of single, living mammalian cells is described in detail. FRAP measurements indicate that over 75% of the wheat germ agglutinin receptor (WGA-receptor) complexes on the surface of human embryo fibroblasts are mobile. These WGA-receptor complexes diffuse laterally (as opposed to flow) on the cell surface with a diffusion coefficient in the range of 2 × 10^?11 to 2 × 10^?10 cm²/sec. Both the percentage of mobile WGA-receptor complexes and the mean diffusion coefficient of these complexes are higher than that obtained from earlier FRAP measurements of the mobility of concanavalin A-receptor (Con A-receptor) complexes in a variety of cell types. The possible reasons for the differing mobilities of WGA and Con A receptors are discussed.     

A Refined Reaction-Diffusion Model of Tau-Microtubule Dynamics and Its Application in FDAP Analysis

Fluorescence decay after photoactivation (FDAP) and fluorescence recovery after photobleaching (FRAP) are well established approaches for studying the interaction of the microtubule (MT)-associated protein tau with MTs in neuronal cells. Previous interpretations of FDAP/FRAP data have revealed dwell times of tau on MTs in the range of several seconds. However, this is difficult to reconcile with a dwell time recently measured by single-molecule analysis in neuronal processes that was shorter by two orders of magnitude. Questioning the validity of previously used phenomenological interpretations of FDAP/FRAP data, we have generalized the standard two-state reaction-diffusion equations by 1), accounting for the parallel and discrete arrangement of MTs in cell processes (i.e., homogeneous versus heterogeneous distribution of tau-binding sites); and 2), explicitly considering both active (diffusion upon MTs) and passive (piggybacking upon MTs at rates of slow axonal transport) motion of bound tau. For some idealized cases, analytical solutions were derived. By comparing them with the full numerical solution and Monte Carlo simulations, the respective validity domains were mapped. Interpretation of our FDAP data (from processes of neuronally differentiated PC12 cells) in light of the heterogeneous formalism yielded independent estimates for the association (∼2 ms) and dwell (∼100 ms) times of tau to/on a single MT rather than in an MT array. The dwell time was shorter by orders of magnitude than that in a previous report where a homogeneous topology of MTs was assumed. We found that the diffusion of bound tau was negligible in vivo, in contrast to an earlier report that tau diffuses along the MT lattice in vitro. Methodologically, our results demonstrate that the heterogeneity of binding sites cannot be ignored when dealing with reaction-diffusion of cytoskeleton-associated proteins. Physiologically, the results reveal the behavior of tau in cellular processes, which is noticeably different from that in vitro.     

A closed-form analytic expression for FRAP formula for the binding diffusion model

One of the most dominant methods cells use for a large class of cellular processes is reaction (or binding) diffusion kinetics, which are controlled by kinetic constants such as diffusion coefficients and on/off binding rate constants. Fluorescence recovery after photobleaching (FRAP) can be used to determine these kinetic constants in living cells. While an analytic expression for FRAP formulae for pure diffusion has been available for some time, an analytic FRAP formula for the binding diffusion model has not been reported yet. Here, we present an analytic FRAP formula for the binding diffusion model in an explicit form allowing for diffusion of the bound complex for either a uniform circle laser profile or a Gaussian laser profile.     

Measuring diffusion and binding kinetics by contact area FRAP

The immunological synapse is a stable intercellular structure that specializes in substance and signal transfer from one immune cell to another. Its formation is regulated in part by the diffusion of adhesion and signaling molecules into, and their binding of countermolecules in the contact area. The stability of immunological synapses allows receptor-ligand interactions to approximate chemical equilibrium despite other dynamic aspects. We have developed a mathematical model that describes the coupled reaction-diffusion process in an established immunological synapse. In this study, we extend a previously described contact area fluorescence recovery after photobleaching (FRAP) experiment to test the validity of the model. The receptor binding activity and lateral mobility of fluorescently labeled, lipid-anchored ligands in the bilayer resulted in their accumulation, as revealed by a much higher fluorescence intensity inside the contact area than outside. After complete photobleaching of the synapse, fluorescence recovery requires ligands to dissociate and rebind, and to diffuse in and out of the contact area. Such a FRAP time course consequently provides information on reaction and diffusion, which can be extracted by fitting the model solution to the data. Surprisingly, reverse rates in the two-dimensional contact area were at least 100-fold slower than in three-dimensional solution. As previously reported in immunological synapses, a significant nonrecoverable fraction of fluorescence was observed with one of two systems studied, suggesting some ligands either dissociated or diffused much more slowly compared with other ligands in the same synapse. The combined theory and experiment thus provides a new method for in situ measurements of kinetic rates, diffusion coefficients, and nonrecoverable fractions of interacting molecules in immunological synapses and other stable cell-bilayer junctions.     

Simplified Equation to Extract Diffusion Coefficients from Confocal FRAP Data

Quantitative measurements of diffusion can provide important information about how proteins and lipids interact with their environment within the cell and the effective size of the diffusing species. Confocal fluorescence recovery after photobleaching (FRAP) is one of the most widely accessible approaches to measure protein and lipid diffusion in living cells. However, straightforward approaches to quantify confocal FRAP measurements in terms of absolute diffusion coefficients are currently lacking. Here, we report a simplified equation that can be used to extract diffusion coefficients from confocal FRAP data using the half time of recovery and effective bleach radius for a circular bleach region, and validate this equation for a series of fluorescently labeled soluble and membrane‐bound proteins and lipids. We show that using this approach, diffusion coefficients ranging over three orders of magnitude can be obtained from confocal FRAP measurements performed under standard imaging conditions, highlighting its broad applicability.     

Regulation of the lateral diffusion of WGA-labeled glycoconjugates in human leukocytes. Comparison between adult granulocytes and differentiating promyelocytic HL60 cells

The regulation of the membrane mobility of glycoconjugates in human polymorphonuclear leukocytes (PMNL) was studied by comparing adult PMNL with promyelocytic HL60 cells before and after stimulation of differentiation in HL60 cells with phorbol-myristate acetate (PMA) with respect to lateral diffusion of wheat germ agglutinin (WGA)-labeled glycoconjugates. For this purpose we developed a novel variant of microscope equipment for the study of fluorescence recovery after photobleaching (FRAP) and continuous fluorescence microphotolysis (CFM) using a mini-computer for handling of shutters, data acquisition, and calculations. This equipment is presented in the report. We found that PMA-induced differentiation in HL60 cells reduced the lateral diffusion coefficient (D) of WGA-labeled membrane entities from about 1.5 to 1.0 x 10(-10) cm2/s, which was close to that found for adult blood PMNL, i.e., 1-1.2 x 10(-10) cm2/s. The lateral mobility (D x 10(10)) of succinylated WGA (S-WGA) was 2.3 and 1.7 cm2/s in undifferentiated and PMA-differentiated HL60 cells, respectively, indicating that WGA might have cross-linked membrane receptors, resulting in the slower diffusion. The results are discussed in relation to the effect of phagocyte maturation on the mobility of membrane components.