The light-induced double-bond isomerization of the visual pigment rhodopsin operates a molecular-level optomechanical energy transduction, which triggers a crucial protein structure change. In fact, rhodopsin isomerization occurs according to a unique, ultrafast mechanism that preserves mode-specific vibrational coherence all the way from the reactant excited state to the primary photoproduct ground state. The engineering of such an energy-funnelling function in synthetic compounds would pave the way towards biomimetic molecular machines capable of achieving optimum light-to-mechanical energy conversion. Here we use resonance and off-resonance vibrational coherence spectroscopy to demonstrate that a rhodopsin-like isomerization operates in a biomimetic molecular switch in solution. Furthermore, by using quantum chemical simulations, we show why the observed coherent nuclear motion critically depends on minor chemical modifications capable to induce specific geometric and electronic effects. This finding provides a strategy for engineering vibrationally coherent motions in other synthetic systems.

The activation of rhodopsin, the light-sensitive G-protein-coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. In this study we used hundreds of quantum-classical trajectories to show that, 15 fs after light absorption, a degeneracy between the reactive excited state and a neighbouring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show here that such splitting is modulated by protein electrostatics, thus linking amino acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.

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.

While the light-induced population dynamics of different photoresponsive proteins has been investigated spectroscopically, systematic computational studies have not yet been possible due to the phenomenally high cost of suitable high level quantum chemical methods and the need of propagating hundreds, if not thousands, of nonadiabatic trajectories. Here we explore the possibility of studying the photodynamics of rhodopsins by constructing and investigating quantum mechanics/molecular mechanics (QM/MM) models featuring reduced retinal chromophores. In order to do so we use the sensory rhodopsin found in the cyanobacterium Anabaena PCC7120 (ASR) as a benchmark system. We find that the basic mechanistic features associated with the excited state dynamics of ASR QM/MM models are reproduced using models incorporating a minimal (i.e., three double-bond) chromophore. Furthermore, we show that ensembles of nonadiabatic ASR trajectories computed using the same abridged models replicate, at both the CASPT2 and CASSCF levels of theory, the trends in spectroscopy and lifetimes estimated using unabridged models and observed experimentally at room temperature. We conclude that a further expansion of these studies may lead to low-cost QM/MM rhodopsin models that may be used as effective tools in high-throughput in silico mutant screening.

While rotary molecular switches based on neutral and cationic organic π-systems have been reported, structurally homologous anionic switches providing complementary properties have not been prepared so far. Here we report the design and preparation of a molecular switch mimicking the anionic p-HBDI chromophore of the green fluorescent protein. The investigation of the mechanism and dynamics of the E/Z switching function is carried out both computationally and experimentally. The data consistently support axial rotary motion occurring on a sub-picosecond time scale. Transient spectroscopy and trajectory simulations show that the nonadiabatic decay process occurs in the vicinity of a conical intersection (CInt) between a charge transfer state and a covalent/diradical state. Comparison of our anionic p-HBDI-like switch with the previously reported cationic N-alkyl indanylidene pyrrolinium switch mimicking visual pigments reveals that these similar systems translocate, upon vertical excitation, a similar net charge in the same axial direction.