Computational evidence in favor of a two-state, two-mode model of the retinal chromophore photoisomerization

Molecular Switching in Confined Spaces: Effects of Encapsulating the DHA/VHF Photo-Switch in Cucurbiturils

Confinement of reactive chemical species uniquely affects chemical reactivity by restricting the physical space available and by restricting access to interactions with the solvent. In Nature, for example, confined protein binding pockets govern processes following photoisomerization reactions and the isomerizations themselves. Here we describe the first example of a dihydroazulene/vinylheptafulvene (DHA/VHF) photo-switch functioning in water, and we show how its switching behavior is strongly influenced by supramolecular interactions with a series of cucurbit[n]uril (CB) host molecules. In CB7 inclusion complexes, the kinetics of the thermal VHF-to-DHA back-reaction is accelerated, while in CB8 inclusion complexes, the kinetics is slowed down as compared to the free photo-switch. The effect of the CB encapsulation of the photo-switch can be effectively canceled by introducing a guest that binds the CB more strongly. According to DFT calculations, a stabilization of the reactive s-cis VHF conformer relative to the s-trans VHF appears to be a contributing factor responsible for the accelerated back-reaction when encapsulated in CB7.

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.

Schiff Base Proton Acceptor Assists Photoisomerization of Retinal Chromophores in Bacteriorhodopsin

In this study, we investigated the ultrafast dynamics of bacteriorhodopsins (BRs) from Haloquadratum walsbyi (HwBR) and Haloarcula marismortui (HmBRI and HmBRII). First, the ultrafast dynamics were studied for three HwBR samples: wild-type, D93N mutation, and D104N mutation. The residues of the D93 and D104 mutants correspond to the control by the Schiff base proton acceptor and donor of the proton translocation subchannels. Measurements indicated that the negative charge from the Schiff base proton acceptor residue D93 interacts with the ultrafast and substantial change of the electrostatic potential associated with chromophore isomerization. By contrast, the Schiff base proton donor assists the restructuring of the chromophore cavity hydrogen-bond network during the thermalization of the vibrational hot state. Second, the ultrafast dynamics of the wild-types of HwBR, HmBRI, and HmBRII were compared. Measurements demonstrated that the hydrogen-bond network in the extracellular region in HwBR and HmBRII slows the photoisomerization of retinal chromophores, and the negatively charged helices on the cytoplasmic side of HwBR and HmBRII accelerate the thermalization of the vibrational hot state of retinal chromophores. The similarity of the correlation spectra of the wild-type HmBRI and D104N mutant of HwBR indicates that inactivation of the Schiff base proton donor induces a positive charge on the helices of the cytoplasmic side.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

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

Counterion sensitive photodynamics of the retinal chromophore in solution.

Femtosecond spectroscopic study of photochromic reactions of bacteriorhodopsin and visual rhodopsin

Photochromic ultrafast reactions of bacteriorhodopsin (H. salinarum) and bovine rhodopsin were conducted with a femtosecond two-pump probe pulse setup with the time resolution of 20-25fs. The dynamics of the forward and reverse photochemical reactions for both retinal-containing proteins was compared. It is demonstrated that when retinal-containing proteins are excited by femtosecond pulses, dynamics pattern of the vibrational coherent wave packets in the course of the reaction is different for bacteriorhodopsin and visual rhodopsin. As shown in these studies, the low-frequencies that form a wave packets experimentally observed in the dynamics of primary products formation as a result of retinal photoisomerization have different intensities and are clearer for bovine rhodopsin. Photo-reversible reactions for both retinal proteins were performed from the stage of the relatively stable photointermediates that appear within 3-5ps after the light pulse impact. It is demonstrated that the efficiency of the reverse phototransition K-form→bacteriorhodopsin is almost five-fold higher than that of the Batho-intermediate→visual rhodopsin phototransition. The results obtained indicate that in the course of evolution the intramolecular mechanism of the chromophore-protein interaction in visual rhodopsin becomes more perfect and specific. The decrease in the probability of the reverse chromophore photoisomerization (all-trans→11-cis retinal) in primary photo-induced rhodopsin products causes an increase in the efficiency of the photoreception process.

Probing the Interstate Coupling near a Conical Intersection by Optical Spectroscopy

