Transient kinetic studies on the interaction of Ras and the Ras-binding domain of c-Raf-1 reveal rapid equilibration of the complex

Signal transduction at point-blank range: analysis of a spatial coupling mechanism for pathway crosstalk

The plasma membrane provides a physical platform for the orchestration of molecular interactions and biochemical conversions involved in the early stages of receptor-mediated signal transduction in living cells. In that context, we introduce here the concept of spatial coupling, wherein simultaneous recruitment of different enzymes to the same receptor scaffold facilitates crosstalk between different signaling pathways through the local release and capture of activated signaling molecules. To study the spatiotemporal dynamics of this mechanism, we have developed a Brownian dynamics modeling approach and applied it to the receptor-mediated activation of Ras and the cooperative recruitment of phosphoinositide 3-kinase (PI3K) by activated receptors and Ras. Various analyses of the model simulations show that cooperative assembly of multimolecular complexes nucleated by activated receptors is facilitated by the local release and capture of membrane-anchored signaling molecules (such as active Ras) from/by receptor-bound signaling proteins. In the case of Ras/PI3K crosstalk, the model predicts that PI3K is more likely to be recruited by activated receptors bound or recently visited by the enzyme that activates Ras. By this mechanism, receptor-bound PI3K is stabilized through short-range, diffusion-controlled capture of active Ras and Ras/PI3K complexes released from the receptor complex. We contend that this mechanism is a means by which signaling pathways are propagated and spatially coordinated for efficient crosstalk between them.

Effector proteins exert an important influence on the signaling-active state of the small GTPase Cdc42

GTP-binding (G) proteins regulate the flow of information in cellular signaling pathways by alternating between a GTP-bound "active" state and a GDP-bound "inactive" state. Cdc42, a member of the Rho family of Ras-related small G-proteins, plays key roles in the regulation of cell shape, motility, and growth. Here we describe the high resolution x-ray crystal structure for Cdc42 bound to the GTP analog guanylyl beta,gamma-methylene-diphosphonate (GMP-PCP) (i.e. the presumed signaling-active state) and show that it is virtually identical to the structures for the signaling-inactive, GDP-bound form of the protein, contrary to what has been reported for Ras and other G-proteins. Especially surprising was that the GMP-PCP- and GDP-bound forms of Cdc42 did not show detectable differences in their Switch I and Switch II loops. Fluorescence studies using a Cdc42 mutant in which a tryptophan residue was introduced at position 32 of Switch I also showed that there was little difference in the Switch I conformation between the GDP- and GMP-PCP-bound states (i.e. <10%), which again differed from Ras where much larger changes in Trp-32 fluorescence were observed when comparing these two nucleotide-bound states (>30%). However, the binding of an effector protein induced significant changes in the Trp-32 emission specifically from GMP-PCP-bound Cdc42, as well as in the phosphate resonances for GTP bound to this G-protein as indicated in NMR studies. An examination of the available structures for Cdc42 complexed to different effector proteins, versus the x-ray crystal structure for GMP-PCP-bound Cdc42, provides a possible explanation for how effectors can distinguish between the GTP- and GDP-bound forms of this G-protein and ensure that the necessary conformational changes for signal propagation occur.

Photophysical properties of popular fluorescent adenosine nucleotide analogs used in enzyme mechanism probing

Fluorescent nucleotide analogs are widely used in mechanistic studies of nucleotide binding and utilizing proteins. We describe here an overview of the photophysical parameters of the most popular nucleotide analogs that have a fluorescent N-methylanthraniloyl-group attached at various positions of the nucleotide. Steady state absorption and fluorescence spectra of free chromophores depend on the type of modification (ribose, base or phosphate moiety) and the addition of proteins suggests that the labeled nucleotides also vary in sensitivity depending upon their local protein environment. Fluorescence lifetime measurements imply two to three lifetimes for each nucleotide with complex changes in dependence on solvent but more importantly also on the protein. The measured quantum yields quantify the increase in fluorescence for (C8)-MABA-ADP, MANT-ATP and (Pgamma)-MABA-ATP as 153%, 93% and 14% when bound to DnaK, ClpB and Trap1, respectively, compared to free in buffer solution.

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.

B- and C-RAF display essential differences in their binding to Ras: the isotype-specific N terminus of B-RAF facilitates Ras binding

