-
新型冠状病毒肺炎（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的药理作用机制提供理论依据。
HTML
-
通过《中华人民共和国药典》（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.1. 柴胡达胸合剂药材的性味归经、活性成分和靶点的筛选
1.2. “药材-活性成分-靶点”网络图的构建
1.3. PPI网络构建
1.4. GO基因注释和KEGG通路分析
-
柴胡达胸合剂“药材-性味归经”网络如图1所示，图中节点的大小代表该节点在网络中的度值。由图可知，度值最大的性味归经分别为寒（度值=8）、苦（度值=13）和肺经（度值=14）。
-
通过TCMSP平台检索到柴胡达胸合剂中的十八味药材共包含1 977个化合物。以OB≥30%，DL≥0.18及HL≥4 h为条件筛选并去重后，共得到221个活性成分。在TCMSP检索活性成分的作用靶点，并到Uniprot数据库中查找其标准名称，最后共得到259个作用靶点。柴胡达胸合剂的“中药-化合物-活性成分-靶点”信息，结果见表1。
中药名称 化合物（个） 活性成分（个） 靶点（个） 柴胡 349 14 165 黄芩 143 32 95 法半夏 116 11 70 全瓜蒌 80 7 8 黄连 48 10 153 枳实 65 14 98 甘草 280 76 198 浙贝母 17 4 29 桔梗 102 7 63 百部 110 18 82 前胡 101 11 157 紫菀 91 13 175 款冬花 148 19 155 槟榔 52 3 10 草果 59 5 136 藿香 94 9 147 佩兰 60 6 73 虎杖 62 8 159 -
通过检索GeneCards和OMIM数据库，共收集得到COVID-19相关基因352个。将柴胡达胸合剂的作用靶点和COVID-19相关基因取交集制作Venn图，共得到51个交集基因，即柴胡达胸合剂治疗COVID-19的作用靶点（图2）。
-
柴胡达胸合剂治疗COVID-19“药材-活性成分-靶点”网络共包含234个节点（药材节点18个，有效成分节点165个，靶点节点51个）。网络中棱形代表药材节点，倒三角代表有效成分节点，圆形代表靶点节点（图3）。节点的颜色越深或节点图形越大，表明该节点在网络中的度值越高。每一圈的最低点为该圈度值最大的节点，度值沿逆时针方向逐渐减小，且3个有效成分节点圈由外往里节点度值依次减小。网络中化合物节点中位数为4，高于2倍中位数的化合物节点共有12个（表2），这些化合物可能是柴胡达胸合剂治疗COVID-19的主要活性成分。
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 
-
在“药材-活性成分-靶点”网络中的靶点节点不包含PPI信息，因此，对柴胡达胸合剂治疗COVID-19靶点进行PPI分析，结果如图4A所示。图中节点形状越大，表明其度值越高，越可能为柴胡达胸合剂治疗COVID-19的核心作用靶点。对网络中度值前30的节点作条形图，节点度值排名前10的蛋白为CASP3、MAPK3、IL-6、MAPK8、IL-10、CXCL8、MAPK1、IL-1B、PTGS2 和CCL2（图4B）。
-
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）。
-
对柴胡达胸合剂治疗COVID-19的51个作用靶点进行 KEGG 信号通路富集分析，筛选出P＜0.05的信号通路156条，选取富集基因最多的10条通路作气泡图。排名前5的信号通路为糖尿病并发症中的AGE-RAGE信号通路、甲型流感、IL-17信号通路、TNF信号通路和乙型肝炎（图6）。
2.1. 柴胡达胸合剂的药材性味归经
2.2. 柴胡达胸合剂的活性成分及其作用靶点
2.3. 柴胡达胸合剂治疗 COVID-19的作用靶点
2.4. “药材-活性成分-靶点”网络
2.5. PPI网络
2.6. GO基因注释
2.7. 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产生治疗作用。
DownLoad:
