Role of glutamine-61 in the hydrolysis of GTP by p21H-ras: an experimental and theoretical study

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Mechanism-Based Redesign of GAP to Activate Oncogenic Ras

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

GAP positions catalytic H-Ras residue Q61 for GTP hydrolysis in molecular dynamics simulations, complicating chemical rescue of Ras deactivation

Functional interaction of Ras signaling proteins with upstream, negative regulatory GTPase activating proteins (GAPs) represents a crucial step in cellular decision making related to growth and survival. Key components of the catalytic transition state for Ras deactivation by GAP-accelerated hydrolysis of Ras-bound guanosine triphosphate (GTP) are thought to include an arginine residue from the GAP (the arginine finger), a glutamine residue from Ras (Q61), and a water molecule that is likely coordinated by Q61 to engage in nucleophilic attack on GTP. Here, we use in-vitro fluorescence experiments to show that 0.1-100 mM concentrations of free arginine, imidazole, and other small nitrogenous molecule fail to accelerate GTP hydrolysis, even in the presence of the catalytic domain of a mutant GAP lacking its arginine finger (R1276A NF1). This result is surprising given that imidazole can chemically rescue enzyme activity in arginine-to-alanine mutant protein tyrosine kinases (PTKs) that share many active site components with Ras/GAP complexes. Complementary all-atom molecular dynamics (MD) simulations reveal that an arginine finger GAP mutant still functions to enhance Ras Q61-GTP interaction, though less extensively than wild-type GAP. This increased Q61-GTP proximity may promote more frequent fluctuations into configurations that enable GTP hydrolysis as a component of the mechanism by which GAPs accelerate Ras deactivation in the face of arginine finger mutations. The failure of small molecule analogs of arginine to chemically rescue catalytic deactivation of Ras is consistent with the idea that the influence of the GAP goes beyond the simple provision of its arginine finger. However, the failure of chemical rescue in the presence of R1276A NF1 suggests that the GAPs arginine finger is either unsusceptible to rescue due to exquisite positioning or that it is involved in complex multivalent interactions. Therefore, in the context of oncogenic Ras proteins with mutations at codons 12 or 13 that inhibit arginine finger penetration toward GTP, drug-based chemical rescue of GTP hydrolysis may have bifunctional chemical/geometric requirements that are more difficult to satisfy than those that result from arginine-to-alanine mutations in other enzymes for which chemical rescue has been demonstrated.

Functional and biological heterogeneity of KRASQ61 mutations

Missense mutations at the three hotspots in the guanosine triphosphatase (GTPase) RAS-Gly12, Gly13, and Gln61 (commonly known as G12, G13, and Q61, respectively)-occur differentially among the three RAS isoforms. Q61 mutations in KRAS are infrequent and differ markedly in occurrence. Q61H is the predominant mutant (at 57%), followed by Q61R/L/K (collectively 40%), and Q61P and Q61E are the rarest (2 and 1%, respectively). Probability analysis suggested that mutational susceptibility to different DNA base changes cannot account for this distribution. Therefore, we investigated whether these frequencies might be explained by differences in the biochemical, structural, and biological properties of KRASQ61 mutants. Expression of KRASQ61 mutants in NIH 3T3 fibroblasts and RIE-1 epithelial cells caused various alterations in morphology, growth transformation, effector signaling, and metabolism. The relatively rare KRASQ61E mutant stimulated actin stress fiber formation, a phenotype distinct from that of KRASQ61H/R/L/P, which disrupted actin cytoskeletal organization. The crystal structure of KRASQ61E was unexpectedly similar to that of wild-type KRAS, a potential basis for its weak oncogenicity. KRASQ61H/L/R-mutant pancreatic ductal adenocarcinoma (PDAC) cell lines exhibited KRAS-dependent growth and, as observed with KRASG12-mutant PDAC, were susceptible to concurrent inhibition of ERK-MAPK signaling and of autophagy. Our results uncover phenotypic heterogeneity among KRASQ61 mutants and support the potential utility of therapeutic strategies that target KRASQ61 mutant-specific signaling and cellular output.

Enhanced BRAF engagement by NRAS mutants capable of promoting melanoma initiation

A distinct profile of NRAS mutants is observed in each tumor type. It is unclear whether these profiles are determined by mutagenic events or functional differences between NRAS oncoproteins. Here, we establish functional hallmarks of NRAS mutants enriched in human melanoma. We generate eight conditional, knock-in mouse models and show that rare melanoma mutants (NRAS G12D, G13D, G13R, Q61H, and Q61P) are poor drivers of spontaneous melanoma formation, whereas common melanoma mutants (NRAS Q61R, Q61K, or Q61L) induce rapid tumor onset with high penetrance. Molecular dynamics simulations, combined with cell-based protein-protein interaction studies, reveal that melanomagenic NRAS mutants form intramolecular contacts that enhance BRAF binding affinity, BRAF-CRAF heterodimer formation, and MAPK > ERK signaling. Along with the allelic series of conditional mouse models we describe, these results establish a mechanistic basis for the enrichment of specific NRAS mutants in human melanoma.

Identification of functional substates of KRas during GTP hydrolysis with enhanced sampling simulations

As the hub of major signaling pathways, Ras proteins are implicated in 19% of tumor-caused cancers due to perturbations in their conformational and/or catalytic properties. Despite numerous studies, the functions of the conformational substates for the most important isoform, KRas, remain elusive. In this work, we perform an extensive simulation analysis on the conformational landscape of KRas in its various chemical states during the GTP hydrolysis cycle: the reactant state KRasGTP·Mg²⁺, the intermediate state KRasGDP·Pi·Mg²⁺ and the product state KRasGDP·Mg²⁺. The results from enhanced sampling simulations reveal that State 1 of KRasGTP·Mg²⁺ has multiple stable substates in solution, one of which might account for interacting with GEFs. State 2 of KRasGTP·Mg²⁺ features two substates "Tyr32_in" and "Tyr32_out", which are poised to interact with effectors and GAPs, respectively. For the intermediate state KRasGDP·Pi·Mg²⁺, Gln61 and Pi are found to assume a broad set of conformations, which might account for the weak oncogenic effect of Gln61 mutations in KRas in contrast to the situation in HRas and NRas. Finally, the product state KRasGDP·Mg²⁺ has more than two stable substates in solution, pointing to a conformation-selection mechanism for complexation with GEFs. Based on these results, some specific inhibition strategies for targeting the binding sites of the high-energy substates of KRas during GTP hydrolysis are discussed.

Not all RAS mutations are equal: A detailed review of the functional diversity of RAS hot spot mutations

