The retinal chromophore/chloride ion pair: structure of the photoisomerization path and interplay of charge transfer and covalent states

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.

Comparative quantum-classical dynamics of natural and synthetic molecular rotors show how vibrational synchronization modulates the photoisomerization quantum efficiency

We use quantum-classical trajectories to investigate the origin of the different photoisomerization quantum efficiency observed in the dim-light visual pigment Rhodopsin and in the light-driven biomimetic molecular rotor para-methoxy N-methyl indanylidene-pyrrolinium (MeO-NAIP) in methanol. Our results reveal that effective light-energy conversion requires, in general, an auxiliary molecular vibration (called promoter) that does not correspond to the rotary motion but synchronizes with it at specific times. They also reveal that Nature has designed Rhodopsin to exploit two mechanisms working in a vibrationally coherent regime. The first uses a wag promoter to ensure that ca. 75% of the absorbed photons lead to unidirectional rotations. The second mechanism ensures that the same process is fast enough to avoid directional randomization. It is found that MeO-NAIP in methanol is incapable of exploiting the above mechanisms resulting into a 50% quantum efficiency loss. However, when the solvent is removed, MeO-NAIP rotation is predicted to synchronize with a ring-inversion promoter leading to a 30% increase in quantum efficiency and, therefore, biomimetic behavior.

Retinal photoisomerization versus counterion protonation in light and dark-adapted bacteriorhodopsin and its primary photoproduct

Discovered over 50 years ago, bacteriorhodopsin is the first recognized and most widely studied microbial retinal protein. Serving as a light-activated proton pump, it represents the archetypal ion-pumping system. Here we compare the photochemical dynamics of bacteriorhodopsin light and dark-adapted forms with that of the first metastable photocycle intermediate known as "K". We observe that following thermal double isomerization of retinal in the dark from bio-active all-trans 15-anti to 13-cis, 15-syn, photochemistry proceeds even faster than the ~0.5 ps decay of the former, exhibiting ballistic wave packet curve crossing to the ground state. In contrast, photoexcitation of K containing a 13-cis, 15-anti chromophore leads to markedly multi-exponential excited state decay including much slower stages. QM/MM calculations, aimed to interpret these results, highlight the crucial role of protonation, showing that the classic quadrupole counterion model poorly reproduces spectral data and dynamics. Single protonation of ASP212 rectifies discrepancies and predicts triple ground state structural heterogeneity aligning with experimental observations. These findings prompt a reevaluation of counter ion protonation in bacteriorhodopsin and contribute to the broader understanding of its photochemical dynamics.

Photochemistry of the Light-Driven Sodium Pump Krokinobacter eikastus Rhodopsin 2 and Its Implications on Microbial Rhodopsin Research: Retrospective and Perspective

The discovery of the light-driven sodium pump Krokinobacter eikastus rhodopsin 2 (KR2) in 2013 has changed the paradigm that cation transport in microbial rhodopsins is restricted to the translocation of protons. Even though this finding is already remarkable by itself, it also reignited more general discussions about the functional mechanism of ion transport. The unique composition of the retinal binding pocket in KR2 with a tight interaction between the retinal Schiff base and its respective counterion D116 also has interesting implications on the photochemical pathway of the chromophore. Here, we discuss the most recent advances in our understanding of the KR2 functionality from the primary event of photon absorption by all-trans retinal up to the actual protein response in the later phases of the photocycle, mainly from the point of view of optical spectroscopy. In this context, we furthermore highlight some of the ongoing debates on the photochemistry of microbial rhodopsins and give some perspectives for promising future directions in this field of research.

Control of Protonated Schiff Base Excited State Decay within Visual Protein Mimics: A Unified Model for Retinal Chromophores

