Biochemical Clustering of Monomeric GTPases of the Ras Superfamily *,S

Marie-Elaine Caruso{ddagger},§, Sarah Jenna{ddagger},||, Simon Beaulne{ddagger}, Eun-Hye Lee,||, Anne Bergeron**, Cédric Chauve**, Philippe Roby{ddagger}{ddagger}, Jean-François Rual§§, David E. Hill§§, Marc Vidal§§, Roger Bossé{ddagger}{ddagger},¶¶ and Eric Chevet{ddagger},§,||||

From the {ddagger} Organelle Signaling Laboratory, Department of Surgery and Montreal Proteomics Network, McGill University, Montreal, Quebec H3A 2B2, Canada, the ** Laboratoire de Combinatoire et d’Informatique Mathématique, University of Quebec at Montreal, Montreal, Quebec H3C 3P8, Canada, {ddagger}{ddagger} PerkinElmer BioSignal, Montreal, Quebec H3J 1R4, Canada, and the §§ Centre for Cancer Systems Biology and Department of Cancer Biology, Dana-Farber Cancer Institute and the Department of Genetics, Harvard Medical School, Boston, Massachusetts 02115


    ABSTRACT
 TOP
 ABSTRACT
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
To date phylogeny has been used to compare entire families of proteins based on their nucleotide or amino acid sequence. Here we developed a novel analytical platform allowing a systematic comparison of protein families based on their biochemical properties. This approach was validated on the Rho subfamily of GTPases. We used two high throughput methods, referred to as AlphaScreenTM and FlashPlate®, to measure nucleotide binding capacity, exchange, and hydrolysis activities of small monomeric GTPases. These two technologies have the characteristics to be very sensitive and to allow homogenous and high throughput assays. To analyze and integrate the data obtained, we developed an algorithm that allows the classification of GTPases according to their enzymatic activities. Integration and hierarchical clustering of these results revealed unexpected features of the small Rho GTPases when compared with primary sequence-based trees. Hence we propose a novel phylobiochemical classification of the Ras superfamily of GTPases.


Ras GTPases are small GTP-binding proteins involved in diverse cellular processes such as apoptosis, cell proliferation/differentiation, cytoskeleton reorganization, and membrane trafficking. These proteins cycle between an inactive GDP-bound form and an active GTP-bound form. Under their GTP-bound form their active conformation allows them to interact with effector molecules and to generate specific biological responses (1). Two intrinsic enzymatic activities regulate the balance between GTP-bound and GDP-bound conformations: 1) GDP/GTP exchange and 2) GTP hydrolysis activities (2). Whereas the intrinsic enzymatic activities have been characterized for a few of these proteins, such studies were carried out in non-systematic manners and by different research groups (36).

The analysis of full genome sequences has allowed large scale biology experimental approaches consisting of systematic characterization of entire families of proteins instead of investigating them individually. Developing such "functional proteomic" approaches requires the design of high throughput and versatile experimental procedures (7, 8). For the study of small G proteins, conventional filtration assays, although powerful at small scale, are not appropriate for systematic studies. Indeed this methodology has a low to very low throughput, is not versatile, and requires large amounts of reagents. To overcome these issues and characterize GTP-binding proteins of the Ras superfamily in a systematic manner, we developed a robust assay platform using two well adopted high throughput technologies, AlphaScreenTM (9) and FlashPlate® (10).

AlphaScreen is a bead-based non-radioactive and homogenous proximity assay used to measure interaction between biological binding partners. The principle of this technology relies on the use of a Donor bead and an Acceptor bead that generate a light signal when brought into proximity (<200 nm). Upon laser excitation at 680 nm, the Donor beads, containing a photosensitizer, will generate short lived singlet oxygen that can diffuse only a short distance before returning to the ground state. The Acceptor beads, containing chemiluminescers and fluorophores, will react with this singlet oxygen and will emit an amplified light signal measurable at ~600 nm. AlphaScreen (PerkinElmer) provides highly versatile, sensitive, and homogeneous assays that allow us to perform studies at a higher throughput and at a lower cost.

FlashPlates are white polystyrene microplates in which the interior of wells is coated with a thin layer of polystyrene-based scintillation reagent. The principle of this homogeneous radiometric technology is based on proximity between a radioligand and the scintillation reagent. A target protein, such as antibodies or GSH, is also labeled on the wall of the plate and will bring the protein of interest close to the scintillation reagent. The interaction between the radioisotopes and the proteins of interest will activate the scintillation reagent and emit a luminescent signal. FlashPlate (PerkinElmer) technology reduces considerably the amount of reagent to use and eliminates the time-consuming washes normally needed to separate bound from unbound radioligands.

In this study we took advantage of these two technologies to characterize in a systematic manner four biochemical properties of the Ras family of small G proteins and we reclassified them according to their activities instead of their primary amino acid sequences. Due to their high level of interspecies conservation, we selected the Rho GTPases as prototype proteins to develop and validate our assays.


    EXPERIMENTAL PROCEDURES
 TOP
 ABSTRACT
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Materials
[{gamma}-35S]GTP{gamma}S and [{gamma}-33P]GTP were purchased form PerkinElmer Life Sciences. ATP, GDP, GTP, and GTP{gamma}S1 were purchased form Sigma. Goat anti-GST antibodies were from Amersham Biosciences.

