Atomistic design of microbial opsin-based blue-shifted optogenetics tools

Structural changes of Natronomonas pharaonis halorhodopsin in its late photocycle revealed by solid-state NMR spectroscopy

Natronomonas pharaonis halorhodopsin (NpHR) is a light-driven Cl- inward pump that is widely used as an optogenetic tool. Although NpHR is previously extensively studied, its Cl- uptake process is not well understood from the protein structure perspective, mainly because in crystalline lattice, it has been difficult to analyze the structural changes associated with the Cl- uptake process. In this study, we used solid-state NMR to analyze NpHR both in the Cl--bound and -free states under near-physiological transmembrane condition. Chemical shift perturbation analysis suggested that while the structural change caused by the Cl- depletion is widespread over the NpHR molecule, residues in the extracellular (EC) part of helix D exhibited significant conformational changes that may be related to the Cl- uptake process. By combining photochemical analysis and dynamic nuclear polarization (DNP)-enhanced solid-state NMR measurement on NpHR point mutants for the suggested residues, we confirmed their importance in the Cl- uptake process. In particular, we found the mutation at Ala165 position, located at the trimer interface, to an amino acid with bulky sidechain (A165V) significantly perturbs the late photocycle and disrupts its trimeric assembly in the Cl--free state as well as during the ion-pumping cycle under the photo-irradiated condition. This strongly suggested an outward movement of helix D at EC part, disrupting the trimer integrity. Together with the spectroscopic data and known NpHR crystal structures, we proposed a model that this helix movement is required for creating the Cl- entrance path on the extracellular surface of the protein and is crucial to the Cl- uptake process.

Robust optogenetic inhibition with red-light-sensitive anion-conducting channelrhodopsins

Channelrhodopsins (ChRs) are light-gated ion channels widely used to optically activate or silence selected electrogenic cells, such as individual brain neurons. Here, we describe identifying and characterizing a set of anion-conducting ChRs (ACRs) from diverse taxa and representing various branches of the ChR phylogenetic tree. The Mantoniella squamata ACR (MsACR1) showed high sensitivity to yellow-green light (λmax at 555 nm) and was further engineered for optogenetic applications. A single amino-acid substitution that mimicked red-light-sensitive rhodopsins like Chrimson shifted the photosensitivity 20 nm toward red light and accelerated photocurrent kinetics. Hence, it was named red and accelerated ACR, raACR. Both wild-type and mutant are capable optical silencers at low light intensities in mouse neurons in vitro and in vivo, while raACR offers a higher temporal resolution.

Exploring Protonation State, Ion Binding, and Photoactivated Channel Opening of an Anion Channelrhodopsin by Molecular Simulations

Channelrhodopsins are light-gated ion channels with a retinal chromophore found in microbes and are widely used in optogenetics, a field of neuroscience that utilizes light to regulate neuronal activity. GtACR1, an anion conducting channelrhodopsin derived from Guillardia theta, has attracted attention for its application as a neuronal silencer in optogenetics because of its high conductivity and selectivity. However, atomistic mechanisms of channel photoactivation and ion conduction have not yet been elucidated. In the present study, we investigated the molecular characteristics of GtACR1 and its photoactivation processes by molecular simulations. The QM/MM RWFE-SCF method which combines highly accurate quantum chemistry calculations with long-time molecular dynamics (MD) simulations were used to model protein structures of the wild-type and mutants with different protonation states of key groups and to calculate absorption energies for verification of the models. The QM/MM modeling together with MD simulations of free-energy calculations favors protonation of a key counterion carboxyl group of Asp234 with a strong binding of a chloride ion in the extracellular pocket in the dark state. A channel open state was also successfully modeled by the QM/MM RWFE-SCF free-energy optimizations, providing atomistic insights into the channel activation mechanism.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

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.

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.

Automatic Rhodopsin Modeling with Multiple Protonation Microstates

Automatic Rhodopsin Modeling (ARM) is a simulation protocol providing QM/MM models of rhodopsins capable of reproducing experimental electronic absorption and emission trends. Currently, ARM is restricted to a single protonation microstate for each rhodopsin model. Herein, we incorporate an extension of the minimal electrostatic model (MEM) into the ARM protocol to account for all relevant protonation microstates at a given pH. The new ARM+MEM protocol determines the most important microstates contributing to the description of the absorption spectrum. As a test case, we have applied this methodology to simulate the pH-dependent absorption spectrum of a toy model, showing that the single-microstate picture breaks down at certain pH values. Subsequently, we applied ARM+MEM toAnabaenasensory rhodopsin, confirming an improved description of its absorption spectrum when the titration of several key residues is considered.

Structural basis for ion selectivity in potassium-selective channelrhodopsins

KCR channelrhodopsins (K+-selective light-gated ion channels) have received attention as potential inhibitory optogenetic tools but more broadly pose a fundamental mystery regarding how their K+ selectivity is achieved. Here, we present 2.5-2.7 Å cryo-electron microscopy structures of HcKCR1 and HcKCR2 and of a structure-guided mutant with enhanced K+ selectivity. Structural, electrophysiological, computational, spectroscopic, and biochemical analyses reveal a distinctive mechanism for K+ selectivity; rather than forming the symmetrical filter of canonical K+ channels achieving both selectivity and dehydration, instead, three extracellular-vestibule residues within each monomer form a flexible asymmetric selectivity gate, while a distinct dehydration pathway extends intracellularly. Structural comparisons reveal a retinal-binding pocket that induces retinal rotation (accounting for HcKCR1/HcKCR2 spectral differences), and design of corresponding KCR variants with increased K+ selectivity (KALI-1/KALI-2) provides key advantages for optogenetic inhibition in vitro and in vivo. Thus, discovery of a mechanism for ion-channel K+ selectivity also provides a framework for next-generation optogenetics.

