|
||
ARTICLE |
Correspondence to Brad S. Rothberg: rothberg{at}uthscsa.edu
|
|
|---|
and ß subunits. In this study, we used perturbation mutagenesis to analyze the role of S0 in channel gating. Single residues in the S0 region of the BK channel were substituted with tryptophan to give a large change in side chain volume; native tryptophans in S0 were substituted with alanine. The effects of the mutations on voltage- and Ca2+-dependent gating were quantified using patch-clamp electrophysiology. Three of the S0 mutants (F25W, L26W, and S29W) showed especially large shifts in their conductance–voltage (G-V) relations along the voltage axis compared to wild type. The G-V shifts for these mutants persisted at nominally 0 Ca2+, suggesting that these effects cannot arise simply from altered Ca2+ sensitivity. The basal open probabilities for these mutants at hyperpolarized voltages (where voltage sensor activation is minimal) were similar to wild type, suggesting that these mutations may primarily perturb voltage sensor function. Further analysis using the dual allosteric model for BK channel gating showed that the major effects of the F25W, L26W, and S29W mutations could be accounted for primarily by decreasing the equilibrium constant for voltage sensor movement. We conclude that S0 may make functional contact with other transmembrane regions of the BK channel to modulate the equilibrium between resting and active states of the channel's voltage sensor.
Dr. Fan's present address is Department of Biomedical Engineer-ing, The Fourth Military Medical University, Xian, China.
| INTRODUCTION |
|---|
|
|
|---|
BK channels are unique, however, in that they contain a Ca2+ sensor that is formed by the cytoplasmic carboxyl ends of the four channel subunits (Schreiber et al., 1999
; Bian et al., 2001
; Bao et al., 2002
; Bao et al., 2004
). Also, BK channels contain an additional transmembrane region, S0, at the amino end of each subunit (Fig. 1) (Wallner et al., 1996
; Meera et al., 1997
).
This segment contains several residues that are highly conserved among BK channel orthologues. This sequence conservation may indicate a role for these side chains in channel structure and function, yet the function of the S0 segment has not yet been studied in detail. It also seems that the S0 region may form an important component of the interacting region between the BK channel and its ß subunits (Wallner et al., 1996
). Thus it will be important to learn the role of the S0 segment in gating if we are to understand the molecular basis of this interaction.
|
-helical S0 segment, then Phe-25, Leu-26, and Ser-29 would be localized to a contiguous patch on the same face of the helix. Together our results are consistent with several S0 residues forming a tightly focused functional contact with other parts of the BK channel to modulate the equilibrium between resting and active states of the channel's transmembrane voltage-sensing domains. | MATERIALS AND METHODS |
|---|
|
|
|---|
Solutions and Electrophysiological Recordings
All experiments were performed with excised inside-out patches from identified transfected cells 1–2 d after transfection. Experiments were done at room temperature (22–24°C). Currents were low-pass filtered at 3–5 kHz and digitized at 20 kHz. Except where noted, solutions bathing both sides of the membrane contained 160 mM KCl and 10 mM HEPES (pH 7.4). For most of the data reported below, the solution at the cytoplasmic face of the patch additionally contained either 2 mM N-(2-hydroxyethyl)-ethylenediamine-triacetic acid (HEDTA), to obtain free [Ca2+]'s ranging from 1 to 15 µM, calculated using the program MaxChelator (www.stanford.edu/
cpatton/maxc.html) for pH 7.4, T = 22°C, and ionic strength = 0.16. For some patches, additional data were obtained at 110 µM free [Ca2+], using 2 mM nitrilotriacetic acid (NTA) to buffer the [Ca2+], or at nominally 0 Ca2+, using 2 mM EGTA with no added Ca2+. CaCl2 was added to these intracellular solutions to bring the free [Ca2+] to the indicated levels. The free [Ca2+] values reported in this paper were measured from recording solutions using a Ca2+-sensitive electrode (Orion), based on a set of standard solutions that contained 160 mM KCl, 10 mM HEPES at pH 7.4, and added CaCl2 ranging up to 10 mM (in the absence of Ca2+ buffers). From these measurements, we estimated the contaminating [Ca2+] with no added CaCl2 to be
20 µM, originating mainly from the KCl and HEPES salts.
For recordings in nominally 0 Ca2+i, the bath solution at the cytoplasmic face of the patch was composed of 50 mM KCl, 10 mM HEPES, and 2 mM EGTA, pH 7.4, with no added Ca2+. Under these conditions, we estimate that the free [Ca2+] should be no greater than 0.0004 µM (assuming 20 µM Ca2+ present as a contaminant), calculated with MaxChelator. This concentration is outside of the range where Ca2+ has any significant effect on BK channel gating (Rothberg and Magleby, 2000
).
Voltage protocols for each dataset were typically repeated at least five times for each patch, and current traces were averaged online. To minimize voltage errors due to series resistance, we analyzed only recordings in which the maximal current was <4 nA. We estimate that the maximal voltage error contributed by series resistance for these recordings was
6 mV. Analyzed data are presented without correction for series resistance.
Data Analysis
We quantified the effects of mutations on gating by comparing parameters obtained from G-V relations after fitting data with a Boltzmann function,
![]() |
To additionally account for the weakly voltage-dependent gating that occurs with voltage sensors deactivated, data in Fig. 4 were fitted with a modified Boltzmann function:
![]() | (1) |
For wild-type data and selected mutants, gating was analyzed further by global fitting of G-V and PO vs. V data obtained at several different [Ca2+] with a dual allosteric model (Horrigan and Aldrich, 2002
):
![]() | (2) |
Data predicted by the dual allosteric model were calculated from a starting set of kinetic parameters; model predictions were then compared with the experimental data, and kinetic parameters were optimized using an iterative
2 minimization routine, with a modified Simplex search algorithm (Magleby and Weiss, 1990
; Rothberg and Magleby, 1998
). Automated fitting was combined with manual parameter adjustment in some cases, in order to better account for low open probability data. Model parameters that were not well constrained by experimental data obtained for these studies were fixed at or near values determined in previous studies on BK channels (Horrigan and Aldrich, 2002
; Ma et al., 2006
; Wang et al., 2006
).
| RESULTS |
|---|
|
|
|---|
Tryptophan Substitution of S0 Residues at Several Positions Can Produce Shifts in the Voltage-activation of the Channel
Fig. 2 A shows representative current traces obtained from patches containing wild-type BK channels or the L26W mutant, at 9.9 µM Ca2+i, illustrating that a single mutation in S0 could produce a large shift in the voltage activation of the channel toward more positive potentials.
For L26W, the magnitude of these shifts were greater than +50 mV for [Ca2+]
0.4 µM (
V1/2 = 88 ± 3 mV at 9.9 µM Ca2+). Two other mutations in S0 produced relatively large positive G-V shifts, at F25W (
V1/2 = 71 ± 4 mV at 9.9 µM Ca2+) and S29W (
V1/2 = 79 ± 7 mV at 9.9 µM Ca2+). We also observed shifts in voltage activation toward more negative potentials with some of the S0 mutants, which ranged to –18 ± 3 mV for G35W at 2.4 µM Ca2+. Despite the large shifts in V1/2 that were sometimes observed, all of the mutants in these experiments retained both voltage-dependent and Ca2+-dependent gating. Fig. 2 B illustrates that an increase in [Ca2+] at the cytoplasmic side of the patch from 0.4 to 9.9 µM produced a leftward shift of 94 mV for wild-type channels, and a leftward shift of 65 mV for L26W channels.
|
Fig. 3 A summarizes the mean
V1/2 values compared with wild type for each mutant, over a range of [Ca2+], and illustrates that the mutation effects on
V1/2 were sensitive mainly to the position of the mutation.
For example, one of the largest shifts we observed was with the F25W mutation (see above); this mutation gave a modest increase in side chain volume of only
40
3. In contrast, the G35W mutation gave a large increase in side chain volume (
167
3), but resulted G-V curves that were nearly indistinguishable from wild type at 9.9 µM Ca2+. These results suggest that the magnitude and the direction of these shifts are sensitive primarily to the position of the mutation rather than the size of the volume change.
|
We estimated the slopes of G-V curves for each S0 mutant by fitting with the Boltzmann equation, and quantified mutation effects as the fractional change in G-V slope compared with wild type. Mean slopes for wild-type channels from these experiments at 0.4, 2.4, and 9.9 µM Ca2+ were 15.7 ± 1.6, 12.6 ± 0.9, and 11.3 ± 0.8 mV/e-fold change. The results (Fig. 3 B) illustrate that tryptophan mutations in S0 often resulted in slight changes in the slopes of G-V curves. Fractional changes in G-V slope ranged from a 24 ± 6% increase for F25W at 0.4 µM Ca2+ (corresponding to an actual mean slope of 12.0 ± 1.0 mV/e-fold change) to slope decreases of 30% for S28W and L37W, both at 9.9 µM Ca2+.
We found that the slope changes observed in some of the S0 mutants were not statistically significant compared with slopes of wild-type channels at the same Ca2+. This observation in itself may not rule out the possibility that mutations that alter the G-V slope by
30% do not alter the effective gating charge of the channel, or other aspects of channel gating. Possible relations between G-V slope and perturbations in the gating mechanism will be addressed below.
Based on these results, the major effects among these S0 mutations appear to be the shifts in mutant channel V1/2, observed with F25W, L26W, and S29W. Mechanistically, these shifts could result from a change in the equilibria of Ca2+ sensor or voltage sensor movement, a change in the coupling between either or both sensors and pore gating, or a direct effect on the channel's basal open–closed equilibrium (Rothberg and Magleby, 1999
; Cui and Aldrich, 2000
; Rothberg and Magleby, 2000
; Horrigan and Aldrich, 2002
).
S0 Mutations Appear To Act through Mechanisms Other than a Change in Ca2+ Sensitivity
The shifts in V1/2 observed with several S0 mutants are consistent with a role for the S0 segment in channel gating. To begin to understand the role of S0 in gating, we focused on three mutations that yielded a large effect on channel gating, F25W, L26W, and S29W.
A simple test of whether the mutations are having an obligatory effect on the channel's Ca2+ sensing mechanism is to measure channel activity in nominally 0 Ca2+. If the shift observed at higher Ca2+ is due to a change in Ca2+ sensing, then the V1/2 of these channels at 0 Ca2+ should be the same as wild type (i.e., the shift should disappear). We found for F25W, L26W, and S29W that the shifts persisted even at nominally 0 Ca2+ (Fig. 4 A), which argues against these mutation effects acting solely through the Ca2+ activation mechanism (WT V1/2 = 194 ± 5.1 mV; F25W V1/2 = 220 ± 3.1 mV; L26W V1/2 = 232 ± 3.4 mV; S29W V1/2 = 211 ± 3.5 mV). We further explore potential mechanistic consequences of altered Ca2+ sensitivity in a later section.
|
50–300) channels in each patch, measured as NPo. We estimated the number of channels in each patch (N) by first measuring the current amplitude at voltages where all channels were opened, and then dividing this number by the single-channel conductance (Wang et al., 2006We further used the mean values estimated for zL for wild-type and mutant channels at 2.4 µM to estimate values for L0 at this Ca2+ obtained from individual patches, by fitting the Po data obtained at voltages between –85 and –125 mV with an exponential voltage dependence extrapolated to 0 mV.
Our results indicate that L0 estimated at nominally 0 Ca2+ and L0 (estimated at 2.4 and 13 µM Ca2+) could be altered two- to nearly fivefold among the F25W, L26W, or S29W mutant channels (Fig. 4). However, the mean value of L0 estimated for wild-type channels at 2.4 µM Ca2+ (2.4 x 10–5 ± 1.2 x 10–5) was not significantly different from the means estimated for the three mutant channels (8.9 x 10–5 ± 7.9 x 10–5, F25W; 1.2 x 10–5 ± 2.9 x 10–6, L26W; 2.2 x 10–5 ± 4.8 x 10–6, S29W). In addition, it will be shown below that, in the context of the dual allosteric model for BK channel gating (Horrigan and Aldrich, 2002
), changes in L0 of this magnitude are, on their own, insufficient to account for the G-V shifts observed with these mutants over the range of [Ca2+] we examined. These results suggest that a direct change in gate equilibrium is not likely to be solely responsible for the G-V shifts observed with these three mutations.
A rightward shift in V1/2 could also be related to either a decrease in the allosteric voltage-coupling factor (D in the dual allosteric model) (Horrigan and Aldrich, 2002
), a decrease in the effective gating valence of the voltage sensor (zJ), or a decrease in the voltage sensor equilibrium constant (J0). A decrease in either D or zJ that would give rise to a large (>25 mV) G-V shift might be apparent in our G-V data as a decrease in the slopes of G-V curves compared with wild type; for a decrease in D, this slope change would be especially apparent at nominally 0 Ca2+ (discussed below). We have not observed a statistically significant decrease in G-V slope for any of our S0 mutants, including F25W, L26W, or S29W (Fig. 3 B). For these three mutants, we estimated Boltzmann parameters additionally at nominally 0 Ca2+ (Fig. 4 A) and 13 µM Ca2+ (Fig. 4 C), and again they were not distinguished from wild type.
To explore the effects of changing these parameters both qualitatively and quantitatively to gain insight toward possible mechanisms underlying the mutation effects in F25W, L26W, and S29W, we used the dual allosteric model for BK channel gating. We first arrived at two similar sets of parameters that described wild-type gating behavior for data obtained at nominally 0, 2.4, and 13 µM Ca2+, including low open probability data obtained using single-channel recording (Fig. 4). These wild-type fit parameters were constrained by the G-V data presented in Fig. 5 A, and by parameter values obtained previously for BK channels (Horrigan and Aldrich, 2002
; Ma et al., 2006
; Wang et al., 2006
). Then, using these parameters (presented in Table I, Fit A and Fit B), we changed each parameter individually and compared the predicted G-V data with the experimental data obtained for F25W, L26W, and S29W.
Our rationale was that if we could describe the major features of the mutation effects by changing a single parameter in the model, then this would indicate the simplest mechanism that may underlie the mutation effect. It is import to note, however, that this would not necessarily be the only mechanism underlying the mutation effect.
|
|
18-fold decrease), zJ = 0.65 e0 (16% increase), and L0 = 7.5 x 10–7 (3.8-fold increase); the dashed line, however, shows the prediction of Fit B parameters with J0 = 0.00825 (an 8.5-fold decrease), and no changes in the other parameters from wild type. In Fig. 5 D, the solid line shows the prediction of Fit A parameters for S29W with J0 = 0.0074 (an
6.8-fold decrease), and zJ = 0.65 e0 and L0 = 7.5 x 10–7 as with L26W. The dashed line shows the prediction of Fit B parameters with the change of only J0 to 0.022 (a 2.2-fold decrease). Quantitatively, the G-V shifts observed for F25W, L26W, and S29W at nominally 0 Ca2+ are smaller than the shifts observed in the presence of Ca2+ (for F25W, 26 mV at 0 Ca2+, 58 mV at 2.4 µM Ca2+; for L26W, 38 mV at 0 Ca2+, 79 mV at 2.4 µM Ca2+; for S29W, 17 mV at 0 Ca2+, 63 mV at 2.4 µM Ca2+). Would this, in principle, demand an effect on Ca2+ sensing (in the context of the model)? Interestingly, the smaller shifts observed with 0 Ca2+ in these mutants are captured remarkably well, as illustrated in Fig. 5. Although these results do not eliminate the possibility that these mutations may alter Ca2+ sensing or Ca2+–sensor/gate coupling, they do demonstrate that the size of a G-V shift may increase with increasing Ca2+ even with a change only in the voltage sensor equilibrium constant (in the context of the dual allosteric model). The effects of changing model parameters that determine Ca2+ sensing is explored in more detail below.
Exploration of Alternative Mechanisms To Describe the G-V Shifts
In the context of the dual allosteric model, it is possible to generate shifts in G-V curves by changing other parameters that define voltage-dependent gating besides J0, such as L0, zJ, and D, as well as the parameters K and C, which in part define the Ca2+ dependence of gating. We have explored other possible mechanisms in some detail, using different parameter sets in two examples.
Using L26W, which gave the largest G-V shifts, we explored whether changing any one individual parameter could approximate these mutation effects. In Fig. 6 A, we fitted the G-V curves using the wild-type Fit B parameters, and altering only L0.
This revealed that decreasing L0 by
100-fold (necessary to generate a +80-mV G-V shift compared with WT at 2.4 µM Ca2+) could not account either for the observed Po at negative voltages or the G-V relation observed at nominally 0 Ca2+ (Fig. 6 A). Similarly, a decrease in D overpredicted the G-V shift observed at nominally 0 Ca2+, while a decrease in zJ both overpredicted the G-V shift for 0 Ca2+ and underpredicted the G-V shift observed for 13 µM Ca2+ for this mutant (Fig. 6 B). Finally, changing only the Ca2+ sensor parameters C or KD by an amount sufficient to produce the large G-V shifts observed with L26W predicted a striking decrease in Ca2+ dependence (Fig. 6 C).
|
In Fig. 7 A, we fitted the G-V curve for F25W at nominally 0 Ca2+ using the wild-type Fit A parameters, by individually altering either D, zJ, or L0 (changing C or K would not shift the G-V curve at nominally 0 Ca2+).
We then used these sets of parameters to predict the G-V curves at 2.4 and 13 µM Ca2+. This approach revealed that changes in any of these individual parameters that would produce a similar rightward G-V shift at 0 Ca2+ compared with the WT Fit A prediction (to
220 mV, and
+30 mV shift) would not account for the sizes of G-V shifts observed at higher Ca2+. For example, decreases in either D or L0 underpredicted V1/2 values at 2.4 µM (V1/2 = 115 mV by decreasing J0 vs. 83 mV by decreasing D or L0) and at 13 µM Ca2+ (V1/2 = 59 mV by decreasing J0 vs. 26 mV by decreasing D or L0). Similarly, a decrease in zJ underpredicted the V1/2 values at 2.4 µM (V1/2 = 90 mV by decreasing zJ) and at 13 µM Ca2+ (V1/2 = 19 mV by decreasing zJ). In a second approach (Fig. 7, B–D), we fitted the G-V curves for F25W as we did for L26W in Fig. 6, except for F25W we used the wild-type Fit A parameters (and individually changed L0, zJ, D, KD, or C). These results illustrate that, consistent with the results in Fig. 6, changing L0, zJ, or D typically overpredicted the G-V shift observed at nominally 0 Ca2+, and in the case of zJ, also underpredicted the G-V shift observed at 13 µM Ca2+.
|
37% decrease compared with changing J0) and the solid cyan curve (D = 20.7) has a slope of 19 mV/e-fold change (
25% decrease compared with changing J0). Also, the sizes of the G-V shifts resulting from decreases in either zJ or D would be Ca2+ dependent, with larger rightward shifts occurring with lower [Ca2+] (also illustrated in Fig. 7 C). Large reductions in G-V slope were not observed with F25W, L26W, or S29W, and although there was some Ca2+ dependence observed with the G-V shifts for these mutants, the larger shifts occurred with higher [Ca2+]. These observations suggest that the large rightward G-V shifts observed with these mutations are not likely to arise from decreases in zJ or D alone. Interestingly, either raising KD to 68 µM or lowering C to 3.4 predicted mutant V1/2 values that were close to those predicted by changing J0 (Fig. 7 D). However, changing these parameters alone predicted V1/2 values at 0 Ca2+ that were identical to those of WT channels, and also underpredicted the observed Po at voltages less than –100 mV. As with any of the model parameters, we cannot entirely rule out a more complicated model in which a change in Ca2+ sensitivity contributes to the mutation effects along with changes in voltage sensitivity and/or basal open probability (i.e., L0).
Together, these results suggest that changing individual parameters other than J0 cannot generate all of the observed G-V shifts and low Po data as well as changing the parameter J0 alone. Thus, the simplest model to account for the gating effects of F25W, L26W, and S29W can do so primarily by decreasing J0.
| DISCUSSION |
|---|
|
|
|---|
subunit and its auxiliary ß1 subunit is critical to the channel's physiological role in the control of electrical excitability in smooth muscle (Brenner et al., 2000
Function of the S0 Segment
In the present study, we have found that single amino acid changes in the BK channel S0 segment to tryptophan at key positions can produce large changes in channel gating. We went on to analyze the mutation effects systematically, in the context of our current thinking on BK channel gating, to identify a possible functional role of S0. Our results have suggested that the large G-V shifts observed with F25W, L26W, and S29W mutations are consistent with decreases in the voltage sensor equilibrium constant (J0), such that larger depolarizations are required to fully activate the voltage sensor in these mutants compared with wild-type channels.
Although a decrease in J0 could account for the major features of these mutations, there are some features of the mutant data that are not captured perfectly by modifying this single parameter (Fig. 5). One must therefore consider the possibility that additional gating parameters (such as L0 and zJ) are affected slightly, as illustrated in the ability of both Fit A and Fit B to provide reasonably good descriptions of both wild-type and mutant data. In previous studies, BK model parameters have been additionally constrained by including data obtained from gating current experiments, which are technically quite challenging for these channels (Horrigan and Aldrich, 2002
; Bao and Cox, 2005
). Such experiments may further resolve possible roles of S0 in modulating voltage sensor movement.
If these mutations reveal a functional role for the S0 segment in modulating the equilibrium between the resting and active states of the BK voltage sensor, then this may also be consistent with a possible role for S0 in transducing the actions of the ß1 subunit, which include modulation of voltage sensor equilibrium (Bao and Cox, 2005
; Wang and Brenner, 2006
).
A Possible Structure–Function Relation
If we assume that the transmembrane S0 segment is helical, as with the other transmembrane domains of Kv channel voltage sensors (Jiang et al., 2003
; Lee et al., 2005
; Long et al., 2005
), then the positions of Phe-25, Leu-26, and Ser-29 would be localized to a contiguous region on the surface of S0. Fig. 8 A illustrates that this surface would be located on one face of the helix, and the disruptive effects of tryptophan substitution at these positions suggests that these residues may approach or contact a portion of the voltage sensor.
Phe-25, Leu-26, and Ser-29 are among the S0 residues that are identical among all known slo1 orthologues (Fig. 1; Fig. 8 B), which is consistent with the idea that they are important in BK channel structure and function. In addition, chimeras containing the dslo1 S0 region on an hslo1 background appear to retain the gating properties of wild-type hslo1 (Wallner et al., 1996
), consistent with the idea that the identities of these highly conserved residues contribute to determination of the channel's functional properties.
|
| ACKNOWLEDGMENTS |
|---|
This work was supported by grants from the American Heart Association, Texas Affiliate (0265124Y), and the National Institutes of Health (GM68523) to B.S. Rothberg.
Olaf S. Andersen served as editor.
Submitted: 12 September 2006
Accepted: 25 January 2007
| REFERENCES |
|---|
|
|
|---|
-Helical structural elements within the voltage-sensing domains of a K+ channel. J. Gen. Physiol. 115:33–50.[Medline]This article has been cited by other articles:
![]() |
G. Liu, S. I. Zakharov, L. Yang, S.-X. Deng, D. W. Landry, A. Karlin, and S. O. Marx Position and Role of the BK Channel {alpha} Subunit S0 Helix Inferred from Disulfide Crosslinking J. Gen. Physiol., May 26, 2008; 131(6): 537 - 548. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. M. Papazian S0, Where Is It? J. Gen. Physiol., May 26, 2008; 131(6): 531 - 536. [Full Text] [PDF] |
||||
![]() |
I. Schroeder and U.-P. Hansen Saturation and Microsecond Gating of Current Indicate Depletion-induced Instability of the MaxiK Selectivity Filter J. Gen. Physiol., July 1, 2007; 130(1): 83 - 97. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||