In this work we address the question of the KCa3.1 channel pore structure in the closed configuration in relation to the contribution of the C-terminal end of the S6 segments to the Ca2+-dependent gating process. Our results based on SCAM (substituted cysteine accessibility method) experiments first demonstrate that the S6 transmembrane segment of the open KCa3.1 channel contains two distinct functional domains delimited by V282 with MTSEA and MTSET binding leading to a total channel inhibition at positions V275, T278, and V282 and to a steep channel activation at positions A283 and A286. The rates of modification by MTSEA (diameter 4.6 Å) of the 275C (central cavity) and 286C residues (S6 C-terminal end) for the closed channel configuration were found to differ by less than sevenfold, whereas experiments performed with the larger MTSET reagent (diameter 5.8 Å) resulted in modification rates 103–104 faster for cysteines at 286 compared with 275. Consistent with these results, the modification rates of the cavity lining 275C residue by MTSEA, Et-Hg+, and Ag+ appeared poorly state dependent, whereas modification rates by MTSET were 103 faster for the open than the closed configuration. A SCAM analysis of the channel inner vestibule in the closed state revealed in addition that cysteine residues at 286 were accessible to MTS reagents as large as MTS-PtrEA, a result supported by the observation that binding of MTSET to cysteines at positions 283 or 286 could neither sterically nor electrostatically block the access of MTSEA to the closed channel cavity (275C). It follows that the closed KCa3.1 structure can hardly be accountable by an inverted teepee-like structure as described for KcsA, but is better represented by a narrow passage centered at V282 (equivalent to V474 in Shaker) connecting the channel central cavity to the cytosolic medium. This passage would not be however restrictive to the diffusion of small reagents such as MTSEA, Et-Hg+, and Ag+, arguing against the C-terminal end of S6 forming an obstructive barrier to the diffusion of K+ ions for the closed channel configuration.
Introduction
Ca2+-activated potassium channels (KCa) are present in most mammalian cell types, where their primary role is to establish a link between the various Ca2+-based second messenger systems and the electrical properties of the cells. Three main classes of KCa to date have been identified based on their permeation properties and pharmacology (Vergara et al., 1998). They include the charybdotoxin- and iberiotoxin-sensitive KCa1.1 channels of large conductance (150–220 pS), the intermediate conductance (20–50 pS) KCa3.1 channels inhibited by clotrimazole (Rittenhouse et al., 1997) and TRAM34 (Wulff et al., 2001), and the apamine-sensitive and -insensitive SK channels of small conductance (KCa2.1, KCa2.2, and KCa2.3) (Kohler et al., 1996; Stocker, 2004). The KCa3.1 channel is a tetrameric protein with each subunit comprising 427 amino acids organized in six transmembrane segments S1–S6 with a pore motif between segments 5 and 6. In contrast to KCa1.1, the gating process of SK and KCa3.1 is voltage insensitive and the Ca2+ sensitivity is conferred by the Ca2+-binding protein calmodulin (CaM), constitutively bound in the C terminus to each of the channel subunits in a 1:1 ratio (Khanna et al., 1999). CaM is also essential for the assembly and trafficking of the SK and KCa3.1 channel subunits (Joiner et al., 2001; Lee et al., 2003).
A 3D homology-based model of the pore-forming S6 transmembrane segment for the closed KCa3.1 configuration was proposed by our laboratory (Simoes et al., 2002) using the bacterial KcsA channel structure as template (Doyle et al., 1998). The resulting radial distribution of the α carbons for residues V275 to N292 along the S6 transmembrane segment is illustrated in Fig. 1 A. As seen, the V275, T278, and V282 residues are presented as lining the channel pore with V275 and T278 contributing to the formation of a central inner cavity ∼10 Å wide. The V284 and V285 residues are predicted in turn to be oriented opposite to the pore lumen with the residue A286 at the C-terminal end of S6 pointing toward the pore central axis. More importantly, the diameter of the KCa3.1 conducting pathway is expected to vary along the channel central axis of diffusion with a minimum van der Waals diameter of 2.0 Å at the level of the V282 residue. It follows therefore that a pore structure for the closed KCa3.1 channel based on a KcsA template would be characterized by a bundle crossing region extending from V282 to A286 with the presence of a tight hydrophobic seal at the level of the V282 residue. Data supporting this model would strongly argue for a KCa3.1 activation gate located at the C-terminal end of the transmembrane S6 segments (for example see LeMasurier et al., 2001; Cordero-Morales et al., 2006).
In this work we address the question of the KCa3.1 channel pore structure in the closed configuration. Our results provide evidence that the pore structure of the closed KCa3.1 channel cannot be accounted for by the inverted teepee-like structure prevailing for KcsA, but rather support a model where the closed KCa3.1 is characterized by a narrow passage centered at V282 connecting the channel central cavity to the cytosolic medium. We conclude that the Ca2+-dependent gating of KCa3.1 involves different sections of the S6 segment with the C-terminal end of S6 not constituting a hydrophobic seal capable to control K+ ion flow.
Materials And Methods
Cloning, Sequencing, and Site-directed Mutagenesis of the KCa3.1 Channel
KCa3.1 channel cDNA was obtained by RT-PCR from HeLa cells as previously described (Simoes et al., 2002). Site-directed mutagenesis was performed using the QuickChange site-directed mutagenesis kit (Stratagene). Point mutations were obtained using 25-mer mutated oligonucleotides with the wild-type KCa3.1 as template and confirmed by sequencing the entire coding region on both strands by automated fluorescently labeled sequencing (Bio S&T Inc.).
Oocytes
Mature oocytes (stage V or VI) were obtained from Xenopus laevis frogs anaesthetized with 3-aminobenzoic acid ethyl ester. The follicular layer was removed by incubating the oocytes in a Ca2+-free Barth's solution containing collagenase (1.6 mg/ml; Sigma-Aldrich) for 45 min. Defolliculated oocytes were stored at 18°C in Barth's solution supplemented with 5% horse serum, 2.5 mM Na-pyruvate, 100 U/ml penicillin, 0.1 mg/ml kanamycin, and 0.1 mg/ml streptomycin. The composition of the Barth's solution was (in mM) 88 NaCl, 3 KCl, 0.82 MgSO4, 0.41 CaCl2, 0.33 Ca(NO3)2, and 5 HEPES (pH 7.6). Oocytes were patched 3–5 d after coinjection of 0.1–1 ng of the cDNA coding for KCa3.1 and 1 ng of cDNA coding for a green fluorescent protein that was used as a marker for nuclear injection.
Prior to patch clamping, defolliculated oocytes were bathed in a hyperosmotic solution containing (in mM) 250 KCl, 1 MgSO4, 1 EGTA, 50 sucrose, and 10 HEPES buffered at pH 7.4 with KOH. The vitelline membrane was then peeled off using fine forceps, and the oocyte was transferred to a superfusion chamber for patch clamp measurements.
Solutions
The bath and patch pipette solutions contained (in mM) 200 K2SO4, 1.8 MgCl2, 0.025 CaCl2, 25 HEPES, buffered at pH 7.4 with KOH. The use of sulfate salts prevented excessive contaminations from endogenous Ca2+-dependent chloride channels while enabling to chelate contaminant divalent cations such as Ba2+ (maximum free Ba2+ concentration: 0.5 nM in 200 mM K2SO4). Calcium-free solutions were prepared by adding 1 mM EGTA to 200 mM K2SO4 solutions without CaCl2. In experiments where Et-Hg+ was used as thiol modifying agent, the bath and patch pipette solution consisted of (in mM) 145 K-gluconate, 5 KCl, 2.5 MgCl2, 0.1 EGTA, and 10 HEPES (pH adjusted to 7.4 with KOH). Ag+ solutions were prepared by adding Ag2SO4 to Cl−-free 200 mM K2SO4 solutions that contained 60 mM EDTA without CaCl2 and where MgCl2 has been substituted by MgSO4. The free Ag+ concentration was calculated using the Eqcal software (Biosoft). Ca2+-free 200 K2SO4 solutions at pH 6.5 and 5.5 were prepared by replacing HEPES by 25 mM MES plus KOH with either 10 mM EGTA (pH 5.5) or 2 mM EGTA (pH 6.5). For Ca2+-free 200 K2SO4 solutions at pH 8.5, pH was adjusted with (in mM) 25 Tris plus HCl in the presence of 1 EGTA. 2-Aminoethyl methanethiosulfonate hydrobromide (MTSEA), sodium (2-sulfanoethyl) methanethiosulfonate (MTSES), [2-(trimethylammonium)ethyl] methanethiosulfonate bromide (MTSET), 2-(aminocarbonyl)ethyl methanethiosulfonate (MTSACE), and [3-(triethylammonium) propyl] methanethiosulfonate (MTS-PtrEA) (Toronto Research Chemicals Inc.) were added directly into the bath solution just before application. Stock solution of MTSACE (500 mM) was prepared in DMSO. Bath solution changes were performed as described previously using a RSC-160 rapid solution changer system (BioLogic). The solution exchange time was estimated at 30 ms (Banderali et al., 2004).
Patch Clamp Recordings
Inside-out single channel recordings were performed using an Axopatch 200A amplifier (Axon Instruments, Inc.). Patch pipettes were pulled from borosilicate capillaries using a Narishige pipette puller (model PP-83) and used uncoated. The resistance of the patch electrodes ranged from 2 to 5 MΩ. Unless specified otherwise, the membrane potential is expressed as −Vp, where Vp is the pipette applied potential. Data acquisition was performed using a Digidata 1320A acquisition system (Axon Instruments, Inc.) at a sampling rate of 2.0 kHz and filtered at 500 Hz. Experiments were performed at room temperature (22°C).
Models To Estimate the Rate of Cysteine Modification by MTS Reagents
The binding of a MTS molecule to a cysteine residue at a channel site Sf can formally be expressed as (Wilson and Karlin, 1998)where n corresponds to the number of MTS-cysteine complexes formed n ∈ {0, 3}, ki, the MTS entry rate to the cysteine binding site, ko, the MTS exit rate from the cysteine binding site, km, the binding rate of the MTS to a single cysteine, and
and the binding of four MTS reagents to a channel can now be modeled as a five state kinetic scheme with
Because the proposed kinetic scheme with Kn given by Eq. 1 leads to a hardly tractable expression and to limit the number of adjustable parameters in our curve fitting procedure, we chose to analyze our experimental data according to two different kinetic models.
(1) Model I.
In conditions where the exit rate from the binding site is much faster than the reaction with the cysteine (km ≪ ko), the channel modification rate at low MTS concentrations corresponds to Kn = (4 − n)kmkifn[MTS]/ko (see Appendix) with
(2) Model II.
In conditions where the reaction rate with the cysteine is considered much faster than the entry and exit rates ki and ko (km ≫ ko + ki[MTS]), the rates of transition Kn can be expressed as
The time-dependent variation in current following the addition of MTS molecules can now be expressed as
where P(Bn,t) is the probability for a channel to be in state Bn at time t, N/(1 + Aξn) the mean number of open channels in state Bn, A the channel open/closed equilibrium constant before the addition of MTS, ΔI the unitary current amplitude, and g, the fractional decrease in unitary current following the binding of a single MTS molecule to the channel. With fn = 1, Eq. 2 reduces for the Model I to the expression proposed by Flynn and Zagotta (2003) to determine the effects of MTS on CNG channel gating and permeation. The parameter ξ was introduced to account for the change in open probability when n MTS molecules are added to the channel structure. The expression Aξn can formally be expressed as exp(−(ΔG + nΔGMTS)/KT) where ΔG is the difference in free energy between the channel closed and open configurations and ΔGMTS the change in free energy between the closed and open configurations due to the binding of a single MTS reagent. If ξ = 1, the channel open probability becomes independent of the number of MTS–cysteine complexes formed, and Eq. 2 describes the time-dependent current inhibition resulting from the binding of MTS molecules to the channel. The time course of the current activation or inhibition was analyzed by adjusting the parameters N, f0, A, ξ, and KC (Model I) or N, A, ξ,
Rates of MTS-induced modification could under these conditions be directly computed as the ratio K0/[MTS], where K0 is the rate of channel block for a given MTS concentration. Curve fitting the experimental data to Eq. 2 and theoretical predictions based on Model I and Model II were performed using Mathematica 5.1 (Wolfram Research Inc.).
MTSEA Hydrolysis
Because the rate of MTSEA hydrolysis increases at alkaline pH (Karlin and Akabas, 1998), the time course of the MTSEA-induced current variations for experiments performed at pH 8.5 was corrected to account for the time-dependent changes in MTSEA concentration. The mean lifetime of MTSEA at pH 8.5 was estimated at 92 s based on absorbance measurements at 245 nm performed in 200 mM K2SO4 conditions.
Computer-based Homology Modeling
The closed KCa3.1 channel model structure was produced by automated homology modeling using Modeller v7.0 and involved the generation of 150 models using KcsA (PDB 1K4C) as template (Sali and Blundell, 1993). The model with the lowest objective function and the lowest RMS deviation between the template and the model was kept and used as model structure for the closed KCa3.1 channel. The overall structural quality of the generated models was evaluated by PROCHECK (Laskowski et al., 1993) and ProQ (Wallner and Elofsson, 2003).
Online Supplemental Materials
The PDB file generated by Modeller containing the atom coordinates of the best model obtained for the closed KCa3.1 has been included as online supplemental material (available at http://www.jgp.org/cgi/content/full/jgp.200609726/DC1).
Results
Evidence for Two Functional Domains along S6
Fig. 2 presents typical macroscopic currents in response to channel modification by internal application of MTSEA for cysteine residues engineered at positions 275–287 along S6. MTSEA is a positively charged compound with dimensions corresponding to a cylinder 8.2 Å long with a diameter of 4.6 Å (space-filling estimation). The inside-out current recordings illustrated in Fig. 2 were performed in the presence of 25 μM Ca2+ and thus reflect the action of internal MTSEA on the open KCa3.1 channel (Po ranging from 0.2 to 0.4; see Banderali et al., 2004). The results in Fig. 2 confirm that the wild-type KCa3.1 is rather insensitive to MTSEA (1 mM), in accordance with the model proposed in Fig. 1 A where the C276 and C277 residues are oriented opposite to the pore lumen and thus poorly exposed to the water-filled conduction pathway. In contrast, a total and irreversible MTSEA-dependent current inhibition was recorded with the V275C and T278C mutants (Fig. 2) and with the V282C channel as well (not depicted), suggesting that the V275, T278, and V282 residues are accessible to water and oriented toward the channel pore. In addition, fitting the current traces obtained for the V275C, T278C, and V282C mutants to Eq. 2 revealed a time course for current inhibition characterized by a two exponential function and a minimal fractional decrease in unitary current per MTSEA of 80% (parameter g = 0.8 in Eq. 2). This analysis thus suggests that the binding of two MTSEA in the channel central cavity may be required for a total current inhibition, with the first MTSEA already accounting for 80% of the current decrease. In contrast to the total current inhibition recorded with the V275C, T278C, and V282C mutants, the modification by MTSEA of cysteine residues engineered at positions 283 and 286 and to a lesser extent 285 (unpublished data) led to a robust current increase. Ca2+ dose–response measurements performed on the A286C mutant showed no change in the Ca2+ concentration for half activation compared with wild-type KCa3.1 ([Ca2+]1/2 = 1.2 μM) (Simoes et al., 2002), but the Ca2+ sensitivity of the A286C and A283C mutants activated by MTSET appeared slightly increased with a [Ca2+]1/2 estimated at 300 nM (unpublished data). The increase in channel activity observed following the binding of MTS reagents cannot therefore be ascribed to a change in Ca2+ sensitivity, as 25 μM Ca2+ constitutes a saturating Ca2+ concentration with and without MTS binding. Examples of the single channel current fluctuations illustrating the action of MTSEA on the A286C mutant are presented in Fig. 3. As seen, the binding of MTSEA to A286C resulted in a ≈50% decrease in unitary current amplitude accompanied by a strong increase in open probability (35-fold). These observations confirm the single-channel results on the MTSET-dependent activation of the A283C and A286C mutants reported in a previous work (Simoes et al., 2002). On the basis of these parameters, the time course of the current variations initiated by MTSEA for A286C and A283C channels was curve fitted to Eq. 2 using g = 0.15 with ξ ranging from 0.2 to 0.7. The current variations were best described by the successive binding of four MTSEA molecules to the channel, leading to either multiple exponential or sigmoid type variation in current amplitude as a function of time depending on the value of the parameter ξ. Altogether, these results indicate that the transmembrane S6 segment has two functional domains delimited by V282, a first domain extending from V275 to V282 constituting the channel central cavity and a second domain extending from A283 to A286 involved in the control of the channel open probability.
The Open Pore Structure as Revealed by MTSEA and MTSET
Fig. 4 summarizes the modification rates derived from Eq. 2 for cysteine residues engineered along S6 using either MTSEA or MTSET as modifying agents. These experiments were performed in the presence of 25 μM free Ca2+ to ensure that KCa3.1 channels were maximally activated. This figure clearly shows that MTSEA did not strongly discriminate between cysteines at 275 predicted to be lining the channel central cavity (Fig. 1 A), and cysteines engineered at the C-terminal end of S6 (A286) (Fig. 4 A). Notably, the slowest rates of modification were observed with the T278C and V282C mutants, respectively, indicating that these residues were less accessible than 275C to MTSEA modification. These conclusions remain valid whether the reaction rates Kn were computed according to Model I where Kn = (4 − n)KC[MTS]fn, with KC the modification rate by cysteine or to Model II with
Accessibility of Pore Residues in the Closed Conformation
A similar approach was used to probe the pore structure of the closed KCa3.1 channel. In these inside-out experiments, pulses of MTS reagents were applied for periods ranging from 0.5 to 5 s in zero internal Ca2+ conditions. The binding of the MTS reagent to the target cysteines was assayed by measuring the magnitude of the currents in 25 μM internal Ca2+ after each MTS pulse following the washout of the thiol modifying agent (see Fig. 5). The action of MTSEA and MTSET on the closed V275C and A286C channels, respectively, is illustrated in Fig. 5. The inside-out recordings presented in this figure provide evidence that cysteine residues engineered at positions 275 and 286 are accessible to MTSEA and MTSET, respectively, in zero internal Ca2+ conditions. The modification rates measured in zero Ca2+ for cysteines substituting for S6 residues predicted on the basis of the KcsA structure to be facing the channel pore are summarized in Fig. 6 (A and B). The rates were computed according to Model I and Model II by fitting to Eq. 2 the time-dependent current variations resulting from the repetitive application of either MTSEA or MTSET. The overall profile presented in Fig. 6 A for MTSEA shows that the modification rate per cysteine KC estimated from Model I varied from 480 ± 22 M−1s−1 (N = 3) for cysteines at the C-terminal end of S6 (286C) to 29 ± 3 M−1s−1 (N = 7) and 7 ± 3 M−1s−1 (N = 2) for cysteines at positions 275 and 278, respectively. In contrast, the modification rates estimated for the larger MTSET were reduced 104-fold for cysteines at 275 compared with cysteines at the cytoplasmic end of S6 (211 ± 25 M−1s−1, N = 4, for 286C). The access of MTS reagents to residues inside the channel central cavity for the closed KCa3.1 appears therefore to strongly depend upon the dimensions of the reagents comparatively to residues located down 282. This conclusion holds independently of the kinetic model used. Experiments performed using the T278C and V282C mutants led to no detectable MTSET-induced current modifications over a period of 5 min for the closed channel (Simoes et al., 2002).
How Does Channel Opening Affect the Channel Pore Structure ?
The results presented in Figs. 4 and 6 were also used to probe the conformational changes between the KCa3.1 pore structure in zero and saturating (25 μM) internal Ca2+, respectively. Fig. 7 A summarizes the rates of modification by MTSEA and MTSET of cysteine residues substituting for amino acids predicted to be facing the pore for the KCa3.1 open and closed configurations. MTSEA was found to poorly discriminate between the channel open and closed conformations with modification rates for the V275C and A286C mutants, respectively, 7 and 1.7 times higher in 25 μM internal Ca2+ relative to zero internal Ca2+ conditions. A sevenfold difference in the V275C modification rate by MTSEA between the open and closed states reflects structural changes affecting the accessibility of residues inside the channel central cavity, but can hardly account for the gating of KCa3.1 by Ca2+. Control experiments were also performed to determine to what extent the presence of two endogenous cysteines at positions 276 and 277 could contribute to the weak state dependence of the V275C accessibility to MTSEA. Experiments performed using the triple mutant V275C-C276G-C277S showed that the open/closed channel accessibility to MTSEA differs by ≈4-fold, indicating that our conclusion on MTSEA discriminating weakly between the open and closed states truly reflects the accessibility of the cavity lining 275C residues. Fig. 7 A also indicates that the modification rates estimated for the T278C and V282C mutants is state independent, suggesting that the region distal to V282 is not obstructive to MTSEA diffusion for the open and closed configurations. In contrast, the modification rates measured with MTSET for cysteines at 275 was reduced at least 103-fold in zero Ca2+ as compared with rates measured in 25 μM Ca2+ conditions. This 103-fold variation is likely to represent a minimum value as V275C channels in zero Ca2+ were on the average inhibited by <20% after a 7-min exposure to MTSET (10 mM). The results presented in Fig. 7 A also indicate that the modification rates estimated for cysteines at positions 283 and 286 are poorly state dependent, confirming that the accessibility of these residues is not influenced by the conformational changes triggered by channel opening. Finally, the effect of the pore dimensions in the closed configuration on the accessibility to the 275C was further investigated using ethyl mercury (Et-Hg+) and Ag+ as probes. Et-Hg+ is smaller than MTSEA with dimensions corresponding to a sphere of a 4.1 Å diameter (space-filling estimation). This molecule is closer to the dimensions of a dehydrated K+ ion (2.66 Å diameter) while being more hydrophilic than MTSEA. With a diameter of 2.55 Å, Ag+ is similar in size to a K+ ion. Because experiments performed using Ag+ as a probe required a strong buffering of the Ag+ concentration with the Ca2+-chelator EDTA, accessibility measurements could not be performed in this case for the open channel configuration. Examples of inside-out current recordings illustrating the respective action of Et-Hg+ and Ag+ (25 nM) on the V275C mutant are presented in Fig. 8. Both reagents caused a strong current inhibition in zero Ca2+ conditions, indicating that they could access the channel cavity with the channel in the closed state. The results of these experiments are shown in Fig. 7 B. The time course of the current inhibition induced by Et-Hg+ and Ag+ could be well fitted by a single exponential, so that the rates of modification were computed for a simple two state model using K0 as described in Materials and methods. The modification rates measured with Et-Hg+ for the closed and open configurations of the V275C mutant were found to be nearly identical with values of 5,000 ± 1,000 s−1M−1 (n = 3) and 4,000 ± 800 s−1M−1 (n = 3) in zero and 25 μM internal Ca2+, respectively. A modification rate of 3.1 ± 0.6 × 107 s−1M−1 (n = 4) was estimated for Ag+ acting on the closed V275C channel. This value is less than 10 times smaller than the diffusion-limited rate expected for this ion (∼108 s−1M−1; del Camino and Yellen, 2001). Taken together, these results support the proposal that small positively charged thiol-reactive reagents such as MTSEA, Et-Hg+, and Ag+ permeate through the water-filled channel pore of the closed KCa3.1 channel and have access to the cavity lining residue V275 in a state-independent manner.
Do Modification Rate Measurements Truly Reflect MTS Accessibility to the Binding Site ?
The expression described in Eq. 1 indicates that the global rate of modification Kn measured experimentally includes the rate of reaction km, the exit rate ko, and the entry rate ki. It follows that the rate measured experimentally cannot always be considered as an estimation of accessibility as it is not a direct evaluation of ki, the entry rate of the MTS reagent to the binding site. Experiments were thus conducted where the accessibility parameter ki was analyzed from the modification rates measured with MTSEA and MTSET on the V275C mutant. The validity of Eq. 1 in Materials and methods was first tested by estimating the modification rate K0 as a function of the MTSEA concentration for the V275C mutant in 25 μM and zero internal Ca2+ conditions. As the current inhibition by MTSEA for the closed and open channel could be well approximated by a single exponential function, modification rates were thus computed using g = 1 in Eq. 2 with K0 given by Eq. 3.
The results of these experiments are summarized in Fig. 9 A. The modification rates K0 measured in zero internal Ca2+ show a strictly linear dependence as a function of [MTSEA] for concentrations ranging from 0.1 to 30 mM. This observation indicates that the cysteine binding site at 275 has a low affinity for MTSEA ([MTSEA]1/2 ≫ 30 mM; see Appendix). The modification rates measured in 25 μM Ca2+ could be fitted to Eq. 1, leading to km = 7.0 s−1 and [MTSEA]1/2 = 35 mM. A [MTSEA]1/2 value of ≈35 mM contrasts with an estimated [MTSEA]1/2 value of 100 μM reported for the Kir 6.2 channel, suggesting important differences in MTSEA entry and exit rates between the two channels (Phillips et al., 2003). Again, as shown in Fig. 7 A, the modification by MTSEA of the cysteines at position 275 appeared poorly state dependent over the entire MTSEA concentration range considered. Clearly, experiments performed at 1–5 mM MTSEA are within the concentration range where Kn is directly proportional to [MTSEA], thus supporting our approach to measure modification rate using either Model I or II (see Materials and methods).
Next we took advantage of the fact that the reaction of MTS reagents with cysteine residues is strongly pH dependent, the reaction of MTS with protonated cysteine being 109-fold slower than the reaction with cysteine in a deprotonated form (Karlin and Akabas, 1998). The reaction rate km can thus formally be expressed as
where km0 is the reaction rate for a totally deprotonated cysteine, and pKa-cys the expected pKa value for cysteine ≈ 8.5. On the basis of the results presented in Fig. 9, equation Kn can now be rewritten as
where
Are MTSEA and K+ Ions Using the Same Diffusion Pathway?
The results obtained with MTSEA suggest that molecules of 4.6 Å in diameter can access the channel cavity in the closed configuration. This is clearly not supported by a model of the closed KCa3.1 based on the KcsA structure. However, with a pKa > 8.5, 6% of the MTSEA molecules are neutral at pH 7.4 and can possibly access cysteines at 275 by diffusing through the lipid phase. A series of experiments was thus performed in which the modification rate of V275C by MTSEA was measured in zero and 25 μM Ca2+ conditions in the presence of the quaternary ammonium ion TBA in the channel central cavity. TBA applied internally has been documented to cause a total block of the V275C channel with half inhibition at 20 μM (Banderali et al., 2004), and the presence of positively charged TBA in the channel cavity should affect the time course of the MTSEA action if MTSEA diffuses in a protonated form. Examples of inside-out recordings obtained in 1 mM TBA conditions are presented in Fig. 11 A. Our results indicate that the presence of TBA decreased the MTSEA modification rates of the V275C mutant computed for K0 from 840 ± 136 M−1s−1 (N = 7) in zero TBA to 48 ± 8 M−1s−1 (N = 2) in 1 mM TBA for the open channel, and from 112 ± 12 M−1s−1 (N = 7) in the absence of TBA to 40 ± 4 M−1s−1 (N = 4) with TBA for the closed channel configuration. Our results thus indicate that MTSEA was sensitive to the presence of TBA in the channel cavity when the channel was maintained in a closed state. More importantly, the rates of modification were identical for the open and closed configurations in the presence of TBA. In conditions where the rate of reaction km with TBA becomes rate limiting, so that K0 ≈ 4kikm[MTSEA]/ko, this observation would be indicative of a ratio ki/ko identical for the open and closed state, arguing for the charged form of MTSEA being in equilibrium between the bulk medium and the central cavity when the channel is closed. If alternatively the presence of TBA leads to km ≫ ko, so that K0 ≈ ki[MTSEA], our results would be indicative of TBA limiting the entry rate ki, independently of the channel state, thus supporting the proposal that MTSEA accesses the 275C residues for the closed channel in a charged form. This proposal is also in agreement with the results presented in curve C (empty squares) of Fig. 10 on the pH dependence of the modification rate by MTSEA of the closed V275C mutant with TBA in the channel central cavity. The resulting titration curve led to a pKa of 6.9. As the modification rate K0 = 0.040 s−1 for 1 mM [MTSEA] in the presence of TBA, Eq. 5 leads to ki ≈ 54 M−1s−1. This analysis indicates that the presence of TBA caused a ≈9-fold decrease of the MTSEA entry rate, thus supporting the proposal that MTSEA accesses the closed channel cavity in a charged form. In addition, Eq. 5 predicts that for low MTSEA concentrations, the pH dependence of the MTSEA modification rate in conditions where only the neutral form of MTSEA has access to the 275C residue is given by
where pKmtsea ≈ 8.5 (Karlin and Akabas, 1998), [MTSEA] total the total MTSEA concentration in the internal solution, and pKa-cys the pKa of the target cysteine residue (8.5). This expression predicts that the pH dependence of K0 should essentially reflect the titration curve of MTSEA with a pKa value ≈ 8.5 independently of the parameter β. The titration curves presented in Fig. 10 for the V275C mutant should have therefore been identical with a pKa close to 8.5. This is not in agreement with our experimental observations as we found a pKa of 7.8 in the absence (Fig. 10, curve B) and 6.9 in the presence (Fig. 10, curve C) of TBA for the modification rate by MTSEA of 275C. The presence of TBA is not expected to affect the pKa of MTSEA, indicating that the observed variation of modification rates as a function of pH are not related to the protonation state of MTSEA. Our results do not therefore support the proposal of MTSEA diffusing in a neutral form to the channel cavity. Finally, because MTSEA was documented to be membrane permeant and to cause transmembrane cysteine modification (Holmgren et al., 1996; Karlin and Akabas, 1998), control experiments were performed in which the modification rate of the V275C mutant was measured in zero Ca2+ conditions with 5 mM cysteine added to the pipette solution. There was no effect due to the presence of cysteine in the patch electrode on the modification rate of the V275C mutant by MTSEA (unpublished data), ruling out a potential contribution of a trans-inhibition coming from the diffusion of MTSEA through the membrane.
Can the Structure of the Closed KCa3.1 Channel distal to V282 be Represented by KcsA ?
The results obtained with MTSET on A283C and A286C mutants indicate that cysteines at 283 and 286 were accessible to MTSET in the closed channel configuration. The radial distribution of the α carbons for residues V275 to N292 presented in Fig. 1 A suggests however that A283 and A286 should participate in the formation of a bundle crossing region. Experiments were thus performed to investigate the topology of the pore region distal to V282 for the closed KCa3.1. In a first series of experiments, the accessibility of the 286C residue to the large MTS reagent MTS-PtrEA (space filling diameter of 9.5 and 12 Å long) was tested in zero Ca2+ conditions. The current recording presented in Fig. 11 B clearly shows that MTS-PtrEA could access the 286C residue with the channel in the closed configuration as the subsequent addition of Ca2+ revealed a strong channel activation as observed with MTSEA and MTSET. In a second series of experiments we took advantage of the fact that MTSET–cysteine complexes can be formed at positions 283 and 286 in zero internal Ca2+, thus adding one to four positively charged residues within the channel inner vestibule. On the basis of the model presented in Fig. 1 A, the addition of positively charged residues at position 283 or 286 is expected to decrease the entry rate of MTSEA within the channel cavity in zero internal Ca2+ conditions and affect the modification rate of the cysteine residues engineered at V275. The experimental protocol used for measuring the rates of modification for the V275C-A283C and V275C-A286C double mutant channels is presented in Fig. 11 C. The time course of the current inhibition could be fitted by a single exponential so that the modification rates were estimated directly from K0. The modification rates of 275C by MTSEA measured with MTSET–cysteine complexes at position 286 and 283 were estimated at 168 ± 8 M−1s−1 (N = 3) and 104 ± 4 M−1s−1 (N = 3), respectively, compared with 112 ± 12 M−1s−1 (N = 7) for the V275C simple mutant. The presence of positively charged residues within the channel inner vestibule for the closed KCa3.1 failed therefore to decrease the accessibility of MTSEA to cysteine residues located inside the channel central cavity. Finally, the KcsA-based structure predicted for the closed KCa3.1 illustrated in Fig. 12 A shows that the channel inner vestibule contains a double ring of histidine residues at positions 297 and 299. An additional histidine ring is also located at the level of the A286 residue coming from the four histidines in S5 at position 203. Because histidine residues are expected to be ∼90% protonated at pH 5.5, substituting histidine to alanine under these conditions should cause the neutralization of ≈11 equivalent positive charges in the channel inner vestibule. The results presented in Fig. 12 B on the modification rates estimated at pH 5.5 for the triple V275C-H297A-H299A and quadruple V275C-H297A-H299A-H203A mutant channels provide evidence that the neutralization of ≈11 equivalent positive charges at the inner entrance of the KCa3.1 does not affect the cysteine accessibility of MTSEA to 275C in the closed channel configuration.
Discussion
Two main conclusions follow from the observations reported in this work. First the S6 segment of the KCa3.1 channel contains two distinct functional domains delimited by the V282 residue with MTSEA and MTSET binding leading to a total channel inhibition for residues V275–V282 (V275, T278, and V282) and to a steep channel activation for residues distal to V282 (A283 and A286). Second, the pore structure of the closed KCa3.1 channel cannot be accounted for by the inverted teepee-like structure prevailing for the KcsA channel with a constriction at the level of V282 tight enough to be impermeable to molecules such as MTSEA, Et-Hg+, and Ag+. Rather, our results support a model where the dimensions of the pore between the cytosolic medium and central cavity for the closed channel range from 4.6 to 5.8 Å in diameter.
Structure of the Closed KCa3.1 Channel and Channel Gating
Our experiments on the closed KCa3.1 showed that the channel modification rates by MTSET were ≈103– 104 times faster for cysteine residues engineered at positions 283 and 286 compared with cysteines substituting the cavity-lining V275 residues. This effect was strongly size dependent as experiments performed with the smaller MTSEA led to modification rates differing by <20-fold between the A286C and V275C channels. These observations strongly argue for a pore structure of the closed KCa3.1 channel with a constriction located at the level of the T278–A282 residues tight enough to impair the passage of MTSET but leaky to MTSEA. Based on the cross section areas predicted for MTSET and MTSEA it follows that pore diameter for the closed channel between the cytosolic medium and the channel cavity should be ranging from 4.6 to 5.8 Å, which is smaller in size than a fully hydrated K+ ion (6 Å × 6 Å).
These conclusions are also in agreement with our observations that the modification rates of the closed V275C mutant by MTSEA, Et-Hg+, and Ag+ were poorly state dependent, in contrast to MTSET, which yielded modification rates at least 103-fold slower in zero than in saturating Ca2+. Notably, the modification rates by MTSEA of cysteines at position 282 were nearly identical for the open and closed states, suggesting a minimal structural change at this site in response to channel opening. Although modification rates with Ag+ could not be estimated for the open V275C mutant, a value of 3 × 107 s−1M−1 for the closed state is in agreement with the results reported for the closed/open CNG and Kir2.1 (IRK1J) channels (Flynn and Zagotta, 2003; Xiao et al., 2003) while being 100 times faster than the rates of modification measured in Shaker for the residue equivalent to V282 in KCa3.1 (del Camino and Yellen, 2001). In addition, a value of 3.1 ± 0.6 × 107 s−1M−1 is close to the diffusion-limited rate of modification expected for Ag+, arguing for an unrestricted diffusion of Ag+ in the channel cavity with V275C in the closed state. Our observations are thus in line with SCAM results reported for the SK2, CNG, and Kir2.1 channels where the helix bundle crossing region was described as being leaky to small sulphydryl reagents such as MTSEA and Ag+ but not MTSET (Flynn and Zagotta, 2003; Xiao et al., 2003). These authors concluded that the gate controlling K+ ion flow in ligand-gated channels should be located at the level or close to the selectivity filter. Their results and the observations presented in this work contrast with the findings reported on the voltage-gated Shaker channel where a huge decrease in the accessibility of the residues above the bundle crossing (V474 corresponding to V282 in KCa3.1, see Fig. 1 B) was measured with MTSET, MTSEA, and Ag+ when the channel was in a closed configuration (Liu et al., 1997; del Camino and Yellen, 2001). The structural changes involved in the opening of the pore thus appear to differ considerably between KCa3.1 and voltage-gated channels such as Shaker.
Can MTSEA Be Trapped in the Channel Cavity?
Because the MTSEA results set the lower limit of the pore diameter for the closed channel to 4.6 Å, which is nearly two times larger than the van der Walls diameter (2.6 Å) of a K+ ion but less than the diameter of a fully hydrated K+ (6 Å), alternative interpretations for the weak state dependence observed with the V275C mutant need to be considered. One of these interpretations refers to the trapping of MTSEA inside the channel cavity (Phillips et al., 2003). The kinetic scheme related to MTSEA trapping can be expressed as (see Phillips et al., 2003)where α and υ correspond to the average channel opening and closing rates, OM and CM the kinetic states associated to the open and closed channel configurations with MTSEA in the binding site, and CB, OB the state following the irreversible binding of MTSEA to the target cysteines. We found from single channel recordings that the mean open time of the V275C mutant is of the order of 16 ms and varies little as a function of the internal Ca2+ concentration (unpublished data). This observation sets the upper limit for the closing rate υ to 63 s−1. In addition, the results presented in Fig. 9 indicate that the reaction rate, km, for the open V275C mutant is approximately equal to 7 s−1 with [MTSAE]1/2 ≈ 35 mM, and our pH analysis showed that in the closed configuration the MTSEA entry rate for the V275C mutant corresponds to ki ≈ 480 M−1s−1. This value constitutes a lower limit, the entry rate being expected to be higher for the open channel configuration. With ki set to 104 M−1s−1 and km = 7 s−1, the expression derived for [MTS]1/2 (Appendix, Eq. A5) predicts an exit rate ko ≈ 320 s−1. The variation in modification rate as a function of the channel open probability can be estimated by solving the kinetic scheme 3 numerically with υ set to 63 s−1. This analysis led to the conclusion that a trapping mechanism would require an open probability Po > 0.05 to account for a measured inhibition rate of 0.1 s−1 in zero Ca2+ for 1 mM MTSEA. An identical result was obtained assuming an entry rate ki = 105 M−1s−1. If ki is presumed to be the same for the open and closed state (ki ≈ 480 M−1s−1), a Po > 0.1 would then be required to account for our results. A Po of 0.05 corresponds to the level of channel activity expected for a 0.6 μM internal Ca2+ solution, which is at least 50 times higher than the expected Ca2+ concentration in our EGTA-buffered Ca2+-free solutions (Simoes et al., 2002). A smaller Po of 0.02 would however be required if the mean open time of the V275C mutant is set to 0.1 or 1 ms with ki > 104 M−1s−1. This level of channel activity still corresponds to an internal Ca2+ concentration of 0.5 μM, arguing also against a trapping mechanism. Finally, Po ranging from 0.02 to 0.1 in the absence of Ca2+ would have resulted in current level for the V275C mutant clearly distinct from the current level after MTSEA block. Such a behavior was not observed experimentally. The absence of state-dependent modifications by MTSEA of the V275C mutant is not likely therefore to be accountable by a trapping mechanism.
An alternative explanation to the observation that MTSEA has access from the cytoplasmic side to the channel cavity with KCa3.1 in a closed configuration would be that MTSEA uses a diffusion pathway not applicable to K+ ions and/or diffuses through the channel pore in a deprotonated form. Evidence arguing against this interpretation comes from the TBA protection experiments and from the pH dependence of the 275C modification rate by MTSEA, which support a mechanism where MTSEA accesses the channel central cavity from the cytosolic medium as a protonated cation. Our MTSEA results are thus compatible with a model where the weak state dependence of the 275C modification rate by MTSEA is likely to reflect a pore structure for the closed KCa3.1 with a pore constriction tight enough to impair the passage of MTSET but leaky to smaller reagents such as MTSEA, Ag+, and Et-Hg+.
Toward a 3D Representation of the Closed KCa3.1 Channel Structure
The results on the 275C accessibility to MTSEA in zero Ca2+ also need to be interpreted in light of our findings on the dimensions of the closed KCa3.1 inner vestibule. Our results show for instance that the large MTS-PtrEA had access to the 286C residue in zero Ca2+. Because a MTS-PtrEA molecule fits into a cylinder 9.5 Å diameter and 12 Å long, the functional pore dimensions for the closed channel below A286 should be >9.5 Å comparatively to 4.2 Å as predicted on the basis of a KcsA-like structure. Furthermore, the addition of positively charged MTSET–cysteine complexes at position 283 did not impair MTSEA from diffusing inside the channel cavity in the closed configuration. With a predicted pore diameter of 6.2 Å at A283, a KcsA-based representation of the closed KCa3.1 can hardly account for these observations. Taken together, our analysis does not support a model for the closed KCa3.1 channel where the cytoplasmic pore structure would be represented by an inverted tepee-like structure as described for KcsA. In accordance with this proposal, the neutralization at pH 5.5 of 11 equivalent positive charges in the V275C-H203A-H297A-H299A mutant did not alter the modification rates by MTSEA of the cysteines at 275. This observation would be compatible with a model where the residues distal to V282 residues are connected to rather wide internal vestibule. Alternatively, the secondary structure of the S6 segment distal to 283 may not be totally represented by an α helix as predicted on the basis of a KcsA structure but may contain some turn/coil regions as observed with MthK channel, so that the residues predicted to be facing the channel pore may be oriented in another direction. This could explain the inefficiency of the MTSET–cysteine complexes at 283 to impair the access of MTSEA to the channel cavity with the channel closed.
Gating of KCa3.1 Channel, a Multilocation Process
The extent of the contribution of the bundle crossing region to the ligand-dependent gating of ion channels is not fully established. Recent results obtained with the pH-dependent KcsA and Ca2+-activated MthK channels argue for a multiple gate system where the stimulus-dependent gate at bundle crossing would be determinant to the duration of individual opening channel bursts with intraburst fluctuations governed by a second gating process located at the level of the selectivity filter or close to it (Cordero-Morales et al., 2006; Zadek and Nimigean, 2006). A recent study where the opening of the bundle crossing in KcsA was monitored using a fluorescence-based approach also concluded in the presence of a second gate at the level of the selectivity filter (Blunck et al., 2006). The finding that the accessibility of cysteine residues at position 275 was poorly state dependent when MTSEA, Et-Hg+, and Ag+ were used as sulphydryl reagents does not favor a model where the constriction region from V282 to A286 would act as an activation gate. This conclusion is in accordance with our results demonstrating that the closed KCa3.1 cannot be represented by a KcsA-like structure. The possibility that the selectivity filter region contributes to the activation of ligand-gated channels has been discussed in several studies (Bruening-Wright et al., 2002) and a similar mechanism could prevail in KCa3.1. The S6 segments would then act as transducers conveying the structural changes from the bundle crossing to the selectivity filter. Our results demonstrate nevertheless that the action of Ca2+ on KCa3.1 is consistent with an opening of the pore at the C-terminal end of S6. This conclusion is mainly supported by the observation that the accessibility of the cysteines at 275 to MTSET is state dependent. It remains to be determined to what extent an increase in pore dimensions corresponding at least to the diameter of an MTSET molecule (5.8 Å) or larger is essential for channel gating. Opening of the pore in the C terminus of S6 may serve other functions such as to modify the geometry of the cavity and/or the selectivity filter region for gating, or favor the diffusion of K+ ions through the V282 constriction region thus increasing the channel conductance (Li and Aldrich, 2006). The observation that MTSEA and MTSET could activate the A283C and A286C mutants also raises the possibility of a contribution of the S6 segment from A283 to A286 to channel gating distinct from a gating process involving the selectivity filter and/or the cavity region. Several single channel analyses have already demonstrated that the KCa3.1 channel open probability is smaller than one at saturating internal Ca2+ concentrations, with typical values of 0.2–0.4 (Sauvé et al., 1986; Keen et al., 1999). The fact that MTSET and MTSEA acting on A283C or A286C caused an increase of the channel maximum open probability at saturating Ca2+ therefore suggests additional gating mechanisms that are distinct from a Ca2+-dependent opening of the channel pore at the C-terminal end of S6. These possibilities need further investigations.
Conclusions
Our results do not support a KcsA-like structure for the pore region of the closed KCa3.1 channel. It is therefore unlikely that a constriction at the C-terminal end of S6 contributes to the gating process of KCa3.1 by Ca2+. Our results point nevertheless toward a pore structure with a diameter d such that 4.6 Å < d < 5.8 Å connecting the channel central cavity and a wide cytoplasmic vestibule.
Appendix
The time course of cysteine modification predicted from the kinetic Scheme 1 presented in Materials and methods is given by
where
and
In conditions where Ks ≫ Kn, the expression (A1) reduces to
with Kn now expressed as
Eq. A5 corresponds to a standard Michaelis Mentel equation with a maximum modification rate given by (4 − n)km and a [MTS] concentration for half maximum rate ([MTS]1/2) equal to
so that Kn is now a linear function of [MTS] with a slope given by
The modification rate per cysteine KC in the kinetic Model I presented in Materials and methods thus formally corresponds to KC = kmki/ko with
Acknowledgments
We thank Drs. Daniel Devor and Benoit Roux for stimulating discussions. We also acknowledge the work of Ms. Julie Verner for expert oocyte preparation.
This work was supported by grants from the Canadian Institutes of Health Research (MOP 7769), from the Canadian Heart and Stroke Foundation, and from the Canadian Cystic Fibrosis Foundation.
Olaf S. Andersen served as editor.
References
Abbreviations used in this paper: MTS, methanethiosulfonate; MTSET, [2-(trimethylammonium)ethyl] methanethiosulfonate bromide; MTSEA, 2-aminoethyl methanethiosulfonate hydrobromide; MTS-PtrEA, [3-(triethylammonium)propyl] methanethiosulfonate; Et-Hg+, ethylmercury; MTSACE, 2-(aminocarbonyl)ethyl methanethiosulfonate; MTSES, sodium (2-sulfanoethyl) methanethiosulfonate.