Allosteric modulation of Ras-GTP is linked to signal transduction through RAF kinase

Conformations of KRAS4B Affected by Its Partner Binding and G12C Mutation: Insights from GaMD Trajectory-Image Transformation-Based Deep Learning

Binding of partners and mutations highly affects the conformational dynamics of KRAS4B, which is of significance for deeply understanding its function. Gaussian accelerated molecular dynamics (GaMD) simulations followed by deep learning (DL) and principal component analysis (PCA) were carried out to probe the effect of G12C and binding of three partners NF1, RAF1, and SOS1 on the conformation alterations of KRAS4B. DL reveals that G12C and binding of partners result in alterations in the contacts of key structure domains, such as the switch domains SW1 and SW2 together with the loops L4, L5, and P-loop. Binding of NF1, RAF1, and SOS1 constrains the structural fluctuation of SW1, SW2, L4, and L5; on the contrary, G12C leads to the instability of these four structure domains. The analyses of free energy landscapes (FELs) and PCA also show that binding of partners maintains the stability of the conformational states of KRAS4B while G12C induces greater mobility of the switch domains SW1 and SW2, which produces significant impacts on the interactions of GTP with SW1, L4, and L5. Our findings suggest that partner binding and G12C play important roles in the activity and allosteric regulation of KRAS4B, which may theoretically aid in further understanding the function of KRAS4B.

KRAS Mutation Subtypes and Their Association with Other Driver Mutations in Oncogenic Pathways

The KRAS mutation stands out as one of the most influential oncogenic mutations, which directly regulates the hallmark features of cancer and interacts with other cancer-causing driver mutations. However, there remains a lack of precise information on their cooccurrence with mutated variants of KRAS and any correlations between KRAS and other driver mutations. To enquire about this issue, we delved into cBioPortal, TCGA, UALCAN, and Uniport studies. We aimed to unravel the complexity of KRAS and its relationships with other driver mutations. We noticed that G12D and G12V are the prevalent mutated variants of KRAS and coexist with the TP53 mutation in PAAD and CRAD, while G12C and G12V coexist with LUAD. We also noticed similar observations in the case of PIK3CA and APC mutations in CRAD. At the transcript level, a positive correlation exists between KRAS and PIK3CA and between APC and KRAS in CRAD. The existence of the co-mutation of KRAS and other driver mutations could influence the signaling pathway in the neoplastic transformation. Moreover, it has immense prognostic and predictive implications, which could help in better therapeutic management to treat cancer.

Molecular Mechanism of Phosphorylation-Mediated Impacts on the Conformation Dynamics of GTP-Bound KRAS Probed by GaMD Trajectory-Based Deep Learning

The phosphorylation of different sites produces a significant effect on the conformational dynamics of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations were combined with deep learning (DL) to explore the molecular mechanism of the phosphorylation-mediated effect on conformational dynamics of the GTP-bound KRAS. The DL finds that the switch domains are involved in obvious differences in conformation contacts and suggests that the switch domains play a key role in the function of KRAS. The analyses of free energy landscapes (FELs) reveal that the phosphorylation of pY32, pY64, and pY137 leads to more disordered states of the switch domains than the wild-type (WT) KRAS and induces conformational transformations between the closed and open states. The results from principal component analysis (PCA) indicate that principal motions PC1 and PC2 are responsible for the closed and open states of the phosphorylated KRAS. Interaction networks were analyzed and the results verify that the phosphorylation alters interactions of GTP and magnesium ion Mg2+ with the switch domains. It is concluded that the phosphorylation pY32, pY64, and pY137 tune the activity of KRAS through changing conformational dynamics and interactions of the switch domains. We anticipated that this work could provide theoretical aids for deeply understanding the function of KRAS.

Unraveling the Significance of Mg2+ Dependency and Nucleotide Binding Specificity of H-RAS

RAS is a small GTPase and acts as a binary molecular switch; the transition from its active to inactive state plays a crucial role in various cell signaling processes. Molecular dynamics simulations at the atomistic level suggest that the absence of cofactor Mg2+ ion generally leads to pronounced structural changes in the Switch-I than Switch-II regions and assists GTP binding. The presence of the Mg2+ ion also restricts the rotation of ϒ phosphate and enhances the hydrolysis rate of GTP. Further, the simulations reveal that the stability of the protein is almost uncompromised when Mg2+ is replaced with Zn2+ and not the Ca2+ ion. The specificity of H-RAS to GTP was evaluated by substituting with ATP and CTP, which indicates that the binding pocket tolerates purine bases over pyrimidine bases. However, the D119 residue specifically interacts with the guanine base and serves as one of the primary interactions that leads to the selectivity of GTP over ATP. The ring displacement of 32Y serves as gate dynamics in H-RAS which are important for its interaction with GAP for the nucleotide exchange and is restricted in the presence of ATP. Finally, the point mutations 61, 16, and 32 influence the structural changes, specifically in the Switch-II region, which are expected to impact the GTP hydrolysis and thus are termed oncogenic mutations.

Robust Strategy for Hit-to-Lead Discovery: NMR for SAR

Establishing robust structure-activity relationships (SARs) is key to successful drug discovery campaigns, yet it often remains elusive due to screening and hit validation artifacts (false positives and false negatives), which frequently result in unproductive downstream expenditures of time and resources. To address this issue, we developed an integrative biophysics-driven strategy that expedites hit-to-lead discovery, mitigates false positives/negatives and common hit validation errors, and provides a robust approach to obtaining accurate binding and affinity measurements. The advantage of this method is that it vastly improves the clarity and reproducibility for affinity-driven SAR by monitoring and eliminating confounding factors. We demonstrate the ease at which high-quality micromolar binders can be generated from the initial millimolar fragment screening hits against an "undruggable" protein target, HRas.

Kinetic and thermodynamic allostery in the Ras protein family

Allostery, the transfer of information between distant parts of a macromolecule, is a fundamental feature of protein function and regulation. However, allosteric mechanisms are usually not explained by protein structure, requiring information on correlated fluctuations uniquely accessible to molecular simulation. Existing work to extract allosteric pathways from molecular dynamics simulations has focused on thermodynamic correlations. Here, we show how kinetic correlations encode complementary information essential to explain observed variations in allosteric regulation. We applied kinetic and thermodynamic correlation analysis on atomistic simulations of H, K, and NRas isoforms in the apo, GTP, and GDP-bound states of Ras protein, with and without complexing to its downstream effector, Raf. We show that switch I and switch II are the primary components of thermodynamic and kinetic allosteric networks, consistent with the key roles of these two motifs. These networks connect the switches to an allosteric loop recently discovered from a crystal structure of HRas. This allosteric loop is inactive in KRas, but is coupled to the hydrolysis arm switch II in NRas and HRas. We find that the mechanism in the latter two isoforms are thermodynamic and kinetic, respectively. Binding of Raf-RBD further activates thermodynamic allostery in HRas and KRas but has limited effect on NRas. These results indicate that kinetic and thermodynamic correlations are both needed to explain protein function and allostery. These two distinct channels of allosteric regulation, and their combinatorial variability, may explain how subtle mutational differences can lead to diverse regulatory profiles among enzymatic proteins.

Allosteric site variants affect GTP hydrolysis on Ras

RAS GTPases are proto-oncoproteins that regulate cell growth, proliferation, and differentiation in response to extracellular signals. The signaling functions of RAS, and other small GTPases, are dependent on their ability to cycle between GDP-bound and GTP-bound states. Structural analyses suggest that GTP hydrolysis catalyzed by HRAS can be regulated by an allosteric site located between helices 3, 4, and loop 7. Here we explore the relationship between intrinsic GTP hydrolysis on HRAS and the position of helix 3 and loop 7 through manipulation of the allosteric site, showing that the two sites are functionally connected. We generated several hydrophobic mutations in the allosteric site of HRAS to promote shifts in helix 3 relative to helix 4. By combining crystallography and enzymology to study these mutants, we show that closure of the allosteric site correlates with increased hydrolysis of GTP on HRAS in solution. Interestingly, binding to the RAS binding domain of RAF kinase (RAF-RBD) inhibits GTP hydrolysis in the mutants. This behavior may be representative of a cluster of mutations found in human tumors, which potentially cooperate with RAF complex formation to stabilize the GTP-bound state of RAS.

RAS-targeted cancer therapy: Advances in drugging specific mutations

Rat sarcoma (RAS), as a frequently mutated oncogene, has been studied as an attractive target for treating RAS-driven cancers for over four decades. However, it is until the recent success of kirsten-RAS (KRAS)^G12C inhibitor that RAS gets rid of the title "undruggable". It is worth noting that the therapeutic effect of KRAS^G12C inhibitors on different RAS allelic mutations or even different cancers with KRAS^G12C varies significantly. Thus, deep understanding of the characteristics of each allelic RAS mutation will be a prerequisite for developing new RAS inhibitors. In this review, the structural and biochemical features of different RAS mutations are summarized and compared. Besides, the pathological characteristics and treatment responses of different cancers carrying RAS mutations are listed based on clinical reports. In addition, the development of RAS inhibitors, either direct or indirect, that target the downstream components in RAS pathway is summarized as well. Hopefully, this review will broaden our knowledge on RAS-targeting strategies and trigger more intensive studies on exploiting new RAS allele-specific inhibitors.

Path Sampling Simulations Reveal How the Q61L Mutation Alters the Dynamics of KRas

Flexibility is essential for many proteins to function, but can be difficult to characterize. Experiments lack resolution in space and time, while the time scales involved are prohibitively long for straightforward molecular dynamics simulations. In this work, we present a multiple state transition path sampling simulation study of a protein that has been notoriously difficult to characterize in its active state. The GTPase enzyme KRas is a signal transduction protein in pathways for cell differentiation, growth, and division. When active, KRas tightly binds guanosine triphosphate (GTP) in a rigid state. The protein-GTP complex can also visit more flexible states, in which it is not active. KRas mutations can affect the conversion between these rigid and flexible states, thus prolonging the activation of signal transduction pathways, which may result in tumor formation. In this work, we apply path sampling simulations to investigate the dynamic behavior of KRas-4B (wild type, WT) and the oncogenic mutant Q61L (Q61L). Our results show that KRas visits several conformational states, which are the same for WT and Q61L. The multiple state transition path sampling (MSTPS) method samples transitions between the different states in a single calculation. Tracking which transitions occur shows large differences between WT and Q61L. The MSTPS results further reveal that for Q61L, a route to a more flexible state is inaccessible, thus shifting the equilibrium to more rigid states. The methodology presented here enables a detailed characterization of protein flexibility on time scales not accessible with brute-force molecular dynamics simulations.

Delineating The RAS Conformational Landscape

Mutations in RAS isoforms (KRAS, NRAS, and HRAS) are among the most frequent oncogenic alterations in many cancers, making these proteins high priority therapeutic targets. Effectively targeting RAS isoforms requires an exact understanding of their active, inactive, and druggable conformations. However, there is no structural catalog of RAS conformations to guide therapeutic targeting or examining the structural impact of RAS mutations. Here we present an expanded classification of RAS conformations based on analyses of the catalytic switch 1 (SW1) and switch 2 (SW2) loops. From 721 human KRAS, NRAS, and HRAS structures available in the Protein Data Bank (206 RAS-protein co-complexes, 190 inhibitor-bound, and 325 unbound, including 204 WT and 517 mutated structures), we created a broad conformational classification based on the spatial positions of Y32 in SW1 and Y71 in SW2. Clustering all well-modeled SW1 and SW2 loops using a density-based machine learning algorithm defined additional conformational subsets, some previously undescribed. Three SW1 conformations and nine SW2 conformations were identified, each associated with different nucleotide states (GTP-bound, nucleotide-free, and GDP-bound) and specific bound proteins or inhibitor sites. The GTP-bound SW1 conformation could be further subdivided based on the hydrogen bond type made between Y32 and the GTP γ-phosphate. Further analysis clarified the catalytic impact of G12D and G12V mutations and the inhibitor chemistries that bind to each druggable RAS conformation. Overall, this study has expanded our understanding of RAS structural biology, which could facilitate future RAS drug discovery.

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.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

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.

Equilibria between conformational states of the Ras oncogene protein revealed by high pressure crystallography

In this work, we experimentally investigate the allosteric transitions between conformational states on the Ras oncogene protein using high pressure crystallography. Ras protein is a small GTPase involved in central regulatory processes occurring in multiple conformational states. Ras acts as a molecular switch between active GTP-bound, and inactive GDP-bound states, controlling essential signal transduction pathways. An allosteric network of interactions between the effector binding regions and the membrane interacting regions is involved in Ras cycling. The conformational states which coexist simultaneously in solution possess higher Gibbs free energy than the ground state. Equilibria between these states can be shifted by applying pressure favouring conformations with lower partial molar volume, and has been previously analyzed by high-pressure NMR spectroscopy. High-pressure macromolecular crystallography (HPMX) is a powerful tool perfectly complementary to high-pressure NMR, allowing characterization at the molecular level with a high resolution the different allosteric states involved in the Ras cycling. We observe a transition above 300 MPa in the crystal leading to more stable conformers. Thus, we compare the crystallographic structures of Ras(wt)·Mg²⁺·GppNHp and Ras(D33K)·Mg²⁺·GppNHp at various high hydrostatic pressures. This gives insight into per-residue descriptions of the structural plasticity involved in allosteric equilibria between conformers. We have mapped out at atomic resolution the different segments of Ras protein which remain in the ground-state conformation or undergo structural changes, adopting excited-energy conformations corresponding to transient intermediate states. Such in crystallo phase transitions induced by pressure open the possibility to finely explore the structural determinants related to switching between Ras allosteric sub-states without any mutations nor exogenous partners.

The equilibria between structural states induced by pressure within the crystal structure of Ras are illustrated with different colors corresponding to different Ras substates.

Conformational Fluctuations in GTP-Bound K-Ras: A Metadynamics Perspective with Harmonic Linear Discriminant Analysis

Biomacromolecules often undergo significant conformational rearrangements during function. In proteins, these motions typically consist in nontrivial, concerted rearrangement of multiple flexible regions. Mechanistic, thermodynamics, and kinetic predictions can be obtained via molecular dynamics simulations, provided that the simulation time is at least comparable to the relevant time scale of the process of interest. Because of the substantial computational cost, however, plain MD simulations often have difficulty in obtaining sufficient statistics for converged estimates, requiring the use of more-advanced techniques. Central in many enhanced sampling methods is the definition of a small set of relevant degrees of freedom (collective variables) that are able to describe the transitions between different metastable states of the system. The harmonic linear discriminant analysis (HLDA) has been shown to be useful for constructing low-dimensional collective variables in various complex systems. Here, we apply HLDA to study the free-energy landscape of a monomeric protein around its native state. More precisely, we study the K-Ras protein bound to GTP, focusing on two flexible loops and on the region associated with oncogenic mutations. We perform microsecond-long biased simulations on the wild type and on G12C, G12D, G12 V mutants, describe the resulting free-energy landscapes, and compare our predictions with previous experimental and computational studies. The fast interconversion between open and closed macroscopic states and their similar thermodynamic stabilities are observed. The mutation-induced effects include the alternations of the relative stabilities of different conformational states and the introduction of many microscopic metastable states. Together, our results demonstrate the applicability of the HLDA-based protocol for the conformational sampling of multiple flexible regions in folded proteins.

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.

KRAS G12C Mutations in NSCLC: From Target to Resistance

Lung cancer represents the most common form of cancer, accounting for 1.8 million deaths globally in 2020. Over the last decade the treatment for advanced and metastatic non-small cell lung cancer have dramatically improved largely thanks to the emergence of two therapeutic breakthroughs: the discovery of immune checkpoint inhibitors and targeting of oncogenic driver alterations. While these therapies hold great promise, they face the same limitation as other inhibitors: the emergence of resistant mechanisms. One such alteration in non-small cell lung cancer is the Kirsten Rat Sarcoma (KRAS) oncogene. KRAS mutations are the most common oncogenic driver in NSCLC, representing roughly 20-25% of cases. The mutation is almost exclusively detected in adenocarcinoma and is found among smokers 90% of the time. Along with the development of new drugs that have been showing promising activity, resistance mechanisms have begun to be clarified. The aim of this review is to unwrap the biology of KRAS in NSCLC with a specific focus on primary and secondary resistance mechanisms and their possible clinical implications.

Mutation-Induced Impacts on the Switch Transformations of the GDP- and GTP-Bound K-Ras: Insights from Multiple Replica Gaussian Accelerated Molecular Dynamics and Free Energy Analysis

Mutations yield significant effect on the structural flexibility of two switch domains, SW1 and SW2, in K-Ras, which is considered as an important target of anticancer drug design. To unveil a molecular mechanism with regard to mutation-mediated tuning on the activity of K-Ras, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations followed by analysis of free energy landscapes (FELs) are performed on the GDP- and GTP-bound wild-type (WT), G12V, and D33E K-Ras. The results suggest that G12V and D33E not only evidently change the flexibility of SW1 and SW2 but also greatly affect correlated motions of SW1 and SW2 separately relative to the P-loop and SW1, which exerts a certain tuning on the activity of K-Ras. The information stemming from the analyses of FELs reveals that the conformations of SW1 and SW2 are in high disorders in the GDP- and GTP-associated WT and mutated K-Ras, possibly producing significant effect on binding of guanine nucleotide exchange factors or effectors to K-Ras. The interaction networks of GDP and GTP with K-Ras are identified and the results uncover that the instability in hydrogen-bonding interactions of SW1 with GDP and GTP is mostly responsible for conformational disorder of SW1 and SW2 as well as tunes the activity of oncogenic K-Ras.

Engineering subtilisin proteases that specifically degrade active RAS

We describe the design, kinetic properties, and structures of engineered subtilisin proteases that degrade the active form of RAS by cleaving a conserved sequence in switch 2. RAS is a signaling protein that, when mutated, drives a third of human cancers. To generate high specificity for the RAS target sequence, the active site was modified to be dependent on a cofactor (imidazole or nitrite) and protease sub-sites were engineered to create a linkage between substrate and cofactor binding. Selective proteolysis of active RAS arises from a 2-step process wherein sub-site interactions promote productive binding of the cofactor, enabling cleavage. Proteases engineered in this way specifically cleave active RAS in vitro, deplete the level of RAS in a bacterial reporter system, and also degrade RAS in human cell culture. Although these proteases target active RAS, the underlying design principles are fundamental and will be adaptable to many target proteins.

How to make an undruggable enzyme druggable: lessons from ras proteins

Significant advances have been made toward discovering allosteric inhibitors for challenging drug targets such as the Ras family of membrane-associated signaling proteins. Malfunction of Ras proteins due to somatic mutations is associated with up to a quarter of all human cancers. Computational techniques have played critical roles in identifying and characterizing allosteric ligand-binding sites on these proteins, and to screen ligand libraries against those sites. These efforts, combined with a wide range of biophysical, structural, biochemical and cell biological experiments, are beginning to yield promising inhibitors to treat malignancies associated with mutated Ras proteins. In this chapter, we discuss some of these developments and how the lessons learned from Ras might be applied to similar other challenging drug targets.

Isoform-Specific Destabilization of the Active Site Reveals a Molecular Mechanism of Intrinsic Activation of KRas G13D

Ras GTPases are mutated at codons 12, 13, and 61, with different frequencies in KRas, HRas, and NRas and in a cancer-specific manner. The G13D mutant appears in 25% of KRas-driven colorectal cancers, while observed only rarely in HRas or NRas. Structures of Ras G13D in the three isoforms show an open active site, with adjustments to the D13 backbone torsion angles and with disconnected switch regions. KRas G13D has unique features that destabilize the nucleotide-binding pocket. In KRas G13D bound to GDP, A59 is placed in the Mg²⁺ binding site, as in the HRas-SOS complex. Structure and biochemistry are consistent with an intermediate level of KRas G13D bound to GTP, relative to wild-type and KRas G12D, observed in genetically engineered mouse models. The results explain in part the elevated frequency of the G13D mutant in KRas over the other isoforms of Ras.

The Frequency of Ras Mutations in Cancer

Ras is frequently mutated in cancer, however, there is a lack of consensus in the literature regarding the cancer mutation frequency of Ras, with quoted values varying from 10%-30%. This variability is at least in part due to the selective aggregation of data from different databases and the dominant influence of particular cancer types and particular Ras isoforms within these datasets. To provide a more definitive figure for Ras mutation frequency in cancer, we cross-referenced the data in all major publicly accessible cancer mutation databases to determine reliable mutation frequency values for each Ras isoform in all major cancer types. These percentages were then applied to current U.S. cancer incidence statistics to estimate the number of new patients each year that have Ras-mutant cancers. We find that approximately 19% of patients with cancer harbor Ras mutations, equivalent to approximately 3.4 million new cases per year worldwide. We discuss the Ras isoform and mutation-specific trends evident within the datasets that are relevant to current Ras-targeted therapies.

Nerve growth factor in metabolic complications and Alzheimer's disease: Physiology and therapeutic potential

As the population ages, obesity and metabolic complications as well as neurological disorders are becoming more prevalent, with huge economic burdens on both societies and families. New therapeutics are urgently needed. Nerve growth factor (NGF), first discovered in 1950s, is a neurotrophic factor involved in regulating cell proliferation, growth, survival, and apoptosis in both central and peripheral nervous systems. NGF and its precursor, proNGF, bind to TrkA and p75 receptors and initiate protein phosphorylation cascades, resulting in changes of cellular functions, and are associated with obesity, diabetes and its complications, and Alzheimer's disease. In this article, we summarize changes in NGF levels in metabolic and neuronal disorders, the signal transduction initiated by NGF and proNGF, the physiological and pathophysiological relevance, and therapeutic potential in treating chronic metabolic diseases and cognitive decline.

Distinct Binding Preferences between Ras and Raf Family Members and the Impact on Oncogenic Ras Signaling

The Ras GTPases are frequently mutated in human cancer, and, although the Raf kinases are essential effectors of Ras signaling, the tumorigenic properties of specific Ras-Raf complexes are not well characterized. Here, we examine the ability of individual Ras and Raf proteins to interact in live cells using bioluminescence resonance energy transfer (BRET) technology. We find that C-Raf binds all mutant Ras proteins with high affinity, whereas B-Raf exhibits a striking preference for mutant K-Ras. This selectivity is mediated by the acidic, N-terminal segment of B-Raf and requires the K-Ras polybasic region for high-affinity binding. In addition, we find that C-Raf is critical for mutant H-Ras-driven signaling and that events stabilizing B-Raf/C-Raf dimerization, such as Raf inhibitor treatment or certain B-Raf mutations, can allow mutant H-Ras to engage B-Raf with increased affinity to promote tumorigenesis, thus revealing a previously unappreciated role for C-Raf in potentiating B-Raf function.

Cosolvent Analysis Toolkit (CAT): a robust hotspot identification platform for cosolvent simulations of proteins to expand the druggable proteome

Cosolvent Molecular Dynamics (MD) simulations are increasingly popular techniques developed for prediction and characterization of allosteric and cryptic binding sites, which can be rendered “druggable” by small molecule ligands. Despite their conceptual simplicity and effectiveness, the analysis of cosolvent MD trajectories relies on pocket volume data, which requires a high level of manual investigation and may introduce a bias. In this work, we present CAT (Cosolvent Analysis Toolkit): an open-source, freely accessible analytical tool, suitable for automated analysis of cosolvent MD trajectories. CAT is compatible with commonly used molecular graphics software packages such as UCSF Chimera and VMD. Using a novel hybrid empirical force field scoring function, CAT accurately ranks the dynamic interactions between the macromolecular target and cosolvent molecules. To benchmark, CAT was used for three validated protein targets with allosteric and orthosteric binding sites, using five chemically distinct cosolvent molecules. For all systems, CAT has accurately identified all known sites. CAT can thus assist in computational studies aiming at identification of protein “hotspots” in a wide range of systems. As an easy-to-use computational tool, we expect that CAT will contribute to an increase in the size of the potentially ‘druggable’ human proteome.

Ras-ERK1/2 Signaling Promotes The Development Of Osteosarcoma By Regulating H2BK12ac Through CBP

Background

H2BK12ac is an important histone acetylation pattern of H2B, which has been reported in several cancers. However, whether H2BK12ac joins in Ras-ERK1/2 activation-induced osteosarcoma (OS) cell behaviors remain unclear. The study explored this peradventure and revealed the underlying mechanism.

Methods

MG-63 cells were transfected with pEGFP-N1, pEGFP-Ras^WT and pEGFP-K-Ras^G12V/T35S, H2BK12ac and ERK1/2 expression levels were analyzed by Western blot. Effects of H2BK12ac on cell viability, migration, colony formation and cell cycle were investigated by MTT, Transwell, soft-agar colony formation and flow cytometry assays. RT-qPCR and ChIP were performed to study the effect of H2BK12ac and CBP on ERK1/2-downstream gene transcriptions.

Results

H2BK12ac was specifically down-regulated by Ras-ERK1/2 activation in MG-63 cells. Down-regulated H2BK12ac participated in regulating cell proliferation and migration of MG-63 cells, meanwhile, affected the transcription of ERK1/2-downstream genes. Additionally, silence of HDAC1 up-regulated H2BK12ac expression, and inhibited the promoting effect of Ras-ERK1/2 on MG-63 cells' proliferation, migration and RNA expression levels of ERK1/2-downstream genes. Further, the degradation of CBP mediated by MDM2 was discovered to be linked to Ras-ERK1/2 activation-induced H2BK12ac down-regulation.

Conclusion

These findings from the study demonstrated that Ras-ERK1/2 signaling could promote the development of OS via regulating H2BK12ac through MDM2-mediated CBP degradation.

Retracted: Ras-ERK1/2 signaling promotes the development of uveal melanoma by downregulating H3K14ac

Ras-extracellular signal-regulated protein kinases 1 and 2 (ERK1/2) signaling has been proposed as the crucial regulators in the development of various cancers. Histone acetylation at H3 lysine 14 (H3K14ac) is closely associated with gene expression and DNA damage. However, whether H3K14ac participates in mediating Ras-ERK1/2-induced cell proliferation and migration in uveal melanoma cells remains unknown. The purpose of this study is to investigate the effect of H3K14ac on Ras-ERK1/2 affected uveal melanoma cell phenotypes. MP65 cells were transfected with Ras ^WT and Ras ^G12V/T35S , the unloaded plasmid of pEGFP-N1 served as a negative control. Protein levels of phosphorylated ERK1/2 ^Thr202 and H3K14ac were assessed by western blot assay. Cell viability, number of colonies, migration, and the downstream genes of ERK1/2 were analyzed by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2-H-tetrazolium bromide, soft-agar colony formation, transwell, and chromatin immunoprecipitation assays. HA-tag vectors of CLR3 and TIP60 and the small interfering RNAs that specific for CLR3 and MDM2 were transfected into MP65 cells to uncover the effects of CLR3, TIP60, and MDM2 on Ras-ERK1/2 mediated H3K14ac expression and MP65 cell phenotypes. We found that, Ras-ERK1/2 decreased H3K14ac expression in MP65 cells, and H3K14ac significantly suppressed Ras-ERK1/2-induced cell viability, colony formation, and migration in MP65 cells. Moreover, the transcription of CYR61, IGFBP3, WNT16B, NT5E, GDF15, and CARD16 was regulated by H3K14ac. Additionally, CLR3 silence recovered H3K14ac expression and reversed the effect of Ras-ERK1/2 on MP65 cell proliferation, migration and the mRNAs of ERK1/2 downstream genes. Besides, Ras-ERK1/2 decreased H3K14ac expression by MDM2-mediated TIP60 degradation. In conclusion, Ras-ERK1/2 promoted uveal melanoma cells growth and migration by downregulating H3K14ac via MDM2-mediated TIP60 degradation.

Oncogenic G12D mutation alters local conformations and dynamics of K-Ras

K-Ras is the most frequently mutated oncoprotein in human cancers, and G12D is its most prevalent mutation. To understand how G12D mutation impacts K-Ras function, we need to understand how it alters the regulation of its dynamics. Here, we present local changes in K-Ras structure, conformation and dynamics upon G12D mutation, from long-timescale Molecular Dynamics simulations of active (GTP-bound) and inactive (GDP-bound) forms of wild-type and mutant K-Ras, with an integrated investigation of atomistic-level changes, local conformational shifts and correlated residue motions. Our results reveal that the local changes in K-Ras are specific to bound nucleotide (GTP or GDP), and we provide a structural basis for this. Specifically, we show that G12D mutation causes a shift in the population of local conformational states of K-Ras, especially in Switch-II (SII) and α3-helix regions, in favor of a conformation that is associated with a catalytically impaired state through structural changes; it also causes SII motions to anti-correlate with other regions. This detailed picture of G12D mutation effects on the local dynamic characteristics of both active and inactive protein helps enhance our understanding of local K-Ras dynamics, and can inform studies on the development of direct inhibitors towards the treatment of K-Ras^G12D-driven cancers.

Uncommon mutational profiles of metastatic colorectal cancer detected during routine genotyping using next generation sequencing

RAS genotyping is mandatory to predict anti-EGFR monoclonal antibodies (mAbs) therapy resistance and BRAF genotyping is a relevant prognosis marker in patients with metastatic colorectal cancer. Although the role of hotspot mutations is well defined, the impact of uncommon mutations is still unknown. In this study, we aimed to discuss the potential utility of detecting uncommon RAS and BRAF mutation profiles with next-generation sequencing. A total of 779 FFPE samples from patients with metastatic colorectal cancer with valid NGS results were screened and 22 uncommon mutational profiles of KRAS, NRAS and BRAF genes were selected. In silico prediction of mutation impact was then assessed by 2 predictive scores and a structural protein modelling. Three samples carry a single KRAS non-hotspot mutation, one a single NRAS non-hotspot mutation, four a single BRAF non-hotspot mutation and fourteen carry several mutations. This in silico study shows that some non-hotspot RAS mutations seem to behave like hotspot mutations and warrant further examination to assess whether they should confer a resistance to anti-EGFR mAbs therapy for patients bearing these non-hotspot RAS mutations. For BRAF gene, non-V600E mutations may characterise a novel subtype of mCRC with better prognosis, potentially implying a modification of therapeutic strategy.

Advances in NMR Methods to Identify Allosteric Sites and Allosteric Ligands

NMR allows assessment of protein structure in solution. Unlike conventional X-ray crystallography that provides snapshots of protein conformations, all conformational states are simultaneously accessible to analysis by NMR. This is a significant advantage for discovery and characterization of allosteric effects. These effects are observed when binding at one site of the protein affects another distinct site through conformational transitions. Allosteric regulation of proteins has been observed in multiple physiological processes in health and disease, providing an opportunity for the development of allosteric inhibitors. These compounds do not directly interact with the orthosteric site of the protein but influence its structure and function. In this book chapter, we provide an overview on how NMR methods are utilized to identify allosteric sites and to discover novel inhibitors, highlighting examples from the field. We also describe how NMR has contributed to understanding of allosteric mechanisms and propose that it is likely to play an important role in clarification and further development of key concepts of allostery.

Structure-Guided Protein Transition Modeling with a Probabilistic Roadmap Algorithm

Proteins are macromolecules in perpetual motion, switching between structural states to modulate their function. A detailed characterization of the precise yet complex relationship between protein structure, dynamics, and function requires elucidating transitions between functionally-relevant states. Doing so challenges both wet and dry laboratories, as protein dynamics involves disparate temporal scales. In this paper, we present a novel, sampling-based algorithm to compute transition paths. The algorithm exploits two main ideas. First, it leverages known structures to initialize its search and define a reduced conformation space for rapid sampling. This is key to address the insufficient sampling issue suffered by sampling-based algorithms. Second, the algorithm embeds samples in a nearest-neighbor graph where transition paths can be efficiently computed via queries. The algorithm adapts the probabilistic roadmap framework that is popular in robot motion planning. In addition to efficiently computing lowest-cost paths between any given structures, the algorithm allows investigating hypotheses regarding the order of experimentally-known structures in a transition event. This novel contribution is likely to open up new venues of research. Detailed analysis is presented on multiple-basin proteins of relevance to human disease. Multiscaling and the AMBER ff14SB force field are used to obtain energetically-credible paths at atomistic detail.

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.

Pharmacological Induction of RAS-GTP Confers RAF Inhibitor Sensitivity in KRAS Mutant Tumors

Targeting KRAS mutant tumors through inhibition of individual downstream pathways has had limited clinical success. Here we report that RAF inhibitors exhibit little efficacy in KRAS mutant tumors. In combination drug screens, MEK and PI3K inhibitors synergized with pan-RAF inhibitors through an RAS-GTP-dependent mechanism. Broad cell line profiling with RAF/MEK inhibitor combinations revealed synergistic efficacy in KRAS mutant and wild-type tumors, with KRAS^G13D mutants exhibiting greater synergy versus KRAS^G12 mutant tumors. Mechanistic studies demonstrate that MEK inhibition induced RAS-GTP levels, RAF dimerization and RAF kinase activity resulting in MEK phosphorylation in synergistic tumor lines regardless of KRAS status. Taken together, our studies uncover a strategy to rewire KRAS mutant tumors to confer sensitivity to RAF kinase inhibition.

The K-Ras, N-Ras, and H-Ras Isoforms: Unique Conformational Preferences and Implications for Targeting Oncogenic Mutants

Ras controls a multitude of cellular signaling processes, including cell proliferation, differentiation, and apoptosis. Deregulation of Ras cycling often promotes tumorigenesis and various other developmental disorders, termed RASopothies. Although the structure of Ras has been known for many decades, it is still one of the most highly sought-after drug targets today, and is often referred to as "undruggable." At the center of this paradoxical protein is a lack of understanding of fundamental differences in the G domains between the highly similar Ras isoforms and common oncogenic mutations, despite the immense wealth of knowledge accumulated about this protein to date. A shift in the field during the past few years toward a high-resolution understanding of the structure confirms the hypothesis that each isoform and oncogenic mutation must be considered individually, and that not all Ras mutations are created equal. For the first time in Ras history, we have the ability to directly compare the structures of each wild-type isoform to construct a "base-line" understanding, which can then be used as a springboard for analyzing the effects of oncogenic mutations on the structure-function relationship in Ras. This is a fundamental and large step toward the goal of developing personalized therapies for patients with Ras-driven cancers and diseases.

Structural insight into the rearrangement of the switch I region in GTP-bound G12A K-Ras

K-Ras, a molecular switch that regulates cell growth, apoptosis and metabolism, is activated when it undergoes a conformation change upon binding GTP and is deactivated following the hydrolysis of GTP to GDP. Hydrolysis of GTP in water is accelerated by coordination to K-Ras, where GTP adopts a high-energy conformation approaching the transition state. The G12A mutation reduces intrinsic K-Ras GTP hydrolysis by an unexplained mechanism. Here, crystal structures of G12A K-Ras in complex with GDP, GTP, GTPγS and GppNHp, and of Q61A K-Ras in complex with GDP, are reported. In the G12A K-Ras-GTP complex, the switch I region undergoes a significant reorganization such that the Tyr32 side chain points towards the GTP-binding pocket and forms a hydrogen bond to the GTP γ-phosphate, effectively stabilizing GTP in its precatalytic state, increasing the activation energy required to reach the transition state and contributing to the reduced intrinsic GTPase activity of G12A K-Ras mutants.

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.

Molecular Dynamics Simulations and Dynamic Network Analysis Reveal the Allosteric Unbinding of Monobody to H-Ras Triggered by R135K Mutation

Ras proteins, as small GTPases, mediate cell proliferation, survival and differentiation. Ras mutations have been associated with a broad spectrum of human cancers and thus targeting Ras represents a potential way forward for cancer therapy. A recently reported monobody NS1 allosterically disrupts the Ras-mediated signaling pathway, but its efficacy is reduced by R135K mutation in H-Ras. However, the detailed mechanism is unresolved. Here, using molecular dynamics (MD) simulations and dynamic network analysis, we explored the molecular mechanism for the unbinding of NS1 to H-Ras and shed light on the underlying allosteric network in H-Ras. MD simulations revealed that the overall structures of the two complexes did not change significantly, but the H-Ras–NS1 interface underwent significant conformational alteration in the mutant Binding free energy analysis showed that NS1 binding was unfavored after R135K mutation, which resulted in the unfavorable binding of NS1. Furthermore, the critical residues on H-Ras responsible for the loss of binding of NS1 were identified. Importantly, the allosteric networks for these important residues were revealed, which yielded a novel insight into the allosteric regulatory mechanism of H-Ras.

The small GTPases K-Ras, N-Ras, and H-Ras have distinct biochemical properties determined by allosteric effects

H-Ras, K-Ras, and N-Ras are small GTPases that are important in the control of cell proliferation, differentiation, and survival, and their mutants occur frequently in human cancers. The G-domain, which catalyzes GTP hydrolysis and mediates downstream signaling, is 95% conserved between the Ras isoforms. Because of their very high sequence identity, biochemical studies done on H-Ras have been considered representative of all three Ras proteins. We show here that this is not a valid assumption. Using enzyme kinetic assays under identical conditions, we observed clear differences between the three isoforms in intrinsic catalysis of GTP by Ras in the absence and presence of the Ras-binding domain (RBD) of the c-Raf kinase protein (Raf-RBD). Given their identical active sites, isoform G-domain differences must be allosteric in origin, due to remote isoform-specific residues that affect conformational states. We present the crystal structure of N-Ras bound to a GTP analogue and interpret the kinetic data in terms of structural features specific for H-, K-, and N-Ras.

A Novel HRAS Mutation Independently Contributes to Left Ventricular Hypertrophy in a Family with a Known MYH7 Mutation

Several genetic conditions can lead to left ventricular hypertrophy (LVH). Among them, hypertrophic cardiomyopathy (HCM), caused by mutations in sarcomere genes, is the most common inherited cardiac disease. Instead, RASopathies, a rare class of disorders characterized by neuro-cardio-facial-cutaneous abnormalities and sometimes presenting with LVH, are caused by mutations in the RAS-MAPK pathway. We report on a 62-years-old male who presented isolated severe obstructive LVH but did not carry the sarcomere mutation previously identified in his affected relatives. By exome sequencing, we detected a novel mutation in HRAS gene (NM_005343.2:p.Arg68Trp), present also in the proband's daughter, who showed mild LVH and severe intellectual disability. The cardiac phenotype was indistinguishable between family members carrying either mutation. In silico studies suggested that the mutated HRAS protein is constitutionally activated. Consistently, functional characterization in vitro confirmed elevated HRAS-GTP accumulation and downstream RAS-MAPK pathway activation that are known to drive cell proliferation in LVH. Our study emphasizes the role of RAS signaling in cardiac hypertrophy and highlights the complexity in differential diagnosis of RASopathies. In fact, the mild features of RASopathy and the recurrence of sarcomeric HCM in this family delayed the correct diagnosis until comprehensive genetic testing was performed.

Drugging Ras GTPase: a comprehensive mechanistic and signaling structural view

Ras proteins are small GTPases, cycling between inactive GDP-bound and active GTP-bound states. Through these switches they regulate signaling that controls cell growth and proliferation. Activating Ras mutations are associated with approximately 30% of human cancers, which are frequently resistant to standard therapies. Over the past few years, structural biology and in silico drug design, coupled with improved screening technology, led to a handful of promising inhibitors, raising the possibility of drugging Ras proteins. At the same time, the invariable emergence of drug resistance argues for the critical importance of additionally honing in on signaling pathways which are likely to be involved. Here we overview current advances in Ras structural knowledge, including the conformational dynamic of full-length Ras in solution and at the membrane, therapeutic inhibition of Ras activity by targeting its active site, allosteric sites, and Ras-effector protein-protein interfaces, Ras dimers, the K-Ras4B/calmodulin/PI3Kα trimer, and targeting Ras with siRNA. To mitigate drug resistance, we propose signaling pathways that can be co-targeted along with Ras and explain why. These include pathways leading to the expression (or activation) of YAP1 and c-Myc. We postulate that these and Ras signaling pathways, MAPK/ERK and PI3K/Akt/mTOR, act independently and in corresponding ways in cell cycle control. The structural data are instrumental in the discovery and development of Ras inhibitors for treating RAS-driven cancers. Together with the signaling blueprints through which drug resistance can evolve, this review provides a comprehensive and innovative master plan for tackling mutant Ras proteins.

Phosphoproteomic Analyses of NRAS(G12) and NRAS(Q61) Mutant Melanocytes Reveal Increased CK2α Kinase Levels in NRAS(Q61) Mutant Cells

In melanoma, mutant and thereby constantly active neuroblastoma rat sarcoma (NRAS) affects 15-20% of tumors, contributing to tumor initiation, growth, invasion, and metastasis. Recent therapeutic approaches aim to mimic RAS extinction by interfering with critical signaling pathways downstream of the mutant protein. This study investigates the phosphoproteome of primary human melanocytes bearing mutations in the two hot spots of NRAS, NRAS(G12) and NRAS(Q61). Stable isotope labeling by amino acids in cell culture followed by mass spectrometry identified 14,155 spectra of 3,371 unique phosphopeptides mapping to 1,159 proteins (false discovery rate < 2%). Data revealed pronounced PI3K/AKT signaling in NRAS(G12V) mutant cells and pronounced mitogen-activated protein kinase (MAPK) signaling in NRAS(Q61L) variants. Computer-based prediction models for kinases involved, revealed that CK2α is significantly overrepresented in primary human melanocytes bearing NRAS(Q61L) mutations. Similar differences were found in human NRAS(Q61) mutant melanoma cell lines that were also more sensitive to pharmacologic CK2α inhibition compared with NRAS(G12) mutant cells. Furthermore, CK2α levels were pronounced in patient samples of NRAS(Q61) mutant melanoma at the mRNA and protein level. The preclinical findings of this study reveal that codon 12 and 61 mutant NRAS cells have distinct signaling characteristics that could allow for the development of more effective, mutation-specific treatment modalities.

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.

Searching for the Chokehold of NRAS Mutant Melanoma

Up to 18% of melanomas harbor mutations in the neuroblastoma rat-sarcoma homolog (NRAS). Yet, decades of research aimed to interfere with oncogenic RAS signaling have been largely disappointing and have not resulted in meaningful clinical outputs. Recent advances in disease modeling, structural biology, and an improved understanding of RAS cycling as well as RAS signaling networks have renewed hope for developing strategies to selectively block hyperactive RAS function. This review discusses direct and indirect blocking of activated RAS with a focus on current and potential future therapeutic approaches for NRAS mutant melanoma.

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.

KRAS insertion mutations are oncogenic and exhibit distinct functional properties

Oncogenic KRAS mutations introduce discrete amino acid substitutions that reduce intrinsic Ras GTPase activity and confer resistance to GTPase-activating proteins (GAPs). Here we discover a partial duplication of the switch 2 domain of K-Ras encoding a tandem repeat of amino acids G60_A66dup in a child with an atypical myeloproliferative neoplasm. K-Ras proteins containing this tandem duplication or a similar five amino acid E62_A66dup mutation identified in lung and colon cancers transform the growth of primary myeloid progenitors and of Ba/F3 cells. Recombinant K-Ras(G60_A66dup) and K-Ras(E62_A66dup) proteins display reduced intrinsic GTP hydrolysis rates, accumulate in the GTP-bound conformation and are resistant to GAP-mediated GTP hydrolysis. Remarkably, K-Ras proteins with switch 2 insertions are impaired for PI3 kinase binding and Akt activation, and are hypersensitive to MEK inhibition. These studies illuminate a new class of oncogenic KRAS mutations and reveal unexpected plasticity in oncogenic Ras proteins that has diagnostic and therapeutic implications.

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.

Direct Attack on RAS: Intramolecular Communication and Mutation-Specific Effects

The crystal structure of RAS was first solved 25 years ago. In spite of tremendous and sustained efforts, there are still no drugs in the clinic that directly target this major driver of human cancers. Recent success in the discovery of compounds that bind RAS and inhibit signaling has fueled renewed enthusiasm, and in-depth understanding of the structure and function of RAS has opened new avenues for direct targeting. To succeed, we must focus on the molecular details of the RAS structure and understand at a high-resolution level how the oncogenic mutants impair function. Structural networks of intramolecular communication between the RAS active site and membrane-interacting regions on the G-domain are disrupted in oncogenic mutants. Although conserved across the isoforms, these networks are near hot spots of protein-ligand interactions with amino acid composition that varies among RAS proteins. These differences could have an effect on stabilization of conformational states of interest in attenuating signaling through RAS. The development of strategies to target these novel sites will add a fresh direction in the quest to conquer RAS-driven cancers. Clin Cancer Res; 21(8); 1810-8. ©2015 AACR. See all articles in this CCR Focus section, "Targeting RAS-Driven Cancers."

Computational allosteric ligand binding site identification on Ras proteins

A number of computational techniques have been proposed to expedite the process of allosteric ligand binding site identification in inherently flexible and hence challenging drug targets. Some of these techniques have been instrumental in the discovery of allosteric ligand binding sites on Ras proteins, a group of elusive anticancer drug targets. This review provides an overview of these techniques and their application to Ras proteins. A summary of molecular docking and binding site identification is provided first, followed by a more detailed discussion of two specific techniques for binding site identification in ensembles of Ras conformations generated by molecular simulations.

Targeting the KRAS variant for treatment of non-small cell lung cancer: potential therapeutic applications

Lung cancer is the leading cause of cancer deaths worldwide, with non-small cell lung cancer (NSCLC) accounting for 80% of all lung cancers. Kirsten rat sarcoma viral oncogene homolog (KRAS) is one of the deadliest cancer-related proteins and plays a pivotal role in the most aggressive and lethal human cancers, including lung adenocarcinoma where it represents one of the most frequently mutated oncogene. Although therapeutic progresses have made an impact over the last decade, median survival for patients with advanced lung cancer remains disappointing, with a 5-year worldwide survival rate of <15%. For more than 20 years it has been recognized that constitutively active signaling downstream of KRAS is a fundamental driver of lung tumorigenesis. However, years of pursuit have failed to yield a drug that can safely curb KRAS activity; up to now no approved therapies exist for KRAS-mutant NSCLC. The aim of this review is to discuss the current knowledge of KRAS-mutated NSCLC, touching upon KRAS clinical relevance as a prognostic and predictive biomarker, with an emphasis on novel therapeutic approaches for the treatment of KRAS-variant NSCLC.