Absorption of Schiff-Base Retinal Chromophores in Vacuo

A Set of Quantum-Mechanically Derived Force Fields for Natural and Synthetic Retinal Photoswitches

The diverse biological functions of rhodopsins are all triggered by the photoexcitation of retinal protonated Schiff base chromophores. This diversity can be traced back not only to variations in protein scaffolds in which the chromophore is embedded, but also to the different isomeric forms of the chromophore itself, whose role is crucial in several processes. Although most computational approaches for these systems often require classical molecular dynamics, efforts in providing a set of parameters able to accurately and consistently model several isomeric chromophores are lacking in the literature. The most recent efforts entail either refinements of general purpose force fields lacking in accuracy, or parametrization strategies that include environmental effects, which makes the resulting parameters not transferable to a different embedding. In this work, we provide accurate intramolecular force fields based on data purposely computed using Møller-Plesset second order perturbation theory, specifically tailored for varied natural retinal protonated Schiff bases and synthetic analogues often employed in retinal-based photoswitches. We demonstrate the quality of our quantum-mechanically derived force fields (QMD-FFs) through a wide set of validation tests. These consistently indicate that QMD-FFs outperform in all cases transferable, general-purpose FFs, delivering an excellent description of each chromophore in terms of equilibrium geometries, conformational landscapes, and optical properties in comparison to literature data, experimental measurements, and reference QM calculations. Our intramolecular QMD-FFs, distributed in electronic format, can be adopted to describe these chromophores in complex environments, exploiting intermolecular parameters compatible with those available in the literature for biological macromolecules.

In Silico Study of a Bacteriorhodopsin/TiO2 Hybrid System at the Molecular Level

Bacteriorhodopsin (bR) is a light-harvesting membrane protein that represents a promising sensitizer of TiO2 for photovoltaic and photoelectrochemical devices. However, despite numerous experimental studies, the molecular-level understanding of the bR/TiO2 hybrid system is still unsatisfactory. In this contribution, we report the construction and analysis of an atomistic model of such a system. To do so, both steered molecular dynamics-molecular dynamics and quantum mechanics/molecular mechanics computations are applied to four different bR orientations on the anatase TiO2 surface. The resulting bR/TiO2 models are then used to compute the light absorption maxima changes relative to those of bR. We show that all four models reproduce the experimentally observed blue-shift value induced by bR binding on TiO2 and could be used to study the binding and binding-induced protein modifications. We conclude that the constructed models could provide a basis for future studies aiming to simulate the complex long-range electron transfer mechanism in bR/TiO2-based solar energy conversion devices as well as in engineering bR to achieve enhanced efficiencies.

Mechanical Modulation of S0-S1 and S0-T1 Energy Gaps of 11-cis and All-trans Retinal Schiff Bases

The retinal Schiff base is a chromophore of significant biological relevance, as it is responsible for capturing sunlight in rhodopsins, which are photoactive proteins found in various living organisms. Additionally, this chromophore is subjected to various mechanical forces in different proteins, which alter its structure and, consequently, its properties. To thoroughly understand the mechanical response limits of the retinal excitation energy, a simple first-order formalism has been developed to quantify the chromophore's optimal mechanical response to applied external forces (on the order of tens of pN). Additionally, the response to larger forces is analyzed by using an algorithm to explore the potential energy surfaces. It can be concluded that the retinal Schiff base exhibits a significant mechanical response and that the optimal forces and displacements involve certain coordinates typically of low frequency, showing differences between the S1 and T1 states, as well as between the 11-cis and all-trans isomers. Additionally, the possibility of mechanically modulating the bond length alternation using mechanical forces is ruled out.

Absorption and Excited-State Coherences of Cryogenically Cold Retinal Protonated Schiff Base in Vacuo

Retinal protonated Schiff base (RPSB), found in its all-trans conformer in Bacteriorhodopsin, undergoes barrier-controlled isomerization upon photoabsorption through polyene chain torsion. The effects of the protein environment on the active vibrations during photoabsorption and their redistribution are still not understood. This paper reports on femtosecond time-resolved action-absorption measurements of cryogenically cooled gas-phase all-trans RPSB, which exhibit two coherent vibrational oscillations, 167(14) cm- 1 ${^{ - 1} }$ and 117(1) cm- 1 ${^{ - 1} }$ , of the first excited state with dephasing times of∼ 1 ${{\rm{ \sim }}1}$ ps. The absence of the high-frequency vibration in solution and the low-frequency vibration in the protein indicates that these vibrations are sensitive to environments. An action-absorption spectrum of cryogenically cold all-trans RPSB, reveals a∼ 310 ${{\rm{ \sim }}310}$  cm- 1 ${^{ - 1} }$ active vibration when using a hole-burning technique and 1500 cm- 1 ${^{ - 1} }$ C=C stretching modes.

