Machine learning-based predictive models for identifying high active compounds against HIV-1 integrase

HIV-integrase is an important drug target because it catalyzes chromosomal integration of proviral DNA towards establishing latent infection. Computer-aided drug design has immensely contributed to identifying and developing novel antiviral drugs. We have developed various machine learning-based predictive models for identifying high activity compounds against HIV-integrase. Multiclass models were built using support vector machine with reasonable accuracy on the test and evaluation sets. The developed models were evaluated by rigorous validation approaches and the best features were selected by Boruta method. As compared to the model developed from all descriptors set, a slight improvement was observed among the selected descriptors. Validated models were further used for virtual screening of potential compounds from ChemBridge library. Of the six high active compounds predicted from selected models, compounds 9103124, 6642917 and 9082952 showed the most reasonable binding-affinity and stable-interaction with HIV-integrase active-site residues Asp64, Glu152 and Asn155. This was in agreement with previous reports on the essentiality of these residues against a wide range of inhibitors. We therefore highlight the rigorosity of validated classification models for accurate prediction and ranking of high active lead drugs against HIV-integrase.

Design, Synthesis, Molecular Modeling and Anti-HIV Assay of Novel Quinazolinone Incorporated Coumarin Derivatives

BACKGROUND

The emergence of drug-resistant viral strains has created the need for the development of novel anti-HIV agents with a diverse structure that targets key enzymes in the HIV lifecycle.

OBJECTIVE

Considering the pharmacophore of integrase inhibitors, one of the validated targets for anti-HIV therapy, we designed a quinazolinone incorporated coumarin scaffold to affect HIV.

METHODS

Coumarin is a beta enol ester and also a well-known drug scaffold. Designed structures were prepared using a one-pot three-component reaction from 3-amino-4-hydroxycoumarin, isatoic anhydride and benzaldehyde derivatives.

RESULTS

In vitro anti-HIV and cytotoxicity assay indicated that more than half of the compounds had EC50 values lower than 50 µM. Unsubstituted phenyl derivative showed the highest activity and selectivity with an EC50 value of 5 µM and a therapeutic index of 7. Compounds were docked into the integrase active site to investigate the probable mechanism of action. Accordingly, the hydroxyl moiety of coumarin along with the carbonyl of the quinazolinone ring could function as the metal chelating group. Quinazolinone and phenyl groups interact with side chains of IN residues, as well.

CONCLUSION

Here, a novel anti-HIV scaffold is represented for further modification and in-vivo studies.

Recent advances in the discovery of small-molecule inhibitors of HIV-1 integrase

AIDS caused by the infection of HIV is a prevalent problem today. Rapid development of drug resistance to existing drug classes has called for the discovery of new targets. Within the three major enzymes (i.e., HIV-1 protease, HIV-1 reverse transcriptase and HIV-1 integrase [IN]) of the viral replication cycle, HIV-1 IN has been of particular interest due to the absence of human cellular homolog. HIV-1 IN catalyzes the integration of viral genetic material with the host genome, a key step in the viral replication process. Several novel classes of HIV IN inhibitors have been explored by targeting different sites on the enzyme. This review strives to provide readers with updates on the recent developments of HIV-1 IN inhibitors.

AIDS is an epidemic disease that endangers the lives of millions of people across the world. The AIDS virus, also known as HIV, has developed resistance to the majority of available drugs on the market, thus requiring the need for new drugs. HIV integrase is one of the key viral enzymes required for viral cell proliferation. Since there is no similar enzyme in the human body, major emphasis is being made to develop therapeutics for this novel target. The drugs that are at various stages of development for this target are reviewed here.

The advancement of multidimensional QSAR for novel drug discovery - where are we headed?

INTRODUCTION

The Multidimensional quantitative structure-activity relationship (multidimensional-QSAR) method is one of the most popular computational methods employed to predict interesting biochemical properties of existing or hypothetical molecules. With continuous progress, the QSAR method has made remarkable success in various fields, such as medicinal chemistry, material science and predictive toxicology. Areas covered: In this review, the authors cover the basic elements of multidimensional -QSAR including model construction, validation and application. It includes and emphasizes the very recent developments of multidimensional -QSAR such as: HQSAR, G-QSAR, MIA-QSAR, multi-target QSAR. The advantages and disadvantages of each method are also discussed and typical examples of their application are detailed. Expert opinion: Although there are defects in multidimensional-QSAR modeling, it is still of enormous help to chemists, biologists and other researchers in various fields. In the authors' opinion, the latest more precise and feasible QSAR models should be further developed by integrating new descriptors, algorithms and other relevant computational techniques. Apart from being applied in traditional fields (e.g. lead optimization and predictive risk assessment), QSAR should be used more widely as a routine method in other emerging research fields including the modeling of nanoparticles(NPs), mixture toxicity and peptides.

