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

Optimizing the use of open-source software applications in drug discovery

Drug discovery is a time consuming and costly process. Recently, a trend towards the use of in silico computational chemistry and molecular modeling for computer-aided drug design has gained significant momentum. This review investigates the application of free and/or open-source software in the drug discovery process. Among the reviewed software programs are applications programmed in JAVA, Perl and Python, as well as resources including software libraries. These programs might be useful for cheminformatics approaches to drug discovery, including QSAR studies, energy minimization and docking studies in drug design endeavors. Furthermore, this review explores options for integrating available computer modeling open-source software applications in drug discovery programs.

Decoys for docking

Molecular docking is widely used to predict novel lead compounds for drug discovery. Success depends on the quality of the docking scoring function, among other factors. An imperfect scoring function can mislead by predicting incorrect ligand geometries or by selecting nonbinding molecules over true ligands. These false-positive hits may be considered "decoys". Although these decoys are frustrating, they potentially provide important tests for a docking algorithm; the more subtle the decoy, the more rigorous the test. Indeed, decoy databases have been used to improve protein structure prediction algorithms and protein-protein docking algorithms. Here, we describe 20 geometric decoys in five enzymes and 166 "hit list" decoys-i.e., molecules predicted to bind by our docking program that were tested and found not to do so-for beta-lactamase and two cavity sites in lysozyme. Especially in the cavity sites, which are very simple, these decoys highlight particular weaknesses in our scoring function. We also consider the performance of five other widely used docking scoring functions against our geometric and hit list decoys. Intriguingly, whereas many of these other scoring functions performed better on the geometric decoys, they typically performed worse on the hit list decoys, often highly ranking molecules that seemed to poorly complement the model sites. Several of these "hits"from the other scoring functions were tested experimentally and found, in fact, to be decoys. Collectively, these decoys provide a tool for the development and improvement of molecular docking scoring functions. Such improvements may, in turn, be rapidly tested experimentally against these and related experimental systems, which are well-behaved in assays and for structure determination.

Virtual screening of human 5-aminoimidazole-4-carboxamide ribonucleotide transformylase against the NCI diversity set by use of AutoDock to identify novel nonfolate inhibitors

AICAR transformylase (5-aminoimidazole-4-carboxamide ribonucleotide transformylase) is a folate-dependent activity of the bifunctional protein ATIC (AICAR transformylase and IMP cyclohydrolase) and is responsible for catalyzing the penultimate step of the de novo purine biosynthetic pathway. As such, AICAR transformylase has been proposed as a potential target for antineoplastic drug design. Virtual screening of the human AICAR transformylase active site by use of AutoDock against the NCI diversity set, a library of compounds with nonredundant pharmacophore profiles, has revealed 44 potential inhibitor candidates. In vitro inhibition assay of 16 soluble compounds from this list revealed that eight compounds with novel scaffolds, relative to the general folate template, had micromolar inhibition. Subsequent extension of docking trials on compounds with similar scaffolds from the entire NCI-3D database has unveiled 11 additional inhibitors that were confirmed by the in vitro inhibition assay. In particular, one compound, NSC30171, had nanomolar inhibition (K(i) = 154 nM, IC(50) = 600 nM) against AICAR transformylase. These 19 inhibitors serve as novel templates/scaffolds for development of more potent and specific non-folate-based AICAR transformylase inhibitors.

Structure-based discovery of human l-xylulose reductase inhibitors from database screening and molecular docking

Human L-xylulose reductase (XR) is an enzyme of the glucuronic acid/uronate cycle of glucose metabolism and is a possible target for treatment of the long-term complications of diabetes. In this study we utilised the molecular modelling program DOCK to analyse the 249,071 compounds of the National Cancer Institute Database and retrieved those compounds with high predicted affinity for XR. Several carboxylic acid-based compounds were tested and shown to inhibit XR. These included nicotinic acid (IC50=100 microM), benzoic acid (IC50=29 microM) and their derivatives. These results extend and improve upon the activities of known, commercially available inhibitors of XR such as the aliphatic fatty acid n-butyric acid (IC50=64 microM). To optimise the interaction between the inhibitor and the holoenzyme, the program GRID was used to design de novo compounds based on the inhibitor benzoic acid. The inclusion of a hydroxy-phenyl group and a phosphate to the benzoic acid molecule increased the net binding energy by 1.3- and 2.4-fold, respectively. The resultant compounds may produce inhibitors with improved specificity for XR.

Efficient method for high-throughput virtual screening based on flexible docking: discovery of novel acetylcholinesterase inhibitors

A method of easily finding ligands, with a variety of core structures, for a given target macromolecule would greatly contribute to the rapid identification of novel lead compounds for drug development. We have developed an efficient method for discovering ligand candidates from a number of flexible compounds included in databases, when the three-dimensional (3D) structure of the drug target is available. The method, named ADAM&EVE, makes use of our automated docking method ADAM, which has already been reported. Like ADAM, ADAM&EVE takes account of the flexibility of each molecule in databases, by exploring the conformational space fully and continuously. Database screening has been made much faster than with ADAM through the tuning of parameters, so that computational screening of several hundred thousand compounds is possible in a practical time. Promising ligand candidates can be selected according to various criteria based on the docking results and characteristics of compounds. Furthermore, we have developed a new tool, EVE-MAKE, for automatically preparing the additional compound data necessary for flexible docking calculation, prior to 3D database screening. Among several successful cases of lead discovery by ADAM&EVE, the finding of novel acetylcholinesterase (AChE) inhibitors is presented here. We performed a virtual screening of about 160 000 commercially available compounds against the X-ray crystallographic structure of AChE. Among 114 compounds that could be purchased and assayed, 35 molecules with various core structures showed inhibitory activities with IC(50) values less than 100 microM. Thirteen compounds had IC(50) values between 0.5 and 10 microM, and almost all their core structures are very different from those of known inhibitors. The results demonstrate the effectiveness and validity of the ADAM&EVE approach and provide a starting point for development of novel drugs to treat Alzheimer's disease.