Plasmids and Constructs
Caenorhabditis elegans CDC-42 (RO7G3.1), RAC-2 (KO3D3.10b), CED-10 (C09G12.8b), MIG-2 (C35C5.4), and RHO-1 (Y51H4A.3) plus one putative uncharacterized GTPase, CRP-1 (Y32F6B.3), were obtained from the C. elegans ORFeome (11). The mRAC-1 and mRAC-1N17 were described previously (12). The above ORFs were recombined in a bacterial expression vector (pGEX-2TK) to produce N-terminally tagged GST fusion proteins using the Gateway technology (Invitrogen).

Nucleotide Biotinylation
GTP{gamma}S and GDP were biotinylated, respectively, using biotin-PEG-maleinate (Pierce) dissolved in MES, pH 6, and biotin-PEG-COO-NHS (Pierce) dissolved in carbonate buffer, pH 8.5, according to the manufacturer’s instructions. For both nucleotides, reaction mixtures were incubated for 2 h at 37 °C and purified using HPLC. Concentrations were measured by optical density at 260 nm. Biotinylated nucleotides were validated for their biological properties toward GTPases (data not shown).

Expression and Purification of Recombinant Proteins
Bacterial recombinant GST-GTPases were expressed and purified as described previously (13).

AlphaScreen Assays
AlphaScreen assays were performed in Costar 384-well microplates in a final reaction volume of 25 µl. Streptavidin Donor beads and protein G Acceptor beads (PerkinElmer BioSignal) were used at a final concentration of 0.02 mg/ml per well. The assays were performed in AlphaScreen buffer (50 mM HEPES, pH 7.4, 100 mM NaCl, 5 mM MgCl2, 1 mM DTT, and 1 mg/ml BSA). All incubations were performed at 23 °C. Laser excitations were carried out at 680 nm, and readings were performed at 520–620 nm using an AlphaQuest reader (Packard).

Competition Assays—
During purification, GST-GTPases were unloaded from their intrinsic nucleotide by adding 5 mM EDTA in lysate before purification. Different concentrations of nucleotides (ranging from 0.4 nM to 3 µM) were incubated with a 30 nM concentration of unloaded GST-GTPases and a 30 nM concentration of a biotinylated derivative of GTP{gamma}S (b-GTP{gamma}S). 25 nM goat anti-GST antibody (Amersham Biosciences) was added and incubated for 20 min. Donor and Acceptor beads were added simultaneously and incubated for 1 h before reading.

Exchange Activity—
Purified GST-GTPases were preloaded with a biotinylated derivative of GDP (b-GDP) as follows. 150 nM proteins were incubated with 100 nM b-GDP in loading buffer (25 mM Tris, pH 7.4, 100 mM NaCl, 0.1 mM DTT, and 5 mM EDTA) for 30 min at 30 °C. Reactions were then placed on ice, and 10 mM MgCl2 was added. 5 µl of preloaded GST-GTPases (final concentration, 30 nM) were then added in microplate wells with different concentration of unlabeled GTP (ranging from 0.4 nM to 3 µM). The plate was incubated for 20 min, and 25 nM goat anti-GST (Amersham Biosciences) was added. After 20 min of incubation, Donor and Acceptor beads were added simultaneously, and the plate was analyzed after 1 h of incubation at 23 °C.

FlashPlate Assays
These assays were performed in glutathione-coated 96-wells FlashPlates (PerkinElmer). Bacterial lysates containing GST-GTPases were diluted with PBS (50 µl) and incubated in the FlashPlate for 45 min at 4 °C in the presence of 5 mM EDTA to allow the binding of GST-GTPases on the GSH plate and to unload the protein from its intrinsic nucleotide.

Binding Assay—
The lysate was removed, and the wells were washed two times with loading buffer (20 mM Tris, pH 7.4, 100 mM NaCl, and 0.1 mM DTT). Then 50 µl of loading buffer containing different concentrations of [{gamma}-35S]GTP{gamma}S (20–1000 nCi) were added. The signal produced was monitored on a ß scintillation counter (Packard) every 5 min until the plateau was reached. The plate was placed back on ice, and 10 mM MgCl2 was added. The radioactive reaction mixture was then removed from the plate, and 50 µl of stopping buffer (20 mM Tris, pH 7.4, 0.1 mM DTT, 5 mM MgCl2, 1 mg/ml BSA, and 100 µM GTP{gamma}S) were added. Signal was monitored on a ß scintillation counter (Packard) every 5 min for another 30 min to measure nucleotide dissociation.

