Search

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 CHF 2 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 CH 3 and CF 3 analogues, respectively, providing insight into the deployment of this unique amino acid., Competing Interests: The authors declare no competing financial interest.

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., (© 2017 The Protein Society.)

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., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

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

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.

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

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.

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.

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.

Optimization of a 5-oxopyrrolopyridine series based upon structure-activity relationships (SARs) developed from our previous efforts on a number of related bicyclic series yielded compound 2s (BMS-767778) with an overall activity, selectivity, efficacy, PK, and developability profile suitable for progression into the clinic. SAR in the series and characterization of 2s are described.

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., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

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., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

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., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

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.

Continued structure-activity relationship (SAR) exploration within our previously disclosed azolopyrimidine containing dipeptidyl peptidase-4 (DPP4) inhibitors led us to focus on an imidazolopyrimidine series in particular. Further study revealed that by replacing the aryl substitution on the imidazole ring with a more polar carboxylic ester or amide, these compounds displayed not only increased DPP4 binding activity but also significantly reduced human ether-a-go-go related gene (hERG) and sodium channel inhibitory activities. Additional incremental adjustment of polarity led to permeable molecules which exhibited favorable pharmacokinetic (PK) profiles in preclinical animal species. The active site binding mode of these compounds was determined by X-ray crystallography as exemplified by amide 24c. A subsequent lead molecule from this series, (+)-6-(aminomethyl)-5-(2,4-dichlorophenyl)-N-(1-ethyl-1H-pyrazol-5-yl)-7-methylimidazo[1,2-a]pyrimidine-2-carboxamide (24s), emerged as a potent, selective DPP4 inhibitor that displayed excellent PK profiles and in vivo efficacy in ob/ob mice.

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.

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.

An indole-based P1 moiety was incorporated into a previously established factor Xa inhibitor series. The indole group was designed to hydrogen-bond with the carbonyl of Gly218, while its 3-methyl or 3-chloro substituent was intended to interact with Tyr228. These interactions were subsequently observed in the X-ray crystal structure of compound 18. SAR studies led to the identification of compound 20 as the most potent FXa inhibitor in this series (IC(50) = 2.4 nM, EC(2xPT) = 1.2 microM). An in-depth energetic analysis suggests that the increased binding energy of 3-chloroindole-versus 3-methylindole-containing compounds in this series is due primarily to (a) the more hydrophobic nature of chloro- versus methyl-containing compounds and (b) an increased interaction of 3-chloroindole versus 3-methylindole with Gly218 backbone. The stronger hydrophobicity of chloro- versus methyl-substituted aromatics may partly explain the general preference for chloro- versus methyl-substituted P1 groups in FXa, which extends beyond the current series.

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.

Atazanavir, which is marketed as REYATAZ, is the first human immunodeficiency virus type 1 (HIV-1) protease inhibitor approved for once-daily administration. As previously reported, atazanavir offers improved inhibitory profiles against several common variants of HIV-1 protease over those of the other peptidomimetic inhibitors currently on the market. This work describes the X-ray crystal structures of complexes of atazanavir with two HIV-1 protease variants, namely, (i) an enzyme optimized for resistance to autolysis and oxidation, referred to as the cleavage-resistant mutant (CRM); and (ii) the M46I/V82F/I84V/L90M mutant of the CRM enzyme, which is resistant to all approved HIV-1 protease inhibitors, referred to as the inhibitor-resistant mutant. In these two complexes, atazanavir adopts distinct bound conformations in response to the V82F substitution, which may explain why this substitution, at least in isolation, has yet to be selected in vitro or in the clinic. Because of its nearly symmetrical chemical structure, atazanavir is able to make several analogous contacts with each monomer of the biological dimer.

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.

Penicillin G acylase is an important enzyme in the commercial production of semisynthetic penicillins used to combat bacterial infections. Mutant strains of Providencia rettgeri were generated from wild-type cultures subjected to nutritional selective pressure. One such mutant, Bro1, was able to use 6-bromohexanamide as its sole nitrogen source. Penicillin acylase from the Bro1 strain exhibited an altered substrate specificity consistent with the ability of the mutant to process 6-bromohexanamide. The X-ray structure determination of this enzyme was undertaken to understand its altered specificity and to help in the design of site-directed mutants with desired specificities. In this paper, the structure of the Bro1 penicillin G acylase has been solved at 2.5 A resolution by molecular replacement. The R-factor after refinement is 0.154 and R-free is 0.165. Of the 758 residues in the Bro1 penicillin acylase heterodimer (alpha-subunit, 205; beta-subunit, 553), all but the eight C-terminal residues of the alpha-subunit have been modeled based on a partial Bro1 sequence and the complete wild-type P. rettgeri sequence. A tightly bound calcium ion coordinated by one residue from the alpha-subunit and five residues from the beta-subunit has been identified. This enzyme belongs to the superfamily of Ntn hydrolases and uses Ogamma of Ser beta1 as the characteristic N-terminal nucleophile. A mutation of the wild-type Met alpha140 to Leu in the Bro1 acylase hydrophobic specificity pocket is evident from the electron density and is consistent with the observed specificity change for Bro1 acylase. The electron density for the N-terminal Gln of the alpha-subunit is best modeled by the cyclized pyroglutamate form. Examination of aligned penicillin acylase and cephalosporin acylase primary sequences, in conjunction with the P. rettgeri and Escherichia coli penicillin acylase crystal structures, suggests several mutations that could potentially allow penicillin acylase to accept charged beta-lactam R-groups and to function as a cephalosporin acylase and thus be used in the manufacture of semi-synthetic cephalosporins.

