Novel allosteric sites on Ras for lead generation

Exploring Binding Sites in Chagas Disease Protein TcP21 Using Integrated Mixed Solvent Molecular Dynamics Approaches

Chagas disease, caused by the protozoan Trypanosoma cruzi, remains a significant global health burden, particularly in Latin America, where millions are at risk. This disease predominantly affects socioeconomically vulnerable populations, aggravating economic inequality, marginalization, and low political visibility. Despite extensive research, effective treatments are still lacking, partly due to the complex biology of the parasite and its infection mechanisms. This study focuses on TcP21, a novel 21 kDa protein secreted by extracellular amastigotes, which has been implicated in T. cruzi infection via an alternative infective pathway. Although the potential of TcP21 for understanding Chagas disease is promising, further exploration is necessary, particularly in identifying potential binding sites on its surface. Computational tools offer a versatile and effective strategy for preliminary binding site assessment, facilitating a more cost-efficient allocation of experimental resources. In this study, we employed three independent computational approaches─mixed solvent molecular dynamics simulations (MSMD), fragment-based molecular docking, and pharmacophore model docking coupled with molecular dynamics simulations─to identify potential binding sites and provide comprehensive insights into TcP21. The three methodologies converged on a common site located on the external surface of the protein, characterized by key residues such as GLU55, ASP52, VAL70, ILE62, and TRP77. The protonated amino, acetamido, and phenyl groups of the pharmacophore probe were consistently observed to interact with the site via a network of salt bridges, hydrogen bonds, charge-charge interactions, and alkyl-π interactions, suggesting these groups play a significant role in ligand binding. This study does not aim to propose specific therapeutic hits but to highlight a still unknown and unexplored protein involved in T. cruzi cell invasion. In this regard, given the strong correlation between the three distinct approaches used for mapping, we consider this study offers valuable insights for further research into P21 and its role in Chagas disease.

Inhibition of GTPase KRASG12D: a review of patent literature

INTRODUCTION

KRAS is a critical oncogenic protein intricately involved in tumor progression, and the difficulty in targeting KRAS has led it to be classified as an 'undruggable target.' Among the various KRAS mutations, KRASG12D is highly prevalent and represents a promising therapeutic target, yet there are currently no approved inhibitors for it.

AREA COVERED

This review summarizes numerous patents and literature featuring inhibitors or degraders of KRASG12D through searching relevant information in PubMed, SciFinder and Web of Science databases from 2021 to February 2024, providing an overview of the research progress on inhibiting KRASG12D in terms of design strategies, chemical structures, biological activities, and clinical advancements.

EXPERT OPINION

Since the approval of AMG510 (Sotorasib), there has been an increasing focus on the inhibition of KRASG12D, leading to numerous reports of related inhibitors and degraders. Among them, MRTX1133, as the first KRASG12D inhibitor to enter clinical trials, has demonstrated excellent tumor suppression in various KRASG12D-bearing human tumor xenograft models. It is important to note, however, that understanding the mechanisms of acquired resistance caused by KRAS inhibition and developing additional combination therapies is crucial. Moreover, seeking covalent inhibition of KRASG12D also holds significant potential.

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.

Progress in RAS-targeted therapeutic strategies: From small molecule inhibitors to proteolysis targeting chimeras

As a widely considerable target in chemical biology and pharmacological research, rat sarcoma (RAS) gene mutations play a critical driving factor in several fatal cancers. Despite the great progress of RAS subtype-specific inhibitors, rapid acquired drug resistance could limit their further clinical applications. Proteolysis targeting chimera (PROTAC) has emerged as a powerful tool to handle "undruggable" targets and exhibited significant therapeutic benefit for the combat of drug resistance. Owing to unique molecular mechanism and binding kinetics, PROTAC is expected to become a feasible strategy to break the bottleneck of classical RAS inhibitors. This review aims to discuss the current advances of RAS inhibitors and especially focus on PROTAC strategy targeting RAS mutations and their downstream effectors for relevant cancer treatment.

The energetic and allosteric landscape for KRAS inhibition

Thousands of proteins have been validated genetically as therapeutic targets for human diseases1. However, very few have been successfully targeted, and many are considered 'undruggable'. This is particularly true for proteins that function via protein-protein interactions-direct inhibition of binding interfaces is difficult and requires the identification of allosteric sites. However, most proteins have no known allosteric sites, and a comprehensive allosteric map does not exist for any protein. Here we address this shortcoming by charting multiple global atlases of inhibitory allosteric communication in KRAS. We quantified the effects of more than 26,000 mutations on the folding of KRAS and its binding to six interaction partners. Genetic interactions in double mutants enabled us to perform biophysical measurements at scale, inferring more than 22,000 causal free energy changes. These energy landscapes quantify how mutations tune the binding specificity of a signalling protein and map the inhibitory allosteric sites for an important therapeutic target. Allosteric propagation is particularly effective across the central β-sheet of KRAS, and multiple surface pockets are genetically validated as allosterically active, including a distal pocket in the C-terminal lobe of the protein. Allosteric mutations typically inhibit binding to all tested effectors, but they can also change the binding specificity, revealing the regulatory, evolutionary and therapeutic potential to tune pathway activation. Using the approach described here, it should be possible to rapidly and comprehensively identify allosteric target sites in many proteins.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Automation of Drug Discovery through Cutting-edge In-silico Research in Pharmaceuticals: Challenges and Future Scope

The rapidity and high-throughput nature of in silico technologies make them advantageous for predicting the properties of a large array of substances. In silico approaches can be used for compounds intended for synthesis at the beginning of drug development when there is either no or very little compound available. In silico approaches can be used for impurities or degradation products. Quantifying drugs and related substances (RS) with pharmaceutical drug analysis (PDA) can also improve drug discovery (DD) by providing additional avenues to pursue. Potential future applications of PDA include combining it with other methods to make insilico predictions about drugs and RS. One possible outcome of this is a determination of the drug potential of nontoxic RS. ADME estimation, QSAR research, molecular docking, bioactivity prediction, and toxicity testing all involve impurity profiling. Before committing to DD, RS with minimal toxicity can be utilised in silico. The efficacy of molecular docking in getting a medication to market is still debated despite its refinement and improvement. Biomedical labs and pharmaceutical companies were hesitant to adopt molecular docking algorithms for drug screening despite their decades of development and improvement. Despite the widespread use of "force fields" to represent the energy exerted within and between molecules, it has been impossible to reliably predict or compute the binding affinities between proteins and potential binding medications.

Discovery of Hit Compounds Targeting the P4 Allosteric Site of K-RAS, Identified through Ensemble-Based Virtual Screening

Mutants of Ras are oncogenic drivers of a large number of human tumors. Despite being recognized as an attractive target for the treatment of cancer, the high affinity for its substrate tagged the protein as undruggable for a few years. The identification of cryptic pockets on the protein surface gave the opportunity to identify molecules capable of acting as allosteric modulators. Several molecules were disclosed in recent years, with sotorasib and adagrasib already approved for clinical use. The present study makes use of computational methods to characterize eight prospective allosteric pockets (P1-P8) in K-Ras, four of which (P1-P4) were previously characterized in the literature. The present study also describes the results of a virtual screening study focused on the discovery of hit compounds, binders of the P4 site that can be considered as peptidomimetics of a fragment of the SOS αI helix, a guanine exchange factor of Ras. After a detailed description of the computational procedure followed, we disclose five hit compounds, prospective binders of the P4 allosteric site that exhibit an inhibitory capability higher than 30% in a cell proliferation assay at 50 μM.

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.

Reduced dynamic complexity allows structure elucidation of an excited state of KRASG13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS^G13D, which suggests an exchange between two conformational states in solution. Methyl and ³¹P NMR spectra of active KRAS^G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major P^γ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS^G13D and KRAS^G13D-RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS^G13D, showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

Eliminating oncogenic RAS: back to the future at the drawing board

RAS drug development has made enormous strides in the past ten years, with the first direct KRAS inhibitor being approved in 2021. However, despite the clinical success of covalent KRAS-G12C inhibitors, we are immediately confronted with resistances as commonly found with targeted drugs. Previously believed to be undruggable due to its lack of obvious druggable pockets, a couple of new approaches to hit this much feared oncogene have now been carved out. We here concisely review these approaches to directly target four druggable sites of RAS from various angles. Our analysis focuses on the lessons learnt during the development of allele-specific covalent and non-covalent RAS inhibitors, the potential of macromolecular binders to facilitate the discovery and validation of targetable sites on RAS and finally an outlook on a future that may engage more small molecule binders to become drugs. We foresee that the latter could happen mainly in two ways: First, non-covalent small molecule inhibitors may be derived from the development of covalent binders. Second, reversible small molecule binders could be utilized for novel targeting modalities, such as degraders of RAS. Provided that degraders eliminate RAS by recruiting differentially expressed E3-ligases, this approach could enable unprecedented tissue- or developmental stage-specific destruction of RAS with potential advantages for on-target toxicity. We conclude that novel creative ideas continue to be important to exterminate RAS in cancer and other RAS pathway-driven diseases, such as RASopathies.

Exploring the state- and allele-specific conformational landscapes of Ras: understanding their respective druggabilities

Recent advances in direct inhibition of Ras benefit from the protein's intrinsic dynamic nature that derives therapeutically vulnerable conformers bearing transiently formed cryptic pockets. Hotspot mutants of Ras are major tumor drivers and are hyperactivated in cells at variable levels, which may require allele-specific strategies for effective targeting. However, it remains unclear how the prevalent oncogenic mutations and activation states perturb the free energy landscape governing the protein dynamics and druggability. Here we characterized the nucleotide state- and allele-dependent alterations of Ras conformational dynamics using a combined NMR experimental and computational approach and constructed quantitative ensembles revealing the conservation of the cryptic SI/II-P and SII-P pockets in different states and alleles. Highly local but critical conformational reorganizations that undermine the SII-P accessibility to residue 12 have been identified as a common mechanism resulting in the low reactivities of Ras·GTP as well as Ras(G12D)·GDP with covalent SII-P inhibitors. Our results strongly support the conformational selection scenario for interactions between Ras and the previously reported binders and offer insights for the future development of state- and allele-specific, as well as pan-Ras, inhibitors.

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

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

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.

