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

Classification of KRAS-Activating Mutations and the Implications for Therapeutic Intervention

Members of the family of RAS proto-oncogenes, discovered just over 40 years ago, were among the first cancer-initiating genes to be discovered. Of the three RAS family members, KRAS is the most frequently mutated in human cancers. Despite intensive biological and biochemical study of RAS proteins over the past four decades, we are only now starting to devise therapeutic strategies to target their oncogenic properties. Here, we highlight the distinct biochemical properties of common and rare KRAS alleles, enabling their classification into functional subtypes. We also discuss the implications of this functional classification for potential therapeutic avenues targeting mutant subtypes.

SIGNIFICANCE

Efforts in the recent past to inhibit KRAS oncogenicity have focused on kinases that function in downstream signal transduction cascades, although preclinical successes have not translated to patients with KRAS-mutant cancer. Recently, clinically effective covalent inhibitors of KRASG12C have been developed, establishing two principles that form a foundation for future efforts. First, KRAS is druggable. Second, each mutant form of KRAS is likely to have properties that make it uniquely druggable.

What Does the Brønsted Slope Measure in the Phosphoryl Transfer Transition State?

The structural and energetic features of phosphate and phosphonate hydrolysis in Protein Phosphatase-1 (PP1) and water are studied using quantum mechanical (QM) cluster models. The calculations are able to reproduce observed kinetic isotope effects and capture several key trends in the experimental Brønsted plots: the β l g values are rather different for phosphate and phosphonate ester hydrolysis in solution but are similar in PP1. Detailed analyses of structure, charge distribution and bond order of computed transition states support the general conclusion from experimental study that phosphoryl transfer transition states are different for the two classes of substrates in solution but similar in PP1. On the other hand, the microscopic models also highlight notable differences between the phosphate and phosphonate transition states, which are manifested in not only structure but also kinetic isotope effects. Overall, we find that while β l g / β E Q , l g generally correlates with the partial charge on leaving group oxygen and the fractional bond order of the breaking P- O l g bond, the precise mapping between β l g / β E Q , l g and P- O l g bond order in the transition state is difficult due largely to the cross talk between breaking and forming P-O bonds. Therefore, further supporting previous analyses of limitations of free energy relations, our results suggest that while free energy relation is a valuable tool for probing the nature of transition state, a quantitative mapping of β l g and β l g / β E Q , l g values to structure or charge in the transition state should be conducted with great care.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Specific Substates of Ras To Interact with GAPs and Effectors: Revealed by Theoretical Simulations and FTIR Experiments

The oncogenic Ras protein adopts various specific conformational states to execute its function in signal transduction. The large number of Ras structures obtained from X-ray and NMR experiments illustrates the diverse conformations that Ras adopts. It is difficult, however, to connect specific structural features with Ras functions. We report the free-energy landscape of Ras·GTP based on extensive explicit solvent simulations. The free-energy map clearly shows that the functional state 2 of Ras·GTP in fact has two distinct substates, denoted here as "Tyr32in" and "Tyr32out". Unbiased MD simulations show that the two substrates interconvert on the submicrosecond scale in solution, pointing to a novel mechanism for Ras·GTP to selectively interact with GAPs and effectors. This proposal is further supported by time-resolved FTIR experiments, which demonstrate that Tyr32 destabilizes the Ras·GAP complex and facilitates an efficient termination of Ras signaling.

The protonation states of GTP and GppNHp in Ras proteins

The small GTPase Ras transmits signals in a variety of cellular signaling pathways, most prominently in cell proliferation. GTP hydrolysis in the active center of Ras acts as a prototype for many GTPases and is the key to the understanding of several diseases, including cancer. Therefore, Ras has been the focus of intense research over the last decades. A recent neutron diffraction crystal structure of Ras indicated a protonated γ-guanylyl imidodiphosphate (γ-GppNHp) group, which has put the protonation state of GTP in question. A possible protonation of GTP was not considered in previously published mechanistic studies. To determine the detailed prehydrolysis state of Ras, we calculated infrared and NMR spectra from quantum mechanics/molecular mechanics (QM/MM) simulations and compared them with those from previous studies. Furthermore, we measured infrared spectra of GTP and several GTP analogs bound to lipidated Ras on a membrane system under near-native conditions. Our findings unify results from previous studies and indicate a structural model confirming the hypothesis that γ-GTP is fully deprotonated in the prehydrolysis state of Ras.

Review: Translational GTPases

Translational GTPases (trGTPases) play key roles in facilitating protein synthesis on the ribosome. Despite the high degree of evolutionary conservation in the sequences of their GTP-binding domains, the rates of GTP hydrolysis and nucleotide exchange vary broadly between different trGTPases. EF-Tu, one of the best-characterized model G proteins, evolved an exceptionally rapid and tightly regulated GTPase activity, which ensures rapid and accurate incorporation of amino acids into the nascent chain. Other trGTPases instead use the energy of GTP hydrolysis to promote movement or to ensure the forward commitment of translation reactions. Recent data suggest the GTPase mechanism of EF-Tu and provide an insight in the catalysis of GTP hydrolysis by its unusual activator, the ribosome. Here we summarize these advances in understanding the functional cycle and the regulation of trGTPases, stimulated by the elucidation of their structures on the ribosome and the progress in dissecting the reaction mechanism of GTPases. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 463-475, 2016.

Neutron Crystal Structure of RAS GTPase Puts in Question the Protonation State of the GTP γ-Phosphate

RAS GTPase is a prototype for nucleotide-binding proteins that function by cycling between GTP and GDP, with hydrogen atoms playing an important role in the GTP hydrolysis mechanism. It is one of the most well studied proteins in the superfamily of small GTPases, which has representatives in a wide range of cellular functions. These proteins share a GTP-binding pocket with highly conserved motifs that promote hydrolysis to GDP. The neutron crystal structure of RAS presented here strongly supports a protonated γ-phosphate at physiological pH. This counters the notion that the phosphate groups of GTP are fully deprotonated at the start of the hydrolysis reaction, which has colored the interpretation of experimental and computational data in studies of the hydrolysis mechanism. The neutron crystal structure presented here puts in question our understanding of the pre-catalytic state associated with the hydrolysis reaction central to the function of RAS and other GTPases.

Overview of simulation studies on the enzymatic activity and conformational dynamics of the GTPase Ras

Over the last 40 years, we have learnt a great deal about the Ras onco-proteins. These intracellular molecular switches are essential for the function of a variety of physiological processes, including signal transduction cascades responsible for cell growth and proliferation. Molecular simulations and free energy calculations have played an essential role in elucidating the conformational dynamics and energetics underlying the GTP hydrolysis reaction catalysed by Ras. Here we present an overview of the main lessons from molecular simulations on the GTPase reaction and conformational dynamics of this important anti-cancer drug target. In the first part, we summarise insights from quantum mechanical and combined quantum mechanical/molecular mechanical simulations as well as other free energy methods and highlight consensus viewpoints as well as remaining controversies. The second part provides a very brief overview of new insights emerging from large-scale molecular dynamics simulations. We conclude with a perspective regarding future studies of Ras where computational approaches will likely play an active role.

Transition state in DNA polymerase β catalysis: rate-limiting chemistry altered by base-pair configuration

Kinetics studies of dNTP analogues having pyrophosphate-mimicking β,γ-pCXYp leaving groups with variable X and Y substitution reveal striking differences in the chemical transition-state energy for DNA polymerase β that depend on all aspects of base-pairing configurations, including whether the incoming dNTP is a purine or pyrimidine and if base-pairings are right (T•A and G•C) or wrong (T•G and G•T). Brønsted plots of the catalytic rate constant (log(k_pol)) versus pK_a4 for the leaving group exhibit linear free energy relationships (LFERs) with negative slopes ranging from −0.6 to −2.0, consistent with chemical rate-determining transition-states in which the active-site adjusts to charge-stabilization demand during chemistry depending on base-pair configuration. The Brønsted slopes as well as the intercepts differ dramatically and provide the first direct evidence that dNTP base recognition by the enzyme–primer–template complex triggers a conformational change in the catalytic region of the active-site that significantly modifies the rate-determining chemical step.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

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

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

Structural characterization of the Mycobacterium tuberculosis biotin biosynthesis enzymes 7,8-diaminopelargonic acid synthase and dethiobiotin synthetase

Mycobacterium tuberculosis (Mtb) depends on biotin synthesis for survival during infection. In the absence of biotin, disruption of the biotin biosynthesis pathway results in cell death rather than growth arrest, an unusual phenotype for an Mtb auxotroph. Humans lack the enzymes for biotin production, making the proteins of this essential Mtb pathway promising drug targets. To this end, we have determined the crystal structures of the second and third enzymes of the Mtb biotin biosynthetic pathway, 7,8-diaminopelargonic acid synthase (DAPAS) and dethiobiotin synthetase (DTBS), at respective resolutions of 2.2 and 1.85 A. Superimposition of the DAPAS structures bound either to the SAM analogue sinefungin or to 7-keto-8-aminopelargonic acid (KAPA) allowed us to map the putative binding site for the substrates and to propose a mechanism by which the enzyme accommodates their disparate structures. Comparison of the DTBS structures bound to the substrate 7,8-diaminopelargonic acid (DAPA) or to ADP and the product dethiobiotin (DTB) permitted derivation of an enzyme mechanism. There are significant differences between the Mtb enzymes and those of other organisms; the Bacillus subtilis DAPAS, presented here at a high resolution of 2.2 A, has active site variations and the Escherichia coli and Helicobacter pylori DTBS have alterations in their overall folds. We have begun to exploit the unique characteristics of the Mtb structures to design specific inhibitors against the biotin biosynthesis pathway in Mtb.

Expression, purification, and characterization of soluble K-Ras4B for structural analysis

A p21 GTPase K-Ras4B plays an important role in human cancer and represents an excellent target for cancer therapeutics. Currently, there are no drugs directly targeting K-Ras4B. In part, this is due to the lack of structural information describing unique features of K-Ras4B. Here we describe a methodology allowing production of soluble, well-folded K-Ras4B for structural analysis. The key points in K-Ras4B preparation are low temperature expression and extraction of K-Ras4B from the insoluble fraction using a nucleotide loading procedure in the presence of Mg(2+) and citrate, a low affinity chelator. Additionally, a significant amount of K-Ras4B could be extracted from the soluble fraction. We show that recombinant K-Ras4B is monomeric in solution. Excellent NMR signal dispersion suggests that the protein is well-folded and is amenable to solution structure determination. In addition, using phospholipid bilayer nanodiscs we show that recombinant K-Ras4B interacts with lipids and that this interaction is mediated by the C-terminal hypervariable region.

The EVB as a quantitative tool for formulating simulations and analyzing biological and chemical reactions

