For nearly a decade, the difficulties associated with both the determination and reproducibility of Ras-dependency indexes (RDIs) have limited their application and further delineation of the biology underlying Ras dependency. In this report, we describe the application of a computational single sample gene set enrichment analysis (ssGSEA) method to derive RDIs with gene expression data. The computationally derived RDIs across the Cancer Cell Line Encyclopedia (CCLE) cell lines show excellent agreement with the experimentally derived values and high correlation with a previous in-house siRNA effector node (siREN) study and external studies. Using EMT signature-derived RDIs and data from cell lines representing the extremes in RAS dependency, we identified enriched pathways distinguishing these classes, including the Fas signaling pathway and a putative Ras-independent pathway first identified in NK cells. Importantly, extension of the method to patient samples from The Cancer Genome Atlas (TCGA) showed the same consensus differential expression patterns for these two pathways across multiple tissue types. Last, the computational RDIs displayed a significant association with TCGA cancer patients' survival outcomes. Together, these lines of evidence confirm that our computationally derived RDIs faithfully represent a measure of Ras dependency in both cancer cell lines and patient samples. The application of such computational RDIs can provide insights into Ras biology and potential clinical applications.

RAS proteins are binary switches, cycling between ON and OFF states during signal transduction. These switches are normally tightly controlled, but in RAS-related diseases, such as cancer, RASopathies, and many psychiatric disorders, mutations in the RAS genes or their regulators render RAS proteins persistently active. The structural basis of the switch and many of the pathways that RAS controls are well known, but the precise mechanisms by which RAS proteins function are less clear. All RAS biology occurs in membranes: a precise understanding of RAS' interaction with membranes is essential to understand RAS action and to intervene in RAS-driven diseases.

In cancer research, RAS biology has been focused on only a handful of tumor types. While RAS genes have long been suspected as common contributors to a wide spectrum of cancer types, robust evidence is required to firmly establish their critical oncogenic significance. We present a data mining study using DepMap genome-wide CRISPR screening data, which provide substantial evidence to support the prominent pervasive oncogenic role and tissue-specific permissiveness of RAS gene mutations. Differential analysis of CRISPR effect scores identifies K- or N-RAS genes as the most differential gene in contrasts of (K-, N-, combined) RAS mutant versus wild-type cell lines across multiple tissue types. The distinguished tissue-specific pattern of KRAS vs. NRAS as top differential genes in subsets of tissue types and evidence from genome data supported the idea of KRAS- and NRAS-engaged tissue types. To our knowledge, this is the first report of prominent pervasive oncogenic role of RAS mutations revealed by gene dependency data that is beyond the current understanding of the oncogenic role of RAS genes and their well-known involved tissue types. Our findings strongly support RAS mutations as primary oncogenic drivers beyond traditionally recognized cancer types and offer insights into their tissue-specific permissiveness.

The 3 human RAS genes play pivotal roles regulating proliferation, differentiation, and survival in normal cells and become mutated in 15% to 20% of all human tumors and amplified in many others. In this report, we examined data from The Cancer Genome Atlas to investigate the relationship between RAS gene mutational status and messenger RNA expression. We show that all 3 RAS genes exhibit increased expression when they are mutated in a context-dependent manner. In the case of KRAS, this increase is manifested by a larger proportional increase in KRAS4A than KRAS4B, although both increase significantly. In addition, the mutational status of RAS genes can be associated with expression changes in other RAS genes, with most of these cases showing decreased expression. The mutational status associations with expression are recapitulated in cancer cell lines. Increases in expression are mediated by both copy number variation and contextual differences, including mutational status of epidermal growth factor receptor (EGFR) and BRAF. These findings potentially reveal an adaptive response during tumor evolution that is dependent on the mutational status of proximal genes in the RAS pathway and cellular context. Cell contextual differences in these adaptations may influence therapeutic responsiveness and alternative resistance mechanisms.

RAS proteins cycling between the active-form (GTP-bound) and inactive-form (GDP-bound) play a key role in cell signaling pathways that control cell survival, proliferation, and differentiation. Mutations at codon 12, 13, and 61 in RAS are known to attenuate its GTPase activity favoring the RAS active state and constitutively active downstream signaling. This hyperactivation accounts for various malignancies including pancreatic, lung, and colorectal cancers. Active KRAS is found to exist in equilibrium between two rapidly interconverting conformational states (State1-State2) in solution. Due to this dynamic feature of the protein, the ¹H-¹⁵N correlation cross-peak signals of several amino acid (AA) residues of KRAS belonging to the flexible loop regions are absent from its 2D ¹H-¹⁵N HSQC spectrum within and near physiological solution pH. A threonine to serine mutation at position 35 (T35S) shifts the interconverting equilibrium to State1 conformation and enables the emergence of such residues in the 2D ¹H-¹⁵N HSQC spectrum due to gained conformational rigidity. We report here the ¹H^N, ¹⁵N, and ¹³C backbone resonance assignments for the 19.2 kDa (AA 1-169) protein constructs of KRAS-GppNHp harboring T35S mutation (KRAS^T35S/C118S-GppNHp) and of its oncogenic counterpart harboring the Q61L mutation (KRAS^{T35S/Q61L/C118S}-GppNHp) using heteronuclear, multidimensional NMR spectroscopy at 298 K. High resolution NMR data allowed the unambiguous assignments of ¹H-¹⁵N correlation cross-peaks for all the residues except for Met1. Furthermore, 2D ¹H-¹⁵N HSQC overlay of two proteins assisted in determination of Q61L mutation-induced chemical shift perturbations for select residues in the regions of P-loop, Switch-II, and helix α3.

Deregulation of KRAS4b signaling pathway has been implicated in 30% of all cancers. Membrane localization of KRAS4b is an essential step for the initiation of the downstream signaling cascades that guide various cellular mechanisms. KRAS4b plasma membrane (PM) binding is mediated by the insertion of a prenylated moiety that is attached to the terminal carboxy-methylated cysteine, in addition to electrostatic interactions of its positively charged hypervariable region with anionic lipids. Calmodulin (CaM) has been suggested to selectively bind KRAS4b to act as a negative regulator of the RAS/mitogen-activated protein kinase (MAPK) signaling pathway by displacing KRAS4b from the membrane. However, the mechanism by which CaM can recognize and displace KRAS4b from the membrane is not well understood. In this study, we employed biophysical and structural techniques to characterize this mechanism in detail. We show that KRAS4b prenylation is required for binding to CaM and that the hydrophobic pockets of CaM can accommodate the prenylated region of KRAS4b, which might represent a novel CaM-binding motif. Remarkably, prenylated KRAS4b forms a 2:1 stoichiometric complex with CaM in a nucleotide-independent manner. The interaction between prenylated KRAS4b and CaM is enthalpically driven, and electrostatic interactions also contribute to the formation of the complex. The prenylated KRAS4b terminal KSKTKC-farnesylation and carboxy-methylation is sufficient for binding and defines the minimal CaM-binding motif. This is the same region implicated in membrane and phosphodiesterase6-δ binding. Finally, we provide a structure-based docking model by which CaM binds to prenylated KRAS4b. Our data provide new insights into the KRAS4b-CaM interaction and suggest a possible mechanism whereby CaM can regulate KRAS4b membrane localization.

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.

RAF kinases are integral to the RAS-MAPK signaling pathway, and proper RAF1 folding relies on its interaction with the chaperone HSP90 and the cochaperone CDC37. Understanding the intricate molecular interactions governing RAF1 folding is crucial for comprehending this process. Here, we present a cryo-EM structure of the closed-state RAF1-HSP90-CDC37 complex, where the C-lobe of the RAF1 kinase domain binds to one side of the HSP90 dimer, and an unfolded N-lobe segment of the RAF1 kinase domain threads through the center of the HSP90 dimer. CDC37 binds to the kinase C-lobe, mimicking the N-lobe with its HxNI motif. We also describe structures of HSP90 dimers without RAF1 and CDC37, displaying only N-terminal and middle domains, which we term the semi-open state. Employing 1 μs atomistic simulations, energetic decomposition, and comparative structural analysis, we elucidate the dynamics and interactions within these complexes. Our quantitative analysis reveals that CDC37 bridges the HSP90-RAF1 interaction, RAF1 binds HSP90 asymmetrically, and that HSP90 structural elements engage RAF1s unfolded region. Additionally, N- and C-terminal interactions stabilize HSP90 dimers, and molecular interactions in HSP90 dimers rearrange between the closed and semi-open states. Our findings provide valuable insight into the contributions of HSP90 and CDC37 in mediating client folding.

SHOC2 acts as a strong synthetic lethal interactor with MEK inhibitors in multiple KRAS cancer cell lines. SHOC2 forms a heterotrimeric complex with MRAS and PP1C that is essential for regulating RAF and MAPK-pathway activation by dephosphorylating a specific phosphoserine on RAF kinases. Here we present the high-resolution crystal structure of the SHOC2-MRAS-PP1C (SMP) complex and apo-SHOC2. Our structures reveal that SHOC2, MRAS, and PP1C form a stable ternary complex in which all three proteins synergistically interact with each other. Our results show that dephosphorylation of RAF substrates by PP1C is enhanced upon interacting with SHOC2 and MRAS. The SMP complex forms only when MRAS is in an active state and is dependent on SHOC2 functioning as a scaffolding protein in the complex by bringing PP1C and MRAS together. Our results provide structural insights into the role of the SMP complex in RAF activation and how mutations found in Noonan syndrome enhance complex formation, and reveal new avenues for therapeutic interventions.

Farnesylation and carboxymethylation of KRAS4b (Kirsten rat sarcoma isoform 4b) are essential for its interaction with the plasma membrane where KRAS-mediated signaling events occur. Phosphodiesterase-δ (PDEδ) binds to KRAS4b and plays an important role in targeting it to cellular membranes. We solved structures of human farnesylated-methylated KRAS4b in complex with PDEδ in two different crystal forms. In these structures, the interaction is driven by the C-terminal amino acids together with the farnesylated and methylated C185 of KRAS4b that binds tightly in the central hydrophobic pocket present in PDEδ. In crystal form II, we see the full-length structure of farnesylated-methylated KRAS4b, including the hypervariable region. Crystal form I reveals structural details of farnesylated-methylated KRAS4b binding to PDEδ, and crystal form II suggests the potential binding mode of geranylgeranylated-methylated KRAS4b to PDEδ. We identified a 5-aa-long sequence motif (Lys-Ser-Lys-Thr-Lys) in KRAS4b that may enable PDEδ to bind both forms of prenylated KRAS4b. Structure and sequence analysis of various prenylated proteins that have been previously tested for binding to PDEδ provides a rationale for why some prenylated proteins, such as KRAS4a, RalA, RalB, and Rac1, do not bind to PDEδ. Comparison of all four available structures of PDEδ complexed with various prenylated proteins/peptides shows the presence of additional interactions due to a larger protein-protein interaction interface in KRAS4b-PDEδ complex. This interface might be exploited for designing an inhibitor with minimal off-target effects.