An evaluation of sotorasib for the treatment of patients with non-small cell lung cancer with KRAS G12C mutations

INTRODUCTION

Improving the clinical outcomes of patients with KRAS^G12C-mutated non-small cell lung cancer (NSCLC), the majority of whom are current or former smokers, has been a barrier to improving population-level outcomes in NSCLC. Novel, effective KRAS^G12C inhibitors are emerging and sotorasib is the first member of that class to achieve commercial availability.

AREAS COVERED

In this review, we survey the epidemiology of KRAS^G12C-mutated NSCLC, as well as sotorasib's chemistry, pharmacology, and clinical trial data.

EXPERT OPINION

While sotorasib's development has been unique and exciting, questions persist regarding its intracranial penetrance, optimal dose, and efficacy relative to standard-of-care therapy. Improvements in the clinical activity of KRAS inhibition will hinge on better understanding of resistance mechanisms, the development of broad-spectrum inhibitors with activity beyond G12C mutations, and combination therapy targeting multiple mediators of KRAS signaling and alternative pathways. From a regulatory perspective, sotorasib's development may, in time, prove to be an instructive example for early-phase clinical trialists and regulators focused on dose optimization.

Computer-Aided Drug Design Boosts RAS Inhibitor Discovery

The Rat Sarcoma (RAS) family (NRAS, HRAS, and KRAS) is endowed with GTPase activity to regulate various signaling pathways in ubiquitous animal cells. As proto-oncogenes, RAS mutations can maintain activation, leading to the growth and proliferation of abnormal cells and the development of a variety of human cancers. For the fight against tumors, the discovery of RAS-targeted drugs is of high significance. On the one hand, the structural properties of the RAS protein make it difficult to find inhibitors specifically targeted to it. On the other hand, targeting other molecules in the RAS signaling pathway often leads to severe tissue toxicities due to the lack of disease specificity. However, computer-aided drug design (CADD) can help solve the above problems. As an interdisciplinary approach that combines computational biology with medicinal chemistry, CADD has brought a variety of advances and numerous benefits to drug design, such as the rapid identification of new targets and discovery of new drugs. Based on an overview of RAS features and the history of inhibitor discovery, this review provides insight into the application of mainstream CADD methods to RAS drug design.

Modeling receptor flexibility in the structure-based design of KRASG12C inhibitors

KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS^G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS^G12C—sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS^G12C. In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS^G12C, we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS^G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.

Therapeutic Targeting of Rab GTPases: Relevance for Alzheimer's Disease

Rab GTPases (Rabs) are small proteins that play crucial roles in vesicle transport and membrane trafficking. Owing to their widespread functions in several steps of vesicle trafficking, Rabs have been implicated in the pathogenesis of several disorders, including cancer, diabetes, and multiple neurodegenerative diseases. As treatments for neurodegenerative conditions are currently rather limited, the identification and validation of novel therapeutic targets, such as Rabs, is of great importance. This review summarises proof-of-concept studies, demonstrating that modulation of Rab GTPases in the context of Alzheimer's disease (AD) can ameliorate disease-related phenotypes, and provides an overview of the current state of the art for the pharmacological targeting of Rabs. Finally, we also discuss the barriers and challenges of therapeutically targeting these small proteins in humans, especially in the context of AD.

Molecular dynamics simulations reveal the selectivity mechanism of structurally similar agonists to TLR7 and TLR8

TLR7 and TLR8 are key members of the Toll-like receptor family, playing crucial roles in the signaling pathways of innate immunity, and thus become attractive therapeutic targets of many diseases including infections and cancer. Although TLR7 and TLR8 show a high degree of sequence homology, their biological response to small molecule binding is very different. Aiming to understand the mechanism of selective profiles of small molecule modulators against TLR7 and TLR8, we carried out molecular dynamic simulations on three imidazoquinoline derivatives bound to the receptors separately. They are Resiquimod (R), Hybrid-2 (H), and Gardiquimod (G), selective agonists of TLR7 and TLR8. Our MD trajectories indicated that in the complex of TLR7-R and TLR7-G, the two chains forming the TLR7 dimer tended to remain “open” conformation, while the rest systems maintained in the closed format. The agonists R, H, and G developed conformational deviation mainly on the aliphatic tail. Furthermore, we attempted to quantify the selectivity between TLR7 and TLR8 by binding free energies via MM-GBSA method. It showed that the three selected modulators were more favorable for TLR7 than TLR8, and the ranking from the strongest to the weakest was H, R and G, aligning well with experimental data. In the TLR7, the flexible and hydrophobic aliphatic side chain of H has stronger van der Waals interactions with V381 and F351 but only pick up interaction with one amino acid residue i.e. Y353 of TLR8. Unsurprisingly, the positively charged side chain of G has less favorable interaction with I585 of TLR7 and V573 of TLR8 explaining G is weak agonist of both TLR7 and TLR8. All three imidazoquinoline derivatives can form stable hydrogen bonds with D555 of TLR7 and the corresponding D543 of TLR8. In brief, the set of total 400ns MD studies sheds light on the potential selectivity mechanisms of agonists towards TLR7 and TLR8, indicating the van der Waals interaction as the driving force for the agonists binding, thus provides us insights for designing more potent and selective modulators to cooperate with the hydrophobic nature of the binding pocket.

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.

