留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

应中央军委要求,2022年9月起,《药学实践杂志》将更名为《药学实践与服务》,双月刊,正文96页;2023年1月起,拟出版月刊,正文64页,数据库收录情况与原《药学实践杂志》相同。欢迎作者踊跃投稿!

基于网络药理学研究柴胡达胸合剂治疗新型冠状病毒肺炎的作用机制

肖志军 刘萃萃 陆赛花 蔡建 徐峰

肖志军, 刘萃萃, 陆赛花, 蔡建, 徐峰. 基于网络药理学研究柴胡达胸合剂治疗新型冠状病毒肺炎的作用机制[J]. 药学实践与服务, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023
引用本文: 肖志军, 刘萃萃, 陆赛花, 蔡建, 徐峰. 基于网络药理学研究柴胡达胸合剂治疗新型冠状病毒肺炎的作用机制[J]. 药学实践与服务, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023
XIAO Zhijun, LIU Cuicui, LU Saihua, CAI Jian, XU Feng. The mechanism study on Chaihudaxiong mixture in the treatment of coronavirus disease 2019 with network pharmacology approach[J]. Journal of Pharmaceutical Practice and Service, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023
Citation: XIAO Zhijun, LIU Cuicui, LU Saihua, CAI Jian, XU Feng. The mechanism study on Chaihudaxiong mixture in the treatment of coronavirus disease 2019 with network pharmacology approach[J]. Journal of Pharmaceutical Practice and Service, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023

基于网络药理学研究柴胡达胸合剂治疗新型冠状病毒肺炎的作用机制

doi: 10.12206/j.issn.1006-0111.202004023
基金项目: 上海市卫生健康委员会科研课题(20154Y0123),上海市卫生健康委员会临床药学重点专科建设项目(2018-2020)
详细信息
    作者简介:

    肖志军,硕士,主管药师,Email:[email protected]

    通讯作者: 徐 峰,博士,主任药师,研究方向:临床合理用药,Email:[email protected]
  • 中图分类号: R285.5

The mechanism study on Chaihudaxiong mixture in the treatment of coronavirus disease 2019 with network pharmacology approach

图(6) / 表(2)
计量
  • 文章访问数:  6896
  • HTML全文浏览量:  2913
  • PDF下载量:  45
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-05
  • 修回日期:  2020-05-08
  • 网络出版日期:  2020-07-27
  • 刊出日期:  2020-07-25

基于网络药理学研究柴胡达胸合剂治疗新型冠状病毒肺炎的作用机制

doi: 10.12206/j.issn.1006-0111.202004023
    基金项目:  上海市卫生健康委员会科研课题(20154Y0123),上海市卫生健康委员会临床药学重点专科建设项目(2018-2020)
    作者简介:

    肖志军,硕士,主管药师,Email:[email protected]

    通讯作者: 徐 峰,博士,主任药师,研究方向:临床合理用药,Email:[email protected]
  • 中图分类号: R285.5

摘要:   目的  基于网络药理学方法探讨柴胡达胸合剂治疗新型冠状病毒肺炎(COVID-19)的潜在药理作用机制。  方法  利用中药系统药理学数据库和分析平台(TCMSP)筛选柴胡达胸合剂活性成分和对应的作用靶点,并通过Uniprot数据库标准化靶点名称;在GeneCards和OMIM数据库检索冠状病毒相关基因,并与柴胡达胸合剂作用靶点取交集,筛选出治疗作用靶点;利用Cytoscape 3.7.2软件,构建和分析“药材-活性成分-靶点”网络图;通过String平台分析靶点蛋白相互作用,并使用R软件的相关包进行GO基因注释和KEGG信号通路分析。  结果  筛选出165个有效成分和51个作用靶点,进一步分析发现主要活性成分为β-谷甾醇和11个黄酮类化合物,核心作用靶点为CASP3、MAPK3、IL-6、MAPK8、IL-10、CXCL8、MAPK1、IL-1B等。GO基因注释得到GO条目共1 722个(P<0.05),其中生物学过程条目1 612个,细胞组成条目30个,分子功能条目80个。KEGG信号通路筛选出信号通路156条(P<0.05),富集基因较多的信号通路为糖尿病并发症中的AGE-RAGE信号通路、甲型流感、IL-17信号通路、TNF信号通路和乙型肝炎。  结论  该研究初步揭示了柴胡达胸合剂多成分、多靶点、多通路对COVID-19发挥治疗作用的特点,为进一步阐明柴胡达胸合剂治疗COVID-19的药理作用机制提供理论依据。

