Published online Feb 27 2006. doi:10.1085/jgp.200509421
The Rockefeller University Press, 0022-1295 $8.00
JGP, Volume 127, Number 3, 309-328
Role of Charged Residues in the S1S4 Voltage Sensor of BK Channels
Zhongming Ma,
Xing Jian Lou, and
Frank T. Horrigan
Department of Physiology, University of Pennsylvania School of Medicine, Philadelphia, PA 19104
Correspondence to Frank T. Horrigan: Horrigan{at}mail.med.upenn.edu
The activation of large conductance Ca2+-activated (BK) potassium channels is weakly voltage dependent compared to Shaker and other voltage-gated K+ (KV) channels. Yet BK and KV channels share many conserved charged residues in transmembrane segments S1S4. We mutated these residues individually in mSlo1 BK channels to determine their role in voltage gating, and characterized the voltage dependence of steady-state activation (Po) and IK kinetics (
(IK)) over an extended voltage range in 050 µM [Ca2+]i. mSlo1 contains several positively charged arginines in S4, but only one (R213) together with residues in S2 (D153, R167) and S3 (D186) are potentially voltage sensing based on the ability of charge-altering mutations to reduce the maximal voltage dependence of PO. The voltage dependence of PO and
(IK) at extreme negative potentials was also reduced, implying that the closedopen conformational change and voltage sensor activation share a common source of gating charge. Although the position of charged residues in the BK and KV channel sequence appears conserved, the distribution of voltage-sensing residues is not. Thus the weak voltage dependence of BK channel activation does not merely reflect a lack of charge but likely differences with respect to KV channels in the position and movement of charged residues within the electric field. Although mutation of most sites in S1S4 did not reduce gating charge, they often altered the equilibrium constant for voltage sensor activation. In particular, neutralization of R207 or R210 in S4 stabilizes the activated state by 37 kcal mol1, indicating a strong contribution of nonvoltage-sensing residues to channel function, consistent with their participation in state-dependent salt bridge interactions. Mutations in S4 and S3 (R210E, D186A, and E180A) also unexpectedly weakened the allosteric coupling of voltage sensor activation to channel opening. The implications of our findings for BK channel voltage gating and general mechanisms of voltage sensor activation are discussed.
Abbreviations used in this paper: BK, large conductance Ca2+-activated; MES, methanesulfonic acid; WT, wild-type.
X.J. Lou's present address is Johnson & Johnson Pharmaceutical Research and Development, Spring House, PA 19477.
 |
INTRODUCTION
|
|---|
To sense voltage, ion channels must contain charged residues whose position in the membrane electric field change upon voltage sensor activation. Identifying these residues and their contributions to voltage sensitivity are critical to understanding the gating of any voltage-dependent channel. Voltage-dependent K+ channels share a common "core" architecture, including six transmembrane segments (S1S6) with an S1S4 voltage sensor and S5S6 pore domain in each subunit (Fig. 1 A). Many conserved basic and acidic amino acids are present in S1S4 that could potentially contribute to voltage sensing (Fig. 1 B). However, in Shaker, only four residues in S4 and one in S2 account for the majority of gating charge (Fig. 1 B, double starred positions) (Aggarwal and MacKinnon, 1996
; Seoh et al., 1996
). That a few charged residues make large contributions to voltage sensitivity is supported by evidence in Shaker (Larsson et al., 1996
; Starace and Bezanilla, 2001
), hERG (Zhang et al., 2004
), and Na+ channels (Yang et al., 1996
) that the accessibility of S4 positions to internal and external solutions can switch sides upon voltage sensor activation, implying that some charges traverse the entire electric field. Although the majority of charged residues in S1S4 appear to be nonvoltage sensing, their conservation suggests they also play important roles in channel structure or function. Salt bridge interactions between oppositely charged amino acid side chains in different transmembrane segments, including nonvoltage-sensing residues, are observed in the crystal structures of KvAP (Jiang et al., 2003
) and Kv1.2 (Long et al., 2005a
), and impact Shaker biogenesis (Papazian et al., 1995
; Tiwari-Woodruff et al., 1997
). Nonvoltage-sensing residues also influence channel gating, suggesting they participate in interactions that depend on the activation state of the voltage sensor (Tiwari-Woodruff et al., 2000
).
Large conductance Ca2+-activated (BK) channels (Slo1) exhibit a pattern of charged residues in S1S4, similar to that of Shaker and other strongly voltage-dependent K (Kv) channels (Fig. 1 B). However, an equivalent of only 2.6 e per channel are estimated to traverse the electric field upon activation of mSlo1 (Horrigan and Aldrich, 1999
), a fivefold smaller gating charge than Shaker (13 e) (Aggarwal and MacKinnon, 1996
; Seoh et al., 1996
). The molecular basis for this marked difference in voltage sensitivity is unclear. Slo1 channels are Ca2+ sensitive; but Ca2+ acts independent of voltage and has no effect on gating charge (Horrigan and Aldrich, 2002
). Slo1 lacks two of the five most important voltage-sensing residues in Shaker (Fig. 1 B), including one in S4; but this difference seems insufficient to account for a fivefold reduction in gating charge. Slo1 retains a series of three regularly spaced arginines in S4 that match key voltage-sensing residues in Shaker even when alternative alignments are considered (Diaz et al., 1998
). Moreover, the pattern of arginines in the NH2-terminal half of S4 in Slo1 appears identical to that of Kv2.1 (Fig. 1 B), which exhibits a gating charge similar to that of Shaker (Islas and Sigworth, 1999
).
Several mechanisms could account for the relatively weak voltage dependence of BK channel activation. First, voltage sensor movement in Slo1 may be limited compared with Shaker, such that the number of voltage-sensing positions is reduced and/or each residue makes a smaller contribution to gating charge. A reduced contribution of individual residues seems likely, since even one S4 arginine traversing the entire electric field, as in Shaker, should exceed the total gating charge of mSlo1. Alternatively, S4 residues could make large contributions to gating charge if additional sites make contributions of opposite polarity such that the net gating charge is small. Finally, fundamental differences in the structure and/or movement of voltage sensors in BK and KV channels may exist, such that voltage-sensing residues in Slo1 correspond to nonvoltage-sensing positions in Shaker, and the pattern of charge residues is not predictive of differences in voltage sensitivity.
Here we examine the effects on voltage-dependent activation of mutating each charged residue in the S1S4 segments of mSlo1 to understand their role in channel function, to identify voltage-sensing residues, and to determine why BK channels are weakly voltage dependent. Our results indicate that while many charged residues influenced BK channel gating, only four are potentially voltage sensing, and only one is in S4. Charge neutralization at any of these sites reduces voltage sensitivity but the decreases in gating charge are small (
0.3 e per subunit). The limited contribution of S4 suggests that voltage sensor movement in BK channels is reduced relative to KV channels and indicates that charged residues outside of S4 make significant contributions to gating charge. That potential voltage-sensing residues in S2 are not in the same location for mSlo1 and Shaker suggests that additional differences exist in the movement or conformation of voltage sensors in these channels. Thus several factors may contribute to the weak voltage dependence of BK channel activation. In addition, our results show that nonvoltage-sensing residues make important contributions to voltage gating, in some cases influencing the energetic coupling of voltage sensor activation to channel opening, in others profoundly influencing the equilibrium constant for voltage sensor activation.
 |