Hydrolysis Activity—
Bacterial lysate containing EDTA was replaced by 50 µl of loading buffer (20 mM Tris, pH 7.4, 100 mM NaCl, and 0.1 mM DTT) containing either 500 nCi of [{gamma}-35S]GTP{gamma}S or 500 nCi of [{gamma}-33P]GTP. The interaction between GTPases and these radionucleotides was carried out for 30 min after which 10 mM MgCl2 was added. The radioactive reaction mixture was then replaced by 50 µl of hydrolysis buffer (20 mM Tris, pH 7.4, 100 mM NaCl, 5 mM MgCl2, 1 mM DTT, 0.01% BSA, and 1 µM GTP), and the radioactivity was monitored for another 60 min using a ß scintillation counter (Packard).

Data Clustering, Tree Representation, and Comparison
C. elegans GTPase amino acid sequences were aligned using ClustalW (Blosum30 (Henikoff) matrix).2 Matrices were computed from the alignment using PHYLIP ProtDist according to the Jones-Taylor-Thornton model. The hierarchical clustering tree was then computed using the PHYLIP Neighbor command based on the UPGMA method and visualized using TreeView. Biochemical values obtained for each GTPase were normalized (between 0 and 1) and used to generate distance matrices based on activity(ies) similarities (alone or in combination). These matrices were calculated based on Euclidian distances, and hierarchical clustering trees were then computed using the PHYLIP Neighbor command based on the UPGMA method and visualized using TreeView. Tree comparison was carried out using PHYLIP Treedist.


    RESULTS
 TOP
 ABSTRACT
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Rho GTPases Identified in the C. elegans Genome—
The ORFeome projects, which consist of the individual cloning of nearly all protein-encoding ORFs from the ATG to the stop codon predicted from genome sequences, are the cornerstone of system biology approaches. The C. elegans ORFeome project is one of the most advanced and had reached more than 84% completion in 2003 (11). In C. elegans, six genes have been predicted to encode potential Rho GTPases, and the corresponding amino acid sequences were aligned and compared with the mouse RAC-1 protein sequence as shown in Fig. 1A. The corresponding ORFs were retrieved from the C. elegans ORFeome (14, 15) and expressed as N-terminally tagged GST fusion proteins in Escherichia coli (DH5{alpha}). The GST fusion proteins were expressed at high levels and purified to homogeneity (Fig. 1B). Using these small G proteins, we developed four distinct assays to measure their affinity for guanine nucleotides as well as their exchange and hydrolysis activities.



View larger version (68K):
[in this window]
[in a new window]
 