Charge transfer characteristics in rhodopsin mimics during photoexcitation

To gain insights into the light-harvesting capabilities of the chromophores, it is essential to understand their molecular and electronic structures within their natural chemical or biological contexts. Rhodopsins display varied absorption characteristics due to the interaction between the chromophore retinal and its surrounding protein environments. In this study, we employed a quantum mechanics/molecular mechanics approach to examine a series of artificially designed rhodopsin mimics based on human cellular retinol acid binding protein 2 (hCRABP II). We elucidated the electron transfer within the all-trans protonated Schiff base upon light excitation, and our calculated absorption spectra show well consistency with the experimental result. Furthermore, the interaction mechanisms between the chromophore and the protein were investigated, and the relationship between the blueshifts and redshifts in the absorption spectra was analyzed. Our calculation demonstrates that the blueshifts and redshifts in the rhodopsin mimics correlate well with attractive (such as the hydrogen bonds or electrostatic interactions) and repulsive interactions (such as the steric effects) between the chromophore and the protein environment, respectively. These findings could provide hints for designing rhodopsin with absorption spectra at different wavelengths.

Low-temperature FTIR spectroscopy of the L/Q switch of proteorhodopsin

Rhodopsins are photoreceptive membrane proteins containing a retinal chromophore, and the color tuning mechanism in rhodopsins is one of the important topics. Color switch is a color-determining residue at the same position, where replacement of red- and blue-shifting amino acids in two wild-type rhodopsins causes spectral blue- and red-shifts, respectively. The first and most famous color switch in microbial rhodopsins is the L/Q switch in proteorhodopsins (PRs). Green- or blue-absorbing PR (GPR or BPR) contains Leu and Gln at position 105 of the C-helix (TM3), respectively, and their replacement converted absorbing colors. The L/Q switch enables bacteria to absorb green or blue light in shallow or deep ocean waters, respectively. Although Gln and Leu are hydrophilic and hydrophobic residues, respectively, a comprehensive mutation study of position 105 in GPR revealed that the λmax correlated with the volume of residues, not the hydropathy index. To gain structural insights into the mechanism, we applied low-temperature FTIR spectroscopy of L105Q GPR, and the obtained spectra were compared with those of GPR and BPR. The difference FTIR spectra of L105Q GPR were similar to those of BPR, not GPR, implying that the L/Q switch converts the GPR structure into a BPR structure in terms of the local environments of the retinal chromophore. It includes retinal skeletal vibration, hydrogen-bonding strength of the protonated Schiff base, amide-A vibration (peptide backbone), and protein-bound water molecules. Consequently color is switched accompanying such structural alterations, and known as the L/Q switch.

Photoisomerization pathway of the microbial rhodopsin chromophore in solution

Photoisomerization is a key photochemical reaction in microbial and animal rhodopsins. It is well established that such photoisomerization is highly selective; all-trans to 13-cis, and 11-cis to all-trans forms in microbial and animal rhodopsins, respectively. Nevertheless, unusual photoisomerization pathways have been discovered recently in microbial rhodopsins. In an enzymerhodopsin NeoR, the all-trans chromophore is isomerized into the 7-cis form exclusively, which is stable at room temperature. Although, the 7-cis form is produced by illumination of retinal, formation of the 7-cis form was never reported for a protonated Schiff base of all-trans retinal in solution. Present HPLC analysis of retinal oximes prepared by hydroxylamine reaction revealed that all-trans and 7-cis forms cannot be separated from the syn peaks under the standard HPLC conditions, while it is possible by the analysis of the anti-peaks. Consequently, we found formation of the 7-cis form by the photoreaction of all-trans chromophore in solution, regardless of the protonation state of the Schiff base. Upon light absorption of all-trans protonated retinal Schiff base in solution, excited-state relaxation accompanies double-bond isomerization, producing 7-cis, 9-cis, 11-cis, or 13-cis form. In contrast, specific chromophore-protein interaction enforces selective isomerization into the 13-cis form in many microbial rhodopsins, but into 7-cis in NeoR.

