整合药效团及共价对接用于虚拟筛选JAK3抑制剂

大量研究表明JAK3与炎症疾病的发生、发展具有密切的关系,使得JAK3成为一个极具潜力的药物靶点。其中,JAK3共价抑制剂因其选择性高、活性强的特点受到广泛关注。但是,JAK3与其家族的其他成员同源性高,使得开发JAK3选择性抑制剂充满挑战。计算机虚拟筛选方法可以在分子水平对JAK3的结构特征进行针对性筛选,但是传统的共价对接方法效率较低、准确度欠佳,因此本文提出了一种结合药效团和共价对接的虚拟筛选策略。该联用方法从DrugBank数据库成功地筛选出已报道的JAK3临床抑制剂,表明了这种虚拟筛选不仅具有较高的效率,同时具备了较强的筛选准确性,为JAK3共价抑制剂的虚拟筛选提供一定的指导作用。     

GPR40激动剂作用方式及药效团模型研究

为阐明GPR40与其激动剂分子之间的相互作用方式,构建可用于GPR40激动剂分子先导化合物筛选的药效团模型。我们利用分子对接技术将GPR40与其激动剂分子进行对接,分析分子与受体之间相互作用的关键氨基酸和结合方式,采用药团模型法分别构建了基于受体-配体复合物(CBP)和基于激动剂分子共同特征(HipHop)的药效团模型,HipHop模型采用测试集法进行验证。结果显示GPR40与小分子相互作用的关键氨基酸主要有ARG183、TYR91、TYR2240、ARG2258、PHE142等,相互作用方式则主要为氢键、盐桥、Pi-Pi Stacking以及疏水作用,以药团模型法构建了10个HipHop模型,其中8号药效团为最优模型,可用于GPR40激动剂分子的虚拟筛选研究,这为GPR40激动剂药物分子设计奠定了理论基础。     

基于药效团模型的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个预测活性较好的化合物,可作为进一步药物研发的候选化合物.     

新型酪氨酸激酶小分子抑制的三维药效团研究

通过CATALYST软件包得到了两类HER2抑制的三维药效团模型。尽管亚苄基丙二腈化合物和3-取代吲哚啉-2-酮系列化合物具有完全不同的骨架结构,但得到的药效团却具有共同的特性,这表明当这两类抑制剂和受体发生相互作用时,采用了相似的结合模式。共同的药效团模型包括一个氢键受体,一个氢键给体,一个脂肪类疏水团以及一个芳香类疏水团。根据药效团模型,我们还进行了三维构效关系的研究,结果表明得到的药效团模型具有很好的预测能力(线性回归系数R≈0.96)。药效团模型对于研究酪氨酸激酶小分子抑制剂的结构与活性关系,以及评估和预测此类未知化合物活性具人重要的意义。     

新型酪氨酸激酶小分子抑制的三维药效团研究

通过CATALYST软件包得到了两类HER2抑制的三维药效团模型。尽管亚苄基丙二腈化合物和3-取代吲哚啉-2-酮系列化合物具有完全不同的骨架结构,但得到的药效团却具有共同的特性,这表明当这两类抑制剂和受体发生相互作用时,采用了相似的结合模式。共同的药效团模型包括一个氢键受体,一个氢键给体,一个脂肪类疏水团以及一个芳香类疏水团。根据药效团模型,我们还进行了三维构效关系的研究,结果表明得到的药效团模型具有很好的预测能力(线性回归系数R≈0.96)。药效团模型对于研究酪氨酸激酶小分子抑制剂的结构与活性关系,以及评估和预测此类未知化合物活性具人重要的意义。     

人类腺苷受体A3亚型拮抗剂的构效关系分析

构建人类腺苷受体A3亚型药效团模型和三维蛋白结构模型用于作用模式研究.以18个来源于文献具有腺苷受体A3亚型拮抗活性的化合物作为训练集,使用HypoGen方法构建药效团模型.通过同源模建和分子动力学模拟构建了人类腺苷受体A3亚型的三维蛋白模型,并利用PROCHECK方法评估该模型的合理性,对所得的结构使用分子对接程序进行作用模式分析,药效团模型和同源模建结果相互匹配较好.使用新药效团模型对MDL药物数据库(MDDR)中包含的约120000个化合物进行虚拟筛选,得到了8个候选化合物,用于进一步的生物学评价和活性测定.本工作对于人类腺苷受体A3亚型拮抗剂的设计和抗哮喘药物的研发具有一定的理论指导和应用价值.     

基于虚拟筛选和分子动力学模拟发现新型ULK1抑制剂

Unc-51样自噬激活激酶1（unc-51-like autophagy activating kinase 1,ULK1）作为自噬启动的重要调控因子,是肿瘤治疗的关键靶点之一。首先,以已知ULK1抑制剂为基础构建药效团模型,通过药效团模型筛选、分子对接以及分子力学广义波恩表面积（Molecular Mechanics/Generalized Born Surface Area,MM/GBSA）结合自由能计算等方法,对含有52万多个类药性小分子的数据库进行虚拟筛选,得到具有较高理论亲和力的化合物。随后,50ns的分子动力学模拟验证了蛋白质-配体复合物结合的稳定性,最后10ns的平均结合自由能的计算研究进一步验证了配体的结合能力。结果表明,6个化合物（F5258-0159、F3407-0428、F0529-1100、F0696-3531、F3222-5280、F6525-5596）具有骨架新颖、分子对接分数和结合自由能数值优异及与ULK1的结合状态稳定等特点,可以作为新型潜在的ULK1抑制剂用于肿瘤治疗的研究,也为新型ULK1抑制剂的设计和研发提供新的研究思路。     

基于药效团模型从中草药数据库中搜索gpIIb/IIIa受体抑制剂

选择20 个3,4-二氢-1(1H)-异喹啉酮类gpIIb/IIIa受体抑制剂作为训练集, 利用Catalyst软件包建立了gpIIb/IIIa受体抑制剂三维药效团模型. 探讨了药效团作用模式. 并通过建立的可靠性最佳的药效团模型(线性回归系数r=0.7715), 从中草药数据库中虚拟筛选了gpIIb/IIIa受体抑制剂, 通过实验活性测定得到了8个抑制ADP活化全血血小板聚集的IC50从40到100 μmol·L-1的化合物, 进一步证明了所建药效团模型的有效性.     

一种新的De Novo生物活性分子设计方法

由药效团进行虚拟活性结构生成与3D-QSAR模型相结合,筛选出有前途的结构多样性的化合物,并从中寻找活性先导化合物,是一种新的分子设计方法。采用这种方法对抗小麦赤霉病类含氟农药进行了研究,共生成了53个虚拟活性结构,通过3D-QSAR模型筛选出其中10个活性较高的结构,在活性最高的化合物基础上进行了结构修饰,得到了活性更高且毒性较低的理想化合物。研究结果表明这种方法能突破原模型化合物结构模式的局限,可以找到结构新颖的活性先导化合物,是一种非常有前途的分子设计方法,而且具有较高的筛选效率。     

基于药效团模型从中草药数据库中搜索gpⅡb/Ⅲa受体抑制剂

选择20个3,4-二氢-1(1H)-异喹啉酮类gpⅡb/Ⅲa受体抑制剂作为训练集,利用Catalyst软件包建立了gpⅡb/Ⅲa受体抑制剂三维药效团模型.探讨了药效团作用模式.并通过建立的可靠性最佳的药效团模型(线性回归系数r=0.7715),从中草药数据库中虚拟筛选了gpⅡb/Ⅲa受体抑制剂,通过实验活性测定得到了8个抑制ADP活化全血血小板聚集的IC50从40到100μmol·L-1的化合物,进一步证明了所建药效团模型的有效性.     

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

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

Pharmacophore generation,atom-based 3D-QSAR and molecular dynamics simulation analyses of pyridine-3-carboxamide-6-yl-urea analogues as potential gyrase B inhibitors

DNA gyrase subunit B (GyrB) is an attractive drug target for the development of antibacterial agents with therapeutic potential. In the present study, computational studies based on pharmacophore modelling, atom-based QSAR, molecular docking, free binding energy calculation and dynamics simulation were performed on a series of pyridine-3-carboxamide-6-yl-urea derivatives. A pharmacophore model using 49 molecules revealed structural and chemical features necessary for these molecules to inhibit GyrB. The best fitted model AADDR.13 was generated with a coefficient of determination (r²) of 0.918. This model was validated using test set molecules and had a good r² of 0.78. 3D contour maps generated by the 3D atom-based QSAR revealed the key structural features responsible for the GyrB inhibitory activity. Extra precision molecular docking showed hydrogen bond interactions with key amino acid residues of ATP-binding pocket, important for inhibitor binding. Further, binding free energy was calculated by the MM-GBSA rescoring approach to validate the binding affinity. A 10 ns MD simulation of inhibitor #47 showed the stability of the predicted binding conformations. We identified 10 virtual hits by in silico high-throughput screening. A few new molecules were also designed as potent GyrB inhibitors. The information obtained from these methodologies may be helpful to design novel inhibitors of GyrB.     

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.     

Refined pharmacophore features for virtual screening of human thromboxane A2 receptor antagonists

For a long time, the structural basis of TXA2 receptor is limited due to the lack of crystal structure information, till the release of the crystal structure of TXA2 receptor, which deepens our understanding about ligand recognition and selectivity mechanisms of this physiologically important receptor. In this research, we report the successful implementation in the discovery of an optimal pharmacophore model of human TXA2 receptor antagonists through virtual screening. Structure-based pharmacophore models were generated based on two crystal structures of human TXA2 receptor (PDB entry 6IIU and 6IIV). Docking simulation revealed interaction modes of the virtual screening hits against TXA2 receptor, which was validated through molecular dynamics simulation and binding free energy calculation. ADMET properties were also analyzed to evaluate the toxicity and physio-chemical characteristics of the hits. The research would provide valuable insight into the binding mechanisms of TXA2 receptor antagonists and thus be helpful for designing novel antagonists.     

Pharmacophore modeling coupled with molecular dynamic simulation approach to identify new leads for meprin-β metalloprotease

Human meprin beta metalloprotease, a small subgroup of the astacin family, is a potent drug target for the treatment of several disorders such as fibrosis, neurodegenerative disease in particular Alzheimer and inflammatory bowel diseases. In this study, a ligand-based pharmacophore approach has been used for the selection of potentially active compounds to understand the inhibitory activities of meprin-β by using the sulfonamide scaffold based inhibitors. Using this dataset, a pharmacophore model (Hypo1) was selected on the basis of a highest correlation coefficient (0.959), lowest total cost (105.89) and lowest root mean square deviation (1.31 Å) values. All the pharmacophore hypotheses generated from the candidate inhibitors comprised four features: two hydrogen-bond acceptor, one hydrogen-bond donor and one zinc binder feature. The best validated pharmacophore model (Hypo1) was used for virtual screening of compounds from several databases. The selective hit compounds were filtered by drug likeness property, acceptable ADMET profile, molecular docking and DFT study. Molecular dynamic simulations with the final 10 hit compounds revealed that a large number of non-covalent interactions were formed with the active site and specificity sub-pockets of the meprin beta metalloprotease. This study assists in the development of the new lead molecules as well as gives a better understanding of their interaction with meprin-β.     

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.     

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.     

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.     

Integrating pharmacophore mapping,virtual screening,density functional theory,molecular simulation towards the discovery of novel apolipoprotein (apoE ε4) inhibitors

AimAn integrated protocol of virtual screening involving molecular docking, pharmacophore probing, and simulations was established to identify small novel molecules targeting crucial residues involved in the variant apoE ε4 to mimic its behavior as apoE2 thereby eliminating the amyloid plaque accumulation and facilitating its clearance.Materials and MethodsAn excellent ligand-based and structure-based approach was made to identify common pharmacophoric features involving structure-based docking with respect to apoE ε4 leading to the development of apoE ε4 inhibitors possessing new scaffolds. An effort was made to design multiple-substituted triazine derivatives series bearing a novel scaffold. A structure-based pharmacophore mapping was developed to explore the binding sites of apoE ε4 which was taken into consideration. Subsequently, virtual screening, ADMET, DFT searches were at work to narrow down the proposed hits to be forwarded as a potential drug likes candidates. Further, the binding patterns of the best-proposed hits were studied and were forwarded for molecular dynamic simulations of 10 ns for its structural optimization.ResultsSelectivity profile for the most promising candidates was studied, revealing significantly C13 and C15 to be the most potent compounds. The proposed hits can be forwarded for further study against apoE ε4 involved in neurological disorder Alzheimer’s.     

Identification of Novel TRPC5 Inhibitors by Pharmacophore-Based and Structure-Based Approaches

Canonical transient receptor potential-5 (TRPC5), which belongs to the subfamily of transient receptor potential (TRP) channels, is a non-selective cation channel mainly expressed in the central nervous system and shows more restricted expression in the periphery. TRPC5 plays a crucial role in human physiology and pathology, for instance, anxiety, depression, epilepsy, pain, memory and chronic kidney disease (CKD). However, due to lack of the effective and selective inhibitors, its physiological and pathological mechanism remains so far unknown. It is therefore pivotal to identify potential TRPC5 inhibitors. We have applied ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS) methods. The pharmacophore models of TRPC5 antagonists generated by using the HypoGen and HipHop algorithms were used as a query model for the screening of potential inhibitors against the Specs database. The resultant hits from LBVS were further screened by SBVS. SBVS was carried out based on the homology model generation of human TRPC5, binding site identification, molecular dynamics optimization and molecular docking studies. In our systematic screening approaches, we have identified 7 hits compounds with comparable dock score after Lipinski and Veber rules, ADMET, PAINS analysis, cluster analysis, and similarity analysis. In conclusion, the current research provides novel backbones for the new-generation of TRPC5 inhibitors.