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

Cooperative genomic lesions in HRAS-mutant cancers predict resistance to farnesyltransferase inhibitors

In the clinical development of farnesyltransferase inhibitors (FTIs) for HRAS-mutant tumors, responses varied by cancer type. Co-occurring mutations may affect responses. We aimed to uncover cooperative genetic events specific to HRAS-mutant tumors and to study their effect on sensitivity to FTIs. Using targeted sequencing data from the MSK-IMPACT and Dana-Farber Cancer Institute Genomic Evidence Neoplasia Information Exchange databases, we identified comutations that were observed predominantly in HRAS-mutant versus KRAS-mutant or NRAS-mutant cancers. HRAS-mutant cancers had a higher frequency of coaltered mutations (48.8%) in the MAPK, PI3K, or RTK pathway genes, compared with KRAS-mutant (41.4%) and NRAS-mutant (38.4%) cancers (p < 0.05). Class 3 BRAF, NF1, PTEN, and PIK3CA mutations were more prevalent in HRAS-mutant lineages. To study the effects of comutations on sensitivity to FTIs, HrasG13R was transfected into "RASless" (Kraslox/lox/Hras-/-/Nras-/-/RERTert/ert) mouse embryonic fibroblasts (MEFs), which sensitized nontransfected MEFs to tipifarnib. Comutation in the form of Pten or Nf1 deletion and Pik3caH1047R transduction led to resistance to tipifarnib in HrasG13R-transfected MEFs in the presence or absence of KrasWT, whereas BrafG466E transduction led to resistance to tipifarnib only in the presence of KrasWT. Combined treatment with tipifarnib and MEK inhibition sensitized cells to tipifarnib in all settings, including in MEFs with PI3K pathway comutations. HRAS-mutant tumors demonstrate lineage-dependent MAPK or PI3K pathway alterations, which confer resistance to tipifarnib. The combined use of FTIs and MEK inhibition is a promising strategy for HRAS-mutant tumors.

Insight into structural dynamics involved in activation mechanism of full length KRAS wild type and P-loop mutants

KRAS protein is known to be frequently mutated in various cancers. The most common mutations being at position 12, 13 and 61. The positions 12 and 13 form part of the phosphate binding region (P-loop) of KRAS. Owing to mutation, the protein remains in continuous active state and affects the normal cellular process. Understanding the structural changes owing to mutations in GDP-bound (inactive state) and GTP-bound (active state) may help in the design of better therapeutics. To understand the structural flexibility due to the mutations specifically located at P-loop regions (G12D, G12V and G13D), extensive molecular dynamics simulations (24 μs) have been carried for both inactive (GDP-bound) and active (GTP-bound) structures for the wild type and these mutants. The study revealed that the local structural changes at the site of mutations allosterically guide changes in distant regions of the protein through hydrogen bond and hydrophobic signalling network. The dynamic cross correlation analysis and the comparison of the correlated motions among different systems manifested that changes in SW-I, SW-II, α3 and the loop preceding α3 affects the interactions of GDP/GTP with different regions of the protein thereby affecting its hydrolysis. Further, the Markov state modelling analysis confirmed that the mutations, especially G13D imparts rigidity to structure compared to wild type and thus limiting its conformational state in either intermediate state or active state. The study suggests that along with SW-I and SW-II regions, the loop region preceding the α3 helix and α3 helix are also involved in affecting the hydrolysis of nucleotides and may be considered while designing therapeutics against KRAS.

Insights into the Anticancer Mechanisms Modulated by Gamma and Delta Tocotrienols in Colorectal Cancers

Colorectal cancer (CRC) is a growing concern all over the world. There has been a concerted effort to identify natural bioactive compounds that can be used to prevent or overcome this condition. Tocotrienols (T3s) are a naturally occurring form of vitamin E known for various therapeutic effects, such as anticancer, antioxidant, neuroprotective, and anti-inflammatory activities. The literature evidence suggests that two T3 analogues, ie, gamma (γ)- and delta (δ)-T3, can modulate cancers via several cancer-related signaling pathways. The aim of this review was to compile and analyze the existing literature on the diverse anticancer mechanisms of γT3 and δT3 exhibited in CRC cells, to showcase the anticancer potential of T3s. Medline was searched for research articles on anticancer effects of γT3 and δT3 in CRC published in the past 2 decades. A total of 38 articles (26 cell-based, 9 animal studies, 2 randomized clinical trials, and 1 scoping review) that report anticancer effects of γT3 and δT3 in CRC were identified. The findings reported in those articles indicate that γT3 and δT3 inhibit the proliferation of CRC cells, induce cell cycle arrest and apoptosis, suppress metastasis, and produce synergistic anticancer effects when combined with well-established anticancer agents. There is preliminary evidence that shows that T3s affect telomerase functions and support anticancer immune responses. γT3 and δT3 have the potential for development as anticancer agents.

High Pressure Promotes Binding of the Allosteric Inhibitor Zn2+-Cyclen in Crystals of Activated H-Ras

