About the intrinsic photochemical properties of the 11-cis retinal chromophore: computational clues for a trap state and a lever effect in Rhodopsin catalysis

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Trajectory surface hopping molecular dynamics simulations for retinal protonated Schiff-base photoisomerization

Global switching trajectory surface hopping molecular dynamics simulations are performed using accurate on-the-fly (TD)CAM-B3LYP/6-31G potential energy surfaces to study retinal protonated Schiff-base photoisomerization up to S₁ excitation. The simulations detected two-layer conical intersection networks: one is at an energy as high as 8 eV and the other is in the energy range around 3-4 eV. Six conical intersections within the low-layer energy region that correspond to active conical intersections under experimental conditions are found via the use of pairwise isomers, within which nonadiabatic molecular dynamics simulations are performed. Eight isomer products are populated with simulated sampling trajectories from which the simulated quantum yield in the gas phase is estimated to be 0.11 (0.08) moving from the all-trans isomer to the 11-cis (11-cis to all-trans) isomer in comparison with an experimental value of 0.09 (0.2) in the solution phase. Each conical intersection is related to one specific twist angle accompanying a related CC double bond motion during photoisomerization. Nonplanar distortion of the entire dynamic process has a significant role in the formation of the relevant photoisomerization products. The present simulation indicates that all hopping points show well-behaved potential energy surface topology, as calculated via the conventional TDDFT method, at conical intersections between S₁ and S₀ states. Therefore, the present nonadiabatic dynamics simulations with the TDDFT method are very encouraging for simulating various large systems related to retinal Schiff-base photoisomerization in the real world.

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.

Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems

A new flavor of the frozen natural orbital complete active space second-order perturbation theory method (FNO-CASPT2, Aquilante et al., J. Chem. Phys. 131, 034113) is proposed herein. In this new implementation, the virtual space in Cholesky decomposition-based CASPT2 computations (CD-CASPT2) is truncated by excluding those orbitals that contribute the least toward preserving a predefined value of the trace of an approximate density matrix, as that represents a measure of the amount of dynamic correlation retained in the model. In this way, the amount of correlation included is practically constant at all nuclear arrangements, thus allowing for the computation of smooth electronic states surfaces and energy gradients-essential requirements for theoretical studies in photochemistry. The method has been benchmarked for a series of relevant biochromophores for which large speed-ups have been recorded while retaining the accuracy achieved in the corresponding CD-CASPT2 calculations. Both vertical excitation energies and gradient calculations have been carried out to establish general guidelines as to how much correlation needs to be retained in the calculation for the results to be consistent with the CD-CASPT2 findings. Our results feature errors within a tenth of an eV for the most difficult cases and have been validated to be used for gradient computations where an up to 3-fold speed-up is observed depending on the size of the system and the basis set employed.

Molecular bases for the selection of the chromophore of animal rhodopsins

The functions of microbial and animal rhodopsins are triggered by the isomerization of their all-trans and 11-cis retinal chromophores, respectively. To lay the molecular basis driving the evolutionary transition from the all-trans to the 11-cis chromophore, multiconfigurational quantum chemistry is used to compare the isomerization mechanisms of the sensory rhodopsin from the cyanobacterium Anabaena PCC 7120 (ASR) and of the bovine rhodopsin (Rh). It is found that, despite their evolutionary distance, these eubacterial and vertebrate rhodopsins start to isomerize via distinct implementations of the same bicycle-pedal mechanism originally proposed by Warshel [Warshel A (1976) Nature 260:678-683]. However, by following the electronic structure changes of ASR (featuring the all-trans chromophore) during the isomerization, we find that ASR enters a region of degeneracy between the first and second excited states not found in Rh (featuring the 11-cis chromophore). We show that such degeneracy is modulated by the preorganized structure of the chromophore and by the position of the reactive double bond. It is argued that the optimization of the electronic properties of the chromophore, which affects the photoisomerization efficiency and the thermal isomerization barrier, provided a key factor for the emergence of the striking amino acid sequence divergence observed between the microbial and animal rhodopsins.

