Development and screening of a polyketide virtual library for drug leads against a motilide pharmacophore.

A virtual library of macrocyclic polyketide molecules was generated and screened to identify novel, conformationally constrained potential motilin receptor agonists ("motilides"). A motilide pharmacophore model was generated from the potent 6,9-enol ether erythromycin and known derivatives from the literature. The pharmacophore for each molecular conformation was a point in a distance-volume space based on presentation of the putative binding moieties. Two methods, one fragment based method and the other reaction based, were explored for constructing the polyketide virtual library. First, a virtual library was assembled from monomeric fragments using the CHORTLES language. Second, the virtual library was assembled by the in silico application of all possible polyketide synthase enzyme reactions to generate the product library. Each library was converted to low-energy 3D conformations by distance geometry and standard minimization methods. The distance-volume metric was calculated for low-energy conformations of the members of the virtual polyketide library and screened against the enol ether pharmacophore. The goal was to identify novel macrocycles that satisfy the pharmacophore. We identified three conformationally constrained, novel polyketide series that have low-energy conformations satisfying the distance-volume constraints of the motilide pharmacophore.

Leishmania major: molecular modeling of cysteine proteases and prediction of new nonpeptide inhibitors.

The crystal structures of papain, cruzain, and human liver cathepsin B were used to build homology-based enzyme models of a cathepsin L-like cysteine protease (cpL) and a cathepsin B-like cysteine protease (cpB) from the protozoan parasite Leishmania major. Although structurally a member of the cathepsin B subfamily, the L. major cpB is not able to cleave synthetic substrates having an arginine in position P2. This biochemical property correlates with the prediction of a glycine instead of a glutamic acid at position 205 (papain numbering). The modeled active sites of the L. major cpB and cpL were used to screen the Available Chemicals Directory (a database of about 150,000 commercially available compounds) for potential cysteine protease inhibitors, using DOCK3.5. Based on both steric and force field considerations, 69 compounds were selected. Of these, 18 showed IC50's between 50 and 100 microM and 3 had IC50's below 50 microM. A secondary library of compounds, originally derived from a structural screen against the homologous protease of Plasmodium falciparum (falcipain), and subsequently expanded by combinatorial chemistry, was also screened. Three inhibitors were identified which were not only effective against the L. major protease but also inhibited parasite growth at 5-50 microM.

Structure-based identification of an inducer of the low-pH conformational change in the influenza virus hemagglutinin: irreversible inhibition of infectivity.

Past efforts to employ a structure-based approach to design an inhibitor of the fusion-inducing conformational change in the influenza virus hemagglutinin (HA) yielded a family of small benzoquinones and hydroquinones. The most potent of these, tert-butyl hydroquinone (TBHQ), inhibits both the conformational change in HA from strain X:31 influenza virus and viral infectivity in tissue culture cells with 50% inhibitory concentrations in the micromolar range (D. L. Bodian, R. B. Yamasaki, R. L. Buswell, J. F. Stearns, J. M. White, and I. D. Kuntz, Biochemistry 32:2967-2978, 1993). A new structure-based inhibitor design search was begun which involved (i) the recently refined crystal structure (2.1-A resolution) of the HA ectodomain, (ii) new insights into the conformational change, and (iii) improvements in the molecular docking program, DOCK. As a result, we identified new inhibitors of HA-mediated membrane fusion. Like TBHQ, most of these molecules inhibit the conformational change. One of the new compounds, however, facilitates rather than inhibits the HA conformational change. Nonetheless, the facilitator, diiodofluorescein, inhibits HA-mediated membrane fusion and, irreversibly, infectivity. We further characterized the effects of inhibitors from both searches on the conformational change and membrane fusion activity of HA as well as on viral infectivity. We also isolated and characterized several mutants resistant to each class of inhibitor. The implications of our results for HA-mediated membrane fusion, anti-influenza virus therapy, and structure-based inhibitor design are discussed.

Matching chemistry and shape in molecular docking.

We have added a chemical filter to the ligand placement algorithm of the molecular docking program DOCK. DOCK places ligands in receptors using local shape features. Here we label these shape features by chemical type and insist on complementary matches. We find fewer physically unrealistic complexes without reducing the number of complexes resembling the known ligand-receptor configurations. Approximately 10-fold fewer complexes are calculated and the new algorithm is correspondingly 10-fold faster than the previous shape-only matching. We tested the new algorithm's ability to reproduce three known ligand-receptor complexes: methotrexate in dihydrofolate reductase, deoxyuridine monophosphate in thymidylate synthase and pancreatic trypsin inhibitor in trypsin. The program found configurations within 1 A of the crystallographic mode, with fewer non-native solutions compared with shape-only matching. We also tested the program's ability to retrieve known inhibitors of thymidylate synthase and dihydrofolate reductase by screening molecular databases against the enzyme structures. Both algorithms retrieved many known inhibitors preferentially to other compounds in the database. The chemical matching algorithm generally ranks known inhibitors better than does matching based on shape alone.

Inhibition of the fusion-inducing conformational change of influenza hemagglutinin by benzoquinones and hydroquinones.

Influenza hemagglutinin (HA) undergoes a conformational change that is required for viral entry. The rearrangement includes exposure of the fusion peptide, a hydrophobic segment buried in the trimer interface of the native protein. Since fusion peptide release triggers the membrane fusion event crucial for viral replication, inhibition of fusion peptide exposure should prevent infection. We reasoned that small molecules that bind to HA and stabilize its nonfusogenic conformation would block viral activity. A computer-assisted method was used to select putative HA ligands. One of the selected compounds, 4A,5,8,8A-tetrahydro-5,8-methano-1,4-naphthoquinone, prevented the conversion of X31 HA to a conformation recognized by alpha-fusion peptide antisera. Several derivatives of this compound, including both benzoquinones and hydroquinones, also showed inhibition. The most effective compounds tested have IC50S between 1 and 20 microM. Representative compounds also inhibited virus-induced syncytia formation, HA-mediated hemolysis, and viral infectivity in vitro. The inhibitors are attractive leads for the development of antiviral drugs and can serve as probes of the mechanism of the conformational change of HA.