1. 引言
胶质瘤起源于神经胶质祖细胞,是中枢神经系统最常见的肿瘤,以弥漫性组织浸润伴周围脑实质破坏为主要特征、手术切除辅以放化疗为主要治疗手段。现代医学模式下中医疗法辨证统一、因人制宜、阴阳平衡的观念优势日益凸显,中西医结合抗肿瘤增强药物敏感性、缓解不良反应、改善精神状态、提高生活质量势在必行。扶正抑瘤汤以红芪、当归、莪术、墓头回3:1:1:3配伍,扶正祛邪、益气养阴,是甘肃省知名中医多年临床应用疗效验证的经验方,课题组前期研究也显示扶正抑瘤汤能协同II型干扰素(Interferon-γ, INF-γ)基因转染的骨髓间充质干细胞(Bone Mesenchymal Stem Cells, BMSCs)抑制胶质瘤裸鼠肿瘤生长,抑瘤的同时减毒增效,维持BMSCs稳态 [1] [2]。为进一步探索扶正抑瘤汤治疗胶质瘤的关键有效成分和关键靶点,挖掘分析生物网络药理学大数据,明确扶正抑瘤汤治疗胶质瘤的作用机制,奠定胶质瘤靶向治疗理论基础。
2. 资料与方法
2.1. 收集资料
2.1.1. 收集扶正抑瘤汤有效成分及靶基因
利用中药系统药理学数据库与分析平台(TCMSP数据库https://old.tcmsp-e.com/tcmsp.php)以药物口服生物利用度(OB)大于30%,成药性(DL)大于0.18筛选收集扶正抑瘤汤有效成分及对应靶蛋白,并在通用蛋白资源数据库(UniProt数据库https://www.uniprot.org/)为靶蛋白匹配人类已验证的唯一靶基因。
2.1.2. 收集胶质瘤靶基因
利用疾病基因关联数据库(DisGeNET数据库https://www.disgenet.org/)以Glioma为关键词、疾病基因分数大于0.3筛选收集胶质瘤靶基因;利用人类基因数据库(GeneCards数据库https://www.genecards.org/)以Glioma和Brain glioma为关键词、疾病基因相关分数大于等于10筛选收集胶质瘤靶基因,合并去重以上两数据,所得基因作为胶质瘤靶基因。
2.1.3. 收集扶正抑瘤汤与胶质瘤交集靶基因
在维恩平台2.1.0 (Venny平台https://bioinfogp.cnb.csic.es/tools/venny/)分别导入步骤2.1.1和2.1.2所得数据获取扶正抑瘤汤与胶质瘤靶基因交集并绘制维恩图。
2.2. 分析数据
2.2.1. 分析获得胶质瘤发病关键靶基因
利用蛋白质相互作用网络功能富集分析数据库11.0 (String数据库https://string-db.org/)分析疾病基因相关分数降序排列前100位的胶质瘤靶基因蛋白互作网络(PPI network),筛选出证据可信、最小互作分数大于0.7、最大互作蛋白数小于50的蛋白质,经网络数据集成、分析和可视化软件Cytoscape3.8.2软件CytoHubba插件计算节点中心度,以中心度值大于中位值降序排列前10位的靶基因作为胶质瘤发病的关键靶基因。
2.2.2. 分析获得扶正抑瘤汤治疗胶质瘤的关键靶基因
将方剂疾病交集基因导入String数据库,筛选证据可信、最小互作分数大于0.7、最大互作蛋白数小于50的蛋白质,利用CytoHubba插件计算节点中心度并输出蛋白互作网络图,以中心度值大于中位值降序排列前11位的基因作为扶正抑瘤汤治疗胶质瘤的关键靶基因。
2.2.3. 分析获得扶正抑瘤汤治疗胶质瘤的关键药物有效成分
利用CytoHubba插件分析扶正抑瘤汤有效成分与中心度值前50顺位的方剂疾病交集基因之间的关系网络并计算节点中心度,中心度值较高者确定为扶正抑瘤汤治疗胶质瘤的关键有效成分。
2.2.4. 扶正抑瘤汤治疗胶质瘤功能富集分析
利用大卫生物信息资源数据库6.8 (DAVID数据库https://david.ncifcrf.gov/)对方剂疾病交集基因进行基因本体分析(GO分析),以P值、FDR值均小于0.05、基因计数降序排列前5位筛选获得扶正抑瘤汤治疗胶质瘤的主要生物功能并使用GraphPad Prism8软件绘图。
利用京都基因和基因组百科全书(KEGG数据库https://www.genome.jp/kegg/)对方剂疾病交集基因作通路分析,找出关键靶基因对应的通路作为扶正抑瘤汤治疗胶质瘤的主要通路。
2.3. 构建扶正抑瘤汤治疗胶质瘤作用机制网络图
深入分析总结前7步获得的数据,将方剂–有效成分、有效成分–靶基因、KEGG通路–靶基因三组对应数据导入Cytoscape,绘制扶正抑瘤汤治疗胶质瘤作用机制网络图。
3. 结果
3.1. 资料收集结果
3.1.1. 获得扶正抑瘤汤有效成分及靶基因
通过TCMSP数据库筛选收集到扶正抑瘤汤有效成分共17个(见表1),分别是红芪13个、当归2个、莪术1个、墓头回1个,其中当归和墓头回的有效成分均与红芪重合,去重后余14个有效成分,对应216个靶基因。
![](Images/Table_Tmp.jpg)
Table 1. Active ingredients list of Fuzheng Yiliu Decoction
表1. 扶正抑瘤汤有效成分表
3.1.2. 获得胶质瘤靶基因
通过DisGeNET数据库筛选到胶质瘤靶基因71个,通过GeneCards数据库以Glioma为关键词筛选到胶质瘤靶基因149个,以Brain Glioma为关键词筛选到靶基因1456个,两个数据库合并去重获得胶质瘤靶基因1456个。
3.1.3. 获得扶正抑瘤汤与胶质瘤交集靶基因
利用Venny平台绘制扶正抑瘤汤与胶质瘤靶基因交集图(见图1),获得141个方剂疾病交集基因。
![](//html.hanspub.org/file/101-1574652x8_hanspub.png?20140216001514100)
Figure 1. Venn diagram of intersection genes of Fuzheng Yiliu Decoction and glioma
图1. 扶正抑瘤汤与胶质瘤交集基因维恩图
3.2. 数据分析结果
3.2.1. 发现胶质瘤发病关键靶基因
Cytoscape绘制胶质瘤发病靶基因蛋白互作网络图(见图2),节点颜色越深中心度值越高,取中心度值大于中位值前10顺位的靶基因作为胶质瘤发病关键靶基因,分别是AKT1、p53、VEGFA、MAPK1、STAT3、MYC、PTEN、CCND1、PIK3CA、EGF。
![](//html.hanspub.org/file/101-1574652x9_hanspub.png?20140216001514100)
Figure 2. Protein-Protein interaction network of target gene in glioma
图2. 胶质瘤发病靶基因蛋白互作网络图
3.2.2. 发现扶正抑瘤汤治疗胶质瘤的关键靶基因
Cytosape绘制方剂疾病交集基因蛋白互作网络图(见图3),节点颜色越深中心度值越高,取中心度值大于中位值前11顺位的基因作为扶正抑瘤汤治疗胶质瘤的关键靶基因(见表2),其中有5个治疗关键靶基因与胶质瘤发病关键基因吻合,标记为绿色,分别为p53、AKT1、MYC、MAPK1,VEGFA。
![](//html.hanspub.org/file/101-1574652x10_hanspub.png?20140216001514100)
Figure 3. Protein-Protein interaction network of intersection gene between Fuzheng Yiliu Decoction and glioma
图3. 扶正抑瘤汤与胶质瘤交集基因蛋白互作网络图
![](Images/Table_Tmp.jpg)
Table 2. Key target genes of Fuzheng Yiliu Decoction in treating glioma
表2. 扶正抑瘤汤治疗胶质瘤关键靶基因
3.2.3. 发现扶正抑瘤汤治疗胶质瘤发挥作用的关键药物有效成分
用CytoHubba分析方剂有效成分与中心度值前50顺位的方剂疾病交集基因之间的网络关系并输出网络分析图(见图4),节点颜色越深,代表方剂组分发挥作用越关键,其中左侧顺位排列的11个节点为扶正抑瘤汤治疗胶质瘤的关键靶基因,绿色是和胶质瘤发病关键靶基因吻合的5个靶基因节点,图中深紫色节点红芪槲皮素(quercetin)与关键靶基因连线最多,中心度值最大,取中心度值降序前3位作为扶正抑瘤汤治疗胶质瘤的关键有效成分,既槲皮素(quercetin)、柚皮素(naringenin)、芒柄花黄素(formononetin),均为红芪所含有效成分。
![](//html.hanspub.org/file/101-1574652x11_hanspub.png?20140216001514100)
Figure 4. Network analysis diagram of the key active components of Fuzheng Yiliu Decoction in treating glioma
图4. 扶正抑瘤汤治疗胶质瘤关键有效成分网络分析图
3.2.4. 扶正抑瘤汤治疗胶质瘤主要功能富集
DAVID数据库方剂疾病交集基因GO分析结果(见图5)显示扶正抑瘤汤主要通过促进RNA聚合酶II启动子转录、DNA转录、信号转导等生物学过程在细胞组分上发挥蛋白结合、酶结合、ATP结合等分子功能治疗胶质瘤。
KEGG数据库分析方剂疾病交集基因获得263条相关通路,筛选与扶正抑瘤汤治疗胶质瘤关键靶基因有关的15条通路列表(见表3),可见扶正抑瘤汤可能通过PI3K/Akt通路、MAPK通路、Ras通路、p53通路、Ca离子通路、VEGF信号通路等多条通路发挥功效。
![](//html.hanspub.org/file/101-1574652x12_hanspub.png?20140216001514100)
Figure 5. Gene Ontology analysis of Fuzheng Yiliu Decoction in treating glioma
图5. 扶正抑瘤汤治疗胶质瘤基因本体分析结果
3.3. 构建扶正抑瘤汤治疗胶质瘤作用机制网络图
综合分析前7项数据,通过Cytoscpae构建4味中药、14个有效成分、141个交集靶基因、15条通路连结而成的扶正抑瘤汤治疗胶质瘤作用机制网络图(见图6),由图可知红芪13条有效成分最多,墓头回有效成分与红芪关键有效成分槲皮素(MOL000098)重合,当归有效成分与红芪有效成分豆甾醇(MOL00049)、植物甾醇(MOL000358)重合,莪术有自己独有的有效成分常春藤素(MOL000296),四种中药组成扶正抑瘤汤,通过多条通路,多个靶点发挥综合治疗胶质瘤功效。
![](//html.hanspub.org/file/101-1574652x13_hanspub.png?20140216001514100)
Figure 6. Network diagram of Fuzheng Yiliu Decoction in treating glioma
图6. 扶正抑瘤汤治疗胶质瘤作用机制网络图
4. 讨论
肿瘤分子靶向疗法较手术治疗、放化疗等副作用更小,疗效更佳。骨髓间充质干细胞(BMSCs)是肿瘤生物治疗的理想载体,在肿瘤组织中具有较好的归巢作用,但其成瘤性、迁移性也是它临床应用安全性的壁垒所在。扶正抑瘤汤协同白细胞介素12 (IL-12)基因修饰的BMSCs具有抗裸鼠胶质瘤作用 [3],且含药血清对胶质瘤微环境下BMSCs增殖具有抑制作用,但不影响其肿瘤定位能力,这些发现为扶正抑瘤汤联合BMSCs靶向治疗胶质瘤提供可能 [4]。
通过此次研究进一步明确扶正抑瘤汤对胶质瘤的作用机制,发现方剂关键有效成分。本研究发现扶正抑瘤汤中红芪有效成分最多,所含槲皮素、柚皮素、芒柄花黄素为方剂治疗胶质瘤的关键有效成分,三者均为红芪所含黄酮类化合物,其中槲皮素几乎对应所有方剂治疗胶质瘤的关键靶基因,位列有效成分之首。实验研究证明槲皮素可通过抑制细胞增殖、诱导细胞凋亡显著抑制大鼠颅内C6胶质瘤细胞生长 [5];桑东平等研究显示槲皮素可以通过抑制热休克蛋白27 (Hsp27)增强胶质母细胞瘤细胞对标准化疗药物替莫唑胺(TMZ)的敏感性 [6];也可以通过槲皮素预处理,提高胶质瘤祖细胞自噬活性,增强瘤细胞对放疗照射的敏感性 [7];槲皮素还具有促进免疫器官功能恢复、提高机体抗氧化水平的功能 [8],这些证据都充分说明槲皮素在治疗胶质瘤中的重要地位。
柚皮素抑瘤机制也有相关研究支持,柚皮苷可以通过抑制磷酸肌醇三激酶(PI3K)的磷酸化抑制胶质瘤细胞生长,大约有70%的胶质瘤发现PI3K的异常激活,而较低浓度的柚皮苷即可抑制PI3K活化水平,并对其下游苏氨酸蛋白激酶(Akt)的磷酸化水平发挥抑制作用 [9]。齐冰丽等 [10] 发现柚皮素通过PI3K/Akt/NF-κB通路可影响下游线粒体途径中抑凋亡蛋白(Bcl-2)和促凋亡蛋白(Bax)的表达,调节卵巢癌细胞增殖与凋亡,这一结果与本研究所得方剂治疗胶质瘤的关键通路和靶基因类似。
应用芒柄花黄素开发抗癌药物也有诸多进展,一些经典肿瘤信号通路均被证明是芒柄花黄素的作用靶点 [11],在乳腺癌裸鼠模型中芒柄花黄素作为成纤维细胞生长因子受体2抑制剂(FGFR2),能阻断FGFR2介导的PI3K/Akt信号通路 [12],降低肿瘤细胞活性;芒柄花黄素具有抗血管生成潜能,抑制血管内皮细胞生长因子活性 [13];芒柄花黄素还可以降低多药耐药(MDR)的发生,显著增强替莫唑胺对胶质瘤C6细胞的化疗疗效 [14]。
当归所含2个有效成分均为植物甾醇,能显著增加荷瘤小鼠巨噬细胞吞噬率,具有一定的免疫调节作用 [15];莪术具有活血化瘀之功效 [16],所含1个有效成分为三萜皂苷类化合物常春藤素;墓头回虽然仅含1个有效成分,但与红芪槲皮素重合,协同增强红芪抑瘤作用。
至此,可以明确红芪在扶正抑瘤汤治疗胶质瘤中的君药地位,当归、墓头回有效成分与红芪重合辅助君药抑瘤应为臣药,莪术调和诸药,针对兼病起扶正功效。四味中药君臣有序、四象平衡共同发挥扶正祛邪,益气养阴之用,这也进一步验证了经验方中红芪、当归、莪术、墓头回3:1:1:3配伍方案的科学性。
本次探索研究还发现扶正抑瘤汤通过PI3K/Akt通路、MAPK通路、p53通路、Ras通路、VEGF通路等多条通路作用于AKT1、MAPK1、p53、VEGFA等多个靶基因发挥治疗作用。证据表明神经纤维蛋白1 (NF1)和肿瘤抑制蛋白p53 (p53)的缺失增强了胶质瘤的形成,NF1可以抑制Ras/MAPK通路活性,降低细胞重现性 [17]。Ras-Raf-MAPK通路在受到细胞外信号刺激后相继激活癌基因Ras、丝/苏氨酸激酶(Raf-1)、MAPK激酶最终激活MAPK,被激活的MAPK转至细胞核内激活转录因子调节细胞增殖、转化和凋亡。辛艳超等 [18] 研究显示p38MAPK在胶质瘤中的异常表达随胶质瘤病理级别升高而逐渐增强,闻公灵等研究也证实,miR-206可靶向结合MAPK1抑制该通路活性调节胶质瘤细胞放射敏感性 [19]。
p53为抑癌基因,编码p53肿瘤蛋白调节细胞周期。孙晓玲等 [20] 对90例原发性星形细胞瘤患者临床病理资料整理分析发现,p53在星形细胞瘤组织中阳性率为58.89%显著高于正常脑组织,且瘤体越大、恶性程度越高阳性率越高,p53阳性患者的生存期也均短于阴性表达患者。p53的突变与染色质重塑调节因子ATRX失活显著相关,ATRX突变激活端粒延长替代机制(ALT),是肿瘤中一种不依赖于端粒酶活性的端粒长度维持机制 [21]。
c-Myc是典型的促癌基因,参与调控细胞周期、细胞分化和凋亡。老年脑胶质瘤组织中c-Myc蛋白的阳性率为59.6%且与病理分化显著相关,较正常脑组织阳性率差异显著 [22]。王宸等研究提示冬凌草可能通过Hippo信号通路促进转录共激活因子相关蛋白(YAP)磷酸化下调进入细胞核的YAP量进而减少下游c-Myc的表达发挥抗瘤作用 [23]。
VEGF信号通路靶向药物肿瘤治疗是近年来研究的热点,脑胶质瘤组织中VEGF表达水平异常升高与肿瘤微血管形成紧密相关 [24],榄香烯联合VEGF多克隆抗体有助于抑制脑胶质瘤体组织生长 [25]。
尽管此次论文结果有诸多理论数据支持,但由于靶基因和相关通路是基于网络药理学预测获得,扶正抑瘤汤在具体应用过程中的实际通路和效应还需进一步利用分子对接、人群差异表达分析、实验研究、临床观察等方法验证。
5. 结论
综上所述,扶正抑瘤汤以红芪为方剂之首,四组分合理配伍;以黄酮类化合物为有效成分之首,多成分协同增效;通过AKT1、MAPK、p53、VEGF、MYC等多条靶基因相关通路综合发挥胶质瘤治疗作用。
基金项目
2022年甘肃中医药大学研究生创新创业基金(2022CX69)。
NOTES
*通讯作者。