Modelling retinal chromophores photoisomerization: from minimal models in vacuo to ultimate bidimensional spectroscopy in rhodopsins

Modelling of retinal photoisomerization in different environments is reviewed and ultimate ultrafast electronic spectroscopy is proposed for obtaining new insights.

Gas-Phase Retinal Spectroscopy: Temperature Effects Are But a Mirage

We employ state-of-the-art first-principle approaches to investigate whether temperature effects are responsible for the unusually broad and flat spectrum of protonated Schiff base retinal observed in photodissociation spectroscopy, as has recently been proposed. We first carefully calibrate how to construct a realistic geometrical model of retinal and show that the exchange-correlation M06-2X functional yields an accurate description while the commonly used complete active space self-consistent field method (CASSCF) is not adequate. Using modern multiconfigurational perturbative methods (NEVPT2) to compute the excitations, we then demonstrate that conformations with different orientations of the β-ionone ring are characterized by similar excitations. Moreover, other degrees of freedom identified as active in room-temperature molecular dynamics simulations do not yield the shift required to explain the anomalous spectral shape. Our findings indicate that photodissociation experiments are not representative of the optical spectrum of retinal in the gas phase and call for further experimental characterization of the dissociation spectra.

Aborted double bicycle-pedal isomerization with hydrogen bond breaking is the primary event of bacteriorhodopsin proton pumping

Quantum mechanics/molecular mechanics calculations based on ab initio multiconfigurational second order perturbation theory are employed to construct a computer model of Bacteriorhodopsin that reproduces the observed static and transient electronic spectra, the dipole moment changes, and the energy stored in the photocycle intermediate K. The computed reaction coordinate indicates that the isomerization of the retinal chromophore occurs via a complex motion accounting for three distinct regimes: (i) production of the excited state intermediate I, (ii) evolution of I toward a conical intersection between the excited state and the ground state, and (iii) formation of K. We show that, during stage ii, a space-saving mechanism dominated by an asynchronous double bicycle-pedal deformation of the C10═C11─C12═C13─C14═N moiety of the chromophore dominates the isomerization. On this same stage a N─H/water hydrogen bond is weakened and initiates a breaking process that is completed during stage iii.

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.

Acceleration of the Z to E photoisomerization of penta-2,4-dieniminium by hydrogen out-of-plane motion: theoretical study on a model system of retinal protonated Schiff base

We report the result of comparison between two reaction coordinates [on the potential energy surface of the first excited state (S(1))] produced by CASSCF and these energies recalculated by MRMP2 in the Z to E photoisomerization of penta-2,4-dieniminium (PDI) as the minimal model of the retinal protonated Schiff base (RPSB). One coordinate is the S(1) state minimum-energy-path (MEP) in mass-weighted coordinates from the S(1) vertically excited point, where a strong hydrogen-out-of plane (HOOP) motion is not exhibited. The energy profile of the S(1) MEP at the MRMP2//CASSCF level shows a barrier for the rotation around the reactive C-C and hits the S(1)/S(0) degeneracy space where the central C-C-C-C dihedral angle is distorted by 65 degrees . The other coordinate is an S(1) coordinate obtained by the relaxed scan strategy. The relaxed coordinate along the central C-C-C-C dihedral angle, which we call the HOOP coordinate, shows strong HOOP motion. According to the MRMP2//CASSCF calculation, there is no barrier on the HOOP coordinate. Furthermore, the S(1) to S(0) transition may be possible without the large skeletal deformation by HOOP motion because the HOOP coordinate encounters the S(1)/S(0) degeneracy space where the central C-C-C-C dihedral angle is distorted by only 40 degrees . Consequently, if PDI is a suitable model molecule for the RPSB as often assumed, the 11-cis to all-trans photoisomerization is predicted to be accelerated by the HOOP motion.