The disordered hypervariable region and the folded catalytic domain of oncogenic K-Ras4B partner in phospholipid binding

Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation

K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.

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.

Recapitulation of cell-like KRAS4b membrane dynamics on complex biomimetic membranes

Understanding the spatiotemporal distribution and dynamics of RAS on the plasma membrane (PM) is the key for elucidating the molecular mechanisms of the RAS signaling pathway. Single particle tracking (SPT) experiments show that in cells, KRAS diffuses in at least three interchanging states on the cellular PM; however, KRAS remains monomeric and always shows homogeneous diffusion on artificial membranes. Here, we show for the first time on a supported lipid bilayer composed of heterogeneous lipid components that we can recapitulate the three-state diffusion of KRAS seen in cells. The use of a biologically relevant eight-lipid system opens a new frontier in the biophysical studies of RAS and other membrane associated proteins on a biomimetic system that recapitulates the complexity of a cellular PM.

•

KRAS4b shows homogeneous diffusion on simple 2-lipids bilayer

•

KRAS4b shows a cell-like, three-state diffusion on a complex 8-lipid bilayer

•

Phase separation in lipids favors the multi-state diffusion of KRAS4b

•

The complex lipid composition favors RAS nanoclustering irrespective of nucleotide state

Ras isoform-specific expression, chromatin accessibility, and signaling

The anchorage of Ras isoforms in the membrane and their nanocluster formations have been studied extensively, including their detailed interactions, sizes, preferred membrane environments, chemistry, and geometry. However, the staggering challenge of their epigenetics and chromatin accessibility in distinct cell states and types, which we propose is a major factor determining their specific expression, still awaits unraveling. Ras isoforms are distinguished by their C-terminal hypervariable region (HVR) which acts in intracellular transport, regulation, and membrane anchorage. Here, we review some isoform-specific activities at the plasma membrane from a structural dynamic standpoint. Inspired by physics and chemistry, we recognize that understanding functional specificity requires insight into how biomolecules can organize themselves in different cellular environments. Within this framework, we suggest that isoform-specific expression may largely be controlled by the chromatin density and physical compaction, which allow (or curb) access to "chromatinized DNA." Genes are preferentially expressed in tissues: proteins expressed in pancreatic cells may not be equally expressed in lung cells. It is the rule-not an exception, and it can be at least partly understood in terms of chromatin organization and accessibility state. Genes are expressed when they can be sufficiently exposed to the transcription machinery, and they are less so when they are persistently buried in dense chromatin. Notably, chromatin accessibility can similarly determine expression of drug resistance genes.

Unveiling the Dynamics of KRAS4b on Lipid Model Membranes

Small GTPase proteins are ubiquitous and responsible for regulating several processes related to cell growth and differentiation. Mutations that stabilize their active state can lead to uncontrolled cell proliferation and cancer. Although these proteins are well characterized at the cellular scale, the molecular mechanisms governing their functions are still poorly understood. In addition, there is limited information about the regulatory function of the cell membrane which supports their activity. Thus, we have studied the dynamics and conformations of the farnesylated KRAS4b in various membrane model systems, ranging from binary fluid mixtures to heterogeneous raft mimics. Our approach combines long time-scale coarse-grained (CG) simulations and Markov state models to dissect the membrane-supported dynamics of KRAS4b. Our simulations reveal that protein dynamics is mainly modulated by the presence of anionic lipids and to some extent by the nucleotide state (activation) of the protein. In addition, our results suggest that both the farnesyl and the polybasic hypervariable region (HVR) are responsible for its preferential partitioning within the liquid-disordered (Ld) domains in membranes, potentially enhancing the formation of membrane-driven signaling platforms.

Mutational landscape of K-Ras substitutions at 12th position-a systematic molecular dynamics approach

K-Ras is a small GTPase and acts as a molecular switch by recruiting GEFs and GAPs, and alternates between the inert GDP-bound and the dynamic GTP-bound forms. The amino acid at position 12 of K-Ras is a hot spot for oncogenic mutations (G12A, G12C, G12D, G12R, G12S, and G12V), disturbing the active fold of the protein, leading to cancer development. This study aimed to investigate the potential conformational changes induced by these oncogenic mutations at the 12^th position, impairing GAP-mediated GTP hydrolysis. Comprehensive computational tools (iStable, FoldX, SNPeffect, DynaMut, and CUPSAT) were used to evaluate the effect of these six mutations on the stability of wild type K-Ras protein. The docking of GTP with K-Ras was carried out using AutoDock4.2, followed by molecular dynamics simulations. Furthermore, on comparison of binding energies between the wild type K-Ras and the six mutants, we have demonstrated that the G12A and G12V mutants exhibited the strongest binding efficiency compared to the other four mutants. Trajectory analyses of these mutations revealed that G12A encountered the least deviation, fluctuation, intermolecular H-bonds, and compactness compared to the wildtype, which was supported by the lower Gibbs free energy value. Our study investigates the molecular dynamics simulations of the mutant K-Ras forms at the 12^th position, which expects to provide insights about the molecular mechanisms involved in cancer development, and may serve as a platform for targeted therapies against cancer. Communicated by Ramaswamy H. Sarma.

Multivalent assembly of KRAS with the RAS-binding and cysteine-rich domains of CRAF on the membrane

