K-Ras4B/calmodulin/PI3Kα: A promising new adenocarcinoma-specific drug target?

Novel approaches for the development of direct KRAS inhibitors: structural insights and drug design

INTRODUCTION

Hyperactivated RAS signaling is reported in 13% of all human cancers, in which ~80% resulted due to KRAS mutations alone. Direct inhibition of KRAS is an important aspect in treating KRAS-related tumors. Despite the efforts of more than four decades, not many KRAS inhibitors have been successful in obtaining clinical approval, except the very recent FDA approval for sotorasib. In recent years, the understanding of structural insights and allosteric pocket identification at catalytic sites of KRAS are likely to provide an excellent opportunity for the development of much more effective clinical candidates.

AREA COVERED

The presented review article mainly summarizes the developments of small molecule KRAS inhibitors as drug candidates and rational approaches that are being utilized for the selective targeting of KRAS signaling in the mutant cancer cells.

EXPERT OPINION

After the initial success in targeting the mutant KRAS G12C variants, the search has been shifted to address the challenges concerning the resistance and efficacy of small molecule KRAS inhibitors. However, the contribution of other KRAS mutations at G12V, G13C, and G13D variants causing cancers is much higher than the mutations at G12C. In view of this aspect, specific attention is required to target all other mutations as well. Accordingly, for the development of KRAS targeted therapies, the design of small molecule inhibitors that can inhibit KRAS signaling and as well as target inhibition of other signaling pathways like RAS-SOS and RAS-PI3K has to be explored extensively.

Advances in Intracellular Calcium Signaling Reveal Untapped Targets for Cancer Therapy

Intracellular Ca2+ distribution is a tightly regulated process. Numerous Ca2+ chelating, storage, and transport mechanisms are required to maintain normal cellular physiology. Ca2+-binding proteins, mainly calmodulin and calbindins, sequester free intracellular Ca2+ ions and apportion or transport them to signaling hubs needing the cations. Ca2+ channels, ATP-driven pumps, and exchangers assist the binding proteins in transferring the ions to and from appropriate cellular compartments. Some, such as the endoplasmic reticulum, mitochondria, and lysosomes, act as Ca2+ repositories. Cellular Ca2+ homeostasis is inefficient without the active contribution of these organelles. Moreover, certain key cellular processes also rely on inter-organellar Ca2+ signaling. This review attempts to encapsulate the structure, function, and regulation of major intracellular Ca2+ buffers, sensors, channels, and signaling molecules before highlighting how cancer cells manipulate them to survive and thrive. The spotlight is then shifted to the slow pace of translating such research findings into anticancer therapeutics. We use the PubMed database to highlight current clinical studies that target intracellular Ca2+ signaling. Drug repurposing and improving the delivery of small molecule therapeutics are further discussed as promising strategies for speeding therapeutic development in this area.

Novel PI3K/Akt/mTOR pathway inhibitors plus radiotherapy: Strategy for non-small cell lung cancer with mutant RAS gene

Non-small cell lung cancer (NSCLC) with RAS -mutant gene has been the most difficult obstacle to overcome. Over 25% of muted lung adenocarcinomas have RAS mutation. The prognosis of NSCLC patients with RAS-mutant genes is always poor because there is no effective drug to suppress RAS-mutant genes. NSCLC patients with RAS-mutant usually develop resistance to radiotherapy and chemotherapy, which in some cases leads to a 5-10% survival rate for non-small cell lung cancer (NSCLC). As little clinical symptom of NSCLC was presented at its early stages, thus it always brings in disappointing treatment outcome. Currently, NSCLC presents the highest morbidity and mortality all over the world. The combination of PI3K/AKT/mTOR pathway inhibitors with radiotherapy is a novel strategy to improve radiosensitivity and therapeutic outcome of NSCLC with a RAS-mutant gene. There have been many preclinical studies and clinical trials on the effect of PI3K/AKT/mTOR pathway inhibitors combined with radiotherapy in NSCLC with a RAS-mutant gene have been reported in the past years. This review provides current knowledge of the combination of PI3K/Akt/mTOR pathway inhibitors with radiotherapy, which prove to be a significant improvement for the treatment of NSCLC patients with RAS mutations and will benefit NSCLC patients with RAS mutations.