Conical intersections are points where adiabatic potential energy surfaces cross. The interstate coupling between the potential energy surfaces plays a crucial role in many processes associated with conical intersections. Still no method exists to measure this coupling driving the chemical reactions between the potential energy surfaces involved. In this Letter, using a generic model for photoisomerization, we propose a novel experimental approach to estimate the coupling that mixes the electronic states near a conical intersection. The approach is based on analyzing the vibrational wavepacket of the reactant in the adiabatic ground and excited electronic states. The nuclear wavepacket dynamics are extracted from linear absorption and two-dimensional electronic spectroscopy. Comparing the frequencies of the coupling mode in the adiabatic ground and excited states from models with and without coupling between the potential energy surfaces suggests an experimental tool to determine the interstate coupling.

Impact of Vibrational Coherence on the Quantum Yield at a Conical Intersection

We study the vibrationally coherent quantum dynamics of an electronic wave packet in the vicinity of a conical intersection within a three-state two-mode model. By transforming the coherent tuning and coupling modes into the bath, the underdamped dynamics of the resulting effective three-state model is solved efficiently by the numerically exact hierarchy equation of motion approach. The transient excited-state absorption and two-dimensional spectra reveal the impact of vibrational coherence on the relaxation pathways of the wave packet. We find that both the quantum yield and the isomerization rate are crucially influenced by the vibrational coherence of the wave packet. A less coherent wave packet can traverse the conical intersection more rapidly, while the resulting quantum yield is smaller. Finally, we show that repeated passages of the wave packet through the conical intersection can lead to measurable interference effects in the form of Stueckelberg oscillations.

Ultrafast Relaxation Dynamics of the Ethylene Cation C(2)H(4)+

We present a combined experimental and computational study of the relaxation dynamics of the ethylene cation. In the experiment, we apply an extreme-ultraviolet-pump/infrared-probe scheme that permits us to resolve time scales on the order of 10 fs. The photoionization of ethylene followed by an infrared (IR) probe pulse leads to a rich structure in the fragment ion yields reflecting the fast response of the molecule and its nuclei. The temporal resolution of our setup enables us to pinpoint an upper bound of the previously defined ethylene-ethylidene isomerization time to 30 ± 3 fs. Time-dependent density functional based trajectory surface hopping simulations show that internal relaxation between the first excited states and the ground state occurs via three different conical intersections. This relaxation unfolds on femtosecond time scales and can be probed by ultrashort IR pulses. Through this probe mechanism, we demonstrate a route to optical control of the important dissociation pathways leading to separation of H or H2.

Opsin Effect on the Electronic Structure of the Retinylidene Chromophore in Rhodopsin

Direct examination of experimental NMR parameters combined with electronic structure analysis was used to provide a first-principle interpretation of NMR experiments and give a precise evaluation of how the electronic perturbation of the protein environment affects the electronic properties of the retinylidene chromophere in rhodopsin. To this end, we pursued a theoretical analysis using a combination of tools including quantum mechanics/molecular mechanics (QM/MM) at the Density Functional Theory (DFT) level, in conjunction with gauge independent atomic orbital (GIAO) calculations of (13)C NMR chemical shieldings and (1)J(CC) spin-spin coupling constants obtained with the Coupled Perturbed DFT (CPDFT) method. The opsin effect on the retinylidene chromophere is interpreted as an inductive effect of Glu-113 which readjusts the weighting factors of resonance substructures of the conjugated chain of the chromophere. These changes give a rationalization to the alternating effect of the (13)C chemical shifts magnitudes when comparing the retinylidene chromophere in the presence and absence of the protein environment. Conversely, perturbation of π orbitals has little to no effect over (1)J (13)C-(13)C spin-spin coupling constants, as they are mainly dominated by the Fermi contact term, and hence the counteraion effect is restricted to the vicinity of the perturbation. Thus, the apparent contradiction between experimental findings based on chemical shifts (deep penetration) and one-bond J-couplings (localized effects of the protonated Schiff base at the chain terminus) is in fact a consequence of different properties responding differently to the same external perturbation.

Mechanisms responsible for the photocurrent in bacteriorhodopsin