In silico methods and tools for drug discovery

In the past, conventional drug discovery strategies have been successfully employed to develop new drugs, but the process from lead identification to clinical trials takes more than 12 years and costs approximately $1.8 billion USD on average. Recently, in silico approaches have been attracting considerable interest because of their potential to accelerate drug discovery in terms of time, labor, and costs. Many new drug compounds have been successfully developed using computational methods. In this review, we briefly introduce computational drug discovery strategies and outline up-to-date tools to perform the strategies as well as available knowledge bases for those who develop their own computational models. Finally, we introduce successful examples of anti-bacterial, anti-viral, and anti-cancer drug discoveries that were made using computational methods.

Computational modeling reveals key molecular properties and dynamic behavior of disruptor of telomeric silencing 1-like (DOT1L) and partnering complexes involved in leukemogenesis

Disruptor of telomeric silencing 1‐like (DOT1L) is the only non‐SET domain histone lysine methyltransferase (KMT) and writer of H3K79 methylation on nucleosomes marked by H2B ubiquitination. DOT1L has elicited significant attention because of its interaction or fusion with members of the AF protein family in blood cell biology and leukemogenic transformation. Here, our goal was to extend previous structural information by performing a robust molecular dynamic study of DOT1L and its leukemogenic partners combined with mutational analysis. We show that statically and dynamically, D161, G163, E186, and F223 make frequent time‐dependent interactions with SAM, while additional residues T139, K187, and N241 interact with SAM only under dynamics. Dynamics models reveal DOT1L, SAM, and H4 moving as one and show that more than twice the number of DOT1L residues interacts with these partners, relative to the static structure. Mutational analyses indicate that six of these residues are intolerant to substitution. We describe the dynamic behavior of DOT1L interacting with AF10 and AF9. Studies on the dynamics of a heterotrimeric complex of DOT1L1‐AF10 illuminated describe coordinated motions that impact the relative position of the DOT1L HMT domain to the nucleosome. The molecular motions of the DOT1L–AF9 complex are less extensive and highly dynamic, resembling a swivel‐like mechanics. Through molecular dynamics and mutational analysis, we extend the knowledge previous provided by static measurements. These results are important to consider when describing the biochemical properties of DOT1L, under normal and in disease conditions, as well as for the development of novel therapeutic agents.

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.

RAS-inhibiting biologics identify and probe druggable pockets including an SII-α3 allosteric site

RAS mutations are the most common oncogenic drivers across human cancers, but there remains a paucity of clinically-validated pharmacological inhibitors of RAS, as druggable pockets have proven difficult to identify. Here, we identify two RAS-binding Affimer proteins, K3 and K6, that inhibit nucleotide exchange and downstream signaling pathways with distinct isoform and mutant profiles. Affimer K6 binds in the SI/SII pocket, whilst Affimer K3 is a non-covalent inhibitor of the SII region that reveals a conformer of wild-type RAS with a large, druggable SII/α3 pocket. Competitive NanoBRET between the RAS-binding Affimers and known RAS binding small-molecules demonstrates the potential to use Affimers as tools to identify pharmacophores. This work highlights the potential of using biologics with small interface surfaces to select unseen, druggable conformations in conjunction with pharmacophore identification for hard-to-drug proteins.

Oncogenic RAS mutants remain difficult to target with small molecules. Here, the authors show that RAS-binding Affimer proteins inhibit RAS signaling while binding diverse regions on the RAS surface, suggesting the potential to use Affimers as tools to identify new binding pockets and pharmacophores.

Normal Mode Analysis of KRas4B Reveals Partner Specific Dynamics

Ras GTPase interacts with its regulators and downstream effectors for its critical function in cellular signaling. Targeting the disrupted mechanisms in Ras-related human cancers requires understanding the distinct dynamics of these protein-protein interactions. We performed normal mode analysis (NMA) of KRas4B in wild-type or mutant monomeric and neurofibromin-1 (NF1), Son of Sevenless 1 (SOS1) or Raf-1 bound dimeric conformational states to reveal partner-specific dynamics of the protein. Gaussian network model (GNM) analysis showed that the known KRas4B lobes further partition into subdomains upon binding to its partners. Furthermore, KRas4B interactions with different partners suppress the flexibility in not only their binding sites but also distant residues in the allosteric lobe in a partner-specific way. The conformational changes can be driven by intrinsic residue fluctuations of the open state KRas4B-GDP, as we illustrated with anisotropic network model (ANM) analysis. The allosteric paths connecting the nucleotide binding residues to the allosteric site at α3-L7 portray differences in the inactive and active states. These findings help in understanding the partner-specific KRas4B dynamics, which could be utilized for therapeutic targeting.

