| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
odarczykDepartment of Biophysics, Institute of Experimental Physics, University of Warsaw, Warsaw, Poland
Correspondence: Address reprint requests to Borys Kierdaszuk, Fax: 48-22-554-0001; E-mail: borys{at}biogeo.uw.edu.pl.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
It is well-known that decays of fluorescence intensity in proteins often exhibit a complex behavior, the origin of which can result from multiple (ground state) conformations, protein dynamics, spectral relaxation, or other interactions between fluorophores and their environment (Lakowicz, 2000
). The decay of protein fluorescence is usually described by a multiexponential model, i.e., the arithmetic sum of a number of exponents characterized by lifetimes and pre-exponential factors. This is based on the interpretation of each discrete component with the aid of a particular protein conformation, including dynamic equilibrium between rotational isomers of tryptophan (Ross et al., 1992a
) and tyrosine residues (Ross et al., 1992b
, 1986
; Laws et al., 1986
). However, there are many cases for which an individual decay component has no physical interpretation, e.g., in the case of a multitude of conformational substates in proteins with possible time-dependent interconversions between them (Alcala et al., 1987
). Furthermore, fluorophores that undergo dipolar relaxation (Ladokhin, 1999
), as well as chemically heterogeneous systems (Gryczynski et al., 1988
), display a complex fluorescence decay and cannot be properly described using the rotamer hypothesis (Ladokhin and White, 2001
). Even for the indole moiety in solution, it is difficult to find unique conformations corresponding to discrete decay components (Szabo et al., 1983
; Gryczynski et al., 1988
). In such cases application of the multiexponential decay model is often arbitrary, and introduction of additional exponential components (few extra parameters) may improve fitting, but often without physical significance. An alternative approach suggests that, for a complex heterogeneous system, continuous lifetime distributions would be more relevant than a sum of discrete terms. We have to consider a better model to describe complex fluorescence relaxation linking the character of the decay law to the distribution of lifetimes.
| THEORY |
|---|
|
|
|---|
![]() | (1) |
(t) is the decay time of the system.
Solution of Eq. 1 in the case of
(t) =
0 = const leads to the classical single-exponential decay function
![]() | (2) |
Traditional analyses of a system of emitting residues treat the interaction with the environment as a small perturbation, and allow description of the excited-state decay of the system by its eigenvalues, i.e., discrete quantities (amplitudes and fluorescence decay lifetimes). These analyses have also shown that, with the assumption of noninteracting resonance transitions (separated by large energy gaps compared to their typical width), each individual fluorophore exhibits an exponential decay behavior. Obviously for a mixture of isolated and noninteracting residues, the decay of I(t) is given by the multiexponential model
![]() | (3) |
i are the amplitudes and decay times of the M exponential components of the fluorescence decay. The mean decay time is given by 
e
=
Ai
i.
On the other hand, decay of excited states may be treated as a result of interaction of emitting residues with their environment. Following this, one should consider many intercoupled states with coupling between them described by the coupling matrix Wij, the elements of which were a collection of transition rates from the ith energy level to the jth one. Such a mixing of configurations is present in proteins, where the energy of the excited state is dissipated among many conformational substates. It means that the component representing a transition between the excited and ground states of the emitting residues, in reference to the whole protein, is decomposed into a large number of resonances corresponding to different conformational substates of the protein in the ground and excited states. Previous investigations (Alcala et al., 1987
; Alcala, 1994
) indicate that proteins can exhibit a distribution of conformations. It is also known that, at room temperature, the timescale of a transition (interconversion) between different energy levels corresponding to different conformations can be of the same order of magnitude as the excited state decay time (Alcala et al., 1987
).
In such instances the conformational equilibria of amino acids in solution or in proteins, and overall protein conformations, as well as numerous quenching phenomena, point to a lifetime distribution, P(
), which better describes the contribution of each component to the fluorescence decay than discrete values of amplitudes Ai for each decay component. Consequently, the fluorescence decay I(t) is determined by
![]() | (4) |
The lifetime distribution approach in fluorescence decay analysis is usually used without a theoretical basis for the lifetime distribution, which is arbitrarily considered as Gaussian or Lorentzian. In this approach a part of the distribution P(
) exists for
< 0 (cutoff problem), and, therefore, needs additional normalization (Lakowicz, 1999
). The cutoff problem is avoided with the geometric distribution function (e.g., top-hat distributed decay times). In this case, P(
) goes to zero in a well-defined manner and is given by
![]() | (5) |
is full width of the top-hat function.
In an alternative approach, P(
) is not described by any particular function. Thus, one can ask about the most expected distribution under certain constraints. Such a distribution can be found using the maximum entropy method (Montroll and Schlesinger, 1983
; Sobczyk and Tr
bicki, 1993
). In the maximum entropy formalism, one seeks the distribution function P(
= 1/
) that maximizes the entropy
![]() | (6) |
), which is subject to the constraints of normalization (with m0 = 1 for k = 0), and L data constraints, given by
![]() | (7) |
![]() | (8) |
k) are determined by L + 1 equations given by Eq. 7. Thus, we have the constraints on normalization (
), the mean value of the decay rate specified as 

