基于不同病理階段生物信息學(xué)分析篩選強直性脊柱炎治療中藥有效成分

王志強，宮彩霞，李振彬

強直性脊柱炎(AS)是一種以骶髂關(guān)節(jié)和脊柱炎癥和結(jié)構(gòu)性損傷為特點的慢性系統(tǒng)性疾病[1]。疾病早期以骶髂關(guān)節(jié)炎、附著點炎和脊柱炎為主要表現(xiàn)，疾病后期因進行性脊柱韌帶和小關(guān)節(jié)骨化，導(dǎo)致固定性胸腰椎后凸畸形[2]。目前，AS的主要治療目標(biāo)是通過控制炎癥和癥狀來維持患者正常軀體功能，防止進行性結(jié)構(gòu)性損傷，改善遠期生活質(zhì)量[3]。目前，臨床指南推薦AS一線治療藥物為非甾體類抗炎藥(NSAIDs)[4-5]。此外，糖皮質(zhì)激素、抗風(fēng)濕藥也是臨床常用的治療藥物[6]。雖然這些藥物有一定療效，但是長期治療可能出現(xiàn)多種不良反應(yīng)[6]。臨床研究證實，中藥具有良好的抗炎、鎮(zhèn)痛效果，能有效緩解AS炎性癥狀、提高生活質(zhì)量[7-8]。本研究基于對AS不同病理階段的生物信息學(xué)分析，探討AS的發(fā)病機制，并篩選出可能有效的中藥成分，為進一步研究提供理論依據(jù)。

1 研究方法

1.1AS不同病理階段靶點的收集和篩選 分別以“ankylosing spondylitis”、“sacroiliac arthritis(骶髂關(guān)節(jié)炎)”、“enthesitis(附著點炎)”、“spinal ligament ossification(脊柱韌帶骨化)”和“facet joint ossification(小關(guān)節(jié)骨化)”為關(guān)鍵詞，在GeneCards數(shù)據(jù)庫檢索相應(yīng)的疾病靶點。“ankylosing spondylitis”、“sacroiliac arthritis”和“enthesitis”靶點取交集，獲得炎癥損傷階段靶點；“ankylosing spondylitis”、“spinal ligament ossification”和“facet joint ossification”靶點取交集，獲得結(jié)構(gòu)性損傷階段靶點。

1.2蛋白質(zhì)相互作用(PPI)網(wǎng)絡(luò)的構(gòu)建及分析 分別把炎癥損傷和結(jié)構(gòu)性損傷階段靶點導(dǎo)入STRING平臺(https://string-db.org/，Version 11.0)[9]，蛋白種屬選擇“人類”，構(gòu)建PPI網(wǎng)絡(luò)，各參數(shù)保持默認(rèn)設(shè)置，用Cytoscape軟件(http://www.cytoscape.org/，Version 3.7.2)[10]實現(xiàn)網(wǎng)絡(luò)可視化。利用Cytoscape軟件中的MCODE插件對網(wǎng)絡(luò)進行拓?fù)浞治?，利用cytoHubba插件獲得PPI網(wǎng)絡(luò)排名前10的核心靶點，各參數(shù)保持默認(rèn)。

1.3功能注釋和通路富集分析 使用R軟件(Version 3.6.1)通過ClusterProfiler包(Version 3.18.0)[11]對PPI網(wǎng)絡(luò)中經(jīng)MCODE分析獲得的評分最高的Cluster靶點進行GO和KEGG通路富集分析，經(jīng)P<0.05、Q<0.05篩選后，取排名前20的生物學(xué)功能和通路進行可視化處理。

1.4AS不同病理階段的中藥有效成分篩選 在中藥系統(tǒng)藥理數(shù)據(jù)庫與分析平臺(TCMSP，http://tcmspw.com/tcmsp.php，Version 2.3)[12]檢索可作用于cytoHubba插件獲得的核心靶點的所有中藥有效成分，并以口服生物利用度(OB)≥30%和類藥性(DL)≥0.18作為進一步篩選條件。

1.5分子對接 將篩選得到的中藥有效成分作為配體，目前指南推薦的NSAIDs、腫瘤壞死因子抑制劑和IL17A單抗的作用靶點PTGS2、TNF和IL17A作為受體，用Discovery Studio軟件(Version 16.1.0)的CDOCKER模式進行分子對接驗證。中藥有效成分和靶點的3D結(jié)構(gòu)分別從PubChem數(shù)據(jù)庫(https://pubchem.ncbi.nlm.nih.gov/)[13]和RCSB數(shù)據(jù)庫(http://www.rcsb.org/)[14]下載。

2 結(jié)果

2.1AS不同病理階段靶點的篩選 在GeneCards數(shù)據(jù)庫獲得相應(yīng)靶點2123、251、241、1174和458個。取交集后，獲得炎癥損傷階段靶點103個，結(jié)構(gòu)性損傷階段靶點109個，見圖1。

圖1 AS不同病理階段靶點的篩選

2.2PPI網(wǎng)絡(luò)的構(gòu)建及分析 炎癥損傷階段的PPI網(wǎng)絡(luò)包括99個節(jié)點和1267條連線，結(jié)構(gòu)性損傷階段的PPI網(wǎng)絡(luò)包括107個節(jié)點和1916條連線，見圖2。經(jīng)MCODE插件進行拓?fù)浞治龊?，炎癥損傷階段的PPI網(wǎng)絡(luò)共獲得7個Cluster，結(jié)構(gòu)性損傷階段的PPI網(wǎng)絡(luò)共獲得3個Cluster，見表1。利用cytoHubba插件獲得2個PPI網(wǎng)絡(luò)排名前10的核心靶點，其中炎癥損傷階段的核心靶點為IL6、TNF、IL10、CSF2、IL17A、IL1B、CXCL8、IFNG、TLR4、IL18，結(jié)構(gòu)性損傷階段的核心靶點為IL6、GAPDH、VEGFA、INS、CXCL8、CCL2、STAT3、TNF、IL1B、PTGS2。

表1 AS 2個病理階段PPI網(wǎng)絡(luò)的MCODE拓?fù)浞治鼋Y(jié)果

圖2 AS 2個病理階段的PPI網(wǎng)絡(luò)和核心靶點

2.3功能注釋和通路富集分析 對AS 2個病理階段PPI網(wǎng)絡(luò)經(jīng)MCODE拓?fù)浞治龊?，獲得分值最高的43、42個靶點進行GO和KEGG通路富集分析。炎癥損傷階段靶點分別富集到生物過程(BP)1020個、細(xì)胞組分(CC)13個、分子功能(MF)30個、通路66條；結(jié)構(gòu)性損傷階段分別富集到BP 1866個、CC 16個、MF 47個、通路139條。2個病理階段排名前20的生物學(xué)功能和通路見圖3。在富集分析得到的信號通路中，AS 2個病理階段有高達65條相同的信號通路，其中排名前20的相同信號通路包括細(xì)胞因子-細(xì)胞因子受體相互作用、類風(fēng)濕關(guān)節(jié)炎、炎性腸病、瘧疾、結(jié)核、甲型流感、百日咳、Chagas病、IL17、Th17細(xì)胞分化等。

圖3 AS 2個病理階段靶點的GO和KEGG通路富集分析

2.4AS不同病理階段核心靶點的中藥有效成分篩選 在TCMSP數(shù)據(jù)庫中檢索，可作用于“2.2”中獲得的AS不同病理階段核心靶點的所有中藥有效成分?？赏瑫r作用于炎癥損傷階段5個以上核心靶點的中藥有效成分為槲皮素(6個)和白藜蘆醇(5個)，可同時作用于結(jié)構(gòu)性損傷階段5個以上核心靶點的中藥有效成分包括白藜蘆醇(8個)、金雀異黃素(8個)、槲皮素(7個)、尼古丁(7個)、熊果酸(6個)、黃體酮(6個)、表沒食子兒茶素沒食子酸酯(5個)、蘆丁(5個)、漢黃芩素(5個)、漆黃素(5個)、楊梅黃酮(5個)、雷公藤甲素(5個)。上述中藥有效成分經(jīng)OB≥30%和DL≥0.18篩選后，槲皮素、表沒食子兒茶素沒食子酸酯、漢黃芩素、漆黃素和雷公藤甲素符合要求。

