Ras Superfamily Of GTPases Research Articles

Mutational analysis of the ribosome assembly GTPase RbgA provides insight into ribosome interaction and ribosome-stimulated GTPase activation

Ribosome biogenesis GTPase A protein (RbgA) is an essential GTPase required for the biogenesis of the 50S subunit in Bacillus subtilis. Homologs of RbgA are widely distributed in bacteria and eukaryotes and are implicated in ribosome assembly in the mitochondria, chloroplast and cytoplasm. Cells depleted of RbgA accumulate an immature large subunit that is missing key ribosomal proteins. RbgA, unlike many members of the Ras superfamily of GTPases, lacks a defined catalytic residue for carrying out guanosine triphosphate (GTP) hydrolysis. To probe RbgA function in ribosome assembly, we used a combined bioinformatics, genetic and biochemical approach. We identified a RNA-binding domain within the C-terminus of RbgA that is structurally similar to AmiR–NasR Transcription Anti-termination Regulator (ANTAR) domains, which are known to bind structured RNA. Mutation of key residues in the ANTAR domain altered RbgA association with the ribosome. We identified a putative catalytic residue within a highly conserved region of RbgA, His9, which is contained within a similar PGH motif found in elongation factor Tu (EF-Tu) that is required for GTP hydrolysis on interaction with the ribosome. Finally, our results support a model in which the GTPase activity of RbgA directly participates in the maturation of the large subunit rather than solely promoting dissociation of RbgA from the 50S subunit.

Ubiquitination and activation of a Rab GTPase promoted by a β2-Adrenergic Receptor/HACE1 complex

We and others have shown that trafficking of G-protein-coupled receptors is regulated by Rab GTPases. Cargo-mediated regulation of vesicular transport has received great attention lately. Rab GTPases, which form the largest branch of the Ras GTPase superfamily, regulate almost every step of vesicle-mediated trafficking. Rab GTPases are well-recognized targets of human diseases but their regulation and the mechanisms connecting them to cargo proteins are still poorly understood. Here, we show by overexpression and depletion studies that HACE1, a HECT-domain-containing ubiquitin ligase, promotes the recycling of the β₂-adrenergic receptor (β₂AR), a prototypical G-protein-coupled receptor, through a Rab11a-dependent mechanism. Interestingly, the β₂AR in conjunction with HACE1 triggered ubiquitylation of Rab11a, as observed by western blot analysis. LC-MS/MS experiments determined that Rab11a is ubiquitylated on Lys145. A Rab11a-K145R mutant failed to undergo β₂AR-HACE1-induced ubiquitylation and inhibited the HACE1-mediated recycling of the β₂AR. Rab11a, but not Rab11a-K145R, was activated by β₂AR-HACE1, indicating that ubiquitylation of Lys145 is involved in activation of Rab11a. Co-expression of β₂AR-HACE1 also potentiated ubiquitylation of Rab6a and Rab8a, but not of other Rab GTPases that were tested. We report a novel regulatory mechanism of Rab GTPases through their ubiquitylation, with associated functional effects demonstrated on Rab11a. This suggests a new pathway whereby a cargo protein, such as a G-protein-coupled receptor, can regulate its own trafficking by inducing the ubiquitylation and activation of a Rab GTPase.

RABL-regulated pathways: a new tale in sperm function

RABL-regulated pathways: a new tale in sperm function

Analysis of rab GTPases.

This unit describes methods to identify proteins interacting with, and regulating, Rab GTPases. Rabs form the largest subgroup of Ras superfamily GTPases, and act as molecular switches controlling the specificity of membrane trafficking. The regulation and the signal readout of Rabs are mediated by four groups of proteins, the GDP-GTP exchange factors (GEFs), GTPase activating proteins (GAPs), Rab chaperones, and effector proteins. Rabs are activated at the membrane surface by specific GEFs that promote the exchange of bound GDP to GTP. Effector proteins then bind to the activated Rab GTPase and mediate the cellular response at this membrane surface. Finally, a Rab-specific GAP-protein inactivates the GTPase by catalyzing the hydrolysis of bound GTP to GDP, thus terminating the cellular response. The methods described here are valuable for characterizing these different types of activity, and assigning which Rab they act on.