Artificial biomimetic chromophore-protein complexes inspired by natural visual pigments can feature color tunability across the full visible spectrum. However, control of excited state dynamics of the retinal chromophore, which is of paramount importance for technological applications, is lacking due to its complex and subtle photophysics/photochemistry. Here, ultrafast transient absorption spectroscopy and quantum mechanics/molecular mechanics simulations are combined for the study of highly tunable rhodopsin mimics, as compared to retinal chromophores in solution. Conical intersections and transient fluorescent intermediates are identified with atomistic resolution, providing unambiguous assignment of their ultrafast excited state absorption features. The results point out that the electrostatic environment of the chromophore, modified by protein point mutations, affects its excited state properties allowing control of its photophysics with same power of chemical modifications of the chromophore. The complex nature of such fine control is a fundamental knowledge for the design of bio-mimetic opto-electronic and photonic devices.

Ionic Atmosphere Effect on the Absorption Spectrum of a Flavoprotein: A Reminder to Consider Solution Ions

This study utilizes the FMN-dependent NADH:quinone oxidoreductase from Pseudomonas aeruginosa PAO1 to investigate the effect of introducing an active site negative charge on the flavin absorption spectrum both in the absence and presence of a long-range electrostatic potential coming from solution ions. There were no observed changes in the flavin UV-visible spectrum when an active site tyrosine (Y277) becomes deprotonated in vitro. These results could only be reproduced computationally using average solvent electrostatic configuration (ASEC) QM/MM simulations that include both positive and negative solution ions. The same calculations performed with minimal ions to neutralize the total protein charge predicted that deprotonating Y277 would significantly alter the flavin absorption spectrum. Analyzing the distribution of solution ions indicated that the ions reorganize around the protein surface upon Y277 deprotonation to cancel the effect of the tyrosinate on the flavin absorption spectrum. Additional biochemical experiments were performed to test this hypothesis.

First-Principles Nonadiabatic Dynamics Simulation of Azobenzene Photodynamics in Solutions

The photoisomerization of azobenzene is a prototypical reaction of various light-activated processes in material and biomedical sciences. However, its reaction mechanism has been under debate for decades, partly due to the challenges in computational simulations to accurately describe the molecule's photodynamics. A recent study (J. Am. Chem. Soc. 2020, 142 (49), 20,680-20,690) addressed the challenges by combining the hole-hole Tamm-Dancoff Approximated (hh-TDA) density functional theory (DFT) method with the ab initio multiple spawning (AIMS) algorithm. The hh-TDA-DFT/AIMS method was applied to first-principles nonadiabatic dynamics simulation of azobenzene's photodynamics in the vacuum. However, it remains necessary to benchmark this new method in realistic molecular environments against experimental data. In the current work, the hh-TDA-DFT/AIMS method was employed in a quantum mechanics/molecular mechanics setting to characterize the trans azobenzene's photodynamics in explicit methanol and n-hexane solvents, following both the S₁ (nπ*) and S₂ (ππ*) excitations. The simulated absorption and fluorescence spectra following the S₂ excitation quantitatively agree with the experiments. However, the hh-TDA-DFT method overestimates the torsional barrier on the S₁ state, leading to an overestimation of the S₁ state lifetime. The excited-state population decays to the ground state through two competing channels. The reactive channel partially yields the cis azobenzene photoproduct, and the unreactive channel exclusively leads to the reactant. The S₂ excitation increases the decay through the unreactive channel and thus decreases the isomerization quantum yield compared to the S₁ excitation. The solvent slows down the azobenzene's torsional dynamics on the S₁ state, but its polarity minimally affects the reaction kinetics and quantum yields. Interestingly, the dynamics of the central torsion and angles of azobenzene play a critical role in determining the final isomer of the azobenzene. This benchmark study validates the hh-TDA-DFT/AIMS method's accuracy for simulating the azobenzene's photodynamics in realistic molecular environments.

Electrostatic Control of Photoisomerization in Channelrhodopsin 2

