Ultrafast excited state dynamics of the green fl uorescent protein chromophore and its kindling fl uorescent protein analogue

Fluorescent proteins exhibit a very diverse range of photochemical behaviour, from efficient fluorescence through photochromism to photochemical reactivity. Remarkably this diverse behaviour arises from chromophores which have very similar structures. Here we describe measurements and modelling of the excited state dynamics in the chromophores of GFP (HBDI) and the kindling fluorescent protein, KFP (FHBMI). The methods are ultrafast fluorescence spectroscopy with sub 50 fs time resolution and the modelling is based on the Smoluchowski equation. The excited state decays of both chromophores are very fast, longer for their anions than for the neutral form and independent of wavelength. Detailed studies show the mean fluorescence wavelength to be independent of time. The excited state decay times are also observed to be a very weak function of solvent polarity and viscosity. These results are modelled utilising recently calculated potential energy surfaces for the ground and excited states as a function of the twist coordinates about the two bridging bonds of the chromophore. For FHBMI and the scarce data on the neutral HBDI the calculations are not successful suggesting the need for refinement of these potential energy surfaces. For HBDI in methanol the simulation is successful provided a strong dependence of the radiationless decay rate on the coordinate is assumed. Such dependence should be included in future calculations of excited state dynamics. When the simulations are extended to more viscous solvents they fail to reproduce the observed weak viscosity dependence. The implications of these results for the nature of the coordinate leading to radiationless decay in the chromophore and for the photodynamics of fluorescent proteins are discussed.


Introduction
2][3][4] GFP, and its more photostable School of Chemistry, University of East Anglia, Norwich Research Park, Norwich NR4 7TJ.E-mail: s.meech@ uea.ac.uk and uorescent mutants, 5 are readily cloned and expressed attached to a target protein, making the expression and subsequent location of that protein in the living cell visible to uorescence imaging.This new technology sparked a revolution in bioimaging which continues to the present time.The chromophore responsible for the green uorescence, 4-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI, Fig. 1) was characterized and isolated by Niwa and co-workers. 6,7The chromophore is formed in vivo in an intramolecular reaction between three adjacent amino acids. 8Mutations in GFP can modify either the chromophore or its immediate environment in a way that tunes the spectrum across the entire visible range, facilitating multicolour imaging and the imaging of protein-protein interactions through uorescence resonance energy transfer. 5ecent advances in protein engineering mean that the chromophore can be modied with unnatural amino acids, expanding the range of possible colours. 9It has been shown that some relatively minor changes to the chromophore (such as methylation in the phenyl group, the location of the OH group) can have a dramatic effect on the electronic spectrum. 10,113][14] Intriguingly, HBDI (or chromophores very similar to it) demonstrate an extraordinarily diverse range of photophysical behaviour in these different protein environments.While the anionic state of the chromophore in wild type GFP is strongly uorescent, the same state in other FPs is only weakly uorescent (such proteins being better labelled as chromoproteins).In other FPs the chromophore undergoes a photochemical reaction with an adjacent amino acid to generate a new chromophore with a different emission spectrum. 15,16Such proteins (e.g.kaede) have been used as 'optical highlighters' to label specic populations of proteins in situ.][22] This unexpectedly diverse behaviour has fuelled interest in the photophysics of HBDI itself.It was established that the strong uorescence of GFP is quenched when the protein is denatured, and that HBDI in uid solution is very weakly uorescent (quantum yield of ca 0.2 Â 10 À4 ).However, uorescence is recovered by freezing the solution to the glass transition temperature. 6,7,23,246][27] A number of derivatives of HBDI have been prepared that modify, in some cases quite dramatically, the electronic spectra of the chromophore.However, in most cases radiationless decay remains ultrafast. 11,28One feature these derivatives have in common is the methine bridge linking the two aromatic rings, so it appears likely that the radiationless decay coordinate involves reorganization about this group.Consistent with this, illumination of HBDI in its stable cis state was found to lead to formation of the metastable trans form. 29The spectra of the two isomers are similar and their excited state decays are indistinguishable. 30Further evidence for the importance of structural reorganization in the methine bonds comes from the observed suppression of radiationless decay when motion about these bonds is restricted, 31,32 and from the observation that a sterically crowded derivative of HBDI, with a twisted ground state, exhibits enhanced radiationless relaxation. 11etailed experimental studies of the excited state decay of HBDI showed that the decay rate is only moderately sensitive to solvent viscosity and solvent polarity, although uorescence is recovered as the glass transition temperature is approached, i.e. in very viscous media. 26,27,33The observed temperature dependence and weak dependence on viscosity suggested that a volume-conserving coordinate couples the excited state to the ground state over a small or zero energy barrier.The weak dependence on solvent polarity suggests that the radiationless decay does not require stabilization of a charge transfer state.
Early theoretical calculations suggested that the most likely volume-conserving coordinate associated with the methine bridge, the "hula twist", is energetically up-hill, and so incompatible with the observed near barrierless ultrafast relaxation. 34,35Thus, the details of the radiationless decay mechanism and the mechanism by which the protein matrix renders the chromophore strongly uorescent are incompletely understood.7][38][39][40][41][42][43][44][45][46] In a number of cases these calculations (which mainly refer to the gas phase) have investigated excited state dynamics, and two possible pathways leading to radiationless decay to the ground state have been suggested.These involve rotations about either of the two methine bridge bonds.In both cases the coordinates are barrierless or nearly so, consistent with the observed ultrafast decay in uid solution.These calculations are also in good agreement with recent experimental observation of the ultrafast dynamics of HBDI in the gas phase, where the initially excited Franck-Condon state relaxed to a twisted intermediate in a few hundred femtoseconds, and from the twisted state to the ground state in about 1 ps. 47The success of these calculations in the gas phase is important, but does not address the mechanism of uorescence enhancement in the protein.There have been a number of attempts to incorporate some representation of the medium into the calculations, either as a continuum or using the actual structure of the protein. 44,48However, none of these include dynamics associated with the environment.
In this work we extend our time-resolved studies of HBDI uorescence to include a detailed analysis of the time-dependent emission spectrum.We also compare data on HBDI with measurements on a newly synthesised derivative of HBDI related to the chromophore of the kindling uorescent protein (KFP), specically 2-formyl-4-hydroxybenzylidene-3-methylimidazolinone, FHBMI, Fig. 1.KFP is one of the photoactive proteins in which a uorescent state can be photochemically generated; 49,50 the excited state potential energy surface of the KFP chromophore (2-acetyl-4-hydroxybenzylidene-3-methylimidazolinone, FHBMI) was investigated by Olsen. 37The time-dependent spectra observed here are modelled using a generalized Smoluchowski equation.6][57][58][59][60][61] Here we replace the model PES with one determined from theoretical calculations, specically the calculations of Olsen for both HBDI and AHBMI in their anion states. 36,37We then incorporate realistic models of the solvent friction experienced by the assumed reaction coordinate for the FP chromophore.We further extend these calculations using the model PESs generated from the much more limited calculations on the neutral chromophore, 45 and thus discuss the neutral chromophore dynamics.In this way we endeavour to combine realistic calculated potential energy surfaces with a realistic model of the medium.

Experimental
The experimental arrangement for uorescence upconversion has been described in detail elsewhere, and is based on a 1 W 100 MHz titanium sapphire laser generating sub 20 fs pulses centered at 800 nm. 55The excitation wavelength was usually the second harmonic at 400 nm with power up to 9 mW.Fluorescence was detected by a reective microscope objective and focused into a thin nonlinear crystal where it was mixed with the residual 800 nm pulses, which were routed via a delay stage that controls the timing with femtosecond accuracy.The resulting sum frequency intensity samples the instantaneous uorescence intensity at the arrival time of the 800 nm pulse.Thus the temporal resolution is limited only by the pulsewidths of the 400 and 800 nm pulses.Careful control of the dispersion in each beam yields a time resolution of 45 fs, determined by measuring the upconverted Raman scattering of heptane.An example of one time-resolved measurement on neutral HBDI in methanol is shown in Fig. 2. As has been noted previously the decay is ultrafast and departs somewhat from a simple single exponential. 27BDI was prepared by standard methods, and the synthesis of FHBMI will be detailed elsewhere. 6All samples were made up to have an OD of roughly 0.5 at the excitation wavelength in a 1 mm pathlength cell (typically a concentration of 0.5 mM).
To construct the time-resolved emission spectra, uorescence decay proles were recorded across the emission spectrum by selecting the phase matching condition in the up-converting nonlinear crystal and tuning the wavelength of the SFG detection monochromator.These individual wavelength-resolved temporal proles are area normalized to the corresponding intensity in the steady state emission spectra to create a three dimensional time-intensity-wavelength surface.In a number of cases the range of detection wavelengths was limited by scattered pump radiation (a particular problem with short pulse broadband radiation).In these cases two-photon excitation was used, as has been described by Joo et al. 62 For these experiments a separate optical pathway for the pump beam was constructed, again with separate dispersion compensation.In this case the time resolution was measured as 60 fs.

Results and discussion
Steady state absorption and emission data for FHBMI are compared with those for HBDI in Fig. 3.In Table 1 the decay parameters for uorescence decays measured at the peak wavelength of the emission spectra are also compared.Both absorption and emission spectra for FHBMI are slightly red-shied from that of HBDI for both the neutral and anionic forms.The small shi in the anion spectrum on replacement of a weakly donating methyl group with the electron-withdrawing CHO group is in marked contrast to the large effect of exchanging a hydrogen for methyl on the phenolic ring reported previously. 11The relatively benign nature of the methyl for formyl substitution is conrmed by the uorescence decay times which are very similar for HBDI and FHBMI in both neutral and anion states.The similarity of HBDI and FHBMI in both charge states contrasts with the very different spectroscopic and photochemical properties of GFP and KFP. 63,64This is further evidence of the critical role played by protein structure and by specic interactions between the chromophore and adjacent amino acid residues in determining the photophysical properties of FPs.The most complete representation of excited state dynamics comes not from the single wavelength measurements reported in Table 1, but from observations of the temporal evolution of the entire emission spectrum.Such time-resolved where a ¼ 2g[ n À n p ]G.The amplitude is h, g the asymmetry parameter, G the width parameter, and n p the peak wavenumber.These log-normal parameters form the basis for further analysis of the time-resolved spectra.In Fig. 5 the time dependence of the mean frequency (rst moment) and integrated intensity of the log-normal ts are displayed.The integrated intensity represents the time-dependent excited state population (assuming there is no temporal evolution in the transition moment) and can be seen to be non-single exponential for both chromophores.The anionic states are longer lived than the neutral for both FHBMI and HBDI, consistent with what was observed in the single wavelength study (Table 1).Equally signicant is the observation that the mean frequency of the spectrum is essentially independent of time in all cases (Fig. 5a).A time-dependent shi in the spectrum is expected in (at least) two eventualities.First, if the excited electronic state has a signicantly different dipole moment to the ground state, then solvation dynamics should stabilise the new excited state, resulting in a time-resolved red shi in the emission spectrum occurring on the timescale of solvent reorganization. 65Second, if the excited state relaxes along a coordinate which causes closer approach of ground and excited states, and the transition moment remains nite, the emission is predicted to shi to the red with time.Thus, the fact that no time-dependent shi is observed places important constraints on the nature of the excited state dynamics.
The starting point for our simulation of the time-dependent spectra are the potential energy surfaces calculated by Olsen for the ground and excited states of HBDI and AHBMI, a molecule closely related to FHBMI, in their anionic states (which are typically the uorescent states in FPs).Calculations were made using two different sets of surfaces for the anion, one based on a state-averaged CASSCF model and the second on a perturbation theory-corrected CASSCF model (MR). 37he calculation was made as a function of the reaction coordinate which includes both double and single bond twist coordinates of the methine bridge (Fig. 1), varying between 0 and 90 , and these together create a single PES.We restrict our attention to chromophores which are initially in their cis ground state; we have shown elsewhere that cis and trans HBDI have indistinguishable excited state decay times, which accords with the very similar shapes of their calculated PES. 66hese PES are input to the generalized Smoluchowski eqn (2), which describes the diffusive evolution of a population density r(z,t) under a diffusion coefficient D (calculated from the Stokes-Einstein-Debye equation) along the calculated PES [S 1 (z)] where z is a generalised twist coordinate. 54g. 4 Time-resolved emission spectra (points) fit to the log-normal function (lines) for (a) neutral HBDI 50 fs to 800 fs (b) anionic HBDI 50 fs to 2ps (c) anionic FHBMI 50 fs to 2 ps.Solvent methanol, excitation wavelength 400 nm.
In previous approaches using (2) to simulate time-dependent spectra the generally unknown coordinate dependence of the excited state PES was calculated assuming coupling between harmonic bright and dark states. 52,54In the present case we digitized published PESs for HBDI and FHBMI anions and t them to polynomials for the calculation of r(z,t) using (2).The last term in eqn (2) describes nonradiative decay from the excited state to the ground state, where G(z) is a coordinate-dependent probability for decay back to the ground state (a 'sink' function).Several forms can be chosen for G(z) but in this work a pair of Gaussians function was chosen, described by: where k G1 and k G2 are the rate constants describing coupling of the ground and excited state, s 1 and s 2 are the widths and z 01 and z 02 are the centres, and we have set k G1 k G2 , s 1 ¼ s 2 and z 01 ¼ 1 and z 02 ¼ À1; the latter choice accords with the calculation of a conical intersection between or close approach of S 0 and S 1 at 90 twist.Consequently the solution to (2) for a particular choice of the diffusion coefficient, D, gives the population dynamics on the upper surface, r(z,t).The rate of population decay back to the ground state, however, depends on the overlap of the population with the sink functions (3).Finally, the excited state populations are converted to spectra for comparison with the measurements (Fig. 4 and 5), which requires the choice of a lineshape function g(n 0 (z),n(z) À n 0 (z)), where n 0 (z) is the coordinate-dependent spacing between the S 0 and S 1 surfaces.The lineshape function is here taken to be the log-normal lineshape.The nal requirement is a model for a coordinate-dependent transition moment, M 2 (z).These two combine with the calculated population density to give the uorescence and absorption spectra respectively: 52,54 I fl (t) f Ð dz g(n 0 (z),n(z) À n 0 (z))M 2 (z)r(z,t)n 3 (4) The M(z) in ( 4) and ( 5) accounts for the possibility that the transition moment may be a function of the twist coordinate.This has already been shown to be the case for a methylated HBDI derivative in which the substituents twist the methine bond away from the most stable planar state through steric crowding, giving rise to a weaker transition moment. 11To model the transition moment in (4) and ( 5) we use an arbitrary function which has the property of going to zero at the 90 twist angles.
The relevant coordinate-dependent parameters for the simulation are plotted in Fig. 6.
Fig. 6 shows that electronic excitation of the planar ground state generates a Franck-Condon excited state, in which the molecule is on a nearly at region of the PES and able to undergo barrierless or near barrierless twisting about one or both bridging bonds.These features are common to a number of other calculated PESs for HBDI.Excited state dynamics are typically obtained from these calculations by following a minimum energy pathway for the excited state, as it relaxes down the twisting and stretching coordinates.A consequence of this twisting is an increase in the energy of the ground state until at some point the two PESs approach one another or merge at a conical intersection.Once the PESs become sufficiently close it is assumed that internal conversion from excited to ground state occurs, consistent with the known ultrafast relaxation of HBDI in solution.
In line with this scenario we have calculated the excited state dynamics for the HBDI anion in methanol according to eqn (2) (the choice of D being discussed further below).Since the theoretical calculations described above typically include neither a coordinate dependence of the transition moment nor a rate for the radiationless decay, we initially set the second term in (2) to zero to focus on the dynamics of the population density.The resulting r(z,t), calculated timeresolved emission spectra and mean emission frequency are plotted in Fig. 7, using the M 2 (z) in Fig. 6.This procedure predicts a time-dependent spectral shi in the emission spectrum aer excitation, arising from the approach of ground and excited state surfaces with increasing z.This does not agree with the experimentally observed time independent spectrum (Fig. 5).It is shown below that the experimental result can only be reproduced when the simulation incorporates a coordinate dependence of both the transition moment and rate of decay to the ground state.
To simulate the spectra we assume a form for M 2 (z) (Fig. 6) and calculate r(z,t) adjusting the parameters of G(z) to obtain a best t to the calculated absorption and time-dependent emission spectra using ( 4) and ( 5).The results for HBDI are shown in Fig. 8.In this case a good t to the spectral shape and amplitude is obtained, but it does require a strong dependence of the rate of internal conversion on the twist coordinate (the actual example recovered from the t is shown in Fig. 6).Within this simple one dimensional representation such a coordinate dependence is essential to reproduce the time-independent mean frequency observed (Fig. 5): a narrower Gaussian (i.e. a 'pinhole sink') and a smaller k lead to a red shi in the spectrum with time.The results in Fig. 8 were obtained with a xed M 2 (z) and varying G(z).Clearly these two parameters can have compensating effects, and xing G(z) and varying M 2 (z) achieves similar (but less good) ts.These two cases do, however, correspond to physically different phenomena, and should be distinguished.In Fig. 8 decay is directly to the ground state, as a consequence of the second term in (2).The alternative corresponds to the build up of a population in the twisted state (analogous to the results seen in Fig. 7) which does not contribute to the emission due to a M 2 (z) which falls to zero rapidly with increasing |z|.These two cases are distinguishable by studying transient absorption spectra, to which a twisted dark state may be expected to contribute.We have measured transient infra-red spectroscopy of HBDI (neutral) in solution and nd that the excited state decays on a sub-picosecond timescale to yield a vibrationally hot ground state which cools in picoseconds.This is consistent with efficient internal conversion, i.e. a dominant role for kG(z) in the decay.
Fig. 8a shows the ts obtained for anionic HBDI using the MR surfaces with the hydrodynamic volume used for D (the choices for D being either a measured hydrodynamic volume or a calculated benzene ring volume, see below). 37,67In S1 the equivalent ts obtained using the volume of benzene for D and the CASSCF surfaces are shown.The time dependence of the integrated intensity and rst moment for all ts are shown in 8b and c.As taken from Olsen's paper, the S 0 CASSCF and MR surfaces differ by approximately 1300 cm À1 at 0 , while on the S 1 surface, the CASSCF surface is approximately 10700 cm À1 higher in energy.Once adjusted vertically to match the experimental absorption and t ¼ 50 fs resolved uorescence spectrum, there is no energy difference between the S 1 states at 0 , where the absorption takes place, but at 90 the atter MR surface is roughly 7000 cm À1 higher in energy.The surfaces provide reasonable ts to the experimental data at early times.Using the benzene volume provides a better t to the integrated intensity aer 1 ps, with the CASSCF surface giving the better t overall.The hydrodynamic volume is, however, better at earlier times, with little difference between the MR and CASSCF surfaces.None of the ts reproduce the small blue shi seen experimentally, but these variations are on the same scale as the precision of the calculated spectra.All ts required the broad (large s) and deep (large k) sinks shown in Fig. 6.In contrast to these results for HBDI, the same calculation failed for FHBMI.Since the experimental spectra do not appear to be signicantly different it is likely that the failure reects a less suitable PES for this molecule.This may arise from one of two causes.First, that there was some shortcoming in calculation for the more complex chromophore.Alternatively, the chromophore we have studied differs from that calculated by the exchange of a methyl for hydrogen (or a formyl for an acetyl substitution); this may have a signicant effect on the PES.
In Fig. 9 we compare the calculated potential surface for the HBDI anion with a recent calculation for the energy of the neutral form at 0 and 90 torsion about the double bond. 37The neutral form of HBDI presents a much more challenging case for accurate calculation than the anion, 45 explaining the relative paucity of theoretical results with which to model the data.To obtain a PES for simulation of the neutral HBDI (Fig. 9a) the anionic surface has been scaled to match the two neutral points at 0 and 90 .A gradient was calculated given by the energy difference between the anionic surface and the two points for the neutral energy.This gradient was applied to the anionic surface for the single bond torsion resulting in a calculated coordinate dependent neutral surface connecting the two calculated points.It can be seen that the calculation apparently predicts a atter PES for the neutral HBDI than for the anion.The constructed PES were vertically shied to simulate the observed absorption spectrum and then used to calculate the time-dependent uorescence spectra, assuming the form of M(z) shown in Fig. 6.In Fig. 9b the calculated time dependence of the integrated area predicts a decay for the neutral which is slower than that of the anion, which is opposite to the experimental observations.This underlines the difficulty of accurately calculating the excited state surface for neutral HBDI.This difficulty sets a signicant new challenge for theoretical calculations, since the neutral state of Fig. 8 (a) Comparison of the spectra calculated using eqn (5) (red) and the log normal fits to the data (black), eqn (1), for the HBDI anion in methanol using the hydrodynamic volume.(b) and (c) Calculated integrated intensity and first moment of the spectra using the volume of benzene (circles) and the hydrodynamic volume (diamonds) for D in (2) (see text for details) for the MR surfaces (black) and the CASSCF surfaces (red).Data are for HBDI (black crosses).
HBDI is the one most oen responsible for the excited state reactions which are critical in second generation uorescent proteins. 68,69inally we turn to the solvent friction effects on the excited state decay.All of the data above have been presented in methanol solution and the diffusion coefficient has been calculated from the Stokes-Einstein-Debye relation: where h is the solvent viscosity and V is the volume, where the volume is calculated from the estimated volume of the isomerizing unit.This was estimated in two ways, rst from the tabulated volume of a benzene ring and second from the hydrodynamic volume of benzonitrile, measured elsewhere. 70The results for both volumes are in qualitative agreement with the experimental data (Fig. 8a, b and c) in methanol.This suggests that eqn ( 7) is an appropriate model for the friction in methanol.However, it is worth noting that the recovery times in methanol are quite close to those observed in the gas phase, 47 suggesting that the rate in this very uid solvent may be determined by the shape of the PES rather than friction of the medium.
To test this friction model further, the excited state decay times of the HBDI in a series of solvent of varying viscosity were measured (Fig. 10 and Table 2).These measurements were made at two or three wavelengths, which was sufficient to check that the more viscous solvents do not introduce a wavelength dependence (which would in turn indicate a time-dependent spectrum).The excited state decays remained independent of wavelength.The measured ultrafast decay is seen to be only slightly dependent on the solvent viscosity (Fig. 10a) in line with previous observations suggesting a volume-conserving (thus friction independent) reaction coordinate. 11,26,27In contrast, when we apply the corresponding diffusion coefficient (7) in the calculation a very strong viscosity dependence is predicted, dramatically stronger than that observed experimentally (Fig. 10b).This suggests a signicant shortcoming in the model for the diffusion in eqn (2), which probably arises because of the signicant friction experienced for motion along the coordinates recovered in the calculations.

