help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Clayton, A. H. A.
Right arrow Articles by Jovin, T. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Clayton, A. H. A.
Right arrow Articles by Jovin, T. M.

Biophys J, September 2002, p. 1631-1649, Vol. 83, No. 3

Dynamic Fluorescence Anisotropy Imaging Microscopy in the Frequency Domain (rFLIM)

Andrew H. A. Clayton, Quentin S. Hanley, Donna J. Arndt-Jovin, Vinod Subramaniam, and Thomas M. Jovin

Department of Molecular Biology, Max Planck Institute for Biophysical Chemistry, D-37077 Göttingen, Germany


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL
RESULTS
EXPERIMENTAL RESULTS
DISCUSSION
APPENDIX
REFERENCES

We describe a novel variant of fluorescence lifetime imaging microscopy (FLIM), denoted anisotropy-FLIM or rFLIM, which enables the wide-field measurement of the anisotropy decay of fluorophores on a pixel-by-pixel basis. We adapted existing frequency-domain FLIM technology for rFLIM by introducing linear polarizers in the excitation and emission paths. The phase delay and intensity ratios (AC and DC) between the polarized components of the fluorescence signal are recorded, leading to estimations of rotational correlation times and limiting anisotropies. Theory is developed that allows all the parameters of the hindered rotator model to be extracted from measurements carried out at a single modulation frequency. Two-dimensional image detection with a sensitive CCD camera provides wide-field imaging of dynamic depolarization with parallel interrogation of different compartments of a complex biological structure such as a cell. The concepts and technique of rFLIM are illustrated with a fluorophore-solvent (fluorescein-glycerol) system as a model for isotropic rotational dynamics and with bacteria expressing enhanced green fluorescent protein (EGFP) exhibiting depolarization due to homotransfer of electronic excitation energy (emFRET). The frequency-domain formalism was extended to cover the phenomenon of emFRET and yielded data consistent with a concentration depolarization mechanism resulting from the high intracellular concentration of EGFP. These investigations establish rFLIM as a powerful tool for cellular imaging based on rotational dynamics and molecular proximity.


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL
RESULTS
EXPERIMENTAL RESULTS
DISCUSSION
APPENDIX
REFERENCES

The fluorescence anisotropy (r) of intrinsic or extrinsic fluorophores provides valuable information about the state and environment of the corresponding biomolecular carrier on or within a cell. The steady-state anisotropy, <A><AC>r</AC><AC>&cjs1171;</AC></A>, of a distinct molecular species undergoing isotropic rotational diffusion is related to the excited state lifetime tau  and the rotational correlation time phi  according to the Jablonski equation (Jablonski, 1960),
<FR><NU>r<SUB><UP>o</UP></SUB></NU><DE><A><AC>r</AC><AC>&cjs1171;</AC></A></DE></FR>=1+<FR><NU>&tgr;</NU><DE>&phgr;</DE></FR>=1+&sfgr;, (1)
where ro is a limiting value (in the absence of rotation) given by the relative orientation of the absorption and emission transition moments, and sigma  is the ratio tau /phi . From Eq. 1, it follows that the anisotropy is a parameter with the capability of revealing changes in orientational distribution and in the excited state lifetime with great sensitivity. Such a circumstance arises readily in cell biology, e.g., from alterations in local microviscosity or other restrictions to diffusional motion, biochemical environment, complex formation, and molecular proximity manifested by hetero- or homo-transfer of electronic energy. Furthermore, the fluorescence anisotropy and lifetime are both intrinsic parameters, unlike the intensity signals used to determine them, and are thereby relatively insensitive to experimental factors such as light path and geometry. This fortunate circumstance permits more reliable comparisons between data obtained from different experiments and from different laboratories.

The measurement of fluorescence anisotropy in a microscope offers the possibility for characterizing very small domains on the surface or in the interior of the cell. For example, the cytoplasmic viscosity of living cells was determined by anisotropy imaging of fluid-phase fluorophores (Dix and Verkman, 1990). In addition, highly oriented samples, such as labeled protein constituents of the plasma membrane, can be readily visualized by anisotropy microscopy, as documented in the pioneering study by Axelrod (1979) of the orientational distribution of dyes bound to erythrocyte ghost membranes.

Measurements of dynamic depolarization or anisotropy decay provide additional and complementary information, such as rotational correlation times sensitive to the rotational volume and shape of the fluorophore. Anisotropy decay can reveal the multiplicity of rotational correlation times reflecting heterogeneity in the molecular population and in the size, shape, and internal motions of the fluorophore-biomolecule conjugate. Real molecules often exhibit molecular asymmetry, anisotropic rotational modes, and multiple lifetimes, and, in the cellular milieu, molecular heterogeneity is the rule rather than the exception. An adequate mathematical formulation of such systems is necessarily complex, as is already evident from Eq. 1 due to the dependence of <A><AC>r</AC><AC>&cjs1171;</AC></A> on three other parameters. The need to image the dynamic depolarization of fluorophores in cellular environments directly prompted the present investigation.

Determinations of rotational correlation times are performed either in the time or frequency domain. For example, using an excitation pulse as a forcing function, Eq. 1 is transformed into a decay process in which the fluorescence lifetime tau  is no longer linked to the correlation time phi  
r(t)=(r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB>)e<SUP><UP>−t/&phgr;</UP></SUP>+r<SUB>∞</SUB>, i(t)=e<SUP><UP>−t/&tgr;</UP></SUP>/&tgr;, (2)
where i(t) is the instantaneous emission intensity, normalized to a unit integrated value. In Eq. 2, the parameter rinfinity reflects the extent to which the rotational reorientation is hindered. An example is the case of a transmembrane protein undergoing molecular rotation in a lipid bilayer (Lipari and Szabo, 1980; Kinosita et al., 1982; Thevenin et al., 1994; Martin-Fernandez et al., 1998). For this case, <A><AC>r</AC><AC>&cjs1171;</AC></A> is obtained by integrating the intensity-weighted r(t) (Eq. 3); the same approach applies to heterogeneous systems and more complex (multiexponential, nonexponential) time courses for the decay of the excited state or rotational diffusion,
<A><AC>r</AC><AC>&cjs1171;</AC></A>=<LIM><OP>∫</OP><LL>0</LL><UL>∞</UL></LIM> r(t) · i(t) <UP>d</UP>t→<FR><NU>r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB></NU><DE><A><AC>r</AC><AC>&cjs1171;</AC></A>−r<SUB>∞</SUB></DE></FR>=1+&sfgr;. (3)
By itself, the fluorescence lifetime is of great value for the generation of contrast in fluorescence microscopy. Fluorescence lifetime imaging microscopy (FLIM) provides pixel-by-pixel discrimination among fluorophores differing in excited state kinetics due to distinct photophysical characteristics or microenvironments. FLIM has been implemented with time-domain (Wang et al., 1991; Ghiggino et al., 1992; Scully et al., 1997; Cole et al., 2000) and frequency-domain (Marriot et al., 1991; Lakowicz et al., 1992; Piston et al., 1992; Gadella et al., 1993; Squire et al., 2000; Verveer et al., 2000; Hanley et al., 2001) techniques.

