Virtual screening of dipeptidyl peptidase-4 inhibitors using quantitative structure–activity relationship-based artificial intelligence and molecular docking of hit compounds

Dipeptidyl peptidase-4 (DPP-4) inhibitors are becoming an essential drug in the treatment of type 2 diabetes mellitus; however, some classes of these drugs exert side effects, including joint pain and pancreatitis. Studies suggest that these side effects might be related to secondary inhibition of DPP-8 and DPP-9. In this study, we identified DPP-4-inhibitor hit compounds selective against DPP-8 and DPP-9. We built a virtual screening workflow using a quantitative structure–activity relationship (QSAR) strategy based on artificial intelligence to allow faster screening of millions of molecules for the DPP-4 target relative to other screening methods. Five regression machine learning algorithms and four classification machine learning algorithms were applied to build virtual screening workflows, with the QSAR model applied using support vector regression (R²_pred 0.78) and the classification QSAR model using the random forest algorithm with 92.2% accuracy. Virtual screening results of > 10 million molecules obtained 2 716 hits compounds with a pIC₅₀ value of > 7.5. Additionally, molecular docking results of several potential hit compounds for DPP-4, DPP-8, and DPP-9 identified CH0002 as showing high inhibitory potential against DPP-4 and low inhibitory potential for DPP-8 and DPP-9 enzymes. These results demonstrated the effectiveness of this technique for identifying DPP-4-inhibitor hit compounds selective for DPP-4 and against DPP-8 and DPP-9 and suggest its potential efficacy for applications to discover hit compounds of other targets.     

Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.     

基于药效团模型的DHODH抑制剂构效关系研究

利用药效团模型研究二氢乳清酸脱氢酶(Dihydroorotate dehydrogenase,DHODH)抑制剂的构效关系,为DHODH抑制剂的虚拟筛选提供新的方法.以31个具有DHODH抑制活性的化合物为训练集化合物,半数抑制浓度(IC50)范围为7~63000 nmol/L,利用Catalyst/HypoGen算法构建DHODH抑制剂药效团模型,通过对训练集化合物多个构象进行叠合,提取药效团特征及三维空间限制构建药效团模型.利用基于CatScramble的交叉验证方法及评价模型对已知活性化合物的活性预测能力,确定较优药效团模型.模型包含1个氢键受体、3个疏水中心,表征了受体配体相互作用时可能发生的氢键相互作用、疏水相互作用和π-π相互作用,4个药效特征在三维空间的排列概括了DHODH抑制剂产生活性的结构特点.所得较优模型对训练集化合物及测试集化合物的计算活性值与实验活性值的相关系数分别为0.8405和0.8788.利用药效团模型对来源于微生物的系列化合物进行虚拟筛选,筛选出59个预测活性较好的化合物,可作为进一步药物研发的候选化合物.     

Pharmacophore modeling,virtual screening and molecular docking of ATPase inhibitors of HSP70

Heat shock protein 70 is an effective anticancer target as it influences many signaling pathways. Hence the study investigated the important pharmacophore feature required for ATPase inhibitors of HSP70 by generating a ligand based pharmacophore model followed by virtual based screening and subsequent validation by molecular docking in Discovery studio V4.0. The most extrapolative pharmacophore model (hypotheses 8) consisted of four hydrogen bond acceptors. Further validation by external test set prediction identified 200 hits from Mini Maybridge, Drug Diverse, SCPDB compounds and Phytochemicals. Consequently, the screened compounds were refined by rule of five, ADMET and molecular docking to retain the best competitive hits. Finally Phytochemical compounds Muricatetrocin B, Diacetylphiladelphicalactone C, Eleutheroside B and 5-(3-{[1-(benzylsulfonyl)piperidin-4-yl]amino}phenyl)- 4-bromo-3-(carboxymethoxy)thiophene-2-carboxylic acid were obtained as leads to inhibit the ATPase activity of HSP70 in our findings and thus can be proposed for further in vitro and in vivo evaluation.     

Design and synthesis of novel diamide derivatives of glycine as antihyperglycemic agents

