-
乳腺癌是发生在乳腺上皮组织的恶性肿瘤，居女性癌症发病率的首位，是引起女性死亡的主要原因之一[1]。乳腺癌最常用的联合化疗可降低复发率，但引起的肿瘤细胞多药耐药和严重的不良反应限制了其临床疗效。中医药已被广泛用于癌症的治疗，可通过提高肿瘤对药物的敏感性、调节机体免疫功能，增强化疗效果，改善患者预后[2]。其优势主要体现在中药作用多靶点、多途径而产生广泛的生物活性和较低的毒性，现已被用作乳腺癌在内的多种癌症的替代疗法[3]。
苦参（Sophora flavescens）味苦性寒，始载于《神农本草经》，具有清热燥湿、杀虫、利尿的功效。苦参的主要化学成分包括生物碱、黄酮类、苯丙素、萜类、甾体、有机酸和挥发油等。其中苦参生物碱占3%，苦参黄酮占1.5%，是苦参的两类特征性成分，具有抗炎和抗肿瘤等药理作用，对急性髓系白血病、鼻咽癌、宫颈癌、乳腺癌、结直肠癌、膀胱癌、胰腺癌等均显示出较好的抗肿瘤活性[4]。
目前关于苦参抗乳腺癌的研究多局限于探索苦参中某种活性成分作用于特定靶点或通路的抗乳腺癌的机制，未从整体角度阐释苦参抗乳腺癌的作用机制。因此，本研究基于网络药理学的方法，从中药作用的整体理论和中医药现代化的研究思路，探究苦参抗乳腺癌的作用，初步构建其主要活性成分作用于多靶点、多通路的作用机制网络，为开发基于苦参的抗乳腺癌新药和临床用药提供参考[5-6]。
-
中药系统药理学分析平台（TCMSP）；中医药百科全书数据库（ETCM）；通用蛋白质数据库（UniProt）；人类基因数据库（Genecards）；生物信息学和化学信息学数据库（Drugbank）；人类孟德尔遗传数据库（OMIM）；治疗靶点数据库（TTD）；蛋白质互作平台（STRING10.5）；生物学信息注释数据库 （DAVID 6.8）。数据库检索日期为2021年12月。
-
应用TCMSP数据库检索关键词“苦参”的所有活性成分，选择“ingredient”，设置阈值标准为：口服生物利用度（OB）≥30% ，且类药性（DL）≥0.18，并结合ETCM数据库和国内外文献，筛选出苦参的活性成分。
-
从PubChem 数据库获得筛选出的苦参活性成分的2D结构式导入Swiss Target Prediction数据库，预测苦参活性成分的靶点，物种选择Homo sapiens，导入结构式，筛选与苦参活性化合物相关的作用靶点。通过 UniProt数据库，将检索得到的所有靶点蛋白校正为其标准名称。
-
将搜索关键词设置为“breast cancer”，使用GeneCards、Drugbank、OMIM、TTD数据库对癌症的靶点进行搜索，将Drugbank、OMIM、TTD数据库中的乳腺癌相关基因进行合并、去重后，得到乳腺癌的靶点[7]。
-
将以上四个数据库获取的乳腺癌靶基因与苦参活性成分靶基因输入Venny 2.1.0软件，取交集作为苦参抗乳腺癌的靶点，筛选出苦参抗乳腺癌的潜在靶点。
-
将上述靶点及主要活性成分导入 Cytoscape 3.7.2 软件，构建苦参抗乳腺癌的苦参-核心活性成分-靶点网络图。用Cytoscape 3.7.2 软件中“Networkanalyzer”功能，对网络拓扑属性进行分析，得到苦参抗乳腺癌的主要活性成分。
-
将上述得到的主要靶基因上传到STRING数据库，设置物种为“Homo sapiens”，进行检索，利用STRING数据库对PPI网络图进行分析，得到相关靶点参数。将该参数上传至Cytoscape 3.7.2 软件，构建苦参蛋白互作网络（PPI）图。
-
将苦参抗乳腺癌靶点输入DAVID 6.8平台，结合Hiplot平台分别进行GO功能富集分析与KEGG通路富集分析，依据P<0.01，筛选苦参抗乳腺癌的靶点富集的生物功能及信号通路。
-
根据苦参抗乳腺癌的苦参-核心活性成分-靶点网络图，选择核心活性成分与主要靶点进行分子对接。从PubChem数据库下载苦参核心活性成分的SDF文件，在PDB数据库寻找各主要靶点的PDB文件，利用Schrodinger 2018软件进行分子对接。
-
人乳腺癌MCF-7细胞购自中国科学院上海细胞库。以含10% FBS的DMEM高糖培养基添加100 U/ml青霉素和100 µg/ml链霉素，5% CO2的37˚C 培养箱培养，取对数生长期的细胞用于后续实验。
-
苦参提取物（原药材生药量=20∶1，批号：20220710，西安瑞迪生物科技有限公司）；胎牛血清（批号：10099-141）、胰蛋白酶（批号：2323491）购自GIBCO；DMEM高糖培养基（批号：8122184，Hyclone）；CCK8试剂盒（批号：C0039）、Hoechst 33342荧光染料（批号：C1029）购自碧云天生物技术有限公司；4%多聚甲醛（批号：CR2110036）、Androgen Receptor（AR）抗体（批号：GB11253）、cy3-山羊抗兔二抗（批号：GB21303）、DAPI染料（批号：G1012）、抗荧光淬灭封片剂（批号：G1401）购自Servicebio。
-
取对数生长期的MCF-7细胞，调整细胞密度至1×105/ml，接种于96孔板中，每孔100 μl，培养过夜待细胞贴壁后，加入不同浓度的苦参提取物（0、1.25、2.5、5、10、20、40 mg/ml），每组设置6个复孔。空白组加入相应的无细胞培养基，相同处理。药物作用24 h后，显微镜下观察并拍照；吸弃原培养基，PBS洗涤2次，每孔加入100 μl含有10 μl的CCK8的新鲜培养基，于5% CO2的37˚C 培养箱中继续培养2 h后，多功能酶标仪测定450 nm波长处OD值，根据OD值计算细胞的存活率：细胞存活率=[（给药组OD值−空白组OD值）/（对照组OD值−空白组OD值）]×100%。
-
取对数生长期的MCF-7细胞，调整细胞密度至1×105/ml，接种于6孔板中，培养过夜待细胞贴壁后，加入不同浓度的苦参提取物处理，药物作用48 h后，吸弃原培养基，PBS洗涤2次，每孔加入1 ml含有100 μl Hoechst 33342染色液的新鲜培养基，室温染色20 min后，PBS洗涤3次，荧光显微镜下观察并拍照，并使用Image J软件进行荧光定量分析。
-
取对数生长期的MCF-7细胞，调整细胞密度至1×105/ml，接种于6孔板细胞爬片上，培养过夜待细胞贴壁后，加入不同浓度的苦参提取物处理，药物作用48 h后，吸弃原培养基，PBS洗涤3次，每孔加 1 ml 4%的多聚甲醛，固定细胞15 min，弃去固定液，PBS 洗涤 3 次，加入细胞破膜液破膜30 min，PBS洗涤3次，每次5 min。每孔加入 1 ml 的 5% BSA溶液，置于摇床上室温封闭 30 min，弃去封闭液，加入 AR一抗（1∶500稀释），4 ℃孵育过夜，PBS洗涤3次，每次5 min。 加入荧光素标记的cy3-山羊抗兔二抗，摇床上室温孵育 1 h， PBS 洗涤 3 次，每次5 min。 加入1 μg/ml DAPI，室温孵育10 min，PBS 洗涤 3 次，每次5 min。 待玻片稍干后，用抗荧光淬灭封片剂封片。荧光显微镜下观察并分别采集AR染色、DAPI染色及两者合并图像，并使用Image J软件进行荧光定量分析，每孔至少随机选取3个不同视野。
-
在TCMSP数据库中检索苦参共得到113个化合物，根据OB≥30%、DL≥0.18，结合ETCM数据库和国内外文献，最终确定结构明确的化合物36种，见表1。
表 1 苦参中36种活性化合物基本信息
分子编号 化合物名称 分子式 OB(%) DL MOL006596 大豆素 C20H18O5 97.27 0.76 MOL000456 菜豆素 C20H18O4 78.20 0.73 MOL001484 高丽槐素 C16H12O5 75.18 0.54 MOL004941 甘草素 C15H12O4 71.12 0.18 MOL000392 刺芒柄花素 C16H12O4 69.67 0.21 MOL006565 异苦参碱 C15H24N2O 68.68 0.25 MOL003648 山槐素 C16H12O5 65.83 0.54 MOL003627 槐果碱 C15H22N2O 64.26 0.25 MOL005944 苦参碱 C15H24N2O 63.77 0.25 MOL006627 12,13-去氢苦参碱 C15H22N2O 62.23 0.25 MOL006571 安纳基林 C15H20N2O 62.01 0.24 MOL006626 勒奇黄烷酮 G C20H20O6 60.97 0.40 MOL003680 槐定碱 C15H24N2O 60.07 0.25 MOL006566 12,13-去氢苦参碱 C15H22N2O 58.34 0.25 MOL006649 5α-羟基苦参碱 C15H24N2O2 55.42 0.28 MOL006630 去甲藤黄素 C15H10O6 54.93 0.24 MOL006628 白羽扇豆碱 C15H24N2O 52.71 0.24 MOL006623 苦参新醇 T C25H30O7 51.28 0.64 MOL006650 三叶豆紫檀苷丙二酸酯 C25H24O13 48.69 0.52 MOL006652 三叶豆紫檀苷 C22H22O11 48.53 0.74 MOL006604 异黄腐醇 C21H22O5 48.09 0.39 MOL005100 橙皮素 C16H14O6 47.74 0.27 MOL006613 苦参素 C16H14O5 47.62 0.38 MOL000098 槲皮素 C15H10O7 46.43 0.28 MOL003347 贯叶金丝桃素 C35H52O4 44.03 0.60 MOL003673 怀特酮 C20H18O5 42.80 0.36 MOL006622 苦参新醇 O C27H30O13 42.41 0.76 MOL001040 柚皮素 C15H12O5 42.36 0.21 MOL003676 槐胺碱 C15H20N2O 42.16 0.25 MOL006582 5α,9α-二羟基苦参碱 C15H24N2O3 40.93 0.32 MOL003542 8-异戊烯醇去甲基淫羊藿苷 C20H18O6 38.04 0.39 MOL006569 14β-羟基苦参碱 C15H24N2O2 37.26 0.29 MOL000006 木犀草素 C15H10O6 36.16 0.25 MOL006561 14α-羟基苦参碱 C15H24N2O2 35.73 0.29 MOL006570 9α-羟基槐胺碱 C15H22N2O2 35.23 0.29 MOL006563 9α-羟基苦参碱 C15H24N2O2 32.04 0.29 -
通过Swiss Target Prediction数据库预测共得到苦参活性成分的药物靶点509个，分别从TTD数据库、Drugbank数据库和OMIM数据库获取乳腺癌疾病靶点137个、244个和506个，去除重复项得到756个疾病靶点。从GeneCards数据库获得乳腺癌疾病靶点15230个，保留分值≥4.38（2倍中位数）的疾病靶点3813个，合并4个数据库检索结果并删去重复项，得到275个乳腺癌相关靶点。采用Venny 2.1软件，绘制苦参活性成分靶点与乳腺癌疾病靶点韦恩图，得到70个共有靶点（图1）。
图 1 苦参抗乳腺癌靶点韦恩图
-
根据Cytoscape 3.7.2软件绘制的苦参抗乳腺癌“中药-活性成分-靶点”网络图的参数，筛选化合物Degree（度）值>5的节点，共筛选出24种核心活性成分，对应67个靶点，构建苦参抗乳腺癌的“中药-核心活性成分-靶点”网络图，如图2所示，每个参与互作的蛋白节点用一个圆圈表示，圆圈面积越大、颜色越深表示其Degree值越大，在药物-核心活性成分-靶点网络中，总共包括92个节点（其中包含1个药物节点、24个化合物节点和67个靶点节点）和 630 条边。按照Degree（度）值>5排列活性成分，从中可见木犀草素、槲皮素、菜豆素、甘草素等的Degree值最大，是苦参抗乳腺癌的核心活性成分。
图 2 苦参-核心活性成分-靶点网络图
为了更好地理解苦参抗乳腺癌的作用机制，按照Degree>5排列主要靶点，筛选出苦参抗乳腺癌的主要疾病靶点25个。用STRING数据库分析该网络参数，导入Cytoscape软件构建苦参蛋白互作网络（PPI）图，以明确生物分子间相互作用相关的重要靶蛋白（图3）。筛选出前12位的主要核心靶蛋白为EGFR、AKT1、ESR1、SRC、CYP19A1、AR、ABCB1、CDK1、ESR2、KDR、IGF1R和TOP1，各核心靶点的拓扑参数见表2。
图 3 苦参抗乳腺癌的PPI网络（筛选后）
表 2 核心靶点及其拓扑参数
靶点 蛋白名称 度 介数 紧密度 EGFR 表皮生长因子受体 16 1.28×10−1 7.86×10−1 AKT1 丝氨酸/苏氨酸蛋白激酶 1 16 1.28×10−1 7.86×10−1 ESR1 雌激素受体 α 15 1.45×10−1 7.59×10−1 SRC 非受体酪氨酸激酶 Src 14 1.58×10−1 7.33×10−1 CYP19A1 细胞色素P450酶 19A1 11 0.95×10−1 6.47×10−1 AR 雄激素受体 10 0.32×10−1 6.47×10−1 ABCB1 ATP结合盒B亚家族成员1 9 0.65×10−1 6.11×10−1 CDK1 周期素依赖性激酶1 8 0.97×10−1 5.95×10−1 ESR2 雌激素受体β 8 0.18×10−1 5.95×10−1 KDR 血管内皮生长因子受体2 8 0.12×10−1 5.95×10−1 IGF1R 胰岛素样生长因子1受体 8 0.08×10−1 5.95×10−1 TOP1 DNA拓扑异构酶1 6 0.08×10−1 5.64×10−1 -
利用DAVID数据库对筛选出的12个关键靶点进行GO 功能富集分析，依据P<0.01，选出GO条目49个，其中生物过程（BP）29个，分子功能（MF）17个，及细胞组成（CC）3个。利用Hiplot在线平台对每个种类的前15个条目进行可视化分析。结果显示，苦参抗乳腺癌靶点的生物过程主要涉及蛋白自磷酸化、凋亡过程的负向调控、信号转导、细胞增殖、胞内雌激素受体信号通路等。分子功能涉及酶结合反应、一氧化氮合成酶调节活性、蛋白酪氨酸激酶活性、染色质结合反应、类固醇结合反应等。细胞过程涉及细胞核、细胞质膜和细胞膜（图4）。
图 4 苦参抗乳腺癌关键靶点的GO富集分析
-
采用Hiplot在线平台对苦参关键作用靶点的KEGG信号通路进行富集分析，依据P<0.01筛选KEGG信号通路，以Gene Ratio代表目标基因中属于这个pathway的基因的数量在背景基因中所占的百分比，P 值以及调控通路上的基因个数来衡量 KEGG 富集程度，其中气泡颜色代表 P 值的大小，Gene Ratio代表富集的程度大小，气泡大小代表该通路中目标基因的数量，对前15个通路进行可视化分析，同时绘制苦参抗乳腺癌的基因-概念网络图。图5所示，按照P值从小到大富集的通路为内分泌抵抗、EGFR酪氨酸激酶抑制剂耐药、肿瘤中的蛋白多糖、雌激素信号通路、乳腺癌、催乳素信号通路、黏着斑、Rap1信号通路、前列腺癌、VEGF信号通路、黏附连接、黑色素瘤、神经胶质瘤、ErbB信号通路、Ras信号通路；各信号通路对应的靶点包括SRC、IGF1R、ESR2、ESR1、EGFR、AKT1、KDR和AR（图6）。
图 5 苦参抗乳腺癌关键靶点的KEGG通路富集分析
图 6 苦参抗乳腺癌关键靶点的基因-概念网络图
-
将24种核心活性成分的分子配体与12个靶点进行分子对接，得到275组分子对接数据。各靶点与核心活性成分对接能力的交叉数据如图7所示，图中红色越深，表示结合自由能越高，对接能力越强。结果显示活性成分与靶点EGFR、ESR1、CYP19A1、AR和ESR2等的对接能力较好。苦参中三种活性成分：甘草素、柚皮素和怀特酮与关键靶点AR结合的对接分数均大于10，结合能力较其他活性成分和靶点更强。
图 7 分子对接分数的热图分析
-
CCK8检测结果显示，MCF-7细胞的存活率随着苦参提取物（KS）浓度的增大而降低，呈现药物浓度依赖性，当KS浓度达到20 mg/ml时，MCF-7细胞的存活率约为49%（图8A）；绘制KS对MCF-7细胞作用的IC50曲线，计算IC50值约为17.3 mg/ml（图8B）。
图 8 CCK8检测MCF-7细胞增殖能力
-
采用Hoechst33342荧光染色观察细胞凋亡形态变化，空白对照组的MCF-7细胞生长状态良好，细胞膜和细胞核形态正常，细胞核呈现均匀的淡蓝色弥散的荧光，未观察到明显凋亡的细胞形态；KS各浓度组Hoechst33342荧光强度较空白对照组显著增强（P<0.01），细胞形态发生变化，细胞数量略有减少，细胞呈不规则膨胀或皱缩，部分细胞膜破裂，细胞核呈碎块状致密浓染，且当KS浓度为20 mg/ml时，细胞凋亡特征最为明显，凋亡数量明显增加，荧光强度最强呈现蓝白色，表明 KS可诱导乳腺癌MCF-7细胞凋亡（图9）。
图 9 Hoechst33342荧光染色检测MCF-7细胞凋亡（×200）
-
为验证苦参作用靶点的分子对接结果，采用免疫荧光染色检测MCF-7细胞内AR蛋白定位和表达水平。采用DAPI对细胞核进行染色，发现AR定位在细胞核内。与对照组相比，20 mg/ml 的KS浓度组，AR荧光强度显著最强（P<0.01），提示KS可显著上调AR蛋白的表达水平，与分子对接的结果一致（图10）。
图 10 免疫荧光染色观察MCF-7细胞AR表达水平（×200）
-
乳腺癌在《妇人大全良方》《本草纲目》等中医古籍属“乳岩”“乳石痈”的范畴，由于肝伤失其条达，气血瘀滞，热毒结滞于乳中而发病。热毒是导致恶性肿瘤发生发展的主要原因，脏腑经络、气血津液受热毒侵犯是恶性肿瘤发生发展的关键[8]。扶正和祛邪是中药在肿瘤综合治疗过程中的两个重要环节，祛邪类药物主要以清热解毒、软坚散结、活血化瘀等功效为主，具有诱导细胞分化、促进细胞凋亡、逆转多药耐药、抗侵袭转移等多方面的作用[9]。因此，清热解毒属祛邪的范畴，可作为乳腺癌的主要治法之一。苦参具有清热燥湿、杀虫、利尿之功效，在传统医学中主要用于治疗湿热黄疸、泻痢、皮肤瘙痒等病症。现代药理学研究表明，苦参具有抗氧化、抗炎、抗心律失常、平喘、抗溃疡和抗肿瘤等多种药理活性[10]。在抗肿瘤方面的研究显示，苦参及其制剂在抗肿瘤、减轻化疗不良反应和逆转肿瘤耐药方面均发挥重要作用。
本研究构建的苦参-核心活性成分-靶点网络发现苦参抗乳腺癌的核心成分有24种，主要为黄酮类和生物碱类。靶点EGFR、AKT1、ESR1、SRC、CYP19A1、AR、ABCB1、CDK1、ESR2、KDR、IGF1R和TOP1在苦参抗乳腺癌中发挥了重要作用。EGFR 属于受体酪氨酸激酶的 ErbB 家族成员，几乎在所有亚型的乳腺癌患者中常过度表达[11]。本研究发现苦参最核心的作用靶点为EGFR。EGFR表达增加与乳腺癌的发生发展有关，可通过激活Ras-Raf-MEK-ERK、PI3K-Akt-mTOR 和 Src-STAT3 通路来发挥作用[12-13]。有研究显示EGFR介导的β-连环蛋白核蓄积对于AKT1抑制诱导的乳腺癌转移至关重要[14]。AKT1属于AKT家族成员，是本研究显示的苦参作用的关键靶点之一。有研究表明AKT1和细胞周期蛋白D1（Cyclin D1）之间存在交叉调节，两者在不同亚型乳腺癌中的表达呈正相关性，且在Luminal A、Luminal B和三阴性乳腺癌中的相关性最强[15-16]。基于二代测序分析显示，AKT1基因突变与原发性乳腺癌相比在转移性乳腺癌中显著富集，并被明确为可调控的靶点[17-18]。目前，已有多种靶向AKT的小分子抑制剂正在进行临床试验，可能为包括三阴性乳腺癌在内的难治性乳腺癌患者带来新的希望[19]。但AKT1抑制乳腺癌转移的具体机制和下游信号通路尚不完全清楚[20-21]。而单纯的AKT抑制剂作为抗癌药物可能带来争议，如药物开发难度大价格昂贵，对靶点的选择特异性问题导致的不良反应包括高血糖、腹泻、皮疹及可能出现的其他不利因素[22-23]。而AKT抑制剂 Capivasertib（AZD5363）在AKT1突变肿瘤患者中表现出显著的临床意义，特别是对于ER阳性的乳腺癌患者[24]；同样，AKT抑制剂GDC-0068在三阴性乳腺癌中与紫杉醇联用也延长了患者的无进展生存期[25]。本研究中显示的苦参中核心活性成分对AKT1的靶向性可能为乳腺癌患者治疗药物的开发带来新的思路。
内分泌治疗在ER阳性乳腺癌患者的治疗中占有举足轻重的地位，该疗法包括选择性雌激素受体调节剂（SERMs），选择性雌激素受体抑制剂（SERDs）和芳香化酶抑制剂（AIs）。ESR1基因是编码ER的基因，在ER阳性乳腺癌中存在ESR1基因突变、扩增和融合，但扩增和融合的频率远低于点突变且与耐药的相关性仍不明确[26]。药物靶点本身的突变通常是肿瘤逃避药物作用的优先机制，ESR1突变已成为内分泌抵抗的关键机制，近年来研究发现ESR1配体结合域（LBD）的突变是 ER阳性转移性乳腺癌耐药的常见驱动因素[27]。这些LBD点突变可通过不依赖激素的ER转录活性，直接导致乳腺癌对AI的抗性并降低对他莫昔芬和氟维司群的敏感性[28]。而雄激素受体（AR）在ER阳性的乳腺癌中的作用一直以来都存在争议，最新的研究显示，AR通过改变ER和关键转录共激活因子的基因组分布，在内分泌抵抗的ER阳性乳腺癌中发挥抑癌作用，并支持AR激动剂作为ER阳性乳腺癌的治疗策略[29]。此外，研究者采用中药方剂治疗三阴性乳腺癌患者，发现AR阳性的三阴性乳腺癌患者预后优于AR阴性患者[30]。但目前尚无AR激动剂被用于乳腺癌的治疗，且苦参作为具有多靶点作用的中药，与AR特异性抑制剂的作用机制存在一定差异。本研究发现内分泌抵抗是苦参抗乳腺癌最核心的通路，分子对接的结果显示主要活性成分与靶点AR的结合能最高。在体外细胞实验中，发现苦参提取物可显著抑制ER表达阳性的人乳腺癌MCF-7细胞增殖，诱导细胞凋亡，上调AR蛋白的表达，进一步验证了分子对接的结果，提示苦参对AR的上调可能成为治疗乳腺癌新的治疗策略，为乳腺癌新的治疗方案提供理论依据。
综上所述，本研究通过网络药理学发现苦参核心成分的主要靶点包括AKT1、ESR1、CYP19A1和AR等，提示苦参可能通过多靶点增强乳腺癌对内分泌治疗的敏感性。分子对接的结果显示苦参核心成分与AR的结合能力最强，体外实验也验证了苦参对AR表达的上调作用，提示苦参对AR的上调可能成为治疗乳腺癌新的治疗策略。但本研究还存在一定的局限性，如本研究以单味中药苦参作为研究对象分析其抗乳腺癌的活性成分和作用机制，而以苦参为主要组分的中药复方或中成药在临床中的应用更为常见，其他辅助中药也可能会发挥协同抗乳腺癌的作用。此外，本研究以网络药理学的方法分析苦参抗乳腺癌的作用机制，这些研究结果还需要更多的基础实验和临床研究来进一步验证。
Potential mechanism of Sophora flavescens against breast cancer via network pharmacology and molecular docking
-
摘要:
目的 采用网络药理学和分子对接的方法分析苦参抗乳腺癌的主要活性成分及潜在分子机制。 方法 采用TCMSP、ETCM数据库和查阅国内外文献收集筛选苦参化学成分，通过Swiss Target Prediction 数据库对有效活性成分的靶点进行预测，通过GeneCards、TTD、Drugbank、OMIM收集乳腺癌相关靶点，根据Venny 2.1.0软件筛选出苦参抗乳腺癌疾病靶点；采用Cytoscape软件构建苦参-核心活性成分-靶点网络图，用 STRING数据库分析共有靶点，构建PPI网络图，通过DAVID数据库和Hiplot平台对关键靶点蛋白进行GO功能富集分析和KEGG通路富集分析，并采用Schrodinger软件对活性成分与靶点进行分子对接，同时采用分子生物学方法对关键靶点进行验证。 结果 从苦参中共筛选出结构明确的活性成分36个，筛选出70个苦参抗乳腺癌的疾病靶点，其中EGFR、AKT1、ESR1、SRC、CYP19A1、AR、ABCB1等12个核心靶点参与乳腺癌中内分泌抵抗、EGFR酪氨酸激酶抑制剂耐药、雌激素等信号通路，且核心成分与关键靶点AR之间的结合能最高。体外细胞实验显示，苦参提取物可抑制乳腺癌细胞的增殖，诱导细胞凋亡，上调AR蛋白的表达。 结论 苦参通过多成分作用于多靶点和多个信号通路，调节复杂的生物过程发挥抗乳腺癌的作用，苦参对AR蛋白的上调可能成为治疗乳腺癌新的治疗策略。 Abstract:Objective To analyze the main active components and potential molecular mechanism of Sophora flavescens against breast cancer based on network pharmacology and molecular docking. Methods The chemical constituents were collected and screened by TCMSP, ETCM database and literature review. The targets of active ingredients were predicted by Swiss Target Prediction database. Breast cancer-related targets were collected by GeneCards, TTD, Drugbank and OMIM. The anti-breast cancer targets of Sophora flavescens were screened by Venny 2.1.0 software. Cytoscape software was used to construct the network diagram of Sophora flavescens-key active ingredients-targets. STRING database was used to analyze the common targets, and PPI network diagram was constructed. GO function enrichment analysis and KEGG pathway enrichment analysis of key target proteins were performed by DAVID database and Hiplot online platform. Schrodinger software was used to calculate the molecular docking between the active ingredients and targets. Molecular biological methods were used to verify the key targets. Results A total of 36 active components with clear structures were screened from Sophora flavescens. 70 anti-breast cancer targets of Sophora flavescens were screened out. 12 core targets including EGFR, AKT1, ESR1, SRC, CYP19A1, AR and ABCB1 participate in endocrine resistance, EGFR tyrosine kinase inhibitors and estrogen signaling pathways in breast cancer. Moreover, the docking score between the core component and the key target AR is the highest. In vitro experiments showed that the extract of Sophora flavescens can inhibit the proliferation of breast cancer cells, induce cell apoptosis and up-regulate AR protein expression. Conclusion It was revealed that Sophora flavescens plays an anti-breast cancer role by regulating complex biological processes through multiple components acting on multiple targets and signaling pathways. The upregulation of AR protein by Sophora flavescens may become a new therapeutic strategy for the treatment of breast cancer. -
Key words:
- Sophora flavescens /
- breast cancer /
- network pharmacology /
- molecular docking /
- androgen receptor
-
表 1 苦参中36种活性化合物基本信息
分子编号 化合物名称 分子式 OB(%) DL MOL006596 大豆素 C20H18O5 97.27 0.76 MOL000456 菜豆素 C20H18O4 78.20 0.73 MOL001484 高丽槐素 C16H12O5 75.18 0.54 MOL004941 甘草素 C15H12O4 71.12 0.18 MOL000392 刺芒柄花素 C16H12O4 69.67 0.21 MOL006565 异苦参碱 C15H24N2O 68.68 0.25 MOL003648 山槐素 C16H12O5 65.83 0.54 MOL003627 槐果碱 C15H22N2O 64.26 0.25 MOL005944 苦参碱 C15H24N2O 63.77 0.25 MOL006627 12,13-去氢苦参碱 C15H22N2O 62.23 0.25 MOL006571 安纳基林 C15H20N2O 62.01 0.24 MOL006626 勒奇黄烷酮 G C20H20O6 60.97 0.40 MOL003680 槐定碱 C15H24N2O 60.07 0.25 MOL006566 12,13-去氢苦参碱 C15H22N2O 58.34 0.25 MOL006649 5α-羟基苦参碱 C15H24N2O2 55.42 0.28 MOL006630 去甲藤黄素 C15H10O6 54.93 0.24 MOL006628 白羽扇豆碱 C15H24N2O 52.71 0.24 MOL006623 苦参新醇 T C25H30O7 51.28 0.64 MOL006650 三叶豆紫檀苷丙二酸酯 C25H24O13 48.69 0.52 MOL006652 三叶豆紫檀苷 C22H22O11 48.53 0.74 MOL006604 异黄腐醇 C21H22O5 48.09 0.39 MOL005100 橙皮素 C16H14O6 47.74 0.27 MOL006613 苦参素 C16H14O5 47.62 0.38 MOL000098 槲皮素 C15H10O7 46.43 0.28 MOL003347 贯叶金丝桃素 C35H52O4 44.03 0.60 MOL003673 怀特酮 C20H18O5 42.80 0.36 MOL006622 苦参新醇 O C27H30O13 42.41 0.76 MOL001040 柚皮素 C15H12O5 42.36 0.21 MOL003676 槐胺碱 C15H20N2O 42.16 0.25 MOL006582 5α,9α-二羟基苦参碱 C15H24N2O3 40.93 0.32 MOL003542 8-异戊烯醇去甲基淫羊藿苷 C20H18O6 38.04 0.39 MOL006569 14β-羟基苦参碱 C15H24N2O2 37.26 0.29 MOL000006 木犀草素 C15H10O6 36.16 0.25 MOL006561 14α-羟基苦参碱 C15H24N2O2 35.73 0.29 MOL006570 9α-羟基槐胺碱 C15H22N2O2 35.23 0.29 MOL006563 9α-羟基苦参碱 C15H24N2O2 32.04 0.29 
下载: 导出CSV
表 2 核心靶点及其拓扑参数
靶点 蛋白名称 度 介数 紧密度 EGFR 表皮生长因子受体 16 1.28×10−1 7.86×10−1 AKT1 丝氨酸/苏氨酸蛋白激酶 1 16 1.28×10−1 7.86×10−1 ESR1 雌激素受体 α 15 1.45×10−1 7.59×10−1 SRC 非受体酪氨酸激酶 Src 14 1.58×10−1 7.33×10−1 CYP19A1 细胞色素P450酶 19A1 11 0.95×10−1 6.47×10−1 AR 雄激素受体 10 0.32×10−1 6.47×10−1 ABCB1 ATP结合盒B亚家族成员1 9 0.65×10−1 6.11×10−1 CDK1 周期素依赖性激酶1 8 0.97×10−1 5.95×10−1 ESR2 雌激素受体β 8 0.18×10−1 5.95×10−1 KDR 血管内皮生长因子受体2 8 0.12×10−1 5.95×10−1 IGF1R 胰岛素样生长因子1受体 8 0.08×10−1 5.95×10−1 TOP1 DNA拓扑异构酶1 6 0.08×10−1 5.64×10−1
下载: 导出CSV
-
[1] SIEGEL R, MILLER K, WAGLE N S, et al. Cancer statistics, 2023[J]. CA Cancer J Clin, 2023, 73: 17-48. doi: 10.3322/caac.21763 [2] ZHANG X Y, QIU H, LI C S, et al. The positive role of traditional Chinese medicine as an adjunctive therapy for cancer[J]. Biosci Trends, 2021, 15(5): 283-298. doi: 10.5582/bst.2021.01318 [3] GREENLEE H, DUPONT-REYES M J, BALNEAVES L G, et al. Clinical practice guidelines on the evidence-based use of integrative therapies during and after breast cancer treatment[J]. CA Cancer J Clin, 2017, 67(3): 194-232. doi: 10.3322/caac.21397 [4] CHEN M H, GU Y Y, ZHANG A L, et al. Biological effects and mechanisms of matrine and other constituents of Sophora flavescens in colorectal cancer[J]. Pharmacol Res, 2021, 171: 105778. doi: 10.1016/j.phrs.2021.105778 [5] 李梢. 基于生物网络调控的方剂研究模式与实践[J]. 中西医结合学报, 2007, 5(5): 489-493. [6] 林玲, 吴勤研, 宋熙薇. 基于网络药理学的康艾注射液抗非小细胞肺癌作用机制分析[J]. 药学与临床研究, 2021, 29(2): 86-91. [7] WANG Y X, ZHANG S, LI F C, et al. Therapeutic target database 2020: enriched resource for facilitating research and early development of targeted therapeutics[J]. Nucleic Acids Res, 2020, 48(D1): D1031-D1041. [8] 李华华, 杨峰, 裴俊文, 等. 杨峰教授在清热解毒理论指导下运用复方苦参注射液治疗恶性肿瘤经验[J]. 中医研究, 2021, 34(7): 74-76. [9] 李佳鑫, 王冰, 于淼. 基于“扶正祛邪”治则的中药抗肿瘤作用机制的研究进展[J]. 中草药, 2021, 52(18): 5751-5757. [10] LI J J, ZHANG X, SHEN X C, et al. Phytochemistry and biological properties of isoprenoid flavonoids from Sophora flavescens Ait[J]. Fitoterapia, 2020, 143: 104556. doi: 10.1016/j.fitote.2020.104556 [11] MASUDA H, ZHANG D W, BARTHOLOMEUSZ C, et al. Role of epidermal growth factor receptor in breast cancer[J]. Breast Cancer Res Treat, 2012, 136(2): 331-345. doi: 10.1007/s10549-012-2289-9 [12] BURNESS M L, GRUSHKO T A, OLOPADE O I. Epidermal growth factor receptor in triple-negative and basal-like breast cancer: promising clinical target or only a marker?[J]. Cancer J, 2010, 16(1): 23-32. doi: 10.1097/PPO.0b013e3181d24fc1 [13] NAKAI K, HUNG M C, YAMAGUCHI H. A perspective on anti-EGFR therapies targeting triple-negative breast cancer[J]. Am J Cancer Res,2016,6(8):1609-1623. [14] LI W, HOU J Z, NIU J, et al. Akt1 inhibition promotes breast cancer metastasis through EGFR-mediated β-catenin nuclear accumulation[J]. Cell Commun Signal, 2018, 16(1): 82. doi: 10.1186/s12964-018-0295-1 [15] GAO F, ALWHAIBI A, SABBINENI H, et al. Suppression of Akt1-β-catenin pathway in advanced prostate cancer promotes TGFβ1-mediated epithelial to mesenchymal transition and metastasis[J]. Cancer Lett, 2017, 402: 177-189. doi: 10.1016/j.canlet.2017.05.028 [16] CHEN K, JIAO X M, ROCCO A D, et al. Endogenous cyclin D1 promotes the rate of onset and magnitude of mitogenic signaling via Akt1 Ser473 phosphorylation[J]. Cell Rep, 2020, 32(11): 108151. doi: 10.1016/j.celrep.2020.108151 [17] GEELEN C T, SAVAS P, TEO Z L, et al. Clinical implications of prospective genomic profiling of metastatic breast cancer patients[J]. Breast Cancer Res, 2020, 22(1): 1-13. doi: 10.1186/s13058-019-1178-0 [18] DOSSOU A S, BASU A. The emerging roles of mTORC1 in macromanaging autophagy[J]. Cancers, 2019, 11(10): 1422. doi: 10.3390/cancers11101422 [19] SCHMID P, ABRAHAM J, CHAN S, et al. Capivasertib plus paclitaxel versus placebo plus paclitaxel as first-line therapy for metastatic triple-negative breast cancer: the PAKT trial[J]. J Clin Oncol, 2020, 38(5): 423-433. doi: 10.1200/JCO.19.00368 [20] LI C W, XIA W Y, LIM S O, et al. AKT1 inhibits epithelial-to-mesenchymal transition in breast cancer through phosphorylation-dependent Twist1 degradation[J]. Cancer Res, 2016, 76(6): 1451-1462. doi: 10.1158/0008-5472.CAN-15-1941 [21] LIU H, RADISKY D C, NELSON C M, et al. Mechanism of Akt1 inhibition of breast cancer cell invasion reveals a protumorigenic role for TSC2[J]. Proc Natl Acad Sci U S A, 2006, 103(11): 4134-4139. doi: 10.1073/pnas.0511342103 [22] HANKER A B, KAKLAMANI V, ARTEAGA C L. Challenges for the clinical development of PI3K inhibitors: strategies to improve their impact in solid tumors[J]. Cancer Discov, 2019, 9(4): 482-491. doi: 10.1158/2159-8290.CD-18-1175 [23] VERNIERI C, CORTI F, NICHETTI F, et al. Everolimus versus alpelisib in advanced hormone receptor-positive HER2-negative breast cancer: targeting different nodes of the PI3K/AKT/mTORC1 pathway with different clinical implications[J]. Breast Cancer Res, 2020, 22(1): 33. doi: 10.1186/s13058-020-01271-0 [24] HYMAN D M, SMYTH L M, DONOGHUE M T A, et al. AKT inhibition in solid tumors with AKT1 mutations[J]. J Clin Oncol, 2017, 35(20): 2251-2259. doi: 10.1200/JCO.2017.73.0143 [25] KIM S B, DENT R, IM S A, et al. Ipatasertib plus paclitaxel versus placebo plus paclitaxel as first-line therapy for metastatic triple-negative breast cancer (LOTUS): a multicentre, randomised, double-blind, placebo-controlled, phase 2 trial[J]. Lancet Oncol, 2017, 18(10): 1360-1372. doi: 10.1016/S1470-2045(17)30450-3 [26] RAO G H, PIEROBON M, KIM I K, et al. Inhibition of AKT1 signaling promotes invasion and metastasis of non-small cell lung cancer cells with K-RAS or EGFR mutations[J]. Sci Rep, 2017, 7(1): 7066. doi: 10.1038/s41598-017-06128-9 [27] LEI J T, GOU X X, SEKER S, et al. ESR1 alterations and metastasis in estrogen receptor positive breast cancer[J]. J Cancer Metastasis Treat, 2019, 5: 38. [28] JESELSOHN R, BUCHWALTER G, ANGELIS C D, et al. ESR1 mutations—a mechanism for acquired endocrine resistance in breast cancer[J]. Nat Rev Clin Oncol, 2015, 12(10): 573-583. doi: 10.1038/nrclinonc.2015.117 [29] HICKEY T E, SELTH L A, CHIA K M, et al. The androgen receptor is a tumor suppressor in estrogen receptor-positive breast cancer[J]. Nat Med, 2021, 27(2): 310-320. doi: 10.1038/s41591-020-01168-7 [30] YU L H, YU Q H, XU C C, et al. The expression of androgen receptor in triple-negative breast cancer and the effect of a traditional Chinese medicine formula on disease-free survival[J]. Gland Surg, 2022, 11(11): 1772-1783. doi: 10.21037/gs-22-508 -
点击查看大图
图(10) / 表(2)
计量
- 文章访问数: 1371
- HTML全文浏览量: 360
- PDF下载量: 34
- 被引次数: 0