We have undertaken the extension of the FLIM technology to include the generation of images based on the dynamic depolarization of fluorescence, and denote the approach featured in this report as anisotropy FLIM (rFLIM). rFLIM complements standard steady-state polarization microscopy (Axelrod, 1979) and corresponding static and dynamic polarization measurements performed in solution or suspension. It adds the imaging capability to previous systems that were limited to single-point anisotropy decay measurements in the microscope (Verkman et al., 1990; Martin-Fernandez et al., 1998; Tramier et al., 2000).

The implementation of rFLIM with two-dimensional (2D) image detection using a sensitive charge-coupled device (CCD) camera provides wide-field imaging of dynamic depolarization, e.g., with the intent of simultaneously interrogating different compartments of a cell. To this end, we modified a frequency-domain FLIM apparatus (Gadella et al., 1993; Schneider and Clegg, 1997; Hanley et al., 2001) by the introduction of linear polarizers in the excitation and emission paths (Fig. 1). In this report, we first describe our experimental implementation of rFLIM and then the theory of frequency-domain rotational dynamics extended to rFLIM. The latter is illustrated with dynamic depolarization measurements of fluorescein in glycerol/water mixtures and of enhanced green fluorescent protein (EGFP) in bacteria. Finally, we discuss other applications, limitations, and future perspectives of rFLIM.



View larger version (28K):
[in this window]
[in a new window]
 
FIGURE 1   Schematic of rFLIM apparatus. The microchannel plate intensifier was modulated at the photocathode. The signal generators driving the image intensifier and the acousto-optical modulator (AOM) were phase locked to a 10-MHz frequency standard. The relative phase (delta  phase) between the two was adjusted under computer control, cycling through a series of n phase steps of 360°/n. An image was recorded at each phase step. The zero-order light exiting the AOM was isolated from the higher orders by an iris and relayed to the microscope illumination system through a multimode optical fiber. The optical modulation frequency was double the AOM driving frequency.


    EXPERIMENTAL
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL
RESULTS
EXPERIMENTAL RESULTS
DISCUSSION
APPENDIX
REFERENCES

rFLIM apparatus

A schematic of the rFLIM instrument is shown in Fig. 1. The frequency domain fluorescence lifetime and rotational dynamics imaging instrument was built as an add-on to a standard research grade fluorescence microscope (E-600, Nikon, Düsseldorf, Germany). The FLIM system used as the basis for the measurement of rotational dynamics was similar to previously described instruments (Gadella et al., 1993; Schneider and Clegg, 1997; Hanley et al., 2001). A modulated image intensifier (model C5825, Hamamatsu Photonics, Herrsching, Germany) was installed on the camera port of the microscope. The modulation frequency and phase were set with a computer-controlled signal generator (model 2030; Marconi Instruments, St. Albans, UK). The intensified image was relayed to a 14-bit cooled CCD camera (model KX-2 containing a Kodak KAF-1600 CCD with 9-µm square pixels in a 1536 × 1024 format; Apogee Instruments, Tucson, AZ) using a tandem pair of f/1.8 (AF Nikkor, Nikon) camera lenses. The 488-nm line of an Ar+ ion laser (Innova 90-5; Coherent, Santa Clara, CA) was used to illuminate the microscope field. The laser light was modulated by an acousto-optic modulator (AOM) (model SWM-10044; IntraAction Corp., Belwood, IL). A second signal generator (model 2023; Marconi) was frequency and phase locked to the intensifier signal generator. Its output signal was amplified with a 37-dB broadband linear power amplifier (model 403LA; Electronic Navigation Industries, Rochester, NY) and used to drive the AOM at a frequency in most experiments of 28.84 MHz, leading to optical modulation at f = 57.68 MHz. The zero-order light from the AOM was selected and relayed through a multimode optical fiber (model CU-3578 UV-vis; Multimode Fiber Optics, East Hanover, NJ) to the microscope illumination port. The optical fiber was vibrated mechanically at audio frequencies to scramble modes and thus reduce speckle. The dichroic filter was a 505 DRLP02 and the emission longpass filter (>515 nm) a 515EFLP (Omega Optical, Brattleboro, VT).

Anisotropy measurements were implemented by the addition of a sheet polarizer (Y-FA, Nikon) to the excitation light source and an analyzer (PW44, Schott, Mainz, Germany; or a Polaroid film polarizer) on the emission path of the microscope oriented either parallel or perpendicular to the excitation beam polarizer. Three data sets were taken; a conventional lifetime image with both excitation and emission polarizers removed; a 0°/0° (excitation/emission) polarization image; and a 0°/90° polarization image. A lifetime and anisotropy standard, 1 µM aqueous rhodamine 6G (Hanley et al., 2001), was also measured to provide a pixel-by-pixel calibration of relative modulation depth, phase shift, and polarization dependent transmissions (G-factor, Eq. 4 below).

Polarization confocal laser scanning microscopy

Confocal polarization-dependent images of the EGFP-expressing bacteria were acquired with a modified LSM 310 confocal laser scanning microscope (model LSM310; Zeiss, Göttingen, Germany), equipped with an external Ar+ ion laser (80 mW, 488 nm) as a photobleaching source. Fluorescence images of the parallel and perpendicular emission of the bacterial embedded in agarose were recorded with a 63× (1.2 NA) C-Apochromat water immersion objective, using excitation at 488 nm and emission at >515 nm. The unpolarized emission of a 0.2 µM fluorescein solution served to calculate the instrument G-factor (1.23) under the same conditions.

The total fluorescence and anisotropy images were generated with Scil-Image (TNO TPD, Delft, The Netherlands) software after background subtraction and pixel registration of the two polarized emission images by cross-correlation analysis. The final anisotropy image was smoothed by a 3 × 3 uniform filter function. Background pixels were excluded by threshold masking. The mean r for each anisotropy image was calculated from the mean intensity-weighted pixel values. Measurements were performed at ambient temperature (23°C).

Data acquisition and processing in rFLIM

