| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
* Raymond and Beverly Sackler Faculty of Exact Sciences, School of Chemistry, Tel-Aviv University, Tel-Aviv, Israel; and Department of Physical Chemistry and the Fritz Haber Research Center, The Hebrew University, Jerusalem, Israel
Correspondence: Address reprint requests to Dr. Noam Agmon, Tel.: 972-2-658-5687; Fax: 972-2-651-3742; E-mail: agmon{at}fh.huji.ac.il.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
–4). Formation of the wild-type chromophore (p-hydroxybenzylidene-imidazolidinone) occurs via post-translational cyclization followed by auto-oxidation within a tripeptide unit (Ser-65, Tyr-66, Gly-67) of the 238 amino-acid sequence. X-ray diffraction (5,6) revealed that the chromophore (Cro) is protectively housed along a coaxial helix threaded through the center of an 11-stranded ß-barrel. It is covalently anchored and effectively secluded from the aqueous solvent surrounding the protein. Additional noncovalent coupling of the Cro to the protein backbone is facilitated via an extended hydrogen-bond (HB) network.
At room temperature (RT), wild-type (wt) GFP exhibits two main absorption peaks with maxima at 398 nm (band A) and 478 nm (band B). The ratio of the absorption bands depends weakly on the pH. As the pH increases, the B-band gets stronger while maintaining an isosbestic point indicative of two-state kinetics. At RT there are two emission peaks in the steady-state fluorescence spectrum, a very weak band at 460 nm and a very strong band at 510 nm. Time-resolved emission showed a decrease in the 460-nm emission concomitant with a rise in the 510-nm emission. This has been interpreted (7,8) as evidence for ESPT from the hydroxyl of Tyr-66 (R*OH form, A-state, 460 nm), resulting in the green fluorescence of the anion (R*O– form, 510 nm).
Previous time-resolved studies of excited GFP focused on the complex decay of the emission at short times up to 150 ps (7–10). Leiderman et al. (11) have extended the fluorescence measurements to long-times (up to 10 ns), because this time regime can reveal important facets of the recombination reaction between the conjugate base and the diffusing proton. Our experience with ESPT in homogeneous solutions (water, alcohols, dimethylsulfoxide) shows that this reaction is reversible in the excited state (ES), and thus it is possible to monitor the probability that the dissociated and solvated proton, which diffuses with a diffusion constant D, will recombine geminately with the deprotonated base (12,13). The theoretical formulation for this process in d-dimensions leads to the asymptotic (
) long-time behavior of the normalized fluorescence signal, If (t;
), at the wavelength characteristic of the protonated Cro:
 | (1) |
f is the fluorescence lifetime of the deprotonated form.
In solution, one consistently finds t–3/2 asymptotics (13), which is commensurate with the three-dimensional character of the homogeneous solution (d = 3 in Eq. 1). A surprising finding of Leiderman et al. (11) was that a similar behavior is observed for wt-GFP at RT. This was unexpected from the prevailing model (14,15) of irreversible proton migration along a short HB network connecting Tyr-66, via a water molecule and the hydroxyl of Ser-205, to the carboxylate of Glu-222. In Leiderman et al. (11), the t–3/2 decay was interpreted as evidence for three-dimensional proton diffusion. Since there was no effect of proton scavengers on the kinetics, it was concluded that the diffusion occurs inside the protein. These proton scavengers can be relatively high concentrations (0.4 M) of fluoride or acetate anions, which are added to the buffer solution.
A more detailed examination of the x-ray data (5,6) was consequently initiated (16). It was not possible to find three-dimensional HB arrays, which could support three-dimensional proton diffusion inside the GFP barrel. However, linear chains of HBs (i.e., proton wires) were found to be more prevalent within the GFP than previously reported. These proton wires involve a few crystallographic water molecules plus mainly oxygen atoms on amino-acid side chains as well as backbone carbonyls.
Two functionally relevant proton wires were observed:
) from RT to a wide temperature range, 87–291 K. Interestingly, we find a transition from the t–3/2 decay at RT to a t–1/2 behavior below 230 K. The implications are that at low temperatures the proton migrates within the protein, possibly along the one-dimensional HB chain of Fig. 1. The transition at higher temperatures is ascribed to a conformational change within the protein, which allows proton exit possibly via the threonine gateway of Fig. 2.
|
|
|
|
MATERIALS AND EXPERIMENTAL METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
Samples were placed in a 1-mm optical path-length quartz cell. The temperature of the irradiated sample was controlled by placing it in a liquid nitrogen cryostat with thermal stability of ±1 K. The absorbance at 397 nm was typically 0.1 OD. Steady-state fluorescence spectra at low temperature were taken using CVI SM-240 spectrometer (CVI Spectral Products, Putnam, CT) with 2 nm resolution.
Time-resolved measurements
Time-resolved fluorescence was acquired using the time-correlated single-photon counting (TCSPC) technique, the method of choice when sensitivity, large dynamic range, and low intensity illumination are important criteria in fluorescence decay measurements. For excitation we used a cavity-dumped Ti:sapphire femtosecond laser (Mira, Coherent, Santa Clara, CA), which provides short pulses (80 fs), operating at the second harmonic generation frequency (spectral range of 380–400 nm) and a relatively low repetition rate (500 kHz). A low rate might be important to ensure that the GFP is fully relaxed before re-excitation. The excitation pulse energy was reduced by neutral density filters to 1 pJ. We checked the sample before and after time-resolved measurements, not finding noticeable irradiation-induced changes in the absorption spectra or the transient emission decay. We conclude that, under our irradiation conditions, no sample deterioration could be detected.
The TCSPC detection system is based on a Hamamatsu 3809U photomultiplier (Hamamatsu City, Japan) and an Edinburgh Instruments TC900 computer module (West Lothian, United Kingdom) for TCSPC. The overall instrumental response function (IRF) was 40 ps (full-width half-maximum). The large dynamic range of the TCSPC system (more than four decades) enables us to accurately determine the non-exponential photoluminescence decay profiles of wt-GFP fluorescence. Measurements were taken at 10-nm spectral width at two frequencies: = 450–455 nm (R*OH kinetics) and = 505–525 nm (R*O– kinetics), as well as at = 485 nm (stabilized RO–, see below). The latter was used to determine the ES lifetime, f = 1/kf. The lifetime-corrected signals, If(t;) exp(t/f), when measured at the frequencies of the acid and the anion and normalized, are denoted below by P(t) and 1 – P(t), respectively.
|
THEORY |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
 | (2) |
