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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 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 conductancevoltage 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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 (
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 subunits surrounding a central pore (
subunits display the essential properties of native BKCa channels: Ca2+ sensitivity, voltage sensitivity, K+ selectivity, and a large single-channel conductance (
subunit (1,197 amino acids) and a smaller auxiliary ß subunit (191 amino acids) (
and ß subunits were shown to be integral membrane proteins. Subsequent cloning and expression of these subunits revealed that while the
subunit forms a functional channel, coexpression of the ß subunit increases the sensitivity of the channel to Ca2+ at many membrane voltages (
Recently, additional BKCa ß subunits have been discovered and characterized (
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 conductancevoltage relation toward more negative voltages as internal Ca2+ is raised (
How then might the BKCa ß1 subunit increase the apparent affinity of the channel for Ca2+? At the single channel level, 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;
![]() |
MATERIALS AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Channel Expression
The BKCa subunit clone (mbr5) was kindly provided by Dr. Larry Salkoff (
0.050.5 ng of total cRNA was injected into Xenopus laevis oocytes 26 d before recording. ß1 and
cRNA were mixed in a ratio of 6:1 (ß/
) 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 (. Typically, 5090% 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 14 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.
|
|
|
|
|
|
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 (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 and mSlo
+ ß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.
|
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.
|
|
We tested whether there was a constant shift between mSlo and mSlo
+ ß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+]. A2A6 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 A0A11 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 (A7A11) 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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
We will refer to channels composed of subunits alone as mSlo
and those composed of both
and ß1 subunits as mSlo
+ ß1. The influence that internal Ca2+ and membrane voltage have on the steady state open probability of mSlo
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
G-V curve moves leftward along the voltage axis. When the ß1 subunit is coexpressed with
, a similar Ca2+-dependent shifting is observed (Fig 1 B). The G-V curves determined with ß1 present, however, differ from those determined with
alone in that each mSlo
+ ß1 G-V curve lies to the left of the corresponding mSlo
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
+ ß1 G-V curves are more widely spaced than are the mSlo
curves. As shown in Fig 1 C, when half-maximal-activationvoltage (V1/2) is plotted as a function of [Ca2+], the wider spacing of the mSlo
+ ß1 curves creates a function with a steeper slope.
Perhaps less apparent, the mSlo + ß1 G-V curves are also somewhat less steep than are the mSlo
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
+ ß1 than it is for mSlo
, 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 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
+ ß1 it is 1.2 µM. At 0 mV, [Ca2+] 1/2 for mSlo
is 32.8 µM, while for mSlo
+ ß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 gating including: the leftward shifting of the mSlo
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
macroscopic currents (
gating properties (see below and
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 (
![]() |
(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 (
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+] (
![]() |
(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 4AE. 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 (
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 and mSlo
+ ß 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
V1/2 vs. [Ca2+] relation and the corresponding mSlo
+ ß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 + ß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
+ ß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
and mSlo
+ ß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 + ß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 ( 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
To investigate this issue further, we examined mSlo and mSlo
+ ß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
+ ß1 I-V curve is more shallow than the mSlo
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 and mSlo
+ ß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
and mSlo
+ ß1, respectively. Similar to the observations of
and +200 mV for mSlo
+ ß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 + ß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 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
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
Gi. This would act to oppose the right shift expected to accompany a decrease in voltage sensitivity, and thus it could maintain the mSlo
+ ß1 V1/2 close to that of mSlo
, while allowing G-V steepness to decrease. Interestingly, however, according to Scheme I, over a range of [Ca2+], a decrease in
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
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.
|
|
|
|
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 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
therefore appears to be close to 100 nM. Interestingly, the same experiment produces a very similar pattern with mSlo
+ ß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
and mSlo
+ ß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
and mSlo
+ ß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 ( and mSlo
+ ß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
and mSlo
+ ß1 in Fig 9 F. Although at 660 nM [Ca2+] mSlo
+ ß1 responds with a larger leftward shift than does mSlo
(likely reflecting mSlo
+ ß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
and mSlo
+ ß1 channels, respectively. Thus, a small (19%) change in KO is suggested.
Analysis in Terms of a More Complex Model
Although Scheme Ibased 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's gating charge moves before channel opening rather with the same time course as opening, as is predicted by Scheme I (
has a complex relationship between its deactivation time constant and membrane voltage that is inconsistent with this simple scheme (
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
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 (
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 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 and mSlo
+ ß1 G-V data with Equation 8. For the mSlo
data, we restricted L(0), Vhc, D, Q, and Z to a small range of values determined by
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
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
G-V relation as a function of [Ca2+] fairly well. From the parameters of the fit (listed in Fig 12 A) KC for mSlo
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.
|
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 + ß1 from [Ca2+]critical analysis. As shown in Fig 12 B, over the [Ca2+] range
175 µ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
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 + ß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
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 + ß1 fit to their value in the mSlo
fit and recalculated the model G-V curves. We did this without redoing the fitting. The effects of such changes are shown in Fig 12CF. 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
+ ß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 fit. Interestingly, this also substantially disturbed the model's fit to the mSlo
+ ß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
fit, a good fit to the mSlo
+ ß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 + ß1 G-V data when either one, but not both, of these parameters were held at their mSlo
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
to
14.8 for mSlo
+ ß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
fit and 0.6 µM for the mSlo
+ ß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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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+ (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 (
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 (
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 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 (
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 ( and mSlo
+ ß1 are more nearly parallel than are the V1/2 vs. [Ca2+] plots, thus suggesting that some of mSlo
+ ß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
and mSlo
+ ß1 z*V1/2 vs. [Ca2+] plots, and thus some of the increased spacing of the mSlo
+ ß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 V1/2/
[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 (
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 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
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 (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 (
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 ( and mSlo
+ ß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 and mSlo
+ ß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 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
's Ca2+-dependent behavior. Remarkably, we have found that we could reproduce the Ca2+-dependent shifting of the mSlo
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
Gºi, a reduction in gating charge, and a subtle but important change in Ca2+ binding affinity.
channel. Without performing a similar analysis, we would not wish to rely on the estimates of these parameters suggested for the mSlo
+ ß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
and mSlo
+ ß1 under the conditions of our experiments. For mSlo
, these estimates are: KO = 0.8 µM and KC = 7.4 µM; and, for mSlo
+ ß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 are somewhat different from those we have published previously at similar [Ca2+] (
and mSlo
+ ß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
and mSlo
+ ß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 (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, 20-fold, and that this effect is evident even in the absence of [Ca2+]. If we suppose, as has been suggested by
Gi, as our work suggests. Because plots of V1/2 vs. [Ca2+] are not parallel for the mSlo
and mSlo
+ ß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 and mSlo
+ ß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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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