Picosecond quantum-classical dynamics reveals that the coexistence of light-induced microbial and animal chromophore rotary motion modulates the isomerization quantum yield of heliorhodopsin.

Rhodopsins are light-responsive proteins forming two vast and evolutionary distinct superfamilies whose functions are invariably triggered by the photoisomerization of a single retinal chromophore. In 2018 a third widespread superfamily of rhodopsins called heliorhodopsins was discovered using functional metagenomics. Heliorhodopsins, with their markedly different structural features with respect to the animal and microbial superfamilies, offer an opportunity to study how evolution has manipulated the chromophore photoisomerization to achieve adaptation. One question is related to the mechanism of such a reaction and how it differs from that of animal and microbial rhodopsins. To address this question, we use hundreds of quantum-classical trajectories to simulate the spectroscopically documented picosecond light-induced dynamics of a heliorhodopsin from the archaea thermoplasmatales archaeon (TaHeR). We show that, consistently with the observations, the trajectories reveal two excited state decay channels. However, inconsistently with previous hypotheses, only one channel is associated with the -C13C14- rotation of microbial rhodopsins while the second channel is characterized by the -C11C12- rotation typical of animal rhodopsins. The fact that such -C11C12- rotation is aborted upon decay and ground state relaxation, explains why illumination of TaHeR only produces the 13-cis isomer with a low quantum efficiency. We argue that the documented lack of regioselectivity in double-bond excited state twisting motion is the result of an "adaptation" that could be completely lost via specific residue substitutions modulating the steric hindrance experienced along the isomerization motion.

Experimental Assessment of the Electronic and Geometrical Structure of a Near-Infrared Absorbing and Highly Fluorescent Microbial Rhodopsin.

The recently discovered Neorhodopsin (NeoR) exhibits absorption and emission maxima in the near-infrared spectral region, which together with the high fluorescence quantum yield makes it an attractive retinal protein for optogenetic applications. The unique optical properties can be rationalized by a theoretical model that predicts a high charge transfer character in the electronic ground state (S0) which is otherwise typical of the excited state S1 in canonical retinal proteins. The present study sets out to assess the electronic structure of the NeoR chromophore by resonance Raman (RR) spectroscopy since frequencies and relative intensities of RR bands are controlled by the ground and excited state's properties. The RR spectra of NeoR differ dramatically from those of canonical rhodopsins but can be reliably reproduced by the calculations carried out within two different structural models. The remarkable agreement between the experimental and calculated spectra confirms the consistency and robustness of the theoretical approach.

Retinal chromophore charge delocalization and confinement explain the extreme photophysics of Neorhodopsin.

The understanding of how the rhodopsin sequence can be modified to exactly modulate the spectroscopic properties of its retinal chromophore, is a prerequisite for the rational design of more effective optogenetic tools. One key problem is that of establishing the rules to be satisfied for achieving highly fluorescent rhodopsins with a near infrared absorption. In the present paper we use multi-configurational quantum chemistry to construct a computer model of a recently discovered natural rhodopsin, Neorhodopsin, displaying exactly such properties. We show that the model, that successfully replicates the relevant experimental observables, unveils a geometrical and electronic structure of the chromophore featuring a highly diffuse charge distribution along its conjugated chain. The same model reveals that a charge confinement process occurring along the chromophore excited state isomerization coordinate, is the primary cause of the observed fluorescence enhancement.