沈 浮,沈樂樂,盧 敏,盧芳國,鄺高艷**,吳泳蓉,趙昱東,夏 雨,鄧 博**
(1. 湖南中醫藥大學研究生院 長沙 410208;2. 寧夏醫科大學研究生院 銀川 750004;3. 湖南中醫藥大學第一附屬醫院 長沙 410007;4. 湖南中醫藥大學醫學院 長沙 410208;5. 成都中醫藥大學研究生院 成都 610075)
骨關節炎(osteoarthritis,OA)和類風濕關節炎(rheumatoid arthritis,RA)分別為退行性關節炎及炎性關節炎的代表性疾病,是目前導致關節功能障礙和致殘的重要原因,其慢性持續的關節疼痛狀態嚴重影響人們生活質量[1,2]。OA 的基本病理特征主要包括軟骨細胞的退變和細胞外基質降解,常累及膝、髖等負重關節,RA 是一種以慢性滑膜炎及血管翳形成為主要病理改變的慢性自身免疫性疾病[3,4]。盡管兩者的發病機制有所不同,但均有 IL-1b、IL-6、TNF-α 等炎癥介質的參與[5]。目前尚無特效藥物,西醫治療以緩解疼痛癥狀、延緩病情進展為主。二者均屬于中醫“痹癥”的范疇,且病因病機也多有相似,素體虛弱,正氣不足,衛外不固,腠理不密,是誘發本病的重要因素[6,7],中醫辯證施治對于痹癥的治療具有獨特的優勢,且療效顯著、不良反應少,適宜長期服用,近年來愈發受到重視[8]。

表1 數據庫及軟件工具
追風透骨膠囊最早出處可追溯到明朝董宿《奇效良方》中的大追風散,“百節酸疼……項背拘急……面上游風,狀若蟲行,一切頭風”記載的大追風散初衷本用于祛風化痰而利頭,而痹癥同為風邪所害,并多兼有氣滯血瘀,后人配伍祛濕散寒、行氣活血之品,并經劑型優化改革,遂成追風透骨膠囊[6,9]。臨床上使用追風透骨膠囊治療OA、RA 療效顯著,在多項臨床研究中顯示,追風透骨膠囊可有效改善RA 患者的關節疼痛、腫脹及晨僵等臨床癥狀,并顯著降低RA 患者紅細胞沉降率、C反應蛋白以及類風濕因子等炎性指標[6,10],追風透骨膠囊亦可降低OA 患者的骨關節指數評分(VOMAC)、視覺模擬量表(VAS)評分而顯著改善患者關節活動度及疼痛反應,降低血清炎性介質水平,下調血液粘稠度,提高患者生活質量[6,11-13],對于痹癥的臨床治療有良好的應用價值。異病同治為中醫獨有的診療思維,最早出自《黃帝內經》,指在中醫理論指導下,不同疾病在其病情進展過程中出現相似的證候,從而采取同樣的的治療方法,可獲同樣理想的療效,反映了中醫辨證論治的特點[14,15]。追風透骨膠囊對于寒濕痹阻型的OA、RA 均具有良好的療效,體現了異病同治的理念[16]。
然而,由于中藥復方中成分多且復雜,疾病和中醫證候涉及到許多分子變化,因此大多數中藥方劑在復雜疾病中的組合規律和作用仍有待闡明[17]。隨著現代生物信息學的發展,海量信息獲取技術的成熟,尤其是計算能力及系統建模的發展,為中藥方劑在“生物網絡”的應用賦予了實質性的內涵[18],高通量信息技術為中藥復方干預疾病靶標的機制研究提供快捷的方法和新的機遇,并與中藥復方復雜成分及多途徑、靶點干預疾病的機理向貼合[19,20]。近年來,不少學者陸續開展了追風透骨膠囊治療OA、RA 的臨床對照試驗研究,但其內在機制尚缺乏完整的現代醫學證據。本研究利用高通量信息技術及分子對接法探討追風透骨膠囊對于OA、RA 異病同治的潛在分子機制進行預測,以期進一步拓展追風透骨膠囊異病同治的科學內涵,為該藥的深入研究和臨床應用提供更多可靠的理論依據。
如表1所示。
本研究依托TCMSP、BATMAN-TCM 數據庫檢索追風透骨膠囊的21 位味中藥(地龍、朱砂未收集在內):制 川 烏(chuanwu)、香 附(xiangfu)、川 芎(chuanxiong)、麻黃(mahuang)、制草烏(caowu)、當歸(danggui)、赤小豆(chixiaodou)、羌活(qianghuo)、赤芍(chishao)、細辛(xixin)、制天南星(tiannanxing)、白芷(baizhi)、甘草(gancao)、白術(baizhu)、沒藥(moyao)、乳香(ruxiang)、茯苓(fuling)、桂枝(guizhi)、天麻(tianma)、甘松(gansong)、防風(fangfeng)的所有化學成分。在TCMSP數據庫的篩選標準為,口服生物利用度(oral bioavailability,OB)≥30%并且類藥性指數(drug-likeness,DL)≥ 0.18[21],在 BATMAN-TCM 數據庫的篩選條件為:靶點蛋白為高可信度(prediction score cutoff > 20)[22],并將上述候選化合物的潛在靶點上傳至UniProt數據庫,進一步獲得對應的基因信息。

