In-silico predicted mouse melanopsins with blue spectral shifts deliver efficient subcellular signaling

Melanopsin is a photopigment belonging to the G Protein-Coupled Receptor (GPCR) family expressed in a subset of intrinsically photosensitive retinal ganglion cells (ipRGCs) and responsible for a variety of processes. The bistability and, thus, the possibility to function under low retinal availability would make melanopsin a powerful optogenetic tool. Here, we aim to utilize mouse melanopsin to trigger macrophage migration by its subcellular optical activation with localized blue light, while simultaneously imaging the migration with red light. To reduce melanopsin's red light sensitivity, we employ a combination of in silico structure prediction and automated quantum mechanics/molecular mechanics modeling to predict minimally invasive mutations to shift its absorption spectrum towards the shorter wavelength region of the visible spectrum without compromising the signaling efficiency. The results demonstrate that it is possible to achieve melanopsin mutants that resist red light-induced activation but are activated by blue light and display properties indicating preserved bistability. Using the A333T mutant, we show that the blue light-induced subcellular melanopsin activation triggers localized PIP3 generation and macrophage migration, which we imaged using red light, demonstrating the optogenetic utility of minimally engineered melanopsins.

Structural Models of the First Molecular Events in the Heliorhodopsin Photocycle

Retinylidene conformations and rearrangements of the hydrogen-bond network in the vicinity of the protonated Schiff base (PSB) play a key role in the proton transfer process in the Heliorhodopsin photocycle. Photoisomerization of the retinylidene chromophore and the formation of photoproducts corresponding to the early intermediates were modeled using a combination of molecular dynamics simulations and quantum mechanical/molecular mechanics calculations. The resulting structures were refined, and the respective excitation energies were calculated. Aided by metadynamics simulations, we constructed a photoisomerized intermediate where the 13-cis retinylidene chromophore is rotated about a parallel pair of double bonds at C13=C14 and C15=NZ double bonds. We demonstrate how the deprotonation of the Schiff base and the concomitant protonation of the Glu107 counterion are only favored because of these rearrangements.

A Portrait of the Chromophore as a Young System-Quantum-Derived Force Field Unraveling Solvent Reorganization upon Optical Excitation of Cyclocurcumin Derivatives

The study of fast non-equilibrium solvent relaxation in organic chromophores is still challenging for molecular modeling and simulation approaches, and is often overlooked, even in the case of non-adiabatic dynamics simulations. Yet, especially in the case of photoswitches, the interaction with the environment can strongly modulate the photophysical outcomes. To unravel such a delicate interplay, in the present contribution we resorted to a mixed quantum-classical approach, based on quantum mechanically derived force fields. The main task is to rationalize the solvent reorganization pathways in chromophores derived from cyclocurcumin, which are suitable for light-activated chemotherapy to destabilize cellular lipid membranes. The accurate and reliable decryption delivered by the quantum-derived force fields points to important differences in the solvent's reorganization, in terms of both structure and time scale evolution.

Convergent mechanism underlying the acquisition of vertebrate scotopic vision

High sensitivity of scotopic vision (vision in dim light conditions) is achieved by the rods' low background noise, which is attributed to a much lower thermal activation rate (kth) of rhodopsin compared with cone pigments. Frogs and nocturnal geckos uniquely possess atypical rods containing noncanonical cone pigments that exhibit low kth, mimicking rhodopsin. Here, we investigated the convergent mechanism underlying the low kth of rhodopsins and noncanonical cone pigments. Our biochemical analysis revealed that the kth of canonical cone pigments depends on their absorption maximum (λmax). However, rhodopsin and noncanonical cone pigments showed a substantially lower kth than predicted from the λmax dependency. Given that the λmax is inversely proportional to the activation energy of the pigments in the Hinshelwood distribution-based model, our findings suggest that rhodopsin and noncanonical cone pigments have convergently acquired low frequency of spontaneous-activation attempts, including thermal fluctuations of the protein moiety, in the molecular evolutionary processes from canonical cone pigments, which contributes to highly sensitive scotopic vision.

Quantum Mechanics/Molecular mechanics calculations predict A1, not A2, is present in melanopsin (Opn4m) of red-eared slider turtles (Trachemys scripta elegans)

