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 engineering of microbial rhodopsins with enhanced fluorescence is of great importance in the expanding field of optogenetics. Here we report the discovery of two mutants (W76S/Y179F and L83Q) of a sensory rhodopsin from the cyanobacterium Anabaena PCC7120 with opposite fluorescence behavior. In fact, while W76S/Y179F displays, with respect to the wild-type protein, a nearly tenfold increase in red-light emission, the second is not emissive. Thus, the W76S/ Y179F, L83Q pair offers an unprecedented opportunity for the investigation of fluorescence enhancement in microbial rhodopsins, which is pursued by combining transient absorption spectroscopy and multi-configurational quantum chemistry. The results of such an investigation point to an isomerization-blocking electronic effect as the direct cause of instantaneous (subpicosecond) fluorescence enhancement.

A novel atomistic methodology to perform free energy geometry optimization of a retinal chromophore covalently bound to any rhodopsin-like protein cavity is presented and benchmarked by computing the absorption maxima wavelengths (λ) of distant rhodopsin systems. The optimization is achieved by computing the Nagaoka's Free Energy Gradient (FEG) within an Average Solvent Electrostatic Configuration (ASEC) atomistic representation of the thermodynamic equilibrium and minimizing such quantity via an iterative procedure based on sequential classical MD and constrained QM/MM geometry optimization steps. The performance of such an ASEC-FEG protocol is assessed at the CASPT2//CASSCF/Amber level by reproducing the λ values observed for 12 mutants of redesigned human cellular retinol binding protein II (hCRBPII) systems; a set of 10 distant wild-type rhodopsins from vertebrates, invertebrates, eubacteria, and archaea organisms; and finally a set of 10 rhodopsin mutants from an eubacterial rhodopsin. The results clearly show that the proposed protocol, which can be easily extended to any protein incorporating a covalently bound ligand, yields correct λ trends with limited absolute errors.