Discovery and SAR of aryl hydroxy pyrimidinones as potent small molecule agonists of the GPCR APJ.

This article describes the discovery of aryl hydroxy pyrimidinones and the medicinal chemistry efforts to optimize this chemotype for potent APJ agonism. APJ is a G-protein coupled receptor whose natural agonist peptide, apelin, displays hemodynamic improvement in the cardiac function of heart failure patients. A high throughput screen was undertaken to identify small molecule hits that could be optimized to mimic the apelin in vitro response. A potent and low molecular weight aryl hydroxy pyrimidinone analog 30 was identified through optimization of an HTS hit and medicinal chemistry efforts to improve its properties.

VISTA is an acidic pH-selective ligand for PSGL-1.

Co-inhibitory immune receptors can contribute to T cell dysfunction in patients with cancer1,2. Blocking antibodies against cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) and programmed cell death 1 (PD-1) partially reverse this effect and are becoming standard of care in an increasing number of malignancies3. However, many of the other axes by which tumours become inhospitable to T cells are not fully understood. Here we report that V-domain immunoglobulin suppressor of T cell activation (VISTA) engages and suppresses T cells selectively at acidic pH such as that found in tumour microenvironments. Multiple histidine residues along the rim of the VISTA extracellular domain mediate binding to the adhesion and co-inhibitory receptor P-selectin glycoprotein ligand-1 (PSGL-1). Antibodies engineered to selectively bind and block this interaction in acidic environments were sufficient to reverse VISTA-mediated immune suppression in vivo. These findings identify a mechanism by which VISTA may engender resistance to anti-tumour immune responses, as well as an unexpectedly determinative role for pH in immune co-receptor engagement.

Evaluating a Multiscale Mechanistic Model of the Immune System to Predict Human Immunogenicity for a Biotherapeutic in Phase 1.

A mechanistic model of the immune response was evaluated for its ability to predict anti-drug antibody (ADA) and their impact on pharmacokinetics (PK) and pharmacodynamics (PD) for a biotherapeutic in a phase 1 clinical trial. Observed ADA incidence ranged from 33 to 67% after single doses and 27-50% after multiple doses. The model captured the single dose incidence well; however, there was overprediction after multiple dosing. The model was updated to include a T-regulatory (Treg) cell mediated tolerance, which reduced the overprediction (relative decrease in predicted incidence rate of 21.5-59.3% across multidose panels) without compromising the single dose predictions (relative decrease in predicted incidence rate of 0.6-13%). The Treg-adjusted model predicted no ADA impact on PK or PD, consistent with the observed data. A prospective phase 2 trial was simulated, including co-medication effects in the form of corticosteroid-induced immunosuppression. Predicted ADA incidences were 0-10%, depending on co-medication dosage. This work demonstrates the utility in applying an integrated, iterative modeling approach to predict ADA during different stages of clinical development.

Discovery of Clinical Candidate BMS-823778 as an Inhibitor of Human 11ß-Hydroxysteroid Dehydrogenase Type 1 (11ß-HSD-1).

BMS-823778 (2), a 1,2,4-triazolopyridinyl-methanol derived analog, was identified as a potent and selective inhibitor of human 11ß-hydroxysteroid dehydrogenase type 1 (11ß-HSD-1) enzyme (IC50 = 2.3 nM) with >10,000-fold selectivity over 11ß-HSD-2. Compound 2 exhibits robust acute pharmacodynamic effects in cynomolgus monkeys (ED50 = 0.6 mg/kg) and in diet-induced obese (DIO) mice (ED50 = 34 mg/kg). Compound 2 also showed excellent inhibition in an ex vivo adipose DIO mouse model (ED50 = 5.2 mg/kg). Oral bioavailability ranges from 44% to 100% in preclinical species. Its favorable development properties, pharmacokinetics, high adipose-to-plasma concentration ratio, and preclinical pharmacology profile have prompted the evaluation of 2 for the treatment of type 2 diabetes and metabolic syndrome in phase 2 clinical trials.

