Excited-state singlet manifold and oscillatory features of a nonatetraeniminium retinal chromophore model

Dual Photoisomerization on Distinct Potential Energy Surfaces in a UV-Absorbing Rhodopsin

UV-absorbing rhodopsins are essential for UV vision and sensing in all kingdoms of life. Unlike the well-known visible-absorbing rhodopsins, which bind a protonated retinal Schiff base for light absorption, UV-absorbing rhodopsins bind an unprotonated retinal Schiff base. Thus far, the photoreaction dynamics and mechanisms of UV-absorbing rhodopsins have remained essentially unknown. Here, we report the complete excited- and ground-state dynamics of the UV form of histidine kinase rhodopsin 1 (HKR1) from eukaryotic algae, using femtosecond stimulated Raman spectroscopy (FSRS) and transient absorption spectroscopy, covering time scales from femtoseconds to milliseconds. We found that energy-level ordering is inverted with respect to visible-absorbing rhodopsins, with an optically forbidden low-lying S₁ excited state that has Ag^– symmetry and a higher-lying UV-absorbing S₂ state of Bu⁺ symmetry. UV-photoexcitation to the S₂ state elicits a unique dual-isomerization reaction: first, C13=C14 cis–trans isomerization occurs during S₂–S₁ evolution in <100 fs. This very fast reaction features the remarkable property that the newly formed isomer appears in the excited state rather than in the ground state. Second, C15=N16 anti–syn isomerization occurs on the S₁–S₀ evolution to the ground state in 4.8 ps. We detected two ground-state unprotonated retinal photoproducts, 13-trans/15-anti (all-trans) and 13-cis/15-syn, after relaxation to the ground state. These isomers become protonated in 58 μs and 3.2 ms, respectively, resulting in formation of the blue-absorbing form of HKR1. Our results constitute a benchmark of UV-induced photochemistry of animal and microbial rhodopsins.

Point Mutation of Anabaena Sensory Rhodopsin Enhances Ground-State Hydrogen Out-of-Plane Wag Raman Activity

The interaction between the retinal protonated Schiff base (RPSB) and surrounding protein residues inside the retinal pocket is believed to play a major role in the ultrafast isomerization of the former. Coherent time-resolved vibrational spectroscopic techniques are applied to reveal the effect of changes in the protein architecture by point mutations (V112N and L83Q) close to the RPSB in Anabaena sensory rhodopsin (ASR). Our study reveals that such point mutations have a minor effect on the low-frequency (<400 cm^-1) torsional modes but dramatically influence the ground-state vibrational Raman activity of the C₁₄-H out-of-plane (HOOP) wag mode (800-820 cm^-1). In mutated ASR, the increase of HOOP Raman activity in the ground state is experimentally observed for the all- trans RPSB, which has shorter excited-state lifetime than in wild-type ASR. This indicates that predistortion of the RPSB inside the mutated retinal pocket is a major factor in the acceleration of the isomerization rate.

Mapping the ultrafast vibrational dynamics of all-trans and 13-cis retinal isomerization in Anabaena Sensory Rhodopsin

Discrepancies in the isomerization dynamics and quantum yields of the trans and cis retinal protonated Schiff base is a well-known issue in the context of retinal photochemistry. Anabaena Sensory Rhodopsin (ASR) is a microbial retinal protein that comprises a retinal chromophore in two ground state (GS) conformations: all-trans, 15-anti (AT) and 13-cis, 15-syn (13C). In this study, we applied impulsive vibrational spectroscopic techniques (DFWM, pump-DFWM and pump-IVS) to ASR to shed more light on how the structural changes take place in the excited state within the same protein environment. Our findings point to distinct features in the ground state structural conformations as well as to drastically different evolutions in the excited state manifold. The ground state vibrational spectra show stronger Raman activity of the C14-H out-of-plane wag (at about 805 cm-1) for the 13C isomer than that for the AT isomer, which hints at a pre-distortion of 13C in the ground state. Evolution of the Raman frequency after interaction with the actinic pulse shows a blue-shift for the C[double bond, length as m-dash]C stretching and CH3 rocking mode for both isomers. For AT, however, the blue-shift is not instantaneous as observed for the 13C isomer, rather it takes more than 200 fs to reach the maximum frequency shift. This frequency blue-shift is rationalized by a decrease in the effective conjugation length during the isomerization reaction, which further confirms a slower formation of the twisted state for the AT isomer and corroborates the presence of a barrier in the excited state trajectory previously predicted by quantum chemical calculations.

