Biophysical and Biochemical Characterization of Structurally Diverse Small Molecule Hits for KRAS Inhibition.

We describe six compounds as early hits for the development of direct inhibitors of KRAS, an important anticancer drug target. We show that these compounds bind to KRAS with affinities in the low micromolar range and exert different effects on its interactions with binding partners. Some of the compounds exhibit selective binding to the activated form of KRAS and inhibit signal transduction through both the MAPK or the phosphatidylinositide 3-kinase PI3K-protein kinase B (AKT) pathway in cells expressing mutant KRAS. Most inhibit intrinsic and/or SOS-mediated KRAS activation while others inhibit RAS-effector interaction. We propose these compounds as starting points for the development of non-covalent allosteric KRAS inhibitors.

Membrane anchoring of the DIRAS3 N-terminal extension permits tumor suppressor function.

DIRAS3 is an imprinted tumor suppressor gene encoding a GTPase that has a distinctive N-terminal extension (NTE) not found in other RAS proteins. This NTE and the prenylated C-terminus are required for DIRAS3-mediated inhibition of RAS/MAP signaling and PI3K activity at the plasma membrane. In this study, we applied biochemical, biophysical, and computational methods to characterize the structure and function of the NTE. The NTE peptide recognizes phosphoinositides PI(3,4,5)P3 and PI(4,5)P2 with rapid kinetics and strong affinity. Lipid binding induces NTE structural change from disorder to amphipathic helix. Mass spectrometry identified N-myristoylation of DIRAS3. All-atom molecular dynamic simulations predict DIRAS3 could adhere to the membrane through both termini, suggesting the NTE is involved in targeting and stabilizing DIRAS3 on the membrane by double anchoring. Overall, our results are consistent with DIRAS3's function as a tumor suppressor, whereby the membrane-bound DIRAS3 can effectively target PI3K and KRAS at the membrane.

A new ferrocene derivative blocks K-Ras localization and function by oxidative modification at His95.

Ras proteins are membrane-bound GTPases that regulate essential cellular processes at the plasma membrane (PM). Constitutively active mutations of K-Ras, one of the three Ras isoforms in mammalian cells, are frequently found in human cancers. Ferrocene derivatives, which elevate cellular reactive oxygen species (ROS), have shown to block the growth of non-small cell lung cancers harboring oncogenic mutant K-Ras. Here, we tested a novel ferrocene derivative on the growth of pancreatic ductal adenocarcinoma and non-small cell lung cancer. Our compound, which elevated cellular ROS levels, inhibited the growth of K-Ras-driven cancers, and abrogated the PM binding and signaling of K-Ras in an isoform-specific manner. These effects were reversed upon antioxidant supplementation, suggesting a ROS-mediated mechanism. We further identified that K-Ras His95 residue plays an important role in this process, and it is putatively oxidized by cellular ROS. Together, our study demonstrates that the redox system directly regulates K-Ras/PM binding and signaling via oxidative modification at the His95, and proposes a role of oncogenic mutant K-Ras in the recently described antioxidant-induced growth and metastasis of K-Ras-driven cancers.

Small-Molecule Inhibition of KRAS through Conformational Selection.

