Synthesis of thiophene-2-carboxamidines containing 2-aminothiazoles and their biological evaluation as urokinase inhibitors.

The serine protease urokinase (uPa) has been implicated in the progression of both breast and prostate cancer. Utilizing structure based design, the synthesis of a series of substituted 4-[2-amino-1,3-thiazolyl]-thiophene-2-carboxamidines is described. Further optimization of this series by substitution of the terminal amine yielded urokinase inhibitors with excellent activities.

Design and synthesis of diaminopyrrolidinone inhibitors of human osteoclast cathepsin K.

The structure-based design and synthesis of lactam-constrained azapeptide inhibitors of human cathepsin K are described. Enhanced stability to proteolytic cleavage over acyclic analogues is discussed.

ASP147 in the third transmembrane helix of the rat mu opioid receptor forms ion-pairing with morphine and naltrexone.

We tested the hypotheses that the carboxylate side chain of Asp147 of the mu opioid receptor interacts with the protonated nitrogen of naltrexone and morphine and that this interaction is important for pharmacological properties of the two compounds. Mutation of Asp147 to Ala or Asn substantially reduced the affinity of naltrexone and the affinity, potency and efficacy of morphine, while the Glu mutant had similar properties as the wildtype, indicating the significant role of the carboxylate group of Asp147 in receptor binding and activation. This role could be due to its direct interaction with ligands or involvement in interhelical interactions. The unprotonated analogs of naltrexone and morphine, cyclopropylcarbonyl noroxymorphone (CPCNOM) and N-formylnormorphine (NFNM), respectively, were used to discriminate between these mechanisms. CPCNOM was much less potent as an antagonist and had substantially lower affinity for the mu receptor than naltrexone. Similarly, NFNM was unable to activate the mu receptor and had much lower affinity than morphine. These results indicate the importance of the protonated nitrogen. Notably, the D147A and D147N mutations did not appreciably affect the binding affinities of CPCNOM and NFNM. In addition, the D147E mutant had similar affinities for CPCNOM and NFNM as the D147A and D147N mu receptors. Thus, the carboxylate group of Asp147 is not important for binding of the two unprotonated compounds. These results indicate that the carboxylate group of Asp147 of the mu receptor interacts directly with the protonated nitrogen of naltrexone and morphine and this interaction is important for binding and receptor activation.

Structure-based design of non-peptide, carbohydrazide-based cathepsin K inhibitors.

Using binding models which were based on the X-ray crystal structure of an amino acid-based active site-spanning inhibitor complexed with cathepsin K, Cbz-leucine mimics have been developed, leading ultimately to the design of a potent cathepsin K inhibitor free of amino acid components. These mimics, which consist of alpha-substituted biphenylacetyl groups in place of Cbz-leucine moieties, effectively mimic all aspects of the Cbz-leucine moieties which are important for inhibitor binding. The predicted directions of binding for the inhibitors were confirmed by mass spectral analysis of their complexes with cathepsin K, which gave results consistent with acylation of the enzyme and loss of the acylhydrazine portion of the inhibitor which binds on the S' side of the active site. The binding models were found to be very predictive of relative inhibitor potency as well as direction of inhibitor binding. These results strengthen the validity of a strategy involving iterative cycles of structure-based design and inhibitor synthesis and evaluation for the discovery of non-peptide inhibitors.

Use of papain as a model for the structure-based design of cathepsin K inhibitors: crystal structures of two papain-inhibitor complexes demonstrate binding to S'-subsites.

Papain has been used as a surrogate enzyme in a drug design effort to obtain potent and selective inhibitors of cathepsin K, a new member of the papain superfamily of cysteine proteases that is selectively and highly expressed in osteoclasts and is implicated in bone resorption. Here we report the crystal structures of two papain-inhibitor complexes and the rational design of novel cathepsin K inhibitors. Unlike previously known crystal structures of papain-inhibitor complexes, our papain structures show ligand binding extending deep within the S'-subsites. The two inhibitor complexes, carbobenzyloxyleucinyl-leucinyl-leucinal and carbobenzyloxy-L-leucinyl-L-leucinyl methoxymethyl ketone, were refined to 2.2- and 2.5-A resolution with R-factors of 0.190 and 0. 217, respectively. The S'-subsite interactions with the inhibitors are dominated by an aromatic-aromatic stacking and an oxygen-aromatic ring edge interaction. The knowledge of S'-subsite interactions led to a design strategy for an inhibitor spanning both subsites and yielded a novel, symmetric inhibitor selective for cathepsin K. Simultaneous exploitation of both S- and S'-sites provides a general strategy for the design of cysteine protease inhibitors having high specificity to their target enzymes.

