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.     

Platymonas subcordiformis Channelrhodopsin-2 (PsChR2) Function: II. RELATIONSHIP OF THE PHOTOCHEMICAL REACTION CYCLE TO CHANNEL CURRENTS*

Channelrhodopsins, such as the algal phototaxis receptor Platymonas subcordiformis channelrhodopsin-2 (PsChR2), are light-gated cation channels used as optogenetic tools for photocontrol of membrane potential in living cells. Channelrhodopsin (ChR)-mediated photocurrent responses are complex and poorly understood, exhibiting alterations in peak current amplitude, extents and kinetics of inactivation, and kinetics of the recovery of the prestimulus dark current that are sensitive to duration and frequency of photostimuli. From the analysis of time-resolved optical absorption data, presented in the accompanying article, we derived a two-cycle model that describes the photocycles of PsChR2. Here, we applied the model to evaluate the transient currents produced by PsChR2 expressed in HEK293 cells under both fast laser excitation and step-like continuous illumination. Interpretation of the photocurrents in terms of the photocycle kinetics indicates that the O states in both cycles are responsible for the channel current and fit the current transients under the different illumination regimes. The peak and plateau currents in response to a single light step, a train of light pulses, and a light step superimposed on a continuous light background observed for ChR2 proteins are explained in terms of contributions from the two parallel photocycles. The analysis shows that the peak current desensitization and recovery phenomena are inherent properties of the photocycles. The light dependence of desensitization is reproduced and explained by the time evolution of the concentration transients in response to step-like illumination. Our data show that photocycle kinetic parameters are sufficient to explain the complex dependence of photocurrent responses to photostimuli.     

Spectral characteristics of the photocycle of channelrhodopsin-2 and its implication for channel function

In 2003, channelrhodopsin-2 (ChR2) from Chlamydomonas reinhardtii was discovered to be a light-gated cation channel, and since that time the channel became an excellent tool to control by light neuronal cells in culture as well as in living animals with high temporal and spatial resolution in a noninvasive manner. However, little is known about the spectral properties and their relation to the channel function. We have expressed ChR2 in the yeast Pichia pastoris and purified the protein. Flash-photolysis data were combined with patch-clamp studies to elucidate the photocycle. The protein absorbs maximally at ∼ 480 nm before light excitation and shows flash-induced absorbance changes with at least two different photointermediates. Four relaxation processes can be extracted from the time course that we have analysed in a linear model for the photocycle leading to the kinetic intermediates P₁ to P₄. A short-lived photointermediate at 400 nm, suggesting a deprotonation of the retinal Schiff base, is followed by a red-shifted (520 nm) species with a millisecond lifetime. The first three kinetic intermediates in the photocycle, P₁ to P₃, are described mainly by the red-shifted 520-nm species. The 400-nm species contributes to a smaller extent to P₁ and P₂. The fourth one, P₄, is spectroscopically almost identical with the ground state and lasts into the seconds time region. We compared the spectroscopic data to current measurements under whole-cell patch-clamp conditions on HEK 293 cells. The lifetimes of the spectroscopically and electrophysiologically determined intermediates are in excellent agreement. The intermediates P₂ and P₃ (absorbing at 520 nm) are identified as the cation permeating states of the channel. Under stationary light, a modulation of the photocurrent by green light (540 nm) was observed. We conclude that the red-shifted spectral species represents the open channel state, and the thermal relaxation of this intermediate, the transition from P₃ to P₄, is coupled to channel closing.     

Characterization of a Highly Efficient Blue-shifted Channelrhodopsin from the Marine Alga Platymonas subcordiformis

