Network Pharmacology-based Analysis of the Anti-atherosclerosis Mechanism of Scutellarin and Experimental Validation in Vivo
-
摘要:
目的 基于网络药理学分析灯盏乙素(Scutellarin) 治疗动脉粥样硬化(atherosclerosis,AS) 潜在的分子生物学机制。 方法 通过PubChem数据库进行灯盏乙素潜在的靶点预测,运用Disgenet,CeneCards筛选与AS关联的靶基因,利用Cytoscape 3.9.0建立疾病-靶点-成份可视化网络,使用String构建蛋白互作网络图,利用David数据库进行GO及KECG的网络构建。另外,高脂喂养复制APOE-/-小鼠AS模型, 用不同浓度的灯盏乙素灌胃治疗,利用RT-PCR方法验证筛选所得核心通路的表达变化。 结果 (1)网络药理学预测结果:灯盏乙素治疗AS的核心靶点共有肿瘤坏死因子-α (TNF-α),AKT丝氨酸/苏氨酸激酶(AKT serine/threonine kinase,MAPK),低密度脂蛋白受体(Low density lipoprotein receptor),载脂蛋白E(apolipoprotein E,APOE)等共106个。灯盏乙素调控PI3K/AKT信号通路(hsa04151)、Ras信号通路(hsa04014)、MAPK信号通路(hsa04010)(P < 0.05)等涉及细胞的增殖,炎症反应,氧化应激,及细胞凋亡等生物学过程来治疗AS;(2) 体内实验验证的结果:和模型组相比,灯盏乙素组呈浓度依耐性地下调PI3K、AKT、mTOR和Bcl-2 的mRNA表达,上调Bax和Caspase-3的mRNA表达,( P < 0.05)。另外,和模型组相比,灯盏乙素则呈浓度依耐性地下调小鼠血清TNF-α和IL-1β的含量( P < 0.05)。 结论 经过网络药理学筛选和RT-PCR验证,灯盏乙素通过调控PI3K/AKT/mTOR信号通路抑制APOE-/-小鼠AS模型炎症反应,促进细胞凋亡而治疗AS。 -
关键词:
- 网络药理学 /
- 灯盏乙素 /
- 动脉粥样硬化 /
- APOE-/- 小鼠AS模型 /
- PI3K/AKT/mTOR
Abstract:Objective To analyze the potential molecular biological mechanism of scutellarin for the treatment of atherosclerosis (AS) based on network pharmacology. Methods The PubChem database was used to predict the potential targets of scutellarin, and the Disgenet and CeneCards databases were used to predict and screen the genes related to atherosclerosis, and the Cytoscape 3.7.0 software was used to establish and construct the component-target-disease network visualization network. The GO and KECG pathways were analyzed using the David database. In addition, an APOE-/- mouse AS model was developed and treated with different concentrations of scutellarin, and the expression changes of the core targets were verified by RT-PCR. Results 1. Network pharmacological prediction results: the core targets of Scutellarin for AS treatment were tumor necrosis factor-α (TNF-α), AKT serine/threonine kinase 1 (MAPK), low density lipoprotein receptor (LDLR), apolipoprotein E (APOE), among other targets, with 106 in total. The key targets of Scutellarin are PI3K/AKT signaling pathway (hsa04151), Ras signaling pathway (hsa04014) and MAPK signaling pathway (hsa04010) (P < 0.05), which are involved in cell proliferation, inflammatory response, oxidative stress and apoptosis, to treat AS. 2. Results of in vivo experiments: compared with the model group, the mRNA expression of PI3K, AKT, mTOR and Bcl-2 was down-regulated and the mRNA expression of Bax and Caspase-3 was up-regulated in the Scutellarin group. In addition, compared with model group, Scutellarin down-regulated the contents of TNF-α and IL-1β in serum concentration-dependent (P < 0.05). Conclusion Through network pharmacology screening and RT-PCR verification, Scutellarin inhibits the process of inflammatory response in a mouse model of AS by regulating the PI3K/AKT/mTOR signaling pathway. -
Key words:
- Network pharmacology /
- Scutellarin /
- Atherosclerosis /
- APOE-/- mice AS model /
- PI3K\AKT\mTOR
-
宫颈癌是世界范围内最常见的癌症之一,在所有女性恶性肿瘤中其发病率和死亡率均排名第四[1]。在我国,宫颈癌的发病率和死亡率分别是16.56/10万和11.78/10万,给社会造成了严重的疾病和经济负担[2]。95%以上的宫颈癌是由于持续感染高风险人乳头状瘤病毒(high risk- human papillomavirus,HR-HPV)造成的[3−4]。宫颈癌的发展可分为宫颈上皮内瘤变(cervical intraepithelial neoplasia,CIN)和宫颈癌两个主要阶段[5],在此过程中,HPV感染和宿主遗传因素相互作用,共同决定了疾病的进展[6]。近年来的研究证据表明,宿主的遗传变异,如微小RNA(microRNAs,miRNAs)的多态性在宫颈癌的发生和进展中起到了重要作用[7−10]。
miRNAs是一类长度约为20~24个核苷酸的非编码单链RNA[11−12],通过调控各种细胞信号通路对肿瘤细胞产生重要影响,在包括宫颈癌在内的多种癌症发展中发挥着重要作用,比如转移[13]、免疫浸润[14]、侵袭[15]等。多项研究发现,miR-146a在癌症中的异常表达,既可以发挥肿瘤抑制作用,也可能促进肿瘤形成。例如,miR-146a在乳腺癌中的高表达与肿瘤细胞增殖有关[16];在基底样乳腺癌细胞中的过度表达,导致p53功能改变从而增强了肿瘤发生的潜力[17]。此外,Yue等[18]研究发现,miR-146a基因的单核苷酸多态性(single nucleotide polymorphisms,SNP)会导致其成熟序列表达差异,从而增加个体罹患宫颈癌的风险。
本研究选取miR-146a基因的两个SNP位点rs57095329、rs6864584,并在云南汉族人群中探讨这两个SNP位点与CIN发生风险的相关性,以期为CIN的预防和诊断提供候选靶点。
1. 资料与方法
1.1 研究对象
本次研究共纳入2017年5月至2019年10月间96例在该院被诊断为CIN的患者作为CIN组,另选取同期在该医院225例进行健康体检的个体作为对照组。CIN组的纳入标准为:(1)根据美国国立综合癌症网络(National Comprehensive Cancer Network,NCCN)肿瘤学临床实践指南以及国际妇产科联盟(International Federation of Gynecology and Obstetrics,FIGO)分期指导诊断[19],确诊为癌前病变;(2)根据宫颈病理诊断标准明确CIN病理分期[20];(3)患者取样前未接受放化疗等抗肿瘤辅助手段治疗;(4)患者不具有其他恶性肿瘤以及心血管疾病、糖尿病等慢性疾病;(5)临床资料完整。排除采样前接受过放化疗等抗肿瘤辅助治疗方法的患者以及其他恶性肿瘤患者。本研究方案已通过昆明医科大学第三附属医院伦理委员会批准(KYCS2021193),所有研究对象均签署知情同意书。
1.2 DNA提取
使用QIAamp DNA Blood Mini Kit血液DNA提取试剂盒从全血样本中提取DNA,检测DNA的浓度和纯度,并将DNA样本保存于-80 ℃冰箱备用。
1.3 基因分型实验
采用TaqMan探针法对SNP位点进行基因分型,所使用的TaqMan探针(TaqMan SNP Genotyping Assays)及分型试剂(TaqMan Genotyping Master Mix)均购于美国ABI公司。其中,rs57095329探针的编号为C__90078480_10,rs6864584由公司进行定制合成。具体反应体系和条件请参照笔者之前的研究[12,21−22]。每块384孔反应板中包含3个阳性对照和至少1个阴性对照。
1.4 统计学分析
采用SPSS 26.0软件进行数据统计分析。CIN组和对照组间的年龄为计量资料,经检验符合正态分布且方差齐,以均数±标准差($ \bar x \pm s $)表示,分布差异通过独立样本t检验完成。使用SHEsis在线软件(http://analysis.bio-x.cn/)进行哈迪-温伯格平衡检验(Hardy-Weinberg Equilibrium,HWE)以评估样本的群体代表性,两组间不同SNP位点基因型和等位基因分布频率为计数资料,以例数和百分比[n (%)]表示,并通过卡方检验分析两组间的分布差异[23]。此外,借助SNPstats(https://snpstats.net/start.htm)在线软件对两个SNP位点进行遗传模式分析。检验水准设定为α = 0.05,对于多重比较,P值使用Bonferroni法校正为0.025。
2. 结果
2.1 受试者基础临床信息
本研究共纳入321例受试者,对照组共225例,平均年龄(45.55±12.88)岁;CIN组共96例,平均年龄(45.34±10.65)岁,其中CIN I 期23例(24.0%),CIN Ⅱ期13例(13.5%),CIN Ⅲ期60例(62.5%)。采用独立样本t检验比较对照组和CIN组之间年龄分布,差异无统计学意义(t = 0.812,P = 0.445)。
2.2 SNP位点与CIN和CC的相关性分析
2.2.1 Hardy-Weinberg平衡检验
对两个SNP位点在各样本组间的分布进行HWE平衡检验,结果表明rs57095329、rs6864584均处于HWE平衡(P > 0.05),说明样本具有群体代表性。
2.2.2 SNP位点在对照组和CC、CIN中的分布频率
rs57095329位点A等位基因在CIN组的频率显著低于对照组(P < 0.001,OR = 0.48,95%CI:0.32~0.70),AA基因型在CIN组中的频率显著低于对照组(P < 0.001)。rs6864584位点的等位基因和基因型分布频率在对照组与CIN组之间差异无统计学意义(P = 0.473;P = 0.602),见表1。
表 1 SNP位点的等位基因和基因型在对照组和CIN组中的分布特征[n (%)]Table 1. Distribution characteristics of allele and genotype at SNP loci in the control group and CIN group [n (%)]SNPs 等位基因/基因型 对照组 CIN组 HWE 对照组 vs CIN组 χ2 P χ2 P OR (95%CI) rs57095329 A 370(82.2) 132(68.8) 0.742 0.389 14.325 < 0.001* 0.48(0.32~0.70) G 80(17.8) 60(31.2) A/A 154(68.4) 43(44.8) 15.888 < 0.001* A/G 62(27.6) 46(47.9) G/G 9(4.0) 7(7.3) rs6864584 C 38(8.4) 13(6.8) 0.116 0.733 0.515 0.473 0.79(0.41~1.51) T 412(91.6) 179(93.2) C/C 2(0.9) 0(0.0) 1.014 0.602 C/T 34(15.1) 13(13.5) T/T 189(84.0) 83(86.5) *P < 0.025。 2.2.3 SNP位点在对照组和CIN组的遗传模式分析
在对照组与CIN组的比较中,显性模式(AIC = 380,BIC = 387.6)是rs57095329位点的最优遗传模式,在该模式下,基因型A/G-G/G相对于A/A增加了CIN的发生风险(P < 0.001;OR = 2.67,95%CI:1.64~4.37)。rs6864584位点的最优遗传模式为隐性模式(AIC = 394.2,BIC = 401.8),该位点基因型与CIN发生的风险无显著相关性(P > 0.05),见表2。
表 2 SNP位点在对照组和CIN中的遗传模式分析[n (%)]Table 2. Genetic model analysis of SNP loci in the control and CIN groups [n (%)]SNPs 模型 基因型 对照组 CIN组 对照组 vs CIN组 χ2 P OR (95%CI) rs57095329 共显性 A/A 154(68.4) 43 (44.8) 14.555 < 0.001* 1.00 A/G 62 (27.6) 46 (47.9) 2.66(1.60~4.42) G/G 9 (4.0) 7 (7.3) 2.79(0.98~7.91) 显性 A/A 154(68.4) 43 (44.8) 15.879 < 0.001* 1.00 A/G-G/G 71 (31.6) 53 (55.2) 2.67(1.64~4.37) 隐性 A/A-A/G 216(96.0) 89 (92.7) 1.539 0.230 1.00 G/G 9 (4.0) 7 (7.3) 1.89(0.68~5.23) 超显性 A/A-G/G 163(72.4) 50 (52.1) 12.496 < 0.001* 1.00 A/G 62 (27.6) 46 (47.9) 2.42(1.47~3.97) 逻辑累加 --- --- --- 17.055 < 0.001* 2.12(1.42~3.16) rs6864584 共显性 T/T 189(84.0) 83 (86.5) 0.155 0.450 1.00 T/C 34 (15.1) 13 (13.5) 0.87(0.44~1.73) C/C 2 (0.9) 0 (0.0) --- 显性 T/T 189(84.0) 83 (86.5) 0.314 0.570 1.00 T/C-C/C 36 (16.0) 13 (13.5) 0.82(0.41~1.63) 隐性 T/T-T/C 223(99.1) 96(100.0) 0.859 0.230 1.00 C/C 2 (0.9) 0 (0.0) --- 超显性 T/T-C/C 191(84.9) 83 (86.5) 0.133 0.710 1.00 T/C 34 (15.1) 13 (13.5) 0.88(0.44~1.75) 逻辑累加 --- --- --- 0.519 0.470 0.79(0.41~1.51) *P < 0.05。 2.3 rs57095329潜在作用机制探讨
在得到rs57095329可能与CIN风险相关的初步结果后,笔者对其潜在作用机制进行了探讨。分析结果显示,rs57095329可能影响与43种蛋白结合而发挥作用,分数为2 b,见图1。在子宫组织中,rs57095329对SPI1的表达有影响(P = 0.050),见图2。
3. 讨论
目前研究已经确定持续感染高危型HPV是导致宫颈癌最主要的原因[24]。病毒、病毒和环境因素之间复杂的相互作用决定了最终的感染结局[25]。宿主癌基因或肿瘤抑制基因的分子改变是癌症发展的关键驱动因素[26]。研究表明miRNA具有抗癌基因或致癌基因的功能[27−28]。miRNA中的SNP与包括癌症在内的多种疾病密切相关,其中的某些SNP与人类癌症的易感性存在相关性[18,29−31]。miRNA中的SNP可改变miRNA的表达水平或结合亲和力,进而调控靶基因的表达,从而显著影响细胞通路调控决定细胞命运[32−33]。本研究对宫颈癌组织与正常组织差异表达的miR-416a基因进行分析,选取位于该基因启动区域上的SNP位点rs57095329、rs6864584进行基因分型。
在本研究中,rs57095329位点的等位基因A可能是CIN发生的保护性因素。在显性模式下,携带G等位基因(A/G-G/G)的个体CIN发生风险显著升高。由于miR-146a通过负反馈回路诱导基因产物下调,对先天和适应性免疫应答会诱导免疫调节作用,该miRNA上的SNP位点与癌症治疗预防相关[34]。SNP位点rs57095329、rs6864584位于miR-146a的启动子区域,rs57095329位于序列上游17 kb处[35−36]。Luo等[35]发现rs57095329在亚洲人中与风险相关的G等位基因与红斑性狼疮患者中的miR-146a的表达降低有关。Hefzy等[37]在慢性乙型肝炎病毒感染患者中发现rs57095329 病例组中AG和GG基因型的频率高于对照组,表明rs57095329可作为埃及人群中慢性乙型肝炎病毒致病的潜在危险因素。笔者研究发现与上述结果一致。对rs57095329相关作用机制的初步探索结果提示,rs57095329对子宫组织中SPI1的表达有影响。转录因子 SPI1通常通过调控细胞分化而发挥作用,也可作为单体与不同的DNA结合蛋白协同作用以激活靶基因转录[38]。2019年,Tao 等[39]研究结果证实,SPI1可以被募集到PARP9启动子并激活其转录来促进宫颈癌细胞的侵袭性。此外,笔者注意到不同研究中对rs57095329的分析结果存在不一致性,这可能是由于rs57095329位点在不同疾病中发挥的作用并不统一。各研究中人群的遗传背景差异以及样本量的不同也是值得考虑的重要影响因素。因此,未来需要对rs57095329位点进行更深入地研究,明确其与不同疾病之间的联系。
本次研究发现rs6864584位点与CIN的发生风险无相关性。部分研究所得结果与笔者的结果是一致的。Zhang等[40]在自身免疫性疾病易感性研究中表明,rs6864584 C等位基因与患病风险的降低相关,但并不显著。此外,部分研究与笔者的结果不一致。例如Salimi等[41]发现rs6864584 T/C降低了TT+TC下自发性反复自然流产的风险。rs6864584在不同研究中存在差异,未来需要在更丰富的队列中进行进一步的考察。
综上所述,本研究发现在CIN组中miR-146a基因SNP位点rs57095329的A等位基因分布频率低于对照组,推测该等位基因可能降低了CIN的发生风险。rs6864584位点与CIN的发生风险无相关性。miR-146a-rs57095329在宫颈癌中发挥的功能有必要进行更深入的研究。
-
图 9 灯盏乙素对PI3K/AKT/mTOR信号通路及炎症因子的影响[(
$\bar x \pm s $ ),n = 6)]A-a:PI3K的mRNA 表达;A-b:AKT的mRNA 表达;A-c:mTOR的mRNA 表达。B-a:Bax的mRNA 表达;B-b:Bcl-2的mRNA 表达;B-c:Caspase-3的mRNA 表达。C: 小鼠血清TNF-α的含量。D:小鼠血清IL-1β的含量。与空白对照组相比较,*P < 0.05,** P < 0.01;与模型组比较,#P < 0.05,#P < 0.01。
Figure 9. Effects of Scutellarin on PI3K/AKT/mTOR signaling pathway-related genes and inflammatory factors[(
$\bar x \pm s $ ),n = 6)]表 1 实验涉及的部分数据库
Table 1. Main databases involved in the experiment
名称 网址 用途 Pubchem https://pubchem.ncbi.nlm.nih.gov/ 化合物相关信息下载 Swiss http://www.swisstargetprediction.ch/ 药物相关靶点检索 Pharm mapper http://www.lilab-ecust.cn/pharmmapper/ 药物相关靶点检索 Genecards https://www.genecards.org/ 疾病相关靶点检索 Disgenet https://www.disgenet.org/ 疾病相关靶点检索 Venny https://bioinfogp.cnb.csic.es/tools/venny/ 疾病-药物靶点交集图 Uniprot https://www.uniprot.org/ 靶点名称翻译 David https://david.ncifcrf.gov/ KEGG、GO分析 String https://string-db.org/ 蛋白质相互作用网络 微生信平台 http://www.bioinformatics.com.cn/ 数据处理 Cytoscape v3.9.0软件 https://cytoscape.org/ 数据处理及绘图 表 2 引物序列
Table 2. The primer sequence
引物名称 长度(bp) 引物序列5′→3′ PI3K 20 Forward primer:CTTATGTCCTTGGCATTGGT Reverse primer:CAAAGTCTATGTGGAAGAGCTG AKT1 20 Forward primer:TCACCCAGTGACAACTCAG Reverse primer:AAACTCGTTCATGGTCACAC Bcl-2 19 Forward primer:TGACTGAGTACCTGAACCG Reverse primer:TAGTTCCACAAAGGCATCC Bax 19 Forward primer:CTGCAGAGGATGATTGCTG Reverse primer:GTCTGCAAACATGTCAGCT Caspase-3 21 Forward primer:AACTCTTCATCATTCAGGCCT Reverse primer:CCATATCATCGTCAGTTCCAC β-actin 20 Forward primer:CTTCCTGGGTATGGAATCCT Reverse primer:TCTTTACGGATGTCAACGTC mTOR 20 Forward primer:GCAGATTTGCCAACTATCTTCGG Reverse primer:CAGCGGTAAAAGTGTCCCCTG -
[1] Sun J J,Yin X W,Liu H H,et al. Rapamycin inhibits ox-LDL-induced inflammation in human endothelial cells in vitro by inhibiting the mTORC2/PKC/c-Fos pathway[J]. Acta Pharmacol Sin,2018,39(3):336-344. doi: 10.1038/aps.2017.102 [2] Troidl K,Schubert C,Vlacil A K,et al. The lipopeptide MALP-2 promotes collateral growth[J]. Cells,2020,9(4):997. doi: 10.3390/cells9040997 [3] Xi J Y,Rong Z. Scutellarin ameliorates high glucose-induced vascular endothelial cells injury by activating PINK1/Parkin-mediated mitophagy[J]. J Ethnopharmacol,2021,271:113855. doi: 10.1016/j.jep.2021.113855 [4] Keeter,W C,Ma S,Stahr N,et al. Atherosclerosis and multi-organ-associated pathologies[J]. Semin Immunopathol,2022,44(3):363-374. doi: 10.1007/s00281-022-00914-y [5] Wu M F,Xu K Z,Guo YG,et al. Lipoprotein(a) and atherosclerotic cardiovascular disease:Current understanding and future perspectives[J]. Cardiovasc Drugs Ther,2019,33(6):739-748. doi: 10.1007/s10557-019-06906-9 [6] Panes O,Gonzalez C,Hidalgo P,et al. Platelet tissue factor activity and membrane cholesterol are increased in hypercholesterolemia and normalized by rosuvastatin,but not by atorvastatin.[J]. Atherosclerosis,2017,257:164-171. doi: 10.1016/j.atherosclerosis.2016.12.019 [7] Lee Y P,Cho Y,Kim E J,et al. Reduced expression of pyruvate kinase in kidney proximal tubule cells is a potential mechanism of pravastatin altered glucose metabolism[J]. Sci Rep,2019,9(1):5318. doi: 10.1038/s41598-019-39461-2 [8] Ott B R,Daiello L A,Dahabreh I J,et al. Do statins impair cognition? A systematic review and meta-analysis of randomized controlled trials[J]. J Gen Intern Med,2015,30(3):348-358. doi: 10.1007/s11606-014-3115-3 [9] Mo J,Yang R H,Li F,et al. Scutellarin protects against vascular endothelial dysfunction and prevents atherosclerosis via antioxidation[J]. Phytomedicine,2018,42:66-74. doi: 10.1016/j.phymed.2018.03.021 [10] Zhou S,Ai Z,Li W,et al. Deciphering the pharmacological mechanisms of taohe-chengqi decoction extract against renal fibrosis through integrating network pharmacology and experimental validation In vitro and in vivo[J]. Front Pharmacol,2020,11(1):425. [11] 何信用,王俊岩,宋囡,等. 基于网络药理学的葛根素治疗动脉粥样硬化潜在分子机制研究[J]. 中华中医药学刊,2020,38(9):116-120. doi: 10.13193/j.issn.1673-7717.2020.09.029 [12] Fernandez-Hernando C,Jozsef L,Jenkins D,et al. Absence of akt1 reduces vascular smooth muscle cell migration and survival and induces features of plaque vulnerability and cardiac dysfunction during atherosclerosis[J]. Arterioscler Thromb Vasc Biol,2009,29(12):2033-2040. doi: 10.1161/ATVBAHA.109.196394 [13] 李敏静,郭莉,陈晔,等. 血府逐瘀汤含药血清对低氧诱导的大鼠肺动脉平滑肌细胞增殖及PI3K/AKT/mTOR信号通路的影响[J]. 中国中西医结合杂志,2020,40(7):836-841. [14] 张莉,白林英,母凯茜,等. 灯盏乙素对小型猪动脉粥样硬化模型颈总动脉斑块组织中基质金属蛋白酶-1、2、9 表达的影响[J]. 中药药理与临床,2020,36(6):91-97. [15] Ashino T,Yamamoto M,Yoshida T,et al. Redox-sensitive transcription factor Nrf2 regulates vascular smooth muscle cell migration and neointimal hyperplasia[J]. Arterioscler Thromb Vasc Biol,2013,33(4):760-768. doi: 10.1161/ATVBAHA.112.300614 -