Picosecond resolution of tyrosine fluorescence and anisotropy decays by 2-GHz frequency-domain fluorometry

We extended the technique of frequency-domain fluorometry to an upper frequency limit of 2000 MHz. This was accomplished by using the harmonic content of a laser pulse train (3.76 MHz, 5 ps) from a synchronously pumped and cavity-dumped dye laser. We used a microchannel plate photomultiplier as the detector to obtain the 2-GHz bandwidth. This new instrument was used to examine tyrosine intensity and anisotropy decays from peptides and proteins. These initial data sets demonstrate that triply exponential tyrosine intensity decays are easily recoverable, even if the mean decay time is less than 1 ns. Importantly, the extended frequency range provides good resolution of rapid and/or multiexponential tyrosine anisotropy decays. Correlation times as short as 15 ps have been recovered for indole, with an uncertainty of +/- 3 ps. We recovered a doubly exponential anisotropy decay of oxytoxin (29 and 454 ps), which probably reflects torsional motions of the phenol ring and overall rotational diffusion, respectively. Also, a 40-ps component was found in the anisotropy decay of bovine pancreatic trypsin inhibitor, which may be due to rapid torsional motions of the tyrosine residues and/or energy transfer among these residues. The rapid component has an amplitude of 0.05, which is about 16% of the total anisotropy. The availability of 2-GHz frequency-domain data extends the measurable time scale for fluorescence to overlap with that of molecular dynamics calculations.     

Lifetime distributions and anisotropy decays of indole fluorescence in cyclohexane/ethanol mixtures by frequency-domain fluorometry

We used frequency-domain fluorometry to measure intensity and anisotropy decay of indole fluorescence in cyclohexane/ethanol mixtures at 20 degrees C. In 100% cyclohexane or 100% ethanol the intensity decay of indole appears to be a single exponential with decay times of 7.66 and 4.10 ns, respectively. In cyclohexane containing a small percentage of ethanol (up to 10%), we observed increased heterogeneity in intensity decay, resulting in a 10-fold increase in chi 2R for the single-exponential fit, as compared with the double-exponential model. We obtained comparable or better fits using unimodal Lorentzian and Gaussian lifetime distributions (two floating parameters) than for the two-exponential model (three floating parameters). We believe that the distribution of decay times reflects a range of indole solvation states in the dominately nonpolar solutions. This result suggests that a variety of hydrogen-bonding configurations could be one origin of the distributions of decay times observed for tryptophan emission from proteins. We also measured rotational diffusion of indole in cyclohexane, ethanol and its mixtures at 20 degrees C. The picosecond correlation times required that the mean decay times be decreased by acrylamide quenching (in ethanol) or energy transfer (in cyclohexane). In ethanol we observed nearly isotropic rotation of indole; in cyclohexane we obtained two correlation times of 17 and 73 ps. The shorter correlation time in cyclohexane appears to be due to the slip boundary condition, which was found to be progressively eliminated by small percentages of ethanol. Hence, hydrogen-bonding interactions appear to have a substantial effect on the rotational dynamics of indole.     

Correction for contaminant fluorescence in frequency-domain fluorometry

We describe a general method to correct for contaminant fluorescence when using the technique of frequency-domain fluorometry. The method can be applied regardless of the origin of the background signal, from scattered light, impurity fluorescence, or both. The procedure requires measurement of the frequency-dependent phase and modulation of the background at enough frequencies to approximate the decay law of the background. We also describe a general method to propagate the uncertainties in the measured phase and modulation values into the corrected values. This propagation is necessary to ensure proper weighting of the frequency-dependent data in the least-squares fitting algorithms. The practical usefulness of this correction method is demonstrated using frequency-domain data for one and two component mixtures which were deliberately contaminated with scattered light and/or other fluorophores.     

Demonstration of an associated anisotropy decay by frequency-domain fluorometry

We used frequency-domain fluorometry to demonstrate the presence of an associated decay of fluorescence anisotropy. In such systems the individual correlation times are associated with distinct emitting species, each with its own characteristic lifetime and rotational correlation times. We obtained an associated system using 1-anilino-8-naphthalenesulfonic acid (ANS) in the presence of increasing amounts of apomyoglobin. When both free and apomyoglobin-bound ANS contributed to the emission the differential polarized phase angles become negative at particular frequencies, even though the fundamental anisotropy (r0) is greater than zero. Additionally, the modulated anisotropy decreases at high frequencies. Both observations appear to be the unique consequence of an associated anisotropy decay, and are not possible for a multiexponential anisotropy decay of a single species.     

