G Protein Activation without a GEF in the Plant Kingdom


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.


Published in the journal: . PLoS Genet 8(6): e32767. doi:10.1371/journal.pgen.1002756
Category: Research Article
doi: 10.1371/journal.pgen.1002756

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.

Phylogeny of G protein components in the plant kingdom.
Fig. 1. Phylogeny of G protein components in the plant kingdom.
(A) Conservation of G protein components in the plant kingdom. *, Genomes have not been sequenced; N/I, Not identified; +/−, 7TM-RGS gene is identified in a date palm, but not in the grass-family. (B–D) Phylogram of the consensus maximum likelihood (ML) tree as determined for Gα, Gβ and multi transmembrane RGS protein sequences. The trees were rooted with Homo sapiensi1, Gαq, Gβ1, Gβ5 or S. moellendorffii 7TM-RGS genes. Bootstrap values above 40 are shown near each branch. White, gray or black circles indicate that the branch was supported more than 90%, 80% or 70%, respectively. Species and gene names are shown in blue, green, yellow or purple colors, indicating eudicots, monocots, gymnosperms or non-vascular plants, respectively. Scale bar represents 0.1 substitutions per sites. See Figure S1, S2, S3 for aligned sequences used for creating the trees.

Tab. 1. G-protein components found in genome or EST databases of the land plant kingdom.
G-protein components found in genome or EST databases of the land plant kingdom.
Homologous genes of A.thaliana AtGPA1, AGB1, AGG1, AGG2, AGG3 and AtRGS1 were identified with translated protein BLAST against plant genome database through Phytozome v7.0 (released on Apr/8/2011; www.phytozome.net), PlantGDB database (http://www.plantgdb.org/), EST database registered in NCBI (www.ncbi.nlm.nih.gov), and EST database of M. polymorpha (More than 2 million sequence tags, http://Marchantia.pmb.lif.kyoto-u.ac.jp). (1) All RGS genes were predicted to have 7- or 5-transmembrane domain, except with non-transmembrane S. italica RGS. (2) DNA sequences with ∼28,000 gene predictions of P. dactylifera were downloaded from Weill Cornell Medical College in Qatar (http://qatar-weill.cornell.edu/research/datepalmGenome/download.html). G protein components were found with local BLAST search. (3) Sequence of G protein components was collected from EST database. Note that there might be unidentified G protein components in their species.

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 (Kobs = 5.80 min−1). 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, Kobs = 0.92 min−1). Our value is similar to that published in other studies [4], [5], [12]. Likewise, Gα from pine (PtGα1, Kobs = 6.85 min−1) and liverwort (MpGα1, Kobs = 1.84 min−1) 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).

Activation of plant Gα subunits.
Fig. 2. Activation of plant Gα subunits.
(A, B) Time course of [35S]GTPγS binding to 1 µM Gα at 20°C. Data are presented as the mean ± SEM of more than triplicates. (C, D) Intrinsic tryptophan fluorescence of Gα was measured at room temperature. 5 µM GTPγS was added to 400 nM Gα in the cuvette at time 0. Data are mean ± SEM of duplicate samples.

Tab. 2. Rates of nucleotide exchange and GTP hydrolysis.
Rates of nucleotide exchange and GTP hydrolysis.
Nucleotide exchange was measured with intrinsic Trp fluorescence and with [35S]GTPγS binding. Rate of GTP hydrolysis was measured with single turn [γ-32P]GTP hydrolysis. The rate-limiting step was determined by comparing the rates of GTP binding and hydrolysis. The percentage of bound GTP was approximated by the following equation: GTP binding rate/(GTP binding rate + GTP hydrolysis rate). Note that GTP/GDP ratio and Gα interacting proteins will affect the state in vivo. Rates were measured at 20°C and are reported as min−1.

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 32PO4 from [γ-32P]GTP in single turnover GTPase experiments (Figure 3D and Table 2). AtGPA1 (Kcat = 0.047 min−1 [4], [5]) and OsRGA1 (Kcat = 0.052 min−1) had slow rates of GTP hydrolysis. In contrast, liverwort MpGα1 (Kcat = 0.87 min−1) 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.

