The Mechanism of Guanosine Nucleotide Hydrolysis by p21 c-Ha-ras

The Ras switch in structural and historical perspective

Since its discovery as an oncogene more than 40 years ago, Ras has been and still is in the focus of many academic and pharmaceutical labs around the world. A huge amount of work has accumulated on its biology. However, many questions about the role of the different Ras isoforms in health and disease still exist and a full understanding will require more intensive work in the future. Here we try to survey some of the structural findings in a historical perspective and how it has influenced our understanding of structure-function and mechanistic relationships of Ras and its interactions. The structures show that Ras is a stable molecular machine that uses the dynamics of its switch regions for the interaction with all regulators and effectors. This conformational flexibility has been used to create small molecule drug candidates against this important oncoprotein.

Mechanism of the intrinsic arginine finger in heterotrimeric G proteins

Heterotrimeric G proteins are crucial molecular switches that maintain a large number of physiological processes in cells. The signal is encoded into surface alterations of the Gα subunit that carries GTP in its active state and GDP in its inactive state. The ability of the Gα subunit to hydrolyze GTP is essential for signal termination. Regulator of G protein signaling (RGS) proteins accelerates this process. A key player in this catalyzed reaction is an arginine residue, Arg178 in Gα_i1, which is already an intrinsic part of the catalytic center in Gα in contrast to small GTPases, at which the corresponding GTPase-activating protein (GAP) provides the arginine "finger." We applied time-resolved FTIR spectroscopy in combination with isotopic labeling and site-directed mutagenesis to reveal the molecular mechanism, especially of the role of Arg178 in the intrinsic Gα_i1 mechanism and the RGS4-catalyzed mechanism. Complementary biomolecular simulations (molecular mechanics with molecular dynamics and coupled quantum mechanics/molecular mechanics) were performed. Our findings show that Arg178 is bound to γ-GTP for the intrinsic Gα_i1 mechanism and pushed toward a bidentate α-γ-GTP coordination for the Gα_i1·RGS4 mechanism. This movement induces a charge shift toward β-GTP, increases the planarity of γ-GTP, and thereby catalyzes the hydrolysis.

RGK family G-domain:GTP analog complex structures and nucleotide-binding properties

The RGK family of small G-proteins, including Rad, Gem, Rem1, and Rem2, is inducibly expressed in various mammalian tissues and interacts with voltage-dependent calcium channels and Rho kinase. Many questions remain regarding their physiological roles and molecular mechanism. Previous crystallographic studies reported RGK G-domain:guanosine di-phosphate structures. To test whether RGK proteins undergo a nucleotide-induced conformational change, we determined the crystallographic structures of Rad:GppNHp and Rem2:GppNHp to 1.7 and 1.8 Å resolutions, respectively. Also, we characterized the nucleotide-binding properties and conformations for Gem, Rad, and several structure-based mutants using fluorescence spectroscopy. The results suggest that RGK G-proteins may not behave as Ras-like canonical nucleotide-induced molecular switches. Further, the RGK proteins have differing structures and nucleotide-binding properties, which may have implications for their varied action on effectors.

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.

Role of the arginine finger in Ras.RasGAP revealed by QM/MM calculations

In the Ras.Ras.GAP complex, hydrolysis of guanosine triphosphate is strongly accelerated GAP as compared to Ras alone. This is largely attributed to the arginine finger R789(GAP) pointing to AlF(x) in the transition state analogue. We performed QM/MM simulations where triphosphate was treated using the quantum mechanical method of density functional theory, while the protein complex and water environment were described classically using MD. Compared to Ras, the crucial electron shift, bond stretching and distortion towards an eclipsed gamma-to-beta orientation are much more pronounced. The arginine finger is shown to act by displacing water out of the binding niche. The resulting enhanced electrostatic field catalyses the cleavage step.

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.

Phosphate Recognition in Structural Biology

Drug-discovery research in the past decade has seen an increased selection of targets with phosphate recognition sites, such as protein kinases and phosphatases, in the past decade. This review attempts, with the help of database-mining tools, to give an overview of the most important principles in molecular recognition of phosphate groups by enzymes. A total of 3003 X-ray crystal structures from the RCSB Protein Data Bank with bound organophosphates has been analyzed individually, in particular for H-bonding interactions between proteins and ligands. The various known binding motifs for phosphate binding are reviewed, and similarities to phosphate complexation by synthetic receptors are highlighted. An analysis of the propensities of amino acids in various classes of phosphate-binding enzymes showed characteristic distributions of amino acids used for phosphate binding. This review demonstrates that structure-based lead development and optimization should carefully address the phosphate-binding-site environment and also proposes new alternatives for filling such sites.

QM/MM modeling the Ras-GAP catalyzed hydrolysis of guanosine triphosphate

The mechanism of the hydrolysis reaction of guanosine triphosphate (GTP) by the protein complex Ras-GAP (p21(ras) - p120(GAP)) has been modeled by the quantum mechanical-molecular mechanical (QM/MM) and ab initio quantum calculations. Initial geometry configurations have been prompted by atomic coordinates of a structural analog (PDBID:1WQ1). It is shown that the minimum energy reaction path is consistent with an assumption of two-step chemical transformations. At the first stage, a unified motion of Arg789 of GAP, Gln61, Thr35 of Ras, and the lytic water molecule 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 transition state TS1. At the second stage, Gln61 abstracts and releases protons within the subsystem including Gln61, the lytic water molecule and the gamma-phosphate group of GTP through the corresponding transition state TS2. Direct quantum calculations show that, in this particular environment, the reaction GTP + H(2)O --> GDP + H(2)PO(4) (-) can proceed with reasonable activation barriers of less than 15 kcal/mol at every stage. This conclusion leads to a better understanding of the anticatalytic effect of cancer-causing mutations of Ras, which has been debated in recent years.

Phosphoryl transfer in Ras proteins, conclusive or elusive?

The chemical mechanism of GTP hydrolysis by GTP-binding proteins of the Ras superfamily continues to inspire both experimental and computational biologists. The debate centres on the nature of the transition state, with arguments for both dissociative and associative, and whether there is a common GTPase mechanism for these proteins. In a recent structural analysis of Rab11, the product P(i) was found in an unusual configuration. This finding indicates that substrate-assisted catalysis might operate as a mechanism to enable nucleophilic attack in the intrinsic GTPase reaction, and would thus favour a pentavalent phosphorane intermediate. Recent findings on the GAP-mediated reaction of different Ras proteins suggest that a common mechanism might not exist and that G proteins probably show a continuum of electronic configurations in the transition state.

Fourier Transform Infrared Spectroscopy on the Rap·RapGAP Reaction, GTPase Activation without an Arginine Finger

GTPase activating proteins (GAPs) down-regulate Ras-like proteins by stimulating their GTP hydrolysis, and a malfunction of this reaction leads to disease formation. In most cases, the molecular mechanism of activation involves stabilization of a catalytic Gln and insertion of a catalytic Arg into the active site by GAP. Rap1 neither possesses a Gln nor does its cognate Rap-GAP employ an Arg. Recently it was proposed that RapGAP provides a catalytic Asn, which substitutes for the Gln found in all other Ras-like proteins (Daumke, O., Weyand, M., Chakrabarti, P. P., Vetter, I. R., and Wittinghofer, A. (2004) Nature 429, 197-201). Here, RapGAP-mediated activation has been investigated by time-resolved Fourier transform infrared spectroscopy. Although the intrinsic hydrolysis reactions of Rap and Ras are very similar, the GAP-catalyzed reaction shows unique features. RapGAP binding induces a GTP(*) conformation in which the three phosphate groups are oriented such that they are vibrationally coupled to each other, in contrast to what was seen in the intrinsic and the Ras.RasGAP reactions. However, the charge shift toward beta-phosphate observed with RasGAP was also observed for RapGAP. A GDP.P(i) intermediate accumulates in the GAP-catalyzed reaction, because the release of P(i) is eight times slower than the cleavage reaction, and significant GTP synthesis from GDP.P(i) was observed. Partial steps of the cleavage reaction are correlated with structural changes of protein side groups and backbone. Thus, the Rap.RapGAP catalytic machinery compensates for the absence of a cis-Gln by a trans-Asn and for the catalytic Arg by inducing a different GTP conformation that is more prone to be attacked by a water molecule.

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.

The mechanism of Ras GTPase activation by neurofibromin

Individual rate constants have been determined for each step of the Ras.GTP hydrolysis mechanism, activated by neurofibromin. Fluorescence intensity and anisotropy stopped-flow measurements used the fluorescent GTP analogue, mantGTP (2'(3')-O-(N-methylanthraniloyl)GTP), to determine rate constants for binding and release of neurofibromin. Quenched flow measurements provided the kinetics of the hydrolytic cleavage step. The fluorescent phosphate sensor, MDCC-PBP was used to measure phosphate release kinetics. Phosphate-water oxygen exchange, using (18)O-substituted GTP and inorganic phosphate (P(i)), was used to determine the extent of reversal of the hydrolysis step and of P(i) binding. The data show that neurofibromin and P(i) dissociate from the NF1.Ras.GDP.P(i) complex with identical kinetics, which are 3-fold slower than the preceding cleavage step. A model is presented in which the P(i) release is associated with the change of Ras from "GTP" to "GDP" conformation. In this model, the conformation change on P(i) release causes the large change in affinity of neurofibromin, which then dissociates rapidly.

IQGAP1 is a component of Cdc42 signaling to the cytoskeleton

The Ras-GAP related protein IQGAP1 binds several proteins, including actin, calmodulin, E-cadherin and the Rho family GTPase Cdc42. To gain insight into its in vivo function, IQGAP1 was overexpressed in mammalian cells. Transfection of IQGAP1 significantly increased the levels of active, GTP-bound Cdc42, resulting in the formation of peripheral actin microspikes. By contrast, transfection of an IQGAP1 mutant lacking part of the GAP-related domain (IQGAP1deltaGRD) substantially decreased the amount of GTP-bound Cdc42 in cell lysates. Consistent with these findings, IQGAP1DeltaGRD blocked Cdc42 function in cells that stably overexpress constitutively active Cdc42 and abrogated the effect of bradykinin on Cdc42. In cells transfected with IQGAP1deltaGRD, bradykinin was unable to activate Cdc42, translocate Cdc42 to the membrane fraction, or induce filopodia production. IQGAP1deltaGRD transfection altered cellular morphology, producing small, round cells that closely resemble Cdc42-/- cells. Some insight into the mechanism was provided by in vitro analysis, which revealed that IQGAP1deltaGRD increased the intrinsic GTPase activity of Cdc42, thereby increasing the amount of inactive, GDP-bound Cdc42. These data imply that IQGAP1 has a crucial role in transducing Cdc42 signaling to the cytoskeleton.

Ras-mediated cleavage of a GTP analogue by a novel mechanism

The small guanosine triphosphate (GTP) binding protein Ras is involved in many cellular signal transduction processes leading to cell growth, differentiation and apoptosis. Mutations in ras genes are found in a large number of human tumours. GTP hydrolysis, the process that normally leads to the transition of the Ras protein from the active (GTP-bound) form to the inactive (GDP-bound) form is impaired due to these oncogenic mutations. In contrast, the GTP analogue 3,4-diaminobenzophenone(DABP)-phosphoramidate-GTP, a substrate for GTP-binding proteins, enables switching to the inactive GDP form in both wild-type and oncogenic Ras. Here we show by HPLC, mass spectrometry and NMR spectroscopy that the mechanism of this DABP-GTPase reaction is different from the physiological GTPase reaction. The gamma-phosphate group is not attacked by a nucleophilic water molecule, but rather by the aromatic amino group of the analogue, which leads to the generation of a stable cyclic diamidate product. These findings have potential implications for the development of anti-Ras drugs.

Monitoring the GAP catalyzed H-Ras GTPase reaction at atomic resolution in real time

The molecular reaction mechanism of the GTPase-activating protein (GAP)-catalyzed GTP hydrolysis by Ras was investigated by time resolved Fourier transform infrared (FTIR) difference spectroscopy using caged GTP (P(3)-1-(2-nitro)phenylethyl guanosine 5'-O-triphosphate) as photolabile trigger. This approach provides the complete GTPase reaction pathway with time resolution of milliseconds at the atomic level. Up to now, one structural model of the GAP x Ras x GDP x AlF(x) transition state analog is known, which represents a "snap shot" along the reaction-pathway. As now revealed, binding of GAP to Ras x GTP shifts negative charge from the gamma to beta phosphate. Such a shift was already identified by FTIR in GTP because of Ras binding and is now shown to be enhanced by GAP binding. Because the charge distribution of the GAP x Ras x GTP complex thus resembles a more dissociative-like transition state and is more like that in GDP, the activation free energy is reduced. An intermediate is observed on the reaction pathway that appears when the bond between beta and gamma phosphate is cleaved. In the intermediate, the released P(i) is strongly bound to the protein and surprisingly shows bands typical of those seen for phosphorylated enzyme intermediates. All these results provide a mechanistic picture that is different from the intrinsic GTPase reaction of Ras. FTIR analysis reveals the release of P(i) from the protein complex as the rate-limiting step for the GAP-catalyzed reaction. The approach presented allows the study not only of single proteins but of protein-protein interactions without intrinsic chromophores, in the non-crystalline state, in real time at the atomic level.

The mechanism of GTP hydrolysis by Ras probed by Fourier transform infrared spectroscopy

Time-resolved Fourier transform infrared spectroscopy (FTIR) in combination with photo-induced release of (18)O-labeled caged nucleotide has been employed to address mechanistic issues of GTP hydrolysis by Ras protein. Infrared spectroscopy of Ras complexes with nitrophenylethyl (NPE)-[alpha-(18)O(2)]GTP, NPE-[beta-(18)O(4)]GTP, or NPE-[gamma-(18)O(3)]GTP upon photolysis or during hydrolysis afforded a substantially improved mode assignment of phosphoryl group absorptions. Photolysis spectra of hydroxyphenylacyl-GTP and hydroxyphenylacyl-GDP bound to Ras and several mutants, Ras(Gly(12))-Mn(2+), Ras(Pro(12)), Ras(Ala(12)), and Ras(Val(12)), were obtained and yielded valuable information about structures of GTP or GDP bound to Ras mutants. IR spectra revealed stronger binding of GDP beta-PO(3)(2-) moiety by Ras mutants with higher activity, suggesting that the transition state is largely GDP-like. Analysis of the photolysis and hydrolysis FTIR spectra of the [beta-nonbridge-(18)O(2), alphabeta-bridge-(18)O]GTP isotopomer allowed us to probe for positional isotope exchange. Such a reaction might signal the existence of metaphosphate as a discrete intermediate, a key species for a dissociative mechanism. No positional isotope exchange was observed. Overall, our results support a concerted mechanism, but the transition state seems to have a considerable amount of dissociative character. This work demonstrates that time-resolved FTIR is highly suitable for monitoring positional isotope exchange and advantageous in many aspects over previously used methods, such as (31)P NMR and mass spectrometry.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Structure of Gialpha1.GppNHp, autoinhibition in a galpha protein-substrate complex

The structure of the G protein Gialpha1 complexed with the nonhydrolyzable GTP analog guanosine-5'-(betagamma-imino)triphosphate (GppNHp) has been determined at a resolution of 1.5 A. In the active site of Gialpha1. GppNHp, a water molecule is hydrogen bonded to the side chain of Glu43 and to an oxygen atom of the gamma-phosphate group. The side chain of the essential catalytic residue Gln204 assumes a conformation which is distinctly different from that observed in complexes with either guanosine 5'-O-3-thiotriphosphate or the transition state analog GDP.AlF4-. Hydrogen bonding and steric interactions position Gln204 such that it interacts with a presumptive nucleophilic water molecule, but cannot interact with the pentacoordinate transition state. Gln204 must be released from this auto-inhibited state to participate in catalysis. RGS proteins may accelerate the rate of GTP hydrolysis by G protein alpha subunits, in part, by inserting an amino acid side chain into the site occupied by Gln204, thereby destabilizing the auto-inhibited state of Galpha.

Raman difference studies of GDP and GTP binding to c-Harvey ras

The vibrational spectra of phosphate modes for GDP and GTP bound to the c-Harvey p21(ras) protein have been determined using 18O isotope edited Raman difference spectroscopy. A number of the phosphate stretch frequencies are changed upon GDP/GTP binding to ras, and the results are analyzed by ab initio calculations and through the use of empirical relationships that relate bond orders and bond lengths to vibrational frequencies. Bound GDP is found to be strongly stabilized by its interactions, mostly electrostatic, with the active site Mg2+. Bound GTP also interacts with the active site Mg2+ via its beta-phosphate group, as expected on the basis of crystallographic studies of bound GppNp. The angle between the nonbridging P&bondDot;O bonds of the gamma-phosphate of bound GTP increase by about 1-2 degrees compared to its solution value, thus bringing about a geometry that is closer to planar for these bonds as expected for the putative pentacoordinated transition state geometry of the phosphotransfer reaction. Modeling of the interactions at the nucleotide binding site suggests that the water molecule in-line with the P-O bond is positioned to bring about the change in bond angle. Moreover, a weak fifth bond (about 0.03 vu) appears to be formed between it and the gamma-phosphorus atom of bound GTP with a concomitant weakening of the O-P bond between the GDP leaving group and the gamma-phosphorus atom. Hence, an important role of the active site structure appears to be the strategic positioning of this in-line water. These structural results are consistent with a reaction pathway for GTP hydrolysis in ras of synchronous bond formation between the gamma-phosphorus of GTP and the attacking nucleophile and bond breaking between the gamma-phosphorus and the GDP leaving group.

Nucleotide hydrolysis-dependent conformational changes in p21(ras) as studied using ESR spectroscopy

We have employed ESR spectroscopy using guanine nucleotides that contain a spin label at the 2',3'-position of the ribose to investigate structural changes in the proto-oncogene product p21(ras) that are dependent on nucleotide hydrolysis. The three nucleotide analogs used were 2',3'-(2,2,5, 5-tetramethyl-3-pyrroline-1-oxyl-3-carboxylic acid ester (SL) GTP, SL-GDP, and the non-hydrolyzable analog SL-guanylylimidodiphosphate. SL-GTP was hydrolyzed by p21 with rates similar to those for GTP hydrolysis and appears to be an excellent substrate analog. The ESR spectra of SL-GTP and SL-GDP in complex with p21 differ significantly when acquired at 0 degrees C or 5 degrees C indicating different environments (conformations) of the protein-bound radicals depending on the phosphorylation state of the bound nucleotide. We calculated the rate constant for the conformational change as deduced from the changes in the corresponding ESR spectra upon incubation of the p21.SL-GTP complex at 25 degrees C and compared it to the rate constant of hydrolysis of SL-GTP at the same temperature. The rate constant deduced from the ESR method was similar to that determined by a high performance liquid chromatography technique. The data are in agreement with the idea that a conformational change during GTP hydrolysis by p21 occurs simultaneously with the actual hydrolysis step.

The interaction of Ras with GTPase-activating proteins

Ras plays a major role as a molecular switch in many signal transduction pathways which lead to cell growth and differentiation. The GTPase reaction of Ras is of central importance in the function of the switch since it terminates Ras-effector interactions. GTPase-activating proteins (GAPs) accelerate the very slow intrinsic hydrolysis reaction of the GTP-bound Ras by several orders of magnitude and thereby act as presumably negative regulators of Ras action. The GTP hydrolysis of oncogenic mutants of Ras remains unaltered. In this review we discuss recent biochemical and structural findings relating to the mechanism of GAP action, which strengthen the hypothesis that GAP accelerates the actual cleavage step by stabilizing the transition state of the phosphoryl transfer reaction.

Conformational and dynamic differences between N-ras P21 bound to GTPgammaS and to GMPPNP as studied by NMR

Heteronuclear-edited proton-detected NMR methods are used to study the nucleotide-dependent conformational changes between the GMPPNP form of human N-ras P21 as compared to GDP and GTPgammaS forms. Full-length N-ras P21 was also compared with protein truncated beyond residue 167, to search for interaction points between the more invariant part of the protein and the variable C-terminal section. In both cases, the reporter was the 15N-H 2D spectrum of aspartate amide groups labeled with 15N. Small truncation-induced changes were seen in the spectrum at the resonances of Asp-54, -108, and -109 which are not far from the C-terminal and, surprisingly, at Asp-57 which is more remote. The spectrum obtained for the GMPPNP-ligated form is similar to that of the GTPgammaS form, except that peaks of several residues are weak at low temperature, and strongly temperature-dependent in their intensity, and a new resonance appears at 15 degrees C and above. The observations are discussed in terms of a two-state model for the GMPPNP-ligated protein, previously proposed by Geyer et al. [(1996) Biochemistry 35, 10308-10320].

The crystal structure of human rac1, a member of the rho-family complexed with a GTP analogue

The crystal structure of human rac1, a member of the rho family of small G-proteins, complexed with the non-hydrolysable GTP analogue, guanosine-5'-(beta gamma-imino)triphosphate (GMPPNP), has been determined by X-ray analysis at a resolution of 1.38 A. Comparison with the structure of H-ras indicates that rac1 has an extra alpha-helical domain that is characteristic of the rho G proteins, and may be involved in the signalling pathway of this family.

G protein mechanisms: insights from structural analysis

This review is concerned with the structures and mechanisms of a superfamily of regulatory GTP hydrolases (G proteins). G proteins include Ras and its close homologs, translation elongation factors, and heterotrimeric G proteins. These proteins share a common structural core, exemplified by that of p21ras (Ras), and significant sequence identity, suggesting a common evolutionary origin. Three-dimensional structures of members of the G protein superfamily are considered in light of other biochemical findings about the function of these proteins. Relationships among G protein structures are discussed, and factors contributing to their low intrinsic rate of GTP hydrolysis are considered. Comparison of GTP- and GDP-bound conformations of G proteins reveals how specific contacts between the gamma-phosphate of GTP and the switch II region stabilize potential effector-binding sites and how GTP hydrolysis results in collapse (or reordering) of these surfaces. A GTPase-activating protein probably binds to and stabilizes the conformation of its cognate G protein that recognizes the transition state for hydrolysis, and may insert a catalytic residue into the G protein active site. Inhibitors of nucleotide release, such as the beta gamma subunit of a heterotrimeric G protein, bind selectively to and stabilize the GDP-bound state. Release factors, such as the translation elongation factor, Ts, also recognize the switch regions and destabilize the Mg(2+)-binding site, thereby promoting GDP release. G protein-coupled receptors are expected to operate by a somewhat different mechanism, given that the GDP-bound form of many G protein alpha subunits does not contain bound Mg2+.

Structure of the GDP-Pi complex of Gly203-->Ala gialpha1: a mimic of the ternary product complex of galpha-catalyzed GTP hydrolysis

BACKGROUND

G proteins play a vital role in transmembrane signalling events. In their inactive form G proteins exist as heterotrimers consisting of an alpha subunit, complexed with GDP and a dimer of betagamma subunits. Upon stimulation by receptors, G protein alpha subunits exchange GDP for GTP and dissociate from betagamma . Thus activated, alphasubunits stimulate or inhibit downstream effectors. The duration of the activated state corresponds to the single turnover rate of GTP hydrolysis, which is typically in the range of seconds. In Gialpha1, the Gly203-->Ala mutation reduces the affinity of the substrate for Mg2+, inhibits a key conformational step that occurs upon GTP binding and consequently inhibits the release of betagamma subunits from the GTP complex. The structure of the Gly203-->Ala mutant of Gialpha1 (G203AGialpha1) bound to the slowly hydrolyzing analog of GTP (GTPgammaS) has been determined in order to elucidate the structural changes that take place during hydrolysis.

RESULTS

We have determined the three dimensional structure of a Gly203-->Ala mutant of Gialpha1 at 2.6 A resolution. Although crystals were grown in the presence of GTPgammaS and Mg2+, the catalytic site contains a molecule of GDP and a phosphate ion, but no Mg2+. The phosphate ion is bound to a site near that occupied by the gamma-phosphate of GTPgammaS in the activated wild-type alpha subunit. A region of the protein, termed the Switch II helix, twists and bends to adopt a conformation that is radically different from that observed in other Gialpha1 subunit complexes.

CONCLUSIONS

Under the conditions of crystallization, the Gly203-->Ala mutation appears to stabilize a conformation that may be similar, although perhaps not identical, to the transient ternary product complex of Gialpha1-catalyzed GTP hydrolysis. The rearrangement of the Switch II helix avoids a potential steric conflict caused by the mutation. However, it appears that dissociation of the gamma-phosphate from the pentacoordinate intermediate also requires a conformational change in Switch II. Thus, a conformational rearrangement of the Switch II helix may be required in Galpha-catalyzed GTP hydrolysis.

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

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

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

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

Site-directed mutagenesis of Thermus thermophilus elongation factor Tu. Replacement of His85, Asp81 and Arg300

His85 in Thermus thermophilus elongation factor Tu (EF-Tu) was replaced by glutamine, leucine and glycine residues, leading to [H85Q]EF-Tu, [H85L] EF-Tu and [H85G]EF-Tu, respectively. Asp81 was replaced by alanine leading to [D81A]EF-Tu, and replacement of Arg300 provided [R300I]EF-Tu. Glycine in position 85 of domain I induces a protease-sensitive site in domain II and causes complete protein degradation in vivo. A similar effect was observed when Asp81 was replaced by alanine or Arg300 by isoleucine. Degradation is probably due to disturbed interactions between the domains of EF-Tu.GTP, inducing a protease-sensitive cleavage site in domain II. [H85Q]EF-Tu, which can be effectively overproduced in Escherichia coli, is slower in poly(U)-dependent poly(Phe) synthesis, has lower affinity to aminoacyl-tRNA but shows only a slightly reduced rate of intrinsic GTP hydrolysis compared to the native protein. The GTPase of this protein variant is not efficiently stimulated by aminoacyl-tRNA and ribosomes. The slow GTPase of [H85Q]EF-Tu increases the fidelity of translation as measured by leucine incorporation into poly(Phe) in in vitro poly(U)-dependent ribosomal translation. Replacement of His85 in T. thermophilus EF-Tu by leucine completely deactivates the GTPase activity but does not substantially influence the aminoacyl-tRNA binding. [H85L]EF-Tu is inactive in poly(U)-dependent poly(Phe)-synthesis. The rate of nucleotide dissociation is highest for [H85L]EF-Tu, followed by [H85Q]EF-Tu and native T. thermophilus EF-Tu. Mutation of His85, a residue which is not directly involved in the nucleotide binding, thus influences the interaction of EF-Tu domains, nucleotide binding and the efficiency and rate of GTPase activity.

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

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

Biochemical characterization of the essential GTP-binding protein Obg of Bacillus subtilis

An essential guanine nucleotide-binding protein, Obg, of Bacillus subtilis has been characterized with respect to its enzymatic activity for GTP. The protein was seen to hydrolyze GTP with a Km of 5.4 microM and a kcat of 0.0061 min-1 at 37 degrees C. GDP was a competitive inhibitor of this hydrolysis, with an inhibition constant of 1.7 microM at 37 degrees C. The dissociation constant for GDP from the Obg protein was 0.5 microM at 4 degrees C and was estimated to be 1.3 microM at 37 degrees C. Approximately 80% of the purified protein was capable of binding GDP. In addition to hydrolysis of GTP, Obg was seen to autophosphorylate with this substrate. Subsequent release of the covalent phosphate proceeds at too slow a rate to account for the overall rate of GTP hydrolysis, indicating that in vitro hydrolysis does not proceed via the observed phosphoamidate intermediate. It was speculated that the phosphorylated form of the enzyme may represent either a switched-on or a switched-off configuration, either of which may be normally induced by an effector molecule. This enzyme from a temperature-sensitive mutant of Obg did not show significantly altered GTPase activity at the nonpermissive temperature.

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

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

Separate GTP binding and GTPase activating domains of a G alpha subunit

Most members of the guanosine triphosphatase (GTPase) superfamily hydrolyze guanosine triphosphate (GTP) quite slowly unless stimulated by a GTPase activating protein or GAP. The alpha subunits (G alpha) of the heterotrimeric G proteins hydrolyze GTP much more rapidly and contain an approximately 120-residue insert not found in other GTPases. Interactions between a G alpha insert domain and a G alpha GTP-binding core domain, both expressed as recombinant proteins, show that the insert acts biochemically as a GAP. The results suggest a general mechanism for GAP-dependent hydrolysis of GTP by other GTPases.

An NMR comparison of the changes produced by different guanosine 5'-triphosphate analogs in wild-type and oncogenic mutant p21ras

We have used nuclear magnetic resonance spectroscopy to compare the conformational changes produced by replacement of bound GDP by the GTP analogs guanosine 5'-O-(3-thiotriphosphate) (GTP gamma S) and guanylyl (beta, gamma-imido)diphosphate (GMPPNP) in wild-type p21ras as well as the oncogenic mutant (G12D)p21ras. We have used isotope-edited nuclear magnetic resonance spectroscopy to observe the amide resonances of selectively [15N]glycine and [15N]isoleucine labeled p21ras-nucleotide complexes. We find that eight of the nine resonances that respond strongly to GTP gamma S and GMPPNP binding are the same but that the nature of the effect appears different. With GTP gamma S, seven new resonances replace the eight resonances specifically associated with GDP-p21ras, but in GMPPNP-p21ras only two resonances replace the GDP-specific resonances that are lost. The resonance of Gly 60 is clearly shown to be responsive to replacement of GDP by GMPPNP, in addition to glycines 10, 12, 13, 15, and 75 and isoleucines 36, 21, and one other, that were found to respond to GTP gamma S by Miller et al. [Miller, A.-F., Papastavros, M. Z., & Redfield, A.G. (1992) Biochemistry 31, 10208-10216). The two GMPPNP-specific resonances observed appear in positions similar to GTP gamma S-specific resonances, and the GTP gamma S-specific resonances, although not lost altogether, are weaker than the GDP-specific resonances they replace. Thus, the two GTP analogs have similar effects on the spectrum of p21ras, suggesting that the effects are due to features common to both analogs. We propose that active site resonance intensities are specifically attenuated when GTP analogs are bound because interactions with the gamma-phosphate of GTP analogs couple the flexible loops 2 and 4 to the rigid loop 1 of the active site. The conformational heterogeneity and dynamics of loops 2 and 4 would be constrained by loop 1 but also transmitted to it. Coupled conformational exchange on a common intermediate time scale could explain the simultaneous loss of resonances from all three loops in the active site. In our comparison of wild-type and (G12D) GDP-p21ras, we find that the resonance of Ile 36 is not visible in (G12D)p21ras. In (G12D)p21ras, replacement of GDP by GTP gamma S causes the resonances of glycines 10, 13, 15, 60, and 75 and isoleucine 21 and four others to shift from their GDP-specific positions. GTP gamma S-specific resonances are observed for all but two of these.(ABSTRACT TRUNCATED AT 400 WORDS)

Mapping the nucleotide-dependent conformational change of human N-ras p21 in solution by heteronuclear-edited proton-observed NMR methods

Heteronuclear-edited proton-detected NMR methods are used to study the nucleotide-dependent conformational change between GDP- and GTP gamma S-bound forms of human N-ras p21. Amide groups of Asp are used as sensitive probes. When GTP gamma S is substituted for GDP in cellular N-ras p21, the chemical shifts of resonances Asp-47, -126, -154, and Asn-172, as well as Gly-77 and -151, are not sensitive to nucleotide exchange, whereas Asp-30, -33, -38, -54, -57, -69, -92, -105, and -119 are affected. Distinct chemical shift changes of Asp-33, -38, and -69 indicate that substantial structure changes occur in the effector-binding region and the switch II region. Crystallographic studies of H-ras p21 have indicated that the conformational differences are confined to residues 32-38 and 60-76. Our observations indicate that the nucleotide-dependent structural transitions of the protein in solution may not be identical to those in the crystal. They suggest that the peptide beyond Glu-76 participates in a conformational switch, and possibly is involved in effector function. We propose that the region roughly from Asp-92 to -105, and the region of guanine base-binding motif(s), e.g., 116NKXD, are candidate sites recognized by either a GDP/GTP release factor or a GTPase-affected protein.

Phosphorylation site of NtrC, a protein phosphatase whose covalent intermediate activates transcription

The NtrC transcription factor is a member of a family of homologous prokaryotic regulatory proteins that participate in the transduction of extracellular and nutritional signals. It has been demonstrated that the phosphate group from a histidine residue of the phosphorylated NtrB protein autokinase is transferred to the NtrC protein. Phosphorylation of the NtrC protein is transient and activates its transcriptional enhancement activity. We have investigated the site of phosphorylation of the Salmonella typhimurium NtrC protein and find that it is an aspartate residue (Asp-54) that is found within a sequence conserved in all of the members of the family of regulatory proteins. We propose that this phosphorylation is an NtrC protein histidine phosphatase catalytic intermediate. This conclusion suggests that the NtrC family should be viewed not as kinase substrates but as enzymes that can catalyze the hydrolysis of their activated forms in a concentration-independent fashion. They are similar in this sense to eukaryotic signal-transducing GTPases.

The structure of Ras protein: a model for a universal molecular switch

X-ray crystallography has revealed the molecular architecture of the cellular and oncogenic forms of p21Ha-ras, the protein encoded by the human Ha-ras gene, in both its active (GTP-bound) and in its inactive (GDP-bound) forms. From comparison of these two structures, a mechanism is suggested for the GTPase hydrolysis reaction that triggers the conformational change necessary for signal transduction. The structures have also allowed identification of the structural consequences of point mutations and the way in which they interfere with the intrinsic GTPase activity of p21ras. The p21ras structure is similar to that of the G-domain of elongation factor Tu (EF-Tu) from Escherichia coli, suggesting that p21ras can serve as a good model for other guanine nucleotide binding proteins.

Crystal structures at 2.2 A resolution of the catalytic domains of normal ras protein and an oncogenic mutant complexed with GDP

The biological functions of ras proteins are controlled by the bound guanine nucleotide GDP or GTP. The GTP-bound conformation is biologically active, and is rapidly deactivated to the GDP-bound conformation through interaction with GAP (GTPase Activating Protein). Most transforming mutants of ras proteins have drastically reduced GTP hydrolysis rates even in the presence of GAP. The crystal structures of the GDP complexes of ras proteins at 2.2 A resolution reveal the detailed interaction between the ras proteins and the GDP molecule. All the currently known transforming mutation positions are clustered around the bound guanine nucleotide molecule. The presumed "effector" region and the GAP recognition region are both highly exposed. No significant structural differences were found between the GDP complexes of normal ras protein and the oncogenic mutant with valine at position 12, except the side-chain of the valine residue. However, comparison with GTP-analog complexes of ras proteins suggests that the valine side-chain may inhibit GTP hydrolysis in two possible ways: (1) interacting directly with the gamma-phosphate and altering its orientation or the conformation of protein residues around the phosphates; and/or (2) preventing either the departure of gamma-phosphate on GTP hydrolysis or the entrance of a nucleophilic group to attack the gamma-phosphate. The structural similarity between ras protein and the bacterial elongation factor Tu suggests that their common structural motif might be conserved for other guanine nucleotide binding proteins.

The GTPase superfamily: conserved structure and molecular mechanism

GTPases are conserved molecular switches, built according to a common structural design. Rapidly accruing knowledge of individual GTPases--crystal structures, biochemical properties, or results of molecular genetic experiments--support and generate hypotheses relating structure to function in other members of the diverse family of GTPases.