网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

药对“虎杖-土茯苓”治疗痛风性关节炎的网络药理学研究  PDF

  • 刘强 1
  • 张开伟 2
1. 贵州中医药大学研究生院(贵州 贵阳 550025); 2. 贵州中医药大学第一附属医院(贵州 贵阳 550001)

最近更新:2022-01-10

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

目的

采用网络药理学、分子对接技术预测药对“虎杖-土茯苓”治疗痛风性关节炎(gouty arthritis,GA)的潜在机理。

方法

利用TCMSP平台检索收集药对化学成分、对应靶点;借助五大疾病数据库获取GA已知靶点;利用Venny 2.1.0得到药对与GA两者间的交集靶点;运用SRTING数据库、CytoHubba插件构建交集靶点PPI网络;采用DAVID平台对交集靶点进行GO、KEGG通路富集分析;最后将筛选出的核心成分、关键靶点进行分子对接。

结果

共获得23个化学成分,199个靶点,其中治疗GA的潜在靶点有76个。预测得知核心成分为槲皮素、木犀草素、β-谷甾醇,关键靶点为JUN、TNF、RELA、IL6。GO功能富集分析显示,药对治疗GA的生物功能:287条(生物过程:172个,细胞组分:42个,分子功能:73个);KEGG通路富集共筛选出信号通路:122条。分子对接结果显示,槲皮素、木犀草素、β-谷甾醇与JUN、TNF有良好的结合活力。

结论

药对“虎杖-土茯苓”通过槲皮素、木犀草素等成分,作用于JUN、TNF等靶点,癌症通路、TNF信号通路等通路参与下,发挥抗炎、抗凋亡等治疗GA的作用。

痛风性关节炎(gouty arthritis,GA)是常见的关节炎性病变,与嘌呤异常代谢/尿酸排泄障碍有关,以受累关节发生急性晶体滑膜炎为典型特[

1]。流行病学调查显示,全球GA的患病率呈逐年增长的趋[2],目前临床一线药物长期服用均有不同程度的不良反应,如不加以控制,则会影响患者的生活质量,甚至致残。为此,“治未病、整体观、辨证审因”的中医理念则显得格外重要。研究显示中医药在治疗GA方面具有效果独特、成本低且安全性高等多项优势,为GA防治做出独特贡[3]

虎杖源自蓼科植物虎杖的根茎及根,味苦、性微冷,入肝、胆、肺经,以解毒镇痛、祛风利湿为功效,现代研究表明虎杖治疗GA具有确切的疗[

4];土茯苓源自百合科植物土茯苓的根茎,味甘淡、性平,入肝、胃经,以除湿止痛、解毒利关节见长,善治筋骨挛痛,为临床防治GA的常用中[5]。一项中医药防治GA的临床用药规律研究显示“虎杖-土茯苓”两味药物为治疗GA的常见核心药物组合之[6],但该药对潜在作用机理及核心成分尚未完全清楚。因此,本研究选取药对“虎杖-土茯苓”,基于网络药理学多中心-多途径的思路,预测该药对治疗GA的潜在作用机理,以期为该药对治疗GA的后期实验研究提出科学理论支撑,丰富中医药治疗GA的现代药理内涵。

1 资料及方法

1.1 数据库及软件工具

表1所示。

表1 数据库及软件工具
数据库名称网址数据分析软件
TCMSP https://lsp.nwu.edu.cn/tcmsp.php Cytoscape 3.7.2
Uniport https://www.uniprot.org CytoHubba插件
GeneCards https://www.genecards.org/ R语言 4.0.2软件
OMIM https://www.omim.org Autodock 4.2.6 分子对接软件
DisGeNET https://www.disgenet.org Pymol 2.0.1 可视化软件
Drugbank https://go.drugbank.com/
CTD https://ctdbase.org/
STRING https://string-db.org
DAVID https://david.ncifcrf.gov/
Venny 2.1.0 https://bioinfogp.cnb.csic.es /tools /venny /index.html
PDB https: //www.rcsb.org/

1.2 化学成分及靶基因的获取

