基于网络药理学和分子对接探讨鸦胆子治疗结直肠癌的作用机制

结直肠癌（ colorectal cancer ， CRC ）是消化系统常见的恶性肿瘤之一，是全球癌症死亡率的第 2 大原因 [1] 。大多数 CRC 是散发的，遗传及环境因素是其重要的致病因素 [2] 。化疗是 CRC 的关键治疗方案之一 [3] 。传统中药制剂作为癌症化疗的辅助治疗已经在全世界范围内被广泛接受 [4] 。大量研究及临床实践向肿瘤的综合治疗不断拓展，已经证明传统中药作为辅助疗法，与化疗或者放疗相结合，可以提高疗效、减少不良反应 [5-7] 。

鸦胆子 Brucea javanica (L.) Merr. 是原生于我国东南部及其他热带、亚热带地区的苦木科灌木。以鸦胆子果实为原料生产的中药抗肿瘤注射液，已广泛在临床被用于肺癌、肺癌脑转移和胃肠道肿瘤的辅助治疗。临床研究证明鸦胆子可以增强临床上晚期非小细胞肺癌的疗效并降低化疗的不良反应 [8-9] 。鸦胆子苦醇是鸦胆子的重要活性成分 [10] 。研究发现，鸦胆子苦醇对乳腺癌、鼻咽癌、肺癌和其他恶性肿瘤均有抑制作用 [11-14] 。鸦胆子果实的乙醇提取物可能通过上调 p53 和抑制核因子 -κB （ nuclear factor-κB ， NF-κB ）抗人结肠癌 HCT116 细胞生长 [15] 。 然而，其潜在机制尚不清楚。

网络药理学是一种结合计算机科学与医学来阐述药物对疾病作用机制的新兴方法，能系统地研究、鉴定中药配方的生物活性化合物并可视化其多靶点、多途径的作用机制，已被广泛应用于预测、分析中药的药理作用及潜在机制 [4,16-17] 。因此，本研究从网络药理学角度筛选鸦胆子的活性成分及针对 CRC 的靶点和信号通路。分子对接是研究小分子药物与靶点蛋白受体之间的相互作用及亲和力的模拟计算方法。分子对接可以作为网络药理学的进一步验证，二者互补地用于中药研究 [18] 。本研究通过分子对接模拟鸦胆子已鉴定组分与其可能靶标蛋白的分子相互作用并计算结合力。鸦胆子苦醇作为被验证的鸦胆子抗肿瘤活性成分 [10] ，本研究进一步在裸鼠荷瘤模型中验证其抗 CRC 的潜力及网络药理学预测的效应机制。

1 材料

1.1 动物

12 只 SPF 级雌性 BALB/c-Nu 小鼠， 6 周龄，体质量 16 ～ 20 g ，购自北京华阜康生物科技股份有限公司，质量合格证编号 No.110322220103372262 ，质量许可证号 SCXK （京） 2019-0008 ，动物使用许可证号 SYXK （吉） 2018-0001 。动物于相对湿度 45% ～ 55% 、温度（ 22 ± 1 ）℃、 12 h 光暗循环的环境下，适应性饲养 1 周，自由进食饮水。动物实验经吉林大学基础学院伦理委员会批准（批准号 202073 ）。

1.2 细胞

HCT116 细胞购自武汉普诺赛生命科技有限公司。

1.3 药品与试剂

