Calcium-driven regulation of voltage-sensing domains in BK channels

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: September 26, 2019 (This version)
Accepted Manuscript published: September 11, 2019 (Go to version)
Accepted: September 10, 2019
Received: January 7, 2019

1. Of interest
The substrate-binding domains of the osmoregulatory ABC importer OpuA transiently interact

Marco van den Noort, Panagiotis Drougkas ... Bert Poolman

Research Article May 2, 2024
Further reading

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

Allosteric interactions between the voltage-sensing domain (VSD), the Ca²⁺-binding sites, and the pore domain govern the mammalian Ca²⁺- and voltage-activated K⁺ (BK) channel opening. However, the functional relevance of the crosstalk between the Ca²⁺- and voltage-sensing mechanisms on BK channel gating is still debated. We examined the energetic interaction between Ca²⁺ binding and VSD activation by investigating the effects of internal Ca²⁺ on BK channel gating currents. Our results indicate that Ca²⁺ sensor occupancy has a strong impact on VSD activation through a coordinated interaction mechanism in which Ca²⁺ binding to a single α-subunit affects all VSDs equally. Moreover, the two distinct high-affinity Ca²⁺-binding sites contained in the C-terminus domains, RCK1 and RCK2, contribute equally to decrease the free energy necessary to activate the VSD. We conclude that voltage-dependent gating and pore opening in BK channels is modulated to a great extent by the interaction between Ca²⁺ sensors and VSDs.

https://doi.org/10.7554/eLife.44934.001

Introduction

Diverse cellular events involve calcium ions as a primary mediator in the signal transduction pathways triggering, among other signaling processes, Ca²⁺-activated conductances. Since the BK channels are regulated by cytosolic Ca²⁺ and depolarizing voltages (Marty, 1981; Pallotta et al., 1981; Latorre et al., 1982), they are integrators of physiological stimuli involving intracellular Ca²⁺ elevation and membrane excitability. BK channels are modular proteins in which each module accomplishes a specific function. Thus, different modules harbor voltage and Ca²⁺ sensors that communicate allosterically with the channel gate (Cox et al., 1997; Horrigan and Aldrich, 1999; Horrigan and Aldrich, 2002; Horrigan et al., 1999; Rothberg and Magleby, 1999; Rothberg and Magleby, 2000; Cui and Aldrich, 2000). Functional BK channels are formed by homotetramers of α-subunits (Shen et al., 1994), each comprising a transmembrane voltage-sensing domain (VSD) and an intracellular Ca²⁺-sensing C-terminal domain (CTD) that can independently modulate the ion conduction gate in the pore domain (PD) (Latorre et al., 2017). The CTDs consist of two non-identical regulators of the conductance of K⁺ domains (RCK1 and RCK2) arranged in a ring-like tetrameric structure dubbed the gating ring (Wu et al., 2010; Yuan et al., 2010; Yuan et al., 2012; Hite et al., 2017; Tao et al., 2017). Each RCK domain contains distinct ligand-binding sites capable of detecting Ca²⁺ in the micromolar range (Schreiber and Salkoff, 1997; Bao et al., 2002; Xia et al., 2002).

In the absence of Ca²⁺, the activation of the VSD decreases the free energy necessary to fully open the BK channels in an allosteric fashion (Horrigan and Aldrich, 1999; Horrigan et al., 1999). With 0 Ca²⁺, very positive membrane potentials are required to drive all voltage sensors to their active conformations (Cui et al., 1997; Stefani et al., 1997; Horrigan et al., 1999; Contreras et al., 2012), ultimately leading to near-maximal BK activation. Thus, in cells such as neurons, an appreciable open probability of BK channels at physiologically relevant voltages necessarily requires the activation of the Ca²⁺ sensors on the gating ring. The allosteric interplays established between the functional and structural modules (VSD-PD, CTD-PD, and CTD-VSD) are key for enabling BK channels to operate over a wide dynamic range of internal Ca²⁺ and voltage, thereby fine-tuning the channel’s gating machinery. Therefore, understanding the structure-functional bases that underlie the Ca²⁺ and voltage activation mechanisms interrelationship becomes essential to unveil how the channel behaves under different physiological conditions.

The voltage dependence of Ca²⁺-dependent gating ring rearrangements (Miranda et al., 2013; Miranda et al., 2018) and RCK1 site occupancy (Sweet and Cox, 2008; Savalli et al., 2012; Miranda et al., 2018) as well as the perturbation of VSD movements by Ca²⁺ binding (Savalli et al., 2012) support the idea that the energetic interaction between both specialized sensors may be crucial for BK channel activation. The physical CTD-VSD interface has been suggested to provide the structure capable of mediating the crosstalk between these sensory modules and their synergy in activating the pore domain (Yang et al., 2007; Sun et al., 2013; Tao et al., 2017; Zhang et al., 2017). However, the strength of the interaction between voltage and Ca²⁺ sensors and their relevance to BK channel activation is still debated (cf. Horrigan and Aldrich, 2002; Carrasquel-Ursulaez et al., 2015). Also, the functional role that each of the high-affinity Ca²⁺-binding sites plays in the CTD-VSD allosteric interaction is an open question. The RCK1 and RCK2 Ca²⁺-binding sites have distinct functional properties conferred by their different molecular structures and relative positions within the gating ring (Wu et al., 2010; Yuan et al., 2010; Yuan et al., 2012; Hite et al., 2017; Tao et al., 2017). Thus, the RCK sites differ in their Ca²⁺ binding affinities (Bao et al., 2002; Xia et al., 2002; Sweet and Cox, 2008), divalent cations selectivity (Oberhauser et al., 1988; Schreiber and Salkoff, 1997; Zeng et al., 2005; Zhou et al., 2012), voltage dependence (Sweet and Cox, 2008; Savalli et al., 2012; Miranda et al., 2018) and in their contribution to allosteric gating mechanisms (Yang et al., 2010; Yang et al., 2015). In particular, only the RCK1 site appears to be involved in communicating the Ca²⁺-dependent conformational changes towards the membrane-spanning VSD (Savalli et al., 2012; Miranda et al., 2018). Recently, the Aplysia BK structure has revealed that the N-lobe of the RCK1 domain is in a non-covalent contact with the VSD and the S4-S5 linker that connects the voltage sensor to the pore domain. This RCK1-VSD interaction surface is rearranged when comparing the liganded and Ca²⁺-free structures (Hite et al., 2017; Tao et al., 2017). In fact, it has been hypothesized that any Ca²⁺-induced rearrangements of the gating ring should ultimately be transmitted to the pore domain via the VSD (Hite et al., 2017; Zhou et al., 2017). Thereby, defining the extent to which Ca²⁺ binding influences to VSD is central for determining the importance of the crosstalk between sensors in decreasing the free energy necessary to open the BK channel.

Here, we examined the Ca²⁺-dependence of the VSD activation by estimating the allosteric coupling between the Ca²⁺ and voltage sensors. By analyzing gating currents under unliganded and Ca²⁺-saturated conditions, we found a strong energetic influence of Ca²⁺-binding on voltage sensor equilibrium in a manner that is independent of channel opening. These findings show that a major component in the synergistic Ca²⁺ and voltage activation of BK channels resides in the strong coupling between Ca²⁺ binding and the voltage sensor activation. We also found that the Ca²⁺-dependence of voltage sensor activation is consistent with a CTD-VSD allosteric coupling that occurs through a concerted interaction scheme in which each Ca²⁺ bound to one subunit affects all voltage sensors in the BK tetramer equally. Notably, we found that the two distinct RCK1 and RCK2 Ca²⁺ sensors contribute equally to the VSD activation via independent allosteric pathways.

Results

Allosteric coupling between Ca²⁺-binding and voltage sensor activation is strong

We characterized the effects of Ca²⁺-binding on voltage sensor activation in BK channels by analyzing the gating currents measured in inside-out patches of Xenopus laevis oocyte membrane. The amount of gating charge displaced ( $Q_{C}$ ) at each Ca²⁺ concentration was obtained by integrating the initial part of the decay of the ON-gating current (I_G-ON), which was fitted to a single exponential (fast ON-gating; see Materials and methods). As we show below, we determined only the gating charge displaced before the opening of the BK channel. Figure 1A–B show representative I_G-ON records in response to 160 mV voltage step in the nominal absence of Ca²⁺ (‘zero’ Ca²⁺) and in saturating Ca²⁺ concentration (100 µM Ca²⁺). In Figure 1A–B, we also show the initial time courses of the corresponding macroscopic K⁺ current (I_K) activation at the same voltage and internal Ca²⁺ conditions. The I_G-ON relaxation exhibits an almost complete decay before the I_K achieves an exponential time course in 0 Ca²⁺ conditions (Figure 1A). Thus, the time constant of the I_K activation (~3.4 ms at 160 mV) following a delay of ~160 µs is consistent with the movement of the voltage sensors preceding channel opening. Under saturating internal Ca²⁺ conditions, the I_G-ON time course develops 20 times faster than the exponential kinetic of the I_K activation ( $τ_{I_{G} - O N}$ = 30 µs and $Δ t_{I_{K}}$ = 660 µs at 160 mV) and is also almost complete within the time interval comprised by the I_K delay ( $Δ t_{I_{K} (100 μ M)}$ = 84 µs) (Figure 1B). Thus, the fast I_G-ON relaxation reflects the movement of the gating charge in the channel’s closed conformation, regardless of internal Ca²⁺ concentration.

Figure 1

Download asset Open asset

Ca²⁺ binding strongly affects the activation of VSD in BK channels.

(**A–B**) Fast component of the ON-gating current (*I_G*-ON) at 0 and 100 μM internal Ca²⁺ concentration, respectively. The representative *I_G*-ON records were evoked by applying a 160 mV voltage step of 1 ms duration. The first 100 µs of the *I_G*-ON were fitted to a single exponential function (red line; $τ_{I_{G} - O N (0 C a^{2 +})}$ = 58 µs and $τ_{I_{G} - O N (100 μ M)}$ = 30 µs). The area under the curve described by the monoexponential fit (gray areas) was integrated to obtain the charge displaced between closed states ( $Q_{C}$ ). For comparison, the initial time course of the macroscopic K⁺ current (*I_K*) activation (solid black line) obtained at 0 and 100 μM Ca²⁺ were superimposed on the *I_G*-ON under the same internal Ca²⁺ conditions. The *I_K* were evoked by a 15 ms pulse to 160 mV. The *I_K* is described by an exponential function with initial delay: $Δ t_{I_{K} (0 C a^{2 +})}$ = 0.16 ms and $τ_{I_{K} (0 {C a}^{2 +})}$ = 3.37 ms at 0 Ca²⁺; and $Δ t_{I_{K} (100 μ M)}$ = 0.08 ms and $τ_{I_{K} (100 μ M)}$ = 0.66 ms at 100 μM Ca²⁺. After 1 ms, the *I_K* increased to about 20% and 80% of its steady-state amplitude in 0 Ca²⁺ and 100 μM Ca²⁺, respectively. (C) Voltage-dependence of the $Q_{C}$ and of the gating current time constants ( $τ_{I_{G} - O N}$ ) at 0 Ca²⁺ (open circles) and 100 μM Ca²⁺ (filled circles). Gating charge-voltage relationships ( $Q_{C} (V) / Q_{C, M A X}$ ) were obtained by integrating the fast component for each ON *I_G* trace (from −90 mV to 350 mV). Boltzmann fitting to the experimental data (mean ± SEM) is indicated by solid line at ‘zero’ Ca²⁺ ( $V_{H}$ = 174.5 ± 2.4 mV, $z_{Q}$ = 0.60 ± 0.01, $n$ = 25) and by a dashed line at 100 µM Ca²⁺ ( $V_{H}$ = 31.9 ± 4.5 mV, $z_{Q}$ = 0.66 ± 0.01, $n$ = 7). Right ordinate shows the time constants data (mean ± SEM) of the exponential decays of *I_G*-ON ( $τ_{I_{G} - O N}$ ) plotted against the voltage. The best fit to a two-state model of the VSD activation where the $z_{J}$ was constrained to the values found for the $Q_{C} (V)$ relation ( $z_{Q}$ ) is indicated by solid lines at ‘zero’ Ca²⁺ ( $α_{0}$ = 3.73 ms⁻¹, $β_{0}$ = 76.10 ms⁻¹ and $δ$ = 0.29) and dashed line at 100 µM Ca²⁺ ( $α_{0}$ = 7.28 ms⁻¹, $β_{0}$ = 6.98 ms⁻¹ and $δ$ = 0.36). The corresponding $V_{H}$ at each internal Ca²⁺ concentration is indicated by a vertical line. (D) Semi-logarithmic plot of the $Q_{C} (V) / Q_{C, M A X}$ curves at 0 Ca²⁺ and saturating Ca²⁺ concentration (100 µM). The allosteric parameter $E$ determines the vertical displacement of the 0 mV intercept (dashed vertical line) of the $Q_{C} (V) / Q_{C, M A X}$ curve at 100 µM Ca²⁺ relative to the 0 Ca²⁺ condition.