Recruitment of RAF kinases to the plasma membrane was initially proposed to be mediated by Ras proteins via interaction with the RAF Ras binding domain (RBD). Data reporting that RAF kinases possess high affinities for particular membrane lipids support a new model in which Ras-RAF interactions may be spatially restricted to the plane of the membrane. Although the coupling features of Ras binding to the isolated RAF RBD were investigated in great detail, little is known about the interactions of the processed Ras with the functional and full-length RAF kinases. Here we present a quantitative analysis of the binding properties of farnesylated and nonfarnesylated H-Ras to both full-length B- and C-RAF in the presence and absence of lipid environment. Although isolated RBD fragments associate with high affinity to both farnesylated and nonfarnesylated H-Ras, the full-length RAF kinases revealed fundamental differences with respect to Ras binding. In contrast to C-RAF that requires farnesylated H-Ras, cytosolic B-RAF associates effectively and with significantly higher affinity with both farnesylated and nonfarnesylated H-Ras. To investigate the potential farnesyl binding site(s) we prepared several N-terminal fragments of C-RAF and found that in the presence of cysteine-rich domain only the farnesylated form of H-Ras binds with high association rates. The extreme N terminus of B-RAF turned out to be responsible for the facilitation of lipid independent Ras binding to B-RAF, since truncation of this region resulted in a protein that changed its kinase properties and resembles C-RAF. In vivo studies using PC12 and COS7 cells support in vitro results. Co-localization measurements using labeled Ras and RAF documented essential differences between B- and C-RAF with respect to association with Ras. Taken together, these data suggest that the activation of B-RAF, in contrast to C-RAF, may take place both at the plasma membrane and in the cytosolic environment.

Tumour prevention by a single antibody domain targeting the interaction of signal transduction proteins with RAS

Many disease-related processes occur via protein complexes that are considered undruggable with small molecules. An example is RAS, which is frequently mutated in cancer and contributes to initiation and maintenance of the disease by constitutive signal transduction through protein interaction with effector proteins, like PI3K, RAF and RALGDS. Such protein interactions are therefore significant targets for therapy. We describe a single immunoglobulin variable region domain that specifically binds to activated GTP-bound RAS and prevents RAS-dependent tumorigenesis in a mouse model. The crystal structure of the immunoglobulin-RAS complex shows that the variable region competitively binds to the conformationally variant regions of RAS, where its signalling effector molecules interact. This allows the plasma membrane targeted single domain intrabody to inhibit signalling by mutant RAS. This mode of action is a novel advance to directly interfere with oncogenic RAS function in human cancer and shows a universally applicable approach to develop macromolecules to combat cancer. In addition, this method illustrates a general means for interfering with protein interactions that are commonly considered intractable as conventional drug targets.

Massive sequence perturbation of the Raf ras binding domain reveals relationships between sequence conservation, secondary structure propensity, hydrophobic core organization and stability

The contributions of specific residues to the delicate balance between function, stability and folding rates could be determined, in part by [corrected] comparing the sequences of structures having identical folds, but insignificant sequence homology. Recently, we have devised an experimental strategy to thoroughly explore residue substitutions consistent with a specific class of structure. Using this approach, the amino acids tolerated at virtually all residues of the c-Raf/Raf1 ras binding domain (Raf RBD), an exemplar of the common beta-grasp ubiquitin-like topology, were obtained and used to define the sequence determinants of this fold. Herein, we present analyses suggesting that more subtle sequence selection pressure, including propensity for secondary structure, the hydrophobic core organization and charge distribution are imposed on the Raf RBD sequence. Secondly, using the Gibbs free energies (DeltaG(F-U)) obtained for 51 mutants of Raf RBD, we demonstrate a strong correlation between amino acid conservation and the destabilization induced by truncating mutants. In addition, four mutants are shown to significantly stabilize Raf RBD native structure. Two of these mutations, including the well-studied R89L, are known to severely compromise binding affinity for ras. Another stabilized mutant consisted of a deletion of amino acid residues E104-K106. This deletion naturally occurs in the homologues a-Raf and b-Raf and could indicate functional divergence. Finally, the combination of mutations affecting five of 78 residues of Raf RBD results in stabilization of the structure by approximately 12 kJ mol(-1) (DeltaG(F-U) is -22 and -34 kJ mol(-1) for wt and mutant, respectively). The sequence perturbation approach combined with sequence/structure analysis of the ubiquitin-like fold provide a basis for the identification of sequence-specific requirements for function, stability and folding rate of the Raf RBD and structural analogues, highlighting the utility of conservation profiles as predictive tools of structural organization.

Dynamic polymorphism of Ras observed by single molecule FRET is the basis for molecular recognition

Ras regulates signal transduction pathway function by dynamically interacting with various effectors. To understand the basis for Ras function, its conformational dynamics were measured in the absence and presence of effectors using single molecule fluorescence resonance energy transfer (FRET) between probes located on the Switch II region and GTP. The time trajectories of FRET efficiency from GTP-bound Ras showed that this conformation spontaneously varies among multiple states. Among them, a low FRET state was identified as an inactive state. The transition involving the inactive conformational state occurred in the time range of seconds. In contrast, fluctuation occurring most probably between multiple active high FRET conformational states lasted approximately 30 ms but converged to a specific conformational state upon binding to an effector. Thus, Ras conformation spontaneously fluctuates to readily interact with various effectors.

The isolated catalytic hairpin of the Ras-specific guanine nucleotide exchange factor Cdc25Mm retains nucleotide dissociation activity but has impaired nucleotide exchange activity

