The discovery and optimization of naphthalene-linked P2-P4 Macrocycles as inhibitors of HCV NS3 protease.

Naphthalene-linked P2-P4 macrocycles within a tri-peptide-based acyl sulfonamide chemotype have been synthesized and found to inhibit HCV NS3 proteases representing genotypes 1a and 1b with single digit nanomolar potency. The pharmacokinetic profile of compounds in this series was optimized through structural modifications along the macrocycle tether as well as the P1 subsite. Ultimately a compound with oral bioavailability of 100% in rat, and a long half-life in plasma was obtained. However, compounds in this macrocyclic series exhibited cardiac effects in an isolated rabbit heart model and for this reason further optimization efforts were discontinued.

Structure-activity relationships of 4-hydroxy-4-biaryl-proline acylsulfonamide tripeptides: A series of potent NS3 protease inhibitors for the treatment of hepatitis C virus.

The design and synthesis of a series of tripeptide acylsulfonamides as potent inhibitors of the HCV NS3/4A serine protease is described. These analogues house a C4 aryl, C4 hydroxy-proline at the S2 position of the tripeptide scaffold. Information relating to structure-activity relationships as well as the pharmacokinetic and cardiovascular profiles of these analogues is provided.

Ligand placement based on prior structures: the guided ligand-replacement method.

The process of iterative structure-based drug design involves the X-ray crystal structure determination of upwards of 100 ligands with the same general scaffold (i.e. chemotype) complexed with very similar, if not identical, protein targets. In conjunction with insights from computational models and assays, this collection of crystal structures is analyzed to improve potency, to achieve better selectivity and to reduce liabilities such as absorption, distribution, metabolism, excretion and toxicology. Current methods for modeling ligands into electron-density maps typically do not utilize information on how similar ligands bound in related structures. Even if the electron density is of sufficient quality and resolution to allow de novo placement, the process can take considerable time as the size, complexity and torsional degrees of freedom of the ligands increase. A new module, Guided Ligand Replacement (GLR), was developed in Phenix to increase the ease and success rate of ligand placement when prior protein-ligand complexes are available. At the heart of GLR is an algorithm based on graph theory that associates atoms in the target ligand with analogous atoms in the reference ligand. Based on this correspondence, a set of coordinates is generated for the target ligand. GLR is especially useful in two situations: (i) modeling a series of large, flexible, complicated or macrocyclic ligands in successive structures and (ii) modeling ligands as part of a refinement pipeline that can automatically select a reference structure. Even in those cases for which no reference structure is available, if there are multiple copies of the bound ligand per asymmetric unit GLR offers an efficient way to complete the model after the first ligand has been placed. In all of these applications, GLR leverages prior knowledge from earlier structures to facilitate ligand placement in the current structure.

Automating crystallographic structure solution and refinement of protein-ligand complexes.

High-throughput drug-discovery and mechanistic studies often require the determination of multiple related crystal structures that only differ in the bound ligands, point mutations in the protein sequence and minor conformational changes. If performed manually, solution and refinement requires extensive repetition of the same tasks for each structure. To accelerate this process and minimize manual effort, a pipeline encompassing all stages of ligand building and refinement, starting from integrated and scaled diffraction intensities, has been implemented in Phenix. The resulting system is able to successfully solve and refine large collections of structures in parallel without extensive user intervention prior to the final stages of model completion and validation.

Design, synthesis, functional and structural characterization of an inhibitor of N-acetylneuraminate-9-phosphate phosphatase: observation of extensive dynamics in an enzyme/inhibitor complex.

The design, synthesis and characterization of a phosphonate inhibitor of N-acetylneuraminate-9-phosphate phosphatase (HDHD4) is described. Compound 3, where the substrate C-9 oxygen was replaced with a nonlabile CH2 group, inhibits HDHD4 with a binding affinity (IC50 11µM) in the range of the native substrate Neu5Ac-9-P (compound 1, Km 47µM). Combined SAR, modeling and NMR studies are consistent with the phosphonate group in inhibitor 3 forming a stable complex with native Mg(2+). In addition to this key interaction, the C-1 carboxylate of the sugar interacts with a cluster of basic residues, K141, R104 and R72. Comparative NMR studies of compounds 3 and 1 with Ca(2+) and Mg(2+) are indicative of a highly dynamic process in the active site for the HDHD4/Mg(2+)/3 complex. Possible explanations for this observation are discussed.

Structural basis for CARM1 inhibition by indole and pyrazole inhibitors.

