Electrostatic contributions to protein-protein binding affinities: application to Rap/Raf interaction.

The challenge of evaluating absolute binding free energies of protein-protein complexes is addressed using the scaled Protein Dipoles Langevin Dipoles (PDLD/S) model in combination with the Linear Response Approximation (LRA). This is done by taking the complex between Rap1A (Rap) and the p21ras binding domain of c-Raf (Raf-RBD) (Nassar et al., Nature 375:554-560, 1995) as a model system. Several formulations and different thermodynamic cycles are explored taking advantage of the LRA method and considering the protein reorganization during complex formation. The performance of different approximations is examined by comparing the calculated and observed absolute binding energies for the native complex and some of its mutants. The evaluation of the contributions of individual residues to the binding free energy, which is referred to here as group contributions is also examined. Special attention is paid to the role of the "dielectric constant," epsilon(in) which is in fact a scaling factor that represents the contributions that are treated implicitly. It is found that explicit consideration of protein relaxation is crucial for obtaining reasonable results with small values of epsilon(in), but it is also found that such a treatment of protein-protein interactions is very challenging and does not always give stable results. This indicates that more advanced explicit calculations should be based on experimentally determined structures of both the complex and the isolated proteins. Nevertheless, it is demonstrated that the qualitative trend of the effect of mutations can be reproduced by considering the effect of protein reorganization implicitly, using epsilon(in) approximately 25 for ionized residues and epsilon(in) approximately 4 for polar residues. Thus, it is concluded that an explicit treatment of solvent relaxation (which is common to current continuum models) does not provide sufficient compensation for turning off the charges of ionized residues on the interaction surface of the Raf-RBD/Rap complex. Representing the missing contribution by large epsilon(in) can, of course, reproduce the observed effect of ionized residues, but now the contribution of uncharged residues will be largely underestimated. Regardless of these conceptual problems, it is established that a very simple nonrelaxed approach, where the relaxation of both the protein and the solvent are considered implicitly, can provide an effective qualitative way for evaluating group contributions, using large and small values for epsilon(in) of ionized and neutral residues, respectively. As much as the actual system studied is concerned we find that more residues than generally assumed play a role in Raf-RBD/Rap interaction. This includes residues that are not located at the protein-protein interaction surface. These residues contribute to the binding energy through direct charge-charge interaction without leading to drastic structural changes. The overall contribution of the surface residues is quite significant since Raf and Rap are positively and negatively charged, respectively, and their charges are distributed along the interaction site between the two proteins.

The role of the metal ion in the p21ras catalysed GTP-hydrolysis: Mn2+ versus Mg2+.

GTP and ATP hydrolysing proteins have an absolute requirement for a divalent cation, which is usually Mg2+, as a cofactor in the enzymatic reaction. Other phosphoryl transfer enzymes employ more than one divalent ion for the enzymatic reaction. It is shown here for p21ras, a well studied example of GTP hydrolysing proteins, that the GTP-hydrolysis rate is significantly faster if Mg2+ is replaced by Mn2+, both in the presence or absence of its GTPase-activating protein Ras-GAP. This effect is not due to a different stoichiometry of metal ion binding, since one metal ion is sufficient for full enzymatic activity. To determine the role of the metal ion, the crystal structure of p21(G12P). GppCp complexed with Mn2+ was determined and shown to be very similar to the corresponding p21(G12P). GppCp.Mg2+ structure. Especially the coordination sphere around the metal ions is very similar, and no second metal ion binding site could be detected, consistent with the assumption that one metal ion is sufficient for GTP hydrolysis. In order to explain the biochemical differences, we analysed the GTPase reaction mechanism with a linear free energy relationships approach. The result suggests that the reaction mechanism is not changed with Mn2+ but that the transition metal ion Mn2+ shifts the pKa of the gamma-phosphate by almost half a unit and increases the reaction rate due to an increase in the basicity of GTP acting as the general base. This suggests that the intrinsic GTPase reaction could be an attractive target for anti-cancer drug design. By using Rap1A and Ran, we show that the acceleration of the GTPase by Mn2+ appears to be a general phenomenon of GTP-binding proteins.

Linear free energy relationships in the intrinsic and GTPase activating protein-stimulated guanosine 5'-triphosphate hydrolysis of p21ras.

