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.

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.

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.

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.

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.

Kinetic mechanism of p53 oncogenic mutant aggregation and its inhibition.

Aggregation of destabilized mutants of the tumor suppressor p53 is a major route for its loss of activity. In order to assay drugs that inhibit aggregation of p53, we established the basic kinetics of aggregation of its core domain, using the mutant Y220C that has a mutation-induced, druggable cavity. Aggregation monitored by light scattering followed lag kinetics. Electron microscopy revealed the formation of small aggregates that subsequently grew to larger amorphous aggregates. The kinetics of aggregation produced surprising results: progress curves followed either by the binding of Thioflavin T or the fluorescence of the protein at 340 nm fitted well to simple two-step sequential first-order lag kinetics with rate constants k(1) and k(2) that were independent of protein concentration, and not to classical nucleation-growth. We suggest a mechanism of first-order formation of an aggregation competent state as being rate determining followed by rapid polymerization with the higher order kinetics. By measuring the inhibition kinetics of k(1) and k(2), we resolved that the process with the higher rate constant followed that of the lower. Further, there was only partial inhibition of k(1) and k(2), which showed two parallel pathways of aggregation, one via a state that requires unfolding of the protein and the other of partial unfolding with the ligand still bound. Inhibition kinetics of ligands provides a useful tool for probing an aggregation mechanism.

Lithocholic acid is an endogenous inhibitor of MDM4 and MDM2.

The proteins MDM2 and MDM4 are key negative regulators of the tumor suppressor protein p53, which are frequently upregulated in cancer cells. They inhibit the transactivation activity of p53 by binding separately or in concert to its transactivation domain. MDM2 is also a ubiquitin ligase that leads to the degradation of p53. Accordingly, MDM2 and MDM4 are important targets for drugs to inhibit their binding to p53. We found from in silico screening and confirmed by experiment that lithocholic acid (LCA) binds to the p53 binding sites of both MDM2 and MDM4 with a fivefold preference for MDM4. LCA is an endogenous steroidal bile acid, variously reported to have both carcinogenic and apoptotic activities. The comparison of LCA effects on apoptosis in HCT116 p53(+/+) vs. p53(-/-) cells shows a predominantly p53-mediated induction of caspase-3/7. The dissociation constants are in the µM region, but only modest inhibition of binding of MDM2 and MDM4 is required to negate their upregulation because they have to compete with transcriptional coactivator p300 for binding to p53. Binding was weakened by structural changes in LCA, and so it may be a natural ligand of MDM2 and MDM4, raising the possibility that MDM proteins may be sensors for specific steroids.

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.

Halogen-enriched fragment libraries as leads for drug rescue of mutant p53.

The destabilizing p53 cancer mutation Y220C creates a druggable surface crevice. We developed a strategy exploiting halogen bonding for lead discovery to stabilize the mutant with small molecules. We designed halogen-enriched fragment libraries (HEFLibs) as starting points to complement classical approaches. From screening of HEFLibs and subsequent structure-guided design, we developed substituted 2-(aminomethyl)-4-ethynyl-6-iodophenols as p53-Y220C stabilizers. Crystal structures of their complexes highlight two key features: (i) a central scaffold with a robust binding mode anchored by halogen bonding of an iodine with a main-chain carbonyl and (ii) an acetylene linker, enabling the targeting of an additional subsite in the crevice. The best binders showed induction of apoptosis in a human cancer cell line with homozygous Y220C mutation. Our structural and biophysical data suggest a more widespread applicability of HEFLibs in drug discovery.

Using halogen bonds to address the protein backbone: a systematic evaluation.

Halogen bonds are specific embodiments of the sigma hole bonding paradigm. They represent directional interactions between the halogens chlorine, bromine, or iodine and an electron donor as binding partner. Using quantum chemical calculations at the MP2 level, we systematically explore how they can be used in molecular design to address the omnipresent carbonyls of the protein backbone. We characterize energetics and directionality and elucidate their spatial variability in sub-optimal geometries that are expected to occur in protein-ligand complexes featuring a multitude of concomitant interactions. By deriving simple rules, we aid medicinal chemists and chemical biologists in easily exploiting them for scaffold decoration and design. Our work shows that carbonyl-halogen bonds may be used to expand the patentable medicinal chemistry space, redefining halogens as key features. Furthermore, this data will be useful for implementing halogen bonds into pharmacophore models or scoring functions making the QM information available for automatic molecular recognition in virtual high throughput screening.

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.

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 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.

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.

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.

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.

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.

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.