Channelrhodopsin 2 (ChR2) is the most commonly used tool in optogenetics. Because of its faster photocycle compared to wild-type (WT) ChR2, the E123T mutant of ChR2 is a useful optogenetic tool when fast neuronal stimulation is needed. Interestingly, in spite of its faster photocycle, the initial step of the photocycle in E123T (photoisomerization of retinal protonated Schiff base or RPSB) was found experimentally to be much slower than that of WT ChR2. The E123T mutant replaces the negatively charged E123 residue with a neutral T123 residue, perturbing the electric field around the RPSB. Understanding the RPSB photoisomerization mechanism in ChR2 mutants will provide molecular-level insights into how ChR2 photochemical reactivity can be controlled, which will lay the foundation for improving the design of optogenetic tools. In this work, we combine ab initio nonadiabatic dynamics simulation, excited state free energy calculation, and reaction path search to comprehensively characterize the RPSB photoisomerization mechanism in the E123T mutant of ChR2. Our simulation agrees with previous experiments in predicting a red-shifted absorption spectrum and significant slowdown of photoisomerization in the E123T mutant. Interestingly, our simulations predict similar photoisomerization quantum yields for the mutant and WT despite the differences in excited-state lifetime and absorption maximum. Upon mutation, the neutralization of the negative charge on the E123 residue increases the isomerization barrier, alters the reaction pathway, and changes the relative stability of two fluorescent states. Our findings provide new insight into the intricate role of the electrostatic environment on the RPSB photoisomerization mechanism in microbial rhodopsins.

Directing excited state dynamics via chemical substitution: A systematic study of π-donors and π-acceptors at a carbon-carbon double bond

Functional group substituents are a ubiquitous tool in ground-state organic chemistry often employed to fine-tune chemical properties and obtain desired chemical reaction outcomes. Their effect on photoexcited electronic states, however, remains poorly understood. To help build an intuition for these effects, we have studied ethylene, substituted with electron acceptor (cyano) and/or electron donor (methoxy) substituents, both theoretically and experimentally: using ab initio quantum molecular dynamics and time-resolved photoelectron spectroscopy. Our results show the consistent trend that photo-induced ethylenic dynamics is primarily localized to the carbon with the greater electron density. For doubly substituted ethylenes, the trend is additive when both substituents are located on opposite carbons, whereas the methoxy group (in concert with steric effects) dominates when both substituents are located on a single carbon atom. These results point to the development of rules for structure-dynamics correlations; in this case, a novel mechanistic ultrafast photochemistry for conjugated carbon chains employing long-established chemical concepts.

Multistate Multiconfiguration Quantum Chemical Computation of the Two-Photon Absorption Spectra of Bovine Rhodopsin

Recently, progress in IR sources has led to the discovery that humans can detect infrared (IR) light. This is hypothesized to be due to the two-photon absorption (TPA) events promoting the retina dim-light rod photoreceptor rhodopsin to the same excited state populated via one-photon absorption (OPA). Here, we combine quantum mechanics/molecular mechanics and extended multiconfiguration quasi-degenerate perturbation theory calculations to simulate the TPA spectrum of bovine rhodopsin (Rh) as a model for the human photoreceptor. The results show that the TPA spectrum of Rh has an intense S₀ → S₁ band but shows also S₀ → S₂ and S₀ → S₃ transitions whose intensities, relative to the S₀ → S₁ band, are significantly increased when compared to the corresponding bands of the OPA spectrum. In conclusion, we show that IR light in the 950 nm region can be perceived by rod photoreceptors, thus supporting the two-photon origin of the IR perception. We also found that the same photoreceptor can perceive red (i.e., close to 680 nm) light provided that TPA induces population of S₂.

Effect of point mutations on the ultrafast photo-isomerization of Anabaena sensory rhodopsin

