Chimeras of channelrhodopsin-1 and -2 from Chlamydomonas reinhardtii exhibit distinctive light-induced structural changes from channelrhodopsin-2

Twisting and Protonation of Retinal Chromophore Regulate Channel Gating of Channelrhodopsin C1C2

Channelrhodopsins (ChRs) are light-gated ion channels and central optogenetic tools that can control neuronal activity with high temporal resolution at the single-cell level. Although their application in optogenetics has rapidly progressed, it is unsolved how their channels open and close. ChRs transport ions through a series of interlocking elementary processes that occur over a broad time scale of subpicoseconds to seconds. During these processes, the retinal chromophore functions as a channel regulatory domain and transfers the optical input as local structural changes to the channel operating domain, the helices, leading to channel gating. Thus, the core question on channel gating dynamics is how the retinal chromophore structure changes throughout the photocycle and what rate-limits the kinetics. Here, we investigated the structural changes in the retinal chromophore of canonical ChR, C1C2, in all photointermediates using time-resolved resonance Raman spectroscopy. Moreover, to reveal the rate-limiting factors of the photocycle and channel gating, we measured the kinetic isotope effect of all photoreaction processes using laser flash photolysis and laser patch clamp, respectively. Spectroscopic and electrophysiological results provided the following understanding of the channel gating: the retinal chromophore highly twists upon the retinal Schiff base (RSB) deprotonation, causing the surrounding helices to move and open the channel. The ion-conducting pathway includes the RSB, where inflowing water mediates the proton to the deprotonated RSB. The twisting of the retinal chromophore relaxes upon the RSB reprotonation, which closes the channel. The RSB reprotonation rate-limits the channel closing.

Concerted primary proton transfer reactions in a thermophilic rhodopsin studied by time-resolved infrared spectroscopy at high temperature

The primary proton transfer reactions of thermophilic rhodopsin, which was first discovered in an extreme thermophile, Thermus thermophilus JL-18, were investigated using time-resolved Fourier transform infrared spectroscopy at various temperatures ranging from 298 to 343 K (25 to 70 °C) and proton transport activity analysis. The analyses were performed using counterion (D95E, D95N, D229E, and D229N) and proton donor mutants (E106D and E106Q) as well. First, the initial proton transfer from the protonated retinal Schiff base (PRSB) to D95 was identified. The temperature dependency showed that the proton transfer reaction in the intermediate states dramatically changed above 318 K (45 °C). In addition, the proton transfer reaction correlated well with the structural change from turn to β-strand in the protein moiety, suggesting that this step may be regulated by the rigidity of the loop region. We also elucidated that the proton transfer reaction from proton donor E106 to the retinal Schiff base occurred synchronously with the primary proton transfer from the PRSB to D95. Surprisingly, we discovered that the direction of proton transfer was regulated by the secondary counterion, D229. Comparative analysis of Gloeobacter rhodopsin from the mesophile, Gloeobacter violaceus, highlighted that the primary proton transfer reactions in thermophilic rhodopsin were optimized at high temperatures partly due to the specific turn to β-strand structural change. This was not observed in Gloeobacter rhodopsin and other related proteins such as bacteriorhodopsin.

Time-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin

Channelrhodopsins (ChRs) are microbial light-gated ion channels utilized in optogenetics to control neural activity with light . Light absorption causes retinal chromophore isomerization and subsequent protein conformational changes visualized as optically distinguished intermediates, coupled with channel opening and closing. However, the detailed molecular events underlying channel gating remain unknown. We performed time-resolved serial femtosecond crystallographic analyses of ChR by using an X-ray free electron laser, which revealed conformational changes following photoactivation. The isomerized retinal adopts a twisted conformation and shifts toward the putative internal proton donor residues, consequently inducing an outward shift of TM3, as well as a local deformation in TM7. These early conformational changes in the pore-forming helices should be the triggers that lead to opening of the ion conducting pore.

Optogenetics: A revolutionary approach for the study of depression

