Re-emerging Aspartic Protease Targets: Examining Cryptococcus neoformans Major Aspartyl Peptidase 1 as a Target for Antifungal Drug Discovery.

Cryptococcosis is an invasive infection that accounts for 15% of AIDS-related fatalities. Still, treating cryptococcosis remains a significant challenge due to the poor availability of effective antifungal therapies and emergence of drug resistance. Interestingly, protease inhibitor components of antiretroviral therapy regimens have shown some clinical benefits in these opportunistic infections. We investigated Major aspartyl peptidase 1 (May1), a secreted Cryptococcus neoformans protease, as a possible target for the development of drugs that act against both fungal and retroviral aspartyl proteases. Here, we describe the biochemical characterization of May1, present its high-resolution X-ray structure, and provide its substrate specificity analysis. Through combinatorial screening of 11,520 compounds, we identified a potent inhibitor of May1 and HIV protease. This dual-specificity inhibitor exhibits antifungal activity in yeast culture, low cytotoxicity, and low off-target activity against host proteases and could thus serve as a lead compound for further development of May1 and HIV protease inhibitors.

Characterization of Nine Compounds Isolated from the Acid Hydrolysate of Lonicera fulvotomentosa Hsu et S. C. Cheng and Evaluation of Their In Vitro Activity towards HIV Protease.

In this study, we isolated nine compounds from the acid hydrolysate of the flower buds of Lonicera fulvotomentosa Hsu et S. C. Cheng and characterized their chemical structures using 1H-NMR, 13C-NMR, and electron ionization mass spectroscopy (EI-MS). These compounds were identified as ß-sitosterol (1), 5,5'-dibutoxy-2,2'-bifuran (2), nonacosane-10-ol (3), ethyl (3ß)-3,23-dihydroxyolean-12-en-28-oate (4), oleanolic acid (5), ethyl caffeate (6), caffeic acid (7), isovanillin (8), and hederagenin (9), with 4 as a new triterpene compound. Inhibitory activity against human immunodeficiency virus (HIV) protease was also evaluated for the compounds, and only ethyl caffeate, caffeic acid, and isovanillin (6, 7, and 8) exhibited inhibitory effects, with IC50 values of 1.0 µM, 1.5 µM, and 3.5 µM, respectively. Molecular docking with energy minimization and subsequent molecular dynamic (MD) simulation showed that ethyl caffeate and caffeic acid bound to the active site of HIV protease, while isovanillin drifted out from the active site and dissociated into bulk water during MD simulations, and most of the binding residues of HIV protease have been previously identified for HIV protease inhibitors. These results suggest that caffeic acid derivatives may possess inhibitory activities towards HIV protease other than previously reported inhibitory activities against HIV integrase, and thus ethyl caffeate and caffeic acid could be used as lead compounds in developing potential HIV protease inhibitors, and possibly even dual-function inhibitors against HIV.

Novel Central Nervous System (CNS)-Targeting Protease Inhibitors for Drug-Resistant HIV Infection and HIV-Associated CNS Complications.

