连花清瘟胶囊提高COVID-19临床治愈率的机制：基于网络药理学和分子对接技术

鄢海燕，邹纯才

皖南医学院药学院，安徽 芜湖 241002

中华人民共和国国家卫生健康委办公厅和中华人民共和国国家中医药管理局办公室发布的《新型冠状病毒肺炎诊疗方案》（试行第七版）（国卫办医函〔2020〕184号）［1］明确新型冠状肺炎（COVID-19）临床主要表现为发热、干咳和乏力。在中医治疗中，明确COVID-19属于中医“疫”病范畴，病因为感受“疫戾”之气，各地可根据病情、当地气候特点以及不同体质情况进行辩证论治，并推荐连花清瘟胶囊用于医学观察期的治疗。

Hu等［2］的最新研究论文阐述了连花清瘟胶囊用于临床治疗COVID-19的疗效，即经过连花清瘟治疗组治疗14 d（4粒/次，3次/d）后，连花清瘟胶囊治疗组的恢复率达91.5%，明显高于对照组（82.4%），且连花清瘟胶囊治疗组对于发热、乏力、咳嗽等症状的治愈时间中位数明显缩短。该文认为，从安全性和有效性角度来看，可以考虑使用连花清瘟胶囊改善COVID-19的临床症状（发热、干咳和乏力）。

以系统生物学理论为基础的网络药理学，通过对生物系统的网络分析来选取特定信号节点，构建活性成分-蛋白靶点-信号通路之间的复杂网络来探讨药物的作用机制。现有网络药理学研究多关注药物作用于靶点蛋白用于疾病的治疗［3-4］，对于疾病与其临床症状共有靶点蛋白间的网络药理学研究少有报道。COVID-19的临床主要表现为发热、干咳和乏力。涉及发热、干咳和乏力等临床症状的靶点蛋白较多，就COVID-19而言，这些临床症状似应与COVID-19间存在共有的靶点蛋白。连花清瘟胶囊改善COVID-19的临床症状提高临床治愈率似应是作用于COVID-19及其临床症状（发热、干咳和乏力）共有靶点蛋白的结果。为此，为更全面地研究连花清瘟胶囊提高COVID-19临床治愈率的机制和物质基础，本文以“fever”、“coughing ”、“fatigue ”等临床表现及“2019-nCoV”、“COVID-19”为疾病关键词，筛选出“fever”、“coughing ”、“fatigue ”、“2019-nCoV”、“COVID-19”中的共同靶点蛋白。同时选取连花清瘟胶囊中连翘、金银花等13味中药，对涉及到的活性成分进行筛选，通过网络药理学方法研究连花清瘟胶囊提高COVID-19临床治愈率的潜在作用靶点和信号通路及生物学过程，并采用分子对接技术研究核心活性成分与主要通路中的靶点蛋白的对接，探讨连花清瘟胶囊提高COVID-19临床治愈率的机制和物质基础。

1 资料和方法

1.1 数据资料来源

1.1.1 药材来源 连花清瘟胶囊（国药准字Z20040063）由石家庄以岭药业股份有限公司生产。成分为：连翘、金银花、炙麻黄、炒苦杏仁、石膏、板蓝根、绵马贯众、鱼腥草、广藿香、大黄、红景天、薄荷脑、甘草。功能主治为清瘟解毒、宣肺泄热。用于治疗流行性感冒属热毒袭肺证，症见：发热或高热，恶寒，肌肉酸痛，鼻塞流涕，咳嗽，头痛，咽干咽痛，舌偏红，苔黄或黄腻等。

