article: G protein-gated K+ channels (GIRK; Kir3), activated by Gβγ subunits derived from Gi/o proteins, regulate heartbeat and neuronal excitability and plasticity. Both neurotransmitter-evoked (Ievoked) and neurotransmitter-independent basal (Ibasal) GIRK activities are physiologically important, but mechanisms of Ibasal and its relation to Ievoked are unclear. We have previously shown for heterologously expressed neuronal GIRK1/2, and now show for native GIRK in hippocampal neurons, that Ibasal and Ievoked are interrelated: the extent of activation by neurotransmitter (activation index, Ra) is inversely related to Ibasal. To unveil the underlying mechanisms, we have developed a quantitative model of GIRK1/2 function. We characterized single-channel and macroscopic GIRK1/2 currents, and surface densities of GIRK1/2 and Gβγ expressed in Xenopus oocytes. Based on experimental results, we constructed a mathematical model of GIRK1/2 activity under steady-state conditions before and after activation by neurotransmitter. Our model accurately recapitulates Ibasal and Ievoked in Xenopus oocytes, HEK293 cells and hippocampal neurons; correctly predicts the dose-dependent activation of GIRK1/2 by coexpressed Gβγ and fully accounts for the inverse Ibasal-Ra correlation. Modeling indicates that, under all conditions and at different channel expression levels, between 3 and 4 Gβγ dimers are available for each GIRK1/2 channel. In contrast, available Gαi/o decreases from ~2 to less than one Gα per channel as GIRK1/2' s density increases. The persistent Gβγ/channel (but not Gα/channel) ratio support a strong association of GIRK1/2 with Gβγ, consistent with recruitment to the cell surface of Gβγ, but not Gα, by GIRK1/2. Our analysis suggests a maximal stoichiometry of 4 Gβγ but only 2 Gαi/o per one GIRK1/2 channel. The unique, unequal association of GIRK1/2 with G protein subunits, and the cooperative nature of GIRK gating by Gβγ, underlie the complex pattern of basal and agonist-evoked activities and allow GIRK1/2 to act as a sensitive bidirectional detector of both Gβγ and Gα.
  G proteins and the linked G protein-coupled receptors (GPCRs) are prominent regulators of excitability, which activate or inhibit ion channels by a variety of mechanisms [1]. This paper focuses on the quantitative analysis of the classical GPCR-initiated signaling cascade that culminates in the activation of GIRK channels (G protein-gated K+ channel; Kir3). GIRKs are important transducers of inhibitory neurotransmitter effects in heart and brain. They regulate heartbeat, neuronal excitability and plasticity, analgesia, alcohol and drug effects, and are implicated in a number of disorders such as epilepsy, Down syndrome, bipolar disorder, atrial fibrillation and primary aldosteronism [2,3, 4,5, 6]. GIRK is also the first-discovered effector of Gβγ [7] and a prototypical model of membrane-delimited G protein signaling. In the now classical scheme, the agonist-bound GPCR catalyzes GDP/GTP exchange at Gα and the separation of Gαi/oGTP from Gβγ; Gβγ directly binds to GIRK and triggers channel opening [8,9, 10,11]. Mammalian GIRKs are usually heterotetramers of GIRK1 with one of the other subunits (GIRK2, GIRK3 and GIRK4). GIRK1/2 is predominant in mammalian brain, but heterotetrameric GIRK1/3, GIRK2/3 and homotetrameric GIRK2 are also abundant in certain brain regions [2]. A GIRK channel is activated by direct binding of up to 4 molecules of Gβγ, but partial activation is achieved by binding of 1–3 Gβγ molecules [12,13,14,15,16]. NMR studies [17], crystal structure [18] and docking models [19] of GIRK-Gβγ complexes have confirmed the 4: 1 Gβγ: GIRK stoichiometry, showing binding of one Gβγ to each interface between adjacent GIRK subunits. Further, a strong association of GIRKs with Gβγ has been suggested by co-immunoprecipitation and Förster/Bioluminescence Resonance Energy Transfer (FRET/BRET, respectively) [20,21,22,23,24]. In support, in Xenopus oocytes, GIRK1-containing channels recruit Gβγ to the plasma membrane (PM) [25]. GIRK also binds Gαi/o subunits which regulate the channel' s basal activity, specificity and kinetics of signaling [26,27,28,29,30,31,32,33,34], but the mechanisms are poorly understood. No FRET between GIRK subunits and Gαi/o could be detected in the PM [24,35,36]; GIRK1 does not recruit Gαi to the PM [25] and binds Gαi/o in vitro less strongly than Gβγ [36]. The stoichiometry of Gα-GIRK interaction is unknown. Traditionally, GIRKs have been regarded as inhibitory devices operated exclusively by inhibitory neurotransmitters which elicit the GIRK' s evoked response (Ievoked). However, recent studies revealed that neuronal GIRKs also have a substantial basal activity, Ibasal [37,38,39]. GIRK' s basal activity and the balance between Ibasal and Ievoked are important determinants of neuronal excitability [39,40], bistability of neuronal networks [41], neuronal plasticity [42,43,44], dendritic integration [45], atrial arrhythmia and remodeling [46], and have recently been proposed to be related to effects of Li+, a drug used in the treatment of bipolar disorder [47]. Thus, changes in Ibasal and its relation to Ievoked are physiologically relevant and need to be understood. The molecular mechanisms of Ibasal and Ievoked have been extensively studied in heterologous model systems, mainly Xenopus oocytes and human embryonic kidney (HEK) cells (e. g. [48,49]). We discovered that, for the neuronal GIRK1/2, Ibasal and Ievoked are coupled. Incremental expression of GIRK1/2 in Xenopus oocytes revealed an inverse correlation between Ibasal and the extent of activation by transmitter. The higher Ibasal, the smaller the index of activation by the transmitter (Ra) and by coexpressed Gβγ (Rβγ) [30]. The Ibasal-Ievoked coupling was regulated by Gαi: coexpression of Gαi3 reduced Ibasal, increased agonist- and Gβγ-induced GIRK currents (a phenomenon we dubbed" priming" by Gαi/o”), and eliminated the inverse correlation between Ibasal and Ra [30,31,50,51]. These findings compelled an unusual explanation of the underlying mechanism. We proposed that Gβγ available for GIRK regulation is in excess over Gαi/o, thus the high Ibasal of GIRK1/2. We suggested that the magnitude of Ibasal and its relation to Ievoked are crucially regulated by the availability of Gαi/o [30,31,50]. Here we demonstrate that cultured hippocampal neurons show the same inverse relation between Ibasal and Ra as previously found in oocytes and HEK cells. This prompted us to further use these heterologous systems to address the coupling between GIRK' s basal and evoked activity. In the present work we have developed a quantitative model for Ibasal and Ievoked of GIRK1/2, which uses experimentally determined micro- and macroscopic parameters of GIRK1/2 currents and surface densities and accurately simulates and predicts macroscopic GIRK1/2 currents under a variety of conditions. Furthermore, modeling allowed to assess the apparent molar ratios of Gα and Gβγ available for GIRK, which we term “functional stoichiometry”. Our analysis reveals that, in Xenopus oocytes, HEK cells, and hippocampal neurons, 3 to 4 Gβγ molecules are available for the activation of GIRK1/2 channel over a wide range of surface densities of the channel, even when no exogenous Gβγ is coexpressed with GIRK. Calculations in Xenopus oocytes suggest a substantial increase in total concentration of Gβγ in the PM when large amounts of GIRK1/2 are expressed, corroborating the proposed mechanism of recruitment of Gβγ by GIRK1 to the PM [25]. In contrast, modeling shows that at most two Gα molecules are available for channel’s activation, even after overexpression of Gαi3. Furthermore, the Gα/GIRK ratio decreases with increasing channel density. The unequal and variable stoichiometry of GIRK1/2-associated Gα and Gβγ qualitatively and quantitatively explains the inverse Ra-Ibasal relation. Our results indicate a significant extent of association between GIRK1/2 and Gβγ, and support the notion that Gα is a non-obligatory partner in the GIRK-G protein signaling complex [50], but GαGDP plays a crucial role in regulating basal activity and, consequently, the magnitude of agonist response.
  First, we wanted to characterize the relation between GIRK’s Ibasal and Ievoked in hippocampal neurons, known to preferentially express GIRK1/2 [2]. Ibasal and Ievoked were measured using standard experimental paradigms ([39,47]; Fig 1A). Baclofen was used to activate the endogenous GABAB receptor and to generate Ievoked [52]. Net GIRK’s Ibasal was revealed as shown in Fig 1A (see also S1 Fig), by adding 100–120 nM tertiapin-Q (TPNQ), which selectively blocks >90% of GIRK currents in hippocampal neurons [39,43,53]. To characterize the relation between Ibasal and Ievoked, we utilized the activation index Ra, defined as Itotal/Ibasal (where Itotal is the total GIRK current; see Fig 1 and S1 Fig) [51]. GIRK currents of cultured hippocampal neurons showed considerable variability: Ibasal, 0. 2–26 pA/pF, Ievoked, 1–65 pA/pF (n = 65). Strikingly, there was a strong inverse correlation between Ra and Ibasal (Fig 1B, closed triangles), which was similar to that observed in oocytes expressing GIRK1/2 (open circles). The strength of the correlation indicates that it may be driven by a distinct molecular mechanism of potential physiological importance. The similarity of this distinctive phenomenon in hippocampal neurons and GIRK1/2-expressing oocytes encouraged us to further investigate it in the Xenopus oocyte expression system. The oocyte is particularly suitable for accurate control of protein expression (by titrating the injected RNA) and for current measurements, which are essential for quantitative modeling of Ibasal and Ievoked of GIRK1/2. Gβγ is well-established as the main gating agent for GIRK’s Ievoked [8,9]. This is also true for Ibasal of heterologously expressed GIRK1/2, which is suppressed by up to 80–90% by the expression of Gβγ-binding proteins such as C-terminus of β-adrenergic kinase, phosducin or Gα, both in Xenopus oocytes [30,31,50] and HEK293 cells [28]. In this work, we did not manipulate cellular levels of phosphatidylinositol diphosphate (PIP2), and used healthy cells which always showed robust GIRK currents, indicating levels of PIP2 sufficient for channel activation [54]. Thus, under the conditions used in this work, Gβγ was the main gating factor determining the steady-state macroscopic GIRK current (I). In a general form, I is described [1] by: I = Isingle·Po·N, (1) where Isingle is the single-channel current, N is the number of functional channels in the PM, and Po is the channel’s open probability. In a heterologous expression system, the channel’s surface density (N/S, where S is the surface area of the cell) can be experimentally manipulated and measured. Isingle of GIRK channels is an activation-independent parameter; Po is the gating parameter that changes as a function of the concentration of Gβγ available for GIRK activation by agonist or added Gβγ [9,11]. We start the development of the model by considering how Gβγ, available for activation of GIRK, can be derived from heterotrimeric Gαβγ (Fig 2A). In the absence of GPCR-activated G protein cycle, a small fraction of G proteins dissociates into free GαGDP and Gβγ due to finite affinity of their interaction [55,56] (the left branch of the reaction in Fig 2A). This free Gβγ can contribute to Ibasal [57,58]. Addition of agonist activates the GPCR and promotes GDP-GTP exchange at Gα and full or partial separation of GαGTP from Gβγ (the right branch of the reaction in Fig 2A; [59,60,61]). In our experiments in Xenopus oocytes and HEK293 cells, we coexpressed the muscarinic receptor 2 (m2R) which couples to Gi/o, and used acetylcholine (ACh) at supramaximal doses [62], 2–10 μM, in order to achieve a complete separation/rearrangement between GαGTP and Gβγ [63,64]. There are two existing models of GIRK gating by Gβγ. The allosteric kinetic model of Kurachi and colleagues, developed for cardiac GIRK1/4 [65,66,67], adequately describes the kinetics and magnitude of agonist- and GTPγS-evoked currents and the effect of RGS proteins. However, this model does not address Ibasal, does not include an explicit Gβγ binding step (Gβγ acts catalytically) and assumes very high surface densities of G proteins, in the order of 56 mM, which are incompatible with our measurements or those of others (see below and Discussion). It would be difficult to adjust this model for our purposes and to adequately describe Ibasal, or to implement the Gβγ recruitment phenomenon. The second model, termed here “concerted model”, was previously developed by us to describe the G protein-dependent activation of GIRK by Na+ (Fig 2B) [58]. The model included a description of Ibasal and an explicit Gβγ-GIRK binding step, but not the GαGTP-Gβγ dissociation step or Ievoked. Further, it assumed opening of the channel only when all four Gβγ binding sites are occupied (Fig 2B), which does not concur with the experimental findings that suggest a graded contribution of each bound Gβγ molecule [13,14,68]. Therefore, in the present work, we have developed a “graded contribution” model (Fig 2B and 2C), where each Gβγ-occupied state can contribute to channel opening and thus to Po [13,14,68]. To date it is not known whether Gβγ binding to GIRK is truly cooperative (i. e. whether Gβγ occupancy at one binding site increases the affinity of Gβγ binding at another site). Therefore, for simplicity, in the graded contribution model we assume that Gβγ binding to GIRK is sequential and the affinity of each Gβγ-binding site is independent of the occupancy of other sites. However, overall the process of gating is cooperative, since occupancy of each additional Gβγ-binding site increases Po in a more-than-additive manner (Fig 2C). The relative contributions of each Gβγ-occupied channel state to Po are adopted from published data for the homologous GIRK1/4 channel [13,14]. Thus, one bound Gβγ causes channel opening with a Po which is 1% of the maximal Po, Po, max; two Gβγ give 6%, three Gβγ 26%, and four Gβγ 100% of Po, max (Fig 2C). This approach is applicable only to steady-state calculations of macroscopic currents since it omits the kinetic details, but it allows to bypass the need to determine (or assume) a large number of unknown parameters: state-dependent changes in channel’s affinity to Gβγ, rates of closed-open transitions from different Gβγ-bound states, and the contributions of several potential open states. Once the channel achieved a state with n Gβγ bound, its fractional Po is known and does not depend on the pathway by which the channel opens. For simplicity, in calculating the steady-state Po for each Gβγ-occupied state, all open states (usually 2 are reported for GIRKs; [69,70,71,72]) were pooled into a single one (see Fig 2B). We have also considered a more general model with 4 separate closed states, in which each closed subunit can open independently of other subunits and the opening is promoted by Gβγ binding, giving rise to four open states (the “separate gating transitions model”, Fig 2B). The scheme shown also describes an alternative case in which subsequent closed states C1-C4 arise from the Gβγ-free closed state C0, the transitions between closed states are driven by Gβγ binding, and there are 5 separate C-O transitions. In both cases, it can be shown that, utilizing the approach based on graded contributions of each Gβγ-occupied channel state to Po, the derivation of Po converges to the same lead equation (eq 6) as for the graded contribution model (see Supplemental Discussion, S2 Text). Therefore, in this study we implemented the graded contribution model to simplify calculations. Throughout this work we also used an extended version of the concerted model, with the inclusion of the GPCR-induced dissociation of GαGTP from Gβγ, to cross-check the conclusions of the graded contribution model. Quantitative description and modeling of signaling cascades require the evaluation of amounts, stoichiometry, and affinities of interactions of participating proteins [73,74,75,76,77]. We took an approach that rests as much as possible on experimentally determined parameters. The data necessary for simulation by the model (eqs 5–13 in Methods) are the surface density of GIRKs and G proteins in Xenopus oocytes, macroscopic parameters of GIRK1/2 gating (whole-cell Ibasal, I evoked and the current induced by coexpression of Gβγ, Iβγ, see S2 Fig) at different channel densities, and Isingle and Po, max. Other parameters were experimentally determined in other works. To estimate single-channel current (Isingle) and Po of Gβγ- and agonist-activated GIRK1/2, we expressed the channels at low density with m2R and recorded channel activity in cell-attached patches (Fig 3A and 3B). Isingle was determined from amplitude distribution histograms (Fig 3A and 3B; right panels). The Gaussian fits to these histograms showed one main conductance level, suggesting that subconductance states, if any, did not significantly contribute to Po. The average Isingle was identical for ACh and Gβγ activation (Fig 3C), ~2. 8 pA. Po was estimated from patches containing 1 to 3 channels (see Methods). When channels were activated by Gβγ expressed at a saturating dose with no agonist present, Po was 0. 105±0. 018 (Fig 3D). The Gβγ RNA dose used (5 ng Gβ RNA, 1 ng Gγ RNA) consistently produced maximal macroscopic activation (see below, for example S8 Fig), and ACh generated negligible whole-cell Ievoked which was ~10% of Ibasal (Ra = 1. 1 ± 0. 02, n = 14). Thus, free endogenous GαGTP produced upon activation of m2R did not substantially affect the GIRK1/2 current evoked by saturating Gβγ. We therefore conclude that the Po measured in oocytes expressing saturating Gβγ is the Po, max of GIRK1/2, within a possible ~10% error. In comparison, when GIRK1/2 was activated via the coexpressed m2R (no Gα or Gβγ were coexpressed) with 2–5 μM ACh in the pipette, Po was 0. 037±0. 008, less than half of Po, max. (Fig 3B and 3D). (The actual Po could be higher because of the desensitization observed with ACh but not with Gβγ; see Methods). For further analysis and modeling of whole-cell Ibasal and Ievoked, we varied the surface density of GIRK1/2. The design was to obtain low, intermediate and high densities of GIRK1/2 by injecting 25,100–200 or 1000–2000 pg RNA of each subunit per oocyte. The cells expressed 1 or 2 ng of m2R RNA which did not affect Ibasal (S2C Fig and ref. [78]) but could always produce the maximal Ievoked [78]. Ievoked was elicited by ACh at 10 μM, a saturating dose. Under these conditions, all Gαi/o should convert to GαGTP, so that all available Gβγ can bind to the channel and activate it. The data are summarized in Table 1; main findings are also briefly highlighted in Fig 3E and 3F. We measured Ibasal, Ievoked and Itotal in each oocyte (set 1 in Table 1, S2A Fig). In separate groups of oocytes expressing saturating Gβγ, where channel’s Po reached Po, max, we measured Iβγ (set 2 in Table 1, S2B Fig). It is noteworthy that in oocytes, at all channel densities, Iβγ was 1. 6–2 fold greater than Itotal, the total GIRK current (Ibasal + Ievoked) without coexpressed Gβγ (Table 1). A similar Iβγ/Itotal ratio of 1. 66 was observed in HEK293 cells (Table 2), where all data have been pooled together (because GIRK1/2 expression levels have not been monitored). Similarly, Iβγ/Itotal ratio of ~2. 2 for GIRK1/2 expressed in HEK cells can be estimated from the data of Wydeven et al. [79] who activated GIRK with baclofen via GABAB receptors (Ievoked ~40 pA/pF, Iβγ ~ 90 pA/pF). The inverse value, Itotal/Iβγ, ranged 0. 5–0. 62 at different channel densities (Fig 3F, left panel). Since GIRK1/2 was maximally activated by the doses of Gβγ used in these experiments, Itotal/Iβγ is equal to the fraction of maximal activation after GPCR activation, Po/Po, max. Note that the single-channel data (Fig 3D) show less than 40% activation with saturating ACh (although, as noted above, this is probably an underestimate because of desensitization). In all, in oocytes and HEK cells, when the channel is activated by an agonist, only 40–60% of maximal Po is achieved. From the data of Fig 3, one can approximately estimate the amount of Gβγ molecules that are bound to the channel after maximal activation by agonist. This is done by comparing between measured values of Itotal/Iβγ (Fig 3F, left panel) and the expected Po/Po, max [13,14] from Fig 2C. To facilitate the comparison, we have redrawn the plot of Fig 2C in a simple graphic form (Fig 3F, right panel), and projected the values of Itotal/Iβγ onto the Po/Po, max plot (red dashed lines). For 3 bound Gβγ, the expected Po/Po, max is 0. 26, and for 4 Gβγ it is 1. Thus, with Itotal/Iβγ of 0. 4–0. 6, we estimate that, even without coexpression of Gβγ, between 3 and 4 Gβγ are available for activation of a single GIRK channel at all channel densities. For saturating Gβγ, Eq 1 for Iβγ takes the form: Iβγ= Isingle·Po, max·N. (2) From here, we calculated the total number of functional channels in the PM (N) and the corresponding channel density per μm2 of the PM. As shown in Table 1, our" low" ," intermediate" and" high" expression levels correspond to approximately 2. 7,9. 7 and 21. 7 channels/μm2, respectively. In the following, the data of Table 1 served as the basis for testing the predictions of the model and for calculating Gα and Gβγ available for channel activation. To obtain an independent estimate of the density of GIRK1/2 in the PM, we used quantitative immunochemistry of GIRK1 in cytosol-free, manually separated plasma membranes of Xenopus oocytes (Fig 4A and 4B) [80,81]. GIRK1 was coexpressed with GIRK2 at high density, and GIRK1 was probed with a C-terminally directed antibody. Western blots of manually separated PM and of the rest of the cells without the nucleus (" cytosol" ) showed the presence of two forms of GIRK1. A double band of about 55–58 KDa was always observed, and an additional higher diffuse band was seen in four out of seven experiments (Fig 4A and 4B). These bands correspond to partly and fully glycosylated channels, respectively [82,83]. Three oocyte batches showed only partly glycosylated bands in the PM. Since oocytes injected with 2 ng RNA always had large GIRK1/2 currents, it is likely that both partly and fully glycosylated channels are functional at the PM, in agreement with [83]. Notably, the main fraction of the channel was found in the cytosolic fraction (most likely endoplasmic reticulum and Golgi), largely in a partly-glycosylated form (Fig 4A). This is not unexpected, because in Xenopus oocytes the PM constitutes only a very small fraction of the cell' s total mass [81]. Next, molar amounts of GIRK in the PM fraction were calculated, taking into account the presence of two GIRK1 subunits in each channel. Calibration of antibody-produced signal was done with known amounts of the GST-fused distal C-terminus of GIRK1 which contains the epitope for the antibody. Note that this method yields channel levels in concentration units (e. g. mole/L). Both GIRK and Gβγ are associated with the PM (Gβγ is membrane-anchored by a lipid moiety [84]), and the interaction between GIRK’s cytosolic domain and Gβγ takes place within the submembrane space. Therefore, to compare data with GIRK1/2 surface densities obtained by channel counting from currents (Table 1), and for further modeling, we have converted two-dimensional protein densities to protein concentrations within the interaction space as previously described (e. g. [58,85,86,87]; see Supplemental Discussion, S2 Text), according to C = N/ (W·S·A), (3) where C is the concentration of protein in the submembrane space, N is the number of protein molecules in the membrane, S is a membrane area, A is Avogadro number and W is the width of the interaction space. For calculations, we used Soocyte = 2×107 μm2 (deduced from an oocyte' s capacitance of 200 nF [88] and specific capacitance of plasma membrane of 1 μF/cm2), and W was assumed to be 10 nm. The latter roughly corresponds to the molecular size of the complex of Gβγ and the cytosolic part of GIRK [18]. The influence of this parameter on the conclusions of the model was tested later (see below, S6 Fig, panels D and E). Consequently, the standard conversion factor between channel density (number of channels/μm2) and channel concentration (nM) is 1 channel/μm2 = 166 nM. Conversion of channel concentrations determined in Fig 4 into surface densities using Eq 3 gave ~12–14 channels/μm2 for both partly and fully glycosylated channels in the PM, and the average total amount of GIRK1/2 in PM (with partly + fully glycosylated GIRK1) was 19. 1±6. 8 channels/μm2 (Fig 4B). This is in good agreement with the independent assessment of ~22 channels/μm2 obtained from measurements of Iβγ for the high GIRK1/2 expression group (Table 1). To note, the latter was calculated using Eq 2 with Po, max measured at low channel densities. If Po, max were different at high GIRK1/2 expression levels, the densities calculated from Iβγ and immunochemistry would not match. The close correspondence between the two independent approaches indicates that Po, max is preserved at the high expression level. We conclude that the total surface density of GIRK1/2 channels in the PM can be satisfactorily estimated from whole-cell currents (Eq 2, Table 1). Such measurements are calibration-independent and accurate [1], and therefore GIRK1/2 can be used as a “molecular ruler”. In this procedure, the fluorescently labeled GIRK1/2, with its surface density calculated from Iβγ, will serve as a reference for estimating PM densities of other fluorescently labeled proteins. To use GIRK1/2 as a molecular ruler, we expressed YFP-GIRK1 (GIRK1 with Yellow Fluorescent Protein (YFP) fused to the N-terminus). Single channel analysis of Gβγ-activated YFP-GIRK1 coexpressed with GIRK2, YFP-GIRK1/2, showed the same Po, max and Isingle as in wild-type GIRK1/2 (S3 Fig), allowing the use of this construct for calibration purposes. To obtain high current levels of YFP-GIRK1/2 we usually had to inject 5 ng/oocyte of channel’s RNA. To estimate the expression levels of Gβγ using YFP-GIRK1/2 as molecular ruler, we expressed YFP-GIRK1/2 and, in separate oocytes of the same batch, Gβγ in which either Gβ or Gγ was labeled with YFP. Expression of YFP was monitored from fluorescence intensity in the PM (Fig 4C). In addition, Iβγ was measured and the surface density of YFP-GIRK1/2 was calculated. The amount of YFP molecules per μm2 was calculated assuming a 2: 2 GIRK1: GIRK2 stoichiometry in a heterotetramer [89] (See also Supplemental Discussion, S2 Text). The surface density of Gβγ-YFP was then calculated based on intensity ratios of YFP-GIRK1 and Gβγ-YFP. A typical experiment is shown and explained in Fig 4D. To validate the estimates of Gβγ expression, in four experiments as in Fig 4D we also measured the levels of Gβ-YFP (coexpressed with unlabeled Gγ) by quantitative immunochemistry in manually separated plasma membranes. We used purified Gβγ to calibrate the signal produced by the Gβ antibody (Fig 4E). We also constructed and expressed an YFP-fused construct corresponding to Xenopus laevis Gβ1, YFP-Gβ-XL (see Methods). Western blots showed a prominent ~36 KDa band of the endogenous Gβ, and ~70 KDa bands corresponding to the expressed YFP-Gβ or YFP-Gβ-XL (Fig 4E). The surface density of the expressed YFP-Gβγ assessed by the quantitative immunochemical method was 28±6 molecules/μm2, close to the estimate of 22. 1±8. 7 molecules/μm2 obtained in the same experiments from measurements of fluorescence using YFP-GIRK1 as" molecular ruler" (Fig 4F; 53 oocytes, n = 4 experiments; P = 0. 295). These results demonstrate the feasibility of the molecular ruler methodology and provide a good estimate of the expressed Gβγ-YFP. In several sets of experiments (see also below) we consistently found that, with 5 ng RNA of Gβγ, its surface density ranged between 20 and 30 molecules/μm2. We next utilized YFP-Gβ-XL as a caliper for the endogenous oocyte' s Gβ. Results of 4 experiments showed that, in Western blots, Gβ antibody used here gave similar signal with Gβ-XL as with bovine Gβ1 (S3 Fig, panels C, D). We then estimated the surface density of the endogenous Gβ (the 37 kDa band in Fig 4E) to be 24±4. 6 molecules/μm2 (n = 4). We have also estimated the concentrations of total and cytosolic endogenous Gβγ from the 4 experiments of Fig 4E, assuming an oocyte’s water volume of 0. 5 μl [88]. The total Gβγ concentration was 173±44 nM, the concentration of Gβγ in the cytosolic fraction was 171±44 nM. We define the molar amounts of proteins physically available for the function of the cascade as functional stoichiometry. It can change depending on availability of a protein, in contrast to limiting (maximal) stoichiometry which reflects the maximal molar ratios of interacting proteins. For example, if one GIRK channel can interact with at most 4 Gβγ and 4 Gα molecules, then the limiting GIRK: Gβγ: Gα stoichiometry is 1: 4: 4. Having determined the Po, max and surface densities of GIRK1/2 and endogenous Gβγ, we were now able to simulate macroscopic GIRK currents in oocytes and to assess the functional stoichiometry of GIRK1/2-Gα-Gβγ. We initially assumed that all of the Ibasal in oocytes was Gβγ-dependent. The affinities of GIRK-Gβγ and Gα-Gβγ interactions were adopted from published work: KD = 1. 86 nM for GαGDP-Gβγ binding [56], and KD = 50 nM for the GIRK-Gβγ interaction, as estimated by biochemical methods [90]. Simulations were done using eqs 5–16 as explained in the Methods section. Simulated data were compared to experimental measurements of GIRK1/2 activity for the three GIRK1/2 surface density groups (Table 1 and Fig 3E). Note that GIRK densities were calculated from Iβγ but simulations were done for the separately measured Ibasal and Ievoked, avoiding circular reasoning. We first tried to simulate the experimental data by assuming that only endogenous Gαβγ is available for the activation of GIRK1/2 (S4A Fig). However, no satisfactory description of data can be obtained under this assumption. Simulations that assumed recruitment of 3–4 Gβγ with GIRK, without Gα, gave a better approximation to the data (S4 Fig, panels B, C). Next, we turned to a more accurate assessment of functional GIRK1/2: Gβγ: Gα stoichiometry. Our model allows to calculate the amounts of Gα and Gβγ available for GIRK1/2 without any prior knowledge or assumption about the G protein concentrations in the cell, directly from experimental data. This idea is illustrated graphically in Fig 5 for the high GIRK1/2 surface density group of Table 1. The procedure consists of two steps. First, because all Gβγ available for GIRK is free to activate the channel after addition of agonist, the total GIRK-available Gβγ can be calculated from Itotal, as shown in Fig 5A. Here, the solid green line presents the simulated Itotal for a range of GIRK-available Gβγ densities, normalized to channel density (Gβγ: GIRK ratio). The intercept of simulated Itotal and the experimentally observed Itotal produces the estimate of the Gβγ: GIRK ratio. The numeric calculation of Gβγ density from experimental Itotal can also be done by substituting the channel density and Itotal into the eqs 5–13 and 16 (see Supplemental Methods (S1 Text) for Matlab routines). Next, the estimate of Gβγ is used to calculate the available Gα from Ibasal. Ibasal is determined by" free" Gβγ that is not bound to GαGDP in the absence of agonist. The calculation is done by substituting channel density, Ibasal and the estimated value of available Gβγ into the eqs 5–16 (see Supplemental Methods, S1 Text); a graphical illustration is shown in Fig 5B. The summary of calculations, made with our standard assumption of KD = 50 nM for the GIRK-Gβγ interaction, is presented in Table 3 and Fig 5C and 5D. These calculations show that there are ~3–4 Gβγ, but less than 2 Gα, available for each channel. Moreover, while Gβγ: GIRK ratio remains relatively constant throughout the range of analyzed channel densities, there is a sharp decrease in Gα: GIRK ratio with the increase in channel density (Table 3; see also below and S6 Fig). The estimate of KD for the GIRK-Gβγ interaction varies depending on the method used, from KD = 4–10 nM determined in excised patches [91,92,93] to ~50 nM in direct biochemical measurements [90]. To check for model’s stability regarding this parameter, we calculated the available Gβγ and Gα for a range of KD from 5 to 100 nM. As shown in Fig 5C and 5D, the estimates of Gβγ: GIRK and Gα: GIRK ratios, and the trend in their changes as a function of channel density, remain highly stable within the examined range of KD. Similar conclusions were attained using the concerted model (Table 3 and S5 Fig). The latter does not involve any assumptions for fractional Gβγ contributions to total Po. Similarity of conclusion of the two models alleviates concerns regarding the use of values for fractional Gβγ contributions, adopted from GIRK1/4 studies, to simulate GIRK1/2. Table 3 also shows the total Gβγ and Gα available for GIRK in PM, in molecules/μm2, calculated for each of the three channel densities. It is easy to see that, at intermediate and high channel densities, the endogenous Gβγ, 24 molecules/μm2 (that was present in the PM before the expression of GIRK) cannot account for the observed GIRK1/2 activation. In contrast, estimates of total Gα available for GIRK1/2 remain within the limits of endogenous Gα. As mentioned before, the calculations have been made without any assumption regarding the presence or amount of endogenous Gαβγ, and made no specific a priori assumptions regarding recruitment of Gβγ or Gα. Hence, modeling independently predicts the necessity for GIRK1/2-related increase in PM density of Gβγ, corroborating the experimentally observed recruitment of Gβγ, but not Gα, by GIRK1/2 [25]. Next, we addressed the possible contribution of intrinsic, Gβγ-independent activity to Ibasal. About 10–20% of Ibasal of GIRK1/2 in Xenopus oocytes [30,31,50] and HEK293 cells [28] persists after expression of Gβγ scavengers or Gα. Whereas the residual Ibasal may reflect incomplete Gβγ chelation, a genuine Gβγ-independent fraction of Ibasal cannot be discarded. In the extreme case it may contribute up to 20% of GIRK1/2 Ibasal. This may account for up to 10% of Po, max (because Ibasal can reach at most half of Iβγ, which is the indicator of Po, max; Tables 1 and 2). We have therefore extended the model to include the contribution of a hypothetical intrinsic Gβγ-independent channel activity. We assume that the intrinsic basal Po of a channel (Po, intrinsic) is an inherent, density-independent property of a single channel, best described as a fraction of Po, max. We thus repeated our calculations of GIRK1/2-available Gα and Gβγ assuming a Po, intrinsic in the range between 1% and 10% of Po, max (S6 Fig, panels A, B). For these calculations, Eq 6 (Methods) was modified in the following way: I=isingle⋅N⋅Po, max⋅ (ϕ⋅∑14fp, x⋅ϕx+ (1-ϕ) ), (4) where Φ is the fraction of Po, max which is Gβγ-dependent (see Eqs 6 and 7 in the Methods for definitions of other parameters). In the whole range of Po, intrinsic tested, the estimation of 3–4 Gβγ per channel remained highly stable (S6 Fig). The estimate of less than 2 Gα per channel also persisted except at the highest Po, intrinsic and low GIRK1/2 density, where Gα: GIRK ratio slightly exceeded 2 (S6B Fig, low surface density, black bar). Interestingly, for a significant Gβγ-independent intrinsic activity (10% of Po, max), up to ~60% of macroscopic Ibasal could be Gβγ-independent, especially at low channel densities which are common in native cells (S6C Fig). This finding may be relevant to some cells. For instance, coexpression of the Gβγ scavenger phosducin did not significantly reduce Ibasal in atrial cardiomyocytes [94], where the predominant channel is GIRK1/4. Finally, we considered the possible impact of variation in the presumed width (W) of the submembrane space within which the GIRK-Gβγ interactions occur. S6D and S6E Fig) shows that the main conclusions regarding the functional stoichiometry of GIRK, Gα and Gβγ remained largely unchanged over a wide range of W, 1–20 nm. We next tested the ability of the model to predict a new result: the dose dependency of activation by Gβγ, using the estimates of available Gα and Gβγ calculated from basal and agonist-evoked currents. We injected increasing amounts of wt-Gβγ mRNA into oocytes expressing GIRK1/2 at a constant density, and monitored both GIRK currents and Gβγ expression. Relative levels of Gβγ in the PM were directly measured in giant membrane patches of the oocytes [30,95] (Fig 6A) using the anti-Gβ antibody. Absolute surface densities of the expressed Gβγ (X axis in Fig 6C) were calculated assuming that 5 ng mRNA of Gβγ gives 30±4 molecules/μm2 (n = 47 oocytes). This density was calculated based on 3 experiments performed during the same period as the experiments of Fig 6 and S8 Fig, with wt-Gβ and YFP-Gγ, and using YFP-GIRK1/2 as the molecular ruler. As shown in Fig 6B, expression levels of Gβγ in the PM (grey bars) reached maximum at 5 ng RNA/oocyte. Iβγ (red circles) reached maximum already at 1 ng Gβγ RNA. Thus, maximal activation of GIRK1/2 has been attained already at submaximal expression levels of Gβγ (see also S8 Fig). Channel density of 13. 75 channels/μm2 was calculated based on Iβγ measured after expression of 5 ng RNA of exogenous Gβγ. Gβγ and Gα available to GIRK without the coexpression of exogenous Gβγ were calculated from Ibasal and Itotal (as in Fig 5), yielding ~43 molecules of Gβγ and 10 molecules Gα per μm2. We remind that the high density of available “endogenous” Gβγ in the presence of GIRK1/2 is due to Gβγ recruitment, explaining the high Ibasal and the relatively low index of activation of GIRK by Gβγ in a given oocyte, Rβγ. (Rβγ was defined as Iβγ/[average Ibasal], where average Ibasal was determined in a group of oocytes of the same experiment which expressed the channel without Gβγ. The definitions are as in [51]. See S2 Fig for definition of Rβγ). Using the calculated GIRK1/2 density and the amounts of available Gα and Gβγ before coexpression of Gβγ, we next calculated the predicted Iβγ and Rβγ for a range of doses (surface densities) of exogenously coexpressed Gβγ (Fig 6C). The predicted dose-dependencies of Iβγ and Rβγ (blue lines) are in agreement with experimental data (red circles). Assuming that 5 ng/oocyte of Gβ RNA gives either less (20 Gβγ/μm2) or more (44 Gβγ/μm2) molecules of coexpressed Gβγ instead of 30 Gβγ/μm2 produced similar predictions, still in good agreement with experiment (S7 Fig). Thus, the results of the simulations are relatively insensitive to a 50% variation in our estimate of coexpressed Gβγ. Further, very similar results were obtained in a separate experiment using a different experimental design, where we expressed increasing doses of Gβγ-YFP and calibrated Gβγ-YFP density using YFP-GIRK1/2 as molecular ruler (S8 Fig). We note that, because channel’s density is estimated from Iβγ obtained with a saturating dose of Gβγ, the good agreement between measured and predicted Iβγ at this RNA dose might be expected. However, the densities of Gα and Gβγ available to GIRK before coexpression of exogenous Gβγ (0 point on X-axis in Fig 6C) are calculated from Ibasal and Ievoked. Therefore, in both experiments (Fig 6 and S8 Fig), the satisfactory predictions of Iβγ and Rβγ at intermediate Gβγ doses, or the shape of the dose-response curves of Iβγ and Rβγ vs. Gβγ density, do not result from circular reasoning and are not trivial. This is illustrated by showing simulations that assume equal amounts of endogenous Gα and Gβγ (1,10 or 24 molecules/μm2) available for GIRK1/2, and no Gβγ recruitment. The use of these “classical” assumptions failed to reproduce the experimental data (black, red and green lines in Fig 6C and S7 Fig). In particular, saturation of Rβγ was predicted to happen at much higher doses of coexpressed Gβγ than in the experiment, obviously because the presumed initial basal level of Gβγ available to the channels was too low, thus requiring expression of more additional Gβγ. We next evaluated the model’s applicability to another expression system (HEK293 cells), and also to hippocampal neurons that natively express GIRK1/2 channels. We re-grouped raw data previously obtained in HEK cells expressing m2R and GIRK1/2 [51] (Table 2), the data from cultured hippocampal neurons (Fig 1), and the data obtained in oocytes, in a uniform manner. To enable direct comparison between the different systems, we arbitrarily segregated all the recordings into 4 groups, based on basal GIRK currents: <3 pA/pF, 3–13 pA/pF, 13–50 pA/pF and > 50 pA/pF (Table 4). For modeling, we needed to estimate channel surface densities, which have not been directly measured in HEK293 cells and neurons. To this end, we used Itotal to indirectly assess the channel densities. The ratio Iβγ/Itotal is fairly consistent in oocytes and HEK cells, ranging between 1. 6 and 2. 2 (Table 1 and ref. [96]). Thus, for each Ibasal range, we calculated Iβγ from Itotal assuming Iβγ/Itotal = 2 (Fig 7). Then we calculated densities using Eq 2 and Po, max of 0. 105, as measured in oocytes. The segregation of GIRK activity by Ibasal yielded relatively similar channel density groups in all cell types; neurons did not appear to express large amounts of GIRK, and the high density group was empty (Fig 7A). Next, using the procedure of Fig 5, we calculated the Gβγ and then Gα available for GIRK in all cases (Fig 7B and 7C). In all three systems, the relationship between G protein subunits and channel density was very similar to that found in oocytes. The most striking result was the persistent, channel density-independent availability of more than 3 Gβγ molecules per channel (Fig 7B). In all cases, the estimated Gα: GIRK ratio was about 2 for the low channel densities, but this number decreased as channel density increased (Fig 7C). Nevertheless, estimates of Gα: GIRK ratio in oocytes were lower than in HEK cells or neurons, indicating that there was a relative shortage of GIRK-associated Gαi/o in oocytes. We have therefore reanalyzed the experiments [30] in which Gα was co-expressed in oocytes in 5–10 fold excess (in terms of RNA quantities) over GIRK1/2 (Fig 7, denoted as “oocytes+Gα”). These doses produced the maximal GIRK-specific" priming" effect of Gαi3: strong reduction in Ibasal and increase in Ievoked without a significant reduction in Itotal (S1 Table). To note, coexpression of Gαi3 also produced Ra of ~10 which was comparable to the neurons with the lowest Ibasal and highest Ra (compare S1 Table and Table 4). Calculation of available Gβγ and Gα showed a robust persistence of Gβγ: GIRK ratio of above 3 (Fig 7B). Expectedly, the estimate of the available Gα greatly increased after Gαi3 overexpression, but, remarkably, Gα: GIRK ratio did not exceed 2 Gα molecules/channel (Fig 7C). Variations in Iβγ/Itotal ratio and Po, max in different cells could bias our estimates of channel density (Eq 2) and thus also estimates of Gβγ and Gα. Therefore, for neurons, we repeated our calculations for a range of Iβγ/Itotal ratios between 1. 5 and 3 (S9A Fig) and Po, max between 0. 05 and 0. 2 (S9C Fig). For comparison, a similar range of Iβγ/Itotal ratios was also tested for the oocyte data (S9B Fig). The exact values of Gβγ: GIRK and Gα: GIRK ratios varied, especially with changes in Iβγ/Itotal ratio. Generally, the lowest channel density is most sensitive to perturbations, and, for Iβγ/Itotal = 1. 5 (the lowest ratio tested), calculated Gβγ: GIRK and Gα: GIRK ratios exceed our usual estimates. However, this ratio is lower than that observed experimentally (Table 1), likely causing an overestimate of the values of Gβγ and Gα. In all, although the absence of direct measurements of channel densities and Po, max in HEK cells and neurons introduces an element of uncertainty, our results support the functional stoichiometry of 3–4 Gβγ and 2 or less Gα molecules per GIRK1/2 channel. Importantly, for a wide range of parameters, the main trends persist: available Gβγ is in excess over Gα; Gβγ: GIRK ratio remains high (>3) whereas Gα: GIRK ratio decreases as Ibasal increases. The systematic study presented above supports our hypothesis [51] that the inverse Ra-Ibasal relationship for GIRK1/2 reflects a progressive decline in GIRK1/2-associated Gα relative to Gβγ. We could now test this hypothesis quantitatively, and establish whether the calculated changes in GIRK: Gα: Gβγ functional stoichiometry can fully account for the Ra-Ibasal relationship shown in Fig 1. To this end, we simulated the changes in Ra as a function of Ibasal for a range of channel densities. We used channel densities, Ibasal values and Gβγ: GIRK and Gα: GIRK ratios calculated above (Tables 1 and 3 for oocytes, Table 4 for neurons). No free parameters were allowed. The results are shown in Fig 8, for oocytes (grey circles) and hippocampal neurons (black triangles). At this point, the channel density estimates and thus the simulations for oocytes are more reliable than for neurons. First, for simplicity, we assumed a constant Gβγ: GIRK ratio at all densities (3. 5 in oocytes and 3. 4 in neurons; see Fig 7). For further simulations, in order to construct continuous curves, we needed to assign numeric values for Gα: GIRK ratios within the full range of channel densities, based on the individual data points calculated for the “density groups” (Fig 8A). Since the reduction in Gα: GIRK ratio as a function of channel density is a process of unknown nature, the data points were arbitrarily fitted to a hyperbolic function (Fig 8A, solid lines). Next, we simulated the relation between Ibasal and channel density, by substituting the obtained values of Gα: GIRK into eqs 5–15 (Fig 8B). The simulation gave a good match to data of Table 1 (oocytes) and Table 4 (neurons), indicating that the fitting procedure of Fig 8A was satisfactory. Finally, values of Gα: GIRK from Fig 8A and Ibasal from Fig 8B were used to simulate the Ra-Ibasal relationship (Fig 8C, solid black and grey lines), matching well the raw data (triangles and circles). Simulations based on a constant relations of GIRK, Gα and Gβγ at all channel densities (allowing Gβγ and Gα recruitment) could not account for the observed trend in Ra changes. This is exemplified for a 1: 4: 2 and 1: 4: 4 GIRK: Gβγ: Gα stoichiometry (Fig 8C, red and green dashed lines, respectively). We conclude that the decrease in Gα available for GIRK activation at higher Ibasal can fully account for the inverse Ra-Ibasal relationship in both experimental systems.
  In this work we have quantitatively analyzed the GPCR-Gi/o-GIRK1/2 cascade, focusing on basal (Ibasal) and agonist-evoked (Ievoked) activities, both of which regulate neuronal excitability. We developed a mathematical model which allows quantification and simulation of macroscopic GIRK1/2 currents under steady-state conditions, before and after activation by neurotransmitter or by Gβγ. Our simulations fully rested on experimental data and parameters obtained in this and previous works. The modeling accurately described basal and evoked GIRK1/2 currents in two expression systems and in hippocampal neurons in a wide range of channel’s surface densities, correctly predicted the dose-dependent activation of GIRK1/2 by coexpressed Gβγ in Xenopus oocytes, and fully accounted for the inverse correlation between Ibasal and agonist activation index (Ra) previously observed in heterologous systems and, as shown here, also in hippocampal neurons. Our experimental findings and the model lay the basis for further analysis of the GPCR-Gi/o-GIRK cascade, for example for GIRKs of different subunit composition, and in different cells. Importantly, the present quantitative analysis provides novel and often unanticipated insights into the mechanisms of GIRK regulation by G protein subunits, Gβγ and GαGDP. It reveals an unequal and, moreover, variable functional GIRK1/2: Gβγ: Gα stoichiometry: 1) Under all conditions tested, between 3 and 4 Gβγ dimers are available for GIRK1/2; 2) Only two or less Gα are available per GIRK1/2 channel; 3) Increase of GIRK1/2 surface density is accompanied by a proportional increase in Gβγ (which is recruited by the channel), but not Gα. The unequal, effector-dependent Gα-Gβγ stoichiometry within the GIRK1/2 signaling cascade is an unexpected departure from classical schemes which usually assume that, prior to GPCR activation, the heterotrimeric G proteins available to the effector exist as stoichiometric complexes of Gα and Gβγ [97]. We propose that the unique functional stoichiometry of GIRK1/2 with G protein subunits, and the cooperative nature of GIRK gating by Gβγ, underlie the complex pattern of basal and agonist-evoked activities and allow GIRK1/2 to act as a sensitive bidirectional detector of both Gβγ and GαGDP. Our conceptual model of GIRK1/2 regulation (Fig 9) rests on the main findings of this study regarding the GIRK1/2: Gβγ: Gα stoichiometry (points 1–3 above) and the notion that, for Gβγ to activate GIRK, it must have its Gα-interacting interface exposed and free to contact GIRK [19,92,98,99]. In the resting state, the channel’s environment is enriched in 3–4 molecules of Gβγ and 1–2 Gα. In this scenario, between one and three Gβγ molecules are not associated with GαGDP and can bind and activate GIRK, resulting in a basal activity that is between 1 and 26% of total Po, max (see Fig 2C). Because of the gating cooperativity, occupancy of the first two Gβγ-binding sites yields low Ibasal. The fewer Gα, the more “free” Gβγ remains to occupy the activation sites at rest, yielding higher Ibasal. After GPCR-induced separation of GαGTP from Gβγ (lower arm of the scheme), due to gating cooperativity, addition of each Gα-free Gβγ ensures a robust 4–6 fold activation (e. g. going from 2 to 3 or 3 to 4 Gβγ-occupied sites). An even stronger activation takes place with a shift from 2 to 4 Gβγ-occupied sites (×16); and so on. Overexpression of Gβγ “sequesters” Gα and allows full occupancy of all Gβγ binding sites (middle arm of the scheme). Finally, overexpression of GIRK1/2 recruits Gβγ but not Gα, increasing Gβγ/Gα ratio (upper arm of the scheme). The balance between available GαGDP and Gβγ yields a continuum of basal activity magnitudes even on the level of a single channel, and sensitively regulates the extent of activation by the agonist. The scheme emphasizes the important fact that, given the relatively fixed amount of Gβγ available for GIRK1/2 activation, it is the availability of Gαi/o that determines the level of basal activity and, consequently, the extent of activation by agonists (as experimentally observed previously; [30,31]). The imbalance between Gβγ and Gα renders GIRK1/2 with a sizeable Ibasal, allowing it to act as a bidirectional, servo-like device [51] where its activity can be regulated not only by positive (Gβγ, Na+, PIP2) but also negative (GαGDP, protein kinase C, GαqGTP, PIP2 depletion) stimuli. Surface levels of endogenous G protein subunits and of heterologously expressed channels and G proteins crucially determine the behavior of the reconstituted signaling cascade, but they have never been quantitatively studied in the past in this common model system. We obtained very close estimates of surface levels of GIRK using two independent approaches: quantitative immunochemistry (which measures all channels in the PM) and electrophysiology (which counts only functional channels) (Fig 4). This indicated that the majority of GIRK1/2 channels in the PM of Xenopus oocytes were functional. Further use of GIRK1/2 as a molecular fluorescent ruler for Gβγ yielded Gβγ density very close to density estimated from quantitative immunochemistry, affirming the “molecular ruler” procedure with GIRK1/2 and lending additional support to the measurements of GIRK1/2 density. We estimated total cellular endogenous Gβγ in the oocyte as ~170 nM, similar to other cell types, 200–800 nM [100] and to the recent high-precision mass spectrometry measurement of ~200 nM Gβ in Xenopus eggs [101]. Notably, examination of data reported in the latter work [101] suggests a total concentration of all Gα in Xenopus eggs of ~350 nM. Thus, total oocyte' s Gβγ is not in excess over Gα, supporting our assumption (S4 Fig) that endogenous heterotrimeric G proteins are in the Gαβγ form before activation by GPCRs or coexpression of GIRK. Our estimate of 24 molecules/μm2 of endogenous PM-associated Gβγ, presumably as Gαβγ, is comparable to the ~40 molecules/μm2 evaluated in HEK cells [102]. In terms of concentration, 24 molecules/μm2 corresponds to ~4 μM, much higher than the cytosolic level of ~0. 2 μM. Such enrichment of G proteins at the PM is expected, because of the lipid modification of both Gα and Gγ [84,103]. It is probable that a substantial fraction of the PM-associated endogenous Gβγ is not available for GIRK activation, being associated with Gαs or Gαq rather than Gαi/o, or located in separate PM compartments, or associated with other effectors such as adenylyl cyclase [104]. Simulations showed that the main conclusions of our study are not affected by assuming a wide range of endogenous Gi/o proteins available for GIRK, from 1 to 24 molecules/μm2. In this work, we varied the levels of heterologously expressed GIRK1/2 and Gβγ. For GIRK1/2, the surface densities ranged from about 1 to 30 molecules/μm2 in oocytes and about 1–60 channels/μm2 in HEK293 cells (Table 4 and Fig 7). Densities of <20 channels/μm2 were estimated in cultured hippocampal neurons (Fig 7A). This is comparable to 9–10 channels/μm2 found in spines of cerebellar Purkinje neurons by quantitative electron microscopy [105]. Thus, the" low" and" intermediate" densities of expressed GIRK1/2 and Gβγ in oocytes and HEK293 cells may be the most physiologically relevant. For our description of GIRK currents, we have considered several models of GIRK channel gating by Gβγ (Fig 2B): the previously developed" concerted activation" model [58], the”graded contribution” model, and a more general model that includes 4 to 5 closed-open transitions. We were able to show that the latter two models converge to the same form of description of macroscopic steady-state Po based on fractional contributions of channels occupied by 1 to 4 Gβγ molecules, experimentally demonstrated for GIRK1/4 [13,14] (see Results and Supplemental Discussion S2 Text for a detailed discussion). We have therefore chosen the graded contribution model as our main tool to simulate and predict GIRK1/2 currents. Despite its relative simplicity, this model incorporates several complex properties of GIRK gating, and provides a strong computational tool for the analysis of the G protein-GIRK signaling. First, in this model we implemented the gating cooperativity of GIRK, by including the graded contribution of each bound Gβγ to channel opening (Fig 2C). Second, to our knowledge, this is the first model to describe both basal and agonist-evoked GIRK activity. Third, the model allows to estimate the amount of G protein subunits available for channel activation without any a priori assumptions regarding the levels of endogenous Gα or Gβγ; the GIRK-available Gα and Gβγ are calculated from experimental data (Fig 5). Finally, our method of calculating the functional stoichiometry of GIRK and Gβγ applies even if Ibasal is partly due to the presence of a GPCR activated by low dose of an ambient neurotransmitter. In this case, the available Gα calculated using the method of Fig 5 will represent only that fraction of GIRK-coupled Gα that is still in its GDP-bound form. In this work we left aside auxiliary/modulatory proteins such as RGS, focusing on the minimal essential composition of the cascade. We also have not addressed the impact of direct GIRK-Gα interaction. We and others did not find significant direct effects of GαiGDP on GIRK1/2 gating; the main function of GαGDP is the prevention of basal activation of the channel by ambient Gβγ and the release of Gβγ for channel activation by agonist/GPCR [31,34,36,106]. This function is fully implemented in our model. As for GαiGTP, it regulates the kinetics of Ievoked of GIRK1/2 but barely affects the steady-state amplitude [29,36]. Our present results suggest that, for GIRK1/2, effects of activated GαGTP on GIRK1/2 current amplitude are negligible, within a possible ~10% error. These considerations justify the omission of Gα-GIRK binding reactions from model’s equations. It remains possible that certain GαGTP may differently regulate GIRK channels of other subunit compositions or under certain conditions [34,107]. Both graded contribution and concerted models showed that Ibasal and Ievoked of the expressed GIRK1/2, and their changes with channel surface density, cannot be accurately described unless there is a recruitment of Gβγ to the PM by the channel. Recruitment of Gα is probably negligible (S4 Fig, Table 3). The prediction of preferential availability of Gβγ over Gα for GIRK1/2 concurs with the experimental findings [25], supporting the model’s validity. Further validation came from predicting the system’s response to perturbation in the form of dose-dependent response to coexpression of Gβγ, producing a satisfactory simulation of both Iβγ and Rβγ (Fig 6, S7 and S8 Figs). Finally, on the basis of experimental data and the calculated molar ratios of GIRK: Gβγ: Gα for different channel densities, the model fully accounted for the inverse correlation between Ra and Ibasal in GIRK1/2 (Fig 8), which was the starting point of this endeavor. We also made the first steps to extend the model to hippocampal neurons. The analysis of GIRK behavior in HEK cells and hippocampal neurons yielded estimates of Gβγ: GIRK and Gα: GIRK ratios close to those obtained in oocytes (3–4 Gβγ and ≤ 2 Gα), and a satisfactory simulation of the inverse Ra-Ibasal relation observed in the neurons (Fig 8). Uncertainties remain, because calculation of channel density in neurons relied on oocyte and partly HEK cell data for Iβγ/Itotal ratios and Po, max. A variety of GIRK compositions and distinct localization and density in cellular compartments further complicate the picture. Future quantitative studies are warranted for a more accurate description of GIRK activity in various neurons. Our analysis provides two major insights into functional stoichiometry of GIRK1/2 vs. the G protein subunits. First, the molar ratios are both uneven and can change as a function of channels’ density in the PM. The stoichiometry of more than 3 Gβγ per GIRK1/2 is practically invariable, whereas the amount of available Gαi/o is lower than Gβγ and further drops sharply as the level of expression of the channel increases. These estimates of Gβγ and Gα availability remained remarkably stable under a wide range of potentially variable parameters, such as KD of GIRK-Gβγ binding, the extent of Gβγ-independent basal activity, the size of submembrane reaction space, the Itotal/Iβγ ratio and Po, max in neurons, etc. Second, unexpectedly, the limiting stoichiometry for Gβγ: Gα: GIRK is 4: 2: 1. The uneven Gβγ: Gα stoichiometry and the decrease in Gβγ: Gα ratio were suggested by previous qualitative findings [50,51]. Our new results support this hypothesis and provide new insights into the underlying mechanism. Assumptions of pre-assembly of GIRK1/2 with 1,2, 3 or 4 Gαβγ heterotrimers failed to recapitulate the observed macroscopic currents and Ra (S4 Fig). The decrease in Gα/GIRK ratio as channel’s levels increase is consistent with total PM concentration of Gαi/o being relatively constant at all GIRK1/2 densities (Table 3). This is in agreement with little or no recruitment of Gα to the PM by GIRK1/2 [25,50]. On the other hand, the total amount of GIRK-available Gβγ in the PM increases as more GIRK1/2 channels are expressed, substantially exceeding the “basal” concentration of Gβγ of the naïve oocytes. The conspicuous persistence of Gβ: GIRK stoichiometry is best demonstrated in the oocytes, where it rests on a full quantitative analysis of experimental data. The identical estimates obtained in HEK cells and neurons (Fig 7B), though based on partial data, provide further support. The ability of GIRK1/2 to sustain a steady Gβγ-enriched environment strongly argues for a strong association between Gβγ and GIRK1/2, in line with the proposed high-affinity “anchoring” and recruitment of Gβγ by GIRK1 [25]. The mechanism of Gβγ recruitment is unknown but may be due to co-trafficking from the endoplasmic reticulum [22] or" kinetic scaffolding" and similar mechanisms [63,108,109], as discussed in [25]. At present we cannot rule out that the participation of another Gβγ-binding PM protein, such as a GPCR, or unknown scaffolding proteins, is important for the enrichment of Gβγ. It is widely accepted that the GIRK signaling cascade occurs within signaling complexes of GIRK channels with subunits of Gi/o heterotrimeric proteins and some GPCRs (reviewed in [4,61,110,111,112]). The uneven and variable stoichiometry within the GIRK1/2-Gβγ-Gαi signaling complex revealed by this study is compatible with a high-affinity, dynamic complex, where the channel and the G protein subunits are allowed to dissociate and reassociate in the PM [113]. These considerations justify our use of standard kinetic formalism for modeling. Finite affinity and reversibility of GIRK-Gβγ interaction is also supported by the demonstration of competition between two Gβγ effectors, GIRK and voltage-gated calcium channels, for available Gβγ in sympathetic neurons [114], and by recent GIRK2 reconstitution studies in lipid bilayers [106]. In contrast to Gβγ, our present analysis supports the notion [36,51,106] that, at least in heterologous or artificial systems, Gαi/oGDP is not an obligatory partner in the complex. Notably, in neurons and HEK cells the calculated Gα: GIRK ratios are higher than in oocytes (Fig 7 and Table 4); thus the higher Ra in these cells. The greater availability of Gα in HEK cells and neurons may reflect the presence of scaffolding or trafficking aids absent in the oocytes. However, the inverse Ra-Ibasal correlation and the (calculated) reduction in Gα at high channel densities are maintained in neurons and HEK cells, supporting the expendability of Gα. The apparent limiting stoichiometry of 4 Gβγ molecules per channel is not surprising, since the model explicitly includes 4 Gβγ-binding sites per channel. However, the calculated availability of two or less Gαi/oGDP, under most conditions examined, was unexpected. The limit of 2 Gα per channel was only slightly exceeded in simulations of lowest channel densities and when allowing substantial deviations from our standard assumptions (S9 Fig). Most conspicuously, overexpression of Gαi3, which reduced Ibasal by 75–80% and elevated the activation index Ra to about 10, increased the calculated Gα: GIRK ratio from less than 0. 5 to 2—but no more—Gα molecules per channel (Fig 7, S1 Table). Taken together, our results point to a limiting stoichiometry of 2 Gα molecules available for a GIRK1/2 channel. It is not clear what limits the amount of available Gα, but it is tempting to speculate that this limit reflects the maximal number of Gα molecules that can interact with GIRK1/2. The actual stoichiometry of this interaction is unknown, but the NMR study of Shimada and colleagues [115] indicates that the interacting surface of GαGTP requires two GIRK1 subunits of a GIRK1 tetramer for full contact. One GIRK1 subunit interacts with the helical domain of Gα and the other one with the GTPase domain [115]. We emphasize that our conclusions are valid for GIRK1/2 but may not be so for other GIRK channels. Thus, homomeric GIRK2 channels, with their low Ibasal and very high response to Gβγ, do not recruit Gβγ to the PM [25] and probably do not show pre-association with 3–4 Gβγ, or may have more Gα. Excess of Gβγ over Gα has been observed in the phototransduction cascade [116,117] though it has not been linked to any effector. We speculate that effector-dependent changes in the balance of Gα and Gβγ may take place with effectors other than GIRK, playing a role in their regulation.
  All experiments were approved by Tel Aviv University Committee for Animal Use and Care (permits M-08-081 and M-13-002 for Xenopus frogs and M-12-061 for mice). Female frogs, maintained at 20±2°C on 10 h light/14 h dark cycle, were anaesthetized in a 0. 17% solution of procaine methanesulphonate (MS222), and portions of ovary were removed through a small incision on the abdomen. The incision was sutured, and the animal was held in a separate tank until it had fully recovered from the anesthesia and then returned to the tank. The animals did not show any signs of postoperational distress and were allowed to recover for at least 3 months until the next surgery. Following the final collection of oocytes, anaesthetized frogs were killed by decapitation and double pithing. Xenopus oocytes were injected with RNA, and incubated in for 3–4 days at 20–22°C in NDE-96 solution (in mM: 96 NaCl, 2 KCl, 1 CaCl2,1 MgCl2,2. 5 Na-pyruvate, 50 μg/ml gentamycin, 5 mM HEPES/NaOH, pH 7. 5). All experiments with nerve cells derived from newborn mice, that have been analyzed in this paper, have been performed previously [47], and no additional mice have been used. The anti-GIRK1 polyclonal antibody was from Alomone Labs (Jerusalem), #APC-005. This antibody was raised against the distal C-terminal residues 437–501 of mouse GIRK1. 3 μg/mL were used for Western blots. The anti-Gβ polyclonal antibody (T-20, from Santa Cruz. #sc-378) is directed against the last 50 amino acids (a. a.) of mouse Gβ. 0. 4 μg/mL were used for Western blots. Most DNA constructs were as reported previously: bovine Gβ1, bovine Gγ2, human muscarinic type 2 receptor (M2R), human Gαi3, rat GIRK1, mouse GIRK2, YFP-GIRK1, GIRK2-HA, Gβ1-YFP, Gγ2-YFP and Gγ2-CFP [36,51,57]. To compare antibody labeling efficiency of bovine vs. Xenopus Gβ, we created the YFP-Gβ-XL construct in which the last 50 amino acids of the C-terminus in bovine Gβ were made identical to those of Xenopus Gβ1, by mutating 7 a. a. : V296I, A299C, A302R, D303E, A305V, A309S and D322S by standard PCR protocols. Preparation and storage of Gβ1γ2 and of the GST-fused distal C-terminus of GIRK1, GST-dCT (a. a. 365-501) used for calibrations of Fig 4A were done as described previously [69,118]. RNA was synthesized in vitro [57]. Amounts of injected RNA are indicated in the text, Tables 1, S1 and in Figure legends. Fluorescence levels of the expressed YFP (yellow fluorescent protein) and cerulean (termed here CFP, cyan fluorescent protein) were measured in intact Xenopus oocytes essentially as described [25,32]. Both YFP and CFP carried mutations that increase stability and reduce dimerization [32]. Briefly, oocytes were imaged in ND96 solution in a 0. 7 mm glass-bottom dish using Zeiss 510META confocal microscope with a 20× air objective. Images were acquired in the spectral mode. CFP was excited at 405 nm and emission was measured at 481–492 nm. YFP was excited at 514 nm and emission was measured at 535–546 nm. Fluorescent signals were averaged from 3 regions of interest (ROI) using Zeiss LSM Image Browser, and averaged background measured at an area outside the cell was subtracted. The average signal from uninjected oocytes was subtracted for final analysis. Saturation of emission measurement was strictly avoided to ensure that the readout of the confocal microscope was linear within the range of measurement. All measurements were made in the linear range of the recording apparatus. In calculating Gβγ surface density through comparing fluorescent signals from YFP-GIRK and Gβγ-YFP, the fluorescent intensities of the two proteins were compared directly from oocytes of the same batch on the same day, as described in Fig 4D. No correction for non-fluorescent (improperly folded) YFP [119] was needed (assuming that the percent of misfolding was similar in all YFP fusion proteins used here), because the ionic current Iβγ (a fluorescence-independent parameter) was used as the basis for GIRK density estimates. Giant membrane patches of oocyte membrane were prepared and imaged as described [25,95]. For imaging, fixated membranes were immunostained with the anti-Gβ antibody at 1: 200 dilution. Cy3-conjugated anti-rabbit secondary antibody (Jackson ImmunoResearch) was used for imaging with the 543 nm laser, and excitation was measured at 566–577 nm. Background fluorescence from an area outside the giant patch was subtracted. For final analysis shown in Fig 6, the signal from uninjected oocytes was subtracted from all groups. Manual separation of plasma membranes from the rest of the oocyte (“cytosol”) was performed as described [81], with modifications. Plasma membranes together with the vitelline membranes (extracellular collagen-like matrix) were removed manually with fine forceps after a 5–15 min incubation in a low osmolarity solution (5 mM NaCl, 5 mM HEPES, and protease inhibitors (Roche Complete Protease Inhibitors Cocktail, 1 tablet/50 ml), pH = 7. 5). The remainder of cell (cytosols) was processed separately. First, the nuclei were separated by centrifuge for 10 min at 700×g at 4°C. Plasma membranes and cytosols were solubilized in 35 μl running buffer (2% SDS, 10% glycerol, 5% β-mercaptoethanol, 0. 05% Bromophenol Blue, 62. 5 mM Tris-HCl pH 6. 8) and heated to 65°C for 5 min. Samples were electrophoresed on 12% polyacrylamide-SDS gel, and transferred to nitrocellulose membranes for Western blotting with the various antisera. The signals were visualized using the SuperSignal kit (Thermo) and quantitated using ImageJ software (National Institutes of Health, USA). In general, the macroscopic GIRK1/2 current, I, can be calculated utilizing a modified (Eq 1): I= Isingle·Po·N/fsc (5) where Isingle is a unitary current and N is the number of channels [1]. fcs is a solution conversion factor between solutions used for whole-cell (24 mM K+) and in cell-attached patches (144 mM K+). fsc was estimated as 4. 63±0. 26 (n = 6) by measuring GIRK currents in the same oocytes in the two solutions, in whole-cell configuration. For the graded contribution model, we define Po as: Po=Po, max⋅∑14fp, x⋅ϕx (6) where Po, max is the maximal open probability, fp, x is the fraction of Po contributed by x Gβγ-occupied channel (x is an integer between 1 and 4), and the ϕx is the fraction of channels in the x Gβγ-occupied state and can be calculated according to: ϕx=[Cx]Ctotal (7) where [Cx] is the concentration of channels in x Gβγ-occupied state and Ctotal is the total channel concentration in membrane. For the concerted model, Eq 6 is reduced to Po=Po, max⋅ϕ4 where ϕ4 is the fraction of channels with four bound molecules of Gβγ [58]. For the graded contribution model, we calculated values of fp, x based on data described by Ivanova-Nikolova et al. (1998) [13] rendering ~ 0. 01,0. 06,0. 26 and 1 values corresponding to 1–4 Gβγ occupied states. Based on mass-action law and Fig 2B, GIRK1/2 channel activity can be described by the following system of eqs (8–12): [C0]⋅[Gβγ]=14KD⋅[C1] (8) [C1]⋅[Gβγ]=23KD⋅[C2] (9) [C2]⋅[Gβγ]=32KD⋅[C3] (10) [C3]⋅[Gβγ]=4⋅KD⋅[C4] (11) [C0]+[C1]+[C2]+[C3]+[C4]=Ctotal (12) where the KD is a dissociation constant of Gβγ and GIRK1/2. For our simulations we routinely used KD = 50 nM as measured in direct biochemical experiments [90], but a range of other values has also been tested as explained in the Results. For both concerted and graded contribution models, G protein dissociation reaction required for modeling of GIRK1/2 basal activity according to the schemes described in Fig 2B is formulated as: kon⋅[Gβγ]⋅[GαGDP]=koff⋅[GαGDPGβγ] (13) Gβγtotal=[Gβγ]+[GαGDPβγ]+[C1]+2⋅[C2]+3⋅[C3]+4⋅[C4] (14) Gαtotal=[GαGDP]+[GαGDPGβγ] (15) where Gβγtotal and Gαtotal are total concentrations of corresponding subunits available for interaction with the channel and kon and koff are association and dissociation constants of Gi protein subunits (0. 7•106 M-1 s-1 and 0. 0013 s-1, respectively [56]). Eqs 8–12, combined with Eqs (13–15) is the most general form of description of a system containing GIRK channel and G proteins. For simulation of agonist-evoked activity, with saturating doses of both GPCR and agonist, we assumed a complete dissociation of G- protein heterotrimer [63,64] and thus Eq 14 is changed to: Gβγtotal= [Gβγ]+[C1]+2[C2]+3[C3]+4[C4] (16) Simulations were performed using Matlab and Berkeley Madonna software. Steady-state simulations used in model development and application, as well as in most Figures, were done utilizing Matlab 6. 5 function “solve” which is a part of Symbolic Math Toolbox. This function first looks for analytical solution, and if the former is absent, switches to numerical iterative algorithm (“trust region algorithm”, “quasi-Newton algorithm”). MATLAB routines for the calculation of Gβγ and Gα available for GIRK with the graded contribution model are shown in Supplemental Methods (S1 Text). In several cases we tested a range of arguments to produce continuous curves range of changes in GIRK1/2 currents or their ratios (Ra, Rβγ) (Fig 6C, S4, S7 and S8 Figs). Here, in order to reduce calculation time, we utilized Berkley Madonna software which implements 4th order Runge-Kutta method for numerical solution of differential equations (see Supplemental Methods, S1 Text). To assure lack of inconsistencies in calculation, we have compared the Matlab and Berkeley Madonna calculation results for a large number of cases and always obtained the same numbers. Imaging data on protein expression have been normalized as described previously [125]. Fluorescence intensity in each oocyte or giant membrane was calculated relative to the average signal in the oocytes of the control group of the same experiment. This procedure yields average normalized intensity as well statistical variability (e. g. SEM) in all treatment groups as well as in the control group. Statistical analysis was performed with SigmaPlot 11 (Systat Software Inc. , San Jose, CA, USA). If the data passed the Shapiro-Wilk normality test and the equal variance test, two-group comparisons were performed using t-test. If not, we performed the Mann-Whitney Rank Sum Test. Multiple group comparison was done with one-way ANOVA if the data were normally distributed. ANOVA on ranks was performed whenever the data did not distribute normally. Tukey’s post-hoc test was performed for normally distributed data and Dunn’s post-hoc test otherwise. Unless specified otherwise, the data in the graphs is presented as mean ± SEM. Correlation between two parameters (such as basal current and Ra) was tested using the Spearman correlation test by running this test on raw data using the statistical module of SigmaPlot 11.