CARM1 (co-activator-associated arginine methyltransferase 1) is a PRMT (protein arginine N-methyltransferase) family member that catalyses the transfer of methyl groups from SAM (S-adenosylmethionine) to the side chain of specific arginine residues of substrate proteins. This post-translational modification of proteins regulates a variety of transcriptional events and other cellular processes. Moreover, CARM1 is a potential oncological target due to its multiple roles in transcription activation by nuclear hormone receptors and other transcription factors such as p53. Here, we present crystal structures of the CARM1 catalytic domain in complex with cofactors [SAH (S-adenosyl-L-homocysteine) or SNF (sinefungin)] and indole or pyazole inhibitors. Analysis of the structures reveals that the inhibitors bind in the arginine-binding cavity and the surrounding pocket that exists at the interface between the N- and C-terminal domains. In addition, we show using ITC (isothermal titration calorimetry) that the inhibitors bind to the CARM1 catalytic domain only in the presence of the cofactor SAH. Furthermore, sequence differences for select residues that interact with the inhibitors may be responsible for the CARM1 selectivity against PRMT1 and PRMT3. Together, the structural and biophysical information should aid in the design of both potent and specific inhibitors of CARM1.

7-Oxopyrrolopyridine-derived DPP4 inhibitors-mitigation of CYP and hERG liabilities via introduction of polar functionalities in the active site.

Design, synthesis, and SAR of 7-oxopyrrolopyridine-derived DPP4 inhibitors are described. The preferred stereochemistry of these atropisomeric biaryl analogs has been identified as Sa. Compound (+)-3t, with a K(i) against DPP4, DPP8, and DPP9 of 0.37 nM, 2.2, and 5.7 µM, respectively, showed a significant improvement in insulin response after single doses of 3 and 10 µmol/kg in ob/ob mice.

Arylsulfonamidopiperidone derivatives as a novel class of factor Xa inhibitors.

The design, synthesis and SAR of a novel class of valerolactam-based arylsulfonamides as potent and selective FXa inhibitors is reported. The arylsulfonamide-valerolactam scaffold was derived based on the proposed bioisosterism to the arylcyanoguanidine-caprolactam core in known FXa inhibitors. The SAR study led to compound 46 as the most potent FXa inhibitor in this series, with an IC(50) of 7 nM and EC(2×PT) of 1.7 µM. The X-ray structure of compound 40 bound to FXa shows that the sulfonamide-valerolactam scaffold anchors the aryl group in the S1 and the novel acylcytisine pharmacophore in the S4 pockets.

Cyanoguanidine-based lactam derivatives as a novel class of orally bioavailable factor Xa inhibitors.

The N,N'-disubstituted cyanoguanidine is an excellent bioisostere of the thiourea and ketene aminal functional groups. We report the design and synthesis of a novel class of cyanoguanidine-based lactam derivatives as potent and orally active FXa inhibitors. The SAR studies led to the discovery of compound 4 (BMS-269223, K(i)=6.5nM, EC(2xPT)=32muM) as a selective, orally bioavailable FXa inhibitor with an excellent in vitro liability profile, favorable pharmacokinetics and pharmacodynamics in animal models. The X-ray crystal structure of 4 bound in FXa is presented and key ligand-protein interactions are discussed.

Aroylguanidine-based factor Xa inhibitors: the discovery of BMS-344577.

We report the design and synthesis of a novel class of N,N'-disubstituted aroylguanidine-based lactam derivatives as potent and orally active FXa inhibitors. The structure-activity relationships (SAR) investigation led to the discovery of the nicotinoyl guanidine 22 as a potent FXa inhibitor (FXa IC(50)=4 nM, EC(2xPT)=7 microM). However, the potent CYP3A4 inhibition activity (IC(50)=0.3 microM) of 22 precluded its further development. Detailed analysis of the X-ray crystal structure of compound 22 bound to FXa indicated that the substituent at the 6-position of the nicotinoyl group of 22 would be solvent-exposed, suggesting that efforts to attenuate the unwanted CYP activity could focus at this position without affecting FXa potency significantly. Further SAR studies on the 6-substituted nicotinoyl guanidines resulted in the discovery of 6-(dimethylcarbamoyl) nicotinoyl guanidine 36 (BMS-344577, IC(50)=9 nM, EC(2xPT)=2.5 microM), which was found to be a selective, orally efficacious FXa inhibitor with an excellent in vitro liability profile, favorable pharmacokinetics and pharmacodynamics in animal models.

Involvement of DPP-IV catalytic residues in enzyme-saxagliptin complex formation.