Crystal structures of a multidrug-resistant human immunodeficiency virus type 1 protease reveal an expanded active-site cavity

The goal of this study was to use X-ray crystallography to investigate the structural basis of resistance to human immunodeficiency virus type 1 (HIV-1) protease inhibitors. We overexpressed, purified, and crystallized a multidrug-resistant (MDR) HIV-1 protease enzyme derived from a patient failing on several protease inhibitor-containing regimens. This HIV-1 variant contained codon mutations at positions 10, 36, 46, 54, 63, 71, 82, 84, and 90 that confer drug resistance to protease inhibitors. The 1.8-angstrom (A) crystal structure of this MDR patient isolate reveals an expanded active-site cavity. The active-site expansion includes position 82 and 84 mutations due to the alterations in the amino acid side chains from longer to shorter (e.g., V82A and I84V). The MDR isolate 769 protease "flaps" stay open wider, and the difference in the flap tip distances in the MDR 769 variant is 12 A. The MDR 769 protease crystal complexes with lopinavir and DMP450 reveal completely different binding modes. The network of interactions between the ligands and the MDR 769 protease is completely different from that seen with the wild-type protease-ligand complexes. The water molecule-forming hydrogen bonds bridging between the two flaps and either the substrate or the peptide-based inhibitor are lacking in the MDR 769 clinical isolate. The S1, S1', S3, and S3' pockets show expansion and conformational change. Surface plasmon resonance measurements with the MDR 769 protease indicate higher k(off) rates, resulting in a change of binding affinity. Surface plasmon resonance measurements provide k(on) and k(off) data (K(d) = k(off)/k(on)) to measure binding of the multidrug-resistant protease to various ligands. This MDR 769 protease represents a new antiviral target, presenting the possibility of designing novel inhibitors with activity against the open and expanded protease forms.

A structure-based strategy to identify new molecular scaffolds targeting the bacterial ribosomal A-site

The need for novel antibiotics is widely recognized. A well validated target of antibiotics is the bacterial ribosome. Recent X-ray structures of the ribosome bound to antibiotics have shed new light on the binding sites of these antibiotics, providing fresh impetus for structure-based strategies aiming at identifying new ribosomal ligands. In that respect, the ribosomal decoding region of the aminoacyl-tRNA acceptor site (A-site) is of particular interest because oligonucleotide model systems of this site are available for crystallography, NMR and compound binding assays. This work presents how these different resources can be combined in a hierarchical screening strategy which has led to the identification of new A-site ligands. The approach exploits an X-ray structure of the A-site against which large and diverse libraries of compounds were computationally docked. The complementarity of the compounds to the A-site was assessed using a scoring function specifically calibrated for RNA targets. Starting from approximately 1 million compounds, the computational selection of candidate ligands allowed us to focus the experimental work on 129 compounds, 34 of which showed affinity for the A-site in a FRET-based binding assay. NMR experiments confirmed binding to the A-site for some compounds. For the most potent compound in the FRET assay, a tentative binding mode is suggested, which is compatible with the NMR data and the limited SAR in this series. Overall, the results validate the screening strategy.

Predicting peptide binding to MHC pockets via molecular modeling, implicit solvation, and global optimization

Development of a computational prediction method based on molecular modeling, global optimization, and implicit solvation has produced accurate structure and relative binding affinity predictions for peptide amino acids binding to five pockets of the MHC molecule HLA-DRB1*0101. Because peptide binding to MHC molecules is essential to many immune responses, development of such a method for understanding and predicting the forces that drive binding is crucial for pharmaceutical design and disease treatment. Underlying the development of this prediction method are two hypotheses. The first is that pockets formed by the peptide binding groove of MHC molecules are independent, separating the prediction of peptide amino acids that bind within individual pockets from those that bind between pockets. The second hypothesis is that the native state of a system composed of an amino acid bound to a protein pocket corresponds to the system's lowest free energy. The prediction method developed from these hypotheses uses atomistic-level modeling, deterministic global optimization, and three methods of implicit solvation: solvent-accessible area, solvent-accessible volume, and Poisson-Boltzmann electrostatics. The method predicts relative binding affinities of peptide amino acids for pockets of HLA-DRB1*0101 by determining computationally an amino acid's global minimum energy conformation. Prediction results from the method are in agreement with X-ray crystallography data and experimental binding assays.

Effect of haloperidol and risperidone on amyloid precursor protein levels in vivo

The neurotoxic beta-amyloid peptide of Alzheimer's disease is formed from the amyloid precursor protein (APP), which is a member of an evolutionarily highly conserved gene family with significant functional importance. Because behavioral and psychiatric symptoms treated with antipsychotics may influence the course of the disease, we have investigated traditional and atypical antipsychotic drugs, administered through the intraperitoneal route, for their effects on rat cortical APP. Western-immunoblotting was utilized for semi-quantitative evaluation of APP levels. Treatment with haloperidol resulted in an acute elevation of cortical APP both in therapeutic and toxic doses, however, it had no significant chronic impact on APP. Atypical antipsychotic risperidone did not change cortical APP concentration. These results indicate that both haloperidol and risperidone are considered to be relatively safe with respect to APP metabolism. Possible mechanisms, including involvement of calcium and APP itself as a receptor, are discussed.

Molecular recognition and docking algorithms

Molecular docking is an invaluable tool in modern drug discovery. This review focuses on methodological developments relevant to the field of molecular docking. The forces important in molecular recognition are reviewed and followed by a discussion of how different scoring functions account for these forces. More recent applications of computational chemistry tools involve library design and database screening. Last, we summarize several critical methodological issues that must be addressed in future developments.

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.

Pharmacophoric features of nucleosidic HIV-1RT inhibitors

