密蒙花揮發(fā)油成分的定量結(jié)構(gòu)與色譜保留關(guān)系研究

秦正龍 冷檢 馮長君

摘要：采用電拓?fù)鋺B(tài)指數(shù)（Ek）及分子價連接性指數(shù)（mXtV）對從密蒙花 （Buddleja officinalis Maxim）揮發(fā)油中分離出的27種化合物進(jìn)行結(jié)構(gòu)表征，通過多元線性回歸及最佳變量子集方法，建立了色譜保留時間的定量構(gòu)效關(guān)系模型。該模型的相關(guān)系數(shù)（R2）為0.988，逐一剔除法交叉驗證系數(shù)（Q2）為0.982，計算值與實驗值頗為吻合，較好地揭示了密蒙花揮發(fā)油成分色譜保留時間的遞變規(guī)律，該研究對中藥化學(xué)成分的分離、純化、檢測、結(jié)構(gòu)表征及其藥用價值的開發(fā)，具有理論指導(dǎo)意義。

關(guān)鍵詞：密蒙花；揮發(fā)油；保留時間；電拓?fù)錉顟B(tài)指數(shù)；分子價連接性指數(shù)

中圖分類號：TS201.2 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2015）03-0674-02

DOI：10.14088/j.cnki.issn0439-8114.2015.03.043

Relationship between Quantitative Structure and Retention Time of Volatile Components from Buddleja officinalis Maxim

QIN Zheng-long1，LENG Jian1，F(xiàn)ENG Chang-jun2

（1.School of Chemistry & Chemical Engineering， Jiangsu Normal University， Xuzhou 221116， Jiangsu， China；

2.School of Chemistry & Chemical Engineering， Xuzhou Institute of Technology， Xuzhou 221008， Jiangsu， China）

Abstract： The atom-type electrotopological state index （Ek） and molecular valence connectivity index （mXtV） were used to characterize the structures of 27 volatile compounds from Buddleja officinalis Maxim. A quantitative structure-retention relationship（QSRR） between Ek，mXtV and chromatographic retention time （RT） of volatile components from Buddleja officinalis Maxim was developed by using multiple linear regression and Leaps-and-Bounds regression. Multiple correlation coefficient （R2） and cross-validation correlation coefficient （Q2） of leave-one-out were 0.988 and 0.982 respectively. The calculated values were in good agreement with experimental data， and the model can elucidate well the change of chromatographic retention time for volatile compounds from Buddleja officinalis Maxim. The method may play an important role in the theoretical guidance， such as separation， identification and confirming the structure of chemical composition from the Chinese medicine and exploiting their medical value.

Key words： Buddleja officinalis Maxim；volatile component；retention time；atom-type electrotopological state index；molecular valence connectivity index

密蒙花（Buddleja officinalis Maxim）又名蒙花、染飯花、小錦花、酒藥花、羊耳朵、斷腸草，是馬錢科醉魚草屬植物密蒙花的干燥花蕾及其花序[1]，廣泛分布于中國華北、華南、西南等地區(qū)[2]。密蒙花是傳統(tǒng)中藥，味甘，性微寒，具有涼血、養(yǎng)肝、祛風(fēng)、明目之功效，主要用于目赤腫痛，多淚羞明，肝虛目暗，視物昏花[3，4]。此外，密蒙花還有多種藥理活性，具有降低血糖、調(diào)節(jié)免疫、抗炎、抗菌、抗氧化及抗腫瘤等功效[5，6]。

目前，有關(guān)密蒙花的研究報道主要集中在它的化學(xué)成分及生物活性[7，8]，而對其分子結(jié)構(gòu)與性質(zhì)之間的相關(guān)關(guān)系的研究報道則較少。為此，本試驗利用電拓?fù)鋺B(tài)指數(shù)（Ek）[9]及分子價連接性指數(shù)（mXtV）[10]對密蒙花中的27種揮發(fā)油成分進(jìn)行結(jié)構(gòu)表征，并將它們與色譜保留時間（RT）進(jìn)行回歸分析，建立了定量結(jié)構(gòu)色譜保留關(guān)系模型，對預(yù)測色譜保留值、選擇色譜分離條件、解釋色譜分配機(jī)理等具有重要意義。