Membrane anchoring of farnesylated KRAS is critical for activation of RAF kinases, yet our understanding of how these proteins interact on the membrane is limited to isolated domains. The RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF engage KRAS and the plasma membrane, unleashing the kinase domain from autoinhibition. Due to experimental challenges, structural insight into this tripartite KRAS:RBD-CRD:membrane complex has relied on molecular dynamics simulations. Here, we report NMR studies of the KRAS:CRAF RBD-CRD complex. We found that the nucleotide-dependent KRAS-RBD interaction results in transient electrostatic interactions between KRAS and CRD, and we mapped the membrane interfaces of the CRD, RBD-CRD, and the KRAS:RBD-CRD complex. RBD-CRD exhibits dynamic interactions with the membrane through the canonical CRD lipid-binding site (CRD β7-8), as well as an alternative interface comprising β6 and the C terminus of CRD and β2 of RBD. Upon complex formation with KRAS, two distinct states were observed by NMR: State A was stabilized by membrane association of CRD β7-8 and KRAS α4-α5 while state B involved the C terminus of CRD, β3-5 of RBD, and part of KRAS α5. Notably, α4-α5, which has been proposed to mediate KRAS dimerization, is accessible only in state B. A cancer-associated mutation on the state B membrane interface of CRAF RBD (E125K) stabilized state B and enhanced kinase activity and cellular MAPK signaling. These studies revealed a dynamic picture of the assembly of the KRAS-CRAF complex via multivalent and dynamic interactions between KRAS, CRAF RBD-CRD, and the membrane.

A Non-Canonical Calmodulin Target Motif Comprising a Polybasic Region and Lipidated Terminal Residue Regulates Localization

Calmodulin (CaM) is a Ca²⁺-sensor that regulates a wide variety of target proteins, many of which interact through short basic helical motifs bearing two hydrophobic ‘anchor’ residues. CaM comprises two globular lobes, each containing a pair of EF-hand Ca²⁺-binding motifs that form a Ca²⁺-induced hydrophobic pocket that binds an anchor residue. A central flexible linker allows CaM to accommodate diverse targets. Several reported CaM interactors lack these anchors but contain Lys/Arg-rich polybasic sequences adjacent to a lipidated N- or C-terminus. Ca²⁺-CaM binds the myristoylated N-terminus of CAP23/NAP22 with intimate interactions between the lipid and a surface comprised of the hydrophobic pockets of both lobes, while the basic residues make electrostatic interactions with the negatively charged surface of CaM. Ca²⁺-CaM binds farnesylcysteine, derived from the farnesylated polybasic C-terminus of KRAS4b, with the lipid inserted into the C-terminal lobe hydrophobic pocket. CaM sequestration of the KRAS4b farnesyl moiety disrupts KRAS4b membrane association and downstream signaling. Phosphorylation of basic regions of N-/C-terminal lipidated CaM targets can reduce affinity for both CaM and the membrane. Since both N-terminal myristoylated and C-terminal prenylated proteins use a Singly Lipidated Polybasic Terminus (SLIPT) for CaM binding, we propose these polybasic lipopeptide elements comprise a non-canonical CaM-binding motif.

Calmodulin disrupts plasma membrane localization of farnesylated KRAS4b by sequestering its lipid moiety

KRAS4b is a small guanosine triphosphatase (GTPase) protein that regulates several signal transduction pathways that underlie cell proliferation, differentiation, and survival. KRAS4b function requires prenylation of its C terminus and recruitment to the plasma membrane, where KRAS4b activates effector proteins including the RAF family of kinases. The Ca²⁺-sensing protein calmodulin (CaM) has been suggested to regulate the localization of KRAS4b through direct, Ca²⁺-dependent interaction, but how CaM and KRAS4b functionally interact is controversial. Here, we determined a crystal structure, which was supported by solution nuclear magnetic resonance (NMR), that revealed the sequestration of the prenyl moiety of KRAS4b in the hydrophobic pocket of the C-terminal lobe of Ca²⁺-bound CaM. Our engineered fluorescence resonance energy transfer (FRET)-based biosensor probes (CaMeRAS) showed that, upon stimulation of Ca²⁺ influx by extracellular ligands, KRAS4b reversibly translocated in a Ca²⁺-CaM-dependent manner from the plasma membrane to the cytoplasm in live HeLa and HEK293 cells. These results reveal a mechanism underlying the inhibition of KRAS4b activity by Ca²⁺ signaling pathways.

The quaternary assembly of KRas4B with Raf-1 at the membrane

Proximally located in the membrane, oncogenic Ras dimers (or nanoclusters) can recruit and promote Raf dimerization and MAPK (Raf/MEK/ERK) signaling. Among Ras isoforms, KRas4B is the most frequently mutated. Recent data on the binary KRas4B–Raf-1 complex suggested that Raf-1 CRD not only executes membrane anchorage, but also supports the high-affinity interaction of Raf-1 RBD with KRas4B catalytic domain. For a detailed mechanistic picture of Raf activation at the membrane, we employ explicit MD simulations of the quaternary KRas4B–Raf-1 complex. The complex contains two active GTP-bound KRas4B proteins forming a dimer through the allosteric lobe interface and two tandem RBD-CRD segments of Raf-1 interacting with the effector lobes at both ends of the KRas4B dimer. We show that Raf-1 RBD-CRD supports stable KRas4B dimer at preferred interface and orientation at the membrane, thereby cooperatively enhancing the affinity of the KRas4B–Raf-1 interaction. We propose that a Ras dimer at the membrane can increase the population of proximal Raf kinase domains, promoting kinase domain dimerization in the cytoplasm. Collectively, the dynamic Ras–Raf assembly promotes Raf activation not by allostery; instead, Ras activates Raf by shifting its ensemble toward kinase domain-accessible states through enhanced affinity at the membrane.