Structure-based design of cathepsin K inhibitors containing a benzyloxy-substituted benzoyl peptidomimetic.

Peptidomimetic cathepsin K inhibitors have been designed using binding models which were based on the X-ray crystal structure of an amino acid-based, active site-spanning inhibitor complexed with cathepsin K. These inhibitors, which contain a benzyloxybenzoyl group in place of a Cbz-leucine moiety, maintained good inhibitory potency relative to the amino acid-based inhibitor, and the binding models were found to be very predictive of relative inhibitor potency. The binding mode of one of the inhibitors was confirmed by X-ray crystallography, and the crystallographically determined structure is in close qualitative agreement with the initial binding model. These results strengthen the validity of a strategy involving iterative cycles of structure-based design, inhibitor synthesis and evaluation, and crystallographic structure determination for the discovery of peptidomimetic inhibitors.

Design of potent and selective human cathepsin K inhibitors that span the active site.

Potent and selective active-site-spanning inhibitors have been designed for cathepsin K, a cysteine protease unique to osteoclasts. They act by mechanisms that involve tight binding intermediates, potentially on a hydrolytic pathway. X-ray crystallographic, MS, NMR spectroscopic, and kinetic studies of the mechanisms of inhibition indicate that different intermediates or transition states are being represented that are dependent on the conditions of measurement and the specific groups flanking the carbonyl in the inhibitor. The species observed crystallographically are most consistent with tetrahedral intermediates that may be close approximations of those that occur during substrate hydrolysis. Initial kinetic studies suggest the possibility of irreversible and reversible active-site modification. Representative inhibitors have demonstrated antiresorptive activity both in vitro and in vivo and therefore are promising leads for therapeutic agents for the treatment of osteoporosis. Expansion of these inhibitor concepts can be envisioned for the many other cysteine proteases implicated for therapeutic intervention.

Determination of the amino acid residue involved in [3H]beta-funaltrexamine covalent binding in the cloned rat mu-opioid receptor.

We previously demonstrated that [3H]beta-funaltrexamine ([3H]beta-FNA) labeled the rat mu opioid receptor expressed in Chinese hamster ovary cells with high specificity, and [3H]beta-FNA-labeled receptors migrated as one broad band with a mass of 80 kDa. In this study, we determined the region and then the amino acid residue of the mu receptor involved in the covalent binding of [3H]beta-FNA. [3H]beta-FNA-labeled receptors were solubilized and purified to approximately 10% purity by immunoaffinity chromatography with antibodies against a C-terminal domain peptide. The site of covalent bond formation was determined to be within Ala206-Met243 by CNBr cleavage of partially purified labeled mu receptors and determinations of sizes of labeled receptor fragments. The amino acid residue of beta-FNA covalent incorporation was then determined by site-directed mutagenesis studies within this region. Mutation of Lys233 to Ala, Arg, His, and Leu completely eliminated covalent binding of [3H]beta-FNA, although these mutants bound beta-FNA with high affinity. Mutations of other amino acid residues did not affect covalent binding of [3H]beta-FNA. These results indicate that [3H]beta-FNA binds covalently to Lys233. Since [3H]beta-FNA is a rigid molecule, the information will be very useful for molecular modeling of interaction between morphinans and the mu receptor.

A paradigm for drug discovery using a conformation from the crystal structure of a presentation scaffold.

We describe a structural validation of the use of presentation scaffolds for control and elucidation of bioactive conformations of peptides. The protein REI-RGD34--produced by inserting the sequence RIPRGDMP into the CDR1 loop region of the immunoglobulin VL domain REI--strongly inhibits fibrinogen binding to the integrins alpha IIb beta 3 and alpha V beta 3. In the X-ray crystal structure of their protein at 2.4 A resolution, the RGD-containing loop exhibits defined electron density that is consistent with models for the bioactive conformations of ligands of these receptors based on previous small-molecule studies. Furthermore, a search of a small-molecule database with conformational information derived from the structure of REI-RGD34 identified constrained peptides and peptidomimetics known to be antagonists of the platelet receptor alpha IIb beta 3.

A check on rational drug design: crystal structure of a complex of human immunodeficiency virus type 1 protease with a novel gamma-turn mimetic inhibitor.