A quantum pharmacological study on nucleosidic inhibitors of HIV-1RT has been performed. The main aim of this study is to discuss the pharmacophoric features (conformational and electrostatic) of nucleosidic inhibitors and compare them with normal substrate dNTP. Present study stresses on the need to refine nucleosidic drugs, as combination therapy to date is still one of the best remedies for AIDS. The results of ab initio HF calculations indicate very little effect of 3' substituent on ring puckering and suggest that potency regulation may be via very intricate phosphate-catalytic triad interactions. Our MESP maps also show charge complementarity between the drug and receptor.

Virtual screening to enrich hit lists from high-throughput screening: a case study on small-molecule inhibitors of angiogenin

"Hit lists" generated by high-throughput screening (HTS) typically contain a large percentage of false positives, making follow-up assays necessary to distinguish active from inactive substances. Here we present a method for improving the accuracy of HTS hit lists by computationally based virtual screening (VS) of the corresponding chemical libraries and selecting hits by HTS/VS consensus. This approach was applied in a case study on the target-enzyme angiogenin, a potent inducer of angiogenesis. In conjunction with HTS of the National Cancer Institute Diversity Set and ChemBridge DIVERSet E (approximately 18,000 compounds total), VS was performed with two flexible library docking/scoring methods, DockVision/Ludi and GOLD. Analysis of the results reveals that dramatic enrichment of the HTS hit rate can be achieved by selecting compounds in consensus with one or both of the VS functions. For example, HTS hits ranked in the top 2% by GOLD included 42% of the true hits, but only 8% of the false positives; this represents a sixfold enrichment over the HTS hit rate. Notably, the HTS/VS method was effective in selecting out inhibitors with midmicromolar dissociation constants typical of leads commonly obtained in primary screens.

A model binding site for testing scoring functions in molecular docking

Prediction of interaction energies between ligands and their receptors remains a major challenge for structure-based inhibitor discovery. Much effort has been devoted to developing scoring schemes that can successfully rank the affinities of a diverse set of possible ligands to a binding site for which the structure is known. To test these scoring functions, well-characterized experimental systems can be very useful. Here, mutation-created binding sites in T4 lysozyme were used to investigate how the quality of atomic charges and solvation energies affects molecular docking. Atomic charges and solvation energies were calculated for 172,118 molecules in the Available Chemicals Directory using a semi-empirical quantum mechanical approach by the program AMSOL. The database was first screened against the apolar cavity site created by the mutation Leu99Ala (L99A). Compared to the electronegativity-based charges that are widely used, the new charges and desolvation energies improved ranking of known apolar ligands, and better distinguished them from more polar isosteres that are not observed to bind. To investigate whether the new charges had predictive value, the non-polar residue Met102, which forms part of the binding site, was changed to the polar residue glutamine. The structure of the resulting Leu99Ala and Met102Gln double mutant of T4 lysozyme (L99A/M102Q) was determined and the docking calculation was repeated for the new site. Seven representative polar molecules that preferentially docked to the polar versus the apolar binding site were tested experimentally. All seven bind to the polar cavity (L99A/M102Q) but do not detectably bind to the apolar cavity (L99A). Five ligand-bound structures of L99A/M102Q were determined by X-ray crystallography. Docking predictions corresponded to the crystallographic results to within 0.4A RMSD. Improved treatment of partial atomic charges and desolvation energies in database docking appears feasible and leads to better distinction of true ligands. Simple model binding sites, such as L99A and its more polar variants, may find broad use in the development and testing of docking algorithms.

A common mechanism underlying promiscuous inhibitors from virtual and high-throughput screening

High-throughput and virtual screening are widely used to discover novel leads for drug design. On examination, many screening hits appear non-drug-like: they act noncompetitively, show little relationship between structure and activity, and have poor selectivity. Attempts to develop these peculiar molecules into viable leads are often futile, and much time can be wasted on the characterization of these "phony" hits. Despite their common occurrence, the mechanism of action of these promiscuous molecules remains unknown. To investigate this problem, 45 diverse screening hits were studied. Fifteen of these were previously reported as inhibitors of various receptors, including beta-lactamase, malarial protease, dihydrofolate reductase, HIV Tar RNA, thymidylate synthase, kinesin, insulin receptor, tyrosine kinases, farnesyltransferase, gyrase, prions, triosephosphate isomerase, nitric oxide synthase, phosphoinositide 3-kinase, and integrase; 30 were from an in-house screening library of a major pharmaceutical company. In addition to their original targets, 35 of these 45 compounds were shown to inhibit several unrelated model enzymes. These 35 screening hits included compounds, such as fullerenes, dyes, and quercetin, that have repeatedly shown activity against diverse targets. When tested against the model enzymes, the compounds showed time-dependent but reversible inhibition that was dramatically attenuated by albumin, guanidinium, or urea. Surprisingly, increasing the concentration of the model enzymes 10-fold largely eliminated inhibition, despite a 1000-fold excess of inhibitor; a well-behaved competitive inhibitor did not show this behavior. One model to explain these observations was that the active form of the promiscuous inhibitors was an aggregate of many individual molecules. To test this hypothesis, light scattering and electron microscopy experiments were performed. The nonspecific inhibitors were observed to form particles of 30-400 nm diameter by both techniques. In control experiments, a well-behaved competitive inhibitor and an inactive dye-like molecule were not observed to form aggregates. Consistent with the hypothesis that the aggregates are the inhibitory species, the particle size and IC(50) values of the promiscuous inhibitors varied monotonically with ionic strength; a competitive inhibitor was unaffected by changes in ionic strength. Unexpectedly, aggregate formation appears to explain the activity of many nonspecific inhibitors and may account for the activity of many promiscuous screening hits. Molecules acting via this mechanism may be widespread in drug discovery screening databases. Recognition of these compounds may improve screening results in many areas of pharmaceutical interest.

Peptide interactions with G-protein coupled receptors

