NMR Data-Driven Docking of HDM2-Inhibitor Complexes.

We present an automated NMR-guided docking workflow that can be used to generate models of protein-ligand complexes based on data from NOE NMR experiments. The first step is to generate a number of intermolecular distance constraints from experimental NOE data. Then, the ligand is docked on an ensemble of receptor structures to account for protein flexibility, and multiple poses are generated. Finally, we use the NOE-based constraints to filter and score docking poses based on the percentage of NOE constraints that are consistent with protein-ligand interatomic distances. This workflow was successfully used during a lead optimization project to generate models of synthetic protein-protein interaction (PPI) inhibitors bound to the HDM2 protein.

Discovery of a new series of PI3K-δ inhibitors from Virtual Screening.

PI3K-δ mediates key immune cell signaling pathways and is a target of interest for treatment of oncological and immunological disorders. Here we describe the discovery and optimization of a novel series of PI3K-δ selective inhibitors. We first identified hits containing an isoindolinone scaffold using a combined ligand- and receptor-based virtual screening workflow, and then improved potency and selectivity guided by structural data and modeling. Careful optimization of molecular properties led to compounds with improved permeability and pharmacokinetic profile, and high potency in a whole blood assay.

Discovery of IDO1 inhibitors containing a decahydroquinoline, decahydro-1,6-naphthyridine, or octahydro-1H-pyrrolo[3,2-c]pyridine scaffold.

A series of IDO1 inhibitors containing a decahydroquinoline, decahydro-1,6-naphthyridine, or octahydro-1H-pyrrolo[3,2-c]pyridine scaffold were identified with good cellular and human whole blood activity against IDO1. These inhibitors contain multiple chiral centers and all diastereomers were separated. The absolute stereochemistry of each isomers were not determined. Compounds 15 and 27 stood out as leads due to their good cellular as well as human whole blood IDO1 inhibition activity, low unbound clearance, and reasonable mean residence time in rat cassette PK studies.

SAR towards indoline and 3-azaindoline classes of IDO1 inhibitors.

A novel series of IDO1 inhibitors have been identified with good IDO1 Hela cell and human whole blood activity. These inhibitors contain an indoline or a 3-azaindoline scaffold. Their structure-activity-relationship studies have been explored. Compounds 37 and 41 stood out as leads due to their good potency in IDO1 Hela assay, good IDO1 unbound hWB IC50s, reasonable unbound clearance, and good MRT in rat and dog PK studies.

Discovery and optimization of heteroaryl piperazines as potent and selective PI3Kδ inhibitors.

A high-throughput screening (HTS) campaign identified a class of heteroaryl piperazines with excellent baseline affinity and selectivity for phosphoinositide 3-kinase δ (PI3Kδ) over closely related isoforms. Rapid evaluation and optimization of structure-activity relationships (SAR) for this class, leveraging the modular nature of this scaffold, facilitated development of this hit class into a series of potent and selective inhibitors of PI3Kδ. This effort culminated in the identification of 29, which displayed excellent potency in enzyme and cell-based assays, as well as favorable pharmacokinetic and off-target profiles.

Design of selective PI3Kδ inhibitors using an iterative scaffold-hopping workflow.

PI3Kδ mediates key immune cell signaling pathways and is a target of interest for multiple indications in immunology and oncology. Here we report a structure-based scaffold-hopping strategy for the design of chemically diverse PI3Kδ inhibitors. Using this strategy, we identified several scaffolds that can be combined to generate new PI3Kδ inhibitors with high potency and isoform selectivity. In particular, an oxindole-based scaffold was found to impart exquisite selectivity when combined with several hinge binding motifs.

Workflows and performances in the ranking prediction of 2016 D3R Grand Challenge 2: lessons learned from a collaborative effort.