圖1 追風透骨膠囊藥物有效活性成分分布情況
OA、RA 相 關 疾 病 靶 點 來 源 于 GeneCards、DisGeNET 數據庫,將化合物作用靶點與疾病靶點取交集,獲取追風透骨膠囊干預兩者的的共有靶標,并導入Cytoscape 3.7.0 構建“藥物-化合物-共有靶標”網絡圖以備后續分析。
利用STRING 數據庫進行共有靶點的蛋白質間相互作用關系(protein-protein interaction,PPI)分析,設置combined score >0.9,并將結果利用Cytoscape 3.7.0繪制PPI網絡圖,并進行拓撲分析,以度值大于中位數的2倍篩選出關鍵共有靶點。
為深入探究對追風透骨膠囊異病同治相關生物途徑及機制,本研究運用R語言調用Bioconductor數據包進行追風透骨膠囊干預OA、RA 靶基因的基因本體功能富集(gene ontology,GO)與京都基因與基因組百科全書(Kyoto Encyclopedia of Genes and Genomics,KEGG)通路富集分析,設定P<0.05具有統計學意義,并對分析結果進行可視化。
為進一步驗證追風透骨膠囊中核心化合物與OA和RA 共同作用靶蛋白的結合情況,利用SYBYL-X 2.1.1 軟件中Surflex-dock 模塊進行分子對接,以追風透骨膠囊中核心化合物為配體,OA 和RA 共同作用靶蛋白為受體,該對接工具匹配方法主要基于結構及形狀的相似性,有高準確性、高陽性率及匹配速度迅速的特征[23]。一般認為,受體與配體的自由結合能越小,其穩定性越強,則Total score 值越大。Total score ≥5.0意味著配體與靶蛋白之結合活性較優,而Total score ≥7.0 則可預測該配體分子與靶蛋白結合活性較為強烈[24]。
通過檢索TCMSP、BATMAN-TCM 數據庫最終篩選收集到追風透骨膠囊中21 味藥物(地龍、朱砂未收集在內)所含的活性化學成分共276 種,其中川烏1個、香附 14 個、川芎 6 個、麻黃 21 個、草烏 4 個、當歸 2個、赤小豆 1 個、羌活 9 個、赤芍 11 個、細辛 6 個、天南星6個、白芷16個、甘草74個、白術4個、沒藥22個、乳香 2 個、茯苓 3 個、桂枝 6 個、天麻 5 個、甘松 4 個、防風17 個,活性化合物在各類中藥的分布比例如圖1A 所示,共獲取有效活性化合物靶標共475 個。由于各類中藥含有相同類型的活性成分,本研究組對各組分的中藥所含成分進歸納(圖1B),β-谷甾醇、谷甾醇、豆甾醇、槲皮素、茚酚、山柰酚在追風透骨膠囊中各組分中藥的分布范圍最廣。

