Next Article in Journal
Clinical, Pathological and Virological Outcomes of Tissue-Homogenate-Derived and Cell-Adapted Strains of Porcine Epidemic Diarrhea Virus (PEDV) in a Neonatal Pig Model
Next Article in Special Issue
Targeting Yezo Virus Structural Proteins for Multi-Epitope Vaccine Design Using Immunoinformatics Approach
Previous Article in Journal
Epidemiology and Diagnostics of Cacao Swollen Shoot Disease in Ghana: Past Research Achievements and Knowledge Gaps to Guide Future Research
Previous Article in Special Issue
Computational Analysis of CD46 Protein Interaction with SARS-CoV-2 Structural Proteins: Elucidating a Putative Viral Entry Mechanism into Human Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anti-Dengue: A Machine Learning-Assisted Prediction of Small Molecule Antivirals against Dengue Virus and Implications in Drug Repurposing

1
Virology Unit, Institute of Microbial Technology, Council of Scientific and Industrial Research (CSIR), Sector 39A, Chandigarh 160036, India
2
Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
*
Author to whom correspondence should be addressed.
Viruses 2024, 16(1), 45; https://doi.org/10.3390/v16010045
Submission received: 5 October 2023 / Revised: 20 December 2023 / Accepted: 21 December 2023 / Published: 27 December 2023
(This article belongs to the Special Issue Computational Drug Discovery for Viral Infections)

Abstract

:
Dengue outbreaks persist in global tropical regions, lacking approved antivirals, necessitating critical therapeutic development against the virus. In this context, we developed the “Anti-Dengue” algorithm that predicts dengue virus inhibitors using a quantitative structure–activity relationship (QSAR) and MLTs. Using the “DrugRepV” database, we extracted chemicals (small molecules) and repurposed drugs targeting the dengue virus with their corresponding IC50 values. Then, molecular descriptors and fingerprints were computed for these molecules using PaDEL software. Further, these molecules were split into training/testing and independent validation datasets. We developed regression-based predictive models employing 10-fold cross-validation using a variety of machine learning approaches, including SVM, ANN, kNN, and RF. The best predictive model yielded a PCC of 0.71 on the training/testing dataset and 0.81 on the independent validation dataset. The created model’s reliability and robustness were assessed using William’s plot, scatter plot, decoy set, and chemical clustering analyses. Predictive models were utilized to identify possible drug candidates that could be repurposed. We identified goserelin, gonadorelin, and nafarelin as potential repurposed drugs with high pIC50 values. “Anti-Dengue” may be beneficial in accelerating antiviral drug development against the dengue virus.

Graphical Abstract

1. Introduction

Dengue, a viral disease transmitted by mosquitoes, exhibits a rapid transmission rate and is particularly common in tropical and subtropical areas. Consequently, it presents a substantial burden in terms of both mortality and morbidity [1]. Dengue was first registered in 1780 in Madras (now Chennai). The initial virology-confirmed outbreak occurred in Calcutta and along India’s eastern coast from 1963 to 1964. DHF was recorded in the Philippines in 1953–1954 [2]. Since 1950, frequent dengue outbreaks have occurred in Southeast Asian countries [3]. The World Health Organization (WHO) has reported a significant increase in the global burden of dengue over the past two decades. Roughly half of the world’s population is at risk of dengue infection, with an estimated 100 to 400 million infections yearly [4,5].
Dengue virus (DENV) is a single positive-stranded RNA virus belonging to the genus Flavivirus and Flaviviridae family. DENV has four serotypes, namely DENV-1, DENV-2, DENV-3, and DENV-4. DENV transmission occurs through the bite of DENV-carrying female Aedes mosquitoes, mainly by Aedes aegypti and, rarely, by Aedes albopictus, which leads to severe health issues known as dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) [5]. DENV comprises 10,723 nucleotides (approximately 11 kb), enciphering larger polyprotein precursors containing ~3391 amino acid residues. DENV polyproteins, after cleavage by host and virus proteases, constitute three structural proteins named C (capsid); prM (pre-membrane); E (envelope); and seven nonstructural proteins called NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5 [6]. The clinical symptoms of dengue infection can range from asymptomatic to severe illnesses that may result in fatalities. Different categories exist in symptomatic cases, including mild acute undifferentiated febrile illness (UF), DF, DHF, DSS, and uncommon dengue (UD) or expanded dengue syndrome (EDS) [7].
Several research groups identified novel potent DENV inhibitors. Low, June Su Yin et al. identified Narasin as a novel antiviral agent with an IC50 of less than 1 μM against all DENV serotypes [8]. In another study, Raekiansyah, Muhareva et al. highlighted Brefeldin A as a promising and novel antiviral compound, displaying an IC50 range of 54.6 to 65.7 nM against all DENV serotypes [9]. Bardiot et al. studied the potential of KU Leuven’s compound library in inhibiting DENV-2 through a CPE reduction assay. They determined a promising inhibitor, 2-((3,4-dimethoxyphenyl) amino)-1-(1H-indol-3-yl)-2-phenylethan-1-one, for DENV [10].
Further, several experimental studies have been performed to determine the activity of repurposed drugs against the DENV. Drug repurposing could be a promising approach to looking for effective antivirals against the DENV. For example, quinine [11], N-Acetylcysteine [12], and Antiemetic Metoclopramide [13] have been used as repurposed drugs against DENV. Likewise, many more antivirals as potential repurposed drug candidates have been explored against the DENV [14]. Still, fewer antivirals are under clinical trial; therefore, we must explore more chemicals/inhibitors to get a highly effective and potent antiviral against DENV.
In this endeavor, computational approaches can be used to predict potent antivirals to reduce the time and cost. It could also be advantageous to accelerate the drug discovery process. In light of this, our group developed various machine learning-based antiviral predictors using the quantitative structure–activity relationship (QSAR) information of molecules/peptides such as AVCpred [15], AVPpred [16], AVP-IC50 Pred [17], HIVprotI [18], Anti-flavi [19], etc. Recently, we have developed a predictive algorithm for SARS-CoV-2, i.e., anti-corona [20], and for Ebola virus, i.e., anti-Ebola [21]. However, the platform is required to predict the repurposed drugs targeting the DENV utilizing machine learning techniques (MLTs).
In this study, we developed the “Anti-Dengue” predictive algorithm using various MLTs like support vector machine (SVM), artificial neural network (ANN), k-nearest neighbor (kNN), and random forest (RF). This algorithm predicts the efficacy of chemicals and drugs against DENV by assessing their inhibition efficiency, measured in terms of pIC50 and IC50 values (μM). Furthermore, we have also identified various effective repurposed drug candidates by scanning the “DrugBank” database through the best predictive model.

2. Materials and Methods

For developing the anti-dengue predictor, the workflow is given in Figure 1.

2.1. Data Collection

The antiviral entries were procured from the “DrugRepV” database to develop an “Anti-Dengue” predictor. The “DrugRepV” database encompasses chemicals (small molecules) and repurposed drugs designed to target epidemic and pandemic viruses, comprising a total of 8485 entries. This dataset provides comprehensive information, including antiviral names, drug types, primary and secondary indications, viral strains, pathways, assay details, clinical status, and more [22].
The steps for fetching out the antiviral entries are given below using the standard method [23]:
  • We obtained 900 antiviral entries for the DENV in the “DrugRepV” database.
  • The antiviral entries were filtered based on IC50/EC50 values, SMILES, molecular weight, etc. to acquire only relevant candidates, i.e., 238.
  • Using the formula pIC50 = −log10(IC50(M)), the IC50 is converted into pIC50, where the IC50 is the dimensionless activity that can be expressed in molar concentrations. Higher values of pIC50 showed greater potency and vice versa.
