李瑞平 雷雨歌 李 健 張 越

沈陽醫(yī)學院公共衛(wèi)生學院，遼寧沈陽 110034

砷作為已知的人類致癌物，廣泛存在于自然界 中[1]。通過環(huán)境污染和職業(yè)接觸進入體內(nèi)，對人體皮膚、神經(jīng)系統(tǒng)、心血管、肝臟及腎臟等器官有明顯的損傷[2]。砷在體內(nèi)代謝過程中產(chǎn)生多種自由基和非自由基產(chǎn)物，引起細胞紊亂，直接攻擊細胞或誘發(fā)脂質(zhì)過氧化物引起機體氧化與抗氧化代謝失調(diào)，造成氧化應激損傷[3]。氧化應激的產(chǎn)生使得大量的ROS 產(chǎn)生，造成細胞的損傷和凋亡[4]。姜黃素（Cur）是姜黃（姜黃科）根莖的主要Cur 類化合物，傳統(tǒng)上被用來治療咳嗽、黃疸、感冒、肝臟疾病和炎癥性疾病，被認為是一種安全性能高、療效好、副作用小的藥物，可作為糖尿病及心血管并發(fā)癥的治療藥物[5]。但Cur 對砷誘導心臟損傷的具體機制尚不明確，值得進一步研究。本研究旨在通過網(wǎng)絡(luò)藥理學方法探索Cur 抗砷誘導心臟損傷的潛在靶點和途徑。以砷處理的H9C2 細胞作為體外心肌細胞模型，驗證Cur 對細胞凋亡的調(diào)控作用及其機制。結(jié)果表明，Cur 至少部分地通過激活PI3K/AKT 通路來減輕亞砷酸鈉（NaAsO2）介導的細胞凋亡。

1 材料與方法

1.1 Cur靶點的預測

從本草組鑒HERB 數(shù)據(jù)庫（http://herb.ac.cn/）中篩選Cur 的潛在靶點，該數(shù)據(jù)庫是中藥高通量實驗和參考數(shù)據(jù)庫。

1.2 砷致心臟損傷疾病靶點預測

從GeneCard、在線人類孟德爾遺傳（OMIM）數(shù)據(jù)庫獲得疾病相關(guān)靶點，關(guān)鍵詞為砷誘導心臟損傷，并去除兩個數(shù)據(jù)庫的重復靶點。

1.3 疾病和藥物共同靶點的篩選

將藥物和疾病的靶點輸入Venny（https://bioinfogp.cnb.csic.es/tools/venny/）平臺，得到藥物和疾病的交集靶點及韋恩圖。

1.4 交集靶點蛋白質(zhì)相互作用（PPI）網(wǎng)絡(luò)構(gòu)建

將藥物與疾病的交集靶點輸入String 數(shù)據(jù)庫，將物種設(shè)定為人“Homo Sapiens”，最低相互作用閾值設(shè)為“medium confidence（＞0.700）”，將其導入Cytoscape 3.7.2 軟件，繪制PPI 網(wǎng)絡(luò)。Cytohubba 是Cytoscape 軟件中尋找核心靶點的一個插件，用于發(fā)現(xiàn)復雜網(wǎng)絡(luò)中的關(guān)鍵目標[6]。利用Cytohubba 中的五種算法，節(jié)點連接度（Deg）、最大集團中心性（MCC）、最大領(lǐng)域分量（MNC）、邊緣滲透分量（EPC）與Closeness（Clo）選取前5 個基因，根據(jù)5 種算法篩選后的基因取交集，如果同一基因同時能在五種算法之中，則認為是核心基因。

1.5 基因功能和信號通路的富集分析

GO 是注釋基因及其表達產(chǎn)物的常用方法，由生物過程（biological process）、細胞組成（cellular component）和分子功能（molecular function）3 部分組成，運用DAVID 數(shù)據(jù)庫對疾病和藥物的交集靶點進行GO 功能分析，說明其在基因功能中的作用。對Cur 和砷導致的心臟損傷疾病的交集靶點進行KEGG 信號通路富集分析，從而了解Cur 治療該疾病發(fā)生改變的關(guān)鍵信號通路。

1.6 關(guān)鍵靶點與Cur的分子對接