Raman Optical Activity of Retinal Chromophore in Sensory Rhodopsin II

Raman optical activity (ROA) spectroscopy was used to study the conformation of the retinal chromophore in sensory rhodopsin II (SRII), which is a blue-green light sensor of microbes. The ROA spectrum consisted of the negative vibrational bands of the chromophore, whose relative intensities are similar to those of the parent Raman spectrum. This spectral feature was explained by the left-handed helical twist of the retinal chromophore on the basis of quantum chemical calculations. On the other hand, we found that the chromophore conformation based on the crystal structures of SRII has a right-handed helical twist, which does not agree with the observation. This specific result suggests that the consistency with chiro-optical properties can be a key criterion for the accurate prediction and/or evaluation of chromophore conformation in retinal-binding proteins.

A blue-shifted anion channelrhodopsin from the Colpodellida alga Vitrella brassicaformis

Microbial rhodopsins, a family of photoreceptive membrane proteins containing the chromophore retinal, show a variety of light-dependent molecular functions. Channelrhodopsins work as light-gated ion channels and are widely utilized for optogenetics, which is a method for controlling neural activities by light. Since two cation channelrhodopsins were identified from the chlorophyte alga Chlamydomonas reinhardtii, recent advances in genomic research have revealed a wide variety of channelrhodopsins including anion channelrhodopsins (ACRs), describing their highly diversified molecular properties (e.g., spectral sensitivity, kinetics and ion selectivity). Here, we report two channelrhodopsin-like rhodopsins from the Colpodellida alga Vitrella brassicaformis, which are phylogenetically distinct from the known channelrhodopsins. Spectroscopic and electrophysiological analyses indicated that these rhodopsins are green- and blue-sensitive pigments (λ_max =  ~ 550 and ~ 440 nm) that exhibit light-dependent ion channeling activities. Detailed electrophysiological analysis revealed that one of them works as a monovalent anion (Cl^-, Br^- and NO₃^-) channel and we named it V. brassicaformis anion channelrhodopsin-2, VbACR2. Importantly, the absorption maximum of VbACR2 (~ 440 nm) is blue-shifted among the known ACRs. Thus, we identified the new blue-shifted ACR, which leads to the expansion of the molecular diversity of ACRs.

Pump-like channelrhodopsins: Not just bridging the gap between ion pumps and ion channels

Channelrhodopsins are microbial rhodopsins that work as light-gated ion channels. Their importance has become increasingly recognized due to their ability to control the membrane potential of specific cells in a light-dependent manner. This technology, termed optogenetics, has revolutionized neuroscience, and numerous channelrhodopsin variants have been isolated or engineered to expand the utility of optogenetics. Pump-like channelrhodopsins (PLCRs), one of the recently discovered channelrhodopsin subfamilies, have attracted broad attention due to their high sequence similarity to ion-pumping rhodopsins and their distinct properties, such as high light sensitivity and ion selectivity. In this review, we summarize the current understanding of the structure-function relationships of PLCRs and discuss the challenges and opportunities of channelrhodopsin research.

Properties of a Single Amino Acid Residue in the Third Transmembrane Domain Determine the Kinetics of Ambient Light-Sensitive Channelrhodopsin

Channelrhodopsins have been utilized in gene therapy to restore vision in patients with retinitis pigmentosa and their channel kinetics are an important factor to consider in such applications. We investigated the channel kinetics of ComV1 variants with different amino acid residues at the 172nd position. Patch clamp methods were used to record the photocurrents induced by stimuli from diodes in HEK293 cells transfected with plasmid vectors. The channel kinetics (τon and τoff) were considerably altered by the replacement of the 172nd amino acid and was dependent on the amino acid characteristics. The size of amino acids at this position correlated with τon and decay, whereas the solubility correlated with τon and τoff. Molecular dynamic simulation indicated that the ion tunnel constructed by H172, E121, and R306 widened due to H172A variant, whereas the interaction between A172 and the surrounding amino acids weakened compared with H172. The bottleneck radius of the ion gate constructed with the 172nd amino acid affected the photocurrent and channel kinetics. The 172nd amino acid in ComV1 is a key residue for determining channel kinetics as its properties alter the radius of the ion gate. Our findings can be used to improve the channel kinetics of channelrhodopsins.

All-optical physiology resolves a synaptic basis for behavioral timescale plasticity

Learning has been associated with modifications of synaptic and circuit properties, but the precise changes storing information in mammals have remained largely unclear. We combined genetically targeted voltage imaging with targeted optogenetic activation and silencing of pre- and post-synaptic neurons to study the mechanisms underlying hippocampal behavioral timescale plasticity. In mice navigating a virtual-reality environment, targeted optogenetic activation of individual CA1 cells at specific places induced stable representations of these places in the targeted cells. Optical elicitation, recording, and modulation of synaptic transmission in behaving mice revealed that activity in presynaptic CA2/3 cells was required for the induction of plasticity in CA1 and, furthermore, that during induction of these place fields in single CA1 cells, synaptic input from CA2/3 onto these same cells was potentiated. These results reveal synaptic implementation of hippocampal behavioral timescale plasticity and define a methodology to resolve synaptic plasticity during learning and memory in behaving mammals.

