Michael A.
Parkes
,
Ciara
Phillips
,
Michael J.
Porter
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 24th March 2016
Understanding how the interactions between a chromophore and its surrounding protein control the function of a photoactive protein remains a challenge. Here, we present the results of photoelectron spectroscopy measurements and quantum chemistry calculations aimed at investigating how substitution at the coumaryl tail of the photoactive yellow protein chromophore controls competing relaxation pathways following photoexcitation of isolated chromophores in the gas phase with ultraviolet light in the range 350–315 nm. The photoelectron spectra are dominated by electrons resulting from direct detachment and fast detachment from the 21ππ* state but also have a low electron kinetic energy component arising from autodetachment from lower lying electronically excited states or thermionic emission from the electronic ground state. We find that substituting the hydrogen atom of the carboxylic acid group with a methyl group lowers the threshold for electron detachment but has very little effect on the competition between the different relaxation pathways, whereas substituting with a thioester group raises the threshold for electron detachment and appears to ‘turn off’ the competing electron emission processes from lower lying electronically excited states. This has potential implications in terms of tuning the light-induced electron donor properties of photoactive yellow protein.
Fig. 1 Left: PYP chromophore in its protein environment (Protein Data Bank). Right: Structures of the deprotonated PYP chromophores employed in this work. |
There have been numerous experimental investigations of the spectroscopy and excited state dynamics of isolated model PYP chromophores, in the gas-phase and in solution, aimed at disentangling the roles of the chromophore and the protein.9–27 Action spectra of the deprotonated coumaric acid anion in the gas-phase have been found to have a maximum around 430 nm, which is reasonably close to the absorption maximum of PYP (446 nm)21 and it is now generally accepted that this similarity is the combined effect of numerous chromophore-protein interactions. Although the absorption spectra of the gas-phase chromophore and PYP are similar, time-resolved measurements of a deprotonated methyl terminated ketone analogue of the PYP chromophore in the gas-phase, excited at 400 nm (3.1 eV), revealed that 80% of the excited state population relaxed back to the ground state on a timescale of 52 ps via a twisted intermediate, similar to isolated chromophores in solution, and that the remaining 20% of the excited state population underwent autodetachment (electron emission).16 From this, the authors concluded that one of the roles of the protein environment was to funnel the excited state population through the S1/S0 conical intersection seam to the cis isomer and to suppress electron emission. It is worth noting that electron emission has also been observed in isolated chromophores in aqueous solution12 and in PYP28 at higher photoexcitation energies. Interestingly, the spectrum of the electron emitted from the chromophore in PYP was found to be similar to that of a solvated electron, suggesting that the protein pocket provides a local environment for an electron that is similar to water.
Studies of numerous chromophore analogues in aqueous solution have found that the substituent on the coumaryl tail of the chromophore plays a key role in the excited state dynamics, highlighting the significance of the thioester link to the Cys69 amino acid residue. These studies found that derivatives with strong electron accepting groups, such as thioesters, return to the initial trans configuration on timescales of ∼1 ps, whereas those with weak electron accepting groups, such as amides or carboxylate, have longer excited state lifetimes, up to 10 ps for the carboxylate, and measurable trans–cis isomerisation yields.9,11,12,18
From a theoretical perspective, there have been numerous studies of PYP and the PYP chromophore, in the gas-phase and in solution, aimed at investigating the role of the environment on its electronic structure and dynamics.29–37 Various strategies have been employed, ranging from studies of a few solvent molecules or the closest protein residues to quantum-mechanical/molecular mechanics and averaged solvent electrostatic potential/molecular dynamics methods. Such studies have highlighted the importance of torsions around the C–C single bond between the phenoxide ring and the CC double bond in controlling the trans–cis isomerisation process in the protein.29 They have also rationalised the blue solvatochromic shift of pCA− and its derivatives in terms of a transfer of electron density from the phenolate end of the chromophore to the rest of the chromophore, resulting in a decrease in dipole moment between the ground and first electronically excited states.34 García-Prieto and coworkers have recently reported a detailed theoretical study of the absorption spectra of different chromophore analogues in the gas-phase and solution and found that the presence of a sulfur atom on the coumaryl tail modulates the solvent effect for the first few excited electronic states37 which could explain the differences observed in the excited state dynamics in aqueous solution discussed above.9,11,12,18
In a previous anion photoelectron spectroscopy study of gas-phase pCA− and its ortho- and meta-isomers, we found that moving the position of the O− group on the phenoxide group of the chromophore controlled the competition between electron emission and internal conversion.27 The focus of this paper is to determine whether the thioester link also plays a role in controlling the competition between electron emission and internal conversion in the gas-phase. We employ a combination of anion photoelectron spectroscopy and quantum chemistry calculations to investigate the effect of substituting the hydrogen atom of the carboxylic acid group with a methyl group and substituting the oxygen atom of the methoxy group with a sulfur atom. Specifically, we excite pCA−, its methyl ester, pCE−, and its methyl thioester, pCT−, (Fig. 1) with ultraviolet light in the range 350–315 nm (3.54–3.94 eV), within the 21ππ* ← S0 absorption band. Improving our understanding of the role of the thioester link in controlling electron emission in PYP may provide a basis for practical applications such as the rational design of photoactive materials with specific redox properties.
Vertical excitation energies (VEEs) of the singlet excited electronic states of the phenolate deprotonated coumaric acid chromophores were calculated using the CAM-B3LYP/6-311++G(3df,3pd) method. The long-range corrected version of B3LYP using the Coulomb Attenuating Method (CAM)50 was chosen for its potential to describe excited states with charge-transfer character and because it has been shown to yield reasonable results for the phenolate form of the pCA− anion.22
Vertical detachment energies (VDEs) were determined using electron propagator theory (EPT)51,52 and a 6-311++G(3df,3pd) basis set. We have already benchmarked this approach against the high-resolution photoelectron spectrum of the phenoxide anion and shown that it yields good results for the deprotonated green fluorescent protein chromophore anion.40 Adiabatic detachment energies (ADEs) were determined using the B3LYP/6-311++G(3df,3pd) method, ADE = Eradicalmin + EradicalZPE − [Eanionmin + EanionZPE], where subscripts ‘min’ and ‘ZPE’ represent the minimum energy and zero-point energies, respectively. Higher VDEs were determined by calculating the VEEs of the first electronically excited doublet states of the neutral radicals at the optimised geometries of the corresponding anions using the CAM-B3LYP/6-311++G(3df,3pd) method and adding these to the VDEs calculated using the EPT method.
Fig. 2 Photoelectron spectra of pCA−, pCE− and pCT− recorded at 350 nm (3.54 eV), 320 nm (3.87 eV) and 315 nm (3.94 eV) (solid black lines). The 350 nm photoelectron spectra (shaded grey) are superimposed on photoelectron spectra recorded at shorter wavelengths, whose intensities have been scaled so the rising edges are aligned. Dashed black lines mark the calculated VDEs (Table 1) and solid red lines mark the calculated VEEs determined by García-Prieto and coworkers for the two lowest energy singlet ππ* states,37 listed in Table 2. |
All the photoelectron spectra are dominated by broad, unresolved features at low eBEs. The rising edges of these broad features remain at constant eBE for all photon energies and are therefore attributed to direct photodetachment (PD). The maxima of the 350 nm (3.54 eV) spectra are in good agreement (∼0.1 eV) with the calculated vertical detachment energies (VDEs) and adiabatic detachment energies (ADEs), which are presented in Table 1. It is worth noting that there is significant photoelectron emission at lower eBEs than the calculated VDEs and ADEs in all the photoelectron spectra. In pCA−, and one of its isomers oCA−, this has been attributed to rotation around the single bond between the CC double bond and the phenoxide group, so we assume these torsional motions play a similar role in pCE− and pCT−. The fact that the calculated VDEs and ADEs are so close to one another suggests that the minimum energy geometries of the anion and the radical are very similar, in agreement with the results of our combined photoelectron spectroscopy and computational study of gas-phase pCA− and its ortho- and meta-isomers.27 Our calculated VDEs are also in good agreement with the EOM-IP-CCSD/6-311+G(df,pd) value reported by Zuev et al. for pCA− (2.92 eV),53 the OVGF/aug-cc-pVDZ value reported by Gromov et al. for pCT− (2.90 eV)30 and the experimental VDE reported for the methyl ketone analogue (2.9 eV).16
Chromophore | VDE | ADE (0–0) | Expt |
---|---|---|---|
pCA− | 3.00 (0.879) | 2.99 | 2.91 ± 0.05 |
pCE− | 2.91 (0.875) | 2.91 | 2.83 ± 0.05 |
pCT− | 3.07 (0.879) | 3.03 | 2.96 ± 0.05 |
It is worth noting that pCA− can be formed during electrospray ionisation as a carboxylate anion or as a phenolate anion, whereas pCE− and pCT− have only one deprotonation site and can only be formed as phenolate anions. It is clear from the photoelectron spectra (Fig. 2) that all three chromophores have similar VDEs (around 3 eV), which supports our earlier suggestion that pCA− is formed in its phenolate form in our instrument when using methanol with a few drops of aqueous ammonia as a solvent.27 This contrasts with the observations of Almasian et al.,25 but perhaps emphasises the importance of instrumental parameters, such as the position of the electrospray head, as well as choice of solvent in determining the deprotonation site.54 The VDE for the carboxylate form of pCA− (4.68 eV)27 is significantly higher than for the phenolate form so it is also possible that both phenolate and carboxylate forms are present in our instrument but that we are only sensitive to the phenolate form in experiments with photons <4 eV.
Although the VDEs of the chromophores are similar, the maxima in the photoelectron spectra and calculated VDEs (Table 1) increase in the order pCE− < pCA− < pCT−. This trend can be understood in terms of the stabilising effect of electron accepting substituents lowering the energy of the resonantly stabilised anions as the electron affinity increases in the order OMe < OH < SMe.
The calculated vertical excitation energies for the first two 1ππ* states and the first 1nπ* states of pCA−, pCE− and pCT− are listed in Table 2, alongside the results of higher-level SA-CASSCF(14,12)-PT2/cc-pVDZ calculations by García-Prieto and coworkers37 and experimental values reported by Rocha-Rinza and coworkers.21 The molecular orbitals involved in the transitions are plotted in Fig. 3. Our calculated value for pCE− is consistent with that calculated using the CAM-B3LYP/aug-cc-pVTZ method22 and our value for pCA− lies within the range of values calculated by Zuev et al.53 and Uppsten and Durbeej,55 although all the values are ∼0.5 eV higher than the experimental values determined from action spectra21 and those calculated using the SA-CASSCF(14,12)-PT2/cc-pVDZ method.37 Nonetheless, the characters of the excited states calculated using the CAM-B3LYP/6-311++G(3df,3pd) method, and the energies of the excited states with respect to one another, are in good agreement with those calculated using the SA-CASSCF(14,12)-PT2/cc-pVDZ method. To guide our interpretation of the photoelectron spectra, we have marked the more accurate VEEs determined by García-Prieto and coworkers using the SA-CASSCF(14,12)-PT2/cc-pVDZ method37 on our experimental spectra (Fig. 2) and use the configurations obtained from the CAM-B3LYP/6-311++G(3df,3pd) to determine the characters of the excited states with respect to the detachment continua.
The effect of substituting the hydrogen on the carboxylic acid group in pCA− for a methyl group (pCE−) has very little effect on any of the VEEs. However, substituting the oxygen atom for a sulfur atom (pCT−) causes the 11ππ* and 11nπ* states to red-shift by 0.1–0.2 eV, although it has very little effect on the 21ππ* state. These observations can be understood in terms of the molecular orbitals involved in the transitions (Fig. 3). Transitions to the 11ππ* and 11nπ* states are to the π1* molecular orbital, which is delocalised across the anion and therefore stabilised by the electron accepting methyl thioester group, whereas the 21ππ* transition is to the π2* molecular orbital, which is localised on the phenoxide group and is barely influenced by changing the substituents on the coumaryl tail.
It is clear from Fig. 2 that the 11ππ* and D0 states are very close to one another in pCA− and pCE−, but that the 11ππ* state is 0.2–0.3 eV lower in energy than the D0 state in pCT−. This suggests that the thioester link between the chromophore and the protein plays a role in ensuring that the 11ππ* state is bound with respect to photodetachment following excitation at the maximum of 1ππ* ← S0 absorption band. This is consistent with the fact that there have not been any reports of electron emission from PYP following excitation within the 1ππ* ← S0 absorption maximum at 446 nm (2.78 eV), even though electron emission has been observed following excitation at higher photon energies.28
As the photon energy increases, the broad features in the photoelectron spectra change shape on the high eBE side, characteristic of an indirect PD process following resonant excitation of a higher lying electronically excited state (autodetachment). In pCA−, and one of its isomers mCA−, this broadening has already been rationalised in terms of resonant excitation of higher-lying 21ππ* excited states followed by autodetachment.27 Similar effects have also been observed in photoelectron spectra of the deprotonated green fluorescent protein (GFP) chromophore following resonant excitation of higher lying excited 1ππ* states.40,56,57
The 350 nm (3.54 eV) photoelectron spectra are least influenced by resonant excitation of the 21ππ* states, which have VEEs around 3.8–3.9 eV for all three chromophores (Fig. 2), so the 350 nm spectra have been superimposed on the 320 nm and 315 nm photoelectron spectra to highlight the contributions from resonant excitation of the 21ππ* states in the photoelectron spectra recorded at shorter wavelengths. The difference is most pronounced for the 320 nm photoelectron spectrum of pCT−, when the photon energy (3.87 eV) is resonant with the 2ππ* ← S0 absorption maximum (3.86 eV37).
In Fig. 3, the possible electronic relaxation and electron emission processes following photoexcitation of pCA−, pCE− and pCT− are illustrated on a Jablonski diagram and the resonance characters of the first three singlet excited electronic states with respect to the D0 and D1 electronic continua are also shown. The 21ππ* state responsible for the broadening on the high eBE edges of the photoelectron spectra (Fig. 2) is an excited shape resonance with respect to the D0 continuum, implying a strong coupling between the 21ππ* state and the D0 continuum and fast electron emission.
The photoelectron spectra of pCA− and pCE− also have features at high eBE (low eKE) that shift towards higher eBEs as the photon energy increases, characteristic of indirect PD processes. Interestingly, this feature is not observed in the photoelectron spectra of pCT−. The first question to ask is whether the D1 continua are accessible energetically. The VDEs to the first electronically excited states of the neutral radicals have been calculated to be 4.38, 4.32 and 4.50 eV for pCA−, pCE− and pCT−, respectively. The calculated VDEs are higher than the photon energies used in our experiments; however, it is possible that they are overestimated by ∼0.5 eV, similar to the VEEs calculated for the anion using the same method (Table 2), in which case the D1 continuum would be accessible energetically in pCA− and pCE− at 320 and 315 nm and at 350 nm if the ADEs were 0.3–0.4 eV lower than the VDEs. The next question to ask is if any of the excited states of pCA− and pCE− are coupled strongly to the D1 continuum. The 21ππ* state is a Feshbach resonance with respect to the D1 continuum, implying a weak coupling between the 21ππ* state and the D1 continuum and slow electron emission. Thus, indirect electron emission from the 21ππ* state to the D1 continuum is unlikely to compete with fast electron emission to the D0 continuum. However, ultrafast internal conversion to a lower lying excited electronic state may compete with the fast electron emission to the D0 continuum. Subsequent electron emission from lower lying electronically excited states could generate low eKE electrons, if the displacements between the minima of these states and the D0 state were sufficiently small that the Franck–Condon factors were largest between the high vibrational levels of the electronically excited states and high vibrational levels of D0 (Fig. 3). The 11nπ* state has shape resonance character with respect to the D1 continuum and Feshbach character with respect to the D0 continuum, whereas the 11ππ* state has shape resonance character with respect to the D0 detachment continuum and Feshbach character with respect to the D1 continuum. Consequently, it seems likely that any population relaxing to the 11nπ* state would undergo fast electron emission to the D1 continuum or internal conversion to 11ππ* or S0 states and any population relaxing to the 11ππ* state would undergo fast electron emission to the D0 continuum or internal conversion to S0. Internal conversion processes populating high vibrational levels of the electronic ground state of the anion would also result in low eKE electrons from thermionic emission (Fig. 3).
It is possible that the low eKE (high eBE) electrons observed in the 350 nm photoelectron spectra of pCA− and pCE− arise from a different relaxation process than those in the 320 nm and 315 nm photoelectron spectra. The action absorption spectrum for pCE− shows that 350 nm lies between the 11ππ* and 21ππ* absorption bands,22 so the non-zero absorption could be the result of populating high lying vibrational levels of the 11ππ* state which could then autodetach to the D0 continuum or undergo internal conversion to S0 followed by thermionic emission.
Thus, the presence of high eBE (low eKE) electrons in the 315 nm and 320 nm photoelectron spectra of pCA− and pCE−, that shift towards higher eBEs as the photon energy increases, suggests that IC to lower lying electronic states or the ground electronic state compete with PD. The presence of high eBE (low eKE) electrons in the 350 nm photoelectron spectra of pCA− and pCE− can either be attributed to a similar process or to resonant excitation of the 11ππ* state followed by autodetachment or internal conversion to the ground state and thermionic emission. The absence of high eBE (low eKE) electrons in the 315 nm and 320 nm photoelectron spectra of pCT− could be the result of raising the threshold for detachment into the D1 continuum and ‘turning off’ the 11nπ* → D1 + e− detachment channel. However, if this were the case, subsequent electronic relaxation from the 11nπ* state to the 11ππ* or S0 states would themselves generate low eKE electrons by fast electron detachment into the D0 continuum or thermionic emission. A more likely explanation seems to be that lowering the 11nπ* and 11ππ* states in the Franck–Condon region makes conical intersections between the 21ππ* and 11nπ* or 11ππ* states inaccessible or less accessible, effectively ‘turning off’ internal conversion to these states and the ground electronic state, leaving indirect electron emission from the 21ππ* state to the D0 continuum as the only relaxation pathway. The absence of high eBE (low eKE) electrons in the 350 nm photoelectron spectrum of pCT− can be explained either as ‘turning off’ internal conversion processes from the low vibrational levels of the 21ππ* state or ‘turning off’ autodetachment and internal conversion from high vibrational levels of the 11ππ* state. This opens the interesting possibility of using changes to the link between the chromophore and the protein in PYP to manipulate the UV induced electron donor properties of the protein and isolated chromophores, with the potential to monitor and manipulate redox processes.58
Footnote |
† Electronic supplementary information (ESI) available: Synthesis of S-methyl (E)-3-(4-hydroxyphenyl)prop-2-enethioate (pCT). Coordinates of the optimised geometries of the deprotonared anions. See DOI: 10.1039/c6cp00565a |
This journal is © the Owner Societies 2016 |