The dataset containing drugs/inhibitors to create the model is given in Supplementary Table S1.

2.2. Descriptor Calculation

The chemical structures of antiviral candidates were used to procure the chemical information, such as the simplified molecular-input line-entry system (SMILES), then reformed into 3D-SDF format utilizing the open Babel v3.1.1 command line tool [24]. Further, these SDF files served as inputs for withdrawing chemical descriptors and fingerprints.

2.3. Compounds/Inhibitors Feature Extraction

The computation of 1D, 2D, and 3D molecular descriptors and fingerprints using 3D-SDF structures was performed using PaDEL software (version 2.21) to calculate 17,968 descriptors [25]. One-dimensional descriptors are substructural molecular fragment-based descriptors (H-Bond acceptor/donor, fingerprints, fragments count, etc.). Two-dimensional descriptors are structural and physicochemical properties-based descriptors (topological and electronic information, topological descriptors, connectivity indices, etc.). Three-dimensional descriptors are derived from the 3D conformation of the molecules (geometrical, as well as spatial, information of molecules, comparative molecular similarity index analysis (CoMSIA), solvent accessible area, comparative molecular field analysis (CoMFA), polar and nonpolar surface areas (PSAs and NPSAs), etc. [26]. Molecular fingerprints are another way of depicting the molecule structure where binary digits (bits) help find or differentiate between the specific substructures in the molecule. The descriptors and fingerprints are essential when studying drugs or chemicals to determine their QSAR [27].

2.4. Feature Selection

Feature selection involves identifying and eliminating redundant and irrelevant features to obtain significant features that can improve the accuracy of the developed models [28]. The feature selection was performed with the help of the perceptron, SVR, and DT methods in the recursive feature elimination (RFE) module available in the scikit-learn library to find the top 50, 100, 150, and 200 relevant features. Among these, the top 100 features of the perceptron method were used as input for implementing the machine learning algorithms in this study [29,30].

2.5. Machine Learning Algorithms

In this current study, we involved the implementation of SVM, ANN, kNN, and RF.
SVM is a supervised machine learning algorithm that can be utilized for regression and classification tasks. It generally creates several hyperplanes but needs to find the best hyperplane with a maximum margin that classifies the data more accurately. There are two categories of SVM, namely linear SVM and nonlinear SVM. Linear SVM is typically used for data that can be separated linearly, while nonlinear SVM is designed for data that cannot be separated linearly. The kernel function is also used to alter the training data with the help of which nonlinear decision surface is converted to a linear equation, i.e., usable form for data processing [19].
RF is also a supervised machine learning algorithm that can be utilized for regression and classification tasks. RF performs functions by forming decision trees using a training dataset, and the outturn it makes is based on the mean prediction [31].
An ANN is an effort to imitate the neuron network that comprises the human brain to make the computer learn things and respond accordingly as humans do. It typically comprises three layers: the input, hidden, and output layers. These layers transform the input into a meaningful output [32].
The kNN algorithm is a MLT that does not assume any specific form for the underlying data distribution and is supervised in nature. It can be applied to perform classification or regression tasks [33]. It is frequently known as memory-based, instance-based, or lazy learning. It is based on the pick out of the nearest neighbor for a query data point based on the distance, which can be calculated by Euclidean distance, Minkowski distance, Manhattan distance, Hamming distance, etc.

2.6. Generation of Random Datasets

To create independent validation datasets, we used a random selection process to choose approximately 10% of the available data, while the remaining 90% was utilized for training and testing purposes of the models. We repeated this procedure five times, resulting in five sets of training/testing and independent validation data, each containing 238 molecules (T214 + V24).

2.7. Ten-Fold Cross-Validation

To assess the performance of the machine learning predictive models, we employed the ten-fold cross-validation method. This technique involved splitting the training/testing dataset into ten equal parts. During each iteration, nine parts were combined for training, while the remaining part was used for testing to assess the model’s performance. All ten parts were used as testing data at least once, and the overall model performance was evaluated based on the average performance of all the testing parts. Additionally, to validate the performance of the developed model, we used an independent/external dataset that was not utilized during the model’s training and testing.

2.8. Model Performance Assessment

The developed model performance was evaluated by calculating the mean absolute error (MAE), mean squared error (MSE), root mean squared error (RMSE), coefficient of determination (R2), and Pearson’s correlation coefficient (PCC or R) using the formulas as given below.
P C C = n n = 1 n   E i a c t E i p r e d n = 1 n   E i a c t n = 1 n   E i p r e d n n = 1 n   ( E i a c t ) 2 ( n = 1 n   E i a c t ) 2 n n = 1 n   ( E i p r e d ) 2 ( n = 1 n   E i p r e d ) 2
M A E = 1 n n = 1 n   E i p r e d E i a c t
M S E = 1 n n = 1 n   E i p r e d E i a c t ²
R M S E = 1 n n = 1 n   ( E i p r e d E i a c t ) 2
where n, Eact, and Epred are the dataset size and actual and predicted values, respectively.

2.9. Applicability Domain Analysis

Moreover, along with the model performance, model accuracy for the new prediction also plays a crucial role. Applicability domain analysis defines the boundary of the developed model for its reliability. For accurate predictions of a new compound using a developed model, it is essential for the chemical properties of the compound to fall within the applicability domain of the compounds employed in training the model [34]. The reliability of these developed models was assessed using the William’s plot based on the distance-based leverage approach. These plots depict the relationship between the leverage and standardized residuals. The formula of the leverage threshold (h*) is
Leverage threshold (h*) = 3(p + 1)/n
where p = number of descriptors utilized in developing the model; n = number of compounds used in the training dataset.
The reliability of the predicted model was observed to be dependent on a majority of the data points falling within the leverage threshold (h*). To confirm the strength and effectiveness of the developed models created using the SVM, RF, kNN, and ANN algorithms, we plotted a scatter plot between the predicted pIC50 values and actual pIC50 values.

2.10. Decoy Sets Analysis

Decoys were generated for these drug candidates using the DecoyFinder 2.0 tool [35]. DecoyFinder 2.0 utilizes a molecular weight-based method to generate decoys. The ZINC20 database was used as a source of a subset containing 4.78 million drug-like molecules to make the decoys [36]. Six decoy datasets were developed, having 238 random decoys of active drug candidates. Further format conversion and molecular descriptors were calculated to determine the pIC50 values. Eventually, a correlation was made in terms of the PCC between the decoy pIC50 and actual pIC50 of each decoy dataset’s equivalent active drug candidates.

2.11. Chemical Clustering Analysis

The chemical diversity of these drug candidates was evaluated by executing chemical clustering using ChemMine tools. We used the multidimensional scaling (MDS) algorithm and Binning clustering with the same similarity cut-off of 0.6 in both methods [37].

2.12. Drug Repurposing

Using the best predictive model based on SVM, we predict the potent repurposed drug candidates by scanning the more than 2000 FDA-approved drugs present within the DrugBank database [38]. We excluded those drugs from our DrugBank scanning approach that were already used in the model development. We converted the file format of these drugs and generated 17,968 molecular descriptors using PaDEL software. Further, we extracted the top 100 perceptron features involved in developing the best model. Subsequently, these DrugBank drugs, along with the 100 features, were employed to predict novel, potentially effective repurposed drug candidates with elevated pIC50 values against DENV.

2.13. Web Server Development

The best-performing SVM predictive model was implemented on the “Anti-Dengue” web server to assess the effectiveness of chemicals and drugs in inhibiting the DENV, as indicated by inhibition efficiencies such as the pIC50 and IC50 values (μM). The “Anti-Dengue” web server was constructed utilizing LAMP software (Ubuntu 12.04.2 LTS), incorporating Linux as the operating system, Apache as the web server, MySQL as the relational database management system, and PHP (Perl or Python) as the object-oriented scripting language. The front end of the “Anti-Dengue” web server was developed using HTML, CSS, and PHP, while the scripting languages, viz., python, perl, and JavaScript, were used at the back end of the web server. The web server predicts the inhibition efficiency in terms of the IC50 and pIC50 on the best-performing SVM model. To enhance user accessibility, we provide dedicated web pages such as “Help” and “Frequently Asked Questions” on the server for user guidance and assistance.

3. Results

3.1. Feature Selection Approach

Among all 17,968 descriptors, the top 100 features of the drugs were selected for developing the models. In the case of the support vector regression (SVR) method, the features are E1i, geomShape, FP258, KRFP320, KRFP307, ExtFP465, KRFPC3056, etc. Similarly, in the decision tree (DT) regression method, the features are SubFPC26, AATSC3m, ATSC1i, ATSC8p, ATSC8e, ATSC6e, ATSC6m, etc. Moreover, the perceptron method’s components are KRFPC52, ExtFP897, E3u, E2m, FP258, ExtFP41, ExtFP953, etc. The complete list of the top 100 features that were extracted using these three methods (SVR, DT, and perceptron) of the recursive feature elimination module is provided in Supplementary Table S2.

3.2. Performance of Developed Machine Learning-Based QSAR Models

To identify inhibitors of the DENV, we developed robust prediction models using four MLTs. These methods included SVM, ANN, kNN, and RF. The predicted models were developed using 100 top features/descriptors selected using the RFE module from the scikit-learn library.
Various statistical measures were utilized to evaluate the effectiveness of the developed QSAR models, including the MAE, MSE, RMSE, R2, and PCC. The MAE, or mean absolute error, is a metric used to measure the average magnitude of errors between the predicted and actual values. It is calculated by taking the average of the absolute differences between each predicted value and its corresponding actual value. The MAE tells about the closeness of the predicted values to the actual values. These values are negative-oriented values; that is, the more negative values, the more superior the developed model.
The MSE, or mean squared error, is a metric commonly used to quantify the average squared difference between predicted and actual values. It involves calculating the squared differences for each data point, averaging these squared differences, and then taking the square root to obtain the final result. The MSE gives more weight to larger errors than smaller ones, making them sensitive to outliers.
The RMSE measures the average magnitude of the errors between the predicted and actual values, with the square root applied to make the result more interpretable in the same units as the original data.
An R2 value of 1 depicts the data perfectly fitting into the model, whereas a value of 0 shows that the data do not fit into the model at all.
PCC values show the correlation between the inhibitors’ predicted and actual pIC50 values. PCC values lie between −1 and +1, where the −1 value shows a negative correlation, 0 values depict no correlation, and the +1 value implies a positive correlation. The R2 values show how well the data can fit in a statistical model.
The training and testing datasets for the DENV prediction models exhibited PCC values of 0.71 for SVM, 0.65 for ANN, 0.34 for kNN, and 0.45 for RF. For an independent dataset, the PCC values were 0.81 for SVM, 0.74 for ANN, 0.68 for kNN, and 0.54 for RF. The performance metrics for the best models developed using SVM, RF, kNN, and ANN for the DENV are presented in Table 1, Table 2, Table 3 and Table 4. Further information about all of the models developed for DENV inhibitors can be found in Supplementary Table S3. Detailed information on the actual and predicted IC50 of the independent validation dataset is available in Supplementary Table S4.

3.3. Applicability Domain Analysis

An applicability domain analysis using a William’s plot showed the leverage threshold (h*) value comes out to be 1.415 for models predicted using algorithms. Out of four predictive algorithms, the SVM model was found to be reliable, as most of the data points lie within the leverage threshold (h*), as given in Figure 2. Figure 3 displays a scatter plot between the actual pIC50 values and predicted pIC50 values for both the training/testing and independent validation datasets, illustrating that most of the points are clustered around the trend line. This indicates that the developed QSAR models are highly reliable. Supplementary Table S5 contains the information used for the William’s plot in the applicability domain analysis. Supplementary Table S6 contains information about the actual and predicted pIC50 values for the scatter plot.

3.4. Validation Using the Decoy Set

Unlike active molecules, decoys refer to molecules that cannot bind to their target. To confirm the predictive model’s reliability, the inhibitory activity in terms of the pIC50 was calculated for all six random decoy sets and then compared in terms of pIC50 with their corresponding active molecules (Supplementary Table S7). Decoy sets 1–6 showed the PCC values 0.117, 0.045, −0.0002, −0.091, −0.043, and −0.028, respectively, and their correlation is displayed using a scatter plot in Figure 4.

3.5. Chemical Diversity Analysis

A chemical diversity analysis was conducted to check the structural heterogeneity of the anti-dengue chemical compounds. A binning clustering analysis revealed that anti-dengue chemical compounds could be sorted into 124 bins or clusters (Supplementary Table S8). A 2D and 3D multidimensional scaling plot was generated to illustrate the dissimilarity of anti-dengue chemical compounds in chemical space, utilizing the same similarity cut-off as the binning clustering analysis, as shown in Figure 5.

3.6. Prediction of Promising Repurposed Anti-Dengue Drug Candidates

The most effective predictive model, based on SVM, was utilized to forecast repurposed drugs from the approved drugs category of the “DrugBank” database. The top 25 predicted candidates are listed in Table 5.

3.7. Anti-Dengue Web Server

To predict the effectiveness of anti-dengue chemicals, users should paste/upload the input in SDF format. The output will be received in a tabular format that includes Query ID, SMILES, the inhibition efficiency as pIC50 and IC50 (μM), 2D structure, and descriptor. The computation time for unknown chemicals typically ranges between 2 and 5 min. Users can keep track of their jobs by noting the job ID and accessing the “check job status” page to retrieve the results at any time. The “Anti-Dengue” web server is freely available at https://bioinfo.imtech.res.in/manojk/antidengue/.

4. Discussion

Dengue is an emerging health problem across the globe. Due to the absence of approved antiviral treatments or a universal vaccine for DENV infection, several research teams are focused on developing inhibitors that target various components, such as structural, nonstructural, host, and non-specific targets. In this concern, focusing on computational approaches for developing antivirals would be a better step to accelerate drug discovery research [39]. Hence, in the present research work, we developed a machine learning-based prediction algorithm, “Anti-Dengue”, to identify new potential repurposed drug candidates targeting DENV.
In this study, we employed multiple machine learning techniques (MLTs): support vector machine (SVM), artificial neural network (ANN), k-nearest neighbor (kNN), and random forest (RF) to develop a better predictive model. Additionally, we explored three feature selection methods: perceptron, SVR, and DT. By combining these MLTs with four feature sets comprising the top 50, 100, 150, and 200 features and considering five random datasets (214 molecules in training/testing and 24 molecules in independent datasets generated from 238 molecules), we developed a total of 240 models. Following an assessment of the performance parameters, such as the mean absolute error (MAE), mean squared error (MSE), root mean squared error (RMSE), coefficient of determination (R2), and Pearson’s correlation coefficient (PCC or R), of these models, we provided 12 predictive models details in Table 1, Table 2, Table 3 and Table 4. Finally, we selected a specific model for further analyses like the applicability domain, scatter plot, decoy dataset, etc. This chosen model is characterized by 100 features utilizing the perceptron feature selection method. Detailed information on all MLTs with 100 feature sets using all three feature selection methods and random sets is provided in Supplementary Table S3. This SVM model was integrated into the web implementation and employed to predict potential repurposed drug candidates against the dengue virus, and the top 25 predicted drug candidates are listed in Table 5.
We utilized four different MLTs, namely SVM, RF, ANN, and kNN, to develop highly effective predictive models. These MLTs have been employed by various researchers in a multitude of studies [40]. For example, Mpropred for the prediction of SARS-CoV-2 main protease antagonists [41], TargIDe for predicting the molecules with antibiofilm activity against Pseudomonas aeruginosa [42], EBOLApred for predicting cell entry inhibitors against the Ebola virus [43], and StackHCV for the identification of inhibitors against the NS5 protein of the Hepatitis C virus [44]. Similarly, we have utilized these techniques to create predictive algorithms such as AVCpred for predicting general effective antiviral compounds [15]: AVPpred, the first algorithm for predicting antiviral peptides [16], AVP-IC50 Pred for predicting antiviral peptides activity in terms of the IC50, i.e., the half-maximal inhibitory concentration [17], HIVprotI for predicting and designing inhibitors targeting Human Immunodeficiency Virus (HIV) proteins [18], and anti-flavi for predicting and designing various novel antiviral compounds, particularly for flaviviruses [19]. Recently, some predictive algorithms were developed for predicting repurposed drugs/inhibitors specifically for a virus, such as anti-Ebola for the Ebola virus [21] and anti-corona for SARS-CoV-2 [20]. To develop the predictor in the context of the DENV, we extracted the most relevant features from the 17,968 molecular descriptors and fingerprints. Out of all the MLTs employed to construct the predictive models, SVM outperformed RF, kNN, and ANN. SVM produced a PCC of 0.71 on the training/testing dataset and 0.81 on the independent validation dataset.
Further model robustness was cross-checked by plotting a William’s plot in the applicability domain analysis and plotting the actual vs. predicted pIC50 values to validate the robustness of the predicted model. We used the decoys of each active drug candidate to further check the reliability of the “Anti-Dengue” predictive models. Then, we compared the pIC50 values of inactive decoy molecules with their corresponding active molecule, which further confirms the reliability and robustness of the developed “Anti-Dengue” predictive models.
Furthermore, a chemical clustering analysis for the 238 molecules was also assessed using the multidimensional scaling (MDS) algorithm and binning clustering methods. Chemical clustering is generally used to identify outliers and understand chemical compounds’ arrangement in a chemical space. The binning clustering method made the chemical clusters based on the user-defined similarity cut-off values. We used a Tanimoto coefficient (Tc) of 0.6 as the similarity coefficient, which is the proportion of the features shared among two compounds divided by their union, i.e., c/(a + b + c), where c is the number of features common in both compounds, while a and b are the number of features that are unique in one or the other compound, respectively [45]. The Tanimoto coefficient value generally lies between 0 and 1, with higher values depicting greater similarity and vice versa. Using a Tc of 0.6 showed that compounds are joining with 0.6 similarity or more to aggregate numerous clusters using the “single linkage” rule. As many clusters are forming in the anti-dengue chemicals, they are well spread in the chemical space. The binning cluster results are represented in tabular form with the compound ID, bin/cluster size, and bin/cluster ID. Multidimensional scaling (MDS) creates a matrix of “item-to-item” distances, and each item is assigned with coordinates and represents these distances in the form of 2D and 3D scatter plots. MDS-generated plots show that anti-dengue chemicals are well distributed in the 2D and 3D chemical space. Binning clustering utilizes internally developed C++ implementation, and MDS uses the “cmdscale” function implemented in R. These methods showed that these chemicals are very dissimilar [20,46].
The developed predictive model identified several potentially effective repurposed drugs for the treatment of DENV from the “approved” drugs category within the DrugBank database. Furthermore, we conducted a literature review to verify the status of the top predicted drugs. We discovered that some hits have been investigated through experimental reports or in silico analysis. For example, Carro, Ana C., Luana E. Piccini, and Elsa B. Damonte tested chlorpromazine as an endocytic inhibitor against DENV-2 entry into myeloid cells in the presence or absence of antibodies [47]. Similarly, Shahen, Mohamed et al. showed that Loratadine (LRD), along with ReDuNing (RDN) and Acetaminophen, decreases the susceptibility, as well as the severity of, DENV by targeting the miRNA interacting with the potential target genes [48]. Likewise, Boonyasuppayakorn, Siwaporn et al. checked Primaquine, along with known FDA-approved antimalarial drugs like chloroquine and amodiaquine, to inhibit the viral proteases and DENV replication using protease, as well as reporter replication-based assays [49]. Malakar, Shilu et al. evaluated the four Food and Drug Administration (FDA)-approved drugs: azelaic acid, quinine sulfate, aminolevullic acid, and mitoxantrone hydrochloride. Quinine had the most potent activity against the DENV-2 virus strain. Quinine was found to inhibit DENV production by 80% compared to the controls. In a dose-dependent manner, it decreased DENV RNA and viral protein synthesis, consequently impeding replication [11]. Therefore, repurposed drug candidates predicted from our method have the potential to work as antiviral agents that could accelerate the drug discovery process for combating DENV infection.
Several researchers have conducted in silico studies aimed at identifying repurposed drugs against the DENV. These studies encompassed techniques like the transcriptomics-based bioinformatics approach, molecular simulations, molecular docking, pharmacophore model-based drug repurposing, and others [50,51]. These studies include datasets like phytocompound databases, natural products, small molecules, and FDA-approved drugs. Nonetheless, our study diverged from these methodologies, as we integrated four distinct MLTs to predict agents with anti-dengue properties. To develop the predictive models, we employed a range of chemically diverse anti-dengue compounds that have been experimentally validated by different research groups. Additionally, our best predictive models have been integrated into the web server, a feature that sets them apart from any previously documented computational studies for the DENV.
Recurring occurrences of DENV outbreaks characterized by significant mortality and fatality rates are causing significant global apprehension, as there is no approved drug or universal vaccine available for the treatment of DENV infection. Therefore, utilizing computational methods could prove highly beneficial in accelerating the discovery of potent inhibitors against the DENV. In this endeavor, “Anti-Dengue” is the first dedicated web server based on MLTs to find novel potential repurposing drug candidates against DENV infection.
The limitations of the current study are primarily associated with the size of the dataset. Specifically, the relatively small number of entries related to the dengue virus poses a constraint, as a larger dataset could enhance the predictive model’s performance. Another limitation is that the Anti-Dengue web server is currently employing a highly effective SVM-based predictive model for the identification of potential inhibitors/repurposed drugs in terms of inhibition efficiency, as indicated by the pIC50 and IC50 values (μM) against the dengue virus. Unfortunately, alternative machine learning models were not integrated due to their inferior performance on the existing dataset. It is our belief that the development of more robust predictive models using machine learning may be achievable in the future with the availability of additional data. A third limitation is that the “Anti-Dengue” web server is designed exclusively for small molecules, as it is trained on chemicals and FDA-approved drugs, and is not applicable to peptides, antibodies, etc.

5. Conclusions

We developed a QSAR-based algorithm, “Anti-Dengue”(https://bioinfo.imtech.res.in/manojk/antidengue/), which utilizes SVM, ANN, kNN, and RF. Predictive models were developed to identify the potent inhibitors against the DENV. The performance of these predictive models was found to be good, with a PCC of up to 0.71 on the training/testing dataset and a PCC of up to 0.81 on the independent validation dataset. Further applicability domain, chemical clustering, and decoy dataset analyses showed that these predictive models are reliable and robust in nature. The “DrugBank database” was scanned to predict the potential repurposed drug candidates against the DENV. As a result, it will facilitate the rapid development of antivirals that are effective against the DENV.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/v16010045/s1. Table S1: Table showing the drugs/chemicals taken from DrugRepV database targeting Dengue virus used for the development of predictive models; Table S2: Table showing the top 100 selected features for dengue virus from 3 different Recurssive feature elimination techniques i.e., Support Vector regression, Decision Tree regression and Perceptron method; Table S3: The statistical measures of performance of the all the predictive models developed for dengue virus using support vector machine (SVM), random forest (RF), k-nearest neighbour (kNN), artificial neural network (ANN) and deep neural network (DNN) machine-learning techniques utilizing support vector regression (SVR), decision tree regression (DTR) and perceptron method (PCT) based selected features during ten-fold cross validation on five random training/testing and independent validation datasets; Table S4: The data of actual versus predicted pIC50 values used as independent validation set in best performing SVM model; Table S5: The input data for applicability domain analysis of predictive models developed for dengue virus; Table S6: The data of actual versus predicted pIC50 values used in models’ development for dengue virus; Table S7: The data of actual versus predicted pIC50 values of decoy datasets for dengue virus; Table S8: The datasets are clustered into bins using a binning clustering algorithm with similarity cut-off of 0.6 for dengue virus.

Author Contributions

M.K. carried out the Conceptualization, Supervision, Writing—review and editing, Software, Project administration, and Funding acquisition. S.G. performed the Data curation, Methodology, Software, Validation, Formal analysis, Visualization, and Writing—original draft. A.T. performed the Formal analysis, Methodology, and Software. A.R. executed the Methodology, Software, Visualization, and Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge funding support from the CSIR- Institute of Microbial Technology (OLP0143 and STS038).

Data Availability Statement

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Harapan, H.; Michie, A.; Sasmono, R.T.; Imrie, A. Dengue: A Minireview. Viruses 2020, 12, 829. [Google Scholar] [CrossRef] [PubMed]
  2. Gupta, N.; Srivastava, S.; Jain, A.; Chaturvedi, U.C. Dengue in India. Indian J. Med. Res. 2012, 136, 373–390. [Google Scholar] [PubMed]
  3. Ooi, E.-E.; Gubler, D.J. Dengue in Southeast Asia: Epidemiological characteristics and strategic challenges in disease prevention. Cad. Saude Publica 2009, 25 (Suppl. 1), S115–S124. [Google Scholar] [CrossRef] [PubMed]
  4. Bhatt, S.; Gething, P.W.; Brady, O.J.; Messina, J.P.; Farlow, A.W.; Moyes, C.L.; Drake, J.M.; Brownstein, J.S.; Hoen, A.G.; Sankoh, O.; et al. The global distribution and burden of dengue. Nature 2013, 496, 504–507. [Google Scholar] [CrossRef] [PubMed]
  5. Brady, O.J.; Gething, P.W.; Bhatt, S.; Messina, J.P.; Brownstein, J.S.; Hoen, A.G.; Moyes, C.L.; Farlow, A.W.; Scott, T.W.; Hay, S.I. Refining the Global Spatial Limits of Dengue Virus Transmission by Evidence-Based Consensus. PLoS Negl. Trop. Dis. 2012, 6, e1760. [Google Scholar] [CrossRef] [PubMed]
  6. Dwivedi, V.D.; Tripathi, I.P.; Tripathi, R.C.; Bharadwaj, S.; Mishra, S.K. Genomics, proteomics and evolution of dengue virus. Brief. Funct. Genom. 2017, 16, 217–227. [Google Scholar] [CrossRef] [PubMed]
  7. Kalayanarooj, S. Clinical Manifestations and Management of Dengue/DHF/DSS. Trop. Med. Health 2011, 39 (Suppl. 4), 83–87. [Google Scholar] [CrossRef]
  8. Low, J.S.Y.; Wu, K.X.; Chen, K.C.; Ng, M.M.-L.; Chu, J.J.H. Narasin, a novel antiviral compound that blocks dengue virus protein expression. Antivir. Ther. 2011, 16, 1203–1218. [Google Scholar] [CrossRef]
  9. Raekiansyah, M.; Mori, M.; Nonaka, K.; Agoh, M.; Shiomi, K.; Matsumoto, A.; Morita, K. Identification of novel antiviral of fungus-derived brefeldin A against dengue viruses. Trop. Med. Health 2017, 45, 32. [Google Scholar] [CrossRef]
  10. Bardiot, D.; Koukni, M.; Smets, W.; Carlens, G.; McNaughton, M.; Kaptein, S.; Dallmeier, K.; Chaltin, P.; Neyts, J.; Marchand, A. Discovery of Indole Derivatives as Novel and Potent Dengue Virus Inhibitors. J. Med. Chem. 2018, 61, 8390–8401. [Google Scholar] [CrossRef]
  11. Malakar, S.; Sreelatha, L.; Dechtawewat, T.; Noisakran, S.; Yenchitsomanus, P.T.; Chu, J.J.H.; Limjindaporn, T. Drug repurposing of quinine as antiviral against dengue virus infection. Virus Res. 2018, 255, 171–178. [Google Scholar] [CrossRef] [PubMed]
  12. Tafere, G.G.; Wondafrash, D.Z.; Demoz, F.B. Repurposing of N-Acetylcysteine for the Treatment of Dengue Virus-Induced Acute Liver Failure. Hepat. Med. 2020, 12, 173–178. [Google Scholar] [CrossRef] [PubMed]
  13. Shen, T.-J.; Hanh, V.T.; Nguyen, T.Q.; Jhan, M.-K.; Ho, M.-R.; Lin, C.-F. Repurposing the Antiemetic Metoclopramide as an Antiviral against Dengue Virus Infection in Neuronal Cells. Front. Cell. Infect. Microbiol. 2020, 10, 606743. [Google Scholar] [CrossRef] [PubMed]
  14. Botta, L.; Rivara, M.; Zuliani, V.; Radi, M. Drug repurposing approaches to fight Dengue virus infection and related diseases. Front. Biosci. (Landmark Ed.) 2018, 23, 997–1019. [Google Scholar] [CrossRef] [PubMed]
  15. Qureshi, A.; Kaur, G.; Kumar, M. AVCpred: An integrated web server for prediction and design of antiviral compounds. Chem. Biol. Drug Des. 2017, 89, 74–83. [Google Scholar] [CrossRef] [PubMed]
  16. Thakur, N.; Qureshi, A.; Kumar, M. AVPpred: Collection and prediction of highly effective antiviral peptides. Nucleic Acids Res. 2012, 40, W199–W204. [Google Scholar] [CrossRef]
  17. Qureshi, A.; Tandon, H.; Kumar, M. AVP-IC50 Pred: Multiple machine learning techniques-based prediction of peptide antiviral activity in terms of half maximal inhibitory concentration (IC50). Biopolymers 2015, 104, 753–763. [Google Scholar] [CrossRef]
  18. Qureshi, A.; Rajput, A.; Kaur, G.; Kumar, M. HIVprotI: An integrated web based platform for prediction and design of HIV proteins inhibitors. J. Cheminform. 2018, 10, 12. [Google Scholar] [CrossRef]
  19. Rajput, A.; Kumar, M. Anti-flavi: A Web Platform to Predict Inhibitors of Flaviviruses Using QSAR and Peptidomimetic Approaches. Front. Microbiol. 2018, 9, 3121. [Google Scholar] [CrossRef]
  20. Rajput, A.; Thakur, A.; Mukhopadhyay, A.; Kamboj, S.; Rastogi, A.; Gautam, S.; Jassal, H.; Kumar, M. Prediction of repurposed drugs for Coronaviruses using artificial intelligence and machine learning. Comput. Struct. Biotechnol. J. 2021, 19, 3133–3148. [Google Scholar] [CrossRef]
  21. Rajput, A.; Kumar, M. Anti-Ebola: An initiative to predict Ebola virus inhibitors through machine learning. Mol. Divers. 2022, 26, 1635–1644. [Google Scholar] [CrossRef] [PubMed]
  22. Rajput, A.; Kumar, A.; Megha, K.; Thakur, A.; Kumar, M. DrugRepV: A compendium of repurposed drugs and chemicals targeting epidemic and pandemic viruses. Brief. Bioinform. 2021, 22, 1076–7084. [Google Scholar] [CrossRef] [PubMed]
  23. Rajput, A.; Kumar, A.; Kumar, M. Computational Identification of Inhibitors Using QSAR Approach against Nipah Virus. Front. Pharmacol. 2019, 10, 71. [Google Scholar] [CrossRef] [PubMed]
  24. O’Boyle, N.M.; Banck, M.; James, C.A.; Morley, C.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An open chemical toolbox. J. Cheminform. 2011, 3, 33. [Google Scholar] [CrossRef] [PubMed]
  25. Yap, C.W. PaDEL-descriptor: An open source software to calculate molecular descriptors and fingerprints. J. Comput. Chem. 2011, 32, 1466–1474. [Google Scholar] [CrossRef] [PubMed]
  26. Kar, S.; Leszczynski, J. Exploration of Computational Approaches to Predict the Toxicity of Chemical Mixtures. Toxics 2019, 7, 15. [Google Scholar] [CrossRef] [PubMed]
  27. Perkins, R.; Fang, H.; Tong, W.; Welsh, W.J. Quantitative structure-activity relationship methods: Perspectives on drug discovery and toxicology. Environ. Toxicol. Chem. 2003, 22, 1666–1679. [Google Scholar] [CrossRef]
  28. Hira, Z.M.; Gillies, D.F. A Review of Feature Selection and Feature Extraction Methods Applied on Microarray Data. Adv. Bioinform. 2015, 2015, 198363. [Google Scholar] [CrossRef]
  29. Lin, X.; Yang, F.; Zhou, L.; Yin, P.; Kong, H.; Xing, W.; Lu, X.; Jia, L.; Wang, Q.; Xu, G. A support vector machine-recursive feature elimination feature selection method based on artificial contrast variables and mutual information. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2012, 910, 149–155. [Google Scholar] [CrossRef]
  30. Gholami, B.; Norton, I.; Tannenbaum, A.R.; Agar, N.Y.R. Recursive feature elimination for brain tumor classification using desorption electrospray ionization mass spectrometry imaging. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. IEEE Eng. Med. Biol. Soc. Annu. Int. Conf. 2012, 2012, 5258–5261. [Google Scholar] [CrossRef]
  31. Petkovic, D.; Altman, R.; Wong, M.; Vigil, A. Improving the explainability of Random Forest classifier—User centered approach. Pac. Symp. Biocomput. 2018, 23, 204–215. [Google Scholar] [PubMed]
  32. Choi, R.Y.; Coyner, A.S.; Kalpathy-Cramer, J.; Chiang, M.F.; Campbell, J.P. Introduction to Machine Learning, Neural Networks, and Deep Learning. Transl. Vis. Sci. Technol. 2020, 9, 14. [Google Scholar] [PubMed]
  33. Greener, J.G.; Kandathil, S.M.; Moffat, L.; Jones, D.T. A guide to machine learning for biologists. Nat. Rev. Mol. Cell Biol. 2022, 23, 40–55. [Google Scholar] [CrossRef] [PubMed]
  34. Kar, S.; Roy, K.; Leszczynski, J. Applicability Domain: A Step Toward Confident Predictions and Decidability for QSAR Modeling. Methods Mol. Biol. 2018, 1800, 141–169. [Google Scholar] [CrossRef] [PubMed]
  35. Cereto-Massagué, A.; Guasch, L.; Valls, C.; Mulero, M.; Pujadas, G.; Garcia-Vallvé, S. DecoyFinder: An easy-to-use python GUI application for building target-specific decoy sets. Bioinformatics 2012, 28, 1661–1662. [Google Scholar] [CrossRef] [PubMed]
  36. Irwin, J.J.; Tang, K.G.; Young, J.; Dandarchuluun, C.; Wong, B.R.; Khurelbaatar, M.; Moroz, Y.S.; Mayfield, J.; Sayle, R.A. ZINC20—A Free Ultralarge-Scale Chemical Database for Ligand Discovery. J. Chem. Inf. Model. 2020, 60, 6065–6073. [Google Scholar] [CrossRef] [PubMed]
  37. Backman, T.W.H.; Cao, Y.; Girke, T. ChemMine tools: An online service for analyzing and clustering small molecules. Nucleic Acids Res. 2011, 39, W486–W491. [Google Scholar] [CrossRef] [PubMed]
  38. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef]
  39. Kumar, S.; Kumar, G.S.; Maitra, S.S.; Malý, P.; Bharadwaj, S.; Sharma, P.; Dwivedi, V.D. Viral informatics: Bioinformatics-based solution for managing viral infections. Brief. Bioinform. 2022, 23, 1–36. [Google Scholar] [CrossRef]
  40. Niazi, S.K.; Mariam, Z. Recent Advances in Machine-Learning-Based Chemoinformatics: A Comprehensive Review. Int. J. Mol. Sci. 2023, 24, 11488. [Google Scholar] [CrossRef]
  41. Ferdous, N.; Reza, M.N.; Hossain, M.U.; Mahmud, S.; Napis, S.; Chowdhury, K.; Mohiuddin, A.K. Mpropred: A machine learning (ML) driven Web-App for bioactivity prediction of SARS-CoV-2 main protease (Mpro) antagonists. PLoS ONE 2023, 18, e0287179. [Google Scholar] [CrossRef] [PubMed]
  42. Carneiro, J.; Magalhães, R.P.; de la Oliva Roque, V.M.; Simões, M.; Pratas, D.; Sousa, S.F. TargIDe: A machine-learning workflow for target identification of molecules with antibiofilm activity against Pseudomonas aeruginosa. J. Comput. Aided. Mol. Des. 2023, 37, 265–278. [Google Scholar] [CrossRef] [PubMed]
  43. Adams, J.; Agyenkwa-Mawuli, K.; Agyapong, O.; Wilson, M.D.; Kwofie, S.K. EBOLApred: A machine learning-based web application for predicting cell entry inhibitors of the Ebola virus. Comput. Biol. Chem. 2022, 101, 107766. [Google Scholar] [CrossRef] [PubMed]
  44. Malik, A.A.; Chotpatiwetchkul, W.; Phanus-umporn, C.; Nantasenamat, C.; Charoenkwan, P.; Shoombuatong, W. StackHCV: A web-based integrative machine-learning framework for large-scale identification of hepatitis C virus NS5B inhibitors. J. Comput. Aided. Mol. Des. 2021, 35, 1037–1053. [Google Scholar] [CrossRef] [PubMed]
  45. Barba, M.; Dutoit, R.; Legrain, C.; Labedan, B. Identifying reaction modules in metabolic pathways: Bioinformatic deduction and experimental validation of a new putative route in purine catabolism. BMC Syst. Biol. 2013, 7, 1–16. [Google Scholar] [CrossRef] [PubMed]
  46. Kamboj, S.; Rajput, A.; Rastogi, A.; Thakur, A.; Kumar, M. Targeting non-structural proteins of Hepatitis C virus for predicting repurposed drugs using QSAR and machine learning approaches. Comput. Struct. Biotechnol. J. 2022, 20, 3422–3438. [Google Scholar] [CrossRef] [PubMed]
  47. Carro, A.C.; Piccini, L.E.; Damonte, E.B. Blockade of dengue virus entry into myeloid cells by endocytic inhibitors in the presence or absence of antibodies. PLoS Negl. Trop. Dis. 2018, 12, e0006685. [Google Scholar] [CrossRef] [PubMed]
  48. Shahen, M.; Guo, Z.; Shar, A.H.; Ebaid, R.; Tao, Q.; Zhang, W.; Wu, Z.; Bai, Y.; Fu, Y.; Zheng, C.; et al. Dengue virus causes changes of MicroRNA-genes regulatory network revealing potential targets for antiviral drugs. BMC Syst. Biol. 2018, 12, 2. [Google Scholar] [CrossRef]
  49. Boonyasuppayakorn, S.; Reichert, E.D.; Manzano, M.; Nagarajan, K.; Padmanabhan, R. Amodiaquine, an antimalarial drug, inhibits dengue virus type 2 replication and infectivity. Antivir. Res. 2014, 106, 125–134. [Google Scholar] [CrossRef]
  50. Punekar, M.; Kasabe, B.; Patil, P.; Kakade, M.B.; Parashar, D.; Alagarasu, K.; Cherian, S. A Transcriptomics-Based Bioinformatics Approach for Identification and In Vitro Screening of FDA-Approved Drugs for Repurposing against Dengue Virus-2. Viruses 2022, 14, 2150. [Google Scholar] [CrossRef]
  51. Kumar, S.; Bajrai, L.H.; Faizo, A.A.; Khateb, A.M.; Alkhaldy, A.A.; Rana, R.; Azhar, E.I.; Dwivedi, V.D. Pharmacophore-Model-Based Drug Repurposing for the Identification of the Potential Inhibitors Targeting the Allosteric Site in Dengue Virus NS5 RNA-Dependent RNA Polymerase. Viruses 2022, 14, 1827. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The workflow includes retrieving dengue inhibitors from DrugRepV and converting SMILES to SDF format. Molecular descriptors/fingerprints are calculated using PaDEL software, followed by the recursive feature elimination (RFE) module for feature selection. SVM, ANN, kNN, and RF MLTs are employed with ten-fold cross-validation for predictive algorithms. The performance is evaluated using MAE, MSE, RMSE, R2, and PCC values. Further, the model’s robustness is analyzed with applicability domain, scatter plots, and decoy sets. Potent repurposed drugs are predicted by scanning the “DrugBank” database.
Figure 1. The workflow includes retrieving dengue inhibitors from DrugRepV and converting SMILES to SDF format. Molecular descriptors/fingerprints are calculated using PaDEL software, followed by the recursive feature elimination (RFE) module for feature selection. SVM, ANN, kNN, and RF MLTs are employed with ten-fold cross-validation for predictive algorithms. The performance is evaluated using MAE, MSE, RMSE, R2, and PCC values. Further, the model’s robustness is analyzed with applicability domain, scatter plots, and decoy sets. Potent repurposed drugs are predicted by scanning the “DrugBank” database.
Viruses 16 00045 g001
Figure 2. The applicability domain analysis of the support vector machine was assessed by a William’s plot between the leverage and standardized residuals of the molecules.
Figure 2. The applicability domain analysis of the support vector machine was assessed by a William’s plot between the leverage and standardized residuals of the molecules.
Viruses 16 00045 g002
Figure 3. The robustness of the (a) support vector machine, (b) artificial neural network, (c) k-nearest neighbor, and (d) random forest-based predicted models was assessed by scatter plots between the actual and predicted pIC50 values of the molecules.
Figure 3. The robustness of the (a) support vector machine, (b) artificial neural network, (c) k-nearest neighbor, and (d) random forest-based predicted models was assessed by scatter plots between the actual and predicted pIC50 values of the molecules.
Viruses 16 00045 g003
Figure 4. To evaluate the reliability of the predicted models based on SVM, a scatter plot was generated to compare the actual and decoy pIC50 values of (a) Decoy Set 1, (b) Decoy Set 2, (c) Decoy Set 3, (d) Decoy Set 4, (e) Decoy Set 5, and (f) Decoy Set 6.
Figure 4. To evaluate the reliability of the predicted models based on SVM, a scatter plot was generated to compare the actual and decoy pIC50 values of (a) Decoy Set 1, (b) Decoy Set 2, (c) Decoy Set 3, (d) Decoy Set 4, (e) Decoy Set 5, and (f) Decoy Set 6.
Viruses 16 00045 g004
Figure 5. Chemical diversity analysis: (a) 2D multidimensional scaling plot and (b) 3D multidimensional scaling plot of the anti-dengue compounds.
Figure 5. Chemical diversity analysis: (a) 2D multidimensional scaling plot and (b) 3D multidimensional scaling plot of the anti-dengue compounds.
Viruses 16 00045 g005
Table 1. “Anti-Dengue” predictive model performances during 10-fold cross-validation using the SVM machine learning technique.
Table 1. “Anti-Dengue” predictive model performances during 10-fold cross-validation using the SVM machine learning technique.
AlgorithmFeature SelectionModel ParametersDatasetRMSEMSEMAER2PCC
SVMPerceptronsvm_param_32_kernel_rbf_gamma_0.005_C_10T2140.690.470.480.470.71
V240.430.190.360.560.81
SVMSVRsvm_param_32_kernel_rbf_gamma_0.005_C_10T2140.720.550.510.390.68
V240.380.150.310.660.84
SVMDTsvm_param_1_kernel_rbf_gamma_0.1_C_0.01T2140.970.990.71−0.070.41
V240.650.420.560.020.36
Table 2. “Anti-Dengue” predictive model performances during 10-fold cross-validation using the ANN machine learning technique.
Table 2. “Anti-Dengue” predictive model performances during 10-fold cross-validation using the ANN machine learning technique.
AlgorithmFeature SelectionModel ParametersDataset RMSEMSEMAER2PCC
ANNPerceptronANN__paras_19_activation_identity_solver_sgd_learning_constantT2140.720.590.520.040.65
V240.580.330.460.220.74
ANNSVRANN__paras_26_activation_identity_solver_lbfgs_learning_invscalingT2140.520.320.360.620.67
V240.340.110.260.740.90
ANNDTANN__paras_14_activation_relu_solver_adam_learning_invscalingT2144.63108.311.77−98.10.45
V240.90.820.63−0.90.43
Table 3. “Anti-Dengue” predictive model performances during 10-fold cross-validation using the kNN machine learning technique.
Table 3. “Anti-Dengue” predictive model performances during 10-fold cross-validation using the kNN machine learning technique.
AlgorithmFeature SelectionModel Parameters Dataset RMSEMSEMAER2PCC
kNNPerceptronknn_k9T2140.890.830.70.00.34
V240.50.250.410.410.68
kNNSVRknn_k7T2140.870.810.670.070.35
V240.460.210.370.510.74
kNNDTknn_k9T2140.90.880.660.020.37
V240.480.230.380.460.72
Table 4. “Anti-Dengue” predictive models performance during 10-fold cross-validation using RF machine learning technique.
Table 4. “Anti-Dengue” predictive models performance during 10-fold cross-validation using RF machine learning technique.
AlgorithmFeature SelectionModel Parameters Dataset RMSEMSEMAER2PCC
RFPerceptronrf__paras_30_n_200_depth_12_split_5_leaf_4T2140.890.820.660.070.45
V240.570.330.470.240.54
RFSVRrf__paras_44_n_300_depth_12_split_2_leaf_2T2140.840.760.630.150.49
V240.450.20.360.540.79
RFDTrf__paras_30_n_200_depth_12_split_5_leaf_4T2140.840.740.610.130.54
V240.450.20.370.530.77
PCC—Pearson’s correlation coefficient, R2—coefficient of determination, MAE—mean absolute error, MSE—mean squared error, and RMSE—root mean squared error.
Table 5. The top hits of the predicted repurposed drug candidates.
Table 5. The top hits of the predicted repurposed drug candidates.
DrugBankIDDrug NamePrimary IndicationPredicted_pIC50Status
DB00014GoserelinBreast cancer and prostate cancer8.42 Not yet tested
DB00644GonadorelinFunction of gonadotropes and
the pituitary
8.19 Not yet tested
DB00666NafarelinCentral precocious puberty in
children of both sexes and
treatment of endometriosis
8.03 Not yet tested
DB11279Brilliant greenTo prevent infections of the
umbilical cord
8.03 Not yet tested
DB01284TetracosactideScreening of patients presumed
to have adrenocortical
insufficiency
7.91 Not yet tested
DB12887TazemetostatMetastatic or locally advanced
epithelioid sarcoma is not
eligible for complete resection.
7.83 Not yet tested
DB00626BacitracinWound infections, pneumonia,
skin and eye infections
7.83 Not yet tested
DB01061AzlocillinPseudomonas aeruginosa,
Haemophilus influenzae and
Escherichia coli infections
7.81 Not yet tested
DB01403MethotrimeprazineFor the treatment of psychosis,
particular those of
schizophrenia, and manic
phases of bipolar disorder
7.8 Not yet tested
DB01621PipotiazineChronic non-agitated
schizophrenic patients
7.67 Not yet tested
DB01147CloxacillinTreatment of beta-hemolytic
streptococcal, pneumococcal,
and staphylococcal infections
7.67 Not yet tested
DB06788HistrelinPalliative treatment of
advanced prostate cancer
7.65 Not yet tested
DB09320Procaine
benzylpenicillin
Local anesthetic and antibiotic
combination for bacterial infections
7.62 Not yet tested
DB00434CyproheptadineAppetite stimulation, allergic
symptoms, and treatment of
serotonin syndrome
7.51 Not yet tested
DB09570IxazomibMultiple myeloma7.51Not yet tested
DB09473Indium In-111
Oxyquinoline
Radiolabeling autologous
leukocytes
7.5 Not yet tested
DB04826ThenalidineNot available7.41Not yet tested
DB00477ChlorpromazinePreoperative anxiety,
nausea, vomiting, bipolar
disorder, and schizophrenia
7.27 Experimental
DB00948MezlocillinLungs, urinary tract, skin
gram-negative infections
7.39 Not yet tested
DB01201RifapentinePulmonary tuberculosis7.39 Not yet tested
DB00455LoratadineManage the symptoms of
allergic rhinitis
6.8 Experimental
DB01087PrimaquineTo prevent relapse of vivax
Malaria
6.69 Experimental
DB00468QuinineUncomplicated
Plasmodium falciparum
Malaria
6.65 Experimental
DB01583LiotrixPrimary, secondary or tertiary hypothyroidism6.63Not yet tested
DB09225ZotepineSchizophrenia6.63Not yet tested
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gautam, S.; Thakur, A.; Rajput, A.; Kumar, M. Anti-Dengue: A Machine Learning-Assisted Prediction of Small Molecule Antivirals against Dengue Virus and Implications in Drug Repurposing. Viruses 2024, 16, 45. https://doi.org/10.3390/v16010045

AMA Style

Gautam S, Thakur A, Rajput A, Kumar M. Anti-Dengue: A Machine Learning-Assisted Prediction of Small Molecule Antivirals against Dengue Virus and Implications in Drug Repurposing. Viruses. 2024; 16(1):45. https://doi.org/10.3390/v16010045

Chicago/Turabian Style

Gautam, Sakshi, Anamika Thakur, Akanksha Rajput, and Manoj Kumar. 2024. "Anti-Dengue: A Machine Learning-Assisted Prediction of Small Molecule Antivirals against Dengue Virus and Implications in Drug Repurposing" Viruses 16, no. 1: 45. https://doi.org/10.3390/v16010045

APA Style

Gautam, S., Thakur, A., Rajput, A., & Kumar, M. (2024). Anti-Dengue: A Machine Learning-Assisted Prediction of Small Molecule Antivirals against Dengue Virus and Implications in Drug Repurposing. Viruses, 16(1), 45. https://doi.org/10.3390/v16010045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop