Solution structure of the Ras binding domain of the protein kinase Byr2 from Schizosaccharomyces pombe

Mutational analyses of the interacting domains of Schizosaccharomyces pombe Byr2 with 14-3-3s

The Byr2 kinase of fission yeast Schizosaccharomyces pombe is recruited to the membrane with the assistance of Ras1. Byr2 is also negatively regulated by 14-3-3 proteins encoded by rad24 and rad25. We conducted domain and mutational analysis of Byr2 to determine which region is critical for its binding to 14-3-3 proteins. Rad24 and Rad25 bound to both the Ras interaction domain in the N-terminus and to the C-terminal catalytic domain of Byr2. When amino acid residues S87 and T94 of the Ras-interacting domain of Byr2 were mutated to alanine, Rad24 could no longer bind to Byr2. S402, S566, S650, and S654 mutations in the C-terminal domain of Byr2 also abolished its interaction with Rad24 and Rad25. More than three mutations in the C-terminal domain were required to abolish completely its interaction with 14-3-3 protein, suggesting that multiple residues are involved in this interaction. Expression of the N-terminal domain of Byr2 in wild-type cells lowered the mating ratio, because it likely blocked the interaction of Byr2 with Ste4 and Ras1, whereas expression of the catalytic domain of Byr2 increased the mating ratio as a result of freeing from intramolecular regulation by the N-terminal domain of Byr2. The S87A and T94A mutations of Byr2 increased the mating ratio and attenuated inhibition of Byr2 by Rad24; therefore, these two amino acids are critical for its regulation by Rad24. S566 of Byr2 is critical for activity of Byr2 but not for its interaction with 14-3-3 proteins. In this study, we show that 14-3-3 proteins interact with two separate domains in Byr2 as negative regulators.

E. coli elongation factor Tu bound to a GTP analogue displays an open conformation equivalent to the GDP-bound form

According to the traditional view, GTPases act as molecular switches, which cycle between distinct ‘on’ and ‘off’ conformations bound to GTP and GDP, respectively. Translation elongation factor EF-Tu is a GTPase essential for prokaryotic protein synthesis. In its GTP-bound form, EF-Tu delivers aminoacylated tRNAs to the ribosome as a ternary complex. GTP hydrolysis is thought to cause the release of EF-Tu from aminoacyl-tRNA and the ribosome due to a dramatic conformational change following P_i release. Here, the crystal structure of Escherichia coli EF-Tu in complex with a non-hydrolysable GTP analogue (GDPNP) has been determined. Remarkably, the overall conformation of EF-Tu·GDPNP displays the classical, open GDP-bound conformation. This is in accordance with an emerging view that the identity of the bound guanine nucleotide is not ‘locking’ the GTPase in a fixed conformation. Using a single-molecule approach, the conformational dynamics of various ligand-bound forms of EF-Tu were probed in solution by fluorescence resonance energy transfer. The results suggest that EF-Tu, free in solution, may sample a wider set of conformations than the structurally well-defined GTP- and GDP-forms known from previous X-ray crystallographic studies. Only upon binding, as a ternary complex, to the mRNA-programmed ribosome, is the well-known, closed GTP-bound conformation, observed.

Extending the Lifetime of Native GTP-Bound Ras for Site-Resolved NMR Measurements: Quantifying the Allosteric Dynamics

Characterization of native GTP-bound Ras is important for an appreciation of its cellular signaling and for the design of inhibitors, which however has been depressed by its intrinsic instability. Herein, an effective approach for extending the lifetime of Ras⋅GTP samples by exploiting the active role of Son of Sevenless (Sos) is demonstrated that sustains the activated state of Ras. This approach, combined with a postprocessing method that suppresses residual Ras⋅GDP signals, is applied to the site-resolved NMR measurement of the allosteric dynamics of Ras⋅GTP. The observed network of concerted motions well covers the recently identified allosteric inhibitor-binding pockets, but the motions are more confined than those of Ras⋅GppNHp, advocating the use of native GTP for development of allosteric inhibitors. The Sos-based approach is anticipated to generally facilitate experiments on active Ras when native GTP is preferred.

