Role of the ß1 Subunit in Large-Conductance Ca2+-activated K+ Channel Gating Energetics: Mechanisms of Enhanced Ca2+ Sensitivity

D.H. Coxa,b and R.W. Aldrichc
a Molecular Cardiology Research Institute, New England Medical Center
b Department of Neuroscience, Tufts University School of Medicine, Boston, Massachusetts 02111
c Department of Molecular and Cellular Physiology and Howard Hughes Medical Institute, Stanford University, Stanford, California 94305

Correspondence to: D.H. Cox, New England Medical Center Hospitals, 750 Washington St., Box 7868, Boston, MA 02111. Fax:617-636-0576 E-mail:dcox01{at}emerald.tufts.edu.


  Abstract
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

Over the past few years, it has become clear that an important mechanism by which large-conductance Ca2+-activated K+ channel (BKCa) activity is regulated is the tissue-specific expression of auxiliary ß subunits. The first of these to be identified, ß1, is expressed predominately in smooth muscle and causes dramatic effects, increasing the apparent affinity of the channel for Ca2+ 10-fold at 0 mV, and shifting the range of voltages over which the channel activates -80 mV at 9.1 µM Ca2+. With this study, we address the question: which aspects of BKCa gating are altered by ß1 to bring about these effects: Ca2+ binding, voltage sensing, or the intrinsic energetics of channel opening? The approach we have taken is to express the ß1 subunit together with the BKCa {alpha} subunit in Xenopus oocytes, and then to compare ß1's steady state effects over a wide range of Ca2+ concentrations and membrane voltages to those predicted by allosteric models whose parameters have been altered to mimic changes in the aspects of gating listed above. The results of our analysis suggest that much of ß1's steady state effects can be accounted for by a reduction in the intrinsic energy the channel must overcome to open and a decrease in its voltage sensitivity, with little change in the affinity of the channel for Ca2+ when it is either open or closed. Interestingly, however, the small changes in Ca2+ binding affinity suggested by our analysis (Kc 7.4 µM -> 9.6 µM; Ko = 0.80 µM -> 0.65 µM) do appear to be functionally important. We also show that ß1 affects the mSlo conductance–voltage relation in the essential absence of Ca2+, shifting it +20 mV and reducing its apparent gating charge 38%, and we develop methods for distinguishing between alterations in Ca2+ binding and other aspects of BKCa channel gating that may be of general use.

Key Words: BKCa channel, mSlo, beta subunit, Ca2+ binding, allosteric model


  INTRODUCTION
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

Several mechanisms have evolved to tune the functional properties of ion channel proteins and thereby the electrical activity of various cell types. Prominent among them is the expression of auxiliary subunits that alter channel gating (Isom et al. 1994 ; Gurnett and Campbell 1996 ; Trimmer 1998 ). How auxiliary subunits accomplish this task is in most cases not understood, but finding answers to this question is important both for gaining a better understanding of the relationship between the structure of ion channels and how they function, and for providing information that may be used to design pharmaceuticals that can modulate the many physiological processes in which ion channels play important roles.

One ion channel whose activity is markedly affected by the expression of an auxiliary subunit is the large conductance Ca2+-activated potassium (BKCa)1 channel. In certain tissues, it appears that this channel is composed simply of four identical {alpha} subunits surrounding a central pore (Shen et al. 1994 ; Rothberg and Magleby 1999 ). Indeed, cloning and heterologous expression studies have revealed that channels composed of only {alpha} subunits display the essential properties of native BKCa channels: Ca2+ sensitivity, voltage sensitivity, K+ selectivity, and a large single-channel conductance (Atkinson et al. 1991 ; Adelman et al. 1992 ; Butler et al. 1993 ; Pallanck and Ganetzky 1994 ; Tseng-Crank et al. 1994 ). When BKCa channels where purified from bovine smooth muscle, however, two different subunits where isolated: a pore-forming {alpha} subunit (1,197 amino acids) and a smaller auxiliary ß subunit (191 amino acids) (Garcia-Calvo et al. 1994 ; Knaus et al. 1994b ; Giangiacomo et al. 1995 ). Both {alpha} and ß subunits were shown to be integral membrane proteins. Subsequent cloning and expression of these subunits revealed that while the {alpha} subunit forms a functional channel, coexpression of the ß subunit increases the sensitivity of the channel to Ca2+ at many membrane voltages (McManus et al. 1995 ; Wallner et al. 1995 ; Tseng-Crank et al. 1996 ). The presence of the ß subunit therefore appears to explain previous observations that indicated that smooth muscle BKCa channels are more sensitive to Ca2+ than are BKCa channels in many other tissues (McManus 1991 ; Tanaka et al. 1997 ).

Recently, additional BKCa ß subunits have been discovered and characterized (Ramanathan et al. 1999 ; Wallner et al. 1999 ; Xia et al. 1999 ; Brenner et al. 2000 ; Meera et al. 2000 ; Weiger et al. 2000 ). Although structurally similar, these novel ßs each alter BKCa channel gating in distinct ways, and they have different tissue distributions. They thus appear to have evolved to provide tissue-specific tuning of BKCa channel function, altering such properties as Ca2+ sensitivity, voltage–range-of-activation, and activation, deactivation, and inactivation kinetics. In this study, we have analyzed the functional effects of the original ß subunit, now designated ß1.

An intriguing property of BKCa channel gating even in the absence of ß1 is that in addition to being sensitive to internal Ca2+, these channels are also sensitive to changes in membrane voltage, and there is a synergy between these stimuli. One manifestation of this synergy is a shift in the channel's conductance–voltage relation toward more negative voltages as internal Ca2+ is raised (Marty 1981 ; Barrett et al. 1982 ; Latorre et al. 1982 ; Methfessel and Boheim 1982 ; Wong et al. 1982 ; Moczydlowski and Latorre 1983 ). Another is an apparent increase in the channel's Ca2+ affinity as the membrane voltage is made more positive (Markwardt and Isenberg 1992 ; Cui et al. 1997 ). Recently, we have developed allosteric models that mimic this synergy (Cox et al. 1997a ). They do so not by supposing that changes in voltage are directly affecting the binding of Ca2+ to the channel in a given conformation, but rather by supposing that the favorable decrease in free energy imparted to the channel's open conformation by one stimulus lessens the amount of work necessary for the other stimulus to activate the channel. The behavior of these models underscores the possibility that changes in the channel's apparent Ca2+ affinity can be brought about by means other than a direct effect on the channel's Ca2+ binding sites.

How then might the BKCa ß1 subunit increase the apparent affinity of the channel for Ca2+? At the single channel level, Nimigean and Magleby 1999 , Nimigean and Magleby 2000 have demonstrated that ß1 increases the time the mSlo channel occupies a set of predominately open bursting states, and that simply an increase in all of the channel's forward Ca2+ binding rates cannot account for this behavior. Thus, perhaps something more than just a change in Ca2+ binding affinity is occurring. On the other hand, their experiments did not rule out ß1-induced changes in Ca2+ unbinding rates. In fact, Meera et al. 1996 have found that substantial Ca2+ (~300 nM) must be present for human ß1 to shift the hSlo G-V relation toward more negative potentials. This would seem to argue for a direct effect of ß1 on Ca2+ binding or unbinding because when Ca2+ is not present, and thus its binding properties are irrelevant, ß1's steady state effects are apparently lost.

More generally, we may list at least three mechanisms by which the ß1 subunit may plausibly increase the BKCa channel's apparent Ca2+ affinity. It may alter the voltage sensitivity of the channel and thereby indirectly alter the channel's apparent Ca2+ affinity. It may change the energetics of conformational changes that are separate from both voltage sensing and Ca2+ binding, or it may in fact change the affinities of Ca2+ binding sites. To distinguish between these possibilities, we have examined the gating properties of the mSlo channel (the BKCa channel cloned from mouse; Butler et al. 1993 ) in the presence and absence of the ß1 subunit (bovine) over a wide range of conditions, and analyzed these changes in terms of allosteric models that represent our current understanding of BKCa channel gating (Cox et al. 1997a ; Cui et al. 1997 ; Horrigan and Aldrich 1999 ; Horrigan et al. 1999 ; Rothberg and Magleby 1999 ). The results of this analysis suggest that ß1 increases the Ca2+ sensitivity of the mSlo channel by both altering aspects of gating that do not involve Ca2+ binding, and by making subtle but important changes to the channel's Ca2+ binding affinity when it is either open or closed.


  MATERIALS AND METHODS
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

Channel Expression
The BKCa {alpha} subunit clone (mbr5) was kindly provided by Dr. Larry Salkoff (Butler et al. 1993 ), and the ß subunit clone (bovine ß1) was kindly provided by Dr. Owen McManus (McManus et al. 1995 ). Both clones were propagated in the Escherichia coli strain Top 10. In vitro transcription was performed with the mMessage mMachine kit with T3 or T7 RNA polymerase (Ambion Inc.). To record macroscopic currents ~0.05–0.5 ng of total cRNA was injected into Xenopus laevis oocytes 2–6 d before recording. ß1 and {alpha} cRNA were mixed in a ratio of 6:1 (ß/{alpha}) before injection. We found this ratio to be well above that necessary for saturation of ß1's effects.

Electrophysiology
Electrophysiological recordings were performed essentially as described previously (Cox et al. 1997b ). All recordings were done in the inside-out patch clamp configuration (Hamill et al. 1981 ). Patch pipettes were made of borosilicate glass (VWR Micropipettes). Their tips were coated with wax (Sticky Wax) and fire polished before use. Data were acquired using an Axopatch 200-A patch-clamp amplifier (Axon Instruments, Inc.) in the resistive feedback mode and a Macintosh-based computer system using "Pulse" acquisition software (HEKA Electronik) and the ITC-16 hardware interface (Instrutech). Records were digitized at 20-µs intervals (50,000 samples/s) and low pass filtered at 10 KHz with the Axopatch's four pole bessel filter. All experiments were carried out at 23°C. Under most conditions, before current records were analyzed and displayed (see, for example, Fig 6), capacity and leak currents were subtracted using a P/5 leak subtraction protocol with a holding potential of -120 mV and voltage steps opposite in polarity to those in the experimental protocol. With the ß1 subunit present, there was significant holding current at -120 mV with 39 or 74 µM [Ca2+]. Under these conditions, capacity and leak currents were measured by lowering [Ca2+] to 0.0005 µM and repeating the experimental protocol. Because of the limited voltage range, no channel activation was observed. The currents recorded with 0.0005 µM [Ca2+] were then used for capacity and leak current subtraction. As estimated from the amplitude and decay time constant of capacity currents, the series resistance (Rs) under our experimental conditions was 5–6 M{Omega}. Typically, 50–90% of Rs was compensated for using the Axopatch 200-A's Rs compensation circuitry. To increase the signal-to-noise ratio, typically three to six current series were taken consecutively under identical conditions and averaged before analysis and display. Maximum current amplitudes in each patch were typically 1–4 nA. In a series of 45 experiments, 23 without and 22 with the ß1 subunit, the means and standard deviations of the current amplitudes recorded with 74 µM [Ca2+] at +150 mV were 2.38 ± 1.19 nA without and 2.44 ± 1.00 nA with the ß1 subunit. These means are not statistically, significantly different (Students t test, P = 0.84); however, the 95% confidence interval calculated from these data for the difference between mean current amplitudes is 0.06 ± 0.66 nA.



View larger version (24K):
[in this window]
[in a new window]
 
Figure 1. The effects of Ca2+ on the mSlo G-V relation with and without the ß1 subunit. G-V relations determined without (A) and with (B) ß1 coexpression. [Ca2+] are as indicated in A. Each curve represents the average of between 6 and 15 experiments, as indicated in Table 1. Error bars represent SEM. Solid curves represent fits to the Boltzmann function {=}.Parameters of the fits are as follows, in order of increasing [Ca2+]. mSlo{alpha} V1/2 (mV) = 115.0, 101.9, 65.7, 51.8, 30.4, and 21.2; z = 1.25, 1.25, 1.27, 1.20, 1.10, and 1.09. mSlo{alpha} + ß 1 V1/2 (mV) = 79.2, 26.4, -9.6, -30.6, -90.1, and -95.6; z = 1.01, 1.18, 1.08, 0.99, 1.04, and 0.90. Dashed curves represent fits to Equation 1. Fit parameters are as follows: mSlo{alpha} L(0) = 1,263, Q = 1.18, KC = 6.00 µM, KO = 1.24 µM; mSlo{alpha} + ß1 L(0) 281, Q = 1.02, KC = 9.82 µM, KO = 0.82 µM. (C) Plots of V1/2 vs. [Ca2+] for mSlo{alpha} (•) and mSlo{alpha} + ß1 ({circ}) channels. Data are from Table 1. Error bars indicate SEM and are often smaller than the plot symbols.