Depression is a severe and chronic mental disorder that affects millions of individuals worldwide. Symptoms include depressed mood, loss of interest, reduced motivation and suicidal thoughts. Even though findings from genetic, molecular and imaging studies have helped provide some clues regarding the mechanisms underlying depression-like behaviors, there are still many unanswered questions that need to be addressed. Optogenetics, a technique developed in the early 2000s, has proved effective in the study and treatment of depression and depression-like behaviors and has revolutionized already known experimental techniques. This technique employs light and genetic tools to either inhibit or excite specific neurons or pathways within the brain. In this review paper, an up-to-date understanding of the use of optogenetics in the study of depression-like behaviors is provided, along with suggestions for future research directions.

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.

The effect on ion channel of different protonation states of E90 in channelrhodopsin-2: a molecular dynamics simulation

Channelrhodopsin-2 (ChR2) is a cationic channel protein that has been extensively studied in optogenetics. The ion channel is opened via a series of proton transfers and H-bond changes during the photocycle but the detailed mechanism is still unknown. Molecular dynamics (MD) simulations with enhanced sampling were performed on the dark-adapted state (i.e., D470) and two photocycle intermediates (P₁⁵⁰⁰ and P₂³⁹⁰) to study the proton transfer path of the Schiff base and the subsequent conformational changes. The results suggest there are two possible proton transfer pathways from the Schiff base to proton acceptors (i.e., E123 or D253), depending on the protonation of E90. If E90 is protonated in the P₁⁵⁰⁰ state, the proton on the Schiff base will transfer to E123. The polyene chain of 13-cis retinal tilts and opens the channel that detours the blocking central gate (CG) and forms a narrow channel through the transmembrane helices (TM) 2, 3, 6 and 7. In contrast, if E90 deprotonates after retinal isomerization, the primary proton acceptor is D253, and an almost-open channel through TM1, 2, 3 and 7 is generated. The channel diameter is very close to the experimental value. The potential mean force (PMF) suggests that the free energy is extremely low for ions passing through this channel.

With E90 protonated, the proton acceptor of RSBH⁺ is E123 with a narrow channel along TM3; while with E90 deprotonated, proton transfer from RSBH⁺ to D253 generates an approximately open channel along TM2.

Gate-keeper of ion transport-a highly conserved helix-3 tryptophan in a channelrhodopsin chimera, C1C2/ChRWR

Microbial rhodopsin is a large family of membrane proteins having seven transmembrane helices (TM1-7) with an all-trans retinal (ATR) chromophore that is covalently bound to Lys in the TM7. The Trp residue in the middle of TM3, which is homologous to W86 of bacteriorhodopsin (BR), is highly conserved among microbial rhodopsins with various light-driven functions. However, the significance of this Trp for the ion transport function of microbial rhodopsins has long remained unknown. Here, we replaced the W163 (BR W86 counterpart) of a channelrhodopsin (ChR), C1C2/ChRWR, which is a chimera between ChR1 and 2, with a smaller aromatic residue, Phe to verify its role in the ion transport. Under whole-cell patch clamp recordings from the ND7/23 cells that were transfected with the DNA plasmid coding human codon optimized C1C2/ChRWR (hWR) or its W163F mutant (hWR-W163F), the photocurrents were evoked by a pulsatile light at 475 nm. The ion-transporting activity of hWR was strongly altered by the W163F mutation in 3 points: (1) the H⁺ leak at positive membrane potential (V_m) and its light-adaptation, (2) the attenuation of cation channel activity and (3) the manifestation of outward H⁺ pump activity. All of these results strongly suggest that W163 has a role in stabilizing the structure involved in the gating-on and -off of the cation channel, the role of “gate keeper”. We can attribute the attenuation of cation channel activity to the incomplete gating-on and the H⁺ leak to the incomplete gating-off.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Regulation of Photocycle Kinetics of Photoactive Yellow Protein by Modulating Flexibility of the β-Turn