We have previously reported (Newlander et al., J. Med. Chem. 1993, 36, 2321-2331) the design of human immunodeficiency virus type 1 (HIV-1) protease inhibitors incorporating C7 mimetics that lock three amino acid residues of a peptide sequence into a gamma-turn. The design of one such compound, SB203238, was based on X-ray structures of reduced amide aspartyl protease inhibitors. It incorporates a gamma-turn mimetic in the P2-P1' position, where the carbonyl of the C7 ring is replaced with an sp3 methylene group yielding a constrained reduced amide. It shows competitive inhibition with Ki = 430 nM at pH 6.0. The three-dimensional structure of SB203238 bound to the active site of HIV-1 protease has been determined at 2.3 A resolution by X-ray diffraction and refined to a crystallographic R-factor (R = sigma magnitude of Fo magnitude of - magnitude of Fc magnitude of /sigma magnitude of Fo magnitude of, where Fo and Fc are the observed and calculated structure factor amplitudes, respectively) of 0.177. The inhibitor lies in an extended conformation in the active site; however, because of the constrained geometry of the C7 ring, it maintains fewer hydrogen bonds with the protein than in most other HIV-1 protease-inhibitor complexes. More importantly, the inhibitor binds to the enzyme differently than predicted in its design, by binding with the P2-P1' alpha-carbon atoms shifted by approximately one-half a residue toward the N-terminus from their presumed positions. This study illustrates the importance of structural information in an approach to rational drug design.

An orally bioavailable HIV-1 protease inhibitor containing an imidazole-derived peptide bond replacement: crystallographic and pharmacokinetic analysis.

(2R,4S,5S,1'S)-2-Phenylmethyl-4-hydroxy-5-(tert-butoxycarbonyl) amino-6-phenylhexanoyl-N-(1'-imidazo-2-yl)-2'-methylpropanamide (compound 2) is a tripeptide analogue inhibitor of HIV-1 protease in which a C-terminal imidazole substituent constitutes an isoelectronic, structural mimic of a carboxamide group. Compound 2 is a potent inhibitor of the protease (K(i) = 18 nM) and inhibits HIV-1 acute infectivity of CD4+ T-lymphocytes (IC50 = 570 nM). Crystallographic analysis of an HIV-1 protease-compound 2 complex demonstrates that the nitrogen atoms of the imidazole ring assume the same hydrogen-bonding interactions with the protease as amide linkages in other peptide analogue inhibitors. The sole substitution of the C-terminal carboxamide of a hydroxyethylene-containing tripeptide analogue with an imidazole group imparts greatly improved pharmacokinetic and oral bioavailability properties on the compound compared to its carboxamide-containing homologue (compound 1). While the oral bioavailability of compound 1 in rats was negligible, compound 2 displayed oral bioavailabilities of 30% and 14%, respectively, in rats and monkeys.

Rational design, synthesis, and crystallographic analysis of a hydroxyethylene-based HIV-1 protease inhibitor containing a heterocyclic P1'--P2' amide bond isostere.

The rational design and synthesis of a highly potent inhibitor of HIV-1 protease have been accomplished. The inhibitor, SB 206343, is based on a model derived from the structure of the MVT-101/HIV-1 protease complex and contains a 4(5)-acylimidazole ring as an isosteric replacement for the P1'--P2' amide bond. It is a competitive inhibitor with an apparent inhibition constant of 0.6 nM at pH 6.0. The three-dimensional structure of SB 206343 bound in the active site of HIV-1 protease has been determined at 2.3 A resolution by X-ray diffraction techniques and refined to a crystallographic discrepancy factor, R (= sigma parallel Fo magnitude of/Fc parallel/sigma magnitude of), of 0.194. The inhibitor is held in the enzyme by a set of hydrophobic and polar interactions. N-3 of the imidazole ring participates in a novel hydrogen-bonding interaction with the bound water molecule, demonstrating the effectiveness of the imidazole ring as an isosteric replacement for the P1'--P2' amide bond in hydroxyethylene-based HIV-1 protease inhibitors. Also present are hydrogen-bonding interactions between N-1 of the imidazole ring and the carbonyl of Gly-127 as well as between the imidazole acyl carbonyl oxygen and the amide nitrogen of Asp-129, exemplifying the peptidomimetic nature of the 4(5)-acylimidazole isostere. All of these interactions are in qualitative agreement with those predicted by the model.

A shape- and chemistry-based docking method and its use in the design of HIV-1 protease inhibitors.