Melanopsin is a photopigment that plays a role in non-visual, light-driven, cellular processes such as modulation of circadian rhythms, retinal vascular development, and the pupillary light reflex (PLR). In this study, computational methods were used to understand which chromophore is harbored by melanopsin in red-eared slider turtles (Trachemys scripta elegans). In mammals, the vitamin A derivative 11-cis-retinal (A1) is the chromophore, which provides functionality for melanopsin. However, in red-eared slider turtles, a member of the reptilian class, the identity of the chromophore remains unclear. Red-eared slider turtles, similar to other freshwater vertebrates, possess visual pigments that harbor a different vitamin A derivative, 11-cis-3,4-didehydroretinal (A2), making their pigments more sensitive to red-light than blue-light, therefore, suggesting the chromophore to be the A2 derivative instead of the A1. To help resolve the chromophore identity, in this work, computational homology models of melanopsin in red-eared slider turtles were first constructed. Next, quantum mechanics/molecular mechanics (QM/MM) calculations were carried out to compare how A1 and A2 derivatives bind to melanopsin. Time dependent density functional theory (TDDFT) calculations were then used to determine the excitation energy of the pigments. Lastly, calculated excitation energies were compared to experimental spectral sensitivity data from responses by the irises of red-eared sliders. Contrary to what was expected, our results suggest that melanopsin in red-eared slider turtles is more likely to harbor the A1 chromophore than the A2. Furthermore, a glutamine (Q622.56) and tyrosine (Y853.28) residue in the chromophore binding pocket are shown to play a role in the spectral tuning of the chromophore.

Extended Single-Electron Transfer Model and Dynamically Associated Energy Transfer Event in a Dual-Functional Catalyst System

Organic photocatalysis has been developed flourishingly to rely on bimolecular energy transfer (EnT) or oxidative/reductive electron transfer (ET), promoting a variety of synthetic transformations. However, there are rare examples to merge EnT and ET processes rationally within one chemical system, of which the mechanistic investigation still remains in its infancy. Herein, the first mechanistic illustration and kinetic assessments of the dynamically associated EnT and ET paths were conducted for realizing the C-H functionalization in a cascade photochemical transformation of isomerization and cyclization by using the dual-functional organic photocatalyst of riboflavin. An extended single-electron transfer model of transition-state-coupled dual-nonadiabatic crossings was explored to analyze the dynamic behaviors in the proton transfer-coupled cyclization. This can also be used to clarify the dynamic correlation with the EnT-driven E → Z photoisomerization that has been kinetically evaluated by using Fermi's golden rule with the Dexter model. The present computational results of electron structures and kinetic data contribute to a fundamental basis for understanding the photocatalytic mechanism of the combined operation of EnT and ET strategies, which will guide the design and manipulation for the implementation of multiple activation modes based on a single photosensitizer.

On the fluorescence enhancement of arch neuronal optogenetic reporters

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter.

Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency

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. Here, we use hundreds of quantum-classical trajectories to show that, 15 femtoseconds after light absorption, a degeneracy between the reactive excited state and a neighboring 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 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.

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.

Torsionally broken symmetry assists infrared excitation of biomimetic charge-coupled nuclear motions in the electronic ground state

The concerted interplay between reactive nuclear and electronic motions in molecules actuates chemistry. Here, we demonstrate that out-of-plane torsional deformation and vibrational excitation of stretching motions in the electronic ground...

Analytical Gradient Theory for Resolvent-Fitted Second-Order Extended Multiconfiguration Perturbation Theory (XMCQDPT2)

We present the formulation and implementation of an analytical gradient algorithm for extended multiconfiguration quasidegenerate perturbation theory (XMCQDPT2) with the resolvent-fitting approximation by Granovsky. This algorithm is powerful when optimizing molecular configurations with a moderate-sized active space and many electronic states. First, we present the powerfulness and accuracy of resolvent-fitting approximations compared to canonical XMCQDPT2 theory. Then, we demonstrate the utility of the current algorithm in frequency analyses, optimizing the minimum energy conical intersection geometries of the retinal chromophore model RPSB6 and evaluating nuclear gradients when there are many electronic states. Furthermore, we parallelize the algorithm using the OpenMP/MPI hybrid approach. Additionally, we report the computational cost and parallel efficiency of the program.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Free Energy Computation for an Isomerizing Chromophore in a Molecular Cavity via the Average Solvent Electrostatic Configuration Model: Applications in Rhodopsin and Rhodopsin-Mimicking Systems

