-
2019年12月，湖北省武汉市爆发新型冠状病毒肺炎（简称“新冠肺炎”，coronavirus disease 2019，COVID-19），由于新型冠状病毒（SARS-CoV-2）传染力强、传播范围广，我国将新冠肺炎定义为急性呼吸道传染病并按甲类传染病进行管理。有效的管控措施使得目前国内疫情得到了一定控制，但中、重症患者的治疗情势依然严峻复杂[1-2]。
金振口服液由羚羊角、平贝母、大黄、黄芩、牛黄、青礞石、生石膏、甘草八味中药构成，具有清热解毒，祛痰止咳的功效，主要用于急性支气管炎符合痰热咳嗽者，表现为发热、咳嗽、咳吐黄痰；现代药理学研究表明，金振口服液中含有芦荟大黄素、大黄酚、黄芩素等成分[3-4]。目前，金振口服液已经成为援助一线的重要物资[5]。此外，在第四期中国新冠疫情防控经验国际分享交流会上，钟南山院士指出金振口服液具有较好的体外抗病毒与抗炎活性，其TC50值为29.44 mg/ml，EC50值为5.85 mg/ml，SI值为5.03。本文旨在通过对金振口服液进行网络药理学分析，以期明确其物质基础及可能的作用机制，为抗疫贡献一份力量。
HTML
-
本研究采用计算系统生物学实验室的中药系统药理学数据库和分析平台(TCMSP) (
http://lsp.nwu.edu.cn/tcmsp.php )检索中药羚羊角、平贝母、大黄、黄芩、牛黄、青礞石、生石膏、甘草的活性成分。因石膏、青礞石和羚羊角作为矿物类和动物类药材的特殊性，故本文方法暂时未将其纳入分析范围[6]。根据化合物口服生物利用度(OB)与类药性(DL)进行筛选，OB阈值设为OB≥30%，DL阈值设为DL≥0.18，最终筛选出154个化合物，入血化合物141个并采用TCMSP数据库获得对应的靶点。 -
以novel coronavirus pneumonia为关键词检索GeneCard（https://www.genecards.org/）数据库，收集与新冠肺炎相关的作用靶点。
-
为明确疾病靶点与药物靶点之间的关系，将两部分靶点进行交集，输入Cytoscape3.6.1软件构建出化合物-靶点网络。Cytoscape软件的核心架构是网络，每个节点是基因或分子，节点与节点之间的连接代表这些生物分子之间的相互作用，节点的度值代表网络中节点与节点相连的数目。
-
为更深层次探究靶点之间的相互作用关系，将交集所得的共同靶点上传至线上软件STRING10.5 (
http://string db.org )，构建蛋白互作网络。物种选择为Homosapiens，参数保持默认设置，获取PPI网络。 -
利用R Project中的org.Hs.eg.db程序包(HYPERLINK
http://www.bioconductor.org/packages/release/data/annotation/html/org.Hs.eg.db.html )将基因symbols转变为Entrez IDs，进而利用clusterProfiler程序包(3.6版)(HYPERLINKhttp://bioconductor.org/packages http://bioconductor.org/packages/release/bioc/html/clusterProfiler.html ）对网络进行GO分析和KEGG通路分析。设定P<0.05，并按照涉及的靶点数目多少进行排序，获取GO富集分析中分子功能、细胞组成和生物过程的前15个条目以及KEGG富集的前20条信号通路。
1.1. 药物靶点筛选
1.2. 疾病靶点的收集
1.3. 构建化合物-靶点网络
1.4. 构建蛋白互作网络
1.5. 生物过程与通路分析
-
通过TCMSP检索到平贝母、大黄、黄芩、牛黄、甘草中化合物共580种：其中46个来自平贝母，92个来自大黄，143个来自黄芩，19个来自牛黄，280个来自甘草。以OB≥30%且DL≥0.18筛选出化合物158种：其中9个来自平贝母，16个来自大黄，36个来自黄芩，5个来自牛黄，92个来自甘草。入血成分144种：7个来自平贝母，10个来自大黄，32个来自黄芩，5个来自牛黄，90个来自甘草。其中，β-谷甾醇为平贝母、黄芩、甘草的共有成分。部分化合物见表1。
化合物编号 化合物名称 口服生物利用度 类药性 MOL008839 脱氧胆酸甲酯 34.63 0.73 MOL008845 脱氧胆酸 40.72 0.68 MOL008846 ZINC01280365 46.38 0.49 MOL000953 CLR 37.87 0.68 MOL001004 天竺葵素 37.99 0.21 MOL010176 堇青碱 36.81 0.71 MOL010193 七环碱 69.46 0.43 MOL010197 去氢贝母碱 52.43 0.67 MOL004440 贝母辛碱 57.40 0.81 MOL009593 去氢浙贝母碱 60.07 0.67 MOL002235 泽兰黄醇素 50.80 0.41 MOL002259 葡萄糖苷 41.65 0.63 MOL002268 大黄酸 47.07 0.28 MOL002280 Torachrysone-8-O-beta-D-(6'-oxayl)-glucoside 43.02 0.74 MOL002281 决明内酯 46.46 0.24 MOL002288 大黄素-1-O-β-D-葡萄糖苷 44.81 0.80 MOL002297 胡萝卜甙 35.89 0.70 MOL000358 β-谷甾醇 36.91 0.75 MOL000471 芦荟大黄素 83.38 0.24 MOL000096 (-)-儿茶素 49.68 0.24 MOL001689 金合欢素 34.97 0.24 MOL000173 汉黄芩素 30.68 0.23 MOL000228 （2R）-7-羟基-5-甲氧基-2-苯基铬-4-酮 55.23 0.20 MOL002714 黄芩素 33.52 0.21 MOL002909 5,7,2,5-四羟基-8,6-二甲氧基黄酮 33.82 0.45 MOL002910 红花素 41.15 0.24 MOL002913 二氢黄芩苷 40.04 0.21 MOL002914 枇杷内酯（黄烷酮） 41.35 0.24 MOL002915 丹参素 49.07 0.33 MOL002917 5,2'，6'-三羟基-7,8-二甲氧基黄酮 45.05 0.33 MOL002925 5,7,2'，6'-四羟基黄酮 37.01 0.24 MOL002927 黄芩黄酮II 69.51 0.44 MOL002928 千层纸黄素A 41.37 0.23 MOL002932 二甲氧基黄酮 76.26 0.29 MOL002933 5,7,4'-三羟基-8-甲氧基黄酮 36.56 0.27 MOL002937 二氢木蝴蝶素 66.06 0.23 MOL000525 去甲汉黄芩素 39.40 0.21 MOL000552 5,2'-二羟基-6,7,8-三甲氧基黄酮 31.71 0.35 MOL000073 表儿茶素 48.96 0.24 -
采用Cytoscape3.6.1软件构建的化合物-靶点网络包含140个化合物节点、47个靶点节点和483条边，用紫色节点表示甘草中的化合物，深绿色节点表示黄芩中的化合物，橙色节点表示大黄中的化合物，红色节点表示平贝母中的化合物，黄色节点表示共有化合物，浅绿色节点表示靶点(图1)。依据网络拓扑学性质可知，有较多节点的化合物或药物靶点在整个网络中可能起到关键作用，因此，本研究筛选节点度较大的节点进行分析。表2列出了化合物-靶点网络中的关键节点及其拓扑学性质。排名前5位的化合物分别为槲皮素（quercetin）、汉黄芩素（wogonin）、山奈酚（kaempferol）、柚皮素（naringrnin）、甘草查尔酮A（licochalcone a）。分别能与38个、13个、13个、11个、9个靶点蛋白发生作用。从靶点的角度看，度值排名前5位的是前列腺素过氧化物合酶2（PTGS2）、一氧化氮合酶（NOS2）、过氧化物酶体增生激活受体γ（PPARG）、前列腺素过氧化物合酶1（PTGS1）、丝裂原活化蛋白激酶14（MAPK14），分别能与114个、84个、72个、70个和50个化合物发生相互作用，体现了中药多靶点的特性。
名称 类别 度值 PTGS2 gene 114 NOS2 gene 76 PPARG gene 71 PTGS1 gene 57 MAPK14 gene 50 quercetin mol 38 wogonin mol 13 kaempferol mol 13 naringenin mol 11 BCL2 gene 10 CASP3 gene 10 licochalcone a mol 9 -
应用STRING软件构建PPI网络(图2)，在此网络图中包含49个节点，1 186条边，同时得到网络中关键靶点的频次(图3)。根据“度值>均值”筛选出关键节点26个，包括：MAPK3、MAPK8、IL6、CASP3、IL10、MAPK1、CXCL8、CCL2、IL1B、PTGS2、IFNG、IL4、FOS、ICAM1、IL2、MAPK14、EGFR、RELA、CASP8、CAT、STAT1、BCL2L1、HMOX1、NOS2、PPARG、CXCL10。度值最高的是MAPK3和MAPK8，度值为40；其次为IL-6，度值为39；CASP3度值为38。度值大的靶点提示在网络调控中起着关键作用，度值大的靶点很可能是金振口服液治疗NCP的关键靶点。
-
GO生物过程条目（图4）主要包括：对脂多糖的反应（response to lipopolysaccharide）、细胞内脂多糖反应（cellular response to lipopolysaccharide）、氧化应激反应（response to oxidative stress）、对机械刺激的反应（response to mechanical stimulus）以及外源性凋亡信号通路（extrinsic apoptotic signaling pathway）等；分子功能包括：细胞因子受体结合（cytokine receptor binding）、细胞因子活性（cytokine activity）、受体-配体活性（receptor ligand activity）、趋化因子受体结合（chemokine receptor binding）以及趋化因子活性（chemokine activity）；并且参与了膜筏（membrane raft）、膜微区（membrane microdomain）、膜区（membrane region）结构的组成。KEGG富集信号通路（图5，图6）主要包括：AGE-RAGE信号通路（AGE-RAGE signaling pathway in diabetic）、IL-17信号通路（IL-17 signaling pathway）、甲型流感信号通路（Influenza A）、TNF信号通路（TNF signaling pathway）、南美锥虫病信号通路（Chagas disease）等。
2.1. 筛选药物活性成分
2.2. 化合物-靶点网络构建
2.3. PPI网络的构建与分析
2.4. GO富集与KEGG富集分析
-
在这场防治疫情保卫战中，中医药发挥了积极作用，大量中药方剂已经应用于临床并显示出其明确的疗效[7-8]。中医学认为新冠肺炎属于疫病的范畴。病因是湿热戾气，病机是由于湿热戾气从口鼻而入，首先犯肺，使肺的宣发肃降功能失常，从而出现咳嗽；因为邪气性质湿热，湿邪易阻遏气机损伤阳气，热邪易伤津耗气，故而患者会出现乏力、发热的表现[9-10]。随着疾病逐渐发展，邪气按卫气营血的规律传变入里，最终逆传心包，出现神昏谵语等危重证候。
现代医学认为新冠肺炎是由SARS-CoV-2导致的，故而抗病毒是治疗新冠肺炎的重要举措，临床研究证实，金振口服液可以治疗呼吸道合胞病毒、严重呼吸窘迫综合征冠状病毒、柯萨奇病毒等引起的上呼吸道感染[11-13]，但物质基础与作用机制尚不明确。本研究发现金振口服液抗炎的活性物质主要为黄酮类化合物，包括槲皮素、汉黄芩素、山奈酚、柚皮素、甘草查尔酮A等。现代研究表明：槲皮素可与病毒表面蛋白相互作用阻断病毒的黏附[14]；汉黄芩素可以通过抑制病毒基因复制和蛋白合成抑制病毒增殖[15]；柚皮素可以抑制神经氨酸酶活性与血凝素活性发挥抗病毒的作用[16]。KEGG富集分析结果表明，金振口服液抗SARS-CoV-2的机制可能是通过作用于甲型流感病毒信号通路中PKC、ERK、STAT1等靶点实现的（图7）。
基于目前的临床调查，现代医学认为新冠肺炎可以通过自身的免疫应答导致肺上皮细胞短时间内分泌大量的炎性细胞因子，这些炎性细胞因子会使得炎症细胞大量渗入肺部，引起免疫激化。新冠肺炎的严重程度与机体炎症反应的强度相关，随着疾病发展，如果患者体内出现过度的炎症反应，即细胞因子风暴，将可能导致肺上皮细胞和内皮细胞大量凋亡，最终引起急性肺损伤与呼吸窘迫综合征等危重证候[17-19]。
现代医学应用糖皮质激素来控制细胞因子风暴，但糖皮质激素的副作用使得抗疫陷入进退两难的境地。中药在抗炎作用上具有明显优势，在中医理论的指导下，发挥整体调节的优势，并根据药物配伍减毒增效，是目前治疗新冠肺炎的重要手段。PPI结果表明：金振口服液中的活性物质主要作用于炎性介质，其可能是通过降低白细胞介素、肿瘤坏死因子等炎症因子的分泌水平，调节巨噬细胞的极化等发挥抗炎作用[20-21]。现代研究表明：槲皮素具有抑制炎症反应、抗呼吸道感染的作用[22-23]；汉黄芩素具有抗氧化活性、抗炎等作用[24]；山奈酚可以通过抑制肺组织促炎因子减轻脂多糖导致的小鼠急性肺损伤[25]。KEGG富集分析结果表明其抗炎机制可能是通过作用于IL-17信号通路和TNF信号通路实现的（图8、图9）。实验研究表明：金振口服液可以通过抑制p65、NF-κB、ERK1/2、p38多个靶蛋白磷酸化，从而阻断核转录因子-κB，丝裂原活化蛋白激酶炎症通路，降低炎症细胞因子的释放，从而改善急性肺损伤[26-27]，证实了网络药理学的结果。
综上所述，金振口服液治疗新冠肺炎的有效成分可能是槲皮素、汉黄芩素、山奈酚、柚皮素、甘草查尔酮A等，这些成分通过炎症相关靶点参与IL-17、TNF、甲型流感等信号通路，从而发挥抗病毒、抗炎的作用。因此，在金振口服液的生产和质量控制过程中，应当重视槲皮素、汉黄芩素等主要活性物质。此外，由于石膏、青礞石和羚羊角作为矿物类和动物类药材的特殊性，笔者未能在数据库中获取其活性成分，故在本研究结论中未能体现其作用。但中医学认为生石膏辛甘大寒可解肌退热；青礞石咸平可坠痰下气；羚羊角咸寒可散血解毒，三者合用共奏清热止咳、宣肺涤痰之功。临床上应用的银翘白虎汤、礞石滚痰丸、羚羊角颗粒在退热、止咳方面疗效显著[28-30]。现代药理学研究表明：石膏具有一定的抗炎效应，可能是通过抑制TLR4/NF-κB通路实现的[31]；青礞石可以通过抑制海马内胶质纤维酸性蛋白的表达发挥镇静作用[32]；羚羊角可降低炎症因子如前列腺素、白细胞介素等的表达起到解热作用[33]。因此，对于矿石药以及动物药的活性成分及作用机制也应深入研究，以期为传统中药的开发和利用提供理论依据；同时，药物活性成分在炮制过程中改变、中药活性成分含量不确定等因素，使得本研究结果可能与临床效果存在一定偏差。但是，网络药理学的研究结果也可为中医药治疗新冠肺炎提供一定的理论基础，为临床上用药选择提供一定依据。
致谢
感谢金永生老师在网络药理学和活性化合物收集等方面给予的技术支持。
DownLoad:
