Allosteric Voltage Gating of Potassium Channels II: mSlo Channel Gating Charge Movement in the Absence of Ca2+

Frank T. Horrigana and Richard W. Aldricha
a From the Department of Molecular and Cellular Physiology, Howard Hughes Medical Institute, Stanford University School of Medicine, Stanford, California 94305

Correspondence to: Richard W. Aldrich, Department of Molecular and Cellular Physiology, Howard Hughes Medical Institute, Stanford University School of Medicine, Stanford, CA 94305., raldrich{at}leland.stanford.edu (E-mail), Fax: 650-725-4463; (fax)


  Abstract
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
Acknowledgements
References

Large-conductance Ca2+-activated K+ channels can be activated by membrane voltage in the absence of Ca2+ binding, indicating that these channels contain an intrinsic voltage sensor. The properties of this voltage sensor and its relationship to channel activation were examined by studying gating charge movement from mSlo Ca2+-activated K+ channels in the virtual absence of Ca2+ (<1 nM). Charge movement was measured in response to voltage steps or sinusoidal voltage commands. The charge–voltage relationship (Q–V) is shallower and shifted to more negative voltages than the voltage-dependent open probability (G–V). Both ON and OFF gating currents evoked by brief (0.5-ms) voltage pulses appear to decay rapidly ({tau}ON = 60 µs at +200 mV, {tau}OFF = 16 µs at -80 mV). However, QOFF increases slowly with pulse duration, indicating that a large fraction of ON charge develops with a time course comparable to that of IK activation. The slow onset of this gating charge prevents its detection as a component of IgON, although it represents ~40% of the total charge moved at +140 mV. The decay of IgOFF is slowed after depolarizations that open mSlo channels. Yet, the majority of open channel charge relaxation is too rapid to be limited by channel closing. These results can be understood in terms of the allosteric voltage-gating scheme developed in the preceding paper (Horrigan, F.T., J. Cui, and R.W. Aldrich. 1999. J. Gen. Physiol. 114:277–304). The model contains five open (O) and five closed (C) states arranged in parallel, and the kinetic and steady-state properties of mSlo gating currents exhibit multiple components associated with C–C, O–O, and C–O transitions.

Key Words: calcium, potassium channel, BK channel, ion channel gating, gating current


  Introduction
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
Acknowledgements
References

Large-conductance Ca2+-activated K+ channels (BK channels)1 are sensitive to membrane potential as well as intracellular calcium. Although the voltage dependence of these channels is weak compared with that of many purely voltage-gated K+ (Kv) channels (Barrett et al. 1982 ; Cui et al. 1997 ), voltage gating is likely to be of central mechanistic importance to BK channel function. Sequence similarity between the Slo family of BK channels and Kv channels, including the presence of a charged 17 domain, suggests that the basic structure of the BK channel may resemble that of Kv channels, and that Ca2+ acts to regulate the function of this voltage-dependent "core" (Cui et al. 1997 ; Meera et al. 1996 ; Wei et al. 1994 ). Consistent with this hypothesis, BK channels can be maximally activated by voltage in the absence of Ca2+ (Cox et al. 1997 ; Cui et al. 1997 ), gating currents are detected under these conditions (Stefani et al. 1997 ), and mutations in the 17 region alter the voltage dependence of activation (Diaz et al. 1998 ; Cui, J., and R.W. Aldrich, manuscript in preparation). Thus, BK channel voltage sensitivity reflects the action of an intrinsic voltage sensor, and Ca2+ is not required for channel opening. Ca2+ shifts many of the voltage-dependent parameters of BK channel gating to more negative voltages (Barrett et al. 1982 ; Cox et al. 1997 ). These results are consistent with a model in which Ca2+ binding allosterically regulates voltage-dependent channel activation (Cox et al. 1997 ; Cui et al. 1997 ). Such a mechanism implies that voltage-dependent gating plays a critical role in the response of BK channels to both Ca2+ and voltage.

In the preceding article (Horrigan et al. 1999 ), we examined the response of mSlo Ca2+-activated K+ channels to voltage by recording K+ current in the absence of Ca2+. The kinetic and steady-state properties of mSlo IK indicate that the mechanism of voltage gating is complex. In response to a voltage step, IK activates with an exponential time course after a brief, voltage-dependent delay. The exponential relaxation of IK suggests that a rate-limiting step dominates channel activation. However, the delay indicates that rapid voltage-dependent transitions also exist in the activation pathway. The time constant of IK relaxation ({tau}(IK)) and steady-state open probability (Po) both exhibit complex voltage dependencies that are inconsistent with many conventional sequential gating schemes. A particularly important finding is that {tau}(IK) and Po become less voltage dependent at very negative voltages. To account for these results, we proposed a voltage-gating scheme based on an allosteric mechanism.

This diagram illustrates the idea that mSlo channels undergo a transition between a closed (C) and open (O) conformation, and that this transition is influenced allosterically by the state of four independent and identical voltage sensors (one for each subunit). We assume each voltage sensor can undergo a transition between a resting (R) and an activated (A) conformation, and the equilibrium constant for the C–O transition (L) increases by a constant factor (D) for each voltage sensor that is activated. Similarly, the equilibrium constant for voltage sensor activation (J) increases D-fold in favor of the activated state, when the channel opens. Therefore, the factor D embodies the allosteric interaction between voltage-sensor activation and channel opening. This mechanism results in a gating scheme (Figure 144) that contains a parallel arrangement of open and closed states.



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 1. mSlo gating capacitance. (A) Voltage-dependent membrane capacitance (Cg) is plotted versus voltage for patches excised from HEK 293 cells that expressed mSlo channels (mSlo) or were not transfected (control). Each trace represents the average response to two voltage ramps (see panel C). The baseline (Cg = 0) was set to the mean of Cm between -150 and -100 mV. The Cg–V relationship for mSlo is fit by the derivative of a Boltzmann function with respect to voltage (dashed line) with z = 0.59 e and Vh = 133 mV. (B) The Qg–V relationship, determined by integrating the Cg–V, is fit by a Boltzmann function. The shape of the average GK–V relationship (symbols; from Horrigan et al. 1999 ) can be approximated by raising the Qg–V fit to the 4th power (Fit4). (C) The voltage protocol used to measure Cg in panel A consists of a 1-s voltage ramp from -160 to +200 mV superimposed with a sinusoidal command (30 mV amplitude, 868 Hz). Capacitance was determined for each period of the sinusoid. (D) Cg–V relationships obtained from channels expressed in a Xenopus oocyte using sinusoidal voltage amplitudes of 3 mV (symbols) or 30 mV (solid line) at 1,736 Hz are superimposable, indicating that sinwave amplitude does not effect the measurement. The 3 and 30 mV traces represent the average response to eight and four voltage ramps, respectively. Data are plotted at 3-mV intervals, representing the mean capacitance over eight periods of the sinusoid. The internal solution did not contain crown ether (18C6TA).



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 2. Cg represents mSlo charge movement. (A) The Ca2+ sensitivity of Cg is illustrated by comparing Cg–V relationships obtained in 0 or 60 µM [Ca2+]i from the same patch at 868 Hz. Peak Cg in 0 Ca was 39% larger than in 60 Ca2+, but traces were normalized to show the shift in the position of peak Cg along the voltage axis. (B) Cg–V relationships for mSlo wild-type (WT, 1,736 Hz) and mutant (R207Q, 1,500 Hz) channels are similar in shape but are shifted by more than 250 mV relative to each other. Fits, representing the derivative of Boltzmann functions, are superimposed on the data (dashed lines; WT: z = 0.58 e, Vh = 132 mV; R207Q: z = 0.74, Vh = -143 mV). (C) Cg–V relationships for wild-type channels decrease in amplitude as the frequency of the sinusoidal voltage command is increased from 200 to 6,944 Hz. (D) The orthogonal component of the admittance signal (Gg) increases with frequency. (E) DC current measured during the voltage ramp (see Materials and Methods) changes linearly with voltage, demonstrating that the G signal does not represent voltage-dependent changes in Gm and indicating a membrane/seal resistance (125 GV) over the entire voltage range. (F) Cg (solid symbols) and Gg (open symbols) measured at +120 mV are plotted versus frequency for two experiments and are fit by Lorenzian functions, described in the text, with a time constant of 70 µs.



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 3. mSlo gating current. (A) mSlo Ig evoked in response to a 0.5-ms pulse to +160 mV from a holding potential of -80 mV. The trace represents the averaged response to eight pulses. IgON and IgOFF are fit by exponential functions (dashed lines). (B) A family of Ig evoked in response to 1-ms pulses to different voltages (0–160 mV in 40-mV steps). (C) QON–V and QOFF–V relationships were obtained by integrating IgON and IgOFF, respectively (from B) over 1-ms intervals. The Qg–V relationship was obtained from Cg measurements at 868 Hz in the same patch. Qfast was determined from an exponential fit to IgON (see below). (D) IgON evoked at +160 mV is compared with the initial time course of IK activation measured at the same voltage from a different experiment. IK is fit with an exponential function (dashed line) from 0.5 to 20 ms after the start of the pulse. The IK scale bar represents 10% of the steady-state amplitude. IgON is also fit with an exponential function, and the shaded area under the fit was used to determine Qfast. (E) IgON and Cg measured from a single patch were integrated to determine Qfast and Qg, respectively, as plotted in F. The Qfast–V relationship is fit with a Boltzmann function (z = 0.57 e, Vh = 136 mV).



View larger version (27K):
[in this window]
[in a new window]
 
