Next Article in Journal
Intermolecular Mechanism and Dynamic Investigation of Avian Influenza H7N9 Virus’ Susceptibility to E119V-Substituted Peramivir–Neuraminidase Complex
Next Article in Special Issue
Rac1: A Regulator of Cell Migration and a Potential Target for Cancer Therapy
Previous Article in Journal
Application of New Energy Thermochromic Composite Thermosensitive Materials of Smart Windows in Recent Years
Previous Article in Special Issue
Identification of New Regulators of Pancreatic Cancer Cell Sensitivity to Oxaliplatin and Cisplatin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Machine Learning and Chemoinformatics-Based Screening of Mycotic Compounds against Kinesin Spindle ProteinEg5 for Lung Cancer Therapy

1
Centre for Environmental Assessment and Climate Change, G.B. Pant National Institute of Himalayan Environment (GBP-NIHE), Kosi-Katarmal, Almora 263643, Uttarakhand, India
2
Department of Botany, DSB Campus, Kumaun University, Nainital 263002, Uttarakhand, India
3
ENVIS Centre on Himalayan Ecology, G.B. Pant National Institute of Himalayan Environment (GBP-NIHE), Kosi-Katarmal, Almora 263643, Uttarakhand, India
4
Centre for Biodiversity Conservation and Management, G.B. Pant National Institute of Himalayan Environment (GBP-NIHE), Kosi-Katarmal, Almora 263643, Uttarakhand, India
5
ICAR-Indian Veterinary Research Institute, Bengaluru 560024, Karnataka, India
6
Department of Biotechnology, Bhimtal Campus, Kumaun University, Nainital 263136, Uttarakhand, India
7
Genetic Engineering and Biotechnology, University of Rajshahi, Rajshahi 6205, Bangladesh
8
Nutrition and Bromatology Group, Department of Analytical Chemistry and Food Science, Faculty of Science, Universidade de Vigo, E-32004 Ourense, Spain
9
Department of Pharamaceutics, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
10
Department of Pharmacy Practice, College of Pharamcy, AlMaarefa University, Ad Diriyah 13713, Saudi Arabia
11
Department of Clinical Laboratory Science, Faculty of Applied Medical Science, Najran University, Najran 61441, Saudi Arabia
12
Department of Botany, Soban Singh Jeena University, Almora 263601, Uttarakhand, India
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(5), 1639; https://doi.org/10.3390/molecules27051639
Submission received: 26 January 2022 / Revised: 11 February 2022 / Accepted: 17 February 2022 / Published: 2 March 2022
(This article belongs to the Special Issue Targeting of Signaling Pathways for Cancer Therapy)

Abstract

:
Among the various types of cancer, lung cancer is the second most-diagnosed cancer worldwide. The kinesin spindle protein, Eg5, is a vital protein behind bipolar mitotic spindle establishment and maintenance during mitosis. Eg5 has been reported to contribute to cancer cell migration and angiogenesis impairment and has no role in resting, non-dividing cells. Thus, it could be considered as a vital target against several cancers, such as renal cancer, lung cancer, urothelial carcinoma, prostate cancer, squamous cell carcinoma, etc. In recent years, fungal secondary metabolites from the Indian Himalayan Region (IHR) have been identified as an important lead source in the drug development pipeline. Therefore, the present study aims to identify potential mycotic secondary metabolites against the Eg5 protein by applying integrated machine learning, chemoinformatics based in silico-screening methods and molecular dynamic simulation targeting lung cancer. Initially, a library of 1830 mycotic secondary metabolites was screened by a predictive machine-learning model developed based on the random forest algorithm with high sensitivity (1) and an ROC area of 0.99. Further, 319 out of 1830 compounds screened with active potential by the model were evaluated for their drug-likeness properties by applying four filters simultaneously, viz., Lipinski’s rule, CMC-50 like rule, Veber rule, and Ghose filter. A total of 13 compounds passed from all the above filters were considered for molecular docking, functional group analysis, and cell line cytotoxicity prediction. Finally, four hit mycotic secondary metabolites found in fungi from the IHR were screened viz., (−)-Cochlactone-A, Phelligridin C, Sterenin E, and Cyathusal A. All compounds have efficient binding potential with Eg5, containing functional groups like aromatic rings, rings, carboxylic acid esters, and carbonyl and with cell line cytotoxicity against lung cancer cell lines, namely, MCF-7, NCI-H226, NCI-H522, A549, and NCI H187. Further, the molecular dynamics simulation study confirms the docked complex rigidity and stability by exploring root mean square deviations, root mean square fluctuations, and radius of gyration analysis from 100 ns simulation trajectories. The screened compounds could be used further to develop effective drugs against lung and other types of cancer.

Graphical Abstract

1. Introduction

Cancer is the second-leading cause of death globally after cardiovascular disease, and it is accountable for a 19.3 million increase in cases with 10 million deaths in 2020. Worldwide, the most commonly diagnosed cancer is female breast (11.7%), followed by lung cancer (~11.4%), colorectal cancers (10.0%), liver cancer (8.3%), stomach cancer (7.7%), and prostate cancer (7.3%) [1]. Being the second-leading cause of cancer-related deaths in developed countries, lung cancer is mounting at alarming rates in developing countries and low-middle-income countries, responsible for ~1.3 million deaths per year [2]. The over-expression of Eg5 has been reported in lung cancer and has an essential role in cell migration and angiogenesis impairment [3,4,5,6]. In this context, the mitotic kinesin, Eg5, is a promising target as it has a critical role in assembling the mitotic spindle during cell division. Kinesin superfamily proteins or KIF11/Eg5 are microtubule-based motor proteins that generate directional movement along microtubules. KIFs are core proteins, and these are not only essential for intracellular transport but are also important for various cellular and morphology functions [7]. KIF11/Eg5 is also identified as a prognostic factor. Overall, it indicates poor survival and the worst progression-free survival in lung cancer patients. Its knockdown is reported for induced G2/M phase arrest and improved apoptosis in lung cancer cells [8]. Various studies reported that the natural compounds inhibit the abnormal Eg5 effect and minimize tumors in patients suffering from lung cancer [9,10]. Therefore, Eg5 is a good target for anti-cancer drug development, and natural products may provide a huge source of compounds with therapeutic potential. With the above context, several fungal-derived natural products, such as lovastatin, echinocandin B, and cyclosporine A, are well reported with excellent therapeutic potential. Currently, several fungal secondary metabolites from the Indian Himalayan Region (IHR) have been identified as a vital source of potent therapeutic compounds. A diverse structural class of secondary metabolites, including aromatic compounds, amino acids, anthacenones, butanolides, butenolides, cytochalasans, macrolides, naphthalenones, pyrones, and terpenes, are produced by fungi [11]. Therefore, in the present study, mycotic secondary metabolites are screened for their anti-lung cancer property targeting the Eg5 protein. Several in silico techniques were used for screening, based on machine learning, chemo-informatics, and molecular dynamics simulation.
Currently, machine learning (ML) techniques are widely used in the cancer drug discovery process, including structure-based molecular docking, deep learning, and ligand-based cheminformatics modeling. ML also applies to proteo-chemometrics modeling, cell phenotype data, transcriptomics, and electronic health records (EHRs) based ML [12]. Drug development using artificial intelligence (Al), such as ML, is less costly and more time effective than the rational drug discovery process. In the drug discovery market, Al is projected to reach USD 1434 million by 2024 from USD 259 million in 2019 at a compound annual growth rate (CAGR) of 40.8% [13]. The machine learning-based quantitative structure-activity relationships (QSAR) models are widely used nowadays for screening huge compounds to predict potential inhibitors against diverse diseases. These types of models codify the chemical 3D structure based on molecular descriptors, such as constitutional, functional groups, geometrical, thermodynamic, topological, and quantum mechanical. However, the development of a QSAR model with high accuracy is a challenging task that depends on the correct analysis and selection of computed descriptors, application of proper classifier, and proper statistical understating of the model [14]. In the present work, machine-learning-based QSAR modeling, molecular docking, and other chemoinformatics-based screening techniques and molecular dynamics simulation were employed for screening natural inhibitors against the Eg5 protein targeting lung cancer.