f, which is assumed to be identical for both R*OH and R*O–. The effect of ES decay in this case is simply to multiply the calculated kinetics without decay by exp(–t/f). Thus, we need to solve the kinetics of Eqs. 2 without the ES decay.
Since the reaction under consideration is inherently fast, it is not possible to neglect the diffusive transport of the proton, which couples to the reactions under consideration. In the diffusion model (12,13), one assumes that the proton can execute spherically symmetric diffusive motion in a d-dimensional space. At the contact distance, r = a, it may associate with the excited base with a rate coefficient ka or dissociate from it with a rate coefficient kd. The competing irreversible quenching reaction (rate coefficient kq) is assumed, for simplicity, to occur at the same contact distance a. The potential V(r) (in units of the thermal energy, kBT) describes the proton-anion interaction as a function of their separation, r. In solution-phase studies (13), a Coulomb (or screened Coulomb) potential was used but, as we shall see below, this is inapplicable for the protein interior. The association (equilibrium) constant is given by K*eq = ka e–V(a)/kd, which includes the effect of the potential at the contact distance.
The mathematical depiction of this problem utilizes the d-dimensional spherically symmetric Smoluchowski equation
 | (3) |
 | (4) |
 | (5) |
) for p(r, t) at r = a couples the reversible and irreversible reactions:
 | (6) |
d = 1, 2
, and 4 for d = 1, 2, and 3, respectively. Some limiting cases of this equation can be solved analytically (see below). However, for arbitrary times and potentials one needs to solve Eq. 3 numerically, and this is conveniently done using the Windows application (18) for solving the Spherically Symmetric Diffusion Problem (SSDP, Ver. 2.65). Its graphical interface allows one to convolute the theoretical solution with the IRF and compare with the time-resolved experimental data.
Analytic solution
A central result of the present analysis is that, below a critical temperature (of 230 K), the fluorescence signal from the excited GFP Cro behaves as if the proton diffuses in one dimension—namely, d = 1 in Eqs. 3 and 6. In this case, the above equations admit an analytic solution in several limiting cases. The simplest case is that of no potential (unbiased diffusion) and no quenching, V(r) = 0 and kq = 0. This problem has been solved two decades ago (17). Interestingly, we will see that this simple case almost suffices to describe our low temperature results.
We present this solution below, mainly to show (19) that from the time dependence of P(t) it is possible to determine only two combinations of the three parameters, kd, ka, and D. These combinations are kd and 
both having units of 1/time. It is not possible to determine ka and D independently, because these two parameters depend also on the units of r, whereas our data depends only on time. Thus, although diffusion is indispensable in explaining the data, the diffusion constant itself can be determined only from measurements that sample the spatial distribution, p(r, t).
The analytic solution can be conveniently obtained (17) by the Laplace transform method. It depends on two roots of a quadratic equation
 | (7) |
 | (8) |
(z) is defined by
 | (9) |
±, and from their product and sum,
 | (10) |
