Search PubMed⌕ Search

Biomedical subjects

A Warshel

Publications and source records attributed to A Warshel.

At least 19 recordsLinked to original sources
PubMed2000

How important are entropic contributions to enzyme catalysis?

The idea that enzymes accelerate their reactions by entropic effects has played a major role in many prominent proposals about the origin of enzyme catalysis. This idea implies that the binding to an enzyme active site freezes the motion of the reacting fragments and eliminates their entropic contributions, (delta S(cat)(double dagger))', to the activation energy. It is also implied that the binding entropy is equal to the activation entropy, (delta S(w)(double dagger))', of the corresponding solution reaction. It is, however, difficult to examine this idea by experimental approaches. The present paper defines the entropic proposal in a rigorous way and develops a computer simulation approach that determines (delta S(double dagger))'. This approach allows us to evaluate the differences between (delta S(double dagger))' of an enzymatic reaction and of the corresponding reference reaction in solution. Our approach is used in a study of the entropic contribution to the catalytic reaction of subtilisin. It is found that this contribution is much smaller than previously thought. This result is due to the following: (i) Many of the motions that are free in the reactants state of the reference solution reaction are also free at the transition state. (ii) The binding to the enzyme does not completely freeze the motion of the reacting fragments so that (delta S(double dagger))' in the enzymes is not zero. (iii) The binding entropy is not necessarily equal to (delta S(w)(double dagger))'.

Catalysis↗
PubMed2000

How does GAP catalyze the GTPase reaction of Ras? A computer simulation study.