Initial excited-state relaxation of locked retinal protonated schiff base chromophore. An insight from coupled cluster and multireference perturbation theory calculations

The initial S₁ excited-state relaxation of retinal protonated Schiff base (RPSB) analog with central C11C12 double bond locked by eight-membered ring (locked-11.8) was investigated by means of multireference perturbation theory methods (XMCQDPT2, XMS-CASPT2, MS-CASPT2) as well as single-reference coupled-cluster CC2 method. The analysis of XMCQDPT2-based geometries reveals rather weak coupling between in-plane and out-of-plane structural evolution and minor energetical relaxation of three locked-11.8 conformers. Therefore, a strong coupling between bonds length inversion and backbone out-of-plane deformation resulting in a very steep S₁ energy profile predicted by CASSCF/CASPT2 calculations is in clear contradiction with the reference XMCQDPT2 results. Even though CC2 method predicts good quality ground-state structures, the excited-state structures display more advanced torsional deformation leading to ca. 0.2 eV exaggerated energy relaxation and significantly red shifted (0.4-0.7 eV) emission maxima. According to our findings, the initial photoisomerization process in locked-11.8, and possibly in other RPSB analogs, studied fully (both geometries and energies) by multireference perturbation theory may be somewhat slower than predicted by CASSCF/CASPT2 or CC2 methods. © 2018 Wiley Periodicals, Inc.

A new twist in the photophysics of the GFP chromophore: a volume-conserving molecular torsion couple

Dynamics of a nonplanar GFP chromophore are studied experimentally and theoretically. Coupled torsional motion is responsible for the ultrafast decay.

The simple structure of the chromophore of the green fluorescent protein (GFP), a phenol and an imidazolone ring linked by a methyne bridge, supports an exceptionally diverse range of excited state phenomena. Here we describe experimentally and theoretically the photochemistry of a novel sterically crowded nonplanar derivative of the GFP chromophore. It undergoes an excited state isomerization reaction accompanied by an exceptionally fast (sub 100 fs) excited state decay. The decay dynamics are essentially independent of solvent polarity and viscosity. Excited state structural dynamics are probed by high level quantum chemical calculations revealing that the fast decay is due to a conical intersection characterized by a twist of the rings and pyramidalization of the methyne bridge carbon. The intersection can be accessed without a barrier from the pre-twisted Franck–Condon structure, and the lack of viscosity dependence is due to the fact that the rings twist in the same direction, giving rise to a volume-conserving decay coordinate. Moreover, the rotation of the phenyl, methyl and imidazolone groups is coupled in the sterically crowded structure, with the methyl group translating the rotation of one ring to the next. As a consequence, the excited state dynamics can be viewed as a torsional couple, where the absorbed photon energy leads to conversion of the out-of-plane orientation from one ring to the other in a volume conserving fashion. A similar modification of the range of methyne dyes may provide a new family of devices for molecular machines, specifically torsional couples.

Excited-state minima and emission energies of retinal chromophore analogues: Performance of CASSCF and CC2 methods as compared with CASPT2

This study provides gas-phase S₁ excited-state geometries along with emission and adiabatic energies for methylated/demethylated and ring-locked analogues of protonated Schiff base retinal models comprising system of five conjugated double bonds (PSB5), using second order multiconfiguration perturbation theory (CASPT2). CASPT2 results serve as reference data to assess the performance of CC2 (second-order approximate coupled cluster singles and doubles) and a commonly used CASSCF/CASPT2 protocol, that is, complete active space self-consistent field (CASSCF) geometry optimization followed by CASPT2 energy calculation. We find that the CASSCF methodology fails to locate planar S₁ minimum energy structures for four out of five investigated planar models in contrast to CC2 and CASPT2 methods. However, for those which were found: one planar and two twisted minima, there is an excellent agreement between CASSCF and CASPT2 results in terms of geometrical parameters, one-electron properties, as well as emission and adiabatic energies. CC2 performs well for in-plane S₁ minima and their spectroscopic and electronic properties. However, this picture deteriorates for twisted minima. As expected, the CC2 description of the S₂ electronic state, with strong multireference and significant double excitation character, is very poor, exhibiting errors in transition energies exceeding 1 eV. They may be substantially diminished by recalculating transition energies with CASPT2 method. Our work shows that CASSCF/CASPT2 and CC2 shortcomings may influence gas-phase retinal analogues' excited state description in a dramatic way. © 2017 Wiley Periodicals, Inc.