在中药药理学分析平台(TCMSP)中,分别以虎杖、土茯苓为检索词,将挑选标准设置为OB(%)≥30,DL≥0.18;获取满足标准的化学成分及其对应的蛋白靶点。在Uniprot数据库中查询并获得蛋白靶点的对应靶基因名称,在Cytoscape 3.7.2中构建“药物-活性成分-靶点”图。

1.3 GA相关靶标预测

选择GeneCards、OMIM、DisGeNET、Drugbank、CTD数据库,以“gouty arthritis”或“GA”为关键词,检索GA的相关靶标基因,整理并删除重复项。通过Venny 2.1.0绘制药对靶点与GA靶点韦恩图,取交集靶点,并将其视为该药对治疗GA的潜在作用靶点。

1.4 构建“化学成分-交集靶点”网络图

整合1.3中交集靶点所映射出的化学成分,选用Cytoscape 3.7.2构建“化学成分-交集靶点”关系图,以便阐明“虎杖-土茯苓”化学成分与GA靶点间的联系。使用Net-work Analyzer功能,以度值(Degree)为衡量标准,评估关系网络中治疗GA的核心化学成分。

1.5 蛋白质相互网络(PPI)的构建

在STRING数据库中输入交集靶点,物种限定为“人类”,筛选阈值条件为>0.900,创建蛋白互作网络(PPI)。将tsv格式数据导入Cytoscape 3.7.2,利用CytoHubba插件分析数据。

1.6 基因富集分析

DAVID平台对交集靶点进行GO、KEGG通路富集分析(P<0.05),将排名前10的GO、排名前20的KEGG导入R语言4.0.2可视化;筛选排名前20通路所对应核心靶点,依托Cytoscape 3.7.2绘制“核心靶点-通路”图。

1.7 分子对接

运用分子对接技术预测分析核心化合物作用于核心靶点的结合活力。在TCMSP平台下载核心化合物mol 2格式3D结构;PDB数据库下载关键靶点蛋白结构,Autodock 4.2.6对两者进行分子对接;利用Pymol 2.0.1进行可视化处理。

2 结果

2.1 “虎杖-土茯苓”活性成分及靶基因

通过TCMSP平台挑选出“虎杖”治疗GA的有效活性成分10个,相应靶点168个;“土茯苓”治疗GA的有效成分15个,相应靶点177个,两者共同拥有有效活性成分为槲皮素和β-谷甾醇,去重后靶基因共199个(见表2)。利用Cytoscape 3.7.2构建包括220个节点、534条边的“药物-活性成分-靶点”网络图(见图1)。

