Thermodynamics of Ras/effector and Cdc42/effector interactions probed by isothermal titration calorimetry

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.

The P34G mutation reduces the transforming activity of K-Ras and N-Ras in NIH 3T3 cells but not of H-Ras

Ras proteins (H-, N-, and K-Ras) operate as molecular switches in signal transduction cascades controlling cell proliferation, differentiation, or apoptosis. The interaction of Ras with its effectors is mediated by the effector-binding loop, but different data about Ras location to plasma membrane subdomains and new roles for some docking/scaffold proteins point to signaling specificities of the different Ras proteins. To investigate the molecular mechanisms for these specificities, we compared an effector loop mutation (P34G) of three Ras isoforms (H-, N-, and K-Ras4B) for their biological and biochemical properties. Although this mutation diminished the capacity of Ras proteins to activate the Raf/ERK and the phosphatidylinositol 3-kinase/AKT pathways, the H-Ras V12G34 mutant retained the ability to cause morphological transformation of NIH 3T3 fibroblasts, whereas both the N-Ras V12G34 and the K-Ras4B V12G34 mutants were defective in this biological activity. On the other hand, although both the N-Ras V12G34 and the K-Ras4B V12G34 mutants failed to promote activation of the Ral-GDS/Ral A/PLD and the Ras/Rac pathways, the H-Ras V12G34 mutant retained the ability to activate these signaling pathways. Interestingly, the P34G mutation reduced specifically the N-Ras and K-Ras4B in vitro binding affinity to Ral-GDS, but not in the case of H-Ras. Thus, independently of Ras location to membrane subdomains, there are marked differences among Ras proteins in the sensitivity to an identical mutation (P34G) affecting the highly conserved effector-binding loop.

Crystal structure of Rab9 complexed to GDP reveals a dimer with an active conformation of switch II

The small GTPase Rab9 is an essential regulator of vesicular transport from the late endosome to the trans-Golgi network, as monitored by the redirection of the mannose-6-phosphate receptors. The crystal structure of Rab9 complexed to GDP, Mg(2+), and Sr(2+) reveals a unique dimer formed by an intermolecular beta-sheet that buries the switch I regions. Surface area and shape complementarity calculations suggest that Rab9 dimers can form an inactive, membrane-bound pool of Rab9 . GDP that is independent of GDI. Mg(2+)-bound Rab9 represents an inactive state, but Sr(2+)-bound Rab9 . GDP displays activated switch region conformations, mimicking those of the GTP state. A hydrophobic tetrad is formed resembling an effector-discriminating epitope found only in GTP-bound Rab proteins.

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.

Comparative functional analysis of the Rac GTPases

Small GTPases of the Rho family including Rac, Rho and Cdc42 regulate different cellular processes like reorganization of the actin cytoskeleton by acting as molecular switches. The three distinct mammalian Rac proteins share very high sequence identity but how their specificity is achieved is hitherto unknown. Here we show that Rac1 and Rac3 are very closely related concerning their biochemical properties, such as effector interaction, nucleotide binding and hydrolysis. In contrast, Rac2 displays a slower nucleotide association and is more efficiently activated by the Rac-GEF Tiam1. Modeling and normal mode analysis support the idea that altered dynamics of Rac2 at the switch I region may be responsible for different biochemical properties. These results provide insight into the individual functionalities of the Rac isoforms.

Differential interaction and activation of Rho family GTPases by p210bcr-abl and p190bcr-abl

The p210(bcr-abl) and p190(bcr-abl) fusion proteins, respectively responsible for chronic myelogenous leukemia and acute lymphoblastic leukemia, present deregulated tyrosine kinase activity and abnormal localization. The Dbl homology domain of Bcr, activating Rho GTPases, is present in p210(bcr-abl), but absent in p190(bcr-abl). We investigated the interaction of Bcr-Abl chimeras and Rho proteins by coimmunoprecipitation, pull-down experiments and GEF activity measurement. RhoA, Rac1 and Cdc42 interact in vivo with p210(bcr-abl) only. Moreover, the three types of GTPases are activated in vitro and in vivo by p210(bcr-abl). Nevertheless, Rac1 and Cdc42, but not RhoA, are activated by p190(bcr-abl) in vitro and in vivo. Part of this GEF activity of p190(bcr-abl) is probably attributable to p95(vav), which is complexed with both p190(bcr-abl) and p210(bcr-abl) in an activated form. p160(bcr), also in complex with Bcr-Abl, presents no GEF activity in p190(bcr-abl)-expressing cells. These results suggest that differential activation of Rho proteins should play a major role in Bcr-Abl-induced leukemogenesis.

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.

Applications of calorimetric methods to drug discovery and the study of protein interactions

Recent studies report the application of isothermal titration calorimetry and differential scanning calorimetry to the study of protein-ligand interactions, allosteric cooperativity and aspects of protein folding. New methods of data analysis compare alternative methods for determining ligand binding enthalpy and analyze potential sources of error in the experimental measurement of other thermodynamic parameters. Several reports examine issues concerning drug design and the correlation of thermodynamic and X-ray structural data. New instruments allow volumetric effects in biochemical systems to be evaluated calorimetrically and to substantially expand the throughput of differential scanning calorimetry measurements in drug discovery and other high-throughput applications.

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.

Converging free energy estimates: MM-PB(GB)SA studies on the protein-protein complex Ras-Raf

Estimating protein-protein interaction energies is a very challenging task for current simulation protocols. Here, absolute binding free energies are reported for the complex H-Ras/C-Raf1 using the MM-PB(GB)SA approach, testing the internal consistency and model dependence of the results. Averaging gas-phase energies (MM), solvation free energies as determined by Generalized Born models (GB/SA), and entropic contributions calculated by normal mode analysis for snapshots obtained from 10 ns explicit-solvent molecular dynamics in general results in an overestimation of the binding affinity when a solvent-accessible surface area-dependent model is used to estimate the nonpolar solvation contribution. Applying the sum of a cavity solvation free energy and explicitly modeled solute-solvent van der Waals interaction energies instead provides less negative estimates for the nonpolar solvation contribution. When the polar contribution to the solvation free energy is determined by solving the Poisson-Boltzmann equation (PB) instead, the calculated binding affinity strongly depends on the atomic radii set chosen. For three GB models investigated, different absolute deviations from PB energies were found for the unbound proteins and the complex. As an alternative to normal-mode calculations, quasiharmonic analyses have been performed to estimate entropic contributions due to changes of solute flexibility upon binding. However, such entropy estimates do not converge after 10 ns of simulation time, indicating that sampling issues may limit the applicability of this approach. Finally, binding free energies estimated from snapshots of the unbound proteins extracted from the complex trajectory result in an underestimate of binding affinity. This points to the need to exercise caution in applying the computationally cheaper "one-trajectory-alternative" to systems where there may be significant changes in flexibility and structure due to binding. The best estimate for the binding free energy of Ras-Raf obtained in this study of -8.3 kcal mol(-1) is in good agreement with the experimental result of -9.6 kcal mol(-1), however, further probing the transferability of the applied protocol that led to this result is necessary.