Rhodopsin photosensors of phototactic algae act as light-gated cation channels when expressed in animal cells. These proteins (channelrhodopsins) are extensively used for millisecond scale photocontrol of cellular functions (optogenetics). We report characterization of PsChR, one of the phototaxis receptors in the alga Platymonas (Tetraselmis) subcordiformis. PsChR exhibited ∼3-fold higher unitary conductance and greater relative permeability for Na⁺ ions, as compared with the most frequently used channelrhodopsin-2 from Chlamydomonas reinhardtii (CrChR2). Photocurrents generated by PsChR in HEK293 cells showed lesser inactivation and faster peak recovery than those by CrChR2. Their maximal spectral sensitivity was at 445 nm, making PsChR the most blue-shifted channelrhodopsin so far identified. The λ_max of detergent-purified PsChR was 437 nm at neutral pH and exhibited red shifts (pK_a values at 6.6 and 3.8) upon acidification. The purified pigment undergoes a photocycle with a prominent red-shifted intermediate whose formation and decay kinetics match the kinetics of channel opening and closing. The rise and decay of an M-like intermediate prior to formation of this putative conductive state were faster than in CrChR2. PsChR mediated sufficient light-induced membrane depolarization in cultured hippocampal neurons to trigger reliable repetitive spiking at the upper threshold frequency of the neurons. At low frequencies spiking probability decreases less with PsChR than with CrChR2 because of the faster recovery of the former. Its blue-shifted absorption enables optogenetics at wavelengths even below 400 nm. A combination of characteristics makes PsChR important for further research on structure-function relationships in ChRs and potentially useful for optogenetics, especially for combinatorial applications when short wavelength excitation is required.     

Computational Optogenetics: Empirically-Derived Voltage- and Light-Sensitive Channelrhodopsin-2 Model

Channelrhodospin-2 (ChR2), a light-sensitive ion channel, and its variants have emerged as new excitatory optogenetic tools not only in neuroscience, but also in other areas, including cardiac electrophysiology. An accurate quantitative model of ChR2 is necessary for in silico prediction of the response to optical stimulation in realistic tissue/organ settings. Such a model can guide the rational design of new ion channel functionality tailored to different cell types/tissues. Focusing on one of the most widely used ChR2 mutants (H134R) with enhanced current, we collected a comprehensive experimental data set of the response of this ion channel to different irradiances and voltages, and used these data to develop a model of ChR2 with empirically-derived voltage- and irradiance- dependence, where parameters were fine-tuned via simulated annealing optimization. This ChR2 model offers: 1) accurate inward rectification in the current-voltage response across irradiances; 2) empirically-derived voltage- and light-dependent kinetics (activation, deactivation and recovery from inactivation); and 3) accurate amplitude and morphology of the response across voltage and irradiance settings. Temperature-scaling factors (Q₁₀) were derived and model kinetics was adjusted to physiological temperatures. Using optical action potential clamp, we experimentally validated model-predicted ChR2 behavior in guinea pig ventricular myocytes. The model was then incorporated in a variety of cardiac myocytes, including human ventricular, atrial and Purkinje cell models. We demonstrate the ability of ChR2 to trigger action potentials in human cardiomyocytes at relatively low light levels, as well as the differential response of these cells to light, with the Purkinje cells being most easily excitable and ventricular cells requiring the highest irradiance at all pulse durations. This new experimentally-validated ChR2 model will facilitate virtual experimentation in neural and cardiac optogenetics at the cell and organ level and provide guidance for the development of in vivo tools.     

Monitoring light-induced structural changes of Channelrhodopsin-2 by UV-visible and Fourier transform infrared spectroscopy

Channelrhodopsin-2 (ChR2) is a microbial type rhodopsin and a light-gated cation channel that controls phototaxis in Chlamydomonas. We expressed ChR2 in COS-cells, purified it, and subsequently investigated this unusual photoreceptor by flash photolysis and UV-visible and Fourier transform infrared difference spectroscopy. Several transient photoproducts of the wild type ChR2 were identified, and their kinetics and molecular properties were compared with those of the ChR2 mutant E90Q. Based on the spectroscopic data we developed a model of the photocycle comprising six distinguishable intermediates. This photocycle shows similarities to the photocycle of the ChR2-related Channelrhodopsin of Volvox but also displays significant differences. We show that molecular changes include retinal isomerization, changes in hydrogen bonding of carboxylic acids, and large alterations of the protein backbone structure. These alterations are stronger than those observed in the photocycle of other microbial rhodopsins like bacteriorhodopsin and are related to those occurring in animal rhodopsins. UV-visible and Fourier transform infrared difference spectroscopy revealed two late intermediates with different time constants of tau = 6 and 40 s that exist during the recovery of the dark state. The carboxylic side chain of Glu(90) is involved in the slow transition. The molecular changes during the ChR2 photocycle are discussed with respect to other members of the rhodopsin family.     