MATERIALS AND METHODS
|
|---|
Channel Expression
Experiments were performed with the mbr5 clone of the mouse homologue of the Slo1 gene (mSlo1) (Butler et al., 1993
) expressed in Xenopus oocytes. The clone was modified to facilitate mutagenesis and was propagated and cRNA transcribed as described previously (Cox et al., 1997
). Oocytes were injected with
0.550 ng of cRNA, incubated at 18°C, and studied 37 d after injection. Site-directed mutagenesis was performed with the QuikChange XL Site-Directed Mutagenesis Kit (Stratagene) and confirmed by sequencing. Some mutants (E180A, D186A, R210N, and E219Q) were provided by J. Cui (Washington University, St Louis, MO).
Electrophysiology and Data Analysis
Currents were recorded using the patch clamp technique in the inside-out or outside-out configuration (Hamill et al., 1981
). For IK measurements, the internal solution contained (in mM) 110 K-methanesulfonic acid (K-MES), 20 HEPES. 40 µM (+)-18-crown-6-tetracarboxylic acid (18C6TA) was added to chelate contaminant Ba2+ (Diaz et al., 1996
; Neyton, 1996
). In addition, "0 Ca2+" solution contained 5 mM EGTA, reducing free Ca2+ to an estimated 0.8 nM in the presence of
10 µM contaminant Ca2+. Solutions containing Ca2+ were buffered with 5 mM HEDTA, and free [Ca2+] was measured with a Ca2+ electrode (Orion Research Inc.). Nominal [Ca2+] reported as 1, 10, and 50 µM corresponded to measured concentrations of 0.87, 8.7, and 49.9 µM, respectively. Ca2+ was added as CaCl2, and [Cl] was adjusted to 10 mM with HCl. The external solution contained 110 K-MES, 2 MgCl2, 6 HCl, 20 HEPES. For gating current measurements, the internal solution contained 135 NMDG-MES, 6 NMDG-Cl, 20 HEPES, 2 EGTA, 40 µM 18C6TA. The external solution contained 125 TEA-MES, 2 TEA-Cl, 2 MgCl2, 20 HEPES. The pH of all solutions was adjusted to 7.2 with MES. Experiments were performed at room temperature (2022°C).
Data were acquired with an Axopatch 200B amplifier (Axon Instruments, Inc.) in patch mode with the Axopatch's filter set at 100 kHz. Currents were subsequently filtered by an 8-pole Bessel filter (Frequency Device, Inc.) at 20 kHz and sampled at 100 kHz with an 18-bit A/D converter (Instrutech ITC-18). A P/4 protocol was used for leak subtraction (Armstrong and Bezanilla, 1974
) with a holding potential of 80 mV. For some mutants, with increased PO, current traces in high [Ca2+] were leak subtracted using a leak trace recorded in 0 Ca2+ at V
140 mV before and after changing [Ca2+]i. Electrodes were made from thick-walled 1010 glass (World Precision Instruments, Inc.) and their tips coated with wax (KERR Sticky Wax). The electrode's resistance in the bath solution (1.02.5 M
) was used as an estimate of series of resistance (RS) for correcting the voltage at which macroscopic IK was recorded. Series resistance error was <15 mV for all data presented. A Macintosh-based computer system was used in combination with Pulse Control acquisition software (Herrington and Bookman, 1995
) and Igor Pro for graphing and data analysis (WaveMetrics, Inc.). A Levenberg-Marquardt algorithm was used to perform nonlinear least-squared fits.
Open probability (PO) was estimated over a wide voltage range by recording macroscopic IK when PO was high (
0.0050.05), and single channel currents in the same patch when PO was low (
0.010.1). Macroscopic conductance (GK) was determined from tail currents at 80 mV, following 30-ms voltage pulses, and was normalized by GKmax measured in 50 µM Ca2+ to estimate PO. At more negative voltages, NPO was determined from steady-state recordings of 160 s duration that were digitally filtered at 5 kHz. NPO was determined from all-points amplitude histograms by measuring the fraction of time spent (PK) at each open level (K) using a half-amplitude criteria and summing their contributions.
. Po was then determined by estimating N from GKmax in 50 µM Ca2+ (N = GKmax/gK, where gK is the single channel conductance at 80 mV).
Patch to patch variation in the half-activation voltage (Vh) of the GK-V relationship is observed for Slo1 channels (Stefani et al., 1997
; Horrigan et al., 1999
) and causes broadening in the averaged voltage-dependent relationships relative to individual experiments. To compensate for such variation, Vh was determined for each patch and Po-V relations were shifted along the voltage axis by
Vh = (
Vh
Vh), where
Vh
is the mean for all experiments at the same [Ca2+], before averaging (Horrigan and Aldrich, 1999
; Horrigan et al., 1999
).
Mean activation charge displacement (qa = kT d(ln(PO)/dV)) was measured from the slope of the ln(PO)-V relation by linear regression over 60-mV intervals (approximately four data points) and plotted against mean voltage. This procedure minimized noise while introducing errors of at most 5% in qa, based on simulations using the HCA model (Horrigan et al., 1999
). Fits of the qa-V relation by gating schemes were similarly determined by simulating the data at 20-mV intervals and determining qa by linear regression.
Gating currents were measured using admittance analysis as previously described (Horrigan and Aldrich, 1999
). In brief, inside-out patches were excised into nominally K+-free solutions containing isotonic TEA in the extracellular solution to block residual ionic currents. Currents were recorded in response to 0.51-s voltage ramps upon which a sinusoidal voltage command (1736 hz, 60 mV peak to peak) was superimposed. Admittance was calculated for each cycle of the sinwave, and capacitance was determined after correcting for phase shifts due to instrumentation. The voltage-dependent component of the capacitance signal due to gating current was integrated over the voltage range to determine the Q-V relation, which was sampled and plotted at 20-mV intervals. The resulting Q-V is a pseudo steady-state measurement that approximates the closed channel charge distribution (QC) at voltages where PO is small (Horrigan and Aldrich, 1999
).
The steady-state data in 0 Ca2+ were fit with the HCA model using a semi-automated procedure to estimate the best values and standard deviation for each parameter based on least-squares criteria. After determining the best fit to the log(Po)-V relation, with all parameters allowed to vary, the parameters were refined by iteratively fitting qa-V and log(PO)-V. First the voltage sensor charge (zJ) was varied to fit qa at voltages where PO < 0.1. Under these conditions, qa is highly dependent on zJ. Log(PO)-V was then fit with all parameters except zJ allowed to vary. The procedure was repeated until zJ approached a constant value. In this way, reasonable fits were obtained to both log(PO) and qa, and the values of charge (zJ) and coupling (D) were better constrained than by fitting either relation alone. In some cases where gating was shifted to very negative or positive voltages (R207Q,E; R210C; R213C,E; and D153C,K), the coupling factor D was poorly constrained and was held to the wild-type (WT) value until zJ was determined; and then D was allowed to vary with zJ held constant. Data in multiple [Ca2+] were fit by eye using the HA model (see RESULTS).
 |