RAB-Like 2 Has an Essential Role in Male Fertility, Sperm Intra-Flagellar Transport, and Tail Assembly

A significant percentage of young men are infertile and, for the majority, the underlying cause remains unknown. Male infertility is, however, frequently associated with defective sperm motility, wherein the sperm tail is a modified flagella/cilia. Conversely, a greater understanding of essential mechanisms involved in tail formation may offer contraceptive opportunities, or more broadly, therapeutic strategies for global cilia defects. Here we have identified Rab-like 2 (RABL2) as an essential requirement for sperm tail assembly and function. RABL2 is a member of a poorly characterized clade of the RAS GTPase superfamily. RABL2 is highly enriched within developing male germ cells, where it localizes to the mid-piece of the sperm tail. Lesser amounts of Rabl2 mRNA were observed in other tissues containing motile cilia. Using a co-immunoprecipitation approach and RABL2 affinity columns followed by immunochemistry, we demonstrated that within developing haploid germ cells RABL2 interacts with intra-flagella transport (IFT) proteins and delivers a specific set of effector (cargo) proteins, including key members of the glycolytic pathway, to the sperm tail. RABL2 binding to effector proteins is regulated by GTP. Perturbed RABL2 function, as exemplified by the Mot mouse line that contains a mutation in a critical protein–protein interaction domain, results in male sterility characterized by reduced sperm output, and sperm with aberrant motility and short tails. Our data demonstrate a novel function for the RABL protein family, an essential role for RABL2 in male fertility and a previously uncharacterised mechanism for protein delivery to the flagellum.

Role of Electrostatics in Differential Binding of RalGDS to Rap Mutations E30D and K31E Investigated by Vibrational Spectroscopy of Thiocyanate Probes

The human protein Rap1A (Rap) is a member of the Ras superfamily of GTPases that binds to the downstream effector Ral guanine nucleotide dissociation stimulator (RalGDS). Although Ras and Rap have nearly identical amino acid sequences and structures along the effector binding surface, the charge reversal mutation Rap K31E has previously been shown to increase the dissociation constant of Rap-RalGDS docking to values similar to that of Ras-RalGDS docking. This indicates that the difference in charge at position 31 could provide a mechanism for Ral to distinguish two structurally similar but functionally distinct GTPases, which would be of vital importance for appropriate biological function. In this report, vibrational Stark effect spectroscopy, dissociation constant measurements, and molecular dynamics simulations were used to investigate the role that electrostatic field differences caused by the charge reversal mutation Rap K31E play in determining the binding specificity of RalGDS to Rap versus Ras. To do this, six variants of RalGDS carrying a thiocyanate electrostatic probe were docked with three Rap mutants, E30D, K31E, and E30D/K31E. The change in absorption energy of the thiocyanate probe caused by RalGDS docking to these Rap variants was then compared to that observed with wild-type Ras. Three trends emerged: the expected reversion behavior, an additive behavior of the two single mutations, and cancelation of the effects of each single mutation in the double mutant. These observations are explained with a physical model of the position of the thiocyanate probe with respect to the mutated residue based on molecular dynamics simulations.

Phosphorylation of the effector complex HOPS by the vacuolar kinase Yck3p confers Rab nucleotide specificity for vacuole docking and fusion

The homotypic fusion of yeast vacuoles requires the Rab-family GTPase Ypt7p and its effector complex, homotypic fusion and vacuole protein sorting complex (HOPS). Although the vacuolar kinase Yck3p is required for the sensitivity of vacuole fusion to proteins that regulate the Rab GTPase cycle-Gdi1p (GDP-dissociation inhibitor [GDI]) or Gyp1p/Gyp7p (GTPase-activating protein)-this kinase phosphorylates HOPS rather than Ypt7p. We addressed this puzzle in reconstituted proteoliposome fusion reactions with all-purified components. In the presence of HOPS and Sec17p/Sec18p, there is comparable fusion of 4-SNARE (soluble N-ethylmaleimide-sensitive factor attachment protein receptor) proteoliposomes when they have Ypt7p bearing either GDP or GTP, a striking exception to the rule that only GTP-bound forms of Ras-superfamily GTPases have active conformations. However, the phosphorylation of HOPS by recombinant Yck3p confers a strict requirement for GTP-bound Ypt7p for binding phosphorylated HOPS, for optimal membrane tethering, and for proteoliposome fusion. Added GTPase-activating protein promotes GTP hydrolysis by Ypt7p, and added GDI captures Ypt7p in its GDP-bound state during nucleotide cycling. In either case, the net conversion of Ypt7:GTP to Ypt7:GDP has no effect on HOPS binding or activity but blocks fusion mediated by phosphorylated HOPS. Thus guanine nucleotide specificity of the vacuolar fusion Rab Ypt7p is conferred through downstream posttranslational modification of its effector complex.

RasGRF1 regulates proliferation and metastatic behavior of human alveolar rhabdomyosarcomas

The involvement of the Ras superfamily of GTPases in the pathogenesis of rhabdomysarcoma (RMS) is not well understood. While mutant H-Ras leads to embryonal RMS (ERMS) formation in experimental animals and in Costello syndrome patients, no data exists on the potential role of Ras GTPases in the pathogenesis of alveolar RMS (ARMS). To address this issue better, we focused on the role of the GTP exchange factor RasGRF1 in this process. We observed that, in comparison to normal skeletal muscle cells, RasGRF1 mRNA is upregulated in the majority of human ARMS cell lines and subsequently confirmed its high expression in patient samples. By employing confocal microscopy analysis, we observed RasGRF1 accumulation in cell filopodia, which suggests its involvement in ARMS cell migration. Furthermore, we observed that RasGRF1 becomes phosphorylated in ARMS after stimulation by several pro-metastatic factors, such as SDF-1 and HGF/SF, as well as after exposure to growth-promoting Igf-2 and insulin. More importantly, activation of RasGRF1 expression correlated with activation of p42/44 MAPK and AKT. When the expression of RasGRF1 was down-regulated in ARMS cells by an shRNA strategy, these RasGRF1-kd RMS cells did not respond to stimulation by SDF-1, HGF/SF, Igf-2 or insulin by phosphorylation of p42/44 MAPK and AKT and lost their chemotactic responsiveness; however, their adhesion was not affected. We also observed that RasGRF1-kd ARMS cells proliferated at a very low rate in vitro, and, more importantly, after inoculation into immunodeficient SCID/beige inbred mice they formed significantly smaller tumors. We conclude that RasGRF1 plays an important role in ARMS pathogenesis and is a new potential therapeutic target to inhibit ARMS growth.

Abstract 2147: The RRas2/TC21 GTPase is essential for breast tumorigenesis and lung metastasis

Abstract Rras2/TC21 is a Ras superfamily GTPase with still unknown functions in vivo. Several observations suggest that this protein may be involved in tumorigenic processes. Thus, it has been shown that the expression of constitutively active RRas2 mutants leads to high levels of cell transformation in vitro. Moreover, this GTPase has been found overexpressed in a number of tumors, including lymphomas, oral and esophageal carcinoma, and breast cancer. To address the role of this protein in epithelial tumors, we have carried out an extensive characterization of this GTPase using RRas2 knockout mice and shRNA knockdown tumor cells. These studies revealed that RRas2 is not important for the development and/or malignant progression of skin and oral squamous cell carcinoma. Instead, we have observed that the RRas2 deficiency severily impairs breast cancer tumorigenesis and the tropism of metastatic cancer cells to the lung. The characterization of the role of RRas2 in primary tumorigenesis indicates that this GTPase is important for the proliferation and angiogenesis of tumors. In the case of metastasis, RRas2 is important for the survival and proliferation of cancer cells in the lung parenchyma. By contrast, this GTPase is not required for the lymph node infiltration, intravasation, or extravasation of cancer cells. Interestingly, the requirement of RRas2 is still maintained in breast cancer cells that express active mutants of K-Ras, indicating that RRas2 plays essential roles downstream of both Her2 and K-Ras. Signaling experiments indicate that RRas2 contributes to those processes by modulating the activation of phosphatidylinositol 3-kinase ≤ and downstream targets such as Akt, FoxO1/FoxO3 and S6 ribosomal protein. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 2147. doi:1538-7445.AM2012-2147

GTPase Networks in Membrane Traffic