Inhibition of Ras activity coordinates cell fusion with cell-cell contact during yeast mating

In the fission yeast Schizosaccharomyces pombe, pheromone signaling engages a signaling pathway composed of a G protein-coupled receptor, Ras, and a mitogen-activated protein kinase (MAPK) cascade that triggers sexual differentiation and gamete fusion. Cell-cell fusion requires local cell wall digestion, which relies on an initially dynamic actin fusion focus that becomes stabilized upon local enrichment of the signaling cascade on the structure. We constructed a live-reporter of active Ras1 (Ras1-guanosine triphosphate [GTP]) that shows Ras activity at polarity sites peaking on the fusion structure before fusion. Remarkably, constitutive Ras1 activation promoted fusion focus stabilization and fusion attempts irrespective of cell pairing, leading to cell lysis. Ras1 activity was restricted by the guanosine triphosphatase-activating protein Gap1, which was itself recruited to sites of Ras1-GTP and was essential to block untimely fusion attempts. We propose that negative feedback control of Ras activity restrains the MAPK signal and couples fusion with cell-cell engagement.

Activation of Mst11 and Feedback Inhibition of Germ Tube Growth in Magnaporthe oryzae

Appressorium formation and invasive growth are two important steps in the infection cycle of Magnaporthe oryzae that are regulated by the Mst11-Mst7-Pmk1 mitogen-activated protein kinase (MAPK) pathway. However, the molecular mechanism involved in the activation of Mst11 MAPK kinase kinase is not clear in the rice blast fungus. In this study, we functionally characterized the regulatory region of Mst11 and its self-inhibitory binding. Deletion of the middle region of Mst11, which contains the Ras-association (RA) domain and two conserved phosphorylation sites (S453 and S458), blocked Pmk1 activation and appressorium formation. However, the MST11(ΔRA) transformant MRD-2 still formed appressoria, although it was reduced in virulence. Interestingly, over 50% of its germ tubes branched and formed two appressoria by 48 h, which was suppressed by treatments with exogenous cAMP. The G18V dominant active mutation enhanced the interaction of Ras2 with Mst11, suggesting that Mst11 has stronger interactions with the activated Ras2. Furthermore, deletion and site-directed mutagenesis analyses indicated that phosphorylation at S453 and S458 of Mst11 is important for appressorium formation and required for the activation of Pmk1. We also showed that the N-terminal region of Mst11 directly interacted with its kinase domain, and the S789G mutation reduced their interactions. Expression of the MST11(S789G) allele rescued the defect of the mst11 mutant in plant infection and resulted in the formation of appressoria on hydrophilic surfaces, suggesting the gain-of-function effect of the S789G mutation. Overall, our results indicate that the interaction of Mst11 with activated Ras2 and phosphorylation of S453 and S458 play regulatory roles in Mst11 activation and infection-related morphogenesis, possibly by relieving its self-inhibitory interaction between its N-terminal region and the C-terminal kinase domain. In addition, binding of Mst11 to Ras2 may be involved in the feedback inhibition of cAMP signaling and further differentiation of germ tubes after appressorium formation.

Elucidating the mode of action of a typical Ras state 1(T) inhibitor

The small GTPase Ras is an essential component of signal transduction pathways within the cell, controlling proliferation, differentiation, and apoptosis. Only in the GTP-bound form does Ras interact strongly with effector molecules such as Raf-kinase, thus acting as a molecular switch. In the GTP-bound form, Ras exists in a dynamic equilibrium between at least two distinct conformational states, 1(T) and 2(T), offering different functional properties of the protein. Zn2+-cyclen is a typical state 1(T) inhibitor; i.e., it interacts selectively with Ras in conformational state 1(T), a weak effector binding state. Here we report that active K-Ras4B, which is prominently found to be mutated in human tumors, exhibits a dynamic equilibrium like H-Ras, which can be modulated by Zn2+-cyclen. The titration experiments of Ras with Zn2+-cyclen indicate a cooperatively coupled binding of the ligands to the two interaction sites on Ras that could be identified for H-Ras previously. Our data further indicate that as in state 2(T) where induced fit produces the substate 2(T)* after effector binding, a corresponding substate 1(T)* can be detected at the state 1(T) mutant Ras(T35A). The interaction of Zn2+-cyclen with Ras not only shifts the equilibrium toward the weak effector binding state 1(T) but also perturbs the formation of substate 1(T)*, thus enhancing the inhibitory effect. Although Zn2+-cyclen shows an affinity for Ras in only the millimolar range, its potency of inhibition corresponds to a competitive state 2 inhibitor with micromolar binding affinity. Thus, the results demonstrate the mode of action and potency of this class of allosteric Ras inhibitors.

