董宏偉，涂清波

(南京中醫(yī)藥大學翰林學院，江蘇 泰州 225300)

肉桂(Cinnamomum cassia Presl)是樟科樟屬植物的一種常綠樹，具有較高的醫(yī)用和食用價值。《本草綱目》記載肉桂具有“補下焦不足，治沉寒痼冷之病，滲泄止渴，去榮衛(wèi)中風寒，表虛自汗”等功效。我國是肉桂原產(chǎn)地和種植大國，國內(nèi)資源豐富及產(chǎn)量巨大。肉桂精油是從桂皮、桂枝、桂葉等部位提取的芳香味油狀液體?，F(xiàn)代藥理表明[1]，其具有抗氧化、抗炎、抗菌和抗癌等多種生物活性。中藥治療疾病存在多成分、多靶點、多通路同時作用的特性[2]。氧化應激是一種生物體內(nèi)部產(chǎn)生的自由基和活性氧物質(zhì)引起的一系列有害生理過程[3]，會導致脂質(zhì)過氧化、蛋白質(zhì)損傷、DNA斷裂等多種細胞損傷，進而引起許多疾病，如心血管疾病、腫瘤等。肉桂精油的抗氧化作用明確，但其具體的作用機制研究還不夠清楚。網(wǎng)絡藥理學適用于中藥復雜成分的解析和開發(fā)[4]，可對肉桂精油抗氧化的作用機制進行系統(tǒng)全面的分析。本研究運用氣質(zhì)聯(lián)用和網(wǎng)絡藥理學對肉桂精油活性成分、作用靶點、信號通路等進行系統(tǒng)分析，初步探究肉桂精油的抗氧化機制，為其進一步藥理研究和產(chǎn)品開發(fā)提供科學依據(jù)。

1 材料與方法

1.1 實驗試劑

桂皮粉為食品級，購自泰州市邦藍食品有限公司；甲醇為色譜級，購自國藥集團化學試劑有限公司。

1.2 肉桂精油化學成分分析

根據(jù)劉彤[5]等提取方法稍加修改，稱取 100 g 肉桂粉置于 1 L 燒瓶中，加入 600 mL 蒸餾水回流提取 2 h，收集得到肉桂精油。運用氣質(zhì)聯(lián)用技術(GC-MS)[6]。簡單地說，用色譜甲醇稀釋肉桂精油，經(jīng) 0.22 μm 微孔濾膜過濾后進氣相色譜，色譜柱為HP-INNOWax polyethylene glycol(30 m×0.32 mm×0.25 μm)，載氣為氦氣，流速為 2.0 mL/min。設置 40 ℃ 開始保溫 2 min，5 ℃/min 升高至 240 ℃，保持 18 min，注射口溫度為 250 ℃，分流比為20∶1。電離源和傳輸線的溫度分別為230 ℃和250 ℃，離子能量為 70 eV，質(zhì)量范圍為40～400 Da。

1.3 肉桂精油活性成分對應靶點的獲取

利用TCMSP數(shù)據(jù)庫檢索肉桂精油主要化學成分的對應靶點[7]，通過Pubchem數(shù)據(jù)庫，獲得主要成分的SMILE Name，并導入SWISS Taget Prediction數(shù)據(jù)庫中，獲取單個主要成分的Gene名。

1.4 肉桂精油抗氧化相關靶點的收集

在Genecards疾病數(shù)據(jù)庫中，以“Anti-oxidant”為關鍵詞，在Genecards數(shù)據(jù)庫中檢索并獲取對應靶點[8]，從而獲得與抗氧化相關的靶點信息，去除重復值，匯總抗氧化的相關靶點。將藥物和疾病靶基因取交集，獲得兩者的共同靶點基因，繪制Venn圖。

1.5 蛋白質(zhì)相互作用網(wǎng)絡的構(gòu)建及可視化分析

將共同靶點數(shù)據(jù)集導入String數(shù)據(jù)庫[9]，得到蛋白質(zhì)相互作用網(wǎng)絡，通過運用Cytoscape軟件構(gòu)建“藥材-成分-靶點”可視化、多元素、多水平間的作用關系。網(wǎng)絡中各節(jié)點分別代表活性成分和關鍵靶點基因；網(wǎng)絡中邊用來連接活性成分與關鍵靶點基因；連接到網(wǎng)絡的節(jié)點以度值為單位進行表示。某節(jié)點與其他節(jié)點連接的邊數(shù)越多，則說明該節(jié)點在網(wǎng)絡中扮演的角色越重要，度值也越大。