The Hypervariable Region of K-Ras4B Governs Molecular Recognition and Function

The flexible C-terminal hypervariable region distinguishes K-Ras4B, an important proto-oncogenic GTPase, from other Ras GTPases. This unique lysine-rich portion of the protein harbors sites for post-translational modification, including cysteine prenylation, carboxymethylation, phosphorylation, and likely many others. The functions of the hypervariable region are diverse, ranging from anchoring K-Ras4B at the plasma membrane to sampling potentially auto-inhibitory binding sites in its GTPase domain and participating in isoform-specific protein-protein interactions and signaling. Despite much research, there are still many questions about the hypervariable region of K-Ras4B. For example, mechanistic details of its interaction with plasma membrane lipids and with the GTPase domain require further clarification. The roles of the hypervariable region in K-Ras4B-specific protein-protein interactions and signaling are incompletely defined. It is also unclear why post-translational modifications frequently found in protein polylysine domains, such as acetylation, glycation, and carbamoylation, have not been observed in K-Ras4B. Expanding knowledge of the hypervariable region will likely drive the development of novel highly-efficient and selective inhibitors of K-Ras4B that are urgently needed by cancer patients.

Novel compound cedrelone inhibits hepatocellular carcinoma progression via PBLD and Ras/Rap1

Although it is known that Phenazine biosynthesis-like domain-containing protein (PBLD) expression is downregulated in hepatocellular carcinoma (HCC), its biological function is unclear. Additionally, no agents capable of upregulating PBLD exist. In the current study, the relationship between PBLD and HCC was analyzed using clinicopathological specimens. A HCC cell model, microarray analysis and an animal model were used to verify the therapeutic effect of cedrelone on HCC. The present study demonstrated that PBLD inhibited HCC progression. Furthermore, the present study revealed that cedrelone possessed treated-HCC capabilities via targeted PBLD overexpression. The epithelial-mesenchymal transition phenotype and growth rate were inhibited and the apoptosis ratio was promoted by cedrelone following PBLD overexpression. The Ras and Ras-proximate-1 signaling pathways were also determined to be regulated by cedrelone via PBLD activation in HCC. PBLD may therefore be an independent predictor of HCC progression and a novel target for HCC treatment. Additionally, the PBLD activator, cedrelone, may be a potential drug for HCC treatment in the future.

Phosphatidic acid in membrane rearrangements

Phosphatidic acid (PA) is the simplest cellular glycerophospholipid characterized by unique biophysical properties: a small headgroup; negative charge; and a phosphomonoester group. Upon interaction with lysine or arginine, PA charge increases from -1 to -2 and this change stabilizes protein-lipid interactions. The biochemical properties of PA also allow interactions with lipids in several subcellular compartments. Based on this feature, PA is involved in the regulation and amplification of many cellular signalling pathways and functions, as well as in membrane rearrangements. Thereby, PA can influence membrane fusion and fission through four main mechanisms: it is a substrate for enzymes producing lipids (lysophosphatidic acid and diacylglycerol) that are involved in fission or fusion; it contributes to membrane rearrangements by generating negative membrane curvature; it interacts with proteins required for membrane fusion and fission; and it activates enzymes whose products are involved in membrane rearrangements. Here, we discuss the biophysical properties of PA in the context of the above four roles of PA in membrane fusion and fission.

The structural basis for Ras activation of PI3Kα lipid kinase

PI3Kα is a principal Ras effector that phosphorylates PIP2 to PIP3 in the PI3K/Akt/mTOR pathway. How Ras activates PI3K has been unclear: is Ras' role confined to PI3K recruitment to the membrane or does Ras activation also involve allostery? Recently, we determined the mechanism of PI3Kα activation at the atomic level. We showed the vital role and significance of conformational change in PI3Kα activation. Here, by a 'best-match for hydrogen-bonding pair' (BMHP) computational protocol and molecular dynamics (MD) simulations, we model the atomic structure of KRas4B in complex with the Ras binding domain (RBD) of PI3Kα, striving to understand the mechanism of PI3Kα activation by Ras. Point mutations T208D, K210E, and K227E disrupt the KRas4B-RBD interface in the models, in line with the experiments. We identify allosteric signaling pathways connecting Ras to RBD in the p110α subunit. However, the observed weak allosteric signals coupled with the detailed mechanism of PI3Kα activation make us conclude that the dominant mechanistic role of Ras is likely to be recruitment and restriction of the PI3Kα population at the membrane. Thus, RTK recruits the PI3Kα to the membrane and activates it by relieving its autoinhibition exerted by the nSH2 domain, leading to exposure of the kinase domain, which permits PIP2 binding. Ras recruitment can shift the PI3Kα ensemble toward a population where the kinase domain surface and the active site position and orientation favor PIP2 insertion. This work helps elucidate Ras-mediated PI3K activation and explores the structural basis for Ras-PI3Kα drug discovery.

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