Recently, there has been growing interest in the electrical properties of bacteriorhodopsin (bR), a protein belonging to the transmembrane protein family. Several experiments pointed out the role of green light in enhancing the current flow in nanolayers of bR, thus confirming potential applications of this protein in the field of optoelectronics. By contrast, the mechanisms underlying the charge transfer and the associated photocurrent are still far from being understood at a microscopic level. To take into account the structure-dependent nature of the current, in a previous set of papers we suggested a mechanism of sequential tunneling among neighboring amino acids. As a matter of fact, when irradiated with green light, bR undergoes a conformational change at a molecular level. Thus, the role played by the protein tertiary-structure in modeling the charge transfer cannot be neglected. The aim of this paper is to go beyond previous models, in the framework of a new branch of electronics we call proteotronics, which exploits the ability of using proteins as reliable, well-understood materials for the development of novel bioelectronic devices. In particular, the present approach assumes that the conformational change is not the unique transformation the protein undergoes when irradiated by light. Instead, the light can also promote an increase of the protein state free energy that, in turn, should modify its internal degree of connectivity. This phenomenon is here described by the change of the value of an interaction radius associated with the physical interactions among amino acids. The implemented model enables us to achieve a better agreement between theory and experiments in the region of a low applied bias by preserving the level of agreement at high values of applied bias. Furthermore, results provide new insights on the mechanisms responsible for bR photoresponse.

The G protein-coupled receptor rhodopsin: a historical perspective

Rhodopsin is a key light-sensitive protein expressed exclusively in rod photoreceptor cells of the retina. Failure to express this transmembrane protein causes a lack of rod outer segment formation and progressive retinal degeneration, including the loss of cone photoreceptor cells. Molecular studies of rhodopsin have paved the way to understanding a large family of cell-surface membrane proteins called G protein-coupled receptors (GPCRs). Work started on rhodopsin over 100 years ago still continues today with substantial progress made every year. These activities underscore the importance of rhodopsin as a prototypical GPCR and receptor required for visual perception-the fundamental process of translating light energy into a biochemical cascade of events culminating in vision.

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.

Retinal shows its true colours: photoisomerization action spectra of mobility-selected isomers of the retinal protonated Schiff base

Isomers of the retinal protonated Schiff base are separated and probed using laser radiation in a tandem ion mobility spectrometer yielding isomer-specific electronic spectra.

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.

Human infrared vision is triggered by two-photon chromophore isomerization

Vision relies on photoactivation of visual pigments in rod and cone photoreceptor cells of the retina. The human eye structure and the absorption spectra of pigments limit our visual perception of light. Our visual perception is most responsive to stimulating light in the 400- to 720-nm (visible) range. First, we demonstrate by psychophysical experiments that humans can perceive infrared laser emission as visible light. Moreover, we show that mammalian photoreceptors can be directly activated by near infrared light with a sensitivity that paradoxically increases at wavelengths above 900 nm, and display quadratic dependence on laser power, indicating a nonlinear optical process. Biochemical experiments with rhodopsin, cone visual pigments, and a chromophore model compound 11-cis-retinyl-propylamine Schiff base demonstrate the direct isomerization of visual chromophore by a two-photon chromophore isomerization. Indeed, quantum mechanics modeling indicates the feasibility of this mechanism. Together, these findings clearly show that human visual perception of near infrared light occurs by two-photon isomerization of visual pigments.

Combination of transient 2D-IR experiments and ab initio computations sheds light on the formation of the charge-transfer state in photoexcited carbonyl carotenoids

The excited state dynamics of carbonyl carotenoids is very complex because of the coupling of single- and doubly excited states and the possible involvement of intramolecular charge-transfer (ICT) states. In this contribution we employ ultrafast infrared spectroscopy and theoretical computations to investigate the relaxation dynamics of trans-8'-apo-β-carotenal occurring on the picosecond time scale, after excitation in the S2 state. In a (slightly) polar solvent like chloroform, one-dimensional (T1D-IR) and two-dimensional (T2D-IR) transient infrared spectroscopy reveal spectral components with characteristic frequencies and lifetimes that are not observed in nonpolar solvents (cyclohexane). Combining experimental evidence with an analysis of CASPT2//CASSCF ground and excited state minima and energy profiles, complemented with TDDFT calculations in gas phase and in solvent, we propose a photochemical decay mechanism for this system where only the bright single-excited 1Bu(+) and the dark double-excited 2Ag(-) states are involved. Specifically, the initially populated 1Bu(+) relaxes toward 2Ag(-) in 200 fs. In a nonpolar solvent 2Ag(-) decays to the ground state (GS) in 25 ps. In polar solvents, distortions along twisting modes of the chain promote a repopulation of the 1Bu(+) state which then quickly relaxes to the GS (18 ps in chloroform). The 1Bu(+) state has a high electric dipole and is the main contributor to the charge-transfer state involved in the dynamics in polar solvents. The 2Ag(-) → 1Bu(+) population transfer is evidenced by a cross peak on the T2D-IR map revealing that the motions along the same stretching of the conjugated chain on the 2Ag(-) and 1Bu(+) states are coupled.

Shape of Multireference, Equation-of-Motion Coupled-Cluster, and Density Functional Theory Potential Energy Surfaces at a Conical Intersection