RESULTS
|
|---|
Mutation of Charged Residues in S1S4 Alter mSlo1 Gating
To study the role of charged residues in the BK channel voltage sensor, we examined the effects on the voltage dependence of steady-state activation (PO) and macroscopic kinetics (
(IK)) of mutating individual sites in the S1, S2, S3, and S4 segments of mSlo1 indicated in Fig. 1 B. Charged residues were replaced by neutral (A, C, N, Q) or charged (E, K) amino acids to determine the effects of reducing or reversing charge. In addition, a neutral residue in S2 (Y163) was replaced by Glu (E), to match a voltage-sensing position in Shaker (E293). IK was recorded over a wide voltage range in the virtual absence of internal Ca2+ (0 Ca2+, see MATERIALS AND METHODS) to isolate effects of mutation on voltage-dependent gating from those on Ca2+-dependent mechanisms or the closed to open conformation change (Horrigan and Aldrich, 2002
). Measurements were also performed in various [Ca2+]i, including a nearly saturating concentration (50 µM) used to determine the maximal conductance (GKmax) and estimate the number of channels (N) in each patch to determine Po (see MATERIALS AND METHODS).
Charge-altering mutations at all sites tested produced functional channels that express macroscopic currents in excised patches. In response to voltage pulses, currents recorded in 0 Ca2+ or 50 µM Ca2+ activate and deactivate with an approximately exponential time course, qualitatively similar to the WT channel (Fig. 2, AC). However most mutations in S2, S3, and S4 altered channel gating significantly, as shown in Fig. 2 D by comparing normalized GK-V relations for WT and mutant channels in 0 Ca2+. Although R207Q, R207E, or R210C in S4 increased PO relative to the WT and shifted G-V curves to more negative voltages, most mutations inhibited activation. Indeed some mutants appeared almost nonfunctional in 0 Ca2+, but we were able to characterize their gating properties owing to the high single channel conductance and Ca2+ sensitivity of BK channels. For instance, a voltage pulse to +180 mV in 0 Ca2+ evoked only single channel currents from a patch expressing R213C channels (Fig. 2 C, left), whereas macroscopic currents are evoked from the same patch in 50 µM Ca2+ (Fig. 2 C, right). This shows that hundreds of channels are present and functional, and that PO in 0 Ca2+ is low but measurable based on the resolution of discrete opening events.

View larger version (50K):
[in this window]
[in a new window]
|
Figure 2. Mutation of charged residues in voltage sensor domain alter mSlo1 gating. Families of IK evoked from (A) WT (B) R207Q, and (C) R213C channels in response to pulses to different voltages at 20-mV intervals from a holding potential of 80 mV. Currents in 0 Ca2+ and 50 µM Ca2+ for each channel were recorded from the same patch. Dashed lines indicate the zero current level. (D) GK-V relations (mean ± SEM) for WT and mutant channels are normalized by GKmax determined in 50 µM Ca2+ and are fit by Boltzmann functions (lines). WT (), R113Q ( ), D133Q ( ), D153K ( ), D153C ( ), Y163E ( ), R167E ( ), R167A ( ), E180A ( ), D186A ( ), R201Q ( ), R207E ( ), R207Q ( ), R210C ( ), R210E ( ), R210N ( ), R213C ( ), R213E ( ), E219Q ( ).
| |
The following two sections and Fig. 3 describe the analytical framework used to identify voltage-sensing residues and determine which aspects of channel gating are altered by various mutations. The results of this analysis are summarized in Fig. 4 and Table I.

View larger version (39K):
[in this window]
[in a new window]
|
Figure 3. Data analysis and modeling. (A) Po-V relation for WT mSlo1 channels in 0 Ca2+ (mean ± SEM) is plotted on a semilog scale and fit by the HCA model (solid line, see be-low). The limiting slope at negative voltages is indicated by a dashed line representing the partial charge associated with channel opening (zL = 0.3 e). (B) qa-V relation determined from the logarithmic slope of mean PO (see MATERIALS AND METHODS). (C) HCA gating mechanism asserts that equilibration of the channel gate between closed and open (CO) is influenced allosterically by four independent and identical voltage sensors that each can be in a resting or activated state (RA). The equilibrium constants for channel opening L and voltage sensor activation J are voltage dependent with partial charges zL and zJ, respectively ( ). The activation of each voltage sensor increases the CO equilibrium constant by an allosteric coupling factor D. (D) Gating scheme shows the closed and open states (Ci, Oi, where i = number of activated voltage sensors) and equilibrium constants specified by the HCA mechanism. Horizontal transitions represent voltage sensor activation, and vertical transitions represent channel opening. Solid lines in A and B represent best fits to the HCA model (zJ = 0.59 e, zL = 0.3 e, J0 = 0.017, Lo = 1.4 x 106, D = 23.1) determined by a semi-automated method (see MATERIALS AND METHODS). Dashed lines in B indicate the charge distribution for open channels (QO, Eq. 5) and the difference between charge distributions for open and closed channels Qo Qc, Eq. 6), defined by the fit parameters. (EH) Predicted effects of changing various parameters in the HCA model. (E) Changes in peak qa, expressed as a fraction of the control value (fqaMAX), produced by changing gating charge (zJ, zL) are directly proportional to the fractional change in total gating charge (fzT, where zT = zL + 4zJ). Dashed line indicates unity slope (b = 1). Solid lines indicate effects of proportional changes in both zL and zJ (b = 0.99) or changing zJ alone (b = 0.94). (F) Effects on fqaMAX of fractional changes in L0, J0, or D. (G) Effects on PO-V relation of changing zJ(0.35 e), D(5.1), J0(1.7), or L0(1.4 x 104) and leaving other parameters the same as the WT (solid line). (H) The decreases in zJ or D both decrease qaMAX relative to the WT (solid line) but have distinct effects on the qa-V relations. In particular, the steepness of the foot of the curve is reduced together with zJ but is not altered by changing D.
| |

