KRAS interaction with RAF1 RAS-binding domain and cysteine-rich domain provides insights into RAS-mediated RAF activation

QM/MM study reveals novel mechanism of KRAS and KRASG12R catalyzed GTP hydrolysis

As a crucial drug target, KRAS can regulate most cellular processes involving guanosine triphosphate (GTP) hydrolysis. However, the mechanism of GTP hydrolysis has remained controversial over the past decades. Here, several different GTP hydrolysis mechanisms catalyzed by wild-type KRAS (WT-KRAS) and KRASG12R mutants were discussed via four QM/MM calculation models. Based on the computational results, a Mg2+-coordinated H2O-mediated GTP hydrolysis mechanism was proposed. In this mechanism, a Mg2+-coordinated H2O first protonates the fully deprotonated GTP, and then the Mg2+ coordinated hydroxyl anion is generated. The Pγ-O bond is formed via the SN2 attack of the second H2O on the Pγ atom of the GTP, leading to the simultaneous cleavage of the Pγ-O bond. Meanwhile, the hydroxyl anion coordinated to Mg2+ and generated in the first step acts as a proton acceptor from water. This Mg2+ coordinated H2O-involved GTP hydrolysis mechanism may also be suitable for Mg2+-catalyzed ATP hydrolysis. Furthermore, the mechanism of GTP hydrolysis catalyzed by the KRASG12R mutant, whose hydrolysis rate was approximately 40-fold slower than WT-KRAS, was also discussed. Our QM/MM calculations reveal that GTP is easily protonated by the residue R12, and the energy barrier of GTP hydrolysis catalyzed by the KRASG12R mutant is lower than the corresponding barrier for WT-KRAS. Nevertheless, molecular dynamics (MD) simulations reveal that R12, a residue characterized by significant steric hindrance, is positioned at the GTP site where the nucleophilic attack by water occurs during Pγ-O bond formation, thereby strongly impeding the approach of water molecules to GTP. As a result, the GTP hydrolysis rate catalyzed by the KRASG12R mutant was severely impaired. Uncovering the GTP hydrolysis mechanism catalyzed by the WT-KRAS and KRASG12R mutant may also give a reasonable explanation for the relationship between the KRASG12R mutation and the occurrence of cancer. We hope this finding will provide useful guidance for drug discovery that targets KRAS.

Allostery in Disease: Anticancer Drugs, Pockets, and the Tumor Heterogeneity Challenge

Charting future innovations is challenging. Yet, allosteric and orthosteric anticancer drugs are undergoing a revolution and taxing unresolved dilemmas await. Among the imaginative innovations, here we discuss cereblon and thalidomide derivatives as a means of recruiting neosubstrates and their degradation, allosteric heterogeneous bifunctional drugs like PROTACs, drugging phosphatases, inducers of targeted posttranslational protein modifications, antibody-drug conjugates, exploiting membrane interactions to increase local concentration, stabilizing the folded state, and more. These couple with harnessing allosteric cryptic pockets whose discovery offers more options to modulate the affinity of orthosteric, active site inhibitors. Added to these are strategies to counter drug resistance through drug combinations co-targeting pathways to bypass signaling blockades. Here, we discuss on the molecular and cellular levels, such inspiring advances, provide examples of their applications, their mechanisms and rational. We start with an overview on difficult to target proteins and their properties-rarely, if ever-conceptualized before, discuss emerging innovative drugs, and proceed to the increasingly popular allosteric cryptic pockets-their advantages-and critically, issues to be aware of. We follow with drug resistance and in-depth discussion of tumor heterogeneity. Heterogeneity is a hallmark of highly aggressive cancers, the core of drug resistance unresolved challenge. We discuss potential ways to target heterogeneity by predicting it. The increase in experimental and clinical data, computed (cell-type specific) interactomes, capturing transient cryptic pockets, learned drug resistance, workings of regulatory mechanisms, heterogeneity, and resistance-based cell signaling drug combinations, assisted by AI-driven reasoning and recognition, couple with creative allosteric drug discovery, charting future innovations.

Pioneer in Molecular Biology: Conformational Ensembles in Molecular Recognition, Allostery, and Cell Function

My focus shifted over the years. In 1978 I developed an efficient O(n3) dynamic programming algorithm for the-then open problem of RNA secondary structure prediction. This algorithm, now dubbed the "Nussinov algorithm", "Nussinov plots", and "Nussinov diagrams", is still taught in classes across Europe and the U.S. As sequences started coming out in the 1980s, I started seeking genome-encoded functional signals, later becoming a bioinformatics trend. In the early 1990s I transited to proteins, co-developing a powerful computer vision-based docking algorithm. In the late 1990s, I proposed the foundational role of conformational ensembles in molecular recognition and allostery. At the time, conformational ensembles and free energy landscapes were viewed as physical properties of protein molecules but were not associated with function. The classical view of molecular recognition and binding was based on only two conformations captured by crystallography: open and closed. I proposed that all conformational states preexist. Proteins always have not one folded form- nor two-but many folded forms. Thus, rather than inducing fit, binding can work by shifting the ensembles between states, and this shifting, or redistributing the ensembles to maintain equilibrium, is the origin of the allosteric effect and protein, thus cell, function. This new paradigm impacted community views in allosteric drug design, catalysis, and regulation. Dynamic conformational ensemble shifts are now acknowledged as the origin of recognition, allostery, and signaling, underscoring that conformational ensembles-not proteins-are the workhorses of the cell, pioneering the fundamental idea that dynamic ensembles are the driving force behind cellular processes.

Discovery of Elironrasib (RMC-6291), a Potent and Orally Bioavailable, RAS(ON) G12C-Selective, Covalent Tricomplex Inhibitor for the Treatment of Patients with RAS G12C-Addicted Cancers

The discovery of elironrasib (RMC-6291) represents a significant breakthrough in targeting the previously deemed undruggable GTP-bound, active KRASG12C. To target the active state of RAS (RAS(ON)) directly, we have employed an innovative tri-complex inhibitor (TCI) modality involving formation of a complex with an inhibitor, the intracellular chaperone protein CypA, and the target protein KRASG12C in its GTP-bound form. The resulting tri-complex inhibits oncogenic signaling, inducing tumor regressions across various preclinical models of KRASG12C mutant human cancers. Here we report structure-guided medicinal chemistry efforts that led to the discovery of elironrasib, a potent, orally bioavailable, RAS(ON) G12C-selective, covalent, tri-complex inhibitor. The investigational agent elironrasib is currently undergoing phase 1 clinical trials (NCT05462717, NCT06128551, NCT06162221), with preliminary data indicating clinical activity in patients who had progressed on first-generation inactive state-selective KRASG12C inhibitors.

Ionic control of small GTPase HRas using calmodulin

HRas is a small GTPase that plays physiologically important roles in various intracellular signal transduction processes, such as cell growth and proliferation. The structure and action mechanisms of HRas have been well characterized, leading to its widespread use as a molecular switch in bionanomachines. Calmodulin (CaM), a calcium ion-binding protein, acts as an ion-binding molecular switch and activates the target enzymes. We previously demonstrated that the fusion protein of HRas (M13-HRas) with the CaM target peptide M13 at the N-terminus of HRas exhibits reversible regulation of GTPase activity and the interaction between M13-HRas and the downstream signalling factor Raf by calcium ions with CaM. In this study, we prepared two new HRas fusion proteins with the M13 peptide at the C-terminus (HRas-M13) and both termini (M13-HRas-M13) of HRas and analysed the calcium-dependent regulation of HRas function. M13-HRas-M13 more efficiently controlled GTPase, interaction with Raf and the HRas regulator GEF by calcium ions with CaM.

Accurate Characterization of the Allosteric Energy Landscapes, Binding Hotspots and Long-Range Communications for KRAS Complexes with Effector Proteins : Integrative Approach Using Microsecond Molecular Dynamics, Deep Mutational Scanning of Binding Energetics and Allosteric Network Modeling

KRAS is a pivotal oncoprotein that regulates cell proliferation and survival through interactions with downstream effectors such as RAF1. Oncogenic mutations in KRAS, including G12V, G13D, and Q61R, drive constitutive activation and hyperactivation of signaling pathways, contributing to cancer progression. Despite significant advances in understanding KRAS biology, the structural and dynamic mechanisms of KRAS binding and allostery by which oncogenic mutations enhance KRAS-RAF1 binding and signaling remain incompletely understood. In this study, we employ microsecond molecular dynamics simulations, Markov State Modeling, mutational scanning and binding free energy calculations together with dynamic network modeling to elucidate the effect of KRAS mutations and characterize the thermodynamic and allosteric drivers and hotspots of KRAS binding and oncogenic activation. Our simulations revealed that oncogenic mutations stabilize the open active conformation of KRAS by differentially modulating the flexibility of the switch I and switch II regions, thereby enhancing RAF1 binding affinity. The G12V mutation rigidifies both switch I and switch II, locking KRAS in a stable, active state. In contrast, the G13D mutation moderately reduces switch I flexibility while increasing switch II dynamics, restoring a balance between stability and flexibility. The Q61R mutation induces a more complex conformational landscape, characterized by the increased switch II flexibility and expansion of functional macrostates, which promotes prolonged RAF1 binding and signaling. Mutational scanning of KRAS-RAF1 complexes identified key binding affinity hotspots, including Y40, E37, D38, and D33, and together with the MM-GBSA analysis revealed the hotspots leverage synergistic electrostatic and hydrophobic binding interactions in stabilizing the KRAS-RAF1 complexes. Network-based analysis of allosteric communication identifies critical KRAS residues (e.g., L6, E37, D57, R97) that mediate long-range interactions between the KRAS core and the RAF1 binding interface. The central β-sheet of KRAS emerges as a hub for transmitting conformational changes, linking distant functional sites and facilitating allosteric regulation. Strikingly, the predicted allosteric hotspots align with experimentally identified allosteric binding hotspots that define the energy landscape of KRAS allostery. This study highlights the power of integrating computational modeling with experimental data to unravel the complex dynamics of KRAS and its mutants. The identification of binding hotspots and allosteric communication routes offers new opportunities for developing targeted therapies to disrupt KRAS-RAF1 interactions and inhibit oncogenic signaling. Our results underscore the potential of computational approaches to guide the design of allosteric inhibitors and mutant-specific therapies for KRAS-driven cancers.

Plasma membrane-associated ARAF condensates fuel RAS-related cancer drug resistance

RAF protein kinases are major RAS effectors that function by phosphorylating MEK. Although all three RAF isoforms share a conserved RAS binding domain and bind to GTP-loaded RAS, only ARAF uniquely enhances RAS activity. Here we uncovered the molecular basis of ARAF in regulating RAS activation. The disordered N-terminal sequence of ARAF drives self-assembly, forming ARAF-RAS condensates tethered to the plasma membrane. These structures concentrate active RAS locally, impeding NF1-mediated negative regulation of RAS, thereby fostering receptor tyrosine kinase (RTK)-triggered RAS activation. In RAS-mutant tumors, loss of the ARAF N terminus sensitizes tumor cells to pan-RAF inhibition. In hormone-sensitive cancers, increased ARAF condensates drive endocrine therapy resistance, whereas ARAF depletion reverses RTK-dependent resistance. Our findings delineate ARAF-RAS protein condensates as distinct subcellular structures sustaining RAS activity and facilitating oncogenic RAS signaling. Targeting ARAF-RAS condensation may offer a strategy to overcome drug resistance in both wild-type and mutant ARAF-mediated scenarios.

Structural insights into isoform-specific RAS-PI3Kα interactions and the role of RAS in PI3Kα activation

Mutations in RAS and PI3Kα are major drivers of human cancer. Their interaction plays a crucial role in activating PI3Kα and amplifying the PI3K-AKT-mTOR pathway. Disrupting RAS-PI3Kα interaction enhances survival in lung and skin cancer models and reduces tumor growth and angiogenesis, although the structural details of this interaction remain unclear. Here, we present structures of KRAS, RRAS2, and MRAS bound to the catalytic subunit (p110α) of PI3Kα, elucidating the interaction interfaces and local conformational changes upon complex formation. Structural and mutational analyses highlighted key residues in RAS and PI3Kα impacting binding affinity and revealed isoform-specific differences at the interaction interface in RAS and PI3K isoforms, providing a rationale for their differential affinities. Notably, in the RAS-p110α complex structures, RAS interaction with p110α is limited to the RAS-binding domain and does not involve the kinase domain. This study underscores the pivotal role of the RAS-PI3Kα interaction in PI3Kα activation and provides a blueprint for designing PI3Kα isoform-specific inhibitors to disrupt this interaction.

Dynamic conformational equilibria in the active states of KRAS and NRAS

The design of potent RAS inhibitors benefits from a molecular understanding of the dynamics in KRAS and NRAS and their oncogenic mutants. Here we characterize switch-1 dynamics in GTP-state KRAS and NRAS by 31P NMR, by 15N relaxation dispersion NMR, hydrogen-deuterium exchange mass spectrometry (HDX-MS), and molecular dynamics simulations. In GMPPNP-bound KRAS and NRAS, we see the co-existence of two conformational states, corresponding to an "inactive" state-1 and an "active" state-2, as previously reported. The KRAS oncogenic mutations G12D, G12C and G12V only slightly affect this equilibrium towards the "inactive" state-1, with rank order wt < G12C < G12D < G12V. In contrast, the NRAS Q61R oncogenic mutation shifts the equilibrium fully towards the "active" state-2. Our molecular dynamics simulations explain this by the observation of a transient hydrogen bond between the Arg61 side chain and the Thr35 backbone carbonyl oxygen. NMR relaxation dispersion experiments with GTP-bound KRAS Q61R confirm a drastic decrease in the population of state-1, but still detect a small residual population (1.8%) of this conformer. HDX-MS indicates that higher populations of state-1 correspond to increased hydrogen-deuterium exchange rates in some regions and increased flexibility, whereas low state-1 populations are associated with KRAS rigidification. We elucidated the mechanism of action of a potent KRAS G12D inhibitor, MRTX1133. Binding of this inhibitor to the switch-2 pocket causes a complete shift of KRAS G12D towards the "inactive" conformation and prevents binding of effector RAS-binding domain (RBD) at physiological concentrations, by signaling through an allosteric network.

Single-cell sensor analyses reveal signaling programs enabling Ras-G12C drug resistance

Clinical resistance to rat sarcoma virus (Ras)-G12C inhibitors is a challenge. A subpopulation of cancer cells has been shown to undergo genomic and transcriptional alterations to facilitate drug resistance but the immediate adaptive effects on Ras signaling in response to these drugs at the single-cell level is not well understood. Here, we used Ras biosensors to profile the activity and signaling environment of endogenous Ras at the single-cell level. We found that a subpopulation of KRas-G12C cells treated with Ras-G12C-guanosine-diphosphate inhibitors underwent adaptive signaling and metabolic changes driven by wild-type Ras at the Golgi and mutant KRas at the mitochondria, respectively. Our Ras biosensors identified major vault protein as a mediator of Ras activation through its scaffolding of Ras signaling pathway components and metabolite channels. Overall, methods including ours that facilitate direct analysis on the single-cell level can report the adaptations that subpopulations of cells adopt in response to cancer therapies, thus providing insight into drug resistance.

Structural Insights Into the Impact of the Glycine-Rich Loop Mutation in Noonan Syndrome on the ATP Binding Pocket of CRAF Kinase

The pathogenic G361A variant of CRAF, associated with increased intrinsic kinase activity in Noonan syndrome (NS), remains poorly understood in terms of its molecular and structural impact on kinase activity. To elucidate the mechanistic implications of the glycine to alanine substitution at residue 361 in CRAF, we employed molecular dynamics simulations. Our findings reveal that this mutation predominantly affects the ATP binding pocket and critical intermolecular interactions within the active cleft that favors the phosphate transfer reaction. Notably, our data highlight significant alterations in key interactions involving Lys470/Asp486 and ATP.Mg2+ in CRAFG361A that are absent in wild-type CRAF. Additionally, we identified a novel interaction mode between Lys431 and γ-phosphate in wild-type CRAF, a residue evolutionarily conserved in CRAFs but not in related kinases such as BRAF, ARAF, and KSR1/2. Furthermore, observed shifts in the αC-helix and G-loop relative to the wild-type correlate with an enlarged ATP-binding cavity in the mutant, reflecting structural adaptations due to these mutations. Overall, these structural insights underscore the elevated intrinsic kinase activity of the CRAFG361A variant and provide crucial mechanistic details that could inform the development of specific inhibitors targeting this variant.

Computationally designed sensors detect endogenous Ras activity and signaling effectors at subcellular resolution

The utility of genetically encoded biosensors for sensing the activity of signaling proteins has been hampered by a lack of strategies for matching sensor sensitivity to the physiological concentration range of the target. Here we used computational protein design to generate intracellular sensors of Ras activity (LOCKR-based Sensor for Ras activity (Ras-LOCKR-S)) and proximity labelers of the Ras signaling environment (LOCKR-based, Ras activity-dependent Proximity Labeler (Ras-LOCKR-PL)). These tools allow the detection of endogenous Ras activity and labeling of the surrounding environment at subcellular resolution. Using these sensors in human cancer cell lines, we identified Ras-interacting proteins in oncogenic EML4-Alk granules and found that Src-Associated in Mitosis 68-kDa (SAM68) protein specifically enhances Ras activity in the granules. The ability to subcellularly localize endogenous Ras activity should deepen our understanding of Ras function in health and disease and may suggest potential therapeutic strategies.

Oxidative protein damage negatively affects protein-protein interaction: The case of KRAS-cRAF

Protein-protein interactions (PPIs) play crucial roles in cellular signaling, transmitting signals from the cell surface to its interior. One of the most important signaling cascades is the RAS-RAF-MEK-ERK pathway. This pathway is initiated by various upstream signaling reactions, including receptor tyrosine kinase (RTK) activation, and it controls many biological functions like cell proliferation, differentiation, and survival. Once RAS is activated, it binds RAF and relays the signal to downstream proteins. The RAS-binding domain (RBD) in RAF protein plays a crucial role in this process, facilitating the RAS-ERK pathway signaling. In this study, we explored the effect of oxidative stress induced by UV radiation on the KRAS-RBD interaction. Using the Split Intein-Mediated Protein Ligation (SIMPL) method, we assessed the impact of different UV doses on KRAS-RBD interactions and observed a disruption of this interaction at higher doses. UV-treated samples exhibited high levels of protein carbonylation, as detected by Oxime Blot and mass spectrometry (MS) analysis, indicating oxidative damage. The MS results provided detailed insights into specific carbonylation modifications on the KRAS protein. Our study demonstrates that protein oxidation and carbonylation can disrupt protein-protein interactions, specifically the KRAS/c-RAF interaction. These findings highlight the impact of oxidative stress on signaling pathways, such as those triggered by UV irradiation. A deeper understanding of these molecular changes may aid in developing therapies targeting diseases linked to oxidative stress, including cancer.

Motif-guided identification of KRAS-interacting proteins

BACKGROUND

For decades, KRAS has always been a huge challenge to the field of drug discovery for its significance in cancer progression as well as its difficulties in being targeted as an "undruggable" protein. KRAS regulates downstream signaling pathways through protein-protein interactions, whereas many interaction partners of KRAS remain unknown.

RESULTS

We developed a workflow to computationally predict and experimentally validate the potential KRAS-interacting proteins based on the interaction mode of KRAS and its known binding partners. We extracted 17 KRAS-interacting motifs from all experimentally determined KRAS-containing protein complexes as queries to identify proteins containing fragments structurally similar to the queries in the human protein structure database using our in-house protein-protein interaction prediction method, PPI-Miner. Finally, out of the 78 predicted potential interacting proteins of KRAS, 10 were selected for experimental validation, including BRAF, a previously reported interacting protein, which served as the positive control in our validation experiments. Additionally, a known peptide that binds to KRAS, KRpep-2d, was also used as a positive control. The predicted interacting motifs of these 10 proteins were synthesized to perform biolayer interferometry assays, with 4 out of 10 exhibiting binding affinities to KRAS, and the strongest, GRB10, was selected for further validation. Additionally, the interaction between GRB10 (RA-PH domain) and KRAS was confirmed via immunofluorescence and co-immunoprecipitation.

CONCLUSIONS

These results demonstrate the effectiveness of our workflow in predicting potential interacting proteins for KRAS and deepen the understanding of KRAS-driven tumor mechanisms and the development of therapeutic strategies.

The Raf/LIN-45 C-terminal distal tail segment negatively regulates signaling in Caenorhabditis elegans

Raf protein kinases act as Ras-GTP sensing components of the ERK signal transduction pathway in animal cells, influencing cell proliferation, differentiation, and survival. In humans, somatic and germline mutations in the genes BRAF and RAF1 are associated with malignancies and developmental disorders. Recent studies shed light on the structure of activated Raf, a heterotetramer consisting of Raf and 14-3-3 dimers, and raised the possibility that a Raf C-terminal distal tail segment (DTS) regulates activation. We investigated the role of the DTS using the Caenorhabditis elegans Raf ortholog lin-45. Truncations removing the DTS strongly enhanced lin-45(S312A), a weak gain-of-function allele equivalent to RAF1 mutations found in patients with Noonan Syndrome. We genetically defined three elements of the LIN-45 DTS, which we termed the active site binding sequence (ASBS), the KTP motif, and the aromatic cluster. In the context of lin-45(S312A), the mutation of each of these elements enhanced activity. We used AlphaFold to predict DTS protein interactions for LIN-45, fly Raf, and human BRAF within the activated heterotetramer complex. We propose the following distinct functions for the LIN-45 DTS elements: (1) the ASBS binds the kinase active site as an inhibitor; (2) phosphorylation of the KTP motif modulates the DTS-kinase domain interaction; and (3) the aromatic cluster anchors the DTS in an inhibitory conformation. Human RASopathy-associated variants in BRAF affect residues of the DTS, consistent with these predictions. This work establishes that the Raf/LIN-45 DTS negatively regulates signaling in C. elegans and provides a model for its function in other Raf proteins.

Two hearts beat as one: the debate over RAS dimers continues

A recent report by Yun et al. describes the detection of RAS dimers using intact mass spectrometry and investigates the role that membrane lipids, nucleotide state, and binding partners have in their formation.

Comprehensive structure-function analysis reveals gain- and loss-of-function mechanisms impacting oncogenic KRAS activity

To dissect variant-function relationships in the KRAS oncoprotein, we performed deep mutational scanning (DMS) screens for both wild-type and KRASG12D mutant alleles. We defined the spectrum of oncogenic potential for nearly all possible KRAS variants, identifying several novel transforming alleles and elucidating a model to describe the frequency of KRAS mutations in human cancer as a function of transforming potential, mutational probability, and tissue-specific mutational signatures. Biochemical and structural analyses of variants identified in a KRASG12D second-site suppressor DMS screen revealed that attenuation of oncogenic KRAS can be mediated by protein instability and conformational rigidity, resulting in reduced binding affinity to effector proteins, such as RAF and PI3-kinases, or reduced SOS-mediated nucleotide exchange activity. These studies define the landscape of single amino acid alterations that modulate the function of KRAS, providing a resource for the clinical interpretation of KRAS variants and elucidating mechanisms of oncogenic KRAS inactivation for therapeutic exploitation.

A structurally informed human protein-protein interactome reveals proteome-wide perturbations caused by disease mutations

To assist the translation of genetic findings to disease pathobiology and therapeutics discovery, we present an ensemble deep learning framework, termed PIONEER (Protein-protein InteractiOn iNtErfacE pRediction), that predicts protein-binding partner-specific interfaces for all known protein interactions in humans and seven other common model organisms to generate comprehensive structurally informed protein interactomes. We demonstrate that PIONEER outperforms existing state-of-the-art methods and experimentally validate its predictions. We show that disease-associated mutations are enriched in PIONEER-predicted protein-protein interfaces and explore their impact on disease prognosis and drug responses. We identify 586 significant protein-protein interactions (PPIs) enriched with PIONEER-predicted interface somatic mutations (termed oncoPPIs) from analysis of approximately 11,000 whole exomes across 33 cancer types and show significant associations of oncoPPIs with patient survival and drug responses. PIONEER, implemented as both a web server platform and a software package, identifies functional consequences of disease-associated alleles and offers a deep learning tool for precision medicine at multiscale interactome network levels.

Generating Protein Structures for Pathway Discovery Using Deep Learning

Resolving the intricate details of biological phenomena at the molecular level is fundamentally limited by both length- and time scales that can be probed experimentally. Molecular dynamics (MD) simulations at various scales are powerful tools frequently employed to offer valuable biological insights beyond experimental resolution. However, while it is relatively simple to observe long-lived, stable configurations of, for example, proteins, at the required spatial resolution, simulating the more interesting rare transitions between such states often takes orders of magnitude longer than what is feasible even on the largest supercomputers available today. One common aspect of this challenge is pathway discovery, where the start and end states of a scientific phenomenon are known or can be approximated, but the mechanistic details in between are unknown. Here, we propose a representation-learning-based solution that uses interpolation and extrapolation in an abstract representation space to synthesize potential transition states, which are automatically validated using MD simulations. The new simulations of the synthesized transition states are subsequently incorporated into the representation learning, leading to an iterative framework for targeted path sampling. Our approach is demonstrated by recovering the transition of a RAS-RAF protein domain (CRD) from membrane-free to interacting with the membrane using coarse-grain MD simulations.

Effector Binding Sequentially Alters KRAS Dimerization on the Membrane: New Insights Into RAS-Mediated RAF Activation

RAS proteins are peripheral membrane GTPases that activate multiple downstream effectors for cell proliferation and differentiation. The formation of a signaling RAS-RAF complex at the plasma membrane is implicated in a quarter of all human cancers; however, the underlying mechanism remains unclear. In this work, nanodisc platforms and paramagnetic relaxation enhancement (PRE) analyses to determine the structure of a hetero-tetrameric complex comprising KRAS and the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of activated RAF1 are employed. The binding of the RBD or RBD-CRD differentially alters the dimerization modes of KRAS on both anionic and neutral membranes, validated by interface-specific mutagenesis. Notably, the RBD binding allosterically generated two distinct KRAS dimer interfaces in equilibrium, favored by KRAS free and in complex with the RBD-CRD, respectively. Additional interactions of the CRD with both KRAS protomers are mutually cooperative to stabilize a new dimer configuration of KRAS bound to the RBD-CRD. The RAF binding sequentially alters KRAS dimerization, providing new insights into RAF activation, including a configurational transition of the KRAS dimer to provide an interaction site for the CRD and release the autoinhibited RAF complex. These methods are applicable to many other signaling protein complexes on the membrane.

Mediating kinase activity in Ras-mutant cancer: potential for an individualised approach?

It is widely acknowledged that there is a considerable number of oncogenic mutations within the Ras superfamily of small GTPases which are the driving force behind a multitude of cancers. Ras proteins mediate a plethora of kinase pathways, including the MAPK, PI3K, and Ral pathways. Since Ras was considered undruggable until recently, pharmacological targeting of pathways downstream of Ras has been attempted to varying success, though drug resistance has often proven an issue. Nuances between kinase pathway activation in the presence of various Ras mutants are thought to contribute to the resistance, however, the reasoning behind activation of different pathways in different Ras mutational contexts is yet to be fully elucidated. Indeed, such disparities often depend on cancer type and disease progression. However, we are in a revolutionary age of Ras mutant targeted therapy, with direct-targeting KRAS-G12C inhibitors revolutionising the field and achieving FDA-approval in recent years. However, these are only beneficial in a subset of patients. Approximately 90% of Ras-mutant cancers are not KRAS-G12C mutant, and therefore raises the question as to whether other distinct amino acid substitutions within Ras may one day be targetable in a similar manner, and indeed whether better understanding of the downstream pathways these various mutants activate could further improve therapy. Here, we discuss the favouring of kinase pathways across an array of Ras-mutant oncogenic contexts and assess recent advances in pharmacological targeting of various Ras mutants. Ultimately, we will examine the utility of individualised pharmacological approaches to Ras-mediated cancer.

Targeting the RAS upstream and downstream signaling pathway for cancer treatment

Cancer often involves the overactivation of RAS/RAF/MEK/ERK (MAPK) and PI3K-Akt-mTOR pathways due to mutations in genes like RAS, RAF, PTEN, and PIK3CA. Various strategies are employed to address the overactivation of these pathways, among which targeted therapy emerges as a promising approach. Directly targeting specific proteins, leads to encouraging results in cancer treatment. For instance, RTK inhibitors such as imatinib and afatinib selectively target these receptors, hindering ligand binding and reducing signaling initiation. These inhibitors have shown potent efficacy against Non-Small Cell Lung Cancer. Other inhibitors, like lonafarnib targeting Farnesyltransferase and GGTI 2418 targeting geranylgeranyl Transferase, disrupt post-translational modifications of proteins. Additionally, inhibition of proteins like SOS, SH2 domain, and Ras demonstrate promising anti-tumor activity both in vivo and in vitro. Targeting downstream components with RAF inhibitors such as vemurafenib, dabrafenib, and sorafenib, along with MEK inhibitors like trametinib and binimetinib, has shown promising outcomes in treating cancers with BRAF-V600E mutations, including myeloma, colorectal, and thyroid cancers. Furthermore, inhibitors of PI3K (e.g., apitolisib, copanlisib), AKT (e.g., ipatasertib, perifosine), and mTOR (e.g., sirolimus, temsirolimus) exhibit promising efficacy against various cancers such as Invasive Breast Cancer, Lymphoma, Neoplasms, and Hematological malignancies. This review offers an overview of small molecule inhibitors targeting specific proteins within the RAS upstream and downstream signaling pathways in cancer.

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.

KRASG12C Inhibitors in Non-Small Cell Lung Cancer: A Review

Rat sarcoma virus (RAS) GTPase is one of the most important drivers of non-small cell lung cancer (NSCLC). RAS has three different isoforms (Harvey rat sarcoma viral oncogene homolog [HRAS], Kirsten rat sarcoma viral oncogene homolog [KRAS] and Neuroblastoma ras viral oncogene homolog [NRAS]), of which KRAS is most commonly mutated in NSCLC. The mutated KRAS protein was historically thought to be "undruggable" until the development of KRASG12C inhibitors. In this review, from the aspect of brain metastasis, we aim to provide an overview of the advances in therapies that target KRASG12C, the limitations of the current treatments, and future prospects in patients with KRAS p.G12C mutant NSCLC.

Capturing RAS oligomerization on a membrane

RAS GTPases associate with the biological membrane where they function as molecular switches to regulate cell growth. Recent studies indicate that RAS proteins oligomerize on membranes, and disrupting these assemblies represents an alternative therapeutic strategy. However, conflicting reports on RAS assemblies, ranging in size from dimers to nanoclusters, have brought to the fore key questions regarding the stoichiometry and parameters that influence oligomerization. Here, we probe three isoforms of RAS [Kirsten Rat Sarcoma viral oncogene (KRAS), Harvey Rat Sarcoma viral oncogene (HRAS), and Neuroblastoma oncogene (NRAS)] directly from membranes using mass spectrometry. We show that KRAS on membranes in the inactive state (GDP-bound) is monomeric but forms dimers in the active state (GTP-bound). We demonstrate that the small molecule BI2852 can induce dimerization of KRAS, whereas the binding of effector proteins disrupts dimerization. We also show that RAS dimerization is dependent on lipid composition and reveal that oligomerization of NRAS is regulated by palmitoylation. By monitoring the intrinsic GTPase activity of RAS, we capture the emergence of a dimer containing either mixed nucleotides or GDP on membranes. We find that the interaction of RAS with the catalytic domain of Son of Sevenless (SOScat) is influenced by membrane composition. We also capture the activation and monomer to dimer conversion of KRAS by SOScat. These results not only reveal the stoichiometry of RAS assemblies on membranes but also uncover the impact of critical factors on oligomerization, encompassing regulation by nucleotides, lipids, and palmitoylation.

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. We propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

Functional and structural insights into RAS effector proteins

RAS proteins are conserved guanosine triphosphate (GTP) hydrolases (GTPases) that act as molecular binary switches and play vital roles in numerous cellular processes. Upon GTP binding, RAS GTPases adopt an active conformation and interact with specific proteins termed RAS effectors that contain a conserved ubiquitin-like domain, thereby facilitating downstream signaling. Over 50 effector proteins have been identified in the human proteome, and many have been studied as potential mediators of RAS-dependent signaling pathways. Biochemical and structural analyses have provided mechanistic insights into these effectors, and studies using model organisms have complemented our understanding of their role in physiology and disease. Yet, many critical aspects regarding the dynamics and biological function of RAS-effector complexes remain to be elucidated. In this review, we discuss the mechanisms and functions of known RAS effector proteins, provide structural perspectives on RAS-effector interactions, evaluate their significance in RAS-mediated signaling, and explore their potential as therapeutic targets.

Signaling from RAS to RAF: The Molecules and Their Mechanisms

RAF family protein kinases are a key node in the RAS/RAF/MAP kinase pathway, the signaling cascade that controls cellular proliferation, differentiation, and survival in response to engagement of growth factor receptors on the cell surface. Over the past few years, structural and biochemical studies have provided new understanding of RAF autoregulation, RAF activation by RAS and the SHOC2 phosphatase complex, and RAF engagement with HSP90-CDC37 chaperone complexes. These studies have important implications for pharmacologic targeting of the pathway. They reveal RAF in distinct regulatory states and show that the functional RAF switch is an integrated complex of RAF with its substrate (MEK) and a 14-3-3 dimer. Here we review these advances, placing them in the context of decades of investigation of RAF regulation. We explore the insights they provide into aberrant activation of the pathway in cancer and RASopathies (developmental syndromes caused by germline mutations in components of the pathway).

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

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

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

The Raf/LIN-45 C-terminal distal tail segment negatively regulates signaling in Caenorhabditis elegans

Raf protein kinases act as Ras-GTP sensing components of the ERK signal transduction pathway in animal cells, influencing cell proliferation, differentiation, and survival. In humans, somatic and germline mutations in the genes BRAF and RAF1 are associated with malignancies and developmental disorders. Recent studies shed light on the structure of activated Raf, a heterotetramer consisting of Raf and 14-3-3 dimers, and raised the possibility that a Raf C-terminal distal tail segment (DTS) regulates activation. We investigated the role of the DTS using the Caenorhabditis elegans, which has a single Raf ortholog termed lin-45 . We discovered that truncations removing the DTS strongly enhanced lin-45(S312A) , a weak gain-of-function allele equivalent to RAF1 mutations found in patients with Noonan Syndrome. We generated mutations to test three elements of the LIN-45 DTS, which we termed the active site binding sequence (ASBS), the KTP motif, and the aromatic cluster. In the context of lin-45(S312A), mutation of either the ASBS, KTP motif, or aromatic cluster enhanced activity. We used AlphaFold to predict DTS protein interactions for LIN-45, fly Raf, and human BRAF, within the activated heterotetramer complex. We propose distinct functions for the LIN-45 DTS elements: i) the ASBS binds the kinase active site as an inhibitor, ii) phosphorylation of the KTP motif modulates DTS-kinase domain interaction, and iii) the aromatic cluster anchors the DTS in an inhibitory conformation. This work establishes that the Raf/LIN-45 DTS negatively regulates signaling in C. elegans and provides a model for its function in other Raf proteins.

Molecular determinants of Ras-mTORC2 signaling

Recent research has identified the mechanistic Target of Rapamycin Complex 2 (mTORC2) as a conserved direct effector of Ras proteins. While previous studies suggested the involvement of the Switch I (SWI) effector domain of Ras in binding mTORC2 components, the regulation of the Ras-mTORC2 pathway is not entirely understood. In Dictyostelium, mTORC2 is selectively activated by the Ras protein RasC, and the RasC-mTORC2 pathway then mediates chemotaxis to cAMP and cellular aggregation by regulating the actin cytoskeleton and promoting cAMP signal relay. Here, we investigated the role of specific residues in RasC's SWI, C-terminal allosteric domain, and hypervariable region (HVR) related to mTORC2 activation. Interestingly, our results suggest that RasC SWI residue A31, which was previously implicated in RasC-mediated aggregation, regulates RasC's specific activation by the Aimless RasGEF. On the other hand, our investigation identified a crucial role for RasC SWI residue T36, with secondary contributions from E38 and allosteric domain residues. Finally, we found that conserved basic residues and the adjacent prenylation site in the HVR, which are crucial for RasC's membrane localization, are essential for RasC-mTORC2 pathway activation by allowing for both RasC's own cAMP-induced activation and its subsequent activation of mTORC2. Therefore, our findings revealed new determinants of RasC-mTORC2 pathway specificity in Dictyostelium, contributing to a deeper understanding of Ras signaling regulation in eukaryotic cells.

Clonal competition assays identify fitness signatures in cancer progression and resistance in multiple myeloma

Multiple myeloma (MM) is a genetically heterogeneous disease and the management of relapses is one of the biggest clinical challenges. TP53 alterations are established high-risk markers and are included in the current disease staging criteria. KRAS is the most frequently mutated gene affecting around 20% of MM patients. Applying Clonal Competition Assays (CCA) by co-culturing color-labeled genetically modified cell models, we recently showed that mono- and biallelic alterations in TP53 transmit a fitness advantage to the cells. Here, we report a similar dynamic for two mutations in KRAS (G12A and A146T), providing a biological rationale for the high frequency of KRAS and TP53 alterations at MM relapse. Resistance mutations, on the other hand, did not endow MM cells with a general fitness advantage but rather presented a disadvantage compared to the wild-type. CUL4B KO and IKZF1 A152T transmit resistance against immunomodulatory agents, PSMB5 A20T to proteasome inhibition. However, MM cells harboring such lesions only outcompete the culture in the presence of the respective drug. To better prevent the selection of clones with the potential of inducing relapse, these results argue in favor of treatment-free breaks or a switch of the drug class given as maintenance therapy. In summary, the fitness benefit of TP53 and KRAS mutations was not treatment-related, unlike patient-derived drug resistance alterations that may only induce an advantage under treatment. CCAs are suitable models for the study of clonal evolution and competitive (dis)advantages conveyed by a specific genetic lesion of interest, and their dependence on external factors such as the treatment.

Regulation of RAF family kinases: new insights from recent structural and biochemical studies

The RAF kinases are required for signal transduction through the RAS-RAF-MEK-ERK pathway, and their activity is frequently up-regulated in human cancer and the RASopathy developmental syndromes. Due to their complex activation process, developing drugs that effectively target RAF function has been a challenging endeavor, highlighting the need for a more detailed understanding of RAF regulation. This review will focus on recent structural and biochemical studies that have provided 'snapshots' into the RAF regulatory cycle, revealing structures of the autoinhibited BRAF monomer, active BRAF and CRAF homodimers, as well as HSP90/CDC37 chaperone complexes containing CRAF or BRAFV600E. In addition, we will describe the insights obtained regarding how BRAF transitions between its regulatory states and examine the roles that various BRAF domains and 14-3-3 dimers play in both maintaining BRAF as an autoinhibited monomer and in facilitating its transition to an active dimer. We will also address the function of the HSP90/CDC37 chaperone complex in stabilizing the protein levels of CRAF and certain oncogenic BRAF mutants, and in serving as a platform for RAF dephosphorylation mediated by the PP5 protein phosphatase. Finally, we will discuss the regulatory differences observed between BRAF and CRAF and how these differences impact the function of BRAF and CRAF as drivers of human disease.

Deciphering the pathogenic role of rare RAF1 heterozygous missense mutation in the late-presenting DDH

Background

Developmental Dysplasia of the Hip (DDH) is a skeletal disorder where late-presenting forms often escape early diagnosis, leading to limb and pain in adults. The genetic basis of DDH is not fully understood despite known genetic predispositions.

Methods

We employed Whole Genome Sequencing (WGS) to explore the genetic factors in late-presenting DDH in two unrelated families, supported by phenotypic analyses and in vitro validation.

Results

In both cases, a novel de novo heterozygous missense mutation in RAF1 (c.193A>G [p.Lys65Glu]) was identified. This mutation impacted RAF1 protein structure and function, altering downstream signaling in the Ras/ERK pathway, as demonstrated by bioinformatics, molecular dynamics simulations, and in vitro validations.

Conclusion

This study contributes to our understanding of the genetic factors involved in DDH by identifying a novel mutation in RAF1. The identification of the RAF1 mutation suggests a possible involvement of the Ras/ERK pathway in the pathogenesis of late-presenting DDH, indicating its potential role in skeletal development.

RAS and SHOC2 Roles in RAF Activation and Therapeutic Considerations

Mutations in RAS proteins play a pivotal role in the development of human cancers, driving persistent RAF activation and deregulating the Mitogen-Activated Protein Kinase (MAPK) signaling pathway. While progress has been made in targeting specific oncogenic RAS proteins, effective drug-based therapies for the majority of RAS mutations remain limited. Recent investigations on RAS-RAF complexes and the SHOC2-MRAS-PP1C holoenzyme complex have provided crucial insights into the structural and functional aspects of RAF activation within the MAPK signaling pathway. Moreover, these studies have also unveiled new blueprints for developing inhibitors allowing us to think beyond the current RAS and MEK inhibitors. In this review, we explore the roles of RAS and SHOC2 in activating RAF and discuss potential therapeutic strategies to target these proteins. A comprehensive understanding of the molecular interactions involved in RAF activation and their therapeutic implications holds the potential to drive innovative approaches in combating RAS/RAF-driven cancers.

Reconstitution and characterization of BRAF in complex with 14-3-3 and KRAS4B on nanodiscs

RAF kinases are key components of the RAS-MAPK signaling pathway, which drives cell growth and is frequently overactivated in cancer. Upstream signaling activates the small GTPase RAS, which recruits RAF to the cell membrane, driving a transition of the latter from an auto-inhibited monomeric conformation to an active dimer. Despite recent progress, mechanistic details underlying RAF activation remain unclear, particularly the role of RAS and the membrane in mediating this conformational rearrangement of RAF together with 14-3-3 to permit RAF kinase domain dimerization. Here, we reconstituted an active complex of dimeric BRAF, a 14-3-3 dimer and two KRAS4B on a nanodisc bilayer and verified that its assembly is GTP-dependent. Biolayer interferometry (BLI) was used to compare the binding affinities of monomeric versus dimeric full-length BRAF:14-3-3 complexes for KRAS4B-conjugated nanodiscs (RAS-ND) and to investigate the effects of membrane lipid composition and spatial density of KRAS4B on binding. 1,2-Dioleoyl-sn-glycero-3-phospho-L-serine (DOPS) and higher KRAS4B density enhanced the interaction of BRAF:14-3-3 with RAS-ND to different degrees depending on BRAF oligomeric state. We utilized our reconstituted system to dissect the effects of KRAS4B and the membrane on the kinase activity of monomeric and dimeric BRAF:14-3-3 complexes, finding that KRAS4B or nanodiscs alone were insufficient to stimulate activity, whereas RAS-ND increased activity of both states of BRAF. The reconstituted assembly of full-length BRAF with 14-3-3 and KRAS on a cell-free, defined lipid bilayer offers a more holistic biophysical perspective to probe regulation of this multimeric signaling complex at the membrane surface.

RAF1 gene fusions are recurrent driver events in infantile fibrosarcoma-like mesenchymal tumors

Infantile fibrosarcomas (IFS) and congenital mesoblastic nephroma (CMN) are rare myofibroblastic tumors of infancy and early childhood commonly harboring the ETV6::NTRK3 gene fusion. IFS/CMN are considered as tumors with an 'intermediate prognosis' as they are locally aggressive, but rarely metastasize, and generally have a favorable outcome. A fraction of IFS/CMN-related neoplasms are negative for the ETV6::NTRK3 gene rearrangement and are characterized by other chimeric proteins promoting MAPK signaling upregulation. In a large proportion of these tumors, which are classified as IFS-like mesenchymal neoplasms, the contributing molecular events remain to be identified. Here, we report three distinct rearrangements involving RAF1 among eight ETV6::NTRK3 gene fusion-negative tumors with an original histological diagnosis of IFS/CMN. The three fusion proteins retain the entire catalytic domain of the kinase. Two chimeric products, GOLGA4::RAF1 and LRRFIP2::RAF1, had previously been reported as driver events in different cancers, whereas the third, CLIP1::RAF1, represents a novel fusion protein. We demonstrate that CLIP1::RAF1 acts as a bona fide oncoprotein promoting cell proliferation and migration through constitutive upregulation of MAPK signaling. We show that the CLIP1::RAF1 hyperactive behavior does not require RAS activation and is mediated by constitutive 14-3-3 protein-independent dimerization of the chimeric protein. As previously reported for the ETV6::NTRK3 fusion protein, CLIP1::RAF1 similarly upregulates PI3K-AKT signaling. Our findings document that RAF1 gene rearrangements represent a recurrent event in ETV6::NTRK3-negative IFS/CMN and provide a rationale for the use of inhibitors directed to suppress MAPK and PI3K-AKT signaling in these cancers. © 2024 The Pathological Society of Great Britain and Ireland.

SIRT4 as a novel interactor and candidate suppressor of C-RAF kinase in MAPK signaling

Cellular responses leading to development, proliferation, and differentiation depend on RAF/MEK/ERK signaling, which integrates and amplifies signals from various stimuli for downstream cellular responses. C-RAF activation has been reported in many types of tumor cell proliferation and developmental disorders, necessitating the discovery of potential C-RAF protein regulators. Here, we identify a novel and specific protein interaction between C-RAF among the RAF kinase paralogs, and SIRT4 among the mitochondrial sirtuin family members SIRT3, SIRT4, and SIRT5. Structurally, C-RAF binds to SIRT4 through the N-terminal cysteine-rich domain, whereas SIRT4 predominantly requires the C-terminus for full interaction with C-RAF. Interestingly, SIRT4 specifically interacts with C-RAF in a pre-signaling inactive (serine 259-phosphorylated) state. Consistent with this finding, the expression of SIRT4 in HEK293 cells results in an up-regulation of pS259-C-RAF levels and a concomitant reduction in MAPK signaling as evidenced by strongly decreased phospho-ERK signals. Thus, we propose an additional extra-mitochondrial function of SIRT4 as a cytosolic tumor suppressor of C-RAF-MAPK signaling, besides its metabolic tumor suppressor role of glutamate dehydrogenase and glutamate levels in mitochondria.

Carbene Footprinting Directs Design of Genetically Encoded Proximity-Reactive Protein Binders

Genetically encoding proximal-reactive unnatural amino acids (PrUaas), such as fluorosulfate-l-tyrosine (FSY), into natural proteins of interest (POI) confer the POI with the ability to covalently bind to its interacting proteins (IPs). The PrUaa-incorporated POIs hold promise for blocking undesirable POI-IP interactions. Selecting appropriate PrUaa anchor sites is crucial, but it remains challenging with the current methodology, which heavily relies on crystallography to identify the proximal residues between the POIs and the IPs for the PrUaa anchorage. To address the challenge, here, we propose a footprinting-directed genetically encoded covalent binder (footprinting-GECB) approach. This approach employs carbene footprinting, a structural mass spectrometry (MS) technique that quantifies the extent of labeling of the POI following the addition of its IP, and thus identifies the responsive residues. By genetically encoding PrUaa into these responsive sites, POI variants with covalent bonding ability to its IP can be produced without the need for crystallography. Using the POI-IP model, KRAS/RAF1, we showed that engineering FSY at the footprint-assigned KRAS residue resulted in a KRAS variant that can bind irreversibly to RAF1. Additionally, we inserted FSY at the responsive residue in RAF1 upon footprinting the oncogenic KRASG12D/RAF1, which lacks crystal structure, and generated a covalent binder to KRASG12D. Together, we demonstrated that by adopting carbene footprinting to direct PrUaa anchorage, we can greatly expand the opportunities for designing covalent protein binders for PPIs without relying on crystallography. This holds promise for creating effective PPI inhibitors and supports both fundamental research and biotherapeutics development.

The RAF cysteine-rich domain: Structure, function, and role in disease

RAF kinases, consisting of ARAF, BRAF and CRAF, are direct effectors of RAS GTPases and critical for signal transduction through the RAS-MAPK pathway. Driver mutations in BRAF are commonplace in human cancer, while germline mutations in BRAF and CRAF cause RASopathy development syndromes. However, there remains a lack of effective drugs that target RAF function, which is partially due to the complexity of the RAF activation cycle. Therefore, greater understanding of RAF regulation is required to identify new approaches that target its function in disease. A key piece of this puzzle is the RAF zinc finger, often referred to as the cysteine-rich domain (CRD). The CRD is a lipid and protein binding domain which plays complex and opposing roles in the RAF activation cycle. Firstly, it supports the RAS-RAF interaction during RAF activation by binding to phosphatidylserine (PS) in the plasma membrane and by making direct RAS contacts. Conversely, under quiescent conditions the CRD also plays a critical role in maintaining RAF in a closed, autoinhibited state. However, the interplay between these activities and their relative importance for RAF activation were not well understood. Recent structural and biochemical studies have contributed greatly to our understanding of these roles and identified functional differences between BRAF CRD and that of CRAF. This chapter provides an in-depth review of the CRDs roles in RAF regulation and how they may inform novel approaches to target RAF function.

PreMode predicts mode-of-action of missense variants by deep graph representation learning of protein sequence and structural context

Accurate prediction of the functional impact of missense variants is important for disease gene discovery, clinical genetic diagnostics, therapeutic strategies, and protein engineering. Previous efforts have focused on predicting a binary pathogenicity classification, but the functional impact of missense variants is multi-dimensional. Pathogenic missense variants in the same gene may act through different modes of action (i.e., gain/loss-of-function) by affecting different aspects of protein function. They may result in distinct clinical conditions that require different treatments. We developed a new method, PreMode, to perform gene-specific mode-of-action predictions. PreMode models effects of coding sequence variants using SE(3)-equivariant graph neural networks on protein sequences and structures. Using the largest-to-date set of missense variants with known modes of action, we showed that PreMode reached state-of-the-art performance in multiple types of mode-of-action predictions by efficient transfer-learning. Additionally, PreMode's prediction of G/LoF variants in a kinase is consistent with inactive-active conformation transition energy changes. Finally, we show that PreMode enables efficient study design of deep mutational scans and optimization in protein engineering.

RAF and MEK Inhibitors in Non-Small Cell Lung Cancer

Lung cancer, despite recent advancements in survival rates, represents a significant global health burden. Non-small cell lung cancer (NSCLC), the most prevalent type, is driven largely by activating mutations in Kirsten rat sarcoma viral oncogene homologue (KRAS) and receptor tyrosine kinases (RTKs), and less in v-RAF murine sarcoma viral oncogene homolog B (BRAF) and mitogen-activated protein-kinase kinase (MEK), all key components of the RTK-RAS-mitogen-activated protein kinase (MAPK) pathway. Learning from melanoma, the identification of BRAFV600E substitution in NSCLC provided the rationale for the investigation of RAF and MEK inhibition as a therapeutic strategy. The regulatory approval of two RAF-MEK inhibitor combinations, dabrafenib-trametinib, in 2017, and encorafenib-binimetinib, in 2023, signifies a breakthrough for the management of BRAFV600E-mutant NSCLC patients. However, the almost universal emergence of acquired resistance limits their clinical benefit. New RAF and MEK inhibitors, with distinct biochemical characteristics, are in preclinical and clinical development. In this review, we aim to provide valuable insights into the current state of RAF and MEK inhibition in the management of NSCLC, fostering a deeper understanding of the potential impact on patient outcomes.

The cell biology of HIV-1 latency and rebound

Transcriptionally latent forms of replication-competent proviruses, present primarily in a small subset of memory CD4+ T cells, pose the primary barrier to a cure for HIV-1 infection because they are the source of the viral rebound that almost inevitably follows the interruption of antiretroviral therapy. Over the last 30 years, many of the factors essential for initiating HIV-1 transcription have been identified in studies performed using transformed cell lines, such as the Jurkat T-cell model. However, as highlighted in this review, several poorly understood mechanisms still need to be elucidated, including the molecular basis for promoter-proximal pausing of the transcribing complex and the detailed mechanism of the delivery of P-TEFb from 7SK snRNP. Furthermore, the central paradox of HIV-1 transcription remains unsolved: how are the initial rounds of transcription achieved in the absence of Tat? A critical limitation of the transformed cell models is that they do not recapitulate the transitions between active effector cells and quiescent memory T cells. Therefore, investigation of the molecular mechanisms of HIV-1 latency reversal and LRA efficacy in a proper physiological context requires the utilization of primary cell models. Recent mechanistic studies of HIV-1 transcription using latently infected cells recovered from donors and ex vivo cellular models of viral latency have demonstrated that the primary blocks to HIV-1 transcription in memory CD4+ T cells are restrictive epigenetic features at the proviral promoter, the cytoplasmic sequestration of key transcription initiation factors such as NFAT and NF-κB, and the vanishingly low expression of the cellular transcription elongation factor P-TEFb. One of the foremost schemes to eliminate the residual reservoir is to deliberately reactivate latent HIV-1 proviruses to enable clearance of persisting latently infected cells-the "Shock and Kill" strategy. For "Shock and Kill" to become efficient, effective, non-toxic latency-reversing agents (LRAs) must be discovered. Since multiple restrictions limit viral reactivation in primary cells, understanding the T-cell signaling mechanisms that are essential for stimulating P-TEFb biogenesis, initiation factor activation, and reversing the proviral epigenetic restrictions have become a prerequisite for the development of more effective LRAs.

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.

Conditional Cooperativity in RAS Assembly Pathways on Nanodiscs and Altered GTPase Cycling

Self-assemblies (i.e., nanoclusters) of the RAS GTPase on the membrane act as scaffolds that activate downstream RAF kinases and drive MAPK signaling for cell proliferation and tumorigenesis. However, the mechanistic details of nanoclustering remain largely unknown. Here, size-tunable nanodisc platforms and paramagnetic relaxation enhancement (PRE) analyses revealed the structural basis of the cooperative assembly processes of fully processed KRAS, mutated in a quarter of human cancers. The cooperativity is modulated by the mutation and nucleotide states of KRAS and the lipid composition of the membrane. Notably, the oncogenic mutants assemble in nonsequential pathways with two mutually cooperative 'α/α' and 'α/β' interfaces, while α/α dimerization of wild-type KRAS promotes the secondary α/β interaction sequentially. Mutation-based interface engineering was used to selectively trap the oligomeric intermediates of KRAS and probe their favorable interface interactions. Transiently exposed interfaces were available for the assembly. Real-time NMR demonstrated that higher-order oligomers retain higher numbers of active GTP-bound protomers in KRAS GTPase cycling. These data provide a deeper understanding of the nanocluster-enhanced signaling in response to the environment. Furthermore, our methodology is applicable to assemblies of many other membrane GTPases and lipid nanoparticle-based formulations of stable protein oligomers with enhanced cooperativity.

Cryo-EM Structures of CRAF2/14-3-32 and CRAF2/14-3-32/MEK12 Complexes

RAF protein kinases are essential effectors in the MAPK pathway and are important cancer drug targets. Structural understanding of RAF activation is so far based on cryo-electron microscopy (cryo-EM) and X-ray structures of BRAF in different conformational states as inactive or active complexes with KRAS, 14-3-3 and MEK1. In this study, we have solved the first cryo-EM structures of CRAF2/14-3-32 at 3.4 Å resolution and CRAF2/14-3-32/MEK12 at 4.2 Å resolution using CRAF kinase domain expressed as constitutively active Y340D/Y341D mutant in insect cells. The overall architecture of our CRAF2/14-3-32 and CRAF2/14-3-32/MEK12 cryo-EM structures is highly similar to corresponding BRAF structures in complex with 14-3-3 or 14-3-3/MEK1 and represent the activated dimeric RAF conformation. Our CRAF cryo-EM structures provide additional insights into structural understanding of the activated CRAF2/14-3-32/MEK12 complex.

BRAF - a tumour-agnostic drug target with lineage-specific dependencies

In June 2022, the FDA granted Accelerated Approval to the BRAF inhibitor dabrafenib in combination with the MEK inhibitor trametinib for the treatment of adult and paediatric patients (≥6 years of age) with unresectable or metastatic BRAFV600E-mutant solid tumours, except for BRAFV600E-mutant colorectal cancers. The histology-agnostic approval of dabrafenib plus trametinib marks the culmination of two decades of research into the landscape of BRAF mutations in human cancers, the biochemical mechanisms underlying BRAF-mediated tumorigenesis, and the clinical development of selective RAF and MEK inhibitors. Although the majority of patients with BRAFV600E-mutant tumours derive clinical benefit from BRAF inhibitor-based combinations, resistance to treatment develops in most. In this Review, we describe the biochemical basis for oncogenic BRAF-induced activation of MAPK signalling and pan-cancer and lineage-specific mechanisms of intrinsic, adaptive and acquired resistance to BRAF inhibitors. We also discuss novel RAF inhibitors and drug combinations designed to delay the emergence of treatment resistance and/or expand the population of patients with BRAF-mutant cancers who benefit from molecularly targeted therapies.

Live-cell target engagement of allosteric MEKi on MEK-RAF/KSR-14-3-3 complexes

The RAS-mitogen-activated protein kinase (MAPK) pathway includes KSR, RAF, MEK and the phospho-regulatory sensor 14-3-3. Specific assemblies among these components drive various diseases and likely dictate efficacy for numerous targeted therapies, including allosteric MEK inhibitors (MEKi). However, directly measuring drug interactions on physiological RAS-MAPK complexes in live cells has been inherently challenging to query and therefore remains poorly understood. Here we present a series of NanoBRET-based assays to quantify direct target engagement of MEKi on MEK1 and higher-order MEK1-bound complexes with ARAF, BRAF, CRAF, KSR1 and KSR2 in the presence and absence of 14-3-3 in living cells. We find distinct MEKi preferences among these complexes that can be compiled to generate inhibitor binding profiles. Further, these assays can report on the influence of the pathogenic BRAF-V600E mutant on MEKi binding. Taken together, these approaches can be used as a platform to screen for compounds intended to target specific complexes in the RAS-MAPK cascade.

Membrane lipids drive formation of KRAS4b-RAF1 RBDCRD nanoclusters on the membrane

The oncogene RAS, extensively studied for decades, presents persistent gaps in understanding, hindering the development of effective therapeutic strategies due to a lack of precise details on how RAS initiates MAPK signaling with RAF effector proteins at the plasma membrane. Recent advances in X-ray crystallography, cryo-EM, and super-resolution fluorescence microscopy offer structural and spatial insights, yet the molecular mechanisms involving protein-protein and protein-lipid interactions in RAS-mediated signaling require further characterization. This study utilizes single-molecule experimental techniques, nuclear magnetic resonance spectroscopy, and the computational Machine-Learned Modeling Infrastructure (MuMMI) to examine KRAS4b and RAF1 on a biologically relevant lipid bilayer. MuMMI captures long-timescale events while preserving detailed atomic descriptions, providing testable models for experimental validation. Both in vitro and computational studies reveal that RBDCRD binding alters KRAS lateral diffusion on the lipid bilayer, increasing cluster size and decreasing diffusion. RAS and membrane binding cause hydrophobic residues in the CRD region to penetrate the bilayer, stabilizing complexes through β-strand elongation. These cooperative interactions among lipids, KRAS4b, and RAF1 are proposed as essential for forming nanoclusters, potentially a critical step in MAP kinase signal activation.