A Method for High Fidelity Optogenetic Control of Individual Pyramidal Neurons In vivo

Optogenetic methods have emerged as a powerful tool for elucidating neural circuit activity underlying a diverse set of behaviors across a broad range of species. Optogenetic tools of microbial origin consist of light-sensitive membrane proteins that are able to activate (e.g., channelrhodopsin-2, ChR2) or silence (e.g., halorhodopsin, NpHR) neural activity ingenetically-defined cell types over behaviorally-relevant timescales. We first demonstrate a simple approach for adeno-associated virus-mediated delivery of ChR2 and NpHR transgenes to the dorsal subiculum and prelimbic region of the prefrontal cortex in rat. Because ChR2 and NpHR are genetically targetable, we describe the use of this technology to control the electrical activity of specific populations of neurons (i.e., pyramidal neurons) embedded in heterogeneous tissue with high temporal precision. We describe herein the hardware, custom software user interface, and procedures that allow for simultaneous light delivery and electrical recording from transduced pyramidal neurons in an anesthetized in vivo preparation. These light-responsive tools provide the opportunity for identifying the causal contributions of different cell types to information processing and behavior.     

The chromophore structure of the long-lived intermediate of the C128T channelrhodopsin-2 variant

The photocycle of the light-activated channel, channelrhodopsin-2 C128T, has been studied by resonance Raman (RR) spectroscopy focussing on the intermediates P380 and P353 that constitute a side pathway in the recovery of the parent state. The P353 species displays a UV–vis absorption spectrum with a fine-structure reminiscent of the reduced-retro form of bacteriorhodopsin, whereas the respective RR spectra differ substantially. Instead, the RR spectra of the P380/P353 intermediate couple are closely related to that of a free retinal in the all-trans configuration. These findings imply that the parent state recovery via P380/P353 involves the transient hydrolysis and re-formation of the retinal–protein linkage.     

光敏感通道(Channelrhodopsin-2)——神经回路功能和神经系统疾病研究的新工具

光敏感通道(channelrhodopsin-2,ChR2)是一种受光脉冲控制的具有7次跨膜结构的非选择性阳离子通道蛋白,自1991年从莱茵衣藻中发现后被许多实验室所关注.依据ChR2可以快速形成光电流,使细胞发生去极化反应的电生理特性,ChR2已被广泛应用于神经系统的研究.与传统的神经系统研究方法如电生理技术、神经药理学方法相比,用光脉冲控制带有ChR2的神经元的活动,具有更高的空间选择性和特异性.ChR2作为光基因技术的核心组成部分,对神经科学是一个崭新的应用前景广泛的研究工具.近年来ChR2不仅应用于视觉、躯体感觉、听觉和嗅觉等多条感觉神经回路的形态和功能研究,还被应用于各种临床神经系统疾病的研究.本文总结了目前ChR2在神经系统中的研究进展,并对ChR2未来的应用前景作了进一步展望.     

Photoactivation of channelrhodopsin

Channelrhodopsins (ChRs) are light-gated ion channels that control photomovement of microalgae. In optogenetics, ChRs are widely applied for light-triggering action potentials in cells, tissues, and living animals, yet the spectral properties and photocycle of ChR remain obscure. In this study, we cloned a ChR from the colonial alga Volvox carteri, VChR. After electrophysiological characterization in Xenopus oocytes, VChR was expressed in COS-1 cells and purified. Time-resolved UV-visible spectroscopy revealed a pH-dependent equilibrium of two dark species, D(470)/D(480). Laser flashes converted both with tau approximately 200 mus into major photointermediates P(510)/P(530), which reverted back to the dark states with tau approximately 15-100 ms. Both intermediates were assigned to conducting states. Three early intermediates P(500)/P(515) and P(390) were detected on a ns to mus time scale. The spectroscopic and electrical data were unified in a photocycle model. The functional expression of VChR we report here paves the way toward a broader structure/function analysis of the recently identified class of light-gated ion channels.     

