Rhodopsin-based light-harvesting system for sustainable synthetic biology.

Rhodopsins, a diverse class of light-sensitive proteins found in various life domains, have attracted considerable interest for their potential applications in sustainable synthetic biology. These proteins exhibit remarkable photochemical properties, undergoing conformational changes upon light absorption that drive a variety of biological processes. Exploiting rhodopsin's natural properties could pave the way for creating sustainable and energy-efficient technologies. Rhodopsin-based light-harvesting systems offer innovative solutions to a few key challenges in sustainable engineering, from bioproduction to renewable energy conversion. In this opinion article, we explore the recent advancements and future possibilities of employing rhodopsins for sustainable engineering, underscoring the transformative potential of these biomolecules.

RhoMax: Computational Prediction of Rhodopsin Absorption Maxima Using Geometric Deep Learning.

Microbial rhodopsins (MRs) are a diverse and abundant family of photoactive membrane proteins that serve as model systems for biophysical techniques. Optogenetics utilizes genetic engineering to insert specialized proteins into specific neurons or brain regions, allowing for manipulation of their activity through light and enabling the mapping and control of specific brain areas in living organisms. The obstacle of optogenetics lies in the fact that light has a limited ability to penetrate biological tissues, particularly blue light in the visible spectrum. Despite this challenge, most optogenetic systems rely on blue light due to the scarcity of red-shifted opsins. Finding additional red-shifted rhodopsins would represent a major breakthrough in overcoming the challenge of limited light penetration in optogenetics. However, determining the wavelength absorption maxima for rhodopsins based on their protein sequence is a significant hurdle. Current experimental methods are time-consuming, while computational methods lack accuracy. The paper introduces a new computational approach called RhoMax that utilizes structure-based geometric deep learning to predict the absorption wavelength of rhodopsins solely based on their sequences. The method takes advantage of AlphaFold2 for accurate modeling of rhodopsin structures. Once trained on a balanced train set, RhoMax rapidly and precisely predicted the maximum absorption wavelength of more than half of the sequences in our test set with an accuracy of 0.03 eV. By leveraging computational methods for absorption maxima determination, we can drastically reduce the time needed for designing new red-shifted microbial rhodopsins, thereby facilitating advances in the field of optogenetics.

Parallel photocycle kinetic model of anion channelrhodopsin GtACR1 function.

The light-gated anion channelrhodopsin GtACR1 is an important optogenetic tool for neuronal silencing. Its photochemistry, including its photointermediates, is poorly understood. The current mechanistic view presumes BR-like kinetics and assigns the open channel to a blue-absorbing L intermediate. Based on time-resolved absorption and electrophysiological data, we recently proposed a red-absorbing spectral form for the open channel state. Here, we report the results of a comprehensive kinetic analysis of the spectroscopic data combined with channel current information. The time evolutions of the spectral forms derived from the spectroscopic data are inconsistent with the single chain mechanism and are analyzed within the concept of parallel photocycles. The spectral forms partitioned into conductive and nonconductive parallel cycles are assigned to intermediate states. Rejecting reversible connections between conductive and nonconductive channel states leads to kinetic schemes with two independent conductive states corresponding to the fast- and slow-decaying current components. The conductive cycle is discussed in terms of a single cycle and two parallel cycles. The reaction mechanisms and reaction rates for the wild-type protein, the A75E, and the low-conductance D234N and S97E protein variants are derived. The parallel cycles of channelrhodopsin kinetics, its relation to BR photocycle, and the role of the M intermediate in channel closure are discussed.

A subgroup of light-driven sodium pumps with an additional Schiff base counterion.

Light-driven sodium pumps (NaRs) are unique ion-transporting microbial rhodopsins. The major group of NaRs is characterized by an NDQ motif and has two aspartic acid residues in the central region essential for sodium transport. Here we identify a subgroup of the NDQ rhodopsins bearing an additional glutamic acid residue in the close vicinity to the retinal Schiff base. We thoroughly characterize a member of this subgroup, namely the protein ErNaR from Erythrobacter sp. HL-111 and show that the additional glutamic acid results in almost complete loss of pH sensitivity for sodium-pumping activity, which is in contrast to previously studied NaRs. ErNaR is capable of transporting sodium efficiently even at acidic pH levels. X-ray crystallography and single particle cryo-electron microscopy reveal that the additional glutamic acid residue mediates the connection between the other two Schiff base counterions and strongly interacts with the aspartic acid of the characteristic NDQ motif. Hence, it reduces its pKa. Our findings shed light on a subgroup of NaRs and might serve as a basis for their rational optimization for optogenetics.