The Structural Basis of the Farnesylated and Methylated KRas4B Interaction with Calmodulin

Ca²⁺-calmodulin (CaM) extracts KRas4B from the plasma membrane, suggesting that KRas4B/CaM interaction plays a role in regulating Ras signaling. To gain mechanistic insight, we provide a computational model, supported by experimental structural data, of farnesylated/methylated KRas4B_1-185 interacting with CaM in solution and at anionic membranes including signaling lipids. Due to multiple interaction modes, we observe diverse conformational ensembles of the KRas4B-CaM complex. A highly populated conformation reveals the catalytic domain interacting with the N-lobe and the hypervariable region (HVR) wrapping around the linker with the farnesyl docking to the extended CaM's C-lobe pocket. Alternatively, KRas4B can interact with collapsed CaM with the farnesyl penetrating CaM's center. At anionic membranes, CaM interacts with the catalytic domain with large fluctuations, drawing the HVR. Signaling lipids establishing strong salt bridges with CaM prevent membrane departure. Membrane-interacting KRas4B-CaM complex can productively recruit phosphatidylinositol 3-kinase α (PI3Kα) to the plasma membrane, serving as a coagent in activating PI3Kα/Akt signaling.

Drugging K-RasG12C through covalent inhibitors: Mission possible?

Ras, whose mutants are present in approximately 30% of human tumours, is one of the most important oncogenes. Drugging Ras is thus regarded as the quest for the Holy Grail in cancer therapeutics development. Despite more than three decades of efforts, drug discovery targeting Ras constantly fails, rendering Ras undruggable, due to its smooth surface and picomolar affinity towards guanosine substrates. The most frequently mutated isoform of Ras is K-Ras, accounting for >85% of Ras-driven cancers, and one majority of them is the G12C mutation. Recent advances in structural biology shed light on drugging Ras, and one of the cutting-edge breakthroughs is the design of covalent G12C-specific inhibitors targeting the mutated cysteine. This type of inhibitor can be classified into substrate-competitive orthosteric inhibitors and non-competitive allosteric inhibitors. They display improved selectivity and enhanced potency due to their G12-specific and irreversible covalent binding nature. Thus, they represent a new hope for revolutionizing the conventional characterization of Ras as "undruggable" and pave a promising avenue for further drug discovery. Here, we provide comprehensive structural and medicinal chemical insights into K-Ras covalent inhibitors specific for the G12C mutant. We first present an in-depth analysis of the conformations of the inhibitor binding pockets. Then, all the latest covalent ligands selectively inhibiting K-Ras^G12C are reviewed. Finally, we examine the current challenges faced by this new class of anti-Ras inhibitors.

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.

The mechanism of PI3Kα activation at the atomic level

PI3K lipid kinases phosphorylate PIP₂ to PIP₃ in the PI3K/Akt/mTOR pathway to regulate cellular processes. They are frequently mutated in cancer. Here we determine the PI3Kα activation mechanism at the atomic level. Unlike protein kinases where the substrate abuts the ATP, crystal structures indicate that in PI3Kα, the distance between the γ phosphate of the ATP and the PIP₂ lipid substrate is over 6 Å, much too far for the phosphoryl transfer, raising the question of how catalysis is executed. PI3Kα has two subunits, the catalytic p110α and the regulatory p85α. Our simulations show that release of the autoinhibition exerted by the nSH2 domain of the p85α triggers significant conformational change in p110α, leading to the exposure of the kinase domain for membrane interaction. Structural rearrangement in the C-lobe of the kinase domain reduces the distance between the ATP γ-phosphate and the substrate, offering an explanation as to how phosphoryl transfer is executed. An alternative mechanism may involve ATP relocation. This mechanism not only explains how oncogenic mutations promote PI3Kα activation by facilitating nSH2 release, or nSH2-release-induced, allosteric motions; it also offers an innovative, PI3K isoform-specific drug discovery principle. Rather than competing with nanomolar range ATP in the ATP-binding pocket and contending with ATP pocket conservation and massive binding targets, this mechanism suggests blocking the PI3Kα sequence-specific cavity between the ATP-binding pocket and the substrate binding site. Targeting isoform-specific residues in the cavity may prevent PIP₂ phosphorylation.