Bryonia dioica (Cucurbitaceae family) produces at least two type I ribosome-inactivating proteins, bryodin 1 (BD1) and bryodin 2 (BD2). A cDNA sequence encoding BD1 was isolated from B. dioica leaf mRNA using degenerative oligonucleotides and codes for a 22 amino acid signal peptide followed by a protein of 267 residues. Expression of two recombinant BD1 (rBD1) forms in Escherichia coli yielded proteins of 267 (to the natural stop codon) and 247 amino acids (to the putative cleavage site yielding the mature protein) that had identical protein synthesis inhibition activity as compared to native BD1. The substitution of Lys for Glu at position 189 near the active site reduced the ability of rBD1 to inhibit protein synthesis by 10-fold. Toxicologic analysis showed that rBD1 was well tolerated in rodents with LD50 values of 40 mg/kg in mice and >25 mg/kg in rats. A crystal of mature rBD1 protein was used to collect X-ray diffraction data to 2.1 A resolution. The rBD1 crystal structure was solved and showed extensive homology with other type I RIPs and A chains of type II RIPs. The studies described here demonstrate that rBD1 retains full biologic activity and serve as a guide for using this potent, yet nontoxic, RIP in the construction of single-chain immunotoxin fusion proteins.

Uridine 5'-diphospho-N-acetylenolpyruvylglucosamine reductase (MurB), the second enzyme in the peptidoglycan synthetic pathway of Escherichia coli, has been crystallized in two previously unreported forms, one orthorhombic and the other monoclinic. MurB (molecular mass 38 kDa) crystallizes in a range of conditions that utilize polyethylene glycol fractions as precipitants, and crystals can be grown with or without the enzyme's substrate, uridine 5'-diphospho-N-acetylenolpyruvylglucosamine. X-ray diffraction from crystals of the orthorhombic form extends to 2 A resolution and shows the symmetry and systematic absences of space group P2(1)2(1)2(1). These crystals show significant variations in cell dimensions at room temperature and at 100 K. A crystal used to collect a 2.0 A resolution data set at a synchrotron source showed cell dimensions at ca 100 K of a = 51.0, b = 79.3 and c = 87.1 A, indicating one molecule peroasymmetric unit. The monoclinic crystals scatter X-rays to 3.0 A resolution consistent with space group P2(1), unit-cell dimensions (ca 100 K) a = 50.7, b = 92.4, c = 85.5 A, and beta = 104 degrees, and two molecules per asymmetric unit. Mercury derivatives have been prepared with both orthorhombic and monoclinic forms, and efforts are underway to exploit these derivatives to determine the structure of this protein.

Two isoforms of the heterodimeric enzyme penicillin G acylase (EC 3.5.1.11) from Providencia rettgeri ATCC 31052 (strain Bro1) were purified to near homogeneity. The isoforms exhibited comparable enzymatic activities but differed slightly in the molecular weight and pI of their respective alpha-subunit. The origin of this difference was traced to the partial conversion of the N-terminal Gln of the alpha-subunit to pyrrolidonecarboxylic acid (pyro-Glu). The boundaries of the mature enzyme within the translated DNA sequence of the wild-type propeptide (GenBank M86533) were determined. The results conclusively identified the length of the signal peptide and the position of the spacer cleaved from the propeptide to form the active heterodimer. The molecular weights of the alpha- and beta-subunits, based on these termini, were 23.7 and 62.2 kDa, respectively. Both isoforms were crystallized independently as hexagonal bipyramids up to 0.60 mm in diameter in either space group P6(1)22 or P6(5)22 (a = b = 140.5 A and c = 209.5 A) from ammonium sulfate solutions buffered by 50 mM potassium phosphate at pH 7.5. The presence of glycerol, although not required, facilitated crystal growth. Native and heavy atom derivative data were collected to 3.0 A resolution, and the calculation of isomorphous replacement phases is under way.

A previously developed immobilization technique involving latex coatings on solid particulate supports was investigated further for penicillin G production by Penicillium chrysogenum. Several modifications were found to decrease the germination lag time, including a higher spore concentration, a thinner latex layer, an increased latex porosity, and a decreased drying time. This approach enabled the development of immobilized mycelial pellets within 2-3 days from the onset of biocatalyst preparation and incubation.A continuous immobilized-cell airlift bioreactor produced penicillin G in a series of runs in which the production phase lasted up to 30 days. The productivity of this system was 3-6 times greater than the productivity of the corresponding free-cell shake flask fermentation.