Glu1022.53-Mediated Early Conformational Changes in the Process of Light-Induced Green Cone Pigment Activation.

Ligand-triggered activation of G protein-coupled receptors (GPCRs) relies on the phenomenon of loose allosteric coupling, which involves conformational alterations spanning from the extracellular ligand-binding domain to the cytoplasmic region, where interactions with G proteins occur. During the GPCR activation process, several intermediate and equilibrium states orchestrate the movement of the flexible and rigid transmembrane (TM) segments of the GPCR. Monitoring early conformational changes is important in unraveling the structural intricacies of the loose allosteric coupling. Here, we focus on the lumi intermediate formed by thermal relaxation from the initial photointermediate, batho in primate green cone pigment (MG), a light-sensitive GPCR responsible for color vision. Our findings from light-induced Fourier transform infrared difference spectroscopy reveal its similarity with rhodopsin, which mediates twilight vision, specifically involving the flip motion of the ß-ionone ring, the relaxation of the torsional structure of the retinal, and local perturbations in the α-helix upon lumi intermediate formation. Conversely, we observe a hydrogen bond modification specific to MG's protonated carboxylic acid, identifying its origin as Glu1022.53 situated in TM2. The weakening of the hydrogen bond strength at Glu1022.53 during the transition from the batho to the lumi intermediates corresponds to a slight outward movement of TM2. Additionally, within the X-ray crystal structure of the rhodopsin lumi intermediate, we note the relocation of the Met862.53 side chain in TM2, expanding the volume of the retinal binding pocket. Consequently, the position of 2.53 emerges as the early step in the conformational shift toward light-induced activation. Moreover, given the prevalence of IR-insensitive hydrophobic amino acids at position 2.53 in many rhodopsin-like GPCRs, including rhodopsin, the hydrogen bond alteration in the C═O stretching band at Glu1022.53 of MG can be used as a probe for tracing conformational changes during the GPCR activation process.

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.

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.

Color Tuning in Bovine Rhodopsin through Polarizable Embedding.

Bovine rhodopsin is among the most studied proteins in the rhodopsin family. Its primary activation mechanism is the photoisomerization of 11-cis retinal, triggered by the absorption of a UV-visible photon. Different mutants of the same rhodopsin show different absorption wavelengths due to the influence of the specific amino acid residues forming the cavity in which the retinal chromophore is embedded, and rhodopsins activated at different wavelengths are, for example, exploited in the field of optogenetics. In this letter, we present a procedure for systematically investigating color tuning in models of bovine rhodopsin and a set of its mutants embedded in a membrane bilayer. Vertical excitation energy calculations were carried out with the polarizable embedding potential for describing the environment surrounding the chromophore. We show that polarizable embedding outperformed regular electrostatic embedding in determining both the vertical excitation energies and associated oscillator strengths of the systems studied.

Configurational Changes of Retinal Schiff Base during Membrane Na+ Transport by a Sodium Pumping Rhodopsin.

Microbial rhodopsins are photoreceptors containing the retinal Schiff base chromophore and are ubiquitous among microorganisms. The Schiff base configuration of the chromophore, 15-anti (C═N trans) or 15-syn (C═N cis), is structurally important for their functions, such as membrane ion transport, because this configuration dictates the orientation of the positively charged NH group that interacts with substrate ions. The 15-anti/syn configuration is thus essential for elucidating the ion-transport mechanisms in microbial rhodopsins. Here, we identified the Schiff base configuration during the photoreaction of a sodium pumping rhodopsin from Indibacter alkaliphilus using Raman spectroscopy. We found that the unique configurational change from the 13-cis, 15-anti to all-trans, 15-syn form occurs between the photointermediates termed O1 and O2, which accomplish the Na+ uptake and release, respectively. This isomerization is considered to give rise to the highly irreversible O1 → O2 step that is crucial for unidirectional Na+ transport.

