Voltage-dependent dynamics of the BK channel cytosolic gating ring are coupled to the membrane-embedded voltage sensor
Abstract
In humans, large conductance voltage- and calcium-dependent potassium (BK) channels are regulated allosterically by transmembrane voltage and intracellular Ca2+. Divalent cation binding sites reside within the gating ring formed by two Regulator of Conductance of Potassium (RCK) domains per subunit. Using patch-clamp fluorometry, we show that Ca2+ binding to the RCK1 domain triggers gating ring rearrangements that depend on transmembrane voltage. Because the gating ring is outside the electric field, this voltage sensitivity must originate from coupling to the voltage-dependent channel opening, the voltage sensor or both. Here we demonstrate that alterations of the voltage sensor, either by mutagenesis or regulation by auxiliary subunits, are paralleled by changes in the voltage dependence of the gating ring movements, whereas modifications of the relative open probability are not. These results strongly suggest that conformational changes of RCK1 domains are specifically coupled to the voltage sensor function during allosteric modulation of BK channels.
https://doi.org/10.7554/eLife.40664.001Introduction
The open probability of large conductance voltage-and Ca2+-activated K+ (BK or slo1) channels is regulated allosterically by voltage and intracellular concentration of divalent ions (Barrett et al., 1982; Moczydlowski and Latorre, 1983; Horrigan and Aldrich, 2002; Latorre et al., 2017). This feature makes BK channels important regulators of physiological processes such as neurotransmission and muscular function, where they couple membrane voltage and the intracellular concentration of Ca2+ (Robitaille and Charlton, 1992; Hu et al., 2001; Wang et al., 2001; Raffaelli et al., 2004). The BK channel is formed in the membrane as tetramers of α subunits, encoded by the KCNMA1 gene (Shen et al., 1994; Quirk and Reinhart, 2001). Each α subunit contains seven transmembrane domains (S0 to S6), a small extracellular N-terminal domain and a large intracellular C-terminal domain (Wallner et al., 1996; Meera et al., 1997; Tao et al., 2017) (Figure 2a). Similar to other voltage-gated channels, the voltage across the membrane is sensed by the voltage sensor domain (VSD), containing charged amino acids within transmembrane segments S2, S3 and S4 (Díaz et al., 1998; Ma et al., 2006; Pantazis and Olcese, 2012; Tao et al., 2017). The sensor for divalent cations is at the C-terminal region and is formed by two Regulator of Conductance for K+ domains (RCK1 and RCK2) per α subunit (Wei et al., 1994; Moss and Magleby, 2001; Xia et al., 2002; Zeng et al., 2005; Wu et al., 2010). In the tetramer, four RCK1-RCK2 tandems pack against each other in a large structure known as the gating ring (Wu et al., 2010; Yuan et al., 2011; Giraldez and Rothberg, 2017; Tao et al., 2017; Zhou et al., 2017). Two high-affinity Ca2+ binding sites are located in the RCK2 (also known as ‘Ca2+ bowl’) and RCK1 domains, respectively. Additionally, a site with low affinity for Mg2+ and Ca2+ is located at the interface between the VSD and the RCK1 domain (Shi and Cui, 2001; Zhang et al., 2001; Bao et al., 2002; Xia et al., 2002; Yang et al., 2007; Yang et al., 2008a; Tao et al., 2017) (Figure 2a). The high-affinity binding sites show structural dissimilarity (Zhang et al., 2010; Tao et al., 2017) and different affinity for divalent ions (Zeng et al., 2005). Apart from Ca2+, it has been described that Cd2+ selectively binds to the RCK1 site, whereas Ba2+ and Mg2+ show higher affinity for the RCK2 site (Xia et al., 2002; Zeng et al., 2005; Yang et al., 2008b; Zhou et al., 2012; Miranda et al., 2016). Thus, intracellular concentrations of Ca2+, Cd2+, Ba2+ or Mg2+ can shift the voltage dependence of BK activation towards more negative potentials. Using patch clamp fluorometry (PCF), we have shown that these cations trigger independent conformational changes of RCK1 and/or RCK2 within the gating ring, measured as large changes in the efficiency of Fluorescence Resonance Energy Transfer (FRET) between fluorophores introduced into specific sites in the BK tetramer. These rearrangements depend on the specific interaction of the divalent ions with their high-affinity binding sites, showing different dependences on cation concentration and membrane voltage (Miranda et al., 2013; Miranda et al., 2016). To date, the proposed transduction mechanism by which divalent ion binding increases channel open probability was a conformational change of the gating ring that leads to a physical pulling of the channel gate, where the linker between the S6 transmembrane domain and the RCK1 region acts like a passive spring (Niu et al., 2004). Such a mechanism would be analogous to channel activation by ligand binding in glutamate receptor or cyclic nucleotide-gated ion channels, also tetramers (Sobolevsky et al., 2009; James et al., 2017). Our previous results do not support this as the sole mechanism underlying coupling of divalent ion binding to channel opening, since the gating ring conformational changes that we have recorded: 1) are not strictly coupled to the opening of the channel’s gate, and 2) show different voltage dependence for each divalent ion. In addition, the recent cryo-EM structure of the full slo1 channel of Aplysia californica (Hite et al., 2017; Tao et al., 2017) shows that the RCK1 domain of the gating ring is in contact with the VSD, predicting that changes in the voltage sensor position could be reflected in the voltage dependent gating ring reorganizations.
Understanding the nature of the voltage dependence associated with individual rearrangements produced by binding of divalent ions to the gating ring is essential to untangle the mechanism underlying the role of such rearrangements in BK channel gating. To this end, we have now performed PCF measurements with human BK channels heterologously expressed in Xenopus oocytes, including a range of VSD mutations or co-expressed with different regulatory subunits. Here we provide evidence for a functional interaction between the gating ring and the voltage sensor in full-length, functional BK channels at the plasma membrane, in agreement with the structural data from Aplysia BK. Moreover, these data support a pathway that couples to divalent ion binding to channel opening through the voltage sensor.
Results
Voltage dependence of gating ring rearrangements is associated to activation of the RCK1 binding site
BK α subunits labeled with fluorescent proteins CFP and YFP in the linker between the RCK1 and RCK2 domains (position 667) retain the functional properties of wild-type BK channels (Miranda et al., 2013; Miranda et al., 2016). This allowed us to use PCF to detect conformational rearrangements of the gating ring measured as changes in FRET efficiency (E) between the fluorophores (Miranda et al., 2013; Miranda et al., 2016). Binding of Ca2+ ions to both high-affinity binding sites (RCK1 and Ca2+ bowl) produces an activation of BK channels, coincident with an increase in E from basal levels reaching saturating values at high Ca2+ concentrations (Miranda et al., 2013 and Figure 1a). In addition, we observed that the E signal has the remarkable property that in intermediate Ca2+ concentrations (from 4 µM to 55 µM) it shows voltage dependence besides its Ca2+ dependence (Miranda et al., 2013 and Figure 1a). As discussed previously (Miranda et al., 2013), these changes in E with voltage are not conformational dynamics of the gating ring that simply follow the voltage dependence of VSD. For instance, at 0 Ca2+ concentrations movements of the VSD occurs between 0 and +300 mV (Stefani et al., 1997; Horrigan et al., 1999; Horrigan and Aldrich, 2002; Zhang et al., 2014; Carrasquel-Ursulaez et al., 2015; Zhang et al., 2017). However, we do not observe changes in E between 0 and +240 mV (Figure 1a). Similarly, at 100 µM Ca2+, charge movement takes place between −100 and +150 mV (Carrasquel-Ursulaez et al., 2015), while our FRET signals at 95 µM Ca2+ do not vary within this voltage range (Figure 1a). Independent activation of high-affinity binding sites by other divalent ions (Ba2+, Cd2+, or Mg2+ (Miranda et al., 2016)) led us to postulate that Ca2+ activation has a site-dependent relation to voltage. To further evaluate the effect of individual high-affinity Ca2+ binding sites on the voltage-dependent component of the gating ring conformational changes we first selectively mutated the binding sites. Mutations D362A and D367A (Xia et al., 2002; Zeng et al., 2005) were introduced in the BK667CY construct (BK667CYD362A/D367A) to remove the high-affinity binding site located in the RCK1 domain. Figure 1b shows the relative conductance and E values for the BK667CYD362A/D367A construct at different membrane voltages for various Ca2+ concentrations. As described previously, the G-V curves show a significantly reduced shift to more negative potentials when Ca2+ is increased, as compared to the non-mutated BK667CY (Figure 1a–b, left panels). Specific activation of the Ca2+ bowl renders a smaller change in E values, which are not voltage-dependent within the voltage range tested (Figure 1b, right panel). To test the effect of eliminating the RCK2 Ca2+ binding site -the Ca2+ bowl- we mutated five aspartates to alanines (5D5A) (Bao et al., 2002). As expected, activation of only the RCK1 domain by Ca2+ reduced the Ca2+-dependent shift in the GV curves (Figure 1c, left panel). Even though the extent to which the E values changed with Ca2+ was reduced (Figure 1c), there was a persistent voltage dependence equivalent to that shown in Figure 1a corresponding to the non-mutated channel (most appreciable at 12 µM and 22 µM Ca2+ concentrations; Figure 1c, right panel) (Miranda et al., 2013). Further, at these two Ca2+ concentrations the changes in E occurred within the same voltage range (+60–120 mV) in channels with the Ca2+ bowl mutated (Figure 1c) or not (Figure 1a). This effect seems not to be attributable to Ca2+ binding to unknown binding sites in the channel, since the double mutation of the RCK1 and RCK2 sites abolishes the change in the FRET signal (Figure 1d). Altogether, these results indicate that the voltage-dependent component of the gating ring conformational changes triggered by Ca2+ in the BK667CY construct depends on activation of the RCK1 binding site. Because the gating ring is not within the transmembrane region, it is not expected to be directly influenced by the transmembrane voltage. Therefore, the voltage-dependent FRET signals must be coupled to the dynamics of the gate region associated with the opening and closing of the channel and/or those of the voltage sensor domain.
![](https://iiif.elifesciences.org/lax/40664%2Felife-40664-fig1-v2.tif/full/617,/0/default.jpg)
Voltage dependence of gating ring rearrangements is associated to activation of the RCK1 binding site.
G-V (left panels) and E-V curves (right panels) obtained simultaneously at several Ca2+concentrations from (a) the BK667CY construct, (b) mutation of the RCK1 high-affinity site (D362A/D367A), (c) mutation of the Ca2+ bowl (5D5A), or (d) both (D362A/D367A 5D5A). Note that the voltage dependence of the E signal is only abolished after mutating the RCK1 high-affinity binding site (b) or both (d). Data corresponding to each Ca2+ concentration are color-coded as indicated in the legend at the bottom. Solid curves in the G-V graphs represent Boltzmann fits. For reference, grey shadows in (a–d) left panels represent the full range of G-V curves corresponding to non-mutated BK667CY channels from 0 μM Ca2+ to 95 µM Ca2+ (indicated with colored dashed lines). Data points and error bars represent average ± SEM (n = 3–14, N = 2–8). Part of the data in (a, b and c) are taken from (Miranda et al., 2013) and (Miranda et al., 2016).
The voltage-dependent conformational changes of the gating ring are not related to the opening and closing of the pore domain
To test whether the voltage-dependent FRET signals relate to the opening and closing of the channel (intrinsic gating) we used two modifications of BK channel function in which the relative probability of opening is shifted in the voltage axis, yet the actual dynamics of voltage sensor are expected to be unaltered (Figure 2b). We reasoned that, if the voltage-dependent FRET signals of the gating ring are coupled to the opening and closing, they should follow a similar displacement with voltage. The first BK channel construct is the α subunit including the single point mutation F315A, which has been described to shift the voltage dependence of the relative conductance of the channel to more positive potentials, by uncoupling the voltage sensor activation from the gate opening (Figure 2c) (Carrasquel-Ursulaez et al., 2015). Figure 2d shows the relative conductance and E vs. voltage for the BK667CYF315A mutant at various Ca2+concentrations. Our results show that the shift of the relative probability of opening to more positive potentials (Figure 2d, left panel) does not lead to changes in the voltage dependence of the gating ring FRET signals (Figure 2d, right panel).
![](https://iiif.elifesciences.org/lax/40664%2Felife-40664-fig2-v2.tif/full/617,/0/default.jpg)
Modification of the voltage dependence of gate opening does not affect the gating ring voltage-dependent conformational changes.
(a) Topology of the BKα subunit where the voltage sensing domain (VSD), Ca2+ sensing domain (gating ring, GR) and pore domain (PD) are indicated by colored dashed lines boxes (see main text for a full description). (b) The three BK functional modules (VSD, PD, GR), schematically represented as colored boxes, interact allosterically. (c) Diagram representing the main effect of the F315A mutation, which is the uncoupling of the VSD to the PD. (d) G-V (left panel) and E-V curves (right panel) obtained simultaneously at several Ca2+concentrations after mutation of the F315 site to alanine (BK667CYF315A). It should be noted that the extent of the shifts induced by the mutation are smaller than previously reported (Carrasquel-Ursulaez et al., 2015), which could arise from the different experimental conditions and/or our fluorescent construct. (e) The interaction with the γ1 subunit favors the VSD-PD coupling mechanism (f) G-V (left) and E-V curves (right) of BK667CY α subunits co-expressed with γ1 subunits. In all panels, data corresponding to each Ca2+ concentration are color-coded as indicated in the bottom legend. Colored dashed lines represent the G-V and E-V curves corresponding to BK667CYα channels (Miranda et al., 2013; Miranda et al., 2016). The solid curves in the G-V graphs represent Boltzmann fits. The full range of G-V curves from 0 μM Ca2+ to 95 µM Ca2+ from BK667CY is represented as a grey shadow in left panels (d and f), for reference. Data points and error bars represent average ± SEM (n = 3–8; N = 2–3).
The second modification of BK function consisted in co-expressing the wild type α subunit with the auxiliary subunit γ1 (Yan and Aldrich, 2010; Yan and Aldrich, 2012; Gonzalez-Perez et al., 2014; Li and Yan, 2016). In this case, the relative probability of opening is shifted to more negative potentials by increasing the coupling between the voltage sensor and the gate of the channel (Figure 2e). This construct adds the advantage of representing a physiologically relevant modification of channel gating. Figure 2f shows the relative conductance and E vs. voltage in oocytes co-expressing the BK667CYα and γ1 at voltages ranging from −160 to +260 mV, with three [Ca2+] concentrations: nominal 0, 12 µM and 22 µM. As expected, the presence of the γ1 subunit drives the relative conductance curves to more negative potentials (Figure 2f, left panel) compared to the values obtained without γ1 (Figure 2f, dashed lines). Remarkably, the change in the voltage dependence of the relative conductance induced by γ1 does not alter the simultaneously recorded FRET signals (Figure 2f, right panel), which remains indistinguishable from that recorded with BK667CYα (Figure 2f, dashed lines).
The dynamics of the VSD are directly reflected in the gating ring conformation
Using the allosteric HA model of BK channel function, Horrigan and Aldrich (2002) proposed that Ca2+ binding to the Ca2+ bowl is coupled to the voltage sensor activation. Yet, the strength of that interaction (allosteric constant E) was smaller than those corresponding to Ca2+- or V-sensors with channel opening (Horrigan and Aldrich, 2002). Interestingly, when E was derived from gating currents data, a larger value was obtained (Carrasquel-Ursulaez et al., 2015). Further, Ca2+ binding to the RCK1 domain (but not to the Ca2+ bowl) is voltage-dependent (Sweet and Cox, 2008), which as the authors hypothesized might originate from physical interactions between the voltage sensors and the RCK1 domains. Additionally, using the cut-open oocyte voltage-clamp fluorometry approach, Savalli et al. (2012) showed that fluorescence emission from reporters within the VSD could change upon uncaged Ca2+ stimuli. This evidence indicates that the VSD is coupled to the gating ring, but none of these approaches directly monitored the conformational changes of the gating ring structure. Therefore, we decided to explore whether the voltage dependence of the gating ring movements is attributable to the voltage sensor activation. To this end we modified the voltage dependence of the VSD activation by co-expression with β auxiliary subunits or by introducing specific mutations in the VSD (Figure 3 and Figure 4). The effects of co-expressing BK α subunit with the four different types of auxiliary β subunits have been extensively studied (Tseng-Crank et al., 1996; Behrens et al., 2000; Brenner et al., 2000; Cox and Aldrich, 2000; Uebele et al., 2000; Lingle et al., 2001; Zeng et al., 2001; Bao and Cox, 2005; Orio and Latorre, 2005; Yang et al., 2008a; Sweet and Cox, 2009; Contreras et al., 2012; Li and Yan, 2016). β1 subunit has been previously proposed to alter the voltage sensor-related voltage dependence, as well as the intrinsic opening of the gate and Ca2+ sensitivity (Figure 3a) (Cox and Aldrich, 2000; Bao and Cox, 2005; Orio and Latorre, 2005; Sweet and Cox, 2009; Contreras et al., 2012; Castillo et al., 2015). Recordings from BK667CYα co-expressed with β1 subunits reveal the expected modifications in the voltage dependence of the relative conductance, that is an increase in the apparent Ca2+ sensitivity (Figure 3b, left panel) (Wallner et al., 1995; Cox and Aldrich, 2000; Bao and Cox, 2005; Orio and Latorre, 2005; Sweet and Cox, 2009; Contreras et al., 2012). In addition, it has been reported that β1 subunit alters the function of the VSD (Orio and Latorre, 2005; Castillo et al., 2015). Notably, the E-V curves are shifted to more negative potentials (Figure 3b, right panel), similarly to the described modification (Castillo et al., 2015). The structural determinants of the β1 subunit influence on the VSD reside within its N-terminus, which has been shown by engineering a chimera between the β3b subunit (which does not influence the VSD) and the N-terminus of the β1 (β3bNβ1) (Castillo et al., 2015). We recapitulated this strategy. First, we co-expressed BK667CY α subunits with β3b and observed the expected inactivation of the ionic currents at positive potentials, yet with different blockade kinetics (see Figure 3—figure supplement 1) (Uebele et al., 2000; Xia et al., 2000; Lingle et al., 2001). The relative open probability of this complex is like BK667CYα alone, except that at extreme positive potentials the values of relative conductance at the tails decrease due to inactivation (Figure 3—figure supplement 1b, left panel). The values of E vs V remained comparable to those observed for BK667CYα (Figure 3—figure supplement 1b, right panel). We then co-expressed the β3bNβ1 chimera (Castillo et al., 2015) with BK667CYα (Figure 3c). This complex did not modify the relative conductance vs. voltage relationship (Figure 3d, left panel) as compared with BK667CYα alone (Figure 3d, grey shadow). On the other hand, while the magnitude of the FRET change is the same as in BK667CYα, the voltage dependence of E values at [Ca2+] of 4 μM, 12 μM and 22 μM shifted to more negative potentials compared to the values of BK667CYα alone (Figure 3d, right panel, compare dashed to solid lines). Altogether, these results indicate that the alteration of the voltage dependence of the voltage sensor induced by the amino terminal of β1within the β3bNβ1 chimera underlies the modification of the voltage dependence of the gating ring conformational changes, reinforcing the hypothesis that this voltage dependence is directly related to VSD function.
![](https://iiif.elifesciences.org/lax/40664%2Felife-40664-fig3-v2.tif/full/617,/0/default.jpg)
Co-expression with β subunits.
(a) β1 subunits have been shown to directly regulate VSD function, shifting Vh(j) to more negative values (b) Left panel, G-V curves obtained at several Ca2+concentrations after co-expression of BK667CY with the β1 subunit, which induces a leftward shift in the E-V curves obtained simultaneously (right). (c) β3bNβ1 chimeras produce similar effects to β1 on VSD function, since they retain the N-terminal region of β1 (Castillo et al., 2015). (d) G-V (left) and E-V curves (right) of BK667CY α subunits co-expressed with the β3bNβ1 chimera. Data corresponding to each Ca2+ concentration are color-coded as indicated in the legend at the bottom. Colored dashed lines represent the G-V and E-V curves corresponding to BK667CYα channels (Miranda et al., 2013; Miranda et al., 2016). The solid curves in the G-V graphs represent Boltzmann fits. The full range of G-V curves from 0 μM Ca2+ to 95 µM Ca2+ from BK667CY is represented as a grey shadow in left panels (b and d), for reference. Data points and error bars represent average ± SEM (n = 3–10; N = 2–4).
![](https://iiif.elifesciences.org/lax/40664%2Felife-40664-fig4-v2.tif/full/617,/0/default.jpg)
Mutation of charged residues of BK VSD.
VSD activation was altered by mutation of charged residues in the VSD that modify its voltage of half activation, Vh(j) (a) The R210E mutation induces a negative shift of Vh(j)(b) G-V (left panel) and E-V curves (right panel) obtained simultaneously from constructs BK667CY containing the R210E mutation at several Ca2+concentrations. (c) The E219R mutation produces a negative shift of Vh(j) (d) G-V (left panel) and E-V curves (right panel) obtained simultaneously from constructs BK667CY containing the E219R mutation at several Ca2+concentrations. (e) The R213E mutation induces a large positive shift of Vh(j) values. (f) G-V (left panel) and E-V curves (right panel) obtained simultaneously from constructs BK667CY containing the R213E mutation at several Ca2+concentrations. Data corresponding to each Ca2+ concentration are color-coded as indicated in the bottom legend. Colored dashed lines represent the G-V and E-V curves corresponding to non-mutated BK667CYα channels (Miranda et al., 2013; Miranda et al., 2016). The solid curves in the G-V graphs represent Boltzmann fits. The full range of G-V curves from 0 μM Ca2+ to 95 µM Ca2+ from BK667CY is represented as a grey shadow in left panels (b), (d and f), for reference. Data points and error bars represent average ± SEM (n = 4–10; N = 3–4).
VSD activation can also be altered by introducing single point mutations that modify the voltage of half activation of the voltage sensor, Vh(j). This parameter is determined by fitting data to the HA allosteric model (Ma et al., 2006) or directly from gating current measurements (Zhang et al., 2014). Mutations of charged amino acids on the VSD have been reported to produce different modifications in the Vh(j) values. In some cases, other parameters related to BK channel activation are additionally affected by the mutations. Mutation R210E shifts the Vh(j) value from +173 mV to +25 mV at 0 Ca2+ in BK channels (Figure 4a) (Ma et al., 2006). Consistent with this, introduction of this mutation in BK667CYα (BK667CYR210E) caused a shift of the relative conductance vs. voltage dependence towards more negative potentials (Figure 4b, left panel) as compared to BK667CY (Figure 4b, left panel, grey shadow). Simultaneously measured E values showed a negative shift in the voltage dependence of the FRET signal at intermediate Ca2+ concentrations (Figure 4b, right panel). Mutation E219R had been previously shown to produce a large negative shift in Vh(j) from +150 mV to +40 mV (ΔVh(j) = −110 mV; Figure 4c), additionally modifying the Ca2+ sensitivity and the coupling between the VSD and channel gate (Zhang et al., 2014). As previously reported, BK667CYE219R showed modified relative conductance vs. voltage relationships at different Ca2+ concentrations (Figure 4d, left panel) (Zhang et al., 2014). In addition, this construct revealed a shift to more negative potentials in the E vs. voltage dependence at intermediate Ca2+ concentrations (12 μM and 22 μM Ca2+; Figure 4d, right panel), paralleling the reported negative shift in Vh(j) (Ma et al., 2006; Zhang et al., 2014). Since mutations displacing the Vh(j) to more negative potentials induce equivalent shifts in the voltage dependence of the gating ring motion (measured as E), we tested if other mutations previously reported to induce positive shifts on Vh(j) (Ma et al., 2006) were also associated with changes of the E-V curves in the same direction. As shown by Ma et al., the largest effect on Vh(j) is induced by the R213E mutation, producing a shift of ΔVh(j)=+337 mV (Figure 4e) (Ma et al., 2006). The BK667CYR213E construct showed a significant shift in the voltage dependence of the relative conductance to more positive potentials (Figure 4f, left panel). Notably, this effect was paralleled by a large displacement in the E vs. voltage dependence towards more positive potentials (Figure 4f, right panel). Taken together, our data show that modifications of the Vh(j) values caused by mutating the VSD charged residues are reflected in equivalent changes in the voltage dependence of the gating ring conformational rearrangements, which occur in analogous directions and with proportional magnitudes at intermediate Ca2+ concentrations.
All these results on the VSD modifications and their corresponding changes in FRET signals support the existence of a direct coupling mechanism between the VSD function and the gating ring conformational changes.
Parallel alterations of the voltage dependence of VSD function and gating ring motions by selective activation of the RCK1 binding site
We have previously shown that specific interaction of Cd2+ with the RCK1 binding site leads to activation of the BK channel, which is accompanied by voltage-dependent changes in the E values at intermediate Cd2+ concentrations of 10 μM and 30 μM (Miranda et al., 2016). To further assess the role of the RCK1 binding site activation in the voltage dependence of the gating ring motions, we studied activation by Cd2+ of selected BK667CY VSD mutants (Figure 5). Addition of Cd2+ to the BK667CYE219R mutant (Figure 5a) shifted the voltage dependence of E towards more negative potentials at intermediate Cd2+ concentrations (10 μM and 30 μM; Figure 5b) when compared to non-mutated BK667CY (Figure 5b; dashed lines). This change in the E-V curves induced by selective activation of the RCK1 binding site with Cd2+ paralleled the large negative shift (ΔVh(j) = −110 mV) previously reported with the E219R mutant BK channels (Ma et al., 2006; Zhang et al., 2014). We also tested Cd2+ activation in the mutant BK667CYR201Q, which shifts the Vh(j) parameter by 47 mV towards positive potentials (Figure 5c) (Ma et al., 2006). Addition of Cd2+ rendered right-shifted E vs. voltage relationships (Figure 5d, right panel), following the direction of the predicted Vh(j) shift described for this mutant BK channel (Ma et al., 2006). Finally, addition of Cd2+ to the BK667CYF315A construct (Figure 5e) (Carrasquel-Ursulaez et al., 2015) did not have any effect on the E-V relationship (Figure 5f). These results are consistent with a mechanism in which specific binding of Cd2+ to the RCK1 binding site allows voltage-dependent conformational changes in the gating ring that are directly related to VSD activation.
![](https://iiif.elifesciences.org/lax/40664%2Felife-40664-fig5-v2.tif/full/617,/0/default.jpg)
Voltage dependence of gating ring rearrangements after specific activation of RCK1 high-affinity binding site by Cd2+.
(a) Effect of the VSD E219R mutation on the selective activation of RCK1 by Cd2+. (b) G-V (left panels) and E-V curves (right panels) obtained simultaneously at several Ca2+concentrations from constructs BK667CYE219R. (c) VSD R201Q mutation induces a positive shift of Vh(j) (d) G-V (left panels) and E-V curves (right panels) obtained simultaneously at several Cd2+concentrations from constructs BK667CYR201Q (e) Effect of the F315A mutation on the selective activation of RCK1 by Cd2+. (f) G-V (left panels) and E-V curves (right panels) obtained simultaneously at several Cd2+ concentrations from constructs BK667CYF315A. Data corresponding to each Cd2+ concentration are color-coded as indicated in the legend at the bottom. Colored dashed lines represent the G-V and E-V curves corresponding to BK667CYα channels (Miranda et al., 2013; Miranda et al., 2016). The solid curves in the G-V graphs represent Boltzmann fits. The full range of G-V curves from 0 μM Cd2+ to 100 µM Cd2+ corresponding to non-mutated BK667CY is represented as a grey shadow in left panels (b), (d), and (f), for reference. Data points and error bars represent average ± SEM (n = 3–4; N = 2).
Voltage dependence of Ba2+-induced gating ring movement is related to function of the channel gate
Ca2+, Mg2+ and Ba2+ bind to the Ca2+ bowl and trigger conformational changes of the gating ring region (Miranda et al., 2016). However, the effects of these ions on BK function and gating ring motions are fundamentally different. Notably, Ba2+ induces a rapid blockade of the BK current after a transient activation that is measurable at low Ba2+ concentrations (Zhou et al., 2012; Miranda et al., 2016) (Figure 6a). In addition, we previously showed that the gating ring conformational motions induced by Ba2+ show a voltage-dependent component, which is not observed when Ca2+ or Mg2+ bind to the Ca2+ bowl (Miranda et al., 2013; Miranda et al., 2016) (Figure 6b). We combined mutagenesis with the cation-specific activation strategy to identify the structural source of the voltage dependence in Ba2+-triggered gating ring motions. In this case, alteration of VSD function by mutating charged residues (Figure 6c and e) was not reflected in any change of the E vs. voltage relationships, as shown in Figure 6d and f for constructs BK667CYR210E and BK667CYR213E, respectively. These results indicate that the voltage dependence of Ba2+-induced gating ring conformational changes, unlike those induced by Ca2+ and Cd2+ through activation of the RCK1 binding site, may not be related to VSD activation. This conclusion is further supported by the lack of changes in Ba2+ responses when mutations in the VSD were made in a RCK1 Ca2+ binding site knockout (D362A D367A) background (Figure 6—figure supplement 1b & c). Next, we studied the effect of Ba2+ on BK667CY channels containing the F315A mutation (Figure 6g) (Carrasquel-Ursulaez et al., 2015). As shown in Figure 6h, the E values reached similar levels to those of non-mutated BK667CY channels at saturating Ba2+ concentrations. However, at intermediate concentrations of Ba2+ the E-V curves were shifted towards more positive potentials when compared with BK667CY channels (Figure 6h, dashed line). These results suggest that the voltage-dependent component of the conformational changes triggered by Ba2+ binding to the Ca2+ bowl are not directly related to VSD activation, but rather to the function of the channel gate.
![](https://iiif.elifesciences.org/lax/40664%2Felife-40664-fig6-v2.tif/full/617,/0/default.jpg)
Voltage dependence of gating ring movements triggered by Ba2+.
(a) The RCK2 site is selectively activated by Ba2+, which additionally induces pore block. (b) FRET efficiency (E) data obtained at several Ba2+ concentrations from BK667CY constructs (Miranda et al., 2016). (c) Effect of the VSD R210E mutation after selective activation of the RCK2 binding site by Ba2+. (d) E-V curves obtained at several Ba2+ concentrations from BK667CYR210E constructs. (e) Effect of the VSD R213E mutation after selective activation of the RCK2 binding site by Ba2+. (f) E-V curves obtained at several Ba2+ concentrations from BK667CYR213E constructs. (g) Effect of the F315A mutation after selective activation of the RCK2 binding site by Ba2+ (h) E-V curves obtained at several Ba2+ concentrations from BK667CYF315A constructs. Data corresponding to each Ba2+ concentration are color-coded according to the legend at the bottom. For reference, the curve corresponding to 100 μM Ba2+ from the BK667CY construct shown in (b) is also shown as a colored dashed line in panels (b, d, f and h). Data points and error bars represent average ± SEM (n = 4–6; N = 2–3).
Discussion
Using fluorescently labeled BKα subunit constructs reporting protein dynamics between the RCK1 and RCK2 domains, we previously demonstrated that the channel high-affinity binding sites can be independently activated by different divalent ions, inducing energetically-additive rearrangements of the gating ring measured as changes in the FRET efficiency values, E (Miranda et al., 2013; Miranda et al., 2016). Further, the effects of Ca2+, Cd2+ and Ba2+ on the E values showed a voltage-dependent component, for which we could not provide an explanation. Voltage dependence of Ca2+-induced rearrangements seemed to be specifically related to RCK1 activation, since only the mutation of that binding site resulted in voltage-independent E signals (Miranda et al., 2016 and Figure 1). One possibility to explain this result is the existence of direct structural interactions of the RCK1 domain and the VSD. Interestingly, the recently obtained cryo-EM full BK structure from Aplysia californica revealed the existence of specific protein-protein interfaces formed by the amino terminal lobes of the RCK1 domains facing the transmembrane domain and the VSD/S4-S5 linkers (Hite et al., 2017). According to the structural data obtained in saturating Mg2+ and Ca2+ concentrations, gating of the channel by Ca2+ was proposed to be mediated, at least partly, by displacement of these interfaces causing the VSD/S4-S5 linkers to move, contributing to pore opening ((Hite et al., 2017; Tao et al., 2017); but see also (Zhou et al., 2017)). Our work provides functional data supporting this mechanism. Our data show that mutations altering the voltage dependence of BK VSD are reflected in the voltage dependence of the gating ring movements triggered by activation of the RCK1 binding site by Ca2+ or Cd2+. Mutations altering VSD function by inducing large leftward shifts in the Vh(j) values (Ma et al., 2006; Zhang et al., 2014) strongly correlate with negative shifts in the voltage dependence of the E signals. Likewise, mutations inducing positive shifts in the VSD voltage dependence of the voltage sensor function are reflected in E-V shifts towards more positive membrane voltages. Interestingly, we also observe a correlation between the changes in the slope of the G-V curves and that of the E-V curves (e.g. Figure 4f; see also Supplementary file 1), suggesting the existence of an interaction between the VSD and the gating ring. This idea is further supported by the effect of β1 which has been proposed to alter the voltage dependence of VSD function (Wallner et al., 1995; Cox and Aldrich, 2000; Nimigean and Magleby, 2000; Bao and Cox, 2005; Orio and Latorre, 2005; Contreras et al., 2012; Castillo et al., 2015). We observed that β1 and β3bNβ1 induce a leftward shift in the E-V curves. Conversely, two experimental strategies known to influence the G-V curves without direct interference with the VSD did not affect the voltage dependence of E. The lack of effect on the E-V curves of the mutation F315A can be explained because the shift in the G-V curves arises from the influence of this mutation in the C↔O transition with minor effects on the voltage dependence of the gating currents (Carrasquel-Ursulaez et al., 2015). Analogously, no change in the voltage dependence of E was observed after co-expression of BKα with the γ1 subunit, which shifts the voltage dependence of pore opening by enhancing its allosteric coupling with the voltage sensor activation (Yan and Aldrich, 2010). As with the mutation F315A, the presence of γ1 subunit produces a minor shift in the Q-V distributions, not paralleling the large shift in the G-V curves (Carrasquel-Ursulaez and Ramon Latorre, personal communication).
A puzzling result from our previous study was the observation that Ba2+ binding to the Ca2+ bowl triggers voltage-dependent conformational changes (Miranda et al., 2016). Even though we still do not know the mechanisms of this unique response to Ba2+, here we learned that it is not related to the dynamics of VSD, but rather influenced by perturbations affecting the opening and closing of the channel at the pore domain. Why Ba2+ but not Ca2+? A possible answer for this question is that Ba2+ has the additional property of blocking the permeation pathway (Miller, 1987; Neyton and Miller, 1988; Zhou et al., 2012), which could somehow be transmitted allosterically to the gating ring. If simply ion permeation blockade is what matters, then we might expect that blocking permeation with the high affinity quaternary ammonium derivative N-(4-[benzoyl]benzyl)-N,N,N-tributylammonium (bb-TBA) (Tang et al., 2009) should produce a voltage dependent FRET signal with Ca2+ activation. But, it does not (Figure 6—figure supplement 1d). Another possibility for the Ba2+ effect could be a direct allosteric interaction between the intrinsic gating in the pore and the divalent binding site in RCK2, which needs to be tested further.
Irrespectively of the fluorescent construct (Miranda et al., 2013) or the divalent ion used to activate the BK channel (Miranda et al., 2016), we have consistently observed that the conformational changes monitored as changes in the FRET efficiency are not strictly coupled to the intrinsic gating of the channel. In this study, we have found that the consequences of the voltage dependence of the intrinsic gating by manipulations of the VSD and the pore region are paralleled by the FRET efficiencies. These results rule out the possibilities that FRET signals derive from conformational changes in an unknown Ca2+ binding site or that they are completely uncoupled to the intrinsic gating.
In conclusion, our functional data show a strong correlation between the VSD function and the RCK1 conformational changes, suggesting a transduction mechanism from ion binding to change the channel activation. This transduction mechanism is in agreement with the existence of structural interactions between the RCK1 domain and the VSD. The correlation between VSD function and the RCK1 conformational changes is not observed between RCK2 and VSD, suggesting the existence of a different transduction mechanism that may include an indirect mechanism through the RCK1 or RCK1-S6 linker.
Materials and methods
Molecular biology and heterologous expression of tagged channels
Request a detailed protocolFluorescent BK α subunits were labelled with CFP or YFP using a transposon-based insertion method (Giraldez et al., 2005). Subunits labelled in the position 667 were subcloned into the pGEMHE oocyte expression vector (Liman et al., 1992). RNA was transcribed in vitro with T7 polymerase (Ambion, Thermo Fisher Scientific, Waltham, USA), and injected at a ratio 3:1 of CFP: YFP into Xenopus laevis oocytes, giving a population enriched in 3CFP:1YFP labelled tetramers (BK667CY) (Miranda et al., 2013; Miranda et al., 2016). Individualized Oocytes were obtained from Xenopus laevis extracted ovaries (Nasco, Fort Anderson, WI, USA). Neutralization of the Ca2+ bowl was achieved by mutating five consecutive aspartate residues to alanines (5D5A: 894–899) (Bao et al., 2002) on the BK667CY background. Elimination of RCK1 high-affinity Ca2+ sensitivity was achieved by double mutation D362A and D367A (Xia et al., 2002; Zeng et al., 2005; Zhang et al., 2010). Mutations were performed using standard procedures (Quickchange, Agilent Technologies, Santa Clara, USA). Auxiliary subunits (β3b, γ1 and chimera β3bNβ1) were co-injected with the BK667CFP/BK667YFP RNA mix at a 5:1 wt ratio, giving molar ratios above 20:1.
Patch-clamp fluorometry and FRET
Request a detailed protocolBorosilicate pipettes with a large tip (0.7–1 MΩ in symmetrical K+) were used to obtain inside-out patches excised from Xenopus laevis oocytes expressing BK667CY. Currents were recorded with the Axopatch 200B amplifier and Clampex software (Axon Instruments, Molecular Devices, Sunnyvale, USA). Recording solutions contained (in mM): pipette, 40 KMeSO3, 100 N-methylglucamine-MeSO3, 20 HEPES, 2 KCl, 2 MgCl2, 100 µM CaCl2 (pH 7.4); bath solution, 40 KMeSO3, 100 N-methylglucamine-MeSO3, 20 HEPES, 2 KCl, 1 EGTA, and MgCl2 or BaCl2 to give the appropriate divalent concentration previously estimated using Maxchelator software (maxchelator.standford.edu) (Bers et al., 1994). Solutions containing Cd2+ were prepared with a bath solution containing KF instead of K-Mes to precipitate the contaminant Ca2+ previously to the administration of the proper concentration of CdCl2 estimated with Maxchelator. Solutions containing different ion concentrations were exchanged using a fast solution-exchange system (BioLogic, Claix, France). All experiments were performed in various batches of oocytes, using different Ca2+ solutions prepared over time.
Simultaneous fluorescent and electrophysiological recordings were obtained as previously described (Miranda et al., 2013; Miranda et al., 2016). Conductance-voltage (G-V) curves were obtained from tail currents using standard procedures. The G-V relations were fit with the Boltzmann function: G/Gmax = 1/(1 + exp (-zF(V-Vhalf)/RT), where Gmax is the maximum tail current, z is the voltage dependence of activation, Vhalf is the half-activation voltage of the ionic current. T is the absolute temperature (295K), F is the Faraday’s constant and R the universal gas constant. Fit parameters are provided in Supplementary file 1. Conformational changes of the gating ring were tracked as intersubunit changes of the FRET efficiency between CFP and YFP as previously reported (Miranda et al., 2013; Miranda et al., 2016). Analysis of the FRET signal was performed using emission spectra ratios. We calculated the FRET efficiency as E=(RatioA-RatioA0)/(RatioA1-RatioA0), where RatioA and RatioA0 are the emission spectra ratios for the FRET signal and the control only in the presence of acceptor respectively (Zheng and Zagotta, 2003); RatioA1 is the maximum emission ratio that we can measure in our system (Miranda et al., 2013; Miranda et al., 2016). This value of E is proportional to FRET efficiency (Zheng and Zagotta, 2003). The E value showed is an average of the E value corresponding to each tetramer present in the membrane patch and represent an estimation of the distance between the fluorophores located in the same position of the four subunits of the tetramer. Where possible, the E-V relations were fit with the Boltzmann function: E = 1/(1 + exp (-zF(V-Vhalf)/RT), where z is the voltage dependence of the gating ring movement (E) and Vhalf is the half-activation voltage of the fluorescent signal. Fit parameters are provided in Supplementary file 1.
Data availability
All data generated and analysed during this study are included in the manuscript.
References
-
Elimination of the BK(Ca) channel's high-affinity Ca(2+) sensitivityThe Journal of General Physiology 120:173–189.https://doi.org/10.1085/jgp.20028627
-
Gating and ionic currents reveal how the BKCa channel's Ca2+ sensitivity is enhanced by its beta1 subunitThe Journal of General Physiology 126:393–412.https://doi.org/10.1085/jgp.200509346
-
Properties of single calcium-activated potassium channels in cultured rat muscleThe Journal of Physiology 331:211–230.https://doi.org/10.1113/jphysiol.1982.sp014370
-
A practical guide to the preparation of Ca2+ buffersMethods in Cell Biology 40:3–29.https://doi.org/10.1016/B978-0-12-374841-6.00001-3
-
Cloning and functional characterization of novel large conductance calcium-activated potassium channel beta subunits, hKCNMB3 and hKCNMB4Journal of Biological Chemistry 275:6453–6461.https://doi.org/10.1074/jbc.275.9.6453
-
Role of the beta1 subunit in large-conductance Ca(2+)-activated K(+) channel gating energetics. Mechanisms of enhanced Ca(2+) sensitivityThe Journal of General Physiology 116:411–432.https://doi.org/10.1085/jgp.116.3.411
-
Role of the S4 segment in a voltage-dependent calcium-sensitive potassium (hSlo) channelJournal of Biological Chemistry 273:32430–32436.https://doi.org/10.1074/jbc.273.49.32430
-
Generation of functional fluorescent BK channels by random insertion of GFP variantsThe Journal of General Physiology 126:429–438.https://doi.org/10.1085/jgp.200509368
-
Understanding the conformational motions of RCK gating ringsThe Journal of General Physiology 149:431–441.https://doi.org/10.1085/jgp.201611726
-
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
-
Coupling between voltage sensor activation, Ca2+ binding and channel opening in large conductance (BK) potassium channelsThe Journal of General Physiology 120:267–305.https://doi.org/10.1085/jgp.20028605
-
Molecular determinants of BK channel functional diversity and functioningPhysiological Reviews 97:39–87.https://doi.org/10.1152/physrev.00001.2016
-
Modulation of BK channel function by auxiliary beta and gamma subunitsInternational Review of Neurobiology 128:51–90.https://doi.org/10.1016/bs.irn.2016.03.015
-
Role of charged residues in the S1-S4 voltage sensor of BK channelsThe Journal of General Physiology 127:309–328.https://doi.org/10.1085/jgp.200509421
-
Trapping single ions inside single ion channelsBiophysical Journal 52:123–126.https://doi.org/10.1016/S0006-3495(87)83196-X
-
Discrete Ba2+ block as a probe of ion occupancy and pore structure in the high-conductance Ca2+ -activated K+ channelThe Journal of General Physiology 92:569–586.https://doi.org/10.1085/jgp.92.5.569
-
Differential effects of beta 1 and beta 2 subunits on BK channel activityThe Journal of General Physiology 125:395–411.https://doi.org/10.1085/jgp.200409236
-
Relative transmembrane segment rearrangements during BK channel activation resolved by structurally assigned fluorophore-quencher pairingThe Journal of General Physiology 140:207–218.https://doi.org/10.1085/jgp.201210807
-
BK potassium channels control transmitter release at CA3-CA3 synapses in the rat hippocampusThe Journal of Physiology 557:147–157.https://doi.org/10.1113/jphysiol.2004.062661
-
The contribution of RCK domains to human BK channel allosteric activationJournal of Biological Chemistry 287:21741–21750.https://doi.org/10.1074/jbc.M112.346171
-
Tetraethylammonium block of Slowpoke calcium-activated potassium channels expressed in Xenopus oocytes: evidence for tetrameric channel formationPflugers Archiv European Journal of Physiology 426:440–445.https://doi.org/10.1007/BF00388308
-
Intracellular Mg(2+) enhances the function of BK-type Ca(2+)-activated K(+) channelsThe Journal of General Physiology 118:589–606.https://doi.org/10.1085/jgp.118.5.589
-
Measurements of the BKCa channel's high-affinity Ca2+ binding constants: effects of membrane voltageThe Journal of General Physiology 132:491–505.https://doi.org/10.1085/jgp.200810094
-
Measuring the influence of the BKCa {beta}1 subunit on Ca2+ binding to the BKCa channelThe Journal of general physiology 133:139–150.https://doi.org/10.1085/jgp.200810129
-
Closed-channel block of BK potassium channels by bbTBA requires partial activationThe Journal of General Physiology 134:409–436.https://doi.org/10.1085/jgp.200910251
-
Cloning and functional expression of two families of beta-subunits of the large conductance calcium-activated K+ channelJournal of Biological Chemistry 275:23211–23218.https://doi.org/10.1074/jbc.M910187199
-
Characterization of and modulation by a beta-subunit of a human maxi KCa channel cloned from myometriumReceptors & Channels 3:185–199.
-
Activation of Slo1 BK channels by Mg2+ coordinated between the voltage sensor and RCK1 domainsNature Structural & Molecular Biology 15:1152–1159.https://doi.org/10.1038/nsmb.1507
-
Gating properties conferred on BK channels by the beta3b auxiliary subunit in the absence of its NH(2)- and COOH terminiThe Journal of General Physiology 117:607–628.https://doi.org/10.1085/jgp.117.6.607
-
Divalent cation sensitivity of BK channel activation supports the existence of three distinct binding sitesThe Journal of General Physiology 125:273–286.https://doi.org/10.1085/jgp.200409239
-
Allosteric regulation of BK channel gating by Ca(2+) and Mg(2+) through a nonselective, low affinity divalent cation siteThe Journal of General Physiology 118:607–636.https://doi.org/10.1085/jgp.118.5.607
-
Deletion of cytosolic gating ring decreases gate and voltage sensor coupling in BK channelsThe Journal of General Physiology 149:373–387.https://doi.org/10.1085/jgp.201611646
-
Threading the biophysics of mammalian Slo1 channels onto structures of an invertebrate Slo1 channelThe Journal of General Physiology 149:985–1007.https://doi.org/10.1085/jgp.201711845
Article and author information
Author details
Funding
National Institute of Neurological Disorders and Stroke (ZIA-NS002993)
- Pablo Miranda
- Miguel Holmgren
H2020 European Research Council (ERC-CoG-2014-648936)
- Teresa Giraldez
Ministerio de Economía y Competitividad (SAF2013-50085-EXP)
- Teresa Giraldez
Ministerio de Economía y Competitividad (RyC-2012-11349)
- Teresa Giraldez
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
MH and PM were supported by the intramural section of the National Institutes of Health (NINDS). TG was funded by the Spanish Ministry of Economy and Competitivity (grants SAF2013-50085-EXP and RyC-2012–11349) and the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement 648936). We thank Deepa Srikumar for technical assistance and Andrew Plested for useful comments on the manuscript. The γ1 clone and the β3bNβ1 chimera were kind gifts from Chris Lingle and Ramon Latorre, respectively.
Copyright
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Metrics
-
- 1,461
- views
-
- 247
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Structural Biology and Molecular Biophysics
The M2 proton channel aids in the exit of mature influenza viral particles from the host plasma membrane through its ability to stabilize regions of high negative gaussian curvature (NGC) that occur at the neck of budding virions. The channels are homo-tetramers that contain a cytoplasm-facing amphipathic helix (AH) that is necessary and sufficient for NGC generation; however, constructs containing the transmembrane spanning helix, which facilitates tetramerization, exhibit enhanced curvature generation. Here we used all-atom molecular dynamics (MD) simulations to explore the conformational dynamics of M2 channels in lipid bilayers revealing that the AH is dynamic, quickly breaking the 4-fold symmetry observed in most structures. Next, we carried out MD simulations with the protein restrained in 4-fold and 2-fold symmetric conformations to determine the impact on the membrane shape. While each pattern was distinct, all configurations induced pronounced curvature in the outer leaflet, while conversely, the inner leaflets showed minimal curvature and significant lipid tilt around the AHs. The MD-generated profiles at the protein-membrane interface were then extracted and used as boundary conditions in a continuum elastic membrane model to calculate the membrane bending energy of each conformation embedded in different membrane surfaces characteristic of a budding virus. The calculations show that all three M2 conformations are stabilized in inward-budding, concave spherical caps and destabilized in outward-budding, convex spherical caps, the latter reminiscent of a budding virus. One of the C2-broken symmetry conformations is stabilized by 4 kT in NGC surfaces with the minimum energy conformation occurring at a curvature corresponding to 33 nm radii. In total, our work provides atomistic insight into the curvature sensing capabilities of M2 channels and how enrichment in the nascent viral particle depends on protein shape and membrane geometry.
-
- Chromosomes and Gene Expression
- Structural Biology and Molecular Biophysics
The intricate structural organization of the human nucleus is fundamental to cellular function and gene regulation. Recent advancements in experimental techniques, including high-throughput sequencing and microscopy, have provided valuable insights into nuclear organization. Computational modeling has played significant roles in interpreting experimental observations by reconstructing high-resolution structural ensembles and uncovering organization principles. However, the absence of standardized modeling tools poses challenges for furthering nuclear investigations. We present OpenNucleome—an open-source software designed for conducting GPU-accelerated molecular dynamics simulations of the human nucleus. OpenNucleome offers particle-based representations of chromosomes at a resolution of 100 KB, encompassing nuclear lamina, nucleoli, and speckles. This software furnishes highly accurate structural models of nuclear architecture, affording the means for dynamic simulations of condensate formation, fusion, and exploration of non-equilibrium effects. We applied OpenNucleome to uncover the mechanisms driving the emergence of ‘fixed points’ within the nucleus—signifying genomic loci robustly anchored in proximity to specific nuclear bodies for functional purposes. This anchoring remains resilient even amidst significant fluctuations in chromosome radial positions and nuclear shapes within individual cells. Our findings lend support to a nuclear zoning model that elucidates genome functionality. We anticipate OpenNucleome to serve as a valuable tool for nuclear investigations, streamlining mechanistic explorations and enhancing the interpretation of experimental observations.