Search

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.

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

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.

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.

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.

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)

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.

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.

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.

We present a method to explore the interaction of flexible ligands with receptors of known geometry on the basis of molecular shape. This method is an extension of that described by Kuntz et al. (J. Mol. Biol. 1982, 161, 269). The shape of a binding site on a macromolecular receptor is represented as a set of overlapping spheres. Each ligand is divided into a small set of large rigid fragments that are docked separately into the binding site and then rejoined later in the calculation. The division of ligands into separate fragments allows a degree of flexibility at the position that joins them. The rejoined fragments are then energy minimized in the receptor site. We illustrate the method with two test cases: dihydrofolate reductase/methotrexate and prealbumin/thyroxine. For each test case, the method finds binding geometries for the ligand near that observed crystallographically as well as others that provide good steric fit with the receptor.

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.