Thiophene antibacterials that allosterically stabilize DNA-cleavage complexes with DNA gyrase.

A paucity of novel acting antibacterials is in development to treat the rising threat of antimicrobial resistance, particularly in Gram-negative hospital pathogens, which has led to renewed efforts in antibiotic drug discovery. Fluoroquinolones are broad-spectrum antibacterials that target DNA gyrase by stabilizing DNA-cleavage complexes, but their clinical utility has been compromised by resistance. We have identified a class of antibacterial thiophenes that target DNA gyrase with a unique mechanism of action and have activity against a range of bacterial pathogens, including strains resistant to fluoroquinolones. Although fluoroquinolones stabilize double-stranded DNA breaks, the antibacterial thiophenes stabilize gyrase-mediated DNA-cleavage complexes in either one DNA strand or both DNA strands. X-ray crystallography of DNA gyrase-DNA complexes shows the compounds binding to a protein pocket between the winged helix domain and topoisomerase-primase domain, remote from the DNA. Mutations of conserved residues around this pocket affect activity of the thiophene inhibitors, consistent with allosteric inhibition of DNA gyrase. This druggable pocket provides potentially complementary opportunities for targeting bacterial topoisomerases for antibiotic development.

Structure-guided design of antibacterials that allosterically inhibit DNA gyrase.

A series of DNA gyrase inhibitors were designed based on the X-ray structure of a parent thiophene scaffold with the objective to improve biochemical and whole-cell antibacterial activity, while reducing cardiac ion channel activity. The binding mode and overall design hypothesis of one series was confirmed with a co-crystal structure with DNA gyrase. Although some analogs retained both biochemical activity and whole-cell antibacterial activity, we were unable to significantly improve the activity of the series and analogs retained activity against the cardiac ion channels, therefore we stopped optimization efforts.

Synthesis and structure-activity relationships of imidazo[1,2-a]pyrimidin-5(1H)-ones as a novel series of beta isoform selective phosphatidylinositol 3-kinase inhibitors.

A series of PI3K-beta selective inhibitors, imidazo[1,2-a]-pyrimidin-5(1H)-ones, has been rationally designed based on the docking model of the more potent R enantiomer of TGX-221, identified by a chiral separation, in a PI3K-beta homology model. Synthesis and SAR of this novel chemotype are described. Several compounds in the series demonstrated potent growth inhibition in a PTEN-deficient breast cancer cell line MDA-MB-468 under anchorage independent conditions.

Synthesis and structure-activity relationships of 1,2,4-triazolo[1,5-a]pyrimidin-7(3H)-ones as novel series of potent ß isoform selective phosphatidylinositol 3-kinase inhibitors.

A series of 1,2,4-triazolo[1,5-a]pyrimidin-7(3H)-ones with excellent enzyme inhibition, improved isoform selectivity, and excellent inhibition of downstream phosphorylation of AKT has been identified. Several compounds in the series demonstrated potent (∼ 0.100 µM IC(50)) growth inhibition in a PTEN deficient cancer cell line.

Discovery of a First-in-Class Gut-Restricted RET Kinase Inhibitor as a Clinical Candidate for the Treatment of IBS.

Abdominal pain and abnormal bowel habits represent major symptoms for irritable bowel syndrome (IBS) patients that are not adequately managed. Although the etiology of IBS is not completely understood, many of the functions of the gastrointestinal (GI) tract are regulated by the enteric nervous system (ENS). Inflammation or stress-induced expression of growth factors or cytokines may lead to hyperinnervation of visceral afferent neurons in GI tract and contribute to the pathophysiology of IBS. Rearranged during transfection (RET) is a neuronal growth factor receptor tyrosine kinase critical for the development of the ENS as exemplified by Hirschsprung patients who carry RET loss-of-function mutations and lack normal colonic innervation leading to colonic obstruction. Similarly, RET signaling in the adult ENS maintains neuronal function by contributing to synaptic formation, signal transmission, and neuronal plasticity. Inhibition of RET in the ENS represents a novel therapeutic strategy for the normalization of neuronal function and the symptoms of IBS patients. Herein, we describe our screening effort and subsequent structure-activity relationships (SARs) in optimizing potency, selectivity, and mutagenicity of the series, which led to the discovery of a first-in-class, gut-restricted RET kinase inhibitor, 2-(4-(4-ethoxy-6-oxo-1,6-dihydropyridin-3-yl)-2-fluorophenyl)-N-(5-(1,1,1-trifluoro-2-methylpropan-2-yl)isoxazol-3-yl)acetamide (15, GSK3179106), as a clinical candidate for the treatment of IBS. GSK3179106 is a potent, selective, and gut-restricted pyridone hinge binder small molecule RET kinase inhibitor with a RET IC50 of 0.3 nM and is efficacious in vivo.

The role of quantum mechanics in structure-based drug design.