A kinetic-optimized CoChR variant with enhanced high-frequency spiking fidelity

Channelrhodopsins are a promising toolset for noninvasive optical manipulation of genetically identifiable neuron populations. Existing channelrhodopsins have generally suffered from a trade-off between two desired properties: fast channel kinetics and large photocurrent. Such a trade-off hinders spatiotemporally precise optogenetic activation during both one-photon and two-photon photostimulation. Furthermore, the simultaneous use of spectrally separated genetically encoded indicators and channelrhodopsins has generally suffered from non-negligible crosstalk in photocurrent or fluorescence. These limitations have hindered crosstalk-free dual-channel experiments needed to establish relationships between multiple neural populations. Recent large-scale transcriptome sequencing revealed one potent optogenetic actuator, the channelrhodopsin from species Chloromonas oogama (CoChR), which possessed high cyan light-driven photocurrent but slow channel kinetics. We rationally designed and engineered a kinetic-optimized CoChR variant that was faster than native CoChR while maintaining large photocurrent amplitude. When expressed in cultured hippocampal pyramidal neurons, our CoChR variant improved high-frequency spiking fidelity under one-photon illumination. Our CoChR variant's blue-shifted excitation spectrum enabled simultaneous cyan photostimulation and red calcium imaging with negligible photocurrent crosstalk.

Structural insights into light-driven anion pumping in cyanobacteria

Transmembrane ion transport is a key process in living cells. Active transport of ions is carried out by various ion transporters including microbial rhodopsins (MRs). MRs perform diverse functions such as active and passive ion transport, photo-sensing, and others. In particular, MRs can pump various monovalent ions like Na+, K+, Cl-, I-, NO3-. The only characterized MR proposed to pump sulfate in addition to halides belongs to the cyanobacterium Synechocystis sp. PCC 7509 and is named Synechocystis halorhodopsin (SyHR). The structural study of SyHR may help to understand what makes an MR pump divalent ions. Here we present the crystal structure of SyHR in the ground state, the structure of its sulfate-bound form as well as two photoreaction intermediates, the K and O states. These data reveal the molecular origin of the unique properties of the protein (exceptionally strong chloride binding and proposed pumping of divalent anions) and sheds light on the mechanism of anion release and uptake in cyanobacterial halorhodopsins. The unique properties of SyHR highlight its potential as an optogenetics tool and may help engineer different types of anion pumps with applications in optogenetics.

Cryo-EM structures of the channelrhodopsin ChRmine in lipid nanodiscs

Microbial channelrhodopsins are light-gated ion channels widely used for optogenetic manipulation of neuronal activity. ChRmine is a bacteriorhodopsin-like cation channelrhodopsin (BCCR) more closely related to ion pump rhodopsins than other channelrhodopsins. ChRmine displays unique properties favorable for optogenetics including high light sensitivity, a broad, red-shifted activation spectrum, cation selectivity, and large photocurrents, while its slow closing kinetics impedes some applications. The structural basis for ChRmine function, or that of any other BCCR, is unknown. Here, we present cryo-EM structures of ChRmine in lipid nanodiscs in apo (opsin) and retinal-bound (rhodopsin) forms. The structures reveal an unprecedented trimeric architecture with a lipid filled central pore. Large electronegative cavities on either side of the membrane facilitate high conductance and selectivity for cations over protons. The retinal binding pocket structure suggests channel properties could be tuned with mutations and we identify ChRmine variants with ten-fold decreased and two-fold increased closing rates. A T119A mutant shows favorable properties relative to wild-type and previously reported ChRmine variants for optogenetics. These results provide insight into structural features that generate an ultra-potent microbial opsin and provide a platform for rational engineering of channelrhodopsins with improved properties that could expand the scale, depth, and precision of optogenetic experiments.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

Archaeal Lipids Regulating the Trimeric Structure Dynamics of Bacteriorhodopsin for Efficient Proton Release and Uptake

S-TGA-1 and PGP-Me are native archaeal lipids associated with the bacteriorhodopsin (bR) trimer and contribute to protein stabilization and native dynamics for proton transfer. However, little is known about the underlying molecular mechanism of how these lipids regulate bR trimerization and efficient photocycling. Here, we explored the specific binding of S-TGA-1 and PGP-Me with the bR trimer and elucidated how specific interactions modulate the bR trimeric structure and proton release and uptake using long-term atomistic molecular dynamic simulations. Our results showed that S-TGA-1 and PGP-Me are essential for stabilizing the bR trimer and maintaining the coherent conformational dynamics necessary for proton transfer. The specific binding of S-TGA-1 with W80 and K129 regulates proton release on the extracellular surface by forming a “Glu-shared” model. The interaction of PGP-Me with K40 ensures proton uptake by accommodating the conformation of the helices to recruit enough water molecules on the cytoplasmic side. The present study results could fill in the theoretical gaps of studies on the functional role of archaeal lipids and could provide a reference for other membrane proteins containing similar archaeal lipids.

Advancements in the Quest to Map, Monitor, and Manipulate Neural Circuitry