Scaling up Functional Analyses of the G Protein-Coupled Receptor Rhodopsin.

Eukaryotic cells use G protein-coupled receptors (GPCRs) to convert external stimuli into internal signals to elicit cellular responses. However, how mutations in GPCR-coding genes affect GPCR activation and downstream signaling pathways remain poorly understood. Approaches such as deep mutational scanning show promise in investigations of GPCRs, but a high-throughput method to measure rhodopsin activation has yet to be achieved. Here, we scale up a fluorescent reporter assay in budding yeast that we engineered to study rhodopsin's light-activated signal transduction. Using this approach, we measured the mutational effects of over 1200 individual human rhodopsin mutants, generated by low-frequency random mutagenesis of the GPCR rhodopsin (RHO) gene. Analysis of the data in the context of rhodopsin's three-dimensional structure reveals that transmembrane helices are generally less tolerant to mutations compared to flanking helices that face the lipid bilayer, which suggest that mutational tolerance is contingent on both the local environment surrounding specific residues and the specific position of these residues in the protein structure. Comparison of functional scores from our screen to clinically identified rhodopsin disease variants found many pathogenic mutants to be loss of function. Lastly, functional scores from our assay were consistent with a complex counterion mechanism involved in ligand-binding and rhodopsin activation. Our results demonstrate that deep mutational scanning is possible for rhodopsin activation and can be an effective method for revealing properties of mutational tolerance that may be generalizable to other transmembrane proteins.

Modeling pH-Dependent Biomolecular Photochemistry.

The tuning mechanism of pH can be extremely challenging to model computationally in complex biological systems, especially with respect to the photochemical properties. This article reports a protocol aimed at modeling pH-dependent photodynamics using a combination of constant-pH molecular dynamics and semiclassical nonadiabatic molecular dynamics simulations. With retinal photoisomerization in Anabaena sensory rhodopsin (ASR) as a testbed, we show that our protocol produces pH-dependent photochemical properties, such as the isomerization quantum yield or decay rates. We decompose our results into single-titrated residue contributions, identifying some key tuning amino acids. Additionally, we assess the validity of the single protonation state picture to represent the system at a given pH and propose the most populated protein charge state as a compromise between cost and accuracy.

Vibrational Study of the Inward Proton Pump Xenorhodopsin NsXeR: Switch Order Determines Vectoriality.

Common proton pumps, e.g. HsBR and PR, transport protons out of the cell. Xenorhodopsins (XeR) were the first discovered microbial rhodopsins which come as natural inward proton pumps. In this work we combine steady-state (cryo-)FTIR and Raman spectroscopy with time-resolved IR and UV/Vis measurements to roadmap the inward proton transport of NsXeR and pinpoint the most important mechanistic features. Through the assignment of characteristic bands of the protein backbone, the retinal chromophore, the retinal Schiff base and D220, we could follow the switching processes for proton accessibility in accordance with the isomerization / switch / transfer model. The corresponding transient IR signatures suggest that the initial assignment of D220 as the proton acceptor needs to be questioned due to the temporal mismatch of the Schiff base and D220 protonation steps. The switching events in the K-L and MCP-MEC transitions are finely tuned by changes of the protein backbone and rearrangements of the Schiff base. This finely tuned mechanism is disrupted at cryogenic temperatures, being reflected in the replacement of the previously reported long-lived intermediate GS* by an actual redshifted (O-like) intermediate.

Cation-π Interactions Greatly Influence Ion Transportability of the Light-Driven Sodium Pump KR2: A Molecular Dynamics Study.

Krokinobacter eikastus rhodopsin 2 (KR2) is a typical light-driven sodium pump. Although wild-type KR2 exhibits high Na+ selectivity, mutagenesis performed on the residues constituting the entrance enables permeation of K+ and Cs+, while the underlying mechanism remains elusive. This study presents a comprehensive molecular dynamics investigation, including force field optimization, metadynamics, and alchemical free energy methods, to explore the N61L/G263F mutant of KR2, which exhibits transportability for K+ and Cs+. The introduced Phe263 residue can directly promote ion binding at the entrance through cation-π interactions, while the N61L mutation can enhance ion binding at Phe46 by relieving steric hindrance. These results suggest that cation-π interactions may significantly influence the ion transportability and selectivity of KR2, which can provide important insights for protein engineering and the design of artificial ion transporters.

