Correlation of HBV DNA load with liver pathology,liver function,and serological markers of liver fibrosis in children with chronic hepatitis B

肝癌作为我国常见的一种恶性程度较高的肿瘤，死亡率较高且呈上升趋势，位居全球癌症的第三位^[1]。现临床多采用手术、放化疗等多种方法联合治疗，但治疗效果尚未达到预期。且相较于毒副反应较大的化学治疗，中医药作为肿瘤治疗常用的辅助治疗方式，其减毒增效的特点有利于提高患者的治疗效果及生存质量。

药对是中医临床医师在治疗过程中常用的一种独具特色的药物配伍方式，其组成相较方剂较为简单，但药对皆由“病、证、症”出发，精简配对，亦可发挥综合效用，增强疗效。本文研究基于国医大师刘尚义教授临床辨治肿瘤的数据挖掘分析。药对“川楝子-延胡索”是刘尚义教授临床常用的一组行气、止痛经方药对^[2]，二药既是药对，亦是名方金铃子散的组成。金铃子散出自于金元四大家之一刘完素所著医书《素问·病机气宜保命集》，其中川楝子可疏肝行气、清肝泻火；延胡索可行气活血、祛瘀止痛，两者配伍具有清肝泻热、活血散瘀定痛之效。此外，现代药理学研究^[3-5]表明川楝子、延胡索均具有良好的抗肿瘤、抗炎镇痛作用。但目前研究尚未有完整、明确的药对“川楝子-延胡索”用于治疗肝癌的药效机制，故开展此次研究提供部分理论基础。

1.   方法

1.1   搜集、筛选药对“川楝子-延胡索”的成分及对应靶点

