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

Real-time structural transitions are coupled to chemical steps in ATP hydrolysis by Eg5 kinesin

At the biochemical level, motor proteins are enzymatic molecules that function by converting chemical energy into mechanical motion. The key element for energy transduction and a major unresolved question common for all motor proteins is the coordination between the chemical and conformational steps in ATP hydrolysis. Here we show time-lapse monitoring of an in vitro ATP hydrolysis reaction by the motor domain of a human Kinesin-5 protein (Eg5) using difference Fourier transform infrared spectroscopy and UV photolysis of caged ATP. In this first continuous observation of a biological reaction coordinate from substrate to product, direct spectral markers for two catalytic events are measured: proton abstraction from nucleophilic water by the catalytic base and formation of the inorganic phosphate leaving group. Simultaneous examination of conformational switching in Eg5, in parallel with catalytic steps, shows structural transitions in solution consistent with published crystal structures of the prehydrolytic and ADP-bound states. In addition, we detect structural modifications in the Eg5 motor domain during bond cleavage between the beta- and gamma-phosphates of ATP. This conclusion challenges mechanochemical models for motor proteins that utilize only two stages of the catalytic cycle to drive force and motion.

Allosteric modulation of Ras positions Q61 for a direct role in catalysis

Ras and its effector Raf are key mediators of the Ras/Raf/MEK/ERK signal transduction pathway. Mutants of residue Q61 impair the GTPase activity of Ras and are found prominently in human cancers. Yet the mechanism through which Q61 contributes to catalysis has been elusive. It is thought to position the catalytic water molecule for nucleophilic attack on the gamma-phosphate of GTP. However, we previously solved the structure of Ras from crystals with symmetry of the space group R32 in which switch II is disordered and found that the catalytic water molecule is present. Here we present a structure of wild-type Ras with calcium acetate from the crystallization mother liquor bound at a site remote from the active site and likely near the membrane. This results in a shift in helix 3/loop 7 and a network of H-bonding interactions that propagates across the molecule, culminating in the ordering of switch II and placement of Q61 in the active site in a previously unobserved conformation. This structure suggests a direct catalytic role for Q61 where it interacts with a water molecule that bridges one of the gamma-phosphate oxygen atoms to the hydroxyl group of Y32 to stabilize the transition state of the hydrolysis reaction. We propose that Raf together with the binding of Ca(2+) and a negatively charged group mimicked in our structure by the acetate molecule induces the ordering of switch I and switch II to complete the active site of Ras.

ATP hydrolysis in Eg5 kinesin involves a catalytic two-water mechanism

Motor proteins couple steps in ATP binding and hydrolysis to conformational switching both in and remote from the active site. In our kinesin.AMPPPNP crystal structure, closure of the active site results in structural transformations appropriate for microtubule binding and organizes an orthosteric two-water cluster. We conclude that a proton is shared between the lytic water, positioned for gamma-phosphate attack, and a second water that serves as a general base. To our knowledge, this is the first experimental detection of the catalytic base for any ATPase. Deprotonation of the second water by switch residues likely triggers subsequent large scale structural rearrangements. Therefore, the catalytic base is responsible for initiating nucleophilic attack of ATP and for relaying the positive charge over long distances to initiate mechanotransduction. Coordination of switch movements via sequential proton transfer along paired water clusters may be universal for nucleotide triphosphatases with conserved active sites, such as myosins and G-proteins.

Lone pair ... pi interactions between water oxygens and aromatic residues: quantum chemical studies based on high-resolution protein structures and model compounds

The pi electron cloud of aromatic centers is known to be involved in several noncovalent interactions such as C--H...pi, O--H...pi, and pi...pi interactions in biomolecules. Lone-pair (lp) ... pi interactions have gained attention recently and their role in biomolecular structures is being recognized. In this article, we have carried out systematic analysis of high-resolution protein structures and identified more than 400 examples in which water oxygen atoms are in close contact (distance < 3.5 A) with the aromatic centers of aromatic residues. Three different methods were used to build hydrogen atoms and we used a consensus approach to find out potential candidates for lp...pi interactions between water oxygen and aromatic residues. Quantum mechanical calculations at MP2/6-311++G(d,p) level on model systems based on protein structures indicate that majority of the identified examples have energetically favorable interactions. The influence of water hydrogen atoms was investigated by sampling water orientations as a function of two parameters: distance from the aromatic center and the angle between the aromatic plane and the plane formed by the three water atoms. Intermolecular potential surfaces were constructed using six model compounds representing the four aromatic amino acids and 510 different water orientations for each model compound. Ab initio molecular orbital calculations at MP2/6-311++G(d,p) level show that the interaction energy is favorable even when hydrogen atoms are farthest from the aromatic plane while water oxygen is pointing toward the aromatic center. The strength of such interaction depends upon the distance of water hydrogen atoms from the aromatic substituents. Our calculations clearly show that the lp...pi interactions due to the close approach of water oxygen and aromatic center are influenced by the positions of water hydrogen atoms and the aromatic substituents.

Mechanism of the chemical step for the guanosine triphosphate (GTP) hydrolysis catalyzed by elongation factor Tu

Elongation factor Tu (EF-Tu), the protein responsible for delivering aminoacyl-tRNAs (aa-tRNAs) to ribosomal A site during translation, belongs to the group of guanosine-nucleotide (GTP/GDP) binding proteins. Its active 'on'-state corresponds to the GTP-bound form, while the inactive 'off'-state corresponds to the GDP-bound form. In this work we focus on the chemical step, GTP+H(2)O-->GDP+Pi, of the hydrolysis mechanism. We apply molecular modeling tools including molecular dynamics simulations and the combined quantum mechanical-molecular mechanical calculations for estimates of reaction energy profiles for two possible arrangements of switch II regions of EF-Tu. In the first case we presumably mimic binding of the ternary complex EF-Tu.GTP.aa-tRNA to the ribosome and allow the histidine (His85) side chain of the protein to approach the reaction active site. In the second case, corresponding to the GTP hydrolysis by EF-Tu alone, the side chain of His85 stays away from the active site, and the chemical reaction GTP+H(2)O-->GDP+Pi proceeds without participation of the histidine but through water molecules. In agreement with the experimental observations which distinguish rate constants for the fast chemical reaction in EF-Tu.GTP.aa-tRNA.ribosome and the slow spontaneous GTP hydrolysis in EF-Tu, we show that the activation energy barrier for the first scenario is considerably lower compared to that of the second case.

Semisynthesis of H-Ras with a glutamic acid methylester at position 61

The mechanism of the guanosine triphosphate (GTP) hydrolysis reaction of small G-proteins such as Ras is generally understood; however, some important molecular details are still missing. One example concerns the role of Gln61 in the catalysis of the GTP hydrolysis reaction. This amino acid is frequently mutated in oncogenic Ras leading to constitutively active variants of the protein. To elucidate the role of Gln61, subtle structural changes were introduced at this position by exchanging the natural occurring glutamine against a glutamic acid methyl ester (GluOme). Thereby the H-bond donor properties of this residue are changed and analysis of the GTP hydrolysis reaction can provide information on the function of the native carboxamide moiety. Using a semisynthetic approach, Ras(1-166)Gln61GluOMe was synthesized by sequential native chemical ligation of three unprotected peptide segments. Peptides Ras(1-50) and Ras(51-79)Gln61GluOMe were synthesized using Boc chemistry. The C-terminal peptide Ras(80-166) was expressed in E. coli. Initial tests of this semisynthetic strategy were performed by synthesis of the N- and C-terminally truncated protein variant Ras(39-101)Gln61GluOMe. The identified optimal reaction conditions were then applied to the synthesis of Ras(1-166)Gln61GluOMe. Refolding of the semisynthetic product in the presence of GTP was successful and revealed intrinsic GTPase activity of Ras(1-166)Gln61GluOMe.

Transformation efficiency of RasQ61 mutants linked to structural features of the switch regions in the presence of Raf

Transformation efficiencies of Ras mutants at residue 61 range over three orders of magnitude, but the in vitro GTPase activity decreases 10-fold for all mutants. We show that Raf impairs the GTPase activity of RasQ61L, suggesting that the Ras/Raf complex differentially modulates transformation. Our crystal structures show that, in transforming mutants, switch II takes part in a network of hydrophobic interactions burying the nucleotide and precatalytic water molecule. Our results suggest that Y32 and a water molecule bridging it to the gamma-phosphate in the wild-type structure play a role in GTP hydrolysis in lieu of the Arg finger in the absence of GAP. The bridging water molecule is absent in the transforming mutants, contributing to the burying of the nucleotide. We propose a mechanism for intrinsic hydrolysis in Raf-bound Ras and elucidate structural features in the Q61 mutants that correlate with their potency to transform cells.

Structure of the ROC domain from the Parkinson's disease-associated leucine-rich repeat kinase 2 reveals a dimeric GTPase

Mutations in leucine-rich repeat kinase 2 (LRRK2) are the most common cause of Parkinson's disease (PD). LRRK2 contains a Ras of complex proteins (ROC) domain that may act as a GTPase to regulate its protein kinase activity. The structure of ROC and the mechanism(s) by which it regulates kinase activity are not known. Here, we report the crystal structure of the LRRK2 ROC domain in complex with GDP-Mg(2+) at 2.0-A resolution. The structure displays a dimeric fold generated by extensive domain-swapping, resulting in a pair of active sites constructed with essential functional groups contributed from both monomers. Two PD-associated pathogenic residues, R1441 and I1371, are located at the interface of two monomers and provide exquisite interactions to stabilize the ROC dimer. The structure demonstrates that loss of stabilizing forces in the ROC dimer is likely related to decreased GTPase activity resulting from mutations at these sites. Our data suggest that the ROC domain may regulate LRRK2 kinase activity as a dimer, possibly via the C-terminal of ROC (COR) domain as a molecular hinge. The structure of the LRRK2 ROC domain also represents a signature from a previously undescribed class of GTPases from complex proteins and results may provide a unique molecular target for therapeutics in PD.

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.

Surface change of Ras enabling effector binding monitored in real time at atomic resolution

Ras, the prototype of the Ras superfamily, acts as a molecular switch for cell growth. External growth signals induce a GDP-to-GTP exchange. This modifies the Ras surface (Ras(on)GTP) and enables effector binding, which then activates signal-transduction pathways. GTP hydrolysis, catalysed by Ras and GAP, returns the signal to "off" (Ras(off)GDP). Oncogenic mutations in Ras prevent this hydrolysis, and thereby cause uncontrolled cell growth. In the Ras(off)-to-Ras(on) transition, the Ras surface is changed by a movement of the switch I loop that controls effector binding. We monitored this surface change at atomic resolution in real time by time-resolved FTIR (trFTIR) spectroscopy. In the transition from Ras(off) to Ras(on) a GTP-bound intermediate is now identified, in which effector binding is still prevented (Ras(off)GTP). The loop movement from Ras(off)GTP to Ras(on)GTP was directly monitored by the C=O vibration of Thr35. The structural change creates a binding site with a rate constant of 5 s(-1) at 260 K. A small molecule that shifted the equilibrium from the Ras(on)GTP state towards the Ras(off)GTP state would prevent effector binding, even if hydrolysis were blocked by oncogenic mutations. We present a spectroscopic fingerprint of both states that can be used as an assay in drug screening for such small molecules.

Slow conformational dynamics of the guanine nucleotide-binding protein Ras complexed with the GTP analogue GTPgammaS

The guanine nucleotide-binding protein Ras occurs in solution in two different conformational states, state 1 and state 2 with an equilibrium constant K(12) of 2.0, when the GTP analogue guanosine-5'-(beta,gamma-imido)triphosphate or guanosine-5'-(beta,gamma-methyleno)triphosphate is bound to the active centre. State 2 is assumed to represent a strong binding state for effectors with a conformation similar to that found for Ras complexed to effectors. In the other state (state 1), the switch regions of Ras are most probably dynamically disordered. Ras variants that exist predominantly in state 1 show a drastically reduced affinity to effectors. In contrast, Ras(wt) bound to the GTP analogue guanosine-5'-O-(3-thiotriphosphate) (GTPgammaS) leads to (31)P NMR spectra that indicate the prevalence of only one conformational state with K(12) > 10. Titration with the Ras-binding domain of Raf-kinase (Raf-RBD) shows that this state corresponds to effector binding state 2. In the GTPgammaS complex of the effector loop mutants Ras(T35S) and Ras(T35A) two conformational states different to state 2 are detected, which interconvert over a millisecond time scale. Binding studies with Raf-RBD suggest that both mutants exist mainly in low-affinity states 1a and 1b. From line-shape analysis of the spectra measured at various temperatures an activation energy DeltaH(|) (1a1b) of 61 kJ.mol(-1) and an activation entropy DeltaS(|) (1a1b) of 65 J.K(-1).mol(-1) are derived. Isothermal titration calorimetry on Ras bound to the different GTP-analogues shows that the effective affinity K(A) for the Raf-RBD to Ras(T35S) is reduced by a factor of about 20 compared to the wild-type with the strongest reduction observed for the GTPgammaS complex.

Structure of the GMPPNP-stabilized NG domain complex of the SRP GTPases Ffh and FtsY

Ffh and FtsY are GTPase components of the signal recognition particle co-translational targeting complex that assemble during the SRP cycle to form a GTP-dependent and pseudo twofold symmetric heterodimer. Previously the SRP GTPase heterodimer has been stabilized and purified for crystallographic studies using both the non-hydrolysable GTP analog GMPPCP and the pseudo-transition state analog GDP:AlF4, revealing in both cases a buried nucleotide pair that bridges and forms a key element of the heterodimer interface. A complex of Ffh and FtsY from Thermus aquaticus formed in the presence of the analog GMPPNP could not be obtained, however. The origin of this failure was previously unclear, and it was thought to have arisen from either instability of the analog, or, alternatively, from differences in its interactions within the tightly conscribed composite active site chamber of the complex. Using insights gained from the previous structure determinations, we have now determined the structure of the SRP GTPase targeting heterodimer stabilized by the non-hydrolysable GTP analog GMPPNP. The structure demonstrates how the different GTP analogs are accommodated within the active site chamber despite slight differences in the geometry of the phosphate chain. It also reveals a K+ coordination site at the highly conserved DARGG loop at the N/G interdomain interface.

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.

Ras oncogenes and their downstream targets

RAS proteins are small GTPases, which serve as master regulators of a myriad of signaling cascades involved in highly diverse cellular processes. RAS oncogenes have been originally discovered as retroviral oncogenes, and ever since constitutively activating RAS mutations have been identified in human tumors, they are in the focus of intense research. In this review, we summarize the biochemical properties of RAS proteins, trace down the evolution of RAS signaling and present an overview of the spatio-temporal activation of major RAS isoforms. We further discuss RAS effector pathways, their role in normal and transformed cell physiology and summarize ongoing attempts to interfere with aberrant RAS signaling. Finally, we comment on the role of micro RNAs in modulating RAS expression, contribution of RAS to stem cell function and on high-throughput analyses of RAS signaling networks.

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.

Hyperquenching for protein cryocrystallography

When samples having volumes characteristic of protein crystals are plunge cooled in liquid nitrogen or propane, most cooling occurs in the cold gas layer above the liquid. By removing this cold gas layer, cooling rates for small samples and modest plunge velocities are increased to 1.5 × 10(4) K s(-1), with increases of a factor of 100 over current best practice possible with 10 μm samples. Glycerol concentrations required to eliminate water crystallization in protein-free aqueous mixtures drop from ∼28% w/v to as low as 6% w/v. These results will allow many crystals to go from crystallization tray to liquid cryogen to X-ray beam without cryoprotectants. By reducing or eliminating the need for cryoprotectants in growth solutions, they may also simplify the search for crystallization conditions and for optimal screens. The results presented here resolve many puzzles, such as why plunge cooling in liquid nitrogen or propane has, until now, not yielded significantly better diffraction quality than gas-stream cooling.

A phosphoryl transfer intermediate in the GTPase reaction of Ras in complex with its GTPase-activating protein

The hydrolysis of nucleoside triphosphates by enzymes is used as a regulation mechanism in key biological processes. Here, the GTP hydrolysis of the protein complex of Ras with its GTPase-activating protein is monitored at atomic resolution in a noncrystalline state by time-resolved FTIR spectroscopy. At 900 ms, after the attack of water at the gamma-phosphate, there appears a H2PO4- intermediate that is shown to be hydrogen-bonded in an eclipsed conformation to the beta-phosphate of GDP. The H2PO4- intermediate is in a position where it can either reform GTP or be released from the protein in 7 s in the rate-limiting step of the GTPase reaction. We propose that such an intermediate also occurs in other GTPases and ATPases.

Class-1 release factor eRF1 promotes GTP binding by class-2 release factor eRF3

In eukaryotes, termination of mRNA translation is triggered by the essential polypeptide chain release factors eRF1, recognizing all three stop codons, and eRF3, a member of the GTPase superfamily with a role that has remained opaque. We have studied the kinetic and thermodynamic parameters of the interactions between eRF3 and GTP, GDP and the non-hydrolysable GTP analogue GDPNP in the presence (K(D)(GDP)=1.3+/-0.2 muM, K(D)(GTP) approximately 200 muM and K(D)(GDPNP)>160 muM) as well as absence (K(D)(GDP)=1.9+/-0.3 muM, K(D)(GTP) 0.7+/-0.2 muM and K(D)(GDPNP) approximately 200 muM) of eRF1. From the present data we propose that (i) free eRF3 has a strong preference to bind GDP compared to GTP (ii) eRF3 in complex with eRF1 has much stronger affinity to GTP than free eRF3 (iii) eRF3 in complex with PABP has weak affinity to GTP (iv) eRF3 in complex with eRF1 does not have strong affinity to GDPNP, implying that GDPNP is a poor analogue of GTP for eRF3 binding.

The role of residue stability in transient protein-protein interactions involved in enzymatic phosphate hydrolysis. A computational study

Finding why protein-protein interactions (PPIs) are so specific can provide a valuable tool in a variety of fields. Statistical surveys of so-called transient complexes (like those relevant for signal transduction mechanisms) have shown a tendency of polar residues to participate in the interaction region. Following this scheme, residues in the unbound partners have to compete between interacting with water or interacting with other residues of the protein. On the other hand, several works have shown that the notion of active site electrostatic preorganization can be used to interpret the high efficiency in enzyme reactions. This preorganization can be related to the instability of the residues important for catalysis. In some enzymes, in addition, conformational changes upon binding to other proteins lead to an increase in the activity of the enzymatic partner. In this article the linear response approximation version of the semimacroscopic protein dipoles Langevin dipoles (PDLD/S-LRA) model is used to evaluate the stability of several residues in two phosphate hydrolysis enzymes upon complexation with their activating partners. In particular, the residues relevant for PPI and for phosphate hydrolysis in the CDK2/Cyclin A and Ras/GAP complexes are analyzed. We find that the evaluation of the stability of residues in these systems can be used to identify not only active site regions but it can also be used as a guide to locate "hot spots" for PPIs. We also show that conformational changes play a major role in positioning interfacing residues in a proper "energetic" orientation, ready to interact with the residues in the partner protein surface. Thus, we extend the preorganization theory to PPIs, extrapolating the results we obtained from the above-mentioned complexes to a more general case. We conclude that the correlation between stability of a residue in the surface and the likelihood that it participates in the interaction can be a general fact for transient PPIs.

A newly designed microspectrofluorometer for kinetic studies on protein crystals in combination with x-ray diffraction

We present a new design for a fluorescence microspectrophotometer for use in kinetic crystallography in combination with x-ray diffraction experiments. The FLUMIX device (Fluorescence spectroscopy to monitor intermediates in x-ray crystallography) is built for 0 degrees fluorescence detection, which has several advantages in comparison to a conventional fluorometer with 90 degrees design. Due to the reduced spatial requirements and the need for only one objective, the system is highly versatile, easy to handle, and can be used for many different applications. In combination with a conventional stereomicroscope, fluorescence measurements or reaction initiation can be performed directly in a hanging drop crystallization setup. The FLUMIX device can be combined with most x-ray sources, normally without the need of a specialized mechanical support. As a biological model system, we have used H-Ras p21 with an artificially introduced photo-labile GTP precursor (caged GTP) and a covalently attached fluorophore (IANBD amide). Using the FLUMIX system, detailed information about the state of photolyzed crystals of the modified H-Ras p21 (p21(mod)) could be obtained. Measurements in combination with a synchrotron beamline showed significant fluorescence changes in p21(mod) crystals even within a few seconds of x-ray exposure at 100 K.

Empirical rules facilitate the search for binding sites on protein surfaces

Computational surface screening of 3D protein structures is a valuable means of finding possible docking sites for substrates, effectors and similar molecules. It can be improved by considering properties of molecules which are known to bind to protein surfaces, and thus reflect the required properties of binding sites. In-depth studies are available on drugs and lead compounds as binding partners with statistically assured properties. Here we present a simple strategy for finding binding sites, which is based on the empirical rule-of-five by Lipinski et al. for oral drugs and the rule-of-three by Congreve et al. for leads. The fast automated search with the new C-code TRIDOCK yields a preliminary set of sites, thus facilitating further investigation by visual, comparative and quantitative work. Possible binding sites are tagged by pseudo-atoms added to the structure file for molecular graphical evaluation. Usually, the strategy yields not just a few single sites, but an accumulation of several sites in known substrate binding pockets. Clusters are also found at known or putative protein-protein docking interfaces. A comparison of the activated and inactivated form of the GTPase Ras reveals clear differences and identifies a niche, which is possibly a suitable new target for compounds that bind specifically to activated Ras.

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.

Zooming in on ATP hydrolysis in F1

We summarize our current view of the reaction mechanism in F(1)-ATPase as it has emerged from experiment, theory, and computational studies over the last several years. ATP catalysis in the catalytic binding pockets of F(1) takes place without the release of any significant free energy and is efficiently driven by the combined action of two water molecules utilizing a so-called protein-relay mechanism. The chemical reaction itself is controlled by the spatial position of a key arginine residue.

Conformational states of human H-Ras detected by high-field EPR, ENDOR, and 31P NMR spectroscopy

Ras is a central constituent of the intracellular signal transduction that switches between its inactive state with GDP bound and its active state with GTP bound. A number of different X-ray structures are available. Different magnetic resonance techniques were used to characterise the conformational states of the protein and are summarised here. 31P NMR spectroscopy was used as probe for the environment of the phosphate groups of the bound nucleotide. It shows that in liquid solution additional conformational states in the GDP as well as in the GTP forms coexist which are not detected by X-ray crystallography. Some of them can also be detected by solid-state NMR in the micro crystalline state. EPR and ENDOR spectroscopy were used to probe the environment of the divalent metal ion (Mg2+ was replaced by Mn2+) bound to the nucleotide in the protein. Here again different states could be observed. Substitution of normal water by 17O-enriched water allowed the determination of the number of water molecules in the first coordination sphere of the metal ion. In liquid solution, they indicate again the existence of different conformational states. At low temperatures in the frozen state ENDOR spectroscopy suggests that only one state exists for the GDP- and GTP-bound form of Ras, respectively.

Theoretical IR spectroscopy based on QM/MM calculations provides changes in charge distribution, bond lengths, and bond angles of the GTP ligand induced by the Ras-protein

The GTPase Ras p21 is a crucial switch in cellular signal transduction. Fourier transform infrared (FTIR) spectra of the substrate guanosine triphosphate (GTP) show remarkable changes when it binds to the enzyme. The reduced band widths indicate that the flexible GTP molecule is guided by the protein into a preferred conformation. The delocalized phosphate vibrations of unbound GTP become localized. The frequency shifts show an electron movement toward beta-phosphate, which probably contributes to catalysis by reducing the free activation energy. To quantify these qualitative observations we performed QM/MM molecular dynamics simulations of Ras.GTP and GTP in water. The triphosphate part of GTP was treated quantum mechanically using density functional theory (DFT). Vibrational spectra were calculated in harmonic approximation with an average deviation of 3% from the experimental frequencies. This provides a high confidence in the computational results as vibrational spectra are highly sensitive to conformation and charge distribution. As compared to GTP in water, Ras-bound GTP shows a shift of negative charge of approximately 0.2 e toward the beta-phosphate from gamma-phosphate and from alpha-phosphate due to the positive charge of the magnesium ion, to a lesser extent of Lys-16, and surprisingly without any effect of the P-loop backbone. Magnesium and Gly-13 twist and bend the gamma-O-beta bonds such that the crucial bond is stretched before cleaving.

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.

Structural flexibility of small GTPases. Can it explain their functional versatility?

Multiple interactions with many different partners are responsible for the amazing functional versatility of proteins, especially those participating in cellular regulation. The structural properties that could facilitate multiple interactions are examined for small GTPases. The role of cellular constraints, compartmentation and scaffolds on protein-protein interactions is considered.

Superoxide anion radical modulates the activity of Ras and Ras-related GTPases by a radical-based mechanism similar to that of nitric oxide

Ras GTPases cycle between inactive GDP-bound and active GTP-bound states to modulate a diverse array of processes involved in cellular growth control. The activity of Ras is up-regulated by cellular agents, including both protein (guanine nucleotide exchange factors) and redox-active agents (nitric oxide (NO) and superoxide anion radical (O2*). We have recently elucidated the mechanism by which NO promotes guanine nucleotide dissociation of redox-active NKCD motif-containing Ras and Ras-related GTPases. In this study, we show that guanine nucleotide dissociation is enhanced upon exposure of the redox-active GTPases, Ras and Rap1A, to O2* and provide evidence for the efficient guanine nucleotide reassociation in the presence of the radical quenching agent ascorbate to complete guanine nucleotide exchange. In vivo, guanine nucleotide reassociation is necessary to populate Ras in its biologically active GTP-bound form after the dissociation of GDP. We further show that treatment of the redox-active GTPases with O2* releases GDP in form of an unstable the oxygenated GDP adduct, putatively assigned as 5-oxo-GDP. 5-Oxo-GDP was not produced from either the C118S or the F28L Ras variants upon the treatment of O2*, supporting the involvement of residues Cys118 and Phe28 in O2*-mediated Ras guanine nucleotide dissociation. These results indicate that the mechanism of O2*-mediated Ras guanine nucleotide dissociation is similar to that of NO/O2-mediated Ras guanine nucleotide dissociation.

Mechanism of free radical nitric oxide-mediated Ras guanine nucleotide dissociation

Ras proteins cycle between GDP-bound and GTP-bound states to modulate a diverse array of cellular growth processes. In this study, we have elucidated a mechanism by which nitric oxide, in the presence of oxygen (NO/O2), regulates Ras activity. We show that treatment of Ras with NO/O2 causes conversion of Ras-bound GDP into a free 463.3 Da nucleotide-nitration product. Mass and UV/visible spectroscopic analyses suggest that this nitration product is 5-guanidino-4-nitroimidazole diphosphate (NIm-DP), a degradation product of 5-nitro-GDP. These results indicate that NO/O2 mediates Ras guanine nucleotide exchange (GNE) by conversion of Ras-bound GDP into an unstable 5-nitro-GDP. 5-Nitro-GDP can be produced by radical-based reaction of the GDP guanine base with nitrogen dioxide (*NO2). We also provide evidence that the Ras Phe28 side-chain plays a key role in the formation of a NO/O2-induced Ras 5-nitro-GDP product. We previously proposed a mechanism of NO/O2-mediated Ras GNE, in which *NO2, formed by the reaction of NO with O2, generates a Ras Cys118 thiyl radical (Ras-S118) intermediate. In the present study, we provide evidence for a radical-based mechanism of NO/O2-mediated Ras GNE. According to this mechanism, reaction of NO with O2 produces *NO2. *NO2 then reacts with Ras to produce Ras-S118, which withdraws an electron from the Ras-bound guanine nucleotide base to produce a guanine nucleotide diphosphate cation radical (G(+)-DP) via the Phe28 side-chain. G(+)-DP is subsequently converted to a neutral radical, and can react with another *NO2 to produce 5-nitro-GDP. This radical-based reaction process disrupts key binding interactions between Ras and the guanine base, resulting in release of GDP from Ras and its conversion to free 5-nitro-GDP. This mechanism is likely to be common to other NKCD motif-containing Ras superfamily GTPases, as NO/O2 also facilitates GNE on the redox-active Rap1A and Rab3A GTPases.

X-ray crystallography of protein-ligand interactions

Crystal structures of protein-ligand complexes provide a detailed view of their spatial arrangement and interactions. In the case of stable, unreactive ligands, such as inhibitors or allosteric regulators, the complexes can be generated by cocrystallization or by soaking the ligand into fully grown crystals. In order to obtain highly occupied stochiometric complexes, the concentration and amount of ligand used needs to be considered. Protein complexes with reactive short-lived species that occur in chemical or binding reactions can be determined using monochromatic X-ray diffraction techniques via kinetic trapping approaches. To this end, the kinetics of the reaction has to be determined in the crystalline state and triggering methods to start the reaction need to be established. To facilitate data interpretation, the experimental conditions are usually chosen such that the peak concentration of the reactive species under investigation is maximized.

UreG, a chaperone in the urease assembly process, is an intrinsically unstructured GTPase that specifically binds Zn2+

Bacillus pasteurii UreG, a chaperone involved in the urease active site assembly, was overexpressed in Escherichia coli BL21(DE3) and purified to homogeneity. The identity of the recombinant protein was confirmed by SDS-PAGE, protein sequencing, and mass spectrometry. A combination of size exclusion chromatography and multiangle and dynamic laser light scattering established that BpUreG is present in solution as a dimer. Analysis of circular dichroism spectra indicated that the protein contains large portions of helices (15%) and strands (29%), whereas NMR spectroscopy indicated the presence of conformational fluxionality of the protein backbone in solution. BpUreG catalyzes the hydrolysis of GTP with a kcat=0.04 min(-1), confirming a role for this class of proteins in coupling energy requirements and nickel incorporation into the urease active site. BpUreG binds two Zn2+ ions per dimer, with a KD=42 +/- 3 microm, and has a 10-fold lower affinity for Ni2+. A structural model for BpUreG was calculated by using threading algorithms. The protein, in the fully folded state, features the typical structural architecture of GTPases, with an open beta-barrel surrounded by alpha-helices and a P-loop at the N terminus. The protein dynamic behavior observed in solution is critically discussed relative to the structural model, using algorithms for disorder predictions. The results suggest that UreG proteins belong to the class of intrinsically unstructured proteins that need the interaction with cofactors or other protein partners to perform their function. It is also proposed that metal ions such as Zn2+ could have important structural roles in the urease activation process.

Rapid evolution in conformational space: a study of loop regions in a ubiquitous GTP binding domain

The rapidly evolving subsets of a protein are often evident in multiple sequence alignments as poorly defined, gap-containing regions. We investigated the 3D context of these regions observed in 28 protein structures containing a GTP-binding domain assumed to be homologous to the transforming factor p21-RAS. The phylogenetic depth of this data set is such that it is possible to observe lineages sharing a common protein core that diverged early in the eukaryotic cell history. The sequence variability among these homolog proteins is directly linked to the structural variability of surface loops. We demonstrate that these regions are self-contained and thus mostly free of the evolutionary constraints imposed by the conserved core of the domain. These intraloop interactions have the property to create stem-like structures. Interestingly, these stem-like structures can be observed in loops of varying size, up to the size of small protein domains. We propose a model under which the diversity of protein topologies observed in these loops can be the product of a stochastic sampling of sequence and conformational space in a near-neutral fashion, while the proximity of the functional features of the domain core allows novel beneficial traits to be fixed. Our comparative observations, limited here to the proteins containing the RAS-like GTP-binding domain, suggest that a stochastic process of insertion/deletion analogous to "budding" of loops is a likely mechanism of structural innovation. Such a framework could be experimentally exploited to investigate the folding of increasingly complex model inserts.

ATP hydrolysis in the betaTP and betaDP catalytic sites of F1-ATPase

The enzyme F1-adenosine triphosphatase (ATPase) is a molecular motor that converts the chemical energy stored in the molecule adenosine triphosphate (ATP) into mechanical rotation of its gamma-subunit. During steady-state catalysis, the three catalytic sites of F1 operate in a cooperative fashion such that at every instant each site is in a different conformation corresponding to a different stage along the catalytic cycle. Notwithstanding a large amount of biochemical and, recently, structural data, we still lack an understanding of how ATP hydrolysis in F1 is coupled to mechanical motion and how the catalytic sites achieve cooperativity during rotatory catalysis. In this publication, we report combined quantum mechanical/molecular mechanical simulations of ATP hydrolysis in the betaTP and betaDP catalytic sites of F1-ATPase. Our simulations reveal a dramatic change in the reaction energetics from strongly endothermic in betaTP to approximately equienergetic in betaDP. The simulations identify the responsible protein residues, the arginine finger alphaR373 being the most important one. Similar to our earlier study of betaTP, we find a multicenter proton relay mechanism to be the energetically most favorable hydrolysis pathway. The results elucidate how cooperativity between catalytic sites might be achieved by this remarkable molecular motor.

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.

High resolution crystal structure of human Rab9 GTPase: a novel antiviral drug target

Rab GTPases and their effectors facilitate vesicular transport by tethering donor vesicles to their respective target membranes. Rab9 mediates late endosome to trans-Golgi transport and has recently been found to be a key cellular component for human immunodeficiency virus-1, Ebola, Marburg, and measles virus replication, suggesting that it may be a novel target in the development of broad spectrum antiviral drugs. As part of our structure-based drug design program, we have determined the crystal structure of a C-terminally truncated human Rab9 (residues 1-177) to 1.25-A resolution. The overall structure shows a characteristic nucleotide binding fold consisting of a six-stranded beta-sheet surrounded by five alpha-helices with a tightly bound GDP molecule in the active site. Structure-based sequence alignment of Rab9 with other Rab proteins reveals that its active site consists of residues highly conserved in the Rab GTPase family, implying a common catalytic mechanism. However, Rab9 contains seven regions that are significantly different in conformation from other Rab proteins. Some of those regions coincide with putative effector-binding sites and switch I and switch II regions identified by structure/sequence alignments. The Rab9 structure at near atomic resolution provides an excellent model for structure-based antiviral drug design.

Analysis of the open and closed conformations of the GTP-binding protein YsxC from Bacillus subtilis

Genetic analysis has suggested that the product of the Bacillus subtilis ysxC gene is essential for survival of the microorganism and hence may represent a target for the development of a novel anti-infective agent. B.subtilis YsxC is a member of the translation factor related class of GTPases and its crystal structure has been determined in an apo form and in complex with GDP and GMPPNP/Mg2+. Analysis of these structures has allowed us to examine the conformational changes that occur during the process of nucleotide binding and GTP hydrolysis. These structural changes particularly affect parts of the switch I and switch II region of YsxC, which become ordered and disordered, respectively in the "closed" or "on" GTP-bound state and disordered and ordered, respectively, in the "open" or "off" GDP-bound conformation. Finally, the binding of the magnesium cation results in subtle shifts of residues in the G3 region, at the start of switch II, which serve to optimize the interaction with a key aspartic acid residue. The structural flexibility observed in YsxC is likely to contribute to the role of the protein, possibly allowing transduction of an essential intracellular signal, which may be mediated via interactions with a conserved patch of surface-exposed, basic residues that lies adjacent to the GTP-binding site.

Biomolecular cryocrystallography: structural changes during flash-cooling

To minimize radiation damage, crystal structures of biological macromolecules are usually determined after rapid cooling to cryogenic temperatures, some 150-200 K below the normal physiological range. The biological relevance of such structures relies on the assumption that flash-cooling is sufficiently fast to kinetically trap the macromolecule and associated solvent in a room-temperature equilibrium state. To test this assumption, we use a two-state model to calculate the structural changes expected during rapid cooling of a typical protein crystal. The analysis indicates that many degrees of freedom in a flash-cooled protein crystal are quenched at temperatures near 200 K, where local conformational and association equilibria may be strongly shifted toward low-enthalpy states. Such cryoartifacts should be most important for strongly solvent-coupled processes, such as hydration of nonpolar cavities and surface regions, conformational switching of solvent-exposed side chains, and weak ligand binding. The dynamic quenching that emerges from the model considered here can also rationalize the glass transition associated with the atomic fluctuations in the protein.

The significance of the free energy of hydrolysis of GTP for signal-transducing and regulatory GTPases

A large number of GTP/GDP binding proteins, which in general have intrinsic and/or stimulatable GTPase activity, have been identified in recent years and are involved in a wide range of cellular regulatory and signal transducing processes. A common property of these proteins is that they exist in what is generally described as an active form when GTP is bound and an inactive (resting) form when GDP is present. Thus, the intrinsic or stimulated GTPase activity of these 'enzymes' serves to turn off a signal or to terminate a regulated process. It has been suggested that these proteins, together with ATPases whose prime function is to convert the free energy of ATP hydrolysis into another form of energy or into energy-requiring chemical reactions should be grouped together under the heading of 'energyases'. In this article, this suggestion is examined from the point of view of identifying the role of the free-energy of hydrolysis of GTP in the signal-transducing or regulatory process of the GTPases. It is concluded that there is a qualitative difference between ATPases and classical GTPases, in the sense that a quantitative relationship between the free-energy of GTP hydrolysis and the appearance of this energy in a different form cannot be directly defined. The significance of the high free energy of hydrolysis is that it allows efficient transition from the active to the inactive state of GTPases in spite of the tendency of the strong interaction of the GTP-bound form with a partner molecule ('effector'), an essential feature of their mode of action, to stabilize the GTP-bound form.

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.

Observation of slow dynamic exchange processes in Ras protein crystals by 31P solid state NMR spectroscopy

The folding, structure and biological function of many proteins are inherently dynamic properties of the protein molecule. Often, the respective molecular processes are preserved upon protein crystallization, leading, in X-ray diffraction experiments, to a blurring of the electron density map and reducing the resolution of the derived structure. Nuclear magnetic resonance (NMR) is known to be an alternative method to study molecular structure and dynamics. We designed and built a probe for phosphorus solid state NMR that allows for the first time to study static properties as well as dynamic processes in single-crystals of a protein by NMR spectroscopy. The sensitivity achieved is sufficient to detect the NMR signal from individual phosphorus sites in a 0.3mm(3) size single-crystal of GTPase Ras bound to the nucleotide GppNHp, that is, the signal from approximately 10(15) phosphorus nuclei. The NMR spectra obtained are discussed in terms of the conformational variability of the active center of the Ras-nucleotide complex. We conclude that, in the crystal, the protein complex exists in three different conformations. Magic angle spinning (MAS) NMR spectra of a powder sample of Ras-GppNHp show a splitting of one of the phosphate resonances and thus confirm this conclusion. The MAS spectra provide, furthermore, evidence of a slow, temperature-dependent dynamic exchange process in the Ras protein crystal.

The structural basis for the transition from Ras-GTP to Ras-GDP

The conformational changes in Ras that accompany the hydrolysis of GTP are critical to its function as a molecular switch in signaling pathways. Understanding how GTP is hydrolyzed by revealing the sequence of intermediary structures in the reaction is essential for understanding Ras signaling. Until now, no structure of an intermediate in GTP hydrolysis has been experimentally determined for Ras alone. We have solved the crystal structure of the Ala-59 to Gly mutant of Ras, (RasA59G), bound to guanosine 5'-imidotriphosphate or GDP to 1.7-A resolution. In the guanosine 5'-imidotriphosphate-bound form, this mutant adopts a conformation that is intermediate between the GTP- and GDP-bound forms of wild-type Ras and that is similar to what has been predicted by molecular dynamics simulation [Ma, J. P. & Karplus, M. (1997) Proc. Natl. Acad. Sci. USA 94, 11905-11910]. This conformation is stabilized by direct and water-mediated interactions between the switch 1 and switch 2 regions and is characterized by an increase in the binding affinity for GTP. We propose that the structural changes promoted by the Ala-59 to Gly mutation exhibit a discrete conformational state assumed by wild-type Ras during GTP hydrolysis.

Protein-water interactions in a dynamic world

Protein-water interactions are key to biological function. They have an underlying dynamic component that pervades the functional roles associated both with particular systems and with the properties of proteins in general. This article focuses on the specific ways in which the dynamics of water are important to protein structure, motion and adaptability to changes in the protein environment.

Alkaline phosphatase revisited: hydrolysis of alkyl phosphates

Escherichia coli alkaline phosphatase (AP) is the prototypical two metal ion catalyst with two divalent zinc ions bound approximately 4 A apart in the active site. Studies spanning half a century have elucidated many structural and mechanistic features of this enzyme, rendering it an attractive model for investigating the potent catalytic power of bimetallic centers. Unfortunately, fundamental mechanistic features have been obscured by limitations with the standard assays. These assays generate concentrations of inorganic phosphate (P(i)) in excess of its inhibition constant (K(i) approximately 1 muM). This tight binding by P(i) has affected the majority of published kinetic constants. Furthermore, binding limits k(cat)/K(m) for reaction of p-nitrophenyl phosphate, the most commonly employed substrate. We describe a sensitive (32)P-based assay for hydrolysis of alkyl phosphates that avoids the complication of product inhibition. We have revisited basic mechanistic features of AP with these alkyl phosphate substrates. The results suggest that the chemical step for phosphorylation of the enzyme limits k(cat)/K(m). The pH-rate profile and additional results suggest that the serine nucleophile is active in its anionic form and has a pK(a) of < or = 5.5 in the free enzyme. An inactivating pK(a) of 8.0 is observed for binding of both substrates and inhibitors, and we suggest that this corresponds to ionization of a zinc-coordinated water molecule. Counter to previous suggestions, inorganic phosphate dianion appears to bind to the highly charged AP active site at least as strongly as the trianion. The dependence of k(cat)/K(m) on the pK(a) of the leaving group follows a Brønsted correlation with a slope of beta(lg) = -0.85 +/- 0.1, differing substantially from the previously reported value of -0.2 obtained from data with a less sensitive assay. This steep leaving group dependence is consistent with a largely dissociative transition state for AP-catalyzed hydrolysis of phosphate monoesters. The new (32)P-based assay employed herein will facilitate continued dissection of the AP reaction by providing a means to readily follow the chemical step for phosphorylation of the enzyme.

The conformation of bound GMPPNP suggests a mechanism for gating the active site of the SRP GTPase

BACKGROUND

The signal recognition particle (SRP) is a phylogenetically conserved ribonucleoprotein that mediates cotranslational targeting of secreted and membrane proteins to the membrane. Targeting is regulated by GTP binding and hydrolysis events that require direct interaction between structurally homologous "NG" GTPase domains of the SRP signal recognition subunit and its membrane-associated receptor, SR alpha. Structures of both the apo and GDP bound NG domains of the prokaryotic SRP54 homolog, Ffh, and the prokaryotic receptor homolog, FtsY, have been determined. The structural basis for the GTP-dependent interaction between the two proteins, however, remains unknown.

RESULTS

We report here two structures of the NG GTPase of Ffh from Thermus aquaticus bound to the nonhydrolyzable GTP analog GMPPNP. Both structures reveal an unexpected binding mode in which the beta-phosphate is kinked away from the binding site and magnesium is not bound. Binding of the GTP analog in the canonical conformation found in other GTPase structures is precluded by constriction of the phosphate binding P loop. The structural difference between the Ffh complex and other GTPases suggests a specific conformational change that must accompany movement of the nucleotide from an "inactive" to an "active" binding mode.

CONCLUSIONS

Conserved side chains of the GTPase sequence motifs unique to the SRP subfamily may function to gate formation of the active GTP bound conformation. Exposed hydrophobic residues provide an interaction surface that may allow regulation of the GTP binding conformation, and thus activation of the GTPase, during the association of SRP with its receptor.

Microsecond exchange of internal water molecules in bacteriorhodopsin

The proton-conducting pathway of bacteriorhodopsin (BR) contains at least nine internal water molecules that are thought to be key players in the proton translocation mechanism. Here, we report the results of a multinuclear (1H, 2H, 17O) magnetic relaxation dispersion (MRD) study with the primary goal of determining the rate of exchange of these internal water molecules with bulk water. This rate is of interest in current attempts to elucidate the molecular details of the proton translocation mechanism. The relevance of water exchange kinetics is underscored by recent crystallographic findings of substantial variations in the number and locations of internal water molecules during the photocycle. Moreover, internal water exchange is believed to be governed by conformational fluctuations in the protein and can therefore provide information about the thermal accessibility of functionally important conformational substates. The present 2H and 17O MRD data show that at least seven water molecules, or more if they are orientationally disordered, in BR have residence times (inverse exchange rate constant) in the range 0.1-10 micros at 277 K. At least five of these water molecules have residence times in the more restrictive range 0.1-0.5 micros. These results show that most or all of the deeply buried water molecules in BR exchange on a time-scale that is short compared to the rate-limiting step in the photocycle. The MRD measurements were performed on BR solubilized in micelles of octyl glucoside. From the MRD data, the rotational correlation time of detergent-solubilized BR was determined to 35 ns at 300 K, consistent with a monomeric protein in complex with about 150 detergent molecules. The solubilized protein was found to be stable in the dark for at least eight months at 277 K.

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.