Is the Photoinduced Isomerization in Retinal Protonated Schiff Bases a Single- or Double-Torsional Process?

Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment

Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.

Recommendations for Velocity Adjustment in Surface Hopping

This study investigates velocity adjustment directions after hopping in surface hopping dynamics. Using fulvene and a protonated Schiff base (PSB4) as case studies, we investigate the population decay and reaction yields of different sets of dynamics with the velocity adjusted in either the nonadiabatic coupling, gradient difference, or momentum directions. For the latter, in addition to the conventional algorithm, we investigated the performance of a reduced kinetic energy reservoir approach recently proposed. Our evaluation also considered velocity adjustment in the directions of approximate nonadiabatic coupling vectors. While results for fulvene are susceptible to the adjustment approach, PSB4 is not. We correlated this dependence to the topography near the conical intersections. When nonadiabatic coupling vectors are unavailable, the gradient difference direction is the best adjustment option. If the gradient difference is also unavailable, a semiempirical vector direction or the momentum direction with a reduced kinetic energy reservoir becomes an excellent option to prevent an artificial excess of back hoppings. The precise velocity adjustment direction is less crucial for describing the nonadiabatic dynamics than the kinetic energy reservoir's size.

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.

Benchmarking two-photon absorption strengths of rhodopsin chromophore models with CC3 and CCSD methodologies: An assessment of popular density functional approximations

This work presents the investigations of the impact of an increasing electron correlation in the hierarchy of coupled-cluster methods, i.e., CC2, CCSD, and CC3, on two-photon absorption (2PA) strengths for the lowest excited state of the minimal rhodopsin's chromophore model-cis-penta-2,4-dieniminium cation (PSB3). For a larger chromophore's model [4-cis-hepta-2,4,6-trieniminium cation (PSB4)], CC2 and CCSD calculations of 2PA strengths were performed. Additionally, 2PA strengths predicted by some popular density functional theory (DFT) functionals differing in HF exchange contribution were assessed against the reference CC3/CCSD data. For PSB3, the accuracy of 2PA strengths increases in the following order: CC2 < CCSD < CC3, with the CC2 deviation from both higher-level methods exceeding 10% at 6-31+G* basis sets and 2% at aug-cc-pVDZ basis set. However, for PSB4, this trend is reversed and CC2-based 2PA strength is larger than the corresponding CCSD value. Among the DFT functionals investigated, CAM-B3LYP and BHandHLYP provide 2PA strengths in best compliance with reference data, however, with the error approaching an order of magnitude.

A study of the photochemical behaviour and relaxation mechanisms of anti–syn isomerisation around quinazolinone –N–N bonds

High-resolution NMR experiments revealed that differently substituted quinazolinone-based Schiff bases undergo anti to syn isomerisation on exposure to ultraviolet light in DMSO solution.

How a linear triazene photoisomerizes in a volume-conserving fashion

Understanding deactivation mechanisms of functional groups is a key step to design novel photo-active devices and molecular imaging agents. Here, we elucidate the photochemistry of linear triazenes, an extended analogue of the photo-switchable azo group, exemplarily for the widely used DNA-minor-groove binder berenil. Combining ultrafast spectroscopy and ab initio calculations unveils that the E-azo,s-trans structure of berenil predominates in the gas phase and in aqueous solution, and ADC(2) intrinsic reaction coordinate calculations disclose that the excited-state relaxation to the S1 minima/conical intersections follows a two-step mechanism: N[double bond, length as m-dash]N bond stretching followed by a bicycle-pedal rotation in the triazene bridge. Furthermore, studying the ground-state pathways shows that a fraction of the molecules relaxes back to the E-azo,s-trans isomer while the other part photoisomerizes to the Z-azo,s-trans via a hula-twist motion, as evidenced by experimental quantum yields of Φ ≈ 0.5 found for berenil in water, ethylene glycol, or bound to β-trypsin. Moreover, our studies show that while the excited-state relaxation is insensitive to the environment, the ground-state dynamics depend on biomolecular binding partners.

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.

Ultrafast photochemistry of anabaena sensory rhodopsin: experiment and theory

Light induced isomerization of the retinal chromophore activates biological function in all retinal protein (RP) driving processes such as ion-pumping, vertebrate vision and phototaxis in organisms as primitive as archea, or as complex as mammals. This process and its consecutive reactions have been the focus of experimental and theoretical research for decades. The aim of this review is to demonstrate how the experimental and theoretical research efforts can now be combined to reach a more comprehensive understanding of the excited state process on the molecular level. Using the Anabaena Sensory Rhodopsin as an example we will show how contemporary time-resolved spectroscopy and recently implemented excited state QM/MM methods consistently describe photochemistry in retinal proteins. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

