Discovery of DPP IV inhibitors by pharmacophore modeling and QSAR analysis followed by in silico screening

Fragment-based QSAR study to explore the structural requirements of DPP-4 inhibitors: a stepping stone towards better type 2 diabetes mellitus management

Dipeptidyl peptidase-4 (DPP-4) inhibitors belong to a prominent group of pharmaceutical agents that are used in the governance of type 2 diabetes mellitus (T2DM). They exert their antidiabetic effects by inhibiting the incretin hormones like glucagon-like peptide-1 and glucose-dependent insulinotropic polypeptide which, play a pivotal role in the regulation of blood glucose homoeostasis in our body. DPP-4 inhibitors have emerged as an important class of oral antidiabetic drugs for the treatment of T2DM. Surprisingly, only a few 2D-QSAR studies have been reported on DPP-4 inhibitors. Here, fragment-based QSAR (Laplacian-modified Bayesian modelling and Recursive partitioning (RP) approaches have been utilized on a dataset of 108 DPP-4 inhibitors to achieve a deeper understanding of the association among their molecular structures. The Bayesian analysis demonstrated satisfactory ROC values for the training as well as the test sets. Meanwhile, the RP analysis resulted in decision tree 3 with 2 leaves (Tree 3: 2 leaves). This present study is an effort to get an insight into the pivotal fragments modulating DPP-4 inhibition.

Ligand-based pharmacophore modeling and machine learning for the discovery of potent aurora A kinase inhibitory leads of novel chemotypes

Aurora-A (AURKA) is serine/threonine protein kinase involved in the regulation of numerous processes of cell division. Numerous studies have demonstrated strong association between AURKA and cancer. AURKA is overexpressed in many cancers, such as colon, breast and prostate cancers. Consequently, AURKA has emerged as promising target for therapeutic intervention in cancer management. Herein, we describe a computational workflow for the discovery of novel anti-AURKA inhibitory leads starting with ligand-based assessment of the pharmacophoric space of six diverse sets of inhibitors. Subsequently, machine learning/QSAR modeling was coupled with genetic function algorithm to search for the best possible combination of machine learner, ligand-based pharmacophore(s) and molecular descriptors capable of explaining variation in anti-AURKA bioactivities within a collected list of inhibitors. Two learners succeeded in achieving acceptable structure/activity correlations, namely, random forests and extreme gradient boosting (XGBoost). Three pharmacophores emerged in the successful ML models. These were then used as 3D search queries to mine the National Cancer Institute database for novel anti-AURKA leads. Top-ranking 38 hits were assessed in vitro for their anti-AURKA bioactivities. Among them, three compounds exhibited promising dose-response curves, demonstrating experimental IC50 values ranging from sub-micromolar to low micromolar values. Remarkably, two of these compounds are of novel chemotypes.

Structure‐Based Approaches for the Prediction of Alzheimer's Disease Inhibitors: Comparative Interactions Analysis, Pharmacophore Modeling and Molecular Dynamics Simulations

Due to its significant role in neurodegeneration, Cyclin‐dependent kinase 5 (CDK5) has emerged as a potential target for addressing neuropathological disorders, including Alzheimer's disease (AD). The application of CDK5 inhibitors has demonstrated promise in the treatment of AD. This prompted us to model this interesting target using a computational workflow named Docking‐based Comparative Intermolecular Contacts Analysis (dbCICA). Approaches including 3D‐QSAR, genetic algorithm, and pharmacophore modeling were employed to discover new CDK inhibitors.

Discovery of new Cdc2-like kinase 4 (CLK4) inhibitors via pharmacophore exploration combined with flexible docking-based ligand/receptor contact fingerprints and machine learning

Cdc2-like kinase 4 (CLK4) inhibitors are of potential therapeutic value in many diseases particularly cancer. In this study, we combined extensive ligand-based pharmacophore exploration, ligand–receptor contact fingerprints generated by flexible docking, physicochemical descriptors and machine learning-quantitative structure–activity relationship (ML-QSAR) analysis to investigate the pharmacophoric/binding requirements for potent CLK4 antagonists. Several ML methods were attempted to tie these properties with anti-CLK4 bioactivities including multiple linear regression (MLR), random forests (RF), extreme gradient boosting (XGBoost), probabilistic neural network (PNN), and support vector regression (SVR). A genetic function algorithm (GFA) was combined with each method for feature selection. Eventually, GFA-SVR was found to produce the best self-consistent and predictive model. The model selected three pharmacophores, three ligand–receptor contacts and two physicochemical descriptors. The GFA-SVR model and associated pharmacophore models were used to screen the National Cancer Institute (NCI) structural database for novel CLK4 antagonists. Three potent hits were identified with the best one showing an anti-CLK4 IC₅₀ value of 57 nM.

Ligand-based pharmacophores, ligand–receptor contact fingerprints, physicochemical descriptors and machine learning were combined to probe binding of potent CLK4 antagonists. GFA-SVR gave the best model. Virtual screening identified 3 nanomolar hits.

QSAR-guided pharmacophoric modeling reveals important structural requirements for Polo kinase 1 (Plk1) inhibitors

Targeting Polo-like kinase 1 (Plk1) by molecular inhibitors is being a promising approach for tumor therapy. Nevertheless, insufficient methodical analyses have been done to characterize the interactions inside the Plk1 binding pocket. In this study, an extensive combined ligand and structure-based drug design workflow was conducted to data-mine the structural requirements for Plk1 inhibition. Consequently, the binding modes of 368 previously known Plk1 inhibitors were investigated by pharmacophore generation technique. The resulted pharmacophores were engaged in the context of Genetic function algorithm (GFA) and Multiple linear regression (MLR) analyses to search for a prognostic QSAR model. The most successful QSAR model was with statistical criteria of (r2277 = 0.76, r2adj = 0.76, r2pred = 0.75, Q2 = 0.73). Our QSAR-selected pharmacophores were validated by Receiver Operating Characteristic (ROC) curve analysis. Later on, the best QSAR model and its associated pharmacophoric hypotheses (HypoB-T4-5, HypoI-T2-7, HypoD-T4-3, and HypoC-T3-3) were used to identify new Plk1 inhibitory hits retrieved from the National Cancer Institute (NCI) database. The most potent hits exhibited experimental anti-Plk1 IC50 of 1.49, 3.79. 5.26 and 6.35 μM. Noticeably, our hits, were found to interact with the Plk1 kinase domain through some important amino acid residues namely, Cys67, Lys82, Cys133, Phe183, and Asp194.

Application of machine learning methods for the development of antidiabetic drugs

Diabetes is a chronic non-communicable disease caused by several different routes, which has attracted increasing attention. In order to speed up the development of new selective drugs, machine learning (ML) technology has been applied in the process of diabetes drug development, which opens up a new blueprint for drug design. This review provides a comprehensive portrayal of the application of ML in antidiabetic drug use.

Ligand Based Pharmacophore Modeling Followed by Biological Screening Lead to Discovery of Novel PDK1 Inhibitors as Anticancer Agents

BACKGROUND

Phosphoinositide-Dependent Kinase-1 (PDK1) is a serine/threonine kinase, which belongs to AGC kinase family required by cancer cells.

METHODS

Pharmacophoric space of 86 PDK1 inhibitors using six diverse sets of inhibitors was explored to identify high-quality pharmacophores. The best combination of pharmacophoric models and physicochemical descriptors was selected by genetic algorithm-based QSAR analysis that can elucidate the variation of bioactivity within the training inhibitors. Two successful orthogonal pharmacophores emerged in the optimum QSAR equation (r2 69 = 0.90, r2 LOO= 0.86, F= 51.92, r2 PRESS against 17 test inhibitors = 0.79). Receiver Operating Characteristic (ROC) curve analyses were used to estimate the QSAR-selected pharmacophores.

RESULTS

5 out of 11 compounds tested had shown potential intracellular PDK1 inhibition with the highest inhibition percent for compounds 92 and 93 as follows; 90 and 92% PDK1 inhibition, respectively.

CONCLUSION

PDK1 inhibitors are potential anticancer agents that can be discovered by combination method of ligand based design with QSAR and ROC analysis.

Atom Type Independent Modeling of the Conformational Energy of Benzylic, Allylic, and Other Bonds Adjacent to Conjugated Systems

Applications of computational methods to predict binding affinities for protein/drug complexes are routinely used in structure-based drug discovery. Applications of these methods often rely on empirical force fields (FFs) and their associated parameter sets and atom types. However, it is widely accepted that FFs cannot accurately cover the entire chemical space of drug-like molecules, due to the restrictive cost of parametrization and the poor transferability of existing parameters. To address these limitations, initiatives have been carried out to develop more transferable methods, in order to allow for more rigorous descriptions of any drug-like molecule. We have previously reported H-TEQ, a method which does not rely on atom types and incorporates well established chemical principles to assign parameters to organic molecules. The previous implementation of H-TEQ (a torsional barrier prediction method) only covered saturated and lone pair containing molecules; here, we report our efforts to incorporate conjugated systems into our model. The next step was the evaluation of the introduction of unsaturations. The developed model (H-TEQ3.0) has been validated on a wide variety of molecules containing heteroaromatic groups, alkyls, and fused ring systems. Our method performs on par with one of the most commonly used FFs (GAFF2), without relying on atom types or any prior parametrization.

Discovery of Novel DPP-IV Inhibitors as Potential Candidates for the Treatment of Type 2 Diabetes mellitus Predicted by 3D QSAR Pharmacophore Models, Molecular Docking and de novo Evolution

Dipeptidyl peptidase-IV (DPP-IV) rapidly breaks down the incretin hormones glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic peptide (GIP). Thus, the use of DPP-IV inhibitors to retard the degradation of endogenous GLP-1 is a possible mode of therapy correcting the defect in incretin-related physiology. The aim of this study is to find a new small molecule and explore the inhibition activity to the DPP-IV enzyme using a computer aided simulation. In this study, the predicted compounds were suggested as potent anti-diabetic candidates. Chosen structures were applied following computational strategies: The generation of the three-dimensional quantitative structure-activity relationship (3D QSAR) pharmacophore models, virtual screening, molecular docking, and de novo Evolution. The method also validated by performing re-docking and cross-docking studies of seven protein systems for which crystal structures were available for all bound ligands. The molecular docking experiments of predicted compounds within the binding pocket of DPP-IV were conducted. By using 25 training set inhibitors, ten pharmacophore models were generated, among which hypo1 was the best pharmacophore model with the best predictive power on account of the highest cost difference (352.03), the lowest root mean squared deviation (RMSD) (2.234), and the best correlation coefficient (0.925). Hypo1 pharmacophore model was used for virtual screening. A total of 161 compounds including 120 from the databases, 25 from the training set, 16 from the test set were selected for molecular docking. Analyzing the amino acid residues of the ligand-receptor interaction, it can be concluded that Arg125, Glu205, Glu206, Tyr547, Tyr662, and Tyr666 are the main amino acid residues. The last step in this study was de novo Evolution that generated 11 novel compounds. The derivative dpp4_45_Evo_1 by all scores CDOCKER_ENERGY (CDOCKER, -41.79), LigScore1 (LScore1, 5.86), LigScore2 (LScore2, 7.07), PLP1 (-112.01), PLP2 (-105.77), PMF (-162.5)—have exceeded the control compound. Thus the most active compound among 11 derivative compounds is dpp4_45_Evo_1. Additionally, for derivatives dpp4_42_Evo_1, dpp4_43_Evo2, dpp4_46_Evo_4, and dpp4_47_Evo_2, significant upward shifts were recorded. The consensus score for the derivatives of dpp4_45_Evo_1 from 1 to 6, dpp4_43_Evo2 from 4 to 6, dpp4_46_Evo_4 from 1 to 6, and dpp4_47_Evo_2 from 0 to 6 were increased. Generally, predicted candidates can act as potent occurring DPP-IV inhibitors given their ability to bind directly to the active sites of DPP-IV. Our result described that the 6 re-docked and 27 cross-docked protein-ligand complexes showed RMSD values of less than 2 Å. Further investigation will result in the development of novel and potential antidiabetic drugs.

Novel 1,3,4-Triaryl Pyrazoles: Synthesis, QSAR Studies and Cytotoxicity against Breast Cancer

Background:

The existence of drug-resistance and lack of selectivity encourages scientists to search for novel and more selective cytotoxic agents.

Objective:

In this work, novel 1,3,4-triarylpyrazole derivatives were synthesized to study their cytotoxicity on MCF7 (human breast Cell Line). In addition, QSAR studies were performed to show the relation between the cytotoxic activity and the structural features of our new synthesized pyrazole derivatives.

Methods:

Pyrazole-4-carbaldehyde derivative 3 was utilized as a starting material for the preparation of the new pyarazole derivatives. These target compounds were screened for their cytotoxic activity against MCF-7 followed by study cell cycle of the most active compounds. Finally, pharmacophore modeling and QSAR Studies was carried out.

Results:

Among these compounds; 5d and 8b showed the highest anti-proliferative activity (IC50 = 4.9 and 2.11 µM, respectively). Flow cytometric analysis showed that, compounds 5d and 8b arrested the cell cycle in addition to induction of apoptosis in MCF7 cells. Moreover, their stimulation effect on caspases 3/7 was examined to explore their mechanism of induction of apoptosis and the results showed that their proapoptotic activity could be due to the activation of caspases 3/7.

Conclusion:

Pyrazole derivatives 5d and 8b displayed potent bioactivities, indicating that these compounds could be considered as a new lead for more investigation in the future

Activity and selectivity cliffs for DPP-IV inhibitors: Lessons we can learn from SAR studies and their application to virtual screening

The inhibition of dipeptidyl peptidase-IV (DPP-IV) has emerged over the last decade as one of the most effective treatments for type 2 diabetes mellitus, and consequently (a) 11 DPP-IV inhibitors have been on the market since 2006 (three in 2015), and (b) 74 noncovalent complexes involving human DPP-IV and drug-like inhibitors are available at the Protein Data Bank (PDB). The present review aims to (a) explain the most important activity cliffs for DPP-IV noncovalent inhibition according to the binding site structure of DPP-IV, (b) explain the most important selectivity cliffs for DPP-IV noncovalent inhibition in comparison with other related enzymes (i.e., DPP8 and DPP9), and (c) use the information deriving from this activity/selectivity cliff analysis to suggest how virtual screening protocols might be improved to favor the early identification of potent and selective DPP-IV inhibitors in molecular databases (because they have not succeeded in identifying selective DPP-IV inhibitors with IC₅₀ ≤ 100 nM). All these goals are achieved with the help of available homology models for DPP8 and DPP9 and an analysis of the structure-activity studies used to develop the noncovalent inhibitors that form part of some of the complexes with human DPP-IV available at the PDB.

Ligand-based computer aided drug design reveals new tropomycin receptor kinase a (TrkA) inhibitors

Targeting tropomycin kinase A (TrkA) by small molecule inhibitors is considered as a promising strategy for treating several human cancers. To achieve this goal, a ligand based QSAR model was applied using the Discovery studio 4.5 (DS 4.5). Hence, a total list of 161 TrkA inhibitors was investigated. The TrkA inhibitors were extensively explored to detect their optimal physicochemical properties and pharmacophoric binding modes, which were converted into numeric descriptors and allowed to compete within the context of the Genetic Function Algorithm (GFA) approximations to find the subset of terms that correlates best with the activity. The resulted successful QSAR equation had statistical criteria of (r²₁₂₉=0.67, r²_LOO=0.61 r²_PRESS against 32 external test inhibitors=0.50). Afterwards, the most successful pharmacophore: HypoB-T5-3, was used to screen compounds within the National Cancer institute (NCI) database. Only 41 compounds were retrieved and 21 of them exhibited anti-TrkA activity. The most potent hit had an IC₅₀ value of 2.4μM. Later, upon docking the active hits into the TrkA binding pocket, important interactions were revealed including hydrogen bonding with the amino acids Asp668 and Lys544 in addition to the cation-π interactions with the sidechain of Arg559.

Ligand-based modeling of Akt3 lead to potent dual Akt1/Akt3 inhibitor

Akt1 and Akt3 are important serine/threonine-specific protein kinases involved in G2 phase required by cancer cells to maintain cell cycle and to prevent cell death. Accordingly, inhibitors of these kinases should have potent anti-cancer properties. This prompted us to use pharmacophore/QSAR modeling to identify optimal binding models and physicochemical descriptors that explain bioactivity variation within a set of 74 diverse Akt3 inhibitors. Two successful orthogonal pharmacophores were identified and further validated using receiver operating characteristic (ROC) curve analyses. The pharmacophoric models and associated QSAR equation were applied to screen the national cancer institute (NCI) list of compounds for new Akt3 inhibitors. Six hits showed significant experimental anti-Akt3 IC₅₀ values, out of which one compound exhibited dual low micromolar anti-Akt1 and anti-Akt3 inhibitory profiles.

Structural Basis of Bioactivity of Food Peptides in Promoting Metabolic Health

Bioactive peptides have many structural features that enable them to become functional in controlling several biological processes in the body, especially those related to metabolic health. This chapter provides an overview of the multiple targets of food-derived peptides against metabolic health problems (e.g., hypertension, dyslipidemia, hyperglycemia, oxidative stress) and discusses the importance of structural chemistry in determining the bioactivities of peptides and protein hydrolysates.

Discovery of novel potent nuclear factor kappa-B inhibitors (IKK-β) via extensive ligand-based modeling and virtual screening

Inhibitor kappa-B kinase-beta (IKK-β) controls the activation of nuclear transcription factor kappa-B and has been linked to inflammation and cancer. Therefore, inhibitors of this kinase should have potent anti-inflammatory and anticancer properties. Accordingly, we explored the pharmacophoric space of 218 IKK-β inhibitors to identify high-quality binding models. Subsequently, genetic algorithm-based quantitative structure activity relationship (QSAR) analysis was employed to select the best possible combination of pharmacophoric models and physicochemical descriptors that explain bioactivity variation among training compounds. Three successful pharmacophores emerged in 2 optimal QSAR equations (r₁²₁₇₅  = 0.733, r₁²_LOO  = 0.52, F1 = 65.62, r₁²_PRESS against 43 test inhibitors = 0.63 and r₂²₁₇₅  = 0.683, r₂²_LOO  = 0.52, F2 = 72.66, r₂²_PRESS against 43 test inhibitors = 0.65). Two pharmacophores were merged in a single binding model. Receiver operating characteristic curve validation proved the excellent qualities of this model. The merged pharmacophore and the associated QSAR equations were applied to screen the National Cancer Institute list of compounds. Ten hits were found to exhibit potent anti-IKK-β bioactivity, out of which, one illustrates IC50 of 11.0nM.

Computer-aided discovery of new FGFR-1 inhibitors followed by in vitro validation

AIM

FGFR-1 is an oncogenic kinase involved in several cancers. FGFR1-specific inhibitors have shown promising results against several human cancers prompting us to model this interesting target. Toward the end, we implemented elaborate ligand-based and structure-based computational workflows to explore the pharmacophoric requirements for potent FGFR-1 inhibitors. Results & methodology: Structure-based and ligand-based modeling applied on 59 diverse FGFR-1 inhibitors yielded novel pharmacophore and quantitative structure-activity relationship models that were used to scan the National Cancer Institute's structural database for novel leads. Four potent hits were captured, with the most active having IC₅₀ of 426 nM. Identities and purities of active hits were established using nuclear magnetic resonance and mass spectroscopy.

CONCLUSION

Elaborate ligand-based (pharmacophore/quantitaive structure-activity relationship) and structure-based (docking-based comparative intermolecular contacts analysis) modeling provided deep understanding of ligand binding within FGFR-1 as evidenced by the virtually captured new potent leads.

Ligand-based modeling of diverse aryalkylamines yields new potent P-glycoprotein inhibitors

The P-glycoprotein (P-gp) efflux pump has an important role as a natural detoxification system in many types of normal and cancer cells. P-gp is implicated in multiple drug resistance (MDR) exhibited by several types of cancer against a multitude of anticancer chemotherapeutic agents, and therefore, it is clinically validated target for cancer therapy. Accordingly, in this study we combined exhaustive pharmacophore modeling and quantitative structure-activity relationship (QSAR) analysis to explore the structural requirements for potent P-gp inhibitors employing 130 known P-gp ligands. Genetic function algorithm (GFA) coupled with k nearest neighbor (kNN) or multiple linear regression (MLR) analyses were employed to build self-consistent and predictive QSAR models based on optimal combinations of pharmacophores and physicochemical descriptors. Successful pharmacophores were complemented with exclusion spheres to optimize their receiver operating characteristic curve (ROC) profiles. Optimal QSAR models and their associated pharmacophore hypotheses were validated by identification and experimental evaluation of new promising P-gp inhibitory leads retrieved from the National Cancer Institute (NCI) structural database. Several potent hits were captured. The most potent hit decreased the IC50 of doxorubicin from 0.906 to 0.190 μM on doxorubicin resistant MCF7 cell-line.

Identification of novel inhibitors for Pim-1 kinase using pharmacophore modeling based on a novel method for selecting pharmacophore generation subsets

Targeting Proviral integration-site of murine Moloney leukemia virus 1 kinase, hereafter called Pim-1 kinase, is a promising strategy for treating different kinds of human cancer. Headed for this a total list of 328 formerly reported Pim-1 kinase inhibitors has been explored and divided based on the pharmacophoric features of the most active molecules into 10 subsets projected to represent potential active binding manners accessible to ligands within the binding pocket of Pim-1 kinase. Discovery Studio 4.1 (DS 4.1) was employed to detect potential pharmacophoric active binding manners anticipated by Pim-1 Kinase inhibitors. The pharmacophoric models were then allowed to compete within Quantitative Structure Activity Relationship (QSAR) framework with other 2D descriptors. Accordingly Genetic algorithm and multiple linear regression investigation were engaged to find the finest QSAR equation that has the best predictive power r262(2) = 0.70, F = 119.14, rLOO(2) = 0.693, rPRESS(2) against 66 external test inhibitors = 0.71 q(2) = 0.55. Three different pharmacophores appeared in the successful QSAR equation this represents three different binding modes for inhibitors within the Pim-1 kinase binding pocket. Pharmacophoric models were later used to screen compounds within the National Cancer Institute database. Several low micromolar Pim-1 Kinase inhibitors were captured. The most potent hits show IC50 values of 0.77 and 1.03 µM. Also, upon analyzing the successful QSAR Equation we found that some polycyclic aromatic electron-rich structures namely 6-Chloro-2-methoxy-acridine can be considered as putative hits for Pim-1 kinase inhibition.

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.

Discovery of nanomolar phosphoinositide 3-kinase gamma (PI3Kγ) inhibitors using ligand-based modeling and virtual screening followed by in vitro analysis

Phosphoinositide 3-kinase gamma (PI3Kγ) is member of a family of enzymes involved in cancer pathogenesis. Accordingly, considerable efforts have been carried out to develop new PI3Kγ inhibitors. Towards this end we explored the pharmacophoric space of PI3Kγ using three diverse sets of inhibitors. Subsequently, we employed genetic algorithm-based QSAR analysis to select optimal combination of pharmacophoric models and physicochemical descriptors that can explain bioactivity variation within training inhibitors. Interestingly, two successful pharmacophores were selected within two statistically consistent QSAR models. The close similarity among the two binding models prompted us to merge them in a hybrid pharmacophore. The resulting model showed superior receiver operator characteristic curve (ROC) and closely resembled binding interactions seen in crystallographic ligand-PI3Kγ complexes. The resulting model was employed to screen the national cancer institute (NCI) list of compounds to search for new PI3Kγ ligands. After testing captured hits in vitro, 19 compounds showed nanomolar IC50 values against PI3Kγ. The chemical structures and purities of most potent hits were validated using NMR and MS experiments.

The use of docking-based comparative intermolecular contacts analysis to identify optimal docking conditions within glucokinase and to discover of new GK activators

Glucokinase (GK) is involved in normal glucose homeostasis and therefore it is a valid target for drug design and discovery efforts. GK activators (GKAs) have excellent potential as treatments of hyperglycemia and diabetes. The combined recent interest in GKAs, together with docking limitations and shortages of docking validation methods prompted us to use our new 3D-QSAR analysis, namely, docking-based comparative intermolecular contacts analysis (dbCICA), to validate docking configurations performed on a group of GKAs within GK binding site. dbCICA assesses the consistency of docking by assessing the correlation between ligands' affinities and their contacts with binding site spots. Optimal dbCICA models were validated by receiver operating characteristic curve analysis and comparative molecular field analysis. dbCICA models were also converted into valid pharmacophores that were used as search queries to mine 3D structural databases for new GKAs. The search yielded several potent bioactivators that experimentally increased GK bioactivity up to 7.5-folds at 10 μM.

Elaborate ligand-based modeling coupled with QSAR analysis and in silico screening reveal new potent acetylcholinesterase inhibitors

Inhibition of the enzyme acetylcholinesterase (AChE) has been shown to alleviate neurodegenerative diseases prompting several attempts to discover and optimize new AChE inhibitors. In this direction, we explored the pharmacophoric space of 85 AChE inhibitors to identify high quality pharmacophores. Subsequently, we implemented genetic algorithm-based quantitative structure-activity relationship (QSAR) modeling to select optimal combination of pharmacophoric models and 2D physicochemical descriptors capable of explaining bioactivity variation among training compounds (r2(68)=0.94, F-statistic=125.8, r2 LOO=0.92, r2 PRESS against 17 external test inhibitors = 0.84). Two orthogonal pharmacophores emerged in the QSAR equation suggesting the existence of at least two binding modes accessible to ligands within AChE binding pocket. The successful pharmacophores were comparable with crystallographically resolved AChE binding pocket. We employed the pharmacophoric models and associated QSAR equation to screen the national cancer institute list of compounds. Twenty-four low micromolar AChE inhibitors were identified. The most potent gave IC50 value of 1.0 μM.

Synthesis of Some 4,5-Dihydrothieno[3,2-e][1,2,4]Triazolo[4,3-a] Pyrimi-dine-2-Carboxamides as Anti-Inflammatory and Analgesic Agents

A new series 4,5-dihydrothieno[3,2-e][1,2,4]triazolo[4,3-a]pyrimidine-2-carboxamide was synthesized. Twenty one newly synthesized compounds were investigated for their anti-inflammatory and analgesic activity using acute and subacute formalin-induced paw edema models and diclofenac Na as a reference. The acute toxicity (ALD50) and ulcerogenic effects of the active compounds were also determined. The thienotriazolopyrimidines 10a, 10c and 11c were found to exhibit remarkable anti-inflammatory activity at both models in addition to good analgesic activity with a delayed onset of action. Moreover, the active compounds showed high GI safety level and are well tolerated by experimental animals with high safety margin (ALD50 > 0.4 g/kg). Docking study using Molecular Operating Environment (MOE) version 2008.10 into COX-2 has been made for derivatives of highest anti-inflammatory activity.

New leads for DPP IV inhibition: structure-based pharmacophore mapping and virtual screening study

Dipeptidyl peptidase IV (DPP IV) is an attractive target for the development of new antidiabetic drugs. DPP IV inhibitors improve glycemic control by preventing the rapid inactivation of the incretin hormones; glucagon-like peptide 1 (GLP-1) and glucose-dependent insulinotropic peptide. In the current study, virtual screening, using 2D and 3D filters implemented in a hierarchical cascade, was employed to identify new DPP IV inhibitors. Co-crystallized ligands, with potent DPP IV-inhibitory activities, were utilized to generate structure-based pharmacophore models using DS Visualizer software. The derived pharmacophore maps were validated using in-house built database containing active and inactive DPP IV inhibitors. Subsequently, the optimum validated pharmacophore model was used as a search query against two 3D-databases (NCI and in-house built drug databases). Further hit filtration was carried out employing 2D virtual filters based on Lipinski's rule of 5; number of rotatable bonds and other physicochemical filters. 3D filter using high-throughput molecular docking was also applied. As a result, 5 novel DPP IV inhibitors were discovered as potential lead compounds and later confirmed via in vitro bioassay.

Elaborate ligand-based modeling reveal new migration inhibitory factor inhibitors

Recent research suggested the involvement of migration inhibitor factor (MIF) in cancer and inflammatory diseases, which prompted several attempts to develop new MIF inhibitors. Accordingly, we investigated the pharmacophoric space of 79 MIF inhibitors using seven diverse subsets of inhibitors to identify plausible binding hypotheses (pharmacophores). Subsequently, we implemented genetic algorithm and multiple linear regression analysis to select optimal combination of pharmacophores and physicochemical descriptors capable of explaining bioactivity variation within the training compounds (QSAR model, r63=0.62, F=42.8, rLOO(2)=0.721,rPRESS(2) against 16 external test inhibitors=0.58). Two orthogonal pharmacophores appeared in the optimal QSAR model suggestive of at least two binding modes available to ligands inside MIF binding pocket. Subsequent validation using receiver operating characteristic (ROC) curves analysis established the validity of these two pharmacophores. We employed these pharmacophoric models and associated QSAR equation to screen the National Cancer Institute (NCI) list of compounds. Eight compounds gave >50% inhibition at 100μM. Two molecules illustrated >75% inhibition at 10μM.

DPP-IV inhibitory potential of naringin: an in silico, in vitro and in vivo study

The incretin based therapies are an emerging class of antidiabetic drugs, with two categories: one is glucagone like peptide-1 (GLP-1) agonists and the other one is dipeptidyl peptidase (CD26; DPP-IV) inhibitors. However, in the DPP-IV inhibitors category only few compounds are commercially available and also have some undesirable effects. Therefore, in the present work we tried to explore a naturally occurring compound naringin for its potential DPP-IV inhibition and antidiabetic potential. It is noteworthy that this compound is abundantly present in orange peel and thus may provide cost effective treatment for diabetes, especially type 2 diabetes mellitus. In the present study, we have conducted virtual docking study and observed tight binding of naringin, as shown by higher negative values of H bond lengths, while in vitro DPP-IV inhibition assay has also shown better inhibition by naringin. In vivo study, in response to 10 days administration of 40 mg/kg of naringin twice daily to Wistar albino rats, inhibited the serum levels of DPP-IV activity, random glucose concentration with concomitant increase in insulin levels. All the comparisons were made with the standard commercially available drug sitagliptin.

Application of docking-based comparative intermolecular contacts analysis to validate Hsp90α docking studies and subsequent in silico screening for inhibitors

Heat shock protein (Hsp90α) has been recently implicated in cancer, prompting several attempts to discover and optimize new Hsp90α inhibitors. Towards this end, we docked 83 diverse Hsp90α inhibitors into the ATP-binding site of this chaperone using several docking-scoring settings. Subsequently, we applied our newly developed computational tool--docking-based comparative intramolecular contacts analysis (dbCICA)--to assess the different docking conditions and select the best settings. dbCICA is based on the number and quality of contacts between docked ligands and amino acid residues within the binding pocket. It assesses a particular docking configuration based on its ability to align a set of ligands within a corresponding binding pocket in such a way that potent ligands come into contact with binding site spots distinct from those approached by low-affinity ligands, and vice versa. The optimal dbCICA models were translated into valid pharmacophore models that were used as 3D search queries to mine the National Cancer Institute's structural database for new inhibitors of Hsp90α that could potentially be used as anticancer agents. The process culminated in 15 micromolar Hsp90α ATPase inhibitors.

Design of novel FLT-3 inhibitors based on dual-layer 3D-QSAR model and fragment-based compounds in silico

FMS-like tyrosine kinase 3 (FLT-3) is strongly correlated with acute myeloid leukemia, but no FLT-3-inhibitor cocomplex structure is available to assist the design of therapeutic inhibitors. Hence, we propose a dual-layer 3D-QSAR model for FLT-3 that integrates the pharmacophore, CoMFA, and CoMSIA. We then coupled the model with the fragment-based design strategy to identify novel FLT-3 inhibitors. In the first layer, the previously established model, Hypo02, was evaluated in terms of its correlation coefficient (r), RMS, cost difference, and configuration cost, with values of 0.930, 1.24, 106.45, and 16.44, respectively. Moreover, Fischer's cross-validation test of data generated by Hypo02 yielded a 98% confidence level, and the validation of the testing set yielded a best r value of 0.87. The features of Hypo02 were separated into two parts and then used to screen the MiniMaybridge fragment compound database. Nine novel FLT-3 inhibitors were generated in this layer. In the second layer, Hypo02 was subjected to an alignment rule to generate CoMFA- and CoMSIA-based models, for which the partial least-squares validation method was utilized. The values of q(2), r(2), and predictive r(2) were 0.58, 0.98, and 0.76, respectively, derived from the CoMFA model with steric and electrostatic fields. The CoMSIA model with five different fields yielded values of 0.54, 0.97, and 0.76 for q(2), r(2), and predictive r(2), respectively. The CoMFA and CoMSIA models were used to constrain 3D structures of the nine novel FLT-3 inhibitors. This dual-layer 3D-QSAR model constitutes a valuable tool to easily and quickly screen and optimize novel potential FLT-3 inhibitors for the treatment of acute myeloid leukemia.

Synthesis and QSAR study of novel cytotoxic spiro[3H-indole-3,2'(1'H)-pyrrolo[3,4-c]pyrrole]-2,3',5'(1H,2'aH,4'H)-triones

1,3-Dipolar cycloaddition reaction of 1-aryl-1H-pyrrole-2,5-diones 1a-e with non-stabilized azomethine ylides, generated in situ via decarboxylative condensation of isatins 2a-c and sarcosine (3) in refluxing ethanol, afforded 4'-aryl-5'a,6'-dihydro-1'-methyl-spiro[3H-indole-3,2'(1'H)-pyrrolo[3,4-c]pyrrole]-2,3',5'(1H,2'aH,4'H)-triones 4a-o in good yields. Compound 4l exhibited high anti-tumor activity against HEPG2 (liver cancer) cell line (IC(50) = 12.16 μM) compared to that of Doxorubicin (IC(50) = 7.36 μM), and the other synthesized compounds revealed moderate anti-tumor properties against HCT116 (colon), MCF7 (breast) and HEPG2 (liver) human tumor cell lines. 3D-Pharmacophore modeling and quantitative structure-activity relationship (QSAR) analysis were combined to explore the structural requirements controlling the observed anti-tumor properties. It was found that the major structural factors affecting potency of these compounds were related to their basic skeleton.

Docking-based comparative intermolecular contacts analysis as new 3-D QSAR concept for validating docking studies and in silico screening: NMT and GP inhibitors as case studies

The significant role played by docking algorithms in drug discovery combined with their serious pitfalls prompted us to envisage a novel concept for validating docking solutions, namely, docking-based comparative intermolecular contacts analysis (dbCICA). This novel approach is based on the number and quality of contacts between docked ligands and amino acid residues within the binding pocket. It assesses a particular docking configuration on the basis of its ability to align a set of ligands within a corresponding binding pocket in such a way that potent ligands come into contact with binding site spots distinct from those approached by low-affinity ligands and vice versa. In other words, dbCICA evaluates the consistency of docking by assessing the correlation between ligands' affinities and their contacts with binding site spots. Optimal dbCICA models can be translated into valid pharmacophore models that can be used as 3-D search queries to mine structural databases for new bioactive compounds. dbCICA was implemented to search for new inhibitors of candida N-myristoyl transferase as potential antifungal agents and glycogen phosphorylase (GP) inhibitors as potential antidiabetic agents. The process culminated in five selective micromolar antifungal leads and nine GP inhibitory leads.

Elaborate ligand-based modeling reveals new nanomolar heat shock protein 90α inhibitors

Heat shock protein (Hsp90α) has been recently implicated in cancer prompting several attempts to discover and optimize new Hsp90α inhibitors. Toward this end, we explored the pharmacophoric space of 83 Hsp90α inhibitors using six diverse sets of inhibitors to identify high-quality pharmacophores. Subsequently, genetic algorithm and multiple linear regression analysis were employed to select an optimal combination of pharmacophoric models and 2D physicochemical descriptors capable of accessing a self-consistent quantitative structure activity relationship (QSAR) of optimal predictive potential (r(67)(2)=0.811, F 42.8, r(LOO)(2)=0.748, r(PRESS)(2) (against 16 external test inhibitors) = 0.619). Three orthogonal pharmacophores emerged in the QSAR equation suggesting the existence of at least three binding modes accessible to ligands within the Hsp90α binding pocket. Receiver operating characteristic (ROC) curves analysis established the validity of QSAR-selected pharmacophores. We employed the pharmacophoric models and associated QSAR equation to screen the national cancer institute (NCI) list of compounds and our in-house-built drugs and agrochemicals database (DAC). Twenty-five nanomolar and low micromolar Hsp90α inhibitors were identified. The most potent were formoterol, amodaquine, primaquine, and midodrine with IC(50) values of 3, 5, 6, and 20 nM, respectively.