View larger version (46K):
[in this window]
[in a new window]
|
Figure 4. Effects of voltage sensor mutations on gating charge and allosteric coupling. (A) qaMAX for each voltage sensor mutant was determined from the mean PO-V relation in 0 Ca2+ (see MATERIALS AND METHODS) and normalized to that of the WT (1.89 e). qaMAX for R210C could not be determined for reasons discussed in the text. (B) The total gating charge (zT = zL + 4zJ) and (C) allosteric coupling factor D were estimated from fits of PO-V and qa-V relations to the allosteric models and normalized to those for the WT (zT = 2.62 e, D = 23.9). Shaded bars were determined by semi-automated fitting of the HCA model to 0 Ca2+ data (see MATERIALS AND METHODS). Solid bars represent fits to data in all [Ca2+] using the HA model (Table I) with parameters adjusted manually. (D and E) Confidence limits for manual fits were estimated for zT or D by adjusting zJ (panel D) or D (panel E) respectively from the best fit values (solid line) to exceed the error range of PO as indicated by dashed lines.
| |
Identification of Voltage-sensing Residues
Voltage-sensing residues contribute to the total gating charge (zT) moved in each channel upon voltage sensor activation and channel opening. Therefore, to identify these residues we must determine the impact of charge-altering mutations on zT. The GK-V relations in Fig. 2 D are fit by Boltzmann functions, characterized by an apparent charge zAPP that was in many cases reduced. However, effects of mutation on zAPP are not sufficient evidence that a residue is voltage sensing (Seoh et al., 1996
). A more direct indication of gating charge is provided by the mean activation charge displacement qa, defined in terms of the logarithmic slope of the PO-V relationship (Sigg and Bezanilla, 1997
):
 | (1) |
For channels that open only when all voltage sensors are activated, qa increases at negative voltages approaching a maximum or limiting slope (qalim) equal to zT when Po is small (Almers, 1978
; Sigg and Bezanilla, 1997
). The limiting slope has been used to estimate zT for KV channels (Seoh et al., 1996
; Islas and Sigworth, 1999
) and Na+ channels (Hirschberg et al., 1995
). However, for BK channels, the relationship between qa and zT is more complex because channels can open when no voltage sensors are activated (Horrigan and Aldrich, 2002
). The mean Po-V relation for mSlo1 in 0 Ca2+ is plotted on a semilog scale in Fig. 3 A and the corresponding qa-V relation in Fig. 3 B. At extreme negative voltages (<100 mV) Po becomes very small (
107106) and weakly voltage dependent, exhibiting a limiting slope (qaLim = 0.3 e; Fig. 3 A, dashed line) that represents the small charge associated with channel opening. At more positive voltages, qa increases to a peak (qaMAX = 1.9 ± 0.1 e) that reflects voltage sensor activation (Horrigan and Aldrich, 2002
).
The relationship between qa and zT can be understood based on gating current measurements and some assumptions about voltage gating embodied in a model that accounts for BK channel gating in 0 Ca2+ (Fig. 3, C and D) (Horrigan et al., 1999
). The HCA model postulates that the transition between the closed (C) and open (O) conformation is influenced allosterically by the state (resting [R] or activated [A]) of four independent and identical voltage sensors (Fig. 3 C). The energetic effect of voltage sensor activation on the CO equilibrium is represented by a factor D and is fully described in terms of a two-tiered gating scheme (Fig. 3 D) representing the closed and open channel with different numbers (04) of activated voltage sensors. The equilibrium constants for the CO and RA transitions are both voltage dependent with partial charges zL and zJ, respectively. Therefore the total gating charge involved in activating all four voltage sensors and opening the channel is zT = (4zJ + zL). zJ can be determined from the QC-V relation, which is measured from gating current when channels are closed and is well described by a Boltzmann function with a mean charge of zJ = 0.58 e (Horrigan and Aldrich, 2002
). zL can be determined from the limiting slope qaLim = 0.3 e. Therefore, the estimated value of zT = 2.62 e exceeds qa at all voltages (Fig. 3 B). However the HCA model reproduces the Po-V and qa-V relations using these estimates of zJ, and zL (Fig. 3, A and B, solid curves) and therefore reproduces the relationship between zT and qa.
In the present study we assessed the impact of mutations on zT by measuring qaMAX, and by fitting the Po-V relations in various [Ca2+] with an expanded version of the HCA model that includes the effects of Ca2+ (Table I, HA model) (Horrigan and Aldrich, 2002
). In principle, zT could be estimated as for the WT channel, from gating currents and limiting slope. However, the expression of mutant channels was usually insufficient for gating current measurements. On the other hand, moderate levels of channel expression are sufficient to routinely measure the PO-V relation and qaMAX even at very low PO, owing to the high single channel conductance of BK channels.
Although qaMAX underestimates the total gating charge, simulations based on the HCA model imply that qaMAX provides a reasonable indicator of changes in zT. That is, the fractional change in qaMAX for a mutant relative to the WT (
) approximates the fractional change in zT (
) under most conditions. Fig. 3 E plots fqaMAX versus fzT as charge in the model is altered either by varying zJ alone or by varying zJ and zL proportionally (solid lines). In both cases, a linear relationship is observed with a slope near unity, indicating that fqaMAX approximates fzT over a wide range of gating charge. qaMAX is also sensitive to changes in the allosteric coupling factor D but is virtually unaffected by large changes in the equilibrium constants for voltage sensor activation (J0) or channel opening (L0) (Fig. 3 F). Examples of the effects on the Po-V and qa-V relations of altering the parameters zJ, L0, J0, or D individually are illustrated in Fig. 3 (G and H), respectively.
Distinguishing between Changes in Gating Charge and Allosteric Coupling
Because changes in qaMAX may reflect changes in coupling (D) as well as changes in gating charge (zJ, zL), it is important that we distinguish between these two possibilities. In general, we were able to determine which parameters in the HA model were altered by fitting the steady-state data (PO and qa) over a wide range of conditions. Many of the parameters are highly constrained by subsets of the data at extremes of voltage, [Ca2+], and PO (Horrigan and Aldrich, 2002
). When fitting these data we also focused on several features that depend specifically on coupling or charge. As described below, the magnitude of the Po increase between extreme positive and negative voltages and the shape of the qa-V relation provide direct indications of changes in coupling and charge respectively.
Allosteric coupling determines the extent to which the CO equilibrium constant and PO is increased by voltage sensor activation (Horrigan et al., 1999
). Therefore changes in coupling should alter the magnitude of the PO increase between extreme negative and positive voltages where voltage sensors are in the resting or activated state, respectively. For the WT channel in 0 Ca2+, PO increases from a weakly voltage-dependent plateau near 107 at extreme negative voltages to a level close to unity at voltages (
+300 mV) where the Q-V relation is saturating and voltage sensors are fully activated (Horrigan and Aldrich, 1999
). A simulated decrease in gating charge (zJ) that reduces qaMAX by 35% has no effect on the magnitude of the PO increase between extreme negative and positive voltages (Fig. 3 G,
zJ). However, a decrease in coupling (D factor) that also reduces qaMAX by 35% greatly inhibits the voltage-dependent PO increase (Fig. 3 G,
D). In this case, PO fails to approach unity when voltage sensors are fully activated and instead achieves a plateau whose weak voltage dependence, like that at negative voltages, represents the charge associated with channel opening (zL). Thus, decreases in coupling or charge should have effects on the PO-V relation that are readily distinguishable, provided the resting and activated states of the voltage sensor are experimentally accessible, and Po does not saturate before all voltage sensors are activated. Saturation could result from enhanced coupling (increased D factor) or stabilization of the open state (increased L0) caused by mutation. Likewise, the presence of Ca2+ can obscure changes in coupling by stabilizing the open state.
Changes in gating charge have effects at negative voltages on the voltage dependence of qa and PO that are independent of changes in coupling. For any channel with multiple open states the qa-V relation can be described as a function of PO and the charge distributions (Q-V relations) for open and closed channels (Horrigan and Aldrich, 2002
):
 | (2) |
