Elimination of the BKCa Channel's High-Affinity Ca2+ Sensitivity

Lin Bao, Anne M. Rapin, Ericka C. Holmstrand and Daniel H. Cox

Molecular Cardiology Research Institute, New England Medical Center, and the Department of Neuroscience, Tufts University School of Medicine, Boston, MA 02111

Address correspondence to Daniel H. Cox, 750 Washington St., NEMC Hospitals, Box 7868, Boston, MA 02111. Fax: 617-636-0576; E-mail: dan.cox{at}tufts.edu


   ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We report here a combination of site-directed mutations that eliminate the high-affinity Ca2+ response of the large-conductance Ca2+-activated K+ channel (BKCa), leaving only a low-affinity response blocked by high concentrations of Mg2+. Mutations at two sites are required, the "Ca2+ bowl," which has been implicated previously in Ca2+ binding, and M513, at the end of the channel's seventh hydrophobic segment. Energetic analyses of mutations at these positions, alone and in combination, argue that the BKCa channel contains three types of Ca2+ binding sites, one of low affinity that is Mg2+ sensitive (as has been suggested previously) and two of higher affinity that have similar binding characteristics and contribute approximately equally to the power of Ca2+ to influence channel opening. Estimates of the binding characteristics of the BKCa channel's high-affinity Ca2+-binding sites are provided.

Key Words: mSlo • Slo • potassium channel • Ca2+ bowl • Ca2+ binding


   INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Large-conductance Ca2+-activated K+ channels (BKCa)* play an important feedback role in a number of physiological processes, most notably smooth-muscle contraction and neurosecretion, which involve changes in both membrane potential and intracellular Ca2+ concentration ([Ca2+]) (Latorre et al., 1989Go; McManus, 1991Go; Robitaille and Charlton, 1992Go; Robitaille et al., 1993Go; Nelson and Quayle, 1995Go). To accomplish this function they have evolved a unique sensitivity to both these signaling modalities that manifests as a Ca2+-dependent, leftward shifting of the channel's conductance-voltage (G-V) relation (Barrett et al., 1982Go; Methfessel and Boheim, 1982Go; Moczydlowski and Latorre, 1983Go). Biophysical studies have shed considerable light on the mechanism by which voltage-sensor activation influences channel opening (Cox et al., 1997Go; Stefani et al., 1997Go; Diaz et al., 1998Go; Rothberg and Magleby, 1998Go; Horrigan and Aldrich, 1999Go; Horrigan et al., 1999Go; Talukder and Aldrich, 2000Go). Less is known, however, about the BKCa channel's Ca2+-sensing mechanism and it is upon this mechanism that the present work is focused.

BKCa channels (Fig. 1) are structurally similar to K+ channels gated only by voltage (KV). Four identical subunits, known by their gene name slo, form a fully functional channel (Shen et al., 1994Go), and the amino-terminal third of the Slo protein contains several membrane-spanning regions, including a pore domain and a S4 helix that likely forms the channel's voltage sensors (Adelman et al., 1992Go; Butler et al., 1993Go; Diaz et al., 1998Go; Cui and Aldrich, 2000Go). Unlike KV channels, however, the Slo subunit contains a large intracellular domain that confers Ca2+ sensitivity on a voltage-gated structure (Wei et al., 1994Go; Schreiber et al., 1999Go; Moss and Magleby, 2001Go). The COOH-terminal part of this intracellular domain is termed the channel's "tail," while the portion of Slo that precedes the tail is termed the "core" (Fig. 1) (Wei et al., 1994Go). Neither core nor tail contain canonical Ca2+-binding motifs, and so the origin of the BKCa channel's Ca2+ sensitivity remains unresolved.



View larger version (14K):
[in this window]
[in a new window]
 
FIGURE 1. Schematic diagrams of the BKCa channel. (A) Diagram of the Slo subunit, four of which form a fully functional BKCa channel. Indicated are the core and tail domains, the pore helix (P), hydrophobic regions (S1-S10), the RCK domain, and the Ca2+-bowl. (B) Schematic diagram of a Slo tetramer.

 
There is reason to believe, however, that the tail is involved in Ca2+ sensing. Transfer of a large portion of the mouse Slo (mSlo) tail to a Ca2+-insensitive Slo homologue (Slo3) confers significant Ca2+-sensitivity (Schreiber et al., 1999Go), and 45Ca2+ binding to a tail polypeptide is reduced 56% by the mutation of five aspartates (Bian et al., 2001Go) that are among ten acidic residues in an unusually well-conserved tail domain termed the "Ca2+ bowl" (Schreiber and Salkoff, 1997Go). This domain contains 28 amino acids, and mutations here have been shown to cause rightward G-V shifts at constant [Ca2+], shifts similar to those observed with the wild-type channel when [Ca2+] is lowered (Schreiber and Salkoff, 1997Go; Bian et al., 2001Go; Braun and Sy, 2001Go). Thus, the Ca2+ bowl may form a Ca2+-binding site.

Arguing against an exclusive or even dominant role for the Ca2+ bowl in Ca2+ sensing, however, are the following observations: when the same aspartates were mutated in the full-length channel, the channel remained Ca2+ sensitive, its G-V relation still shifting well over 100 mV in response to changes in [Ca2+] (Bian et al., 2001Go); and this has been the case for all Ca2+ bowl mutations so far reported (Schreiber and Salkoff, 1997Go; Bian et al., 2001Go; Braun and Sy, 2001Go). No mutation in the Ca2+ bowl has been shown to substantially alter the concentration range over which the channel responds to Ca2+, and mutations outside the Ca2+ bowl have also been shown to cause rightward G-V shifts at constant [Ca2+] (Diaz et al., 1998Go; Cui and Aldrich, 2000Go; Braun and Sy, 2001Go). Thus, if there is a Ca2+-binding site in the Ca2+ bowl, then there are likely to be others outside this region as well (Schreiber and Salkoff, 1997Go; Schreiber et al., 1999Go; Bian et al., 2001Go).

To investigate further the role of the Ca2+ bowl in the Ca2+-sensing mechanism of the BKCa channel, we made a series of mutations in the mSlo Ca2+ bowl and analyzed their effects on the energetics of Ca2+-dependent gating. Our data argue that Ca2+ bowl mutations can eliminate the ability of a high-affinity Ca2+-binding site to promote channel opening, and therefore suggest that a Ca2+-binding site is present in this region. In addition, however, we have identified another site, at the end of the channel's seventh hydrophobic domain, that when mutated, also impairs Ca2+ sensing. And when the two sites are mutated together, all high-affinity Ca2+ response is lost. Thus, we propose that the BKCa channel contains two types of high-affinity Ca2+-binding sites, and in the experiments detailed below we estimate their binding characteristics.


   MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Channel Expression
The BKCa {alpha} subunit clone (mbr5) was propagated in the E. coli strains Top 10 or XL1-blue. In vitro transcription was performed with the "mMessage mMachine" kit with T3 RNA polymerase (Ambion). To record macroscopic currents ~0.5–50 ng of cRNA was injected into Xenopus laevis oocytes (stage IV-V) 2–6 d before recording.

Mutagenesis
All mutations were made with the Quick Change Site–directed mutagenesis kit (Stratagene), and mutations were identified by sequencing around the point of the mutation. Further, the coding regions of the following channels were sequenced in their entirety (wild-type mSlo [mbr5], {Delta}896–903, M513I, {Delta}896–903+M513I, {Delta}899–903+M513I).

Electrophysiology
Electrophysiological recordings were performed essentially as described previously (Cox and Aldrich, 2000Go). All recordings were done in the inside-out patch-clamp configuration (Hamill et al., 1981Go). Patch pipettes were made of borosilicate glass (VWR micropipettes), and had resistances of typically 1–2 Mohms in our recording solutions. Their tips were coated with wax (Sticky Wax) and fire polished before use. Data were acquired using an Axopatch 200B patch clamp amplifier (Axon Instruments, Inc.) or a list EPC-7 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 and low pass filtered at 10 KHz. All experiments were performed at room temperature, 22–24°C. Before current records were analyzed and displayed, 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.

