2.6 共同差异表达蛋白统计统计模型组/假手术组差异蛋白在骨疏康组的表达情况。模型组表达上调，而在骨疏康组表达下调的蛋白有25 个；模型组表达下调，而在骨疏康组表达上调的蛋白有2 个，即模型组共有27 个异常表达的蛋白在骨疏康组得到纠正，见表2。

2.7 共同差异表达蛋白PPI 分析为获得蛋白质的更多功能信息，通过STRING 蛋白相互作用网络分析工具分析，发现共同差异蛋白中Gstz1[10]、Cxxc1[11]、Mecp2、Aldh3b1、Tspan14[12]、Ryr1、Myoz1 位于功能网络节点，其余蛋白均在网络边缘，见图7。图中Gstz1、Cxxc1、Mecp2 之间存在相互作用，Aldh3b1 与Tspan14、Ryr1 与Myoz1 之间存在相互作用。

3 讨论 Discussion

3.1 既往他人的研究贡献及存在不足随着骨质疏松症发病率逐年攀升，成为世界性健康问题和社会公共问题[13]。探索骨疏康等中医药治疗骨质疏松症的作用机制，可为临床应用提供更多实验依据，为骨质疏松症的中医药治疗及预防提供更大的优势和理论补充。临床研究表明，骨疏康治疗原发性骨质疏松症效果显著[14-15]，从调整全身功能入手防治骨质疏松症为其特点。药理学研究显示骨疏康对骨吸收、骨形成、修复骨微结构等方面具有一定作用[16]：WANG 等[17]实验表明骨疏康可抑制去卵巢小鼠破骨细胞的生成和刺激成骨细胞的生成；CHAI 等[18]研究发现骨疏康干预通过上调相关成骨因子而显著增强BMP-2/Smads 信号通路，且改善骨质疏松大鼠的骨微结构；LI 等[19]用不同剂量的骨疏康处理卵巢切除术小鼠，能显著提高维生素D 和钙的代谢水平。随着蛋白组学的迅速发展，成为近年来的研究热点，部分文献报道骨疏康防治骨质疏松症的具体作用机制尚未阐释其相关作用靶点，其分子机制尚未明确。

3.2 研究的意义利用新兴的蛋白组学技术和生物信息学分析，从蛋白质组学角度挖掘骨疏康治疗骨质疏松症的新靶点和信号通路，为临床应用提供实验基础。

此次研究中骨疏康组大鼠骨密度显著提高。从差异表达蛋白GO 功能富集分析来看，骨疏康组/模型组差异蛋白主要参与免疫球蛋白受体结合、血管内皮生长因子受体结合等生物过程。最新研究表明免疫球蛋白家族23(IGSF23)介导破骨细胞分化，抑制其表达可逆转去卵巢小鼠的骨流失[20]；并且研究也显示血管内皮生长因子与血管内皮生长因子受体结合诱导血管生成，而血管侵入参与骨重建过程[21-22]。从差异表达蛋白KEGG 通路分析来看，骨疏康组/模型组差异蛋白主要富集于酪氨酸代谢、乙醇代谢、JAK-STAT 信号通路等，同时与模型组/假手术组差异蛋白共同富集于酪氨酸代谢信号通路。已有研究表明将四环素与一个酪氨酸结合形成一个新化合物抑制了去卵巢后所致骨小梁的丢失，改善骨小梁空间结构，增加骨小梁板的厚度，其机制与抑制破骨细胞介导的骨吸收有关[23]。从共同差异表达蛋白来看，模型组Zc3hc1[24-25]、Svs3b、Hbe2、Mecp2、Myoz1 等27 个 异 常 表达蛋白在骨疏康组得以纠正，提示了骨疏康的治疗作用。结合骨质疏松症病理生理特点以及GO 功能注释、KEGG 通路富集分析、PPI 互作分析结果筛选出与骨代谢相关的Mecp2、Aldh3b1、Ryr1、Myoz14 个共同差异表达蛋白和重要通路进行分析。

