Metadynamics simulations reveal mechanisms of Na+ and Ca2+ transport in two open states of the channelrhodopsin chimera, C1C2.

Cation conducting channelrhodopsins (ChRs) are a popular tool used in optogenetics to control the activity of excitable cells and tissues using light. ChRs with altered ion selectivity are in high demand for use in different cell types and for other specialized applications. However, a detailed mechanism of ion permeation in ChRs is not fully resolved. Here, we use complementary experimental and computational methods to uncover the mechanisms of cation transport and valence selectivity through the channelrhodopsin chimera, C1C2, in the high- and low-conducting open states. Electrophysiology measurements identified a single-residue substitution within the central gate, N297D, that increased Ca2+ permeability vs. Na+ by nearly two-fold at peak current, but less so at stationary current. We then developed molecular models of dimeric wild-type C1C2 and N297D mutant channels in both open states and calculated the PMF profiles for Na+ and Ca2+ permeation through each protein using well-tempered/multiple-walker metadynamics. Results of these studies agree well with experimental measurements and demonstrate that the pore entrance on the extracellular side differs from original predictions and is actually located in a gap between helices I and II. Cation transport occurs via a relay mechanism where cations are passed between flexible carboxylate sidechains lining the full length of the pore by sidechain swinging, like a monkey swinging on vines. In the mutant channel, residue D297 enhances Ca2+ permeability by mediating the handoff between the central and cytosolic binding sites via direct coordination and sidechain swinging. We also found that altered cation binding affinities at both the extracellular entrance and central binding sites underly the distinct transport properties of the low-conducting open state. This work significantly advances our understanding of ion selectivity and permeation in cation channelrhodopsins and provides the insights needed for successful development of new ion-selective optogenetic tools.

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.

Channelrhodopsins with distinct chromophores and binding patterns.

Channelrhodopsins are popular optogenetic tools in neuroscience, but remain poorly understood mechanistically. Here we report the cryo-EM structures of channelrhodopsin-2 (ChR2) from Chlamydomonas reinhardtii and H. catenoides kalium channelrhodopsin (KCR1). We show that ChR2 recruits an endogenous N-retinylidene-PE-like molecule to a previously unidentified lateral retinal binding pocket, exhibiting a reduced light response in HEK293 cells. In contrast, H. catenoides kalium channelrhodopsin (KCR1) binds an endogenous retinal in its canonical retinal binding pocket under identical condition. However, exogenous ATR reduces the photocurrent magnitude of wild type KCR1 and also inhibits its leaky mutant C110T. Our results uncover diverse retinal chromophores with distinct binding patterns for channelrhodopsins in mammalian cells, which may further inspire next generation optogenetics for complex tasks such as cell fate control.

Unique hydrogen-bonding network in a viral channelrhodopsin.

Channelrhodopsins (CRs) are used as key tools in optogenetics, and novel CRs, either found from nature or engineered by mutation, have greatly contributed to the development of optogenetics. Recently CRs were discovered from viruses, and crystal structure of a viral CR, OLPVR1, reported a very similar water-containing hydrogen-bonding network near the retinal Schiff base to that of a light-driven proton-pump bacteriorhodopsin (BR). In both OLPVR1 and BR, nearly planar pentagonal cluster structures are comprised of five oxygen atoms, three oxygens from water molecules and two oxygens from the Schiff base counterions. The planar pentagonal cluster stabilizes a quadrupole, two positive charges at the Schiff base and an arginine, and two negative charges at the counterions, and thus plays important roles in light-gated channel function of OLPVR1 and light-driven proton pump function of BR. Despite similar pentagonal cluster structures, present FTIR analysis revealed different hydrogen-bonding networks between OLPVR1 and BR. The hydrogen bond between the protonated Schiff base and a water is stronger in OLPVR1 than in BR, and internal water molecules donate hydrogen bonds much weaker in OLPVR1 than in BR. In OLPVR1, the bridged water molecule between the Schiff base and counterions forms hydrogen bonds to D76 and D200 equally, while the hydrogen-bonding interaction is much stronger to D85 than to D212 in BR. The present interpretation is supported by the mutation results, where D76 and D200 equally work as the Schiff base counterions in OLPVR1, but D85 is the primary counterion in BR. This work reports highly sensitive hydrogen-bonding network in the Schiff base region, which would be closely related to each function through light-induced alterations of the network.