This function can be approximated by the difference of QO and QC when Po is small:
 | (3) |
when PO << 1.
The relationship between gating charge and qa defined by Eq. 3 is model independent and analogous to the limiting-slope analysis that has been applied to Shaker and other voltage-gated channels. In general, the shape and magnitude of QO-V and QC-V relations depend on gating charge, whereas coupling influences only their relative position along the voltage axis (Horrigan and Aldrich, 2002
). Thus, Eq. 3 describes a bell-shaped function whose steepness is determined by charge and whose width is influenced by coupling. These properties can be defined exactly in the case of the HCA model, where both QO and QC are described by Boltzmann functions with identical voltage dependence but different half-activation voltages (VHC and VHO).
where
 | (4) |
where
 | (5) |
Combining Eqs. 35 yields an expression for qa based on the HCA model:
 | (6) |
when PO << 1.
At extreme negative voltages qa can be further approximated by QO (Eq. 5) (Horrigan and Aldrich, 2002
). Fig. 3 B compares the voltage dependence of qa from the HCA model (solid line) to that of Eqs. 5 and 6 (dashed lines). Eq. 6 provides a good approximation to the qa-V relation near its peak and at more negative voltages, and therefore helps describe several properties of qa. First, qaMAX derived from Eq. 6 depends only on gating charge and coupling, as suggested by the simulations in Fig. 3 (E and F):
 | (7) |
Second, the width of the qa-V relation is related to the difference in half-activation voltages (
VH) between QO and QC and will therefore depend on both voltage sensor charge (zJ) and coupling (D):
 | (8) |
Importantly,
VH has a reciprocal dependence on zJ and D that is distinct from that of qaMAX. Therefore, decreases in charge or coupling that reduce qaMAX will increase or decrease
VH, respectively, having opposite effects on the width of the qa-V relation.
Simulations in Fig. 3 H demonstrate that decreases in D or zJ that have identical effects on qaMAX have distinct effects on qa-V shape. In addition to the opposing effects of D and zJ on qa-V width, the steepness of qa-V is altered by gating charge but is independent of coupling. In particular, the foot of the qa-V relation is well approximated by QO (Fig. 3 B) and therefore has a voltage sensitivity that reflects only zJ and a voltage-independent component equal to zL. In practice, measurements at the foot of the qa-V relation for many mutants were noisy owing to the small amplitude of the corresponding PO data. Therefore, to determine zJ, we took the approach of fitting qa with the HCA model over a broader voltage range that includes qaMAX. However, the fit was restricted to voltages where PO < 0.1, satisfying the simplifying assumption in Eq. 6. In this way, fits depend on coupling as well as charge, but free parameters were limited and their effects could be distinguished based on the above analysis.
Effects of Voltage Sensor Mutations on Gating Charge and Coupling
The effects of individual S1S4 mutations on gating charge were assessed by comparing fractional changes in the maximal logarithmic slope of the PO-V relation (qaMAX) (Fig. 4 A) to those in total gating charge (zT) (Fig. 4 B) and coupling (D) (Fig. 4 C) estimated from fits to the allosteric models. Mutations at six positions in S2 (D153, R167), S3 (E180, D186), and S4 (R210, R213) reduced qaMAX by significantly more than 10%. Four of these (D153, R167, D186, and R213) also reduced zT and are potentially voltage sensing. Decreases in qaMAX by R210E and E180A appear entirely accounted for by a reduction in the coupling factor D, whereas D186A reduces both charge and coupling.
The total gating charge (zT = 4zJ + zL) and coupling factor D were estimated by manual and semi-automated fitting protocols with similar results (solid and shaded bars in Fig. 4, B and C). The mean PO-V, log(PO)-V, and qa-V relations in the presence and absence of Ca2+ were fit simultaneously using the HA model, with parameters adjusted manually (Table I) to provide the best fit by eye. The charge associated with the closedopen transition (zL) was additionally constrained by the voltage dependence of
(IK), as described below. The dissociation constant for Ca2+ (KD) and a weak interaction between voltage sensor activation and Ca2+ binding (allosteric factor E) have little impact on the determination of zT or D and were held constant as a simplifying assumption (see Table I legend). Confidence intervals for zT or D were estimated from the smallest change in zJ (Fig. 4 D) or D (Fig. 4 E) required to exceed the error range (SEM) of Po. Because data was obtained over a wide range of PO, V, and [Ca2+], using macroscopic and single channel methods, the parameters were well constrained, but we could not determine a satisfactory method of weighting the different data to allow automated fitting. We were able to fit the 0 Ca2+ data using a semi-automated procedure (see MATERIALS AND METHODS) that yielded values of zT and D (Fig. 4, B and C, shaded bars) consistent with the manual estimates. Based on comparison of the mean and errors bars determined by these two methods a 10% change in zT or qaMAX or a 20% change in D were deemed significant.
Potential Voltage-sensing Residues in S2, S3, and S4
Steady-state data for mutants at the four potential voltage-sensing sites are compared with that of the WT in Figs. 5 and 6. Mean log(PO)-V and qa-V relations for the WT in 050 µM Ca2+ are plotted in Fig. 5 A and Fig. 6 A, respectively, together with fits to the HA model (solid lines). The corresponding data for R213C, D153K, R167E, and D186A in Fig. 5 (BE) and Fig. 6 (BE) are also fit with the HA model (solid lines) and compared with the WT fit (dashed lines). For the purpose of illustration, qa is shown only in 0 Ca2+ and 50 µM Ca2+. However, it is evident from Fig. 5 that the maximal slope of the PO-V relations (i.e., qaMAX) for mutants were similarly reduced at all [Ca2+]. In addition, mutations increased the width and decreased the steepness of qa-V relations, consistent with decreases in gating charge (Fig. 6).