We present a novel technique for computing the free energy differences between two chromophore "isomers" hosted in a molecular environment (a generalized solvent). Such an environment may range from a relatively rigid protein cavity to a flexible solvent environment. The technique is characterized by the application of the previously reported "average electrostatic solvent configuration" method, and it is based on the idea of using the free energy perturbation theory along with a chromophore annihilation procedure in thermodynamic cycle calculations. The method is benchmarked by computing the ground-state room-temperature relative stabilities between (i) the cis and trans isomers of prototypal animal and microbial rhodopsins and (ii) the analogue isomers of a rhodopsin-like light-driven molecular switch in methanol. Furthermore, we show that the same technology can be used to estimate the activation free energy for the thermal isomerization of systems i-ii by replacing one isomer with a transition state. The results show that the computed relative stability and isomerization barrier magnitudes for the selected systems are in line with the available experimental observation in spite of their widely diverse complexity.

Synthesis and Photoswitching Properties of Bioinspired Dissymmetric γ-Pyrone, an Analogue of Cyclocurcumin

Cyclocurcumin (CC), a turmeric curcuminoid with potential therapeutic properties, is also a natural photoswitch that may undergo E/Z photoisomerization under UV light. To be further exploited in relevant biological applications, photoactivation under near-infrared (NIR) irradiation is required. Such requirement can be met through opportune chemical modifications, by favoring two-photon absorption (TPA) probability. Herein, a general and efficient synthesis of a biomimetic 2,6-disubstituted-γ-pyrone analogue of CC is described, motivated by the fact that molecular modeling previews an order of magnitude increase of its NIR TPA compared to CC. Three retrosynthetic pathways have been identified (i) via an aryl-oxazole intermediate or via aryl-diynone through (ii) a bottom-up or (iii) a top-down approach. While avoiding the passage through unstable synthons or low-yield intermediate reactions, only the latest approach could conveniently afford the 2,6-disubstituted-γ-pyrone analogue of CC, in ten steps and with an overall yield of 18%. The photophysical properties of our biomimetic analogue have also been characterized showing an improved photoisomerization yield over the parent natural compound. The potentially improved nonlinear optical properties, as well as enhanced stability, may be correlated to the enforcement of the planarity of the pyrone moiety leading to a quadrupolar D-π-A-π-D system.

Light-induced difference FTIR spectroscopy of primate blue-sensitive visual pigment at 163 K

Structural studies of color visual pigments lag far behind those of rhodopsin for scotopic vision. Using difference FTIR spectroscopy at 77 K, we report the first structural data of three primate color visual pig­ments, monkey red (MR), green (MG), and blue (MB), where the batho-intermediate (Batho) exhibits photo­equilibrium with the unphotolyzed state. This photo­chromic property is highly advantageous for limited samples since the signal-to-noise ratio is improved, but may not be applicable to late intermediates, because of large structural changes to proteins. Here we report the photochromic property of MB at 163 K, where the BL intermediate, formed by the relaxation of Batho, is in photoequilibrium with the initial MB state. A comparison of the difference FTIR spectra at 77 and 163 K provided information on what happens in the process of transition from Batho to BL in MB. The coupled C₁₁=C₁₂ HOOP vibration in the planer structure in MB is decoupled by distortion in Batho after retinal photoisomerization, but returns to the coupled C₁₁=C₁₂ HOOP vibration in the all-trans chromophore in BL. The Batho formation accompanies helical structural perturbation, which is relaxed in BL. Protein-bound water molecules that form an extended water cluster near the retinal chromophore change hydrogen bonds differently for Batho and BL, being stronger in the latter than in the initial state. In addition to structural dynamics, the present FTIR spectra show no signals of protonated carboxylic acids at 77 and 163 K, suggesting that E181 is deprotonated in MB, Batho and BL.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

QM/MM Investigation of the Spectroscopic Properties of the Fluorophore of Bacterial Luciferase