Recent years have seen dramatic improvements in computer power, allowing ever more challenging problems to be approached. In light of this, it is imperative to have a quantitative model for examining chemical reactivity, both in the condensed phase and in solution, as well as to accurately quantify physical organic chemistry (particularly as experimental approaches can often be inconclusive). Similarly, computational approaches allow for great progress in studying enzyme catalysis, as they allow for the separation of the relevant energy contributions to catalysis. Due to the complexity of the problems that need addressing, there is a need for an approach that can combine reliability with an ability to capture complex systems in order to resolve long-standing controversies in a unique way. Herein, we will demonstrate that the empirical valence bond (EVB) approach provides a powerful way to connect the classical concepts of physical organic chemistry to the actual energies of enzymatic reactions by means of computation. Additionally, we will discuss the proliferation of this approach, as well as attempts to capture its basic chemistry and repackage it under different names. We believe that the EVB approach is the most powerful tool that is currently available for studies of chemical processes in the condensed phase in general and enzymes in particular, particularly when trying to explore the different proposals about the origin of the catalytic power of enzymes.

Transition state structures and the roles of catalytic residues in GAP-facilitated GTPase of Ras as elucidated by (18)O kinetic isotope effects

Ras-catalyzed guanosine 5' triphosphate (GTP) hydrolysis proceeds through a loose transition state as suggested in our previous study of (18)O kinetic isotope effects (KIE) [ Du , X. et al. ( 2004 ) Proc. Natl. Acad. Sci. U.S.A. 101 , 8858 - 8863 ]. To probe the mechanisms of GTPase activation protein (GAP)-facilitated GTP hydrolysis reactions, we measured the (18)O KIEs in GTP hydrolysis catalyzed by Ras in the presence of GAP(334) or NF1(333), the catalytic fragment of p120GAP or NF1. The KIEs in the leaving group oxygens (the beta nonbridge and the beta-gamma bridge oxygens) reveal that chemistry is rate-limiting in GAP(334)-facilitated GTP hydrolysis but only partially rate-limiting in the NF1(333)-facilitated GTP hydrolysis reaction. The KIEs in the gamma nonbridge oxygens and the leaving group oxygens reveal that the GAP(334) or NF1(333)-facilitated GTP hydrolysis reaction proceeds through a loose transition state that is similar in nature to the transition state of the GTP hydrolysis catalyzed by Ras alone. However, the KIEs in the pro-S beta, pro-R beta, and beta-gamma oxygens suggest that charge increase on the beta-gamma bridge oxygen is more prominent in the transition states of GAP(334)- and NF1(333)-facilitated reactions than that catalyzed by the intrinsic GTPase activity of Ras. The charge distribution on the two beta nonbridge oxygens is also very asymmetric. The catalytic roles of active site residues were inferred from the effect of mutations on the reaction rate and KIEs. Our results suggest that the arginine finger of GAP and amide protons in the P-loop of Ras stabilize the negative charge on the beta-gamma bridge oxygen and the pro-S beta nonbridge oxygen of a loose transition state, whereas Lys-16 of Ras and Mg(2+) are only involved in substrate binding.

Dineopentyl phosphate hydrolysis: evidence for stepwise water attack

Phosphate ester hydrolysis is ubiquitous in biology, playing a central role in energy production, signaling, biosynthesis, and the regulation of protein function among other things. Although the mechanism of action of the enzymes regulating this reaction has been the focus of intensive research in the past few decades, the correct description of this apparently simple reaction remains controversial. A clear understanding of the mechanism that takes place in solution is crucial to be able to evaluate whether proposals for the enzyme-catalyzed mechanisms are reasonable. For the pH-independent hydrolysis of phosphate diesters, several kinetically equivalent mechanisms are plausible, including hydroxide attack on the neutral phosphate. However, it is very difficult to measure the rate of this reaction directly by experimental methods, so it has been evaluated by examining the rate of hydrolysis of neutral phosphate triesters, where a methyl group has replaced a proton. This may not be an accurate model of the neutral phosphate diester and does not provide information about a reaction pathway that is concerted with nucleophilic attack to generate a similar phosphorane. We have carefully mapped out free energy surfaces for both hydroxide and water attack on the dineopentyl phosphate anion and for water attack on the neutral diester. In doing so, we have accurately reproduced existing experimental data and demonstrate that water attack proceeds through an associative mechanism with proton transfer to the phosphate to generate a phosphorane intermediate. Our data show that the substrate-as-base mechanism is viable for phosphate ester hydrolysis, which may have important implications for the studies of phosphate ester hydrolysis by enzymes.

On the interpretation of the observed linear free energy relationship in phosphate hydrolysis: a thorough computational study of phosphate diester hydrolysis in solution

The hydrolysis of phosphate esters is crucially important to biological systems, being involved in, among other things, signaling, energy transduction, biosynthesis, and the regulation of protein function. Despite this, there are many questions that remain unanswered in this important field, particularly with regard to the preferred mechanism of hydrolysis of phosphate esters, which can proceed through any of multiple pathways that are either associative or dissociative in nature. Previous comparisons of calculated and observed linear free energy relationships (LFERs) for phosphate monoester dianions with different leaving groups showed that the TS character gradually changes from associative to dissociative with the increasing acidity of the leaving group, while reproducing the experimental LFER. Here, we have generated ab initio potential energy surfaces for the hydrolysis of phosphate diesters in solution, with a variety of leaving groups. Once again, the reaction changes from a compact concerted pathway to one that is more expansive in character when the acidity of the leaving group increases. When such systems are examined in solution, it is essential to take into consideration the contribution of solute to the overall activation entropy, which remains a major computational challenge. The popular method of calculating the entropy using a quasi-harmonic approximation appears to markedly overestimate the configurational entropy for systems with multiple occupied energy wells. We introduce an improved restraint release approach for evaluating configurational entropies and apply this approach to our systems. We demonstrate that when this factor is taken into account, it is possible to reproduce the experimental LFER for this system with reasonable accuracy.