Controlling the hydrolysis rate of GTP bound to guanine nucleotide binding proteins is crucial for the right timing of many biological processes. Theoretical, structural, and functional studies have demonstrated that in p21ras the substrate of the reaction, GTP itself, plays a central role by acting as the base catalyst. This substrate-assisted reaction mechanism was analyzed with the help of linear free energy relationships (LFERs). Here we present experimental data that further support the proposed mechanism. We extend the LFER analysis to a wide range of oncogenic as well as nontransforming Ras mutants. It is illustrated that almost all Ras variants follow the observed LFER and thus also the same reaction path. Further, the reduced GTPase reaction rate that characterizes the oncogenic effect of many of the p21 mutants found in human tumors seems to be a consequence of a slightly reduced pKa of the gamma-phosphate group of bound GTP. Factors causing a pKa deviation of just 0.5 unit are enough to slow the intrinsic GTPase reaction rate significantly, and the system may exhibit as a consequence of this an oncogenic potential. Interestingly, we also found oncogenic mutations that do not follow the regular LFER. This suggests that the oncogenic effect of distinct Ras mutants has a different physical origin. The results presented might aid in the design of drugs aimed at reactivating the GTPase reaction of many oncogenic p21ras mutants. We also analyzed the stimulated GTPase reaction of p21ras by the GTPase activating protein (GAP) and the GTPase reaction of Rap1A, a Ras-related GTP binding protein, with similar approaches. The corresponding results indicate that the GAP-stimulated GTPase as well as the Rap1A-catalyzed reaction seem to follow the same substrate-assisted reaction mechanism. However, the correlation coefficient for the GAP-catalyzed reaction is different from the corresponding coefficient for the intrinsic reaction. While the intrinsic reaction exhibits a Brønsted slope of beta = 2.1, the corresponding value for the GAP-activated reaction is beta = 4.9.

Mechanistic analysis of the observed linear free energy relationships in p21ras and related systems.

Previous studies of the GTPase reaction catalyzed by p21ras have indicated that the logarithm of the observed reaction rate and the pKa of the bound GTP are correlated by the Brønsted relationship log(kcat) = beta pKa + A. While most of the Ras mutants display a Brønsted slope beta of 2.1, a small set of oncogenic mutants exhibit a beta of > > 1. On the other hand, it was found that the corresponding Brønsted slope for the GTPase reaction of p21ras in the presence of GTPase Activating Protein (GAP) is about beta = 4.9. The present work explores the basis for such linear free energy relationships (LFERs) in general and applies these concepts to p21ras and related systems. It is demonstrated that the optimal way to analyze LFER is by using Marcus type parabolas that represent the reactant, intermediate, and product state of the reaction in a relevant energy diagram. The observed LFER is used to analyze the actual free energy surface and reaction path of the intrinsic GTPase reaction in p21ras. From this, a model reaction profile can be constructed that explains how a LFER can arise and also how the different observed Brønsted coefficients can be rationalized. This analysis is augmented by solvent isotope effect studies. It is pointed out that the overall activation barrier reflects the energy of the proton transfer (PT) step, although this step does not include the actual transition state of the hydrolysis reaction. The proposed GTP as a base mechanism is compared to a recently proposed reaction scheme where Gln61 serves as a proton shuttle in a concerted mechanism. It is shown by unique energy considerations that the concerted mechanism is unlikely. Other alternative mechanisms are also considered, and their consistency with the observed LFER and other factors is discussed. Finally, we analyze the observed LFER for the GTPase reaction of p21ras in the presence of GAP and discuss its relevance for the mechanism of GAP activation.

Conformational transitions in p21ras and in its complexes with the effector protein Raf-RBD and the GTPase activating protein GAP.

31P NMR revealed that the complex of p21ras with the GTP analog GppNHp.Mg2+ exists in two conformational states, states 1 and 2. In wild-type p21ras the equilibrium constant K1(12) between the two states is 1.09. The population of these states is different for various mutants but independent of temperature. The activation enthalpy delta H ++ and activation entropy delta S ++ for the conformational transitions were determined by full-exchange matrix analysis for wild-type p21ras and p21ras(S65P). For the wild-type protein one obtains delta H ++ = 89 +/- 2 kJ mol-1 and delta S ++ = 102 +/- 20 J mol-1 K-1 and for the mutant protein delta H ++ = 93 +/- 7 kJ mol-1 and delta S ++ = 138 +/- 30 J mol-1 K-1. The study of various p21ras mutants suggests that the two states correspond to different conformations of loop L2, with Tyr-32 in two different positions relative to the bound nucleotide. High-field EPR at 95 GHz suggest that the observed conformational transition does not directly influence the coordination sphere of the protein-bound metal ion. The influence of this transition on loop L4 was studied by 1H NMR with mutants E62H and E63H. There was no indication that L4 takes part in the transition described in L2, although a reversible conformational change could be induced by decreasing the pH value. The exchange between the two states is slow on the NMR time scale (< 10 s-1): at approximately pH 5 the population of the two states is equal. The interaction of p21ras-triphosphate complexes with the Ras-binding domain (RBD) of the effector protein c-Raf-1, Raf-RBD, and with the GTPase activating protein GAP was studied by 31P NMR spectroscopy. In complex with Raf-RBD the second conformation of p21ras (state 2) is stabilized. In this conformation Tyr-32 is located in close proximity to the phosphate groups of the nucleotide, and the beta-phosphate resonance is shifted upfield by 0.7 ppm. Spectra obtained in the presence of GAP suggest that in the ground state GAP does not interact directly with the nucleotide bound to p21ras and does not induce larger conformational changes in the neighborhood of the nucleotide. The experimental data are consistent with a picture where GAP accelerates the exchange process between the two states and simultaneously increases the population of state 1 at higher temperature.