FIG. 1. Analysis of the six C. elegans Rho GTPases. A, amino acid sequence alignment (ClustalW) of the six C. elegans RHO GTPases. B, SDS-PAGE analysis of the six GST-tagged RHO GTPases identified in the C. elegans genome. GST-GTPases were expressed in bacteria and purified to homogeneity using glutathione-Sepharose 4B resin.

 
Small G Protein Affinity for Nucleotides—
The first two assays, developed using both FlashPlates and AlphaScreen, aimed at measuring the affinity of small G proteins for GTP and other nucleotides. GSH-coated FlashPlates were used to perform saturation binding assays where increasing concentrations of non-hydrolyzable [{gamma}-35S]GTP{gamma}S were incubated in the presence of a constant amount of nucleotide-free GST-GTPase. GST alone was also used as a negative control (see "Experimental Procedures"; Fig. 2A). When [{gamma}-35S]GTP{gamma}S interacts with the proteins associated with the plate, ß particles emitted from 35S activate the scintillation reagent embedded in the walls of the microplate. A luminescent signal is then emitted and measured. As shown in Fig. 2B, the signal intensity increased proportionally to the radionucleotide concentration and reached a plateau. An apparent Kd (Kd(app)) corresponding to the concentration of [{gamma}-35S]GTP{gamma}S necessary to reach 50% of the maximum signal was measured (Fig. 2B) in a prototypical experiment using CDC-42 (Kd(app) = 7.2 nM). Kd(app) values obtained with the six C. elegans putative Rho GTPases (CDC-42, CED-10, RAC-2, RHO-1, CRP-1, and MIG-2) were then compared together to set a rank order of GTPase GTP binding potency. In addition, the well characterized mRAC-1 (12, 16) and the GTP binding-deficient mutant mRAC-1N17 (17) were used as controls. Our results show that these proteins clustered in two groups according to their affinity toward GTP{gamma}S: on one hand CDC-42, RAC-2, CED-10, RHO-1, and mRAC-1 showed a high [{gamma}-35S]GTP{gamma}S binding with Kd(app) values ranging from 2.8 to 7.4 nM, and on the other hand MIG-2, CRP-1, and mRAC-1N17 showed a low affinity for [{gamma}-35S]GTP{gamma}S with Kd(app) ranging from 11.7 to 14.5 nM (Fig. 2C). The results obtained with mRAC-1N17 were consistent with previous reports showing that this mutant displays low GTP affinity (17).



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 2. [{gamma}-35S]GTP{gamma}S saturation binding assay. A, FlashPlate experimental design for [{gamma}-35S]GTP{gamma}S binding assay. Nucleotide-free GST-GTPases were fixed on a GSH FlashPlate and incubated with different concentrations of [{gamma}-35S]GTP{gamma}S. Radioligand interacting with the GST-GTPase will activate the FlashPlate and generate a light signal as described under "Experimental Procedures." B, the saturation curve obtained with the binding assay shows the kinetic of association between the [{gamma}-35S]GTP{gamma}S and GST-CDC-42 (•) or GST alone ({square}). The Kd(app) value corresponds to the concentration at which 50% of the maximum binding is reached. This graph represents a single experiment done in triplicate. C, comparison of the Kd(app) values obtained in the [{gamma}-35S]GTP{gamma}S loading assay for six C. elegans GTPases (CDC-42, RAC-1, CED-10, RHO-1, MIG-2, and CRP-1) and two mouse GTPases (mRAC-1 and the GTP binding-deficient mRAC-1N17). The Kd(app) values represent the average of at least three independent experiments performed in triplicate.

 
Our next objective was to develop competition assays to determine specific nucleotide affinities. These assays were carried out using both AlphaScreen and FlashPlate. Both platforms gave comparable results for CDC-42 affinities for GTP or GDP (Supplemental Fig. 1). To compare the six potential C. elegans Rho GTPases, we used AlphaScreen because of its higher throughput potential and its non-radioactive applications. In this assay, b-GTP{gamma}S was produced and used as a tracer binding to GTP-unloaded GST-GTPases as described under "Experimental Procedures." In our system, the binding of b-GTP{gamma}S to GST-GTPases brings anti-GST Acceptor beads into proximity of streptavidin-coated Donor beads thus allowing the generation of an AlphaScreen signal (Fig. 3A). Competition isotherms are generated by adding increasing concentrations of unlabeled nucleotides, which progressively prevent b-GTP{gamma}S binding to GST-GTPase. A concentration-dependent signal decrease is then observed (Fig. 3B). Competition curves with GTP were used to generate IC50 values, which are apparent affinity constants: an IC50 value corresponds to the GTP concentration inducing 50% of signal loss. As shown in Fig. 3B, GTP competed the binding of b-GTP{gamma}S to CDC-42 with an IC50 value of 45 nM. IC50 values specific for each GTPases were then measured except for mRAC-1N17 for which an outranged signal was detected due to its expected low affinity for b-GTP{gamma}S. When the IC50 values were compared together, Rho GTPases clustered in the same two groups previously reported with saturation binding studies using the FlashPlate. High affinity Rho GTPases showed IC50 values ranging from 11.8 to 33.5 nM, whereas the lowest affinity Rho GTPases generated IC50 values varying from 97.6 to 294.5 nM (Fig. 3C). ATP was used as a nonspecific nucleotide control. As expected, ATP was not able to compete the binding of b-GTP{gamma}S to the GTPases. In all the assays, GST was included as another negative control to measure background level.



View larger version (18K):
[in this window]
[in a new window]
 