2. Methods

2.1. Machine Learning Model Development and Screening

In the present work, the complete PubChem Human A549 Lung Tumor Cell Growth Inhibition Assay (AID: 317) data set of percentage inhibition for the cancer line was considered for model development [15]. The assay consisted of 3317 compounds in total, out of which for model development, 278 compounds (50% growth inhibition at 10 uM) were considered as active, and 294 compounds with negative or not having growth inhibition capacity were considered as inactive (Supplementary Materials, Table S1).
For screening, a structural library of 1830 mycotic secondary metabolites from 184 medicinal fungi was retrieved from the Medicinal Fungi Secondary Metabolite And Therapeutics (MeFSAT) database [16]. All compounds within the model and screening data set were further converted from the three-dimensional (3D) standard data format (SDF) to the molecular-input line-entry system (SMILES) format with the help of the O’babel software [17]. Thereafter, molecular descriptors were calculated for the compounds. Molecular descriptors generally encode various information about molecular fragmentation. A total of 2633 molecular descriptors were obtained, including one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) fingerprint features using the PaDEL-descriptor software [18]. Thereafter, both datasets, i.e., the CHEMBL dataset (AID 371) and the MeFSAT dataset were subjected to the Waikato Environment for Knowledge Analysis (WEKA) software [19] for building the predictive machine-learning model and screening. The CfsSubsetEval module implemented in WEKA was used for the significant descriptor selection. The descriptors with zero values were excluded followed by the removal of highly correlated descriptors. A total 16 descriptors from 8 descriptor types selected by the CFS Subset Evaluator were considered for the model development. Selected descriptor types include autocorrelation (ATSC6i, AATSC8s, GATS7c, GATS1e, GATS1p, GATS2i), basic group count (nBase), burden-modified eigenvalues (SpMin3_Bhe, SpMin1_Bhs), crippenLogP, atom type electrotopological state (nHsOH, SssNH, maxHCsatu), information content (BIC2), ring count (nFRing), and XLogP. Among them are the autocorrelation descriptors pool chemical information given by property values in specified molecule regions and structural information. They are based on a conceptual segmentation of the molecular structure and the application of an autocorrelation function to molecular properties of different molecular regions [20]. The second descriptor class is nBase, which counts the number of basic groups in the compounds that tend to be protonated in the gastrointestinal tract. Basic compounds have higher polarity and low lipophilicity that limits passive absorption across bio-membranes [21]. Further, the third descriptor class, Burden modified eigenvalues, derived from the eigenvalues of a modified adjacency matrix, is chemically intuitive in that its elements relate to atomic and bonding properties of the compounds [22]. The fourth, CrippenLogP descriptor class, simply determines the log P of the compounds by fitting an extensive training set of 9920 molecules with r2 of 0.918 and σ of 0.677 [23]. The fifth descriptor class, atom type electrotopological state descriptors, characterized by atoms in molecules, is introduced as the electrotopological state index, which combines both the electronic character and the topological environment of each skeletal atom in a molecule. They help to identify atoms and regions in the molecule which are important for activity [24]. The BIC2 descriptors are vital for the count of bond information content index (neighborhood symmetry of 2-order). The second to last descriptor, nFRing, gives information about the number of fused rings in a compound. The last descriptor class, XLogP, calculates partition coefficients of solutes in octanol/water responsible for lipophilicity of compounds, a vital parameter responsible for its ability to cross the cell membrane [25].
To develop a machine-learning model the PubChem dataset was divided into 70:30 as a training and test set. The training set was made up of 374 compounds while the test set had 161 compounds (Supplementary Materials, Table S1). Five machine-learning algorithms, namely, random forest, J48, decision stump, random tree, and REPTree classifier with bagging (Bootstrap Aggregation) ensemble algorithm were applied individually to get the best screening model. Random forest is one of the most frequently used ensemble-learning algorithms in the field of in silico drug discovery. It generates multiple decision trees and for one test instance, and the classification of each tree is regarded as one vote, and the test descriptors were confirmed by integrating the prediction results of each tree (sub-classifier) [26]. In J48, the main root of the decision tree is established based on the gain and gain ratio values of the descriptors. In the decision stump algorithm, the classification process is performed by considering only one key feature in the sample set and has root node and leaves. In a random tree algorithm, trees are constructed using randomly selected samples from the descriptor data set and a random one is selected. Whereas, the REPTree creates rapid multiple decision trees, using information gain, prunes them using reduced error pruning and selects the best tree among all generated trees [27]. The accuracy of the models was evaluated by using the following statistical parameters [28]:
Binary   classification   accuracy = TN + TP TN + TP + FN + FP ,   Sensitivity = TP FN + TP ,   Specificity = TN FP + TN ,
Accuracy _ random = TP + FN × TP + FP + TN + FN × TN + FP N ^ 2
Matthews   correlation   coefficient   MCC = TN × TP FN × FP FP + TP FN + TP FP + TN FN + TN ,  
TP = true positives; FP = false positives; TN = true negatives; FN = false negatives.
The best machine model was used as a predictive model to screen anti-lung cancer compounds from the library of 1830 mycotic secondary metabolites.

