付先軍，王振國，李學博

（山東中醫(yī)藥大學中醫(yī)藥經(jīng)典理論教育部實驗室/中醫(yī)文獻研究所/信息管理學院濟南250355）

基于多維結(jié)構(gòu)描述符的中藥“性-構(gòu)”關系研究＊

付先軍＊＊，王振國＊＊，李學博

（山東中醫(yī)藥大學中醫(yī)藥經(jīng)典理論教育部實驗室/中醫(yī)文獻研究所/信息管理學院濟南250355）

目的：為了解析中藥寒熱藥性與其成分之間的相關關系，從物質(zhì)基礎層面詮釋中藥寒熱藥性的科學內(nèi)涵。方法：篩選《中國藥典》第一部明確有主要有效成分記載的植物類中藥共計334種，中藥主要成分280種，計算每種成分的分子描述符（包括一維、二維、三維及半經(jīng)驗量化計算分子描述符），通過廣義規(guī)則歸納算法（Generalized Rule Induction，GRI）篩選強相關的分子描述符。結(jié)果：與中藥寒熱藥性相關性較強的主要成分分子結(jié)構(gòu)描述符多與分子的能量狀態(tài)有關系，分子能量狀態(tài)越高，與中藥熱性關聯(lián)關系越大；而分子能量越低，與中藥寒性關聯(lián)關系越大；而平性中藥不明顯。除分子能量狀態(tài)外，中藥藥性還可能與其主要成分的脂水分配系數(shù)、拓撲分子極性表面積，碳、氫、氧原子數(shù)，氫鍵供體數(shù)目等分子結(jié)構(gòu)描述符具有一定的關聯(lián)關系。結(jié)論：中藥藥性確實與中藥主要成分的分子結(jié)構(gòu)具有一定的關聯(lián)關系，其中最重要的分子描述符是三維分子描述符中的分子能量，提示中藥藥性作用可能與中藥主要有效成分的分子能量狀態(tài)以及對機體能量代謝的影響而產(chǎn)生的綜合效應有關。

中藥寒熱藥性結(jié)構(gòu)描述符關聯(lián)規(guī)則“性-構(gòu)”關系

中藥藥性理論是中醫(yī)藥理論體系中的重要組成部分，中藥寒熱藥性理論是中藥藥性理論的核心，反映藥物對人體陰陽盛衰，寒熱變化方面的作用傾向，是說明藥物作用性質(zhì)的重要概念之一，對指導中藥的臨床應用具有重要意義，《本草經(jīng)集注》記載“其甘苦之味可略，有毒無毒易知，惟冷熱須明”[1]，突出寒熱藥性的重要性，以現(xiàn)代科學手段、方法揭示和闡明其機理有著重要意義，中藥寒熱藥性理論的現(xiàn)代研究已成為中醫(yī)藥科學研究的熱點與前沿，從20世紀中葉以來，國內(nèi)外眾多學者從中藥寒熱藥性的作用機理、物質(zhì)基礎等方面對中藥四性理論進行了一系列的研究和探討，取得了一定進展。寒熱藥性是中藥有效成分作用于機體的客觀反應，對中藥化學成分與其藥性作用之間相互關系的研究，一直是寒熱藥性物質(zhì)基礎研究的熱點。盛良[2,3]提出化學元素具有寒涼溫熱四性，認為中藥的寒熱溫涼藥性反映了物質(zhì)在化學反應中電子得失（包括偏移）的能力；蔣俊等[4]認為中藥中的成分/組分可表現(xiàn)相應的藥性屬性（寒、熱、溫、涼）；李和光、王曉燕[5,6]總結(jié)出寒性、熱性中藥紫外光譜吸收峰的特征規(guī)律，發(fā)現(xiàn)寒藥、熱藥的主要差別在氯仿和石油醚提取液。

上述研究成果揭示，決定中藥藥性的物質(zhì)基礎是中藥所含物質(zhì)成分，寒涼、溫熱藥性的物質(zhì)基礎可能有共性，既有反映同一藥性的類似結(jié)構(gòu)化合物或成分，而中藥的物質(zhì)成分是由多種不同含量的化合物所組成的，分子是構(gòu)成物質(zhì)的基本單位，化合物內(nèi)部分子的結(jié)構(gòu)特征及原子間的組合方式等結(jié)構(gòu)信息（分子描述符）決定了化合物所表現(xiàn)的性質(zhì)[7-9]。因此可以通過研究化合物的分子結(jié)構(gòu)與其活性、性質(zhì)之間的關系而達到對中藥生物活性、整體藥性的預測和掌握，建立起中藥主要成分的結(jié)構(gòu)與活性、藥性之間的關系紐帶[10]。然而目前報道多從化學成分、微量元素研究了寒熱藥性的物質(zhì)基礎，還沒有以數(shù)據(jù)挖掘方法研究中藥主要成分多維結(jié)構(gòu)描述符與寒熱藥性之間關系的報道。

1 材料與方法

1.1 中藥的篩選

研究藥物篩選依據(jù)：①限于植物藥，其他類別的藥物暫不涉及；②四性記載明確、臨床使用公認、學術無爭議的藥物；③2015版《中華人民共和國藥典》[11]收錄的藥物；④主要有效成分清楚。

1.2 信息來源及標準化處理

所篩選中藥藥性、主要成分等基本信息全部來源于《中國藥典》2015版第一部。所有信息均經(jīng)過標準化處理，中藥藥名以《中國藥典》2015版第一部記載的正名為準。《中國藥典》2015版第一部記載的中藥藥性有大寒、寒、微寒、涼、平、大熱、熱、溫、微溫等，由于中藥藥性中寒與大寒、微寒、涼，熱與大熱、溫、微溫只是程度上的不同，沒有本質(zhì)區(qū)別，因此，中藥四性實質(zhì)上可以看作是寒、熱二性，將《中國藥典》記載中大寒、寒、微寒、涼統(tǒng)一歸為寒性藥，將大熱、熱、溫、微溫統(tǒng)一歸為熱性藥，寒熱不明顯，文獻記載為平性藥的依然歸為平性藥。主要成分信息主要來源于《中國藥典》2015版第一部，同時參考來源于清華同方全文期刊庫、維普全文期刊庫以及美國化學文摘（Chemical Abstracts，CA）數(shù)據(jù)庫的文獻資料，選取在中藥中含量較高、藥效作用與中藥臨床療效具有一定相關性的成分。

1.3 研究方法

1.3.1 結(jié)構(gòu)優(yōu)化和分子描述符的計算

將所有的中藥有效成分的結(jié)構(gòu)進行結(jié)構(gòu)優(yōu)化，優(yōu)化好的分子進行分子描述符計算，如組成描述符、拓撲描述符、電荷描述符、連接指數(shù)描述符等多種描述符，用這些描述符來定量的表達分子的結(jié)構(gòu)特征如原子之間的聯(lián)系、分支、分子形狀、對稱性、分子中電荷分布的特征以及分子的量子化學結(jié)構(gòu)特征等[12,13]。目前約有1 650個分子描述符，剔除相關度大于0.8的分子描述符以及沒有比較明確物理含義的描述符[14]，剩余546個分子描述符，其中一維分子描述符114個，二維分子描述符288個，三維分子描述符55個，半經(jīng)驗量化分子描述符（SEMIEMPIRICAL QM PROPERTY）89個。

1.3.2 廣義規(guī)則歸納算法（Generalized Rule Induction，GRI）

GRI從數(shù)據(jù)中提取一組規(guī)則，找出信息容量最高的規(guī)則，信息容量的度量采用指數(shù)衡量，該指數(shù)把規(guī)則的普適性（Support）和精確性（Confidence）都考慮在內(nèi)[15,16]。

在本研究中主要是通過廣義關聯(lián)規(guī)則挖掘中藥藥性與其主要成分描述符之間的關聯(lián)關系，即發(fā)現(xiàn)中藥主要成分與其藥性的關聯(lián)關系強弱，其中支持度是中藥中主要成分的結(jié)構(gòu)描述符（前項）與中藥藥性（后項）同時出現(xiàn)的頻率。置信度是中藥主要成分結(jié)構(gòu)描述符（前項）出現(xiàn)時，后項中藥藥性出現(xiàn)的概率。

2 結(jié)果

2.1 篩選的中藥及主要成分