鸦胆子苦醇（批号 B26016 ，质量分数≥ 98% ）、 β- 环糊精（批号 T66291 ）购自上海源叶生物科技有限公司；二甲基亚砜（ dimethyl sulfoxide ， DMSO ，批号 D8418 ）购自美国 Sigma 公司； RIPA 裂解液、 BCA 蛋白浓度测定试剂盒（批号 P0006C ）购自上海碧云天生物技术有限公司； 4% ～ 15% 连续梯度聚丙烯酰胺预制胶（批号 DG101-01 ）购自北京全式金生物技术股份有限公司；表皮因子生长受体（ epidermal growth factor receptor ， EGFR ）抗体（批号 AF6043 ）、磷脂酰肌醇 3- 激酶（ phosphatidylinositol 3-kinase ， PI3K ）抗体（批号 AF6241 ）、蛋白激酶 B （ protein kinase B ， Akt ）抗体（批号 AF6261 ）、细胞周期蛋白依赖性激酶 4 （ cyclin-dependent kinase 4 ， CDK4 ）抗体（批号 DF6102 ）、细胞周期蛋白 B1 （ cyclin B1 ）抗体（批号 AF6168 ）、 cyclin D1 抗体（批号 AF0931 ）、剪切型半胱氨酸天冬氨酸蛋白酶 -3 （ cleaved cystein-asparate protease-3 ， cleaved Caspase-3 ）抗体（批号 AF7022 ）、 B 淋巴细胞瘤 -2 （ B-cell lymphoma-2 ， Bcl-2 ）抗体（批号 AF6139 ）、 Bcl-2 相关 X 蛋白（ Bcl-2 associated X protein ， Bax ）抗体（批号 AF0120 ）、蛋白基质金属蛋白酶 7 （ matrix metalloproteinase 7 ， MMP7 ）抗体（批号 AF0218 ）、 MMP2 抗体（批号 AF5330 ）、甘油醛 -3- 磷酸脱氢酶（ glyceraldehyde-3-phosphate dehydrogenase ， GAPDH ）抗体（批号 AF7021 ）、电化学发光试剂盒（批号 KF001 ）均购自江苏亲科生物研究中心有限公司； HRP 标记的羊抗兔二抗（批号 SA00001-2 ）购自武汉三鹰生物技术有限公司。

1.4 仪器

ChemiScope 6000 Touch 型凝胶成像仪（上海勤翔科学仪器有限公司）； SC-412 型 4 ℃立式冷藏柜（青岛市海尔股份有限公司）； PowerPac3000 型电泳仪（美国 Bio-Rad 公司）； MOV-112F 型细胞培养箱（日本三洋公司）； V18R 型高速冷冻离心机（瑞士 Dynamica 公司）； BETS-100 型摇床振荡器（海门市其林贝尔公司）； Hemo 3600V 型 Shinova 血液学分析仪（上海麦本医疗科技有限公司）。

2 方法

2.1 鸦胆子的网络药理学分析及分子对接

2.1.1 鸦胆子的活性成分及靶点 利用中医系统药理学数据库（ TCMSP ， https://old.tcmsp-e.com/tcmsp.php ）、中医药资料库（ http://tcm.cmu.edu.tw/ ）以及文献检索的方法收集鸦胆子的主要活性成分，将能在 PubChem 中检索对应的 CAS 号和 SMILES 号的成分进行收集。同时，对收集的成分进行药动学信息检索，以口服生物利用度（ oral bioavailability ， OB ） ≥ 30% 、类药性（ drug-likeness ， DL ）≥ 0.18 为筛选条件。

2.1.2 CRC 的靶点及交集靶点的获取 在 GeneCards （ https://www.genecards.org ）中，以“ CRC ”为检索词，物种限定为“ Homo sapiens ”，对其进行注释和预测。收集相关性值大于 1 的靶点作为疾病靶点。在 Venn 2.1.0 制图平台（ https://bioinfogp.cnb. csic.es/tools/venny/index.html ）绘制 CRC 与鸦胆子的交集靶点。

2.1.3 “活性成分 -CRC ”网络的构建 通过 Cytoscape 3.7.2 （ https://cytoscape.org ）可视化“活性成分 -CRC ”网络，连接的节点越密集，受该成分调控的靶点越多，表明该成分是鸦胆子抗 CRC 的重要活性成分。

2.1.4 蛋白质 - 蛋白质相互作用（ protein-protein interaction ， PPI ）网络的构建及关键靶点筛选 将 CRC 靶点导入 STRING 11.0 数据库（ https://cn.string-db.org ）进行分析，模式设置为 “ 多种蛋白 ” ，物种限定为“ Homo sapiens ”。预读数据后，其余蛋白的置信度≥ 0.99 ，与其余蛋白无联系的蛋白被隐藏。然后构建 PPI 网络，分数设置为 0.4 ，剔除不符合得分的靶标。利用 Barplot 绘图函数，对蛋白质连接节点的个数进行技术并绘制直方图，寻找核心靶点。

2.1.5 基因本体（ gene ontology ， GO ）功能及京都基因与基因组百科全书（ Kyoto encyclopedia of genes and genomes ， KEGG ）通路富集分析 利用 R 3.6.0 软件进行 GO 功能富集分析，从靶细胞的分子功能、细胞成分和生物过程确定基因富集的功能， P 值与富集程度呈正相关。用 R 将结果绘制成相应的直方图和气泡图，选取筛选校正后 P ＜ 0.05 的富集结果。利用 KEGG 途径富集分析获得枢纽靶点作用的途径。