The inhibition of DPP-IV by saxagliptin has been proposed to occur through formation of a covalent but reversible complex. To evaluate further the mechanism of inhibition, we determined the X-ray crystal structure of the DPP-IV:saxagliptin complex. This structure reveals covalent attachment between S630 and the inhibitor nitrile carbon (C-O distance <1.3 A). To investigate whether this serine addition is assisted by the catalytic His-Asp dyad, we generated two mutants of DPP-IV, S630A and H740Q, and assayed them for ability to bind inhibitor. DPP-IV H740Q bound saxagliptin with an approximately 1000-fold reduction in affinity relative to DPP-IV WT, while DPP-IV S630A showed no evidence for binding inhibitor. An analog of saxagliptin lacking the nitrile group showed unchanged binding properties to the both mutant proteins, highlighting the essential role S630 and H740 play in covalent bond formation between S630 and saxagliptin. Further supporting mechanism-based inhibition by saxagliptin, NMR spectra of enzyme-saxagliptin complexes revealed the presence of three downfield resonances with low fractionation factors characteristic of short and strong hydrogen bonds (SSHB). Comparison of the NMR spectra of various wild-type and mutant DPP-IV:ligand complexes enabled assignment of a resonance at approximately 14 ppm to H740. Two additional DPP-IV mutants, Y547F and Y547Q, generated to probe potential stabilization of the enzyme-inhibitor complex by this residue, did not show any differences in inhibitor binding either by ITC or NMR. Together with the previously published enzymatic data, the structural and binding data presented here strongly support a histidine-assisted covalent bond formation between S630 hydroxyl oxygen and the nitrile group of saxagliptin.

The structure of Dasatinib (BMS-354825) bound to activated ABL kinase domain elucidates its inhibitory activity against imatinib-resistant ABL mutants.

Chronic myeloid leukemia (CML) is caused by the constitutively activated tyrosine kinase breakpoint cluster (BCR)-ABL. Current frontline therapy for CML is imatinib, an inhibitor of BCR-ABL. Although imatinib has a high rate of clinical success in early phase CML, treatment resistance is problematic, particularly in later stages of the disease, and is frequently mediated by mutations in BCR-ABL. Dasatinib (BMS-354825) is a multitargeted tyrosine kinase inhibitor that targets oncogenic pathways and is a more potent inhibitor than imatinib against wild-type BCR-ABL. It has also shown preclinical activity against all but one of the imatinib-resistant BCR-ABL mutants tested to date. Analysis of the crystal structure of dasatinib-bound ABL kinase suggests that the increased binding affinity of dasatinib over imatinib is at least partially due to its ability to recognize multiple states of BCR-ABL. The structure also provides an explanation for the activity of dasatinib against imatinib-resistant BCR-ABL mutants.

Interactive comparison and remediation of collections of macromolecular structures.

Often similar structures need to be compared to reveal local differences throughout the entire model or between related copies within the model. Therefore, a program to compare multiple structures and enable correction any differences not supported by the density map was written within the Phenix framework (Adams et al., Acta Cryst 2010; D66:213-221). This program, called Structure Comparison, can also be used for structures with multiple copies of the same protein chain in the asymmetric unit, that is, as a result of non-crystallographic symmetry (NCS). Structure Comparison was designed to interface with Coot(Emsley et al., Acta Cryst 2010; D66:486-501) and PyMOL(DeLano, PyMOL 0.99; 2002) to facilitate comparison of large numbers of related structures. Structure Comparison analyzes collections of protein structures using several metrics, such as the rotamer conformation of equivalent residues, displays the results in tabular form and allows superimposed protein chains and density maps to be quickly inspected and edited (via the tools in Coot) for consistency, completeness and correctness.

Potent Inhibitors of Hepatitis C Virus NS3 Protease: Employment of a Difluoromethyl Group as a Hydrogen-Bond Donor.

The design and synthesis of potent, tripeptidic acylsulfonamide inhibitors of HCV NS3 protease that contain a difluoromethyl cyclopropyl amino acid at P1 are described. A cocrystal structure of 18 with a NS3/4A protease complex suggests the presence of a H-bond between the polarized C-H of the CHF2 moiety and the backbone carbonyl of Leu135 of the enzyme. Structure-activity relationship studies indicate that this H-bond enhances enzyme inhibitory potency by 13- and 17-fold compared to the CH3 and CF3 analogues, respectively, providing insight into the deployment of this unique amino acid.

Discovery of a Potent Acyclic, Tripeptidic, Acyl Sulfonamide Inhibitor of Hepatitis C Virus NS3 Protease as a Back-up to Asunaprevir with the Potential for Once-Daily Dosing.

The discovery of a back-up to the hepatitis C virus NS3 protease inhibitor asunaprevir (2) is described. The objective of this work was the identification of a drug with antiviral properties and toxicology parameters similar to 2, but with a preclinical pharmacokinetic (PK) profile that was predictive of once-daily dosing. Critical to this discovery process was the employment of an ex vivo cardiovascular (CV) model which served to identify compounds that, like 2, were free of the CV liabilities that resulted in the discontinuation of BMS-605339 (1) from clinical trials. Structure-activity relationships (SARs) at each of the structural subsites in 2 were explored with substantial improvement in PK through modifications at the P1 site, while potency gains were found with small, but rationally designed structural changes to P4. Additional modifications at P3 were required to optimize the CV profile, and these combined SARs led to the discovery of BMS-890068 (29).

