HIVprotI: an integrated web based platform for prediction and design of HIV proteins inhibitors
Journal of Cheminformatics volume 10, Article number: 12 (2018)
A number of anti-retroviral drugs are being used for treating Human Immunodeficiency Virus (HIV) infection. Due to emergence of drug resistant strains, there is a constant quest to discover more effective anti-HIV compounds. In this endeavor, computational tools have proven useful in accelerating drug discovery. Although methods were published to design a class of compounds against a specific HIV protein, but an integrated web server for the same is lacking. Therefore, we have developed support vector machine based regression models using experimentally validated data from ChEMBL repository. Quantitative structure activity relationship based features were selected for predicting inhibition activity of a compound against HIV proteins namely protease (PR), reverse transcriptase (RT) and integrase (IN). The models presented a maximum Pearson correlation coefficient of 0.78, 0.76, 0.74 and 0.76, 0.68, 0.72 during tenfold cross-validation on IC50 and percent inhibition datasets of PR, RT, IN respectively. These models performed equally well on the independent datasets. Chemical space mapping, applicability domain analyses and other statistical tests further support robustness of the predictive models. Currently, we have identified a number of chemical descriptors that are imperative in predicting the compound inhibition potential. HIVprotI platform (http://bioinfo.imtech.res.in/manojk/hivproti) would be useful in virtual screening of inhibitors as well as designing of new molecules against the important HIV proteins for therapeutics development.
Human Immunodeficiency Virus (HIV) is one of the reasons for human death and suffering worldwide. It causes Acquired Immunodeficiency Syndrome (AIDS) in which gradual breakdown of the immune system allows critical opportunistic diseases to flourish . As per the UNAIDS report, around 78 million people have become infected with HIV and 35 million people have died of AIDS-related illnesses since the start of the epidemic. In 2015 alone there were about 36.9 million people living with HIV of which 1.1 million died (http://www.unaids.org/en/resources/campaigns/HowAIDSchangedeverything/factsheet). Due to the high genetic variability and mutation rate of HIV, vaccines are not available to curb the HIV infection .
Researchers have put a considerable focus on HIV therapy and a lot of compounds have been tested against this pathogen [3, 4]. However, a few antiretroviral drugs have been able to slow the disease progression. These drugs blocked the function of proteins implicated in certain stages of the HIV life-cycle . Different HIV enzymes are needed for the development of the retrovirus including reverse transcriptase (RT), protease (PR) and integrase (IN) . RT creates complementary DNA from an RNA template which can integrate into the host genome and its inhibitors are widely used as antiretroviral drugs . For example, the first anti-HIV drug zidovudine or azidothymidine (a nucleoside analog) was approved by the Food and Drug Administration (FDA) in 1987. It inhibits HIV reverse transcriptase, hence thwarting viral replication . PR slices the newly synthesized polyproteins at the relevant positions to form the mature protein apparatus and is a major drug-target for treatment of HIV . In 1995, saquinavir (invirase) became the first approved protease inhibitor. It blocks the enzyme’s active site, thus restricting the processing of HIV poly-proteins [10, 11]. The IN enzyme enables the virus to integrate its genetic material into the DNA of the host cell for a long-term infection. Compounds that inhibit the IN enzyme have demonstrated potent anti-HIV activity . For example, raltegravir (isentress), the first integrase inhibitor was approved by FDA in 2007 . Presently about 30 antiretroviral drugs are prescribed for the clinical treatment of AIDS . An improved knowledge of the structure and function viral proteins has led antiviral drug developers to design better antivirals to treat HIV infections .
To conserve capital and time for finding novel drugs, scientists have extensively used different computational approaches to scan virtual compound libraries prior to the wet lab experiments . The preferred targeted region should be off-target free and conserved across many strains of a virus for broad activity. Once the target is chosen, candidate antivirals can be selected by predicting the potential inhibitor using bioinformatics approaches [17, 18]. Amongst the diverse methods, quantitative structure activity relationship (QSAR) is being regularly used [19,20,21,22]. In QSAR, associations involving chemical descriptors and activity are employed to envisage the properties of other compounds . The chemical descriptors present the structural information of a compound as numerical values . Virtual screening employing QSAR is a valuable bioinformatics approach which helps to identify and devise of new antiviral drugs .
Several attempts have been made for predicting specific types of compounds against different HIV proteins (discussed later). Nevertheless, till date there no web server/software, which can regressively estimate the IC50/percentage inhibition activity of diverse types of, compounds against different HIV proteins. To accommodate this requirement, we created HIVprotI, a web based algorithm for prediction and design of protein specific anti-HIV compounds. In this approach, we employed experimentally validated inhibitors against RT, PR, IN (with IC50/percentage inhibition) from ChEMBL . We calculated molecular descriptors and performed feature selection to pick the best performing descriptors, which were employed to build support vector machine (SVM) based QSAR models for the prediction of inhibitors against HIV proteins. We further incorporated the models in the HIVprotI web server, which will be useful for virtual screening and scheming novel inhibitors directed against HIV.
In the present study, we have employed diverse datasets of inhibitors with experimentally validated IC50/percent inhibition activity against PR, RT and IN. The data was collected from ChEMBL resource (https://www.ebi.ac.uk/chembl/) by target browser (taxonomy tree) as well as target search using keywords like ‘Human Immunodeficiency Virus’, ‘HIV’, ‘protease’, ‘reverse transcriptase’ and ‘integrase’ etc. Initially among the inhibitors, majority of data belonged to RT, PR and IN with 3882, 3180, 2732 (IC50) and 1000, 740, 406 (percent inhibition) compounds respectively. After filtering entries with required information and eliminating redundant entries, we were left with 2126, 1895, 1240 (IC50) and 563, 518, 186 (percent inhibition) molecules correspondingly for the above mentioned proteins (hence 06 datasets) (Tables 1, 2). We have three times randomly picked ~ 10% of data as independent/validation dataset from each of six datasets. In each case this ~ 10% of the compounds were used for validation of the QSAR predictive models developed using the remaining 90% data during the training/testing . This process is iterated three times for each of the six datasets and performances were comparable as detailed in the Additional file 1: Table S1. These datasets were employed for descriptor calculation and development of the models. The datasets can be accessed from the URLs: http://bioinfo.imtech.res.in/manojk/hivproti/ic50_datasets.php and http://bioinfo.imtech.res.in/manojk/hivproti/datasets.php.
To develop protein specific prediction models, we calculated around 18,000 molecular descriptors which include geometric, electrostatic, structural, constitutional, path and graph fingerprints etc. utilizing the open source PaDEL software .
To increase the speed of computation and eliminate unrelated features we chose the most necessary molecular descriptors employing the filter ‘RemoveUseless’ and attribute evaluator ‘ClassifierSubsetEval’ with ‘BestFirst’ as the search method in Waikato Environment for Knowledge Analysis (Weka) suite . ClassifierSubsetEval estimates feature subsets on training/testing data utilizing a classifier to evaluate the worth of different feature sets.
We created protein specific QSAR models for each of the 3-inhibitor classes (RT, PR and IN) employing SMOreg machine learning algorithm in Weka package. SMOreg executes the SVM in regression approach . Chosen chemical descriptors and fingerprints were utilized for building the prediction models. The models were assessed by means of tenfold cross validation as well as independent validation . The overall methodology for model development is depicted in Fig. 1. However, we used Tropsha’s validation test to statistically validate the prediction ability of developed model [31, 32].
To calculate the performance of the QSAR models, we used various statistical measures such as Pearson’s correlation coefficient (PCC), Coefficient of Determination, Mean absolute error and Root-mean-square error described as follows.
Pearson’s correlation coefficient (R) calculates the correlation between two variables.
Here n is the size of test set while Eipred and Eiact are the predicted and actual efficacies correspondingly.
A PCC value of 1 implies full positive correlation, 0 implies no correlation while − 1 implies full negative correlation.
Coefficient of Determination (R2) signifies how well a data fits the statistical model. An R2 value of 1 states that the model totally fits the data. On the other hand, a value of 0 implies that the model does not fit the data in any way.
Mean absolute error (MAE) calculates the closeness of predictions to the actual results.
Here Eipred is the predicted value, Eiact the true and |Eipred–Eiact| the absolute error.
Root-mean-square error (RMSE) calculates the mean magnitude of the error.
MAEs and RMSEs are negatively-oriented values, that is, lower values are superior.
Chemical space mapping
Chemical space mapping of the datasets were done employing Ches-Mapper, a Java based application . It involves creation of 3-dimensional structure (mmff94 force field), feature extraction (Chemical Development Kit (CDK) and hashed finger prints), clustering (k-means cascade method), embedding in 3-D space (or Dimensionality reduction) and alignment of compounds using maximum common subgraph (MCS). IC50 and percentage inhibition datasets of RT, PR and IN were individually mapped in the chemical space to understand the relationship between structure, physicochemical properties and biological aspects. All the clusters were provided as superimposed images with information of the number of sequences, 3-D embedding quality, and respective MCS.
Applicability domain (AD) of a QSAR model helps to measure its certainty in prediction . We used Model Disturbance Index (MDI) v/s Prediction error (PE) method to calculate AD of all the prediction models . It is calculated through Java-based Applicability Domain-Model Disturbance Index (AD-MDI) software (http://nanobridges.eu/software/). Validation datasets of RT, PR and IN for both IC50 and percentage inhibition were used to check their reliability on the respective training data sets. AD of the models was provided in form of scatter plots between MDI and PE.
Performance of QSAR models
To facilitate the identification of the most efficient descriptors of anti HIV drugs against the three proteins, we calculated the correlation between chemical descriptors of the anti-HIV compounds and their IC50/percent inhibition. We have used wide-ranging datasets from the ChEMBL bioactivity resource . After feature selection, the relevant descriptors that remained were 45, 61, 55 (IC50) and 42 49, 23 (percent inhibition) for PR, RT and IN respectively.
During tenfold cross validation, we attained maximum correlation (PCC) of 0.78, 0.76, 0.74 (IC50) and 0.76, 0.68, 0.72 (percent inhibition) in case of PR, RT and IN respectively. Further we reached a maximum PCC of 0.73, 0.72, 0.70 (IC50) and 0.70, 0.63, 0.65 (percent inhibition) on independent validation datasets for PR, RT, and IN correspondingly (Tables 3, 4).
Furthermore, all the models are statistically significant with p value < 0.001. Other statistical parameters used in the creation of prediction models are listed in Additional file 1: Tables S2 and S3. Scatter plots of predicted and actual activities are depicted in Figs. 2 and 3.
A compilation of chosen molecular descriptors such as atom type electrotopological state, partial charge, extended topochemical atom and several path/graph fingerprints were helpful in developing the algorithm (Additional file 1: Tables S4 and S5). The atom type electrotopological state holds information of electronic state of the bonded atom in a compound and its topological nature in the milieu of the entire molecular structure . Likewise, the extended topochemical atom indices provide details of the electronic environment of the atoms, bonds, functional groups and branching  The details of other molecular descriptors have been discussed by Yap .
The prediction models have been incorporated into an open source and simple to web application, ‘HIVprotI’. Here one can predict the inhibition activity of query compounds against the different HIV proteins in terms of IC50/percent inhibition. The web server components include:
Using this module one can submit/sketch the query molecule [38, 39]. Users can select the proteins on which they desire to virtually screen the query compound. Following the submission of an input molecule by the user, its structure is optimized using Obminimize programme (https://openbabel.org/wiki/Obminimize) to optimize the geometry and minimize the energy for a molecule before descriptor calculation and prediction . On submission, it predicts IC50/percent inhibition activity against the HIV proteins. Users can also analyze the various properties of the query compound (Fig. 4).
Batch mode submission
Users can also submit multiple molecules simultaneously to check their inhibition efficiency against the desired HIV proteins. Clickable example molecules are given on the web server to help the users for easily getting started. This module will facilitate the researchers to virtually screen large number of compounds and select the ones with desired efficacy value. In addition, this component also enables the users to choose drug-like compounds by calculating the requisite properties. The batch mode can be accessed through the url: http://bioinfo.imtech.res.in/manojk/hivproti/batch.php.
Using this module, one can create analogs of their query structures based on user-defined components to evaluate the efficacy of the modified compounds on the selected HIV proteins. The structures are generated using SmiLib -a Java-based tool for rapid combinatorial library enumeration .
The result output displays predicted IC50/percent inhibition activity against the chosen HIV proteins. In addition one can view the various chemical attributes of the query compound like structure, Hydrogen/Lipinski bond donors/acceptors, rotatable/rigid bonds, logP value etc. to recognize drug-like compounds (Fig. 4).
HIVprotI also gives the users a search tool to find the data used in our study. In this component, the compounds targeting the PR, RT and IN proteins along with their structures are available in a database. The inhibitor entries can be easily searched and filtered from the web-application.
Comparison with existing approaches
There are several QSAR methods exist for predicting various categories of HIV protein inhibitors, which are compared with HIVprotI algorithm as shown in Table 5. Nonetheless, these approaches are very specific and deal with a selected group of inhibitors such as quinolines , pyrimidones , processing inhibitors  etc. Owing to this rationale they envisage the compounds which are similar to the inhibitor class with a good correlation, but do not perform well on other structurally dissimilar inhibitors for the same HIV protein. Majority of such studies are based on a restricted quantity of inhibitors. Moreover, none of them till date have provided any web server/software to enable the researchers to screen AVCs or compare the output from different studies. To check the performance of existing methods we used their datasets and developed QSAR models using HIVprotI approach. Performance of such developed models during training/testing was similar to the reported one, but they did not perform well on our independent validation datasets. It can be due to the low number of compounds used in those studies and it is seen that the performance increased with the rise in the number and types of compounds (Table 6). Moreover, there are certain online resources like PASS Online (www.way2drug.com/passonline), SuperPred  and ChemProt server  for deciphering the biological potential and targets of different chemical compounds using information of existing drugs, environmental chemicals and natural products databases. In addition, AVCpred algorithm developed by our group helps in the prediction of generic antiviral compounds . However, none of them is an HIV protein specific predictor. To address these limitations, HIVprotI has been developed employing more and varied inhibitors besides providing an open source integrated website for prediction and screening of protein specific anti-HIV compounds.
We used multidimensional scaling (MDS) to visualize the similitude of compounds in the different datasets. MDS calculates distance matrices by all-against-all evaluation of molecules from atom pair similarity values. The location or coordinates for each compound correspond to the distances graphically in a scatter plots. The plots were generated using Chemmine clustering workbench available at http://chemmine.ucr.edu/ . Compounds that are similar are close to one another while dissimilar compounds are placed far apart. The clustering analysis revealed that the datasets employed in the development of IC50 based QSAR models are chemically more diverse compared to the percent inhibition ones (Fig. 5). Further the compounds targeting RT and IN were comparatively more dispersed in the chemical space than those directed against PR. Since the QSAR models developed in this study are more in number as well as type in comparison to other studies (Additional file 1: Figure S1), hence the algorithm will be better in predicting diverse types of HIV protein inhibitors (https://figshare.com/articles/Additional_file_2/5607103).
Chemical space mapping
Chemical space mapping uses 3-D embedding and clustering to display the similarity among compounds i.e. spatial proximity between two compounds based on their feature similarity. Chemical space mapping was done for each dataset of IC50 and percentage inhibition. Clustering was performed through k-means cascade algorithm. The datasets were grouped into subgroups, with similar feature values in each cluster. For example, in case of IN inhibitors (IC50) with 1238 compounds, there were 10 clusters with 45, 113, 329, 25, 323, 185, 85, 12, 71 and 50 molecules respectively. Figures integrating mapped 3-D superimposed clusters, outline of individual cluster, MCS of each cluster along with information of the number of molecules in each cluster is shown in Fig. 6 (IC50 datasets) and Additional file 1: Figure S2 (percentage inhibition datasets).
To hinder HIV proliferation, the anti-HIV compounds aim at important proteins of HIV that are involved in various steps of its life cycle such as replication, transcription, maturation, integration, etc. . The drugs ought to be reasonably non-toxic to humans [49, 50]. The multiple stages of HIV life cycle can be inhibited employing compounds that can restrain viral enzymes like PR, RT, IN etc. that are requisite for HIV survival in the host cell. A growing list of these inhibitors are in clinical use and novel ones are also under trials .
Discovering innovative and enhanced drugs is a key objective in the management of HIV. Nevertheless, finding new inhibitors or compounds is a tedious procedure . To accelerate the development of new inhibitors, computational methods employing QSAR approaches are widely employed to optimize the research budget prior to experimentation . QSAR based algorithms have been extensively utilized in the selection of lead compounds and designing novel drugs .
In the present study, we created protein specific QSAR models to spot the probability of a given compound being an HIV protein inhibitor utilizing selected molecular descriptors of experimentally proven inhibitors against the HIV proteins. The open source PaDEL software was employed to compute numerous types of chemical descriptors followed by attribute selection approach to eliminate the irrelevant descriptors. We used machine learning to build the QSAR based models with good performance on various data sets of experimentally proven data from ChEMBL resource for specific HIV proteins. The developed models also displayed high performance while validated through independent data sets. Further, good predicting ability of the produced models was also observed by applying the statistical tests (Tropsha’s validation tests) for the continuous predictive models calculating Rext, k, k′ parameters as reported by (Additional file 1: Table S6) Golbraikh and Tropsha  and Vrontaki et al. . Simultaneously, the robustness of the QSAR models was also examined using Y-randomization test by comparing their performance (Q2 and R2) to the models generated using randomized inhibition values. We noticed that models with high Q2 and R2 values developed with actual inhibition values were consistent whereas there was no model with high Q2 and R2 values developed with randomized inhibition values (Additional file 1: Table S7) [52, 53]. However, the algorithm can be enhanced with the availability of more high-throughput data on these enzymes in future.
Chemical clustering for each dataset were performed through multidimensional scaling (MDS) as well as k-means cascade algorithm. The clusters comprised of aromatic/ringed compounds like toluene, ethenyl benzene etc. Although each class of inhibitors has characteristic compounds but we found some common MCS between different data sets e.g. 1,2,4-pentatriene between IC50 datasets of IN and RT, benzene derivatives (e.g. ethenyl benzene) in IN and PR inhibitor datasets etc. This indicates that some inhibitors can have multiple targets. Hence, chemical space mapping would help to fetch the information about characteristic as well as common compounds among each class of inhibitors and further assist in finding broad-spectrum anti-HIV drugs.
The applicability domain of the prediction models was verified by means of Williams plot (Figs. 7, 8) wherein standardized residuals are graphed against leverages . If the standardized residual of a molecule is more than thrice the standard deviation (± 3σ), the molecule is considered to be an outlier. The caution value of leverage (h*) is computed as 3p/n, where p represents the number of selected descriptors plus one and n represents the number of compounds used in training [55, 56]. If the leverage of a compound is more than h*, it is labeled as an outlier. The plots reveal that the leverages of bulk of the molecules do not go beyond the caution value (h*) in the QSAR models and thus the applicability of the models is reasonable. AD was also evaluated on independent test data sets by model population approach i.e. by checking relationship between MDI and PE (Additional file 1: Figures S3 and S4). The scatter plots further confirm the reliability of the QSAR models with low outliers.
We also compared our prediction method with earlier algorithms and found that the latter are optimized to a specific class of inhibitors and do not perform well when tested with other types of inhibitors for the same target protein. Besides none of the earlier published methods provide a software or web server for the researchers. On the other hand, HIVprotI web server has useful services like sketching new compounds and estimate their inhibition activity against multiple HIV proteins. Users can also screen several molecules concurrently using batch mode module on the web server. In addition searchable databases of both experimental and predicted datasets are also provided. The HIVprotI algorithm will aid the researchers in envisaging new anti-HIV compounds and virtually analyze the outcome of alterations on current drugs.
The HIVprotI is the first integrated web algorithm to predict anti-HIV compounds using experimentally verified data sets. Three QSAR prediction models for PR, RT and IN were developed to make all-inclusive predictions as well as screen compounds for their inhibition potential in a high throughput manner. The HIVprotI would be useful for scientists working in the field of anti-HIV therapeutics.
Acquired Immunodeficiency Syndrome
Food and Drug Administration
Human Immunodeficiency Virus
mean absolute error
quantitative structure-activity relationship
simplified molecular-input line-entry system
support vector machine
Maartens G, Celum C, Lewin SR (2014) HIV infection: epidemiology, pathogenesis, treatment, and prevention. Lancet 384(9939):258–271
Arts EJ, Hazuda DJ (2012) HIV-1 antiretroviral drug therapy. Cold Spring Harb Perspect Med 2(4):a007161
Hung TC, Lee WY, Chen KB, Chan YC, Chen CY (2014) Lead screening for HIV-1 integrase (IN) inhibited by traditional Chinese medicine. Biomed Res Int 2014:479367
Darnag R, Schmitzer A, Belmiloud Y, Villemin D, Jarid A, Chait A et al (2010) Quantitative structure-activity relationship studies of TIBO derivatives using support vector machines. SAR QSAR Environ Res 21(3–4):231–246
Flexner C (2007) HIV drug development: the next 25 years. Nat Rev Drug Discov 6(12):959–966
Barre-Sinoussi F, Ross AL, Delfraissy JF (2013) Past, present and future: 30 years of HIV research. Nat Rev Microbiol 11(12):877–883
Sarafianos SG, Marchand B, Das K, Himmel DM, Parniak MA, Hughes SH et al (2009) Structure and function of HIV-1 reverse transcriptase: molecular mechanisms of polymerization and inhibition. J Mol Biol 385(3):693–713
Brook I (1987) Approval of zidovudine (AZT) for acquired immunodeficiency syndrome. A challenge to the medical and pharmaceutical communities. JAMA 258(11):1517
Abbenante G, Fairlie DP (2005) Protease inhibitors in the clinic. Med Chem 1(1):71–104
James JS (1995) Saquinavir (Invirase): first protease inhibitor approved–reimbursement, information hotline numbers. AIDS Treat News 237:1–2
Bragman K (1996) Saquinavir: an HIV proteinase inhibitor. Adv Exp Med Biol 394:305–317
Prada N, Markowitz M (2010) Novel integrase inhibitors for HIV. Expert Opin Investig Drugs 19(9):1087–1098
Gandhi RT (2007) Raltegravir, a new HIV integrase inhibitor. AIDS Clin Care 19(12):105
Gu WG, Zhang X, Yuan JF (2014) Anti-HIV drug development through computational methods. AAPS J 16(4):674–680
De Clercq E (2009) Anti-HIV drugs: 25 compounds approved within 25 years after the discovery of HIV. Int J Antimicrob Agents 33(4):307–320
Plewczynski D, von Grotthuss M, Spieser SA, Rychlewski L, Wyrwicz LS, Ginalski K et al (2007) Target specific compound identification using a support vector machine. Comb Chem High Throughput Screen 10(3):189–196
De Clercq E (2002) Strategies in the design of antiviral drugs. Nat Rev Drug Discov 1(1):13–25
De Clercq E (2005) Antiviral drug discovery and development: where chemistry meets with biomedicine. Antiviral Res 67(2):56–75
Pissurlenkar RR, Khedkar VM, Iyer RP, Coutinho EC (2011) Ensemble QSAR: a QSAR method based on conformational ensembles and metric descriptors. J Comput Chem 32(10):2204–2218
Demchuk E, Ruiz P, Chou S, Fowler BA (2011) SAR/QSAR methods in public health practice. Toxicol Appl Pharmacol 254(2):192–197
Ruusmann V, Sild S, Maran U (2014) QSAR DataBank—an approach for the digital organization and archiving of QSAR model information. J Cheminform 6:25
Pourbasheer E, Riahi S, Ganjali MR, Norouzi P (2011) QSAR study of C allosteric binding site of HCV NS5B polymerase inhibitors by support vector machine. Mol Divers 15(3):645–653
Gonzalez-Diaz H, Romaris F, Duardo-Sanchez A, Perez-Montoto LG, Prado-Prado F, Patlewicz G et al (2010) Predicting drugs and proteins in parasite infections with topological indices of complex networks: theoretical backgrounds, applications, and legal issues. Curr Pharm Des 16(24):2737–2764
Cruz-Monteagudo M, Cordeiro MN, Tejera E, Dominguez ER, Borges F (2012) Desirability-based multi-objective QSAR in drug discovery. Mini Rev Med Chem 12(10):920–935
Prado-Prado FJ, Borges F, Uriarte E, Perez-Montoto LG, Gonzalez-Diaz H (2009) Multi-target spectral moment: QSAR for antiviral drugs vs. different viral species. Anal Chim Acta 651(2):159–164
Bento AP, Gaulton A, Hersey A, Bellis LJ, Chambers J, Davies M et al (2014) The ChEMBL bioactivity database: an update. Nucleic Acids Res. 42:D1083–D1090
Rajput A, Gupta AK, Kumar M (2015) Prediction and analysis of quorum sensing peptides based on sequence features. PLoS ONE 10(3):e0120066
Yap CW (2011) PaDEL-descriptor: an open source software to calculate molecular descriptors and fingerprints. J Comput Chem 32(7):1466–1474
Frank E, Hall M, Trigg L, Holmes G, Witten IH (2004) Data mining in bioinformatics using Weka. Bioinformatics 20(15):2479–2481
Qureshi A, Thakur N, Kumar M (2013) VIRsiRNApred: a web server for predicting inhibition efficacy of siRNAs targeting human viruses. J Transl Med 11:305
Golbraikh A, Tropsha A (2002) Beware of q2! J Mol Graph Model 20(4):269–276
Vrontaki E, Melagraki G, Mavromoustakos T, Afantitis A (2015) Exploiting ChEMBL database to identify indole analogs as HCV replication inhibitors. Methods 71:4–13
Gutlein M, Karwath A, Kramer S (2012) CheS-Mapper-chemical space mapping and visualization in 3D. J Cheminform 4(1):7
Weaver S, Gleeson MP (2008) The importance of the domain of applicability in QSAR modeling. J Mol Graph Model 26(8):1315–1326
Yan J, Zhu WW, Kong B, Lu HB, Yun YH, Huang JH et al (2014) A combinational strategy of model disturbance and outlier comparison to define applicability domain in quantitative structural activity relationship. Mol Inform 33(8):503–513
Hall LH, Kier LB (1995) Electrotopological state indices for atom types: a novel combination of electronic, topological, and valence state information. J Chem Inf Comput Sci 35(6):1039–1045
Das RN, Roy K (2013) QSPR with extended topochemical atom (ETA) indices. 4. Modeling aqueous solubility of drug like molecules and agrochemicals following OECD guidelines. Struct Chem 24(1):303–331
Melagraki G, Ntougkos E, Rinotas V, Papaneophytou C, Leonis G, Mavromoustakos T et al (2017) Cheminformatics-aided discovery of small-molecule Protein-Protein Interaction (PPI) dual inhibitors of Tumor Necrosis Factor (TNF) and Receptor Activator of NF-kappaB Ligand (RANKL). PLoS Comput Biol 13(4):e1005372
Braga RC, Alves VM, Muratov EN, Strickland J, Kleinstreuer N, Trospsha A et al (2017) Pred-skin: a fast and reliable web application to assess skin sensitization effect of chemicals. J Chem Inf Model 57(5):1013–1017
O’Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR (2011) Open Babel: an open chemical toolbox. J Cheminform 3:33
Schüller A, Hähnke V, Schneider G (2007) SmiLib v2.0: a Java-based tool for rapid combinatorial library enumeration. QSAR Comb Sci 26(3):407–410
Sun XH, Guan JQ, Tan JJ, Liu C, Wang CX (2012) 3D-QSAR studies of quinoline ring derivatives as HIV-1 integrase inhibitors. SAR QSAR Environ Res 23(7–8):683–703
Kaushik S, Gupta SP, Sharma PK, Anwer Z (2011) A QSAR study on a series of N-methyl pyrimidones acting as HIV integrase inhibitors. Indian J Biochem Biophys 48(6):427–434
Kong Y, Xuan S, Yan A (2014) Computational models on quantitative prediction of bioactivity of HIV-1 integrase 3′ processing inhibitors. SAR QSAR Environ Res 25(9):729–746
Nickel J, Gohlke BO, Erehman J, Banerjee P, Rong WW, Goede A et al (2014) SuperPred: update on drug classification and target prediction. Nucleic Acids Res 42:W26–W31
Taboureau O, Nielsen SK, Audouze K, Weinhold N, Edsgard D, Roque FS et al (2011) ChemProt: a disease chemical biology database. Nucleic Acids Res 39:D367–D372
Qureshi A, Kaur G, Kumar M (2017) AVCpred: an integrated web server for prediction and design of antiviral compounds. Chem Biol Drug Des 89(1):74–83
Backman TW, Cao Y, Girke T (2011) ChemMine tools: an online service for analyzing and clustering small molecules. Nucleic Acids Res 39:W486–W491
Littler E, Oberg B (2005) Achievements and challenges in antiviral drug discovery. Antivir Chem Chemother 16(3):155–168
Reusser P (2000) Antiviral therapy: current options and challenges. Schweiz Med Wochenschr 130(4):101–112
De Clercq E (2012) Human viral diseases: what is next for antiviral drug discovery? Curr Opin Virol 2(5):572–579
Melagraki G, Afantitis A (2015) A risk assessment tool for the virtual screening of metal oxide nanoparticles through Enalos InSilicoNano platform. Curr Top Med Chem 15(18):1827–1836
Zhang S, Golbraikh A, Oloff S, Kohn H, Tropsha A (2006) A novel automated lazy learning QSAR (ALL-QSAR) approach: method development, applications, and virtual screening of chemical databases using validated ALL-QSAR models. J Chem Inf Model 46(5):1984–1995
Minovski N, Zuperl S, Drgan V, Novic M (2013) Assessment of applicability domain for multivariate counter-propagation artificial neural network predictive models by minimum euclidean distance space analysis: a case study. Anal Chim Acta 759:28–42
Fatemi MH, Heidari A, Gharaghani S (2015) QSAR prediction of HIV-1 protease inhibitory activities using docking derived molecular descriptors. J Theor Biol 369:13–22
Zhao Q, Yang K, Li W, Xing B (2014) Concentration-dependent polyparameter linear free energy relationships to predict organic compound sorption on carbon nanotubes. Sci Rep 4:3888
Deeb O, Goodarzi M (2010) Exploring QSARs for inhibitory activity of non-peptide HIV-1 protease inhibitors by GA-PLS and GA-SVM. Chem Biol Drug Des 75(5):506–514
Basak SC, Mills D, Garg R, Bhhatarai B (2010) A QSAR study of HIV protease inhibitors using theoretical descriptors. Curr Comput Aided Drug Des 6(4):269–282
Silva DG, Freitas MP (2012) MIA-QSAR modeling of the anti-HIV-1 protease activities and cytotoxicities of ritonavir analogs. Curr Drug Saf 7(4):269–281
Ferreira LG, Andricopulo AD (2015) Fragment-based QSAR and structural analysis of a series of hydroxyethylamine derivatives as HIV-1 protease inhibitors. Comb Chem High Throughput Screen 18(5):464–475
Hu R, Doucet JP, Delamar M, Zhang R (2009) QSAR models for 2-amino-6-arylsulfonylbenzonitriles and congeners HIV-1 reverse transcriptase inhibitors based on linear and nonlinear regression methods. Eur J Med Chem 44(5):2158–2171
Mandal AS, Roy K (2009) Predictive QSAR modeling of HIV reverse transcriptase inhibitor TIBO derivatives. Eur J Med Chem 44(4):1509–1524
Sabet R, Fassihi A, Moeinifard B (2009) QSAR study of PETT derivatives as potent HIV-1 reverse transcriptase inhibitors. J Mol Graph Model 28(2):146–155
Latha RS, Vijayaraj R, Singam ER, Chitra K, Subramanian V (2011) 3D-QSAR and docking studies on the HEPT derivatives of HIV-1 reverse transcriptase. Chem Biol Drug Des 78(3):418–426
Veerasamy R, Subramaniam DK, Chean OC, Ying NM (2012) Designing hypothesis of substituted benzoxazinones as HIV-1 reverse transcriptase inhibitors: QSAR approach. J Enzyme Inhib Med Chem 27(5):693–707
Jain Pancholi N, Gupta S, Sapre N, Sapre NS (2014) Design of novel leads: ligand based computational modeling studies on non-nucleoside reverse transcriptase inhibitors (NNRTIs) of HIV-1. Mol BioSyst 10(2):313–325
Cheng Z, Zhang Y, Fu W (2010) QSAR study of carboxylic acid derivatives as HIV-1 Integrase inhibitors. Eur J Med Chem 45(9):3970–3980
Gupta P, Sharma A, Garg P, Roy N (2013) QSAR study of curcumine derivatives as HIV-1 integrase inhibitors. Curr Comput Aided Drug Des 9(1):141–150
MK conceived, designed and provided valuable guidance throughout the implementation. AQ collected the data. AQ and GK developed the predictive models. AR and AQ performed the chemical space mapping and applicability domain analyses. AQ and AR designed the web interface. AQ, AR and MK wrote the manuscript. All authors read and approved the final manuscript.
The project was funded by Council of Scientific and Industrial Research (GENESIS-BSC0121) and Department of Biotechnology (GAP0001).
The authors declare that they have no competing interests.
Availability and requirements
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supporting information including Table S1. Performance of QSAR predictive models on three times randomly picked ~ 10% independent/validation data. These models were developed using remaining ~ 90% data during training/testing respectively for each of the six datasets; Table S2. Details of statistical parameters used for the development of IC50 based QSAR models; Table S3. Details of statistical parameters used for the development of percent inhibition based QSAR models; Table S4. Details of chemical descriptors used in the development of IC50 based QSAR models; Table S5. Details of chemical descriptors used in the development of percent inhibition based QSAR models; Table S6. Details of slopes k (predicted vs. observed inhibition) and k’ (observed vs. predicted inhibition) of the regression lines for the QSAR models; Table S7. Details of Y-randomization test performed on the QSAR models; Figure S1. Chemical space analysis of QSAR studies (Table 5) for Protease (PR) (a, b), Reverse Transcriptase (RT) (c, d) and Integrase (IN) (e, f) respectively; Figure S2. Chemical space mapping outline of (a) Integrase (IN), (b) Protease (PR) and (c) Reverse Transcriptase (RT) inhibitors (percentage inhibition) with internal circle showing clustering and 3-D embedding of compounds, middle circle with exact (zoomed) superimposed cluster and outermost circle with specific MCS of each cluster; Figure S3. Scatter plot depicting the applicability domain for IC50 datasets of (a) Integrase (IN), (b) Protease (PR) and (c) Reverse Transcriptase (RT); Figure S4. Scatter plot depicting the applicability domain for percentage inhibition datasets of (a) Integrase (IN), (b) Protease (PR) and (c) Reverse Transcriptase (RT).
Source code of HIVProtI web server.
About this article
Cite this article
Qureshi, A., Rajput, A., Kaur, G. et al. HIVprotI: an integrated web based platform for prediction and design of HIV proteins inhibitors. J Cheminform 10, 12 (2018). https://doi.org/10.1186/s13321-018-0266-y