《中國藥典》2015版第一部明確有主要有效成分記載的植物類中藥共計334種，占總的植物類中藥65.62%，剔除重復出現(xiàn)以及寒熱藥性中藥都出現(xiàn)的主要有效成分，剩余280種成分，其中出現(xiàn)頻數(shù)較高的化學成有綠原酸、金絲桃苷、毛蕊花糖苷、槲皮素等（表1）。這些成分在植物類中藥中的含量都非常高，最低有0.02%，最高達96%，而且這些成分大多具有明顯的藥理作用，在中藥成分中具有一定代表性。

2.2 基于GRI研究中藥藥性及其主要成分結(jié)構(gòu)描述符關聯(lián)關系

2.2.1 主要成分一維分子描述符與中藥藥性之間的相關性分析

與中藥藥性關聯(lián)程度較高的主要成分結(jié)構(gòu)描述符分析（支持度≥10%，置信度≥60%）。從表2可以看出，中藥的寒熱藥性可能與其主要成分的碳、氫、氧原子數(shù)有一定關聯(lián)關系。

2.2.2 主要成分二維分子描述符與中藥藥性之間的相關性分析

與中藥藥性關聯(lián)程度較高的主要成分結(jié)構(gòu)描述符分析（單項支持度≥10%，置信度≥60%；二項支持度≥10%，置信度≥80%）。從表3可以看出，中藥寒性與其主要成分的氫鍵供體數(shù)目及6元環(huán)的數(shù)量有關，而中藥的熱性與其主要成分原子數(shù)，環(huán)鍵數(shù)及鍵數(shù)具有關聯(lián)關系。

2.2.3 主要成分三維分子描述符與中藥藥性之間的相關性分析

與中藥藥性關聯(lián)程度較高的主要成分結(jié)構(gòu)描述符分析（單項支持度≥5%，置信度≥60%；二項支持度≥10%，置信度≥80%）。中藥寒性主要與其主要成分的Jurs DPSA 1分子描述符（帶電表面局部區(qū)域的差異，即部分陽性溶劑可及表面面積減去部分負溶劑可及表面面積）相關性比較大；中藥熱性與Jurs FPSA 2（描述局部電荷在分子表面上的分布情況）和Jurs DPSA 1分子描述符相關性較大；當中藥主要有效成分Jurs DPSA 1分子描述符大于565.19時，多為寒性中藥；成分Jurs DPSA 1分子描述符小于502.77多為熱性中藥（表4）。

2.2.4 主要成分半經(jīng)驗量化分子描述符（SEMIEMPIRICALQM PROPERTY）與中藥藥性之間的相關性分析

與中藥藥性關聯(lián)程度較高的主要成分結(jié)構(gòu)描述符分析（單項支持度≥10%，置信度≥60%；二項支持度≥10%，置信度≥70%）。從表5結(jié)果來看，寒性中藥主要成分生成熱（Heat of Formation）多小于-107.86，而熱性中藥主要成分的總能量（Total Energy）多大于-2 860.58。

3 討論

在上述結(jié)果中可以看出，與中藥寒熱藥性相關性較強的主要成分分子結(jié)構(gòu)描述符多與分子的能量狀態(tài)有關系，而且分子能量狀態(tài)越高，與中藥熱性關聯(lián)關系越大；而分子能量越低，與中藥寒性關聯(lián)關系越大。

而早在20世紀60-70年代，就有學者從藥物對機體能量代謝方面的影響來研究中藥寒熱藥性，如侯燦[17]認為“寒”或“熱”是以熱量不足或熱量過剩為其共同發(fā)病學原因的一種機體典型反應狀態(tài)，寒性熱性藥可改變這種狀態(tài)；70年代，日本學者桑木崇秀[18]證明不同藥性中藥對整體能量代謝的影響是不同的；這些研究結(jié)果證明了寒涼、溫熱性中藥對機體能量代謝的影響是有差異的，溫熱藥多能刺激代謝，增加體內(nèi)熱量產(chǎn)生，而寒涼藥能減少體內(nèi)產(chǎn)熱；而肖小河課題組[19,20]提出了“中醫(yī)藥（藥性）熱力學觀”（thermodynamic outlook of Chinese medicine），從體外到體內(nèi)、從實驗到臨床，較系統(tǒng)地考察了機體新陳代謝過程中的能量轉(zhuǎn)移和熱變化，及其不同寒熱藥性中藥的干預效果，揭示中藥寒熱藥性差異的客觀性以及“寒者熱之、熱者寒之”的科學性，為審視和研究中藥藥性理論提供了新的視角和方法；樊冬麗等[21]則對麻黃湯和麻杏石甘湯進行了微量熱學比較，結(jié)果發(fā)現(xiàn)溫熱藥方麻黃湯作用于細菌生長代謝時，所釋放的熱量要高于麻杏石甘湯。