Peptide recognition by G-protein coupled receptors (GPCRs) is reviewed with an emphasis on the indirect approach used to determine the receptor-bound conformation of peptide ligands. This approach was developed in response to the lack of detailed structural information available for these receptors. Recent advances in the structural determination of rhodopsin (the GPCR of the visual system) by crystallography have provided a scaffold for homology modeling of the inactive state of a wide variety of GPCRs that interact with peptide messages. Additionally, the ability to mutate GPCRs and assay compounds of similar chemical structure to test a common binding site on the receptor provides a firm experimental basis for structure-activity studies. Recognition motifs, common in other well-studied systems such as proteolytic enzymes and major histocompatibility class receptors (MHC) are reviewed briefly to provide a basis of comparison. Finally, the development of true peptidomimetics is contrasted with nonpeptide ligands, discovered through combinatorial chemistry. In many systems, the evidence suggests that the peptide ligands bind at the interface between the transmembrane segments and the extracellular loops, while nonpeptide antagonists bind within the transmembrane segments. Plausible models of GPCRs and the mechanism by which they activate G-proteins on binding peptides are beginning to emerge.

Anti-HIV activity of medicinal herbs: usage and potential development

The acquired immunodeficiency syndrome (AIDS) is a result of human immunodeficiency virus (HIV) infection which subsequently leads to significant suppression of immune functions. AIDS is a significant threat to the health of mankind, and the search for effective therapies to treat AIDS is of paramount importance. Several chemical anti-HIV agents have been developed. However, besides the high cost, there are adverse effects and limitations associated with using chemotherapy for the treatment of HIV infection. Thus, herbal medicines have frequently been used as an alternative medical therapy by HIV positive individuals and AIDS patients. The aim of this review is to summarize research findings for herbal medicines, which are endowed with the ability to inhibit HIV. In this article, we will emphasize a Chinese herbal medicine, Scutellaria baicalensis Georgi and its identified components (i.e., baicalein and baicalin), which have been shown to inhibit infectivity and replication of HIV. Potential development of anti-AIDS compounds using molecular modeling methods will also be discussed.

Novel inhibitors of Trypanosoma cruzi dihydrofolate reductase

There is an urgent need for the development of new drugs to treat Chagas' disease, which is caused by the protozoan parasite Trypanosoma cruzi. The enzyme dihydrofolate reductase (DHFR) has been a very successful drug target in a number of diseases and we decided to investigate it as a potential drug target for Chagas' disease. A homology model of the enzyme was used to search the Cambridge Structural Database using the program DOCK 3.5. Compounds were then tested against the enzyme and the whole parasite. Compounds were also screened against the related parasite, Trypanosoma brucei.

Polyoxometalate HIV-1 protease inhibitors. A new mode of protease inhibition

Nb-containing polyoxometalates (POMs) of the Wells-Dawson class inhibit HIV-1 protease (HIV-1P) by a new mode based on kinetics, binding, and molecular modeling studies. Reaction of alpha(1)-K(9)Li[P(2)W(17)O(61)] or alpha(2)-K(10)[P(2)W(17)O(61)] with aqueous H(2)O(2) solutions of K(7)H[Nb(6)O(19)] followed by treatment with HCl and KCl and then crystallization affords the complexes alpha(1)-K(7)[P(2)W(17)(NbO(2))O(61)] (alpha(1)()1) and alpha(2)-K(7)[P(2)W(17)(NbO(2))O(61)] (alpha(2)()1) in 63 and 86% isolated yields, respectively. Thermolysis of the crude peroxoniobium compounds (72-96 h in refluxing H(2)O) prior to treatment with KCl converts the peroxoniobium compounds to the corresponding polyoxometalates (POMs), alpha(1)-K(7)[P(2)W(17)NbO(62)] (alpha(1)()2) and alpha(2)-K(7)[P(2)W(17)NbO(62)] (alpha(2)()2), in moderate yields (66 and 52%, respectively). The identity and high purity of all four compounds were confirmed by (31)P NMR and (183)W NMR. The acid-induced dimerization of the oxo complexes differentiates sterically between the cap (alpha(2)) site and the belt (alpha(1)) site in the Wells-Dawson structure (alpha(2)()2 dimerizes in high yield; alpha(1)()2 does not). All four POMs exhibit high activity in cell culture against HIV-1 (EC(50) values of 0.17-0.83 microM), are minimally toxic (IC(50) values of 50 to >100 microM), and selectively inhibit purified HIV-1 protease (HIV-1P) (IC(50) values for alpha(1)()1, alpha(2)()1, alpha(1)()2, and alpha(2)()2 of 2.0, 1.2, 1.5, and 1.8 microM, respectively). Thus, theoretical, binding, and kinetics studies of the POM/HIV-1P interaction(s) were conducted. Parameters for [P(2)W(17)NbO(62)](7)(-) were determined for the Kollman all-atom (KAA) force field in Sybyl 6.2. Charges for the POM were obtained from natural population analysis (NPA) at the HF/LANL2DZ level of theory. AutoDock 2.2 was used to explore possible binding locations for the POM with HIV-1P. These computational studies strongly suggest that the POMs function not by binding to the active site of HIV-1P, the mode of inhibition of all other HIV-1P protease inhibitors, but by binding to a cationic pocket on the "hinge" region of the flaps covering the active site (2 POMs and cationic pockets per active homodimer of HIV-1P). The kinetics and binding studies, conducted after the molecular modeling, are both in remarkable agreement with the modeling results: 2 POMs bind per HIV-1P homodimer with high affinities (K(i) = 1.1 +/- 0.5 and 4.1 +/- 1.8 nM in 0.1 and 1.0 M NaCl, respectively) and inhibition is noncompetitive (k(cat) but not K(m) is affected by the POM concentration).

Inhibition of HIV-1 protease by a boron-modified polypeptide