In our implementation of rFLIM, the excitation light and the detector (intensifier) gain were modulated by the same radial frequency omega . A series of images was acquired while adjusting the relative phase between the AOM and the intensifier (Fig. 1) over one period (2pi radians, 360°). The images were taken in sets of n (n = 4, 8, or 16) such that each successive image was measured with an incremental phase shift of 360°/n. Each series of images was accompanied by a corresponding blank image acquired with all shutters open except for the laser source. Acquisition proceeded with a reversal of the sequence of phase shifts for achieving a first-order photobleaching compensation; the second data series was processed as a second period of the sampled sinusoidal waveform. Acquisition was the same for sample, standard, and all polarization states. Two sets of measurements were made with the emission polarizer oriented either parallel or perpendicular. Each data set was processed (Hanley et al., 2001) to yield frequency-dependent amplitudes, and a phase shift between the excitation and the emitted light. The two resulting polarized emission components are modulated at the same frequency but phase shifted relative to each other; I||,AC leads Iperp ,AC by Delta Phi . In addition, the AC amplitude of the perpendicular component is reduced relative to that of the parallel component. These signals are characterized by the ratio YAC = I||,AC/Iperp ,AC, and the corresponding ratio of DC (omega  = 0) magnitudes by YDC = I||,DC/Iperp ,DC.

Processing of the data involved five steps. 1) The lifetime images were computed from data obtained for the sample and lifetime reference in the absence of polarization optics. 2) The Delta Phi , YAC, and YDC images were calculated for the sample data sets. 3) The equivalent standardization images were computed for use as diagnostics and for determining G-factor corrections. 4) The latter were applied to the YAC, and YDC images according to the procedures given below. 5) Images of the anisotropy parameters (phi rorinfinity ) were generated according to the formalism developed in the Results section (Eqs. 14, 18-20, 24-27).

The pixel-by-pixel corrections to the amplitude images (YAC, YDC) were based on measurements of a reference solution with known anisotropy. G-factors (Gi) accounting for differences in system response for the two polarized emission components were computed at each (ith) image pixel using a reference solution of known anisotropy (usually ~0), according to
G<SUB><UP>i</UP></SUB>=<FENCE><FR><NU>1+2<A><AC>r</AC><AC>&cjs1171;</AC></A><SUB><UP>ref</UP></SUB></NU><DE>1−<A><AC>r</AC><AC>&cjs1171;</AC></A><SUB><UP>ref</UP></SUB></DE></FR></FENCE> <FR><NU>1</NU><DE>Y<SUP>ref</SUP><SUB>DC,i</SUB></DE></FR>, (4)
where <A><AC>r</AC><AC>&cjs1171;</AC></A>ref is the steady-state anisotropy of the reference solution measured in a spectrofluorometer and Y<UP><SUB>DC,i</SUB><SUP>ref</SUP></UP> is the DC intensity ratio (I||,i/Iperp ,i) of the polarized emission measured in the microscope (see also Dix and Verkman, 1990). The reference measurements were made either directly before or after the sample measurement with identical optical configuration (objective, dichroic filter, etc.). The values of Gi were 1.15 ± 0.08 [Plan 2× (NA 0.06) air objective] and differed by less than 0.02 among the different objectives tested [Plan 2×, PlanFluor 20× (NA 0.75) air, PlanFluor 40× (NA 0.75) air, Plan Apo 60× (NA 1.2) water]. The corrected AC and DC ratios of the sample were obtained from
Y<SUP><UP>corr</UP></SUP><SUB>(<UP>AC,DC</UP>)<UP>,i</UP></SUB>=G<SUB><UP>i</UP></SUB> · Y<SUP><UP>meas</UP></SUP><SUB>(<UP>AC,DC</UP>)<UP>,i</UP></SUB>. (5)
We observed a zero phase shift (within experimental error) for excitation with unpolarized light. Consequently, no corrections for instrumental effects on the phase shift were required. The programs Scil-Image (TNO TPD) and Mathematica 4.0 (Wolfram Research, Champaign, IL), including the Digital Imaging Processing add-on, were used for further image processing (segmentation and masking), generation of frequency histograms, and statistical analyses. Regions of interest (corresponding to individual cuvettes in a dual cuvette image or a particular cellular compartment in the cells) were isolated using an image intensity threshold mask.

Reference measurements in a spectrofluorometer

The fluorescence anisotropies of fluorescein solutions in glycerol-buffer mixtures and of rhodamine 6G in water were determined at 23°C with a commercial spectrofluorometer (enhanced model C-60, Photon Technology International, Monmouth Junction, NJ). The excitation wavelength was 488 nm, the emission wavelength 540 nm, and the slit spectral bandwidth 5 nm. For the rhodamine 6G solution, <A><AC>r</AC><AC>&cjs1171;</AC></A>ref was 0.012.

The fluorescence anisotropy and relative peak intensities of EGFP solutions were determined at 20°C with a Cary Eclipse spectrofluorometer (Varian Australia Pty, Melbourne) equipped with motorized excitation and emission polarizers. The excitation wavelength varied between 450 and 480 nm, and the emission spectra were measured up to 600 nm with a 5-nm bandwidth. A Suprasil cuvette with a 0.25-mm square cross section (<1 µL fill volume) was used to minimize inner filter and reabsorption effects. The intensity values were corrected for absorption A along the lightpath by the factor A ln[10]/(1 - 10-A). Anisotropy values determined in the two spectrofluorimeters agreed within 5%.

Materials and EGFP-expressing bacteria

Glycerol-buffer solutions over a range of glycerol concentrations from 10 to 80% (w/w) were prepared gravimetrically using spectroscopic grade glycerol (Aldrich, Milwaukee, WI) and 10 mM Tris-HCl, pH 8.0. Solutions for rFLIM rotational correlation time measurements contained 10, 61, and 70% glycerol by weight. The viscosities (eta , 23°C) of these solutions (1.2, 10.6, 21.3 cP, respectively) were estimated from published tables (Weast, 1979). Solutions for viscosity-dependent fluorescence measurements were prepared by addition of a small volume (<50 µL) of ~2 mM fluorescein in ethanol to a cuvette containing the glycerol/buffer mixture. The solutions were not deoxygenated. EGFP solutions were prepared volumetrically in 0.1 M Tris-HCl, 0.1 M NaCl, pH 8.5 buffer by dilution from a stock solution of known concentration, based on the extinction coefficient at 490 nm of 55 mM-1cm-1 (Patterson et al., 1997); we used varepsilon 480nm = 50 mM-1.

Escherichia coli bacteria expressing EGFP were cultured at 37°C as described previously (Jakobs et al., 2000). Before rFLIM measurements, the bacteria were suspended in 1% agarose to minimize movement and mounted between a microscope slide and a coverslip. Some of the EGFP leaked into the extracellular medium. A sample from the same bacterial culture was centrifuged and the supernatant was found to contain free EGFP, the spectral properties of which were identical to those obtained by deliberate lysis of the bacteria via repeated freeze-thaw cycles. The mean concentration of intracellular EGFP was estimated by filtration of the bacterial suspension over a Millipore HA 0.45 µ nitrocellulose filter. The weight of the packed cells corresponding to a given volume of culture medium was determined and the EGFP fluorescence of the suspension and filtrate compared to a known (2.8 µM EGFP) reference solution.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL
RESULTS
EXPERIMENTAL RESULTS
DISCUSSION
APPENDIX
REFERENCES

