Development and NMR validation of minimal pharmacophore hypotheses for the generation of fragment libraries enriched in heparanase inhibitors

A combined strategy based on the development of pharmacophore hypotheses and NMR approaches is reported for the identification of novel inhibitors of heparanase, a key enzyme involved in tumor metastasis through the remodeling of the subepithelial and subendothelial basement membranes, resulting in the dissemination of metastatic cancer cells. Several pharmacophore hypotheses were initially developed from the most active heparanase inhibitors known to date and, after their application to a pool of 27 known heparanase inhibitors and a database of 1,120 compounds approved by the FDA, a four-point pharmacophore model was selected as the most predictive. This model was subsequently applied to a database of 686 chemical fragments, and a subset of 100 fragments accomplishing completely or partially the four-point model was selected to perform nuclear magnetic resonance experiments to validate the hypothesis. The experimental studies confirmed the reliability of our pharmacophore model, its applicability to in silico databases in order to reduce the number of compounds to be experimentally screened, and the possibility of generating fragment libraries enriched in heparanase inhibitors. Rafael Gozalbes and Silvia Mosulén contributed equally to this work.     

Pharmacophore based 3D-QSAR modeling,virtual screening and docking for identification of potential inhibitors of β-secretase

The enzyme β-secretase-1 is responsible for the cleavage of the amyloid precursor protein, a vital step in the process of the formation of amyloid-β peptides which are known to lead to neurodegeneration causing Alzheimer’s disease. Challenges associated with toxicity and blood brain permeation inability of potential inhibitors, continue to evade a successful therapy, thus demanding the search and development of highly active and effective inhibitors. Towards these efforts, we used a ligand based pharmacophore model generation from a dataset of known inhibitors whose activities against β-secretase hovered in the nano molar range. The identified 5 feature pharmacophore model, AHHPR, was validated via three dimensional quantitative structure activity relationship as indicated by r², q² and Pearson R values of 0.9013, 0.7726 and 0.9041 respectively. For a dataset of compounds with nano molar activity, the important pharmacophore features present in the current model appear to be similar with those observed in the models resulting from much wider activity range of inhibitors. Virtual screening of the ChemBridge CNS-Set™, a database having compounds with a better suitability for central nervous system based disorders followed by docking and analysis of the ligand protein interactions resulted in the identification of eight prospective compounds with considerable diversity. The current pharmacophore model can thus be useful for the identification, design and development of potent β-secretase inhibitors which by optimization can be potential therapeutics for Alzheimer’s disease.     

表皮生长因子受体酪氨酸激酶抑制剂的药效团研究

根据一系列表皮生长因子受体酪氨酸激酶抑制剂的三维定量构效关系研究，得 到了该类抑制剂的药效团，研究结果与Novartis的药效团模型相当类似．药效团包 括一个氢键受体，一个氢键给体，一个疏水区和一个带有氯或溴原子药效团对于研 究表皮生长因子受体酪氨酸激酶抑制剂结构与活性的关系具有重要的意义．通过三 维数据库搜索可能会得到新的先导化合物．     

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.     

Development of predictive pharmacophore model for in silico screening,and 3D QSAR CoMFA and CoMSIA studies for lead optimization,for designing of potent tumor necrosis factor alpha converting enzyme inhibitors

A chemical feature-based pharmacophore model was developed for Tumor Necrosis Factor-α converting enzyme (TACE) inhibitors. A five point pharmacophore model having two hydrogen bond acceptors (A), one hydrogen bond donor (D) and two aromatic rings (R) with discrete geometries as pharmacophoric features was developed. The pharmacophore model so generated was then utilized for in silico screening of a database. The pharmacophore model so developed was validated by using four compounds having proven TACE inhibitory activity which were grafted into the database. These compounds mapped well onto the five listed pharmacophoric features. This validated pharmacophore model was also used for alignment of molecules in CoMFA and CoMSIA analysis. The contour maps of the CoMFA/CoMSIA models were utilized to provide structural insight for activity improvement of potential novel TACE inhibitors. The pharmacophore model so developed could be used for in silico screening of any commercial/in house database for identification of TACE inhibiting lead compounds, and the leads so identified could be optimized using the developed CoMSIA model. The present work highlights the tremendous potential of the two mutually complementary ligand-based drug designing techniques (i.e. pharmacophore mapping and 3D-QSAR analysis) using TACE inhibitors as prototype biologically active molecules.     

Study of Aldo-keto Reductase 1C3 Inhibitor with Novel Framework for Treating Leukaemia Based on Virtual Screening and In vitro Biological Activity Testing

