Nucleotide inhibition of the pancreatic ATP-sensitive K+ channel explored with patch-clamp fluorometry
Abstract
Pancreatic ATP-sensitive K+ channels (KATP) comprise four inward rectifier subunits (Kir6.2), each associated with a sulphonylurea receptor (SUR1). ATP/ADP binding to Kir6.2 shuts KATP. Mg-nucleotide binding to SUR1 stimulates KATP. In the absence of Mg2+, SUR1 increases the apparent affinity for nucleotide inhibition at Kir6.2 by an unknown mechanism. We simultaneously measured channel currents and nucleotide binding to Kir6.2. Fits to combined data sets suggest that KATP closes with only one nucleotide molecule bound. A Kir6.2 mutation (C166S) that increases channel activity did not affect nucleotide binding, but greatly perturbed the ability of bound nucleotide to inhibit KATP. Mutations at position K205 in SUR1 affected both nucleotide affinity and the ability of bound nucleotide to inhibit KATP. This suggests a dual role for SUR1 in KATP inhibition, both in directly contributing to nucleotide binding and in stabilising the nucleotide-bound closed state.
Introduction
ATP-sensitive K+ channels (KATP) couple the metabolic state of a cell to its electrical activity (Ashcroft and Rorsman, 2013). In pancreatic β-cells, closure of KATP in response to glucose uptake triggers insulin secretion. As such, mutations in KATP that affect its response to changes in cellular metabolism cause diseases of insulin secretion, for example neonatal diabetes and persistent hyperinsulinemic hypoglycaemia of infancy (PHHI; Quan et al., 2011; Ashcroft et al., 2017). KATP is composed of four inwardly rectifying K+ channel subunits (Kir6.2 in pancreatic β-cells), which form the channel pore, and four modulatory sulphonylurea receptor subunits (SUR1 in β-cells; Figure 1A; Aguilar-Bryan et al., 1995; Inagaki et al., 1995; Sakura et al., 1995; Inagaki et al., 1997). SUR1 is a member of the ABC transporter family but lacks any transport activity (Aguilar-Bryan et al., 1995; Tusnády et al., 1997). KATP responds to metabolism via adenine nucleotide binding to three distinct classes of intracellular nucleotide-binding site (one on each Kir6.2 subunit and two on each SUR1 subunit—making twelve sites in total [Vedovato et al., 2015]). Binding of ATP or ADP to Kir6.2 inhibits KATP channel activity (Tucker et al., 1997; Proks et al., 2010), whereas binding of nucleotides to SUR1 stimulates KATP (Nichols et al., 1996; Tucker et al., 1997). The stimulatory activity of nucleotides on KATP depends on Mg2+ (Gribble et al., 1998), whereas their inhibitory effect on Kir6.2 does not (Tucker et al., 1997).
In addition to nucleotide-dependent activation, SUR1 confers several other properties on Kir6.2. First, association with SUR1 increases the open probability () of Kir6.2 (Babenko and Bryan, 2003; Chan, 2003; Fang et al., 2006). Despite this increase in , SUR1 also paradoxically increases the apparent affinity for nucleotide inhibition at Kir6.2 by an unknown mechanism (Tucker et al., 1997). SUR1 is also responsible for high-affinity inhibition of KATP by antidiabetic sulphonylureas and glinides as well as activation by KATP-specific K+ channel openers (Tucker et al., 1997). Finally, SUR1 and Kir6.2 must co-assemble to ensure mutual exit from the endoplasmic reticulum and correct trafficking to the plasma membrane (Zerangue et al., 1999).
To date, the primary means of studying nucleotide-dependent effects on KATP channel function has been with electrophysiological approaches, which measure the summed activity of all three classes of binding site acting in concert. Thus, it can be difficult to separate the contributions of each class of site to the opening and closing of the channel pore and to properly distinguish between changes in nucleotide binding and channel gating. To overcome these limitations, we have applied a novel approach to directly measure nucleotide binding to each individual class of site in KATP (Puljung et al., 2019). This method utilizes Förster resonance energy transfer (FRET) between channels labelled with the fluorescent unnatural amino acid 3-(6-acetylnaphthalen-2-ylamino)−2-aminopropanoic acid (ANAP) and fluorescent trinitrophenyl (TNP) analogues of adenine nucleotides (Figure 1B). As we show here, this method is readily combined with patch-clamp electrophysiology so that nucleotide binding and regulation of current can be measured simultaneously. This has enabled us to quantitatively assess nucleotide binding to Kir6.2 and explore how this is coupled to channel inhibition in both wild-type KATP and KATP carrying mutations that impair ATP inhibition.
Results
Measuring nucleotide binding to Kir6.2
We previously used this FRET-based binding assay to measure nucleotide binding to the second nucleotide-binding site of SUR1 (Puljung et al., 2019). To measure binding to Kir6.2 in the complete KATP complex (four full-length Kir6.2 subunits co-expressed with four full-length SUR1 subunits), we replaced a tryptophan at position 311 (W311) of Kir6.2 that is 26 Å from the location of the inhibitory nucleotide-binding site with ANAP (Figure 1C) such that each subunit is labelled with one ANAP molecule. We designate this construct Kir6.2*. Based on the theoretical FRET efficiency calculated from the Förster equation and available cryo-EM structures (Martin et al., 2017; Martin et al., 2019), we expect 91% FRET efficiency between ANAP at position 311 and a TNP-ATP molecule bound to Kir6.2, and only 4% FRET efficiency to TNP-ATP bound to the closest nucleotide-binding site on SUR1 (nucleotide binding site 1; Figure 1C; Figure 1D). We also expect very little FRET between ANAP at position 311 and TNP-ATP bound to neighbouring Kir6.2 subunits (see Materials and methods for more details on our calculations).
ANAP incorporation into Kir6.2 was achieved as described previously (Chatterjee et al., 2013; Zagotta et al., 2016; Puljung et al., 2019). Briefly, HEK-293T cells were co-transfected with a plasmid encoding a Kir6.2 construct with a C-terminal GFP tag and an amber stop codon (TAG) replacing the codon corresponding to amino acid position 311 (W311TAG-GFP) and a plasmid encoding an ANAP-specific tRNA/tRNA synthetase pair (pANAP). We also included a dominant negative eukaryotic ribosomal release factor (eRF-E55D) in our transfections, which has been shown to increase the yield of full-length, ANAP-labelled protein (Schmied et al., 2014; Puljung et al., 2019). When cultured in the presence of ANAP, full length, fully ANAP-labelled Kir6.2 protein was produced and successfully trafficked to the membrane in the presence of SUR1 (Figure 1—figure supplement 1; see Materials and methods). We used GFP-tagged Kir6.2 constructs throughout this study unless otherwise indicated, to help identify cells or membranes expressing KATP.
Currents were measured from patches from cells expressing KATP that were excised in the absence of Mg2+. Under such conditions, nucleotides can bind to both sites on SUR1, but no activation occurs, allowing inhibitory currents to be measured in isolation (Gribble et al., 1998; Ueda et al., 1999; Puljung et al., 2019). Kir6.2*-GFP + SUR1 exhibited nearly identical sensitivity to ATP inhibition as Kir6.2-GFP + SUR1 (Figure 1—figure supplement 2A), indicating that replacement of W311 with ANAP did not affect inhibition of KATP. Whereas both constructs were inhibited by TNP-ATP with a higher apparent affinity relative to ATP, incorporation of ANAP resulted in channels with a slightly lower TNP-ATP sensitivity relative to wild-type ( of 6.2 µM compared to an of 1.2 µM, Figure 1—figure supplement 2B,C).
Kir6.2-GFP has been demonstrated to traffic to the plasma membrane in the absence of SUR1 and form functional channels (John et al., 1998; Makhina and Nichols, 1998). In a luminescence-based, surface-expression assay, we did not detect HA-tagged Kir6.2*-GFP at the plasma membrane in the absence of SUR1 (Figure 1—figure supplement 1E). To verify that the currents measured in our experiments in which Kir6.2*-GFP was co-transfected with SUR1 were the result of Kir6.2*-GFP + SUR1 and not Kir6.2*-GFP alone, we measured the sensitivity of currents to inhibition by the sulphonylurea tolbutamide, a property conferred by the SUR1 subunit. Whereas currents from unlabelled wild-type Kir6.2-GFP expressed in the absence of SUR1 were not affected by 100 µM tolbutamide, both wild-type Kir6.2-GFP and Kir6.2*-GFP currents were inhibited to a similar extent when expressed with SUR1 (46.5% ±0.04% and 57.7% ±0.02%, respectively; Figure 1—figure supplement 2D). The extent of inhibition was similar to previous measurements of tolbutamide inhibition (Tucker et al., 1997), confirming that Kir6.2*-GFP was co-assembled with SUR1 at the plasma membrane.
To measure nucleotide binding, cells transfected with Kir6.2*-GFP + SUR1 were briefly sonicated, leaving behind unroofed plasma membrane fragments (Heuser, 2000; Zagotta et al., 2016; Puljung et al., 2019) containing ANAP-labelled KATP channels with the intracellular nucleotide-binding sites exposed to the bath solution. The sample was excited with a 385 nm LED and emitted fluorescence from the membrane fragments was passed through a spectrometer, allowing us to separate ANAP, GFP, and TNP-ATP fluorescence by peak wavelength (Figure 1E,F). As expected from FRET, increasing the concentration of TNP-ATP caused a decrement in the ANAP peak at 472 nm and a concomitant increase in the TNP-ATP peak at 561 nm (Figure 1F). We used the quenching of the ANAP peak as a direct measure of TNP-ATP binding as this signal was specific to KATP. In contrast, the peak TNP-ATP fluorescence may include contributions from both specific and non-specific nucleotide binding, as well as direct excitation of TNP-ATP in solution by the 385 nm excitation light. Due to the sharp cut-off of the GFP emission spectrum at shorter wavelengths, our measurements of peak ANAP fluorescence were unaffected by the presence of the GFP tag on Kir6.2.
We fit concentration-response data for ANAP quenching by TNP-ATP with the Hill equation, to produce estimates of apparent affinity (, the half maximal effective concentration) and (ANAP quenching at saturating concentrations of TNP-ATP; Figure 1G). was 93%, in good agreement with the 91% predicted by the Förster equation and theoretical distance measurements (Figure 1D), suggesting that we were able to measure binding directly to the inhibitory site at Kir6.2. To confirm this, we introduced a well-studied neonatal diabetes mutation (G334D) into the Kir6.2 binding site, which drastically reduces the sensitivity of the channel to inhibition by nucleotides (Drain et al., 1998; Masia et al., 2007; Proks et al., 2010). Based on the cryo-electron microscopy structures of KATP, this mutation is expected to interfere with nucleotide binding directly (Figure 1H inset, Martin et al., 2017). The resulting construct Kir6.2*,G334D-GFP + SUR1 displayed drastically reduced ANAP quenching over the range of TNP-ATP concentrations tested. We therefore conclude that our binding measurements were specific for the inhibitory nucleotide-binding site on Kir6.2. This observation is consistent with the interpretation that the G334D mutation causes neonatal diabetes by preventing nucleotide binding.
Measuring current inhibition and nucleotide binding simultaneously
The apparent affinity of Kir6.2*-GFP + SUR1 for TNP-ATP in unroofed membranes was 25.6 µM (Figure 1G and Table 1). This value is higher than the apparent affinity for nucleotide inhibition (6.2 µM) measured using patch-clamp (Figure 1—figure supplement 2C). However, both binding and current measurement are a function of the intrinsic binding affinity, the channel , and the ability of agonist, once bound, to close the channel. Furthermore, the functional state of KATP in unroofed membranes is unclear. This is a particular problem with KATP channels, which run down due to slow dissociation of phosphatidylinositol 4,5-bisphosphate (PIP2), reducing the over time even in the absence of nucleotides (Proks et al., 2016).
As measuring either nucleotide binding or ionic currents in isolation only offers limited mechanistic insight into inhibition of KATP, we turned to patch-clamp fluorometry (PCF, Zheng and Zagotta, 2003). Using PCF, we can measure TNP-ATP binding to Kir6.2 and channel activity simultaneously (Figure 2), providing us with direct access to the relationship between nucleotide binding and channel function. We simultaneously measured fluorescence emission spectra and ionic currents for Kir6.2*-GFP + SUR1 in inside-out, excised membrane patches. The apparent negative fluorescence intensities at high TNP-ATP concentrations are due to imperfect background subtraction, and do not affect our measurements of ANAP intensities (see Materials and methods). As before, all measurements were performed in the presence of Mg2+ chelators, such that nucleotide inhibition could be measured in the absence of activation (Tucker et al., 1997; Gribble et al., 1998). Strikingly, current inhibition occurred at a lower range of concentrations compared to nucleotide binding (Figure 2C,D). The apparent for inhibition calculated from Hill fits was an order of magnitude lower than the for binding measured in the same patches (Figure 2C, Table 2). We considered several different gating models to explain this observation. In each model, we assumed the channel pore was able to open and close in the absence of ligand with an equilibrium constant , where and , reflecting the ability of KATP to open and close in the absence of nucleotides. This excludes the possibility of induced-fit models which would not predict unliganded channel closings. Induced-fit models also cannot account for separation between the binding and gating curves which we observe in Figure 2C,D (Changeux and Edelstein, 2011). Each model also had parameters representing the intrinsic binding affinity to the closed state (, where ) and the factor by which nucleotide binding favours channel closure (, where ).
Our simultaneous binding and current measurements were well fit with a Monod-Wyman-Changeux (MWC)-type model (Figure 2D,E; Monod et al., 1965) which has been previously proposed to explain KATP channel inhibition (Enkvetchakul and Nichols, 2003; Craig et al., 2008; Vedovato et al., 2015). In our MWC-type model, each ligand binding event () is independent and each bound ligand favours the closed state by the same factor (). Simultaneous measurement of binding (fluorescence) and gating (current) allowed us to obtain well constrained fits to our model. To obtain free parameter (, , ) estimates and verify that each parameter was well and uniquely determined, we employed a Bayesian Markov chain Monte Carlo (MCMC) method previously employed by Hines et al. (2014). Using this approach, we constructed posterior probability distributions for the free parameters of our MWC-type model (Figure 2F, Table 3). Based on these distributions, we estimated (), (), and . The very low value indicates that nucleotide binding was tightly coupled to channel closure; that is nucleotides have a very strong preference for the closed state of the channel. The low value for also explains why the channels were nearly completely inhibited at TNP-ATP concentrations at which not all the binding sites were occupied, as well as the degree to which channel inhibition is complete at saturating concentrations of TNP-ATP. Our estimate of was quite low and broadly distributed. We repeated our fits with fixed to a value consistent with previous single channel measurements (0.8, ; John et al., 1998; Enkvetchakul et al., 2000; Ribalet et al., 2006). This had only a very small effect on our estimates of and (Figure 2—figure supplement 1). The broad distribution of in our fit may represent current rundown which occurs during our patch-clamp recordings and is expected to affect the open-closed equilibrium. Cross-correlation plots (in parameter space) of the values derived from our fits produced well bounded ellipsoids, indicating that our parameters were uniquely determined (Figure 2—figure supplement 1A).
In addition to the full MWC-type model we considered alternate models (Figure 2—figure supplement 2). These included a model in which only the first binding event influences the open-closed equilibrium of the channel (single-binding model; Figure 2—figure supplement 2B, Table 3), and an MWC-style model with an additional parameter to allow for direct negative cooperativity between binding sites (negative cooperativity model; Figure 2—figure supplement 2C, Table 3). The single-binding model yielded very similar parameter estimates to our full MWC-type model (Figure 2—figure supplement 2D, Table 3). This is a consequence of being so low that even in the MWC-type model most channels are closed when only a single nucleotide is bound. The cooperative model improved our fits, but not enough to justify the inclusion of an additional free parameter (see Discussion).
Kir6.2-C166S affects the ability of bound nucleotides to close KATP
To provide a rigorous test as to whether our experimental system was capable of separating nucleotide binding from subsequent channel gating, we introduced a mutation (Kir6.2-C166S) which increases of KATP and decreases sensitivity of the channel to inhibition by nucleotides (Trapp et al., 1998). C166 is located near the bundle-crossing gate of Kir6.2 (Figure 3A). Other mutations at this site cause neonatal diabetes (Flanagan et al., 2006; Gloyn et al., 2006).
In unroofed membranes, Kir6.2*,C166S-GFP + SUR1 bound TNP-ATP with an very similar to that of Kir6.2*-GFP + SUR1 (Figure 3B, 32.0 µM and 25.6 µM, respectively), which suggests only a small change in apparent nucleotide affinity. This is an unexpected finding, as one might expect that an increase in would allosterically cause a decrease in the apparent affinity for inhibitory nucleotide binding. To resolve this conflict, we again turned to PCF (Figure 3C,D). Rundown was much slower for Kir6.2*,C166S-GFP + SUR1, which may reflect the increased of this construct. Measuring current inhibition in combination with nucleotide binding confirmed that whereas the apparent nucleotide binding affinity was unchanged by the C166S mutation, current inhibition occurred at much higher concentrations compared to binding and was incomplete (Figure 3D). How can we explain this paradox? Fits of the data with our MWC-type model (Figure 3D,E) suggest that, in addition to the expected effect on , the C166S mutation profoundly affects the ability of bound ligand to stabilise the closed state of the channel () without affecting (Figure 3E, Table 3). We propose that, in addition to increasing the of the channel, C166 is also important in the transduction pathway from the inhibitory nucleotide binding site on Kir6.2 to the channel gate.
Our MWC-type model predicts a nucleotide-insensitive current plateau at high concentrations, with the height of the plateau at saturating nucleotide concentrations given by . For example, when we see a current plateau of just under 5%, and as increases (as we see with Kir6.2*,C166S-GFP + SUR1) so does the plateau (Figure 3—figure supplement 1C,D). We attempted to confirm the existence of a nucleotide-insentitive current plateau for Kir6.2*,C166S-GFP + SUR1 at high concentrations of TNP-ATP, but were unable to test inhibition at concentrations of over 1 mM as our stocks of TNP-ATP were prepared from triethylammonium salts. Triethylamine concentrations of over 1 mM inhibited KATP and influenced our results (Figure 3—figure supplement 1A,B). However, we saw only partial inhibition of Kir6.2*,C166S-GFP + SUR1 by 10 mM ATP, which supports the existence of a nucleotide-insensitive current fraction at high concentrations. This observation has been previously reported for mutations at Kir6.2-C166 in some constructs (Trapp et al., 1998; Ribalet et al., 2006) but not others (Enkvetchakul et al., 2001).
Exploring the effect of SUR1 on nucleotide inhibition of KATP
SUR1 plays a complex role in the regulation of Kir6.2. It increases the of the channel and allows for the activation of the channel by Mg-nucleotides (Nichols et al., 1996; Tucker et al., 1997; Babenko and Bryan, 2003; Chan, 2003; Fang et al., 2006). However, it also increases the sensitivity of Kir6.2 to nucleotide inhibition (Babenko and Bryan, 2003; Chan, 2003; Fang et al., 2006). To understand the effect of SUR1 on nucleotide inhibition of KATP, we expressed Kir6.2*-GFP in the absence of SUR1 in unroofed membranes and measured TNP-ATP binding (Figure 4—figure supplement 1A). We found only a small increase (approximately 1.5-fold) in apparent compared to the same construct in the presence of SUR1 (37.6 µM and 25.6 µM respectively). Unfortunately, we were unable to achieve high enough expression of Kir6.2*-GFP alone to carry out PCF experiments in the absence of SUR1. However, we were able to measure currents from unlabelled Kir6.2-GFP alone (Figure 4—figure supplement 1B). As expected Kir6.2-GFP alone was much less sensitive to inhibition by TNP-ATP than Kir6.2-GFP + SUR1.
As Kir6.2*-GFP expression in the absence of SUR1 was not sufficient for PCF recordings, we took a mutational approach to better understand the role of SUR1 in inhibitory nucleotide binding. SUR1-K205 is located in the L0 linker of SUR1, which connects the first set of transmembrane domains (TMD0) to the ABC core structure (Figure 1A, Figure 4A; Martin et al., 2017; Puljung, 2018). This loop is adjacent to the inhibitory nucleotide-binding site on Kir6.2 and the interface between neighbouring Kir6.2 subunits. Mutations at SUR1-K205 were previously shown to reduce sensitivity of KATP to nucleotide-dependent inhibition (Pratt et al., 2012; Ding et al., 2019), and a recent cryo-EM stucture suggests that SUR1-K205 may directly coordinate the phosphates of ATP bound to Kir6.2 (Ding et al., 2019). Other mutations in L0 are associated with neonatal diabetes (Ashcroft et al., 2017) and PHHI (Snider et al., 2013).
We introduced a charge neutralization (alanine, K205A) and a charge reversal (glutamate, K205E) mutation at this position and measured simultaneous nucleotide binding and current inhibition with PCF (Figure 4B,C,D). The binding and inhibition curves for TNP-ATP almost perfectly overlaid for the SUR1-K205A mutant (Figure 4C). The same was also true for SUR1-K205E (Figure 4D). Data were fit with the MWC-type model as before. Mutating K205 to an alanine or a glutamate resulted in an apparent decrease in nucleotide binding affinity (Figure 4C,D,E). This was reflected by a decrease in the estimated for TNP-ATP, which correlated with the degree of conservation of the mutation, that is we observed a larger effect for the charge reversal compared to the charge neutralization mutation (Figure 4E). However, in addition to direct effects of K205 on nucleotide binding, we also observed a shift in for both mutations (Figure 4E). This suggests a dual role for SUR1 in KATP inhibition, both in contributing to nucleotide binding and in stabilising the nucleotide-bound closed state.
Discussion
We have developed a novel approach that allows for site-specific measurement of nucleotide binding to KATP and concomitant measurements of channel current. Performing these measurements simultaneously allowed us to examine nucleotide regulation of KATP function in great detail. We used a Bayesian approach to fit models to our combined fluorescence/current data sets to extract meaningful functional parameters with a minimum of prior assumptions. Such insights would not be possible from experiments in which macroscopic currents or binding were measured in isolation.
PCF has been used successfully by other labs to simultaneously measure ligand binding and gating in HCN channels (Biskup et al., 2007; Kusch et al., 2010; Wu et al., 2011). These groups measured fluorescence from a cyclic nucleotide analogue that increased its quantum yield when bound, minimizing background fluorescence from unbound ligand. Additional background subtraction could be performed by imaging the patches using confocal microscopy such that a region corresponding to the patch membrane could be computationally selected, thus omitting background fluorescence from the surrounding solution (Biskup et al., 2007; Kusch et al., 2010). In our PCF experiments, we used a FRET-based approach to measure ligand binding. We acquired fluorescence emission spectra, such that donor fluorescence could be separated from acceptor fluorescence by wavelength. This allowed us to directly assess binding from the quenching of donor fluorescence, which was specific to KATP. FRET also provided the spatial sensitivity necessary to discriminate between nucleotide binding directly to Kir6.2 and to the nucleotide-binding sites of SUR1. We assume that any TNP-ATP bound non-specifically to our membranes would be too far from Kir6.2 to cause appreciable FRET. This assumption was confirmed by the lack of FRET between TNP-ATP and a Kir6.2*-GFP mutant (G334D), in which nucleotide binding was severely disrupted (Figure 1H).
Previous studies have suggested that KATP inhibition follows an MWC-type model (Trapp et al., 1998; Enkvetchakul and Nichols, 2003; Drain et al., 2004; Craig et al., 2008; Vedovato et al., 2015). The majority of this earlier work was performed using single-channel measurements of mutated and/or concatenated channel subunits. In this study, we confirm these results using minimally perturbed channels with nucleotide sensitivity similar to that of wild-type KATP (Figure 1—figure supplement 2A). By using an MCMC approach to model fitting, we can also evaluate our models to assess how well the derived parameters were determined by the data. MCMC fits provide a basis for determining credible intervals for our parameter estimates. This allows for direct comparison of values derived from wild-type and different mutant constructs.
Although we did not explicitly include the effects of PIP2 on KATP gating in our model formulations, we assumed that the effects of PIP2 on were implicitly modelled in our parameter ; rundown due to dissociation of PIP2 manifests as a decrease in rather than a change in the number of channels. Although we were able to extract identifiable parameter estimates for , and , our estimates of for each model we considered were appreciably less well constrained than for the other parameters. We expect that this uncertainty arises from measuring a heterogeneous population of channels with regard to PIP2 binding. Fixing to values obtained from the literature (Figure 2—figure supplement 1, Figure 3—figure supplement 2, Figure 4—figure supplement 2, Figure 4—figure supplement 3) allowed us to extract estimates for and that were functionally identical to those derived from unconstrained fits, suggesting that the uncertainty of does not affect our inferences for these other parameters. Although the value for that we obtained from the literature for Kir6.2*-GFP + SUR1 (0.8) is on the extreme end of the posterior probability distribution for from our fits, we would not expect it to change our estimates of or as these values were stable over a broad range of values of (Figure 2—figure supplement 1) for this construct. Therefore, PCF represents a robust means to compare and between different mutated KATP constructs without worrying about the confounding effects of rundown.
Previous studies suggest that, whereas KATP closure occurs via a concerted mechanism, individual nucleotide binding events at Kir6.2 are not equivalent (Markworth et al., 2000). Earlier attempts to determine the stoichiometry of inhibitory nucleotide binding to Kir6.2 (i.e. how many ATPs must bind to induce channel closure) have produced models ranging from those in which binding of a single nucleotide completely shuts KATP to an MWC-type model in which each binding event is independent and contributes equally to channel closure (Trapp et al., 1998; Markworth et al., 2000; Enkvetchakul and Nichols, 2003; Drain et al., 2004; Wang et al., 2007; Craig et al., 2008; Vedovato et al., 2015). To resolve this controversy, we fit our data with both single-binding and MWC-type models. At very low values for , such as we derived from our experiments, the predictions of both models are functionally very similar. Even in our MWC-type model, we expect most KATP channels to be closed when just one molecule of nucleotide is bound.
It has been proposed that there is direct negative cooperativity between binding events at different subunits on Kir6.2 (Wang et al., 2007). We fit our data to an extended MWC-type model including an additional free parameter (), representing negative binding cooperativity between subunits (Figure 2—figure supplement 2). Not surprisingly this model improved the fit to our data as assessed by the Bayes factor, which represents the marginal likelihood of one model over another to explain our observations (Wagenmakers, 2007; Gronau et al., 2017). We also tested the cooperative model using approximate leave-one-out cross validation, which assesses the ability of a model to predict new or out-of-sample data using in-sample fits (Vehtari et al., 2017). Although in this work, we are primarily concerned with the inferences made from our fits, the ability of a model to make predictions is a good measure of its usefulness. Based on this criterion, the cooperative model has no more predictive accuracy than either the MWC-type model or the single-binding model. Therefore, the inclusion of an additional free parameter is not justified. Furthermore, whereas the cooperative model yielded good fits with identifiable parameters for Kir6.2*-GFP + SUR1 channels, it failed to yield identifiable parameters for all the mutants considered. Thus, this model did not allow for direct comparison between constructs. However, it remains a possibility that these mutations function in part by abolishing binding cooperativity between subunits.
We performed all our experiments on mutated, tagged channels using a fluorescent derivative of ATP. This allowed us to fit mechanistic models and readily compare between mutated constructs that affect nucleotide inhibition of KATP. This raises an obvious question: how relevant are our findings to inhibition of wild-type KATP by ATP? In a previous paper, we estimated and from an MWC-type model based on fits to published data for ATP inhibition of wild-type Kir6.2 + SUR1 (Proks et al., 2010; Vedovato et al., 2015). The value we obtained for (0.03) was quite similar to that we report here from our PCF measurements (0.04). We also obtained a similar estimate for in our previous model (3.0 × 104 M-1 vs 2.1 × 104 M-1 from our PCF experiments). Despite obtaining similar parameters, past experiments in which only ionic currents were measured, did not allow us to distinguish between competing gating models. Measuring currents and fluorescence simultaneously allowed for better model selection and aided in our ability to identify constrained parameters.
We compared the parameters derived for inhibitory nucleotide binding to those estimated for nucleotide activation of KATP based on experiments in which currents and binding were measured in separate preparations (Puljung et al., 2019). In those experiments, we estimated a value for , the factor by which binding of MgTNP-ADP to SUR1 stabilied channel opening, of 2.2. Although this value was derived using a different nucleotide, it still provides an approximate basis for comparing the coupling of nucleotide stimulation through SUR1 to nucleotide inhibition via binding to Kir6.2. If both activation and inhibition proceed via MWC-type models, the open closed equilibrium at saturating nucleotide concentrations is given by multiplied by or , respectively. The degree of stabiliation of the open state of KATP can be calculated as for activation. Stabiliation of the closed state is given by . Based on our observations, saturating concentrations of MgTNP-ADP stabilied the open state by −1.9 kcal mol-1 (−7.9 kJ mol-1). At saturating concentrations, TNP-ATP stabilied the closed state of KATP by −7.9 kcal mol-1 (31.8 kJ mol-1). Thus, assuming excitatory and inhibitory processes are independent, inhibition would be expected to dominate under conditions at which all the nucleotide binding sites are occupied. This is consistent with published measurements of wild-type KATP in the presence of Mg2+(Proks et al., 2010). In our previous study, we estimated for MgTNP-ADP binding to the stimulatory second nucleotide binding site of SUR1 to be 5.8 × 104 M-1 (), higher affinity than the we report here for TNP-ATP binding to the inhibitory site on Kir6.2 (2.1 × 104 M-1, ). Higher affinity binding to the stimulatory site may explain the ability of MgADP to increase KATP currents in the presence of ATP (Gribble et al., 1998). This phenomenon may also explain the bell-shaped MgADP concentration-response curve for KATP, which shows an increase in current at low concentrations, followed by inhibition at higher concentrations (Proks et al., 2010; Vedovato et al., 2015). Future experiments in which activation and inhibition are measured by PCF for the same ligand will allow us to model the complex response of KATP under conditions where all three nucleotide binding sites simultaneously affect channel gating (i.e. in the presence of Mg2+).
Mutations that cause neonatal diabetes reduce the sensitivity of KATP to nucleotide inhibition, and reduction in nucleotide sensitivity is broadly correlated with disease severity (McTaggart et al., 2010). We studied two residues on Kir6.2 that have been implicated in diabetes and have been proposed to affect nucleotide sensitivity via different mechanisms. We find that G334D drastically reduced the apparent affinity for nucleotide binding to KATP in unroofed membranes. In our MWC-type models, this could only be explained by a dramatic decrease in . This corroborates earlier hypotheses that mutating G334 directly disrupts inhibitory nucleotide binding to Kir6.2 (Drain et al., 1998). Due to poor expression, we were unable to test this construct using PCF. Therefore, we could not obtain accurate estimates of and .
In contrast to G334D, the C166S mutation does not directly affect nucleotide binding to Kir6.2, but rather disrupts the ability of bound nucleotide to close the channel. This contributes to the decreased nucleotide sensitivity which was previously attributed solely to an increased . In the future, we hope to use this rigorous approach to assess a whole panel of neonatal diabetes mutations in Kir6.2 to better understand the mechanism by which they cause disease.
Using PCF allowed us to probe more deeply into the role of SUR1 in regulating nucleotide inhibition of KATP. The cytoplasmic L0 loop of SUR1 was previously implicated in modulation of and nucleotide sensitivity of Kir6.2 (Babenko and Bryan, 2003; Chan, 2003; Pratt et al., 2012). We find that, in addition to directly contributing to tighter nucleotide binding at Kir6.2, SUR1 plays a critical role in preferentially stabilising the closed state of the channel when nucleotides are bound. Whereas a single nucleotide-binding event is sufficient for channel closure when Kir6.2 is associated with wild-type SUR1, mutating residue K205 reduced the ability of a single nucleotide to close the channel. This difference manifests in both our MWC-type and single-binding models (Figure 4—figure supplement 3 and Table 3).
In addition to providing mechanistic insights into disease-associated mutations in Kir6.2, our PCF-based approach allows us to probe the interactions between Kir6.2 and SUR1 on two different levels. As we show here, we can use this method to examine the effects of SUR1 on inhibitory nucleotide binding to Kir6.2. We can also adapt this method to study activation of Kir6.2 by nucleotides bound to the stimulatory sites on SUR1. Mutations in SUR1 that cause neonatal diabetes may do so by disrupting inhibitory binding/gating or enhancing the stimulatory effects of nucleotides. The formalism developed in this study provides a rigorous way to mechanistically assess the effects of these mutations. Our approach should be readily adaptable to the study of other nucleotide-gated channels including the cystic fibrosis transmembrane conductance regulator (CFTR, also an ABC-family protein) and purinergic P2X receptors.
Materials and methods
Molecular biology
Request a detailed protocolHuman Kir6.2 and SUR1 were subcloned into pcDNA4/TO and pCGFP_EU vectors for expression of wild-type and GFP-tagged constructs, respectively. pcDNA4/TO and pANAP were obtained from Addgene. peRF1-E55D and pCGFP_EU were kind gifts from the Chin Laboratory (MRC Laboratory of Molecular Biology, Cambridge, UK) and the Gouaux Laboratory (Vollum Institute, Oregon, USA) respectively. Amber stop codons and point mutations were introduced using the QuikChange XL system (Stratagene; San Diego, CA). All constructs were confirmed by DNA sequencing (DNA Sequencing and Services, University of Dundee, Scotland).
Cell culture and channel expression
Request a detailed protocolHEK-293T cells were obtained from and verified/tested for mycoplasma by LGC standards (ATTC CRL-3216, Middlesex, UK). Our working stock tested negative for mycoplasma contamination using the MycoAlert Mycoplasma Detection Kit (Lonza Bioscience; Burton on Trent, UK). Cells were plated onto either poly-L-lysine coated borosilicate glass coverslips (VWR International; Radnor, PA) or poly-D-lysine coated glass-bottomed FluoroDishes (FD35-PDL-100, World Precision Instruments). ANAP-tagged Kir6.2 constructs were labelled using amber stop codon suppression as described by Chatterjee et al. (2013). Transfections were carried out 24 hr after plating using TransIT-LT1 (Mirus Bio LLC; Madison, WI) at a ratio of 3 µl per µg of DNA. Unless specified otherwise, all transfections included a Kir6.2 construct with an amber stop codon (TAG) at position 311 (Kir6.2-W311TAG), SUR1, pANAP and eRF1-E55D in the ratio 0.5:1.5:1:1. Transfected cells were cultured in Dulbecco’s Modified Eagle Medium (Sigma; St. Louis, MO) + 10% foetal bovine serum, 100 U ml-1 penicillin and 100 µg ml-1 streptomycin (Thermo Fisher Scientific; Waltham, MA) supplemented with 20 μM ANAP (free acid, AsisChem; Waltham, MA). Cells were incubated at 33 °C and in the presence of 300 µM tolbutamide to enhance protein expression and channel trafficking to the plasma membrane (Yan et al., 2007; Lin et al., 2015). eRF1-E55D was included to increase efficiency of ANAP incorporation (Schmied et al., 2014). Experiments were carried out 2–4 days after transfection. We also expressed constructs labelled with ANAP at positions I182, F183, F198, and I210. Kir6.2-F183*, Kir6.2-F198*, and Kir6.2-I210* co-expressed with SUR1 did not produce sufficient currents for subsequent experimentation. Mutations at I182 are known to produce profound effects on nucleotide inhibition of KATP (Li et al., 2000). Thus, we did not consider this site for further experimentation.
Western blots
Request a detailed protocolTo confirm our ability to express full-length Kir6.2*-GFP, we performed western blots for HA-tagged Kir6.2 constructs in detergent-solubilized HEK-293T cells (Figure 1—figure supplement 1C). The HA tag plus a short linker (YAYMEKGITDLAYPYDVPDY) was inserted in the extracellular region following helix M1 of Kir6.2 between L100 and A101. Transfected HEK-293T cells grown in 6-well plates were harvested in cold PBS (Life Technologies Limited; Paisley, UK), pelleted at 0.2 x g for 2.5 min and resuspended in lysis buffer containing 0.5% Triton X-100, 100 mM potassium acetate, and a cOmplete protease inhibitor tablet (one tablet/50 ml, Roche; Basel, Switzerland), buffered to pH 7.4. After a 30 min benzonase (Sigma) treatment at room temperature, samples were mixed with a DTT containing reducing agent and loading buffer (NuPAGE, Invitrogen; Carlsbad, CA) and run on a precast Bis-Tris 4–12% poly-acrylamide gel at 200 V for 40 min. Proteins were wet transferred overnight onto polyvinylidene difluoride (PVDF) membranes (Immobilon P, Merck Millipore; Burlington, VT) in 25 mM Tris, 192 mM glycine, 20% methanol, and 0.1% SDS at 10 V on ice. Membranes were blocked with 5% milk in TBS-Tw (150 mM NaCl, 0.05% Tween 20, 25 mM Tris, pH 7.2) before staining for 30 min with a 1:1000 dilution of rat anti-HA monoclonal antibody in TBS-Tw (clone 3F10, Roche). After washing with TBS-Tw, membranes were incubated for 30 min with a 1:20,000 dilution of HRP-conjugated goat anti-rat polyclonal antibodies in TBS-Tw (Jackson ImmunoResearch; Ely, UK). Detection was performed using the SuperSignal West Pico Chemiluminescent Substrate (Thermo Fisher) and a C-DiGit Blot Scanner (Licor Biosciences; Lincoln, NE). Analysis was performed using custom code written in Python.
Transfection of wild-type Kir6.2-HA or Kir6.2-HA-GFP resulted in two bands on the western blots. The upper bands were close to the expected sizes for full-length Kir6.2-HA and Kir6.2-HA-GFP (46 kDa and 77kDa, respectively). We consistently observed a lower molecular weight band as well. This band must correspond to an N-terminally truncated Kir6.2 product, as the apparent molecular weight shifted with addition of the C-terminal GFP tag. Based on the molecular weight, we predict that the truncated protein product initiated from a start codon in the first transmembrane domain. Therefore, we believe it is unlikely that this protein would form functional channels or traffic to the plasma membrane. When Kir6.2-W311TAG-HA or Kir6.2-W311TAG-HA-GFP were co-transfected with SUR1, pANAP, and eRF1-E55D, and cells were cultured in the presence of ANAP, the western blots were similar to wild-type Kir6.2-HA or Kir6.2-HA-GFP. Over 90% full-length Kir6.2*-HA-GFP was produced under these conditions (Figure 1—figure supplement 1D). We were unable to quantify the percentage of full-length Kir6.2*-HA produced as the C-terminally truncated band resulting from termination at the TAG codon was very similar in size to the N-terminally truncated band. Co-expression with SUR1 increased the percentage of full-length Kir6.2*-HA-GFP produced (Figure 1—figure supplement 1D). In the absence of ANAP, we did not observe any full-length Kir6.2, indicating that there was no read-through of the amber (TAG) stop codon (Figure 1—figure supplement 1D).
Confocal microscopy
Request a detailed protocolConfocal imaging was performed using a spinning-disk system (Ultra-VIEW VoX, PerkinElmer; Waltham, MA) mounted on an IX81 microscope (Olympus; Southend-on-Sea, UK) with a Plan Apo 60x oil immersion objective (NA = 1.4), provided by the Micron Advanced Bioimaging Unit, Oxford. Transfected HEK-293T cells were incubated for 15 min with 1 nM CellMask Deep Red (Thermo Fisher) to stain plasma membranes before washing with PBS and imaging. ANAP was excited with a solid-state laser at 405 nm. GFP and CellMask were excited with an argon laser at 488 nm and 633 nm respectively. Images were captured on an EMCCD camera (ImagEM; Hamamatsu Photonics; Welwyn Garden City, UK) binned at 2 × 2 pixels and analysed using Python. A median filter with a box size of 32 × 32 pixels was applied to improve the signal-to-noise ratio by reducing background fluorescence.
We examined the surface expression of our ANAP-labelled constructs using confocal microscopy (Figure 1—figure supplement 1A,B). When Kir6.2-W311TAG-GFP was co-transfected with SUR1 along with pANAP and eRF1-E55D in the presence of ANAP, the ANAP and GFP fluorescence were co-localized at the plasma membrane. When wild-type Kir6.2-GFP was transfected under the same conditions, only GFP fluorescence was observed at the plasma membrane. ANAP fluorescence was diffuse and confined to the cytoplasm or intracellular structures. Thus, the plasma-membrane ANAP signal was specific for Kir6.2*-GFP.
Surface expression assays
Request a detailed protocolWe measured surface expression of HA-tagged Kir6.2 subunits using an approach outlined by Zerangue et al. (1999) ( Puljung et al., 2019). Cells were plated on 19 mm coverslips coated with poly-L-lysine and transfected as described above. Following incubation, cells were rinsed with PBS before fixation with 10% formalin for 30 min at room temperature. After washing again, cells were blocked with 1% BSA in PBS for 30 min at 4 °C before a 1 hr incubation at 4 °C with a 1:1000 dilution (in PBS) of rat anti-HA monoclonal antibodies. Cells were then washed 5 times on ice with 1% BSA in PBS followed by a 30 min incubation at 4 °C with a 1:2000 dilution of HRP-conjugated goat anti-rat polyclonal antibodies. Cells were washed 5 times in PBS + 1% BSA and 4 times in PBS. Coverslips were removed from the culture dishes and placed in clean, untreated dishes for measurement. 300 µl of SuperSignal ELISA Femto Maximum Sensitivity Substrate (Thermo Fisher) was added to each sample and the luminescence was measured using a Glomax 20/20 Luminometer (Promega; Madison, WI) after a 10 s incubation.
HEK-293T cells were transfected with Kir6.2 constructs with or without a TAG stop codon corresponding to position 311. Cells were co-transfected with pANAP and eRF1-E55D in the presence or absence of SUR1 and cultured with or without ANAP. Wild-type Kir6.2-HA and Kir6.2-HA-GFP in the presence of SUR1 were included as positive controls. Kir6.2 constructs with no HA tag served as negative controls. In the presence of ANAP, we observed strong trafficking of Kir6.2*-HA-GFP to the plasma membrane, but much less trafficking of Kir6.2*-HA (Figure 1—figure supplement 1E). When cells were cultured in the absence of ANAP, we observed little to no Kir6.2 surface expression from cells that were transfected with Kir6.2-W311TAG-HA or Kir6.2-W311TAG-HA-GFP, suggesting that prematurely truncated constructs did not traffic to the plasma membrane. In the absence of SUR1, surface expression was weak for both wild-type and tagged constructs, despite the reported ability of Kir6.2-GFP to traffic to the plasma membrane in the absence of SUR1 (John et al., 1998; Makhina and Nichols, 1998).
Epifluorescence imaging and spectroscopy
Request a detailed protocolEpifluorescence imaging and spectroscopy were performed using a Nikon Eclipse TE2000-U microscope with a 60x water immersion objective (Plan Apo VC, NA = 1.2, Nikon; Kingston upon Thames, UK) or a 100x oil immersion objective (Nikon, Apo TIRF, NA = 1.49). Imaging of ANAP was performed using a 385 nm LED source (ThorLabs; Newton, NJ) with a 390/18 nm band-pass excitation filter, an MD416 dichroic and a 479/40 nm band-pass emission filter (all from ThorLabs). GFP was imaged using a 490 nm LED source (ThorLabs) with a 480/40 nm band-pass excitation filter, a DM505 dichroic, and a 510 nm long-pass emission filter (all from Chroma; Bellows Falls, VT). Fluorescence spectra were collected by exciting ANAP as above but using a 400 nm long-pass emission filter (ThorLabs), then passing emitted light through an IsoPlane 160 Spectrometer (Princeton Instruments; Trenton, NJ) with a 300 grooves mm-1 grating. Images were collected with 1 s exposures on a Pixis 400BR_eXcelon CCD (Princeton Instruments).
Electrophysiology
Request a detailed protocolPatch pipettes were pulled from thick-walled borosilicate glass capillaries (GC150F-15, Harvard Apparatus; Holliston, MA) to a resistance of 1.5 MΩ to 2.5 MΩ when filled with pipette solution. Currents were recorded at −60 mV from excised inside-out patches using an Axopatch 200B amplifier equipped with a Digidata 1322A digitizer and using pClamp 10 software (Molecular Devices; San Jose, CA). Currents were low-pass filtered at 5 kHz and digitized at 20 kHz. The bath solution (intracellular) contained 140 mM KCl, 10 mM HEPES, 1 mM EDTA and 1 mM EGTA (pH 7.3 with KOH). The pipette solution (extracellular) contained 140 mM KCl, 10 mM HEPES and 1 mM EDTA (pH 7.4 with KOH). All experiments were carried out in Mg2+-free conditions. Currents were leak corrected using the current remaining in bath solution containing 5 mM barium acetate at 60 mV, assuming a linear leak with a reversal potential of 0 mV. Inhibition was calculated and corrected for rundown by alternating test concentrations of nucleotide solution with nucleotide-free solution, then expressing the test currents as a fraction of the average of the control currents before and after the test solution as described previously (Proks et al., 2010).
FRET calculations
Request a detailed protocolWe calculated the expected FRET efficiency between ANAP incorporated at amino acid position 311 and a docked TNP-ATP molecule as described previously (Puljung et al., 2019). The equivalency between FRET efficiency (measured as ANAP quenching) and nucleotide binding is based on two main assumptions. Firstly, we assume that the observed quenching from a bound nucleotide does not differ dramatically between open and closed states of the channel. As there is no open-state structure of KATP, we do not know exactly how much relative movement would occur between a bound TNP-ATP and Kir6.2-W311. However, based on cryo-EM structures of apo and nucleotide-bound Kir6.2 we do not expect to see a change in the distance between these two positions (Martin et al., 2019).
Secondly, we assume that the ANAP and TNP-ATP molecules on each subunit do not undergo energy transfer with those on other subunits to an extent which would dramatically change the observed quenching. At saturating TNP-ATP concentrations, where each ANAP-labelled site on Kir6.2 is occupied, FRET between ANAP and the closest acceptor will be kinetically favoured and the overall FRET efficiency will not be affected by cross-talk between neighbouring sites (Corry et al., 2005). In the limiting case, at low TNP-ATP concentrations, one would expect a large proportion of Kir6.2 tetramers (with four ANAP-labelled binding sites) bound to only a single TNP-ATP molecule. In this case, we expect a 4% overestimation of nucleotide binding as calculated using a numerical method to simulate a single TNP-ATP acceptor with multiple ANAP donors based on the distances calculated from our docking (Figure 1C; Deplazes et al., 2012). This may have resulted in our binding curves becoming artifically shallow at low concentrations. However, this difference is not significant in the context of our measurements as it is smaller than the observed error of our measurements at low TNP-ATP concentrations.
Unroofed binding measurements
Request a detailed protocolUnroofed membranes were prepared as described previously (Heuser, 2000; Zagotta et al., 2016; Puljung et al., 2019). A coverslip plated with transfected HEK-293T cells was removed from the culture media and rinsed with PBS. The coverslip was then briefly sonicated using a probe sonicator (Vibra-cell; Newtown, CT) leaving behind adherent plasma membrane fragments. Cells cultured on FluoroDishes were rinsed and sonicated directly in the dish. Unroofed membrane fragments were nearly invisible in bright-field images and identified by their GFP and ANAP fluorescence. Fluorescent TNP-nucleotides (Jena Bioscience; Jena, Germany) were diluted in bath solution and perfused onto unroofed membranes using a valve- controlled microvolume superfusion system (µFlow, ALA Scientific Instruments; Farmingdale, NY).
Fluorescence spectra were collected as described above. A region of interest corresponding to the membrane fragment was manually selected and line-averaged for each wavelength. A similarly sized region of background was selected and averaged, then subtracted from the spectrum of interest. After subtraction, ANAP intensity was calculated by averaging the fluorescence intensity measured between 469.5 nm and 474.5 nm. Bleaching was corrected by fitting the normalised ANAP intensity of exposures taken during perfusion with nucleotide-free solution to a single exponential decay of the form
then using the fit to correct the intensity of exposures taken during perfusion with test nucleotide solutions.
Some experiments were excluded from further analysis due to obvious cross-contamination between different solutions within the µFlow superfusion system. These were identified by noticeable colour changes in the solution in the delivery tubes.
Patch-clamp fluorometry
Request a detailed protocolThe tip of the patch pipette was centred on the slit of the spectrometer immediately after patch excision. Currents were measured as described above. Fluorescence emission spectra from the excised patch were acquired concurrently with current measurements, both during test solution application as well as nucleotide-free solution. Background subtraction was slightly imperfect due to the exclusion of TNP-ATP from the volume of the glass of the pipette, resulting in spectra that have negative intensities at the TNP-ATP peak at high nucleotide concentrations. However, this over-subtraction does not affect the size of the ANAP peak, which we used to quantify nucleotide binding.
ANAP bleaching was corrected as for the unroofed binding experiments with Equation 1 (Figure 2—figure supplement 3A). Due to the lower signal-to-noise ratio for PCF compared to the unroofed membranes, we performed experiments from both high-to-low and low-to-high TNP-ATP concentrations to minimise artifacts from our bleaching corrections. Kir6.2*-GFP + SUR1 showed consistent bleaching time courses (Figure 2—figure supplement 3B) and an average of 34% of the initial ANAP fluorescence intensity remained at the end of each experiment (Figure 2—figure supplement 3C).
Some experiments were excluded from further analysis due to low fluorescence intensity, as we were concerned about a low signal to noise ratio influencing our results.
Data processing and presentation
Request a detailed protocolRaw spectrographic images and current traces were pre-processed in Python and Clampfit (Axon) before analysis with R (R Development Core Team, 2019). Where applicable, all experimental data points are displayed in each figure. The number of experiments is reported in the figure legends and tables. To help visualise uncertainty and prevent some data points being hidden, points are arranged with a small amount of horizontal jitter; vertical position remains unaffected. Unless otherwise stated, summary statistics are overlaid as the mean with error bars representing the standard error of the mean. Where these error bars are not visible, they are smaller than the size of the point used for the mean.
Hill fits to fluorescence quenching were nonlinear least-squares fits to the following equation:
where represents normalised fluorescence intensity and and are values. Current inhibition data were fit to the same equation but with representing normalised current magnitude, instead of , and instead of .
Bayesian model fitting
Request a detailed protocolThe MWC-type models considered (Figure 2 and Figure 2—figure supplement 2) were formulated as follows:
When no ligand is present (i.e. when ), Equation 4 becomes:
We can use this to normalise the predicted changes in the open fraction to an observed change in current as:
Two variations on the full MWC model were also considered, and diagrammatic formulations are shown in Figure 2—figure supplement 2. The first was similar to the MWC-type model, except that the channels close after one molecule of TNP-ATP binding with subsequent binding events having no effect.
The second alternate model was the same as the full MWC model, but with an additional term describing binding cooperativity between Kir6.2 subunits.
Each model was fit to the combined patch-clamp fluorometry datasets using the brms package (Gelman et al., 2015; Burkner, 2017) in R. Prior probability distributions for each parameter were supplied as:
so that all priors are flat apart from L, which is weakly informative with 99% of its density falling between unliganded open probabilities of 0.01 and 0.99, and 85% falling between 0.1 and 0.9.
We considered two alternative sets of priors for fitting the MWC-type model to our mutant constructs (Figure 4—figure supplement 4). We generated a narrow informative prior by fitting normal distributions to the posterior probability density of our fits to Kir6.2*-GFP + SUR1, and a broad informative prior by increasing the standard deviation of the fitted normal distribution by a factor of ten (Figure 4—figure supplement 4A). Using narrow informative priors results in poorer fits as it does not allow for high enough values of to explain the data, whereas the broad informative priors result in fits which do not visibly differ much from the neutral priors (Figure 4—figure supplement 4B). Whereas the fits using the broad informative priors were visually similar to those using neutral priors, the posterior probability distributions for the parameters were slightly different (Figure 4—figure supplement 4C). Notably, due to the broad prior distribution supplied for , the posterior probability distribution for is also very broad for each construct, with a large amount of probability density over unlikely values, that is over half the probability density for Kir6.2*-GFP + SUR1-K205A and Kir6.2*-GFP + SUR1-K205E was for values lower than 0.01 (corresponding to a biologically implausible unliganded ).
Each model was run with four independent chains for 10,000 iterations each after a burn-in period of 20,000 iterations, saving every 10th sample for a total of 4000 samples per model. Each model parameter achieved a minimum effective sample size of 3500 and a potential scale reduction statistic () of 1.00. Where applicable, the posterior probabilities of each parameter are reported as the median and the 95% equal-tailed interval. Bayes factors were calculated using bridge-sampling (Gronau et al., 2017), and leave-one-out cross-validation (LOO-CV) was performed using the loo package (Vehtari et al., 2017).
Docking
Request a detailed protocolComputational docking of TNP-ATP into the nucleotide binding site of Kir6.2 was performed using AutoDock-Vina (Trott and Olson, 2010) and Pymol (Schrödinger, LLC; New York, NY). 11 TNP-ATP structures from the Protein Data Bank (PDB accession #s 1I5D, 3AR7, 5NCQ, 5SVQ, 5XW6, 2GVD, 5A3S, 2PMK, and 3B5J) were used as starting poses and a 15 Å × 11.25 Å × 15 Å box was centred on the ATP bound to Kir6.2 in PDB accession #6BAA (Martin et al., 2017). Protonation states for each residue were assigned using PDB2PQR and PROPKA 3.0 (Dolinsky et al., 2004). The modal highest-scoring pose from the docking run was selected (PDB accession #5XW6, Kasuya et al., 2017) and distances were measured from a pseudo atom at the centre of the fluorescent moiety. TNP-ATP (PDB #3AR7, Toyoshima et al., 2011) was positioned into the first nucleotide binding domain of SUR1 (PDB #6PZI, Martin et al., 2019) using the alignment tool in Pymol.
Chemicals and stock solutions
Request a detailed protocolUnless otherwise noted, all chemicals were obtained from Sigma. TNP-ATP was obtained as a 10 mM aqueous stock from Jena Bioscience and stored at −20° C aqueous stocks of ANAP-TFA were prepared by dissolving the free acid in 30 mM NaOH, and were stored at −20° C. Tolbutamide stocks (50 mM) were prepared in 100 mM KOH and stored at −20° C.
Data availability
All data and associated code are available on github (https://github.com/smusher/KATP_paper_2019; copy archived at https://github.com/elifesciences-publications/KATP_paper_2019) and have also been uploaded to Dryad (https://doi.org/10.5061/dryad.0vt4b8gtv).
-
Dryad Digital RepositoryNucleotide inhibition of the pancreatic ATP-sensitive K+ channel explored with patch-clamp fluorometry.https://doi.org/10.5061/dryad.0vt4b8gtv
References
-
Neonatal diabetes and the KATP channel: from mutation to therapyTrends in Endocrinology & Metabolism 28:377–387.https://doi.org/10.1016/j.tem.2017.02.003
-
K(ATP) channels and islet hormone secretion: new insights and controversiesNature Reviews Endocrinology 9:660–669.https://doi.org/10.1038/nrendo.2013.166
-
Sur domains that associate with and gate KATP pores define a novel gatekeeperThe Journal of Biological Chemistry 278:41577–41580.https://doi.org/10.1074/jbc.C300363200
-
Brms: an r package for bayesian multilevel models using stanJournal of Statistical Software 80:1–28.https://doi.org/10.18637/jss.v080.i01
-
Conformational selection or induced fit? 50 years of debate resolvedF1000 Biology Reports 3:19.https://doi.org/10.3410/B3-19
-
A genetically encoded fluorescent probe in mammalian cellsJournal of the American Chemical Society 135:12540–12543.https://doi.org/10.1021/ja4059553
-
How ATP inhibits the open KATP channelThe Journal of General Physiology 132:131–144.https://doi.org/10.1085/jgp.200709874
-
ExiFRET: flexible tool for understanding FRET in complex geometriesJournal of Biomedical Optics 17:011005.https://doi.org/10.1117/1.JBO.17.1.011005
-
PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculationsNucleic Acids Research 32:W665–W667.https://doi.org/10.1093/nar/gkh381
-
Concerted gating mechanism underlying KATP channel inhibition by ATPBiophysical Journal 86:2101–2112.https://doi.org/10.1016/S0006-3495(04)74269-1
-
ATP interaction with the open state of the K(ATP) channelBiophysical Journal 80:719–728.https://doi.org/10.1016/S0006-3495(01)76051-1
-
Gating mechanism of KATP channels: function fits formThe Journal of General Physiology 122:471–480.https://doi.org/10.1085/jgp.200308878
-
Stan: a probabilistic programming language for bayesian inference and optimizationJournal of Educational and Behavioral Statistics 40:530–543.https://doi.org/10.3102/1076998615606113
-
KCNJ11 activating mutations are associated with developmental delay, epilepsy and neonatal diabetes syndrome and other neurological featuresEuropean Journal of Human Genetics 14:824–830.https://doi.org/10.1038/sj.ejhg.5201629
-
Determination of parameter identifiability in nonlinear biophysical models: a bayesian approachThe Journal of General Physiology 143:401–416.https://doi.org/10.1085/jgp.201311116
-
Independent trafficking of KATP channel subunits to the plasma membraneThe Journal of Biological Chemistry 273:3369–3374.https://doi.org/10.1074/jbc.273.6.3369
-
The role of the KATP channel in glucose homeostasis in health and disease: more than meets the isletThe Journal of Physiology 588:3201–3209.https://doi.org/10.1113/jphysiol.2010.191767
-
On the nature of allosteric transitions: a plausible modelJournal of Molecular Biology 12:88–118.https://doi.org/10.1016/S0022-2836(65)80285-6
-
Engineered interaction between SUR1 and Kir6.2 that enhances ATP sensitivity in KATP channelsThe Journal of General Physiology 140:175–187.https://doi.org/10.1085/jgp.201210803
-
Activation of the K(ATP) channel by Mg-nucleotide interaction with SUR1The Journal of General Physiology 136:389–405.https://doi.org/10.1085/jgp.201010475
-
Running out of time: the decline of channel activity and nucleotide activation in Adenosine triphosphate-sensitive K-channelsPhilosophical Transactions of the Royal Society B: Biological Sciences 371:20150426.https://doi.org/10.1098/rstb.2015.0426
-
Cryo-electron microscopy structures and progress toward a dynamic understanding of KATP channelsThe Journal of General Physiology 150:653–669.https://doi.org/10.1085/jgp.201711978
-
Current understanding of K ATP channels in neonatal diseases: focus on insulin secretion disordersActa Pharmacologica Sinica 32:765–780.https://doi.org/10.1038/aps.2011.57
-
SoftwareR: A Language and Environment for Statistical ComputingR Foundation for Statistical Computing, Vienna, Austria.
-
Efficient multisite unnatural amino acid incorporation in mammalian cells via optimized pyrrolysyl tRNA synthetase/tRNA expression and engineered eRF1Journal of the American Chemical Society 136:15577–15583.https://doi.org/10.1021/ja5069728
-
Genotype and phenotype correlations in 417 children with congenital hyperinsulinismThe Journal of Clinical Endocrinology & Metabolism 98:E355–E363.https://doi.org/10.1210/jc.2012-2169
-
Molecular analysis of ATP-sensitive K channel gating and implications for channel inhibition by ATPThe Journal of General Physiology 112:333–349.https://doi.org/10.1085/jgp.112.3.333
-
Autodock vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreadingJournal of Computational Chemistry 31:455–461.https://doi.org/10.1002/jcc.21334
-
The Nucleotide-Binding sites of SUR1: a mechanistic modelBiophysical Journal 109:2452–2460.https://doi.org/10.1016/j.bpj.2015.10.026
-
Practical bayesian model evaluation using leave-one-out cross-validation and WAICStatistics and Computing 27:1413–1432.https://doi.org/10.1007/s11222-016-9696-4
-
A practical solution to the pervasive problems ofp valuesPsychonomic Bulletin & Review 14:779–804.https://doi.org/10.3758/BF03194105
-
Subunit-stoichiometric evidence for kir6.2 channel gating, ATP binding, and binding-gating couplingMolecular Pharmacology 71:1646–1656.https://doi.org/10.1124/mol.106.030528
-
State-dependent cAMP binding to functioning HCN channels studied by patch-clamp fluorometryBiophysical Journal 100:1226–1232.https://doi.org/10.1016/j.bpj.2011.01.034
-
Measuring distances between TRPV1 and the plasma membrane using a noncanonical amino acid and transition metal ion FRETThe Journal of General Physiology 147:201–216.https://doi.org/10.1085/jgp.201511531
Article and author information
Author details
Funding
Biotechnology and Biological Sciences Research Council (BB/R002517/1)
- Frances M Ashcroft
- Michael C Puljung
Wellcome (203731/Z/16/A)
- Samuel G Usher
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We wish to thank Raul Terron Exposito for technical assistance and Dr. Natascia Vedovato for helpful discussions. Dr. James Cantley provided access to the Licor scanner for western blots. This work was supported by the Biotechnology and Biological Science Research Council (BB/R002517/1) and the Wellcome Trust Oxion graduate program.
Copyright
© 2020, Usher et al.
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,682
- views
-
- 233
- downloads
-
- 19
- 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
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.
-
- Structural Biology and Molecular Biophysics
Osmotic stress and chloride regulate the autophosphorylation and activity of the WNK1 and WNK3 kinase domains. The kinase domain of unphosphorylated WNK1 (uWNK1) is an asymmetric dimer possessing water molecules conserved in multiple uWNK1 crystal structures. Conserved waters are present in two networks, referred to here as conserved water networks 1 and 2 (CWN1 and CWN2). Here, we show that PEG400 applied to crystals of dimeric uWNK1 induces de-dimerization. Both the WNK1 the water networks and the chloride-binding site are disrupted by PEG400. CWN1 is surrounded by a cluster of pan-WNK-conserved charged residues. Here, we mutagenized these charges in WNK3, a highly active WNK isoform kinase domain, and WNK1, the isoform best studied crystallographically. Mutation of E314 in the Activation Loop of WNK3 (WNK3/E314Q and WNK3/E314A, and the homologous WNK1/E388A) enhanced the rate of autophosphorylation, and reduced chloride sensitivity. Other WNK3 mutants reduced the rate of autophosphorylation activity coupled with greater chloride sensitivity than wild-type. The water and chloride regulation thus appear linked. The lower activity of some mutants may reflect effects on catalysis. Crystallography showed that activating mutants introduced conformational changes in similar parts of the structure to those induced by PEG400. WNK activating mutations and crystallography support a role for CWN1 in WNK inhibition consistent with water functioning as an allosteric ligand.