There is currently no specific therapeutics for the HIV-1-related central nervous system (CNS) complications. Here we report that three newly designed CNS-targeting HIV-1 protease inhibitors (PIs), GRL-083-13, GRL-084-13, and GRL-087-13, which contain a P1-3,5-bis-fluorophenyl or P1-para-monofluorophenyl ring, and P2-bis-tetrahydrofuran (bis-THF) or P2-tetrahydropyrano-tetrahydrofuran (Tp-THF), with a sulfonamide isostere, are highly active against wild-type HIV-1 strains and primary clinical isolates (50% effective concentration [EC50], 0.0002 to ∼0.003 µM), with minimal cytotoxicity. These CNS-targeting PIs efficiently suppressed the replication of HIV-1 variants (EC50, 0.002 to ∼0.047 µM) that had been selected to propagate at high concentrations of conventional HIV-1 PIs. Such CNS-targeting PIs maintained their antiviral activity against HIV-2ROD as well as multidrug-resistant clinical HIV-1 variants isolated from AIDS patients who no longer responded to existing antiviral regimens after long-term therapy. Long-term drug selection experiments revealed that the emergence of resistant-HIV-1 against these CNS-targeting PIs was substantially delayed. In addition, the CNS-targeting PIs showed the most favorable CNS penetration properties among the tested compounds, including various FDA-approved anti-HIV-1 drugs, as assessed with the in vitro blood-brain barrier reconstruction system. Crystallographic analysis demonstrated that the bicyclic rings at the P2 moiety of the CNS-targeting PIs form strong hydrogen-bond interactions with HIV-1 protease (PR) active site. Moreover, both the P1-3,5-bis-fluorophenyl and P1-para-monofluorophenyl rings sustain greater van der Waals contacts with PR than in the case of darunavir (DRV). The data suggest that the present CNS-targeting PIs have desirable features for treating patients infected with wild-type and/or multidrug-resistant HIV-1 strains and might serve as promising preventive and/or therapeutic candidates for HIV-1-associated neurocognitive disorders (HAND) and other CNS complications.

HIV-1 Protease Uses Bi-Specific S2/S2' Subsites to Optimize Cleavage of Two Classes of Target Sites.

Retroviral proteases (PRs) have a unique specificity that allows cleavage of sites with or without a P1' proline. A P1' proline is required at the MA/CA cleavage site due to its role in a post-cleavage conformational change in the capsid protein. However, the HIV-1 PR prefers to have large hydrophobic amino acids flanking the scissile bond, suggesting that PR recognizes two different classes of substrate sequences. We analyzed the cleavage rate of over 150 combinations of six different HIV-1 cleavage sites to explore rate determinants of cleavage. We found that cleavage rates are strongly influenced by the two amino acids flanking the amino acids at the scissile bond (P2-P1/P1'-P2'), with two complementary sets of rules. When P1' is proline, the P2 side chain interacts with a polar region in the S2 subsite of the PR, while the P2' amino acid interacts with a hydrophobic region of the S2' subsite. When P1' is not proline, the orientations of the P2 and P2' side chains with respect to the scissile bond are reversed; P2 residues interact with a hydrophobic face of the S2 subsite, while the P2' amino acid usually engages hydrophilic amino acids in the S2' subsite. These results reveal that the HIV-1 PR has evolved bi-functional S2 and S2' subsites to accommodate the steric effects imposed by a P1' proline on the orientation of P2 and P2' substrate side chains. These results also suggest a new strategy for inhibitor design to engage the multiple specificities in these subsites.

The inhibition effect of garlic-derived compounds on human immunodeficiency virus type 1 and saquinavir.

Garlic has been used as a traditional medicine to treat various diseases. Garlic reduces the risk of some diseases. This protective effect is due to the organosulfur compounds of garlic. The aim of this study was to investigate the inhibition effects of garlic-derived compounds on human immunodeficiency virus type 1 (HIV-1) and as the most important anti-HIV-1 medicine. The activation of saquinavir is believed to be the principal mechanism behind the protective effects of HIV-1. Our theoretical calculations are performed for blood phase by using the density functional theory for the main compounds of garlic. The chemical activity and solubility of ajoene and the mainly derived compounds of garlic as theoretical calculations are important for the medical research comparing with the other compounds of the garlic. The theoretical calculations have helped us to determine which active ingredient of the garlic having inhibition effects on HIV-1 and saquinavir.

Polygonumins A, a newly isolated compound from the stem of Polygonum minus Huds with potential medicinal activities.

