蔡孟成，刘益群，俞超芹，金永生 （海军军医大学：. 基础医学院, . 长海医院; . 药学院，上海 200433）

清开灵是由《温病条辨》安宫牛黄丸加减化裁而来，由珍珠母、栀子、水牛角、板蓝根、金银花五味中药以及胆酸、猪去氧胆酸、黄芩苷3 种活性成分组成。具有清热解毒、镇静安神的功效，临床上常用于病毒性感冒、上呼吸道感染、急性扁桃体炎等疾病。现代研究表明，清开灵具有良好的抗病毒作用，对流感病毒H1N1、H5N1 和H7N9 以及2 型登革热病毒等均有一定的抑制作用[1-2]。冠状病毒感染引起的肺炎是由病毒感染引起，与流感的起因及症状极其相似[3]，因此，清开灵也有一定的干预冠状病毒的作用。网络药理学是近年来整合系统生物学和药理学的一门新兴学科，通过大量的生物学信息，对已有的分子生物学数据进行系统全面的分析，从蛋白、分子、基因层面揭示机体与疾病的相互作用机制，可以为多靶点药物的开发和机制研究提供依据；分子对接是通过受体的特征以及受体和药物分子之间的相互作用方式来进行药物设计的方法，通过结合能的大小判定受体与药物分子结合的能力。本文旨在通过网络药理学与分子对接技术探讨清开灵干预病毒性感冒的物质基础与作用机制，以期为今后清开灵用于临床预防和治疗病毒性感冒提供理论依据。

1 材料与方法

1.1 药物靶点的筛选