Aldo-keto reductase 1C3(AKR1C3) is a potential target for the treatment of acute myeloid leukaemia and T-cell acute lymphoblastic leukaemia. In this study, pharmacophore models, molecular docking and virtual screening of target prediction were used to find a potential AKR1C3 inhibitor. Firstly, eight bacteriocin derivatives(Z1-Z8) were selected as training sets to construct 20 pharmacophore models. The best pharmacophore model MODEL_016 was obtained by Decoy test(the enrichment degree was 21.5117, and the fitting optimisation degree was 0.9668). Secondly, MODEL_016 was used for the virtual screening of ZINC database. Thirdly, the hit 83256 molecules were docked into the AKR1C3 protein. Compared to the total scores and interactions between compounds and protein, 16532 candidate compounds with higher docking scores and interactions with important residues PHE306 and TRP227 were screened. Lastly, eight compounds(A1-A8) that had good absorption, distribution, metabolism, excretion and toxicity(ADMET) properties were obtained by target prediction. Compounds A3 and A7 with high total score and good target prediction results were selected for in vitro biological activity test, whose IC₅₀ values were 268.3 and 88.94 μmol/L, respectively. The results provide an important foundation for the discovery of novel AKR1C3 inhibitors. The research methods used in this study can also provide important references for the research and development of new drugs.     

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 general strategy for creating "inactive-conformation" abl inhibitors

Kinase inhibitors that bind to the ATP cleft can be broadly classified into two groups: those that bind exclusively to the ATP site with the kinase assuming a conformation otherwise conducive to phosphotransfer (type I), and those that exploit a hydrophobic site immediately adjacent to the ATP pocket made accessible by a conformational rearrangement of the activation loop (type II). To date, all type II inhibitors were discovered by using structure-activity-guided optimization strategies. Here, we describe a general pharmacophore model of type II inhibition that enables a rational "hybrid-design" approach whereby a 3-trifluoromethylbenzamide functionality is appended to four distinct type I scaffolds in order to convert them into their corresponding type II counterparts. We demonstrate that the designed compounds function as type II inhibitors by using biochemical and cellular kinase assays and by cocrystallography with Abl.     

Ensemble pharmacophore meets ensemble docking: a novel screening strategy for the identification of RIPK1 inhibitors

Programmed cell death has been a fascinating area of research since it throws new challenges and questions in spite of the tremendous ongoing research in this field. Recently, necroptosis, a programmed form of necrotic cell death, has been implicated in many diseases including neurological disorders. Receptor interacting serine/threonine protein kinase 1 (RIPK1) is an important regulatory protein involved in the necroptosis and inhibition of this protein is essential to stop necroptotic process and eventually cell death. Current structure-based virtual screening methods involve a wide range of strategies and recently, considering the multiple protein structures for pharmacophore extraction has been emphasized as a way to improve the outcome. However, using the pharmacophoric information completely during docking is very important. Further, in such methods, using the appropriate protein structures for docking is desirable. If not, potential compound hits, obtained through pharmacophore-based screening, may not have correct ranks and scores after docking. Therefore, a comprehensive integration of different ensemble methods is essential, which may provide better virtual screening results. In this study, dual ensemble screening, a novel computational strategy was used to identify diverse and potent inhibitors against RIPK1. All the pharmacophore features present in the binding site were captured using both the apo and holo protein structures and an ensemble pharmacophore was built by combining these features. This ensemble pharmacophore was employed in pharmacophore-based screening of ZINC database. The compound hits, thus obtained, were subjected to ensemble docking. The leads acquired through docking were further validated through feature evaluation and molecular dynamics simulation.     

Chemical genomic and proteomic methods for determining kinase inhibitor selectivity

The clinical success of the Bcr-Abl tyrosine kinase inhibitor Gleevec((R)) and the recent clinical approval of a number of small molecule drugs that target protein kinases have intensified the search for novel protein kinase inhibitors. Since most small molecule kinase inhibitors target the highly conserved ATP-binding pocket of this enzyme family, the target selectivity of these molecules is a major concern. Due to the large size of the human kinome, it is a formidable challenge to determine the absolute specificity of a given protein kinase inhibitor, but recent technological developments have made substantial progress in achieving this goal. This review summarizes some of the most recent experimental techniques that have been developed for the determination of protein kinase inhibitor selectivity. Special emphasis is placed on the results of these screens and the general insights that they provide into kinase inhibitor target selectivity.     

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