2.2. Drug-Likeness Screening and Molecular Docking

The DruLiTo-screening tool was further used to screen the ligands with drug-likeness properties [29]. Chemical descriptors such as atom molar refractivity (AMR), AlogP, logP, H-bond acceptor (HBA), H-bond donor (HBD), molecular weight (MW), number of rotatable bonds (nRB), number of atoms, number of acidic groups, rotatable bond count (RC), number of rigid bonds (nRigidB), nAtomRing, nHB, and total polar surface area (TPSA) for all of the ligands were calculated using the chemistry development kit (CDK) [30].
DruLiTo screened the ligands based on the various drug-likeness rules, such as Lipinski’s rule, CMC-50-like rule, Veber rule, and Ghose filter [31]. After screening through the machine-learning model, screened compounds were subjected to target-specific virtual screening. For that, the three-dimensional structure of the target protein, Eg5 (PDB id: 4BXN), was downloaded from an online protein data bank portal [32]. The target protein and machine-learning screened compounds were docked using the PyRx graphic user interface (GUI), which utilizes the AutoDockVina virtual screening program [31]. All screened ligands were docked in the active site of the Eg5 protein with coordinates of X = 11.82, Y = 37.46, and Z = −9.79 to obtain protein–ligand complexes. The compound, N-(3-aminopropyl)-N-[(1R)-1-(3-benzyl-7-chloro-4-oxo-4H-chromen-2-yl)-2-methylpropyl]-4, methylbenzamide, was considered as a reference compound for examining the necessary binding score. Further, the docked protein and ligand complexes were visualized using UCSF chimera [33]. The in-depth interactions were investigated using the protein–ligand interaction profiler server (PLIP) [34].

2.3. Functional Group Analysis and Cell Line Cytotoxicity Prediction