Critical role of Asp227 in the photocycle of proteorhodopsin

The photocycle of the proton acceptor complex mutant D227N of the bacterial retinal protein proteorhodopsin is investigated employing steady state pH-titration experiments in the UV-visible range as well as femtosecond-pump-probe spectroscopy and flash photolysis in the visible spectral range. The evaluation of the pH-dependent spectra showed that the neutralization of the charge at position 227 has a remarkable influence on the ground state properties of the protein. Both the pK(a) values of the primary proton acceptor and of the Schiff base are considerably decreased. Femtosecond-time-resolved measurements demonstrate that the general S(1) deactivation pathway; that is, the K-state formation is preserved in the D227N mutant. However, the pH-dependence of the reaction rate is lost by the substitution of Asp227 with an asparagine. Also no significant kinetic differences are observed upon deuteration. This is explained by the lack of a strongly hydrogen-bonded water in the vicinity of Asp97, Asp227, and the Schiff base or a change in the hydrogen bonding of it (Ikeda et al. (2007) Biochemistry 46, 5365-5373). The flash photolysis measurements prove a considerably elongated photocycle with pronounced pH-dependence. Interestingly, at pH 9 the M-state is visible until the end of the reaction cycle, leading to the conclusion that the mutation does not only lower the pK(a) of the Schiff base in the unphotolyzed ground state but also prevents an efficient reprotonation reaction.

Electrostatic control of the photoisomerization efficiency and optical properties in visual pigments: on the role of counterion quenching

Hybrid QM(CASPT2//CASSCF/6-31G*)/MM(Amber) computations have been used to map the photoisomerization path of the retinal chromophore in Rhodopsin and explore the reasons behind the photoactivity efficiency and spectral control in the visual pigments. It is shown that while the electrostatic environment plays a central role in properly tuning the optical properties of the chromophore, it is also critical in biasing the ultrafast photochemical event: it controls the slope of the photoisomerization channel as well as the accessibility of the S(1)/S(0) crossing space triggering the ultrafast decay. The roles of the E113 counterion, the E181 residue, and the other amino acids of the protein pocket are explicitly analyzed: it appears that counterion quenching by the protein environment plays a key role in setting up the chromophore's optical properties and its photochemical efficiency. A unified scenario is presented that discloses the relationship between spectroscopic and mechanistic properties in rhodopsins and allows us to draw a solid mechanism for spectral tuning in color vision pigments: a tunable counterion shielding appears as the elective mechanism for L<-->M spectral modulation, while a retinal conformational control must dictate S absorption. Finally, it is suggested that this model may contribute to shed new light into mutations-related vision deficiencies that opens innovative perspectives for experimental biomolecular investigations in this field.

Functional electric field changes in photoactivated proteins revealed by ultrafast Stark spectroscopy of the Trp residues

Ultrafast transient absorption spectroscopy of wild-type bacteriorhodopsin (WT bR) and 2 tryptophan mutants (W86F and W182F) is performed with visible light excitation (pump) and UV probe. The aim is to investigate the photoinduced change in the charge distribution with 50-fs time resolution by probing the effects on the tryptophan absorption bands. A systematic, quantitative comparison of the transient absorption of the 3 samples is carried out. The main result is the absence in the W86F mutant of a transient induced absorption band observed at approximately 300-310 nm in WT bR and W182F. A simple model describing the dipolar interaction of the retinal moiety with the 2 tryptophan residues of interest allows us to reproduce the dominant features of the transient signals observed in the 3 samples at ultrashort pump-probe delays. In particular, we show that Trp(86) undergoes a significant Stark shift induced by the transient retinal dipole moment. The corresponding transient signal can be isolated by direct subtraction of experimental data obtained for WT bR and W86F. It shows an instantaneous rise, followed by a decay over approximately 500 fs corresponding to the isomerization time. Interestingly, it does not decay back to zero, thus revealing a change in the local electrostatic environment that remains long after isomerization, in the K intermediate state of the protein cycle. The comparison of WT bR and W86F also leads to a revised interpretation of the overall transient UV absorption of bR.

Color tuning in rhodopsins: the mechanism for the spectral shift between bacteriorhodopsin and sensory rhodopsin II

The mechanism of color tuning in the rhodopsin family of proteins has been studied by comparing the optical properties of the light-driven proton pump bacteriorhodopsin (bR) and the light detector sensory rhodopsin II (sRII). Despite a high structural similarity, the maximal absorption is blue-shifted from 568 nm in bR to 497 nm in sRII. The molecular mechanism of this shift is still a matter of debate, and its clarification sheds light onto the general mechanisms of color tuning in retinal proteins. The calculations employ a combined quantum mechanical/molecular mechanical (QM/MM) technique, using a DFT-based method for ground state properties and the semiempirical OM2/MRCI method and ab initio SORCI method for excited state calculations. The high efficiency of the methodology has allowed us to study a wide variety of aspects including dynamical effects. The absorption shift as well as various mutation experiments and vibrational properties have been successfully reproduced. Our results indicate that several sources contribute to the spectral shift between bR and sRII. The main factors are the counterion region at the extracellular side of retinal and the amino acid composition of the binding pocket. Our analysis allows a distinction and identification of the different effects in detail and leads to a clear picture of the mechanism of color tuning, which is in good agreement with available experimental data.

Photochemistry of a Retinal Protonated Schiff-Base Analogue Mimicking the Opsin Shift of Bacteriorhodopsin

A retinal Schiff base analogue which artificially mimics the protein-induced red shifting of absorption in bacteriorhodopsin (BR) has been investigated with femtosecond multichannel pump probe spectroscopy. The objective is to determine if the catalysis of retinal internal conversion in the native protein BR, which absorbs at 570 nm, is directly correlated with the protein-induced Stokes shifting of this absorption band otherwise known as the "opsin shift". Results demonstrate that the red shift afforded in the model system does not hasten internal conversion relative to that taking place in a free retinal-protonated Schiff base (RPSB) in methanol solution, and stimulated emission takes place with biexponential kinetics and characteristic timescales of approximately 2 and 10.5 ps. This shows that interactions between the prosthetic group and the protein that lead to the opsin shift in BR are not directly involved in reducing the excited-state lifetime by nearly an order of magnitude. A sub-picosecond phase of spectral evolution, analogues of which are detected in photoexcited retinal proteins and RPSBs in solution, is observed after excitation anywhere within the intense visible absorption band. It consists of a large and discontinuous spectral shift in excited-state absorption and is assigned to electronic relaxation between excited states, a scenario which might also be relevant to those systems as well. Finally, a transient excess bleach component that tunes with the excitation wavelength is detected in the data and tentatively assigned to inhomogeneous broadening in the ground state absorption band. Possible sources of such inhomogeneity and its relevance to native RPSB photochemistry are discussed.

Following Photoinduced Dynamics in Bacteriorhodopsin with 7-fs Impulsive Vibrational Spectroscopy

Sub-10-fs laser pulses are used to impulsively photoexcite bacteriorhodopsin (BR) suspensions and probe the evolution of the resulting vibrational wave packets. Fourier analysis of the spectral modulations induced by transform-limited as well as linearly chirped excitation pulses allows the delineation of excited- and ground-state contributions to the data. On the basis of amplitude and phase variations of the modulations as a function of the dispersed probe wavelength, periodic modulations in absorption above 540 nm are assigned to ground-state vibrational coherences induced by resonance impulsive Raman spectral activity (RISRS). Probing at wavelengths below 540 nm-the red edge of the intense excited-state absorption band-uncovers new vibrational features which are accordingly assigned to wave packet motions along bound coordinates on the short-lived reactive electronic surface. They consist of high- and low-frequency shoulders adjacent to the strong C=C stretching and methyl rock modes, respectively, which have ground-state frequencies of 1008 and 1530 cm-1. Brief activity centered at approximately 900 cm-1, which is characteristic of ground-state HOOP modes, and strong modulations in the torsional frequency range appear as well. Possible assignments of the bands and their implication to photoinduced reaction dynamics in BR are discussed. Reasons for the absence of similar signatures in the pump-probe spectral modulations at longer probing wavelengths are considered as well.

Computational photochemistry of retinal proteins

High spectral tunability and quantum yield are the striking features of rhodopsin photochemistry. They rely on a strong and complex interaction of their chromophore, the protonated Schiff base of retinal, with its protein environment. In this article, we review the progress in the computational modeling of these systems, focusing on the optical properties and the excited state dynamics. While the earlier success of atomistic theoretical models was based on the breakthrough in X-ray crystallography and combined quantum mechanical molecular mechanical (QM/MM) methodology, recent advances point out the importance of high-level QM methods and the incorporation of effects that are neglected in conventional QM/MM or ONIOM schemes, like polarization and charge transfer.

Ultrafast excited state dynamics of the protonated Schiff base of all-trans retinal in solvents

We present a comparative study of the ultrafast photophysics of all-trans retinal in the protonated Schiff base form in solvents with different polarities and viscosities. Steady-state spectra of retinal in the protonated Schiff base form show large absorption-emission Stokes shifts (6500-8100 cm(-1)) for both polar and nonpolar solvents. Using a broadband fluorescence up-conversion experiment, the relaxation kinetics of fluorescence is investigated with 120 fs time resolution. The time-zero spectra already exhibit a Stokes-shift of approximately 6000 cm(-1), indicating depopulation of the Franck-Condon region in < or =100 fs. We attribute it to relaxation along skeletal stretching. A dramatic spectral narrowing is observed on a 150 fs timescale, which we assign to relaxation from the S(2) to the S(1) state. Along with the direct excitation of S(1), this relaxation populates different quasistationary states in S(1), as suggested from the existence of three distinct fluorescence decay times with different decay associated spectra. A 0.5-0.65 ps decay component is observed, which may reflect the direct repopulation of the ground state, in line with the small isomerization yield in solvents. Two longer decay components are observed and are attributed to torsional motion leading to photo-isomerization. The various decay channels show little or no dependence with respect to the viscosity or dielectric constant of the solvents. This suggests that in the protein, the bond selectivity of isomerization is mainly governed by steric effects.

Structure, Spectroscopy, and Spectral Tuning of the Gas-Phase Retinal Chromophore: The β-Ionone “Handle” and Alkyl Group Effect

The low-lying singlet states (i.e. S0, S1, and S2) of the chromophore of rhodopsin, the protonated Schiff base of 11-cis-retinal (PSB11), and of its all-trans photoproduct have been studied in isolated conditions by using ab initio multiconfigurational second-order perturbation theory. The computed spectroscopic features include the vertical excitation, the band origin, and the fluorescence maximum of both isomers. On the basis of the S0-->S1 vertical excitation, the gas-phase absorption maximum of PSB11 is predicted to be 545 nm (2.28 eV). Thus, the predicted absorption maximum appears to be closer to that of the rhodopsin pigment (2.48 eV) and considerably red-shifted with respect to that measured in solution (2.82 eV in methanol). In addition, the absorption maxima associated with the blue, green, and red cone visual pigments are tentatively rationalized in terms of the spectral changes computed for PSB11 structures featuring differently twisted beta-ionone rings. More specifically, a blue-shifted absorption maximum is explained in terms of a large twisting of the beta-ionone ring (with respect to the main conjugated chain) in the visual S-cone (blue) pigment chromophore. In contrast, the chromophore of the visual L-cone (red) pigment is expected to have a nearly coplanar beta-ionone ring yielding a six double bond fully conjugated framework. Finally, the M-cone (green) chromophore is expected to feature a twisting angle between 10 and 60 degrees. The spectroscopic effects of the alkyl substituents on the PSB11 spectroscopic properties have also been investigated. It is found that they have a not negligible stabilizing effect on the S1-S0 energy gap (and, thus, cause a red shift of the absorption maximum) only when the double bond of the beta-ionone ring conjugates significantly with the rest of the conjugated chain.

Contrasting the excited-state dynamics of the photoactive yellow protein chromophore: protein versus solvent environments

Wavelength- and time-resolved fluorescence experiments have been performed on the photoactive yellow protein, the E46Q mutant, the hybrids of these proteins containing a nonisomerizing "locked" chromophore, and the native and locked chromophores in aqueous solution. The ultrafast dynamics of these six systems is compared and spectral signatures of isomerization and solvation are discussed. We find that the ultrafast red-shifting of fluorescence is associated mostly with solvation dynamics, whereas isomerization manifests itself as quenching of fluorescence. The observed multiexponential quenching of the protein samples differs from the single-exponential lifetimes of the chromophores in solution. The locked chromophore in the protein environment decays faster than in solution. This is due to additional channels of excited-state energy dissipation via the covalent and hydrogen bonds with the protein environment. The observed large dispersion of quenching timescales observed in the protein samples that contain the native pigment favors both an inhomogeneous model and an excited-state barrier for isomerization.

Calculating Absorption Shifts for Retinal Proteins: Computational Challenges

Rhodopsins can modulate the optical properties of their chromophores over a wide range of wavelengths. The mechanism for this spectral tuning is based on the response of the retinal chromophore to external stress and the interaction with the charged, polar, and polarizable amino acids of the protein environment and is connected to its large change in dipole moment upon excitation, its large electronic polarizability, and its structural flexibility. In this work, we investigate the accuracy of computational approaches for modeling changes in absorption energies with respect to changes in geometry and applied external electric fields. We illustrate the high sensitivity of absorption energies on the ground-state structure of retinal, which varies significantly with the computational method used for geometry optimization. The response to external fields, in particular to point charges which model the protein environment in combined quantum mechanical/molecular mechanical (QM/MM) applications, is a crucial feature, which is not properly represented by previously used methods, such as time-dependent density functional theory (TDDFT), complete active space self-consistent field (CASSCF), and Hartree-Fock (HF) or semiempirical configuration interaction singles (CIS). This is discussed in detail for bacteriorhodopsin (bR), a protein which blue-shifts retinal gas-phase excitation energy by about 0.5 eV. As a result of this study, we propose a procedure which combines structure optimization or molecular dynamics simulation using DFT methods with a semiempirical or ab initio multireference configuration interaction treatment of the excitation energies. Using a conventional QM/MM point charge representation of the protein environment, we obtain an absorption energy for bR of 2.34 eV. This result is already close to the experimental value of 2.18 eV, even without considering the effects of protein polarization, differential dispersion, and conformational sampling.

Counterion Controlled Photoisomerization of Retinal Chromophore Models: a Computational Investigation

CASPT2//CASSCF photoisomerization path computations have been used to unveil the effects of an acetate counterion on the photochemistry of two retinal protonated Schiff base (PSB) models: the 2-cis-penta-2,4-dieniminium and the all-trans-epta-2,4,6-trieniminium cations. Different positions/orientations of the counterion have been investigated and related to (i) the spectral tuning and relative stability of the S0, S1, and S2 singlet states; (ii) the selection of the photochemically relevant excited state; (iii) the control of the radiationless decay and photoisomerization rates; and, finally, (iv) the control of the photoisomerization stereospecificity. A rationale for the results is given on the basis of a simple (electrostatic) qualitative model. We show that the model readily explains the computational results providing a qualitative explanation for different aspects of the experimentally observed "environment" dependent PSB photochemistry. Electrostatic effects likely involved in controlling retinal photoisomerization stereoselectivity in the protein are also discussed under the light of these results, and clues for a stereocontrolled electrostatically driven photochemical process are presented. These computations provide a rational basis for the formulation of a mechanistic model for photoisomerization electrostatic catalysis.

A global investigation of excited state surfaces within time-dependent density-functional response theory

This work investigates the capability of time-dependent density functional response theory to describe excited state potential energy surfaces of conjugated organic molecules. Applications to linear polyenes, aromatic systems, and the protonated Schiff base of retinal demonstrate the scope of currently used exchange-correlation functionals as local, adiabatic approximations to time-dependent Kohn-Sham theory. The results are compared to experimental and ab initio data of various kinds to attain a critical analysis of common problems concerning charge transfer and long range (nondynamic) correlation effects. This analysis goes beyond a local investigation of electronic properties and incorporates a global view of the excited state potential energy surfaces.