2.1.6 分子对接 将小分子药物的 mol2 格式的文件导入 AutoDockTools （ Version 1.5.7 ），对小分子药物进行平衡电荷、将非极性氢原子与相对应的碳原子合并等修饰获得 3D 结构，然后转化为 PDBQT 格式的文件。从蛋白质数据库（ https://www.rcsb.org/ ， Protein Data Bank ， PDB ）获取 EGFR （ ID ： 1M17 ）、 Caspase-3 （ ID ： 4QUB ）、 Cyclin D1 （ ID ： 5VZU ）的晶体结构。通过 AutoDockTools 对受体蛋白进行去水、加氢、平衡电荷等预处理后转化为 PDBQT 格式的文件。利用 AutoDockTools 对小分子药物和配体蛋白进行分子对接，将受体和配体的 PDBQT 文件导入 AutoDockTools ，构建对接口袋。在大分子上设置中心，并设置 X 、 Y 和 Z 的参数以保证蛋白质完全被覆盖。将输出的最大结合模式数设置为 20 。分子对接后生成活性位点位置，并计算结合能和氢键数。结合能为负值说明小分子药物配体可以与受体靶点蛋白自发结合，结合能越小说明结合越稳定，结合能＜ −17.782 kJ/mol 即为能较好结合。用 OpenBabel 将对接好的文件转换为 PDB 格式，通过 PyMOL （ 3D ， Version 2.2.0 ）和 LigPlot （ 2D ）对偶联模型进行可视化，并展示氢键，这 2 种模型都被广泛应用于可视化对接结果，包括配体与蛋白质主链或侧链元件之间的对接位点和氢键相互作用模式。

2.2 鸦胆子苦醇治疗CRC的作用研究

2.2.1细胞培养 HCT116 细胞用含 10% 胎牛血清、 100 U/mL 青霉素和 100 μg/mL 链霉素的 DMEM 培养基，于 37 ℃ 、 5% CO2 的培养箱中培养。

2.2.2 动物模型的制备 BALB/c-Nu 小鼠右背侧 sc 100 μL HCT116 细胞（ 2 × 106 个 ），待肿瘤生长体积为 80 ～ 100 mm3 ， 小鼠随机分为对照组和鸦胆子苦醇（ 2 mg/kg ）组，每组 6 只。鸦胆子苦醇溶于 1% DMSO 和 20% β- 环糊精中，给药组 ip 药物，对照组 ip 等体积的 DMSO 和 β- 环糊精，每隔 1 d 给药 1 次。每隔 1 d 记录肿瘤大小和小鼠体质量，当肿瘤大小达到 2000 mm3 时，采用颈椎脱位法处死小鼠，收集肿瘤，测量肿瘤体积和质量，采集小鼠眶后静脉窦全血。

2.2.3 Western blotting 检测肿瘤组织细胞周期、侵袭、迁移、增殖及 EGFR/PI3K/Akt 通路相关蛋白表达 取各组小鼠肿瘤组织，加入含蛋白酶抑制剂和磷酸酶抑制剂的预冷 RIPA 裂解液提取蛋白， 12 000 r/min 离心 10 min ，取上清，采用 BCA 蛋白浓度测定试剂盒测定蛋白浓度。蛋白样品经十二烷基硫酸钠 - 聚丙烯酰胺凝胶电泳，转至 PVDF 膜，用 5% 脱脂牛奶室温封闭 30 min ，分别加入 EGFR （ 1 ∶ 1000 ）、 PI3K （ 1 ∶ 1000 ）、 Akt （ 1 ∶ 1000 ）、 CDK4 （ 1 ∶ 1000 ）、 cyclin B1 （ 1 ∶ 1000 ）、 cyclin D1 （ 1 ∶ 1000 ）、 cleaved Caspase-3 （ 1 ∶ 1000 ）、 Bax （ 1 ∶ 2000 ）、 Bcl-2 （ 1 ∶ 1000 ）、 MMP7 （ 1 ∶ 1500 ）、 MMP2 （ 1 ∶ 1000 ）、 GAPDH （ 1 ∶ 20 000 ）抗体， 4 ℃孵育过夜；加入二抗（ 1 ∶ 20 000 ），室温孵育 1 h 。使用电化学发光试剂盒显影，采用成像仪及 Image J 2.1 软件分析条带灰度值。

