2020年《全球癌症统计报告》表明，原发性肝癌的发病率排名第六，死亡率排名第三，严重影响了人们的生活质量^[1-2]。中医药治疗肿瘤历史悠久、疗效独特，是我国肝癌患者综合治疗方案的重要组成部分之一。传统中医药学与现代术疗、化疗等西医疗法相结合，能够起到减轻症状，增效减毒，增强免疫的效果，提高患者的生存率及生活质量^[3-4]。

百合地黄汤为中药经典名方之一，出自汉代张仲景《金匮要略方论•百合狐惑阴阳毒病证治篇》^[5-6]，含“百合七枚(擘)，生地黄汁一升”，始载用于“百合病，不经吐、下、发汗，病形如初者”。百合归心、肺经，具有养阴润肺，清心安神之功效，地黄归心经、肝经、肾经，有清热凉血，生津可治阴虚发热等功效。现代药理学研究发现，百合地黄汤高剂量组对肝癌H22荷瘤小鼠有抑瘤作用，且抑瘤作用呈剂量依赖关系^[7]。百合地黄汤加黄芪能降低肝癌小鼠PCNA蛋白及Cyclin D1蛋白的表达，可抑制肝肿瘤和Lewis肺细胞增殖，有较好的抗肺肿瘤转移作用^[8]。目前对百合地黄汤抗癌作用报道较少，有待深入研究。

现代信息的发展有助于对中国传统中医药领域进行抗肿瘤药物的数据挖掘。网络药理学^[9-13]是多向药理学、计算生物学、大数据分析等学科融合而成的，体现了中医整体观、辨证论治、方剂配伍的理论，为从生物学分子网络的结构和功能来认知中药药性提供了方法学契机和有力工具。本文通过网络药理学构建多层网络，预测百合地黄汤抗肝癌信号通路和作用机制，以期推动其在在肿瘤防治方面的认识。

