Molecular modeling studies focused on 5-HT7 versus 5-HT1A selectivity. Discovery of novel phenylpyrrole derivatives with high affinity for 5-HT7 receptors

The present study discusses the well-known 5-HT7/5-HT1A selectivity issue through a new series of phenylpyrrole derivatives. The first hits emerged from a virtual screening performed on a chemolibrary. Further study led to an optimization of a preliminary 5-HT7 pharmacophore model. The importance of each pharmacophoric feature is confirmed, but these characteristics have to be coupled to geometric constraints in order to achieve a 5-HT7 selectivity. Indeed, 5-HT1A affinity probably arises from extended conformations, whereas a bent one appears to be best suited for 5-HT7 selectivity.     

On the value of homology models for virtual screening: discovering hCXCR3 antagonists by pharmacophore-based and structure-based approaches

Human chemokine receptor CXCR3 (hCXCR3) antagonists have potential therapeutic applications as antivirus, antitumor, and anti-inflammatory agents. A novel virtual screening protocol, which combines pharmacophore-based and structure-based approaches, was proposed. A three-dimensional QSAR pharmacophore model and a structure-based docking model were built to virtually screen for hCXCR3 antagonists. The hCXCR3 antagonist binding site was constructed by homology modeling and molecular dynamics (MD) simulation. By combining the structure-based and ligand-based screenings results, 95% of the compounds satisfied either pharmacophore or docking score criteria and would be chosen as hits if the union of the two searches was taken. The false negative rates were 15% for the pharmacophore model, 14% for the homology model, and 5% for the combined model. Therefore, the consistency of the pharmacophore model and the structural binding model is 219/273 = 80%. The hit rate for the virtual screening protocol is 273/286 = 95%. This work demonstrated that the quality of both the pharmacophore model and homology model can be measured by the consistency of the two models, and the false negatives in virtual screening can be reduced by combining two virtual screening approaches.     

5-HT3受体拮抗剂药效团模型的构建

以31个来源于MDDR数据库中具有抑制鼠Bezold-Jarisch反射作用的5-HT3受体拮抗剂作为训练集化合物, 构建5-HT3受体拮抗剂药效团模型. 训练集化合物具备结构多样性, 来源于相同药理模型, 活性值ED50范围为0.05~320 μg/kg i.v.. 利用Catalyst计算5-HT3受体拮抗剂的最优药效团由一个氢键受体、一个疏水基团、一个正电离子化基团、一个芳香环特征和6个排除体积组成; Fixed cost值、Null cost 值、Δcost值和Configuration cost值分别为112.6, 172.0, 59.4和7.248. 训练集化合物活性的计算值与实测值相关系数为0.9031, 偏差值为0.8976, 基于Fischer的交叉验证结果表明药效团模型具有较高的置信度, 所得药效团对训练集化合物活性值的预测结果显示有较好的预测能力, 可用于数据库搜索指导发现新的具有该活性的先导化合物, 也可用于中药或天然产物药物研究开发.     

Dihydropyrazole and dihydropyrrole structures based design of Kif15 inhibitors as novel therapeutic agents for cancer

Mitotic Kinesin motors, Eg5 and Kif15, have recently emerged as good targets for cancer as they play an inevitable role during mitosis. But, most of the Eg5 inhibitors were found ineffective when the cancer cells develop resistance to them by escalating the expression of Kif15 as alternative to Eg5. Therefore, the drugs that target Kif15 became necessary to be used either as a single or in combination with Eg5 inhibitors. The present study used 39 dihydropyrazole and 13 dihydropyrrole derivatives that were having in vitro inhibitory potential against kinesin motors to develop a common pharmacophore hypothesis AHRR and atom-based QSAR model. The model was used for virtual screening of ZINC database and the resultant hits were docked against Kif15. The four drug candidates with high docking score were examined for their activity and pharmacokinetic behaviour. Based on the results these drugs could be considered as lead candidates in further drug development for cancer.     

Molecular design based on 3D pharmacophores. Applications to 5-HT7 receptors

A definition of a pharmacophore for the 5-HT7 antagonists was carried out by searching the common chemical features of selective antagonists from the literature. A molecular design is described by analyzing the differences between this new pharmacophore and three other 3D serotonin pharmacophores previously described. This comparison led to the synthesis of a new series of potent 5-HT7 antagonists.     

Combining ethnopharmacology and virtual screening for lead structure discovery: COX-inhibitors as application example