Using the constrained DFT approach in generating diabatic surfaces and off diagonal empirical valence bond terms for modeling reactions in condensed phases

The empirical valence bond (EVB) model provides an extremely powerful way for modeling and analyzing chemical reactions in solutions and proteins. However, this model is based on the unverified assumption that the off diagonal elements of the EVB Hamiltonian do not change significantly upon transfer of the reacting system from one phase to another. This ad hoc assumption has been rationalized by its consistency with empirically observed linear free energy relationships, as well as by other qualitative considerations. Nevertheless, this assumption has not been rigorously established. The present work explores the validity of the above EVB key assumption by a rigorous numerical approach. This is done by exploiting the ability of the frozen density functional theory (FDFT) and the constrained density functional theory (CDFT) models to generate convenient diabatic states for QM/MM treatments, and thus to examine the relationship between the diabatic and adiabatic surfaces, as well as the corresponding effective off diagonal elements. It is found that, at least for the test case of S(N)()2 reactions, the off diagonal element does not change significantly upon moving from the gas phase to solutions and thus the EVB assumption is valid and extremely useful.

H and other transfers in enzymes and in solution: theory and computations, a unified view. 2. Applications to experiment and computations

Equations obtained in part I for the free-energy barrier to one-step enzymatic reactions between bound reactants are discussed. The rate is expressed in terms of lambdao (protein reorganization energy), DeltaG(o) (standard free energy of reaction of the H-transfer step), bond breaking/bond forming term, w (work terms), and H-transmission property. Two alternative approximations for the coupling of the bond breaking/bond forming and protein are distinguished experimentally in favorable cases by the DeltaG(o) where the maximum deuterium kinetic isotope effect occurs. Plots of log rate versus DeltaG(o) and properties such as DeltaS* and DeltaS(o) are discussed. The weak or zero T-dependence of the kinetic isotope effect for wild-type enzymes operating under physiological conditions is interpreted in terms of vanishing (or isotopically insensitive) w plus transfer from the lowest H-state. Static and dynamic protein flexibility is discussed. While the many correlations accessible for electron transfers are not available for H-transfers in enzymes, a combination of experiment, computation, and analytical approaches can assist in evaluating the utility of the present equations and in suggesting further experiments and computations. A protein reorganization energy lambdao is obtained in the literature from the extended valence bond formalism where diabatic electronic states are used. A method is suggested for extracting it when instead a bond distance difference coordinate is used. The results may provide a bridge between the two approaches.

Mechanisms of guanosine triphosphate hydrolysis by Ras and Ras-GAP proteins as rationalized by ab initio QM/MM simulations

The hydrolysis reaction of guanosine triphosphate (GTP) by p21(ras) (Ras) has been modeled by using the ab initio type quantum mechanical-molecular mechanical simulations. Initial geometry configurations have been prompted by atomic coordinates of the crystal structure (PDBID: 1QRA) corresponding to the prehydrolysis state of Ras in complex with GTP. Multiple searches of minimum energy geometry configurations consistent with the hydrogen bond networks have been performed, resulting in a series of stationary points on the potential energy surface for reaction intermediates and transition states. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step mechanism of GTP hydrolysis. At the first stage, a unified action of the nearest residues of Ras and the nearest water molecules results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low barrier (16.7 kcal/mol) transition state TS1. At the second stage, the inorganic phosphate is formed in consequence of proton transfers mediated by two water molecules and assisted by the Gln61 residue from Ras. The highest transition state at this segment, TS3, is estimated to have an energy 7.5 kcal/mol above the enzyme-substrate complex. The results of simulations are compared to the previous findings for the GTP hydrolysis in the Ras-GAP (p21(ras)-p120(GAP)) protein complex. Conclusions of the modeling lead to a better understanding of the anticatalytic effect of cancer causing mutation of Gln61 from Ras, which has been debated in recent years.

Modifying the beta,gamma leaving-group bridging oxygen alters nucleotide incorporation efficiency, fidelity, and the catalytic mechanism of DNA polymerase beta

DNA polymerase catalysis and fidelity studies typically compare incorporation of "right" versus "wrong" nucleotide bases where the leaving group is pyrophosphate. Here we use dGTP analogues replacing the beta,gamma-bridging O with CH2, CHF, CF2, or CCl2 to explore leaving-group effects on the nucleotidyl transfer mechanism and fidelity of DNA polymerase (pol) beta. T.G mismatches occur with fidelities similar to dGTP with the exception of the CH2 analogue, which is incorporated with 5-fold higher fidelity. All analogues are observed to bind opposite template C with Kds between 1 and 4 microM, and structural evidence suggests that the analogues bind in essentially the native conformation, making them suitable substrates for probing linear free energy relationships (LFERs) in transient-kinetics experiments. Importantly, Brnsted correlations of log(kpol) versus leaving-group pKa for both right and wrong base incorporation reveal similar sensitivities (betalg approximately -0.8) followed by departures from linearity, suggesting that a chemical step rather than enzyme conformational change is rate-limiting for either process. The location of the breaks relative to pKas of CF2, O, and the sterically bulky CCl2-bridging compounds suggests a modification-induced change in the mechanism by stabilization of leaving-group elimination. The results are addressed theoretically in terms of the energetics of successive primer 3'-O addition (bond forming) and pyrophosphate analogue elimination (bond breaking) reaction energy barriers.

Structural basis of effector regulation and signal termination in heterotrimeric Galpha proteins