The role of the significant flexibility of the β-turn in photoactive yellow protein (PYP) due to Gly115 was studied. G115A and G115P mutations were observed to accelerate the photocycle and shift the equilibrium between the late photocycle intermediate (pB) and its precursor (pR) toward pR. Thermodynamic analysis of dark-state recovery from pB demonstrated that the transition state (pB^⧧) has a negative change in transition heat capacity, suggesting that an exposed hydrophobic surface of pB is buried in pB^⧧. Fourier transform infrared spectroscopy showed that the structural ensemble of pB is populated by the compact structure in G115P. Taken together, the rigid structure induced by mutation of Gly115 facilitates its transition to pB^⧧, which adopts a substantially more compact structure as opposed to the ensemble-averaged structure of pB. The photocycle kinetics of PYP may be fine-tuned by modulating the flexibility of the 115 loop to activate an appropriate number of transducer proteins.

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Spectroscopic study of the transmembrane domain of a rhodopsin-phosphodiesterase fusion protein from a unicellular eukaryote

The choanoflagellate Salpingoeca rosetta contains a chimeric rhodopsin protein composed of an N-terminal rhodopsin (Rh) domain and a C-terminal cyclic nucleotide phosphodiesterase (PDE) domain. The Rh-PDE enzyme light-dependently decreases the concentrations of cyclic nucleotides such as cGMP and cAMP. Photoexcitation of purified full-length Rh-PDE yields an "M" intermediate with a deprotonated Schiff base, and its recovery is much faster than that of the enzyme domain. To gain structural and mechanistic insights into the Rh domain, here we expressed and purified the transmembrane domain of Rh-PDE, Rh-PDE(TMD), and analyzed it with transient absorption, light-induced difference UV-visible, and FTIR spectroscopy methods. These analyses revealed that the "K" intermediate forms within 0.005 ms and converts into the M intermediate with a time constant of 4 ms, with the latter returning to the original state within 4 s. FTIR spectroscopy revealed that all-trans to 13-cis photoisomerization occurs as the primary event during which chromophore distortion is located at the middle of the polyene chain, allowing the Schiff base to form a stronger hydrogen bond. We also noted that the peptide backbone of the α-helix becomes deformed upon M intermediate formation. Results from site-directed mutagenesis suggested that Glu-164 is protonated and that Asp-292 acts as the only Schiff base counterion in Rh-PDE. A strong reduction of enzymatic activity in a D292N variant, but not in an E164Q variant, indicated an important catalytic role of the negative charge at Asp-292. Our findings provide further mechanistic insights into rhodopsin-mediated, light-dependent regulation of second-messenger levels in eukaryotic microbes.

An Atomistic Model of a Precursor State of Light-Induced Channel Opening of Channelrhodopsin

Channelrhodopsins (ChRs) are microbial light-gated ion channels with a retinal chromophore and are widely utilized in optogenetics to precisely control neuronal activity with light. Despite increasing understanding of their structures and photoactivation kinetics, the atomistic mechanism of light gating and ion conduction remains elusive. Here, we present an atomic structural model of a chimeric ChR in a precursor state of the channel opening determined by an accurate hybrid molecular simulation technique and a statistical theory of internal water distribution. The photoactivated structure features extensive tilt of the chromophore accompanied by redistribution of water molecules in its binding pocket, which is absent in previously known photoactivated structures of analogous photoreceptors, and widely agrees with structural and spectroscopic experimental evidence of ChRs. The atomistic model manifests a photoactivated ion-conduction pathway that is markedly different from a previously proposed one and successfully explains experimentally observed mutagenic effects on key channel properties.

Structural insights into ion conduction by channelrhodopsin 2