We report and characterize ground-state and excited-state potential energy profiles using a variety of electronic structure methods along a loop lying on the branching plane associated with a conical intersection (CI) of a reduced retinal model, the penta-2,4-dieniminium cation (PSB3). Whereas the performance of the equation-of-motion coupled-cluster, density functional theory, and multireference methods had been tested along the excited- and ground-state paths of PSB3 in our earlier work, the ability of these methods to correctly describe the potential energy surface shape along a CI branching plane has not yet been investigated. This is the focus of the present contribution. We find, in agreement with earlier studies by others, that standard time-dependent DFT (TDDFT) does not yield the correct two-dimensional (i.e., conical) crossing along the branching plane but rather a one-dimensional (i.e., linear) crossing along the same plane. The same type of behavior is found for SS-CASPT2(IPEA=0), SS-CASPT2(IPEA=0.25), spin-projected SF-TDDFT, EOM-SF-CCSD, and, finally, for the reference MRCISD+Q method. In contrast, we found that MRCISD, CASSCF, MS-CASPT2(IPEA=0), MS-CASPT2(IPEA=0.25), XMCQDPT2, QD-NEVPT2, non-spin-projected SF-TDDFT, and SI-SA-REKS yield the expected conical crossing. To assess the effect of the different crossing topologies (i.e., linear or conical) on the PSB3 photoisomerization efficiency, we discuss the results of 100 semiclassical trajectories computed by CASSCF and SS-CASPT2(IPEA=0.25) for a PSB3 derivative. We show that for the same initial conditions, the two methods yield similar dynamics leading to isomerization quantum yields that differ by only a few percent.

Modelling time-resolved two-dimensional electronic spectroscopy of the primary photoisomerization event in rhodopsin

Time-resolved two-dimensional (2D) electronic spectra (ES) tracking the evolution of the excited state manifolds of the retinal chromophore have been simulated along the photoisomerization pathway in bovine rhodopsin, using a state-of-the-art hybrid QM/MM approach based on multiconfigurational methods. Simulations of broadband 2D spectra provide a useful picture of the overall detectable 2D signals from the near-infrared (NIR) to the near-ultraviolet (UV). Evolution of the stimulated emission (SE) and excited state absorption (ESA) 2D signals indicates that the S1 → SN (with N ≥ 2) ESAs feature a substantial blue-shift only after bond inversion and partial rotation along the cis → trans isomerization angle, while the SE rapidly red-shifts during the photoinduced skeletal relaxation of the polyene chain. Different combinations of pulse frequencies are proposed in order to follow the evolution of specific ESA signals. These include a two-color 2DVis/NIR setup especially suited for tracking the evolution of the S1 → S2 transitions that can be used to discriminate between different photochemical mechanisms of retinal photoisomerization as a function of the environment. The reported results are consistent with the available time-resolved pump-probe experimental data, and may be used for the design of more elaborate transient 2D electronic spectroscopy techniques.

The EF loop in green proteorhodopsin affects conformation and photocycle dynamics

The proteorhodopsin family consists of retinal proteins of marine bacterial origin with optical properties adjusted to their local environments. For green proteorhodopsin, a highly specific mutation in the EF loop, A178R, has been found to cause a surprisingly large redshift of 20 nm despite its distance from the chromophore. Here, we analyze structural and functional consequences of this EF loop mutation by time-resolved optical spectroscopy and solid-state NMR. We found that the primary photoreaction and the formation of the K-like photo intermediate is almost pH-independent and slower compared to the wild-type, whereas the decay of the K-intermediate is accelerated, suggesting structural changes within the counterion complex upon mutation. The photocycle is significantly elongated mainly due to an enlarged lifetime of late photo intermediates. Multidimensional MAS-NMR reveals mutation-induced chemical shift changes propagating from the EF loop to the chromophore binding pocket, whereas dynamic nuclear polarization-enhanced (13)C-double quantum MAS-NMR has been used to probe directly the retinylidene conformation. Our data show a modified interaction network between chromophore, Schiff base, and counterion complex explaining the altered optical and kinetic properties. In particular, the mutation-induced distorted structure in the EF loop weakens interactions, which help reorienting helix F during the reprotonation step explaining the slower photocycle. These data lead to the conclusion that the EF loop plays an important role in proton uptake from the cytoplasm but our data also reveal a clear interaction pathway between the EF loop and retinal binding pocket, which might be an evolutionary conserved communication pathway in retinal proteins.

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.

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.