K-Ras G-domain binding with signaling lipid phosphatidylinositol (4,5)-phosphate (PIP2): membrane association, protein orientation, and function

Ras genes potently drive human cancers, with mutated proto-oncogene GTPase KRAS4B (K-Ras4B) being the most abundant isoform. Targeted inhibition of oncogenic gene products is considered the "holy grail" of present-day cancer therapy, and recent discoveries of small-molecule KRas4B inhibitors were made thanks to a deeper understanding of the structure and dynamics of this GTPase. Because interactions with biological membranes are key for Ras function, Ras-lipid interactions have become a major focus, especially because such interactions evidently involve both the Ras C terminus for lipid anchoring and its G-protein domain. Here, using NMR spectroscopy and molecular dynamics simulations complemented by biophysical- and cell-biology assays, we investigated the interaction between K-Ras4B with the signaling lipid phosphatidylinositol (4,5)-phosphate (PIP2). We discovered that the β2 and β3 strands as well as helices 4 and 5 of the GTPase G-domain bind to PIP2 and identified the specific residues in these structural elements employed in these interactions, likely occurring in two K-Ras4B orientation states relative to the membrane. Importantly, we found that some of these residues known to be oncogenic when mutated (D47K, D92N, K104M, and D126N) are critical for K-Ras-mediated transformation of fibroblast cells, but do not substantially affect basal and assisted nucleotide hydrolysis and exchange. Moreover, the K104M substitution abolished localization of K-Ras to the plasma membrane. The findings suggest that specific G-domain residues can critically regulate Ras function by mediating interactions with membrane-associated PIP2 lipids; these insights that may inform the future design of therapeutic reagents targeting Ras activity.

Calmodulin (CaM) Activates PI3Kα by Targeting the "Soft" CaM-Binding Motifs in Both the nSH2 and cSH2 Domains of p85α

PI3Kα is a key lipid kinase in the PI3K/Akt pathway. Its frequent oncogenic mutations make it a primary drug target. Calmodulin (CaM) activates PI3Kα independently of extracellular signals, indicating a significant role in oncogenic PI3Kα activation. Here, we reveal the atomic-scale structures of CaM in complexes with the nSH2 and cSH2 domains of the regulatory p85α subunit of PI3Kα, and illustrate how CaM activates PI3Kα by targeting the "soft 1-5-10" CaM-binding motifs in both nSH2 and cSH2 domains. Experiment observed CaM binding cSH2 first, followed by nSH2 binding hours later. CaM typically prefers binding helical peptides. Here we observe that, unlike in cSH2, the CaM-binding motif in nSH2 populates a mixed β-sheet/α-helix/random coil structure. The population shift from a β-sheet toward CaM's favored α-helical conformation explains why the nSH2 domain needs a longer time for CaM binding in the experiments. The "soft" CaM-binding motifs in both nSH2 and cSH2 domains establish strong CaM-PI3Kα interactions, collectively facilitating PI3Kα activation. This work uncovers the structural basis for CaM-driven PI3Kα activation.

Autoinhibition in Ras effectors Raf, PI3Kα, and RASSF5: a comprehensive review underscoring the challenges in pharmacological intervention

Autoinhibition is an effective mechanism that guards proteins against spurious activation. Despite its ubiquity, the distinct organizations of the autoinhibited states and their release mechanisms differ. Signaling is most responsive to the cell environment only if a small shift in the equilibrium is required to switch the system from an inactive (occluded) to an active (exposed) state. Ras signaling follows this paradigm. This underscores the challenge in pharmacological intervention to exploit and enhance autoinhibited states. Here, we review autoinhibition and release mechanisms at the membrane focusing on three representative Ras effectors, Raf protein kinase, PI3Kα lipid kinase, and NORE1A (RASSF5) tumor suppressor, and point to the ramifications to drug discovery. We further touch on Ras upstream and downstream signaling, Ras activation, and the Ras superfamily in this light, altogether providing a broad outlook of the principles and complexities of autoinhibition.

Cyclase-associated protein 1 (CAP1) is a prenyl-binding partner of Rap1 GTPase

Rap1 proteins are members of the Ras subfamily of small GTPases involved in many biological responses, including adhesion, cell proliferation, and differentiation. Like all small GTPases, they work as molecular allosteric units that are active in signaling only when associated with the proper membrane compartment. Prenylation, occurring in the cytosol, is an enzymatic posttranslational event that anchors small GTPases at the membrane, and prenyl-binding proteins are needed to mask the cytoplasm-exposed lipid during transit to the target membrane. However, several of these proteins still await discovery. In this study, we report that cyclase-associated protein 1 (CAP1) binds Rap1. We found that this binding is GTP-independent, does not involve Rap1's effector domain, and is fully contained in its C-terminal hypervariable region (HVR). Furthermore, Rap1 prenylation was required for high-affinity interactions with CAP1 in a geranylgeranyl-specific manner. The prenyl binding specifically involved CAP1's C-terminal hydrophobic β-sheet domain. We present a combination of experimental and computational approaches, yielding a model whereby the high-affinity binding between Rap1 and CAP1 involves electrostatic and nonpolar side-chain interactions between Rap1's HVR residues, lipid, and CAP1 β-sheet domain. The binding was stabilized by the lipid insertion into the β-solenoid whose interior was occupied by nonpolar side chains. This model was reminiscent of the recently solved structure of the PDEδ-K-Ras complex; accordingly, disruptors of this complex, e.g. deltarasin, blocked the Rap1-CAP1 interaction. These findings indicate that CAP1 is a geranylgeranyl-binding partner of Rap1.

Raf-1 Cysteine-Rich Domain Increases the Affinity of K-Ras/Raf at the Membrane, Promoting MAPK Signaling

K-Ras4B preferentially activates Raf-1. The high-affinity interaction of Ras-binding domain (RBD) of Raf with Ras was solved, but the relative position of Raf's cysteine-rich domain (CRD) in the Ras/Raf complex at the membrane and key question of exactly how it affects Raf signaling are daunting. We show that CRD stably binds anionic membranes inserting a positively charged loop into the amphipathic interface. Importantly, when in complex with Ras/RBD, covalently connected CRD presents the same membrane interaction mechanism, with CRD locating at the space between the RBD and membrane. To date, CRD's role was viewed in terms of stabilizing Raf-membrane interaction. Our observations argue for a key role in reducing Ras/RBD fluctuations at the membrane, thereby increasing Ras/RBD affinity. Even without K-Ras, via CRD, Raf-1 can recruit to the membrane; however, by reducing the Ras/RBD fluctuations and enhancing Ras/RBD affinity at the membrane, CRD promotes Raf's activation and MAPK signaling over other pathways.

Oncogenic Ras Isoforms Signaling Specificity at the Membrane

How do Ras isoforms attain oncogenic specificity at the membrane? Oncogenic KRas, HRas, and NRas (K-Ras, H-Ras, and N-Ras) differentially populate distinct cancers. How they selectively activate effectors and why is KRas4B the most prevalent are highly significant questions. Here, we consider determinants that may bias isoform-specific effector activation and signaling at the membrane. We merge functional data with a conformational view to provide mechanistic insight. Cell-specific expression levels, pathway cross-talk, and distinct interactions are the key, but conformational trends can modulate selectivity. There are two major pathways in oncogenic Ras-driven proliferation: MAPK (Raf/MEK/ERK) and PI3Kα/Akt/mTOR. All membrane-anchored, proximally located, oncogenic Ras isoforms can promote Raf dimerization and fully activate MAPK signaling. So why the differential statistics of oncogenic isoforms in distinct cancers and what makes KRas so highly oncogenic? Many cell-specific factors may be at play, including higher KRAS mRNA levels. As a key factor, we suggest that because only KRas4B binds calmodulin, only KRas can fully activate PI3Kα/Akt signaling. We propose that full activation of both MAPK and PI3Kα/Akt proliferative pathways by oncogenic KRas4B-but not by HRas or NRas-may help explain why the KRas4B isoform is especially highly populated in certain cancers. We further discuss pharmacologic implications. Cancer Res; 78(3); 593-602. ©2017 AACR.

Oncogenic Ras Isoforms Signaling Specificity at the Membrane

How do Ras isoforms attain oncogenic specificity at the membrane? Oncogenic KRas, HRas, and NRas (K-Ras, H-Ras, and N-Ras) differentially populate distinct cancers. How they selectively activate effectors and why is KRas4B the most prevalent are highly significant questions. Here, we consider determinants that may bias isoform-specific effector activation and signaling at the membrane. We merge functional data with a conformational view to provide mechanistic insight. Cell-specific expression levels, pathway cross-talk, and distinct interactions are the key, but conformational trends can modulate selectivity. There are two major pathways in oncogenic Ras-driven proliferation: MAPK (Raf/MEK/ERK) and PI3Kα/Akt/mTOR. All membrane-anchored, proximally located, oncogenic Ras isoforms can promote Raf dimerization and fully activate MAPK signaling. So why the differential statistics of oncogenic isoforms in distinct cancers and what makes KRas so highly oncogenic? Many cell-specific factors may be at play, including higher KRAS mRNA levels. As a key factor, we suggest that because only KRas4B binds calmodulin, only KRas can fully activate PI3Kα/Akt signaling. We propose that full activation of both MAPK and PI3Kα/Akt proliferative pathways by oncogenic KRas4B-but not by HRas or NRas-may help explain why the KRas4B isoform is especially highly populated in certain cancers. We further discuss pharmacologic implications. Cancer Res; 78(3); 593-602. ©2017 AACR.

Is Nanoclustering essential for all oncogenic KRas pathways? Can it explain why wild-type KRas can inhibit its oncogenic variant?