https://doi.org/10.7554/eLife.44934.002

The $Q_{C}$ was determined over a wide range of membrane potentials at low and high internal (100 μM) Ca²⁺ concentration. Increasing the Ca²⁺ concentration promoted a large leftward shift ( ${∆ V}_{H}$ = -142.6 ± 4.5 mV) of the normalized $Q_{C}$ ( $Q_{C} (V) / Q_{C, M A X}$ curves) (Figure 1C). In spite of with the large leftward shift of the $Q_{C}$ at high Ca²⁺ concentrations, we found no appreciable slow component in I_G-ON (cf. Horrigan and Aldrich, 2002; see also Figure 2—figure supplement 1A). As expected, a large Ca²⁺-dependent leftward shift was also observed in the time constants of the exponential decays of I_G-ON ( $τ_{I_{G} - O N}$ ). The $τ_{I_{G} - O N} (V)$ curves were fitted to a two state model ( $τ (V) = 1 / (α (V) + β (V)$ ) where the forward ( $α$ ) and backward ( $β$ ) rate constants represent the resting-active (R-A) transitions of the voltage sensors, that determine the equilibrium constant of the VSD activation ( $J (V) = α (V) / β (V)$ ). The predominant effect of Ca²⁺-binding on VSD activation appears to cause a decrease in the backward rate constant at zero voltage ( $β_{0}$ ; see Materials and methods). $β_{0}$ decreases from 76 ms⁻¹ at 0 Ca²⁺ to 7 ms⁻¹ at 100 μM Ca²⁺, which results in a shift in the equilibrium of the voltage sensors towards their active conformation. Such a large shift ( $Δ V_{H}$ = -142.6 ± 4.5 mV) implies that Ca²⁺ binding to the RCK Ca²⁺-binding sites alters the VSD equilibrium, which in consequence causes a decrease in the free energy ( $Δ Δ G_{V}^{C a}$ ) that defines the voltage sensor R-A equilibrium ( $J$ ) by ~8 kJ/mol ( $Δ Δ G_{V}^{C a}$ = -7.98 ± 0.27 kJ/mol). In terms of the allosteric gating scheme (Horrigan and Aldrich, 2002), this result means that the $J$ equilibrium constant of VSD becomes amplified by an allosteric factor $E$ equal to 26.4 at Ca²⁺-saturated conditions (Figure 1D), revealing a strong allosteric coupling between Ca²⁺-binding sites and voltage sensors.

Families of gating currents (I_G) were evoked at different intracellular Ca²⁺ concentrations ([Ca²⁺]_i) ranging from 0.1 to 100 µM in K⁺-free solution (Figure 2A). For all experiments, we first measured I_G in 'zero' Ca²⁺ condition and then perfused the internal side with solutions containing different concentrations of Ca²⁺. The increase in internal Ca²⁺ promoted a leftward shift of the $Q_{C}$ versus voltage ( $Q_{C} (V)$ ) curves (Figure 2B-C), which indicates that Ca²⁺-binding facilitates the activation of the voltage sensor, being more prominent as Ca²⁺-binding site occupancy increases.

Figure 2 with 2 supplements see all

Download asset Open asset

Ca²⁺-dependent effects on VSD activation in BK channels.

(A) Representative gating current (*I_G*) recordings at different internal Ca²⁺ concentrations (from 0 to 100 μM). *I_G* was evoked by the indicated voltage protocol of 1 ms duration. Representative gating current records are from different patches with the exception of 0 and 100 μM Ca²⁺. (B) Gating charge-voltage relationships ( $Q_{C} (V)$ ) were obtained by integrating the fast component for each ON *I_G* trace. Normalized gating charge data $(Q_{C} (V) / Q_{C, M A X})$ (mean ± SEM) were fitted using a single Boltzmann function (solid lines). The fit parameters are as follows: ‘zero’ Ca²⁺ ( $V_{H}$ = 174.5 ± 2.4 mV, $z_{Q}$ = 0.60 ± 0.01, $n$ = 25); 0.1 µM Ca²⁺ ( $V_{H}$ = 162.4 ± 4.2 mV, $z_{Q}$ = 0.59 ± 0.01, $n$ = 5); 0.5 µM Ca²⁺ ( $V_{H}$ = 151.6 ± 1.3 mV, $z_{Q}$ = 0.60 ± 0.01, $n$ = 5); 1 µM Ca²⁺ ( $V_{H}$ = 137.1 ± 5.1 mV, $z_{Q}$ = 0.61 ± 0.01, $n$ = 5); 5 µM Ca²⁺ ( $V_{H}$ = 121.9 ± 3.8 mV, $z_{Q}$ = 0.63 ± 0.01, $n$ = 5); 10 µM Ca²⁺ ( $V_{H}$ = 67.3 ± 17.1 mV, $z_{Q}$ = 0.64 ± 0.08, $n$ = 4); 100 µM Ca²⁺ ( $V_{H}$ = 31.9 ± 4.5 mV, $z_{Q}$ = 0.66 ± 0.01, $n$ = 7). (C) $V_{H}$ obtained from the $Q_{C} (V)$ curves as a function of Ca²⁺ concentration (mean ± SEM). Ca²⁺ binding produces a leftward shift in $V_{H} (Δ V_{H})$ : 0.1 µM Ca²⁺ ( $Δ V_{H}$ = −12.1 ± 3.5 mV); 0.5 µM Ca²⁺ ( $Δ V_{H}$ = −22.9 ± 1.8 mV); 1 µM Ca²⁺ ( $Δ V_{H}$ = −37.1 ± 3.5 mV); 5 µM Ca²⁺ ( $Δ V_{H}$ = −50.3 ± 4.7 mV); 10 µM Ca²⁺ ( $Δ V_{H}$ = −107.1 ± 17.1 mV); 100 µM Ca²⁺ ( $Δ V_{H}$ = −142.6 ± 4.5 mV). One-way ANOVA followed by Dunnett’s post-hoc test analysis was used to assess statistical significance of the Ca²⁺-induced shifts in $V_{H}$ (***p<0.001).

https://doi.org/10.7554/eLife.44934.003

Figure 2A shows that the OFF-gating current was dramatically modified becoming smaller in amplitude and with slower kinetics as the internal Ca²⁺ concentration was increased (see also Figure 2—figure supplement 1A and Figure 2—figure supplement 2A). At least two components (fast and slow components) could be resolved in the OFF gating current decay in 'zero' and 10 µM internal Ca²⁺ conditions (Figure 2—figure supplement 1B-C). In Figure 2—figure supplement 1D, we fitted $τ_{I_{G} - O N} (V)$ data for 0 and 10 µM internal Ca²⁺ to a two-state R-A model. Figure 2—figure supplement 1D shows that the time constants of the fast relaxation of the I_G-OFF evoked at -90 mV (open and closed orange circles) are in reasonable agreement with the time constants values extrapolated to the same voltage from the two-state model of the I_G-ON kinetic. The change in the time constant of the fast component of OFF-gating current detected at 10 µM Ca²⁺ relative to 'zero' Ca²⁺ condition agrees with a Ca²⁺-induced effect primarily on the backward rate of the R-A transition ( $β_{0}$ decreases from 76 ms⁻¹ at 0 Ca²⁺ to 10.6 ms⁻¹ at 10 μM Ca²⁺, see Figure 2—figure supplement 1D), which is consistent with the large $Q_{C} (V)$ shift observed ( $Δ V_{H (10 μ M)}$ = -107.1 ± 17.1 mV). On the other hand, the relative contribution of the slower component increasing as internal Ca²⁺ is increased, reflecting an increase in the open probability of the channel (Figure 2—figure supplement 1E-F). This kinetic behavior recapitulates the effect described on gating charge displacement as a function of the depolarizing pulse duration (Horrigan and Aldrich, 2002; Carrasquel-Ursulaez et al., 2015; Contreras et al., 2012), and confirms that this phenomenon is associated with the time course of channel opening revealing the allosteric interaction between voltage sensors and the pore gate (Horrigan and Aldrich, 2002).

Since the OFF gating charge cannot be accurately estimated from the OFF gating current recorded at -90 mV, we performed two experiments using 100 µM internal Ca²⁺ concentration and recorded the OFF gating currents at -150 mV and -200 mV (Figure 2—figure supplement 2A). At these applied voltages, the charge displaced in the ON is recovered in the OFF (Figure 2—figure supplement 2B). As expected from the large shift to the left along the voltage axis of $Q_{C}$ at saturating internal Ca²⁺ concentrations, there was no difference between the $Q_{C}$ curve and the gating charge-voltage curves obtained from the OFF gating currents 2 ms after the onset of the voltage pulse, for applied voltages of -150 or -200 mV (Figure 2—figure supplement 2C).

Ca²⁺ binding to a single α-subunit affects the R-A voltage sensor equilibrium of all four subunits equally

Taking advantage of the dose-dependent effect of Ca²⁺ on voltage sensor activation we investigated the underlying mechanism of the communication between the Ca²⁺ binding and voltage sensors in the context of the well-established Horrigan-Aldrich (HA) allosteric gating model (Horrigan and Aldrich, 2002). Two different mechanisms were proposed by Horrigan and Aldrich for the interaction between the Ca²⁺-binding sites and voltage sensors. The first mechanism proposes that Ca²⁺ binding to one α-subunit only affects the VSD in the same subunit (Scheme I) (Figure 3A), whereas the second mechanism proposes that the Ca²⁺ binding affects all four VSD equally (Scheme II) (Figure 3B). It should be noted that the standard HA model makes two simplifying assumptions. First, the model considers that there is a single Ca²⁺-binding site per α-subunit. Second, the model assumes the Scheme I as the Ca²⁺ binding-VSD interaction scheme underlying the general gating mechanism of BK channel (Horrigan and Aldrich, 2002).

Figure 3 with 2 supplements see all

Download asset Open asset

Model-dependent behavior of the $Q_{C} (V)$ curves based on the CTD-VSD interaction mechanisms according to the fractional occupancy of Ca²⁺-binding sites.

(**A–B**) Cartoons representing two interaction schemes between voltage sensors and Ca²⁺-binding sites (modified from Horrigan and Aldrich, 2002). Scheme I (A) assumes that Ca²⁺-binding only affects the voltage sensor of one α-subunit ( $E_{M 1}$ ), whereas Scheme II (B) predicts that binding of Ca²⁺ to one α-subunit will affect VSD in all subunits equally, increasing the voltage sensor equilibrium constants ( $J$ ) $E_{M 2}$ -fold in all four subunits ( $E_{M 2}^{4}$ , when the four Ca²⁺ sites are occupied). In both schemes, a single Ca²⁺-binding site is considered in each α-subunit. (**C–D**) Predictions of $Q_{C} (V)$ relationships at different internal calcium concentration (from 0 to 10 mM) by two distinctive interaction mechanisms between Ca²⁺-binding sites and voltage sensors (Scheme I and Scheme II), respectively. $Q_{C} (V)$ curves were generated using Equation 4 (blue: Model I) or Equation 6 (red: Model II), and the following set of parameters: $z_{J}$ = 0.58, $J_{0}$ = 0.018, $K_{D}$ = 11 µM and $E_{M 1}$ = $E_{M 2}^{4}$ = 25.

https://doi.org/10.7554/eLife.44934.006

We simulated the normalized $Q_{C} (V)$ curves over a wide range of Ca²⁺ concentrations (from 0 to 10 mM) for each Ca²⁺-VSD interaction scheme (Figure 3C-D). We assumed that the fast gating currents measured correspond to the charge displaced by the R-A transitions and does not include the charge associated with the transition between the activated states (see Figure 3—figure supplement 1B-C). In other words, Equations 4 and 6 given in the Appendix describing the $Q_{C} (V)$ relationships for Schemes I and II are based on the assumption that channels equilibrate between Ca²⁺-bound states while the channel is closed and voltage sensors are not activated. We further assume that this distribution is not altered when gating current is measured because Ca²⁺ binding is slower than voltage-sensor activation. These assumptions are reasonable since the Ca²⁺-binding rate constant estimated for BK is 1.8 x 10⁸ M⁻¹s⁻¹ (Hou et al., 2016) implying that at 10 µM internal Ca²⁺ (the highest non-saturating Ca²⁺ concentration tested) the time constant of the Ca²⁺ binding is ~300 µs, while the VSD activates with a time constant of ~60 μs at $V$ = $V_{H}$ (see the Appendix for details of the simulations).

At extreme conditions of low (0.03 to 0.1 µM) and high (≥100 µM) internal Ca²⁺, VSD activation behaves in a mechanism-independent manner since all voltage sensors are in the same functional state (unliganded or saturated). However, the distinctive effects on $Q_{C} (V)$ curves at intermediate Ca²⁺ concentrations (1-10 µM) provide useful signatures to distinguish between the two mechanisms. Indeed, Scheme I predicts two functional states of the VSD depending on the occupancy status of the Ca²⁺ site (Ca²⁺ bound and unbound) such that the $Q_{C} (V)$ curve is described by the fractional distribution of the unliganded and Ca²⁺-saturated functional states like an all-or-none allosteric effect (Figure 3C; Figure 3—figure supplement 1B and Equation 4 in Appendix). By contrast, the Ca²⁺-binding effect on the VSD activation according to Scheme II is characterized by a five-component Boltzmann function (Figure 3—figure supplement 1C and Equation 6 in Appendix). Each component represents a single functional state determined by the number of Ca²⁺ bound to the channel (from 0 to 4). In such a case, the $Q_{C} (V)$ curves resulting from a distribution of functional states are equivalent to a single Boltzmann function, leftward shifted by an incremental allosteric effect (from $E$ to $E^{4}$ ) as the number of Ca²⁺ bound to the channel increases (Figure 3D). The experimental leftward shift of the $Q_{C} (V)$ curves occur with constant slope ( $z_{Q}$ ) (Figure 2B), which is consistent with Model II. It can be argued that Scheme I can produce single Boltzmann curves without a change in $z_{Q}$ if Ca²⁺ re-equilibration is fast enough. However, Figure 3—figure supplement 2 demonstrates that to recover the experimental $z_{Q}$ requires a Ca²⁺ binding rate constant that is 100-fold faster than the one reported by Hou et al. (2016) exceeding the diffusion limit constraint. Therefore, we can conclude that Ca²⁺-binding is slow enough to safely ignore calcium re-equilibration during the 100 μs it takes to measure the gating charges.

To further elucidate the mechanism by which Ca²⁺ and voltage sensors interact, we performed fits of the $Q_{C} (V)$ data using the two different models represented by Scheme I and Scheme II (Figure 4A-B). The allosteric factor $E$ that accounts for the coupling between the Ca²⁺-binding sites and the voltage sensors was constrained to values calculated from the experimental data for the $Q_{C} (V)$ shift at Ca²⁺ saturating conditions (100 µM) in relation to the same curve in the absence of Ca²⁺. The $z_{J}$ , $J_{0}$ and $K_{D}$ parameters obtained during the fitting procedure were very similar for each model (Table 1). The fitted values for the affinity constant ( $K_{D}$ = 6 µM) agree with previous reports (Horrigan and Aldrich, 2002; Cox et al., 1997) although they are slightly smaller than those estimated for the closed conformation of the channel ( $K_{D}$ = 11 µM). Based on the Akaike Information Criterion (AIC) (Akaike, 1974), the model that fits to the $Q_{C} (V, [{C a}^{2 +}])$ data best is Model II, whereas the probability of the Model I being the best model is negligible ( $w_{i}$ = 10⁻⁴⁷) (Table 1). Both models generates a $V_{H}$ - $\log [{C a}^{2 +}]$ curve that accounts reasonably well for the dose-response experimental data (Figure 4C). However, Model I predicts a pronounced decrease in the $z_{Q}$ parameter of the $Q_{C} (V)$ curves at intermediate Ca²⁺ conditions which differs markedly from the experimental values (Figure 4D). This prediction is a consequence of the fractional distribution of two very distinctive functional states of the voltage sensors (unliganded and Ca²⁺-saturated, see above). Additionally, the behavior of $Q_{C} (V)$ curves at intermediate Ca²⁺ concentrations (1-10 µM) is qualitatively consistent with the phenotype exhibited by the Ca²⁺-VSD scheme II (Figure 3D and Figure 4B). Thus, the experimental dose-dependent effect of Ca²⁺ on voltage sensor activation suggests that Ca²⁺-binding to a single α-subunit of BK channels increases $E$ -fold the equilibrium constant $J$ that defines the equilibrium between resting and active conformations of the voltage sensors in all four subunits.

Figure 4 with 1 supplement see all

Download asset Open asset

Dose-dependent effect of Ca²⁺ on voltage sensor activation is predicted by a Ca²⁺-VSD interaction mechanism in which Ca²⁺-binding affects equally the VSD in all four α-subunits.

(**A–B**) The experimental $Q_{C} (V)$ data were fitted using the two possible allosteric interaction mechanisms between voltage and calcium sensors described by Scheme I and Scheme II. The blue and red lines represent the global fits by Model I and Model II, respectively. The allosteric factor $E$ ( $E_{M 1}$ and $E_{M 2}$ ) was constrained to the value obtained from the individual fitting of the $Q_{C} (V) / Q_{C, M A X}$ curves at 0 and 100 µM Ca²⁺ (experimental $E$ ( $E_{e x p}$ ) equal to 26.4, see Table 1). The $z_{J}$ , $J_{0}$ and $K_{D}$ parameters were allowed to vary freely. Note that the allosteric factor $E$ for Model I ( $E_{M 1}$ ) and Model II ( $E_{M 2}$ ) have different interpretations, since $E_{M 1} = E_{e x p}$ whereas $E_{M 2} = \sqrt[4]{E_{e x p}}$ given that the four voltage sensor will be altered in 2.3-fold ( $E_{M 2}$ = 2.27) with each additional Ca²⁺ bound. (**C–D**) The Ca²⁺-dependence of $z_{Q}$ - $Q_{C} (V)$ (C) and $V_{H}$ - $Q_{C} (V)$ (D) curves are superimposed with the $z_{Q}$ and $V_{H}$ values predicted by Model I (blue line) and Model II (red line).

https://doi.org/10.7554/eLife.44934.009

High-affinity Ca²⁺-binding sites in RCK1 and RCK2 domains contribute equally to the allosteric coupling between Ca²⁺ and voltage sensors

Under physiological conditions, the RCK1 and RCK2 high-affinity Ca²⁺-binding sites are responsible for all the calcium sensitivity of the activation of BK channel (Schreiber and Salkoff, 1997; Bao et al., 2002; Bao et al., 2004; Xia et al., 2002). However, distinct physiological roles of the RCK1 Ca²⁺-sensor and Ca²⁺ bowl may be based on their functionally and structurally distinctive properties (Zeng et al., 2005; Sweet and Cox, 2008; Yang et al., 2010; Savalli et al., 2012; Tao et al., 2017). We next asked what the energetic contribution to VSD equilibrium was of the two high-affinity Ca²⁺-binding sites contained in the RCK1 and RCK2 domains.

To elucidate the effect of each Ca²⁺-sensor on the VSD activation we used mutations that selectively and separately abolish the function of the two different RCK Ca²⁺-sites. Disruption of the RCK1 Ca²⁺-sensor by the double mutant D362A/D367A (Xia et al., 2002) significantly reduced (48%, $Δ V_{H}$ _{(D362A/D367A)} = -74.9 ± 4.7 mV) the leftward shift of the $Q_{C} (V)$ curves at 100 μM Ca²⁺ compared with that for the wild-type (WT) BK channel (Figure 5A, C). We also examined the effect of the mutant M513I (Bao et al., 2002) which has been shown to eliminate the Ca²⁺ sensitivity mediated by the RCK1 site (Bao et al., 2002; Bao et al., 2004; Zhang et al., 2010). In this mutant, the 100 µM Ca²⁺-induced shift in $V_{H}$ of the VSD activation curve was also considerably smaller compared to that for the WT channel (about 54%, $Δ V_{H}$ _(M513I) = -65.4 ± 2.6 mV) (Figure 6). Therefore, both mutations affect the Ca²⁺-induced enhancement of the activation of the voltage sensor very similarly through the RCK1 site (Figure 6C), although their mechanisms of action could be quite different. The M513 residue appears to participate in the stabilization of the proper conformation of the RCK1 Ca²⁺-site whereas D367 is a key residue in the coordination of Ca²⁺ ions (Wu et al., 2010; Zhang et al., 2010; Tao et al., 2017). On the other hand, neutralization of the residues forming part of the Ca²⁺ bowl (Schreiber and Salkoff, 1997) (5D5A mutant, see Materials and methods) on the RCK2 domain decreased the leftward shift of the $Q_{C} (V)$ curve by approximately 54% ( $Δ V_{H}$ _(5D5A) = -65.7 ± 4.7 mV) when Ca²⁺ was increased up to 100 μM (Figure 5B, D). The effect of Ca²⁺ binding on ${∆ V}_{H}$ contributed by each high-affinity Ca²⁺ site is roughly half that for WT channels with both sites intact (Figure 5E). Therefore, both high-affinity Ca²⁺-binding sites contribute approximately equally to the decrease in free energy that is necessary to activate the VSD. Indeed, the change of free energy of the resting-active equilibrium of the voltage sensor in response to Ca²⁺-binding at RCK2 site is ~ -4 kJ/mol ( $Δ Δ G_{V}^{C a}$ _{(D362A/D367A)} = -4.2 ± 0.3 kJ/mol and $Δ Δ G_{V}^{C a}$ _(M513I) = -3.6 ± 0.5 kJ/mol) (Figure 5C and Figure 6C). Similarly, the occupation of the RCK1 Ca²⁺-binding site decreased the free energy necessary to activate the VSD in -3.8 ± 0.4 kJ/mol ( $Δ Δ G_{V}^{C a}$ _(5D5A)). Remarkably, these findings reveal an additive effect of Ca²⁺-binding to the RCK1 and Ca²⁺ bowl sites on the VSD activation, which suggest independent allosteric pathways through which they exert their modulation on the VSD.

Figure 5

Download asset Open asset

The high-affinity Ca²⁺-binding sites contribute equally to the allosteric coupling between calcium and voltage sensors in BK channels.

(**A–B**) Representative gating current (*I_G*) recordings at 0 and 100 μM of [Ca²⁺]_i for the RCK1 site mutant (D362A/D367A) and the RCK2 site mutant (5D5A), respectively. (**C–D**) Gating charge-voltage curves ( $Q_{C} (V)$ ) were obtained at 0 Ca²⁺ (open symbols) and 100 μM Ca²⁺ (filled symbols) for D362A/D367A and 5D5A mutants, respectively. Boltzmann fitting to the experimental data (mean ± SEM) is indicated by solid lines ( $V_{H}$ _{(D362A/D367A)}= 178.0 ± 2.7 mV, $z_{Q}$ = 0.58 ± 0.01, $n$ = 12 and $V_{H}$ _(5D5A)= 176.4 ± 4.6 mV, $z_{Q}$ = 0.58 ± 0.01, $n$ = 17 at ‘zero’ Ca²⁺; $V_{H}$ _{(D362A/D367A)}= 104.2 ± 7.3 mV, $z_{Q}$ = 0.56 ± 0.02, $n$ = 7 and $V_{H}$ _(5D5A)= 110.8 ± 6.7 mV, $z_{Q}$ = 0.58 ± 0.02, $n$ = 6 at 100 µM Ca²⁺). For comparison, all $Q_{C} (V)$ plots include the Boltzmann fit of the $Q_{C} (V)$ curves for WT at 0 Ca²⁺ (dashed black line) and 100 μM Ca²⁺ (solid black line). (E) Quantification of the $V_{H}$ shift $(Δ V_{H})$ in the $Q_{C} (V)$ curves and the free energy change $(Δ Δ G_{V}^{C a})$ induced by 100 μM Ca²⁺. The non-parametric t-test was used to evaluate statistical significances between WT BK channel and the RCK sites mutants (***p<0.001).

https://doi.org/10.7554/eLife.44934.011

Figure 6

Download asset Open asset

Mutations abolishing Ca²⁺-sensing by the RCK1 binding-site reduce the Ca²⁺-induced effect on voltage sensors activation similarly.

(A) Representative gating current (I_G) recordings at 0 and 100 μM [Ca²⁺]_i for the RCK1 site mutant M513I. (B) Gating charge-voltage curves $Q_{C} (V)$ were obtained at 0 Ca²⁺ and 100 μM Ca²⁺ (open and filled symbols) for the M513I mutant. Boltzmann fit to the experimental data (mean ± SEM) is indicated by solid lines ( $V_{H}$ _(M513I)= 170.4 ± 4.4 mV, $z_{Q}$ = 0.58 ± 0.01, $n$ = 17 at ‘zero’ Ca²⁺ and $V_{H}$ _(M513I)= 105.0 ± 6.3 mV, $z_{Q}$ = 0.62 ± 0.03, $n$ = 4 at 100 µM Ca²⁺). For comparison, the $Q_{C} (V)$ plot includes the Boltzmann fit of the $Q_{C} (V)$ curves for WT at 0 Ca²⁺ and 100 μM Ca²⁺ (dashed and solid black line). (C) Quantification of the $V_{H}$ shift $(Δ V_{H})$ in the $Q_{C} (V)$ curves and the free energy change induced by 100 μM Ca²⁺ $(Δ Δ G_{V}^{C a})$ . A non-parametric t-test was used to compare WT and RCK1-site mutants BK channels (***p<0.001).

https://doi.org/10.7554/eLife.44934.012

Taking these results into account, we expanded the Ca²⁺-VSD interaction models considering the energetic contribution of the two kinds of Ca²⁺ sensors on the VSD per α-subunit ( $E_{W T} = E_{S 1} * E_{S 2}$ ) (Figure 3—figure supplement 1D). As described in the above model fittings, the allosteric factors $E$ for each one RCK1 and RCK2 sites ( $E_{S 1}$ and $E_{S 2}$ ) were constrained to values equivalent to the Ca²⁺-induced energetic perturbations of the voltage sensor equilibrium for the 5D5A and D362A/D367A mutants, respectively. The inclusion of two Ca²⁺ sites per subunit in the Ca²⁺-VSD interaction schemes increases the functional states of the VSD for both models (Figure 3—figure supplement 1E). Therefore, it becomes difficult to discriminate between different models based on the phenotype of the $Q_{C} (V)$ curves at intermediate Ca²⁺ concentrations (Figure 4—figure supplement 1A-B). However, according to the AIC criteria (compare AICs for the models in Table 1), Model II - Two Ca²⁺ sites is the model that best describes the $Q_{C} (V, [{C a}^{2 +}])$ data ( $w_{i}$ = 0.99). As mentioned when describing the one Ca²⁺ site models, the goodness of the different models can be better appreciated by the Ca²⁺-dependence of the $z_{Q}$ parameter (Figure 4—figure supplement 1C-D). The fractional distribution of the distinct functional states of the VSDs defined by Model I -Two Ca²⁺ sites (unliganded, RCK1 site occupied, RCK2 site occupied, and the two sites occupied; see Figure 3—figure supplement 1E) produce a $z_{Q}$ - $\log [{C a}^{2 +}]$ curve with a pronounced minimum when the Ca²⁺ concentration is close to the $K_{D}$ value (Figure 4—figure supplement 1C). Thus, Model I -Two Ca²⁺ sites predicts a prominent Ca²⁺- dependence of the $z_{Q}$ whereas as we found experimentally, the $z_{Q}$ predicted by Modell II is essentially independent of the Ca²⁺ concentration. The slight apparent increase in $z_{Q}$ data observed in Figure 4—figure supplement 1C as the Ca²⁺ concentration is increased is not statistically significant.

We note here that when the intrasubunit cooperativity factor ( $G$ ) is allowed to vary freely, Model II produces estimates of the $K_{D}$ parameters that are out the range relative to the apparent Ca²⁺ affinities previously reported in the literature for the Ca²⁺-binding sites ( $K_{D (R C K 1)}$ = 13 - 24 µM and $K_{D (R C K 2)}$ = 3 - 5 µM) (Sweet and Cox, 2008; Bao et al., 2002; Xia et al., 2002). In addition, Model II reaches better estimates of the $K_{D}$ parameters (Table 1) if we consider non cooperative interactions between the Ca²⁺-binding sites. However, although the $K_{D}$ value for the RCK2 Ca²⁺-binding site is in agreement with previous reports, Model II underestimates the value for $K_{D (R C K 1)}$ .

Discussion

Recent functional and structural studies have revealed the existence of a major interplay between voltage- and Ca²⁺-sensing modules in the BK channel (Yuan et al., 2010; Savalli et al., 2012; Miranda et al., 2013; Miranda et al., 2016; Miranda et al., 2018; Carrasquel-Ursulaez et al., 2015; Hite et al., 2017; Tao et al., 2017; Zhang et al., 2017), which offer a new perspective on our understanding of its multimodal gating mechanism. However, the details of the CTD-VSD allosteric coupling as well its molecular nature have yet to be firmly established since their direct assessment is experimentally challenging. Based on the functional independence of the distinct structural domains involved (PD, CTD, and VSD), the energetic relationship between the sensory modules can be directly defined by comparing the change in the equilibrium of the voltage sensor under two extreme Ca²⁺ conditions: unliganded and saturated (Horrigan and Aldrich, 2002).

Thus, by characterizing the voltage dependence of charge movement in the virtual absence of internal Ca²⁺ and at Ca²⁺ concentrations that saturate the Ca²⁺ high-affinity sites, this work directly establishes that Ca²⁺-binding significantly facilitates VSD activation through direct energetic contribution to the R-A equilibrium ( $Δ V_{H}$ = -142.6 ± 4.5 mV and $Δ Δ G_{V}^{C a}$ = -7.98 ± 0.27 kJ/mol). Although this result is in agreement with previous reports from our laboratory ( $Δ V_{H}$ = -140 mV and $Δ Δ G_{V}^{C a}$ = -7.9 kJ/mol; at [Ca²⁺]_i = 100 µM) (Carrasquel-Ursulaez et al., 2015), it is at odds with the smaller leftward shift obtained by Horrigan and Aldrich (2002) at saturating Ca²⁺ concentration ( $Δ V_{H}$ = -33 mV and $Δ Δ G_{V}^{C a}$ = -1.9 kJ/mol; at [Ca²⁺]_i = 70 µM). The reason for the contradictory findings is not clear to us. At this stage, we can only identify two differences between the experimental procedures followed by Horrigan and Aldrich (2002) and by us. First, there is a difference in the BK channel clones used (hSlo1 vs mSlo1). Since mBK channels share a high degree of identity with hBKs (96%), it seems unlikely that the differences in Ca²⁺ binding and voltage sensor coupling are due to different BK clones. Second, Horrigan and Aldrich did all their experiments at 0 and 70 μM Ca²⁺ in HEK cells. Here, one can argue that the strength of coupling between Ca²⁺ binding and the voltage sensor may be different in distinct expression systems (e.g., differential modulation) but this is something that is outside the scope of the present paper. Moreover, even if we assume that the calcium effect on VSD is underestimated at 70 µM Ca²⁺ (Horrigan and Aldrich, 2002) compared to 100 µM (our work), we observed a significantly greater effect of Ca²⁺ concentrations (1, 5 and 10 µM, Figure 2) when less than 50% of the Ca²⁺ sensors are occupied ( $K_{D}$ = 11 µM; Horrigan and Aldrich, 2002; Cox et al., 1997).

Fluorescence studies that optically track the motion of the voltage sensor or of the gating ring provide two lines of evidence that support the present findings. First, conformational rearrangements of the voltage sensors detected using voltage-clamp fluorometry can be provoked by Ca²⁺-binding to the high-affinity sites. A sudden rise in intracellular [Ca²⁺] caused by a UV flash induced-photolysis of caged Ca²⁺ prompts a leftward shift in both conductance-voltage ( $G (V)$ ) and fluorescence-voltage ( $F (V)$ ) relationships. These results suggest that functional activation of the gating ring is propagated to the VSD, leading to structural perturbations of voltages sensors, thereby favoring its active conformation (Savalli et al., 2012). Second, the structural rearrangement of the gating ring in response to Ca²⁺ has a voltage dependence (Miranda et al., 2013; Miranda et al., 2018) attributable to the operation of the voltage sensor. The origin of these voltage-dependent motions has recently been established via modifications of the voltage-sensing function of the BK channel using the patch-clamp fluorometry technique (Miranda et al., 2018). Both mutations of the charged residue on the S4 transmembrane segment (R210, R213, and E219) and the co-expression of the β1-subunit with BKα channel, modify the conformational changes of the gating ring triggered by depolarization in correspondence to the observed $G (V)$ shift measured for these channel constructs. In contrast, perturbations of pore opening equilibrium (e.g. through the F315A mutation or the assembly of BKα channel with γ1-subunit) does not modify the voltage-dependent reorganization of the gating ring (Miranda et al., 2018).

Mechanistically, how might the CTD-VSD coupling occur in manner that is independent of channel opening? Taking into account the homotetrameric configuration of the BK channel, Horrigan and Aldrich (2002) defined the general gating scheme of BK channel considering the simplest CTD-VSD interaction model in which voltage sensors and Ca²⁺-binding sites interact solely within the same subunit. However, the VSD movement at non-saturating Ca²⁺ conditions observed here, which entail distinct functional states of the Ca²⁺-binding sites (unliganded and liganded), unveiled that the standard HA model is unable to explain the mechanistic interaction governing the allosteric coupling between the Ca²⁺ and voltage sensors. Assuming that Ca²⁺-sensors are independent and they modify the voltage sensor in the same subunit only, Scheme I would predict $Q_{C} (V)$ curves characteristic of an all-or-none model showing two well-distinguishable Boltzmann components corresponding to the fractions of unliganded and Ca²⁺-saturated BK channels (Figure 4A). Conversely, an energetic effect of each Ca²⁺-site on all the voltage sensors of the tetramer would lead to an equivalent functional status of each VSD, so that the $Q_{C} (V)$ curves would shift as the occupancy of the Ca²⁺ sites increased. Proposing that the VSD and Ca²⁺ sites interact in the manner described by Scheme II reproduces reasonably well the behavior of the Ca²⁺-dependent gating charge movement observed in our experiments (Figure 4B–D). This concerted CTD-VSD communication may underlie a mechanism analogous to the mechanical strategy of interaction between the homooctameric ring of RCK domains and the pore module described for bacterial K⁺ channels (Jiang et al., 2002; Ye et al., 2006; Lingle, 2007; Pau et al., 2011; Smith et al., 2012; Smith et al., 2013). In both MthK and BK channels the Ca²⁺-site occupancy triggers a conformational change corresponding to a symmetric overall rearrangement of the cytosolic tetrameric structure that is ultimately propagated to the transmembrane regions (TMD) via C-linker and in the BK channel, also via the protein-protein interfaces between the gating ring and the TMD (Jiang et al., 2002; Jiang et al., 2003; Ye et al., 2006; Yuan et al., 2010; Yuan et al., 2012; Pau et al., 2011; Smith et al., 2012; Tao et al., 2017). Thus, we can speculate that each Ca²⁺-binding event produces a gradual conformational expansion of the gating ring affecting the four voltage sensors in each step through the progressive perturbations within the protein-protein interfaces.

As mentioned above, the communication pathway through which the Ca²⁺-driven conformational changes are propagated to the voltage sensors appears to reside on the CTD-VSD interface that involves non-covalent interactions between RCK1 N-lobe and S0-S4 transmembrane segments (Yang et al., 2007; Yang et al., 2008; Yang et al., 2010; Sun et al., 2013; Hite et al., 2017; Tao et al., 2017). Although it is not possible at present to dismiss the possibility that the Ca²⁺ binding effect on the VSD workings is, at least in part, mediated by the covalent pulling of the C-linker, we recall that the VSDs are domain-swapped with RCK domains in the gating ring (Tao et al., 2017). Therefore, binding of Ca²⁺ makes the RCK1 N-lobe pull on the S6 helix from its subunit whereas the modification of the contact surface between the RCK1 N-lobe with the voltage sensor of an adjacent subunit induces an outward displacement of the voltage sensor (Hite et al., 2017). These structural arguments make us favor the non-covalent interaction between the CTD and the VSD as the source of the coupling between these two structures. Scanning mutagenesis of RCK1-N terminal subdomain indicates that residues on the βA-αC region are involved in the allosteric connection of the Ca²⁺-dependent activation mediated by RCK1 site occupancy (Yang et al., 2010). In line with this study, the selective activation of the RCK1 domain was identified as being responsible for the Ca²⁺-induced VSD rearrangement (Savalli et al., 2012) and the voltage dependence of the Ca²⁺-driven motions of gating ring (Miranda et al., 2016; Miranda et al., 2018), suggesting that CTD-VSD allosteric coupling is primarily determined by the RCK1 site (Sweet and Cox, 2008). However, our results are inconsistent with this picture. The constructs D362A/D367A and 5D5A (D894A-D898A) selectively impair the Ca²⁺-sensitivity of the RCK1- and RCK2-sensors, respectively, by neutralizing the residues that are involved in contributing to Ca²⁺-coordination (Zhang et al., 2010; Tao et al., 2017). Comparing the fast gating charge movement at 0 Ca²⁺ and saturating Ca²⁺ conditions reveals that the energetic effect of Ca²⁺-binding on voltage sensor equilibrium is practically identical (~ −4 kJ/mol) for both the D362A/D367A and the 5D5A mutations (Figure 5). Thus, our findings establish that the RCK2-driven contribution to CTD-VSD energetic coupling is quite similar to the RCK1-driven contribution. The functional role of the RCK2-sensor on Ca²⁺-sensitivity of VSD activation was further corroborated using the M513I mutation (Figure 6). This point mutation hinders the Ca²⁺-dependent activation associated with the RCK1-sensor, presumably by disrupting the structural integrity of the binding site and the transduction pathway through the βA-αC region (Zhang et al., 2010). Thus, another residue involved in the BK Ca²⁺-dependent activation mediated by the RCK1 Ca²⁺-binding site but not forming part of the site itself, decreases the $Q_{C} (V)$ leftward shift almost in the same amount as does the D362A/D367A mutant.

We found that the energetic contribution of each RCK site to the voltage sensor equilibrium is the same and its combination mimics the VSD Ca²⁺-sensitivity of the fully occupied sites. These findings remind us of early reports showing that mutations in each RCK site shift the Ca²⁺-dependent $G (V)$ by approximately one-half relative to the effect seen for WT channels (Bao et al., 2002; Xia et al., 2002). Thus, our results suggest that the two RCK-sensors contribute independently to the modulation on the VSD, although we cannot eliminate the possibility of some cooperativity between them. Indeed, various lines of evidence indicate that there is some, albeit modest, cooperativity between the two high-affinity Ca²⁺-binding sites although its nature is still unclear (Qian et al., 2006; Savalli et al., 2012; Sweet and Cox, 2008). Intra and intersubunit structural connectivity supports the putative cooperative interactions between the Ca²⁺ sensors at the gating ring (Hite et al., 2017; Yuan et al., 2012). In fact, a recent functional study of the intrasubunit connections between the RCK1 site and Ca²⁺ bowl (R514-Y904/E902 interactions) has shown that such connections are potential candidates for the structural determinants underlying a cooperative mechanism between the RCK1- and RCK2-sensor. These interactions are involved in either the preservation of the integrity of RCK1 Ca²⁺-binding site or define the allosteric propagation pathway of the chemical energy induced by Ca²⁺ binding towards transmembrane domains (Kshatri et al., 2018). On the basis of the cryo-EM structure of Aplysia californica BK channel, (Hite et al., 2017) proposed that there should be a positive cooperativity between the Ca²⁺-binding at RCK1 site and the Ca²⁺ bowl since the Ca²⁺-induced conformational change of the RCK1-N lobes from closed to open configuration depends on the functional state (unliganded and liganded) of both RCK sites.

Our analysis based on the CTD-VSD interaction model can not specify likely cooperative relations among the two high-affinity Ca²⁺ sites within the same α-subunit. In fact, the analysis of our gating current data when using cooperativity between the two high-affinity sites (Table 1) produced a result at odds with previous reports. Based on the structural information of the BK channels, it could be considered that cooperative interactions between the Ca²⁺ sensors of the different α-subunits (Hite et al., 2017) can also account, in part, for the Ca²⁺-dependent behavior of the VSDs. However, functional studies point out a more relevant role of the intrasubunit cooperativity (albeit modest) than the intersubunit cooperativity between the RCK Ca²⁺ sites (Niu and Magleby, 2002; Qian et al., 2006). Thus, although a concerted CTD-VSD model (Scheme II) gives better explanation than an independent CTD-VSD model (Scheme I) to the allosteric communication of the calcium and voltage sensors, more work is required to explore improved models able to reproduce more accurately the properties of the interaction CTD-VSD mechanism in BK channel.

Table 1

Parameters for the best fits of the $Q_{C} (V)$ data using different Ca²⁺-VSD interaction models.

https://doi.org/10.7554/eLife.44934.013

One Ca²⁺-site per α-subunit			Two Ca²⁺-sites per α-subunit
	Model I	Model II		Model I		Model II
Parameters	Model I	Model II	Parameters	With cooperativity	Without cooperativity	With cooperativity	Without cooperativity
*z_J*	0.61	0.61	z_J	0.61	0.61	0.61	0.61
J₀	0.020	0.018	J₀	0.019	0.019	0.021	0.019
E	26.4^*	2.27^*	E_S₁	4.57^*		1.46^*
E	26.4^*	2.27^*	E_S₂	5.35^*		1.52^*
K_D (µM)	6.4	6.1	K_D₁ (µM)	3.2	4.9	837.7	5.9
K_D (µM)	6.4	6.1	K_D₂ (µM)	631.7	6.9	6.6	5.9
			G	56.1	1^*	120.4	1^*
AIC	−948.4	−1150.9		−1088.9	−1090.3	−1162.0	−1147.5
L_i	4^*10⁻⁴⁷	0.004		1^*10⁻¹⁶	3^*10⁻¹⁶	1	7^*10⁻⁴
w_i	4^*10⁻⁴⁷	0.004		1^*10⁻¹⁶	3^*10⁻¹⁶	0.995	0.001

^*Fixed parameters in the model fitting. AIC values correspond to Akaike Information Criterion to select the best fit model. ℒ_i and w_i are the relative likelihood and the weight of each model within the set of candidate models.

In conclusion, our results depict a remarkable, and direct energetic interplay between the specialized sensory modules (VSD and CTD) of the BK channel. Our findings together with the emerging structural-functional information establish a new paradigm about how stimuli integration (depolarization and intracellular Ca²⁺) modulates this channel’s activation and its relevance within a physiological context. Notable and unexpected is the equivalent contribution of the distinct ligand-binding sites in the cytosolic domain to the allosteric regulation of voltage sensing. Additional studies to discern the molecular bases underlying the Ca²⁺ and voltage propagation pathways and the cooperative interactions of the RCK1 and RCK2 regulatory domains may provide new clues about the dual gating mechanism of BK channel.

Materials and methods

Channel expression

Request a detailed protocol

Xenopus laevis oocytes were used as a heterologous system to express BK channels. The cDNA coding for the human BK α-subunit (U11058) was provided by L. Toro (University of California, Los Angeles, CA). The cDNA coding for independent mutants of each two high-affinity Ca²⁺ site from BK channel, the double mutant D362A/D367A (Xia et al., 2002), the mutant M513I (Bao et al., 2002) in the RCK1 Ca²⁺-binding site, and the mutant 5D5A (Schreiber and Salkoff, 1997) (D894A/D895A/D896A/D897A/D898A) in the RCK2 Ca²⁺-binding site or calcium bowl, were kindly provided by M. Holmgren (National Institutes of Health, Bethesda, MD). The cRNA was prepared using mMESSAGE mMACHINE (Ambion) for in vitro transcription. Xenopus laevis oocytes were injected with 50 ng of cRNA and incubated in an ND96 solution (in mM: 96 NaCl, 2 KCl, 1.8 CaCl₂, 1 MgCl₂, 5 HEPES, pH 7.4) at 18°C for 4–8 days before electrophysiological recordings.

Electrophysiological recordings

Request a detailed protocol

All recordings were made by using the patch-clamp technique in the inside-out configuration. Data were acquired with an Axopatch 200B (Molecular Devices) amplifier and the Clampex 10 (Molecular Devices) acquisition software. Gating current (I_G) records were elicited by 1 ms voltage steps from −90 to 350 mV in increments of 10 mV. Both the voltage command and current output were filtered at 20 kHz using an 8-pole Bessel low-pass filter (Frequency Devices). Current signals were sampled with a 16-bit A/D converter (Digidata 1550B; Molecular Devices), using a sampling rate of 500 kHz. Linear membrane capacitance and leak subtraction were performed based on a P/4 protocol (Armstrong and Bezanilla, 1974).

Borosilicate capillary glasses (1B150F-4, World Precision Instruments) were pulled in a horizontal pipette puller (Sutter Instruments). After fire-polishing, pipette resistance was 0.5–1 MΩ. The external (pipette) solution contained (in mM): 110 tetraethylammonium (TEA)-MeSO₃, 10 HEPES, 2 MgCl₂; pH was adjusted to 7.0. The internal solution (bath) contained (in mM): N-methyl-D-glucamine (NMDG)-MeSO₃, 10 HEPES, and 5 EGTA for ‘zero Ca²⁺’ solution (∼0.8 nM, based on the presence of ∼10 μM contaminant [Ca²⁺] (Cui et al., 1997). An agar bridge containing 1 M NaMES connected the internal solution to a pool of the external solution grounded with an Ag/AgCl electrode. The calculated bridge/bath junction potential was ~0.8 mV. For test solutions at different Ca²⁺ concentrations (0.1–100 μM), CaCl₂ was added to reach the desired free [Ca²⁺], and 5 mM EGTA (0.1–0.5 μM) or HEDTA (1–10 μM) was used as calcium buffer. No Ca²⁺ chelator was used in 100 µM Ca²⁺ solutions. Free calcium concentration was estimated using the WinMaxChelator Software and checked with a Ca²⁺-electrode (Hanna Instruments). All experiments were performed at room temperature (20–22°C). To measure I_G at different Ca²⁺ concentrations in the same oocyte, the patch was excised and washed with an appropriate internal solution using at least 10 times the chamber volume.

Data analysis

Request a detailed protocol

All data analysis was performed using Clampfit 10 (Molecular Devices, RRID:SCR_011323), Matlab (MathWorks, RRID:SCR_001622) and Excel 2007 (Microsoft, RRID:SCR_016137). The first 50-100 µs of the ON-gating currents were fitted to a single exponential function and the area under the curve (Figure 1A-B) was integrated to obtain the charge displaced between closed states ( $Q_{C}$ ) (Horrigan and Aldrich, 1999; Horrigan and Aldrich, 2002; Carrasquel-Ursulaez et al., 2015; Contreras et al., 2012). $Q_{C} (V)$ data for each [Ca²⁺]_i were fitted using a Boltzmann function:

Q_{C} (V) = \frac{Q_{C, M A X}}{1 + e^{(\frac{- z_{Q} F (V - V_{H})}{R T})}}

where $Q_{C, M A X}$ is the maximum charge, $z_{Q}$ is the voltage dependence of activation, $V_{H}$ is the half-activation voltage, $T$ is the absolute temperature (typically 295 K), $F$ is the Faraday’s constant, and $R$ is the universal gas constant. $Q_{C, M A X}$ , $V_{H}$ , and $z_{Q}$ were determined using least square minimization. $Q_{C} (V)$ curves were aligned by shifting them along the voltage axis by the mean $Δ V = (⟨ V_{H} ⟩ - V_{H})$ to generate a mean curve that did not alter the voltage dependence (Horrigan et al., 1999). All error estimates are SEM.

For each experiment, the time constants obtained from exponential fits to ON-gating currents were shifted along the voltage axis by $Δ V$ to determine the mean $τ_{I_{G} - O N} (V)$ relationships. $τ_{I_{G} - O N} (V)$ data were fitted to a two-states process described by

τ (V) = 1 / (α (V) + β (V))

where $α (V) = α_{0} e^{(z_{J} δ F V / R T)}$ and $β (V) = β_{0} e^{(z_{J} (δ - 1) F V / R T)}$ are, respectively, the forward and backward rate constants which determine the equilibrium constant $J$ of the voltage sensor (see below). The parameter $δ$ is the electrical distance at which the peak of the energy barrier that separates the two resting (R)-active (A) states of the voltage sensor is located, and $z_{J}$ is the number of gating charges displaced during the R-A transition.

The Ca²⁺-induced effect on VSD activation was quantified as the $V_{H}$ shift relative to 'zero' Ca²⁺ condition: $Δ V_{H} = V_{H ({[C a^{2 +}]}_{i})} - V_{H (0 {[C a^{2 +}]}_{i})}$ . For wild-type (WT) BK channel and the RCK Ca²⁺-sensor mutants (D362A/D367A, M513I and 5D5A), the energetic contribution of Ca²⁺-binding on resting-active (R-A) equilibrium of the voltage sensor was calculated as changes in Gibbs free energy of VSD activation induced by 100 µM Ca²⁺:

Δ Δ G_{V}^{C a} = F (z_{Q (100 μ M {[C a^{2 +}]}_{i})} V_{H (100 μ M {[C a^{2 +}]}_{i})} - z_{Q (0 {[C a^{2 +}]}_{i})} V_{H (0 {[C a^{2 +}]}_{i})})

Model fitting

Request a detailed protocol

We fit the $Q_{C} (V, [{C a}^{2 +}])$ experimental data using two distinct interaction mechanisms between Ca²⁺-binding sites and voltage sensor (see Scheme I and Scheme II in the Figure 3A-B) within the framework of Horrigan-Aldrich (HA) general allosteric model (Horrigan and Aldrich, 2002). Assumptions and considerations for the equations that describe each one of the Ca²⁺-VSD interaction schemes are given in the Appendix. In terms of the HA allosteric mechanisms, the voltage sensor R-A equilibrium is defined by the equilibrium constant $J$ according to the relation:

J = e^{\frac{z_{J} F (V - V_{H})}{R T}} = J_{0} e^{\frac{z_{J} F V}{R T}}

where $J_{0}$ is the zero voltage equilibrium constant and $z_{J}$ the gating charge displacement per voltage sensor. Thus, the fraction of the total charge displaced essentially between closed states, $(Q_{C} (V) / Q_{C, M A X})$ in the absence of calcium can be written as:

Q_{C} (V) / Q_{C, M A X} (C a^{2 +} ≪ K_{D}) = \frac{1}{1 + J^{- 1}}

where $K_{D}$ is the dissociation constant of the high-affinity calcium-binding site with all voltage sensors at rest and the channel closed. In the presence of saturating Ca²⁺ (100 µM), the equilibrium of the R-A transition $J$ becomes amplified by the allosteric factor $E$ , which defines the coupling between Ca²⁺-binding sites and voltage sensors, being

Q_{C} (V) / Q_{C, M A X} (C a^{2 +} ≫ K_{D}) = \frac{1}{1 + {(J E)}^{- 1}}

and

J E = J_{0} e^{\frac{R T \ln (E) + z_{J} F V}{R T}}

The $Q_{C} (V) / Q_{C, M A X}$ measured in the presence of high [Ca²⁺] and 'zero Ca²⁺' condition at the same voltage (so that $J$ be canceled out) but in the limit where $({J E)}^{- 1} ≫ 1$ is:

{\frac{Q_{C} (V) / Q_{C, M A X} (C a^{2 +} ≫ K_{D})}{Q_{C} (V) / Q_{C, M A X} (C a^{2 +} ≪ K_{D})}}_{(lim J E^{- 1} ≫ 1)} = \frac{J E}{J} = E

Thus, the Gibbs free energy perturbation of the voltage sensor R-A equilibrium when the high-affinity binding sites are approximately 100% occupied by Ca²⁺ (100 µM) is a straightforward measure of the allosteric factor: $E = e^{- Δ Δ G_{V}^{C a} / R T}$ .

Based on these conditions, the values of the allosteric parameter $E$ were calculated and introduced in each of the two Ca²⁺-VSD interaction models as a fixed parameter. Once $E$ was obtained, the families of $Q_{C} (V, [{C a}^{2 +}])$ curves were simultaneously fitted to Equation 4 and Equation 6 (see Appendix) and estimating the $z_{J}$ , $J_{0}$ and $K_{D}$ parameters for each model by minimizing least-square values. To select the best Ca²⁺-VSD interaction scheme that describes the experimental data, the fits provided by each model were compared according to their Akaike Information Criterion ( $A I C$ ) values (Akaike, 1974), calculated as ${A I C}_{i} = 2 p_{i} - 2 \ln (L_{i})$ , where $p_{i}$ is the number of free parameters and $\ln (L_{i})$ is the maximum log-likelihood of the model $i$ . The best fit being the one that achieves the lowest ${A I C}_{i}$ value. Minimum ${A I C}_{i}$ ( ${A I C}_{M I N}$ ) values were used as model selection criteria. Using the ${A I C}_{i}$ weights ( $w_{i}$ ), we estimated the probability of model $i$ is the best model given the data and the set of candidate models. $w_{i}$ are based on the relative likelihood of each tested model $i$ which is a function of the difference in ${A I C}_{i}$ score and the best model: $Δ A I C_{i} = A I C_{i} - A I C_{M I N}$ (Anderson and Burnham, 2002). From the $Δ A I C_{i}$ we obtained an estimate of the relative likelihood of model $i$ ( $ℒ_{i}$ ) by the simple transform: $L_{i} = e x p (- \frac{1}{2} Δ A I C_{i})$ . $w_{i}$ is calculated normalizing the $ℒ_{i}$ for each model: $w_{i} = ℒ_{i} / \sum_{k = 1}^{K} ℒ_{k}$ , where $K$ is the number of candidate models.

The models for the Ca²⁺-VSD interaction schemes were extended including two high-affinity Ca²⁺-binding sites per α-subunit (Figure 3—figure supplement 1D-E). The contribution of each Ca²⁺-binding site to the free energy of the voltage sensor equilibrium may be split in two, such as $E = E_{S 1} * E_{S 2} = e^{- (Δ Δ G_{V}^{C a} (S 1) + Δ Δ G_{V}^{C a} (S 2)) / R T}$ , where $E_{S 1}$ and $E_{S 2}$ are the allosteric factors $E$ for the RCK1 and RCK2 sites. Thus, for the global fit of the $Q_{C} (V, [{C a}^{2 +}])$ curves, we constrained the allosteric parameter $E_{S 1}$ and $E_{S 2}$ obtained experimentally for the RCK2 Ca²⁺-sensor mutant (5D5A) and RCK1 Ca²⁺-sensor mutant (D362A/D367A), respectively, as described above. The rest of the parameters $z_{J}$ , $J_{0}$ , $K_{D 1}$ , $K_{D 2}$ , and $G$ , where $K_{D 1}$ and $K_{D 2}$ are the dissociation constants of the RCK1 and RCK2 sites and $G$ is a cooperativity factor between the two sites within the same α-subunit of the BK channel, were allowed to vary freely.

Appendix 1

Assumptions and model testing

We assumed that the four voltage sensors act independently, transiting between two states, resting (R) and active (A), governed by the voltage-dependent equilibrium constant $J$ . The R–A equilibrium is displaced toward the active state by membrane depolarization, generating a fast gating charge movement ( $Q_{C}$ ) before channel opening. Additionally, the Ca²⁺-binding to high-affinity sites shifts the voltage sensor equilibrium toward its active configuration through an allosteric coupling described by the factor $E$ (Figure 3—figure supplement 1A). By assuming the simplified standard model for the BK channels (Horrigan and Aldrich, 2002), where each α-subunit has a single Ca²⁺-binding site, we established the possible states and their connections through which each voltage sensor transits in the presence of Ca²⁺ (Figure 3—figure supplement 1B–C) following the CTD-VSD interaction mechanisms described by Scheme I and II (Figure 3A–B).

For Scheme I, in which Ca²⁺-binding sites and voltage sensors interact within the same α-subunit, the activation of each VSD can occur through the R₀-A₀ or R₁-A₁ transitions according to the functional state of the Ca²⁺ site (unbound or Ca²⁺ bound). The equilibrium of such transitions is governed by $J$ and $J E_{M 1}$ , respectively (Figure 3—figure supplement 1B). In the case of Scheme II, in which binding of Ca²⁺ to a single α-subunit affects all four voltage sensors equally, the R-A equilibrium for each VSD would be affected by the number of Ca²⁺ bound in the channel (0-4) depicted in the model (Model II) as five possible R-A transitions. According to this model, the $J$ constant increase $E_{M 2}$ -fold for each Ca²⁺ site occupied (Figure 3—figure supplement 1C). For both schemes, the horizontal transitions R-R and A-A represent the Ca²⁺-binding equilibrium ( $K$ or $K E$ ) when the VSD is in the resting and active conformation, respectively. The equilibrium constant $K$ is defined as the bound/unbound probability ratio for each Ca²⁺-binding site and depends on Ca²⁺ concentration ([Ca²⁺]) and the Ca²⁺ dissociation constant ( $K_{D}$ ): $K = [{C a}^{2 +}] / K_{D}$ .

Here, we assumed that the voltage sensor movement at ON-gating currents is in equilibrium relative with the binding of Ca²⁺. The assumption is reasonable since the Ca²⁺-binding rate constant estimated for BK channel is 1.8 x 10⁸ M⁻¹s⁻¹ (Hou et al., 2016) implying that at 10 µM internal Ca²⁺ the association time constant is 340 µs. Thus, Ca²⁺ binding at this Ca²⁺ concentration proceeds with a time constant about six-fold longer than does the movement of the voltage sensor (~60 μs). Based on this consideration, the R-A transitions in the models would be predominant transitions whose proportion would be determined by the [Ca²⁺] and $K_{D}$ . Therefore, predictions of the $Q_{C} (V)$ curves at different Ca²⁺ concentrations for Model I and Model II were based on a given fractional occupancy of Ca²⁺ sites established by the probability of Ca²⁺ bound $(b)$ and unbound $(1 - b)$ for each Ca²⁺-sensor, and the energetic contribution to VSD equilibrium.

Simulations of the $Q_{C} (V)$ curves using the Scheme I (Model I) were obtained using the equation

\frac{Q_{C} (V)}{Q_{C, M A X}} = (1 - b) (\frac{1}{1 + J^{- 1}}) + b (\frac{1}{1 + {(J E_{M 1})}^{- 1}});

where

b = \frac{1}{1 + K^{- 1}} = \frac{1}{1 + \frac{K_{D}}{[C a^{2 +}]}} = \frac{[C a^{2 +}]}{[C a^{2 +}] + K_{D}};

and

J = J_{0} e^{\frac{z_{J} F V}{R T}}

Substituting $b$ and $J$ into Equation 1, the Ca²⁺-dependent voltage sensor activation for Model I is given by the equation

\frac{Q_{C} (V)}{Q_{C, M A X}} = (\frac{K_{D}}{[C a^{2 +}] + K_{D}}) (\frac{1}{1 + \frac{e^{\frac{- z_{J} F V}{R T}}}{J_{0}}}) + (\frac{[C a^{2 +}]}{[C a^{2 +}] + K_{D}}) (\frac{1}{1 + \frac{e^{\frac{- z_{J} F V}{R T}}}{J_{0} E_{M 1}}})

Thus, the $Q_{C} (V)$ curves are determined by the proportion of two functional VSD populations with a distinctive effect (unliganded effect or Ca²⁺-saturated effect). Consequently, the $Q_{C} (V)$ curves are represented by a weighted sum of two Boltzmann functions.

Meanwhile, for the concerted CTD-VSD interaction Scheme II (Model II), the $Q_{C} (V, [{C a}^{2 +}])$ curves would be determined using the general equation:

\frac{Q_{C} (V)}{Q_{C, M A X}} = \sum_{x = 0}^{n} (\begin{matrix} n \\ x \end{matrix}) {(1 - b)}^{n - x} b^{x} (\frac{1}{1 + {(J E_{M 2}^{x})}^{- 1}})

The expression in the first bracket represents the fraction of VSD belonging to a channel with $x$ (0 to 4) Ca²⁺ bound, according to a binomial probability distribution. Thus, the $Q_{C} (V)$ curves result in a weighted sum of five distinct Boltzmann functions corresponding to the five possible R-A transitions (Figure 3—figure supplement 1C). By setting $n =$ 4 because the tetrameric symmetry of the channels, and substituting $b$ and $J$ into the previous equation (Equation 5) we have

\frac{Q_{C} (V)}{Q_{C, M A X}} = \sum_{x = 0}^{4} (\begin{matrix} 4 \\ x \end{matrix}) {(\frac{K_{D}}{[C a^{2 +}] + K_{D}})}^{4 - x} {(\frac{[C a^{2 +}]}{[C a^{2 +}] + K_{D}})}^{x} (\frac{1}{1 + \frac{e^{\frac{- z_{J} F V}{R T}}}{J_{0} E_{M 2}^{x}}})

It should be noted that at limiting Ca²⁺ conditions, both schemes become equivalent where the VSD activation is characterized by a single Boltzmann function. At zero Ca²⁺, the $Q_{C} (V)$ curves are described by

\frac{Q_{C} (V)}{Q_{C, M A X}} = (\frac{1}{1 + \frac{e^{\frac{- z_{J} F V}{R T}}}{J_{0}}}),

whereas Ca²⁺ saturating concentration $J$ is multiplied by the allosteric factor $E$ , where $E = E_{M 1} = E_{M 2}^{4}$ depending on the model (Model I or Model II):

\frac{Q_{C} (V)}{Q_{C, M A X}} = (\frac{1}{1 + \frac{e^{\frac{- z_{J} F V}{R T}}}{J_{0} E}})

Given that each α-subunit has two Ca²⁺-binding sites, we expanded the CTD-VSD interaction schemes (Figure 3—figure supplement 1C) considering the existence of two Ca²⁺-binding sites (Figure 3—figure supplement 1D-E). The Model I and II include the energetic contribution of RCK1 and RCK2 Ca²⁺-sites to the VSD activation. The factor $E = E_{S 1} {* E}_{S 2}$ where $E_{S 1}$ and $E_{S 1}$ are the allosteric coupling between the VSD and the RCK1 Ca²⁺-site and RCK2 Ca²⁺-site, respectively. The $K_{1}$ and $K_{2}$ constants define the bound/unbound transition for each RCK1 and RCK2 sites with $K_{1} = [{C a}^{2 +}] / K_{D 1}$ and $K_{2} = [{C a}^{2 +}] / K_{D 2}$ . Assuming that the Ca²⁺ sensors of distinct α-subunit do not interact, we only consider intrasubunit cooperativity between the RCK1 and RCK2 sites defined by the factor $G$ . Thus, the occupancy of one RCK site will affect Ca²⁺-binding equilibrium to the other RCK site in the α-subunit ( ${G K}_{1}$ and ${G K}_{2}$ ) (Figure 3—figure supplement 1E). Note that for Model II the equilibrium $J$ of the VSD increase $E_{S 1}$ -fold and $E_{S 2}$ -fold for each Ca²⁺ bound to the RCK1 and RCK2 sites, respectively, reaching $J E_{S 1}^{4} E_{S 2}^{4}$ when the eight Ca²⁺ sites are occupied.

Data availability

All data generated or analysed during this study are included in the manuscript and supporting files.

References

1. Akaike H
(1974) A new look at the statistical model identification
IEEE Transactions on Automatic Control 19:716–723.
https://doi.org/10.1109/TAC.1974.1100705
- Google Scholar
Book
1. Anderson DR
2. Burnham K
(2002) Model Selection and Multi-Model Inference: A Practical Information-Theoretic Approach
New York: Springer-Verlag.
https://doi.org/10.1007/b97636
- Google Scholar
1. Armstrong CM
2. Bezanilla F
(1974) Charge movement associated with the opening and closing of the activation gates of the na channels
The Journal of General Physiology 63:533–552.
https://doi.org/10.1085/jgp.63.5.533
- PubMed
- Google Scholar
1. Bao L
2. Rapin AM
3. Holmstrand EC
4. Cox DH
(2002) Elimination of the BK(Ca) channel's high-affinity Ca(2+) sensitivity
The Journal of General Physiology 120:173–189.
https://doi.org/10.1085/jgp.20028627
- PubMed
- Google Scholar
1. Bao L
2. Kaldany C
3. Holmstrand EC
4. Cox DH
(2004) Mapping the BKCa channel's "Ca2+ bowl": side-chains essential for Ca2+ sensing
The Journal of General Physiology 123:475–489.
https://doi.org/10.1085/jgp.200409052
- PubMed
- Google Scholar
(2015) Hydrophobic interaction between contiguous residues in the S6 transmembrane segment acts as a stimuli integration node in the BK channel
The Journal of General Physiology 145:61–74.
https://doi.org/10.1085/jgp.201411194
- PubMed
- Google Scholar
(2012) Modulation of BK channel voltage gating by different auxiliary β subunits
PNAS 109:18991–18996.
https://doi.org/10.1073/pnas.1216953109
- PubMed
- Google Scholar
1. Cox DH
2. Cui J
3. Aldrich RW
(1997) Allosteric gating of a large conductance Ca-activated K+ channel
The Journal of General Physiology 110:257–281.
https://doi.org/10.1085/jgp.110.3.257
- PubMed
- Google Scholar
1. Cui J
2. Cox DH
3. Aldrich RW
(1997) Intrinsic voltage dependence and Ca2+ regulation of mslo large conductance Ca-activated K+ channels
The Journal of General Physiology 109:647–673.
https://doi.org/10.1085/jgp.109.5.647
- PubMed
- Google Scholar
1. Cui J
2. Aldrich RW
(2000) Allosteric linkage between voltage and ca(2+)-dependent activation of BK-type mslo1 K(+) channels
Biochemistry 39:15612–15619.
https://doi.org/10.1021/bi001509+
- PubMed
- Google Scholar
(2017) Structural basis for gating the high-conductance Ca²⁺-activated K⁺ channel
Nature 541:52–57.
https://doi.org/10.1038/nature20775
- PubMed
- Google Scholar
(1999) Allosteric voltage gating of potassium channels I. mslo ionic currents in the absence of ca(2+)
The Journal of General Physiology 114:277–304.
https://doi.org/10.1085/jgp.114.2.277
- PubMed
- Google Scholar
1. Horrigan FT
2. Aldrich RW
(1999) Allosteric voltage gating of potassium channels II. mslo channel gating charge movement in the absence of ca(2+)
The Journal of General Physiology 114:305–336.
https://doi.org/10.1085/jgp.114.2.305
- PubMed
- Google Scholar
1. Horrigan FT
2. Aldrich RW
(2002) Coupling between voltage sensor activation, Ca2+ binding and channel opening in large conductance (BK) potassium channels
The Journal of General Physiology 120:267–305.
https://doi.org/10.1085/jgp.20028605
- PubMed
- Google Scholar
1. Hou P
2. Xiao F
3. Liu H
4. Yuchi M
5. Zhang G
6. Wu Y
7. Wang W
8. Zeng W
9. Ding M
10. Cui J
11. Wu Z
12. Wang LY
13. Ding J
(2016) Extrapolating microdomain ca(2+) dynamics using BK channels as a ca(2+) sensor
Scientific Reports 6:17343.
https://doi.org/10.1038/srep17343
- PubMed
- Google Scholar
1. Jiang Y
2. Lee A
3. Chen J
4. Cadene M
5. Chait BT
6. MacKinnon R
(2002) Crystal structure and mechanism of a calcium-gated potassium channel
Nature 417:515–522.
https://doi.org/10.1038/417515a
- PubMed
- Google Scholar
1. Jiang Y
2. Ruta V
3. Chen J
4. Lee A
5. MacKinnon R
(2003) The principle of gating charge movement in a voltage-dependent K+ channel
Nature 423:42–48.
https://doi.org/10.1038/nature01581
- PubMed
- Google Scholar
(2018) Functional validation of Ca2+-binding residues from the crystal structure of the BK ion channel
Biochimica Et Biophysica Acta (BBA) - Biomembranes 1860:943–952.
https://doi.org/10.1016/j.bbamem.2017.09.023
- Google Scholar
(1982) Reconstitution in planar lipid bilayers of a Ca2+-dependent K+ channel from transverse tubule membranes isolated from rabbit skeletal muscle
PNAS 79:805–809.
https://doi.org/10.1073/pnas.79.3.805
- PubMed
- Google Scholar
(2017) Molecular determinants of BK channel functional diversity and functioning
Physiological Reviews 97:39–87.
https://doi.org/10.1152/physrev.00001.2016
- PubMed
- Google Scholar
1. Lingle CJ
(2007) Gating rings formed by RCK domains: keys to gate opening
The Journal of General Physiology 129:101–107.
https://doi.org/10.1085/jgp.200709739
- PubMed
- Google Scholar
1. Marty A
(1981) Ca-dependent K channels with large unitary conductance in chromaffin cell membranes
Nature 291:497–500.
https://doi.org/10.1038/291497a0
- PubMed
- Google Scholar
(2013) State-dependent FRET reports calcium- and voltage-dependent gating-ring motions in BK channels
PNAS 110:5217–5222.
https://doi.org/10.1073/pnas.1219611110
- PubMed
- Google Scholar
(2016) Interactions of divalent cations with calcium binding sites of BK channels reveal independent motions within the gating ring
PNAS 113:14055–14060.
https://doi.org/10.1073/pnas.1611415113
- PubMed
- Google Scholar
(2018) Voltage-dependent dynamics of the BK channel cytosolic gating ring are coupled to the membrane-embedded voltage sensor
eLife 7:e40664.
https://doi.org/10.7554/eLife.40664
- PubMed
- Google Scholar
1. Niu X
2. Magleby KL
(2002) Stepwise contribution of each subunit to the cooperative activation of BK channels by Ca2+
PNAS 99:11441–11446.
https://doi.org/10.1073/pnas.172254699
- PubMed
- Google Scholar
(1988) Activation by divalent cations of a Ca2+-activated K+ channel from skeletal muscle membrane
The Journal of General Physiology 92:67–86.
https://doi.org/10.1085/jgp.92.1.67
- PubMed
- Google Scholar
(1981) Single channel recordings of Ca2+-activated K+ currents in rat muscle cell culture
Nature 293:471–474.
https://doi.org/10.1038/293471a0
- PubMed
- Google Scholar
(2011) Structure and function of multiple Ca2+-binding sites in a K+ channel regulator of K+ conductance (RCK) domain
PNAS 108:17684–17689.
https://doi.org/10.1073/pnas.1107229108
- PubMed
- Google Scholar
1. Qian X
2. Niu X
3. Magleby KL
(2006) Intra- and intersubunit cooperativity in activation of BK channels by ^Ca2+
The Journal of General Physiology 128:389–404.
https://doi.org/10.1085/jgp.200609486
- PubMed
- Google Scholar
1. Rothberg BS
2. Magleby KL
(1999) Gating kinetics of single large-conductance Ca2+-activated K+ channels in high Ca2+ suggest a two-tiered allosteric gating mechanism
The Journal of General Physiology 114:93–124.
https://doi.org/10.1085/jgp.114.1.93
- PubMed
- Google Scholar
1. Rothberg BS
2. Magleby KL
(2000) Voltage and Ca2+ activation of single large-conductance Ca2+-activated K+ channels described by a two-tiered allosteric gating mechanism
The Journal of General Physiology 116:75–100.
https://doi.org/10.1085/jgp.116.1.75
- PubMed
- Google Scholar
1. Savalli N
2. Pantazis A
3. Yusifov T
4. Sigg D
5. Olcese R
(2012) The contribution of RCK domains to human BK channel allosteric activation
Journal of Biological Chemistry 287:21741–21750.
https://doi.org/10.1074/jbc.M112.346171
- PubMed
- Google Scholar
1. Schreiber M
2. Salkoff L
(1997) A novel calcium-sensing domain in the BK channel
Biophysical Journal 73:1355–1363.
https://doi.org/10.1016/S0006-3495(97)78168-2
- PubMed
- Google Scholar
1. Shen K-Z
2. Lagrutta A
3. Davies NW
4. Standen NB
5. Adelman JP
6. North RA
(1994) Tetraethylammonium block of slowpoke calcium-activated potassium channels expressed in xenopus oocytes: evidence for tetrameric channel formation
Pflügers Archiv European Journal of Physiology 426:440–445.
https://doi.org/10.1007/BF00388308
- Google Scholar
(2012) Crystal structure of a ba(2+)-bound gating ring reveals elementary steps in RCK domain activation
Structure 20:2038–2047.
https://doi.org/10.1016/j.str.2012.09.014
- PubMed
- Google Scholar
(2013) Structural basis of allosteric interactions among Ca2+-binding sites in a K+ channel RCK domain
Nature Communications 4:1–10.
https://doi.org/10.1038/ncomms3621
- Google Scholar
1. Stefani E
2. Ottolia M
3. Noceti F
4. Olcese R
5. Wallner M
6. Latorre R
7. Toro L
(1997) Voltage-controlled gating in a large conductance Ca2+-sensitive K+channel (hslo)
PNAS 94:5427–5431.
https://doi.org/10.1073/pnas.94.10.5427
- PubMed
- Google Scholar
1. Sun X
2. Shi J
3. Delaloye K
4. Yang X
5. Yang H
6. Zhang G
7. Cui J
(2013) The interface between membrane-spanning and cytosolic domains in Ca²+-dependent K+ channels is involved in β subunit modulation of gating
Journal of Neuroscience 33:11253–11261.
https://doi.org/10.1523/JNEUROSCI.0620-13.2013
- PubMed
- Google Scholar
1. Sweet TB
2. Cox DH
(2008) Measurements of the BKCa channel's high-affinity ^Ca2+binding constants: effects of membrane voltage
The Journal of General Physiology 132:491–505.
https://doi.org/10.1085/jgp.200810094
- PubMed
- Google Scholar
(2017) Cryo-EM structure of the open high-conductance Ca2+-activated K+ channel
Nature 541:46–51.
https://doi.org/10.1038/nature20608
- Google Scholar
1. Wu Y
2. Yang Y
3. Ye S
4. Jiang Y
(2010) Structure of the gating ring from the human large-conductance ca(2+)-gated K(+) channel
Nature 466:393–397.
https://doi.org/10.1038/nature09252
- PubMed
- Google Scholar
1. Xia XM
2. Zeng X
3. Lingle CJ
(2002) Multiple regulatory sites in large-conductance calcium-activated potassium channels
Nature 418:880–884.
https://doi.org/10.1038/nature00956
- PubMed
- Google Scholar
1. Yang H
2. Hu L
3. Shi J
4. Delaloye K
5. Horrigan FT
6. Cui J
(2007) Mg2+ mediates interaction between the voltage sensor and cytosolic domain to activate BK channels
PNAS 104:18270–18275.
https://doi.org/10.1073/pnas.0705873104
- PubMed
- Google Scholar
1. Yang H
2. Shi J
3. Zhang G
4. Yang J
5. Delaloye K
6. Cui J
(2008) Activation of Slo1 BK channels by Mg2+ coordinated between the voltage sensor and RCK1 domains
Nature Structural & Molecular Biology 15:1152–1159.
https://doi.org/10.1038/nsmb.1507
- PubMed
- Google Scholar
1. Yang J
2. Krishnamoorthy G
3. Saxena A
4. Zhang G
5. Shi J
6. Yang H
7. Delaloye K
8. Sept D
9. Cui J
(2010) An epilepsy/dyskinesia-associated mutation enhances BK channel activation by potentiating Ca2+ sensing
Neuron 66:871–883.
https://doi.org/10.1016/j.neuron.2010.05.009
- PubMed
- Google Scholar
1. Yang H
2. Zhang G
3. Cui J
(2015) BK channels: multiple sensors, one activation gate
Frontiers in Physiology 6:1–16.
https://doi.org/10.3389/fphys.2015.00029
- Google Scholar
1. Ye S
2. Li Y
3. Chen L
4. Jiang Y
(2006) Crystal structures of a ligand-free MthK gating ring: insights into the ligand gating mechanism of K+ channels
Cell 126:1161–1173.
https://doi.org/10.1016/j.cell.2006.08.029
- PubMed
- Google Scholar
1. Yuan P
2. Leonetti MD
3. Pico AR
4. Hsiung Y
5. MacKinnon R
(2010) Structure of the human BK channel Ca2+-activation apparatus at 3.0 A resolution
Science 329:182–186.
https://doi.org/10.1126/science.1190414
- PubMed
- Google Scholar
(2012) Open structure of the Ca2+ gating ring in the high-conductance Ca2+-activated K+ channel
Nature 481:94–97.
https://doi.org/10.1038/nature10670
- Google Scholar
1. Zeng XH
2. Xia XM
3. Lingle CJ
(2005) Divalent cation sensitivity of BK channel activation supports the existence of three distinct binding sites
The Journal of General Physiology 125:273–286.
https://doi.org/10.1085/jgp.200409239
- PubMed
- Google Scholar
1. Zhang G
2. Huang SY
3. Yang J
4. Shi J
5. Yang X
6. Moller A
7. Zou X
8. Cui J
(2010) Ion sensing in the RCK1 domain of BK channels
PNAS 107:18700–18705.
https://doi.org/10.1073/pnas.1010124107
- PubMed
- Google Scholar
1. Zhang G
2. Geng Y
3. Jin Y
4. Shi J
5. McFarland K
6. Magleby KL
7. Salkoff L
8. Cui J
(2017) Deletion of cytosolic gating ring decreases gate and voltage sensor coupling in BK channels
The Journal of General Physiology 149:373–387.
https://doi.org/10.1085/jgp.201611646
- PubMed
- Google Scholar
1. Zhou Y
2. Zeng XH
3. Lingle CJ
(2012) Barium ions selectively activate BK channels via the Ca2+-bowl site
PNAS 109:11413–11418.
https://doi.org/10.1073/pnas.1204444109
- PubMed
- Google Scholar
1. Zhou Y
2. Yang H
3. Cui J
4. Lingle CJ
(2017) Threading the biophysics of mammalian Slo1 channels onto structures of an invertebrate Slo1 channel
The Journal of General Physiology 149:985–1007.
https://doi.org/10.1085/jgp.201711845
- PubMed
- Google Scholar

Article and author information

Author details

Yenisleidy Lorenzo-Ceballos
1. Doctorado en Ciencias Mención Neurociencia, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
2. Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
Contribution
Conceptualization, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing—original draft

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4309-9314
Willy Carrasquel-Ursulaez

Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile

Contribution
Conceptualization, Software, Formal analysis, Investigation, Writing—review and editing

Competing interests
No competing interests declared
Karen Castillo

Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile

Contribution
Validation, Investigation, Methodology, Writing—review and editing

Competing interests
No competing interests declared
Osvaldo Alvarez
1. Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
2. Departamento de Biología, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
Contribution
Conceptualization, Software, Formal analysis, Writing—review and editing

Competing interests
No competing interests declared
Ramon Latorre

Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile

Contribution
Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing

For correspondence
ramon.latorre@uv.cl

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-6044-5795

Funding

AFOSR (No. FA9550-16-1-0384)

Ramon Latorre

FONDECYT (Grant No. 1180999)

Karen Castillo

CONICYT-PFCHA (Doctoral fellowship No. 63140149)

Yenisleidy Lorenzo-Ceballos

Chilean Ministry of Economy, Development, and Tourism (Millennium Scientific Initiative P029-022-F)

Yenisleidy Lorenzo-Ceballos
Willy Carrasquel-Ursulaez
Karen Castillo
Osvaldo Alvarez
Ramon Latorre

FONDECYT (Grant No. 1150273)

Ramon Latorre

FONDECYT (Grant No. 1190203)

Ramon Latorre

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dr. John Ewer for his comments on the manuscript and to Mrs. Luisa Soto (University of Valparaiso) for excellent technical assistance. This research was supported by FONDECYT Grant No. 1150273, FONDECYT Grant No. 1190203 and AFOSR No. FA9550-16-1-0384 to RL; CONICYT-PFCHA Doctoral fellowships No. 63140149 to YLC; FONDECYT Grant No. 1180999 to KC The Centro Interdisciplinario de Neurociencia de Valparaiso is a Millennium Institute supported by the Millennium Scientific Initiative of the Chilean Ministry of Economy, Development, and Tourism (P029-022-F).

Ethics

Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Ethics Committee for Animal Experimentation of the University of Valparaíso. All of the animals were handled according to approved institutional animal care and use committee protocols (BEA031-14) of the University of Valparaiso. All surgery was performed under tricaine anesthesia, and every effort was made to minimize suffering.

Version history

Received: January 7, 2019
Accepted: September 10, 2019
Accepted Manuscript published: September 11, 2019 (version 1)
Version of Record published: September 26, 2019 (version 2)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,555

views
321

downloads
11

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Yenisleidy Lorenzo-Ceballos
Willy Carrasquel-Ursulaez
Karen Castillo
Osvaldo Alvarez
Ramon Latorre

(2019)

Calcium-driven regulation of voltage-sensing domains in BK channels

eLife 8:e44934.

https://doi.org/10.7554/eLife.44934

Share this article

Cite this article

Ca2+ binding strongly affects the activation of VSD in BK channels.

Ca2+-dependent effects on VSD activation in BK channels.

Model-dependent behavior of the QC(V) curves based on the CTD-VSD interaction mechanisms according to the fractional occupancy of Ca2+-binding sites.

Dose-dependent effect of Ca2+ on voltage sensor activation is predicted by a Ca2+-VSD interaction mechanism in which Ca2+-binding affects equally the VSD in all four α-subunits.

The high-affinity Ca2+-binding sites contribute equally to the allosteric coupling between calcium and voltage sensors in BK channels.

Mutations abolishing Ca2+-sensing by the RCK1 binding-site reduce the Ca2+-induced effect on voltage sensors activation similarly.

Parameters for the best fits of the QC(V) data using different Ca2+-VSD interaction models.

Author details

Yenisleidy Lorenzo-Ceballos

Contribution

Competing interests

Willy Carrasquel-Ursulaez

Contribution

Competing interests

Karen Castillo

Contribution

Competing interests

Osvaldo Alvarez

Contribution

Competing interests

Ramon Latorre

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organisms

Further reading

Ca²⁺ binding strongly affects the activation of VSD in BK channels.

Ca²⁺-dependent effects on VSD activation in BK channels.

Model-dependent behavior of the $Q_{C} (V)$ curves based on the CTD-VSD interaction mechanisms according to the fractional occupancy of Ca²⁺-binding sites.

Dose-dependent effect of Ca²⁺ on voltage sensor activation is predicted by a Ca²⁺-VSD interaction mechanism in which Ca²⁺-binding affects equally the VSD in all four α-subunits.

The high-affinity Ca²⁺-binding sites contribute equally to the allosteric coupling between calcium and voltage sensors in BK channels.

Mutations abolishing Ca²⁺-sensing by the RCK1 binding-site reduce the Ca²⁺-induced effect on voltage sensors activation similarly.

Parameters for the best fits of the $Q_{C} (V)$ data using different Ca²⁺-VSD interaction models.