Cdc25Mm is a mammalian Ras-specific guanine nucleotide exchange factor (GEF). By homology modeling we show that it shares with Sos-GEF the structure of the putative catalytic HI hairpin where the dominant negative T1184E mutation is located. Similarly to Cdc25MmT1184E, the isolated wild-type and mutant hairpins retain the ability to displace Ras-bound nucleotide, originate a stable Ras/GEF complex and downregulate the Ras pathway in vivo. These results indicate that nucleotide re-entry and Ras/GEF dissociation--final steps in the GEF catalytic cycle--require GEF regions different from the HI hairpin. GEF down-sizing could lead to development of novel Ras inhibitors.

The ubiquitin domain superfold: structure-based sequence alignments and characterization of binding epitopes

Ubiquitin-like domains are present, apart from ubiquitin-like proteins themselves, in many multidomain proteins involved in different signal transduction processes. The sequence conservation for all ubiquitin superfold family members is rather poor, even between subfamily members, leading to mistakes in sequence alignments using conventional sequence alignment methods. However, a correct alignment is essential, especially for in silico methods that predict binding partners on the basis of sequence and structure. In this study, using 3D-structural information we have generated and manually corrected sequence alignments for proteins of the five ubiquitin superfold subfamilies. On the basis of this alignment, we suggest domains for which structural information will be useful to allow homology modelling. In addition, we have analysed the energetic and electrostatic properties of ubiquitin-like domains in complex with various functional binding proteins using the protein design algorithm FoldX. On the basis of an in silico alanine-scanning mutagenesis, we provide a detailed binding epitope mapping of the hotspots of the ubiquitin domain fold, involved in the interaction with different domains and proteins. Finally, we provide a consensus fingerprint sequence that identifies all sequences described to belong to the ubiquitin superfold family. It is possible that the method that we describe may be applied to other domain families sharing a similar fold but having low levels of sequence homology.

Current challenges in quantitative modeling of epidermal growth factor signaling

Over the last decade, epidermal growth factor (EGF) signaling has been used repeatedly as a test-bed for pioneering computational systems biology. Recent breakthroughs in our molecular understanding of EGF signaling pose new challenges for mathematical modeling strategies. Three key areas emerge as particularly relevant: the pervasive importance of compartmentalization and endosomal trafficking; the complexity of signalosome complexes; and the regulatory influence of diffusion and spatiality. Each one of them demands a drastic change in current computational approaches. We discuss recent developments in the field that address these emerging aspects in a new generation of more realistic - and potential more useful - models of EGF signaling.

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.

Structural flexibility of small GTPases. Can it explain their functional versatility?

Multiple interactions with many different partners are responsible for the amazing functional versatility of proteins, especially those participating in cellular regulation. The structural properties that could facilitate multiple interactions are examined for small GTPases. The role of cellular constraints, compartmentation and scaffolds on protein-protein interactions is considered.

Relevance of the kinetic equilibrium of forces to the control of the cell cycle by Ras proteins

In higher organisms, the replacement of GDP bound to Ras proteins with GTP, under the participation of an exchange factor, is an important step in the initiation of cell division. Ras-GTP activates kinases and other effectors, which pass signals to the cell nucleus and to the cytoskeleton. The active state of Ras is terminated by hydrolysis of the bound GTP with the assistance of an activating protein (GAP). Knowledge of these regulatory events is based on extensive experimental data, but many aspects of their interpretation are still controversial. It is assumed here that a significant part of the free energy released when two partners associate is stored in a 'kinetic equilibrium of forces' (KEF), and used to facilitate the separation from a third partner. The activation of the Raf kinase is explained primarily in terms of an allosteric effect of Ras-GTP on the phosphate transfer in the catalytic region of the kinase. A mechanism is proposed for the modification of GAP by Ras-GTP, which is believed to be a prerequisite for the well-known crosstalk between the Ras- and Rho-dependent signalling pathways. The cell, by meeting the requirements for KEF, manages to reduce activation barriers, thus significantly accelerating the regulatory events and other complex biological reaction sequences.

Electrostatically optimized Ras-binding Ral guanine dissociation stimulator mutants increase the rate of association by stabilizing the encounter complex

Association of two proteins can be described as a two-step process, with the formation of an encounter complex followed by desolvation and establishment of a tight complex. Here, by using the computer algorithm PARE, we designed a set of mutants of the Ras effector protein Ral guanine nucleotide dissociation stimulator (RalGDS) with optimized electrostatic steering. The fastest binding RalGDS mutant, M26K,D47K,E54K, binds Ras 14-fold faster and 25-fold tighter compared with WT. A linear correlation was found between the calculated and experimental data, with a correlation coefficient of 0.97 and a slope of 0.65 for the 24 mutants produced. The data suggest that increased electrostatic steering specifically stabilizes the encounter complex and transition state. This conclusion is backed up by Phi analysis of the encounter complex and transition state of the RalGDS(M26K,D47K,E54K)/Ras complex, with both values being close to 1. Upon further formation of the final complex, the increased Coulombic interactions are probably counterbalanced by the cost of desolvation of charges, keeping the dissociation rate constant almost unchanged. This mechanism is also reflected by the mutual compensation of enthalpy and entropy changes quantified by isothermal titration calorimetry. The binding constants of the faster binding RalGDS mutants toward Ras are similar to those of Raf, the most prominent Ras effector, suggesting that the design methodology may be used to switch between signal transduction pathways.