Computational Structural Biology: Successes, Future Directions, and Challenges

Computational biology has made powerful advances. Among these, trends in human health have been uncovered through heterogeneous 'big data' integration, and disease-associated genes were identified and classified. Along a different front, the dynamic organization of chromatin is being elucidated to gain insight into the fundamental question of genome regulation. Powerful conformational sampling methods have also been developed to yield a detailed molecular view of cellular processes. when combining these methods with the advancements in the modeling of supramolecular assemblies, including those at the membrane, we are finally able to get a glimpse into how cells' actions are regulated. Perhaps most intriguingly, a major thrust is on to decipher the mystery of how the brain is coded. Here, we aim to provide a broad, yet concise, sketch of modern aspects of computational biology, with a special focus on computational structural biology. We attempt to forecast the areas that computational structural biology will embrace in the future and the challenges that it may face. We skirt details, highlight successes, note failures, and map directions.

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.

Interaction of Calmodulin with the cSH2 Domain of the p85 Regulatory Subunit

Calmodulin (CaM) is a calcium sensor protein that directly interacts with the dual-specificity (lipid and protein) kinase PI3Kα through the SH2 domains of the p85 regulatory subunit. In adenocarcinomas, the CaM interaction removes the autoinhibition of the p110 catalytic subunit of PI3Kα, leading to activation of PI3Kα and promoting cell proliferation, survival, and migration. Here we demonstrate that the cSH2 domain of p85α engages its two CaM-binding motifs in the interaction with the N- and C-lobes of CaM as well as the flexible central linker, and our nuclear magnetic resonance experiments provide structural details. We show that in response to binding CaM, cSH2 exposes its tryptophan residue at the N-terminal region to the solvent. Because of the flexible nature of both CaM and cSH2, multiple binding modes of the interactions are possible. Binding of CaM to the cSH2 domain can help release the inhibition imposed on the p110 subunit, similar to the binding of the phosphorylated motif of RTK, or phosphorylated CaM (pCaM), to the SH2 domains. Amino acid sequence analysis shows that CaM-binding motifs are common in SH2 domains of non-RTKs. We speculate that CaM can also activate these kinases through similar mechanisms.

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.

Computational Insights into the Interactions between Calmodulin and the c/nSH2 Domains of p85α Regulatory Subunit of PI3Kα: Implication for PI3Kα Activation by Calmodulin

Calmodulin (CaM) and phosphatidylinositide-3 kinase (PI3Kα) are well known for their multiple roles in a series of intracellular signaling pathways and in the progression of several human cancers. Crosstalk between CaM and PI3Kα has been an area of intensive research. Recent experiments have shown that in adenocarcinoma, K-Ras4B is involved in the CaM-PI3Kα crosstalk. Based on experimental results, we have recently put forward a hypothesis that the coordination of CaM and PI3Kα with K-Ras4B forms a CaM-PI3Kα-K-Ras4B ternary complex, which leads to the formation of pancreatic ductal adenocarcinoma. However, the mechanism for the CaM-PI3Kα crosstalk is unresolved. Based on molecular modeling and molecular dynamics simulations, here we explored the potential interactions between CaM and the c/nSH2 domains of p85α subunit of PI3Kα. We demonstrated that CaM can interact with the c/nSH2 domains and the interaction details were unraveled. Moreover, the possible modes for the CaM-cSH2 and CaM-nSH2 interactions were uncovered and we used them to construct a complete CaM-PI3Kα complex model. The structural model of CaM-PI3Kα interaction not only offers a support for our previous ternary complex hypothesis, but also is useful for drug design targeted at CaM-PI3Kα protein-protein interactions.

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.

Modeling of RAS complexes supports roles in cancer for less studied partners

Background

RAS protein interactions have predominantly been studied in the context of the RAF and PI3kinase oncogenic pathways. Structural modeling and X-ray crystallography have demonstrated that RAS isoforms bind to canonical downstream effector proteins in these pathways using the highly conserved switch I and II regions. Other non-canonical RAS protein interactions have been experimentally identified, however it is not clear whether these proteins also interact with RAS via the switch regions.