View larger version (41K):
[in this window]
[in a new window]
|
Figure 5. Effect of potential voltage-sensing residue mutations on PO. Mean log(PO)-V relations for (A) WT, (B) R213C, (C) D153K, (D) R167A, and (E) D186A in 050 µM Ca2+. Solid lines are the fits to the HA model (Table I parameters). Dashed lines represent the WT fit. (F) Log(PO)-V relations in 0 Ca2+ for WT ( ), D153C( ), D153K ( ), R213C ( ), R213E ( ), R167A ( ), R167E ( ), D186A(), and R113Q ( ) are superimposed to show that mutations reduce voltage dependence at the foot of the curve, consistent with a change in gating charge. Mutant curves were shifted along both axes as described below to match the WT at 90 mV. According to analysis presented in the text, the qa-V relation can be approximated by QO (Eq. 5) at extreme negative voltages. Therefore ln(PO) can be approximated as the sum of a linear and exponential function: ln(PO) = [ln(L0/(1 + L0))+(zL/kT)V] + [4DJ0exp(zJV/kT)]. Curves were shifted along the voltage axis such that the exponential component of this equation for all curves equaled 0.1 log unit at 90 mV, representing the transition between a linear and nonlinear function (i.e., the foot of the curve). Thus, the voltage shift for mutants was V = V0.1(WT) V0.1(mutant) where V0.1= (kT/zJ)ln[0.1/(4DJ0log(e))], and the parameters zJ, D, and J0 were taken from Table I.
| |
The data in Fig. 5 (BE) and Fig. 6 (BE) show the effects of a single mutation at each potential voltage-sensing site. At three of these sites (R213, D153, and R167), the effects of both charge neutralization and charge reversal were compared and found to be qualitatively similar (Fig. 4, AC; Table I). However the effects on zJ of charge reversal were always greater than neutralization (Fig. 5 F; Table I), consistent with these residues acting as voltage sensors. The log(PO)-V relations in 0 Ca2+ for these mutants and the WT are superimposed in Fig. 5 F. Because the qa-V relation at extreme negative voltages can be approximated by QO (Eq. 5; Fig. 3 B), the foot of the log(PO)-V relation is expected to change shape only if gating charge is altered. In all cases, the mutants in Fig. 5 F show a reduced voltage sensitivity at the most negative voltages, consistent with decreases in gating charge.
The Location of Voltage-sensing Positions in S2
Two of the four potential voltage-sensing residues in mSlo1 are in S2 (D153, R167). KV channels also contain conserved acidic and basic residues at these positions (Fig. 1 B). However, the residue corresponding to D153 in Shaker (E283) is nonvoltage sensing (Seoh et al., 1996
), and the functional role of that corresponding to R167 has not to our knowledge been reported. Shaker does have a potential voltage-sensing residue in S2 (E293) (Seoh et al., 1996
) that corresponds to an uncharged residue in Slo1 (Y163) (Fig. 1 B). To determine if this position can sense voltage in mSlo1 we replaced it with the corresponding residue in Shaker (Y163E), but no significant change in qaMAX or zT was observed (Fig. 4, A and B; Table I). Thus the distribution of voltage-sensing and nonvoltage-sensing positions in the S2 of BK and KV channels appears to be different.
The Voltage Sensor Contributes to Charge Associated with Channel Opening (zL)
The limiting slope of Po at negative voltages (qalim) represents the gating charge that moves when channels open (zL) while voltage sensors are not activated (i.e., C0O0 transition in HCA model; Fig. 3 D). The value of zL for the WT (0.30 e, Table I) was determined as the slope of mean log(Po) from 140 to 280 mV (50 µM Ca2+: zL = 0.295 ± 0.012 [n = 10]; 5 µM Ca2+: zL = 0.303 ± 0.021 [n = 3]) and confirms a previous estimate over a more limited voltage range (Horrigan and Aldrich, 2002
). zL is small and could result from subtle changes in the pore domain associated with channel opening, including redistribution of ions and the electric field within the pore (Jiang et al., 2002
) or movement of charged residues that are present in S6. However, mutation of potential voltage-sensing residues in S2 (D153, R167) and S4 (R213) reduced zL (Table I), suggesting that channel opening and voltage sensor activation share a common source of gating charge.
That zL is reduced by voltage sensor mutation is evident by comparing the qa-V relations of D153K, R167A, and R213C to that of the WT in 50 µM Ca2+ (Fig. 6 F). Over a range of negative voltages, qa for these mutants is significantly less than the WT value of zL (dashed line). A decrease in zL is also obvious by comparing the limiting slope of the Po-V relations for R167A (Fig. 6 D) to that of the WT fit (dashed lines). By contrast, mutation of a nonvoltage-sensing residue (Y163E) shifted the qa-V relation to more positive voltages relative to the WT but had no effect on qalim (Fig. 6 F), indicating that zL was unchanged. In the case of WT and Y163E channels, a limiting value of qa was clearly achieved over a wide range of negative voltages, allowing zL to be determined precisely. In general, however, zL is difficult to determine because at the extreme voltages necessary to measure qalim, steady-state recording for prolonged periods is problematic and PO is small. Measurement error is minimized in 50 µM Ca2+ (e.g., Fig. 6 F) since PO is increased more than 1,000-fold relative to 0 Ca2+ (Horrigan and Aldrich, 2002
). However, in some cases, PO could not be measured at sufficiently negative voltages to achieve a limiting slope (e.g., D186A; Fig. 6 E). More commonly, a limiting slope was achieved or closely approached but could not be measured over a wide enough voltage range to determine zL within 0.1 e. By contrast, IK deactivation kinetics can be measured over a much wider range of negative voltages (Horrigan and Aldrich, 2002
). Therefore, the voltage dependence of
(IK) (Fig. 7) was used to estimate the values of zL in Table I.
Macroscopic currents activate or deactivate in response to a voltage step with an exponential time course following a brief delay (Horrigan and Aldrich, 2002
). The mean time constants of IK relaxation (
(IK)) for the WT in 0 Ca2+ and 50 µM Ca2+ are plotted over an extended voltage range in Fig. 7 A. At extreme negative voltages,
(IK) is limited by the rate of channel closing (O0 to C0 in the HCA model; Fig. 3 D) and is well fit by an exponential function from 400 to 150 mV (solid lines) whose slope indicates the partial charge (zN) for closing, and is virtually unaffected by Ca2+. Because zL represents the sum of zN and the partial charge for opening, it is likely that zN will change if zL is altered. The normalized
(IK)-V relations for the WT and several mutants in 50 µM Ca2+ are compared in Fig. 7 B at extreme negative voltages. The slopes of such relations were in many cases indistinguishable from that of the WT (zN; Fig. 7 C). However, mutation of potential voltage-sensing residues (R213, D153, and R167) significantly reduce zN (P < 0.0001; Fig. 7 C), consistent with decreases in zL.
Because zN is easier to determine than the limiting slope of PO, zL for mutants (Table I) were usually estimated based on the fractional change in zN relative to the WT:
 | (9) |
