基于网络药理学探讨扶正化瘀方治疗肝癌的作用机制

齐婧姝， 平大冰， 孙 鑫， 薛静波， 陶艳艳， 彭 渊， 刘成海,,3,4

1 上海中医药大学附属曙光医院 肝病研究所， 上海 201203； 2 上海市中医临床重点实验室， 上海 201203；3 肝肾疾病病证教育部重点实验室， 上海 201203； 4 上海中医健康服务协同创新中心， 上海 201203

原发性肝癌(简称肝癌)包括肝细胞癌和肝内胆管细胞癌[1]，其病死率居全世界恶性肿瘤第三位[2]，我国每年肝癌死亡人数占全球肝癌死亡的51%[3]。肝癌往往起病隐匿、进展迅速、复发快、预后差[4]，是危害人类健康的重大疾病。因此，对肝癌的防治是临床研究的重点。

扶正化瘀方(fuzheng huayu，FZHY)由丹参、冬虫夏草、桃仁、松花粉、绞股蓝、五味子6味药组成，适用于乙型肝炎肝纤维化属瘀血阻络，肝肾不足证者。本课题组前期工作[5]表明，该方可有效改善肝炎后肝硬化患者肝功能、改善肝组织炎症，并逆转肝纤维化；方中多味药及其不同化学成分在肝纤维化各个病理环节融对抗、保护、调节于一体，标本兼治。已有临床和基础研究表明，FZHY可以降低肝癌的恶性程度[6]，防治实验性肝癌的发生率[7-9]，然而其防治肝癌的活性成分以及潜在的作用机制尚不清楚。因此，基于中药复方多成分、多靶点、多途径发挥药理机制的作用特点，本研究以FZHY为载体，运用网络药理学方法，从分子网络水平系统探索该方多成分、多靶点与肝癌疾病间复杂网状关系，为后续实验研究提供线索。通过整合多个数据库信息收集整理FZHY的主要化学成分及其靶标信息、肝癌相关基因信息，构建FZHY治疗肝癌的成分-靶标网络，旨在通过网络分析解析FZHY治疗肝癌的作用机制，为临床和基础科研提供理论依据。

1 材料与方法

1.1 搜索FZHY主要活性成分的作用靶点 FZHY中各味药的成分信息从TCMSP(https://tcmspw.com/tcmsp.php，数据更新至2020-12-20)、BATMAN(http://bionet.ncpsb.org/batman-tcm/，数据更新至2020-12-20)、Drugbank(https://go.drugbank.com/，数据更新至2020-12-20)等数据库中检索获得；各成分的靶标信息从TCMSP、Drugbank中获取，构建FZHY有效成分潜在的靶标群。方中各味药的归经信息从《中国药典》(2015年版)中获取。TCMSP数据库筛选过滤标准为：口服生物利用度(oral bioavailable,OB)≥30%且类药性(drug likeness,DL)≥0.18，获得FZHY多成分靶标基因。

1.2 提取肝癌相关基因 通过检索“Liver cancer”“Hepatic carcinoma”“Hepatoma”“Hepatocellular carcinoma”等关键词，从GeneCards人类基因数据库(https://www.genecards.org，数据更新至2020-12-20)、OMIM在线人类孟德尔遗传数据库(Online Mendelian Inheritance in Man,OMIM, https://omim.org，数据更新至2020-12-20)中搜索肝癌相关的基因，构建肝癌致病基因的靶点数据集，获得肝癌疾病靶标基因。

1.3 FZHY抗肝癌靶标网络的构建与分析 对FZHY多成分靶标基因和肝癌疾病靶标基因进行基因去重，绘制韦恩图，获得FZHY靶标群与肝癌的共有基因，并在STRING数据库(http://sring-db.org/)构建蛋白-蛋白互作关系(protein-protein interaction，PPI)分析。PPI分析模式限定为“Multiple proteins”，物种设定为“Homo Sapiens”，设置可信度≥0.7，隐藏孤立的蛋白，对数据进行分析后导出PPI网络图，再利用Cytoscape软件建立多成分-基因-疾病网络模型。

1.4 GO功能富集分析和KEGG通路富集分析 利用R语言中的“colorspace”“stringi”“ggplot2”“clusterProfiler”“enrichplot”“org.Hs.EG.db”包对FZHY有效成分潜在作用靶点进行功能注释(gene ontology,GO)，即主要进行分子功能(molecular function，MF)、生物过程(biological process，BP)和细胞组成(cellular component，CC)3个部分的富集分析，并对相关基因进行KEGG通路富集分析，设定显著性为P<0.05。

1.5 成分-靶点-通路网络 结合FZHY抗肝癌的主要成分、核心靶点与KEGG分析所得通路，运用Cytoscape软件绘制成分-靶点-通路网络图，分析其整体关联。

2 结果

2.1 FZHY包含的成分及对应靶标 FZHY由6味中药组成，根据OB≥30%和DL≥0.18筛选出中药包含的主要成分共142个。合并各成分的靶标基因，除去被多个成分共同调控的重复基因，获得197个药物靶标基因(表1)。

2.2 FZHY的靶标网络分析 分别以检索词“Liver cancer”“Hepatic carcinoma”“Hepatoma”“Hepatocellular carcinoma”在GeneCards和OMIM数据库中检索肝癌靶标基因，分别获得相应疾病相关基因数(表2)。合并去除重复基因，共获得疾病相关基因17 869个。

表2 在线数据库检索关键词及疾病相关基因数

2.3 FZHY成分-基因-肝癌疾病靶标基因网络分析 对FZHY多成分靶标基因和肝癌疾病靶标基因进行基因去重，绘制韦恩图，获得FZHY靶标群与肝癌靶标群的共有基因192个(药物-疾病靶标基因)(图1)。利用Cytoscape软件构建FZHY有效成分-肝癌靶标基因网络图(附录1)，得到FZHY治疗肝癌的潜在化合物有95种，其中与肝癌靶标基因连线最多为槲皮素和木犀草素(Mol ID: MOL000098、MOL000006)，提示槲皮素和木犀草素是FZHY调控肝癌的重要活性成分。PPI网络分析药物-疾病靶标基因间的关系(附录2)，根据基因-基因连线数分析PPI网络核心基因，显示主要为INS、IL-6、EGFR等基因(图2)，提示这些靶标能被FZHY多味中药和多个成分共同调控，因此在靶标网络中可能具有更重要的影响。

表1 FZHY包含化学成分的对应靶标数

图1 基于GeneCards和OMIM数据库预测的肝癌与FZHY的靶点基因交集

图2 FZHY治疗肝癌的关键靶标基因

2.4 网络模型中重要靶标生物功能分析 针对核心靶点进行GO-BP分析，确定靶蛋白具有显著性意义(P<0.05)的生物学过程主要涉及：对药物的反应、对氧水平的反应、对类固醇激素的反应、对缺氧的反应、对金属离子的反应、对氧水平降低的反应、对氧化应激的反应、细胞对化学应激的反应、对酮的反应、对营养物质的反应。GO-CC分析结果提示涉及的细胞组分主要有：突触后膜、膜筏、膜微域、膜区、突触膜、囊腔、分泌颗粒腔、细胞质囊腔、突触膜整体成分、突触后膜整体成分等。GO-MF分析结果提示涉及的分子功能主要有：DNA-结合转录因子结合、羧酸结合、有机酸结合、类固醇结合、神经递质受体活性、核受体活性、配体-激活转录因子活性、突触后神经递质受体活性、半胱氨酸-类型内肽酶活性参与凋亡过程、半胱氨酸-类型内肽酶活性参与凋亡信号通路等(附录3)。

对FZHY靶标网络中的关键基因进行KEGG信号通路富集分析，发现前50位信号通路有22条与肝癌密切相关，具体见表3和附录4。

用Cytoscape软件构建FZHY治疗肝癌关键靶标与关键信号通路的网络图(附录5)，发现肝癌22条关键信号通路与FZHY治疗肝癌的80个潜在靶点关系密切。

表3 FZHY治疗肝癌关键靶标的KEGG富集通路

3 讨论

在我国，肝癌是肝炎后肝硬化最常见的并发症[10]，给患者和社会带来沉重的负担。尽管外科手术切除及肝移植、射频消融、介入治疗以及个性化化疗方案和靶向药物的应用使得肝癌患者获益良多，但复发仍是临床所面临的重大难题。对于难以手术切除肝癌的综合治疗方案，中西结合治疗是提高肝癌疗效的一个重要途径[11-12]。FZHY(制剂为胶囊或片剂)由丹参、桃仁、虫草菌丝、松花粉、五味子、绞股蓝组成，对肝炎后肝硬化具有良好的治疗效果[13]，已有国内同行发表FZHY抗肝癌的相关临床与基础研究，但对于其抗肝癌的活性成分、作用靶点以及分子机制尚不完全明确。

本研究采用网络药理学方法从分子网络角度系统分析了FZHY对肝癌的潜在治疗机制，通过TCMSP平台，筛选得到FZHY活性成分142种，药物靶点基因197个；将FZHY对应的靶点与肝癌的疾病靶点基因取交集，获得交集基因192个，筛选得到潜在化合物95种，其中部分化合物抗肝癌的有效性及其药理机制已被证实，起到重要调控作用的活性成分有：槲皮素、木犀草素等。据报道，槲皮素抗肝癌的作用可能与其灭活caspase/Cyto-c通路、抑制AP-2β/hTERT、NF-κB/COX-2、阻断Akt/ERK1/2信号通路及竞争性抑制葡萄糖转运蛋白(GLUT-1)密切相关[14-15];木犀草素可通过TGFβ1、p53和Fas/FasL三条信号通路抑制肝癌细胞生长[16]。上述化合物位于FZHY治疗肝癌疾病-活性成分-靶基因网络的核心，可能在FZHY抗肝癌过程中起到重要作用。

进一步对活性成分靶点-信号通路网络进行分析，得到80个抗肝癌靶点，涉及22条信号通路。FZHY潜在作用的关键靶点有INS、IL-6、EGFR等；主要通过调控细胞凋亡、TNF等信号通路而抑制肝癌的发生和发展，并促进肝癌细胞凋亡等多途径、多靶点发挥作用。研究显示，胰岛素可促进肝癌细胞增殖和侵袭，抑制细胞的凋亡，其功能与 Cleaved caspase 3、MMP-2、MMP-9 蛋白表达及 PI3K/AKT 信号通路密切相关[17]；IL-6在人肝细胞癌中的表达显著上调，与肝癌恶性进展呈正相关[18]，其对肝癌的正向调节作用可能与其抑制miRNA-197的表达来调节STAT3总蛋白质的表达，继而促进p-STAT3的表达关系密切[19]；文献报道，在29%～85%的肝癌组织中有EGFR过表达的现象[20]，甚至有专家建议将EGFR作为肝癌治疗的靶点[21]，其机制可能为EGF/EGFR通过PI3K/AKT信号通路磷酸化p300 Ser1834位点，促进组蛋白H3K9位点乙酰化，使DKK1的转录增加，进而导致肝癌细胞的转移[22]。

通过GO功能富集分析发现，FZHY对于肝癌的调控主要集中在DNA-结合转录因子结合、核受体结合等过程，提示FZHY可能通过抑制癌细胞增殖，诱导癌细胞凋亡而发挥作用。通过KEGG通路富集分析追踪FZHY治疗肝癌的代谢通路，并作FZHY治疗肝癌关键靶标的关键信号通路-基因网络图分析整体联系。结果发现，FZHY主要作用于细胞凋亡、TNF等信号通路。其中凋亡信号通路可以影响胰岛素的产生和释放[23-24]，从而抑制肿瘤细胞增殖和侵袭能力，促进其凋亡。由此可知，FZHY可能可以通过上述途径发挥对肝癌的调控作用。

综上所述，本研究运用网络药理学技术手段探讨了FZHY活性成分与肝癌疾病及其相关基因的网状关联，通过研究该方多成分、多靶点、多途径作用机制，为进一步深入开展的实验研究提供了参考和线索。本课题组研究将围绕其抗肝癌药效物质及作用机制开展进一步基础与临床转化研究，以期开发出更适用于临床治疗肝癌的中药制剂。

利益冲突声明：本研究不存在研究者、伦理委员会成员、受试者监护人以及与公开研究成果有关的利益冲突。

作者贡献声明：齐婧姝负责撰写修改论文；平大冰、孙鑫、薛静波参与文献的收集整理；彭渊、陶艳艳负责拟定写作思路及制作网络图；刘成海负责指导撰写文章并最后定稿。

附录1～5见二维码