1 材料與方法

1.1 數(shù)據(jù)來源

密蒙花揮發(fā)油中27種化學(xué)成分的保留時間取自文獻(xiàn)[11]。

1.2 兩類拓?fù)渲笖?shù)的計算

應(yīng)用ChemDraw Ultra 9.0構(gòu)建密蒙花中27種揮發(fā)油成分的分子結(jié)構(gòu)，保存為mol文件，在Matlab7.1.0下調(diào)用其分子結(jié)構(gòu)，應(yīng)用修正的程序[12]計算得到46種電拓?fù)鋺B(tài)指數(shù)（Ek），12種分子價連接性指數(shù)（mXtV）。

1.3 方法

將計算得到的兩類拓?fù)渲笖?shù)作為自變量，相應(yīng)的RT作為因變量，利用最佳變量子集回歸選擇最優(yōu)變量組合，建立QSRR模型，采用逐一剔除法對模型的預(yù)測能力及魯棒性進(jìn)行檢驗，以交叉驗證系數(shù)（Q2）予以評價。如果Q2≥0.5，說明所建模型具有良好的魯棒性及預(yù)測能力[13，14]。

2 結(jié)果與討論

2.1 密蒙花揮發(fā)油成分QSRR模型的構(gòu)建

將密蒙花中27種揮發(fā)油成分的RT與計算得到的兩類拓?fù)渲笖?shù)一起輸入SPSS 13.0，在95%的置信區(qū)間內(nèi)進(jìn)行回歸分析，應(yīng)用最佳變量子集建立的QSRR模型見表1。其中R，R2，R2adj，Q2，S，F(xiàn)分別為相關(guān)系數(shù)、判定系數(shù)、校正判定系數(shù)、逐一剔除法的交叉驗證系數(shù)、估計標(biāo)準(zhǔn)誤差、Fischer檢驗值。

由表1可知，隨著模型中變量數(shù)目的增多，其R、R2、R2adj均逐漸增大，S逐漸減小，NO.4后R、R2、R2adj增幅極小，而Q2先逐漸增大至0.982，而后又下降。此外，為了使構(gòu)建的模型具有良好的可信度，樣本數(shù)與變量數(shù)的比值不能太小，所建模型才有意義[15]。綜合考慮，確定四元最佳QSRR模型為：

RT=-10.727+7.4491XVP+4.2414XVc-1.299E1-1.022E2 （1）

式中，1XVP、4XVc分別為1階路徑分子價連接性指數(shù)和4階簇分子價連接性指數(shù)；E1、E2依次為第1類和第2類原子的電拓?fù)鋺B(tài)指數(shù)。

n=27，R=0.994，R2=0.988，R2adj=0.986，Q2=0.982 S=1.102，F(xiàn)=462.966

用上述模型給出的計算值與實驗值頗為吻合。

2.2 QSRR模型的魯棒性檢驗

由于所研究的是小樣本（n<30），故選用逐一剔除法[16]。每次只剔除一個化合物，用余下的化合物作為訓(xùn)練集，按模型（1）的形式回歸，重復(fù)27次，得到27個Jackknifed相關(guān)系數(shù)（R2）。根據(jù)R2的雷達(dá)圖評價模型的穩(wěn)健性。如圖1，以0.980為圓心，0.003為間距，27個Jackknife的R2全部落在0.983～0.990之間，波動很小。表明模型（1）中不存在異常的RT數(shù)據(jù)及隨機(jī)相關(guān)，具有總體可接受的穩(wěn)健性及良好的預(yù)測能力。另外，模型（1）的R2adj與Q2相差0.004，遠(yuǎn)遠(yuǎn)小于0.3，說明模型沒有過擬合，不存在不相關(guān)的其他變量，數(shù)據(jù)中不存在離域點。

2.3 模型的構(gòu)效關(guān)系分析

色譜保留時間（RT）主要是由組分和固定相之間的作用力決定[17]，作用力越大，則RT越大。分子間的作用力包括色散力、誘導(dǎo)力、取向力和氫鍵。分子價連接性指數(shù)揭示了分子的大小、形狀、分枝，反映了色散力的大??；電拓?fù)鋺B(tài)指數(shù)蘊(yùn)含了分子中各原子的電子狀態(tài)及電性作用，可表征誘導(dǎo)力、取向力的強(qiáng)弱，它們共同揭示了影響RT的本質(zhì)因素。