The 2016 D3R Grand Challenge 2 includes both pose and affinity or ranking predictions. This article is focused exclusively on affinity predictions submitted to the D3R challenge from a collaborative effort of the modeling and informatics group. Our submissions include ranking of 102 ligands covering 4 different chemotypes against the FXR ligand binding domain structure, and the relative binding affinity predictions of the two designated free energy subsets of 15 and 18 compounds. Using all the complex structures prepared in the same way allowed us to cover many types of workflows and compare their performances effectively. We evaluated typical workflows used in our daily structure-based design modeling support, which include docking scores, force field-based scores, QM/MM, MMGBSA, MD-MMGBSA, and MacroModel interaction energy estimations. The best performing methods for the two free energy subsets are discussed. Our results suggest that affinity ranking still remains very challenging; that the knowledge of more structural information does not necessarily yield more accurate predictions; and that visual inspection and human intervention are considerably important for ranking. Knowledge of the mode of action and protein flexibility along with visualization tools that depict polar and hydrophobic maps are very useful for visual inspection. QM/MM-based workflows were found to be powerful in affinity ranking and are encouraged to be applied more often. The standardized input and output enable systematic analysis and support methodology development and improvement for high level blinded predictions.

Performance of multiple docking and refinement methods in the pose prediction D3R prospective Grand Challenge 2016.

We describe the performance of multiple pose prediction methods for the D3R 2016 Grand Challenge. The pose prediction challenge includes 36 ligands, which represent 4 chemotypes and some miscellaneous structures against the FXR ligand binding domain. In this study we use a mix of fully automated methods as well as human-guided methods with considerations of both the challenge data and publicly available data. The methods include ensemble docking, colony entropy pose prediction, target selection by molecular similarity, molecular dynamics guided pose refinement, and pose selection by visual inspection. We evaluated the success of our predictions by method, chemotype, and relevance of publicly available data. For the overall data set, ensemble docking, visual inspection, and molecular dynamics guided pose prediction performed the best with overall mean RMSDs of 2.4, 2.2, and 2.2 Å respectively. For several individual challenge molecules, the best performing method is evaluated in light of that particular ligand. We also describe the protein, ligand, and public information data preparations that are typical of our binding mode prediction workflow.

High-resolution crystal structures of factor XIa coagulation factor in complex with nonbasic high-affinity synthetic inhibitors.

Factor XI (FXI) is a key enzyme in the coagulation pathway and an attractive target for the development of anticoagulant drugs. A small number of high-resolution crystal structures of FXIa in complex with small synthetic inhibitors have been published to date. All of these ligands have a basic P1 group and bind exclusively in the nonprime side of the active site of FXIa. Here, two structures of FXIa in complex with nonbasic inhibitors that occupy both the prime and nonprime sides of the active site are presented. These new structures could be valuable in the design and optimization of new FXIa synthethic inhibitors.

Diminishing GSH-Adduct Formation of Tricyclic Diazepine-based Mutant IDH1 Inhibitors.

Mutant isocitrate dehydrogenase 1 (IDH1) has been identified as an attractive oncology target for which >70% of grade II and III gliomas and ∼10% of acute myeloid leukemia (AML) harbor somatic IDH1 mutations. These mutations confer a neomorphic gain of function, leading to the production of the oncometabolite (R)-2-hydroxyglutarate (2-HG). We identified and developed a potent, selective, and orally bioavailable brain-penetrant tricyclic diazepine scaffold that inhibits mutant IDH1. During the course of in vitro metabolism studies, GSH-adduct metabolites were observed. The hypothesis for GSH-adduct formation was driven by the electron-rich nature of the tricyclic core. Herein, we describe our efforts to reduce the electron-rich nature of the core. Ultimately, a strategy focused on core modifications to block metabolic hot spots coupled with substitution pattern changes (C8 N → C linked) led to the identification of new tricyclic analogues with minimal GSH-adduct formation across species while maintaining an overall balanced profile.

Oxetane Promise Delivered: Discovery of Long-Acting IDO1 Inhibitors Suitable for Q3W Oral or Parenteral Dosing.