Targeting HIV-1 integrase with strand transfer inhibitors

HIV-1 integrase (IN) is a retroviral enzyme essential for integration of genetic material into the DNA of the host cell and hence for viral replication. The absence of an equivalent enzyme in humans makes IN an interesting target for anti-HIV drug design. This review briefly overviews the structural and functional properties of HIV-1 IN. We analyze the binding modes of the established drugs, clinical candidates and a comprehensive library of leads based on innovative chemical scaffolds of HIV-1 IN strand transfer inhibitors (INSTIs). Computational clustering techniques are applied for identifying structural features relating to bioactivity. From bio- and chemo-informatics analyses, we provide novel insights into structure-activity relationships of INSTIs and elaborate new strategies for design of innovative inhibitors.

The Influence of Bioisosteres in Drug Design: Tactical Applications to Address Developability Problems

The application of bioisosteres in drug discovery is a well-established design concept that has demonstrated utility as an approach to solving a range of problems that affect candidate optimization, progression, and durability. In this chapter, the application of isosteric substitution is explored in a fashion that focuses on the development of practical solutions to problems that are encountered in typical optimization campaigns. The role of bioisosteres to affect intrinsic potency and selectivity, influence conformation, solve problems associated with drug developability, including P-glycoprotein recognition, modulating basicity, solubility, and lipophilicity, and to address issues associated with metabolism and toxicity is used as the underlying theme to capture a spectrum of creative applications of structural emulation in the design of drug candidates.

Docking-based CoMFA and CoMSIA study of azaindole carboxylic acid derivatives as promising HIV-1 integrase inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were performed based on a series of azaindole carboxylic acid derivatives that had previously been reported as promising HIV-1 integrase inhibitors. Docking studies to explore the binding mode were performed based on the highly active molecule 36. The best docked conformation of molecule 36 was used as template for alignment. The comparative molecular field analysis (CoMFA) model (including steric and electrostatic fields) yielded the cross validation q (2) = 0.655, non-cross validation r (2) = 0.989 and predictive r (2) pred = 0.979. The best comparative molecular similarity indices analysis (CoMSIA) model (including steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields) yielded the cross validation q (2) = 0.719, non-cross validation r (2) = 0.992 and predictive r (2) pred = 0.953. A series of new azaindole carboxylic acid derivatives were designed and the HIV-1 integrase inhibitory activities of these designed compounds were predicted based on the CoMFA and CoMSIA models.

Design and synthesis of novel pyrimidone analogues as HIV-1 integrase inhibitors

A series of novel pyrimidone analogues have been designed and synthesized as HIV-1 integrase (IN) inhibitors. This study demonstrated that introducing a substituent in the N1-position of the pyrimidone scaffold does not significantly influence IN inhibitory activity. Molecular docking studies showed these compounds could occupy the IN active site and form pi-pi interactions with viral DNA nucleotides DC16 and DA17 to displace reactive viral DNA 3'OH and block intasome activity.

New scaffolds of natural origin as Integrase-LEDGF/p75 interaction inhibitors: virtual screening and activity assays

The disruption of crucial interactions between HIV-1 Integrase and cellular cofactor LEDGF/p75 represents an emerging approach for the design and development of new antiretroviral agents. In this study we report the successful application of a structure-based virtual screening strategy for the discovery of natural hit structures able to inhibit Integrase-LEDGF/p75 interaction. The application of sequential filters (drug-likeness, 3D-pharmacophore mapping, docking, molecular dynamics simulations) yielded a hit list of compounds, out of which 9 were tested in the in vitro AlphaScreen assays and 8 exhibited a detectable inhibition of the interaction between the two proteins. The best inhibitors belong to different chemical classes and could be represent a good starting point for further optimization and structure-activity relationship studies.

4-Substituted 2-Hydroxyisoquinoline-1,3(2H,4H)-diones as a Novel Class of HIV-1 Integrase Inhibitors

A series of 2-hydroxy-1,3-dioxoisoquinoline-4-carboxamides featuring an N-hydroxyimide chelating functionality was evaluated for their inhibitory properties against human immunodeficiency virus type 1 integrase (HIV-1 IN). Several derivatives displayed low nanomolar IC50 values comparable to that of the clinically used raltegravir. A marked effect of one compound on both primary IN-catalyzed reactions, strand transfer (ST), and 3' processing (3'-P), emphasizes a novel IN inhibition mechanism establishing it as a potential new generation IN inhibitor. Substitution of the 2-hydroxyisoquinoline-1,3-dione scaffold at position 4 by carboxamido chains was beneficial for antiviral activity since reproducible low micromolar anti-HIV activities were obtained for the first time within this scaffold.

Identification of potential HIV-1 integrase strand transfer inhibitors: in silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski's rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.

A symmetric region of the HIV-1 integrase dimerization interface is essential for viral replication

HIV-1 integrase (IN) is an important target for contemporary antiretroviral drug design research. Historically, efforts at inactivating the enzyme have focused upon blocking its active site. However, it has become apparent that new classes of allosteric inhibitors will be necessary to advance the antiretroviral field in light of the emergence of viral strains resistant to contemporary clinically used IN drugs. In this study we have characterized the importance of a close network of IN residues, distant from the active site, as important for the obligatory multimerization of the enzyme and viral replication as a whole. Specifically, we have determined that the configuration of six residues within a highly symmetrical region at the IN dimerization interface, composed of a four-tiered aromatic interaction flanked by two salt bridges, significantly contributes to proper HIV-1 replication. Additionally, we have utilized a quantitative luminescence assay to examine IN oligomerization and have determined that there is a very low tolerance for amino acid substitutions along this region. Even conservative residue substitutions negatively impacted IN multimerization, resulting in an inactive viral enzyme and a non-replicative virus. We have shown that there is a very low tolerance for amino acid variation at the symmetrical dimeric interface region characterized in this study, and therefore drugs designed to target the amino acid network detailed here could be expected to yield a significantly reduced number of drug-resistant escape mutations compared to contemporary clinically-evaluated antiretrovirals.

Resistance to integrase inhibitors

Integrase (IN) is a clinically validated target for the treatment of human immunodeficiency virus infections and raltegravir exhibits remarkable clinical activity. The next most advanced IN inhibitor is elvitegravir. However, mutant viruses lead to treatment failure and mutations within the IN coding sequence appear to confer cross-resistance. The characterization of those mutations is critical for the development of second generation IN inhibitors to overcome resistance. This review focuses on IN resistance based on structural and biochemical data, and on the role of the IN flexible loop i.e., between residues G140-G149 in drug action and resistance.

Synthesis and anti-HSV-1 evaluation of new 3H-benzo[b]pyrazolo[3,4-h]-1,6-naphthyridines and 3H-pyrido[2,3-b]pyrazolo[3,4-h]-1,6-naphthyridines

BACKGROUND

Herpes simplex virus type-1 (HSV-1) is the primary cause of facial lesions (mouth, lips, and eyes) in humans. The widespread use of acyclovir and nucleoside analogues has led to emergence of HSV strains that are resistant to these drugs. Recently, non-nucleoside anti-HSV compounds have received considerable attention. 1,6-Naphthyridines are a class of heterocyclic compounds that exhibit a broad spectrum of biological activities such as inhibitor of HIV-1 integrase, HCMV, FGF receptor-1 tyrosine kinase, and the enzyme acetylcholinesterase. We previously reported the synthesis, SAR studies, and evaluation anti-HSV-1 activity of 3H-benzo[b]pyrazolo[3,4-h]-1,6-naphthyridines. In the course of our search for new 1,6-naphthyridines derivatives with potential activity against HSV-1, we have synthesized and evaluated new 3H-benzo[b]pyrazolo[3,4-h]-1,6-naphthyridines (1a-k) and 3H-pyrido[2,3-b]pyrazolo[3,4-h]-1,6-naphthyridines (2a-c).

RESULTS

A known synthetic approach was used for preparing new 3H-benzo[b]pyrazolo[3,4-h]-1,6-naphthyridines (1a-k) and 3H-pyrido[2,3-b]pyrazolo[3,4-h]-1,6-naphthyridines (2a-c), starting from ethyl 4-chloro-1-phenyl-1H-pyrazolo[3,4-b]pyridine-5-carboxylate (7). All compounds were identified by FTIR, 1H NMR, and mass spectrometry. The antiviral effect on HSV-1 virus replication was determined.