The long-time asymptotic behavior of P(t) can be obtained directly from that of erfc(z); see Abramowitz and Stegun (20). This can be worked out even when V(r) 
0, provided that V(r) 
0 as r 
(21, 22). It leads to Eq. 1 with K*eq = ka e–V(a)/kd. Thus, when only the long-time asymptotics is measured, its slope allows us to determine only a single parameter, which is the ratio of the two parameters in Eq. 10.
More general cases with kq 0 and/or a linear potential, V(r) = 
r, can also be worked out analytically (23, 24, 25), though the solution becomes more complicated. In these cases the additional parameters that could possibly be determined from fits to the complete time-dependence of P(t) are 
and 
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
510 nm at RT and blue-shifting by 5 nm at the lowest T. This is attributed to R*O– formed by the photodissociation of R*OH. Separated from it by an iso-emissive point, a smaller band at 485 nm grows with decreasing temperature. This band has previously been observed at 77 K (7). Now we see that it appears rather abruptly below 230 K.
|
|
)), one concludes that it is due to direct excitation of ground-state RO–. Possibly the protein conformation of this species, or the absence of an available proton within the protein, prevent its conversion to ground-state ROH, which characterizes the majority of the protein molecules in the sample. ). In contrast, the protein molecules emitting at 485 nm do not contain an available proton within their barrel structure that could interact with the Cro.
Interestingly, the 485-nm peak appears rather abruptly as T is lowered below 230 K. This transition could indicate the onset of a protein conformational change. The blue shift, and its enhancement at low T, suggests ground-state stabilization of the anionic species. Because the transient signal at 485 nm is a pure exponential, the population emitting at 510 nm does not convert to the 485-nm population during the ES lifetime, so that the conformational change must occur on a longer timescale in the ground-state. Hence, this population is decoupled from the low-T ES acid-base reaction, and could thus be ignored in the treatment of the ES kinetics. In contrast, the main 510-nm band continues shifting and widening to the red also above 230 K. Thus its broadening could involve less localized protein modes with a wider distribution of barrier heights.
Acid fluorescence
The ROH emission occurs around 460 nm. It is barely seen in the steady-state spectrum (Fig. 3) because it is 500 times weaker than the anion emission. Nevertheless, its intensity is sufficient for TCSPC measurements. The central result of the present work is the temperature-dependence of the time-resolved fluorescence for wt-GFP (Fig. 5). Here we present results in 30% glycerol-water (v/v) solutions, although we have obtained qualitatively similar data for other glycerol contents (12% and 60%). The data are multiplied by exp(kft), where kf = 0.35 ns–1 was determined from the pure exponential decay of the emission from the directly excited anion at 485 nm. For the acid form it is not possible to measure kf directly due to its rapid decay, so we apply the same value, assuming that (like the 485-nm band) it has no temperature-dependence.
|
230 K. (Note that this temperature is close to the transition temperature of the 485-nm band in Fig. 3). The data below this temperature are presented in Fig. 5, A and B. Fig. 5 A gives the short time behavior on a linear scale, whereas Fig. 5 B shows the long-time behavior on a log-log scale (with data shifted vertically for clarity). It exhibits a power law asymptotic decay as in solution-phase ESPT (13), except that the power is t–1/2 rather than t–3/2. This strongly suggests that the proton is migrating along one-dimensional proton wires within the protein.
To check this hypothesis, we have fitted the data to the one-dimensional diffusion model, solving Eq. 3 numerically for d = 1 to obtain P(t) (lines). The same calculation is used to compare 1 – P(t) with the RO– kinetics, as measured in the 505–510 nm range (Fig. 5 C). The parameters, collected in Table 1, correspond in a first approximation to free diffusion (V(r) 
0) with the simplest reversible boundary condition (i.e., kq 0). As shown in the Theory section and in Agmon (19), the transient traces in this limit allow one to determine kd and It is not possible, from an experiment containing only temporal information, to determine independently the two spatially dependent parameters, ka and D. In liquid solutions D is known from mobility measurements, but within the protein D is unknown.
To achieve a global fit, which agrees simultaneously with the kinetics at both acid and base wavelengths, we have introduced two additional effects. First, we included a small proton quenching effect with a rate coefficient, kq, which is typically 4% of ka. It accounts for the decrease in the R*O– signal at long-times (Fig. 5 C, inset). Such a decrease can also be explained by increasing kf, but in our case kf has been determined independently from the decay of the 485-nm peak, and was found to be temperature-independent (below 230 K). Note that our depiction of the quenching reaction is least accurate at short times, because Eq. 6 assumes that it occurs at the same contact distance as dissociation and recombination, whereas it probably takes place at a site on the back ring of the Cro (on atom N2 in Fig. 1). However, at longer times the proton migrates sufficiently far that the difference in distance to the two sites becomes relatively small.
The second effect is a weak linear attractive potential (in one dimension), V(r) = r. It enhances the tail of the R*OH signal, which is otherwise diminished by the quenching effect. We find that the attracting Coulomb potential traditionally used for ESPT in solution (13) is inadequate, as it leads to excessive recombination at short and intermediate times, even if a much larger dielectric constant is used for the interior of the protein than commonly assumed. The success of electrostatic considerations on nanometric scales is contrasted with their inadequacy on the molecular scale.
By fitting the time-resolved kinetics, one can determine at most five parameters: kf, kd, 
and The radiative decay rate constant, kf, was kept constant throughout the whole temperature range at the value obtained for the 485-nm band below 230 K. This leaves us with four adjustable parameters, 
kq/ka, and ka (Table 1). The last two parameters were found to be temperature-independent below 200 K. Above it, they decrease to zero over a rather small temperature range. Thus, below 200 K we adjust only two parameters at each temperature. The good fits obtained simultaneously to both ROH and RO– kinetics on all timescales lend credibility to the one-dimensional diffusion model below 230 K.
Above this temperature, the one-dimensional model deteriorates with increasing T (Fig. 5 D). There is a switchover to steeper decay occurring at earlier times with increasing T. Above 270 K, the long-time tail seems to approach the t–3/2 asymptotic behavior (dashed line) as previously reported (11). This could indicate that, on these timescales (say, t > 300 ps at RT), the proton is exiting from the protein (16) either irreversibly or reversibly, after executing three-dimensional diffusion in the external solution. Commensurate with the proton exit event, the proton-quenching process is eliminated (see red curve in Fig. 5 C, inset), and kq decreases to zero.
Another signal of the transition temperature is seen in the temperature-dependence of the dissociation rate constant in Fig. 6. At very low T, kd is nearly constant, suggesting the involvement of a quantum tunneling effect. Then kd increases with increasing T up to 230 K, where it reaches a maximum. Possibly the conformational change occurring near the transition temperature modifies the HB network around the Cro in a way that slows the dissociation process at higher T. At this stage, however, we do not analyze this temperature-dependence quantitatively, because an IRF of 40 ps can distort the fast dissociation process.
|
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
). The weaker R*OH band was studied at short times (100 ps) to demonstrate the ESPT phenomenon (7–10), but its long-time behavior remained unexplored until recently (11). In addition, there were no systematic studies of the temperature-dependence. Here we presented a first report of the temperature-dependence of the transient fluorescence from wt-GFP, with particular emphasis on the R*OH band up to t > 10 ns. We find that the transient fluorescence kinetics exhibit a transition around 230 K. Below it, there is t–1/2 asymptotic decay, with the whole time-dependence fitting quantitatively to a one-dimensional diffusion model for a geminate reversible reaction, indicating that the proton is migrating along one-dimensional proton wires within the protein. Above 230 K, the long-time behavior switches to a steeper t–3/2 decay, suggesting a conformational change that allows the proton to escape to a region from which its probability of return to the active site is much more limited.
Proton wires
It is commonly assumed that proton conduction within proteins occurs along proton wires, which are specialized HB networks for transporting protons within proteins. Such wires are used to conduct protons to, from, or between active sites in enzymes or across the cell membrane in proton pumps (26,27). Evidence for proton wires is typically obtained from x-ray structures and mutagenesis work. However, we are not aware that proton motion along such networks has ever been recorded in real-time.
Below 230 K, wt-GFP seems to be frozen into two major conformations in its ground electronic state. The minority is locked in the anionic state, giving rise to the 485-nm fluorescence band. When excited, it decays mono-exponentially with a rate constant kf = 0.35 ns–1, which is independent of temperature in this range, allowing us to obtain a reliable value for the ES lifetime, f = 1/kf.
The majority of the population is in the ROH state, in which the OH of Tyr-66 is connected to a rather rigid HB network. Its fluorescence shows a remarkably slow decay with the ubiquitous t–1/2 tail characteristic of reversible proton migration in one dimension. The anion shows a small excess decay (over that of the directly-excited minority form at 485 nm), suggesting a small degree of quenching by the dissociated proton. This interpretation is corroborated by fits to a one-dimensional diffusion model for reversible dissociation-association with quenching (Fig. 5 B). Rigorously speaking, the t–1/2 asymptotics should hold only in the absence of quenching and for potentials that decay to zero at infinity. For a linear potential (as employed here), the asymptotic behavior shows transitions between various forms, such as t–3/2, t–1/2 and exponential (23–25). In the present case, however, the competing quenching reaction is minimal (kq/ka 0.04). It affects the long-time behavior only marginally in our observation time-window, causing minor deviations from the t–1/2 behavior.
Thus, in accordance with the prevailing mechanism for proton migration in GFP (14,15), the proton is transferred in the ES from Tyr-66 to Glu-222. However, this transfer is reversible and probably continues beyond Glu-222. The pathway up to Glu-222 involves only two intermediate proton acceptors: water and Ser-205. Irreversible transfer along this route should result in biexponential kinetics. Conceivably, longer HB chains are required to allow the diffusional power-law to evolve, such as the one continuing from Glu-222 to Asp-82 (Fig. 1 and (16)).
Our measurements, devoid of a characteristic length-scale, cannot determine how far within the protein the proton migrates. This fundamental question might, however, be addressed in combination with mutagenetic work. If the proton-wire could be bisected at points progressively closer to the proton release point, the power law should terminate by a flat plateau at progressively shorter times. Alternatively, blocking of proton entry might lock the Cro in its RO– state indefinitely, precluding the observation of ESPT altogether.
The possibility that proton migration is not limited to a single linear pathway follows from the observation of proton quenching. Quantum chemistry calculations (28) suggest that the GFP Cro is actually a bifunctional compound (29), with the N2 nitrogen on the imidazolidinone enhancing its basicity in the ES. Attack of the geminate proton on this site will form the zwitterion (30), which decays rapidly via conical intersection to the ground-state (31). Indeed, Fig. 1 shows an HB connecting a water molecule along the Asp-82–Glu-222 pathway to N2 (yellow), and this could allow the formation of the dark zwitterion as manifested in the quenching kinetics.
The possibility of seeing proton-wires using ESPT could potentially be utilized to map out conjectured pathways in other proteins, by artificially tailoring an ESPT Cro at a location where connection to the HB pathway could be made. The formation of the desired connection could be verified by x ray, and proton migration along the pathway could then be studied by time-resolved fluorescence.
Proton exit
As the temperature is increased above 230 K, the probability of the geminate proton to return to the dissociation site on Tyr-66 and rebind diminishes, and then the fluorescence signal from the protonated Cro drops below the t–1/2 behavior observed for T < 230 K (Fig. 5 D). This could result from a conformational transition occurring within the protein, one which is also connected to the disappearance of the 485-nm band in Fig. 3. This conformational change does not occur instantly upon excitation, but rather with some delay, which is manifested in the residual t–1/2 kinetics observed at intermediate times before the conformational change. The delay time increases with decreasing T, until it eventually becomes infinite below 230 K. Clearly, in a complex protein structure many different conformational changes may occur. Here we focus on one such change that has been noted in comparing the x-ray structures of the protonated and deprotonated forms of wt-GFP (14).
It is known that the side chain of Thr-203 can assume two conformations in wt-GFP. In the ROH form, Thr-203 points its OH group away from that of Tyr-66. In the RO– form, the Thr-203 side chain is rotated by 120°, so that its OH forms a HB to the phenoxide, RO–···HO. This is demonstrated in Fig. 2. If this rotational motion freezes below 230 K, there would be a fraction of the population which is trapped indefinitely in the RO– state, giving rise to the 485-nm band. Moreover, it has been noticed (16) that the rotated Thr-203 forms a short HB-network leading from Tyr-66 to the backbone carbonyl of His-148, as seen in Fig. 2. This group is already on the surface of the protein, which should allow for facile proton escape to solution. In the ground-state, the majority of the population is in the ROH state. Below 230 K, the threonine switch does not open after photoexcitation, so that the ROH releases its proton to the long proton wire of Fig. 1. Above 230 K, the switch opens (with a delay of a few 100 ps), and then the proton hops across the short HB connection seen in Fig. 2 to the external solution.
As noted before (11), the steeper decay at the higher temperatures shows a t–3/2 dependence, which could be indicative of proton diffusion in three-dimensional space with reversible recombination (13). Diffusion is now taking place outside the protein, with occasional reentry and rebinding to the excited Cro (16). However, this reversible exit scenario is at odds with the experimental finding of essentially zero effect for relatively high concentrations of external proton scavengers (11).
An alternative interpretation may be that proton exit is irreversible, because it is difficult for it to find its way back into the small hole on the protein surface. Indeed, pH jump experiments (32,33) suggest that proton entry from bulk solution occurs on the 100-µs timescale. Thus, in this scenario, after the proton escapes it can no longer affect the Cro fluorescence. The subsequent non-exponential decay might then arise from a different effect, such as an inhomogeneous distribution in the Thr-203 side-chain angle which determines the exit rate.
An independent evidence for proton exit comes from the observation of cessation of geminate quenching above 230 K (Fig. 5 C, inset). The dissociated proton outside the protein cannot reach the N2 site at the back side of the Cro, where the quenching occurs. The lifetime of the R*O– formed by protic dissociation (510-nm band) then increases, coinciding with that of the directly excited RO– population (485-nm band).
These conclusions do not contradict the transient IR measurements of Stoner-Ma et al. (34), which indicate that Glu-222 gets (partially) protonated in deuterated GFP at RT on a 300-ps timescale. Until the Thr-203 switch opens, the proton wanders along the internal proton wire(s), in accordance with the fact that the one-dimensional diffusion model still fits the data at short times in Fig. 5 D. At 281 K we observe the transition to the t–3/2 decay at 300 ps, while for deuterated GFP this should occur even at later times. Thus we expect the transient IR signal to start decreasing at longer times.
The 230 K transition
An interesting observation is the transition in GFP steady-state and transient fluorescence around 230 K. It is manifested at three different characteristic wavelengths:
485 nm
The emission at 485 nm is attributed to anionic GFP, which does not get reprotonated even after long residence times in the ground-state. This small peak in the emission spectrum (Fig. 3) disappears abruptly near 230 K.
510 nm
Emission at this wavelength is attributed to the GFP anion, which is formed in the ES from the acidic R*OH form. Unlike the 485 K population, its fluorescence signal at short-times grows with a delay. At long-times and above 230 K, the decay at both wavelengths is exponential with a time constant of 2.9 ns. Below 230 K, the emission at this wavelength decays faster, which we attribute to quenching by the geminate proton. Thus the quenching reaction in Eq. 2 bears evidence to the 230 K transition.
450 nm
The weak emission at 460 nm is attributed to excited GFP in its R*OH form. The transient kinetics show two signatures of the 230 K transition:
The accumulation of several different experimental evidences for the transition temperature indicates that there might be an observable change in the protein structure around T = 230 K. It may therefore be interesting to re-measure the x-ray diffraction of a wt-GFP sample above and below this temperature to check this possibility.
Proton re-entry
Our fluorescence measurements, which follow the ES on a 20-ns timescale, do not probe the proton reentry process. Here, two types of measurements are relevant. By monitoring the transient absorption of the ROH band, one can follow its ground-state recovery after the ES deprotonation process. It was found to proceed nearly in parallel to ES decay, with a time-constant of 3–4 ns (35). In contrast, pH-jump experiments (33) indicate that proton entry into anionic GFP mutants occurs on the much slower 100-µs timescale. It could be that the two kinds of experiments track different GFP populations. If the R*O– decays before proton exit has occurred, reprotonation by the geminate proton can be very fast. However, for those GFP molecules that have lost their proton, reprotonation requires proton entry from solution, which is a much slower process.
Two proton wires connecting the Tyr-66 OH with the protein surface are depicted in Figs. 1 and 2. As discussed in Agmon (16), the protein surfaces near their end-points have completely different characteristics. The short pathway in Fig. 2 ends on a hydrophobic surface patch, which should not be attractive to protons. The long pathway in Fig. 1 ends at Glu-5 on the bottom surface of the GFP barrel. This face contains many Glu and Asp residues, which could serve to collect protons from solution. The structure of this pathway is analogous to that of the D-pathway in cytochrome-c oxidase, which serves to conduct protons into the active site of that protein. By homology, proton reentry into GFP from solution at RT could proceed mainly along the long pathway in Fig. 1.
|
CONCLUSION |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
), or become distributed over so many proton-wires that its behavior appears three-dimensional (11). The unique combination of proton wires with an ESPT Cro within the same protein allows one to follow these proton migration processes in real-time using transient fluorescence, which was not achieved in studies of other proteins thus far.
Our findings are consistent with the mechanism suggested in Agmon (16), where GFP under normal conditions acts as a proton pump. In this scenario, photo-excitation induces rotation of the Thr-203 side chain leading to proton ejection to solution (Fig. 2). The proton is then replenished in the ground-state from another face of the protein via the long internal proton wire of Fig. 1. From this viewpoint, GFP has similar features to those of bacteriorhodopsin (br), for example, which makes it an efficient proton pump. Both proteins have a Cro that utilizes light energy to induce proton transfer: ESPT in GFP and ground-state PT in br. Both proteins couple PT and conformational change to ensure the directionality of the transfer process. The isomerization occurs within the Cro in br or right next to it (Thr-203) in GFP. Both have postulated exit and entry paths for the proton along intraprotein proton wires. These wires are not clearly visible in br, which has many labile water molecules. Remarkably, in the rigid GFP barrel these proton conduction pathways are much more stable, so that they could be inferred (16) directly from the x-ray structure (Figs. 1 and 2). The combination of more stable proton wires, with the faster PT occurring in the ES, could actually make the proton-pumping cycle in GFP faster than that of br.
If these conjectures can be corroborated, GFP mutants might have a new line of applications, because proton pumping proteins are candidates for devices capable of direct light conversion into proton gradients (36). GFP is considerably more rigid than known membrane proton pumping proteins, so that it may be expected to exhibit higher stability in such devices.
In marine biology, bioluminescence (which derives its excitation energy from a chemical reaction) is actually a much more prevalent phenomenon than biofluorescence. It has recently been suggested that ESPT plays a role also in bioluminescence, e.g., from aequorin and coelenterazine-based photoproteins in general (37). Thus, the detailed study of ESPT in photoproteins may be essential also for understanding bioluminescence phenomena on the molecular level.
|
ACKNOWLEDGEMENTS |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
This research was supported in part by The Israel Science Foundation (N.A., grant No. 191/03), the U.S.-Israel Binational Science Foundation (D.H.), and the James-Franck German-Israel Program in Laser-Matter Interaction (D.H.).
Submitted on June 24, 2005; accepted for publication October 24, 2005.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMATERIALS AND EXPERIMENTAL...THEORYRESULTSDISCUSSIONCONCLUSIONACKNOWLEDGEMENTSREFERENCES |
|---|
2. Tsien, R. Y. 1998. The green fluorescent protein. Annu. Rev. Biochem. 67:509–544.[CrossRef][Medline]
3. Remington, S. J. 2000. Structural basis for understanding spectral variations in green fluorescent protein. Meth. Enzymol. 305:196–210.[CrossRef][Medline]
4. Zimmer, M. 2002. Green fluorescent protein (GFP): applications, structure, and related photophysical behavior. Chem. Rev. 102:759–781.[CrossRef][Medline]
5. Yang, F., L. G. Moss, and G. N. Phillips, Jr. 1996. The molecular structure of green fluorescent protein. Nat. Biotechnol. 14:1246–1251.[CrossRef][Medline]
6. Örmo, M., A. B. Cubitt, K. Kallio, L. A. Gross, R. Y. Tsien, and S. J. Remington. 1996. Crystal structure of the Aequorea victoria green fluorescent protein. Science. 273:1392–1395.[Abstract]
7. Chattoraj, M., B. A. King, G. U. Bublitz, and S. G. Boxer. 1996. Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer. Proc. Natl. Acad. Sci. USA. 93:8362–8367.
8. Lossau, H., A. Kummer, R. Heinecke, F. Pöllinger-Dammer, C. Kompa, G. Bieser, T. Jonsson, C. M. Silva, M. M. Yang, D. C. Youvan, and M. E. Michel-Beyerle. 1996. Time-resolved spectroscopy of wild-type and mutant green fluorescent proteins reveals excited state deprotonation consistent with fluorophore-protein interactions. Chem. Phys. 213:l–16.
9. Striker, G., V. Subramaniam, C. A. M. Seidel, and A. Volkmer. 1999. Photochromicity and fluorescence lifetimes of green fluorescent protein. J. Phys. Chem. B. 103:8612–8617.
10. Winkler, K., J. Lindner, H. Bürsing, and P. Vöhringer. 2002. Ultrafast Raman-induced Kerr effect of water: single molecule versus collective motions. J. Chem. Phys. 113:4674–4682.[CrossRef]
11. Leiderman, P., M. Ben-Ziv, L. Genosar, D. Huppert, K. M. Solntsev, and L. M. Tolbert. 2004. Study of the long-time fluorescence tail of the green fluorescent protein. J. Phys. Chem. B. 108:8043–8053.
12. Pines, E., D. Huppert, and N. Agmon. 1988. Geminate recombination in excited-state proton transfer reactions: numerical solution of the Debye-Smoluchowski equation with back-reaction boundary conditions. J. Chem. Phys. 88:5620–5630.[CrossRef]
13. Agmon, N. 2005. Elementary steps in excited-state proton transfer. J. Phys. Chem. A. 109:13–35.[CrossRef][Medline]
14. Brejc, K., T. K. Sixma, P. A. Kitts, S. R. Kain, R. Y. Tsien, M. Ormö, and S. J. Remington. 1997. Structural basis for dual excitation and photoisomerization of the Aequorea victoria green fluorescent protein. Proc. Natl. Acad. Sci. USA. 94:2306–2311.
15. Palm, G. J., A. Zdanov, G. A. Gaitanaris, R. Stauber, G. N. Pavlakis, and A. Wlodawer. 1997. The structural basis for spectral variations in green fluorescent protein. Nat. Struct. Biol. 4:361–365.[CrossRef][Medline]
16. Agmon, N. 2005. Proton pathways in green fluorescence protein. Biophys. J. 88:2452–2461.
17. Agmon, N. 1984. Diffusion with back-reaction. J. Chem. Phys. 81:2811–2817.[CrossRef]
18. Krissinel', E. B. and N. Agmon. 1996. Spherical symmetric diffusion problem. J. Comput. Chem. 17:1085–1098.[CrossRef]
19. Agmon, N. 2005. Reduced parameter set determinable from geminate kinetics. Chem. Phys. Lett. 417:530–534.[CrossRef]
20. Abramowitz, M., and I. A. Stegun, editors. 1970. Handbook of Mathematical Functions. Dover, New York.
21. Agmon, N., E. Pines, and D. Huppert. 1988. Geminate recombination in proton-transfer reactions. II. Comparison of diffusional and kinetic schemes. J. Chem. Phys. 88:5631–5638.[CrossRef]
22. Agmon, N., and A. Szabo. 1990. Theory of reversible diffusion-influenced reactions. J. Chem. Phys. 92:5270–5284.[CrossRef]
23. Kim, H., K. J. Shin, and N. Agmon. 1999. Excited-state reversible geminate recombination with quenching in one dimension. J. Chem. Phys. 111:3791–3799.[CrossRef]
24. Kim, H., K. J. Shin, and N. Agmon. 2001. Diffusion-influenced reversible geminate recombination in one dimension. II. Effect of a constant field. J. Chem. Phys. 114:3905–3912.[CrossRef]
25. Kim, H., and K. J. Shin. 2004. Diffusion-influenced reversible geminate recombination in one dimension. III. Field effect on the excited-state reaction. J. Chem. Phys. 120:9142–9150.[CrossRef][Medline]
26. Decoursey, T. E. 2003. Voltage-gated proton channels and other proton transfer pathways. Physiol. Rev. 83:475–579.
27. Frank, R. A. W., C. M. Titman, J. V. Pratap, B. F. Luisi, and R. N. Perham. 2004. A molecular switch and proton wire synchronize the active sites in thiamine enzymes. Science. 306:872–876.
28. Scharnagl, C., and R. A. Raupp-Kossmann. 2004. Solution pKa values of the green fluorescent protein chromophore from hybrid quantum-classical calculations. J. Phys. Chem. B. 108:477–489.
29. Bardez, E. 1999. Excited-state proton transfer in bifunctional compounds. Isr. J. Chem. 39:319–332.
30. Voityuk, A. A., M.-E. Michel-Beyerle, and N. Rösch. 1998. Quantum chemical modeling of structure and absorption spectra of the chromophore in green fluorescent proteins. Chem. Phys. 231:13–25.[CrossRef]
31. Weber, W., V. Helms, J. A. McCammon, and P. W. Langhoff. 1999. Shedding light on the dark and weakly fluorescent states of green fluorescent proteins. Proc. Natl. Acad. Sci. USA. 96:6177–6182.
32. Mallik, R., J. B. Udgaonkar, and G. Krishnamoorthy. 2003. Kinetics of proton transfer in a green fluorescent protein: a laser-induced pH jump study. Proc. Indian Acad. Sci. (Chem. Sci.) 115:307–317.
33. Saxena, A. M., J. B. Udgaonkar, and G. Krishnamoorthy. 2005. Protein dynamics control proton transfer from bulk solvent to protein interior: a case study with a green fluorescent protein. Protein Sci. 14:1787–1799.
34. Stoner-Ma, D., A. A. Jaye, P. Matousek, M. Towrie, S. R. Meech, and P. J. Tonge. 2005. Observation of excited-state proton transfer in green fluorescent protein using ultrafast vibrational spectroscopy. J. Am. Chem. Soc. 127:2864–2865.[CrossRef][Medline]
35. Kummer, A. D., J. Wiehler, H. Rehaber, C. Kompa, B. Steipe, and M. E. Michel-Beyerle. 2000. Effects of threonine 203 replacements on excited-state dynamics and fluorescence properties of the green fluorescent protein (GFP). J. Phys. Chem. B. 104:4791–4798.
36. Luo, T.-J. M., R. Soong, E. Lan, B. Dunn, and C. Montemagno. 2005. Photo-induced proton gradients and ATP biosynthesis produced by vesicles encapsulated in a silica matrix. Nature Mater. 4:220–224.[CrossRef][Medline]
37. Vysotski, E. S., and J. Lee. 2004. Ca2+-regulated photoproteins: structural insight into the bioluminescence mechanism. Acc. Chem. Res. 37:405–415.[CrossRef][Medline]
This article has been cited by other articles:
 |
|
 |
|
  A. B. Copty, Y. Neve-Oz, I. Barak, M. Golosovsky, and D. Davidov Evidence for a Specific Microwave Radiation Effect on the Green Fluorescent Protein Biophys. J., August 15, 2006; 91(4): 1413 - 1423. [Abstract] [Full Text] [PDF]
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||