將排名前三的靶點與Cur 進行分子對接，通過PubChem 數(shù)據(jù)庫下載成分的結(jié)構(gòu)，通過PDB 數(shù)據(jù)庫下載排名靠前的核心靶點結(jié)構(gòu)，將成分及靶點結(jié)構(gòu)通過AutoDock 軟件進行加氫、計算電荷處理后進行分子對接，結(jié)合能最低的結(jié)構(gòu)為最優(yōu)結(jié)構(gòu)，通過Open Babel GUI 軟件將對接結(jié)果的.pdbqt 格式轉(zhuǎn)化為.pdb 格式，最后使用PyMOL 軟件對圖像進行處理分析。

1.7 細胞培養(yǎng)

H9C2 細胞購自中國科學院，NaAsO2來自美國Sigma 公司。

1.8 細胞活力

將H9C2 心肌細胞接種與96 孔板中，接種密度為5×104個/ml，鋪板24 h 后對細胞進行NaAsO2染毒實驗。實驗的分組設(shè)置為：對照（Con）組、NaAsO2組（6、8、10 μmol/L），NaAsO2（6 μmol/L）與Cur組（2.5、5、7.5、10、15 μmol/L）。

1.9 熒光定量聚合酶鏈式反應

取大鼠心肌細胞，提取RNA 并按照反轉(zhuǎn)錄試劑盒（Takara）將RNA 逆轉(zhuǎn)錄成cDNA，以GAPDH為內(nèi)參，RT-PCR 檢測各因子的表達，引物序列由上海生工公司合成，以GAPDH 為參照基因。

1.10 統(tǒng)計學分析

采用SPSS 22.0 統(tǒng)計學軟件進行數(shù)據(jù)分析，計量資料用均數(shù)±標準差（± s）表示，采用t檢驗；P＜ 0.05 為差異有統(tǒng)計學意義。

2 結(jié)果

2.1 靶點的預測

從HERB 數(shù)據(jù)庫中篩選出Cur 靶點152 個。以“Arsenic induces heart damage”為關(guān)鍵詞，GeneCards數(shù)據(jù)庫以“Relevance score”大于1 進行篩選，得出疾病靶點2207 個。OMIN 數(shù)據(jù)庫中篩選出197 個疾病靶點。

2.2 PPI網(wǎng)絡(luò)的構(gòu)建

將篩選出來的Cur 靶點與疾病靶點，運用VENNY 2.1 繪制韋恩圖，見圖1A。將兩者的共同靶點導入String 數(shù)據(jù)庫，根據(jù)設(shè)定條件得到PPI 網(wǎng)絡(luò)。運用Cytoscape 3.7.2 軟件導入節(jié)點信息，見圖1B。其中，根據(jù)5 種算法獲得5 個核心基因分別是AKT1、AMPK3、JUN、TP53、SRC。

2.3 GO和KEGG通路分析

DAVID 數(shù)據(jù)庫以FDR＜0.05 對結(jié)果進行篩選，對119 個交集靶點進行GO 功能富集分析和KEGG通路富集分析。共有973 個GO 術(shù)語條目，主要的GO 生物過程包括RNA 聚合酶Ⅱ啟動子轉(zhuǎn)錄的正調(diào)控、凋亡過程的負調(diào)控等，表明Cur 對砷誘導的心臟損傷有抗凋亡作用。細胞組分主要與細胞質(zhì)、細胞核有關(guān)。分子功能與蛋白結(jié)合與酶結(jié)合有關(guān)，揭示了蛋白質(zhì)在疾病的重要性，見圖2。主要KEGG 信號通路有癌癥通路、PI3K/AKT 信號通路，PI3K/AKT 通路富集程度較高，有34 個靶基因，見圖3。

圖2 Cur 治療砷誘導心臟損傷GO 富集分析

圖3 Cur 治療砷誘導心臟損傷KEGG 富集分析

2.4 分子對接

將Cur 與排名前3 位靶點進行分子對接，值越小表示結(jié)合活性越高，成分越容易與受體結(jié)合，所示對接結(jié)果見表1。分子與標靶蛋白的結(jié)合能均小于0，說明配體與受體均可以自發(fā)結(jié)合。分子對接圖見圖4。

表1 分子對接結(jié)果圖

圖4 Cur 治療砷誘導心臟損傷的分子對接示意圖

2.5 細胞活力