This chapter addresses, from a molecular structural perspective gained from examination of x-ray crystallographic and biochemical data, the mechanisms by which GTP-bound Galpha subunits of heterotrimeric G proteins recognize and regulate effectors. The mechanism of GTP hydrolysis by Galpha and rate acceleration by GAPs are also considered. The effector recognition site in all Galpha homologues is formed almost entirely of the residues extending from the C-terminal half of alpha2 (Switch II) together with the alpha3 helix and its junction with the beta5 strand. Effector binding does not induce substantial changes in the structure of Galpha*GTP. Effectors are structurally diverse. Different effectors may recognize distinct subsets of effector-binding residues of the same Galpha protein. Specificity may also be conferred by differences in the main chain conformation of effector-binding regions of Galpha subunits. Several Galpha regulatory mechanisms are operative. In the regulation of GMP phospodiesterase, Galphat sequesters an inhibitory subunit. Galphas is an allosteric activator and inhibitor of adenylyl cyclase, and Galphai is an allosteric inhibitor. Galphaq does not appear to regulate GRK, but is rather sequestered by it. GTP hydrolysis terminates the signaling state of Galpha. The binding energy of GTP that is used to stabilize the Galpha:effector complex is dissipated in this reaction. Chemical steps of GTP hydrolysis, specifically, formation of a dissociative transition state, is rate limiting in Ras, a model G protein GTPase, even in the presence of a GAP; however, the energy of enzyme reorganization to produce a catalytically active conformation appears to be substantial. It is possible that the collapse of the switch regions, associated with Galpha deactivation, also encounters a kinetic barrier, and is coupled to product (Pi) release or an event preceding formation of the GDP*Pi complex. Evidence for a catalytic intermediate, possibly metaphosphate, is discussed. Galpha GAPs, whether exogenous proteins or effector-linked domains, bind to a discrete locus of Galpha that is composed of Switch I and the N-terminus of Switch II. This site is immediately adjacent to, but does not substantially overlap, the Galpha effector binding site. Interactions of effectors and exogenous GAPs with Galpha proteins can be synergistic or antagonistic, mediated by allosteric interactions among the three molecules. Unlike GAPs for small GTPases, Galpha GAPs supply no catalytic residues, but rather appear to reduce the activation energy for catalytic activation of the Galpha catalytic site.

Dissection of the GTPase mechanism of Ras protein by MD analysis of Ras mutants

Controlling the hydrolysis rate of GTP bound to the p21ras protein is crucial for the delicate timing of many biological processes. A few mechanisms were suggested for the hydrolysis of GTP. To gain more insight into the individual elementary events of GTP hydrolysis, we carried out molecular dynamic analysis of wild-type p21ras and some of its mutants. It was recently shown that Ras-related proteins and mutants generally follow a linear free energy relationship (LFER) relating the rate of reaction to the pK(a) of the gamma-phosphate group of the bound GTP, indicating that proton transfer from the attacking water to the GTP is the first elementary event in the GTPase mechanism. However, some exceptions were observed. Thus, the Gly12 --> Aspartic p21ras (G12D) mutant had a very low GTPase activity although its pK(a) was very close to that of the wild-type ras. Here we compared the molecular dynamics (MD) of wild-type Ras and G12D, showing that in the mutant the catalytic water molecule is displaced to a position where proton transfer to GTP is unfavorable. These results suggest that the mechanism of GTPase is indeed composed of an initial proton abstraction from water by the GTP, followed by a nucleophilic attack of the hydroxide ion on the gamma-phosphorus of GTP.

Crystallographic evidence for substrate-assisted GTP hydrolysis by a small GTP binding protein

GTP hydrolysis by small GTP binding proteins of the Ras superfamily is a universal reaction that controls multiple cellular regulations. Its enzymic mechanism has been the subject of long-standing debates as to the existence/identity of the general base and the electronic nature of its transition state. Here we report the high-resolution crystal structure of a small GTP binding protein, Rab11, solved in complex with GDP and Pi. Unexpectedly, a Pi oxygen and the GDP-cleaved oxygen are located less than 2.5 A apart, suggesting that they share a proton, likely in the form of a low-barrier hydrogen bond. This implies that the gamma-phosphate of GTP was protonated; hence, that GTP acts as a general base. Furthermore, this interaction should establish at, and stabilize, the transition state. Altogether, we propose a revised model for the GTPase reaction that should reconcile earlier models into a unique substrate-assisted mechanism.

Conformational states of Ras complexed with the GTP analogue GppNHp or GppCH2p: implications for the interaction with effector proteins

The guanine nucleotide-binding protein Ras occurs in solution in two different states, state 1 and state 2, when the GTP analogue GppNHp is bound to the active center as detected by (31)P NMR spectroscopy. Here we show that Ras(wt).Mg(2+).GppCH(2)p also exists in two conformational states in dynamic equilibrium. The activation enthalpy DeltaH(++)(12) and the activation entropy DeltaS(++)(12) for the transition from state 1 to state 2 are 70 kJ mol(-1) and 102 J mol(-1) K(-1), within the limits of error identical to those determined for the Ras(wt).Mg(2+).GppNHp complex. The same is true for the equilibrium constants K(12) = [2]/[1] of 2.0 and the corresponding DeltaG(12) of -1.7 kJ mol(-1) at 278 K. This excludes a suggested specific effect of the NH group of GppNHp on the equilibrium. The assignment of the phosphorus resonance lines of the bound analogues has been done by two-dimensional (31)P-(31)P NOESY experiments which lead to a correction of the already reported assignments of bound GppNHp. Mutation of Thr35 in Ras.Mg(2+).GppCH(2)p to serine leads to a shift of the conformational equilibrium toward state 1. Interaction of the Ras binding domain (RBD) of Raf kinase or RalGDS with Ras(wt) or Ras(T35S) shifts the equilibrium completely to state 2. The (31)P NMR experiments suggest that, besides the type of the side chain of residue 35, a main contribution to the conformational equilibrium in Ras complexes with GTP and GTP analogues is the effective acidity of the gamma-phosphate group of the bound nucleotide. A reaction scheme for the Ras-effector interaction is presented which includes the existence of two conformations of the effector loop and a weak binding state.