Redshift of the purple membrane absorption band and the deprotonation of tyrosine residues at high pH: Origin of the parallel photocycles of trans-bacteriorhodopsin

At high pH (> 8) the 570 nm absorption band of all-trans bacteriorhodopsin (bR) in purple membrane undergoes a small (1.5 nm) shift to longer wavelengths, which causes a maximal increase in absorption at 615 nm. The pK of the shift is 9.0 in the presence of 167 mM KCl, and its intrinsic pK is ~8.3. The red shift of the trans-bR absorption spectrum correlates with the appearance of the fast component in the light-induced L to M transition, and absorption increases at 238 and 297 nm which are apparently caused by the deprotonation of a tyrosine residue and red shift of the absorption of tryptophan residues. This suggests that the deprotonation of a tyrosine residue with an exceptionally low pK (pK_a ≈ 8.3) is responsible for the absorption shift of the chromophore band and fast M formation. The pH and salt dependent equilibrium between the two forms of bR, “neutral” and “alkaline,” bR ↔ bR_a, results in two parallel photocycles of trans-bR at high pH, differing in the rate of the L to M transition. In the pH range 10-11.8 deprotonation of two more tyrosine residues is observed with pK's ~ 10.3 and 11.3 (in 167 mM KCL). Two simple models discussing the role of the pH induced tyrosine deprotonation in the photocycle and proton pumping are presented.

It is suggested that the shifts of the absorption bands at high pH are due to the appearance of a negatively charged group inside the protein (tyrosinate) which causes electrochromic shifts of the chromophore and protein absorption bands due to the interaction with the dipole moments in the ground and excited states of bR (Stark effect). This effect gives evidence for a significant change in the dipole moment of the chromophore of bR upon excitation.

Under illumination alkaline bR forms, besides the usual photocycle intermediates, a long-lived species with absorption maximum at 500 nm (P500). P500 slowly converts into bR_a in the dark. Upon illumination P500 is transformed into an intermediate having an absorption maximum at 380 nm (P380). P380 can be reconverted to P500 by blue light illumination or by incubation in the dark.

     

Resonance Raman and FTIR spectroscopic characterization of the closed and open states of channelrhodopsin-1

Channelrhodopsin-1 from Chlamydomonas augustae (CaChR1) is a light-activated cation channel, which is a promising optogenetic tool. We show by resonance Raman spectroscopy and retinal extraction followed by high pressure liquid chromatography (HPLC) that the isomeric ratio of all-trans to 13-cis of solubilized channelrhodopsin-1 is with 70:30 identical to channelrhodopsin-2 from Chlamydomonas reinhardtii (CrChR2). Critical frequency shifts in the retinal vibrations are identified in the Raman spectrum upon transition to the open (conductive P₂³⁸⁰) state. Fourier transform infrared spectroscopy (FTIR) spectra indicate different structures of the open states in the two channelrhodopsins as reflected by the amide I bands and the protonation pattern of acidic amino acids.     

Platymonas subcordiformis Channelrhodopsin-2 Function: I. THE PHOTOCHEMICAL REACTION CYCLE*

