JDQ443, a Structurally Novel, Pyrazole-Based, Covalent Inhibitor of KRASG12C for the Treatment of Solid Tumors.

Rapid emergence of tumor resistance via RAS pathway reactivation has been reported from clinical studies of covalent KRASG12C inhibitors. Thus, inhibitors with broad potential for combination treatment and distinct binding modes to overcome resistance mutations may prove beneficial. JDQ443 is an investigational covalent KRASG12C inhibitor derived from structure-based drug design followed by extensive optimization of two dissimilar prototypes. JDQ443 is a stable atropisomer containing a unique 5-methylpyrazole core and a spiro-azetidine linker designed to position the electrophilic acrylamide for optimal engagement with KRASG12C C12. A substituted indazole at pyrazole position 3 results in novel interactions with the binding pocket that do not involve residue H95. JDQ443 showed PK/PD activity in vivo and dose-dependent antitumor activity in mouse xenograft models. JDQ443 is now in clinical development, with encouraging early phase data reported from an ongoing Phase Ib/II clinical trial (NCT04699188).

Discovery, Preclinical Characterization, and Early Clinical Activity of JDQ443, a Structurally Novel, Potent, and Selective Covalent Oral Inhibitor of KRASG12C.

Covalent inhibitors of KRASG12C have shown antitumor activity against advanced/metastatic KRASG12C-mutated cancers, though resistance emerges and additional strategies are needed to improve outcomes. JDQ443 is a structurally unique covalent inhibitor of GDP-bound KRASG12C that forms novel interactions with the switch II pocket. JDQ443 potently inhibits KRASG12C-driven cellular signaling and demonstrates selective antiproliferative activity in KRASG12C-mutated cell lines, including those with G12C/H95 double mutations. In vivo, JDQ443 induces AUC exposure-driven antitumor efficacy in KRASG12C-mutated cell-derived (CDX) and patient-derived (PDX) tumor xenografts. In PDX models, single-agent JDQ443 activity is enhanced by combination with inhibitors of SHP2, MEK, or CDK4/6. Notably, the benefit of JDQ443 plus the SHP2 inhibitor TNO155 is maintained at reduced doses of either agent in CDX models, consistent with mechanistic synergy. JDQ443 is in clinical development as monotherapy and in combination with TNO155, with both strategies showing antitumor activity in patients with KRASG12C-mutated tumors. SIGNIFICANCE: JDQ443 is a structurally novel covalent KRASG12C inhibitor with a unique binding mode that demonstrates potent and selective antitumor activity in cell lines and in vivo models. In preclinical models and patients with KRASG12C-mutated malignancies, JDQ443 shows potent antitumor activity as monotherapy and in combination with the SHP2 inhibitor TNO155. This article is highlighted in the In This Issue feature, p. 1397.

ReSCoSS: a flexible quantum chemistry workflow identifying relevant solution conformers of drug-like molecules.

Conformational equilibria are at the heart of drug design, yet their energetic description is often hampered by the insufficient accuracy of low-cost methods. Here we present a flexible and semi-automatic workflow based on quantum chemistry, ReSCoSS, designed to identify relevant conformers and predict their equilibria across different solvent environments in the Conductor-like Screening Model for Real Solvents (COSMO-RS) framework. We demonstrate the utility and accuracy of the workflow through conformational case studies on several drug-like molecules from literature where relevant conformations are known. We further show that including ReSCoSS conformers significantly improves COSMO-RS based predictions of physicochemical properties over single-conformation approaches. ReSCoSS has found broad adoption in the in-house drug discovery and development work streams and has contributed to establishing quantum-chemistry methods as a strategic pillar in ligand discovery.

Discovery and Optimization of Novel SUCNR1 Inhibitors: Design of Zwitterionic Derivatives with a Salt Bridge for the Improvement of Oral Exposure.

G-protein-coupled receptor SUCNR1 (succinate receptor 1 or GPR91) senses the citric cycle intermediate succinate and is implicated in various pathological conditions such as rheumatoid arthritis, liver fibrosis, or obesity. Here, we describe a novel SUCNR1 antagonist scaffold discovered by high-throughput screening. The poor permeation and absorption properties of the most potent compounds, which were zwitterionic in nature, could be improved by the formation of an internal salt bridge, which helped in shielding the two opposite charges and thus also the high polarity of zwitterions with separated charges. The designed compounds containing such a salt bridge reached high oral bioavailability and oral exposure. We believe that this principle could find a broad interest in the medicinal chemistry field as it can be useful not only for the modulation of properties in zwitterionic compounds but also in acidic or basic compounds with poor permeation.