View larger version (15K):
[in this window]
[in a new window]
 
Figure 2. Ca2+ dose–response curves for mSlo{alpha} (•) and mSlo{alpha} + ß1 ({circ}) at (A) +60 mV and (B) 0 mV. Data are from those shown in Fig 1A and Fig B. Smooth curves represent fits to the Hill equation G/Gmax = Amplitude/{1 + (Kd/[Ca])n}. Fit parameters are as follows: +60 mV mSlo{alpha} Kd = 3.84 µM, n = 1.5, Amplitude = 0.82; mSlo{alpha} + ß1 Kd = 1.20 µM, n = 3.5, Amplitude = 1.00: 0 mV mSlo{alpha} Kd = 32.8 µM, n = 1.3, Amplitude = 0.39; mSlo{alpha} + ß1 Kd = 3.43 µM, n = 1.9, Amplitude = 0.97.



View larger version (23K):
[in this window]
[in a new window]
 
Figure 3. Scheme I, two-tiered gating scheme. Those states in the top tier are designated closed. Those in the bottom tier are designated open. The central conformational change is voltage dependent with gating charge Q, and equilibrium constant L(0) = [closed]/[open]. KC1, KC2, KC3, and KC4 represent Ca2+ dissociation constants in the closed conformation. KO1, KO2, KO3, and KO4 represent Ca2+ dissociation constants in the open conformation. When KC1 = KC2 = KC3 = KC4 and KO1 = KO2 = KO3 = KO4, Scheme I represents a voltage-dependent version of the Monod-Wyman-Changeux model of allosteric proteins (Monod et al. 1965 ).



View larger version (37K):
[in this window]
[in a new window]
 
Figure 4. Effects of changes in various aspects of gating according to the voltage-dependent MWC model. (A) Control G-V simulations, with (from left to right) 0, 1, 10, and 100 µM [Ca2+]. Model parameters are as indicated. (B) Simulated effect of decreasing L(0). (C) Simulated effect of decreasing KC. (D) Simulated effect of decreasing Q. (E) Simulated effect of decreasing both Q and L(0). In A–E, the control curves of A are shown in gray for comparison. (F) Simulated V1/2 vs. [Ca2+] and (G) Q*V1/2 vs. [Ca2+] plots for the various conditions shown in A–E. Symbols are as indicated in A–E. Notice in G that only when the actual affinity of the model for Ca2+ is altered does the Q*V1/2 vs. [Ca2+] plot change shape.



View larger version (10K):
[in this window]
[in a new window]
 
Figure 5. (A) z*V1/2 vs. [Ca2+] plots for mSlo{alpha} ({blacksquare}) and mSlo{alpha} + ß1 ({square}). Each point represents the product of z and V1/2 values listed in Table 1. In B, the mSlo{alpha} + ß1 curve has been shifted vertically so that the two curves overlap at 9.1 µM. Errors bars indicate SEM and are often smaller than the plot symbols.



View larger version (21K):
[in this window]
[in a new window]
 
Figure 6. The ß1 subunit has effects on mSlo macroscopic current kinetics at subnanomolar [Ca2+]. Shown are families of current traces recorded in the absence (A) and presence (B) of ß1. For each trace, the membrane voltage was held at -50 mV and stepped to the indicated voltages. [Ca2+] equaled 0.5 nM. Superimposed on each trace are monoexponential fits. Time constants from these fits are plotted in C ({circ}). Also plotted are deactivation time constants (•) determined from monoexponential fits to currents elicited with a protocol that activated the channels with a depolarizing step, and then deactivated with steps to the indicated potentials (traces not shown).

Solutions
Recording solutions were composed of the following (mM): pipette solution: 80 KMeSO3, 60 N-methyl-glucamine-MeSO3, 20 HEPES, 2 KCl, 2 MgCl2, pH 7.20; internal solution: 80 KMeSO3, 60 N-methyl-glucamine MeSO3, 20 HEPES, 2 KCl, 1 HEDTA or 5 EGTA, and CaCl2 sufficient to give the appropriate free Ca2+ concentration, pH 7.20. EGTA (Sigma-Aldrich) was used as the Ca2+ buffer for solutions intended to contain <1 µM free Ca2+. HEDTA (Sigma-Aldrich) was used as the Ca2+ buffer for solutions intended to contain between 1 and 50 µM free Ca2+, and no Ca2+ chelator was used in solutions intended to contain >50 µM free Ca2+. 50 µM (+)-18-crown-6-tetracarboxylic acid (18C6TA) was added to all internal solutions to prevent Ba2+ block at high voltages.

The appropriate amount of total Ca2+ (100 mM CaCl2 standard solution; Orion Research Inc.) to add to the base internal solution containing 1 mM HEDTA or 5 mM EGTA to yield the desired free Ca2+ concentration was calculated using the program Max Chelator (Bers et al. 1994 ; available online at www.stanford.edu/~cpatton/maxc.html), and the proton and Ca2+ binding constants of Bers (supplied with the program) for pH 7.20, T = 23°C, and ionic strength = 0.15. The ability of 18C6TA to chelate Ca2+, as well as K+ and Ba2+, was also considered in these calculations using the following dissociation constants: Ca2+ 10-8 M (Dietrich 1985 ), K+ 3.3 x 10-6 M (Dietrich 1985 ), and Ba2+ 1.6 x 10-10 M (Diaz et al. 1996 ). For solutions intended to contain 0.5 µM free Ca2+ or more, free Ca2+ was measured with a Ca-sensitive electrode (Orion Research Inc.), and the measured value was reported. Free Ca2+ measurements were precise to within ~8%. Standard Ca2+ solutions for calibration of the Ca-sensitive electrode were prepared by adding known amounts (1, 10, 100, 1,000, and 10,000 µM) of standard CaCl2 solution to a solution identical to our base internal solution, except that no chelator was added. Endogenous Ca2+ in this solution was estimated from the deviation from linearity of the Ca-sensitive electrode's response at 10 µM added Ca2+, and was typically ~10 µM. Endogenous Ca2+ was then compensated for when making Ca2+-buffered solutions. When preparing solutions intended to contain <1 µM free Ca2+, the base internal solution was passed through a chelex 100 column before addition of chelator and Ca2+. This lowered the endogenous free Ca2+ to ~3 µM.

The concentration of Ca2+ in the solution bathing the cytoplasmic face of the patch was exchanged using a sewer pipe flow system (DAD 12) purchased from Adams & List Associates, Ltd. To minimize potential systematic variation, all comparisons between mSlo{alpha} and mSlo{alpha} + ß1 currents were done with the same experimental set up and the same solutions. The same batches of oocytes were also often but not always used.

Determination of G-V Curves and Fitting
At [Ca2+] >=600 nM, G-V relations were determined from the amplitude of tail currents 200 µs after repolarization to a fixed membrane potential (-80 mV) after voltage steps to the indicated test voltages. Each G-V relation was fitted with a Boltzmann function [G=]and normalized to the peak of the fit. At [Ca2+] < 600 nM, G-V relations were determined from tail current measurements 200 µs after repolarization to +40 mV. Since these curves did not reach saturation, the percentage of channels active at +200 mV in a given patch was estimated by comparison to current amplitudes recorded from the same patch with 9.1 µM [Ca2+], as described in the text (see Fig 7). Boltzmann fits were performed with "Igor Pro" graphing and curve fitting software (WaveMetrics Inc.) using the Levenberg-Marquardt algorithm to perform nonlinear least squares fits.



View larger version (20K):
[in this window]
[in a new window]
 
Figure 7. Effects of the ß1 subunit at subnanomolar [Ca2+]. (A) Current–voltage curves for mSlo{alpha} (•, n = 11) and mSlo{alpha} + ß1 ({circ}, n = 8) determined from patches superfused with a solution buffered to 0.5 nM [Ca2+]. Each current measurement was made at the end of a 50-ms voltage step to the indicated test voltage. Each curve has been normalized to their values at +200 mV. (B) G-V relations were determined as described in the text from the same data as used to generate the I-V plots in A. Smooth curves represent Boltzmann fits. Parameters of the fits are: mSlo{alpha} V1/2 = 180 mV, z = 0.98, mSlo{alpha} + ß1 V1/2 = 200 mV, z = 0.61. (C) In each experiment, data was acquired at both 9.1 µM (squares) and 0.5 nM [Ca2+] (circles). G-V plots for both concentrations are shown in C. In A–C, error bars indicate SEM.

Model fits to many G-V curves simultaneously were performed with Table-Curve 3D software, which employs an 80-bit Levenberg-Marquardt algorithm (Jandel Scientific) to minimize the sum of squares. When fitting to the 50-state model, several sets of initial parameters were tested. The resulting fits were usually similar, but not always identical. The best fit was then chosen as that which minimized the sum of squares.

Statistical Tests
Unpaired, two-tailed Student's t tests were used to test for significant differences between the mean V1/2 and z values listed in Table 1 and Table 2.


 
View this table:
[in this window]
[in a new window]
 
Table 1. G-V Parameters


 
View this table:
[in this window]
[in a new window]
 
Table 2. G-V Parameters

We tested whether there was a constant shift between mSlo{alpha} and mSlo{alpha} + ß1 mean z*V1/2 values over all [Ca2+] by fitting the data used to calculate the values plotted in Fig 5 A (below) with two statistical models. Model 1 assumed that the ß1 subunit had a constant affect on z*V1/2 regardless of [Ca2+]. Model 2 allowed the effects of ß1 on z*V1/2 to vary with [Ca2+]. The two models were fitted to the data using a maximum likelihood criteria. The linear models employed were as follows.

Model 1:(means assumed parallel). For this model, we define the following parameters. A0 represents the expected z*V1/2 value in the absence of the ß1 subunit at 0.99 µM Ca2+. A1 represents the effect of ß1 on z*V1/2, which in this model is assumed to be the same for all [Ca2+]. A2–A6 represent the additional effect of [Ca2+] over that observed at 0.99 µM for each measurement depending on the [Ca2+] at which a particular measurement was taken. Under this model, the expected value for a particular data point recorded, for example, with the ß1 subunit present and 9.1 µM [Ca2+], would be A0 (the expected value at 0.99 µM) + A1 (the effect of the ß1 subunit) + A4 (the additional effect of 9.1 µM [Ca2+] over that caused by 0.99 µM). Using statistical nomenclature, we may write the model then as: expected value at every point i = A0 + A1 I(ab_i = b) + A2 I(ca_i = 1.86) + A3 I(ca_i = 4.63) + A4 I(ca_i = 9.1) + A5 I(ca_i = 39) + A6 I(ca_i = 74).

Here, as follows from the definitions below, the terms not related to the conditions of a particular measurement evaluate to 0. The following definitions are applied. I() is the incremental function defined as [I(x = y) = 1, if x = y and I(x = y) = 0 if x is not equal to y]: ab_i = b if for the ith observation the ß1 subunit was present, ab_i = a if for the ith observation the ß1 subunit was absent, ca_i = 0.99 if for the ith observation [Ca2+] = 0.99, ca_i = 1.86 if for the ith observation [Ca2+] = 1.86, ca_i = 4.63 if for the ith observation [Ca2+] = 4.63, and so on.

