劉文靜 楊娟 黃美雅,3 張楚天 孫攀 黃云梅,3*

1.福建中醫(yī)藥大學(xué)藥學(xué)院，福建 福州 350122 2.福建中醫(yī)藥大學(xué)中西醫(yī)結(jié)合研究院，福建 福州350122 3.福建省中西醫(yī)結(jié)合老年性疾病重點(diǎn)實(shí)驗(yàn)室，福建 福州 350122

絕經(jīng)后骨質(zhì)疏松癥(postmenopausal osteoporosis，PMOP)是因機(jī)體雌激素水平下降，骨質(zhì)流失增加，骨微結(jié)構(gòu)破壞，易致骨折的一種代謝性疾病[1]。中醫(yī)理論認(rèn)為絕經(jīng)后骨質(zhì)疏松癥的主要發(fā)病基礎(chǔ)是機(jī)體腎虧脾虛，故臨床上以補(bǔ)腎健脾、強(qiáng)筋壯骨為用藥原則。健骨顆粒是臨床常用經(jīng)驗(yàn)方，由煅狗骨、淫羊藿、山茱萸、黨參、山藥等中藥組成，具有補(bǔ)腎健脾、強(qiáng)筋壯骨之功效。前期實(shí)驗(yàn)[2-5]表明健骨顆粒可通過降低miR-14及CylinD1等的表達(dá)促進(jìn)成骨細(xì)胞分化，通過抑制RANKL等的表達(dá)從而抑制破骨細(xì)胞形成，發(fā)揮治療絕經(jīng)后骨質(zhì)疏松癥的作用。

代謝組學(xué)主要研究疾病、藥物效應(yīng)或其他刺激反應(yīng)所引起機(jī)體代謝物波動(dòng)的詳細(xì)情況[6]。中藥方劑發(fā)揮作用通常是方劑中有效成分網(wǎng)絡(luò)對機(jī)體疾病狀態(tài)下生物分子網(wǎng)絡(luò)調(diào)節(jié)的綜合結(jié)果[7]，適合以“組學(xué)”的手段進(jìn)一步開展系統(tǒng)性的研究。本文在前期健骨顆粒實(shí)驗(yàn)基礎(chǔ)上通過對大鼠血清代謝組學(xué)研究，探討健骨顆粒對去卵巢骨質(zhì)疏松大鼠的作用機(jī)理，以期為臨床防治絕經(jīng)后骨質(zhì)疏松癥提供新的實(shí)驗(yàn)依據(jù)。

1 材料與方法

1.1 材料

1.1.1藥材：健骨顆粒組方為煅狗骨、淫羊藿、山茱萸、山藥、黨參等。原藥材購于福建省醫(yī)藥公司，由福建中醫(yī)藥研究院中試車間負(fù)責(zé)加工制備，每克健骨顆粒含原生藥2.9 g。

1.1.2儀器：Discovery Wi雙能X線骨密度儀(美國 Hologic公司)；IG-A1000 N萬能材料試驗(yàn)機(jī)(日本島津儀器公司)；Q Exact質(zhì)譜儀(美國Thermo Scientific公司)；UltiMate 3000高效液相系統(tǒng)(美國Thermo Scientific公司 )；HYPERSIL GOLD AQ色譜柱(100 mm × 2.1 mm，1.9 μm，美國Thermo Scientific公司)。

1.1.3動(dòng)物：60只3月齡SPF級雌性SD大鼠，質(zhì)量(280±20)g，許可證號(hào)：SCXK(浙)2014-0001，購于浙江省實(shí)驗(yàn)動(dòng)物中心，飼養(yǎng)于福建中醫(yī)藥大學(xué)實(shí)驗(yàn)動(dòng)物中心鼠類實(shí)驗(yàn)室，環(huán)境設(shè)施為SPF級，合格證號(hào)：SYXK(閩)2014-0005。在溫度 21 ℃、相對濕度 70 %～80 %、通風(fēng)良好、自然光照的環(huán)境下，給予大鼠標(biāo)準(zhǔn)飼料飼養(yǎng)，自由進(jìn)食、飲水、活動(dòng)。