DPP-4 inhibition is one of the most extensively explored approaches for the management of type 2 diabetes (T2D). Most DPP-4 inhibitors in the market contain a proline mimetic active pharmacophore. Herein, we report the design, synthesis, and preliminary evaluation of a series of novel diamide derivatives of glycine, devoid of the proline mimic, for the treatment of T2D. As predicted from in silico studies, the diamide derivatives of glycine exhibited comparable DPP-4 inhibition with the standard as confirmed by the preliminary in vitro studies. Compound 6b was found to be the most potent (IC₅₀ 94.82 nM) DPP-4 inhibitor among all the molecules synthesized in the series.     

Pharmacophore modelling and electronic feature analysis of hydroxamic acid derivatives,the HIV integrase inhibitors

Hydroxamic acid derivatives with metal ion binding properties were collected from the literature to generate a pharmacophore and 3D-QSAR model for HIV strand transfer inhibition. The derived pharmacophore model (AAAHRR) recognizes both metal ion binding site and hydrophobic group. The QSAR model generated using this hypothesis expressed statistical significance (r ² = 0.971 for the training set and q ² = 0.913 for the test set). The ability of this pharmacophore model to retrieve other metal ion binding inhibitors was examined by screening the ChemBank database (ligandinfo) incorporated with 10 known strand transfer inhibitors. The studied favourable and unfavourable contours of chemical features (H-bond donor, acceptor and hydrophobic sites) revealed the role of hydrophobic substitution at the fluorobenzene ring and cyclization of the metal ion binding hydroxamic acid in effective integrase inhibition. Analysis of the frontier orbitals, HOMO and LUMO revealed that the nucleophilic / electrophilic interactions depend on the significant overlapping observed at the azaindole and hydroxamic acid groups. In essence, the generated pharmacophore model is competent enough to disclose the essential site-specific interactions involved in the inhibition of HIV integrase, and hence can be used in virtual screening to identify novel scaffolds as leads with increased anti-viral potency.     

Multiple receptor-ligand based pharmacophore modeling and molecular docking to screen the selective inhibitors of matrix metalloproteinase-9 from natural products

Matrix metalloproteinase-9 (MMP-9) is an attractive target for cancer therapy. In this study, the pharmacophore model of MMP-9 inhibitors is built based on the experimental binding structures of multiple receptor-ligand complexes. It is found that the pharmacophore model consists of six chemical features, including two hydrogen bond acceptors, one hydrogen bond donor, one ring aromatic regions, and two hydrophobic (HY) features. Among them, the two HY features are especially important because they can enter the S1′ pocket of MMP-9 which determines the selectivity of MMP-9 inhibitors. The reliability of pharmacophore model is validated based on the two different decoy sets and relevant experimental data. The virtual screening, combining pharmacophore model with molecular docking, is performed to identify the selective MMP-9 inhibitors from a database of natural products. The four novel MMP-9 inhibitors of natural products, NP-000686, NP-001752, NP-014331, and NP-015905, are found; one of them, NP-000686, is used to perform the experiment of in vitro bioassay inhibiting MMP-9, and the IC₅₀ value was estimated to be only 13.4 µM, showing the strongly inhibitory activity of NP-000686 against MMP-9, which suggests that our screening results should be reliable. The binding modes of screened inhibitors with MMP-9 active sites were discussed. In addition, the ADMET properties and physicochemical properties of screened four compounds were assessed. The found MMP-9 inhibitors of natural products could serve as the lead compounds for designing the new MMP-9 inhibitors by carrying out structural modifications in the future.     

Construction of the Pharmacophore Model of Glycogen Synthase Kinase-3 Inhibitors

A three dimensional pharmacophore model has been generated for glycogen synthase kinase-3 （GSK-3） inhibitots. A dataset consisting of 89 compounds was selected on the basis of the information content of the structures and activity data as required by the CATALYST system. The optimum model with four features （one hydrogen-bond acceptor unit, one ring aromatic unit, and two hydrophobic aromatic units） was selected with a good correlation coefficient （0.95）. This model is able to predict the activity of other known GSK-3 inhibitors not included in the model generation, and can be further used to identify structurally diverse compounds with desired biological activity by virtual screening.     

