A structural model of a Ras-Raf signalosome

Structures, dynamics, complexes, and functions: From classic computation to artificial intelligence

Computational approaches can provide highly detailed insight into the molecular recognition processes that underlie drug binding, the assembly of protein complexes, and the regulation of biological functional processes. Classical simulation methods can bridge a wide range of length- and time-scales typically involved in such processes. Lately, automated learning and artificial intelligence methods have shown the potential to expand the reach of physics-based approaches, ushering in the possibility to model and even design complex protein architectures. The synergy between atomistic simulations and AI methods is an emerging frontier with a huge potential for advances in structural biology. Herein, we explore various examples and frameworks for these approaches, providing select instances and applications that illustrate their impact on fundamental biomolecular problems.

RAS G-domains allosterically contribute to the recognition of lipid headgroups and acyl chains

Mutant RAS are major contributors to cancer and signal primarily from nanoclusters on the plasma membrane (PM). Their C-terminal membrane anchors are main features of membrane association. However, the same RAS isoform bound to different guanine nucleotides spatially segregate. Different RAS nanoclusters all enrich a phospholipid, phosphatidylserine (PS). These findings suggest more complex membrane interactions. Our electron microscopy-spatial analysis shows that wild-types, G12V mutants, and membrane anchors of isoforms HRAS, KRAS4A, and KRAS4B prefer distinct PS species. Mechanistically, reorientation of KRAS4B G-domain exposes distinct residues, such as Arg 135 in orientation state 1 (OS1) and Arg 73/Arg 102 in OS2, to the PM and differentially facilitates the recognition of PS acyl chains. Allele-specific oncogenic mutations of KRAS4B also shift G-domain reorientation equilibrium. Indeed, KRAS4BG12V, KRAS4BG12D, KRAS4BG12C, KRAS4BG13D, and KRAS4BQ61H associate with PM lipids with headgroup and acyl chain specificities. Distribution of these KRAS4B oncogenic mutants favors different nanoscale membrane topography. Thus, RAS G-domains allosterically facilitate membrane lateral distribution.

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.

Development of a Data-Driven Integrative Model of a Bacterial Chromosome

The chromosome of archetypal bacteria E. coli is known for a complex topology with a 4.6 × 106 base pairs (bp) long sequence of nucleotides packed within a micrometer-sized cellular confinement. The inherent organization underlying this chromosome eludes general consensus due to the lack of a high-resolution picture of its conformation. Here we present our development of an integrative model of E. coli at a 500 bp resolution (https://github.com/JMLab-tifrh/ecoli_finer), which optimally combines a set of multiresolution genome-wide experimentally measured data within a framework of polymer based architecture. In particular the model is informed with an intragenome contact probability map at 5000 bp resolution derived via the Hi-C experiment and RNA-sequencing data at 500 bp resolution. Via dynamical simulations, this data-driven polymer based model generates an appropriate conformational ensemble commensurate with chromosome architectures that E. coli adopts. As a key hallmark of the E. coli chromosome the model spontaneously self-organizes into a set of nonoverlapping macrodomains and suitably locates plectonemic loops near the cell membrane. As novel extensions, it predicts a contact probability map simulated at a higher resolution than precedent experiments and can demonstrate segregation of chromosomes in a partially replicating cell. Finally, the modular nature of the model helps us devise control simulations to quantify the individual role of key features in hierarchical organization of the bacterial chromosome.

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications

Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.

Crystal Packing Reveals a Potential Autoinhibited KRAS Dimer Interface and a Strategy for Small-Molecule Inhibition of RAS Signaling

KRAS GTPases harbor oncogenic mutations in more than 25% of human tumors. KRAS is considered to be largely undruggable due to the lack of a suitable small-molecule binding site. Here, we report a unique crystal structure of His-tagged KRASG12D that reveals a remarkable conformational change. The Switch I loop of one His-KRASG12D structure extends into the Switch I/II pocket of another His-KRASG12D in an adjacent unit cell to create an elaborate interface that is reminiscent of high-affinity protein-protein complexes. We explore the contributions of amino acids at this interface using alanine-scanning studies with alchemical free energy perturbation calculations based on explicit-solvent molecular dynamics simulations. Several interface amino acids were found to be hot spots as they contributed more than 1.5 kcal/mol to the protein-protein interaction. Computational analysis of the complex revealed the presence of two large binding pockets that possess physicochemical features typically found in pockets considered druggable. Small-molecule binding to these pockets may stabilize this autoinhibited structure of KRAS if it exists in cells to provide a new strategy to inhibit RAS signaling.

Emodin derivative E35 and its combination with autophagy inhibitor against acute myeloid leukemia cells in vitro and in vivo

Acute myeloid leukemia (AML) is an aggressive hematopoietic malignancy with poor prognosis and high recurrence rate. The discovery of more effective therapeutic strategies for AML plays a crucial role. The present work showed that E35, a novel derivative of emodin, significantly inhibited cell proliferation and induced autophagy and apoptosis in AML cells. Treatment with E35 markedly induced Beclin-1, LC3-II, cleaved Caspase-9 and PARP, and suppressed mitogen-activated protein kinase (MAPK) pathway. E35 exposure evoked autophagic activity prior to apoptosis induction, and autophagy inhibition by 3-methyladenine (3-MA) dramatically increased E35-induced apoptosis in both AML cell lines and patient-derived AML cells. Nevertheless, study on AML xenograft model showed that the combination E35 with 3-MA exhibited much more inhibitory effects on leukemia cell growth in vivo. No obvious adverse reactions occurred in the xenograft animals administered E35 alone or its cotreatment with 3-MA. These findings suggest that E35 could exert anti-leukemia effects, and that the combination of E35 and autophagy inhibitor might prove a more highly efficient strategy for AML treatment.

GLUT3 promotes macrophage signaling and function via RAS-mediated endocytosis in atopic dermatitis and wound healing

The facilitative GLUT1 and GLUT3 hexose transporters are expressed abundantly in macrophages, but whether they have distinct functions remains unclear. We confirmed that GLUT1 expression increased after M1 polarization stimuli and found that GLUT3 expression increased after M2 stimulation in macrophages. Conditional deletion of Glut3 (LysM-Cre Glut3fl/fl) impaired M2 polarization of bone marrow-derived macrophages. Alternatively activated macrophages from the skin of patients with atopic dermatitis showed increased GLUT3 expression, and a calcipotriol-induced model of atopic dermatitis was rescued in LysM-Cre Glut3fl/fl mice. M2-like macrophages expressed GLUT3 in human wound tissues as assessed by transcriptomics and costaining, and GLUT3 expression was significantly decreased in nonhealing, compared with healing, diabetic foot ulcers. In an excisional wound healing model, LysM-Cre Glut3fl/fl mice showed significantly impaired M2 macrophage polarization and delayed wound healing. GLUT3 promoted IL-4/STAT6 signaling, independently of its glucose transport activity. Unlike plasma membrane-localized GLUT1, GLUT3 was localized primarily to endosomes and was required for the efficient endocytosis of IL-4Rα subunits. GLUT3 interacted directly with GTP-bound RAS in vitro and in vivo through its intracytoplasmic loop domain, and this interaction was required for efficient STAT6 activation and M2 polarization. PAK activation and macropinocytosis were also impaired without GLUT3, suggesting broader roles for GLUT3 in the regulation of endocytosis. Thus, GLUT3 is required for efficient alternative macrophage polarization and function, through a glucose transport-independent, RAS-mediated role in the regulation of endocytosis and IL-4/STAT6 activation.

Mechanism-Based Redesign of GAP to Activate Oncogenic Ras

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

Structural and functional analyses of a germline KRAS T50I mutation provide insights into Raf activation

A T50I substitution in the K-Ras interswitch domain causes Noonan syndrome and emerged as a third-site mutation that restored the in vivo transforming activity and constitutive MAPK pathway activation by an attenuated KrasG12D,E37G oncogene in a mouse leukemia model. Biochemical and crystallographic data suggested that K-RasT50I increases MAPK signal output through a non-GTPase mechanism, potentially by promoting asymmetric Ras:Ras interactions between T50 and E162. We generated a "switchable" system in which K-Ras mutant proteins expressed at physiologic levels supplant the fms like tyrosine kinase 3 (FLT3) dependency of MOLM-13 leukemia cells lacking endogenous KRAS and used this system to interrogate single or compound G12D, T50I, D154Q, and E162L mutations. These studies support a key role for the asymmetric lateral assembly of K-Ras in a plasma membrane-distal orientation that promotes the formation of active Ras:Raf complexes in a membrane-proximal conformation. Disease-causing mutations such as T50I are a valuable starting point for illuminating normal Ras function, elucidating mechanisms of disease, and identifying potential therapeutic opportunities for Rasopathy disorders and cancer.

The E3/E4 ubiquitin ligase UFD-2 suppresses normal and oncogenic signaling mediated by a Raf ortholog in Caenorhabditis elegans

Signaling by the kinase cascade composed of Raf, MEK, and ERK is critical for animal development and is often inappropriately activated in human malignancies. We sought to identify factors that control signaling mediated by the Caenorhabditis elegans Raf ortholog LIN-45. A genetic screen showed that the degradation of LIN-45 required the E3/E4 ubiquitin ligase UFD-2. Both UFD-2 and its partner, the ATP-dependent segregase CDC-48, were required for the developmental regulation of LIN-45 protein abundance. We showed that UFD-2 acted in the same pathway as the E3 ubiquitin ligase SCFSEL-10 to decrease LIN-45 abundance in cells in which Raf-MEK-ERK signaling was most highly active. UFD-2 also reduced the protein abundance of activated LIN-45 carrying a mutation equivalent to the cancer-associated BRAF(V600E) variant. Our structure-function studies showed that the disruption of LIN-45 domains that mediate protein-protein interactions, including the conserved cysteine-rich domain and 14-3-3 binding motifs, were required for UFD-2-independent degradation of LIN-45. We propose a model in which UFD-2 and CDC-48 act downstream of SCFSEL-10 to remove LIN-45 from its protein interaction partners and facilitate proteasomal targeting and degradation. These findings imply that UFD-2 and CDC-48 may be important for Raf degradation during normal and oncogenic Ras and MAPK signaling in mammalian cells.

The interrelation of galectins and autophagy

Autophagy is a vital physiological process that maintains intracellular homeostasis by removing damaged organelles and senescent or misfolded molecules. However, excessive autophagy results in cell death and apoptosis, which will lead to a variety of diseases. Galectins are a type of animal lectin that binds to β-galactosides and can bind to the cell surface or extracellular matrix glycans, affecting a variety of immune processes in vivo and being linked to the development of many diseases. In many cases, galectins and autophagy both play important regulatory roles in the cellular life course, yet our understanding of the relationship between them is still incomplete. Galectins and autophagy may share common etiological cofactors for some diseases. Hence, we summarize the relationship between galectins and autophagy, aiming to draw attention to the existence of multiple associations between galectins and autophagy in a variety of physiological and pathological processes, which provide new ideas for etiological diagnosis, drug development, and therapeutic targets for related diseases.

Chimera and Tandem-Repeat Type Galectins: The New Targets for Cancer Immunotherapy

In humans, a total of 12 galectins have been identified. Their intracellular and extracellular biological functions are explored and discussed in this review. These galectins play important roles in controlling immune responses within the tumour microenvironment (TME) and the infiltration of immune cells, including different subsets of T cells, macrophages, and neutrophils, to fight against cancer cells. However, these infiltrating cells also have repair roles and are hijacked by cancer cells for pro-tumorigenic activities. Upon a better understanding of the immunomodulating functions of galectin-3 and -9, their inhibitors, namely, GB1211 and LYT-200, have been selected as candidates for clinical trials. The use of these galectin inhibitors as combined treatments with current immune checkpoint inhibitors (ICIs) is also undergoing clinical trial investigations. Through their network of binding partners, inhibition of galectin have broad downstream effects acting on CD8⁺ cytotoxic T cells, regulatory T cells (Tregs), Natural Killer (NK) cells, and macrophages as well as playing pro-inflammatory roles, inhibiting T-cell exhaustion to support the fight against cancer cells. Other galectin members are also included in this review to provide insight into potential candidates for future treatment(s). The pitfalls and limitations of using galectins and their inhibitors are also discussed to cognise their clinical application.

The Self-Association of the KRAS4b Protein is Altered by Lipid-Bilayer Composition and Electrostatics

KRAS is a peripheral membrane protein that regulates multiple signaling pathways, and is mutated in ≈30 % of cancers. Transient self-association of KRAS is essential for activation of the downstream effector RAF and oncogenicity. The presence of anionic phosphatidylserine (PS) lipids in the membrane was shown to promote KRAS self-assembly, however, the structural mechanisms remain elusive. Here, we employed nanodisc bilayers of defined lipid compositions, and probed the impact of PS concentration on KRAS self-association. Paramagnetic NMR experiments demonstrated the existence of two transient dimer conformations involving alternate electrostatic contacts between R135 and either D153 or E168 on the "α4/5-α4/5" interface, and revealed that lipid composition and salt modulate their dynamic equilibrium. These dimer interfaces were validated by charge-reversal mutants. This plasticity demonstrates how the dynamic KRAS dimerization interface responds to the environment, and likely extends to the assembly of other signaling complexes on the membrane.

Consensus on the RAS dimerization hypothesis: Strong evidence for lipid-mediated clustering but not for G-domain-mediated interactions

One of the open questions in RAS biology is the existence of RAS dimers and their role in RAF dimerization and activation. The idea of RAS dimers arose from the discovery that RAF kinases function as obligate dimers, which generated the hypothesis that RAF dimer formation might be nucleated by G-domain-mediated RAS dimerization. Here, we review the evidence for RAS dimerization and describe a recent discussion among RAS researchers that led to a consensus that the clustering of two or more RAS proteins is not due to the stable association of G-domains but, instead, is a consequence of RAS C-terminal membrane anchors and the membrane phospholipids with which they interact.

Genetic variants of SOS2, MAP2K1 and RASGRF2 in the RAS pathway genes predict survival of HBV-related hepatocellular carcinoma patients

The RAS pathway participates in the cascade of proliferation and cell division process, and the activated RAS pathway can lead to tumorigenesis including hepatocellular carcinoma (HCC). However, few studies have explored the effects of genetic variants in the RAS pathway-related genes on the survival of patients with HBV-related HCC. In the present study, we assessed the associations between 11,658 single-nucleotide polymorphisms (SNPs) in 62 RAS pathway genes and the overall survival (OS) of 866 HBV-related HCC individuals, which were randomly split (1:1) into discovery and validation datasets. As a result, three potentially functional SNPs were identified, based on multivariable cox proportional hazards regression analyses, in SOS Ras/Rho guanine nucleotide exchange factor 2 (SOS2, rs4632055 A > G), Ras protein-specific guanine nucleotide releasing factor 2 (RASGRF2, rs26418A > G) and mitogen-activated protein kinase 1 (MAP2K1,rs57120695 C > T), which were significantly and independently associated with OS of HBV-related HCC patients [adjusted hazards ratios (HRs) of 1.42, 1.32 and 1.50, respectively; 95% confidence intervals (CI), 1.14 to 1.76, 1.15 to 1.53 and 1.15 to 1.97, respectively; P = 0.001, < 0.001 and 0.003, respectively]. Additionally, the joint effects as the unfavorable genotypes of these three SNPs showed a significant association with the poor survival of HCC (trend test P < 0.001). The expression quantitative trait loci (eQTL) analysis further revealed that the rs4632055 G allele and the rs26418 A allele were associated with lower mRNA expression levels of SOS2 and RASGRF2, respectively. Collectively, these potentially functional SNPs of RASGRF2, SOS2 and M2PAK1 may become potential prognostic biomarkers for HBV-related HCC after hepatectomy.

SOS1-inspired hydrocarbon-stapled peptide as a pan-Ras inhibitor

Blocking the interaction between Ras and Son of Sevenless homolog 1 (SOS1) has been an attractive therapeutic strategy for treating cancers involving oncogenic Ras mutations. K-Ras mutation is the most common in Ras-driven cancers, accounting for 86%, with N-Ras mutation and H-Ras mutation accounting for 11% and 3%, respectively. Here, we report the design and synthesis of a series of hydrocarbon-stapled peptides to mimic the alpha-helix of SOS1 as pan-Ras inhibitors. Among these stapled peptides, SSOSH-5 was identified to maintain a well-constrained alpha-helical structure and bind to H-Ras with high affinity. SSOSH-5 was furthermore validated to bind with Ras similarly to the parent linear peptide through structural modeling analysis. This optimized stapled peptide was proven to be capable of effectively inhibiting the proliferation of pan-Ras-mutated cancer cells and inducing apoptosis in a dose-dependent manner by modulating downstream kinase signaling. Of note, SSOSH-5 exhibited a high capability of crossing cell membranes and strong proteolytic resistance. We demonstrated that the peptide stapling strategy is a feasible approach for developing peptide-based pan-Ras inhibitors. Furthermore, we expect that SSOSH-5 can be further characterized and optimized for the treatment of Ras-driven cancers.

Classification of GTP-dependent K-Ras4B active and inactive conformational states

Classifying reliably active and inactive molecular conformations of wildtype (WT) and mutated oncogenic proteins is a key, ongoing challenge in molecular cancer studies. Here, we probe the GTP-bound K-Ras4B conformational dynamics using long-time atomistic molecular dynamics (MD) simulations. We extract and analyze the detailed underlying free energy landscape of WT K-Ras4B. We use two key reaction coordinates, labeled d₁ and d₂ (i.e., distances coordinating the P_β atom of the GTP ligand with two key residues, T35 and G60), shown to correlate closely with activities of WT and mutated K-Ras4B. However, our new K-Ras4B conformational kinetics study reveals a more complex network of equilibrium Markovian states. We show that a new reaction coordinate is required to account for the orientation of acidic K-Ras4B sidechains such as D38 with respect to the interface with binding effector RAF1 and rationalize the activation/inactivation propensities and the corresponding molecular binding mechanisms. We use this understanding to unveil how a relatively conservative mutation (i.e., D33E, in the switch I region) can lead to significantly different activation propensities compared with WT K-Ras4B. Our study sheds new light on the ability of residues near the K-Ras4B-RAF1 interface to modulate the network of salt bridges at the binding interface with the RAF1 downstream effector and, thus, to influence the underlying GTP-dependent activation/inactivation mechanism. Altogether, our hybrid MD-docking modeling approach enables the development of new in silico methods for quantitative assessment of activation propensity changes (e.g., due to mutations or local binding environment). It also unveils the underlying molecular mechanisms and facilitates the rational design of new cancer drugs.

Lifetime of actin-dependent protein nanoclusters

Protein nanoclusters (PNCs) are dynamic collections of a few proteins that spatially organize in nanometer-length clusters. PNCs are one of the principal forms of spatial organization of membrane proteins, and they have been shown or hypothesized to be important in various cellular processes, including cell signaling. PNCs show remarkable diversity in size, shape, and lifetime. In particular, the lifetime of PNCs can vary over a wide range of timescales. The diversity in size and shape can be explained by the interaction of the clustering proteins with the actin cytoskeleton or the lipid membrane, but very little is known about the processes that determine the lifetime of the nanoclusters. In this paper, using mathematical modeling of the cluster dynamics, we model the biophysical processes that determine the lifetime of actin-dependent PNCs. In particular, we investigated the role of actin aster fragmentation, which had been suggested to be a key determinant of the PNC lifetime, and we found that it is important only for a small class of PNCs. A simple extension of our model allowed us to investigate the kinetics of protein-ligand interaction near PNCs. We found an anomalous increase in the lifetime of ligands near PNCs, which agrees remarkably well with experimental data on RAS-RAF kinetics. In particular, analysis of the RAS-RAF data through our model provides falsifiable predictions and novel hypotheses that will not only shed light on the role of RAS-RAF kinetics in various cancers, but also will be useful in studying membrane protein clustering in general.

Dynamic regulation of RAS and RAS signaling

RAS proteins regulate most aspects of cellular physiology. They are mutated in 30% of human cancers and 4% of developmental disorders termed Rasopathies. They cycle between active GTP-bound and inactive GDP-bound states. When active, they can interact with a wide range of effectors that control fundamental biochemical and biological processes. Emerging evidence suggests that RAS proteins are not simple on/off switches but sophisticated information processing devices that compute cell fate decisions by integrating external and internal cues. A critical component of this compute function is the dynamic regulation of RAS activation and downstream signaling that allows RAS to produce a rich and nuanced spectrum of biological outputs. We discuss recent findings how the dynamics of RAS and its downstream signaling is regulated. Starting from the structural and biochemical properties of wild-type and mutant RAS proteins and their activation cycle, we examine higher molecular assemblies, effector interactions and downstream signaling outputs, all under the aspect of dynamic regulation. We also consider how computational and mathematical modeling approaches contribute to analyze and understand the pleiotropic functions of RAS in health and disease.

α4-α5 Helices on Surface of KRAS Can Accommodate Small Compounds That Increase KRAS Signaling While Inducing CRC Cell Death

KRAS is the most frequently mutated oncogene associated with the genesis and progress of pancreatic, lung and colorectal (CRC) tumors. KRAS has always been considered as a therapeutic target in cancer but until now only two compounds that inhibit one specific KRAS mutation have been approved for clinical use. In this work, by molecular dynamics and a docking process, we describe a new compound (P14B) that stably binds to a druggable pocket near the α4-α5 helices of the allosteric domain of KRAS. This region had previously been identified as the binding site for calmodulin (CaM). Using surface plasmon resonance and pulldown analyses, we prove that P14B binds directly to oncogenic KRAS thus competing with CaM. Interestingly, P14B favors oncogenic KRAS interaction with BRAF and phosphorylated C-RAF, and increases downstream Ras signaling in CRC cells expressing oncogenic KRAS. The viability of these cells, but not that of the normal cells, is impaired by P14B treatment. These data support the significance of the α4-α5 helices region of KRAS in the regulation of oncogenic KRAS signaling, and demonstrate that drugs interacting with this site may destine CRC cells to death by increasing oncogenic KRAS downstream signaling.

Mutations in the α4-α5 allosteric lobe of RAS do not significantly impair RAS signaling or self-association

Mutations in one of the three RAS genes (HRAS, KRAS, and NRAS) are present in nearly 20% of all human cancers. These mutations shift RAS to the GTP-loaded active state due to impairment in the intrinsic GTPase activity and disruption of GAP-mediated GTP hydrolysis, resulting in constitutive activation of effectors such as RAF. Because activation of RAF involves dimerization, RAS dimerization has been proposed as an important step in RAS-mediated activation of effectors. The α4-α5 allosteric lobe of RAS has been proposed as a RAS dimerization interface. Indeed, the NS1 monobody, which binds the α4-α5 region within the RAS G domain, inhibits RAS-dependent signaling and transformation as well as RAS nanoclustering at the plasma membrane. Although these results are consistent with a model in which the G domain dimerizes through the α4-α5 region, the isolated G domain of RAS lacks intrinsic dimerization capacity. Furthermore, prior studies analyzing α4-α5 point mutations have reported mixed effects on RAS function. Here, we evaluated the activity of a panel of single amino acid substitutions in the α4-α5 region implicated in RAS dimerization. We found that these proposed "dimerization-disrupting" mutations do not significantly impair self-association, signaling, or transformation of oncogenic RAS. These results are consistent with a model in which activated RAS protomers cluster in close proximity to promote the dimerization of their associated effector proteins (e.g., RAF) without physically associating into dimers mediated by specific molecular interactions. Our findings suggest the need for a nonconventional approach to developing therapeutics targeting the α4-α5 region.

NRas activity is regulated by dynamic interactions with nanoscale signaling clusters at the plasma membrane

NRas is a key mediator of the mitogenic pathway in normal cells and in cancer cells. Its dynamics and nanoscale organization at the plasma membrane (PM) facilitate its signaling. Here, we used two-color photoactivated localization microscopy to resolve the organization of individual NRas and associated signaling proteins in live melanoma cells, with resolution down to ∼20 nm. Upon EGF activation, a fraction of NRas and BRAF (dis)assembled synchronously at the PM in co-clusters. NRas and BRAF clusters associated with GPI-enriched domains, serving as possible nucleation sites for these clusters. NRas and BRAF association in mutual clusters was reduced by the NRas farnesylation inhibitor lonafarnib, yet enhanced by the BRAF inhibitor vemurafenib. Surprisingly, dispersed NRas molecules associated with the periphery of self-clusters of either Grb2 or NF1. Thus, NRas-mediated signaling, which is critical in health and disease, is regulated by dynamic interactions with functional clusters of BRAF or other related proteins at the PM.

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

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

Structural and biophysical properties of farnesylated KRas interacting with the chaperone SmgGDS-558

KRas is a small GTPase and membrane-bound signaling protein. Newly synthesized KRas is post-translationally modified with a membrane-anchoring prenyl group. KRas chaperones are therapeutic targets in cancer due to their participation in trafficking oncogenic KRas to membranes. SmgGDS splice variants are chaperones for small GTPases with basic residues in their hypervariable domain (HVR), including KRas. SmgGDS-607 escorts pre-prenylated small GTPases, while SmgGDS-558 escorts prenylated small GTPases. We provide a structural description of farnesylated and fully processed KRas (KRas-FMe) in complex with SmgGDS-558 and define biophysical properties of this interaction. Surface plasmon resonance measurements on biomimetic model membranes quantified the thermodynamics of the interaction of SmgGDS with KRas, and small-angle x-ray scattering was used to characterize complexes of SmgGDS-558 and KRas-FMe structurally. Structural models were refined using Monte Carlo and molecular dynamics simulations. Our results indicate that SmgGDS-558 interacts with the HVR and the farnesylated C-terminus of KRas-FMe, but not its G-domain. Therefore, SmgGDS-558 interacts differently with prenylated KRas than prenylated RhoA, whose G-domain was found in close contact with SmgGDS-558 in a recent crystal structure. Using immunoprecipitation assays, we show that SmgGDS-558 binds the GTP-bound, GDP-bound, and nucleotide-free forms of farnesylated and fully processed KRas in cells, consistent with SmgGDS-558 not engaging the G-domain of KRas. We found that the dissociation constant, Kd, for KRas-FMe binding to SmgGDS-558 is comparable with that for the KRas complex with PDEδ, a well-characterized KRas chaperone that also does not interact with the KRas G-domain. These results suggest that KRas interacts in similar ways with the two chaperones SmgGDS-558 and PDEδ. Therapeutic targeting of the SmgGDS-558/KRas complex might prove as useful as targeting the PDEδ/KRas complex in KRas-driven cancers.

Structure-function analysis of the SHOC2-MRAS-PP1C holophosphatase complex

Receptor tyrosine kinase (RTK)-RAS signalling through the downstream mitogen-activated protein kinase (MAPK) cascade regulates cell proliferation and survival. The SHOC2-MRAS-PP1C holophosphatase complex functions as a key regulator of RTK-RAS signalling by removing an inhibitory phosphorylation event on the RAF family of proteins to potentiate MAPK signalling1. SHOC2 forms a ternary complex with MRAS and PP1C, and human germline gain-of-function mutations in this complex result in congenital RASopathy syndromes2-5. However, the structure and assembly of this complex are poorly understood. Here we use cryo-electron microscopy to resolve the structure of the SHOC2-MRAS-PP1C complex. We define the biophysical principles of holoenzyme interactions, elucidate the assembly order of the complex, and systematically interrogate the functional consequence of nearly all of the possible missense variants of SHOC2 through deep mutational scanning. We show that SHOC2 binds PP1C and MRAS through the concave surface of the leucine-rich repeat region and further engages PP1C through the N-terminal disordered region that contains a cryptic RVXF motif. Complex formation is initially mediated by interactions between SHOC2 and PP1C and is stabilized by the binding of GTP-loaded MRAS. These observations explain how mutant versions of SHOC2 in RASopathies and cancer stabilize the interactions of complex members to enhance holophosphatase activity. Together, this integrative structure-function model comprehensively defines key binding interactions within the SHOC2-MRAS-PP1C holophosphatase complex and will inform therapeutic development .

Identification of a BRAF/PA28γ/MEK1 signaling axis and its role in epithelial-mesenchymal transition in oral submucous fibrosis

Oral submucous fibrosis (OSF) is a chronic and insidious oral potentially malignant disorder associated with a 4-17% risk of oral squamous cell carcinoma (OSCC). Our previous study found that proteasomal activator 28 gamma (PA28γ) is frequently overexpressed in oral squamous cell carcinoma and negatively correlated with poor patient prognosis. However, the role of PA28γ in the occurrence and development of OSF remains unclear. Here, we screened PA28γ-related genes and investigated their function in OSF. We demonstrated that the expression of PA28γ was positively associated with MEK1 and gradually elevated from normal to progressive stages of OSF tissue. Arecoline, a pathogenic component of OSF, could upregulate the protein levels of PA28γ and phosphorylated MEK1 and contribute to epithelial to mesenchymal transition (EMT) in epithelial cells. Notably, PA28γ could interact with MEK1 and upregulate its phosphorylation level. Furthermore, arecoline upregulated BRAF, which can interact with PA28γ and upregulate its protein level. Additionally, BRAF, PA28γ, and MEK1 could form protein complexes and then enhance the MEK1/ERK signaling pathways. The concrete mechanism of the protein stability of PA28γ is that BRAF mediates its degradation by inhibiting its ubiquitination. These findings underscore the instrumental role of PA28γ in the BRAF/MEK1 pathway and enhanced EMT through MEK1/ERK activation in OSF.

Galectin-3 Is a Crucial Immunological Disease Marker in Patients with Fungal Keratitis

Fungal keratitis, one of the most common infectious eye diseases in China, often results in a poor prognosis due to a delayed diagnosis and the insufficiency of effective therapy. There is an urgent need to identify specific biomarkers for the disease. In this study, we screened out tear proteins in patients with fungal keratitis by microsphere-based immunoassay analysis. Levels of cytokine expression were determined in both human corneal epithelial cell models in vitro and the corneas of patients by western blot, quantitative polymerase chain reaction (qPCR), and immunofluorescence analysis. Neutrophil activation was examined by flow cytometry analysis. The relationship between the cytokine expression and neutrophils was evaluated by immunofluorescence costaining and correlation analysis. These results demonstrated that the galectin-3 expression level was increased in both cell model and patient samples at the early and late stages of fungal keratitis. The neutrophils were significantly activated during the disease course of fungal keratitis. Meanwhile, colocalization and a positive correlation between galectin-3 and neutrophils were observed, suggesting that galectin-3 may play a crucial role in the recruitment of neutrophils and immune regulation of fungal keratitis. In conclusion, galectin-3 could be a key disease marker implying a beneficial immune response in the pathogenesis of fungal keratitis, which might be a target of therapeutic strategy in the future.

Therapeutic Targeting the Allosteric Cysteinome of RAS and Kinase Families

Allosteric mechanisms are pervasive in nature, but human-designed allosteric perturbagens are rare. The history of KRAS^G12C inhibitor development suggests that covalent chemistry may be a key to expanding the armamentarium of allosteric inhibitors. In that effort, irreversible targeting of a cysteine converted a non-deal allosteric binding pocket and low affinity ligands into a tractable drugging strategy. Here we examine the feasibility of expanding this approach to other allosteric pockets of RAS and kinase family members, given that both protein families are regulators of vital cellular processes that are often dysregulated in cancer and other human diseases. Moreover, these heavily studied families are the subject of numerous drug development campaigns that have resulted, sometimes serendipitously, in the discovery of allosteric inhibitors. We consequently conducted a comprehensive search for cysteines, a commonly targeted amino acid for covalent drugs, using AlphaFold-generated structures of those families. This new analysis presents potential opportunities for allosteric targeting of validated and understudied drug targets, with an emphasis on cancer therapy.

Structural insights into Ras regulation by SIN1

Both the mTORC2 and Ras-ERK pathways respond to growth factor stimulation and play critical roles in cell growth and proliferation, disarray of these pathways leads to many diseases, especially cancer. These two signaling pathways crosstalk at many levels; recently it's become clear that the SIN1 component of mTORC2 could interact with Ras family small GTPases, but how these two proteins interact at the molecular level and the functional outcomes of this interaction remain to be addressed. In this work we determined the high-resolution structure of Ras-SIN1 complexes and revealed the detailed interaction mechanism. We also showed that Ras-SIN1 association inhibits insulin-induced ERK activation. Insights from this work could improve our understanding of the disease-causing mechanism of errant mTORC2 or Ras proteins.

Over the years it has been established that SIN1, a key component of mTORC2, could interact with Ras family small GTPases through its Ras-binding domain (RBD). The physical association of Ras and SIN1/mTORC2 could potentially affect both mTORC2 and Ras-ERK pathways. To decipher the precise molecular mechanism of this interaction, we determined the high-resolution structures of HRas/KRas-SIN1 RBD complexes, showing the detailed interaction interface. Mutation of critical interface residues abolished Ras-SIN1 interaction and in SIN1 knockout cells we demonstrated that Ras-SIN1 association promotes SGK1 activity but inhibits insulin-induced ERK activation. With structural comparison and competition fluorescence resonance energy transfer (FRET) assays we showed that HRas-SIN1 RBD association is much weaker than HRas-Raf1 RBD but is slightly stronger than HRas-PI3K RBD interaction, providing a possible explanation for the different outcome of insulin or EGF stimulation. We also found that SIN1 isoform lacking the PH domain binds stronger to Ras than other longer isoforms and the PH domain appears to have an inhibitory effect on Ras-SIN1 binding. In addition, we uncovered a Ras dimerization interface that could be critical for Ras oligomerization. Our results advance our understanding of Ras-SIN1 association and crosstalk between growth factor-stimulated pathways.

Conformational control and regulation of the pseudokinase KSR via small molecule binding interactions

Pseudokinases often operate through functionally related enzymes and receptors. A prime example is the pseudokinase KSR (Kinase Suppressor of RAS), which can act as both an amplifier and inhibitor of members in the RAS-MAPK (Mitogen Activated Protein Kinase) signaling pathway. KSR is structurally related to the active RAF kinases over multiple domains; moreover, the pseudokinase domain of KSR forms physical and regulatory complexes with both RAF and MEK through distinct interfaces. Characterization of small molecule interactions on KSR has been used to uncover novel chemical tools and understand the mechanism of action of clinical drugs. Here, we elaborate on assays and structural methods for measuring binding at orthosteric and interfacial binding sites on KSR. These distinct small molecule pockets provide therapeutic paths for targeting KSR1 and KSR2 pseudokinases in disease, including in RAS and RAF mutant cancers.

Comparative Analysis and Isoform-Specific Therapeutic Vulnerabilities of KRAS Mutations in Non-Small Cell Lung Cancer

PURPOSE

Activating missense mutations of KRAS are the most frequent oncogenic driver events in lung adenocarcinoma (LUAD). However, KRAS isoforms are highly heterogeneous, and data on the potential isoform-dependent therapeutic vulnerabilities are still lacking.

EXPERIMENTAL DESIGN

We developed an isogenic cell-based platform to compare the oncogenic properties and specific therapeutic actionability of KRAS-mutant isoforms. In parallel, we analyzed clinicopathologic and genomic data from 3,560 patients with non-small cell lung cancer (NSCLC) to survey allele-specific features associated with oncogenic KRAS mutations.

RESULTS

In isogenic cell lines expressing different mutant KRAS isoforms, we identified isoform-specific biochemical, biological, and oncogenic properties both in vitro and in vivo. These exclusive features correlated with different therapeutic responses to MEK inhibitors, with KRAS G12C and Q61H mutants being more sensitive compared with other isoforms. In vivo, combined KRAS G12C and MEK inhibition was more effective than either drug alone. Among patients with NSCLCs that underwent comprehensive tumor genomic profiling, STK11 and ATM mutations were significantly enriched among tumors harboring KRAS G12C, G12A, and G12V mutations. KEAP1 mutation was significantly enriched among KRAS G12C and KRAS G13X LUADs. KRAS G13X-mutated tumors had the highest frequency of concurrent STK11 and KEAP1 mutations. Transcriptomic profiling revealed unique patterns of gene expression in each KRAS isoform, compared with KRAS wild-type tumors.

CONCLUSIONS

This study demonstrates that KRAS isoforms are highly heterogeneous in terms of concurrent genomic alterations and gene-expression profiles, and that stratification based on KRAS alleles should be considered in the design of future clinical trials.

Signalling dynamics, cell decisions, and homeostatic control in health and disease

Cell signalling engenders cells with the capability to receive and process information from the intracellular and extracellular environments, trigger and execute biological responses, and communicate with each other. Ultimately, cell signalling is responsible for maintaining homeostasis at the cellular, tissue and systemic level. For this reason, cell signalling is a topic of intense research efforts aimed to elucidate how cells coordinate transitions between states in developing and adult organisms in physiological and pathological conditions. Here, we review current knowledge of how cell signalling operates at multiple spatial and temporal scales, focusing on how single-cell analytical techniques reveal mechanisms underpinning cell-to-cell variability, signalling plasticity, and collective cellular responses.

A Structure is Worth a Thousand Words: New Insights for RAS and RAF Regulation

The RAS GTPases are frequently mutated in human cancer, with KRAS being the predominant tumor driver. For many years, it has been known that the structure and function of RAS are integrally linked, as structural changes induced by GTP binding or mutational events determine the ability of RAS to interact with regulators and effectors. Recently, a wealth of information has emerged from structures of specific KRAS mutants and from structures of multiprotein complexes containing RAS and/or RAF, an essential effector of RAS. These structures provide key insights regarding RAS and RAF regulation as well as promising new strategies for therapeutic intervention.

SIGNIFICANCE

The RAS GTPases are major drivers of tumorigenesis, and for RAS proteins to exert their full oncogenic potential, they must interact with the RAF kinases to initiate ERK cascade signaling. Although binding to RAS is typically a prerequisite for RAF to become an activated kinase, determining the molecular mechanisms by which this interaction results in RAF activation has been a challenging task. A major advance in understanding this process and RAF regulation has come from recent structural studies of various RAS and RAF multiprotein signaling complexes, revealing new avenues for drug discovery.

Identification of the nucleotide-free state as a therapeutic vulnerability for inhibition of selected oncogenic RAS mutants

RAS guanosine triphosphatases (GTPases) are mutated in nearly 20% of human tumors, making them an attractive therapeutic target. Following our discovery that nucleotide-free RAS (apo RAS) regulates cell signaling, we selectively target this state as an approach to inhibit RAS function. Here, we describe the R15 monobody that exclusively binds the apo state of all three RAS isoforms in vitro, regardless of the mutation status, and captures RAS in the apo state in cells. R15 inhibits the signaling and transforming activity of a subset of RAS mutants with elevated intrinsic nucleotide exchange rates (i.e., fast exchange mutants). Intracellular expression of R15 reduces the tumor-forming capacity of cancer cell lines driven by select RAS mutants and KRAS(G12D)-mutant patient-derived xenografts (PDXs). Thus, our approach establishes an opportunity to selectively inhibit a subset of RAS mutants by targeting the apo state with drug-like molecules.

Khan et al. develop a high-affinity monobody to nucleotide-free RAS that, when expressed intracellularly, inhibits oncogenic RAS-mediated signaling and tumorigenesis. This study reveals the feasibility of targeting the nucleotide-free state to inhibit tumors driven by oncogenic RAS mutants that possess elevated nucleotide exchange activity.

Ras Multimers on the Membrane: Many Ways for a Heart-to-Heart Conversation

Formation of Ras multimers, including dimers and nanoclusters, has emerged as an exciting, new front of research in the ‘old’ field of Ras biomedicine. With significant advances made in the past few years, we are beginning to understand the structure of Ras multimers and, albeit preliminary, mechanisms that regulate their formation in vitro and in cells. Here we aim to synthesize the knowledge accrued thus far on Ras multimers, particularly the presence of multiple globular (G-) domain interfaces, and discuss how membrane nanodomain composition and structure would influence Ras multimer formation. We end with some general thoughts on the potential implications of Ras multimers in basic and translational biology.