Conclusions
Ultrafast uorescence spectroscopy of the chromophore of GFP (HBDI) and a molecule closely related to the chromophore of KFP (FHBMI) have been measured in their neutral and anionic states.The spectra and dynamics of the two chromophores in solution are similar, in contrast to the behaviour of the two proteins.This highlights the importance of including the protein matrix in determining the photophysical properties of FPs.The time-resolved uorescence spectra were measured with a 50 fs time resolution.For both molecules and for neutral and anion states the spectral mean frequency is independent of time.
The excited state dynamics of the anionic form of HBDI have been simulated by solving the Smoluchowski equation for diffusion on a potential energy surface.The inputs for the simulations were potential energy surfaces calculated for the chromophore. 37,67The simulation for HBDI accurately reproduced the experimentally observed dynamics in methanol provided a strong coordinate dependence of the radiationless decay rate and transition moment were taken into account.This suggests that the simplest models for excited state dynamics in FP chromophores, relaxation on an excited state PES to a twisted geometry followed by internal conversion to the ground state at that geometry (i.e. at a 'pinhole sink' or conical intersection), are inadequate.To reproduce the experimental behaviour theoretical models need to include the possibility of coordinate-dependent rates of radiationless decay, with decay occurring rather early on along the coordinate.An attempt to simulate the faster dynamics observed in the neutral form of the chromophore on the basis of the very limited data for its excited state PES were not successful, suggesting the need for more detailed calculation on the neutral chromophore, which is of particular importance in photochemically active FPs.
Although the model calculations using a realistic PES were successful for the HBDI anion in methanol they completely fail to correctly predict the observed very weak effect of solvent viscosity on lifetime.As suggested elsewhere the observed weak dependence of the HBDI decay time on solvent viscosity is consistent with a volume-conserving isomerization coordinate, i.e. one which displaces a much smaller volume of solvent than either of the single bond rotations which generated the PESs employed here. 26,27The possibility of a volume-conserving 'hula twist' coordinate promoting radiationless decay has been considered, but calculations show that this is energetically up-hill, which is inconsistent with the observed ultrafast barrierless nature of the internal conversion. 34An alternative volume conserving radiationless decay coordinate, which was reported to be important in at least some calculations of the excited state dynamics of HBDI, is pyramidalization at the central bridging carbon atom. 35,38,40Such a shape change appears to be sufficiently small to be consistent with the observed weak viscosity dependence.Relaxation along a pyramidalization coordinate has been invoked to account for the ultrafast excited state double bond rotation observed in a family of synthetic molecular motors. 71,72hese considerations of the very efficient radiationless decay in HBDI and FHBMI still leave open the question of how the protein tunes the quantum yield of the chromophore between >0.8 and <0.01 depending on the environment.As has been discussed extensively, the electronic spectra of HBDI are shied signicantly on incorporation into GFP. 68This suggests the possibility that the interaction between HBDI and the protein leads to changes in electronic structure, such that the surfaces calculated for the gas phase may not be an accurate representation of the chromophore in the protein.Alternatively protein -chromophore interactions may provide the friction required to resist relaxation along the barrierless PES.It is difficult to see how simple steric crowding in the protein could resist a single bond rotation coordinate, given that there is sufficient room in the cavity of the bbarrel to allow at least one bond to rotate. 73This is especially true as full rotation is evidently not required to promote radiationless decay (Fig. 6).It is yet more difficult to see how such steric crowding could prevent isomerization along a volume-conserving coordinate, which is the nature of the coordinate most consistent with the observed viscosity dependence.We speculate that a motion which may be restricted in the protein, but not in free solution, is the length of the chromophore.Both twisting and pyramidalization will lead to some modulation of the length of the chromophore.At least in the highly uorescent anionic state of wild type GFP the chromophore is strongly H-bonded at both the phenolic and imidazole rings.The stretching of these H-bonds as a consequence of the excited state dynamics may be sufficient to generate the friction required (or equivalently generate a sufficiently high energy barrier) to suppress the excited state structural relaxation leading to rapid internal conversion.This further highlights the importance of including proteinchromophore interactions when simulating the excited state dynamics of FPs.