FIG. 3. b-GTP{gamma}S competition assay. A, AlphaScreen experimental design for the b-GTP{gamma}S competition assay. Nucleotide-free GST-GTPases were incubated with a constant concentration of b-GTP{gamma}S and increasing concentrations of unlabeled competitors. The streptavidin Donor bead (D) recognizes b-GTP{gamma}S, and the protein G Acceptor bead (A) binds to the GST-GTPase via an anti-GST antibody. When the b-GTP{gamma}S interacts with the protein, the two beads are in close proximity, and a light signal is detectable upon laser excitation of the Donor bead. B, b-GTP{gamma}S competition curves obtained with either ATP ({blacktriangleup}) or GTP (•) on CDC-42 or GST ({square}). The curves are the result of a single experiment done in triplicate. The IC50 value corresponds to the competitor concentration required for the loss of 50% of the maximum signal (obtained in the absence of competitor). C, comparison of the IC50 values obtained for six C. elegans GTPases (CDC-42, RAC-1, CED-10, RHO-1, MIG-2, and CRP-1) and two mouse GTPases (mRAC-1 and the GTP binding-deficient mRAC-1N17) using the b-GTP{gamma}S competition assay. The IC50 values for the competition assay represent the average of at least three independent experiments performed in triplicate. NA, not applicable; CPS, counts per second.

 
Small G Protein Exchange Activity—
A third assay aiming to measure the GDP/GTP exchange activity was developed using AlphaScreen. GST-GTPases were loaded with b-GDP and incubated with increasing concentrations of GTP. Donor and Acceptor beads were added, and the exchange of b-GDP to GTP was measured as described under "Experimental Procedures" (Fig. 4A). As described above for the b-GTP{gamma}S binding assay, exchange of b-GDP to GTP led to a signal decrease. The exchange activity was also characterized by an IC50 value obtained when increasing concentrations of GTP were used to compete the binding of b-GDP to a Rho GTPase such as CDC-42 (Fig. 4B). The IC50 value corresponds to the concentration of GTP necessary for dissociating 50% of the b-GDP initially bound. To compare exchange activities for various GTPases, the specific activity of the GTPases was measured, and a similar amount of GTPase·b-GDP complexes was used in the assay. Parallel analyses of the exchange activities displayed by the Rho GTPases showed that CDC-42, RAC-2, CED-10, RHO-1, and mRAC-1 have high exchange activities with IC50 values ranging from 0.65 to 13.3 nM, whereas MIG-2 and mRAC-1N17 showed lower exchange activities (IC50 values of 49 and 2000 nM, respectively) consistent with their lower affinity for GTP (Fig. 4C). No significant exchange activity was detected for CRP-1. This phenomenon is attributable to the lower affinity of this protein to GTP. Additional experiments to measure the relative affinity of the GTPases for nucleotides using competition assay (as in Fig. 3A) revealed that CRP-1 had a lower affinity for the GTP than GDP (Supplemental Fig. 2B). These results were consistent with the data previously obtained using filtration assay and showing that CRP-1 had a very low exchange activity and atypical characteristics (18).



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 4. b-GDP/GTP exchange activity assay. A, AlphaScreen experimental design for the measure of the exchange activity. GST-GTPases were loaded with b-GDP and incubated with different concentrations of GTP. The same AlphaScreen beads as in Fig. 2 were used. When b-GDP is exchanged to GTP, the distance between the Donor and Acceptor beads increases, and a loss of signal is measured. B, curves obtained in a single assay done in triplicate for GST-CDC-42 ({blacksquare}) and GST alone ({triangledown}) using this b-GDP/GTP exchange assay. The IC50 value corresponds to the GTP concentration where 50% of the b-GDP is exchanged to a GTP. C, comparison of the IC50 values obtained for six C. elegans GTPases (CDC-42, RAC-1, CED-10, RHO-1, MIG-2, and CRP-1) and two mouse GTPases (mRAC-1 and the GTP binding-deficient mRAC-1N17) using the b-GDP/GTP exchange assay. The IC50 values for the exchange activity correspond to the average of at least three independent experiments performed in triplicate. CPS, counts per second.

 
Small G Protein GTP Hydrolysis Activity—
A fourth assay was developed with GSH-coated FlashPlates to measure GTP hydrolysis activity. GTPases were loaded with [{gamma}-33P]GTP, and the level of radionucleotide associated with the GTPases was measured as a function of time (Fig. 5A). The signal decrease observed is produced from a combination of hydrolysis of the [{gamma}-33P]GTP and dissociation of the nucleotide. To account for the latter, GTP dissociation from the GTPase was measured independently after loading with a non-hydrolyzable GTP analog, [{gamma}-35S]GTP{gamma}S. The pharmacodynamic parameter measured is the half-life of radionucleotide binding to the Rho GTPases. This parameter corresponds to the time required to achieve 50% hydrolysis after signal correction to account for nucleotide dissociation. A prototypical experiment using CDC-42 revealed that this protein has a high hydrolysis activity (half-life = 10.4 min; Fig. 5B). In our assay, CDC-42, RAC-2, CED-10, MIG-2, and mRAC-1 showed comparable hydrolysis activities with half-lives ranging from 12.0 to 14.1 min, whereas RHO-1 and CRP-1 clustered in another group with half-lives of 41.4 and 95.8 min, respectively (Fig. 5C). Our results were consistent with those reported in the literature for several Rho GTPases using conventional filtration assays (Supplemental Table 1 and Ref. 21).



View larger version (21K):
[in this window]
[in a new window]
 