English Abstract

肖志军, 刘萃萃, 陆赛花, 蔡建, 徐峰. 基于网络药理学研究柴胡达胸合剂治疗新型冠状病毒肺炎的作用机制[J]. 药学实践与服务, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023
引用本文: 肖志军, 刘萃萃, 陆赛花, 蔡建, 徐峰. 基于网络药理学研究柴胡达胸合剂治疗新型冠状病毒肺炎的作用机制[J]. 药学实践与服务, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023
XIAO Zhijun, LIU Cuicui, LU Saihua, CAI Jian, XU Feng. The mechanism study on Chaihudaxiong mixture in the treatment of coronavirus disease 2019 with network pharmacology approach[J]. Journal of Pharmaceutical Practice and Service, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023
Citation: XIAO Zhijun, LIU Cuicui, LU Saihua, CAI Jian, XU Feng. The mechanism study on Chaihudaxiong mixture in the treatment of coronavirus disease 2019 with network pharmacology approach[J]. Journal of Pharmaceutical Practice and Service, 2020, 38(4): 289-295. doi: 10.12206/j.issn.1006-0111.202004023
  • 新型冠状病毒肺炎(coronavirus disease 2019,COVID-19),是指2019年始发、由严重急性呼吸道综合征冠状病毒2型(severe acute respiratory syn-drome coronavirus 2,SARS-CoV-2)引起的肺炎。截至2020年4月5日,全球共确诊COVID-19患者1 093 349例,死亡58 620例[1]。目前尚无针对COVID-19的特异性治疗药物,一些化学药物包括氯喹/羟氯喹、洛匹那韦/利托那韦、瑞德西韦等正在临床开展随机对照研究。临床实践表明,清肺排毒汤和连花清瘟胶囊等多种中药方剂和制剂对COVID-19有良好的治疗效果。据国家卫生健康委员会报道,在我国确诊的COVID-19病例中,有74 187人使用了中医药,占91.5%;中医药能够缓解症状,减少轻型、普通型疾病向重型发展,提高治愈率、降低病死率,总有效率达90%以上[2]

    柴胡达胸合剂,曾用名为“强力肺炎1号”,是国医大师梅国强教授为痰热壅肺证COVID-19患者制定的中药处方[3]。柴胡达胸合剂由小柴胡汤、小陷胸汤、达原饮、止嗽散共同组方,包含柴胡、黄芩、法半夏、全瓜蒌、黄连、枳实、甘草、浙贝母、桔梗、百部、前胡、紫苑、款冬花、槟榔、草果、藿香、佩兰、虎杖共十八味中药。由于临床使用疗效显著,2020年2月23日,湖北省药品监督管理局下发制剂备案批件,包括柴胡达胸合剂在内的2个由湖北省中医院研制的医院制剂获批用于防治COVID-19[4]

    网络药理学是基于系统生物学和多向药理学技术和方法,通过构建“药物-基因-疾病”网络,分析药物在网络中与特定节点相互作用的关系,从整体角度探索药物与机体相互作用的一门学科[5]。2007年,Hopkins首次在Nature Biotechnology杂志上发表述评,提出网络药理学这一概念,并认为其为发现新药的新范式[6]。中药通过多成分、多靶点、多通路对疾病产生治疗作用,利用网络药理学方法,可系统阐明中药治疗疾病的药理作用机制[7]

    因此,本研究运用网络药理学方法,筛选柴胡达胸合剂治疗COVID-19的活性成分和作用靶点,构建“药材-活性成分-靶点”网络图,然后对靶点进行蛋白质-蛋白质相互作用(protein-protein interac-tion, PPI)、GO基因注释和KEGG信号通路分析,为进一步阐明柴胡达胸合剂治疗COVID-19的药理作用机制提供理论依据。

    • 通过《中华人民共和国药典》(2015年版)手工检索柴胡达胸合剂的十八味中药材的性味归经,利用Cytoscape 3.7.2软件制作“药材-性味归经”网络图。在中药系统药理学数据库和分析平台(TCMSP)以柴胡达胸合剂中的十八味中药材名为关键词检索得到所有中药的化学成分。生物利用度(oral bioavailability,OB)和半衰期(half life,HL)是影响药动学的重要参数,而类药性(drug-likeness,DL)可以反映化合物的理化性质与已上市的药物是否类似。根据TCMSP数据库推荐的筛选标准,本研究以OB≥30%,DL≥0.18及HL≥4 h为条件,筛选得到每味中药材的活性成分。同时,通过 TCMSP 数据库查找各活性成分的作用靶点,归纳整理后利用Uniprot 数据库(https://www.uniprot.org/)标准化靶点名称。

    • 以coronavirus为关键词,检索GeneCards(https://www.genecards.org/)和OMIM(https://omim.org/)数据库获得COVID-19潜在相关基因;利用Venn图在线绘制工具(http://bioinformatics.psb.ugent.be/webtools/Venn/),将柴胡达胸合剂活性成分的作用靶点和COVID-19相关基因取交集,得到柴胡达胸合剂治疗 COVID-19作用靶点;最后,将得到的药材-活性成分-靶点关联性文本导入Cytoscape 3.7.2软件,构建并分析“药材-活性成分-靶点”网络。

    • 将柴胡达胸合剂治疗 COVID-19作用靶点导入STRING 蛋白相互作用数据库(https://string-db.org/),物种选定为Homo sapiens,获得PPI信息并导入Cytoscape 3.7.2软件,利用NetworkAnalyzer功能,分析网络中每个靶点的度值(dgree value),使用R软件Graphics包,绘制条形图展示度值排名前20位的靶点。

    • 利用R软件的org.Hs.eg.db和clusterProfiler包,对柴胡达胸合剂治疗COVID-19的作用靶点进行GO基因注释和KEGG通路分析,以P<0.05进行筛选,得到柴胡达胸合剂对COVID-19发挥治疗作用参与的生物学过程和信号通路,并绘制气泡图展示结果。

    • 柴胡达胸合剂“药材-性味归经”网络如图1所示,图中节点的大小代表该节点在网络中的度值。由图可知,度值最大的性味归经分别为寒(度值=8)、苦(度值=13)和肺经(度值=14)。

      图  1  柴胡达胸合剂的“药材-性味归经”网络图

    • 通过TCMSP平台检索到柴胡达胸合剂中的十八味药材共包含1 977个化合物。以OB≥30%,DL≥0.18及HL≥4 h为条件筛选并去重后,共得到221个活性成分。在TCMSP检索活性成分的作用靶点,并到Uniprot数据库中查找其标准名称,最后共得到259个作用靶点。柴胡达胸合剂的“中药-化合物-活性成分-靶点”信息,结果见表1

      表 1  柴胡达胸合剂的“中药-化合物-活性成分-靶点”信息表

      中药名称化合物(个)活性成分(个)靶点(个)
      柴胡34914165
      黄芩1433295
      法半夏1161170
      全瓜蒌8078
      黄连4810153
      枳实651498
      甘草28076198
      浙贝母17429
      桔梗102763
      百部1101882
      前胡10111157
      紫菀9113175
      款冬花14819155
      槟榔52310
      草果595136
      藿香949147
      佩兰60673
      虎杖628159
    • 通过检索GeneCards和OMIM数据库,共收集得到COVID-19相关基因352个。将柴胡达胸合剂的作用靶点和COVID-19相关基因取交集制作Venn图,共得到51个交集基因,即柴胡达胸合剂治疗COVID-19的作用靶点(图2)。

      图  2  柴胡达胸合剂作用靶点和COVID-19相关基因Venn图

    • 柴胡达胸合剂治疗COVID-19“药材-活性成分-靶点”网络共包含234个节点(药材节点18个,有效成分节点165个,靶点节点51个)。网络中棱形代表药材节点,倒三角代表有效成分节点,圆形代表靶点节点(图3)。节点的颜色越深或节点图形越大,表明该节点在网络中的度值越高。每一圈的最低点为该圈度值最大的节点,度值沿逆时针方向逐渐减小,且3个有效成分节点圈由外往里节点度值依次减小。网络中化合物节点中位数为4,高于2倍中位数的化合物节点共有12个(表2),这些化合物可能是柴胡达胸合剂治疗COVID-19的主要活性成分。

      表 2  高于2倍度值中位数的化合物信息表

      TCMSP数据库ID化合物名称对应药材度值
      MOL000098槲皮素(quercetin)草果、柴胡、甘草、虎杖、黄连、藿香、款冬花、前胡、紫菀46
      MOL000006木犀草素(luteolin)虎杖、桔梗、佩兰、枳实、紫菀23
      MOL000422山奈酚(kaempferol)柴胡、甘草、款冬花、紫菀16
      MOL000358β-谷甾醇(beta-sitosterol)百部、半夏、浙贝母、虎杖、黄芩、款冬花、前胡、紫菀15
      MOL000173汉黄芩素(wogonin)黄芩14
      MOL004328柚皮素(naringenin)甘草、枳实13
      MOL002714黄芩素(baicalein)半夏、黄芩10
      MOL000497甘草查尔酮A(licochalcone A)甘草10
      MOL001689刺槐素(acacetin)黄芩、桔梗10
      MOL005828川陈皮素(nobiletin)枳实10
      MOL000354异鼠李素(isorhamnetin)柴胡、甘草、紫菀9
      MOL005916葛花苷元(irisolidone)藿香9

      图  3  柴胡达胸合剂治疗COVID-19“药材-活性成分-靶点”网络图

    • 在“药材-活性成分-靶点”网络中的靶点节点不包含PPI信息,因此,对柴胡达胸合剂治疗COVID-19靶点进行PPI分析,结果如图4A所示。图中节点形状越大,表明其度值越高,越可能为柴胡达胸合剂治疗COVID-19的核心作用靶点。对网络中度值前30的节点作条形图,节点度值排名前10的蛋白为CASP3、MAPK3、IL-6、MAPK8、IL-10、CXCL8、MAPK1、IL-1B、PTGS2 和CCL2(图4B)。

      图  4  靶点PPI网络图(A)和度值条形图(B)

    • GO基因注释将基因的功能分为3个部分:参与的生物学过程(biological process,BP),所处的细胞组分(cellular component,CC),执行的分子功能(molecular function,MF)。对柴胡达胸合剂治疗COVID-19的51个作用靶点进行GO基因注释,以P<0.05为条件进行筛选,结果得到GO条目共1 722个,其中BP条目1 612个,CC条目30个,MF条目80个。选取每个部分的前5个条目作气泡图,富集最多基因且P值最小的BP、CC和MF条目分别为脂多糖反应、膜筏和细胞因子受体结合(图5)。

      图  5  GO基因注释气泡图

    • 对柴胡达胸合剂治疗COVID-19的51个作用靶点进行 KEGG 信号通路富集分析,筛选出P<0.05的信号通路156条,选取富集基因最多的10条通路作气泡图。排名前5的信号通路为糖尿病并发症中的AGE-RAGE信号通路、甲型流感、IL-17信号通路、TNF信号通路和乙型肝炎(图6)。

      图  6  KEGG信号通路分析气泡图

    • COVID-19隶属于中医的“温疫”、“疫病”范畴,病因为感受“异气”、“疠气”,疠气夹湿,病位在肺、脾[3]。柴胡达胸合剂用于痰热壅肺证患者,其病因为痰热互结,壅闭于肺,致使肺失宣降而表现的肺经实热证候。本研究首先对柴胡达胸合剂的十八味中药材进行性味归经的网络分析,结果发现柴胡达胸合剂组方的性味以“苦寒”最多且主归肺经。“苦寒”药能清热泻火,消除热症,多用于具有实热特征病证[8]。因此,柴胡达胸合剂与COVID-19的病机、病位相符。

      利用中药、疾病相关数据库,本研究筛选出柴胡达胸合剂治疗COVID-19的165个活性成分和51个作用靶点,表明柴胡达胸合剂治疗COVID-19具有多成分、多靶点的特点。通过构建和分析“药材-活性成分-靶点”网络图,发现网络中度值较高的12个活性成分:槲皮素、木犀草素、山奈酚、β-谷甾醇、汉黄芩素、柚皮素、黄芩素、甘草查尔酮A、刺槐素、川陈皮素、异鼠李素和葛花苷元。除β-谷甾醇外,其余11个成分均属于黄酮类化合物。

      黄酮类化合物广泛存在于自然界的多种植物中,具有包括抗炎、抗动脉粥样硬化和抗肿瘤等多种药理作用[9]。黄酮类化合物还有良好的抗病毒作用,对流感病毒、呼吸道合胞病毒、单纯疱疹病毒和柯萨奇病毒等都有抑制作用[10]。体外实验研究表明,槲皮素和木犀草素均能够抑制SARS-CoV 3CL蛋白酶活性,对SARS-CoV 产生抑制作用[11-12]。β-谷甾醇也能抑制SARS-CoV 3CL蛋白酶活性[13]。SARS-CoV-2与SARS-CoV基因序列有约80%同源性,两者3CL蛋白酶结构有相似性[14]。此外,分子对接结果发现山柰酚、槲皮素、 黄芩素、 木犀草素、 汉黄芩素、β-谷甾醇与SARS-CoV-2 3CL蛋白酶均有较高的结合活性[15]。因此,通过直接抑制SARS-CoV-2 3CL蛋白酶活性,可能是柴胡达胸合剂治疗COVID-19的药理作用机制之一。

      通过分析PPI网络发现,CASP3、MAPK3、IL-6、MAPK8、IL-10、CXCL8、MAPK1、IL-1B、PTGS2、CCL2等在网络中有较高的度值。CASP3基因编码的caspase-3蛋白是细胞凋亡过程中重要的终末剪切酶,研究发现caspase-3蛋白在SARS-CoV病毒引起的组织细胞凋亡过程中发挥了重要作用[16]。MAPK基因编码的丝裂原活化蛋白激酶参与细胞的增殖、分化、迁移及凋亡等多个生物学过程。此外,炎症因子风暴被认为是重症COVID-19患者组织损伤的病理机制之一。SARS-CoV-2病毒在体内激活T细胞,产生粒细胞-巨噬细胞集落刺激因子(GM-CSF)和IL-6等细胞因子,随后GM-CSF会进一步激活CD14+CD16+炎性单核细胞,进一步升高IL-6等炎性因子,形成炎症因子风暴,导致严重的肺部和其他器官免疫损伤[17]。细胞因子根据其在炎性反应中的作用不同可分为促炎性细胞因子(如 IL-1、IL-6、IFN-α、IFN-γ、TNF-α 等)和抑炎性细胞因子(如 IL-4、IL-10 等)两类。SARS患者的肺部炎症和肺损伤与患者血浆中的IL-1B、IL-6、IL-12等促炎性细胞因子水平升高引起炎症因子风暴有关[18]。同样地,COVID-19患者血浆IL-1B、IFN-γ、CXC趋化因子-10(CXCL-10)等促炎性细胞因子水平也升高[19]。本研究中,IL-6、IL-10和IL-1B均为核心作用靶点。因此,减少组织细胞凋亡、降低促炎性细胞因子和升高抑炎性细胞因子水平可能也是柴胡达胸合剂治疗COVID-19的药理作用机制。

      对靶点进行GO基因注释的结果表明,柴胡达胸合剂治疗COVID-19的主要生物学过程为脂多糖反应、对源于细菌的分子的反应和氧化应激反应等。脂多糖是G-菌细胞壁的组成成分,可诱导细胞产生炎性反应,其中促炎性细胞因子发挥了重要介导作用。IL-6可促进B细胞分化,并活化MAPK,激活STAT转录因子,从而加重炎性反应[20]。KEGG分析富集到156条信号通路,主要涉及糖尿病并发症中的AGE-RAGE信号通路、甲型流感、IL-17信号通路、TNF信号通路和乙型肝炎,表明柴胡达胸合剂治疗COVID-19多通路的特点。

      综上所述,本研究采用网络药理学方法,初步揭示了柴胡达胸合剂可能一方面通过多种黄酮类化合物和β-谷甾醇直接抑制SARS-CoV-2 3CL蛋白酶活性,另一方面通过多成分、多靶点、多通路减少组织细胞凋亡、降低促炎性细胞因子和升高抑炎性细胞因子水平,从而对COVID-19产生治疗作用。

参考文献 (20)

目录

    /

    返回文章
    返回