Members of the Rab or ARF/Sar branches of the Ras GTPase superfamily regulate almost every step of intracellular membrane traffic. A rapidly growing body of evidence indicates that these GTPases do not act as lone agents but are networked to one another through a variety of mechanisms to coordinate the individual events of one stage of transport and to link together the different stages of an entire transport pathway. These mechanisms include guanine nucleotide exchange factor (GEF) cascades, GTPase-activating protein (GAP) cascades, effectors that bind to multiple GTPases, and positive-feedback loops generated by exchange factor-effector interactions. Together these mechanisms can lead to an ordered series of transitions from one GTPase to the next. As each GTPase recruits a unique set of effectors, these transitions help to define changes in the functionality of the membrane compartments with which they are associated.

The GTPase TcRjl of the human pathogen Trypanosoma cruzi is involved in the cell growth and differentiation

The protozoan parasite Trypanosoma cruzi, the etiological agent of Chagas Disease, undergoes through a complex life cycle where rounds of cell division and differentiation occur initially in the gut of triatominae vectors and, after transmission, inside of infected cells in vertebrate hosts. Members of the Ras superfamily of GTPases are molecular switches which play pivotal regulatory functions in cell growth and differentiation. We have previously described a novel GTPase in T. cruzi, TcRjl, which belongs to the RJL family of Ras-related GTP binding proteins. Here we show that most of TcRjl protein is found bound to GTP nucleotides and may be locked in this stage. In addition, we show that TcRjl is located close to the kinetoplast, in a region corresponding possibly to flagellar pocket of the parasite and the expression of a dominant-negative TcRjl construct (TcRjlS37N) displays a significative growth phenotype in reduced serum medium. Remarkably, overexpression of TcRjl inhibits differentiation of epimastigotes to trypomastigote forms and promotes the accumulation of intermediate differentiation stages. Our data suggest that TcRjl might play a role in the control of the parasite growth and differentiation.

Chemical Biology of Lipidated Proteins

Many signaling proteins such as the members of the Ras superfamily of GTPases are posttranslationally modified by covalent attachment of lipid groups, which is crucial for the correct localization and function of these proteins. Numerous lipidated proteins are oncogens often found mutated in several human cancers. Therefore, several therapeutic strategies have been developed based on the inhibition of the enzymes involved in these lipidation steps. Here, we will summarize the results on protein lipidation inhibition, mainly focusing on the small molecules targeting the isoprenylation and acylation of proteins.

Mechanism of activation and inactivation of Gq/phospholipase C-β signaling nodes.

The phospholipase C (PLC) isozymes catalyze conversion of phosphatidylinositol 4,5-bisphosphate (PtdIns(4,5)P2) into the Ca2+-mobilizing second messenger inositol 1,4,5-trisphosphate (inositol (1,4,5)P3) and the protein kinase-activating second messenger diacylglycerol (DAG).1 Many PLC-dependent cellular responses occur in addition to those mediated by these classical second messengers since the activities of a broad range of membrane, cytosolic, and cytoskeletal proteins are regulated by PtdIns(4,5)P2 binding.2–4 Mammals express 13 different PLCs (Figure 1), which differ markedly in their modes of upstream regulation and physiological functions.5 Many growth factors, antigens, and other extracellular stimuli signal through tyrosine phosphorylation of the PLC-γ isozymes, whereas an even larger group of hormones, neurotransmitters, chemoattractant and chemosensory molecules, and other extracellular stimuli promote physiological effects through heterotrimeric G protein-dependent activation of the PLC-β isozymes. These and other forms of receptor-promoted signaling pathways also communicate to Ras superfamily GTPases, which in turn directly activate certain PLC isoyzmes. Thus, PLC-dependent signal transduction provides one of the major fabrics for communication of cells, and delineation of its function at all levels from the intact animal to the atomic resolution of mechanism is fundamental to understanding mammalian biology. Figure 1 Conserved domain structure of the mammalian PLC isozymes. The 13 functional human PLC isozymes were aligned on the basis of the conservation of the protein sequence, and a dendrogram was constructed to cluster similar sequences into shared branches. The ... Many excellent reviews on PLC-dependent signaling are available that focus, for example, on early aspects of the discovery and function of receptor-promoted formation of Ins(1,4,5)P3/diacylglycerol and mobilization of Ca2+,1,6,7 classification and regulation of the PLC isozymes,8–10 tyrosine phosphorylation-dependent activation,11,12 activation through Ras superfamily GTPases,5,13,14 physiology,15,16 and structure/function.17–19 We have limited this review to a general introduction of the domain and structural features of the PLC isozymes and then consider in detail recent advances made in understanding the mechanisms through which Gα-subunits of the Gq family bind to and activate PLC-β isozymes and how the PLC-β isozymes in turn promote inactivation of this signaling complex by stimulating GTP hydrolysis by the GTP-activated G protein.