The formation of a complex between p21(ras) and GAP accelerates the GTPase reaction of p21(ras) and terminates the signal for cell proliferation. The understanding of this rate acceleration is important for the elucidation of the role of Ras mutants in tumor formation. In principle there are two main options for the origin of the effect of GAP. One is a direct electrostatic interaction between the residues of GAP and the transition state of the Ras-GAP complex and the other is a GAP-induced shift of the structure of Ras to a configuration that increases the stabilization of the transition state. This work examines the relative importance of these options by computer simulations of the catalytic effect of Ras. The simulations use the empirical valence bond (EVB) method to study the GTPase reaction along the alternative associative and dissociative paths. This approach reproduces the trend in the overall experimentally observed catalytic effect of GAP: the calculated effect is 7 +/- 3 kcal/mol as compared to the observed effect of approximately 6.6 kcal/mol. Furthermore, the calculated effect of mutating Arg789 to a nonpolar residue is 3-4 kcal/mol as compared to the observed effect of 4.5 kcal/mol for the Arg789Ala mutation. It is concluded, in agreement with previous proposals, that the effect of Arg789 is associated with its direct interaction with the transition state charge distribution. However, calculations that use the coordinates of Ras from the Ras-GAP complex (referred to here as Ras') reproduce a significant catalytic effect relative to the Ras coordinates. This indicates that part of the effect of GAP involves a stabilization of a catalytic configuration of Ras. This configuration increases the positive electrostatic potential on the beta-phosphate (relative to the corresponding situation in the free Ras). In other words, GAP stabilizes the GDP bound configuration of Ras relative to that of the GTP-bound conformation. The elusive oncogenic effect of mutating Gln61 is also explored. The calculated effect of such mutations in the Ras-GAP complex are found to be small, while the observed effect is very large (8.7 kcal/mol). Since the Ras is locked in its Ras-GAP configuration in our simulations, we conclude that the oncogenic effect of mutation of Gln61 is indirect and is associated most probably with the structural changes of Ras upon forming the Ras-GAP complex. In view of these and the results for the Ras' we conclude that GAP activates Ras by both direct electrostatic stabilization of the transition state and an indirect allosteric effect that stabilizes the GDP-bound form. The present study also explored the feasibility of the associative and dissociative mechanism in the GTPase reaction of Ras. It is concluded that the reaction is most likely to involve an associative mechanism.

Computer Simulation↗
PubMed2000

Examining methods for calculations of binding free energies: LRA, LIE, PDLD-LRA, and PDLD/S-LRA calculations of ligands binding to an HIV protease.

Several strategies for evaluation of the protein-ligand binding free energies are examined. Particular emphasis is placed on the Linear Response Approximation (LRA) (Lee et. al., Prot Eng 1992;5:215-228) and the Linear Interaction Energy (LIE) method (Aqvist et. al., Prot Eng 1994;7:385-391). The performance of the Protein Dipoles Langevin Dipoles (PDLD) method and its semi-microscopic version (the PDLD/S method) is also considered. The examination is done by using these methods in the evaluating of the binding free energies of neutral C2-symmetric cyclic urea-based molecules to Human Immunodeficiency Virus (HIV) protease. Our starting point is the introduction of a thermodynamic cycle that decomposes the total binding free energy to electrostatic and non-electrostatic contributions. This cycle is closely related to the cycle introduced in our original LRA study (Lee et. al., Prot Eng 1992;5:215-228). The electrostatic contribution is evaluated within the LRA formulation by averaging the protein-ligand (and/or solvent-ligand) electrostatic energy over trajectories that are propagated on the potentials of both the polar and non-polar (where all residual charges are set to zero) states of the ligand. This average involves a scaling factor of 0.5 for the contributions from each state and this factor is being used in both the LRA and LIE methods. The difference is, however, that the LIE method neglects the contribution from trajectories over the potential of the non-polar state. This approximation is entirely valid in studies of ligands in water but not necessarily in active sites of proteins. It is found in the present case that the contribution from the non-polar states to the protein-ligand binding energy is rather small. Nevertheless, it is clearly expected that this term is not negligible in cases where the protein provides preorganized environment to stabilize the residual charges of the ligand. This contribution can be particularly important in cases of charged ligands. The analysis of the non-electrostatic term is much more complex. It is concluded that within the LRA method one has to complete the relevant thermodynamic cycle by evaluating the binding free energy of the "non-polar" ligand, l;, where all the residual charges are set to zero. It is shown that the LIE term, which involves the scaling of the van der Waals interaction by a constant beta (usually in the order of 0.15 to 0.25), corresponds to this part of the cycle. In order to elucidate the nature of this non-electrostatic term and the origin of the scaling constant beta, it is important to evaluate explicitly the different contributions to the binding energy of the non-polar ligand, DeltaG(bind,l;). Since this cannot be done at present (for relatively large ligands) by rigorous free energy perturbation approaches, we evaluate DeltaG(bind,l;) by the PDLD approach, augmented by microscopic calculations of the change in configurational entropy upon binding. This evaluation takes into account the van der Waals, hydrophobic, water penetration and entropic contributions, which are the most important free energy contributions that make up the total DeltaG(bind,l;). The sum of these contributions is scaled by a factor straight theta and it is argued that obtaining a quantitative balance between these contributions should result in straight theta = 1. By doing so we should have a reliable estimate of the value of the LIE beta and a way to understand its origin. The present approach gives straight theta values between 0.5 and 0.73, depending on the approximation used. This is encouraging but still not satisfying. Nevertheless, one might be able to use our PDLD approach to estimate the change of the LIE straight theta between different protein active sites. It is pointed out that the LIE method is quite similar to our original approach where the electrostatic term was evaluated by the LRA method and the non-electrostatic term by the PDLD method (with its vdw, solvation,

Energy Transfer↗
PubMed1999

Simulating proton translocations in proteins: probing proton transfer pathways in the Rhodobacter sphaeroides reaction center.

A general method for simulating proton translocations in proteins and for exploring the role of different proton transfer pathways is developed and examined. The method evaluates the rate constants for proton transfer processes using the energetics of the relevant proton configurations. The energies (DeltaG((m))) of the different protonation states are evaluated in two steps. First, the semimicroscopic version of the protein dipole Langevin dipole (PDLD/S) method is used to evaluate the intrinsic energy of bringing the protons to their protein sites, when the charges of all protein ionized residues are set to zero. Second, the interactions between the charged groups are evaluated by using a Coulomb's Law with an effective dielectric constant. This approach, which was introduced in an earlier study by one of the authors of the current report, allows for a very fast determination of any DeltaG((m)) and for practical evaluation of the time-dependent proton population: That is, the rate constants for proton transfer processes are evaluated by using the corresponding DeltaG((m)) values and a Marcus type relationship. These rate constants are then used to construct a master equation, the integration of which by a fourth-order Runge-Kutta method yields the proton population as a function of time. The integration evaluates, 'on the fly,' the changes of the rate constants as a result of the time-dependent changes in charge-charge interaction, and this feature benefits from the fast determination of DeltaG((m)). The method is demonstrated in a preliminary study of proton translocation processes in the reaction center of Rhodobacter sphaeroides. It is found that proton transfer across water chains involves significant activation barriers and that ionized protein residues probably are involved in the proton transfer pathways. The potential of the present method in analyzing mutation experiments is discussed briefly and illustrated. The present study also examines different views of the nature of proton translocations in proteins. It is shown that such processes are controlled mainly by the electrostatic interaction between the proton site and its surroundings rather than by the local bond rearrangements of water molecules that are involved in the proton pathways. Thus, the overall rate of proton transport frequently is controlled by the highest barrier along the conduction pathway. Proteins 1999;36:484-500.

Computer Simulation↗
PubMed1999

Role of active site residues in the glycosylase step of T4 endonuclease V. Computer simulation studies on ionization states.

T4 Endonuclease V (EndoV) is a base excision repair enzyme that removes thymine dimers (TD) from damaged DNA. To elucidate the role of the active site residues in catalysis, their pK(a)'s were evaluated using the semimicroscopic version of the protein dipoles-Langevin dipoles method (PDLD/S). Contributions of different effects to the pK(a) such as charge-charge interactions, conformational rearrangement, protein relaxation, and DNA binding were analyzed in detail. Charging of the active site residues was found to be less favorable in the complex than in the free enzyme. The pK(a) of the N-terminus decreased from 8.01 in the free enzyme to 6.52 in the complex, while the pK(a) of Glu-23 increased from 1. 52 to 7.82, which indicates that the key residues are neutral in the reactant state of the glycosylase step. These pK(a)'s are in agreement with the optimal pH range of the reaction and support the N-terminus acting as a nucleophile. The Glu-23 in its protonated form is hydrogen bonded to O4' of the sugar of 5' TD and can play a role in increasing the positive charge of C1' and, hence, accelerating the nucleophilic substitution. Furthermore, the neutral Glu-23 is a likely candidate to protonate O4' to induce ring opening required to complete the glycosylase step of EndoV. The positively charged Arg-22 and Arg-26 provide an electrostatically favorable environment for the leaving base. To distinguish between S(N)1 and S(N)2 mechanisms of the glycosylase step the energetics of protonating O2 of 5' TD was calculated. The enzyme was found to stabilize the neutral thymine by approximately 3.6 kcal/mol, whereas it destabilizes the protonated thymine by approximately 6.6 kcal/mol with respect to an aqueous environment. Consequently, the formation of a protonated thymine intermediate is unlikely, indicating an S(N)2 reaction mechanism for the glycosylase step.

Bacteriophage T4↗
PubMed1999

Mechanistic alternatives in phosphate monoester hydrolysis: what conclusions can be drawn from available experimental data?

Phosphate monoester hydrolysis reactions in enzymes and solution are often discussed in terms of whether the reaction pathway is associative or dissociative. Although experimental results for solution reactions have usually been considered as evidence for the second alternative, a closer thermodynamic analysis of observed linear free energy relationships shows that experimental information is consistent with the associative, concerted and dissociative alternatives.

Animals↗
PubMed1998

Computer simulations of enzyme catalysis: finding out what has been optimized by evolution.

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Animals↗
PubMed1998

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.

Computer Simulation↗
PubMed1998

Reorganization energy of the initial electron-transfer step in photosynthetic bacterial reaction centers.

The reorganization energy (lambda) for electron transfer from the primary electron donor (P*) to the adjacent bacteriochlorophyll (B) in photosynthetic bacterial reaction centers is explored by molecular-dynamics simulations. Relatively long (40 ps) molecular-dynamics trajectories are used, rather than free energy perturbation techniques. When the surroundings of the reaction center are modeled as a membrane, lambda for P* B --> P+ B- is found to be approximately 1.6 kcal/mol. The results are not sensitive to the treatment of the protein's ionizable groups, but surrounding the reaction center with water gives higher values of lambda (approximately 6.5 kcal/mol). In light of the evidence that P+ B- lies slightly below P* in energy, the small lambda obtained with the membrane model is consistent with the speed and temperature independence of photochemical charge separation. The calculated reorganization energy is smaller than would be expected if the molecular dynamics trajectories had sampled the full conformational space of the system. Because the system does not relax completely on the time scale of electron transfer, the lambda obtained here probably is more pertinent than the larger value that would be obtained for a fully equilibrated system.

Calorimetry↗
PubMed1998

The effect of protein relaxation on charge-charge interactions and dielectric constants of proteins.

The effect of the reorganization of the protein polar groups on charge-charge interaction and the corresponding effective dielectric constant (epsilon(eff)) is examined by the semimicroscopic version of the Protein Dipole Langevin Dipoles (PDLD/S) method within the framework of the Linear Response Approximation (LRA). This is done by evaluating the interactions between ionized residues in the reaction center of Rhodobacter sphaeroides, while taking into account the protein reorganization energy. It is found that an explicit consideration of the protein relaxation leads to a significant increase in epsilon(eff) and that semimicroscopic models that do not take this relaxation into account force one to use a large value for the so-called "protein dielectric constant," epsilon(p), of the Poisson-Boltzmann model or for the corresponding epsilon(in) in the PDLD/S model. An additional increase in epsilon(eff) is expected from the reorganization of ionized residues and from changes in the degree of water penetration. This finding provides further support for the idea that epsilon(in) (or epsilon(p)) represents contributions that are not considered explicitly. The present study also provides a systematic illustration of the nature of epsilon(eff), supporting our previously reported view that charge-charge interactions correspond to a large value of this "dielectric constant," even in protein interiors. It is also pointed out that epsilon(eff) for the interaction between ionizable groups in proteins is very different from the effective dielectric constant, epsilon'(eff), that determines the free energy of ion pairs in proteins (epsilon'(eff) reflects the effect of preoriented protein dipoles). Finally, the problems associated with the search for a general epsilon(in) are discussed. It is clarified that the epsilon(in) that reproduces the effect of protein relaxation on charge-charge interaction is not equal to the epsilon(in) that reproduces the corresponding effect upon formation of individual charges. This reflects fundamental inconsistencies in attempts to cast microscopic concepts in a macroscopic model. Thus one should either use a large epsilon(in) for charge-charge interactions and a small epsilon(in) for charge-dipole interactions or consider the protein relaxation microscopically.

Binding Sites↗
PubMed1998

Electrostatic effects in macromolecules: fundamental concepts and practical modeling.

The past few years have seen an exponential growth in the calculations of electrostatic energies of macromolecules and an increased recognition of the crucial role of electrostatic effects. This review considers the current state of the field. Focus is placed on calculations of pKas, redox potentials and binding energies in macromolecules and clarification of the fact that the value of the dielectric 'constant' of a protein depends on its definition and that small dielectric constants should not be used in describing charge-charge interactions by current continuum models.

Enzymes↗
PubMed1997

A fast estimate of electrostatic group contributions to the free energy of protein-inhibitor binding.

Dissecting ligand-protein binding free energies in individual contributions of protein residues (which are referred to here as 'group contributions') is of significant importance. For example, such contributions could help in estimating the corresponding mutational effects and in studies of drug resistance problems. However, the meaning of group contributions is not always uniquely defined and the approximations for rapid estimates of such contributions are not well developed. In this paper, the nature of group contributions to binding free energy is examined, focusing particularly on electrostatic contributions which are expected to be well behaved. This analysis examines different definitions of group contributions; the 'relaxed' group contributions that represent the change in binding energy upon mutation of the given residue to glycine, and the 'non-relaxed' group contributions that represent the scaled Coulomb interaction between the given residue and the ligand. Both contributions are defined and evaluated by the linear response approximation (LRA) of the PDLD/ S method. The present analysis considers the binding of pepstatin to endothiapepsin and 23 of its mutants as a test case for a neutral ligand. The 'non-relaxed' group contributions of 15 endothiapepsin residues show significant peaks in the 'electrostatic fingerprint'. The residues that contribute to the electrostatic fingerprint are located in the binding site of endothiapepsin. They include the aspartic dyad (Asp32, Asp215) with adjacent residues and the flap region. Twelve of these 15 residues have a heavy atom distance of <3.75 A to pepstatin. The contributions of 8 (10) of these 12 residues can be reconciled with the calculated 'relaxed' group contributions where one allows the protein and solvent (solvent only) to relax upon mutation of the given residue to glycine. On the other hand, it was found that residues at the second 'solvation shell' can have relaxed contributions that are not captured by the non-relaxed approach. Hence, whereas residues with significant non-relaxed electrostatic contributions are likely to contribute to binding, residues with small non-relaxed contributions may still affect the binding energy. At any rate, it is established here that even in the case of uncharged inhibitors it is possible to use the non-relaxed electrostatic fingerprint to detect 'hot' residues that are responsible for binding. This is significant since some versions of the non-relaxed approximation are faster by several orders of magnitude than more rigorous approaches. The general applicability of this approach is outlined, emphasizing its potential in studies of drug resistance where it is crucial to have a rapid way of anticipating the effect of mutation on both drug binding and catalysis.

Aspartic Acid↗
PubMed1996

Energy considerations show that low-barrier hydrogen bonds do not offer a catalytic advantage over ordinary hydrogen bonds.

Low-barrier hydrogen bonds have recently been proposed as a major factor in enzyme catalysis. Here we evaluate the feasibility of transition state (TS) stabilization by low-barrier hydrogen bonds in enzymes. Our analysis focuses on the facts that (i) a low-barrier hydrogen bond is less stable than a regular hydrogen bond in water, (ii) TSs are more stable in the enzyme active sites than in water, and (iii) a nonpolar active site would destabilize the TS relative to its energy in water. Combining these points and other experimental and theoretical facts in a physically consistent frame-work shows that a low-barrier hydrogen bond cannot stabilize the TS more than an ordinary hydrogen bond. The reason for the large catalytic effect of active site hydrogen bonds is that their formation entails a lower reorganization energy than their solution counterparts, due to the preorganized enzyme environment.

Binding Sites↗
PubMed1996

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.

Binding Sites↗
PubMed1996

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.

Catalysis↗
PubMed1996

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.

Binding Sites↗