Virtual screening of large libraries of organic compounds combined with pharmacological high throughput screening is widely used for drug discovery in the pharmaceutical industry. Our aim was to explore the efficiency of using a biased 3D database comprising secondary metabolites from antiinflammatory medicinal plants as a source for the virtual screening. For this study pharmacophore models of cyclooxygenase I and II (COX-1, COX-2), key enzymes in the inflammation process, were generated with structure-based as well as common feature based modeling, resulting in three COX hypotheses. Four different multiconfomational 3D databases limited in molecular weight between 300 and 700 Da were applied to the screening in order to compare and analyze the obtained hit rates. Two of them were created in-house (DIOS, NPD). The database DIOS consists of 2752 compounds from phytochemical reports of antiinflammatory medicinal plants described by the ethnopharmacological source 'De material medica' of Pedanius Dioscorides, whereas NPD contains almost 80,000 compounds gathered arbitrarily from natural sources. In addition, two available multiconformational 3D libraries comprising marketed and development drug substances (DWI and NCI), mainly originating from synthesis, were used for comparison. As a test of the pharmacophore models' capability in natural sources, the models were used to search for known COX inhibitory natural products. This was achieved with some exceptions, which are discussed in the paper. Depending on the hypothesis used, DWI and NCI library searches produced hit rates in the range of 6.6% to 13.7%. A slight increase of the number of molecules assessed for binding was achieved with the database of natural products (NPD). Using the biased 3D database DIOS, however, the average increase of efficiency reached 77% to 133% compared to the hit rates resulting from WDI and NCI. The statistical benefit of a combination of an ethnopharmacological approach with the potential of computer aided drug discovery by in silico screening was demonstrated exemplified on the applied targets COX-1 and COX-2.     

Rapid shape-based ligand alignment and virtual screening method based on atom/feature-pair similarities and volume overlap scoring

Shape-based methods for aligning and scoring ligands have proven to be valuable in the field of computer-aided drug design. Here, we describe a new shape-based flexible ligand superposition and virtual screening method, Phase Shape, which is shown to rapidly produce accurate 3D ligand alignments and efficiently enrich actives in virtual screening. We describe the methodology, which is based on the principle of atom distribution triplets to rapidly define trial alignments, followed by refinement of top alignments to maximize the volume overlap. The method can be run in a shape-only mode or it can include atom types or pharmacophore feature encoding, the latter consistently producing the best results for database screening. We apply Phase Shape to flexibly align molecules that bind to the same target and show that the method consistently produces correct alignments when compared with crystal structures. We then illustrate the effectiveness of the method for identifying active compounds in virtual screening of eleven diverse targets. Multiple parameters are explored, including atom typing, query structure conformation, and the database conformer generation protocol. We show that Phase Shape performs well in database screening calculations when compared with other shape-based methods using a common set of actives and decoys from the literature.     

Hierarchical virtual screening: identification of potential high-affinity and selective β(3)-adrenergic receptor agonists

The hierarchical virtual screening (HVS) study, consisting of pharmacophore modelling, docking and VS of the generated focussed virtual library, has been carried out to identify novel high-affinity and selective β(3)-adrenergic receptor (β-AR) agonists. The best pharmacophore model, comprising one H-bond donor, two hydrophobes, one positive ionizable and one negative ionizable feature, was developed based on a training set of 51 β(3)-AR agonists using the pharmacophore generation protocol implemented in Discovery Studio. The model was further validated with the test set, external set and ability of the pharmacophoric features to complement the active site amino acids of the homology modelled β(3)-AR developed using MODELLER software. The focussed virtual library was generated using the structure-based insights gained from our earlier reported comprehensive study focussing on the structural basis of β-AR subtype selectivity of representative agonists and antagonists. The HVS with the sequential use of the best pharmacophore model and homology modelled β(3)-AR in the screening of the generated focussed library has led to the identification of potential virtual leads as novel high-affinity and selective β(3)-AR agonists.     

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.     

Ligand-based molecular modeling study on a chemically diverse series of cholecystokinin-B/gastrin receptor antagonists: generation of predictive model

Pharmacophore hypotheses were developed for six structurally diverse series of cholecystokinin-B/gastrin receptor (CCK-BR) antagonists. A training set consisting of 33 compounds was carefully selected. The activity spread of the training set molecules was from 0.1 to 2100 nM. The most predictive pharmacophore model (hypothesis 1), consisting of four features, namely, two hydrogen bond donors, one hydrophobic aliphatic, and one hydrophobic aromatic feature, had a correlation (r) of 0.884 and a root-mean-square deviation of 1.1526, and the cost difference between null cost and fixed cost was 81.5 bits. The model was validated on a test set consisting of six different series of 27 structurally diverse compounds and performed well in classifying active and inactive molecules correctly. This validation approach provides confidence in the utility of the predictive pharmacophore model developed in this work as a 3D query tool in the virtual screening of drug-like molecules to retrieve new chemical entities as potent CCK-BR antagonists. The model can also be used to predict the biological activities of compounds prior to their costly and time-consuming synthesis.     

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.     

基于药效团的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小分子抑制剂。     

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

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

In Silico Screening of Natural Compounds for Candidates 5HT6 Receptor Antagonists against Alzheimer’s Disease

Alzheimer’s disease (AD), a devastating neurodegenerative disease, is the focus of pharmacological research. One of the targets that attract the most attention for the potential therapy of AD is the serotonin 5HT6 receptor, which is the receptor situated exclusively in CNS on glutamatergic and GABAergic neurons. The neurochemical impact of this receptor supports the hypothesis about its role in cognitive, learning, and memory systems, which are of critical importance for AD. Natural products are a promising source of novel bioactive compounds with potential therapeutic potential as a 5HT6 receptor antagonist in the treatment of AD dementia. The ZINC—natural product database was in silico screened in order to find the candidate antagonists of 5-HT6 receptor against AD. A virtual screening protocol that includes both short-and long-range interactions between interacting molecules was employed. First, the EIIP/AQVN filter was applied for in silico screening of the ZINC database followed by 3D QSAR and molecular docking. Ten best candidate compounds were selected from the ZINC Natural Product database as potential 5HT6 Receptor antagonists and were proposed for further evaluation. The best candidate was evaluated by molecular dynamics simulations and free energy calculations.     

SHAFTS: a hybrid approach for 3D molecular similarity calculation. 1. Method and assessment of virtual screening

We developed a novel approach called SHAFTS (SHApe-FeaTure Similarity) for 3D molecular similarity calculation and ligand-based virtual screening. SHAFTS adopts a hybrid similarity metric combined with molecular shape and colored (labeled) chemistry groups annotated by pharmacophore features for 3D similarity calculation and ranking, which is designed to integrate the strength of pharmacophore matching and volumetric overlay approaches. A feature triplet hashing method is used for fast molecular alignment poses enumeration, and the optimal superposition between the target and the query molecules can be prioritized by calculating corresponding "hybrid similarities". SHAFTS is suitable for large-scale virtual screening with single or multiple bioactive compounds as the query "templates" regardless of whether corresponding experimentally determined conformations are available. Two public test sets (DUD and Jain's sets) including active and decoy molecules from a panel of useful drug targets were adopted to evaluate the virtual screening performance. SHAFTS outperformed several other widely used virtual screening methods in terms of enrichment of known active compounds as well as novel chemotypes, thereby indicating its robustness in hit compounds identification and potential of scaffold hopping in virtual screening.     

Pharmacophore modeling as an efficient tool in the discovery of novel noncompetitive AMPA receptor antagonists

A three-dimensional pharmacophore model for the binding of noncompetitive AMPA receptor antagonists was developed in order to map common structural features of highly active compounds. This hypothesis, which consists of two hydrophobic regions, one hydrogen bond acceptor and one aromatic region, was successfully used as framework for the design of a new class of allosteric modulators containing a tetrahydroisoquinoline skeleton and for in silico screening. The promising biological results suggested that the identified molecules might be useful "lead compounds" for future drug development.     

Receptor-Based Pharmacophore Modeling in the Search for Natural Products for COVID-19 Mpro

Considering the urgency of the COVID-19 pandemic, we developed a receptor-based pharmacophore model for identifying FDA-approved drugs and hits from natural products. The COVID-19 main protease (M^pro) was selected for the development of the pharmacophore model. The model consisted of a hydrogen bond acceptor, donor, and hydrophobic features. These features demonstrated good corroboration with a previously reported model that was used to validate the present model, showing an RMSD value of 0.32. The virtual screening was carried out using the ZINC database. A set of 208,000 hits was extracted and filtered using the ligand pharmacophore mapping, applying the lead-like properties. Lipinski’s filter and the fit value filter were used to minimize hits to the top 2000. Simultaneous docking was carried out for 200 hits for natural drugs belonging to the FDA-approved drug database. The top 28 hits from these experiments, with promising predicted pharmacodynamic and pharmacokinetic properties, are reported here. To optimize these hits as M^pro inhibitors and potential treatment options for COVID-19, bench work investigations are needed.