The IQGAP-related protein DGAP1 mediates signaling to the actin cytoskeleton as an effector and a sequestrator of Rac1 GTPases

Proteins are typically categorized into protein families based on their domain organization. Yet, evolutionarily unrelated proteins can also be grouped together according to their common functional roles. Sequestering proteins constitute one such functional class, acting as macromolecular buffers and serving as an intracellular reservoir ready to release large quantities of bound proteins or other molecules upon appropriate stimulation. Another functional protein class comprises effector proteins, which constitute essential components of many intracellular signal transduction pathways. For instance, effectors of small GTP-hydrolases are activated upon binding a GTP-bound GTPase and thereupon participate in downstream interactions. Here we describe a member of the IQGAP family of scaffolding proteins, DGAP1 from Dictyostelium, which unifies the roles of an effector and a sequestrator in regard to the small GTPase Rac1. Unlike classical effectors, which bind their activators transiently leading to short-lived signaling complexes, interaction between DGAP1 and Rac1-GTP is stable and induces formation of a complex with actin-bundling proteins cortexillins at the back end of the cell. An oppositely localized Rac1 effector, the Scar/WAVE complex, promotes actin polymerization at the cell front. Competition between DGAP1 and Scar/WAVE for the common activator Rac1-GTP might provide the basis for the oscillatory re-polarization typically seen in randomly migrating Dictyostelium cells. We discuss the consequences of the dual roles exerted by DGAP1 and Rac1 in the regulation of cell motility and polarity, and propose that similar signaling mechanisms may be of general importance in regulating spatiotemporal dynamics of the actin cytoskeleton by small GTPases.

Structurally unique interaction of RBD-like and PH domains is crucial for yeast pheromone signaling

The Ras-binding domain is conserved among fungal Ste11 MAPKKKs and is critical for mating in fungi. Its interaction with Ras1 is critical for Schizosaccharomyces pombe mating, whereas in Saccharomyces cerevisiae its interaction with the Ste5 PH domain plays the crucial role. The binding partner of RBD for fungal mating is shifted from Ras to a PH domain in fungi in which Ste5 exists.

The Ste5 protein forms a scaffold that associates and regulates the components of the mitogen-activated protein (MAP) kinase cascade that controls mating-pheromone-mediated signaling in the yeast Saccharomyces cerevisiae. Although it is known that the MEK kinase of the pathway, Ste11, associates with Ste5, details of this interaction have not been established. We identified a Ras-binding-domain-like (RBL) region in the Ste11 protein that is required specifically for the kinase to function in the mating pathway. This module is structurally related to domains in other proteins that mediate Ras-MAP kinase kinase kinase associations; however, this RBL module does not interact with Ras, but instead binds the PH domain of the Ste5 scaffold. Structural and functional studies suggest that the key role of this PH domain is to mediate the Ste5–Ste11 interaction. Overall these two evolutionarily conserved modules interact with each other through a unique interface, and thus in the pheromone pathway the structural context of the RBL domain contribution to kinase activation has been shifted through a change of its interaction partner from Ras to a PH domain.

Specific Conformational States of Ras GTPase upon Effector Binding