The Functionality of the DC Pair in a Rhodopsin Guanylyl Cyclase from Catenaria anguillulae.

Rhodopsin guanylyl cyclases (RGCs) belong to the class of enzymerhodopsins catalyzing the transition from GTP into the second messenger cGMP, whereas light-regulation of enzyme activity is mediated by a membrane-bound microbial rhodopsin domain, that holds the catalytic center inactive in the dark. Structural determinants for activation of the rhodopsin moiety eventually leading to catalytic activity are largely unknown. Here, we investigate the mechanistic role of the D283-C259 (DC) pair that is hydrogen bonded via a water molecule as a crucial functional motif in the homodimeric C. anguillulae RGC. Based on a structural model of the DC pair in the retinal binding pocket obtained by MD simulation, we analyzed formation and kinetics of early and late photocycle intermediates of the rhodopsin domain wild type and specific DC pair mutants by combined UV-Vis and FTIR spectroscopy at ambient and cryo-temperatures. By assigning specific infrared bands to S-H vibrations of C259 we are able to show that the DC pair residues are tightly coupled. We show that deprotonation of D283 occurs already in the inactive L state as a prerequisite for M state formation, whereas structural changes of C259 occur in the active M state and early cryo-trapped intermediates. We propose a comprehensive molecular model for formation of the M state that activates the catalytic moiety. It involves light induced changes in bond strength and hydrogen bonding of the DC pair residues from the early J state to the active M state and explains the retarding effect of C259 mutants.

Osmotic stress studies of G-protein-coupled receptor rhodopsin activation.

We summarize and critically review osmotic stress studies of the G-protein-coupled receptor rhodopsin. Although small amounts of structural water are present in these receptors, the effect of bulk water on their function remains uncertain. Studies of the influences of osmotic stress on the GPCR archetype rhodopsin have given insights into the functional role of water in receptor activation. Experimental work has discovered that osmolytes shift the metarhodopsin equilibrium after photoactivation, either to the active or inactive conformations according to their molar mass. At least 80 water molecules are found to enter rhodopsin in the transition to the photoreceptor active state. We infer that this movement of water is both necessary and sufficient for receptor activation. If the water influx is prevented, e.g., by large polymer osmolytes or by dehydration, then the receptor functional transition is back shifted. These findings imply a new paradigm in which rhodopsin becomes solvent swollen in the activation mechanism. Water thus acts as an allosteric modulator of function for rhodopsin-like receptors in lipid membranes.

A rapid, tag-free way to purify functional GPCRs.

G protein-coupled receptors (GPCRs) play diverse signaling roles and represent major pharmaceutical targets. Consequently, they are the focus of intense study, and numerous advances have been made in their handling and analysis. However, a universal way to purify GPCRs has remained elusive, in part because of their inherent instability when isolated from cells. To address this, we have developed a general, rapid, and tag-free way to purify GPCRs. The method uses short peptide analogs of the Gα subunit C terminus (Gα-CT) that are attached to chromatography beads (Gα-CT resin). Because the Gα-CT peptides bind active GPCRs with high affinity, the Gα-CT resin selectively purifies only active functional receptors. We use this method to purify both rhodopsin and the ß2-adrenergic receptor and show they can be purified in either active conformations or inactive conformations, simply by varying elution conditions. While simple in concept-leveraging the conserved GPCR-Gα-CT binding interaction for the purpose of GPCR purification-we think this approach holds excellent potential to isolate functional receptors for a myriad of uses, from structural biology to proteomics.

Retinal-Based Anion Pump from the Cyanobacterium Tolypothrix campylonemoides.

In this work, TcaR rhodopsin from the cyanobacterium Tolypothrix campylonemoides was characterized. Analysis of the amino acid sequence of TcaR revealed that this protein possesses a TSD motif that differs by only one amino acid from the TSA motif of the known halorhodopsin chloride pump. The TcaR protein was expressed in E. coli, purified, and incorporated into proteoliposomes and nanodiscs. Functional activity was measured by electric current generation through the planar bilayer lipid membranes (BLMs) with proteoliposomes adsorbed on one side of the membrane surface, as well as by fluorescence using the voltage-dependent dye oxonol VI. We have shown that TcaR rhodopsin functions as a powerful anion pump. Our results show that the novel microbial anion transporter, TcaR, deserves deeper investigation and may be of interest both for fundamental studies of membrane proteins and as a tool for optogenetics.

