COMBINE analysis of the specificity of binding of Ras proteins to their effectors

The small guanosine triphosphate (GTP)-binding proteins of the Ras family are involved in many cellular pathways leading to cell growth, differentiation, and apoptosis. Understanding the interaction of Ras with other proteins is of importance not only for studying signalling mechanisms but also, because of their medical relevance as targets, for anticancer therapy. To study their selectivity and specificity, which are essential to their signal transfer function, we performed COMparative BINding Energy (COMBINE) analysis for 122 different wild-type and mutant complexes between the Ras proteins, Ras and Rap, and their effectors, Raf and RalGDS. The COMBINE models highlighted the amino acid residues responsible for subtle differences in binding of the same effector to the two different Ras proteins, as well as more significant differences in the binding of the two different effectors (RalGDS and Raf) to Ras. The study revealed that E37, D38, and D57 in Ras are nonspecific hot spots at its effector interface, important for stabilization of both the RalGDS-Ras and Raf-Ras complexes. The electrostatic interaction between a GTP analogue and the effector, either Raf or RalGDS, also stabilizes these complexes. The Raf-Ras complexes are specifically stabilized by S39, Y40, and D54, and RalGDS-Ras complexes by E31 and D33. Binding of a small molecule in the vicinity of one of these groups of amino acid residues could increase discrimination between the Raf-Ras and RalGDS-Ras complexes. Despite the different size of the RalGDS-Ras and Raf-Ras complexes, we succeeded in building COMBINE models for one type of complex that were also predictive for the other type of protein complex. Further, using system-specific models trained with only five complexes selected according to the results of principal component analysis, we were able to predict binding affinities for the other mutants of the particular Ras-effector complex. As the COMBINE analysis method is able to explicitly reveal the amino acid residues that have most influence on binding affinity, it is a valuable aid for protein design.     

Structural and biochemical analysis of Ras-effector signaling via RalGDS.

The structure of the complex of Ras with the Ras-binding domain of its effector RalGDS (RGS-RBD), the first genuine Ras-effector complex, has been solved by X-ray crystallography. As with the Rap-RafRBD complex (Nasser et al., 1995), the interaction is via an inter-protein beta-sheet between the switch I region of Ras and the second strand of the RGS-RBD sheet, but the details of the interactions in the interface are remarkably different. Mutational studies were performed to investigate the contribution of selected interface residues to the binding affinity. Gel filtration experiments show that the Ras x RGS-RBD complex is a monomer. The results are compared to a recently determined structure of a similar complex using a Ras mutant (Huang et al., 1998) and are discussed in relation to partial loss-of-function mutations and the specificity of Ras versus Rap binding.     

Electrostatic fields near the active site of human aldose reductase: 1. New inhibitors and vibrational stark effect measurements

Vibrational Stark effect spectroscopy was used to measure electrostatic fields in the hydrophobic region of the active site of human aldose reductase (hALR2). A new hALR2 inhibitor was designed and synthesized that contains a nitrile probe with a Stark tuning rate of 0.77 cm-1/(MV/cm). Mutations to amino acid residues in the vicinity of the nitrile functional group were selected based on electrostatics calculations, possible complications from hydrogen bonds near the nitrile, and comparison with the active site of human aldehyde reductase, whose structure is very similar. Changes in the absorption energy of the nitrile probe when bound to those mutated proteins were then used to quantify perturbations to the protein's electrostatic field. Electrostatic field changes as large as -10 MV/cm were observed. Measured electrostatic fields were compared to predictions based on continuum electrostatics calculations, revealing that substantial modifications to the calculation strategy are necessary. The effects of hydrogen bonding of amino acid side chains to the nitrile probe are considered, and applications of vibrational Stark effect spectroscopy to investigations of ligand binding and biological function are discussed.     

The activation of RalGDS can be achieved independently of its Ras binding domain. Implications of an activation mechanism in Ras effector specificity and signal distribution.