ATP-Magnesium Coordination: Protein Structure-Based Force Field Evaluation and Corrections

In the numerous molecular recognition and catalytic processes across biochemistry involving adenosine triphosphate (ATP), the common bioactive form is its magnesium chelate, ATP·Mg2+. In aqueous solution, two chelation geometries predominate, distinguished by bidentate and tridentate Mg2+-phosphate coordination. These are approximately isoenergetic but separated by a high energy barrier. Force field-based atomistic simulation studies of this complex require an accurate representation of its structure and energetics. Here we focused on the energetics of ATP·Mg2+ coordination. Applying an enhanced sampling scheme to circumvent prohibitively slow sampling of transitions between coordination modes, we observed striking contradictions between Amber and CHARMM force field descriptions, most prominently in opposing predictions of the favored coordination mode. Through further configurational free energy calculations, conducted against a diverse set of ATP·Mg2+-protein complex structures to supplement otherwise limited experimental data, we quantified systematic biases for each force field. The force field calculations were strongly predictive of experimentally observed coordination modes, enabling additive corrections to the coordination free energy that deliver close agreement with experiment. We reassessed the applicability of the thus corrected force field descriptions of ATP·Mg2+ for biomolecular simulation and observed that, while the CHARMM parameters display an erroneous preference for overextended triphosphate configurations that will affect many common biomolecular simulation applications involving ATP, the force field energy landscapes broadly agree with experimental measurements of solution geometry and the distribution of ATP·Mg2+ structures found in the Protein Data Bank. Our force field evaluation and correction approach, based on maximizing consistency with the large and heterogeneous collection of structural information encoded in the PDB, should be broadly applicable to many other systems.

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.

Silencing of Oncogenic KRAS by Mutant-Selective Small Interfering RNA

Oncogenic mutations in the KRAS gene are well-established drivers of cancer. While the recently developed KRAS^G12C inhibitors offer a targeted KRAS therapy and have shown success in the clinic, KRAS^G12C represents only 11% of all KRAS mutations. Current therapeutic approaches for all other KRAS mutations are both indirect and nonmutant-selective, largely focusing on inhibition of downstream KRAS effectors such as MAP kinases. Inhibition of KRAS downstream signaling results in a system-wide down-modulation of the respective targets, raising concerns about systemic cell toxicity. Here, we describe a custom short interfering RNA oligonucleotide (EFTX-D1) designed to preferentially bind mRNA of the most commonly occurring KRAS missense mutations in codons 12 and 13. We determined that EFTX-D1 preferentially reduced the mutant KRAS sequence versus wild-type at the levels of both transcription and translation and reversed oncogenic KRAS-induced morphologic and growth transformation. Furthermore, EFTX-D1 significantly impaired the proliferation of several KRAS mutant cancer cell lines in 2-D as well as 3-D assays. Taken together, our data indicate a novel use of RNA interference to target oncogenic KRAS-driven cancers specifically.

Inhibition of RAS: proven and potential vulnerabilities

RAS is a membrane localized small GTPase frequently mutated in human cancer. As such, RAS has been a focal target for developing cancer therapeutics since its discovery nearly four decades ago. However, efforts to directly target RAS have been challenging due to the apparent lack of readily discernable deep pockets for binding small molecule inhibitors leading many to consider RAS as undruggable. An important milestone in direct RAS inhibition was achieved recently with the groundbreaking discovery of covalent inhibitors that target the mutant Cys residue in KRAS(G12C). Surprisingly, these G12C-reactive compounds only target mutant RAS in the GDP-bound state thereby locking it in the inactive conformation and blocking its ability to couple with downstream effector pathways. Building on this success, several groups have developed similar compounds that selectively target KRAS(G12C), with AMG510 and MRTX849 the first to advance to clinical trials. Both have shown early promising results. Though the success with these compounds has reignited the possibility of direct pharmacological inhibition of RAS, these covalent inhibitors are limited to treating KRAS(G12C) tumors which account for <15% of all RAS mutants in human tumors. Thus, there remains an unmet need to identify more broadly efficacious RAS inhibitors. Here, we will discuss the current state of RAS(G12C) inhibitors and the potential for inhibiting additional RAS mutants through targeting RAS dimerization which has emerged as an important step in the allosteric regulation of RAS function.

Comparative MD simulations and advanced analytics based studies on wild-type and hot-spot mutant A59G HRas

