Combinatorial support vector machines approach for virtual screening of selective multi-target serotonin reuptake inhibitors from large compound libraries

A multiple-step screening protocol to identify norepinephrine and dopamine reuptake inhibitors for depression

Depression severely impairs the health of people all over the world. Cognitive dysfunction due to depression has resulted in a severe economic burden to family and society induced by the reduction of social functioning of patients. Norepinephrine-dopamine reuptake inhibitors (NDRIs) targeted with the human norepinephrine transporter (hNET) and distributed with the human dopamine transporter (hDAT) simultaneously treat depression and improve cognitive function, and they effectively prevent sexual dysfunction and other side effects. Because many patients continue to poorly respond to NDRIs, it is urgent to discover novel NDRI antidepressants that do not interfere with cognitive function. The aim of this work was to selectively identify novel NDRI candidates acting against hNET and hDAT from large compound libraries by a comprehensive strategy integrating support vector machine (SVM) models, ADMET, molecular docking, in vitro binding assays, molecular dynamics simulation, and binding energy calculation. First, 6522 compounds that do not inhibit the human serotonin transporter (hSERT) were obtained by SVM models of hNET, hDAT, and non-target hSERT with similarity analyses from compound libraries. ADMET and molecular docking were then used to identify compounds that could potently bind to the hNET and hDAT with satisfactory ADMET, and 4 compounds were successfully identified. According to their docking scores and ADMET information, 3719810 was advanced for profiling by in vitro assays as a novel NDRI lead compound due to its strongest druggability and balancing activities. Encouragingly, 3719810 performed comparative activities on two targets, with Ki values of 7.32 μM for hNET and 5.23 μM for hDAT. To obtain candidates with additional activities and balance the activities of 2 targets, 5 analogs were optimized, and 2 novel scaffold compounds were successively designed. By assessment of molecular docking, molecular dynamics simulations, and binding energy calculations, 5 compounds were validated as NDRI candidates with high activities, and 4 of them performed acceptable balancing activities acting on hNET and hDAT. This work supplied promising novel NDRIs for treatment of depression with cognitive dysfunction or other related neurodegenerative disorders, and also provided a strategy for highly efficient and cost-effective identification of inhibitors for dual targets with homologous non-targets.

Evolution of Support Vector Machine and Regression Modeling in Chemoinformatics and Drug Discovery

The support vector machine (SVM) algorithm is one of the most widely used machine learning (ML) methods for predicting active compounds and molecular properties. In chemoinformatics and drug discovery, SVM has been a state-of-the-art ML approach for more than a decade. A unique attribute of SVM is that it operates in feature spaces of increasing dimensionality. Hence, SVM conceptually departs from the paradigm of low dimensionality that applies to many other methods for chemical space navigation. The SVM approach is applicable to compound classification, and ranking, multi-class predictions, and –in algorithmically modified form– regression modeling. In the emerging era of deep learning (DL), SVM retains its relevance as one of the premier ML methods in chemoinformatics, for reasons discussed herein. We describe the SVM methodology including strengths and weaknesses and discuss selected applications that have contributed to the evolution of SVM as a premier approach for compound classification, property predictions, and virtual compound screening.

Exploring the binding mechanism of positive allosteric modulators in human metabotropic glutamate receptor 2 using molecular dynamics simulations

Positive allosteric modulators (PAMs) of human metabotropic glutamate receptor 2 (hmGlu₂) are well-known in the treatment of psychiatric disorders for their higher selectivity and lower tolerance risk. A variety of PAMs have been reported over the last decade and two compounds were in Phase II clinical trials for schizophrenia and anxiety. These trials were discontinued on account of the unsatisfactory therapeutic efficacy, but PAMs were explored as novel treatments for addiction and epilepsy. Thus, it is still important to explore novel hmGlu₂ PAMs in the near future. Nowadays, the challenges in optimizing drug potency and improving scaffold diversity for PAMs are the noncomprehensive character analyses of multiple scaffolds; the exploration of the binding modes of PAMs in the allosteric binding site have been proposed to reduce this difficulty. However, there has been no comprehensive research about the binding profiles of PAMs in the hmGlu₂ receptor. To address this issue, this work explores the binding characters of eight PAMs representing five chemical series by multiple computational methods. As a result, the shared binding modes of the eight studied PAMs interacting with 15 residues in the allosteric binding site were defined. In addition, the reduced hydrophobicity with low electronegativity of R₁, increased hydrophobicity with low negative electron density of R₂ and the electronegativity of the linker were identified as indicators that regulate the affinity of PAMs. This finding agrees well with the physicochemical properties of reported multiple series PAMs. This comprehensive work sheds additional light on the binding mechanism and physicochemical regularity underlining PAMs affinity and could be further utilized as a structural and energetic blueprint for discovering and assessing novel PAMs for hmGlu₂.

Discovery of novel dual VEGFR2 and Src inhibitors using a multistep virtual screening approach

AIM

Simultaneous inhibition of VEGFR2 and Src may enhance the efficacy of VEGFR2-targeted cancer therapeutics. Hence, development of dual inhibitors on VEGFR2 and Src can be a useful strategy for such treatments.