Small GTPases of the Ras family are major players of signal transduction in eukaryotic cells. They receive signals from a number of receptors and transmit them to a variety of effectors. The distribution of signals to different effector molecules allows for the generation of opposing effects like proliferation and differentiation. To understand the specificity of Ras signaling, we investigated the activation of RalGDS, one of the Ras effector proteins with guanine-nucleotide exchange factor activity for Ral. We determined the GTP level on RalA and showed that the highly conserved Ras binding domain (RBD) of RalGDS, which mediates association with Ras, is important but not sufficient to explain the stimulation of the exchange factor. Although a point mutation in the RBD of RalGDS, which abrogates binding to Ras, renders RalGDS independent to activated Ras, an artificially membrane-targeted version of RalGDS lacking its RBD could still be activated by Ras. The switch II region of Ras is involved in the activation, because the mutant Y64W in this region is impaired in the RalGDS activation. Furthermore, it is shown that Rap1, which was originally identified as a Ras antagonist, can block Ras-mediated RalGDS signaling only when RalGDS contains an intact RBD. In addition, kinetic studies of the complex formation between RalGDS-RBD and Ras suggest that the fast association between RalGDS and Ras, which is analogous to the Ras/Raf case, achieves signaling specificity. Conversely, the Ras x RalGDS complex has a short lifetime of 0.1 s and Rap1 forms a long-lived complex with RalGDS, possibly explaining its antagonistic effect on Ras.     

Electrostatic effects of mutations of Ras glutamine 61 measured using vibrational spectroscopy of a thiocyanate probe

Mutations of human oncoprotein p21(Ras) (hereafter Ras) at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that electrostatically stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. We have tested the interactions between amino acids at this position and water by measuring changes in the electrostatic field experienced by a nitrile probe positioned near Ras Q61 using vibrational Stark effect (VSE) spectroscopy. We mutated this glutamine to every amino acid except cysteine and proline and then incubated these mutants with a Ral guanine nucleotide dissociation stimulator (Ral) containing the I18C mutation that was chemically labeled with a thiocyanate vibrational spectroscopic probe. The formation of the docked Ras Q61X-labeled Ral complex was confirmed by measurement of the dissociation constant of the interaction. We measured the absorption energy of this nitrile to determine any differences in electrostatic environment in the immediate vicinity of the thiocyanate probe between wild type and mutants of Ras. For each Ras Q61X mutant, we correlate the change in electrostatic field at position 61 with the solvent accessible surface area of polar components of the mutant side chain determined from a Boltzmann-weighted ensemble of structures, as well as the residue's hydration potential. These results support the hypothesis that the role of Ras Q61 is to stabilize water in or near the active site during GTP hydrolysis. The substantial effect that nonpolar side chains of Ras Q61X have on the absorption energy of the thiocyanate must be investigated with further experiments.     

A detailed thermodynamic analysis of ras/effector complex interfaces

Many cellular functions are based on the interaction and crosstalk of various signaling proteins. Among these, members of the Ras family of small GTP-binding proteins are important for communicating signals into different pathways. In order to answer the question of how binding affinity and specificity is achieved, we analyzed binding energetics on the molecular level, with reference to the available structural data. The interaction of two members of the Ras subfamily with two different effector proteins, namely Raf and RalGDS, were investigated using isothermal titration calorimetry and a fluorescence-based method. Experiments with alanine mutants, located in the complex interfaces, yielded an energy map for the contact areas of the Ras/effector complexes, which could be differentiated into enthalpy and entropy contributions. In addition, by using double mutant cycle analysis, we probed the energetic contribution of selected pairs of amino acid residues. The resulting energy landscapes of the Ras/effector interface areas show a highly different topology when comparing the two effectors, Raf and RalGDS, demonstrating the specificity of the respective interactions. Particularly, we observe a high degree of compensating effects between enthalpy and entropy; differences between these components are much greater than the overall free energy differences. This is observed also when using the software FOLD-X to predict the effect of point mutations on the crystal structures of the different complexes. Prediction of the free energy changes shows a very good correlation with the experimentally observed energies. Furthermore, in line with experimental data, energy decomposition indicates that many different components of large magnitude counteract each other to produce a smaller change in overall free energy, illustrating the importance of long-range electrostatic forces in complex formation.     

Electrostatic fields near the active site of human aldose reductase: 2. New inhibitors and complications caused by hydrogen bonds

Vibrational Stark effect spectroscopy was used to measure electrostatic fields in the hydrophobic region of the active site of human aldose reductase (hALR2). A new nitrile-containing inhibitor was designed and synthesized, and the X-ray structure of its complex, along with cofactor NADP(+), with wild-type hALR2 was determined at 1.3 ? resolution. The nitrile is found to be in the proximity of T113, consistent with a hydrogen bond interaction. Two vibrational absorption peaks were observed at room temperature in the nitrile region when the inhibitor binds to wild-type hALR2, indicating that the nitrile probe experiences two different microenvironments, and these could be empirically separated into a hydrogen-bonded and non-hydrogen-bonded population by comparison with the T113A mutant, in which a hydrogen bond to the nitrile is not present. Classical molecular dynamics simulations based on the structure predict a double-peak distribution in protein electric fields projected along the nitrile probe. The interpretation of these two peaks as a hydrogen bond formation-dissociation process between the probe nitrile group and a nearby amino acid side chain is used to explain the observation of two IR bands, and the simulations were used to investigate the molecular details of this conformational change. Hydrogen bonding complicates the simplest analysis of vibrational frequency shifts as being due solely to electrostatic interactions through the vibrational Stark effect, and the consequences of this complication are discussed.     

Recognizing and defining true Ras binding domains I: biochemical analysis

Common domain databases contain sequence motifs which belong to the ubiquitin fold family and are called Ras binding (RB) and Ras association (RalGDS/AF6 Ras associating) (RA) domains. The name implies that they bind to Ras (or Ras-like) GTP-binding proteins, and a few of them have been documented to qualify as true Ras effectors, defined as binding only to the activated GTP-bound form of Ras. Here we have expressed a large number of these domains and investigated their interaction with Ras, Rap and M-Ras. While their (albeit weak) sequence homology suggest that the domains adopt a common fold, not all of them bind to Ras proteins, irrespective of whether they are called RB or RA domains. We used fluorescence spectroscopy and isothermal titration calorimetry to show that the binding affinities vary over a large range, and are usually specific for either Ras or Rap. Moreover, the specificity is dictated by a set of key residues in the interface. Stopped-flow kinetic analysis showed that the association rate constants determine the different affinities of effector binding, while the dissociation rate constants are in a similar range. Manual sequence analysis allowed us to define positively charged sequence epitopes in certain secondary structure elements of the ubiquitin fold (beta1, beta2 and alpha1) which are located at similar positions and comprise the hot spots of the binding interface. These residues are important to qualify an RA/RB domain as a true candidate Ras or Rap effector.     

Vibrational Stark effects calibrate the sensitivity of vibrational probes for electric fields in proteins

Infrared spectroscopy is widely used to probe local environments and dynamics in proteins. The introduction of a unique vibration at a specific site of a protein or more complex assembly offers many advantages over observing the spectra of an unmodified protein. We have previously shown that infrared frequency shifts in proteins can arise from differences in the local electric field at the probe vibration. Thus, vibrational frequencies can be used to map electric fields in proteins at many sites or to measure the change in electric field due to a perturbation. The Stark tuning rate gives the sensitivity of a vibrational frequency to an electric field, and for it to be useful, the Stark tuning rate should be as large as possible. Vibrational Stark effect spectroscopy provides a direct measurement of the Stark tuning rate and allows a quantitative interpretation of frequency shifts. We present vibrational Stark spectra of several bond types, extending our work on nitriles and carbonyls and characterizing four additional bond types (carbon-fluorine, carbon-deuterium, azide, and nitro bonds) that are potential probes for electric fields in proteins. The measured Stark tuning rates, peak positions, and extinction coefficients provide the primary information needed to design amino acid analogues or labels to act as probes of local environments in proteins.     

Zinc(II)cyclen-peptide conjugates interacting with the weak effector binding state of Ras

Zinc(II)cyclen-peptide hybrid compounds and bis-zinc(II)cyclen complexes are prepared as potential binders of the guanine nucleotide binding protein Ras, an important molecular switch in cellular signal transduction. The design of the compounds is based on the previous observation that zinc(II)cyclen complexes could serve as lead compounds for inhibitors of Ras-effector interaction and thus be able to interrupt Ras induced signal transduction. Zinc(II)cyclen selectively stabilizes conformational state 1 of active Ras, a conformational state with drastically decreased affinity to effector proteins like Raf-kinase. To achieve higher binding affinities of such Ras-Raf interaction inhibitors, zinc(II)cyclen conjugates with short peptides, derived from the sequence of the Ras-activator SOS, were prepared by solid phase synthesis protocols. Dinuclear bis-zinc(II)cyclen complexes were obtained from alkyne-azide cycloaddition reactions. NMR investigations of the prepared compounds revealed that the peptide conjugates do not lead to an increase in Ras binding affinity of the metal complex-peptide conjugates. The dinuclear zinc complexes lead to an immediate precipitation of the protein prohibiting spectroscopic investigations of their binding.     

Double-mutant analysis of the interaction of Ras with the Ras-binding domain of RGL