Polygonumins A, a new compound, was isolated from the stem of Polygonum minus. Based on NMR results, the compound's structure is identical to that of vanicoside A, comprising four phenylpropanoid ester units and a sucrose unit. The structure differences were located at C-3″″'. The cytotoxic activity of polygonumins A was evaluated on several cancer cell lines by a cell viability assay using tetrazolium dye 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT). The compound showed the highest antiproliferative (p < 0.05) activities against K562 (Human Leukaemia Cell Line), MCF7 (Human breast adenocarcinoma cell line), and HCT116 (Colorectal cancer cells) cells. Cytotoxic studies against V79-4 cells were carried out and showed that polygonumins A was toxic at 50 µg/ml, suggesting that this compound may be used as an anticancer drug without affecting normal cells. Polygonumins A also showed promising activity as an HIV-1 protease inhibitor with 56% relative inhibition. Molecular docking results indicated that the compound possesses high binding affinity towards the HIV protease over the low binding free energy range of -10.5 to -11.3 kcal/mol. P. minus is used in Malaysian traditional medicine for the treatment of tumour cells. This is the first report on the use of P. minus as an HIV-1 protease inhibitor.

A Modified P1 Moiety Enhances In Vitro Antiviral Activity against Various Multidrug-Resistant HIV-1 Variants and In Vitro Central Nervous System Penetration Properties of a Novel Nonpeptidic Protease Inhibitor, GRL-10413.

We report here that GRL-10413, a novel nonpeptidic HIV-1 protease inhibitor (PI) containing a modified P1 moiety and a hydroxyethylamine sulfonamide isostere, is highly active against laboratory HIV-1 strains and primary clinical isolates (50% effective concentration [EC50] of 0.00035 to 0.0018 µM), with minimal cytotoxicity (50% cytotoxic concentration [CC50] = 35.7 µM). GRL-10413 blocked the infectivity and replication of HIV-1NL4-3 variants selected by use of atazanavir, lopinavir, or amprenavir (APV) at concentrations of up to 5 µM (EC50 = 0.0021 to 0.0023 µM). GRL-10413 also maintained its strong antiviral activity against multidrug-resistant clinical HIV-1 variants isolated from patients who no longer responded to various antiviral regimens after long-term antiretroviral therapy. The development of resistance against GRL-10413 was significantly delayed compared to that against APV. In addition, GRL-10413 showed favorable central nervous system (CNS) penetration properties as assessed with an in vitro blood-brain barrier (BBB) reconstruction system. Analysis of the crystal structure of HIV-1 protease in complex with GRL-10413 demonstrated that the modified P1 moiety of GRL-10413 has a greater hydrophobic surface area and makes greater van der Waals contacts with active site amino acids of protease than in the case of darunavir. Moreover, the chlorine substituent in the P1 moiety interacts with protease in two distinct configurations. The present data demonstrate that GRL-10413 has desirable features for treating patients infected with wild-type and/or multidrug-resistant HIV-1 variants, with favorable CNS penetration capability, and that the newly modified P1 moiety may confer desirable features in designing novel anti-HIV-1 PIs.

Enthalpy screen of drug candidates.

The enthalpic and entropic contributions to the binding affinity of drug candidates have been acknowledged to be important determinants of the quality of a drug molecule. These quantities, usually summarized in the thermodynamic signature, provide a rapid assessment of the forces that drive the binding of a ligand. Having access to the thermodynamic signature in the early stages of the drug discovery process will provide critical information towards the selection of the best drug candidates for development. In this paper, the Enthalpy Screen technique is presented. The enthalpy screen allows fast and accurate determination of the binding enthalpy for hundreds of ligands. As such, it appears to be ideally suited to aid in the ranking of the hundreds of hits that are usually identified after standard high throughput screening.

Computational Prediction of HIV-1 Resistance to Protease Inhibitors.

