Intracellular blockade by quaternary ammonium (QA) molecules of many potassium channels is state dependent, where the requirement for channel opening is evidenced by a time-dependent component of block in the macroscopic record. Whether this is the case for Ca2+- and voltage-activated potassium (BK) channels, however, remains unclear. Previous work (Li, W., and R.W. Aldrich. 2004. J. Gen. Physiol. 124:43–57) tentatively proposed a state-dependent, trapping model, but left open the possibility of state-independent block. Here, we found BK channel blockade by a novel QA derivative, bbTBA, was time dependent, raising the possibility of state-dependent, open channel block. Alternatively, the observed voltage dependence of block could be sufficient to explain time-dependent block. We have used steady-state and kinetic measurements of bbTBA blockade in order to discriminate between these two possibilities. bbTBA did not significantly slow deactivation kinetics at potentials between −200 and −100 mV, suggesting that channels can close unhindered by bound bbTBA. We further find no evidence that bbTBA is trapped inside BK channels after closing. Measurements of steady state fractional block at +40 mV revealed a 1.3-fold change in apparent affinity for a 33-fold change in Po, in striking contrast to the 31-fold change predicted by state-dependent block. Finally, the appearance of a third kinetic component of bbTBA blockade at high concentrations is incompatible with state-dependent block. Our results suggest that access of intracellular bbTBA to the BK channel cavity is not strictly gated by channel opening and closing, and imply that the permeation gate for BK channels may not be intracellular.
The mechanism of ion channel opening and closing is widely thought of in terms of a gate that permits or restricts access to the inner pore via a reversible conformational change. Much attention has been focused on locating the gate of K+ channels, as a first step toward understanding the conformational changes that underlie the gating process. Clay Armstrong's innovative use of intracellular organic blockers that could block K+ channels only when they were open (Armstrong, 1966, 1969, 1971) produced the first strong evidence for an intracellular activation gate. The classic features of state-dependent (open channel) block defined in this work (Armstrong, 1971) include (a) time-dependent block of macroscopic current, where steady-state block reaches equilibrium only after open probability is sufficiently high, (b) a “foot in the door” effect of blocker on channel closing, where internally bound quaternary ammonium (QA) ions of a sufficient size were found to delay closure of the activation gate, and (c) “trapping” of a blocker molecule inside the channel pore after closure of the activation gate. These observations strongly suggested that intracellular blockers required channel opening in order to access (and exit) the pore, and implied that the gating structure regulating access of the blocker to the pore must be located on the intracellular face of the channel. Numerous functional studies of K+ channels (Armstrong and Hille, 1972; MacKinnon and Miller, 1989; MacKinnon and Yellen, 1990; MacKinnon et al., 1990; Yellen et al., 1991; Heginbotham et al., 1992; Liu et al., 1997; Shin et al., 2001) and structural studies of bacterial (Doyle et al., 1998; Jiang et al., 2002a; Jiang et al., 2003; Kuo et al., 2003) and mammalian (Long et al., 2005) K+ channels have confirmed Armstrong's original hypothesis.
The structure of the mammalian voltage-gated K+ channel Kv 1.2 (Long et al., 2005) exemplifies the detailed organization of an archetypal K+ channel pore. The S5 and S6 transmembrane domains from each of four subunits define the central pore, a large water-filled cavity constrained at the extracellular end by the narrow selectivity filter, and at the intracellular end by the “bundle crossing,” where the S6 helices come into close apposition. For the majority of K+ channels, permeation of K+ ions appears to be governed by constriction and enlargement of the S6 aperture. On the other hand, it is clear that some related ion channels employ a different gating mechanism. The case for cyclic nucleotide-gated (CNG) channels is particularly compelling, where Flynn and Zagotta (2001) found that although access of internally applied MTSET to pore residues depends upon the state of the channel, those residues are equally accessible to internal Ag+ (roughly equivalent to a K+ ion in size) in the open or closed state. The authors propose that although the larger MTS reagent is gated by conformational changes of the S6 bundle crossing, permeant ions may be gated at the selectivity filter for CNG channels. Similarly, small-conductance Ca2+-activated K+ (SK) channels (Bruening-Wright et al., 2002), ATP-sensitive inward rectifier K+ (KATP) channels (Proks et al., 2003), and KcsA K+ channels (Cordero-Morales et al., 2006) may also be gated at the selectivity filter.
For BK Ca2+- and voltage-activated K+ channels, the location of the permeation gate remains an open question. Although the structure of a bacterial Ca2+-gated K+ channel (MthK) is known (Jiang et al., 2002a), and functional evidence from BK channels is consistent with motion of the intracellular S6 domains during gating (Niu et al., 2004), the mechanism of BK channel gating is far from fully understood. In some ways, intracellular block of BK channels appears similar to classical open channel block. Most notably, block of BK channels by intracellular Shaker ball peptide is state dependent; the peptide blocks only open channels and must exit before channels can close (unpublished data). Several other observations could also be consistent with an intracellular gate for BK channels: mutations in S6 disrupt activation gating (Piskorowski, 2005), BK channels can close when the selectivity filter is occupied by either Na+ or Ba2+ (Miller et al., 1987; Neyton and Pelleschi, 1991), and high concentrations of external K+ have no obvious slowing effect on deactivation of BK channels (Neyton and Pelleschi, 1991; Demo and Yellen, 1992).
In spite of these observations, however, there is equally compelling evidence suggesting that the BK channel gate may lie elsewhere. First, QA molecules, which have been shown functionally to promote C-type inactivation (a form of selectivity filter gating in Kv channels; Choi et al., 1991) and are proposed to induce constriction of the Kcsa selectivity filter in structural studies (Lenaeus et al., 2005), also speed closing of BK channels (Li and Aldrich, 2004), raising the possibility that BK channels gate at the selectivity filter via a C-type collapse mechanism. Also, permeant thallium favors Kcsa selectivity filter collapse (Zhou and MacKinnon, 2003) and speeds BK channel deactivation (Piskorowski and Aldrich, 2006). Second, although binding of the N-terminal inactivation peptide to the intracellular mouth of the pore of many K+ channels follows predictions of open channel block (Demo and Yellen, 1991), the case for BK channels is far more complex and intriguing. Most importantly, binding and unbinding of the BK inactivation particle (provided by the β subunit) does not require channels to visit the fully conducting open state (Benzinger et al., 2006). Third, no evidence has been obtained thus far for trapping of large intracellular blocker molecules like leupeptin, curare, QX-314, and chlorisondamine (Solaro et al., 1997). Fourth, external Rb+ has been shown to destabilize the closed state, presumably due to occupancy of the selectivity filter (slowing deactivation and speeding channel opening; Demo and Yellen, 1992). Finally, a novel conotoxin, K-BtX, appears to stabilize the open state of BK channels from the external side (Fan et al., 2003). It is thus not unreasonable to consider the possibility that BK channels are gated at the selectivity filter, particularly in light of the fact that three other channels activated by intracellular ligands, the CNG, KATP, and SK channels, appear to use selectivity filter gates (Flynn and Zagotta, 2001; Bruening-Wright et al., 2002; Proks et al., 2003).
The classic QA blocker experiments that have worked so well for other K+ channels have been less informative for BK channels because of the relatively large size of the BK channel intracellular cavity. In order for a time-dependent component of block to be visible, the kinetics of block must be slow with respect to the intrinsic rate of channel opening, and the kinetics of block by even the largest molecule tested thus far, decyltriethylammonium (C10), have been too fast to resolve (Li and Aldrich, 2004). We sought to find a QA blocker with sufficiently slow kinetics to allow resolution of time-dependent block and to discern whether intracellular block of BK channels requires channel opening. Here we provide evidence that the QA ion bbTBA (N-(4-[benzoyl]benyl)-N,N,N-tributylammonium bromide) can block both closed and open channels. We conclude that access of bbTBA to the BK channel pore does not appear to be strictly regulated by an intracellular gating structure, and speculate that permeant K+ ions may be gated in the selectivity filter itself.
Materials And Methods
All experiments were conducted using the mSlo gene encoding the mouse BK channel (Butler et al., 1993). Freshly isolated Xenopus oocytes were injected with ∼5 ng or ∼0.05 ng cRNA for macroscopic or microscopic experiments, respectively. Patch clamp recordings were made 1–4 d after injection.
All recordings were made from oocyte macropatches using the inside-out patch configuration (Hamill et al., 1981). Pipettes were pulled from borosilicate glass (VWR Scientific) to yield resistances of 1–3 MΔ, and coated with sticky wax (Kerr Corporation). All experiments were performed at 22°C. Data were acquired with an Axopatch 1B or Axopatch 200A amplifier (Axon Instruments, Inc.), low pass filtered at 10 kHz, and digitized using an ITC-16 converter (Instrutech). Sampling interval was set at 10 or 20 μs using Pulse software (HEKA Electronik) and a Macintosh G3 computer. Data were leak subtracted using a −P/4 protocol stepping from leak holding potential of −120 mV. Macroscopic ionic currents were usually averaged from three or more consecutive sweeps. Idealization of microscopic recordings was accomplished using the TAC software suite (Bruxton Corporation); all further analysis of macrosopic and microscopic data was performed with Igor (WaveMetrics, Inc.).
Analysis and Curve Fitting
Dose–response data were fitted with the Hill equation:
where IbbTBA and Icontrol = current with and without bbTBA, respectively, Kd = apparent dissociation constant for bbTBA, and n = Hill coefficient. The voltage dependence of kOFF, kON, or Kd was determined by fitting with the Woodhull equation (Woodhull, 1973):
where V = test potential, k(0) = the rate at 0 mV, z = valence of bbTBA (+1), and δ = effective electrical distance; together, zδ roughly approximates the voltage dependence of blockade. Voltage dependence of block was also measured by fitting fractional unblock as a function of voltage with a Woodhull equation of the form
Conductance–voltage relationships were calculated from isochronal tail current measurements at −80 mV after variable test steps, and normalized to maximal tail current amplitude. This procedure avoids complications induced by fast Ca2+ block of steady-state currents, and thus accurately reflects relative macroscopic open probability (Po) as a function of voltage (Cox, Cui and Aldrich, 1997). The data were fitted with a Boltzman function:
where G = conductance, V1/2 = potential for half-maximal activation, and z = apparent equivalent gating charge. Current traces depicted in the figures were zero subtracted for display purposes. In general, the SEM is given with averaged values, while the 95% confidence interval is provided for fitted parameters.
Predictions of Percent Block (Fig. 3 C)
The apparent Kds (Kapp) for bbTBA or tetrabutylammonium (TBA) were estimated from the average fractional block (F) observed for each blocker alone using FQA = [QA]/[QA] + KappQA. This yielded KappbbTBA = 5.9 μM and KappTBA = 0.628 mM. The KappbbTBA in the presence of TBA was then calculated: KappbbTBA/TBA = KappbbTBA(1 + [TBA]/KappTBA = 15.3 μM; the fraction of current blocked by bbTBA in the presence of TBA is thus FbbTBA/TBA = [bbTBA]/[bbTBA]+KappbbTBA/TBA = 0.25. Similarly, the KappTBA/bbTBA = 1.16 mM, and the FTBA/bbTBA = 0.46. Together, both blockers block 0.25 + 0.46 = 71% of total current.
Percent block was calculated by summing the fraction of current blocked by TBA, FTBA = 0.62, and the fraction of current remaining in the presence of TBA that is blocked by bbTBA, FbbTBA/TBA = (1 − FTBA)FbbTBA = 0.18, yielding 80% block.
Final Model for bbTBA Blockade
A final model for state-independent blockade of BK channels by bbTBA was constructed (Scheme 1), which accurately reproduced steady-state and kinetic features of block (see Figs. 1, 3, 5, 8, 10, and 11).
This model, based on the simulation protocol developed by Horrigan et al. (1999), consisted of four connected states, with parameters taken from our average macroscopic and microscopic data (see Table I). The final zδ value of 0.15 was calculated from the average values of zδ derived from six different measurements, and closely matches that observed for C10 and TBA (∼0.16 and ∼0.21; Li and Aldrich 2004). Kd(0 mV) was allowed to vary somewhat (8–10 μM) in order to compensate for the observed shift in Kd over time and variability in individual patches. Because microscopic measurements were made earlier in time, the average on and off rates are biased toward the lower Kd(0) observed at that time (∼5 μM; Table I). Assuming that the decline in apparent affinity of bbTBA over time is the result of a degradative process, we expect that the increased Kd reflects a decrease in kON, with relatively little change in kOFF. Thus, in order for our final model to accommodate measurements made later in time, we constrained the model to the macroscopic steady-state parameters Kd(0) = 8–10 μM and zδ = 0.15, and also required that kOFF (0) = 562 s−1, while kON(0) was scaled to match the Kd (kON = 0.625 × 108 M−1s−1). Microscopic measurements of the voltage dependence of kON(0.05) and kOFF(0.7) were assumed to be accurate, and kOFF was adjusted slightly (to 0.1) to yield zδtotal = 0.15. Note that without these adjustments, the kinetic features of block, including the triphasic time dependence observed at high concentrations, were still reproducible. Our final model makes two assumptions: (1) bbTBA binding does not alter the intrinsic rates of channel activation and deactivation (e.g., the same values for α and β were used for the closed/open transition for both unblocked and blocked channels), and (2) intrinsic affinity of bbTBA for the open and closed states is equal (which our data suggest is not entirely accurate; see Figs. 8–10,9).
The pipette (external) solution contained (in mM) 140 KMeSO3, 20 HEPES, 2 KCl, and 2 Mg2+ (pH 7.2). The bath (internal) solution contained (in mM) 136 KMeSO3, 20 HEPES, 6 KCl, CaCl2, and HEDTA. The amount of Ca2+ and HEDTA added to the internal solution was varied by estimating the amount of CaCl2 and buffer required to yield a given final Ca2+ concentration with MAX CHELATOR software (Bers, 1982). Thus, solutions of 0–4 μM free Ca2+ were buffered with 2 mM HEDTA, and those containing 10–20 μM free Ca2+ contained 1 mM HEDTA, while 100 μM free Ca2+ solutions were not buffered. To minimize voltage-dependent block by contaminant barium, (+) −18-crown-6-tetracarbolxylic acid (18C6TA) was added to the internal solution just before recording at a final concentration of 40 μM (Cox et al., 1997). Both 18C6TA and bbTBA were prepared as millimolar stocks in internal solution and stored at −20°C before diluting into fresh solutions immediately before beginning experiments. Each patch was constantly bathed with internal solution using a sewer pipe system (DAD-12; Adams and List Assoc., Ltd.), which achieves solution changes in <1 s.
N-(4-[benzoyl]benyl)-N,N,N-tributylammonium bromide (bbTBA) was acquired from Spectra Group Limited, Inc. All experiments were performed using dilutions of an initial 50 mM stock prepared from newly obtained salt, which was stored frozen in small aliquots at −20°C. We did notice a slow decline in efficacy of block over a period of ∼9 mo. This apparent degradation of bbTBA was independent of storage in solution (50 mM stocks frozen at −20°C) or as a salt, based on comparison of Kd at +100 mV of (a) the original stock solution, 3.2 μM, (b) the same stock solution 9 mo later, 4.4 μM, and (c) a freshly made solution prepared from 9-mo-old salt, 4.4 μM. The latter measurements (b and c) were made from Hill fits of dose–response data obtained from single patches, n = 3 patches. The absorption spectra (measured after 9 mo) for the latter two solutions were also comparable (Amax = 0.20 or 0.18, respectively); no measurement had been made for the original stock solution 9 mo prior. This slow decline in efficacy of block adds some variability to our averages, but will otherwise not affect the conclusions of this work, as the fractional block experiments were all performed in an ∼3-mo period, and all comparisons were made to control experiments performed concurrently. Tetrabutylammonium (TBA) was purchased from Alfa Aesar. All other chemicals were purchased from Sigma-Aldrich.
Macroscopic bbTBA Blockade of BK Channels Is Time Dependent
To characterize blockade of BK channels by bbTBA, we first performed dose–response experiments, where macroscopic K+ currents were elicited from inside-out macropatches using a series of voltage steps in the absence and presence of varying concentrations of bbTBA (Fig. 1). Fig. 1 A illustrates blockade of BK currents measured at +140 mV. bbTBA block was dose dependent and fully reversible (see Figs. 3 and 9). In Fig. 1 B, average fractional block, calculated from steady-state ionic currents in the presence and absence of bbTBA, was plotted as a function of bbTBA concentration for each of six test voltages in order to determine the dosage for half-maximal block (Kd). Hill fits of these data yielded Kd values in the low micromolar range, with Hill coefficients near 1. The resulting Kd values were then plotted as a function of test voltage, for voltages ≥+40 mV (Fig. 1 C), where open probability (Po) in 100 μM Ca2+ is expected to saturate. As expected for a charged molecule (see structure in Fig. 1 B), block was modestly voltage dependent, with a Kd at 0 mV of 9.5 μM, and a total voltage dependence zδ of 0.15, similar to other QA blockers in BK channels (Li and Aldrich, 2004). The lesser degree of voltage dependence compared with QA block of Kv channels (zδ ∼ 0.4) may reflect the larger diameter, and presumed lower resistance, of the vestibule in BK channels. Fig. 1 demonstrates that bbTBA behaves as a typical QA blocker, where dose-dependent blockade results from binding of bbTBA to a single voltage-dependent binding site.
A fast time-dependent component of blockade was visible under certain recording conditions where activation is sufficiently fast to allow resolution of block kinetics (see Fig. 1 A). This component of block is examined in more detail in Fig. 2. Representative current families elicited by test potentials ranging from +20 to +120 mV in the absence (Fig. 2 A) and presence (Fig. 2 B) of 5 μM bbTBA illustrate the voltage range over which the time-dependent component is visible in 100 μM Ca2+ internal solution. Because BK channel activation is slowed in the absence of Ca2+ (Cox et al., 1997), the lowest concentration of Ca2+ for which time-dependent block was reliably observed (at extreme positive test potentials) was ∼4 μM Ca2+ (see Fig. 7). Fig. 2 C shows, on a faster time scale, time-dependent block at bbTBA concentrations ranging from 1 to 20 μM at +180 mV in 100 μM Ca2+ solutions. Note that the time constant of the current relaxation is well fitted by a single exponential function (dashed lines) and is linearly related to bbTBA concentration in this range (Pearson's R = 0.998; unpublished data), consistent with a bimolecular blocking reaction.
Time-dependent blockade by QA ions often signifies open channel block (Scheme 2); however, it is important to note that the implication is not limited to this possibility.
The approximately twofold change in fractional block (estimated from the peak and steady-state levels of current in bbTBA) during the time-dependent component is well within the range that could be explained by voltage-dependent, state-independent block (Scheme 3) alone (assuming zδ ∼ 0.15, which predicts an approximately threefold change in block for the 260-mV step shown in Fig. 2 C), without a need to invoke state-dependent block.
The bbTBA Binding Site Is Located in the Pore
Before using bbTBA to evaluate accessibility of the pore region of the BK channel, it is important to confirm that the bbTBA binding site is located at the pore. We expect that bbTBA is a pore blocker since block is voltage dependent and bbTBA is structurally analogous to organic cations such as TBA, for which evidence of pore block is strong. To provide direct evidence of pore block, we evaluated the ability of the known pore blocker TBA to compete with bbTBA for binding to the BK channel. The experiment was conducted as depicted in Fig. 3 A####; ionic currents at +120 mV were measured from channels before, during, and after exposure to 5 μM bbTBA, 1 mM TBA, or a solution containing both blockers (5 μM bbTBA + 1 mM TBA). Traces averaged from five consecutive steps (Fig. 3 B) were then used to calculate the percent block observed in each blocker solution. The average percent block from three patches is shown in Fig. 3 C. If both blockers can bind and block channels simultaneously (noncompetitive block; Scheme 4), we expect the degree of block by bbTBA + TBA to be 80%; however, we observe that block is weaker (71%), exactly what one can predict from a model of competitive block (Scheme 5) where bbTBA and TBA compete for a single binding site.
The data in Fig. 3 B were compared with the predictions of our final model (described in Materials and methods and Table I) to assess its adequacy in describing steady-state and kinetic properties of bbTBA blockade, and to illustrate the predictions made in Fig. 3 C. Blockade of ionic current by bbTBA was well described by our final model, as were currents in TBA (using estimates of block kinetics from Li and Aldrich ; solid lines). The degree of blockade observed in the presence of both blockers matched predictions of the competitive model (Scheme 5, solid line), but was not compatible with noncompetitive block (Scheme 4, dashed line). Note that the virtual absence of time dependence can be accounted for by both models. These data, along with the bimolecular kinetics and voltage dependence, provide good evidence that bbTBA is a pore blocker, as expected.
To evaluate more carefully the kinetics and voltage dependence of bbTBA blockade, we performed microscopic experiments (Fig. 4). Fig. 4 A illustrates representative single channel records from inside-out patches at +100 mV in 100 μM Ca2+ internal solution, in the absence and presence of either 1 or 3 μM bbTBA. While the BK channel activity is characterized by high open probability (Po) and short flickery closings, Po decreases and open times become progressively shorter as the concentration of bbTBA is increased. Dwell time histograms constructed from idealized records were fitted with single exponential functions in order to yield the mean open and shut dwell times (Fig. 4 B). Block rates calculated from these mean dwell times were then plotted as a function of bbTBA concentration in order to yield the second order kON at +100 mV, the average kOFF, and the resultant Kd (Fig. 4 C). As an alternative to the kinetic measurements, we also determined the microscopic Kd using steady-state measurements of open probability in various bbTBA concentrations (Fig. 4 D). The results from Hill fits to normalized Po plotted as a function of bbTBA concentration were roughly similar to those of the kinetic measurements in Fig. 4 C.
The measurements illustrated in Fig. 4 were performed over a range of test potentials in order to yield the average microscopic block parameters shown in Fig. 5. Values of kON (Fig. 5 A), kOFF (Fig. 5 B), and Kd (Fig. 5 C) derived from fits of these data to a Woodhull model are shown in Table I. The general agreement between block parameters derived from steady-state and kinetic measurements validates our final model of bbTBA block, which assumes a single bbTBA binding site and single step-blocking reaction, and adequate correction for native channel behavior in the calculation of kON. Indeed, despite the discrepancy in Kd (discussed in Materials and methods), our final model describes the data in Fig. 5 C quite well (solid lines). Furthermore, the voltage dependence of bbTBA blockade from the microscopic and macroscopic measurements also match closely (Table I), confirming that the macroscopic measurements are not contaminated by leak or nonspecific ionic currents. Finally, we performed kinetic measurements of bbTBA block in 2 μM Ca2+ as well (unpublished data) and confirmed that bbTBA block does not appear to depend on Ca2+ intrinsically. Together, these single channel measurements give us a more reliable idea of the kinetics and voltage dependence of bbTBA blockade, which are important in considering whether time-dependent block can be accounted for by voltage-dependent block alone.
bbTBA Does Not Behave as a Classical Open Channel Blocker for BK Channels
If BK channels posses an intracellular gate that occludes the bbTBA binding site when closed, one might predict that blocker binding and gate closing are mutually exclusive events, and that bound blocker would prevent or hinder closing of the activation gate. This is exactly what Armstrong originally observed for squid K+ channels (Armstrong, 1971), and is commonly observed for state-dependent blockers. In Fig. 6, we ask whether bbTBA alters the rate of channel closing for BK channels.
Tail currents measured from channels in the presence of 5 μM bbTBA were well fitted by single exponential functions (Fig. 6 A), although fits did deviate slightly from single exponential kinetics at potentials positive to −80 mV (discussed below). In contrast to expectations for state-dependent block, where bbTBA is predicted to slow the macroscopic rate of deactivation by ∼2.3-fold, Fig. 6 B illustrates that bbTBA had no obvious effect on deactivation time constants over a wide range of potentials, although a progressive, mild slowing is apparent at positive potentials above −120 mV (1.1–1.4-fold). Although initially counterintuitive, our final state-independent model does in fact predict a progressive 1.2–1.3-fold slowing of deactivation at potentials between −200 and −50 mV, caused by the tendency for some open channels (either blocked or unblocked) to either unblock (OB-O-C) or block (O-OB-CB) before closing, rather than close directly in a single step (OB-CB or O-C). Thus, the kinetics of bbTBA block are such that although unblock does not require opening (in this model), a significant fraction of channels tend to traverse the same two-step closing pathways that state-dependent block would require, especially at positive potentials. TBA is known to speed deactivation (Li and Aldrich, 2004), an effect that was recapitulated in a single patch shown in Fig. 6 C, illustrating that the approximately twofold speeding of deactivation by TBA reported previously was easily distinguishable in our measurements. The results of Fig. 6 suggest that closure of the activation gate is not physically hindered by the presence of bbTBA, which is consistent with at least two alternative interpretations: either closure of the gate does not occlude the bbTBA binding site (state-independent block) or the bbTBA binding site is gated, but channels close with a blocker molecule trapped inside the pore (an alternative model of state-dependent blockade).
Some channels are capable of trapping small molecules inside the pore cavity after closure of the activation gate (Armstrong, 1971; Strichartz, 1973; Neely and Lingle, 1986; Holmgren, Smith and Yellen, 1997). This would be represented by extension of the state- dependent model (Scheme 2) into C-O-OB-CB. Modeling confirms that this trapping model could also explain the mild slowing of deactivation by bbTBA and still assume that block is state dependent. Indeed, there is evidence that the internal cavity of BK channels is larger than that of most K+ channels (Li and Aldrich, 2004), and thus may be large enough to accommodate bbTBA with a closed (cytoplasmic) gate. If BK channels can trap bbTBA, then one would expect that a substantial fraction of open channels that are blocked during the test pulse will remain blocked immediately after rapid closure and reopening, and current will simply rise to steady-state block levels with minimal or no time dependence. However, we found no evidence for trapping of bbTBA inside BK channels (Fig. 7). Fig. 7 A shows a trapping experiment performed in 20 μM internal Ca2+, where three control pulses to +150 mV were followed by three pulses in the presence of 5 μM bbTBA. The interpulse interval was set to 5 ms at −200 mV, in order to close channels rapidly and minimize the possibility of bbTBA escaping from open channels, either before closing or during rare reopening events (Po at −200 mV ∼ 0.0001; Horrigan and Aldrich, 2002). Fig. 7 A shows that time-dependent block remained intact for each test pulse. In the context of a state-dependent model of block, the persistence of the time dependence would suggest that channels do not remain blocked during the interpulse interval, i.e., channels cannot trap bbTBA and have to wait for blocker to exit before closing (in which case one would expect to see slowing of deactivation). Obviously, this interpretation is difficult to reconcile with Fig. 6. Alternatively, a state-independent model of block posits that time dependence simply reflects equilibration of channels to the change in test potential; one would not expect this to change for three identical test pulses.
The trapping experiment was repeated in 4 μM internal Ca2+ (Fig. 7 B) where probability of channel opening at −200 mV is expected to be 10-fold lower than in 20 μM Ca2+ (∼0.00001; Horrigan and Aldrich, 2002). The apparent inability of channels to trap bbTBA (and the lack of slowed deactivation at negative potentials in Fig. 6) cannot be explained by rapid escape of bbTBA before channels close, since the expected rate of unblock at −200 mV (1234 s−1; extrapolated from Fig. 5 B) is roughly four times slower than the rate of channel closing in Fig. 7 B (4680 s−1). This was confirmed by demonstrating that block of tail currents at negative potentials was much greater than predicted by extrapolation of steady-state block at positive potentials (Fig. 7 C), indicating that channels are still blocked at the time of closing. Note further that deactivation in 4 μM internal Ca2+ (Fig. 7 B) is even faster than the measurements made in Fig. 7 C, using 100 μM Ca2+ (Cui et al., 1997). The recovery from block of channels between test pulses in Fig. 7 (A and B) is consistent with the expected rate of unblock at −200 mV noted above, and, given the exceedingly low probability of channel opening, implies that bbTBA dissociates from closed channels. Together, the results of Figs. 6 and 7 are most easily understood in terms of a state-independent model of blockade, where the gate is located in a position that does not interfere with the deactivation process and does not permit trapping of intracellular bbTBA inside the cavity.
bbTBA Blockade Does Not Depend Strongly upon Open Probability
One can directly assay the state dependence of channel blockade by comparing the efficacy of block under conditions where Po is high and low, while voltage is kept constant. The BK channel is particularly amenable to such experiments by virtue of its activation by both voltage and Ca2+. Thus, one can measure ionic current (and blockade) at a single test potential while activating channels to different levels of open probability with Ca2+. If bbTBA block requires channel opening, steady-state block should be maximal at high Po and minimal at low Po where only a few channels have reached the open state. Conversely, block that is independent of channel conformation should be equally effective at all levels of Po. To discriminate between these possibilities, macroscopic steady-state fractional block was measured at several test potentials under a variety of different Ca2+ conditions representing a wide range of open probabilities (Fig. 8). Fig. 8 A illustrates two representative measurements of ionic current blockade at +90 mV from the same patch exposed to either 100 or 1 μM internal Ca2+, ±5 μM bbTBA. The relative Po vs. voltage relationships for this patch (Fig. 8 B) show that open probability at this voltage is ∼1 or ∼0.25 under these conditions, respectively. It is apparent in Fig. 8 A that despite the large difference in open probability, fractional block appears roughly equivalent (although block at low Po is ∼5% less effective). When fractional unblock is plotted as a function of test potential together with relative Po for both Ca2+ conditions (Fig. 8 B), it is evident that fractional block changes modestly with test voltage, in accordance with a Woodhull model of voltage-dependent block, even between +90 and +180 mV where Po changes dramatically. The voltage dependence of blockade is characterized further in Fig. 8 C.
Measurements of the kind illustrated in Fig. 8 B were made under a variety of different Ca2+ conditions ranging from 1 to 100 μM Ca2+, and fractional unblock was then plotted as a function of test potential (Fig. 8 C) or relative Po (Fig. 8 D). Fig. 8 C illustrates the average voltage dependence of block at low, medium, and high Po. The apparent affinity of block (derived from fitting of each dataset to Eq. 3) at low Po was slightly lower (∼10 μM) than at high Po (∼8 μM), although our final model of bbTBA block (Kd(0) = 8.9 μM and zδ = 0.15, solid line) describes all three datasets quite well. The similarity of the voltage dependence of block at high and low Po (Fig. 8 C: high Po, zδ = 0.14; low Po, zδ = 0.17) suggests that block of open and closed states occurs at the same site across the transmembrane field; however the degree of voltage dependence could also be affected by ion occupancy. In Fig. 8 D, the same measurements of fractional unblock were plotted individually as a function of open probability, for each of three different voltages. Individual data points at each voltage were then fit with a linear function. No obvious dependence of block on Po is apparent from these data. Fig. 8 suggests that bbTBA can bind to either closed or open BK channels, with no strong preference for the open conformation.
Because the macroscopic measurements shown in Fig. 8 are subject to limitations such as voltage-dependent Ba2+ block (Cox et al., 1997) and difficulty in measuring small current amplitudes at low Po, as well as potential errors in estimating open probability from macroscopic Po-V curves, we performed similar experiments using microscopic measurements of fractional block (Fig. 9). The time course of a typical experiment is shown in Fig. 9 A. 10-s sweeps at +40 mV were recorded from an excised patch perfused with 2 μM internal Ca2+, which was alternated with 5 μM bbTBA several times before repeating the process with 4, 10, and 100 μM Ca2+ solutions ± 5 μM bbTBA. This procedure allowed us to monitor directly the stability of Po over the duration of the experiment and ensure accurate measurement of Po changes induced by blocker. The first two sweeps in the absence and presence of bbTBA for 2 μM Ca2+ (sweeps 1 and 2 in Fig. 9 A) are illustrated in Fig. 9 B. As shown in the top left panel of Fig. 9 B, Po in 2 μM Ca2+ was 0.09; 100 ms of this record shown on an expanded time scale below illustrates that channel activity mainly consisted of short openings punctuated by flickery closing events. When 5 μM bbTBA was applied, Po was reduced to 0.05, primarily due to a reduction in mean open time (Fig. 9 B, right). Fractional block at Po ∼ 0.09 computed from the raw amplitude histograms in Fig. 9 C was 45%.
Sweeps 3 and 4 recorded in 100 μM internal Ca2+ ± 5 μM bbTBA are illustrated in Fig. 9 D. In this case, control Po (0.89) was reduced to 0.44 in the presence of bbTBA. Note from the data on a faster time scale that channel activity is qualitatively similar to that observed in 2 μM Ca2+; control openings are still interrupted by flickery closings, and mean open time in blocker is similarly reduced. Fractional block computed from the amplitude histograms in Fig. 9 E at Po ∼ 0.9 was 51%. This experiment shows that, as the macroscopic data suggested, block does appear modestly reduced at low Po. To see whether this ∼6% decrease in block at low Po is compatible with state-dependent block, we compared data from several patches to model predictions in Fig. 10.
Individual measurements of microscopic fractional block from several experiments were plotted as a function of Po (Fig. 10). Predictions of fractional block made using state-dependent (Scheme 2) and state-independent (Scheme 3) models are shown as dashed lines. The data are best fit by a line that matches closely the predictions of a state-independent model of block. However, it is apparent from these data that, on average, a 7% decline in block is observed over the range of Po from 0.97 to 0.03. This translates into a modest 1.3-fold change in apparent affinity at +40 mV, from 4.5 to 6 μM. Although this change in apparent affinity deviates from a model of purely state-independent block, it is clearly not compatible with the large change in apparent affinity predicted by the state-dependent model for the same change in Po (which predicts a 31-fold lowering of apparent affinity, giving a Kd of 144 μM at Po = 0.03).
The mechanistic picture of bbTBA blockade that emerges from Figs. 8–10,9 is one where binding of bbTBA is enhanced to a small degree in the open state, but not gated in an absolute fashion by channel conformation. The difference in apparent affinity between the open and closed conformations results in a modest 7% change in the fraction of blocked channels over a range of Po encompassing nearly two orders of magnitude. The final piece of evidence that we present against open channel block comes from returning again to the initial observation that bbTBA blockade is time dependent.
At high concentrations of bbTBA (∼75–100 μM), a third component of time-dependent block becomes apparent in the macroscopic record. Fig. 11 A illustrates the three kinetic phases of time-dependent block for currents at +90 or +150 mV in the presence of 75 μM bbTBA; current initially increases to a maximum, and then rapidly declines, and finally increases again to steady state. These complex kinetics are predicted by our final model of state-independent block, and, importantly, are not compatible with state-dependent block (Fig. 11 B). The data at +90 mV were roughly approximated by our final model (dashed black line). Only minor adjustment of those parameters was required to fit the steady-state level of blockade in this patch; we also added the assumption that bbTBA slows activation in order to more accurately describe the kinetics (solid red line). (It is important to note that, while we have no evidence to justify this particular adjustment, this slowing of activation was not required in order to reproduce the triphasic kinetics of block.) The data require that some fraction of channels start in the closed, blocked state, consistent with expectations of voltage-dependent, state-independent block. The smaller fraction of closed, blocked channels required to fit our data (35–50%, compared with predictions of ∼80%) may be due to the preference of bbTBA for the open state. Finally, the state-dependent version of our final adjusted model is grossly incompatible with the observed data, predicting a much larger change in fractional block during the time-dependent component, and biphasic kinetics (solid blue line).
Fig. 11 C uses the occupancy of each state over time to explain how these complex kinetics arise in the ionic currents, and why they cannot be observed for state-dependent block. The final adjusted model is illustrated in the inset, including the rate constants used for modeling block at +90 mV in 75 μM bbTBA. Each state in the model is represented by a different color, and the fraction of channels populating that state is plotted over time in the same color. As in Fig. 11 B, we assume that 35% of channels are initially closed and blocked (CB), while the remaining 65% are closed (C). At the start of the voltage step, occupancy of the closed state falls as many channels quickly become blocked (C-CB), and the remaining channels more slowly enter the open state (C-O). Thus, occupancy of CB rises and then begins to decline slowly as channels continue on to the open blocked state (CB-OB); this in turn is reflected by the steady rise in occupancy of OB. Meanwhile, those channels that entered the open state initially rapidly become blocked and contribute to the population of OB. Finally, the slow transition of channels from C to O comes to a halt, and occupancy reaches a steady state. Thus, occupancy of the open state (which shapes the kinetics of current blockade) is dominated in the first phase by channels entering the open state (C-O); during the second phase, channels leave the open state to become blocked; and finally in the third phase, the slow kinetics of the closed to open transition again dominate as occupancy of O and OB reaches steady-state equilibrium.
Triphasic block kinetics arise from the simultaneous combination of the relatively fast blocking process and the comparably slower activation process. The rapid population of the blocked states requires direct connectivity between the closed and blocked states (C and CB). This is in direct contrast with state-dependent models, where channel opening is the rate-limiting step for blockade, and explains why triphasic block kinetics cannot be obtained from state-dependent models. Although triphasic kinetics do not require some fraction of channels starting in CB, the correct amplitude of the peak current does. No permutation of state-dependent models, including a trapping model, was able to account for this complex time dependence. Fig. 11 provides strong support for a state-independent mechanism of bbTBA blockade of BK channels.
Time-dependent blockade, although often suggestive of open channel block, is equally compatible with models that do not require the open channel conformation for blocker to access its binding site. We found that the time dependence of bbTBA block is a direct result of the intrinsic voltage dependence of block, with little dependence on channel opening. The present work provides the following evidence of state independent block. (a) Macroscopic deactivation was not significantly affected by bbTBA at negative potentials, and only mildly slowed at positive potentials; this is in direct contrast to expectations for state-dependent block, which predicts slowing at all potentials, especially negative voltages, and suggests that bound bbTBA does not prevent channel closing. (b) Trapping experiments suggest that bbTBA cannot be confined within the cavity after channel closure, indicating that closure of the activation gate does not prohibit bbTBA dissociation. (c) Macroscopic and microscopic measurements of steady-state block as a function of open probability suggest that block does in fact have a mild preference for the open channel conformation, but to a degree that deviates substantially from predictions of a purely state-dependent model. The observation that ∼45% of channels are blocked when fewer than 10% of channels are open further supports the idea that channels can be blocked in the closed state. (d) Finally, our final state-independent model of bbTBA block was able to accurately describe the kinetics of blockade at high concentrations of bbTBA, where the time dependence of block was triphasic. These unusual triphasic kinetics, which result from simultaneous blocking and opening of channels, could not be reproduced by any state-dependent model we tested. Together, our results strongly suggest that closure of the activation gate does not present a significant barrier to bbTBA binding site access.
Implications for Location of the BK Channel Permeation Gate
Relatively unrestricted access of intracellular bbTBA in the open and closed channel conformations implies that K+ flux is gated at a location that is extracellular to the bbTBA binding site. Since QA ions are known to dock at the internal entrance to the selectivity filter, deep within the channel's intracellular cavity (Choi et al., 1993; Zhou et al., 2001a; Lenaeus et al., 2005), this places the permeation gate at or above the selectivity filter. Our findings suggest that current assumptions about BK channel architecture and gating may need to be reevaluated.
Jiang et al. (2001, 2002a,b) postulate that mammalian BK channels share the same gating mechanism proposed for the bacterial Ca2+-gated K+ channel MthK. Extrapolation of their gating model from the prokaryotic MthK to the eukaryotic BK channel was based on <20% amino acid conservation between the two proteins, and a shift in the voltage dependence of activation after mutation of a putative salt bridge site in the BK channel (Jiang et al., 2001). Their gating model involves an intracellular gating ring which, upon ligand binding, tugs the S6 bundle helices apart. Niu et al. (2004) have expanded upon this gating model by proposing that both Ca2+ and voltage activation of BK channels involves passive coupling of either the ligand-binding gating ring or the voltage-sensing apparatus to the S6 gate via separate hinges, based on the observation that deletions and insertions into the lower S6 region alter gating. However, mutations in the S6 regions of KATP and CNG channels, which do not appear to use intracellular S6 gates, also disrupt gating (Trapp et al., 1998; Flynn and Zagotta, 2003). In addition, Ca2+-dependent gating of the MthK channel differs in significant ways from that of the BK channel, suggesting the possibility of different gating mechanisms (Zadek and Nimigean, 2006). Thus, in the absence of any direct evidence in support of an intracellular gate (or gating ring), the degree to which this model represents BK channel gating remains unclear.
Recent work (unpublished data) provides strong evidence that conformational rearrangements associated with gating regulate access of intracellular Shaker ball peptides to the BK channel pore. The support for this comes from time-dependent block, slowing of macroscopic deactivation and a hook in the tail currents, and correlation of the degree of fractional block with macroscopic Po. Although clearly demonstrating state-dependent block by ball peptides, that study is not conclusive about the gate for K+ ions, since these 20 amino acid peptides are orders of magnitude larger than permeant K+ ions. How can we reconcile gated access of a bulky peptide with relatively unrestricted access of a smaller quaternary ammonium derivative? One could speculate that permeation gating of BK channels involves movements at or above the selectivity filter that are obligatorily coupled to rearrangements of an intracellular protein structure (for example, the putative gating ring). Reorganization of intracellular regions, including the S6 helices, may indeed gate access of ball peptide to the pore, while for the smaller blocker bbTBA, these same movements only slightly alter access and/or binding to a constantly available site. The reverse conformational change in S6, while permitted by bound bbTBA, could be sterically prohibited until ball peptide dissociates, resulting in slower macroscopic deactivation.
Might BK Channels Have a Selectivity Filter Gate?
Gating of K+ flux at the selectivity filter of BK channels is an attractive possibility, when one considers the classic experiments of Miller et al. (1987), which showed not only that occupancy of BK channel selectivity filter by Ba2+ destabilizes the closed state (Miller, 1987; Miller et al., 1987), but also that, in the closed state, the selectivity filter is inaccessible to Ba2+ from both the internal and external face of the membrane, as if both ends of the selectivity filter pinch shut (Miller et al., 1987). This possibility is made all the more attractive by recent work illustrating that BK channel gating is unusually sensitive to selectivity filter occupancy; permeant thallium favors flickery gating on the submillisecond time scale, which may be a direct result of rapid conformational changes of the selectivity filter (Piskorowski and Aldrich, 2006), and also prevents TBA from speeding deactivation (Piskorowski, 2005).
One of the best studied selectivity filter gates is the C-type inactivation gate of Kv channels. Functional evidence strongly suggests that C-type inactivation involves conformational changes of the outer mouth and selectivity filter, as well as the S4 voltage-sensing region, of Kv channels (Timpe et al., 1988; Choi, Aldrich and Yellen, 1991; Lopez-Barneo et al., 1993; Yellen et al., 1994; Baukrowitz and Yellen, 1996; Loots and Isacoff, 1998). A similar mechanism may account for voltage-dependent gating of Kcsa (Cordero-Morales et al., 2006). A rather striking feature of BK channels, by contrast, is their lack of intrinsic inactivation, even after long periods of recording. A simple explanation for this might be that the C-type inactivation gate for BK channels is in fact the permeation gate itself.
Comparisons with CNG Channel Gating
Gating of CNG channels, which are also activated by intracellular ligand and are evolutionary relatives of BK channels, embodies many of the basic principles we propose for BK channels. In general, activation appears to involve a substantial rearrangement of the lower S6 bundle crossing, as evidenced by state-dependent block by intracellular tetracaine (Fodor et al., 1997a,b), dequalinium (Rosenbaum et al., 2004), and state-dependent QA block and accessibility (Contreras and Holmgren, 2006). These motions, which widen the constriction from ∼5 to 12 Å (Flynn and Zagotta, 2001; Flynn and Zagotta, 2003), are required to accommodate both spontaneous salt bridge formation at lower S6 residues only in the closed state, and the large difference in reactivity of the relatively large MTSET molecule with pore-lining residues in the open and closed states (Flynn and Zagotta, 2001). However, constriction of the internal pore entrance does not constitute permeation gating for CNG channels: smaller Ag+ and MTSEA ions react roughly equivalently with pore residues regardless of channel conformation (Flynn and Zagotta, 2001), and pore block by Ni2+ is state independent (Karpen et al., 1993). Rather, it is thought that S6 conformational changes are coupled to intrinsic permeation gating of CNG channels at the selectivity filter, a proposal that draws support from state-dependent conformational changes observed in the outer pore region (Becchetti and Roncaglia, 2000; Liu and Siegelbaum, 2000).
CNG channels also exemplify that the apparent affinity of blockers can depend, often dramatically, on channel conformation, without being gated, e.g., enhanced coordination of Ni2+ by lower S6 histidine residues in the open state (Gordon and Zagotta, 1995a,b), and increased apparent affinity of several large blockers in the closed state (Fodor et al., 1997a,b; Rosenbaum et al., 2004). Similarly, block of BK channels by C10 showed a preference for the open state, which was far less than one would expect based on pure open channel block (Li and Aldrich, 2004). The observed slight preference of bbTBA for the open conformation in the present study appears similar.
Widening of the intracellular bundle crossing for both CNG and BK channels gates access of larger molecules, while smaller reagents and permeant ions are able to pass through the aperture regardless of channel state. If the largest molecules able to pass the bundle crossing freely are any indication of size, the diameter of the closed state bundle crossing is larger for BK channels than for CNG channels, consistent with the larger dimensions of the BK channel cavity (Li and Aldrich, 2004). We speculate that BK channels may share with their CNG cousins a conserved gating mechanism, which involves conformational changes at the bundle crossing and permeation gating at the selectivity filter.
Conclusions and Future Directions
Our work presents evidence for a voltage-gated channel that is not gated at the S6 bundle crossing. Although the exact location of the permeation gate for BK channels remains an open question, the most obvious candidate position is the selectivity filter itself. Despite numerous potential pitfalls (including 29 native cysteines, functional consequences of MTS modification, activation of BK channels by Cd2+, and buffering requirements for Ag+, which would stipulate minimal Ca2+ conditions), future experiments using cysteine accessibility and small probes (e.g., Ag+, Cd2+, or MTSEA) may be useful to convincingly establish the location of the permeation gate.
Thanks to Jon Sack for planting the seed that started the story, technical help with programming and data analysis, and unbounded enthusiasm for channel biophysics. Thanks to Weiyan Li for guidance and support at every stage of this project. We also thank Justin Dubois for assistance in locating novel QA molecules for electrophysiological assays.
This work was supported by the Howard Hughes Medical Institute and the Mathers Foundation.
David C. Gadsby served as editor.
Abbreviations used in this paper: bbTBA, N-(4-[benzoyl]benyl)-N,N,N-tributylammonium bromide; BK, large-conductance Ca2+- and voltage-activated potassium; C10, decyltriethylammonium; CNG, cyclic nucleotide-gated; Po, open probability; QA, quaternary ammonium; TBA, tetrabutylammonium.