Theory

We first examine the theory of anisotropy and dynamic depolarization as applied to standard cuvette-based measurements of randomly orientated fluorophores using a homodyne measurement scheme. In conventional instruments, the sample is excited with linearly polarized light. The emission is detected perpendicular to the excitation propagation direction and viewed through an analyzing polarizer oriented either parallel (I||) or perpendicular (Iperp ) to the excitation beam polarizer. These two signals differ in general because the polarized incident beam generates an anisotropic distribution of excited state fluorophores, i.e., by photoselection (Albrecht, 1961). Rotational diffusion of the fluorophores during the excited state lifetime reduces the disparity between I|| and Iperp , resulting in dynamic depolarization of the emitted light according to Eq. 2. The fluorescence anisotropy r(t) at any time point is defined as (Jablonski, 1960)
r(t)=<FR><NU>I<SUB>∥</SUB>(t)−I<SUB>⊥</SUB>(t)</NU><DE>I<SUB>∥</SUB>(t)+2I<SUB>⊥</SUB>(t)</DE></FR>. (6)
The theory underlying the determination of rotational correlation times, differential polarized phase fluorometry, a method introduced by Weber and coworkers, is found in several excellent papers and reviews (Weber, 1977, 1978; Mantulin and Weber, 1977; Lakowicz et al., 1979; Lakowicz 1983; Jameson and Hazlett, 1986). The three measured parameters (Delta Phi , YAC, YDC) reflect the excited state lifetime and/or the intrinsic rotational diffusion properties of the fluorophore (and its conjugate) (Eq. 2), and an experimental variable (the modulation frequency omega ) that can be manipulated at will.

For the hindered rotator represented by Eq. 2 and using our nomenclature, one obtains
&Dgr;&PHgr;=<UP>tan<SUP>−1</SUP></UP><FENCE><FR><NU>3 · &ohgr;&tgr; · &sfgr; · (r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB>)</NU><DE><AR><R><C>(1−r<SUB><UP>o</UP></SUB>)(1+2r<SUB><UP>o</UP></SUB>)[<UP>1+</UP>(<UP>&ohgr;&tgr;</UP>)<SUP><UP>2</UP></SUP>]</C></R><R><C><UP>+</UP>[<UP>2+r<SUB>o</SUB></UP>−r<SUB>∞</SUB>(4r<SUB><UP>o</UP></SUB>−1)]&sfgr;</C></R><R><C>   +(1−r<SUB>∞</SUB>)(1+2r<SUB>∞</SUB>)&sfgr;<SUP>2</SUP></C></R></AR></DE></FR></FENCE>, (7)

Y<SUB><UP>AC</UP></SUB>=<RAD><RCD><FR><NU>(1+2r<SUB><UP>o</UP></SUB>)<SUP>2</SUP>(&ohgr;&tgr;)<SUP>2</SUP>+[(1+2r<SUB><UP>o</UP></SUB>)+(1+2r<SUB>∞</SUB>)&sfgr;]<SUP>2</SUP></NU><DE>(1−r<SUB><UP>o</UP></SUB>)<SUP>2</SUP>(&ohgr;&tgr;)<SUP>2</SUP>+[(1−r<SUB><UP>o</UP></SUB>)+(1−r<SUB>∞</SUB>)&sfgr;]<SUP>2</SUP></DE></FR></RCD></RAD>, (8)

Y<SUB><UP>DC</UP></SUB>=<FR><NU>1+2r<SUB><UP>o</UP></SUB>+(1+2r<SUB>∞</SUB>)&sfgr;</NU><DE>1−r<SUB><UP>o</UP></SUB>+(1−r<SUB>∞</SUB>)&sfgr;</DE></FR>; <A><AC>r</AC><AC>&cjs1171;</AC></A>=<FR><NU>Y<SUB><UP>DC</UP></SUB>−1</NU><DE>Y<SUB><UP>DC</UP></SUB>+2</DE></FR>. (9)
Note that we have used the dimensionless variables omega tau and sigma , incorporating tau  and phi . That is, the rotational correlation time is expressed relative to the lifetime and the latter serves to scale the radial frequency. Our parameters differ from those used in the early publications (R = 1/(6phi ), Weber, 1977, 1978; Lakowicz et al., 1979). The relationships of Eqs. 7-9 are depicted in Fig. 2 as a function of omega tau and sigma . The phase difference Delta Phi approaches 0 at both extremes of frequency (omega  = 0, infinity ) and passes through a maximum, whereas YAC increases monotonically with omega , albeit to a degree dependent on sigma ; YDC is not a function of frequency. The maximal value of Delta Phi , for any given parameter set (ro, rinfinity , sigma ), is given by (Weber, 1978)


&Dgr;&PHgr;<SUB><UP>max,&ohgr;&tgr;</UP></SUB><UP>=tan<SUP>−1</SUP></UP><FENCE><FR><NU>3(r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB>)&sfgr;</NU><DE>2<RAD><RCD>(1−r<SUB><UP>o</UP></SUB>)(1+2r<SUB><UP>o</UP></SUB>)[1−r<SUB><UP>o</UP></SUB>+(1−r<SUB>∞</SUB>)&sfgr;][1+2r<SUB><UP>o</UP></SUB>+(1+2r<SUB>∞</SUB>)&sfgr;]</RCD></RAD></DE></FR></FENCE>, (10)



View larger version (56K):
[in this window]
[in a new window]
 
FIGURE 2   Dependence of rFLIM measured parameters (Delta Phi , YAC, YDC) and their partial derivatives on omega tau and sigma . See text for discussion.

with
&ohgr;&tgr;<SUB><UP>&Dgr;&PHgr;max</UP></SUB>=<RAD><RCD><FENCE>1+<FR><NU>(1−r<SUB>∞</SUB>)</NU><DE>(1−r<SUB><UP>o</UP></SUB>)</DE></FR> &sfgr;</FENCE> <FENCE>1+<FR><NU>(1+2r<SUB>∞</SUB>)</NU><DE>(1+2r<SUB><UP>o</UP></SUB>)</DE></FR> &sfgr;</FENCE></RCD></RAD>≈1+&sfgr;, (11)
and, for any given parameter set (ro, rinfinity , omega tau ), by


&Dgr;&PHgr;<SUB><UP>max,&sfgr;</UP></SUB>=<UP>tan<SUP>−1</SUP></UP><FENCE><FR><NU>3(r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB>)&ohgr;&tgr;</NU><DE>2+r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB>(4r<SUB><UP>o</UP></SUB>−1)+2<RAD><RCD>[(1−r<SUB><UP>o</UP></SUB>)(1+2r<SUB><UP>o</UP></SUB>)(1−r<SUB>∞</SUB>)(1+2r<SUB>∞</SUB>)(1+(&ohgr;&tgr;)<SUP>2</SUP>)]</RCD></RAD></DE></FR></FENCE>, (12)