Ultrafast terahertz Stark spectroscopy reveals the excited-state dipole moments of retinal in bacteriorhodopsin

The photoinduced all-trans to 13-cis isomerization of the retinal Schiff base represents the ultrafast first step in the reaction cycle of bacteriorhodopsin (BR). Extensive experimental and theoretical work has addressed excited-state dynamics and isomerization via a conical intersection with the ground state. In conflicting molecular pictures, the excited state potential energy surface has been modeled as a pure S[Formula: see text] state that intersects with the ground state, or in a 3-state picture involving the S[Formula: see text] and S[Formula: see text] states. Here, the photoexcited system passes two crossing regions to return to the ground state. The electric dipole moment of the Schiff base in the S[Formula: see text] and S[Formula: see text] state differs strongly and, thus, its measurement allows for assessing the character of the excited-state potential. We apply the method of ultrafast terahertz (THz) Stark spectroscopy to measure electric dipole changes of wild-type BR and a BR D85T mutant upon electronic excitation. A fully reversible transient broadening and spectral shift of electronic absorption is induced by a picosecond THz field of several megavolts/cm and mapped by a 120-fs optical probe pulse. For both BR variants, we derive a moderate electric dipole change of 5 [Formula: see text] 1 Debye, which is markedly smaller than predicted for a neat S[Formula: see text]-character of the excited state. In contrast, S[Formula: see text]-admixture and temporal averaging of excited-state dynamics over the probe pulse duration gives a dipole change in line with experiment. Our results support a picture of electronic and nuclear dynamics governed by the interaction of S[Formula: see text] and S[Formula: see text] states in a 3-state model.

Mutational analyses identify a single amino acid critical for color tuning in proteorhodopsins

Microbial rhodopsins are light-activated proteins that contain seven transmembrane alpha-helices. Spectral tuning in microbial rhodopsins is a useful optogenetic tool. In this study, we report a new site that controls spectral tuning. In the proteorhodopsins ISR34 and ISR36, a single amino-acid substitution at Cys189 caused an absorption maximum shift of 44 nm, indicating spectral tuning at a specific site. Comparison of single amino acid substitutions was conducted using photochemical and photobiological approaches. The maximum absorption for red-shift was measured for mutations at positions 189 and 105 in ISR34, both residues being equally important. Structural changes resulting from amino acid substitutions are related to pK_a  values, pumping activity, and spectral tuning.

Far-Red Absorbing Rhodopsins, Insights From Heterodimeric Rhodopsin-Cyclases

The recently discovered Rhodopsin-cyclases from Chytridiomycota fungi show completely unexpected properties for microbial rhodopsins. These photoreceptors function exclusively as heterodimers, with the two subunits that have very different retinal chromophores. Among them is the bimodal photoswitchable Neorhodopsin (NeoR), which exhibits a near-infrared absorbing, highly fluorescent state. These are features that have never been described for any retinal photoreceptor. Here these properties are discussed in the context of color-tuning approaches of retinal chromophores, which have been extensively studied since the discovery of the first microbial rhodopsin, bacteriorhodopsin, in 1971 (Oesterhelt et al., Nature New Biology, 1971, 233 (39), 149-152). Further a brief review about the concept of heterodimerization is given, which is widely present in class III cyclases but is unknown for rhodopsins. NIR-sensitive retinal chromophores have greatly expanded our understanding of the spectral range of natural retinal photoreceptors and provide a novel perspective for the development of optogenetic tools.

Spectroscopy and photoisomerization of protonated Schiff-base retinal derivatives in vacuo

The protonated Schiff-base retinal acts as the chromophore in bacteriorhodopsin as well as in rhodopsin. In both cases, photoexcitation initializes fast isomerization which eventually results in storage of chemical energy or signaling. The details of the photophysics for this important chromophore is still not fully understood. In this study, action-absorption spectra and photoisomerization dynamics of three retinal derivatives are measured in the gas phase and compared to that of the protonated Schiff-base retinal. The retinal derivatives include C₉C₁₀trans-locked, C₁₃C₁₄trans-locked and a retinal derivative without the β-ionone ring. The spectroscopy as well as the isomerization speed of the chromophores are altered significantly as a consequence of the steric constraints.

Effect of Freezing out Vibrational Modes on Gas-Phase Fluorescence Spectra of Small Ionic Dyes