2.2.4 肝肾功及血常规指标检测 采用 Shinova 血液学分析仪分析小鼠的肝功能指数、肾功能指数和血常规指标，包括丙氨酸氨基转移酶（ alanine aminotransferase ， ALT ）、天冬氨酸氨基转移酶（ aspartate aminotransferase ， AST ）、 γ- 谷氨酰转肽酶（ γ-glutamyl transpeptidase ， γ-GT ）、总胆汁酸（ total bile acid ， T-BIL ）、白蛋白（ albumin ， ALB ）、血尿素氮（ blood urea nitrogen ， BUN ）、肌酐（ creatinine ， CREA ）、尿酸（ uric acid ， UA ）、白细胞（ white blood cell ， WBC ）、红细胞（ red blood cell ， RBC ）、血小板（ platelet ， PLT ）和血红蛋白（ haemoglobin ， HGB ）。

2.2.5 统计学分析 数据以 表示，采用 SPSS 26.0 软件进行单因素方差分析，使用 GraphPad Prism 8.0.1 软件可视化结果。

3 结果

3.1 网络药理学分析及分子对接

3.1.1 鸦胆子活性成分及对应靶点的获取 通过 TCMSP 检索共获得 67 个活性成分，以 OB ≥ 30% 、 DL ≥ 0.18 作为补充条件筛选，共获得 15 个鸦胆子主要活性成分（表 1 ）。收集活性成分对应的靶点，去除重复值后共获得 34 个靶点。

3.1.2 CRC 相关的靶点及与药物交集靶点的获取 Genecard 数据库中检索到人类 CRC 相关的靶点共 9410 个，经过用 Venn 2.1.0 得到 32 个鸦胆子与 CRC 的交集靶点（图 1 ），包括 EGFR 、 Caspase-3 、 cyclin D1 、 CDK4 、 MMP-9 和血管内皮生长因子 A （ vascular endothelial growth factor A ， VEGFA ）等。

3.1.3 “活性成分 -CRC ”网络的构建 通过 Cytoscape 3.7.2 可视化“活性成分 -CRC ”网络（图 2 ），包含 37 个节点和 97 个连接。椭圆节点表示靶点，方框节点表示活性成分，灰色线表示节点间相互作用。

3.1.4 PPI 网络分析 将 32 个鸦胆子与 CRC 交集靶点输入 STRING ，构建 PPI 网络（图 3 ），共获得 32 个节点和 152 条节点之间的连线。节点是鸦胆子与 CRC 交集靶点编码的蛋白，节点的连线越多，表示相互作用越多，靶点的作用越重要。根据节点的连线数量，通过 Barplot 计算并统计前 30 个靶点（图 4 ），前 3 名为 EGFR 、 Caspase-3 和 cyclin D1 。

3.1.5 GO 功能和 KEGG 通路富集分析 如图 5 所示，鸦胆子治疗 CRC 主要涉及细胞凋亡、自噬、基因表达负调控、蛋白磷酸化以及蛋白结合、转录因子结合等，主要作用于 PI3K/Akt 信号通路、肿瘤坏死因子（ tumor necrosis factor ， TNF ）信号通路和 p53 信号通路等。

3.1.6 分子对接 选取 PPI 网络中排名前 3 名的靶点 EGFR 、 Caspase-3 及 cyclin D1 ，分别与鸦胆子治疗 CRC 的 4 个关键活性成分进行分子对接（图 6 ）， LigPlot 将结果二维可视化，标明氢键的个数、长度及小分子药物配体连接的氨基酸残基（图 7 ）。一般 认为，小于 −17.782 、 −20.92 或 −29.288 kJ/mol 的结合能分别表明配体与受体之间有一定的、好的或较强的结合活性。结合能反映了受体和配体之间结合的可能性。结合能越低，受体与配体的亲和力越高，构象越稳定。其中，鸦胆子苦醇与 EGFR 的结合能最低（ −37.530 5 kJ/mol ）。与鸦胆子苦醇和木犀草素相比，鸦胆子苷 B 和 β- 谷甾醇与 3 个靶点蛋白的结合能均较高，亲和力较低。鸦胆子苦醇与 EGFR 、 Caspase-3 及 cyclin D1 的结合能分别为 −37.530 5 、 −33.137 3 、 −30.710 6 kJ/mol ，平均结合能为 −33.792 8 kJ/mol 。木犀草素与 EGFR 、 Caspase-3 及 cyclin D1 的结合能为 −26.652 1 、 −30.585 0 、 −25.689 8 kJ/mol ，平均结合能为 −27.642 3 kJ/mol 。