Fragment-Based Discovery of Non-bisphosphonate Binders of Trypanosoma brucei Farnesyl Pyrophosphate Synthase.

Trypanosoma brucei is the causative agent of human African trypanosomiasis (HAT). Nitrogen-containing bisphosphonates, a current treatment for bone diseases, have been shown to block the growth of the T. brucei parasites by inhibiting farnesyl pyrophosphate synthase (FPPS); however, due to their poor pharmacokinetic properties, they are not well suited for antiparasitic therapy. Recently, an allosteric binding pocket was discovered on human FPPS, but its existence on trypanosomal FPPS was unclear. We applied NMR and X-ray fragment screening to T. brucei FPPS and report herein on four fragments bound to this previously unknown allosteric site. Surprisingly, non-bisphosphonate active-site binders were also identified. Moreover, fragment screening revealed a number of additional binding sites. In an early structure-activity relationship (SAR) study, an analogue of an active-site binder was unexpectedly shown to bind to the allosteric site. Overlaying identified fragment binders of a parallel T. cruzi FPPS fragment screen with the T. brucei FPPS structure, and medicinal chemistry optimisation based on two binders revealed another example of fragment "pocket hopping". The discovery of binders with new chemotypes sets the framework for developing advanced compounds with pharmacokinetic properties suitable for the treatment of parasitic infections by inhibition of FPPS in T. brucei parasites.

Benchmarking Different QM Levels for Usage with COSMO-RS.

The COSMO-RS method is an established method for the prediction of fluid phase properties such as activity coefficients, liquid-liquid equilibria, and free energy of solvation. It is also frequently used in quantum chemistry-based chemical reaction modeling to predict the solvation contribution to the reactions. The COSMOtherm software, which features the currently most advanced implementation of COSMO-RS, is based on quantum chemical COSMO calculations using the BP functional with the def2-TZVPD basis set. As the accuracy of COSMO-RS depends on the accuracy of the underlying quantum chemical (QC) calculation, it is important to validate the currently used level against other common, presumably superior, approaches such as the more recently developed M06-2x hybrid density functional or wave function methods such as MP2. As compared to other applications where the electronic energy is the most important result of the QC calculation, the COSMO-RS method has a much higher dependence on the molecular polarity and thus the electron density distribution. We find that MP2, PBE0, and M06-2x perform slightly worse in fully reparametrized COSMO-RS with respect to the prediction of experimentally measured properties like pKa or logP. Although MP2 was reported to yield better polarities than most DFT functionals for spin unpolarized molecules, this theoretical advantage does not manifest in a practical benefit for the prediction of thermodynamic properties with a refitted COSMO-RS parameter set. Other pure DFT functionals such as PBE or TPSS can be used instead of BP, but again, no practical advantage is expected as they yield extremely similar polarities to the original BP calculations.

High accuracy quantum-chemistry-based calculation and blind prediction of macroscopic pKa values in the context of the SAMPL6 challenge.

Recent advances in the development of low-cost quantum chemical methods have made the prediction of conformational preferences and physicochemical properties of medium-sized drug-like molecules routinely feasible, with significant potential to advance drug discovery. In the context of the SAMPL6 challenge, macroscopic pKa values were blindly predicted for a set of 24 of such molecules. In this paper we present two similar quantum chemical based approaches based on the high accuracy calculation of standard reaction free energies and the subsequent determination of those pKa values via a linear free energy relationship. Both approaches use extensive conformational sampling and apply hybrid and double-hybrid density functional theory with continuum solvation to calculate free energies. The blindly calculated macroscopic pKa values were in excellent agreement with the experiment.

Mutations in sphingosine-1-phosphate lyase cause nephrosis with ichthyosis and adrenal insufficiency.