表1 來源于不同藥性的中藥主要成分

表2 主要成分一維分子描述符與中藥藥性之間的相關性分析

表3 主要成分二維分子描述符與中藥藥性之間的相關性分析

續(xù)表

表4 主要成分三維分子描述符與中藥藥性之間的相關性分析

表5 主要成分半經(jīng)驗量化分子描述符與中藥藥性之間的相關性分析

根據(jù)前人研究的成果，結(jié)合我們的研究結(jié)果，可以推斷，部分中藥不同藥性的作用，可能是通過其不同分子結(jié)構(gòu)（分子能量狀態(tài)）的主要有效成分，進入機體以后，由藥物分子的物理內(nèi)能以及藥物與生物機體之間的相互作用，影響對機體能量代謝而產(chǎn)生的綜合效應。

除了分子能量狀態(tài)之外，中藥藥性可能與其主要成分的脂水分配系數(shù)、拓撲分子極性表面積，碳、氫、氧原子數(shù)，氫鍵供體數(shù)目等分子結(jié)構(gòu)描述符具有一定的關聯(lián)關系。

4 結(jié)論

通過研究結(jié)果可以推斷，中藥藥性確實與中藥主要成分的分子結(jié)構(gòu)具有一定的關聯(lián)關系，其中最重要的分子描述符是三維分子描述符中的分子能量，提示中藥藥性作用可能與中藥主要有效成分的分子能量狀態(tài)以及對機體能量代謝的影響而產(chǎn)生的綜合效應有關。

1王振國,王鵬,李峰,等.中藥藥性理論現(xiàn)代研究回顧與展望.山東中醫(yī)藥大學學報,2008,32(2):94-97.

2盛良.論化學元素具有寒熱溫涼四性-論中藥四性與現(xiàn)代化學的統(tǒng)一.中國中醫(yī)基礎醫(yī)學雜志,2003,9(1):14-16.

3盛良.論中藥礦物藥四性與無機化學的結(jié)合—二論中藥四性與現(xiàn)代化學的統(tǒng)一.中國中醫(yī)基礎醫(yī)學雜志,2004,10(3):24-26.

4蔣俊,賈曉斌,陳斌,等.基于組分結(jié)構(gòu)假說構(gòu)建中藥藥性矢量結(jié)構(gòu)方程.中草藥,2012,43(3):422-426.

5李和光.中藥寒熱藥性物質(zhì)基礎紫外光譜指紋圖譜數(shù)字化表征體系的構(gòu)建.濟南：山東中醫(yī)藥大學碩士學位論文,2007:22-24.

6王曉燕,李峰.60種中藥脂類成分GC-MS圖譜與寒熱藥性關系的統(tǒng)計識別模式研究.遼寧中醫(yī)雜志,2015,42(6):1303-1305.

7常美佳,羅盛,楊旭曙,等.2,4-噻唑烷二酮類醛糖還原酶抑制劑的三維定量結(jié)構(gòu)-活性關系及分子對接研究.南京師大學報(自然科學版),2015,38(3):31-40.

8Yu J,Li B,Lee S,et al.Relationship between physical properties and chemical structures of poly(ethylene terephthalate-co-ethylene isophthalate).J Appl Polym Sci,1999,73(7):1191-1195.

9Kotritzky A R,Maran U,Lobanov V S,et al.Structurally diverse quantitative structure-property relationship correlations of technologically relevant physical properties.J Chem Inform Comput Sci,2000,40(1): 1-18.

10付先軍,王鵬,王振國.從中藥“性-構(gòu)關系”探索構(gòu)建寒熱藥性成分要素表征體系的研究構(gòu)想.世界科學技術-中醫(yī)藥現(xiàn)代化,2011,13 (5):919-924.

11中華人民共和國藥典委員會.中華人民共和國藥典2015版(一部)北京:中國醫(yī)藥出版社,2015.

12 Gasteiger J,Engle T.化學信息學教程.北京：化學工業(yè)出版社,2005.