Identification of potential Tpx inhibitors against pathogen-host interactions

Yersinia organisms cause many infectious diseases by invading human cells and delivering their virulence factors via the type three secretion system (T3SS). One alternative strategy in the fight against these pathogenic organisms is to interfere with their T3SS. Previous studies demonstrated that thiol peroxidase, Tpx is functional in the assembly of T3SS and its inhibition by salicylidene acylhydrazides prevents the secretion of pathogenic effectors. In this study, the aim was to identify potential inhibitors of Tpx using an integrated approach starting with high throughput virtual screening and ending with molecular dynamics simulations of selected ligands. Virtual screening of ZINC database of 500,000 compounds via ligand-based and structure-based pharmacophore models retrieved 10,000 hits. The structure-based pharmacophore model was validated using high-throughput virtual screening (HTVS). After multistep docking (SP and XP), common scaffolds were used to find common substructures and the ligand binding poses were optimized using induced fit docking. The stability of the protein–ligand complex was examined with molecular dynamics simulations and the binding free energy of the complex was calculated. As a final outcome eight compounds with different chemotypes were proposed as potential inhibitors for Tpx. The eight ligands identified by a detailed virtual screening protocol can serve as leads in future drug design efforts against the destructive actions of pathogenic bacteria.     

Pharmacophore Modeling and in Silico Screening Studies to Design Potential KDR Kinase Inhibitors

A novel ligand‐based pharmacophore model for KDR kinase was generated on the basis of chemical features of 30 KDR kinase inhibitors. This pharmacophore model consists of one hydrogen‐bond acceptor, one hydrogen‐bond donor and two hydrophobic groups. Several methods have been used to validate the model, suggesting that it can serve as a reliable tool for virtual screening to facilitate the discovery of novel KDR inhibitors. The model was then used as database search query from the National Cancer Institute (NCI) database for the rational design to identify new hit compound.     

Identification of potential influenza virus endonuclease inhibitors through virtual screening based on the 3D-QSAR model

Influenza endonucleases have appeared as an attractive target of antiviral therapy for influenza infection. With the purpose of designing a novel antiviral agent with enhanced biological activities against influenza endonuclease, a three-dimensional quantitative structure-activity relationships (3D-QSAR) model was generated based on 34 influenza endonuclease inhibitors. The comparative molecular similarity index analysis (CoMSIA) with a steric, electrostatic and hydrophobic (SEH) model showed the best correlative and predictive capability (q ² = 0.763, r ² = 0.969 and F = 174.785), which provided a pharmacophore composed of the electronegative moiety as well as the bulky hydrophobic group. The CoMSIA model was used as a pharmacophore query in the UNITY search of the ChemDiv compound library to give virtual active compounds. The 3D-QSAR model was then used to predict the activity of the selected compounds, which identified three compounds as the most likely inhibitor candidates.     

基于药效团的BRD4靶向小分子抑制剂虚拟筛选

BRD4靶点和多种肿瘤密切相关,是具有良好成药性的热门靶点。本文选取活性较好且结构差异较大的BRD4小分子抑制剂作为训练集分子,基于配体小分子共同特征(HipHop)方法使用Discovery Studio 3.0分子模拟软件构建了药效团。药效团通过测试集验证、ROC曲线验证(SE(sensitivity)=0.93765、SP(specificity)=0.89500、(AUC)=0.956),结果表明构建得到的药效团具有较强的可靠性和较高的可信度。药效团模型含有1个芳环中心、1个疏水基团、2个氢键受体四个药效特征元素。此药效团被用于ZINC数据库进行虚拟筛选,共筛选了861203个分子,命中率为0.782%。再对筛选得到的分子经过分子对接、ADMET成药性预测、构象分析并讨论分子-蛋白相互作用模式,最终得到了21个有潜力的BRD4小分子抑制剂。     

Virtual screening of B-Raf kinase inhibitors: A combination of pharmacophore modelling,molecular docking, 3D-QSAR model and binding free energy calculation studies