Solutions
In general, recording solutions were composed of the following (in 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 1 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 containing <0.8 µM free [Ca2+]. HEDTA (Sigma-Aldrich) was used as the Ca2+ buffer for solutions containing between 0.8 and 10 µM free [Ca2+], and no Ca2+ chelator was used in solutions containing 100 µ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 1 mM EGTA to yield the desired free Ca2+ concentration was calculated using the program Max Chelator (Bers et al., 1994Go), which was downloaded from 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, 1985Go), K+ 3.3 x 10-6 M (Dietrich, 1985Go), Ba2+ 1.6 x 10-10 M (Diaz et al., 1996Go). 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 purchased from World Precision Instruments (10-8–10-2 M). Endogenous [Ca2+] in our internal solution before addition of Ca2+ chelator was estimated from the deviation from linearity of the Ca-sensitive electrode's response at 10 µM added [Ca2+], and was 16–20 µM. Endogenous [Ca2+] was then compensated for when making Ca2+-buffered solutions.

During our experiments, the solution bathing the cytoplasmic face of the patch was exchanged using a sewer pipe flow system (DAD 12) purchased from Adams and List Assoc. Ltd.

In the experiments of Fig. 5 , Mg2+ was often added to the internal solution. When this was the case, the following solutions were used. For the solution containing nominally 0 [Ca2+] and 10 mM [Mg2+], our base internal solution described above was used with the addition of 5 mM rather than 1 mM EGTA, and 12.1 mM MgCl2 (Sigma-Aldrich 4.9 M stock). For the solution containing 100 µM [Ca2+] and 10 mM [Mg2+] our base internal solution was used without added chelator, and 90 µM CaCl2 and 10 mM MgCl2 were added. For solutions containing 100 mM [Mg2+] our base internal was also used, except N-methyl-glucamine was not included, and then either 15 mM EGTA and 112 mM MgCl2 were added to create a solution containing nominally 0 [Ca2+] and 100 mM free [Mg2+], or 0 EGTA, 100 mM MgCl2, and 90 µM CaCl2 were added to create a solution containing 100 µM free [Ca2+] and 100 mM free [Mg2+].



View larger version (29K):
[in this window]
[in a new window]
 
FIGURE 5. The remaining Ca2+ response of M513I+{Delta}896–903 is blocked by Mg2+. (A) Top, G-V curves for the double mutant M513I+{Delta}896–903. Circles represent <3 nM [Ca2+]. Triangles represent 100 µM [Ca2+]. Open symbols indicate the absence of Mg2+. Grayed symbols indicate the presence of 10 mM [Mg2+]. Darkened symbols indicate the presence of 100 mM [Mg2+]. Each curve has been fitted with a Boltzmann function. Bottom, plots of the change in V1/2 as [Ca2+] is increased from <3 nM to 100 µM in the presence of either 0, 10, or 100 mM [Mg2+]. (B) Same experiment as in A except the wild-type channel was used. Mean V1/2 values were determined from Boltzmann fits to the data from each experiment individually.

 
Determination of G-V Curves and Fitting
Conductance-voltage (G-V) relations were determined from the amplitude of tail currents measured 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 = Gmax/[1 + e-zF{V - V1/2}/RT]) and normalized to the peak of the fit. All curve fitting was done with "Igor Pro" graphing and curve fitting software (WaveMetrics, Inc.) using the Levenberg-Marquardt algorithm to perform nonlinear least squares fits. Values in the text are given ± the standard error of the mean.


   RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Mutations at the Ca2+ Bowl Alter Ca2+ Binding
We began our investigation by mutating to alanine five contiguous residues (899–903) in the mSlo Ca2+ bowl, four potentially Ca2+-coordinating aspartates, and a proline (commonly a helix-breaking residue that might reasonably be important for maintaining the secondary structure of the surrounding region). Despite this fairly dramatic manipulation, however, the 899–903A mutant shares many functional characteristics with the wild-type channel (Fig. 2) . Except at 0.8 µM where the mutant's G-V curve is somewhat shallower, at most [Ca2+] the shapes of the wild-type and mutant curves are indistinguishable (Fig. 2 F and Table I) . The mutant and wild-type G-V curves are also in similar positions at 3 nM [Ca2+] (although the mutant's curve is somewhat right-shifted; Fig. 2 D and Table I). And perhaps most importantly, the mutant channel is still quite Ca2+ sensitive. Its G-V relation shifts over 150 mV as [Ca2+] is increased from 3 nM, a concentration too low to be sensed by the channel (Meera et al., 1996Go; Cui et al., 1997Go; Cox and Aldrich, 2000Go), to 100 µM, near saturating for the channel's high-affinity Ca2+-binding sites (Cox et al., 1997Go; Shi and Cui, 2001Go; Zhang et al., 2001Go). Thus, the 899–903A mutation does not eliminate the Ca2+ sensitivity of the mSlo channel, but it does reduce the power of Ca2+ to shift the mSlo G-V curve leftward along the voltage axis such that by 100 µM the mutant's G-V curve stands 52 mV to the right of the corresponding wild-type curve (Fig. 2, A, C, and D; Table I).



View larger version (34K):
[in this window]
[in a new window]
 
FIGURE 2. Ca2+ bowl mutations reduce mSlo's Ca2+ response. (A) G-V relations determined from inside-out Xenopus oocyte macropatches expressing the mSlo protein. [Ca2+] are as indicated on the figure. (B) Wild-type (top) and mutant (895–903A; bottom) mSlo current families recorded with voltage steps to between -40 and 130 mV with 100 µM internal [Ca2+]. Tail potentials are -80 mV; holding potentials are -120 mV (top) and -100 mV (bottom). (C) G-V relations determined from patches expressing the Ca2+ bowl mutant 899–903A. The residues mutated to alanine are indicated above the plot. (D) Half maximal activation voltage (V1/2) vs. [Ca2+] plots for the data in panels A, C, and E. The points plotted are average parameter values determined from experiments fitted individually with a Boltzmann function (see Table I for values). Symbols represent channel type as indicated on the figure. (E) G-V relations determined from patches expressing the Ca2+ bowl mutant 895–903A. (F) Apparent gating valence (z), determined from Boltzmann fits, plotted as a function of [Ca2+] (see Table I for values). Symbols represent channel type as in D. Error bars in this and subsequent figures represent standard error of the mean.

 

View this table:
[in this window]
[in a new window]
 
TABLE I Boltzmann-Fit Parameters for Wild-type and Mutant mSlo G-V Curvesa

 
These data indicate that the 899–903A mutant is similar to a class of Ca2+-bowl mutants termed by Schreiber and Salkoff (1997)Go to be "(+)-shifted," as at 100 µM, 10 µM, and 0.8 µM [Ca2+], the mutant's G-V relation is between 45 and 65 mV right or (+)-shifted relative to wild-type (see Fig. 2 D and Table I). However, because of the similarity of the mutant and wild-type curves at very low [Ca2+] (something also described for the [+]-shifted deletion mutant {Delta}897–899 by Schreiber and Salkoff [1997]Go), one might also consider the 899–903A channel's G-V curves, examined over a series of [Ca2+], as compressed along the voltage axis.

Because G-V compression can be attributed to an alteration in the voltage-sensing aspect of gating (Cox et al., 1997Go; Cox and Aldrich, 2000Go; Cui and Aldrich, 2000Go), it is reasonable to question whether the Ca2+-binding properties of the mSlo channel have been altered at all by the 899–903A mutation. However, the mutation's relatively small effects on G-V steepness and position at very low [Ca2+] and its large effect on G-V position at high [Ca2+] make a compelling case for a true change in the channel's Ca2+-binding properties. This can be evaluated quantitatively by interpreting these data in terms of the voltage-dependent Monod-Wyman-Changeaux model (VD-MWC), which supposes a central voltage-dependent closed-to-open conformational change that is allosterically regulated by Ca2+ binding (Cox et al., 1997Go; Cox and Aldrich, 2000Go) (see diagram in Fig. 6 B). According to this model, the position of the channel's G-V relation as indexed by its half-activation voltage (V1/2) is given by the following function of [Ca2+]