The light-gated ion channel channelrhodopsin 2 (ChR2) from Chlamydomonas reinhardtii is a major optogenetic tool. Photon absorption starts a well-characterized photocycle, but the structural basis for the regulation of channel opening remains unclear. We present high-resolution structures of ChR2 and the C128T mutant, which has a markedly increased open-state lifetime. The structure reveals two cavities on the intracellular side and two cavities on the extracellular side. They are connected by extended hydrogen-bonding networks involving water molecules and side-chain residues. Central is the retinal Schiff base that controls and synchronizes three gates that separate the cavities. Separate from this network is the DC gate that comprises a water-mediated bond between C128 and D156 and interacts directly with the retinal Schiff base. Comparison with the C128T structure reveals a direct connection of the DC gate to the central gate and suggests how the gating mechanism is affected by subtle tuning of the Schiff base's interactions.

Retinal isomerization and water-pore formation in channelrhodopsin-2

Channelrhodopsin-2 (ChR2) is a light-sensitive ion channel widely used in optogenetics. Photoactivation triggers a trans-to-cis isomerization of a covalently bound retinal. Ensuing conformational changes open a cation-selective channel. We explore the structural dynamics in the early photocycle leading to channel opening by classical (MM) and quantum mechanical (QM) molecular simulations. With QM/MM simulations, we generated a protein-adapted force field for the retinal chromophore, which we validated against absorption spectra. In a 4-µs MM simulation of a dark-adapted ChR2 dimer, water entered the vestibules of the closed channel. Retinal all-trans to 13-cis isomerization, simulated with metadynamics, triggered a major restructuring of the charge cluster forming the channel gate. On a microsecond time scale, water penetrated the gate to form a membrane-spanning preopen pore between helices H1, H2, H3, and H7. This influx of water into an ion-impermeable preopen pore is consistent with time-resolved infrared spectroscopy and electrophysiology experiments. In the retinal 13-cis state, D253 emerged as the proton acceptor of the Schiff base. Upon proton transfer from the Schiff base to D253, modeled by QM/MM simulations, we obtained an early-M/P2390-like intermediate. Rapid rotation of the unprotonated Schiff base toward the cytosolic side effectively prevents its reprotonation from the extracellular side. From MM and QM simulations, we gained detailed insight into the mechanism of ChR2 photoactivation and early events in pore formation. By rearranging the network of charges and hydrogen bonds forming the gate, water emerges as a key player in light-driven ChR2 channel opening.

Spatiotemporal Control of GPR37 Signaling and Its Behavioral Effects by Optogenetics

Despite the progress in deorphanization of G Protein-Coupled Receptors (GPCRs), ≈100 GPCRs are still classified as orphan receptors without identified endogenous ligands and with unknown physiological functions. The lack of endogenous ligands triggering GPCR signaling has hampered the study of orphan GPCR functions. Using GPR37 as an example, we provide here the first demonstration of the channelrhodopsin 2 (ChR2)-GPCR approach to bypass the endogenous ligand and selectively activate the orphan GPCR signal by optogenetics. Inspired by the opto-XR approach, we designed the ChR2-GPR37 chimera, in which the corresponding parts of GPR37 replaced the intracellular portions of ChR2. We showed that optogenetic activation of ChR2/opto-GPR37 elicited specific GPR37 signaling, as evidenced by reduced cAMP level, enhanced ERK phosphorylation and increased motor activity, confirming the specificity of opto-GPR37 signaling. Besides, optogenetic activation of opto-GPR37 uncovered novel aspects of GPR37 signaling (such as IP-3 signaling) and anxiety-related behavior. Optogenetic activation of opto-GPR37 permits the causal analysis of GPR37 activity in the defined cells and behavioral responses of freely moving animals. Importantly, given the evolutionarily conserved seven-helix transmembrane structures of ChR2 and orphan GPCRs, we propose that opto-GPR37 approach can be readily applied to other orphan GPCRs for their deorphanization in freely moving animals.

Structural properties determining low K+ affinity of the selectivity filter in the TWIK1 K+ channel