B-Raf kinase has been identified as an important target in recent cancer treatment. In order to discover structurally diverse and novel B-Raf inhibitors (BRIs), a virtual screening of BRIs against ZINC database was performed by using a combination of pharmacophore modelling, molecular docking, 3D-QSAR model and binding free energy (ΔG_bind) calculation studies in this work. After the virtual screening, six promising hit compounds were obtained, which were then tested for inhibitory activities of A375 cell lines. In the result, five hit compounds show good biological activities (IC₅₀ < 50 μM). The present method of virtual screening can be applied to find structurally diverse inhibitors, and the obtained five structurally diverse compounds are expected to develop novel BRIs.     

Identification of Novel β3-Adrenoceptor Agonists Using Energetic Analysis, Structure Based Pharmacophores and Virtual Screening

β3 Adrenergic receptor (β3-AR), is a potential therapeutic target for the treatment of type II diabetes and obesity. We report the identification of novel compounds as β3-AR agonists by integrating different approaches of energetic analysis, structure based pharmacophore designing and virtual screening. In a step wise filtering protocol, structure based virtual screening of 2,33,450 compounds was done. These molecules were docked into the active site of the receptor utilizing three levels of accuracy; ligands passing the HTVS (high throughput virtual screening) step were subsequently analyzed in Glide SP (Standard Precision) and finally in Glide XP (Extra Precision) to estimate the receptor ligand binding affinities. In the second step a total of 300 pharmacophore hypotheses were generated from a set of known and diverse β3-AR agonists. The best hypothesis showed six features: three hydrogen bond acceptors, one positively charged group, and two aromatic rings. To cross validate, pharmacophore filtering was done on the set of shortlisted compounds from structure based VS (virtual screening). The different screening techniques employed were validated using enrichment factor calculations. The energetic based Pharmacophore performed fairly well at distinguishing active from the inactive compounds and yielded a greater diversity of active molecules whereas the number of actives retrieved in the case of structure based screening was the highest.     

Pharmacophore identification, in silico screening, and virtual library design for inhibitors of the human factor Xa

Factor Xa inhibitors are innovative anticoagulant agents that provide a better safety/efficacy profile compared to other anticoagulative drugs. A chemical feature-based modeling approach was applied to identify crucial pharmacophore patterns from 3D crystal structures of inhibitors bound to human factor Xa (Pdb entries 1fjs, 1kns, 1eqz) using the software LIGANDSCOUT and CATALYST. The complex structures were selected regarding the criteria of high inhibitory potency (i.e. all ligands show K(i) values against factor Xa in the subnanomolar range) and good resolution (i.e. at least 2.2 A) in order to generate selective and high quality pharmacophore models. The resulting chemical-feature based hypotheses were used for virtual screening of commercial molecular databases such as the WDI database. Furthermore, a ligand-based molecular modeling approach was performed to obtain common-feature hypotheses that represent the relevant chemical interactions between 10 bioactive factor Xa inhibitors and the protein, respectively. In a next step a virtual combinatorial library was designed in order to generate new compounds with similar chemical and spatial properties as known inhibitors. The software tool ILIB DIVERSE was used for this procedure in order to provide new scaffolds of this group of anticoagulants. Finally we present the combination of these two techniques, hence virtual screening was performed with selective pharmacophore models in a focused virtual combinatorial database. De novo derived molecular scaffolds that were able to adequately satisfy the pharmacophore criteria are revealed and are promising templates for candidates for further development.     

A pseudo-ligand approach to virtual screening

A virtual screening method is presented that is grounded on a receptor-derived pharmacophore model termed "virtual ligand" or "pseudo-ligand". The model represents an idealized constellation of potential ligand sites that interact with residues of the binding pocket. For rapid virtual screening of compound libraries the potential pharmacophore points of the virtual ligand are encoded as an alignment-free correlation vector, avoiding spatial alignment of pharmacophore features between the pharmacophore query (i.e., the virtual ligand) and the candidate molecule. The method was successfully applied to retrieving factor Xa inhibitors from a Ugi three-component combinatorial library, and yielded high enrichment of actives in a retrospective search for cyclooxygenase-2 (COX-2) inhibitors. The approach provides a concept for "de-orphanizing" potential drug targets and identifying ligands for hitherto unexplored or allosteric binding pockets.     