Nonadiabatic photodynamics of a retinal model in polar and nonpolar environment

The nonadiabatic photodynamics of the all-trans-2,4-pentadiene-iminium cation (protonated Schiff base 3, PSB3) and the all-trans-3-methyl-2,4-pentadiene-iminium cation (MePSB3) were investigated in the gas phase and in polar (aqueous) and nonpolar (n-hexane) solutions by means of surface hopping using a multireference configuration-interaction (MRCI) quantum mechanical/molecular mechanics (QM/MM) level. Spectra, lifetimes for radiationless deactivation to the ground state, and structural and electronic parameters are compared. A strong influence of the polar solvent on the location of the crossing seam, in particular in the bond length alternation (BLA) coordinate, is found. Additionally, inclusion of the polar solvent changes the orientation of the intersection cone from sloped in the gas phase to peaked, thus enhancing considerably its efficiency for deactivation of the molecular system to the ground state. These factors cause, especially for MePSB3, a substantial decrease in the lifetime of the excited state despite the steric inhibition by the solvent.

Control Mechanisms of Photoisomerization in Protonated Schiff Bases

We performed ab initio excited-state molecular dynamics simulations of a gas-phase photoexcited protonated Schiff base (C1-N2═C3-C4═C5-C6) to search for control mechanisms of its photoisomerization. The excited molecule twists by ∼90° around either the N2C3 bond or the C4C5 bond and relaxes to the ground electronic state through a conical intersection with either a trans or cis outcome. We show that a large initial distortion of several dihedral angles and a specific normal vibrational mode combining pyramidalization and double-bond twisting can lead to a preferential rotation of atoms around the C4C5 bond. We also show that selective pretwisting of several dihedral angles in the initial ground state thermal ensemble (by analogy to a protein pocket) can significantly increase the fraction of photoreactive (cis → trans) trajectories. We demonstrate that new ensembles with higher degrees of control over the photoisomerization reaction can be obtained by a computational directed evolution approach on the ensembles of molecules with the pretwisted geometries.

The effect of protein environment on photoexcitation properties of retinal

Retinal is the photon absorbing chromophore of rhodopsin and other visual pigments, enabling the vertebrate vision process. The effects of the protein environment on the primary photoexcitation process of retinal were studied by time-dependent density functional theory (TDDFT) and the algebraic diagrammatic construction through second order (ADC(2)) combined with our recently introduced reduction of virtual space (RVS) approximation method. The calculations were performed on large full quantum chemical cluster models of the bluecone (BC) and rhodopsin (Rh) pigments with 165-171 atoms. Absorption wavelengths of 441 and 491 nm were obtained at the B3LYP level of theory for the respective models, which agree well with the experimental values of 414 and 498 nm. Electrostatic rather than structural strain effects were shown to dominate the spectral tuning properties of the surrounding protein. The Schiff base retinal and a neighboring Glu-113 residue were found to have comparable proton affinities in the ground state of the BC model, whereas in the excited state, the proton affinity of the Schiff base is 5.9 kcal/mol (0.26 eV) higher. For the ground and excited states of the Rh model, the proton affinity of the Schiff base is 3.2 kcal/mol (0.14 eV) and 7.9 kcal/mol (0.34 eV) higher than for Glu-113, respectively. The protein environment was found to enhance the bond length alternation (BLA) of the retinyl chain and blueshift the first absorption maxima of the protonated Schiff base in the BC and Rh models relative to the chromophore in the gas phase. The protein environment was also found to decrease the intensity of the second excited state, thus improving the quantum yield of the photoexcitation process. Relaxation of the BC model on the excited state potential energy surface led to a vanishing BLA around the isomerization center of the conjugated retinyl chain, rendering the retinal accessible for cis-trans isomerization. The energy of the relaxed excited state was found to be 30 kcal/mol (1.3 eV) above the minimum ground state energy, and might be related to the transition state of the thermal activation process.

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.

Studies of the ground and excited-state surfaces of the retinal chromophore using CAM-B3LYP

The isomerization of the 11-cis isomer (PSB11) of the retinal chromophore to its all-trans isomer (PSBT) is examined. Optimized structures on both the ground state and the excited state are calculated, and the dependence on torsional angles in the carbon chain is investigated. Time-dependent density functional theory is used to produce excitation energies and the excited-state surface. To avoid problems with the description of excited states that can arise with standard DFT methods, the CAM-B3LYP functional was used. Comparing CAM-B3LYP with B3LYP results indicates that the former is significantly more accurate, as a consequence of which detailed cross sections of the retinal excited-state surface are obtained.