Intensity and anisotropy decays of the Wye base of yeast tRNA(Phe) as measured by frequency-domain fluorometry

The intensity and anisotropy decays of Wye base fluorescence from yeast tRNA(Phe) were determined by frequency-domain fluorometry. The intensity decay is at least a double exponential in the presence and absence of Mg2+, but the multi-exponential character of the decay is more pronounced in the absence of Mg2+. The anisotropy decay displays components due to overall tRNA rotational diffusion and to local torsional motions. The amplitude of the local motion is decreased 2-fold by the presence of Mg2+. The results are broadly consistent with a more homogeneous environment for the Wye base in the presence of Mg2+.     

Intensity and anisotropy decays of the tyrosine calmodulin proteolytic fragments, as studied by GHz frequency-domain fluorescence

Frequency-domain fluorescence measurements to 2 GHz were able to recover and account for essentially all of the intrinsic tyrosine anisotropy of calmodulin and its proteolytic fragments containing one or two tyrosine residues. Low-temperature measurements have detected a very rapid initial anisotropy decay in the 2-tyrosine species which may be attributed to radiationless energy transfer between the two tyrosines. The observed values of the rotational correlation times indicate that both tyrosines of calmodulin possess considerable mobility, which decreases in the presence of Ca2+ and at low temperatures.     

Picosecond resolution of oxytocin tyrosyl fluorescence by 2 GHz frequency-domain fluorometry

The technique of frequency-domain fluorometry has been extended to 2000 MHz using the harmonic content of a picosecond laser source and a microchannel plate photomultiplier tube. This new instrument was used to resolve complex subnanosecond intensity and anisotropy decays of the tyrosyl emission of oxytocin. The intensity decay was found to contain at least three exponential components, 80, 359 and 927 ps. The anisotropy analysis revealed a 29 ps torsional motion of the tyrosine residue as well as a 454 ps overall rotational correlation time. The time resolution of this method should permit the comparison of experimental results with theoretical models for motions of proteins.     

Analysis of time-resolved fluorescence anisotropy decays.

We discuss the analysis of time-correlated single photon counting measurements of fluorescence anisotropy. Particular attention was paid to the statistical properties of the data. The methods used previously to analyze these experiments were examined and a new method was proposed in which parallel- and perpendicular-polarized fluorescence curves were fit simultaneously. The new method takes full advantage of the statistical properties of the measured curves; and, in some cases, it is shown to be more sensitive than other methods to systematic errors present in the data. Examples were presented using experimental and simulated data. The influence of fitting range on extracted parameters and statistical criteria for evaluating the quality of fits are also discussed.     

Anisotropy decays of single tryptophan proteins measured by GHz frequency-domain fluorometry with collisional quenching

We used harmonic-content frequency-domain fluorometry to determine the anisotropy decays of a variety of single tryptophan peptides and proteins. Resolution of the rapid and complex anisotropy decays was enhanced by global analysis of the data measured in the presence of quenching by either oxygen or acrylamide. For each protein, and for each quencher, data were obtained at four to six quencher concentrations, and the data analyzed globally to recover the anisotropy decay. The decrease in decay times produced by quenching allows measurements to an upper frequency limit of 2 GHz. The chosen proteins provided a range of exposures of the tryptophan residues to the aqueous phase, these being ACTH, monellin, Staphylococcus nuclease and ribonuclease T ₁, in order of decreasing exposure. Examination of indole and several small peptides demonstrates the resolution limitations of the measurements; a correlation time of 12 ps was measured for indole in methanol at 40°C. Comparison of the anisotropy decays of gly-trp-gly with leu-trp-leu revealed stearic effects of the larger leucine side chains on the indole ring. The anisotropy decay of gly-trp-gly revealed a 40 ps component for the indole side chain, which was resolved from the overall 150 ps correlation time of the tripeptide. Only the longer correlation time was observed for leu-trp-leu. With the exception of ribonuclease T ₁, each of the proteins displayed a subnanosecond component in the anisotropy decay which we assign to independent motions of the tryptophan residues. For example, Staphylococcus nuclease and monellin displayed segmental tryptophan motions with correlation times of 80 and 275 ps, respectively. The amplitudes of the rapid components increased with increasing exposure to the aqueous phase. These highly resolved anisotropy decays for proteins of known structure are suitable for comparison with molecular dynamic simulations.Abbreviations Ac acrylamide - ACTH adrenocorticotropin hormone (1–24) - BPTI bovine pancreatic trypsin inhibitor - NATA N-acetyl-L-tryptophanamide - RNase T ₁ ribonuclease T ₁ - S. Nuclease staphylococcus aureus nuclease Supported by grants DMB-8804931 and DIR-8710401 from the National Science Foundation, and GM-39617 from the National Institutes of Health. J. R. Lakowicz acknowledges support from the Medical Biotechnology Center at the University of Maryland. I. Gryczynski was on leave from University of Gáansk, Institute of Experimental Physics, Gdansk, Poland, with partial support from CPBP 01.06.2.01 (Poland). H. Cherek was on leave from Nicholas Copernicus University, Torun, Poland, with partial support from CPBP 01.06.2.03 Offprint requests to: J. R. Lakowicz     