利用文献和TCMSP中药系统药理分析平台^[14]( https://old.tcmsp-e.com/tcmsp.php)检索百合和地黄中的化学成分及相关参数。主要通过文献收集地黄中化学成分，通过TCMSP数据库收集百合中化学成分，借助PubChem数据库^[15]( https://pubchem.ncbi.nlm.nih.gov)获取其分子结构，分子结构统一保存为“sdf”格式，进一步导入SwissADME数据库^[16]( http://www.swissadme.ch/index.php)中获得百合、地黄的GI值(血糖生成指数)、Lipinski值(类药五原则)、OB值(口服生物利用度)。

借助SwissTargetPrediction数据库^[17]( http://www.swisstargetprediction.ch)对收集的成分进行靶点匹配和预测，得到百合地黄汤化学成分的作用靶点。在OMIM^[18]( https://omim.org)、GeneCards数据库^[19]( https://www.genecards.org)，以“Hepatocellular Carcinoma”为关键词进行检索，得到肝癌相关基因并导出结果，借助Uniprot数据库^[20]( https://www.uniprot.org)将靶点蛋白名称标准化。将百合、地黄化学成分作用靶点、肝癌相关靶点导入在线韦恩图网站( http://bioinformatics.psb.ugent.be/webtools/Venn)，获得百合地黄汤作用靶点与肝癌的交集靶点，即百合地黄汤抗肝癌的潜在作用靶点。

为了在蛋白水平上进一步了解百合地黄汤与疾病靶点的作用机制，将百合地黄汤抗肝癌的潜在作用靶点输入到STRING数据库( https://www.string-db.org)中，物种限定为“homo sapiens”，置信度得分高于0.7，收集获得的蛋白−蛋白相互作用信息，将其导入Cytoscape3.7.2软件^[21]中。构建蛋白相互作用图，利用network analysis计算每个节点的打分，筛选打分较高的节点，从而得出PPI网络中的核心基因。

以百合地黄汤化学成分、作用靶点为数据来源，采用Cytoscape 3.7.2软件构建“药物−化学成分−作用靶点”网络。网络中的节点为药物、化学成分、作用靶点。将百合地黄汤的潜在抗肝癌靶点与化合物做映射，获得活性化合物。在Cytoscape 3.7.2软件中建立“药物−活性成分−抗肝癌靶点”网络图，借助network analysis进行拓扑学分析，获得每个化学成分的自由度，从而得出百合地黄汤抗肝癌的核心成分。

为了进一步了解上述选择的核心基因的功能及其在信号通路中的作用，将筛选的百合地黄汤抗肝癌的潜在作用靶点导入David数据库^[22]( https://david.ncifcrf.gov/)，并将物种限制在人类，所有的基因名称都被更正为它们的官方名称。在上述数据库搜索和转换操作，设定阈值P<0.05，进行GO生物过程富集分析和KEGG信号通路富集分析，选取前15的通路通过在线绘图网站OMishare Tools( https://www.omicshare.com/tools/)将结果可视化。

在Pubchem数据库中下载百合地黄汤抗肝癌的核心成分的3D分子结构，在PDB数据库( https://www.rcsb.org/)中下载PPI网络中的核心基因的靶点蛋白。通过AutoDock Tool 1.5.6软件进行分子对接前处理，如去水、确定分子对接坐标等，再利用AutoDock Vina 4.2软件^[23]对PPI网络中的核心基因和百合地黄汤抗肝癌的核心成分进行分子对接。运用PyMol 2.4软件展示得到的分子对接结果，立体表现出化合物配体在蛋白靶点中的情况。

经过数据库、文献查询，查找出百合地黄汤中共含有化合物131个。为了更全面地筛选活性成分，将收集到的百合地黄汤中所有化合物进行分析考察。表1展示了部分百合地黄汤化合物的GI值、Lipinski值、OB值。

表 1(Table 1)

表 1 百合地黄汤部分化合物信息 Table 1 Information of some compounds of Baihe Dihuang Decoction 化合物 GI值 Lipinski值 OB值 2-甲氧基-4-乙烯基苯酚 高 满足 0.55 豆甾醇糖苷 高 MW>500 0.55 橙花基乙酸酯 高 满足 0.55 正丁基-β-D-果糖苷 高 满足 0.55 邻苯二甲酸二丁酯 高 满足 0.55 地黄甙C 高 满足 0.55 氟哌啶醇苷B 高 满足 0.56 2-丁基硼酸 高 满足 0.56 8-表庚酸 高 满足 0.56 阿魏酸 高 满足 0.56 香草酸 高 满足 0.56 山奈酚 高 满足 0.55 26-O-β-D-glucopyranosyl-nuatigenin-3-O-α-L- rhamnopyranosyl (1→2)-β-D-glucopyranoside_qt 高 满足 0.55 26-O-β-D-Glucopyranosyl-3β,26-dihydroxy- choleslen-16,22-dioxo-3-O-α-L-rhamnopy- ranosyl-(1-2)-β-D-glucopyranoside_qt 高 满足 0.55 百合皂苷 低 MW>500, NorO>10, NHorOH>5 0.11 百合苷C 低 NHorOH>5 0.55 梓醇 低 NHorOH>5 0.55 毛蕊花糖苷 低 MW>500, NorO>10, NHorOH>5 0.17 齐敦果酸 低 MLOGP>4.15 0.56

表 1 百合地黄汤部分化合物信息 Table 1 Information of some compounds of Baihe Dihuang Decoction

经SwissTargetPrediction数据库预测，得到百合作用靶点711个，地黄作用靶点454个，核查数据后删除重复项，共得到药物靶点360个。在OMIM、GeneCards数据库中，经筛选共得到肝癌相关靶点488个。将360个药物靶点与488个疾病靶点导入在线韦恩图网站，得到58个交集基因(见图1)。这58个交集基因即为百合地黄汤抗肝癌作用靶点，其中百合、地黄、肝癌都对应的靶点有MMP2、KDR、EGFR、ALK、IL2、LGALS3、CDK1、PPARG、ABCB1、STAT3、F2、CA9、ESR1、HRAS、AR、ESR2和MMP9。

图 1(Figure 1)

图 1 百合地黄汤与肝癌对应靶点匹配韦恩图 Figure 1 Venn diagram of Baihe Dihuang Decoction and corresponding targets of HCC

由肝癌相关基因与药物靶点基因映射筛选出共同靶点58个，构建的PPI网络包括58个节点，821条边，见图2。图中节点的大小依自由度而设定，节点之间的关系强弱由边的线条粗细表示，边的线条越粗则两蛋白的Combine score(连接分数)值越高。采用Cytoscape 3.7.2软件进行网络拓扑学分析，结果表明，间距中心度、自由度及紧密中心度大于均值(间距中心度为0.009，自由度为28.31，紧密中心度为0.68)的关键节点有15个，占总节点数的25%。表2展示了百合地黄汤抗肝癌靶点蛋白质相互作用靶点自由度值排名，其中CASP3、VEGFA、STAT3与其他靶点的连接最多，节点最大，这说明它们与其他靶点蛋白相互作用更强，在该网络中发挥了关键作用。

图 2(Figure 2)

图 2 百合地黄汤抗肝癌靶点蛋白质相互作用图 Figure 2 Anti-HCC target protein interaction diagram of Baihe Dihuang Decoction

表 2(Table 2)

表 2 百合地黄汤抗肝癌靶点蛋白质相互作用靶点自由度值排名 Table 2 Ranking of the degree of target protein interaction of Baihe Dihuang Decoction against HCC 排名 靶点 自由度值 排名 靶点 自由度值 1 VEGFA 52 6 ALB 49 2 STAT3 51 7 SRC 48 3 CASP3 51 8 HRAS 48 4 EGFR 50 9 ESR1 47 5 AKT1 50 10 JUN 46

表 2 百合地黄汤抗肝癌靶点蛋白质相互作用靶点自由度值排名 Table 2 Ranking of the degree of target protein interaction of Baihe Dihuang Decoction against HCC

百合地黄汤药物−化学成分−作用靶点网络图见图3，共有499个节点和1319条边，包括131个化合物节点，平均自由度是5.286。其中“BH”红色三角节点代表百合药物，“DH”红色三角节点代表地黄药物，“DH+数字”蓝色椭圆节点代表地黄化学成分，“BH+数字”紫色椭圆节点代表百合化学成分，绿色椭圆节点代表作用靶点。由图3可知，同一活性成分可以作用于不同靶标，同一靶标也可对应不同活性成分，体现出百合地黄汤抗肝癌具有多成分、多靶标的特点。进一步分析发现，百合地黄汤抗肝癌的靶点可以映射到87个化合物，这87个化合物即为百合地黄汤抗肝癌活性化合物。“药物−活性成分−抗肝癌靶点”网络图见图4，“BH”红色椭圆节点代表百合药物，“DH”红色椭圆节点代表地黄药物，“BH+数字”椭圆蓝色节点代表百合活性成分，“DH+数字”紫色椭圆节点代表地黄活性成分，绿色椭圆节点代表抗肝癌靶点。边线代表药物−活性成分、活性成分−抗肝癌靶点的相互关联。网络拓扑学分析得到了各个节点之间的自由度，表3展示了“药物−活性成分−抗肝癌靶点”网络图中活性成分与靶点的自由度排名前10的活性成分，经统计自由度值较高的活性成分有山奈酚、百合皂苷、百合苷C、26-O-β-D-glucopyranosylnuatigenin-3-O-α-L-rhamnopyranosyl(1→2)-o-[β-D-glucopyranosyl(1→4)]-β-D-glucopyranoside_qt。

图 3(Figure 3)

图 3 百合地黄汤药物−化学成分−作用靶点网络图 Figure 3 Network diagram of drug-chemical composition-targets of Baihe Dihuang Decoction

表 3(Table 3)

表 3 百合地黄汤药物−活性成分−抗肝癌靶点网络自由度值排名 Table 3 Ranking of degree value of drug-active ingredient-anti-HCC target network of Baihe Dihuang Decoction 排名 代称 化合物 自由度值 1 DH40 kaempferol 19 2 BH70 26-O-β-D-glucopyranosylnuatigenin-3-O-α-L-rhamnopyranosyl(1→2)-o- [β-D-glucopyranosyl(1→4)]-β-D-glucopyranoside_qt 6 3 BH65 brownioside 5 4 BH15 lilioside C 5 5 BH68 27-O-(3-hydroxy-3-methylglutaroyl)-isonarthogenin-3-O-α-L-rhamnopyranosyl- (1→2)-O-[β-D-glucopyranosyl-(1→4)]-β-D-glucopyranosid 5 6 BH6 Stigmasta-5,22-dien-3-O-β-D-glucopyranoside 4 7 BH62 darutoside 4 8 BH67 27-O-(3-hydroxy-3-methylglutaroyl)-isonarthogenin-3-O-α-L-rhamnopyranosyl- (1→2)-O-[β-D-glucopyranosyl-(1→4)]-β-D-glucopyranosid_qt 4 9 BH75 26-O-β-D-glucopyranosyl-3β,26-dihydroxy-5-cholesten-16,22- dioxo-3-O-α-L-rhamnopyranosyl-(1→2)-β-D-glucopyranoside 4 10 BH13 n-butyl-β-D-fructoufranoside 4

表 3 百合地黄汤药物−活性成分−抗肝癌靶点网络自由度值排名 Table 3 Ranking of degree value of drug-active ingredient-anti-HCC target network of Baihe Dihuang Decoction

图 4(Figure 4)

图 4 百合地黄汤药物−活性成分−抗肝癌靶点网络图 Figure 4 Network diagram of drug-active ingredient-anti-HCC targets of Baihe Dihuang Decoction

对百合地黄汤抗肝癌作用所涉及的58个蛋白靶点进行GO功能富集分析和KEGG信号通路富集分析，以P<0.05为筛选条件，得到富集基因数量较多的GO功能和信号通路。在KEGG富集分析中共涉及91条富集通路，包括了癌症蛋白多糖、癌症通路、PI3K-Akt信号通路、VEGF信号通路、癌症的中心代谢、胰细胞、肾癌、Rap1信号通路、H1F信号通路等。GO功能富集分析中确定了388个GO条目。生物过程相关的条目286个，涉及药物反应、细胞增殖及调控、蛋白质磷酸化、基因表达及调控、细胞迁移调控、信号传导等。分子功能方面相关的条目68个，涉及蛋白质结合、蛋白络氨酸激酶活性、蛋白激酶活性、ATP结合酶活性、相同蛋白质结合、类固醇激素受体、转录因子结合等。细胞组成相关的条目34个，涉及细胞器、质膜、突触等方面。图5列出了GO功能富集分析的前24个条目，图6列出了信号通路富集分析的前15个条目。

图 5(Figure 5)

图 5 百合地黄汤百合地黄汤抗肝癌GO功能富集分析 Figure 5 GO function enrichment analysis of anti-HCC targets of Baihe Dihuang Decoction

图 6(Figure 6)

图 6 百合地黄汤抗肝癌靶点的信号通路富集 Figure 6 Enrichment diagram of anti-HCC signal pathway of Baihe Dihuang Decoction

根据上述筛选得到的87个活性成分，选取活性成分中自由度值前4的化合物(山奈酚、百合皂苷、百合苷C、26-O-β-D-glucopyranosylnuatigenin-3-O-α-L-rhamnopyranosyl(1→4)-o- [β-D-glucopyranosyl(1→4)]-β-D-glucopyranoside_qt)与位于前3位的核心基因(STAT3、CASP3、VEGFA)使用AutoDock Vina软件进行分子对接，获得活性化合物与靶标相互作用的最小结合能。分子对接结果见表4，结合能都小于−5 kJ·mol⁻¹，可见百合地黄汤与这3个肝癌蛋白靶点具有较好的亲和性，图7展示了活性化合物与对接结合能最低的靶点的分子对接情况。

表 4(Table 4)

表 4 分子对接结果 Table 4 Molecular docking results 化合物 结合能/(kJ·mol−1) CASP3 STAT3 VEGFA 山奈酚 −7.4 −7.4 −6.4 26-O-β-D-glucopyranosylnuatigenin-3-O-α-L-rhamnopyranosyl (1→2)-o-[β-D-glucopyranosyl(1→4)]-β-D-glucopyranoside_qt −8.4 −8.3 −7.9 百合皂苷 −8.2 −8.5 −9.1 百合苷C −5.5 −5.3 −5

表 4 分子对接结果 Table 4 Molecular docking results

图 7(Figure 7)

图 7 百合地黄汤中活性成分与关键靶点的分子对接图 Figure 7 Molecular docking diagram of active components and key targets in Baihe Dihuang Decoction

本文进行了百合地黄汤中抗肝癌主要活性成分及潜在作用机制的探讨。百合地黄汤中含有131种化合物成分。在前期摸索中，采用文献报道常用筛选标准(OB≥30%，DL≥0.18)对化学成分进行筛选^[24]，仅得到8个化合物，其中大部分含有较高药理活性的化合物如百合皂苷、百合苷C、梓醇、毛蕊花糖苷、梓醇、齐墩果酸等被筛除掉了。鉴于该标准的局限性，笔者综合考量，将百合地黄汤的化合物都进行分析。本文经网络药理学的方法筛选得到了87个活性成分及58个抗肝癌作用靶点，其中山奈酚、百合皂苷、百合苷C、26-O-β-D-glucopyranosylnuatigenin-3-O-α-L-rhamnopyranosyl(1→2)-o-[β-D-glucopyranosyl(1→4)]-β-D-glucopyranoside_qt等活性成分对应靶点较多，大多为皂苷类化合物。中药皂苷类成分可通过抑制肿瘤细胞增殖、诱导肿瘤细胞凋亡、抑制侵袭转移、抑制肿瘤血管新生、逆转肿瘤多药耐药、调节机体免疫、诱导肿瘤细胞自噬死亡等机制发挥抗肿瘤作用^[25-27]。现代药理研究表明，山奈酚通过调控环氧合酶-2的转录、抑制P-糖蛋白的功能和表达、抑制白细胞介素-4抗细胞增殖来诱导细胞凋亡、干预细胞信号转导和增强抑癌基因活性及抑制癌基因表达^[28]。

PPI网络显示百合地黄汤潜在作用的关键靶点有STAT3、CASP3、VEGFA。STAT3在细胞恶变过程中起关键性作用，肿瘤的增殖、细胞凋亡抑制、侵袭和转移、新血管生成以及肿瘤干细胞的自我更新和分化会导致其高表达从而减缓细胞恶变^[29-31]。CASP3是与免疫细胞凋亡密切相关的半胱氨酸天蛋白酶，能促进免疫细胞发挥作用^[32]。Folkman^[33]提出了肿瘤性血管生成理论，肿瘤的发生和发展依赖于肿瘤血管的生成。VEGFA是血管内皮生长因子的一种，可以通过与VEGFA2(血管内皮生长因子受体)结合，影响血管内皮细胞增殖、存活、迁移及血管通透性，进而抑制肿瘤细胞的增殖^[34]。

信号通路富集表明百合地黄汤主要是调控癌症蛋白聚糖、VEGF信号通路和PI3K-Akt通路而抑制肿瘤细胞增殖、促进肿瘤细胞凋亡等发挥抗肝癌作用。现代药理研究表明^[27]，肝癌发展进程受多条信号通路传导调控，中医药可同时干预多条信号通路，多靶点的调控转导因子表达，从而抑制、杀伤肿瘤细胞，改善症状和体征，发挥防治肝癌的作用。肿瘤微环境中的许多蛋白聚糖(PGs)已被证明是关键的大分子，这些大分子有助于各种类型的癌症的生物学发展，包括增殖、黏附、血管生成和转移，从而影响肿瘤的进展，通过调节细胞增殖和凋亡、调节转录等功能来治疗癌症^[35]。Napoleone Ferrara等^[36]研究发现，VEGF通过与肝癌细胞表面的结合，激活Akt激酶的相关信号通路而促进肿瘤细胞的迁移。Zou Y等^[37]研究证明在多种肿瘤干细胞中，PI3K-Akt信号通路的激活发挥重要作用，抑制PI3K-Akt信号的激活可以抑制肿瘤干细胞的自我更新，并减少肿瘤干细胞的比例。由此预测，百合地黄汤通过以上通路抑制肿瘤细胞增殖，防治肿瘤细胞复发转移，同时调节机体免疫功能，提高机体免疫细胞杀伤肿瘤细胞的能力。分子对接结果显示，百合地黄汤的4个核心活性成分与肝癌的3个核心靶点(STAT3、CASP3、VEGFA)具有良好亲和活性。

综上所述，本文以百合地黄汤主要活性成分为研究对象，通过网络药理学技术分析了其潜在的抗肝癌靶点及信号通路，并结合分子对接方法验证了百合地黄汤主要活性成分能与肝癌靶点发生相互作用。初步揭示了百合地黄汤抗肝癌作用具有多靶点、多通路的治疗特点，为后续深入验证百合地黄汤抗肝癌的作用机制提供了方向与参考。