Abstract
The Ca^{2+}dependent gating mechanism of largeconductance calciumactivated K^{+} (BK) channels from cultured rat skeletal muscle was examined from low (4 μM) to high (1,024 μM) intracellular concentrations of calcium (Ca^{2+}_{i}) using singlechannel recording. Open probability (P_{o}) increased with increasing Ca^{2+}_{i} (K_{0.5} 11.2 ± 0.3 μM at +30 mV, Hill coefficient of 3.5 ± 0.3), reaching a maximum of ∼0.97 for Ca^{2+}_{i} ∼ 100 μM. Increasing Ca^{2+}_{i} further to 1,024 μM had little additional effect on either P_{o} or the singlechannel kinetics. The channels gated among at least three to four open and four to five closed states at high levels of Ca^{2+}_{i} (>100 μM), compared with three to four open and five to seven closed states at lower Ca^{2+}_{i}. The ability of kinetic schemes to account for the singlechannel kinetics was examined with simultaneous maximum likelihood fitting of twodimensional (2D) dwelltime distributions obtained from low to high Ca^{2+}_{i}. Kinetic schemes drawn from the 10state MonodWymanChangeux model could not describe the dwelltime distributions from low to high Ca^{2+}_{i}. Kinetic schemes drawn from Eigen's general model for a ligandactivated tetrameric protein could approximate the dwelltime distributions but not the dependency (correlations) between adjacent intervals at high Ca^{2+}_{i}. However, models drawn from a general 50 state twotiered scheme, in which there were 25 closed states on the upper tier and 25 open states on the lower tier, could approximate both the dwelltime distributions and the dependency from low to high Ca^{2+}_{i}. In the twotiered model, the BK channel can open directly from each closed state, and a minimum of five open and five closed states are available for gating at any given Ca^{2+}_{i}. A model that assumed that the apparent Ca^{2+}binding steps can reach a maximum rate at high Ca^{2+}_{i} could also approximate the gating from low to high Ca^{2+}_{i}. The considered models can serve as working hypotheses for the gating of BK channels.
Introduction
Largeconductance Ca^{2+}activated K^{+} (BK)^{1} channels, which are activated by micromolar concentrations of intracellular Ca^{2+} (Ca^{2+}_{i}) and by depolarization (Barrett et al. 1982) are present in a wide variety of tissues (reviewed by Rudy 1988; Marty 1989; McManus 1991; Latorre 1994; Conly, 1996). Activated BK channels reduce membrane excitability by allowing K^{+} efflux through their opened pores, which drives the membrane potential more negative. Hence, BK channels form a link in a negative feedback loop that decreases excitability in response to both increased Ca^{2+}_{i} and depolarization.
Information about the gating mechanism of BK channels has accumulated from kinetic studies on native and heterologously expressed channels. BK channels are homotetramers, formed from four alpha subunits (Shen et al. 1994), with one or more Ca^{2+}binding sites per subunit (Schreiber and Salkoff 1997). BK channels typically display a steep relation between Ca^{2+}_{i} and open probability, P_{o}, with Hill coefficients usually in the range of 2–4 (Barrett et al. 1982; Wong et al. 1982; Golowasch et al. 1986; Oberhauser et al. 1988; Reinhart et al. 1989; McManus and Magleby 1991; Art et al. 1995; DiChiara and Reinhart 1995; Cui et al. 1997). Hill coefficients >1 suggest a cooperative action of Ca^{2+} on P_{o}, and are often associated with allosteric modulation of activity (Colquhoun 1973; Fersht 1985).
Models drawn from the MonodWymanChangeux (MWC) model (Monod et al. 1965) describing the transitions of tetrameric allosteric proteins have often been used to describe the Ca^{2+}dependent gating of BK channels. For example, McManus and Magleby 1991 found that a model with three open and five closed states drawn from the MWC model could account for many of the basic features of the Ca^{2+} dependence of the singlechannel gating kinetics from low to intermediate levels of activity, and Cox et al. 1997b found that the full 10state MWC model (and also an extended 12state model) could account for macroscopic conductance–voltage relations over a wide range of Ca^{2+}_{i}.
Despite these successes, MWCtype models appear too simple to account for certain details of BK channel gating. Cox et al. 1997b found that such a model would not account for the details of the macroscopic activation and deactivation kinetics at large positive or negative potentials, and a number of studies have suggested that closed states in addition to those in the MWC model may contribute to the gating (Wu et al. 1995; Cox et al. 1997b). For example, Rothberg and Magleby 1998 found that gating mechanisms drawn from the MWC model underpredicted the numbers of brief closed intervals (flickers) adjacent to the various open intervals. Singlechannel analysis then suggested that there were additional closed states, either intermediate (in the activation pathway) or secondary (beyond the activation pathway) that contributed to the generation of the flickers.
We now test further whether the above extensions of the MWC model are consistent with the underlying gating mechanism. Our approach is to examine the activity of single BK channels from low (4 μM) to very high (∼1 mM) intracellular concentrations of Ca^{2+} (Ca^{2+}_{i}). By driving the channel toward its fully liganded states, we estimated the numbers of open and closed states involved in the gating at high Ca^{2+}_{i}, and through analysis of twodimensional dwelltime distributions and dependency plots, we determined the minimum number of independent transition pathways connecting the fully liganded open and closed states.
We found that at high Ca^{2+}_{i}, the BK channel gates among three to four open and four to five closed states, with two or more effective transition pathways (gateway states) connecting the open and closed states. Neither the general model of Eigen 1968 nor the 55state extension of Eigen's model by Cox et al. 1997b appear consistent with these findings for the gating at high Ca^{2+}_{i}. Extensions of Eigen's model to a twotiered gating mechanism, in which the upper tier was comprised of closed states and the lower tier was comprised of open states, provided a model that would allow gating in the necessary numbers of open and closed states at high Ca^{2+}_{i}, together with the required multiple transition pathways between open and closed states. In the twotiered model, each of the closed conformations of the channel can undergo a direct (possibly concerted) transition to an open conformation. Models drawn from the general twotiered model could describe the singlechannel kinetics from low to high Ca^{2+}_{i}. Since the general twotiered model (and simple extensions of this model) contain previous models that have been used to account for gating under more limited conditions (McManus and Magleby 1991; Rothberg and Magleby 1998), it can serve to unify descriptions of gating. A simpler model that assumed the Ca^{2+}dependent transition rates approached a maximum in high Ca^{2+}_{i} could also approximate many features of the singlechannel kinetics from low to high Ca^{2+}_{i}. The considered models can serve as working hypotheses to further study the gating mechanisms of BK channels.
Methods
Preparation
Currents flowing through single BK channels in patches of surface membrane excised from primary cultures of rat skeletal muscle (myotubes) were recorded using the patch clamp technique (Hamill et al. 1981). Cultures of rat myotubes were prepared from fetal skeletal muscle as described previously (Barrett et al. 1982; Bello and Magleby 1998). All recordings were made at +30 mV using the excised insideout configuration of the patch clamp technique in which the intracellular surface of the patch was exposed to the bathing solution. Analysis was restricted to patches containing a single BK channel. Singlechannel patches were identified by observing openings to only a single open channel conductance level during several minutes of recording in which the P_{o} was >0.4. Experiments were performed at room temperature (22–24°C). Data are presented as mean ± SD.
Solutions
The solutions bathing both sides of the membrane contained 150 mM KCl and 5 mM TES (Ntris(hydroxymethyl)methyl2aminoethane sulfonate) pH buffer, with the pH of the solutions adjusted to 7.0. Contaminant Ca^{2+}_{i} was determined by atomic absorption spectrometry. The solution at the intracellular side of the membrane also contained added Ca^{2+} (as CaCl_{2}), to bring the Ca^{2+} concentration at the intracellular surface (Ca^{2+}_{i}) to the indicated levels. (The solutions did not contain Ca^{2+} buffers.) No Ca^{2+} was added to the extracellular (pipette) solution. Solutions were changed through the use of a microchamber (Barrett et al. 1982).
Recording and Measuring Interval Durations and Identifying Normal Activity
Singlechannel currents were recorded on a digital data recorder (DC37 kHz; Instrutech Corp.), lowpass filtered with a fourpole Bessel filter to give a final effective filtering of 6–10 kHz (−3 dB), and sampled by computer at a rate of 125–200 kHz. The effective filtering is expressed in terms of dead time, which is the duration of an underlying interval before filtering that would just reach 50% of the singlechannel current amplitude with filtering (Colquhoun and Sigworth 1995). Five channels in which stable singlechannel data were obtained over a wide range of Ca^{2+}_{i} were analyzed in detail. The channels and their dead times were: B04, 16 μs; B06 and B16, 28.5 μs; B12, 22.9 μs, and B14, 17.9 μs. The sampled currents were analyzed using custom programs written in the laboratory. The methods used to set the level of filtering to exclude false events that could arise from noise, measure interval durations with halfamplitude threshold analysis, test for stability, and identify activity in the various modes using stability plots, have been described previously, including the precautions taken to prevent artifacts in the analysis (McManus et al. 1987; McManus and Magleby 1988, McManus and Magleby 1989; Magleby 1992).
Restriction of Analysis to Activity in the Normal Mode
The analysis in the present study was restricted to channel activity in the normal mode, which typically involves ∼96% of the detected intervals (McManus and Magleby 1988). Activity in modes other than normal, including the low activity mode and transitions to isolated long shut intervals, were removed before analysis. The low activity mode is readily distinguished from normal activity over the range of low to high Ca^{2+}_{i}, as detailed in Rothberg et al. 1996. The isolated long shut intervals are also readily identified at high Ca^{2+}_{i}, and appear to contribute an insignificant number of the long shut intervals at lower Ca^{2+}_{i} (Rothberg et al. 1996). The mean frequency and duration of the isolated long shut intervals are Ca^{2+} independent, indicating that isolated long shut intervals do not arise from discrete Ca^{2+} block (Rothberg et al. 1996). In spite of their low frequency, appreciable time can be spent in the isolated long shut intervals at high Ca^{2+}_{i} because most of the other closed intervals are so brief. For higher Ca^{2+}_{i}s that give P_{o}s > 0.8, BK channels can spend >30% of their time in the low activity mode that has a P_{o} of ∼0.001 (Rothberg et al. 1996). Because of the potential effects of the low activity mode and isolated long shut intervals on the singlechannel record, it is essential to identify and exclude these (and other modes), as we have done, when studying normal activity over a wide range of Ca^{2+}_{i}.
Log Binning and Plotting of One and TwoDimensional DwellTime Distributions
Both one (1D) and twodimensional (2D) dwelltime distributions were analyzed. The 1D distributions of open and closed interval durations were logbinned as described previously (McManus et al. 1987) at a resolution of 25 bins per log unit for fitting with mixtures (sums) of exponential components. Details of estimating the numbers of significant exponential components are in McManus and Magleby 1988, including the use of the likelihood ratio test for significance. The 2D distributions of adjacent open–closed interval pairs (Fredkin et al. 1985) were logbinned at a resolution of 10 bins per log unit, as described previously (Magleby and Weiss 1990b; Rothberg et al. 1997). The 2D distributions were fitted with mixtures (sums) of 2D exponential components to estimate the numbers of underlying 2D exponential components, and hence the minimum number of kinetic states, as detailed in Rothberg et al. 1997.
With filtering, detected intervals with durations less than about twice the dead time are narrowed (McManus et al. 1987; Colquhoun and Sigworth 1995). For the plotting of 2D dwelltime distributions, the plotted intervals have not been corrected for this narrowing. For the fitting of kinetic models using 2D dwelltime distributions, the measured durations of these intervals were corrected to their estimated true durations before binning and fitting, using the numerical method in Colquhoun and Sigworth 1995. Even with correction, the data were typically fit starting at ∼1.5 dead times for two reasons: (a) intervals whose underlying (true) durations are less than the dead time can be detected when data are filtered, placing an excess of intervals in the first few bins (Magleby and Weiss 1990a), and (b) intervals that are undetected due to filtering (missed events) can introduce extra (phantom) components with time constants typically less than onehalf the dead time (Roux and Sauve 1985; Blatz and Magleby 1986; Hawkes et al. 1992; Magleby and Weiss 1990a), which can also add excess intervals to the first few bins. Depending on the gating mechanism, noise, and filtering, the excess intervals arising from either phenomena can lead to the detection of excess brief exponential components if the fitting includes intervals less than ∼1.5–2 dead times.
Dependency Plots
Dependency plots were constructed from the 2D dwelltime distributions as detailed in Magleby and Song 1992 and Rothberg and Magleby 1998. In brief, the dependency for each bin of open–closed interval pairs with mean durations t_{O} and t_{C} is: 1 where N_{Obs}(t_{O},t_{C}) is the experimentally observed number of interval pairs in bin (t_{O},t_{C}), and N_{Ind}(t_{O},t_{C}) is the calculated number of interval pairs in bin (t_{O},t_{C}) if adjacent open and closed intervals pair independently (at random). The method of calculating expected frequencies for observations that are independent (contingency tables) is a common statistical procedure (compare Mendenhall et al. 1990). The expected number of interval pairs in bin (t_{O},t_{C}) for independent pairing is: 2 where P(t_{O}) is the probability of an open interval falling in the row of bins with a mean open duration of t_{O}, and P(t_{C}) is the probability of a closed interval falling in the column of bins with a mean closed duration of t_{C}. P(t_{O}) is given by the number of open intervals in row t_{O} divided by the total number of open intervals in all rows, and P(t_{C}) is given by the number of closed intervals in column t_{C} divided by the total number of closed intervals in all columns. Since open and closed intervals are paired, the total number of open intervals is equal to the total number of closed intervals. Each open and closed interval forms two pairs: one with the preceding interval and one with the following interval. Hence, the number of interval pairs in a 2D dwelltime distribution is equal to the number of open plus closed intervals minus one.
Estimating the Most Likely Rate Constants for Kinetic Schemes
The most likely rate constants for the examined kinetic schemes were estimated from the simultaneous fitting of the 2D dwelltime distributions obtained at six different Ca^{2+}_{i} using the iterative maximum likelihood fitting procedure described in Rothberg and Magleby 1998 with Qmatrix methods to calculate the predicted 2D dwelltime distributions (Fredkin et al. 1985; Colquhoun and Hawkes 1995b). Corrections for missed events were applied during the fitting using the method described in Crouzy and Sigworth 1990. Since the gating of native BK channels in rat skeletal muscle appears consistent with microscopic reversibility (Song and Magleby 1994), microscopic reversibility was maintained during the fitting.
The number of underlying 2D exponential components that sum to form a 2D dwelltime distribution is given by the product of the numbers of open and closed states (Fredkin et al. 1985; Rothberg et al. 1997). Thus, a single 2D dwelltime distribution fitted with, for example, three open and five closed 2D exponential components can potentially define 22 parameters: three open and five closed time constants and the volumes of 14 of 15 (3 × 5) underlying 2D components. One of the volumes is not free since the volumes must sum to 1.0. Simultaneous fitting of three to six 2D dwelltime distributions obtained at different Ca^{2+}_{i} further increased the information available to estimate the fixed number of rate constants that defined the models. In spite of all this information, rate constants can remain illdefined in complex models (Fredkin et al. 1985; Bauer et al. 1987), which is why the analysis in this paper has been restricted to models with the fewest numbers of states that can approximate the examined kinetic features of the data. Rather than estimating the uncertainty in the rate constants for any single experiment, the uncertainty was estimated by comparing rate constants obtained from fitting data from different channels.
For a given model and rate constants, the equilibrium occupancies of the states could be calculated as described in Colquhoun and Hawkes 1995b. The frequencies of entry into each state could then be calculated by dividing the equilibrium occupancy of a state by its mean lifetime.
Evaluating Kinetic Schemes
Normalized likelihood ratios (NLR_{1000}) have been used to indicate how well any given kinetic scheme describes the 2D dwelltime distributions when compared with a theoretical best description of the data. Normalization corrects for the differences in numbers of interval pairs among experiments so that comparisons can be made among channels. The normalized likelihood ratio per 1,000 interval pairs is defined as 3 where lnS is the natural logarithm of the maximum likelihood estimate for the observed 2D dwelltime distributions given the kinetic scheme, lnT is the natural logarithm of the likelihood of the theoretical best description of the observed distributions, and n is the total number of fitted interval pairs (events) in the observed dwelltime distributions (McManus and Magleby 1991; Weiss and Magleby 1992; Rothberg and Magleby 1998).
The NLR gives a measure of how well different kinetic schemes describe the distributions, but it cannot be used to directly rank schemes since no penalty is applied for the numbers of free parameters. To overcome this difficulty, models were ranked using an information criteria approach (Akaike 1974; Horn 1987), which has the limitation that the significance level is not known. If 4 then model g is ranked above model f, where m_{g} and m_{f} are the maximum likelihood estimates for models g and f, and k_{g} and k_{f} are the number of free parameters for each scheme. We also compared the Akaike rankings to those of the Schwarz 1978 criterion and found that the Akaike rankings for different schemes were more consistent among channels and also were in better agreement with rankings based on the visual perception of the ability of models to describe major features of the data.
Estimating the Theoretical Best Descriptions of the 2D DwellTime Distributions
To evaluate models, it was useful to have an estimate of the theoretical best descriptions of the dwelltime distributions for comparison to the distributions predicted by the various examined gating mechanisms. For gating consistent with a discrete state Markov process (the rate constants remain constant in time for constant experimental conditions), as appears to be the case for BK channels (McManus and Magleby 1989; Petracchi et al. 1991), the log likelihood for a theoretical best description of the 2D dwelltime distribution could be obtained by fitting the distribution with an uncoupled (generic) scheme (Rothberg and Magleby 1998). Uncoupled schemes have direct transition pathways between every open and closed state, and no transition pathways between open states or between closed states (Kienker 1989). The loglikelihood value for the theoretical best fit for the simultaneous fitting of dwelltime distributions obtained at different Ca^{2+}_{i} was given by the sum of the loglikelihood values obtained by fitting data at each Ca^{2+} separately with the uncoupled scheme.
Using Simulation to Generate Predicted SingleChannel Data
To make comparisons between the observed distributions and those predicted by the kinetic models, simulated singlechannel current records were generated with filtering equivalent to that used in the analysis of the experimental current records and with noise similar to that in the experimental current record. The simulated singlechannel currents were then analyzed in the same way that the experimental currents were in order to obtain the predicted 2D distributions, dependency plots, and numbers of exponential components observed in the predicted dwelltime distributions. The method used to simulate singlechannel currents with filtering and noise is detailed in Magleby and Weiss 1990a. Visual comparison between observed and predicted distributions and dependency plots served to guide the analysis, indicate where models were inadequate, and give a sense for how well the models described the data. The critical assessment of the models was based on quantitative comparisons using maximum likelihood and Akaike rankings (Formula 3 and Formula 4).
Online Supplemental Material
Figure S1 presents the equilibrium occupancy of the states in Scheme VII at low (5.5 μM) and high (1,024 μM) Ca^{2}+_{i}. At low Ca^{2}+_{i}, the channel spent most of its time (91.7%) in closed states C8–C11, considerably less time (7.5%) in open states O1–O6, and even less time (0.8%) in the intermediate closed states C12–C17. At high Ca^{2}+_{i}, the channel spent essentially all (99.4%) of its time in the fully liganded open and closed states, with most (96.6%) of the time spent in open states O1 and O4, considerably less time (2.8%) in closed states C7, C12, and C15. A very small amount of time (0.5%) was spent in open states O2 and O5, with negligible time (0.1%) spent in all of the remaining open and closed states. Thus, at high Ca^{2}+_{i}, Scheme VII gates as though there were one effective transition pathway between open and closed states (C15–O4), inconsistent with the observed significant dependencies at high Ca^{2}+_{i}. Available at http://www.jgp.org/cgi/content/full/114/1/93/DC1
Figure S2 presents the estimated rate constants for the three channels examined in detail for Scheme IV, Scheme V, Scheme VI, Scheme VII, Scheme XII, and Scheme XIII. Available at http://www.jgp.org/cgi/content/full/114/1/93/DC1
Results
The findings presented in this paper are based on a complete singlechannel analysis of data obtained from three different BK channels, each studied at six different Ca^{2+}_{i}, from low to high (4–1,024 μM). In addition to these three channels, two additional channels were analyzed at three Ca^{2+}_{i}, including high Ca^{2+}_{i}, and an additional four channels were also examined, with findings consistent with the channels analyzed in more detail. All of the data, analysis, and figures presented in this paper are restricted to data collected during normal (mode) activity, which typically includes ∼96% of the intervals (McManus and Magleby 1988). All experiments were performed at +30 mV.
Saturation of P_{o} and the Mean Open and Closed Interval Durations at High Ca^{2+}_{i}
Currents flowing through a single BK channel in an insideout patch of membrane excised from a cultured rat skeletal muscle cell are shown in Fig. 1 A. These traces are representative of BK channel gating during normal activity over a range of Ca^{2+}_{i}, and illustrate a small fraction of the ∼25 min of continuous recording from this channel. The current traces are excerpts of data obtained with 5.5, 12.3, 132, and 1,024 μM Ca^{2+}_{i} that gave mean open probabilities (P_{o}) during normal activity of 0.061, 0.50, 0.97, and 0.97, respectively. Thus, while the 2.2fold increase in Ca^{2+}_{i} from 5.5 to 12.3 μM gave an 8.2fold increase in P_{o}, the 7.8fold increase in Ca^{2+}_{i} from 132 to 1,024 μM resulted in no increase in P_{o}, indicating that the channel had reached a maximum level of activation. The saturation in P_{o} at high Ca^{2+}_{i} is readily apparent in the current traces in Fig. 1.
The effect of Ca^{2+} on P_{o} for five different single BK channels during normal activity is illustrated in Fig. 2 A. For these channels, increasing Ca^{2+}_{i} in the 4–20 μM range increased P_{o} from ∼0.02 to 0.85, while further increases in Ca^{2+}_{i} to either 132 or 1,024 μM led to a maximum P_{o} ranging from 0.93 to 0.98. Maximum P_{o}s within this range have been observed previously (Barrett et al. 1982; Moczydlowski and Latorre 1983; Rothberg et al. 1996; Cox et al. 1997a).
Fitting P_{o} versus Ca^{2+}_{i} with the Hill equation for data from these five channels gave a K_{0.5} (the Ca^{2+}_{i} for a P_{o} of 0.5) of 11.1 ± 0.7 μM, with a Hill coefficient of 3.5 ± 0.6, consistent with at least four Ca^{2+}binding steps contributing to maximal channel activation. The mean of the maximal fitted P_{o} was 0.95 ± 0.03, and the fitted line indicated that a Ca^{2+} of ∼40 μM was sufficient to drive the channel to within 1% of the maximum P_{o}. Hill coefficients ranging from 2–4 are a common feature of BK channels (see introduction), and higher slopes have been reported (Golowasch et al. 1986; Munoz et al. 1998), consistent with the channel binding at least four Ca^{2+} to become fully activated.
The P_{o} did not reach 1.0 during normal activity at high Ca^{2+}_{i} because of frequent sojourns to brief closed states that generated flickers, and also because of much less frequent sojourns to longer closed intervals with durations of typically 1–10 ms. Examples of the frequent flickers and of the less frequent longer closed intervals during normal activity at high Ca^{2+}_{i} (>100 μM) are presented in Fig. 1 B for data obtained at 132 and 1,024 μM Ca^{2+}_{i}. Analysis of the current records obtained at 132 and 1,024 μM Ca^{2+}_{i} indicated that only 2.1 ± 1.0% of the closed intervals had durations >1 ms, and even fewer, 0.026 ± 0.023%, had durations >10 ms. (The estimates obtained at both 132 and 1,024 μM Ca^{2+}_{i} from each of five channels were combined, as the percentages of longer closed intervals at the two different high Ca^{2+}_{i} were not significantly different.)
The increase in P_{o} with increasing Ca^{2+}_{i} was associated with increases in the mean open interval durations and decreases in the mean closed interval durations (Fig. 2B and Fig. C), consistent with previous reports for BK channels from skeletal muscle (Barrett et al. 1982; Magleby and Pallotta 1983; McManus and Magleby 1991) and mSlo BK channels (Nimigean and Magleby 1999).
The increase in mean open interval duration (Fig. 2 B, thick line; Hill coefficient of 1.02) was much less Ca^{2+} sensitive than the decrease in mean closed interval duration (Fig. 2 C, thick line; Hill coefficient of 3.48), suggesting that the major effect of Ca^{2+}_{i} is to decrease the durations of the closed intervals. The decreased filtering for channel B04 (Fig. 2 B, ⋄) would contribute to the briefer observed open times for this channel, as more of the flickers would be captured.
Just as P_{o} saturated at high Ca^{2+}_{i}, the mean durations of the open and closed intervals also saturated in high Ca^{2+}_{i} at ∼3.6 and ∼0.15 ms, respectively (Fig. 2B and Fig. C, thick line). Thus, any viable mechanism for the gating of the BK channel must account for a saturation in mean open times, mean closed times, and in P_{o} at high Ca^{2+}_{i}.
The Gating Kinetics Remain Unchanged at High Ca^{2+}_{i}
To further characterize the effect of high Ca^{2+}_{i} on gating, the open and closed interval durations during normal activity were measured and plotted as 1D dwelltime distributions in Fig. 3. The thick lines in the open and closed distributions in Fig. 3 are the fits with mixtures (sums) of three open and five closed exponential components, respectively. Increasing Ca^{2+}_{i} from 5.5 to 132 μM shifted the open intervals to longer durations (note rightward shift of the major peak for the open times) and the closed intervals to briefer durations (note leftward shift and decreasing amplitude of the peak describing the longer closed times.). A further increase in Ca^{2+}_{i} from 132 to 1,024 μM had little additional affect on either the open or closed distributions, as seen in Fig. 3D and Fig. H, where the fits to the distributions at 132 μM Ca^{2+}_{i} are plotted as thin lines on the fits to the distributions at 1,024 μM Ca^{2+}_{i}. The thin lines essentially superimpose the thick lines, indicating essentially unchanged gating kinetics at high Ca^{2+}_{i}. Fig. 3 H, inset, where the ordinate is plotted on a log scale to present the tails of the distributions at high gain, shows that the thin line also superimposes the thick line (within the range of the data) at the longer intervals where the frequency of occurrence of intervals is low. Results consistent with those in Fig. 3 were observed for four additional channels.
Table 1 presents the time constants and areas of the three open and five closed exponential components fitted to the 1D distributions in Fig. 3. For increases in Ca^{2+}_{i} up to 132 μM, the time constant of the briefest open component remained relatively unchanged, the time constants of the two longer open components increased with increasing Ca^{2+}_{i}, and the area from the briefest open component shifted into the longer open components. In contrast, the time constants of the longer closed components became briefer and their areas shifted from the longer to the briefer closed components. Increasing the Ca^{2+}_{i} from 132 to 1,024 μM then had little additional effect on the open and closed components, when compared with the large effects observed for lower concentrations of Ca^{2+}_{i}.
Table 2 presents the mean ± SD of the time constants and areas of the exponential components for data obtained from five channels at two levels of high Ca^{2+}_{i}. For four of the channels, the data were obtained at 132 and 1,024 μM Ca^{2+}_{i}, and for the fifth channel the data were obtained at 100 and 1,000 μM Ca^{2+}_{i}. The data from the fifth channel was pooled with the data from the other four since there was no apparent difference in the findings. Increasing the Ca^{2+}_{i} from 132 to 1,024 μM had no significant effect on either the time constants or areas of any of the open or closed components (Table 2, P > 0.05; paired t test). The measured increase in the time constant of the slowest closed component in Table 1 for channel B06 would be consistent with stochastic variation, as there were <20 intervals contributing to this component at each of the two levels of high Ca^{2+}_{i}.
The apparent lack of effect of an eightfold increase in Ca^{2+}_{i} on the dwelltime distributions at high Ca^{2+}_{i} suggest that the rates for the transitions between states that dominate the gating at high Ca^{2+} are either Ca^{2+} independent or saturate at high Ca^{2+}. These possibilities will be considered later.
Lack of Discrete Ca Block with High Ca^{2+}_{i}
It is well established that Ba^{2+} produces discrete block of BK channels (Vergara and Latorre 1983; Miller 1987). If Ca^{2+} produces a similar block, then a potential difficulty with conducting experiments in high Ca^{2+}_{i} is the possibility of discrete (slow) channel block by Ca^{2+}. Closed intervals arising from discrete Ca^{2+} block could then be mistaken for closed intervals arising from sojourns to closed states during the gating. The observations in Fig. 1 and Fig. 3 at high Ca^{2+}_{i} are consistent with our previous findings of a lack of evidence for discrete Ca^{2+} block of native BK channels from rat skeletal muscle (Rothberg et al. 1996). Increasing Ca^{2+}_{i} eightfold from 132 to 1,024 μM had little effect on the closed dwelltime distributions (Fig. 3C, Fig. D, Fig. G, and Fig. H). If appreciable closed intervals arose from discrete Ca^{2+} block, then increasing Ca^{2+}_{i} eightfold might have been expected to have a noticeable effect on the distributions, which was not observed.
As a more critical test for discrete block, we tabulated the frequency of occurrence of closed intervals >1 and >10 ms at the two different Ca^{2+}_{i}. For the five examined BK channels, the frequency of closed intervals with durations >1 ms was 7.4 ± 3.6 s^{−1} for 132 μM Ca^{2+}_{i} and 5.6 ± 2.5 s^{−1} for 1,024 μM Ca^{2+}_{i}, values that were not significantly different . The frequency of closed intervals with durations >10 ms was 0.12 ± 0.11 s^{−1} with 132 μM Ca^{2+}_{i} and 0.076 ± 0.086 s^{−1} with 1,024 μM Ca^{2+}_{i}, values that were also not significantly different .
Thus, the lack of effect of Ca^{2+} on the closed dwelltime distributions from brief to long times, and also on all closed intervals >1 and >10 ms, suggests that discrete Ca^{2+} block did not contribute to the closed intervals. As expected, increasing Ca^{2+}_{i} to 1,024 μM did decrease the conductance of the channel by ∼10%, presumably due to a screening (fast blocking) effect (Ferguson 1991). Cox et al. 1997a studying mSlo BK channels also found that high Ca^{2+}_{i} (1,000 μM) reduced singlechannel conductance through a fast block, but did not induce discrete block. Since high Ca^{2+}_{i} does not introduce closed intervals arising from discrete block, it is possible to estimate the numbers of kinetic states that contribute to the gating at high Ca^{2+}_{i}.
At Least Three to Four Open and Four to Five Closed States Are Entered during Gating in High Ca^{2+}_{i}
Assuming that the gating is consistent with a discretestate Markov model (McManus and Magleby 1989; Petracchi et al. 1991), the number of significant exponential components required to describe the dwelltime distributions gives an estimate of the minimum number of states entered during the gating (Colquhoun and Hawkes 1981, Colquhoun and Hawkes 1995a). Estimates of the numbers of significant exponential components were made from fitting both 1D and 2D dwelltime distributions (see methods). (Examples of 1D distributions were presented in Fig. 3 and examples of 2D distributions will be presented in a later section.) These estimates were then plotted against Ca^{2+}_{i} in Fig. 4.
In general, the 1D and 2D methods gave similar estimates of the numbers of exponential components, but there were some differences. 2D fitting can have an increased ability to detect components over 1D fitting when the same numbers of intervals are analyzed (Rothberg et al. 1997). However, with 2D fitting, when very brief intervals pair with longer intervals, both intervals in the pair are excluded from the fitting to avoid potential artifacts arising from fitting intervals with durations <1.5–2 dead times (see methods). The combination of the greater ability to detect components together with the fitting of fewer intervals for 2D fitting, when compared with 1D fitting, would contribute to the minor differences in the numbers of significant components estimated by the 1 and 2D fitting methods.
At the lower concentrations of Ca^{2+}_{i}, the distributions were typically described by three to four open and four to seven closed components for fitting with either 1 or 2D distributions (Fig. 4), consistent with previous observations (McManus and Magleby 1988; Wu et al. 1995; Rothberg and Magleby 1998). At the high concentrations of 132 and 1,024 μM Ca^{2+}_{i}, the open distributions were still typically described by three to four open components, and the closed distributions were typically described by four to five closed components. Thus, at kinetically saturating levels of Ca^{2+}_{i}, the gating typically involved transitions among at least three to four open and four to five closed kinetic states.
Theoretical Models for the Gating of Tetrameric Ligandgated Channels
Functional BK channels can be formed by four alpha subunits (Shen et al. 1994), and the BK channels in skeletal muscle, as studied here, are unlikely to be associated with the auxiliary beta subunit that increases the Ca^{2+} sensitivity of the channel (TsengCrank et al. 1996; Chang et al. 1997; Nimigean and Magleby 1999). If one or more Ca^{2+} ions bind to each alpha subunit to fully activate the channel (Schreiber and Salkoff 1997), then this would give a basis for the high Hill coefficients of 2–4 that are typically observed for activation (see introduction).
Although relatively simple in concept, with four subunits and at least one Ca^{2+}binding site per subunit, theoretical models for the gating of a ligandactivated homotetramer, such as the BK channel, can be highly complex, with 35–55 potential states (Eigen, 1969; Cox et al. 1997b). The 55state model (Cox et al. 1997b) is given by Scheme I, where each subunit can assume either of two conformations, indicated by squares or circles, and each subunit in either conformation can either be free of Ca^{2+} (open symbols) or bound with a Ca^{2+} (shaded symbols). Many of the states are isoforms, in which subunits with diagonal and adjacent conformational changes and Ca^{2+} bindings have potentially different functional properties. (Scheme I)
For the gating of the channel, the states in the top row of Scheme I are assumed to represent closed states of the channel, and the states in the bottom row are assumed to represent open states. The conductance of the states in the middle three rows is less clear, but may be open, closed, or partially conducting (Cox et al. 1997a). Scheme I reduces to the 35state model of Eigen 1968 if it is assumed that diagonal and adjacent subunits in the same conformation have identical properties. Because of the complexity of the 55 and 35state schemes, Scheme I is often reduced further to the 25state model described by Scheme II (Fersht 1985).
If it is further assumed that each collection of isoforms has the same properties, and that conformational changes are concerted so that all four subunits undergo conformational changes simultaneously or that the lifetimes of the middle three rows of states in Scheme I and Scheme II are very brief, then Scheme I and Scheme II reduce to the 10state MonodWymanChangeux model for allosteric proteins (Monod et al. 1965) described by Scheme III. (Scheme II) (Scheme III)
Simple Models Fail to Predict the SingleChannel Kinetics at High Ca^{2+}_{i}
It has been shown previously for low to intermediate levels of Ca^{2+}_{i} that the Ca^{2+} dependence of the singlechannel kinetics of BK channels in rat skeletal muscle can be approximated by the gating mechanisms described by Fig. 4, Scheme V, Scheme VI (McManus and Magleby 1991; Rothberg and Magleby 1998). Thus, we first examined whether these schemes might also account for the singlechannel kinetics through high Ca^{2+}_{i}.
Scheme IV is drawn from the MWC model. Scheme VI is an expansion of Scheme IV because it includes closed states beyond the activation pathway. Scheme V can be viewed as a condensed version of Scheme I and Scheme II if it is assumed that the intermediate states in Scheme I and Scheme II are too brief to be detected or that only a subset can be detected. It is the additional brief closed states C9, C10, and C11 that generate most of the flickers (brief closings) in Scheme V and Scheme VI, whereas flickers for Scheme IV are generated mainly by sojourns between states O2 and C5. Flickers are highly characteristic of singlechannel currents and can be seen in Fig. 1. (Scheme IV) (Scheme V) (Scheme VI)
Before examining whether Scheme IV, Scheme V, Scheme VI could account for the gating from low through high Ca^{2+}_{i}, we first examined whether they could describe the 1D dwelltime distributions from low to intermediate levels of Ca^{2+}_{i}, as reported previously (McManus and Magleby 1991; Rothberg and Magleby 1998). Rate constants for each scheme were estimated by the simultaneous fitting of 2D dwelltime distributions obtained at three different Ca^{2+}_{i} of 5.5, 8.3, and 12.3 μM (see methods). The most likely rate constants were then used to simulate singlechannel data with noise and filtering equivalent to that of the experimental data. The simulated current records were then analyzed in the same manner as the experimental data to determine the predicted 1D dwelltime distributions. As expected, the predicted distributions gave excellent descriptions of the observed 1D dwelltime distributions from low to intermediate Ca^{2+}_{i} for the three channels examined in detail. For example, for the data in Fig. 3A, Fig. B, Fig. E, and Fig. F, at 5.5 and 12.3 μM Ca^{2+}_{i}, the predicted distributions typically overlapped or were within a line width of the thick lines describing the dwelltime distributions in these figures (not shown).
We next examined whether Scheme IV, Scheme V, Scheme VI could predict the distributions from low to high Ca^{2+}_{i}. Rate constants for each scheme were estimated by the simultaneous fitting of 2D dwelltime distributions obtained at six different Ca^{2+}_{i} of 5.5, 8.3, 12.3, 20.3, 132, and 1,024 μM. These rate constants were then used to predict the observed distributions over the full range of Ca^{2+}_{i}. As shown in Fig. 5, Fig. 4 could not simultaneously describe the dwelltime distributions from low to high Ca^{2+}_{i}. Scheme V gave a better description than Scheme IV, but still could not describe the distributions. The predictions of Scheme VI were very slightly better than those of Scheme V, and are not shown. Table 3 shows the schemes rankings: Scheme VI > Scheme V > Scheme IV (Akaike rankings, Formula 4).
The inability of these schemes to account for the singlechannel kinetics from low to high Ca^{2+}_{i} indicates that the models described by Scheme IV, Fig. 5, Scheme VI are too simple. Such a finding is, perhaps, not surprising since Fig. 4 with eight states and Scheme V and Scheme VI with 11 states include only a small subset of the minimal 55 potential states for the gating of a ligandactivated homotetrameric channel based on theoretical considerations (Scheme I). Nevertheless, it is possible that the channel does not gate as described by the theoretical 55state model or, if it does, that only a fraction of the potential 55 states actually contribute to the gating. Consequently, we examined the differences between the observed and predicted responses for these schemes to determine how the minimal schemes might be expanded to be more consistent with the gating of the channel.
For high Ca^{2+}_{i}, Scheme IV, Scheme V, Scheme VI predicted fewer brief openings than were observed in the singlechannel data (Fig. 5). Brief openings in Scheme V typically arose from sojourns such as C11O3C11, and in Scheme VI from sojourns such as C11O3C11 or C6O3C6, as O3 was the open state with the briefest lifetime in both schemes. For these schemes, the channel would be unlikely to reach O3 or the associated closed states at high Ca^{2+}_{i}, as the high Ca^{2+}_{i} would drive the gating towards the fully liganded states. Thus, to generate more brief open intervals at high Ca^{2+}_{i}, there needs to be a means for the gating to reach directly one or more brief open states from the fully liganded closed states. In addition, the underprediction of the long closed intervals at high Ca^{2+}_{i} suggests that the channel may also gate among additional closed states at high Ca^{2+}_{i} not included in Scheme IV, Scheme V, Scheme VI. Further evidence that the channel gates among additional states at high Ca^{2+}_{i} arises from the observation in Fig. 4 that the channels typically entered at least three to four open and four to five closed states at high Ca^{2+}_{i}, while analysis of the dwelltime distributions predicted by Scheme IV, Scheme V, Scheme VI at high Ca^{2+}_{i} indicated that the distributions were described by only one open and two to three closed components.
Scheme VII Can Describe the Open and Closed 1D DwellTime Distributions from Low to High Ca^{2+}_{i}
We next examined how Scheme V might be expanded to provide the required access to more open and closed states at high Ca^{2+}_{i}. Scheme V is contained within the general Scheme I and Scheme II if it is assumed that the three rows of intermediate states in the general schemes are closed states that can be collapsed into one row of brief lifetime intermediate closed states. However, since it is possible that the last two rows of states in Scheme I and Scheme II are open states and that the remaining intermediate closed states do not collapse into one row of closed states, we investigated whether simplified schemes drawn from Scheme II, with an assumption of three rows of closed states and two rows of open states, could account for the data.
Scheme VII presents a gating mechanism of this type, where open states O4O5O6 (with brief lifetimes) and closed states C12C13C14 are the additional rows of open and closed states when compared with Scheme V. With high Ca^{2+}_{i}, the channel could now make sojourns from closed states to brief open states, such as C15O4C15, giving brief open intervals, and also make sojourns from closed states to longer open (and compound open) states, such as C15O4O1O4C15, giving longer open intervals. The extra row of closed states should also allow the generation of longer closed intervals at high Ca^{2+}_{i} by allowing more sojourns among closed states between openings. (Scheme VII)
The most likely rate constants for Scheme VII were determined from the simultaneous fitting of 2D dwelltime distributions at six different Ca^{2+}_{i} (5.5, 8.3, 12.3, 20.3, 132, and 1,024 μM). The thick lines in Fig. 5 show that the additional states allowed Scheme VII to describe the 1D open and closed dwelltime distributions from low to high Ca^{2+}_{i}. In this scheme, the mean lifetimes at 5.5 μM Ca^{2+}_{i} of states O4, O5, and O6, of 0.11, 0.09, and 0.03 ms, tend to be brief compared with the lifetimes of the final row of open states O1, O2, and O3 of 0.50, 0.18, and 0.06 ms. The improved ability of Scheme VII to describe the singlechannel gating when compared with Scheme IV, Scheme V, Scheme VI (compare thick line to dotted and thin lines in Fig. 5) is also reflected in the greatly improved likelihood ratios in Table 3.
We also examined whether a scheme like Scheme VII, but with one less row of intermediate closed states, could account for the data. For the three channels examined, the likelihood estimates were two to four orders of magnitude less than for Scheme VII, and the reduced scheme ranked below Scheme VII for all three channels (not shown).
Using Kinetic Structure to Evaluate the Gating Mechanism
Scheme VII accounts for the 1D open and closed dwelltime distributions from low to very high Ca^{2+}_{i} (Fig. 5). Such a description would be sufficient to predict P_{o} as a function of Ca^{2+}_{i} over a wide range of activity, but the 1D distributions do not take into account the correlation information between adjacent intervals, which can give insight into the connections (transition pathways) among the various states (Magleby and Song 1992; Rothberg and Magleby 1998). To examine the correlation information, we determined whether Scheme VII could account for the kinetic structure of the singlechannel data. The kinetic structure is described by 2D dwelltime distributions and dependency plots. The 2D distributions indicate the frequency of occurrence of pairs of adjacent open and closed intervals (Fredkin et al. 1985; Magleby and Weiss 1990b), and the dependency plots convey information about the correlations of adjacent interval durations (Magleby and Song 1992).
The kinetic structure for the same channel featured in the previous figures (channel B06) is shown in Fig. 6 at four different Ca^{2+}_{i}. The 2D dwelltime distributions (Fig. 6, left) are plotted on log–log coordinates with the logs of the durations of adjacent open and closed intervals locating the position of the bin on the x and y axis, respectively. The z axis plots the square root of the numbers of intervals in each bin. These 2D dwelltime distributions thus extend the Sigworth and Sine 1987 transform used in the previous figures to two dimensions. From the 2D dwelltime distributions it can be seen that pairs of long open intervals adjacent to brief closed intervals (flickers) occur most frequently of all the interval pairs, and this is the case from low to high Ca^{2+}_{i} (Fig. 6, left, position 4). It is these interval pairs that give rise to the characteristic longer openings separated by flickers in the experimental data (Fig. 1). At the lowest Ca^{2+}_{i} of 5.5 μM, there were also large numbers of longer open intervals adjacent to longer closed intervals (Fig. 6, left, position 6) and brief open intervals adjacent to longer closed intervals (left, position 3). As the Ca^{2+}_{i} was raised, the longer closed intervals shifted to briefer durations. The peak at position 6 with 5.5 μM Ca^{2+}_{i} shifted towards position 5 at 12.3 μM Ca^{2+}_{i}, and position 4 at higher Ca^{2+}_{i}.
The dependency plots in Fig. 6 (right) present the fractional excess or deficit of interval pairs of specified durations over that expected if the intervals paired at random. Dependencies of +0.5 or −0.5 would indicate a 50% excess or 50% deficit of interval pairs over the number expected if open and closed intervals paired independently. The thick lines indicate a dependency of zero. Because the dependency plots present magnified representations of excesses and deficits in the numbers of observed interval pairs relative to the numbers expected for independent pairing, they must be interpreted with some caution, as the estimates of dependency can be unreliable where the numbers of observed interval pairs per bin in the 2D dwelltime distributions are small. Consequently, references to dependency will only be made when the referenced dependencies are known to be significantly different from zero. Such dependencies will be referred to by numbers on the dependency plots. Examples of which areas of the dependency plots are significantly different from zero are presented in Rothberg and Magleby 1998 for low to intermediate levels of Ca^{2+}_{i} and will be presented in a later section for data obtained at high Ca^{2+}_{i}.
The dependency plots in Fig. 6 (right) indicate that over the wide range of examined Ca^{2+}_{i} there was a deficit of brief open intervals adjacent to brief closed intervals (position 1), an excess of brief open intervals adjacent to both intermediate and long closed intervals (positions 2 and 3), and a deficit of long open intervals adjacent to long closed intervals (position 6; not clearly visible in the presented orientation, but visible when the plots were rotated). At the lower Ca^{2+}_{i} of 5.5 and 12.3 μM, there was also an excess of longer open intervals adjacent to brief closed intervals (position 4). These specific excesses and deficits of interval pairs give rise to the characteristic saddle shape of the dependency plots for BK channels (Rothberg and Magleby 1998), and were consistently seen for all the examined channels.
While the kinetic structure at intermediate levels of Ca^{2+}_{i} for six different BK channels has been presented previously (Rothberg and Magleby 1998), it was important to determine whether the kinetic structure at high Ca^{2+}, as shown in Fig. 6, was consistently observed. (It will be shown in a later section that the presence or absence of significant dependencies at high Ca^{2+}_{i} is a key factor in distinguishing gating mechanism.) Fig. 7 presents such plots for data obtained at 132 and 1,024 μM Ca^{2+}_{i} for two additional channels. Although there were some obvious differences in the magnitudes of the dependencies, the general shape of the kinetic structure in these plots was the same as for the channel shown in Fig. 6G and Fig. H, in that there were obvious deficits of brief open intervals adjacent to brief closed intervals (Fig. 7, right, position 1) and obvious excesses of brief open intervals adjacent to intermediate duration closed intervals (right, position 2).
Dependency significance plots were made for the data obtained at high Ca^{2+}_{i} to estimate which dependencies were significantly different from zero. A paired t test was used to compare the number of interval pairs in each bin of the observed 2D dwelltime distribution with the number expected if adjacent open and closed intervals paired independently, by using a moving 3 × 3 bin array as detailed in Rothberg and Magleby 1998. Results are shown in Fig. 8 for channel B12. The plots in Fig. 8A and C, present the dependency significance for the front and back views, respectively, of the dependency plot in Fig. 7 E obtained at 132 μM Ca^{2+}_{i}, and the plots in Fig. 8B and Fig. D, present the dependency significance for the front and back views of the dependency plot in Fig. 7 F obtained at 1,024 μM Ca^{2+}_{i}. The dependency significance plots present the logarithm of the estimated P value, which is multiplied by the sign of the dependency to indicate whether the paired intervals are in excess or deficit. The thick lines on the plots at −1.3 and 1.3 indicate a significance level of . Absolute values of dependency significance > 1.3, 2, 3, and 4 would indicate P < 0.05, 0.01, 0.001, and 0.0001, respectively.
From Fig. 8 it can be seen that the dependencies at the numbered positions 1, 2, 4, and 5 were significantly different from zero at both 132 and 1,024 μM Ca^{2+}_{i}. Importantly, there was a significant excess of long open intervals adjacent to brief closed intervals (position 4), even though the fractional excess of these interval pairs was small in the dependency plots (Fig. 7E and Fig. F). The reason for this apparent discrepancy is that most of the interval pairs fall at position 4 in high Ca^{2+}_{i}, as can be seen from the 2D dwelltime distributions (Fig. 7, A–D), so that even an appreciable excess of interval pairs at position 4 would still appear small when plotted as dependency, which plots the fractional excess of intervals (Formula 1). In six of six dependency significance plots that were examined at high Ca^{2+}_{i}, the dependencies at positions 1, 2, and 5 were significant, and the dependency at position 4 was significant in five of six plots. The one plot where significance was not observed at position 4 had fewer numbers of analyzed intervals.
The observations that the numbers of detected kinetic states remained relatively unchanged from low to high Ca^{2+}_{i} (Fig. 4) and that the general shapes of the dependency plots also remained relatively unchanged from low to high Ca^{2+}_{i} (Fig. 6 and Fig. 7) raise the possibility that the basic gating mechanism remains relatively unchanged from low to high Ca^{2+}_{i}.
Scheme VII Describes the Kinetic Structure at Low but not High Ca^{2+}_{i}
To determine whether Scheme VII could account for the kinetic structure, the most likely rate constants for Scheme VII, determined from the simultaneous fitting of 2D dwelltime distributions at six different Ca^{2+}_{i} (5.5, 8.3, 12.3, 20.3, 132, and 1,024 μM), were used to simulate singlechannel data for Scheme VII with noise and filtering equivalent to that of the experimental data. The simulated current records were then analyzed in the same manner as the experimental data to determine the predicted kinetic structure shown in Fig. 9.
Scheme VII captured the basic features of the kinetic structure at 5.5 and 12.3 μM Ca^{2+}_{i} (compare the predicted 2D dwelltime distributions and dependency plots in Fig. 9A, Fig. B, Fig. E and Fig. F, to the observed plots in Fig. 6, A, B, E, and F). Scheme VII also captured the basic features of the 2D dwelltime distributions at the high Ca^{2+}_{i} of 132 and 1,024 μM (compare Fig. 9C and Fig. D with Fig. 6C and Fig. D). However, Scheme VII predicted that little or no dependence would be observed at high Ca^{2+}_{i} (Fig. 9G and Fig. H), in contrast to the significant dependencies observed in the experimental data (Fig. 6G and Fig. H, and Fig. 7 and Fig. 8). Thus, Scheme VII predicted that open and closed intervals would pair randomly at high Ca^{2+}_{i}, in contrast to the dependent pairing observed in the data. These observations indicate that Scheme VII is too simple to capture the features of the gating at high Ca^{2+}_{i}.
The Gating Is Effectively Restricted to the Fully Liganded States of Scheme VII at High Ca^{2+}_{i}
To explore why Scheme VII did not predict the dependency at high Ca^{2+}_{i}, we calculated the equilibrium occupancies of the various open and closed states for this scheme. At low to intermediate Ca^{2+}_{i} (5.5–12.3 μM), the channel readily entered all the states in Scheme VII, with occupancy biased towards the closed states C8–C11 at the low Ca^{2+}_{i}. At high Ca^{2+}_{i} (132 and 1,025 μM), the channel spent 99.6% of its time in the fully liganded states, with 96.6% in open states 1 and 4 and 2.8% in closed states 7, 12, and 15. (The calculated equilibrium occupancies for Fig. 7 can be found in online supplemental Figure S1 [http://www.jgp. org/cgi/full/114/1/93/DC1].)
The reason that Scheme VII predicted a lack of dependence between adjacent open and closed intervals at high Ca^{2+}_{i} is now readily apparent. Because the channel spent 99.4% of its time in the fully liganded column of states at high Ca^{2+}_{i}, the channel would essentially gate as Scheme VIII. Scheme VIII has a single effective transition pathway between the open and closed states, given by C15O4. A single effective transition pathway gives a single gateway state, which would lead to independent pairing of open and closed intervals and lack of significant dependencies (McManus et al. 1985; Colquhoun and Hawkes 1987; Song and Magleby 1994). (Scheme VIII)
If Scheme VII does effectively gate in the fully liganded column of states at high Ca^{2+}_{i}, as indicated by Scheme VIII, then a maximum of two open and three closed components would be detected in the dwelltime distributions predicted by Scheme VII at high Ca^{2+}_{i}. This was found to be the case. Fitting exponentials to dwelltime distributions simulated with Scheme VII at high Ca^{2+}_{i} gave two open and three closed components, compared with the typically three to four open and four to five closed states detected in the experimental data. Scheme VII also predicted a lack of dependence and too few components at high Ca^{2+}_{i} for the two other channels analyzed in detail.
Scheme I, Fig. 2, and Scheme VII Can Be Rejected
The above findings indicate that Scheme VII can thus be rejected as a model for gating, as it cannot describe the dependencies at high Ca^{2+}_{i}. By analogy, Fig. 1 and Scheme II and all schemes based on subsets of states drawn from these schemes, such as Scheme II, Scheme III, Scheme IV, Scheme V, Scheme VI, can also be rejected as such schemes would also not describe the dependencies at high Ca^{2+}_{i}, provided that the forward rate constants for binding of Ca^{2+} are sufficiently rapid at high Ca^{2+} to effectively keep the gating in the fully liganded column of states.
The Kinetic Structure at High Ca^{2+}_{i} Suggests a TwoTiered Gating Mechanism
The above results suggest that models for gating at high Ca^{2+}_{i} must: (a) allow effective transitions among at least three to four open and four to five closed states at high Ca^{2+}_{i} to generate the required numbers of exponential components in the dwelltime distributions (Fig. 4), (b) have two or more independent transition pathways between open and closed states (two or more gateway states) at high Ca^{2+}_{i} to generate a dependent relationship between the durations of adjacent intervals (Fig. 6, Fig. 7, Fig. 8), and (c) gate among the three to four open and four to five closed states at high Ca^{2+}_{i} in a manner that is essentially independent of Ca^{2+}_{i} for Ca^{2+}_{i} > 100 μM to account for the observed lack of effect of Ca^{2+} on the gating at high Ca^{2+}_{i} (Fig. 1, Fig. 2, Fig. 3 and Table 1 and Table 2).
Since Scheme I and Scheme II are theoretical schemes based on a ligandactivated homotetramer, it might be useful if these schemes could be modified to account for gating at high Ca^{2+}_{i}. As pointed out by Cox et al. 1997b, it is not necessarily clear which of the states in Scheme I are open and which are closed. Perhaps all 55 (Scheme I) or 25 (Scheme II) states are closed, each with the potential of opening through a concerted conformational change of the subunits. Concerted conformational changes leading to opening have been suggested previously for the gating of Shaker channels (Bezanilla et al. 1994; Zagotta et al. 1994; Zheng and Sigworth 1998). If each closed state can open for BK channels, then Scheme I and Scheme II would become twotiered 110 and 50state models, respectively, with the upper tiers composed entirely of closed states and the lower tiers of an equal number of open states. Preliminary data using macroscopic ionic and gating currents from mSlo are consistent with such models (Horrigan and Aldrich 1998). The extension of Scheme II to a twotiered model is given by Scheme IX, where in the graphic representation all 25 closed states in the upper tier are visible and only the first two rows of the 25 open states in the lower tier are visible
Scheme X Can Approximate the Kinetic Structure From Low to High Ca^{2+}_{i}
A difficulty with Scheme IX, as was also the case for Scheme I and Scheme II, is that Scheme IX has so many rate constants that it would be difficult if not impossible to determine unique rate constants for this scheme, even by the simultaneous fitting of 2D dwelltime distributions over a range of Ca^{2+}_{i} (see methods). Consequently, we examined whether a reduced version of Scheme IX might be sufficient to approximate the kinetic structure. The number of rows of closed states in the upper tier was reduced from five to three, and the number of rows of open states in the lower tier was reduced from five to two to obtain the reduced Scheme X. Some of the open and closed states with zero and one bound Ca^{2+} were also omitted in Scheme X to reduce further the number of rate constants because, over the range of Ca^{2+}_{i} examined in this paper, it might be expected that these states would contribute little to the gating (Cui et al. 1997; Rothberg and Magleby 1998). (Scheme IX) (Scheme X)
Scheme X is like Scheme VII, except that Scheme X has a total of six independent transition pathways between the open and closed states, compared with three in Fig. 7. Two of the independent transition pathways connect fully liganded open and closed states, which would then allow at least two functional gateway states for gating with high Ca^{2+}_{i}. Two or more gateway states would be needed to generate the dependencies observed in high Ca^{2+}_{i}.
Scheme X was tested by determining the most likely rate constants from the simultaneous fitting of 2D dwelltime distributions obtained at six different Ca^{2+}_{i} ranging from 5.5 to 1,024 μM. The most likely rate constants were then used with Scheme X to obtain the predicted kinetic structure in Fig. 10. The predicted kinetic structure (compare to the observed in Fig. 6) shows that Scheme X captured the major features from low to high Ca^{2+}_{i}, including the general shapes of the dependency plots at high Ca^{2+}_{i}. Scheme X predicted the excess of brief open intervals adjacent to the longer closed intervals (position 2) and the deficits of brief open intervals adjacent to brief closed intervals (position 1) that were not predicted by Scheme VII at high Ca^{2+}_{i} (compare Fig. 10G and Fig. H, with Fig. 9 and Fig. 6G and Fig. H). For intervals with durations less than ∼0.05 ms, Scheme X predicted too great of a deficit at position 1. This could reflect an inadequacy of the model or it could reflect the fact that dwell times <0.05 ms were not fitted, so that the predictions of the model were not constrained below this time.
As would be expected from the reasonable descriptions of the kinetic structure, Scheme X also described the Ca^{2+} dependence of the 1D dwelltime distributions. The distributions predicted by Scheme X essentially superimposed the thick lines in Fig. 5. Scheme X also described the Ca^{2+} dependence of P_{o} and of the mean open and closed interval durations (Fig. 2, thin, dotted, and dashed lines).
Assessing the Various Gating Mechanisms
Scheme X could also describe the kinetic structure obtained from the two other channels examined over a wide range of Ca^{2+}_{i}. The rankings of the various kinetic schemes for the three channels together with the NLR_{1000} are presented in Table 3. The NLR_{1000}, which gives a measure of how well the schemes describe the data (see below and methods), indicated that Scheme X was more likely than the other examined schemes for all three channels (larger values of NLR_{1000} are more likely).
While the NLR_{1000} can indicate which schemes are most likely, it does not apply any penalties for additional free parameters. Consequently, the schemes were ranked by the Akaike criteria, which applies a penalty for additional free parameters (Formula 4). The general rankings were: Scheme X > Scheme VII > Scheme IV, Scheme V, Scheme VI (Table ). The Akaike test ranks schemes, but does not give the significance of the rankings. The likelihood ratio test can be used to estimate the significance of rankings for nested models (Horn and Lange 1983; see examples in McManus and Magleby 1988). Scheme VII and Scheme X ranked significantly above Fig. 4, Scheme V, Scheme VI and Scheme Vsat (P < 0.001) for all three channels. Interestingly, Scheme X ranked above Scheme VII for only two of three channels, and this ranking was significant for only one channel (B06, P < 0.001) in spite of the fact that Scheme X gave better likelihoods than Scheme VII for all three channels.
The apparent discrepancy between visual observations and the significance of some of the rankings may reflect the necessarily conservative nature of statistical tests. Alternatively, the discrepancy may reflect that Scheme X is still too simple, so that obvious improvements in some aspects of the gating, such as in the dependency plots, are countered by minimal improvements or even small detrimental changes in other aspects of the gating. Thus, Scheme X may have to be expanded into Scheme IX to obtain sufficiently improved descriptions of the data to outweigh the heavy penalty imposed by the ranking tests. Consistent with this possibility, Scheme X predicts only two open and three closed exponential components at high Ca^{2+}_{i}. It will be shown in a later section that the gating of the fully liganded channel is described better by models that more closely approximate Scheme IX.
The NLR_{1000} values in Table 3 give a numerical measure of how well the various schemes describe the experimental data. A NLR of 1.0 indicates that a kinetic scheme describes the 2D dwelltime distributions as well as the theoretical best description for a discrete state Markov model (see methods). The thick lines in Fig. 3 show the theoretical best description of the 1D distributions. For channel B06, the values of NLR_{1000} (normalized to 1,000 interval pairs) ranged from 3.15 × 10^{−33} for Scheme IV to 3.25 × 10^{−4} for Scheme X. These values give likelihood ratios per interval pair of 0.93 [(3.15 × 10^{−33})^{0.001}] for Scheme IV and 0.99 [(3.25 × 10^{−4})^{0.001}] for Scheme X. Such values suggest an average likelihood difference per interval pair between the predicted and theoretical best fits of 7% for Scheme IV and only 1% for Scheme X. The 7% difference per interval pair is readily seen (Fig. 5, dotted line), while the small 1% difference is still visually apparent as less than perfect descriptions of the data (Fig. 5, thick line, and compare Fig. 10 to Fig. 6).
Scheme X Predicts too Few Components at High Ca^{2+}_{i}
In spite of its relative success, Scheme X is still too simple. Analysis of simulated data indicated that Scheme X predicted only two significant open and four significant closed components at high Ca^{2+}_{i}, compared with the three to four open and five closed components in the experimental data. This underprediction is not surprising, since for Scheme X the high Ca^{2+}_{i} would effectively drive the gating towards the two open and three closed fully liganded states.
Models for the Gating of the Fully Liganded Channel
Since Scheme X predicts too few components at high Ca^{2+}_{i}, we explored what types of gating mechanisms might be consistent with the gating at high Ca^{2+}_{i}. In theory, Scheme IX could be examined directly, but the data would be insufficient to constrain the large numbers of rate constants for the 50state model. Consequently, we explored the reduced models given by Scheme XI, Scheme XII, Scheme XIII, which are all composed of fully liganded states.
Scheme XI is drawn from the 10state model describing the fully liganded states in Scheme IX (the rightmost column of states) with two open states excluded because they had no significant effect on the likelihood estimates. Scheme XII expands Scheme XI so that transitions to the open states pass through intermediate states, and Scheme XIII has closed states beyond the activation pathway. Scheme XII and Scheme XIII were examined because BK channels gate with large numbers of brief closings (flickers) at high Ca^{2+}_{i}, just as they do at lower Ca^{2+}_{i} (Fig. 1, Fig. 3, and Fig. 6). Intermediate and/or secondary states are associated with the generation of flickers at lower Ca^{2+}_{i} (Rothberg and Magleby 1998). Additional support for possible intermediate closed states between the closed and open states comes from the observation of Cui et al. 1997 that there is a 50–150μs delay in the activation of BK channels by voltage steps. Indirect support for possible secondary closed states comes from observations on Shaker K^{+} channels, where secondary states appear to contribute to the gating (Hoshi et al. 1994; Zagotta et al. 1994; Schoppa and Sigworth 1998). Homology between Shaker channels and the core region of BK channels (Wei et al., 1994; Toro et al. 1998) then suggests the possibility of considering such secondary states for BK channels. (Scheme XI) (Scheme XII) (Scheme XIII) (Scheme XIV)
Scheme XI, Scheme XII, Scheme XIII were fitted to the 2D dwelltime distributions obtained at a single high Ca^{2+}_{i} of 1,024 μM for each of three channels. All three of these schemes gave reasonable descriptions of the kinetic structure at high Ca^{2+}_{i} (not shown, but slightly better visually when compared with the experimental data than that shown in Fig. 10D and Fig. H), and all three schemes gave detected numbers of open and closed states within the range observed in the experimental data at high Ca^{2+}_{i}. Scheme XI gave three open and four closed components, and Scheme XII and Scheme XIII each gave three open and five closed components. For purposes of comparison, a simpler Fig. 14 with only two open and three closed states, which describes the fully liganded states in Scheme X, was also examined.
Both the NLR_{1000} and Akaike criteria ranked the schemes in the order: Scheme XII ∼ Scheme XIII > Scheme XI >> Scheme XIV (Table 4). The likelihood ratio test was applied to the nested Scheme XIII, Scheme XI, and Fig. 14, and gave highly significant (P < 0.001) rankings of: Scheme XIII > Scheme XI > Scheme XIV.
The values of the NLR_{1000} for Scheme XI, Scheme XII, and Scheme XIII (Table 4) indicated that the fits given by these schemes were considerably better than for the simpler Scheme XIV, and approached the theoretical best descriptions of the single data sets for discrete state Markov models. The values of the NLR_{1000} for these schemes ranged from 0.061 to 0.959, giving likelihood ratios per interval pair ranging from 0.9972 (0.061^{0.001}) to 0.9999 (0.899^{0.001}), suggesting little difference in likelihood per interval pair between the observed and theoretical best descriptions of the data.
These findings indicate that the gating of BK channels at high Ca^{2+}_{i} can be approximated by models based on the fully liganded states in Scheme IX. Scheme XI was drawn from the fully liganded states in Scheme IX. Adding three additional brief states as either intermediate states (Scheme XII) or secondary states (Scheme XIII) to generate additional flickers significantly improved the description of the data (Table 4). Whether the intermediate or secondary states are needed, or whether these additional states simply provide a means to compensate for the fact that fitting only the fully liganded states excludes potential contributions to the gating from transitions back to the states with fewer than four bound Ca^{2+} is not yet clear. What is clear, however, is that at least three open and five to eight closed states, as described by Scheme XI, Scheme XII, Scheme XIII, are required to describe the gating at high Ca^{2+}_{i} equivalent to the theoretical best description. These schemes lack Ca^{2+}dependent rate constants and apply only for Ca^{2+}_{i} > ∼100 μM, where the gating kinetics are little affected by Ca^{2+}_{i}.
Scheme V with an Assumption of Saturation in the Ca^{2+}dependent Rate Constants Could Approximate the Gating from Low to High Ca^{2+}_{i}
We also explored an alternative explanation to account for the lack of effect of Ca^{2+}_{i} on the gating at high Ca^{2+}_{i}. In all of the above considered schemes, the binding rate of Ca^{2+}_{i} was assumed to be a first order reaction, increasing linearly with Ca^{2+}_{i}. Thus, the effective rate constants for binding are given by the product of Ca^{2+}_{i} times the rate constants expressed per micromole per second. An upper limit for the rate constant for such a diffusioncontrolled process is ∼10^{9} M^{−1} s^{−1} (Fersht 1985; Cui et al. 1997). There is, however, no a priori reason to think that the effective binding rate would necessarily increase linearly with Ca^{2+}_{i} at high Ca^{2+}_{i}. The physical structure of the Ca^{2+}binding sites is not yet known, but if the binding sites are in a vestibule with some additional negative charged groups, then the local concentration of Ca^{2+} at the binding sites at lower Ca^{2+} could be greater than that in the bulk solution (Van der Kloot and Cohen 1979; Green and Andersen 1991; Nonner and Eisenberg 1998), so that the local concentration could reach a maximum as the concentration of Ca^{2+} in the bulk solution is raised. This could give an apparent saturation in the Ca^{2+}dependent rate constants.
Alternatively, if the binding, which is represented by a onestep process in the kinetic schemes, is actually a twostep process that involves binding followed by a conformational change, then the apparent binding rate would saturate if the second step becomes rate limiting at high Ca^{2+}_{i} (Fersht 1985).
Since the physical details involving Ca^{2+} binding and action are not known, we explored these two saturation models by using an approach that was independent of a detailed physical model. As both processes would have the effect of reducing the effective concentration of Ca^{2+}_{i} at high Ca^{2+}_{i}, we examined whether the kinetic structure from low to high Ca^{2+}_{i} could be described by letting the effective Ca^{2+}_{i} at high Ca^{2+}_{i} be less than the actual Ca^{2+}_{i}, to mimic apparent saturation of the binding step. Since it is not known what the value of the effective Ca^{2+}_{i} would be at high Ca^{2+}_{i}, this value was estimated by iterative fitting. The 2D dwelltime distributions obtained at six different Ca^{2+}_{i} (5.5, 8.3, 12.3, 20.3, 132, and 1,024 μM) were simultaneously fitted to estimate the most likely rate constants for Scheme V, and also the most likely effective concentrations of Ca^{2+}_{i} for the data obtained at 132 and 1,024 μM Ca^{2+}_{i}. When fitting, the Ca^{2+}_{i} used for the data obtained at the four lower Ca^{2+}_{i} was fixed to the experimental values.
Scheme V with effective Ca^{2+}_{i} of 56.9 and 60.0 μM for the data obtained at 132 and 1,024 μM Ca^{2+}_{i}, respectively, and the actual Ca^{2+}_{i} for the other four data sets obtained at lower Ca^{2+}_{i} could approximate the basic features of the kinetic structure from low to high Ca^{2+}_{i}. The predicted kinetic structure at 1,024 μM Ca^{2+}_{i} is shown in Fig. 11 and was visually indistinguishable from the predicted kinetic structure at 132 μM Ca^{2+}_{i}. Comparison of the predicted kinetic structure in Fig. 11 to that in Fig. 6C, Fig. D, Fig. G, and Fig. H, showed that Scheme V with an assumption of saturation could approximate the data at high Ca^{2+}_{i}. Scheme V with saturation also described the data at lower Ca^{2+}_{i}, with the predicted structure similar to that in Fig. 10 (not shown). For the other two channels studied in a similar manner, the effective Ca^{2+}_{i}s at 132 and 1,024 μM Ca^{2+}_{i} were 67.8 and 114 μM (channel B12) and 35.2 and 39.4 μM (channel B14). The likelihoods indicated that Scheme V with an effective saturation in the binding rate (Scheme Vsat) described the kinetic structure from low to high Ca^{2+}_{i} slightly less well than Scheme VII and Scheme X (Table 3). The observation that the effective Ca^{2+}_{i} was greater for all three channels for the data collected at 1,024 μM than for the data collected at 132 μM suggests that increasing Ca^{2+}_{i} from 132 to 1,024 μM may have some additional effects on the gating, but any effects would be small since the exponential components describing the dwelltime distributions at 132 and 1,024 μM Ca^{2+}_{i} were not significantly different. Cox et al. 1997b have also found (for mSlo) that high Ca^{2+}_{i} may have additional effects on gating.
The results in this section show that a relatively simple gating mechanism (Scheme V) with the added assumption of apparent saturation in the Ca^{2+} binding steps at high Ca^{2+}_{i} can approximate the gating from low to high Ca^{2+}_{i}. It will be discussed later that saturating models may be less appropriate than twotiered models.
Estimated Rate Constants for the Examined Kinetic Schemes
Fig. 12 presents the estimated rate constants for Scheme X and Scheme XI for the three channels examined in detail. Online supplemental Figure S2 (http://www. jgp.org/cgi/full/114/1/93/DC1) presents estimated rate constants for these same three channels for most of the other examined schemes. The rate constants for the examined schemes typically ranged from ∼0 to 45,000/s, indicating a large range in the height of the energy barriers between the various states. Rate constants were limited so as not to exceed 45,000/s, as letting them go higher gave little improvement in the fits. Estimated rate constants for the simpler models (Schemes IV–VI and XI–XIII) were relatively consistent from channel to channel. For the most complex gating mechanism examined (Scheme X), there could be considerable variability in the estimates, depending on the specific rate constants. In those cases where there was considerable variability in estimates of the rate constants among channels, the variability was typically associated with poorly defined rate constants, as these rate constants could be fixed to various values with little effect on the likelihood values after refitting. The rate constants for Scheme X will be used in the discussion to describe how the channel gates at low and high Ca^{2+}.
Discussion
This study used detailed singlechannel analysis to examine the Ca^{2+}dependent gating of native BK channels in cultured rat skeletal muscle. We have extended previous studies by examining the effects of high Ca^{2+}_{i} to obtain critical information about mechanism when the gating is driven towards the fully liganded states. Maximum likelihood fitting together with comparisons of the observed and predicted 2D dwelltime distributions and dependency plots (the kinetic structure) were used to evaluate gating mechanisms.
Properties of Gating in High Ca^{2+}_{i}
For low to intermediate levels of Ca^{2+}_{i}, the gating was highly Ca^{2+}_{i} dependent, with a Hill coefficient of ∼3.5, within the range of 2–4 typically observed for BK channels (see introduction). In contrast, the gating was Ca^{2+} independent for high levels of Ca^{2+}_{i} (>100 μM). Increasing Ca^{2+}_{i} 8–10fold to 1,024 μM had little effect on P_{o} (∼0.97), the mean open and closed times, the 1 and 2D dwelltime distributions, and the dependency plots (Fig. 1, Fig. 2, Fig. 3, Fig. 4, Fig. 6, Fig. 7, Fig. 8, and Table 1 and Table 2). Estimates of the numbers of significant exponential components in the dwelltime distributions indicated that the channel entered at least three to four open and four to five closed states during normal activity at high Ca^{2+}_{i}, and that the estimated numbers of states did not change when Ca^{2+}_{i} was increased from 100 to 1,024 μM (Fig. 4). Significant dependencies (correlations) between the durations of adjacent open and closed intervals in the dependency plots (Fig. 6, Fig. 7, Fig. 8) indicated that transitions between the open and closed states occurred over at least two independent transition pathways (two or more gateway states) at high Ca^{2+}_{i}.
Thus, as a first approximation, gating at high Ca^{2+}_{i} involves Ca^{2+}independent transitions among at least three to four open and four to five closed states, with two or more independent transition pathways among the open and closed states. Models with these characteristics (Schemes XI–XIII) gave excellent descriptions of the kinetic structure for data limited to high Ca^{2+}_{i} (Table 4).
The MWC Model Is Inconsistent With Gating from Low through High Ca^{2+}_{i}
The MWC model (Monod et al. 1965) for allosteric proteins (Scheme III) predicts that high Ca^{2+}_{i} would drive the gating towards the two fully liganded open and closed states, resulting in a simple twostate gating mechanism at high Ca^{2+}_{i} with a single transition pathway between the states. These predictions of the MWC model are inconsistent with our observations at high Ca^{2+}_{i} of multiple open and closed states connected by two or more independent transition pathways. The MWC model or extensions of the MWC model (Schemes IV–VI) could not describe the gating from low to high Ca^{2+}_{i} (Fig. 5 and Table 3). Clearly, the MWC model can be rejected for the gating of BK channels in skeletal muscle.
Models with a Single Effective Gateway State at High Ca^{2+}_{i} Can Be Rejected
Since models based on the MWC model were inconsistent with the data at high Ca^{2+}_{i}, we examined more complex models based on Eigen 1968 general 35state allosteric model for tetrameric proteins and the 55state extension of Eigen's model (Cox et al. 1997a) shown in Scheme I. Since both the 35 and the 55state models have too many rate constants to estimate practically, we condensed these models to a 25state model (Scheme II) by assuming that the isoforms of each state were kinetically indistinguishable (as in Fersht 1985), and then examined models drawn from the 25state model. Scheme VII could describe the 1D dwelltime distributions from low to high Ca^{2+}_{i} (Fig. 5), but predicted no dependence between adjacent intervals at high Ca^{2+}_{i} (Fig. 9), in contrast to the significant dependencies observed in the experimental data (Fig. 6, Fig. 7, Fig. 8).
The lack of dependence for Scheme VII at high Ca^{2+}_{i} results because the high Ca^{2+}_{i} drives the gating towards the fully liganded states, as summarized by Scheme VIII, where there is a single effective transition pathway between the open and closed states. The observation that Scheme VII predicted no dependence at high Ca^{2+}_{i} suggests by analogy that the more complex 25–55state models (Scheme I and Scheme II), and all schemes in which there would be a single effective transition pathway between open and closed states at high Ca^{2+}_{i}, can be rejected.
A General TwoTiered Model Can Account for the Gating from Low to High Ca^{2+}_{i}
Schemes IIV can be classified as one tiered because the open and closed states can be contained within a single plane. Since the examined onetiered models were inadequate to describe the data at high Ca^{2+}_{i}, we considered what type of model would be required. If all the states in Scheme I and Scheme II are closed states, then there would be a sufficient number of fully liganded closed states at high Ca^{2+} to generate the four to five observed closed exponential components observed at high Ca^{2+}. If each of the closed states in Scheme I and Scheme II can make a direct transition to an open state, then this would give enough fully liganded open states to generate the three to four observed open components at high Ca^{2+}. If the gating is effectively confined to the fully liganded states at high Ca^{2+}_{i}, then the gating at high Ca^{2+}_{i} would be described by Scheme XV, where all the subunits are bound with Ca^{2+}. Each subunit in this scheme can exist in two conformational states, and a concerted conformational change of all subunits is required for opening. Scheme XV has five independent transition pathways between the fully liganded closed and open states (upper and lower tiers, respectively) that would allow dependence to be generated between open and closed intervals at high Ca^{2+}_{i}. Scheme XI, consistent with Scheme XV, gave excellent descriptions of the gating in high Ca^{2+}_{i} (Table 4). (Scheme XV)
Extending Scheme XV to include states with zero to three bound Ca^{2+} would give the general twotiered model described by Scheme IX, with 25 closed states on the upper tier and 25 open states on the lower tier. Scheme X, a reduced version of Fig. 9, could describe the kinetic structure of the channel from low to high Ca^{2+}_{i} (Fig. 10). Scheme X also generated singlechannel current records that closely mimicked, except for stochastic variation, the experimental current records, as can be seen by comparing the simulated (predicted) records in Fig. 13 to the experimental records in Fig. 1. The Ca^{2+}independent gating kinetics at high Ca^{2+}_{i} together with the few longer closed intervals are present in the simulated records.
Although the more complex Scheme IX was not tested directly, this general scheme should give an even better description of the gating than Scheme X, as Scheme X is contained within Scheme IX. Hence, Fig. 9 can serve as a working hypothesis for the Ca^{2+}dependent gating of BK channels. Horrigan and Aldrich (personal communication), based on analysis of macroscopic ionic and gating currents from mSlo, have also found evidence for twotiered gating mechanisms.
Gating in Scheme IX Occurs Among Five Subschemes that Differ in the Number of Bound Ca^{2+}
In the context of Scheme IX, it can be seen that, for any fixed number (0–4) of Ca^{2+} bound to the channel, the channel could gate among at least five open and five closed states. Thus, Scheme IX can be viewed as being comprised of five subschemes (one of which is shown in Scheme XV). The subunit conformations of the analogous states in each of the subschemes are the same, as are the transition pathways among the various open and closed states. Hence, both the conformations of the subunits and the connections among states are identical for each subscheme, independent of the number of bound Ca. Ca^{2+}_{i} acts by driving the gating from the subscheme comprised of the states with zero bound Ca^{2+} towards the subscheme comprised of the states with four bound Ca^{2+}. The binding of Ca^{2+} stabilizes the open states. The dynamics of this Ca^{2+}dependent shift will be presented in a later section.
Additional support for Scheme IX comes from the observation that BK channels gate at very low Ca^{2+}_{i} (Barrett et al. 1982; Pallotta 1985; Meera et al. 1996; Rothberg and Magleby 1996; Cui et al. 1997), among multiple closed (Talukder and Aldrich 1998) and open states (Nimigean, Rothberg, and Magleby, unpublished observations). The unliganded states in Scheme IX have five gateway states, suggesting that there would be a dependent relationship among the durations of open and closed intervals at zero Ca^{2+}_{i}, consistent with unpublished observations (Nimigean, Rothberg, and Magleby).
Why the Dependency Plots Appear Similar from Low to High Ca^{2+}_{i}
The saddle shape of the dependency plots (Fig. 6 and Fig. 7) reflects the inverse relationship between the durations of adjacent open and closed intervals. This inverse relationship suggests that more stable (longer duration) open states are effectively connected to less stable (briefer duration) closed states (McManus et al. 1985). In terms of Scheme IX, the general saddle shape of the dependency plots remains the same from low through high Ca^{2+}_{i} because the subschemes that are entered from low to high Ca^{2+}_{i} are comprised of the same numbers of open and closed states with the same subunit conformations and connections among the states. If the relative stability of the connected open and closed states depends mainly on the subunit conformations rather than on the numbers of bound Ca^{2+}_{i}, then the general shape of the dependency plots would remain the same from low through high Ca^{2+}_{i}.
Gating Dynamics at Low and High Ca^{2+}_{i}
Since Scheme X (drawn from Scheme IX) could account for the kinetic structure from low to high Ca^{2+}_{i}, Scheme X was examined to gain insight into how Ca^{2+} activates the channel. Fig. 14 presents the equilibrium occupancy, mean lifetime, and frequency of entry for each state for gating at low Ca^{2+}_{i} and high Ca^{2+}_{i} for Scheme X. Notice that bars of considerable height in these plots can have small values due to the logarithmic ordinates.
At low Ca^{2+}_{i} (Fig. 14, A–C), the channel readily gates among all of the states in Scheme X, as indicated by the frequencies of entry into each state (C), but spends most of its time in the closed states with zero (82.9%) and one (8.2%) bound Ca^{2+}_{i} (A). At high Ca^{2+}_{i}, the gating of the channel is effectively confined to the fully liganded states with four bound Ca^{2+}, where it spends 99.2% of its time. Little time (0.76%) is spent in the states with three bound Ca^{2+}, and negligible time (<0.04%) is spent in the states with 0, 1, or 2 bound Ca^{2+}_{i} (Fig. 14, D–F). The high P_{o} in high Ca^{2+} reflects that the channel spends 96.6% of its time in the open states with four bound Ca^{2+}, where Ca^{2+} stabilizes the open states. The flickers (brief closings) in the singlechannel current record arise mainly from transitions to the brief closed states C12– C16 at lower Ca^{2+}_{i} and C12 and C15 at high Ca^{2+}_{i}.
The Considered Schemes Are Only Approximations of the Gating Mechanism
As emphasized by Schoppa and Sigworth 1998, all kinetic modeling is by nature approximate. Although Scheme X could give good descriptions of the singlechannel kinetics from low to high Ca^{2+}_{i}, this model must be considered as only an approximation of the actual underlying gating mechanism. Scheme X with 17 states was drawn from the general 50state model described by Scheme IX. Additional open and closed states with zero and one bound Ca^{2+}_{i} would have to be added to Scheme X to allow gating in very low Ca^{2+}_{i}, and additional rows of open and closed states would have to be added to Scheme X to generate the observed numbers of open and closed components at high Ca^{2+}_{i}. The addition of these states would bring Scheme X to the general 50state model described by Scheme IX.
However, even the general 50state Scheme IX is a reduced model compared with what the actual gating mechanism is likely to be. Scheme IX excludes the isoforms of the various open and closed states. Including all the isoforms of each state (see Scheme I) would expand Scheme IX from 50 to 110 states.
BK channels may have additional closed states beyond the activation pathway (secondary states). Rothberg and Magleby 1998 found that models with such secondary states were consistent with the gating from low to intermediate Ca^{2+}_{i}, and Scheme XIII with such secondary states gave excellent descriptions of gating in high Ca^{2+}_{i} (Table 4). Such secondary states, if present, would add a third tier to Scheme IX and Scheme X consisting of closed states beyond the open states. Secondary states may also contribute to the gating of other K^{+} channels (Hoshi et al. 1994; Schoppa and Sigworth 1998).
The skeletal muscle BK channel often passes through a brief lifetime subconductance state upon opening and closing (Ferguson et al. 1993). It is not clear whether these subconductance states would arise from some of the states in Scheme IX and Scheme X, from isoforms of the states (see Scheme I) that were not included in Scheme IX and Scheme X, or whether additional states would have to be added to Scheme IX and Scheme X to account for the subconductance levels. For example, the concerted conformational changes that occur between closing and opening may occur in two steps, rather than the one indicated in Scheme IX and Scheme X (Schoppa and Sigworth 1998.
Changing Ca^{2+}_{i} more than two orders of magnitude, as was done in our experiments, would be expected to alter surface charge (Hille et al. 1975; Green and Andersen 1991) and hence gating (Moczydlowski et al. 1985). Since the considered gating mechanisms did not take surface charge into consideration, the question arises as to what effects this omission might have on the conclusions of our study. Changing the membrane potential ±20 mV (to mimic possible surface charge effects) did not alter the observed numbers of exponential components or the general saddle shape of the dependency plots (our unpublished observations). Since the rejection of previous models in favor of twotiered gating mechanisms in our study was based on factors that relate to the observed numbers of states and dependency at high Ca^{2+}_{i}, surface charge effects would be unlikely to alter the conclusion of twotiered gating mechanisms reached in this study, but could alter some of the rate constants.
Although Scheme X could describe many features of the data, this does not exclude the possibility that other rather different mechanisms might also account for the data. A model with fewer states than Scheme X, but with apparent saturation in Ca^{2+} binding rate could also give reasonable descriptions of the Ca^{2+} dependence of the singlechannel kinetics from low to high Ca^{2+}_{i} (Fig. 11 and Table 3). Nevertheless, we prefer the general twotiered approach based on Scheme IX to the more ad hoc saturation models, as twotiered models provide a means to account for the complexity of the gating in zero Ca^{2+}_{i} and are consistent with a tetrameric protein.
BK channels can gate in a number of different modes, with 96% of the intervals occurring during activity in the normal mode (McManus and Magleby 1988). The gating mechanisms developed in this study apply only to normal mode activity and would have to be expanded to account for activity in other modes. The gating mechanisms would also have to be expanded or modified to account for the effects of permeant ions on channel activity (Demo and Yellen 1992; Mienville and Clay 1996), the different gating properties of other BK channels such as dSlo (from Drosophila) where mean open times are relatively Ca^{2+} independent (Lagrutta et al. 1994; Moss et al. 1999), and the gating effects of various beta subunits not present in skeletal muscle (McManus et al. 1995; Dworetzky et al. 1996; Solaro et al. 1997; Tanaka et al. 1997; Nimigean and Magleby 1999; Waller et al. 1999).
Voltage Dependence of Gating
The voltage dependence of BK channels is an intrinsic property of the channel and does not appear to arise through voltagedependent increases in Ca^{2+} binding (Cox et al. 1997b; Cui et al. 1997; Rothberg and Magleby 1996, Rothberg and Magleby 1999; Stefani et al. 1997; Diaz et al. 1998). BK channels show homology to the superfamily of voltagedependent K^{+} channels, including an S4 voltage sensor (Atkinson et al. 1991; Adelman et al. 1992; Butler et al. 1993; Pallanck and Ganetzky 1994; Diaz et al. 1998). In the context of the general twotiered Scheme IX, depolarization could increase P_{o} through voltagedependent transitions of two general types. Depolarization could drive the concerted conformational changes that occur when states on the closed tier open to states on the open tier, or depolarization could drive the conformational changes of the individual subunits for transitions among states on each tier.
Recent preliminary observations on large multistate models using analysis of macroscopic ionic and gating currents (Horrigan and Aldrich 1998) and singlechannel recordings (Rothberg and Magleby, unpublished observations) suggest that the voltage dependence lies mainly in the Ca^{2+}independent closed–closed and open–open steps with the concerted closed–open step being less voltage dependent. Consistent with this hypothesis, gating charge movement can precede channel opening, and charge movement can also occur after the pore is open (Stefani et al. 1997). These observations suggest that it may only be necessary to add voltage dependence to the rate constants in the models that we have considered in order to account for the major features of the voltage dependence of the singlechannel kinetics. Nevertheless, Shaker K^{+} channels gate as if each subunit undergoes three transitions in sequence followed by two final concerted transitions for opening (Schoppa and Sigworth 1998), and movement of the S2 segment may precede that of the S4 segment (Cha and Bezanilla 1997), consistent with multiple conformations of each subunit. If the subunits in BK channels also have multiple conformations, then the potential numbers of states would be greatly increased over the models considered here.
Conclusion
This study develops gating mechanisms that can describe the Ca^{2+} dependence of the kinetic structure of BK channels from low to kinetically saturating levels of Ca^{2+}_{i}. These models are drawn from a general 50state twotiered model in which each closed state in the upper tier can make a direct transition to an open state in the lower tier. Our previous models that describe the Ca^{2+}dependent gating over more limited conditions are contained within the general 50state model. Thus, the 50state model serves to unify previous studies, and can provide a framework for further studies on mechanism through singlechannel analysis of gating at very low Ca^{2+}_{i} and of the voltage dependence of the gating.
Acknowledgments
This work was supported in part by grants from the National Institutes of Health (NS007044 and AR32805) and the Muscular Dystrophy Association.
Note added in Proof. Forthcoming kinetic studies lend additional support to 9. Gating currents and macroscopic ionic currents recorded from mSlo BK channels in the absence of Ca^{2}+ are consistent with a model that is described by the unliganded states in 9. (Horrigan, F.T., J. Cui, and R.W. Aldrich. 1999. Allosteric voltagegating of potassium channels I: mSlo ionic currents in the absence of Ca^{2+}. J. Gen. Physiol. In press. Horrigan, F.T., and R.W. Aldrich. 1999. Allosteric voltagegating of potassium channels II: mSlo channel gating charge movement in the absence of Ca^{2+}. J. Gen. Physiol. In press.)
Footnotes

1used in this paper: 1D, 1dimensional; BK channel, largeconductance calciumactivated K+ channel; MWC, MonodWymanChangeux; NLR1000, normalized likelihood ratio

The online version of this article contains supplemental material.
 Submitted: 5 March 1999
 Revision requested 6 May 1999
 Accepted: 7 May 1999