Fig. 1
Fig. 1 The chemical structures of HBDI and FBHMI in their neutral states.The anion states are deprotonated at the phenol OH.

Fig. 2
Fig. 2 Example of the measured fluorescence decay of neutral HBDI (black).The up-conversion of the solvent Raman signal scattered from the pump pulse is shown (red).This determines the time resolution of the experiment.

Fig. 5
Fig. 5 Time dependence of the (a) first moment and (b) integrated intensity for HBDI (black) and FBHMI (blue) recovered from the log normal plots.Open symbols neutral, closed anion.

Fig. 6
Fig. 6 Parameters used in the simulations.Solid and dashed black lines are the ground and excited state PESs calculated by Olsen.The red circles are polynomial fits to the PES used in the simulation.The green line is the Gaussian decay functions and the dashed blue line is the coordinate dependent transition moment.

Fig. 7
Fig. 7 (a) the r(z,t) calculated from (2) with the final term omitted (no decay function) The population decays at Franck-Condon position z ¼ 0 and builds up at z ¼ AE1.(b) The spectra calculated for the data in (a) using (4) and (5).

Fig. 9
Fig. 9 (a) Anionic HBDI S 0 and S 1 surfaces calculated by Olsen (black).Calculated S 0 and S 1 neutral HBDI points are shown (black crosses).Neutral surfaces were created by scaling the anionic surfaces to the neutral points (red) (see text for details).(b) Time dependence of the integrated intensity of HBDI predicted by the neutral surface (open circles) and Olsen's anionic surface (filled circles).

Fig. 10 (
Fig. 10 (a) Time-resolved fluorescence decay for HBDI neutral in a range of solvent.(b) Calculated integrated intensity from (2) where D for the appropriate solvent was used and other parameters fixed as for methanol.

Table 1
Results of fits of the fluorescence decay to a sum of two exponentials.The weights are a i , the decay times s i and the weighted mean decay time <s> are readily obtained from wavelength-intensity slices of the intensitytime -energy surface described above.Examples are shown for HBDI and FHBMI in methanol in Fig.4.The data points are well represented by a t to the lognormal function: Paper Faraday Discussions Open Access Article.Published on 01 March 2013.Downloaded on 08/06/2017 10:42:29.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.spectra

Table 2
Solvent dependence of fits of the fluorescence decay to a sum of two exponentials.The weights are a i , the decay times s i and the weighted mean decay time <s>