潘丹陽，何敬愉，陳慧芳，陳英杰，劉孟華，龔盛昭，

（1. 廣州環(huán)亞化妝品科技股份有限公司，廣東 廣州 510663；2. 南方醫(yī)科大學 藥學院，廣東 廣州 510515；3. 廣東輕工職業(yè)技術(shù)學院，廣東 廣州 510300）

脫發(fā)（Alopecia） 或稱禿發(fā)，是一種以毛發(fā)減少為特征的皮膚疾?。?］。它的種類較多，如先天性脫發(fā)、斑禿以及雄激素源性脫發(fā)等，但病因尚不完全清楚。一般認為與遺傳、內(nèi)分泌功能、感染病灶、自身免疫功能、精神因素以及營養(yǎng)狀況等都有一定關(guān)系［2］。生姜為姜科植物姜（Zingiber officinale Rosc.）的新鮮根莖，具有解表散寒，溫中止嘔，化痰止咳，解魚蟹毒的功效［3］?，F(xiàn)代藥理學研究發(fā)現(xiàn)，生姜具有抗菌、抗腫瘤、抗炎、解熱、抗氧化、降血糖、保肝、利尿、降膽固醇等作用［4］。揮發(fā)油是生姜中的主要成分，不僅應用到醫(yī)藥產(chǎn)品中，還應用到功能食品、化妝品和香料中?；瘜W研究表明生姜揮發(fā)油主要含有萜類化合物，如單萜烯類、單萜烯類氧化物、倍半萜烯類、倍半萜烯類氧化物，還含有酚類和鏈烴類化物［5-6］。研究表明：生姜提取物可增加人毛囊真皮乳頭細胞，從而改善脫發(fā)癥狀［7］。亦有報道證實生姜提取物具有明顯的促毛發(fā)生長作用，其對毛發(fā)的生長速度影響呈現(xiàn)劑量依賴性［8］，但其作用機制尚未明確。本文旨在明確生姜揮發(fā)油化學成分的基礎(chǔ)上，采用網(wǎng)絡藥理學方法探究其緩解斑禿癥狀的可能分子機制，為其應用提供科學數(shù)據(jù)。

1 材料與方法

1.1 材料

利用中藥系統(tǒng)藥理學分析平臺TCMSP和BATMAN-TCM平臺得到生姜揮發(fā)油相關(guān)的化學成分。

1.2 方法

1.2.1 蛋白靶點的收集、篩選和處理

利 用GeneCard數(shù) 據(jù) 庫（https：//www.genecards.org/）對“alopecia”（Score≥10）進行檢索，對檢索結(jié)果進行保存。另外，UniProt數(shù)據(jù)庫（https：//sparql.uniprot.org/）以“alopecia+reviewed +Homosapiens”進行檢索，并對檢索結(jié)果進行保存。利用Venn分析工具，將病理靶點進行疊加，除去重復靶點。在RCSB Protein Data Bank（RCSB PDB）數(shù)據(jù)庫（http：//www.rcsb.org/）篩選靶蛋白，并根據(jù)是否存在蛋白-配體復合物、晶體結(jié)構(gòu)以及分辨率，最終確定了57個靶蛋白。在進行分子對接前，應用Sybyl7.3軟件對靶蛋白進行修飾處理，包括去除原配體、去除水分子、加氫以及對氨基酸殘基和力場進行優(yōu)化等處理。

1.2.2 分子對接

使用Sybyl7.3軟件中的Surflex-Dock插件將生姜揮發(fā)性成分和人類脫發(fā)癥相關(guān)靶蛋白進行分子半柔性對接。

1.2.3 分子網(wǎng)絡的構(gòu)建

