G Protein Activation without a GEF in the Plant Kingdom

Summary

Animal heterotrimeric G proteins are activated by guanine nucleotide exchange factors (GEF), typically seven transmembrane receptors that trigger GDP release and subsequent GTP binding. In contrast, the Arabidopsis thaliana G protein (AtGPA1) rapidly activates itself without a GEF and is instead regulated by a seven transmembrane Regulator of G protein Signaling (7TM-RGS) protein that promotes GTP hydrolysis to reset the inactive (GDP-bound) state. It is not known if this unusual activation is a major and constraining part of the evolutionary history of G signaling in eukaryotes. In particular, it is not known if this is an ancestral form or if this mechanism is maintained, and therefore constrained, within the plant kingdom. To determine if this mode of signal regulation is conserved throughout the plant kingdom, we analyzed available plant genomes for G protein signaling components, and we purified individually the plant components encoded in an informative set of plant genomes in order to determine their activation properties in vitro. While the subunits of the heterotrimeric G protein complex are encoded in vascular plant genomes, the 7TM-RGS genes were lost in all investigated grasses. Despite the absence of a Gα-inactivating protein in grasses, all vascular plant Gα proteins examined rapidly released GDP without a receptor and slowly hydrolyzed GTP, indicating that these Gα are self-activating. We showed further that a single amino acid substitution found naturally in grass Gα proteins reduced the Gα-RGS interaction, and this amino acid substitution occurred before the loss of the RGS gene in the grass lineage. Like grasses, non-vascular plants also appear to lack RGS proteins. However, unlike grasses, one representative non-vascular plant Gα showed rapid GTP hydrolysis, likely compensating for the loss of the RGS gene. Our findings, the loss of a regulatory gene and the retention of the “self-activating” trait, indicate the existence of divergent Gα regulatory mechanisms in the plant kingdom. In the grasses, purifying selection on the regulatory gene was lost after the physical decoupling of the RGS protein and its cognate Gα partner. More broadly these findings show extreme divergence in Gα activation and regulation that played a critical role in the evolution of G protein signaling pathways.

Introduction

There are few well-understood examples of how signaling pathways evolved. In particular, it is not known how extant signaling molecules evolve characteristics including intrinsic activity, regulatory mechanisms and binding partners. Neutral selection theory proposes that genes released from constraints are gradually deleted from a genome. However, the processes whereby signaling genes are freed from constraints are not known and uninvestigated.

Heterotrimeric G proteins are well characterized molecular switches that are activated in response to extracellular stimuli [1], [2]. The G protein activation state is determined by the balance between rates of GDP-release (nucleotide exchange) and intrinsic GTP-hydrolysis by the Gα subunit of the heterotrimer [1], [2]. For all metazoan and yeast Gα proteins, GDP-release is slower than GTP-hydrolysis, and thus the G protein is predominantly GDP-bound in its resting state. However, both nucleotide exchange and hydrolysis are conditionally controlled by regulatory proteins in cells [3]. In animals and yeast, G protein-coupled receptors (GPCR) accelerate GDP release to favor the active GTP-bound state. Regulator of G Signaling (RGS) proteins accelerate GTP hydrolysis to favor the inactive GDP-bound state.

In contrast to this paradigm found in animals, Arabidopsis thaliana (Arabidopsis) Gα (AtGPA1) spontaneously self-activates without the aid of a GPCR or non-receptor GEF [4], [5]. Thus, in the absence of regulatory proteins, AtGPA1 would be predominantly GTP-bound [4], [5]. Instead, AtGPA1 inactivation is regulated in vivo by the single Arabidopsis RGS protein, AtRGS1 [4], [6], [7], which accelerates the intrinsically slow GTPase activity of AtGPA1 [4], [7]. AtRGS1 is the first identified protein with a domain architecture consisting of an N-terminal 7TM domain fused to an RGS domain [7], [8].