Steroid-resistant nephrotic syndrome (SRNS) causes 15% of chronic kidney disease cases. A mutation in 1 of over 40 monogenic genes can be detected in approximately 30% of individuals with SRNS whose symptoms manifest before 25 years of age. However, in many patients, the genetic etiology remains unknown. Here, we have performed whole exome sequencing to identify recessive causes of SRNS. In 7 families with SRNS and facultative ichthyosis, adrenal insufficiency, immunodeficiency, and neurological defects, we identified 9 different recessive mutations in SGPL1, which encodes sphingosine-1-phosphate (S1P) lyase. All mutations resulted in reduced or absent SGPL1 protein and/or enzyme activity. Overexpression of cDNA representing SGPL1 mutations resulted in subcellular mislocalization of SGPL1. Furthermore, expression of WT human SGPL1 rescued growth of SGPL1-deficient dpl1Δ yeast strains, whereas expression of disease-associated variants did not. Immunofluorescence revealed SGPL1 expression in mouse podocytes and mesangial cells. Knockdown of Sgpl1 in rat mesangial cells inhibited cell migration, which was partially rescued by VPC23109, an S1P receptor antagonist. In Drosophila, Sply mutants, which lack SGPL1, displayed a phenotype reminiscent of nephrotic syndrome in nephrocytes. WT Sply, but not the disease-associated variants, rescued this phenotype. Together, these results indicate that SGPL1 mutations cause a syndromic form of SRNS.

Discovery of novel pyrrolidineoxy-substituted heteroaromatics as potent and selective PI3K delta inhibitors with improved physicochemical properties.

In the recent years, PI3Kδ has emerged as a promising target for the treatment of B- and T-cell mediated inflammatory diseases. We present a cellular assay activity analysis for our previously reported 4,6-diaryl quinazoline PI3Kδ inhibitor series that suggests an optimal logP range between 2 and 3. We discovered novel analogues in this lipophilicity space that feature a chiral pyrrolidineoxy-group as a replacement for the position-4 aromatic ring of 4,6-diaryl quinazolines. These Fsp3 enriched derivatives retain potency and selectivity towards PI3Kδ. Compared to 4,6-diaryl quinazolines, their permeability profile is improved and molecular weight as well as PSA are reduced. These modifications offer additional possibilities for derivative generation in a favorable physicochemical property space and thus increase the chances to identify a clinical candidate.

Harnessing Fluorine-Sulfur Contacts and Multipolar Interactions for the Design of p53 Mutant Y220C Rescue Drugs.

Many oncogenic mutants of the tumor suppressor p53 are conformationally unstable, including the frequently occurring Y220C mutant. We have previously developed several small-molecule stabilizers of this mutant. One of these molecules, PhiKan083, 1-(9-ethyl-9H-carbazole-3-yl)-N-methylmethanamine, binds to a mutation-induced surface crevice with a KD = 150 µM, thereby increasing the melting temperature of the protein and slowing its rate of aggregation. Incorporation of fluorine atoms into small molecule ligands can substantially improve binding affinity to their protein targets. We have, therefore, harnessed fluorine-protein interactions to improve the affinity of this ligand. Step-wise introduction of fluorines at the carbazole ethyl anchor, which is deeply buried within the binding site in the Y220C-PhiKan083 complex, led to a 5-fold increase in affinity for a 2,2,2-trifluoroethyl anchor (ligand efficiency of 0.3 kcal mol(-1) atom(-1)). High-resolution crystal structures of the Y220C-ligand complexes combined with quantum chemical calculations revealed favorable interactions of the fluorines with protein backbone carbonyl groups (Leu145 and Trp146) and the sulfur of Cys220 at the mutation site. Affinity gains were, however, only achieved upon trifluorination, despite favorable interactions of the mono- and difluorinated anchors with the binding pocket, indicating a trade-off between energetically favorable protein-fluorine interactions and increased desolvation penalties. Taken together, the optimized carbazole scaffold provides a promising starting point for the development of high-affinity ligands to reactivate the tumor suppressor function of the p53 mutant Y220C in cancer cells.

Arylpiperidines as a new class of oxidosqualene cyclase inhibitors.

The cyclization of oxidosqualene to lanosterol, catalyzed by the enzyme oxidosqualene cyclase (OSC), goes through a number of carbocationic high energy intermediates (HEI), and mimicking these intermediates is a promising approach for the development of OSC inhibitors. 3-Arylpiperidines (or tetrahydropyridines) were designed as steroidomimetic rings A + C equivalents containing two protonable amino groups for mimicking both the pro-C4 HEI and the pro-C20 HEI of the OSC-mediated cyclization cascade. Inhibitory activity is strongly dependent on the nature of the lipophilic substituent representing an equivalent of the sterol side chain. Here aromatic residues (substituted benzyl, cinnamyl, naphthylmethyl) were found to be most suitable. Docking experiments on a first optimized 3-arylpiperidine compound led to an isomeric 4-arylpiperidine with submicromolar activity on human OSC. This inhibitor reduced total cholesterol biosynthesis in a cellular assay with an IC50 value of 0.26 µM.

