Time-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin
Abstract
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.
Introduction
Channelrhodopsins (ChRs) are light-gated cation channels first identified as sensory photoreceptors in green algae (Nagel et al., 2003). ChRs consist of seven transmembrane helices and a retinal chromophore covalently attached to a conserved lysine residue via a Schiff base linkage. ChRs allow fast ion flux across cell membranes upon illumination and thus are widely utilized as optogenetic tools to precisely control neural activity (Schneider et al., 2015; Ernst et al., 2014; Zhang et al., 2011). As in the well-studied bacterial proton pump bacteriorhodopsin (BR), light absorption triggers all-trans to 13-cis isomerization of the retinal, and induces the transitions of optically distinguishable intermediates that lead to channel opening and closing (photocycle process) (Kuhne et al., 2019; Bamann et al., 2010; Bamann et al., 2008; Figure 1). Spectroscopic studies have provided detailed characterizations of the respective intermediates in the photocycle of ChRs. In the best studied variant, Chlamydomonas reinhardtii ChR2 (CrChR2) (Bamann et al., 2008; Stehfest and Hegemann, 2010; Lórenz-Fonfría et al., 2013), retinal isomerization leads to an intermediate with a red-shiftedabsorption peak, P1520 (K or L in BR), and sequential deprotonation and reprotonation of the Schiff base generate P2390 and P3520 (M and N in BR), respectively (Figure 1). Photocycle transitions are determined in most ChRs by an aspartate–cysteine pair in the retinal binding pocket (DC-pair) (Berndt et al., 2009; Yizhar et al., 2011; Figure 1—figure supplement 1), but the mechanism by which ion pore formation is associated with these residues is still under debate. Recent studies have suggested that the ion conducting pore is formed during the late P2390 intermediate Kuhne et al., 2019, involving the rearrangement of the pore-forming TM helices (Müller et al., 2011), and a two-dimensional crystal study indeed revealed low resolution projection maps of the light-induced helix movement in ChR (Müller et al., 2011; Müller et al., 2015). A detailed mechanism for how retinal isomerization causes conformational changes and how these changes induce ion pore opening has remained elusive.
The recent development of the X-ray free electron laser (XFEL) technology has enabled high-resolution visualization of the time-resolved molecular conformational changes in crystals, thus profoundly advancing our understanding of biological reactions (Tosha et al., 2017; Nango et al., 2016; Suga et al., 2017; Shimada et al., 2017). In this study, we performed the time-resolved serial femtosecond crystallography (TR-SFX) analysis of ChR, which revealed critical conformational changes occurring early in the photocycle that probably induce the cation permeation in the later stages.
Results
Time-resolved spectroscopy of C1C2
We used a chimeric construct between Chlamydomonas reinhardtii ChR1 (CrChR1) and ChR2 (CrChR2) for the TR-SFX experiments (hereafter referred to as C1C2) (Kato et al., 2012; Figure 1—figure supplement 1), which has normal ChR1 photocycle properties (Hontani et al., 2017; Inaguma et al., 2015; Figure 2a,b). C1C2 yielded high-density microcrystals that diffracted XFEL to about 2.3 Å at SACLA (Figure 3a,b). The dark state structure of C1C2 determined at SACLA is almost the same as the previous structure solved by the synchrotron analysis, including the chromophore retinal, the interhelical interactions of the inner and central gates, and the extracellular water access channel in the dark state (Figure 3c,d, and Figure 3—figure supplement 1; Hontani et al., 2017). The C1C2 photocycle was reported in detail recently (Hontani et al., 2017), revealing different accumulations of the P1 and P2 photo-intermediates, as compared to CrChR2. Preceding the TR-SFX experiment, we investigated the photocycle and ion pore formation processes of C1C2. In the flash-photolysis experiments, C1C2 showed similar photo-intermediates as observed in the previous study (Hontani et al., 2017), in detergent micelles or reconstituted lipid membranes (Figure 2a,b and Figure 2—figure supplement 1). The results indicate a rapid rise of 520 nm absorption after photo-excitation, occurring before 10 μs, which corresponds to the formation of the P1 intermediate with an isomerized retinal (Hontani et al., 2017; Nogly et al., 2018), and the subsequent increases at 390 nm and 520 nm show formation of the P2 and P3 intermediates with a deprotonated and reprotonated retinal Schiff base, respectively (Figure 2a). The absorption spectra of P1 and P3 resemble each other and cannot be distinguished in C1C2, but the delayed rise of the 530 nm absorption after 360 μs shows the formation of the P3 intermediate at the later time point. To correlate the different photocycle intermediates to formation of the ion conducting pore, we next measured the C1C2 photocurrents evoked by 7-ns laser flash (Figure 2c and Figure 2—figure supplement 2). The C1C2 photocurrents are initiated by a fast, inward directed charge transfer that likely coincides with the retinal deprotonation (P2 formation), which has not been fully time-resolved in the whole cell patch clamp recordings. Afterwards, passive photocurrents rise with a time constant of 70–140 µs (rather faster at acidic pH) and change direction depending on the membrane voltage. The photocurrent remains after reprotonation of the retinal Schiff base and declines in a biexponential manner, with two voltage dependent time constants of 8–40 ms and 60–130 ms and a kinetic amplitude ratio that depends on the pH, correlating neither precisley with the P2nor the P3 intermediate. The absorption spectrum of ChR depends mostly on the microenvironment of the chromophore retinal, especially on the Schiff base environment (Ritter et al., 2008), and does not need to reflect changes in the ion pore structure . As recently shown for CrChR2 (Kuhne et al., 2019), also C1C2 channel opening takes place only after deprotonation of the retinal Schiff base, during further confirmational changes of the P2 intermediate.
To determine the time points for the TR-SFX study, we next investigated the photocycle of C1C2 in the lipidic cubic phase (LCP) crystals. The flash illumination induced a rapid elevation at 530 nm, confirming retinal isomerization (P1 formation) (Figure 2d and Figure 2—figure supplement 3a, c). The elevation of the 390 nm absorption, which is most prominent in 2.1–21 ms, showed the P2 formation, but the subsequent increase at the 530 nm absorption was not observed in the crystal (Figure 2d and Figure 2—figure supplement 1 and Figure 2—figure supplement 3a,c), suggesting that the transition from P2 to P3 is hindered in the crystal, although the C1C2 crystal has only minimal interactions at the extracellular membrane region (Figure 3—figure supplement 2). Nonetheless, the accumulation of the P2 intermediate is more prominent in the crystals, as the spectral increase at 390 nm is relatively higher, as compared to the similar increase detected in micelles and reconstituted membranes (Figure 2a,d and Figure 2—figure supplement 1, Figure 2—figure supplement 3b,d). In addition, the flash photo-activation experiment using HEK cells showed that the photocurrent rose in between 100–1000 μs, regardless of pH, which corresponds to a second late P2 intermediate (Figure 2c and Figure 2—figure supplement 2), in a good agreement with the previous studies in ChR2 (Kuhne et al., 2015; Kuhne et al., 2019). Therefore, we selected the time points corresponding to the P1 to P2 transition in the crystals for the TR-SFX study, to understand the detailed mechanism of the light-dependent channel opening in ChR.
Structural change in ChR
The LCP matrix containing C1C2 microcrystals was introduced to the XFEL chamber in a continuous flow. Since the crystallized sample was swollen and almost in the liquid sponge phase, we added paraffin and solid LCP matrix to enable stable and constant sample flow to the FEL beam (detailed in Materials and methods). The diffraction images were collected at 30 Hz, with the 470 nm pump-pulse laser at 15 Hz and the delayed time points at 1 μs, 50 μs, 250 μs, 1 ms, and 4 ms (Figure 3—figure supplement 3a,c and Supplementry file 1). The absence of light contamination was checked by the second dark data set, obtained with a 10 Hz pump laser pulse (Figure 3—figure supplement 3c,d). The difference-Fourier maps were calculated at 2.5 Å for each delayed time point, with the dark state SFX structure as the reference (Figure 4), and they clearly revealed the structural changes in the protein core region of C1C2 (Figure 4—figure supplement 1; Figure 4—figure supplement 2 and Video 1). Due to the low quantum efficiency of the retinal isomerization, which is estimated as only about 30% in C1C2 (Hontani et al., 2017), we calculated the extrapolated structure factor amplitudes and the models were refined against the extrapolated maps (detailed in Materials and methods and Figure 4—figure supplement 3). We also calculated the Fclight–Fcdark maps using the models, which shows a similar appearance with the difference-Fourier map and thus confirms the conformational change (Figure 4—figure supplement 4).
Since the retinal isomerization reaction is generally completed within a few picoseconds (Hontani et al., 2017; Neumann-Verhoefen et al., 2013; Bühl et al., 2018), the observed changes correspond to the events following the retinal isomerization. The isomerization of the retinal induces conformational changes of the retinal-linked residue, Lys296, and TM7 (Figure 4, left panels). The positive densities along Lys296 suggest its downward shift (toward the extracellular side), which is detected from 1 μs and becomes most prominent at 4 ms (Figure 4, left panels). In accordance with this shift, positive and negative difference densities are observed around the proximal residues in TM7, such as Trp299, indicating an overall downward shift of the middle portion of TM7. The similar downward shift of TM7 (Helix G) was the most distinct feature in the L-like intermediate of BR, observed in both the cold-trapped synchrotron structure (Neutze et al., 2002) and the recent TR-SFX study (Nango et al., 2016). In the TR-SFX experiment of BR, the primary retinal isomerization reaction induced local distortion around the C13–C14 double bond of the retinal, which was subsequently propagated as conformational changes to other protein regions, also involving a downward shift of TM7. Therefore, the retinal-attached residue in C1C2 probably undergoes similar conformational changes, although the primary retinal isomerization was not visualized, due to the limited resolution of the current TR-SFX experiment. Notably, at 50 μs we observed positive and negative difference densities at the putative internal proton acceptor residue, Asp292 (Figure 4, left panels), indicating a slight downward shift of this residue, which might be associated with the proton accepting reaction.
Retinal kink-induced outward shift of TM3
The most prominent change occurs in TM3 (Figure 4, right panels). A strong negative density appears at Cys167, accompanied by a complementary positive density at the distal site from the retinal, suggesting an outward shift of Cys167. In addition, a strong positive density appears near Pro168 on the outside of the TM3 segment, which indicates an outward shift of the cytoplasmic segment of TM3, originating at the helix kink at the Pro168 residue. Along with these changes, negative and positive densities appear at the C13 methyl group and on the lateral side of the polyene chain of the retinal, respectively. This kink and lateral shift of retinal allows lateral movement of TM6 as well as TM3 (Figure 4, right panels). The current difference Fourier map is in excellent agreement with the recent MD simulation combined with QM/MM calculations (Cheng et al., 2018), in which the retinal adopts a kinked conformation toward TM3 upon the Schiff base deprotonation, probably due to the destabilized resonance structure of the polyene chain (Figure 5), although the QM/MM simulation focuses only on the molecular events after the Schiff base deprotonation and thus does not exactly match the experimental time scale. Together, the sequential changes in the difference-Fourier maps suggest that the retinal kink occurs, followed by the consequent outward shift of TM3. Consistently, the positive/negative peaks at the retinal appear slightly earlier, whereas the change in TM3 becomes slower and reaches the maximum at the later time points of 1 ms and 4 ms (Figure 5).
The current TR-SFX experiment revealed two critical conformational changes occurring early in the photocycle process: the outward shift of TM3 and the downward shift of TM7, both of which constitute gate constriction sites along the ion permeation pathway in the dark state (Figure 3). TM7 contributes essentially to the ‘central gate’, in which Asn297 forms hydrogen bond interactions with Ser102 (TM1) and Glu129 (TM2). Therefore, the TM7 shift, as well as the outward shift of Asp292 upon proton acceptance, should lead to rearrangements of the central gate residues. The intracellular end of TM3 shapes the ‘inner-gate’, in which His173 forms ionic interactions with Arg307 in TM7 and Glu121 and Glu122 in TM2. The recent MD simulation showed a similar retinal kink and subsequent water influx from the intracellular side (Figure 5; Cheng et al., 2018). Therefore, the conformational change in the middle of TM3 might propagate to the intracellular segment and induce the opening of the inner gate.
Cys167 and Asp195 have been identified as important residues for channel gating and are thus termed DC-gate or DC-pair, which is conserved among most ChRs. FTIR and structural studies have suggested that Asp195 functions as the internal proton donor for the deprotonated Schiff base, via an internal hydrated water molecule and Cys167 (Nack et al., 2010; Bamann et al., 2010; Lórenz-Fonfría et al., 2013). A previous FTIR study showed an increase of the S-H stretching frequency of the Cys167 side chain upon illumination (Ito et al., 2014). The DC-pair residues are critical for the fast channel kinetics of ChRs, as mutations of any of them drastically slowed the open/close kinetics of the channel, resulting in step function opsin (SFO) mutants, with long-lasting photocurrents that can be switched on and off by different wavelengths of light (Berndt et al., 2009; Yizhar et al., 2011). Previous 2D crystallographic studies of the CrChR2 C128T mutant (corresponding to Cys167 in C1C2 [Figure 1—figure supplement 1]) suggested rearrangements of TM2, 6, and 7 in the open channel conformation, that was here induced by the stationary light (Müller et al., 2011; Müller et al., 2015). These conformational changes are partly consistent with the current TR-SFX analysis. However, the TM3 shift is only observed in the current study, indicating the critical role of the Cys residue in the TM3 shift. The movement of TM7, on the other hand, is similarly observed in both experiments, although to different degrees. We therefore propose that the conformational changes in TM7 and TM3 are the two important triggers that cooperatively induce channel opening by eliciting sequential conformational changes in the surrounding helices, such as TM2 and TM6 (Figure 6), as observed in the ‘open’ conformation of the SFO mutant of CrChRs (Müller et al., 2015). This model is also in good agreement with the results of the vibrational spectroscopic study showing a two-phasic helix hydration occurring during the early and late stages of the P2 intermediate, in which only the latter leads to the fully-open channel pore (Lórenz-Fonfría et al., 2015). Furthermore, site-directed labeling with infrared sensitive, unnatural amino acids subsequently showed more specificly, that espacially hydration changes in the intracellular gate are correlated with channel opening (Krause et al., 2019). In summary, our TR-SFX experiment has revealed two critical conformational changes in TM7 and TM3 that occur during the early part of the photocycle, and these triggers would induce water and proton influx and subsequent pore dilation, although only slight structural perturbations in the gate constrictions were observed in the current study, as they were probably affected by the crystal packing interactions (Figure 4—figure supplement 5).
Discussion
Channel pore gating in ChR
The current results highlight the diversity and similarity of rhodopsin proteins. TM7 is covalently attached to the retinal and thus undergoes a similar shift in other rhodopsin proteins, such as the proton pump BR (Nango et al., 2016), and the Na+ pump KR2 (Skopintsev et al., 2020), albeit with some variations. Therefore, this structural movement is probably conserved in rhodopsin proteins. Most notably, the retinal kink and TM3 shift were also observed in KR2, even though it is evolutionally and functionally distant from ChRs, while such a conformational change was not observed in BR. In BR, the retinal is transiently kinked toward the opposite direction, but only during the very early time after photoexcitation (around 10 ps to 16 ns), and causes totally different changes in the protein part. Therefore, the retinal kink and subsequent conformational changes are likely to be affected by the residues constituting the retinal pocket. Indeed, the outward shift of TM3 was not observed in the Cys167 mutant of C1C2, as elucidated by the light-activated projection maps on the 2D crystallographic study (Müller et al., 2015). A recent MD simulation demonstrated the hydrogen bonding interactions between the deprotonated Schiff base and the hydroxyl group of Thr166, located beside the Cys167 residue (Cheng et al., 2018), which are likely to affect the retinal conformation. The same threonine is highly conserved among microbial rhodopsins, and in CrChR2 it was also proposed to couple the protonation changes in the counterion region to the hydrogen bonding changes in the DC-pair and the conformational changes of TM3 (Watanabe et al., 2012). Together, the hydrogen-bonding network beside the retinal, involving the DC-pair, may affect the retinal conformation after isomerization, in good agreement with the experimental evidence that the residues lining the retinal pocket affect channel kinetics (Oda et al., 2018). These results would pave the way for the development of novel optogenetic tools with different kinetics.
Materials and methods
Laser flash photolysis experiment
Request a detailed protocolThe time-evolution of the transient absorption changes of photo-excited C1C2 was measured as previously described (Inoue et al., 2013; Inoue et al., 2016). The purified protein sample was solubilized in buffer containing 150 mM NaCl, 50 mM Tris-HCl (pH 8.0), 0.05% n-dodecyl-β-D-maltoside (DDM), 5% glycerol, and 0.01% cholesteryl hemisuccinate (CHS). For the measurement of C1C2 in liposomes, the purified C1C2 protein was mixed with a lipid mixture of POPE and POPG (molar ratio = 3: 1) at a protein-to-lipid molar ratio of 1: 50 (the final concentrations of C1C2 and lipids were adjusted to 5 μM and 250 μM, respectively), and the DDM was removed by BioBeads (SM-2, Bio-Rad). The sample solution was placed in a quartz cuvette with a 1 cm optical path length and illuminated with a 480 nm beam generated by an OPO system, pumped by the third harmonics generation of a nanosecond Nd3+: YAG laser (LT-2214, LOTIS TII, Minsk, Republic of Belarus). The intensities of the transmitted probe light from a Xe arc lamp (L8004, Hamamatsu Photonics, Japan) were measured before and after laser excitation with an ICCD linear array detector (C8808-01, Hamamatsu, Japan) at 24°C, and transient absorption spectra were obtained by calculating the ratios between them. Ninety spectra were averaged to improve the signal-to-noise (S/N) ratio. To obtain the detailed time evolution of the transient absorption change after photoexcitation, the change in the intensity of the monochromated output of a Xe arc lamp (L9289-01, Hamamatsu Photonics, Japan), passed through the sample at 24°C, was observed by a photomultiplier tube (R10699, Hamamatsu Photonics, Japan) at 375, 445, and 520 nm. The signal from the photomultiplier tube was averaged to improve the S/N ratio and stored in a digital-storage-oscilloscope (DPO7104, Tektronix, Japan). The obtained kinetics data were fitted with a multi-exponential function.
Dark adaptation kinetics of C1C2 was measured for the 0.05% DDM-solubilized sample at 24°C, as previously described (Inoue et al., 2016). The dark-adapted sample was illuminated for 1 min by using the output from a 1-kW tungsten–halogen projector lamp (Master HILUX-HR, Rikagaku, Japan), through an interference filter at 460 nm. After turning off the light, the spectra or the absorption at a specific wavelength were measured by a UV-visible spectrometer (V-730, JASCO).
Flash photo-activation current measurement
Request a detailed protocolHEK293 cells were cultured in Dulbecco’s Modified Medium (Biochrom), supplemented with 10% (v/v) fetal bovine serum (FBS Superior; Biochrom), 1 μM all-trans retinal, and 100 µg/ml penicillin/streptomycin (Biochrom, Berlin, Germany). Cells were seeded on poly-lysine coated glass coverslips at a concentration of 1 × 105 cells/ml, and transiently transfected using the FuGENE HD Transfection Reagent (Promega) 28–48 hr before measurement. C1C2 was fused to an mCherry-fluorophor and expressed under the control of the CMV-promotor.
Patch pipettes were prepared from borosilicate glass capillaries (G150F-3; Warner Instruments) using a P-1000 micropipette puller (Sutter Instruments), and fire polished with a final pipette resistance of 1.5–3.0 MΩ. Fluorescent cells were identified using an Axiovert 100 inverted microscope (Carl Zeiss). Single ns laser pulses were generated by an Opolette Nd:YAG laser/OPO system (OPOTEK), selected by an LS6ZM2 shutter system (Vincent Associates), and delivered to the cells through the optics of the microscope and a 90/10 beamsplitter (Chroma). The laser intensity was set to 5% using the built-in motorized variable attenuator, resulting in a pulse energy of 100 ± 20 µJ/mm2. Photocurrents were recorded with an AxoPatch 200B amplifier (Molecular Devices) at a maximal bandwidth of 100 kHz and digitized using a DigiData 1440A digitizer (Molecular Devices) at a sampling rate of 250 kHz. Membrane resistance was typically >1 GΩ, while access resistance was below 10 MΩ. Pipette capacity, series resistance, and cell capacity compensations were applied. A 140 mM NaCl agar bridge was employed for the reference bath electrode. Intracellular pipette and extracellular bath solutions contained 110 mM NaCl, 1 mM KCl, 1 mM CsCl, 2 mM CaCl2, and 2 mM MgCl2, and were buffered to pH 6.0, 7.2, or 8.0 by 10 mM MES, HEPES, or TRIS (with glucose added to 290 mOsm for intracellular and 310 mOsm for extracellular solutions). Measurements were controlled by the pCLAMP software (Molecular Devices). Post measurement photocurrent recordings were first baseline corrected and time shifted to align with the rising edges of the Q-switch signals of the laser pulse, using the Clampfit 10.4 software (Molecular Devices). They were then binned to 50 logarithmically spaced data points per temporal decade, normalized to peak photocurrents at −60 mV, and averaged for up to 10 individual repetitions for each cell and voltage condition with a custom-written Matlab script (MathWorks). Kinetic analysis was performed in Origin 9.1 (OriginLab). Passive photocurrents were calculated by subtracting the photocurrents at 0 mV, which were recorded in the absence of an electrochemical gradient and represent the inner protein charge transfer reactions. For sufficient statistical significance, each measurement was repeated multiple times on different biological replicates, in at least three independent experiments. The exact number of biological replicates for each measurement is provided in the figure legend.
Time-resolved visible absorption spectroscopy
Request a detailed protocolThe photocycle reaction was induced by a 6-ns, 470 nm pump pulse from an optical parametric oscillator (OPO) (NT230, EKSPLA), and the spectral changes during the reaction were measured by using a fiber-coupled spectrometer (Flame-S, Ocean Optics), with a microsecond white-light probe pulse from a Xe-flash lamp (L11316-11-11, Hamamatsu Photonics). A portion of the probe light (5%) was used to correct the pulse-to-pulse fluctuation of the probe light intensity. The pump pulse and the probe light were focused on the sample point, with beam diameters of 300 μm and 40 μm, respectively. The pump energy was reduced to 15 μJ (0.21 nJ/μm [Schneider et al., 2015]), to avoid sample damage by the repetitive pump illuminations during the data accumulation. The pump pulse, the probe light source, and the spectrometer were synchronized by pulse generators (DG645, Stanford Research Systems), and the pump-probe delay time was adjusted with a timing jitter of ±20 ns. The repetition rates of the pump and the probe were 0.5 Hz and 1 Hz, respectively. In the solution measurement, the sample was packed in a spinning quartz cell with an optical path length of 1 mm and spun at 50 rpm during the measurement. The sample concentration was ~10 mg/ml, dissolved in 50 mM Tris-HCl, pH 8.0, containing 150 mM NaCl and 5% glycerol. In the microcrystal measurement, we used the C1C2 crystals obtained under similar conditions as in the TR-SFX (below). Specifically, a pillar of 30 μl protein-laden LCP, approximately 2 cm in length, was extruded from the syringe and soaked in 900 μl of the crystallization buffer, containing 100 mM HEPES, pH 7.0, 280 mM KH2PO4, and 30% PEG500DM, and incubated for 1 week in the dark at 20°C, which yielded crystals in a liquid sponge phase floating on the buffer solution. The sponge phase layers were transferred by pipetting and packed in two quartz windows with a 50 μm-thick spacer. Only a few absorption signals from the non-crystallized area were observed, and we used the transmitted light intensity from the non-crystallized area as the reference to calculate the absorption spectra of C1C2. The sample temperature was kept at 20°C for both the solution and microcrystal measurements. To decrease the noise level of the time-resolved spectra, singular value decomposition (SVD) (Shrager and Hendler, 1998) was used with Igor Pro (Wave Matrix). In the SVD, matrix A comprising the intensity changes at each wavelength and time (column and row, respectively) was decomposed according to A = USVT, where S represents the diagonal matrix with the singular values in decreasing order. This allows an evaluation of the number of significant orthonormal basis vectors for the wavelength (columns of U) and the time (columns of V). By reconstructing the matrix A’ from only the most significant components, random noise was eliminated from the raw spectra. The time-resolved spectra for the microcrystal and solution samples were well reconstructed by two and three components, respectively. Based on the SVD, global fitting was then performed for the time-resolved reconstructed spectra with Igor Pro (Wave Matrix). In the global fitting, the spectra were fitted with a sum of exponentials with apparent time constants τi as global parameters, where i is the number of components used for the spectral reconstruction in the SVD analysis.
Crystallization of C1C2
Request a detailed protocolThe vector construction and protein production and purification were performed according to the previous report (Kato et al., 2012). C1C2 was finally purified by size-exclusion chromatography, in 150 mM NaCl, 50 mM Tris-HCl, pH 8.0, 5% glycerol, 0.05% DDM, and 0.01% CHS, and concentrated to ~15 mg/ml for crystallization. C1C2 was mixed with monoolein (Sigma) in a 2:3 protein to lipid ratio (w/w), using a 100 μl volume Gastight syringe (Hamilton). A pillar of 30 μl protein-laden LCP, approximately 2 cm in length, was extruded from the syringe and soaked in 900 μl of the crystallization buffer, containing 100 mM MES, pH 6.9, 100 mM Na formate, and 30% PEG500DM, and incubated for 2–3 weeks in the dark at 20°C. Crystals were observed with a Hirox KH8700 digital microscope.
Sample preparation for LCP-SFX
Request a detailed protocol100 μl of LCP pillars with grown C1C2 microcrystals were collected and loaded into a Hamilton syringe and mixed with an equivalent volume of buffer-containing LCP (2:3 ratio of the buffer containing 150 mM NaCl and 50 mM Tris-HCl, pH 8.0, to monoolein). A 10 μl aliquot of liquid paraffin (Wako) was further added to the LCP material and homogenously mixed to enable smooth flow in the SFX experiment. All of these procedures were performed in the dark under red LED lights (over 600 nm wavelength).
Diffraction experiment at SACLA-FEL
Request a detailed protocolSingle-shot XFEL data collection was performed, using femtosecond X-ray pulses from the SACLA at BL3. The pulse parameters of SACLA were as follows: pulse duration, <10 fs; X-ray energy, 7 keV; energy bandwidth, 33.5013 eV (FWHM); pulse flux, 2.0 × 1011 photon/pulse photons per pulse; beam size 1.5 μm (H) × 1.5 μm (W); repetition rate, 30 Hz. The C1C2 crystals were mixed with protein-less LCP and paraffin oil, loaded into an injector with a 75 μm inner diameter nozzle, and set in a diffraction chamber filled with helium gas, in a setup called Diverse Application Platform for Hard X-ray Diffraction in SACLA (DAPHNIS). The flow rate was set to 2.5 μl/min (9.4 mm/s). The diffraction patterns were recorded on a multiport CCD detector with an eight sensor module (Kameshima et al., 2014). The excitation laser pulses were provided at 15 Hz, and the XFEL pulses had a repetition rate of 30 Hz. Each ‘pump-on’ image was followed by ‘pump-off’ images, which were recorded separately. The pump-on images were used to analyze the structural differences, whereas the diffraction data for the dark-state were collected in a separate run to avoid light contamination.
Pump laser setup
Request a detailed protocolTo advance the C1C2 samples to the excited state, a 6 ns visible light pump pulse with a wavelength of 470 nm was provided to the samples at Δt = 1 μs, 50 μs, 250 μs, 1 ms, and 4 ms, from an OPO (NT230, EKSPLA), which was Q-switched at 15 Hz using a pulse generator (DG645, Stanford Research Systems). The pulse generator was also used to synchronize the pump pulse with the XFEL pulse and control the delay time (Δt) with a timing jitter of <1 ns. Samples were thus irradiated from both sides of the LCP microjet, a geometry that was chosen to optimize the crystal excitation efficiency. The pump focal size was set to 40 μm (FWHM) and the pump energy was 20 μJ (10 μJ from each direction) at the sample point. The pump beam center was aligned 6 μm upstream from the XFEL spot. The pump light was scattered, as the LCP microjet is semi-transparent. To avoid the interference from the pump scattering, the sample flow rate (9.4 mm/s) was set to be the same as that in the TR-SFX analysis of BR, in which each pump-illuminated crystal for the pump-on XFEL images was separated by 630 µm. A portion of the pump beam (5%) was isolated with a beam sampler, and its photodiode-detected signal was used to monitor the laser status and to tag the diffraction image with ‘light-on’. Since the repetition rate of the pump laser was half of that of the XFEL, the ‘light-on’ data and the ‘dark’ data were collected alternately.
Data processing and structural analysis
Request a detailed protocolData collection was monitored by a real-time data processing pipeline (Nakane et al., 2016), based on Cheetah (Barty et al., 2014). Diffraction images with more than 20 spots were considered as hits and processed in CrystFEL, version 0.6.3 (White et al., 2012), with the parameters of hitfinderAlgorithm = 6, hitfinderADC = 40, hitfinderMinSNR = 5. Diffraction spots for indexing were detected by the built-in spot finding algorithms ‘zaef’ (Zaefferer, 2000) or ‘peakfinder8’ (Barty et al., 2014) in CrystFEL, and indexing was performed with DirAx (Duisenberg, 1992). Intensities were merged by Monte Carlo integration with the process_hkl command in the CrystFEL suite, with the per-image resolution cutoff (--push-res 1.2). Data collection statistics are summarized in Supplementary file 1. As the cell parameters do not vary among the different delayed-time points, the difference Fourier electron density maps (|Fobs|light–|Fobs|dark)·exp[iΦcalc] were calculated by FFT (Ten Eyck, 1973) in the CCP4 suite (Winn et al., 2011), with phases calculated from the refined dark state model.
Structural refinement
Request a detailed protocolThe structure was solved by the molecular replacement method implemented in Phenix Phaser-MR, using the model of C1C2 (PDB: 3UG9). The structure was manually modified to fit into the electron density maps, using the program Coot (Emsley et al., 2010), and then refined with the programs Phenix (Adams et al., 2010) and Refmac5 in the CCP4 suite (Vagin et al., 2004). The model structures of the respective delayed time points were generated by real-space refinement against the difference-Fourier maps by COOT, with the restraints of protein geometry and retinal planarity between C5and C13. Extrapolated data were calculated by a linear approximation, as follows: FExtra = (Folight − Fodark)/R + Fodark, where R represents the activation ratio of molecules and FExtra represents the extrapolated structure factor amplitude. Folight was scaled to Fodark by using the scaleit program in the CCP4 suite (Howell and Smith, 1992) before calculating FExtra. The activation ratio R (the ratio of molecules in the excited state after light illumination) was determined by calculating extrapolated maps with phases of the dark state at different activation ratios in steps of 0.01, until a feature of the dark state at Cys167 disappeared. Next, 4 ms model was manually adjusted to fit the extrapolated map and refined by Refmac5 against FExtra and its estimated error. Then, each time point model was refined by Refmac5 with external structure restraints using 4 ms structure by the Geman–McClure function (Murshudov et al., 2011) the refinement statistics are. Figures were prepared with Cuemol (http://www.cuemol.org).
Data availability
Coordinates and structure factors for the TR-SFX structures have been deposited in the PDB with the accession codes 7C86 (dark), 7E6Y (1 μs), 7E6Z (50 μs), 7E70 (250 μs), 7E71 (1 ms), and 7E6X(4 ms). Raw diffraction images have been deposited in the Coherent X-ray Imaging Data Bank (accession ID 150, https://doi.org/10.11577/1768368).
-
RCSB Protein Data BankID 7C86. Time-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin: Dark state structure.
-
Coherent X-ray Imaging Data BankTime-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin.https://doi.org/10.11577/1768368
-
Worldwide Protein Data BankTime-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin: 1 µs structure.https://doi.org/10.2210/pdb7E6Y/pdb
-
Worldwide Protein Data BankTime-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin: 50 µs structure.https://doi.org/10.2210/pdb7E6Z/pdb
-
Worldwide Protein Data BankTime-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin: 250 µs structure.https://doi.org/10.2210/pdb7E70/pdb
-
Worldwide Protein Data BankTime-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin: 1 ms structure.https://doi.org/10.2210/pdb7E71/pdb
-
Worldwide Protein Data BankTime-resolved serial femtosecond crystallography reveals early structural changes in channelrhodopsin: 4 ms structure.https://doi.org/10.2210/pdb7E6X/pdb
References
-
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica Section D Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
Spectral characteristics of the photocycle of channelrhodopsin-2 and its implication for channel functionJournal of Molecular Biology 375:686–694.https://doi.org/10.1016/j.jmb.2007.10.072
-
Cheetah: software for high-throughput reduction and analysis of serial femtosecond X-ray diffraction dataJournal of Applied Crystallography 47:1118–1131.https://doi.org/10.1107/S1600576714007626
-
Ultrafast Protein Response in Channelrhodopsin-2 Studied by Time-Resolved Infrared SpectroscopyThe Journal of Physical Chemistry Letters 9:7180–7184.https://doi.org/10.1021/acs.jpclett.8b03382
-
Indexing in single-crystal diffractometry with an obstinate list of reflectionsJournal of Applied Crystallography 25:92–96.https://doi.org/10.1107/S0021889891010634
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Microbial and animal rhodopsins: structures, functions, and molecular mechanismsChemical Reviews 114:126–163.https://doi.org/10.1021/cr4003769
-
Reaction dynamics of the chimeric channelrhodopsin C1C2Scientific Reports 7:7217.https://doi.org/10.1038/s41598-017-07363-w
-
Identification of heavy-atom derivatives by normal probability methodsJournal of Applied Crystallography 25:81–86.https://doi.org/10.1107/S0021889891010385
-
Chimeras of Channelrhodopsin-1 and -2 from Chlamydomonas reinhardtii exhibit distinctive Light-induced structural changes from Channelrhodopsin-2Journal of Biological Chemistry 290:11623–11634.https://doi.org/10.1074/jbc.M115.642256
-
A light-driven sodium ion pump in marine BacteriaNature Communications 4:1678.https://doi.org/10.1038/ncomms2689
-
A natural light-driven inward proton pumpNature Communications 7:13415.https://doi.org/10.1038/ncomms13415
-
Water-containing hydrogen-bonding network in the active center of channelrhodopsinJournal of the American Chemical Society 136:3475–3482.https://doi.org/10.1021/ja410836g
-
Development of an X-ray pixel detector with multi-port charge-coupled device for X-ray free-electron laser experimentsReview of Scientific Instruments 85:033110.https://doi.org/10.1063/1.4867668
-
Early formation of the ion-conducting pore in channelrhodopsin-2Angewandte Chemie International Edition 54:4953–4957.https://doi.org/10.1002/anie.201410180
-
Projection structure of channelrhodopsin-2 at 6 å resolution by electron crystallographyJournal of Molecular Biology 414:86–95.https://doi.org/10.1016/j.jmb.2011.09.049
-
Light-induced Helix movements in channelrhodopsin-2Journal of Molecular Biology 427:341–349.https://doi.org/10.1016/j.jmb.2014.11.004
-
REFMAC5 for the refinement of macromolecular crystal structuresActa Crystallographica Section D Biological Crystallography 67:355–367.https://doi.org/10.1107/S0907444911001314
-
The DC gate in Channelrhodopsin-2: crucial hydrogen bonding interaction between C128 and D156Photochemical & Photobiological Sciences 9:194–198.https://doi.org/10.1039/b9pp00157c
-
Data processing pipeline for serial femtosecond crystallography at SACLAJournal of Applied Crystallography 49:1035–1041.https://doi.org/10.1107/S1600576716005720
-
Ultrafast infrared spectroscopy on channelrhodopsin-2 reveals efficient energy transfer from the retinal chromophore to the proteinJournal of the American Chemical Society 135:6968–6976.https://doi.org/10.1021/ja400554y
-
Bacteriorhodopsin: a high-resolution structural view of vectorial proton transportBiochimica Et Biophysica Acta (BBA) - Biomembranes 1565:144–167.https://doi.org/10.1016/S0005-2736(02)00566-7
-
Crystal structure of the red light-activated channelrhodopsin chrimsonNature Communications 9:3949.https://doi.org/10.1038/s41467-018-06421-9
-
Monitoring Light-induced structural changes of Channelrhodopsin-2 by UV-visible and fourier transform infrared spectroscopyJournal of Biological Chemistry 283:35033–35041.https://doi.org/10.1074/jbc.M806353200
-
Biophysics of channelrhodopsinAnnual Review of Biophysics 44:167–186.https://doi.org/10.1146/annurev-biophys-060414-034014
-
Some pitfalls in curve-fitting and how to avoid them: a case in pointJournal of Biochemical and Biophysical Methods 36:157–173.https://doi.org/10.1016/S0165-022X(98)00007-4
-
Evolution of the channelrhodopsin photocycle modelChemPhysChem 11:1120–1126.https://doi.org/10.1002/cphc.200900980
-
Crystallographic fast fourier transformsActa Crystallographica Section A 29:183–191.https://doi.org/10.1107/S0567739473000458
-
REFMAC5 dictionary: organization of prior chemical knowledge and guidelines for its useActa Crystallographica Section D Biological Crystallography 60:2184–2195.https://doi.org/10.1107/S0907444904023510
-
Structural model of channelrhodopsinJournal of Biological Chemistry 287:7456–7466.https://doi.org/10.1074/jbc.M111.320309
-
CrystFEL : a software suite for snapshot serial crystallographyJournal of Applied Crystallography 45:335–341.https://doi.org/10.1107/S0021889812002312
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
New developments of computer-aided crystallographic analysis in transmission electron microscopyJournal of Applied Crystallography 33:10–25.https://doi.org/10.1107/S0021889899010894
Article and author information
Author details
Funding
Japan Society for the Promotion of Science (16H06294)
- Osamu Nureki
Japan Society for the Promotion of Science (18J21256)
- Kazumasa Oda
Japan Society for the Promotion of Science (17H05000)
- Tomohiro Nishizawa
Japan Science and Technology Agency (JPMJPR14L9)
- Kunio Hirata
Japan Science and Technology Agency (JPMJPR14L8)
- Tomohiro Nishizawa
Japan Agency for Medical Research and Development (Platform Project for Supporting Drug Discovery andLife Science Research (Basis for Supporting Innovative Drug Discovery and Life Science Research(BINDS)))
- Kunio Hirata
Deutsche Forschungsgemeinschaft (SFB1078 (B2))
- Johannes Vierock
Cluster of Excellence in Inflammation Research (UniCat)
- Johannes Vierock
Cluster of Excellence in Inflammation Research (BIG-NSE)
- Johannes Vierock
Cluster of Excellence in Inflammation Research (E4)
- Peter Hegemann
Cluster of Excellence in Inflammation Research (SPP 1926)
- Johannes Vierock
European Research Council (ERC-2016-StG 714762)
- Johannes Vierock
Hertie Foundation (Senior Research Professor)
- Peter Hegemann
Ministry of Education, Culture, Sports, Science and Technology (X-ray Free-Electron Laser Priority Strategy Program)
- So Iwata
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Sae Okazaki and Masae Miyazaki for assistance in obtaining the C1C2 crystals, and Rie Tanaka, Toshi Arima, Ayumi Yamashita, and Jun Kobayashi for assistance in the collection of the X-ray data diffraction data at SACLA. We also thank the beamline staff scientists at SACLA and BL32XU of SPring-8 (Hyogo, Japan). The XFEL experiments were performed at the BL2 of SACLA, with the approval of the Japan Synchrotron Radiation Research Institute (JASRI) (proposal nos. 2015B8053, 2016A8054, 2016B8067, 2017A8022, and 2017B8035), and SPring-8 BL32XU (20170511, 20170721, and 20170728), with the approval of RIKEN. We acknowledge computational support from the SACLA HPC system and Mini-K super computer system.
Copyright
© 2021, Oda et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,343
- views
-
- 553
- downloads
-
- 48
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.
-
- Structural Biology and Molecular Biophysics
Osmotic stress and chloride regulate the autophosphorylation and activity of the WNK1 and WNK3 kinase domains. The kinase domain of unphosphorylated WNK1 (uWNK1) is an asymmetric dimer possessing water molecules conserved in multiple uWNK1 crystal structures. Conserved waters are present in two networks, referred to here as conserved water networks 1 and 2 (CWN1 and CWN2). Here, we show that PEG400 applied to crystals of dimeric uWNK1 induces de-dimerization. Both the WNK1 the water networks and the chloride-binding site are disrupted by PEG400. CWN1 is surrounded by a cluster of pan-WNK-conserved charged residues. Here, we mutagenized these charges in WNK3, a highly active WNK isoform kinase domain, and WNK1, the isoform best studied crystallographically. Mutation of E314 in the Activation Loop of WNK3 (WNK3/E314Q and WNK3/E314A, and the homologous WNK1/E388A) enhanced the rate of autophosphorylation, and reduced chloride sensitivity. Other WNK3 mutants reduced the rate of autophosphorylation activity coupled with greater chloride sensitivity than wild-type. The water and chloride regulation thus appear linked. The lower activity of some mutants may reflect effects on catalysis. Crystallography showed that activating mutants introduced conformational changes in similar parts of the structure to those induced by PEG400. WNK activating mutations and crystallography support a role for CWN1 in WNK inhibition consistent with water functioning as an allosteric ligand.