Plants and animals diverged from each other early in eukaryotic history. Based on recent evolutionary findings [9], [10], the plant kingdom is the most distinct group from animal lineages that are within Unikonts [10]. Our recent finding that G protein signaling in Arabidopsis differs greatly from that of animals raised the question of how these distinct signaling modules evolved in eukaryotes. Whether or not Bikonta other than Arabidopsis possess self-activating Gα proteins was unknown and, in cases, controversial. One group reported slow nucleotide exchange for the rice Gα [11], while another group reported relatively fast nucleotide exchange [12], albeit slower than the well-characterized Arabidopsis Gα protein [4], [5]. The Gα protein from Glycine max (soybean) may also possess relatively rapid GDP release [13] although there is no direct biochemical evidence supporting this idea.

Here we show that the plant kingdom employs G protein activation mechanisms distinct from those found in the animal kingdom. We analyzed plant genomes for G protein signaling components, and purified an informative subset of these components for biochemical analysis. We found that the trait of self-activating Gα was conserved throughout the plant kingdom. However, mechanisms that regulate G protein signal initiation differed throughout the plant kingdom, with some species lacking RGS proteins. We also provide evidence for the evolutionary route from one signal regulation mechanism to another. Specifically, we found in monocots that a single amino acid mutation in Gα disrupted the RGS-Gα interface and may have resulted in subsequent loss of the RGS genes from the genome. Collectively, these characteristics distinguish plant G protein signal regulation from the well-known paradigm from the animal kingdom. More broadly, this study illustrates the mechanism for how a strict functional pair (i.e. a signaling component and its regulator), commonly found in eukaryotes, was disrupted and resulted in rewiring of a cellular signaling network.

Results

G protein signaling components in plants

To identify signaling modules in the plant kingdom, homologous sequences of Gα, Gβ and Gγ genes were collected from genomic or expression sequence tags (EST) databases as described in Materials and Methods (Table 1 and Table S1, Figures S1, S2, S3). For reference, an evolutionary tree is provided in Figure 1A and includes the species described from here on. Typically vascular plants had one or two Gα genes, but G. max, a partially diploidized tetraploid, had four Gα genes. Physcomitrella patens (moss), a non-vascular land plant, lacked Gα homologs, although another non-vascular plant, Marchantia polymorpha (liverwort), possessed one Gα gene. One or two Gβ genes were encoded in all land plants analyzed, with the exception of soybean, which had four Gβ homologs. Likewise, Gγ genes resembling Arabidopsis Gγ genes [14], [15] were conserved in all land plants, with a few gene duplications. Notably, moss contained genes encoding the Gβγ dimer, but lacked a canonical Gα protein (Figure 1A). The moss genome encoded a gene (XP_001772174.1) homologous to Arabidopsis extra large GTP-binding protein (XLGA), although it should be noted that the moss gene lack a sequence for phosphate-binding loop (P-loop) and a glutamate residue in switch II region, each critical for G protein function.

Chlamydomonas reinhardtii and Volvox carteri, (unicellular and multicellular green algae, respectively) had no homologous genes for Gα, Gβ, Gγ and RGS, but a partial sequence of a Gα homologue was found in the EST database of Coleochaete scutata (a green alga, JG445935), a descendant of the most probable immediate ancestral group to land plants. These results suggest that non-vascular plant and chlorophyta lost some elements of the heterotrimeric G protein complex in their lineages.

Next, we searched for G protein regulatory elements. Previous analysis showed that plants lack canonical G protein coupled receptors [16], [17], and our analysis of new plant genomes/ESTs supported this finding. We discovered that genes encoding RGS proteins were not present in the most studied monocots, the cereals, even though RGS genes were present in all other vascular plants (eudicots, gymnosperms and a spikemoss). Although all grasses lacked a standard 7TM RGS protein, the grass, Setaria italica (foxtail millet) and the non-grass monocot Phoenix dactylifera (date palm) each possessed a gene that appears to encode an RGS protein. Unlike eudicots, however, the S. italica RGS lacked the transmembrane domains (Figure 1A and Table 1 and Table S1). Two eudicot RGS genes (Ricinus comunis and G. max RGS2) were predicted to have five transmembranes instead of seven transmembranes predicted for the founding member and prototype of the multi-transmembrane domain RGS family, AtRGS1. No RGS-homologous genes were found in non-vascular plants (liverwort and moss). Together, these results suggest that RGS proteins arose in an ancestor of vascular plants, but RGS genes were subsequently lost in many monocots.