1.2 方法

1.2.1動(dòng)物實(shí)驗(yàn)及樣本采集：將60只雌性SD大鼠隨機(jī)分為正常組、模型組和健骨顆粒組，每組20只。模型組和健骨顆粒組經(jīng)腹腔注射0.3 mL/100 g 的10%水合氯醛麻醉，摘除雙側(cè)卵巢，術(shù)后每只肌注8萬單位青霉素抗感染，正常組僅切除卵巢周圍少量脂肪組織，其余操作同卵巢切除組。術(shù)后常規(guī)喂養(yǎng)一周后，健骨顆粒組大鼠按生藥5.8 g/(kg·d)給予健骨顆粒藥液，其余組給予等量生理鹽水灌胃。分別于第6 周、12 周各取材一半，取右側(cè)股骨、脛骨及第四腰椎進(jìn)行骨密度及生物力學(xué)檢測，并隨機(jī)每組選4只大鼠取腹主動(dòng)脈血分離血清。本次動(dòng)物實(shí)驗(yàn)經(jīng)福建中醫(yī)藥大學(xué)動(dòng)物倫理委員會(huì)批準(zhǔn)，并遵守所有與動(dòng)物倫理使用相關(guān)機(jī)構(gòu)和政府的法規(guī)。

1.2.2骨密度檢測：將大鼠右側(cè)股骨放置在帶有小動(dòng)物骨密度測定軟件的骨密度(bone mineral density, BMD)掃描儀探頭下檢測，經(jīng)計(jì)算機(jī)自動(dòng)分析得出大鼠右側(cè)股骨 BMD 相關(guān)數(shù)據(jù)。

1.2.3骨骼生物力學(xué)性能檢驗(yàn):將右側(cè)脛骨及第四腰椎去除周圍軟組織，分別運(yùn)用三點(diǎn)抗彎法和壓縮法測定脛骨與腰椎的最大載荷。

1.2.4色譜及質(zhì)譜檢測：正離子流動(dòng)相：緩沖液 A=0.1 %甲酸、水、10 mmol/L 醋酸銨；緩沖液 B=0.1 %甲酸、水、90 %乙腈，10 mmol/L 醋酸銨。負(fù)離子流動(dòng)相：緩沖液 A=5 %(v/v)乙腈-10 mmol/L 醋酸銨，pH 9；緩沖液 B=90 %(v/v)乙腈-10 mmol/L 醋酸銨，pH 9。梯度洗脫：0～1 min，100 % A，0 % B；9～12 min，0 % A，0 % B；12.1～15 min，100 % A，0 % B；流速0.35 mL/min。質(zhì)量控制(quanlity contral，QC)樣品進(jìn)行正離子掃描及負(fù)離子掃描(全掃描范圍m/z 70～1 050)，離子強(qiáng)度TOP15的母離子進(jìn)行二級質(zhì)譜鑒定。母離子使用HCD方法碎裂，進(jìn)行二級質(zhì)譜序列鑒定，獲得質(zhì)譜檢測原始數(shù)據(jù)。

1.3 數(shù)據(jù)處理

通過Compound Discover(CD)軟件進(jìn)行數(shù)據(jù)提取處理，基于QC樣品TIC圖疊加，查看LC-MS檢測系統(tǒng)的穩(wěn)定性。通過mzCloud和綜合數(shù)據(jù)庫ChemSpider在線搜索鑒定化合物，Pathway Analysis數(shù)據(jù)庫檢索和確認(rèn)差異代謝物，綜合單變量統(tǒng)計(jì)和多變量統(tǒng)計(jì)分析篩選具有潛在生物學(xué)意義的差異代謝物。樣品間歸一化處理及樣品代謝物間標(biāo)準(zhǔn)化處理后，采用正交偏最小二乘判別分析( OPLS-DA)制火山圖。t檢驗(yàn)考察代謝物組間差異的顯著性，運(yùn)用SPSS17.0對差異性代謝物取Log后進(jìn)行單因素方差分析。通過差異代謝物的代謝通路和顯著差異代謝物富集分析篩選潛在的靶標(biāo)路徑。