Membrane-anchored oncogenic KRas can dimerize, form nanoclusters, and signal through the MAPK (Raf/MEK/ERK) and PI3Kα/Akt/mTOR. Both pathways are needed in KRAS-driven proliferation. Here we ask: Is oncogenic KRas nanoclustering (or dimerization) essential for all KRas signaling pathways? Raf kinase domain dimerization, thus MAPK activation, requires KRas nanoclusters. By contrast, the PI3Kα heterodimer acts as a monomeric unit; thus, does PI3Kα activation and PI3Kα/Akt/mTOR signaling require nanoclustering? Further, calmodulin binds only to oncogenic KRas4B. Here we ask: Does calmodulin downregulate KRas4B cancer development as suggested early on, or promote it? We also ask: Why is oncogenic KRas4B the most abundant isoform? Does wild-type Ras indeed inhibit its oncogenic variants as data appeared to suggest? And related to the last question, why is wild-type KRas a more potent inhibitor of its oncogenic form than wild-type NRas of its oncogenic form? Resolving these cardinal questions, and others, such as how exactly does RASSF5 (NORE1A) act as tumor suppressor, and why Ras isoforms tend to occur in distinct cancer types are crucial for effective pharmacology. In this review, we take a nanoclustering/dimerization-centric outlook and show that many questions can be explained by simply considering Ras nanoclustering.

Phosphorylated Calmodulin Promotes PI3K Activation by Binding to the SH2 Domains

How calmodulin (CaM) acts in KRAS-driven cancers is a vastly important question. CaM binds to and stimulates PI3Kα/Akt signaling, promoting cell growth and proliferation. Phosphorylation of CaM at Tyr99 (pY99) enhances PI3Kα activation. PI3Kα is a lipid kinase. It phosphorylates PIP2 to produce PIP3, to which Akt binds. PI3Kα has two subunits: the regulatory p85 and the catalytic p110. Here, exploiting explicit-solvent MD simulations we unveil key interactions between phosphorylated CaM (pCaM) and the two SH2 domains in the p85 subunit, confirm experimental observations, and uncover PI3Kα's mechanism of activation. pCaMs form strong and stable interactions with both nSH2 and cSH2 domains, with pY99 being the dominant contributor. Despite the high structural similarity between the two SH2 domains, we observe that nSH2 prefers an extended CaM conformation, whereas cSH2 prefers a collapsed conformation. Notably, collapsed CaM is observed after binding of an extended CaM to K-Ras4B. Thus, the more populated extended pCaM conformation targets nSH2 to release its autoinhibition of p110 catalytic sites. This executes the key activation step of PI3Kα. Independently, K-Ras4B allosterically activates p110. These events are at the cell membrane, which contributes to tighten the PI3Kα Ras binding domain/K-Ras4B interaction, to accomplish K-Ras4B allosteric activation, with a minor contribution from cSH2.

Molecular Dynamics Simulations and Dynamic Network Analysis Reveal the Allosteric Unbinding of Monobody to H-Ras Triggered by R135K Mutation

Ras proteins, as small GTPases, mediate cell proliferation, survival and differentiation. Ras mutations have been associated with a broad spectrum of human cancers and thus targeting Ras represents a potential way forward for cancer therapy. A recently reported monobody NS1 allosterically disrupts the Ras-mediated signaling pathway, but its efficacy is reduced by R135K mutation in H-Ras. However, the detailed mechanism is unresolved. Here, using molecular dynamics (MD) simulations and dynamic network analysis, we explored the molecular mechanism for the unbinding of NS1 to H-Ras and shed light on the underlying allosteric network in H-Ras. MD simulations revealed that the overall structures of the two complexes did not change significantly, but the H-Ras–NS1 interface underwent significant conformational alteration in the mutant Binding free energy analysis showed that NS1 binding was unfavored after R135K mutation, which resulted in the unfavorable binding of NS1. Furthermore, the critical residues on H-Ras responsible for the loss of binding of NS1 were identified. Importantly, the allosteric networks for these important residues were revealed, which yielded a novel insight into the allosteric regulatory mechanism of H-Ras.

Intrinsic protein disorder in oncogenic KRAS signaling

How Ras, and in particular its most abundant oncogenic isoform K-Ras4B, is activated and signals in proliferating cells, poses some of the most challenging questions in cancer cell biology. In this paper, we ask how intrinsically disordered regions in K-Ras4B and its effectors help promote proliferative signaling. Conformational disorder allows spanning long distances, supports hinge motions, promotes anchoring in membranes, permits segments to fulfil multiple roles, and broadly is crucial for activation mechanisms and intensified oncogenic signaling. Here, we provide an overview illustrating some of the key mechanisms through which conformational disorder can promote oncogenesis, with K-Ras4B signaling serving as an example. We discuss (1) GTP-bound KRas4B activation through membrane attachment; (2) how farnesylation and palmitoylation can promote isoform functional specificity; (3) calmodulin binding and PI3K activation; (4) how Ras activates its RASSF5 cofactor, thereby stimulating signaling of the Hippo pathway and repressing proliferation; and (5) how intrinsically disordered segments in Raf help its attachment to the membrane and activation. Collectively, we provide the first inclusive review of the roles of intrinsic protein disorder in oncogenic Ras-driven signaling. We believe that a broad picture helps to grasp and formulate key mechanisms in Ras cancer biology and assists in therapeutic intervention.

PDEδ Binding to Ras Isoforms Provides a Route to Proper Membrane Localization