Exploiting Transient Protein States for the Design of Small-Molecule Stabilizers of Mutant p53.

The destabilizing p53 cancer mutation Y220C creates an extended crevice on the surface of the protein that can be targeted by small-molecule stabilizers. Here, we identify different classes of small molecules that bind to this crevice and determine their binding modes by X-ray crystallography. These structures reveal two major conformational states of the pocket and a cryptic, transiently open hydrophobic subpocket that is modulated by Cys220. In one instance, specifically targeting this transient protein state by a pyrrole moiety resulted in a 40-fold increase in binding affinity. Molecular dynamics simulations showed that both open and closed states of this subsite were populated at comparable frequencies along the trajectories. Our data extend the framework for the design of high-affinity Y220C mutant binders for use in personalized anticancer therapy and, more generally, highlight the importance of implementing protein dynamics and hydration patterns in the drug-discovery process.

Experimental and Theoretical Evaluation of the Ethynyl Moiety as a Halogen Bioisostere.

Bioisosteric replacements are widely used in medicinal chemistry to improve physicochemical and ADME properties of molecules while retaining or improving affinity. Here, using the p53 cancer mutant Y220C as a test case, we investigate both computationally and experimentally whether an ethynyl moiety is a suitable bioisostere to replace iodine in ligands that form halogen bonds with the protein backbone. This bioisosteric transformation is synthetically feasible via Sonogashira cross-coupling. In our test case of a particularly strong halogen bond, replacement of the iodine with an ethynyl group resulted in a 13-fold affinity loss. High-resolution crystal structures of the two analogues in complex with the p53-Y220C mutant enabled us to correlate the different affinities with particular features of the binding site and subtle changes in ligand binding mode. In addition, using QM calculations and analyzing the PDB, we provide general guidelines for identifying cases where such a transformation is likely to improve ligand recognition.

Tracing the evolution of the p53 tetramerization domain.

The tetrameric transcription factors p53, p63, and p73 evolved from a common ancestor and play key roles in tumor suppression and development. Surprisingly, p63 and p73 require a second helix in their tetramerization domain for the formation of stable tetramers that is absent in human p53, raising questions about the evolutionary processes leading to diversification. Here we determined the crystal structure of the zebrafish p53 tetramerization domain, which contains a second helix, reminiscent of p63 and p73, combined with p53-like features. Through comprehensive phylogenetic analyses, we systematically traced the evolution of vertebrate p53 family oligomerization domains back to the beginning of multicellular life. We provide evidence that their last common ancestor also had an extended p63/p73-like domain and pinpoint evolutionary events that shaped this domain during vertebrate radiation. Domain compaction and transformation of a structured into a flexible, intrinsically disordered region may have contributed to the expansion of the human p53 interactome.

Halogen-enriched fragment libraries as chemical probes for harnessing halogen bonding in fragment-based lead discovery.

Halogen bonding has recently experienced a renaissance, gaining increased recognition as a useful molecular interaction in the life sciences. Halogen bonds are favorable, fairly directional interactions between an electropositive region on the halogen (the σ-hole) and a number of different nucleophilic interaction partners. Some aspects of halogen bonding are not yet understood well enough to take full advantage of its potential in drug discovery. We describe and present the concept of halogen-enriched fragment libraries. These libraries consist of unique chemical probes, facilitating the identification of favorable halogen bonds by sharing the advantages of classical fragment-based screening. Besides providing insights into the nature and applicability of halogen bonding, halogen-enriched fragment libraries provide smart starting points for hit-to-lead evolution.

Structure of the SAS-6 cartwheel hub from Leishmania major.