To uncover the structural and dynamical determinants involved in the highly specific binding of Ras GTPase to its effectors, the conformational states of Ras in uncomplexed form and complexed to the downstream effectors Byr2, PI3Kγ, PLCε, and RalGDS were investigated using molecular dynamics and cross-comparison of the trajectories. The subtle changes in the dynamics and conformations of Ras upon effector binding require an analysis that targets local changes independent of global motions. Using a structural alphabet, a computational procedure is proposed to quantify local conformational changes. Positions detected by this approach were characterized as either specific for a particular effector, specific for an effector domain type, or as effector unspecific. A set of nine structurally connected residues (Ras residues 5–8, 32–35, 39–42, 55–59, 73–78, and 161–165), which link the effector binding site to the distant C-terminus, changed dynamics upon effector binding, indicating a potential effector-unspecific signaling route within the Ras structure. Additional conformational changes were detected along the N-terminus of the central β-sheet. Besides the Ras residues at the effector interface (e.g., D33, E37, D38, and Y40), which adopt effector-specific local conformations, the binding signal propagates from the interface to distant hot-spot residues, in particular to Y5 and D57. The results of this study reveal possible conformational mechanisms for the stabilization of the active state of Ras upon downstream effector binding and for the structural determinants responsible for effector specificity.

Ras signaling in yeast

Since the study of yeast RAS and adenylate cyclase in the early 1980s, yeasts including budding and fission yeasts contributed significantly to the study of Ras signaling. First, yeast studies provided insights into how Ras activates downstream signaling pathways. Second, yeast studies contributed to the identification and characterization of GAP and GEF proteins, key regulators of Ras. Finally, the study of yeast provided many important insights into the understanding of C-terminal processing and membrane association of Ras proteins.

Novel RasGRF1-derived Tat-fused peptides inhibiting Ras-dependent proliferation and migration in mouse and human cancer cells

Mutations of RAS genes are critical events in the pathogenesis of different human tumors and Ras proteins represent a major clinical target for the development of specific inhibitors to use as anticancer agents. Here we present RasGRF1-derived peptides displaying both in vitro and in vivo Ras inhibitory properties. These peptides were designed on the basis of the down-sizing of dominant negative full-length RasGRF1 mutants. The over-expression of these peptides can revert the phenotype of K-RAS transformed mouse fibroblasts to wild type, as monitored by several independent biological readouts, including Ras-GTP intracellular levels, ERK activity, morphology, proliferative potential and anchorage independent growth. Fusion of the RasGRF1-derived peptides with the Tat protein transduction domain allows their uptake into mammalian cells. Chemically synthesized Tat-fused peptides, reduced to as small as 30 residues on the basis of structural constraints, retain Ras inhibitory activity. These small peptides interfere in vitro with the GEF catalyzed nucleotide dissociation and exchange on Ras, reduce cell proliferation of K-RAS transformed mouse fibroblasts, and strongly reduce Ras-dependent IGF-I-induced migration and invasion of human bladder cancer cells. These results support the use of RasGRF1-derived peptides as model compounds for the development of Ras inhibitory anticancer agents.

Cu2+-cyclen as probe to identify conformational states in guanine nucleotide binding proteins

(31)P NMR spectroscopy is a suitable method for identifying conformational states in the active site of guanine nucleotide binding proteins detecting the nucleotide placed there. Because there is no labeling necessary, this method is gaining increasing interest. By (31)P NMR spectroscopy two major conformational states, namely state 1(T) and state 2(T), can be detected in active Ras protein characterized by different chemical shifts. Depending on the conformational state Ras shows clearly different physiological properties. Meanwhile analogous conformational equilibria could also be shown for other members of the Ras superfamily. It is often difficult to determine the conformational states of the proteins on the basis of chemical shift alone; therefore, direct detection would be a great advantage. With the use of Cu(2+)-cyclen which selectively interacts only with one of the major conformational states (state 1) one has a probe to distinguish between the two states, because only proteins existing in conformational state 1 interact with the Cu(2+)-cyclen at low millimolar concentrations. The suitability was proven using Ras(wt) and Ras mutants, Ras complexed with GTP, GppNHp, or GTPγS, as well as two further members of the Ras superfamily namely Arf1 and Ran.