Biochemical characterization of the Ran-RanBP1-RanGAP system: are RanBP proteins and the acidic tail of RanGAP required for the Ran-RanGAP GTPase reaction?

RanBP type proteins have been reported to increase the catalytic efficiency of the RanGAP-mediated GTPase reaction on Ran. Since the structure of the Ran-RanBP1-RanGAP complex showed RanBP1 to be located away from the active site, we reinvestigated the reaction using fluorescence spectroscopy under pre-steady-state conditions. We can show that RanBP1 indeed does not influence the rate-limiting step of the reaction, which is the cleavage of GTP and/or the release of product P(i). It does, however, influence the dynamics of the Ran-RanGAP interaction, its most dramatic effect being the 20-fold stimulation of the already very fast association reaction such that it is under diffusion control (4.5 x 10(8) M(-1) s(-1)). Having established a valuable kinetic system for the interaction analysis, we also found, in contrast to previous findings, that the highly conserved acidic C-terminal end of RanGAP is not required for the switch-off reaction. Rather, genetic experiments in Saccharomyces cerevisiae demonstrate a profound effect of the acidic tail on microtubule organization during mitosis. We propose that the acidic tail of RanGAP is required for a process during mitosis.

Probing the sensory rhodopsin II binding domain of its cognate transducer by calorimetry and electrophysiology

Sensory rhodopsin II, a repellent phototaxis receptor from Natronobacterium pharaonis (NpSRII) forms a tight complex with its cognate transducer (NpHtrII). Light excitation of the receptor triggers conformational changes in both proteins, thereby activating the cellular two-component signalling cascade. In membranes, the two proteins form a 2:2 complex, which dissociates to a 1:1 heterodimer in micelles. Complexed to the transducer sensory rhodopsin II is no longer capable of light-driven proton pumping. In order to elucidate the dimerisation and the size of the receptor-binding domain of the transducer, isothermal titration calorimetry and electrophysiological experiments have been carried out. It is shown, that an N-terminal sequence of 114 amino acid residues is sufficient for tight binding (K(d)=240nM; DeltaH=-17.6kJmol(-1)) and for inhibiting the proton transfer. These data and results obtained from selected site-directed mutants indicate a synergistic interplay of transducer transmembrane domain (1-82) and cytoplasmic peptide (83-114) leading to an optimal and specific interaction between receptor and transducer.

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.

Total chemical synthesis of a functional interacting protein pair: the protooncogene H-Ras and the Ras-binding domain of its effector c-Raf1

Generation of biological function by chemical methods is potentially of great importance for the understanding and targeting of physiological processes. Chemical synthesis of proteins offers the ability to alter the properties of target protein molecules in a tailor-made fashion. In the present work it is demonstrated that this methodology can be expanded to the elucidation of protein-protein interactions as exemplified by the complete chemical synthesis of the protooncogene product H-Ras as well as of the Ras-binding domain (RBD) of its effector c-Raf1. The 166-aa polypeptide chain of H-Ras was synthesized by native chemical ligation of three unprotected peptide segments. Similarly, the 81-aa RBD was prepared by ligation of two peptide segments. Both RBD and Ras displayed functional and spectroscopic properties indistinguishable from their recombinant forms as judged by CD spectroscopy and from transient kinetic measurements of the Ras-RBD interaction as well as from nucleotide replacement reactions in Ras. An unnatural amino acid bearing a nitrobenzofurazan side chain was introduced into position 91 of the RBD, providing unique fluorescence properties. The association transient of nitrobenzofurazan labeled with Ras.guanosine 5'-beta,gamma-imidotriphosphate showed a slow phase that had not been detected in earlier work by using other signals.

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.

Observation of slow dynamic exchange processes in Ras protein crystals by 31P solid state NMR spectroscopy

The folding, structure and biological function of many proteins are inherently dynamic properties of the protein molecule. Often, the respective molecular processes are preserved upon protein crystallization, leading, in X-ray diffraction experiments, to a blurring of the electron density map and reducing the resolution of the derived structure. Nuclear magnetic resonance (NMR) is known to be an alternative method to study molecular structure and dynamics. We designed and built a probe for phosphorus solid state NMR that allows for the first time to study static properties as well as dynamic processes in single-crystals of a protein by NMR spectroscopy. The sensitivity achieved is sufficient to detect the NMR signal from individual phosphorus sites in a 0.3mm(3) size single-crystal of GTPase Ras bound to the nucleotide GppNHp, that is, the signal from approximately 10(15) phosphorus nuclei. The NMR spectra obtained are discussed in terms of the conformational variability of the active center of the Ras-nucleotide complex. We conclude that, in the crystal, the protein complex exists in three different conformations. Magic angle spinning (MAS) NMR spectra of a powder sample of Ras-GppNHp show a splitting of one of the phosphate resonances and thus confirm this conclusion. The MAS spectra provide, furthermore, evidence of a slow, temperature-dependent dynamic exchange process in the Ras protein crystal.

