Exploration of Cryptic Pockets Using Enhanced Sampling Along Normal Modes: A Case Study of KRAS G12D.

Identification of cryptic pockets has the potential to open new therapeutic opportunities by discovering ligand binding sites that remain hidden in static apo structures of a target protein. Moreover, allosteric cryptic pockets can become valuable for designing target-selective ligands when the natural ligand binding sites are conserved in variants of a protein. For example, before an allosteric cryptic pocket was discovered, KRAS was considered undruggable due to its smooth surface and conservation of the GDP/GTP binding pocket across the wild type and oncogenic isoforms. Recent identification of the Switch-II cryptic pocket in the KRASG12C mutant and FDA approval of anticancer drugs targeting this site underscores the importance of cryptic pockets in solving pharmaceutical challenges. Here, we present a newly developed approach for the exploration of cryptic pockets using weighted ensemble molecular dynamics simulations with inherent normal modes as progress coordinates applied to the wild type KRAS and the G12D mutant. We performed extensive all-atomic simulations (>400 µs) with and without several cosolvents (xenon, ethanol, benzene), and analyzed trajectories using three distinct methods to search for potential binding pockets. These methods have been applied as a proof-of-concept to KRAS and have shown they can predict known cryptic binding sites. Furthermore, we performed ligand-binding simulations of a known inhibitor (MRTX1133) to shed light on the nature of cryptic pockets in KRASG12D and the role of conformational selection vs induced-fit mechanism in the formation of these cryptic pockets.

Discovery of mobocertinib, a potent, oral inhibitor of EGFR exon 20 insertion mutations in non-small cell lung cancer.

In the treatment of non-small cell lung cancer (NSCLC), patients harboring exon 20 insertion mutations in the epidermal growth factor receptor (EGFR) gene (EGFR) have few effective therapies because this subset of mutants is generally resistant to most currently approved EGFR inhibitors. This report describes the structure-guided design of a novel series of potent, irreversible inhibitors of EGFR exon 20 insertion mutations, including the V769_D770insASV and D770_N771insSVD mutants. Extensive structure-activity relationship (SAR) studies led to the discovery of mobocertinib (compound 21c), which inhibited growth of Ba/F3 cells expressing the ASV insertion with a half-maximal inhibitory concentration of 11 nM and with selectivity over wild-type EGFR. Daily oral administration of mobocertinib induced tumor regression in a Ba/F3 ASV xenograft mouse model at well-tolerated doses. Mobocertinib was approved in September 2021 for the treatment of adult patients with advanced NSCLC with EGFR exon 20 insertion mutations with progression on or after platinum-based chemotherapy.

Computational and Functional Analyses of HER2 Mutations Reveal Allosteric Activation Mechanisms and Altered Pharmacologic Effects.

Amplification of HER2 can drive the proliferation of cancer cells, and several inhibitors of HER2 have been successfully developed. Recent advances in next-generation sequencing now reveal that HER2 is subject to mutation, with over 2,000 unique variants observed in human cancers. Several examples of oncogenic HER2 mutations have been described, and these primarily occur at allosteric sites outside the ATP-binding site. To identify the full spectrum of oncogenic HER2 driver mutations aside from a few well-studied mutations, we developed mutation-allostery-pharmacology (MAP), an in silico prediction algorithm based on machine learning. By applying this computational approach to 820 single-nucleotide variants, a list of 222 known and potential driver mutations was produced. Of these 222 mutations, 111 were screened by Ba/F3-retrovirus proliferation assays; 37 HER2 mutations were experimentally determined to be driver mutations, comprising 15 previously characterized and 22 newly identified oncogenic mutations. These oncogenic mutations mostly affected allosteric sites in the extracellular domain (ECD), transmembrane domain, and kinase domain of HER2, with only a single mutation in the HER2 orthosteric ATP site. Covalent homodimerization was established as a common mechanism of activation among HER2 ECD allosteric mutations, including the most prevalent HER2 mutation, S310F. Furthermore, HER2 allosteric mutants with enhanced covalent homodimerization were characterized by altered pharmacology that reduces the activity of existing anti-HER2 agents, including the mAb trastuzumab and the tyrosine kinase inhibitor lapatinib. Overall, the MAP-scoring and functional validation analyses provided new insights into the oncogenic activity and therapeutic targeting of HER2 mutations in cancer. SIGNIFICANCE: This study identified new oncogenic HER2 allosteric mutations, including ECD mutations that share covalent dimerization as a mechanism of oncogenicity, suggesting the need for novel inhibitors to treat HER2-mutant cancers.

Discovery of Novel, Potent Inhibitors of Hydroxy Acid Oxidase 1 (HAO1) Using DNA-Encoded Chemical Library Screening.