CCK8 結(jié) 果 表 明 從6 μmol/L 時 開 始，隨 著NaAsO2染毒劑量增加，細胞活性下降（P＜ 0.05）。在Cur 與NaAsO2聯(lián)合處理實驗中Cur 濃度為2.5 μmol/L 時，細胞活性上升最顯著，因此，在后續(xù)的實驗選用2.5 μmol/L Cur 和6 μmol/L NaAsO2進行研究。見圖5。

圖5 Cur 和NaAsO2 對H9C2 細胞存活率的影響（n=6）

2.6 Cur通過PI3K/AKT信號通路在砷誘導心肌細胞中發(fā)揮抗凋亡作用

與Con 組相比，NaAsO2組Caspase-9 的mRNA水 平 升 高，AKT、BCL-2 的mRNA 水 平 降 低; 與NaAsO2組相比，Cur+NaAsO2聯(lián)合組中PI3K、AKT mRNA 水平升高，Caspase-9、BAX 的mRNA 水平降低，見圖6。

圖6 H9C2 細胞中PI3K（A）、AKT（B）、Caspase-3（C）、Caspase-9（D）、Bax（E）、BCL-2（F）的mRNA 相對表達水平（n=3）

3 討論

全世界有超過3 億人口面臨高砷暴露的威脅，我國因砷暴露受累人群近2000 萬[7]。本研究采用網(wǎng)絡(luò)藥理學、分子對接和體外實驗驗證的方法，揭示Cur 對砷誘導心臟的作用機制。從數(shù)據(jù)庫中選取Cur 152 個靶點，PPI 網(wǎng)絡(luò)構(gòu)建中根據(jù)5 種算法得出AKT1、MAPK3、TP53、JUN、SRC5 個核心基因發(fā)揮著重要作用。分子對接結(jié)果顯示AKT1 與Cur結(jié)合活性最好。AKT1 是AKT 重要的亞型，AKT 一般指蛋白激酶B（protein kinase B，PKB），是一種絲氨酸/蘇氨酸蛋白激酶，也是PI3K 下游的主要靶點之一。PI3K/AKT 信號通路具有廣泛的生物學效應，參與調(diào)節(jié)細胞分裂、增殖、凋亡、代謝等多個環(huán)節(jié)，與心血管疾病密切相關(guān)[8]。MAPK3、JUN 與Cur 結(jié)合活性較好。在人類中，MAPK3 在細胞周期進展和凋亡中起關(guān)鍵作用[9]。c-jun 通過促進細胞存活和抑制細胞凋亡，在低氧心肌反應中起到保護作用。GO 富集分析的生物學過程涉及細胞凋亡過程。KEGG 信號通路分析表明PI3K/AKT 信號通路富集度較高，體外實驗也證實了Cur 通過激活PI3K/AKT 信號通路抑制砷處理的心肌細胞的凋亡。

PI3K/AKT 是調(diào)節(jié)心臟功能的一條重要的信號轉(zhuǎn)導通路，主要通過血管內(nèi)皮生長因子介導血管生成、抑制心肌細胞凋亡、重構(gòu)心室、促進細胞能量代謝等方面調(diào)節(jié)心臟功能[10]，該信號在心肌細胞存活和程序性死亡的過程中發(fā)揮關(guān)鍵的調(diào)控作用[11]。此通路與線粒體途徑細胞凋亡過程中Bcl-2 家族蛋白的活性密切相關(guān)，Bcl-2 蛋白家族是線包括抗凋亡蛋白Bcl-2 和促凋亡蛋白Bax 等[12]。Caspase-9作為線粒體途徑細胞凋亡的起始者，活化的Caspase-9 進一步裂解Caspase-3 和Caspase-7，凋亡過程開始[13]。研究表明，Cur 可以通過抑制心肌細胞凋亡，進而降低冠心病的發(fā)病風險[14]。Cur 可以改善糖尿病心肌病大鼠心室重構(gòu)[15]。

綜上所述，本研究在網(wǎng)絡(luò)藥理學的基礎(chǔ)上，將PI3K/AKT 通路與Cur 對砷誘導心臟的保護作用聯(lián)系起來，并利用體外實驗驗證了網(wǎng)絡(luò)藥理學的分析。通絡(luò)網(wǎng)絡(luò)藥理學挖掘了潛在通路，這些通路也值得研究，但本研究未對其他通路進行驗證。本研究只利用體外實驗進行驗證，具有一定的局限性。