Model 2: (means not assumed parallel). Following the above discussion, Model 2 may be written as: expected value of p_i = A0 + A1 I(ab_i = b) + A2 I(ca_i = 1.86) + A3 I(ca_i = 4.63) + A4 I(ca_i = 9.1) + A5 I(ca_i = 39) + A6 I(ca_i = 74) + A7 I(ab_i = b) I(ca_i = 1.86) + A8 I(ab_i = b) I(ca_i = 4.63) + A9 I(ab_i = b) I(ca_i = 9.1) + A10 I(ab_i = b) I(ca_i = 39) + A11 I(ab_i = b) I(ca_i = 74), where A0–A11 are parameters of the model. Here the effect of the ß1 subunit is not assumed to be the same at each [Ca2+], and so there is a separate parameter that represents the additional effect of the ß1 subunit (A7–A11) at each [Ca2+] above 0.99 µM. Under this model, the expected value for a particular data point recorded, for example, with the ß1 subunit present at 9.1 µM [Ca2+] would be A0 (the expected value at 0.99 µM) + A1 (the effect of the ß1 subunit at 0.99 µM [Ca2+]) + A4 (the additional effect of 9.1 µM [Ca2+] over that caused by 0.99 µM) + A9 (the additional effect of the ß1 subunit observed at 9.1 µM over that observed at 0.99 µM).

The likelihood of each individual observation was calculated from the expected value under the model at the appropriate condition and the estimated variance of the mean at each point. Model performance was compared using a likelihood ratio test that assumed normal errors, and allowed for heterogeneity of variance. The test was implemented with S-PLUS 2000 software (Mathsoft), and it gave a chi-square value of 38.236 on five degrees of freedom.


  RESULTS
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

We will refer to channels composed of {alpha} subunits alone as mSlo{alpha} and those composed of both {alpha} and ß1 subunits as mSlo{alpha} + ß1. The influence that internal Ca2+ and membrane voltage have on the steady state open probability of mSlo{alpha} is often displayed as a series of conductance-voltage curves determined at several internal Ca2+ concentrations ([Ca2+]) (Fig 1 A). The curves in Fig 1 A correspond to [Ca2+] ranging from ~1 to ~75 µM. As [Ca2+] is raised, the mSlo{alpha} G-V curve moves leftward along the voltage axis. When the ß1 subunit is coexpressed with {alpha}, a similar Ca2+-dependent shifting is observed (Fig 1 B). The G-V curves determined with ß1 present, however, differ from those determined with {alpha} alone in that each mSlo{alpha} + ß1 G-V curve lies to the left of the corresponding mSlo{alpha} curve. Over this [Ca2+] range, at every combination of [Ca2+] and voltage, the channels are more active when ß1 is present. This effect was first demonstrated by McManus et al. 1995 . Notice also that the mSlo{alpha} + ß1 G-V curves are more widely spaced than are the mSlo{alpha} curves. As shown in Fig 1 C, when half-maximal-activation–voltage (V1/2) is plotted as a function of [Ca2+], the wider spacing of the mSlo{alpha} + ß1 curves creates a function with a steeper slope.

Perhaps less apparent, the mSlo{alpha} + ß1 G-V curves are also somewhat less steep than are the mSlo{alpha} curves. In Fig 1A and Fig B, each G-V curve represents the average of many experiments (for n values, see Table 1). These curves were fitted with Boltzmann functions (solid curves), and the parameters of these fits are listed in the legend. At each [Ca2+], the Boltzmann parameter z, which is a measure of the voltage sensitivity of the channel, is smaller for mSlo{alpha} + ß1 than it is for mSlo{alpha}, on average by 14%. When mean parameter values from each experiment fitted individually are examined (Table 1), a similar decline in voltage sensitivity is observed at all [Ca2+] except 1.9 µM, although the small decline at 39 µM is not clearly statistically significant. Over all, however, looked at in this way, the average decline in z is 11%.

The effects of the ß1 subunit can also be seen from the point of view of changing [Ca2+] at constant voltage (Fig 2). Looked at in this way, at many voltages ß1 increases the channel's apparent Ca2+ affinity. For example at +60 mV the apparent Ca2+ affinity of mSlo{alpha} as judged by the concentration at which the channel's Ca2+ dose-response curve is half maximal ([Ca2+] 1/2), is 3.8 µM, while for mSlo{alpha} + ß1 it is 1.2 µM. At 0 mV, [Ca2+] 1/2 for mSlo{alpha} is 32.8 µM, while for mSlo{alpha} + ß1 it is 3.4 µM. Since the G-V curves in Fig 1 and the dose-response curves in Fig 2 represent different ways of looking at the same steady-state gating behavior, the effects of the ß1 subunit illustrated in these figures are manifestations of the same underlying process.

An Interpretation Using a Simple Model
What aspects of gating does the ß1 subunit alter to bring about the effects described above? Because a primary effect of ß1 is to alter the position of the mSlo G-V curve at many [Ca2+], to answer this question it seems important to understand what determines this curve's position. In general, the position of the G-V curve of any voltage-gated channel is determined by the free energy difference between open and closed states in the absence of an applied voltage and the degree to which an applied voltage alters this energy difference. This later factor is itself influenced by many factors including: the number of gating charges the channel has and how far each can move through the membrane's electric field, the degree to which these charges move in concert, and the degree to which the change in energy associated with each moving charge affects the energy difference between open and closed. Without some hypothesis about the physical situation, then, it is difficult to relate changes in the channel's G-V relation to changes in specific aspects of gating.

Recently, however, we have shown that models of the form of Scheme I (Fig 3) can mimic many aspects of mSlo{alpha} gating including: the leftward shifting of the mSlo{alpha} G-V curve as [Ca2+] is raised, an apparent increase in Ca2+ affinity as the membrane voltage is made more positive, near maximal activation with strong depolarizations at very low [Ca2+], and many properties of the kinetics of mSlo{alpha} macroscopic currents (Cui et al. 1997 ; Cox et al. 1997a ). The dashed lines in Fig 1 A represent a fit with a model of this form. While Scheme I cannot account for all mSlo{alpha} gating properties (see below and Cox et al. 1997a ; Stefani et al. 1997 ; Horrigan and Aldrich 1999 ; Horrigan et al. 1999 ; Rothberg and Magleby 1999 ), the success of this scheme in many respects suggests that it can be useful, at least initially, as a guide for interpreting the effects of the ß1 subunit in terms of changes in specific aspects of mSlo gating. We have discussed the properties of models of the form of Scheme I previously (Cox et al. 1997a ). We review them here in so far as they relate to the analysis that follows.

Scheme I supposes the following. The channel exists in two conformations, open and closed. The equilibrium between open and closed conformations is voltage dependent, with depolarization favoring opening. The channel is a homotetramer containing four Ca2+ binding sites, one on each subunit. Ca2+ can bind to both the open and closed conformations of the channel but on average binds better to the open conformation, and thereby promotes opening. The simplest form of Scheme I (wherein in a given conformation each binding site binds Ca2+ equivalently and independently) corresponds to a voltage-dependent version of the well known Monod-Wyman-Changeux (MWC) model of allosteric proteins (Monod et al. 1965 ). Its open probability as a function of [Ca2+] and voltage (V) is given by:

(1)

where L(0) represents the open-to-closed equilibrium constant in the absence of an applied voltage, a measure of the intrinsic energetics of channel opening. Q represents the equivalent gating charge associated with this equilibrium, KC the closed-conformation Ca2+ dissociation constant, and KO the open-conformation Ca2+ dissociation constant (F, R, and T have their usual meanings). Thus, four parameters determine the behavior of the model at equilibrium and, as illustrated in Fig 4, changes in any of these parameters can affect leftward G-V shifts at many [Ca2+]. Of interest to us, however, is which of these changes mimic the effects of the BKCa ß1 subunit.

The position of the voltage-dependent (VD) MWC model's G-V relation on the voltage axis, as indexed by its half-maximal activation voltage (V1/2) is given by:

(2)

and the derivative of Equation 2 with respect to [Ca2+] is:

(3)

Interestingly, this derivative is not a function of L(0). This means that according to the VD-MWC model alterations in the intrinsic energy difference between open and closed states will affect G-V position to the same extent at each [Ca2+], and if the ß1 subunit were only affecting this parameter, there would be no change in the slope of a plot of V1/2 vs. [Ca2+]. This is illustrated in Fig 4 B. Here, a simulated decrease in L(0) causes equivalent leftward G-V shifts at each [Ca2+] (compare Fig 4B with A), with the spacing between the curves maintained. On a plot of V1/2 vs. [Ca2+], this appears as a downward translation along the V1/2 axis with no change in slope (see Fig 4 F). In contrast to this prediction, we see in Fig 1 C that ß1 not only shifts the mSlo channel's V1/2 vs. [Ca2+] relation downward, but it also increases its steepness. That is, the magnitude of the ß1-induced G-V shift increases with increasing [Ca2+]. This suggests that the effects of ß1 are more complicated than just an alteration in the intrinsic energy difference between open and closed. The general downward movement of the channel's V1/2 vs. [Ca2+] relation, however, does suggest that this energy difference may be decreased by ß1 in combination with changes in other aspects of gating.

According to Equation 3, an increase in the steepness of a plot of V1/2 vs. [Ca2+] can come about either by changes in the channel's Ca2+ binding constants or by a decrease in its voltage sensitivity. An example of the former possibility is shown in Fig 4 C. Here, raising KC in our simulation produces leftward G-V shifts when [Ca2+] is present, with larger shifts at higher [Ca2+] (compare Fig 4C with A). At many voltages, this would appear as an increase in apparent Ca2+ affinity. At 0 mV, for example, the model channel is only 28% activated at 10 µM [Ca2+] when KC equals 10 µM, but 78% activated at the same [Ca2+] when KC equals 100 µM. Thus, by decreasing the channel's true Ca2+ binding affinity in the closed conformation, we can, in fact, create an apparent increase in Ca2+ affinity. This underscores the sometimes counterintuitive nature of even simple allosteric systems and highlights the need for a mechanistic and mathematical framework for interpreting alterations in their behavior. In this case, the increase in the model's apparent Ca2+ affinity with increasing KC arises because as KC and KO become farther apart, Ca2+ binding has a more powerful effect on the model's closed-to-open equilibrium. This can be seen from Equation 3, where dV1/2/d[Ca] depends always on the difference between KC and KO such that as these constants are separated, larger effects of Ca2+ are expected. Furthermore, because dV1/2/d[Ca] is a function of both KC and KO, any change in either of these constants will affect the steepness of the model's V1/2 vs. [Ca2+] relation.

Equation 3 also states that the steepness of the model's V1/2 vs. [Ca2+] relation is inversely related to its equivalent gating charge Q. This relation is therefore expected to become steeper as Q decreases and the model's G-V curves become more shallow (Cox et al. 1997a ; Cui, J., and R.W. Aldrich, manuscript in preparation). This is illustrated in Fig 4D and Fig F, where we have decreased Q in our simulation from 1.4 to 0.7 and observe an increase in G-V spacing (Fig 4 D), and thus an increase in V1/2 vs. [Ca2+] slope (Fig 4 F). Thus, some of the increased steepness of the mSlo channel's V1/2 vs. [Ca2+] relation may be due to the general decrease in voltage sensitivity evident upon ß1 coexpression (see Fig 1).

According to Scheme I (Fig 3), however, we can eliminate changes in voltage sensitivity (and thus G-V steepness) from consideration if we plot Q*V1/2 vs. [Ca2+] rather than V1/2 vs. [Ca2+] (Cox et al. 1997a ; Cui, J., and R.W. Aldrich, manuscript in preparation). This can be seen by rearranging Equation 2 to:

(4)

and again taking the derivative with respect to [Ca2+]:

(5)

This derivative is independent of both L(0) and Q and depends only on the model's Ca2+ binding constants. According to the VD-MWC model, then, a plot of Q*V1/2 vs. [Ca2+] will only change shape upon ß1 coexpression, if ß1 is actually altering Ca2+ binding site affinities. This is illustrated in Fig 4 G, where Q*V1/2 vs. [Ca2+] is plotted for the simulations in Fig 4A–E. Furthermore, although we have come to this conclusion based on the mathematics of the VD-MWC model (Equation 1 HREF="#FD2">Equation 2Equation 3Equation 4Equation 5), it is possible to show that this conclusion holds for all two-tiered models similar in form to Scheme I (Cox et al. 1997a ). It therefore does not depend on the simplifying, MWC assumption that in either channel conformation all Ca2+ binding steps are equivalent. Nor is it necessary to suppose that the channel has any particular number of Ca2+ binding sites. All that is required is a central voltage-dependent conformational change allosterically modulated by Ca2+ binding (Cox et al. 1997a ).