where KO and KC represent the model's Ca2+ dissociation constants when the channel is open or closed, respectively, L(0) represents the equilibrium constant of the central closed-to-open conformational change in the absence of an applied voltage, z represents the gating charge associated with this conformational change, n represents the number of Ca2+-binding sites that influence opening, and R, T, and F have their usual meanings. The amount that V1/2 will change with a given increase in [Ca2+] ({Delta}V1/2) is then

and the change in the product zFV1/2 is given by

(1)



View larger version (37K):
[in this window]
[in a new window]
 
FIGURE 6. Estimating {Delta}G(0)o-c based on the VD-MWC model. (A, C, and E) Wild-type and mutant G-V curves for the indicated mutants determined over an expanded series of [Ca2+] as indicated in C. (B) Diagram of the voltage-dependent MWC model (Cox et al., 1997Go). Horizontal transitions represent Ca2+ binding. Vertical transitions represent the conformational change by which the channel opens. KC and KO represent the model's closed and open-state Ca2+ dissociation constants, respectively. (D) V1/2 vs. [Ca2+] plots for the data in A, C, and E (see Table I for values). (F) zFV1/2 vs. [Ca2+] plots for the data in A, C, and E. Mean values of z and V1/2 were determined from fits to individual experiments separately and are listed in Table I. For the VD-MWC model, zFV1/2 = {Delta}G(0)o-c. Symbols indicate channel type as in D.

 
Eq. 1 is important here because for the VD-MWC and related models (see for example Cox et al., 1997Go), zFV1/2 is equal to the free energy difference between open and closed at 0 mV ({Delta}G(0)o-c) (Cui and Aldrich, 2000Go), and thus {Delta}zFV1/2 is equal to the change in this free-energy difference ({Delta}{Delta}G(0)o-c) as [Ca2+] is increased, a quantity, unlike {Delta}V1/2, that depends only on the Ca2+-binding properties of the model (Cox et al., 1997Go; Cox and Aldrich, 2000Go; Cui and Aldrich, 2000Go). For wild-type mSlo, {Delta}zFV1/2,(3 nM to 100 µM), evaluated from Boltzmann fits (see Table I), is 18.5 ± 0.66 kJ/mol, whereas for the 899–903A mutant it is 11.1 ± 0.60 kJ/mol. Thus, the 899–903A mutation appears to reduce the degree to which 100 µM [Ca2+] alters the channel's free energy difference between open and closed by 7.4 kJ/mol or 40%, and Eq. 1 suggests it must do so by reducing either the change in affinity that occurs at each Ca2+-binding site as the channel opens (that is the ratio KC/KO), or the number of binding sites the channel contains (n). This latter possibility is particularly interesting, as it would suggest that the 899–903A mutation knocks out a subset of the channels Ca2+ binding sites. Indeed, based on similar findings, Schreiber and Salkoff (1997)Go have argued that "(+)-shifted" mutants represent a complete loss of Ca2+-bowl function and that the channel's remaining Ca2+ sensitivity arises from Ca2+ binding elsewhere.

If this idea is correct, then mutations in the Ca2+ bowl more severe than 899–903A would not be expected to display a more severe phenotype, as under this hypothesis Ca2+ binding at the Ca2+ bowl is completely lost with 899–903A. To test this prediction, we mutated to alanine four additional residues in the heart of the Ca2+ bowl (895–903), three aspartates and one glutamine. Remarkably, despite the loss of 7 of 10 Ca2+ bowl acidic side chains, the 895–903A channel's equilibrium gating behavior is still appreciably Ca2+ sensitive, and very similar to that of the less severe 899–903A mutant (Fig. 2 E and Table I). Other than a small leftward (rather than rightward) G-V shift at 3 nM [Ca2+], neutralizing an additional three acidic and therefore potentially Ca2+-coordinating side chains had no further effect on G-V position as a function of [Ca2+] (Fig. 2 D and Table I). Thus, the smaller 899–903A mutation appears to be sufficient to completely eliminate Ca2+ binding at the Ca2+ bowl.

To examine, however, whether deleting rather than altering Ca2+-bowl residues might create a more severe phenotype, we also made the following deletion mutants: {Delta}895–903, {Delta}896–903, {Delta}898–903, and {Delta}899–903 deletions of 9, 8, 6, and 5 residues, respectively. Unfortunately, the {Delta}895–903 mutant did not produce stable currents. The others, however, did express well and displayed phenotypes that are similar to that of the alanine-substitution mutants described above: little effect at 3 nM and a fairly large reduction in the power of 100 µM [Ca2+] to shift the mSlo G-V relation leftward (Fig. 3) . In fact, the phenotype of the smaller deletion ({Delta}899–903) is virtually identical to that of 895–903A, and that of the larger deletions differs primarily at 100 µM [Ca2+], where the G-V curves for {Delta}896–903 and {Delta}898–903 are 74 and 79 mV right-shifted relative to wild-type, rather than 56 mV as is observed for 895–903A (Fig. 3 C and Table I). Despite this difference, the similar phenotypes of the five Ca2+-bowl mutants described above suggest the disruption of a common function, and we might tentatively conclude that the alanine-substitution mutants (899–903A, 895–903A) and the smaller deletion ({Delta}899–903) cause a near complete loss of function at the Ca2+ bowl, whereas the larger deletions ({Delta}896–903, {Delta}898–903) eliminate Ca2+-bowl function.



View larger version (32K):
[in this window]
[in a new window]
 
FIGURE 3. Many Ca2+ bowl mutants behave similarly. (A) Wild-type G-V relations determined at a series of [Ca2+] as in Fig. 2. (B) G-V relations for the mutant {Delta}899–903. The residues deleted are indicated above the plot with a line drawn through them. (C) V1/2 vs. [Ca2+] plots for the channels in A, B, D, and F (see Table I for values). (D) G-V relations for the mutant {Delta}898–903. (E) z vs. [Ca2+] plots for the channels in A, B, D, and F (see Table I for values). Symbols represent mutants as indicated in C. (F) G-V relations for the mutant {Delta}896–903. In A, B, D, and F symbols represent [Ca2+] as indicated in B.

 
Another Ca2+-binding Site
The above conclusion implies, as has been suggested (Schreiber and Salkoff, 1997Go; Bian et al., 2001Go), that high-affinity Ca2+-binding sites, not associated with the Ca2+ bowl, are responsible for the remaining Ca2+ sensitivity of the mutants described above. To search for these sites we examined a string of three acidic residues (EED 520–522) just past the channel's seventh hydrophobic segment (S7) in the intracellular part of the core (Wei et al., 1994Go; Jiang et al., 2001Go). We mutated all three to alanine, but the resulting mutant showed little deviation from wild-type behavior (unpublished data). However, when we mutated a nearby residue, M513, as shown in Fig. 4 B, the mutation created a channel with a reduced Ca2+ sensitivity similar to that of the Ca2+-bowl mutants described above (compare Fig. 4 B to Fig. 3 B, D, and F). Interestingly, M513 is at the end of S7, which places it 370 amino acids away from the Ca2+ bowl in the channel's primary sequence. Thus, it could be that this mutation affects a different Ca2+ binding site from that affected by mutations in the Ca2+ bowl but with perhaps similar binding characteristics. As we know little, however, about the tertiary structure of the BKCa channel, it could also be that mutations at the two sites alter the same binding site.



View larger version (29K):
[in this window]
[in a new window]
 