Neural circuits and the cells that comprise them represent the functional units of the brain. Circuits relay and process sensory information, maintain homeostasis, drive behaviors, and facilitate cognitive functions such as learning and memory. Creating a functionally-precise map of the mammalian brain requires anatomically tracing neural circuits, monitoring their activity patterns, and manipulating their activity to infer function. Advancements in cell-type-specific genetic tools allow interrogation of neural circuits with increased precision. This review provides a broad overview of recombination-based and activity-driven genetic targeting approaches, contemporary viral tracing strategies, electrophysiological recording methods, newly developed calcium, and voltage indicators, and neurotransmitter/neuropeptide biosensors currently being used to investigate circuit architecture and function. Finally, it discusses methods for acute or chronic manipulation of neural activity, including genetically-targeted cellular ablation, optogenetics, chemogenetics, and over-expression of ion channels. With this ever-evolving genetic toolbox, scientists are continuing to probe neural circuits with increasing resolution, elucidating the structure and function of the incredibly complex mammalian brain.

Engineering rhodopsins' activation spectra using a FRET-based approach

In the past decade, optogenetics has become a nearly ubiquitous tool in neuroscience because it enables researchers to manipulate neural activity with high temporal resolution and genetic specificity. Rational engineering of optogenetic tools has produced channelrhodopsins with a wide range of kinetics and photocurrent magnitude. Genome mining for previously unidentified species of rhodopsin has uncovered optogenetic tools with diverse spectral sensitivities. However, rational engineering of a rhodopsin has thus far been unable to re-engineer spectral sensitivity while preserving full photocurrent. Here, we developed and characterized ChroME-mTFP, a rhodopsin-fluorescent protein fusion that drives photocurrent through Förster resonance energy transfer (FRET). This FRET-opsin mechanism artificially broadened the activation spectrum of the blue-green-light-activated rhodopsin ChroME by approximately 50 nm, driving higher photocurrent at blue-shifted excitation wavelengths without sacrificing kinetics. The excitation spectra's increase at short wavelengths enabled us to optogenetically excite neurons at lower excitation powers with shorter wavelengths of light. Increasing this rhodopsin's sensitivity to shorter, bluer wavelengths pushes it toward dual-channel, crosstalk-free optogenetic stimulation and imaging with green-light-activated sensors. However, this iteration of FRET-opsin suffers from some imaging-light-induced photocurrent crosstalk from green or yellow light due to maintained, low-efficiency excitation at longer wavelengths.

Hybrid QM/MM Free-Energy Evaluation of Drug-Resistant Mutational Effect on the Binding of an Inhibitor Indinavir to HIV-1 Protease

A human immunodeficiency virus-1 (HIV-1) protease is a homodimeric aspartic protease essential for the replication of HIV. The HIV-1 protease is a target protein in drug discovery for antiretroviral therapy, and various inhibitor molecules of transition state analogues have been developed. However, serious drug-resistant mutants have emerged. For understanding the molecular mechanism of the drug resistance, an accurate examination of the impacts of the mutations on ligand binding and enzymatic activity is necessary. Here, we present a molecular simulation study on the ligand binding of indinavir, a potent transition state analogue inhibitor, to the wild-type protein and a V82T/I84V drug-resistant mutant of the HIV-1 protease. We employed a hybrid ab initio quantum mechanical/molecular mechanical (QM/MM) free-energy optimization technique which combines a highly accurate QM description of the ligand molecule and its interaction with statistically ample conformational sampling of the MM protein environment by long-time molecular dynamics simulations. Through the free-energy calculations of protonation states of catalytic groups at the binding pocket and of the ligand-binding affinity changes upon the mutations, we successfully reproduced the experimentally observed significant reduction of the binding affinity upon the drug-resistant mutations and elucidated the underlying molecular mechanism. The present study opens the way for understanding the molecular mechanism of drug resistance through the direct quantitative comparison of ligand binding and enzymatic reaction with the same accuracy.

Structural basis for channel conduction in the pump-like channelrhodopsin ChRmine

ChRmine, a recently discovered pump-like cation-conducting channelrhodopsin, exhibits puzzling properties (large photocurrents, red-shifted spectrum, and extreme light sensitivity) that have created new opportunities in optogenetics. ChRmine and its homologs function as ion channels but, by primary sequence, more closely resemble ion pump rhodopsins; mechanisms for passive channel conduction in this family have remained mysterious. Here, we present the 2.0 Å resolution cryo-EM structure of ChRmine, revealing architectural features atypical for channelrhodopsins: trimeric assembly, a short transmembrane-helix 3, a twisting extracellular-loop 1, large vestibules within the monomer, and an opening at the trimer interface. We applied this structure to design three proteins (rsChRmine and hsChRmine, conferring further red-shifted and high-speed properties, respectively, and frChRmine, combining faster and more red-shifted performance) suitable for fundamental neuroscience opportunities. These results illuminate the conduction and gating of pump-like channelrhodopsins and point the way toward further structure-guided creation of channelrhodopsins for applications across biology.

Mechanisms of long-distance allosteric couplings in proton-binding membrane transporters

Membrane transporters that use proton binding and proton transfer for function couple local protonation change with changes in protein conformation and water dynamics. Changes of protein conformation might be required to allow transient formation of hydrogen-bond networks that bridge proton donor and acceptor pairs separated by long distances. Inter-helical hydrogen-bond networks adjust rapidly to protonation change, and ensure rapid response of the protein structure and dynamics. Membrane transporters with known three-dimensional structures and proton-binding groups inform on general principles of protonation-coupled protein conformational dynamics. Inter-helical hydrogen bond motifs between proton-binding carboxylate groups and a polar sidechain are observed in unrelated membrane transporters, suggesting common principles of coupling protonation change with protein conformational dynamics.

Mutational analyses identify a single amino acid critical for color tuning in proteorhodopsins