Quantum chemical modeling of the GTP hydrolysis by the RAS-GAP protein complex

We present results of the modeling for the hydrolysis reaction of guanosine triphosphate (GTP) in the RAS-GAP protein complex using essentially ab initio quantum chemistry methods. One of the approaches considers a supermolecular cluster composed of 150 atoms at a consistent quantum level. Another is a hybrid QM/MM method based on the effective fragment potential technique, which describes interactions between quantum and molecular mechanical subsystems at the ab initio level of the theory. Our results show that the GTP hydrolysis in the RAS-GAP protein complex can be modeled by a substrate-assisted catalytic mechanism. We can locate a configuration on the top of the barrier corresponding to the transition state of the hydrolysis reaction such that the straightforward descents from this point lead either to reactants GTP+H(2)O or to products guanosine diphosphate (GDP)+H(2)PO(4)(-). However, in all calculations such a single-step process is characterized by an activation barrier that is too high. Another possibility is a two-step reaction consistent with formation of an intermediate. Here the Pgamma-O(Pbeta) bond is already broken, but the lytic water molecule is still in the pre-reactive state. We present arguments favoring the assumption that the first step of the GTP hydrolysis reaction in the RAS-GAP protein complex may be assigned to the breaking of the Pgamma-O(Pbeta) bond prior to the creation of the inorganic phosphate.

Kinetic isotope effects in Ras-catalyzed GTP hydrolysis: evidence for a loose transition state

A remote labeling method has been developed to determine (18)O kinetic isotope effects (KIEs) in Ras-catalyzed GTP hydrolysis. Substrate mixtures consist of (13)C-depleted GTP and [(18)O,(13)C]GTP that contains (18)O at phosphoryl positions of mechanistic interest and (13)C at all carbon positions of the guanosine moiety. Isotope ratios of the nonvolatile substrates and products are measured by using a chemical reaction interface/isotope ratio mass spectrometer. The isotope effects are 1.0012 (0.0026) in the gamma nonbridge oxygens, 1.0194 (0.0025) in the leaving group oxygens (the beta-gamma oxygen and the two beta nonbridge oxygens), and 1.0105 (0.0016) in the two beta nonbridge oxygens. The KIE in the beta-gamma bridge oxygen was computed to be 1.0116 or 1.0088 by two different methods. The significant KIE in the leaving group reveals that chemistry is largely rate-limiting whereas the KIEs in the gamma nonbridge oxygens and the leaving group indicate a loose transition state that approaches a metaphosphate. The KIE in the two beta nonbridge oxygens is roughly equal to that in the beta-gamma bridge oxygen. This indicates that, in the transition state, Ras shifts one-half of the negative charge that arises from P(gamma)-O(beta-gamma) fission from the beta-gamma bridge oxygen to the two beta nonbridge oxygens. The KIE effects, interpreted in light of structural and spectroscopic data, suggest that Ras promotes a loose transition state by stabilizing negative charge in the beta-gamma bridge and beta nonbridge oxygens of GTP.

Why does the Ras switch "break" by oncogenic mutations?

The elucidation of the structure of the RasGAP complex provides what is perhaps the most detailed link between protein structure and cancer causing mutations. In particular, it is known that mutations of Gln 61 destroy the GTPase activity of the complex, locks the cell in its ON state and thus, can cause cancer. It is entirely unclear however, why this specific mutation is so important. The present work uncovers the elusive role of Gln 61 by computer simulation of the GTPase reaction in Ras, RasGAP and of their mutants. Simulations of the effects of mutations of Gln 61 reproduce the corresponding observed changes in activation energies and allow us to analyze the energy contributions to these effects. It is found that Gln 61 does not operate in a direct chemical way nor by a direct electrostatic or steric interaction with the transition state (TS). Instead, oncogenic mutations of Gln 61 lead to the destruction of the exquisitely preorganized catalytic configuration of the active site of the RasGAP complex. This "allosteric" effect causes a major reduction in the electrostatic stabilization of the TS. Our findings have general relevance to other proteins that control signal transduction processes.

Hydrolysis of cyclic phosphates by ribonuclease A: a computational study using a simplified ab initio quantum model

The second step in the enzyme-catalyzed hydrolysis of phosphate esters by ribonuclease A (RNase A) was studied using an ab initio quantum-based model of the active site including constrained parts of three critical residues, His-12, His-119, and Lys-41, and a small substrate. The competition between release of the cyclic phosphate intermediate and subsequent hydrolysis following transphosphorylation was explored to determine the electronic factors that contribute to preferential intermediate product release observed experimentally. The structural and energetic results obtained at both the RHF and MP2 levels reveal several contributing factors consistent with experimental observation. Although the intrinsic electronic effects tend to favor hydrolysis slightly with an overall activation free energy of approximately 70 kJ mol(-1), entropic and environmental effects favor release of the cyclic phosphate intermediate over hydrolysis. Exploration of the second, hydrolysis step also revealed interesting similarity with the transphosphorylation step, including the observation of autocatalysis by the substrate. Moreover, both steps of the overall RNase A reaction reveal multiple pathways involving proton transfers to sites of similar proton affinities. The anionic phosphate in both steps can act as a stable proton binding site as protons are moved around the active site throughout the progress of the reaction. These results suggest autocatalysis may be representative of more general behavior in enzymes containing highly charged substrates, especially phosphates.