The photocycle kinetics of Platymonas subcordiformis channelrhodopsin-2 (PsChR2), among the most highly efficient light-gated cation channels and the most blue-shifted channelrhodopsin, was studied by time-resolved absorption spectroscopy in the 340–650-nm range and in the 100-ns to 3-s time window. Global exponential fitting of the time dependence of spectral changes revealed six lifetimes: 0.60 μs, 5.3 μs, 170 μs, 1.4 ms, 6.7 ms, and 1.4 s. The sequential intermediates derived for a single unidirectional cycle scheme based on these lifetimes were found to contain mixtures of K, L, M, O, and P molecular states, named in analogy to photointermediates in the bacteriorhodopsin photocycle. The photochemistry is described by the superposition of two independent parallel photocycles. The analysis revealed that 30% of the photoexcited receptor molecules followed Cycle 1 through the K, M, O, and P states, whereas 70% followed Cycle 2 through the K, L, M, and O states. The recovered state, R, is spectrally close, but not identical, to the dark state on the seconds time scale. The two-cycle model of this high efficiency channelrhodopsin-2 (ChR) opens new perspectives in understanding the mechanism of channelrhodopsin function.     

Kinetics of proton release and uptake by channelrhodopsin-2

Electrophysiological experiments showed that the light-activated cation channel channelrhodopsin-2 (ChR2) pumps protons in the absence of a membrane potential. We determined here the kinetics of transient pH change using a water-soluble pH-indicator. It is shown that ChR2 released protons prior to uptake with a stoichiometry of 0.3 protons per ChR2. Comparison to the photocycle kinetics revealed that proton release and uptake match rise and decay of the P(3)(520) intermediate. As the P(3)(520) state also represents the conductive state of cation channeling, the concurrence of proton pumping and channel gating implies an intimate mechanistic link of the two functional modes. Studies on the E123T and S245E mutants show that these residues are not critically involved in proton translocation.     

Optimizing the spatial resolution of Channelrhodopsin-2 activation

Over the past few years, the light-gated cation channel Channelrhodopsin-2 (ChR2) has seen a remarkable diversity of applications in neuroscience. However, commonly used wide-field illumination provides poor spatial selectivity for cell stimulation. We explored the potential of focal laser illumination to map photocurrents of individual neurons in sparsely transfected hippocampal slice cultures. Interestingly, the best spatial resolution of photocurrent induction was obtained at the lowest laser power. By adjusting the light intensity to a neuron's spike threshold, we were able to trigger action potentials with a spatial selectivity of less than 30 microm. Experiments with dissociated hippocampal cells suggested that the main factor limiting the spatial resolution was ChR2 current density rather than scattering of the excitation light. We conclude that subcellular resolution can be achieved only in cells with a high ChR2 expression level and that future improved variants of ChR2 are likely to extend the spatial resolution of photocurrent induction to the level of single dendrites.     

Bimodal Activation of Different Neuron Classes with the Spectrally Red-Shifted Channelrhodopsin Chimera C1V1 in Caenorhabditis elegans

The C. elegans nervous system is particularly well suited for optogenetic analyses of circuit function: Essentially all connections have been mapped, and light can be directed at the neuron of interest in the freely moving, transparent animals, while behavior is observed. Thus, different nodes of a neuronal network can be probed for their role in controlling a particular behavior, using different optogenetic tools for photo-activation or –inhibition, which respond to different colors of light. As neurons may act in concert or in opposing ways to affect a behavior, one would further like to excite these neurons concomitantly, yet independent of each other. In addition to the blue-light activated Channelrhodopsin-2 (ChR2), spectrally red-shifted ChR variants have been explored recently. Here, we establish the green-light activated ChR chimera C1V1 (from Chlamydomonas and Volvox ChR1′s) for use in C. elegans. We surveyed a number of red-shifted ChRs, and found that C1V1-ET/ET (E122T; E162T) works most reliable in C. elegans, with 540–580 nm excitation, which leaves ChR2 silent. However, as C1V1-ET/ET is very light sensitive, it still becomes activated when ChR2 is stimulated, even at 400 nm. Thus, we generated a highly efficient blue ChR2, the H134R; T159C double mutant (ChR2-HR/TC). Both proteins can be used in the same animal, in different neurons, to independently control each cell type with light, enabling a further level of complexity in circuit analyses.     

Characterization of Engineered Channelrhodopsin Variants with Improved Properties and Kinetics