Microbial rhodopsins are light-activated proteins that contain seven transmembrane alpha-helices. Spectral tuning in microbial rhodopsins is a useful optogenetic tool. In this study, we report a new site that controls spectral tuning. In the proteorhodopsins ISR34 and ISR36, a single amino-acid substitution at Cys189 caused an absorption maximum shift of 44 nm, indicating spectral tuning at a specific site. Comparison of single amino acid substitutions was conducted using photochemical and photobiological approaches. The maximum absorption for red-shift was measured for mutations at positions 189 and 105 in ISR34, both residues being equally important. Structural changes resulting from amino acid substitutions are related to pK_a  values, pumping activity, and spectral tuning.

Emerging Diversity of Channelrhodopsins and Their Structure-Function Relationships

Cation and anion channelrhodopsins (CCRs and ACRs, respectively) from phototactic algae have become widely used as genetically encoded molecular tools to control cell membrane potential with light. Recent advances in polynucleotide sequencing, especially in environmental samples, have led to identification of hundreds of channelrhodopsin homologs in many phylogenetic lineages, including non-photosynthetic protists. Only a few CCRs and ACRs have been characterized in detail, but there are indications that ion channel function has evolved within the rhodopsin superfamily by convergent routes. The diversity of channelrhodopsins provides an exceptional platform for the study of structure-function evolution in membrane proteins. Here we review the current state of channelrhodopsin research and outline perspectives for its further development.

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.

Rhodopsins at a glance

Rhodopsins are photoreceptive membrane proteins consisting of a common heptahelical transmembrane architecture that contains a retinal chromophore. Rhodopsin was first discovered in the animal retina in 1876, but a different type of rhodopsin, bacteriorhodopsin, was reported to be present in the cell membrane of an extreme halophilic archaeon, Halobacterium salinarum, 95 years later. Although these findings were made by physiological observation of pigmented tissue and cell bodies, recent progress in genomic and metagenomic analyses has revealed that there are more than 10,000 microbial rhodopsins and 9000 animal rhodopsins with large diversity and tremendous new functionality. In this Cell Science at a Glance article and accompanying poster, we provide an overview of the diversity of functions, structures, color discrimination mechanisms and optogenetic applications of these two rhodopsin families, and will also highlight the third distinctive rhodopsin family, heliorhodopsin.

Pro219 is an electrostatic color determinant in the light-driven sodium pump KR2

Color tuning in animal and microbial rhodopsins has attracted the interest of many researchers, as the color of their common retinal chromophores is modulated by the amino acid residues forming the chromophore cavity. Critical cavity amino acid residues are often called “color switches”, as the rhodopsin color is effectively tuned through their substitution. Well-known color switches are the L/Q and A/TS switches located in the C and G helices of the microbial rhodopsin structure respectively. Recently, we reported on a third G/P switch located in the F helix of the light-driven sodium pumps of KR2 and JsNaR causing substantial spectral red-shifts in the latter with respect to the former. In order to investigate the molecular-level mechanism driving such switching function, here we present an exhaustive mutation, spectroscopic and computational investigation of the P219X mutant set of KR2. To do so, we study the changes in the absorption band of the 19 possible mutants and construct, semi-automatically, the corresponding hybrid quantum mechanics/molecular mechanics models. We found that the P219X feature a red-shifted light absorption with the only exception of P219R. The analysis of the corresponding models indicate that the G/P switch induces red-shifting variations via electrostatic interactions, while replacement-induced chromophore geometrical (steric) distortions play a minor role. However, the same analysis indicates that the P219R blue-shifted variant has a more complex origin involving both electrostatic and steric changes accompanied by protonation state and hydrogen bond networks modifications. These results make it difficult to extract simple rules or formulate theories for predicting how a switch operates without considering the atomistic details and environmental consequences of the side chain replacement.

Nakajima, Pedraza-González et al. provide a comprehensive investigation of amino acid mutations at position 219 of the sodium pump rhodopsin, KR2, and their role in the color tuning of the retinal chromophore. They prepared P219X (X= A, C, D, E, F, G, H, I, K, L, M, N, Q, R, S, T, V, W, Y) mutants of KR2, and find that all mutants are red-shifted, except for P219R, highlighting its role as a color determinant in the light-driven pump KR2.

Structure of a retinal chromophore of dark-adapted middle rhodopsin as studied by solid-state nuclear magnetic resonance spectroscopy

Middle rhodopsin (MR) found from the archaeon Haloquadratum walsbyi is evolutionarily located between two different types of rhodopsins, bacteriorhodopsin (BR) and sensory rhodopsin II (SRII). Some isomers of the chromophore retinal and the photochemical reaction of MR are markedly different from those of BR and SRII. In this study, to obtain the structural information regarding its active center (i.e., retinal), we subjected MR embedded in lipid bilayers to solid-state magic-angle spinning nuclear magnetic resonance (NMR) spectroscopy. The analysis of the isotropic ¹³C chemical shifts of the retinal chromophore revealed the presence of three types of retinal configurations of dark-adapted MR: (13-trans, 15-anti (all-trans)), (13-cis, 15-syn), and 11-cis isomers. The higher field resonance of the 20-C methyl carbon in the all-trans retinal suggested that Trp182 in MR has an orientation that is different from that in other microbial rhodopsins, owing to the changes in steric hindrance associated with the 20-C methyl group in retinal. ¹³Cζ signals of Tyr185 in MR for all-trans and 13-cis, 15-syn isomers were discretely observed, representing the difference in the hydrogen bond strength of Tyr185. Further, ¹⁵N NMR analysis of the protonated Schiff base corresponding to the all-trans and 13-cis, 15-syn isomers in MR showed a strong electrostatic interaction with the counter ion. Therefore, the resulting structural information exhibited the property of stable retinal conformations of dark-adapted MR.