We then phylogenetically analyzed the evolution of G protein signaling components. Generally, phylogenies of genes encoding plant Gα, Gβ and 7TM-RGS matched those generated with other genes used for phylogeny construction [18] (Figure 1B–1D). Near the end of angiosperm evolution, monocot and eudicot Gβ had approximately the same branch length from the common ancestor (Figure 1D). However, Gα evolution was accelerated in the monocot branch: the branch length of Poaceae (grass family) Gα from the common ancestor with P. dactylifera was almost twice as long as that of date palm Gα (Figure 1C). We hypothesize that this accelerated evolution of monocot Gα subunits compensated for the loss of RGS genes and/or was the result of the loss (discussed below).

The “self-activating” trait of the Gα protein is conserved in the plant kingdom

We included representatives from a eudicot (A. thaliana AtGPA1), a grass (Oryza sativa OsRGA1), a gymnosperm (Pinus taeda PtGα1), and a nonvascular plant (M. polymorpha MpGα1). First, we characterized the nucleotide exchange rates of these proteins using the non-hydrolysable GTP analog, GTPγS (Figure 2B and Table 2). Consistent with previous results [4], [5], we found that AtGPA1 had fast nucleotide exchange (K_obs = 5.80 min⁻¹). In contrast to an early report that suggested OsRGA1 had slow nucleotide exchange [12], we found that OsRGA1 exchange nearly matched that of AtGPA1 (OsRGA1, K_obs = 0.92 min⁻¹). Our value is similar to that published in other studies [4], [5], [12]. Likewise, Gα from pine (PtGα1, K_obs = 6.85 min⁻¹) and liverwort (MpGα1, K_obs = 1.84 min⁻¹) also had fast nucleotide exchange. These nucleotide exchange rates were corroborated by measuring the rate of the activation-dependent change in intrinsic Gα fluorescence [19], [20] (Figure 2C and 2D). Together these data suggest that the trait of fast GDP release is conserved in land plants, and likely arose in a common ancestor of this super group (Table 2).

For a Gα protein to be called “self-activating,” it must release GDP/bind GTP faster than it hydrolyzes GTP. In other words, the Gα should accumulate in its active form without a regulatory protein. Thus, we measured the rate of Gα-GTP accumulation in the presence of hydrolysable GTP (Figure 3A–3C). In this reaction, activated Gα would only be observed if the rate of nucleotide exchange was faster than the rate of GTP hydrolysis (i.e. when the Gα protein is “self-activating”) [5]. All tested plant Gα subunits accumulated in the active state with GTP (Figure 3A–C). Grass Gα (OsRGA1) and eudicot Gα (AtGPA1) displayed sustained activation in the presence of GTP. However, the Gα from liverwort (MpGα1) quickly returned to the inactive form, even in the presence of a 10-fold molar excess of GTP (Figure 3C), suggesting that the liverwort Gα had a fast rate of GTP hydrolysis. To test this hypothesis, we directly measured the intrinsic rates of inactivation of the selected Gα subunits by quantifying release of ³²PO₄ from [γ-³²P]GTP in single turnover GTPase experiments (Figure 3D and Table 2). AtGPA1 (K_cat = 0.047 min⁻¹ [4], [5]) and OsRGA1 (K_cat = 0.052 min⁻¹) had slow rates of GTP hydrolysis. In contrast, liverwort MpGα1 (K_cat = 0.87 min⁻¹) had a 16-times faster GTP hydrolysis rate than AtGPA1 and OsRGA1, indicating that MpGα1 efficiently inactivates itself without an RGS protein, yet hydrolysis remains the rate-limiting step. Together, these results suggest that land plant Gα subunits are all “self-activating” due to rapid nucleotide exchange relative to GTP hydrolysis and that the controlled step for activating G signaling is at GTP hydrolysis. The ideal element for this control is a 7TM-RGS protein, represented by the prototype AtRGS1. However, the absence of 7TM-RGS proteins in grasses indicates an alternative regulatory mechanism must exist in this class.