2 結(jié)果

2.1 骨密度檢測結(jié)果

大鼠右側(cè)股骨骨密度如表1所示，與正常組相比，同期模型組骨密度下降(P<0.01)；與模型組相比，同期健骨顆粒組大鼠骨密度均升高(12周，P<0.01)。

表1 股骨近端骨密度檢測結(jié)果Table 1 The test result of bone mineral density(g/cm2,

2.2 骨生物力學(xué)檢測結(jié)果

大鼠脛骨與腰椎的最大載荷結(jié)果如表2所示。與正常組相比，同期模型組骨組織生物力學(xué)最大載荷均下降(P<0.05或P<0.01)；與模型組比，同期健骨顆粒組大鼠脛骨及第四腰椎的最大載荷均升高(6周腰椎最大載荷P<0.05，12周均為P<0.01)。

表2 各組大鼠脛骨及第四腰椎最大載荷檢測結(jié)果

2.3 基于LC-MS/MS平臺(tái)的大鼠血清代謝組學(xué)檢測

QC樣品的穩(wěn)定性試驗(yàn)考察基于HPLC-MS/MS平臺(tái)正離子模式(positive ion mode，Pos)和負(fù)離子模式(negative ion mode，Neg)兩種電離方式同時(shí)檢測，大大提高代謝物覆蓋率。QC樣品TIC圖疊加顯示LC-MS系統(tǒng)穩(wěn)定性結(jié)果如圖1，各血清樣品中的離子峰分離較明顯。

圖1 大鼠血清樣品TIC圖Fig.1 TIC picturefrom LC-MS/MS spectra of serum注：A：正離子模式；B：負(fù)離子模式。

2.4 分層聚類分析血液數(shù)據(jù)結(jié)果

分層聚類分析所得熱圖(圖2)直觀顯示出正負(fù)離子通道鑒定的所有血清代謝物在各組之間的含量變化趨勢。圖中水平軸和垂直軸分別代表樣品與變量信息，正常組與其他各組顯著分開。模型組與正常組組間差異提示去卵巢大鼠機(jī)體代謝發(fā)生變化；健骨顆粒組與模型組組間差異明顯，提示中藥對去卵巢大鼠的代謝有調(diào)節(jié)作用。

圖2 分層聚類對比分析各組大鼠血清內(nèi)源性物質(zhì)的相關(guān)性Fig.2 Clustering result of endogenous materials of rats serum metabolites shown as heatmap注：A、E：6周正常組/模型組；B、F: 6周模型組/健骨顆粒組；C、G:12周正常組/模型組；D、H: 12周模型組/健骨顆粒組；A～D為正離子模式；E～H為負(fù)離子模式；-正常組；-模型組；-健骨顆粒組；顏色深淺反映變量值大小，紅色：含量最高，綠色/藍(lán)色：含量最低。

2.5 絕經(jīng)后骨質(zhì)疏松大鼠機(jī)體代謝表型的變化

正常組、模型組及健骨顆粒組組間代謝物譜差異的OPLS-DA結(jié)果如圖3所示，正常組和模型組的血清樣本點(diǎn)在PCI維可以完全分離，且各樣本點(diǎn)在一定范圍內(nèi)聚集較好，說明去卵巢后骨質(zhì)疏松大鼠機(jī)體生理及物質(zhì)代謝情況已經(jīng)發(fā)生了明顯變化。采用 7 倍交叉驗(yàn)證法驗(yàn)證OPLS-DA模型可靠性，結(jié)果見表3。從生理代謝物的變化層面可認(rèn)為本研究建立的模型成功，可靠性較高。

表3 OPLS-DA 模型的主要參數(shù)Table 3 Main parameters of the OPLS-DA model

圖3 大鼠血清OPLS-DA得分圖Fig.3 OPLS-DA score plots derived from LC-MS/MS spectra of serum注：A、E：6周正常組/模型組；B、F: 6周模型組/健骨顆粒組；C、G:12周正常組/模型組；D、H: 12周模型組/健骨顆粒組；A～D為正離子模式；E～H為負(fù)離子模式。

2.6 絕經(jīng)后骨質(zhì)疏松癥相關(guān)潛在生物標(biāo)志物的鑒定及差異代謝物分析

本研究使用OPLS-DA模型中變量重要性投影(VIP)參數(shù)評價(jià)且經(jīng)t檢驗(yàn)篩選差異性代謝物，結(jié)果見圖4。經(jīng)t檢驗(yàn)篩選出27種差異代謝物作為潛在的生物標(biāo)志物并進(jìn)行單因素方差分析。研究發(fā)現(xiàn)，與正常組相比，模型組中脂質(zhì)代謝及嘧啶代謝的相關(guān)代謝產(chǎn)物乙酰左旋肉堿、花生四稀酸、硬脂酸、膽汁酸、皮質(zhì)酮及尿苷、胞嘧啶核苷、尿嘧啶等含量均下降，與模型組相比，健骨顆粒組其含量均上升；模型組中氨基酸代謝相關(guān)的纈氨酸及L-精氨酸的含量均上升，健骨顆粒組含量與模型組相比均下降(表4)。其中，16-羥基十六烷酸、硬脂酸、泛酸等組間對比差異顯著(P<0.05)，具體結(jié)果見圖5。

表4 歸一化處理血清與絕經(jīng)后骨質(zhì)疏松模型相關(guān)的生物標(biāo)志物的信息Table 4 Normalized serum biomarker information related to postmenopausal osteoporosis model

圖4 差異代謝物篩選火山圖Fig.4 Volcano picture of DEM about potentially different metabolites注：A、E：6周正常組/模型組；B、F: 6周模型組/健骨顆粒組；C、G:12周正常組/模型組；D、H: 12周模型組/健骨顆粒組；A～D為正離子模式；E～H為負(fù)離子模式；紅色：顯著上調(diào)的代謝物；藍(lán)色：顯著下調(diào)的代謝物；黑色：閾值外的代謝物。

圖5 具有顯著差異性的代謝產(chǎn)物Fig.5 Metabolites with significant difference注：組間對比顯著差異性，*P<0.05；組間對比顯著差異性，#P<0.005。

2.7 血清中差異代謝物相關(guān)代謝通路分析

經(jīng)路徑分析，結(jié)合“2.6”項(xiàng)結(jié)果對所富集的代謝通路篩選見表5。泛酸和輔酶A、初級膽汁酸、類固醇激素及不飽和脂肪酸的生物合成、嘧啶代謝和其他次生代謝產(chǎn)物的生物合成影響值大于0.10，提示OPLS-DA所篩選出的差異代謝物與這些代謝途徑密切相關(guān)，說明健骨顆粒發(fā)揮治療作用與這些代謝途徑有關(guān)。

表5 構(gòu)建分析通路結(jié)果(Impact > 0.10)Table 5 Result from pathway analysis (Impact>0.10)

3 討論

卵巢切除術(shù)是經(jīng)典的絕經(jīng)后骨質(zhì)疏松癥動(dòng)物模型制作方法，結(jié)果顯示，本次實(shí)驗(yàn)造模成功，藥物干預(yù)后骨密度與骨生物力學(xué)指標(biāo)均有所改善，表明健骨顆粒對去卵巢后骨質(zhì)疏松癥有治療作用。

絕經(jīng)后體內(nèi)雌激素水平下降，炎性因子濃度及氧化應(yīng)激標(biāo)志物均升高[8-9]。動(dòng)物實(shí)驗(yàn)表明，去卵巢后大鼠體內(nèi)雌激素水平下降，從而引起代謝機(jī)制紊亂及內(nèi)源性代謝物含量變化，尤其是體內(nèi)脂代謝。此外，雌激素可減少膽汁酸排出[10-11]，緩解蛋白消耗，調(diào)節(jié)多種炎性疾病[12-15]。代謝組學(xué)研究[16-18]發(fā)現(xiàn)，去卵巢后骨質(zhì)疏松大鼠體內(nèi)雌激素水平下降，可引起脂代謝、核酸代謝及氨基酸代謝紊亂，導(dǎo)致肥胖、脂肪肝、骨質(zhì)疏松等一系列疾病。中藥干預(yù)可調(diào)節(jié)大鼠機(jī)體的脂代謝、氨基酸代謝及核酸代謝等[19-20]，使機(jī)體恢復(fù)正常。本研究結(jié)果與文獻(xiàn)報(bào)道大體一致，去卵巢大鼠體內(nèi)脂代謝、核酸代謝和氨基酸代謝異常，健骨顆粒干預(yù)后大鼠體內(nèi)脂代謝、核酸代謝及氨基酸代謝得到調(diào)節(jié)，機(jī)體有恢復(fù)正常的趨勢，其中脂類代謝的變化尤其顯著。

脂類代謝包括脂肪酸代謝、初級膽汁酸合成、類固醇生物合成、花生四稀酸代謝、神經(jīng)鞘磷脂的合成等。脂肪酸類具有抗炎功能[21-22]。膽汁酸可影響維生素D的腸道吸收，對成骨細(xì)胞有破壞作用,影響進(jìn)鈣、磷吸收，新骨形成及舊骨溶解，致低骨密度及骨流失，引起膽汁淤積性骨質(zhì)疏松癥[23-24]。與模型組相比，健骨顆粒組中這些代謝產(chǎn)物的含量有恢復(fù)正常的趨勢，但6周膽汁酸含量變化除外。推測其因卵巢切除后體內(nèi)雌激素含量下降致膽汁酸代謝增強(qiáng)，所以模型組膽汁酸含量明顯升高。此外，乙酰左旋肉堿具有運(yùn)輸脂肪的能力，與正常組相比，模型組中其含量下降，與模型組相比，健骨顆粒組含量上升，表明藥物干預(yù)可調(diào)節(jié)機(jī)體脂代謝。脂代謝相關(guān)物質(zhì)的含量變化說明去卵巢骨質(zhì)疏松大鼠機(jī)體脂代謝紊亂，健骨顆?？赏ㄟ^改善脂代謝對絕經(jīng)后骨質(zhì)疏松起一定的干預(yù)調(diào)節(jié)作用。

本次篩選的相關(guān)代謝產(chǎn)物還有皮質(zhì)酮、纈氨酸、泛酸、大豆苷元以及與核酸代謝相關(guān)的尿嘧啶、尿苷及胞嘧啶核苷。其中，皮質(zhì)酮可誘導(dǎo)氧化應(yīng)激，與花生四稀酸相互作用干擾正常生理行為[21]，致成骨細(xì)胞凋亡增加,分化能力減弱及雌激素護(hù)骨作用降低[8,25]。纈氨酸是人體必需氨基酸之一及生糖氨基酸。其缺乏會(huì)抑制嗜中性與嗜酸性白細(xì)胞增殖[26]，影響免疫功能及抑制炎癥調(diào)節(jié)功能。泛酸是脂肪酸合成類固醇所必需的物質(zhì)，是合成輔酶A(CoA)的重要成分，CoA參與生物體內(nèi)糖類、脂肪、蛋白質(zhì)的代謝[27]。大豆苷元可促進(jìn)人成骨樣細(xì)胞分化,使骨量增加,骨結(jié)構(gòu)改善,可改善骨質(zhì)疏松癥狀[28-29]。核酸類物質(zhì)具有遺傳、介導(dǎo)和催化生化反應(yīng)、提供或轉(zhuǎn)移能量等功能。其含量變化說明模型組機(jī)體核酸代謝及氨基酸代謝異常，而健骨顆粒可改善去卵巢后骨質(zhì)疏松機(jī)體核酸代謝及氨基酸代謝。

綜上，基于LC-MS技術(shù)的代謝組學(xué)研究發(fā)現(xiàn)，去卵巢后大鼠體內(nèi)脂代謝、嘧啶代謝及氨基酸代謝等紊亂，健骨顆粒干預(yù)可改善其代謝情況。因此，健骨顆粒可能通過調(diào)節(jié)多種代謝通路從而改善絕經(jīng)后骨質(zhì)疏松癥。