Structural Factors Determining the Absorption Spectrum of Channelrhodopsins: A Case Study of the Chimera C1C2

Channelrhodopsins are photosensitive proteins that trigger flagella motion in single-cell algae and have been successfully utilized in optogenetic applications. In optogenetics, light is used to activate neural cells in living organisms, which can be achieved by exploiting the ion channel signaling of channelrhodopsins. Tailoring channelrhodopsins for such applications includes the tuning of the absorption maximum. In order to establish rational design and to obtain a desired spectral shift, a basic understanding of the absorption spectrum is required. We have studied the chimera C1C2 as a representative of this protein family and the first member with an available crystal structure. For this purpose, we sampled the conformations of C1C2 using quantum mechanical/molecular mechanical molecular dynamics and subjected the resulting snapshots of the trajectory to excitation energy calculations using ADC(2) and simplified time-dependent density functional theory. In contrast to previous reports, we found that different hydrogen-bonding networks-involving the retinal protonated Schiff base, the putative counterions E162 and D292, and water molecules-had only a small impact on the absorption spectrum. However, in the case of deprotonated E162, increasing the distance to the Schiff base hydrogen-bonding partner led to a systematic blue shift. The β-ionone ring rotation was identified as another important contributor. Yet the most important factors were found to be the bond length alternation and bond order alternation that were linearly correlated to the absorption maximum by up to 62 and 82%, respectively. We ascribe this novel insight into the structural basis of the absorption spectrum to our enhanced protein setup that includes membrane embedding as well as long and extensive sampling.

Crystal structure of schizorhodopsin reveals mechanism of inward proton pumping

Schizorhodopsins (SzRs), a new rhodopsin family identified in Asgard archaea, are phylogenetically located at an intermediate position between type-1 microbial rhodopsins and heliorhodopsins. SzRs work as light-driven inward H⁺ pumps as xenorhodopsins in bacteria. Although E81 plays an essential role in inward H⁺ release, the H⁺ is not metastably trapped in such a putative H⁺ acceptor, unlike the other H⁺ pumps. It remains elusive why SzR exhibits different kinetic behaviors in H⁺ release. Here, we report the crystal structure of SzR AM_5_00977 at 2.1 Å resolution. The SzR structure superimposes well on that of bacteriorhodopsin rather than heliorhodopsin, suggesting that SzRs are classified with type-1 rhodopsins. The structure-based mutagenesis study demonstrated that the residues N100 and V103 around the β-ionone ring are essential for color tuning in SzRs. The cytoplasmic parts of transmembrane helices 2, 6, and 7 are shorter than those in the other microbial rhodopsins, and thus E81 is located near the cytosol and easily exposed to the solvent by light-induced structural change. We propose a model of untrapped inward H⁺ release; H⁺ is released through the water-mediated transport network from the retinal Schiff base to the cytosol by the side of E81. Moreover, most residues on the H⁺ transport pathway are not conserved between SzRs and xenorhodopsins, suggesting that they have entirely different inward H⁺ release mechanisms.

Specific residues in the cytoplasmic domain modulate photocurrent kinetics of channelrhodopsin from Klebsormidium nitens

Channelrhodopsins (ChRs) are light-gated ion channels extensively applied as optogenetics tools for manipulating neuronal activity. All currently known ChRs comprise a large cytoplasmic domain, whose function is elusive. Here, we report the cation channel properties of KnChR, one of the photoreceptors from a filamentous terrestrial alga Klebsormidium nitens, and demonstrate that the cytoplasmic domain of KnChR modulates the ion channel properties. KnChR is constituted of a 7-transmembrane domain forming a channel pore, followed by a C-terminus moiety encoding a peptidoglycan binding domain (FimV). Notably, the channel closure rate was affected by the C-terminus moiety. Truncation of the moiety to various lengths prolonged the channel open lifetime by more than 10-fold. Two Arginine residues (R287 and R291) are crucial for altering the photocurrent kinetics. We propose that electrostatic interaction between the rhodopsin domain and the C-terminus domain accelerates the channel kinetics. Additionally, maximal sensitivity was exhibited at 430 and 460 nm, the former making KnChR one of the most blue-shifted ChRs characterized thus far, serving as a novel prototype for studying the molecular mechanism of color tuning of the ChRs. Furthermore, KnChR would expand the optogenetics tool kit, especially for dual light applications when short-wavelength excitation is required.

Tashiro et al. describe a new channelrhodopsin variant from a terrestrial algal species and the role of the C-terminal domain in regulatory function. This far-blue-shifted channelrhodopsin may contribute to optogenetic tool research in the future.

Ab Initio Evaluation of the Redox Potential of Cytochrome c

