The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

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.

Common Patterns of Hydrolysis Initiation in P-loop Fold Nucleoside Triphosphatases

The P-loop fold nucleoside triphosphate (NTP) hydrolases (also known as Walker NTPases) function as ATPases, GTPases, and ATP synthases, are often of medical importance, and represent one of the largest and evolutionarily oldest families of enzymes. There is still no consensus on their catalytic mechanism. To clarify this, we performed the first comparative structural analysis of more than 3100 structures of P-loop NTPases that contain bound substrate Mg-NTPs or their analogues. We proceeded on the assumption that structural features common to these P-loop NTPases may be essential for catalysis. Our results are presented in two articles. Here, in the first, we consider the structural elements that stimulate hydrolysis. Upon interaction of P-loop NTPases with their cognate activating partners (RNA/DNA/protein domains), specific stimulatory moieties, usually Arg or Lys residues, are inserted into the catalytic site and initiate the cleavage of gamma phosphate. By analyzing a plethora of structures, we found that the only shared feature was the mechanistic interaction of stimulators with the oxygen atoms of gamma-phosphate group, capable of causing its rotation. One of the oxygen atoms of gamma phosphate coordinates the cofactor Mg ion. The rotation must pull this oxygen atom away from the Mg ion. This rearrangement should affect the properties of the other Mg ligands and may initiate hydrolysis according to the mechanism elaborated in the second article.

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.

Structure-based inhibitor design of mutant RAS proteins-a paradigm shift

As a member of small GTPase family, KRAS protein is a key physiological modulator of various cellular activities including proliferation. However, mutations of KRAS present in numerous cancer types, most frequently in pancreatic (> 60%), colorectal (> 40%), and lung cancers, drive oncogenic processes through overactivation of proliferation. The G12C mutation of KRAS protein is especially abundant in the case of these types of malignancies. Despite its key importance in human disease, KRAS was assumed to be non-druggable for a long time since the protein seemingly lacks potential drug-binding pockets except the nucleotide-binding site, which is difficult to be targeted due to the high affinity of KRAS for both GDP and GTP. Recently, a new approach broke the ice and provided evidence that upon covalent targeting of the G12C mutant KRAS, a highly dynamic pocket was revealed. This novel targeting is especially important since it serves with an inherent solution for drug selectivity. Based on these results, various structure-based drug design projects have been launched to develop selective KRAS mutant inhibitors. In addition to the covalent modification strategy mostly applicable for G12C mutation, different innovative solutions have been suggested for the other frequently occurring oncogenic G12 mutants. Here we summarize the latest advances of this field, provide perspectives for novel approaches, and highlight the special properties of KRAS, which might issue some new challenges.

Residue interaction networks of K-Ras protein with water molecules identifies the potential role of switch II and P-loop

The mutant K-Ras with aberrant signaling is the primary cause of several cancers. The proposed study investigated the influence of water molecules in K-Ras crystal structure, where they have a significant function by understanding their residue interaction networks (RINs). We analyzed the RINs of K-Ras with and without water molecules and determined their interaction properties. RINs were developed with the help of StructureViz2 and RINspector; further, the changes in K-Ras backbone flexibility were predicted with the DynaMine. We found that the residues K42, I142, and L159 are the hotspots from water, including the K-Ras-GTP complex with the highest residue centrality analysis (RCA) Z-score. The DynaMine prediction calculated the NMR S² value for the frequently mutated positions G12, G13, and Q61 showing a minor shift in flexibility, which make up the P-Loop and switch II of the K-Ras protein. This flexibility shift can account for changes in conformational activity and the protein's GTPase activity, making it difficult to recognize by the effectors and exchange factors. Taken together, our study helps in understanding the functional importance of the water molecules in K-Ras protein and the impact of mutation that modulate the conformational state of the protein.

Ras Isoforms from Lab Benches to Lives-What Are We Missing and How Far Are We?

The central protein in the oncogenic circuitry is the Ras GTPase that has been under intense scrutiny for the last four decades. From its discovery as a viral oncogene and its non-oncogenic contribution to crucial cellular functioning, an elaborate genetic, structural, and functional map of Ras is being created for its therapeutic targeting. Despite decades of research, there still exist lacunae in our understanding of Ras. The complexity of the Ras functioning is further exemplified by the fact that the three canonical Ras genes encode for four protein isoforms (H-Ras, K-Ras4A, K-Ras4B, and N-Ras). Contrary to the initial assessment that the H-, K-, and N-Ras isoforms are functionally similar, emerging data are uncovering crucial differences between them. These Ras isoforms exhibit not only cell-type and context-dependent functions but also activator and effector specificities on activation by the same receptor. Preferential localization of H-, K-, and N-Ras in different microdomains of the plasma membrane and cellular organelles like Golgi, endoplasmic reticulum, mitochondria, and endosome adds a new dimension to isoform-specific signaling and diverse functions. Herein, we review isoform-specific properties of Ras GTPase and highlight the importance of considering these towards generating effective isoform-specific therapies in the future.

Mechanistic insights into the effect of phosphorylation on Ras conformational dynamics and its interactions with cell signaling proteins

Ras undergoes interconversion between the active GTP-bound state and the inactive GDP-bound state. This GTPase cycle, which controls the activities of Ras, is accelerated by Ras GTPase-activating proteins (GAPs) and guanine nucleotide exchange factors (SOS). Oncogenic Ras mutations could affect the GTPase cycle and impair Ras functions. Additionally, Src-induced K-Ras Y32/64 dual phosphorylation has been reported to disrupt GTPase cycle and hinder Ras downstream signaling. However, the underlying mechanisms remain unclear. To address this, we performed molecular dynamics simulations (~30 μs in total) on unphosphorylated and phosphorylated K-Ras4B in GTP- and GDP-bound states, and on their complexes with GTPase cycle regulators (GAP and SOS) and the effector protein Raf. We found that K-Ras4B dual phosphorylation mainly alters the conformation at the nucleotide binding site and creates disorder at the catalytic site, resulting in the enlargement of GDP binding pocket and the retard of Ras-GTP intrinsic hydrolysis. We observed phosphorylation-induced shift in the distribution of Ras-GTP inactive-active sub-states and recognized potential druggable pockets in the phosphorylated Ras-GTP. Moreover, decreased catalytic competence or signal delivery abilities due to reduced binding affinities and/or distorted catalytic conformations of GAP, SOS and Raf were observed. In addition, the allosteric pathway from Ras/Raf interface to the distal Raf L4 loop was compromised by Ras phosphorylation. These results reveal the mechanisms by which phosphorylation influences the intrinsic or GAP/SOS catalyzed transformations between GTP- and GDP-bound states of Ras and its signal transduction to Raf. Our findings project Ras phosphorylation as a target for cancer drug discovery.

Understanding the Biological Activities of Vitamin D in Type 1 Neurofibromatosis: New Insights into Disease Pathogenesis and Therapeutic Design

Vitamin D is a fat-soluble steroid hormone playing a pivotal role in calcium and phosphate homeostasis as well as in bone health. Vitamin D levels are not exclusively dependent on food intake. Indeed, the endogenous production-occurring in the skin and dependent on sun exposure-contributes to the majority amount of vitamin D present in the body. Since vitamin D receptors (VDRs) are ubiquitous and drive the expression of hundreds of genes, the interest in vitamin D has tremendously grown and its role in different diseases has been extensively studied. Several investigations indicated that vitamin D action extends far beyond bone health and calcium metabolism, showing broad effects on a variety of critical illnesses, including cancer, infections, cardiovascular and autoimmune diseases. Epidemiological studies indicated that low circulating vitamin D levels inversely correlate with cutaneous manifestations and bone abnormalities, clinical hallmarks of neurofibromatosis type 1 (NF1). NF1 is an autosomal dominant tumour predisposition syndrome causing significant pain and morbidity, for which limited treatment options are available. In this context, vitamin D or its analogues have been used to treat both skin and bone lesions in NF1 patients, alone or combined with other therapeutic agents. Here we provide an overview of vitamin D, its characteristic nutritional properties relevant for health benefits and its role in NF1 disorder. We focus on preclinical and clinical studies that demonstrated the clinical correlation between vitamin D status and NF1 disease, thus providing important insights into disease pathogenesis and new opportunities for targeted therapy.

Pumping the brakes on RAS - negative regulators and death effectors of RAS

Mutations that activate the RAS oncoproteins are common in cancer. However, aberrant upregulation of RAS activity often occurs in the absence of activating mutations in the RAS genes due to defects in RAS regulators. It is now clear that loss of function of Ras GTPase-activating proteins (RasGAPs) is common in tumors, and germline mutations in certain RasGAP genes are responsible for some clinical syndromes. Although regulation of RAS is central to their activity, RasGAPs exhibit great diversity in their binding partners and therefore affect signaling by multiple mechanisms that are independent of RAS. The RASSF family of tumor suppressors are essential to RAS-induced apoptosis and senescence, and constitute a barrier to RAS-mediated transformation. Suppression of RASSF protein expression can also promote the development of excessive RAS signaling by uncoupling RAS from growth inhibitory pathways. Here, we will examine how these effectors of RAS contribute to tumor suppression, through both RAS-dependent and RAS-independent mechanisms.

The road to ERK activation: Do neurons take alternate routes?

The ERK cascade is a central signaling pathway that regulates a wide variety of cellular processes including proliferation, differentiation, learning and memory, development, and synaptic plasticity. A wide range of inputs travel from the membrane through different signaling pathway routes to reach activation of one set of output kinases, ERK1&2. The classical ERK activation pathway beings with growth factor activation of receptor tyrosine kinases. Numerous G-protein coupled receptors and ionotropic receptors also lead to ERK through increases in the second messengers calcium and cAMP. Though both types of pathways are present in diverse cell types, a key difference is that most stimuli to neurons, e.g. synaptic inputs, are transient, on the order of milliseconds to seconds, whereas many stimuli acting on non-neural tissue, e.g. growth factors, are longer duration. The ability to consolidate these inputs to regulate the activation of ERK in response to diverse signals raises the question of which factors influence the difference in ERK activation pathways. This review presents both experimental studies and computational models aimed at understanding the control of ERK activation and whether there are fundamental differences between neurons and other cells. Our main conclusion is that differences between cell types are quite subtle, often related to differences in expression pattern and quantity of some molecules such as Raf isoforms. In addition, the spatial location of ERK is critical, with regulation by scaffolding proteins producing differences due to colocalization of upstream molecules that may differ between neurons and other cells.

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.

Controlling the Revolving and Rotating Motion Direction of Asymmetric Hexameric Nanomotor by Arginine Finger and Channel Chirality

Nanomotors in nanotechnology are as important as engines in daily life. Many ATPases are nanoscale biomotors classified into three categories based on the motion mechanisms in transporting substrates: linear, rotating, and the recently discovered revolving motion. Most biomotors adopt a multisubunit ring-shaped structure that hydrolyzes ATP to generate force. How these biomotors control the motion direction and regulate the sequential action of their multiple subunits is intriguing. Many ATPases are hexameric with each monomer containing a conserved arginine finger. This review focuses on recent findings on how the arginine finger controls motion direction and coordinates adjacent subunit interactions in both revolving and rotating biomotors. Mechanisms of intersubunit interactions and sequential movements of individual subunits are evidenced by the asymmetrical appearance of one dimer and four monomers in high-resolution structural complexes. The arginine finger is situated at the interface of two subunits and extends into the ATP binding pocket of the downstream subunit. An arginine finger mutation results in deficiency in ATP binding/hydrolysis, substrate binding, and transport, highlighting the importance of the arginine finger in regulating energy transduction and motor function. Additionally, the roles of channel chirality and channel size are discussed as related to controlling one-way trafficking and differentiating the revolving and rotating mechanisms. Finally, the review concludes by discussing the conformational changes and entropy conversion triggered by ATP binding/hydrolysis, offering a view different from the traditional concept of ATP-mediated mechanochemical energy coupling. The elucidation of the motion mechanism and direction control in ATPases could facilitate nanomotor fabrication in nanotechnology.

Hybrid QM/MM vs Pure MM Molecular Dynamics for Evaluating Water Distribution within p21N-ras and the Resulting GTP Electronic Density

p21^ras protein activity, regulated by GTP hydrolysis, constitutes an active field of research for the development of cancer targeted therapies that would concern ∼30% of human tumors to which specific mutations have been associated. Indeed, the catalyzing mechanisms provided by the protein environment during GTP hydrolysis and how they are impaired by specific mutations remain to be fully elucidated. In this article, we present results from molecular mechanics (MM) molecular dynamics (MD) simulations and density functional theory (DFT) calculations carried out for wild-type p21 ^N-ras and six Gln 61 mutants. In the first part, we present the water distribution within the active site of the wild-type protein according to MM MD. Significant differences are observed when comparing the results to the previous distribution assessed through quantum mechanics/molecular mechanics (QM/MM) MD. Such method-dependent results highlight the importance of accounting for the electrostatic coupling between the protein complex and the solvent molecules in identifying hydration sites. In the second part, we present the results from DFT calculations performed to determine the electronic distribution of the GTP ligand, considering the wild-type active site arrangement according to both classical and hybrid approaches. Only in the QM/MM-based configuration is the ligand electronic density similar to that of a GDP-like state observed experimentally. For this reason, in the last set of calculations carried out for p21 ^N-ras Gln 61 mutants, only the active site structural conformations obtained through hybrid MD are considered. Through the analysis of the GTP electronic density, we conclude that the wild-type active site arrangement according to QM/MM MD is closer to a catalytically efficient conformation of the protein than the arrangement according to MM MD. Hence, water distribution according to the hybrid approach must correspond to the optimal placement of solvent in the active site. Within all of the studied Gln 61 substituted proteins, p21^ras major catalyzing effect, which consists of stabilizing a more GDP-like state, is lost.

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

Computational insights of K1444N substitution in GAP-related domain of NF1 gene associated with neurofibromatosis type 1 disease: a molecular modeling and dynamics approach

The NF1 gene encodes for neurofibromin protein, which is ubiquitously expressed, but most highly in the central nervous system. Non-synonymous SNPs (nsSNPs) in the NF1 gene were found to be associated with Neurofibromatosis Type 1 disease, which is characterized by the growth of tumors along nerves in the skin, brain, and other parts of the body. In this study, we used several in silico predictions tools to analyze 16 nsSNPs in the RAS-GAP domain of neurofibromin, the K1444N (K1423N) mutation was predicted as the most pathogenic. The comparative molecular dynamic simulation (MDS; 50 ns) between the wild type and the K1444N (K1423N) mutant suggested a significant change in the electrostatic potential. In addition, the RMSD, RMSF, Rg, hydrogen bonds, and PCA analysis confirmed the loss of flexibility and increase in compactness of the mutant protein. Further, SASA analysis revealed exchange between hydrophobic and hydrophilic residues from the core of the RAS-GAP domain to the surface of the mutant domain, consistent with the secondary structure analysis that showed significant alteration in the mutant protein conformation. Our data concludes that the K1444N (K1423N) mutant lead to increasing the rigidity and compactness of the protein. This study provides evidence of the benefits of the computational tools in predicting the pathogenicity of genetic mutations and suggests the application of MDS and different in silico prediction tools for variant assessment and classification in genetic clinics.

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.

Review: Ras GTPases and myosin: Qualitative conservation and quantitative diversification in signal and energy transduction

Most GTPases and many ATPases belong to the P‐loop class of proteins with significant structural and mechanistic similarities. Here we compare and contrast the basic properties of the Ras family GTPases and myosin, and conclude that there are fundamental similarities but also distinct differences related to their specific roles. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 422–430, 2016.

The Structural Basis of Oncogenic Mutations G12, G13 and Q61 in Small GTPase K-Ras4B

Ras mediates cell proliferation, survival and differentiation. Mutations in K-Ras4B are predominant at residues G12, G13 and Q61. Even though all impair GAP-assisted GTP → GDP hydrolysis, the mutation frequencies of K-Ras4B in human cancers vary. Here we aim to figure out their mechanisms and differential oncogenicity. In total, we performed 6.4 μs molecular dynamics simulations on the wild-type K-Ras4B (K-Ras4B^WT-GTP/GDP) catalytic domain, the K-Ras4B^WT-GTP–GAP complex, and the mutants (K-Ras4B^{G12C/G12D/G12V}-GTP/GDP, K-Ras4B^G13D-GTP/GDP, K-Ras4B^Q61H-GTP/GDP) and their complexes with GAP. In addition, we simulated ‘exchanged’ nucleotide states. These comprehensive simulations reveal that in solution K-Ras4B^WT-GTP exists in two, active and inactive, conformations. Oncogenic mutations differentially elicit an inactive-to-active conformational transition in K-Ras4B-GTP; in K-Ras4B^G12C/G12D-GDP they expose the bound nucleotide which facilitates the GDP-to-GTP exchange. These mechanisms may help elucidate the differential mutational statistics in K-Ras4B-driven cancers. Exchanged nucleotide simulations reveal that the conformational transition is more accessible in the GTP-to-GDP than in the GDP-to-GTP exchange. Importantly, GAP not only donates its R789 arginine finger, but stabilizes the catalytically-competent conformation and pre-organizes catalytic residue Q61; mutations disturb the R789/Q61 organization, impairing GAP-mediated GTP hydrolysis. Together, our simulations help provide a mechanistic explanation of key mutational events in one of the most oncogenic proteins in cancer.

Ras Conformational Ensembles, Allostery, and Signaling

Ras proteins are classical members of small GTPases that function as molecular switches by alternating between inactive GDP-bound and active GTP-bound states. Ras activation is regulated by guanine nucleotide exchange factors that catalyze the exchange of GDP by GTP, and inactivation is terminated by GTPase-activating proteins that accelerate the intrinsic GTP hydrolysis rate by orders of magnitude. In this review, we focus on data that have accumulated over the past few years pertaining to the conformational ensembles and the allosteric regulation of Ras proteins and their interpretation from our conformational landscape standpoint. The Ras ensemble embodies all states, including the ligand-bound conformations, the activated (or inactivated) allosteric modulated states, post-translationally modified states, mutational states, transition states, and nonfunctional states serving as a reservoir for emerging functions. The ensemble is shifted by distinct mutational events, cofactors, post-translational modifications, and different membrane compositions. A better understanding of Ras biology can contribute to therapeutic strategies.

Computational characterization of the chemical step in the GTP hydrolysis by Ras-GAP for the wild-type and G13V mutated Ras

The free energy profiles for the chemical reaction of the guanosine triphosphate hydrolysis GTP + H2O → GDP + Pi by Ras-GAP for the wild-type and G13V mutated Ras were computed by using molecular dynamics protocols with the QM(ab initio)/MM potentials. The results are consistent with the recent measurements of reaction kinetics in Ras-GAP showing about two-order reduction of the rate constant upon G13V mutation in Ras: the computed activation barrier on the free energy profile is increased by 3 kcal/mol upon the G13V replacement. The major reason for a higher energy barrier is a shift of the "arginine finger" (R789 from GAP) from the favorable position in the active site. The results of simulations provide support for the mechanism of the reference reaction according to which the Q61 side chain directly participates in chemical transformations at the proton transfer stage.

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

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

Lessons from computer simulations of Ras proteins in solution and in membrane

BACKGROUND

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

SCOPE OF REVIEW

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

MAJOR CONCLUSIONS

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

GENERAL SIGNIFICANCE

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.

Assessment of the potential pathogenicity of missense mutations identified in the GTPase-activating protein (GAP)-related domain of the neurofibromatosis type-1 (NF1) gene

Neurofibromatosis type-1 (NF1) is caused by constitutional mutations of the NF1 tumor-suppressor gene. Although ∼85% of inherited NF1 microlesions constitute truncating mutations, the remaining ∼15% are missense mutations whose pathological relevance is often unclear. The GTPase-activating protein-related domain (GRD) of the NF1-encoded protein, neurofibromin, serves to define its major function as a negative regulator of the Ras-MAPK (mitogen-activated protein kinase) signaling pathway. We have established a functional assay to assess the potential pathogenicity of 15 constitutional nonsynonymous NF1 missense mutations (11 novel and 4 previously reported but not functionally characterized) identified in the NF1-GRD (p.R1204G, p.R1204W, p.R1276Q, p.L1301R, p.I1307V, p.T1324N, p.E1327G, p.Q1336R, p.E1356G, p.R1391G, p.V1398D, p.K1409E, p.P1412R, p.K1436Q, p.S1463F). Individual mutations were introduced into an NF1-GRD expression vector and activated Ras was assayed by an enzyme-linked immunosorbent assay (ELISA). Ten NF1-GRD variants were deemed to be potentially pathogenic by virtue of significantly elevated levels of activated GTP-bound Ras in comparison to wild-type NF1 protein. The remaining five NF1-GRD variants were deemed less likely to be of pathological significance as they exhibited similar levels of activated Ras to the wild-type protein. These conclusions received broad support from both bioinformatic analysis and molecular modeling and serve to improve our understanding of NF1-GRD structure and function.

Klebsiella aerogenes UreF: identification of the UreG binding site and role in enhancing the fidelity of urease activation

The Ni-containing active site of Klebsiella aerogenes urease is assembled through the concerted action of the UreD, UreE, UreF, and UreG accessory proteins. UreE functions as a metallochaperone that delivers Ni to a UreD-UreF-UreG complex bound to urease apoprotein, with UreG serving as a GTPase during enzyme activation. This study focuses on the role of UreF, previously proposed to act as a GTPase activating protein (GAP) of UreG. Sixteen conserved UreF surface residues that may play roles in protein-protein interactions were independently changed to Ala. When produced in the context of the entire urease gene cluster, cell-free extracts of nine site-directed mutants had less than 10% of the wild-type urease activity. Enrichment of the variant forms of UreF, as the UreE-F fusion proteins, uniformly resulted in copurification of UreD and urease apoprotein, whereas UreG bound to only a subset of the species. Notably, weakened interaction with UreG correlated with the low-activity mutants. The affected residues in UreF map to a distinct surface on the crystal structure, defining the UreG binding site. In contrast to the hypothesis that UreF is a GAP, the UreD-UreF-UreG-urease apoprotein complex containing K165A UreF exhibited significantly greater levels of GTPase activity than that containing the wild-type protein. Additional studies demonstrated the UreG GTPase activity was largely uncoupled from urease activation for the complex containing this UreF variant. Further experiments with these complexes provided evidence that UreF gates the GTPase activity of UreG to enhance the fidelity of urease metallocenter assembly, especially in the presence of the noncognate metal Zn.

GDP-bound and nucleotide-free intermediates of the guanine nucleotide exchange in the Rab5·Vps9 system

Many GTPases regulate intracellular transport and signaling in eukaryotes. Guanine nucleotide exchange factors (GEFs) activate GTPases by catalyzing the exchange of their GDP for GTP. Here we present crystallographic and biochemical studies of a GEF reaction with four crystal structures of Arabidopsis thaliana ARA7, a plant homolog of Rab5 GTPase, in complex with its GEF, VPS9a, in the nucleotide-free and GDP-bound forms, as well as a complex with aminophosphonic acid-guanylate ester and ARA7·VPS9a(D185N) with GDP. Upon complex formation with ARA7, VPS9 wedges into the interswitch region of ARA7, inhibiting the coordination of Mg(2+) and decreasing the stability of GDP binding. The aspartate finger of VPS9a recognizes GDP β-phosphate directly and pulls the P-loop lysine of ARA7 away from GDP β-phosphate toward switch II to further destabilize GDP for its release during the transition from the GDP-bound to nucleotide-free intermediates in the nucleotide exchange reaction.

Neurofibromin: a general outlook

Neurofibromin is a cytoplasmic protein that is predominantly expressed in neurons, Schwann cells, oligodendrocytes, astrocytes and leukocytes. It is encoded by the gene NF1, located on chromosome 17, at q11.2, and has different biochemical functions, including association to microtubules and participation in several signaling pathways. Alterations in this protein are responsible for a phacomatosis named neurofibromatosis type 1.

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

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

Relevance of the kinetic equilibrium of forces to the control of the cell cycle by Ras proteins

In higher organisms, the replacement of GDP bound to Ras proteins with GTP, under the participation of an exchange factor, is an important step in the initiation of cell division. Ras-GTP activates kinases and other effectors, which pass signals to the cell nucleus and to the cytoskeleton. The active state of Ras is terminated by hydrolysis of the bound GTP with the assistance of an activating protein (GAP). Knowledge of these regulatory events is based on extensive experimental data, but many aspects of their interpretation are still controversial. It is assumed here that a significant part of the free energy released when two partners associate is stored in a 'kinetic equilibrium of forces' (KEF), and used to facilitate the separation from a third partner. The activation of the Raf kinase is explained primarily in terms of an allosteric effect of Ras-GTP on the phosphate transfer in the catalytic region of the kinase. A mechanism is proposed for the modification of GAP by Ras-GTP, which is believed to be a prerequisite for the well-known crosstalk between the Ras- and Rho-dependent signalling pathways. The cell, by meeting the requirements for KEF, manages to reduce activation barriers, thus significantly accelerating the regulatory events and other complex biological reaction sequences.

GTP hydrolysis mechanism of Ras-like GTPases

The Ras-like GTPases regulate diverse cellular functions via the chemical cycle of binding and hydrolyzing GTP molecules. They alternate between GTP- and GDP-bound conformations. The GTP-bound conformation is biologically active and promotes a cellular function, such as signal transduction, cytoskeleton organization, protein synthesis/translocation, or a membrane budding/fusion event. GTP hydrolysis turns off the GTPase switch by converting it to the inactive GDP-bound conformation. The fundamental GTP hydrolysis mechanism by these GTPases has generated considerable interest over the last two decades but remained to be firmly established. This review provides an update on the catalytic mechanism with discussions on recent developments from kinetic, structural, and model studies in the context of the various GTP hydrolysis models proposed over the years.

A novel testicular RhoGAP-domain protein induces apoptosis

The GTPase-activating proteins (GAPs) accelerate the hydrolysis of GTP to GDP by small GTPases. The GTPases play diverse roles in many cellular processes, including proliferation, cell motility, endocytosis, nuclear import/export, and nuclear membrane formation. Little is known about GAP-domain proteins in spermatogenesis. We isolated a novel RhoGAP domain-containing tGAP1 protein from male germ cells that exhibits unusual properties. The tGAP1 is expressed at low levels in early spermatogonia. Robust transcription initiates in midpachytene spermatocytes and continues after meiosis. The 175-kDa tGAP1 protein localizes to the cytoplasm of spermatocytes and to the cytoplasm and nucleus in spermatids. The protein contains four GAP domain-related sequences, in contrast to all other GAP proteins that harbor one such domain. No activity toward RhoA, Rac1, or Cdc42 could be detected. Results of transfection studies in various somatic cells indicated that low-level tGAP1 expression significantly slows down the cell cycle. Expression of higher levels of tGAP1 by infection of somatic cells with recombinant adenoviruses demonstrated that tGAP1 efficiently induces apoptosis, which to our knowledge is the first such demonstration for a RhoGAP protein. Based on its subcellular location in spermatids and its activity, tGAP1 may play a role in nuclear import/export.

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.

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

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

Theoretical studies on the hydrolysis of mono-phosphate and tri-phosphate in gas phase and aqueous solution

Phosphate hydrolysis by GTPases plays an important role as a molecular switch in signal transduction and as an initiator of many other biological processes. Despite the centrality of this ubiquitous reaction, the mechanism is still poorly understood. As a first step to understand the mechanisms of this process, the nonenzymatic hydrolysis of mono-phosphate and tri-phosphate esters were systematically studied in gas phase and aqueous solution using hybrid density functional methods. The dielectric effect of the environment on the energetics of these processes was also explored. Theoretical results show that for mono-phosphate ester, the dissociative pathway is much more favorable than the associative pathway. However, the reaction barriers for the dissociative and associative pathways of tri-phosphate hydrolysis are very close in aqueous solution, though the dissociative pathway is more favorable in the gas phase. High dielectric solvents, such as water, significantly lower the activation barrier of the associative pathway due to the greater solvation energy of the associative transition states than that of the reactant complex. By contrast, the barrier of the dissociative pathway, with respect to the gas phase, is less sensitive to the surrounding dielectric. In the associative hydrolysis pathway of the tri-phosphate ester, negative charge is transferred from the gamma-phosphate to beta-phosphate through the bridging ester oxygen and results in Pgamma-O bond dissociation. No analogous charge transfer was observed in the dissociative pathway, where Pgamma-O bond dissociation resulted from proton transfer from the gamma-phosphate to the bridge oxygen. Finally, the active participation of local water molecules can significantly lower the activation energy of the dissociative pathway for both mono-phosphate and tri-phosphate.

Oncogenic mutations and packing defects in protein structure

Oncogenic mutations in expressed proteins are of primary interest to understand tumor formation but their structural consequences bearing on protein function are not clearly understood. In this contribution I report on two illustrative examples, p21ras and p57, revealing that such mutations have an effect on specific structural deficiencies in the packing of the protein structure, i. e., on backbone hydrogen bonds insufficiently shielded from water attack. These structural deficiencies in the wild type are typically "corrected intermolecularly" by protein complexation or protein-ligand association. However, in the oncogenic mutants, these binding signals are partially or completely suppressed: the mutated residues properly wrap or desolvate the hydrogen bonds intramolecularly. Thus, the interactivity of the proteins becomes impaired: their binding affinity decreases sharply, as there is no thermodynamic benefit from removing water surrounding properly desolvated hydrogen bonds. The results, specialized for p21ras and p53, reveal how oncogenic mutations determine a hindrance to GAP-induced hydrolysis (p21) and decrease binding affinity for DNA (p53). Furthermore, the oncogenic potential of mutations in residues not directly engaged in the interface electrostatics is assessed. The results suggest that a high sensitivity of structural defects to genetic accident might be a necessary condition to establish the existence of a proto-oncogene, an angle that merits a systematic study.

Structural fingerprints of the Ras-GTPase activating proteins neurofibromin and p120GAP

Ras specific GTPase activating proteins (GAPs), neurofibromin and p120GAP, bind GTP bound Ras and efficiently complement its active site. Here we present comparative data from mutations and fluorescence-based assays of the catalytic domains of both RasGAPs and interpret them using the crystal structures. Three prominent regions in RasGAPs, the arginine-finger loop, the phenylalanine-leucine-arginine (FLR) region and alpha7/variable loop contain structural fingerprints governing the GAP function. The finger loop is crucial for the stabilization of the transition state of the GTPase reaction. This function is controlled by residues proximal to the catalytic arginine, which are strikingly different between the two RasGAPs. These residues specifically determine the orientation and therefore the positioning of the arginine finger in the Ras active site. The invariant FLR region, a hallmark for RasGAPs, indirectly contributes to GTPase stimulation by forming a scaffold, which stabilizes Ras switch regions. We show that a long hydrophobic side-chain in the FLR region is crucial for this function. The alpha7/variable loop uses several conserved residues including two lysine residues, which are involved in numerous interactions with the switch I region of Ras. This region determines the specificity of the Ras-RasGAP interaction.

Relative molecular similarity in selected chemical carcinogens and the nucleoside triphosphate chain

Several markers of cell toxicity are useful as screening tests for epigenetic carcinogens. The direct effects of chemicals on ATPase and GTPase function are pertinent to the early stages of carcinogenesis. Interference with triphosphate-diphosphate exchange mechanisms may result from the interaction of carcinogens with the substrate triphosphate chain. To investigate this hypothesis, a computational chemistry programme is used in this study to investigate molecular similarity in ATPase inhibitors, carcinogens and tumour promoters, in relation to the nucleoside triphosphate chain. The results show that atoms in the investigated molecular structures superimpose on sets of oxygen atoms in the triphosphate chain with interatomic distances < 0.3A. Relative molecular similarity to the substrate triphosphate chain is discussed in terms of the established inhibitory properties of carcinogens/tumour promoters on ATPase function, the carcinogen/ tumour promoting properties of ATPase inhibitors and the prediction of carcinogenic activity from chemical structure.