Six boronated tetrapeptides with the carboxy moiety of phenylalanine replaced by dihydroxyboron were synthesized, and their activities against human immunodeficiency virus 1 (HIV-1) protease subsequently investigated. The sequences of these peptides were derived from HIV-1 protease substrates, which included the C-terminal part of the scissile bond (Phe-Pro) within the gag-pol polyprotein. Enzymatic studies showed that these compounds were competitive inhibitors of HIV-1 protease with K(i) values ranging from 5 to 18 microM when experiments were performed at high enzyme concentrations (above 5 x 10(-8) M); however, at low protease concentrations inhibition was due in part to an increase of the association constants of the protease subunits. Ac-Thr-Leu-Asn-PheB inhibited HIV-1 protease with a K(i) of 5 microM, whereas the non-boronated parental compound was inactive at concentrations up to 400 microM, which indicates the significance of boronation in enzyme inhibition. The boronated tetrapeptides were inhibitory to an HIV-1 protease variant that is resistant to several HIV-1 protease inhibitors. Finally, fluorescence analysis showed that the interactions between the boronated peptide Ac-Thr-Leu-Asn-PheB and HIV-1 protease resulted in a rapid decrease of fluorescence emission at 360 nm, which suggests the formation of a compound/enzyme complex. Boronated peptides may provide useful reagents for studying protease biochemistry and yield valuable information toward the development of protease dimerization inhibitors.

Identification of HIV-1 integrase inhibitors via three-dimensional database searching using ASV and HIV-1 integrases as targets

Integration of viral DNA into the host cell genome is a critical step in the life cycle of HIV. This essential reaction is catalyzed by integrase (IN) through two steps, 3'-processing and DNA strand transfer. Integrase is an attractive target for drug design because there is no known cellular analogue and integration is essential for successful replication of HIV. A computational three-dimensional (3-D) database search was used to identify novel HIV-1 integrase inhibitors. Starting from the previously identified Y3 (4-acetylamino-5-hydroxynaphthalene-2,7-disulfonic acid) binding site on the avian sarcoma virus integrase (ASV IN), a preliminary search of all compounds in the nonproprietary, open part of the National Cancer Institute 3-D database yielded a collection of 3100 compounds. A more rigorous scoring method was used to rescreen the 3100 compounds against both ASV IN and HIV-1 IN. Twenty-two of those compounds were selected for inhibition assays against HIV-1 IN. Thirteen of the 22 showed inhibitory activity against HIV-1 IN at concentrations less than 200 microM and three of them showed antiviral activities in HIV-1 infected CEM cells with effective concentrations (EC50) ranging from 0.8 to 200 microM. Analysis of the computer-generated binding modes of the active compounds to HIV-1 IN showed that simultaneous interaction with the Y3 site and the catalytic site is possible. In addition, interactions between the active compounds and the flexible loop involved in the binding of DNA by IN are indicated to occur. The structural details and the unique binding motif between the HIV-1 IN and its inhibitors identified in the present work may contribute to the future development of IN inhibitors.

HIV protease as a target for retrovirus vector-mediated gene therapy

The dimeric aspartyl protease of HIV has been the subject of intense research for almost a decade. Knowledge of the substrate specificity and catalytic mechanism of this enzyme initially guided the development of several potent peptidomimetic small molecule inhibitors. More recently, the solution of the HIV protease structure led to the structure-based design of improved peptidomimetic and non-peptidomimetic antiviral compounds. Despite the qualified success of these inhibitors, the high mutation rate associated with RNA viruses continues to hamper the long-term clinical efficacy of HIV protease inhibitors. The dimeric nature of the viral protease has been conducive to the investigation of dominant-negative inhibitors of the enzyme. Some of these inhibitors are defective protease monomers that interact with functional monomers to form inactive protease heterodimers. An advantage of macromolecular inhibitors as compared to small-molecule inhibitors is the increased surface area of interaction between the inhibitor and the target gene product. Point mutations that preserve enzyme activity but confer resistance to small-molecule inhibitors are less likely to have an adverse effect on macromolecular interactions. The use of efficient retrovirus vectors has facilitated the delivery of these macromolecular inhibitors to primary human lymphocytes. The vector-transduced cells were less susceptible to HIV infection in vitro, and showed similar levels of protection compared to other macromolecular inhibitors of HIV replication, such as RevM10. These preliminary results encourage the further development of dominant-negative HIV protease inhibitors as a gene therapy-based antiviral strategy.

Isolation of acein-2, a novel angiotensin-I-converting enzyme inhibitory peptide derived from a tryptic hydrolysate of human plasma

We previously described a novel angiotensin-I-converting enzyme (ACE) inhibitory peptide, designated Acein-1, that was isolated from a tryptic hydrolysate of human plasma. We now report a second such inhibitory peptide, Acein-2 obtained from the same hydrolysate. The peptide was purified by gel filtration and cation exchange chromatography followed by reversed-phase gradient and isocratic high performance liquid chromatography. Acein-2 was found to be a tripeptide, Leu-Ile-Tyr, which is thought to correspond to f(518-520) of human alpha2-macroglobulin. The synthetic tripeptide showed a potent dose-dependent inhibition of ACE, with an IC(50) value of 0.82 micromol/l. Lineweaver-Burk plots suggested that Acein-2 as well as the previously described Acein-1 are non-competitive inhibitors.

Structure-based discovery of small molecule inhibitors targeted to protein tyrosine phosphatase 1B

Protein tyrosine phosphatases (PTPases) are involved in the control of tyrosine phosphorylation levels in the cell and are believed to be crucial for the regulation of a multitude of cellular functions. A detailed understanding of the role played by PTPases in various signaling pathways has not yet been achieved, and potent and selective PTPase inhibitors are essential in the quest to determine the functionality of individual PTPases. Using the DOCK methodology, we have carried out a structure-based, computer-assisted search of an available chemical database in order to identify low molecular weight, nonpeptidic PTP1B inhibitors. We have identified several organic molecules that not only possess inhibitory activity against PTP1B but which also display significant selectivity for PTP1B. This indicates that although structural features important for pTyr recognition are conserved among different PTPases, it is possible to generate selective inhibitors targeted primarily to the catalytic site. Kinetic analysis and molecular modeling experiments suggest that the PTP1B active site possesses significant plasticity such that substituted and extended aromatic systems can be accommodated. The newly identified molecules provide a molecular framework upon which therapeutically useful compounds can ultimately be based, and systematic optimization of these lead compounds is likely to further enhance their potency and selectivity.

