PLS and shape-based similarity analysis of maleimides--GSK-3 inhibitors

New Insights on the Activity and Selectivity of MAO-B Inhibitors through In Silico Methods

To facilitate the identification of novel MAO-B inhibitors, we elaborated a consolidated computational approach, including a pharmacophoric atom-based 3D quantitative structure-activity relationship (QSAR) model, activity cliffs, fingerprint, and molecular docking analysis on a dataset of 126 molecules. An AAHR.2 hypothesis with two hydrogen bond acceptors (A), one hydrophobic (H), and one aromatic ring (R) supplied a statistically significant 3D QSAR model reflected by the parameters: R² = 0.900 (training set); Q² = 0.774 and Pearson's R = 0.884 (test set), stability s = 0.736. Hydrophobic and electron-withdrawing fields portrayed the relationships between structural characteristics and inhibitory activity. The quinolin-2-one scaffold has a key role in selectivity towards MAO-B with an AUC of 0.962, as retrieved by ECFP4 analysis. Two activity cliffs showing meaningful potency variation in the MAO-B chemical space were observed. The docking study revealed interactions with crucial residues TYR:435, TYR:326, CYS:172, and GLN:206 responsible for MAO-B activity. Molecular docking is in consensus with and complementary to pharmacophoric 3D QSAR, ECFP4, and MM-GBSA analysis. The computational scenario provided here will assist chemists in quickly designing and predicting new potent and selective candidates as MAO-B inhibitors for MAO-B-driven diseases. This approach can also be used to identify MAO-B inhibitors from other libraries or screen top molecules for other targets involved in suitable diseases.

Natural Compounds as DPP-4 Inhibitors: 3D-Similarity Search, ADME Toxicity, and Molecular Docking Approaches

Type 2 diabetes mellitus is one of the most common diseases of the 21st century, caused by a sedentary lifestyle, poor diet, high blood pressure, family history, and obesity. To date, there are no known complete cures for type 2 diabetes. To identify bioactive natural products (NPs) to manage type 2 diabetes, the NPs from the ZINC15 database (ZINC-NPs DB) were screened using a 3D shape similarity search, molecular docking approaches, and ADMETox approaches. Frequently, in silico studies result in asymmetric structures as “hit” molecules. Therefore, the asymmetrical FDA-approved diabetes drugs linagliptin (8-[(3R)-3-aminopiperidin-1-yl]-7-but-2-ynyl-3-methyl-1-[(4-methylquinazolin-2-yl)methyl]purine-2,6-dione), sitagliptin ((3R)-3-amino-1-[3-(trifluoromethyl)-6,8-dihydro-5H-[1,2,4]triazolo[4,3-a]pyrazin-7-yl]-4-(2,4,5-trifluorophenyl)butan-1-one), and alogliptin (2-[[6-[(3R)-3-aminopiperidin-1-yl]-3-methyl-2,4-dioxopyrimidin-1-yl]methyl]benzonitrile) were used as queries to virtually screen the ZINC-NPs DB and detect novel potential dipeptidyl peptidase-4 (DPP-4) inhibitors. The most promising NPs, characterized by the best sets of similarity and ADMETox features, were used during the molecular docking stage. The results highlight that 11 asymmetrical NPs out of 224,205 NPs are potential DPP-4 candidates from natural sources and deserve consideration for further in vitro/in vivo tests.

Small Molecules of Natural Origin as Potential Anti-HIV Agents: A Computational Approach

The human immunodeficiency virus type 1 (HIV-1), one of the leading causes of infectious death globally, generates severe damages to people's immune systems and makes them susceptible to serious diseases. To date, there are no drugs that completely remove HIV from the body. This paper focuses on screening 224,205 natural compounds of ZINC15 NPs subset to identify those with bioactivity similar to non-nucleoside reverse transcriptase inhibitors (NNRTIs) as promising candidates to treat HIV-1. To reach the goal, an in silico approach involving 3D-similarity search, ADMETox, HIV protein-inhibitor prediction, docking, and MM-GBSA free-binding energies was trained. The FDA-approved HIV drugs, efavirenz, etravirine, rilpivirine, and doravirine, were used as queries. The prioritized compounds were subjected to ADMETox, docking, and MM-GBSA studies against HIV-1 reverse transcriptase (RT). Lys101, Tyr181, Tyr188, Trp229, and Tyr318 residues and free-binding energies have proved that ligands can stably bind to HIV-1 RT. Three natural products (ZINC37538901, ZINC38321654, and ZINC67912677) containing oxan and oxolan rings with hydroxyl substituents and one (ZINC2103242) having 3,6,7,8-tetrahydro-2H-pyrido[1,2-a]pyrazine-1,4-dione core exhibited comparable profiles to etravirine and doravirine, with ZINC2103242 being the most promising anti-HIV candidate in terms of drug metabolism and safety profile. These findings may open new avenues to guide the rational design of novel HIV-1 NNRTIs.