Parallel photocycle kinetic model of anion channelrhodopsin GtACR1 function.

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

Molecular Mechanisms behind Circular Dichroism Spectral Variations between Channelrhodopsin and Heliorhodopsin Dimers.

Channelrhodopsin (ChR) and heliorhodopsin (HeR) are microbial rhodopsins with similar structures but different circular dichroism (CD) spectra: ChR shows biphasic negative and positive bands, whereas HeR shows a single positive band. We explored the physicochemical factors underlying these differences through computational methods. Using the exciton model based on first-principles computations, we obtained the CD spectra of ChR and HeR. The obtained spectra indicate that the protein dimer structures and the quantum mechanical treatment of the retinal chromophore and its interacting amino acids are crucial for accurately reproducing the experimental spectra. Further calculations revealed that the sign of the excitonic coupling was opposite between the ChR and HeR dimers, which was attributed to the contrasting second term of the orientation factor between the two retinal chromophores. These findings demonstrate that slight variations in the intermolecular orientation of the two chromophores can result in significant differences in the CD spectral shape.

Control of polymers' amorphous-crystalline transition enables miniaturization and multifunctional integration for hydrogel bioelectronics.

Soft bioelectronic devices exhibit motion-adaptive properties for neural interfaces to investigate complex neural circuits. Here, we develop a fabrication approach through the control of metamorphic polymers' amorphous-crystalline transition to miniaturize and integrate multiple components into hydrogel bioelectronics. We attain an about 80% diameter reduction in chemically cross-linked polyvinyl alcohol hydrogel fibers in a fully hydrated state. This strategy allows regulation of hydrogel properties, including refractive index (1.37-1.40 at 480 nm), light transmission (>96%), stretchability (139-169%), bending stiffness (4.6 ± 1.4 N/m), and elastic modulus (2.8-9.3 MPa). To exploit the applications, we apply step-index hydrogel optical probes in the mouse ventral tegmental area, coupled with fiber photometry recordings and social behavioral assays. Additionally, we fabricate carbon nanotubes-PVA hydrogel microelectrodes by incorporating conductive nanomaterials in hydrogel for spontaneous neural activities recording. We enable simultaneous optogenetic stimulation and electrophysiological recordings of light-triggered neural activities in Channelrhodopsin-2 transgenic mice.

Channel Gating in Kalium Channelrhodopsin Slow Mutants.

Kalium channelrhodopsin 1 from Hyphochytrium catenoides (HcKCR1) is the first discovered natural light-gated ion channel that shows higher selectivity to K+ than to Na+ and therefore is used to silence neurons with light (optogenetics). Replacement of the conserved cysteine residue in the transmembrane helix 3 (Cys110) with alanine or threonine results in a >1,000-fold decrease in the channel closing rate. The phenotype of the corresponding mutants in channelrhodopsin 2 is attributed to breaking of a specific interhelical hydrogen bond (the "DC gate"). Unlike CrChR2 and other ChRs with long distance "DC gates", the HcKCR1 structure does not reveal any hydrogen bonding partners to Cys110, indicating that the mutant phenotype is likely caused by disruption of direct interaction between this residue and the chromophore. In HcKCR1_C110A, fast photochemical conversions corresponding to channel gating were followed by dramatically slower absorption changes. Full recovery of the unphotolyzed state in HcKCR1_C110A was extremely slow with two time constants 5.2 and 70 min. Analysis of the light-minus-dark difference spectra during these slow processes revealed accumulation of at least four spectrally distinct blue light-absorbing photocycle intermediates, L, M1 and M2, and a UV light-absorbing form, typical of bacteriorhodopsin-like channelrhodopsins from cryptophytes. Our results contribute to better understanding of the mechanistic links between the chromophore photochemistry and channel conductance, and provide the basis for using HcKCR1_C110A as an optogenetic tool.

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.

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.

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.

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.

The crystal structure of bromide-bound GtACR1 reveals a pre-activated state in the transmembrane anion tunnel.