3.2 鸦胆子苦醇抗CRC的作用验证

3.2.1 鸦胆子苦醇对荷瘤裸鼠模型肿瘤生长的影响 如图 8 所示，与对照组比较，鸦胆子苦醇组裸鼠肿瘤生长速度被明显抑制（ P ＜ 0.01 ），证明了鸦胆子苦醇对 CRC 具有明显的抑制作用。

3.2.2 鸦胆子苦醇对荷瘤裸鼠模型肿瘤组织细胞周期、侵袭、迁移、增殖及EGFR/PI3K/Akt 通路相关蛋白表达的影响 如图 9 所示，与对照组比较，鸦胆子苦醇组肿瘤组织中 EGFR 、 PI3K 、 Akt 、 CDK4 、 cyclin D1 、 MMP9 、 MMP7 和 Bcl-2 蛋白表达水平均显著降低（ P ＜ 0.01 ）， cyclin B1 、 cleaved Caspase-3 、 Bax 蛋白表达水平均明显升高（ P ＜ 0.01 ）， Bax/Bcl-2 值显著升高（ P ＜ 0.01 ）。

3.2.3 鸦胆子苦醇对荷瘤裸鼠生物安全性 评价 为研究鸦胆子苦醇对荷瘤裸鼠代谢器官的影响，检测对照组和鸦胆子苦醇组荷瘤裸鼠的肝功能指标（ ALT 、 AST 、 ALB 、 γ-GT 、 T-BIL ）及肾功能指标（ BUN 、 CREA 、 UA ）和血常规指标（ WBC 、 RBC 、 PLT 、 HGB ），如图 10 所示，与对照组比较，鸦胆子苦醇组荷瘤裸鼠肝、肾功能和血常规指标均无显著差异。

4 讨论

本研究通过网络药理学分析鸦胆子的主要活性成分及其抗 CRC 的潜在机制，共获得鸦胆子苦醇、木犀草素、鸦胆子苷 B 、 β- 谷甾醇 4 个小分子药物活性成分和 EGFR 、 Caspase-3 、 cyclin D1 等 32 个核心靶点。 KEGG 通路富集分析表明鸦胆子抑制 CRC 可能与 PI3K/Akt 通路、 TNF 信号通路和 p53 信号通路等相关。在 4 个活性成分中，鸦胆子苦醇和木犀草素对应的靶点明显多于鸦胆子苷 B 和 β- 谷甾醇。将网络药理学预测的核心成分鸦胆子苦醇及木犀草素与关键靶点蛋白进行分子对接， 2 个核心成分与 EGFR 、 Caspase-3 和 cyclin D1 蛋白的活性结合位点的结合能均小于 −16.736 kJ/mol ，表明核心成分与关键的靶点蛋白之间可以自发结合且具有较强的结合活性。在筛选的 4 种鸦胆子活性成分中，鸦胆子苦醇和木犀草素均可与 EGFR 相互作用， EGFR 是多种靶向药物的治疗靶点，如吉非替尼和厄洛替尼用于治疗非小细胞肺癌，拉帕替尼用于治疗晚期或转移性乳腺癌 [19] 。本研究发现，木犀草素与 EGFR 的 Lys721 、 Glu738 、 Ala731 形成 3 个氢键，鸦胆子苦醇与 EGFR 的 Ile917 形成 1 个氢键。木犀草素和鸦胆子苦醇对 EGFR 的亲和力表明木犀草素 -EGFR 、鸦胆子苦醇 -EGFR 相互作用的有效性。