MATERIALS & METHODS

A multistep virtual screening protocol, comprising ligand-based support vector machines method, drug-likeness rules filter and structure-based molecular docking, was developed and employed to identify dual inhibitors of VEGFR2 and Src from a large commercial chemical library. Kinase inhibitory assays and cell viability assays were then used for experimental validation.

RESULTS

A set of compounds belonging to six different molecular scaffolds was identified and sent for biological evaluation. Compound 3c belonging to the 2-amino-3-cyanopyridine scaffold exhibited good antiproliferative effect and dual-target activities against VEGFR2 and Src.

CONCLUSION

This study demonstrated the ability of the multistep virtual screening approach to identify novel multitarget agents.

A lazy learning-based QSAR classification study for screening potential histone deacetylase 8 (HDAC8) inhibitors

Histone deacetylases 8 (HDAC8) is an enzyme repressing the transcription of various genes including tumour suppressor gene and has already become a target of human cancer treatment. In an effort to facilitate the discovery of HDAC8 inhibitors, two quantitative structure-activity relationship (QSAR) classification models were developed using K nearest neighbours (KNN) and neighbourhood classifier (NEC). Molecular descriptors were calculated for the data set and database compounds using ADRIANA.Code of Molecular Networks. Principal components analysis (PCA) was used to select the descriptors. The developed models were validated by leave-one-out cross validation (LOO CV). The performances of the developed models were evaluated with an external test set. Highly predictive models were used for database virtual screening. Furthermore, hit compounds were subsequently subject to molecular docking. Five hits were obtained based on consensus scoring function and binding affinity as potential HDAC8 inhibitors. Finally, HDAC8 structures in complex with five hits were also subjected to 5 ns molecular dynamics (MD) simulations to evaluate the complex structure stability. To the best of our knowledge, the NEC classification model used in this study is the first application of NEC to virtual screening for drug discovery.

In silico de novo design of novel NNRTIs: a bio-molecular modelling approach

Six novel NNRTIs (DABO) with high efficacy are designed by assessing the interaction potential and structural requirements using chemometric analyses (SVM, BPNN and MLR) on structural descriptors.

Design of novel leads: ligand based computational modeling studies on non-nucleoside reverse transcriptase inhibitors (NNRTIs) of HIV-1

Researchers are on the constant lookout for new antiviral agents for the treatment of AIDS. In the present work, ligand based modeling studies are performed on analogues of substituted phenyl-thio-thymines, which act as non-nucleoside reverse transcriptase inhibitors (NNRTIs) and novel leads are extracted. Using alignment-dependent descriptors, based on group center overlap (SALL, HDALL, HAALL and RALL), an alignment-independent descriptor (S log P), a topological descriptor (Balaban index (J)) and a 3D descriptor dipole moment (μ) and shape based descriptors (Kappa 2 index ((2)κ)), a correlation is derived with inhibitory activity. Linear and non-linear techniques have been used to achieve the goal. Support Vector Machine (SVM, R = 0.929, R(2) = 0.863) and Back Propagation Neural Network (BPNN, R = 0.928, R(2) = 0.861) methods yielded near similar results and outperformed Multiple Linear Regression (MLR, R = 0.915, R(2) = 0.837). The predictive ability of the models are cross-validated using a test dataset (SVM: R = 0.846, R(2) = 0.716, BPNN: R = 0.841, R(2) = 0.707 and MLR: R = 0.833, R(2) = 0.694). It is concluded that the hydrophobicity (S log P) and the polarity (μ) of a ligand and the presence of hydrogen donor (HDALL) moieties are the deciding factors in improving antiviral activity and pharmaco-therapeutic properties. Based on the above findings, a virtual dataset is created to extract probable leads with reasonable antiviral activity as well as better pharmacophoric properties.

A two-step target binding and selectivity support vector machines approach for virtual screening of dopamine receptor subtype-selective ligands

Target selective drugs, such as dopamine receptor (DR) subtype selective ligands, are developed for enhanced therapeutics and reduced side effects. In silico methods have been explored for searching DR selective ligands, but encountered difficulties associated with high subtype similarity and ligand structural diversity. Machine learning methods have shown promising potential in searching target selective compounds. Their target selective capability can be further enhanced. In this work, we introduced a new two-step support vector machines target-binding and selectivity screening method for searching DR subtype-selective ligands, which was tested together with three previously-used machine learning methods for searching D1, D2, D3 and D4 selective ligands. It correctly identified 50.6%–88.0% of the 21–408 subtype selective and 71.7%–81.0% of the 39–147 multi-subtype ligands. Its subtype selective ligand identification rates are significantly better than, and its multi-subtype ligand identification rates are comparable to the best rates of the previously used methods. Our method produced low false-hit rates in screening 13.56 M PubChem, 168,016 MDDR and 657,736 ChEMBLdb compounds. Molecular features important for subtype selectivity were extracted by using the recursive feature elimination feature selection method. These features are consistent with literature-reported features. Our method showed similar performance in searching estrogen receptor subtype selective ligands. Our study demonstrated the usefulness of the two-step target binding and selectivity screening method in searching subtype selective ligands from large compound libraries.