FIGURE 4. Double mutants that eliminate mSlo's high-affinity Ca2+ response. (A) Wild-type G-V relations. (B) G-V relations over a series of [Ca2+] for the mutant M513I. (C) V1/2 vs. [Ca2+] plots for the channels in A, B, D, F (see Table I for values). Symbols represent mutants as indicated on the figure. (D) G-V relations for the double mutant M513I+{Delta}896–903. (E) z vs. [Ca2+] plots (see Table I for values). Symbols represent mutants as indicated in C. (F) G-V relations for the double mutant M513I+{Delta}899–903. In A, B, D, and F symbols represent [Ca2+] as indicated in B.

 
To distinguish between these possibilities we made the double mutant M513I+{Delta}896–903. If each mutation acts independently, then their effects should be additive; if they act on the same process, this would not be expected (Goldstein et al., 1994Go; Stocker and Miller, 1994Go; Hidalgo and MacKinnon, 1995Go). G-V curves from the double mutant are shown in Fig. 4 D, and, as is evident, to a first approximation the mutations are additive. The two mutations reduce the mSlo channel's Ca2+ sensitivity considerably more than either mutation alone. In fact, the mutant channel shows almost no response to [Ca2+] below 10 µM, and it shows only a 45 ± 3.6 mV leftward G-V shift as [Ca2+] is increased from 3 nM to 100 µM (as compared with the 202 ± 5.8 mV shift of the wild-type channel, Table I). Thus, the two mutations together create a channel with only a low-affinity Ca2+ response, and interestingly the double mutant M513I+{Delta}899–903, which involves the smaller Ca2+-bowl deletion, displays an identical phenotype (Fig. 4 F).

Mg2+ Eliminates the Double Mutant's Ca2+ Sensitivity
BKCa channel gating can be modulated by intracellular Mg2+ (Golowasch et al., 1986Go) and millimolar concentrations of Mg2+ have recently been shown to activate the mSlo channel via a site in its core domain that also has a low affinity for Ca2+ (Shi and Cui, 2001Go; Zhang et al., 2001Go). As the double mutants described above show only a low-affinity Ca2+ response, we tested whether the remaining Ca2+ sensitivity of the M513I+{Delta}896–903 mutant is due to Ca2+ binding at this low-affinity, Mg2+-sensitive site. We examined the shift of the double mutant's G-V relation in response to 100 µM [Ca2+] in the absence and presence of 10 and 100 mM [Mg2+]. If Ca2+ binding at the Mg2+-sensitive site is responsible for the double mutant's low-affinity Ca2+ response, then 100 mM [Mg2+] would be expected to occupy this site and eliminate the effect of Ca2+, whereas 10 mM [Mg2+] would be expected to inhibit but not eliminate the effect of Ca2+ (Shi and Cui, 2001Go; Zhang et al., 2001Go).

As shown in Fig. 5 A, these predictions were born out in our data. In the presence of 10 mM [Mg2+] the {Delta}896–903+M513I channel's G-V shift, when exposed to 100 µM [Ca2+], is reduced from 45 ± 3.6 mV to 24 ± 3.4 mV, and in the presence of 100 mM [Mg2+] it is completely eliminated ({Delta}V1/2 = -2 ± 6.8 mV, see Fig. 5 A, bottom). Furthermore, when the same experiment is performed with the wild-type channel (Fig. 5 B), 100 mM [Mg2+] does not completely eliminate the effect of 100 µM [Ca2+], but rather reduces {Delta}V1/2 by an amount similar to its effect on the double mutant (Fig. 5 B bottom). These results argue that Mg2+ only effectively blocks the channel's low-affinity Ca2+ binding sites (as has been reported, Shi and Cui, 2001Go; Zhang et al., 2001Go), and that the weak Ca2+ response of the double mutant is due to low-affinity Ca2+ binding at Mg2+-sensitive sites. Furthermore, the high-affinity response of the wild-type channel must therefore arise from Ca2+ binding at sites whose functional effects are completely eliminated by the double mutant. This is remarkable, because it suggests the presence of three types of Ca2+-binding sites on the mSlo channel, one of low affinity whose effects can be blocked by high concentrations of Mg2+, and two of higher affinity, one whose effects are eliminated by mutations in the Ca2+ bowl, and another whose effects are eliminated by mutation at M513.

Estimating the Ca2+-binding Characteristics of the Channel's High-Affinity Sites
The mutations described above must cause complete loss of function at their respective sites; otherwise, the double mutant would respond to Ca2+ even after its low-affinity sites are occupied by Mg2+. Thus, they provide an opportunity to estimate the binding characteristics of the sites they disrupt. To do this, however, requires extracting from the wild-type and mutant G-V data the degree to which Ca2+ binding is affecting {Delta}G(0)o-c, the free energy difference between open and closed in the absence of an applied voltage. That is, we would like to plot {Delta}G(0)o-c rather than V1/2 as a function of [Ca2+], so as to examine a physically meaningful Ca2+-dependent quantity that is independent of the mechanism by which the channel responds to voltage, and yet readily interpretable in terms of the binding properties of the channel's Ca2+-binding sites.

A strategy which has been used to do this in the past is to fit each G-V relation with a Boltzmann function and from the parameters of the fit estimate {Delta}G(0)o-c as {Delta}G(0)o-c = zFV1/2 (Cox and Aldrich, 2000Go; Cui and Aldrich, 2000Go; Jiang et al., 2001Go; Shi and Cui, 2001Go), and in fact this was the strategy we used above in examining the 899–903A mutant. In effect, it assumes a two-state model for the channel's voltage-sensing mechanism, and then on the basis of this model estimates the effect of Ca2+ binding on the free energy difference that the applied voltage must overcome to open the channel. As discussed above, for the voltage-dependent MWC model (Fig. 6 B) this strategy is strictly correct (Cox et al., 1997Go; Cox and Aldrich, 2000Go; Cui and Aldrich, 2000Go).

Results of such an analysis for the wild-type channel and the {Delta}896–903 and M513I mutants are plotted in Fig. 6 F, together with standard V1/2 vs. [Ca2+] plots (Fig. 6 D). New data have been added for each channel at 71, 130, and 360 nM [Ca2+] so as to get a better sense of the shapes of the resulting {Delta}G(0)o-c vs. [Ca2+] relations. Upon examining these relations, we see that the zFV1/2 values, or equivalently the {Delta}G(0)o-c values, of the three channels are very similar at 3 nM [Ca2+], indicating that the intrinsic free energy difference between open and closed is unaltered by the mutations. Also, 100 µM [Ca2+] causes a smaller change in the {Delta}G(0)o-c values of the two mutants than it does that of the wild-type channel, reflecting a weaker influence of Ca2+ binding on the channel's closed-to-open conformational change. More interesting, however, are the shapes of the three {Delta}G(0)o-c vs. [Ca2+] curves. Over the [Ca2+] range 3 nM to 0.8 µM, where V1/2 for all three channels is changing appreciably (see Fig. 6 D), {Delta}G(0)o-c for all three channels remains nearly constant (Fig. 6 F). This is unexpected, because, according to the current view, it is only through altering {Delta}G(0)o-c that Ca2+ binding alters V1/2. That is, Ca2+ binding is thought to shift the mSlo G-V relation, not because it alters the channel's voltage-sensing mechanism, but because it decreases the energy barrier this mechanism must work against. Thus, V1/2 is not expected to change without a concomitant change in {Delta}G(0)o-c, and the incongruity we observe between Fig. 6, D and F, suggests that either this view is not correct, or the method we have used to estimate {Delta}G(0)o-c is inaccurate, at least at [Ca2+] less than ~1 µM.

To test the latter possibility, we used a model of BKCa channel gating developed by Horrigan et al. (1999)Go that is more sophisticated in its voltage-sensing mechanism than the VD-MWC model. It supposes four voltage sensors, one in each subunit, that move independently to promote channel opening. Voltage-sensor movement is not required for channel opening, however, but rather promotes opening in a manner similar to Ca2+ binding in an MWC system. It also supposes a weak voltage dependence associated with its central closed-to-open conformational change, and not including Ca2+ binding, it consists of the 10 states illustrated in Fig. 7 A. Adding then four identical Ca2+-binding sites that also allosterically promote opening leads to the 50-state model depicted in Fig. 7 B (Cox and Aldrich, 2000Go; Rothberg and Magleby, 2000Go). Although apparently complex, this model contains only seven free parameters, and yet it has been shown to account for mSlo macroscopic currents over wide ranges of [Ca2+] and voltage both with and without the BKCa ß1 subunit (Cox and Aldrich, 2000Go) and, as well, gating currents at nominally 0 [Ca2+] (Horrigan and Aldrich, 1999Go; Horrigan et al., 1999Go). Furthermore, a scheme of the same form was developed simultaneously to account for mSlo's single-channel gating behavior (Rothberg and Magleby, 2000Go). Thus, the model illustrated in Fig. 7 B, although overly simple in its Ca2+-sensing mechanism, represents the latest view of the relationship between Ca2+ binding and voltage sensing in BKCa-channel activation.