While action spectroscopy of cold molecular ions is a well-established technique to provide vibrationally resolved absorption features, fluorescence experiments are still challenging. Here we report the fluorescence spectra of pyronin-Y and resorufin ions at 100 K using a newly constructed setup. Spectra narrow upon cooling, and the emission maxima blueshift. Temperature effects are attributed to the population of vibrational excited levels in S₁, and that frequencies are lower in S₁ than in S₀. This picture is supported by calculated spectra based on a Franck-Condon model that not only predicts the observed change in maximum, but also assigns Franck-Condon active vibrations. In-plane vibrational modes that preserve the mirror plane present in both S₀ and S₁ of resorufin and pyronin Y account for most of the observed vibrational bands. Finally, at low temperatures, it is important to pick an excitation wavelength as far to the red as possible to not reheat the ions.

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.

DECAY DYNAMICS IN MOLECULAR BEAMS

The phenomenon of power law decays in molecular beams is reviewed. The transition from a canonical to a microcanonical description of the decay is analyzed, and the appearance of the power law decay derived. Deviations from a power law often contain information on parallel competing processes. This is illustrated with examples where thermal radiation or dark unimolecular channels are the competing processes. Also corrections to the power law due to finite heat capacities and from nonideal energy distributions are derived. Finally, the consequences for the interpretation of action spectroscopy data are reviewed. © 2020 The Authors. Mass Spectrometry Reviews published by Wiley Periodicals, Inc. Mass Spec Rev.

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

Benchmark and performance of long-range corrected time-dependent density functional tight binding (LC-TD-DFTB) on rhodopsins and light-harvesting complexes

The chromophores of rhodopsins (Rh) and light-harvesting (LH) complexes still represent a major challenge for a quantum chemical description due to their size and complex electronic structure. Since gradient corrected and hybrid density functional approaches have been shown to fail for these systems, only range-separated functionals seem to be a promising alternative to the more time consuming post-Hartree-Fock approaches. For extended sampling of optical properties, however, even more approximate approaches are required. Recently, a long-range corrected (LC) functional has been implemented into the efficient density functional tight binding (DFTB) method, allowing to sample the excited states properties of chromophores embedded into proteins using quantum mechanical/molecular mechanical (QM/MM) with the time-dependent (TD) DFTB approach. In the present study, we assess the accuracy of LC-TD-DFT and LC-TD-DFTB for rhodopsins (bacteriorhodopsin (bR) and pharaonis phoborhodopsin (ppR)) and LH complexes (light-harvesting complex II (LH2) and Fenna-Matthews-Olson (FMO) complex). This benchmark study shows the improved description of the color tuning parameters compared to standard DFT functionals. In general, LC-TD-DFTB can exhibit a similar performance as the corresponding LC functionals, allowing a reliable description of excited states properties at significantly reduced cost. The two chromophores investigated here pose complementary challenges: while huge sensitivity to external field perturbation (color tuning) and charge transfer excitations are characteristic for the retinal chromophore, the multi-chromophoric character of the LH complexes emphasizes a correct description of inter-chromophore couplings, giving less importance to color tuning. None of the investigated functionals masters both systems simultaneously with satisfactory accuracy. LC-TD-DFTB, at the current stage, although showing a systematic improvement compared to TD-DFTB cannot be recommended for studying color tuning in retinal proteins, similar to some of the LC-DFT functionals, because the response to external fields is still too weak. For sampling of LH-spectra, however, LC-TD-DFTB is a viable tool, allowing to efficiently sample absorption energies, as shown for three different LH complexes. As the calculations indicate, geometry optimization may overestimate the importance of local minima, which may be averaged over when using trajectories. Fast quantum chemical approaches therefore may allow for a direct sampling of spectra in the near future.

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.

Ultrafast photoisomerisation of an isolated retinoid

The photoinduced excited state dynamics of gas-phase trans-retinoate (deprotonated trans-retinoic acid, trans-RA^-) are studied using tandem ion mobility spectrometry coupled with laser spectroscopy, and frequency-, angle- and time-resolved photoelectron imaging. Photoexcitation of the bright S₃(ππ*) ← S₀ transition leads to internal conversion to the S₁(ππ*) state on a ≈80 fs timescale followed by recovery of S₀ and concomitant isomerisation to give the 13-cis (major) and 9-cis (minor) photoisomers on a ≈180 fs timescale. The sub-200 fs stereoselective photoisomerisation parallels that for the retinal protonated Schiff base chromophore in bacteriorhodopsin. Measurements on trans-RA^- in methanol using the solution photoisomerisation action spectroscopy technique show that 13-cis-RA^- is also the principal photoisomer, although the 13-cis and 9-cis photoisomers are formed with an inverted branching ratio with photon energy in methanol when compared with the gas phase, presumably due to solvent-induced modification of potential energy surfaces and inhibition of electron detachment processes. Comparison of the gas-phase time-resolved data with transient absorption spectroscopy measurements on retinoic acid in methanol suggest that photoisomerisation is roughly six times slower in solution. This work provides clear evidence that solvation significantly affects the photoisomerisation dynamics of retinoid molecules.

Red-shifting mutation of light-driven sodium-pump rhodopsin

Microbial rhodopsins are photoreceptive membrane proteins that transport various ions using light energy. While they are widely used in optogenetics to optically control neuronal activity, rhodopsins that function with longer-wavelength light are highly demanded because of their low phototoxicity and high tissue penetration. Here, we achieve a 40-nm red-shift in the absorption wavelength of a sodium-pump rhodopsin (KR2) by altering dipole moment of residues around the retinal chromophore (KR2 P219T/S254A) without impairing its ion-transport activity. Structural differences in the chromophore of the red-shifted protein from that of the wildtype are observed by Fourier transform infrared spectroscopy. QM/MM models generated with an automated protocol show that the changes in the electrostatic interaction between protein and chromophore induced by the amino-acid replacements, lowered the energy gap between the ground and the first electronically excited state. Based on these insights, a natural sodium pump with red-shifted absorption is identified from Jannaschia seosinensis.

Microbial rhodopsins are photoreceptive and widely used in optogenetics for which they should preferable function with longer-wavelength light. Here, authors achieve a 40-nm red-shift in the absorption wavelength of a sodium-pump rhodopsin (KR2) by altering the distribution of the retinal chromophore.

Intrinsic photoisomerization dynamics of protonated Schiff-base retinal

The retinal protonated Schiff-base (RPSB) in its all-trans form is found in bacterial rhodopsins, whereas visual rhodopsin proteins host 11-cis RPSB. In both cases, photoexcitation initiates fast isomerization of the retinal chromophore, leading to proton transport, storage of chemical energy or signaling. It is an unsolved problem, to which degree this is due to protein interactions or intrinsic RPSB quantum properties. Here, we report on time-resolved action-spectroscopy studies, which show, that upon photoexcitation, cis isomers of RPSB have an almost barrierless fast 400 fs decay, whereas all-trans isomers exhibit a barrier-controlled slow 3 ps decay. Moreover, formation of the 11-cis isomer is greatly favored for all-trans RPSB when isolated. The very fast photoresponse of visual photoreceptors is thus directly related to intrinsic retinal properties, whereas bacterial rhodopsins tune the excited state potential-energy surface to lower the barrier for particular double-bond isomerization, thus changing both the timescale and specificity of the photoisomerization.

The primary photoresponse of protonated Schiff-base retinal in visual and bacterial rhodopsins is fast sub-ps isomerisation. Here, the authors show that the fast photoisomerization of rhodopsin is related to an intrinsic retinal property, whereas bacterial rhodopsins tune the excited-state potential-energy surface and improve the isomerization timescale and specificity.

Action-spectroscopy studies of positively charge-tagged azobenzene in solution and in the gas-phase

The absorption of a positively charge-tagged azobenzene molecule is studied in the gas-phase by measuring photoinduced fragmentation of ions as a function of time. This technique provides information on prompt as well as delayed fragmentation, and a single dissociation channel after one-photon absorption is identified. The spectra in solution, as well as in the gas-phase, show a weak S₀ → S₁, a strong S₀ → S₂, and a broad absorption band in the UV regime. The bands are assigned through time dependent density functional theory calculations. The ratio of the various absorption bands depends on the trans to cis isomerization fraction and may be tuned by light irradiation. Gas-phase absorption spectra are presented and discussed in terms of trans and cis isomers.

Action and Ion Mobility Spectroscopy of a Shortened Retinal Derivative

The development of tandem ion mobility spectroscopy (IMS) known as IMS-IMS has led to extensive research into isomerizations of isolated molecules. Many recent works have focused on the retinal chromophore which is the optical switch used in animal vision. Here, we study a shortened derivative of the chromophore, which exhibits a rich IM spectrum allowing for a detailed analysis of its isomerization pathways, and show that the longer the chromophore is, the lower the barrier energies for isomerization are. Graphical Abstract.

Is the choice of a standard zeroth-order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal?

To account for systematic error of CASPT2 method empirical modification of the zeroth-order Hamiltonian with Ionization Potential-Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S-IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S-IPEA, is whether the introduction of the modified zeroth-order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0-0.35 a.u., on vertical excitation energies to low-lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S-IPEA significantly overestimates both S₀ →S₁ and S₀ →S₂ energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth-order Hamiltonian or small value of the IPEA shift in a range of 0.05-0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.

Shapes, Dynamics, and Stability of β-Ionone and Its Two Mutants Evidenced by High-Resolution Spectroscopy in the Gas Phase

The conformational landscapes of β-ionone and two mutants (α-ionone and β-damascone) have been analyzed by means of state-of-the-art rotational spectroscopy and quantum-chemical calculations. The experiments performed at high resolution and sensitivity have provided a deep insight into their conformational spaces, assigning more than 8000 transitions corresponding to the rotational structures of 54 different species (3 isomers, 14 conformers, and 40 isotopologues). Methyl internal rotation dynamics were also observed and analyzed. The work proved the great flexibility of β-ionone due to its flatter potential energy surface. This feature confers on β-ionone a wider ability to interconvert between conformers with rather similar energies with respect to its mutants, allowing the retinal ligand to better adapt inside the binding pocket.

Excitation spectra of retinal by multiconfiguration pair-density functional theory

Retinal is the chromophore in proteins responsible for vision. The absorption maximum of retinal is sensitive to mutations of the protein. However, it is not easy to predict the absorption spectrum of retinal accurately, and questions remain even after intensive investigation. Retinal poses a challenge for Kohn-Sham density functional theory (KS-DFT) because of the charge transfer character in its excitations, and it poses a challenge for wave function theory because the large size of the molecule makes multiconfigurational perturbation theory methods expensive. In this study, we demonstrate that multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to predict the vertical excitation energies of 11-Z retinal, and it reproduces the experimentally determined absorption band widths and peak positions better than complete active space second-order perturbation theory (CASPT2). The consistency between complete active space self-consistent field (CASSCF) and KS-DFT dipole moments is demonstrated to be a useful criterion in selecting the active space. We also found that the nature of the terminal groups and the conformations of retinal play a significant role in the absorption spectrum. By considering a thermal distribution of conformations, we predict an absorption spectrum of retinal that is consistent with the experimental gas-phase spectrum. The location of the absorption peak and the spectral broadening based on MC-PDFT calculations agree better with experiments than those of CASPT2.

The UV-visible action-absorption spectrum of all-trans and 11-cis protonated Schiff base retinal in the gas phase

The UV-visible absorption of retinal in its protonated Schiff-base form is studied in the gas phase. In particular, transitions to highly-excited electronic states, S_n, in the all-trans and 11-cis forms are considered, and several new states are discovered. Their positions and strengths are compared to state of the art quantum calculations. The location of these states are particularly important when new fs pump-probe experiments are designed to investigate the fast excited-state dynamics of retinal chromophores.

Elucidation of the intrinsic optical properties of hydrogen-bonded and protonated flavin chromophores by photodissociation action spectroscopy

The intrinsic optical properties of the flavin chromophore when engaged in hydrogen bonding or being protonated were elucidated by photo-induced action spectroscopy and computations.

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.

The femtosecond-to-second photochemistry of red-shifted fast-closing anion channelrhodopsin PsACR1

Femtosecond-to-second complete photocycle model of anion channelrhodopsin PsACR1.

Polarizable Embedding Density Matrix Renormalization Group

The polarizable embedding (PE) approach is a flexible embedding model where a preselected region out of a larger system is described quantum mechanically, while the interaction with the surrounding environment is modeled through an effective operator. This effective operator represents the environment by atom-centered multipoles and polarizabilities derived from quantum mechanical calculations on (fragments of) the environment. Thereby, the polarization of the environment is explicitly accounted for. Here, we present the coupling of the PE approach with the density matrix renormalization group (DMRG). This PE-DMRG method is particularly suitable for embedded subsystems that feature a dense manifold of frontier orbitals which requires large active spaces. Recovering such static electron-correlation effects in multiconfigurational electronic structure problems, while accounting for both electrostatics and polarization of a surrounding environment, allows us to describe strongly correlated electronic structures in complex molecular environments. We investigate various embedding potentials for the well-studied first excited state of water with active spaces that correspond to a full configuration-interaction treatment. Moreover, we study the environment effect on the first excited state of a retinylidene Schiff base within a channelrhodopsin protein. For this system, we also investigate the effect of dynamical correlation included through short-range density functional theory.

