Proton transfer pathway in anion channelrhodopsin-1
Abstract
Anion channelrhodopsin from Guillardia theta (GtACR1) has Asp234 (3.2 Å) and Glu68 (5.3 Å) near the protonated Schiff base. Here, we investigate mutant GtACR1s (e.g., E68Q/D234N) expressed in HEK293 cells. The influence of the acidic residues on the absorption wavelengths was also analyzed using a quantum mechanical/molecular mechanical approach. The calculated protonation pattern indicates that Asp234 is deprotonated and Glu68 is protonated in the original crystal structures. The D234E mutation and the E68Q/D234N mutation shorten and lengthen the measured and calculated absorption wavelengths, respectively, which suggests that Asp234 is deprotonated in the wild-type GtACR1. Molecular dynamics simulations show that upon mutation of deprotonated Asp234 to asparagine, deprotonated Glu68 reorients toward the Schiff base and the calculated absorption wavelength remains unchanged. The formation of the proton transfer pathway via Asp234 toward Glu68 and the disconnection of the anion conducting channel are likely a basis of the gating mechanism.
Introduction
Anion channelrhodopsins (ACRs) are light-gated anion channels that undergo photoisomerization at the retinal chromophore, which is covalently attached to a conserved lysine residue via the protonated Schiff base, from all-trans to 13-cis. Natural ACRs were identified in the cryptophyte Guillardia theta (GtACR1 and GtACR2) (Govorunova et al., 2015). ACRs hyperpolarize the membrane through anion import and can widely be used as neural silencing tools in optogenetics (Wiegert et al., 2017; Miyazaki et al., 2019). Microbial rhodopsins have acidic residues or Cl– at the Schiff base moiety to stabilize the protonated Schiff base as counterions. Counterions play a major role in determining the absorption wavelength and the function of the protein (Tsujimura and Ishikita, 2020). The X-ray crystal structures of GtACR1 show that two acidic residues, Glu68 and Asp234, exist at the corresponding positions (Figure 1; Kim et al., 2018; Li et al., 2019).
It was proposed that both Glu68 and Asp234 were protonated in GtACR1 (Kim et al., 2018; Sineshchekov et al., 2016; Yi et al., 2016; Kandori, 2020) in contrast to other microbial rhodopsins because the absorption wavelengths remain unchanged upon the E68Q and D234N mutations (Kim et al., 2018; Sineshchekov et al., 2016; Yi et al., 2016). Indeed, the C=C stretching frequency of the retinal is not significantly affected upon the E68Q and D234N mutations in resonance Raman spectroscopy, which implies that the electrostatic interaction between the retinal and protein environment remains unchanged (Yi et al., 2016). In addition, the C=O stretching frequency for a protonated carboxylate, which is observed in the wild-type GtACR1, disappears in the E68Q (Yi et al., 2017; Dreier et al., 2021) and D234N (Kim et al., 2018) GtACR1s according to Fourier transform infrared (FTIR) spectroscopy analysis.
Nevertheless, it is an open question whether Asp234 is protonated. Kim et al. pointed out that the loss of photocurrent in the D234N GtACR1 cannot be easily understood if Asp234 is protonated as the influence of the mutation of protonated aspartate to asparagine on the protein function is often small (Kim et al., 2018). GtACR1 crystal structures show that the residues at the Schiff base moiety are highly conserved between GtACR1 and bacteriorhodopsin (BR). Tyr57, Arg82, Tyr185, and Lys216, which are responsible for the low pKa of –2 for the counterion Asp212 in BR (Saito et al., 2012), are fully conserved as Tyr72, Arg94, Tyr207, and Lys238 in GtACR1. Note that the counterion Asp85 in BR, which increases pKa(Asp212) by 6 (Saito et al., 2012), is replaced with Ser97 in GtACR1. This suggests that the pKa of Asp234 in GtACR1 is even lower than the low pKa of –2 for Asp212 in BR. According to Li et al., Tyr72 and Tyr207 donate H-bonds to Asp234 in GtACR1 (Li et al., 2019): this suggests that deprotonated Asp234 is stabilized, decreasing pKa(Asp234), as observed for deprotonated Asp212 in BR. In addition, resonance Raman spectroscopy analysis indicates that the Schiff base Lys238 is also protonated (Yi et al., 2016) as observed in other microbial rhodopsins. The presence of the positively charged Schiff base needs to have an adjacent negative charge (e.g., deprotonated acidic residue) to effectively decrease the energy in GtACR1. To the best of our knowledge, microbial rhodopsins have more than one deprotonated acidic residue adjacent to the retinal Schiff base (e.g., Tsujimura and Ishikita, 2020). This also holds true for channelrhodopsin from Chlamydomonas noctigama (Chrimson) and rhodopsin phosphodiesterase (Rh-PDE), which have both deprotonated and protonated acidic residues near the Schiff base (Vierock et al., 2017; Watari et al., 2019). That is, either Glu68 or Asp234 may be deprotonated in GtACR1. Deprotonation of Asp234 is energetically more favorable than deprotonation of Glu68 in the presence of protonated Schiff base as the electrostatic interaction with Asp234 (3.2 Å) is larger than with Glu68 (5.3 Å). Alternatively, Cl– may exist and act as a counterion, as observed in Cl– pumping rhodopsins (Kim et al., 2016; Kolbe et al., 2000). However, the corresponding electron density is not observed in GtACR1 (Kim et al., 2018; Li et al., 2019). In addition, no spectral changes are reported upon deionization of the sample or exchange from Cl– to SO42– buffer (Sineshchekov et al., 2016). So far, the counterion of GtACR1 remains unknown.
Recently, Dreier et al. proposed that Asp234 is deprotonated in the dark and acts as a counterion according to FTIR measurements and molecular dynamics (MD) simulations (Dreier et al., 2021). The C=O stretching frequencies of 1740 (–)/1732 (+) cm–1 for protonated Asp234 at 77 K observed by Kim et al., 2018 were not observed at 293 K by Dreier et al., 2021. In addition, MD simulations indicated that the H-bond network that involves Tyr72, Tyr207, and Asp234 was stable with deprotonated Asp234 but unstable with protonated Asp234 (Dreier et al., 2021). Indeed, the presence of deprotonated Asp234 was already suggested based on the homology modeling of GtACR2 (Kojima et al., 2018) before the crystal structures of GtACR1 were reported.
GtACR1 undergoes a photocycle including K, L, M, N, and O intermediates (Figure 2; Sineshchekov et al., 2016). The L-state represents the anion conducting state. The L- to M-state transition involves the deprotonation of the Schiff base and the photocurrent decay. The fast photocurrent decay (fast channel closing) corresponds to the M-state formation (i.e., proton release from the Schiff base), and the slow photocurrent decay corresponds to the M-state decay (Sineshchekov et al., 2016; Sineshchekov et al., 2015; Figure 2). Glu68 is likely to accept a proton from the Schiff base upon the M-state formation as a decrease in the accumulation of the M-state was observed in the E68Q GtACR1 (Sineshchekov et al., 2016). However, it remains unclear whether Glu68 is the initial proton acceptor in the wild-type GtACR1. The GtACR1 crystal structures show that Asp234 is closer to the Schiff base (Figure 1; Kim et al., 2018; Li et al., 2019; Li et al., 2021). In addition, the E68Q mutation did not completely inhibit the M-state formation, which indicates that an alternative proton accepter exists (Sineshchekov et al., 2016).
To identify the counterion and clarify the proton-mediated gating mechanism of GtACR1, we investigate the Glu68 and Asp234 mutant proteins (E68Q, E68D, D234N, D234E, and E68Q/D234N) expressed in HEK293 cells. The protonation states are calculated by solving the Poisson–Boltzmann equation and evaluated by conducting MD simulations. Using a quantum mechanical/molecular mechanical (QM/MM) approach, the absorption wavelengths are calculated and the microscopic origin of the wavelength shifts upon the mutations is analyzed.
Results
Protonation states of Glu68 and Asp234
The protonation pattern (Table 1) and pKa values (Table 2 and Table 3) calculated solving the linear Poisson-Boltzmann equation show that Asp234 is deprotonated, whereas Glu68 is protonated in the GtACR1 crystal structures (Kim et al., 2018; Li et al., 2019; Table 1; Supplementary file 1A). The calculated protonation pattern shows that Asp234 is deprotonated in the wild-type GtACR1 even using the MD-generated conformations with protonated Glu68/protonated Asp234 or deprotonated Glu68/protonated Asp234 (Table 1, Supplementary file 1B, Figure 2—figure supplement 1). These results suggest that deprotonation of Asp234, ‘the only residue directly interacting with the protonated Schiff base (Li et al., 2019)’, is a prerequisite to stabilize the protonated Schiff base, as suggested by Dreier et al., 2021.
pKa(Asp234) = –5 (Table 2) is significantly low and even lower than pKa(Asp212) = –2 in BR (Saito et al., 2012). The crystal structures show that the residues at the Schiff base moiety are highly conserved between GtACR1 and BR. Tyr72 and Tyr207 donate H-bonds to each carboxyl O site of Asp234 in GtACR1 (Li et al., 2019; Dreier et al., 2021), whereas Tyr57 and Tyr185 donate H-bonds to each carboxyl O site of Asp212 in BR (Saito et al., 2012). Thus, each tyrosine residue stabilizes the deprotonated state of Asp234, decreasing pKa(Asp234) in GtACR1 by ~3 (Table 2), as observed in BR (Saito et al., 2012). The tendency is also observed for the conserved residue pairs, Arg94/Arg82 and Lys238/Lys216, in GtACR1/BR (Table 2). Asp85, which increases pKa(Asp212) in BR by 6, is replaced with Ser97, which has no influence on pKa(Asp234) in GtACR1 (Table 2). This discrepancy contributes to the low pKa(Asp234) in GtACR1, which is lower than pKa(Asp212) in BR. As far as the original geometry of the GtACR1 crystal structure is analyzed, no residue that increases pKa(Asp234) significantly is identified (Table 2).
Recently, Li et al. reported the GtACR1 conformation (pre-activating state), where Arg94 forms a salt-bridge with Asp234 (Li et al., 2021). The influence of Arg94 on pKa(Asp234) (~3) indicates that the electrostatic link between Arg94 and Asp234 exists even in the ground state (Table 2). It seems possible that the electrostatic interaction between deprotonated Asp234 and channel-gating Arg94 (Table 2) is absent in the D234N GtACR1, leading to the loss of the photocurrent (Kim et al., 2018).
In contrast, pKa(Glu68) is high, 12 (Table 3), which is consistent with the reported protonation state of Glu68 (Yi et al., 2016; Yi et al., 2017; Dreier et al., 2021). The high pKa(Glu68) value can be primarily due to the presence of anionic Asp234 whose deprotonated state is stabilized by the protonated Schiff base (Table 3). Charge neutral Ala53 exists at the corresponding position in BR, which is also consistent with the protonation of Glu68 (Table 3).
Exceptionally, Glu68 is deprotonated only in the D234N GtACR1. The influence of the protonated Schiff base is weaker on Glu68 than on Asp234 (Table 2 and Table 3), which allows to stabilize the putative protonated Glu68 conformation in MD simulations (Table 1). However, the experimentally measured absorption wavelength cannot be reproduced unless Glu68 is deprotonated in the D234N GtACR1 (Table 4, Supplementary file 1C). This is consistent with the absence of the 1708 cm–1 band in the D234N GtACR1, which is assigned to protonated Glu68 in FTIR measurements (Dreier et al., 2021). These results confirm that the presence of a negative charge at the protonated Schiff base moiety is a prerequisite to stabilize the protonated Schiff base, as observed in other microbial rhodopsins (Tsujimura and Ishikita, 2020), including Chrimson and Rh-PDE (Vierock et al., 2017; Watari et al., 2019).
Absorption wavelengths of the Glu68 and Asp234 mutant GtACR1s
The experimentally measured absorption wavelengths for the Glu68 and Asp234 mutant proteins (E68Q, E68D, D234N, D234E, and E68Q/D234N) expressed in HEK293 cells are shown in Figure 3.
In the present study, the experimentally measured absorption wavelengths are the same for the wild-type and D234N GtACR1s (Table 4), which is consistent with results reported previously (Kim et al., 2018; Sineshchekov et al., 2016; Yi et al., 2016). Notably, the calculated absorption wavelengths are also the same for the wild-type and D234N GtACR1s irrespective of deprotonated Asp234 in the wild-type GtACR1 (Table 4, Supplementary file 1C, Figure 3—figure supplements 1 and 2). This can be explained as, in the D234N GtACR1, deprotonated Glu68 moves toward the positively charged Schiff base, which fully substitutes a role of deprotonated Asp234 in stabilizing the positively charged Schiff base (Figure 4a and b, Table 5). A similar conformation of Glu68, which orients toward the Schiff base, was previously reported for the corresponding residues of GtACR2 (Glu64) (Kojima et al., 2018) and channelrhodopsin from Chlamydomonas reinhardtii (Glu90) (Volkov et al., 2017). Thus, the absence of the change in the absorption wavelength upon D234N mutation (Kim et al., 2018) does not necessarily indicate that Asp234 is protonated in the wild-type GtACR1.
The existence of deprotonated Asp234 in the wild-type GtACR1 can also be understood from the absorption wavelength in the D234E GtACR1. The distance between Glu234 and the Schiff base (2.9 Å in MD-generated conformations) in the D234E GtACR1 is shorter than that between Asp234 and the Schiff base (3.4 Å in MD-generated conformations) in the wild-type GtACR1 because glutamate is longer than aspartate (Figure 4a and c, Figure 4—figure supplement 1). The absorption wavelength is short as the electrostatic interaction between the deprotonated counterion and the protonated Schiff base is strong (Tsujimura and Ishikita, 2020). Remarkably, the D234E mutation leads to a decrease in the absorption wavelength (Table 4, Figure 3, Figure 3—figure supplements 1 and 2), which suggests that Asp234 is deprotonated in the wild-type GtACR1. The decrease in the measured absorption wavelength of 10 nm could not be reproduced when we forced Asp234 in wild-type and Glu234 in D234E GtACR1s to protonate (Supplementary file 1D). The electrostatic contributions of charge-neutral protonated Asp234 in wild-type and protonated Glu234 in D234E GtACR1s to the absorption wavelengths are small (Supplementary file 1D).
As far as we are aware, the absorption wavelength of the isolated E68Q/D234N GtACR1 is not reported (e.g., Sineshchekov et al., 2016; Dreier et al., 2021). We successfully isolated a photoactive form of E68Q/D234N GtACR1 using the HEK293 cell expression system, which has been widely used for the functional expression in animal rhodopsins (Kojima et al., 2017; Yamashita et al., 2010). The experimentally measured absorption wavelength in the isolated E68Q/D234N protein is 7 nm longer than that in the wild-type protein (Figure 3), which indicates that Glu68 or Asp234 must be deprotonated in the wild-type GtACR1. As Asp234 is closer to the Schiff base (3.2 Å) than Glu68 (5.3 Å) (Li et al., 2019), it seems more likely that Asp234 is deprotonated in the wild-type GtACR1.
Microbial rhodopsins, including Chrimson and Rh-PDE (Vierock et al., 2017; Watari et al., 2019), have more than one deprotonated acidic residue adjacent to the Schiff base (Tsujimura and Ishikita, 2020). The loss of two acidic residues upon the E68Q/D234N mutation requires an additional negative charge as far as the Schiff base remains protonated. Thus, it seems possible that Cl– exist to stabilize the protonated Schiff base specifically in the E68Q/D234N GtACR1 because the next closest acidic residue, Glu60, is 10 Å away from the Schiff base. The presence of Cl– in the E68Q/D234N GtACR1 is not reported. To investigate the existence of Cl–, isolated E68Q/D234N samples were solubilized in Cl–-free buffer. However, denaturation of the samples did not allow us to conclude the existence of Cl–. In QM/MM calculations and MD simulations, the binding of Cl– at Thr71/Asn234 or Ser97/Lys238 is more stable in the E68Q/D234N GtACR1 than in the wild-type protein (Figure 4—figure supplement 2, Supplementary file 1E). The increase in the calculated absorption wavelength upon the E68Q/D234N mutation (33 nm) is overestimated in the absence of Cl–, whereas the corresponding increase (4 nm) is at the same level as that measured experimentally in the presence of Cl– (Table 4). Thus, Cl– is likely to exist near the protonated Schiff base to compensate for the loss of two acidic residues in the E68Q/D234N GtACR1.
The E68Q mutation does not alter the absorption wavelength (Table 4) as reported previously (Sineshchekov et al., 2016; Yi et al., 2016), thereby suggesting that Glu68 is protonated in the presence of deprotonated Asp234 (e.g., wild-type GtACR1) (Table 1).
In general, blue light-sensitive microbial rhodopsins (e.g., Sensory rhodopsin II and Chlamydomonas channelrhodopsins) show the main absorbance peak with spectral shoulder at shorter wavelength region (e.g., Takahashi et al., 1990). Based on these, it seems likely that the wide band of E68D is due to the existence of the spectral shoulder of this blue-shifted mutant.
Discussion
Our finding of deprotonated Asp234 in the ground state of GtACR1 can explain the following observations: loss of photocurrent upon the D234N mutation (Kim et al., 2018) can be due to loss of Asp234, which is deprotonated in the wild-type GtACR1. It seems possible that the electrostatic interaction between deprotonated Asp234 and channel-gating Arg94 (Li et al., 2021; Table 2) is absent in the D234N GtACR1, leading to loss of the photocurrent (Kim et al., 2018). Intriguingly, MD simulations show that upon the D234N mutation, deprotonated Glu68 reorients toward and interferes with the channel bottle neck (Figure 5). It seems likely that Glu68 acts as a proton acceptor, forming the M-state (i.e., fast channel closing), in the D234N GtACR1 (Sineshchekov et al., 2016), as deprotonated Glu68 is sufficiently close to the protonated Schiff base (Figure 4b). The reorientation of deprotonated negatively charged Glu68 toward the protonated Schiff base and the interference with the channel bottle neck may also explain why the photocurrent owing to the anion conduction is abolished (Kim et al., 2018) irrespective of the accumulation of the M-state (Sineshchekov et al., 2016) (i.e., with deprotonated Schiff base) in the D234N GtACR1. It seems likely that the anion conduction is inhibited in the anion conducting L-state in the D234N GtACR1 because Glu68 already interferes with the channel bottle neck in the ground state. Thus, not only the gating (Arg94) but also the conduction (Glu68) can be inhibited in the D234N GtACR1.
The mechanism is also likely to hold true for the M-state in the wild-type GtACR1, although the MD simulations were conducted based on the dark state structures. As far as we are aware, no intermediate structures of GtACR1 have been reported. The recent time-resolved X-ray free electron laser (XFEL) structures of cation channelrhodopsin C1C2 show that the distance between the Schiff base and Glu129 (Glu68 in GtACR1) remains unaffected during the early part of the photocycle irrespective of the isomerization of the retinal (Figure 5—figure supplement 1; Oda et al., 2021). In C1C2, not only the isomerization of the retinal but also a protein conformational change is required for the conducting-channel formation during the photocycle, as no continuous channel exists in the ground state (Kato et al., 2012). In contrast, a continuous channel spanning through the protein already exists in the ground state of GtACR1 (Figure 5; Li et al., 2019). From the analogy, it seems plausible that the Schiff base interacts electrostatically with Glu68 in the M-state, inhibiting the anion conduction.
The formation of an H-bond between Asn234 and deprotonated Glu68 in the D234N GtACR1 (Figure 5) also suggests that Glu68 accepts the proton from protonated Asp234 in the M-state of the wild-type GtACR1 (Figure 6). Based on the observation of the ground state structure, it seems possible that the proton transfer pathway that proceeds from the protonated Schiff base via deprotonated Asp234 toward protonated Glu68 can also form in the M-state. Protonated Glu68 can accept the proton from transiently protonated Asp234 and simultaneously donates the proton to the adjacent acceptor group in the H-bond network e.g., Asp-L213 in the bacterial photosynthetic reaction center (Sugo et al., 2021).
Then, the absence of Glu68 as a proton acceptor of Asp234 may affect the release of the proton from the Schiff base in the E68Q GtACR1. Indeed, it has been reported that the E68Q mutation affects the channel-gating mechanism (Sineshchekov et al., 2015), leading to a decrease in the M-state (i.e., deprotonated Schiff base) accumulation (Sineshchekov et al., 2016). The conformation of Glu68 as a proton acceptor of Asp234 interferes with the channel bottle neck (Figure 5). This may explain why the M-state formation (i.e., release of the proton from the Schiff base via Asp234 to Glu68) corresponds to the fast channel closing (Figures 2 and 6). This may also explain why the mutation of Glu68 to glutamine leads to a suppression of the fast channel closing at a physiological pH (Sineshchekov et al., 2015).
Conclusions
It was proposed that Asp234 was protonated in the wild-type GtACR1 (Kim et al., 2018; Sineshchekov et al., 2016; Yi et al., 2016; Kandori, 2020) from the following results: (i) the absorption wavelength remains unchanged upon the D234N mutation (Kim et al., 2018; Sineshchekov et al., 2016; Yi et al., 2016); (ii) the C=C stretching frequency of the retinal is not significantly affected upon the D234N mutation in resonance Raman spectroscopy, which implies that the electrostatic interaction between the retinal and the protein environment remains unchanged (Yi et al., 2016); and (iii) the C=O stretching frequencies of 1740 (–)/1732 (+) cm–1 for a protonated carboxylate, which is observed in the wild-type GtACR1, disappear in the D234N GtACR1 at 77 K (Kim et al., 2018). However, the C=O stretching frequencies of 1740 (–)/1732 (+) cm–1 for protonated Asp234 at 77 K were not observed at 293 K by Dreier et al., 2021.
In contrast, the present results show that Asp234 is deprotonated in the wild-type GtACR1, as indicated by the following findings. (i) The E68Q/D234N mutation leads to an increase in the absorption wavelength (Table 4), which indicates that Glu68 or Asp234 is deprotonated in the wild-type GtACR1 (Table 1). (ii) The absorption wavelength in the D234E GtACR1 is shorter than in the wild-type protein (Table 4), which can be explained only by the presence of a deprotonated acidic residue (Table 1, Supplementary file 1D). (iii) The calculated pKa value of –5 for Asp234 is lower than that of –2 for the corresponding residue Asp212 in BR (Saito et al., 2012; Table 2). The significantly low pKa value can be understood as Asp85 in BR, which increases pKa(Asp212) by 6 (Saito et al., 2012), being replaced with Ser97 in GtACR1 (Table 2). The calculated protonation pattern shows that Asp234 is deprotonated in the wild-type GtACR1 even using the MD-generated conformations with protonated Asp234 (Table 1). (iv) Glu68, which is protonated in the wild-type GtACR1, is deprotonated in the D234N GtACR1 (Table 1). If Glu68 remained protonated in the D234N GtACR1, the absorption wavelength would be significantly longer as compared with the wild-type GtACR1 (Table 4). This is consistent with the FTIR measurements, which show that Glu68 is deprotonated in the D234N GtACR1 (Dreier et al., 2021). In any GtACR1, a negative charge needs to be present as far as the Schiff base is protonated. (v) Mutation of deprotonated Asp234 to uncharged asparagine does not alter the calculated absorption wavelength because deprotonated Glu68 reorients and interacts with the Schiff base in the D234N GtACR1, compensating for the change in the charge at the 234 site (Table 4, Table 5, Figure 4). Thus, the absence of changes in the absorption wavelength upon the D234N mutation (Kim et al., 2018; Sineshchekov et al., 2016; Yi et al., 2016) does not serve as a basis of the presence of protonated Asp234 in the wild-type GtACR1. The charge compensation by Glu68 can also explain why the C=C stretching frequency of the retinal, which reflects the electrostatic interaction between the retinal and the protein environment, does not significantly change upon the D234N mutation in resonance Raman spectroscopy (Yi et al., 2016).
The following mechanism can be deduced from the present findings: in D234N GtACR1, anionic Glu68 reorients toward the Schiff base to interact electrostatically. If Asp234 accepts a proton from the Schiff base in the M-state of the wild-type GtACR1, Glu68 is likely to reorient toward the channel, decreasing in the channel radius and inhibiting the anion conduction structurally. Simultaneously, the approach of anionic Glu68 toward the channel pore inhibits anion conduction electrostatically (Figure 5). The mechanism presented here explains why (i) the loss of photocurrent occurs upon the D234N mutation (Kim et al., 2018), (ii) the M-state formation corresponds to the fast channel closing (Sineshchekov et al., 2016), and (iii) the Glu68 to Gln mutation leads to a suppression of the fast channel closing at a physiological pH (Sineshchekov et al., 2015). The formation of the proton transfer pathway in the M-state, which proceeds from the Schiff base via Asp234 and Glu68 toward the protein bulk surface (Figure 6), can explain (iv) the accumulation of the M-state (i.e., deprotonation of the Schiff base) in the wild-type, D234N, and E68Q GtACR1s.
When the properties of a protein (e.g., absorption wavelength) remain unchanged upon the mutation of aspartate to asparagine, one may assume that the aspartate is protonated. However, this does not hold true for the following cases: (i) when the aspartate is adjacent to the focusing site (e.g., forming an H-bond), because the H-bond character (e.g., polarity and pattern) of asparagine is not identical to that of protonated aspartate irrespective of the same net charge; (ii) when another titratable residue exists near the aspartate (e.g., Glu68 in GtACR1) and the protonation states of the two residues are linked. The present example shows that asparagine mutation is not always equivalent to protonated aspartate especially when it is directly involved in the H-bond with the focusing site.
Materials and methods
Coordinates and atomic partial charges
Request a detailed protocolThe atomic coordinates were taken from the X-ray structure of GtACR1 monomer unit ‘A’ (PDB code 6EDQ; Li et al., 2019). All crystal water molecules were included explicitly in calculations if not otherwise specified. During the optimization of hydrogen atom positions with CHARMM (Brooks et al., 1983), the positions of all heavy atoms were fixed, and all titratable groups (e.g., acidic and basic groups) were ionized. The Schiff base was considered protonated. Atomic partial charges of the amino acids and retinal were obtained from the CHARMM22 (MacKerell et al., 1998) parameter set.
Protonation pattern
Request a detailed protocolThe computation was based on the electrostatic continuum model, solving the linear Poisson–Boltzmann equation with the MEAD program (Bashford and Karplus, 1990). The difference in electrostatic energy between the two protonation states, protonated and deprotonated, in a reference model system was calculated using a known experimentally measured pKa value (e.g., 4.0 for Asp; Nozaki and Tanford, 1967). The difference in the pKa value of the protein relative to the reference system was added to the known reference pKa value. The experimentally measured pKa values employed as references were 12.0 for Arg, 4.0 for Asp, 9.5 for Cys, 4.4 for Glu, 10.4 for Lys, 9.6 for Tyr, (Nozaki and Tanford, 1967), and 7.0 and 6.6 for the Nε and Nδ atoms of His, respectively (Tanokura, 1983a; Tanokura, 1983b; Tanokura, 1983c). All other titratable sites were fully equilibrated to the protonation state of the target site during titration. The dielectric constants were set to 4 inside the protein and 80 for water. All water molecules were considered implicitly. All computations were performed at 300 K, pH 7.0, and with an ionic strength of 100 mM. The linear Poisson–Boltzmann equation was solved using a three-step grid-focusing procedure at resolutions of 2.5, 1.0, and 0.3 Å. The ensemble of the protonation patterns was sampled by the Monte Carlo (MC) method with the Karlsberg program (Rabenstein and Knapp, 2001). The MC sampling yielded the probabilities [protonated] and [deprotonated] of the two protonation states of the molecule.
MD simulations
Request a detailed protocolThe GtACR1 assembly was embedded in a lipid bilayer consisting of 258 1-palmitoyl-2-oleyl-sn-glycero-3-phosphocholine (POPC) molecules using CHARMM-GUI (Jo et al., 2008), and soaked in 29070–29072 TIP3P water models, and 5–7 chloride ions were added to neutralize the system using the VMD plug-ins (Humphrey et al., 1996). After structural optimization with position restraints on heavy atoms of the GtACR1 assembly, the system was heated from 0.1 to 300 K over 5.5 ps with time step of 0.01 fs, equilibrated at 300 K for 1 ns with time step of 0.5 fs, and annealed from 300 to 0 K over 5.5 ps with time step of 0.01 fs. The heating and annealing processes to energetically relax the positions of POPC and TIP3 water molecules were performed with time step of 0.01 fs, as done in previous studies (Kurisaki et al., 2015a; Kurisaki et al., 2015b). To avoid the influence of changes in the retinal Schiff base structure on the excitation energy, the position restraints on heavy atoms of side chains were released and MD simulations were performed; the system was heated from 0.1 K to 300 K over 5.5 ps with time step of 0.01 fs and equilibrated at 300 K for 1 ns with time step of 0.5 fs. The system was equilibrated at 300 K for 5 ns with time step of 1.0 fs, and a production run was conducted over 1 ns with 1.0 fs step for sampling of side-chain orientations. 10 conformations were sampled at 0.1 ns intervals during the 1 ns production run. All MD simulations were conducted with the CHARMM22 (MacKerell et al., 1998) force field parameter set using the MD engine NAMD version 2.11 (Phillips et al., 2005). For MD simulations with time step of 1.0 fs, the SHAKE algorithm for hydrogen constraints was employed (Ryckaert et al., 1977). For temperature and pressure control, the Langevin thermostat and piston were used (Feller et al., 1995; Kubo et al., 1991).
POPC molecules have little effect on the calculated pKa values of Asp234 and Glu68 nor the absorption wavelengths (Supplementary file 1F and G). Using the resulting coordinates, the protonation state of the titratable residues was finally determined with the MEAD (Bashford and Karplus, 1990) and Karlsberg (Rabenstein and Knapp, 2001) programs in the absence of POPC molecules.
QM/MM calculations
Request a detailed protocolUsing 10 MD-generated protein conformations, the geometry was optimized using a QM/MM approach in the absence of POPC molecules. The restricted density functional theory (DFT) method was employed with the B3LYP functional and LACVP* basis sets using the QSite (Schrödinger, 2012) program. The QM region was defined as the retinal and Schiff base (Lys238). All atomic coordinates were fully relaxed in the QM region, and the protonation pattern of titratable residues was implemented in the atomic partial charges of the corresponding MM region. In the MM region, the positions of H atoms were optimized using the OPLS2005 force field (Jorgensen et al., 1996), while the positions of the heavy atoms were fixed.
The absorption energy of microbial rhodopsins is highly correlated with the energy difference between highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) (ΔEHOMO-LUMO) or the lowest excitation energy calculated using time-dependent (TD) DFT (ETD-DFT) of the retinal Schiff base (Tsujimura and Ishikita, 2020; Tsujimura et al., 2021). To calculate absorption energies and corresponding wavelengths, the energy levels of HOMO and LUMO and the lowest excitation energies were calculated in the absence of POPC molecules. The absorption energy (Eabs in eV) was calculated using the following equations, which are obtained for wild-type and six mutant GtACR1s (coefficients of determination R2 = 0.93 for Equation 1 and 0.73 for Equation 2; Tsujimura et al., 2021):
The HOMO-LUMO energy gap and the lowest excitation energy were calculated based on 10 MD-generated/QM/MM-optimized protein conformations (see Source data 1 for the atomic coordinates). Empirically, the correlation between the calculated and experimentally measured absorption energies is higher in ΔEHOMO-LUMO than in ETD-DFT among 13 microbial rhodopsin crystal structures (Tsujimura and Ishikita, 2020) and GtACR1 mutants (Tsujimura et al., 2021). In the present study, we analyze the absorption wavelengths of microbial rhodopsin proteins using Equation 1 based on the empirically corrected Eabs (Zhang and Musgrave, 2007) (in eV) or the corresponding wavelength (in nm).
A QM/MM approach utilizing the polarizable continuum model (PCM) method with a dielectric constant of 78 for the bulk region, in which electrostatic and steric effects created by a protein environment were explicitly considered in the presence of bulk water, was employed. In the PCM method, the polarization points were placed on the spheres with a radius of 2.8 Å from the center of each atom to describe possible water molecules in the cavity. The radii of 2.8–3.0 Å from each atom center and the dielectric constant values of ~80 are likely to be optimal to reproduce the excitation energetics, as evaluated for the polarizable QM/MM/PCM approach (Tamura et al., 2020). The TD-DFT method with the B3LYP functional and 6-31G* basis sets was employed using the GAMESS program (Schmidt et al., 1993). The trends in the shifts of absorption wavelength with respect to wild-type GtACR1 remain unchanged when the functional/basis set is replaced (e.g., the CAM-B3LYP functional; Yanai et al., 2004; Figure 3—figure supplement 2).
The electrostatic contribution of the side chain in the MM region to the absorption wavelength of the retinal Schiff base was obtained from the shift in the HOMO-LUMO energy gap upon the removal of the atomic charges of the focusing side chain.
Gene preparation
Request a detailed protocolThe cDNA of GtACR1 (GenBank accession number: KP171708) was optimized for human codon usage and fused to a C-terminal sequence encoding a hexahistidines-tag. The fusion product was inserted into the pCAGGS mammalian expression vector, as previously described (Kojima et al., 2017; Kojima et al., 2020). GtACR1 cDNAs containing mutations were constructed using the In-Fusion Cloning Kit according to the manufacturer’s instructions (Kojima et al., 2017; Kojima et al., 2020).
Protein expression and purification of GtACR1
Request a detailed protocolThe expression plasmids were transfected into HEK293T cells using the calcium-phosphate method (Kojima et al., 2017; Kojima et al., 2020). HEK 293T cells were a gift from Dr. Satoshi Koike (Tokyo Metropolitan Organization for Medical Research). We have confirmed that the identity has been authenticated by STR profiling and the cell lines tested negative for mycoplasma contamination. We have not used any cell lines from the list of commonly misidentified cell lines maintained by the International Cell Line Authentication Committee. After 1-day incubation, all-trans-retinal (final concentration = 5 μM) was added to transfected cells. After another day incubation, the cells were collected by centrifugation (7510 × g for 10 min) at 4°C and suspended in Buffer-A (50 mM HEPES [pH 7.0] and 140 mM NaCl). All-trans-retinal (final concentration = 0.31 μM) was added to the cell suspension to reconstitute the photoactive pigments by shaking rotatory for more than 12 hr at 4°C. Then, the cells were collected by centrifugation (12900 × g for 30 min) at 4°C and suspended in Buffer-A and solubilized in Buffer-B (20 mM HEPES [pH 7.4], 300 mM NaCl, 5% glycerol, and 1% dodecyl maltoside [DDM]). The solubilized fraction was collected by ultracentrifugation (169,800 × g for 20 min) at 4°C, and the supernatant was applied to a Ni2+ affinity column to purify the pigments. After the column was washed with Buffer-C (20 mM HEPES [pH 7.4], 300 mM NaCl, 5% glycerol, 0.02% DDM, and 20 mM imidazole), the pigment was eluted with a linear gradient of imidazole by Buffer-D (20 mM HEPES [pH 7.4], 300 mM NaCl, 5% glycerol, 0.02% DDM, and 1 M imidazole). Purified samples were concentrated by centrifugation using an Amicon Ultra filter (30,000 Mw cut-off; Millipore, USA) and the buffer was exchanged using PD-10 column (GE Healthcare, USA) to Buffer-E (20 mM HEPES [pH 7.4], 300 mM NaCl, 5% glycerol, and 0.02% DDM).
Spectroscopic analysis
Request a detailed protocolAbsorption spectra of the purified proteins were recorded with a UV–visible spectrophotometer (Shimadzu, UV-2450, UV-2600) in Buffer-E. The samples were kept at 15°C using a thermostat.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
An integral equation to describe the solvation of polar molecules in liquid waterThe Journal of Physical Chemistry B 101:7821–7826.https://doi.org/10.1021/jp971083h
-
CHARMM: a program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry 4:187–217.https://doi.org/10.1002/jcc.540040211
-
CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structuresPLOS Computational Biology 8:e1002708.https://doi.org/10.1371/journal.pcbi.1002708
-
Constant pressure molecular dynamics simulation: the Langevin piston methodThe Journal of Chemical Physics 103:4613–4621.https://doi.org/10.1063/1.470648
-
VMD: visual molecular dynamicsJournal of Molecular Graphics 14:33–38.https://doi.org/10.1016/0263-7855(96)00018-5
-
CHARMM-GUI: a web-based graphical user interface for CHARMMJournal of Computational Chemistry 29:1859–1865.https://doi.org/10.1002/jcc.20945
-
Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquidsJournal of the American Chemical Society 118:11225–11236.https://doi.org/10.1021/ja9621760
-
Structure/function study of photoreceptive proteins by FTIR spectroscopyBulletin of the Chemical Society of Japan 93:904–926.https://doi.org/10.1246/bcsj.20200109
-
Green-sensitive, long-lived, step-functional anion channelrhodopsin-2 variant as a high-potential neural silencing toolThe Journal of Physical Chemistry Letters 11:6214–6218.https://doi.org/10.1021/acs.jpclett.0c01406
-
Potential of mean force between two molecular ions in a polar molecular solvent: a study by the three-dimensional reference interaction site modelThe Journal of Physical Chemistry B 103:7942–7957.https://doi.org/10.1021/jp991300+
-
Dewetting of S1-pocket via water channel upon thrombin-substrate association reactionThe Journal of Physical Chemistry. B 119:15807–15812.https://doi.org/10.1021/acs.jpcb.5b09581
-
Toward understanding allosteric activation of thrombin: a conjecture for important roles of unbound Na+ molecules around thrombinThe Journal of Physical Chemistry. B 119:3635–3642.https://doi.org/10.1021/jp510657n
-
All-atom empirical potential for molecular modeling and dynamics studies of proteinsThe Journal of Physical Chemistry. B 102:3586–3616.https://doi.org/10.1021/jp973084f
-
Book2019.01; Chemical Computing Group ULC, 1010 Sherbooke St. West, Suite #910, Montreal, QC, Canada, H3A 2R7Chemical Computing Group ULC.
-
Acid-base titrations in concentrated guanidine hydrochloride. Dissociation constants of the guamidinium ion and of some amino acidsJournal of the American Chemical Society 89:736–742.https://doi.org/10.1021/ja00980a002
-
Scalable molecular dynamics with NAMDJournal of Computational Chemistry 26:1781–1802.https://doi.org/10.1002/jcc.20289
-
Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanesJournal of Computational Physics 23:327–341.https://doi.org/10.1016/0021-9991(77)90098-5
-
Factors that differentiate the H-bond strengths of water near the Schiff bases in bacteriorhodopsin and Anabaena sensory rhodopsinThe Journal of Biological Chemistry 287:34009–34018.https://doi.org/10.1074/jbc.M112.388348
-
General atomic and molecular electronic structure systemJournal of Computational Chemistry 14:1347–1363.https://doi.org/10.1002/jcc.540141112
-
Placevent: an algorithm for prediction of explicit solvent atom distribution-application to HIV-1 protease and F-ATP synthaseJournal of Computational Chemistry 33:1536–1543.https://doi.org/10.1002/jcc.22984
-
1H-NMR study on the tautomerism of the imidazole ring of histidine residues: I. Microscopic pK values and molar ratios of tautomers in histidine-containing peptidesBiochimica et Biophysica Acta (BBA) - Protein Structure and Molecular Enzymology 742:576–585.https://doi.org/10.1016/0167-4838(83)90276-5
-
1H-NMR study on the tautomerism of the imidazole ring of histidine residues: II. Microenvironments of histidine-12 and histidine-119 of bovine pancreatic ribonuclease ABiochimica et Biophysica Acta (BBA) - Protein Structure and Molecular Enzymology 742:586–596.https://doi.org/10.1016/0167-4838(83)90277-7
-
Insights into the protein functions and absorption wavelengths of microbial rhodopsinsThe Journal of Physical Chemistry. B 124:11819–11826.https://doi.org/10.1021/acs.jpcb.0c08910
-
Mechanism of absorption wavelength shifts in anion channelrhodopsin-1 mutantsBiochimica et Biophysica Acta (BBA) - Bioenergetics 1862:148349.https://doi.org/10.1016/j.bbabio.2020.148349
-
Spectroscopic study of the transmembrane domain of a rhodopsin-phosphodiesterase fusion protein from a unicellular eukaryoteThe Journal of Biological Chemistry 294:3432–3443.https://doi.org/10.1074/jbc.RA118.006277
-
A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP)Chemical Physics Letters 393:51–57.https://doi.org/10.1016/j.cplett.2004.06.011
-
Comparison of DFT methods for molecular orbital eigenvalue calculationsThe Journal of Physical Chemistry. A 111:1554–1561.https://doi.org/10.1021/jp061633o
Article and author information
Author details
Funding
Japan Agency for Medical Research and Development (20dm0207060h0004)
- Yuki Sudo
Core Research for Evolutional Science and Technology (JPMJCR1656)
- Yuki Sudo
- Hiroshi Ishikita
Japan Society for the Promotion of Science (JP21K15054)
- Keiichi Kojima
Japan Society for the Promotion of Science (JP20K21482)
- Yuki Sudo
Japan Society for the Promotion of Science (JP21H0040413)
- Yuki Sudo
Japan Society for the Promotion of Science (JP21H0244613)
- Yuki Sudo
Japan Society for the Promotion of Science (JP18H05155)
- Hiroshi Ishikita
Japan Society for the Promotion of Science (JP18H01937)
- Hiroshi Ishikita
Japan Society for the Promotion of Science (JP20H03217)
- Hiroshi Ishikita
Japan Society for the Promotion of Science (JP20H05090)
- Hiroshi Ishikita
University of Tsukuba (Interdisciplinary Computational Science Program in CCS)
- Hiroshi Ishikita
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This research was supported by AMED (20dm0207060h0004) to YS and JST CREST (JPMJCR1656 to YS and HI), JSPS KAKENHI (JP21K15054 to KK, JP20K21482, JP21H0040413, and JP21H0244613 to YS, and JP18H05155, JP18H01937, JP20H03217, and JP20H05090 to HI), and the Interdisciplinary Computational Science Program in CCS, University of Tsukuba.
Copyright
© 2021, Tsujimura 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
-
- 968
- views
-
- 138
- downloads
-
- 7
- 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
Ubiquitination typically involves covalent linking of ubiquitin (Ub) to a lysine residue on a protein substrate. Recently, new facets of this process have emerged, including Ub modification of non-proteinaceous substrates like ADP-ribose by the DELTEX E3 ligase family. Here, we show that the DELTEX family member DTX3L expands this non-proteinaceous substrate repertoire to include single-stranded DNA and RNA. Although the N-terminal region of DTX3L contains single-stranded nucleic acid binding domains and motifs, the minimal catalytically competent fragment comprises the C-terminal RING and DTC domains (RD). DTX3L-RD catalyses ubiquitination of the 3’-end of single-stranded DNA and RNA, as well as double-stranded DNA with a 3’ overhang of two or more nucleotides. This modification is reversibly cleaved by deubiquitinases. NMR and biochemical analyses reveal that the DTC domain binds single-stranded DNA and facilitates the catalysis of Ub transfer from RING-bound E2-conjugated Ub. Our study unveils the direct ubiquitination of nucleic acids by DTX3L, laying the groundwork for understanding its functional implications.
-
- Biochemistry and Chemical Biology
Escherichia coli ClpB and Saccharomyces cerevisiae Hsp104 are AAA+ motor proteins essential for proteome maintenance and thermal tolerance. ClpB and Hsp104 have been proposed to extract a polypeptide from an aggregate and processively translocate the chain through the axial channel of its hexameric ring structure. However, the mechanism of translocation and if this reaction is processive remains disputed. We reported that Hsp104 and ClpB are non-processive on unfolded model substrates. Others have reported that ClpB is able to processively translocate a mechanically unfolded polypeptide chain at rates over 240 amino acids (aa) per second. Here, we report the development of a single turnover stopped-flow fluorescence strategy that reports on processive protein unfolding catalyzed by ClpB. We show that when translocation catalyzed by ClpB is challenged by stably folded protein structure, the motor enzymatically unfolds the substrate at a rate of ~0.9 aa s−1 with a kinetic step-size of ~60 amino acids at sub-saturating [ATP]. We reconcile the apparent controversy by defining enzyme catalyzed protein unfolding and translocation as two distinct reactions with different mechanisms of action. We propose a model where slow unfolding followed by fast translocation represents an important mechanistic feature that allows the motor to rapidly translocate up to the next folded region or rapidly dissociate if no additional fold is encountered.