with
&sfgr;<SUB><UP>&Dgr;&PHgr;max</UP></SUB>=<RAD><RCD><FR><NU>(1−r<SUB><UP>o</UP></SUB>)(1+2r<SUB><UP>o</UP></SUB>)[1+(&ohgr;&tgr;)<SUP>2</SUP>]</NU><DE>(1−r<SUB>∞</SUB>)(1+2r<SUB>∞</SUB>)</DE></FR></RCD></RAD> (13)

≈<RAD><RCD>1+(&ohgr;&tgr;)<SUP>2</SUP></RCD></RAD>.
The approximations of Eqs. 11 and 13 hold well (±10%) over a large range of ro and rinfinity values. Thus, a frequency sweep of Delta Phi can provide an estimate of sigma  (Eq. 11) and, thereby, of phi .

The dependence of Delta Phi , YAC, and YDC on sigma  is also featured in the three-dimensional (3D) representations of Fig. 2. Delta Phi is a bell-shaped function of sigma  (a limited region of parameter space is shown in Fig. A1). YDC and YAC decrease monotonically with sigma . (These properties apply for ro > rinfinity  > 0.) Important features of these parametric functions are revealed by the partial derivatives of Delta Phi , YAC, and YDC with respect to omega tau and sigma , which are also functions of the same two parameters (Fig. 2). The greatest measurement sensitivity is achieved for the largest magnitudes of the derivatives, i.e., at low values of sigma  and omega . However, the optimal conditions for determining phi  are not necessarily the same as those for tau . It follows that a judicious choice of experimental conditions is required for the most sensitive determination of phi .

The expressions corresponding to an isotropic rotator are given by setting rinfinity to 0 in Eqs. 7-9 (and in Eqs. 14-19 below). In the practical application of Eqs. 7-9, tau  is obtained from independent measurements, and phi  is derived from various strategies for achieving solutions for sigma  (and thus phi  via the definition sigma  = tau /phi ). The approaches depend on the model and differ according to the degree of prior knowledge about the other depolarization parameters (ro and rinfinity ) and the selection of signals (Eqs. 7-9) to be used alone or in combination.

The general solutions of Eqs. 7-9 for sigma  (and thus for phi ) consist of the following expressions (Eqs. 14-19). The program Mathematica was used to perform extensive symbolic manipulations and simulations of these functions. The solution derived from difference phase measurements is given by
&phgr;<SUB>&Dgr;&PHgr;</SUB>=<FR><NU>a</NU><DE>1+<UP>sign</UP><RAD><RCD>b</RCD></RAD></DE></FR> · &tgr; (14)
with
a=<FR><NU>2(1−r<SUB>∞</SUB>)(1+2r<SUB>∞</SUB>)<UP>tan</UP> &Dgr;&PHgr;</NU><DE>3(r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB>)&ohgr;&tgr;−[2+r<SUB><UP>o</UP></SUB>(1−4r<SUB>∞</SUB>)+r<SUB>∞</SUB>]<UP>tan</UP> &Dgr;&PHgr;</DE></FR>, (15)

b=1−<FR><NU>(1−r<SUB><UP>o</UP></SUB>)(1+2r<SUB><UP>o</UP></SUB>)</NU><DE>(1−r<SUB>∞</SUB>)(1+2r<SUB>∞</SUB>)</DE></FR> [1+(&ohgr;&tgr;)<SUP>2</SUP>]a<SUP>2</SUP> (16)

=1−&sfgr;<SUP><UP>2</UP></SUP><SUB><UP>&Dgr;&PHgr;max</UP></SUB>a<SUP>2</SUP>,

<UP>sign</UP>=<UP>− if</UP> &sfgr;<&sfgr;<SUB><UP>&Dgr;&PHgr;max</UP></SUB> (<UP>Eq. 13</UP>); (17)

