Michael A.
Parkes
,
Jonathan
Crellin
,
Alice
Henley
and
Helen H.
Fielding
*
Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, UK. E-mail: h.h.fielding@ucl.ac.uk
First published on 21st May 2018
Indole is an important molecular motif in many biological molecules and exists in its deprotonated anionic form in the cyan fluorescent protein, an analogue of green fluorescent protein. However, the electronic structure of the deprotonated indole anion has been relatively unexplored. Here, we use a combination of anion photoelectron velocity-map imaging measurements and quantum chemistry calculations to probe the electronic structure of the deprotonated indole anion. We report vertical detachment energies (VDEs) of 2.45 ± 0.05 eV and 3.20 ± 0.05 eV, respectively. The value for D0 is in agreement with recent high-resolution measurements whereas the value for D1 is a new measurement. We find that the first electronically excited singlet state of the anion, S1(ππ*), lies above the VDE and has shape resonance character with respect to the D0 detachment continuum and Feshbach resonance character with respect to the D1 continuum.
![]() | ||
Fig. 1 Structure of the deprotonated cyan fluorescent protein chromophore with the deprotonated indole moiety in black. |
Experimentally, the most direct way of probing the electronic structure and non-radiative relaxation dynamics of molecules is photoelectron spectroscopy. In earlier work, we proposed that the electronic structure of GFP and CFP chromophores were governed by the deprotonated phenol and indole molecular units respectively.5 Recently, Anstöter et al. developed this idea further by suggesting that the basic electronic structure of photoactive chromophores can be understood by using Hückel theory to construct the molecular orbitals of a photoactive protein chromophore from linear combinations of the molecular orbitals of their constituent molecular building blocks.6 This implies that to understand the electronic structure and non-radiative relaxation dynamics of CFP requires a detailed understanding of the electronic structure of the deprotonated indole anion building block. The only gas-phase studies of the deprotonated indole anion reported to date are a low-resolution photoelectron spectroscopy study that determined the vertical detachment energy (VDE) to be 2.31 ± 0.15 eV5 and a recent high-resolution slow velocity-map imaging (SEVI) study that reported a very accurate value for the adiabatic detachment energy (ADE) of 2.4315 ± 0.0017 eV.7 Photodetachment spectroscopy has also been employed to study deprotonated tryptophan-containing peptides and their corresponding neutral radicals.8,9 Here, we report the results of a new, combined photoelectron imaging and computational chemistry study of the deprotonated indole anion. We report values for the VDEs to both the ground and first electronically excited states of the neutral deprotonated indole radical and we determine the vertical excitation energy and resonance character of the first electronically excited state of the deprotonated indole anion.
The pBasex inversion method was used to obtain photoelectron velocity and angular distributions.11 Electron kinetic energy (eKE) spectra were obtained by calibrating the radial photoelectron velocity distribution against the photoelectron spectrum of iodide12 and the energy resolution, ΔE/E < 5%. Qualitative insight into the nature of the molecular orbitals from which the electrons are detached can be obtained from the photoelectron angular distributions (PADs) which, for a one-photon process, is given by
I(θ) ∝ 1 + β2P2(cos![]() | (1) |
Vertical detachment energies (VDEs) were calculated using various methods. We used electron propagator theory (EPT),20 a method that has been found to yield reliable VDEs in previous work.21–24 We also used the equation-of-motion coupled-cluster method with single and double excitations for the calculation of ionisation potentials (EOM-IP-CCSD).25 Both of these methods allowed us to calculate the VDE to the first electronically excited state of the neutral radical, D1, as well as that to the electronic ground state of the neutral radical, D0. We also used DFT to calculate the VDE as the energy difference between the anion and neutral radical at the optimised geometry of the anion. Calculating the VDE this way allowed us to make a direct comparison with the adiabatic detachment energy (ADE) calculated as the difference between the lowest vibrational state of the anion and the lowest vibrational state of the neutral radical (0–0 transition).
Vertical excitation energies (VEEs) of the excited singlet states of the anion were calculated using time-dependent DFT (TD-DFT)26 with the CAM-B3LYP functional and 6-311++G(3df,3dp) and aug-cc-pVTZ basis sets. They were also calculated using the ADC(2)/6-311++G(2d,2p) method.27,28 Due to the number of diffuse functions in the basis set, a huge number of continuum states were calculated. Thus, to assist with our interpretation of the ADC(2) calculations, we first used the 6-311G(2d,2p) basis set, which did not give accurate VEEs because it neglected the interaction of the resonance states with the continuum but it allowed us to determine the nature of the molecular orbitals involved in the transitions. Having identified the key molecular orbitals, we then added the diffuse functions to obtain more accurate VEEs. We chose to use the ADC(2) method having recently found it compared favourably with high-level SA-CASSCF(14,12)-PT2/cc-pVDZ calculations,29 for a series of photoactive yellow protein chromophore anions,21 and high-level XMCQDPT2/(aug)-cc-pVTZ level calculations,30 for a series of GFP chromophore anions.22
Geometry optimisations, vibrational frequencies, TD-DFT and EPT calculations were performed using the Gaussian09 software suite,31 EOM-IP-CCSD and ADC(2) calculations were carried out using the Q-Chem software package.32
The photoelectron spectrum was calculated using ezSpectrum (version 3.0).33 These calculations require the equilibrium geometries, harmonic frequencies and normal mode vectors of the deprotonated indole anions and corresponding neutral radicals as input and we used those obtained using the CAM-B3LYP/aug-cc-pVTZ method. The Franck–Condon overlap integrals were calculated in ezSpectrum using the parallel normal mode approximation. The vibrational temperature of the anions was assumed to be 300 K and the minimum intensity threshold was set to 0.001. The maximum number of vibrational quanta in the anion and neutral radical were limited to 3 and 6, respectively. The resulting stick spectra were convoluted with Gaussian instrumental profiles with full-width at half maxima (FWHM) equivalent to the instrumental resolution, ΔE/E = 4% at 0.4 eV eKE.
PADs for direct detachment to D0 and D1 were calculated using ezDyson (version 3.3).34,35 The Dyson orbitals were obtained from the EOM-IP-CCSD/aug-cc-pVDZ calculation and represent the overlap between the N electron molecular wavefunction of the anion in its ground electronic state and the N − 1 electron molecular wavefunction of the corresponding neutral radical. The PADs were calculated in ezDyson using the plane wave approximation.
All the photoelectron spectra have a relatively narrow peak centered at 2.45 ± 0.05 eV eBE, a broad unresolved feature at higher eBE and a peak centered around 3.20 ± 0.05 eV eBE. The peaks around 2.45 eV and 3.20 eV eBE remain at constant eBE for all photon energies, characteristic of direct photodetachment (PD). Between the two peaks at 2.45 eV and 3.20 eV is a feature that shifts to higher eBE with increasing photon energy, characteristic of an indirect PD process. More specifically, this feature is centered around 0.9–1.0 eV eKE and remains at constant eKE with increasing photon energy. The relative intensities of the first two peaks in the photoelectron spectra vary with increasing photon energy, another indication that indirect PD is contributing to the photoelectron spectrum.
The calculated VDEs for D0–S0 and D1–S0 processes are presented in Table 1. The calculated values for detachment to the D0 continuum all lie in the range 2.4–2.6 eV and are in excellent agreement with our measured value of 2.45 ± 0.05 eV eBE. The VDEs are also very close to the CAM-B3LYP/aug-cc-pVTZ calculated 0–0 transition at 2.46 eV, suggesting that the geometry of the anion is similar to the geometry of the neutral radical in its ground electronic state. The calculated values for detachment to the D1 continuum are very close to our experimental measurement of 3.20 eV eBE indicating that this feature arises from direct detachment to the continuum associated with the first electronically excited state of the neutral radical of deprotonated indole.
Method | D0 | D1 |
---|---|---|
B3LYP/6-311++G(3df,3dp) | 2.36 | |
B3LYP/aug-cc-pVDZ | 2.36 | |
B3LYP/aug-cc-pVTZ | 2.44 | |
CAM-B3LYP/6-311++G(3df,3dp) | 2.47 | |
CAM-B3LYP/aug-cc-pVDZ | 2.45 | |
CAM-B3LYP/aug-cc-pVTZ | 2.46 | |
EPT/6-311++G(3df,3dp) | 2.56 | 3.11 |
EOM-IP-CCSD/aug-cc-pVDZ | 2.37 | 3.00 |
Experiment | 2.45 ± 0.05 | 3.20 ± 0.05 |
To assist with our interpretation of the structure in the photoelectron spectra, we have calculated the photoelectron spectrum corresponding to direct photodetachment from the S0 state of the deprotonated indole anion to the D0 state of the neutral radical. The photoelectron stick spectrum is convoluted with a Gaussian instrument function with FWHM ΔE/E = 4% and presented in Fig. 3 together with the experimental 359 nm and 346 nm photoelectron spectra. The most intense line in the stick spectrum is the 0–0 transition at 2.46 eV. The calculated spectrum allows us to identify the dominant features in the photoelectron spectra as a progression in the ν29 mode (1468 cm−1) of the neutral radical, an in-plane vibration (see inset of Fig. 3). Despite the excellent agreement between experiment and theory for the sharp peak at low eBE, there is additional photoelectron signal in the experimental spectra that is not accounted for in the calculation and the profile of this additional photoelectron signal varies with photon energy. Such variation in photoelectron signal is observed at all wavelengths (Fig. 2) and is consistent with indirect PD following photoexcitation of excited electronic states of the anion.
The β2-parameters plotted in Fig. 2 are fairly constant across the photoelectron spectra and slightly negative, as might be expected for photodetachment from a π-orbital. We have calculated the β2-parameters that characterise the PADs corresponding to direct detachment from the S0 state of the deprotonated indole anion to the D0 and D1 states of the neutral radical. The calculated β2-parameters are plotted in Fig. 4, together with the experimental values. The experimental values are determined as intensity-weighted averages from regions of the photoelectron spectra corresponding to eBEs of 2.4–2.5 eV for D0 and 3.0–3.4 eV for D1, with appropriate intensity-averaged standard deviations. The experimental values are plotted as a function of eKE above the experimentally determined D0 and D1 VDEs (2.45 eV and 3.20 eV). The calculated and measured PADs are very similar to one another for both D0 and D1 detachment channels.
To identify the excited electronic state of the anion embedded in the detachment continuum, we calculated the VEEs of electronically excited singlet states of the deprotonated indole anion. The VEEs are summarised in Table 2 and suggest that there is a weak, but not insignificant, πHOMO → π* transition with VEE around 3.8–4.0 eV, which is accessible with the photon energies employed in this work (3.5–4.0 eV). The ADC(2)/6-311G(2d,2p) molecular orbitals (MOs) involved in this transition are presented in Fig. 5. Both the highest occupied molecular orbital (HOMO) and π* state are delocalised across the whole molecular anion.
Method | VEE (f) | Main configuration |
---|---|---|
CAM-B3LYP/aug-cc-pVTZ | 3.97 (0.029) | 0.55 (πHOMO → π*) |
CAM-B3LYP/6-311++G(3df,3dp) | 4.05 (0.032) | 0.67 (πHOMO → π*) |
ADC(2)/6-311++G(2d,2p) | 3.84 (0.029) | 0.44 (πHOMO → π*) |
In order to understand the indirect PD processes contributing to the broadening of the photoelectron spectra of deprotonated indole following photoexcitation in the range 359–310 nm (3.45–4.00 eV), it is useful to consider the main electronic configurations and molecular orbitals of the states involved; these are illustrated in a Jablonski diagram (Fig. 5). Direct PD from S0 to D0 and D1 involve removal of an electron from the HOMO and the HOMO−1, respectively. Photoexcitation of S1 involves excitation of a single electron from the HOMO to the π* molecular orbital and our calculations suggest that the VEE is around 3.8–4.0 eV. For the S1 state of the deprotonated indole anion to be accessible in our experiments, the relaxation energy in S1 would need to be 0.35–0.55 eV. The relaxation energy of the lowest lying 1ππ* state of indole has been calculated to be around 0.5 eV, so this is not at all unreasonable.36 The S1 state has a hole in the HOMO so has shape resonance character with respect to the D0 continuum but has Feshbach resonance character with respect to the D1 continuum. Thus, we expect the S1 state to couple strongly to the D0 continuum and more weakly to the D1 continuum and, therefore, autodetachment from S1 to occur predominantly to D0. The broadening of the D0–S0 photoelectron spectrum can be attributed to this autodetachment from S1 to D0. Despite the clear evidence for a resonance in the continuum, the experimental and calculated PADs are very similar (Fig. 4); this is interesting, since resonances can cause changes in PADs.6 We do not see any evidence of thermionic emission from the ground electronic state of the anion, characterised by a low eKE (high eBE) signal with exponential profile, I(eKE) ∝ exp(−eKE/kBT),37 which suggests that autodetachment from S1 to D0 is faster than internal conversion from S1 to S0.
Our original motivation for this work was to improve our understanding of the role of the deprotonated indole anion as a molecular building block of the CFP chromophore. In this respect, an important measurement is that of the VEE of deprotonated indole, 3.8–4.0 eV, which is higher than that of deprotonated phenol, 3.5 eV.6 If we consider the MOs of the excited states of the photoactive protein chromophore to be linear combinations of the excited states of their constituent molecular building blocks,6 replacing the deprotonated phenol anion with a deprotonated indole anion raises the energy of the π* MO of the substituent relative to the π* MO of the 2,3-dimethyl-5-methylene-3,5-dihydro-4H-imidazol-4-one molecular unit (grey in Fig. 1). Simple Hückel theory then predicts that the first electronically excited state of the deprotonated CFP chromophore will be higher than that of the deprotonated GFP chromophore, which is in agreement with gas-phase action absorption spectroscopy measurements of the deprotonated GFP (2.53 eV)38 and CFP (2.70 eV)4 chromophores. Moreover, the higher lying excited state of the deprotonated CFP chromophore will be localised on the indole moiety and be lower in energy than that of the deprotonated GFP chromophore. Following photoexcitation of S1, GFP has the ability to act as a light-induced electron donor in the presence of an oxidant,39,40 which is valuable for monitoring redox processes in cells. Following UV photoexcitation of higher lying electronically excited states, GFP has been shown to undergo decarboxylation41 and it has been proposed that the primary step might involve electron emission from excited electronic states.42 Thus, since simple Hückel theory predicts that the excited states of the CFP chromophore are similar to those of the GFP chromophore, it seems reasonable to suppose that similar processes involving electron emission from electronically excited states will occur in CFP and that changing the substituent provides a means of tuning these photoinduced processes.
Footnote |
† Electronic supplementary information (ESI) available: Optimised anion and radical geometries; vibrational frequencies and FCFs; ezSpectrum calculations; vertical excited energies and characters; molecular orbitals. See DOI: 10.1039/c8cp01902a |
This journal is © the Owner Societies 2018 |