Inactivation of plant Gα subunits.
Fig. 3. Inactivation of plant Gα subunits.
(A–C) The activation and inactivation rates of Gα were monitored at room temperature by measuring the intrinsic tryptophan fluorescence. 800 nM GTP or 5 µM GTPγS was added to 400 nM Gα in the cuvette at time 0. Data are mean ± SEM of duplicate samples. Note that the turnover rate of Gα depends on stoichimometry of the concentrations of active Gα protein and GTP. The specific activity of the Gα subunits estimated by [35S]GTPγS binding assay were as follows: AtGPA1, 0.46 mol/mol protein; RGA1, 0.55 mol/mol protein; MpGPA1, 0.69 mol/mol protein. (D) Time course of single turnover [γ32P]GTP hydrolysis by 800 nM Gα. The mean ± SEM of duplicate samples is shown.

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α12 and Gα13, and this class of Gα subunits has dedicated RGS Homology (RH) proteins of Rho-family GEFs. Gα12 and Gα13 are not substrates for RGS proteins, which are dedicated GAPs of Gαi and Gαq. A mutation of this lysine of human Gα13 abolishes interaction with the RGS domain of a Gα13 effector, leukemia-associated Rho GEF (LARG) [22].

Evolution of monocot Gα and monocot specific loss of a 7TM-RGS gene.
Fig. 4. Evolution of monocot Gα and monocot specific loss of a 7TM-RGS gene.
(A) Evolutionary process to delete a 7TM-RGS gene in monocots. 1. The angiosperm ancestor possessed a 7TM-RGS gene. 2. After separation from the palm tree lineage, the Gα in the grass lineage lost a functional and physical coupling with its partner RGS protein. This uncoupling occurred by a single amino acid mutation, thus releasing the RGS gene from evolutionary constraint. 3. RGS genes were gradually deleted from grass-family genomes. 4. The extant S.italica RGS gene may still be functional, despite the loss of the 7TM region. (B) Switch I region of Gα. Conserved residues are highlighted with orange. Substituted Asn residue in monocots (Asn195 of OsRGA1) is shown with green. See Table S1 for species names. (C) Predicted structures of OsRGA1 (green) overlayed on Rattus norvegicusi1 and RGS4 (cyan and yellow, PDB: 1AGR [21]). Residues in purple (Gαi1) or red (OsRGA1) are at the binding surface to RGS4. Asn195 of OsRGA1 and Thr182 of Gαi1, discussed in this paper, are illustrated by a stick model. (D) The binding surface between Gαi1 and RGS4. Thr182 of Gαi1 and the interacting residues of RGS4 are shown in a stick model.

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−28) 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 (KD = 17.4 nM), and OsRGA1 had a relatively lower affinity for AtRGS1 (KD = 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 (KD = 115 nM). Reciprocally, the N191T mutation increased OsRGA1 affinity for AtRGS1 by 12-fold (KD = 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−1) 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.

Intrinsic properties of Gα mutants.
Fig. 5. Intrinsic properties of Gα mutants.
(A, B) Time course of [35S]GTPγS binding to 1 µM Gα at 20°C. Data are presented as the mean ± SEM of triplicates. (C, D) Time course of single turnover [γ32P]GTP hydrolysis by 800 nM Gα. The mean ± SEM of duplicate samples is shown.

Affinity of plant Gα to AtRGS1 immobilized on the SPR biosensor.
Fig. 6. Affinity of plant Gα to AtRGS1 immobilized on the SPR biosensor.
Recombinant AtRGS1 (284–459 aa) was immobilized on sensor chip CM5. AlF4-bound Gα subunits or the Gβγ dimer control (A) flowed over the chip at seven different concentrations (6.25, 12.5, 25, 50, 100, 200, 400 nM). Kinetics determined with 1∶1 (Langmular) binding model is shown in Table 3. Wild type Arabidopsis Ga subunit (B), T194N mutant Arabidopsis Gα subunit (C), wild type rice Gα subunit (D), and 195T mutant rice Gα subunit (E).

GAP activity of RGS toward plant Gα subunits.
Fig. 7. GAP activity of RGS toward plant Gα subunits.
(A–C) Time course of steady-state [γ32P]GTP hydrolysis by 500 nM Gα in the presence or absence of 750 nM AtRGS1 were measured over time after incubation at 20°C. Rate of Pi production (mol/mol Gα protein) were shown. Data are mean ± SEM of duplicate samples. (D–H) Single turnover [γ32P]GTP hydrolysis by 500 nM Gα proteins with AtRGS1 (0 nM (black), 5 nM (purple), 12.5 nM (yellow), 50 nM (orange), 125 nM (blue), 500 nM (red) and 750 nM (green)). Data are mean +/− SEM for more than two individual experiments, except OsRGA1 with 50 nM RGS at 1 min and 5 min. Dose dependency of single time point (2 min) was shown in (D).

Tab. 3. Affinity of plant Gα to AtRGS1 immobilized on the SPR biosensor.
Affinity of plant Gα to AtRGS1 immobilized on the SPR biosensor.
Recombinant AtRGS1 (284–459 aa) was immobilized on sensor chip CM5. Active form of Gα subunits was injected at seven different concentrations (6.25, 12.5, 25, 50, 100, 200, 400 nM). Kinetics shown in table was determined with 1∶1 (Langmular) binding model. ka, association rate constant; kd, dissociation rate constant; KD = kd/ka, equilibrium constant.

Tab. 4. Rates of steady state GTP hydrolysis.
Rates of steady state GTP hydrolysis.
Steady state GTP hydrolysis [min−1]. Steady-state [γ-32P]GTP hydrolysis was measured as described in Figure 7 at 20°C and reported as min−1.

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.

Model of G protein activation in the plant kingdom.
Fig. 8. Model of G protein activation in the plant kingdom.
Slow rate of GDP release and GTP hydrolysis is indicated by a thin arrow. A rapid rate is indicated by a thick arrow. In animals, a rate of GDP release from Gα is much slower than that of GTP hydrolysis. Thus, acceleration of the GDP release by GPCR changes the G protein from inactive to active. In eudicots and monocots, GDP release is rapid, and GTP hydrolysis is much slower than the GDP release. Thus, G protein can self-activate without the aid of a GPCR or other GEF. Instead, the eudicot G protein is regulated by a 7TM-RGS protein, which constitutively promotes GTP hydrolysis step on plasma membrane. However, some monocot genomes lack the 7TM-RGS gene, thus some monocot G protein must use an unknown mechanism to regulate activation. In addition, a 7TM-RGS gene is not expressed in a liverwort. However, a liverwort G protein has a rapid rate of both GDP release and GTP hydrolysis, which is likely to compensate for the loss of the 7TM-RGS gene.

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 sapiensi1, Gαq, Gβ1 and Gβ5 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 MgCl2, 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 MgCl2 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 [35S]GTPγS. The rate of GTP hydrolysis was determined with [γ32P]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 MgCl2) was mixed with an equal volume of TEDM buffer containing 5 µM [35S]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 MgCl2) 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 35S bound to the filter was quantified by scintillation counting.