Results

To address this question we constructed a RAS isoform-specific protein-protein interaction network and predicted 3D complexes involving RAS isoforms and interaction partners to identify the most probable interaction interfaces. The resulting models correctly captured the binding interfaces for well-studied effectors, and additionally implicated residues in the allosteric and hyper-variable regions of RAS proteins as the predominant binding site for non-canonical effectors. Several partners binding to this new interface (SRC, LGALS1, RABGEF1, CALM and RARRES3) have been implicated as important regulators of oncogenic RAS signaling. We further used these models to investigate competitive binding and multi-protein complexes compatible with RAS surface occupancy and the putative effects of somatic mutations on RAS protein interactions.

Conclusions

We discuss our findings in the context of RAS localization to the plasma membrane versus within the cytoplasm and provide a list of RAS protein interactions with possible cancer-related consequences, which could help guide future therapeutic strategies to target RAS proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s13628-017-0037-6) contains supplementary material, which is available to authorized users.

Flexible-body motions of calmodulin and the farnesylated hypervariable region yield a high-affinity interaction enabling K-Ras4B membrane extraction

In calmodulin (CaM)-rich environments, oncogenic KRAS plays a critical role in adenocarcinomas by promoting PI3K/Akt signaling. We previously proposed that at elevated calcium levels in cancer, CaM recruits PI3Kα to the membrane and extracts K-Ras4B from the membrane, organizing a K-Ras4B-CaM-PI3Kα ternary complex. CaM can thereby replace a missing receptor-tyrosine kinase signal to fully activate PI3Kα. Recent experimental data show that CaM selectively promotes K-Ras signaling but not of N-Ras or H-Ras. How CaM specifically targets K-Ras and how it extracts it from the membrane in KRAS-driven cancer is unclear. Obtaining detailed structural information for a CaM-K-Ras complex is still challenging. Here, using molecular dynamics simulations and fluorescence experiments, we observed that CaM preferentially binds unfolded K-Ras4B hypervariable regions (HVRs) and not α-helical HVRs. The interaction involved all three CaM domains including the central linker and both lobes. CaM specifically targeted the highly polybasic anchor region of the K-Ras4B HVR that stably wraps around CaM's acidic linker. The docking of the farnesyl group to the hydrophobic pockets located at both CaM lobes further enhanced CaM-HVR complex stability. Both CaM and K-Ras4B HVR are highly flexible molecules, suggesting that their interactions permit highly dynamic flexible-body motions. We, therefore, anticipate that the flexible-body interaction is required to extract K-Ras4B from the membrane, as conformational plasticity enables CaM to orient efficiently to the polybasic HVR anchor, which is partially diffused into the liquid-phase membrane. Our structural model of the CaM-K-Ras4B HVR association provides plausible clues to CaM's regulatory action in PI3Kα activation involving the ternary complex in cell proliferation signaling by oncogenic K-Ras.

CBP501 inhibits EGF-dependent cell migration, invasion and epithelial-to-mesenchymal transition of non-small cell lung cancer cells by blocking KRas to calmodulin binding

The anti-cancer agent CBP501 binds to calmodulin (CaM). Recent studies showed that migration and metastasis are inhibited by several CaM antagonists. However, there is no available evidence that CBP501 has similar effects. Here we found that CBP501 inhibits migration of non-small cell lung cancer (NSCLC) cells in vitro, even in the presence of migration inducing factors such as WNT, IL-6, and several growth factors. CBP501 also inhibited epidermal growth factor (EGF) enhanced invasion and the epithelial-to-mesenchymal transition (EMT), and this inhibition was accompanied by (i) suppression of Akt and ERK1/2 phosphorylation, and (ii) suppression of expression of transcription factor Zeb1 and the mesenchymal marker Vimentin. A pull down analysis performed using sepharose-immobilized CaM showed that CBP501 blocks the interaction between CaM and KRas. Furthermore, EGF induced Akt activation and cell migration was effectively suppressed by KRas down-regulation in NSCLC cells. Stable knockdown of KRas also made cells insensitive to CBP501's inhibition of growth factor-induced migration. Taken together, these results indicate that CBP501 inhibits binding of CaM with KRas and thereby suppresses the PI3K/AKT pathway, migration, invasion and EMT. These findings have identified a previously unrecognized effect of CBP501 on downstream KRas signaling mechanisms involving EMT and invasion, and provide support for the further clinical development of this agent.