FIG. 5. [{gamma}-33P]GTP hydrolysis activity assay. A, FlashPlate experimental design for the [{gamma}-33P]GTP hydrolysis activity. GST-GTPases were fixed on a GSH FlashPlate and were loaded with either 500 nCi of [{gamma}-35S]GTP{gamma}S ({square}) or 500 nCi of [{gamma}-33P]GTP (•). When the GTPase hydrolyzes the {gamma} phosphate of [{gamma}-33P]GTP, the radioactivity is liberated in the reaction mixture, and ß particles are quenched by the aqueous solution leading to the inactivation of the scintillation reaction. A decrease of signal is then measured. [{gamma}-35S]GTP{gamma}S, a non-hydrolyzable analog of the GTP, is used as a control to measure the nucleotide dissociation from the GST fusion protein. B, the half-life value is the time needed for a GTPase to hydrolyze 50% of the initially bound [{gamma}-33P]GTP with the [{gamma}-35S]GTP{gamma}S dissociation parameter excluded. The graph shows the hydrolysis and dissociation curves obtained for CDC-42 in a single experiment done in triplicate. C, comparison of the half-life values obtained for six C. elegans GST-GTPases (CDC-42, RAC-1, CED-10, RHO-1, MIG-2, and CRP-1) and two mouse GST-GTPases (mRAC-1 and the GTP binding-deficient mRAC-1N17) using the hydrolysis activity assay. The half-life values represent the average of at least three independent experiments performed in triplicate.

 
Computation of Hierarchical Trees Based on Biochemical Properties of Small G Proteins—
To compare and integrate our data, we computed hierarchical clustering trees based on the biochemical characteristic of the proteins tested (Fig. 6). As predicted from our results, the enzymatic tree based on GTP binding and GTP{gamma}S competition underlined the lower affinity of MIG-2 and CRP-1 for GTP{gamma}S. These two GTPases clustered in a separate group from the rest of the RHO GTPases tested (Fig. 6, A and B). Similarly the exchange activity-based tree (Fig. 6C) demonstrated that MIG-2 segregated from the others GTPases, and the hydrolysis activity-based tree showed the partitioning of CRP-1 and RHO-1 from the rest of the GTPases (Fig. 6D). Interestingly these observations correlate with the fact that CRP-1 and RHO-1 are the RHO GTPases with the most deviant P-loop, which is critical for GTP hydrolysis (Fig. 1B).



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 6. Comparison and integration of the enzymatic trees of Rho GTPases. Hierarchical clustering trees based on GTPase biochemical activities are shown. The trees show the distance between each GTPase according to their GTP{gamma}S binding (A), their b-GTP{gamma}S competition (B), their kinetic of exchange (C), their speed of GTP hydrolysis (D), and according to the integration of these four activities (E). These trees were then compared with the amino acid sequence-based tree (F). The bar graph (G) shows the distance, calculated using Euclidean matrices, between the SB tree and different hierarchical trees integrating one, two, three, or four activities (A, B, C, and D). The smaller the distance is, the closer the biochemical tree is from the SB tree.

 
We then hypothesized that integrative trees combining different biochemical characteristics should be representative of the amino acid sequence-based tree in the context of conserved families of proteins, and the sequence-based trees should be proportional to the number of biochemical characteristics included in the integrative trees. When the hierarchical clustering tree combining the four biochemical properties of Rho GTPases (Fig. 6E) was compared with the amino acid sequence-based (SB) tree (Fig. 6F), we found that, as expected, the distance was smaller than when the latter was compared with the trees based on individual biochemical characteristics (Supplemental Fig. 3). In all the analyses, GST was used as a control for tree calculation. For each activity, GST clustered out of the Rho GTPases family validating the process of tree calculation. We then calculated the distance between the SB tree and various combinatorial trees based on biochemical properties of the Rho GTPases (combining one, two, three, or four activities; Fig. 6G). Unexpectedly the smallest distances observed between the SB and biochemical trees were when only two biochemical characteristics were combined (Fig. 6G, B+D and A+D). More specifically, the minimal distance to the SB tree was obtained when the affinity for GTP was combined with the hydrolysis activity, thus reflecting the large number of amino acid residues involved in those two processes. In addition, the fact that the exchange activity parameter led to an increase in distance to the SB tree indicates that the small number of residues involved in this process may not be conserved within the Rho family. Together our results indicate that sequence-based trees only provide a partial classification of the Rho GTPases.


    DISCUSSION
 TOP
 ABSTRACT
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
In this study, we developed four different high throughput enzymatic assays to characterize the biochemical properties of small G proteins. To fulfill the needs created by ORFeome scale studies, we used two well accepted high throughput detection technologies, and we developed the algorithms allowing the integrated analysis of the combined data. The first technology, AlphaScreen, allowed us to measure GTP affinity and exchange activity for the small G proteins with a very high reproducibility, with homogeneity, and at low cost. In addition to allowing the characterization of the affinity of GTPases for different nucleotides, this technology represents a powerful tool to characterize GTPases, such as CRP-1, that have unusual properties toward GTP. However, the AlphaScreen platform was not best suited to measure GTP hydrolysis activity. Therefore, we opted for the FlashPlate technology to measure [{gamma}-35S]GTP{gamma}S binding as well as the [{gamma}-33P]GTP hydrolysis activity. When compared with more traditional filtration methods, FlashPlate provided a high throughput tool to follow in a homogenous manner the hydrolysis activity of small G proteins.

We used AlphaScreen and FlashPlate in combination to cover the entire spectrum of small G protein enzymatic activities. Significant results were obtained for six C. elegans and two mouse GTPases of the Rho family. In addition, comparing the biochemical activities of CED-10 and its mouse ortholog mRAC-1 revealed that these two proteins displayed similar affinity for GTP as well as GDP/GTP exchange and GTP hydrolysis activities (Figs. 2–5) even if they share only 83% identity in their amino acid sequence (Fig. 1B). Finally activities obtained for mRAC-1, CED-10, and CDC-42 were comparable to those described in the literature (5, 16, 1820). These observations taken together confirmed that (i) our novel approach provides results comparable in significance to the data previously reported from other studies, (ii) these results are comparable between orthologous proteins, (iii) our experimental platform allows a high coverage of the enzymatic activities of the small G protein, and (iv) our assays are sensitive and reproducible enough to identify subtle enzymatic differences between members of the Ras superfamily of proteins.

This analytical platform would have been incomplete without the tools enabling an integrated analysis of our experimental data. We computed hierarchical trees reporting amino acid sequence proximity between groups of proteins. Using a similar concept, we built hierarchical clustering trees reporting the distance between protein groups based on their biochemical characteristics. The comparison of the amino acid sequence-based and the biochemistry-based trees allowed us to demonstrate that sequence similarity only partially represents the functional characteristics of small G proteins. In addition, we were able to indicate that the amino acid sequence-based trees reflect mostly the combination of two biochemical characteristics, namely the affinity for GTP and the GTP hydrolysis activity. Thus, our analytical platform allowed us to propose an alternative classification of the Rho GTPases based on their enzymatic characteristics.

In this study, we established and designed a platform combining AlphaScreen and FlashPlate technologies to assess the enzymatic characteristics of small GTPases of the Ras superfamily. The proof of concept validation of these assays was made possible by using the subfamily of Rho GTPases from both murine and C. elegans origins. This platform may present several advantages applicable to academic research, for example to study the relationships between phylogeny and functional characteristics by a parallel study of many different proteins from the same phyla. On the other hand this type of platform may be profitable to the industrial/pharmaceutical sector for the screening of GTPase inhibitors/activators. The targeting of specific characteristics such as exchange or hydrolysis activities may therefore represent a specific application of this platform.


    ACKNOWLEDGMENTS
 
We thank the Chevet laboratory for critical advice.


    FOOTNOTES
 
Received, January 25, 2005

Published, MCP Papers in Press, April 6, 2005, DOI 10.1074/mcp.M500025-MCP200

1 The abbreviations used are: GTP{gamma}S, guanosine 5'-3-O-(thio)triphosphate; b-GDP, biotinylated GDP; SB, sequence-based; b-GTP{gamma}S, biotinylated GTP{gamma}S; PEG, polyethylene glycol; NHS, N-hydroxysuccinimide; UPGMA, unweighted pair-group method arithmetic averaged; mRAC-1, mouse RAC-1. Back

2 See www.ebi.ac.uk/clustalw/; EMBL-EBI, ClustalW submission form. Back

