The Ultrafast Photoisomerizations of Rhodopsin and Bathorhodopsin Are Modulated by Bond Length Alternation and HOOP Driven Electronic Effects

Including Tunneling in Non-Born-Oppenheimer Simulations

For electronically nonadiabatic processes in all but the simplest systems, the most practical multidimensional simulation method is a semiclassical approximation in which a trajectory or the center of a wave packet follows a classical path governed by an effective potential energy function. Here, we show how such simulations can be made more realistic by including tunneling by the army ants tunneling method. We illustrate the theory by calculations with model potential energy surfaces; one model study is in the adiabatic limit, and the other one has nonadiabatic transitions between two electronic states during the tunneling event. The army ants tunneling algorithm is used to efficiently sample tunneling events in the trajectories in both cases. This work makes it possible to simulate complex nonadiabatic chemical processes by efficiently including the important quantum effect of tunneling.

Comparison of the isomerization mechanisms of human melanopsin and invertebrate and vertebrate rhodopsins

Comparative modeling and ab initio multiconfigurational quantum chemistry are combined to investigate the reactivity of the human nonvisual photoreceptor melanopsin. It is found that both the thermal and photochemical isomerization of the melanopsin 11-cis retinal chromophore occur via a space-saving mechanism involving the unidirectional, counterclockwise twisting of the =C11H-C12H= moiety with respect to its Lys340-linked frame as proposed by Warshel for visual pigments [Warshel A (1976) Nature 260(5553):679-683]. A comparison with the mechanisms documented for vertebrate (bovine) and invertebrate (squid) visual photoreceptors shows that such a mechanism is not affected by the diversity of the three chromophore cavities. Despite such invariance, trajectory computations indicate that although all receptors display less than 100 fs excited state dynamics, human melanopsin decays from the excited state ∼40 fs earlier than bovine rhodopsin. Some diversity is also found in the energy barriers controlling thermal isomerization. Human melanopsin features the highest computed barrier which appears to be ∼2.5 kcal mol(-1) higher than that of bovine rhodopsin. When assuming the validity of both the reaction speed/quantum yield correlation discussed by Warshel, Mathies and coworkers [Weiss RM, Warshel A (1979) J Am Chem Soc 101:6131-6133; Schoenlein RW, Peteanu LA, Mathies RA, Shank CV (1991) Science 254(5030):412-415] and of a relationship between thermal isomerization rate and thermal activation of the photocycle, melanopsin turns out to be a highly sensitive pigment consistent with the low number of melanopsin-containing cells found in the retina and with the extraretina location of melanopsin in nonmammalian vertebrates.

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.

Retinal conformation governs pKa of protonated Schiff base in rhodopsin activation

We have explored the relationship between conformational energetics and the protonation state of the Schiff base in retinal, the covalently bound ligand responsible for activating the G protein-coupled receptor rhodopsin, using quantum chemical calculations. Guided by experimental structural determinations and large-scale molecular simulations on this system, we examined rotation about each bond in the retinal polyene chain, for both the protonated and deprotonated states that represent the dark and photoactivated states, respectively. Particular attention was paid to the torsional degrees of freedom that determine the shape of the molecule, and hence its interactions with the protein binding pocket. While most torsional degrees of freedom in retinal are characterized by large energetic barriers that minimize structural fluctuations under physiological temperatures, the C6-C7 dihedral defining the relative orientation of the β-ionone ring to the polyene chain has both modest barrier heights and a torsional energy surface that changes dramatically with protonation of the Schiff base. This surprising coupling between conformational degrees of freedom and protonation state is further quantified by calculations of the pKa as a function of the C6-C7 dihedral angle. Notably, pKa shifts of greater than two units arise from torsional fluctuations observed in molecular dynamics simulations of the full ligand-protein-membrane system. It follows that fluctuations in the protonation state of the Schiff base occur prior to forming the activated MII state. These new results shed light on important mechanistic aspects of retinal conformational changes that are involved in the activation of rhodopsin in the visual process.

Photocycle and vectorial proton transfer in a rhodopsin from the eukaryote Oxyrrhis marina

Retinylidene photoreceptors are ubiquitously present in marine protists as first documented by the identification of green proteorhodopsin (GPR). We present a detailed investigation of a rhodopsin from the protist Oxyrrhis marina (OR1) with respect to its spectroscopic properties and to its vectorial proton transport. Despite its homology to GPR, OR1's features differ markedly in its pH dependence. Protonation of the proton acceptor starts at pH below 4 and is sensitive to the ionic conditions. The mutation of a conserved histidine H62 did not influence the pK(a) value in a similar manner as in other proteorhodopsins where the charged histidine interacts with the proton acceptor forming the so-called His-Asp cluster. Mutational and pH-induced effects were further reflected in the temporal behavior upon light excitation ranging from femtoseconds to seconds. The primary photodynamics exhibits a high sensitivity to the environment of the proton acceptor D100 that are correlated to the different initial states. The mutation of the H62 does not affect photoisomerization at neutral pH. This is in agreement with NMR data indicating the absence of the His-Asp cluster. The subsequent steps in the photocycle revealed protonation reactions at the Schiff base coupled to proton pumping even at low pH. The main electrogenic steps are associated with the reprotonation of the Schiff base and internal proton donor. Hence, OR1 shows a different theme of the His-Asp organization where the low pK(a) of the proton acceptor is not dominated by this interaction, but by other electrostatic factors.

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.