Anabaena sensory rhodopsin (ASR) is a particular microbial retinal protein for which light-adaptation leads to the ability to bind both the all-trans, 15-anti (AT) and the 13-cis, 15-syn (13C) isomers of the protonated Schiff base of retinal (PSBR). In the context of obtaining insight into the mechanisms by which retinal proteins catalyse the PSBR photo-isomerization reaction, ASR is a model system allowing to study, within the same protein, the protein-PSBR interactions for two different PSBR conformers at the same time. A detailed analysis of the vibrational spectra of AT and 13C, and their photo-products in wild-type ASR obtained through femtosecond (pump-) four-wave-mixing is reported for the first time, and compared to bacterio- and channelrhodopsin. As part of an extensive study of ASR mutants with blue-shifted absorption spectra, we present here a detailed computational analysis of the origin of the mutation-induced blue-shift of the absorption spectra, and identify electrostatic interactions as dominating steric effects that would entail a red-shift. The excited state lifetimes and isomerization reaction times (IRT) for the three mutants V112N, W76F, and L83Q are studied experimentally by femtosecond broadband transient absorption spectroscopy. Interestingly, in all three mutants, isomerization is accelerated for AT with respect to wild-type ASR, and this the more, the shorter the wavelength of maximum absorption. On the contrary, the 13C photo-reaction is slightly slowed down, leading to an inversion of the ESLs of AT and 13C, with respect to wt-ASR, in the blue-most absorbing mutant L83Q. Possible mechanisms for these mutation effects, and their steric and electrostatic origins are discussed.

Direct evidence for hula twist and single-bond rotation photoproducts

Photoisomerization reactions are quintessential processes driving molecular machines and motors, govern smart materials, catalytic processes, and photopharmacology, and lie at the heart of vision, phototaxis, or vitamin production. Despite this plethora of applications fundamental photoisomerization mechanisms are not well understood at present. The famous hula-twist motion—a coupled single and double-bond rotation—was proposed to explain proficient photoswitching in restricted environments but fast thermal follow-up reactions hamper identification of primary photo products. Herein we describe an asymmetric chromophore possessing four geometrically distinct diastereomeric states that do not interconvert thermally and can be crystallized separately. Employing this molecular setup direct and unequivocal evidence for the hula-twist photoreaction and for photoinduced single-bond rotation is obtained. The influences of the surrounding medium and temperature are quantified and used to favor unusual photoreactions. Based on our findings molecular engineers will be able to implement photo control of complex molecular motions more consciously.

Photoisomerization mechanisms govern important (bio)catalytic reactions and lie at the core of many functional materials. Here, the authors report a molecular setup that allows for the direct and separate observation of three fundamental photoreactions, namely the hula twist, single-bond rotation, as well as double-bond isomerization.

How Methylation Modifies the Photophysics of the Native All- trans-Retinal Protonated Schiff Base: A CASPT2/MD Study in Gas Phase and in Methanol

A comparison between the free-energy surfaces of the all- trans-retinal protonated Schiff base (RPSB) and its 10-methylated derivative in gas phase and methanol solution is performed at CASSCF//CASSCF and CASPT2//CASSCF levels. Solvent effects were included using the average solvent electrostatic potential from molecular dynamics method. This is a QM/MM (quantum mechanics/molecular mechanics) method that makes use of the mean field approximation. It is found that the methyl group bonded to C10 produces noticeable changes in the solution free-energy profile of the S₁ excited state, mainly in the relative stability of the minimum energy conical intersections (MECIs) with respect to the Franck-Condon (FC) point. The conical intersections yielding the 9- cis and 11- cis isomers are stabilized while that yielding the 13- cis isomer is destabilized; in fact, it becomes inaccessible by excitation to S₁. Furthermore, the planar S₁ minimum is not present in the methylated compound. The solvent notably stabilizes the S₂ excited state at the FC geometry. Therefore, if the S₂ state has an effect on the photoisomerization dynamics, it must be because it permits the RPSB population to branch around the FC point. All these changes combine to speed up the photoisomerization in the 10-methylated compound with respect to the native compound.

Impact of Electronic State Mixing on the Photoisomerization Time Scale of the Retinal Chromophore

Spectral data show that the photoisomerization of retinal protonated Schiff base (rPSB) chromophores occurs on a 100 fs time scale or less in vertebrate rhodopsins, it is several times slower in microbial rhodopsins and it is between one and 2 orders of magnitude slower in solution. These time scale variations have been attributed to specific modifications of the topography of the first excited state potential energy surface of the chromophore. However, it is presently not clear which specific environment effects (e.g., electrostatic, electronic, or steric) are responsible for changing the surface topography. Here, we use QM/MM models and excited state trajectory computations to provide evidence for an increase in electronic mixing between the first and the second excited state of the chromophore when going from vertebrate rhodopsin to the solution environments. Ultimately, we argue that a correlation between the lifetime of the first excited state and electronic mixing between such state and its higher neighbor, may have been exploited to evolve rhodopsins toward faster isomerization and, possibly, light-sensitivity.