View larger version (33K):
[in this window]
[in a new window]
 
FIGURE 7. A more complex model of mSlo gating suggests errors in our {Delta}G(0)o-c estimates at low [Ca2+]. (A). Diagram of the model Horrigan et al. (1999)Go used to described the voltage-dependent gating of the mSlo channel. Horizontal transitions represent voltage sensor activation. Vertical transitions represent the conformational change by which the channel opens. JC(V) and JO(V) are the equilibrium constants for voltage-sensor activation when the channel is either closed or open where: JC(V) = JC(0)exp(zFV/RT) and JO(V) = JO(0)exp(zFV/RT); and L(V) is the equilibrium constant between open and closed when no voltage sensors are active and no Ca2+ is bound where: L(V) = L(0)exp(qFV/RT). (B) 50-state mSlo model. Here Ca2+ binding occurs along the long horizontal axis, voltage sensor movement along the short horizontal axis, and transitions from the upper to the lower tier indicate channel opening. The front face of this scheme corresponds to the model in A. (C) Simulated G-V curves from the model in B at the [Ca2+] indicated on the figure. For more discussion of this model see (Cox and Aldrich, 2000Go; Cui and Aldrich, 2000Go; Rothberg and Magleby, 2000Go). Parameters were as follows: JC(0) = 0.066; JO(0) = 1.10; z = 0.51; L(0) = 2e-6; q = 0.4; KC = 10 µM, KO = 1 µM. Each curve is fitted with a Boltzmann function (dashed lines). (D) Plots of {Delta}G(0)o-c as a function of [Ca2+] for the model in C. Darkened circles indicate the true {Delta}G(0)o-c of the model calculated as 4RTln[(1+ [Ca]/KC)/(1 + [Ca]/KO)] + 4RTln[(1 + JC(0))/(1 + JO(0))] - RTln[L(0)], while the open circles represent estimates obtained from Boltzmann-fit parameters as zFV1/2.

 
Fig. 7 C shows simulated G-V curves from this model. The parameters used to generate these curves were similar to those used recently to fit wild-type mSlo data (Cox and Aldrich, 2000Go). Adjacent to these curves (Fig. 7 D) the true {Delta}G(0)o-c of the model is plotted as a function of [Ca2+] (closed circles). These values were determined mathematically from the probability of open-conformation occupancy at 0 mV. Also plotted (open circles) are estimates of {Delta}G(0)o-c obtained by fitting each simulated curve with a Boltzmann function and calculating {Delta}G(0)o-c from the fit parameters.

Interestingly, although the 50-state model assumes no direct effect of Ca2+ binding on voltage sensing, the Boltzmann method of estimating {Delta}G(0)o-c severely underestimates the model's true {Delta}G(0)o-c at [Ca2+] below ~1 µM, where its G-V curves become both increasingly shallower and poorly fitted by a Boltzmann function, characteristics also observed in the data (see Figs. 2 4). Furthermore, the underestimation occurs in such a way that the estimated {Delta}G(0)o-c vs. [Ca2+] relation is almost constant between 0 and 1 µM [Ca2+], when the true relation is changing by ~23%. Thus, the unexpected, flat characteristic we observed in our analysis of the data in Fig. 6 is also observed here, and it thus seems likely that in Fig. 6 F we are indeed underestimating {Delta}G(0)o-c at low [Ca2+] when using zFV1/2 as an estimator.

To attempt to better estimate {Delta}G(0)o-c, we took advantage of the understanding of the voltage-dependent mechanism of mSlo gating developed by Horrigan et al. (1999)Go. Instead of using a Boltzmann function, we fit each G-V curve with the following function, which describes the G-V relation of the scheme diagramed in Fig. 7 A (Horrigan and Aldrich, 1999Go; Horrigan et al., 1999Go).


(2)

Here, JO(0) and JC(0) represent the equilibrium constants for voltage-sensor movement in each subunit at 0 mV when the channel is open or closed, respectively, z represents the gating charge associated with each voltage sensor, q represents the gating charge associated with the model's central conformational change, and ACa is a Ca2+-dependent factor that is logarithmically related to {Delta}G(0)o-c. That is, assuming only that Ca2+ binding promotes opening, and the voltage-dependent gating mechanism described by the scheme in Fig. 7 A, for any [Ca2+], we may write:

(3)
where the first term in Eq. 3 does not vary with [Ca2+] and the second one does.

Thus, we fit the wild-type and mutant G-V curves with Eq. 2, fixing the voltage-dependent gating parameters to values we have used previously to describe mSlo gating (Cox and Aldrich, 2000Go), and from the ACa parameter and Eq. 3 calculated {Delta}G(0)o-c. This method of estimating {Delta}G(0)o-c is precisely accurate for the 50-state model, and precisely accurate as well for any model with a voltage-sensing mechanism described by the scheme in Fig. 7 A, regardless of the number and type of Ca2+-binding sites it contains, as long as those sites do not directly interact with the channel's voltage sensors.

Estimates of {Delta}G(0)o-c as a function of [Ca2+] for the wild-type channel and three of the mutants described above are plotted in Fig. 8 A. Supporting the idea that we have improved our estimates of {Delta}G(0)o-c at low [Ca2+], the {Delta}G(0)o-c vs. [Ca2+] relations for all four channels are no longer flat in the 3 nM to 0.8 µM [Ca2+] range, but now better follow changes in V1/2. In Fig. 8 B we have shifted each curve in A vertically so as to have a value of zero at 3 nM [Ca2+], and thereby created plots of {Delta}{Delta}G(0)o-c, the change in {Delta}G(0)o-c as a function of [Ca2+], a measure of the influence of [Ca2+] on the energetics of channel opening. And in Fig. 8 C we have plotted difference curves obtained by subtracting each mutant's curve in Fig. 8 B from the wild-type curve. These curves measure the Ca2+ influence lost by each mutation and thus are reflective of the properties of the binding site each mutation disrupts.



View larger version (14K):
[in this window]
[in a new window]
 
FIGURE 8. Energetics of Ca2+ binding to mutant and wild-type mSlo channels. (A) Estimates of {Delta}G(0)o-c as a function of [Ca2+] obtained by fitting each experiment with Eq. 2 and then calculating {Delta}G(0)o-c from Eq. 3. Constant parameters were: JC(0) = 0.059, JO(0) = 1.020, z = 0.51, q = 0.4. Symbols indicate channel type as indicated on the plot. (B) The curves in A have been shifted so each curve has a value of 0 at 3 nM [Ca2+]. (C) WT-mutant difference curves. Each mutant's curve in B has been subtracted from the wild-type curve in B to create difference curves that indicate the properties of the site disrupted by each mutation. Error bars, which represent standard error of the mean, are included in A and B; they are often smaller than the plot symbols.

 
To estimate the Ca2+-binding parameters of each site we then fit each curve in Fig. 8 C with the following function, which describes the allosteric effect of four independent ligands influencing a central conformational change (Cui and Aldrich, 2000Go).


(4)

This relation has two free parameters, KC and KO, and listed in Table II are estimates of these parameters for the sites disrupted by each mutation. Examining the results from the {Delta}896–903 and M513I mutations, we see that both types of high-affinity Ca2+-binding site appear to function similarly, each having an affinity for Ca2+ of 3.5–3.8 µM when the channel is closed and 0.80–0.94 µM when the channel is open, and each maximally altering {Delta}G(0)o-c by approximately -3.4 kJ/mol as Ca2+ binds, for a total of approximately -13.5 kJ/mol for each type of site. Also listed under site 1 in Table II, however, are affinity-constant estimates derived from the WT {Delta}899–903 difference curve in Fig. 8 C, because, as discussed below, it is not clear which mutant, {Delta}896–903 or {Delta}899–903, represents elimination of just Ca2+ bowl function. The values derived from the WT-{Delta}899–903 difference curve suggest a somewhat higher Ca2+ affinity for site 1, both in the closed and open conformation (KC = 1.8 µM, KO = 0.6 µM), with a maximum {Delta}{Delta}GCa(0)o-c per site of -2.4 kJ/mol for a total of -9.5 kJ/mol.


View this table:
[in this window]
[in a new window]
 
TABLE II Estimates of the Binding Characteristics of mSlo's High-Affinity Ca2+-binding Sites

 
Are the Mutants Quantitatively Additive?
The data in Fig. 8, A and B, also allow us to determine if mutations in the Ca2+-bowl and M513I act truly independently, as is expected if each affects a distinct binding site, and binding at one site does not influence binding at the other. To test this we plotted, in Fig. 9 A, the effects of 100 µM [Ca2+] on {Delta}G(0)o-c for the mutant and wild-type channels, and, in Fig. 9 B, the degree to which each mutation inhibits the wild-type channel's response to 100 µM [Ca2+]. Consistent with Fig. 3, the {Delta}899–903 mutant is somewhat more responsive to [Ca2+] than {Delta}896–903 (Fig. 9 A), and thus it reduces the channel's response to 100 µM [Ca2+] ({Delta}{Delta}G(0)o-c; 3 nM to 100 µM) somewhat less (Fig. 9 B). Perhaps surprisingly, however, the double mutants involving M513I and either {Delta}899–903 or {Delta}896–903 behave almost identically. Thus, the extra effect of the larger deletion is not evident when combined with M513I, suggesting one pair of mutations is not completely additive. Indeed, the predictions of strict additivity are indicated on Fig. 9 B (solid black lines), and consistent with this observation, the effects of M513I and {Delta}899–903 are truly additive, whereas those of M513I and {Delta}896–903 are less than quantitatively additive by 20%.



View larger version (19K):
[in this window]
[in a new window]
 
FIGURE 9. Additivity of mutations. (A) -{Delta}{Delta}G(0)o-c (3 nM to 100 µM [Ca2+]) is plotted for each channel type. (B) For each channel type is plotted ({Delta}{Delta}G(0)o-c (3 nM to 100 µM [Ca2+])wild-type - {Delta}{Delta}G(0)o-c (3 nM to 100 µM [Ca2+])mutant)/{Delta}{Delta}G(0)o-c (3 nM to 100 µM [Ca2+])wild-type,, which is the percent reduction in response to 100 µM [Ca2+] caused by each mutation. Vertical black lines indicate the predictions of strict additivity, and are thus the sums of the bars for each individual mutant that contributes to a double mutant.

 
The meaning of this result is open to a number of interpretations. One possibility is that the two binding sites do not behave independently, in which case mutations that knockout each site separately would not be expected to be quantitatively additive when combined. Alternatively, it could be that the M513I and {Delta}896–903 mutations have partially overlapping spheres of influence. This would also be expected to produce a less than additive effect. Thus, we cannot be certain which mutation {Delta}896–903 or {Delta}899–903 best represents loss of Ca2+-bowl function, and therefore in Fig. 8 C we have analyzed both.


   DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We have described mutations at two positions in the mSlo channel's intracellular COOH-terminal domain, the Ca2+ bowl and M513, that when mutated individually decrease the channels response to Ca2+ by approximately half, as judged by G-V shifts in response to 100 µM [Ca2+]. Furthermore, the effects of mutations at each site are largely Ca2+ dependent, suggesting they are primarily altering Ca2+ binding, rather than other aspects of gating such as voltage sensing or the intrinsic energetics of channel opening. Reinforcing this conclusion, we have analyzed the behavior of mutants at both sites in terms of allosteric models that provide a means of estimating the degree to which a given mutation alters the power of Ca2+ to influence the energetics of channel opening ({Delta}{Delta}G(0)o-c), a property that depends solely on the Ca2+-binding characteristics of the channel. Results of such analyses suggest that mutations at each site decrease the energy that 100 µM [Ca2+] imparts to the channel's closed-to-open conformational change also by approximately half, as if each site were responsible for approximately half of the channel's Ca2+ sensitivity.

Are There Two High-Affinity Sites?
It is natural to question whether mutations that create similar phenotypes alter the same or different processes. One could imagine that although separated by 370 amino acids in the channel's primary sequence, the Ca2+ bowl and M513 are oriented such that mutations at both sites alter the same process by which Ca2+ binding promotes opening. The observation that the combined effect of mutations at both sites is considerably stronger than mutations at either site alone, however, argues against this hypothesis. In fact, the double mutants, {Delta}896–903+M513I and {Delta}899–903+M513I, do not respond appreciably to [Ca2+] below ~10 µM, whereas the single mutants, {Delta}896–903, {Delta}899–903, and M513I, each start to respond to [Ca2+] (as does the wild-type channel) at ~0.1 µM. Thus, mutations at the two sites are qualitatively additive, and likely alter separate binding sites. Indeed, we find that the energetic effects of {Delta}899–903 and M513I, as shown in Fig. 9 B, are very close to quantitatively additive. The difference between the measured {Delta}{Delta}Go-c value (3 nM to 100 µM) of the double mutant M513I+{Delta}899–903 and the value predicted by strict additivity is 2.25 kJ/mol, a value that is <10% of the {Delta}{Delta}Go-c (3 nM to 100 µM) of the wild-type channel, and less than is typically taken to indicate nonindependence in studies of interactions between mutations (1.5 kT or 3.69 kJ/mol; Schreiber and Fersht, 1995Go; Naranjo and Miller, 1996Go; Ranganathan et al., 1996Go). Thus, it would appear that the {Delta}899–903 and the M513I mutations are acting independently and on separate Ca2+-binding sites.

We cannot rigorously exclude the possibility, however, that the {Delta}899–903 and M513I mutations affect a single binding site in such a way that they happen to behave additively. In fact, the combined effects of the {Delta}896–903 and M513I mutations are less than additive by 5.47 kJ/mol, or 20% of the of the wild-type channel's {Delta}{Delta}Go-C value (3 nM to 100 µM). Here, one mutation partially obscures the effects of the other, suggesting a common process is altered. This result could also be explained, however, by supposing that the binding of Ca2+ at one site influences binding at the other, or more simply that in addition to disrupting the Ca2+ bowl, the {Delta}896–903 mutation has a partial effect at the second site, a conclusion that also requires the presence of two types of high-affinity sites. Thus, we think two distinct high-affinity Ca2+-binding sites are likely.

Supporting this conclusion, Bian et al. (2001)Go have found that the Ca2+-bowl mutation 897–901N reduces the dSlo channel's Hill coefficient for Ca2+ activation by a factor of ~2 at many voltages, a result consistent with the idea that the Ca2+-bowl mutation disrupts one type of high-affinity Ca2+-binding site while leaving another intact. This same mutation also reduces Ca2+ binding to a fusion protein consisting of the last 280 amino acids of dSlo by approximately half. Furthermore, Schreiber and Salkoff (1997)Go have argued for a second type of mSlo Ca2+-binding site based on the observation that mutations at the Ca2+ bowl alter the effects of Ca2+, but leave the activating effects of Cd2+ unaltered. And in another series of experiments they divided the mSlo tail domain into four regions (A–D) and found that both region B, which contains the Ca2+ bowl, and region C which does not, can each confer partial Ca2+ sensitivity upon an insensitive mSlocore–mSlo3tail chimeric channel. Thus, the two high-affinity site interpretation we have given our results is not at odds with previous work.

