目的:运用分子对接从生物碱数据库中筛选具有磷酸乙醇胺甲基转移酶(phosphoethanolamine methyltransferase, PEAMT)抑制活性的成分。方法:以PEAMT为药物靶标,利用相似性筛选、分子对接技术、类药性评估等多轮筛选策略,从生物碱化合物数据库中挖掘命中分子,并运用AutoDock Vina、PyMOL及LigPlus软件探索命中分子与靶点蛋白的结合方式和亲和力。结果:通过相似性虚拟筛选,共获得100个小分子,进行分子对接保留了结合能低于阈值的5个分子,其中4个分子符合Lipinski’s Rule,分别为匹格列酮、5’-腺苷酸、异鸟苷和腺苷,其中吡格列酮和5’-腺苷酸与PEAMT具有较高的结合亲和力,且作用方式与原晶体复合物相似,为最佳命中分子。结论:相似性筛选和分子对接的虚拟筛选策略可为松材线虫PEAMT抑制剂的开发提供研究方向。 Objective: To screen the components with phosphoethanolamine methyltransferase (PEAMT) inhibitory activity from alkaloid database using molecular docking. Methods: The PEAMT was taken as the drug target, and multiple rounds of virtual screening strategies, such as similarity screening, molecular docking technology and drug-like evaluation, were used to explore the hit molecules from alkaloid compound database. AutoDock Vina, PyMOL and LigPlus software were used to explore the binding mode and affinity of the hit molecules with the target protein. Result: 100 small molecules were obtained through similarity virtual screening, and 5 molecules with binding energy lower than threshold were retained by molecular docking, 4 of them were in line with Lipinski’s Rule, including pioglitazone, 5’-adenylic acid, isoguanosine and adenosine, respectively. Pioglitazone and 5’-adenylate were the best hit molecules because of high binding affinity with PEAMT and similar modes of action to that of the original crystal complex. Conclusion: The virtual screening strategy of similarity screening and molecular docking can provide research direction for the development of Bursaphelenchus xylophilus PEAMT inhibitors.
目的:运用分子对接从生物碱数据库中筛选具有磷酸乙醇胺甲基转移酶(phosphoethanolamine methyltransferase, PEAMT)抑制活性的成分。方法:以PEAMT为药物靶标,利用相似性筛选、分子对接技术、类药性评估等多轮筛选策略,从生物碱化合物数据库中挖掘命中分子,并运用AutoDock Vina、PyMOL及LigPlus软件探索命中分子与靶点蛋白的结合方式和亲和力。结果:通过相似性虚拟筛选,共获得100个小分子,进行分子对接保留了结合能低于阈值的5个分子,其中4个分子符合Lipinski’s Rule,分别为匹格列酮、5’-腺苷酸、异鸟苷和腺苷,其中吡格列酮和5’-腺苷酸与PEAMT具有较高的结合亲和力,且作用方式与原晶体复合物相似,为最佳命中分子。结论:相似性筛选和分子对接的虚拟筛选策略可为松材线虫PEAMT抑制剂的开发提供研究方向。
松材线虫,磷酸乙醇胺甲基转移酶(PEAMT),分子对接,虚拟筛选
Si Tang1, Shuangkou Chen1*, Xiaoqing Tan1, Xi Xu1, Zhenning Yang2, Fengming Ren3, Ruijia Wang1, Yaqing Zhang1
1Department of Chemistry and Chemical Engineering, Chongqing University of Science and Technology, Chongqing
2Chongqing Changfeng Chemical Industry Co. LTD, Chongqing
3Institute of Medicinal Plant Cultivation, Chongqing
Received: Sep. 22nd, 2022; accepted: Oct. 8th, 2022; published: Nov. 1st, 2022
Objective: To screen the components with phosphoethanolamine methyltransferase (PEAMT) inhibitory activity from alkaloid database using molecular docking. Methods: The PEAMT was taken as the drug target, and multiple rounds of virtual screening strategies, such as similarity screening, molecular docking technology and drug-like evaluation, were used to explore the hit molecules from alkaloid compound database. AutoDock Vina, PyMOL and LigPlus software were used to explore the binding mode and affinity of the hit molecules with the target protein. Result: 100 small molecules were obtained through similarity virtual screening, and 5 molecules with binding energy lower than threshold were retained by molecular docking, 4 of them were in line with Lipinski’s Rule, including pioglitazone, 5’-adenylic acid, isoguanosine and adenosine, respectively. Pioglitazone and 5’-adenylate were the best hit molecules because of high binding affinity with PEAMT and similar modes of action to that of the original crystal complex. Conclusion: The virtual screening strategy of similarity screening and molecular docking can provide research direction for the development of Bursaphelenchus xylophilus PEAMT inhibitors.
Keywords:Bursaphelenchus xylophilus, Phosphoethanolamine Methyltransferase (PEAMT), Molecular Docking, Virtual Screening
Copyright © 2022 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
松材线虫(Bursaphelenchus xylophilus)为美国研究者Steiner和Bührer首次在蓝变的长叶松(Pinus palustris)木材中发现 [
据报道,利用树干注射法向松树中注射鱼藤酮、盐酸左旋咪唑、磷酸左旋咪唑等植物源触杀型靶向农药,可有效杀死松材线虫,其主要通过抑制松材线虫的呼吸系统发挥灭杀作用 [
基于已知活性配体结构的相似性虚拟筛选采用Open babel 3.1.1软件的fingerprints模块;天然产物数据库中筛选的小分子的类药性评估采用DruLiTo软件 [
从世界蛋白数据库RCSB PDB中搜索磷酸乙醇胺甲基转移酶相关的靶点蛋白,物种来源选择“Bursaphelenchus xylophilus”,获取并下载PDB ID为6WLF的蛋白质结构及其结合的原配体s-腺苷同型半胱氨酸(SAH),其中受体蛋白的三维晶体结构如图1所示。为获取分子对接的对接参数及网格大小,首先利用PyRx对蛋白质结构进行删除水分子、提取配体、加氢及转换为pdbqt格式等系列处理,并运用obabel 命令将提取的原结晶配体转换为pdbqt格式,然后通过设置对接参数及网格大小将原配体对接回6WLF受体蛋白的活性位点中。为验证对接参数的合理性,通常采用对接前后两配体重叠的均方根偏差值RMSD ≤ 2.0 Å [
以6WLF原结晶配体小分子SAH作为模板分子,利用babel 命令调用Open babel的fingerprints模板分别计算生物碱化合物数据库中的412个小分子在FP2,FP3及FP4的计算方法下与模板分子间的相似性系数,并由大到小进行排序,然后选取3种分子指纹下的相似性系数排名靠前的200个分子,于Bioinformatics & Evolutionary Genomics (http://bioinformatics.psb.ugent.be/webtools/Venn)网站在线绘制venn图,获得交集分子,即筛选获得一批与模板分子高度相似的小分子。
图1. PEAMT晶体结构(PDB ID: 6WLF)
PDB ID | 重叠程度 | RMSD (Å) | 中心坐标 (x、y、z)/nm | 网格大小 (x × y × z)/nm | Ki (nmol∙L−1) |
---|---|---|---|---|---|
6WLF | 0.5788 | −3.9097、36.4618、−15.8729 | 17.2187 × 12.9292 × 16.6990 | 19.93 |
表1. 靶点蛋白6WLF的对接参数
以受体蛋白6WLF的原配体SAH为筛选中心,设置合理的对接参数,应用于分子对接以筛选活性化合物。通过重对接确定的6WLF的活性位点如图2所示。运用obabel命令将经过相似性筛选获得的分子批量进行加氢、格式转换(sdf转换为pdbqt格式)等预处理,然后在linux系统下采用AutoDock vina软件将分子与受体蛋白6WLF进行分子对接,计算结合能。以原结晶配体SAH与受体蛋白6WLF结合能的90%为阈值 [
图2. 靶点蛋白6WLF的活性位点
以6WLF原结晶配体结构为模板分子,根据FP2、FP3及FP4三种相似性方法分别计算的前200相似性系数绘制的Venn图如图3所示,发现共有100个交集分子,其中前20个分子的相似性系数列于表2中。从表中数据发现,生物碱数据库中筛选的相似性系数最高的分子为S7868 (S-Adenosyl-L-homocysteine, SAH),即受体蛋白6WLF的原配体,进一步说明相似性虚拟筛选的可靠性和合理性。另相似性系数最高的前三位分子分别为:腺苷(S1647)、阿糖腺苷(S1784)和异鸟苷(S9275)。
图3. FP2、FP3和FP4三种相似性计算方法中前200个分子的Venn图
MOL_ID | 名称 | FP2 | FP3 | FP4 |
---|---|---|---|---|
S7868 | S-Adenosyl-L-homocysteine (SAH) | 1.00 | 0.8125 | 0.8000 |
S1647 | Adenosine | 0.7086 | 0.5882 | 0.5161 |
S1784 | VIRA-A (vidarabine) | 0.7086 | 0.5882 | 0.5161 |
S9275 | Isoguanosine | 0.6899 | 0.5882 | 0.4571 |
S3610 | Cordycepin | 0.6887 | 0.5882 | 0.5161 |
S3190 | N6-methyladenosine (M6A) | 0.6887 | 0.5882 | 0.5161 |
S9351 | 2’-Deoxyadenosine monohydrate | 0.6667 | 0.5882 | 0.5161 |
S9366 | 5’-Adenylic acid | 0.6102 | 0.6250 | 0.5000 |
S2439 | Guanosine | 0.5966 | 0.5882 | 0.4571 |
S2442 | Inosine | 0.5928 | 0.5294 | 0.4054 |
S8095 | Tubercidin | 0.5926 | 0.5882 | 0.5161 |
S5678 | Trans-Zeatin-riboside | 0.5870 | 0.5882 | 0.4545 |
S9369 | Xanthosine Dihydrate | 0.5756 | 0.5294 | 0.3846 |
S9368 | ADP | 0.5714 | 0.6250 | 0.5000 |
S2053 | Cytidine | 0.5610 | 0.5882 | 0.4412 |
S9007 | 2’-Deoxyguanosine monohydrate | 0.5600 | 0.5882 | 0.4324 |
S1985 | ATP disodium | 0.5561 | 0.8125 | 0.6250 |
S3757 | 2’-Deoxyinosine | 0.5542 | 0.5294 | 0.4054 |
S9345 | Disodium 5’-Inosinate | 0.5542 | 0.5294 | 0.4054 |
S2518 | NAD+ | 0.5114 | 0.8750 | 0.5526 |
S1876 | Valaciclovir hydrochloride | 0.4972 | 0.4706 | 0.3590 |
表2. FP2、FP3和FP4三种相似性计算方法中前20个分子的相似性系数
利用Autoduck vina软件将经过相似性筛选的100个分子与受体蛋白6WLF进行分子对接,保留结合能低于阈值的5个生物碱小分子,接着进行类药性评估,最终获得4个命中分子,其结合亲和力Ki及其类药性相关参数列于表3中。其中匹格列酮(Pioglitazone)与PEAMT的结合亲和力和原结晶配体SAH的接近,为最佳命中分子。为探究所筛选配体小分子与PEAMT晶体蛋白6WLF间的作用方式,采用PyMOL和LigPlus软件分别可视了四个命中分子与6WLF的作用模式,如图4所示。由图可知,匹格列酮、5’-腺苷酸、异鸟苷和腺苷4个命中分子均位于由受体蛋白6WLF的PHE81、TYR127、MET26、TRP123、GLY61、ASP80、LEU124、ALA106、ASP105、ASN122、ARG64、PHE65、ASP57和GLY59等氨基酸残基形成的活性口袋中。根据2D相互作用图(图4左),发现4个命中分子主要与活性口袋中的氨基酸残基形成疏水相互作用和氢键相互作用,其中匹格列酮与ASP57、ARG64和PHE65氨基酸残基形成3个氢键,5’-腺苷酸与ASP105和ALA106氨基酸残基形成2个氢键,异鸟苷与ALA106形成1个氢键,而腺苷则与ALA106、ASP105氨基酸残基形成两个氢键,说明氢键相互作用和疏水相互作用极有利于配体小分子与靶点蛋白的结合,两种作用力促进了配体–蛋白复合物的稳定性。同时,从靶点蛋白与其原结晶配体SAH的作用模式图发现,SAH与ALA106、ASN122、ARG64及GLY59氨基酸残基形成了6个氢键,并与活性口袋中的氨基酸残基以疏水作用紧密结合。因此,4个命中分子均能与PEAMT受体蛋白发生强氢键作用和疏水相互作用,具有较高的结合活性,且作用方式与原结晶配体相似,可成为PEAMT抑制剂的候选化合物,而由于匹格列酮和5’-腺苷酸的结合亲和力与原结晶配体SAH接近,为PEAMT抑制剂的最佳获选化合物。
MOL_ID | 名称 | CAS号 | Ki (nmol∙L−1) | MW | LogP | HBA | HBD | nRB |
---|---|---|---|---|---|---|---|---|
S7868 | SAH | 979-92-0 | 19.93 | 384.41 | −2.30 | 9 | 5 | 7 |
S2590 | Pioglitazone | 105355-27-9 | 19.93 | 356.44 | 3.09 | 4 | 1 | 7 |
S1647 | Adenosine | 58-61-7 | 91.11 | 267.24 | −1.57 | 7 | 4 | 2 |
S9275 | Isoguanosine | 1818-71-9 | 46.37 | 283.24 | −2.09 | 7 | 5 | 2 |
S9366 | 5'-Adenylic acid | 61-19-8 | 39.16 | 347.22 | −2.48 | 10 | 5 | 4 |
表3. 4个命中分子的亲和力及其类药性相关参数
6WLF-SAH
6WLF-Pioglitazone
6WLF-5’-Adenylic acid
图4. 命中分子与PEAMT的3D (左图)和2D (右图)结合模式
本研究通过对生物碱化合物数据库进行基于已知活性配体结构的相似性筛选、分子对接筛选及类药性评估等多轮筛选,共获得与松材线虫磷酸乙醇胺甲基转移酶有较强结合活性的四个生物碱类小分子,分别为匹格列酮、5’-腺苷酸、异鸟苷和腺苷,均为具有环状结构的含氮有机化合物。此外,分子对接结果显示,筛选的4种小分子与磷酸乙醇胺甲基转移酶蛋白的原结晶配体具有相似的作用方式,且分子结构中的N原子均能与活性口袋中的氨基酸残基形成多个氢键,进一步表明首先基于结构的相似性筛选策略具有合理性。匹格列酮、5’-腺苷酸、异鸟苷和腺苷4种生物碱化合物均可通过疏水作用和氢键作用与PEAMT晶体蛋白形成稳定的复合物,具有较高的结合活性。匹格列酮和5’-腺苷酸对PEAMT的亲和力与酶蛋白的原配体接近,为PEAMT抑制剂的最佳候选化合物。但上述化合物的酶抑制活性尚需进行生物实验进一步验证。本研究可为松材线虫的磷酸乙醇胺甲基转移酶抑制剂的开发提供研究思路及物质参考。
重庆英才计划–创新创业示范团队项目(CQYC20210309793)资助;重庆英才计划–基础研究与前沿探索项目(cstc2022ycjh-bgzxm0029)资助。
唐 思,陈双扣,谭小庆,徐 曦,杨振宁,任风鸣,王瑞嘉,张亚晴. 基于分子对接的虚拟筛选发现磷酸乙醇胺甲基转移酶抑制剂Virtual Screening of Phosphoethanolamine Methyltransferase, PEAMT Inhibitors Based on Molecular Docking[J]. 药物化学, 2022, 10(04): 318-327. https://doi.org/10.12677/HJMCe.2022.104033