The Ras family of proteins is known to play an important role in cellular signal transduction. The oncoprotein Ras is also found to be mutated in ~90% of the pancreatic cancers, of which G12V, G13V, A59G and Q61L are the known hot-spot mutants. These ubiquitous proteins fall in the family of G-proteins, and hence switches between active GTP bound and inactive GDP bound states, which is hindered in most of its oncogenic mutant counterparts. Moreover, Ras being a GTPase has an intrinsic property to hydrolyze GTP to GDP, which is obstructed due to mutations and lends the mutants stuck in constitutively active state leading to oncogenic behavior. In this regard, the present study aims to understand the dynamics involved in the hot-spot mutant A59G-Ras using long 10μs classical MD simulations (5μs for each of the wild-type and mutant systems) and comparing the same with its wild-type counterpart. Advanced analytics using Markov State Model (MSM) based approach has been deployed to comparatively understand the transition path for the wild-type and mutant systems. Roles of crucial residues like Tyr32, Gln61 and Tyr64 have also been established using multivariate PCA analyses. Furthermore, this multivariate PCA analysis also provides crucial features which may be used as reaction coordinates for biased simulations for further studies. The absence of formation of pre-hydrolysis network is also reported for the mutant conformation, using the distance-based analyses (between crucial residues) of the conserved regions. The implications of this study strengthen the hypothesis that the disruption of the pre-hydrolysis network in the mutant A59G ensemble might lead to permanently active oncogenic conformation in the mutant conformers.

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.

In silico and saturation transfer difference NMR approaches to unravel the binding mode of an andrographolide derivative to K-Ras oncoprotein

Background: Andrographolide and its benzylidene derivatives, SRJ09 and SRJ23, potentially bind oncogenic K-Ras to exert anticancer activity. Their molecular interactions with K-Ras oncoproteins that lead to effective biological activity are of major interest. Methods & results: In silico docking and molecular dynamics simulation were performed using Glide and Desmond, respectively; while saturation transfer difference NMR was performed using GDP-bound K-RasG12V. SRJ23 was found to bind strongly and selectively to K-RasG12V, by anchoring to a binding pocket (namely p2) principally via hydrogen bond and hydrophobic interactions. The saturation transfer difference NMR analysis revealed the proximity of protons of functional moieties in SRJ23 to K-RasG12V, suggesting positive binding. Conclusion: SRJ23 binds strongly and interacts stably with K-RasG12V to exhibit its inhibitory activity.

RAS Function in cancer cells: translating membrane biology and biochemistry into new therapeutics

The three human RAS proteins are mutated and constitutively activated in ∼20% of cancers leading to cell growth and proliferation. For the past three decades, many attempts have been made to inhibit these proteins with little success. Recently; however, multiple methods have emerged to inhibit KRAS, the most prevalently mutated isoform. These methods and the underlying biology will be discussed in this review with a special focus on KRAS-plasma membrane interactions.

Biology, pathology, and therapeutic targeting of RAS

RAS was identified as a human oncogene in the early 1980s and subsequently found to be mutated in nearly 30% of all human cancers. More importantly, RAS plays a central role in driving tumor development and maintenance. Despite decades of effort, there remain no FDA approved drugs that directly inhibit RAS. The prevalence of RAS mutations in cancer and the lack of effective anti-RAS therapies stem from RAS' core role in growth factor signaling, unique structural features, and biochemistry. However, recent advances have brought promising new drugs to clinical trials and shone a ray of hope in the field. Here, we will exposit the details of RAS biology that illustrate its key role in cell signaling and shed light on the difficulties in therapeutically targeting RAS. Furthermore, past and current efforts to develop RAS inhibitors will be discussed in depth.

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.

Computationally Empowered Workflow Identifies Novel Covalent Allosteric Binders for KRASG12C

Due to its frequent mutations in multiple lethal cancers, KRAS is one of the most-studied anticancer targets nowadays. Since the discovery of the druggable allosteric binding site containing a G12C mutation, KRASG12C has been the focus of attention in oncology research. We report here a computationally driven approach aimed at identifying novel and selective KRASG12C covalent inhibitors. The workflow involved initial enumeration of virtual molecules tailored for the KRAS allosteric binding site. Tools such as pharmacophore modeling, docking, and free-energy perturbations were deployed to prioritize the compounds with the best profiles. The synthesized naphthyridinone scaffold showed the ability to react with G12C and inhibit KRASG12C . Analogues were prepared to establish structure-activity relationships, while molecular dynamics simulations and crystallization of the inhibitor-KRASG12C complex highlighted an unprecedented binding mode.

Targeting the Small GTPase Superfamily through Their Regulatory Proteins

The Ras superfamily of small GTPases are guanine-nucleotide-dependent switches essential for numerous cellular processes. Mutations or dysregulation of these proteins are associated with many diseases, but unsuccessful attempts to target the small GTPases directly have resulted in them being classed as "undruggable". The GTP-dependent signaling of these proteins is controlled by their regulators; guanine nucleotide exchange factors (GEFs), GTPase activating proteins (GAPs), and in the Rho and Rab subfamilies, guanine nucleotide dissociation inhibitors (GDIs). This review covers the recent small molecule and biologics strategies to target the small GTPases through their regulators. It seeks to critically re-evaluate recent chemical biology practice, such as the presence of PAINs motifs and the cell-based readout using compounds that are weakly potent or of unknown specificity. It highlights the vast scope of potential approaches for targeting the small GTPases in the future through their regulatory proteins.

Structural And Computational Perspectives of Selectively Targeting Mutant Proteins