* This work was partially supported by Canadian Institutes for Health Research Grant MOP53357(to E. C.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Back

S The on-line version of this manuscript (available at http://www.mcponline.org) contains supplemental material. Back

§ Recipients of salary awards from the Fonds de la Recherche en Santé du Québec. Back

|| Supported by a grant from the Genome Quebec project to The Cell Map Project (Montreal Proteomics Network). Back

¶¶ To whom correspondence may be addressed. Tel.: 514-937-1010 (ext. 230); E-mail: roger.bosse{at}perkinelmer.com

|||| To whom correspondence may be addressed. Tel.: 514-934-1934 (ext. 35468); Fax: 514-843-1411; E-mail: eric.chevet{at}mcgill.ca


    REFERENCES
 TOP
 ABSTRACT
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Takai, Y., Sasaki, T., and Matozaki, T. (2001) Small GTP-binding proteins.Physiol. Rev. 81, 153 –208

  2. Colicelli, J. (2004) Human RAS superfamily proteins and related GTPases.Sci. STKE 250, RE13

  3. Kontani, K., Tada, M., Ogawa, T., Okai, T., Saito, K., Araki, Y., and Katada, T. (2002) Di-Ras, a distinct subgroup of ras family GTPases with unique biochemical properties.J. Biol. Chem. 277, 41070 –41078

  4. Trahey, M., Milley, R. J., Cole, G. E., Innis, M., Paterson, H., Marshall, C. J., Hall, A., and McCormick, F. (1987) Biochemical and biological properties of the human N-ras p21 protein.Mol. Cell. Biol. 7, 541 –544

  5. Shao, H., Kadono-Okuda, K., Finlin, B. S., and Andres, D. A. (1999) Biochemical characterization of the Ras-related GTPases Rit and Rin.Arch. Biochem. Biophys. 371, 207 –219

  6. Haeusler, L. C., Blumenstein, L., Stege, P. Dvorsky, R., and Ahmadian, M. R. (2003) Comparative functional analysis of the Rac GTPases.FEBS Lett. 555, 556 –560

  7. Luan, C. H., Qiu, S., Finley, J. B., Carson, M., Gray, R. J., Huang, W., Johnson, D., Tsao, J., Reboul, J., Vaglio, P., Hill, D. E., Vidal, M., Delucas, L. J., and Luo, M. (2004) High-throughput expression of C. elegans proteins.Genome Res. 14, 2102 –2110

  8. Brasch, M. A., Hartley, J. L., and Vidal, M. (2004) ORFeome cloning and systems biology: standardized mass production of the parts from the parts-list.Genome Res. 14, 2001 –2009

  9. Beaudet, L., Bedard, J., Breton, B., Mercuri, R. J., and Budarf, M. L. (2001) Homogeneous assays for single-nucleotide polymorphism typing using AlphaScreen.Genome Res. 11, 600 –608

  10. Mullinax, T. R., Henrich, G., Kasila, P., Ahern, D. G., Wenske, E. A., Hou, C., Argentieri, D., and Bembenek, M. E. (1999) Monitoring inositol-specific phospholipase C activity using a phospholipid FlashPlate®.J. Biomol. Screen. 4, 151 –155

  11. Reboul, J., Vaglio, P., Rual, J. F., Lamesch, P., Martinez, M., Armstrong, C. M., Li, S., Jacotot, L., Bertin, N., Janky, R., Moore, T., Hudson, J. R., Jr., Hartley, J. L., Brasch, M. A., Vandenhaute, J., Boulton, S., Endress, G. A., Jenna, S., Chevet, E., Papasotiropoulos, V., Tolias, P. P., Ptacek, J., Snyder, M., Huang, R., Chance, M. R., Lee, H., Doucette-Stamm, L., Hill, D. E., and Vidal, M. (2003) C. elegans ORFeome version 1.1: experimental verification of the genome annotation and resource for proteome-scale protein expression.Nat. Genet. 34, 35 –41

  12. Liu, J. F., Chevet, E., Kebache, S., Lemaitre, G., Barritault, D., Larose, L., and Crepin, M. (1999) Functional Rac-1 and Nck signaling networks are required for FGF-2-induced DNA synthesis in MCF-7 cells.Oncogene 18, 6425 –6433

  13. Braun, P., Hu, Y., Shen, B., Halleck, A., Koundinya, M., Harlow, E., and LaBaer, J. (2002) Proteome-scale purification of human proteins from bacteria.Proc. Natl. Acad. Sci. U. S. A. 99, 2654 –2659

  14. Lamesch, P., Milstein, S., Hao, T., Rosenberg, J., Li, N., Sequerra, R., Bosak, S., Doucette-Stamm, L., Vandenhaute, J., Hill, D. E., and Vidal, M. (2004) C. elegans ORFeome version 3.1: increasing the coverage of ORFeome resources with improved gene predictions.Genome Res. 14, 2064 –2069

  15. Reboul, J., Vaglio, P., Tzellas, N., Thierry-Mieg, N., Moore, T., Jackson, C., Shin-i, T., Kohara, Y., Thierry-Mieg, D., Thierry-Mieg, J., Lee, H., Hitti, J., Doucette-Stamm, L., Hartley, J. L., Temple, G. F., Brasch, M. A., Vandenhaute, J., Lamesch, P. E., Hill, D. E., and Vidal, M. (2001) Open-reading-frame sequence tags (OSTs) support the existence of at least 17,300 genes in C. elegans.Nat. Genet. 27, 332 –336

  16. Menard, L., Tomhave, E., Casey, P. J., Uhing, R. J., Snyderman, R., and Didsbury, J. R. (1992) Rac1, a low-molecular-mass GTP-binding-protein with high intrinsic GTPase activity and distinct biochemical properties.Eur. J. Biochem. 206, 537 –546

  17. Feig, L. A., and Cooper, G. M. (1988) Inhibition of NIH 3T3 cell proliferation by a mutant ras protein with preferential affinity for GDP.Mol. Cell. Biol. 8, 3235 –3243

  18. Jenna, S., Caruso, M.-E., Emadali, A., Nguyen, D. T., Dominguez, M., Li, S., Roy, R., Reboul, R., Vidal, M., Bosse, R., and Chevet, E. (2005) Regulation of membrane trafficking by a novel Cdc42-related protein in C. elegans epithelial cells.Mol. Biol. Cell 16, 1629 –1639

  19. Chen, W., Lim, H. H., and Lim, L. (1993) The CDC42 homologue from Caenorhabditis elegans. Complementation of yeast mutation.J. Biol. Chem. 268, 13280 –13285

  20. Knaus, U. G., Heyworth, P. G., Kinsella, B. T., Curnutte, J. T., and Bokoch, G. M. (1992) Purification and characterization of Rac2. A cytosolic GTP-binding protein that regulates human neutrophil NADPH oxidase.J. Biol. Chem. 267, 23575 –23582

  21. Chen, W., Lim, H. H., and Lim, L. (1993) A new member of the ras superfamily, the rac1 homologue from Caenorhabditis elegans. Cloning and sequence analysis of cDNA, pattern of developmental expression, and biochemical characterization of the protein.J. Biol. Chem. 268, 320 –324





This Article
Abstract
Full Text (PDF)
Supplemental Data
All Versions of this Article:
M500025-MCP200v1
4/7/936    most recent
Submit a response
Purchase Article
View Shopping Cart
Alert me when this article is cited
Alert me when eLetters are posted
Alert me if a correction is posted
Services
Email this article to a friend
Similar articles in this journal
Similar articles in PubMed
Alert me to new issues of the journal
Download to citation manager
Glossary
Copyright Permissions
Google Scholar
Articles by Caruso, M.-E.
Articles by Chevet, E.
Articles citing this Article
PubMed
PubMed Citation
Articles by Caruso, M.-E.
Articles by Chevet, E.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 All ASBMB Journals   Journal of Biological Chemistry 
 Journal of Lipid Research   Biochemistry and Molecular Biology Education