為了更準確地找到有效成分，識別分子的主要靶標，利用Cytoscape 3.7.1 軟件進行化合物-脫發(fā)癥靶蛋白網(wǎng)絡圖的繪制。首先，根據(jù)上述分子對接結(jié)果，以靶蛋白與原配體的對接分數(shù)（Total Scores）為閾值，打分高于此閾值的化合物即認為具有較好的活性。本次實驗選擇對接分數(shù)大于或等于6的分子-靶蛋白創(chuàng)建一個Microsoft Excel文件，將Microsoft Excel作為輸入源，導入到Cytoscape 3.7.1軟件中，構(gòu)建分子-靶蛋白網(wǎng)絡圖，并利用Cytoscape 3.7.1中的 “Network analyzer”對生姜分子-靶蛋白網(wǎng)絡特征參數(shù)進行了分析。在構(gòu)建的網(wǎng)絡中，化合物與靶蛋白用節(jié)點表示，當化合物與靶蛋白有相互作用時用邊緣連接。通過得到的拓撲異構(gòu)學參數(shù)來預測可能的活性成分群和關(guān)鍵的治療脫發(fā)靶蛋白。

1.2.4 蛋白-蛋白相互作用（PPI）關(guān)系的構(gòu)建

為了更好地探索生姜揮發(fā)性成分改善脫發(fā)癥的作用機制，利用String數(shù)據(jù)庫建立了靶點與靶點之間的相互作用關(guān)系，將與脫發(fā)癥靶點的基因輸入到String數(shù)據(jù)庫（https：//string-db.org/），限定物種設(shè)置為“Homo Sapiens”，以置信度（0.4）為條件篩選，獲取蛋白相互作用關(guān)系。分析結(jié)果保存為Tsv格式文件，導入 Cytoscape 3.7.1軟件，構(gòu)建PPI網(wǎng)絡圖，并利用Cytoscape 3.7.1中的 “Network analyzer”對生姜分子-靶蛋白網(wǎng)絡特征參數(shù)進行了分析。結(jié)合節(jié)點度值大小，構(gòu)建PPI網(wǎng)絡。

1.2.5 靶點通路分析

借助于DAVID 6.8數(shù)據(jù)庫（https：//david.ncifcrf.gov/）對疾病相關(guān)的藥物活性成分靶點進行GO 功能富集分析與KEGG 通路富集分析。GO功能富集分析結(jié)果根據(jù)錯誤發(fā)現(xiàn)率（FDR）＜0.05標準進行歸類展示。KEGG通路富集分析獲得的信號通路結(jié)果根據(jù)FDR＜0.05標準進行篩選、排序、展示。

2 結(jié)果

2.1 揮發(fā)油成分分析

通過對TCMSP平臺和BATMAN-TCM平臺進行檢索，除去烷烴類化合物后，共獲取生姜揮發(fā)性成分105個。其化合物名稱見表1。

表1 生姜揮發(fā)油中化合物信息

2.2 蛋白靶點的分析

GeneCard數(shù)據(jù)庫檢索到與人類脫發(fā)癥相關(guān)的基因有182個（Score≥10），通過Uniprot數(shù)據(jù)庫檢索到的人類脫發(fā)癥相關(guān)的基因有46個，除去重復基因，共獲取196個與人類脫發(fā)癥相關(guān)的基因。在RCSB PDB數(shù)據(jù)庫篩選蛋白靶點，并根據(jù)是否存在蛋白-配體復合物、晶體結(jié)構(gòu)以及分辨率，最終確定了57個靶蛋白。

2.3 分子-靶點相互作用網(wǎng)絡的構(gòu)建