+<UP>for</UP> &sfgr;>&sfgr;<SUB><UP>&Dgr;&PHgr;max</UP></SUB>.
In experimental practice, it is best to select a frequency for which sigma Delta Phi max is not close to sigma  because b, the argument of the square root function of Eq. 14, vanishes for this condition. Experimental error in the determination of Delta Phi can lead to fluctuations of b about 0, including negative values, and thus render evaluations of Eq. 14 problematical due to the presence of complex numbers. However, a good approximation holds for |bapprox  0: phi Delta Phi  approx  tau /sigma Delta Phi max. For sigma Delta Phi max not equal  sigma , b can still be driven negative in the event that experimental noise generates a recorded Delta Phi exceeding the value of Delta Phi exact (given by Eq. 7) by a factor > [(Delta Phi max,sigma /Delta Phi exact- 1]. In such cases, phi Delta Phi is mathematically and physically indeterminate, inasmuch as all solutions are imaginary.

The solutions derived from the AC and DC components of the relative polarized signals are given by


&phgr;<SUB><UP>Y<SUB>AC</SUB></UP></SUB>=<FR><NU>Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>(1−r<SUB>∞</SUB>)<SUP>2</SUP>−(1+2r<SUB>∞</SUB>)<SUP>2</SUP></NU><DE><AR><R><C>(1+2r<SUB><UP>o</UP></SUB>)(1+2r<SUB>∞</SUB>)−Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>(1−r<SUB><UP>o</UP></SUB>)(1−r<SUB>∞</SUB>)</C></R><R><C> +<RAD><RCD><AR><R><C>Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB><FENCE>9(r<SUB><UP>o</UP></SUB>−r<SUB>∞</SUB>)<SUP>2</SUP>+<FENCE><AR><R><C>2+r<SUB>∞</SUB>(2+5r<SUB>∞</SUB>)+r<SUB><UP>o</UP></SUB>[2−4r<SUB>∞</SUB>(4+r<SUB>∞</SUB>)]</C></R><R><C>+r<SUP><UP>2</UP></SUP><SUB><UP>o</UP></SUB>[5−4r<SUB>∞</SUB>(1−2r<SUB>∞</SUB>)]−Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>(1−r<SUB><UP>o</UP></SUB>)<SUP>2</SUP>(1−r<SUB>∞</SUB>)<SUP>2</SUP></C></R></AR>
</FENCE>(&ohgr;&tgr;)<SUP>2</SUP></FENCE></C></R><R><C> −(1+2r<SUB><UP>o</UP></SUB>)<SUP>2</SUP>(1+2r<SUB>∞</SUB>)<SUP>2</SUP>(&ohgr;&tgr;)<SUP>2</SUP></C></R></AR>
</RCD></RAD></C></R></AR></DE></FR>·&tgr; (18)


&phgr;<SUB><UP>Y<SUB>DC</SUB></UP></SUB>=<FENCE><FR><NU>Y<SUB><UP>DC</UP></SUB>(1−r<SUB>∞</SUB>)−(1+2r<SUB>∞</SUB>)</NU><DE>(1+2r<SUB><UP>o</UP></SUB>)−Y<SUB><UP>DC</UP></SUB>(1−r<SUB><UP>o</UP></SUB>)</DE></FR></FENCE> · &tgr;. (19)
In the simplest application of Eqs. 14-19, rinfinity is set to zero and ro is fixed at an independently-determined value. Next, the rotational correlation times derived from the three experimental signals are compared. In general, agreement among the three values indicates that the system behaves as an isotropic rotator. Disagreement implies more complex modes of anisotropy decay. Values determined in this manner are apparent rotational correlation times and serve as a very simple measure of rotational correlation time heterogeneity. In the event that either ro or rinfinity are known, different solutions can be derived by combination of any two of Eqs. 7-9. For example, for the hindered rotator model, knowledge of ro (a relatively invariant photophysical quantity) provides a solution for variable Delta Phi (i.e., phi ) and rinfinity using the Delta Phi and YDC (or as an equivalent, <A><AC>r</AC><AC>&cjs1171;</AC></A>, in Eq. 9) signals.
&phgr;<SUB><UP>&Dgr;&PHgr;,<A><AC>r</AC><AC>&cjs1171;</AC></A></UP></SUB>=<FR><NU>a</NU><DE>1+<UP>sign</UP><RAD><RCD>b</RCD></RAD></DE></FR> · &tgr;, (20)
with
a=<FENCE><FR><NU>3(r<SUB><UP>o</UP></SUB>−<A><AC>r</AC><AC>&cjs1171;</AC></A>)&ohgr;&tgr;</NU><DE>2(1−<A><AC>r</AC><AC>&cjs1171;</AC></A>)(1+2<A><AC>r</AC><AC>&cjs1171;</AC></A>)<UP>tan</UP> &Dgr;&PHgr;</DE></FR>−1</FENCE><SUP>−1</SUP>, (21)

b= (22)

<FR><NU>9(<A><AC>r</AC><AC>&cjs1171;</AC></A>−r<SUB><UP>o</UP></SUB>)<SUP>2</SUP>−4(1−<A><AC>r</AC><AC>&cjs1171;</AC></A>)(1+2<A><AC>r</AC><AC>&cjs1171;</AC></A>)(1−r<SUB><UP>o</UP></SUB>)(1+2r<SUB><UP>o</UP></SUB>)<UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></NU><DE>[2(1−<A><AC>r</AC><AC>&cjs1171;</AC></A>)(1+2<A><AC>r</AC><AC>&cjs1171;</AC></A>)<UP>tan</UP> &Dgr;&PHgr;−3(r<SUB><UP>o</UP></SUB>−<A><AC>r</AC><AC>&cjs1171;</AC></A>)&ohgr;&tgr;]<SUP>2</SUP></DE></FR> (&ohgr;&tgr;)<SUP>2</SUP>,

<UP>sign</UP>=+ <UP>if</UP> <FR><NU><FENCE>1+<RAD><RCD>b</RCD></RAD></FENCE>&tgr;</NU><DE>(r<SUB><UP>o</UP></SUB>/<A><AC>r</AC><AC>&cjs1171;</AC></A>−1)</DE></FR>>a><FENCE><FR><NU>1+<RAD><RCD>b</RCD></RAD></NU><DE>1−<RAD><RCD>b</RCD></RAD></DE></FR></FENCE> a (23)

− <UP>otherwise</UP>,

r<SUB>∞</SUB>=<FR><NU>(<A><AC>r</AC><AC>&cjs1171;</AC></A>−r<SUB><UP>o</UP></SUB>)</NU><DE>&sfgr;<SUB><UP>&Dgr;&PHgr;,<A><AC>r</AC><AC>&cjs1171;</AC></A></UP></SUB></DE></FR>+<A><AC>r</AC><AC>&cjs1171;</AC></A>=(<A><AC>r</AC><AC>&cjs1171;</AC></A>−r<SUB><UP>o</UP></SUB>)<FENCE><FR><NU>&phgr;<SUB><UP>&Dgr;&PHgr;,<A><AC>r</AC><AC>&cjs1171;</AC></A></UP></SUB></NU><DE>&tgr;</DE></FR></FENCE>+<A><AC>r</AC><AC>&cjs1171;</AC></A>. (24)
A similar approach was advocated by Lakowicz et al. (1979) to determine hindered rotations of a membrane probe in lipid bilayers.

Composite expressions for rotational diffusion parameters

Equations related to those presented above have appeared in the literature in different forms. They share the disadvantage of requiring some prior knowledge of ro and/or rinfinity , information that may not be readily available in microscope-based studies. However, combining all three parameters (Delta Phi , YAC, YDC; Eqs. 7-9) and tau  leads to independent "composite" expressions (Eqs. 25-27) for the three rotational diffusion parameters.


&phgr;<SUP><UP>comp</UP></SUP>=<FENCE><FR><NU>Y<SUB><UP>AC</UP></SUB><FENCE>(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−2Y<SUB><UP>DC</UP></SUB>)(Y<SUB><UP>DC</UP></SUB>+2)<RAD><RCD>1+<UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></RCD></RAD>+Y<SUB><UP>AC</UP></SUB>(Y<SUB><UP>DC</UP></SUB>−4)</FENCE><UP>tan</UP> &Dgr;&PHgr; · &ohgr;&tgr;</NU><DE>(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−4)(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>)+(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−2Y<SUB><UP>DC</UP></SUB>)<SUP>2</SUP> <UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></DE></FR>−1</FENCE><SUP>−1</SUP> · &tgr;, (25)

r<SUP><UP>comp</UP></SUP><SUB><UP>o</UP></SUB>=<FR><NU>(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>+2)(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>)+(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−2Y<SUB><UP>DC</UP></SUB>)(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>+Y<SUB><UP>DC</UP></SUB>)<UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP>−3Y<SUB><UP>AC</UP></SUB>(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>)<RAD><RCD>1+<UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></RCD></RAD></NU><DE>(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−4)(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>)+(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−2Y<SUB><UP>DC</UP></SUB>)<SUP>2</SUP> <UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></DE></FR>, (26)

r<SUP><UP>comp</UP></SUP><SUB><UP>∞</UP></SUB>=<FR><NU><FENCE><AR><R><C>(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>+2)(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>)−3Y<SUB><UP>AC</UP></SUB>[Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>−(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>+Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>)<UP>tan</UP> &Dgr;&PHgr; · &ohgr;&tgr;]<RAD><RCD>1+<UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></RCD></RAD></C></R><R><C>−Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>[2(Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>+2)+(Y<SUB><UP>DC</UP></SUB>−1)(Y<SUB><UP>DC</UP></SUB>+2)<UP>tan</UP> &Dgr;&PHgr; · &ohgr;&tgr;]<UP>tan</UP> &Dgr;&PHgr; · &ohgr;&tgr;+(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−2Y<SUB><UP>DC</UP></SUB>)(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>+Y<SUB><UP>DC</UP></SUB>)<UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></C></R></AR></FENCE></NU><DE>(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−4)(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>)−Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>[2(Y<SUP><UP>2</UP></SUP><SUB><UP>DC</UP></SUB>−4)+(Y<SUB><UP>DC</UP></SUB>+2)<SUP>2</SUP> <UP>tan</UP> &Dgr;&PHgr; · &ohgr;&tgr;]<UP>tan</UP> &Dgr;&PHgr; · &ohgr;&tgr;+(Y<SUP><UP>2</UP></SUP><SUB><UP>AC</UP></SUB>−2Y<SUB><UP>DC</UP></SUB>)<SUP>2</SUP> <UP>tan</UP> &Dgr;&PHgr;<SUP>2</SUP></DE></FR>. (27)

Note that r<UP><SUB>o</SUB><SUP>comp</SUP></UP> does not depend explicitly on either the fluorescence lifetime or the modulation frequency; an implicit dependence is via Delta Phi , YAC, and YDC (Eqs. 7-9). This parameter may well be invariant in a given system under investigation, suggesting that a global analysis over an entire image would provide a constant mean value for use in Eq. 20. This procedure also minimizes quantitative discrepancies by ensuring that the same optical configuration applies to all parameters. Others have noted that distortions to the anisotropy decay induced by high NA objectives are mainly manifested in changes in the ro and not phi  (Tramier et al., 2000). Thus, independent self-consistent measures of ro and phi  are desirable.

In this study, the above formalism was applied to images acquired in the fluorescence microscope and with pixel-by-pixel resolution. The apparent correlation times computed according to the above formalism were used to characterize the system and detect heterogeneity with respect to the correlation time and degree of hindered rotation. In the presence of heterogeneity, including that of the intensity decay function (multiple lifetimes), one must resort to more complex models and multi-frequency measurements, as has been described for FLIM (Squire et al., 2000). If some, but not all, parameters are constant over the sample, a global analysis can also be implemented, as in the case of single-frequency FLIM determinations (Verveer et al., 2000). It is evident that the influence of propagated measurement error will vary among the various approaches represented in Eqs. 14-27.

The above equations are valid provided a unique or appropriate mean fluorescence lifetime can be defined. Fluorophores with monoexponential decays yield equivalent phase and modulation lifetimes. This condition does not hold for more complex decay functions, requiring a more elaborate formalism for the selection of an appropriate lifetime average (e.g., a combination of phase and modulation values) for use in the above equations. However, the latter approach fails if molecular association(s) leads to a coupling of multicomponent intensity and anisotropy decay processes (Szmacinski et al., 1987; Jameson and Sawyer, 1995).

Concentration depolarization due to homotransfer Fluorescence Resonance Energy Transfer (emFRET)

In addition to depolarization by rotational diffusion, the emission anisotropy decreases in the event that excitation energy is transferred between nearby molecules during the excited state lifetime. Concentration depolarization due to FRET between identical molecules is a well-documented phenomenon (Bojarski and Sienicki, 1989) that occurs in concentrated solutions of fluorophores exhibiting a finite excitation-emission overlap integral, i.e., a relatively small Stokes shift. In cell biological studies, homotransfer FRET (which we designate here as emFRET) has been applied to the analysis of protein oligomerization on the plasma membrane (Varma and Mayor, 1998; Blackman et al., 1998; Bene et al., 2000) and within the cell (Gautier et al., 2001). In the latter case, emFRET of EGFP fused to viral thymidine kinase was measured by time-correlated single photon anisotropy decay.

We sought to implement emFRET by the wide-field rFLIM imaging technique. To do so, we extended the phase-modulation frequency domain formalism presented above to include provision for the excited state process (Lakowicz and Balter, 1982) of energy migration. For the specific experimental case treated here, that of emFRET involving EGFP in free solution, the extent of rotational depolarization during the excited state lifetime is limited (sigma  < 1) and a good approximation for the apparent composite anisotropy decay (diffusion + energy migration) is given by the product of the terms representing the two processes (see Fig. 4 of Engstrom et al., 1992). In this case, the impulse response functions (Eq. 1) corresponding to the parallel and perpendicular emission components and for the case of an unhindered spherical rotator are given by
I<SUB>∥,⊥</SUB>(t)=<FR><NU>e<SUP><UP>−t/&tgr;</UP></SUP></NU><DE>3&tgr;</DE></FR> [1+&ggr;<SUB>∥,⊥</SUB>r<SUB><UP>o</UP></SUB>e<SUP><UP>−t/&phgr;</UP></SUP>e<SUP><UP>−&agr;</UP><IT>c</IT>(<UP>t/&tgr;</UP>)<SUP><UP>1/2</UP></SUP></SUP>], (28)
where gamma || = 2 and gamma perp  = -1, c is the fluorophore concentration (in mM units), and alpha , the coefficient of the c(t/tau )1/2 term accounting for energy migration, is given (in mM-1 units) by
&agr;=2 · 10<SUP>−27</SUP>&pgr;<RAD><RCD>&pgr;/3</RCD></RAD>N<SUB><UP>av</UP></SUB>⟨‖&kgr;‖⟩R<SUP><UP>3</UP></SUP><SUB><UP>o</UP></SUB>≈(R<SUP><UP>3</UP></SUP><SUB><UP>o</UP></SUB>/375), (29)
where Ro is the characteristic Förster transfer distance (50% transfer efficiency for a single donor-acceptor FRET pair; units, nm), Nav is Avogadro's constant, and < |kappa |> (= 0.69) is the orientation-averaged square root of the Förster orientation factor kappa 2. This formulation assumes a random initial static distribution of molecular transition moments and a completely depolarized emission from any but the initially excited molecule. The ensemble excited state lifetime remains unaltered in emFRET.

The steady-state anisotropy is obtained by integration of Eq. 28, yielding
<A><AC>r</AC><AC>&cjs1171;</AC></A>=<FR><NU>r<SUB><UP>o</UP></SUB></NU><DE>1+&sfgr;</DE></FR> [1−&bgr;e<SUP>&bgr;<SUP>2</SUP></SUP>&pgr;<SUP>1/2</SUP><UP>erfc</UP>(&bgr;)], (30)
where
&bgr;=<FR><NU>&agr;c</NU><DE>2<RAD><RCD>1+&sfgr;</RCD></RAD></DE></FR>. (31)
The convolution of Eq. 28 with a sinusoidal modulating function leads to analytical equations for the difference phase and modulation corresponding to Eqs. 7-9 but including the effect of energy migration, represented by the parameter alpha . A detailed analysis will be given elsewhere. However, the effect of emFRET on Delta Phi and YAC can be summarized as follows: Delta Phi increases and then decreases with c, particularly for high modulation frequencies and low values of sigma ; and YAC decreases at all frequencies.

Determination of fluorescence anisotropy in the microscope

Conventional descriptions of fluorescence anisotropy apply most readily to standard solution measurements of fluorophores with random orientations in the ground state and detection of the emission with the conventional 90° narrow-aperture configuration, i.e., orthogonal to the propagation direction of the excitation beam. In the microscope, one has to deal with the particular excitation/detection configurations including epi-illumination, high numerical apertures, nonrandom orientations in the ground state (Axelrod, 1979; Florine-Casteel, 1990; Fushimi et al., 1990), polarization distortions produced by biases or birefringence in the optical components (i.e., objectives, apertures, reflectors, filters, prisms: Axelrod, 1979; Florine-Casteel, 1990; Dix and Verkman, 1990; Bahlmann and Hell, 2000; Tramier et al., 2000) and detectors, and photobleaching under conditions of high irradiance.

The perturbation by high numerical aperture (NA) illumination/collection optics of polarization states characteristic of plane waves is manifested by mixing of polarization components along the three optical axes (Jovin, 1979; Axelrod, 1979, 1989; Florine-Casteel, 1990; Dix and Verkman, 1990; Sheppard and Torok, 1997; Bahlmann and Hell, 2000; Tramier et al., 2000). According to Axelrod (1989), the contribution of the observed polarization components of the emission can be formulated as linear sums of products of axial distributions dependent on sample properties and weighting factors that are functions of the NA. We have examined these effects empirically by assessing the influence of NA on the steady-state anisotropies and phase-shifts of known samples.


    EXPERIMENTAL RESULTS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL
RESULTS
EXPERIMENTAL RESULTS
DISCUSSION
APPENDIX
REFERENCES

Microscope validation

Steady-state anisotropies of fluorescein solutions of varying viscosity (achieved by addition of glycerol) were measured in the imaging microscope and compared to the corresponding values obtained in a spectroflurometer. In general, the two sets of measurement were in good agreement and could be accounted for by a single G-factor (Eq. 4). Tests for a systematic depolarization bias due to wide-aperture excitation-detection via the objective lens were also carried out. Using solutions of fluorescein in 80% glycerol, a small depolarizing effect, manifested as slight decreases in YAC and <A><AC>r</AC><AC>&cjs1171;</AC></A>, was observed as the magnification and NA of the objective lens increased (Table 1). Similar results have been reported by others (Dix and Verkman, 1990; Verkman et al., 1990; Tramier et al., 2000). According to the theoretical results of Axelrod (Fig. 2 in Axelrod, 1989), the steady-state anisotropy corresponding to an isotropic rotator with ro = 0.4 and sigma  = 0.6 should decrease from 0.25 to ~0.21 as the NA/eta value varies from 0.06 to 0.88. The variation documented in Table 1 was somewhat smaller, probably due to a contribution from out-of-focus light from the relatively thick solutions. Recent rFLIM data from cultured cells expressing GFP-tagged receptors and calibrated microspheres have revealed more pronounced NA effects (Subramaniam et al., 2002).


                              
View this table:
[in this window]
[in a new window]
 
TABLE 1   Effect of objective lens numerical aperture on rFLIM parameters (Delta Phi , YAC, <A><AC>r</AC><AC>&cjs1171;</AC></A>) for fluorescein in 80% (w/w) glycerol/buffer

It is interesting that the corresponding values of Delta Phi in Table 1 were independent of the aperture, in accordance with the observations of Verkman et al. (1990). As expected, removal of the excitation polarizer yielded a Delta Phi  = 0 within experimental error, implying the absence of an instrumental phase lag between the parallel and perpendicular components of emission. Thus, no correction for instrumental phase shift was required.

Single-frequency rFLIM measurements of fluorescein-glycerol solutions

As a model system for demonstrating the ability of rFLIM to image regions differing in <A><AC>r</AC><AC>&cjs1171;</AC></A> and phi , we imaged two adjacent cuvettes containing fluorescein solutions of different viscosity (10% glycerol/90% buffer, eta  = 1.2 cP; 61% glycerol/39% buffer eta  = 10.6 cP) and FLIM and rFLIM measurements. Figure 3 illustrates the acquired polarization data (Delta Phi , YAC, and YDC), the fluorescence lifetime values (tau phase and tau mod), and the corresponding derived quantities: the rotational correlation times [(phi Delta Phi , phi YAC, phi YDC; Eqs. 14-19, with ro = 0.35, rinfinity  = 0) and (phi Delta Phi ,<A><AC>r</AC><AC>&cjs1171;</AC></A>, rinfinity ; Eqs. 20-24, with ro = 0.35)] determined at a single optical modulation frequency of 58 MHz. The mean values and standard deviations of the lifetime and anisotropy parameter distributions are collected in Table 2. Standard errors have been included as an indication of the errors in the means. The large number of pixels in the distribution led to reasonably precise estimates of the mean, and to standard errors that were typically two orders of magnitude smaller than the standard deviations. Consideration of the standard deviations and standard errors also allows an estimate of the number of pixels required to discriminate between given differences in rotational parameters. The minimal detectable difference between the means of two distributed parameters x and y is given by |x - <A><AC>y</AC><AC>&cjs1171;</AC></A>| = st<RAD><RCD><IT>(n</IT><SUB>x</SUB><IT> + n</IT><SUB>y</SUB><IT>)/(n</IT><SUB>x</SUB><IT>n</IT><SUB>y</SUB>)</RCD></RAD>, where t is the Student's test, s is the pooled standard deviation and nx, ny are the population sizes. Anisotropy decay parameters computed from the means of the entire distributions (Delta Phi , YAC, and YDC) for each condition are also given in Table 2.



View larger version (59K):
[in this window]
[in a new window]
 
FIGURE 3   Spatially-resolved FLIM and rFLIM of fluorescein solutions containing 10% and 61% glycerol. Solutions in adjacent cuvettes: 10% glycerol, left in images and red in histograms; 61% glycerol, right in images and blue in histograms. Top and bottom rows, histograms corresponding to the images in the two center rows. Symbols, all of the binned data points; solid lines, Gaussian fits. Both original measured signals (Delta Phi , YAC, YDC) and derived anisotropy parameters are shown. The objective was a Plan 2× air (NA 0.06). See text for nomenclature and discussion.


                              
View this table:
[in this window]
[in a new window]
 
TABLE 2   rFLIM parameters for fluorescein in glycerol/buffer*

Several points emerge from the images and histograms displayed in Fig. 3. First, spatial resolution was achieved with both acquired data and derived par