Distribution of distances in thiopeptides by fluorescence energy transfer and frequency-domain fluorometry

Frequency-domain fluorescence spectroscopy was employed to examine the decays of tryptophan in Boc-Trp-Met-Asp-Phe-NH2 (donor) and (Formula: see text) (donor-acceptor pair). The efficiency of energy transfer in the thiopeptide amounted to 60%. The measured dispersion of fluorescence decay times was used to recover the donor-acceptor distance distribution. The parameters of the Gaussian distance distribution obtained for this peptide (r, the mean distance (9 A); hw, the halfwidth (25 A)) indicate the lack of a distinct favorable conformation.     

Determination of time-resolved fluorescence emission spectra and anisotropies of a fluorophore-protein complex using frequency-domain phase-modulation fluorometry

We report the first time-resolved fluorescence emission spectra and time-resolved fluorescence anisotropies obtained using frequency-domain fluorescence spectroscopy. We examined the fluorophore p-2-toluidinyl-6-naphthalenesulfonic acid (TNS) in viscous solvents and bound to the heme site of apomyoglobin using multifrequency phase fluorometers. Fluorescence phase shift and modulation data were obtained at modulation frequencies ranging from 1 to 200 MHz. For time-resolved emission spectra, the impulse response for the decay of intensity at each emission wavelength was obtained from the frequency response of the sample at the same emission wavelength. The decays have negative pre-exponential factors, consistent with a time-dependent spectral shift to longer wavelengths. These multiexponential decays were used to construct the time-resolved emission spectra, which were found to be in good agreement with earlier spectra obtained from time-domain measurements. Additionally, time-resolved anisotropies were obtained from the frequency-dependent phase angle differences between the parallel and perpendicularly polarized components of the emission. The rotational correlation times of TNS bound to apomyoglobin are consistent with those expected for this probe rigidly bound to the protein. TNS in propylene glycol also displayed a single exponential decay of anisotropy. These results, in conjunction with the previous successful resolution of multiexponential decays of fluorescence intensity (Lakowicz, J. R., Gratton, E., Laczko, G., Cherek, H., and Limkeman, M. (1984) Biophys. J., in press; Gratton, E., Lakowicz, J. R., Maliwal, B. P., Cherek, H., Laczko, G., and Limkeman, M. (1984) Biophys. J., in press) demonstrate that frequency-domain measurements provide information which is, at a minimum, equivalent to that obtainable from time-domain measurements.     

Resolution of complex anisotropy decays by variable frequency phase-modulation fluorometry: a stimulation study

We used simulations to determine the resolution of complex anisotropy decay laws which is obtainable by frequency-domain fluorometry. The simulations include the effects of torsional and segmental motions of tryptophan residues in proteins, the multiple correlation times of asymmetric molecules, and three-component anisotropy decays. For a protein with a global correlation time of 10 ns it should be possible to resolve torsional motions with correlation times as short as 10 ps if the amplitude of the rapid motion is at least 20% of the total anisotropy decay with r0 = 0.4. Correlation times which differ by only 1.4-fold can be resolved, making this method useful for determination of the shape of proteins and other asymmetric molecules. It is possible to resolve three-component anisotropy decays if the overall difference among the correlation times is 30-fold. Such resolution will be useful for understanding of internal motions of proteins and membranes. The validity of these predictions is demonstrated in the subsequent paper using experimental data for melittin in solution and when bound to membranes (Maliwal, B.P., Hermetter, A. and Lakowicz, J.R. (1986) Biochim. Biophys. Acta 873, 173-181).     

Anisotropy decays of indole, melittin monomer and melittin tetramer by frequency-domain fluorometry and multi-wavelength global analysis