A New View of Pathway-Driven Drug Resistance in Tumor Proliferation

Defeating drug resistance in tumor cell proliferation is challenging. We propose that signaling in cell proliferation takes place via two core pathways, each embodying multiple alternative pathways. We consider drug resistance through an alternative proliferation pathway - within the same or within the other core pathway. Most drug combinations target only one core pathway; blocking both can restrain proliferation. We define core pathways as independent and acting similarly in cell-cycle control, which can explain why their products (e.g., ERK and YAP1) can substitute for each other in resistance. Core pathways can forecast possible resistance because acquired resistance frequently occurs through alternative proliferation pathways. This concept may help to predict the efficacy of drug combinations. The selection of distinct combinations for specific mutated pathways would be guided by clinical diagnosis.

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.

The dynamic mechanism of RASSF5 and MST kinase activation by Ras

As a tumor suppressor, RASSF5 (NORE1A) activates MST1/2 thereby modulating the Hippo pathway. Structurally, activation involves RASSF5 and MST1/2 swapping their SARAH domains to form a SARAH heterodimer. This exposes the MST1/2 kinase domain which homodimerizes, leading to trans-autophosphorylation. The SARAH-SARAH interaction shifts RASSF5 away from its autoinhibited state and relieves MST1/2 autoinhibition. Separate crystal structures are available for the RA (Ras association) domain and SARAH dimer, where SARAH is a long straight α-helix. Using all-atom molecular dynamics simulations, we modeled the RASSF5 RA with a covalently connected SARAH to elucidate the dynamic mechanism of how SARAH mediates between autoinhibition and Ras triggered-activation. Our results show that in inactive RASSF5 the RA domain retains SARAH, yielding a self-associated conformation in which SARAH is in a kinked α-helical motif that increases the binding interface. When RASSF5 binds K-Ras4B-GTP, the equilibrium shifts toward SARAH's interacting with MST. Since the RA/SARAH affinity is relatively low, whereas that of the SARAH heterodimer is in the nM range, we suggest that RASSF5 exerts its tumor suppressor action through competition with other Ras effectors for Ras effector binding site, as well as coincidentally its recruitment to the membrane to help MST activation. Thus, SARAH plays a key role in RASSF5's tumor suppression action by linking the two major pathways in tumor cell proliferation: Ras and the MAPK (tumor cell proliferation-promoting) pathway, and the Hippo (tumor cell proliferation-suppressing) pathway.

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.

Oncogenic KRAS signaling and YAP1/β-catenin: Similar cell cycle control in tumor initiation

Why are YAP1 and c-Myc often overexpressed (or activated) in KRAS-driven cancers and drug resistance? Here, we propose that there are two independent pathways in tumor proliferation: one includes MAPK/ERK and PI3K/A kt/mTOR; and the other consists of pathways leading to the expression (or activation) of YAP1 and c-Myc. KRAS contributes through the first. MYC is regulated by e.g. β-catenin, Notch and Hedgehog. We propose that YAP1 and ERK accomplish similar roles in cell cycle control, as do β-catenin and PI3K. This point is compelling, since the question of how YAP1 rescues K-Ras or B-Raf ablation has recently captured much attention, as well as the mechanism of resistance to PI3K inhibitors. The similarity in cell cycle actions of β-catenin and PI3K can also clarify the increased aggressiveness of lung cancer when both K-Ras and β-catenin operate. Thus, we propose that the two pathways can substitute one another - or together amplify each other - in promoting proliferation. This new understanding of the independence and correspondence of the two pathways in cancer - MAPK/ERK and PI3K/Akt/mTOR; and YAP1 and c-Myc - provide a coherent and significant picture of signaling-driven oncogenic proliferation and may help in judicious, pathway-based drug discovery.

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.