The functional group analysis method was used to functionalize the physical and chemical properties of small groups and atoms that exhibit a distinctive chemical nature when in their original form. The functional group’s frequencies of the screened compounds were analyzed in ‘R’ (version 3.4.3) [35] using the ChemMineR [36] package followed by standard deviation analysis. Five established inhibitors against Eg5 (CID: 132472242, 53494981, 11609157, 11553595, and 11521919) were also subjected to function group analysis for comparison. Functional groups, such as aromatic rings, rings, carboxylic acid esters, and carbonyl, were evaluated for both screened inhibitors and reference inhibitor groups [37]. Thereafter, virtually screened compounds from machine learning, docking, and functional group analysis were additionally executed to evaluate their cell line cytotoxicity prediction. This analysis was carried out by using CLC-Pred web facilities (CLC-Pred: in silico prediction of cytotoxicity for tumor and non-tumor cell lines (way2drug.com, accessed on 8 June 2021)). This estimation of web services was conducted based on the QSAR model build on the prediction of activity spectra for substances (PASS) tools (http://www.way2drug.com/PASSonline, accessed on 8 June 2021) and the training dataset produced based on ChEMBLdb cytotoxicity data [38,39].

2.4. Molecular Dynamics Simulation

The docked complexes of the final screened compounds were further checked for their binding stability to MD simulations using GROMACS 5.0 [40]. The CGenFF server was used to generate topologies and coordinate files. The CHARMM 36 force field was used for topologies of the complex structures of the ligand molecules [41]. The SPC/E water mode l40 was used for solvation of the complexes followed by the addition of counterions for neutralization. The simulation was carried out up to 100 ns time in NVT (constant volume), as well as NPT (constant pressure), maintaining the temperature at 300 K and pressure 1 bar. Finally, the stability of the protein–ligand complexes was analyzed by considering parameters, such as RMSD (root mean square deviations) and radius of gyration (Rg).

3. Results

3.1. Performance Comparison of Different ML Classifiers and Data Set Screening

The performance of the ML algorithms varies with the type of training data set of the model. Therefore, for identification of the best-performing algorithm on the lung cancer inhibitor data set, five classifiers were evaluated for their classification accuracy, namely, random forest, J48, decision stump, random tree, and REP tree classifier with bagging. The performance of the models was assessed with the evaluation of various statistical parameters, viz., correctly classified instances, incorrectly classified instances, Kappa statistic, mean absolute error, root mean square error, Matthews correlation coefficient (MCC), and receiver operating characteristic curve (ROC) area (Figure 1, Table 1). Among all the classifiers, the random forest classifier with the highest sensitivity (1), Kappa statistics of 0.94, and ROC of 0.99 area in the training set was selected for the final screening model (Supplementary Materials, Table S2). The model also performed best in the test set (Supplementary Materials, Table S3). It was able to classify 97.05% of instances of the supplied training set data correctly with a mean out-of-bag (OOB) error of 0.07 (MAE). In the bagging (REPTree) model out-of-bag estimates, 96.52% instances were correctly classified with 0.18 root mean squared error (RMSE). However, bagging (REPTree) and decision stump showed lower sensitivity than random forest. The random tree classifier showed the lowest performance with sensitivity of 0.90, specificity of 0.92, and accuracy of 0.09 (K = 0–4, randomly chosen attributes at each node). The random forest classifier consists of a collection of tree-structured classifiers {h(x, Θ k), k = 1…}, where each tree creates a unit vote against the most popular class at input xand the {Θ k} are independent identically distributed random vectors [42]. The model of the current work was obtained based on the default mtry (number of randomly drawn candidate variables) for classification, square root of the number of descriptors, p = 16, and 100 trees in the forest (ntree). The “forest” built by the model, is an ensemble of decision trees, usually trained with the “bagging” method. The bagging algorithm simulates a defined method before applying a certain training set to remove the instability of learning techniques. The initial training set data is improved by removing some samples and duplicating others in its place of sampling a new training dataset each time by random sampling. The process of sampling necessarily rejects some of the instances and duplicates the others. Bagging simply resamples the native training data instead of creating autonomous data sets from the domain. It makes a combined model that mostly performs better than the single model made from the native training data [43]. The selected model classifies 319 mycotic secondary metabolites with active potential out of 1830 compounds.

3.2. Drug Likeness and Molecular Docking

Screened compounds from ML were further checked for their drug-likeness property by applying four pharmacological filters, namely, Lipinski’s rule, CMC-50-like rule, Veber rule, and Ghose filter. According to Lipinski’s rule, a drug-like compound should have LogP ≤ 5, molecular weight ≤ 500, number of hydrogen bond acceptors ≤ 10, and number of hydrogen bond donors ≤ 5. Whereas, according to the CMC-50-like rule, a drug-like compound must have logP (ALOGP) between 1.3 and 4.1, molar refractivity (AMR) between 70 and 110, molecular weight (MW) between 230 and 390, and the number of atoms (nAT) between 30 and 55. To qualify the Veber rule and Ghose filter, compounds should have rotatable bonds < 12 and polar surface area < 140, logP in the −0.4 to +5.6 range, molar refractivity from 40 to 130, molecular weight from 180 to 480, and the number of atoms from 20 to 70, respectively. A total of 13 compounds out of 319 (Supplementary Materials, Table S4) compounds were able to ensure all the above filters and were considered for molecular docking to check their binding potential with the Eg5 protein (Table 2 and Figure 2). Of those 12, compounds having the lowest binding energy (highest inhibitory potency) were screened further based on the binding energy of the reference compound.

3.3. Functional Group Cell Line Cytotoxicity and Binding Pattern of Screened Molecules

During the analysis of the functional group frequency among some of the established inhibitors of Eg5, viz., the presence of aromatic rings, rings, carboxylic acid esters, and carbonyl groups were shown. A similar frequency pattern was observed in the screened compounds (Figure 3), which reflects their excellent Eg5 inhibitory property. The high frequency of alcohol (ROH), ketone (RCOR), esters (RCOOR), and ether (ROR) groups were analyzed in screened inhibitors as compared to reference inhibitors (Figure 3). The groups, such as alcohol (ROH), ketone (RCOR), and esters (RCOOR), indicate more H-bonding capacity of the compounds or more water solubility. Simultaneously, screened compounds also exhibited a high presence of rings and aromatic rings that prove their strong drug-likeness potential. Aromatic rings commonly promote Van der Waals interactions with the binding site atoms of the protein [44]. Further, the inhibition potential of the compounds in cell lines was checked, and four compounds, namely, (-)-Cochlactone A, Phelligridin C, Sterenin E, and Cyathusal A, were predicted to be active against the lung cancer cell lines, NCI-H226, NCI-H522, A549, and NCI-H187 (Supplementary Materials, Table S5), respectively. The binding profiles of the screened compounds were also checked and all of the above compounds were found with efficient binding potential and can form diverse bonds, including hydrophobic interactions, H-bonds, salt bridges, and π–cation interactions (Table 3, Figure 4). Hydrophobic interactions are vital as they have an energetic contribution to the stability of the protein–ligand complex [45], whereas H-bonds are necessary for enzymatic catalysis to stabilize ligands in the binding pocket [46]. All the screened compounds were found to be able to bond within the switch II (Leu266–Asn289) and α-helix 4 (Gln290–Val303) regions of the allosteric site of the Eg5 protein. ATP-competitive inhibitors of Eg5, including the biphenyl derivatives, have been reported to bind with the protein between the α4 and α6 helices. In Eg5, switch I forms repeated contacts with the γ-phosphate of ATP while switch II offered a β-hairpin conformation [47].

3.4. Binding Stability Analysis of the Screened Compounds

Finally, the structural stability of the screened ligands with the protein was studied by molecular dynamics (MD) simulations. MD simulation was used to predict the movement of every atom in a protein and protein ligand complex over time, based on the general physics model of the governing interatomic interactions. The structural compactness of the screened compounds was assessed by calculating the radius of gyration (Rg) defined as the distribution of atoms of a protein around its axis. All the systems, such as protein and the four complexes, showed stability in the MD simulation approach towards 100 ns time scale between 1.7 and 1.8 nm, reflecting the high structural compactness of the protein, as well as the inhibitor. The lower degree of fluctuation with its consistency throughout the simulation designates more compactness and rigidity of the protein ligand complex. The Rg of Cochlactone A is found to be almost stable in terms of consistency of fluctuations throughout the simulation. Deviations of the conformational stability of macromolecules from the backbone structure to the initial structure were evaluated to check their stability by calculating the root mean square deviation (RMSD). All complexes showed RMSD values between 0.31 and 0.37 nm, and the value starts to increase after 10 ns. The compound, Sterenin E, showed a decreasing pattern after 80 ns. Further, the flexibility of proteins along with all the complexes throughout the simulation process was evaluated by calculating the root mean square fluctuations (RMSF). It indicates the occurrence of local changes along with the protein structure at the selected temperature and pressure. The Cyathusal A_EG5 complex showed the highest RMSF (0.14 nm) values followed by other screened compounds. It can also be understood from Figure 5 that the binding of Cyathusal A makes the protein most flexible in all areas in contrast to the other complexes (Supplementary Materials Table S6, Figure 5).

4. Discussion

The IHR is well known in the field of drug discovery as a treasure house for numerous natural compounds with extensive therapeutic potential. However, the endophytic microbial biodiversity of the area is less explored to date. The endophytic fungus isolated from the medicinal plants of IHR contain several anti-cancer compounds, e.g., Taxol found in Annulohypoxylon sp. Isolated from the Himalayan yew (Taxus wallichiana Zucc.), Vinblastine produced by Curvularia verruculosa, and Vincristine by Fusarium oxysporum isolated from Catharanthus roseus, Podophyllotoxin from Phialocephala fortinii, Camptothecin from Entrophospora infrequens from the plant Nothapodytes nimmoniana etc. [48,49].
Based on the present work, four mycotic secondary metabolites were screened, namely, (−)-Cochlactone A, Phelligridin C, Sterenin E, and Cyathusal A. (−)-Cochlactone A is an alkyl-phenyl ketone, an aromatic compound containing a ketone replaced by one alkyl group and a phenyl group. In 2018, it was found in the fungus Ganoderma cochlear, an edible and medicinal fungus. The compound has a bicycle [4.4.0] decane ring system with a γ-lactone fragment and possesses anti-inflammatory activity [50]. Other species, such as Ganoderma lucidum from the same genus, is well known as a medicinal fungus from the IHR [51]. The second screened compound, Phelligridin C, is a furo [3, 2-c] pyran-4-one derivative reported from the fungus Phellinus igniarius [52]. It is recorded to have cytotoxicity against a human lung cancer cell line. The genus is also reported from the IHR of the Himachal Pradesh and Uttarakhand states of India [53]. The third screened compound, Sterenin E, is a polyketide that has α-glucosidase inhibitory activity from the mushroom, Stereum hirsutum [54]. Naturally, the fungi are known from the Hmuifang forest and Tanhril forest of the IHR from Mizoram, Northeast India [55]. The fourth screened compound, Cyathusal A, is recorded from the mushroom, Cyathus stercoreus, and has free-radical-scavenging activities on the 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical and is thus a potent antioxidant [56]. The fungus is naturally found in the humicolous soil in the Khasi hills of Meghalaya, the upper Shillong area of the IHR [57].

5. Conclusions

The Kinesin protein, Eg5, is a well-established anti-lung cancer target for its role in the assembling of the mitotic spindle directional movement during cell division. Therefore, in the current work, four natural anti-lung cancer inhibitors, viz., (−)-Cochlactone A, Phelligridin C, Sterenin E, and Cyathusal A, have been identified against lung cancer targeting the Eg5 protein by applying the machine-learning and chemoinformatics approaches. All the screened compounds were of mycotic origin and found in the fungi from the IHR. Being of a natural origin, the screened mycotic secondary metabolites may have fewer side effects than other chemotherapeutic drugs used against lung cancer and can be obtained easily by culturing the fungi in the laboratory. They are potential anti-inhibitors as screened by the machine-learning model designed using a lung cancer inhibitor bioassay and have efficient binding potential with the Eg5 protein. Additionally, they have passed all the selected drug-likeness filters and have an excellent presentation of diverse functional groups with pharmacological importance and a stable-binding profile during the MD simulation process, which proves their potential as an anti-cancer drug candidate.

Supplementary Materials

Table S1: Percent growth inhibition data of active and inactive compounds used for ML, Table S2: Training set performance matrix, Table S3: Test set performance matrix, Table S4: Pharmacological properties of screened ligands by machine learning, Table S5: Percent growth inhibition data of active and inactive compounds used for ML, Table S6: Binding stability parameters of the screened compounds during molecular dynamics simulation.

Author Contributions

P.M. designed and wrote the draft manuscript; P.S. performed simulation studies; M.N., M.A.R., S.C., and I.D.B. formulated, critically reviewed, and finalized the entire manuscript; T.J., S.M. (Shalini Mathpal), J.S.-G., R.P., S.M. (Shafi Mahmud), S.A., M.M.G., M.A., A.A.A.A., and M.M.A. helped with the identification of the target and helped in collecting information. All authors have read and agreed to the published version of the manuscript.

Funding

This research work has been carried out in view of the public interest and without receiving any financial support from any source.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are thankful to the Director, G.B. Pant National Institute of Himalayan Environment, Kosi-Katarmal, Almora, Uttarakhand, and Head, Botany Department, S.S.J. University, Almora for providing the necessary facilities to carry out the present study.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

Abbreviations

ALArtificial Intelligence
AMRAtom Molar Refractivity
CAGRCompound Annual Growth Rate
CDKChemistry Development Kit
CSFCorrelation-Based Feature Selection
EREstrogen Receptor
GUIGraphic User Interface
HERElectronic Health Records
IHR Indian Himalayan Region
KIF-11Kinesin Family Member-11
MCF-7Michigan Cancer Foundation-7
MDMolecular Dynamics
MeFSATMedicinal Fungi Secondary Metabolite And Therapeutics
MLMachine Learning
NICNational Cancer Institute
PDBProtein Data Bank
PLIPProtein Legend Interaction Profiler
QEDQuantitative Estimate of Drug-Likeness
QSARQuantitative Structure-Activity Relationship
RGRadius of Gyration
RMSDRoot Mean Square Deviations
ROCReceiver Operating Characteristic Curve
SMILESStandard Data Format (SDF) To Molecular-Input Line-Entry System
TPSATotal Polar Surface Area
USDUnited States Dollar
WEKAWaikato Environment for Knowledge Analysis

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  2. Chand, B.; Kuniyal, J.C.; Chand, R. Ambient Air Quality and Its Sources Surrounding to Hydropower Projects in the Satluj Basin, Northwestern Himalaya, India. MAPAN 2019, 34, 495–510. [Google Scholar] [CrossRef]
  3. Garcia-Saez, I.; Skoufias, D.A. Eg5 Targeting Agents: From new anti-mitotic based inhibitor discovery to cancer therapy and resistance. Biochem. Pharmacol. 2020, 1, 114364. [Google Scholar] [CrossRef] [PubMed]
  4. Bertran, M.T.; Sdelci, S.; Regué, L.; Avruch, J.; Caelles, C.; Roig, J. Nek9 is a Plk1-activated kinase that controls early centrosome separation through Nek6/7 and Eg5. EMBO J. 2011, 30, 2634–2647. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sun, D.; Lu, J.; Ding, K.; Bi, D.; Niu, Z.; Cao, Q.; Zhang, J.; Ding, S. The expression of Eg5 predicts a poor outcome for patients with renal cell carcinoma. Med. Oncol. 2013, 30, 476. [Google Scholar] [CrossRef] [PubMed]
  6. Nand, M.; Maiti, P.; Pant, R.; Kumari, M.; Chandra, S.; Pande, V. Virtual screening of natural compounds as inhibitors of EGFR 696-1022 T790M associated with non-small cell lung cancer. Bioinformation 2016, 12, 311–317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Kim, C.D.; Kim, E.D.; Liu, L.; Buckley, R.S.; Parameswaran, S.; Kim, S.; Wojcik, E.J. Small molecule allosteric uncoupling of microtubule depolymerase activity from motility in human Kinesin-5 during mitotic spindle assembly. Sci. Rep. 2019, 9, 19900. [Google Scholar] [CrossRef]
  8. Li, Z.; Yu, B.; Qi, F.; Li, F. KIF11 Serves as an Independent Prognostic Factor and Therapeutic Target for Patients with Lung Adenocarcinoma. Front. Oncol. 2021, 11, 660218. [Google Scholar] [CrossRef]
  9. Saijo, T.; Ishii, G.; Ochiai, A.; Yoh, K.; Goto, K.; Nagai, K.; Kato, H.; Nishiwaki, Y.; Saijo, N. Eg5 expression is closely correlated with the response of advanced non-small cell lung cancer to antimitotic agents combined with platinum chemotherapy. Lung Cancer 2006, 54, 217–225. [Google Scholar] [CrossRef]
  10. Paier, C.R.K.; Maranhão, S.S.; Carneiro, T.R.; Lima, L.M.; Rocha, D.D.; Santos, R.D.S.; de Farias, K.M.; de Moraes-Filho, M.O.; Pessoa, C. Natural products as new antimitotic compounds for anticancer drug development. Clinics 2018, 73 (Suppl. S1). [Google Scholar] [CrossRef]
  11. Evidente, A.; Kornienko, A.; Cimmino, A.; Andolfi, A.; Lefranc, F.; Mathieu, V.; Kiss, R. Fungal metabolites with anticancer activity. Nat. Prod. Rep. 2014, 31, 617–627. [Google Scholar] [CrossRef] [PubMed]
  12. Issa, N.T.; Stathias, V.; Schürer, S.; Dakshanamurthy, S. Machine and Deep Learning Approaches for Cancer Drug Repurposing. Semin. Cancer Biol. 2021, 68, 132–142. [Google Scholar] [CrossRef] [PubMed]
  13. Global Artificial Intelligence in Drug Discovery Market Size Analysis 2018–2028. Available online: https://Bekryl.com (accessed on 18 December 2021).
  14. Ponzoni, I.; Sebastián-Pérez, V.; Martínez, M.J.; Roca, C.; Pérez, C.D.L.C.; Cravero, F.; Vazquez, G.; Páez, J.A.; Díaz, M.F.; Campillo, N.E. QSAR Classification Models for Predicting the Activity of Inhibitors of Beta-Secretase (BACE1) Associated with Alzheimer’s Disease. Sci. Rep. 2019, 9, 9102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Mendez, D.; Gaulton, A.; Bento, A.P.; Chambers, J.; de Veij, M.; Félix, E.; Magariños, M.P.; Mosquera, J.F.; Mutowo, P.; Nowotka, M.; et al. ChEMBL: Towards direct deposition of bioassay data. Nucleic Acids Res. 2019, 47, 930–940. [Google Scholar] [CrossRef]
  16. Vivek-Ananth, R.P.; Sahoo, A.K.; Kumaravel, K.; Mohanraj, K.; Samal, A. MeFSAT: A curated natural product database specific to secondary metabolites of medicinal fungi. RSC Adv. 2021, 11, 2596–2607. [Google Scholar] [CrossRef]
  17. 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] [Green Version]
  18. 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]
  19. Frank, E.; Hall, M.; Trigg, L.; Holmes, G.; Witten, I.H. Data mining in bioinformatics using Weka. Bioinformatics 2004, 20, 2479–2481. [Google Scholar] [CrossRef] [Green Version]
  20. Puzyn, T.; Leszczynski, J.; Cronin, M.T. (Eds.) Recent Advances in QSAR Studies: Methods and Applications; Springer: Berlin/Heidelberg, Germany, 2010; Volume 8. [Google Scholar] [CrossRef]
  21. Manallack, D.T.; Prankerd, R.J.; Yuriev, E.; Oprea, T.I.; Chalmers, D.K. The Significance of Acid/Base Properties in Drug Discovery. Chem. Soc. Rev. 2013, 42, 485–496. [Google Scholar] [CrossRef] [Green Version]
  22. Burden, F.R. A Chemically Intuitive Molecular Index Based on the Eigenvalues of a Modified Adjacency Matrix. Quant. Struct.-Act. Relatsh. 1997, 16, 309–314. [Google Scholar] [CrossRef]
  23. Wildman, S.A.; Crippen, G.M. Prediction of Physicochemical Parameters by Atomic Contributions. J. Chem. Inf. Comput. Sci. 1999, 39, 868–873. [Google Scholar] [CrossRef]
  24. Hall, L.H.; Mohney, B.; Kier, L.B. The Electrotopological State: An Atom Index for QSAR. Quant. Struct. Act. Relatsh. 1991, 10, 43–51. [Google Scholar] [CrossRef]
  25. Wang, R.; Fu, Y.; Lai, L. A New Atom-Additive Method for Calculating Partition Coefficients. J. Chem. Inf. Comput. Sci. 1997, 37, 615–621. [Google Scholar] [CrossRef]
  26. He, S.; Ye, T.; Wang, R.; Zhang, C.; Zhang, X.; Sun, G.; Sun, X. An In Silico Model for Predicting Drug-Induced Hepatotoxicity. Int. J. Mol. Sci. 2019, 20, 1897. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Taser, P.Y. Application of Bagging and Boosting Approaches Using Decision Tree-Based Algorithms in Diabetes Risk Prediction. Proceedings 2021, 74, 6. [Google Scholar] [CrossRef]
  28. Lučić, B.; Batista, J.; Bojović, V.; Lovrić, M.; Kržić, A.S.; Bešlo, D.; Nadramija, D.; Vikić-Topić, D. Estimation of Random Accuracy and its Use in Validation of Predictive Quality of Classification Models within Predictive Challenges. Croat. Chem. Acta 2019, 92, 379–391. [Google Scholar] [CrossRef] [Green Version]
  29. Geete, A.; Damre, M.; Kokkula, A. Drug Likeness Tool (DruLiTo) Chemistry Development Kit (CDK), Department of Pharmacoinformatics NIPER, Mohali. Available online: http://www.niper.gov.in (accessed on 22 July 2021).
  30. Steinbeck, C.; Han, Y.Q.; Kuhn, S.; Horlacher, O.; Luttmann, E.; Willighagen, E.L. The Chemistry Development Kit (CDK): An open-source Java library for chemo- and bioinformatics. J. Chem. Inf. Comput. Sci. 2003, 43, 493–500. [Google Scholar] [CrossRef] [Green Version]
  31. Nand, M.; Maiti, P.; Joshi, T.; Chandra, S.; Pande, V.; Kuniyal, J.C.; Ramakrishnan, M.A. Virtual screening of anti-HIV1 compounds against SARS-CoV-2: Machine learning modeling, chemoinformatics and molecular dynamics simulation based analysis. Sci. Rep. 2021, 10, 20397. [Google Scholar] [CrossRef]
  32. Talapatra, S.K.; Anthony, N.G.H.; Mackay, S.P.; Kozielski, F. Mitotic Kinesin Eg5 Overcomes Inhibition to the Phase I/II Clinical Candidate SB743921 by an Allosteric Resistance Mechanism. J. Med. Chem. 2013, 56, 6317–6329. [Google Scholar] [CrossRef]
  33. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Meng, E.C.; Couch, G.S.; Croll, T.I.; Morris, J.H.; Ferrin, T.E. UCSF ChimeraX: Structure visualization for researchers, educators, and developers. Protein Sci. Publ. Protein Soc. 2021, 30, 70–82. [Google Scholar] [CrossRef]
  34. Adasme, M.F.; Linnemann, K.L.; Bolz, S.N.; Kaiser, F.; Salentin, S.; Haupt, V.J.; Schroeder, M. PLIP 2021: Expanding the scope of the protein–ligand interaction profiler to DNA and RNA. Nucleic Acids Res. 2021, 49, 530–534. [Google Scholar] [CrossRef] [PubMed]
  35. R Development Core Team. A language and environment for statistical computing: Reference index. Vienna: R Foundation for Statistical Computing. 2010. Available online: http://www.polsci.wvu.edu (accessed on 26 July 2021).
  36. Cao, Y.; Charisi, A.; Cheng, L.-C.; Jiang, T.; Girke, T. ChemmineR: Acompound mining framework for R. Bioinformatics 2008, 24, 1733–1734. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Maiti, P.; Nand, M.; Joshi, T.; Ramakrishnan, M.A.; Chandra, S. Identification of luteolin -7-glucoside and epicatechin gallate from Vernonia cinerea, as novel EGFR L858R kinase inhibitors against lung cancer: Docking and simulation-based study. J. Biomol. Struct. Dyn. 2020, 5048–5057. [Google Scholar] [CrossRef] [PubMed]
  38. Lagunin, A.A.; Dubovskaja, V.I.; Rudik, A.; Pogodin, P.V.; Druzhilovskiy, D.; Gloriozova, T.A.; Filimonov, D.; Sastry, N.G.; Poroikov, V.V. CLC-Pred: A freely available web-service for in silico prediction of human cell line cytotoxicity for drug-like compounds. PLoS ONE 2018, 25, e0191838. [Google Scholar] [CrossRef] [Green Version]
  39. Filimonov, D.; Lagunin, A.A.; Gloriozova, T.A.; Rudik, A.; Druzhilovskii, D.S.; Pogodin, P.V.; Poroikov, V.V. Prediction of the biological activity spectra of organic compounds using the PASS online web resource. Chem. Heterocycl. Compd. 2014, 50, 444–457. Available online: http://www.way2drug.com/PASSonline (accessed on 8 June 2021). [CrossRef]
  40. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef] [Green Version]
  41. Shahbaaz, M.; Nkaule, A.; Christoffels, A. Designing novel possible kinase inhibitor derivatives as therapeutics against Mycobacterium tuberculosis: An in silico study. Sci. Rep. 2019, 9, 4405. [Google Scholar] [CrossRef]
  42. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  43. Jukic, S.; Saracevic, M.; Subasi, A.; Kevric, J. Comparison of Ensemble Machine Learning Methods for Automated Classification of Focal and Non-Focal Epileptic EEG Signals. Mathematics 2020, 8, 1481. [Google Scholar] [CrossRef]
  44. Patrick, G.L. An Introduction to Medicinal Chemistry; Oxford University Press: New York, NY, USA, 2017; pp. 1–832. Available online: https://global.oup.com (accessed on 26 August 2021).
  45. Connelly, P.R.; Thomson, J.A. Heat capacity changes and hydrophobic interactions in the binding of FK506 and rapamycin to the FK506 binding protein. Proc. Natl. Acad. Sci. USA 1992, 89, 4781–4785. [Google Scholar] [CrossRef] [Green Version]
  46. Kostal, J. Chapter Four - Computational Chemistry in Predictive Toxicology. Adv. Mol. Toxicol. 2016, 10, 139–186. [Google Scholar] [CrossRef]
  47. Ogunwa, T.H.; Laudadio, E.; Galeazzi, R.; Miyanishi, T. Insights into the Molecular Mechanisms of Eg5 Inhibition by (+)-Morelloflavone. Pharmaceuticals 2019, 12, 58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Banyal, A.; Thakur, V.; Thakur, R.; Kumar, P. Endophytic Microbial Diversity: A New Hope for the Production of Novel Anti-tumor and Anti-HIV Agents as Future Therapeutics. Curr. Microbiol. 2021, 78, 1699–1717. [Google Scholar] [CrossRef] [PubMed]
  49. Gauchan, D.P.; Vélëz, H.; Acharya, A.; Östman, J.R.; Lundén, K.; Elfstrand, M.; García-Gil, M.R. Annulohypoxylon sp. Strain MUS1, an endophytic fungus isolated from Taxus wallichiana Zucc., produces taxol and other bioactive metabolites. 3 Biotech 2021, 11, 152. [Google Scholar] [CrossRef]
  50. Peng, X.-R.; Lu, S.-Y.; Shao, L.-D.; Zhou, L.; Qiu, M.-H. Structural Elucidation and Biomimetic Synthesis of (±)-Cochlactone A with Anti-Inflammatory Activity. J. Org. Chem. 2018, 83, 5516–5522. [Google Scholar] [CrossRef]
  51. Bijalwan, A.; Bahuguna, K.; Vasishth, A.; Singh, A.; Chaudhary, S.; Dongariyal, A.; Thakur, T.K.; Kaushik, S.; Ansari, M.J.; Alfarraj, S.; et al. Growth performance of Ganoderma lucidum using billet method in Garhwal Himalaya, India. Saudi Biol. Sci. 2021, 28, 2709–2717. [Google Scholar] [CrossRef]
  52. Mo, S.; Wang, S.; Zhou, G.; Yang, Y.; Li, Y.; Chen, A.X.; Shi, J. Phelligridins C-F: Cytotoxic pyrano [4,3-c][2]benzopyran-1,6-dione and furo[3,2-c]pyran-4-one derivatives from the fungus Phellinus igniarius. J. Nat. Prod. 2004, 67, 823–828. [Google Scholar] [CrossRef]
  53. Azeem, U.; Dhingra, G.S.; Shri, R. Pharmacological potential of wood inhabiting fungi of genus Phellinus quél.: An overview. J. Pharmacogn. Phytochem. 2018, 7, 1161–1171. [Google Scholar]
  54. Wang, B.T.; Qi, Q.Y.; Ma, K.; Pei, Y.F.; Han, J.J.; Xu, W.; Li, E.W.; Liu, H.W. Depside α-glucosidase inhibitors from a culture of the mushroom. Stereum hirsutum. Planta Med. 2014, 80, 918–924. [Google Scholar] [CrossRef] [Green Version]
  55. Vabeikhokhei, J.M.; Mangaiha, Z.; Zothanzama, J.; Lalrinawmi, H. Diversity Study of Wood Rotting Fungi from Two different Forests in Mizoram, India. Int. J. Curr. Microbiol. Appl. Sci. 2019, 8, 2775–2785. [Google Scholar] [CrossRef]
  56. Kang, H.-S.; Jun, E.-M.; Park, S.-H.; Heo, S.-J.; Lee, T.-S.; Yoo, I.-D.; Kim, J.-P. Cyathusals A, B, and C, Antioxidants from the Fermented Mushroom Cyathus stercoreus. J. Nat. Prod. 2007, 70, 1043–1045. [Google Scholar] [CrossRef] [PubMed]
  57. Sharma, B.M. Genus Cyathus Haller ex Pers. (Agaricomycetes) from Eastern Himalaya. KAVAKA 2016, 47, 20–26. Available online: http://www.fungiindia.co.in (accessed on 23 October 2021).