To signal, Ras isoforms must be enriched at the plasma membrane (PM). It was suggested that phosphodiesterase-δ (PDEδ) can bind and shuttle some farnesylated Ras isoforms to the PM, but not all. Among these, interest focused on K-Ras4B, the most abundant oncogenic Ras isoform. To study PDEδ/Ras interactions, we modeled and simulated the PDEδ/K-Ras4B complex. We obtained structures, which were similar to two subsequently determined crystal structures. We next modeled and simulated complexes of PDEδ with the farnesylated hypervariable regions of K-Ras4A and N-Ras. Earlier data suggested that PDEδ extracts K-Ras4B and N-Ras from the PM, but surprisingly not K-Ras4A. Earlier analysis of the crystal structures advanced that the presence of large/charged residues adjacent to the farnesylated site precludes the PDEδ interaction. Here, we show that PDEδ can bind to farnesylated K-Ras4A and N-Ras like K-Ras4B, albeit not as strongly. This weaker binding, coupled with the stronger anchoring of K-Ras4A in the membrane (but not of electrostatically neutral N-Ras), can explain the observation why PDEδ is unable to effectively extract K-Ras4A. We thus propose that farnesylated Ras isoforms can bind PDEδ to fulfill the required PM enrichment, and argue that the different environments, PM versus solution, can resolve apparently puzzling Ras observations. These are novel insights that would not be expected based on the crystal structures alone, which provide an elegant rationale for previously puzzling observations of the differential effects of PDEδ on farnesylated Ras family proteins.

Calmodulin and PI3K Signaling in KRAS Cancers

Calmodulin (CaM) uniquely promotes signaling of oncogenic K-Ras; but not N-Ras or H-Ras. How CaM interacts with K-Ras and how this stimulates cell proliferation are among the most challenging questions in KRAS-driven cancers. Earlier data pointed to formation of a ternary complex consisting of K-Ras, PI3Kα and CaM. Recent data point to phosphorylated CaM binding to the SH2 domains of the p85 subunit of PI3Kα and activating it. Modeling suggests that the high affinity interaction between the phosphorylated CaM tyrosine motif and PI3Kα, can promote full PI3Kα activation by oncogenic K-Ras. Our up-to-date review discusses CaM's role in PI3K signaling at the membrane in KRAS-driven cancers. This is significant since it may help development of K-Ras-specific pharmacology.

Intrinsic K-Ras dynamics: A novel molecular dynamics data analysis method shows causality between residue pair motions

K-Ras is the most frequently mutated oncogene in human cancers, but there are still no drugs that directly target it in the clinic. Recent studies utilizing dynamics information show promising results for selectively targeting mutant K-Ras. However, despite extensive characterization, the mechanisms by which K-Ras residue fluctuations transfer allosteric regulatory information remain unknown. Understanding the direction of information flow can provide new mechanistic insights for K-Ras targeting. Here, we present a novel approach –conditional time-delayed correlations (CTC) – using the motions of all residue pairs of a protein to predict directionality in the allosteric regulation of the protein fluctuations. Analyzing nucleotide-dependent intrinsic K-Ras motions with the new approach yields predictions that agree with the literature, showing that GTP-binding stabilizes K-Ras motions and leads to residue correlations with relatively long characteristic decay times. Furthermore, our study is the first to identify driver-follower relationships in correlated motions of K-Ras residue pairs, revealing the direction of information flow during allosteric modulation of its nucleotide-dependent intrinsic activity: active K-Ras Switch-II region motions drive Switch-I region motions, while α-helix-3L7 motions control both. Our results provide novel insights for strategies that directly target mutant K-Ras.

Drugging Ras GTPase: a comprehensive mechanistic and signaling structural view

Ras proteins are small GTPases, cycling between inactive GDP-bound and active GTP-bound states. Through these switches they regulate signaling that controls cell growth and proliferation. Activating Ras mutations are associated with approximately 30% of human cancers, which are frequently resistant to standard therapies. Over the past few years, structural biology and in silico drug design, coupled with improved screening technology, led to a handful of promising inhibitors, raising the possibility of drugging Ras proteins. At the same time, the invariable emergence of drug resistance argues for the critical importance of additionally honing in on signaling pathways which are likely to be involved. Here we overview current advances in Ras structural knowledge, including the conformational dynamic of full-length Ras in solution and at the membrane, therapeutic inhibition of Ras activity by targeting its active site, allosteric sites, and Ras-effector protein-protein interfaces, Ras dimers, the K-Ras4B/calmodulin/PI3Kα trimer, and targeting Ras with siRNA. To mitigate drug resistance, we propose signaling pathways that can be co-targeted along with Ras and explain why. These include pathways leading to the expression (or activation) of YAP1 and c-Myc. We postulate that these and Ras signaling pathways, MAPK/ERK and PI3K/Akt/mTOR, act independently and in corresponding ways in cell cycle control. The structural data are instrumental in the discovery and development of Ras inhibitors for treating RAS-driven cancers. Together with the signaling blueprints through which drug resistance can evolve, this review provides a comprehensive and innovative master plan for tackling mutant Ras proteins.

Membrane-associated Ras dimers are isoform-specific: K-Ras dimers differ from H-Ras dimers