A priori, the ratio of zN to zL need not remain constant. However, zN and qaLim were well correlated (r = 0.96) in cases where both were determined (Fig. 7 C), implying that Eq. 9 is a reasonable approximation. In one case (R167A), the fractional decrease in zN was greater than the decrease in qaLim; but this difference could reflect a failure to achieve a limiting value of qa. For some mutants that behaved similar to the WT (R113Q, D133Q, and R210Q), zN was not measured and zL was set equal to the WT value. For mutations of R207 or R210, zL was constrained by the voltage dependence of Po at positive voltages as described below.
Effects of S4 Mutation on the Equilibrium Constant for Voltage Sensor Activation
Most mutations in S1S4 shifted the Po-V relation to more positive or negative voltages relative to the WT (Fig. 2 D), reflecting changes in the voltage sensor half-activation voltage (VHC(J), Table I). VHC(J) is a function of both voltage sensor charge (zJ) and the intrinsic stability of the activated state, represented by the zero-voltage equilibrium constant J0 (i.e., VHC(J) = (kT/zJ) ln(J0)). Thus, increases in VHC(J) of up to several hundred millivolts produced by mutation of potential voltage-sensing residues R213 and D153 were largely accounted for by decreases in zJ, with little change in J0 (Table I). By contrast, mutation of either R207 or R210 in S4 increased J0 by several orders of magnitude without changing zJ, indicating that the stability of the activated state relative to the resting state was greatly enhanced. These effects on the energetics of voltage sensor activation show that nonvoltage-sensing residues in S4 have an important impact on channel function.
Log(Po)-V relations in 0 Ca2+ for the WT and several mutants that increase Jo (R207Q, R207E, and R210C) are compared in Fig. 8 A. An increase in J0 shifts the steepest part of the curve to more negative voltages without altering its slope, as indicated by the HCA model fits (lines). This effect is most clearly demonstrated for R207Q and R207E because the steepest part of their log(PO)-V relations can be measured, despite shifts of more than 200 mV in their position along the voltage axis relative to the WT. That charge neutralization (R207Q) and charge reversal (R207E) have virtually identical effects suggests that the removal of a positive charge at position 207 was critical for the increase in J0. Gating current measurements confirm that mutation of R207Q enhances voltage sensor activation (Fig. 8 B) (Horrigan and Aldrich, 1999
). The Q-V relation for R207Q, measured with admittance analysis (see MATERIALS AND METHODS) was shifted by 278 mV relative to the WT but is similar in shape, consistent with an increase in JO with no change in zJ.