Figure 1. Statistical performance of different classifiers used for the development of screening model in the training set.
Figure 1. Statistical performance of different classifiers used for the development of screening model in the training set.
Molecules 27 01639 g001
Figure 2. (a) Drug-likeness filters used for screening compounds. (b,c) Active site prediction of Eg5 protein for molecular docking.
Figure 2. (a) Drug-likeness filters used for screening compounds. (b,c) Active site prediction of Eg5 protein for molecular docking.
Molecules 27 01639 g002
Figure 3. (a) Binding free energy of the screened ligands through molecular docking. (b) Functional group frequency comparison between established inhibitors and screened compounds. (c) Binding insights of the screened ligand, Phelligridin-C, with Eg5 protein. (d) 2D structure of screened compound.
Figure 3. (a) Binding free energy of the screened ligands through molecular docking. (b) Functional group frequency comparison between established inhibitors and screened compounds. (c) Binding insights of the screened ligand, Phelligridin-C, with Eg5 protein. (d) 2D structure of screened compound.
Molecules 27 01639 g003
Figure 4. Interaction profile of screened ligands with Eg5 protein: (a) (-)-Cochlactone A, (b) Phelligridin-C, (c) Sterenin E, and (d) Cyathusal A.
Figure 4. Interaction profile of screened ligands with Eg5 protein: (a) (-)-Cochlactone A, (b) Phelligridin-C, (c) Sterenin E, and (d) Cyathusal A.
Molecules 27 01639 g004
Figure 5. Curves illustrating the behavior of the interactions of screened compounds with EG5 protein in the form of Rg, RMSF, and RMSD during MD simulation.
Figure 5. Curves illustrating the behavior of the interactions of screened compounds with EG5 protein in the form of Rg, RMSF, and RMSD during MD simulation.
Molecules 27 01639 g005
Table 1. Comparison of performance of different classifiers for development of screening model in the training set.
Table 1. Comparison of performance of different classifiers for development of screening model in the training set.
Classifier
Name
Correctly Classified Instances % (Value)Kappa StatisticMean Absolute ErrorRoot Mean Square ErrorMCCROC Area
Random forest97.05880.94010.080.17310.9420.989
J4896.79140.93460.050.1750.9370.964
Decision stump96.79140.93460.060.1750.9370.947
Random tree92.78070.85440.070.26870.8550.928
Bagging (REP tree)96.52410.92920.060.18440.9310.96
Table 2. Pharmacological indices of the screened ligands by Lipinski’s rule, CMC-50 like rule, Veber rule, and Ghose filter.
Table 2. Pharmacological indices of the screened ligands by Lipinski’s rule, CMC-50 like rule, Veber rule, and Ghose filter.
Title *12345678910111213
Pharmacological Indices
MW358358372372330346378378380364386370358
logp4422221122442
Alogp11110−11101330
HBA5566786687767
HBD2233234443433
TPSA84849696102123107107134113124104113
AMR9898105105909111311310610410910887
nRB3344333322665
nAtom52525151383946464039504955
nAcidicGroup0000000000100
RC4433334444222
nRigidB26262525232428282928232221
nAromRing1111113322220
nHB77999111010121011910
SAlerts4455550044342
* Compound CID: 1-138970029, 2-139591442, 3-139590281, 4-139590280, 5-16737098, 6-16737097, 7-12085445, 8-54586497, 9-10339712, 10-10248188, 11-77461063, 12-77461065, and 13-10203477.
Table 3. Interaction profile of screened ligands with Eg5 protein.
Table 3. Interaction profile of screened ligands with Eg5 protein.
Ligand NameHydrophobic InteractionsHydrogen BondOther
ResidueAADistanceResidueAADistanceResidueAADistance
(−)-Cochlactone-A79BILE3.85286AGLY1.85Salt Bridges
---286AGLY2.16138BARG4.45
131BPRO3.53297AARG3.13141BHIS5.14
---------
285AALA3.42------
---------
Phelligridin-C79BILE3.9283BARG3.27π–Cation Interactions
125BTYR3.72142BGLN2.5583BARG4.97
---286AGLY1.92
131BPRO3.83290AGLN2.5Salt Bridges
285AALA3.29297AARG2.2138BARG4.1
Sterenin-E79BILE3.4383BARG3.09π–Cation Interactions
82BTYR3.52141BHIS2.6983BARG5.19
---142BGLN2.72Salt Bridges
---142BGLN2.66138BARG4.74
293ALEU3.38287AASN2.55141BHIS5.08
---290AGLN2.66---
---297AARG3.36---
Cyathusal-A79BILE3.7383BARG3.15π–Cation Interactions
---138BARG2.3783BARG5.11
285AALA3.92142BGLN2.51---
---290AGLN2.5---
---297AARG2.45---
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maiti, P.; Sharma, P.; Nand, M.; Bhatt, I.D.; Ramakrishnan, M.A.; Mathpal, S.; Joshi, T.; Pant, R.; Mahmud, S.; Simal-Gandara, J.; et al. Integrated Machine Learning and Chemoinformatics-Based Screening of Mycotic Compounds against Kinesin Spindle ProteinEg5 for Lung Cancer Therapy. Molecules 2022, 27, 1639. https://doi.org/10.3390/molecules27051639