Structure-based identification of small molecule antiviral compounds targeted to the gp41 core structure of the human immunodeficiency virus type 1

Recent X-ray crystallographic determination of the HIV-1 envelope glycoprotein gp41 core structure opened up a new avenue to discover antiviral agents for chemotherapy of HIV-1 infection and AIDS. We have undertaken a systematic study to search for anti-HIV-1 lead compounds targeted to gp41. Using molecular docking techniques to screen a database of 20 000 organic molecules, we found 16 compounds with the best fit for docking into the hydrophobic cavity within the gp41 core and with maximum possible interactions with the target site. Further testing of these compounds by an enzyme-linked immunosorbent assay and virus inhibition assays discerned two compounds (ADS-J1 and ADS-J2) having inhibitory activity at micromolar concentrations on the formation of the gp41 core structure and on HIV-1 infection. These two compounds will be used as leads to design more effective HIV-1 inhibitors targeted to the HIV-1 gp41 core structure.

A method for biomolecular structural recognition and docking allowing conformational flexibility

In this work, we present an algorithm developed to handle biomolecular structural recognition problems, as part of an interdisciplinary research endeavor of the Computer Vision and Molecular Biology fields. A key problem in rational drug design and in biomolecular structural recognition is the generation of binding modes between two molecules, also known as molecular docking. Geometrical fitness is a necessary condition for molecular interaction. Hence, docking a ligand (e.g., a drug molecule or a protein molecule), to a protein receptor (e.g., enzyme), involves recognition of molecular surfaces. Conformational transitions by "hinge-bending" involves rotational movements of relatively rigid parts with respect to each other. The generation of docked binding modes between two associating molecules depends on their three dimensional structures (3-D) and their conformational flexibility. In comparison to the particular case of rigid-body docking, the computational difficulty grows considerably when taking into account the additional degrees of freedom intrinsic to the flexible molecular docking problem. Previous docking techniques have enabled hinge movements only within small ligands. Partial flexibility in the receptor molecule is enabled by a few techniques. Hinge-bending motions of protein receptors domains are not addressed by these methods, although these types of transitions are significant, e.g., in enzymes activity. Our approach allows hinge induced motions to exist in either the receptor or the ligand molecules of diverse sizes. We allow domains/subdomains/group of atoms movements in either of the associating molecules. We achieve this by adapting a technique developed in Computer Vision and Robotics for the efficient recognition of partially occluded articulated objects. These types of objects consist of rigid parts which are connected by rotary joints (hinges). Our method is based on an extension and generalization of the Hough transform and the Geometric Hashing paradigms for rigid object recognition. We show experimental results obtained by the successful application of the algorithm to cases of bound and unbound molecular complexes, yielding fast matching times. While the "correct" molecular conformations of the known complexes are obtained with small RMS distances, additional, predictive good-fitting binding modes are generated as well. We conclude by discussing the algorithm's implications and extensions, as well as its application to investigations of protein structures in Molecular Biology and recognition problems in Computer Vision.

Predicting relative binding affinities of non-peptide HIV protease inhibitors with free energy perturbation calculations

The relative binding free energies in HIV protease of haloperidol thioketal (THK) and three of its derivatives were examined with free energy calculations. THK is a weak inhibitor (IC50 = 15 microM) for which two cocrystal structures with HIV type 1 proteases have been solved [Rutenber, E. et al., J. Biol. Chem., 268 (1993) 15343]. A THK derivative with a phenyl group on C2 of the piperidine ring was expected to be a poor inhibitor based on experiments with haloperidol ketal and its 2-phenyl derivative (Caldera, P., personal communication). Our calculations predict that a 5-phenyl THK derivative, suggested based on examination of the crystal structure, will bind significantly better than THK. Although there are large error bars as estimated from hysteresis, the calculations predict that the 5-phenyl substituent is clearly favored over the 2-phenyl derivative as well as the parent compound. The unfavorable free energies of solvation of both phenyl THK derivatives relative to the parent compound contributed to their predicted binding free energies. In a third simulation, the change in binding free energy for 5-benzyl THK relative to THK was calculated. Although this derivative has a lower free energy in the protein, its decreased free energy of solvation increases the predicted delta delta G (bind) to the same range as that of the 2-phenyl derivative.

Ligand solvation in molecular docking

Solvation plays an important role in ligand-protein association and has a strong impact on comparisons of binding energies for dissimilar molecules. When databases of such molecules are screened for complementarity to receptors of known structure, as often occurs in structure-based inhibitor discovery, failure to consider ligand solvation often leads to putative ligands that are too highly charged or too large. To correct for the different charge states and sizes of the ligands, we calculated electrostatic and non-polar solvation free energies for molecules in a widely used molecular database, the Available Chemicals Directory (ACD). A modified Born equation treatment was used to calculate the electrostatic component of ligand solvation. The non-polar component of ligand solvation was calculated based on the surface area of the ligand and parameters derived from the hydration energies of apolar ligands. These solvation energies were subtracted from the ligand-receptor interaction energies. We tested the usefulness of these corrections by screening the ACD for molecules that complemented three proteins of known structure, using a molecular docking program. Correcting for ligand solvation improved the rankings of known ligands and discriminated against molecules with inappropriate charge states and sizes.

Acein-1, a novel angiotensin-I-converting enzyme inhibitory peptide isolated from tryptic hydrolysate of human plasma

A novel angiotensin-I-converting enzyme (ACE) inhibitory peptide, designated acein-1, was isolated from the tryptic hydrolysate of human plasma. Gel filtration and cation exchange chromatography were performed to purify this peptide, followed by reversed-phase gradient and isocratic high-performance liquid chromatography. Acein-1 was found to be a heptapeptide, Tyr-Leu-Tyr-Glu-Ile-Ala-Arg, corresponding to f(138-144) of human serum albumin. The synthetic heptapeptide, hexapeptide (Tyr-Leu-Tyr-Glu-Ile-Ala, des-7R acein-1) and octapeptide (Tyr-Leu-Tyr-Glu-Ile-Ala-Arg-Arg, acein-1R) showed dose-dependent inhibitions of ACE, and their IC50 values were 16 micromol/l, 500 micromol/l and 86 micromol/l, respectively. Acein-1 might be a non-competitive inhibitor, while acein-1R may be an uncompetitive inhibitor, as shown by Lineweaver-Burk plots.