Protein Field Effect on the Dark State of 11-cis Retinal in Rhodopsin by Quantum Monte Carlo/Molecular Mechanics

The accurate determination of the geometrical details of the dark state of 11-cis retinal in rhodopsin represents a fundamental step for the rationalization of the protein role in the optical spectral tuning in the vision mechanism. We have calculated geometries of the full retinal protonated Schiff base chromophore in the gas phase and in the protein environment using the correlated variational Monte Carlo method. The bond length alternation of the conjugated carbon chain of the chromophore in the gas phase shows a significant reduction when moving from the β-ionone ring to the nitrogen, whereas, as expected, the protein environment reduces the electronic conjugation. The proposed dark state structure is fully compatible with solid-state NMR data reported by Carravetta et al. [J. Am. Chem. Soc. 2004, 126, 3948-3953]. TDDFT/B3LYP calculations on such geometries show a blue opsin shift of 0.28 and 0.24 eV induced by the protein for S1 and S2 states, consistently with literature spectroscopic data. The effect of the geometrical distortion alone is a red shift of 0.21 and 0.16 eV with respect to the optimized gas phase chromophore. Our results open new perspectives for the study of the properties of chromophores in their biological environment using correlated methods.

Computational Photochemistry and Photophysics: the state of the art

This review starts with the most basic concepts in photochemistry and photophysics, followed by a chronological introduction of theoretical methods and relevant applications in the history of computational photochemistry, along with the authors’ comments on the methodologies currently available for photochemical studies. Recent advances in the field are next summarized and discussed, focusing separately on methodology and computational techniques and some highlighted applied works carried out during the last two years on the topics of photodissociations, photostability, photodimerizations, photoisomerizations, proton/hydrogen transfer, photodecarboxylations, charge transport, bioexcimers, chemiluminescence and bioluminescence. We finish this review by conclusions and an outlook of the future.

A large geometric distortion in the first photointermediate of rhodopsin, determined by double-quantum solid-state NMR

Double-quantum magic-angle-spinning NMR experiments were performed on 11,12-(13)C(2)-retinylidene-rhodopsin under illumination at low temperature, in order to characterize torsional angle changes at the C11-C12 photoisomerization site. The sample was illuminated in the NMR rotor at low temperature (~120 K) in order to trap the primary photointermediate, bathorhodopsin. The NMR data are consistent with a strong torsional twist of the HCCH moiety at the isomerization site. Although the HCCH torsional twist was determined to be at least 40°, it was not possible to quantify it more closely. The presence of a strong twist is in agreement with previous Raman observations. The energetic implications of this geometric distortion are discussed.

Tuning the primary reaction of channelrhodopsin-2 by imidazole, pH, and site-specific mutations

Femtosecond time-resolved absorption measurements were performed to investigate the influence of the pH, imidazole concentration, and point mutations on the isomerization process of Channelrhodopsin-2. Apart from the typical spectral characteristics of retinal isomerization, an additional absorption feature rises for the wild-type (wt) on a timescale from tens of ps to 1 ns within the spectral range of the photoproduct and is attributed to an equilibration between different K-intermediates. Remarkably, this absorption feature vanishes upon addition of imidazole or lowering the pH. In the latter case, the isomerization is dramatically slowed down, due to protonation of negatively charged amino acids within the retinal binding pocket, e.g., E123 and D253. Moreover, we investigated the influence of several point mutations within the retinal binding pocket E123T, E123D, C128T, and D156C. For E123T, the isomerization is retarded compared to wt and E123D, indicating that a negatively charged residue at this position functions as an effective catalyst in the isomerization process. In the case of the C128T mutant, all primary processes are slightly accelerated compared to the wt, whereas the isomerization dynamics for the D156C mutant is similar to wt after addition of imidazole.

Backbone modification of retinal induces protein-like excited state dynamics in solution

The drastically different reactivity of the retinal chromophore in solution compared to the protein environment is poorly understood. Here, we show that the addition of a methyl group to the C═C backbone of all-trans retinal protonated Schiff base accelerates the electronic decay in solution making it comparable to the proton pump bacteriorhodopsin. Contrary to the notion that reaction speed and efficiency are linked, we observe a concomitant 50% reduction in the isomerization yield. Our results demonstrate that minimal synthetic engineering of potential energy surfaces based on theoretical predictions can induce drastic changes in electronic dynamics toward those observed in an evolution-optimized protein pocket.