Herein we will focus on the use of quantum mechanics (QM) in drug design (DD) to solve disparate problems from scoring protein-ligand poses to building QM QSAR models. Through the variational principle of QM we know that we can obtain a more accurate representation of molecular systems than classical models, and while this is not a matter of debate, it still has not been shown that the expense of QM approaches is offset by improved accuracy in DD applications. Objectively validating the improved applicability and performance of QM over classical-based models in DD will be the focus of research in the coming years along with research on the conformational sampling problem as it relates to protein-ligand complexes.

A Perspective on Water Site Prediction Methods for Structure Based Drug Design.

Over the last decade, a number of computational methods have been developed, which attempt to evaluate the thermodynamic properties of individual water molecules at the solute-solvent interface, in order to assess contributions to protein-ligand binding. In some cases, these tools tell us what we already know, e.g. that hydrophobic pockets prefer lipophilic substituents, and in other cases the methods only seem to add clarity when retrospectively applied. Hence we have grappled with how to utilize such approaches to understand non-intuitive results and to generate chemistry ideas that otherwise would not have been developed. Here we provide our perspective on these methods and describe how results have been interpreted and applied. We include examples from GSK and elsewhere that highlight how water methods have been (1) utilized retrospectively to explain non-intuitive structure- activity relationships and (2) applied prospectively for chemistry design. Finally, we discuss where this field of study could lead to maximal impact in drug discovery research.

Quantum mechanics in structure-based drug design.

In principle, quantum mechanics provides a more accurate representation of molecular systems than other modeling approaches. While this notion is not a matter of dispute, it has not yet been definitively demonstrated within the realm of structure-based drug design that the use of quantum mechanical methods over the use of classical modeling approaches is justified in consideration of the increase in expense associated with quantum mechanical methods. Demonstrating that quantum mechanics-based methods can be superior to simpler models, and resolving problems relating to estimating the effects of conformational entropy, will provide key areas of interest in the coming years for in silico structure-based drug design. Recent applications using quantum mechanical methods in structure-based drug design are reviewed herein, and applications ranging from scoring receptor-ligand interactions using quantum mechanics to the generation of quantitative structure-activity relationships using quantum mechanics-derived descriptors are discussed.

Discovery of a Novel 2,6-Disubstituted Glucosamine Series of Potent and Selective Hexokinase 2 Inhibitors.

A novel series of potent and selective hexokinase 2 (HK2) inhibitors, 2,6-disubstituted glucosamines, has been identified based on HTS hits, exemplified by compound 1. Inhibitor-bound crystal structures revealed that the HK2 enzyme could adopt an "induced-fit" conformation. The SAR study led to the identification of potent HK2 inhibitors, such as compound 34 with greater than 100-fold selectivity over HK1. Compound 25 inhibits in situ glycolysis in a UM-UC-3 bladder tumor cell line via (13)CNMR measurement of [3-(13)C]lactate produced from [1,6-(13)C2]glucose added to the cell culture.

Large-scale validation of a quantum mechanics based scoring function: predicting the binding affinity and the binding mode of a diverse set of protein-ligand complexes.

Computational methods to calculate binding affinity in protein-ligand interaction are of immense interest because of obvious practical applications in structure-based drug design. Scoring functions attempt to calculate the variation in binding affinity of ligands-inhibitors bound to protein targets at various levels of theory. In this study we use semiempirical quantum mechanics to design a scoring function that can calculate the electrostatic interactions and solvation free energy expected during complexation. This physically based approach has the ability to capture binding affinity trends in a diverse range of protein-ligand complexes. We also show the predictive power of this scoring function within protein targets and its ability to score ligand poses docked to a protein target. We also demonstrate the ability of this scoring function to discriminate between native and decoy poses and highlight the crucial role played by electrostatic interactions in molecular recognition. Finally we compare the performance of our scoring function with other available scoring functions in the literature.

Discovery of a Potent Class of PI3Kα Inhibitors with Unique Binding Mode via Encoded Library Technology (ELT).

In the search of PI3K p110α wild type and H1047R mutant selective small molecule leads, an encoded library technology (ELT) campaign against the desired target proteins was performed which led to the discovery of a selective chemotype for PI3K isoforms from a three-cycle DNA encoded library. An X-ray crystal structure of a representative inhibitor from this chemotype demonstrated a unique binding mode in the p110α protein.

[3a,4]-Dihydropyrazolo[1,5a]pyrimidines: Novel, Potent, and Selective Phosphatidylinositol-3-kinase ß Inhibitors.

A series of novel [3a,4]dihydropyrazolo[1,5a]pyrimidines were identified, which were highly potent and selective inhibitors of PI3Kß. The template afforded the opportunity to develop novel SAR for both the hinge-binding (R3) and back-pocket (R4) substitutents. While cellular potency was relatively modest due to high protein binding, the series displayed low clearance in rat, mouse, and monkey.

Rational Design, Synthesis, and SAR of a Novel Thiazolopyrimidinone Series of Selective PI3K-beta Inhibitors.

A novel thiazolopyrimidinone series of PI3K-beta selective inhibitors has been identified. This chemotype has provided an excellent tool compound, 18, that showed potent growth inhibition in the PTEN-deficient breast cancer cell line MDA-MB-468 under anchorage-independent conditions, and it also demonstrated pharmacodynamic effects and efficacy in a PTEN-deficient prostate cancer PC-3 xenograft mouse model.

Discovery of GSK2126458, a Highly Potent Inhibitor of PI3K and the Mammalian Target of Rapamycin.

Phosphoinositide 3-kinase α (PI3Kα) is a critical regulator of cell growth and transformation, and its signaling pathway is the most commonly mutated pathway in human cancers. The mammalian target of rapamycin (mTOR), a class IV PI3K protein kinase, is also a central regulator of cell growth, and mTOR inhibitors are believed to augment the antiproliferative efficacy of PI3K/AKT pathway inhibition. 2,4-Difluoro-N-{2-(methyloxy)-5-[4-(4-pyridazinyl)-6-quinolinyl]-3-pyridinyl}benzenesulfonamide (GSK2126458, 1) has been identified as a highly potent, orally bioavailable inhibitor of PI3Kα and mTOR with in vivo activity in both pharmacodynamic and tumor growth efficacy models. Compound 1 is currently being evaluated in human clinical trials for the treatment of cancer.

Blind docking of pharmaceutically relevant compounds using RosettaLigand.

It is difficult to properly validate algorithms that dock a small molecule ligand into its protein receptor using data from the public domain: the predictions are not blind because the correct binding mode is already known, and public test cases may not be representative of compounds of interest such as drug leads. Here, we use private data from a real drug discovery program to carry out a blind evaluation of the RosettaLigand docking methodology and find that its performance is on average comparable with that of the best commercially available current small molecule docking programs. The strength of RosettaLigand is the use of the Rosetta sampling methodology to simultaneously optimize protein sidechain, protein backbone and ligand degrees of freedom; the extensive benchmark test described here identifies shortcomings in other aspects of the protocol and suggests clear routes to improving the method.

Pairwise decomposition of residue interaction energies using semiempirical quantum mechanical methods in studies of protein-ligand interaction.

Pairwise decomposition of the interaction energy between molecules is shown to be a powerful tool that can increase our understanding of macromolecular recognition processes. Herein we calculate the pairwise decomposition of the interaction energy between the protein human carbonic anhydrase II (HCAII) and the fluorine-substituted ligand N-(4-sulfamylbenzoyl)benzylamine (SBB) using semiempirical quantum mechanics based methods. We dissect the interaction between the ligand and the protein by dividing the ligand and the protein into subsystems to understand the structure-activity relationships as a result of fluorine substitution. In particular, the off-diagonal elements of the Fock matrix that is composed of the interaction between the ionic core and the valence electrons and the exchange energy between the subsystems or atoms of interest is examined in detail. Our analysis reveals that the fluorine-substituted benzylamine group of SBB does not directly affect the binding energy. Rather, we find that the strength of the interaction between Thr199 of HCAII and the sulfamylbenzoyl group of SBB affects the binding affinity between the protein and the ligand. These observations underline the importance of the sulfonamide group in binding affinity as shown by previous experiments (Maren, T. H.; Wiley: C. E. J. Med. Chem. 1968, 11, 228-232). Moreover, our calculations qualitatively agree with the structural aspects of these protein-ligand complexes as determined by X-ray crystallography.

A quantum mechanics-based scoring function: study of zinc ion-mediated ligand binding.

In this communication, we report the development of a novel quantum mechanics-based scoring function to predict free energy of ligand binding in the zinc metalloenzymes carbonic anhydrase (CA) and carboxypeptidase A (CPA). In particular, the AM1 method is used in conjunction with solvation modeling to predict the relative binding affinities of 18 CA and 5 CPA inhibitors. The effect of metal-ligand charge transfer is also discussed and shown to be different in CPA and CA, providing a further challenge to computing metalloenzyme binding affinities.

Pose scoring by NMR.

Recently, we have developed a fast approach to calculate NMR chemical shifts using the divide and conquer method at the semiempirical level. To demonstrate the utility of this approach for characterizing protein-ligand interactions, we used the deviation of calculated chemical shift perturbations from experiment to determine the orientation of a ligand (GPI-1046) in the binding pocket of the FK506 binding protein (FKBP12). Moreover, we were able to select the native state of the ligand from a collection of decoy poses. A key hydrogen bond between O1 and HN in Ile56 was also identified. Our results suggest that ligand-induced chemical shift perturbations can be used to refine protein/ligand structures.