Centrioles are cylindrical cell organelles with a ninefold symmetric peripheral microtubule array that is essential to template cilia and flagella. They are built around a central cartwheel assembly that is organized through homo-oligomerization of the centriolar protein SAS-6, but whether SAS-6 self-assembly can dictate cartwheel and thereby centriole symmetry is unclear. Here we show that Leishmania major SAS-6 crystallizes as a 9-fold symmetric cartwheel and provide the X-ray structure of this assembly at a resolution of 3.5 Å. We furthermore demonstrate that oligomerization of Leishmania SAS-6 can be inhibited by a small molecule in vitro and provide indications for its binding site. Our results firmly establish that SAS-6 can impose cartwheel symmetry on its own and indicate how this process might occur mechanistically in vivo. Importantly, our data also provide a proof-of-principle that inhibition of SAS-6 oligomerization by small molecules is feasible. DOI: http://dx.doi.org/10.7554/eLife.01812.001.

Machine learning estimates of natural product conformational energies.

Machine learning has been used for estimation of potential energy surfaces to speed up molecular dynamics simulations of small systems. We demonstrate that this approach is feasible for significantly larger, structurally complex molecules, taking the natural product Archazolid A, a potent inhibitor of vacuolar-type ATPase, from the myxobacterium Archangium gephyra as an example. Our model estimates energies of new conformations by exploiting information from previous calculations via Gaussian process regression. Predictive variance is used to assess whether a conformation is in the interpolation region, allowing a controlled trade-off between prediction accuracy and computational speed-up. For energies of relaxed conformations at the density functional level of theory (implicit solvent, DFT/BLYP-disp3/def2-TZVP), mean absolute errors of less than 1 kcal/mol were achieved. The study demonstrates that predictive machine learning models can be developed for structurally complex, pharmaceutically relevant compounds, potentially enabling considerable speed-ups in simulations of larger molecular structures.

Targeting histidine side chains in molecular design through nitrogen-halogen bonds.

Halogen bonds are directional noncovalent interactions that can be used to target electron donors in a protein binding site. In this study, we employ quantum chemical calculations to explore halogen···nitrogen contacts involving histidine side chains. We characterize the energetics on the MP2 level of theory using SCS-MP2 and CCSD(T)/CBS as reference calculations and elucidate their energy profile in suboptimal geometries. We derive simple rules allowing medicinal chemists and chemical biologists to easily determine preferred areas of interaction in a binding site and exploit them for scaffold decoration and design. Our work shows that nitrogen-halogen bonds are valuable interactions that are this far underexploited in patent applications, lead structure, and clinical candidate selection. We highlight their potential to increase binding affinities and suggest that they can significantly contribute to inducing and tuning subtype selectivities.

Small molecule induced reactivation of mutant p53 in cancer cells.

The p53 cancer mutant Y220C is an excellent paradigm for rescuing the function of conformationally unstable p53 mutants because it has a unique surface crevice that can be targeted by small-molecule stabilizers. Here, we have identified a compound, PK7088, which is active in vitro: PK7088 bound to the mutant with a dissociation constant of 140 µM and raised its melting temperature, and we have determined the binding mode of a close structural analogue by X-ray crystallography. We showed that PK7088 is biologically active in cancer cells carrying the Y220C mutant by a battery of tests. PK7088 increased the amount of folded mutant protein with wild-type conformation, as monitored by immunofluorescence, and restored its transcriptional functions. It induced p53-Y220C-dependent growth inhibition, cell-cycle arrest and apoptosis. Most notably, PK7088 increased the expression levels of p21 and the proapoptotic NOXA protein. PK7088 worked synergistically with Nutlin-3 on up-regulating p21 expression, whereas Nutlin-3 on its own had no effect, consistent with its mechanism of action. PK7088 also restored non-transcriptional apoptotic functions of p53 by triggering nuclear export of BAX to the mitochondria. We suggest a set of criteria for assigning activation of p53.

Principles and applications of halogen bonding in medicinal chemistry and chemical biology.

Halogen bonding has been known in material science for decades, but until recently, halogen bonds in protein-ligand interactions were largely the result of serendipitous discovery rather than rational design. In this Perspective, we provide insights into the phenomenon of halogen bonding, with special focus on its role in drug discovery. We summarize the theoretical background defining its strength and directionality, provide a systematic analysis of its occurrence and interaction geometries in protein-ligand complexes, and give recent examples where halogen bonding has been successfully harnessed for lead identification and optimization. In light of these data, we discuss the potential and limitations of exploiting halogen bonds for molecular recognition and rational drug design.