We used frequency-domain fluorescence spectroscopy to measure the fluorescence lifetime and anisotropy decays of indole in propylene glycol, and of the tryptophan emission of melittin monomer and tetramer in water solutions at 5 degrees C. We obtained an increase in resolution of the anisotropy decays by using multiple excitation wavelengths, chosen to provide a range of fundamental anisotropy values. The multi-excitation wavelength anisotropy decays were analyzed globally to recover a single set of correlation times with wavelength-dependent anisotropy amplitudes. Simulated data and kappaR2 surfaces are shown to reveal the effect of multi-wavelength data on the resolution of complex anisotropy decays. For both indole and melittin, the anisotropy decays are heterogeneous and require two correlation times to fit the frequency-domain data. For indole in propylene glycol at 5 degrees C we recovered correlation times of 0.59 and 4.10 ns, which appear to be characteristic of the rigid and asymmetric indole molecule. For melittin monomer the correlation times were 0.13 and 1.75 ns, and for melittin tetramer 0.12 and 3.96 ns. The shorter and longer correlation times of melittin are due to segmental motions and overall rotational diffusion of the polypeptide.     

A study of protein dynamics from anisotropy decays obtained by variable frequency phase-modulation fluorometry: internal motions of N-methylanthraniloyl melittin

Internal motions of melittin and its lipid complexes were studied by anisotropy decays determined by frequency-domain fluorometry. A covalent anthraniloyl probe was attached, probably to lysine-21. The emission spectra indicate that the anthraniloyl moiety is exposed to solvent in both monomeric and tetrameric forms and is present at the lipid-water interfacial region in the lipid complexes. The fluorescence intensity decay of melittin in solution and its lipid complexes was characterized by three lifetimes. The lifetimes were near 1-2 ns, 6-7 ns and 10 ns. At increased temperatures there was an increase in the amplitude of the intermediate lifetime and a decrease in that of the longer lifetime. For all the melittin systems, at least three correlation times were required to fit the anisotropy data. Of the three correlation times, the shortest correlation time represents the local motions of the probe, while the longest represents global motions of the whole molecule. The intermediate correlation time probably represents the dynamics of domains/helices within the molecule. The melittin monomer is highly flexible, with greater than 90% of its anisotropy being lost by the local motions. Even though it is well organized (greater than 75% helical), the tetramer is still a highly flexible molecule, with 70% of its anisotropy being lost by the local motions. The internal motions of melittin decrease upon binding to lipids and are sensitive to the phase state of the lipid complexes.     

Steady-state fluorescence anisotropy and multifrequency phase fluorometry on oxidized phosphatidylcholine vesicles

Multilamellar liposomes, from mixtures of unoxidized (control) and singlet oxygen oxidized phosphatidylcholine, were studied by steady-state fluorescence anisotropy and multifrequency phase fluorometry using 1,6-diphenyl-1,3,5-hexatriene (DPH) as fluorescent probe. Lifetime fluorescence decay of the DPH-labeled liposomes was analyzed either by a model of discrete exponential components and a model that assumes a continuous distribution of lifetime values. Increasing the oxidized phosphatidylcholine content in the liposomes, an increase of the membrane interior polarity and a decrease of membrane fluidity occurs which can be related to the hydroperoxide-lipids and double bonds conjugation, respectively.     

Analysis of fluorescence anisotropy decays by a least square method

A method of fluorescence anisotropy decay analysis is described in this work. The transient anisotropy r(ex)(t) measured in a photocounting pulsefluorimeter is fitted by a non linear least square procedure to the ratio of convolutions of the apparatus response function g(t) by sums of appropriate exponential functions. This method takes rigorously into account the apparatus response function and is applicable to any shape of the later as well as to any values of fluorescence decay times and correlation times. The performances of the method have been tested with data simulated from measured response functions corresponding to an air lamp and a high pressure nitrogen lamp. The statistical standard errors of the anisotropy deca parameters have been found to be smaller than the standard errors previously calculated for the moment method. A systematic error delta in the fluorescence decay time entailed an error deltatheta in the correlation time such as Deltatheta/theta < deltatau/tau. By this method, good fitting of experimental data have been achieved very conveniently and accurately.     

Nanosecond fluorescence anisotropy decays of 1,6-diphenyl-1,3,5-hexatriene in membranes