Do the Mutations Create Complete Loss of Function?
Does each mutation completely eliminate the effect of Ca2+ binding at the site it disrupts? The answer to this question appears to be yes, at least for the {Delta}896–903 and the M513I mutations, as when both mutations are made together all high-affinity Ca2+ response is lost, leaving only a low-affinity response blocked by high concentrations of Mg2+. Indeed, Shi and Cui (2001)Go and Zhang et al. (2001)Go have extensively examined the ability of Mg2+ to activate the mSlo channel. Both groups concluded that Mg2+ activates the channel via a site separate from the channel's high-affinity Ca2+-binding sites. This site was estimated to have an affinity for Ca2+ of 0.6–1 mM when the channel is open and 2–3 mM when the channel is closed, and an affinity for Mg2+ of 2–6 mM when the channel is open and 10–20 mM when the channel is closed (Zhang et al., 2001Go). For a channel with four such sites, a simple VD-MWC type model predicts that this low-affinity site would start to influence the position of the mSlo G-V relation at just above 10 µM [Ca2+], just as we observed with our double mutants (simulations not shown). And it would be ~90% occupied by 100 mM [Mg2+] and therefore almost fully activated. Thus, our double mutant {Delta}896–903+M513I acts exactly as one would predict if it creates a channel whose high-affinity Ca2+-binding sites have lost their function, but whose low-affinity Mg2+-sensitive, Ca2+-binding sites are still intact (Figs. 4 and 5). These results therefore lend further credence to the notion that the BKCa channel has low-affinity, Mg2+-sensitive Ca2+-binding sites that are separate from its high-affinity sites, and they pave the way for the study of these low-affinity sites in isolation.

No Other Binding Sites
The results of our double mutant experiments also argue that the mSlo channel contains no other high-affinity Ca2+-binding sites besides those disrupted by mutations at the Ca2+ bowl and M513 and no low-affinity sites that are not Mg2+sensitive. If either of these suppositions were not correct, then the {Delta}896–903+M513I channel would be expected to respond to increases in [Ca2+] in the presence high concentrations of Mg2+ via these as yet unaccounted for sites. Generally, no such response is observed. It should be mentioned, however, that although we saw no significant change in the G-V position of the {Delta}896–903+M513I channel as [Ca2+] was increased from <3 nM to 100 µM in the presence of 100 mM [Mg2+], we did see a small but consistent change in the shape of the channel's G-V relation. The relation became more shallow by 27% (Fig. 5) as [Ca2+] was raised. Whether this is due to weaker binding of Mg2+ at low voltages, or some real remaining Ca2+ effect will require further investigation.

Estimates of Ca2+-binding Characteristics
Identifying mutations that cause a complete loss of function at mSlo's high-affinity Ca2+-binding sites affords a unique opportunity to estimate the binding characteristics of these sites. Results of our analysis suggest that the site disrupted by M513I, which we refer as site 2, accounts for ~50% of the -26.7 kJ/mol that 100 µM [Ca2+] imparts to the channel closed-to-open conformational change, whereas the Ca2+ bowl site (site 1) imparts 35–50%, depending on which mutation ({Delta}899–90 or {Delta}896–903) is analyzed. Thus, type 1 and 2 sites appear to contribute roughly equally to the high-affinity effects of Ca2+. Assuming one binding site of each type per subunit (four per channel), and taking our estimates from the {Delta}896–903 mutant, we estimate KO for both sites to be 0.8–0.9 µM and KC to be 3.5–3.8 µM. Thus, both sites appear to have very similar binding characteristics. If, however, the channel has an unequal number of sites of each type, as might arise if one type of site is formed at the interface between subunits while the other is not, this may not be the case. But our estimates of the energetic contributions of each type of site (-9.5 to -13.6 kJ/mol for site 1 and -13.4 kJ/mol for site 2) would remain valid.

How Do the Mutations Create a Loss of Function?
How do the mutations we have described exert their effects? Must they act directly at the binding site? Indeed, one might suppose that mutations at either position affect Ca2+ sensitivity, not because they affect Ca2+ binding, but because they affect the ability of Ca2+ binding to influence opening. Suppose, as suggested by Schreiber et al. (1999)Go, that there is a Ca2+-binding domain that swings as if on a hinge after Ca2+ binds, and in so doing promotes opening. One might then suppose that a given mutation breaks the hinge without affecting Ca2+ binding. Thermodynamic constraints dictate, however, that in order for Ca2+ binding to promote opening, the affinity of the binding site must increase as the channel opens, and conversely, if a binding site's affinity increases as the channel opens, binding must promote opening. Thus, for a mutation to eliminate the effectiveness of a binding site, it must either prevent binding altogether, or prevent the change in affinity that occurs at the binding site as the channel opens. Thus, binding and the transduction of binding to opening are inextricably linked, and the mutations we have described must therefore affect Ca2+ binding.

One might still ask, however, whether the mutations we have studied disrupt residues that actually make contact with Ca2+? Our data do not rule out the possibility that some protein endogenous to the oocyte associates with mSlo to form its Ca2+ sensors, as is the case for the small-conductance Ca2+-activated K+ channel, which uses constitutively associated calmodulin as its Ca2+ sensor (Xia et al., 1998Go; Wissmann et al., 2002Go). It is conceivable, then, that mutations at the Ca2+ bowl or M513 disrupt such an association. For the Ca2+ bowl mutations, however, there are good reasons to think this is not the case. First, these mutations remove acidic side chains that could provide the negative charge density normally associated with Ca2+ coordination. Second, two laboratories have found that Ca2+ will bind to recombinant Slo fusion proteins overexpressed and purified from bacteria, and will do so in gel overlay assays where auxiliary proteins would not likely remain associated (Bian et al., 2001Go; Braun and Sy, 2001Go). And third, as mentioned above, the Ca2+-bowl mutation 897–901N inhibits this binding by approximately half (Bian et al., 2001Go). Thus, it seems likely that the Ca2+ bowl does in fact form a Ca2+-binding site, although this is not directly demonstrated by our data.

The situation for the M513I mutation is less clear. Methionines are not commonly Ca2+ coordinators, so it may be that M513I interferes with the linkage between binding and opening (by linkage we mean that which causes the binding site to increase its affinity as the channel opens). This could be the movement of a key atom closer to the Ca2+ ion in its site, as is the case with cGMP when the cyclic nucleotide–gated channel opens (Varnum et al., 1995Go). Or perhaps the M513I mutation inhibits an intersubunit association of intracellular domains that may occur after Ca2+ binds. Such an association has been shown to occur after Ca2+ binds to the SK channel (Schumacher et al., 2001Go), and M513 lies at the very end of mSlo's seventh hydrophobic domain (S7), a region that might reasonably bury its hydrophobic side chains by dimerization, perhaps with its counterpart on a neighboring subunit, when prodded by Ca2+ binding.

Pursuing this idea further, M513 also lies at the very end of mSlo's RCK domain, which spans from just after S6 to the end of S7 (Jiang et al., 2001Go). In the recent crystal structure of the E. coli K+ channel's RCK domain (Jiang et al., 2001Go), the {alpha}G helix, where M513 would be predicted to lie, participates in a dimer interface. Thus, although the homology in this region between the E. coli K+ channel and mSlo is low, and the E. coli channel does not contain a methionine at the position corresponding to mSlo's 513, by structural analogy we might still speculate that M513 is important for the formation of dimers between mSlo's intracellular domains.

Alternatively, it could also be that M513I disrupts a Ca2+-binding site formed by the RCK domain itself, since much of this domain has a structure commonly associated with ligand-binding proteins known as a Rossman fold (Branden and Tooze, 1991Go) and Rossman folds have been shown to form binding sites for a variety of small molecules, including ions (Lee et al., 1995Go). So, although M513 is not in mSlo's predicted Rossman fold, its close proximity suggests its mutation may destabilize a Ca2+-binding site formed therein (Jiang et al., 2001Go).

It is interesting to note that the mSlo3 channel, a Ca2+-insensitive mSlo homologue, contains neither a Ca2+ bowl nor a methionine at position 513 (Schreiber et al., 1998Go) and thus its lack of Ca2+ sensitivity is easily explained in terms of our results. Less easy to explain, however, is the observation that a channel containing the mSlo core and the mSlo3 tail, and thus one that has lost its Ca2+ bowl but retains M513, is in large-part Ca2+ insensitive (Schreiber et al., 1999Go; Moss and Magleby, 2001Go). This result could mean that there are elements of the mSlo tail besides those we have mutated that are required for the proper functioning of a binding site in the core of the channel, or conversely, that M513 interacts with a Ca2+-binding site in the mSlo tail, besides the Ca2+ bowl, that is not present in the mSlo3 tail. Clearly, confirmation of any of these hypotheses must await new physiological and structural data.


   FOOTNOTES
 