Grass-specific loss of the 7TM-RGS gene was preceded by a single amino acid mutation on Gα

Under neutral selection, genes freed from evolutionary constraint are rapidly deleted from the genome. This implies that 7TM-RGS was released from the strict functional constraint with Gα early in grass family history. To determine how this release may have occurred, we modeled the putative RGS - grass Gα protein interaction interface (Figure 4C and 4D) and found that a threonine residue in switch I (Thr194 of AtGPA1) that is critical for interaction with RGS proteins [21] was changed to asparagine in most monocot Gα subunits (Asn195 of OsRGA1, Figure 4B). This threonine residue is conserved in the RGS-sensitive human Gα_i and Gα_q family (Figure 4B) and is located at the center of the interface with RGS protein [21]. The threonine residue is substituted to lysine in Gα₁₂ and Gα₁₃, and this class of Gα subunits has dedicated RGS Homology (RH) proteins of Rho-family GEFs. Gα₁₂ and Gα₁₃ are not substrates for RGS proteins, which are dedicated GAPs of Gα_i and Gα_q. A mutation of this lysine of human Gα₁₃ abolishes interaction with the RGS domain of a Gα₁₃ effector, leukemia-associated Rho GEF (LARG) [22].

Two monocots were atypical in that they retained RGS-encoding genes. Examination of Gα and RGS sequences from these monocots provided insight into how other monocots may have lost RGS genes. First, the RGS protein in the monocot P. dactylifera has the typical (i.e. the Arabidopsis prototype) 7TM-RGS topology and its gene transcription is supported by EST data (Figure S4). The P. dactylifera Gα has the threonine critical for RGS interaction (Figure 4A). Notably, the P. dactylifera gene is longer than the Arabidopsis AtRGS1 gene by 19 kb, primarily due to the dramatic expansion of intron between the 7TM and RGS domains (Figure S4A) suggesting that this part of the gene was subjected to DNA insertion, possibly through transposon activity. The monocot S. italica also encodes a single soluble RGS gene. However, the S. italica Gα has Asn instead of Thr in the switch I region (Figure 4A). These analyses indicate that grass Gα subunits lost the ability to couple with RGS, thus releasing the genetic linkage between the Gα and the RGS protein, although it is also possible that deletion of RGS genes in grasses preceded the Gα mutation.

To trace the evolutionary process leading to the deletion of the 7TM-RGS gene in grasses, we surveyed S. italica genomic sequences surrounding the coding region of the single RGS gene (SiPROV019851m), and we found a hypothetical gene upstream of the RGS gene (SiPROV032159m) with sequence homologous (E value = 1e⁻²⁸) to transmembrane helices 2 through 4 of AtRGS1 (Figure S5 and Table S2). We also found nine ESTs that were homologous to the RGS domain. However, we found no sequence homologous to the region upstream of the 7TM domain (Figure S6). Moreover, the ratio of change in synonomous vs. nonsynonous residues (dn/ds) in comparing the S. italica homologies to their orthologs in date palm and Arabidopsis were higher for the 7TM region than for the RGS region. This suggests that the S. italica RGS domain continued with a function that was under strong purifying selection while the S. italica 7TM domain, although under purifying selection for most of the last 100 million years, has been under neutral or diversifying selection for the last few thousands or millions of years (Table S3).