Parallel screening and activity profiling with HIV protease inhibitor pharmacophore models

Parallel Screening has been introduced as an in silico method to predict the potential biological activities of compounds by screening them with a multitude of pharmacophore models. This study presents an early application example employing a Pipeline Pilot-based screening platform for automatic large-scale virtual activity profiling. An extensive set of HIV protease inhibitor pharmacophore models was used to screen a selection of active and inactive compounds. Furthermore, we aimed to address the usually critically eyed point, whether it is possible in a parallel screening system to differentiate between similar molecules/molecules acting on closely related proteins, and therefore we incorporated a collection of other protease inhibitors including aspartic protease inhibitors. The results of the screening experiments show a clear trend toward most extensive retrieval of known active ligands, followed by the general protease inhibitors and lowest recovery of inactive compounds.     

Structure-Based Discovery of ABCG2 Inhibitors: A Homology Protein-Based Pharmacophore Modeling and Molecular Docking Approach

ABCG2 is an ABC membrane protein reverse transport pump, which removes toxic substances such as medicines out of cells. As a result, drug bioavailability is an unexpected change and negatively influences the ADMET (absorption, distribution, metabolism, excretion, and toxicity), leading to multi-drug resistance (MDR). Currently, in spite of promising studies, screening for ABCG2 inhibitors showed modest results. The aim of this study was to search for small molecules that could inhibit the ABCG2 pump. We first used the WISS MODEL automatic server to build up ABCG2 homology protein from 655 amino acids. Pharmacophore models, which were con-structed based on strong ABCG2 inhibitors (IC₅₀ < 1 μM), consist of two hydrophobic (Hyd) groups, two hydrogen bonding acceptors (Acc2), and an aromatic or conjugated ring (Aro|PiR). Using molecular docking method, 714 substances from the DrugBank and 837 substances from the TCM with potential to inhibit the ABCG2 were obtained. These chemicals maybe favor synthesized or extracted and bioactivity testing.     

Pharmacophore modelling,virtual screening and molecular docking studies on PLD1 inhibitors

Lipid metabolism plays a significant role in influenza virus replication and subsequent infection. The regulatory mechanism governing lipid metabolism and viral replication is not properly understood to date, but both Phospholipase D (PLD1 and PLD2) activities are stimulated in viral infection. In vitro studies indicate that chemical inhibition of PLD1 delays viral entry and reduction of viral loads. The current study reports a three-dimensional pharmacophore model based on 35 known PLD1 inhibitors. A sub-set of 25 compounds was selected as the training set and the remaining 10 compounds were kept in the test set. One hundred and twelve pharmacophore models were generated; a six-featured pharmacophore model (AADDHR.57) with survival score (2.69) produced a statistically significant three-dimensional quantitative structure–activity relationship model with r² = 0.97 (internal training set), r² = 0.71 (internal test set) and Q² = 0.64. The predictive power of the pharmacophore model was validated with an external test set (r² = 0.73) and a systematic virtual screening work-flow was employed showing an enrichment factor of 23.68 at the top 2% of the dataset (active and decoys). Finally, the model was used for screening of the filtered PubChem database to fetch molecules which can be proposed as potential PLD1 inhibitors for blocking influenza infection.     

β分泌酶抑制剂的药效团模型研究

选择活性跨越0.002至25 μmol&#8226;L－1的4类共25个β分泌酶抑制剂作为训练集, 使用Catalyst软件包构建出药效团模型, 并通过对药效团的有效性分析, 筛选得到的最佳模型(correlCorrel＝0.969, Config＝16.32, Δcost＝62.422)由一个环芳香性、一个疏水中心、一个正电荷中心和一个氢键供体组成. 并用其它209个抑制剂组成测试集对模型进行验证, 结果表明该模型显示出较强的预测能力, 能够为进一步的数据库搜索, 寻找新型的β分泌酶抑制剂先导物提供依据.