Conformational states of human rat sarcoma (Ras) protein complexed with its natural ligand GTP and their role for effector interaction and GTP hydrolysis

The guanine nucleotide-binding protein Ras exists in solution in two different conformational states when complexed with different GTP analogs such as GppNHp or GppCH(2)p. State 1 has only a very low affinity to effectors and seems to be recognized by guanine nucleotide exchange factors, whereas state 2 represents the high affinity effector binding state. In this work we investigate Ras in complex with the physiological nucleoside triphosphate GTP. By polarization transfer (31)P NMR experiments and effector binding studies we show that Ras(wt)·Mg(2+)·GTP also exists in a dynamical equilibrium between the weakly populated conformational state 1 and the dominant state 2. At 278 K the equilibrium constant between state 1 and state 2 of C-terminal truncated wild-type Ras(1-166) K(12) is 11.3. K(12) of full-length Ras is >20, suggesting that the C terminus may also have a regulatory effect on the conformational equilibrium. The exchange rate (k(ex)) for Ras(wt)·Mg(2+)·GTP is 7 s(-1) and thus 18-fold lower compared with that found for the Ras·GppNHp complex. The intrinsic GTPase activity substantially increases after effector binding for the switch I mutants Ras(Y32F), (Y32R), (Y32W), (Y32C/C118S), (T35S), and the switch II mutant Ras(G60A) by stabilizing state 2, with the largest effect on Ras(Y32R) with a 13-fold increase compared with wild-type. In contrast, no acceleration was observed in Ras(T35A). Thus Ras in conformational state 2 has a higher affinity to effectors as well as a higher GTPase activity. These observations can be used to explain why many mutants have a low GTPase activity but are not oncogenic.

Protein structure calculation with data imputation: the use of substitute restraints

The amount of experimental restraints e.g., NOEs is often too small for calculating high quality three-dimensional structures by restrained molecular dynamics. Considering this as a typical missing value problem we propose here a model based data imputation technique that should lead to an improved estimation of the correct structure. The novel automated method implemented in AUREMOL makes a more efficient use of the experimental information to obtain NMR structures with higher accuracy. It creates a large set of substitute restraints that are used either alone or together with the experimental restraints. The new approach was successfully tested on three examples: firstly, the Ras-binding domain of Byr2 from Schizosaccharomyces pombe, the mutant HPr (H15A) from Staphylococcus aureus, and a X-ray structure of human ubiquitin. In all three examples, the quality of the resulting final bundles was improved considerably by the use of additional substitute restraints, as assessed quantitatively by the calculation of RMSD values to the "true" structure and NMR R-factors directly calculated from the original NOESY spectra or the published diffraction data.

Ca2+-dependent conformational changes in a C-terminal cytosolic domain of polycystin-2

The PKD1 and PKD2 genes are the genes that are mutated in patients suffering from autosomal dominant polycystic kidney disease. The human PKD2 gene codes for a 968-amino acid long membrane protein called polycystin-2 that represents a cation channel whose activity can be regulated by Ca(2+) ions. By CD, fluorescence, and NMR spectroscopy, we have studied a 117-amino acid-long fragment of the cytoplasmic domain of polycystin-2, polycystin-2-(680-796) that was proposed to contain a Ca(2+)-binding site. NMR structure determination reveals the existence of two Ca(2+)-binding sites in polycystin-2-(680-796) arranged in a typical and an atypical EF-hand motif. In the absence of Ca(2+) the protein forms a dimer that is dissociated by Ca(2+) binding. This dissociation may be related to the Ca(2+) inactivation observed earlier. The calcium affinity of the protein was determined by fluorescence and NMR spectroscopy. At 293 K, the K(D) values for the high and low affinity sites are 55 mum and 179 mum, respectively.

NMR-assignments of a cytosolic domain of the C-terminus of polycystin-2