Essential role of histidine 84 in elongation factor Tu for the chemical step of GTP hydrolysis on the ribosome

Elongation factor Tu (EF-Tu) is a GTP-binding protein that delivers aminoacyl-tRNA to the A site of the ribosome during protein synthesis. The mechanism of GTP hydrolysis in EF-Tu on the ribosome is poorly understood. It is known that mutations of a conserved histidine residue in the switch II region of the factor, His84 in Escherichia coli EF-Tu, impair GTP hydrolysis. However, the partial reaction which is directly affected by mutations of His84 was not identified and the effect on GTP hydrolysis was not quantified. Here, we show that the replacement of His84 with Ala reduces the rate constant of GTP hydrolysis more than 10(6)-fold, whereas the preceding steps of ternary complex binding to the ribosome, codon recognition and, most importantly, the GTPase activation step are affected only slightly. These results show that His84 plays a key role in the chemical step of GTP hydrolysis. Rate constants of GTP hydrolysis by wild-type EF-Tu, measured using the slowly hydrolyzable GTP analog, GTPgammaS, showed no dependence on pH, indicating that His84 does not act as a general base. We propose that the catalytic role of His84 is to stabilize the transition state of GTP hydrolysis by hydrogen bonding to the attacking water molecule or, possibly, the gamma-phosphate group of GTP.

GTPase catalysis by Ras and other G-proteins: insights from Substrate Directed SuperImposition

Comparisons of different protein structures are commonly carried out by superimposing the coordinates of the protein backbones or selected parts of the proteins. When the objective is analysis of similarities and differences in the enzyme's active site, there is an inherent problem in using the same domains for the superimposition. In this work we use a comparative approach termed here "Substrate Directed SuperImposition" (SDSI). It entails the superimposition of multiple protein-substrate structures using exclusively the coordinates of the comparable substrates. SDSI has the advantage of unbiased comparison of the active-site environment from the substrate's point of view. Our analysis extends previous usage of similar approaches to comparison of enzyme catalytic machineries. We applied SDSI to various G-protein structures for dissecting the mechanism of the GTPase reaction that controls the signaling activity of this important family. SDSI indicates that dissimilar G-proteins stabilize the transition state of the GTPase reaction similarly and supports the commonality of the critical step in this reaction, the reorientation of the critical arginine and glutamine. Additionally, we ascribe the catalytic inefficiency of the small G-protein Ras to the great flexibility of its active site and downplay the possible catalytic roles of the Lys16 residue in Ras GTPase. SDSI demonstrated that in contrast to all other Gly12 Ras mutants, which are oncogenic, the Gly12-->Pro mutant does not interfere with the catalytic orientation of the critical glutamine. This suggests why this mutant has a higher rate of GTP hydrolysis and is non-transforming. Remarkably, SDSI also revealed similarities in the divergent catalytic machineries of G-proteins and UMP/CMP kinase. Taken together, our results promote the use of SDSI to compare the catalytic machineries of both similar and different classes of enzymes.

Significance of GTP hydrolysis in Ypt1p-regulated endoplasmic reticulum to Golgi transport revealed by the analysis of two novel Ypt1-GAPs

We here report on the identification and detailed biochemical characterization of two novel GTPase-activating proteins, Gyp5p and Gyp8p, whose efficient substrate is Ypt1p, a Ypt/Rab-GTPase essential for endoplasmic reticulum-to-Golgi trafficking in yeast. Gyp5p accelerated the intrinsic GTPase activity of Ypt1p 4.2 x 10(4)-fold and, surprisingly, the 40-fold reduced GTP hydrolysis rate of Ypt1(Q67L)p 1.5 x 10(4)-fold. At steady state, the two newly discovered GTPase-activating proteins (GAPs) as well as the previously described Gyp1p, which also uses Ypt1p as the preferred substrate, display different subcellular localization. To add to an understanding of the significance of Ypt1p-bound GTP hydrolysis in vivo, yeast strains expressing the GTPase-deficient Ypt1(Q67L)p and having different Ypt1-GAP genes deleted were created. Depending on the genetic background, different mutants exhibited growth defects at low temperature and, already at permissive temperature, various morphological alterations resembling autophagy. Transport of proteins was not significantly impaired. Growth defects of Ypt1(Q67L)-expressing cells could be suppressed on high expression of all three Ypt1-GAPs. We propose that permanently active Ypt1p leads to increased vesicle fusion, which might induce previously unnoticed autophagic degradation of exaggerated membrane-enclosed structures. The data indicate that hydrolysis of Ypt1p-bound GTP is a prerequisite for a balanced vesicle flow between endoplasmic reticulum and Golgi compartments.

Conformational changes during the catalytic cycle of gluconate kinase as revealed by X-ray crystallography