1.1.2 活性成分来源 通过中药系统药理学数据库与分析平台(TCMSP，Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform，https://tcmspw.com/tcmsp.php)及有关文献［5］，分别以连花清瘟胶囊中的13味药来检索活性成分，其中炙麻黄、炒苦杏仁、绵马贯众分别按麻黄、苦杏仁、贯众进行检索。选择符合口服生物利用度（oral bioavailability,OB）≥30%、类药性（drug like,DL）≥0.18并合并各药味中相同的活性成分，绘制结构式并以MDL SDfile（*.sdf）格式保存（检索日期：2020年10月28日）。将这些活性成分的MDL SDfile（*.sdf）格式结构式导入SwissTargetPrediction数据库（http://www.swisstargetprediction.ch/），检索符合Probability≥0.1的活性成分靶点蛋白（检索日期：2020年11月2日）。

1.1.3 候选靶标来源 利用CooLGeN(http://ci.smu.edu.cn/CooLGeN/Home.php)、GeneCards（https://www.genecards.org/） 、TTD（http://bidd.nus.edu.sg/group/cjttd/）数据库，以“fever”、“coughing ”、“fatigue ”、“2019-nCoV”、“COVID-19”为关键词，收集临床表现及疾病靶点蛋白（选择标示符为All human genes），筛选临床表现与疾病共有靶点蛋白，再进一步选择与药材中活性成分相匹配的靶点蛋白。

1.2 靶点通路及生物学功能注释分析

将靶点蛋白复制至DAVID数据库（https://david.ncifcrf.gov/tools.jsp）的列表中，选择标示符为official_gene_symbol，物种注释为Homo sapiens，进行靶点蛋白的GOTERM_BP_DIRECT、GOTERM_CC_DIRECT、GOTERM_MF_DIRECT富集分析和Kyoto Encyclopedia of Genes and Genomes（KEGG）_PATHWAY通路注释分析。

1.3 活性成分-蛋白靶点-生物学功能网络分析和构建

根据预测结果，使用Gephi0.9.2软件将筛选出来的活性成分、蛋白靶点和生物学功能构建活性成分-蛋白靶点-生物学功能网络。

1.4 分子对接分析

本文以SYBYL-X 2.1.1对接软件对靶点蛋白和化学成分进行分子对接分析，获得Total Score值。

2 结果

2.1 连花清瘟胶囊各药味中活性成分

通过TCMSP及有关文献［5］共检索到连花清瘟胶囊各药味中符合OB≥30%、DL≥0.18的活性成分（合并各药味中相同的活性成分）210个。

2.2 靶点的预测

将“2.1”中符合OB≥30%、DL≥0.18的活性成分以MDL SDfile（*.sdf）格式导入SwissTargetPrediction数据库，获得蛋白靶点信息24215条，涉及1320个化学靶点蛋白。筛选符合Probability≥0.1的活性成分靶点蛋白，获得蛋白靶点信息12274条，涉及1159个化学靶点蛋白。

利用CooLGeN、GeneCards、TTD数据库，以“2019-nCoV”、“COVID-19”为关键词，共收集疾病靶点蛋白（All human genes）261个（检索日期：2020年10月28日）。以“fever”、“coughing ”、“fatigue”为关键词，分别收集疾病临床表现靶点蛋白7237、4157、7127个。“fever”、“coughing”、“fatigue”与COVID-19的共有靶点蛋白分别为154、176、106个。筛选各药材中与“fever”、“coughing”、“fatigue”与COVID-19的共有靶点蛋白相同的靶点蛋白及对应的化学成分。并进一步确定药材、临床表现及COVID-19中共同的靶点蛋白共计57个（图1）。

活性成分与COVID-19及其临床表现（发热、咳嗽、乏力）的共有靶点蛋白共计57个：ABL1、ACE、ACE2、AGTR1、AKT1、ALK、AR、BTK、CALCA、CASP1、CASP3、CCL2、CCR5、CTSL、CXCL8、DHODH、DPP4、EGFR、FLT3、G6PD、HPSE、HSP90AA1、HSPA5、IDH1、IL2、IL6、JAK1、JAK2、LNPEP、MAPK1、MCL1、MGMT、MPO、MTOR、NFKB1、NLRP3、NR3C1、P2RX7、PARP1、PIK3CG、PLAT、PLAUR、PLG、PPARA、REN、ROS1、TK1、TLR4、TLR7、TLR9、TNF、TRPV1、TRPV4、TTR、VCP、VEGFA、XPO1，涉及11味药、160个活性成分（注：本文活性成分均以Mol ID表示，具体成分名称可检索TCMSP数据库）。药材、化学成分及靶点蛋白信息（以 MAPK1、NLRP3、HSP90AA1、TLR9、AKT为代表）：

AKT1靶点蛋白：MOL000006（金银花、连翘、炙麻黄），MOL000098（甘草、广藿香、金银花、连翘、鱼腥草、炙麻黄），MOL000239、MOL000354、MOL004828、MOL004863、 MOL004866、 MOL004891、MOL004904、MOL004949、MOL004961（甘 草），MOL000422（甘草、金银花、连翘、绵马贯众、鱼腥草、炙麻黄），MOL001495、MOL003044、MOL003117（金银花），MOL001735（板蓝根），MOL002235（大黄），MOL002823、MOL002881、MOL005842（炙麻黄）；

HSP90AA1靶点蛋白：MOL000500、MOL004806、MOL004838、 MOL004849、 MOL004855、MOL004856、 MOL004857、 MOL004863、MOL004864、 MOL004883、 MOL004884、MOL004904、 MOL004905、 MOL004910、MOL004911、 MOL004935、 MOL004945、MOL004949、 MOL004959、 MOL004980、MOL004988、 MOL004989、 MOL005001、MOL005008、（甘草），MOL001040（绵马贯众），MOL002311、MOL004908（甘 草 、苦 杏 仁），MOL003283、MOL003370（连翘），MOL005016（甘草、广藿香）；

MAPK1靶点蛋白：MOL000522、MOL003283、MOL003330、MOL003370（连 翘），MOL001728、MOL001774、MOL001803（板蓝根），MOL004806、MOL004808、 MOL004838、 MOL004891、MOL004905、MOL004959、MOL004988、MOL004989（甘草），MOL004908（甘草、苦杏仁），MOL005911（广藿香）；

NLRP3靶点蛋白：MOL001495（金银花），MOL001494（金银花、炙麻黄）；

TLR9 靶 点 蛋 白 ：MOL000392、MOL000417、MOL002844、 MOL004815、 MOL004833、MOL004879、 MOL004885、 MOL004891、MOL004907、MOL004957、MOL004990、MOL005016（甘草），MOL001728、MOL001749（板蓝根），MOL004908（甘草、苦杏仁）。

2.3 靶点信号通路与Gene Ontology（GO）分析

将连花清瘟胶囊中各药味活性成分与COVID-19及其临床表现（发热、咳嗽、乏力）共有的57个靶点蛋白导入DAVID数据库，进行靶点蛋白的GOTERM_BP_DIRECT、 GOTERM_CC_DIRECT、GOTERM_MF_DIRECT富集分析和KEGG_PATHWAY通路注释分析。信号通路主要涉及NOD-like receptor signaling pathway、Toll-like receptor signaling pathway、Renin-angiotensin system等35条信号通路。其中Benjamini correction＜0.01的信号通路共有10条（图2）。

图2 连花清瘟胶囊中活性成分潜在靶点的KEGG信号通路富集分析Fig.2 Enrichment analysis of KEGG signal pathway of potential targets of the active components in Lianhua Qingwen capsule.

连花清瘟胶囊中的160个活性成分通过MAPK1，IL6，HSP90AA1，TNF，CCL2，NFKB1，CASP1，NLRP3等靶点蛋白参与COVID-19及其临床表现（发热、咳嗽、乏力）的生物过程、细胞组成和分子功能。生物过程主要与regulation of apoptosis、regulation of programmed cell death等有关，细胞组成主要与cell surface、extracellular space等有关，分子功能主要与protein tyrosine kinase activity、ATP binding等有关（图3、4）。

图3 连花清瘟胶囊中活性成分潜在靶点与COVID-19及其临床表现（发热、咳嗽、乏力）相关的GOTERM_BP_FAT富集分析Fig.3 GOTERM_BP_FAT enrichment analysis of the potential targets of active components in Lianhua Qingwen capsule and the biological functions related to COVID-19 and its clinical manifestations(fever,cough,and fatigue)(Benjamini correction≤0.0001).

图4 连花清瘟胶囊中活性成分潜在靶点与COVID-19及其临床表现（发热、咳嗽、乏力）相关的富集分析Fig.4 GOTERM_CC_FAT(A)and GOTERM_MF_FAT(B)enrichment analysis of the potential target of active components in Lianhua Qingwen capsule and the biological function related to COVID-19 and its clinical manifestations(fever,cough,fatigue)(Benjamini correction≤0.01).

2.4 活性成分-靶点蛋白-信号通路网络的构建

使用Gephi 0.9.2软件构建化学成分-靶点蛋白-信号通路网络。连花清瘟胶囊中的160个活性成分通过57个靶点蛋白、35条信号通路发挥治疗COVID-19或干预COVID-19病变过程的作用，网络关系复杂、蛋白靶点与信号通路多（图5）。

图5 连花清瘟胶囊中药味-部分活性成分-靶点蛋白-作用通路网络图Fig.5 Network diagram of medicinal ingredients()-active component()-target protein()-action pathway()in Lianhua Qingwen capsule.

2.5 部分活性成分与靶点蛋白分子对接结果分析

鉴于 NOD-like receptor signaling pathway、Tolllike receptor signaling pathway的Benjamini correction较小，为5.8E-5、1.5E-6。为此，以NOD-like receptor signaling pathway、Toll-like receptor signaling pathway涉及的靶点蛋白MAPK1（PDB ID:6SLG）、IL6（PDB ID:1ALU）、HSP90AA1（PDB ID:3O0I）、TNF（PDB ID:6M95）、CCL2（PDB ID:5J5Y）、NFKB1（PDB ID:1LE9）、CASP1（PDB ID:1RWV）、NLRP3（PDB ID:3VWE）、PIK3CG（PDB ID:2A4Z）、AKT1（PDB ID:4GAH）、TLR4（PDB ID:3FXI）、TLR7（PDB ID:5T1S）、TLR9（PDB ID:5WYX）为例，对靶点蛋白和化学成分进行分子对接结果分析，验证连花清瘟胶囊化学成分与靶点蛋白的作用。配体与受体间低能量的稳定构象预示着二者间有较大的作用可能性。一般以结合能≤-5kJ/mol或Total Score值≥5.0作为筛选标准。本文以SYBYL-X 2.1.1对接软件进行分子对接，获得Total Score值（表1）。分子对接结果表明，MOL000522、MOL001495、MOL001494等83个化学成分与MAPK1、IL6、HSP90AA1等 12个靶点蛋白的 Total Score值≥5.0。Total Score值≥9.0的代表性核心活性成分与靶点蛋白分子对接模式（图6）。

图6 代表性核心活性成分与靶点蛋白分子对接模式Fig.6 Molecular docking mode of representative core active components and target proteins(Total Score≥9.0).

表1 NOD-like and Toll-like receptor signaling pathway各化学成分与靶点蛋白作用的Total Score值Tab.1 Total score values of each chemical component interacting with target protein in NOD like and total like receptor signaling pathway

续表1-2

续表1-3

3 讨论

通过网络药理学研究连花清瘟胶囊与COVID-19的临床表现（“fever”、“coughing”、“fatigue”）关系结果可知，在涉及到的35条信号通路中，Benjamini correction最小的2个信号通路为NOD-like receptor signaling pathway、Toll-like receptor signaling pathway，Benjamini correction分别为5.8E-5、1.5E-6。NOD受体可识别进入细胞内的病原微生物及其产物，快速启动信号传递，激活天然免疫效应机制；Toll受体可识别细胞外病原体，并将信号传递至细胞内激发机体天然免疫系统。二种模式识别受体可独立进行自我辨识，同时又相互联系协调，在启动天然免疫反应中共同发挥重要作用，可激活免疫系统，破坏免疫耐受状态，纠正自身免疫性疾病［6-7］。连花清瘟胶囊提高COVID-19的临床治愈率可能主要与提高患者自身免疫能力、抑制炎症反应有关。为此，以NOD-like receptor signaling pathway、Toll-like receptor signaling pathway涉及的靶点蛋白及其化学成分进行了分子对接结果分析，验证连花清瘟胶囊化学成分与靶点蛋白的作用。MOL000522、MOL001495、MOL001494等83个化学成分与MAPK1、IL6、HSP90AA1等12个靶点蛋白的Total Score值≥5.0，Total Score值≥9.0涉及到的主要靶点蛋白有MAPK1、NLRP3、HSP90AA1、TLR9、AKT1。

MAPK可调节细胞增殖、分化、死亡、应激反应和凋亡等关键细胞活动［8］，在细胞生命中起着重要作用［9］。MAPK1是多种生化信号的整合点，参与细胞增殖分化、转录调控和细胞的有丝分裂自噬等多种生物学过程，在治疗各种炎症性疾病中发挥重要作用，也是急性肺损伤（ALI）炎症反应和LPS诱导的细胞损伤的重要调节因子［10-14］。由本文网络药理学和分子对接结果可知，来自连翘、甘草、板蓝根、苦杏仁、广藿香的14个化学成分与MAPK1靶点蛋白的分子对接结果Total Score值＞5.0，其中 MOL000522（连翘）、MOL004989（甘草）、MOL003330（连翘）与MAPK1靶点蛋白的分子对接结果Total Score值＞9.0。连花清瘟胶囊作用于MAPK1靶点蛋白，可能通过细胞凋亡调控、细胞程序性死亡的调控等生物学途径，提高患者自身免疫能力、抑制炎症反应发挥提高COVID-19临床治愈率作用。

NLRP3炎症小体作为机体固有免疫及应激系统的重要防御成分，不但可以识别细菌、病毒等病原体引发机体发生固有免疫应答，还负责炎症反应的激活［15-17］，参与了多种疾病发生和进展［18］，在先天免疫中起着重要作用［19］。近年来，愈来愈多的研究表明NLRP3炎性体在肺部感染性疾病中参与并发挥多种作用［20］。如在各种类型ALI中介导了炎症介质的生成和炎症细胞的浸润，增加肺泡上皮细胞的通透性，促进肺水肿的形成［21］。抑制或缺失NLRP3可特异性减轻辐射和脂多糖治疗引起的小鼠肺部炎症［22］。中药成分也可通过抑制NLRP3炎症小体活化，继而减轻脂多糖诱导的ALI［23-25］。由本文网络药理学和分子对接结果可知，MOL001494（金银花、炙麻黄）、MOL001495（金银花）与NLRP3靶点蛋白的分子对接结果Total Score值＞9.0。连花清瘟胶囊作用于NLRP3靶点蛋白，可能通过细胞凋亡调控、细胞程序性死亡的调控等生物学途径，提高机体固有免疫及应激功能，激活并参与炎症反应，在提高COVID-19临床治愈率中发挥重要作用。

作为鸟流感病毒结合受体的细胞表面蛋白HSP90AA1在感染早期通过HSP90AA1-AKT-MTOR途径诱导自噬。一旦病毒识别，HSP90AA1和AKTMTOR通路的直接连接会触发自噬，这是控制感染的关键步骤［26］。由本文网络药理学和分子对接结果可知，来自连翘、甘草、苦杏仁、广藿香、绵马贯众的30个化学成分与HSP90AA1靶点蛋白的分子对接结果Total Score值＞5.0，其中MOL004908（甘草、苦杏仁）、MOL004863、MOL004806（甘草）与HSP90AA1靶点蛋白的分子对接结果Total Score值＞9.0。连花清瘟胶囊作用于HSP90AA1靶点蛋白，可能通过对有机物的反应、细胞生物合成过程的正调控等生物学途径，触发自噬控制病毒感染，在提高COVID-19临床治愈率中发挥重要作用。

TLR9可激活固有免疫系统，参与传染病检测［27-28］。TLR9和磷脂酰肌醇-3-激酶γ（PI3Kγ）是免疫应答中非常重要的效应因子。免疫应答需要PI3Kγ，其中TLR9是相关的触发因子［29-30］。在稳态条件下，不适当的TLR9反应可导致严重的自体炎症性疾病，阻断TLR9炎症通路被激活，控制疾病进展和炎症并发症［31］。在肺缺血再灌注过程中释放的线粒体DNA（mtDNA）会触发TLR9依赖性网络的形成，并导致肺损伤［32］。由本文网络药理学和分子对接结果可知，来自板蓝根、甘草、苦杏仁的14个化学成分与TLR9靶点蛋白的分子对接结果Total Score值＞5.0，其中MOL001749（板蓝根）与TLR9靶点蛋白的分子对接结果Total Score值＞9.0。连花清瘟胶囊作用于TLR9靶点蛋白，可能通过对伤害的反应、细胞因子产生的调节、细胞生物合成过程的正调控等生物学途径，激活固有免疫系统，在预防疾病进展和炎症并发症的发展中起到作用［31］，提高COVID-19临床治愈率。

AKT1是调节细胞存活的信号通路中的一个中心节点。AKT1调控的多种途径在细胞内通过100多种细胞底物的磷酸化进行沟通［33］。如肺纤维化和高碳酸血症。特发性肺纤维化是一种进行性间质性肺炎，以成纤维细胞聚集、胶原沉积和细胞外基质重塑为特征。AKT1通过诱导巨噬细胞产生IL-13来调节肺纤维化，表明靶向AKT1可能同时阻断特发性肺纤维化的纤维化过程［34］。高碳酸血症是严重急性和慢性肺部疾病患者死亡的危险因素。靶向AKT1或升高CO2信号的下游通路可以增强巨噬细胞抗病毒宿主的防御能力，并改善晚期肺部疾病高碳酸血症患者的临床预后［35］。由本文网络药理学和分子对接结果可知，来自甘草、金银花、连翘、绵马贯众、鱼腥草、炙麻黄、广藿香、金银花、板蓝根、大黄的18个化学成分与AKT1靶点蛋白的分子对接结果Total Score值＞5.0，其中MOL001495（金银花）与AKT1靶点蛋白的分子对接结果Total Score值＞9.0。连花清瘟胶囊作用于AKT1靶点蛋白，可能通过细胞凋亡调控、细胞程序性死亡的调控、细胞死亡调控、多细胞生物过程的正调控等生物学途径，在提高COVID-19临床治愈率中发挥重要作用。

COVID-19属中医肺络病，其临床病理特征是免疫系统功能障碍和炎症反应引起的深气道、肺泡的损害，还可见重度肺充血。连花清温胶囊以络病理论为指导，对许多病毒如sars-CoV、mers-CoV和炎症反应有抑制作用，已纳入国家的COVID-19诊断和治疗计划，其理论组方特色和临床基础研究已得到广泛认可，是治疗呼吸系统疾病的代表性药物［36］。最新研究表明［36］，连花清瘟胶囊可抑制体外培养的COVID-19活性，显著缓解COVID-19患者的发热、咳嗽、乏力等症状，为连花清瘟胶囊治疗COVID-19提供了理论和临床依据。