We employ replica-exchange molecular dynamics (REMD) and a hybrid ab initio multiconfigurational quantum mechanics/molecular mechanics (QM/MM) approach to model the absorption and fluorescence properties of bacterial luciferin-luciferase. Specifically, we employ complete active space perturbation theory (CASPT2) and study the effect of active space, basis set, and IPEA shift on the computed energies. We discuss the effect of the protein environment on the fluorophore's excited-state potential energy surface and the role that the protein plays in enhancing the fluorescence quantum yield in bacterial bioluminescence.

Accelerated evolution and positive selection of rhodopsin in Tibetan loaches living in high altitude

Rhodopsin (RH1), the temperature-sensitive visual pigment, attained cold adaptation by functional trade-offs between protein stability and activity. Recent studies suggested convergent selection pressures drove cold adaptation of rhodopsin in high altitude catfishes through nonparallel molecular mechanisms. Here, we tested whether the similar shift occurred in RH1 of Tibetan loaches on the Qinghai-Tibet Plateau (QTP) by investigating the molecular evolution and potential effect on function of RH1. We sequenced RH1 from 27 Triplophysa species, and four lowland loaches and combined these data with published sequences. Tests using a series of models of molecular evolution resulted in strong evidence for accelerated evolution and positive selection in Triplophysa RH1. Three positively selected sites were near key functional domains modulating nonspectral properties of rhodopsin, substitutions of which were likely to compensate for cold-induced decrease in rhodopsin kinetics in cold environments. Moreover, although accelerated evolutionary rates in Tibetan loaches was convergent with those in high altitude catfishes, the sites under positive selection were nonoverlapping. Our findings provide evidence for convergent shift in selection pressures of RH1 in high altitude fish during the ecological transition to cold environment of the QTP.

Electrostatic control of photoisomerization pathways in proteins

Rotation around a specific bond after photoexcitation is central to vision and emerging opportunities in optogenetics, super-resolution microscopy, and photoactive molecular devices. Competing roles for steric and electrostatic effects that govern bond-specific photoisomerization have been widely discussed, the latter originating from chromophore charge transfer upon excitation. We systematically altered the electrostatic properties of the green fluorescent protein chromophore in a photoswitchable variant, Dronpa2, using amber suppression to introduce electron-donating and electron-withdrawing groups to the phenolate ring. Through analysis of the absorption (color), fluorescence quantum yield, and energy barriers to ground- and excited-state isomerization, we evaluate the contributions of sterics and electrostatics quantitatively and demonstrate how electrostatic effects bias the pathway of chromophore photoisomerization, leading to a generalized framework to guide protein design.

Web-ARM: A Web-Based Interface for the Automatic Construction of QM/MM Models of Rhodopsins

This article introduces Web-ARM, a specialized tool, online available, designed to build quantum mechanical/molecular mechanical models of rhodopsins, a widely spread family of light-responsive proteins. Web-ARM allows the rapidly building of models of rhodopsins with a documented quality and the prediction of trends in UV-vis absorption maximum wavelengths, based on their excitation energies computed at the CASPT2//CASSCF/Amber level of theory. Web-ARM builds upon the recently reported, python-based a-ARM protocol [J. Chem. Theory Comput., 2019, 15, 3134-3152] and, as such, necessitates only a crystallographic structure or a comparative model in PDB format and a very basic knowledge of the studied rhodopsin system. The user-friendly web interface uses such input to generate congruous, gas-phase models of rhodopsins and, if requested, their mutants. We present two possible applications of Web-ARM, which showcase how the interface can be employed to assist both research and educational activities in fields at the interface between chemistry and biology. The first application shows how, through Web-ARM, research projects (e.g., rhodopsin and rhodopsin mutant screening) can be carried out in significantly less time with respect to using the required computational photochemistry tools via a command line. The second application documents the use of Web-ARM in a real-life educational/training activity, through a hands-on experience illustrating the concepts of rhodopsin color tuning.

Four resonance structures elucidate double-bond isomerisation of a biological chromophore

Four resonance structures determining the electronic structure of the chromophore’s ground and first excited states. Changing the relative energies of the structures by hydrogen-bonding interactions tunes all chromophore’s photochemical properties.