The RAS family of small GTPases are among the most frequently mutated oncogenes in human cancer. Approximately 20% of cancers harbor a RAS mutation, and >150 different missense mutations have been detected. Many of these mutations have mutant-specific biochemical defects that alter nucleotide binding and hydrolysis, effector interactions and cell signaling, prompting renewed efforts in the development of anti-RAS therapies, including the mutation-specific strategies. Previously viewed as undruggable, the recent FDA approval of a KRAS^G12C-selective inhibitor has offered real promise to the development of allele-specific RAS therapies. A broader understanding of the mutational consequences on RAS function must be developed to exploit additional allele-specific vulnerabilities. Approximately 94% of RAS mutations occur at one of three mutational "hot spots" at Gly¹², Gly¹³ and Gln⁶¹. Further, the single-nucleotide substitutions represent >99% of these mutations. Within this scope, we discuss the mutational frequencies of RAS isoforms in cancer, mutant-specific effector interactions and biochemical properties. By limiting our analysis to this mutational subset, we simplify the analysis while only excluding a small percentage of total mutations. Combined, these data suggest that the presence or absence of select RAS mutations in human cancers can be linked to their biochemical properties. Continuing to examine the biochemical differences in each RAS-mutant protein will continue to provide additional breakthroughs in allele-specific therapeutic strategies.

The Q61H mutation decouples KRAS from upstream regulation and renders cancer cells resistant to SHP2 inhibitors

Cancer cells bearing distinct KRAS mutations exhibit variable sensitivity to SHP2 inhibitors (SHP2i). Here we show that cells harboring KRAS Q61H are uniquely resistant to SHP2i, and investigate the underlying mechanisms using biophysics, molecular dynamics, and cell-based approaches. Q61H mutation impairs intrinsic and GAP-mediated GTP hydrolysis, and impedes activation by SOS1, but does not alter tyrosyl phosphorylation. Wild-type and Q61H-mutant KRAS are both phosphorylated by Src on Tyr32 and Tyr64 and dephosphorylated by SHP2, however, SHP2i does not reduce ERK phosphorylation in KRAS Q61H cells. Phosphorylation of wild-type and Gly12-mutant KRAS, which are associated with sensitivity to SHP2i, confers resistance to regulation by GAP and GEF activities and impairs binding to RAF, whereas the near-complete GAP/GEF-resistance of KRAS Q61H remains unaltered, and high-affinity RAF interaction is retained. SHP2 can stimulate KRAS signaling by modulating GEF/GAP activities and dephosphorylating KRAS, processes that fail to regulate signaling of the Q61H mutant.

SHP2 promotes RAS-driven MAPK signalling, but it is unclear why cancer cells with distinct KRAS mutations exhibit differential sensitivity to SHP2 inhibition. Here the authors show that KRAS Q61H is decoupled from SHP2- mediated upstream regulation, thus Q61H pancreatic cancer cells maintain MAPK signalling and are refractory to SHP2 inhibitors.

The retrospective molecular analysis of large or giant congenital melanocytic nevi in a group of Polish children

Background

Large and giant congenital melanocytic nevi (CMN), benign naevomelanocytic proliferations derived from neural crests, with a projected adult size (PAS) ≥ 20 cm, are connected to a high risk of melanoma and neurocutaneous melanosis. Among several factors, genetic alterations seem to be involved in tumorigenesis. The aim of the present study was to analyse the mutation status of NRAS and BRAF genes in resection specimens from large or giant CMN in a group of Polish patients.

Material and methods

The formalin-fixed, paraffin-embedded resection specimens from 18 patients, fixed in the years of 2006 to 2017, were included in the study. The regions containing the highest load of melanocytes were macrodissected prior to DNA isolation. The NRAS and BRAF mutation status was evaluated using qPCR.

Results

We detected activating mutations in NRAS gene (codons: 12 and 61) in 7 out of the 18 (38.9%) patients. No BRAF mutations were found.

Conclusion

Our study, the first molecular analysis of large/giant CMN in Polish patients, supports the hypothesis that NRAS mutation in codon 61 are frequent, recurrent mutations in large/giant CMN. Moreover, we show, for the first time, that NRAS mutations in codon 12 (p.Gly12Asp) can be also detected in giant CMN. The exact role of these genetic alterations in CMN formation remains to be elucidated.

Conformational Features of Ras: Key Hydrogen-Bonding Interactions of Gln61 in the Intermediate State during GTP Hydrolysis

The Ras protein is one of the most important drug targets for battling cancers. To effectively design novel drugs of Ras, we characterize here its conformational ensembles for the hydrolysis intermediate state RasGDP·Pi and the product state RasGDP by extensive replica-exchange molecular dynamics simulations. Several substates for RasGDP·Pi have been identified, while structural analyses have revealed an unrecognized hydrogen-bonding network that stabilizes the hydrolysis intermediate state. More interestingly, Gln61, which is involved in numerous oncogenic mutations, was found to be engaged in this hydrogen-bonding network, adopting a specific conformation that always points to Pi in contrast to that in the RasGTP state. The simulations also reveal that RasGDP has more than one substate, suggesting a conformational selection mechanism for the interaction between Ras and the guanine nucleotide exchange factors (GEFs). These findings offer new opportunities for the drug design of Ras by stabilizing the hydrolysis intermediate or disrupting its interaction with the GEFs.

The RAS GTPase RIT1 compromises mitotic fidelity through spindle assembly checkpoint suppression

The spindle assembly checkpoint (SAC) functions as a sensor of unattached kinetochores that delays mitotic progression into anaphase until proper chromosome segregation is guaranteed.^1,2 Disruptions to this safety mechanism lead to genomic instability and aneuploidy, which serve as the genetic cause of embryonic demise, congenital birth defects, intellectual disability, and cancer.^3,4 However, despite the understanding of the fundamental mechanisms that control the SAC, it remains unknown how signaling pathways directly interact with and regulate the mitotic checkpoint activity. In response to extracellular stimuli, a diverse network of signaling pathways involved in cell growth, survival, and differentiation are activated, and this process is prominently regulated by the Ras family of small guanosine triphosphatases (GTPases).⁵ Here we show that RIT1, a Ras-related GTPase that regulates cell survival and stress response,⁶ is essential for timely progression through mitosis and proper chromosome segregation. RIT1 dissociates from the plasma membrane (PM) during mitosis and interacts directly with SAC proteins MAD2 and p31^comet in a process that is regulated by cyclin-dependent kinase 1 (CDK1) activity. Furthermore, pathogenic levels of RIT1 silence the SAC and accelerate transit through mitosis by sequestering MAD2 from the mitotic checkpoint complex (MCC). Moreover, SAC suppression by pathogenic RIT1 promotes chromosome segregation errors and aneuploidy. Our results highlight a unique function of RIT1 compared to other Ras GTPases and elucidate a direct link between a signaling pathway and the SAC through a novel regulatory mechanism.

Differential Outcomes in Codon 12/13 and Codon 61 NRAS-Mutated Cancers in the Phase II NCI-MATCH Trial of Binimetinib in Patients with NRAS-Mutated Tumors

PURPOSE

Preclinical and clinical data suggest that downstream inhibition with an MEK inhibitor, such as binimetinib, might be efficacious for NRAS-mutated cancers.

PATIENTS AND METHODS

Patients enrolled in the NCI-MATCH trial master protocol underwent tumor biopsy and molecular profiling by targeted next-generation sequencing. Patients with NRAS-mutated tumors, except melanoma, were enrolled in subprotocol Z1A, a single-arm study evaluating binimetinib 45 mg twice daily. The primary endpoint was objective response rate (ORR). Secondary endpoints included progression-free survival (PFS) and overall survival (OS). A post hoc analysis examined the association of NRAS mutation type with outcome.