Discovery of Clinical Candidate 2-((2S,6S)-2-Phenyl-6-hydroxyadamantan-2-yl)-1-(3'-hydroxyazetidin-1-yl)ethanone [BMS-816336], an Orally Active Novel Selective 11ß-Hydroxysteroid Dehydrogenase Type 1 Inhibitor.

BMS-816336 (6n-2), a hydroxy-substituted adamantyl acetamide, has been identified as a novel, potent inhibitor against human 11ß-hydroxysteroid dehydrogenase type 1 (11ß-HSD1) enzyme (IC50 3.0 nM) with >10000-fold selectivity over human 11ß-hydroxysteroid dehydrogenase type 2 (11ß-HSD2). 6n-2 exhibits a robust acute pharmacodynamic effect in cynomolgus monkeys (ED50 0.12 mg/kg) and in DIO mice. It is orally bioavailable (%F ranges from 20 to 72% in preclinical species) and has a predicted pharmacokinetic profile of a high peak to trough ratio and short half-life in humans. This ADME profile met our selection criteria for once daily administration, targeting robust inhibition of 11ß-HSD1 enzyme for the first 12 h period after dosing followed by an "inhibition holiday" so that the potential for hypothalamic-pituitary-adrenal (HPA) axis activation might be mitigated. 6n-2 was found to be well-tolerated in phase 1 clinical studies and represents a potential new treatment for type 2 diabetes, metabolic syndrome, and other human diseases modulated by glucocorticoid control.

Discovery of pyridyl sulfonamide 11-beta-hydroxysteroid dehydrogenase type 1 (11ß-HSD1) inhibitors for the treatment of metabolic disorders.

A previous disclosure from this lab highlighted the discovery of pyridyl amides as potent 11ß-HSD1 inhibitors. In order to build additional novelty and polarity into this chemotype, replacement of the hydrogen-bonding carbonyl (CO) pharmacophore with the bioisosteric sulfonyl (SO2) group was examined. Despite initial comparisons suggesting the corresponding sulfonamides exhibited weaker activity versus their carbonyl counterparts, further optimization was performed in an effort to identify various potent and unique leads for the program. Judicious incorporation of polar moieties resulted in the identification of compounds with enhanced potency and lipophilicity profiles, resulting in leads with superior aqueous solubility and liver microsomal stability.

Optimization of 1,2,4-Triazolopyridines as Inhibitors of Human 11ß-Hydroxysteroid Dehydrogenase Type 1 (11ß-HSD-1).

Small alkyl groups and spirocyclic-aromatic rings directly attached to the left side and right side of the 1,2,4-triazolopyridines (TZP), respectively, were found to be potent and selective inhibitors of human 11ß-hydroxysteroid dehydrogenase-type 1 (11ß-HSD-1) enzyme. 3-(1-(4-Chlorophenyl)cyclopropyl)-8-cyclopropyl-[1,2,4]triazolo[4,3-a]pyridine (9f) was identified as a potent inhibitor of the 11ß-HSD-1 enzyme with reduced Pregnane-X receptor (PXR) transactivation activity. The binding orientation of this TZP series was revealed by X-ray crystallography structure studies.

Synthesis and structure-activity relationship of 2-adamantylmethyl tetrazoles as potent and selective inhibitors of human 11ß-hydroxysteroid dehydrogenase type 1 (11ß-HSD1).

A series of 2-adamantylmethyl tetrazoles bearing a quaternary carbon at the 2-position of the adamantane ring (i.e. structure A) have been designed and synthesized as novel, potent, and selective inhibitors of human 11ß-HSD1 enzyme. Based on the SAR and the docking experiment, we report for the first time a tetrazole moiety serving as the active pharmacophore for inhibitory activity of 11ß-HSD1 enzyme. Optimization of two regions of A, R(1) and R(2) respectively, was explored with a focus on improving the inhibitory activity (IC50) and the microsomal stability in both human and mouse species. These efforts led to the identification of 26, an orally bioavailable inhibitor of human 11ß-HSD1 with a favorable development profile.

Design, synthesis, and SAR studies of novel polycyclic acids as potent and selective inhibitors of human 11ß-hydroxysteroid dehydrogenase type 1 (11ß-HSD-1).

Starting from high throughput screening hit 2-adamantyl acetic acid 3, a series of polycyclic acids have been designed and synthesized as novel, potent, and selective inhibitors of human 11ß-HSD-1. Structure-activity relationships of two different regions of the chemotype (polycyclic ring and substituents on quaternary carbon) are discussed.

Discovery of 3-hydroxy-4-cyano-isoquinolines as novel, potent, and selective inhibitors of human 11ß-hydroxydehydrogenase 1 (11ß-HSD1).

Derived from the HTS hit 1, a series of hydroxyisoquinolines was discovered as potent and selective 11ß-HSD1 inhibitors with good cross species activity. Optimization of substituents at the 1 and 4 positions of the isoquinoline group in addition to the core modifications, with a special focus on enhancing metabolic stability and aqueous solubility, resulted in the identification of several compounds as potent advanced leads.

Closing the side-chain gap in protein loop modeling.

The success of structure-based drug design relies on accurate protein modeling where one of the key issues is the modeling and refinement of loops. This study takes a critical look at modeled loops, determining the effect of re-sampling side-chains after the loop conformation has been generated. The results are evaluated in terms of backbone and side-chain conformations with respect to the native loop. While models can contain loops with high quality backbone conformations, the side-chain orientations could be poor, and therefore unsuitable for ligand docking and structure-based design. In this study, we report on the ability to model loop side-chains accurately using a variety of commercially available algorithms that include rotamer libraries, systematic torsion scans and knowledge-based methods.

Engineering enzymes for improved catalytic efficiency: a computational study of site mutagenesis in epothilone-B hydroxylase.

Epothilone F, 21-hydroxyl-epothilone B, is an intermediate in the synthesis of BMS-310705, an antitumor compound that has been evaluated in Phase I clinical trials. A bioconversion process utilizing the Gram-positive bacterium Amycolatopsis orientalis was used to prepare epothilone F from epothilone B. In order to improve the yield of epothilone F, a mutagenesis program was performed with the goal of engineering the epothilone-B hydroxylase (EBH) enzyme to improve the yield of epothilone F through oxidative biotransformation. The mutations in EBH increased the yield of epothilone F from 21% in the recombinant expression system to higher than 80% utilizing the best EBH mutants. The studies described here show how a homology model of EBH was used to obtain an understanding of the possible mechanism that led to improved yield of epothilone F in the mutated enzymes. A novel aspect of this study is that it provides some insight into how mutations distant from the binding site can affect enzyme activity.

Pyridine amides as potent and selective inhibitors of 11beta-hydroxysteroid dehydrogenase type 1.

Several series of pyridine amides were identified as selective and potent 11beta-HSD1 inhibitors. The most potent inhibitors feature 2,6- or 3,5-disubstitution on the pyridine core. Various linkers (CH(2)SO(2), CH(2)S, CH(2)O, S, O, N, bond) between the distal aryl and central pyridyl groups are tolerated, and lipophilic amide groups are generally favored. On the distal aryl group, a number of substitutions are well tolerated. A crystal structure was obtained for a complex between 11beta-HSD1 and the most potent inhibitor in this series.

Cofactor-specific modulation of 11beta-hydroxysteroid dehydrogenase 1 inhibitor potency.

11beta-hydroxysteroid dehydrogenase 1 regulates the tissue availability of cortisol by interconverting cortisone and cortisol. It is capable of functioning as both a reductase and a dehydrogenase depending upon the surrounding milieu. In this work, we have studied the reaction mechanism of a soluble form of human 11beta-hydroxysteroid dehydrogenase 1 and its mode of inhibition by potent and selective inhibitors belonging to three different structural classes. We found that catalysis follows an ordered addition with NADP(H) binding preceding the binding of the steroid. While all three inhibitors tested bound to the steroid binding pocket, they differed in their interactions with the cofactor NADP(H). Compound A, a pyridyl amide bound more efficiently to the NADPH-bound form of 11beta-hydroxysteroid dehydrogenase 1. Compound B, an adamantyl triazole, was unaffected by NADP(H) binding and the sulfonamide, Compound C, showed preferential binding to the NADP+ -bound form of 11beta-hydroxysteroid dehydrogenase 1. These differences were found to augment significant selectivity towards inhibition of the reductase reaction versus the dehydrogenase reaction. This selectivity may translate to differences in the in vivo effects of 11beta-hydroxysteroid dehydrogenase 1 inhibitors.

Loopholes and missing links in protein modeling.

This paper provides an unbiased comparison of four commercially available programs for loop sampling, Prime, Modeler, ICM, and Sybyl, each of which uses a different modeling protocol. The study assesses the quality of results and examines the relative strengths and weaknesses of each method. The set of loops to be modeled varied in length from 4-12 amino acids. The approaches used for loop modeling can be classified into two methodologies: ab initio loop generation (Modeler and Prime) and database searches (Sybyl and ICM). Comparison of the modeled loops to the native structures was used to determine the accuracy of each method. All of the protocols returned similar results for short loop lengths (four to six residues), but as loop length increased, the quality of the results varied among the programs. Prime generated loops with RMSDs <2.5 A for loops up to 10 residues, while the other three methods met the 2.5 A criteria at seven-residue loops. Additionally, the ability of the software to utilize disulfide bonds and X-ray crystal packing influenced the quality of the results. In the final analysis, the top-ranking loop from each program was rarely the loop with the lowest RMSD with respect to the native template, revealing a weakness in all programs to correctly rank the modeled loops.

A comparative study of available software for high-accuracy homology modeling: from sequence alignments to structural models.

An open question in protein homology modeling is, how well do current modeling packages satisfy the dual criteria of quality of results and practical ease of use? To address this question objectively, we examined homology-built models of a variety of therapeutically relevant proteins. The sequence identities across these proteins range from 19% to 76%. A novel metric, the difference alignment index (DAI), is developed to aid in quantifying the quality of local sequence alignments. The DAI is also used to construct the relative sequence alignment (RSA), a new representation of global sequence alignment that facilitates comparison of sequence alignments from different methods. Comparisons of the sequence alignments in terms of the RSA and alignment methodologies are made to better understand the advantages and caveats of each method. All sequence alignments and corresponding 3D models are compared to their respective structure-based alignments and crystal structures. A variety of protein modeling software was used. We find that at sequence identities >40%, all packages give similar (and satisfactory) results; at lower sequence identities (<25%), the sequence alignments generated by Profit and Prime, which incorporate structural information in their sequence alignment, stand out from the rest. Moreover, the model generated by Prime in this low sequence identity region is noted to be superior to the rest. Additionally, we note that DSModeler and MOE, which generate reasonable models for sequence identities >25%, are significantly more functional and easier to use when compared with the other structure-building software.

A 3D-QSAR model for CYP2D6 inhibition in the aryloxypropanolamine series.

A comparative molecular similarity index analysis (CoMSiA) has been performed for cytochrome P450 2D6 inhibition on a series of aryloxypropanolamines to determine the factors contributing to this activity. The model is in agreement with a CYP2D6 homology model constructed on the basis of the mammalian CYP2C5 crystal structure. The energy minimized conformations were generated using the systematic search methodology in Sybyl 6.7. The model not only elucidated the relationship between structure and biological activity but, more importantly, provided useful strategies to modulate CYP2D6 affinity in the aryloxypropanolamine series.

An assessment of protein-ligand binding site polarizability.

Electronic polarizability, an important physical property of biomolecules, is currently ignored in most biomolecular calculations. Yet, it is widely believed that polarization could account for a substantial fraction of the total nonbonded energy of a system. This belief is supported by studies of small complexes in vacuum. This perception is driving the development of a new class of polarizable force fields for biomolecular calculations. However, the quantification of this term for protein-ligand complexes has never been attempted. Here we explore the polarizable nature of protein-ligand complexes in order to evaluate the importance of this effect. We introduce two indexes describing the polarizability of protein binding sites. These we apply to a large range of pharmaceutically relevant complexes. We offer a recommendation of particular complexes as test systems with which to determine the effects of polarizability and as test cases with which to test the new generation of force fields. Additionally, we provide a tabulation of the amino acid composition of these binding sites and show that composition can be specific for certain classes of proteins. We also show that the relative abundance of some amino acids is different in binding sites than elsewhere in a protein's structure.