Abbreviations used in this paper: BKCa, large-conductance Ca2+-activated K+ channel; KV, voltage-gated K+ channel; VD-MWC, voltage-dependent Monod-Wyman-Changeux.


   ACKNOWLEDGMENTS
 
We gratefully acknowledge Dr. Kathleen Dunlap, Dr. Robert Blaustein, and Christina Kaldany for helpful comments on the manuscript.

This work was supported by grant R01HL64831 from the National Institutes of Health and by a grant from The Jessie B. Cox Charitable Trust and The Medical Foundation.

Submitted: 14 May 2002
Revised: 11 June 2002
Accepted: 11 June 2002


   REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 


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

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

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

    Bian, S., I. Favre, and E. Moczydlowski. 2001. Ca2+-binding activity of a COOH-terminal fragment of the Drosophila BK channel involved in Ca2+-dependent activation. Proc. Natl. Acad. Sci. USA. 98:4776–4781.[Abstract/Free Full Text]

    Branden, C., and J. Tooze. 1991. Enzymes that bind nucleotides. In Introduction to Protein Structure. New York, Garland Publishing. 141–159.

    Braun, A.F., and L. Sy. 2001. Contribution of potential EF hand motifs to the calcium-dependent gating of a mouse brain large conductance, calcium-sensitive K(+) channel. J. Physiol. 533:681–695.[Abstract/Free Full Text]

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

    Cox, D.H., and R.W. Aldrich. 2000. Role of the ß1 subunit in large-conductance Ca(2+)-activated K(+) channel gating energetics. Mechanisms of enhanced Ca(2+) sensitivity. J. Gen. Physiol. 116:411–432.[Abstract/Free Full Text]

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

    Cui, J., and R.W. Aldrich. 2000. Allosteric linkage between voltage and Ca(2+)-dependent activation of BK-type mslo1 K(+) channels. Biochemistry. 39:15612–15619.[CrossRef][Medline]

    Cui, J., D.H. Cox, and R.W. Aldrich. 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., M. Wallner, E. Stefani, L. Toro, and R. Latorre. 1996. Interaction of internal Ba2+ with a cloned Ca2+-dependent K+ (hslo) channel from smooth muscle. J. Gen. Physiol. 107:399–407.[Abstract]

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

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

    Goldstein, S.A., D.J. Pheasant, and C. Miller. 1994. The charybdotoxin receptor of a Shaker K+ channel: peptide and channel residues mediating molecular recognition. Neuron. 12:1377–1388.[Medline]

    Golowasch, J., A. Kirkwood, and C. Miller. 1986. Allosteric effects of Mg2+ on the gating of Ca2+-activated K+ channels from mammalian skeletal muscle. J. Exp. Biol. 124:5–13.[Abstract]

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

    Hidalgo, P., and R. MacKinnon. 1995. Revealing the architecture of a K+ channel pore through mutant cycles with a peptide inhibitor. Science. 268:307–310.[Medline]

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

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

    Jiang, Y., A. Pico, M. Cadene, B.T. Chait, and R. MacKinnon. 2001. Structure of the RCK domain from the E. coli K+ channel and demonstration of its presence in the human BK channel. Neuron. 29:593–601.[CrossRef][Medline]

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

    Lee, J.O., P. Rieu, M.A. Arnaout, and R. Liddington. 1995. Crystal structure of the A domain from the alpha subunit of integrin CR3 (CD11b/CD18). Cell. 80:631–638.[Medline]

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

    Meera, P., M. Wallner, Z. Jiang, and L. Toro. 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.[CrossRef][Medline]

    Methfessel, C., and G. Boheim. 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 R. Latorre. 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]

    Moss, B.L., and K.L. Magleby. 2001. Gating and conductance properties of BK channels are modulated by the S9-S10 tail domain of the alpha subunit. A study of mSlo1 and mSlo3 wild-type and chimeric channels. J. Gen. Physiol. 118:711–734.[Abstract/Free Full Text]

    Naranjo, D., and C. Miller. 1996. A strongly interacting pair of residues on the contact surface of charybdotoxin and a Shaker K+ channel. Neuron. 16:123–130.[Medline]

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

    Ranganathan, R., J.H. Lewis, and R. MacKinnon. 1996. Spatial localization of the K+ channel selectivity filter by mutant cycle-based structure analysis. Neuron. 16:131–139.[Medline]

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

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

    Rothberg, B.S., and K.L. Magleby. 1998. Investigating single-channel gating mechanisms through analysis of two-dimensional dwell-time distributions. Methods Enzymol. 293:437–456.[Medline]

    Rothberg, B.S., and K.L. Magleby. 2000. Voltage and Ca2+ activation of single large-conductance Ca2+-activated K+ channels described by a two-tiered allosteric gating mechanism. J. Gen. Physiol. 116:75–99.[Abstract/Free Full Text]

    Schreiber, G., and A.R. Fersht. 1995. Energetics of protein-protein interactions: analysis of the barnase-barstar interface by single mutations and double mutant cycles. J. Mol. Biol. 248:478–486.[CrossRef][Medline]

    Schreiber, M., and L. Salkoff. 1997. A novel calcium-sensing domain in the BK channel. Biophys. J. 73:1355–1363.[Abstract]

    Schreiber, M., A. Wei, A. Yuan, J. Gaut, M. Saito, and L. Salkoff. 1998. Slo3, a novel pH-sensitive K+ channel from mammalian spermatocytes. J. Biol. Chem. 273:3509–3516.[Abstract/Free Full Text]

    Schreiber, M., A. Yuan, and L. Salkoff. 1999. Transplantable sites confer calcium sensitivity to BK channels. Nat. Neurosci. 2:416–421.[CrossRef][Medline]

    Schumacher, M.A., A.F. Rivard, H.P. Bachinger, and J.P. Adelman. 2001. Structure of the gating domain of a Ca2+-activated K+ channel complexed with Ca2+/calmodulin. Nature. 410:1120–1124.[CrossRef][Medline]

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

    Shi, J., and J. Cui. 2001. Intracellular Mg(2+) enhances the function of BK-type Ca(2+)-activated K(+) channels. J. Gen. Physiol. 118:589–606.[Abstract/Free Full Text]

    Stefani, E., M. Ottolia, F. Noceti, R. Olcese, M. Wallner, R. Latorre, and L. Toro. 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]

    Stocker, M., and C. Miller. 1994. Electrostatic distance geometry in a K+ channel vestibule. Proc. Natl. Acad. Sci. USA. 91:9509–9513.[Abstract/Free Full Text]

    Talukder, G., and R.W. Aldrich. 2000. Complex voltage-dependent behavior of single unliganded calcium-sensitive potassium channels. Biophys. J. 78:761–772.[Abstract/Free Full Text]

    Varnum, M.D., K.D. Black, and W.N. Zagotta. 1995. Molecular mechanism for ligand discrimination of cyclic nucleotide-gated channels. Neuron. 15:619–625.[Medline]

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

    Wissmann, R., W. Bildl, H. Neumann, A.F. Rivard, N. Klocker, D. Weitz, U. Schulte, J.P. Adelman, D. Bentrop, and B. Fakler. 2002. A helical region in the C terminus of small-conductance Ca2+-activated K+ channels controls assembly with apo-calmodulin. J. Biol. Chem. 277:4558–4564.[Abstract/Free Full Text]

    Xia, X.M., B. Fakler, A. Rivard, G. Wayman, T. Johnson-Pais, J.E. Keen, T. Ishii, B. Hirschberg, C.T. Bond, S. Lutsenko, et al. 1998. Mechanism of calcium gating in small-conductance calcium-activated potassium channels. Nature. 395:503–507.[CrossRef][Medline]

    Zhang, X., C.R. Solaro, and C.J. Lingle. 2001. Allosteric regulation of BK channel gating by Ca(2+) and Mg(2+) through a nonselective, low affinity divalent cation site. J. Gen. Physiol. 118:607–636.[Abstract/Free Full Text]