New Perspectives on an Old Issue: A Comparative MS-CASPT2 and OM2-MRCI Study of Polyenes and Protonated Schiff Bases

Polyenic systems are involved in light perception of numerous living organisms. Although a π-conjugated backbone is a common feature of all polyenes, their photophysics may vary. We provide a comparative quantum mechanical study of low-lying S₁ and S₂ excited states in short (3-5 double bonds) symmetric all-trans linear polyenes and corresponding protonated Schiff bases. In our investigation, we use the well-established ab initio multireference CASPT2 approach and benchmark the efficient semiempirical OM2-MRCI approach against it. For all protonated Schiff bases, MS-CASPT2 results in two distinct S₁ minima with inverted and noninverted bond length pattern, respectively. We find that OM2-MRCI is a computationally affordable and reliable alternative to MS-CASPT2 for investigations of polyenic systems, particularly when highly demanding calculations (e.g. excited-state dynamics) need to be performed.

Fine Tuning of Retinal Photoinduced Decay in Solution

Single methylation at position C₁₀ of the all-trans retinal protonated Schiff base switches its excited-state decay in methanol from a slower picosecond into an ultrafast, protein-like subpicosecond process. QM/MM modeling in conjunction with on-the-fly excited-state dynamics provides fundamental understanding of the fine-tuning mechanics that "catalyzes" the photoinduced decay of solvated retinals. Methylation alters the interplay between the ionic S₁ and covalent S₂ states, reducing the excited-state lifetime by favoring the formation of a S₁ transient fluorescent state with fully inverted bond lengths that accounts for the recorded transient spectroscopy and from which a space-saving conical intersection seam is quickly (<1 ps) reached. Minimal and apparently innocent chemical modifications thus affect the characteristic intramolecular charge-transfer of the S₁ state as well as the interaction with the covalent S₂ excited state, eventually providing the high tunability of retinal photophysics and photochemistry and delivering a new concept for the rational design of retinal-based photoactive molecular devices.

Mechanism of ultrafast non-reactive deactivation of the retinal chromophore in non-polar solvents

Counterion sensitive photodynamics of the retinal chromophore in solution.

Exploring the isomerization paths of push–pull hexatrienes

The first photo-excited singlet state of all-trans donor (amine), acceptor (cyano)-substituted hexatriene forms donor side-twisted minima and its passage towards the cis–trans–trans isomer is characterized by an S₀/S₁ conical intersection.

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.

Excited States and Photochemistry of Chromophores in the Photoactive Proteins Explored by the Combined Quantum Mechanical and Molecular Mechanical Calculations

A photoactive protein usually contains a unique chromophore that is responsible for the initial photoresponse and functions of the photoactive protein are determined by the interaction between the chromophore and its protein surroundings. The combined quantum mechanical and molecular mechanical (QM/MM) approach is demonstrated to be a very useful tool for exploring structures and functions of a photoactive protein with the chromophore and its protein surroundings treated by the QM and MM methods, respectively. In this review, we summarize the basic formulas of the QM/MM approach and emphasize its applications to excited states and photoreactions of chromophores in rhodopsin protein, photoactive yellow protein, and green fluorescent protein.

100 fs photo-isomerization with vibrational coherences but low quantum yield in Anabaena Sensory Rhodopsin

Counter-intuitive photochemistry: in Anabaena Sensory Rhodopsin, the retinal 13-cis isomer isomerizes much faster than all-trans ASR, but with a 3-times lower quantum yield.

Tracking the primary photoconversion events in rhodopsins by ultrafast optical spectroscopy

We review the most recent experimental and computational efforts aimed at exposing the very early phases of the ultrafast isomerization in visual Rhodopsins and we discuss future advanced experiments and calculations.

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.

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.