Molecular basis for functional connectivity between the voltage sensor and the selectivity filter gate in Shaker K+ channels
Abstract
In Shaker K+ channels, the S4-S5 linker couples the voltage sensor (VSD) and pore domain (PD). Another coupling mechanism is revealed using two W434F-containing channels: L361R:W434F and L366H:W434F. In L361R:W434F, W434F affects the L361R VSD seen as a shallower charge-voltage (Q-V) curve that crosses the conductance-voltage (G-V) curve. In L366H:W434F, L366H relieves the W434F effect converting a non-conductive channel in a conductive one. We report a chain of residues connecting the VSD (S4) to the selectivity filter (SF) in the PD of an adjacent subunit as the molecular basis for voltage sensor selectivity filter gate (VS-SF) coupling. Single alanine substitutions in this region (L409A, S411A, S412A, or F433A) are enough to disrupt the VS-SF coupling, shown by the absence of Q-V and G-V crossing in L361R:W434F mutant and by the lack of ionic conduction in the L366H:W434F mutant. This residue chain defines a new coupling between the VSD and the PD in voltage-gated channels.
Introduction
Voltage-gated ion channels play fundamental roles in many physiological processes, such as generating the nerve impulse and regulating neuronal excitability, shaping the pacemaker in the heart, or controlling muscular contractility. Voltage-gated potassium channels (KV) are members of the voltage-gated ion channels superfamily, responsible for the repolarization phase of the action potential, for maximal action potential firing frequency and for simply keeping the membrane potential negative in non-excitable cells. In general, KV channels are assembled as homotetramers. Each subunit monomer contains two functional domains: the voltage sensor domain (VSD), formed by transmembrane segments S1-S4, and the pore domain (PD), composed by transmembrane segments S5, P-loop (a reentrant region in the protein), and S6. The permeation pathway, located in the central axis of the channel, is formed by four PDs flanked by one VSD each.
Presently, there are two tetramerized architectures identified among KV channels: domain-swapped channels (KV1.2 [Long et al., 2005a] and KV7.1 [Sun and MacKinnon, 2017]) and non-domain-swapped channels (EAG1 [Whicher and MacKinnon, 2016], HERG [Wang and MacKinnon, 2017], HCN [Lee and MacKinnon, 2017], BK [Tao et al., 2017], KVAP [Tao and MacKinnon, 2019], and KAT1 [Clark et al., 2020]). In domain-swapped channels, the VSD from one subunit is in close contact with the PD from a neighboring subunit, forming an extensive inter-subunit non-covalent interface. In non-domain-swapped channels, the VSD and the PD from the same subunit forms a similar non-covalent interface and the domains are connected by a shorter version of the S4-S5 linker.
According to a well-studied electromechanical coupling mechanism, henceforth canonical coupling for simplicity, S4 movements driven by changes in the membrane potential lead to changes in S4-S5 linker position, which results in the opening or closing of channel inner gate (bundle crossing), at the intracellular end of the permeation pathway (Kalstrup and Blunck, 2018; Long et al., 2005a; Long et al., 2005b; Lu et al., 2002). The VSD movements, generated by changes in membrane potential, produce gating currents and the charge-voltage (Q-V) curve can be estimated from the current time integral at different voltages. The maximal ionic conductance at each membrane potential serves to estimate the conductance-voltage (G-V) curve. Since VSD movements are the underlying trigger of voltage-dependent inner gate, the gating currents normally precede the pore opening, in voltage and time. Therefore, for a channel that opens as the voltage is made more positive, Q-V curves typically appear to the left of the G-V curves in the voltage axis.
However, the canonical coupling mechanism is not the only functional connection between the VSD and the PD. This can be inferred from an unexpected interaction between L361R, a mutation in S4, and W434F, a mutation in the PD (Carvalho-de-Souza and Bezanilla, 2018). W434F is a well-studied mutation, known to abrogate ionic K+ currents by increasing the speed of C-type inactivation (Perozo et al., 1993; Yang et al., 1997) and it is also widely used to record gating currents in Shaker channels (Bezanilla, 2018). The L361R-W434F interaction was revealed by an unusual crossing between the Q-V curve (from L361R:W434F channels) and the G-V curve (from L361R channels) (Carvalho-de-Souza and Bezanilla, 2018). However, this did not occur when the Q-V curve was estimated from L361R channels, in the absence of W434F and after depleting K+ ions. Therefore, the W434F mutation generated a channel with less voltage sensitivity and with a displaced Q-V curve. This suggests that the VSD and the PD are likely communicating by an alternative pathway, other than through the S4-S5 linker. Such an alternative pathway has been observed in Shaker (Carvalho-de-Souza and Bezanilla, 2019; Conti et al., 2016; Fernández-Mariño et al., 2018) and in KV7.1 (Hou et al., 2017).
In Shaker, a domain-swapped channel (Carvalho-de-Souza and Bezanilla, 2019) showed that this alternative pathway operates independent of the canonical coupling. In fact, it was found to be mediated by the interface of contact between the VSD and the PD from different subunits. They classified this mechanism as noncanonical and demonstrated that it endows voltage dependence to the SF gate. They also hypothesized a series of amino acid residues as the molecular basis for this mechanism. These residues, shown in Figure 1, are arranged parallel to the plane of the membrane and span from the VSD (S4) to an adjacent PD (SF) (Carvalho-de-Souza and Bezanilla, 2019). At the VSD end, depending on the relative position of the S4, the series of residues starts with V363, L366, or V369. Next, comes L409, S411, and S412 in S5, F433, and W434 in the P-loop; at last, Y445 in the SF. From the KV1.2 structure (Chen et al., 2010; Long et al., 2005a), it is possible to appreciate the close contact between these residues.
Here, we studied in detail the role of these residues in the coupling between the voltage sensor and the selectivity filter gate (VS-SF coupling). We assessed it by carrying out individual alanine substitutions at positions L409, S411, S412, and F433 on the L361R:W434F channel. These substitutions were chosen to decrease the van der Waals volume of a particular residue in an attempt to disrupt the functional connection between contiguous residues. All alanine substitutions restored the normal position of the Q-V curve with respect to the G-V curve. We used the same strategy on a mutant channel exhibiting evidence of a VS-SF coupling, the L366H:W434F channel. In this channel, L366H mutation in the VSD unexpectedly restored a fraction of the K+ conduction abrogated by the W434F mutation. Interestingly, all alanine substitutions on L366H:W434F channels were sufficient to decrease the K+ conduction to levels comparable to the W434F phenotype. This indicates that the VS-SF coupling has been disrupted. Together, these results demonstrate that the VS-SF coupling (1) comprises the proposed chain of residues, in domain-swapped channels; (2) is dependent on the volume of these residues; and (3) reveals that the VSD affects the dynamics of the SF gate and that the SF gate affects the movement of the VSD.
Results
Single alanine substitutions disrupt the VS-SF coupling in L361R:W434F channels
One at a time, we introduced the mutations L409A, S411A, S412A, and F433A in a Shaker-IR L361R:W434F background. Shaker-IR L361R:W434F channel studies indicate the mutation W434F actively changes VSD activation dynamics when bearing L361R mutation, evidenced by crossing Q-V and G-V curves. We found that Q-V and G-V curves do not cross each other when any of the alanine substitutions are performed in the presence of the W434F mutation (Figure 2A-F, Table 1). This indicates that the side chain volume in residues L409, S411, S412, and F433 is the base for the communication between the SF and the VSD. The Q-V curves from all mutants containing L361R were satisfactorily fitted with a two-state model equation and the best-fitted values for all curves used in this study are summarized in Table 1. A two-state model equation can satisfactory fit Q-V curves from sensors that appear to populate three or more conformational states during activation, but it does so by underestimating the total gating charges per sensor (Bezanilla and Villalba-Galea, 2013). Another possibility is that the W434F stabilizes the resting state of L361R VSD channel, affecting the VSD activation and/or deactivation, seen as a right shift in the Q-V curve. Altogether, it may be possible that channels bearing L361R mutations have more closed states than WT channels (an issue we did not pursue further in the present study).
Effects by L366H mutation in S4 is consistent with a VS-SF coupling
The single amino acid mutation L366H added to a non-conductive Shaker-IR W434F channel background (L366H:W434F) yields a conductive one. Currents recorded from L366H:W434F channels appear to be highly expressed, as they clearly show gating currents followed by a slower rising current consistent with a K+ conductance 1–2 times larger (Figure 3a–c). This is, however, a much smaller macroscopic conductance to gating currents ratio than in wild-type channels (typically two orders of magnitude or more). These data indicate that, although L366H:W434F channels are conductive, their unitary conductance and/or open probability are dramatically decreased.
Several features of the ionic currents carried by the L366H:W434F channels indicate that they represent K+ moving through the permeation pathway of the PD and not through an eventual omega pore as in some Shaker mutants (Tombola et al., 2005) or being carried by protons as it happens with some other His mutations in the VSD (Starace and Bezanilla, 2004). First, the currents in the L366H:W434F mutant show a decay under sustained voltage-clamp depolarization reminiscent of C-type inactivation, suggesting this process has not been eliminated but rather mitigated (Figure 3a,b). Second, the currents were blocked by 1 mM of 4-aminopyridine (4-AP) (Figure 3c,d), a broad-spectrum K+ channel pore blocker (Kirsch et al., 1993). Third, the currents through L366H:W434F channels are K+-selective as revealed by their reversal potential in different K+ gradients (Figure 3a,b,d). Moreover, the maximal conductance in symmetrical K+ solutions was approximately threefold larger than the conductance in asymmetrical K+ solutions (Figure 3e). It has been reported that under symmetrical K+ concentrations (140 mM), the open probability (Po) of Shaker-IR-W434F is a 100-fold higher as compared to the absence of external K+ (10−7 vs. 10−5) with no expected changes in the single channel conductance (Yang et al., 1997). In a different study, an increase in the number of channels available to open was observed in C-type inactivation enhanced mutants of Shaker (e.g. T449A, T449E, T449K) (López-Barneo et al., 1993). Therefore, we interpret the increased conductance in higher external K+ in L366H:W434F channels as a combination of the increase in the open probability and the increase of the number of channels that open under depolarization, which appears to be a common feature in C-type inactivation enhanced mutants of Shaker.
Moreover, G-V curves from L366H mutants in the absence and in the presence of W434F mutation have very similar voltage dependency and voltage sensitivity (Figure 3f, Table 1). This observation further reinforces the notion that the PD of L366H:W434F channels is conductive and suggests that the canonical coupling underlying channel gating seems to be unaltered in this double mutant. Q-V curves measured in L366H channels with V478W mutation, that prevents inner gate (bundle crossing) opening (Kitaguchi et al., 2004) or W434F (after conductance blockade by 4-AP), are also practically identical (Figure 3g – Table 1). This implies that the W434F mutation does not affect VSD movements in the particular case of L366H-mutant VSD.
Single alanine substitutions in the chained residues eliminate VS-SF coupling in L366H:W434F channels
If L366H and W434F are interacting and this interaction is carried out by the same coupling pathway as demonstrated for L361R and W434F case (Figures 1 and 2), then the ionic K+ currents through L366H:W434F channels should be decreased to levels comparable to the W434F phenotype when a similar alanine mutation strategy is used. Consistently, mutating L409, S411, or S412 to alanine in L366H:W434F channels, the K+ currents were greatly depressed leaving only gating currents to be observed (Figure 4). Once again, the magnitude of gating current recorded attests to the high level of protein expression. Peak of ionic currents from all triple mutants were remarkably small (to ~0.3 µA at 60 mV), indicating L409, S411, and S412 are also involved in the VS-SF coupling in L366H:W434F channels. Due to low expression levels, we were unable to record gating currents from L366H:F433A:W434F triple-mutant channels.
Effects of single alanine mutations on C-type inactivation and VSD to pore conduction coupling
Even though alanine substitutions of amino acids within the connecting chain have been an effective strategy to delineate the VS-SF coupling, it is important to understand the functional consequences of such mutations in regular Shaker-IR-WT channels. We assessed the influence of those alanine mutations on the slow inactivation process and on Q-V and G-V curves (Figure 5). The weighted inactivation time constants (τInac) for WT, L409A, S411A, S412A, and F433A were assessed by fitting a double exponential equation to the decaying phase of K+ currents elicited by long (19 s) depolarized voltage pulses (Figure 5a–e). While the τInac was essentially identical for F433A and WT, it was considerably slower for L409A, S411A, and S412A (Figure 5f). The fraction of the fast component of the inactivation for L409A and S411A was smaller when compared to WT, but similar for S412A and F433A (Figure 5g).
Some of the alanine substitutions (L409A, S411A, and F433A) led to significant right shifts both in the G-V curve and in the available current after several different inactivating voltages pulses (Inac-V curve) (see Figure 5h–l and Table 1). The corresponding Q-V curves from those mutants, on the other hand, were not right-shifted, suggesting that those replaced residues exert a significant role in the VSD to pore conduction coupling in regular channels. We evaluated the change in the energy involved in both the activation and the inactivation processes, as induced by L409A, S411A, S412A, and F433A. To that end, we evaluated (see 'Materials and methods'): (1) the second component of the Q-V curves, taken from fitting a three-state model (Equation 3, see 'Materials and methods') (EQ2, Figure 5m); (2) the conductance (EG, Figure 5n); and (3) the C-type inactivation process (EInac, Figure 5o). We evaluated the second component of the charge movement because this transition is associated with the opening of the pore (Ledwell and Aldrich, 1999). A threshold of 1 kcal/mol, relative to WT, was set as a threshold for significance when considering the energetics of each process. Among all alanine substitutions, only F433A led to a significant change in the energy of the second transition of the Q-V curve (~1.7 kcal/mol – Figure 5m). For the energy associated to the conductance process, only mutation S412A did not lead to a significant change (Figure 5n). Whereas only mutation S411A did not lead to a significant change in the energy related to the inactivation process (Figure 5o).
Altogether, from these results we learned that, apart from F433A, alanine substitutions in the VS-SF coupling have little or no effect in the dynamics of WT VSD. Nevertheless, those same residues are crucially important for the channel VSD to pore conduction coupling. It involves both activation and inactivation processes even though the overall effects were small. Notably, the available current at very positive voltages was considerably increased when mutations S412A and F433A are present in the channels (Figure 5h,k and l), showing the importance of the residues to stabilize the inactivation conformation at the SF.
Interaction between members of the chain does not depend on W434F mutation
We further explored the role of the VS-SF coupling by assessing a possible interaction between residues S411 and F433. The rationale for using those residues lays on their proximity in the protein, according to molecular models available. Another motivation to study the interaction between these two residues is that they have not been previously implicated with a noncanonical electromechanical coupling as shown for S412 (Fernández-Mariño et al., 2018). These two residues also do not affect the inactivation of the channels as W434F mutation does (Perozo et al., 1993). To record gating currents by a mechanism that does not involve the SF and W434F, we used the V478W mutation at the bundle crossing region. This mutation prevents ionic conduction (Kitaguchi et al., 2004; Pau et al., 2017) and does not affect the movement of the VSD (Kitaguchi et al., 2004). Therefore, we can assess the interaction between S411 and F433, in the absence of W434F. We used the thermodynamic mutant cycle and the generalized interaction-energy analysis to calculate this interaction (Chowdhury et al., 2014; Chowdhury and Chanda, 2012; Figure 6). This method uses the median voltage (VMedian) from the Q-V curves to estimate the free energy of activation. The VMedian values calculated are shown in Table 2. The interaction between S411 and F433 is +2.00 ± 0.76 kcal/mol and is above the cutoff value of 1.8 kcal/mol, as previously established for this kind of analysis (Chowdhury et al., 2014; Fernández-Mariño et al., 2018). Therefore, this data supports that the connecting chain operates in the WT channels, independent of the presence of W434F mutation.
Discussion
Recently, several lines of evidence have pointed to the existence of different coupling mechanisms between the VSD and the PD in Shaker channels (Fernández-Mariño et al., 2018; Conti et al., 2016; Carvalho-de-Souza and Bezanilla, 2019; Petitjean et al., 2015). These mechanisms extend beyond the canonical electromechanical S4-S5 linker-based coupling. In this study, we demonstrate that a chain of amino acid residues is determinant for an alternative coupling in Shaker. Three residues in the S5 segment (L409, S411, S412) and two in the P-loop (F433, W434) are key for this functional VS-SF coupling. Crucially, the VSD and the PD involved in this chain are conformationally coupled through adjacent subunits. Yet, we must consider the fact that these aforementioned residues may not be the only ones involved in these processes, as it is possible that other paths may also exist (Conti et al., 2016; Lee et al., 2009).
We studied two functionally distinct abnormal W434F construct channels: L361R:W434F and L366H:W434F. In the first case, W434F affects the movement of L361R VSD seen as a shallower Q-V that crosses the G-V curve from channels not bearing W434F mutation. In the second case, L366H partially relieves the inactivation promoted by W434F, turning a non-conductive channel into a conductive one. The observation of these cases is quite remarkable because it shows the bidirectionality of the putative VS-SF coupling between the VSD and the PD. Irrespective of which abnormal W434F construct we used, single alanine substitutions to any of the residues of the chain are sufficient to disrupt or at least mitigate the abnormal phenotype of the W434F-bearing channels.
It is important to notice that all our observations showing an exacerbated VS-SF coupling were done in the presence of the W434F mutation. Under those conditions, one could argue that this mechanism might operate only after the channel populates the energy well equivalent to a C-type inactivated channel (the W434F phenotype). Our data suggest that the mechanism is in fact always present (Figures 5 and 6) and it is independent of the presence of W434F mutation; however, the W434F mutation makes it more conspicuous because this mutation significantly enhances C-type inactivation (Yang et al., 1997). Therefore, the roles of these chained residues on the overall VSD to pore opening coupling mechanism in regular Shaker channels are validated. Indeed, these results when combined demonstrate that there is a detectable coupling path between the VSD and the PD, as outlined in Figure 1.
Recently, a different study has shown that F244C, a mutation in S1 segment of Shaker channels, also produces a similar relief in the W434F-induced inactivation: the double mutant F244C:W434F also conducts K+ and the current is blocked by 4-AP (Petitjean et al., 2015). F244, one of the residues forming the hydrophobic plug in Shaker channels (Lacroix et al., 2014), is in close proximity to residue S412 (Chen et al., 2010) from a PD of a neighboring subunit. Our data is consistent with these results, where mutations outside the PD, that is, L361R and L366H are able to interact with the W434F. Thus, under our proposal, the interaction between F244C and W434F would take place via the VS-SF coupling as well. Consistently, the F244C homolog in KV1.1, F184C, causes episodic ataxia in human patients (Browne et al., 1995). Mutations in human KV channels at the S412 equivalent have neurological (in KV1.1 [Lee et al., 2004] and KV7.2 [Dedek et al., 2003]) and cardiac (in KV7.1 [Liu et al., 2002]) consequences in humans.
In Shaker channel, the residue I470 is a key participant in the coupling between the activation gate, at the bundle crossing region, and C-type inactivation (Peters et al., 2013). It has been proposed that I470 interacts with T442, in the signature sequence of K+-selective channels, TVGYG, at the SF (Labro et al., 2018). This mechanism constitutes an indirect pathway that couples the voltage sensor with the inner gate and the SF. Several studies have shown that mutations in the outer mouth of the pore, S5 segment and P-loop, affect C-type inactivation, highlighting the importance of these motifs for the inactivation process (De Biasi et al., 1993; Larsson and Elinder, 2000; Molina et al., 1997; Ortega-Sáenz et al., 2000; Perozo et al., 1993; Yang et al., 1997). These studies and our data demonstrate the existence of a direct interaction between the voltage sensor and the SF mediated by the mechanism proposed here.
The above discussion prompts two important questions: what is the contribution of the VS-SF coupling to the overall VSD to pore opening coupling? Furthermore, if the canonical mechanism is disrupted, can the channel be gated by the VS-SF coupling? Presently, we do not have data to support a definite answer for these questions, but it seems that the canonical coupling mechanism is the main mechanism whereby the channels’ pore is gated by the VSD. Studies of channels with a stabilized open pore conformation may help evaluate the physiological role of the VS-SF coupling. Inter-subunit metal bridges between residues V476C and H486 tend to stabilize the channels in the open conformation, displaying K+ currents without apparent activation kinetics, as expected in the absence of a voltage-dependent gating process (Holmgren et al., 1998). A recent study, using T449A:V476C Shaker background and metal bridges, has concluded that the closure of the activation gate (where the canonical coupling mechanism operates) is essential for the recovery from C-type inactivation (Szanto et al., 2020). This suggests that the VS-SF coupling, by itself, is not sufficient to gate the channel’s inner gate and, consequently, the channel. Furthermore, mutations in the bundle crossing region have shown that the pore stays open and does not close even when very negative potentials are applied (Kitaguchi et al., 2004). Regardless, an experimental strategy is yet to be developed to determine the role of VS-SF coupling in the absence of the canonical one.
It is reasonable to believe that the VS-SF coupling cannot, by itself, transfer enough energy to the PD in order to gate the channel, pointing to the notion that it is an accessory or modulatory gating mechanism. Indeed, this alternative mechanism might be more relevant in non-swapped domain channels such as from ether-a-go-go (EAG) K+ channels family (EAG and hERG), which can open even without covalent interaction between the VSD and the PD (without the classical canonical coupling) (Lörinczi et al., 2015). This result suggests that the multiple non-covalent interactions between amino acid residues from the VSD and the PD can endow voltage dependence to the open probability of the PD. Furthermore, in hERG a mutation at W568 has been previously proposed to be involved in a pathway that couples activation and inactivation and when mutated to Leu (W568L), the channel does not exhibit C-type inactivation (Ferrer et al., 2011). Interestingly, residue W568 in hERG is homolog to residue S412 in Shaker demonstrating the physiological relevance of such alternative coupling mechanism in other members of K+ channels family.
In summary, we have shown that residues L409, S411, S412, F433, and W434F comprise the molecular chain of residues that directly connects the voltage sensor to the SF gate, in Shaker K+ channels. Our data also corroborated and clarified the idea of an alternative coupling mechanism. It is also possible that this mechanism is not only important in domain-swapped and non-swapped K+ channels, but also may play a role in other members of the voltage-gated ion channels superfamily, for example, NaV channels.
Materials and methods
Site-directed mutagenesis
Request a detailed protocolShaker zH4 K+ channel with fast inactivation removed (Hoshi et al., 1990), ∆6–46, and cloned into pBSTA vector was used. Mutations were performed using Quick-change II technology (Stratagene, La Jolla, CA), together with custom primers from Integrated DNA Technologies (Integrated DNA Technologies, Inc, Coralville, IA). Shaker cDNA and its mutants were sequenced, linearized by endonuclease NotI (New England Biolabs, Ipswich, MA), and cleaned up with a NucleoSpin Gel and PCR Clean-up kit (Macherey-Nagel, Bethlehem, PA). In vitro transcription kits were used to transcribe cDNA and generate cRNA (T7 RNA expression kit; Ambion Invitrogen, Thermo Fisher Scientific, Waltham, MA).
Oocytes preparation
Request a detailed protocolOocytes were harvested from Xenopus laevis in accordance with experimental protocols #71475 approved by the University of Chicago Institutional Animal Care and Use Committee (IACUC). Following the follicular membrane digestion by collagenase, oocytes were incubated in standard oocytes solution (SOS) containing in mM: 96 NaCl, 2 KCl, 1.8 CaCl2, 1 MgCl2, 0.1 EDTA, 10 HEPES, and pH set to 7.4 with NaOH. SOS was supplemented with 50 µg/ml gentamycin to avoid contamination during incubation. After 6-24 hr of harvesting, defolliculated oocytes stage V-VI were injected with cRNA (5-100 ng diluted in 50 nl of RNAse free water) and incubated for 1-3 days at 12°C or 18°C prior to recording.
Electrophysiological recordings
Request a detailed protocolIonic and gating currents were recorded using cut-open oocyte voltage-clamp method (Stefani and Bezanilla, 1998). Voltage-measuring pipettes were pulled using a horizontal puller (P-87 Model, Sutter Instruments, Novato, CA) to a tip resistance ranging from 0.2 to 0.5 MΩ. Currents were acquired by a setup comprising a Dagan CA-1B amplifier (Dagan, Minneapolis, MN) with a built-in low-pass four-pole Bessel filter for a cutoff frequency of 20–50 kHz. Using a 16-bit A/D converter (USB-1604, Measurement Computing, Norton, MA) for acquisition and controlled by an in-house software (GPatch64MC), data were sampled at 1 MHz, digitally filtered at Nyquist frequency and decimated for a storage acquisition rate of 50–200 kHz. Capacitive transient currents were compensated using a dedicated circuit. The voltage clamp was also controlled using the same in-house software and the 16-bit D/A converter of the same USB-1604 device.
For ionic current measurements, the external solution was composed by (mM): K-methanesulfonic acid (MES) 12, Ca-MES 2, HEPES 10, EDTA 0.1, N-methyl-D-glucamine (NMDG)-MES 108, pH 7.4. For 120 mM external K+ solution, we replace NMDG with K+. The cut-open oocyte internal solution was composed by (mM): K-MES 120, EGTA 2, HEPES 10, pH 7.4. For gating current recordings, we replaced K+ by NMDG in both internal and external solutions. A conditioning hyperpolarized pre-pulse was used with different voltages to ensure that most of the channels were closed prior to the depolarizing pulses. The holding potential was set to −80 or −100 mV according to the mutant. Recordings were performed at room temperature (~17–18°C). All chemicals used were purchased from Sigma-Aldrich (St. Louis, MO).
Data analysis
Request a detailed protocolThe peaks of macroscopic K+ currents were converted into conductance (), using the equation below:
where is the K+ current, is the membrane voltage, is the gas constant, is the temperature in Kelvin, is the Faraday constant. and are the intracellular and extracellular K+ concentrations, respectively.
Each value was normalized by the maximum of the experiment for averaging among experiments and plotting against in order to obtain a G-V curve. The G-V curves were fitted using a two-state model with the following equation:
where is the apparent charge of the transition expressed in units of elementary charge () and is the voltage for 50% of the maximal conductance.
The charge () was calculated from the integral in time of the gating currents. They were normalized, averaged, and plotted against to generate a Q-V curve. The Q-V curves were fitted using a three-state model (Lacroix et al., 2012):
where is proportional to the number of VSDs, and are apparent charges for the first and the second lumped sequential transitions of the VSD movement, respectively, and are the membrane voltages that split the VSD into two states equally during the first and second transitions, respectively.
In some cases, a two-state model was used to best fit the Q-V curves, as shown below:
where is the apparent charge of the transition expressed in units of elementary charge () and is the voltage for 50% of the maximal charge.
Inactivation time constants were taken from the decaying phase of K+ currents after sustained activation by depolarizing voltage pulses. They were fitted using a double exponential shown below:
where and are the amplitudes from the first and the second exponential, respectively, and are the time constants for the first and second exponentials, respectively, and is a baseline adjustment.
For simplicity and analysis purposes, we considered a weighted time constant, calculated by the following equation:
where is the weighted time constant.
In order to build Inac-voltage (Inac-V) curves for each mutant, the peak K+ currents from five to eight oocytes were taken from the test pulse as shown in Figure 5A. They were normalized by their maxima, averaged and plotted against . A two-state model was used to fit the Inac-V curves by the equation below:
where is the apparent charge of the transition expressed in units of elementary charge () and is the voltage dependence of the inactivation process. A1 and A2 are the minimum and the maximum value from the data, respectively.
We calculated the free energy () of the studied processes (second component of charge movement [Q2], conductance activation [G], and inactivation [Inac]) by using and from the Q-V, G-V, and Inac curves. We used the following equation to calculate the free energies:
where Gprocess is the free energy involved in the process, zprocess is the apparent charge involved in the process, F is the Faraday constant, and Vprocess is the voltage necessary to complete half of the total process (the V1/2).
The shift of the studied process was calculated using the following equation:
The standard error (SE) associated to the energies were calculated using the following equation:
To test the interaction between two residues, we used the thermodynamic mutant cycle, calculating the free energy using the generalized interaction-energy analysis that uses the from the Q–V curves (Chowdhury and Chanda, 2012; Chowdhury et al., 2014). The free energy of activation (G) for a channel is calculated using the equation:
The uncertainty associated in G (δG) is:
The number of charges (z) was considered the same for all channels: 13.6 elementary charges (Schoppa et al., 1992; Aggarwal and MacKinnon, 1996; Seoh et al., 1996).
The interaction free energy between residues S411 and F433 is:
The subscripts refer to channel used to calculate the free energy. The uncertainty associated to is given by:
, , , and are, respectively, the standard error of the mean associated with the measurement of for V478W, S411A:V478W, F433A:V478W, and S411A:F433A:V478W channels.
Matlab (The MathWorks, Inc, Natick, MA) and Origin9.0 (Origin Lab Corporation, Northampton, MA) were used for calculating G-Vs, Q-Vs, plotting and fitting the data. Data is shown as mean ± SEM.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
-
Fig 2BID JGP. Nonsensing residues in S3-S4 linker's C terminus affect the voltage sensor set point in K+ channels.
References
-
The gating charge should not be estimated by fitting a two-state model to a Q-V curveJournal of General Physiology 142:575–578.https://doi.org/10.1085/jgp.201311056
-
Identification of two new KCNA1 mutations in episodic ataxia/myokymia familiesHuman Molecular Genetics 4:1671–1672.https://doi.org/10.1093/hmg/4.9.1671
-
Nonsensing residues in S3–S4 linker’s C terminus affect the voltage sensor set point in K+ channelsJournal of General Physiology 150:307–321.https://doi.org/10.1085/jgp.201711882
-
A self-consistent approach for determining pairwise interactions that underlie channel activationJournal of General Physiology 144:441–455.https://doi.org/10.1085/jgp.201411184
-
Estimating the voltage-dependent free energy change of ion channels using the median voltage for activationJournal of General Physiology 139:3–17.https://doi.org/10.1085/jgp.201110722
-
Inactivation determined by a single site in K+ poresPflügers Archiv 422:354–363.https://doi.org/10.1007/BF00374291
-
Gating interaction maps reveal a noncanonical electromechanical coupling mode in the shaker K+ channelNature Structural & Molecular Biology 25:320–326.https://doi.org/10.1038/s41594-018-0047-3
-
Molecular coupling in the human ether-a-go-go-related gene-1 (hERG1) K+ channel inactivation pathwayJournal of Biological Chemistry 286:39091–39099.https://doi.org/10.1074/jbc.M111.292060
-
Stabilizing the closed S6 gate in the shaker kv channel through modification of a hydrophobic sealJournal of General Physiology 124:319–332.https://doi.org/10.1085/jgp.200409098
-
Intermediate state trapping of a voltage sensorJournal of General Physiology 140:635–652.https://doi.org/10.1085/jgp.201210827
-
Mutations in the S4 region isolate the final voltage-dependent cooperative step in potassium channel activationJournal of General Physiology 113:389–414.https://doi.org/10.1085/jgp.113.3.389
-
Effects of external cations and mutations in the pore region on C-type inactivation of shaker potassium channelsReceptors & Channels 1:61–71.
-
Coupling between voltage sensors and activation gate in voltage-gated K+ channelsJournal of General Physiology 120:663–676.https://doi.org/10.1085/jgp.20028696
-
Collapse of conductance is prevented by a glutamate residue conserved in voltage-dependent K(+) channelsJournal of General Physiology 116:181–190.https://doi.org/10.1085/jgp.116.2.181
-
Crystal structure of an inactivated mutant mammalian voltage-gated K+ channelNature Structural & Molecular Biology 24:857–865.https://doi.org/10.1038/nsmb.3457
-
Cut-open oocyte voltage-clamp techniqueMethods in Enzymology 293:300–318.https://doi.org/10.1016/s0076-6879(98)93020-8
-
The activation gate controls steady-state inactivation and recovery from inactivation in shakerJournal of General Physiology 152:e202012591.https://doi.org/10.1085/jgp.202012591
-
An emerging consensus on voltage-dependent gating from computational modeling and molecular dynamics simulationsJournal of General Physiology 140:587–594.https://doi.org/10.1085/jgp.201210873
-
How does the W434F mutation block current in shaker potassium channels?Journal of General Physiology 109:779–789.https://doi.org/10.1085/jgp.109.6.779
Article and author information
Author details
Funding
National Institutes of Health (R01-GM030376)
- Francisco Bezanilla
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are deeply in debt with Drs Miguel Holmgren and Eduardo Perozo, and with Fraol Galan for their careful reading, comments, and suggestions on the manuscript. We thank Li Tang for helping with molecular biology. Supported by NIH grant R01-GM030376.
Ethics
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All of the animals were handled according to approved Institution of Animal Care and Use Committee (IACUC) protocols (#71745) of the University of Chicago.
Copyright
© 2021, Bassetto 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
-
- 1,350
- views
-
- 268
- downloads
-
- 21
- 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
-
- Neuroscience
Time estimation is an essential prerequisite underlying various cognitive functions. Previous studies identified ‘sequential firing’ and ‘activity ramps’ as the primary neuron activity patterns in the medial frontal cortex (mPFC) that could convey information regarding time. However, the relationship between these patterns and the timing behavior has not been fully understood. In this study, we utilized in vivo calcium imaging of mPFC in rats performing a timing task. We observed cells that showed selective activation at trial start, end, or during the timing interval. By aligning long-term time-lapse datasets, we discovered that sequential patterns of time coding were stable over weeks, while cells coding for trial start or end showed constant dynamism. Furthermore, with a novel behavior design that allowed the animal to determine individual trial interval, we were able to demonstrate that real-time adjustment in the sequence procession speed closely tracked the trial-to-trial interval variations. And errors in the rats’ timing behavior can be primarily attributed to the premature ending of the time sequence. Together, our data suggest that sequential activity maybe a stable neural substrate that represents time under physiological conditions. Furthermore, our results imply the existence of a unique cell type in the mPFC that participates in the time-related sequences. Future characterization of this cell type could provide important insights in the neural mechanism of timing and related cognitive functions.
-
- Neuroscience
Granule cells of the cerebellum make up to 175,000 excitatory synapses on a single Purkinje cell, encoding the wide variety of information from the mossy fibre inputs into the cerebellar cortex. The granule cell axon is made of an ascending portion and a long parallel fibre extending at right angles, an architecture suggesting that synapses formed by the two segments of the axon could encode different information. There are controversial indications that ascending axon (AA) and parallel fibre (PF) synapse properties and modalities of plasticity are different. We tested the hypothesis that AA and PF synapses encode different information, and that the association of these distinct inputs to Purkinje cells might be relevant to the circuit and trigger plasticity, similar to the coincident activation of PF and climbing fibre inputs. Here, by recording synaptic currents in Purkinje cells from either proximal or distal granule cells (mostly AA and PF synapses, respectively), we describe a new form of associative plasticity between these two distinct granule cell inputs. We show for the first time that synchronous AA and PF repetitive train stimulation, with inhibition intact, triggers long-term potentiation (LTP) at AA synapses specifically. Furthermore, the timing of the presentation of the two inputs controls the outcome of plasticity and induction requires NMDAR and mGluR1 activation. The long length of the PFs allows us to preferentially activate the two inputs independently, and despite a lack of morphological reconstruction of the connections, these observations reinforce the suggestion that AA and PF synapses have different coding capabilities and plasticity that is associative, enabling effective association of information transmitted via granule cells.