The program DOCK [1,2] has been used successfully to identify molecules which will bind to a specified receptor [3]. The original method ranks molecules based on their shape complementarity to the receptor site and relies on the chemist to bring the appropriate electrostatic or hydrogen bond properties into the molecular skeletons obtained in the search. This is useful when screening a small database of compounds, where it is not likely that molecules with both the correct shape and electrostatic properties will be found. As large databases are more likely to have redundant molecular shapes with a variety of functionality (e.g., members of a congeneric series), it would be useful to have a method which identifies molecules with both the correct shape and functionality. To this end we have modified the DOCK 1.0 method to target user-specified atom types to selected positions in the receptor site. The target sites can be chosen based on structural evidence, calculation or inspection. Targeted-DOCK improves the ability of the DOCK method to find the crystallographically determined binding mode of a ligand. Additionally, targeted-DOCK searches a database of small molecules at 100-1000 times the rate of DOCK 1.0, allowing more molecules to be screened and more sophisticated scoring schemes to be employed. Targeted-DOCK has been used successfully in the design of a novel non-peptide inhibitor of HIV-1 protease.

Strong inhibition of fibrinogen binding to platelet receptor alpha IIb beta 3 by RGD sequences installed into a presentation scaffold.

In order to probe the structural constraints on binding of RGD sequences to the platelet receptor alpha IIb beta 3 we have used recombinant DNA techniques to install the RGD sequence into 'presentation scaffolds', small proteins of known 3-D structure chosen to present guest sequences in constrained orientations. Using Escherichia coli expression systems we made sequence variants in which loop residues of the immunoglobulin VL domain REI and of human interleukin-1 beta were replaced (without changing polypeptide length) by the RGD sequence at positions predicted, based on small molecule studies, to orient the RGD moiety into an active conformation. These variants do not compete for fibrinogen binding to alpha IIb beta 3 up to almost 1 mM concentration. Unfolded or proteolytically fragmented forms of these same proteins do compete, however, showing that the RGD sequences in the mutants must be prohibited from binding by constraints imposed by scaffold structure. To suppress the effects of such structural constraints we constructed two sequence variants in which RGD-containing sequences 42-57 or 44-55 from the snake venom platelet antagonist kistrin were inserted (this increasing the length of the loop) into the third complementarity determining loop of REI. Both of these variants compete strongly for fibrinogen binding with IC50s in the nM range. These results, plus data on kistrin-related peptides also presented here, suggest that the molecular scaffold REI is capable of providing to an installed sequence a structural context and conformation beneficial to binding. The results also suggest that in order to bind well to alpha IIb beta 3, RGD sequences in protein ligands must either project significantly from the surface of the scaffold and/or retain a degree of conformational flexibility within the scaffold. Molecular scaffolds like REI should prove useful in the elucidation of structure-function relationships and the discovery of new active sequences, and may also serve as the basis for novel therapeutic agents.

Inhibition of human immunodeficiency virus-1 protease by a C2-symmetric phosphinate. Synthesis and crystallographic analysis.

The human immunodeficiency virus type 1 (HIV-1) protease is a potential target of acquired immune deficiency syndrome (AIDS) therapy. A highly potent, perfectly symmetrical phosphinate inhibitor of this enzyme, SB204144, has been synthesized. It is a competitive inhibitor of HIV-1 protease, with an apparent inhibition constant of 2.8 nM at pH 6.0. The three-dimensional structure of SB204144 bound to the enzyme has been determined at 2.3-A resolution by X-ray diffraction techniques and refined to a crystallographic discrepancy factor, R (= sigma parallel F(o) magnitude to - Fc parallel/sigma magnitude of F(o)), of 0.178. The inhibitor is held in the enzyme active site by a set of hydrophobic and hydrophilic interactions, including an interaction between Arg8 and the center of the terminal benzene rings of the inhibitor. The phosphinate establishes a novel interaction with the two catalytic aspartates; each oxygen of the central phosphinic acid moiety interacts with a single oxygen of one aspartic acid, establishing a very short (2.2-2.4 A) oxygen-oxygen contact. As with the structures of penicillopepsin bound to phosphinate and phosphonate inhibitors [Fraser, M. E., Strynadka, N. C., Bartlett, P. A., Hanson, J. E., & James, M. N. (1992) Biochemistry 31, 5201-14], we interpret this short distance and the stereochemical environment of each pair of oxygens in terms of a hydrogen bond that has a symmetric single-well potential energy curve with the proton located midway between the two atoms.(ABSTRACT TRUNCATED AT 250 WORDS)

A symmetric inhibitor binds HIV-1 protease asymmetrically.