Closer examination of the assembled sequence provided a clue to the partial gene loss. In the Setaria RGS region, we found two transposons inserted between the conserved and transcribed RGS domain and the apparent 7TM domain (Figure S6). Lack of EST support, suggest that the 7TM domain became a pseudo gene. One insertion is a partial sequence of a LINE transposon, likely resulting from deletion after insertion because the polyA and target site duplication (TSD) are missing. The second insertion was of a previously unknown long terminal repeat (LTR) retrotransposon that we named Alubu. Unequal homologous recombination subsequently converted this insertion into a solo LTR with intact TSD. There are 2 additional Alubu solo LTRs (but no complete elements) in the current Setaria sequence assembly (phytozome 7 http://www.phytozome.net/).

To identify other possible remnants of the 7TM-RGS gene in other grasses that lack RGS genes, we performed a tBLASTx search using the genomic sequences of the S. italica RGS-homologous region (segment 13, bases 1356001–1363646) against other monocot genomes. No homologous sequence of S. italica RGS was found in the other grasses analyzed, although the possibility of highly divergent RGS genes in plant genomes is not excluded.

These results indicate that a vascular plant ancestor had the 7TM-RGS gene. Furthermore, these analyses suggest that grasses gradually lost the RGS gene once it was uncoupled from the Gα protein by mutation of the RGS-Gα interaction interface. More broadly, these analyses point to the mechanism whereby a single amino acid substitution can lead to rewiring of a new signaling network In this case, the mutation led to neutral selection and loss of a regulatory element from the signaling pathway.

Physical uncoupling between RGS and monocot Gα

Our bioinformatics analyses suggested that the single amino acid substitution in the Gα protein-RGS interface was sufficient to release the RGS protein from evolutionary linkage to the Gα protein. To test this hypothesis experimentally, we substituted the threonine with an asparagine in the extant Arabidopsis Gα protein (AtGPA1-T194N) to recapitulate the monocot RGS interaction interface. We also made the comparable reverse substitution in OsRGA1, a representative monocot Gα protein. These mutations did not affect intrinsic nucleotide exchange and GTP hydrolysis rates (Figure 5A–5D). Next, we quantified interaction between these Gα proteins and the RGS protein from Arabidopsis (Figure 6A–6E and Table 3). As shown by surface plasmon resonance (SPR) analysis, AtGPA1 had high affinity for AtRGS1 (K_D = 17.4 nM), and OsRGA1 had a relatively lower affinity for AtRGS1 (K_D = 56.7 nM). We next tested two mutated Gα subunits, AtGPA1-T194N and OsRGA1-N195T. Although these mutations did not affect intrinsic nucleotide exchange and GTP hydrolysis rates (Figure 5A–5D), the T194N mutation reduced AtGPA1 affinity for AtRGS1 by 7-fold (K_D = 115 nM). Reciprocally, the N191T mutation increased OsRGA1 affinity for AtRGS1 by 12-fold (K_D = 4.83 nM). As a second measure of RGS-Gα interaction, we quantified GTPase acceleration by AtRGS1 in a steady-state GTP hydrolysis experiment (Figure 7A–7C and Table 4). Consistent with the affinities from SPR analysis, the T194N mutation of AtGPA1 reduced the GTPase acceleration by AtRGS1, and the N195T mutation of OsRGA1 increased GTPase acceleration by AtRGS1. This change in RGS1 sensitivity conferred by single amino acid substitution was confirmed using enzyme titration assays (Figure 7D–7H). Notably, intrinsic GTP hydrolysis by liverwort MpGα1 was fast (1.1±0.1 min⁻¹) and was not further stimulated by AtRGS1. These results suggest evolution of distinct regulatory systems of plant G proteins in the eudicots, grasses and liverworts.

Collectively, our phylogenetic and biochemical analyses suggest that the grass Gα lost the ability to interact with the regulatory molecule early in the evolutional lineage by substituting one critical residue (Figure 4A). The substitution of threonine to asparagine likely occurred in grasses before the loss of the RGS protein, a stage represented in S. italica, which contains the asparagine substitution, yet still encodes a remnant trace of the 7TM-RGS gene. This suggests that the physical uncoupling of Gα with RGS by single amino acid mutation broke the signaling pathway linkage permitting the subsequent deletion of RGS genes in grasses (Figure 1B).

Discussion

Evolution of “self-activating” Gα and 7TM-RGS

GDP release and GTP hydrolysis by Gα proteins are balanced to establish the steady-state level of the activated Gα subunit of the heterotrimeric G protein complex. In animals, G protein coupled receptors alter this balance to favor the GTP-bound state and relay signals from the outside of the cell to the inside of the cell. Likewise, RGS proteins accelerate GTP hydrolysis to favor the GDP-bound state and terminate intracellular signaling. Our recent discovery that these reactions are differently balanced in animals and Arabidopsis prompted us to examine divergence throughout the lineage and evolution of the G proteins and 7TM-RGS proteins within the plant kingdom. To complement our phylogenetic analyses of plant signaling components, we purified an informative set of plant Gα proteins that spanned the plant kingdom (Figure 8), and also investigated an amino acid substitution that was deduced to have occurred in the grass ancestral Gα protein. All tested Gα subunits were able to release GDP quickly without any other regulatory protein such as a GPCR or other guanine nucleotide exchange factor (Table 2) (i.e. they were “self-activating”). This finding is consistent with the fact that no unequivocal homologous GPCR gene has been characterized in the plant kingdom [16], [17]. These results provide powerful evidence that plant G proteins use different regulatory mechanisms than vertebrates to activate and terminate G protein signaling.

We previously proposed that AtRGS1 functions as a sugar receptor GAP of AtGPA1, operating by a sugar-dependent GAP activity [6], [7]. Here we found that liverworts, representing non-vascular plants, lack RGS genes altogether. In compensation, the liverwort Gα hydrolyzes GTP to GDP quickly without the aid of an RGS protein (Figure 7C and Table 3). The rates of liverwort GDP-release and GTP-hydrolysis were each fast and similar in value (Table 2), suggesting that liverwort Gα activity is equally balanced between the two reactions. Thus, liverwort Gα protein is likely regulated by other proteins yet to be identified. In contrast to liverwort, OsRGA1 (representing monocots), shared similar intrinsic activation/inactivation properties with AtGPA1 (representing eudicots). This self-activating property of OsRGA1 was surprising given that all studied grass genomes lost the standard 7TM RGS gene. As with the nonvascular plants, this finding points toward alternative regulatory mechanisms in grasses that were not identified based on homology to known G protein regulators from animals. These results indicate that plant G proteins use at least three different regulatory mechanisms, not only to activate, but also to terminate G protein signaling.

In addition to GEFs, mammalian G proteins are also regulated by GDP dissociation inhibitor proteins (GDI), which inhibit GDP release from the G protein and stabilize the GDP-bound state [23]. Since all plant Gα subunits spontaneously release GDP (Table 2), and some lack RGS proteins, plant G proteins are likely regulated by molecules having GDI activity. While several proteins and chemicals have GDI activity [23], [24], to date no GDI has been found in the plant kingdom other than the Gβγ dimer, AGB1/AGG1, although this was shown to be insufficient to maintain Gα in the inactive state [25].

Our analyses also identified 7TM RGS gene loss in progress in the S. italica genome. It is not possible to determine whether the insertion of transposon-like elements found in the S. italica RGS gene actually caused loss of the 7TM domain function, or whether this functional loss predated the insertion events. Transposable elements are the most abundant DNA sequences coded in plant genomes, and confer rapid rearrangement of plant genome structure [26], [27]. It is interesting, however, that the S. italica RGS domain continues to be expressed and under purifying selection. The S. italica RGS protein could be specifically coupled with the S. italica Gα with the Asn substitution. However, whether it is still involved in G protein signaling, without the need for the 7TM domain, is not known.

The grass genomes that have been extensively sequenced are dominated by species that have been cultivated for centuries, including rice, sorghum, maize, barley, wheat and S. italica (aka, foxtail millet). However, the observed amino acid change in the Gα protein and the loss of a normal 7TM RGS is not an outcome of domestication per se, as we see the same Gα protein and 7TM genome structure in S. viridis (the wild ancestor of S. italica) (unpub. obs.). We have shown that deletion of AtRGS1 from Arabidopsis results in increased cell growth and proliferation [7], [8], [28], [29]. Also, Gα mutants in Arabidopsis and rice show defects in their development [28], [30]. Our analyses raise the intriguing possibility that G protein signaling regulates growth and development with different regulatory mechanisms in grasses and eudicots. Regulators other than 7TM-RGS await discovery, or the grass family could have divergent RH proteins not identified by BLAST search. In mammals, the Gs-class of α subunits lacks a known RGS protein. For the Gq-class of Gα subunits, phospholipase Cβ functions as a GAP of Gαq [31]. Therefore, it follows that G protein activity in grasses or the other plants may be regulated by divergent effecters or the other binding proteins yet to be identified in plants.

Materials and Methods

Collection of Gα, Gβ, Gγ, and 7TM-RGS sequences

The sequences of G protein signaling components were found using BLASTp (E value<0.1) against protein database and the translated BLAST (tBLASTn, E value<0.1) against genomic DNA sequences registered in Phytozome v7.0 (released on Apr/8/2011; www.phytozome.net) by using A. thaliana genes as queries. Full-length or partial DNA sequences of Gα, Gβ, Gγ, and 7TM-RGS for Triticum aestivum, Hordeum vulgare, P. taeda, Picea glauca, and M.polymorpha were identified with tBLASTn in the nucleotide collection (nr/nt) database or the expressed sequence tags (EST) database at National Center for Biotechnology Information (NCBI) or the species-specific EST database (E value<10). The partial DNA sequences were combined to determine the full cDNA sequences. G protein components of P. dactylifera were found using assembled-gene sequences downloaded from Weill Cornell Medical College in Qatar (http://qatar-weill.cornell.edu/research/datepalmGenome/download.html). 7TM-RGS gene of P. dactylifera was assembled manually. Gα genes of P. taeda and P. glauca were cloned from the cDNA libraries and their sequences were determined, because information from the databases was insufficient to define the full length sequence. To screen all the possible RGS-like genes, P. dactylifera RGS and S. italica RGS and RGS domain sequences from H. sapiens RGS4, G protein-coupled receptor kinase,, LARG and sorting nexin 13 were also used as query sequences (E value<10).

Phylogenetic analysis

Phylogenetic trees were constructed with MEGA5.0 [32]. Full length Gα, Gβ and 7TM-RGS protein sequences were aligned with ClustalW using the following parameters; Gap opening penalty and gap extension penalty for initial pairwise alignment, 10 and 0.1; Gap opening penalty and gap extension penalty for multiple alignment, 10 and 0.2; Gonnet protein weight matrix; Residue-specific penalties, ON; Hydrophilic penalties, ON; Gap separation distance, 4; End gap separation, OFF; Use negative matrix, OFF. The maximum likelihood (ML) trees were made using the Complete-Deletion option of gaps and the JTT (Jones-Taylor-Thornton) substitution model [33] with gamma distributed rate variation, which was estimated as the best-fitting substitution model using MEGA5.0. The consensus phylogenetic trees were shown with the bootstrap values from 1000 repetitions. Homo sapiens Gα_i1, Gα_q, Gβ₁ and Gβ₅ were included as out groups.

Plasmids and proteins

cDNAs of P. taeda RGS and Gα were amplified from the cDNA library and cloned into pENTR-D/TOPO vector. cDNAs corresponding to O. sativa or M. polymorpha Gα were synthesized with optimization of codon usage for E.coli. AtGPA1-T194N and OsRGA1-N191T mutants were created by site directed mutagenesis. The Gα cDNAs were subcloned into pDEST17 (N-terminal 6×His). Recombinant His-tagged Gα proteins were expressed in ArcticExpress RP (Agilent Technologies) or Rosetta(DE3) (Novagen, used only for PtGα1) with 0.5 mM IPTG at 12°C, solubilized in buffer A (50 mM Tris-HCl (pH 7.5), 100 mM NaCl, 5 mM MgCl₂, 5 mM 2-Mercaptoethanol, 10 µM GDP, 1 mM PMSF and 1 µg/ml leupeptin) with 0.25 mg/ml lysozyme and 0.2% NP-40, captured from the soluble fraction with TALON Metal Affinity Resin (Clonetech), washed with buffer A containing 500 mM NaCl and 0.1% sodium cholate, and eluted with buffer A including 500 mM imidazole. 5 mM imidazole was added in crude extracts to reduce nonspecific binding. The purified proteins were dialyzed with buffer A and stored in 20% glycerol at −80°C. Recombinant His-AtRGS1 (284–459 aa) protein was prepared with the same method, except that MgCl₂ and GDP were removed from buffer A.

Gα activity

The rate of GTPγS binding was determined indirectly using intrinsic tryptophan fluorescence of Gα [20] and directly with [³⁵S]GTPγS. The rate of GTP hydrolysis was determined with [γ³²P]GTP.

For GTPγS binding, GDP-loaded Gα (1 µM) in TEDM buffer (50 mM Tris-HCl (pH 7.0), 1 mM EDTA, 1 mM DTT and 5 mM MgCl₂) was mixed with an equal volume of TEDM buffer containing 5 µM [³⁵S]GTPγS (about 5000 cpm/pmol) to start the binding reaction. At a given time points, 100 µl aliquots were quenched in 1 ml of ice-cold wash buffer (20 mM Tris-HCl (pH 7.5), 100 mM NaCl and 25 mM MgCl₂) containing 50 µM GTP and immediately vacuum-filtered onto nitrocellulose. Filters were quickly washed three times with 3 ml of ice-cold wash buffer. The total amount of ³⁵S bound to the filter was quantified by scintillation counting.

For single-turnover GTP hydrolysis reactions, Gα subunit (800 nM) was preloaded with radioactive [γ-³²P]GTP in TEDL (50 mM Tris-HCl (pH 7.5), 10 mM EDTA, 1 mM DTT, and 0.05% lubrol-PX) for 30 min on ice. The hydrolysis reaction was then started by adding 450 µl of TMDL+GTPγS (50 mM Tris-HCl (pH 7.5), 40 mM MgCl₂, 1 mM DTT, 0.05% lubrol-PX, and 400 µM GTPγS) into 1.2 ml of preloaded Gα. At a given time point, duplicate 100 µl aliquots were taken into 1 ml of charcoal (25% (w/v) in 50 mM phosphoric acid (pH 2.0)) to remove non-hydrolyzed [γ-³²P]GTP and proteins. The charcoal tubes were centrifuged, and amount of ³²PO₄ hydrolyzed was measured by scintillation counting of the centrifuged supernatants.

GTP or GTPγS binding with Trp fluorescence and steady state GTP hydrolysis were performed as described previously [5], [25]. Briefly, 400 nM Gα protein was incubated in a cuvette with 1 ml of TEMNG buffer (25 mM Tris-HCl (pH 8.0), 1 mM EDTA, 5 mM MgCl₂, 100 mM NaCl, and 5% glycerol). 800 nM GTP or 5 µM GTPγS was added to the cuvette and the change in the intrinsic fluorescence of Gα protein (excitation at 284 nm, emission at 340 nm) was monitored.

In vitro binding

Affinity between 2 different proteins was measured by Surface Plasmon resonance technology using BIAcore 2000 (GE Healthcare). His-tagged AtRGS1 (284–459aa) was immobilized on sensor chip CM5 with ammine coupling. Temperature, flow rate or running buffer were 25°C, 10 µl/min, or 10 mM Hepes, 150 mM NaCl, 3 mM EDTA, 0.005% Tween-20, 100 µM GDP and 10 mM MgCl₂, respectively. Seven different concentrations (6.25, 12.5, 25, 50, 100, 200 and 400 nM) of His-AtGPA1, AtGPA1-T194N, RGA1, RGA1-N195T or Gβγ (AGB1/AGG1) prepared in running buffer with 20 mM NaF and 100 µM AlCl₃ were flowed onto the sensor chip for 3 min. Dissociation was monitored for 5 min, and the sensor chip was washed with the same running buffer for 10 min at a flow rate of 30 µl/ml. The association (k_a) and dissociation (k_d) rate constants were obtained by fitting the original sensorgrams with a 1∶1 Langmuir binding model.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9