3,3-Disubstituted oxetanes have been utilized as bioisosteres for gem-dimethyl and cyclobutane functionalities. We report the discovery of a novel class of oxetane indole-amine 2,3-dioxygenase (IDO1) inhibitors suitable for Q3W (once every 3 weeks) oral and parenteral dosing. A diamide class of IDO inhibitors was discovered through an automated ligand identification system (ALIS). Installation of an oxetane and fluorophenyl dramatically improved the potency. Identification of a biaryl moiety as an unconventional amide isostere addressed the metabolic liability of amide hydrolysis. Metabolism identification (Met-ID)-guided target design and the introduction of polarity resulted in the discovery of potent IDO inhibitors with excellent pharmacokinetic (PK) profiles in multiple species. To enable rapid synthesis of the key oxetane intermediate, a novel oxetane ring cyclization was also developed, as well as optimization of a literature route on kg scale. These IDO inhibitors may enable unambiguous proof-of-concept testing for the IDO1 inhibition mechanism for oncology.

1H-imidazo[4,5-c]pyridine-4-carbonitrile as cathepsin S inhibitors: separation of desired cellular activity from undesired tissue accumulation through optimization of basic nitrogen pka.

Based on the theoretical understanding of the in vivo lysosomotropism, by adjusting the pk(a) of basic nitrogen containing cathepsin S inhibitors, a set of compounds with pk(a) 6-8 were identified to have excellent cell based Lip10 activity, yet avoiding undesired sequestration in spleen.

Structural insights on ligand recognition at the human leukotriene B4 receptor 1.

The leukotriene B4 receptor 1 (BLT1) regulates the recruitment and chemotaxis of different cell types and plays a role in the pathophysiology of infectious, allergic, metabolic, and tumorigenic human diseases. Here we present a crystal structure of human BLT1 (hBLT1) in complex with a selective antagonist MK-D-046, developed for the treatment of type 2 diabetes and other inflammatory conditions. Comprehensive analysis of the structure and structure-activity relationship data, reinforced by site-directed mutagenesis and docking studies, reveals molecular determinants of ligand binding and selectivity toward different BLT receptor subtypes and across species. The structure helps to identify a putative membrane-buried ligand access channel as well as potential receptor binding modes of endogenous agonists. These structural insights of hBLT1 enrich our understanding of its ligand recognition and open up future avenues in structure-based drug design.

Identification of Potent Reverse Indazole Inhibitors for HPK1.

Hematopoietic progenitor kinase (HPK1), a negative regulator of TCR-mediated T-cell activation, has been recognized as a novel antitumor immunotherapy target. Structural optimization of kinase inhibitor 4 through a systematic two-dimensional diversity screen of pyrazolopyridines led to the identification of potent and selective compounds. Crystallographic studies with HPK1 revealed a favorable water-mediated interaction with Asp155 and a salt bridge to Asp101 with optimized heterocyclic solvent fronts that were critical for enhanced potency and selectivity. Computational studies of model systems revealed differences in torsional profiles that allowed for these beneficial protein-ligand interactions. Further optimization of molecular properties led to identification of potent and selective reverse indazole inhibitor 36 that inhibited phosphorylation of adaptor protein SLP76 in human PBMC and exhibited low clearance with notable bioavailability in in vivo rat studies.

Discovery of Diaminopyrimidine Carboxamide HPK1 Inhibitors as Preclinical Immunotherapy Tool Compounds.

Hematopoietic progenitor kinase 1 (HPK1), a serine/threonine kinase, is a negative immune regulator of T cell receptor (TCR) and B cell signaling that is primarily expressed in hematopoietic cells. Accordingly, it has been reported that HPK1 loss-of-function in HPK1 kinase-dead syngeneic mouse models shows enhanced T cell signaling and cytokine production as well as tumor growth inhibition in vivo, supporting its value as an immunotherapeutic target. Herein, we present the structurally enabled discovery of novel, potent, and selective diaminopyrimidine carboxamide HPK1 inhibitors. The key discovery of a carboxamide moiety was essential for enhanced enzyme inhibitory potency and kinome selectivity as well as sustained elevation of cellular IL-2 production across a titration range in human peripheral blood mononuclear cells. The elucidation of structure-activity relationships using various pendant amino ring systems allowed for the identification of several small molecule type-I inhibitors with promising in vitro profiles.

Utilization of Metabolite Identification and Structural Data to Guide Design of Low-Dose IDO1 Inhibitors.