Mathematical models of protein kinase signal transduction

We have developed a mathematical theory that describes the regulation of signaling pathways as a function of a limited number of key parameters. Our analysis includes linear kinase-phosphatase cascades, as well as systems containing feedback interactions, crosstalk with other signaling pathways, and/or scaffolding and G proteins. We find that phosphatases have a more pronounced effect than kinases on the rate and duration of signaling, whereas signal amplitude is controlled primarily by kinases. The simplest model pathways allow amplified signaling only at the expense of slow signal propagation. More complex and realistic pathways can combine high amplification and signaling rates with maintenance of a stable off-state. Our models also explain how different agonists can evoke transient or sustained signaling of the same pathway and provide a rationale for signaling pathway design.

Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors

We present a computational model that offers an integrated quantitative, dynamic, and topological representation of intracellular signal networks, based on known components of epidermal growth factor (EGF) receptor signal pathways. The model provides insight into signal-response relationships between the binding of EGF to its receptor at the cell surface and the activation of downstream proteins in the signaling cascade. It shows that EGF-induced responses are remarkably stable over a 100-fold range of ligand concentration and that the critical parameter in determining signal efficacy is the initial velocity of receptor activation. The predictions of the model agree well with experimental analysis of the effect of EGF on two downstream responses, phosphorylation of ERK-1/2 and expression of the target gene, c-fos.

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.

Kinetic studies of protein-protein interactions

The structure of a protein-protein interaction, its affinity and thermodynamic characteristics depict a 'frozen' state of a complex. This picture ignores the kinetic nature of complex formation and dissociation, which are of major biological and biophysical interest. This review highlights recent advances in deciphering the kinetic pathway of protein-protein complexation, the nature of the encounter complex, transition state and intermediate along the reaction, and the effects of mutation, viscosity, pH and salt on association.

New insights into the mechanism of protein-protein association

Association of a protein complex follows a two-step mechanism, with the first step being the formation of an encounter complex that evolves into the final complex. Here, we analyze recent experimental data of the association of TEM1-beta-lactamase with BLIP using theoretical calculations and simulation. We show that the calculated Debye-Hückel energy of interaction for a pair of proteins during association resembles an energy funnel, with the final complex at the minima. All attraction is lost at inter-protein distances of 20 A, or rotation angles of >60 degrees from the orientation of the final complex. For faster-associating protein complexes, the energy funnel deepens and its volume increases. Mutations with the largest impact on association (hotspots for association) have the largest effect on the size and depth of the energy funnel. Analyzing existing evidence, we suggest that the transition state along the association pathway is the formation of the final complex from the encounter complex. Consequently, pairs of proteins forming an encounter complex will tend to dissociate more readily than to evolve into the final complex. Increasing directional diffusion by increasing favorable electrostatic attraction results in a faster forming and slower dissociating encounter complex. The possible applicability of electrostatic calculations for protein-protein docking is discussed.

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

BACKGROUND

After activation, small GTPases such as Ras transfer the incoming signal to effectors by specifically interacting with the binding domain of these proteins. Structural details of the binding domain of different effectors determine which pathway is predominantly activated. Byr2 from fission yeast is a functional homolog of Raf, which is the direct downstream target of Ras in mammalians that initiates a protein kinase cascade. The amino acid sequence of Byr2's Ras binding domain is only weakly related to that of Raf, and Byr2's three-dimensional structure is unknown.

RESULTS

We have solved the 3D structure of the Ras binding domain of Byr2 (Byr2RBD) from Schizosaccharomyces pombe in solution. The structure consists of three alpha helices and a mixed five-stranded beta pleated sheet arranged in the topology betabetaalphabetabetaalphabetaalpha with the first seven canonic secondary structure elements forming a ubiquitin superfold. 15N-(1)H-TROSY-HSQC spectroscopy of the complex of Byr2RBD with Ras*Mg(2+)*GppNHp reveals that the first and second beta strands and the first alpha helix of Byr2 are mainly involved in the protein-protein interaction as observed in other Ras binding domains. Although the putative interaction site of H-Ras from human and Ras1 from S. pombe are identical in sequence, binding to Byr2 leads to small but significant differences in the NMR spectra, indicating a slightly different binding mode.

CONCLUSIONS

The ubiquitin superfold appears to be the general structural motif for Ras binding domains even in cases with vanishing sequence identity. However, details of the 3D structure and the interacting interface are different, thereby determining the specifity of the recognition of Ras and Ras-related proteins.

A method to measure the interaction of Rac/Cdc42 with their binding partners using fluorescence resonance energy transfer between mutants of green fluorescent protein