The development of mutations in HIV-1 protease (PR) hinders the activity of antiretroviral drugs, forcing changes in drug prescription. Most resistance assessments used to date rely on expert-based rules on predefined sets of stereotypical mutations; such an information-driven approach cannot capture new polymorphisms or be applied for new drugs. Computational modeling could provide a more general assessment of drug resistance and could be made available to clinicians through the Internet. We have created a protocol involving sequence comparison and all-atom protein-ligand induced fit simulations to predict resistance at the molecular level. We first compared our predictions with the experimentally determined IC50 values of darunavir, amprenavir, ritonavir, and indinavir from reference PR mutants displaying different resistance levels. We then performed analyses on a large set of variants harboring more than 10 mutations. Finally, several sequences from real patients were analyzed for amprenavir and darunavir. Our computational approach detected all of the genotype changes triggering high-level resistance, even those involving a large number of mutations.

Using Hierarchical Virtual Screening To Combat Drug Resistance of the HIV-1 Protease.

Human immunodeficiency virus (HIV) protease inhibitors (PIs) are important components of highly active anti-retroviral therapy (HAART) that block the catalytic site of HIV protease, thus preventing maturation of the HIV virion. However, with two decades of PI prescriptions in clinical practice, drug-resistant HIV mutants have now been found for all of the PI drugs. Therefore, the continuous development of new PI drugs is crucial both to combat the existing drug-resistant HIV strains and to provide treatments for future patients. Here we purpose an HIV PI drug design strategy to select candidate PIs with binding energy distributions dominated by interactions with conserved protease residues in both wild-type and various drug-resistant mutants. On the basis of this strategy, we have constructed a virtual screening pipeline including combinatorial library construction, combinatorial docking, MM/GBSA-based rescoring, and reranking on the basis of the binding energy distribution. We have tested our strategy on lopinavir by modifying its two functional groups. From an initial 751 689 candidate molecules, 18 candidate inhibitors were selected using the pipeline for experimental validation. IC50 measurements and drug resistance predictions successfully identified two ligands with both HIV protease inhibitor activity and an improved drug resistance profile on 2382 HIV mutants. This study provides a proof of concept for the integration of MM/GBSA energy analysis and drug resistance information at the stage of virtual screening and sheds light on future HIV drug design and the use of virtual screening to combat drug resistance.

Design of gem-difluoro-bis-tetrahydrofuran as P2 ligand for HIV-1 protease inhibitors to improve brain penetration: synthesis, X-ray studies, and biological evaluation.

The structure-based design, synthesis, biological evaluation, and X-ray structural studies of fluorine-containing HIV-1 protease inhibitors are described. The synthesis of both enantiomers of the gem-difluoro-bis-THF ligands was carried out in a stereoselective manner using a Reformatskii-Claisen reaction as the key step. Optically active ligands were converted into protease inhibitors. Two of these inhibitors, (3R,3aS,6aS)-4,4-difluorohexahydrofuro[2,3-b]furan-3-yl(2S,3R)-3-hydroxy-4-((N-isobutyl-4-methoxyphenyl)sulfonamido)-1-phenylbutan-2-yl) carbamate (3) and (3R,3aS,6aS)-4,4-difluorohexahydrofuro[2,3-b]furan-3-yl(2S,3R)-3-hydroxy-4-((N-isobutyl-4-aminophenyl)sulfonamido)phenylbutan-2-yl) carbamate (4), exhibited HIV-1 protease inhibitory Ki values in the picomolar range. Both 3 and 4 showed very potent antiviral activity, with respective EC50 values of 0.8 and 3.1 nM against the laboratory strain HIV-1LAI . The two inhibitors exhibited better lipophilicity profiles than darunavir, and also showed much improved blood-brain barrier permeability in an in vitro model. A high-resolution X-ray structure of inhibitor 4 in complex with HIV-1 protease was determined, revealing that the fluorinated ligand makes extensive interactions with the S2 subsite of HIV-1 protease, including hydrogen bonding interactions with the protease backbone atoms. Moreover, both fluorine atoms on the bis-THF ligand formed strong interactions with the flap Gly 48 carbonyl oxygen atom.

Identifying putative drug targets and potential drug leads: starting points for virtual screening and docking.