Herein the discovery of potent IDO1 inhibitors with low predicted human dose is discussed. Metabolite identification (MetID) and structural data were used to strategically incorporate cyclopropane rings into this tetrahydronaphthyridine series of IDO1 inhibitors to improve their metabolic stability and potency. Enabling synthetic chemistry was developed to construct these unique fused cyclopropyl compounds, leading to inhibitors with improved pharmacokinetics and human whole blood potency and a predicted human oral dose as low as 9 mg once daily (QD).

Development of High-Throughput Assays for Evaluation of Hematopoietic Progenitor Kinase 1 Inhibitors.

Hematopoietic progenitor kinase 1 (HPK1), also referred to as mitogen-activated protein kinase kinase kinase kinase 1 (MAP4K1), is a serine/threonine kinase that negatively regulates T-cell signaling by phosphorylating Ser376 of Src homology 2 (SH2) domain-containing leukocyte protein of 76 kDa (SLP-76), a critical mediator of T-cell receptor activation. HPK1 loss of function mouse models demonstrated enhanced immune cell activation and beneficial antitumor activity. To enable discovery and functional characterization of high-affinity small-molecule HPK1 inhibitors, we have established high-throughput biochemical, cell-based, and novel pharmacodynamic (PD) assays. Kinase activity-based time-resolved fluorescence energy transfer (TR-FRET) assays were established as the primary biochemical approach to screen for potent inhibitors and assess selectivity against members of MAP4K and other closely related kinases. A proximal target engagement (TE) assay quantifying pSLP-76 levels as a readout and a distal assay measuring IL-2 secretion as a functional response were established using human peripheral blood mononuclear cells (PBMCs) from two healthy donors. Significant correlations between biochemical and cellular assays as well as excellent correlation between the two donors for the cellular assays were observed. pSLP-76 levels were further used as a PD marker in the preclinical murine model. This effort required the development of a novel ultrasensitive single-molecule array (SiMoA) assay to monitor pSLP-76 changes in mouse spleen.

Discovery of MK-4688: an Efficient Inhibitor of the HDM2-p53 Protein-Protein Interaction.

Identification of low-dose, low-molecular-weight, drug-like inhibitors of protein-protein interactions (PPIs) is a challenging area of research. Despite the challenges, the therapeutic potential of PPI inhibition has driven significant efforts toward this goal. Adding to recent success in this area, we describe herein our efforts to optimize a novel purine carboxylic acid-derived inhibitor of the HDM2-p53 PPI into a series of low-projected dose inhibitors with overall favorable pharmacokinetic and physical properties. Ultimately, a strategy focused on leveraging known binding hot spots coupled with biostructural information to guide the design of conformationally constrained analogs and a focus on efficiency metrics led to the discovery of MK-4688 (compound 56), a highly potent, selective, and low-molecular-weight inhibitor suitable for clinical investigation.

Projected Dose Optimization of Amino- and Hydroxypyrrolidine Purine PI3Kδ Immunomodulators.

The approvals of idelalisib and duvelisib have validated PI3Kδ inhibitors for the treatment for hematological malignancies driven by the PI3K/AKT pathway. Our program led to the identification of structurally distinct heterocycloalkyl purine inhibitors with excellent isoform and kinome selectivity; however, they had high projected human doses. Improved ligand contacts gave potency enhancements, while replacement of metabolic liabilities led to extended half-lives in preclinical species, affording PI3Kδ inhibitors with low once-daily predicted human doses. Treatment of C57BL/6-Foxp3-GDL reporter mice with 30 and 100 mg/kg/day of 3c (MSD-496486311) led to a 70% reduction in Foxp3-expressing regulatory T cells as observed through bioluminescence imaging with luciferin, consistent with the role of PI3K/AKT signaling in Treg cell proliferation. As a model for allergic rhinitis and asthma, treatment of ovalbumin-challenged Brown Norway rats with 0.3 to 30 mg/kg/day of 3c gave a dose-dependent reduction in pulmonary bronchoalveolar lavage inflammation eosinophil cell count.