Molecular mechanism for thermal denaturation of thermophilic rhodopsin

Understanding the factors affecting the stability and function of proteins at the molecular level is of fundamental importance. In spite of their use in bioelectronics and optogenetics, factors influencing thermal stability of microbial rhodopsins, a class of photoreceptor protein ubiquitous in nature are not yet well-understood. Here we report on the molecular mechanism for thermal denaturation of microbial retinal proteins, including, a highly thermostable protein, thermophilic rhodopsin (TR). External stimuli-dependent thermal denaturation of TR, the proton pumping rhodopsin of Thermus thermophilus bacterium, and other microbial rhodopsins are spectroscopically studied to decipher the common factors guiding their thermal stability. The thermal denaturation process of the studied proteins is light-catalyzed and the apo-protein is thermally less stable than the corresponding retinal-covalently bound opsin. In addition, changes in structure of the retinal chromophore affect the thermal stability of TR. Our results indicate that the hydrolysis of the retinal protonated Schiff base (PSB) is the rate-determining step for denaturation of the TR as well as other retinal proteins. Unusually high thermal stability of TR multilayers, in which PSB hydrolysis is restricted due to lack of bulk water, strongly supports this proposal. Our results also show that the protonation state of the PSB counter-ion does not affect the thermal stability of the studied proteins. Thermal photo-bleaching of an artificial TR pigment derived from non-isomerizable trans-locked retinal suggests, rather counterintuitively, that the photoinduced retinal trans-cis isomerization is not a pre-requisite for light catalyzed thermal denaturation of TR. Protein conformation alteration triggered by light-induced retinal excited state formation is likely to facilitate the PSB hydrolysis.

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.

Rejection of the biophoton hypothesis on the origin of photoreceptor dark noise

It has been suggested that retinal “dark light” is caused by photons emitted by the retina itself. The authors show that the “biophoton” radiation from the retina can be detected, but its intensity is ≥100 times lower than necessary to produce the measured physiological noise understood to arise from the spontaneous activation of rhodopsin.

Rod photoreceptors of the vertebrate retina produce, in darkness, spontaneous discrete current waves virtually identical to responses to single photons. The waves comprise an irreducible source of noise (discrete dark noise) that may limit the threshold sensitivity of vision. The waves obviously originate from acts of random activation of single rhodopsin molecules. Until recently, it was generally accepted that the activation occurs due to the rhodopsin thermal motion. Yet, a few years ago it was proposed that rhodopsin molecules are activated not by heat but rather by real photons generated within the retina by chemiluminescence. Using a high-sensitive photomultiplier, we measured intensities of biophoton emission from isolated retinas and eyecups of frogs (Rana ridibunda) and fish (sterlet, Acipenser ruthenus). Retinal samples were placed in a perfusion chamber and emitted photons collected by a high-aperture quartz lens. The collected light was sent to the photomultiplier cathode through a rotating chopper so that a long-lasting synchronous accumulation of the light signal was possible. The absolute intensity of bio-emission was estimated by the response of the measuring system to a calibrated light source. The intensity of the source, in turn, was quantified by measuring rhodopsin bleaching with single-rod microspectrophotometry. We also measured the frequency of discrete dark waves in rods of the two species with suction pipette recordings. Expressed as the rate constant of rhodopsin activation, it was 1.2 × 10⁻¹¹/s in frogs and 7.6 × 10⁻¹¹/s in sterlets. Approximately two thirds of retinal samples of each species produced reliably measurable biophoton emissions. However, its intensity was ≥100 times lower than necessary to produce the discrete dark noise. We argue that this is just a lower estimate of the discrepancy between the hypothesis and experiment. We conclude that the biophoton hypothesis on the origin of discrete dark noise in photoreceptors must be rejected.

Shifts in Selective Pressures on Snake Phototransduction Genes Associated with Photoreceptor Transmutation and Dim-Light Ancestry