Enhanced blue fluorescent protein (EBFP) and enhanced green fluorescent protein (EGFP) mutants of GFP in close proximity to one another can act as a fluorescence resonance energy transfer (FRET) pair. Unstructured amino acid linkers of varying length were inserted between EBFP and EGFP, revealing that linkers even as long as 50 amino acids can be accommodated and still allow FRET to occur. This led to the development of a novel biosensor for Rac/Cdc42 binding to their effector proteins based on the insertion of amino acids 75-118 of p21-activated kinase (PAK) between the GFP mutants. We demonstrate that this protein construct allows significant FRET between EBFP and EGFP and retains the ability to bind to Rac in its GTP-bound form with a binding affinity similar to the uncomplexed PAK fragment, and furthermore, on binding to Rac or Cdc42 a marked change in FRET takes place. This forms the basis for a simple, sensitive, and rapid method to measure binding of Rac/Cdc42 to their effector proteins. Since the signal is dependent upon the interaction with active GTP-bound forms it acts as a biosensor for the activation of Rac/Cdc42. It has the potential for use in live cells and for identifying localization of Rac/Cdc42 within subcellular compartments.

Rim1 and rabphilin-3 bind Rab3-GTP by composite determinants partially related through N-terminal alpha -helix motifs

Rim1 is a protein of the presynaptic active zone, the area of the plasma membrane specialized for neurotransmitter exocytosis, and interacts with Rab3, a small GTPase implicated in neurotransmitter vesicle dynamics. Here, we have studied the molecular determinants of Rim1 that are responsible for Rab3 binding, employing surface plasmon resonance and recombinant, bacterially expressed Rab3 and Rim1 proteins. A site that binds GTP- but not GDP-saturated Rab3 was localized to a short alpha-helical sequence near the Rim1 N terminus (amino acids 19-55). Rab3 isoforms A, C, and D were bound with similar affinities (K(d) = 1-2 microm). Low affinity binding of Rab6A-GTP was also observed (K(d) = 16 microm), whereas Rab1B, -5, -7, -8, or -11A did not bind. Adjacent sequences up to amino acid 387, encompassing differentially spliced sequences, the zinc finger module, and the SGAWFF motif of Rim1, did not significantly contribute to the strength or the specificity of Rab3 binding, whereas a point mutation within the helix (R33G) abolished binding. This Rab3 binding site of Rim1 is reminiscent of the N-terminal alpha-helix that is part of the Rab3-binding region of rabphilin-3, and indeed we observed low affinity, specific binding of Rab3A (K(d) on the order of magnitude of 10-100 microm) to this region of rabphilin-3 alone (amino acids 40-88), whereas additional sequences up to amino acid 178 are needed for high affinity Rab3A binding to rabphilin-3 (K(d) = 10-20 nm). In contrast, an N-terminal alpha-helix motif in aczonin, with sequence similarity to the Rab3-binding site of Rim1, did not bind Rab3A, -C, or -D or several other Rab proteins. These results were qualitatively confirmed in pull-down experiments with native, prenylated Rab3 from brain lysate in Triton X-100. Munc13 bound to the zinc finger domain of Rim1 but not to the rabphilin-3 or aczonin zinc fingers. Pull-down experiments from brain lysate in the presence of cholate as detergent detected binding to downstream Rim1 sequences, between amino acids 56 and 387, of syntaxin and of Rab3. The latter, however, was inhibited rather than stimulated by GTP.

Dynamic properties of the Ras switch I region and its importance for binding to effectors

We have investigated the dynamic properties of the switch I region of the GTP-binding protein Ras by using mutants of Thr-35, an invariant residue necessary for the switch function. Here we show that these mutants, previously used as partial loss-of-function mutations in cell-based assays, have a reduced affinity to Ras effector proteins without Thr-35 being involved in any interaction. The structure of Ras(T35S)(.)GppNHp was determined by x-ray crystallography. Whereas the overall structure is very similar to wildtype, residues from switch I are completely invisible, indicating that the effector loop region is highly mobile. (31)P-NMR data had indicated an equilibrium between two rapidly interconverting conformations, one of which (state 2) corresponds to the structure found in the complex with the effectors. (31)P-NMR spectra of Ras mutants (T35S) and (T35A) in the GppNHp form show that the equilibrium is shifted such that they occur predominantly in the nonbinding conformation (state 1). On addition of Ras effectors, Ras(T35S) but not Ras(T35A) shift to positions corresponding to the binding conformation. The structural data were correlated with kinetic experiments that show two-step binding reaction of wild-type and (T35S)Ras with effectors requires the existence of a rate-limiting isomerization step, which is not observed with T35A. The results indicate that minor changes in the switch region, such as removing the side chain methyl group of Thr-35, drastically affect dynamic behavior and, in turn, interaction with effectors. The dynamics of the switch I region appear to be responsible for the conservation of this threonine residue in GTP-binding proteins.

A sensitive fluorescence monitor for the detection of activated Ras: total chemical synthesis of site-specifically labeled Ras binding domain of c-Raf1 immobilized on a surface