圖2 追風透骨膠囊與疾病靶點韋恩圖
通過 GeneCards、DisGeNET 數據庫搜集 OA、RA 相關疾病靶點,去重后獲得RA 疾病靶點共4856 個,OA疾病靶點共3145 個,將以上2 組疾病靶點與追風透骨膠囊475 個作用靶點進行交集,共獲得交集靶標207個(圖2)。
通過Cytoscape 3.7.0 軟件構建“藥物-化合物-共有靶標”網絡圖(圖3),圖中共有節點410 個,邊1185條,該網絡的平均度值為5.78,大于其平均值2 倍的有效活性成分有12 個,分別是槲皮素(quercetin)、木犀草素(luteolin)、山奈酚(kaempferol)、間羥基苯甲酸(m-hydroxybenzoic acid)、黃芩苷(baicalein)、異鼠李素(isorhamnetin)、刺 槐 素(acacetin)、漢 黃 芩 素(wogonin)、β- 谷 甾 醇(beta-sitosterol)、楊 梅 酮(myricanone)、隱丹參酮(cryptotanshinone)、8-異戊烯-山奈酚(8-Isopentenyl-kaempferol),提示其為該網絡的關鍵性有效活性成分,在追風透骨膠囊干預OA 以及RA 的過程中其重要作用。從網絡圖可以看出,部分活性成分具有相同的作用靶點,說明各組分的藥效具有協同作用。HSP90AA1、ESR1、CALM1、NOS2、AR為追風透骨膠囊干預OA、RA的關鍵靶基因。
通過STRING 數據庫,進行PPI,并通過Cytoscape軟件進行可視化分析(圖4A),共有節點數173 個,邊數863 條,度值中位數為12,根據度值大于中位數2 倍篩選出核心網絡蛋白12 個(圖4B),分別為STAT3、AKT1、APP、JUN、MAPK1、TNF、HSP90AA1、RELA、IL6、MAPK14、VEGFA、MAPK8。這12 個靶蛋白為PPI網絡中的核心節點,很可能是追風透骨膠囊干預OA以及RA的關鍵靶蛋白。

圖3 追風透骨膠囊干預OA、RA的“藥物-化合物-共有靶標”網絡圖
根據P值的大小,本研究用條形圖列舉了各類排名靠前的GO 富集分析結果(圖5),其中生物過程(biological process,BP)2606個,主要涉及脂多糖反應、細菌來源分子反應、氧化應激反應、細胞氧化應激反應、金屬離子反應等(圖5A);細胞組成(cellular component,CC)85 個,主要涉及膜筏、膜微域、膜區、囊泡腔、細胞質囊泡腔等多種細胞結構(圖5B);分子功能(molecular function,MF)170 個,主要涉及細胞核受體活性、轉錄因子活性,直接配體調控特異性DNA 結合等(圖5C)。
在KEGG 富集結果篩選得到153 條主要信號通路,根據通路上富集的基因數目對排序前20的通路及對應的基因導入Cytoscape中繪制“靶基因-信號通路”網絡并進行可視化(圖6)。該網絡圖主要涉及代謝相關通路(AGE-RAGE signaling pathway in diabetic complications、Fluid shear stress and atherosclerosis)、病毒 感 染 相 關 通 路(Kaposi sarcoma-associated herpesvirus infection、Hepatitis B、Hepatitis C、Human cytomegalovirus infection)、炎性反應相關通路(IL-17 signaling pathway)、腫瘤相關通路(TNF signaling pathway、Prostate cancer)等多種信號途徑。