RalGDS is a guanine nucleotide dissociation stimulator for Ral, and one of its homologues is RGL (RalGDS-like). In this study, the effects of mutations of Ras and the Ras-binding domains (RBDs) of RalGDS and RGL on their binding have been systematically examined. The D33A mutation of Ras reduces the abilities to bind RGL-RBD and RalGDS-RBD. To identify the RGL residue interacting with Asp33 of Ras, double-mutant analyses between Ras and RGL-RBD were conducted. For example, the K685A mutation of RGL-RBD has a much smaller effect on the RGL-RBD binding ability of the D33A mutant than on those of other mutants of Ras. Accordingly, it is indicated that the attractive interaction of Asp33 in Ras with Lys685 in RGL-RBD (Lys816 in RalGDS-RBD) contributes to the Ras.RBD association. This interaction is consistent with the crystal structure of the complex of RalGDS-RBD and the E31K Ras mutant [Huang, L., Hofer, F., Martin, G. S., and Kim, S.-H. (1998) Nat. Struct. Biol. 5, 422-426]. This crystal structure exhibits interactions of the mutation-derived Lys31 side chain with three RalGDS residues. Glu31 of Ras discriminates Ras from a Ras-homologue, Rap1, with Lys31, with respect to RalGDS and RGL binding; the E31K mutation of Ras potentiates the abilities to bind RGL-RBD and RalGDS-RBD. To examine the role of Glu31 of the wild-type Ras in the interaction with RGL and RalGDS, double-mutant analyses were conducted. The Ras binding ability of the E689A mutant of RGL-RBD is much stronger than that of the wild-type RGL-RBD, and the E31K mutation of Ras no longer potentiates the Ras binding ability of the E689A mutant. Therefore, the repulsive interaction between Glu31 in Ras and Glu689 in RGL-RBD (Asp820 in RalGDS-RBD) may keep the Ras.RBD association weaker than the Rap1.RBD association, which might be relevant to the regulation of the signaling network.     

Conformational states of Ras complexed with the GTP analogue GppNHp or GppCH2p: implications for the interaction with effector proteins

The guanine nucleotide-binding protein Ras occurs in solution in two different states, state 1 and state 2, when the GTP analogue GppNHp is bound to the active center as detected by (31)P NMR spectroscopy. Here we show that Ras(wt).Mg(2+).GppCH(2)p also exists in two conformational states in dynamic equilibrium. The activation enthalpy DeltaH(++)(12) and the activation entropy DeltaS(++)(12) for the transition from state 1 to state 2 are 70 kJ mol(-1) and 102 J mol(-1) K(-1), within the limits of error identical to those determined for the Ras(wt).Mg(2+).GppNHp complex. The same is true for the equilibrium constants K(12) = [2]/[1] of 2.0 and the corresponding DeltaG(12) of -1.7 kJ mol(-1) at 278 K. This excludes a suggested specific effect of the NH group of GppNHp on the equilibrium. The assignment of the phosphorus resonance lines of the bound analogues has been done by two-dimensional (31)P-(31)P NOESY experiments which lead to a correction of the already reported assignments of bound GppNHp. Mutation of Thr35 in Ras.Mg(2+).GppCH(2)p to serine leads to a shift of the conformational equilibrium toward state 1. Interaction of the Ras binding domain (RBD) of Raf kinase or RalGDS with Ras(wt) or Ras(T35S) shifts the equilibrium completely to state 2. The (31)P NMR experiments suggest that, besides the type of the side chain of residue 35, a main contribution to the conformational equilibrium in Ras complexes with GTP and GTP analogues is the effective acidity of the gamma-phosphate group of the bound nucleotide. A reaction scheme for the Ras-effector interaction is presented which includes the existence of two conformations of the effector loop and a weak binding state.     

Biochemical characterization of the Ras-related GTPases Rit and Rin.