BRD4靶点和多种肿瘤密切相关,是具有良好成药性的热门靶点.本文选取活性较好且结构差异较大的BRD4小分子抑制剂作为训练集分子,基于配体小分子共同特征(HipHop)方法使用Discovery Studio 3.0分子模拟软件构建了药效团.药效团通过测试集验证、ROC曲线验证(SE(sensitivity)=0.937...     

Pharmacophore fingerprint-based approach to binding site subpocket similarity and its application to bioisostere replacement

Bioisosteres have been defined as structurally different molecules or substructures that can form comparable intermolecular interactions, and therefore, fragments that bind to similar protein structures exhibit a degree of bioisosterism. We present KRIPO (Key Representation of Interaction in POckets): a new method for quantifying the similarities of binding site subpockets based on pharmacophore fingerprints. The binding site fingerprints have been optimized to improve their performance for both intra- and interprotein family comparisons. A range of attributes of the fingerprints was considered in the optimization, including the placement of pharmacophore features, whether or not the fingerprints are fuzzified, and the resolution and complexity of the pharmacophore fingerprints (2-, 3-, and 4-point fingerprints). Fuzzy 3-point pharmacophore fingerprints were found to represent the optimal balance between computational resource requirements and the identification of potential replacements. The complete PDB was converted into a database comprising almost 300?000 optimized fingerprints of local binding sites together with their associated ligand fragments. The value of the approach is demonstrated by application to two crystal structures from the Protein Data Bank: (1) a MAP kinase P38 structure in complex with a pyridinylimidazole inhibitor ( 1A9U ) and (2) a complex of thrombin with melagatran ( 1K22 ). Potentially valuable bioisosteric replacements for all subpockets of the two studied protein are identified.     

Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system

In structure-based drug discovery, researchers would like to identify all possible scaffolds for a given target. However, techniques that push the boundaries of chemical space could lead to many false positives or inhibitors that lack specificity for the target. Is it possible to broadly identify the appropriate chemical space for the inhibitors and yet maintain target specificity? To address this question, we have turned to dihydrofolate reductase (DHFR), a well-studied metabolic enzyme of pharmacological relevance. We have extended our multiple protein structure (MPS) method for receptor-based pharmacophore models to use multiple X-ray crystallographic structures. Models were created for DHFR from human and Pneumocystis carinii. These models incorporate a fair degree of protein flexibility and are highly selective for known DHFR inhibitors over drug-like non-inhibitors. Despite sharing a highly conserved active site, the pharmacophore models reflect subtle differences between the human and P. carinii forms, which identify species-specific, high-affinity inhibitors. We also use structures of DHFR from Candida albicans as a counter example. The available crystal structures show little flexibility, and the resulting models give poorer performance in identifying species-specific inhibitors. Therapeutic success for this system may depend on achieving species specificity between the related human host and these key fungal targets. The MPS technique is a promising advance for structure-based drug discovery for DHFR and other proteins of biomedical interest.     

基于虚拟筛选策略发现新型潜在JAK3小分子抑制剂

利用已知活性的分子采用基于配体的策略构建药效团模型,通过基于类药规则、药效团模型、多种精度的分子对接算法、MM/GBSA结合能预测以及ADMET筛选手段对含约250万个分子的数据库进行虚拟筛选。发现5种JAK3抑制剂的新型骨架,其中6个以1-苯基咪唑烷-2-酮为骨架的分子在与JAK3激酶的结合能以及分子的ADMET性质评价方面均表现优异,具备高JAK3抑制剂潜力,被认为是虚拟筛选的命中分子。     

基于分子对接和QSAR方法预测B-Raf II型抑制剂活性

B-Raf激酶在促分裂素原活化蛋白激酶(MAPK)信号转导通路中起着重要作用,已被确定为癌症治疗非常有吸引力的靶标.新型高效B-Raf抑制剂的开发成为癌症治疗的一个热门研究领域.本文以结构多样的B-Raf II型抑制剂为研究对象,联合应用分子对接和定量构效关系(QSAR)模型研究其定量构效关系去探讨抑制活性的起源.两个主题作为研究重点:生物活性构象和描述符.首先对分子对接方法(Glide、Gold、LigandFit、Cdocker和Libdock)进行准确性评价,后将研究的对象分子对接到B-Raf活性位点并获得生物活性构象.基于准确的对接结果,计算得到16个打分评价函数和21个能量描述符,以此构建定量构效关系模型. QSAR结果表明模型具有高度精确的拟合和强的预测能力(模型M1: r² = 0.852, r_(CV)² = 0.790, r_pre² = 0.864;模型M2: r² = 0.738, r_(CV)² = 0.812, r_pre² = 0.8605).同时探讨了对抑制活性有重要影响的描述符,结果表明打分评价函数(G_Score, -ECD, Dock_Score, PMF)与能量描述符(S(hb_ext), DE(int), Emodel)对抑制活性影响非常大.通过虚拟筛选和QSAR模型理论预测,一些新的具有潜在抑制活性的化合物作为B-Raf II型抑制剂被获得.上述信息对于进一步设计新颖高效的B-Raf II型抑制剂提供了有用的指导.     

Comparison of luminescence ADP production assay and radiometric scintillation proximity assay for Cdc7 kinase

Several assay technologies have been successfully adapted and used in HTS to screen for protein kinase inhibitors; however, emerging comparative analysis studies report very low hit overlap between the different technologies, which challenges the working assumption that hit identification is not dependent on the assay method of choice. To help address this issue, we performed two screens on the cancer target, Cdc7-Dbf4 heterodimeric protein kinase, using a direct assay detection method measuring [(33)P]-phosphate incorporation into the substrate and an indirect method measuring residual ADP production using luminescence. We conducted the two screens under similar conditions, where in one, we measured [(33)P]-phosphate incorporation using scintillation proximity assay (SPA), and in the other, we detected luminescence signal of the ATP-dependent luciferase after regenerating ATP from residual ADP (LUM). Surprisingly, little or no correlation were observed between the positives identified by the two methods; at a threshold of 30% inhibition, 25 positives were identified in the LUM screen whereas the SPA screen only identified two positives, Tannic acid and Gentian violet, with Tannic acid being common to both. We tested 20 out of the 25 positive compounds in secondary confirmatory study and confirmed 12 compounds including Tannic acid as Cdc7-Dbf4 kinase inhibitors. Gentian violet, which was only positive in the SPA screen, inhibited luminescence detection and categorized as a false positive. This report demonstrates the strong impact in detection format on the success of a screening campaign and the importance of carefully designed confirmatory assays to eliminate those compounds that target the detection part of the assay.     

Identification of novel PfDHODH inhibitors as antimalarial agents via pharmacophore-based virtual screening followed by molecular docking and in vivo antimalarial activity

Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) catalyses the fourth reaction of de novo pyrimidine biosynthesis in parasites, and represents an important target for the treatment of malaria. In this study, we describe pharmacophore-based virtual screening combined with docking study and biological evaluation as a rational strategy for identification of novel hits as antimalarial agents. Pharmacophore models were established from known PfDHODH inhibitors using the GALAHAD module with IC₅₀ values ranging from 0.033 μM to 142 μM. The best pharmacophore model consisted of three hydrogen bond acceptor, one hydrogen bond donor and one hydrophobic features. The pharmacophore models were validated through receiver operating characteristic and Günere–Henry scoring methods. The best pharmacophore model as a 3D search query was searched against the IBS database. Several compounds with different structures (scaffolds) were retrieved as hit molecules. Among these compounds, those with a Q_FIT value of more than 81 were docked in the PfDHODH enzyme to further explore the binding modes of these compounds. In silico pharmacokinetic and toxicities were predicted for the best docked molecules. Finally, the identified hits were evaluated in vivo for their antimalarial activity in a parasite inhibition assay. The hits reported here showed good potential to become novel antimalarial agents.     