圖4 追風透骨膠囊作用OA、RA共同靶標的PPI網絡及核心靶蛋白
為進一步驗證追風透骨膠囊干預OA 以及RA 的分子機制,從PubChem 數據庫中提取追風透骨膠囊中的關鍵有效活性成分作為配體,從Protein Data Bank數據庫提取PPI中的核心靶蛋白作為受體進行分子對接,對接Total score 結果如圖7 所示,在144 個對接結果中,59個提示有較優的(Total score ≥5.0)對接活性,9 個有強烈的匹配活性(Total score ≥ 7.0),其中β-谷甾醇與核心靶蛋白結合情況最佳,HSP90AA1 是追風透骨膠囊中關鍵化學成分結合活性最優的靶蛋白。
Total score 值最高的4 組對接情況在如圖8 所示,8-異戊烯-山奈酚通過氨基酸殘基TYR139、LEU103、ASN51 形成3 條氫鍵(圖8A);β-谷甾醇通過氨基酸殘基 LYS58 形成 1 條氫鍵(圖 8C);兩者均能與HSP90AA1密切結合。β-谷甾醇與APP通過氨基酸殘基ASN84 形成1 條氫鍵(圖8B);異鼠李素通過氨基酸殘 基 ARG69、GLN37、LYS55、GLU109、LEU110、MET111、SER34 形成 7 條氫鍵與 MAPK-8 受體靶蛋白緊密結合(圖8D)。并且上述配體化合物均能良好的包埋在直徑為5A的受體靶蛋白的活性口袋中。

圖5 追風透骨膠囊治療OA、RA共同潛在靶點的GO富集分析

圖6 追風透骨膠囊干預OA、RA共同潛在靶點的“靶基因-信號通路”網絡

圖7 追風透骨膠囊中的關鍵活性化合物與PPI核心靶蛋白分子對接Total score值熱圖