Similarities and Differences in Photochemistry of Type I and Type II Rhodopsins.

The diversity of the retinal-containing proteins (rhodopsins) in nature is extremely large. Fundamental similarity of the structure and photochemical properties unites them into one family. However, there is still a debate about the origin of retinal-containing proteins: divergent or convergent evolution? In this review, based on the results of our own and literature data, a comparative analysis of the similarities and differences in the photoconversion of the rhodopsin of types I and II is carried out. The results of experimental studies of the forward and reverse photoreactions of the bacteriorhodopsin (type I) and visual rhodopsin (type II) rhodopsins in the femto- and picosecond time scale, photo-reversible reaction of the octopus rhodopsin (type II), photovoltaic reactions, as well as quantum chemical calculations of the forward photoreactions of bacteriorhodopsin and visual rhodopsin are presented. The issue of probable convergent evolution of type I and type II rhodopsins is discussed.

Unique Vibrational Characteristics and Structures of the Photoexcited Retinal Chromophore in Ion-Pumping Rhodopsins.

Photoisomerization of an all-trans-retinal chromophore triggers ion transport in microbial ion-pumping rhodopsins. Understanding chromophore structures in the electronically excited (S1) state provides insights into the structural evolution on the potential energy surface of the photoexcited state. In this study, we examined the structure of the S1-state chromophore in Natronomonas pharaonis halorhodopsin (NpHR), a chloride ion-pumping rhodopsin, using time-resolved resonance Raman spectroscopy. The spectral patterns of the S1-state chromophore were completely different from those of the ground-state chromophore, resulting from unique vibrational characteristics and the structure of the S1 state. Mode assignments were based on a combination of deuteration shifts of the Raman bands and hybrid quantum mechanics-molecular mechanics calculations. The present observations suggest a weakened bond alternation in the π conjugation system. A strong hydrogen-out-of-plane bending band was observed in the Raman spectra of the S1-state chromophore in NpHR, indicating a twisted polyene structure. Similar frequency shifts for the C═N/C═C and C-C stretching modes of the S1-state chromophore in NpHR were observed in the Raman spectra of sodium ion-pumping and proton-pumping rhodopsins, suggesting that these unique features are common to the S1 states of ion-pumping rhodopsins.

Hydrophobic interaction between the TM1 and H8 is essential for rhodopsin trafficking to vertebrate photoreceptor outer segments.

A major unsolved question in vertebrate photoreceptor biology is the mechanism of rhodopsin transport to the outer segment. In rhodopsin-like class A G protein-coupled receptors, hydrophobic interactions between C-terminal α-helix 8 (H8), and transmembrane α-helix-1 (TM1) have been shown to be important for transport to the plasma membrane, however whether this interaction is important for rhodopsin transport to ciliary rod outer segments is not known. We examined the crystal structures of vertebrate rhodopsins and class A G protein-coupled receptors and found a conserved network of predicted hydrophobic interactions. In Xenopus rhodopsin (xRho), this interaction corresponds to F313, L317, and L321 in H8 and M57, V61, and L68 in TM1. To evaluate the role of H8-TM1 hydrophobic interactions in rhodopsin transport, we expressed xRho-EGFP where hydrophobic residues were mutated in Xenopus rods and evaluated the efficiency of outer segment enrichment. We found that substituting L317 and M57 with hydrophilic residues had the strongest impact on xRho mislocalization. Substituting hydrophilic amino acids at positions L68, F313, and L321 also had a significant impact. Replacing L317 with M resulted in significant mislocalization, indicating that the hydrophobic interaction between residues 317 and 57 is exquisitely sensitive. The corresponding experiment in bovine rhodopsin expressed in HEK293 cells had a similar effect, showing that the H8-TM1 hydrophobic network is essential for rhodopsin transport in mammalian species. Thus, for the first time, we show that a hydrophobic interaction between H8 and TM1 is critical for efficient rhodopsin transport to the vertebrate photoreceptor ciliary outer segment.