摘要 目的： 通过网络药理学及分子对接技术， 探讨肉苁蓉苯乙醇总苷(Phenylethanoid Glycosides from Cistanche, CPhGs)治疗肝细胞癌(HepatocellularCarcinoma, HCC)的潜在分子机制。方法：应用TCMSP等数据库及献收集肉苁蓉苯乙醇总苷的活性成分及潜在作用靶点。通过TTD，GeneCards等数据库获得与HCC相关的疾病基因，得到药物靶点与疾病靶点的交集为肉苁蓉苯乙醇总苷治疗HCC的潜在作用靶点，将交集靶点导入PPI网络，筛选关键靶点。将关键靶点导入DAVID数据库，进行GO功能、KEGG通路富集分析。最终，利用Autodock vina软件对筛选出的蛋白与肉苁蓉苯乙醇总苷活性成分进行分子对接。结果：肉苁蓉苯乙醇总苷中有29种活性成分及410个靶基因，涉及PI3K-Akt信号通路、癌症信号通路、癌症蛋白聚糖通路、VEGF信号通路等以发挥治疗HCC的作用。分子对接结果显示，筛选的靶点受体蛋白与活性成分具有较好的结合活性。结论：肉苁蓉苯乙醇总苷可能通过多组分、多靶点和多途径的方式调节癌细胞、血管生成等相关靶标和途径，从而对HCC产生治疗作用。

关键词 肉苁蓉苯乙醇总苷；肝细胞癌；网络药理学；分子对接

原发性肝癌是全球高发病率和致死率的恶性肿瘤之一，五年内的存活率仅为7%[1]。肝细胞癌（Hepatocellular Carcinoma，HCC）是其中最主要的病理类型，约占85%~90%以上[2]。最新数据表明，我国HCC的病死率和发病率分别排在所有恶性肿瘤的第3位和第4位，具有发展较快、易复发的特点[3-4]。肝癌的病因和发病机制尚不清楚，其发生可能与宿主基因多态性、环境因素(包括代谢综合征和黄曲霉毒素B1)、饮酒和病毒因素(包括乙型肝炎病毒和丙型肝炎病毒感染)的相互作用有关[5]。中药可用于晚期HCC患者，具有抗肿瘤、减少放疗化疗的不良反应、提高肿瘤患者免疫力、改善生活质量等作用[6]。

肉苁蓉为列当科植物肉苁蓉Cistanche deserticola Y.C.Ma或德赢vwin登录 Cistanche tubulosa（Schenk）Wight的干燥带鳞叶的肉质茎[7]，具有补肾、益精、润燥、通便的功效，是我国西北干旱地区特有的名贵补益类药材，应用历史悠久。肉苁蓉始载于东汉《神农本草经》：“味甘、微温。主五劳七伤，补中，除茎中寒热痛，养五脏，强阴，益精气，多子，妇人症瘕，久服轻身”，被列为上品[8]。肉苁蓉苯乙醇总苷(PhenylethanoidGlycosides from Cistanche, CPhGs)作为肉苁蓉中的主要有效成分，包括毛蕊花糖苷、松果菊苷、红景天苷、2’-乙酰基毛蕊花糖苷、管花苷B等[9]，具有肝脏保护、抗痴呆、抗疲劳、缺血保护、抗骨质疏松、润肠通便、免疫调节、抗衰老及改善生殖等作用[10]。已有研究证明，CPhGs具有良好的防治肝癌作用。Yuan等人[11]的研究表明CPhGs通过外源性和内源性凋亡途径抑制H22细胞的生长。胡琼等人[12]发现CPhGs能够抑制肝癌细胞的生长，这可能与CPhGs降低荷瘤小鼠血清中AFP的含量，提高荷瘤小鼠的免疫能力有关。目前，尚未发现针对CPhGs防治HCC作用及其分子作用机制的深入研究，其活性成分、作用靶点及分子机制尚未完全阐明。网络药理学将现有的“一靶一药”方式更新为“网络靶点-多成分”的研究方式，即在此基础上对中药的疗效进行综合分析，从系统视角和分子水平探索中医药，构建中医药网络药理学新思路[13-14]。本研究通过网络药理学策略，基于中医药整体观，对CPhGs防治HCC的作用机制进行系统探讨，结合“中药-活性成分-靶点-信号通路”的多层次网络模式，预测CPhGs防治HCC作用的关键靶点及作用机制，并利用分子对接技术进行验证，阐述CPhGs多成分、多靶点、多通路的协同作用机制，为显著提升CPhGs的药用资源开发领域，为后续研究提供参考。