In this work, we experimentally investigate the potency of high pressure to drive a protein toward an excited state where an inhibitor targeted for this state can bind. Ras proteins are small GTPases cycling between active GTP-bound and inactive GDP-bound states. Various states of GTP-bound Ras in active conformation coexist in solution, amongst them, state 2 which binds to effectors, and state 1, weakly populated at ambient conditions, which has a low affinity for effectors. Zn2+-cyclen is an allosteric inhibitor of Ras protein, designed to bind specifically to the state 1. In H-Ras(wt).Mg2+.GppNHp crystals soaked with Zn2+-cyclen, no binding could be observed, as expected in the state 2 conformation which is the dominant state at ambient pressure. Interestingly, Zn2+-cyclen binding is observed at 500 MPa pressure, close to the nucleotide, in Ras protein that is driven by pressure to a state 1 conformer. The unknown binding mode of Zn2+-cyclen to H-Ras can thus be fully characterized in atomic details. As a more general conjunction from our study, high pressure x-ray crystallography turns out to be a powerful method to induce transitions allowing drug binding in proteins that are in low-populated conformations at ambient conditions, enabling the design of specific inhibitors.

Identification and analysis of oncogenic non-synonymous single nucleotide polymorphisms in the human NRAS gene: An exclusive in silico study

BACKGROUND

N-ras protein is encoded by the NRAS gene and operates as GDP-GTP-controlled on/off switching. N-ras interacts with cellular signaling networks that regulate various cellular activities including cell proliferation and survival. The nonsynonymous single nucleotide polymorphism (nsSNPs)-mediated alteration can substantially disrupt the structure and activity of the corresponding protein. N-ras has been reported to be associated with numerous diseases including cancers due to the nsSNPs. A comprehensive study on the NRAS gene to unveil the potentially damaging and oncogenic nsSNPs is yet to be accomplished. Hence, this extensive in silico study is intended to identify the disease-associated, specifically oncogenic nsSNPs of the NRAS gene.

RESULTS

Out of 140 missense variants, 7 nsSNPs (I55R, G60E, G60R, Y64D, L79F, D119G, and V152F) were identified to be damaging utilizing 10 computational tools that works based on different algorithms with high accuracy. Among those, G60E, G60R, and D119G variants were further filtered considering their location in the highly conserved region and later identified as oncogenic variants. Interestingly, G60E and G60R variants were revealed to be particularly associated with lung adenocarcinoma, rhabdomyosarcoma, and prostate adenocarcinoma. Therefore, D119G could be subjected to detailed investigation for identifying its association with specific cancer.

CONCLUSION

This in silico study identified the deleterious and oncogenic missense variants of the human NRAS gene that could be utilized for designing further experimental investigation. The outcomes of this study would be worthwhile in future research for developing personalized medicine.

Crystal structure of NRAS Q61K with a ligand-induced pocket near switch II

The RAS isoforms (KRAS, HRAS and NRAS) have distinct cancer type-specific profiles. NRAS mutations are the second most prevalent RAS mutations in skin and hematological malignancies. Although RAS proteins were considered undruggable for decades, isoform and mutation-specific investigations have produced successful RAS inhibitors that are either specific to certain mutants, isoforms (pan-KRAS) or target all RAS proteins (pan-RAS). While extensive structural and biochemical investigations have focused mainly on K- and H-RAS mutations, NRAS mutations have received less attention, and the most prevalent NRAS mutations in human cancers, Q61K and Q61R, are rare in K- and H-RAS. This manuscript presents a crystal structure of the NRAS Q61K mutant in the GTP-bound form. Our structure reveals a previously unseen pocket near switch II induced by the binding of a ligand to the active form of the protein. This observation reveals a binding site that can potentially be exploited for development of inhibitors against mutant NRAS. Furthermore, the well-resolved catalytic site of this GTPase bound to native GTP provides insight into the stalled GTP hydrolysis observed for NRAS-Q61K.

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.

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

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

Revealing the mechanism of action of a first-in-class covalent inhibitor of KRASG12C (ON) and other functional properties of oncogenic KRAS by 31P NMR

