Accurately measuring compound binding affinities is key to driving the pharmaceutical development process. Rigorous physics-based in silico approaches, particularly alchemical free energy methods, have become a gold-standard tool for estimating compound affinity changes. Here we present the results of a large-scale precompetitive collaborative assessment of relative binding free energy (RBFE) calculations generated by 15 pharmaceutical companies. We evaluate an open-source and MIT-licensed RBFE protocol from the Open Free Energy (OpenFE) ecosystem across both public and blinded private datasets, encompassing over 1,700 ligands in total. For the public dataset, the weighted RMSE across the 58 systems was 1.73(1.53)(1.96) kcal/mol, with 10 of the systems reaching sub-kcal/mol accuracy. For the private dataset, the weighted RMSE across the 37 systems was 2.44(1.94)(3.06) kcal/mol, with only 2 of the systems reaching sub-kcal/mol accuracy, reflecting the increased complexity of real-world drug discovery. The protocol's performance was system-dependent, with no single dominant error source, indicating that accuracy is primarily influenced by input quality and transformation type. Overall, these benchmark results are encouraging and indicate that OpenFE is ready for large-scale industrial applications, with an "out-of-the-box" accuracy that approaches that of commercial solutions. While comparison against published FEP+ results, which were obtained after manual parameter optimization, shows comparable ranking statistics, a gap remains in error statistics, for which we outline possible paths toward improvement. The protocol meets key criteria required for production use in an industrial setting: it shows robust performance, generates reproducible results, and achieves both sufficient throughput and rapid convergence.
KRAS G12C is one of the most prevalent oncogenic mutations in nonsmall cell lung cancer. Herein we describe the discovery and optimization of divarasib (GDC-6036), an orally available, highly potent, and selective covalent KRAS G12C inhibitor. We demonstrate a significant noncovalent binding component of divarasib that contributes to its potency and rapid kinetics. Divarasib has greater potency and kinetics of alkylation compared with other KRAS G12C inhibitors in vitro and shows robust tumor growth inhibition in multiple KRAS G12C-positive cell lines.
KRAS mutations occur frequently in advanced non-small cell lung cancer (aNSCLC); the G12C mutation is the most prevalent. Alterations in STK11 or KEAP1 commonly co-occur with KRAS mutations in aNSCLC. Using real-world data, we assessed the effect of KRAS G12C mutation with or without STK11 and/or KEAP1 mutations on overall survival (OS) in patients with aNSCLC receiving cancer immunotherapy (CIT), chemotherapy, or both in first line (1L) and second line (2L). Patients diagnosed with aNSCLC between January 2011 and March 2020 in a clinico-genomic database were included. Cox proportional hazards models adjusted for left truncation, baseline demographics and clinical characteristics were used to analyze the effect of STK11 and/or KEAP1 co-mutational status on OS in patients with KRAS wild-type (WT) or G12C mutation. Of 2715 patients with aNSCLC without other actionable driver mutations, 1344 (49.5
The growing opportunities recognized for covalent drug inhibitors, like KRAS G12C inhibitors, are driving the need for mass spectrometry methods that can quickly and robustly measure therapeutic drug activity in vivo for drug discovery research and development. Effective front-end sample preparation is critical for proteins extracted from tumors but is generally labor intensive and impractical for large sample numbers typical in pharmacodynamic (PD) studies. Herein, we describe an automated and integrated sample preparation method for the measurement of activity levels of KRAS G12C drug inhibitor alkylation from complex tumor samples involving high throughput detergent removal and preconcentration followed by quantitation using mass spectrometry. We introduce a robust assay with an average intra-assay coefficient of variation (CV) of 4% and an interassay CV of 6% obtained from seven studies, enabling us to understand the relationship between KRAS G12C target occupancy and the therapeutic PD effect from mouse tumor samples. Further, the data demonstrated that the drug candidate GDC-6036, a KRAS G12C covalent inhibitor, shows dose-dependent target inhibition (KRAS G12C alkylation) and MAPK pathway inhibition, which correlate with high antitumor potency in the MIA PaCa-2 pancreatic xenograft model.
KRAS is one of the most frequently mutated oncogenes, with KRAS G12C recently becoming an actionable target for small molecule intervention. GDC-6036 is an investigational KRAS G12C inhibitor that acts by irreversibly binding to the switch II pocket of KRAS G12C when in the inactive GDP-bound state, thereby blocking GTP binding and activation. Assessing target engagement is an essential component of clinical drug development, helping to demonstrate mechanistic activity, guide dose selection, understand pharmacodynamics as it relates to clinical response, and explore resistance. Here, we report the development of an ultra-sensitive approach for assessing KRAS G12C engagement. Immunoaffinity enrichment with a commercially available anti-RAS antibody was combined with a targeted 2D-LC-MS/MS technique to quantify both free and GDC-6036-bound KRAS G12C proteins. A KRAS G12C-positive non-small cell lung cancer xenograft model was dosed with GDC-6036 to assess the feasibility of this assay for analyzing small core needle biopsies. As predicted, dose-dependent KRAS G12C engagement was observed. To date, a sensitivity of 0.08 fmol/μg of total protein has been achieved for both free and GDC-6036-bound KRAS G12C with as little as 4 μg of total protein extracted from human tumor samples. This sub-fmol/μg level of sensitivity provides a powerful potential approach to assess covalent inhibitor target engagement at the site of action using core needle tumor biopsies from clinical studies.
Small molecule inhibitors that target the phosphatidylinositol 3-kinase (PI3K) signaling pathway have received significant interest for the treatment of cancers. The class I isoform PI3Kα is most commonly associated with solid tumors via gene amplification or activating mutations. However, inhibitors demonstrating both PI3K isoform and mutant specificity have remained elusive. Herein, we describe the optimization and characterization of a series of benzoxazepin-oxazolidinone ATP-competitive inhibitors of PI3Kα which also induce the selective degradation of the mutant p110α protein, the catalytic subunit of PI3Kα. Structure-based design informed isoform-specific interactions within the binding site, leading to potent inhibitors with greater than 300-fold selectivity over the other Class I PI3K isoforms. Further optimization of pharmacokinetic properties led to excellent in vivo exposure and efficacy and the identification of clinical candidate GDC-0077 (inavolisib, 32), which is now under evaluation in a Phase III clinical trial as a treatment for patients with PIK3CA-mutant breast cancer.
Small molecules that stabilize inactive protein conformations are an underutilized strategy for drugging dynamic or otherwise intractable proteins. To facilitate the discovery and characterization of such inhibitors, we created a screening platform to identify conformation-locking antibodies for molecular probes (CLAMPs) that distinguish and induce rare protein conformational states. Applying the approach to KRAS, we discovered CLAMPs that recognize the open conformation of KRASG12C stabilized by covalent inhibitors. One CLAMP enables the visualization of KRASG12C covalent modification in vivo and can be used to investigate response heterogeneity to KRASG12C inhibitors in patient tumors. A second CLAMP enhances the affinity of weak ligands binding to the KRASG12C switch II region (SWII) by stabilizing a specific conformation of KRASG12C, thereby enabling the discovery of such ligands that could serve as leads for the development of drugs in a high-throughput screen. We show that combining the complementary properties of antibodies and small molecules facilitates the study and drugging of dynamic proteins.
Structure-based design was utilized to optimize 6,6-diaryl substituted dihydropyrone and hydroxylactam to obtain inhibitors of lactate dehydrogenase (LDH) with low nanomolar biochemical and single-digit micromolar cellular potencies. Surprisingly the replacement of a phenyl with a pyridyl moiety in the chemical structure revealed a new binding mode for the inhibitors with subtle conformational change of the LDHA active site. This led to the identification of a potent, cell-active hydroxylactam inhibitor exhibiting an in vivo pharmacokinetic profile suitable for mouse tumor xenograft study.
Mutations in KRAS are the most common oncogenic driver mutations in human cancers. The KRAS G12C mutation is one of the most prevalent KRAS mutations, present in approximately 12% of non-small cell lung cancer, 4% of colorectal cancer, and up to 4% of other cancer types. The recent discovery and development of covalent KRAS G12C inhibitors provides an opportunity to inhibit what has historically been considered to be an “undruggable” target. GDC-6036 is an orally bioavailable, highly potent and selective KRAS G12C inhibitor, with a median IC50 in the sub-nanomolar range and greater than 18,000-fold selectivity for G12C versus non-G12C cell lines. GDC-6036 demonstrates greater potency and selectivity compared with other KRAS G12C inhibitors in vitro, and complete tumor growth inhibition in multiple KRAS G12C-positive cell lines and in xenograft mouse models. We will highlight the research program that led to the discovery and optimization of GDC-6036, which is currently in clinical development. Citation Format: Hans Purkey. Discovery of GDC-6036, a clinical stage treatment for KRAS G12C-positive cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr ND11.
Abstract The discovery of covalent inhibitors binding the switch II (SWII) pocket has enabled therapeutic intervention in KRASG12C driven tumors and represents a milestone in targeting KRAS-driven cancers. However, the transient nature and high energetic barrier required for binding this pocket has been an obstacle in successfully targeting other KRAS mutant oncoproteins. We report the discovery of KRAS Conformation Locking Antibodies for Molecular Probe discovery (CLAMP)s that specifically recognize the unique conformation of KRASG12C induced by covalent inhibitors. KRAS CLAMPs enable single cell resolution of covalent inhibitor-bound KRASG12C in cells and in vivo tumor models, providing a biomarker for direct target engagement of KRASG12C inhibition. KRAS CLAMPs bind multiple KRAS mutants and stabilize an open conformation of the SWII pocket increasing the affinity of weak non-covalent SWII pocket ligands. This work provides new insights into KRASG12C upon treatment with covalent inhibitors and offers a path towards targeting the SWII pocket in other RAS mutants.
Bruton's tyrosine kinase (Btk) is thought to play a pathogenic role in chronic immune diseases such as rheumatoid arthritis and lupus. While covalent, irreversible Btk inhibitors are approved for treatment of hematologic malignancies, they are not approved for autoimmune indications. In efforts to develop additional series of reversible Btk inhibitors for chronic immune diseases, we sought to differentiate from our clinical stage inhibitor fenebrutinib using cyclopropyl amide isosteres of the 2-amino-pyridyl group to occupy the flat, lipophilic H2 pocket. While drug-like properties were retained-and in some cases improved-a safety liability in the form of hERG inhibition was observed. When a fluorocyclopropyl amide was incorporated, Btk and off-target activity was found to be stereodependent and a lead compound was identified in the form of the (R,R)- stereoisomer.
Abstract Activating mutations in PIK3CA are among the most significant oncogenic events across all cancers, making it an important target for drug development. Yet the application of PI3K inhibitors in the clinic has been limited by the difficulty of achieving an adequate therapeutic window, due to the critical role that PI3K signaling plays in normal physiologic processes, such as glucose homeostasis. In theory, the therapeutic window could be improved if it were possible to design mutant selective inhibitors, as has been demonstrated with other oncogenes such as EGFR. However, unlike EGFR, the most predominant PIK3CA activating mutations do not reside in the kinase active site, presenting a major challenge for rational structure-based design. Nevertheless, it was recently shown that the PI3K inhibitor taselisib is able to achieve modest levels of mutant selectivity both across cancer lines as well as in cell lines that were engineered to express mutant or wild-type PIK3CA. Taselisib was also shown to selectively induce degradation of mutant versus wild-type PIK3CA, leading to the speculation that this degradation may be responsible for the observed selectivity. In order to better understand the origins of mutant selectivity for taselisib and several other PIK3CA inhibitors, we assessed these inhibitors in a variety of biophysical and biochemical assays under conditions designed to mimic physiologic settings. In parallel, we also investigated the mechanistic basis of this selectivity in our engineered cell lines. Our results are consistent with the hypothesis that selective degradation of mutant PIK3CA is the predominant mechanism underlying mutant selectivity for this class of PIK3CA active site inhibitors. This abstract is also being presented as Poster B03. Citation Format: Lan Nguyen, Kyle Edgar, Kyung Song, Stephen Schmidt, Victorai Schutz, Noriko Ishisoko, Eric Torres, Akash Das, Divya Murali, Steve Sideris, Timothy Wendorff, Matt Saabye, Hans Purkey, Jawahar Sudhamsu, Steven Staben, Emily Hanan, Georgia Hatzivassiliou, Lori Friedman, Nicholas F. Endres. Selective degradation of mutant PIK3CA promotes increased mutant specificity in a subset of PI3K ATP-competitive inhibitors [abstract]. In: Proceedings of the AACR Special Conference on Targeting PI3K/mTOR Signaling; 2018 Nov 30-Dec 8; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Res 2020;18(10_Suppl):Abstract nr PR03.
CDK4 and CDK6 are kinases with similar sequences that regulate cell cycle progression and are validated targets in the treatment of cancer. Glioblastoma is characterized by a high frequency of CDKN2A/CCND2/CDK4/CDK6 pathway dysregulation, making dual inhibition of CDK4 and CDK6 an attractive therapeutic approach for this disease. Abemaciclib, ribociclib, and palbociclib are approved CDK4/6 inhibitors for the treatment of HR+/HER2- breast cancer, but these drugs are not expected to show strong activity in brain tumors due to poor blood brain barrier penetration. Herein, we report the identification of a brain-penetrant CDK4/6 inhibitor derived from a literature molecule with low molecular weight and topological polar surface area (MW = 285 and TPSA = 66 angstrom(2)), but lacking the CDK2/1 selectivity profile due to the absence of a basic amine. Removal of a hydrogen bond donor via cyclization of the pyrazole allowed for the introduction of basic and semi-basic amines, while maintaining in many cases efflux ratios reasonable for a CNS program. Ultimately, a basic spiroazetidine (cpK(a) = 8.8) was identified that afforded acceptable selectivity over anti-target CDK1 while maintaining brain-penetration in vivo (mouse K-p,K-uu = 0.20-0.59). To probe the potency and selectivity, our lead compound was evaluated in a panel of glioblastoma cell lines. Potency comparable to abemaciclib was observed in Rb-wild type lines U87MG, DBTRG-05MG, A172, and T98G, while Rb-deficient cell lines SF539 and M059J exhibited a lack of sensitivity.
The pan-proteasome inhibitor bortezomib demonstrated clinical efficacy in off-label trials of Systemic Lupus Erythematosus. One potential mechanism of this clinical benefit is from the depletion of pathogenic immune cells (plasmablasts and plasmacytoid dendritic cells). However, bortezomib is cytotoxic against nonimmune cells, which limits its use for autoimmune diseases. An attractive alternative is to selectively inhibit the immune cell-specific immunoproteasome to deplete pathogenic immune cells and spare nonhematopoietic cells. Here, we disclose the development of highly subunit-selective immunoproteasome inhibitors using insights obtained from the first bona fide human immunoproteasome cocrystal structures. Evaluation of these inhibitors revealed that immunoproteasome-specific inhibition does not lead to immune cell death as anticipated and that targeting viability requires inhibition of both immuno- and constitutive proteasomes. CRISPR/Cas9-mediated knockout experiments confirmed upregulation of the constitutive proteasome upon disruption of the immunoproteasome, protecting cells from death. Thus, immunoproteasome inhibition alone is not a suitable approach to deplete immune cells.
There is a critical need for new antibacterial strategies to counter the growing problem of antibiotic resistance. In Gram-negative bacteria, the outer membrane (OM) provides a protective barrier against antibiotics and other environmental insults. The outer leaflet of the outer membrane is primarily composed of lipopolysaccharide (LPS). Outer membrane biogenesis presents many potentially compelling drug targets as this pathway is absent in higher eukaryotes. Most proteins involved in LPS biosynthesis and transport are essential; however, few compounds have been identified that inhibit these proteins. The inner membrane ABC transporter MsbA carries out the first essential step in the trafficking of LPS to the outer membrane. We conducted a biochemical screen for inhibitors of MsbA and identified a series of quinoline compounds that kill Escherichia coli through inhibition of its ATPase and transport activity, with no loss of activity against clinical multidrug-resistant strains. Identification of these selective inhibitors indicates that MsbA is a viable target for new antibiotics, and the compounds we identified serve as useful tools to further probe the LPS transport pathway in Gram-negative bacteria.
The movement of core-lipopolysaccharide across the inner membrane of Gram-negative bacteria is catalysed by an essential ATP-binding cassette transporter, MsbA. Recent structures of MsbA and related transporters have provided insights into the molecular basis of active lipid transport; however, structural information about their pharmacological modulation remains limited. Here we report the 2.9 Å resolution structure of MsbA in complex with G907, a selective small-molecule antagonist with bactericidal activity, revealing an unprecedented mechanism of ABC transporter inhibition. G907 traps MsbA in an inward-facing, lipopolysaccharide-bound conformation by wedging into an architecturally conserved transmembrane pocket. A second allosteric mechanism of antagonism occurs through structural and functional uncoupling of the nucleotide-binding domains. This study establishes a framework for the selective modulation of ABC transporters and provides rational avenues for the design of new antibiotics and other therapeutics targeting this protein family.
The pan-proteasome inhibitor bortezomib has demonstrated sustained clinical efficacy in off label trials in SLE patients by depleting pathogenic immune cells including highly secretory plasmablasts and plasmacytoid dendritic cells. But bortezomib also effects non-immune cells including neurons which raises significant safety concerns thus limiting the feasibility of pan-proteasome inhibitors as a treatment for autoimmune disease. As an alternative, targeting the immune cell specific immunoproteasome has been proposed as a promising therapy to deplete pathogenic immune cells while avoiding toxicity issues. Using highly specific immunoproteasome inhibitors, we show that this strategy does not lead to immune cell death as anticipated. We generated a large panel of small molecule inhibitors with varying specificity to the β5 subunit of the conventional and immuno- proteasome, and found that effects on viability of plasmablasts and pDCs correlate with the inhibition of the conventional proteasome, not the immunoproteasome. Our data indicates that the immune cells upregulate the conventional proteasomal subunits and prevent the accumulation of ubiquitinated protein. Widely used immunoproteasome inhibitors also block the conventional proteasome which, our data would indicate, is key to their ability to deplete immune cells.