Nanosecond decays of the fluorescence anisotropy, $\text{r}$, were studied for the emission of 1,6-diphenyl-l,3,5-hexatriene (DPH) embedded in a series of mixed multilamellar liposomes containing egg yolk phosphatidylcholine, phosphatidylethanolamine and cholesterol in varying molar ratios, as well as in membranes of intact cells and in virus envelopes.The relative contributions of the fast and the infinitely slow decaying component to the steady-state value, $\text{r}$, of the fluorescence anisotropy were very similar for artifical and biological membranes.Angles, θ, of the cone, by which the motion of the fluorescent molecule is limited, were calculated from the intensity of the infinitely slow decaying anisotropy component and compared with steady-state fluorescence anisotropies and with ‘microviscosities’, 〈η〉. An increase in 〈η〉 from 1.5 to 5.2 P in our systems was accompanied by a decrease in θ from 49° to 30° while the decrease in the mean motional relaxation times, $\text{φ}{}_{\mathrm{<mtext>f</mtext>}}$, of the label molecule was not more than 1 ns and due mainly to changes in the potential, by which the diffusion of DPH in the membrane is restricted. From these observations we conclude that differences in the steady-state fluorescence anisotropy and in ‘microviscosities’ of cholesterol-containing membranes ($\text{r}\text{> 0.15}$) represent changes in the degree of static orientational constraint rather than changes in diffusion rates of the label.     

Resolution of the lifetimes and correlation times of the intrinsic tryptophan fluorescence of human hemoglobin solutions using 2 GHz frequency-domain fluorometry

We used 2 GHz harmonic content frequency-domain fluorescence to measure the intensity and the anisotropy decays from the intrinsic tryptophan fluorescence from human hemoglobin (Hb). The tryptophan intensity decays are dominated by a short-lived component which accounts for 35-60% of the total steady state intensity. The decay time of this short component varies from 9 to 27 ps and this component is sensitive to the ligation state of Hb. Our error analyses indicate the uncertainty is about +/- 3 ps. The intensity decays also show two longer lived components near 0.7 and 8 ns, which are probably due either to impurities or to Hb molecules in conformations which do not permit energy transfer. The anisotropy decays indicate the tryptophan residues in Hb are highly mobile, with apparent correlation times near 55 ps.     

Complex homogeneous and heterogeneous fluorescence anisotropy decays: enhancing analysis accuracy

In biological macromolecules, fluorophores often exhibit multiple depolarizing motions that require multiple lifetimes and rotational relaxation times to define fluorescence intensity and anisotropy decays. The related analysis of time-correlated single-photon counting data becomes uncertain due to the multitude of decay parameters and numerical sensitivity to deconvolution of the instrument response function (IRF) via discretization of integrals. By using simulations we show that improved discretizations based on quadratic and cubic local approximations of the IRF yield more accurate estimation of short rotational relaxation times and lifetimes than the commonly used Grinvald-Steinberg discretization, which in turn appears more reliable than two discretizations based on linear local approximations of the IRF. In addition, our simulation suggests that cubic approximation is the most advantageous in discriminating complex heterogeneous and homogeneous anisotropy decay. We show that among three different information criteria, the Akaike information criterion is best suited for detection of heterogeneity in rotational relaxation times. It is capable of detecting heterogeneity even when anisotropy decay appears homogeneous within statistical errors of estimation.     

Resolvability of fluorescence lifetime distributions using phase fluorometry.

The analysis of the fluorescence decay using discrete exponential components assumes that a small number of species is present. In the absence of a definite kinetic model or when a large number of species is present, the exponential analysis underestimates the uncertainty of the recovered lifetime values. A different approach to determine the lifetime of a population of molecules is the use of probability density functions and lifetime distributions. Fluorescence decay data from continuous distributions of exponentially decaying components were generated. Different magnitudes of error were added to the data to simulate experimental conditions. The resolvability of the distributional model was studied by fitting the simulated data to one and two exponentials. The maximum width of symmetric distributions (uniform, gaussian, and lorentzian), which cannot be distinguished from single and double exponential fits for statistical errors of 1 and 0.1%, were determined. The width limits are determined by the statistical error of the data. It is also shown that, in the frequency domain, the discrete exponential analysis does not uniformly weights all the components of a distribution. This systematic error is less important when probability and distribution functions are used to recover the decay. Finally, it is shown that real lifetime distributions can be proved using multimodal probability density functions. In the companion paper that follows we propose a physical approach, which provides lifetime distribution functions for the tryptophan decay in proteins. In the third companion paper (Alcala, J.R., E. Gratton, and F.J. Prendergast, 1987, Biophys. J., in press) we use the distribution functions obtained to fit data from the fluorescence decay of single tryptophan proteins.