To apply this test to the effects of the BKCa ß1 subunit, we have multiplied each V1/2 value in Table 1 by its corresponding z value and plotted z*V1/2 vs. [Ca2+] for both mSlo{alpha} and mSlo{alpha} + ß 1 in Fig 5 A. Comparing the resulting relations to those in Fig 1 C, we see that, in fact, a good deal of the difference in steepness between the mSlo{alpha} V1/2 vs. [Ca2+] relation and the corresponding mSlo{alpha} + ß1 relation can be accounted for by ß1-induced changes in voltage sensitivity. That is, the two relations in Fig 5 A appear close to parallel. This result suggests that much of ß1's steady state effects can be accounted for by alterations in aspects of gating separate from Ca2+ binding. In particular, it appears that ß1 decreases the intrinsic energy difference between open and closed, thus shifting all G-V curves leftward and, in addition, ß1 decreases the channel's equivalent gating charge, thus expanding their spacing.

To facilitate comparison of the relations plotted in Fig 5 A, in B we have shifted the mSlo{alpha} + ß1 relation upward along the z*V1/2 axis, so that at 9.1 µM [Ca2+] the two relations have the same value. Comparing the two relations in this way we see their similar shape. It is also clear, however, that the mSlo{alpha} + ß1 relation has a somewhat higher value than expected at 0.99 µM [Ca2+], if the ß1 subunit where to have no effect on Ca2+ binding, and it has somewhat lower values than expected at 39 and 74 µM. Also, as the data are plotted in Fig 5 B, the standard errors of the mean for the z*V1/2 values at 0.99 and 39 µM [Ca2+] (indicated on the figure) are small compared with the differences between mSlo{alpha} and mSlo{alpha} + ß1 z*V1/2 values at these [Ca2+]. Thus, the differences in shape between the two curves in Fig 5 A may not be due to random variation in our experiments.

To examine this issue more rigorously, we performed a statistical test on the data used to calculate the values plotted in Fig 5 A. We used a likelihood ratio test and two statistical models. One assumed that ß1 has a constant effect on z*V1/2 as a function of [Ca2+]; that is, that the two sets of data in Fig 5 A are parallel. The other allowed for variation in the magnitude of ß1's effect with changing [Ca2+]. Comparison of the performance of the two models indicates that the two curves in Fig 5 A are not likely to be parallel (P < 0.0001), even when the variation about each mean is considered (for details of the statistical method used see MATERIALS AND METHODS). Thus, in addition to its effects on the energetics and voltage dependence of channel opening, the ß1 subunit may also have direct effects on Ca2+ binding.

To gain some quantitative estimate of how large these effects may be, we fit the mSlo{alpha} + ß1 G-V data in Fig 1 B with the VD-MWC model. The dashed lines in this figure represent the least squares best fit. Examining the parameters of this fit, and comparing them with those that best fit the G-V curves in Fig 1 A, we see that to accommodate the changes produced by ß1 the model, as expected, decreases L(0) (1,263 -> 281.3) and Q (1.2 -> 1.0). In addition, however, small changes in KO (1.24 -> 0.82 µM) and KC (6.0 -> 9.82 µM) are also observed, changes that together increase the ratio (KC/KO) and thereby also tend to expand G-V spacing. Interestingly, similar conclusions are also arrived at when these data are analyzed in terms of a more complex model of BKCa gating discussed below.

The Effects of the ß1 Subunit at Subnanomolar [Ca2+]
The above analysis suggests that although the ß1 subunit may alter Ca2+ binding-site affinities to a small degree (less than twofold in either conformation), much of its steady state effects are on aspects of gating that are separate from Ca2+ binding. If this is in fact the case, then some ß1 effects should be apparent at subnanomolar [Ca2+], where channel gating is Ca2+ independent (Meera et al. 1996 ; Cui et al. 1997 ). Recently, however, Meera et al. 1996 have found that at subnanomolar [Ca2+] there is essentially no change in the hSlo G-V relation upon ß1 coexpression. They concluded that at very low [Ca2+] there is "functional uncoupling" between {alpha} and ß subunits. They called this functional uncoupling a "Ca2+ switch." What this means in a physical sense is unclear. One possibility is that the two subunits become physically disassociated at low [Ca2+]. We know this is not the case, however, because the ß1 subunit has effects on the kinetics of mSlo macroscopic currents that are still evident at subnanomolar [Ca2+]. This is illustrated in Fig 6, and was also pointed out by Meera et al. 1996 (see their Figure 2 C). The other possibility seems to be that ß1's steady state effects are solely on Ca2+ binding such that at very low [Ca2+], where these sites are irrelevant, no change in channel gating is observed. This explanation, however, contradicts the above analysis, which suggested that, at micromolar [Ca2+], a large part of the steady state effects of ß1 were on aspects of gating separate from Ca2+ binding. Furthermore, a recent report indicates that the ß1 subunit does have effects on mSlo gating at subnanomolar [Ca2+] when examined at the single channel level (Nimigean and Magleby 2000 ).

To investigate this issue further, we examined mSlo{alpha} and mSlo{alpha} + ß1 currents at subnanomolar [Ca2+]. In Fig 7 A is plotted mean I-V relations determined from current families like those shown in Fig 6 recorded with 0.5 nM [Ca2+]. Each current measurement was made at the end of a 50-ms voltage step to the indicated test voltage. Each I-V curve was then normalized to its value at +200 mV, and the resulting curves were averaged. Two properties of these curves are of interest. First, even in the presence of ß1, it takes high voltages to activate the mSlo channel, at least +80 mV. And second, the two curves do not superimpose. The mSlo{alpha} + ß1 I-V curve is more shallow than the mSlo{alpha} curve. If ß1 were having no effect on the steady state activation of mSlo, then we would expect the shapes of their I-V curves to be the same. This is true regardless of the absolute degree to which the channels are activated by +200 mV. Thus, there does appear to be an effect of the ß1 subunit on mSlo steady state gating at subnanomolar Ca2+.

To examine this effect further, we determined mSlo{alpha} and mSlo{alpha} + ß1 G-V curves under these low [Ca2+] conditions. This was done by plotting tail current amplitude vs. test voltage for several current families like those shown in Fig 6. Each plot was then normalized to its value at +200 mV and then multiplied by the following ratio: the current recorded with 0.5 nM [Ca2+] at the end of a step to +200 mV, divided by the current recorded with 9.1 µM [Ca2+] under the same conditions. This procedure estimates the percentage of channels that are active at +200 mV by comparing currents recorded from the same patch at 9.1 µM and 0.5 nM [Ca2+]. It was necessary because at ~0.5 nM [Ca2+] neither G-V curve reached saturation by +200 mV, and higher voltages often caused the patches to break. Mean G-V curves determined in this way are shown in Fig 7 B. Most evident is a ß1-induced reduction in voltage sensitivity. Both curves have been fitted with Boltzmann functions with z values of 0.98 and 0.61 for mSlo{alpha} and mSlo{alpha} + ß1, respectively. Similar to the observations of Meera et al. 1996 , however, the V1/2 values of these fits are roughly similar (+180 mV for mSlo{propto} and +200 mV for mSlo{alpha} + ß1), although, statistically, the observed 20-mV difference is significant (P < 0.02).

At first glance, the relatively small and in fact right-shifting effect of the ß1 subunit on mSlo's half-maximal activation voltage at subnanomolar [Ca2+] seems to suggest that much of its ability to cause leftward G-V shifts at higher [Ca2+] is due to its effects on Ca2+ binding, and thus not evident in the essential absence of Ca2+. Upon reflection, however, it becomes clear that the reduction in apparent gating charge induced by ß1 at subnanomolar [Ca2+] would, if left unchecked, shift the mSlo{alpha} + ß1 V1/2 far rightward. Consider, for example, a channel behaving according to Scheme I with a single voltage-dependent step between closed and open. In this case, in the absence of Ca2+, the system reduces to two states, and we may write:

(6)

where {Delta}Gi represents the intrinsic standard free energy difference between open and closed. Thus, according to Equation 6, a 50% reduction in z would be expected to cause a twofold increase in V1/2, and a reduction in z from 0.98 to 0.61, as is observed upon ß1 coexpression, would be expected to shift V1/2 from its mSlo{alpha} value of +180 mV to +294 mV, a difference of 114 mV. We, however, see only a +20-mV rightward shift at subnanomolar [Ca2+] upon ß1 coexpression. This suggests that some compensating effect must be present. In fact, it suggests that in addition to reducing the channel's voltage sensitivity, ß1 decreases {Delta}Gi. This would act to oppose the right shift expected to accompany a decrease in voltage sensitivity, and thus it could maintain the mSlo{alpha} + ß1 V1/2 close to that of mSlo{alpha}, while allowing G-V steepness to decrease. Interestingly, however, according to Scheme I, over a range of [Ca2+], a decrease in {Delta}Gi [which amounts to a decrease in the parameter L(0) in Scheme I] would be expected to cause equivalent leftward G-V shifts at each [Ca2+], and a decrease in voltage sensitivity would be expected to expand G-V spacing. Thus, together these effects could in fact produce little change in V1/2 at subnanomolar [Ca2+], but large leftward G-V shifts at higher [Ca2+]. This is illustrated for the VD-MWC model in Fig 4 E, and it is qualitatively similar to what is observed in the data (Fig 1). Thus, our recordings at subnanomolar [Ca2+] demonstrate an effect of the ß1 subunit on the mSlo G-V relation under this condition, and, in fact, they qualitatively support the analysis above based on recordings with higher [Ca2+].

The ß1 Subunit Has Little Effect on the Affinity of the Channel for Ca2+ when It Is Open
Because Ca2+ binding promotes BKCa channel opening, thermodynamic principles dictate that the average affinity of the channel's binding sites for Ca2+ will be higher when it is open than when it is closed. This suggests that if we make the channel's open probability most energetically sensitive in the absence of Ca2+ (that is, bring it to V1/2) and then see how much Ca2+ it takes to just begin to activate the channels further, this critical [Ca2+] ([Ca2+]critical) will depend primarily on the affinity of the channel for Ca2+ when it is open (KO in the VD-MWC model). In fact, the mathematics of the VD-MWC model support this conclusion. For the VD-MWC model, the relationship between V1/2 and [Ca2+] is given by Equation 2. Because KC is necessarily significantly larger than KO, when [Ca2+] << KC, Equation 2 may be approximated by:

(7)

Equation 7 states that at [Ca2+] substantially less than KO (less than ~5% of KO, depending on Q), small changes in [Ca2+] will not affect V1/2. This is in fact what is observed experimentally (see below and Meera et al. 1996 ). As [Ca2+] is raised, however, at some point it will become appreciable relative to KO and V1/2 will start to move leftward. Furthermore, this [Ca2+]critical will not depend on L(0), which adds a constant term to V1/2, and it will depend only weakly on Q, which influences the magnitude of the shift but not the ratio of [Ca2+] to KO. This is illustrated graphically in Fig 8. Here we have calculated V1/2 vs. [Ca2+] plots for the VD-MWC model using parameters similar to those used to fit the mSlo{alpha} data in Fig 1 A. In Fig 8 A, the model's open-conformation Ca2+ dissociation constant (KO) is varied and we see large changes in [Ca2+]critical. In Fig 8 B, however, when KC is varied, we see very little effect on [Ca2+]critical. Similarly, in Fig 8 C we also see very little effect on [Ca2+]critical when L(0) or Q are varied. Thus, the [Ca2+] at which the VD-MWC model channel's G-V relation just starts to move leftward is determined almost exclusively by its affinity for Ca2+ when it is open. As shown in Fig 11 (below), we have also found this to be true for a more complex BKCa gating model.



View larger version (22K):
[in this window]
[in a new window]
 
Figure 8. The [Ca2+] at which the voltage-dependent MWC model's V1/2 begins to shift is determined primarily by its affinity in the open conformation. (A) V1/2 vs. [Ca2+] plots calculated from Equation 2 with L(0) = 1,000, Q = 1.0, KC = 10 µM, and KO varying as indicated. (B) V1/2 vs. [Ca2+] plots determined from Equation 2 with L(0) = 1,000, Q = 1.0, KO = 1 µM, and KC varying as indicated. (C) V1/2 vs. [Ca2+] plots determined from Equation 2 with KO = 1 µM and KC = 10 µM, and L(0) or Q as above, or varying as indicated.