Mutations in the PKD2 gene lead to the development of polycystic kidney disease (PKD). The PKD2 gene codes for polycystin-2, a cation channel with unknown function. The cytoplasmic, C-terminal domain interacts with a large number of proteins including mDia1, alpha-actinin, PIGEA-14, troponin, and tropomyosin. The C-terminal fragment polycystin-2 (680-796) consisting of 117 amino acids contains a putative calcium binding EF-hand. It was produced in Escherichia coli and enriched uniformly with (13)C and (15)N. The backbone and side chain resonances were assigned by multidimensional NMR methods, the obtained chemical shifts are typical for a partially folded protein. The chemical shifts obtained are in line with the existence of two paired helix-loop-helix (HLH) motifs.

Association rate constants of ras-effector interactions are evolutionarily conserved

Evolutionary conservation of protein interaction properties has been shown to be a valuable indication for functional importance. Here we use homology interface modeling of 10 Ras-effector complexes by selecting ortholog proteins from 12 organisms representing the major eukaryotic branches, except plants. We find that with increasing divergence time the sequence similarity decreases with respect to the human protein, but the affinities and association rate constants are conserved as predicted by the protein design algorithm, FoldX. In parallel we have done computer simulations on a minimal network based on Ras-effector interactions, and our results indicate that in the absence of negative feedback, changes in kinetics that result in similar binding constants have strong consequences on network behavior. This, together with the previous results, suggests an important biological role, not only for equilibrium binding constants but also for kinetics in signaling processes involving Ras-effector interactions. Our findings are important to take into consideration in system biology approaches and simulations of biological networks.

Cellular signal transductions processes are based on protein interactions. Proteins can either associate transiently with each other or form stable complexes, and the strength of the interaction is described by the affinity (the affinity is the ratio between the rate of dissociation and association). Protein complexes with similar affinities can bind and dissociate with different rates, and these rates describe the kinetic properties of protein binding. These kinetic rates are important for signaling; however, to what extent individual changes in such rate constants are biologically important or whether the affinity is more crucial might be different in different signaling processes. In this study we analyze whether association rates are conserved during evolution, because evolutionary conservation of protein biochemical properties is usually a valuable indication of its importance. We analyzed the binding of Ras proteins to effector domains, which are central proteins in many signal transduction pathways, in different organisms. On the basis of homology modeling and energy calculations we find that association rates are conserved, although the sequence similarity decreases compared to the human protein. Our finding should encourage further analysis of the importance of kinetics for cellular signal transduction.

Linking Ras to myosin function: RasGEF Q, a Dictyostelium exchange factor for RasB, affects myosin II functions

Ras guanine nucleotide exchange factor (GEF) Q, a nucleotide exchange factor from Dictyostelium discoideum, is a 143-kD protein containing RasGEF domains and a DEP domain. We show that RasGEF Q can bind to F-actin, has the potential to form complexes with myosin heavy chain kinase (MHCK) A that contain active RasB, and is the predominant exchange factor for RasB. Overexpression of the RasGEF Q GEF domain activates RasB, causes enhanced recruitment of MHCK A to the cortex, and leads to cytokinesis defects in suspension, phenocopying cells expressing constitutively active RasB, and myosin-null mutants. RasGEF Q⁻ mutants have defects in cell sorting and slug migration during later stages of development, in addition to cell polarity defects. Furthermore, RasGEF Q⁻ mutants have increased levels of unphosphorylated myosin II, resulting in myosin II overassembly. Collectively, our results suggest that starvation signals through RasGEF Q to activate RasB, which then regulates processes requiring myosin II.

Choose your own path: specificity in Ras GTPase signaling