BACKGROUND

The Ras.GDP-Ras.GTP cycle plays a central role in eukaryotic signaling cascades. Mutations in Ras which stabilize activated Ras.GTP lead to a continuous stimulation of downstream effectors and ultimately to cell proliferation. Ras mutants which increase the steady-state concentration of Ras.GTP are involved in about 30% of all human cancers. It is therefore of great interest to develop a biosensor which is sensitive to Ras.GTP but not to Ras.GDP.

RESULTS

The Ras binding domain (RBD) of c-Raf1 was synthesized from two unprotected peptide segments by native chemical ligation. Two fluorescent amino acids with structures based on the nitrobenz-2-oxa-1,3-diazole and coumaryl chromophores were incorporated at a site which is close to the RBD/Ras.GTP binding surface. Additionally, a C-terminal tag consisting of His6 was introduced. The Kd values for binding of the site-specifically modified proteins to Ras.GTP are comparable to that of wild-type RBD. Immobilization of C-terminal His6 tag-modified fluorescent RBD onto Ni-NTA-coated surfaces allowed the detection of Ras.GTP in the 100 nM range. Likewise, Ras.GTP/Q61L (an oncogenic mutant of Ras with very low intrinsic GTP hydrolysis activity) can also be detected in this assay system. Ras.GDP does not bind to the immobilized RBD, thus allowing discrimination between inactive and activated Ras.

CONCLUSIONS

The site-specific incorporation of a fluorescent group at a strategic position in a Ras effector protein allows the detection of activated Ras with high sensitivity. This example illustrates the fact that the chemical synthesis of proteins or protein domains makes it possible to incorporate any kind of natural or unnatural amino acid at the position of choice, thereby enabling the facile preparation of specific biosensors, enhanced detection systems for drug screening, or the synthesis of activated proteins, e.g. phosphorylated proteins involved in signaling pathways, as defined molecular species.

Design of inhibitors of Ras--Raf interaction using a computational combinatorial algorithm

Drugs that inhibit important protein-protein interactions are hard to find either by screening or rational design, at least so far. Most drugs on the market that target proteins today are therefore aimed at well-defined binding pockets in proteins. While computer-aided design is widely used to facilitate the drug discovery process for binding pockets, its application to the design of inhibitors that target the protein surface initially seems to be limited because of the increased complexity of the task. Previously, we had started to develop a computational combinatorial design approach based on the well-known 'multiple copy simultaneous search' (MCSS) procedure to tackle this problem. In order to identify sequence patterns of potential inhibitor peptides, a three-step procedure is employed: first, using MCSS, the locations of specific functional groups on the protein surface are identified; second, after constructing the peptide main chain based on the location of favorite locations of N-methylacetamide groups, functional groups corresponding to amino acid side chains are selected and connected to the main chain C(alpha) atoms; finally, the peptides generated in the second step are aligned and probabilities of amino acids at each position are calculated from the alignment scheme. Sequence patterns of potential inhibitors are determined based on the propensities of amino acids at each C(alpha) position. Here we report the optimization of inhibitor peptides using the sequence patterns determined by our method. Several short peptides derived from our prediction inhibit the Ras--Raf association in vitro in ELISA competition assays, radioassays and biosensor-based assays, demonstrating the feasibility of our approach. Consequently, our method provides an important step towards the development of novel anti-Ras agents and the structure-based design of inhibitors of protein--protein interactions.

Crystal structure and functional analysis of Ras binding to its effector phosphoinositide 3-kinase gamma

Ras activation of phosphoinositide 3-kinase (PI3K) is important for survival of transformed cells. We find that PI3Kgamma is strongly and directly activated by H-Ras G12V in vivo or by GTPgammaS-loaded H-Ras in vitro. We have determined a crystal structure of a PI3Kgamma/Ras.GMPPNP complex. A critical loop in the Ras binding domain positions Ras so that it uses its switch I and switch II regions to bind PI3Kgamma. Mutagenesis shows that interactions with both regions are essential for binding PI3Kgamma. Ras also forms a direct contact with the PI3Kgamma catalytic domain. These unique Ras/PI3Kgamma interactions are likely to be shared by PI3Kalpha. The complex with Ras shows a change in the PI3K conformation that may represent an allosteric component of Ras activation.

Biochemical characterization of Rab3-GTPase-activating protein reveals a mechanism similar to that of Ras-GAP