summary: Many neurotransmitters and hormones inhibit the electric activity of excitable cells (such as cardiac cells and neurons) by activating a K+ channel, GIRK (G protein-gated Inwardly Rectifying K+ channel). GIRK channels also possess constitutive “basal” activity which contributes to regulation of neuronal and cardiac excitability and certain disorders, but the mechanism of this activity and its interrelation with the neurotransmitter-evoked activity are poorly understood. In this work we show that key features of basal and neurotransmitter-evoked activities are similar in cultured hippocampal neurons and in two model systems (mammalian HEK293 cells and Xenopus oocytes). Using experimental data of the neuronal GIRK1/2 channel function upon changes in GIRK and G protein concentrations, we constructed a mathematical model that quantitatively accounts for basal and evoked activity, and for the inverse correlation between the two. Our analysis suggests a novel and unexpected mechanism of interaction of GIRK1/2 with the G protein subunits, where the tetrameric GIRK channel can assemble with 4 molecules of the Gβγ subunits but only 2 molecules of Gα. GIRK is a prototypical effector of Gβγ, and the unequal stoichiometry of interaction with G protein subunits may have general implications for G protein signaling.
title: A Quantitative Model of the GIRK1/2 Channel Reveals That Its Basal and Evoked Activities Are Controlled by Unequal Stoichiometry of Gα and Gβγ