RESULTS

In total, 47 eligible patients with a refractory solid tumor harboring a codon 12, 13, or 61 NRAS mutation were treated. Observed toxicity was moderate, and 30% of patients discontinued treatment because of binimetinib-associated toxicity. The ORR was 2.1% (1/47 patients). A patient with malignant ameloblastoma harboring a codon 61 NRAS mutation achieved a durable partial response (PR). A patient with NRAS codon 61-mutated colorectal cancer had an unconfirmed PR, and two other patients with NRAS codon 61-mutated colorectal had stable disease for at least 12 months. In an exploratory analysis, patients with colorectal cancer bearing a NRAS codon 61 mutation (n = 8) had a significantly longer OS (P = 0.03) and PFS (P = 0.007) than those with codon 12 or 13 mutations (n = 16).

CONCLUSIONS

Single-agent binimetinib did not show promising efficacy in NRAS-mutated cancers. The observation of increased OS and PFS in patients with codon 61 NRAS-mutated colorectal cancer merits further investigation.

Antiparasitic dibenzalacetone inhibits the GTPase activity of Rab6 protein of Leishmania donovani (LdRab6), a potential target for its antileishmanial effect

Visceral leishmaniasis (VL, also known as kala-azar) is a vector borne disease caused by obligate intracellular protozoan parasite Leishmania donovani. To overcome the limitations of currently available drugs for VL, molecular target-based study is a promising tool to develop new drugs to treat this neglected tropical disease. One such target we recently identified from L. donovani (Ld) genome (WGS, clinical Indian isolate; BHU 1220, AVPQ01000001) is a small GTP-binding protein, Rab6 protein. We now report a specific inhibitor of the GTPase activity of Rab6 protein of L. donovani (LdRab6) without restricting host enzyme activity. First, to understand the nature of LdRab6 protein, we generated recombinant LdRab6 mutant proteins (rLdRab6) by systematically introducing deletion (two cysteine residues at C-terminal) and mutations [single amino acid substitutions in the conserved region of GTP (Q84L)/GDP(T38N) coding sequence]. The GTPase activity of rLdRab6:GTP and rLdRab6:GDP locked mutant proteins showed ~ 8-fold and ~ 1.5-fold decreases in enzyme activity, respectively, compared to the wild type enzyme activity. The mutant protein rLdRab6:ΔC inhibited the GTPase activity. Sequence alignment analysis of Rab6 protein of L. donovani with Homo sapiens showed identical amino acids in the G conserved region (GTP/GDP-binding sites) but it differed in the C-terminal region. We then evaluated the inhibitory activity of trans-dibenzalacetone (DBA, a synthetic analog of curcumin with strong antileishmanial activity reported earlier by us) in the GTPase activity of LdRab6 protein. Comparative molecular docking analysis of DBA and specific inhibitors of Rab proteins (Lovastatin, BFA, Zoledronate, and NE10790) indicated that DBA had optimum binding affinity with LdRab6 protein. This was further confirmed by the GTPase activity of DBA-treated LdRab6 which showed a basal GTP level significantly lower than that of the wild-type rLdRab6. The results confirm that DBA inhibits the GTPase activity of LdRab6 protein from L. donovani (LdRab6), a potential target for its antileishmanial effect.

Nutrient Signaling and Lysosome Positioning Crosstalk Through a Multifunctional Protein, Folliculin

FLCN was identified as the gene responsible for Birt-Hogg-Dubé (BHD) syndrome, a hereditary syndrome associated with the appearance of familiar renal oncocytomas. Most mutations affecting FLCN result in the truncation of the protein, and therefore loss of its associated functions, as typical for a tumor suppressor. FLCN encodes the protein folliculin (FLCN), which is involved in numerous biological processes; mutations affecting this protein thus lead to different phenotypes depending on the cellular context. FLCN forms complexes with two large interacting proteins, FNIP1 and FNIP2. Structural studies have shown that both FLCN and FNIPs contain longin and differentially expressed in normal versus neoplastic cells (DENN) domains, typically involved in the regulation of small GTPases. Accordingly, functional studies show that FLCN regulates both the Rag and the Rab GTPases depending on nutrient availability, which are respectively involved in the mTORC1 pathway and lysosomal positioning. Although recent structural studies shed light on the precise mechanism by which FLCN regulates the Rag GTPases, which in turn regulate mTORC1, how FLCN regulates membrane trafficking through the Rab GTPases or the significance of the intriguing FLCN-FNIP-AMPK complex formation are questions that still remain unanswered. We discuss the recent progress in our understanding of FLCN regulation of both growth signaling and lysosomal positioning, as well as future approaches to establish detailed mechanisms to explain the disparate phenotypes caused by the loss of FLCN function and the development of BHD-associated and other tumors.

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

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

Multi-target, ensemble-based virtual screening yields novel allosteric KRAS inhibitors at high success rate

RAS mutations account for >15% of all human tumors, and of these ~85% are due to mutations in a particular RAS gene: KRAS. Recent studies revealed that KRAS harbors four druggable allosteric sites. Here, we have (a) used molecular simulations to generate ensembles of wild type and four major oncogenic KRAS mutants (G12V, G12D, G13D, and Q61H); (b) characterized the druggability of each allosteric pocket in each protein; (c) conducted extensive ensemble-based virtual screening using pocket-tailored ligand libraries; (d) prioritized hits through hierarchical postdocking analysis; and (e) validated predicted hits with NMR. Of the 785 diverse potential hits identified by our in silico analysis, we tested 90 for their ability to bind KRAS using NMR and found that nine cause backbone amide chemical shift perturbations of residues near the functionally responsive switch loops, suggesting potential binding. We conducted detailed biophysical analyses on a novel indole-based compound to demonstrate the potential of our workflow to yield lead compounds. We believe the detailed information documented in this work regarding the druggability profile of each allosteric site and the chemical fingerprints of compounds that target them will serve as vital resources for future structure-based drug design efforts against KRAS, a high-value target for cancer therapy.

Macroenvironment-gene-microenvironment interactions in ultraviolet radiation-induced melanomagenesis

Cutaneous malignant melanoma is one of the few major cancers that continue to exhibit a positive rate of increase in the developed world. A wealth of epidemiological data has undisputedly implicated ultraviolet radiation (UVR) from sunlight and artificial sources as the major risk factor for melanomagenesis. However, the molecular mechanisms of this cause-and-effect relationship remain murky and understudied. Recent efforts on multiple fronts have brought unprecedented expansion of our knowledge base on this subject and it is now clear that melanoma is caused by a complex interaction between genetic predisposition and environmental exposure, primarily to UVR. Here we provide an overview of the effects of the macroenvironment (UVR) on the skin microenvironment and melanocyte-specific intrinsic (mostly genetic) landscape, which conspire to produce one of the deadliest malignancies.

Tyrosyl phosphorylation of KRAS stalls GTPase cycle via alteration of switch I and II conformation

Deregulation of the RAS GTPase cycle due to mutations in the three RAS genes is commonly associated with cancer development. Protein tyrosine phosphatase SHP2 promotes RAF-to-MAPK signaling pathway and is an essential factor in RAS-driven oncogenesis. Despite the emergence of SHP2 inhibitors for the treatment of cancers harbouring mutant KRAS, the mechanism underlying SHP2 activation of KRAS signaling remains unclear. Here we report tyrosyl-phosphorylation of endogenous RAS and demonstrate that KRAS phosphorylation via Src on Tyr32 and Tyr64 alters the conformation of switch I and II regions, which stalls multiple steps of the GTPase cycle and impairs binding to effectors. In contrast, SHP2 dephosphorylates KRAS, a process that is required to maintain dynamic canonical KRAS GTPase cycle. Notably, Src- and SHP2-mediated regulation of KRAS activity extends to oncogenic KRAS and the inhibition of SHP2 disrupts the phosphorylation cycle, shifting the equilibrium of the GTPase cycle towards the stalled ‘dark state’.

Deregulation of the RAS GTPase cycle due to mutations in RAS genes is commonly associated with cancer development. Here authors use NMR and mass spectrometry to shows that KRAS phosphorylation via Src alters the conformation of switch I and II regions and thereby impacts the GTPase cycle.

Quantitative Measurement of Intrinsic GTP Hydrolysis for Carcinogenic Glutamine 61 Mutants in H-Ras

Mutations of human oncoprotein p²¹H-Ras (hereafter "Ras") at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. However, there is no comprehensive data set of the effects of mutations to Q61 on the protein's intrinsic catalytic rate, structure, or interactions with water at the active site. Here, we present the first comprehensive and quantitative set of initial rates of intrinsic hydrolysis for all stable variants of RasQ61X. We further conducted enhanced molecular dynamics (MD) simulations of each construct to determine the solvent accessible surface area (SASA) of the side chain at position 61 and compared these results to previously measured changes in electric fields caused by RasQ61X mutations. For polar and negatively charged residues, we found that the rates are normally distributed about an optimal electrostatic contribution, close to that of the native Q61 residue, and the rates are strongly correlated to the number of waters in the active site. Together, these results support a mechanism of GTP hydrolysis in which Q61 stabilizes a transient hydronium ion, which then stabilizes the transition state while the γ-phosphate is undergoing nucleophilic attack by a second, catalytically active water molecule. We discuss the implications of such a mechanism on future strategies for combating Ras-based cancers.

Water Distribution within Wild-Type NRas Protein and Q61 Mutants during Unrestrained QM/MM Dynamics

Point mutations in p21^ras are associated with ∼30% of human tumors by disrupting its GTP hydrolysis cycle, which is critical to its molecular switch function in cellular signaling pathways. In this work, we investigate the impact of Gln 61 substitutions in the structure of the p21^N-ras active site and particularly focus on water reorganization around GTP, which appears to be crucial to evaluate favorable and unfavorable hydration sites for hydrolysis. The NRas-GTP complex is analyzed using a hybrid quantum mechanics/molecular mechanics approach, treating for the first time to our knowledge transient water molecules at the ab initio level and leading to results that account for the electrostatic coupling between the protein complex and the solvent. We show that for the wild-type protein, water molecules are found around the GTP γ-phosphate group, forming an arch extended from residues 12 to 35. Two density peaks are observed, supporting previous results that suggest the presence of two water molecules in the active site, one in the vicinity of residue 35 and a second one stabilized by hydrogen bonds formed with nitrogen backbone atoms of residues 12 and 60. The structural changes observed in NRas Gln 61 mutants result in the drastic delocalization of water molecules that we discuss. In mutants Q61H and Q61K, for which water distribution is overlocalized next to residue 60, the second density peak supports the hypothesis of a second water molecule. We also conclude that Gly 60 indirectly participates in GTP hydrolysis by correctly positioning transient water molecules in the protein complex and that Gln 61 has an indirect steric effect in stabilizing the preorganized catalytic site.

Intersubunit Crosstalk in the Rag GTPase Heterodimer Enables mTORC1 to Respond Rapidly to Amino Acid Availability

mTOR complex I (mTORC1) is a central growth regulator that senses amino acids through a pathway that converges on the Rag GTPases, an obligate heterodimer of two related GTPases. Despite their central role in amino acid sensing, it is unknown why the Rag GTPases are heterodimeric and whether their subunits communicate with each other. Here, we find that the binding of guanosine triphosphate (GTP) to one subunit inhibits the binding and induces the hydrolysis of GTP by the other. This intersubunit communication pushes the Rag GTPases into either of two stable configurations, which represent active "on" or "off" states that interconvert via transient intermediates. Subunit coupling confers on the mTORC1 pathway its capacity to respond rapidly to the amino acid level. Thus, the dynamic response of mTORC1 requires intersubunit communication by the Rag GTPases, providing a rationale for why they exist as a dimer and revealing a distinct mode of control for a GTP-binding protein.

Mechanistic Insights into the Differential Catalysis by RheB and Its Mutants: Y35A and Y35A-D65A

RheB GTPase is a Ras-related molecular switch, which regulates the mTOR signaling pathway by cycling between the active [guanosine triphosphate (GTP)] state and inactive [guanine diphosphate (GDP)] state. Impairment of GTPase activity because of mutations in several small GTPases is known to be associated with several cancers. The conventional GTPase mechanism such as in H-Ras requires a conserved glutamine (Q64) in the switch-II region of RheB to align the catalytic water molecule for efficient GTP hydrolysis. The conformation of this conserved glutamine is different in RheB, resulting in an altered conformation of the entire switch-II region. Studies on the atypical switch-II conformation in RheB revealed a distinct, noncanonical mode of GTP hydrolysis. An RheB mutant Y35A was previously shown to exclusively enhance the intrinsic GTPase activity of RheB, whereas the Y35A-D65A double mutant was shown to reduce the elevated GTPase activity. Here, we have used all-atom molecular dynamics (MD) simulations for comprehensive understanding of the conformational dynamics associated with the fast (Y35A) and slow (Y35A-D65A) hydrolyzing mutants of RheB. Using a combination of starting models from PDB structures and in-silico generated mutant structures, we discuss the observed conformational deviations in wild type (WT) versus mutants. Our results show that a number of interactions of RheB with phosphates of GTP as well as Mg²⁺ are destabilized in Y35A mutant in the switch-I region. We report distinct water dynamics at the active site of WT and mutants. Furthermore, principal component analysis showed significant differences in the conformational space sampled by the WT and mutants. Our observations provide improved understanding of the noncanonical GTP hydrolysis mechanism adopted by RheB and its modulation by Y35A and Y35A-D65A mutants.

Invited review: Small GTPases and their GAPs

Widespread utilization of small GTPases as major regulatory hubs in many different biological systems derives from a conserved conformational switch mechanism that facilitates cycling between GTP-bound active and GDP-bound inactive states under control of guanine nucleotide exchange factors (GEFs) and GTPase activating proteins (GAPs), which accelerate slow intrinsic rates of activation by nucleotide exchange and deactivation by GTP hydrolysis, respectively. Here we review developments leading to current understanding of intrinsic and GAP catalyzed GTP hydrolytic reactions in small GTPases from structural, molecular and chemical mechanistic perspectives. Despite the apparent simplicity of the GTPase cycle, the structural bases underlying the hallmark hydrolytic reaction and catalytic acceleration by GAPs are considerably more diverse than originally anticipated. Even the most fundamental aspects of the reaction mechanism have been challenging to decipher. Through a combination of experimental and in silico approaches, the outlines of a consensus view have begun to emerge for the best studied paradigms. Nevertheless, recent observations indicate that there is still much to be learned. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 431-448, 2016.

Biochemical Classification of Disease-associated Mutants of RAS-like Protein Expressed in Many Tissues (RIT1)

RAS-like protein expressed in many tissues 1 (RIT1) is a disease-associated RAS subfamily small guanosine triphosphatase (GTPase). Recent studies revealed that germ-line and somatic RIT1 mutations can cause Noonan syndrome (NS), and drive proliferation of lung adenocarcinomas, respectively, akin to RAS mutations in these diseases. However, the locations of these RIT1 mutations differ significantly from those found in RAS, and do not affect the three mutational "hot spots" of RAS. Moreover, few studies have characterized the GTPase cycle of RIT1 and its disease-associated mutants. Here we developed a real-time NMR-based GTPase assay for RIT1 and investigated the effect of disease-associated mutations on GTPase cycle. RIT1 exhibits an intrinsic GTP hydrolysis rate similar to that of H-RAS, but its intrinsic nucleotide exchange rate is ∼4-fold faster, likely as a result of divergent residues near the nucleotide binding site. All of the disease-associated mutations investigated increased the GTP-loaded, activated state of RIT1 in vitro, but they could be classified into two groups with different intrinsic GTPase properties. The S35T, A57G, and Y89H mutants exhibited more rapid nucleotide exchange, whereas F82V and T83P impaired GTP hydrolysis. A RAS-binding domain pulldown assay indicated that RIT1 A57G and Y89H were highly activated in HEK293T cells, whereas T83P and F82V exhibited more modest activation. All five mutations are associated with NS, whereas two (A57G and F82V) have also been identified in urinary tract cancers and myeloid malignancies. Characterization of the effects on the GTPase cycle of RIT1 disease-associated mutations should enable better understanding of their role in disease processes.

Molecular Mechanism for Conformational Dynamics of Ras·GTP Elucidated from In-Situ Structural Transition in Crystal

Ras•GTP adopts two interconverting conformational states, state 1 and state 2, corresponding to inactive and active forms, respectively. However, analysis of the mechanism for state transition was hampered by the lack of the structural information on wild-type Ras state 1 despite its fundamental nature conserved in the Ras superfamily. Here we solve two new crystal structures of wild-type H-Ras, corresponding to state 1 and state 2. The state 2 structure seems to represent an intermediate of state transition and, intriguingly, the state 1 crystal is successfully derived from this state 2 crystal by regulating the surrounding humidity. Structural comparison enables us to infer the molecular mechanism for state transition, during which a wide range of hydrogen-bonding networks across Switch I, Switch II and the α3-helix interdependently undergo gross rearrangements, where fluctuation of Tyr32, translocation of Gln61, loss of the functional water molecules and positional shift of GTP play major roles. The NMR-based hydrogen/deuterium exchange experiments also support this transition mechanism. Moreover, the unveiled structural features together with the results of the biochemical study provide a new insight into the physiological role of state 1 as a stable pool of Ras•GTP in the GDP/GTP cycle of Ras.

NRAS mutant melanoma: an overview for the clinician for melanoma management

Melanoma is the deadliest form of skin cancer and the incidence continues to rise in the United States and worldwide. Activating mutations in RAS oncogenes are found in roughly a third of all human cancers. Mutations in NRAS occur in approximately a fifth of cutaneous melanomas and are associated with aggressive clinical behavior. Cells harboring oncogenic NRAS mutations exhibit activation of multiple signaling cascades, including PI3K/Akt, MEK-ERK and RAL, which collectively stimulate cancer growth. While strategies to target N-Ras itself have proven ineffective, targeting one or more N-Ras effector pathways has shown promise in preclinical models. Despite promising preclinical data, current therapies for NRAS mutant melanoma remain limited. Immune checkpoint inhibitors and targeted therapies for BRAF mutant melanoma are transforming the treatment of metastatic melanoma, but the ideal treatment for NRAS mutant melanoma remains unknown. Improved understanding of NRAS mutant melanoma and relevant N-Ras effector signaling modules will be essential to develop new treatment strategies.

Immunohistochemical Detection of NRASQ61R Mutation in Diverse Tumor Types

OBJECTIVES

Testing for NRAS mutation at codon Q61 is of therapeutic, prognostic, and diagnostic importance for metastatic melanoma, thyroid carcinoma, and colorectal carcinoma. Immunohistochemistry for NRASQ61R, the most common NRAS mutation, offers several practical advantages over current molecular diagnostic techniques.

METHODS

We investigated the sensitivity and specificity of NRASQ61R in a series of 149 tumors with known NRAS genotype (72 malignant melanomas, 13 melanocytic nevi, 28 thyroid carcinomas, 25 gastrointestinal carcinomas, and 11 other malignancies).

RESULTS

Thirty-five cases harbored the NRASQ61R mutation (19 malignant melanomas, one melanocytic nevus, 10 thyroid carcinomas, two gastrointestinal carcinomas, and three other malignancies). In this series, the concordance rate between immunohistochemistry and mutational analyses was 100%. The sensitivity and specificity were 100% and 100%, respectively. However, lower staining intensity was observed for thyroid carcinomas in comparison to melanomas and other tumors.

CONCLUSIONS

Our studies confirmed that immunohistochemistry provides excellent sensitivity and specificity for detecting the NRASQ61R mutation in a variety of tumor types in a clinical setting.

Simultaneous and independent tuning of RhoA and Rac1 activity with orthogonally inducible promoters

The GTPases RhoA and Rac1 are key regulators of cell spreading, adhesion, and migration, and they exert distinct effects on the actin cytoskeleton. While RhoA classically stimulates stress fiber assembly and contraction, Rac1 promotes branched actin polymerization and membrane protrusion. These competing influences are reinforced by antagonistic crosstalk between RhoA and Rac1, which has complicated efforts to identify the specific mechanisms by which each GTPase regulates cell behavior. We therefore wondered whether RhoA and Rac1 are intrinsically coupled or whether they can be manipulated independently. To address this question, we placed constitutively active (CA) RhoA under a doxycycline-inducible promoter and CA Rac1 under an orthogonal cumate-inducible promoter, and we stably introduced both constructs into glioblastoma cells. We found that doxycycline addition increased RhoA activity without altering Rac1, and similarly cumate addition increased Rac1 activity without altering RhoA. Furthermore, co-expression of both mutants enabled high activation of RhoA and Rac1 simultaneously. When cells were cultured on collagen hydrogels, RhoA activation prevented cell spreading and motility, whereas Rac1 activation stimulated migration and dynamic cell protrusions. Interestingly, high activation of both GTPases induced a third phenotype, in which cells migrated at intermediate speeds similar to control cells but also aggregated into large, contractile clusters. In addition, we demonstrate dynamic and reversible switching between high RhoA and high Rac1 phenotypes. Overall, this approach represents a unique way to access different combinations of RhoA and Rac1 activity levels in a single cell and may serve as a valuable tool for multiplexed dissection and control of mechanobiological signals.

Modeling the mechanisms of biological GTP hydrolysis

Enzymes that hydrolyze GTP are currently in the spotlight, due to their molecular switch mechanism that controls many cellular processes. One of the best-known classes of these enzymes are small GTPases such as members of the Ras superfamily, which catalyze the hydrolysis of the γ-phosphate bond in GTP. In addition, the availability of an increasing number of crystal structures of translational GTPases such as EF-Tu and EF-G have made it possible to probe the molecular details of GTP hydrolysis on the ribosome. However, despite a wealth of biochemical, structural and computational data, the way in which GTP hydrolysis is activated and regulated is still a controversial topic and well-designed simulations can play an important role in resolving and rationalizing the experimental data. In this review, we discuss the contributions of computational biology to our understanding of GTP hydrolysis on the ribosome and in small GTPases.

Mutation-specific RAS oncogenicity explains NRAS codon 61 selection in melanoma

NRAS mutation at codons 12, 13, or 61 is associated with transformation; yet, in melanoma, such alterations are nearly exclusive to codon 61. Here, we compared the melanoma susceptibility of an NrasQ61R knock-in allele to similarly designed KrasG12D and NrasG12D alleles. With concomitant p16INK4a inactivation, KrasG12D or NrasQ61R expression efficiently promoted melanoma in vivo, whereas NrasG12D did not. In addition, NrasQ61R mutation potently cooperated with Lkb1/Stk11 loss to induce highly metastatic disease. Functional comparisons of NrasQ61R and NrasG12D revealed little difference in the ability of these proteins to engage PI3K or RAF. Instead, NrasQ61R showed enhanced nucleotide binding, decreased intrinsic GTPase activity, and increased stability when compared with NrasG12D. This work identifies a faithful model of human NRAS-mutant melanoma, and suggests that the increased melanomagenecity of NrasQ61R over NrasG12D is due to heightened abundance of the active, GTP-bound form rather than differences in the engagement of downstream effector pathways.

SIGNIFICANCE

This work explains the curious predominance in human melanoma of mutations of codon 61 of NRAS over other oncogenic NRAS mutations. Using conditional "knock-in" mouse models, we show that physiologic expression of NRASQ61R, but not NRASG12D, drives melanoma formation.

Structure-guided mutation of the conserved G3-box glycine in Rheb generates a constitutively activated regulator of mammalian target of rapamycin (mTOR)

Constitutively activated variants of small GTPases, which provide valuable functional probes of their role in cellular signaling pathways, can often be generated by mutating the canonical catalytic residue (e.g. Ras Q61L) to impair GTP hydrolysis. However, this general approach is ineffective for a substantial fraction of the small GTPase family in which this residue is not conserved (e.g. Rap) or not catalytic (e.g. Rheb). Using a novel engineering approach, we have manipulated nucleotide binding through structure-guided substitutions of an ultraconserved glycine residue in the G3-box motif (DXXG). Substitution of Rheb Gly-63 with alanine impaired both intrinsic and TSC2 GTPase-activating protein (GAP)-mediated GTP hydrolysis by displacing the hydrolytic water molecule, whereas introduction of a bulkier valine side chain selectively blocked GTP binding by steric occlusion of the γ-phosphate. Rheb G63A stimulated phosphorylation of the mTORC1 substrate p70S6 kinase more strongly than wild-type, thus offering a new tool for mammalian target of rapamycin (mTOR) signaling.

The DrrAB efflux system of Streptomyces peucetius is a multidrug transporter of broad substrate specificity

The soil bacterium Streptomyces peucetius produces two widely used anticancer antibiotics, doxorubicin and daunorubicin. Present within the biosynthesis gene cluster in S. peucetius is the drrAB operon, which codes for a dedicated ABC (ATP binding cassette)-type transporter for the export of these two closely related antibiotics. Because of its dedicated nature, the DrrAB system is believed to belong to the category of single-drug transporters. However, whether it also contains specificity for other known substrates of multidrug transporters has never been tested. In this study we demonstrate under both in vivo and in vitro conditions that the DrrAB system can transport not only doxorubicin but is also able to export two most commonly studied MDR substrates, Hoechst 33342 and ethidium bromide. Moreover, we demonstrate that many other substrates (including verapamil, vinblastine, and rifampicin) of the well studied multidrug transporters inhibit DrrAB-mediated Dox transport with high efficiency, indicating that they are also substrates of the DrrAB pump. Kinetic studies show that inhibition of doxorubicin transport by Hoechst 33342 and rifampicin occurs by a competitive mechanism, whereas verapamil inhibits transport by a non-competitive mechanism, thus suggesting the possibility of more than one drug binding site in the DrrAB system. This is the first in-depth study of a drug resistance system from a producer organism, and it shows that a dedicated efflux system like DrrAB contains specificity for multiple drugs. The significance of these findings in evolution of poly-specificity in drug resistance systems is discussed.

The folliculin tumor suppressor is a GAP for the RagC/D GTPases that signal amino acid levels to mTORC1

The mTORC1 kinase is a master growth regulator that senses numerous environmental cues, including amino acids. The Rag GTPases interact with mTORC1 and signal amino acid sufficiency by promoting the translocation of mTORC1 to the lysosomal surface, its site of activation. The Rags are unusual GTPases in that they function as obligate heterodimers, which consist of RagA or B bound to RagC or D. While the loading of RagA/B with GTP initiates amino acid signaling to mTORC1, the role of RagC/D is unknown. Here, we show that RagC/D is a key regulator of the interaction of mTORC1 with the Rag heterodimer and that, unexpectedly, RagC/D must be GDP bound for the interaction to occur. We identify FLCN and its binding partners, FNIP1/2, as Rag-interacting proteins with GAP activity for RagC/D, but not RagA/B. Thus, we reveal a role for RagC/D in mTORC1 activation and a molecular function for the FLCN tumor suppressor.

Ragulator is a GEF for the rag GTPases that signal amino acid levels to mTORC1

The mTOR Complex 1 (mTORC1) pathway regulates cell growth in response to numerous cues, including amino acids, which promote mTORC1 translocation to the lysosomal surface, its site of activation. The heterodimeric RagA/B-RagC/D GTPases, the Ragulator complex that tethers the Rags to the lysosome, and the v-ATPase form a signaling system that is necessary for amino acid sensing by mTORC1. Amino acids stimulate the binding of guanosine triphosphate to RagA and RagB but the factors that regulate Rag nucleotide loading are unknown. Here, we identify HBXIP and C7orf59 as two additional Ragulator components that are required for mTORC1 activation by amino acids. The expanded Ragulator has nucleotide exchange activity toward RagA and RagB and interacts with the Rag heterodimers in an amino acid- and v-ATPase-dependent fashion. Thus, we provide mechanistic insight into how mTORC1 senses amino acids by identifying Ragulator as a guanine nucleotide exchange factor (GEF) for the Rag GTPases.

An autoinhibited noncanonical mechanism of GTP hydrolysis by Rheb maintains mTORC1 homeostasis

Rheb, an activator of mammalian target of rapamycin (mTOR), displays low intrinsic GTPase activity favoring the biologically activated, GTP-bound state. We identified a Rheb mutation (Y35A) that increases its intrinsic nucleotide hydrolysis activity ∼10-fold, and solved structures of both its active and inactive forms, revealing an unexpected mechanism of GTP hydrolysis involving Asp65 in switch II and Thr38 in switch I. In the wild-type protein this noncanonical mechanism is markedly inhibited by Tyr35, which constrains the active site conformation, restricting the access of the catalytic Asp65 to the nucleotide-binding pocket. Rheb Y35A mimics the enthalpic and entropic changes associated with GTP hydrolysis elicited by the GTPase-activating protein (GAP) TSC2, and is insensitive to further TSC2 stimulation. Overexpression of Rheb Y35A impaired the regulation of mTORC1 signaling by growth factor availability. We demonstrate that the opposing functions of Tyr35 in the intrinsic and GAP-stimulated GTP catalysis are critical for optimal mTORC1 regulation.

Electrostatic effects of mutations of Ras glutamine 61 measured using vibrational spectroscopy of a thiocyanate probe

Mutations of human oncoprotein p21(Ras) (hereafter Ras) at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that electrostatically stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. We have tested the interactions between amino acids at this position and water by measuring changes in the electrostatic field experienced by a nitrile probe positioned near Ras Q61 using vibrational Stark effect (VSE) spectroscopy. We mutated this glutamine to every amino acid except cysteine and proline and then incubated these mutants with a Ral guanine nucleotide dissociation stimulator (Ral) containing the I18C mutation that was chemically labeled with a thiocyanate vibrational spectroscopic probe. The formation of the docked Ras Q61X-labeled Ral complex was confirmed by measurement of the dissociation constant of the interaction. We measured the absorption energy of this nitrile to determine any differences in electrostatic environment in the immediate vicinity of the thiocyanate probe between wild type and mutants of Ras. For each Ras Q61X mutant, we correlate the change in electrostatic field at position 61 with the solvent accessible surface area of polar components of the mutant side chain determined from a Boltzmann-weighted ensemble of structures, as well as the residue's hydration potential. These results support the hypothesis that the role of Ras Q61 is to stabilize water in or near the active site during GTP hydrolysis. The substantial effect that nonpolar side chains of Ras Q61X have on the absorption energy of the thiocyanate must be investigated with further experiments.

Ca2+-dependent GTPase, extra-large G protein 2 (XLG2), promotes activation of DNA-binding protein related to vernalization 1 (RTV1), leading to activation of floral integrator genes and early flowering in Arabidopsis

Heterotrimeric G proteins, consisting of Gα, Gβ, and Gγ subunits, play important roles in plant development and cell signaling. In Arabidopsis, in addition to one prototypical G protein α subunit, GPA1, there are three extra-large G proteins, XLG1, XLG2, and XLG3, of largely unknown function. Each extra-large G (XLG) protein has a C-terminal Gα-like region and a ∼400 amino acid N-terminal extension. Here we show that the three XLG proteins specifically bind and hydrolyze GTP, despite the fact that these plant-specific proteins lack key conserved amino acid residues important for GTP binding and hydrolysis of GTP in mammalian Gα proteins. Moreover, unlike other known Gα proteins, these activities require Ca(2+) instead of Mg(2+) as a cofactor. Yeast two-hybrid library screening and in vitro protein pull-down assays revealed that XLG2 interacts with the nuclear protein RTV1 (related to vernalization 1). Electrophoretic mobility shift assays show that RTV1 binds to DNA in vitro in a non-sequence-specific manner and that GTP-bound XLG2 promotes the DNA binding activity of RTV1. Overexpression of RTV1 results in early flowering. Combined overexpression of XLG2 and RTV1 enhances this early flowering phenotype and elevates expression of the floral pathway integrator genes, FT and SOC1, but does not repress expression of the floral repressor, FLC. Chromatin immunoprecipitation assays show that XLG2 increases RTV1 binding to FT and SOC1 promoters. Thus, a Ca(2+)-dependent G protein, XLG2, promotes RTV1 DNA binding activity for a subset of floral integrator genes and contributes to floral transition.

The Role of Gln61 in HRas GTP hydrolysis: a quantum mechanics/molecular mechanics study

Activation of the water molecule involved in GTP hydrolysis within the HRas·RasGAP system is analyzed using a tailored approach based on hybrid quantum mechanics/molecular mechanics (QM/MM) simulation. A new path emerges: transfer of a proton from the attacking water molecule to a second water molecule, then a different proton is transferred from this second water molecule to the GTP. Gln(61) will stabilize the transient OH(-) and H(3)O(+) molecules thus generated. This newly proposed mechanism was generated by using, for the first time to our knowledge, the entire HRas-RasGAP protein complex in a QM/MM simulation context. It also offers a rational explanation for previous experimental results regarding the decrease of GTPase rate found in the HRas Q61A mutant and the increase exhibited by the HRas Q61E mutant.

Analysis of Rho GTPase function in axon pathfinding using Caenorhabditis elegans

We provide information and protocols for the analysis of Rho GTPase function in axon pathfinding in Caenorhabditis elegans. The powerful molecular, genetic, imaging, and transgenic tools available in C. elegans make it an excellent system in which to study the in vivo roles of Rho GTPases. Methods for imaging of axon morphology in Rho GTPase single and double mutants are provided, as well as methods for the construction of transgenic C. elegans strains carrying exogenously introduced transgenes that drive the expression of constitutively active and dominant negative mutants.

New insight into the dynamic properties and the active site architecture of H-Ras p21 revealed by X-ray crystallography at very high resolution

Background

In kinetic crystallography, the usually static method of X-ray diffraction is expanded to allow time-resolved analysis of conformational rearrangements in protein structures. To achieve this, reactions have to be triggered within the protein crystals of interest, and optical spectroscopy can be used to monitor the reaction state. For this approach, a modified form of H-Ras p21 was designed which allows reaction initiation and fluorescence readout of the initiated GTPase reaction within the crystalline state. Rearrangements within the crystallized protein due to the progressing reaction and associated heterogeneity in the protein conformations have to be considered in the subsequent refinement processes.

Results

X-ray diffraction experiments on H-Ras p21 in different states along the reaction pathway provide detailed information about the kinetics and mechanism of the GTPase reaction. In addition, a very high data quality of up to 1.0 Å resolution allowed distinguishing two discrete subconformations of H-Ras p21, expanding the knowledge about the intrinsic flexibility of Ras-like proteins, which is important for their function. In a complex of H-Ras•GppNHp (guanosine-5'-(β,γ-imido)-triphosphate), a second Mg²⁺ ion was found to be coordinated to the γ-phosphate group of GppNHp, which positions the hydrolytically active water molecule very close to the attacked γ-phosphorous atom.

Conclusion

For the structural analysis of very high-resolution data we have used a new 'two-chain-isotropic-refinement' strategy. This refinement provides an alternative and easy to interpret strategy to reflect the conformational variability within crystal structures of biological macromolecules. The presented fluorescent form of H-Ras p21 will be advantageous for fluorescence studies on H-Ras p21 in which the use of fluorescent nucleotides is not feasible.

The distinct conformational dynamics of K-Ras and H-Ras A59G

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species.

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.

Oncogenic HRAS mutations cause prolonged PI3K signaling in response to epidermal growth factor in fibroblasts of patients with Costello syndrome

Costello syndrome (CS) is a rare congenital disorder characterized by failure to thrive, craniofacial dysmorphisms, cardiac and skin abnormalities, mental retardation, and predisposition to malignancies. CS is caused by heterozygous gain-of-function mutations in HRAS that also occur as somatic alterations in human tumors. HRAS is one of the three classical RAS proteins and cycles between an active, GTP- and an inactive, GDP-bound conformation. We used primary human skin fibroblasts from patients with CS as a model system to study the functional consequences of HRAS mutations on endogenous signaling pathways. The GTP-bound form of HRAS was significantly enriched in CS compared with normal fibroblasts. Active HRAS is known to stimulate both the RAF-MEK-ERK and the PI3K-AKT signaling cascade. Phosphorylation of MEK and ERK was normal in CS fibroblasts under basal conditions and slightly prolonged after epidermal growth factor (EGF) stimulation. Interestingly, basal phosphorylation of AKT was increased yet more in CS fibroblasts. Moreover, AKT phosphorylation was diminished in the early and enhanced in the late phase of EGF stimulation. Taken together, these results document that CS-associated HRAS mutations result in prolonged signal flux in a ligand-dependent manner. Our data suggest that altered cellular response to growth factors rather than constitutive activation of HRAS downstream signaling molecules may contribute to some of the clinical features in patients with CS.

Targeting the purinome

Purines are critical cofactors in the enzymatic reactions that create and maintain living organisms. In humans, there are approximately 3,266 proteins that utilize purine cofactors and these proteins constitute the so-called purinome. The human purinome encompasses a wide-ranging functional repertoire and many of these proteins are attractive drug targets. For example, it is estimated that 30% of modern drug discovery projects target protein kinases and that modulators of small G-proteins comprise more than 50% of currently marketed drugs. Given the importance of purine-binding proteins to drug discovery, the following review will discuss the forces that mediate protein:purine recognition, the factors that determine druggability of a protein target, and the process of structure-based drug design. A review of purine recognition in representatives of the various purine-binding protein families, as well as the challenges faced in targeting members of the purinome in drug discovery campaigns will also be given.

AtNOS/AtNOA1 is a functional Arabidopsis thaliana cGTPase and not a nitric-oxide synthase

AtNOS1 was previously identified as a potential nitric-oxide synthase (NOS) in Arabidopsis thaliana, despite lack of sequence similarity to animal NOSs. Although the dwarf and yellowish leaf phenotype of Atnos1 knock-out mutant plants can be rescued by treatment with exogenous NO, doubts have recently been raised as to whether AtNOS1 is a true NOS. Moreover, depending on the type of physiological responses studied, Atnos1 is not always deficient in NO induction and/or detection, as previously reported. Here, we present experimental evidence showing that AtNOS1 is unable to bind and oxidize arginine to NO. These results support the argument that AtNOS1 is not a NOS. We also show that the renamed NO-associated protein 1 (AtNOA1) is a member of the circularly permuted GTPase family (cGTPase). AtNOA1 specifically binds GTP and hydrolyzes it. Complementation experiments of Atnoa1 mutant plants with different constructs of AtNOA1 show that GTP hydrolysis is necessary but not sufficient for the physiological function of AtNOA1. Mutant AtNOA1 lacking the C-terminal domain, although retaining GTPase activity, failed to complement Atnoa1, suggesting that this domain plays a crucial role in planta. cGTPases appear to be RNA-binding proteins, and the closest homolog of AtNOA1, the Bacillus subtilis YqeH, has been shown to participate in ribosome assembly and stability. We propose a similar function for AtNOA1 and discuss it in the light of its potential role in NO accumulation and plant development.

Salmonella enterica serovar Typhimurium BipA exhibits two distinct ribosome binding modes

BipA is a highly conserved prokaryotic GTPase that functions to influence numerous cellular processes in bacteria. In Escherichia coli and Salmonella enterica serovar Typhimurium, BipA has been implicated in controlling bacterial motility, modulating attachment and effacement processes, and upregulating the expression of virulence genes and is also responsible for avoidance of host defense mechanisms. In addition, BipA is thought to be involved in bacterial stress responses, such as those associated with virulence, temperature, and symbiosis. Thus, BipA is necessary for securing bacterial survival and successful invasion of the host. Steady-state kinetic analysis and pelleting assays were used to assess the GTPase and ribosome-binding properties of S. enterica BipA. Under normal bacterial growth, BipA associates with the ribosome in the GTP-bound state. However, using sucrose density gradients, we demonstrate that the association of BipA and the ribosome is altered under stress conditions in bacteria similar to those experienced during virulence. The data show that this differential binding is brought about by the presence of ppGpp, an alarmone that signals the onset of stress-related events in bacteria.

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.

Insight into catalysis of a unique GTPase reaction by a combined biochemical and FTIR approach

Rap1 and Rap2 are the only small guanine nucleotide-binding proteins of the Ras superfamily that do not use glutamine for GTP hydrolysis. Moreover, Rap1GAP, which stimulates the GTPase reaction of Rap1 10(5)-fold, does not have the classical "arginine finger" like RasGAP but presumably, introduces an asparagine residue into the active site. Here, we address the requirements of this unique reaction in detail by combining various biochemical methods, such as fluorescence spectroscopy, stopped-flow and time-resolved Fourier transform infrared spectroscopy (FTIR). The fluorescence spectroscopic assay monitors primarily protein-protein interaction steps, while FTIR resolves simultaneously the elementary steps of functional groups labor-free, but it is less sensitive and needs higher concentrations. Combining both methods allows us to distinguish weather mechanistic defects caused by mutation are due to affinity or due to functionality. We show that several mutations of Asn290 block catalysis. Some of the mutants, however, still form a complex with Rap1*GDP in the presence of BeF(x) but not AlF(x), supporting the notion that fluoride complexes are indicators of the ground versus transition state. Mutational analysis also shows that Thr61 is not required for catalysis. While replacement of Thr61 of Rap1 by Leu eliminates GTPase activation by Rap1GAP, the T61A and T61Q mutants have only a minor effect on catalysis, but change the relative rates of cleavage and (P(i)(-)) release. While Rap1GAP(N290A) is completely inactive on wild-type Rap1, it can act on Rap1(T61Q), arguing that Asn290 in trans has a role in catalysis similar to that of the intrinsic Gln in Ras and Rho. Finally, since FTIR works at high, and thus mostly saturating, concentrations, it can clearly separate effects on affinity from purely catalytic modifications, showing that Arg388, conserved between RapGAPs and mutated in the homologous RheBGAP Tuberin, affects binding affinity severely but has no effect on the cleavage reaction itself.