Small G proteins of the Rab family are regulators of intracellular vesicle traffic. Their intrinsic rate of GTP hydrolysis is very low but is enhanced by specific GTPase-activating proteins (GAPs) that switch G proteins to their inactive form. We have characterized the activity of recombinant Rab3-GAP on Rab3A in solution. The K(m) and K(d) values (75 microm) indicate a low affinity of Rab3-GAP for its substrate. The affinity is higher for the transition state analog Rab3A:GDP:AlF(x) (15 microm). The k(cat) (1 s(-)(1)) is within the range of values reported for other GAPs. A mutation in the switch I region of Rab3A disrupted the interaction with Rab3-GAP. Furthermore, Rabphilin, a putative target of Rab3, inhibited the activity of Rab3-GAP on Rab3. Therefore, the Rab3-GAP-binding site involves the switch I region of Rab3 and overlaps with the Rabphilin-binding domain. Substitution of a single arginine residue (Arg-728) of Rab3-GAP disrupted its catalytic activity but not its interaction with Rab3A. We propose that Rab3-GAP, like Ras- and Rho-GAPs, stabilizes the transition state of Rab3 and provides a critical arginine residue to accelerate the GTPase reaction.

Point mutants of c-raf-1 RBD with elevated binding to v-Ha-Ras

A mutational analysis of the Ras-binding domain (RBD) of c-Raf-1 identified three amino acid positions (Asn(64), Ala(85), and Val(88)) where amino acid substitution with basic residues increases the binding of RBD to recombinant v-Ha-Ras. The greatest increase in binding (6-9-fold) was observed with the A85K-RBD mutant. The elevated binding for the A85K-RBD and V88R-RBD mutants was also detected with Ras expressed in cultured mammalian cells, namely NIH-3T3 and BAF cells. None of the wild type residues in RBD positions Asn(64), Ala(85), and Val(88) have been previously implicated in the interaction with Ras (Block, C., Janknecht, R., Herrmann, C., Nassar, N., and Wittinghofer, A. (1996) Nat. Struct. Biol. 3, 244-251; Nassar, N., Horn, G., Herrmann, C., Scherer, A., McCormick, F., and Wittinghofer, A. (1995) Nature 375, 554-560). The discovery of elevated binding among the mutants in these positions implies that additional RBD residues can be used to generate the Ras. RBD complex. These findings are of particular significance in the design of Ras antagonists based on the RBD prototype. The A85K-RBD mutant can be used to develop an assay for measuring the level of activated Ras in cultured cells; Sepharose-linked A85K-RBD.GST fusion protein served as an activation-specific probe to precipitate Ras.GTP but not Ras.GDP from epidermal growth factor-stimulated cells. A85K-RBD precipitates up to 5-fold more Ras.GTP from mammalian cells than wild type RBD.

Kinetic analysis of peptide binding to the TAP transport complex: evidence for structural rearrangements induced by substrate binding

The transporter associated with antigen processing (TAP) plays a key role in the class I major histocompatibility complex (MHC) mediated immune surveillance. It translocates peptides generated by the proteasome complex into the endoplasmic reticulum (ER) for loading onto MHC class I molecules. At the cell surface these MHC complexes are monitored for their antigenic cargo by cytotoxic T-lymphocytes. Peptide binding to TAP is the essential step for peptide selection and for subsequent ATP-dependent translocation into the ER lumen. To examine the pathway of substrate recognition by TAP, we employed peptide epitopes, which were labeled with an environmentally sensitive fluorophore. Upon binding to TAP, a drastic fluorescence quenching of the fluorescent substrate was detected. This allowed us to analyze TAP function in real-time by using a homogeneous assay. Formation of the peptide-TAP complex is composed of a fast association step followed by a slow isomerization of the transport complex. Proton donor groups moving in proximity to the fluorescence label cause fluorescence quenching. Taken together, this peptide-induced structural reorganization may reflect the crosstalk of structural information between the peptide binding site and both nucleotide-binding domains within the TAP complex.

The pre-hydrolysis state of p21(ras) in complex with GTP: new insights into the role of water molecules in the GTP hydrolysis reaction of ras-like proteins

BACKGROUND

In numerous biological events the hydrolysis of guanine triphosphate (GTP) is a trigger to switch from the active to the inactive protein form. In spite of the availability of several high-resolution crystal structures, the details of the mechanism of nucleotide hydrolysis by GTPases are still unclear. This is partly because the structures of the proteins in their active states had to be determined in the presence of non-hydrolyzable GTP analogues (e.g. GppNHp). Knowledge of the structure of the true Michaelis complex might provide additional insights into the intrinsic protein hydrolysis mechanism of GTP and related nucleotides.

RESULTS

The structure of the complex formed between p21(ras) and GTP has been determined by X-ray diffraction at 1.6 A using a combination of photolysis of an inactive GTP precursor (caged GTP) and rapid freezing (100K). The structure of this complex differs from that of p21(ras)-GppNHp (determined at 277K) with respect to the degree of order and conformation of the catalytic loop (loop 4 of the switch II region) and the positioning of water molecules around the gamma-phosphate group. The changes in the arrangement of water molecules were induced by the cryo-temperature technique.

CONCLUSIONS

The results shed light on the function of Gln61 in the intrinsic GTP hydrolysis reaction. Furthermore, the possibility of a proton shuffling mechanism between two attacking water molecules and an oxygen of the gamma-phosphate group can be proposed for the basal GTPase mechanism, but arguments are presented that render this protonation mechanism unlikely for the GTPase activating protein (GAP)-activated GTPase.