AMA Style

Maiti P, Sharma P, Nand M, Bhatt ID, Ramakrishnan MA, Mathpal S, Joshi T, Pant R, Mahmud S, Simal-Gandara J, et al. Integrated Machine Learning and Chemoinformatics-Based Screening of Mycotic Compounds against Kinesin Spindle ProteinEg5 for Lung Cancer Therapy. Molecules. 2022; 27(5):1639. https://doi.org/10.3390/molecules27051639

Chicago/Turabian Style

Maiti, Priyanka, Priyanka Sharma, Mahesha Nand, Indra D. Bhatt, Muthannan Andavar Ramakrishnan, Shalini Mathpal, Tushar Joshi, Ragini Pant, Shafi Mahmud, Jesus Simal-Gandara, and et al. 2022. "Integrated Machine Learning and Chemoinformatics-Based Screening of Mycotic Compounds against Kinesin Spindle ProteinEg5 for Lung Cancer Therapy" Molecules 27, no. 5: 1639. https://doi.org/10.3390/molecules27051639

APA Style

Maiti, P., Sharma, P., Nand, M., Bhatt, I. D., Ramakrishnan, M. A., Mathpal, S., Joshi, T., Pant, R., Mahmud, S., Simal-Gandara, J., Alshehri, S., Ghoneim, M. M., Alruwaily, M., Awadh, A. A. A., Alshahrani, M. M., & Chandra, S. (2022). Integrated Machine Learning and Chemoinformatics-Based Screening of Mycotic Compounds against Kinesin Spindle ProteinEg5 for Lung Cancer Therapy. Molecules, 27(5), 1639. https://doi.org/10.3390/molecules27051639

Article Metrics

Back to TopTop