不同部位獨一味的NMR代謝組模式識別研究

潘 正， 毛 慶， 江 生， 范 剛， 張 藝

（1.成都中醫(yī)藥大學民族藥學院，四川　成都　611130；2.重慶市食品藥品檢驗所，重慶　401121）

不同部位獨一味的NMR代謝組模式識別研究

潘 正1， 毛 慶2， 江 生2， 范 剛1， 張 藝1

（1.成都中醫(yī)藥大學民族藥學院，四川成都611130；2.重慶市食品藥品檢驗所，重慶401121）

目的 運用1H-NMR代謝組學方法研究藏藥獨一味地上部分與地下部分的化學成分差異。方法 獨一味加入CD3OD后超聲提取，將1H-NMR測定的信息進行主成分分析（PCA）和PLS-DA法分析。結(jié)果 數(shù)據(jù)分析表明獨一味地下部分的代謝產(chǎn)物中葡萄糖、蔗糖、果糖的含有量高于地上部分，而脂肪酸低于地上部分。結(jié)論 獨一味的PCA載荷圖中部分的信號仍然難以歸屬，須進一步結(jié)合MS圖譜進行分析。

1H-NMR；獨一味；地上部分；地下部分；葡萄糖；蔗糖；果糖；脂肪酸；主成分分析

唇形科獨一味屬植物獨一味Lamiophlomis rotata（Benth.）Kudo，生長于海拔3 000 m以上的高山草原、碎石灘中，原是藏醫(yī)藥學臨床常用藥材，別名達布巴、達巴巴、吉布孜等，自《中國藥典》1995年版［1］收錄其制劑以來，以其良好的止血鎮(zhèn)痛效果，迅速成為次于云南白藥，年產(chǎn)值超億元的中成藥品種。2005年，獨一味全草作為藥材首次收錄于《中國藥典》［2］，在此以前，獨一味主要以地下部分入藥，由于僅采挖地下部分破壞高原草地植被，為了更好地利用資源，將全草收錄為藥材，即便如此，采挖地下部分對高原草地生態(tài)環(huán)境的破壞，長時間難以恢復(fù)，因此《中國藥典》2010年版［3］將獨一味的干燥地上部分列為藥用部位，基本禁止采挖地下部分。

但不同的藥用部位其化學成分組成和藥理活性不同，迄今為止，獨一味藥材地上部分與地下部分的化學成分組成差異未見報道。為了全面反映兩個不同部位化學成分組成差異，采集西藏、四川、甘肅等12個產(chǎn)地獨一味樣本，對其不同部位提取物進行1H-NMR譜研究，通過化學模式識別法進行主成分分析［4-6］，尋找獨一味不同部位之間的相似性和差異性，為獨一味開發(fā)應(yīng)用等研究提供理論依據(jù)。

1　儀器與試藥

Varian Mercury Plus 400 MHz核磁共振波譜儀（美國Varian公司，第三軍醫(yī)大學-重慶化醫(yī)集團控股聯(lián)合實驗室）；Sartorius BP221S電子天平（北京賽多利斯科學儀器有效公司）；KQ-300B型超聲波清洗器（昆山市超聲儀器有限公司）；氘代甲醇（CD3OD，99.8%）、重水（D2O，99.9%）均購自北京京云重輕科技有限公司，3-三甲基甲硅烷基-2，2，3，3-四氘代丙酸鈉（TSP，美國Sigma-Aldrich公司）。

本實驗所用藥材經(jīng)成都中醫(yī)藥大學張藝教授鑒定為唇形科獨一味屬植物獨一味Lamiophlomis rotata（Benth.）Kudo.的干燥全草。

2　方法

2.1供試液的制備 取獨一味地上部分和地下部分，粉碎至60目。精密稱取約100.0 mg，分別置于2 mL EP管中，分別加入0.8 mL CD3OD，超聲提取30 min，離心，吸取上清液，轉(zhuǎn)移至核磁管中，備用。

2.21H-NMR測定及圖譜處理［7-8］在25℃、35%環(huán)境濕度，采用標準脈沖程序采集圖譜，進行基線校正和相位校正。質(zhì)子觀察頻率400.64 MHz，譜寬6 410.3 Hz，弛豫時間為2 s，掃描次數(shù)64，采集時間為2.556 s，調(diào)用水峰壓制序列壓制水峰信號，總測定時間為6 min 54 s。將核磁文件導入MestReNova軟件，調(diào)整基線后，從δ8.00～0.00以0.04為間隔進行積分，對數(shù)據(jù)進行歸一化處理，再對數(shù)據(jù)進行標準化，最后去掉CD3OD溶劑峰和水峰數(shù)據(jù)，生成一個24×187的矩陣，矩陣的列代表不同的樣品。

2.3數(shù)據(jù)分析 將數(shù)據(jù)預(yù)處理結(jié)果通過SIMCA-P 11.5進行多維數(shù)據(jù)的分析。數(shù)據(jù)的多元統(tǒng)計分析采用主成分分析（PCA）方法進行模式識別。建立模型后，分別產(chǎn)生得分矩陣圖和載荷矩陣圖，通過得分矩陣圖來獲得各組數(shù)據(jù)的分型信息。而載荷矩陣圖能反映各個變量對樣品分型的影響。

采用偏最小二乘法判別分析（partial least squares discriminant analysis，PLS-DA）［9］法進行模式識別。

3　結(jié)果與討論

3.1不同部位化學成分PCA分析 為了區(qū)分獨一味不同部位化學成分的差異、找到獨一味地上部分代替地下部分使用的理論依據(jù)，采用了PCA模式識別技術(shù)對其進行分析，所得PCA得分圖見圖1，結(jié)果表明，前2個主成分的累積貢獻率達到94.7%，說明這2個主成分已經(jīng)可以對94.7%以上的原變量信息進行解釋，故選取這2個主成分進行得分圖和載荷圖的分析。根據(jù)圖中兩組樣本的空間投影距離和混雜程度來看，獨一味地上部分與地下部分代謝產(chǎn)物在PC1軸上均存在顯著差別，而在PC2軸上差異不大。

圖1　獨一味藥材不同部位1H-NMR的PCA分析得分Fig.1 PCA score p lot in different parts of Lamiophlomis rotata

圖2　獨一味藥材不同部位1H-NMR的PCA分析載荷圖Fig.2 PCA loading p lot in different parts of Lam iophlom is rotata

為了區(qū)分這些不同部位的代謝產(chǎn)物，采用PCA載荷圖進行篩查分析。根據(jù)載荷圖（圖2）對差異性代謝成分進行指認，發(fā)現(xiàn)主要為初生代謝產(chǎn)物信號，葡萄糖的信號峰δ3.78（α-anomer）、3.62、3.46，蔗糖的H-1信號峰δ5.42（H，d，J=3.8 Hz），果糖的H-3信號峰δ4.14（H，d，J=8.6 Hz），脂肪酸的信號峰δ1.31（br.s）。結(jié)果顯示獨一味地上地下部位化學成分差異主要表現(xiàn)為，地下部分的代謝產(chǎn)物中葡萄糖、蔗糖、果糖的量高于地上部分，而脂肪酸低于地上部分。地下部分為植物的營養(yǎng)器官，貯藏的糖類物質(zhì)較多，地上部分為植物的莖葉、花序及種子等，具有生長運輸養(yǎng)料和生殖繁殖功能，因此含有的脂肪酸含量較高。

前期研究發(fā)現(xiàn)獨一味中次生代謝產(chǎn)物1H-NMR的指紋區(qū)域主要集中在δ7.80～4.00，對該區(qū)域篩查分析，發(fā)現(xiàn)與PC1正相關(guān)的還有δ7.59、6.69、6.61、4.91、4.74、4.34，與PC1負相關(guān)的主要有δ7.56、7.45、6.65、5.34、5.14和4.51。結(jié)合前期研究，進一步對其進行信號歸屬，從中發(fā)現(xiàn)毛蕊花糖苷中咖啡?；|(zhì)子信號δ7.59，以及苯環(huán)質(zhì)子信號δ6.69、6.61；δ7.56和7.45分別為胡麻屬苷和8-O-乙酰山梔苷甲酯的H-3質(zhì)子信號。雖然上述信號與葡萄糖、蔗糖、果糖和脂肪酸等化合物信號相比，在整個載荷圖中幾乎可以忽略不計，但上述信號的分析，對考察藥材中次生代謝產(chǎn)物差異有重要意義，研究發(fā)現(xiàn)獨一味地上地下部位次生代謝產(chǎn)物差異主要表現(xiàn)為，地下部分的毛蕊花糖苷含量高于地上部分，而胡麻屬苷和8-O-乙酰山梔苷甲酯低于地上部分。

PCA為無監(jiān)督的分析方法，反映了數(shù)據(jù)的原始狀態(tài)，用來觀察實驗樣品的自然分布和組別關(guān)系。進一步區(qū)分不同組別差異，將組內(nèi)的差異降至最低，更能把握多維數(shù)據(jù)的整體特征和變異規(guī)律，采用有監(jiān)督的PLS-DA分析對數(shù)據(jù)進行分析［10-11］。

3.2PLS-DA分析 在樣品的PLS-DA統(tǒng)計分析載荷圖（圖3）中，各個成分距離原點的距離表明了其對判別貢獻的大小，距原點越遠，對樣品分類的貢獻越大，即潛在的代謝標記物。在作為變量的色譜峰中，δ3.78、3.82、3.74、3.62、1.31、1.35、1.11、4.11各峰距離原點較遠，對樣品分類貢獻較大，與PCA分析結(jié)果趨勢相符。

4　結(jié)論

本研究以1H-NMR為識別變量，運用PCA和PLS-DA方法對不同部位的獨一味樣品進行分析，均可實現(xiàn)獨一味的部位分類。從代謝組整體化學信息來看，地上部分與地下部分獨一味藥材間的差異客觀存在，植物中初生代謝產(chǎn)物如葡萄糖、蔗糖、果糖在獨一味地下部位中質(zhì)量分數(shù)較高，而脂肪酸在獨一味地上部位中質(zhì)量分數(shù)較高。通過對獨一味地上地下部位1H-NMR光譜的PCA和PLS-DA分析，結(jié)果提示，獨一味不同部位的化學成分組成主要表現(xiàn)為初生代謝產(chǎn)物的差異，為獨一味地上部位代替地下部位入藥提供了初步的理論支持。

圖3　獨一味藥材1H-NMR的PLS-DA載荷圖Fig.3 PLS-DA loading p lot in different parts of Lam iophlom is rotata

獨一味中含有的次生代謝產(chǎn)物［12-13］如環(huán)烯醚萜、黃酮、苯乙醇苷等在不同的部位肯定存在差異，但在PCA載荷圖譜中相應(yīng)區(qū)域信號對樣品分類貢獻很小，雖然發(fā)現(xiàn)毛蕊花糖苷、胡麻屬苷和8-O-乙酰山梔苷甲酯等成分在藥材地上地下部位含量存在差異，其可靠性仍需結(jié)合其他分析方法進一步驗證；此外，前期雖然對獨一味地上部分的8個次生代謝產(chǎn)物進行了1H-NMR信號歸屬（已整理，另文發(fā)表），但對于PCA載荷圖中δ6.65、5.34、5.14、4.91、4.74、4.51和4.34的信號仍然難以歸屬，下一步課題組將針對上述問題，擬結(jié)合HPLC和MS指紋圖譜綜合分析，進一步發(fā)現(xiàn)其差異性成分，對獨一味的代謝產(chǎn)物差異進行更為全面的研究。

［1］ 國家藥典委員會.中華人民共和國藥典：1995年版一部［S］.北京：化學工業(yè)出版社，1995：555.

［2］ 國家藥典委員會.中華人民共和國藥典：2005年版一部［S］.北京：化學工業(yè)科技出版社，2005：184-185.

［3］ 國家藥典委員會.中華人民共和國藥典：2010年版一部［S］.北京：中國醫(yī)藥科技出版社，2010：245-246.

［4］ Simmler C，Napolitano J G，McAlpine J B，et al.Universal quantitative NMR analysis of complex natural samples［J］.Curr Opin Biotech，2014，25（4）：51-59.

［5］ Kima J，Jung Y，Song B，et al.Discrimination of cabbage（Brassica rapa ssp.pekinensis）cultivars grown in different geographical areas using1H NMR-based metabolomics［J］.FoodChem，2013，137（1-4）：68-75.

［6］ Caligiani A，Palla L，Acquotti D，et al.Application of1H NMR for the characterisation of cocoa beans of different geographical origins and fermentation levels［J］.Food Chem，2014，157（1）：94-99.

［7］ Lee JE，Lee B J，Chung JO，et al.1H NMR-basedmetabolomic characterization during green tea（Camellia sinensis）fermentation［J］.Food Res Int，2011，44（2）：597-604.

［8］ Yanga SO，Shin Y S，Hyun SH，etal.NMR-based metabolic profiling and differentiation of ginseng roots according to cultivation ages［J］.JPharmaceut Biomed，2012，58（1）：19-26.

［9］ Lin W S，Yang C M，Kuo B J.Classifying cultivars of rice（Oryza satiυa L.）based on corrected canopy reflectance spectra data using the orthogonal projections to latent structures（OPLS）method［J］.Chemometr Intell Lab，2012，115（6）：25-36.

［10］ Godoy JL，Vega JR，et al.Relationships between PCA and PLS-regression［J］.Chemometr Intell Lab，2014，130（15）：182-191.

［11］ Lavine B K，RayensW S.Referencemodule in chemistry，molecular sciences and chemical engineering，from Comprehensiυe Chemometrics：chemical and biochemical data analysis［M］. Waltham，MA：Elsevier，2009：517-540.

［12］ LiM X，Shang X F，Zhang R X，etal.Antinociceptive and anti-inflammatory activities of iridoid glycosides extract of Lamiophlomis rotata（Benth.）Kudo［J］.Fitoterapia.2010，81（3）：167-172.

［13］ 張 鳳，孫連娜，陳萬生.獨一味的化學成分及藥理作用［J］.藥學實踐雜志，2008，26（3）：169-171.

M etabonom ics and pattern recognition study on the different parts of Lamiophlom is rotata（Benth.）Kudo

PAN Zheng1， MAO Qing2， JIANG Sheng2， FAN Gong1， ZHANG Yi1
（1.College of Ethnic Medicine，Chengdu Uniυersity of Traditional Chinese Medicine，Chengdu 611130，China；2.Chongqing Institute for Food and Drug Control，Chongqing 401121，China）

AIM To study the distribution characteristics of components from above-and underground parts of Lamiophlomis rotata（Benth.）Kudo.METHODS The powdered materialswere extracted with CD3OD under ultrasonication.Principal componentanalysis（PCA）and partial least squares discriminantanalysis（PLS-DA）were applied to the analysis of the1H-NMR spectroscopy data of L.rotata.RESULTS The experimental results showed that the levels of glucose，fructose，sucrose in underground parts were much higher than these in aboveground parts.In contrast，the fatty acid level of underground parts was lower than that of above-ground ones. CONCLUSION The loading plot shows that some signals of L.Rotata are difficult to classify and further study combined with MS spectra is necessary.

1H-NMR；Lamiophlomis rotata；above-part；underground part；glucose；fructose；sucrose；fatty acid；principal component analysis（PCA）

R284.1

A

1001-1528（2015）03-0567-04

10.3969/j.issn.1001-1528.2015.03.023

2014-05-26

重慶市自然科學基金項目（CSTC2014jcyjA10004）；重慶市教委科學技術(shù)研究項目（KJ1400433）

潘 正（1976—），男，博士生，副教授，從事中藥化學成分分析。Tel：（023）62460536，E-mail：panzhengbusiness@ aliyun.com