本研究采用计算系统生物学实验室的中药系统药理学数据库和分析平台 (TCMSP，http://lsp.nwu.edu.cn/tcmsp.php) 检索中药栀子、板蓝根、金银花的活性成分。因珍珠母和水牛角作为动物类药材在TCMSP 数据库中未收录，笔者应用Traditional Chinese Medicines Integrated Database (TCMID，http://119.3.41.228:8 000/tcmid/search/) 对珍珠母与水牛角进行活性成分检索。根据化合物口服生物利用度 (OB) 与类药性 (DL) 进行筛选，OB 阈值设为 OB≥30%，DL 阈值设为 DL≥0.18[4-5]。此外，由于黄芩苷易在体内代谢为黄芩素，胆酸易代谢为去氧胆酸，猪去氧胆酸易代谢为猪胆酸[6-8]，故黄芩素、去氧胆酸、猪胆酸也被纳入活性成分范围内。采用 PubChem BioActivity Analysis Service 数据库(https://pubchem.ncbi.nlm.nih.gov/assay/assay.cgi)来获取猪胆酸和猪去氧胆酸的靶点，其余活性成分靶点则通过TCMSP 获得[9]。

1.2 疾病靶点的收集

以 “influenza”为关键词检索 GeneCard（https://www.genecards.org/）数据库，收集与病毒性感冒相关的作用靶点[3]。

1.3 构建化合物-靶点网络

为明确疾病靶点与药物靶点之间的关系，将两部分靶点进行交集，输入 Cytoscape 软件构建化合物-靶点网络。Cytoscape 软件的核心架构是网络，每个节点是基因或活性成分，节点与节点之间的连接代表这些生物分子之间的相互作用，节点的度值代表网络中节点与节点相连的数目，度值越大说明该活性成分或靶点处于核心地位。

1.4 构建蛋白相互作用网络

为更深层次探究靶点之间的相互作用关系，将交集所得的共同靶点上传至线上软件STRING10.5(http://string db.org)，构建蛋白相互作用网络。把物种选择为 Homosapiens，置信度设置为＞0.9，其余参数保持默认设置，获取 PPI 网络。

1.5 生物过程与通路分析

利用 R Project 中的 org.Hs.eg.db 程序包（http://www.bioconductor.org/packages/release/data/annotati on/html/org.Hs.eg.db.html）将基因 symbols 转变为Entrez IDs，进而利用 clusterProfiler 程序包（3.6 版）（http://bioconductor.org/packages http://bioconductor.org/packages/release/bioc/html/clusterProfiler.html）对网络进行 GO 分析和 KEGG 通路分析。设定P＜0.05，并按照涉及的靶点数目多少进行排序，获取 GO 富集分析中分子功能、细胞组成和生物过程的前10 个条目以及 KEGG 富集的前 20 条信号通路。

1.6 分子对接

为进一步筛选清开灵中可能干预病毒性感冒的潜在活性成分，将清开灵中的关键化合物胆酸、猪去氧胆酸、黄芩苷与“1.4”项下筛选出的度值较高的蛋白进行分子对接。从PDB 数据库（https://www.rcsb.org/）下载靶蛋白RELA 与JUN 的三维晶体结构（PDB ID 为1VJ7 和2H7H，PDB 格式），从Pubchem 数 据 库（https://pubchem.ncbi.nlm.nih.gov/）下载胆酸、猪去氧胆酸与黄芩苷的三维结构（CAS 号为81-25-4、83-49-8、21 967-41-9，PDB 格式），采用AutoDock Tools 1.5.6 软件删除靶蛋白的水分子，分离配体和受体，添加非极性氢，计算Gasteiger 电荷，分别保存为pdbqt 格式的文件。将已经保存为PDB 格式的化合物导入AutoDock Tools 1.5.6，添加原子电荷，分配原子类型，所有柔性键均默认可旋转，保存为pdbqt 格式，作为对接配体。运行Auto Dock 软件进行对接，采用Pymol软件可视化对接结果[10]。

2 结果

2.1 潜在活性成分的筛选与对应靶点的收集

通过TCMSP 数据库检索到板蓝根活性成分169 个，金银花活性成分236 个，栀子活性成分98 个，通过TCMID 数据库检索到珍珠母活性成分9 个，水牛角活性成分6 个。根据OB≥30%与DL≥0.18筛选出板蓝根活性成分39 个，金银花活性成分23 个，栀子活性成分15 个，其中MOL000098（槲皮素）、MOL000422（山奈酚）、MOL001194（亚油酸乙酯）和MOL003095（甘露醇）为栀子与金银花的共有成分，MOL000358（β-谷甾醇）和MOL000449（豆甾醇）为栀子、金银花和板蓝根的共有成分（图1 中绿色节点表示活性成分，红色节点表示药物，橙色圈内节点为药物共有成分）。通过TCMSP 数据库获取上述中药活性成分以及黄芩苷、黄芩素、胆酸、去氧胆酸的靶点，通过PubChem BioActivity Analysis Service 获取猪去氧胆酸以及猪胆酸的靶点，除去重复靶点，共获得靶点225 个。

2.2 疾病相关靶点的筛选

从GeneCard 数据库中获得病毒性感冒靶点2 506 个。

2.3 化合物-靶点网络构建

将“2.1”中的225 个清开灵活性成分靶点与“2.2”中的2 506 个病毒性感冒靶点输入Venny2.1软件，绘制韦恩图，两者进行交集后获得疾病-药物成分共同靶点100 个，见图2。采用Cytoscape 软件构建的活性成分-靶点网络，该网络共包含147 个节点(47 个化合物节点，100 个靶点节点)和315 条边，其中，红色节点表示化合物，绿色节点表示药物靶点，连接化合物与靶点的边表示它们之间有作用关系，见图3。表1 列出了化合物-靶点网络中的关键节点、类别及其度值。活性成分度值排名前5 位的是槲皮素、木犀草素、山奈酚、β-胡萝卜素和黄芩素，这5 种化合物的节点数分别是71、31、23、18 和15。靶点度值排名前5 位的是前列腺素过氧化物合酶2（PTGS2）、前列腺素过氧化物合酶1（PTGS1）、丝氨酸蛋白酶1（PRSS1）、一氧化氮合酶（NOS2）和雌激素受体1（ESR1），分别能与44、31、15、14 和13 个化合物发生相互作用，体现了中药多成分、多靶点的特性。

2.4 PPI 网络的构建与分析

图2 清开灵化学成分作用靶点与病毒性感冒相关靶点交集分析

应用 STRING 软件构建PPI 网络(图4)，在此网络图中包含88 个节点，718 条边，同时得到网络中关键靶点的频次 (图5)。根据“度值＞均值”筛选出关键节点34 个，包括：JUN、RELA、MAPK1、IL6、AKT1、MAPK14、CXCL8、APP、IL1B、MAPK8、VEGFA、FOS、EGFR、IL-4、ESR1、IL2、PRKCA、PRKCB、CASP8、CCL2、MYC、NFKBIA、CASP3、CTNNB1、 IFNG、 IL-10、 CXCL2、 RB1、 ARG1、CHUK、CXCL10、MMP9、NOS2、STAT1。度值最高的是JUN，度值为29；其次为RELA，度值为27；MAPK14 度值为25，IL-6 度值为23。度值大的靶点提示在网络调控中起着关键作用，度值大的靶点很可能是清开灵治疗 病毒性感冒的关键靶点。

图3 清开灵活性成分-靶点网络

表1 化合物-靶点网络的关键节点及其度值

2.5 GO 与KEGG 富集分析

GO 富集条目（图6）生物过程主要包括：脂多糖反应、 氧化应激反应、金属离子反应、细胞脂多糖反应，分子功能包括：细胞因子受体结合、细胞因子活性、血红素结合、 磷酸酶结合，其在膜筏、膜微区、膜区、囊腔等处富集。KEGG 富集信号通路（图7）主要包括：卡波济肉瘤相关疱疹病毒感染、乙型肝炎、人巨细胞病毒感染、流体剪切应力与动脉粥样硬化、糖尿病并发症的AGE-RAGE 信号通路等。其中，卡波济肉瘤相关疱疹病毒感染的靶点主要包括IL-6、STAT1、p38、CASP8、Bax、CASP3、CASP9、 CyrC、 IKK、 RB、 E2F、 ICAM1、 Raf、ERK、 AKT、 IKBa、 NFKB、 JNK、 AP-1、 IL-8、COX2、VEGF 和Myc；乙型肝炎的靶点主要包括：STAT1、JNK、p38、AP-1、ERK、IKK、IKB、NFKB、 PKC、 Raf、 STAT、 IL-8、 PCNA、 MMP-9、Bax、Bcl2、CASP9、CASP3、Akt、SAPK、pRb 和E2F（图8、图9）。

图4 清开灵蛋白质-蛋白质相互作用网络

2.6 分子对接

一般认为配体与受体结合的稳态构象能量越低，产生的效应越强。笔者对清开灵中的胆酸、猪去氧胆酸合黄芩苷进行分子对接，对接模式见图10，其中，蓝紫色部分代表蛋白质，红黄色部分代表活性成分。胆酸、猪去氧胆酸合黄芩苷与RELA 和JUN 对接的结合能见表2。其中，HQG表示黄芩苷，ZQYDS 表示猪去氧胆酸，DS 表示胆酸。

3 讨论

清开灵由珍珠母、栀子、水牛角、板蓝根、金银花五味中药以及胆酸、猪去氧胆酸、黄芩苷3 种活性成分组成。金银花中槲皮素、木犀草素、绿原酸等是抗流感病毒的主要活性成分，对甲型流感病毒、疱疹病毒及腺病毒均有一定的抑制作用[11]。栀子提取物对流感病毒感染所致小鼠的肺部炎症有明显抑制作用，在体外对流感病毒等6 种病毒引起的细胞病变有明显抑制作用[12]。板蓝根乙酸乙酯部位的生物碱、有机酸和氨基酸衍生物与抗病毒活性密切相关[13]。现代网络药理学研究表明水牛角-珍珠母药对的解热作用可能与水牛角内的氨基酸代谢与珍珠母内金属离子代谢有关[14]；胆酸-猪去氧胆酸药对可能通过作用于MR 影响水、电解质代谢，从而发挥抗炎作用[9]。黄芩苷可以通过抑制H1N1 流感病毒的中后期mRNA 的转录和复制，抑制神经氨酸苷酶活性，以及抑制病毒诱导的巨噬细胞和肺上皮细胞中促炎基因的表达，从而影响 IL-6和 IL-8 等的合成[15-16]。综上所述，现代研究表明清开灵干预病毒性感冒的机制可能是通过抗病毒、抗炎实现的。

图5 清开灵作用的关键靶点

图6 清开灵GO 富集分析

本研究运用网络药理学研究方法共筛选出清开灵中活性成分90 个，作用靶点225 个，体现了中药多成分、多靶点联合作用的整体性和关联性的特点。根据“活性成分-靶点”网络结果，进一步得到排名前5 的关键化合物为槲皮素、木犀草素、山奈酚、β-胡萝卜素和黄芩素。其中，槲皮素、山柰酚、木犀草素、黄芩素经研究发现也是目前抗击新冠疫情“三药三方”中连花清瘟胶囊、金花清感颗粒、清肺排毒汤等的核心化合物[17-19]。现代研究表明，槲皮素可以与F 蛋白相互作用阻断病毒的黏附，发挥抗人呼吸道合胞病毒的作用[20]；木犀草素可通过干扰外壳蛋白I 复合物的表达，降低甲型流感病毒的含量[21]；β-胡萝卜素对腺病毒转化的293 细胞和人类疱疹病毒转化的 Raji 细胞的生长增殖具有明显的抑制作用，对293 细胞内整合的腺病毒早期基因的表达具有下调作用[22]。可见，实验研究进一步证实了清开灵的活性成分具有一定的抗病毒作用。此外，山奈酚可以减轻脂多糖导致的小鼠急性肺损伤，其机制可能是通过抑制肺组织促炎因子NFκB、IL-1β、HO-1、PTGE、iNOS、TNF-α 和COX-2mRNA 的表达和促进抗炎因子HO-1mRNA 的表达而发挥保护作用[23]；槲皮素、木犀草素、β-胡萝卜素也具有抗炎活性[24-25]。综上可知，抗炎作用也是清开灵干预病毒性感冒的重要机制。

图7 清开灵KEGG 富集分析（气泡图）

图8 卡波济肉瘤相关疱疹病毒感染通路图

图9 乙型肝炎通路图

根据PPI 网络分析结果和“活性成分-靶点”网络结果，清开灵干预病毒性感冒的主要作用靶点有前列腺素过氧化物合酶2（PTGS2）、前列腺素过氧化物合酶1（PTGS1）、丝氨酸蛋白酶1（PRSS1）、一氧化氮合酶（NOS2）、雌激素受体1（ESR1）、细胞因子（IL-6）和丝裂原活化蛋白激酶家族成员（MAPK8）等。GO 功能富集分析结果发现，清开灵活性成分主要涉及氧化应激、细胞应答、炎性反应、免疫应答等生物学过程；KEGG 通路富集主要涉及卡波济肉瘤相关疱疹病毒感染，乙型肝炎，人巨细胞病毒感染，流体剪切应力与动脉粥样硬化，糖尿病并发症的AGE-RAGE 信号通路等通路，靶点富集数目最多的通路为病毒感染相关通路，如卡波济肉瘤相关疱疹病毒感染、乙型肝炎、人巨细胞病毒感染等，推测清开灵对病毒引起的感染具有一定的干预作用。现代研究表明，金银花中的绿原酸可以抑制流感病毒结构蛋白神经氨酸酶，从而阻断子代病毒的释放与传播；板蓝根中的橙黄胡椒酰胺酯可降低RNA 聚合酶活性，同时可抑制NF-κB 因子的转录水平；黄芩中的黄芩苷不仅可以抑制神经氨酸酶还可抑制流感病毒RNA 聚合酶[26]。此外，以上通路均涉及炎症基因，如COX-2、IL-1β、IL-6、IFN-γ 等，同时，TNF 信号通路与IL-17 信号通路在KEGG 富集分析中亦有较高的相关度：TNF 信号通路在抗病毒、介导能量代谢、炎症反应、血管生成以及维持干细胞特性等方面发挥重要作用，IL-17 信号通路与免疫反应及炎症反应的发生与发展密切相关[10]。综上所述，清开灵可能通过抗病毒、抗炎途径起到干预病毒性感冒的作用。

图10 清开灵主要成分与RELA 和JUN 分子对接模式

表2 清开灵核心成分与RELA 和JUN 的结合能

胆酸、猪去氧胆酸以及黄芩苷作为清开灵中的主要成分，JUN 与RELA 为PPI 网络中的关键蛋白。分子对接结果显示，这3 种活性成分与JUN 和RELA 均有一定的结合能力。受体内代谢因素影响，如黄芩苷会代谢为黄芩素，胆酸、猪去氧胆酸亦会在体内代谢为去氧胆酸和猪胆酸。以黄芩苷为例，从数据库的靶点数目来看，黄芩苷的靶点数目不如黄芩素，笔者推测黄芩素多靶点的原因可能是其结构上更易与靶点蛋白结合。而清开灵中应用黄芩苷而不用黄芩素可能是因为黄芩苷不易产生消化道副作用、成本低或者与药动学有关，但实际情况仍需后续探究。

急性病毒性感冒属于中医瘟疫的范畴，具有强烈的传染性与流行性，如SARS 与COVID-19 均属于瘟疫，治疗上大多应用清热解毒药。清开灵中金银花、栀子、板蓝根、水牛角、珍珠母均为寒凉之品，黄芩与猪胆汁在中医学中亦属寒凉，故笔者推测黄芩中提取出的黄芩苷与猪胆汁中提取出的胆酸和猪去氧胆酸亦属寒凉，根据中医理论推测其干预病毒性感冒具有一定作用。临床研究数据显示，清开灵联合奥司他韦可以提高小儿咽部甲型流感病毒的清除率[27]，与儿童清肺口服液联用可治疗小儿呼吸道合胞病毒性肺炎[28]，联合利巴韦林可治疗流行性感冒[29]。综上分析，清开灵用于病毒性感冒的治疗可能与其抗病毒、抑制炎症反应相关，进一步证实了网络药理学分析结果的合理性。但近些年对清开灵注射液的安全性讨论层出不穷，使得人们对该药物望而却步[30]；而且中医学认为该药苦寒败胃，不能够长期应用，故此，在临床上我们仍需辨证论治。