Inhibition of hydroxy acid oxidase 1 (HAO1) is a strategy to mitigate the accumulation of toxic oxalate that results from reduced activity of alanine-glyoxylate aminotransferase (AGXT) in primary hyperoxaluria 1 (PH1) patients. DNA-Encoded Chemical Library (DECL) screening provided two novel chemical series of potent HAO1 inhibitors, represented by compounds 3-6. Compound 5 was further optimized via various structure-activity relationship (SAR) exploration methods to 29, a compound with improved potency and absorption, distribution, metabolism, and excretion (ADME)/pharmacokinetic (PK) properties. Since carboxylic acid-containing compounds are often poorly permeable and have potential active glucuronide metabolites, we undertook a brief, initial exploration of acid replacements with the aim of identifying non-acid-containing HAO1 inhibitors. Structure-based drug design initiated with Compound 5 led to the identification of a nonacid inhibitor of HAO1, 31, which has weaker potency and increased permeability.

Bispecific Estrogen Receptor α Degraders Incorporating Novel Binders Identified Using DNA-Encoded Chemical Library Screening.

Bispecific degraders (PROTACs) of ERα are expected to be advantageous over current inhibitors of ERα signaling (aromatase inhibitors/SERMs/SERDs) used to treat ER+ breast cancer. Information from DNA-encoded chemical library (DECL) screening provides a method to identify novel PROTAC binding features as the linker positioning, and binding elements are determined directly from the screen. After screening ∼120 billion DNA-encoded molecules with ERα WT and 3 gain-of-function (GOF) mutants, with and without estradiol to identify features that enrich ERα competitively, the off-DNA synthesized small molecule exemplar 7 exhibited nanomolar ERα binding, antagonism, and degradation. Click chemistry synthesis on an alkyne E3 ligase engagers panel and an azide variant of 7 rapidly generated bispecific nanomolar degraders of ERα, with PROTACs 18 and 21 inhibiting ER+ MCF7 tumor growth in a mouse xenograft model of breast cancer. This study validates this approach toward identifying novel bispecific degrader leads from DECL screening with minimal optimization.

Novel Autotaxin Inhibitor for the Treatment of Idiopathic Pulmonary Fibrosis: A Clinical Candidate Discovered Using DNA-Encoded Chemistry.

The activity of the secreted phosphodiesterase autotaxin produces the inflammatory signaling molecule LPA and has been associated with a number of human diseases including idiopathic pulmonary fibrosis (IPF). We screened a single DNA-encoded chemical library (DECL) of 225 million compounds and identified a series of potent inhibitors. Optimization of this series led to the discovery of compound 1 (X-165), a highly potent, selective, and bioavailable small molecule. Cocrystallization of compound 1 with human autotaxin demonstrated that it has a novel binding mode occupying both the hydrophobic pocket and a channel near the autotaxin active site. Compound 1 inhibited the production of LPA in human and mouse plasma at nanomolar levels and showed efficacy in a mouse model of human lung fibrosis. After successfully completing IND-enabling studies, compound 1 was approved by the FDA for a Phase I clinical trial. These results demonstrate that DECL hits can be readily optimized into clinical candidates.

Novel Small Molecule Inhibitors of Choline Kinase Identified by Fragment-Based Drug Discovery.

Choline kinase α (ChoKα) is an enzyme involved in the synthesis of phospholipids and thereby plays key roles in regulation of cell proliferation, oncogenic transformation, and human carcinogenesis. Since several inhibitors of ChoKα display antiproliferative activity in both cellular and animal models, this novel oncogene has recently gained interest as a promising small molecule target for cancer therapy. Here we summarize our efforts to further validate ChoKα as an oncogenic target and explore the activity of novel small molecule inhibitors of ChoKα. Starting from weakly binding fragments, we describe a structure based lead discovery approach, which resulted in novel highly potent inhibitors of ChoKα. In cancer cell lines, our lead compounds exhibit a dose-dependent decrease of phosphocholine, inhibition of cell growth, and induction of apoptosis at low micromolar concentrations. The druglike lead series presented here is optimizable for improvements in cellular potency, drug target residence time, and pharmacokinetic parameters. These inhibitors may be utilized not only to further validate ChoKα as antioncogenic target but also as novel chemical matter that may lead to antitumor agents that specifically interfere with cancer cell metabolism.

Discovery of Brigatinib (AP26113), a Phosphine Oxide-Containing, Potent, Orally Active Inhibitor of Anaplastic Lymphoma Kinase.