CONCLUSIONS

The compounds 1d, 1f, 1g, and 1h exhibited the highest anti-HSV-1 activity. In general, 3H-benzo[b]pyrazolo[3,4-h]-1,6-naphthyridines were more effective inhibitors than their corresponding 3H-pyrido[2,3-b]pyrazolo[3,4-h]-1,6-naphthyridines. The compound 1h reduced the virus yield in 91% at 50 μM and exhibited a low cytotoxicity (CC50 600 μM).

Design of HIV-1 integrase inhibitors targeting the catalytic domain as well as its interaction with LEDGF/p75: a scaffold hopping approach using salicylate and catechol groups

HIV-1 integrase (IN) is a validated therapeutic target for antiviral drug design. However, the emergence of viral strains resistant to clinically studied IN inhibitors demands the discovery of novel inhibitors that are structurally as well mechanistically different. Herein, we describe the design and discovery of novel IN inhibitors targeting the catalytic domain as well as its interaction with LEDGF/p75, which is essential for the HIV-1 integration as an IN cofactor. By merging the pharmacophores of salicylate and catechol, the 2,3-dihydroxybenzamide (5a) was identified as a new scaffold to inhibit the strand transfer reaction efficiently. Further structural modifications on the 2,3-dihydroxybenzamide scaffold revealed that the heteroaromatic functionality attached on the carboxamide portion and the piperidin-1-ylsulfonyl substituted at the phenyl ring are beneficial for the activity, resulting in a low micromolar IN inhibitor (5p, IC(50)=5 μM) with more than 40-fold selectivity for the strand transfer over the 3'-processing reaction. More significantly, this active scaffold remarkably inhibited the interaction between IN and LEDGF/p75 cofactor. The prototype example, N-(cyclohexylmethyl)-2,3-dihydroxy-5-(piperidin-1-ylsulfonyl) benzamide (5u) inhibited the IN-LEDGF/p75 interaction with an IC(50) value of 8 μM. Using molecular modeling, the mechanism of action was hypothesized to involve the chelation of the divalent metal ions inside the IN active site. Furthermore, the inhibitor of IN-LEDGF/p75 interaction was properly bound to the LEDGF/p75 binding site on IN. This work provides a new and efficient approach to evolve novel HIV-1 IN inhibitors from rational integration and optimization of previously reported inhibitors.

Synthesis, biological evaluation and 3D-QSAR studies of 3-keto salicylic acid chalcones and related amides as novel HIV-1 integrase inhibitors

HIV-1 integrase is one of the three most important enzymes required for viral replication and is therefore an attractive target for anti retroviral therapy. We herein report the design and synthesis of 3-keto salicylic acid chalcone derivatives as novel HIV-1 integrase inhibitors. The most active compound, 5-bromo-2-hydroxy-3-[3-(2,3,6-trichlorophenyl)acryloyl]benzoic acid (25) was selectively active against integrase strand transfer, with an IC(50) of 3.7 μM. While most of the compounds exhibited strand transfer selectivity, a few were nonselective, such as 5-bromo-3-[3-(4-bromophenyl)acryloyl]-2-hydroxybenzoic acid (15), which was active against both 3'-processing and strand transfer with IC(50) values of 11±4 and 5±2 μM, respectively. The compounds also inhibited HIV replication with potencies comparable with their integrase inhibitory potencies. Thus, 5-bromo-2-hydroxy-3-[3-(2,3,6-trichlorophenyl)acryloyl]benzoic acid (25) and 5-bromo-3-[3-(4-bromophenyl)acryloyl]-2-hydroxybenzoic acid (15) inhibited HIV-1 replication with EC(50) values of 7.3 and 8.7 μM, respectively. A PHASE pharmacophore hypothesis was developed and validated by 3D-QSAR, which gave a predictive r(2) of 0.57 for an external test set of ten compounds. Phamacophore derived molecular alignments were used for CoMFA and CoMSIA 3D-QSAR modeling. CoMSIA afforded the best model with q(2) and r(2) values of 0.54 and 0.94, respectively. This model predicted all the ten compounds of the test set within 0.56 log units of the actual pIC(50) values; and can be used to guide the rational design of more potent novel 3-keto salicylic acid integrase inhibitors.