2.5中藥有效成分與3種靶點蛋白的分子對接 將“2.4”中篩選得到的中藥有效成分作為配體，PTGS2、TNF和IL17A作為受體，Discovery Studio軟件分析結(jié)果顯示，所有配體均可穩(wěn)定地對接到受體的對接口袋中，見圖4。

圖4 篩選到的治療強直性脊柱炎有效中藥成分與關(guān)鍵靶點的對接模式圖

3 討論

在疾病后期，AS患者通常會出現(xiàn)脊柱固定性畸形，隨著進行性胸椎后凸增大、腰椎前凸減小，最終導(dǎo)致全脊柱矢狀面失衡。對于AS患者來說，嚴(yán)重的脊柱畸形不僅是一種外表和心理困擾，而且嚴(yán)重影響患者日常生活和工作質(zhì)量。近年研究發(fā)現(xiàn)，炎癥和結(jié)構(gòu)性損傷是AS脊柱活動障礙的2個重要的生物學(xué)機制[1]。OASIS研究顯示，AS功能損傷受到患者疾病活動性和腰頸椎結(jié)構(gòu)性損傷程度的影響[1]。脊柱的可逆性炎癥和不可逆性結(jié)構(gòu)性損傷導(dǎo)致AS患者脊柱活動障礙。AS早期，脊柱活動障礙更多歸因于炎癥，而結(jié)構(gòu)性損傷是后期病變的主要原因，脊柱韌帶和小關(guān)節(jié)進行性骨化最終導(dǎo)致嚴(yán)重的脊柱后凸畸形[1]。然而，脊柱炎癥和結(jié)構(gòu)性損傷對AS最終脊柱畸形的“貢獻度”仍不明確[1]。

最近研究表明，IL23/IL17信號通路激活與附著點炎和附著點新骨形成有關(guān)[15-16]，IL17可以增強T細(xì)胞活化并刺激多種細(xì)胞產(chǎn)生促炎介質(zhì)[17]。NOD樣受體家族和NLRP3炎癥小體活性在多種炎癥性疾病的發(fā)生發(fā)展中起重要作用，且NLRP3炎癥小體與AS易感性相關(guān)[18]。干擾素-γ(IFN-γ)是已知的Th1型細(xì)胞因子，是唯一的Ⅱ型干擾素。IFN-γ可拮抗Th2型細(xì)胞因子IL4和IL10的分泌，從而抑制Th2細(xì)胞反應(yīng)。有研究表明，IFN-γ表達與rs2430561基因多態(tài)性明顯相關(guān)，提示rs2430561等位基因可能調(diào)節(jié)IFN-γ表達，從而導(dǎo)致AS發(fā)生[19]。IL4通過調(diào)節(jié)巨噬細(xì)胞極化抑制AS發(fā)病[20]。因此，Th1/Th2細(xì)胞失衡也參與了AS的發(fā)病。骨代謝高度依賴于JAK2/STAT3介導(dǎo)的成骨細(xì)胞調(diào)控[21]。核因子-κB(NF-κB)參與了機體炎癥、免疫、細(xì)胞增殖和存活等多種重要過程[22]。NF-κB信號通路高活性會增加AS的患病風(fēng)險[23]。TNF在AS發(fā)病中的作用尚未完全闡明，這是因其在AS炎癥損傷和結(jié)構(gòu)性損傷中的作用并不完全一致[24]。在炎癥損傷階段，TNF與NF-κB受體活化因子配體(RANKL)協(xié)同刺激破骨細(xì)胞形成，增強調(diào)節(jié)成骨細(xì)胞分化和活性的蛋白酶降解從而抑制骨形成，同時誘導(dǎo)細(xì)胞外消化酶表達。MMPs和蛋白聚糖酶可降解膠原蛋白、聚集蛋白聚糖及其他軟骨和椎間盤成分，促進局部骨質(zhì)侵蝕。在結(jié)構(gòu)性損傷階段，TNF可同時激活NF-κB及其下游的BMPs、osterix、矮小相關(guān)基因2、骨鈣素和堿性磷酸酶促進成骨分化。TNF對成骨細(xì)胞形成的作用取決于反應(yīng)細(xì)胞的濃度、暴露時間和分化狀態(tài)。此外，TNF還能通過炎癥級聯(lián)反應(yīng)直接或間接激活Wnt/β-catenin信號通路。Wnt/β-catenin信號通路是一條調(diào)控成骨細(xì)胞分化的骨穩(wěn)態(tài)調(diào)控通路[24]。本研究生物信息學(xué)分析顯示，炎癥損傷階段和結(jié)構(gòu)性損傷階段均涉及上述重要信號通路，且2個階段富集通路有高達65條相同的通路，提示炎癥機制貫穿著AS的始終，炎癥長期控制似乎是避免因結(jié)構(gòu)性損傷而導(dǎo)致殘疾的關(guān)鍵[25]。