電拓?fù)鋺B(tài)指數(shù)（Ek）及分子價連接性指數(shù)（mXtV）較好地表征密蒙花揮發(fā)油組分的結(jié)構(gòu)特征，所建的QSRR模型具有良好的穩(wěn)健性及預(yù)測能力，揭示了密蒙花揮發(fā)油成分RT的遞變規(guī)律，對天然植物化學(xué)成分的分離、分析、結(jié)構(gòu)表征、質(zhì)量檢驗及藥用價值的開發(fā)，具有理論指導(dǎo)意義。

參考文獻(xiàn)：

[1] 王軍憲，李教社，楊廣德，等.高效液相色譜法測定密蒙花中蒙花苷的含量[J].藥物分析雜志，2001，21（2）：103-104.

[2] 熊 勇，熊揚(yáng)波，楊青松.藥用植物密蒙花總黃酮提取及抗氧化性研究[J].生物技術(shù)，2011，21（3）：85-87.

[3] 李懷鳳，彭清華，姚小磊，等.密蒙花總黃酮對去勢雄鼠干眼癥模型角膜和淚腺組織的保護(hù)作用[J].中國中醫(yī)眼科雜志，2010，20（1）：1-6.

[4] 張小曼，彭永芳，馬銀海.正交法微波輻射提取密蒙花黃色素[J].食品科學(xué)，2002，23（10）：91-93.

[5] 王 方，彭清華，姚小磊，等.密蒙花總黃酮含藥血漿干預(yù)干眼癥細(xì)胞凋亡模型STAT1磷酸化蛋白表達(dá)[J]. 國際眼科雜志，2010， 10（1）：5-8.

[6] 崔 穎，張永旺.密蒙花研究進(jìn)展[J]. 甘肅中醫(yī)學(xué)院學(xué)報，2010， 27（2）：65-68.

[7] 韓 澎，崔亞君，郭洪祝，等.密蒙花化學(xué)成分及其活性研究[J].中草藥，2004，35（10）：1086-1090.

[8] 郭 雷，朱文成，劉 超.密蒙花化學(xué)成分及生物活性研究進(jìn)展[J].食品研究與開發(fā)，2012，33（7）：222-225.

[9] HALL L H， KIER L B， BROWN B B. Molecular similarity based on novel atom-type electrotopological state indices[J].Journal of Chemical Information and Computer Sciences，1995，35（6）：1074-1080.

[10] KIER L B， HALL L H. Molecular connectivity in structure-activity analysis[M].New York：John Wiley & Sons，1986.69-75.

[11] 張?zhí)m勝，董光平，劉光明.密蒙花揮發(fā)油化學(xué)成分的研究[J].安徽農(nóng)業(yè)科學(xué)，2010，38（9）：4585-4586.

[12] 黔 南，梁逸曾，王亞麗.直觀隊列命名法的基本原理及其在矩陣與拓?fù)渲笖?shù)計算中的應(yīng)用[J].計算機(jī)與應(yīng)用化學(xué)，2003，20（4）：386-389.

[13] DOUGLAS M H， SUBHASH C B， DENISE M. Assessing model fit by cross-validation[J]. Journal Chemical Information and Computer Sciences， 2003， 43（2）：579-586.

[14] WINKLER D A， BURDEN F R. Holographic QSAR of benzodiazepines[J].Quantitative Structure-Activity Relationships，1998，17（3）：224-231.

[15] 劉 東，章文軍，許 祿.手性羥酸和氨基酸類化合物的構(gòu)效關(guān)系研究[J].化學(xué)學(xué)報，2009，67（2）：145-150.

[16] DIETRICH W S， DREYER N D， HANSCH C， et al.Confidence interval estimators for parameters associated with quantitative structure/activity relationship[J]. Journal of Medicinal Chemical，1980，23（11）：1201-1205.

[17] 秦正龍，馮長君.脂肪醇?xì)庀嗌V保留指數(shù)與結(jié)構(gòu)的相關(guān)性研究[J].色譜，2004，22（4）：452-455.