The availability of 3D models of both drug leads (small molecule ligands) and drug targets (proteins) is essential to molecular docking and computational drug discovery. This chapter describes a simple approach that can be used to identify both drug leads and drug targets using two popular Web-accessible databases: (1) DrugBank and (2) The Human Metabolome Database. First, it is illustrated how putative drug targets and drug leads for exogenous diseases (i.e., infectious diseases) can be readily identified and their 3D structures selected using only the genomic sequences from pathogenic bacteria or viruses as input. The second part illustrates how putative drug targets and drug leads for endogenous diseases (i.e., noninfectious diseases or chronic conditions) can be identified using similar databases and similar sequence input. This chapter is intended to illustrate how bioinformatics and cheminformatics can work synergistically to help provide the necessary inputs for computer-aided drug design.

Distinguishing binders from false positives by free energy calculations: fragment screening against the flap site of HIV protease.

Molecular docking is a powerful tool used in drug discovery and structural biology for predicting the structures of ligand-receptor complexes. However, the accuracy of docking calculations can be limited by factors such as the neglect of protein reorganization in the scoring function; as a result, ligand screening can produce a high rate of false positive hits. Although absolute binding free energy methods still have difficulty in accurately rank-ordering binders, we believe that they can be fruitfully employed to distinguish binders from nonbinders and reduce the false positive rate. Here we study a set of ligands that dock favorably to a newly discovered, potentially allosteric site on the flap of HIV-1 protease. Fragment binding to this site stabilizes a closed form of protease, which could be exploited for the design of allosteric inhibitors. Twenty-three top-ranked protein-ligand complexes from AutoDock were subject to the free energy screening using two methods, the recently developed binding energy analysis method (BEDAM) and the standard double decoupling method (DDM). Free energy calculations correctly identified most of the false positives (≥83%) and recovered all the confirmed binders. The results show a gap averaging ≥3.7 kcal/mol, separating the binders and the false positives. We present a formula that decomposes the binding free energy into contributions from the receptor conformational macrostates, which provides insights into the roles of different binding modes. Our binding free energy component analysis further suggests that improving the treatment for the desolvation penalty associated with the unfulfilled polar groups could reduce the rate of false positive hits in docking. The current study demonstrates that the combination of docking with free energy methods can be very useful for more accurate ligand screening against valuable drug targets.

An infrastructure to mine molecular descriptors for ligand selection on virtual screening.

The receptor-ligand interaction evaluation is one important step in rational drug design. The databases that provide the structures of the ligands are growing on a daily basis. This makes it impossible to test all the ligands for a target receptor. Hence, a ligand selection before testing the ligands is needed. One possible approach is to evaluate a set of molecular descriptors. With the aim of describing the characteristics of promising compounds for a specific receptor we introduce a data warehouse-based infrastructure to mine molecular descriptors for virtual screening (VS). We performed experiments that consider as target the receptor HIV-1 protease and different compounds for this protein. A set of 9 molecular descriptors are taken as the predictive attributes and the free energy of binding is taken as a target attribute. By applying the J48 algorithm over the data we obtain decision tree models that achieved up to 84% of accuracy. The models indicate which molecular descriptors and their respective values are relevant to influence good FEB results. Using their rules we performed ligand selection on ZINC database. Our results show important reduction in ligands selection to be applied in VS experiments; for instance, the best selection model picked only 0.21% of the total amount of drug-like ligands.

Exploring the antiviral activity of juglone by computational method.

Nature has been the best source of medicines for a long time. Many plant extracts have been used as drugs. Juglone occurs in all parts of the Juglandaceae family and is found extensively in the black walnut plants. It possesses antifungal, antimalarial, antibacterial and antiviral properties besides exhibiting cytotoxic effects. Juglone has gained interest by the researchers for its anticancer properties. This article elucidates the antiviral activity of the Juglone by the computational method.

Crystallographic fragment-based drug discovery: use of a brominated fragment library targeting HIV protease.