1.6 GO與KEGG富集分析

對篩選得到的化合物對應的潛在靶點進行GO富集分析和KEGG通路富集分析[10]，GO富集分析選擇生物過程(BP)、分子功能(MF)和細胞組成(CC)3個模塊，通路分析選擇KEGG，設定p<0.05，選取物種為人類，通過KEGG分析可以得到肉桂抗氧化的潛在作用通路。

1.7 “活性成分-核心靶點-通路”網(wǎng)絡圖的構(gòu)建

將涉及的活性成分與通路及其相關靶點導入Cytoscape軟件得到“活性成分-核心靶點-通路”網(wǎng)絡[11]，計算網(wǎng)絡中成分節(jié)點進行拓撲分析，反向篩選肉桂精油抗氧化的有效成分及關鍵基因。

2 結(jié)果與分析

2.1 肉桂精油化學成分分析

經(jīng)GC-MS和NIST標準數(shù)據(jù)庫對比分析，肉桂精油的化學成分包含36個化學物質(zhì)(α-蒎烯、莰烯、1，2-二甲基-3-(1-甲基乙烯基)環(huán)戊烷、3，7，7-三甲基雙環(huán)[4.1.0]庚烷、鄰薄荷-8-烯、3-蒈烯、1，4-桉葉素、o-環(huán)己烯、D-檸檬、桉樹醇、苯甲醇、雙戊烯、1，3，3-三甲基雙環(huán)[2.2.1]庚烷-2-氧、2-(4-甲基苯基)丙-2-醇、反式-1-乙烯基-2-甲基環(huán)己烷、水楊酸甲酯、反式肉桂醛、9-甲氧基雙環(huán)[6.1.0]壬-2，4，6-三烯、巴倫西亞橘烯、2，3，4，7，8，8a-六氫-3，6，8，8-四甲基-1H-3a，7-亞甲基薁-5-甲醇、苯甲酸芐酯、N-(1-甲基乙基)-苯甲胺、鄰苯二甲酸二異丁酯、鄰苯二甲酸二丁酯、肉桂酸芐酯、[2-甲基-2-(4-甲基-3-戊烯基)環(huán)丙]-甲醇、α-紫羅蘭酮、甲位異丁酸松油酯、鳶尾酮、2，4，5-三甲基苯乙醛、(+)-3-蒈烯、環(huán)氧石竹烯、1R，4R-對-甲基-2，8-二烯-1-過氧化氫、1，2-二甲基-1，2-乙二醇、順-7-十四烯醛、(E)-6-十六碳烯-4-炔)，其中峰面積比大于1%的主要成分有5個，分別是反式肉桂醛(74.42%)、鄰苯二甲酸二丁酯(8.28%)、水楊酸甲酯(4.51%)和鄰苯二甲酸二異丁酯(4.05%)和鳶尾酮(1%)。

2.2 活性成分靶點網(wǎng)絡的構(gòu)建

肉桂精油36個化學成分對應的潛在藥物靶點為185個。構(gòu)建的活性成分靶點網(wǎng)絡如圖1所示，共有108個節(jié)點(91個基因靶點，1個中藥名稱節(jié)點和16個成分節(jié)點)150條邊。根據(jù)自由度對活性成分進行排序，主要活性成分分別是甲位異丁酸松油酯、鄰苯二甲酸二異丁酯、肉桂酸芐酯、α-紫羅蘭酮、反式肉桂醛。

圖1 藥物-成分-靶點相互作用圖(其中，黃色代表活性成分，藍色代表靶點)

2.3 抗氧化靶點和核心靶點的獲取

將Genecards數(shù)據(jù)庫收集到的靶點合并去除重復后，得到肉桂精油成分靶點和抗氧化相關靶點的維恩圖如圖2所示，抗氧化的靶點有2038個，肉桂精油中活性成分對應的靶點和抗氧化對應的靶點的交集靶點有91個。

圖2 肉桂精油成分靶點和抗氧化相關靶點的維恩圖

2.4 蛋白質(zhì)相互作用網(wǎng)絡的構(gòu)建