We report the biochemical characterization of Rit and Rin, two members of the Ras superfamily identified by expression cloning. Recombinant Rit and Rin bind GTP and exhibit intrinsic GTPase activity. Conversion of Gln to Leu at position 79 (for Rit) or 78 (for Rin) (equivalent to position 61 in Ras) resulted in a complete loss of GTPase activity. Surprisingly, significant differences were found when the guanine nucleotide dissociation constants of Rit and Rin were compared with the majority of Ras-related GTPases. Both proteins display higher k(off) values for GTP than GDP in the presence of 10 mM Mg(2+). These GTP dissociation rates are 5- to 10-fold faster than most Ras-like GTPases. Despite these unique biochemical properties, our data support the notion that both Rit and Rin function as nucleotide-dependent molecular switches. To begin to address whether these proteins act as regulators of distinct signaling pathways, we examined their interaction with a series of known Ras-binding proteins by yeast two-hybrid analysis. Although Rit, Rin, and Ras have highly related effector domain sequences, Rit and Rin were found to interact with the known Ras binding proteins RalGDS, Rlf, and AF-6/Canoe but not with the Raf kinases, RIN1, or the p110 subunit of phosphatidylinositol 3-kinase. These interactions were GTP and effector domain dependent and suggest that RalGDS, Rlf, and AF-6 are Rit and Rin effectors. Their biochemical properties and interaction with a subset of known Ras effector proteins suggest that Rit and Rin may play important roles in the regulation of signaling pathways and cellular processes distinct from those controlled by Ras.     

Mechanism of Interaction between the General Anesthetic Halothane and a Model Ion Channel Protein, III: Molecular Dynamics Simulation Incorporating a Cyanophenylalanine Spectroscopic Probe

A nitrile-derived amino acid, Phe_CN, has been used as an internal spectroscopic probe to study the binding of an inhalational anesthetic to a model membrane protein. The infrared spectra from experiment showed a blue-shift of the nitrile vibrational frequency in the presence of the anesthetic halothane. To interpret the infrared results and explore the nature of the interaction between halothane and the model protein, all-atom molecular dynamics (MD) simulations have been used to probe the structural and dynamic properties of the protein in the presence and absence of one halothane molecule. The frequency shift analyzed from MD simulations agrees well with the experimental infrared results. Decomposition of the forces acting on the nitrile probes demonstrates an indirect impact on the probes from halothane, namely a change of the protein's electrostatic local environment around the probes induced by halothane. Although the halothane remains localized within the designed hydrophobic binding cavity, it undergoes a significant amount of translational and rotational motion, modulated by the interaction of the trifluorine end of halothane with backbone hydrogens of the residues forming the cavity. This dominant interaction between halothane and backbone hydrogens outweighs the direct interaction between halothane and the nitrile groups, making it a good “spectator” probe of the halothane-protein interaction. These MD simulations provide insight into action of anesthetic molecules on the model membrane protein, and also support the further development of nitrile-labeled amino acids as spectroscopic probes within the designed binding cavity.     

The Ras-Byr2RBD complex: structural basis for Ras effector recognition in yeast.

BACKGROUND: The small GTP binding protein Ras has important roles in cellular growth and differentiation. Mutant Ras is permanently active and contributes to cancer development. In its activated form, Ras interacts with effector proteins, frequently initiating a kinase cascade. In the lower eukaryotic Schizosaccharomyces pombe, Byr2 kinase represents a Ras target that in terms of signal-transduction hierarchy can be considered a homolog of mammalian Raf-kinase. The activation mechanism of protein kinases by Ras is not understood, and there is no detailed structural information about Ras binding domains (RBDs) in nonmammalian organisms. RESULTS: The crystal structure of the Ras-Byr2RBD complex at 3 A resolution shows a complex architecture similar to that observed in mammalian homologous systems, with an interprotein beta sheet stabilized by predominantly polar interactions between the interacting components. The C-terminal half of the Ras switch I region contains most of the contact anchors, while on the Byr2 side, a number of residues from topologically distinct regions are involved in complex stabilization. A C-terminal helical segment, which is not present in the known mammalian homologous systems and which is part of the auto-inhibitory region, has an additional binding site outside the switch I region. CONCLUSIONS: The structure of the Ras-Byr2 complex confirms the Ras binding module as a communication element mediating Ras-effector interactions; the Ras-Byr2 complex is also conserved in a lower eukaryotic system like yeast, which is in contrast to other small GTPase families. The extra helical segment might be involved in kinase activation.     

JAK/STAT3-dependent activation of the RalGDS/Ral pathway in M1 mouse myeloid leukemia cells.

The Ras-related GTPase (Ral) is converted to the GTP-bound form by Ral guanine nucleotide dissociation stimulator (RalGDS), a putative effector protein of Ras. Recently, it was proven that Ral regulates c-Src activity and subsequent phosphorylation of its substrate, STAT3. Here, we show that STAT3 inversely regulates activation of Ral through induction of expression of RalGDS. To identify new leukemia inhibitory factor-induced genes, we have performed representational difference analysis using M1 mouse myeloid leukemia cells and cloned RalGDS. The expression of RalGDS and subsequent activation of RalA were clearly suppressed by a dominant negative form of STAT3 and a JAK inhibitor, JAB/SOCS1/SSI-1, indicating that RalGDS/RalA signaling requires the activation of the JAK/STAT3 pathway. An experiment using a Ras inhibitor demonstrated that full activation of RalA also requires activation of Ras. These results suggest a novel cross-talk between JAK/STAT3 and the Ras/RalGDS/Ral signaling pathways through gp130.     

Display of Ras on filamentous phage through cysteine replacement

Phage display technology has been used in a variety of contexts to understand and manipulate biomolecular interactions between proteins and other biomolecules. In this paper we describe the establishment of a phage display system for elucidation of the interactions between the GTPase Ras and its panel of effectors. It is shown how technical problems associated with phage display of a protein with unpaired cysteines, likely to be caused by the oxidizing environment of the bacterial periplasm into which the protein is directed, can be overcome by cysteine replacement based on functional and structural studies. First, the catalytic domain (residues 1-166) of mammalian H-Ras (Ras) was observed to be displayed on phage in an incorrect conformation not detectable by antibodies recognizing conformational epitopes on Ras. Although truncation of the phage coat protein used as fusion partner (g3p) resulted in minor improvements in the display, Ras was tailored for phage display by cysteine replacement. By replacing the three cysteines at positions 51, 80 and 118 of Ras with the corresponding residues in Saccharomyces cerevisiae RAS1, the resulting fusion-phage is recognized by the conformation-dependent anti-Ras antibodies. Furthermore, display of cysteine-free Ras is demonstrated by GTP-analogue dependent binding to the Ras-binding domain of the Ras-effector Raf1. These data pave the way for analysis of Ras-effector interactions using phage display technology yet demonstrate that phage display of proteins with normally reduced cysteines should be approached with caution.     

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.     

Prediction of Ras-effector interactions using position energy matrices

MOTIVATION: One of the more challenging problems in biology is to determine the cellular protein interaction network. Progress has been made to predict protein-protein interactions based on structural information, assuming that structural similar proteins interact in a similar way. In a previous publication, we have determined a genome-wide Ras-effector interaction network based on homology models, with a high accuracy of predicting binding and non-binding domains. However, for a prediction on a genome-wide scale, homology modelling is a time-consuming process. Therefore, we here successfully developed a faster method using position energy matrices, where based on different Ras-effector X-ray template structures, all amino acids in the effector binding domain are sequentially mutated to all other amino acid residues and the effect on binding energy is calculated. Those pre-calculated matrices can then be used to score for binding any Ras or effector sequences. RESULTS: Based on position energy matrices, the sequences of putative Ras-binding domains can be scanned quickly to calculate an energy sum value. By calibrating energy sum values using quantitative experimental binding data, thresholds can be defined and thus non-binding domains can be excluded quickly. Sequences which have energy sum values above this threshold are considered to be potential binding domains, and could be further analysed using homology modelling. This prediction method could be applied to other protein families sharing conserved interaction types, in order to determine in a fast way large scale cellular protein interaction networks. Thus, it could have an important impact on future in silico structural genomics approaches, in particular with regard to increasing structural proteomics efforts, aiming to determine all possible domain folds and interaction types. AVAILABILITY: All matrices are deposited in the ADAN database (http://adan-embl.ibmc.umh.es/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

RalGDS couples growth factor signaling to Akt activation

The Akt kinase is a key regulator of cell proliferation and survival. It is activated in part by PDK1-induced phosphorylation. Here we show that RalGDS, a Ras effector protein that activates Ral GTPases, has a second function that promotes Akt phosphorylation by PDK1 by bringing these two kinases together. In support of this conclusion is our finding that suppression of RalGDS expression in cells inhibits both epidermal growth factor and insulin-induced phosphorylation of Akt. Moreover, while PDK1 complexes with N-GDS, Akt complexes with the central region of RalGDS through an intermediary, JIP1. The biological significance of this newly discovered RalGDS function is highlighted by the observation that an N-terminally deleted mutant of RalGDS that retains the ability to activate Ral proteins but loses the ability to activate Akt also fails to promote cell proliferation. Thus, RalGDS forms a nexus that transduces growth factor signaling to both Ral GTPase and Akt-mediated signaling cascades.