表2 “虎杖-土茯苓”潜在活性化合物
来源Mol IDMolecule NameOB (%)DLDegree
土茯苓 MOL013118 Neoastilbin 40.54 0.74 1
土茯苓 MOL013119 Enhydrin 40.56 0.74 1
土茯苓 MOL004567 isoengelitin 34.65 0.7 1
土茯苓 MOL013129 (2R,3R)-2-(3,5-dihydroxyphenyl)-3,5,7-trihydroxychroman-4-one 63.17 0.27 3
土茯苓 MOL001736 (-)-taxifolin 60.51 0.27 3
土茯苓 MOL000359 sitosterol 36.91 0.75 3
土茯苓 MOL004575 astilbin 36.46 0.74 3
土茯苓 MOL004580 cis-Dihydroquercetin 66.44 0.27 4
虎杖 MOL013288 Picralinal 58.01 0.75 4
土茯苓 MOL013117 4,7-Dihydroxy-5-methoxyl-6-methyl-8-formyl-flavan 37.03 0.28 5
虎杖 MOL002268 rhein 47.07 0.28 5
虎杖 MOL000492 (+)-catechin 54.83 0.24 7
土茯苓 MOL004576 taxifolin 57.84 0.27 9
虎杖 MOL013281 6,8-Dihydroxy-7-methoxyxanthone 35.83 0.21 10
土茯苓 MOL000546 diosgenin 80.88 0.81 16
土茯苓 MOL000449 Stigmasterol 43.83 0.76 26
虎杖 MOL013287 Physovenine 106.21 0.19 30
土茯苓 MOL004328 naringenin 59.29 0.21 32
虎杖 MOL000006 luteolin 36.16 0.25 51
两者共有 MOL000358 beta-sitosterol 36.91 0.75 58
两者共有 MOL000098 quercetin 46.43 0.28 262
虎杖 MOL002259 Physciondiglucoside 41.65 0.63
虎杖 MOL002280 Torachrysone-8-O-beta-D-(6'-oxayl)-glucoside 43.02 0.74

图1 “虎杖-土茯苓”药物-活性成分-靶点网络图

2.2 GA相关靶点

五大疾病数据库共筛选出GA已知疾病靶点823个,取交集后获得交集靶点76个(见图2),这些交集靶点可能是“虎杖-土茯苓”治疗GA的潜在靶点。构建“化学成分-交集靶点”网络图(见图3),左侧菱形代表交集靶点,右侧圆形代表对应化合物,本研究将“化学成分-交集靶点”网络图中Deree值排名前三的化合物视为核心化合物,分别是槲皮素(MOL000098度值125)、木犀草素(MOL000006度值27)、β-谷甾醇(MOL000358度值22),研究结果表明了药对治疗GA“多成分-多靶点”的特点。

图2 “虎杖-土茯苓”与GA的交集靶点韦恩图

图3 “虎杖-土茯苓”化学成分-交集靶点”网络图

2.3 靶点PPI网络分析

76个交集靶点输入STRING平台,获取tsv格式的蛋白互作网络,利用Cytoscape 3.7.2的CytoHubba插件获得核心靶点,Hubba nodes设置为Degree值排名前10的节点,得到核心靶点:JUN(度值27)、TNF(度值27)、RELA(度值24)、IL6(度值22)、MAPK1(度值19)、MAPK3(度值19)、MAPK14(度值17)、IL1B(度值17)、FOS(度值15)、CXCL8(度值15)(见图4)。图中颜色深浅与节点度值呈正比,中心位置节点代表度值前10的核心基因。

图4 交集基因蛋白互作关系图

2.4 GO、KEGG富集分析结果

DAVID平台对交集靶点进行GO、KEGG分析,挑选出富集度靠前的分子功能(MF)、生物过程(BP)、细胞组分(CC)以及信号通路,使用R语言4.0.2可视化(见图5)。MF共73个,靠前的为蛋白结合(61个靶点/80.2%),相同蛋白质结合(16个靶点/21.1%),酶结合(14个靶点/18.4%)等分子功能;BP共172个,靠前的为RNA聚合酶Ⅱ启动子转录的正调控(27个靶点/35.5%),对药物的反应(23个靶点/30.2%),正调控转录-DNA模板化(20个靶点/26.3%)等生物过程;CC共42个,靠前的为细胞溶质(36个靶点/47.4%),细胞质(33个靶点/43.4%),细胞核(33个靶点/43.4%)等细胞组分。“虎杖-土茯苓”治疗GA的主要信号通路包括122条,如癌症通路(31个靶点/40.8%),TNF信号通路(22个靶点/28.9%),乙肝(22个靶点/28.9%),PI3K-Akt信号通路(22个靶点/28.9%),MAPK信号通路(17个靶点/22.4%)等通路。由此可推测出,“虎杖-土茯苓”治疗GA可能是通过上述通路完成。基于显著富集通路与调控基因的关联性,构建“核心靶点-通路”图(见图6),这些显著富集于核心通路的基因可能是药对治疗GA的关键作用点。

图5 GO、KEGG富集分析图

图6 “核心靶点-信号通路”相互作用网络图

2.5 分子对接结果

主要核心化合物为槲皮素、木犀草素、β-谷甾醇与关键靶基因JUN(5t01)、TNF(4g3y)进行分子对接。结合能越低,两者之间的作用越好。结果显示,所有对接结合能均<0 KJ/moL,表示配体与受体均可自发结合。其中,β-谷甾醇与JUN(23.3)结合最好,其它依次为木犀草素与JUN(19.4),β-谷甾醇与TNF(19.0),木犀草素与TNF(17.4),槲皮素与JUN(17.4),槲皮素与TNF(13.6)(见图7)。

图7 分子对接结果

3 讨论

本研究以筛选药对“虎杖-土茯苓”治疗GA的潜在核心成分为切入点,构建“药物-化合物-靶点-疾病”网络图,探寻潜在作用机理。研究结果显示“虎杖-土茯苓”有效化合物共23个,治疗GA的药效化合物有17个,将槲皮素、木犀草素及β-谷甾醇视为治疗GA核心化合物。槲皮素已被证明具有良好的消炎活性且作用长[

7]。槲皮素对COX-2、NF-κB、AP-1、MAPK等多种炎性因子具有抑制表达的作[8]。既往研究已表明槲皮素在治疗GA方面药效显[9],其表现出很强的抗炎镇痛作[10]。槲皮素的肝肾保护作用可能对阻止GA的肝肾损害进程有[11]。木犀草素与槲皮素均为黄酮类代表化合物,具有理想的抗炎效[12]。木犀草素可以通过调控炎性因子控制GA炎性反[13]。β-谷甾醇属于甾醇类化合物,具有理想的抗氧化、抗炎作[14]。由此可见,“虎杖-土茯苓”核心化合物能从多个方面发挥治疗GA作用。

本研究共获得76个“虎杖-土茯苓”治疗GA的潜在作用靶点,根据PPI网络分析结果提示JUN、TNF、RELA、IL6为关键靶点。JUN作为MAPK通路重要成员之一,与炎症及细胞凋亡关系紧密。研究表明当GA发作时ERK1/2通路可以降低JUN激活缓解MSU晶体引起的炎[

15]。TNF与IL6均是炎症细胞因子家族中关键成员,常被作为衡量GA炎性程度的检测指标,近期研究发现IL6的表达水平与GA的疼痛程度具有较好的相关[16]。既往实验发现TNF与体内尿酸水平、GA的疾病进展及预后关系紧[17]。RELA是炎性通路NF-kB家族关键成员,具有调节炎症、中性粒细胞凋亡的重要作[18]。本研发现的靶点可以考虑为未来研究“虎杖-土茯苓”治疗GA机制的主要关注靶点。

通过对“虎杖-土茯苓”治疗GA的GO及KEGG分析结果可知,癌症通路、TNF信号通路、乙肝、PI3K-Akt信号通路为药对治疗GA的关键信号通路。癌症通路是条与多条信号通路均存在联系的代谢通路,同时本研究中排名前20通路中的PI3K-Akt信号通路、MAPK通路、HIF-1信号通路与之皆有关联,由此可说明此通路可能是“虎杖-土茯苓”治疗GA机制中起到枢纽作用的关键通路。TNF信号通路与细胞凋亡、存活,以及免疫和炎症相关,TNF-α参与关节滑膜炎症反应,可以控制GA受累关节的急性炎[

19]。乙肝病毒可干扰细胞的凋[20]。PI3K-Akt信号通路可被多种类型的细胞刺激激活,具有抗细胞凋亡、增殖及分化功[21]。研究发现通过抑制该通路可以改善大鼠关节软骨细胞破坏,降低关节炎[22]。提示“虎杖-土茯苓”可依托上述关键通路抗击GA关节炎症的作用。

综上所述,本研究利用网络药理学和分子对接技术初步预测了“虎杖-土茯苓”治疗GA的机理。从结果上可知,“虎杖-土茯苓”通过槲皮素、木犀草素等成分,作用于JUN、TNF等靶点,在癌症通路、TNF信号通路等通路参与下,发挥抗炎、抗凋亡等治疗GA的作用,为后期实验进一步研究打下基础。

参考文献

1

ABHISHEK A,RODDY E,DOHERTY M.Gout-a guide for the general and acute physicians[J].Clinical medicine (London, England)201717(1):54-59. [百度学术

2

SINGH JA,GAFFO A.Gout epidemiology and comorbidities[J].Seminars in arthritis and rheumatism202050(3):s11-s6. [百度学术

3

赵欣禹,李可大.痛风性关节炎急性期的中医药治疗进展[J].实用中医内科杂志2021,35(7):11-13. [百度学术

4

谢招虎,解 静,李兆福,.虎杖治疗痛风的作用机制及临床应用研究进展[J].辽宁中医杂志.2019;46(5):1111-1113. [百度学术

5

王 特,张晓宇,张 薇,.土茯苓防治痛风的作用机制研究进展[J].辽宁中医杂志.202148(1):215-217. [百度学术

6

孙宇洁,付书璠,李 慧,.近30年中医药治疗痛风临床用药规律数据挖掘研究[J].中药药理与临床202036(4):208-213. [百度学术

7

BATIHA GE,BESHBISHY AM,IKRAM M,et al.The Pharmacological Activity,Biochemical Properties,and Pharmacokinetics of the Major Natural Polyphenolic Flavonoid:Quercetin[J].Foods(Basel,Switzerland)20209(3):111. [百度学术

8

GARCÍA-MEDIAVILLA V,CRESPO I,COLLADO PS,et al.The anti-inflammatory flavones quercetin and kaempferol cause inhibition of inducible nitric oxide synthase,cyclooxygenase-2 and reactive C-protein,and down-regulation of the nuclear factor kappa B pathway in Chang Liver cells[J].European journal of pharmacology2007557(2-3):221-229. [百度学术

9

谢凯莉,李昭华,董先智,.槲皮素抑制黄嘌呤氧化酶活性的研究进展[J].时珍国医国药201930(9):2223-2225. [百度学术

10

RUIZ-MIYAZAWA KW,STAURENGO-FERRARI L,MIZOKAMI SS,et al.Quercetin inhibits gout arthritis in mice:induction of an opioid-dependent regulation of inflammasome[J].Inflammopharmacology2017,5(15):555-570. [百度学术

11

LEKIĆ N,CANOVÁ NK,HOŘÍNEK A,et al.The involvement of heme oxygenase 1 but not nitric oxide synthase 2 in a hepatoprotective action of quercetin in lipopolysaccharide-induced hepatotoxicity of D-galactosamine sensitized rats[J].Fitoterapia20136(87):20-26. [百度学术

12

周霄楠,韩 超,宋鹏琰,.木犀草素和槲皮素体外抗炎作用研究[J].动物医学进展201738(10):56-61. [百度学术

13

沈瑞明,马丽辉, 郑颜萍.木犀草素通过TLR/MyD88/NF-κB通路参与急性痛风性关节炎大鼠的抗炎作用[J].中南大学学报(医学版)202045(2):115-122. [百度学术

14

刘威良,姬 昱,黄艾祥.β-谷甾醇的研究及开发进展[J].农产品加工2019(1):77,7982. [百度学术

15

JARAMILLO M,GODBOUT M, NACCACHE PH,et al.Signaling events involved in macrophage chemokine expression in response to monosodium urate crystals[J].The Journal of biological chemistry2004279(50):52797-805. [百度学术

16

王思齐,虞 斌, 张晓海, .急性痛风患者疼痛与炎症因子相关性分析[J].医学理论与实践202134(3):371-373. [百度学术

17

赖爱云,徐 健, 陶 丽.TNF-α在痛风性关节炎患者炎性反应中的变化及意义[J].河北医药201941(3):388-391. [百度学术

18

耿云峰,杜鸿斌,刘琳琳,.NF-κB家族成员RelA的翻译后修饰及其生理病理作用的研究进展[J].生命科学202032(5):431,438. [百度学术

19

URECH DM,FEIGE U,EWERT S,et al.Anti-inflammatory and cartilage-protecting effects of an intra-articularly injected anti-TNF{alpha}single-chain Fv antibody (ESBA105) designed for local therapeutic use[J].Annals of the rheumatic diseases201069(2):443-449. [百度学术

20

LIN S,ZHANG YJ.Interference of Apoptosis by Hepatitis B Virus[J].Viruses20179(8):230. [百度学术

21

XIE Y,SHI X,SHENG K,et al.PI3K/Akt signaling transduction pathway,erythropoiesis and glycolysis in hypoxia (Review)[J].Molecular medicine reports201919(2):783-791. [百度学术

22

XUE JF, SHI ZM, ZOU J,et al.Inhibition of PI3K/AKT/mTOR signaling pathway promotes autophagy of articular chondrocytes and attenuates inflammatory response in rats with osteoarthritis[J].Biomedicine & pharmacotherapy=Biomedecine & pharmacotherapie201789:1252-1261. [百度学术