View larger version (37K):
[in this window]
[in a new window]
 
Figure 9. [Ca2+]critical changes very little with ß1 coexpression. (A and B) Current traces recorded from macropatches in response to 400-ms ramps from -150 to +150 mV at the indicated [Ca2+] for mSlo{alpha} (A) and mSlo{alpha} + ß1 (B). All traces in A are from the same patch, as are all traces in B. In E, current amplitude at +150 mV recorded with each [Ca2+] relative to that recorded at 9.1 µM are plotted as a function of [Ca2+]. Plots from 6 mSlo{alpha} (•) and 11 mSlo{alpha} + ß1 ({circ}) experiments are shown. (C and D) G-V relations determined as described with reference to Fig 7 B for mSlo{alpha} (C) and mSlo{alpha} + ß1 (D). [Ca2+] are as indicated. Each curve represents the average of between 2 and 11 experiments as indicated in Table 2. Each curve has been fitted with a Boltzmann function with the following parameters, in order of increasing [Ca2+]. mSlo{alpha} V1/2 (mV) = 177, 176, 182, 174, 171, 110, 105, and 47; z = 0.90, 0.88, 0.93, 0.94, 0.93, 1.47, 1.57, and 1.28; mSlo{alpha} + ß1 V1/2 (mV) = 212, 211, 181, 197, 181, 71, 68, and -33; z = 0.62, 0.60, 0.61, 0.62, 0.55, 1.00, 1.00, and 0.98. In F are shown V1/2 vs. [Ca2+] plots for mSlo{alpha} (•) and mSlo{alpha} + ß1 ({circ}). Each point represents the mean of between 2 and 11 experiments, as indicated in Table 2. Error bars indicate SEM. Smooth curves represent fits to Equation 7. For the fits, the Q parameter was constrained to the mean value from analysis of all of our 0.5 nM G-V data. Fit parameters were as follows: mSlo{alpha} L(0) = 1,034, Q = 0.99, KO = 0.84 µM; mSlo{alpha} + ß1 L(0) = 141, Q = 0.64, KO = 0.68 µM.



View larger version (47K):
[in this window]
[in a new window]
 
Figure 10. (A) Scheme II, allosteric model of Horrigan et al. 1999 . Horizontal transitions represent voltage sensor motion in each of four subunits with KVC and KVO representing the forward microscopic equilibrium constants for these transitions for the closed and open channel, respectively. Vertical transitions represent the conformation change by which the channel opens. All transitions are hypothesized to be voltage dependent. (B) Allosteric model to account for both the Ca2+- and voltage-dependent properties of mSlo gating. This model represents the simplest combination of Schemes I and II. Transitions along the long horizontal axis represent Ca2+ binding and unbinding. Transitions along the short horizontal axis represent voltage sensor movement. Transitions from top to bottom represent channel opening. Implicit in this scheme and Equation 8 are the assumptions that voltage sensors and Ca2+ binding sites in each subunit are identical and act independently, and that voltage-sensor movement does not directly influence Ca2+ binding, and vice versa.



View larger version (24K):
[in this window]
[in a new window]
 
Figure 11. Simulations that demonstrate that [Ca2+]critical for Scheme III, as for Scheme I, depend primarily on the affinity of the model channel when it is open. (A) 50-state-model V1/2 vs. [Ca2+] plots as a function of KO. A series of G-V relations were simulated with Equation 8 at [Ca2+] ranging from 0 to 1,000 µM. V1/2 for each curve was then plotted as a function of [Ca2+]. KO was varied as indicated and the simulations were repeated. (B) Simulated V1/2 vs. [Ca2+] plots were generated as in A, except KC rather than KO was varied. Except where indicated, model parameters were: L(0) = 5e5, Q = 0.40, KC = 10 µM, KO = 1 µM, VhC = 135.9, D = 16.7, Z = 0.51.

While it is difficult to prove this result for all classes of reasonable models, it seems intuitively clear that the concentration at which a ligand begins to allosterically affect a central conformational change will be determined primarily by the affinity of the system when it is most sensitive to ligand given that two requirements are met. (a) In the absence of ligand, the central equilibrium constant is close to 1. For mSlo, we can force this to be true by adjusting the voltage to V1/2. And (b), there is a substantial difference in affinity between active and inactive forms. This requirement, however, is simply necessary for the ligand to have any appreciable effect at all. With the mSlo channel, therefore, if the ß1 subunit is altering the affinity of the channel for Ca2+ when it is open, we would expect to see this as a change in [Ca2+]critical.

We have taken two approaches to examining the effects of ß1 on [Ca2+]critical. To rapidly determine the effects of [Ca2+] on channel activity over a range of voltages, we applied voltage ramps from -150 to +150 mV at the following series of [Ca2+]: 0.5, 2, 10, 50, 100, and 660 nM, and 9.1 µM. Data from a patch expressing mSlo{alpha} channels are shown in Fig 9 A. As is evident, there is very little consistent effect of Ca2+ on channel activity at concentrations below 100 nM, a small effect at 100 nM, and then a large effect at 660 nM. The [Ca2+]critical for mSlo{alpha} therefore appears to be close to 100 nM. Interestingly, the same experiment produces a very similar pattern with mSlo{alpha} + ß1 (Fig 9 B), again no consistent effect below 100 nM, a small effect at 100 nM, and then a large effect at 660 nM. To quantify these results, we measured current amplitudes at +150 mV over the same [Ca2+] range for a series of mSlo{alpha} and mSlo{alpha} + ß1 patches, divided each patch's amplitudes by that recorded at 9.1 µM [Ca2+], and then plotted this ratio as a function of [Ca2+]. The results are shown in Fig 9 E. For a large number of both mSlo{alpha} and mSlo{alpha} + ß1 patches, [Ca2+]critical appears close to 100 nM, with a large effect of Ca2+ at 660 nM. (Fig 9)

To examine more rigorously the effects of the ß1 subunit on [Ca2+]critical, we recorded current families at each [Ca2+] and determined G-V relations as described above with reference to Fig 7. The results of Boltzmann fits to these data are listed in Table 2, and averaged G-V curves are shown in Fig 9C and Fig D. Again, as is consistent with published data (Meera et al. 1996 ), we see very little consistent effect of Ca2+ on both mSlo{alpha} and mSlo{alpha} + ß1 G-V relations at low nanomolar Ca2+, perhaps a small effect at 100 nM, and then a large effect at 660 nM [Ca2+]. Mean V1/2 values from a series of experiments are plotted for both mSlo{alpha} and mSlo{alpha} + ß1 in Fig 9 F. Although at 660 nM [Ca2+] mSlo{alpha} + ß1 responds with a larger leftward shift than does mSlo{alpha} (likely reflecting mSlo{alpha} + ß1's weaker voltage dependence), clearly the Ca2+ concentrations at which the two channels begin to respond to Ca2+ are essentially the same. This suggests that, when open, the affinity of the channel for Ca2+ is minimally affected by ß1, or perhaps not affected at all. To estimate KO for the two channel types, we fit the plots in Fig 9 F with Equation 7. These fits yielded KO estimates of 0.84 and 0.68 µM for mSlo{alpha} and mSlo{alpha} + ß1 channels, respectively. Thus, a small (19%) change in KO is suggested.

Analysis in Terms of a More Complex Model
Although Scheme I–based models mimic many properties of BKCa channel gating, there are some properties that they do not reproduce. In particular, gating current measurements indicate that the majority of mSlo{alpha}'s gating charge moves before channel opening rather with the same time course as opening, as is predicted by Scheme I (Stefani et al. 1997 ; Horrigan and Aldrich 1999 ). Also, mSlo{alpha} has a complex relationship between its deactivation time constant and membrane voltage that is inconsistent with this simple scheme (Cox et al. 1997a ; Horrigan et al. 1999 ). Recently, Horrigan et al. 1999 performed a detailed analysis of the macroscopic gating behavior of BKCa channels in the essential absence of Ca2+. Based on their results, they proposed that, rather than there being just a single voltage-dependent step between opened and closed, voltage sensors in each subunit move rapidly in response to changes in membrane voltage regardless of whether the channel is open or closed, and that this movement, in an allosteric manor, favors but is not required for channel opening (Horrigan et al. 1999 ; Horrigan and Aldrich 1999 ). This idea, coupled with the assumption of a weak voltage dependence associated with the central opening conformational change, is represented by Scheme II (Fig 10 A). Here, horizontal transitions represent voltage sensor movement with an X indicating a voltage sensor in its active conformation, and vertical transitions represent channel opening. Interestingly, Horrigan et al. 1999 showed that Scheme II can account well for the properties of mSlo{alpha} macroscopic gating and ionic currents in the essential absence of Ca2+, and they proposed that this scheme, coupled with the allosteric mechanism by which Ca2+ affects channel gating in Scheme I, may also account for BKCa gating at higher [Ca2+]. This idea in its simplest form is represented by Scheme III (Fig 10 B). In Scheme III, the subset of states represented by the 2 x 5 face at the front is equivalent to Scheme I, and the subset of states represented by the 2 x 5 face at the left edge is equivalent to Scheme II. Scheme III has also been recently proposed by Rothberg and Magleby 1999 to account for the gating of single skeletal muscle BKCa channels at high [Ca2+].

Although Scheme III has many states (50), it is based on the straightforward idea that both Ca2+ and voltage allosterically influence a central closed-to-open conformational change. The large number of states arises naturally as a property of gating systems regulated by two stimuli. Despite its many states, the equilibrium behavior of Scheme III is governed by only seven parameters: L(0), the open-to-closed equilibrium constant when no voltage sensors are active and no Ca2+ binding sites are occupied; Q, the gating charge associated with this equilibrium; Vhc, the voltage at which a single voltage sensor is half the time active when the channel is closed; Vho, the same for when the channel is open; Z, the equivalent gating charge associated with each voltage sensor's movement; and KC and KO, the open and closed Ca2+ dissociation constants, respectively. Popen for Scheme III is given by the following function of [Ca2+] and voltage:

(8)

In keeping with other work (Cox et al. 1997a ; Horrigan and Aldrich 1999 ), we may also define the ratio of Ca2+ dissociation constants as C (C = KC/KO), and the ratio of voltage sensor equilibrium constants when no voltage is applied as D [D = e(ZFVhc/RT)/e(ZFVho/RT)].

Interestingly, although it is not possible to write V1/2 as an explicit function of [Ca2+] for Scheme III, simulations indicate that, as with the VD-MWC model, the concentration at which Scheme III starts to respond to Ca2+ is also determined primarily by the affinity of the channel for Ca2+ when it is open. This is illustrated in Fig 11. Here we have constructed V1/2 vs. [Ca2+] plots from G-V curves simulated with a model of the form of Scheme III. The parameters used in these simulations were similar to those used below to fit the mSlo{alpha} data (see below). Varying KO (Fig 11 A) has a large effect on [Ca2+]critical, while varying KC (B) has very little effect. Thus, according to Schemes I and III, the fact that the mSlo channel's [Ca2+]critical changes very little upon ß1 coexpression suggests very little or no effect of ß1 on the open channel's Ca2+ binding affinity.

To further examine the effects of the ß1 subunit in terms of Scheme III, we fitted mSlo{alpha} and mSlo{alpha} + ß1 G-V data with Equation 8. For the mSlo{alpha} data, we restricted L(0), Vhc, D, Q, and Z to a small range of values determined by Horrigan et al. 1999 to reasonably fit mSlo{alpha} gating and ionic current data in the absence of Ca2+. We also restricted KO to be close (within 5%) to the value (0.84 µM) estimated for mSlo{alpha} from the [Ca2+]critical analysis described with reference to Fig 9. As shown in Fig 12 A, even with these restrictions, Scheme III can reproduce the shifting of the mSlo{alpha} G-V relation as a function of [Ca2+] fairly well. From the parameters of the fit (listed in Fig 12 A) KC for mSlo{alpha} is estimated to be 7.4 µM. When both Ca2+ dissociation constants were free to vary, a very similar best fit was found (not shown) with KC equal to 7.3 µM and KO equal to 0.72 µM.



View larger version (48K):
[in this window]
[in a new window]
 
Figure 12. 50-state model analysis. (A) mSlo{alpha} and (B) mSlo{alpha} + ß1 G-V relations over a series of [Ca2+] were fit with Equation 8 as described in the text. The resulting parameters are as indicated on the figure. Those values indicated in parentheses were not parameters of the fit, but rather calculated from the fit parameters. In A, the parameters other than those related to Ca2+ binding were constrained to be within the following ranges, which might reasonably fit the data of Horrigan et al. 1999 at subnanomolar [Ca2+] (Frank Horrigan, personal communication). L(0) 0.25–1 x 106, Q 0.35–0.45, Vhc 135–155, D = 10–20, Z = 0.51–0.59. K0 was constrained to be within 5% of the value estimated from [Ca2+]critical analysis 0.8–0.88 µM and KC was unconstrained. In B, KO was constrained to be within 5% of the value estimated from [Ca2+]critical analysis 0.65–0.72 µM, all other parameters were free to vary. In C–F, various parameters from the fit in B were restored to their values in A. In C, L(0) was restored. In D, the voltage-sensing parameters Q, Z, Vhc, and D were restored. In E, L(0), Q, Z, Vhc, and D were restored. In F, KO and KC were restored. In each panel, G-V relations at the following [Ca2+] are shown: from right to left, 0.0005, 0.99, 1.86, 4.63, 9.1, 39, and 74 µM. Error bars indicate SEM.

We next fit the G-V data acquired with the ß1 subunit present. In doing so, we constrained only KO, again keeping it close to (within 5%) the value (0.68 µM) estimated for mSlo{alpha} + ß1 from [Ca2+]critical analysis. As shown in Fig 12 B, over the [Ca2+] range ~1–75 µM, the 50-state model and the VD-MWC model fits are quite similar (compare with Fig 1 B, dashed curves). This is perhaps not surprising as the mechanism by which Ca2+ binding affects channel opening is the same for the two schemes. The more complex voltage-dependent gating mechanism of the 50-state model, however, also allows it to account for the prominent change in G-V steepness evident in Fig 12 B between 0.0005 and 0.99 µM [Ca2+]. This is something that is also evident without ß1 coexpression (see Fig 12 A and Cui et al. 1997 ), but is most pronounced when ß1 is present. Scheme I–based models cannot account for changes in G-V steepness with changing [Ca2+] because they contain all of their gating charge in the central close-to-open conformational change.

Most interesting, however, is the mechanism by which the 50-state model accounts for the ß1 subunit's steady state effects. As expected based on the preceding analysis, it decreases L(0), which reflects in part the intrinsic free energy difference between open and closed,2 and it decreases its voltage sensitivity (decreasing its gating charge z and altering Vho and Vhc), while changing its Ca2+ dissociation constants very little (KO moves from 0.8 to 0.65 µM, and KC moves from 7.4 to 9.6 µM). Also, we found that it was not possible to achieve any sort of reasonable fit to the mSlo{alpha} + ß1 G-V data if only KO and KC were allowed to vary freely while the parameters not related to Ca2+ binding were constrained to those used to fit the mSlo{alpha} data. Thus, analysis of our data in terms of Scheme III, as with Scheme I, suggests that although the ß1 subunit increases the mSlo channel's apparent affinity for Ca2+, much of this effect is via its action on aspects of gating that are separate from, but energetically coupled to, Ca2+ binding.

To examine the extent that changes in various aspects of gating account for ß1's effects, we restored various 50-state model parameters from the mSlo{alpha} + ß1 fit to their value in the mSlo{alpha} fit and recalculated the model G-V curves. We did this without redoing the fitting. The effects of such changes are shown in Fig 12C–F. Clearly, restoration of both the model's L(0) value (Fig 12 C) and its voltage-sensing parameters (D), as well as both together (E), severely disturbed the original mSlo{alpha} + ß1 fit. In fact, qualitatively, the effects of changing these aspects of gating are similar to what is predicted by Scheme I (Fig 4). Increasing L(0) shifts all G-V curves rightward (Fig 12 C), while increasing the model's voltage sensitivity increases G-V steepness and collapses their spacing (D). Thus, in these respects, Schemes I and III behave similarly.

We also restored just the model's Ca2+ dissociation constants to their values in the mSlo{alpha} fit. Interestingly, this also substantially disturbed the model's fit to the mSlo{alpha} + ß1 G-V curves (Fig 12 F). Not at very low [Ca2+], but as [Ca2+] is increased, the model's G-V spacing becomes too narrow and the fit becomes very poor. Also, when we did redo the fitting, but restricted KC and KO to the values used in the mSlo{alpha} fit, a good fit to the mSlo{alpha} + ß1 G-V data was not found. Thus, while both Schemes I and III suggest only small changes in the channel's Ca2+ dissociation constants with ß1 coexpression, the small changes that do occur appear to be important.

It is not clear, however, whether both KC and KO need to change to account for ß1's effects, as we were able to find reasonable fits to the mSlo{alpha} + ß1 G-V data when either one, but not both, of these parameters were held at their mSlo{alpha} G-V fit values. What appeared necessary to fit the data, however, was an increase in the ratio KC/KO (Fig 12 C) from ~9.3 for mSlo{alpha} to ~14.8 for mSlo{alpha} + ß1. Also, the standard errors of the fit parameter KC for the fits shown in Fig 12A and Fig B, (0.5 µM for the mSlo{alpha} fit and 0.6 µM for the mSlo{alpha} + ß1 fit, as determined by the curve fitting software) are small, relative to the KC estimates: <10%. This also suggests that the change in Fig 12 C that Scheme III uses to account for ß1's effects is required. Interestingly, this change corresponds to a change in the standard free energy difference between open and closed at saturating ligand of only 4.6 kJ/mol.


  DISCUSSION
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

BKCa channels are found in a wide variety of tissues where they are most often thought to provide feedback control over processes that involve membrane depolarization and a rise in intracellular Ca2+ (Latorre et al. 1989 ). One place they do this is at nerve terminals. Here, Ca2+ entry through nearby Ca2+ channels, in combination with membrane depolarization, results in rapid activation of BKCa currents that then contribute to the repolarization of the action potential (Robitaille and Charlton 1992 ; Robitaille et al. 1993a , Robitaille et al. 1993b ). Thus, neuronal BKCa channels are likely optimized for speed rather than sensitivity as Ca2+ concentrations in close proximity to Ca2+ channels are thought to rise as high as 100 µM (Roberts 1994 ), a concentration that according to Fig 1 is sufficient to achieve a Popen of ~0.5 at the peak of the action potential (approximately +20 mV). BKCa channels also play a central role in the control of smooth muscle contractility. In vascular smooth muscle, brief puffs or "sparks" of Ca2+ are released periodically from the sarcoplasmic reticulum. These sparks, which are thought to achieve Ca2+ concentrations near the plasma membrane in the tens of micromolar (Perez et al. 1999 ), activate nearby BKCa channels on the plasma membrane, and the ensuing hyperpolarization leads to muscle relaxation (Nelson et al. 1995 ; Nelson and Quayle 1995 ; Perez et al. 1999 ). In this tissue, however, the membrane voltage seldom exceeds -20 mV and typically rests near -40 mV (Nelson and Quayle 1995 ). Thus, for vascular smooth muscle BKCa channels to be similarly activated, they must be sensitive to Ca2+ at lower voltages than their neuronal counterparts.

To achieve this increase in sensitivity, nature apparently has evolved an auxiliary BKCa ß subunit (ß1) that shifts the BKCa channel's G-V relation -70 to -80 mV, into the physiological voltage range at concentrations of 5 or 10 µM. This shift increases the apparent Ca2+ affinity of these channels at many voltages. It has been the goal of this study to elucidate in a biophysical sense the mechanism by which the ß1 subunit exerts this powerful effect, and in particular to determine whether the apparent increase in Ca2+ affinity evident upon ß1 coexpression is due to a real increase in Ca2+ binding site affinity.

Perhaps the most straightforward approach to this question would be to measure the binding of Ca2+ to the channel. The high nonspecific binding associated with attempting to measure 45Ca2+ binding to a low affinity integral membrane protein, however, renders this approach, at present, technically unfeasible. We have therefore relied on electrophysiological measurements. These measurements, however, are necessarily indirect as they measure channel opening rather than Ca2+ binding, and therefore require interpretation. Binding measurements, if achieved, however, would also require some mechanistic framework for their interpretation, as such measurements may also be influenced by aspects of gating that are separate from, but linked to, Ca2+ binding (Wyman and Gill 1990 ).

To aid us in the interpretation of our electrophysiological data, we have relied on two models of BKCa channel gating that have been developed after extensive examination of the macroscopic gating behavior of the mSlo{alpha} channel. The first is Scheme I and the special case of Scheme I that corresponds to the voltage-dependent MWC model. While, as we have discussed (Cox et al. 1997a ; Horrigan and Aldrich 1999 ; Horrigan et al. 1999 ), this system is too simple to account for some of the voltage-dependent properties of BKCa channel gating, it can recreate the synergistic relationship between Ca2+ and voltage as they act together to activate the BKCa channel over a wide range of [Ca2+] and voltage. Furthermore, the VD-MWC model is sufficiently simple that it allows one to derive general principles about the ligand-dependent allosteric regulation of voltage-sensitive equilibria. The first of these principles is that an alteration in the intrinsic energy difference between open and closed ({Delta}Gi) can indeed cause large G-V shifts. This is perhaps intuitively clear, and readily understandable as simply due to a change in the energy the voltage-sensing mechanism needs to apply to the central conformational change to affect opening. An important characteristic of these shifts reflected by the VD-MWC model, however, is that as long as the channel's voltage sensitivity remains constant with [Ca2+](as is necessarily the case for Scheme I) the magnitude of the shift will not vary with Ca2+. That is, the spacing of the G-V curves will be maintained. Thus, on a V1/2 vs. [Ca2+] plot, a simple V1/2 offset will be observed. We, however, see a change in V1/2 vs. [Ca2+] slope as well as position in Fig 1 upon ß1 coexpression. This suggests that the ß1 subunit is doing more than just changing the equilibrium constant between open and closed.

The second general principle we may derive from the VD-MWC model is that the spacing of a BKCa channel's G-V relations determined over a series of [Ca2+] is likely to be inversely proportional to the voltage sensitivity of the channel (Cox et al. 1997a ; Cui, J., and R.W. Aldrich, manuscript in preparation). The more shallow the channel's voltage dependence, all other things being equal, the wider the spacing. We saw this in Fig 4 for Scheme I and in Fig 12 for Scheme III. This phenomenon is readily understandable as being due to the fact that as the channel's voltage sensitivity decreases, it takes a larger change in the transmembrane electrical field to apply the same amount of energy to the closed-to-open conformational change. Thus, an increase in the steepness of a channel's V1/2 vs. [Ca2+] plot is not necessarily indicative of an increase in Ca2+ binding affinity. For this reason, in our analysis we have multiplied the channel's effective gating valence at each [Ca2+] by its V1/2 and plotted this value as a function of [Ca2+]. As indicated by Equation 5, this results in a plot whose shape, for models of the form of Scheme I, is only sensitive to changes in Ca2+ binding. Applying this test to the effects of the ß1 subunit, we have found that the z*V1/2 vs. [Ca2+] plots for mSlo{alpha} and mSlo{alpha} + ß1 are more nearly parallel than are the V1/2 vs. [Ca2+] plots, thus suggesting that some of mSlo{alpha} + ß1's expanded G-V spacing is due to its general decrease in voltage sensitivity. Analysis of the two curves in Fig 5, however, indicates that there are differences in shape between the mSlo{alpha} and mSlo{alpha} + ß1 z*V1/2 vs. [Ca2+] plots, and thus some of the increased spacing of the mSlo{alpha} + ß1 G-V curves is likely due to effects of ß1 on Ca2+ binding.

A third principle we may derive from the VD-MWC model is that the spacing of the BKCa channel's G-V relation as a function of [Ca2+] will likely change if the model's Ca2+ dissociation constants in the open or closed configuration are altered. This is evident from Equation 3, which shows that {delta}V1/2/{delta}[Ca2+] depends in a complex way on both KC and KO, and this derivative is an even more complex function of KC and KO for Scheme III. In either case, however, the complexity of this relationship makes it difficult to make general predictions about how G-V spacing is likely to change in response to changes in Ca2+ binding constants without specifying the exact nature of the change. One important general principle of such ligand-gated voltage-sensitive systems, however, is that the power of a ligand to shift the channel's G-V relation along the voltage axis [that is, the maximum change in V1/2 that a ligand can produce ({Delta}V1/2max)], is dependent on the change in affinity that occurs upon channel opening, rather than the absolute affinities of the closed or open channel. In fact, for the VD-MWC model in the limit of saturating [Ca2+], we may write:

(9)

which states that {Delta}V1/2max varies as the logarithm of the ratio KC/KO. An important result of this fact is that large leftward G-V shifts at high [Ca2+] can occur if either KO is reduced or KC is increased. In the later case, the channel would in fact bind Ca2+ with lower affinity in the closed conformation and yet affect leftward G-V shifts. Furthermore, such shifts will increase the channel's Ca2+ sensitivity, particularly under conditions where multiple Ca2+ ions normally need to bind before the channel is appreciably open. Thus, an increase in a channel's apparent Ca2+ binding affinity can be the result of a true decrease in binding affinity. In fact, analysis of our data suggests that, among other things, the BKCa channel's KC/KO ratio must change from ~9.3 to ~14.9 to account for ß1's effects, and that simply a small increase in KC may well be responsible. Interestingly, this change in KC/KO corresponds to a change in the standard free energy difference between open and closed at saturating ligand of only 4.6 kJ/mol, less than the energy of a single hydrogen bond. For a weakly voltage-dependent channel such as mSlo (Q ~ 1.2), however, this small energy change converts to a 40-mV change in {Delta}V1/2max. Thus, even quite small changes in the KC/KO ratio can have substantial effects on the power of Ca2+ to shift the channel's G-V relation leftward and thus affect an apparent increase in Ca2+ affinity.

The general principles just described illustrate the complexity of systems containing binding coupled to conformational changes, and they underscore the futility of interpreting results of structural, modulatory, or pharmacological alterations in such cases without the aide of an appropriate and well-tested mechanistic model.

We have also examined the effects of the ß1 subunit on the mSlo G-V relation in the essential absence of Ca2+. A previous report suggested that there is a functional uncoupling between the hSlo channel and the human ß1 subunit under these conditions such that ß1 has no effect on the channel's G-V relation (Meera et al. 1996 ). Similar to this observation, we have found that the ability of ß1 to shift the mSlo G-V relation is greatly reduced at subnanomolar [Ca2+]. Inconsistent with this report, however, we have observed a 20-mV rightward G-V shift induced by the ß1 subunit at subnanomolar [Ca2+] and a 38% reduction in apparent gating charge. Furthermore, we show that these results can be accounted for if we suppose that coupled with the decrease in voltage sensitivity is a decrease in {Delta}Gºi (see Fig 4 E). Since this explanation is consistent with our overall analysis, we greatly favor it over an uncoupling mechanism, particularly as the kinetic effects of ß1 remain quite evident at subnanomolar [Ca2+], and effects of ß1 on mSlo gating have been observed under this condition at the single channel level (Nimigean and Magleby 2000 ). While species differences may account for the differences between our results and those of Meera et al. 1996 , it may also be that the limited voltage range used in their study obscured the effects we have observed. The kinetic effects of ß1 at subnanomolar [Ca2+], however, were clearly evident in both studies.

One advantage in studying the regulation of BKCa channel gating over other ligand-gated channels is the additional sensitivity of the BKCa channel to voltage. This allows the experimenter to manipulate the equilibrium between open and closed to a point where it is most sensitive to ligand binding, its V1/2, and then determine how much ligand it takes to begin to activate the channels further. We have argued that this critical [Ca2+] is reflective of the affinity of the channel when it is open. To examine [Ca2+]critical, we have looked for changes in BKCa currents at high voltages as we increased [Ca2+], and have found that, consistent with published data (Meera et al. 1996 ), both mSlo{alpha} and mSlo{alpha} + ß1 channels begin to respond to [Ca2+] at the same concentration, close to 100 nM. Thus, the ß1 subunit does not appear to substantially alter the mSlo channel's Ca2+ binding affinity when it is open.

In a similar experiment, we have also looked at the point where the mSlo{alpha} and mSlo{alpha} + ß1 channels' G-V relations just start to shift with [Ca2+]. We show, mathematically (Equation 7) for the VD-MWC model and graphically for the VD-MWC and the 50-state models (see Fig 8 and Fig 11), that this Ca2+ concentration depends almost exclusively on the affinity of the channel for Ca2+ when it is open. To be strictly correct, however, it should be stated that the channel's voltage sensitivity can also affect this critical concentration. This is because the size of the shift observed once the channel does start to respond to Ca2+ will be inversely related to voltage sensitivity, and this may influence the ability of the experimenter to detect the shift. We have accounted for changes in voltage sensitivity in our analysis, however, by restricting the value of z in our fits to the channel's V1/2 vs. [Ca2+] plot around [Ca2+]crtical (Fig 9 F) to the value that best fit the G-V at subnanomolar [Ca2+] (the data displayed in Fig 7 B). An alternative approach would be to look at the point where the channel's z*V1/2 vs. [Ca2+] just starts to respond to [Ca2+]. For the VD-MWC model, this point is independent of Q. Clearly, this type of analysis could also be used to determine whether alterations in BKCa channel function brought about by mutation, phosphorylation, or chemical modulation alter the affinity of the BKCa channel for Ca2+ when it is open.

To confirm and extend our analysis based on Scheme I, we have also fit our data with Scheme III. This scheme is identical in general form to that which Rothberg and Magleby 1999 proposed to best describe the Ca2+-dependent properties of single skeletal muscle BKCa channels at +30 mV. It is also a straightforward extension of Scheme II, which, after extensive investigation, Horrigan et al. 1999 found could account well for the gating and ionic current properties of mSlo{alpha} at subnanolmolar [Ca2+]. They proposed that combining Scheme II with the mechanism by which Ca2+ affects channel gating in the MWC model to produce Scheme III might account also for mSlo{alpha}'s Ca2+-dependent behavior. Remarkably, we have found that we could reproduce the Ca2+-dependent shifting of the mSlo{alpha} G-V relation over a wide range of [Ca2+], and as well a change in G-V slope from very low to moderate [Ca2+], while maintaining the voltage-sensing parameters (Q, Z, Vhc, and D) and L(0) of the 50-state model essentially as described by Horrigan et al. 1999 . In fact, this was also the case if we restricted KO to the value estimated from our [Ca2+]critical analysis such that only KC was free to vary. We consider this confirmation of the general theory that Scheme III represents, although clearly further complexities may occur that are not represented by Equation 8 and therefore are not included in the form of this scheme. Such complexities might include Ca2+ binding steps of differing affinities in either the open or closed conformation, or direct interactions between voltage sensors and Ca2+ binding sites. Even without such considerations, however, the added complexity of Scheme III's voltage-sensing mechanism makes it more difficult to derive general principles for its behavior than it is for Scheme I. Simulated changes in various aspects of gating, however, indicate that this system behaves qualitatively like Scheme I (Fig 12), thus validating our use of Scheme I in our preliminary analysis. This is perhaps not surprising since, although Scheme III's voltage-sensing mechanism is more complex than Scheme I's, Scheme III still represents a voltage-sensitive central conformational change allosterically regulated by ligand binding. Furthermore, we have found that both Schemes I and III suggest the same mechanism by which the ß1 subunit alters the mSlo channel's steady state gating behavior: a decrease in {Delta}Gºi, a reduction in gating charge, and a subtle but important change in Ca2+ binding affinity. Horrigan et al. 1999 performed a detailed experimental analysis of gating and ionic currents at subnanomolar [Ca2+], which allowed them to estimate Q, Z, VhC, D, and L(0) for the mSlo{alpha} channel. Without performing a similar analysis, we would not wish to rely on the estimates of these parameters suggested for the mSlo{alpha} + ß1 channel by the fit shown in Fig 12 B. We do think, however, that the broad [Ca2+] range we used in our experiments is sufficient to provide reasonable estimates of the mean Ca2+ dissociation constants for both mSlo{alpha} and mSlo{alpha} + ß1 under the conditions of our experiments. For mSlo{alpha}, these estimates are: KO = 0.8 µM and KC = 7.4 µM; and, for mSlo{alpha} + ß1, they are: KO = 0.65 µM and KC = 9.6 µM.

It should be pointed out, however, that the mean V1/2 and z values we report here for mSlo{alpha} are somewhat different from those we have published previously at similar [Ca2+] (Cui et al. 1997 ). The mean V1/2 values from our current data are on average 5–20 mV right-shifted compared with our previous report, and the mean z values are 10–30% smaller. We do not know the origin of these differences. We believe, however, that they are not due to random sampling, as the large number of experiments included in both studies render the estimates of mean V1/2 and z reported here, and previously, precise to a standard error of <5 mV (V1/2) and a charge of 0.1 (z). More likely, the differences are systematic, due to differences between the solutions in the two studies (we have reduced the K+ concentration in the internal and external solutions from 142 to 82 mM for the current study), the source or condition of the oocytes employed, or some other unknown aspects of the experimental set up. For the comparisons we have made in the present study, however, we do not believe such systematic variation to be an issue, as the mSlo{alpha} and mSlo{alpha} + ß1 data we report here were acquired with the same solutions, using the same experimental setup, in many cases from the same batches of oocytes. Thus, the experiments described here were specifically designed to compare mSlo{alpha} and mSlo{alpha} + ß1 currents, and their results are better compared with one another than with previous work.

The first ß1 subunit to be isolated and characterized was bovine (Knaus et al. 1994a ), and it is this subunit that we used in our experiments. During the course of this study, however, additional mammalian ß1 subunits were cloned, including those from mouse and human (Dworetzky et al. 1996 ; Tseng-Crank et al. 1996 ; Jiang et al. 1999 ). Mouse and human ß1 are ~83% identical to one another and between 81%(mouse) and 86% (human) identical to bovine ß1. Preliminary experiments we have done with mouse ß1 indicate that its effects on mSlo gating are very similar to bovine ß1 at micromolar [Ca2+], and that it produces a somewhat larger rightward G-V shift at subnanomolar [Ca2+]. This difference, however, is not likely to affect the main conclusions of this study.

Recently, Nimigean and Magleby 1999 , Nimigean and Magleby 2000 examined the effects of bovine ß1 on mSlo gating at the single-channel level. They concluded that the ß1 subunit exerts its effects mainly by increasing the time the channel spends in bursting states ~20-fold, and that this effect is evident even in the absence of [Ca2+]. If we suppose, as has been suggested by Rothberg and Magleby 1998 , that the flickers closed that occur during a burst are to states outside the main gating structure of the channel, then the observations of Nimigean and Magleby 1999 , Nimigean and Magleby 2000 fit with our observations. In this case, retention of the channel in bursting states could be accounted for by a reduction in {Delta}Gi, as our work suggests. Because plots of V1/2 vs. [Ca2+] are not parallel for the mSlo{alpha} and mSlo{alpha} + ß1 channels, however, our results suggests that the ß1 subunit also affects other aspects of gating.

In conclusion, results of our analysis of the mSlo{alpha} and mSlo{alpha} + ß1 G-V relation as a function of [Ca2+] suggests that the ß1 subunit exerts its steady state effects by decreasing the intrinsic energy difference between open and closed conformations, thereby shifting the mSlo G-V curve at all [Ca2+] leftward, and by decreasing the channel's voltage sensitivity, which expands G-V spacing as a function of [Ca2+]. In addition, subtle but important changes to one or more of the channels' Ca2+ dissociation constants are required, changes that increase the power of Ca2+ to shift the G-V relation upon binding, an effect that is most pronounced at high [Ca2+]. This work should help provide a framework by which to interpret experiments designed to determine in a molecular sense how the ß1 subunit makes these biophysical changes.


  Footnotes

1 Abbreviations used in this paper: BKCa, large-conductance Ca2+-activated K+ channel; MWC, Monod-Wyman-Changeux; VD-MWC, voltage-dependent MWC.
2 The parameter L(0) is defined differently in Schemes I and III (see Fig 3 and Fig 10 B, respectively). In Scheme I, L(0) represents the equilibrium constant between open and closed at 0 mV in the absence of [Ca2+]; whereas, in Scheme III, L(0) represents the equilibrium constant between open and closed at 0 mV in the absence of [Ca2+] if no voltage sensors are active. For Scheme III, the open-to-closed equilibrium constant in the absence of [Ca2+] at 0 mV is given by:


which thus depends on Z, Vhc, Vho, and L(0).


  Acknowledgements
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

We gratefully acknowledge Dr. Owen McManus for the generous gift of the bovine ß1 subunit clone, Dr. Frank Horrigan for helpful comments on the manuscript, Dr. Jianmin Cui for helpful discussions, and Dr. Michael Fay for his statistical analysis of the data in Fig 5.

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

Submitted: 7 March 2000
Revised: 25 July 2000
Accepted: 26 July 2000


  References
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

    Adelman, J.P., Shen, K.Z., Kavanaugh, M.P., Warren, R.A., Wu, Y.N., Lagrutta, A., Bond, C.T., and North, R.A. 1992. Calcium-activated potassium channels expressed from cloned complementary DNAs. Neuron. 9:209-216[Medline].

    Atkinson, N.S., Robertson, G.A., and Ganetzky, B. 1991. A component of calcium-activated potassium channels encoded by the Drosophila slo locus. Science. 253:551-555[Medline].

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

    Bers, D., Patton, C., and Nuccitelli, R. 1994. A practical guide to the preparation of Ca buffers. Methods Cell Biol. 40:3-29[Medline].

    Brenner, R., Jegla, T.J., Wickenden, A., Liu, Y., and Aldrich, R.W. 2000. Cloning and functional characterization of novel potassium channel beta subunits, hKCNMB3 and hKCNMB4. J. Biol. Chem. 275:6453-6461[Abstract/Free Full Text].

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

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

    Cox, D.H., Cui, J., and Aldrich, R.W. 1997b. Separation of gating properties from permeation and block in mslo large conductance Ca-activated K+ channels. J. Gen. Physiol. 109:633-646[Abstract/Free Full Text].

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

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

    Dietrich, B. 1985. Coordination chemistry of alkali and alkali-earth cations with macrocyclic ligands. J. Chem. Ed. 62:954-964.

    Dworetzky, S.I., Boissard, C.G., Lum, R.J., McKay, M.C., Post, M.D., Trojnacki, J.T., Chang, C.P., and Gribkoff, V.K. 1996. Phenotypic alteration of a human BK (hSlo) channel by hSlobeta subunit coexpression: changes in blocker sensitivity, activation/relaxation and inactivation kinetics, and protein kinase A modulation. J. Neurosci. 16:4543-4550[Abstract/Free Full Text].

    Garcia-Calvo, M., Knaus, H.G., McManus, O.B., Giangiacomo, K.M., Kaczorowski, G.J., and Garcia, M.L. 1994. Purification and reconstitution of the high-conductance, calcium-activated potassium channel from tracheal smooth muscle. J. Biol. Chem. 269:676-682[Abstract/Free Full Text].

    Giangiacomo, K.M., Garcia-Calvo, M., Hans-Gunther, K., Mullmann, T.J., Garcia, M.L., and McManus, O. 1995. Functional reconstitution of the large-conductance, calcium-activated potassium channel purified from bovine aortic smooth muscle. Biochemistry. 34:15849-15862[Medline].

    Gurnett, C.A., and Campbell, K.P. 1996. Transmembrane auxiliary subunits of voltage-dependent ion channels. J. Biol. Chem. 271:27975-27978[Free Full Text].

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

    Horrigan, F.T., and Aldrich, R.W. 1999. Allosteric voltage gating of potassium channels II. Mslo channel gating charge movement in the absence of Ca2+. J. Gen. Physiol. 114:305-336[Abstract/Free Full Text].

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

    Isom, L.L., De Jongh, K.S., and Catterall, W.A. 1994. Auxiliary subunits of voltage-gated ion channels. Neuron 12:1183-1194[Medline].

    Jiang, Z., Wallner, M., Meera, P., and Toro, L. 1999. Human and rodent MaxiK channel beta-subunit genes: cloning and characterization. Genomics. 55:57-67[Medline].

    Knaus, H.G., Folander, K., Garcia-Calvo, M., Garcia, M.L., Kaczorowski, G.J., Smith, M., and Swanson, R. 1994a. Primary sequence and immunological characterization of beta-subunit of high conductance Ca(2+)-activated K+ channel from smooth muscle. J. Biol. Chem. 269:17274-17278[Abstract/Free Full Text].

    Knaus, H.G., Garcia-Calvo, M., Kaczorowski, G.J., and Garcia, M.L. 1994b. Subunit composition of the high conductance calcium-activated potassium channel from smooth muscle, a representative of the mSlo and slowpoke family of potassium channels. J. Biol. Chem. 269:3921-3924[Abstract/Free Full Text].

    Latorre, R., Oberhauser, A., Labarca, P., and Alvarez, O. 1989. Varieties of calcium-activated potassium channels. Annu. Rev. Physiol. 51:385-399[Medline].

    Latorre, R., Vergara, C., and Hidalgo, C. 1982. Reconstitution in planar lipid bilayers of a Ca2+-dependent K+ channel from transverse tubule membranes isolated from rabbit skeletal muscle. Proc. Natl. Acad. Sci. USA. 79:805-809[Abstract].

    Markwardt, F., and Isenberg, G. 1992. Gating of maxi K+ channels studied by Ca2+ concentration jumps in excised inside-out multi-channel patches (myocytes from guinea pig urinary bladder). J. Gen. Physiol. 99:841-862[Abstract].

    Marty, A. 1981. Ca-dependent K channels with large unitary conductance in chromaffin cell membranes. Nature. 291:497-500[Medline].

    McManus, O.B. 1991. Calcium-activated potassium channels: regulation by calcium. J. Bioenerg. Biomembr. 23:537-560[Medline].

    McManus, O.B., Helms, L.M., Pallanck, L., Ganetzky, B., Swanson, R., and Leonard, R.J. 1995. Functional role of the beta subunit of high conductance calcium-activated potassium channels. Neuron. 14:645-650[Medline].

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

    Meera, P., Wallner, M., and Toro, L. 2000. A neuronal subunit (KCNMB4) makes the large conductance, voltage- and Ca2+-activated K+ channel resistant to charybdotoxin and iberiotoxin. Proc. Natl. Acad. Sci. USA. 97:5562-5567[Abstract/Free Full Text].

    Methfessel, C., and Boheim, G. 1982. The gating of single calcium-dependent potassium channels is described by an activation/blockade mechanism. Biophys. Struct. Mech. 9:35-60[Medline].

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

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

    Nelson, M.T., Cheng, H., Rubart, M., Santana, L.F., Bonev, A.D., Knot, H.J., and Lederer, W.J. 1995. Relaxation of arterial smooth muscle by calcium sparks. Science. 270:633-637[Abstract].

    Nelson, M.T., and Quayle, J.M. 1995. Physiological roles and properties of potassium channels in arterial smooth muscle. Am. J. Physiol. Cell Physiol. 268:C799-C822[Abstract/Free Full Text].

    Nimigean, C.M., and Magleby, K.L. 1999. The beta subunit increases the Ca2+ sensitivity of large conductance Ca2+-activated potassium channels by retaining the gating in the bursting states. J. Gen. Physiol. 113:425-440[Abstract/Free Full Text].

    Nimigean, C.M., and Magleby, K.L. 2000. Functional coupling of the ß1 subunit to the large conductance Ca2+-activated K+ channel in the absence of Ca2+: increased Ca2+ sensitivity from a Ca2+-independent mechanism. J. Gen. Physiol. 115:719-736[Abstract/Free Full Text].

    Pallanck, L., and Ganetzky, B. 1994. Cloning and characterization of human and mouse homologs of the Drosophila calcium-activated potassium channel gene, slowpoke. Hum. Mol. Genet. 3:1239-1243[Abstract].

    Perez, G.J., Bonev, A.D., Patlak, J.B., and Nelson, M.T. 1999. Functional coupling of ryanodine receptors to KCa channels in smooth muscle cells from rat cerebral arteries. J. Gen. Physiol. 113:229-238[Abstract/Free Full Text].

    Ramanathan, K., Michael, T.H., Jiang, G.J., Hiel, H., and Fuchs, P.A. 1999. A molecular mechanism for electrical tuning of cochlear hair cells. Science. 283:215-217[Abstract/Free Full Text].

    Roberts, W.M. 1994. Localization of calcium signals by a mobile calcium buffer in frog saccular hair cells. J. Neurosci. 14:3246-3262[Abstract].

    Robitaille, R., Adler, E.M., and Charlton, M.P. 1993a. Calcium channels and calcium-gated potassium channels at the frog neuromuscular junction. J. Physiol. (Paris). 87:15-24[Medline].

    Robitaille, R., and Charlton, M.P. 1992. Presynaptic calcium signals and transmitter release are modulated by calcium-activated potassium channels. J. Neurosci. 12:297-305[Abstract].

    Robitaille, R., Garcia, M.L., Kaczorowski, G.J., and Charlton, M.P. 1993b. Functional colocalization of calcium and calcium-gated potassium channels in control of transmitter release. Neuron. 11:645-655[Medline].

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

    Rothberg, B.S., and Magleby, K.L. 1999. Gating kinetics of single large-conductance Ca2+-activated K+ channels in high Ca2+ suggest a two-tiered allosteric gating mechanism. J. Gen. Physiol. 114:93-124[Abstract/Free Full Text].

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

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

    Tanaka, Y., Meera, P., Song, M., Knaus, H.G., and Toro, L. 1997. Molecular constituents of maxi KCa channels in human coronary smooth muscle: predominant alpha + beta subunit complexes. J. Physiol. 502:545-557[Abstract].

    Trimmer, J.S. 1998. Regulation of ion channel expression by cytoplasmic subunits. Curr. Opin. Neurobiol. 8:370-374[Medline].

    Tseng-Crank, J., Foster, C.D., Krause, J.D., Mertz, R., Godinot, N., DiChiara, T.J., and Reinhart, P.H. 1994. Cloning, expression, and distribution of functionally distinct Ca(2+)-activated K+ channel isoforms from human brain. Neuron. 13:1315-1330[Medline].

    Tseng-Crank, J., Godinot, N., Johansen, T.E., Ahring, P.K., Strobaek, D., Mertz, R., Foster, C.D., Olesen, S.P., and Reinhart, P.H. 1996. Cloning, expression, and distribution of a Ca(2+)-activated K+ channel beta-subunit from human brain. Proc. Natl. Acad. Sci. USA. 93:9200-9205[Abstract/Free Full Text].

    Wallner, M., Meera, P., Ottolia, M., Kaczorowski, G.J., Latorre, R., Garcia, M.L., Stefani, E., and Toro, L. 1995. Characterization of and modulation by a beta-subunit of a human maxi KCa channel cloned from myometrium. Receptors Channels. 3:185-199[Medline].

    Wallner, M., Meera, P., and Toro, L. 1999. Molecular basis of fast inactivation in voltage and Ca2+-activated K+ channels: a transmembrane beta-subunit homolog. Proc. Natl. Acad. Sci. USA. 96:4137-4142[Abstract/Free Full Text].

    Weiger, T.M., Holmqvist, M.H., Levitan, I.B., Clark, F.T., Sprague, S., Huang, W.J., Ge, P., Wang, C.C., Lawson, D., and Jurman, M.E. et al. 2000. A novel nervous system beta subunit that downregulates human large conductance calcium-dependent potassium channels. J. Neurosci. 20:3563-3570[Abstract/Free Full Text].

    Wong, B.S., Lecar, H., and Adler, M. 1982. Single calcium-dependent potassium channels in clonal anterior pituitary cells. Biophys. J. 39:313-317[Abstract].

    Wyman, J., and Gill, S.J. 1990. Binding and Linkage. Mill Valley, University Science Books, pp. 330 pp.

    Xia, X.M., Ding, J.P., and Lingle, C.J. 1999. Molecular basis for the inactivation of Ca2+- and voltage-dependent BK channels in adrenal chromaffin cells and rat insulinoma tumor cells. J. Neurosci. 19:5255-5264[Abstract/Free Full Text].