View larger version (42K):
[in this window]
[in a new window]
|
Figure 8. Effects of S4 mutation on the equilibrium for voltage sensor activation. (A) Mean log(PO)-V relations in 0 Ca2+ for WT( ), R207Q ( ), R207E ( ), and R210C ( ). Mutant curves are shifted to more negative voltages, consistent with an increase in the equilibrium constant for voltage sensor activation (J0) in the HA model (lines, Table I parameters). (B) Normalized Q-V relations for WT and R207Q in 0 Ca2+ measured with admittance analysis from single experiments (see MATERIALS AND METHODS). Curves are fit with Boltzmann functions (WT: z = 0.59 e, Vh = 132 mV; R207Q: z = 0.65 e, Vh = 146 mV). (C) Po-V relations and fits from A are plotted on a linear scale to show that mutants have a similar PO and voltage dependence when voltage sensors are fully activated (i.e., V > 0). (D) qa-V relations for WT and R210C in 0 Ca2+ together with HA model fits showing that the voltage dependence of R210C is consistent with a gating charge similar to the WT. (E) Log(PO)-V relations for R210 mutants in 0 Ca2+ (R210C [ ], R210N [ ], R210E [ ]). R210E was best fit at positive voltages by a Boltzmann function with a partial charge of 0.47 e (solid line). Dashed lines are Boltzmann functions with charge of 0.3 e, representing the value of zL for the WT. (F) (IK)-V relations for WT, R207Q, and R210C. Mutant curves are fit (lines) with a bell-shaped function (IK) = [ 4 + 4]1 at potentials where voltage sensors are fully activated. 4 and 4 represent forward and backward rates for the C4O4 transition in the HCA model (Fig. 3 D) and are exponentially dependent on voltage with partial charges z( 4) and z( 4), respectively.
| |
An important feature of the R207Q data that is reproduced by the model is the weak voltage dependence of Po at positive voltages, evident when the Po-V relation is plotted on a linear scale (Fig. 8 C). For R207Q at V > 0, voltage sensors are fully activated (Fig. 8 B) and PO is limited by the weak voltage dependence of channel opening (C4O4 transition in the HCA model; Fig. 3 D) determined by zL. Therefore the Po-V relation for R207Q is shallower than that of the WT over this voltage range.
In the case of R210C, PO is weakly voltage dependent over a wide voltage range down to 200 mV (Fig. 8 A). A similar observation for R210N led Diaz et. al. (1998) to conclude that R210 is a voltage-sensing residue. However several lines of evidence indicate this is not the case. First, the Po-V relations for R210C and R207Q are virtually identical at V > 0 (Fig. 8 C), implying that the voltage sensors of R210C, like R207Q, are fully activated. Second, when the voltage range of R210C was extended to 300 mV, the voltage dependence of PO increased (Fig. 8 A). Although we could not make measurements at sufficiently negative voltages to determine qaMAX, the qa-V relation for R210C exhibits a steepness similar to that of the WT (Fig. 8 D), consistent with the assumption that zJ is unchanged (dashed fit). These observations suggest that R210C represents an extreme example of the effect represented by R207Q. That is, the voltage sensor equilibrium (J0) for R210C is greatly increased such that the range of voltage sensor activation is shifted to extreme negative voltages, as indicated by the model fits in Fig. 8 (A and D). In this way, R210C appears weakly voltage dependent for V > 200 mV because voltage sensors are constitutively activated.
To confirm that R210 is nonvoltage sensing, we examined the effects of additional mutations at this site (Fig. 8 E). Another charge-neutralizing mutation R210N exhibited reduced Po relative to R210C, but a similar voltage dependence suggested that enhancement of voltage sensor activation does not require a cysteine residue. However a charge reversal (R210E) produced a qualitatively different effect than neutralization. First, the PO-V relation of R210E appears less shifted than R210C and is strongly voltage dependent from 100 to +50 mV. Second, the magnitude of the PO increase between extreme negative and positive voltages is greatly reduced relative to the WT, indicating a decrease in coupling between voltage sensor activation and channel opening (allosteric factor D). As discussed below, the Po-V relation of R210E is consistent with a normal gating charge once changes in D are accounted for. Thus R210 is nonvoltage sensing.
Voltage Sensor Activation May Increase the Charge Associated with Channel Opening
Based on the above results, we conclude that R207 and R210 are nonvoltage sensing and do not contribute to zJ. However, the values of zL determined for mutants at these positions appeared consistently larger than the WT (Table I). We do not think this represents a change in gating charge but rather a difference in the conditions under which zL was determined. For mutants of R207 and R210, zL was adjusted to fit the voltage dependence of PO at positive voltages (e.g., Fig. 8 C) representing the charge to open the channel when voltage sensors are fully activated. For the WT and most mutants, zL was determined from the voltage dependence of Po and
(IK) at extreme negative voltages when voltage sensors are in the resting state. The difference in zL might therefore be explained if the charge associated with channel opening is different when voltage sensors are in the resting or activated state. This hypothesis appears consistent with the Po-V relation for R210E in 0 Ca (Fig. 8 E), representing the only mutant where the weak voltage dependence of PO could be observed at both extreme positive and negative voltages. Boltzmann functions representing WT zL (0.3 e, dashed lines) fit the PO-V relation of R210E at negative voltages but not at positive voltages where the data are best fit with a charge of 0.47 e (solid line). Thus it appears that zL increases when voltage sensors are activated.
To confirm the value of zL for mutants of R207 and R210 we also examined the voltage dependence of
(IK) (Fig. 8 F). When voltage sensors are fully activated, channel gating in 0 Ca2+ should be approximated by a two-state model (C4O4) and the
(IK)-V relation is described by bell-shaped function
(IK) = [
4+
4]1, where
4 and
4 represent the voltage-dependent forward and backward rate constants (Horrigan et al., 1999
). The
(IK)-V relations for R207Q and R210C are well fit by this function (Fig. 8 F); and the summed partial charges associated with the forward and backward rates yield values from 0.34 to 0.58 e for different mutants (R207Q, 0.51 e; R207E, 0.58 e; R210C, 0.43 e; R210E, 0.34 e; R210N, 0.48 e). These estimates are comparable to values of zL determine from Po (Table I), but are probably not as accurate because the data must be fit at voltages near peak
(IK) where neither z(
4) nor z(
4) are well constrained.
Effects on Coupling between Voltage Sensor Activation and Channel Opening
Because decreases in the voltage dependence of Po can be caused by decreases in allosteric coupling, we were careful to examine the effects of mutation on the coupling factor D. In most cases, large voltage-dependent increases in PO were observed in 0 Ca2+, indicating that D was unchanged. However, three mutations (R210E, E180A, and D186A) reduced D significantly (Table I). Two of these in S3 (E180A, D186A) represent relatively small changes (twofold) that have a small impact on qamax (<15%). By contrast, R210E exhibited a 30% reduction in qamax that we attribute entirely to a sevenfold decrease in D.
Both R210E and E180A decreased the coupling factor D without altering gating charge. This is evident by comparing the log(PO)-V relations of these mutants to the WT in 0 Ca2+ (Fig. 9 A). The foot of the curve for mutants and WT is superimposable, unlike mutants of potential voltage-sensing residues (Fig. 5 F), indicating that gating charge is not altered. However the magnitude of the voltage-dependent Po increase is reduced relative to the WT, especially in the case of R210E. PO-V relations for R210E in different [Ca2+] are plotted in Fig. 9 B together with fits to the WT data (dashed lines). The weak voltage dependence of R210E at both positive and negative voltages is consistent with complete activation and deactivation of voltage sensors, respectively. Between these extremes, Po increases sharply over a 150 mV range like WT, consistent with a normal voltage sensor charge, but the magnitude of the Po increase is roughly 1,000-fold smaller than WT. In 0 Ca2+, Po for R210E is less than the WT at positive voltages but is
100-fold greater than the WT at extreme negative voltages. These differences can be reproduced with the HA model (solid lines) by a sevenfold reduction in D and a 100-fold increase in LO. The reduction in D is necessary to account for the reduced magnitude of the voltage-dependent Po increase and also accounts for the decrease in qamax (Fig. 9 C). The qa-V relation for R210E is narrower than that of the WT but has similar steepness at negative voltages, supporting the conclusion that coupling is reduced without altering gating charge.
The ability of R210C and R210N to enhance voltage sensor activation prevented us from determining whether these mutants reduce coupling; we could not measure Po when voltage sensors were in the resting state nor determine qaMAX. The qa-V relation for R210C can be fit with a WT value of D (Fig. 8 D) but is also consistent with a function like that describing R210E (Fig. 9 C). One reason to suspect that R210C does not alter coupling is that its PO in 0 Ca2+ closely matches that of R207Q and R207E when voltage sensors are fully activated (Fig. 8 A). R207 mutants exhibit a large decrease in Po at negative voltages, indicating a value of D similar to the WT. By the same logic, R210N may decrease coupling because Po is reduced relative to the WT, like R210E at