Mutations in KRAS account for about 20% of human cancers. Despite the major progress in recent years toward the development of KRAS inhibitors, including the discovery of covalent inhibitors of the G12C KRAS variant for the treatment of non-small-cell lung cancer, much work remains to be done to discover broad-acting inhibitors to treat many other KRAS-driven cancers. In a previous report, we showed that a 308.4 Da small-molecule ligand [(2R)-2-(N'-(1H-indole-3-carbonyl)hydrazino)-2-phenyl-acetamide] binds to KRAS with low micro-molar affinity [Chem. Biol. Drug Des.2019; 94(2):1441-1456]. Binding of this ligand, which we call ACA22, to the p1 pocket of KRAS and its interactions with residues at beta-strand 1 and the switch loops have been supported by data from nuclear magnetic resonance spectroscopy and microscale thermophoresis experiments. However, the inhibitory potential of the compound was not demonstrated. Here, we show that ACA22 inhibits KRAS-mediated signal transduction in cells expressing wild type (WT) and G12D mutant KRAS and reduces levels of guanosine triphosphate-loaded WT KRAS more effectively than G12D KRAS. We ruled out the direct effect on nucleotide exchange or effector binding as possible mechanisms of inhibition using a variety of biophysical assays. Combining these observations with binding data that show comparable affinities of the compound for the active and inactive forms of the mutant but not the WT, we propose conformational selection as a possible mechanism of action of ACA22.

A new ferrocene derivative blocks KRAS localization and function by oxidative modification at His95.

Ras proteins are membrane-bound GTPases that regulate essential cellular processes at the plasma membrane (PM). Constitutively active mutations of K-Ras, one of the three Ras isoforms in mammalian cells, are frequently found in human cancers. Ferrocene derivatives, which elevate cellular reactive oxygen species (ROS), have shown to block the growth of non-small cell lung cancers (NSCLCs) harboring oncogenic mutant K-Ras. Here, we developed and tested a novel ferrocene derivative on the growth of human pancreatic ductal adenocarcinoma (PDAC) and NSCLC. Our compound inhibited the growth of K-Ras-dependent PDAC and NSCLC and abrogated the PM binding and signaling of K-Ras, but not other Ras isoforms. These effects were reversed upon antioxidant supplementation, suggesting a ROS-mediated mechanism. We further identified K-Ras His95 residue in the G-domain as being involved in the ferrocene-induced K-Ras PM dissociation via oxidative modification. Together, our studies demonstrate that the redox system directly regulates K-Ras PM binding and signaling via oxidative modification at the His95, and proposes a role of oncogenic mutant K-Ras in the recently described antioxidant-induced metastasis in K-Ras-driven lung cancers.

Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor.

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 µs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.

KRAS Inhibitor that Simultaneously Inhibits Nucleotide Exchange Activity and Effector Engagement.

We describe a small molecule ligand ACA-14 (2-hydroxy-5-{[(2-phenylcyclopropyl) carbonyl] amino} benzoic acid) as an initial lead for the development of direct inhibitors of KRAS, a notoriously difficult anticancer drug target. We show that the compound binds to KRAS near the switch regions with affinities in the low micromolar range and exerts different effects on KRAS interactions with binding partners. Specifically, ACA-14 impedes the interaction of KRAS with its effector Raf and reduces both intrinsic and SOS-mediated nucleotide exchange rates. Likely as a result of these effects, ACA-14 inhibits signal transduction through the MAPK pathway in cells expressing mutant KRAS and inhibits the growth of pancreatic and colon cancer cells harboring mutant KRAS. We thus propose compound ACA-14 as a useful initial lead for the development of broad-acting inhibitors that target multiple KRAS mutants and simultaneously deplete the fraction of GTP-loaded KRAS while abrogating the effector-binding ability of the already GTP-loaded fraction.

RAS Nanoclusters Selectively Sort Distinct Lipid Headgroups and Acyl Chains.

RAS proteins are lipid-anchored small GTPases that switch between the GTP-bound active and GDP-bound inactive states. RAS isoforms, including HRAS, NRAS and splice variants KRAS4A and KRAS4B, are some of the most frequently mutated proteins in cancer. In particular, constitutively active mutants of KRAS comprise ∼80% of all RAS oncogenic mutations and are found in 98% of pancreatic, 45% of colorectal and 31% of lung tumors. Plasma membrane (PM) is the primary location of RAS signaling in biology and pathology. Thus, a better understanding of how RAS proteins localize to and distribute on the PM is critical to better comprehend RAS biology and to develop new strategies to treat RAS pathology. In this review, we discuss recent findings on how RAS proteins sort lipids as they undergo macromolecular assembly on the PM. We also discuss how RAS/lipid nanoclusters serve as signaling platforms for the efficient recruitment of effectors and signal transduction, and how perturbing the PM biophysical properties affect the spatial distribution of RAS isoforms and their functions.

Antipsychotic phenothiazine drugs bind to KRAS in vitro.

We used NMR to show that the antipsychotic phenothiazine drugs promazine and promethazine bind to GDP-KRAS. Promazine also binds to oncogenic GDP-KRAS(G12D), and to wild type GppNHp-KRAS. A panel of additional phenothiazines bind to GDP-KRAS but with lower affinity than promazine or promethazine. Binding is most dependent on substitutions at C-2 of the tricyclic phenothiazine ring. Promazine was used to generate an NMR-driven HADDOCK model of the drug/GDP-KRAS complex. The structural model shows the tricyclic phenothiazine ring of promazine associates with the hydrophobic pocket p1 that is bordered by the central ß sheet and Switch II in KRAS. Binding appears to stabilize helix 2 in a conformation that is similar to that seen in KRAS bound to other small molecules. Association of phenothiazines with KRAS may affect normal KRAS signaling that could contribute to multiple biological activities of these antipsychotic drugs. Moreover, the phenothiazine ring represents a new core scaffold on which to design modulators of KRAS activity.

Monoubiquitination of KRAS at Lysine104 and Lysine147 Modulates Its Dynamics and Interaction with Partner Proteins.

KRAS, a 21 kDa guanine nucleotide-binding protein that functions as a molecular switch, plays a key role in regulating cellular growth. Dysregulation of this key signaling node leads to uncontrolled cell growth, a hallmark of cancer cells. KRAS undergoes post-translational modification by monoubiquitination at various locations, including at lysine104 (K104) and lysine147 (K147). Previous studies have suggested that K104 stabilizes helix-2/helix-3 interactions and K147 is involved in nucleotide binding. However, the impact of monoubiquitination at these residues on the overall structure, dynamics, or function of KRAS is not fully understood. In this study, we examined KRAS monoubiquitination at these sites using data from extensive (12 µs aggregate time) molecular dynamics simulations complemented by nuclear magnetic resonance spectroscopy data. We found that ubiquitin forms dynamic nonspecific interactions with various regions of KRAS and that ubiquitination at both sites modulates conformational fluctuations. In both cases, ubiquitin samples a broad range of conformational space and does not form long-lasting noncovalent contacts with KRAS but it adopts several preferred orientations relative to KRAS. To examine the functional impact of these preferred orientations, we performed a systematic comparison of the dominant configurations of the ubiquitin/KRAS simulated complex with experimental structures of KRAS bound to regulatory and effector proteins as well as a model membrane. Results from these analyses suggest that conformational selection and population shift may minimize the deleterious effects of KRAS ubiquitination at K104 and K147 on binding to some but not all interaction partners. Our findings thus provide new insights into the steric effects of ubiquitin and suggest a potential avenue for therapeutic targeting.

Approaches to inhibiting oncogenic K-Ras.

Activating somatic K-Ras mutations are associated with >15% all human tumors and up to 90% of specific tumor types such as pancreatic cancer. Successfully inhibiting abnormal K-Ras signaling would therefore be a game changer in cancer therapy. However, K-Ras has long been considered an undruggable target for various reasons. This view is now changing by the discovery of allosteric inhibitors that directly target K-Ras and inhibit its functions, and by the identification of new mechanisms to dislodge it from the plasma membrane and thereby abrogate its cellular activities. In this review, we will discuss recent progresses and challenges to inhibiting aberrant K-Ras functions by these two approaches. We will also provide a broad overview of other approaches such as inhibition of K-Ras effectors, and offer a brief perspective on the way forward.

Mechanisms of Ras Membrane Organization and Signaling: Ras Rocks Again.

Ras is the most frequently mutated oncogene and recent drug development efforts have spurred significant new research interest. Here we review progress toward understanding how Ras functions in nanoscale, proteo-lipid signaling complexes on the plasma membrane, called nanoclusters. We discuss how G-domain reorientation is plausibly linked to Ras-nanoclustering and -dimerization. We then look at how these mechanistic features could cooperate in the engagement and activation of RAF by Ras. Moreover, we show how this structural information can be integrated with microscopy data that provide nanoscale resolution in cell biological experiments. Synthesizing the available data, we propose to distinguish between two types of Ras nanoclusters, an active, immobile RAF-dependent type and an inactive/neutral membrane anchor-dependent. We conclude that it is possible that Ras reorientation enables dynamic Ras dimerization while the whole Ras/RAF complex transits into an active state. These transient di/oligomer interfaces of Ras may be amenable to pharmacological intervention. We close by highlighting a number of open questions including whether all effectors form active nanoclusters and whether there is an isoform specific composition of Ras nanocluster.

How to make an undruggable enzyme druggable: lessons from ras proteins.

Significant advances have been made toward discovering allosteric inhibitors for challenging drug targets such as the Ras family of membrane-associated signaling proteins. Malfunction of Ras proteins due to somatic mutations is associated with up to a quarter of all human cancers. Computational techniques have played critical roles in identifying and characterizing allosteric ligand-binding sites on these proteins, and to screen ligand libraries against those sites. These efforts, combined with a wide range of biophysical, structural, biochemical and cell biological experiments, are beginning to yield promising inhibitors to treat malignancies associated with mutated Ras proteins. In this chapter, we discuss some of these developments and how the lessons learned from Ras might be applied to similar other challenging drug targets.

Transmembrane potential of physiologically relevant model membranes: Effects of membrane asymmetry.

Transmembrane potential difference (Vm) plays important roles in regulating various biological processes. At the macro level, Vm can be experimentally measured or calculated using the Nernst or Goldman-Hodgkin-Katz equation. However, the atomic details responsible for its generation and impact on protein and lipid dynamics still need to be further elucidated. In this work, we performed a series of all-atom molecular dynamics (MD) simulations of symmetric model membranes of various lipid compositions and cation contents to evaluate the relationship between membrane asymmetry and Vm. Specifically, we studied the impact of the asymmetric distribution of POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine), PIP2 (phosphatidylinositol 4,5-bisphosphate), as well as Na+ and K+ on Vm using atomically detailed MD simulations of symmetric model membranes. The results suggest that, for an asymmetric POPC-POPC/POPS bilayer in the presence of NaCl, the presence of the monovalent anionic lipid POPS in the inner leaflet polarizes the membrane (ΔVm < 0). Intriguingly, replacing a third of the POPS lipids by the polyvalent anionic signaling lipid PIP2 counteracts this effect, resulting in a smaller negative membrane potential. We also found that replacing Na+ ions in the inner region by K+ depolarizes the membrane (ΔVm > 0). These divergent effects arise from variations in the strength of cation-lipid interactions and are correlated with changes in lipid chain order and head-group orientation.

Small Molecule KRAS Inhibitors: The Future for Targeted Pancreatic Cancer Therapy?

Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest solid tumors in the world. Currently, there are no approved targeted therapies for PDAC. Mutations in Kirsten rat sarcoma viral oncogene homologue (KRAS) are known to be a major driver of PDAC progression, but it was considered an undruggable target until recently. Moreover, PDAC also suffers from drug delivery issues due to the highly fibrotic tumor microenvironment. In this perspective, we provide an overview of recent developments in targeting mutant KRAS and strategies to overcome drug delivery issues (e.g., nanoparticle delivery). Overall, we propose that the antitumor effects from novel KRAS inhibitors along with strategies to overcome drug delivery issues could be a new therapeutic way forward in PDAC.

Probing the Conformational and Energy Landscapes of KRAS Membrane Orientation.

Membrane reorientation of oncogenic RAS proteins is emerging as an important modulator of their functions. Previous studies have shown that the most common orientations include those with either the three C-terminal α-helices (OS1) or N-terminal ß-strands (OS2) of the catalytic domain facing the membrane. OS1 and OS2 differ by the degree to which the effector-interacting surface is occluded by the membrane. However, the relative stability of these states and the rates of transition between them remained undetermined. How mutations might modulate preferences for specific orientation states is also far from clear. The current work attempted to address these questions through a comprehensive analysis of two 20 µs-long atomistic molecular dynamics simulations. The simulations were conducted on the oncogenic G12D and Q61H KRAS mutants bound to an anionic lipid bilayer. G12D and Q61H are among the most prevalent cancer-causing mutations at the P-loop and switch 2 regions of KRAS, respectively. We found that both mutants fluctuate in a similar manner between OS1 and OS2 via an intermediate orientation OS0, and both favor the signaling competent OS1 and OS0 over the occluded OS2. However, they differ in the details, such as in the extent to which they sample OS1. Analysis of the orientation free-energy landscapes estimated from the simulations indicate that OS1 and OS2 are the most stable states. However, the overall free energy surface is rugged, indicating a large diversity of conformations including at least two substates in each orientation state that differ in stability only by about 0.5-1.0 kcal/mol. Reversible transitions between OS1 and OS2 occur via two well-defined pathways that traverse OS0. In the minimum energy path, helix 4 remains close to the membrane as the angle of the catalytic domain from the membrane plane changes, resulting in a barrier of ∼1 kcal/mol for OS1/OS2 interconversions. Estimation of the rates of the various transitions based on survival probabilities yielded two rate constants in the order of 107 and 106 s-1, which we attribute to intrinsic protein conformational dynamics and transient protein-lipid interactions, respectively. The faster process dominates every transition, confirming a previous suggestion that RAS membrane reorientation is driven by conformational fluctuations rather than protein-lipid interactions.

Membrane curvature sensing of the lipid-anchored K-Ras small GTPase.

Plasma membrane (PM) curvature defines cell shape and intracellular organelle morphologies and is a fundamental cell property. Growth/proliferation is more stimulated in flatter cells than the same cells in elongated shapes. PM-anchored K-Ras small GTPase regulates cell growth/proliferation and plays key roles in cancer. The lipid-anchored K-Ras form nanoclusters selectively enriched with specific phospholipids, such as phosphatidylserine (PS), for efficient effector recruitment and activation. K-Ras function may, thus, be sensitive to changing lipid distribution at membranes with different curvatures. Here, we used complementary methods to manipulate membrane curvature of intact/live cells, native PM blebs, and synthetic liposomes. We show that the spatiotemporal organization and signaling of an oncogenic mutant K-Ras G12V favor flatter membranes with low curvature. Our findings are consistent with the more stimulated growth/proliferation in flatter cells. Depletion of endogenous PS abolishes K-Ras G12V PM curvature sensing. In cells and synthetic bilayers, only mixed-chain PS species, but not other PS species tested, mediate K-Ras G12V membrane curvature sensing. Thus, K-Ras nanoclusters act as relay stations to convert mechanical perturbations to mitogenic signaling.

Discovery of High-Affinity Noncovalent Allosteric KRAS Inhibitors That Disrupt Effector Binding.

Approximately 15% of all human tumors harbor mutant KRAS, a membrane-associated small GTPase and notorious oncogene. Mutations that render KRAS constitutively active will lead to uncontrolled cell growth and cancer. However, despite aggressive efforts in recent years, there are no drugs on the market that directly target KRAS and inhibit its aberrant functions. In the current work, we combined structure-based design with a battery of cell and biophysical assays to discover a novel pyrazolopyrimidine-based allosteric KRAS inhibitor that binds to activated KRAS with sub-micromolar affinity and disrupts effector binding, thereby inhibiting KRAS signaling and cancer cell growth. These results show that pyrazolopyrimidine-based compounds may represent a first-in-class allosteric noncovalent inhibitors of KRAS. Moreover, by studying two of its analogues, we identified key chemical features of the compound that interact with a set of specific residues at the switch regions of KRAS and play critical roles for its high-affinity binding and unique mode of action, thus providing a blueprint for future optimization efforts.