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

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.

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.

KRAS silencing alters chromatin physical organization and transcriptional activity in colorectal cancer cells

Clinical data revealed that KRAS mutant tumors, while initially sensitive to treatment, rapidly bypass KRAS dependence to acquire a drug-tolerant phenotype. However, the mechanisms underlying the transition from a drug-sensitive to a drug-tolerant state still elude us. Here, we show that global chromatin reorganization is a recurrent and specific feature of KRAS-dependent cells that tolerated KRAS silencing. We show that KRAS-dependent cells undergo G0/G1 cell cycle arrest after KRAS silencing, presenting a transcriptomic signature of quiescence. Proteomic analysis showed upregulated chromatin-associated proteins and transcription-associated biological processes. Accordingly, these cells shifted euchromatin/heterochromatin states, gained topologically associating domains, and altered the nanoscale physical organization of chromatin, more precisely by downregulating chromatin packing domains, a feature associated with the induction of quiescence. In addition, they also accumulated transcriptional alterations over time leading to a diversification of biological processes, linking chromatin alterations to transcriptional performance. Overall, our observations pinpoint a novel molecular mechanism of tolerance to KRAS oncogenic loss driven not by specific gene alterations but by global reorganization of genomic information, in which cells transition chromatin domain structure towards a more quiescent state and gain transcriptional reprogramming capacity.

The switch states of the GDP-bound HRAS affected by point mutations: a study from Gaussian accelerated molecular dynamics simulations and free energy landscapes

Point mutations play a vital role in the conformational transformation of HRAS. In this work, Gaussian accelerated molecular dynamics (GaMD) simulations followed by constructions of free energy landscapes (FELs) were adopted to explore the effect of mutations D33K, A59T and L120A on conformation states of the GDP-bound HRAS. The results from the post-processing analyses on GaMD trajectories suggest that mutations alter the flexibility and motion modes of the switch domains from HRAS. The analyses from FELs show that mutations induce more disordered states of the switch domains and affect interactions of GDP with HRAS, implying that mutations yield a vital effect on the binding of HRAS to effectors. The GDP-residue interaction network revealed by our current work indicates that salt bridges and hydrogen bonding interactions (HBIs) play key roles in the binding of GDP to HRAS. Furthermore, instability in the interactions of magnesium ions and GDP with the switch SI leads to the extreme disorder of the switch domains. This study is expected to provide the energetic basis and molecular mechanism for further understanding the function of HRAS.Communicated by Ramaswamy H. Sarma.

FLEXR GUI: a graphical user interface for multi-conformer modeling of proteins

Proteins are well known 'shapeshifters' which change conformation to function. In crystallography, multiple conformational states are often present within the crystal and the resulting electron-density map. Yet, explicitly incorporating alternative states into models to disentangle multi-conformer ensembles is challenging. We previously reported the tool FLEXR, which, within a few minutes, automatically separates conformational signal from noise and builds the corresponding, often missing, structural features into a multi-conformer model. To make the method widely accessible for routine multi-conformer building as part of the computational toolkit for macromolecular crystallography, we present a graphical user interface (GUI) for FLEXR, designed as a plugin for Coot 1. The GUI implementation seamlessly connects FLEXR models with the existing suite of validation and modeling tools available in Coot. We envision that FLEXR will aid crystallographers by increasing access to a multi-conformer modeling method that will ultimately lead to a better representation of protein conformational heterogeneity in the Protein Data Bank. In turn, deeper insights into the protein conformational landscape may inform biology or provide new opportunities for ligand design. The code is open source and freely available on GitHub at https://github.com/TheFischerLab/FLEXR-GUI.

Mutations influence the conformational dynamics of the GDP/KRAS complex

Mutations near allosteric sites can have a significant impact on the function of KRAS. Three specific mutations, K104Q, G12D/K104Q, and G12D/G75A, which are located near allosteric positions, were selected to investigate the molecular mechanisms behind mutation-induced influences on the activity of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations followed by the principal component analysis (PCA) were performed to improve the sampling of conformational states. The results revealed that these mutations significantly alter the structural flexibility, correlated motions, and dynamic behavior of the switch regions that are essential for KRAS binding to effectors or regulators. Furthermore, the mutations have a significant impact on the hydrogen bonding interactions between GDP and the switch regions, as well as on the electrostatic interactions of magnesium ions (Mg2+) with these regions. Our results verified that these mutations strongly influence the binding of KRAS to its effectors or regulators and allosterically regulate the activity. We believe that this work can provide valuable theoretical insights into a deeper understanding of KRAS function.Communicated by Ramaswamy H. Sarma.

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.

Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics

The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.

Structure-Based Design and Evaluation of Reversible KRAS G13D Inhibitors

Oncogenic KRAS mutations were identified decades ago, yet the selective inhibition of specific KRAS mutant proteins represents an ongoing challenge. Recent progress has been made in targeting certain P-loop mutant proteins, in particular KRAS G12C, for which the covalent inhibition of the GDP state via the Switch II pocket is now a clinically validated strategy. Inhibition of other KRAS mutant proteins such as KRAS G13D, on the other hand, still requires clinical validation. The remoteness of the D13 residue relative to the Switch II pocket in combination with the solvent exposure and conformational flexibility of the D13 side chain, as well as the difficulties of targeting carboxylate residues covalently, renders this specific protein particularly challenging to target selectively. In this report, we describe the design and evaluation of potent and KRAS G13D-selective reversible inhibitors. Subnanomolar binding to the GDP state Switch II pocket and biochemical selectivity over WT KRAS are achieved by leveraging a salt bridge with D13.

Effects of KRAS Genetic Interactions on Outcomes in Cancers of the Lung, Pancreas, and Colorectum

BACKGROUND

KRAS is among the most commonly mutated oncogenes in cancer, and previous studies have shown associations with survival in many cancer contexts. Evidence from both clinical observations and mouse experiments further suggests that these associations are allele- and tissue-specific. These findings motivate using clinical data to understand gene interactions and clinical covariates within different alleles and tissues.

METHODS

We analyze genomic and clinical data from the AACR Project GENIE Biopharma Collaborative for samples from lung, colorectal, and pancreatic cancers. For each of these cancer types, we report epidemiological associations for different KRAS alleles, apply principal component analysis (PCA) to discover groups of genes co-mutated with KRAS, and identify distinct clusters of patient profiles with implications for survival.

RESULTS

KRAS mutations were associated with inferior survival in lung, colon, and pancreas, although the specific mutations implicated varied by disease. Tissue- and allele-specific associations with smoking, sex, age, and race were found. Tissue-specific genetic interactions with KRAS were identified by PCA, which were clustered to produce five, four, and two patient profiles in lung, colon, and pancreas. Membership in these profiles was associated with survival in all three cancer types.

CONCLUSIONS

KRAS mutations have tissue- and allele-specific associations with inferior survival, clinical covariates, and genetic interactions.

IMPACT

Our results provide greater insight into the tissue- and allele-specific associations with KRAS mutations and identify clusters of patients that are associated with survival and clinical attributes from combinations of genetic interactions with KRAS mutations.

Next batter up! Targeting cancers with KRAS-G12D mutations

KRAS is the most frequently mutated oncogene in cancer. Activating mutations in codon 12, especially G12D, have the highest prevalence across a range of carcinomas and adenocarcinomas. With inhibitors to KRAS-G12D now entering clinical trials, understanding the biology of KRAS-G12D cancers, and identifying biomarkers that predict therapeutic response is crucial. In this Review, we discuss the genomics and biology of KRAS-G12D adenocarcinomas, including histological features, transcriptional landscape, the immune microenvironment, and how these factors influence response to therapy. Moreover, we explore potential therapeutic strategies using novel G12D inhibitors, leveraging knowledge gained from clinical trials using G12C inhibitors.

Beyond structural bioinformatics for genomics with dynamics characterization of an expanded KRAS mutational landscape

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations, which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

Allosteric Regulation of Switch-II Domain Controls KRAS Oncogenicity

RAS proteins are GTPases that regulate a wide range of cellular processes. RAS activity is dependent on its nucleotide-binding status, which is modulated by guanine nucleotide exchange factors (GEF) and GTPase-activating proteins (GAP). KRAS can be acetylated at lysine 104 (K104), and an acetylation-mimetic mutation of K104 to glutamine (K104Q) attenuates the in vitro-transforming capacity of oncogenic KRAS by interrupting GEF-induced nucleotide exchange. To assess the effect of this mutation in vivo, we used CRISPR-Cas9 to generate mouse models carrying the K104Q point mutation in wild-type and conditional KrasLSL-G12D alleles. Homozygous animals for K104Q were viable, fertile, and arose at the expected Mendelian frequency, indicating that K104Q is not a complete loss-of-function mutation. Consistent with our previous findings from in vitro studies, however, the oncogenic activity of KRASG12D was significantly attenuated by mutation at K104. Biochemical and structural analysis indicated that the G12D and K104Q mutations cooperate to suppress GEF-mediated nucleotide exchange, explaining the preferential effect of K104Q on oncogenic KRAS. Furthermore, K104 functioned in an allosteric network with M72, R73, and G75 on the α2 helix of the switch-II region. Intriguingly, point mutation of glycine 75 to alanine (G75A) also showed a strong negative regulatory effect on KRASG12D. These data demonstrate that lysine at position 104 is critical for the full oncogenic activity of mutant KRAS and suggest that modulating the sites in its allosteric network may provide a unique therapeutic approach in cancers expressing mutant KRAS.

SIGNIFICANCE

An allosteric network formed by interaction between lysine 104 and residues in the switch-II domain is required for KRAS oncogenicity, which could be exploited for developing inhibitors of the activated oncoprotein.

Reduced dynamic complexity allows structure elucidation of an excited state of KRASG13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS^G13D, which suggests an exchange between two conformational states in solution. Methyl and ³¹P NMR spectra of active KRAS^G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major P^γ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS^G13D and KRAS^G13D-RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS^G13D, showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

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.

Beyond Structural Bioinformatics for Genomics with Dynamics Characterization of an Expanded KRAS Mutational Landscape

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations

G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.

Tumor heterogeneity: preclinical models, emerging technologies, and future applications

Heterogeneity describes the differences among cancer cells within and between tumors. It refers to cancer cells describing variations in morphology, transcriptional profiles, metabolism, and metastatic potential. More recently, the field has included the characterization of the tumor immune microenvironment and the depiction of the dynamics underlying the cellular interactions promoting the tumor ecosystem evolution. Heterogeneity has been found in most tumors representing one of the most challenging behaviors in cancer ecosystems. As one of the critical factors impairing the long-term efficacy of solid tumor therapy, heterogeneity leads to tumor resistance, more aggressive metastasizing, and recurrence. We review the role of the main models and the emerging single-cell and spatial genomic technologies in our understanding of tumor heterogeneity, its contribution to lethal cancer outcomes, and the physiological challenges to consider in designing cancer therapies. We highlight how tumor cells dynamically evolve because of the interactions within the tumor immune microenvironment and how to leverage this to unleash immune recognition through immunotherapy. A multidisciplinary approach grounded in novel bioinformatic and computational tools will allow reaching the integrated, multilayered knowledge of tumor heterogeneity required to implement personalized, more efficient therapies urgently required for cancer patients.

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.

Comprehensive pan-cancer genomic landscape of KRAS altered cancers and real-world outcomes in solid tumors

Recent clinical development of KRAS inhibitors has heightened interest in the genomic landscape of KRAS-altered cancers. We performed a pan-cancer analysis of KRAS-altered samples from 426,706 adult patients with solid or hematologic malignancies using comprehensive genomic profiling; additional analyses included 62,369 liquid biopsy and 7241 pediatric samples. 23% of adult pan-cancer samples had KRAS alterations; 88% were mutations, most commonly G12D/G12V/G12C/G13D/G12R, and prevalence was similar in liquid biopsies. Co-alteration landscapes were largely similar across KRAS mutations but distinct from KRAS wild-type, though differences were observed in some tumor types for tumor mutational burden, PD-L1 expression, microsatellite instability, and other mutational signatures. Prognosis of KRAS-mutant versus other genomic cohorts of lung, pancreatic, and colorectal cancer were assessed using a real-world clinicogenomic database. As specific KRAS inhibitors and combination therapeutic strategies are being developed, genomic profiling to understand co-alterations and other biomarkers that may modulate response to targeted or immunotherapies will be imperative.

Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations

The conformational changes in switch domains significantly affect the activity of NRAS. Gaussian-accelerated molecular dynamics (GaMD) simulations of three separate replicas were performed to decipher the effects of G13D, Q16R, and C118S on the conformational transformation of the GDP-bound NRAS. The analyses of root-mean-square fluctuations and dynamics cross-correlation maps indicated that the structural flexibility and motion modes of the switch domains involved in the binding of NRAS to effectors are highly altered by the G13D, Q61R, and C118Smutations. The free energy landscapes (FELs) suggested that mutations induce more energetic states in NRAS than the GDP-bound WT NRAS and lead to high disorder in the switch domains. The FELs also indicated that the different numbers of sodium ions entering the GDP binding regions compensate for the changes in electrostatic environments caused by mutations, especially for G13D. The GDP-residue interactions revealed that the disorder in the switch domains was attributable to the unstable hydrogen bonds between GDP and two residues, V29 and D30. This work is expected to provide information on the energetic basis and dynamics of conformational changes in switch domains that can aid in deeply understanding the target roles of NRAS in anticancer treatment.

Heterogeneity in RAS mutations: One size does not fit all

Although oncogenic driver mutations in RAS occur in 20% of cancers, heterogeneity in the biologic outputs of different RAS mutants has hampered efforts to develop effective treatments for RAS-mutated cancers. In this issue of Science Signaling, Huynh et al. show that even among KRAS^Q61 mutants, the specific amino acid that is substituted substantially affects mutant KRAS biologic activity and oncogenicity.

Functional and biological heterogeneity of KRASQ61 mutations

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

Targeting KRAS mutant cancers: from druggable therapy to drug resistance

Kirsten Rat Sarcoma Viral Oncogene Homolog (KRAS) is the most frequently mutated oncogene, occurring in a variety of tumor types. Targeting KRAS mutations with drugs is challenging because KRAS is considered undruggable due to the lack of classic drug binding sites. Over the past 40 years, great efforts have been made to explore routes for indirect targeting of KRAS mutant cancers, including KRAS expression, processing, upstream regulators, or downstream effectors. With the advent of KRAS (G12C) inhibitors, KRAS mutations are now druggable. Despite such inhibitors showing remarkable clinical responses, resistance to monotherapy of KRAS inhibitors is eventually developed. Significant progress has been made in understanding the mechanisms of drug resistance to KRAS-mutant inhibitors. Here we review the most recent advances in therapeutic approaches and resistance mechanisms targeting KRAS mutations and discuss opportunities for combination therapy.

Chemical acylation of an acquired serine suppresses oncogenic signaling of K-Ras(G12S)

Drugs that directly impede the function of driver oncogenes offer exceptional efficacy and a therapeutic window. The recently approved mutant selective small-molecule cysteine-reactive covalent inhibitor of the G12C mutant of K-Ras, sotorasib, provides a case in point. KRAS is the most frequently mutated proto-oncogene in human cancer, yet despite success targeting the G12C allele, targeted therapy for other hotspot mutants of KRAS has not been described. Here we report the discovery of small molecules that covalently target a G12S somatic mutation in K-Ras and suppress its oncogenic signaling. We show that these molecules are active in cells expressing K-Ras(G12S) but spare the wild-type protein. Our results provide a path to targeting a second somatic mutation in the oncogene KRAS by overcoming the weak nucleophilicity of an acquired serine residue. The chemistry we describe may serve as a basis for the selective targeting of other unactivated serines.

Identification of potential inhibitors, conformational dynamics, and mechanistic insights into mutant Kirsten rat sarcoma virus (G13D) driven cancers

The mutations at the hotspot region of K-Ras result in the progression of cancer types. Our study aimed to explore the small molecule inhibitors against the G13D mutant K-Ras model with anti-cancerous activity from food and drug administration (FDA)-approved drug compounds. We implemented several computational strategies such as pharmacophore-based virtual screening, molecular docking, absorption, distribution, metabolism and excretion features, and molecular simulation to ensure the identified hit compounds have potential efficacy against G13D K-Ras. We found that the FDA-approved compounds, namely, azelastine, dihydrocodeine, paroxetine, and tramadol, are potential candidates to inhibit the action of G13D mutant K-Ras. All four compounds exhibited similar binding patterns of sotorasib, and a structural binding mechanism with significant hydrophobic contacts. The descriptor features from the QikProp of all four compounds are within allowable limits compared to sotorasib drug. Consequently, a molecular simulation result emphasized that the dihydrocodeine and tramadol exhibited less fluctuation, minimal basin, significant h-bonds, and potent inhibition against G13D K-Ras. As a result, the current research identifies prospective K-Ras inhibitors that could be further improved with biochemical analysis for precision medicine against K-Ras-driven cancers.

Thermal Shift Assay for Small GTPase Stability Screening: Evaluation and Suitability

Thermal unfolding methods are commonly used as a predictive technique by tracking the protein's physical properties. Inherent protein thermal stability and unfolding profiles of biotherapeutics can help to screen or study potential drugs and to find stabilizing or destabilizing conditions. Differential scanning calorimetry (DSC) is a 'Gold Standard' for thermal stability assays (TSA), but there are also a multitude of other methodologies, such as differential scanning fluorimetry (DSF). The use of an external probe increases the assay throughput, making it more suitable for screening studies, but the current methodologies suffer from relatively low sensitivity. While DSF is an effective tool for screening, interpretation and comparison of the results is often complicated. To overcome these challenges, we compared three thermal stability probes in small GTPase stability studies: SYPRO Orange, 8-anilino-1-naphthalenesulfonic acid (ANS), and the Protein-Probe. We studied mainly KRAS, as a proof of principle to obtain biochemical knowledge through TSA profiles. We showed that the Protein-Probe can work at lower concentration than the other dyes, and its sensitivity enables effective studies with non-covalent and covalent drugs at the nanomolar level. Using examples, we describe the parameters, which must be taken into account when characterizing the effect of drug candidates, of both small molecules and Designed Ankyrin Repeat Proteins.

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.

CRISPR in cancer biology and therapy

Over the past decade, CRISPR has become as much a verb as it is an acronym, transforming biomedical research and providing entirely new approaches for dissecting all facets of cell biology. In cancer research, CRISPR and related tools have offered a window into previously intractable problems in our understanding of cancer genetics, the noncoding genome and tumour heterogeneity, and provided new insights into therapeutic vulnerabilities. Here, we review the progress made in the development of CRISPR systems as a tool to study cancer, and the emerging adaptation of these technologies to improve diagnosis and treatment.

Classification of KRAS-Activating Mutations and the Implications for Therapeutic Intervention

Members of the family of RAS proto-oncogenes, discovered just over 40 years ago, were among the first cancer-initiating genes to be discovered. Of the three RAS family members, KRAS is the most frequently mutated in human cancers. Despite intensive biological and biochemical study of RAS proteins over the past four decades, we are only now starting to devise therapeutic strategies to target their oncogenic properties. Here, we highlight the distinct biochemical properties of common and rare KRAS alleles, enabling their classification into functional subtypes. We also discuss the implications of this functional classification for potential therapeutic avenues targeting mutant subtypes.

SIGNIFICANCE

Efforts in the recent past to inhibit KRAS oncogenicity have focused on kinases that function in downstream signal transduction cascades, although preclinical successes have not translated to patients with KRAS-mutant cancer. Recently, clinically effective covalent inhibitors of KRASG12C have been developed, establishing two principles that form a foundation for future efforts. First, KRAS is druggable. Second, each mutant form of KRAS is likely to have properties that make it uniquely druggable.

Regulation of GTPase function by autophosphorylation

A unifying feature of the RAS superfamily is a conserved GTPase cycle by which these proteins transition between active and inactive states. We demonstrate that autophosphorylation of some GTPases is an intrinsic regulatory mechanism that reduces nucleotide hydrolysis and enhances nucleotide exchange, altering the on/off switch that forms the basis for their signaling functions. Using X-ray crystallography, nuclear magnetic resonance spectroscopy, binding assays, and molecular dynamics on autophosphorylated mutants of H-RAS and K-RAS, we show that phosphoryl transfer from GTP requires dynamic movement of the switch II region and that autophosphorylation promotes nucleotide exchange by opening the active site and extracting the stabilizing Mg2+. Finally, we demonstrate that autophosphorylated K-RAS exhibits altered effector interactions, including a reduced affinity for RAF proteins in mammalian cells. Thus, autophosphorylation leads to altered active site dynamics and effector interaction properties, creating a pool of GTPases that are functionally distinct from their non-phosphorylated counterparts.

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.

Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.

Targeting RAS phosphorylation in cancer therapy: Mechanisms and modulators

RAS, a member of the small GTPase family, functions as a binary switch by shifting between inactive GDP-loaded and active GTP-loaded state. RAS gain-of-function mutations are one of the leading causes in human oncogenesis, accounting for ∼19% of the global cancer burden. As a well-recognized target in malignancy, RAS has been intensively studied in the past decades. Despite the sustained efforts, many failures occurred in the earlier exploration and resulted in an ‘undruggable’ feature of RAS proteins. Phosphorylation at several residues has been recently determined as regulators for wild-type and mutated RAS proteins. Therefore, the development of RAS inhibitors directly targeting the RAS mutants or towards upstream regulatory kinases supplies a novel direction for tackling the anti-RAS difficulties. A better understanding of RAS phosphorylation can contribute to future therapeutic strategies. In this review, we comprehensively summarized the current advances in RAS phosphorylation and provided mechanistic insights into the signaling transduction of associated pathways. Importantly, the preclinical and clinical success in developing anti-RAS drugs targeting the upstream kinases and potential directions of harnessing allostery to target RAS phosphorylation sites were also discussed.

Conformational transformation of switch domains in GDP/K-Ras induced by G13 mutants: An investigation through Gaussian accelerated molecular dynamics simulations and principal component analysis

Mutations in K-Ras are involved in a large number of all human cancers, thus, K-Ras is regarded as a promising target for anticancer drug design. Understanding the target roles of K-Ras is important for providing insights on the molecular mechanism underlying the conformational transformation of the switch domains in K-Ras due to mutations. In this study, multiple replica Gaussian accelerated molecular (MR-GaMD) simulations and principal component analysis (PCA) were applied to probe the effect of G13A, G13D and G13I mutations on conformational transformations of the switch domains in GDP-associated K-Ras. The results suggest that G13A, G13D and G13I enhance the structural flexibility of the switch domains, change the correlated motion modes of the switch domains and strengthen the total motion strength of K-Ras compared with the wild-type (WT) K-Ras. Free energy landscape analyses not only show that the switch domains of the GDP-bound inactive K-Ras mainly exist as a closed state but also indicate that mutations evidently alter the free energy profile of K-Ras and affect the conformational transformation of the switch domains between the closed and open states. Analyses of hydrophobic interaction contacts and hydrogen bonding interactions show that the mutations scarcely change the interaction network of GDP with K-Ras and only disturb the interaction of GDP with the switch (SW1). In summary, two newly introduced mutations, G13A and G13I, play similar adjustment roles in the conformational transformations of two switch domains to G13D and are possibly utilized to tune the activity of K-Ras and the binding of guanine nucleotide exchange factors.

Mathematical Modeling to Study KRAS Mutant-Specific Responses to Pathway Inhibition

This chapter will describe how mathematical modeling allows the RAS pathway to be studied with computational experiments. The mathematical model utilized simulates the biochemical reactions that regulate RAS signaling. This type of model incorporates knowledge of reaction mechanisms, including measured quantitative parameters that characterize these reactions for both wild-type and mutant RAS proteins. For an illustrative example, this chapter focuses on how modeling provided new insights that helped solve a problem that challenged the RAS community for nearly a decade: why do colorectal cancers with the KRAS G13D mutation, but not the other common KRAS mutations, benefit from EGFR inhibition? The methods described include computational dose-response experiments and the use of "computational chimeric" RAS mutants.

Parallel Rap1>RalGEF>Ral and Ras signals sculpt the C. elegans nervous system

Ras is the most commonly mutated oncogene in humans and uses three oncogenic effectors: Raf, PI3K, and RalGEF activation of Ral. Understanding the importance of RalGEF>Ral signaling in cancer is hampered by the paucity of knowledge about their function in animal development, particularly in cell movements. We found that mutations that disrupt function of RalGEF or Ral enhance migration phenotypes of mutants for genes with established roles in cell migration. We used as a model the migration of the canal associated neurons (CANs), and validated our results in HSN cell migration, neurite guidance, and general animal locomotion. These functions of RalGEF and Ral are specific to their control of Ral signaling output rather than other published functions of these proteins. In this capacity Ral functions cell autonomously as a permissive developmental signal. In contrast, we observed Ras, the canonical activator of RalGEF>Ral signaling in cancer, to function as an instructive signal. Furthermore, we unexpectedly identified a function for the close Ras relative, Rap1, consistent with activation of RalGEF>Ral. These studies define functions of RalGEF>Ral, Rap1 and Ras signaling in morphogenetic processes that fashion the nervous system. We have also defined a model for studying how small GTPases partner with downstream effectors. Taken together, this analysis defines novel molecules and relationships in signaling networks that control cell movements during development of the nervous system.

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.

Anticancer Activity of Triazolo-Thiadiazole Derivatives and Inhibition of AKT1 and AKT2 Activation

The fusion of 1,2,4-triazole and 1,3,4-thiadiazole rings results in a class of heterocycles compounds with an extensive range of pharmacological properties. A series of 1,2,4-triazolo[3,4-b]-1,2,4-thiadiazoles was synthesized and tested for its enzyme inhibition potential and anticancer activity. The results show that 1,2,4-triazolo[3,4-b]-1,2,4-thiadiazoles display potent anticancer properties in vitro against a panel of cancer cells and in vivo efficacy in HT-29 human colon tumor xenograft in CB17 severe combined immunodeficient (SCID) mice. Preliminary mechanistic studies revealed that KA25 and KA39 exhibit time- and concentration-dependent inhibition of Akt Ser-473 phosphorylation. Molecular modeling experiments indicated that 1,2,4-triazolo[3,4-b]-1,2,4-thiadiazoles bind well to the ATP binding site in Akt1 and Akt2. The low acute toxicity combined with in vitro and in vivo anticancer activity render triazolo[3,4-b]thiadiazoles KA25, KA26, and KA39 promising cancer therapeutic agents.

Molecular assemblies of the catalytic domain of SOS with KRas and oncogenic mutants

Ras is regulated by a specific guanine nucleotide exchange factor Son of Sevenless (SOS), which facilitates the exchange of inactive, GDP-bound Ras with GTP. The catalytic activity of SOS is also allosterically modulated by an active Ras (Ras-GTP). However, it remains poorly understood how oncogenic Ras mutants interact with SOS and modulate its activity. Here, native ion mobility-mass spectrometry is employed to monitor the assembly of the catalytic domain of SOS (SOScat) with KRas and three cancer-associated mutants (G12C, G13D, and Q61H), leading to the discovery of different molecular assemblies and distinct conformers of SOScat engaging KRas. We also find KRasG13D exhibits high affinity for SOScat and is a potent allosteric modulator of its activity. A structure of the KRasG13D•SOScat complex was determined using cryogenic electron microscopy providing insight into the enhanced affinity of the mutant protein. In addition, we find that KRasG13D-GTP can allosterically increase the nucleotide exchange rate of KRas at the active site more than twofold compared to KRas-GTP. Furthermore, small-molecule Ras•SOS disruptors fail to dissociate KRasG13D•SOScat complexes, underscoring the need for more potent disruptors. Taken together, a better understanding of the interaction between oncogenic Ras mutants and SOS will provide avenues for improved therapeutic interventions.

The origins and genetic interactions of KRAS mutations are allele- and tissue-specific

Mutational activation of KRAS promotes the initiation and progression of cancers, especially in the colorectum, pancreas, lung, and blood plasma, with varying prevalence of specific activating missense mutations. Although epidemiological studies connect specific alleles to clinical outcomes, the mechanisms underlying the distinct clinical characteristics of mutant KRAS alleles are unclear. Here, we analyze 13,492 samples from these four tumor types to examine allele- and tissue-specific genetic properties associated with oncogenic KRAS mutations. The prevalence of known mutagenic mechanisms partially explains the observed spectrum of KRAS activating mutations. However, there are substantial differences between the observed and predicted frequencies for many alleles, suggesting that biological selection underlies the tissue-specific frequencies of mutant alleles. Consistent with experimental studies that have identified distinct signaling properties associated with each mutant form of KRAS, our genetic analysis reveals that each KRAS allele is associated with a distinct tissue-specific comutation network. Moreover, we identify tissue-specific genetic dependencies associated with specific mutant KRAS alleles. Overall, this analysis demonstrates that the genetic interactions of oncogenic KRAS mutations are allele- and tissue-specific, underscoring the complexity that drives their clinical consequences.

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.

The Multi-Level Mechanism of Action of a Pan-Ras Inhibitor Explains its Antiproliferative Activity on Cetuximab-Resistant Cancer Cells

Ras oncoproteins play a crucial role in the onset, maintenance, and progression of the most common and deadly human cancers. Despite extensive research efforts, only a few mutant-specific Ras inhibitors have been reported. We show that cmp4–previously identified as a water-soluble Ras inhibitor– targets multiple steps in the activation and downstream signaling of different Ras mutants and isoforms. Binding of this pan-Ras inhibitor to an extended Switch II pocket on HRas and KRas proteins induces a conformational change that down-regulates intrinsic and GEF-mediated nucleotide dissociation and exchange and effector binding. A mathematical model of the Ras activation cycle predicts that the inhibitor severely reduces the proliferation of different Ras-driven cancer cells, effectively cooperating with Cetuximab to reduce proliferation even of Cetuximab-resistant cancer cell lines. Experimental data confirm the model prediction, indicating that the pan-Ras inhibitor is an appropriate candidate for medicinal chemistry efforts tailored at improving its currently unsatisfactory affinity.

KRAS or BRAF mutations cause hepatic vascular cavernomas treatable with MAP2K-MAPK1 inhibition

Human hepatic vascular cavernomas, the most common benign tumor of the liver, were described in the mid-1800s, yet the mechanisms for their formation and effective treatments remain unknown. Here, we demonstrate gain-of-function mutations in KRAS or BRAF genes within liver endothelial cells as a causal mechanism for hepatic vascular cavernomas. We identified gain-of-function mutations in KRAS or BRAF genes in pathological liver tissue samples from patients with hepatic vascular cavernomas. Mice expressing these human KRAS^G12D or BRAF^V600E mutations in hepatic endothelial cells recapitulated the human hepatic vascular cavernoma phenotype of dilated sinusoidal capillaries with defective branching patterns. KRAS^G12D or BRAF^V600E induced “zipper-like” contiguous expression of junctional proteins at sinusoidal endothelial cell–cell contacts, switching capillaries from branching to cavernous expansion. Pharmacological or genetic inhibition of the endothelial RAS–MAPK1 signaling pathway rescued hepatic vascular cavernoma formation in endothelial KRAS^G12D- or BRAF^V600E-expressing mice. These results uncover a major cause of hepatic vascular cavernomas and provide a road map for their personalized treatment.

Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics

Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.

Structural basis of ubiquitination mediated by protein splicing in early Eukarya

BACKGROUND

Inteins are intervening proteins, which are known to perform protein splicing. The reaction results in the production of an intein domain and an inteinless protein, which shows no trace of the insertion. BIL2 is part of the polyubiquitin locus of Tetrahymena thermophila (BUBL), where two bacterial-intein-like (BIL) domains lacking the C + 1 nucleophile, are flanked by two independent ubiquitin-like domains (ubl4/ubl5).

METHODS

We solved the X-ray structures of BIL2 in both the inactive and unprecedented, zinc-induced active, forms. Then, we characterized by mass spectrometry the BUBL splicing products in the absence and in the presence of T.thRas-GTPase. Finally, we investigated the effect of ubiquitination on T.thRas-GTPase by molecular dynamics simulations.

RESULTS

The structural analysis demonstrated that zinc-induced conformational change activates protein splicing. Moreover, mass spectrometry characterization of the splicing products shed light on the possible function of BIL2, which operates as a "single-ubiquitin-dispensing-platform", allowing the conjugation, via isopeptide bond formation (K(εNH₂)-C-ter), of ubl4 to either ubl5 or T.thRas-GTPase. Lastly, we demonstrated that T.thRas-GTPase ubiquitination occurs in proximity of the nucleotide binding pocket and stabilizes the protein active state.

CONCLUSIONS

We demonstrated that BIL2 is activated by zinc and that protein splicing induced by this intein does not take place through classical or aminolysis mechanisms but via formation of a covalent isopeptide bond, causing the ubiquitination of endogenous substrates such as T.thRas-GTPase.

GENERAL SIGNIFICANCE

In this "enzyme-free" ubiquitination mechanism the isopeptide formation, which canonically requires E1-E2-E3 enzymatic cascade and constitutes the alphabet of ubiquitin biology, is achieved in a single, concerted step without energy consumption.

Natural Products Attenuating Biosynthesis, Processing, and Activity of Ras Oncoproteins: State of the Art and Future Perspectives

RAS genes encode signaling proteins, which, in mammalian cells, act as molecular switches regulating critical cellular processes as proliferation, growth, differentiation, survival, motility, and metabolism in response to specific stimuli. Deregulation of Ras functions has a high impact on human health: gain-of-function point mutations in RAS genes are found in some developmental disorders and thirty percent of all human cancers, including the deadliest. For this reason, the pathogenic Ras variants represent important clinical targets against which to develop novel, effective, and possibly selective pharmacological inhibitors. Natural products represent a virtually unlimited resource of structurally different compounds from which one could draw on for this purpose, given the improvements in isolation and screening of active molecules from complex sources. After a summary of Ras proteins molecular and regulatory features and Ras-dependent pathways relevant for drug development, we point out the most promising inhibitory approaches, the known druggable sites of wild-type and oncogenic Ras mutants, and describe the known natural compounds capable of attenuating Ras signaling. Finally, we highlight critical issues and perspectives for the future selection of potential Ras inhibitors from natural sources.

KRAS Ubiquitination at Lysine 104 Retains Exchange Factor Regulation by Dynamically Modulating the Conformation of the Interface

RAS proteins function as highly regulated molecular switches that control cellular growth. In addition to regulatory proteins, RAS undergoes a number of posttranslational modifications (PTMs) that regulate its activity. Lysine 104, a hot spot for multiple PTMs, is a highly conserved residue that forms key interactions that stabilize the RAS helix-2(H2)/helix-3(H3) interface. Mutation at 104 attenuates interaction with guanine nucleotide exchange factors (GEFs), whereas ubiquitination at lysine 104 retains GEF regulation. To elucidate how ubiquitination modulates RAS function, we generated monoubiquitinated KRAS at 104 using chemical biology approaches and conducted biochemical, NMR, and computational analyses. We find that ubiquitination promotes a new dynamic interaction network and alters RAS conformational dynamics to retain GEF function. These findings reveal a mechanism by which ubiquitination can regulate protein function.