Electrostatic control of GTP and GDP binding in the oncoprotein p21ras.

BACKGROUND: p21ras is one of the GTP-binding proteins that act as intercellular molecular switches. The GTP-bound form of p21ras sends a growth-promoting signal that is terminated once the protein is cycled back into its GDP-bound form. The interaction of guanine-nucleotide-exchange factors (GEFs) with p21ras leads to activation of the protein by promoting GDP --> GTP exchange. Oncogenic mutations of p21ras trap the protein in its biological active GTP-bound form. Other mutations interfere with the activity of GEF. Thus, it is important to explore the structural basis for the action of different mutations. RESULTS: The crystal structures of p21ras are correlated with the binding affinities of GTP and GDP by calculating the relevant electrostatic energies. It is demonstrated that such calculations can provide a road map to the location of 'hot' residues whose mutations are likely to change functional properties of the protein. Furthermore, calculations of the effect of specific mutations on GTP and GDP binding are consistent with those observed. This helps to analyze and locate functionally important parts of the protein. CONCLUSIONS: Our calculations indicate that the protein main chain provides a major contribution to the binding energies of nucleotides and probably plays a key role in relaying the effect of GEF action. Analysis of p21ras mutations in residues that are important for the proper function of GEFs suggests that the region comprising residues 62-67 in p21ras is the major GEF-binding site. This analysis and our computer simulations indicate that the effect of GEF is probably propagated to the P-loop (residues 10-17) through interaction between Gly60 and Gly12. This then reduces the interaction between the main-chain dipoles of the P-loop and the nucleotide. Finally, the results also suggest a possible relationship between the GTP --> GDP structural transition and the catalytic effect of the GTPase-activating protein.

Substrate-assisted catalysis as a mechanism for GTP hydrolysis of p21ras and other GTP-binding proteins.

Despite many advances in understanding the structure and function of GTP-binding proteins the mechanism by which these molecules switch from the GTP-bound on-state to the GDP-bound off-state is still poorly understood. Theoretical studies suggest that the activation of the nucleophilic water which hydrolyzes GTP needs a general base. Such a base could not be located in any of the many GTP-binding proteins. Here we present a unique type of linear free energy relationships that not only supports a mechanism for p21ras in which the substrate GTP itself acts as the catalytic base driving the GTPase reaction but can also help to explain why certain mutants of p21ras are oncogenic and others are not.

Why have mutagenesis studies not located the general base in ras p21.

Ras p21 plays a major role in the control of cell growth, and oncogenic mutations of this protein have been found in human cancers. Unfortunately, the detailed mode of action of Ras p21 is still unclear, in spite of the great interest in this protein and the availability of its X-ray crystal structure. In particular, mutagenesis studies of different active site residues could not identify the general base for GTP hydrolysis. Here we tackle this question using a computer simulation approach with clear and reliable energy considerations and conclude that the most likely general base is the bound GTP itself. Obviously, the identification of such a general base cannot be easily accomplished by mutagenesis experiments.

GTP-binding proteins. Structures, interactions and relationships.

Recently available crystal structures show that some, though not all, GTP-binding proteins have a common 'G-domain' topology, variations on which confer distinct functional properties.

On the mechanism of guanosine triphosphate hydrolysis in ras p21 proteins.

The residue Gln61 is assumed to play a major role in the mechanism of ras p21, and mutations of this residue are often found in human tumors. Such mutations lead to a major reduction in the rate of GTP hydrolysis by the complex of ras p21 and the GTPase activating protein (GAP) and lock the protein in a growth-promoting state. This work examines the role of Gln61 in ras p21 by using computer simulation approaches to correlate the structure and energetics of this system. Free energy perturbation calculations and simpler electrostatic considerations demonstrate that Gln61 is unlikely to serve as the general base in the intrinsic GAP-independent reaction of p21. Glutamine is already a very weak base in water, and surprisingly the GlnH+ OH-reaction intermediate is even less stable in the protein active site than in the corresponding reaction in water. The electrostatic field of Glu63, which could in principle stabilize the protonated Gln61, is found to be largely shielded by the surrounding solvent. However, it is still possible that Gln61 is a general base in the GAP/ras p21 complex since this system could enhance the electrostatic effect of Glu63. It is also possible that the gamma-phosphate acts as general base and that Gln61 accelerates the reaction by stabilizing the OH- nucleophile. If such a mechanism is operative, then GAP may enhance the effect of Gln61 by preorienting its hydrogen bonds in the transition-state configuration.