- Research article
Multiple domain insertions and losses in the evolution of the Rab prenylation complex
BMC Evolutionary Biologyvolume 7, Article number: 140 (2007)
Rab proteins are regulators of vesicular trafficking, requiring a lipid modification for proper function, prenylation of C-terminal cysteines. This is catalysed by a complex of a catalytic heterodimer (Rab Geranylgeranyl Transferase – RabGGTase) and an accessory protein (Rab Escort Protein. REP). Components of this complex display domain insertions relative to paralogous proteins. The function of these inserted domains is unclear.
We profiled the domain architecture of the components of the Rab prenylation complex in evolution. We identified the orthologues of the components of the Rab prenylation machinery in 43 organisms, representing the crown eukaryotic groups. We characterize in detail the domain structure of all these components and the phylogenetic relationships between the individual domains.
We found different domain insertions in different taxa, in α-subunits of RGGTase and REP. Our results suggest that there were multiple insertions, expansions and contractions in the evolution of this prenylation complex.
Protein prenyl transferases are a family of protein complexes that catalyze the lipid modification of proteins with isoprenoid groups. There are covalently attached to cysteine residues near or at the C-termini of intracellular proteins via tioether linkages (reviewed in [1, 2]). The family includes Protein Farnesyl Transferase (FTase), Geranylgeranyl Transferase I (GGTase I) and RabGeranylgeranyl Transferase (RabGTTase). FTase and GGTase I modify C-terminal cystein residues in the context of a CAAX motif with a 15 or 20 carbon isoprenoid, respectively, and are thus termed CAAX prenyltransferases (A stands for aliphatic residue, X for any residue). Substrates of the FTase include Ras family small GTPases, nuclear lamins, centromeric proteins among others, whereas substrates of the latter include the Rho family of small GTPases and heterotrimeric G protein γ subunits (reviewed in reference ). RabGGTase in contrast is specific to the Rab family of small GTPases, and catalyses the modification of these with two or sometimes one 20 carbon isoprenoid. It also differs from the CAXX prenyltransferases by the absolute requirement for an accessory protein, termed Rab Escort Protein (REP) for proper catalysis (reviewed in ). Protein prenylation affects proteins involved in a multitude of cellular processes, is involved in a variety of human diseases and therepaeutic approaches and this is extensively reviewed, for example in references [3–6].
Protein prenyl transferases are heterodimeric complexes of a α and a β chain, and in the case of RGGTase a third subunit, the REP protein(s). The genes coding for these enzymes have been cloned in a variety of species and shown to be essential for life [7–9]. FTase and GGTaseI share the α-subunit, but have distinct β-subunits, whereas the α- and β-subunits of RabGGTase are coded by distinct genes. The three-dimensional structure of the three enzymes has been solved, and revealed that the three enzymes are structurally homologous (see Figure 1). The α-subunit is a right-handed, crescent shaped, super-helix composed by 15 α-helices wrapped around the α-α barrel of the β-subunit [10, 11]. RabGGTase is a similar αβ heterodimer but the rat protein displays the insertion of two additional domains in the α-subunit relative to the other prenyltransferases, and to the yeast protein, a Ig-like domain and a Leucine Rich Repeat domain . In contrast with the CAXX prenyltransferases that recognize the C-terminal motif of the substrate proteins via the β, catalytic subunit, RabGGTase recognizes the substrate by interactions of the α-subunit with an extra protein, REP. The rat REP, like the α-subunit, displays an insertion relative to the yeast orthologue and the paralogous RabGDI [13, 14] (Figure 1). Protein prenyltransferases are believed to have evolved from an ancestral heterodimer, which by gene duplication gave rise to the current constellation of subunits .
The role of gene duplication in the evolution of protein complexes has been subject of recent interest [16–19]. Individual gene duplications can contribute to specialization of protein complexes, and can also accumulate to the point where two or more fully independent protein complexes exist [16, 19]. Less well understood is the role of domain insertions and deletions in the context of protein complexes. Protein domains are the basic evolutionary units in protein evolution [20, 21]. Proteins can gain, rearrange and loose domains in the course of evolution. Domain gains result mostly from fusions of independent genes that contain one or more domains and it is perhaps the most frequent mechanism of protein evolution after gene duplication . It can also be a consequence of duplication followed by in-frame fusion, resulting in tandem duplications. Less frequently, domain insertions may happen, situations where one or more domains are inserted into another domain, usually in loop regions [23, 24]. Domain losses are less well understood, but recent evidence suggests that they are most abundant at C-termini and due to the premature stop codons . Rearrangements or circular permutations can be achieved by complete duplication with in-frame fusion followed by partial deletion of domains at the termini, the most frequent mechanism , but also by independent fusion  and also by a mechanism termed "cut and paste" where a gene is fragmented, for example by the action of endonucleases, and then reassembles the resulting fragments in a different order .
Here we characterize the domain architectures of the components of the Rab prenylation complex, i.e. the heterodimer αβ of RGGTase and the accessory protein REP. We find that orthologous proteins display diverse domain architectures, suggestive of multiple independent events involving gains and losses of domains.
Results and discussion
α subunit – tetratricopeptide repeats
The rat α-subunit of RGGTase is very similar to the corresponding α-subunit of FTase, containing 15 α-helices arranged in a crescent-shaped, double layered right-handed superhelix, enveloping the β-subunit [11, 12]. Its structural architecture, together with statistically significant sequence similarity with the tetratricopeptide repeat (TPR) motif indicates that it belongs of the TPR superfamily . An individual TPR is a pair of anti-parallel α-helices, with consensus residues mediating the packing of these helices . In the rat RGGTase α-subunit these 7 TPR are helices 2 to 15. The crystal structure of the Rat RGGTase revealed two additional domains, a C-terminal Leucine Rich Repeat domain (LRR), and an Ig-like domain inserted between helices 11 and 12, i.e. between TPR 5 and 6. This is the exact same loop where in the FTase α-subunit there is a short 310 helix between α-helices 11 and 12 .
We investigated the complete and partial genome sequences of 43 organisms, covering the crown eukaryotic groups discicristates, excavates, alveolata, heterokonts, plants, amoebozoa and the opisthokonts (metazoa, fungi and choanoflagelates) . In all the organisms with a complete genome sequence available, we detect a putative RGGTase α-subunit (Figure S1). Thus we can conclude that the last eukaryotic common ancestor already possessed this subunit. FTase α-subunits were also found in the majority of the organisms studied (Figure S1), suggesting that the duplication that gave rise to the two distinct α-subunits preceded the split of the eukaryotic crown groups.
Detailed phylogenetic reconstruction of the evolutionary relationships between α-subunits based solely on the TPRs shows that the α-subunits of each enzyme are monophyletic (Figure S1). It is then most likely that there was one and only one gene duplication event that created the two paralogous subunits, and that this happened at the base of the eukaryotic tree. Furthermore, the fact that the secondary structure composed of 15 repeating helices with TPRs is conserved suggests that the ancestral of prenyltransferases already had this configuration.
We find that the secondary structure of this subunit varies in evolution, displaying multiple domain insertions in different organisms, which we characterize below.
α subunit – C2-like domain
Rat RGGTase a-subunit displays an inserted globular domain between helices 11 and 12 . The function of this domain is unclear, but it is clearly not involved in contacts with REP nor the Rab substrates . This globular domain is a β sandwich composed of eight strands in two sheets. It is a domain broadly related to the Ig fold, and in the two papers describing the structure of the complex it was termed an Ig-like domain. This is consistent with the CATH hierarchy, a fully automated classification of protein structures . Here we will consider instead the SCOP classification of evolutionary relationships between proteins structures, as it complements automated classification with manual curation [21, 31]. SCOP 1.71 classifies this domain as an independent superfamily, included in the C2 domain-like fold. From now on we will refer to this domain as the C2 domain-like.
The full SCOP hierarchy for this domain is shown in additional file 1 (Figure S2). This fold includes superfamilies such as the C2 domain. This is a domain that is found in multiple eukaryotic proteins and is involved in signaling, vesicular transport, modification of lipids, among other functions . C2 domains usually regulate their respective protein function by establishing Ca2+-dependent and Ca2+-independent phospholipids complexes. One class of C2 domains can bind Ca2+ without binding phospholipids . Recently, the C2 domain of PKCd was shown to mediate protein-protein interactions by binding directly to phosphotyrosine peptides in a sequence-specific manner . It is unclear if the RabGGTase C2-lke domain displays any of these functions. Since the role of Rab isoprenylation is to allow hydrophilic Rab proteins to associated with cellular membranes, it is plausible to think that such modification should occur in proximity to those membranes. If this C2 domain-like is a phospholipid-binding domain, then it could play a role in bringing the prenylation reaction next to membranes. Ultra-structural studies could be used to test this hypothesis.
It was previously observed that this domain was not present in the yeast orthologue of α-subunit RGGTase, but that it would be present in worm also . We analyzed in detail the sequences between helices 11 and 12 of RGGTase, where the Rat C2-like domain is found , and also where the FTase a-subunit displays an inserted 310 helix . We found that only a restricted number of branches on the eukaryotic tree display insertions in this region (Figure 2). These include metazoa, plants and alveolata. All other branches of the tree have a predicted secondary structure similar to that of Bet2 in S. cerevisae, which does no display any insertion between the two TPRs.
Next we investigated whether the insert we found on the three distinct branches of the tree of life are similar, and hence likely to represent one single insertion event in the ancestor of all RGGTases that was subsequently lost in other branches. Or instead whether they are distinct domains, resulting from independent insertion or expansion events. Our results favor the second hypothesis (Figure 2). We find that the inserted domains are similar within taxonomical groups, but different beyond recognition across taxa. By difference beyond recognition we mean BLAST sequence searches , Pfam  and Superfamily [37, 38] domain assignments and secondary structure predictions [39, 40]. Thus, in Metazoa, all Deuterostomes have a clearly defined C2-like domain, but the insects and the nematodes have a distinct insertion conserved solely in their taxonomical group. In the plants, angiosperms have a conserved domain in the same region of the α-subunit, but different beyond recognition from any other insert of the orthologous α-subunits. The same is true for alveolata. The size of the insertion is also conserved with taxonomic group but not across groups. For example, whereas deuterostomes have an insert of near 125 amino residues long, plants have a larger one, nearing 200 residues (Figure S3). Our results are thus compatible at least with up to five distinct insertion and/or expansion events in the same position of the α-subunit of RGGTase.
One exciting possibility is that all these insertions and or expansions represent the same function accomplished by different sequences, as this would expand our repertoire of sequence-function relationships.
Multiple sequence insertions/expansions in the same site suggest that this site is capable of accommodating structural variations more easily than others. Thus there seems to be a structural constraint in place. The fact that the paralogous α-subunit of FTase contains an inserted region in between these helices, in the form of a 310 helix adds support to this hypothesis. However, if the function of these distinct insertions and/or expansions is the same, then this could be the only place where this function is possible, and thus the recurrent use of the same site reflect functional rather than structural constraint. In the absence of information regarding the role played by these insertions/expansions it is not possible to resolve this question.
α subunit – LRR domain
At the C-terminus of rat RGGTase α-subunit there is a Leucine Rich Repeat domain (Figure 1). This domain in not present in FTase (Figure 1). LRRs belong to the SCOP fold of the same name, which groups proteins forming a right-handed β-α superhelix [21, 31]. It is formed by three superfamilies, the RNI-like, Outer arm dynein light chain and the L domain-like which includes the RGGTase LRR domain. LRRs are involved in a variety of biological processes, both in eukaryotes and prokaryotes. Their common role is the establishment of complexes with other proteins [41, 42].
We now investigate whether this LRR domain is a feature of all RGGTase α-subunits, which would suggest that it was present in the ancestral eukaryote, or if it is a recent acquisition restricted to a specific taxon or set of taxa. It is clear from the results shown in Figure 2 that the LRR domain is not universal, as we can only detect it in some animals, in angiosperms and in alveolata. This phylogenetic profile is consistent with two evolutionary scenarios – independent domain fusions or a single domain fusion at the base of the eukaryotic tree followed by a specific domain loss in multiple branches of the eukaryotic tree. Both scenarios seem equally unlikely, so we investigated this further using phylogenetic reconstruction based on the LRR domain sequences only. Our hypotheses is that there is enough phylogenetic signal in these sequences to solve this puzzle. In Figure 3 we show three phylogenies of the LRRs, including sequences from other superfamilies as reference, reconstructed by three distinct methods. Plant and vertebrate LRR are consistently monophyletic, suggesting a common origin. In contrast, the LRR of Ciona, of nematodes and of alveolata segregate with different reference sequences. This varies according to the method used to reconstruct phylogeny. It is thus impossible at this stage to resolve the question regarding the common or independent origin of the LRR sequences in the α-subunit. However, the recent observation that convergent evolution of domain architectures is very rare, with an estimated frequency of 0.4% to 4% , is more supportive of the first scenario.
β-subunit domain architecture
The β-subunit of rat RGGTase is a α-α barrel composed of 12 α helices. It is very similar to the α-α barrel in the β-subunit of FTase . The β-subunits of prenyltransferases are more conserved than the α subunits, but identification and classification of RGGTase β-subunit was simple using a combination of BLAST searches of sequences databases followed by phylogenetic analysis (Figure S4).
We discussed above that the α-subunit display multiple sequence insertions in different species. Below we will show that in some species REP also display inserted sequences. The β-subunit in contrast appears to have an invariable domain architecture throughout evolution. We investigated its structural relatives in order to gain insight whether this is due to structural constraints. In order to do so, we investigated the SCOP hierarchical classification of protein structures. Prenyl transferase's β-subunits belong to the "α-α toroid" fold (SCOP: 48207), Terpenoid cyclases/Protein prenyltransferases superfamily (SCOP: 48239). This fold consists of multihelical proteins displaying up to seven alpha-hairpins arranged in a closed circular array. Therpene synthases are classified into the same superfamily. These proteins are characterized by two α-α domains in the same peptide chain; the first is an α6-α6 barrel of two concentric rings, whereas the second is a barrel with 10 α helices, one 310 helix and at least two β-strands inserted in between helical elements . It appears then that this superfamily can have varying number of helical elements as well as accommodate extra structural elements. It is then plausible that the invariable nature of the domain architecture of β-subunits is due to functional rather than structural constraints.
In animals and yeast it is clear that RGGTase is composed of independent α and β-subunits. In contrast GGTase I and FTase have distinct β-subunits, but share a α-subunit. Thus the GGTase I and FTase are related by a single gene duplication. In contrast, RGGTase is separated from the other prenyltrasnsferases by two gene duplications. Although the number of duplication steps that separates the different enzymes is clear, the order of duplication is not. In other words, we don't know which subunits emerged first and which resulted from these by duplication. We find that all organisms with a complete genome sequence that we investigated display a β-subunit of RGGTase (Figure 2). Considering we also always find a RGGTase α-subunit in the same organisms, it follows that the ancestor of all eukaryotes already had a distinct and fully separated RGGTase.
Rab escort proteins (REP) belong to the same protein family as Rab GDP dissociation inhibitors (RabGDI). They are both classified in the SCOP hierarchy  as FAD/NAD(P)-binding domain fold and superfamily, which suggests a common ancestry. Their structure comprises two domains: domain I include the Rab binding platform, whereas domain II in REP mediates binding to the alpha-subunit of RGGTase [29, 45]. REP and RabGDI share conserved regions, termed SCRs (sequence conserved regions), which are highlighted in Figure 4 as brown boxes. Multiple sequence alignments of REPs and RabGDI reveal that mammalian REPs display an insertion between domain I and domain II, absent in RabGDI, which maps roughly to a sequence region delimited by the conserved regions SCR1B and SCR2 [13, 45]. This insert appears to be larger in vertebrate sequences than in S. cerevisae . The function of this insert is unclear, particularly at the light of the recent structure of REP1 in complex with RGGT, which shows that it is not involved in contacts with the RGGTase subunits nor with Rab substrates .
We identified REPs in all species studied here by searching GeneBank and the other genomic databases (see methods) with known REP sequences. We used phylogenetic reconstruction to classify the obtained sequences into the REP and RabGDI subfamilies (Figure 4B). We identified at least one distinct REP sequence in all species for which a complete genome sequence was available. Mammals have two paralogous REPS (REP-1 or Choroideremia and REP-2 or CHM-like). X. laevis also displays two paralogous REPs, but our phylogenetic analysis suggests that this is an independent and species-specific duplication (Figure 4b). It is interesting to note that in A. thaliana, the duplication of both enzyme subunits was not accompanied by REP duplication.
The most striking observation is that few branches of the tree of life are characterized by the presence of a larger insert region between SCR1B and SCR2 in REPs. Their size indicates that few independent taxa display inserts that are larger than those of yeast MRS6 (yeast REP) and RabGDIs (Figure 4C). At the sequence level there is no real conservation – the inserts are similar within taxa but very different across taxa. Systematic database searches using solely the insert regions from several specie can only find closely related REP sequences. For example, a BLAST search using the insert region of rat or human REP1 will find only deuterostomes REP inserts, but not plant inserts, and vice versa. They are also not similar to any other protein other than REPs, which means that we cannot use this approach to define hypothetical functions for this region.
In the crystal structure of rat REP1, this sequence insert corresponds to a region in the crystal with no clear electron density . Such regions are typically labeled natively unfolded [46, 47]. Natively unfolded proteins are involved in a variety of cellular functions, namely transcriptional and translational regulation, signaling and regulation of the self-assembly of large multi-subunit complexes such as the ribosome and the bacterial flagellum [47, 48]. They are also expected to be involved in a variety of human diseases . Although they can perform their function in the unfolded state, the majority of unfolded regions undergo a process termed induced-folding, in which upon binding to their physiological partners they undergo a transition to a structured form . We tested whether the different inserts were also natively unfolded, using the predictive algorithm Globplot2 . We observed that angiosperm insert sequences are predicted to include a disordered, region, so are all vertebrate sequences and nematodes. In contrast, insects never show a predicted disordered region within the insert region (Figure 4C and S5). Thus, the presence of disordered regions is not restricted to the larger inserts of vertebrates, nor does it seem to correlate with insert size.
In most species studied, we find a single copy of the α- and the β-subunit. In contrast we found two RGGTase α- and β-subunits in A. thaliana. We have previously observed that duplication of whole protein complexes is frequent . We argued for the prevalence of stepwise duplications leading to complete duplication of subunits of complexes, like what is observed with the adaptin tetrameric complexes . The duplication of the two subunits of RGGTase in A. thaliana however appears to have occurred simultaneously, as a result of a whole genome duplication (WGD) estimated to have happened around 38 million years ago . The paralogous pairs localize to distinct chromosomes, to duplicated segments that were mapped to that WGD (α-subunit: At4g2424490 (chr.4) and At5g41820 (chr.5); β-subunit: At3g12070 (chr.3) and At5g12210 (chr.5)).
The functional relevance of maintaining these two copies of RGGTase is unclear. It is clear from all the phylogenetic trees that the subunits are very closely related (e.g. Figure S1), suggesting some selective pressure. A. thaliana may require high levels of this enzyme and concerted gene duplication is one way of boosting the levels of a given gene product . Alternatively, the large number and diversity of Rab GTPases in A. thaliana [55, 56] may require specific regulation of Rab prenylation. This could be achieved by having multiple copies of the enzyme subunits under differential regulation or displaying distinct substrate specificities.
Animals also display an expansion of the Rab family. They don't have multiple copies of the enzyme subunits; they have instead two paralogous REPs that appear to form complexes with distinct substrates and have distinct specificities to the RGGTase . This lends some support to the second hypothesis. It would be thus interesting to investigate whether there is specificity in the pairing of the A. thaliana subunits, and whether the different versions of the enzyme have distinct substrate specificities. In contrast, some animals display two or more copies of the β-subunit for a single α-subunit (Figure 2). Since substrate-binding specificity is not defined by the β-subunit in RabGGTase, it seems plausible that expression levels are at the root of these duplications. Amoebozoa like D. discoideum, which represents another example of independent expansions of the Rab family, have in excess of 50 Rab proteins (unpublished observations). D. discoideum only has one copy of each component (Figure 2), indicating that large Rab families are compatible with a single copy of each RabGGTase component.
The duplication of components of RabGGTase appear thus to be restricted to specific branches of the tree of life (Figure 2). Our analysis also indicates that it is also restricted to RabGGTase. We did not observe duplications of FTase nor GGTase I subunits. This is despite the fact that their substrates show expansions akin to those of the Rab family. For example, the Ras and Rho families expanded from 3 and 6 members respectively in S. cerevisiae, to 22 and 34, respectively, in H. sapiens . It is possible that dosage balance may place a barrier to duplications of subunits [18, 59].
In summary, we identified the orthologues of the components of the Rab prenylation complex in 43 species – the αβ heterodimer RabGGTase, and the accessory proteins REP. We characterized their domain architectures and found that it varies considerably in evolution. The α-subunit can have distinct inserts in two positions, between the tetratricopeptide repeats 5 and 6 and at the C-terminal. These domain insertions are specific to RabGGTase, as the paralogous α-subunit of Farnesyl and Geranylgeranyl (type I) Transferases displays conserved domain architectures. The β subunit has a conserved domain architecture but the REP proteins also have variable inserts between SCR1B and SCR2.
We found the full constellation of protein prenyltransferases in the organisms we analyzed as found in model organisms [1, 2]. This suggests that the Last Common Eukaryotic Ancestor (L.E.C.A.) had the same constellation of enzymes and subunits. The three enzymes are related by gene duplication events, which suggests that there was an ancestral heterodimeric protein prenyl transferase at the base of the eukaryotic tree. The presence of structural homologues of α and β-subunits in all the branches of the tree of life indicates that the eukaryotic innovation required solely novel interactions, rather than novel folds. The precise order of duplication events is however unclear. FTase and GGTase I are separated by a single gene duplication event, whereas RGGTase is separated by two duplication events from either of the other enzymes. The simpler architecture of FTase and GGTase I, and their independence from further components for proper activity suggests that they preceded RGGTase. However, substrates of all three enzymes are ubiquitous in the eukaryotic tree, which further emphasizes that the L.E.C.A. is likely to have already the three enzymes.
In conclusion, the components of the Rab prenylation complex display varied domain architectures in evolution, which are more consistent with multiple independent events in the first insertion of the α-subunit, multiple losses in the second insert, and independent expansions in REP. There is no known function for the variable domains. There is no obvious correlation between Rab family size, and the presence or size of any of the RGGTase inserts, which suggests that these insertions are independent of Rab family expansion. Furthermore, the inserted domains do not parallel each other within or across subunits. This indicates that it is unlikely that they are involved in direct physical contacts, and that they are likely to have independent functions. The few inserted sequences that are similar to known sequences (C2-like and LRR) suggest that the role of these domains is regulatory, possibly involving lipid and protein binding. It is tempting to speculate that they are involved in taxon-specific regulatory interactions not yet described. One future avenue of research will be the identification of putative binding partners, for example by bioinformatic analysis, searching for genes with similar phylogenetic profiles as those of the inserted domains.
Finally, this case study suggests that domain gains and losses may be an important force driving the evolution and diversification of protein complexes. A future avenue of research is the quantification of this contribution.
The analyses were performed on sequenced eukaryotic genomes that were downloaded from NCBI and other eukaryotic genomes databases. A full listing of the organisms investigated is provided in additional file 1. It aims for maximum coverage of eukaryotic diversity within the species with a complete genome sequence. Sequence searches combined BLAST  searches using known RGGTase and REP sequences as query, as well as protein families defined in the Pfam  and superfamily [37, 38] databases. Multiple sequence alignments were performed with ClustalW 1.83 , pairwise alignments used the Smith-Waterman algorithm  (Water in EMBOSS ), in both cases with the a Blosum 62 matrix  and default GAP and extension penalties. Sequence manipulation was done with Jalview 2.1.1 . Domain assignments were done using Superfamily  and Pfam . Secondary structure predictions was performed using Jpred [39, 40]. Phylogenetic reconstruction was done using the Neighbor-Joining clustering algorithm as implemented in ClustalW 1.83 , using 1000 bootstraps, as well as the cladistic methods Maximum Likelihood in the Phylip 3.61 package (ProML)  (Jones-Taylor-Thorthon (JTT) matrix; 100 boostraps) and Bayesian method implemented in MrBayes v. 3.1.2 [68, 69] using Blosum62 as a fixed rate aa model, ran until average standard deviation of split frequencies was lower than 0.01. Trees were drawn using in FigTree v. 1.0 .
All sequences used in this study, as well as alignment and tree files can be found in the supplementary website .
Maurer-Stroh S, Washietl S, Eisenhaber F: Protein prenyltransferases. Genome Biol. 2003, 4 (4): 212-10.1186/gb-2003-4-4-212.
Casey PJ, Seabra MC: Protein prenyltransferases. J Biol Chem. 1996, 271 (10): 5289-5292. 10.1074/jbc.271.10.5289.
Lane KT, Beese LS: Thematic review series: lipid posttranslational modifications. Structural biology of protein farnesyltransferase and geranylgeranyltransferase type I. J Lipid Res. 2006, 47 (4): 681-699. 10.1194/jlr.R600002-JLR200.
Leung KF, Baron R, Seabra MC: Thematic review series: lipid posttranslational modifications. geranylgeranylation of Rab GTPases. J Lipid Res. 2006, 47 (3): 467-475. 10.1194/jlr.R500017-JLR200.
McTaggart SJ: Isoprenylated proteins. Cell Mol Life Sci. 2006, 63 (3): 255-267. 10.1007/s00018-005-5298-6.
Pereira-Leal JB, Hume AN, Seabra MC: Prenylation of Rab GTPases: molecular mechanisms and involvement in genetic disease. FEBS Lett. 2001, 498 (2-3): 197-200. 10.1016/S0014-5793(01)02483-8.
Jiang Y, Rossi G, Ferro-Novick S: Bet2p and Mad2p are components of a prenyltransferase that adds geranylgeranyl onto Ypt1p and Sec4p. Nature. 1993, 366 (6450): 84-86. 10.1038/366084a0.
Jiang Y, Ferro-Novick S: Identification of yeast component A: reconstitution of the geranylgeranyltransferase that modifies Ypt1p and Sec4p. Proc Natl Acad Sci U S A. 1994, 91 (10): 4377-4381. 10.1073/pnas.91.10.4377.
He B, Chen P, Chen SY, Vancura KL, Michaelis S, Powers S: RAM2, an essential gene of yeast, and RAM1 encode the two polypeptide components of the farnesyltransferase that prenylates a-factor and Ras proteins. Proc Natl Acad Sci U S A. 1991, 88 (24): 11373-11377. 10.1073/pnas.88.24.11373.
Taylor JS, Reid TS, Terry KL, Casey PJ, Beese LS: Structure of mammalian protein geranylgeranyltransferase type-I. Embo J. 2003, 22 (22): 5963-5974. 10.1093/emboj/cdg571.
Park HW, Boduluri SR, Moomaw JF, Casey PJ, Beese LS: Crystal structure of protein farnesyltransferase at 2.25 angstrom resolution. Science. 1997, 275 (5307): 1800-1804. 10.1126/science.275.5307.1800.
Zhang H, Seabra MC, Deisenhofer J: Crystal structure of Rab geranylgeranyltransferase at 2.0 A resolution. Structure. 2000, 8 (3): 241-251. 10.1016/S0969-2126(00)00102-7.
Alory C, Balch WE: Organization of the Rab-GDI/CHM superfamily: the functional basis for choroideremia disease. Traffic. 2001, 2 (8): 532-543. 10.1034/j.1600-0854.2001.20803.x.
Wu SK, Zeng K, Wilson IA, Balch WE: Structural insights into the function of the Rab GDI superfamily. Trends Biochem Sci. 1996, 21 (12): 472-476. 10.1016/S0968-0004(96)10062-1.
Zhang H, Grishin NV: The alpha-subunit of protein prenyltransferases is a member of the tetratricopeptide repeat family. Protein Sci. 1999, 8 (8): 1658-1667.
Pereira-Leal JB, Teichmann SA: Novel specificities emerge by stepwise duplication of functional modules. Genome Res. 2005, 15 (4): 552-559. 10.1101/gr.3102105.
Pereira-Leal JB, Levy ED, Kamp C, Teichmann SA: Evolution of protein complexes by duplication of homomeric interactions. Genome Biol. 2007, 8 (4): R51-10.1186/gb-2007-8-4-r51.
Papp B, Pal C, Hurst LD: Dosage sensitivity and the evolution of gene families in yeast. Nature. 2003, 424 (6945): 194-197. 10.1038/nature01771.
Li H, Pellegrini M, Eisenberg D: Detection of parallel functional modules by comparative analysis of genome sequences. Nat Biotechnol. 2005, 23 (2): 253-260. 10.1038/nbt1065.
Vogel C, Teichmann SA, Pereira-Leal J: The relationship between domain duplication and recombination. J Mol Biol. 2005, 346 (1): 355-365. 10.1016/j.jmb.2004.11.050.
Murzin AG, Brenner SE, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. J Mol Biol. 1995, 247 (4): 536-540. 10.1006/jmbi.1995.0159.
Kummerfeld SK, Teichmann SA: Relative rates of gene fusion and fission in multi-domain proteins. Trends Genet. 2005, 21 (1): 25-30. 10.1016/j.tig.2004.11.007.
Russell RB: Domain insertion. Protein Eng. 1994, 7 (12): 1407-1410. 10.1093/protein/7.12.1407.
Aroul-Selvam R, Hubbard T, Sasidharan R: Domain insertions in protein structures. J Mol Biol. 2004, 338 (4): 633-641. 10.1016/j.jmb.2004.03.039.
Weiner J, Beaussart F, Bornberg-Bauer E: Domain deletions and substitutions in the modular protein evolution. Febs J. 2006, 273 (9): 2037-2047. 10.1111/j.1742-4658.2006.05220.x.
Weiner J, Bornberg-Bauer E: Evolution of circular permutations in multidomain proteins. Mol Biol Evol. 2006, 23 (4): 734-743. 10.1093/molbev/msj091.
Bujnicki JM: Sequence permutations in the molecular evolution of DNA methyltransferases. BMC Evol Biol. 2002, 2: 3-10.1186/1471-2148-2-3.
Baldauf SL: The deep roots of eukaryotes. Science. 2003, 300 (5626): 1703-1706. 10.1126/science.1085544.
Pylypenko O, Rak A, Reents R, Niculae A, Sidorovitch V, Cioaca MD, Bessolitsyna E, Thoma NH, Waldmann H, Schlichting I, Goody RS, Alexandrov K: Structure of Rab escort protein-1 in complex with Rab geranylgeranyltransferase. Mol Cell. 2003, 11 (2): 483-494. 10.1016/S1097-2765(03)00044-3.
Pearl F, Todd A, Sillitoe I, Dibley M, Redfern O, Lewis T, Bennett C, Marsden R, Grant A, Lee D, Akpor A, Maibaum M, Harrison A, Dallman T, Reeves G, Diboun I, Addou S, Lise S, Johnston C, Sillero A, Thornton J, Orengo C: The CATH Domain Structure Database and related resources Gene3D and DHS provide comprehensive domain family information for genome analysis. Nucleic Acids Res. 2005, 33 (Database issue): D247-51. 10.1093/nar/gki024.
Andreeva A, Howorth D, Brenner SE, Hubbard TJ, Chothia C, Murzin AG: SCOP database in 2004: refinements integrate structure and sequence family data. Nucleic Acids Res. 2004, 32 (Database issue): D226-9. 10.1093/nar/gkh039.
Nalefski EA, Falke JJ: The C2 domain calcium-binding motif: structural and functional diversity. Protein Sci. 1996, 5 (12): 2375-2390.
Shin OH, Han W, Wang Y, Sudhof TC: Evolutionarily conserved multiple C2 domain proteins with two transmembrane regions (MCTPs) and unusual Ca2+ binding properties. J Biol Chem. 2005, 280 (2): 1641-1651. 10.1074/jbc.M407305200.
Benes CH, Wu N, Elia AE, Dharia T, Cantley LC, Soltoff SP: The C2 domain of PKCdelta is a phosphotyrosine binding domain. Cell. 2005, 121 (2): 271-280. 10.1016/j.cell.2005.02.019.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.
Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR: The Pfam protein families database. Nucleic Acids Res. 2004, 32 (Database issue): D138-41. 10.1093/nar/gkh121.
Wilson D, Madera M, Vogel C, Chothia C, Gough J: The SUPERFAMILY database in 2007: families and functions. Nucleic Acids Res. 2007, 35 (Database issue): D308-13. 10.1093/nar/gkl910.
Gough J, Chothia C: SUPERFAMILY: HMMs representing all proteins of known structure. SCOP sequence searches, alignments and genome assignments. Nucleic Acids Res. 2002, 30 (1): 268-272. 10.1093/nar/30.1.268.
Cuff JA, Clamp ME, Siddiqui AS, Finlay M, Barton GJ: JPred: a consensus secondary structure prediction server. Bioinformatics. 1998, 14 (10): 892-893. 10.1093/bioinformatics/14.10.892.
Cuff JA, Barton GJ: Application of multiple sequence alignment profiles to improve protein secondary structure prediction. Proteins. 2000, 40 (3): 502-511. 10.1002/1097-0134(20000815)40:3<502::AID-PROT170>3.0.CO;2-Q.
Andrade MA, Perez-Iratxeta C, Ponting CP: Protein repeats: structures, functions, and evolution. J Struct Biol. 2001, 134 (2-3): 117-131. 10.1006/jsbi.2001.4392.
Kobe B, Kajava AV: The leucine-rich repeat as a protein recognition motif. Curr Opin Struct Biol. 2001, 11 (6): 725-732. 10.1016/S0959-440X(01)00266-4.
Gough J: Convergent evolution of domain architectures (is rare). Bioinformatics. 2005, 21 (8): 1464-1471. 10.1093/bioinformatics/bti204.
Wendt KU, Poralla K, Schulz GE: Structure and function of a squalene cyclase. Science. 1997, 277 (5333): 1811-1815. 10.1126/science.277.5333.1811.
Alory C, Balch WE: Molecular basis for Rab prenylation. J Cell Biol. 2000, 150 (1): 89-103. 10.1083/jcb.150.1.89.
Uversky VN: Natively unfolded proteins: a point where biology waits for physics. Protein Sci. 2002, 11 (4): 739-756. 10.1110/ps.4210102.
Dyson HJ, Wright PE: Intrinsically unstructured proteins and their functions. Nat Rev Mol Cell Biol. 2005, 6 (3): 197-208. 10.1038/nrm1589.
Xie H, Vucetic S, Iakoucheva LM, Oldfield CJ, Dunker AK, Obradovic Z, Uversky VN: Functional anthology of intrinsic disorder. 1. Biological processes and functions of proteins with long disordered regions. Journal of Proteome Reseach. 2007,
Xie H, Vucetic S, Iakoucheva LM, Oldfield CJ, Dunker AK, Obradovic Z, Uversky VN: Functional anthology of intrinsic disorder. 3. Ligands, Post-Translational modifications,and diseases associated with intrinsically disordered proteins. Journal of Proteome Research. 2007,
Receveur-Brechot V, Bourhis JM, Uversky VN, Canard B, Longhi S: Assessing protein disorder and induced folding. Proteins. 2006, 62 (1): 24-45. 10.1002/prot.20750.
Linding R, Russell RB, Neduva V, Gibson TJ: GlobPlot: Exploring protein sequences for globularity and disorder. Nucleic Acids Res. 2003, 31 (13): 3701-3708. 10.1093/nar/gkg519.
Boehm M, Bonifacino JS: Adaptins: the final recount. Mol Biol Cell. 2001, 12 (10): 2907-2920.
Ermolaeva MD, Wu M, Eisen JA, Salzberg SL: The age of the Arabidopsis thaliana genome duplication. Plant Mol Biol. 2003, 51 (6): 859-866. 10.1023/A:1023001130337.
Sugino RP, Innan H: Selection for more of the same product as a force to enhance concerted evolution of duplicated genes. Trends Genet. 2006, 22 (12): 642-644. 10.1016/j.tig.2006.09.014.
Pereira-Leal JB, Seabra MC: Evolution of the Rab family of small GTP-binding proteins. J Mol Biol. 2001, 313 (4): 889-901. 10.1006/jmbi.2001.5072.
Rutherford S, Moore I: The Arabidopsis Rab GTPase family: another enigma variation. Curr Opin Plant Biol. 2002, 5 (6): 518-528. 10.1016/S1369-5266(02)00307-2.
Larijani B, Hume AN, Tarafder AK, Seabra MC: Multiple factors contribute to inefficient prenylation of Rab27a in Rab prenylation diseases. J Biol Chem. 2003, 278 (47): 46798-46804. 10.1074/jbc.M307799200.
Jiang SY, Ramachandran S: Comparative and evolutionary analysis of genes encoding small GTPases and their activating proteins in eukaryotic genomes. Physiol Genomics. 2006, 24 (3): 235-251.
Veitia RA: Nonlinear effects in macromolecular assembly and dosage sensitivity. J Theor Biol. 2003, 220 (1): 19-25. 10.1006/jtbi.2003.3105.
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.
Smith TF, Waterman MS: Identification of common molecular subsequences. J Mol Biol. 1981, 147 (1): 195-197. 10.1016/0022-2836(81)90087-5.
Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.
Henikoff S, Henikoff JG: Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci U S A. 1992, 89 (22): 10915-10919. 10.1073/pnas.89.22.10915.
Clamp M, Cuff J, Searle SM, Barton GJ: The Jalview Java alignment editor. Bioinformatics. 2004, 20 (3): 426-427. 10.1093/bioinformatics/btg430.
Madera M, Vogel C, Kummerfeld SK, Chothia C, Gough J: The SUPERFAMILY database in 2004: additions and improvements. Nucleic Acids Res. 2004, 32 (Database issue): D235-9. 10.1093/nar/gkh117.
Finn RD, Marshall M, Bateman A: iPfam: visualization of protein-protein interactions in PDB at domain and amino acid resolutions. Bioinformatics. 2005, 21 (3): 410-412. 10.1093/bioinformatics/bti011.
Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.
Computational Genomics Laboratory . [http://eao.igc.gulbenkian.pt/CGL/]
Schalk I, Zeng K, Wu SK, Stura EA, Matteson J, Huang M, Tandon A, Wilson IA, Balch WE: Structure and mutational analysis of Rab GDP-dissociation inhibitor. Nature. 1996, 381 (6577): 42-48. 10.1038/381042a0.
Wheeler DL, Chappey C, Lash AE, Leipe DD, Madden TL, Schuler GD, Tatusova TA, Rapp BA: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2000, 28 (1): 10-14. 10.1093/nar/28.1.10.
We thank members of the Computational Genomics and Cell Cycle Regulation laboratories at the Instituto Gulbenkian de Ciencia for helpful discussions.
JBPL conceived the study and wrote the manuscript. RR and JBPL performed the analysis.