Design,synthesis, biological evaluation and molecular dynamics of LAR inhibitors

In this study, firstly, the pharmacophore model was established based on LAR inhibitors. ZINC database and drug-like database were screened by Hypo-1-LAR model, and the embryonic compound ZINC71414996 was obtained. Based on this compound, we designed 9 compounds. Secondly, the synthetic route of the compound was determined by consulting Reaxys and Scifinder databases, and 9 compounds (1a-1i) were synthesized by nucleophilic substitution, Suzuki reaction and so on. Meanwhile, their structures were confirmed by ¹H NMR and ¹³C NMR. Thirdly, the Enzymatic assays was carried out, the biological evaluation of compounds 1a-1i led to the identification of a novel PTP-LAR inhibitor 1c, which displayed an IC₅₀ value of 4.8 μM. At last, molecular dynamics simulation showed that compounds could interact strongly with the key amino acids LYS1350, LYS1352, ARG1354, TYR1355, LYS1433, ASP1435, TRP1488, ASP1490, VAL1493, SER1523, ARG1528, ARG1561, GLN1570, LYS1681, thereby inhibiting the protein activity. This study constructed the pharmacophore model of LAR protein, designed small-molecule inhibitors, conducted compound synthesis and enzyme activity screening, so as to provide a basis for searching for drug-capable lead compounds.