根據(jù)AS 2個病理階段PPI網(wǎng)絡(luò)的核心靶點，本研究初步篩選出可用于炎癥損傷階段的中藥有效成分是槲皮素，可用于結(jié)構(gòu)性損傷階段的中藥有效成分包括槲皮素、表沒食子兒茶素沒食子酸酯、漢黃芩素、漆黃素和雷公藤甲素。槲皮素能減輕關(guān)節(jié)炎小鼠的關(guān)節(jié)癥狀并降低TNF-α、IL1β、IL17和單核細(xì)胞趨化蛋白-1水平，還可通過調(diào)控NF-κB信號通路減輕炎癥反應(yīng)[26]。此外，槲皮素通過對RANKL/NF-κB受體活化因子/骨保護素信號通路、絲裂原活化蛋白激酶(MAPK)信號通路、Wnt/β-Catenin信號通路、BMPs和轉(zhuǎn)化生長因子-β(TGF-β)信號通路的調(diào)控來影響成骨細(xì)胞和破骨細(xì)胞的發(fā)育和功能[27]。表沒食子兒茶素沒食子酸酯可通過TLR4/NF-κB信號通路抑制脂多糖誘導(dǎo)的RAW264.7巨噬細(xì)胞炎癥反應(yīng)[28]，并能糾正AS模型大鼠Th1/Th2免疫失衡和促進軟骨細(xì)胞凋亡[29]；其還可能通過與雌激素受體相互作用進而激活Wnt/β-catenin信號通路，促進成骨細(xì)胞的增殖和骨向分化[30]。漢黃芩素可通過調(diào)控NF-κB、MAPK、TGF-β1/Smad3信號通路抑制炎癥反應(yīng)[31-32]，同時對破骨細(xì)胞分化有一定影響[33]。漆黃素通過上調(diào)RUNX2轉(zhuǎn)錄活性促進成骨細(xì)胞分化[34]，通過阻斷RANKL誘導(dǎo)活性氧簇產(chǎn)生抑制破骨細(xì)胞分化[35]。雷公藤甲素是中藥雷公藤的主要活性成分，雷公藤甲素的抗炎和免疫抑制特性使其成為治療炎癥性疾病非常有前途的藥物。雷公藤甲素可減少炎性介質(zhì)的產(chǎn)生，減輕動物模型及相關(guān)細(xì)胞系受損器官的氧化應(yīng)激反應(yīng)，抑制T細(xì)胞、巨噬細(xì)胞等炎性細(xì)胞活性，還可通過對RANKL/RANK/OPG和NF-κB信號通路的調(diào)控抑制破骨細(xì)胞炎癥活性[36]。因此，篩選出這些中藥有效成分對炎癥和骨形成過程有干預(yù)作用。經(jīng)初步的分子對接驗證，上述中藥有效成分均可與目前指南推薦藥物的靶點穩(wěn)定對接。

綜上，根據(jù)AS不同病理階段生物信息學(xué)分析，初步篩選出一些潛在的中藥有效成分，為進一步的藥物研發(fā)和疾病治療提供有意義的參考。