13 Cao L,Koo S,Sun J,et al.Structure descriptor for surface passivation in the simulation of atomistic models.Sci China Inform Sci,2017,60：032103.

14李秉軻,叢湧,田之悅,等.基于分子描述符和機器學習方法預測和虛擬篩選MMP-13對MMP-1的選擇性抑制劑.物理化學學報, 2014,30(1):171-182.

15李天瑞,楊寧,馬駿.廣義關聯(lián)規(guī)則基及其挖掘算法.西南交通大學學報,2004,39(1):86-89.

16 Leon C,Biscarri F,Monedero I,et al.Variability and Trend-Based Generalized Rule Induction Model to NTL Detection in Power Companies. IEEE T Power Syst,2011,26(4):1798-1807.

17侯燦.“八綱”病理生理學基礎初步探討.中醫(yī)雜志,1964(12):46-47.

18桑木崇秀.寒熱藥與證候.にほんとぅよぅぃがくざっし,1978,29 (1):7.

19肖小河,王伽伯,趙艷玲,等.藥性熱力學觀及實踐.中國中藥雜志, 2010,35(16):2207-2213.

20趙艷玲,肖小河,劉義.熱力學理論在中醫(yī)藥研究中的應用與發(fā)展.解放軍藥學學報,2005,21(4):291-293.

21樊冬麗,廖慶文,鄢丹,等.基于生物熱力學表達的麻黃湯和麻杏石甘湯的寒熱藥性比較.中國中藥雜志,2007,32(5):421-424.

Study on Relationship between Nature and Multidimensional Structure Descriptor of Main Compounds from Chinese Medicinal Herbs

Fu Xianjun,Wang Zhenguo,Li Xuebo
(Key Laboratory of Classical Theory of Traditional Chinese Medicine of the Ministry of Education of China,Institute of Literature in Chinese Medicine,Taishan Scholar Studio of Specialty of History and Literature of Traditional Chinese Medicine,Shandong University of Traditional Chinese Medicine,Jinan 250355,China)

This study was aimed to analyze the cold and heat nature of Chinese herbal medicine as well as its relation with the major constituent part in order to illustrate the scientific connotation of the cold and heat nature of Chinese herbal medicine from the level of material basis.There were a total of 334 kinds of Chinese herbal medicine with 280 effective ingredients recorded in the first volume of the Chinese Pharmacopoeia.The molecular descriptors of each ingredient(including one-dimensional,two-dimensional,three-dimensional,and semi-empirical calculated molecular descriptors)were calculated.The generalized rule induction(GRI)algorithm was used to mine molecular descriptors with strong correlation to the cold or heat nature of the screened Chinese herbal medicine.The results showed that molecular descriptors of molecular energy state of main components had close relationship with the nature properties of herbs. Compounds with high-energy condition were most from herbs with heat nature;while compounds with low-energy condition were most from herbs with cold nature.In addition to the energy condition of the molecule,topological molecular water partition coefficient,lipid drug property but also the main component of the polar surface area,carbon, hydrogen,oxygen atoms,the number of hydrogen bond donors and other molecular structure descriptors had certain association relationships.It was concluded that the nature of Chinese herbal medicine had certain association relationships with its main ingredients.Among them,the most important molecular descriptors of molecular energy are three-dimensional molecular descriptors.It suggested that the comprehensive effect of Chinese herbal medicine may be relevant to the energy condition of the molecule in main active components of Chinese herbal medicine and the effect on the energy metabolism of the body.

Cold-heat nature of traditional Chinese medicine,structural descriptors,association rules mining,“Nature-Structure Relationship”of TCM

10.11842/wst.2017.04.002

R286.0

A

（責任編輯：王涵，責任譯審：王晶）

2017-03-14

修回日期：2017-03-27

＊國家自然科學基金委面上項目（81473369）：基于“性-構(gòu)”關系的中藥成分寒熱藥性評價，負責人：付先軍；科學技術部國家重點基礎研究發(fā)展計劃“973計劃”課題（2007CB512601）：基于四性的中藥性-效-物質(zhì)關系研究，負責人：王振國。

＊＊通訊作者：付先軍，副教授，主要研究方向：中藥藥性理論指導下的海洋中藥研究；王振國，教授，主要研究方向：中藥藥性理論研究。