Various biochemical activities of metabolism and biosynthesis are fulfilled by redox processes with explicit electron exchange, which furnish redox enzymes with high chemical reactivity. However, theoretical investigation of a redox process, which simultaneously involves a complex electronic change at a redox metal center and conformational reorganization of the surrounding protein environment coupled to the electronic change, requires computationally conflicting approaches, highly accurate quantum chemical calculations, and long-time molecular dynamics (MD) simulations, limiting the physicochemical understanding of biological redox processes. Here, we theoretically examined a redox process of cytochrome c by means of a hybrid molecular simulation technique, which enables one to consistently treat the redox center at the ab initio quantum chemistry level of theory and the protein reorganization with long-time MD simulations on the microsecond timescale. The calculations successfully evaluated a large absolute redox potential, 4.34 eV, with errors of only 0.03 to 0.34 eV to the experimental ones without any problem-specific empirical parameters. Through the long-time MD sampling, large and nonlinear reorganization of the protein environment was unveiled and the molecular determinants for the redox potential were identified. The present ab initio approach significantly expands the applicability of theoretical investigation to biological redox systems with more electronically complicated redox centers such as polynuclear transition metal complexes.

Structure-Function Relationship of Channelrhodopsins

Ion-translocating rhodopsins, especially channelrhodopsins (ChRs), have attracted broad attention as a powerful tool to modulate the membrane potential of cells with light (optogenetics). Because of recent biophysical, spectroscopic, and computational studies, including the structural determination of cation and anion ChRs, our understanding of the molecular mechanism underlying light-gated ion conduction has been greatly advanced. In this chapter, I first describe the background of rhodopsin family proteins including ChR, and how the optogenetics technology has been established from the discovery of first ChR in 2002. I later introduce the recent findings of the structure-function relationship of ChR by comparing the crystal structures of cation and anion ChRs. I further discuss the future goal in the fields of ChR research and optogenetic tool development.

Detection of neuropeptides in vivo and open questions for current and upcoming fluorescent sensors for neuropeptides

During a stress response, various neuropeptides are secreted in a spatiotemporally coordinated way in the brain. For a precise understanding of peptide functions in a stress response, it is important to investigate when and where they are released, how they diffuse, and how they are broken down in the brain. In the past two decades, genetically encoded fluorescent calcium indicators have greatly advanced our knowledge of the functions of specific neuronal activity in regulation of behavioral changes and physiological responses during stress. In addition, various kinds of structural information on G-protein-coupled receptors (GPCRs) for neuropeptides have been revealed. Recently, genetically encoded fluorescent sensors have been developed for detection of neurotransmitters by making use of conformational changes induced by ligand binding. In this review, we summarize the recent and upcoming advances of techniques for detection of neuropeptides and then present several open questions that will be solved by application of recent or upcoming technical advances in detection of neuropeptides in vivo.

Diversity, Mechanism, and Optogenetic Application of Light-Driven Ion Pump Rhodopsins

Ion-transporting microbial rhodopsins are widely used as major molecular tools in optogenetics. They are categorized into light-gated ion channels and light-driven ion pumps. While the former passively transport various types of cations and anions in a light-dependent manner, light-driven ion pumps actively transport specific ions, such as H⁺, Na⁺, Cl^-, against electrophysiological potential by using light energy. Since the ion transport by these pumps induces hyperpolarization of membrane potential and inhibit neural firing, light-driven ion-pumping rhodopsins are mostly applied as inhibitory optogenetics tools. Recent progress in genome and metagenome sequencing identified more than several thousands of ion-pumping rhodopsins from a wide variety of microbes, and functional characterization studies has been revealing many new types of light-driven ion pumps one after another. Since light-gated channels were reviewed in other chapters in this book, here the rapid progress in functional characterization, molecular mechanism study, and optogenetic application of ion-pumping rhodopsins were reviewed.

Cardiac optogenetics: a decade of enlightenment

The electromechanical function of the heart involves complex, coordinated activity over time and space. Life-threatening cardiac arrhythmias arise from asynchrony in these space-time events; therefore, therapies for prevention and treatment require fundamental understanding and the ability to visualize, perturb and control cardiac activity. Optogenetics combines optical and molecular biology (genetic) approaches for light-enabled sensing and actuation of electrical activity with unprecedented spatiotemporal resolution and parallelism. The year 2020 marks a decade of developments in cardiac optogenetics since this technology was adopted from neuroscience and applied to the heart. In this Review, we appraise a decade of advances that define near-term (immediate) translation based on all-optical electrophysiology, including high-throughput screening, cardiotoxicity testing and personalized medicine assays, and long-term (aspirational) prospects for clinical translation of cardiac optogenetics, including new optical therapies for rhythm control. The main translational opportunities and challenges for optogenetics to be fully embraced in cardiology are also discussed.

Mechanism of absorption wavelength shifts in anion channelrhodopsin-1 mutants

Using a quantum mechanical/molecular mechanical approach, we show the mechanisms of how the protein environment of Guillardia theta anion channelrhodopsin-1 (GtACR1) can shift the absorption wavelength. The calculated absorption wavelengths for GtACR1 mutants, M105A, C133A, and C237A are in agreement with experimentally measured wavelengths. Among 192 mutant structures investigated, mutations at Thr101, Cys133, Pro208, and Cys237 are likely to increase the absorption wavelength. In particular, T101A GtACR1 was expressed in HEK293T cells. The measured absorption wavelength is 10 nm higher than that of wild type, consistent with the calculated wavelength. (i) Removal of a polar residue from the Schiff base moiety, (ii) addition of a polar or acidic residue to the β-ionone ring moiety, and (iii) addition of a bulky residue to increase the planarity of the β-ionone and Schiff base moieties are the basis of increasing the absorption wavelength.

Insights into the Protein Functions and Absorption Wavelengths of Microbial Rhodopsins