使用Sybyl7.3軟件中的Surflex-Dock 插件將生姜揮發(fā)油中105個化合物和57個人類脫發(fā)癥相關(guān)蛋白靶點進行分子半柔性對接，如表2所示。根據(jù)設(shè)定的閾值，有48個化合物對48個靶蛋白表現(xiàn)出較好的活性，并構(gòu)建分子-靶蛋白網(wǎng)絡關(guān)系（見圖1，綠色表示靶點；紫色表示化合物；節(jié)點的大小與節(jié)點的度呈正比例關(guān)系）。構(gòu)建的化合物-靶點網(wǎng)絡包含節(jié)點數(shù)為96個，邊數(shù)為470條，平均相鄰節(jié)點數(shù)9.771個，最短路徑為9 120，特征路徑長度2.600，網(wǎng)絡密度0.103，網(wǎng)絡中心度0.293，網(wǎng)絡異質(zhì)性0.825。分析獲得生物網(wǎng)絡和節(jié)點拓撲學特征，挖掘潛在的生物學信息。在網(wǎng)絡中，一個節(jié)點的度值（Degree）表示網(wǎng)絡中和節(jié)點相連的路線的條數(shù)。根據(jù)拓撲學性質(zhì)篩選出的度值較大的節(jié)點進行分析，這些連接化合物或靶點較多的節(jié)點在整個網(wǎng)絡中起到樞紐作用，可能為關(guān)鍵的化合物或者靶點。有36.7%的化合物的作用靶點≥10個，作用靶點≥20個的化合物有10個。從化合物角度來看，度值排名前5位 的 是6-Methylgingediacetate（Com10）、10-Gingediol （Com73）、Gingerenone B（Com72）、10-Gingerol（Com29）、Gingerdione（Com87）。從蛋白靶點角度來看，度值排名前5位為FLG、SLITRK1、IL1B、CASR、ANTR1。

圖1 生姜揮發(fā)油-靶點網(wǎng)絡

表2 用于分子對接的靶點信息

2.4 蛋白-蛋白相互作用（PPI）網(wǎng)絡的構(gòu)建

構(gòu)建的蛋白-蛋白網(wǎng)絡包含節(jié)點數(shù)為40個，邊數(shù)為269條，其中節(jié)點表示蛋白，每條邊則表示蛋白與蛋白之間的相互作用關(guān)系。構(gòu)建的PPI網(wǎng)絡（見圖2，藍色表示Degree≥20的蛋白；粉色表示10≤Degree ＜ 20的蛋白；綠色表示Degree ＜ 10的蛋白；節(jié)點的大小與節(jié)點的度呈正比例關(guān)系）的平均相鄰節(jié)點數(shù)13.45個，最短路徑為1560，特征路徑長度1.822，網(wǎng)絡密度0. 345，網(wǎng)絡中心度0.420，網(wǎng)絡異質(zhì)性0.652。Degree≥20的節(jié)點有12個，包括TNF、IL6、ALB、IL10、IL17A、IL2、PTPRC、IFNG、IL4、ICAM1、IL13、CXCL10。其中度值前3名為TNF、IL6、ALB，分別能與29、28、27個蛋白發(fā)生相互作用。

圖2 生姜揮發(fā)油相關(guān)靶點的PPI網(wǎng)絡

2.5 GO功能富集分析

對生姜揮發(fā)油中化合物有作用的48個蛋白進行GO生物學過程富集分析（見表3）。根據(jù)錯誤發(fā)現(xiàn)率（FDR）＜0.05確定了687個生物過程（Biological Process，BP）條目，21個分子功能（Molecular Function，MF）條目，14個細胞成分（Cellular Component，CC）條目，分別占95.15%、2.91%、1.94%，表3根 據(jù)FDR列出前20個BP條目，前5個MF條目，前5個CC條目。

表3 PPI網(wǎng)絡中的GO條目

2.6 KEGG通路富集分析

通過KEGG通路富集分析，對生姜揮發(fā)油中化合物有作用的48個蛋白在信號通路中的作用進行了研究，得到118條信號通路，根據(jù)FDR＜ 0.01篩選出20條（見表4）與炎癥及免疫相關(guān)的信號通路，如細胞因子-細胞因子受體相互作 用（Cytokine-cytokine receptor interaction）、IL-17信 號 通 路（IL-17 signaling pathway）、T細胞受體信號通路（T cell receptor signaling pathway）、Janus激酶/信號轉(zhuǎn)導與轉(zhuǎn)錄激活子信號通路（Jak-STAT signaling pathway）、Th17細胞分化（Th17 cell differentiation）、磷脂酰肌醇3激酶-絲氨酸/蘇氨酸激酶信號通路（PI3KAkt signaling pathway）。

表4 PPI網(wǎng)絡中的KEGG通路

3 討論