The Ras superfamily of small G proteins contributes importantly to numerous cellular and physiological processes (M. F. Olsen and R. Marais, Semin. Immunol., 2000, 12, 63). This family comprises a large class of proteins (more than 150) which all share a common enzymatic function: hydrolysis of the gamma-phosphate of guanosine triphosphate (GTP) to create the products guanosine diphosphate (GDP) and inorganic phosphate (Y. Takai, T. Sasaki and T. Matozaki, Physiol. Rev., 2001, 81, 153). For this reason Ras family proteins, which include the Ras, Rho, Arf/Sara, Ran and Rab subfamilies, are classified as GTPases (G. W. Reuther and C. J. Der, Curr. Opin. Cell Biol., 2000, 12, 157). Guanine nucleotide coupling is a key regulator of enzymatic function; thus, Ras family GTPases participate in signal transduction. Ras signaling depends on binding to effectors. Many of the known effectors can bind to multiple Ras isotypes, often leading to common cellular outcomes, but each Ras isotype also engages specific effector pathways to mediate unique functions. Further, each Ras isotype can propagate multiple signaling pathways, indicating the presence of cellular determinants which allow for promiscuity in Ras-effector interactions while also maintaining specificity. Small distinctions in sequence, structure, and/or cellular regulation contribute to these differences in Ras-effector binding and subsequent cellular effects. A major focus of investigation in the Ras signaling field is identifying the determinants of these individualized functions. This review will attempt to summarize the current state of understanding of this question (with a particular focus on the Ras subfamily) and the approaches being taken to address it, and will discuss prospective areas for future investigation.

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.

GTP-Ras disrupts the intramolecular complex of C1 and RA domains of Nore1

The novel Ras effector mNore1, capable of inducing apoptosis, is a multidomain protein. It comprises a C1 domain homologous to PKC and an RA domain similar to the Ras effectors AF-6 and RalGDS. Here, we determine the affinity of these two domains to the active forms of Ras and Rap1 using isothermal calorimetric titration. The interaction of Ras/Rap1-GTP with the RA domain of mNore1 is weakened significantly by direct binding of the C1 domain to the RA domain. In order to analyze this observation in atomic detail, we solved the C1 solution structure by NMR. By determining chemical shifts and relaxation rates, we can show an intramolecular complex of C1-RA. GTP-Ras titration and binding to RA disrupts this complex and displaces the C1 domain. Once the C1 domain tumbles freely in solution, a lipid binding interface becomes accessible. Furthermore, we provide evidence of phosphatidylinositol 3-phosphate binding of the free C1 domain.

A general method for the unbiased improvement of solution NMR structures by the use of related X-ray data, the AUREMOL-ISIC algorithm

Background

Rapid and accurate three-dimensional structure determination of biological macromolecules is mandatory to keep up with the vast progress made in the identification of primary sequence information. During the last few years the amount of data deposited in the protein data bank has substantially increased providing additional information for novel structure determination projects. The key question is how to combine the available database information with the experimental data of the current project ensuring that only relevant information is used and a correct structural bias is produced. For this purpose a novel fully automated algorithm based on Bayesian reasoning has been developed. It allows the combination of structural information from different sources in a consistent way to obtain high quality structures with a limited set of experimental data. The new ISIC (Intelligent Structural Information Combination) algorithm is part of the larger AUREMOL software package.

Results

Our new approach was successfully tested on the improvement of the solution NMR structures of the Ras-binding domain of Byr2 from Schizosaccharomyces pombe, the Ras-binding domain of RalGDS from human calculated from a limited set of NMR data, and the immunoglobulin binding domain from protein G from Streptococcus by their corresponding X-ray structures. In all test cases clearly improved structures were obtained. The largest danger in using data from other sources is a possible bias towards the added structure. In the worst case instead of a refined target structure the structure from the additional source is essentially reproduced. We could clearly show that the ISIC algorithm treats these difficulties properly.

Conclusion

In summary, we present a novel fully automated method to combine strongly coupled knowledge from different sources. The combination with validation tools such as the calculation of NMR R-factors strengthens the impact of the method considerably since the improvement of the structures can be assessed quantitatively. The ISIC method can be applied to a large number of similar problems where the quality of the obtained three-dimensional structures is limited by the available experimental data like the improvement of large NMR structures calculated from sparse experimental data or the refinement of low resolution X-ray structures. Also structures may be refined using other available structural information such as homology models.

Conformational states of human H-Ras detected by high-field EPR, ENDOR, and 31P NMR spectroscopy