For single-turnover GTP hydrolysis reactions, Gα subunit (800 nM) was preloaded with radioactive [γ-32P]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 MgCl2, 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 [γ-32P]GTP and proteins. The charcoal tubes were centrifuged, and amount of 32PO4 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 MgCl2, 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 MgCl2, 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 AlCl3 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 (ka) and dissociation (kd) 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


Zdroje

1. GilmanAG 1987 G proteins: transducers of receptor-generated signals. Annu Rev Biochem 56 615 649

2. SprangSR 1997 G protein mechanisms: insights from structural analysis. Annu Rev Biochem 66 639 678

3. SiderovskiDPWillardFS 2005 The GAPs, GEFs, and GDIs of heterotrimeric G-protein alpha subunits. Int J Biol Sci 1 51 66

4. JohnstonCATaylorJPGaoYKimpleAJGrigstonJC 2007 GTPase acceleration as the rate-limiting step in Arabidopsis G protein-coupled sugar signaling. Proc Natl Acad Sci U S A 104 17317 17322

5. JonesJCDuffyJWMachiusMTempleBRDohlmanHG 2011 The crystal structure of a self-activating G protein alpha subunit reveals its distinct mechanism of signal initiation. Sci Signal 4 ra8

6. ChenJGJonesAM 2004 AtRGS1 function in Arabidopsis thaliana. Methods Enzymol 389 338 350

7. ChenJGWillardFSHuangJLiangJChasseSA 2003 A seven-transmembrane RGS protein that modulates plant cell proliferation. Science 301 1728 1731

8. ChenJGGaoYJonesAM 2006 Differential roles of Arabidopsis heterotrimeric G-protein subunits in modulating cell division in roots. Plant Physiol 141 887 897

9. RichardsTACavalier-SmithT 2005 Myosin domain evolution and the primary divergence of eukaryotes. Nature 436 1113 1118

10. RogerAJSimpsonAG 2009 Evolution: revisiting the root of the eukaryote tree. Curr Biol 19 R165 167

11. SeoHSChoiCHLeeSYChoMJBahkJD 1997 Biochemical characteristics of a rice (Oryza sativa L., IR36) G-protein alpha-subunit expressed in Escherichia coli. Biochem J 324 Pt 1 273 281

12. IwasakiYKatoTKaidohTIshikawaAAsahiT 1997 Characterization of the putative alpha subunit of a heterotrimeric G protein in rice. Plant Mol Biol 34 563 572

13. BishtNCJezJMPandeyS 2011 An elaborate heterotrimeric G-protein family from soybean expands the diversity of plant G-protein networks. New Phytol 190 35 48