The crystal structure of the light-gated anion channel GtACR1 reported in our previous Research Article (Li et al., 2019) revealed a continuous tunnel traversing the protein from extracellular to intracellular pores. We proposed the tunnel as the conductance channel closed by three constrictions: C1 in the extracellular half, mid-membrane C2 containing the photoactive site, and C3 on the cytoplasmic side. Reported here, the crystal structure of bromide-bound GtACR1 reveals structural changes that relax the C1 and C3 constrictions, including a novel salt-bridge switch mechanism involving C1 and the photoactive site. These findings indicate that substrate binding induces a transition from an inactivated state to a pre-activated state in the dark that facilitates channel opening by reducing free energy in the tunnel constrictions. The results provide direct evidence that the tunnel is the closed form of the channel of GtACR1 and shed light on the light-gated channel activation mechanism.

The Desensitized Channelrhodopsin-2 Photointermediate Contains 13 -cis, 15 -syn†Retinal Schiff Base.

Channelrhodopsin-2 (ChR2) is a light-gated cation channel and was used to lay the foundations of optogenetics. Its dark state X-ray structure has been determined in 2017 for the wild-type, which is the prototype for all other ChR variants. However, the mechanistic understanding of the channel function is still incomplete in terms of structural changes after photon absorption by the retinal chromophore and in the framework of functional models. Hence, detailed information needs to be collected on the dark state as well as on the different photointermediates. For ChR2 detailed knowledge on the chromophore configuration in the different states is still missing and a consensus has not been achieved. Using DNP-enhanced solid-state MAS NMR spectroscopy on proteoliposome samples, we unambiguously determined the chromophore configuration in the desensitized state, and we show that this state occurs towards the end of the photocycle.

Electrostatic Control of Photoisomerization in Channelrhodopsin 2.

Channelrhodopsin 2 (ChR2) is the most commonly used tool in optogenetics. Because of its faster photocycle compared to wild-type (WT) ChR2, the E123T mutant of ChR2 is a useful optogenetic tool when fast neuronal stimulation is needed. Interestingly, in spite of its faster photocycle, the initial step of the photocycle in E123T (photoisomerization of retinal protonated Schiff base or RPSB) was found experimentally to be much slower than that of WT ChR2. The E123T mutant replaces the negatively charged E123 residue with a neutral T123 residue, perturbing the electric field around the RPSB. Understanding the RPSB photoisomerization mechanism in ChR2 mutants will provide molecular-level insights into how ChR2 photochemical reactivity can be controlled, which will lay the foundation for improving the design of optogenetic tools. In this work, we combine ab initio nonadiabatic dynamics simulation, excited state free energy calculation, and reaction path search to comprehensively characterize the RPSB photoisomerization mechanism in the E123T mutant of ChR2. Our simulation agrees with previous experiments in predicting a red-shifted absorption spectrum and significant slowdown of photoisomerization in the E123T mutant. Interestingly, our simulations predict similar photoisomerization quantum yields for the mutant and WT despite the differences in excited-state lifetime and absorption maximum. Upon mutation, the neutralization of the negative charge on the E123 residue increases the isomerization barrier, alters the reaction pathway, and changes the relative stability of two fluorescent states. Our findings provide new insight into the intricate role of the electrostatic environment on the RPSB photoisomerization mechanism in microbial rhodopsins.

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.

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.

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.

Molecular Dynamics Simulations of Channelrhodopsin Chimera, C1C2.

Molecular dynamics (MD) simulations have been successfully used for modeling dynamic behavior of biologically relevant systems, such as ion channels in representative environments to decode protein structure-function relationships. Protocol presented here describes steps for generating input files and modeling a monomer of transmembrane cation channel, channelrhodopsin chimera (C1C2), in representative environment of 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) planar lipid bilayer, TIP3P water and ions (Na+ and Cl-) using molecular dynamics package NAMD, molecular graphics/analysis tool VMD, and other relevant tools. MD simulations of C1C2 were performed at 303.15 K and in constant particle number, isothermal-isobaric (NpT) ensemble. The results of modeling have helped understand how key interactions in the center of the C1C2 channel contribute to channel gating and subsequent solvent transport across the membrane.

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2.

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.