Canonical K⁺ channels are tetrameric and highly K⁺-selective, whereas two-pore-domain K⁺ (K2P) channels form dimers, but with a similar pore architecture. A two-pore-domain potassium channel TWIK1 (KCNK1 or K2P1) allows permeation of Na⁺ and other monovalent ions, resulting mainly from the presence of Thr-118 in the P1 domain. However, the mechanistic basis for this reduced selectivity is unclear. Using ion-exchange-induced difference IR spectroscopy, we analyzed WT TWIK1 and T118I (highly K⁺-selective) and L228F (substitution in the P2 domain) TWIK1 variants and found that in the presence of K⁺ ions, WT and both variants exhibit an amide-I band at 1680 cm^-1 This band corresponds to interactions of the backbone carbonyls in the selectivity filter with K⁺, a feature very similar to that of the canonical K⁺ channel KcsA. Computational analysis indicated that the relatively high frequency for the amide-I band is well explained by impairment of hydrogen bond formation with water molecules. Moreover, concentration-dependent spectral changes indicated that the K⁺ affinity of the WT selectivity filter was much lower than those of the variants. Furthermore, only the variants displayed a higher frequency shift of the 1680-cm^-1 band upon changes from K⁺ to Rb⁺ or Cs⁺ conditions. High-speed atomic force microscopy disclosed that TWIK1's surface morphology largely does not change in K⁺ and Na⁺ solutions. Our results reveal the local conformational changes of the TWIK1 selectivity filter and suggest that the amide-I bands may be useful "molecular fingerprints" for assessing the properties of other K⁺ channels.

Proton transfer reactions in the red light-activatable channelrhodopsin variant ReaChR and their relevance for its function

Channelrhodopsins (ChRs) are light-gated ion channels widely used for activating selected cells in large cellular networks. ChR variants with a red-shifted absorption maximum, such as the modified Volvox carteri ChR1 red-activatable channelrhodopsin ("ReaChR," λmax = 527 nm), are of particular interest because longer wavelengths allow optical excitation of cells in deeper layers of organic tissue. In all ChRs investigated so far, proton transfer reactions and hydrogen bond changes are crucial for the formation of the ion-conducting pore and the selectivity for protons versus cations, such as Na+, K+, and Ca2+ (1). By using a combination of electrophysiological measurements and UV-visible and FTIR spectroscopy, we characterized the proton transfer events in the photocycle of ReaChR and describe their relevance for its function. 1) The central gate residue Glu130 (Glu90 in Chlamydomonas reinhardtii (Cr) ChR2) (i) undergoes a hydrogen bond change in D → K transition and (ii) deprotonates in K → M transition. Its negative charge in the open state is decisive for proton selectivity. 2) The counter-ion Asp293 (Asp253 in CrChR2) receives the retinal Schiff base proton during M-state formation. Starting from M, a photocycle branching occurs involving (i) a direct M → D transition and (ii) formation of late photointermediates N and O. 3) The DC pair residue Asp196 (Asp156 in CrChR2) deprotonates in N → O transition. Interestingly, the D196N mutation increases 15-syn-retinal at the expense of 15-anti, which is the predominant isomer in the wild type, and abolishes the peak current in electrophysiological measurements. This suggests that the peak current is formed by 15-anti species, whereas 15-syn species contribute only to the stationary current.

Molecular properties of a DTD channelrhodopsin from Guillardia theta

Microbial rhodopsins are membrane proteins found widely in archaea, eubacteria and eukaryotes (fungal and algal species). They have various functions, such as light-driven ion pumps, light-gated ion channels, light sensors and light-activated enzymes. A light-driven proton pump bacteriorhodopsin (BR) contains a DTD motif at positions 85, 89, and 96, which is unique to archaeal proton pumps. Recently, channelrhodopsins (ChRs) containing the DTD motif, whose sequential identity is ~20% similar to BR and to cation ChRs in Chlamydomonas reinhardtii (CrCCRs), were found. While extensive studies on ChRs have been performed with CrCCR2, the molecular properties of DTD ChRs remain an intrigue. In this paper, we studied a DTD rhodopsin from G. theta (GtCCR4) using electrophysiological measurements, flash photolysis, and low-temperature difference FTIR spectroscopy. Electrophysiological measurements clearly showed that GtCCR4 functions as a light-gated cation channel, similar to other G. theta DTD ChRs (GtCCR1-3). Light-driven proton pump activity was also suggested for GtCCR4. Both electrophysiological and flash photolysis experiments showed that channel closing occurs upon reprotonation of the Schiff base, suggesting that the dynamics of retinal and channels are tightly coupled in GtCCR4. From Fourier transform infrared (FTIR) spectroscopy at 77 K, we found that the primary reaction is an all-trans to a 13-cis photoisomerization, like other microbial rhodopsins, although perturbations in the secondary structure were much smaller in GtCCR4 than in CrCCR2.

