1 材料与方法
1.1 橙皮苷潜在作用靶点预测及收集
1.2 橙皮苷治疗肠道黏膜损伤潜在靶点收集
1.3 药物-疾病交集靶点蛋白质-蛋白质相互作用(protein-protein interaction,PPI)网络图的构建及核心靶点筛选
1.4 GO功能和KEGG通路富集分析
1.5 橙皮苷与核心靶点分子对接验证
1.6 橙皮苷与核心靶点分子动力学模拟分析
1.7 体外细胞试验验证
1.7.1 细胞培养
1.7.2 细胞增殖检测
1.7.3 靶基因mRNA相对表达水平测定
表1 引物序列Table 1 Primer sequences |
| 基因 Genes | 引物序列 Primer sequences (5'—3') | 登录号 Accession number | 产物长度 Product length/bp |
|---|---|---|---|
| β-肌动蛋白 β-actin | F: AGTTGAAGGTGGTCTCGTGG R: TGCGGGACATCAAGGAGAAG | NM_001444420.1 | 216 |
| 表皮生长因子受体 EGFR | F: GGACAACCCAGACTACCAGC R: GGCTGTCGGGACTATGTAGC | NM_214007.1 | 217 |
| 纤维连接蛋白1 FN1 | F: GAAGTGACGGACGTGAAGGT R: GGTCACCATGACGGTTGAGT | XM_003133642.5 | 300 |
| B淋巴细胞瘤2 BCL2 | F: GGTGAACTGGGGGAGGATTG R: GTGCCGGTTCAGGTACTCAG | XM_021099593.1 | 130 |
| 前列腺素内过氧化物合酶2 PTGS2 | F: AGACAGCATAAACTGCGCCT R: CATCATCAGACCAGGCACCA | NM_214321.1 | 183 |
| 胱天蛋白酶3 CASP3 | F: TTGAGACGGACAGTGGGACT R: CCGTCCTTTGAATTTCGCCAG | NM_214131.1 | 121 |
| 雌激素受体1 ESR1 | F: GGCCAAGCCCTCTCTTGATT R: AGGATCTCTAGCCAGGCACA | NM_214220.1 | 286 |
1.8 统计分析
2 结果与分析
2.1 橙皮苷和断奶仔猪肠道黏膜损伤靶点预测及筛选
图1 橙皮苷缓解断奶仔猪肠道黏膜损伤的靶点可视化筛选A:橙皮苷靶点基因并集韦恩图 Venn diagram of hesperidin target gene concatenation;B:断奶仔猪肠道黏膜损伤靶点基因并集韦恩图 Venn diagram of intestinal mucosal damage in weaned piglets target gene concatenation;C:橙皮苷与断奶仔猪肠道黏膜损伤共同靶点韦恩图 Venn diagram of common target of hesperidin and intestinal mucosal damage in weaned piglets。 Fig.1 Visual screen of targets for hesperidin alleviating intestinal mucosal damage in weaned piglets |
2.2 橙皮苷-断奶仔猪肠道黏膜损伤交集靶点PPI网络
图2 橙皮苷缓解断奶仔猪肠道黏膜损伤靶点PPI网络图仅注释核心靶点 only annotated core targets。EGFR:表皮生长因子受体 epidermal growth factor receptor;FN1:纤维连接蛋白1 fibronectin;BCL2:B淋巴细胞瘤2 B-cell lymphoma 2;PTGS2:前列腺素内过氧化物合酶2 prostaglandin G/H synthase 2;CASP3:胱天蛋白酶3 caspase 3;ESR1:雌激素受体1 estrogen receptor 1。图4同the same as Fig.4。 Fig.2 PPI network diagram of targets for hesperidin alleviating intestinal mucosal damage in weaned piglets |
表2 橙皮苷缓解断奶仔猪肠道黏膜损伤的核心靶点及其拓扑学特征Table 2 Core targets for hesperidin alleviating intestinal mucosal damage in weaned piglets and its topological characteristics |
| 核心靶点 Core targets | Uniprot编号 Uniprot ID | 中间中心度 Betweenness | 接近中心度 Closeness | 度值 Degree |
|---|---|---|---|---|
| 表皮生长因子受体 EGFR | P00533 | 0.165 6 | 0.830 2 | 35 |
| 纤维连接蛋白1 FN1 | P02751 | 0.145 5 | 0.771 9 | 31 |
| B淋巴细胞瘤2 BCL2 | P10415 | 0.077 3 | 0.771 9 | 31 |
| 前列腺素内过氧化物合酶2 PTGS2 | P35354 | 0.094 3 | 0.758 6 | 30 |
| 胱天蛋白酶3 CASP3 | P42574 | 0.054 6 | 0.745 8 | 29 |
| 雌激素受体1 ESR1 | P03372 | 0.052 1 | 0.721 3 | 27 |
2.3 GO功能和KEGG通路富集分析
图3 橙皮苷缓解断奶仔猪肠道黏膜损伤靶点GO条目柱形图仅注释主要的GO功能 only annotated main GO function。Count:计数;Biological processes:生物过程;Cellular components:细胞组分;Molecular functions:分子功能;proteolysis:蛋白水解;collagen catabolic process:胶原分解代谢过程;apoptotic process:细胞凋亡过程;cell-cell adhesion:细胞间黏附;extracellular matrix organization:细胞外基质组织;cellular response to hypoxia:缺氧的细胞反应;positive regulation of cell migration:细胞迁移的正向调控;kidney development:肾脏发育;positive regulation of apoptotic process:凋亡过程的正向调控;regulation of cell population proliferation:细胞群体增殖调控;negative regulation of neuron apoptotic process:神经元凋亡过程的负调控;extracellular matrix:细胞外基质;early endosome:早期内体;extracellular space:细胞外间隙;mitochondrion:线粒体;plasma membrane:质膜;cytoplasm:细胞质;endoplasmic reticulum:内质网;metalloendopeptidase activity:金属内肽酶活性;zinc ion binding:锌离子结合 protein homodimerization activity:蛋白质同源二聚体活性;integrin binding:整合素结合。 Fig.3 Bar chart of targets GO entries for hesperidin alleviating intestinal mucosal damage in weaned piglets |
图4 橙皮苷缓解肠道黏膜损伤靶点的KEGG通路富集分析结果AGE-RAGE signaling pathway in diabetic complications:糖尿病并发症晚期糖基化产物-晚期糖基化终末产物受体信号通路;Pathways in cancer:癌症通路;Endocrine resistance:内分泌抵抗;Kaposi sarcoma-associated herpesvirus infection:Kaposi肉瘤相关疱疹病毒感染;EGFR tyrosine kinase inhibitor resistance: 表皮生长因子受体酪氨酸激酶抑制剂耐药性;Human papillomavirus infection:人乳头瘤病毒感染;Apoptosis:细胞凋亡;MicroRNAs in cancer:癌症微小RNA;Colorectal cancer:结直肠癌;Chemical carcinogenesis-receptor activation:化学致癌-受体激活;Gastric cancer:胃癌;Lipid and atherosclerosis:脂质和动脉粥样硬化;C-type lectin receptor signaling pathway:C型凝集素受体信号通路;Hepatitis B:乙型肝炎;TNF signaling pathway:肿瘤坏死因子信号通路;Prolactin signaling pathway:催乳素信号通路;Apoptosis-multiple species:多物种细胞凋亡;Pancreatic cancer:胰腺癌;Proteoglycans in cancer:癌症蛋白多糖;Epstein-Barr virus infection:EB病毒感染。 Fig.4 Results of targets KEGG pathway enrichment analysis for hesperidin alleviating intestinal mucosal damage in weaned piglets |
2.4 靶蛋白与橙皮苷分子对接分析
表3 橙皮苷的基本信息和化学结构Table 3 Basic information and chemical structure of hesperidin |
| 名称 Name | PubChem编号 PubChem ID | 分子式 Molecular formula | 分子质量 Molecular weight/u | SMILES分子表达式 SMILES molecular expression formula | 分子结构 Molecular structure |
|---|---|---|---|---|---|
| 橙皮苷 Hesperidin | 10621 | C28H34O15 | 610.56 | OC1=C(OC)C=CC([C@@H]2CC(C3= C(C=C(O[C@@H]4O[C@H](CO [C@@H]5O[C@@H](C)[C@H](O) [C@@H](O)[C@H]5O)[C@@H](O) [C@H](O)[C@H]4O)C=C3O)O2)= O)=C1 | ![]() |
图5 橙皮苷与核心靶蛋白的分子对接可视化图谱Interactions:相互作用;van der Waals:范德华力;Conventional Hydrogen Bond:经典氢键;Carbon Hydrogen Bond:碳氢键;Unfavorable Donor-Donor:不利氢键供体;Pi-Anion:阴离子-π相互作用;Pi-Pi Stacked:π-π堆积;Pi-Pi T-shaped:π-π T型堆积;Pi-Akyl:π-烷基;MET:蛋氨酸 methionine;GLU:谷氨酸 glutamate;ASP:天冬氨酸 aspartate;LYS:赖氨酸 lysine;PRO:脯氨酸 proline;THR:苏氨酸 threonine;SER:丝氨酸 serine;ILE:异亮氨酸 isoleucine;VAL:缬氨酸 valine;ARG:精氨酸 arginine;GLN:谷氨酰胺;GLY:甘氨酸 Glycine;PHE:苯丙氨酸 phenylalanine;ASN:天冬酰胺 asparagine;TRP:色氨酸 tryptophan。 A:橙皮苷-EGFR复合物 hesperidin-EGFR complex;B:橙皮苷-FN1复合物 hesperidin-FN1 complex;C:橙皮苷-BCL2复合物 hesperidin-BCL2 complex;D:橙皮苷-PTGS2复合物 hesperidin-PTGS2 complex;E:橙皮苷-CASP3复合物 hesperidin-CASP3 complex;F:橙皮苷-ESR1复合物 hesperidin-ESR1 complex。 Fig.5 Molecular docking visualization of hesperidin and core target proteins |
表4 橙皮苷与核心靶蛋白的分子对接结果Table 4 Molecular docking results of hesperidin and core target proteins |
| 项目 Items | PDB编号 PDB ID | 结合能 Binding energy/ (kcal/mol) | 作用位点 Combination site |
|---|---|---|---|
| 表皮生长因子受体 EGFR | 8S1M | -5.10 | MET-987、GLU-985、ASP-984、GLN-976、ASP-974、GLY-810、 ARG-973、PRO-975、GLN-812、TYR-813、 ASN808、HIS-805、TYR-801 |
| 纤维连接蛋白1 FN1 | 2CK2 | -4.15 | LYS-86、PRO-87、THR-71、SER-89、ILE-8、VAL-1、ARG-6、 ALA-88、ASN-91、ASP-7、ILE-90、THR-69、CAL-4、SER-2 |
| B淋巴细胞瘤2 BCL2 | 4IEH | -5.63 | ARG-65、ARG-66、ASP-62、GLN-58、GLY-104、VAL-107、TRP-103、 PHE-157、THR-55、TYR-161、ALA-59、GLY-162、PRO-163 |
| 前列腺素内过氧化物合酶2 PTGS2 | 5F1A | -3.70 | GLY-167、PHE-52、ASN-34、VAL-165、PRO-160、LYS-166、 ASP-497、ASP-53、PRO-35、CYS-36、CYS-37、SER-38 |
| 半胱天冬酶3 CASP3 | 4JQY | -2.00 | GLU-248、ASN-208、ASP-990、PHE-252、PHE-247、GLU-246、 TRP-214、ASP-253、PHE-250 |
| 雌激素受体1 ESR1 | 3CBM | -2.60 | AGR-394、GLU-443、GLN-441、PHE-445、TRP-393、 GLU-323、GLY-442、LEU-320、LYS-449、GLY-390、 PRO-324、ILE-326、VAL-446 |
MET:蛋氨酸 methionine;GLU:谷氨酸 glutamate;ASP:天冬氨酸 aspartate;LYS:赖氨酸 lysine;PRO:脯氨酸 proline;THR:苏氨酸 threonine;SER:丝氨酸 serine;ILE:异亮氨酸 isoleucine;VAL:缬氨酸 valine;ARG:精氨酸 arginine;GLN:谷氨酰胺;GLY:甘氨酸 Glycine;PHE:苯丙氨酸 phenylalanine;ASN:天冬酰胺 asparagine;TRP:色氨酸 tryptophan。 |
2.5 靶蛋白-橙皮苷复合物分子动力学模拟结果
图6 BCL2与橙皮苷相互作用的分子动力学模拟分析RMSD:根均方偏差 root mean square deviation;RMSF:均方根波动 root mean square fluctuation;RG:回旋半径 radiuses of gyration;Hydrogen Bonds:氢键;Area:面积;Relative Free Energy:相对自由能。Time:时间。 A:受体-配体复合物0~100 ns的RMSD RMSD of the receptor-ligand complex from 0 to 100 ns;B:受体-配体复合物0~100 ns的RMSF RMSF of the receptor-ligand complex from 0 to 100 ns;C.回旋半径 radius of gyration;D:氢键数量 number of hydrogen bonds;E:溶剂可及表面积分析 solvent accessible surface area analysis;F:自由能景观图 Free energy landscape。图8同 the same as Fig.8。 Fig.6 Molecular dynamics simulation analysis of interaction between BCL2 and hesperidin |
2.6 橙皮苷对氧化应激IPEC-J2细胞增殖的影响
图8 橙皮苷对IPEC-J2细胞活力的影响Blank:空白组 Blank group;HSD-0:0 μmol/L橙皮苷组 0 μmol/L hesperidin group;HSD-100:100 μmol/L橙皮苷组 100 μmol/L hesperidin group;HSD-200:200 μmol/L橙皮苷组 200 μmol/L hesperidin group;HSD-400:400 μmol/L橙皮苷组 400 μmol/L hesperidin group;HSD-800:800 μmol/L橙皮苷组 800 μmol/L hesperidin group。 数据柱标相同小写字母表示差异不显著(P>0.05),不同小写字母表示差异显著(P<0.05)。Value columns with the same small letter mean no significant difference (P>0.05), while with different small letters mean significant difference (P<0.05). Fig.8 Effects of hesperidin on viability of IPEC-J2 cells |
2.7 橙皮苷对氧化应激IPEC-J2细胞核心靶基因表达的影响
表5 橙皮苷对氧化应激IPEC-J2细胞核心靶基因mRNA表达的影响Table 5 Effects of hesperidin on mRNA expression of core target genes in oxidative stress-induced IPEC-J2 cells |
| 基因 Genes | 组别 Groups | P值 P-value | ||
|---|---|---|---|---|
| 空白 Blank | HSD-0 | HSD-400 | ||
| 表皮生长因子受体 EGFR | 1.02±0.25a | 0.41±0.11b | 0.83±0.13a | 0.019 |
| 纤维连接蛋白1 FN1 | 1.01±0.19a | 0.42±0.17b | 0.76±0.10a | 0.021 |
| B淋巴细胞瘤2 BCL2 | 1.01±0.16a | 0.48±0.14b | 0.79±0.11a | 0.033 |
| 前列腺素内过氧化物合酶2 PTGS2 | 1.00±0.12b | 1.66±0.25a | 1.27±0.15b | 0.027 |
| 胱天蛋白酶3 CASP3 | 1.01±0.18b | 1.54±0.28a | 1.25±0.13a | 0.039 |
| 雌激素受体1 ESR1 | 1.01±0.20 | 0.80±0.17 | 0.85±0.15 | 0.353 |
同行数据肩标不同小写字母表示差异显著(P<0.05),相同或无字母表示差异不显著(P>0.05)。 | |
In the same row, values with different small letter superscripts mean significant difference (P<0.05), while with the same or no letter superscripts mean no significant difference (P>0.05). |