中藥是一個由多類別化合物組成的復雜體系，同時通過多分子、多靶點、多途徑作用于疾病，而起到一定的治療效果。脫發(fā)是臨床上常見的皮膚病，由諸多因素引起，存在多種發(fā)病機制［1］。本研究通過分子對接手段分析活性化合物與靶點蛋白間的相互作用，網(wǎng)絡藥理學構(gòu)建化合物-靶點網(wǎng)絡和PPI網(wǎng)絡，探究生姜揮發(fā)油改善和治療脫發(fā)癥的作用機制。當前研究表明生姜提取物具有明顯的促毛發(fā)生長作用，其對毛發(fā)的生長速度影響呈現(xiàn)劑量依賴性［8］。但有一報道提示，涂抹1 mg/ml高濃度的6-姜辣素可顯著減少C57BL/6小鼠的毛囊數(shù)量［9］。生姜揮發(fā)油的使用量是一個影響其治療功效的重要因素。

在構(gòu)建的PPI網(wǎng)絡中白細胞介素（IL）和腫瘤壞死因子-α（TNF-α）在脫發(fā)癥中是重要的靶點。IL是在傳遞信息，激活與調(diào)節(jié)免疫細胞，介導T、B細胞活化、增殖與分化及在炎癥反應中起重要作用［10］。TNF-α是一種重要的促炎因子，通過引起機體失衡和炎癥反應，破壞體內(nèi)穩(wěn)態(tài)維持過程［11］。生姜水提物、6-姜辣素和山姜醇有抗炎作用，通過抑制IL-1β、TNF-α、 IL-6和 IL-8的產(chǎn)生，對UVB誘導的人角質(zhì)形成細胞和小鼠皮膚損傷有保護作用［12］。6-姜烯酚以濃度依賴的方式抑制細胞中TNF-α、IL-1β和IL-6的產(chǎn)生［13］。另外，8-姜辣素、10-姜辣素、8-姜烯酚和10-姜烯酚可以抑制TNF-α的產(chǎn)生［14］。細胞間粘附分子-1（ICAM-1）是一種屬于免疫球蛋白超家族的細胞表面糖蛋白。ICAM-1在許多類型的細胞中低水平表達，其表達主要在NF-κB轉(zhuǎn)錄水平上調(diào)以響應促炎性細胞因子［15］。IFN-γ是由特定抗原或有絲分裂原激活的淋巴細胞產(chǎn)生的，具有重要的免疫調(diào)節(jié)功能［16］。全基因組關(guān)聯(lián)研究表明，毛囊的喪失是由先天免疫系統(tǒng)和適應性免疫系統(tǒng)共同介導的，Th1細胞因子和CD8+T細胞是主要的效應因子［17-18］。脫發(fā)皮損中IFNγ水平的升高可誘導趨化因子CXCL9和CXCL10，這些趨化因子在放大毛囊周圍的T細胞反應中起重要作用，從而使其靶向應對脫發(fā)的特異性自身抗原［19-20］。目前，對生姜揮發(fā)油中免疫和抗炎成分研究主要集中于姜酚和姜烯酚，對于本研究挖掘到的潛在活性化合物6-Methylgingediacetate、10-Gingediol、Gingerenone B和Gingerdione，在免疫和炎癥方面的研究非常少。這也提示生姜揮發(fā)油中姜酮和姜烯酮是兩類重要的活性組分。從分子-靶蛋白網(wǎng)絡和PPI網(wǎng)絡結(jié)果分析，對于脫發(fā)癥改善和治療，生姜中存在多分子與單靶蛋白作用、單分子與多靶蛋白作用的現(xiàn)象，符合中藥多組分、多靶點協(xié)同作用的特征。