Kinetic characteristics of chimeric channelrhodopsins implicate the molecular identity involved in desensitization

Channelrhodopsin (ChR)-1 and ChR2 were the first-identified members of ChRs which are a growing subfamily of microbial-type rhodopsins. Light absorption drives the generation of a photocurrent in cell membranes expressing ChR2. However, the photocurrent amplitude attenuates and becomes steady-state during prolonged irradiation. This process, called desensitization or inactivation, has been attributed to the accumulation of intermediates less conductive to cations. Here we provided evidence that the dark-adapted (DA) photocurrent before desensitization is kinetically different from the light-adapted (LA) one after desensitization, that is, the deceleration of both basal-to-conductive and conductive-to-basal transitions. When the kinetics were compared between the DA and LA photocurrents for the ChR1/2 chimeras, the transmembrane helices, TM1 and TM2, were the determinants of both basal-to-conductive and conductive-to-basal transitions, whereas TM4 may contribute to the basal-to-conductive transitions and TM5 may contribute to the conductive-to-basal transitions, respectively. The fact that the desensitization-dependent decrease of the basal-to-conductive and conductive-to-basal transitions was facilitated by the TM1 exchange from ChR2 to ChR1 and reversed by the further TM2 exchange suggests that the conformation change for the channel gating is predominantly regulated by the interaction between TM1 and TM2. Although the exchange of TM1 from ChR2 to ChR1 showed no obvious influence on the spectral sensitivity, this exchange significantly induced the desensitization-dependent blue shift. Therefore, the TM1 and 2 are the main structures involved in two features of the desensitization, the stabilization of protein conformation and the charge distribution around the retinal-Schiff base (RSB⁺).

Resonance Raman Study of an Anion Channelrhodopsin: Effects of Mutations near the Retinylidene Schiff Base

Optogenetics relies on the expression of specific microbial rhodopsins in the neuronal plasma membrane. Most notably, this includes channelrhodopsins, which when heterologously expressed in neurons function as light-gated cation channels. Recently, a new class of microbial rhodopsins, termed anion channel rhodopsins (ACRs), has been discovered. These proteins function as efficient light-activated channels strictly selective for anions. They exclude the flow of protons and other cations and cause hyperpolarization of the membrane potential in neurons by allowing the inward flow of chloride ions. In this study, confocal near-infrared resonance Raman spectroscopy (RRS) along with hydrogen/deuterium exchange, retinal analogue substitution, and site-directed mutagenesis were used to study the retinal structure as well as its interactions with the protein in the unphotolyzed state of an ACR from Guillardia theta (GtACR1). These measurements reveal that (i) the retinal chromophore exists as an all-trans configuration with a protonated Schiff base (PSB) very similar to that of bacteriorhodopsin (BR), (ii) the chromophore RRS spectrum is insensitive to changes in pH from 3 to 11, whereas above this pH the Schiff base (SB) is deprotonated, (iii) when Ser97, the homologue to Asp85 in BR, is replaced with a Glu, it remains in a neutral form (i.e., as a carboxylic acid) but is deprotonated at higher pH to form a blue-shifted species, (iv) Asp234, the homologue of the protonated retinylidene SB counterion Asp212 in BR, does not serve as the primary counteranion for the protonated SB, and (v) substitution of Glu68 with an Gln increases the pH at which SB deprotonation is observed. These results suggest that Glu68 and Asp234 located near the SB exist in a neutral state in unphotolyzed GtACR1 and indicate that other unidentified negative charges stabilize the protonated state of the GtACR1 SB.

Polarizable Embedded RI-CC2 Method for Two-Photon Absorption Calculations

We present a novel polarizable embedded resolution-of-identity coupled cluster singles and approximate doubles (PERI-CC2) method for calculation of two-photon absorption (TPA) spectra of large molecular systems. The method was benchmarked for three types of systems: a water-solvated molecule of formamide, a uracil molecule in aqueous solution, and a set of mutants of the channelrhodopsin (ChR) protein. The first test case shows that the PERI-CC2 method is in excellent agreement with the PE-CC2 method and in good agreement with the PE-CCSD method. The uracil test case indicates that the effects of hydrogen bonding on the TPA of a chromophore with the nearest environment is well-described with the PERI-CC2 method. Finally, the ChR calculation shows that the PERI-CC2 method is well-suited and efficient for calculations on proteins with medium-sized chromophores.

Kinetic and vibrational isotope effects of proton transfer reactions in channelrhodopsin-2

Channelrhodopsins (ChRs) are light-gated cation channels. After blue-light excitation, the protein undergoes a photocycle with different intermediates. Here, we have recorded transient absorbance changes of ChR2 from Chlamydomonas reinhardtii in the visible and infrared regions with nanosecond time resolution, the latter being accomplished using tunable quantum cascade lasers. Because proton transfer reactions play a key role in channel gating, we determined vibrational as well as kinetic isotope effects (VIEs and KIEs) of carboxylic groups of various key aspartic and glutamic acid residues by monitoring their C=O stretching vibrations in H2O and in D2O. D156 exhibits a substantial KIE (>2) in its deprotonation and reprotonation, which substantiates its role as the internal proton donor to the retinal Schiff base. The unusual VIE of D156, upshifted from 1736 cm(-1) to 1738 cm(-1) in D2O, was scrutinized by studying the D156E variant. The C=O stretch of E156 shifted down by 8 cm(-1) in D2O, providing evidence for the accessibility of the carboxylic group. The C=O stretching band of E90 exhibits a VIE of 9 cm(-1) and a KIE of ∼2 for the de- and the reprotonation reactions during the lifetime of the late desensitized state. The KIE of 1 determined in the time range from 20 ns to 5 ms is incompatible with early deprotonation of E90.

Molecular Dynamics of Channelrhodopsin at the Early Stages of Channel Opening

Channelrhodopsin (ChR) is a light-gated cation channel that responds to blue light. Since ChR can be readily expressed in specific neurons to precisely control their activities by light, it has become a powerful tool in neuroscience. Although the recently solved crystal structure of a chimeric ChR, C1C2, provided the structural basis for ChR, our understanding of the molecular mechanism of ChR still remains limited. Here we performed electrophysiological analyses and all-atom molecular dynamics (MD) simulations, to investigate the importance of the intracellular and central constrictions of the ion conducting pore observed in the crystal structure of C1C2. Our electrophysiological analysis revealed that two glutamate residues, Glu122 and Glu129, in the intracellular and central constrictions, respectively, should be deprotonated in the photocycle. The simulation results suggested that the deprotonation of Glu129 in the central constriction leads to ion leakage in the ground state, and implied that the protonation of Glu129 is important for preventing ion leakage in the ground state. Moreover, we modeled the 13-cis retinal bound; i.e., activated C1C2, and performed MD simulations to investigate the conformational changes in the early stage of the photocycle. Our simulations suggested that retinal photoisomerization induces the conformational change toward channel opening, including the movements of TM6, TM7 and TM2. These insights into the dynamics of the ground states and the early photocycle stages enhance our understanding of the channel function of ChR.

Spectral properties and isomerisation path of retinal in C1C2 channelrhodopsin

Computed torsion profiles along the reactive coordinate in S₁reveal a two-path deactivation mechanism for retinal in C1C2 channelrhodopsin.