鸦胆子苦醇是鸦胆子的主要活性成分。已有研究证实，鸦胆子苦醇可以阻断 PI3K/Akt 通路的信号传导，抑制胃癌等多种癌症的发生发展 [20] 。网络药理学分析结果表明，鸦胆子苦醇可通过 EGFR/PI3K/ Akt 通路抗 CRC 。进一步通过 Western blotting 实验验证了鸦胆子苦醇可以使该级联通路明显失活从而达到抑制 CRC ，验证了网络药理学预测结果。本研究通过构建裸鼠移植瘤模型，研究了鸦胆子苦醇对 CRC 的治疗作用。鸦胆子苦醇明显抑制荷瘤裸鼠模型体内肿瘤的生长，同时与代谢相关的肝脏功能和肾脏功能未受到明显的影响，血常规指标也没有观察到显著差异。表明鸦胆子苦醇具有很好的抑制 CRC 肿瘤细胞生长的作用。细胞周期与肿瘤的发展密切相关，主要由 G1/G0 和 G2/M 期组成，由 cyclin/CDK 复合物和 CDK 抑制剂协调。研究表明，鸦胆子苦醇通过作用于 CDK4 和 cyclin D1 调节非小细胞肺癌的 G1/G0 期， cyclin D1 表达的降低导致 G1/G0 期细胞的阻滞 [21] ，与本研究的结果一致，表明鸦胆子苦醇能够诱导肿瘤细胞的 G1/G0 期阻滞，从而抑制肿瘤细胞增殖。诱导细胞凋亡是抑制肿瘤发展的关键环节，受 Caspase 和 Bcl-2 家族蛋白的调控， Bcl-2 家族蛋白包括促凋亡效应蛋白如 Bax 和抗凋亡效应蛋白如 Bcl-2 。促凋亡效应蛋白在凋亡应激源的作用下导致线粒体外膜通透化，释放细胞色素 C 来切割 Caspase-3 ，从而导致不可逆的细胞凋亡。鸦胆子苦醇可通过上调 cleaved Caspase-3 表达和 Bax/Bcl-2 值，促进多种肿瘤细胞凋亡 [22] ，与本研究结果一致，表明鸦胆子苦醇可诱导肿瘤细胞凋亡，可能涉及线粒体凋亡途径。肿瘤转移是肿瘤患者死亡的相关因素。细胞外基质被 MMP 修饰和降解，从而导致肿瘤细胞的分离和迁移。本研究发现，鸦胆子苦醇能够降低肿瘤组织中 MMP7 和 MMP9 的表达，具有抑制肿瘤细胞侵袭和迁移的能力，与鸦胆子苦醇对胃癌 SGC-7901 细胞的作用一致 [19] 。

药物的积累和清除可能损害器官。因此，通过对小鼠的血液生化指标来评估鸦胆子苦醇是否具有潜在的毒性作用。结果显示，肝、肾功能和血常规指标未见药物引起的显著性差异。因此，鸦胆子苦醇是一种安全、有效、有前景的抗肿瘤佐剂。然而，本研究有几个局限性： ① 来自在线数据库的信息是基于审查和预测的数据，因此，未经证实和未记录的化合物或靶点可能未被列入本研究中； ② 目前对鸦胆子 15 种化合物的定量测定研究尚不完全，因此，未来应进行内容确定的研究； ③ 鸦胆子苦醇虽然为鸦胆子抗 CRC 的最重要的生物活性成分，但不能完全代表鸦胆子，因此，需要进一步的研究来探索鸦胆子体内外治疗 CRC 的潜在分子机制。

本研究先通过网络药理学阐明了鸦胆子抗 CRC 最重要的 4 个生物活性成分（鸦胆子苦醇、木犀草素、鸦胆子苷 B 、 β- 谷甾醇）和 3 个重要靶点（ EGFR 、 Caspase-3 、 cyclin D1 ），主要作用于 EGFR/ PI3K/Akt 信号通路；通过分子对接技术验证了活性成分与靶点的结合能力，发现鸦胆子苦醇是鸦胆子抗 CRC 的关键成分；体内实验进一步证明了鸦胆子苦醇抗 CRC 的作用及分子机制，为鸦胆子治疗 CRC 提供了实验依据，并为从中药中探究活性成分、核心靶点和潜在机制提供了方法。

利益冲突所有作者均声明不存在利益冲突

参考文献（略）

来 源：赵韦欣，王 晴，王梦齐，季安璇，郑晶莹，赵淑华． 基于网络药理学和分子对接探讨鸦胆子治疗结直肠癌的作用机制[J]. 中草药, 2023, 54(6):1850-1859.

【福利时刻】购买细胞及相关实验，扫码或加微信（1278317307）入群。

本公众号已与浙江慧通测评动物实验中心搭建紧密的合作关系：该动物中心已获得动物使用许可并在建动物行为学实验中心，如您需要含药血清制备-成分分析、灌胃给药、造模、药理、毒理以及相关实验技术支持，欢迎联系我们（1278317307）！

【福利时刻】毕业季找工作点击链接，多家企业30W+年薪工作招聘任你选（赶紧点击进入注册，上传简历吧~ ） 。