Comprehensive investigation of selectivity landscape of glycogen synthase kinase-3 inhibitors

Interaction signatures of drug candidates are characteristic to off-target (neutral) and antitarget (negative) effects, inferring reduced efficiency, side-effects and high attrition rate. Today's retroactive scaled-down virtual screening (VS) experiments relying on benchmarking datasets are extensively involved to assess ligand enrichment in the real-world problem. In recent years, unbiased benchmarking sets turned into a tremendous need to assist virtual screening methodologies for emerging drug targets. To date, the benchmarking datasets are quite limited, whereas glycogen synthase kinase-3 (GSK-3) is not included into directories of benchmarking datasets such as DUD-e, MUV, etc. Herein we introduced our in-house algorithm to build an unbiased benchmarking dataset, including highly selective, moderately selective and nonselective inhibitors for a significant therapeutic target - GSK-3, suitable for both ligand-based and structure-based VS approaches. These datasets are unbiased in terms of physico-chemical properties and topological descriptors, as resulted from mean(ROC-AUC) leave-one-out cross-validation (LOO CV). and additional 2 D similarity search. Moreover, we investigated the gradual selectivity dataset by application of multiple 2 D similarity coefficients and distances, 3 D similarity and docking. Besides the resulted links between the enrichment of selective GSK-3 inhibitors and their chemical structures, a database of compounds and their 3 D similarity signatures including cut-off thresholds for enhanced selectivity was generated. 2 D similarity space analysis revealed that selectivity problem cannot be evaluated appropriately with 2 D similarity searching alone. The current analysis provided useful, comprehensive insights, which may facilitate the knowledge-based identification of novel selective GSK-3 inhibitors.Communicated by Ramaswamy H. Sarma.

QSAR studies in the discovery of novel type-II diabetic therapies

INTRODUCTION

Type-II diabetes mellitus (T2DM) is a complex chronic disease that represents a major therapeutic challenge. Despite extensive efforts in T2DM drug development, therapies remain unsatisfactory. Currently, there are many novel and important antidiabetic drug targets under investigation by many research groups worldwide. One of the main challenges to develop effective orally active hypoglycemic agents is off-target effects. Computational tools have impacted drug discovery at many levels. One of the earliest methods is quantitative structure-activity relationship (QSAR) studies. QSAR strategies help medicinal chemists understand the relationship between hypoglycemic activity and molecular properties. Hence, QSAR may hold promise in guiding the synthesis of specifically designed novel ligands that demonstrate high potency and target selectivity.

AREAS COVERED

This review aims to provide an overview of the QSAR strategies used to model antidiabetic agents. In particular, this review focuses on drug targets that raised recent scientific interest and/or led to successful antidiabetic agents in the market. Special emphasis has been made on studies that led to the identification of novel antidiabetic scaffolds.

EXPERT OPINION

Computer-aided molecular design and discovery techniques like QSAR have a great potential in designing leads against complex diseases such as T2DM. Combined with other in silico techniques, QSAR can provide more useful and rational insights to facilitate the discovery of novel compounds. However, since T2DM is a complex disease that includes several faulty biological targets, multi-target QSAR studies are recommended in the future to achieve efficient antidiabetic therapies.

Application of Shape Similarity in Pose Selection and Virtual Screening in CSARdock2014 Exercise

To evaluate the applicability of shape similarity in docking-based pose selection and virtual screening, we participated in the CSARdock2014 benchmark exercise for identifying the correct docking pose of inhibitors targeting factor XA, spleen tyrosine kinase, and tRNA methyltransferase. This exercise provides a valuable opportunity for researchers to test their docking programs, methods, and protocols in a blind testing environment. In the CSARdock2014 benchmark exercise, we have implemented an approach that uses ligand 3D shape similarity to facilitate docking-based pose selection and virtual screening. We showed here that ligand 3D shape similarity between bound poses could be used to identify the native-like pose from an ensemble of docking-generated poses. Our method correctly identified the native pose as the top-ranking pose for 73% of test cases in a blind testing environment. Moreover, the pose selection results also revealed an excellent correlation between ligand 3D shape similarity scores and RMSD to X-ray crystal structure ligand. In the virtual screening exercise, the average RMSD for our pose prediction was found to be 1.02 Å, and it was one of the top performances achieved in CSARdock2014 benchmark exercise. Furthermore, the inclusion of shape similarity improved virtual screening performance of docking-based scoring and ranking. The coefficient of determination (r(2)) between experimental activities and docking scores for 276 spleen tyrosine kinase inhibitors was found to be 0.365 but reached 0.614 when the ligand 3D shape similarity was included.