A library of 68 brominated fragments was screened against a new crystal form of inhibited HIV-1 protease in order to probe surface sites in soaking experiments. Often, fragments are weak binders with partial occupancy, resulting in weak, difficult-to-fit electron density. The use of a brominated fragment library addresses this challenge, as bromine can be located unequivocally via anomalous scattering. Data collection was carried out in an automated fashion using AutoDrug at SSRL. Novel hits were identified in the known surface sites: 3-bromo-2,6-dimethoxybenzoic acid (Br6) in the flap site and 1-bromo-2-naphthoic acid (Br27) in the exosite, expanding the chemistry of known fragments for development of higher affinity potential allosteric inhibitors. At the same time, mapping the binding sites of a number of weaker binding Br-fragments provides further insight into the nature of these surface pockets.

Characterization of small molecule binding. I. Accurate identification of strong inhibitors in virtual screening.

Accurately ranking docking poses remains a great challenge in computer-aided drug design. In this study, we present an integrated approach called MIEC-SVM that combines structure modeling and statistical learning to characterize protein-ligand binding based on the complex structure generated from docking. Using the HIV-1 protease as a model system, we showed that MIEC-SVM can successfully rank the docking poses and consistently outperformed the state-of-art scoring functions when the true positives only account for 1% or 0.5% of all the compounds under consideration. More excitingly, we found that MIEC-SVM can achieve a significant enrichment in virtual screening even when trained on a set of known inhibitors as small as 50, especially when enhanced by a model average approach. Given these features of MIEC-SVM, we believe it provides a powerful tool for searching for and designing new drugs.

Virtual screening for HIV protease inhibitors: a comparison of AutoDock 4 and Vina.

BACKGROUND: The AutoDock family of software has been widely used in protein-ligand docking research. This study compares AutoDock 4 and AutoDock Vina in the context of virtual screening by using these programs to select compounds active against HIV protease. METHODOLOGY/PRINCIPAL FINDINGS: Both programs were used to rank the members of two chemical libraries, each containing experimentally verified binders to HIV protease. In the case of the NCI Diversity Set II, both AutoDock 4 and Vina were able to select active compounds significantly better than random (AUC = 0.69 and 0.68, respectively; p<0.001). The binding energy predictions were highly correlated in this case, with r = 0.63 and iota = 0.82. For a set of larger, more flexible compounds from the Directory of Universal Decoys, the binding energy predictions were not correlated, and only Vina was able to rank compounds significantly better than random. CONCLUSIONS/SIGNIFICANCE: In ranking smaller molecules with few rotatable bonds, AutoDock 4 and Vina were equally capable, though both exhibited a size-related bias in scoring. However, as Vina executes more quickly and is able to more accurately rank larger molecules, researchers should look to it first when undertaking a virtual screen.

Drug candidates from traditional chinese medicines.

Good progress has been made to modernize traditional Chinese medicines by obtaining active components from natural herbs. In this review, some recent works on procuring active components and modernizing traditional Chinese medicines will be covered. In addition, some recent works on drug design using modern drug design tools have been described. With some well defined targets, the traditional Chinese medicine databases have been screened so as to identify those compounds for which the potential as a drug candidate was not known before. Among these studies, two have been selected as examples to be discussed in details. First, new anti-HIV candidates have been detected, namely leucovorin and agaritine derivatives. Subsequently, GTS-21 is proved to be a good candidate for Alzheimer's disease. All these findings may provide useful information for finding effective drug candidates with lower cost.

Inorganic polyhedral metallacarborane inhibitors of HIV protease: a new approach to overcoming antiviral resistance.

HIV protease (PR) is a prime target for rational anti-HIV drug design. We have previously identified icosahedral metallacarboranes as a novel class of nonpeptidic protease inhibitors. Now we show that substituted metallacarboranes are potent and specific competitive inhibitors of drug-resistant HIV PRs prepared either by site-directed mutagenesis or cloned from HIV-positive patients. Molecular modeling explains the inhibition profile of metallacarboranes by their unconventional binding mode.