Arl13b regulates ciliogenesis and the dynamic localization of Shh signaling proteins

Arl13b, a ciliary protein within the ADP-ribosylation factor family and Ras superfamily of GTPases, is required for ciliary structure but has poorly defined ciliary functions. In this paper, we further characterize the role of Arl13b in cilia by examining mutant cilia in vitro and determining the localization and dynamics of Arl13b within the cilium. Previously, we showed that mice lacking Arl13b have abnormal Sonic hedgehog (Shh) signaling; in this study, we show the dynamics of Shh signaling component localization to the cilium are disrupted in the absence of Arl13b. Significantly, we found Smoothened (Smo) is enriched in Arl13b-null cilia regardless of Shh pathway stimulation, indicating Arl13b regulates the ciliary entry of Smo. Furthermore, our analysis defines a role for Arl13b in regulating the distribution of Smo within the cilium. These results suggest that abnormal Shh signaling in Arl13b mutant embryos may result from defects in protein localization and distribution within the cilium.

A Role for Rac3 GTPase in the Regulation of Autophagy

The process of autophagy is situated at the intersection of multiple cell signaling pathways, including cell metabolism, growth, and death, and hence is subject to multiple forms of regulation. We previously reported that inhibition of isoprenylcysteine carboxylmethyltransferase (Icmt), which catalyzes the final step in the post-translational prenylation of so-called CAAX proteins, results in the induction of autophagy which enhances cell death in some cancer cells. In this study, using siRNA-mediated knockdown of a group of small GTPases that are predicted Icmt substrates, we identify Rac3 GTPase as a negative regulator of the process of autophagy. Knockdown of Rac3, but not the closely related isoforms Rac1 and Rac2, results in induction of autophagy. Ectopic expression of Rac3, significantly rescues cells from autophagy and cell death induced by Icmt inhibition, strengthening the notion of an isoform-specific autophagy regulatory function of Rac3. This role of Rac3 was observed in multiple cell lines with varying Rac subtype expression profiles, suggesting its broad involvement in the process. The identification of this less-studied Rac member as a novel regulator provides new insight into autophagy and opens opportunities in identifying additional regulatory inputs of the process.

The human Rgr oncogene is overexpressed in T-cell malignancies and induces transformation by acting as a GEF for Ras and Ral

The Ras superfamily of GTPases is involved in the modification of many cellular processes including cellular motility, proliferation and differentiation. Our laboratory has previously identified the RalGDS related (Rgr) oncogene in a DMBA-induced rabbit squamous cell carcinoma and its human orthologue, hRgr. In the present study, we analyzed the expression levels of the human hRgr transcript in a panel of human hematopoietic malignancies and found that a truncated form (diseased-truncated; Dtr-hrgr) was significantly overexpressed in many T-cell derived neoplasms. Although the Rgr proto-oncogene belongs to the RalGDS family of guanine nucleotide exchange factors (GEFs), we show that upon the introduction of hRgr into fibroblast cell lines it is able to elicit the activation of both Ral and Ras GTPases. Moreover, in vitro guanine nucleotide exchange assays confirm that hRgr promotes Ral and Ras activation through GDP dissociation, which is a critical characteristic of GEF proteins. hRgr has guanine nucleotide exchange activity for both small GTPases and this activity was reduced when a point mutation within the catalytic domain (CDC25) of the protein, (cd) Dtr-hRgr, was utilized. These observations prompted the analysis of the biological effects of hRgr and (cd) hRgr expression in cultured cells. Here, we show that hRgr increases proliferation in low serum, increases invasion, reduces anchorage dependence, and promotes the progression into S phase of the cell cycle; properties that are abolished or severely reduced in the presence of the catalytic dead mutant. We conclude that the ability of hRgr to activate both Ral and Ras is responsible for its transformation-inducing phenotype and it could be an important contributor in the development of some T cell malignancies.