GO功能富集分析發(fā)現(xiàn)，生姜揮發(fā)油在改善脫發(fā)時的基因功能主要體現(xiàn)在對生物過程、分子功能、細胞成分等過程，特別是對免疫系統(tǒng)調(diào)控方面。對基因的KEGG通路富集分析中，大部分基因富集到了與炎癥反應和免疫反應相關(guān)的通路上，主要有細胞因子受體相互作用、IL-17信號通路、T細胞受體信號通路、Jak-STAT信號通路、PI3K-Akt信號通路，分別富集了18，8，8，9，11個基因。細胞因子是可溶的細胞外蛋白或糖蛋白，是細胞的重要細胞間調(diào)節(jié)因子和動員劑，參與固有和適應性炎癥宿主防御、細胞生長、分化、細胞死亡、血管生成以及旨在恢復內(nèi)穩(wěn)態(tài)的發(fā)育和修復過程［21］。白細胞介素17（IL-17）家族通過其相應的受體發(fā)出信號，激活包括NF-κB、MAPKs和C/EBPs在內(nèi)的下游通路，誘導抗菌肽、細胞因子和趨化因子的表達，在急性和慢性炎癥反應中起著至關(guān)重要的作用［22］。生姜提取物可以抑制絲裂原激活蛋白激酶（MAPK）、細胞外信號調(diào)節(jié)激酶1和2（ERK1/2）、p38和c-Jun N-末端激酶（JNK）的磷酸化以及NF-κB的激活，通過MAPK和NF-κB信號通路抑制炎癥介質(zhì)基因的轉(zhuǎn)錄而發(fā)揮 抗炎 作 用［23］。T細 胞 受體（T Cell Receptor，TCR）在T細胞的功能和免疫突觸的形成中起著關(guān)鍵作用。它提供T細胞和抗原呈遞細胞之間的聯(lián)系。TCRs的激活促進了許多信號級聯(lián)，這些級聯(lián)最終通過調(diào)節(jié)細胞因子的產(chǎn)生、細胞存活、增殖和分化來決定細胞的命運［24］。生姜提取物能顯著抑制Th2細胞介導的免疫反應，這表現(xiàn)為IL-4和IL-5的mRNA表達水平降低［25］。6-姜辣素抑制了Th1細胞和Th2細胞從初始T細胞的體外分化，還影響T細胞受體介導的信號轉(zhuǎn)導［26］。Jak-STAT信號通路和PI3K-Akt信號通路是脫發(fā)癥狀中不可忽視的；Jak/STAT途徑是為數(shù)不多的多效性級聯(lián)之一，是一系列細胞因子和生長因子的主要信號機制，用于在動物中傳遞大量的發(fā)育和穩(wěn)態(tài)信號。除了激活STATs外，Jaks還介導其他分子的募集，如MAPK、PI3K等［27］。在PI3K-Akt 信號通路中，PI3K催化細胞膜上磷脂酰肌醇-3，4，5-三磷酸（PIP3）的生成。PIP3反過來充當?shù)诙€信使，幫助激活Akt。激活后的Akt可以通過磷酸化參與凋亡、蛋白質(zhì)合成、代謝和細胞周期的底物來控制關(guān)鍵的細胞過程［28］。研究表明姜烯酮A能夠抑制Jak2；與6-姜烯酚、姜黃素、四氫姜黃素比較，它還能抑制STAT3的磷酸化［29］。此外，6-姜烯酚通過抑制NF-κBp65的磷酸化和核轉(zhuǎn)位來抑制LPS誘導的NF-κB活化［13］。這些說明由炎癥及自身免疫引起的脫發(fā)癥，生姜揮發(fā)油可以通過對炎癥和免疫等多個生物途徑信號通路進行調(diào)節(jié)，發(fā)揮改善和治療作用。

綜上所述，本研究基于網(wǎng)絡藥理學與分子對接技術(shù)，分析了生姜揮發(fā)油的多成分、多靶點以及多通路的復雜網(wǎng)狀關(guān)系，揭示了生姜揮發(fā)油改善脫發(fā)癥狀的潛在分子機制。為了對篩選出來的相關(guān)基因及蛋白進行驗證，需要進一步從細胞及分子層面闡明生姜揮發(fā)油中各個化合物，特別是在本研究中挖掘到的潛在化合物，在免疫和抗炎等信號通路的調(diào)控作用機制，明確化合物的有效性、安全性及起效用量，為臨床使用及產(chǎn)品開發(fā)提供科學依據(jù)。