The visual systems of snakes are heavily modified relative to other squamates, a condition often thought to reflect their fossorial origins. Further modifications are seen in caenophidian snakes, where evolutionary transitions between rod and cone photoreceptors, termed photoreceptor transmutations, have occurred in many lineages. Little previous work, however, has focused on the molecular evolutionary underpinnings of these morphological changes. To address this, we sequenced seven snake eye transcriptomes and utilized new whole-genome and targeted capture sequencing data. We used these data to analyze gene loss and shifts in selection pressures in phototransduction genes that may be associated with snake evolutionary origins and photoreceptor transmutation. We identified the surprising loss of rhodopsin kinase (GRK1), despite a low degree of gene loss overall and a lack of relaxed selection early during snake evolution. These results provide some of the first evolutionary genomic corroboration for a dim-light ancestor that lacks strong fossorial adaptations. Our results also indicate that snakes with photoreceptor transmutation experienced significantly different selection pressures from other reptiles. Significant positive selection was found primarily in cone-specific genes, but not rod-specific genes, contrary to our expectations. These results reveal potential molecular adaptations associated with photoreceptor transmutation and also highlight unappreciated functional differences between rod- and cone-specific phototransduction proteins. This intriguing example of snake visual system evolution illustrates how the underlying molecular components of a complex system can be reshaped in response to changing selection pressures.

Temporal resolution of single-photon responses in primate rod photoreceptors and limits imposed by cellular noise

Sensory receptor noise corrupts sensory signals, contributing to imperfect perception and dictating central processing strategies. For example, noise in rod phototransduction limits our ability to detect light, and minimizing the impact of this noise requires precisely tuned nonlinear processing by the retina. But detection sensitivity is only one aspect of night vision: prompt and accurate behavior also requires that rods reliably encode the timing of photon arrivals. We show here that the temporal resolution of responses of primate rods is much finer than the duration of the light response and identify the key limiting sources of transduction noise. We also find that the thermal activation rate of rhodopsin is lower than previous estimates, implying that other noise sources are more important than previously appreciated. A model of rod single-photon responses reveals that the limiting noise relevant for behavior depends critically on how rod signals are pooled by downstream neurons. NEW & NOTEWORTHY Many studies have focused on the visual system's ability to detect photons, but not on its ability to encode the relative timing of detected photons. Timing is essential for computations such as determining the velocity of moving objects. Here we examine the timing precision of primate rod photoreceptor responses and show that it is more precise than previously appreciated. This motivates an evaluation of whether scotopic vision approaches limits imposed by rod temporal resolution.

Functional trade-offs and environmental variation shaped ancient trajectories in the evolution of dim-light vision

Trade-offs between protein stability and activity can restrict access to evolutionary trajectories, but widespread epistasis may facilitate indirect routes to adaptation. This may be enhanced by natural environmental variation, but in multicellular organisms this process is poorly understood. We investigated a paradoxical trajectory taken during the evolution of tetrapod dim-light vision, where in the rod visual pigment rhodopsin, E122 was fixed 350 million years ago, a residue associated with increased active-state (MII) stability but greatly diminished rod photosensitivity. Here, we demonstrate that high MII stability could have likely evolved without E122, but instead, selection appears to have entrenched E122 in tetrapods via epistatic interactions with nearby coevolving sites. In fishes by contrast, selection may have exploited these epistatic effects to explore alternative trajectories, but via indirect routes with low MII stability. Our results suggest that within tetrapods, E122 and high MII stability cannot be sacrificed-not even for improvements to rod photosensitivity.

Fine Tuning of Retinal Photoinduced Decay in Solution

Single methylation at position C₁₀ of the all-trans retinal protonated Schiff base switches its excited-state decay in methanol from a slower picosecond into an ultrafast, protein-like subpicosecond process. QM/MM modeling in conjunction with on-the-fly excited-state dynamics provides fundamental understanding of the fine-tuning mechanics that "catalyzes" the photoinduced decay of solvated retinals. Methylation alters the interplay between the ionic S₁ and covalent S₂ states, reducing the excited-state lifetime by favoring the formation of a S₁ transient fluorescent state with fully inverted bond lengths that accounts for the recorded transient spectroscopy and from which a space-saving conical intersection seam is quickly (<1 ps) reached. Minimal and apparently innocent chemical modifications thus affect the characteristic intramolecular charge-transfer of the S₁ state as well as the interaction with the covalent S₂ excited state, eventually providing the high tunability of retinal photophysics and photochemistry and delivering a new concept for the rational design of retinal-based photoactive molecular devices.