Using a quantum mechanical/molecular mechanical approach, the absorption wavelength of the retinal Schiff base was calculated based on 13 microbial rhodopsin crystal structures. The results showed that the protein electrostatic environment decreases the absorption wavelength significantly in the cation-conducting rhodopsin but only slightly in the sensory rhodopsin. Among the microbial rhodopsins with different functions, the differences in the absorption wavelengths are caused by differences in the arrangement of the charged residues at the retinal Schiff base binding moiety, namely, one or two counterions at the three common positions. Among the microbial rhodopsins with similar functions, the differences in the polar residues at the retinal Schiff base binding site are responsible for the differences in the absorption wavelengths. Counterions contribute to an absorption wavelength shift of 50-120 nm, whereas polar groups contribute to a shift of up to ∼10 nm. It seems likely that protein function is directly associated with the absorption wavelength in microbial rhodopsins.

The Unlimited Potential of Microbial Rhodopsins as Optical Tools

Microbial rhodopsins, a photoactive membrane protein family, serve as fundamental tools for optogenetics, an innovative technology for controlling biological activities with light. Microbial rhodopsins are widely distributed in nature and have a wide variety of biological functions. Regardless of the many different known types of microbial rhodopsins, only a few of them have been used in optogenetics to control neural activity to understand neural networks. The efforts of our group have been aimed at identifying and characterizing novel rhodopsins from nature and also at engineering novel variant rhodopsins by rational design. On the basis of the molecular and functional characteristics of those novel rhodopsins, we have proposed new rhodopsin-based optogenetics tools to control not only neural activities but also "non-neural" activities. In this Perspective, we introduce the achievements and summarize future challenges in creating optogenetics tools using rhodopsins. The implementation of optogenetics deep inside an in vivo brain is the well-known challenge for existing rhodopsins. As a perspective to address this challenge, we introduce innovative optical illumination techniques using wavefront shaping that can reinforce the low light sensitivity of the rhodopsins and realize deep-brain optogenetics. The applications of our optogenetics tools could be extended to manipulate non-neural biological activities such as gene expression, apoptosis, energy production, and muscle contraction. We also discuss the potentially unlimited biotechnological applications of microbial rhodopsins in the future such as in photovoltaic devices and in drug delivery systems. We believe that advances in the field will greatly expand the potential uses of microbial rhodopsins as optical tools.

Salt bridge impact on global rigidity and thermostability in thermophilic citrate synthase

It has been suggested that structural rigidity is connected to thermostability, e.g. in enzymes from thermophilic microorganisms. We examine the importance of correctly handling salt bridges, and interactions which we term 'strong polars', when constructing the constraint network for global rigidity analysis in these systems. Through a comparison of rigidity in citrate synthases, we clarify the relationship between rigidity and thermostability. In particular, with our corrected handling of strong polar interactions, the difference in rigidity between mesophilic and thermophilic structures is detected more clearly than in previous studies. The increase in rigidity did not detract from the functional flexibility of the active site in all systems once their respective temperature range had been reached. We then examine the distribution of salt bridges in thermophiles that were previously unaccounted for in flexibility studies. We show that in hyperthermophiles these have stabilising roles in the active site; occuring in close proximity to key residues involved in catalysis and binding of the protein.

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.

Assessment of MC-PDFT Excitation Energies for a Set of QM/MM Models of Rhodopsins

A methodology for the automatic production of quantum mechanical/molecular mechanical (QM/MM) models of retinal-binding rhodopsin proteins and subsequent prediction of their spectroscopic properties has been proposed recently by some of the authors. The technology employed for the evaluation of the excitation energies is called Automatic Rhodopsin Modeling (ARM), and it involves the use of the complete active space self-consistent field (CASSCF) method followed by a multiconfiguration second-order perturbation theory (in particular, CASPT2) calculation of external correlation energies. Although it was shown that ARM is capable of successfully reproducing and predicting spectroscopic property trends in chromophore-embedding protein sets, practical applications of such technology are limited by the high computational costs of the multiconfiguration perturbation theory calculations. In the present work we benchmark the more affordable multiconfiguration pair-density functional theory (MC-PDFT) method whose accuracy has been recently validated for retinal chromophores in the gas phase, indicating that MC-PDFT could potentially be used to analyze large (e.g., few hundreds) sets of rhodopsin proteins. Here, we test this theory for a set of rhodopsin QM/MM models whose experimental absorption maxima (λ a max) have been measured. The results indicate that MC-PDFT may be employed to calculate λ a max values for this important class of photoresponsive proteins.

Understanding Colour Tuning Rules and Predicting Absorption Wavelengths of Microbial Rhodopsins by Data-Driven Machine-Learning Approach

The light-dependent ion-transport function of microbial rhodopsin has been widely used in optogenetics for optical control of neural activity. In order to increase the variety of rhodopsin proteins having a wide range of absorption wavelengths, the light absorption properties of various wild-type rhodopsins and their artificially mutated variants were investigated in the literature. Here, we demonstrate that a machine-learning-based (ML-based) data-driven approach is useful for understanding and predicting the light-absorption properties of microbial rhodopsin proteins. We constructed a database of 796 proteins consisting of microbial rhodopsin wildtypes and their variants. We then proposed an ML method that produces a statistical model describing the relationship between amino-acid sequences and absorption wavelengths and demonstrated that the fitted statistical model is useful for understanding colour tuning rules and predicting absorption wavelengths. By applying the ML method to the database, two residues that were not considered in previous studies are newly identified to be important to colour shift.