The crystal structure of gluconate kinase from Escherichia coli has been determined to 2.0 A resolution by X-ray crystallography. The three-dimensional structure was solved by multi-wavelength anomalous dispersion, using a crystal of selenomethionine-substituted enzyme. Gluconate kinase is an alpha/beta structure consisting of a twisted parallel beta-sheet surrounded by alpha-helices with overall topology similar to nucleoside monophosphate (NMP) kinases, such as adenylate kinase. In order to identify residues involved in substrate binding and catalysis, structures of binary complexes with ATP, the ATP analogue adenosine 5'-(beta,gamma-methylene) triphosphate and the product, gluconate-6-phosphate have been determined. Significant conformational changes are induced upon binding of ATP to the enzyme. The largest changes involve a hinge-bending motion of the NMP(bind) part and a motion of the LID with adjacent helices, which opens the cavity to the second substrate, gluconate. Opening of the active site cleft upon ATP binding is the opposite of what has been observed in the NMP kinase family so far, which usually close their active site to prevent fortuitous hydrolysis of ATP. The conformational change positions the side-chain of Arg120 to stack with the purine ring of ATP and the side-chain of Arg124 is shifted to interact with the alpha-phosphate in ATP, at the same time protecting ATP from solvent water. The beta and gamma-phosphate groups of ATP bind in the predicted P-loop. A conserved lysine side-chain interacts with the gamma-phosphate group, and might promote phosphoryl transfer. Gluconate-6-phosphate binds with its phosphate group in a similar position as the gamma-phosphate of ATP, consistent with inline phosphoryl transfer. The gluconate binding-pocket in GntK is located in a different position than the nucleoside binding-site usually found in NMP kinases.

Revisiting the structural flexibility of the complex p21(ras)-GTP: the catalytic conformation of the molecular switch II

The hydrolysis of GTP in p21(ras) triggers conformational changes that regulate the ras/ERK signaling pathway. An important active site residue is Gln61, which has been found to be mutated in 30% of human tumors. The dynamics of the active site conformation is studied by using molecular dynamics simulation of two independent structures of the GTP-bound uncomplexed enzyme. Two distinct conformations of the enzyme are observed, in which the side-chain residue Gln61 is in different orientations. Essential dynamics analysis is used to describe the essential motions in the transition between the two conformations. Results are compared with earlier simulations of p21(ras) and its complex with GTPase activating protein p21-GAP.

gamma-phosphate protonation and pH-dependent unfolding of the Ras.GTP.Mg2+ complex: a vibrational spectroscopy study

The interdependence of GTP hydrolysis and the second messenger functions of virtually all GTPases has stimulated intensive study of the chemical mechanism of the hydrolysis. Despite numerous mutagenesis studies, the presumed general base, whose role is to activate hydrolysis by abstracting a proton from the nucleophilic water, has not been identified. Recent theoretical and experimental work suggest that the gamma-phosphate of GTP could be the general base. The current study investigates this possibility by studying the pH dependence of the vibrational spectrum of the Ras.GTP.Mg(2+) and Ras.GDP.Mg(2+) complexes. Isotope-edited IR studies of the Ras.GTP.Mg(2+) complex show that GTP remains bound to Ras at pH as low as 2.0 and that the gamma-phosphate is not protonated at pH > or = 3.3, indicating that the active site decreases the gamma-phosphate pK(a) by at least 1.1 pK(a) units compared with solution. Amide I studies show that the Ras.GTP.Mg(2+) and Ras.GDP.Mg(2+) complexes partially unfold in what appear to be two transitions. The first occurs in the pH range 5.4-2.6 and is readily reversible. Differences in the pH-unfolding midpoints for the Ras.GTP.Mg(2+) and Ras.GDP.Mg(2+) complexes (3.7 and 4.8, respectively) reveal that the enzyme-gamma-phosphoryl interactions stabilize the structure. The second transition, pH 2.6-1.7, is not readily reversed. The pH-dependent unfolding of the Ras.GTP.Mg(2+) complex provides an alternative interpretation of the data that had been used to support the gamma-phosphate mechanism, thereby raising the issue of whether this mechanism is operative in GTPase-catalyzed GTP hydrolysis reactions.

Substrate assisted catalysis -- application to G proteins

The idea that both the substrate and the enzyme contribute to catalysis (substrate assisted catalysis; SAC) is applicable to guanine nucleotide-binding proteins (G proteins). Naturally occurring SAC uses GTP as a general base in the GTPase reaction catalyzed by G proteins. Engineered SAC has identified a putative rate-limiting step for the GTPase reaction and shown that GTPase-deficient oncogenic Ras mutants are not irreversibly impaired. Thus, anti-cancer drugs could potentially be designed to restore the blocked GTPase reaction.

Phosphoryl transfer by a concerted reaction mechanism in UMP/CMP-kinase

The reaction mechanism of phosphoryl transfer catalyzed by UMP/CMP-kinase from Dictyostelium discoideum was investigated by semiempirical AM1 molecular orbital computations of an active site model system derived from crystal structures that contain a transition state analog or a bisubstrate inhibitor. The computational results suggest that the nucleoside monophosphate must be protonated for the forward reaction while it is unprotonated in the presence of aluminium fluoride, a popular transition state analog for phosphoryl transfer reactions. Furthermore, a compactification of the active site model system during the reaction and for the corresponding complex containing AlF3 was observed. For the active site residues that are part of the LID domain, conformational flexibility during the reaction proved to be crucial. On the basis of the calculations, a concerted phosphoryl transfer mechanism is suggested that involves the synchronous shift of a proton from the monophosphate to the transferred PO3-group. The proposed mechanism is thus analogous to the phosphoryl transfer mechanism in cAMP-dependent protein kinase that phosphorylates the hydroxyl groups of serine residues.

Application of linear free energy relationships to the serpin-proteinase inhibition mechanism

Linear free energy relationships can be used to link the changes in rate constant for a reaction to changes in the equilibrium caused by alterations in structure. While they have most often been used in the analysis of chemical reactions, they have also been employed to resolve questions in enzymology and protein folding. Here we analyze the reaction of a serpin with a panel of six serine proteinases, and observe that a linear free energy relationship exists between the true second-order rate constant for reaction, k(inh), and the inhibition constant, K(I), indicating that formation of the covalent serpin-enzyme complex may be reversible.