圖8 追風透骨膠囊中關鍵活性化合物與PPI核心靶蛋白分子對接對接情況
OA 與RA 均屬于中醫“痹癥”范疇,其病機多為本虛標實,外感風、寒、濕之邪易誘發本病,追風透骨膠囊由制川烏、香附、川芎、麻黃等23 位藥物組成,其中細辛、川烏散寒溫經通絡,羌活、防風祛風除濕,并配伍赤芍、當歸等養血活血之品,共效活血通絡、散寒祛濕之功。現代研究證明,追風透骨膠囊具有顯鎮痛抗炎及抗凝的功效,能加速血液循環、改善血液粘稠度,調節機體免疫平衡[25],其改善關節炎癥反應可能與下調滑液及滑膜組織中miR-155、TGF-β1 的表達有關[26]。川烏在追風透骨膠囊為君藥,其不僅能顯著抑制前列腺素E、組胺等炎性介質的釋放,在體外還可抑制軟骨細胞基質降解,降低單碘乙酸鹽誘導的軟骨細胞損傷,對OA 模型具有一定的軟骨保護作用[27,28]。川芎與當歸的配伍增強了追風透骨膠囊的抗炎鎮痛效應,有研究顯示兩者配伍比例為1:1 時,其活血化瘀之功最優,這與追風透骨膠囊中二者的配比相吻合[29]。芍藥苷為追風透骨膠囊臣藥赤芍的主要成分,有報道顯示,芍藥苷能顯著改善RA 模型大鼠的足腫脹,降低關節內的 TNF-α、IL-1β、IL-6 等促炎因子的釋放,抑制Rho 激酶在滑膜組織內的激活,以及下調p-NF-κB、p65和p-MYPT1的表達[30]。
本研究利用網絡藥理學及分子對接的方法對追風透骨膠囊對于OA、RA 異病同治的機制進行探討,在“藥物-化合物-共有靶標”網絡圖中對追風透骨膠囊干預OA 及RA 共有靶標的關鍵有效活性成分進行了篩選,并通過分子對接技術對篩選得到的活性成分及進行驗證,其中槲皮素、木犀草素、山奈酚的度值排名靠前,而在分子對接結果中β-谷甾醇與核心靶蛋白結合情況最佳,提示它們在追風透骨膠囊干預OA、RA發揮的作用最大。追風透骨膠囊中關鍵活性成分多歸類于黃酮類單體化合物,現代研究表明槲皮素可減少 TNF-α、IL-1β、IL-6 等促炎性細胞因子生成,調控NF-κB、Nrf-2/HO-1 等多種炎癥相關信號通路,降解弱炎癥環境下軟骨細胞內MMP-13 產生、基質降解和細胞凋亡,從而保護關節軟骨,具有治療OA 以及RA的潛在藥物價值[31-34]。木犀草素是一種具有多種藥理特性的天然黃酮,由于其具有較強的抗炎作用,木犀草素已被用于治療多種炎性相關疾病,有報道顯示,木犀草素可通過下調的OA 軟骨細胞JNK、p38 的表達以及抑制NF-κB磷酸化,來抑制NO、TNF-α和IL-6的釋放,并逆轉II 型膠原蛋白的降解,從而減緩軟骨的退化[35,36]。山奈酚通過抑制bFGF-FGFR3–RSK2信號傳導抑制RA-FLSs 的增殖和遷移,從而減少IL-17、IL-21 和TNF-α等活化的T 細胞介導的炎性細胞因子的釋放,進而降低體外和體內的破骨細胞分化,對于RA 的誘導和進展發揮關鍵作用[37,38]。β-谷甾醇具有抗炎效應及免疫調節的功效,其可提高脂多糖(LPS)誘導的J774A1 巨噬細胞的IL-10 活性,減小趨化因子和促炎因子活性,亦可以增加Th-1 型淋巴細胞增殖,調節體內免疫從而抑制IL-6 和TNF 發揮抗炎作用[39]。本研究預測得出在追風透骨膠囊干預OA、RA 的過程主要涉及 HSP90AA1、ESR1、CALM1、NOS2、AR 等靶基因,其中Ⅱ型一氧化氮合酶(NOS2)相關基因與關節內滑膜炎癥密切,普遍認為高濃度的NO 抑制與線粒體代謝有關的多種酶,此外,在病理條件下,迅速分解為具有很強毒性的NO2 相關自由基從而對滑膜、軟骨細胞有直接殺傷作用[40]。CALM1 基因位于編碼Ca M染色體上,有研究顯示Ca M 與細胞內Ca2+結合的情況較為多見,并參與關節軟骨中軟骨細胞的分化及維持軟骨正常形態和功能[41]。本研究通過網絡藥理學對于追風透骨膠囊干預OA、RA 的關鍵活性化合物及潛在靶基因進行預測,部分結果干預OA、RA 的機制已在體內外實驗中得到體現,表明通過此法篩選的化學成分雖有一定的參考價值,另外,仍有部分預測結果在目前發表的文獻中尚未得到驗證,其可能是追風透骨膠囊對于OA、RA 異病同治的潛在活性成分,還需要結合進一步實驗研究。
根據PPI 核心蛋白的網絡分析中的度值排序,靠前的核心靶標主要與炎癥、免疫、增殖以及分化相關聯,分子對接結果提示追風透骨膠囊中關鍵活性成分與HSP90AA1 靶蛋白結合活性最優。HSP90AA1 誘導NF-κB 信號通路組分的抑制后降解,阻滯該信號通路的激活,并導致炎性細胞的細胞因子產生減少,并且越來越多的證據表明,HSP90AA1 相關的靶向藥物減少滑膜炎癥和防止軟骨破壞,有助于RA 的治療[42]。STAT3 是信號轉導與轉錄激活因子(STAT)家族成員之一,其激活可導致促炎細胞因子的產生和免疫反應增加,通過抑制STAT3 的表達可以限制免疫炎癥反應,改善自身免疫性相關的關節炎癥狀[43,44]。AKTI 通過Thr308 和Ser473 的磷酸化激活,活化的AKT1 通過多種信號途徑促進細胞生長,并使凋亡級聯的各種成分失活來參與OA 的進展[45]。IL-6 是重要的促炎因子,在RA、OA 患者的血清以及關節滑液中均明顯升高[46],其可增強 RA 發病過程中 IL-1 和 TNF-α的效應可加速B 淋巴細胞增殖,活化T淋巴細胞,破壞關節軟骨,并通過IL-6 通過JAK/STAT 信號通路的介導影響滑膜細胞纖維化以及骨的再吸收加速病情進展[47-48]。TNF-α亦扮演了重要的角色,可誘導效應T 細胞引起滑膜組織損傷[49,50]。MAPK8、MAPK14 均屬于絲裂原活化蛋白激酶家族,它們參與調節許多細胞功能,如分化、增殖、凋亡和炎癥反應,其參與RA 的疾病進程中得到廣泛驗證[51]。血管生成活動對于OA 與RA 的誘發及進展有著密切的聯系,血管內皮生長因子(VEGF)與內皮細胞(EC)上的同源受體結合,并激活這些細胞產生蛋白水解酶,其降解基底膜導致內皮細胞遷移并進一步增殖,這是RA 滑膜血管生成的重要促炎過程[52]。由此推測追風透骨膠囊對于OA、RA 的干預涉及多種靶蛋白,這些可能追風透骨膠囊對于兩者異病同治的潛在治療靶點,本研究預測的相關靶標基于高通量信息技術進行科學篩選,可為后續相關靶點的深入研究提供有意義的參考。
為了深入了解追風透骨膠囊對于OA 和RA 異病同治的作用機制,本課題組基于R 語言調用相關程序包對追風透骨膠囊干預OA、RA 的交集靶標進行GO及KEGG富集分析,GO富集分析提示追風透骨膠囊異病同治的生物學途徑主要涉及脂多糖生物反應、細菌分子應激反應以及氧化應激等生物學過程,以及核受體、細胞因子活性等多種分子功能。KEGG 分析結果發現炎性反應、代謝、病毒感染、腫瘤相關通路富集較多的靶向基因,其中AGE-RAGE 排名最靠前,對于晚期糖基化終產物(AGEs)信號傳導途徑的深入研究,使人們認識到AGEs 可以誘導人OA 滑膜細胞COX-2 的表達和前列腺素 E2(PGE2)、IL-6 和 MMP-13 的產生[53],另外AGEs 受體中和抗體(RAGE)能有效逆轉人滑膜細胞的炎癥反應和VEGF 生成,提示RAGE 在滑膜細胞活化中起重要作用,從而促進OA 的進展[54]。IL-17通路在RA進程中被MMP13被激活,miR-19a可以通過介導IL-17 信號通路介導HFLS-RA 的增殖和侵襲,從而調節RA 人成纖維樣滑膜細胞中MMP13 的表達,從而參與軟骨細胞的降解[55]。可以發現,追風透骨膠囊作用對于OA、RA 異病同治作用途徑的范圍十分寬泛,并可以通過多種途徑對不同的信號通路及生物學過程進行干預,由于通路的篩選基于相關靶基因在備選通路的富集數目進行篩選,不排除部分通路的預測出現假陽性的結果,因此,后續科研驗證工作仍需結合OA、RA 相關前期文獻研究進行綜合分析,本課題組認為追風透骨膠囊發揮異病同治的機制非常復雜,本次預測篩選的炎癥、代謝及腫瘤相關的通路具有后續研究的參考意義。
本研究仍存在一些不足之處,首先,由于數據庫的局限性,追風透骨膠囊中地龍、朱砂尚未納入本次研究范圍內,且本研究依賴于TCMSP、BATMAN-TCM等數據庫對于追風透骨膠囊各活性成分的收集,由于數據庫本身可能存在對于藥物活性成分及對應靶標的欠完整性及準確性不高等局限,對于各組分中藥在體內的相互作用及代謝情況尚不明確,這與追風透骨膠囊實際干預OA、RA 的機制可能存在偏倚,后期研究可以通過液相色譜等方法來獲取活性成分來提高準確性。其次,根據中醫辯證論治的基本原則,追風透骨膠囊干預寒濕痹阻證的痹癥更符合中醫的病理病機,由于目前暫缺乏中醫證候相關的基因數據庫,因此本研究尚無法做到方證對應,后續研究可以通過基因芯片方法得到該證候相關的差異基因并加以實驗驗證。綜上所述,通過網絡藥理學方法及分子對接技術對追風透骨膠囊干預OA、RA 異病同治作用的機制進行預測,篩選得到的關鍵活性化合物、核心靶標以及相關通路與目前研究OA、RA 的機制具有一致性,并且在分子對接驗證結果中多數活性成分與靶蛋白結合匹配度較高,從一定程度上表明預測結果的可靠性與準確性,體現了追風透骨膠囊具有多成分、多靶點、多途徑的特點,并豐富了追風透骨膠囊異病同治的科學內涵,可為下階段的研究提供參考。