14. ChakravortyDTrusovYZhangWAcharyaBRSheahanMB 2011 An atypical heterotrimeric G-protein gamma-subunit is involved in guard cell K(+) -channel regulation and morphological development in Arabidopsis thaliana. Plant J 67 840 851

15. TrusovYRookesJETilbrookKChakravortyDMasonMG 2007 Heterotrimeric G protein gamma subunits provide functional selectivity in Gbetagamma dimer signaling in Arabidopsis. Plant Cell 19 1235 1250

16. GookinTEKimJAssmannSM 2008 Whole proteome identification of plant candidate G-protein coupled receptors in Arabidopsis, rice, and poplar: computational prediction and in-vivo protein coupling. Genome Biol 9 R120

17. MoriyamaENStropePKOpiyoSOChenZJonesAM 2006 Mining the Arabidopsis thaliana genome for highly-divergent seven transmembrane receptors. Genome Biol 7 R96

18. MooreMJSoltisPSBellCDBurleighJGSoltisDE 2010 Phylogenetic analysis of 83 plastid genes further resolves the early diversification of eudicots. Proc Natl Acad Sci U S A 107 4623 4628

19. HigashijimaTFergusonKMSmigelMDGilmanAG 1987 The effect of GTP and Mg2+ on the GTPase activity and the fluorescent properties of Go. J Biol Chem 262 757 761

20. HigashijimaTFergusonKMSternweisPCRossEMSmigelMD 1987 The effect of activating ligands on the intrinsic fluorescence of guanine nucleotide-binding regulatory proteins. J Biol Chem 262 752 756

21. TesmerJJBermanDMGilmanAGSprangSR 1997 Structure of RGS4 bound to AlF4–activated G(i alpha1): stabilization of the transition state for GTP hydrolysis. Cell 89 251 261

22. SuzukiNTsumotoKHajicekNDaigoKTokitaR 2009 Activation of leukemia-associated RhoGEF by Galpha13 with significant conformational rearrangements in the interface. J Biol Chem 284 5000 5009

23. WillardFSKimpleRJSiderovskiDP 2004 Return of the GDI: the GoLoco motif in cell division. Annu Rev Biochem 73 925 951

24. TakasakiJSaitoTTaniguchiMKawasakiTMoritaniY 2004 A novel Galphaq/11-selective inhibitor. J Biol Chem 279 47438 47445

25. JonesJCTempleBRJonesAMDohlmanHG 2011 Functional reconstitution of an atypical G protein heterotrimer and regulator of G protein signaling protein (RGS1) from Arabidopsis thaliana. J Biol Chem 286 13143 13150

26. BaucomRSEstillJCChaparroCUpshawNJogiA 2009 Exceptional diversity, non-random distribution, and rapid evolution of retroelements in the B73 maize genome. PLoS Genet 5 e1000732

27. BennetzenJL 2005 Transposable elements, gene creation and genome rearrangement in flowering plants. Curr Opin Genet Dev 15 621 627

28. UllahHChenJGYoungJCImKHSussmanMR 2001 Modulation of cell proliferation by heterotrimeric G protein in Arabidopsis. Science 292 2066 2069

29. UllahHChenJGTempleBBoyesDCAlonsoJM 2003 The beta-subunit of the Arabidopsis G protein negatively regulates auxin-induced cell division and affects multiple developmental processes. Plant Cell 15 393 409

30. FujisawaYKatoTOhkiSIshikawaAKitanoH 1999 Suppression of the heterotrimeric G protein causes abnormal morphology, including dwarfism, in rice. Proc Natl Acad Sci U S A 96 7575 7580

31. BersteinGBlankJLJhonDYExtonJHRheeSG 1992 Phospholipase C-beta 1 is a GTPase-activating protein for Gq/11, its physiologic regulator. Cell 70 411 418

32. TamuraKPetersonDPetersonNStecherGNeiM 2011 MEGA5: Molecular evolutionary genetic analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28 2731 2739

33. JonesDTTaylorWRThorntonJM 1992 The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci 8 275 282

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2012 Číslo 6

Nejčtenější v tomto čísle

Tomuto tématu se dále věnují…


Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Pacient na antikoagulační léčbě v akutní situaci
nový kurz
Autoři: MUDr. Jana Michalcová

Kopřivka a její terapie
Autoři: MUDr. Petra Brodská

Uroinfekce v primární péči
Autoři: MUDr. Marek Štefan

Roztroušená skleróza a plánování těhotenství
Autoři: MUDr. Radek Ampapa

Alergenová imunoterapie v léčbě inhalačních alergií
Autoři:

Všechny kurzy
Kurzy Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Nemáte účet?  Registrujte se

Zapomenuté heslo

Zadejte e-mailovou adresu se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se