分别以“川楝子”“延胡索”为关键词在中药系统药理学数据库与分析平台(Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform，TCMSP)中检索搜集中药成分，以“口服生物利用度(OB)≥30%，化合物的类药性(DL)≥0.18”为条件筛选药物有效成分，查找有效成分的对应靶点，并通过Uniprot数据库(http://www.uniprot.org/)将所有作用靶点统一为对应基因名称。

1.2   查找、筛选药对“川楝子-延胡索”治疗肝癌的靶点基因

分别以“hepatic carcinoma”“hepatic cancer”“hepatocellular carcinomas”“liver cancer”为关键词在疾病靶点Genecards数据库(https://www.genecards.org/)中检索查找肝癌的相关基因并导出结果。在Excel表格中选用两倍中位数公式并删除重复靶点后筛选出与肝癌关系更为紧密的基因。最后将药对的药物基因、肝癌基因相映射得出交集基因，该交集基因即为药对治疗肝癌的靶点基因。

1.3   构建“药对-成分-靶点”网络

将药对中“川楝子”“延胡索”筛选后的药物活性成分及靶点基因导入Cytoscape 3.7.1软件中，从而构建药对药物-成分-靶点网络图。

1.4   构建靶点蛋白相互作用网络

将交集靶点基因上传至String数据库(https://string-db.org/)，限定物种为“Homo sapiens”，分析药对“川楝子-延胡索”治疗肝癌的潜在靶点信息及靶点蛋白之间的相互作用关系，并借助Cytoscape 3.7.1软件将所得靶点蛋白信息进行可视化分析。

1.5   分析靶点通路

为了更进一步的分析、明确药对“川楝子-延胡索”的靶点基因在治疗肝癌过程中发挥的生物学功能及相关重要通路，将靶点基因输入至Metascape数据库(https://metascape.org/gp/index.html#/main/step1)进行GO富集分析(含生物过程、细胞组分、分子功能三部分)、KEGG信号通路富集分析，并通过易汉博生物信息在线做图网站Image GP(http://www.ehbio.com/ImageGP/)将分析结果进行可视化处理。

2.   结果

2.1   药对“川楝子-延胡索”的成分及对应作用靶点

在TCMSP数据库中以前述方法进行检索、筛选，得出川楝子有效活性成分9个，其中melianone(MOL002047)、Nimbolidin D(MOL002048)、Nimbolin A(MOL002053)无对应靶点，予以剔除，故其有效成分为6个，对应靶点为181个；延胡索有效活性成分49个，对应靶点1097个；两味中药共同成分1个，为quercetin(MOL000098)，共同靶点基因154个。成分具体情况见表 1。

表  1  药对“川楝子-延胡索”中中药活性成分详情

中药	简称	MOL ID	分子名称	OB(%)	DL(%)
川楝子(简称：CLZ)	CLZ1	MOL001494	Mandenol	42.00	0.19
	CLZ2	MOL001495	Ethyl linolenate	46.10	0.20
	CLZ3	MOL002045	Stigmasterol	43.41	0.76
	CLZ4	MOL002056	(E)-3-[(2S, 3R)-2-(4-hydroxy-3-methoxy-phenyl)-7-methoxy-3-methylol-2, 3-dihydrobenzofuran-5-yl]acrolein	54.74	0.40
	CLZ5	MOL002058	Medioresil	57.20	0.62
延胡索(简称：YHS)	YHS1	MOL001454	Berberine	36.86	0.78
	YHS2	MOL001458	Coptisine	30.67	0.86
	YHS3	MOL001460	Cryptopin	78.74	0.72
	YHS4	MOL001461	Dihydrochelerythrine	32.73	0.81
	YHS5	MOL001463	Dihydrosanguinarine	59.31	0.86
	YHS6	MOL001474	Sanguinarine	37.81	0.86
	YHS7	MOL000217	(S)-Scoulerine	32.28	0.54
	YHS8	MOL002670	Cavidine	35.64	0.81
	YHS9	MOL002903	(R)-Canadine	55.37	0.77
	YHS10	MOL000359	Sitosterol	36.91	0.75
	YHS11	MOL004071	Hyndarin	73.94	0.64
	YHS12	MOL004190	(-)-alpha-N-methylcanadine	45.06	0.80
	YHS13	MOL004191	Capaurine	62.91	0.69
	YHS14	MOL004193	Clarkeanidine	86.65	0.54
	YHS15	MOL004195	Corydaline	65.84	0.68
	YHS16	MOL004196	Corydalmine	52.50	0.59
	YHS17	MOL004197	Corydine	37.16	0.55
	YHS18	MOL004198	Corynoline	46.06	0.85
	YHS19	MOL004199	Corynoloxine	38.12	0.60
	YHS20	MOL004200	Methyl-[2-(3, 4, 6, 7-tetramethoxy-1-phenanthryl)ethyl]amine	61.15	0.44
	YHS21	MOL004202	Dehydrocavidine	38.99	0.81
	YHS22	MOL004203	Dehydrocorybulbine	46.97	0.63
	YHS23	MOL004204	Dehydrocorydaline	41.98	0.68
	YHS24	MOL004205	Dehydrocorydalmine	43.90	0.59
	YHS25	MOL004208	Demethylcorydalmatine	38.99	0.54
	YHS26	MOL004209	13-methyldehydrocorydalmine	35.94	0.63
	YHS27	MOL004210	(1S, 8′R)-6, 7-dimethoxy-2-methylspiro[3, 4-dihydroisoquinoline-1, 7′-6, 8-dihydrocyclopenta [g][1, 3]benzodioxole]-8′-ol	43.95	0.72
	YHS28	MOL004763	Izoteolin	39.53	0.51
	YHS29	MOL004214	Isocorybulbine	40.18	0.66
	YHS30	MOL004215	leonticine	45.79	0.26
	YHS31	MOL004216	13-methylpalmatrubine	40.97	0.63
	YHS32	MOL004220	N-methyllaurotetanine	41.62	0.56
	YHS33	MOL004221	Norglaucing	30.35	0.56
	YHS34	MOL004224	Pontevedrine	30.28	0.71
	YHS35	MOL004225	Pseudocoptisine	38.97	0.86
	YHS36	MOL004226	Pseudoprotopine	53.75	0.83
	YHS37	MOL004228	Saulatine	42.74	0.79
	YHS38	MOL004230	Stylopine	48.25	0.85
	YHS39	MOL004231	Tetrahydrocorysamine	34.17	0.86
	YHS40	MOL004232	Tetrahydroprotopapaverine	57.28	0.33
	YHS41	MOL004233	Thaliporphine	31.87	0.56
	YHS42	MOL004234	2, 3, 9, 10-tetramethoxy-13-methyl-5, 6- dihydroisoquinolino[2, 1-b]isoquinolin-8-one	76.77	0.73
	YHS43	MOL000449	Stigmasterol	43.83	0.76
	YHS44	MOL000785	Palmatine	64.60	0.65
	YHS45	MOL000787	Fumarine	59.26	0.83
	YHS46	MOL000790	Isocorypalmine	35.77	0.59
	YHS47	MOL000791	Bicuculline	69.67	0.88
	YHS48	MOL000793	Bulbocapnine	47.54	0.69
共同成分	A1	MOL000098	Quercetin	46.43	0.28
注：简称用于Cytocape软件制图。

下载: 导出CSV 

| 显示表格

将药对“川楝子-延胡索”对应靶点基因与Uniprot数据库中对应人类基因进行匹配，筛选后得出川楝子对应靶点164个，延胡索对应靶点961个，共同靶点143个，去重后药对“川楝子-延胡索”对应靶点基因共199个。

2.2   药对“川楝子-延胡索”治疗肝癌的靶点基因

检索肝癌的相关基因共42 992个，经两倍中位数及去重筛选后得出肝癌的相关基因共4340个；其中与药对“川楝子-延胡索”对应靶点基因相交集的基因共162个。

2.3   “药对-成分-靶点”网络

对川楝子、延胡索筛选药物活性成分，交集靶点基因导入至Cytoscape 3.7.1软件中构建“药对-成分-靶点”网络图(附录1)。该网络图中共含有254个节点(包括2味中药、54个化合物、162个潜在靶点基因)，1180条边。分析该网络中相关节点的度值可知，中药活性成分中度值最大的是槲皮素(quercetin，A1)，其次为(S)-金黄紫堇碱[(S)-Scoulerine，YHS7]、异紫堇杷明碱(Isocorypalmine，YHS46)、刺罂粟碱(stylopine、YHS30)、延胡索乙素(Hyndarin、YHS11)，说明该药对在治疗肝癌的过程中发挥主要作用的可能是以上活性成分。

2.4   药对“川楝子-延胡索”治疗肝癌的靶点蛋白相互作用网络

将162个交集靶点基因上传至String数据库以分析蛋白之间的相互作用并构建PPI网络图。此外将该分析数据以TSV格式导出并保存。通过Cytoscape 3.7.1软件将该分析数据进行可视化分析，得出药对“川楝子-延胡索”治疗肝癌的靶点蛋白互作网络图(附录2)。分析数据可知，AKT1靶点蛋白中度值最高(118)，其次为IL6(103)、TP53(100)、VEGFA(95)、TNF(90)、CASP3(88)、JUN(88)、MAPK1(88)、MYC(86)、EGFR(83)，说明药对在治疗肝癌时以上基因发挥了重要作用。

2.5   药对“川楝子-延胡索”治疗肝癌的靶点基因的GO功能、KEGG信号通路分析

借助Metascape数据库分析162个交集靶点基因，筛选后得出GO富集条目共2606个(P值均 < 0.01)，包括2321个生物学进程(BP)，主要涉及对有毒物质、无机物、创伤愈合的反应；细胞对有机环状化合物、有机氮化合物及外界刺激的反应；细胞增殖的负调控、细胞死亡的正调控；血管生长发育、凋亡的信号通路等；186个分子功能(MF)，主要涉及蛋白激酶结合、蛋白酶结合、磷酸酶结合、氧化还原酶活性、转录因子结合、细胞因子受体结合、蛋白质结构域特异性结合、蛋白质同聚活性、核受体活性、激素受体结合等；99类细胞组成(CC)，主要涉及膜筏、质膜蛋白复合物、神经细胞体、转录因子复合物、细胞质核周区、细胞外基质、细胞顶部、黏着连接、溶质液泡、囊泡腔等(图 1)。此外，KEGG分析得出相关信号通路共351条，其中主要涉及癌症通路(Pathways in cancer)、AGE-RAGE信号通路(AGE-RAGE signaling pathway in diabetic complications)、流体剪切应力与动脉粥样硬化(Fluid shear stress and atherosclerosis)、TNF信号通路(TNF signaling pathway)、NF-κB信号通路(NF-kappa B signaling pathway)、结直肠癌(Colorectal cancer)、甲状腺激素信号通路等(图 2)。

图  1  药对“川楝子-延胡索”治疗肝癌的靶点基因的GO功能分析(含BP、MF、CC)

下载: 全尺寸图片 幻灯片

图  2  药对“川楝子-延胡索”治疗肝癌的靶点基因的KEGG信号通路分析

下载: 全尺寸图片 幻灯片

3.   讨论

随着网络药理学的发展，许多研究人员采用网络药理学方法分析中药成分、评估药理作用，并通过可视化分析药物分子间的网络关系及相关靶蛋白、通路数据，以探讨中药、方剂发挥疗效的可能分子机制^[6-8]。如巩晓杰等^[9]通过网络药理学研究薏苡仁治疗肝癌的有效活性成分；欧阳思露等^[10]通过网络药理学分析莪术的抗癌机理；李菁等^[11]应用网络药理学、分子对接探析四逆散治疗肝癌的潜在分子机制。因此，网络药理学对于中药的药效开发以及预测潜在的作用机制具有重要意义。

药对“川楝子-延胡索”现临床广泛应用于慢性肝炎、肝硬化、胆囊炎等肝胆疾病，此外亦用于癌症的治疗^[12]。研究^[13-14]表明，通过半仿生提取技术处理金铃子散并将其作用于人宫颈癌HeLa细胞、人肝癌HepG2细胞后，均表达了一定的抗癌作用，且对肝癌细胞作用更为突出；川楝子粗提取物可通过调节人结肠癌细胞SW480、小鼠大肠腺癌细胞CT26的凋亡进程以抑制癌细胞的增殖^[15]；延胡索粉末对小鼠H22肝癌具有明显抑制作用且与其浓度呈一定依赖性^[16]。

此次研究通过网络药理学方法分析药对“川楝子-延胡索”中中药活性成分，构建成分、肝癌靶点相关网络，说明药对治疗肝癌是多个成分、多个靶点共同作用的综合过程。

通过上述分析发现药对中发挥重要作用的成分可能是槲皮素、(S)-金黄紫堇碱、异紫堇杷明碱、刺罂粟碱、延胡索乙素。槲皮素是多种中药的共有成分，也是肿瘤研究的热点成分。槲皮素属于黄酮类药物，其可通过调控细胞生长周期、降低细胞代谢活性、诱导细胞凋亡以抑制肝癌细胞增殖，从而发挥抗肝癌作用^[17-20]。(S)-金黄紫堇碱、异紫堇杷明碱、刺罂粟碱、延胡索乙素均属于异喹啉类生物碱，此类生物碱为延胡索中的主要活性成分，药理研究^[21-23]表明其具有抗肿瘤、镇静、抗炎等多种作用。在抗肿瘤作用方面，异喹啉类生物碱可发挥明显的抗肝癌效用，其主要机制在于通过不同通路促进肿瘤细胞发生凋亡与自噬，抑制细胞的增殖与侵袭转移，控制肿瘤血管生成^[24-25]。

通过药对“川楝子-延胡索”治疗肝癌的靶点蛋白相互作用网络可知该网络关系中的主要作用包括AKT1、IL6、TP53、VEGFA、TNF、CASP3等靶点。分析以上主要发挥治疗肝癌作用的靶点及作用，可将其总结为以下几个方面。(1)炎症相关因子介导肿瘤进程(IL6、TNF)：IL6是具有多种效应的的炎症细胞因子，当机体出现急性损伤时呈现高表达状态^[26]，且参与和调控了肺癌、肝癌等多种恶性肿瘤的发生、发展^[27-28]。TNF包括TNFα和TNFβ，是一类能直接杀伤并造成肿瘤细胞死亡的细胞因子。其中TNFα在调控肿瘤的发展进程中具有双向调节作用，既具有抗肿瘤作用亦可促进肿瘤生长，而具体表现为何种作用取决于TNFα的浓度。当TNFα浓度较高时，可刺激其他相关炎症因子的释放、抑制肿瘤血管新生以最大程度发挥抗肿瘤效应；当TNFα浓度较低时可通过NF-κB信号通路诱导肿瘤细胞侵袭、转移发生以发挥促肿瘤作用^[29-30]。(2)抑制血管新生：VEGFA属于血管内皮生长因子(VEGF)，可升高血管内皮细胞通透性，影响血管内皮细胞增殖、存活、迁移以促进血管生成^[31]。肝癌血供丰富，新生血管的生成可直接影响其病程发展与预后，因此降低VEGF的表达、抑制血管生成、改善肝癌微环境可有效控制肝癌的生长、转移^[32]。(3)抑制肿瘤细胞增殖、促进肿瘤细胞凋亡(AKT1、TP53、CASP3)：AKT1(即AKT)在PI3K/Akt信号通路中发挥重要作用，是通路中的关键因子。研究^[33]表明，在多种肿瘤中均可发现PI3K/Akt信号通路对于促进肿瘤细胞的增殖、生长具有调控作用，这与AKT1的磷酸化被异常活化有紧密关系。此外，研究^[34]表明该信号通路的活化与肝癌的生长有关，且与肝癌患者临床病理特征密切相关。TP53(即P53)是常见的抑癌基因，可抑制肿瘤细胞生长、诱导细胞凋亡^[35]。研究^[36]发现P53容易发生突变，人类约50%的癌症均存在P53突变。突变后P53丧失抑癌活性，在肿瘤细胞增殖中起促进作用，且抑制细胞凋亡，从而加速肿瘤发展^[37]。CASP3(即Caspase 3)是Caspase家族中重要的一员，该家族在细胞凋亡的执行过程中具有不可或缺的作用^[38-39]，且Caspase 3的激活也标志着细胞的凋亡必然开始^[40-41]。而肿瘤的发生发展与肿瘤细胞的凋亡密切相关，研究^[42]表明可通过诱导细胞凋亡从而抑制肝癌细胞的增殖。

GO分子功能阐明药对发挥作用主要与细胞的增殖与凋亡、血管的生长发育有关；KEGG通路分析结果表明药对治疗肝癌过程涉及多个通路，与之相关的通路包括癌症通路、AGE-RAGE信号通路、TNF信号通路、NF-κB信号通路等。多项研究^[43-48]表明NF-κB信号通路参与调控肿瘤的增殖、凋亡过程。NF-κB是细胞中的重要转录调节因子，在大多数细胞中，NF-κB表现为静息状态^[49]。在炎症等刺激下，上游信号通路活化后可激活NF-κB，而激活后的NF-κB可与对应的靶序列在细胞核内结合并调控基因转录、调节细胞黏附分子等以发挥生物学效应^[50-51]。AGE-RAGE信号通路、TNF信号通路均为其上游通路。在AGE-RAGE信号通路中，AGE与其受体RAGE结合后可活化NF-κB信号通路，使NF-κB呈现高表达状态，进而调控其下游基因IL6、TNFα的表达^[52]。亦有研究^[53-55]表明TNF与TNF受体1结合激活TNF信号通路后，可刺激细胞生长因子的表达以促进肝癌细胞的侵袭，诱导VEGF的表达上调以促进肝癌的血管新生从而加速肝癌的恶化，而此过程NF-κB信号通路具有中心主导作用。通路结果与重要靶点分析相吻合，互相验证了文中靶点、通路治疗肝癌的可能机制。

上述分析说明了药对“川楝子-延胡索”治疗肝癌是多种成分、多个靶点基因与信号通路之间相互影响、共同发挥作用的过程，其治疗过程主要是靶点通过介导相关通路调控肿瘤的增殖、凋亡。