將交集靶點導入STRING據(jù)庫后，構(gòu)建蛋白質(zhì)相互作用網(wǎng)絡如圖3a所示，其中抗氧化交集靶點網(wǎng)絡結(jié)構(gòu)圖中有91個節(jié)點，591條邊線組成。通過Cytoscape軟件中的Betweenness、Degree和Closeness算法依次取前50%篩選出肉桂抗氧化的核心靶點，如圖3b所示，肉桂活性成分中抗氧化的核心靶點有10個，分別為CASP3、PTGS2、EGFR、HIF1A、STAT3、TLR4、ESR1、APP、RELA、PPARA。

圖3 a：肉桂精油抗氧化核心靶點PPI網(wǎng)絡；b：肉桂精油活性成分靶點與靶點互相關系

2.5 關鍵靶點的GO功能注釋和KEGG通路富集分析

利用DAVID數(shù)據(jù)庫，對10個肉桂精油活性成分抗氧化靶點基因進行GO和KEGG功能富集分析。GO注釋如圖4所示，富集BP、CC、MF、的富集條目一共有546個。其中BP的條目主要有379個，富集主要部位在細胞膜上，主要涉及細菌的免疫原性和致病性。CC的主要條目有60條，主要涉及細胞的代謝和信號轉(zhuǎn)導，蛋白質(zhì)合成與折疊，物質(zhì)運輸。MF主要條目有107條，主要涉及酶催化的速率、底物選擇性以及酶的穩(wěn)定性。經(jīng)KEGG通路富集分析，得到肉桂精油抗氧化的潛在作用通路共獲得132條通路，按P值按照降序排列，排名前十的信號通路制作氣泡圖，如圖5所示。氣泡點的大小表示相應功能通路中靶基因的數(shù)量，富集度表示屬于該通路中所有被注釋基因的靶基因數(shù)量的比率。P值越低、顏色越紅代表KEGG富集程度更高。表明TNF、FoxO、C-type lectin receptor、PD-L1等信號通路在肉桂精油抗氧化中起主要作用，糖尿病性心肌病通路被認為是一組關鍵通路。肉桂精油抗氧化的關鍵靶點富集在：多巴胺能突觸、糖尿病性心肌病、細胞凋亡、腫瘤壞死、生長激素缺乏癥等相關途徑和通路中。

圖4 肉桂精油抗氧化核心靶點的GO富集分析

圖5 肉桂精油抗氧化核心靶點的KEGG富集分析

2.6 成分-靶點-通路的網(wǎng)絡構(gòu)建

將核心靶點及其相關的活性成分與前20條通路及其相關靶點導入Cytoscape軟件。刪除16條重復得到“活性成分-核心靶點-通路”網(wǎng)絡如圖6所示，該網(wǎng)絡包括192個節(jié)點(其中包括化合物16個，91個靶點，10個通路網(wǎng)絡)267條邊。從通路角度來看，連接度較高的為hsa04728、hsa05415。從靶點角度分析篩選出6個DC≥10的靶點，分別為MAPK11、PIK3CA、CYP19A1、PIK3CA、MAPK14、MAPK10。從化合物角度來看，共篩選出5個主要化合物，分別為甲位異丁酸松油酯、鄰苯二甲酸二異丁酯、肉桂酸芐酯、α-紫羅蘭酮、反式肉桂醛。肉桂精油中存在多成分作用于一個靶點，成一個靶點作用于多個成分，也存在一條通路多個靶點，或不同通路同屬于一個靶點的現(xiàn)象，表明肉桂精油通過多成分、多靶點、多通路方式的抗氧化作用。

圖6 成分-核心靶點-通路網(wǎng)絡

3 結(jié)論

植物精油因其良好的生物活性、安全性、滲透性和生物降解性，被廣泛應用于醫(yī)藥化工、食品工業(yè)和化妝品行業(yè)中。本文運用氣相色譜-質(zhì)譜聯(lián)用和網(wǎng)絡藥理學的方法，發(fā)現(xiàn)肉桂精油主要成分為36個，篩選出對應靶點185個，疾病靶點2038個，共有靶點91個，涉及546個基因功能條目和132條通路。經(jīng)“活性成分-核心靶點-通路”網(wǎng)絡分析，篩選出5種核心成分、6個核心靶點以及2個連接度較高的信號通路，說明肉桂精油可通過多種成分、多個靶點以及不同通路的來實現(xiàn)其抗氧化作用，初步探討了肉桂精油的抗氧化作用機制，為后續(xù)藥理研究和產(chǎn)品開發(fā)提供科學依據(jù)。