Discovery of N-(2-chloro-6-methyl- phenyl)-2-(6-(4-(2-hydroxyethyl)- piperazin-1-yl)-2-methylpyrimidin-4- ylamino)thiazole-5-carboxamide (BMS-354825), a dual Src/Abl kinase inhibitor with potent antitumor activity in preclinical assays.

A series of substituted 2-(aminopyridyl)- and 2-(aminopyrimidinyl)thiazole-5-carboxamides was identified as potent Src/Abl kinase inhibitors with excellent antiproliferative activity against hematological and solid tumor cell lines. Compound 13 was orally active in a K562 xenograft model of chronic myelogenous leukemia (CML), demonstrating complete tumor regressions and low toxicity at multiple dose levels. On the basis of its robust in vivo activity and favorable pharmacokinetic profile, 13 was selected for additional characterization for oncology indications.

Dynamically optimizing experiment schedules of a laboratory robot system with simulated annealing.

A scheduler has been developed for an integrated laboratory robot system that operates in an always-on mode. The integrated system is designed for imaging plates containing protein crystallization experiments, and it allows crystallographers to enter plates at any time and request that they be imaged at multiple time points in the future. The scheduler must rearrange tasks within the time it takes to image one plate, trading off the quality of the schedule for the speed of the computation. For this reason, the scheduler was based on a simulated annealing algorithm with an objective function that makes use of a linear programming solver. To optimize the scheduler, extensive computational simulations were performed involving a difficult but representative scheduling problem. The simulations explore multiple configurations of the simulated annealing algorithm, including both geometric and adaptive annealing schedules, 3 neighborhood functions, and 20 neighborhood diameters. An optimal configuration was found that produced the best results in less than 60 seconds, well within the window necessary to dynamically reschedule imaging tasks as new plates are entered into the system.

The discovery of asunaprevir (BMS-650032), an orally efficacious NS3 protease inhibitor for the treatment of hepatitis C virus infection.

The discovery of asunaprevir (BMS-650032, 24) is described. This tripeptidic acylsulfonamide inhibitor of the NS3/4A enzyme is currently in phase III clinical trials for the treatment of hepatitis C virus infection. The discovery of 24 was enabled by employing an isolated rabbit heart model to screen for the cardiovascular (CV) liabilities (changes to HR and SNRT) that were responsible for the discontinuation of an earlier lead from this chemical series, BMS-605339 (1), from clinical trials. The structure-activity relationships (SARs) developed with respect to CV effects established that small structural changes to the P2* subsite of the molecule had a significant impact on the CV profile of a given compound. The antiviral activity, preclincial PK profile, and toxicology studies in rat and dog supported clinical development of BMS-650032 (24).

Discovery and early clinical evaluation of BMS-605339, a potent and orally efficacious tripeptidic acylsulfonamide NS3 protease inhibitor for the treatment of hepatitis C virus infection.

The discovery of BMS-605339 (35), a tripeptidic inhibitor of the NS3/4A enzyme, is described. This compound incorporates a cyclopropylacylsulfonamide moiety that was designed to improve the potency of carboxylic acid prototypes through the introduction of favorable nonbonding interactions within the S1' site of the protease. The identification of 35 was enabled through the optimization and balance of critical properties including potency and pharmacokinetics (PK). This was achieved through modulation of the P2* subsite of the inhibitor which identified the isoquinoline ring system as a key template for improving PK properties with further optimization achieved through functionalization. A methoxy moiety at the C6 position of this isoquinoline ring system proved to be optimal with respect to potency and PK, thus providing the clinical compound 35 which demonstrated antiviral activity in HCV-infected patients.

High resolution crystal structure of human ß-glucuronidase reveals structural basis of lysosome targeting.

Human ß-glucuronidase (GUS) cleaves ß-D-glucuronic acid residues from the non-reducing termini of glycosaminoglycan and its deficiency leads to mucopolysaccharidosis type VII (MPSVII). Here we report a high resolution crystal structure of human GUS at 1.7 Å resolution and present an extensive analysis of the structural features, unifying recent findings in the field of lysosome targeting and glycosyl hydrolases. The structure revealed several new details including a new glycan chain at Asn272, in addition to that previously observed at Asn173, and coordination of the glycan chain at Asn173 with Lys197 of the lysosomal targeting motif which is essential for phosphotransferase recognition. Analysis of the high resolution structure not only provided new insights into the structural basis for lysosomal targeting but showed significant differences between human GUS, which is medically important in its own right, and E. coli GUS, which can be selectively inhibited in the human gut to prevent prodrug activation and is also widely used as a reporter gene by plant biologists. Despite these differences, both human and E. coli GUS share a high structure homology in all three domains with most of the glycosyl hydrolases, suggesting that they all evolved from a common ancestral gene.