Are the dimer structures of active Ras isoforms similar? This question is significant since Ras can activate its effectors as a monomer; however, as a dimer, it promotes Raf's activation and MAPK (mitogen-activated protein kinase) cell signalling. In the present study, we model possible catalytic domain dimer interfaces of membrane-anchored GTP-bound K-Ras4B and H-Ras, and compare their conformations. The active helical dimers formed by the allosteric lobe are isoform-specific: K-Ras4B-GTP favours the α3 and α4 interface; H-Ras-GTP favours α4 and α5. Both isoforms also populate a stable β-sheet dimer interface formed by the effector lobe; a less stable β-sandwich interface is sustained by salt bridges of the β-sheet side chains. Raf's high-affinity β-sheet interaction is promoted by the active helical interface. Collectively, Ras isoforms' dimer conformations are not uniform; instead, the isoform-specific dimers reflect the favoured interactions of the HVRs (hypervariable regions) with cell membrane microdomains, biasing the effector-binding site orientations, thus isoform binding selectivity.

The Structural Basis of Oncogenic Mutations G12, G13 and Q61 in Small GTPase K-Ras4B

Ras mediates cell proliferation, survival and differentiation. Mutations in K-Ras4B are predominant at residues G12, G13 and Q61. Even though all impair GAP-assisted GTP → GDP hydrolysis, the mutation frequencies of K-Ras4B in human cancers vary. Here we aim to figure out their mechanisms and differential oncogenicity. In total, we performed 6.4 μs molecular dynamics simulations on the wild-type K-Ras4B (K-Ras4B(WT)-GTP/GDP) catalytic domain, the K-Ras4B(WT)-GTP-GAP complex, and the mutants (K-Ras4B(G12C/G12D/G12V)-GTP/GDP, K-Ras4B(G13D)-GTP/GDP, K-Ras4B(Q61H)-GTP/GDP) and their complexes with GAP. In addition, we simulated 'exchanged' nucleotide states. These comprehensive simulations reveal that in solution K-Ras4B(WT)-GTP exists in two, active and inactive, conformations. Oncogenic mutations differentially elicit an inactive-to-active conformational transition in K-Ras4B-GTP; in K-Ras4B(G12C/G12D)-GDP they expose the bound nucleotide which facilitates the GDP-to-GTP exchange. These mechanisms may help elucidate the differential mutational statistics in K-Ras4B-driven cancers. Exchanged nucleotide simulations reveal that the conformational transition is more accessible in the GTP-to-GDP than in the GDP-to-GTP exchange. Importantly, GAP not only donates its R789 arginine finger, but stabilizes the catalytically-competent conformation and pre-organizes catalytic residue Q61; mutations disturb the R789/Q61 organization, impairing GAP-mediated GTP hydrolysis. Together, our simulations help provide a mechanistic explanation of key mutational events in one of the most oncogenic proteins in cancer.

Comparison of the Conformations of KRAS Isoforms, K-Ras4A and K-Ras4B, Points to Similarities and Significant Differences

Human HRAS, KRAS, and NRAS genes encode four isoforms of Ras, a p21 GTPase. Mutations in KRAS account for the majority of RAS-driven cancers. The KRAS has two splice variants, K-Ras4A and K-Ras4B. Due to their reversible palmitoylation, K-Ras4A and N-Ras have bimodal signaling states. K-Ras4A and K-Ras4B differ in four catalytic domain residues (G151R/D153E/K165Q/H166Y) and in their disordered C-terminal hypervariable region (HVR). In K-Ras4A, the HVR is not as strongly positively charged as in K-Ras4B (+6e vs +9e). Here, we performed all-atom molecular dynamics simulations to elucidate isoform-specific differences between the two splice variants. We observe that the catalytic domain of GDP-bound K-Ras4A has a more exposed nucleotide binding pocket than K-Ras4B, and the dynamic fluctuations in switch I and II regions also differ; both factors may influence guanine-nucleotide exchange. We further observe that like K-Kas4B, full-length K-Ras4A exhibits nucleotide-dependent HVR fluctuations; however, these fluctuations differ between the GDP-bound forms of K-Ras4A and K-Ras4B. Unlike K-Ras4B where the HVR tends to cover the effector binding region, in K-Ras4A, autoinhibited states are unstable. With lesser charge, the K-Ras4A HVR collapses on itself, making it less available for binding the catalytic domain. Since the HVRs of N- and H-Ras are weakly charged (+1e and +2e, respectively), autoinhibition may be a unique feature of K-Ras4B.

Ras Conformational Ensembles, Allostery, and Signaling

Ras proteins are classical members of small GTPases that function as molecular switches by alternating between inactive GDP-bound and active GTP-bound states. Ras activation is regulated by guanine nucleotide exchange factors that catalyze the exchange of GDP by GTP, and inactivation is terminated by GTPase-activating proteins that accelerate the intrinsic GTP hydrolysis rate by orders of magnitude. In this review, we focus on data that have accumulated over the past few years pertaining to the conformational ensembles and the allosteric regulation of Ras proteins and their interpretation from our conformational landscape standpoint. The Ras ensemble embodies all states, including the ligand-bound conformations, the activated (or inactivated) allosteric modulated states, post-translationally modified states, mutational states, transition states, and nonfunctional states serving as a reservoir for emerging functions. The ensemble is shifted by distinct mutational events, cofactors, post-translational modifications, and different membrane compositions. A better understanding of Ras biology can contribute to therapeutic strategies.