Evolution of nonspectral rhodopsin function at high altitudes

High-altitude environments present a range of biochemical and physiological challenges for organisms through decreases in oxygen, pressure, and temperature relative to lowland habitats. Protein-level adaptations to hypoxic high-altitude conditions have been identified in multiple terrestrial endotherms; however, comparable adaptations in aquatic ectotherms, such as fishes, have not been as extensively characterized. In enzyme proteins, cold adaptation is attained through functional trade-offs between stability and activity, often mediated by substitutions outside the active site. Little is known whether signaling proteins [e.g., G protein-coupled receptors (GPCRs)] exhibit natural variation in response to cold temperatures. Rhodopsin (RH1), the temperature-sensitive visual pigment mediating dim-light vision, offers an opportunity to enhance our understanding of thermal adaptation in a model GPCR. Here, we investigate the evolution of rhodopsin function in an Andean mountain catfish system spanning a range of elevations. Using molecular evolutionary analyses and site-directed mutagenesis experiments, we provide evidence for cold adaptation in RH1. We find that unique amino acid substitutions occur at sites under positive selection in high-altitude catfishes, located at opposite ends of the RH1 intramolecular hydrogen-bonding network. Natural high-altitude variants introduced into these sites via mutagenesis have limited effects on spectral tuning, yet decrease the stability of dark-state and light-activated rhodopsin, accelerating the decay of ligand-bound forms. As found in cold-adapted enzymes, this phenotype likely compensates for a cold-induced decrease in kinetic rates-properties of rhodopsin that mediate rod sensitivity and visual performance. Our results support a role for natural variation in enhancing the performance of GPCRs in response to cold temperatures.

An experimental comparison of human and bovine rhodopsin provides insight into the molecular basis of retinal disease

Rhodopsin is the visual pigment that mediates dim-light vision in vertebrates and is a model system for the study of retinal disease. The majority of rhodopsin experiments are performed using bovine rhodopsin; however, recent evidence suggests that significant functional differences exist among mammalian rhodopsins. In this study, we identify differences in both thermal decay and light-activated retinal release rates between bovine and human rhodopsin and perform mutagenesis studies to highlight two clusters of substitutions that contribute to these differences. We also demonstrate that the retinitis pigmentosa-associated mutation G51A behaves differently in human rhodopsin compared to bovine rhodopsin and determine that the thermal decay rate of an ancestrally reconstructed mammalian rhodopsin displays an intermediate phenotype compared to the two extant pigments.

Mechanism and bottlenecks in strand photodissociation of split green fluorescent proteins (GFPs)

Split GFPs have been widely applied for monitoring protein-protein interactions by expressing GFPs as two or more constituent parts linked to separate proteins that only fluoresce on complementing with one another. Although this complementation is typically irreversible, it has been shown previously that light accelerates dissociation of a noncovalently attached β-strand from a circularly permuted split GFP, allowing the interaction to be reversible. Reversible complementation is desirable, but photodissociation has too low of an efficiency (quantum yield <1%) to be useful as an optogenetic tool. Understanding the physical origins of this low efficiency can provide strategies to improve it. We elucidated the mechanism of strand photodissociation by measuring the dependence of its rate on light intensity and point mutations. The results show that strand photodissociation is a two-step process involving light-activated cis-trans isomerization of the chromophore followed by light-independent strand dissociation. The dependence of the rate on temperature was then used to establish a potential energy surface (PES) diagram along the photodissociation reaction coordinate. The resulting energetics-function model reveals the rate-limiting process to be the transition from the electronic excited-state to the ground-state PES accompanying cis-trans isomerization. Comparisons between split GFPs and other photosensory proteins, like photoactive yellow protein and rhodopsin, provide potential strategies for improving the photodissociation quantum yield.

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.

Measurements of the stabilities of isolated retinal chromophores

The barrier energies for isomerization and fragmentation were measured for a series of retinal chromophore derivatives using a tandem ion mobility spectrometry approach. These measurements allow us to quantify the effect of charge delocalization on the rigidity of chromophores. We find that the role of the methyl group on the C13 position is pivotal regarding the ground state dynamics of the chromophore. Additionally, a correlation between quasi-equilibrium isomer distribution and fragmentation pathways is observed.