, and the mean value of the logarithm
ln(
)
. The latter expresses the fact that the distribution is determined only for
> 0. Then, taking into account that P(
)d
= P(
)d
, and
ln(
)
-ln(N
0), i.e., the case of a complex system with N decay paths affecting the excited-state lifetime, the distribution of lifetimes P(
), which maximizes the entropy function with the above constraints, is given in the form of the gamma distribution
![]() | (9) |
) can also be obtained in analogy to the case of one decay path, where P(
) is known as a Porter-Thomas distribution, which is, in fact, a gamma distribution with 1 degree of freedom (N = 1; see Dittes, 2000
As shown above, the distribution P(
) in gamma form is predicted mathematically, and is proven not only by a maximum-entropy method but also by a random-matrix-theory approach, describing mixing between allowed degrees of freedom (configurations). In the formulation of random matrix theory, the coupling matrix (Wij) elements, describing the system and its environment, follow some random probability distribution. Only recently has it been proven that, as a result of the composition of contributions from N decay paths, the distribution of eigenvalues of Wij is described by P(
) in the form of Eq. 9 (Wilk and W
odarczyk, 2001
), where
-values are equal to reciprocal eigenvalues.
Furthermore, introduction of the distribution P(
) in the gamma form of Eq. 9 into Eq. 4 leads directly to the decay function
![]() | (10) |
0) and one new parameter of heterogeneity (q), which is defined as
![]() | (11) |
) of fluctuations of
= 1/
around the 1/
0 value (Wilk and W
odarczyk, 1999
o in Eq. 10 results from normalization, and can be substituted by a fitted parameter in the case of non-normalized decay data. The mean decay time
tp
is obtained directly from the integration of Eq. 10, and is given by
![]() | (12) |
In the classical limit, when the number of decay channels goes to infinity, i.e., q
1, the gamma distribution (Eq. 9) becomes a delta function
(
-
0), and the decay function converges from a power-like form (Eq. 10) to a single-exponential form.
Note that the power-like decay function results directly from fluctuations of the parameter
= 1/
in the exponential formula. Inasmuch as protein structural fluctuations occur in the nanosecond-picosecond timescale (Careri et al., 1979
; Karplus and McCammon, 1983
), it follows that these rapid fluctuations would affect fluorescence lifetime values.
It is worth noting that the stretched exponential function,
![]() | (13) |
![]() | (14) |
![]() | (15) |
(t = 0) =
0/q, which means that that system instantly (at the moment of excitation) recognizes that the decay will be characterized by a mean lifetime
0. However, Eq. 14 for t = 0 leads to
(t = 0) = 0, and decay rate 1/
goes to infinity. | EXPERIMENTAL |
|---|
|
|
|---|
Ultraviolet absorption was monitored with a Cary 50 recording instrument (Varian, Palo Alto, CA) fitted with a thermostatically controlled cell compartment, using 5-mm pathlength cuvettes.
Steady-state fluorescence emission and excitation spectra were measured with a FluoroMax spectrofluorimeter (Spex, Metuchen, NJ), with 2-nm spectral resolution for excitation and emission in 5 x 5-mm Suprasil cuvettes, as previously described (Kierdaszuk et al., 2000
).
Measurements of pH (±0.05) were carried out with a CP315m (Elmetron, Zabrze, Poland) pH-meter equipped with a combination semimicro electrode (Orion, Whetstone, Leicester, UK) and temperature sensor.
Time-resolved fluorescence measurements
Fluorescence intensity decays were obtained by time-correlated single photon-counting (measurements of fluorescence, performed on an IBH time-resolved spectrofluorimeter (IBH Consultants, Glasgow, UK) equipped with a laser system consisting of a Ti-Sapphire laser (Mira-900) pumped by an Inova 310 argon ion laser (Coherent, Santa Clara, CA), and two excitation and emission monochromators operated at 4-nm spectral resolution (see Kierdaszuk et al., 2000
, for further details). Typical fluorescence data were convoluted with the instrument response, and fitted to sums of one, two, or three exponentials (Eq. 3), and power-like decay function (Eq. 10) using IBH software and homemade procedures under Matlab, respectively. Two nonlinear fitting procedures, in the form of the Nelder-Meade simplex algorithm (Nelder and Meade, 1965
) or the Marquardt-Levenberg method (Marquardt, 1963
; Levenberg, 1944
), were applied, and led to the same results. The quality of fits was evaluated by the structure observed in the plots of residuals (Figs. 1, 2, and 3, bottom panels) and by the reduced chi-square values (Table 1). All samples were prepared in 50 mM HEPES buffer (pH = 7) at 25°C, unless otherwise stated.
|
|
|
|
| RESULTS |
|---|
|
|
|---|
At neutral pH values, where L-Tyr exists in the zwitterionic form with the phenol residue in the neutral form, its fluorescence decay is well-characterized by a single-exponential function (Laws et al., 1986
). Thus, the typical fluorescence decay of L-tyrosine was analyzed by single-exponential and power-like decay functions (Table 1), with corresponding residuals shown in Fig. 1 (bottom panel). It should be noted that the power-like model provides a good fit to the intensity decay data, with heterogeneity parameter q
1, and with the same mean decay time (
tp
= 3.23 ± 0.06 ns) as that obtained by fit of the single-exponential function (
e
= 3.21 ± 0.02 ns). However, the decay function given by Eq. 10 provides a little better fit, as compared to that obtained for the single-exponential model (Table 1, Fig. 1). These results confirm that the decay function in form of Eq. 10 converge to exponent, when number of decay paths N
(thus q
1), and may fit single-exponential decays.
The acetyl-amide form of L-Tyr (NATyrA), which exists exclusively in the neutral form in aqueous solution at neutral pH, is a better model compound for phenolic residues in a polypeptide chain than parent L-Tyr. Unexpectedly, NATyrA exhibited a more complex fluorescence decay than L-Tyr, usually described by a sum of two exponential terms (Ross et al., 1992a
,b
). It was suggested that this may reflect a conformational equilibrium between two main rotamers identified by 1H NMR spectroscopy (Laws et al., 1986
), but the reason for absence of this effect in the case of L-Tyr fluorescence is so far unknown. Our results for NATyrA fluorescence decay at neutral pH at 25°C (Fig. 2, Table 1) are in line with these observations. A bad fit was obtained for the single-exponential model, as judged by residuals (
)much better, but not fully satisfactory, with the double-exponential model (
). In contrast, an excellent fit was obtained with the power-like decay function (
) and q = 1.07, which led to the number of decay paths (N = 2/(q - 1)), N
29, and may reflect free fluorophore in solution. Although applicability of the multiexponential models is not physically justified, the resultant mean decay time (
e
= 1.63 ± 0.05 ns) for NATyrA fluorescence is in agreement with that obtained with the power-like model (
tp
= 1.53 ± 0.08 ns).
The ternary complex of formycin A (inhibitor) and phosphate (natural co-substrate) with E. coli PNP was chosen as a good example of highly complex fluorescence intensity decay resulting from excitation of tyrosine residues in the enzyme, the N(1)-H and N(2)-H tautomeric forms of formycin A, free in solution, and the latter bound in the active site of the enzyme (Kierdaszuk et al., 2000
). Fluorescence decay of such a complex system cannot be properly described using a single-exponential model (Kierdaszuk et al., 2000
) due to many interacting fluorophore residues, which prevent consideration of the individual decay times (Lakowicz, 1999
). In addition, fluorescence resonance energy transfer between protein tyrosine residues and the N(2)-H form bound by the enzyme (Kierdaszuk et al., 2000
) affects the decay kinetics. Therefore, a continuous lifetime distribution seems to be better than a model that considers a sum of the individual decay components. Indeed, we observed that the power-like decay function given by Eq. 10 provides a good fit of the data (Fig. 3) with
, and high accuracy of the fitted parameters (Table 1). In contrast, a poor fit was obtained with the single-exponential model (
), somewhat improved with the double-exponential model (
), but the values of residuals (Fig. 3, bottom panel) confirm that the power-like decay function is still better than the double-exponential function. Furthermore, the increasing number of exponential terms in the multiexponential model is justified only by the resultant improvement of the goodness of fit, without physical interpretation of the decay components. The value of the mean decay time from the multiexponential model (Kierdaszuk et al., 2000
) roughly agrees with that obtained from Eq. 12 (Table 1).
Tryptophan fluorescence decays were measured for L-Trp (Fig. 4), NATrpA (Fig. 5) and horse liver alcohol dehydrogenase (LADH) (Fig. 6). In contrast to L-Tyr and NATyrA (see above), fluorescence decay of the zwitterionic form of L-Trp in neutral aqueous solution (Fig. 4) was best fitted with the biexponential model (
e
= 2.88 ± 0.05 ns), while neutral NATrpA (Fig. 5) with the monoexponential model (
e
= 2.72 ± 0.01 ns) (Table 1), similar to the results reported previously (Lakowicz, 2000
; and references cited therein). In both cases, the fits of the power-like model were better, as judged on
-values, and led to the q-value, much higher in the case of biexponential than monoexponential decay, as in the case of NATyrA and L-Tyr (Table 1), respectively. It means that the number of decay paths (N) is one range-of-value lower for biexponential decays (
27) than for monoexponential decays, where it is
200, and reflects the fact the power-like function is transformed into monoexponential function, when q
1 (N
).
|
|
|
200 to 20. Additionally, stretched the exponential model and the top-hat distribution function were also fitted to the fluorescence decay of NATrpA in dioxan (Fig. 7). Both models led to the similar fluorescence decay meantime values, i.e., 4.36 ± 0.04 ns (ß = 0.05 ± 0.02), and 4.44 ± 0.01 ns (
= 0.003 ± 0.001), respectively. However, the best fit was obtained with the power-like model (
), as judged on
-values for the stretched exponential model (4.40) and the top-hat distribution function (2.32), as well as on the appropriate residuals (Fig. 7).
|
7 Å), the experimental data argue against significant energy homotransfer (Eftink, 1992
-value and residuals, and show fluorescence lifetimes of components (3.37 ± 0.04 and 6.64 ± 0.02 ns) and decay meantime of 5.78 ± 0.09 ns (Table 1), typical for this protein (Lakowicz et al., 1996
30), i.e., in the range of values observed for L-Trp in aqueous medium, and NATrpA in dioxan (Table 1). Additionally, the Anderson-Darling test (Stephens, 1976
Finally, we should emphasize that the power-like decay model provides additional information about the heterogeneity of the studied fluorophore system, which affects its fluorescence decay. This is based on the relative variance of fluctuations (
) around mean values of lifetime distributions, and the number of decay paths (N = 2/(q - 1)), which can be easily obtained from the fitted parameter q (Eq. 11). With the PNP-FA-Pi ternary complex, q = 1.105 (Table 1), leading to
= 0.105 and to a limited number of decay paths (N
19), which may be related to the numerous phenomena. Some of them should be considered as more important, e.g., 1), conformational equilibrium of each emitting residue; 2), stoichiometry of ligand binding to the six active sites of the enzyme; 3), tautomeric equilibrium of FA, which is affected by enzyme-ligand binding; and 4), fluorescence resonance energy transfer from tyrosine residue(s) to the base moiety of FA in the enzyme-ligand complex (Kierdaszuk et al., 2000
).
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This research was supported by the State Committee for Scientific Research (KBN, grant 6P04A03812, and partially from grants BW 1565/BF and BST 763/BF), and by the Foundation for Polish Science.
Submitted on November 15, 2003; accepted for publication March 19, 2003.
| REFERENCES |
|---|
|
|
|---|
Alcala, J. R. 1994. The effect of harmonic conformational trajectories on protein fluorescence and lifetime distributions. J. Chem. Phys. 101:45784584.
Beechem, J. M., J. R. Knutson, J. B. A. Ross, B. V. Turner, and L. Brand. 1988. Global resolution of heterogeneous decay by phase/modulation fluorometry: mixtures and proteins. Biochemistry. 22:60546058.
Careri, G., P. Fasella, and E. Gratton. 1979. Enzyme dynamics: the statistical physics approach. Annu. Rev. Biophys. Bioeng. 8:6997.[Medline]
Dittes, F. M. 2000. The decay of quantum systems with a small number of open channels. Phys. Rep. 339:215216.
Eftink, M. R., and K. A. Hagman. 1986. Fluorescence lifetime and anisotropy studies with liver alcohol dehydrogenase and its complexes. Biochemistry. 25:66316637.[Medline]
Eftink, M. R. 1992. Luminescence studies with horse liver alcohol dehydrogenase: information on the structure, dynamics, transitions and interactions of this enzyme. Adv. Biophys. Chem. 2:81114.
Gryczynski, I., W. Wiczk, M. L. Johnson, and J. R. Lakowicz. 1988. Lifetime distributions and anisotropy decays of indole fluorescence in cyclohexane/ethanol mixtures by frequency-domain fluorometry. Biophys. Chem. 32:173185.[Medline]
Hirshfield, M. S., S. Chaffe, L. Koro-Johnson, A. Mary, A. A. Smith, and S. A. Short. 1991. Use of site-directed mutagenesis to enhance the epitope-shielding effect of covalent modification of proteins with polyethylene glycol. Proc. Natl. Acad. Sci. USA. 88:71857189.
Karplus, M., and J. A. McCammon. 1983. Dynamics of proteins, elements and function. Annu. Rev. Biochem. 53:263300.
Kierdaszuk, B., A. Modrak-Wójcik, J. Wierzchowski, and D. Shugar. 2000. Formycin A and its n-methyl analogues, specific inhibitor of E. coli purine nucleoside phosphorylase (PNP): induced tautomeric shifts on binding to enzyme, and enzyme
ligand fluorescence resonance energy transfer. Biochim. Biophys. Acta. 1476:109128.[Medline]
Ladokhin, A. S. 1999. Red-edge excitation study of nonexponential fluorescence decay of indole in solution and in a protein. J. Fluorescence. 9:19.
Ladokhin, A. S., and S. H. White. 2001. Alphas and taus of tryptophan fluorescence in membranes. Biophys. J. 81:18251827.
Lakowicz, J. R., B. Kierdaszuk, I. Gryczynski, and H. Malak. 1996. Fluorescence of horse liver alcohol dehydrogenase using one- and two-photon excitation. J. Fluorescence. 6:5159.
Lakowicz, J. R. 1999. Principles of Fluorescence Spectroscopy, 2nd Ed. Kluwer Academic/Plenum Publishers, New York. 130131.
Lakowicz, J. R. 2000. On spectral relaxation in proteins. Photochem. Photobiol. 72:421437.[Medline]
Laws, W. R., A. J. B. Ross, H. R. Wyssbrod, J. M. Beechem, L. Brand, and J. C. Sutherland. 1986. Time-resolved fluorescence and 1H NMR studies of tyrosine and tyrosine analogues: correlation of NMR-determined rotamer populations and fluorescence kinetics. Biochemistry. 25:599607.[Medline]
Lee, B. K. C., J. Siegel, S. E. D. Webb, S. Léveque-Fort, M. J. Cole, R. Jones, K. Dowling, M. J. Lever, and P. M. W. French. 2001. Application of the stretched exponential function to fluorescence lifetime imaging. Biophys. J. 81:12651274.
Levenberg, K. A. 1944. A method for the solution of certain non-linear problems in least squares. Quart. Appl. Math. II:164168.
Mao, C., W. J. Cook, M. Zhou, G. Koszalka, T. A. Krenitsky, and S. E. Ealick. 1997. The crystal structure of E. coli purine nucleoside phosphorylase: a comparison with human enzyme reveals a conserved topology. Structure. 5:13731383.[Medline]
Marquardt, D. W. 1963. An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Ind. Appl. Math. 11:431441.
Montroll, E. W., and M. F. Schlesinger. 1983. Maximum entropy formalism, fractals, scaling phenomena, and 1/f noise: a tale of tails. J. Stat. Phys. 32:209230.
Nelder, J. A., and R. Meade. 1965. A simplex method for function minimization. Comput. J. 7:308313.
Petterson, G. 1987. Liver alcohol dehydrogenase. Crit. Rev. Biochem. 21:349389.[Medline]
Ross, J. B. A., W. R. Laws, A. Buku, J. C. Sutherland, and H. R. Wyssbrod. 1986. Time-resolved fluorescence and 1H NMR studies of tyrosyl residues in oxytocin and small peptides: correlation of NMR-determined conformations of tyrosyl residues and fluorescence decay kinetics. Biochemistry. 25:599607.[Medline]
Ross, J. B. A., H. R. Wyssbrod, R. A. Porter, G. P. Schwartz, C. A. Michaels, and W. R. Laws. 1992a. Correlation of tryptophan intensity decay parameters with 1H NMR-determined rotamer conformations: [tryptophan2]oxytocin. Biochemistry. 31:15851594.[Medline]
Ross, J. B. A., W. R. Laws, K. W. Rousslang, and H. R. Wyssbrod. 1992b. Tyrosine fluorescence and phosphorescence from proteins and polypeptides. In Topics in Fluorescence Spectroscopy. Vol. 3, Biochemical Applications. J. R. Lakowicz, editor. Plenum Press, New York. 163.
Sobczyk, K., and J. Tr
bicki. 1993. Maximum entropy principle and stochastic nonlinear oscillators. Physica A. 193:448468.
Stephens, M. A. 1976. EDF Statistics for goodness of fit and some comparisons. J. Am. Stat. Assoc. 69:730737.
Szabo, A. G., T. M. Stepanik, D. M. Wayner, and N. M. Young. 1983. Conformational heterogeneity of the cooper binding site in azurin. Biophys. J. 41:233244.
Wilk, G., and Z. W
odarczyk. 2001. Nonexponential decays and nonextensivity. Phys. Lett. A. 290:5558.
Wilk, G., and Z. W
odarczyk. 1999. Interpretation of the nonextensivity parameter q in some applications of Tsallis statistics and Lévy distributions. Phys. Rev. Lett. 84:27702773.
Zimmermann, Th., H. Köppel, L. S. Cederbaum, G. Persch, and W. Demtröder. 1988. Confirmation of random-matrix fluctuations in molecular spectra. Phys. Rev. Lett. 61:36.[Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |