Significance: Here we present an unprecedented multiscale simulation platform that enables modeling, hypothesis generation, and discovery across biologically relevant length and time scales to predict mechanisms that can be tested experimentally. We demonstrate that our predictive simulation-experimental validation loop generates accurate insights into RAS-membrane biology. Evaluating over 100,000 correlated simulations, we show that RAS–lipid interactions are dynamic and evolving, resulting in: 1) a reordering and selection of lipid domains in realistic eight-lipid bilayers, 2) clustering of RAS into multimers correlating with specific lipid fingerprints, 3) changes in the orientation of the RAS G-domain impacting its ability to interact with effectors, and 4) demonstration that RAS–RAS G-domain interfaces are nonspecific in these putative signaling domains.

Approved inhibitors of KRASG12C prevent oncogenic activation by sequestering the inactive, GDP-bound (OFF) form rather than directly binding and inhibiting the active, GTP-bound (ON) form. This approach provides no direct target coverage of the active protein. Expectedly, adaptive resistance to KRASG12C (OFF)-only inhibitors is observed in association with increased expression and activity of KRASG12C(ON). To provide optimal KRASG12C target coverage, we have developed BBO-8520, a first-in-class, direct dual inhibitor of KRASG12C(ON) and (OFF) forms. BBO-8520 binds in the Switch-II/Helix3 pocket, covalently modifies the target cysteine, and disables effector binding to KRASG12C(ON). BBO-8520 exhibits potent signaling inhibition in growth factor–activated states, in which current (OFF)-only inhibitors demonstrate little measurable activity. In vivo, BBO-8520 demonstrates rapid target engagement and inhibition of signaling, resulting in durable tumor regression in multiple models, including those resistant to KRASG12C(OFF)-only inhibitors. BBO-8520 is in phase 1 clinical trials in patients with KRASG12C non–small cell lung cancer. Significance: BBO-8520 is a first-in-class direct, small molecule covalent dual inhibitor that engages KRASG12C in the active (ON) and inactive (OFF) conformations. BBO-8520 represents a novel mechanism of action that allows for optimal target coverage and delays the emergence of adaptive resistance seen with (OFF)-only inhibitors in the clinic. See related commentary by Zhou and Westover, p. 455

Significance: The proto-oncogene KRAS , a small GTPase, is frequently mutated in pancreatic, colorectal, and lung cancer. These mutations result in elevated levels of the activated guanosine triphosphate-bound form of KRAS. Localized at the plasma membrane, KRAS functions to recruit effectors, predominantly RAF kinase for activation and initiation of the MAPK signaling cascade. Combining computational and biophysical methods we identify a membrane-distal state of the KRAS G-domain that alternates with two previously described membrane-proximal states through dynamic reorganization of the hypervariable region. Comprising about 90% of the ensemble, this membrane-distal state of the G-domain dominates the proximal states and may facilitate KRAS to recruit cytosolic RAF kinase to the membrane by a fly-casting mechanism.

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

KRAS , the most frequently mutated oncogene in human cancer, produces two isoforms, KRAS4a and KRAS4b, through alternative splicing. These isoforms differ in exon 4, which encodes the final 15 residues of the G-domain and hypervariable regions (HVRs), vital for trafficking and membrane localization. While KRAS4b has been extensively studied, KRAS4a has been largely overlooked. Our multidisciplinary study compared the structural and functional characteristics of KRAS4a and KRAS4b, revealing distinct structural properties and thermal stability. Position 151 influences KRAS4a’s thermal stability, while position 153 affects binding to RAF1 CRD protein. Nuclear magnetic resonance analysis identified localized structural differences near sequence variations and provided a solution-state conformational ensemble. Notably, KRAS4a exhibits substantial transcript abundance in bile ducts, liver, and stomach, with transcript levels approaching KRAS4b in the colon and rectum. Functional disparities were observed in full-length KRAS variants, highlighting the impact of HVR variations on interaction with trafficking proteins and downstream effectors like RAF and PI3K within cells.