Identification of HIV-1 integrase inhibitors based on a four-point pharmacophore

The rapid emergence of human immunodeficiency virus (HIV) strains resistant to available drugs implies that effective treatment modalities will require the use of a combination of drugs targeting different sites of the HIV life cycle. Because the virus cannot replicate without integration into a host chromosome, HIV-1 integrase (IN) is an attractive therapeutic target. Thus, an effective IN inhibitor should provide additional benefit in combination chemotherapy. A four-point pharmacophore has been identified based on the structures of quinalizarin and purpurin, which were found to be potent IN inhibitors using both a preintegration complex assay and a purified enzyme assay in vitro. Searching with this four-point pharmacophore in the 'open' part of the National Cancer Institute three-dimensional structure database produced 234 compounds containing the pharmacophore. Sixty of these compounds were tested for their inhibitory activity against IN using the purified enzyme; 19 were found to be active against IN with IC50 values of less than 100 microM, among which 10 had IC50 values of less than 10 microM. These inhibitors can further serve as leads, and studies are in progress to design novel inhibitors based on the results presented in this study.

Inhibitors of HIV-1 protease: a major success of structure-assisted drug design

Retroviral protease (PR) from the human immunodeficiency virus type 1 (HIV-1) was identified over a decade ago as a potential target for structure-based drug design. This effort was very successful. Four drugs are already approved, and others are undergoing clinical trials. The techniques utilized in this remarkable example of structure-assisted drug design included crystallography, NMR, computational studies, and advanced chemical synthesis. The development of these drugs is discussed in detail. Other approaches to designing HIV-1 PR inhibitors, based on the concepts of symmetry and on the replacement of a water molecule that had been found tetrahedrally coordinated between the enzyme and the inhibitors, are also discussed. The emergence of drug-induced mutations of HIV-1 PR leads to rapid loss of potency of the existing drugs and to the need to continue the development process. The structural basis of drug resistance and the ways of overcoming this phenomenon are mentioned.

Flexible ligand docking using conformational ensembles

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.

Automated docking and the search for HIV protease inhibitors

This article will discuss the motivations, technologies, and future directions of computational automated docking in the context of the structure-based rational design of HIV-1 protease inhibitors. Docking simulations are widely used for screening of compound libraries to identify new drug leads, employing a simple model for rapid testing of thousands of compounds. Docking simulations are also useful for lead enhancement, using more detailed models to analyze the atomic interactions between inhibitors and target macromolecules. Major advances have been reported in the development of empirical force fields, which now allow assessment of relative binding strength and drug specificity, and extensions of automated docking techniques allow de novo drug design.

Entamoeba histolytica: computer-assisted modeling of phosphofructokinase for the prediction of broad-spectrum antiparasitic agents

Pyrophosphate-dependent phosphofructokinase (PPi-PFK) is the rate-limiting glycolytic enzyme found in the pathogenic protists Entamoeba histolytica, Giardia lamblia, Toxoplasma gondii, Trichomonas vaginalis, and Naegleria fowleri. The enzyme differs significantly from ATP-dependent phosphofructokinases found in humans and as such represents an important drug target. Current therapy for infections caused by these pathogens is inadequate, especially for children, pregnant women, and the immune compromised. The development of more selective, safer agents in imperative, as parasitic infections are currently a significant health threat worldwide and will likely become increasingly common agents of disease in the future. For the purpose of designing drugs to treat parasitic infections, we have constructed a model of PPi-PFK from E. histolytica based on the three-dimensional structure of the ATP-dependent PFK from Bacillus stearothermophilus. The model was used with the computer program Dock 3.5 (University of California, San Francisco) to predict the binding of pyrophosphate and selected bisphosphonates to the enzyme. The predicted drug-enzyme interactions suggested that two of these compounds would be competitive inhibitors of pyrophosphate. These drugs were tested against E. histolytica and inhibited the growth of amebae in vitro. This class of compounds may have broad-spectrum antiparasitic activity and, in the future, may facilitate the treatment of serious parasitic infections.

Discovery of HIV-1 integrase inhibitors by pharmacophore searching

Based upon a class of known HIV-1 integrase inhibitors, several pharmacophore models were proposed from molecular modeling studies and validated using a 3D database of 152, compounds for which integrase assay data are known. Using the most probable pharmacophore model as the query, the NCI 3D database of 206,876 compounds was searched, and 340 compounds that contain the pharmacophore query were identified. Twenty-nine of these compounds were selected and tested in the HIV-1 integrase assay. This led to the discovery of 10 novel, structurally diverse HIV-1 integrase inhibitors, four of which have an IC50 value less than 30 microM and are promising lead compounds for further HIV-1 integrase inhibitor development.

Structure-based identification of a ricin inhibitor

Ricin is a potent cytotoxin which has been used widely in the construction of therapeutic agents such as immunotoxins. Recently it has been used by governments and underground groups as a poison. There is interest in identifying and designing effective inhibitors of the ricin A chain (RTA). In this study computer-assisted searches indicated that pterins might bind in the RTA active site which normally recognizes a specific adenine base on rRNA. Kinetic assays showed that pteroic acid could inhibit RTA activity with an apparent Ki of 0.6 mM. A 2.3 A crystal structure of the complex revealed the mode of binding. The pterin ring displaces Tyr80 and binds in the adenine pocket making specific hydrogen bonds to active site residues. The benzoate moiety of pteroic acid binds on the opposite side of Tyr80 making van der Waals contact with the Tyr ring and forming a hydrogen bond with Asn78. Neopterin, a propane triol derivative of pterin, also binds to RTA as revealed by the X-ray structure of its complex with RTA. Neither pterin-6-carboxylic acid nor folic acid bind to the crystal or act as inhibitors. The models observed suggest alterations to the pterin moiety which may produce more potent and specific RTA inhibitors.

Use of the intestinal bile acid transporter for the uptake of cholic acid conjugates with HIV-1 protease inhibitory activity

PURPOSE

To investigate the ability of the human intestinal bile acid transporter to transport cholic acid conjugates with potential HIV-1 protease inhibitory activity.

METHODS

Cholic acid was conjugated at the 24 position of the sterol nucleus with various amino acids and amino acid analogs. The CaCo-2 cell line was used as a model to investigate the interaction of these bile acid conjugates with the human intestinal bile acid transporter. Interaction between the carrier and the conjugates was quantified by inhibition of taurocholic acid transport and confirmed by transport of radiolabelled conjugates in this cell line.

RESULTS

The highest interaction with the transporter, as quantified by inhibition of taurocholic acid transport, occurred when a single negative charge was present around the 24 to 29 region of the sterol nucleus. A second negative charge or a positive charge significantly reduced the interaction. Transport of radiolabelled cholyl-L-Lys-epsilon-tBOC ester and cholyl-D-Asp-beta-benzyl ester was inhibited by taurocholic acid. Of all tested compounds, only cholyl-D-Asp-beta-benzyl ester showed modest HIV-1 protease inhibitory activity with an IC50 of 125 microM.

CONCLUSIONS

Cholic acid-amino acid conjugates with appropriate stereochemistry are recognized and transported by the human bile acid transporter and show modest HIV-1 protease inhibitory activity. Transport of these conjugates by the bile acid carrier is influenced by charge and hydrophobicity around the 24 position of the sterol nucleus.

Analysis of the structure of chemically synthesized HIV-1 protease complexed with a hexapeptide inhibitor. Part I: Crystallographic refinement of 2 A data

The structure of a complex between a hexapeptide-based inhibitor, MVT-101, and the chemically synthesized (Aba 67,95,167,195; Aba: L-alpha-amino-n-butyric acid) protease from the human immunodeficiency virus (HIV-1), reported previously at 2.3 A has now been refined to a crystallographic R factor of 15.4% at 2.0 A resolution. Root mean square deviations from ideality are 0.18 A for bond lengths and 2.4 degrees for the angles. The inhibitor can be fitted to the difference electron density map in two alternative orientations. Drastic differences are observed for positions and interactions at P3/S3 and P3'/S3' subsites of the two orientations due to different crystallographic environments.

Structure-based drug design: computational advances

Structure-based computational methods continue to enhance progress in the discovery and refinement of therapeutic agents. Several such methods and their applications are described. These include molecular visualization and molecular modeling, docking, fragment methods, 3-D database techniques, and free-energy perturbation. Related issues that are discussed include the use of simplified potential energy functions and the determination of the positions of tightly bound waters. Strengths and weaknesses of the various methods are described.

Bile pigments as HIV-1 protease inhibitors and their effects on HIV-1 viral maturation and infectivity in vitro

Using recently developed molecular-shape description algorithms, we searched the Available Chemical Directory for known compounds similar in shape to the potent HIV-1 protease inhibitor Merck L-700,417; 15 compounds most similar in shape to the inhibitor were selected for testing in vitro. Four of these inhibited the protease at 100 microM or less and the most active of the four were the naturally occurring pigments biliverdin and bilirubin. Biliverdin and bilirubin inhibited recombinant HIV-1 protease in vitro at pH 7.8 with K1 values of approx. 1 microM, and also inhibited HIV-2 and simian immunodeficiency virus proteases. The related pyrrolic pigments stercobilin, urobilin, biliverdin dimethyl ester and xanthobilirubic acid showed similar inhibitory activity at low micromolar concentrations. Biliverdin, bilirubin and xanthobilirubic acid did not inhibit viral polyprotein processing in cultured cells, but they reduced viral infectivity significantly. At 100 microM, xanthobilirubic acid affected viral assembly, resulting in a 50% decrease in the generation of infectious particles. In contrast, at the same concentrations biliverdin and bilirubin exerted little or no effect on viral assembly but blocked infection of HeLaT4 cells by 50%. These results suggest that bile pigments might be a new class of potential lead compounds for developing protease inhibitors and they raise the question of whether hyperbilirubinaemia can influence the course of HIV infection.

A new class of HIV-1 protease inhibitor: the crystallographic structure, inhibition and chemical synthesis of an aminimide peptide isostere

The essential role of HIV-1 protease (HIV-1 PR) in the viral life cycle makes it an attractive target for the development of substrate-based inhibitors that may find efficacy as anti-AIDS drugs. However, resistance has arisen to potent peptidomimetic drugs necessitating the further development of novel chemical backbones for diversity based chemistry focused on probing the active site for inhibitor interactions and binding modes that evade protease resistance. AQ148 is a potent inhibitor of HIV-1 PR and represents a new class of transition state analogues incorporating an aminimide peptide isostere. A 3-D crystallographic structure of AQ148, a tetrapeptide isostere, has been determined in complex with its target HIV-1 PR to a resolution of 2.5 A and used to evaluate the specific structural determinants of AQ148 potency and to correlate structure-activity relationships within the class of related compounds. AQ148 is a competitive inhibitor of HIV-1 PR with a Ki value of 137 nM. Twenty-nine derivatives have been synthesized and chemical modifications have been made at the P1, P2, P1', and P2' sites. The atomic resolution structure of AQ148 bound to HIV-1 PR reveals both an inhibitor binding mode that closely resembles that of other peptidomimetic inhibitors and specific protein/inhibitor interactions that correlate with structure-activity relationships. The structure provides the basis for the design, synthesis and evaluation of the next generation of hydroxyethyl aminimide inhibitors. The aminimide peptide isostere is a scaffold with favorable biological properties well suited to both the combinatorial methods of peptidomimesis and the rational design of potent and specific substrate-based analogues.