Statins inhibit T-acute lymphoblastic leukemia cell adhesion and migration through Rap1b

Statins are known to inhibit signaling of Ras superfamily GTPases and reduce T cell adhesion to ICAM-1. Here, we address the hypothesis that statins affect T cell adhesion and migration by modulating the function of specific GTPases. Statins inhibit the synthesis of mevalonic acid, which is required for farnesyl and geranylgeranyl isoprenoid synthesis. Ras superfamily GTPases are post-translationally isoprenylated to facilitate their anchorage to membranes, where they function to stimulate signal transduction processes. We demonstrate that 1 μM statin inhibits the adhesion, migration, and chemotaxis of the T-ALL cell line CCRF-CEM and TEM of CCRF-CEM and PEER T-ALL cells, but higher statin concentrations are needed to inhibit adhesion of primary T cells. Similar effects are observed following treatment with GGTI-298 or RNA interference-mediated knockdown of Rap1b but not Rap1a, Rac1, Rac2, RhoA, or Cdc42. Statins also alter Rap1 activity and Rap1b localization. Rap1 levels are higher in primary T cells than T-ALL cells, which could explain their reduced sensitivity to statins. These results demonstrate for the first time that the closely related Rap1a and Rap1b isoforms have different functions and suggest that statins or Rap1b depletion could be used to reduce tissue invasion in T-ALL.

SUMO wrestling in cell movement

SUMO wrestling in cell movement

Exosomes account for vesicle-mediated transcellular transport of activatable phospholipases and prostaglandins

Exosomes are bioactive vesicles released from multivesicular bodies (MVB) by intact cells and participate in intercellular signaling. We investigated the presence of lipid-related proteins and bioactive lipids in RBL-2H3 exosomes. Besides a phospholipid scramblase and a fatty acid binding protein, the exosomes contained the whole set of phospholipases (A2, C, and D) together with interacting proteins such as aldolase A and Hsp 70. They also contained the phospholipase D (PLD) / phosphatidate phosphatase 1 (PAP1) pathway leading to the formation of diglycerides. RBL-2H3 exosomes also carried members of the three phospholipase A2 classes: the calcium-dependent cPLA(2)-IVA, the calcium-independent iPLA(2)-VIA, and the secreted sPLA(2)-IIA and V. Remarkably, almost all members of the Ras GTPase superfamily were present, and incubation of exosomes with GTPgammaS triggered activation of phospholipase A(2) (PLA(2))and PLD(2). A large panel of free fatty acids, including arachidonic acid (AA) and derivatives such as prostaglandin E(2) (PGE(2)) and 15-deoxy-Delta(12,14)-prostaglandinJ(2) (15-d PGJ(2)), were detected. We observed that the exosomes were internalized by resting and activated RBL cells and that they accumulated in an endosomal compartment. Endosomal concentrations were in the micromolar range for prostaglandins; i.e., concentrations able to trigger prostaglandin-dependent biological responses. Therefore exosomes are carriers of GTP-activatable phospholipases and lipid mediators from cell to cell.

T-RMSD: A Fine-grained, Structure-based Classification Method and its Application to the Functional Characterization of TNF Receptors

This study addresses the relation between structural and functional similarity in proteins. We introduce a novel method named tree based on root mean square deviation (T-RMSD), which uses distance RMSD (dRMSD) variations to build fine-grained structure-based classifications of proteins. The main improvement of the T-RMSD over similar methods, such as Dali, is its capacity to produce the equivalent of a bootstrap value for each cluster node. We validated our approach on two domain families studied extensively for their role in many biological and pathological pathways: the small GTPase RAS superfamily and the cysteine-rich domains (CRDs) associated with the tumor necrosis factor receptors (TNFRs) family. Our analysis showed that T-RMSD is able to automatically recover and refine existing classifications. In the case of the small GTPase ARF subfamily, T-RMSD can distinguish GTP- from GDP-bound states, while in the case of CRDs it can identify two new subgroups associated with well defined functional features (ligand binding and formation of ligand pre-assembly complex). We show how hidden Markov models (HMMs) can be built on these new groups and propose a methodology to use these models simultaneously in order to do fine-grained functional genomic annotation without known 3D structures. T-RMSD, an open source freeware incorporated in the T-Coffee package, is available online.