In the treatment of echinoderm microtubule-associated protein-like 4 (EML4)-anaplastic lymphoma kinase positive (ALK+) non-small-cell lung cancer (NSCLC), secondary mutations within the ALK kinase domain have emerged as a major resistance mechanism to both first- and second-generation ALK inhibitors. This report describes the design and synthesis of a series of 2,4-diarylaminopyrimidine-based potent and selective ALK inhibitors culminating in identification of the investigational clinical candidate brigatinib. A unique structural feature of brigatinib is a phosphine oxide, an overlooked but novel hydrogen-bond acceptor that drives potency and selectivity in addition to favorable ADME properties. Brigatinib displayed low nanomolar IC50s against native ALK and all tested clinically relevant ALK mutants in both enzyme-based biochemical and cell-based viability assays and demonstrated efficacy in multiple ALK+ xenografts in mice, including Karpas-299 (anaplastic large-cell lymphomas [ALCL]) and H3122 (NSCLC). Brigatinib represents the most clinically advanced phosphine oxide-containing drug candidate to date and is currently being evaluated in a global phase 2 registration trial.

Ponatinib inhibits polyclonal drug-resistant KIT oncoproteins and shows therapeutic potential in heavily pretreated gastrointestinal stromal tumor (GIST) patients.

PURPOSE: KIT is the major oncogenic driver of gastrointestinal stromal tumors (GIST). Imatinib, sunitinib, and regorafenib are approved therapies; however, efficacy is often limited by the acquisition of polyclonal secondary resistance mutations in KIT, with those located in the activation (A) loop (exons 17/18) being particularly problematic. Here, we explore the KIT-inhibitory activity of ponatinib in preclinical models and describe initial characterization of its activity in patients with GIST. EXPERIMENTAL DESIGN: The cellular and in vivo activities of ponatinib, imatinib, sunitinib, and regorafenib against mutant KIT were evaluated using an accelerated mutagenesis assay and a panel of engineered and GIST-derived cell lines. The ponatinib-KIT costructure was also determined. The clinical activity of ponatinib was examined in three patients with GIST previously treated with all three FDA-approved agents. RESULTS: In engineered and GIST-derived cell lines, ponatinib potently inhibited KIT exon 11 primary mutants and a range of secondary mutants, including those within the A-loop. Ponatinib also induced regression in engineered and GIST-derived tumor models containing these secondary mutations. In a mutagenesis screen, 40 nmol/L ponatinib was sufficient to suppress outgrowth of all secondary mutants except V654A, which was suppressed at 80 nmol/L. This inhibitory profile could be rationalized on the basis of structural analyses. Ponatinib (30 mg daily) displayed encouraging clinical activity in two of three patients with GIST. CONCLUSION: Ponatinib possesses potent activity against most major clinically relevant KIT mutants and has demonstrated preliminary evidence of activity in patients with refractory GIST. These data strongly support further evaluation of ponatinib in patients with GIST.

Fragment growing and linking lead to novel nanomolar lactate dehydrogenase inhibitors.

Lactate dehydrogenase A (LDH-A) catalyzes the interconversion of lactate and pyruvate in the glycolysis pathway. Cancer cells rely heavily on glycolysis instead of oxidative phosphorylation to generate ATP, a phenomenon known as the Warburg effect. The inhibition of LDH-A by small molecules is therefore of interest for potential cancer treatments. We describe the identification and optimization of LDH-A inhibitors by fragment-based drug discovery. We applied ligand based NMR screening to identify low affinity fragments binding to LDH-A. The dissociation constants (K(d)) and enzyme inhibition (IC(50)) of fragment hits were measured by surface plasmon resonance (SPR) and enzyme assays, respectively. The binding modes of selected fragments were investigated by X-ray crystallography. Fragment growing and linking, followed by chemical optimization, resulted in nanomolar LDH-A inhibitors that demonstrated stoichiometric binding to LDH-A. Selected molecules inhibited lactate production in cells, suggesting target-specific inhibition in cancer cell lines.

Application of MM-GB/SA and WaterMap to SRC Kinase Inhibitor Potency Prediction.

WaterMap and MM-GB/SA scoring methods were applied to an extensive congeneric series of small-molecule SRC inhibitors with high-quality enzyme data and well characterized binding modes to compare the performance of these scoring methods in this data set and to provide insight into the relative strengths of each method. Only minor conformational changes in SRC bound with representative DFG-in class of inhibitors were demonstrated in previous studies; thus, the protein flexibility that normally presents a challenge to pose and potency predictions was minimized in this model system. While WaterMap correctly recognized major trends in the SAR of this series, MM-GB/SA performed better in ranking the relative ligand affinities. The different scoring methods were further analyzed to determine which aspects of series SAR were more amenable to MM-GB/SA than WaterMap scoring.