Theoretical Studies on the Color-Tuning Mechanism in Retinal Proteins

The excited states of the three retinal proteins, bovine rhodopsin (Rh), bacteriorhodopsin (bR), and sensory rhodopsin II (sRII) were studied using the symmetry-adapted cluster-configuration interaction (SAC-CI) and combined quantum mechanical and molecular mechanical (QM/MM) methods. The computed absorption energies are in good agreement with the experimental ones for all three proteins. The spectral tuning mechanism was analyzed in terms of three contributions:  molecular structures of the chromophore in the binding pockets, electrostatic (ES) interaction of the chromophore with the surrounding protein environment, and quantum-mechanical effect between the chromophore and the counterion group. This analysis provided an insight into the mechanism of the large blue-shifts in the absorption peak position of Rh and sRII from that of bR. Protein ES effect is primarily important both in Rh and in sRII, and the structure effect is secondary important in Rh. The quantum-mechanical interaction between the chromophore and the counterion is very important for quantitative reproduction of the excitation energy. These results indicate that the present approach is useful for studying the absorption spectra and the mechanism of the color tuning in the retinal proteins.

Molecular Mechanism of Wide Photoabsorption Spectral Shifts of Color Variants of Human Cellular Retinol Binding Protein II

Color variants of human cellular retinol binding protein II (hCRBPII) created by protein engineering were recently shown to exhibit anomalously wide photoabsorption spectral shifts over ∼200 nm across the visible region. The remarkable phenomenon provides a unique opportunity to gain insight into the molecular basis of the color tuning of retinal binding proteins for understanding of color vision as well as for engineering of novel color variants of retinal binding photoreceptor proteins employed in optogenetics. Here, we report a theoretical investigation of the molecular mechanism underlying the anomalously wide spectral shifts of the color variants of hCRBPII. Computational modeling of the color variants with hybrid molecular simulations of free energy geometry optimization succeeded in reproducing the experimentally observed wide spectral shifts, and revealed that protein flexibility, through which the active site structure of the protein and bound water molecules is altered by remote mutations, plays a significant role in inducing the large spectral shifts.

Photoisomerization action spectrum of retinal protonated Schiff base in the gas phase

The photophysical behaviour of the isolated retinal protonated n-butylamine Schiff base (RPSB) is investigated in the gas phase using a combination of ion mobility spectrometry and laser spectroscopy. The RPSB cations are introduced by electrospray ionisation into an ion mobility mass spectrometer where they are exposed to tunable laser radiation in the region of the S1 ← S0 transition (420-680 nm range). Four peaks are observed in the arrival time distribution of the RPSB ions. On the basis of predicted collision cross sections with nitrogen gas, the dominant peak is assigned to the all-trans isomer, whereas the subsidiary peaks are assigned to various single, double and triple cis geometric isomers. RPSB ions that absorb laser radiation undergo photoisomerization, leading to a detectable change in their drift speed. By monitoring the photoisomer signal as a function of laser wavelength an action spectrum, extending from 480 to 660 nm with a clear peak at 615 ± 5 nm, is obtained. The photoisomerization action spectrum is related to the absorption spectrum of isolated retinal RPSB molecules and should help benchmark future electronic structure calculations.

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.

Ion Mobility Unlocks the Photofragmentation Mechanism of Retinal Protonated Schiff Base

Retinal protonated Schiff base (RPSB) is a key molecular component of biological photoreceptors and bacterial photosynthetic structures, where its action involves photoisomerization around bonds in the polyene chain. In a vacuum environment, collisional activation or exposure to visible light causes the RPSB molecule to disintegrate, producing charged molecular fragments with m/z = 248 Da that cannot be formed by simple cleavage of the polyene chain. Photofragments resulting from laser excitation of RPSB at a wavelength of 532 nm are analyzed in an ion mobility mass spectrometer (IMMS) and found to be the protonated Schiff base of β-ionone. Density functional theory calculations at the M06-2X/cc-pVDZ level support a fragmentation mechanism in which RPSB undergoes an electrocyclization/fragmentation cascade with the production of protonated Schiff base of β-ionone and toluene.

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.