表2 ｜模型组与骨疏康组的共同差异表达蛋白统计Table 2 ｜Statistics of co-differentially expressed proteins betweenGushukanggroup and model group表注：*为模型组表达下调的差异蛋白序号 蛋白名称 中文名称 基因名称 模型组vs. 假手术组 骨疏康组vs. 模型组P 值 差异倍数 P 值 差异倍数1 RCG - Svs3b 0.000 6.329 0.000 0.507 2 Epsilon 2 globin ε2 珠蛋白 Hbe2 0.000 5.373 0.001 0.526 3 Zc3hc1 protein Zc3hc1 蛋白 Zc3hc1 0.010 2.493 0.009 0.568 4 Glyceraldehyde-3-phosphate dehydrogenase 甘油醛-3-磷酸脱氢酶 - 0.002 2.071 0.021 0.760 5 Signal-induced proliferation-associated 1 信号诱导增殖相关1 Sipa1 0.046 1.883 0.048 0.606 6 Soluble calcium-activated nucleotidase 1 可溶性钙激活核苷酸酶1 Cant1 0.005 1.730 0.006 0.757 7 Uncharacterized protein 未鉴定蛋白 - 0.005 1.727 0.038 0.792 8 Maleylacetoacetate isomerase 顺丁烯二酰乙酰乙酸酯异构酶 Gstz1 0.016 1.600 0.021 0.832 9 Lipocalin 7, isoform CRA_a 脂蛋白7，CRA_a 亚型 Tinagl1 0.004 1.595 0.024 0.739 10 Protein unc-45 homolog A 蛋白质unc-45 同源物A Unc45a 0.004 1.539 0.014 0.683 11 Ubiquitin-associated and SH3 domain-containing, A 泛素相关和SH3 结构域，A Ubash3a 0.006 1.538 0.026 0.806 12 Churchill domain-containing 1 包含丘吉尔域1 Churc1 0.014 1.533 0.007 0.611 13 CCR4-NOT transcription complex, subunit 8 CCR4-NOT 转录复合物，亚基8 Cnot8 0.019 1.504 0.011 0.771 14 Aldehyde dehydrogenase family 3 member B1 乙醛脱氢酶家族3 成员B1 Aldh3b1 0.002 1.498 0.003 0.766 15 Leucine-rich repeat flightless-interacting protein 2 富含亮氨酸重复序列的非飞行相互作用蛋白2 Lrrfip2 0.002 1.486 0.001 0.633 16 Coiled-coil domain-containing 115 含线圈域115 Ccdc115 0.001 1.398 0.005 0.765 17 Ryanodine receptor 1 兰尼碱受体1 Ryr1 0.043 1.317 0.003 0.761 18 CXXC finger 1 (PHD domain) CXXC Finger 1(PHD 域) Cxxc1 0.019 1.317 0.025 0.784 19 Tetraspanin 川芎嗪 Tspan14 0.011 1.293 0.006 0.817 20 Lipin 1 脂蛋白1 Lpin1 0.027 1.286 0.036 0.826 21 RNA-binding motif protein 27 RNA 结合基序蛋白27 Rbm27 0.035 1.265 0.015 0.771 22 Ribosome assembly factor mrt4 核糖体装配因子MRT4 Mrto4 0.003 1.259 0.004 0.804 23 BH3-interacting domain death agonist BH3 相互作用域死亡激动剂 Bid 0.003 1.259 0.023 0.807 24 Phosducin-like protein 3 磷素样蛋白3 Pdcl3 0.029 1.257 0.007 0.713 25 Amino acid transporter 氨基酸转运蛋白 Slc1a3 0.009 1.205 0.002 0.726 26 *Methyl-CpG-binding protein 2 甲基-CpG 结合蛋白2 Mecp2 0.040 0.807 0.010 1.227 27 *Myozenin 1 肌红蛋白1 Myoz1 0.023 0.802 0.041 1.383

文章来源：《中国中医药图书情报杂志》 网址: http://www.zgzyytsqbzz.cn/qikandaodu/2021/0613/734.html