Diseases are often caused by mutant proteins. Many drugs have limited effectiveness and/or toxic side effects because of a failure to selectively target the disease-causing mutant variant, rather than the functional wild type protein. Otherwise, the drugs may even target different proteins with similar structural features. Designing drugs that successfully target mutant proteins selectively represents a major challenge. Decades of cancer research have led to an abundance of potential therapeutic targets, often touted to be "master regulators". For many of these proteins, there are no FDA-approved drugs available; for others, off-target effects result in dose-limiting toxicity. Cancer-related proteins are an excellent medium to carry the story of mutant-specific targeting, as the disease is both initiated and sustained by mutant proteins; furthermore, current chemotherapies generally fail at adequate selective distinction. This review discusses some of the challenges associated with selective targeting from a structural biology perspective, as well as some of the developments in algorithm approach and computational workflow that can be applied to address those issues. One of the most widely researched proteins in cancer biology is p53, a tumor suppressor. Here, p53 is discussed as a specific example of a challenging target, with contemporary drugs and methodologies used as examples of burgeoning successes. The oncogene KRAS, which has been described as "undruggable", is another extensively investigated protein in cancer biology. This review also examines KRAS to exemplify progress made towards selective targeting of diseasecausing mutant proteins. Finally, possible future directions relevant to the topic are discussed.

Therapeutic targeting of RAS: New hope for drugging the "undruggable"

RAS is the most frequently mutated oncogene in cancer and a critical driver of oncogenesis. Therapeutic targeting of RAS has been a goal of cancer research for more than 30 years due to its essential role in tumor formation and maintenance. Yet the quest to inhibit this challenging foe has been elusive. Although once considered "undruggable", the struggle to directly inhibit RAS has seen recent success with the development of pharmacological agents that specifically target the KRAS(G12C) mutant protein, which include the first direct RAS inhibitor to gain entry to clinical trials. However, the limited applicability of these inhibitors to G12C-mutant tumors demands further efforts to identify more broadly efficacious RAS inhibitors. Understanding allosteric influences on RAS may open new avenues to inhibit RAS. Here, we provide a brief overview of RAS biology and biochemistry, discuss the allosteric regulation of RAS, and summarize the various approaches to develop RAS inhibitors.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Structural dynamics and allostery of Rab proteins: strategies for drug discovery and design

Rab proteins represent the largest family of the Rab superfamily guanosine triphosphatase (GTPase). Aberrant human Rab proteins are associated with multiple diseases, including cancers and neurological disorders. Rab subfamily members display subtle conformational variations that render specificity in their physiological functions and can be targeted for subfamily-specific drug design. However, drug discovery efforts have not focused much on targeting Rab allosteric non-nucleotide binding sites which are subjected to less evolutionary pressures to be conserved, hence are likely to offer subfamily specificity and may be less prone to undesirable off-target interactions and side effects. To discover druggable allosteric binding sites, Rab structural dynamics need to be first incorporated using multiple experimentally and computationally obtained structures. The high-dimensional structural data may necessitate feature extraction methods to identify manageable representative structures for subsequent analyses. We have detailed state-of-the-art computational methods to (i) identify binding sites using data on sequence, shape, energy, etc., (ii) determine the allosteric nature of these binding sites based on structural ensembles, residue networks and correlated motions and (iii) identify small molecule binders through structure- and ligand-based virtual screening. To benefit future studies for targeting Rab allosteric sites, we herein detail a refined workflow comprising multiple available computational methods, which have been successfully used alone or in combinations. This workflow is also applicable for drug discovery efforts targeting other medically important proteins. Depending on the structural dynamics of proteins of interest, researchers can select suitable strategies for allosteric drug discovery and design, from the resources of computational methods and tools enlisted in the workflow.

Multigram scale synthesis of polycyclic lactones and evaluation of antitumor and other biological properties

An efficient four-step synthesis of tetracyclic lactones from 1,4-benzodioxine-2-carboxylic acid was developed. Ellipticine derivatives exhibit antitumor activity however only a few derivatives without carbazole subunit have been studied to date. Herein, several tetracyclic lactones were synthesized and biologically evaluated. Several compounds (2a, 3a, 4a and 5a) were found to be inhibitors of the Kras-Wnt pathway. The lactone 2a also exerted a potent inhibition of Tau protein translation and was shown to have capacity for CYP1A1-bioactivation. The results obtained are further evidence of the therapeutic potential of tetracyclic lactones related to ellipticine. Molecular modeling studies showed that compound 2a is inserted between helix α3 and α4 of the KRas protein making interactions with the hydrophobic residues Phe90, Glu91, Ile9364, Hie94, Leu133 and Tyr137and a hydrogen bond with residue Arg97.

Structure-Based Analysis of Cryptic-Site Opening

Many proteins in their unbound structures have cryptic sites that are not appropriately sized for drug binding. We consider here 32 proteins from the recently published CryptoSite set with validated cryptic sites, and study whether the sites remain cryptic in all available X-ray structures of the proteins solved without any ligand bound near the sites. It was shown that only few of these proteins have binding pockets that never form without ligand binding. Sites that are cryptic in some structures but spontaneously form in others are also rare. In most proteins the forming of pockets is affected by mutations or ligand binding at locations far from the cryptic site. To further explore these mechanisms, we applied adiabatic biased molecular dynamics simulations to guide the proteins from their ligand-free structures to ligand-bound conformations, and studied the distribution of druggability scores of the pockets located at the cryptic sites.

Discovery of N-(1-Acryloylazetidin-3-yl)-2-(1H-indol-1-yl)acetamides as Covalent Inhibitors of KRASG12C

KRAS regulates many cellular processes including proliferation, survival, and differentiation. Point mutants of KRAS have long been known to be molecular drivers of cancer. KRAS p.G12C, which occurs in approximately 14% of lung adenocarcinomas, 3-5% of colorectal cancers, and low levels in other solid tumors, represents an attractive therapeutic target for covalent inhibitors. Herein, we disclose the discovery of a class of novel, potent, and selective covalent inhibitors of KRAS^G12C identified through a custom library synthesis and screening platform called Chemotype Evolution and structure-based design. Identification of a hidden surface groove bordered by H95/Y96/Q99 side chains was key to the optimization of this class of molecules. Best-in-series exemplars exhibit a rapid covalent reaction with cysteine 12 of GDP-KRAS^G12C with submicromolar inhibition of downstream signaling in a KRAS^G12C-specific manner.

Exploring the Lapse in Druggability: Sequence Analysis, Structural Dynamics and Binding Site Characterization of K-RasG12C Variant, a Feasible Oncotherapeutics Target

BACKGROUND

The difficulty in druggability of K-Ras variant has presented a challenge in the treatment of cancer diseases associated with its dysfunctionality. Despite the identification of different binding sites, limited information exists in the literature about their characteristics. Therefore, identification, crossvalidation and characterization of its druggable sites would aid the design of chemical compounds that will arrest its dysfunctionality related oncogenesis.

OBJECTIVE

This study entails the identification, cross-validation and characterization of K-Ras G12C variant's binding sites for potential druggability, coupled with the elucidation of alterations in 3D conformations and dynamics.

METHOD

Molecular dynamics simulation was carried out on the inactive, the active and the hyperactive K-RasG12Cvariant using the amber software package. The SiteMap software was employed in identifying and characterizing the druggable binding sites while the validation of the binding sites was carried out with the SiteHound and MetaPocket servers.

RESULTS

Four druggable binding sites were identified, validated and characterized based on physicochemical attributes such as size, volume, degree of enclosure or exposure, degree of contact, hydrophobic/hydrophilic character, hydrophobic/hydrophilic balance and hydrogen-bonding features. Conformational studies also revealed that the K-Ras variant exhibited notable structural instability, increased flexibility and a strongly anticorrelated movement compared to the inactive and active wildtype forms.

CONCLUSION

The attributes of the characterized druggable sites will be useful in designing site-specific K-Ras inhibitors for the treatment of K-Ras variant associated cancer diseases.

KRAS-specific inhibition using a DARPin binding to a site in the allosteric lobe

Inhibiting the RAS oncogenic protein has largely been through targeting the switch regions that interact with signalling effector proteins. Here, we report designed ankyrin repeat proteins (DARPins) macromolecules that specifically inhibit the KRAS isoform by binding to an allosteric site encompassing the region around KRAS-specific residue histidine 95 at the helix α3/loop 7/helix α4 interface. We show that these DARPins specifically inhibit KRAS/effector interactions and the dependent downstream signalling pathways in cancer cells. Binding by the DARPins at that region influences KRAS/effector interactions in different ways, including KRAS nucleotide exchange and inhibiting KRAS dimerization at the plasma membrane. These results highlight the importance of targeting the α3/loop 7/α4 interface, a previously untargeted site in RAS, for specifically inhibiting KRAS function.

Multi-target, ensemble-based virtual screening yields novel allosteric KRAS inhibitors at high success rate

RAS mutations account for >15% of all human tumors, and of these ~85% are due to mutations in a particular RAS gene: KRAS. Recent studies revealed that KRAS harbors four druggable allosteric sites. Here, we have (a) used molecular simulations to generate ensembles of wild type and four major oncogenic KRAS mutants (G12V, G12D, G13D, and Q61H); (b) characterized the druggability of each allosteric pocket in each protein; (c) conducted extensive ensemble-based virtual screening using pocket-tailored ligand libraries; (d) prioritized hits through hierarchical postdocking analysis; and (e) validated predicted hits with NMR. Of the 785 diverse potential hits identified by our in silico analysis, we tested 90 for their ability to bind KRAS using NMR and found that nine cause backbone amide chemical shift perturbations of residues near the functionally responsive switch loops, suggesting potential binding. We conducted detailed biophysical analyses on a novel indole-based compound to demonstrate the potential of our workflow to yield lead compounds. We believe the detailed information documented in this work regarding the druggability profile of each allosteric site and the chemical fingerprints of compounds that target them will serve as vital resources for future structure-based drug design efforts against KRAS, a high-value target for cancer therapy.

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.