Channelrhodopsin 2 (ChR2), a light-activated nonselective cationic channel from Chlamydomonas reinhardtii, has become a useful tool to excite neurons into which it is transfected. The other ChR from Chlamydomonas, ChR1, has attracted less attention because of its proton-selective permeability. By making chimeras of the transmembrane domains of ChR1 and ChR2, combined with site-directed mutagenesis, we developed a ChR variant, named ChEF, that exhibits significantly less inactivation during persistent light stimulation. ChEF undergoes only 33% inactivation, compared with 77% for ChR2. Point mutation of Ile¹⁷⁰ of ChEF to Val (yielding “ChIEF”) accelerates the rate of channel closure while retaining reduced inactivation, leading to more consistent responses when stimulated above 25 Hz in both HEK293 cells and cultured hippocampal neurons. In addition, these variants have altered spectral responses, light sensitivity, and channel selectivity. ChEF and ChIEF allow more precise temporal control of depolarization, and can induce action potential trains that more closely resemble natural spiking patterns.     

Chimeras of Channelrhodopsin-1 and -2 from Chlamydomonas reinhardtii Exhibit Distinctive Light-induced Structural Changes from Channelrhodopsin-2

Channelrhodopsin-2 (ChR2) from the green alga Chlamydomonas reinhardtii functions as a light-gated cation channel that has been developed as an optogenetic tool to stimulate specific nerve cells in animals and control their behavior by illumination. The molecular mechanism of ChR2 has been extensively studied by a variety of spectroscopic methods, including light-induced difference Fourier transform infrared (FTIR) spectroscopy, which is sensitive to structural changes in the protein upon light activation. An atomic structure of channelrhodopsin was recently determined by x-ray crystallography using a chimera of channelrhodopsin-1 (ChR1) and ChR2. Electrophysiological studies have shown that ChR1/ChR2 chimeras are less desensitized upon continuous illumination than native ChR2, implying that there are some structural differences between ChR2 and chimeras. In this study, we applied light-induced difference FTIR spectroscopy to ChR2 and ChR1/ChR2 chimeras to determine the molecular basis underlying these functional differences. Upon continuous illumination, ChR1/ChR2 chimeras exhibited structural changes distinct from those in ChR2. In particular, the protonation state of a glutamate residue, Glu-129 (Glu-90 in ChR2 numbering), in the ChR chimeras is not changed as dramatically as in ChR2. Moreover, using mutants stabilizing particular photointermediates as well as time-resolved measurements, we identified some differences between the major photointermediates of ChR2 and ChR1/ChR2 chimeras. Taken together, our data indicate that the gating and desensitizing processes in ChR1/ChR2 chimeras are different from those in ChR2 and that these differences should be considered in the rational design of new optogenetic tools based on channelrhodopsins.     

Flash spectroscopic studies of the kinetics of the halorhodopsin photocycle

The photoreactions of halorhodopsin are complicated by the fact that the parent pigment and its photoproducts interact with chloride. Thus, in any photoreaction scheme at least four species have to be accounted for: HR565 and HR578 Cl-, as well as HR640 and HR520 Cl-. A photocycle scheme proposed earlier places the two main photointermediates of halorhodopsin, HR520 Cl- and HR640, into a single photocycle, with a chloride-dependent equilibrium between them [Oesterhelt, D., Hegemann, P., & Tittor, J. (1985) EMBO J. 4, 2351-2356]. This scheme, with the additional feature of direct photoproduction of HR640 from HR565, was tested in this work by using numerical solutions of the appropriate differential equations to simulate flash-induced absorption changes at 500 nm (production of HR520 Cl-) and at 660 nm (production of HR640). The time scale of the simulation was ms following the flash. Comparison of the simulated curves with experimental traces yielded a unique set of three rate constants. The proposed photocycle scheme and these rate constants predict well the shapes and amplitudes of flash traces at various chloride concentrations. It appears from the photocycle scheme, and the numerical values of rate constants, that chloride is bound with high affinity to the parent halorhodopsin molecule, but with much lower affinity to its main photointermediate. This may be the consequence of the fact that in the parent halorhodopsin in the retinal configuration is all-trans, but in the two photointermediates it is 13-cis.