1 资料与方法

1.1 CPhGs 类成分的收集

通过检索中药系统药理学数据库(TCMSP)、中医药百科全书数据库(ETCM）、中医药研究综合数据库(TCMID)等数据库并结合CNKI 数据库、PubMed 数据库中的相关文献，收集CPhGs 类的化学成分及靶点信息。

1.2 药物潜在治疗靶点的预测及“中药-活性成分-靶标”网络图的构建将检索得到成分的SDF 格式文件导入PharmMapper 数据库，进行潜在靶点的预测分析，并采用UniProt KB 搜索功能对靶点名称进行校正。运用Cytoscape3.8.0 软件构建“中药-活性成分-靶点”网络图。

1.3 HCC 相关基因的筛选

将“Hepatocellular Carcinoma” “HCC”作为关键词在以下4 个数据库中进行检索，获得HCC 靶点的信息：(1)OMIM 数据库；(2)DrugBank 数据库；(3)GeneCards数据库；(4)TTD 数据库。筛选后将靶点信息进行汇总整理。

1.4 CPhGs 治疗HCC 的靶点映射

通过在线绘图软件制作韦恩图，将CPhGs 类成分作用靶点与HCC 靶标进行映射，交集靶点即为CPhGs 类成分治疗HCC 潜在作用靶点，通过韦恩图确定CPhGs 类成分与HCC 的药物-疾病交集靶点。

1.5 PPI 网络的构建

利用STRING数据库（https://string-db.org/）进行共有靶点的蛋白质间相互作用关系分析，设置Combined Score > 0.4，并将结果利用Cytoscape 3.8.0进行拓扑分析，以度值大于中位数的2倍筛选出关键共有靶点。1.6 GO 富集分析和KEGG 通路富集分析利用DAVID 数据库对CPhGs 和HCC 共同作用靶点蛋白进行GO 分析和KEGG 通路富集分析，设定阈值P ≤ 0.05，筛选出具有显著性差异的生物过程和通路，并通过Cytoscape 3.8.0 软件构建“中药-活性成分-靶点-通路”网络。

1.7 核心靶点的分子对接验证

选择关键靶点蛋白用于分子对接，从RCSBPDB 数据库获取这些靶点蛋白的晶体结构文件，从PubChem 数据库获取这些靶点的活性成分的结构文件，利用Autodock_vina 软件计算靶点蛋白与活性成分间的最低结合能并将分子对接模拟验证结果进行可视化。

2 结果

2.1 肉苁蓉苯乙醇总苷的活性成分及其预测靶点通过在TCMSP、ETCM、TCMID数据库及文献[15-16]内检索，最终确定29个化合物作为候选化合物。包括Echinacoside 、Acteoside 、Isoacteoside 、2’-Acetylacteoside、Tubuloside A等，基本信息，见表1。通过数据库及PharmMapper预测出成分的靶标并去重后，获得410个可能的靶标。

2.2 “中药-活性成分-靶标”网络构建与分析

将“2.1”项下的29 个CPhGs 类成分和410 个潜在靶标的相关信息导入Cytoscape 3.8.0 软件中，结果见图1。“中药-活性成分-靶标”网络图涉及439 个节点，7515 条边。图中白色节点为中药；图中黄色节点中的缩写为成分名称，如PhGs9 为Acteoside，其余略；图中蓝色节点为靶标名称，如CASP3 为Caspase-3蛋白。

2.3 HCC 疾病靶点的获取及CPhGs 类成分靶点映射以 “Hepatocellular Carcinoma”“HCC”为关键词检索OMIM、DrugBank、GeneCards、TTD 数据库，通过标准化去重及筛选后共获得HCC 相关靶点1944个。通过韦恩图将CPhGs 成分预测的靶点与HCC 的靶点进行映射，发现共有173 个交集靶点，即CPhGs 治疗HCC 的作用靶点，见图2。

2.4 PPI 网络构建及分析

CPhGs 治疗HCC 的靶点蛋白PPI 网络，见图3A，该网络涉及173 个节点，2202 条边。度值中位数为20，根据度值大于中位数2 倍筛选出核心网络蛋白29个，见图3B，包括ALB、AKT1、EGFR、MAPK1、IGF1、SRC、CASP3、MMP9、HRAS、HSP90AA1、MAPK8、ESR1、MMP2、RHOA、MAPK14、ANXA5、CAT、MDM2、KDR、PIK3R1、CDC42、MAP2K1、STAT1、GRB2、AR、IL2、NOS3、IGF1R、JAK2。以上靶点在PPI 网络中的位置，表明其在CPhGs 治疗HCC 中可能发挥重要作用。

2.5 GO 分类富集分析和KEGG 通路富集分析

将筛选出的 29 个关键靶点进行GO 和KEGG 富集分析，得到生物过程（Biological Process, BP）65 条，主要包括蛋白质磷酸化，细胞增殖、迁移、凋亡，肾上腺素信号通路等；细胞组成（Cellular Components, CC）18 条，主要与细胞质，细胞间连接，细胞内外基质等有关；分子功能（Molecular Function, MF）23 条，主要与蛋白激酶活性、蛋白质结合、酶结合等有关，见图4，以上数据P

值均小于0.05。通过 KEGG 富集得到92 条信号通路（P < 0.05），主要涵盖癌症信号通路（hsa05200: Pathways in cancer）、癌症蛋白聚糖通路（hsa05205: Proteoglycans incancer）、PI3K-Akt 信号通路（hsa04151: PI3K-Akt signaling pathway）、RAS 信号通路(hsa04014: Ras signaling pathway)、FoxO 信号通路（hsa04068: FoxO signalingpathway）、血管内皮生长因子信号通路（hsa04370: VEGF signaling pathway）等，排名前20 的通路，见图5。

2.6 “中药－活性成分－靶点－信号通路”网络构建通过Cytoscape 3.8.0软件构建CPhGs治疗HCC的“中药-活性成分-靶点-通路”网络，见图6，共有79个节点和901条边，涉及29种活性成分、29个关键靶点、20条KEGG通路，体现了CPhGs可能通过多成分、多靶点、多途径发挥多效应的特点。

2.7 分子对接验证

根据“中药-活性成分--靶点”网络中靶点的度值选择CPhGs活性成分5个（毛蕊花糖苷、松果菊苷、2-乙酰毛蕊花糖苷、肉苁蓉苷A、管花苷A）与PPI分析中度值排名前10的靶点蛋白（ALB、AKT1、EGFR、MAPK1、IGF1、SRC、CASP3、MMP9、HRAS、HSP90AA1）进行分子对接，筛选对接结果中结合能低的分子构象即为最优对接结果，其中活性成分与蛋白间的结合能等信息，见表2，具体对接过程，见图7。结果显示，50个对接结果均有较优的（Total score ≤ -5.0）对接活性，其中，43个对接结果有强烈的匹配活性（Total score ≤ -7.0）。以上对接结果充分表明，CPhGs活性组分与HCC关键核心靶基因的结合作用良好。

3 讨论

HCC 是严重威胁人类健康的十大恶性肿瘤之一，具有发病隐匿，进展迅速，治疗困难、转移概率高的特点[17]。HCC 是各种信号通路及细胞因子共同参与的病理生理过程，因此单靶点、单通路治疗很难起效。中药具有多成分、多作用靶点和作用途径复杂等特点，并以“人”为核心，将整体与局部相关联，通过中药有效成分在体内的代谢过程，进而从系统上实现药物对机体的调控。通过“中药-活性成分-靶点”的网络结果显示，CPhGs类成分毛蕊花糖苷、松果菊苷、红景天等均具有较高的防治肝癌活性，在治疗过程中可能发挥重要作用。毛蕊花糖苷对环境、饮食及二乙基亚硝胺（DEN）诱导的大鼠肝癌具有预防作用，其发挥作用可能是由STAT-3依赖性的氧化应激和凋亡调控介导的[18]。另外，毛蕊花糖苷可通过干扰JNK信号通路的激活和调控自噬蛋白的表达来抑制肝癌细胞的增殖[19]。研究表明，松果菊苷（ECH）能抑制二乙基亚硝胺（DEN）诱导的小鼠肝癌，并对HepG2肝癌细胞株具有抗增殖和促凋亡作用，其作用机制可能是通过降低TREM2的表达和PI3K/AKT信号通路发挥其抗肿瘤活性[20]。红景天苷联合奥沙利铂给药可以提高HCC患者的化疗疗效，通过抑制HIF-1α信号通路，显著提高HCC患者对铂类药物的敏感性，抑制缺氧诱导的EMT[21]。综上所述，目前针对CPhGs类成分防治HCC的作用机制研究较少，尚未见有关CPhGs防治HCC作用机制的报道。基于此，本课题应用网络药理学及分子对接技术对CPhGs防治HCC的作用机制进行系统研究，丰富和完善中药、天然药物治疗HCC的研究资料，拓展新疆地产药材肉苁蓉的开发应用价值。

本研究的PPI网络拓扑学分析显示，CPhGs治疗HCC的核心靶点有29个，其中包括AKT1，PIK3R1，STAT1，SRC，MAPK1，CASP3，MAPK8，MAPK14，IGF1R等，均涉及细胞凋亡和细胞增殖过程。采用DAVID数据库对上述29个关键靶点进行GO功能和KEGG富集分析。GO功能富集分析结果显示CPhGs可以对细胞信号转导、细胞迁移、细胞增殖、细胞免疫等生物过程起调控作用，涉及细胞膜、核、质、细胞器膜及细胞空间等细胞组成过程，并通过与蛋白、酶各类受体结合等分子功能产生活性。KEGG功能富集分析结果显示，CPhGs防治HCC可能是通过调控癌症信号通路、癌症蛋白聚糖通路、磷脂酰肌醇3-激酶（PI3K-Akt）信号通路、GnRH信号通路等多条信号转导通路，其中均涉及PI3K/Akt信号通路，其作用尤为关键。

PI3K/Akt信号通路在肝癌的起始过程中起着至关重要的作用，涉及多种生物学过程，如增殖、转移、耐药、放射抗性、能量代谢和自噬等[22-23]。PI3K是由调节亚单位p85和催化亚单位p110组成的异源性二聚体。根据底物和序列同源性，将PI3K分为3个亚型，都能磷酸化磷脂酰肌醇第3位羟基，从而产生特定的磷酸肌醇形式而发挥作用[24]。蛋白激酶B是位于PI3K下游的靶蛋白，是细胞增殖和存活的重要调节因子。目前，已经鉴定出3种Akt基因，分别称为Akt1、Akt2和Akt3。抑制Akt活性可诱导一系列哺乳动物的细胞凋亡，因此，将Akt抑制剂与抗癌药物联合治疗具有突变的肿瘤有一定的疗效[25]。

本研究将“中药-活性成分--靶点”网络中排名前5的核心活性化合物与PPI分析中度值排名前10的蛋白靶点进行分子对接，进一步明确其发挥效用的治疗机制。采用Autodock vina进行对接，选取亲和力最佳的构象，其中最低结合能 ≤ -5.0kcal/mol即表明药效分子与蛋白对接效果良好,结合能越低，表明分子与蛋白结合能力越佳。本研究结果提示肉苁蓉中的苯乙醇苷类活性化合物与AKT1、HSP90AA1等靶标结合后最低结合能均 ≤ -5.0 kcal/mol，表明其具有较强的结合能力，充分验证了网络药理学的预测结果。蛋白激酶B（AKT）是一种人类丝氨酸-苏氨酸激酶，是AGC蛋白激酶家族的成员[26]。AKT蛋白在大多数癌症类型中都高度表达，包括HCC、肺癌、乳腺癌和以及结肠癌，其中AKT1是治疗癌症的主要药物靶点，在细胞增殖和抗凋亡的过程中发挥作用。有研究表明沉默AKT1能够抑制肿瘤生长和诱导细胞凋亡，并可有效上调肿瘤抑制基因PTEN的表达、减少Notch1的表达，靶向AKT1可以抑制HCC的发展[27]。热休克蛋白90α(Heatshock protein 90α，Hsp90α)由HSP90AA1基因编码，是分子伴侣Hsp90的诱导亚型，是真核生物中高度保守且必不可少的伴侣蛋白[28]。HSP90已成为治疗癌症中的关键目标，Hsp90与多种促进肿瘤发生的蛋白相互作用，因此，Hsp90是一个重要的致癌因子。研究表明，HCC患者血清HSP90AA1水平明显高于健康对照组[29]，分泌型HSP90AA1可能是HCC良好的生物标志物，被批准用于中国HCC的临床诊断[30]。

综上所述，本文以CPhGs主要活性成分为研究对象，通过网络药理学技术分析其潜在的防治HCC靶点及信号通路，并结合分子对接方法初步验证了CPhGs主要活性成分能与HCC靶点发生相互作用，为进一步深入探讨其作用机制奠定了基础。