Ras is a central constituent of the intracellular signal transduction that switches between its inactive state with GDP bound and its active state with GTP bound. A number of different X-ray structures are available. Different magnetic resonance techniques were used to characterise the conformational states of the protein and are summarised here. 31P NMR spectroscopy was used as probe for the environment of the phosphate groups of the bound nucleotide. It shows that in liquid solution additional conformational states in the GDP as well as in the GTP forms coexist which are not detected by X-ray crystallography. Some of them can also be detected by solid-state NMR in the micro crystalline state. EPR and ENDOR spectroscopy were used to probe the environment of the divalent metal ion (Mg2+ was replaced by Mn2+) bound to the nucleotide in the protein. Here again different states could be observed. Substitution of normal water by 17O-enriched water allowed the determination of the number of water molecules in the first coordination sphere of the metal ion. In liquid solution, they indicate again the existence of different conformational states. At low temperatures in the frozen state ENDOR spectroscopy suggests that only one state exists for the GDP- and GTP-bound form of Ras, respectively.

A novel mechanism for the modulation of the Ras-effector interaction by small molecules

When proteins require different conformations for their biological function, all these functional states have to coexist simultaneously in solution. However, the corresponding Gibbs free energy differences are usually rather high and thus the conformation with lowest energy predominates in solution whereas the populations of the states with higher energy (excited states) are very small. A stabilization of these excited states can be used as a novel principle to influence the activity of proteins by small molecules. For a proof of this principle, we selected the Ras protein that was shown by (31)P NMR spectroscopy to exist in solution in at least two different conformational states in its GTP form. One of these states shows a drastically reduced affinity to effectors. With Zn(2+)-cyclen we found a small molecule which selectively stabilizes the weak-binding state. It may serve as lead compound for the development of a new type of Ras-inhibitors.

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.

Recognizing and Defining True Ras Binding Domains II: In Silico Prediction Based on Homology Modelling and Energy Calculations

Considering the large number of putative Ras effector proteins, it is highly desirable to develop computational methods to be able to identify true Ras binding molecules. Based on a limited sequence homology among members of the Ras association (RA) and Ras binding (RB) sub-domain families of the ubiquitin super-family, we have built structural homology models of Ras proteins in complex with different RA and RB domains, using the FOLD-X software. A critical step in our approach is to use different templates of Ras complexes, in order to account for the structural variation among the RA and RB domains. The homology models are validated by predicting the effect of mutating hot spot residues in the interface, and residues important for the specificity of interaction with different Ras proteins. The FOLD-X calculated energies of the best-modelled complexes are in good agreement with previously published experimental data and with new data reported here. Based on these results, we can establish energy thresholds above, or below which, we can predict with 96% confidence that a RA/RB domain will or will not interact with Ras. This study shows the importance of in depth structural analysis, high quality force-fields and modelling for correct prediction. Our work opens the possibility of genome-wide prediction for this protein family and for others, where there is enough structural information.

Perturbation of the conformational equilibria in Ras by selective mutations as studied by 31P NMR spectroscopy

Ras regulates a variety of different signal transduction pathways acting as molecular switch. It was shown by liquid and solid-state (31)P NMR spectroscopy that Ras exists in the guanosine-5'-(beta,gamma-imido)triphosphate bound form in at least two conformational states interconverting in millisecond time scale. The relative population between the two conformational states affects drastically the affinity of Ras to its effectors. (31)P NMR spectroscopy shows that the conformational equilibrium can be shifted specifically by point mutations, including mutations with oncogenic potential, thus modifying the effector interactions and their coupling to dynamic properties of the protein.

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.

Ras-effector interactions: after one decade

Ras effectors have convergently developed a common subdomain in their otherwise unrelated protein body for their interaction with Ras. Structural analysis revealed that the mode of interaction is highly similar for all Ras effectors, but is completely different from that of effectors of other subfamilies of small GTPases. Whereas the molecular mechanism of effector activation is still elusive, detailed knowledge about the thermodynamics and dynamics of the interaction with Ras has accumulated.

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.