Individual oncogenic KRAS mutants confer distinct differences in biochemical properties and signaling for reasons that are not well understood. KRAS activity is closely coupled to protein dynamics and is regulated through two interconverting conformations: state 1 (inactive, effector binding deficient) and state 2 (active, effector binding enabled). Here, we use 31P NMR to delineate the differences in state 1 and state 2 populations present in WT and common KRAS oncogenic mutants (G12C, G12D, G12V, G13D, and Q61L) bound to its natural substrate GTP or a commonly used nonhydrolyzable analog GppNHp (guanosine-5'-[(β,γ)-imido] triphosphate). Our results show that GppNHp-bound proteins exhibit significant state 1 population, whereas GTP-bound KRAS is primarily (90% or more) in state 2 conformation. This observation suggests that the predominance of state 1 shown here and in other studies is related to GppNHp and is most likely nonexistent in cells. We characterize the impact of this differential conformational equilibrium of oncogenic KRAS on RAF1 kinase effector RAS-binding domain and intrinsic hydrolysis. Through a KRAS G12C drug discovery, we have identified a novel small-molecule inhibitor, BBO-8956, which is effective against both GDP- and GTP-bound KRAS G12C. We show that binding of this inhibitor significantly perturbs state 1-state 2 equilibrium and induces an inactive state 1 conformation in GTP-bound KRAS G12C. In the presence of BBO-8956, RAF1-RAS-binding domain is unable to induce a signaling competent state 2 conformation within the ternary complex, demonstrating the mechanism of action for this novel and active-conformation inhibitor.

GTP-Bound N-Ras Conformational States and Substates Are Modulated by Membrane and Point Mutation

Oncogenic Ras proteins are known to present multiple conformational states, as reported by the great variety of crystallographic structures. The GTP-bound states are grouped into two main states: the "inactive" state 1 and the "active" state 2. Recent reports on H-Ras have shown that state 2 exhibits two substates, directly related to the orientation of Tyr32: toward the GTP-bound pocket and outwards. In this paper, we show that N-Ras exhibits another substate of state 2, related to a third orientation of Tyr32, toward Ala18 and parallel to the GTP-bound pocket. We also show that this substate is highly sampled in the G12V mutation of N-Ras and barely present in its wild-type form, and that the G12V mutation prohibits the sampling of the GTPase-activating protein (GAP) binding substate, rendering this mutation oncogenic. Furthermore, using molecular dynamics simulations, we explore the importance of the membrane on N-Ras' conformational state dynamics and its strong influence on Ras protein stability. Moreover, the membrane has a significant influence on the conformational (sub)states sampling of Ras. This, in turn, is of crucial importance in the activation/deactivation cycle of Ras, due to the binding of guanine nucleotide exchange factor proteins (GEFs)/GTPase-activating proteins (GAPs).

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

NLRP3 inflammasome activation and symptom burden in KRAS-mutated CMML patients is reverted by IL-1 blocking therapy

Chronic myelomonocytic leukemia (CMML) is frequently associated with mutations in the rat sarcoma gene (RAS), leading to worse prognosis. RAS mutations result in active RAS-GTP proteins, favoring myeloid cell proliferation and survival and inducing the NLRP3 inflammasome together with the apoptosis-associated speck-like protein containing a caspase recruitment domain (ASC), which promote caspase-1 activation and interleukin (IL)-1β release. Here, we report, in a cohort of CMML patients with mutations in KRAS, a constitutive activation of the NLRP3 inflammasome in monocytes, evidenced by ASC oligomerization and IL-1β release, as well as a specific inflammatory cytokine signature. Treatment of a CMML patient with a KRASG12D mutation using the IL-1 receptor blocker anakinra inhibits NLRP3 inflammasome activation, reduces monocyte count, and improves the patient's clinical status, enabling a stem cell transplant. This reveals a basal inflammasome activation in RAS-mutated CMML patients and suggests potential therapeutic applications of NLRP3 and IL-1 blockers.

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.

Neuroblastoma RAS viral oncogene homolog (N-RAS) deficiency aggravates liver injury and fibrosis

Progressive hepatic damage and fibrosis are major features of chronic liver diseases of different etiology, yet the underlying molecular mechanisms remain to be fully defined. N-RAS, a member of the RAS family of small guanine nucleotide-binding proteins also encompassing the highly homologous H-RAS and K-RAS isoforms, was previously reported to modulate cell death and renal fibrosis; however, its role in liver damage and fibrogenesis remains unknown. Here, we approached this question by using N-RAS deficient (N-RAS-/-) mice and two experimental models of liver injury and fibrosis, namely carbon tetrachloride (CCl4) intoxication and bile duct ligation (BDL). In wild-type (N-RAS+/+) mice both hepatotoxic procedures augmented N-RAS expression in the liver. Compared to N-RAS+/+ counterparts, N-RAS-/- mice subjected to either CCl4 or BDL showed exacerbated liver injury and fibrosis, which was associated with enhanced hepatic stellate cell (HSC) activation and leukocyte infiltration in the damaged liver. At the molecular level, after CCl4 or BDL, N-RAS-/- livers exhibited augmented expression of necroptotic death markers along with JNK1/2 hyperactivation. In line with this, N-RAS ablation in a human hepatocytic cell line resulted in enhanced activation of JNK and necroptosis mediators in response to cell death stimuli. Of note, loss of hepatic N-RAS expression was characteristic of chronic liver disease patients with fibrosis. Collectively, our study unveils a novel role for N-RAS as a negative controller of the progression of liver injury and fibrogenesis, by critically downregulating signaling pathways leading to hepatocyte necroptosis. Furthermore, it suggests that N-RAS may be of potential clinical value as prognostic biomarker of progressive fibrotic liver damage, or as a novel therapeutic target for the treatment of chronic liver disease.

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 KRASG12C inhibitors on different RAS allelic mutations or even different cancers with KRASG12C 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.

Direct GDP-KRASG12C inhibitors and mechanisms of resistance: the tip of the iceberg

Kirsten rat sarcoma viral oncogene homolog mutations are observed in 25% of lung adenocarcinoma and 40% of these are G12C mutations. Historically, no approved targeted agents were available for patients with any KRAS mutation, and response rates to standard-of-care therapies were suboptimal. Newly developed inhibitors directed toward KRAS^G12C have been successful in clinical trials with overall response rates ranging between 32% and 46%, and two FDA approvals were granted in May 2021 and December 2022 as second-line or later monotherapies. However, rapid tumor resistance complicates their use as a monotherapy. With the rapid development of this novel class of inhibitors, it is important to discern the different types of tumor resistance that may arise and how each can differently contribute to tumor growth and survival. G12C inhibitor resistance is under investigation and combinations of therapies with G12C inhibitors have been proposed. Much of this insight is gleaned from preclinical investigations, as our knowledge of clinical resistance is in its infancy. In this review, we summarize the preclinical development of KRAS^G12C inhibitors, their clinical evaluations, different types of resistance mechanisms to these compounds, and ways of overcoming them. Finally, we underscore the importance of basic and translational investigations of these molecules in a landscape where their clinical evaluations garner the most attention, and we set the stage for what is to come.

Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations

Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.

Redox-sensitive small GTPase H-Ras in murine astrocytes, an in vitro study

BACKGROUND

Although the protooncogenes small GTPases Ras are redox-sensitive proteins, how they are regulated by redox signaling in the central nervous system (CNS) is still poorly understood. Alteration in redox-sensitive targets by redox signaling may have myriad effects on Ras stability, activity and localization. Redox-mediated changes in astrocytic RAS may contribute to the control of redox homeostasis in the CNS that is connected to the pathogenesis of many diseases.

RESULTS AND METHODS

Here, we investigated the transient physiological induction, at both transcriptional and translational levels, of small GTPases Ras in response to redox stimulation. Cultured astrocytes were treated with hydrogen peroxide as in bolus addition and relative mRNA levels of murine hras and kras genes were detected by qRT-PCR. We found that de novo transcription of hras mRNA in reactive astrocytes is redox-sensitive and mimics the prototypical redox-sensitive gene iNOS. Protein abundance in combination with protein turnover measurements by cycloheximide-chase experiments revealed distinct translation efficiency, GTP-bound enrichment, and protein turnover rates between the two isoforms H-Ras and K-Ras.

CONCLUSION

Reports from recent years support a significant role of H-Ras in driving redox processes. Beyond its canonical functions, Ras may impact on the core astrocytic cellular machinery that operates during redox stimulation.

The Importance of Mg2+ -Free State in Nucleotide Exchange of Oncogenic K-Ras Mutants

For efficient targeting of oncogenic K-Ras interaction sites, a mechanistic picture of the Ras-cycle is necessary. Herein, we used NMR relaxation techniques and molecular dynamics simulations to decipher the role of slow dynamics in wild-type and three oncogenic P-loop mutants of K-Ras. Our measurements reveal a dominant two-state conformational exchange on the ms timescale in both GDP- and GTP-bound K-Ras. The identified low-populated higher energy state in GDP-loaded K-Ras has a conformation reminiscent of a nucleotide-bound/Mg²⁺ -free state characterized by shortened β2/β3-strands and a partially released switch-I region preparing K-Ras for the interaction with the incoming nucleotide exchange factor and subsequent reactivation. By providing insight into mutation-specific differences in K-Ras structural dynamics, our systematic analysis improves our understanding of prolonged K-Ras signaling and may aid the development of allosteric inhibitors targeting nucleotide exchange in K-Ras.

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

Ras dimers have been proposed as building blocks for initiating the extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) cellular signaling pathway. To better examine the structure of possible dimer interfaces, the dynamics of Ras dimerization, and its potential signaling consequences, we performed molecular dynamics simulations totaling 1 ms of sampling, using an all-atom model of two full-length, farnesylated, guanosine triphosphate (GTP)-bound, wild-type KRas4b proteins diffusing on 29%POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine)-mixed POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Our simulations unveil an ensemble of thermodynamically weak KRas dimers spanning multiple conformations. The most stable conformations, having the largest interface areas, involve helix α2 and a hypervariable region (HVR). Among the dimer conformations, we found that the HVR of each KRas has frequent interactions with various parts of the dimer, thus potentially mediating the dimerization. Some dimer configurations have one KRas G-domain elevated above the lipid bilayer surface by residing on top of the other G-domain, thus likely contributing to the recruitment of cytosolic Raf kinases in the context of a stably formed multi-protein complex. We identified a variant of the α4-α5 KRas-dimer interface that is similar to the interfaces obtained with fluorescence resonance energy transfer (FRET) data of HRas on lipid bilayers. Interestingly, we found two arginine fingers, R68 and R149, that directly interact with the beta-phosphate of the GTP bound in KRas, in a manner similar to what is observed in a crystal structure of GAP-HRas complex, which can facilitate the GTP hydrolysis via the arginine finger of GTPase-activating protein (GAP).

Mechanisms of isoform-specific residue influence on GTP-bound HRas, KRas, and NRas

HRas, KRas, and NRas are GTPases with a common set of effectors that control many cell-signaling pathways, including proliferation through Raf kinase. Their G-domains are nearly identical in sequence, with a few isoform-specific residues that have an effect on dynamics and biochemical properties. Here, we use accelerated molecular dynamics (aMD) simulations consistent with solution x-ray scattering experiments to elucidate mechanisms through which isoform-specific residues associated with each Ras isoform affects functionally important regions connected to the active site. HRas-specific residues cluster in loop 8 to stabilize the nucleotide-binding pocket, while NRas-specific residues on helix 3 directly affect the conformations of switch I and switch II. KRas, the most globally flexible of the isoforms, shows greatest fluctuations in the switch regions enhanced by a KRas-specific residue in loop 7 and a highly dynamic loop 8 region. The analysis of isoform-specific residue effects on Ras proteins is supported by NMR experiments and is consistent with previously published biochemical data.

Structure of the MRAS-SHOC2-PP1C phosphatase complex

RAS-MAPK signalling is fundamental for cell proliferation and is altered in most human cancers1-3. However, our mechanistic understanding of how RAS signals through RAF is still incomplete. Although studies revealed snapshots for autoinhibited and active RAF-MEK1-14-3-3 complexes4, the intermediate steps that lead to RAF activation remain unclear. The MRAS-SHOC2-PP1C holophosphatase dephosphorylates RAF at serine 259, resulting in the partial displacement of 14-3-3 and RAF-RAS association3,5,6. MRAS, SHOC2 and PP1C are mutated in rasopathies-developmental syndromes caused by aberrant MAPK pathway activation6-14-and SHOC2 itself has emerged as potential target in receptor tyrosine kinase (RTK)-RAS-driven tumours15-18. Despite its importance, structural understanding of the SHOC2 holophosphatase is lacking. Here we determine, using X-ray crystallography, the structure of the MRAS-SHOC2-PP1C complex. SHOC2 bridges PP1C and MRAS through its concave surface and enables reciprocal interactions between all three subunits. Biophysical characterization indicates a cooperative assembly driven by the MRAS GTP-bound active state, an observation that is extendible to other RAS isoforms. Our findings support the concept of a RAS-driven and multi-molecular model for RAF activation in which individual RAS-GTP molecules recruit RAF-14-3-3 and SHOC2-PP1C to produce downstream pathway activation. Importantly, we find that rasopathy and cancer mutations reside at protein-protein interfaces within the holophosphatase, resulting in enhanced affinities and function. Collectively, our findings shed light on a fundamental mechanism of RAS biology and on mechanisms of clinically observed enhanced RAS-MAPK signalling, therefore providing the structural basis for therapeutic interventions.

The current state of the art and future trends in RAS-targeted cancer therapies

Despite being the most frequently altered oncogenic protein in solid tumours, KRAS has historically been considered ‘undruggable’ owing to a lack of pharmacologically targetable pockets within the mutant isoforms. However, improvements in drug design have culminated in the development of inhibitors that are selective for mutant KRAS in its active or inactive state. Some of these inhibitors have proven efficacy in patients with KRAS^G12C-mutant cancers and have become practice changing. The excitement associated with these advances has been tempered by drug resistance, which limits the depth and/or duration of responses to these agents. Improvements in our understanding of RAS signalling in cancer cells and in the tumour microenvironment suggest the potential for several novel combination therapies, which are now being explored in clinical trials. Herein, we provide an overview of the RAS pathway and review the development and current status of therapeutic strategies for targeting oncogenic RAS, as well as their potential to improve outcomes in patients with RAS-mutant malignancies. We then discuss challenges presented by resistance mechanisms and strategies by which they could potentially be overcome.

The RAS oncogenes are among the most common drivers of tumour development and progression but have historically been considered undruggable. The development of direct KRAS inhibitors has changed this paradigm, although currently clinical use of these novel therapeutics is limited to a select subset of patients, and intrinsic or acquired resistance presents an inevitable challenge to cure. Herein, the authors provide an overview of the RAS pathway in cancer and review the ongoing efforts to develop effective therapeutic strategies for RAS-mutant cancers. They also discuss the current understanding of mechanisms of resistance to direct KRAS inhibitors and strategies by which they might be overcome.

Owing to intrinsic and extrinsic factors, KRAS and other RAS isoforms have until recently been impervious to targeting with small-molecule inhibitors.

Inhibitors of the KRAS^G12C variant constitute a potential breakthrough in the treatment of many cancer types, particularly non-small-cell lung cancer, for which such an agent has been approved by the FDA.

Several forms of resistance to KRAS inhibitors have been defined, including primary, adaptive and acquired resistance; these resistance mechanisms are being targeted in studies that combine KRAS inhibitors with inhibitors of horizontal or vertical signalling pathways.

Mutant KRAS has important effects on the tumour microenvironment, including the immunological milieu; these effects must be considered to fully understand resistance to KRAS inhibitors and when designing novel treatment strategies.

Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation

K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.

Computational Analysis of Deleterious SNPs in NRAS to Assess Their Potential Correlation With Carcinogenesis

The NRAS gene is a well-known oncogene that acts as a major player in carcinogenesis. Mutations in the NRAS gene have been linked to multiple types of human tumors. Therefore, the identification of the most deleterious single nucleotide polymorphisms (SNPs) in the NRAS gene is necessary to understand the key factors of tumor pathogenesis and therapy. We aimed to retrieve NRAS missense SNPs and analyze them comprehensively using sequence and structure approaches to determine the most deleterious SNPs that could increase the risk of carcinogenesis. We also adopted structural biology methods and docking tools to investigate the behavior of the filtered SNPs. After retrieving missense SNPs and analyzing them using six in silico tools, 17 mutations were found to be the most deleterious mutations in NRAS. All SNPs except S145L were found to decrease NRAS stability, and all SNPs were found on highly conserved residues and important functional domains, except R164C. In addition, all mutations except G60E and S145L showed a higher binding affinity to GTP, implicating an increase in malignancy tendency. As a consequence, all other 14 mutations were expected to increase the risk of carcinogenesis, with 5 mutations (G13R, G13C, G13V, P34R, and V152F) expected to have the highest risk. Thermodynamic stability was ensured for these SNP models through molecular dynamics simulation based on trajectory analysis. Free binding affinity toward the natural substrate, GTP, was higher for these models as compared to the native NRAS protein. The Gly13 SNP proteins depict a differential conformational state that could favor nucleotide exchange and catalytic potentiality. A further application of experimental methods with all these 14 mutations could reveal new insights into the pathogenesis and management of different types of tumors.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Prediction of Allosteric Communication Pathways in Proteins

MOTIVATION

Allostery in proteins is an essential phenomenon in biological processes. In this paper, we present a computational model to predict paths of maximum information transfer between active and allosteric sites. In this information theoretic study, we use mutual information as the measure of information transfer, where transition probability of information from one residue to its contacting neighbors is proportional to the magnitude of mutual information between the two residues. Starting from a given residue and using a Hidden Markov Model, we successively determine the neighboring residues that eventually lead to a path of optimum information transfer. The Gaussian approximation of mutual information between residue pairs is adopted. The limits of validity of this approximation are discussed in terms of a nonlinear theory of mutual information and its reduction to the Gaussian form.

RESULTS

Predictions of the model are tested on six widely studied cases, CheY Bacterial Chemotaxis, B-cell Lymphoma extra-large Bcl-xL, Human proline isomerase cyclophilin A (CypA), Dihydrofolate reductase DHFR, HRas GTPase, and Caspase-1. The communication transmission rendering the propagation of local fluctuations from the active sites throughout the structure in multiple paths correlate well with the known experimental data. Distinct paths originating from the active site may likely represent a multi functionality such as involving more than one allosteric site and/or preexistence of some other functional states. Our model is computationally fast and simple, and can give allosteric communication pathways, which are crucial for the understanding and control of protein functionality.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Enhanced BRAF engagement by NRAS mutants capable of promoting melanoma initiation

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

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

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

Insights into the Cross Talk between Effector and Allosteric Lobes of KRAS from Methyl Conformational Dynamics

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-μs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

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.

Extracting phylogenetic dimensions of coevolution reveals hidden functional signals

Despite the structural and functional information contained in the statistical coupling between pairs of residues in a protein, coevolution associated with function is often obscured by artifactual signals such as genetic drift, which shapes a protein's phylogenetic history and gives rise to concurrent variation between protein sequences that is not driven by selection for function. Here, we introduce a background model for phylogenetic contributions of statistical coupling that separates the coevolution signal due to inter-clade and intra-clade sequence comparisons and demonstrate that coevolution can be measured on multiple phylogenetic timescales within a single protein. Our method, nested coevolution (NC), can be applied as an extension to any coevolution metric. We use NC to demonstrate that poorly conserved residues can nonetheless have important roles in protein function. Moreover, NC improved the structural-contact predictions of several coevolution-based methods, particularly in subsampled alignments with fewer sequences. NC also lowered the noise in detecting functional sectors of collectively coevolving residues. Sectors of coevolving residues identified after application of NC were more spatially compact and phylogenetically distinct from the rest of the protein, and strongly enriched for mutations that disrupt protein activity. Thus, our conceptualization of the phylogenetic separation of coevolution provides the potential to further elucidate relationships among protein evolution, function, and genetic diseases.

Enhanced interpretation of 935 hotspot and non-hotspot RAS variants using evidence-based structural bioinformatics

In the current study, we report computational scores for advancing genomic interpretation of disease-associated genomic variation in members of the RAS family of genes. For this purpose, we applied 31 sequence- and 3D structure-based computational scores, chosen by their breadth of biophysical properties. We parametrized our data by assembling a numerically homogenized experimentally-derived dataset, which when use in our calculations reveal that computational scores using 3D structure highly correlate with experimental measures (e.g., GAP-mediated hydrolysis R_Spearman = 0.80 and RAF affinity R_spearman = 0.82), while sequence-based scores are discordant with this data. Performing all-against-all comparisons, we applied this parametrized modeling approach to the study of 935 RAS variants from 7 RAS genes, which led us to identify 4 groups of mutations according to distinct biochemical scores within each group. Each group was comprised of hotspot and non-hotspot KRAS variants, indicating that poorly characterized variants could functionally behave like pathogenic mutations. Combining computational scores using dimensionality reduction indicated that changes to local unfolding propensity associate with changes in enzyme activity by genomic variants. Hence, our systematic approach, combining methodologies from both clinical genomics and 3D structural bioinformatics, represents an expansion for interpreting genomic data, provides information of mechanistic value, and that is transferable to other proteins.

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.

Inhibition of mutant RAS-RAF interaction by mimicking structural and dynamic properties of phosphorylated RAS

Undruggability of RAS proteins has necessitated alternative strategies for the development of effective inhibitors. In this respect, phosphorylation has recently come into prominence as this reversible post-translational modification attenuates sensitivity of RAS towards RAF. As such, in this study, we set out to unveil the impact of phosphorylation on dynamics of HRAS^WT and aim to invoke similar behavior in HRAS^G12D mutant by means of small therapeutic molecules. To this end, we performed molecular dynamics (MD) simulations using phosphorylated HRAS and showed that phosphorylation of Y32 distorted Switch I, hence the RAS/RAF interface. Consequently, we targeted Switch I in HRAS^G12D by means of approved therapeutic molecules and showed that the ligands enabled detachment of Switch I from the nucleotide-binding pocket. Moreover, we demonstrated that displacement of Switch I from the nucleotide-binding pocket was energetically more favorable in the presence of the ligand. Importantly, we verified computational findings in vitro where HRAS^G12D/RAF interaction was prevented by the ligand in HEK293T cells that expressed HRAS^G12D mutant protein. Therefore, these findings suggest that targeting Switch I, hence making Y32 accessible might open up new avenues in future drug discovery strategies that target mutant RAS proteins.

Targeted Long-Read Sequencing Decodes the Transcriptional Atlas of the Founding RAS Gene Family Members

The complicity of human RAS proteins in cancer is a well-documented fact, both due to the mutational hyperactivation of these GTPases and the overexpression of the genes encoding these proteins. Thus, it can be easily assumed that the study of RAS genes at the transcriptional and post-transcriptional level is of the utmost importance. Although previous research has shed some light on the basic mechanisms by which GTPases are involved in tumorigenesis, limited information is known regarding the transcriptional profile of the genes encoding these proteins. The present study highlights for the first time the wide spectrum of the mRNAs generated by the three most significant RAS genes (KRAS, NRAS and HRAS), providing an in-depth analysis of the splicing events and exon/intron boundaries. The implementation of a versatile, targeted nanopore-sequencing approach led to the identification of 39 novel RAS mRNA transcript variants and to the elucidation of their expression profiles in a broad panel of human cell lines. Although the present work unveiled multiple hidden aspects of the RAS gene family, further study is required to unravel the biological function of all the novel alternative transcript variants, as well as the putative protein isoforms.

SHP2 inhibition diminishes KRASG12C cycling and promotes tumor microenvironment remodeling

KRAS is the most frequently mutated human oncogene, and KRAS inhibition has been a longtime goal. Recently, inhibitors were developed that bind KRASG12C-GDP and react with Cys-12 (G12C-Is). Using new affinity reagents to monitor KRASG12C activation and inhibitor engagement, we found that an SHP2 inhibitor (SHP2-I) increases KRAS-GDP occupancy, enhancing G12C-I efficacy. The SHP2-I abrogated RTK feedback signaling and adaptive resistance to G12C-Is in vitro, in xenografts, and in syngeneic KRASG12C-mutant pancreatic ductal adenocarcinoma (PDAC) and non-small cell lung cancer (NSCLC). SHP2-I/G12C-I combination evoked favorable but tumor site-specific changes in the immune microenvironment, decreasing myeloid suppressor cells, increasing CD8+ T cells, and sensitizing tumors to PD-1 blockade. Experiments using cells expressing inhibitor-resistant SHP2 showed that SHP2 inhibition in PDAC cells is required for PDAC regression and remodeling of the immune microenvironment but revealed direct inhibitory effects on tumor angiogenesis and vascularity. Our results demonstrate that SHP2-I/G12C-I combinations confer a substantial survival benefit in PDAC and NSCLC and identify additional potential combination strategies.

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.

Inhibitory activities of 20(R, S)-protopanaxatriol against epidermal growth factor receptor tyrosine kinase

As major metabolites of protopanaxatriol-type ginsenosides, 20(R, S)-protopanaxatriol [20(R, S)-PPT] display multiple bioactivities. This work aimed to investigate the inhibitory activities of 20(R, S)-PPT against epidermal growth factor receptor tyrosine kinase and the potential mechanism. 20(R, S)-PPT inhibited the proliferation of HepG2 cells in a dose-dependent manner and blocked cell cycle progression at G1/G0 phase. Then 20(R, S)-PPT were found to influence the protein expressions involved in epidermal growth factor receptor (EGFR)-mitogen-activated protein kinase (MAPK) signaling pathway. Molecular docking suggested that 20(R, S)-PPT could bind to the active sites of all target proteins in EGFR-MAPK pathway. It is worth noting that 20(R, S)-PPT showed stronger binding capacities with EGFR, compared with other proteins. Hence, this work further investigated the binding interactions and binding stabilities between 20(R, S)-PPT and EGFR. Both hydrophobic interactions and hydrogen bonds contributed to the 20(R, S)-PPT-EGFR binding. In addition, the in vitro inhibitory activities of 20(R, S)-PPT against EGFR tyrosine kinase were observed in a homogeneous time-resolved fluorescence assay, with the IC₅₀ values of 24.10 ± 0.17 and 33.19 ± 0.19 μM respectively. Taken together with the above results, both of 20(R)-PPT and 20(S)-PPT might serve as potential EGFR tyrosine kinase inhibitors.

Crystal Structure Reveals the Full Ras-Raf Interface and Advances Mechanistic Understanding of Raf Activation

Ras and Raf-kinase interact through the Ras-binding (RBD) and cysteine-rich domains (CRD) of Raf to signal through the mitogen-activated protein kinase pathway, yet the molecular mechanism leading to Raf activation has remained elusive. We present the 2.8 Å crystal structure of the HRas-CRaf-RBD_CRD complex showing the Ras-Raf interface as a continuous surface on Ras, as seen in the KRas-CRaf-RBD_CRD structure. In molecular dynamics simulations of a Ras dimer model formed through the α4-α5 interface, the CRD is dynamic and located between the two Ras protomers, poised for direct or allosteric modulation of functionally relevant regions of Ras and Raf. We propose a molecular model in which Ras binding is involved in the release of Raf autoinhibition while the Ras-Raf complex dimerizes to promote a platform for signal amplification, with Raf-CRD centrally located to impact regulation and function.

Antipsychotic phenothiazine drugs bind to KRAS in vitro

We used NMR to show that the antipsychotic phenothiazine drugs promazine and promethazine bind to GDP-KRAS. Promazine also binds to oncogenic GDP-KRAS(G12D), and to wild type GppNHp-KRAS. A panel of additional phenothiazines bind to GDP-KRAS but with lower affinity than promazine or promethazine. Binding is most dependent on substitutions at C-2 of the tricyclic phenothiazine ring. Promazine was used to generate an NMR-driven HADDOCK model of the drug/GDP-KRAS complex. The structural model shows the tricyclic phenothiazine ring of promazine associates with the hydrophobic pocket p1 that is bordered by the central β sheet and Switch II in KRAS. Binding appears to stabilize helix 2 in a conformation that is similar to that seen in KRAS bound to other small molecules. Association of phenothiazines with KRAS may affect normal KRAS signaling that could contribute to multiple biological activities of these antipsychotic drugs. Moreover, the phenothiazine ring represents a new core scaffold on which to design modulators of KRAS activity.

Structural bioinformatics enhances mechanistic interpretation of genomic variation, demonstrated through the analyses of 935 distinct RAS family mutations

MOTIVATION

Protein-coding genetic alterations are frequently observed in Clinical Genetics, but the high yield of variants of uncertain significance remains a limitation in decision making. RAS-family GTPases are cancer drivers, but only 54 variants, across all family members, fall within well-known hotspots. However, extensive sequencing has identified 881 non-hotspot variants for which significance remains to be investigated.

RESULTS

Here, we evaluate 935 missense variants from seven RAS genes, observed in cancer, RASopathies and the healthy adult population. We characterized hotspot variants, previously studied experimentally, using 63 sequence- and 3D structure-based scores, chosen by their breadth of biophysical properties. Applying scores that display best correlation with experimental measures, we report new valuable mechanistic inferences for both hot-spot and non-hotspot variants. Moreover, we demonstrate that 3D scores have little-to-no correlation with those based on DNA sequence, which are commonly used in Clinical Genetics. Thus, combined, these new knowledge bear significant relevance.

AVAILABILITY AND IMPLEMENTATION

All genomic and 3D scores, and markdown for generating figures, are provided in our supplemental data.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Ras subfamily GTPases regulate development, aflatoxin biosynthesis and pathogenicity in the fungus Aspergillus flavus

Ras subfamily proteins are molecular switches in signal transduction pathways of many eukaryotes that regulate a variety of cellular processes. Here, the Ras subfamily, encoded by six genes, was identified in Aspergillus flavus: rasA, rasB, rasC, rab-33, rheb and rsr1. The rsr1 deletion mutant (∆rsr1), rheb deletion mutant (∆rheb) and double deletion mutant (∆rheb/rsr1) displayed significantly decreased growth and sporulation. Sclerotia formation was significantly decreased for ∆rheb or ∆rheb/rsr1 but increased for ∆rsr1. Aflatoxin production was significantly increased in ∆rheb but decreased in ∆rsr1 and ∆rheb/rsr1. We found that rsr1 and rheb are crucial for the pathogenicity of A. flavus. Quantitative proteomics identified 520 differentially expressed proteins (DEPs) for the ∆rsr1 mutant and 133 DEPs for the ∆rheb mutant. These DEPs were annotated in multiple biological processes and KEGG pathways in A. flavus. Importantly, we identified the cytokinesis protein SepA in the protein-protein interaction network of rsr1, and deletion mutants showed that SepA has pleiotropic effects on growth and AF biosynthesis, which may depend on Rsr1 for regulation in A. flavus. Our results indicated that these Ras subfamily proteins exhibited functional redundancy with each other but there were also differences in A. flavus.

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.

Monoubiquitination of KRAS at Lysine104 and Lysine147 Modulates Its Dynamics and Interaction with Partner Proteins

KRAS, a 21 kDa guanine nucleotide-binding protein that functions as a molecular switch, plays a key role in regulating cellular growth. Dysregulation of this key signaling node leads to uncontrolled cell growth, a hallmark of cancer cells. KRAS undergoes post-translational modification by monoubiquitination at various locations, including at lysine104 (K104) and lysine147 (K147). Previous studies have suggested that K104 stabilizes helix-2/helix-3 interactions and K147 is involved in nucleotide binding. However, the impact of monoubiquitination at these residues on the overall structure, dynamics, or function of KRAS is not fully understood. In this study, we examined KRAS monoubiquitination at these sites using data from extensive (12 μs aggregate time) molecular dynamics simulations complemented by nuclear magnetic resonance spectroscopy data. We found that ubiquitin forms dynamic nonspecific interactions with various regions of KRAS and that ubiquitination at both sites modulates conformational fluctuations. In both cases, ubiquitin samples a broad range of conformational space and does not form long-lasting noncovalent contacts with KRAS but it adopts several preferred orientations relative to KRAS. To examine the functional impact of these preferred orientations, we performed a systematic comparison of the dominant configurations of the ubiquitin/KRAS simulated complex with experimental structures of KRAS bound to regulatory and effector proteins as well as a model membrane. Results from these analyses suggest that conformational selection and population shift may minimize the deleterious effects of KRAS ubiquitination at K104 and K147 on binding to some but not all interaction partners. Our findings thus provide new insights into the steric effects of ubiquitin and suggest a potential avenue for therapeutic targeting.

Biophysical and Structural Characterization of Novel RAS-Binding Domains (RBDs) of PI3Kα and PI3Kγ

Phosphatidylinositol-3-kinases (PI3Ks) are lipid kinases that phosphorylate phosphatidylinositol 4,5-bisphosphate to generate a key lipid second messenger, phosphatidylinositol 3,4,5-bisphosphate. PI3Kα and PI3Kγ require activation by RAS proteins to stimulate signaling pathways that control cellular growth, differentiation, motility and survival. Intriguingly, RAS binding to PI3K isoforms likely differ, as RAS mutations have been identified that discriminate between PI3Kα and PI3Kγ, consistent with low sequence homology (23%) between their RAS binding domains (RBDs). As disruption of the RAS/PI3Kα interaction reduces tumor growth in mice with RAS- and epidermal growth factor receptor driven skin and lung cancers, compounds that interfere with this key interaction may prove useful as anti-cancer agents. However, a structure of PI3Kα bound to RAS is lacking, limiting drug discovery efforts. Expression of full-length PI3K isoforms in insect cells has resulted in low yield and variable activity, limiting biophysical and structural studies of RAS/PI3K interactions. This led us to generate the first RBDs from PI3Kα and PI3Kγ that can be expressed at high yield in bacteria and bind to RAS with similar affinity to full-length PI3K. We also solved a 2.31 Å X-ray crystal structure of the PI3Kα-RBD, which aligns well to full-length PI3Kα. Structural differences between the PI3Kα and PI3Kγ RBDs are consistent with differences in thermal stability and may underly differential RAS recognition and RAS-mediated PI3K activation. These high expression, functional PI3K RBDs will aid in interrogating RAS interactions and could aid in identifying inhibitors of this key interaction.

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.