Fig. 4. Voltage dependence and kinetics of fast charge movement. (A1) The normalized Qfast–V relationships for many experiments are fit with Boltzmann functions (z = 0.59 e, dashed lines). The solid line is a Boltzmann function indicating the mean half-activation voltage (<Vh> = 155 mV, z = 0.59 e). (A2) The data from A1 (open symbols) are aligned by shifting them along the voltage axis by {Delta}Vh = (<Vh> - Vh). The mean Qfast–V (filled circles, mean ± SEM) is superimposed on the data together with two Boltzmann fits (Vh = 155 mV; solid line: z = 0.59e, dashed line: z = 0.55 e) and was determined by averaging the shifted data in 15-mV bins. (B1) Time constants of fast Ig relaxation ({tau}gFast) were determined from exponential fits to ON and OFF currents for the experiments in A and are plotted on a log scale versus voltage. (C1) Three {tau}gFast–V relationships from B1 that cover a large voltage range are compared. (B2 and C2) Data from B1 and C1 were shifted along the voltage axis by {Delta}Vh (determined from A) and then normalized to the mean {tau}gFast measured from +100 to +180 mV (59 µs). The solid line in B2 indicates the best fit of a two-state model of voltage-sensor activation where the relationship between the forward ({alpha}) and backward (ß) rates are constrained such that J = {alpha}/ß = 1 at +155 mV (z{alpha} = +0.30 e, zß = -0.21 e, {alpha}(0) = 1,310 s-1, ß(0) = 30,160 s-1). Dashed lines in A2, B2, and C2 represent the parameters ultimately used in the allosteric model to describe closed-state charge movement (z{alpha} = +0.33 e, zß = -0.22 e, {alpha}(0) = 1,100 s-1, ß(0) = 32,120 s-1).



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 5. Slow component of gating-charge movement. A family of Ig was evoked at +140 mV in response to pulses of different duration (0.06–20 ms). (A) Plots the records for 0.06–2-ms pulse duration. The remaining records are shown in Figure 6 A. (B) QOFF was determined by integrating IgOFF for 3 ms after each voltage pulse and is plotted versus pulse duration (Qp). Qp(t) is fit by a double-exponential function with time constants {tau}gFast = 63 µs and {tau}gSlow = 4.22 ms. {tau}gFast was determined by fitting IgON, and QpFast was set equal to Qfast (11.67 fF) determined as in Figure 3 D. (C) The time derivative of the fit to Qp(t) (Qp', dashed line) superimposes on the time course of IgON at +140 mV.



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 6. Slow charge movement is limited by channel activation. (A) Families of Ig evoked by pulses of different duration (0.06–20 ms) to the indicated voltages (HP = -80). Records obtained at +184 and +224 mV were from a different patch than those obtained at +100 and +140 mV. Scale bars represent 100 pA. (B) Qp(t) curves at different voltages from three experiments are fit with double-exponential functions where the time constant and amplitude of the fast component were determined by fitting IgON with an exponential function. Plots from different experiments (different symbols) were normalized to the total fast charge (QTfast) for each patch (see text). The indicated pulse voltages have been corrected based on the Q–V shifts ({Delta}Vh) determined in Figure 4 A for each experiment. (C) The time constants of QpSlow relaxation ({tau}gSlow) (mean ± SD) are plotted versus voltage and compared with the time constants for IK relaxation (mean ± SD) from Horrigan et al., 1999. {tau}gSlow at +224 mV represents a single measurement. (D) Normalized steady-state Q–Vs (open circles) from four experiments were determined with 20-ms pulses and are fit by a Boltzmann function (dashed line, Vh = 143 mV, z = 0.65 e). Qss–V relationships were corrected for {Delta}Vh determined for these experiments in Figure 4 A. Averaged Qfast–V and GK–V relationships are plotted for comparison (filled symbols) and are fit by the allosteric voltage-gating scheme (solid lines, zJ = 0.55 e, Vh(J) = 155 mV, L = 2 x 10-6, zL = 0.4 e, D = 17).



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 7. Changes in OFF kinetics with channel activation. (A) A family of IgOFF evoked at 2100 mV after pulses to +140 mV of 0.06–20 ms duration (from Figure 6 A). Current amplitude is maximal after a 0.5-ms pulse, but IgOFF decays more slowly as pulse duration increases. The baseline for each record is set to the mean current during an interval 4–5 ms after the pulse. (B) The decay of OFF currents are fit by double-exponential functions with {tau}F = 15.5 µs and {tau}M = 59 µs. (C) QOFF obtained by integrating IgOFF from A achieves a steady state within 300 µs after a 0.06-ms pulse (arrow) but relaxes more slowly after longer pulses. (D) The kinetics of QOFF relaxation after a brief (0.06 ms) or prolonged (10–20 ms) pulses are compared by plotting QOFFQss on a semilog scale. Qss is the steady-state value of QOFF measured 3 ms after the pulse. The 0.06-ms trace is fit by a single-exponential function ({tau}F = 15.5 µs). The 10–20-ms trace, representing an average of 10-, 15-, and 20-ms records, is fit by a triple exponential (solid line, {tau}F = 15.5 µs, {tau}M = 59 µs, {tau}S = 448 µs) where the individual components are indicated by dashed lines. (E) A family of QOFFQss for the data in C. Traces are fit with triple exponential functions with the time constants determined from D. (F) QOFF component amplitudes from these fits are plotted versus pulse duration. The relaxation of all three components is fit by exponential functions (solid lines) with a time constant of 4.22 ms. Error bars represent the component amplitudes obtained when {tau}M is changed by ±10% (with {tau}F and {tau}S held constant). The fast component of ON charge (Qfast) is indicated by an arrow. (G) Fast and Medium IgOFF component amplitudes determined from B are plotted versus pulse duration. Solid lines represent exponential fits with a time constant of 4.22 ms. (H) The allosteric model predicts three components of QOFF relaxation corresponding to the indicated transitions in the gating scheme.



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 8. Predictions of the allosteric model. (A) QOFF component amplitudes determined after pulses to +240 mV in 0 Ca2+ are plotted versus pulse duration. The fast component is reduced to <10% of the total OFF charge after a 20-ms pulse. The relaxation of all three components is fit by exponential functions (solid lines) with a {tau} = 0.91 ms. (B) The decay of QOFFQss is plotted on a semilog scale after 0.1- or 20-ms pulses to +160 mV in 60 µM Ca2+ (HP = -80). The 0.1-ms trace is fit by a triple exponential function (solid line, {tau}F = 23.8 µs, {tau}M = 150 µs, {tau}S = 822 µs) and the 20-ms trace is fit with a double-exponential ({tau}M = 150 µs, {tau}S = 822 µs), indicating that the fast component is eliminated when most channels are opened. Dashed lines represent the two components of the 20-ms fit and the fast component of the 0.1-ms fit. (C) Normalized QOFF component amplitudes and total OFF charge (Qp) are plotted versus pulse duration for pulses to +160 mV in 0 Ca2+. OFF components were measured upon repolarization to -80 mV and are normalized to the fast component of ON charge (Qfast) at +160 mV. (D) When QOFF is measured upon repolarization to 0 mV, the Fast component and Qp are unchanged. However, the Medium component decreases and the Slow component increases in a complementary manner. (E) The charge distributions predicted by the allosteric model for Closed (QC) and Open channels (QO) are plotted versus voltage (zJ = 0.55 e, Vh(J) = 155 mV, L = 2 x 10-6, zL = 0.4 e, D = 17). Arrows indicate the predicted amplitudes of Medium and Slow OFF components at repolarization voltages of -80 and 0 mV after a pulse to +160 mV (VP).



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 9. Simulations of Fast Ig. (A) A family of IgON evoked at different voltages (0 to +140 mV) is compared with the prediction of the allosteric scheme (solid lines). Data and simulated traces were both evoked in response to filtered voltage pulses (20 kHz) and then filtered at 20 kHz. (B) A family of gating currents evoked at +140 mV in response to pulses of different duration (from Figure 5 A) is fit by the allosteric model (solid lines). Model parameters for panels A and B are as shown in Table 1 (Case A) with the exception that {alpha} and ß were decreased by 2% ({alpha}(0) = 1,080 s-1, ß(0) = 31,681 s-1) to match this experiment. (C) {tau}gFast measured from simulated traces at different voltages is plotted versus voltage (filled circles) and compared with the {tau}F–V relationship predicted from the parameters assigned to the R–A transition in the model (solid line, {tau} = 1/({alpha}+ß); Case A in Table 1). Open symbols indicate the time constant of the Medium OFF component ({tau}M) measured from several patches. Lines through these data represent predictions of the allosteric scheme (see text). (D) The Qfast–V relationship measured from simulated currents (symbols) is compared with the QC–V relationship specified by the model (line). (E) The time course of Qp predicted by the allosteric model (lines) accounts for the fast component of ON charge but underestimates the magnitude of the slow component.



View larger version (32K):
[in this window]
[in a new window]
 
Fig. 10. Estimating open probability from charge movement. The allosteric model predicts a close relationship between Po and the various QOFF components. A, B, and C plot the voltage dependence of these components measured after 20-ms pulses for three experiments. Solid lines indicate predictions of three models (Cases A, B, and C) described in the text. (A) The Fast OFF component should be proportional to the number of closed channels at the end of the pulse. Therefore, [1 - [QOFFfast(VP)/Qfast(VP)] is plotted as an estimate of steady-state Po, where Qfast is the fast component of ON charge. (B) The Slow OFF component should be directly proportional to Po. The quantity (QOFFslow(VP)/QTfast) is plotted where QTfast is the total fast charge estimated by fitting the Qfast–V relationship with a Boltzmann function. (C) The Medium OFF component is normalized by QTfast and plotted versus voltage. (D) The voltage dependence of the Medium OFF component was also examine by fitting IgOFF with a double-exponential function ({tau}F, {tau}M) and plotting the normalized amplitude of IgOFFmed against voltage. IgOFFmed was normalized by fitting the IgOFF–V relationship with a Boltzmann function corresponding to Case C (z = 0.98 e). (E) The Slow component of ON charge (QpSlow) is expected to exhibit a complex voltage dependence (dashed curves) that is highly sensitive to Po. The data, normalized by QTfast, indicate that the slow component is too large to be accounted for by the initial allosteric model parameters (Case A) but a shift in the Po–V relationship (Cases B and C) produces a better fit. A dashed line indicates the charge assigned to the C–O transition (zL) (F) {tau}gSlow determined from the time course of QpSlow and QOFFslow for many experiments are plotted versus voltage. Solid symbols represent mean ± SEM. Dashed and solid lines represent predictions of Case A and B, respectively (Table 1).



View larger version (35K):
[in this window]
[in a new window]
 
Fig. 11. Ig Simulations using modified parameters. (A) Ig evoked at +140 mV in response to pulses of different duration (HP = -100 mV) are plotted on two different time scales and compared with the predictions of the allosteric scheme using modified parameters corresponding to Case B in Figure 10 (solid lines, Table 1: Case B). (B) Ig evoked at +224 mV is also fit by the model. (C) QOFFQss corresponding to the records in A are plotted on a semilog scale, demonstrating that the model reproduces all three components of QOFF. (D) Qp(t) from Figure 6 B is reproduced by the model (QON(t)) for V >= +140 mV but the slow component amplitude is underestimated at lower voltages. (E) The Cg–V relationship measured from the same patch as in A at 868 Hz is compared with simulations of the allosteric scheme corresponding to Case A and Case B parameters (solid lines). Dashed lines indicate the derivatives of the Qc–V and Qo–V relationships (Q 'C, Q 'O). The scale bar represents 50 fF.



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 12. Predictions of a sequential gating scheme. Predictions of Figure 17 (dashed lines) are compared with data (symbols) and predictions of the allosteric scheme (solid lines). (A) The GK–V relationship can be fit by either model (except at very low Po [ Horrigan et al. 1999 ]), while assigning identical parameters for voltage-sensor movement (zJ = 0.55 e, Vh(J) = 155 mV) and similar charge to the CO transition (z{isin} = 0.315 e: Scheme IV; zL = 0.40 e: Allosteric Scheme). Figure 17 can also approximate the G–V relationship predicted by the allosteric scheme corresponding to Case B by increasing the C–O equilibrium constant ({isin}(0) = 0.42 Case A; {isin}(0) = 4.34 Case B). (B) Both models reproduce the fast and slow components of Qp(t) at +140 mV (from Figure 5B). Rate constants for the allosteric model correspond to those used in Figure 11 (Table 1, Case B). Figure 17 used the same rates ({alpha}, ß) and charge (z{alpha}, zß) to specify closed-state transitions. The forward and backward rate constants for the C–O transition were {delta}(0) = 1,392 s-1 and {gamma}(0) = 322 s-1, respectively, and were assumed symmetrically voltage dependent (z{delta} = +0.158 e, zg = 20.158 e). (C) Both models predict similar bell-shaped QpSlow–V relationships, implying that a large fraction of slow charge arises from voltage-sensor movement rather than the C–O transitions. (D) The instantaneous Ics–V relationship for mSlo was measured in symmetrical 110 mM Cs+ solutions containing no added K+ by activating channels in response to a 50-ms pulse to +200 mV and then stepping to various voltages. ICs was measured 100 µs after the pulse to avoid contamination by IgOFF.



View larger version (38K):
[in this window]
[in a new window]
 
Fig. 14. Scheme I

The horizontal transitions (C–C and O–O) reflect the activation or movement of voltage sensors while vertical (C–O) transitions represent channel opening. The closed and open conformations are each represented by five states, with subscripts (0–4) denoting the number of activated voltage sensors.

For this scheme to reproduce IK, it is necessary that voltage-sensor activation is fast and accounts for most of the channel's voltage dependence while C–O transitions are slow and weakly voltage dependent (Horrigan et al. 1999 ). Closed-state transitions (C–C) must be fast and voltage dependent (zJ = 0.55 e per voltage sensor) to describe the delay in IK activation. C–O transitions must be slow to limit the exponential relaxation of IK. The weak voltage dependence of {tau}(IK) and Po at negative voltages implies that the charge associated with channel opening is small (zL = 0.4 e). Finally, the equilibrium constant L is small (~10-6) and the allosteric factor large (D = 17), equivalent to an interaction energy of 2.8 kT, to account for the shape of the Po–V relationship in 0 Ca2+.

Such a model provides mechanistic insight and places constraints on the possible molecular events that link voltage-sensor movement and channel opening (Horrigan et al. 1999 ). The allosteric relationship between voltage-sensor activation and channel opening requires that the channels can open with any number of voltage sensors activated, including none. Furthermore, it requires that the allosteric transitions from closed to open alter the energetics of voltage-sensor movement such that voltage sensors, present in each subunit, are more easily activated when the channel is open. This effect is analogous to the change in ligand affinity that occurs between the T and R states in an allosteric ligand binding model (Monod et al. 1965 ).

Although the properties of mSlo IK are consistent with the allosteric model, several aspects of the gating scheme are not tightly constrained by the ionic current data (Horrigan et al. 1999 ). Transitions among closed and open states (C–C, O–O) do not immediately alter the number of open channels and, therefore, are not observed directly as a change in IK. Instead, these transitions contribute to the delay in IK activation and to the complex voltage dependence of IK kinetics and steady-state activation. However, any voltage-dependent transition must produce a movement of gating charge that can be detected as gating current (Ig). Gating current provides a direct assay of voltage-sensor movement (C–C, O–O transitions) and, therefore, constrains any voltage-dependent gating scheme. The allosteric model makes specific predictions about the kinetic and steady-state properties of gating charge movement and their relationship to IK. Our experiments examine these predictions and provide a critical test of the model.

Our results are consistent with the assumption that mSlo voltage sensors move rapidly and independently while channels are open or closed. Measurements of the charge associated with voltage-sensor movement are in line with previous estimates based on the ionic current data. Our results also support the prediction that channel opening alters the kinetics of voltage-sensor movement. Finally, we show that some complex kinetic and steady-state properties of mSlo charge movement are reproduced by the proposed gating scheme. These include a large slow component of ON charge that is limited by the speed of channel opening, and three components of OFF charge reflecting C–C, O–O, and C–O transitions. The relationships between these components are consistent with the allosteric model and rule out many alternative schemes.


  Materials and Methods
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
Acknowledgements
References

Channel Expression
Experiments were performed with the mbr5 clone of the mouse homologue of the Slo gene (mSlo), kindly provided by Dr. Larry Salkoff (Washington University School of Medicine, St. Louis, MO). The clone was modified to facilitate mutagenesis and was propagated and cRNA transcribed as described previously by Cox et al. 1997 . Xenopus oocytes were injected with ~50 ng of cRNA (50 nl, 1 ng/nl) 3–7 d before recording. mSlo was also subcloned into a mammalian expression vector (SR{alpha}; kindly provided by Dr. A.P. Braun, University of Calgary, Calgary, Alberta, Canada) containing the SV40 promoter. HEK 293 cells expressing the large T-antigen of the SV40 virus were cotransfected with mSlo and green fluorescent protein (GFP, as a marker) using LipofectAMINE (GIBCO BRL) 3 d before recording.

Electrophysiology
Currents were recorded using the patch clamp technique in the inside-out configuration (Hamill et al. 1981 ). Upon excision, patches were transferred to a separate chamber and washed with at least 20 vol of solution. The internal solution contained (in mM) 135 N-methyl-D-glucamine (NMDG)-MeSO3, 6 NMDG-Cl, 20 HEPES. 40 µM (+)-18-crown-6-tetracarboxylic acid (18C6TA) was added to chelate contaminant Ba2+ (Diaz et al. 1996 ; Neyton 1996 ; Cox et al. 1997 ) unless otherwise indicated. In addition, "0 Ca2+" solutions contained 2 mM EGTA, reducing free Ca2+ to an estimated 0.8 nM in the presence of ~10 µM contaminant Ca2+ (Cox et al. 1997 ). These solutions are considered Ca2+-free for the purposes of this study since [Ca2+]i < 50 nM does not affect Slo channel activation (Cui et al. 1997 ; Meera et al. 1996 ). Solutions containing 60 µM Ca2+ were buffered with 1 mM HEDTA, and free Ca2+ was measured with a Ca2+ electrode (Orion Research, Inc.). The external (pipette) solution contained 125 tetraethylammonium (TEA)-MeSO3, 2 TEA-Cl, 2 MgCl2, 20 HEPES. pH was adjusted to 7.2. Solutions containing 110 mM K+ were as described in the preceding article (Horrigan et al. 1999 ). Experiments were performed at room temperature (20–22°C).

Measurement of rapid gating current in response to voltage pulses requires accurate subtraction of linear capacitive currents due to the electrode and cell membrane. Electrodes were pulled from thick-walled 1010 glass (World Precision Instruments) and coated with wax (sticky wax; Kerr) to minimize electrode capacitance (~1 pF). Pipette access resistance (Rs) ranged from 0.7 to 1.5 MV in K-free solutions. Membrane capacitance ranged from 0.25 to 1 pF as determined by the responses to a -10 mV voltage step from -80 mV before and after sealing the electrode tip onto Sylgard (Dow Corning). Data were acquired with an Axopatch 200B amplifier (Axon Instruments, Inc.) in patch mode at a relatively low gain (1–2 mV/pA) to avoid saturation of capacitive transients in response to voltage steps that often exceeded 300 mV. Both the voltage command and current output were filtered at 20 kHz with 8-pole bessel filters (Frequency Devices, Inc.) to limit the speed of fast capacitive transients so that they could be accurately sampled and subtracted. The Axopatch's internal filter was set at 100 kHz. Currents were sampled at 100 kHz with a 16 bit A/D converter (ITC-16; Instrutech). Ig records were typically signal-averaged in response to at least eight voltage pulses, and a P/-4 protocol was used for leak subtraction (Armstrong and Bezanilla 1974 ) from a holding potential of -80 mV. 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 squares fits.

Simulations
Simulations were calculated at 1-µs intervals using a fifth order Runga-Kutta algorithm with adaptive step size (Press et al. 1992 ) implemented in Igor Pro (Wavemetrics, Inc.). Voltage commands and simulated currents were convolved with the impulse response of a 20 kHz 8-pole bessel filter to reproduce the experimental condition (see Horrigan et al. 1999 ).

Admittance Analysis
Admittance (Y) is defined by the expression Y = I/V where V and I represent the amplitude of the sinusoidal voltage command and resultant current at a specific frequency (f). The admittance of a membrane (Ym) is:

(1)

where Gm and Cm are membrane conductance and capacitance, respectively (j=, {omega} = 2{pi}f). Gm and Cm each represent the sum of a contribution from the lipid bilayer (Gb, Cb) and gating charge movement (Gg(V), Cg (V)) (see RESULTS):


The total admittance of the patch equivalent circuit is:

(2)

where Cs is the stray capacitance of the electrode and holder, and Rs is the series resistance. Combining (1) and (2):

(3)

where

Under typical experimental conditions [RS {cong} 106 , Cm {cong} 1 pF, Gm < 1 nS, {omega} = 5,451 (f = 868 Hz)], T approaches unity, and this expression can be approximated:

(4)

Therefore, Cg(V) can be determined directly as the voltage-dependent component of Yp/{omega} appearing at a phase angle of 90° relative to the command voltage.

For admittance measurements, the membrane was clamped with a sinusoidal voltage command (60 mV peak to peak) generated by the D/A converter of the ITC-16 interface at 18-µs intervals (at least eight samples per cycle of the sinusoid). The voltage command and current signal were both filtered at 20 kHz. Admittance was determined for each cycle of the sinusoid at 0° and 90° after correcting for phase shifts ({Delta}{Phi}) due to the filters and amplifier. These were determined at each frequency by measuring the admittance of an electrode in solution (Y = 1/RS), which should appear at an angle of {Delta}{Phi} relative to the command voltage. DC current was determined as the mean current over each cycle of the sinusoid.


  Results
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
Acknowledgements
References

Gating charge movement was examined in excised macropatches from Xenopus oocytes and HEK 293 cells expressing the pore-forming {alpha} subunit of mSlo Ca2+-activated K channels. Several factors combine to make gating currents more difficult to measure for mSlo than for Kv channels such as Shaker. First, mSlo is less voltage dependent than Shaker, and the gating charge is correspondingly smaller. The steady-state GK–V relationships for mSlo in 0 Ca2+ and Shaker can be approximately fit by Boltzmann functions with equivalent charges of 1 e (Cui et al. 1997 ) and 5.3e (Zagotta et al. 1994b ), respectively. A similar difference is observed based on more sophisticated estimates of total gating charge per channel, ranging from 2.6 to 4.4 e for Slo channels (Horrigan et al. 1999 ; Stefani et al. 1997 ) and from 12 to 16e for Shaker (Aggarwal and MacKinnon 1996 ; Schoppa et al. 1992 ; Zagotta et al. 1994b ). In addition to this 5-fold difference in charge, mSlo exhibits a single channel conductance roughly 10-fold greater than that of Shaker. Thus, the ratio of gating charge to ionic current is almost two orders of magnitude smaller for mSlo than for Shaker. Owing to this relationship, measurement of mSlo gating current (Ig) requires stringent conditions for eliminating IK. For example, in this study, Ig was recorded in the presence of isotonic (125 mM) external TEA, to block the channel pore, even though internal and external solutions were nominally K+-free. Experiments performed in the absence of TEA appear to produce similar results (Stefani et al. 1997 ) but require extensive washing to assure that residual IK is eliminated. Even when TEA was present at a concentration almost 1,000-fold higher than the Ki for mSlo block (Ki = 0.14 mM; Butler et al. 1993 ), Ig could not be recorded in the presence of normal internal K+ (110 mM) because IK was not abolished. Finally, mSlo channels activate only at very positive voltages in 0 Ca2+ with a half-activation voltage of +190 mV for GK (Cui et al. 1997 ; Horrigan et al. 1999 ). The high voltages and large voltage steps needed to activate these channels often proved problematic for leak-free recording of small gating currents.

Two approaches were used to measure mSlo gating charge movement. The first involved clamping the membrane with a sinusoidal voltage command and measuring gating charge as a voltage-dependent component of membrane capacitance using admittance analysis (Fernandez et al. 1982 ). The second involved conventional measurement of gating currents in response to voltage steps. Although the bulk of the analysis was performed using voltage steps, the admittance analysis is presented briefly first (see Figure 1 and Figure 2) to provide an initial characterization of mSlo charge movement and to demonstrate several necessary controls.

Gating Capacitance Measurements and Admittance Analysis
Membrane capacitance (Cm) represents the ability of charge to redistribute across or within the cell membrane in response to a change in voltage. Therefore, Cm includes a nonlinear voltage-dependent contribution from gating charge movement (Cg) as well as a voltage-independent component due to the lipid bilayer. One of the most sensitive methods for measuring capacitance is admittance analysis. The membrane is driven with a sinusoidal voltage command, and the resulting current is analyzed with a phase-sensitive detector to determine Cm as well as other parameters in the membrane equivalent circuit (Gillis 1995 ; Lindau and Neher 1988 ). One advantage of this technique is that residual ionic and "leak" currents can be separated from capacitive (gating) currents based on their phase relative to the command voltage. A related advantage is that measurements can be acquired rapidly without the need for leak-subtraction protocols.

Gating capacitance (Cg) represents the amount of gating charge that moves ({Delta}Qg) in response to a small change in voltage ({Delta}V) and therefore reflects the slope of the Qg–V relationship ({Delta}Qg/{Delta}V). Cg is also dependent on the kinetics of charge movement and is therefore sensitive to the frequency (f) of the sinusoidal voltage command. Thus, capacitance measurements provide an assay of gating charge mobility reflecting both voltage-dependent and kinetic properties. When measured in response to a small amplitude, low frequency voltage perturbation Cg approximates the derivative of the steady-state Q–V relationship (Cg(V) = dQss/dV) (Fernandez et al. 1982 ; Fernandez et al. 1983 ; Taylor and Bezanilla 1979 ). Thus, if the Q–V can be described by a Boltzmann function Qss = [1 + exp(-ze(V - Vh)/kT)]-1 then Cg should exhibit a bell-shaped voltage dependence described by the derivative of a Boltzmann function. This relationship between Cg(V) and Qss(V) is strictly valid only when Cg is measured at a frequency approaching zero (Cg0). However, as discussed below, useful information about mSlo charge movement can be obtained using relatively large amplitude (30 mV) sinusoidal voltage commands at frequencies of hundreds or thousands of Hz.

Figure 1 A shows the Cg–V relationship for mSlo measured at 868 Hz in 0 Ca2+ (see Materials and Methods) from channels expressed in an excised macropatch. Cg exhibits a bell-shaped voltage dependence and is well fit by the derivative of a Boltzmann function (Figure 1 A, Fit). In nontransfected cells, the C–V relationship is flat (Figure 1 A, Control), representing only the uncompensated linear capacitance of the lipid bilayer and electrode (Co = Cb + Cs; see Equation 4). These contributions to the record in Figure 1 A were effectively eliminated by setting the baseline equal to zero at negative voltages (less than -100 mV) where C is voltage independent and presumed equal to Co. The Qg–V relationship was obtained by integrating the Cg–V trace with respect to voltage, and is plotted in Figure 1 B together with the normalized conductance–voltage (GK–V) relationship for mSlo in 0 Ca2+ (mean ± SEM, n = 23 [ Horrigan et al. 1999 ]). The Qg–V relationship is fit by a Boltzmann function (dashed line) characterized by a half-activation voltage (Vh) of 133 mV, corresponding to the peak voltage of the Cg–V curve, and an equivalent charge (z) of 0.59 e. Fits to Cg–V relationships obtained at 868 Hz from many experiments yielded values of Vh = 127.4 ± 3.4 mV and z = 0.61 ± 0.014 e (mean ± SEM, n = 15) in 0 Ca2+. These parameters are similar to those estimated in the preceding paper to characterize the charge and voltage dependence of mSlo voltage sensors (Vh(J) = 145 mV and zJ = 0.55 e [Horrigan et al., 1999]). The GK–V relationship is steeper and shifted to more positive voltages than the Qg–V, and can be approximated by a Boltzmann function (from the Qg–V fit) raised to the 4th power as shown in Figure 1 B.

The Cg–V relationship in Figure 1 A was measured during a 1-s voltage ramp from -160 to +200 mV. A sinusoidal command (f = 868 Hz) was superimposed on the ramp as indicated in Figure 1 C, and patch admittance (Yp) was measured for each cycle of the sinusoid using a phase-sensitive detector implemented in software (Herrington and Bookman 1995 ). Patch capacitance (Cp) was determined based on the expression:

(5)

where the real and imaginary terms represent the components of Yp appearing at phase angles of 0° and 90°, respectively, relative to the command voltage, and {omega} = 2{pi}f. Cg was then defined as the voltage-dependent component of Cp [Cg(V) = Cp(V) - Co]. The Cg–V relationship was unaffected by the polarity of the voltage ramp (data not shown), indicating that a pseudo steady-state condition was achieved at each voltage. Because the amount of gating charge detected was small (1–30 fC), admittance was typically measured using a relatively large amplitude 30-mV sinusoidal voltage command (60 mV peak to peak) to increase the signal to noise ratio. We were concerned that such a perturbation might alter the shape of the Cg–V relationship relative to that obtained with a small amplitude command. However, reduction of the sinwave amplitude from 30 to 3 mV had no detectable effect on the Cg–V relationship (Figure 1 D). This result suggests that the Cg–V was not distorted by the size of the sinusoidal command and is consistent with the weak voltage dependence of mSlo channel gating.

Gating Capacitance Represents mSlo Charge Movement
Although a voltage-dependent component of Cm was not detected in uninjected oocytes, it is important to verify that Cg arises from mSlo channels. High levels of heterologous expression of many membrane proteins in Xenopus oocytes have been shown to upregulate expression of endogenous ion channels (Tzounopoulos et al. 1995 ). It is conceivable that such endogenous channels could contribute to gating charge movement in cells expressing mSlo. However, several lines of evidence argue against such a contribution. First, similar Cg signals were observed using two different expression systems, Xenopus oocytes (Figure 1 D and 2 C) and HEK 293 cells (Figure 1 A and 2 A). Furthermore, Cg is Ca2+-sensitive and can be altered by mutating the mSlo channel.

The Ca2+ sensitivity of Cg is examined in Figure 2 A. Cg–V traces obtained in 0 or 60 µM Ca2+ from the same patch were normalized to peak capacitance and superimposed. The Cg–V relationships are similar in shape but shift to more negative voltages with increasing [Ca2+]i. The GK–V relationship for mSlo also exhibits a negative voltage shift upon application of Ca2+ in this concentration range (Cui et al. 1997 ).

The Ca2+ sensitivity of Cg suggests that this signal represents mSlo charge movement but does not rule out contributions from endogenous Ca2+-sensitive channels. To eliminate this possibility, we examined the properties of an mSlo mutant. Neutralization of a charged residue in the 17 domain of mSlo (R207Q) shifts the GK–V relationship to more negative voltages and reduces its slope relative to that of the wild-type (Diaz et al. 1998 ; Horrigan et al. 1999 ; Cui, J., and R.W. Aldrich, manuscript in preparation). We showed in the preceding paper that these shifts in the G–V can be accounted for by the allosteric voltage-gating scheme if the mutation allows voltage sensors to activate at more negative voltages without altering their charge. Consistent with this hypothesis, the Cg–V relationship for R207Q is approximately the same shape as that for mSlo but is shifted by almost -250 mV (Figure 2 B). This result also confirms that the Cg signal reflects the gating of mSlo.

The Kinetics of mSlo Gating Charge Movement
To assess the speed of charge movement, we examined the frequency dependence of Cg. In the simplest case, where gating charge movement can be represented by a two-state process, such as voltage-sensor activation from R to A, gating admittance

(6)

can be represented by an equivalent circuit consisting of a capacitor Cg0 in series with a resistor Rg0 where {tau}g = Cg0Rg0 is the time constant of gating charge relaxation at a particular voltage (Fernandez et al. 1982 ; Taylor and Bezanilla 1979 ). Where

(7)


(8)

When the frequency of the sinusoidal voltage command is low ({omega} << 1/{tau}), Yg reduces to a purely capacitive signal (Yg {cong} j{omega}Cg0), where Cg0 = dQss/dV. As {omega} increases, Cg should be attenuated, since the gating charge effectively cannot move fast enough to keep up with the voltage command. Cg({omega}) is a Lorenzian function that describes the frequency dependence of Cg. At higher frequencies, the gating current should also change phase with respect to the voltage command such that a component of Yg, described by the function Gg({omega}), appears in phase with the membrane voltage.

The frequency dependence of mSlo charge movement is shown in Figure 2C–F. Consistent with the above predictions, the Cg–V relationship is attenuated as the frequency of the sinusoidal voltage command is increased from 200 to 6,944 Hz (Figure 2 C). At the same time, a voltage-dependent signal appears in the orthogonal G trace and increases at higher frequencies (Figure 2 D). The DC current level during the voltage ramp is small and increases in a roughly linear manner with voltage, indicating a constant membrane resistance (Rm) of ~125 GV (Figure 2 E). Thus, the G signal in Figure 2 D represents a component of gating charge movement (Gg) and not a voltage-dependent change in membrane conductance. Gg is almost eliminated at 200 Hz, consistent with the prediction that Yg will reduce to a purely capacitive signal at low frequencies. Figure 2 F plots the amplitudes of Cg and Gg measured at +120 mV versus frequency for two experiments. Cg and Gg are well fit by Equation 7 and Equation 8, respectively, with a time constant ({tau}) of 70 µs. The relative amplitudes of the admittance components are also consistent with a 70-µs time constant, since Cg and Gg were normalized by Cg0 and Cg0/{tau}, respectively. Thus, a component of mSlo gating charge appears to move much faster than IK activation, which is described by a mean time constant of 1.63 ms at +120 mV (Horrigan et al. 1999 ). We will demonstrate below that an additional component of gating charge moves with the time course of channel activation but is too slow to be detected with admittance analysis.

Conclusions from Capacitance Measurements
Admittance analysis reveals several important properties of mSlo charge movement. Comparison of the Qg–V and GK–V relationships (Figure 1 B) suggests that charge movement can occur at voltages where most channels are closed. The frequency dependence of Cg shows that charge relaxes with a time constant that is much faster than that of IK activation (Figure 2 F). Together, these results suggest that admittance analysis detects charge movement associated with rapid closed-state transitions that precede channel opening. In terms of the allosteric voltage-gating scheme, such transitions result from voltage-sensor movement. That the Qg–V relationship can be fit by a Boltzmann function is consistent with the movement of each voltage sensor being described by a two-state model with a single transition between a resting (R) and an activated state (A). The simple voltage dependence of Qg also supports the notion that the voltage sensors, in different subunits of the mSlo homotetramer, behave identically and act independently. The approximate 4th power relationship between Qg–V and G–V is consistent with the assumption that channel opening is linked to the activation of four voltage sensors. However, as discussed below, this relationship may be affected by the different ionic conditions under which IK and Ig are measured.

Advantages and Limitations of Admittance Analysis
Our results show that admittance analysis provides a sensitive method for detecting and characterizing some aspects of mSlo gating charge movement. By using a large amplitude sinusoidal voltage command in combination with a voltage ramp, we were able to acquire the Cg–V relationship rapidly, and to determine Qg(V) at submillivolt intervals. The speed of mSlo charge movement is advantageous for admittance analysis because it allows measurements to be performed at hundreds or thousands of Hz where the signal to noise ratio is high (Gillis 1995 ; Lindau and Neher 1988 ; Lollike et al. 1995 ). By the same token, as discussed below, this technique is not well suited to detecting slow components of charge movement and may present difficulties in dissecting complex kinetic behavior.

The charge movement detected with capacitance measurements is much faster than IK activation. However, any scheme that assumes the C–O conformational change is voltage dependent or that channel opening affects the ability of voltage sensors to move requires that a component of gating charge will relax with the kinetics of IK activation. The frequency dependence of Cg (Figure 2 F) can be adequately fit by a single Lorenzian function between 200 and 7,000 Hz and therefore provides no evidence for a slow component of gating charge, which should appear as an additional Lorenzian component at low frequencies. However, the frequency range of our measurements may limit our ability to detect such components. For example, charge that moves with a time constant of 2 ms would produce a Cg component that is attenuated by ~85% at frequencies >200 Hz.

Admittance analysis is also not an ideal method for dissecting a model as complex as the one we have proposed for mSlo. The allosteric scheme predicts that multiple kinetic components of charge movement will result from C–C, C–O, and O–O transitions. Admittance analysis detects charge movement associated with perturbations about an equilibrium distribution of channel states, and will therefore contain contributions from all of these sources. Slow transitions associated with channel opening should contribute little to Cg at the frequencies used in our experiments. However, fast transitions among closed or open states (C–C, O–O) should be detected. At voltages less than +100 mV, most channels are closed in 0 Ca2+, and Cg will mainly reflect C–C transitions. However, at more positive voltages, Cg should represent a combination of open- and closed-state charge movement. For this reason, gating currents measured in response to step depolarizations provide a better method for isolating the various transitions predicted by the model.

A Fast Component of Gating Current
Figure 3 A shows Ig evoked in response to a 0.5-ms pulse to +160 mV from a holding potential of -80 mV in 0 Ca2+. The ON current decays rapidly with a time course that is well fit by an exponential function (dashed line) with a time constant of 59 µs, similar to that determined with admittance analysis at +120 mV (70 µs). The OFF current measured at -80 mV decays more quickly, with a time constant of 17 µs. A family of Ig evoked at different voltages (0 to +160 mV) in response to 1-ms pulses is shown in Figure 3 B. The QON–V and QOFF–V relationships obtained by integrating IgON and IgOFF are plotted in Figure 3 C (open symbols) together with the Qg–V relationship obtained from capacitance measurements at 868 Hz in the same patch (solid line). At all voltages, QON and QOFF are equivalent, as expected for gating charge. The gating current and capacitance measurements superimpose from 0 to +120 mV but diverge at +160 mV.

Similar results were obtained with brief voltage pulses and capacitance measurements because both methods mainly detect fast charge movement. Figure 3 D compares the time course of Ig evoked at +160 mV to the initial activation of IK measured at the same voltage from a different experiment. Ig decays, to a large extent, before IK begins to increase. After 1 ms, IK increases to 31% of its steady-state amplitude, representing only 7% of maximum Po. Thus, the channel does not achieve a steady state during a 1-ms pulse, and Ig should reflect little if any slow charge movement that might be associated with channel opening. An important difference between the gating current and capacitance measurements is that the initial decay of IgON represents charge moved when most channels are closed, while Cg is measured after Po has reached a steady state and therefore reflects the behavior of both open and closed channels. Thus, Ig measurements allow better isolation of closed-state transitions owing to the large kinetic difference between IgON and IK.

Fast Ig: Isolation and Voltage Dependence
According to the allosteric model, the initial decay of Ig represents activation of voltage sensors from a resting (R) to an activated state (A) while channels are closed (i.e., C–C transitions). The exponential decay of IgON is consistent with such a two-state model. Moreover, in Figure 3 D, IK achieves an exponential time course (dashed line) at a time (arrows) when the gating current has almost completely decayed. This correlation between Ig and the delay in IK activation is consistent with Ig reflecting closed-state transitions in the activation pathway. However, QON measured during a 1-ms pulse is not only an assay of closed-state charge movement, as some channels do open during this time (Figure 3 D). QON measurements can also be contaminated by outward leak currents that often are observed at voltages greater than +200 mV. To better characterize closed-state transitions, the fast component of ON charge was isolated by fitting an exponential function to the decay of Ig during the first 200 µs of the voltage pulse when most channels are closed. The area under the fit (Qfast), as indicated by the shaded region in Figure 3 D, was used as an estimate of closed-state charge movement (QC). The Qfast–V relation in Figure 3 C (filled symbols) is similar to the Qg–V relation. When data were acquired over a larger voltage range (Figure 3 E), Qg and Qfast diverge at more positive voltages (Figure 3 F). The Qfast–V relationship in Figure 3 F is well fit by a Boltzmann function (z = 0.57 e, Vh = 136 mV). The difference between Qg and Qfast is expected, as it occurs at voltages (>100 mV) where channels begin to open, and Qg therefore cannot be equivalent to QC.

Figure 4 A1 plots the normalized Qfast–V relationships for many experiments. The data were initially fit with Boltzmann functions where all parameters were allowed to vary, yielding a mean equivalent charge <z> = 0.59 ± 0.03 e (mean ± SEM, n = 10). The Q–Vs were then refit with z = <z> and normalized as shown in Figure 4 A1. Although the individual plots are reasonably fit using identical values of z, they are scattered in their position along the voltage axis, similar to the mSlo GK–V relationships (Horrigan et al. 1999 ). To compare the shapes of the Q–Vs, the individual records were aligned as shown in Figure 4 A2 (open symbols) by shifting them along the voltage axis by {Delta}V = <Vh> -Vh where Vh is the half-activation voltage of an individual Q–V and <Vh> is the mean (155 ± 6.5 mV, n = 10) determined from Figure 4 A1. These voltage-shifted plots were then used to determine the average Q–V (Figure 4 A2, filled symbols). A Boltzmann function with z = 0.59 e and Vh = 155 is superimposed on the data (solid line).

Voltage Dependence of Fast Ig Kinetics
To further characterize the properties of closed-channel charge movement, we examined the voltage dependence of fast Ig kinetics. Time constants of fast Ig relaxation ({tau}gFast) were determined from exponential fits to ON and OFF currents for the experiments in Figure 4 A1 and are plotted in Figure 4 B1. OFF currents, measured at voltages less than +40 mV, were evoked after very brief pulses (0.05–0.25 ms) to +160 or +200 mV and therefore should represent the relaxation of closed channels. {tau}gFast exhibits a bell-shaped voltage dependence, consistent with a two-state model of voltage-sensor activation where forward and backward rate constants are voltage dependent. {tau}gFast–V relationships from three experiments that covered a large voltage range are compared in Figure 4 C1. The individual plots are similar in shape but shifted relative to each other along both axes. The amplitude differences resemble those described previously for the delay in IK activation (Horrigan et al. 1999 ) and may reflect temperature variation between experiments conducted at room temperature. To better compare the shape of the {tau}gFast–Vs, the plots were first shifted along the voltage axis based on the Q–V shifts determined in Figure 4 A. The data were then normalized to the mean {tau}gFast determined over an interval around the peak of the {tau}gFast–V (59.0 ± 2.2 µs, n = 10, from +100 to +180 mV). The resulting records, corresponding to Figure 4, B1 and C1, are plotted in Figure 4, B2 and C2, respectively, and exhibit improved alignment of the {tau}gFast–V relationships.

The data in Figure 4 B2 were fit with a function {tau}gFast = 1/({alpha}+ß), representing the predicted {tau}gFast–V relationship for a two-state process where the forward ({alpha}) and backward (ß) rate constants are exponential functions of voltage [{alpha}={alpha}0e()]. Fits were constrained such that the equilibrium constant J = {alpha}/ß equals one at the half-activation voltage of the Qfast–V (Vh(J) = 155 mV). The solid line in Figure 4 B2 represents the best fit and is characterized by a total equivalent charge of zJ = 0.51e (z{alpha} = +0.30 e, zß = -0.21 e).

Estimates of the charge associated with voltage-sensor activation (zJ) based on fits to the Qfast–V and {tau}gFast–V relationships (0.59 and 0.51 e, respectively) apparently differ. However, both relationships can be reasonably fit using the average of these two estimates (0.55 e) (Figure 4, A2, B2, and C2; dashed lines). This value of zJ was also used in the preceding article to reproduce the ionic current data using the allosteric voltage-gating scheme (Horrigan et al. 1999 ). One difference is that the value of Vh(J) used to fit the Qfast–V (155 mV) is 10 mV greater than that previously used to fit IK. In addition, the values of z{alpha} and zß used to fit the {tau}gFast–V relationship (z{alpha} = +0.33 e, zß = -0.22 e) indicate that the R–A transition in the allosteric model is not symmetrically voltage dependent as previously assumed.

A Slow Component of ON Charge Movement
Although the ON currents in Figure 3 appear to decay with a single-exponential time course, there is a significant slow component of charge movement. Figure 5 A plots a family of Ig evoked at +140 mV in response to voltage pulses of different duration (see also Figure 6 A). The peak amplitude of IgOFF increases rapidly with pulse duration, paralleling the rapid decay of IgON, and then remains relatively constant for pulses longer than 0.5 ms. The total gating charge moved during the pulse (Qp) was determined by integrating IgOFF and is plotted versus pulse duration in Figure 5 B. Qp increases with a time course that can be fit by a double-exponential function (solid line) with a fast phase (QpFast) corresponding to the rapid decay of IgON, and an additional phase that is roughly 100-fold slower. The slow component (QpSlow) relaxes with a time constant ({tau}gSlow) of 4.22 ms and represents a significant fraction of the total gating charge movement at +140 mV (43%) but is too slow to be observed as a component of IgON. This point is illustrated in Figure 5 C, which compares IgON evoked at +140 mV to Qp'(t) (dashed line). Qp'(t) is the time derivative of the double-exponential fit to Qp(t) and should represent the time course of IgON (Qp'(t) = dQON/dt = IgON). These two relationships superimpose, demonstrating that observed IgON kinetics are consistent with the presence of a large slow component of ON charge movement.

The predicted amplitude of the slow component of IgON, determined from Qp'(t), is small (2.1 pA) because it decays slowly. For similar reasons, the slow component of ON charge could not be reliably measured from IgON. Small sustained outward currents on the order of a few pA were often observed at high voltages, presumably representing residual ionic or leak current. For example, the current trace in Figure 5 C decays to a steady-state level of 2.2 pA at the end of the pulse. While such small currents have little effect on measurement of Qfast they can contaminate estimates of slow charge determined by integrating IgON over a 20-ms pulse. The slow component of Qp from Figure 5 B is only 8 fC, equivalent to a 0.4 pA current for 20 ms. Measurements of OFF charge (Qp) provide a more reliable estimate of slow charge movement because leak is constant at the holding potential.

The voltage dependence of Qp(t) is examined in Figure 6. Families of Ig evoked at different voltages in response to pulses of 0.06–20 ms duration are shown in Figure 6 A. At each voltage, Qp was plotted versus pulse duration (Figure 6 B) and fit with a double-exponential function as in Figure 5 B. The plots represent data from three experiments and were normalized to the total fast charge movement QTfast estimated from the amplitude of a Boltzmann fit to the Qfast–V relationship for each experiment. The indicated voltages were corrected for shifts in the Qfast–V relationship as determined in Figure 4 A.

A slow component of Qp is observed in Figure 6 B for V >= +100 mV. The time constant of QpSlow ({tau}gSlow) is comparable to that for IK activation ({tau}(IK)) measured from +140 to +240 mV (Figure 6 C). The similar magnitude and voltage dependence of {tau}gSlow and {tau}(IK) suggest that slow charge movement is limited by channel opening. These kinetics also show that gating charge and open probability equilibrate on a similar time scale. Therefore, QOFF determined with a 1-ms voltage pulse, as in Figure 3 C, can underestimate steady-state QOFF (Qss), determined with a 20-ms pulse, by as much as 40%. Despite this difference, the Qss–V and Qfast–V relationships are similar in shape. Figure 6 D compares normalized Qss–Vs from four experiments to the normalized Qfast–V and GK–V relationships. Qss–V almost superimposes with Qfast–V, and the steady-state data were fit with Boltzmann functions with an equivalent charge z = 0.65 ± 0.03 e (mean ± SEM, n = 4), indicating a slightly steeper voltage dependence than Qfast.

The Relationship between Slow Charge Movement and Channel Activation
The predominantly exponential time course of mSlo IK suggests that the kinetics of voltage-dependent activation are dominated by a rate-limiting transition (Horrigan et al. 1999 ). The similar kinetics of QpSlow and IK relaxation implies that slow gating charge movement also reflects this rate-limiting step. It is important to distinguish between two possible sources of slow charge movement. First, the rate-limiting step may represent a voltage-dependent conformational change and therefore contribute directly to QpSlow. Second, the rate-limiting step may contribute indirectly to QpSlow by limiting the speed of other voltage-dependent transitions in the activation pathway. The data suggest that both of these mechanisms contribute to slow charge movement in mSlo.

We have previously concluded that the transition from a closed to open conformation represents the rate-limiting step in mSlo activation and is weakly voltage dependent (Horrigan et al. 1999 ). Hence, the rate-limiting step should contribute directly to slow charge movement. However, the charge associated with the C–O transition (zL = 0.4 e) was estimated to represent only 15% of the total charge per channel. In contrast, slow charge movement in mSlo can represent >40% of the total ON charge (Figure 6 B). These results are inconsistent with the idea that QpSlow merely represents the charge moved during the C–O transition, but they can be understood in terms of the allosteric voltage-gating scheme (Figure 14).

The allosteric model predicts that the majority of charge movement can be attributed to voltage-sensor activation. Fast Ig is evoked in response to a voltage step as sensors initially equilibrate between resting (R) and activated (A) conformations while the channel is closed. Qfast is determined by the voltage-dependent equilibrium constant (J) that characterizes the R–A transition. In addition, a slow component of charge movement should be produced as channels open, representing the C–O transition. However, voltage-sensor movement can also contribute to QpSlow. When a channel opens, the equilibrium constant for voltage-sensor activation increases by the allosteric factor D, causing sensors to reequilibrate between R and A and produce additional charge movement. This charge movement will be slow because the voltage-sensor reequilibration is limited by the speed of channel opening.

The amplitude of QpSlow should depend upon the number of channels that open as well as the fraction of voltage sensors that are initially activated before channels open. For example, at very positive voltages (approximately +300 mV), the model predicts that voltage sensors can be completely activated with channels closed. In this case, channel opening cannot cause additional voltage sensors to be activated so QpSlow will represent only the charge associated with the C–O transition (zL). At less positive voltages, however, QpSlow will represent a combination of channel opening and voltage-sensor reequilibration and may therefore be larger than zL. We will demonstrate later that the magnitude and voltage dependence of QpSlow are consistent with the allosteric gating scheme (Figure 14). The notion that the C–O transition limits slow charge movement is also important in understanding the properties of IgOFF as discussed below.

Three Components of OFF Gating Charge Movement
The large slow component of Qp(t) observed at V >= +140 mV in Figure 6 B indicates that QOFF increases with pulse duration. In contrast, the peak amplitude of IgOFF remains roughly constant or decreases with pulse duration at the same voltages (Figure 6 A). That IgOFF can decrease or remain constant while its integral (QOFF) increases implies that the kinetics of OFF current change with pulse duration. This change is obvious in Figure 7 A, which compares OFF currents evoked at -100 mV after pulses to +140 mV of different duration (0.06–20 ms). Two components of IgOFF are evident from these records. After brief pulses (0.06 or 0.11 ms), OFF current decays with a rapid exponential time course, but an additional slower component appears as pulse duration is increased. The decay of IgOFF at all pulse durations can be well fit by double-exponential functions with time constants of 15.5 and 59 µs (Figure 7 B). Both components decay within 300 µs and therefore appear to be fast relative to the time course of channel closing. Potassium tail currents decay with a time constant of 172 ± 15 µs at -80 mV (Horrigan et al. 1999 ) and therefore require approximately 5{tau}(IK) = 900 µs to decay completely. However, a slower component of OFF charge movement can be detected by plotting the integral of IgOFF (QOFF(t); Figure 7 C). QOFF(t) measured after a brief (0.06 ms) voltage pulse achieves a steady state within 300 µs (Figure 7 C, arrow), consistent with the rapid decay of IgOFF. In contrast, QOFF(t) measured after a 20-ms pulse requires >1 ms to reach a steady state, indicating a slow component of charge relaxation. This component of QOFF is not evident in the corresponding IgOFF trace because it is slow and represents <20% of the total OFF charge.

The components of QOFF(t) relaxation were further analyzed by plotting the quantity (QOFF(t) - QOFFss) where QOFFss is the steady-state value of QOFF(t) measured 3 ms after the voltage pulse (Figure 7 D). The relaxation of (QOFF(t)- QOFFss) after a brief pulse (0.06 ms) can be fit by a single-exponential function as indicated by a linear relationship on this semilog plot ({tau}F = 15.5 µs). The relaxation of (QOFF(t) - QOFFss) after a prolonged pulse (average of 10–20-ms records) is more complicated and was best fit by three exponential components ({tau}F = 15.5 µs, {tau}M = 59 µs, {tau}S = 448 µs), indicated by dashed lines in Figure 7 D, where {tau}F was constrained to that used to fit the 0.06-ms record. On average, time constants of 15.7 ± 1.3, 64.7 ± 10.6, and 580 ± 50 µs were measured at -80 mV (mean ± SEM, n = 6).

The time course of development of the OFF charge components were examined by fitting (QOFF(t)- QOFFss) with triple-exponential functions for all pulse durations (Figure 7 E). The time constants (termed Fast, Medium, and Slow) were determined from the 0.06- and 10–20-ms traces as in Figure 7 D, and component amplitudes were varied to fit the other records. The QOFF component amplitudes (QOFFfast, QOFFmed, and QOFFslow) are plotted versus pulse duration in Figure 7 F. The Fast component develops rapidly and then slowly decreases in amplitude as pulse duration is increased. At the same time, a parallel increase in the Medium and Slow components is observed. The slow relaxations in the development of all three components were fit by exponential functions (solid lines) with a time constant of 4.2 ms. This time constant is identical to that used to fit QpSlow (Figure 5 B) and is therefore assumed to represent the time course of channel opening. As discussed below, the results in Figure 7 F suggest that the Fast component of OFF charge movement represents the relaxation of closed channels, while the Medium and Slow components represent the relaxation of open channels.

Factors Influencing OFF Component Characterization
Accurate separation of QOFF components depends on several factors, including the estimation of their time constants. {tau}F is most easily determined because the fast component is large and can be examined in isolation using brief voltage pulses. The Slow component can also be effectively isolated because {tau}S is almost 10-fold larger than {tau}M. However, the small amplitude of the Slow component and its sensitivity to baseline drift make {tau}S more difficult to determine than {tau}F. The relaxation of QOFF(t) to a steady state in Figure 7 C indicates that IgOFF decays to the baseline level after ~1 ms. A small offset or drift in baseline current can prevent QOFF(t) from achieving such a steady state and affects determination of {tau}S and QOFFslow. To minimize such artifacts, the Ig baseline was typically set equal to the mean current measured during an interval 4–5 ms after the end of the pulse. Despite this precaution, drift in QOFFss was observed in some experiments (data not shown) and contributes to variability in the estimate of {tau}S.

The medium time constant ({tau}M) was also difficult to determine because it is only fourfold slower than {tau}F and cannot be studied under conditions where the Fast and Slow components are absent. Thus, estimates of {tau}M from triple exponential fits to QOFF relaxation were sensitive to the estimates of {tau}F and {tau}S. Error bars in Figure 7 F indicate the effect of ±10% changes in {tau}M on the estimated amplitudes of the different OFF components (with {tau}F and {tau}S held constant). Such variation still allows reasonable fits to QOFF(t) (data not shown); however, an increase in {tau}M results in a decrease in the measured QOFFmed and a complimentary increase in QOFFfast. Larger changes in {tau}M produce inadequate fits to QOFF(t), and the time course of the Medium component development becomes biphasic as the separation of Fast and Medium components is compromised.

Measurements of QOFFmed can be affected by baseline drift or variation in {tau}S. Therefore, the development of Fast and Medium components were also studied by fitting IgOFF with double-exponential functions (Figure 7 B), a procedure that is less sensitive to the slow component. Figure 7 G plots the amplitude of the IgOFF components versus pulse duration, indicating a time course of Fast and Medium component development similar to that determined from QOFF (Figure 7 F). Exponential fits in Figure 7F and Figure G, used identical values of {tau}gSlow. However, in experiments where baseline drift was a problem, fits to IgOFF produced more consistent results and were used to determine {tau}gSlow.

mSlo Charge Movement and Allosteric Voltage Gating
The presence of three components of OFF gating charge movement, their kinetics, and development with pulse duration can be understood in terms of the allosteric voltage-gating scheme (Figure 14). As indicated in Figure 7 H, the allosteric model predicts that OFF charge relaxation will be characterized by Fast, Medium, and Slow components that reflect C–C, O–O, and O–C transitions, respectively. When mSlo channels are closed, OFF currents should represent the relaxation of voltage sensors from an activated to a resting state, corresponding to C–C transitions in the gating scheme. Since brief voltage pulses allow few channels to open, the fast relaxation of IgOFF after such a pulse ({tau}F) mainly reflects the kinetics of this closed-state relaxation pathway. As pulse duration is increased, channels open and their deactivation after the pulse reflects a more complex relaxation pathway involving O–O and O–C transitions. The model predicts that voltage sensors can move even when channels are open. Therefore, the OFF current should exhibit a component that reflects relaxation of voltage sensors from an activated to a resting state, corresponding to O–O transitions in the gating scheme. If these open-state transitions account for the Medium QOFF component, to account for the difference between {tau}M and {tau}F, we must assume that voltage-sensor relaxation is slower when the channel is open than when it is closed. This is a reasonable assumption because the allosteric mechanism requires that channel opening increase the equilibrium constant for voltage-sensor movement D-fold, stabilizing the activated state (A) relative to the resting state (R). Finally, the model predicts that there will be a slow component of OFF charge movement associated with the transition of open channels back to the closed state. Therefore, the slow component should have the same time course as channel deactivation. We will argue later that differences in the observed time course of IK deactivation and slow charge movement ({tau}S) may reflect effects of ionic conditions on channel gating.

If the Fast component of OFF charge movement represents the relaxation of closed channels while Medium and Slow components represent the relaxation of open channels, the effect of pulse duration on the relative amplitude of these components can be understood in terms of the kinetics of channel activation. QOFFfast increases initially because voltage sensors can be activated rapidly during brief pulses while channels are closed. As pulse duration increases, the number of closed channels is reduced and QOFFfast decreases with the time course of channel activation. At the same time, both Medium and Slow components increase, reflecting an increase in the number of open channels.

Charge Movement Measurements Are Not Contaminated by Ionic Currents
An important conclusion from the above analysis is that the slow components of ON and OFF charge movement are limited by channel opening and closing. Since the kinetics of these components are similar to those of IK, it is critical to establish that they do not represent contamination of Ig by residual ionic currents. The slow component of ON charge movement was detected as an increase in QOFF measured after pulses of different duration, whereas the slow OFF charge was seen as a component of QOFF relaxation. Thus, the presence of an inward potassium tail current could potentially contribute to both measurements. This possibility appears unlikely because gating current records that give rise to large slow components of ON charge movement (e.g., Figure 5 and Figure 6) do not exhibit appreciable sustained (ionic) current during the voltage pulse. In addition, the slow increase in Qp with pulse duration involves simultaneous changes in the amplitudes of all three components of QOFF relaxation. The Fast component decreases while the Medium and Slow components increase (Figure 7 F). Although the Slow component relaxes with kinetics similar to that of ionic tail currents, it accounts for only a small fraction of QpSlow. Finally, as discussed below, the relative amplitudes and voltage dependence of the different QOFF components are consistent with previous estimates of the charge and equilibrium properties of C–C, O–O, and C–O transitions in the allosteric scheme (Figure 14).

Testing the Allosteric Voltage-gating Scheme
The relationship between Po and QOFFfast. The allosteric model predicts that the fast component of OFF charge movement should be eliminated after voltage pulses that open all channels. One way to increase Po is by stepping to more positive voltages. Figure 8 A plots the time course of QOFF component development at +240 mV. The decay of QOFFfast is more complete than at +140 mV (Figure 7 F), consistent with a voltage-dependent increase in Po. It is likely that the fast component was not eliminated because, in the absence of Ca2+, mSlo channels are maximally activated only at very positive voltages (greater than +300 mV) (Horrigan et al. 1999 ). However, in the presence of 60 µM Ca2+, channels can be fully activated at +160 mV. Figure 8 B compares the relaxation of QOFFQOFFss after a 0.1- or 20-ms pulse under these conditions. The 0.1-ms trace decays rapidly and is fit by a triple exponential function ({tau}F = 23.8 µs, {tau}M = 150 µs, {tau}S = 822 µs), with the Fast component representing the majority of OFF charge (91%). However, the 20-ms record is well fit by a double-exponential function using only {tau}M and {tau}S. This confirms that the Fast component can be eliminated and that the relaxation of open channels back to the closed state contributes only to the Medium and Slow components of QOFF.

Voltage dependence of QOFF component amplitudes. To further test the allosteric model, we examined the effect of repolarization voltage on the relative amplitudes of QOFF components. Ig was evoked in response to pulses of different duration to +160 mV (0.1–20 ms). After each pulse, the membrane was repolarized to either -80 or 0 mV, and OFF currents were analyzed as in Figure 7. The amplitudes of the three QOFF components are plotted versus pulse duration in Figure 8C and Figure D, for -80 and 0 mV, respectively. The component amplitudes were normalized to Qfast measured in response to a pulse from -80 to +160 mV because a 20% increase in this quantity was observed during the course of the experiment. In the absence of Ca2+, steady-state open probability at 0 mV is expected to be small (<10-4) (Horrigan et al. 1999 ). Therefore, both repolarization voltages should be sufficiently negative to close most channels. The Qss–V relationship (Figure 6 D) indicates that there is also little change in the steady-state charge distribution between -80 and 0 mV, so total QOFF is similar at -80 or 0 mV. Figure 8C and Figure D, shows that the time course of total OFF charge development (Qp(t), open symbols) is also unaffected by repolarization voltage. This is expected, since Qp(t) represents the time course of ON charge movement and should depend only on the voltage during the pulse (+160 mV). Similarly, the development time course of the three QOFF components and the amplitude of the fast component are unaffected by repolarization voltage. However, a change in the relative amplitudes of the Medium and Slow components is observed. QOFFslow increased 2.4-fold at 0 mV while QOFFmed decreased, such that total QOFF remained constant. This complementary change in QOFFmed and QOFFslow supports the idea that both represent charge movement and are not contaminated by ionic currents.

The effect of voltage on the relative amplitude of Slow and Medium components of QOFF can be understood in terms of the allosteric gating scheme (Figure 14). According to the model, the Medium component represents open state (O–O) transitions while the Slow component is limited by channel closing (O–C). Therefore, QOFFmed reflects the voltage-dependent reequilibration of channels among open states. If the membrane is repolarized to a sufficiently negative voltage, QOFFmed will be maximal because open channels will rapidly occupy the leftmost open state (O0) before closing. Under these conditions QOFFslow will be small, representing only the charge moved during the transition from O0 to C0 (zL). However, if the membrane is repolarized to a less negative voltage, the open-channel equilibrium may favor occupancy of intermediate open states (Oi) rather than O0, and QOFFmed will be reduced. At the same time, QOFFslow will increase to reflect relaxation from Oi to the resting closed state (C0).

To examine the quantitative predictions of the allosteric scheme (Figure 14), it is convenient to compare the charge distributions predicted for Closed and Open channels (Figure 8 E, QC(V) and QO(V)). QC can be expressed in terms of the voltage-sensor equilibrium constant J(V) and charge zJ.


(9)

Therefore, QC(V) has the same shape as the Qfast–V relation, with a maximum amplitude of 2.2 e (4 zJ) when zJ = 0.55 e. QO(V) is determined by the open-channel voltage-sensor equilibrium constant (DJ), the voltage-sensor charge zJ, and the charge for the C–O transition (zL = 0.4 e):

(10)

When D is assigned a value of 17, as in the preceding paper, the half-activation voltage for QO(V) is shifted by -130 mV relative to that of QC(V), indicating that voltage sensors are easier to activate when channels are open ({Delta}{Delta}G2.83 kT). The relative amplitudes of QOFFmed and QOFFslow predicted by the model are indicated by arrows in Figure 8 E at repolarization voltages of -80 and 0 mV. If voltage sensors equilibrate before channels close, the Medium OFF component evoked from an open channel can be expressed in terms of QO:

(11)

where VP is the pulse voltage and VR is the repolarization voltage. The Slow OFF component is determined by the difference of QO and QC.


(12)

As illustrated in Figure 8 E, the model predicts that QOFFslow will increase 1.93-fold when OFF charge is measured at 0 mV rather than -80 mV, similar to the 2.38-fold change observed in Figure 8C and Figure D.

Simulations of the allosteric model. The results discussed thus far are qualitatively consistent with the behavior of the allosteric gating scheme (Figure 14). Simulations based on the model as shown in Figure 9, Figure 10, and Figure 11 also reproduce the major features of the data. However, the parameters that were ultimately used to fit Ig differ from those used to describe ionic currents (Horrigan et al. 1999 ). Some of these differences are small and may simply reflect a greater accuracy in characterizing fast voltage-sensor movement with gating currents. Other differences, relating to the slow charge movement, suggest that ionic conditions alter mSlo channel gating.

In the preceding article, parameters for the allosteric scheme (Figure 14) were estimated based on fits to the GK–V relationship, and the voltage dependence of both IK relaxation kinetics ({tau}(IK)-V) and the delay in IK activation ({Delta}t(IK)-V). The charge assigned to voltage-sensor movement (zJ = 0.55 e) was identical to that used here to fit the Qfast–V (Figure 4 A2) and {tau}gFast–V (Figure 4 B2) relationships. However, the half-activation voltage of the Qfast–V (Vh(J) = 155 mV) determined from gating current measurements is 10 mV more positive than previously estimated. Although this discrepancy is small, it is useful to consider several factors that could potentially contribute to such a difference. First, patch to patch variability is observed for both ionic and gating current data in the position of relationships such as the GK–V and Qfast–V along the voltage axis (Figure 4; see also Horrigan et al. 1999 ; Stefani et al. 1997 ). We have attempted to minimize the effects of such shifts by averaging results from many experiments. Nonetheless, such variation could contribute to differences between ionic and gating current data. In addition, the estimate of Vh(J) based on IK recordings is less direct and therefore likely to be less accurate than that determined from gating currents. The previous estimate of Vh(J) was based, in part, on the ability of the allosteric scheme (Figure 14) to fit the {Delta}t(IK)–V relationship. Parameters were assigned to the model with the simplifying assumption that the rate constants for voltage-sensor movement ({alpha}, ß) are symmetrically voltage dependent (z{alpha} = -zß). Under this condition, with Vh(J) = 145 mV, the model reproduces the observation that the maximum delay is observed at approximately +153 mV (Vmax({Delta}t)). However, the {tau}gFast–V relationship (Figure 4 B2) indicates that z{alpha} (0.33 e) is greater than -zß (0.22 e). Under this condition the predicted relationship between Q–V and {Delta}t(IK)–V changes such that Vh(J) > Vmax({Delta}t). Thus, Vh(J) is not merely determined by the {Delta}t(IK)–V relationship but is also influenced by z{alpha} and zß. Finally, experimental conditions were different for IK and Ig measurements and might contribute to a real difference in channel gating. For example, {Delta}t(IK) was measured at a lower temperature (5°C) than Ig (20–22°C). In Shaker K+ channels, decreased temperature has been shown to shift the Q–V relationship to more negative voltages (Rodriguez et al. 1998 ), consistent with the difference in Vh(J) estimated for mSlo from IK and Ig data. In addition, IK was recorded in symmetrical 110 mM K+ while Ig was recorded with NMDG and TEA replacing internal and external K+, respectively. Permeant and blocking ions are known to alter the gating of many K channels (Armstrong 1971 ; Chen et al. 1997 ; Fedida et al. 1999 ; Matteson and Swenson 1986 ; Sala and Matteson 1991 ; Swenson and Armstrong 1981 ; Wang et al. 1999 ; Yeh and Armstrong 1978 ), including BK channels (Demo and Yellen 1992 ; Miller et al. 1987 ; Neyton and Pelleschi 1991 ) (see DISCUSSION).

Initial Ig simulations (Figure 9) were generated using parameters determined from a combination of ionic and gating current measurements. The parameters describing the R–A transition for closed channels (zJ = 0.55, Vh(J) = 155 mV, z{alpha} = 0.33, zß = -0.22) were determined from Qfast–V and {tau}gFast–V relationships as described above. The R–A equilibrium for open channels was assumed to differ from that for closed channels by the allosteric factor D = 17, estimated in the preceding article. The rate constants for this transition were assumed to be symmetrically affected by channel opening such that the forward rate is increased f-fold (f=) and the backward rate is decreased by the same factor. Rate constants for the C–O transitions were identical to those used to fit the IK data at 20°C (Horrigan et al. 1999 ). Finally, simulated Ig was scaled to experimental records by estimating the number of channels (N) based on the expression N = QTfast/4zJ.

Figure 9 A plots a family of IgON evoked at different voltages and compares them to simulations of the allosteric scheme (solid lines). The model reproduces the fast decay and relative amplitudes of these ON currents. The amplitudes of fast gating currents are sensitive to filtering; therefore, the voltage command used in the simulation and the resulting current were filtered at 20 kHz to reproduce experimental conditions (see Materials and Methods). Figure 9 B plots a family of gating currents evoked at +140 mV in response to pulses of different duration (from Figure 5 A). The model (solid lines) reproduces the time course and relative amplitudes of ON and OFF currents in response to brief pulses.

The time constants of Fast and Medium charge movement ({tau}F and {tau}M), predicted by the model, are plotted in Figure 9 C (solid lines). The {tau}F–V relationship is defined ({tau}F = [{alpha} + ß]-1) by the parameters assigned to the R–A transition when the channel is closed, and is identical to the fit of the {tau}gFast data in Figure 4, B2 and C2 (dashed lines). {tau}gFast measured from simulated currents (Figure 9 C, solid symbols) superimposes on {tau}F, confirming that exponential fits to fast Ig can be used to estimate the properties of closed-channel voltage-sensor movement. Similarly, the Qfast–V relationship, determined from these fits, superimposes on the QC–V relationship (Figure 9 D) defined by the model (Equation 9).

The {tau}M–V relationship predicted by the model is the same shape as the {tau}F –V but is shifted to more negative voltages owing to the allosteric interaction between channel opening and voltage-sensor movement ({tau}M = ({alpha}f + ßf/D)-1). Measurements of {tau}M from several experiments (Figure 9 C, open symbols) are similar to those predicted by the model, consistent with the assumption that the forward and backward rate constants for voltage-sensor activation are symmetrically affected by channel opening (i.e., f==4.13). A better fit to the data is obtained if f is increased to 4.8 (dashed line) but, given the limited number and voltage range of {tau}M measurements, we continue to assume f = 4.13 in the following simulations. The similar voltage dependence of the {tau}M and {tau}F data supports the conclusion that both Fast and Medium components of OFF charge represent voltage-sensor movement.

In addition to reproducing Ig in response to brief pulses, the model exhibits a slowing of IgOFF with increased pulse duration (Figure 9 B). However, this effect is more prominent in the data, suggesting that the model underestimates the amount of slow charge movement. To examine the time course and magnitude of slow charge predicted by the model, IgON was simulated in response to 20-ms pulses to different voltages and then integrated to obtain QON(t) (Figure 9 E, solid lines). The time course of QON is biphasic and the fast component matches the data (Qp(t); Figure 9 E, symbols), but the slow component is too small, especially at lower voltages. One possible explanation for this underestimate of QpSlow is that the model underestimates the number of channels that open at different voltages. In other words, the shape of the Po–V relationship and/or its position along the voltage axis may not be accurately reproduced. Since the GK–V relationship was well fit by the allosteric scheme in the preceding paper, this situation could occur if channel opening is enhanced under the ionic conditions where Ig is measured. To test this possibility, we further analyzed the voltage dependence and kinetics of the different charge movement components.

Estimating Po from charge movement. The allosteric model predicts that slow changes in both ON and OFF charge movement components are related to channel opening and closing. Therefore, the amplitudes of these components are related to open probability as well the charge distribution for open (QO) and closed (QC) channels. For example, the fast component of OFF charge depends on QC and the number of closed channels at the end of a voltage pulse (1 - Po):

(13)

For a particular pulse voltage (VP) and holding potential (HP), the second term in this expression can be determined by measuring the fast component of ON charge:

(14)

Therefore, Po can be estimated by comparing fast components of ON and OFF charge:

(15)

Figure 10 A plots the steady-state Po–V relationship estimated in this way for three experiments where QOFFfast(VP) was measured after a 20-ms pulse and QONfast(VP) was determined from an exponential fit to IgON (i.e., Qfast). Although measurements are scattered, reflecting, in part, the previously noted difficulties in separating QOFF components, the data generally follow the shape of the Po–V relationship predicted by the original model parameters (Figure 10 A, Case A) but are shifted to more negative voltages. Two additional Po–V relationships (Cases B and C) are superimposed on the data and will be used throughout this analysis. Case B indicates the prediction of the allosteric scheme (Figure 14) when the equilibrium constant L is increased 12-fold (equivalent to {Delta}{Delta}G = 2.5 kT) while leaving the other parameters unchanged. The Po–V relationship indicated by Case C is roughly the same shape as Case A but is shifted along the voltage axis. Case C was not generated by a gating scheme but can be used in combination with the QC and QO relationships defined by the original model to make predictions about the voltage dependence of different charge movement components. As discussed below, various aspects of the data are consistent with these altered Po–V relationships.

The predicted amplitude of the slow component of OFF charge is directly proportional to Po:

(16)

Figure 10 B plots normalized QOFFslow versus voltage for the same experiments as in Figure 10 A. Again, the data follow the general shape of the Po–V relationship predicted by Case A but appear shifted to more negative voltages. The model relationships were generated from the above expression for QOFFslow where Po was specified by Case A, B, or C in Figure 10 A and QO and QC were determined from the parameters assigned to the original model as illustrated in Figure 8 E. The data and model traces were normalized to the total fast charge movement QTfast for each experiment. According to the model QTfast = 4zJ, therefore, the maximum amplitude of the normalized data should be [QO(HP) - QC(HP)]/4zJ. That the data fall within the amplitude range predicted by the model is therefore consistent with the relative amplitudes of QO(HP), QC(HP), and zJ specified in the model.

The Medium component of OFF charge is larger and therefore easier to measure than QOFFslow but its voltage dependence is determined by QO(V) as well as Po(V):

(17)

Figure 10 C compares the normalized QOFFmed–V relationships to the model predictions. Again, the data plots are similar in shape and magnitude to the prediction of Case A but are shifted to more negative voltages. Both data and model predictions were normalized to QTfast as in Figure 10 B such that the maximum amplitude should be [QO(VP) - QO(HP)]/4zJ. Therefore, the magnitude of QOFFmed is consistent with QO(V) and zJ specified in the model.

The amplitude of the data and model predictions in Figure 10A, Figure B, and Figure C, as noted above, are influenced by several factors in addition to Po. These include model parameters (QO(V), QC(V), zJ, zL) as well as our ability to separate QOFF components and determine QTfast. To better examine the voltage dependence of the data, IgOFFmed–V relationships from several experiments were normalized together with the model predictions to a maximum amplitude of one (Figure 10 D). IgOFFmed is proportional to QOFFmed, so the model relationships represent normalized versions of those used in Figure 10 C. IgOFFmed was normalized based on a Boltzmann fit to the IgOFFmed–V relationship for each experiment (z = 0.98 e). When scaled in this way, the data from different experiments superimpose. Case C represents a Boltzmann fit to these normalized data (z = 0.98 e, Vh = 121 mV). The Po–V relationship for Case C (Figure 10 A) was derived from this fit and the expression for QOFFmed (Equation 17).

The slow component of ON charge movement (QpSlow) should exhibit a complex voltage dependence that is determined by Po(V), QC(V), and QO(V):

(18)

The QpSlow–V relationships plotted in Figure 10 E were normalized by QTfast and exhibit amplitudes that are larger than predicted by Case A, but are similar to those specified by Cases B and C. The model predicts that QpSlow will have a bell-shaped voltage dependence and that QpSlow/QTfast approaches a limiting value of zL/4zJ at positive voltages (Figure 10 E, dashed line). Our measurements do not extend to high enough voltages to test these predictions. However, the data fall close to the relationships determined by Cases B and C over the voltage range tested, and appear to trend downward at the highest voltages. Importantly, the comparison of QpSlow–V relationships for Cases A, B, and C demonstrate that the amount of slow charge movement is highly sensitive to Po and that QpSlow can be considerably larger than the charge associated with the C–O transition.

Finally, we examined the ability of the model to reproduce slow charge movement kinetics. Figure 10 F plots {tau}gSlow over a large voltage range. At positive voltages, {tau}gSlow was measured from the time course of development of the Medium component of IgOFF as in Figure 7 G. At negative voltages, {tau}gSlow was determined from the relaxation of Slow QOFF(t) ({tau}S). The dashed line in Figure 10 F represents a fit of the allosteric scheme (Figure 14) to the time constants of IK relaxation ({tau}(IK)), measured in the preceding paper (equivalent to Case A). {tau}gSlow is faster than {tau}(IK) for V > ~+100 mV and is slower than {tau}(IK) at negative voltages. However, the voltage dependence of {tau}gSlow can be fit (Figure 10 F, solid line) by adjusting the model parameters as specified for Case B where the equilibrium constant L is increased 12-fold. Both {tau}gSlow and {tau}(IK) are weakly voltage dependent from -80 to 0 mV, consistent with the idea that the slow relaxation of QOFF is limited by channel closing.

Taken together, the data in Figure 10 support the hypothesis that the properties of slow charge movement can be accounted for by the allosteric voltage-gating scheme (Figure 14), provided we assume that Po is increased under the conditions where gating currents are measured. Coordinated changes in all three components of QOFF are observed with pulse voltage, consistent with the assumption that their amplitudes depend upon the Po–V relationship. The relative amplitudes of these components are also consistent with their proposed source in terms of the allosteric scheme and with the charges assigned to various transitions in the model. The voltage dependence of the Medium OFF component suggests that the Po–V relationship may be similar in shape to that measured with ionic currents (Case A) but is shifted to more negative voltages (Case C). The Fast and Slow component data are consistent with this hypothesis but are inadequate to test the precise voltage dependence of Po. The data are also insufficient to specify how the model parameters should be altered to account for a change in Po. Case B, assuming a 12-fold increase in the equilibrium constant L, provides a reasonable first approximation that can account for both a shift in the Po–V relationship as well as the observed kinetics of slow charge movement.

Simulations of the modified allosteric scheme. When the C–O transition rates in the allosteric scheme (Figure 14) are modified, as specified by Case B, improved fits to the gating currents are generated. Figure 11A and Figure B, compares simulated currents to Ig evoked at +140 and +224 mV in response to pulses of different duration. The model accurately reproduces the amplitudes of ON and OFF currents, including the decrease in IgOFF amplitude that occurs with increased pulse duration at +224 mV (Figure 11 B). The model also fits the time course of IgOFF and accounts for the slowing of decay kinetics that accompanies increased pulse duration (Figure 11A and Figure B). The time course of OFF charge relaxation (QOFFQOFFss) after +140 mV pulses are plotted on a semilog scale in Figure 11 C, and are well fit at all pulse durations. Thus, the model accurately reproduces the kinetics and amplitudes of the three OFF components. The model can account for the slow time constants of both ON and OFF charge movement at all voltages (Figure 10 F); however, the amplitude of QpSlow is underestimated at low voltages (Figure 10 E, Case B). This point is illustrated in Figure 11 D, which compares Qp(t) at different voltages to QON(t) generated by the model. Both the time course and amplitude of Qp are well fit at V >= +140 mV; however, the slow component predicted by the model at lower voltages is reduced in comparison to the data.

Simulation of gating admittance. To further test the above conclusions, gating currents were simulated in response to a sinusoidal voltage command and compared with admittance analysis results. The CgV relationship is compared with the simulations for Cases A and B (solid lines) in Figure 11 E. Dashed lines indicate the QO'V and QC'–V relationships specified by the model. These relationships are the main determinants of Cg–V since they reflect fast voltage-sensor movement. At voltages below +100 mV where channels are closed, Cg approximates QC'. At higher voltages, Cg represents an average of QC' and QO' weighted by Po. Thus, Cg decreases at positive voltages (approaching QO') when Po is increased (compare Cases A and B). Case A overestimates Cg, suggesting that it underestimates Po. However, as the Po–V relationship is shifted (Case B), the model better approximates the peak amplitude and peak voltage of Cg. The effect of Po on the shape of the Cg–V relationship explains why the mean peak voltage of Cg (+127 mV) is more negative than the half-activation voltage of the Qfast–V relationship (+155 mV).


  Discussion
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
Acknowledgements
References

Examination of gating currents evoked from mSlo Ca2+-activated K+ channels in the absence of Ca2+ has revealed several components of charge movement associated with voltage-dependent gating. We have shown that these results are consistent with an allosteric voltage-gating scheme (Figure 14) that was proposed in the preceding article to account for the kinetic and steady-state properties of mSlo IK in 0 Ca2+ (Horrigan et al. 1999 ). Indeed, many of our experiments were designed to test this model. But before discussing these conclusions concerning the allosteric scheme, it is useful to review our results from a more general perspective. The gating current data lead to several model-independent conclusions and allow many alternative gating schemes to be ruled out.

Previous Models of BK Channel Gating
BK channel gating has been extensively studied at the single channel level (Barrett et al. 1982 ; Magleby and Pallotta 1983a ; Magleby and Pallotta 1983b ; McManus and Magleby 1991 ; Moczydlowski and Latorre 1983 ; Rothberg and Magleby 1998 ). Kinetic analysis reveals complex dwell-time distributions indicating the presence of multiple open and closed states. Based on such analysis, gating schemes have been proposed that contain a parallel arrangement of open and closed states (McManus and Magleby 1991 ), superficially resembling the architecture of our allosteric voltage-gating scheme (Figure 14). However, it is important to recognize that these previous studies were performed in the presence of Ca2+, and that the gating schemes used to describe these data therefore contain Ca2+-bound states and Ca2+-dependent transitions. Thus, the kinetic complexity revealed by the single channel data isn't necessarily related to the mechanism of voltage-dependent gating. Indeed, most schemes derived from single-channel analysis fail to account for the ability of BK channels to open in the absence of Ca2+ binding. By examining mSlo channel gating in the absence of Ca2+, we have characterized this voltage-dependent pathway, thereby defining a boundary condition that must be satisfied by any complete model of BK channel gating and representing a subset of the states that are accessible in the presence of Ca2+.

A model of BK channel gating has been proposed by Cox, Cui, and Aldrich 1997 to account for the effects of voltage and Ca2+ on macroscopic mSlo ionic currents, including their ability to activate in the absence of Ca2+. The model assumes that mSlo channels undergo a single voltage-dependent transition between a closed and open conformation and that Ca2+ binding regulates this transition allosterically. This scheme is essentially a version of the MWC model (Monod et al. 1965 ; Figure 15) where channel opening represents an allosteric transition that alters the affinity of Ca2+-binding sites and is also voltage dependent.



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 15. Scheme II

As in the McManus and Magleby model (1991), Figure 15 contains many states representing different Ca2+-bound versions of the closed and open conformation. However, in the absence of Ca2+, Figure 15 reduces to a two-state model with a single voltage-dependent transition between a unliganded closed and open state (highlighted above). By assuming that voltage-dependent activation can be described by a two-state mechanism, Figure 15 implies that channel opening, voltage-sensor movement, and changes in Ca2+ binding-site affinity all occur during a concerted allosteric transition. Our results demonstrate that a more complicated scheme is required to explain voltage-dependent gating and therefore imply that the interaction of Ca2+ with the channel may also be more complicated than proposed in Figure 15 (see Horrigan et al. 1999 ). In particular, the voltage-dependent C to O transition in Figure 15 does not consist of a completely cooperative (concerted) step, although considerable cooperativity, as formulated by the allosteric voltage-gating model, is present.

The preceding paper examines several properties of mSlo IK that are inconsistent with a two-state model of voltage gating (Horrigan et al. 1999 ). Single channel analysis of mSlo in 0 Ca2+ also provides evidence for multiple closed and open states (Talukder and Aldrich 1998 ). The gating current analysis presented here supports the conclusion that mSlo gating is a multistate process even in the absence of Ca2+. A two-state model of voltage-dependent activation requires that charge movement and channel opening occur simultaneously and therefore exhibit identical kinetic and steady-state properties. In other words, Ig should relax with the same near-exponential kinetics of IK, and the voltage dependence of steady-state charge movement (Q) and open probability (Po) should be identical. Instead, we observe multiple kinetic components of ON and OFF charge movement with major components of both preceding the relaxation of IK. In addition, the normalized Q–V and Po–V relationships are not superimposable (see also Stefani et al. 1997 ). These results indicate that mSlo channel opening cannot be represented by a concerted transition, and that the MWC model (Figure 15) is therefore an oversimplification in this regard, although it captures many of the major features of mSlo behavior.

Fast Ig: Evidence for a Two-state Model of Voltage-sensor Movement
Although the overall response of mSlo channels to voltage is complex, gating currents suggest that the movement of individual voltage sensors can be described by a simple two-state model when channels are closed.

Ig evoked during a voltage step exhibits a prominent fast component (IgFast) representing a majority of ON charge. This fast charge is also detected as a voltage-dependent component of membrane capacitance measured in response to a sinusoidal voltage command, thereby ruling out the possibility that leak subtraction or voltage clamp artifacts contribute to rapid current transients measured in response to large voltage steps. Both admittance analysis and the response to voltage steps indicate that fast gating charge can move at voltages where Po is normally low, and relaxes roughly 100-fold faster than the time constant of IK activation. IgFast decays with exponential kinetics during a time when few channels have opened. The relaxation of OFF current is also fast and single-exponential after brief pulses that open few channels. These results demonstrate that closed unliganded mSlo channels can undergo rapid voltage-dependent transitions.

Because the majority of ON charge moves rapidly, we assume that IgFast can be attributed to voltage-sensor movement. The exponential kinetics of IgFast and lack of a rising phase are consistent with a two-state model in which voltage sensors undergo a transition between a resting (R) and activated (A) conformation. The observation that the Qfast–V relationship is fit by a single Boltzmann function also supports a two-state model. In addition, the time constant of fast Ig relaxation ({tau}gFast) exhibits a bell-shaped voltage dependence that can be fit by the inverse sum of two exponential functions, as predicted for a two-state model in which forward and backward rate constants are voltage dependent.

mSlo channels assemble as homotetramers (Shen et al. 1994 ) and are therefore presumed to contain identical voltage sensors in each subunit. Thus, the simple behavior of IgFast is consistent not only with a two-state model of voltage-sensor movement but also with the idea that voltage sensors act independently. However, interactions between voltage sensors cannot be ruled out simply based on the kinetic and steady-state properties of fast charge movement. While it is true that such interactions could lead to multiexponential IgFast kinetics and a non-Boltzmann Q–V, more subtle effects are also possible. For instance, a model that assumes four voltage sensors move in a concerted manner would also predict two-state behavior, the difference being that the Qfast–V would be fit by a Boltzmann function with equivalent charge (zFast) of 4zJ for a concerted model versus zJ for an independent scheme. To distinguish these two possibilities requires an independent estimate of the fast charge per channel (qfast). Stefani et al. 1997 have reported a total charge (qT) of 4.4 ± 0.8e per channel (mean ± SD, n = 3) for hSlo based on measurements of Ig and ionic current density in different patches from the same oocyte. Although this estimate is not precise and includes both fast and slow charge, its magnitude argues against a concerted model, since zFast = 0.59 e determined for mSlo is much smaller than qT. An independent model would predict a fast charge of 4zFast = 2.36 e, much closer to qT. The relationship between fast charge movement and channel activation, discussed below, also argues against a concerted model of voltage sensor movement and is consistent with an independent scheme. However, uncertainty in some of these measurements, such as the estimate of qfast, prevents us from completely ruling out interaction between voltage sensors.

The Coupling of Voltage-sensor Movement to Channel Activation
Since the decay of IgFast is much faster than the activation of IK, we considered the possibility that fast charge movement might be unrelated to channel activation. An early component of charge movement has been described in Shaker K channels (Sigg et al. 1999 ; Stefani and Bezanilla 1996 ) and squid Na channels (Forster and Greeff 1992 ) that relaxes rapidly (Shaker: {tau} < 10 µs, Na channel: {tau} < 25 µs) and represent <10% of the total gating charge. The speed and small magnitude of this early charge movement suggest it could represent transitions that are not important for channel activation. IgFast described for mSlo is only severalfold slower than these early components and exhibits a similar equivalent charge. However, in the case of mSlo, several lines of evidence support the idea that fast charge movement is coupled to channel activation.

In contrast to the "early" charge movement in Shaker and Na channels, IgFast represents a majority of ON charge (Figure 6 B). In addition, the estimated fast charge per channel Qfast = 4zJ = 2.36 e (assuming independent voltage sensors) is similar to the equivalent charge that characterizes the maximum voltage dependence of Po in 0 Ca (z(Po) = 2.0 e; Horrigan et al., 1999). Thus, the magnitude of IgFast is consistent with the idea that fast charge movement is important for mSlo channel activation.

The kinetic relationship between IgFast and IK also argues that fast charge movement reflects transitions in the activation pathway. Fast charge movement and the delay in IK activation occur on similar time scales. An example in Figure 3 D shows that IgON decays at the same time that IK achieves an exponential time course. Thus, the achievement of a maximal rate of IK activation appears correlated with equilibration of fast gating charge. IK also exhibits a multiexponential rate of increase during the delay (Horrigan et al. 1999 ), supporting the idea that voltage-sensor transitions are not concerted. If the delay in IK depends only on the transitions that give rise to IgFast, then we have previously argued that the delay duration ({Delta}t) should be roughly proportional to {tau}gFast (Horrigan et al. 1999 ). Consistent with this prediction, {tau}gFast (Figure 4) and {Delta}t (Horrigan et al. 1999 ) exhibit similar bell-shaped voltage dependencies that can be characterized by an equivalent charge of 0.55e and peak voltages of 136 and 153 mV, respectively.

Finally, the Qfast–V and Po–V relationships, defining the voltage dependence of fast charge movement and IK activation, respectively, activate over a similar voltage range, consistent with the idea that these two processes are coupled. We have also shown that the normalized GK–V relationship can be approximated by raising the Qg–V relationship to the 4th power (Figure 1 B). As discussed below, an approximate 4th power relationship between Qfast–V and Po–V is predicted by many schemes that assume Po is enhanced by the activation of four voltage sensors. The relationship between Q–V and G–V is an important test of any voltage-dependent model, but experimental factors limit the interpretation of these data in the case of mSlo. The precise relationship because Q–V and G–V is unclear, owing to the likelihood that gating is altered under the conditions where gating currents are measured.

Sequential Voltage-gating Schemes
Taken together, the above observations indicate that the conformational changes underlying fast charge movement are involved in mSlo channel activation. Therefore, any plausible gating scheme must include a pathway that allows rapid voltage-dependent transitions to occur before channels open. The properties of fast charge movement are consistent with these closed-state transitions, arising from the activation of four independent and identical voltage sensors. Two sequential gating schemes incorporating such a mechanism are considered below (Figure 14 and Figure 15) and can reproduce many features of fast charge movement, but can be ruled out based on their failure to account for slow charge movement. These arguments parallel those in the preceding paper based on IK measurements (Horrigan et al. 1999 ), and lead to similar conclusions as to the requirement for an allosteric model.

One of the simplest schemes that can account for the properties of IgFast is the Hodgkin-Huxley (HH) model (Figure 16).



View larger version (6K):
[in this window]
[in a new window]
 
Fig. 16. Scheme III

The HH scheme assumes channels are open when all four voltage sensors are activated and predicts a 4th power relationship between the Q–V and G–V relationships:

(19)

where Q represents the normalized charge distribution defined by the equilibrium constant for voltage-sensor activation (Q = J/(1 + J)). As noted previously, the observed relationship between Qg–V and GK–V appears consistent with this prediction. However, Figure 16 is inadequate because it cannot account for the presence of both fast and slow components of mSlo charge movement. Similarly, the HH scheme cannot reproduce both the brief delay and slow exponential relaxation that characterize IK activation kinetics (Horrigan et al. 1999 ).

Models that assume voltage-sensor activation is followed by a distinct opening transition have proven useful in describing the behavior of channels that deviate from the predictions of the HH scheme (Koren et al. 1990 ; Perozo et al. 1992 ; Schoppa and Sigworth 1998 ; Sigworth 1994 ; Zagotta and Aldrich 1990 ; Zagotta et al. 1994a ). Such models can account for the presence of fast and slow components of gating current as well as nonsigmoidal IK activation kinetics (Horrigan et al. 1999 ; Ledwell and Aldrich 1999 ; Smith-Maxwell et al. 1998a ; Smith-Maxwell et al. 1998b ). Figure 17 (below) assumes that channels can undergo a rate-limiting C–O transition after four independent and identical voltage sensors are activated.



View larger version (6K):
[in this window]
[in a new window]
 
Fig. 17. Scheme IV

Figure 17 predicts an approximate 4th power relationship between G–V and Q–V described by the expression

(20)

where QC, the closed channel charge distribution, defines the voltage dependence of fast charge movement (QC = J/(1 + J)). As illustrated in Figure 12 A, this model can approximate the observed relationship between the Qfast–V and GK–V for mSlo. Figure 17 can also account for a slow component of ON charge movement (Figure 12 B) but, as discussed below, cannot reproduce some important aspects of gating current behavior. Similarly, Figure 17 can approximate the time course of mSlo IK but does not account for the complex voltage dependence of IK relaxation kinetics and open probability (Horrigan et al. 1999 ).

The Slow Component of ON Charge Movement
A slow component of ON charge movement (QpSlow) is detected as an increase in QOFF after pulses of increasing duration. QpSlow develops with the exponential kinetics of IK activation at depolarized voltages, suggesting that activation and slow charge movement are limited by the same transitions. We have shown that the allosteric voltage-gating scheme (Figure 14) can reproduce both the kinetics and voltage dependence of QpSlow (Figure 10E and Figure F). An important conclusion of this analysis is that QpSlow represents not only charge moved during the C–O conformational change but a reequilibration of voltage sensors that is limited by channel opening. The allosteric model predicts that transitions among open states (O–O) can contribute to slow charge movement, since voltage sensors reequilibrate after channels have opened. However, a contribution of voltage-sensor activation to slow charge movement does not require a model with multiple open states.

Figure 17 provides an example of a mechanism by which closed-state transitions contribute to both fast and slow charge movement. Fast Ig is evoked as voltage sensors initially equilibrate between resting (R) and activated (A) while the channel is closed. As channels open, this equilibrium is perturbed because channels can only open when all four voltage sensors are activated. In other words, opening stabilizes the activated voltage sensor, as in the allosteric model. However, in the case of Figure 17, the establishment of a new voltage-sensor equilibrium can only involve transitions between R and A while the channel is closed.

Figure 17 can reproduce the time course of Qp measured at +140 mV, including a large slow component (Figure 12 B). In addition, Figure 17 predicts a QpSlow–V relationship (Figure 12 C, solid lines) similar to that produced by the allosteric model (Figure 12 C, dashed lines). As with the allosteric scheme (Figure 14), a portion of QpSlow represents the charge assigned to the C–O transition (z{isin} = 0.32 e). The bell-shaped voltage dependence of QpSlow predicted by Figure 17 demonstrates that closed-state transitions also contribute to slow charge movement.

In summary, the presence of fast and slow components of ON charge, and their relationship to the time course of IK activation, suggest that the activation pathway must, at minimum, contain a rate-limiting step that is preceded by one or more rapid voltage-dependent transitions. The kinetics and voltage dependence of IgFast and the delay in IK activation, the voltage dependence of Qfast and GK together with the tetrameric structure of the channel further suggest that the rapid transitions may be described by the movement of four independent and identical voltage sensors. Therefore, Figure 17 provides the simplest model that can account for these basic features of the ionic and ON gating current data. However, as discussed below, the properties of OFF charge movement are inconsistent with Figure 17 and indicate that the activation pathway must include multiple open states.

OFF Charge Movement: Evidence for Multiple Open States
OFF currents, recorded after brief voltage pulses, decay with a single-exponential time course. Such a response is predicted by Figure 17 and is consistent with activated voltage sensors relaxing rapidly back to a resting state when channels are closed. However, Figure 17 also predicts that, once channels open, the decay of IgOFF will be limited by the speed of channel closing (Zagotta et al. 1994b ). Therefore, as pulse duration increases, a slow component of OFF charge relaxation should be observed that develops with the time course of IK activation and decays with the kinetics of IK deactivation. At the same time that the slow component increases, the fast component of OFF current should decrease as the number of channels that are closed at the end of the pulse is reduced.

Contrary to the prediction of Figure 17, we observed three components of OFF charge movement. The Fast and Slow components are analogous to those predicted by Figure 17. However, the Medium component, representing a majority of QOFF when channels are maximally activated, provides evidence that channels can undergo transitions among open states. In response to pulses of increasing duration, QOFFfast decreases with approximately the same time course as IK activation while the two slower components, QOFFmed and QOFFslow, increase in parallel. QOFFfast is essentially eliminated under conditions that maximally activate mSlo channels (20 ms at +160 mV in 60 µM Ca2+; Figure 8 B), implying that the Slow and Medium components reflect the relaxation of open channels back to the resting closed state.

We have proposed that the relaxation of the Slow OFF component is limited by the speed of channel closing and, at -80 mV, primarily represents charge moved during the O–C conformational change. QOFFslow represents a minority of the total OFF charge (Figure 10 B), consistent with the notion that the O–C transition is weakly voltage dependent. Similarly, the time constants of slow charge movement ({tau}gSlow) and IK deactivation ({tau}(IK)) are weakly voltage dependent at negative voltages (Figure 10 F). However, the decay of the Slow component is approximately threefold slower than that of potassium tail currents. To account for this difference, we have suggested that channel closing is slowed under the ionic conditions that are used to measure gating currents. The {tau}gSlow–V relationship is similar in shape to the {tau}(IK)–V relationship and can be fit by the allosteric model if the forward rate constants from C to O are increased while the backward rates are decreased relative to those used to describe IK. Such a change requires a 12-fold increase in the C–O equilibrium constant ({Delta}{Delta}G = 2.48 kT), producing a change in the Po–V relationship that appears consistent with the observed voltage dependence of QOFF components and QpSlow (Figure 10, Case B).

The Medium component of OFF charge relaxes ninefold faster than the Slow component and threefold faster than IK tail currents. Thus, regardless of the effect of ionic conditions on channel gating, the Medium component appears to relax faster than channel closing, implying that voltage sensors can move when channels are open. The similar voltage dependence of {tau}M and {tau}F (Figure 4 B2 and 9 C) supports the idea that the Medium component represents voltage-sensor movement. Thus, any plausible gating scheme must include multiple open states with rapid voltage-dependent transitions between them. The voltage dependence of Po leads to the same conclusion (Horrigan et al. 1999 ). The parallel development of Slow and Medium components indicate that once a channel is open, OFF charge relaxation can be described by a constant ratio of QOFFmed and QOFFslow. This behavior is consistent with the idea that equilibration of channels among different open states is fast relative to the speed of IK activation.

A Sequential Scheme with Multiple Open States
A sequential scheme, represented in general form below (Figure 18), could account for Medium and Slow components of OFF charge relaxation, provided transitions among open states are fast compared with the transition from O to C.



View larger version (4K):
[in this window]
[in a new window]
 
Fig. 18. Scheme V

However, such a model is inconsistent with the voltage dependence of steady-state activation. Po is weakly voltage dependent at limiting negative voltages, consistent with a charge of 0.4e assigned the Co–Oo transition in the allosteric scheme (Horrigan et al. 1999 ). For Figure 18 to reproduce this limiting voltage dependence, a total charge of 0.4e must be assigned to the transition between Co and Oo, inconsistent with the assumption that closed-state transitions from Co to Cm represent the activation of four voltage sensors carrying a total charge of 4zfast = 2.36 e.

Figure 18 also appears inconsistent with the relative amplitudes of various ON and OFF charge movement components. For example, we observe that the Medium component of OFF charge is two- to threefold larger than the Slow component measured at -80 mV (Figure 7 E and 8 A). Thus,

(21)

In addition, the fast component of ON charge is larger than the slow component at all voltages studied:

(22)

It can be shown (below) that Figure 18 cannot account for these observations if Equation 21 and Equation 22 are valid at voltages where Po >= 1/2. This last condition cannot be verified directly, but appears reasonable since Equation 21 and Equation 22 are true at +224 to +240 mV (Figure 6 B and 8 A), whereas the half-activation voltage of the GK–V relationship is +190 mV (Horrigan et al. 1999 ). Moreover, we have argued that Po may increase under the conditions where gating currents are measured.

The amplitudes of the different charge movement components for either the allosteric model or Figure 18 can be expressed in terms of QC, QO, and Po as specified by Equation 11Equation 12Equation 13Equation 14Equation 15Equation 16Equation 17Equation 18. Therefore, by substituting Equation 11 and Equation 12, Equation 21 can be rewritten:

(23)

where V is the pulse voltage and HP is the holding potential (-80 mV). Solving for QO(V), we obtain:

(24)

Equation 22 can also be rewritten by substituting Equation 14 and Equation 18.


(25)

Combining Equation 24 and Equation 25:

(26)

For a sequential model like Figure 18, we can further assume:

(27)

Finally, combining Equation 26 and Equation 27:

(28)

Equation 28 reduces to Po(V) < 1/2, indicating that Figure 18 cannot account for the relative amplitude of ON and OFF charge components while also assuming Po >= 1/2.

The Allosteric Voltage-gating Scheme
Figure 18 assumes that closed- and open-state transitions occur sequentially and must therefore represent distinct conformational events. An alternative, represented by the allosteric model, is that C–C and O–O transitions represent the same conformational events, i.e., voltage-sensor movement. The kinetics and voltage dependence of the Fast and Medium components of OFF charge movement (Figure 9 C) are consistent with both C–C and O–O transitions representing voltage-sensor movement, differing only in that the equilibrium constant for voltage-sensor activation is increased when channels open. As demonstrated in this study and in the preceding article, the allosteric model can account for many other properties of mSlo gating in 0 Ca2+.

The allosteric model is mechanistically similar to Figure 17 in that it assumes channels undergo only two types of conformational change: voltage-sensor activation and channel opening. Voltage sensors are assumed to move rapidly and independently in each subunit. Channel opening is relatively slow, weakly voltage dependent, and assumed to represent a concerted transition. Like Figure 17, the allosteric model assumes channel opening stabilizes the activated voltage sensor. Thus, opening results in a slow component of charge movement that is limited by the speed of channel opening but largely represents voltage-sensor charge movement. Unlike Figure 17, the coupling of voltage-sensor activation to channel opening is not an obligatory process but rather an allosteric interaction. Therefore, voltage sensors can move when channels are open, accounting for the Medium component of IgOFF, and channels can open when voltage sensors are not activated, accounting for the weak voltage dependence of Po measured at negative voltages (Horrigan et al. 1999 ).

Although the allosteric model allows channels to open when voltage sensors are in a resting state, it predicts that they are most likely to open when all four are activated. Consequently, channels pass through multiple closed states before opening, consistent with the presence of a delay in IK activation (Horrigan et al. 1999 ). Similarly, the allosteric scheme can account for an approximate 4th power relationship between Q–V and G–V. The model predicts the following relationship between Po and QC.


(29)

When L is small and D >> 1, as determined in the preceding article (L = 2 x 10-6, D = 17, 0 Ca2+), this expression can be approximated as:

(30)

Finally, the allosteric scheme can account for the presence of three components of OFF charge movement as well as the relative amplitudes of various ON and OFF components. In contrast to Figure 18, the allosteric model predicts that pathways traversed during channel activation and deactivation are different. Activation involves fast voltage-sensor movement as channels undergo closed-state transitions before opening. Deactivation involves movement of the same voltage sensors as channels undergo open-state transitions before closing. Because the same voltage sensors are moved during open- and closed-state transitions, the rapid components of ON (Qfast) and OFF charge (QOFFfast, QOFFmed) are of similar amplitude while the slow components are smaller. Furthermore, the relative amplitudes of QOFFmed and QOFFslow change with repolarization voltage in a manner specifically predicted by the allosteric scheme (Figure 8C and Figure D).

Comparison with Previous Studies of BK Channel Gating
mSlo ionic currents. The gating current data support the conclusion from the preceding paper that mSlo channel voltage gating in the absence of Ca2+ can be described by an allosteric scheme. Indeed, many of the model parameters that were derived to fit IK data required little or no adjustment to describe the gating currents (e.g., zJ, Vh(J), zL, D). One feature of the model that could not be determined accurately from IK measurements was the speed of transitions among open states. The Medium component of IgOFF provides a direct assay of these transitions and demonstrates that channel opening slows the relaxation of voltage sensors from A to R. The magnitude of this effect is consistent with our previous estimate of the allosteric factor D = 17, provided we assume channel opening almost symmetrically affects the forward and backward rate constants for the R–A transition (Figure 9 C). However the effect of channel opening on the forward rate was not measured; therefore, the value of D cannot be directly determined from the gating current data.

Single channel studies. Many of the properties of IK that implicate a model with multiple open states are observed only at extreme voltages (Horrigan et al. 1999 ). However, the Medium component of IgOFF demonstrates that the channels can undergo open-state transitions in response to moderate voltage stimuli, such as repolarization from +160 to 0 mV (Figure 8 D). This is important because it suggests that complex open-time distributions described previously for single BK channels may to some extent reflect the occupancy of multiple open states in the voltage-dependent activation pathway in addition to different Ca2+-bound open states.

Gating currents. BK channel gating currents have been described previously for hSlo by Stefani et al. 1997 , and many of their results, obtained in 0 Ca2+, are similar to those presented here for mSlo. For example, the major component of ON current decayed with a rapid exponential time course ({tau} = 57 ± 10 µs, at +200 mV). A slow component of charge movement was also detected but not examined in detail. The Q–V curve determined with brief 1-ms depolarizations was well fit by a Boltzmann function (z = 0.6 e, Vh = 190 ± 15 mV). In addition, the G–V curve was reported to be steeper than the Q–V, and charge movement was observed at voltages where most channels should be closed.

In addition to these similarities, there are important differences between our results and conclusions and those of Stefani et al. 1997 . The normalized G–V and Q–V for hSlo were reported to cross such that the G–V is negative to the Q–V at positive voltages. Such a cross-over, as pointed out by Stefani et al., implies that gating charge can move when channels are open. Although our results lead to a similar conclusion, we do not observe such a relationship between Q–V and G–V in the absence of Ca2+. And the allosteric model used to fit our data does not predict a cross-over even though it allows voltage-dependent open-state transitions. In addition, Stefani et al. conclude that only a small fraction of total gating charge must move before channels can open because (a) the Q–V crosses the foot of the G–V where Po is small, and (b) the Cole-Moore shift is weakly voltage dependent. We disagree with these observations and, for reasons discussed below, conclude that in the absence of Ca2+, the majority of gating charge moves before channels open. Stefani et al. also note that the onset of ionic current overlaps the decay of Ig (see Figure 3) and suggest this provides evidence for open-state charge movement. But most models, including those with a single open state, predict some overlap in the time course of macroscopic Ig and IK. Finally, Stefani et al. conclude that open-state transitions must be weakly voltage dependent to account for the observation that the major component of charge movement is fast. In contrast, we have shown that the allosteric scheme predicts a small slow component of charge movement even when C–C and O–O transitions are assumed to be equally voltage dependent. And the large Medium component of OFF charge movement provides evidence for significant open-state charge movement.

Some of the discrepancies between our results and those of Stefani et al. 1997 probably reflect real differences in the gating of mSlo and hSlo, but we suggest that the cross-over of Q–V and G–V reported for hSlo may also be affected by the conditions used to measure ionic currents. The Q–Vs for mSlo and hSlo are similar in shape based on Boltzmann fits (z = 0.59 e mSlo; z = 0.6 e hSlo), and the half-activation voltage for the hSlo Q–V is 35 mV more positive than that for mSlo (155 mV mSlo). In line with this difference, the GK–V for hSlo, measured in symmetrical K+ and 0 Ca, is roughly 30 mV more positive than that for mSlo (hSlo: Vh = 220 mV [ Diaz et al. 1998 ; Meera et al. 1996 ]; mSlo: Vh = 190 mV V [Horrigan et al. 1999 ). In contrast, the G–V described by Stefani et al. is shifted by -70 mV (Vh = 150 mV) relative to that previously reported for hSlo and is much steeper than that of mSlo based on Boltzmann fits (hSlo: z = 1.3–1.6 e [ Stefani et al. 1997 ]; mSlo: z = 1.0 e [Cui et al., 1997; Horrigan et al. 1999 ]). These differences may reflect the fact that the G–V reported by Stefani et al. was obtained in symmetrical Cs+ rather than K+.

Cs+ permeates BK channels poorly, allowing Cs+ currents to be recorded in patches where IK would be immeasurably large (Stefani et al. 1997 ). However, Cs+ is also known to alter BK channel gating and is likely to affect the G–V. Demo and Yellen 1992 studied Cs+ block of single BK channels and concluded that Cs+ occupancy of the pore destabilizes the closed state. This destabilization shifts the Po–V relationship to more negative voltages and changes its shape because Cs+ block is voltage dependent. Similarly the GCs–V relationship for hSlo is shifted to negative voltages relative to the GK–V for hSlo and is steeper than that for mSlo, consistent with an effect of Cs+ on channel gating (Demo and Yellen 1992 ). An increase in steepness of the G–V could occur if Cs+ occupancy of the pore is voltage dependent. Consistent with this possibility, the instantaneous Ics–V relationship recorded for mSlo (Figure 12 D) is highly nonlinear, indicating that Cs+ permeates more readily at very positive or negative voltages. Cs+ currents also activate slowly, failing to achieve a steady state during a 25-ms pulse (Stefani et al. 1997 ). These kinetics differ from those of IK measured for hSlo (Meera et al. 1996 ) or mSlo (Horrigan et al. 1999 ) in 0 Ca2+, providing additional evidence for an effect of Cs+ on gating and implying that the Gcs–V recorded with 25-ms pulses doesn't represent the steady-state G–V.

hSlo gating currents were also measured in the presence of internal Cs+, using isotonic external TEA to block the channel. Therefore, the use of Cs+ has the apparent advantage of allowing gating and ionic currents to be recorded with the same internal solution. However, the presence of internal Cs+ does not guarantee that channels gate identically when blocked by TEA or conducting Cs+. Aside from the possibility that TEA directly affects gating (see below), TEA may inhibit the effect of Cs+. Demo and Yellen 1992 found that BK channel block by either internal Cs+ or external TEA had no effect on Po, and they concluded that Cs+ could occupy at least two sites in the pore, only one of which affects gating. Thus, it is possible that internal Cs+ cannot occupy the critical gating site when the pore is blocked by TEA. This could explain why hSlo gating currents resemble those recorded for mSlo, and may help account for the observed crossing of Q–V and G–V for hSlo.

As discussed previously, mSlo gating may also be affected by the presence of external TEA, internal NMDG, or the absence of K+ in experiments measuring gating currents. Several properties of slow charge movement summarized in Figure 10 suggest that mSlo channels open more readily under the conditions where gating currents were measured. A 10-fold reduction in internal and external K+ had no appreciable effect on the GK–V (data not shown) but we cannot rule out the possibility that gating is altered by the complete removal of K+ or its replacement with NMDG. Stefani et al. 1997 found that fast Ig, evoked at voltages where channels do not open, were unaltered by the presence of TEA. We also saw no effect of TEA on fast charge movement measured with admittance analysis in 0 K+ (data not shown). However, these experiments do not rule out an effect of TEA on slow charge movement and channel opening.

Another factor that could contribute to an apparent cross-over between Q–V and G–V is the duration of the voltage pulses used to measure gating currents. The Q–V for hSlo was determined using 1-ms pulses and is therefore similar to Qfast and not a steady-state measurement. Stefani et al. 1997 state that the Q–V determined with longer pulses (10–20 ms) in 0 Ca2+ is shifted by -20 mV relative to Q1ms, an effect that is insufficient to account for the apparent cross-over of Q–V and G–V. We also observed little difference between the normalized Qss–V and Qfast–V curves in 0 Ca2+ (Figure 6 D). However, in the presence of 60 µM Ca2+, we observe a large 50 mV shift between Q–Vs determined with 1- or 20-ms pulses (Horrigan, F.T, and R.W. Aldrich, manuscript in preparation). Thus, a cross-over of Q–V and G–V reported by Stefani in 85 µM Ca2+ should be strongly influenced by the use of 1-ms voltage pulses.

Despite uncertainties as to the precise relationship between the Q–V and G–V, gating currents recorded for both mSlo and hSlo show that most charge movement in 0 Ca2+ is fast, indicating that most charge moves before channels open. The kinetics and voltage dependence of the delay in IK activation are also consistent with the idea the multiple voltage-dependent closed-state transitions, accounting for the bulk of charge movement, occur before channels open. Finally, the change in delay duration with prepulse voltage (Cole-Moore shift [ Cole and Moore 1960 ]) is, as stated by Stefani et al., less voltage dependent than that observed for Shaker. However, this is consistent with the overall weak voltage dependence and reduced gating charge of BK channels as compared with Shaker and does not imply that a small proportion of total charge moves before Slo channels open. The Cole-Moore shift for mSlo is well described by the allosteric model, which predicts most voltage sensors activate before channels open (Horrigan et al. 1999 ). The lack of a rising phase in Ig also shows that the earliest closed-state transitions are not weakly voltage dependent.

Detection of Allosteric Voltage Gating
In the preceding article (Horrigan et al. 1999 ), we discussed the possibility that the behavior of voltage-gated ion channels such as Shaker may be consistent with an allosteric voltage-gating scheme (Figure 14) like that used to describe mSlo. One reason to consider this possibility is that many, but not all, features of mSlo ionic currents can be adequately described by sequential gating schemes that have been used to describe a variety of other channels. The behaviors of mSlo that deviate from these conventional schemes and implicate the allosteric model are mainly observed under conditions of extreme voltage or low open probability. Therefore, it is possible that other channels operate through an allosteric mechanism but have not been studied under conditions that are necessary to test this model, which may be even more extreme in those channels than for mSlo.

Many of the gating current properties described here for mSlo can also be accounted for by sequential gating schemes containing a single open state such as Figure 17. However, deviations from the prediction of Figure 17 are more obvious for gating current than for ionic current. The Medium component of OFF charge movement, in particular, provides a direct indication of open-state transitions. Many voltage-gated K+ channels such as Shaker exhibit OFF currents that become slower as pulse duration is increased and channels open (Bezanilla et al. 1991 ; Bezanilla et al. 1994 ; Chen et al. 1997 ; Kanevsky and Aldrich 1999 ; Zagotta et al. 1994b ). However, unlike mSlo, the decay of OFF currents for open Shaker channels appear to be limited by the speed of channel closing as predicted by models like Figure 17 (Zagotta et al. 1994b ). This observation seems to argue against a model with multiple open states. However, as discussed below, the allosteric scheme can account for such results when the speed of voltage-sensor movement or the voltage dependence of open-state transitions is altered.

Two factors allow open-channel charge movement to be detected for mSlo. First, voltage-sensor movement is much faster than channel closing. This difference allows the Medium and Slow components of OFF charge movement to be distinguished and allows open-state transitions to occur before channels close. As discussed below, the relative speed of Ig and IK in channels such as Shaker might prevent detection of open-state transitions. Another factor that is important for detecting open-state transitions is the open-state charge distribution (QO). That is, the voltage dependence of open-state transitions must be such that repolarization to the holding potential causes a redistribution of channels among open states. The effect of repolarization voltage on QOFF components in Figure 8 shows that the Medium component is sensitive to the open-channel charge distribution. Therefore, a change in the voltage dependence of QO might alter the ability to detect open-state transitions. For example, if the allosteric factor D is increased, QO will be shifted to more negative voltages such that QOFFmed measured at -80 mV is reduced.

Figure 13 shows that a slowing of voltage-sensor kinetics reduces the ability to detect open-channel charge movement. IK and Ig were simulated in response to a 20-ms pulse to +240 mV as the forward and backward rates for voltage-sensor movement ({alpha}, ß) were both slowed 10-fold (10x) or 30-fold (30x) relative to those describing mSlo (1x). C–O transition rates and all equilibrium constants were unchanged (relative to Case B parameters). As voltage-sensor movement is slowed, the delay in IK activation increases (Figure 13 A) and gating currents are slowed (Figure 13 B). Under these conditions, IK and Ig resemble those evoked from a channel such as Shaker where IK activation kinetics are sigmoidal and ionic and gating currents relax on a similar time scale. Interestingly, a 30-fold slowing of voltage sensor movement also produces a "hook" in IgOFF, a feature that is also observed in Shaker Ig (Bezanilla et al. 1991 ; Chen et al. 1997 ; Perozo et al. 1992 ). The simulation predicts that Po approaches unity at the end of a pulse to +240 mV. Therefore, IgOFF represents the relaxation of open channels. The time course of QOFF(t) plotted on a semilog scale in Figure 13 C is biphasic when voltage-sensor movement is fast (1x) representing the Medium and Slow components of QOFF. However, kinetically distinct components of QOFF are not evident when voltage-sensor movement is slow (Figure 13 C, 10x and 30x). Similarly, the decay of IgOFF is much faster than that of IK when voltage sensors are fast (Figure 13 D, 1x). However, IgOFF and IK decay with similar kinetics when voltage sensors are slow (Figure 13 D, 10x and 30x). These results demonstrate that OFF charge movement can be limited by the speed of channel deactivation even when multiple open states are present.


 
View this table:
[in this window]
[in a new window]
 
Table 1. Allosteric Scheme Parameters


  Acknowledgements
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
Acknowledgements
References

This work was supported by a grant from the National Institutes of Health (N153294) and by a National Institute of Mental Health Silvio Conte Center for Neuroscience Research grant (MH48108). R.W. Aldrich is an investigator with the Howard Hughes Medical Institute.

Submitted: March 16, 1999; Revised: June 1, 1999; Accepted: June 7, 1999.

1used in this paper: BK channel, large-conductance Ca2+-activated K+ channel; NMDG, N-methyl-D-glucamine; TEA, tetraethylammonium
  References
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
Acknowledgements
References

Aggarwal, S.K., MacKinnon, R. (1996) Contribution of the 17 segment to gating charge in the Shaker K+ channel. Neuron. 16:1169-1177[Medline].

Armstrong, C.M. (1971) Interaction of tetraethylammonium ion derivatives with the potassium channels of giant axons. J. Gen. Physiol. 58:413-437[Abstract/Free Full Text].

Armstrong, C.M., Bezanilla, F. (1974) Charge movement associated with the opening and closing of the activation gates of the Na channels. J. Gen. Physiol. 63:533-552[Abstract/Free Full Text].

Barrett, J.N., Magleby, K.L., Pallotta, B.S. (1982) Properties of single calcium-activated potassium channels in cultured rat muscle. J. Physiol. (Lond.). 331:211-230[Medline].

Bezanilla, F., Perozo, E., Papazian, D.M., Stefani, E. (1991) Molecular basis of gating charge immobilization in Shaker potassium channels. Science. 254:679-683[Medline].

Bezanilla, F., Perozo, E., Stefani, E. (1994) Gating of Shaker K+ channels: II. The components of gating currents and a model of channel activation. Biophys. J. 66:1011-1021[Abstract].

Butler, A., Tsunoda, S., McCobb, D.P., Wei, A., Salkoff, L. (1993) mSlo, a complex mouse gene encoding "maxi" calcium-activated potassium channels. Science. 261:221-224[Medline].

Chen, F.S., Steele, D., Fedida, D. (1997) Allosteric effects of permeating cations on gating currents during K+ channel deactivation. J. Gen. Physiol 110:87-100[Abstract/Free Full Text].

Cole, K.S., Moore, J.W. (1960) Potassium ion current in the squid giant axon: dynamic characteristic. Biophys. J. 1:1-14[Medline].

Cox, D.H., Cui, J., Aldrich, R.W. (1997) Allosteric gating of a large conductance Ca-activated K+ channel. J. Gen. Physiol. 110:257-281[Abstract/Free Full Text].

Cui, J., Cox, D.H., Aldrich, R.W. (1997) Intrinsic voltage dependence and Ca2+ regulation of mslo large conductance Ca-activated K+ channels. J. Gen. Physiol 109:647-673[Abstract/Free Full Text].

Demo, S.D., Yellen, G. (1992) Ion effects on gating of the Ca(2+)-activated K+ channel correlate with occupancy of the pore. Biophys. J. 61:639-648[Abstract].

Diaz, L., Meera, P., Amigo, J., Stefani, E., Alvarez, O., Toro, L., Latorre, R. (1998) Role of the 17 segment in a voltage-dependent calcium-sensitive potassium (hSlo) channel. J. Biol. Chem. 273:32430-32436[Abstract/Free Full Text].

Diaz, F., Wallner, M., Stefani, E., Toro, L., Latorre, R. (1996) Interaction of internal Ba2+ with a cloned Ca2+-dependent K+ (hslo) channel from smooth muscle. J. Gen. Physiol. 107:399-407[Abstract].

Fedida, D., Maruoka, N.D., Lin, S. (1999) Modulation of slow inactivation in human cardiac Kv1.5 channels by extra- and intracellular permeant cations. J. Physiol. (Lond.). 515:315-329[Abstract/Free Full Text].

Fernandez, J.M., Bezanilla, F., Taylor, R.E. (1982) Distribution and kinetics of membrane dielectric polarization. II. Frequency domain studies of gating currents. J. Gen. Physiol. 79:41-67[Abstract].

Fernandez, J.M., Taylor, R.E., Bezanilla, F. (1983) Induced capacitance in the squid giant axon. Lipophilic ion displacement currents. J. Gen. Physiol. 82:331-346[Abstract].

Forster, I.C., Greeff, N.G. (1992) The early phase of sodium channel gating current in the squid giant axon. Characteristics of a fast component of displacement charge movement. Eur. Biophys. J. 21:99-116[Medline].

Gillis, K.D. (1995) Techniques for membrane capacitance measurements. In Sakmann B., Neher E., eds. In Single Channel Recording. NY, NY, Plenum Publishing Corp, 155-198.

Hamill, O.P., Marty, A., Neher, E., Sakmann, B., Sigworth, F.J. (1981) Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pflugers Arch. 391:85-100[Medline].

Herrington, J., Bookman, R.J. (1995) Pulse Control. Miami, University of Miami Press.

Horrigan, F.T., Cui, J., Aldrich, R.W. (1999) Allosteric voltage gating of potassium channels I: mSlo channel ionic currents in the absence of Ca2+. J. Gen. Physiol. 114:277-304[Abstract/Free Full Text].

Kanevsky, M., Aldrich, R.W. (1999) Determinants of voltage-dependent gating and open-state stability in the 18 segment of Shaker potassium channels. J. Gen. Physiol. 114:215-242[Abstract/Free Full Text].

Koren, G., Liman, E.R., Logothetis, D.E., Nadal, G.B., Hess, P. (1990) Gating mechanism of a cloned potassium channel expressed in frog oocytes and mammalian cells. Neuron. 4:39-51[Medline].

Ledwell, J.L., Aldrich, R.W. (1999) Mutations in the 17 region isolate the final voltage-dependent cooperative step in potassium channel activation. J. Gen. Physiol. 113:389-414[Abstract/Free Full Text].

Lindau, M., Neher, E. (1988) Patch-clamp techniques for time-resolved capacitance measurements in single cells. Pflugers Arch. 411:137-146[Medline].

Lollike, K., Borregaard, N., Lindau, M. (1995) The exocytotic fusion pore of small granules has a conductance similar to an ion channel. J. Cell Biol 129:99-104[Abstract].

Magleby, K.L., Pallotta, B.S. (1983a) Burst kinetics of single calcium-activated potassium channels in cultured rat muscle. J. Physiol. (Lond.). 344:605-623[Abstract].

Magleby, K.L., Pallotta, B.S. (1983b) Calcium dependence of open and shut interval distributions from calcium-activated potassium channels in cultured rat muscle. J. Physiol. (Lond.). 344:585-604[Abstract].

Matteson, D.R., Swenson, R.P., Jr. (1986) External monovalent cations that impede the closing of K channels. J. Gen. Physiol. 87:795-816[Abstract].

McManus, O.B., Magleby, K.L. (1991) Accounting for the Ca(2+)-dependent kinetics of single large-conductance Ca(2+)-activated K+ channels in rat skeletal muscle. J. Physiol. (Lond.). 443:739-777[Abstract].

Meera, P., Wallner, M., Jiang, Z., Toro, L. (1996) A calcium switch for the functional coupling between alpha (hslo) and beta subunits (KV, Ca beta) of maxi K channels. FEBS Lett. 382:84-88[Medline].

Miller, C., Latorre, R., Reisin, I. (1987) Coupling of voltage-dependent gating and Ba++ block in the high-conductance, Ca++-activated K+ channel. J. Gen. Physiol. 90:427-449[Abstract].

Moczydlowski, E., Latorre, R. (1983) Gating kinetics of Ca2+-activated K+ channels from rat muscle incorporated into planar lipid bilayers. Evidence for two voltage-dependent Ca2+ binding reactions. J. Gen. Physiol. 82:511-542[Abstract].

Monod, J., Wyman, J., Changeux, J.P. (1965) On the nature of allosteric transitions: a plausible model. J. Mol. Biol. 12:88-118[Medline].

Neyton, J. (1996) A Ba2+ chelator suppresses long shut events in fully activated high-conductance Ca2+-dependent K+ channels. Biophys. J. 71:220-226[Abstract].

Neyton, J., Pelleschi, M. (1991) Multi-ion occupancy alters gating in high-conductance, Ca2+-activated K+ channels. J. Gen. Physiol. 97:641-665[Abstract].

Perozo, E., Papazian, D.M., Stefani, E., Bezanilla, F. (1992) Gating currents in Shaker K+ channels. Implications for activation and inactivation models. Biophys. J. 62:160-168, discussion 169–171[Abstract].

Press, W.H., Teukolsky, S.A., Vetterling, W.T., Flannery, B.P. (1992) Numerical Recipes in C: the Art of Scientific Computing, Vol. 2. Cambridge, Cambridge University Press, pp. 717–722.

Rodriguez, B.M., Sigg, D., Bezanilla, F. (1998) Voltage gating of Shaker K+ channels. The effect of temperature on ionic and gating currents. J. Gen. Physiol 112:223-242[Abstract/Free Full Text].

Rothberg, B.S., Magleby, K.L. (1998) Kinetic structure of large-conductance Ca2+-activated K+ channels suggests that the gating includes transitions through intermediate or secondary states. A mechanism for flickers. J. Gen. Physiol. 111:751-780[Abstract/Free Full Text].

Sala, S., Matteson, D.R. (1991) Voltage-dependent slowing of K channel closing kinetics by Rb+. J. Gen. Physiol. 98:535-554[Abstract].

Schoppa, N.E., McCormack, K., Tanouye, M.A., Sigworth, F.J. (1992) The size of gating charge in wild-type and mutant Shaker potassium channels. Science. 255:1712-1715[Medline].

Schoppa, N.E., Sigworth, F.J. (1998) Activation of Shaker potassium channels. III. An activation gating model for wild-type and V2 mutant channels. J. Gen. Physiol. 111:313-342[Abstract/Free Full Text].

Shen, K.Z., Lagrutta, A., Davies, N.W., Standen, N.B., Adelman, J.P., North, R.A. (1994) Tetraethylammonium block of Slowpoke calcium-activated potassium channels expressed in Xenopus oocytes: evidence for tetrameric channel formation. Pflugers Arch. 426:440-445[Medline].

Sigg, D., Qian, H., Bezanilla, F. (1999) Kramers' diffusion theory applied to gating kinetics of voltage-dependent ion channels. Biophys. J 76:782-803[Abstract/Free Full Text].

Sigworth, F.J. (1994) Voltage gating of ion channels. Q. Rev. Biophys. 27:1-40[Medline].

Smith-Maxwell, C.J., Ledwell, J.L., Aldrich, R.W. (1998a) Role of the 17 in cooperativity of voltage-dependent potassium channel activation. J. Gen. Physiol. 111:399-420[Abstract/Free Full Text].

Smith-Maxwell, C.J., Ledwell, J.L., Aldrich, R.W. (1998b) Uncharged 17 residues and cooperativity in voltage-dependent potassium channel activation. J. Gen. Physiol. 111:421-439[Abstract/Free Full Text].

Stefani, E., Bezanilla, F. (1996) Early events in voltage gating. Biophys. J. 70:143. (Abstr.).

Stefani, E., Ottolia, M., Noceti, F., Olcese, R., Wallner, M., Latorre, R., Toro, L. (1997) Voltage-controlled gating in a large conductance Ca2+-sensitive K+channel (hslo). Proc. Natl. Acad. Sci. USA. 94:5427-5431[Abstract/Free Full Text].

Swenson, R.J., Armstrong, C.M. (1981) K+ channels close more slowly in the presence of external K+ and Rb+. Nature. 291:427-429[Medline].

Talukder, G., Aldrich, R.W. (1998) Intermediate states in the gating of single mSlo calcium activated potassium channels in the absence of calcium. Biophys. J. 75:A217. (Abstr.).

Taylor, R.E., Bezanilla, F. (1979) Comments on the measurement of gating currents in the frequency domain. Biophys. J. 26:338-340[Medline].

Tzounopoulos, T., Maylie, J., Adelman, J.P. (1995) Induction of endogenous channels by high levels of heterologous membrane proteins in Xenopus oocytes. Biophys. J. 69:904-908[Abstract].

Wang, Z., Zhang, X., Fedida, D. (1999) Gating current studies reveal both intra- and extracellular cation modulation of K+ channel deactivation. J. Physiol. (Lond.). 515:331-339[Abstract/Free Full Text].

Wei, A., Solaro, C., Lingle, C., Salkoff, L. (1994) Calcium sensitivity of BK-type KCa channels determined by a separable domain. Neuron. 13:671-681[Medline].

Yeh, J.Z., Armstrong, C.M. (1978) Immobilisation of gating charge by a substance that simulates inactivation. Nature. 273:387-389[Medline].

Zagotta, W.N., Aldrich, R.W. (1990) Voltage-dependent gating of Shaker A–type potassium channels in Drosophila muscle. J. Gen. Physiol. 95:29-60[Abstract].

Zagotta, W.N., Hoshi, T., Aldrich, R.W. (1994a) Shaker potassium channel gating. III: Evaluation of kinetic models for activation. J. Gen. Physiol. 103:321-362[Abstract].

Zagotta, W.N., Hoshi, T., Dittman, J., Aldrich, R.W. (1994b) Shaker potassium channel gating. II: Transitions in the activation pathway. J. Gen. Physiol 103:279-319[Abstract].