Potential advantages of C2-symmetric inhibitors designed for the symmetric HIV-1 protease include high selectivity, potency, stability, and bioavailability. Pseudo-C2-symmetric monools and C2-symmetric diols, containing central hydroxymethylene and (R,R)-dihydroxyethylene moieties flanked by a variety of hydrophobic P1/P1' side chains, were studied as HIV-1 protease inhibitors. The monools and diols were synthesized in 8-10 steps from D-(+)-arabitol and D-(+)-mannitol, respectively. Monools with ethyl or isobutyl P1/P1' side chains were weak inhibitors of recombinant HIV-1 protease (Ki > 10 microM), while benzyl P1/P1' side chains afforded a moderately potent inhibitor (apparent Ki = 230 nM). Diols were 100-10,000x more potent than analogous monools, and a wider range of P1/P1' side chains led to potent inhibition. Both classes of compounds exhibited lower apparent Ki values under high-salt conditions. Surprisingly, monool and diol HIV-1 protease inhibitors were potent inhibitors of porcine pepsin, a prototypical asymmetric monomeric aspartic protease. These results were evaluated in the context of the pseudosymmetric structure of monomeric aspartic proteases and their evolutionary kinship with the retroviral proteases. The X-ray crystal structure of HIV-1 protease complexed with a symmetric diol was determined at 2.6 A. Contrary to expectations, the diol binds the protease asymmetrically and exhibits 2-fold disorder in the electron density map. Molecular dynamics simulations were conducted beginning with asymmetric and symmetric HIV-1 protease/inhibitor model complexes. A more stable trajectory resulted from the asymmetric complex, in agreement with the observed asymmetric binding mode. A simple four-point model was used to argue more generally that van der Waals and electrostatic force fields can commonly lead to an asymmetric association between symmetric molecules.

Structure-based design of nonpeptide inhibitors specific for the human immunodeficiency virus 1 protease.

By using a structure-based computer-assisted search, we have found a butyrophenone derivative that is a selective inhibitor of the human immunodeficiency virus 1 (HIV-1) protease. The computer program creates a negative image of the active site cavity using the crystal structure of the HIV-1 protease. This image was compared for steric complementarity with 10,000 molecules of the Cambridge Crystallographic Database. One of the most interesting candidates identified was bromperidol. Haloperidol, a closely related compound and known antipsychotic agent, was chosen for testing. Haloperidol inhibits the HIV-1 and HIV-2 proteases in a concentration-dependent fashion with a Ki of approximately 100 microM. It is highly selective, having little inhibitory effect on pepsin activity and no effect on renin at concentrations as high as 5 mM. The hydroxy derivative of haloperidol has a similar effect on HIV-1 protease but a lower potency against the HIV-2 enzyme. Both haloperidol and its hydroxy derivative showed activity against maturation of viral polypeptides in a cell assay system. Although this discovery holds promise for the generation of nonpeptide protease inhibitors, we caution that the serum concentrations of haloperidol in normal use as an antipsychotic agent are less than 10 ng/ml (0.03 microM). Thus, concentrations required to inhibit the HIV-1 protease are greater than 1000 times higher than the concentrations normally used. Haloperidol is highly toxic at elevated doses and can be life-threatening. Haloperidol is not useful as a treatment for AIDS but may be a useful lead compound for the development of an antiviral pharmaceutical.

Computerized selection of potential DNA binding compounds.

Using a general shape-search docking algorithm, potential DNA minor groove binding compounds were selected from a subset from the Cambridge Crystallographic Database consisting of almost 10,000 molecules. The crystal structure of the DNA dodecamer as observed in the d(CGCGAATTCGCG)2.netropsin complex served as the target receptor. Surprisingly, the highest scoring compound turned out to be CC-1065, a potent anti-tumour agent. Netropsin itself was number 6 on the list of highest scoring compounds. A number of the top-10 scoring compounds may serve as a source of inspiration for further drug design.

Using shape complementarity as an initial screen in designing ligands for a receptor binding site of known three-dimensional structure.

Finding novel leads from which to design drug molecules has traditionally been a matter of screening and serendipity. We present a method for finding a wide assortment of chemical structures that are complementary to the shape of a macromoleculer receptor site whose X-ray crystallographic structure is known. Each of a set of small molecules from the Cambridge Crystallographic Database (Allen; et al. J. Chem. Doc. 1973, 13, 119) is individually docked to the receptor in a number of geometrically permissible orientations with use of the docking algorithm developed by Kuntz et al. (J. Mol. Biol. 1982, 161, 269). The orientations are evaluated for goodness-of-fit, and the best are kept for further examination using the molecular mechanics program AMBER (Weiner; Kollman J. Comput. Chem. 1981, 106, 765). The shape-search algorithm finds known ligands as well as novel molecules that fit the binding site being studied. The highest scoring orientations of known ligands resemble binding modes generated by interactive modeling or determined crystallographically. We describe the application of this procedure to the binding sites of papain and carbonic anhydrase. While the compounds recovered from the Cambridge Crystallographic Database are not, themselves, likely to be inhibitors or substrates of these enzymes, we expect that the structures from such searches will be useful in the design of active compounds.