基于真實世界及靶點網(wǎng)絡(luò)的鹿瓜多肽注射液治療骨關(guān)節(jié)炎的人群特征、療效、聯(lián)合用藥及機制預(yù)測分析

劉 毅，謝雁鳴*，黎元元*，高明慧

·數(shù)據(jù)挖掘與循證醫(yī)學(xué)·

基于真實世界及靶點網(wǎng)絡(luò)的鹿瓜多肽注射液治療骨關(guān)節(jié)炎的人群特征、療效、聯(lián)合用藥及機制預(yù)測分析

劉 毅1，謝雁鳴1*，黎元元1*，高明慧2

1. 中國中醫(yī)科學(xué)院中醫(yī)臨床基礎(chǔ)醫(yī)學(xué)研究所，北京 100700 2. 中國人民大學(xué)統(tǒng)計學(xué)院，北京 100872

探討在真實世界中鹿瓜多肽注射液治療骨關(guān)節(jié)炎的人群特征及聯(lián)合用藥情況，并通過構(gòu)建靶點網(wǎng)絡(luò)預(yù)測其作用機制。通過調(diào)取醫(yī)院信息系統(tǒng)（hospital information systeam，HIS）數(shù)據(jù)庫中使用過鹿瓜多肽注射液治療骨關(guān)節(jié)炎患者信息，使用SAS9.2及R軟件對其使用人群特征及用藥情況進(jìn)行統(tǒng)計；通過分析應(yīng)用鹿瓜多肽注射液的骨關(guān)節(jié)炎患者的治療結(jié)局以判斷其有效性，并與未使用鹿瓜多肽的骨關(guān)節(jié)炎患者的治療結(jié)局進(jìn)行對比分析；Clementine 12.0對聯(lián)合用藥情況進(jìn)行關(guān)聯(lián)規(guī)則分析，并按照置信度排序；利用CNKI、TCMSP數(shù)據(jù)庫及Cytoscape軟件構(gòu)建“疾病-藥物-成分-靶點”網(wǎng)絡(luò)及蛋白相互作用網(wǎng)絡(luò)，并進(jìn)行拓?fù)鋵W(xué)分析篩選核心成分與核心靶點；將藥物靶點與疾病靶點的交集靶點進(jìn)行基因本體（gene ontology，GO）及京都基因與基因組百科全書（Kyoto encyclopedia of genes and genomes，KEGG）富集分析；最后將核心成分與核心靶點進(jìn)行分子對接。4050例骨關(guān)節(jié)炎患者中，46～65歲年齡段最多（2248例，55.51%），其次為66～75歲年齡段（850例，20.99%），其中女性患者占多數(shù)（2657例，70.35%）。4050例應(yīng)用鹿瓜多肽注射液治療的骨關(guān)節(jié)炎患者的總有效率為97.88%，對比未應(yīng)用鹿瓜多肽注射液的骨關(guān)節(jié)炎患者的治療效果，顯示鹿瓜多肽注射液療效更好。關(guān)聯(lián)規(guī)則顯示鹿瓜多肽注射液治療骨關(guān)節(jié)炎常與激素（倍他米松）、抗生素（頭孢唑啉）、抗炎鎮(zhèn)痛藥（骨肽注射劑）、補腎壯骨劑（金天格膠囊）、活血化瘀劑（脈血康膠囊）合并用藥。靶點網(wǎng)絡(luò)分析得到核心成分棕櫚酸、肉豆蔻酸、色氨酸、維生素A等10種；核心靶點過氧化物酶體增殖物激活受體（peroxisome proliferator-activated receptor gamma，PPARG）、髓過氧物酶（myeloperoxidase，MPO）、受體型酪氨酸蛋白磷酸酶C（receptor-type tyrosine-protein phosphatase C，PTPRC）等9個，其中分子對接較好的組合：維生素A與絲裂原活化蛋白激酶14（mitogen-activated protein kinase 14，MAPK14）、色氨酸與基質(zhì)金屬蛋白酶9（matrix metalloproteinase 9，MMP9）、維生素A與PTPRC。GO及KEGG分析結(jié)果顯示，主要涉及免疫炎癥、激素分泌、血液循環(huán)及MAPK通路等方面。鹿瓜多肽治療骨關(guān)節(jié)炎的患者主要為女性中老年人，其常見證型為肝腎虧虛、濕熱及氣滯血瘀證；其中藥治療多聯(lián)合補腎壯骨劑、活血化瘀劑以求標(biāo)本同治，而化學(xué)藥治療常對癥治療以緩解癥狀。通過靶點網(wǎng)絡(luò)構(gòu)建與分析則得出鹿瓜多肽主要通過維生素A及氨基酸對MMP9、MPO等靶點進(jìn)行調(diào)控從而發(fā)揮作用。

真實世界；靶點網(wǎng)絡(luò)；鹿瓜多肽注射液；人群特征；聯(lián)合用藥；骨關(guān)節(jié)炎；醫(yī)院信息系統(tǒng)

骨關(guān)節(jié)炎（osteoarthritis）是一種全身復(fù)雜因素導(dǎo)致的以軟骨脫失為特征的慢性退行性疾病，世界范圍內(nèi)發(fā)病人口約有3億[1-2]。骨關(guān)節(jié)炎是老年人殘疾的主要原因并可導(dǎo)致疼痛、功能障礙及生活質(zhì)量下降，因此對個人、家庭乃至社會都有較大的危害性[3-4]。目前骨關(guān)節(jié)炎的發(fā)病機制尚未得到充分研究，其藥物治療也缺乏行之有效的改善方案，故對其進(jìn)行研究并探索安全、有效的藥物具有較大的臨床意義[5]。

鹿瓜多肽注射液是由梅花鹿Temminck骨骼及甜瓜子提取而成的注射劑[6]，其含有骨誘導(dǎo)多肽類生物因子、有機鈣、磷等成分，具有調(diào)節(jié)骨代謝、刺激成骨細(xì)胞增殖、調(diào)節(jié)鈣磷代謝、抗炎、鎮(zhèn)痛等作用[7]，廣泛應(yīng)用于骨關(guān)節(jié)炎的治療[8]。為了更深入了解鹿瓜多肽注射液（簡稱鹿瓜多肽）在真實世界中治療骨關(guān)節(jié)炎的人群特征、療效、聯(lián)合用藥情況及可能的機制，本研究對中醫(yī)臨床基礎(chǔ)醫(yī)學(xué)研究所在全國26家三級甲等醫(yī)院建立的信息系統(tǒng)（hospital information systeam，HIS）中使用鹿瓜多肽治療骨關(guān)節(jié)炎患者的一般情況、治療結(jié)局及聯(lián)合用藥情況進(jìn)行分析，同時結(jié)合文獻(xiàn)及相關(guān)數(shù)據(jù)庫以“靶點”為基礎(chǔ)對鹿瓜多肽治療骨關(guān)節(jié)炎的作用機制進(jìn)行分析，以期為其臨床及進(jìn)一步研究提供參考。

1 真實世界研究

1.1 資料與方法

1.1.1 數(shù)據(jù)來源 真實世界數(shù)據(jù)來源于中國中醫(yī)科學(xué)院中醫(yī)臨床基礎(chǔ)醫(yī)學(xué)研究所建立的HIS數(shù)據(jù)倉庫[9]，納入時間為2003年7月—2019年12月，內(nèi)容主要包括患者的一般情況、用藥情況、診斷信息、醫(yī)囑信息、實驗室檢查信息等內(nèi)容[10]。本研究已通過中國中醫(yī)科學(xué)院中醫(yī)臨床基礎(chǔ)醫(yī)學(xué)研究所倫理委員會審批（編號2011-11）。

1.1.2 納入及排除標(biāo)準(zhǔn)

（1）納入標(biāo)準(zhǔn)：主要診斷為骨關(guān)節(jié)炎的患者；其中鹿瓜多肽組為醫(yī)囑單使用鹿瓜多肽次數(shù)不小于3次的骨關(guān)節(jié)炎患者；對照組為住院期間未使用鹿瓜多肽進(jìn)行治療的骨關(guān)節(jié)炎患者。

（2）排除標(biāo)準(zhǔn)：合并用藥為溶媒者（包括5%葡萄糖注射液、0.9%氯化鈉注射液、葡萄糖氯化鈉注射液3種）；合并用藥為氯化鉀和維生素者。

1.1.3 數(shù)據(jù)的標(biāo)準(zhǔn)化和規(guī)范化 為保證數(shù)據(jù)的標(biāo)準(zhǔn)及規(guī)范化，本研究對來自26家不同醫(yī)院的數(shù)據(jù)進(jìn)行了包括患者身份、診斷、藥物名稱、藥物分類及劑量等的標(biāo)準(zhǔn)化處理以保證患者的唯一性、各表之間的相關(guān)性、醫(yī)囑和診斷名稱的一致性、用藥劑量單位和理化指標(biāo)結(jié)果的有效性。其中，西醫(yī)診斷信息規(guī)范化參照《疾病和有關(guān)健康問題的國際統(tǒng)計分類（第10次修訂本）》（ICD-10）[11]；中醫(yī)診斷信息參照《中醫(yī)診斷學(xué)》進(jìn)行規(guī)范化[12]；藥物作用機制分類參照《新編藥物學(xué)》[13]。

1.1.4 數(shù)據(jù)的分析方法 采用SPSS v.18.0、SAS v.9.2對基本信息、中西醫(yī)診斷信息、醫(yī)囑信息進(jìn)行描述性統(tǒng)計，并用Excel 2007輔助做圖，聯(lián)合用藥分析采用Apriori算法建立模型[14]；使用Clementine v.12.0對合并藥物、并發(fā)癥進(jìn)行關(guān)聯(lián)分析；利用2檢驗、Logistic回歸比較應(yīng)用與未應(yīng)用鹿瓜多肽治療骨關(guān)節(jié)的療效差異，以＜0.05為差異有統(tǒng)計學(xué)意義。

1.2 結(jié)果

1.2.1 使用鹿瓜多肽的骨關(guān)節(jié)炎患者基本特征 ①性別：納入使用鹿瓜多肽的骨關(guān)節(jié)炎患者4050例，其中男性1200例（29.65%），女性2850例，（70.35%）；②年齡：46～65歲年齡段患者最多（2248例，55.51%），其次為66～75歲（850例，20.99%）；③職業(yè)：職業(yè)中不便分類的其他從業(yè)人員居多（2631例，65.24%）；④入院情況：入院科室主要集中在骨科（2085例，51.48%），病情多為一般（3678例，90.80%），入院方式主要是門診（3908例，96.50%），住院天數(shù)多為15～28 d（46.05%），住院患者多為醫(yī)?；颊咔移滟M用多處于10 001～20 000元。

1.2.2 使用鹿瓜多肽的骨關(guān)節(jié)炎患者合并的西醫(yī)疾病及中醫(yī)證候 4050例使用鹿瓜多肽的骨關(guān)節(jié)炎患者中，高血壓病是最常見的合并疾?。?88例，9.58%）；其次為骨折（281例，6.94%）、骨質(zhì)疏松癥（247例，6.10%）、腰椎間盤突出癥（227例，5.60%）、2型糖尿?。?82例，4.49%）。排除未錄入中醫(yī)證候信息的多數(shù)西醫(yī)院，從而對567例有中醫(yī)證候的骨關(guān)節(jié)炎患者進(jìn)行頻數(shù)分析，其中最常見者為肝腎虧虛證（220例，38.80%），其次為濕熱證（188例，33.16%）和氣滯血瘀證（11.99%）。

1.2.3 使用鹿瓜多肽的骨關(guān)節(jié)炎患者的用藥劑量、療程 與鹿瓜多肽說明書中標(biāo)明的劑量（靜脈滴注8～12 mL/d，10～15 d為一療程或遵醫(yī)囑）相比，4050例骨關(guān)節(jié)炎患者所使用的最小給藥量為4 mL/d，最高為36 mL/d，劑量范圍在8～12 mL最多（占比50.70%），超說明書劑量用藥1426例，占比38.40%；患者用藥天數(shù)最短1 d，最長28 d，平均值為3.47 d，用藥天數(shù)主要集中在1～3 d（占比79.70%）。

1.2.4 使用鹿瓜多肽的骨關(guān)節(jié)炎患者的療效 將4050例骨關(guān)節(jié)炎患者出院結(jié)局為治愈及好轉(zhuǎn)認(rèn)為有效，則有3964例使用鹿瓜多肽治療骨關(guān)節(jié)炎的患者為有效，總有效率為97.88%；選定其中病程記錄及隨訪完善，有明確記錄治療結(jié)局的病例3071例，同時按照與病例組年齡差不大于5歲、性別相同、入院病情近似進(jìn)行等量匹配而納入未使用鹿瓜多肽治療的骨關(guān)節(jié)炎患者3071例，以結(jié)局為治愈、好轉(zhuǎn)認(rèn)定為有效，其他、無效、死亡認(rèn)定為無效。按照重新劃分的2組進(jìn)行鹿瓜多肽組與對照組的組間比較，若不考慮混雜因素進(jìn)行分析，結(jié)果如表1所示，2組治療結(jié)局差異有統(tǒng)計學(xué)意義，即鹿瓜多肽組療效優(yōu)于對照組（＜0.05）。最后采用3種不同的Logistic回歸分析方法分析鹿瓜多肽組與對照組的治療結(jié)果差異，首先使用單變量的Logistic回歸進(jìn)行分析；然后通過對年齡、性別等人口學(xué)信息及入院病情、醫(yī)療費用、用藥、合并病癥情況等使用傾向評分對對照組病例進(jìn)行加權(quán)處理后，再使用單變量Logistic回歸進(jìn)行分析；最后選取病重天數(shù)為協(xié)變量，進(jìn)行傾向評分加權(quán)的Logistic回歸分析，結(jié)果如表2所示，3種Logistic回歸分析方法的結(jié)果均為＜0.05，綜合可知鹿瓜多肽組與對照組的治療結(jié)果存在差異，并且鹿瓜多肽組的治療效果更佳。

表1 鹿瓜多肽組與對照組骨關(guān)節(jié)炎患者結(jié)局比較

Table 1 Comparison of outcomes of osteoarthritis patients in Cervus and Cucumis Polypeptide Injection group and control group

結(jié)局鹿瓜多肽組（n = 3071）對照組（n = 3071）χ2P 無效/例 48 298189.8910.000 有效/例30232773

1.2.5 使用鹿瓜多肽的骨關(guān)節(jié)炎患者常用合并藥及類別 鹿瓜多肽的聯(lián)用藥物中頻數(shù)前5位的藥物均為化學(xué)藥，包括利多卡因注射液、玻璃酸鈉等；聯(lián)合應(yīng)用最多的中藥則有金天格膠囊、強骨膠囊等5種，如表3所示。聯(lián)合應(yīng)用的化學(xué)藥按照其藥理作用進(jìn)行分類，可知前10類為抗生素類、解熱鎮(zhèn)痛藥等；聯(lián)合應(yīng)用的中藥類別則包括活血化瘀劑、補腎壯骨劑等，具體見表4。

1.2.6 鹿瓜多肽的聯(lián)合用藥關(guān)聯(lián)規(guī)則分析 根據(jù)支持度分別選取前5位鹿瓜多肽最常聯(lián)用的2種及3種的化學(xué)藥和中藥進(jìn)行關(guān)聯(lián)分析，可見其最常聯(lián)合使用的化學(xué)藥是苯巴比妥鈉＋頭孢唑啉（支持度15.97%）及苯巴比妥鈉、骨肽注射劑＋頭孢唑啉（支持度7.40%）；中藥是脈血康膠囊＋金天格膠囊（支持度7.03%）及接骨七厘片、脈血康膠囊＋金天格膠囊（支持度3.72%），見表5。

表2 3種Logistic回歸方法對鹿瓜多肽組及對照組治療骨關(guān)節(jié)炎效果的分析

Table 2 Analysis of effect of three Logistic regression methods on treatment of osteoarthritis in Cervus and Cucumis Polypeptide Injection group and control group

分析方法回歸系數(shù)P值 單變量Logistic回歸1.912 112.00×10?16 傾向評分加權(quán)的單變量Logistic回歸1.663 182.00×10?16 帶協(xié)變量（病重天數(shù)）的傾向評分加權(quán)的Logistic回歸1.663 092.00×10?16

表3 使用鹿瓜多肽的骨關(guān)節(jié)炎患者聯(lián)合使用頻率前5的化學(xué)藥和中藥

Table 3 Top 5 western and traditional Chinese medicines used in combination with Cervus and Cucumis Polypeptide Injection in osteoarthritis patients

化學(xué)藥中藥 名稱n/例占比/%名稱n/例占比/% 利多卡因注射液113928.12金天格膠囊54413.43 玻璃酸鈉 97224.00強骨膠囊45911.33 骨肽注射劑 80119.78白芍總苷膠囊42010.37 倍他米松 77019.01紅花注射液397 9.80 七葉皂苷鈉注射液 68616.94血栓通注射液372 9.19

表4 使用鹿瓜多肽的骨關(guān)節(jié)炎患者聯(lián)合應(yīng)用的中藥和化學(xué)藥類別(前10位)

Table 4 Types of traditional Chinese and western medicines used in combination with Cervus and Cucumis Polypeptide injection in osteoarthritis patients (top 10)

化學(xué)藥中藥 類別n/例占比/%類別n/例占比/% 抗生素類259964.17活血化瘀劑256563.33 抗炎鎮(zhèn)痛藥193947.88補腎壯骨劑126431.21 營養(yǎng)藥157838.96祛風(fēng)除濕劑101825.14 激素151937.51清熱解毒劑 78219.31 鎮(zhèn)痛藥140334.64清熱利濕劑 67216.59 麻醉藥140234.62清熱化痰劑 340 8.40 抗風(fēng)濕藥137834.02益氣養(yǎng)陰劑 218 5.38 抗酸及抗?jié)兯?35833.53瀉熱通便劑 187 4.62 降壓藥110227.21潤腸通便劑 179 4.42 抗骨質(zhì)疏松藥 94023.21益氣養(yǎng)血劑 172 4.25

1.2.7 鹿瓜多肽合并中藥和化學(xué)藥類別關(guān)聯(lián)規(guī)則 按照支持度順序分別選取鹿瓜多肽聯(lián)合用藥前5位的藥物種類組合進(jìn)行分析，可得鹿瓜多肽最常合并化學(xué)藥種類為鎮(zhèn)痛藥＋抗生素（支持度50.40%）；中藥為補腎壯骨劑＋活血化瘀劑（支持度21.053%），見表6。

表5 鹿瓜多肽聯(lián)合使用中藥和化學(xué)藥的關(guān)聯(lián)分析(前5位)

表6 鹿瓜多肽合并中藥和化學(xué)藥藥物類別關(guān)聯(lián)分析(前5位)

Table 6 Association analysis of drug class of Cervus and Cucumis Polypeptide injection combined with traditional Chinese and western medicine (top 5)

合并藥物類別支持度/%置信度/%提升度 中藥類別補腎壯骨劑＋活血化瘀劑21.05369.61.320 清熱解毒劑＋活血化瘀劑 6.85065.71.240 瀉熱通便劑＋活血化瘀劑 4.98660.11.140 補腎壯骨劑＋清熱解毒劑 4.37914.51.390 補腎壯骨劑＋瀉熱通便劑 4.23914.01.690 化學(xué)藥類別鎮(zhèn)痛藥＋抗生素50.40094.81.037 高滲脫水藥＋抗生素43.40094.01.029 營養(yǎng)藥＋抗生素33.40095.61.046 抗酸及抗?jié)兯帲股?2.50096.51.055 高滲脫水藥＋鎮(zhèn)痛藥28.50061.91.164

1.2.8 與鹿瓜多肽聯(lián)合應(yīng)用的中藥、化學(xué)藥關(guān)聯(lián)分析 采用關(guān)聯(lián)分析的方法并利用Clemetine 12.0 web軟件做圖，對所有與鹿瓜多肽合并用藥進(jìn)行關(guān)聯(lián)分析，以求更直觀地表現(xiàn)鹿瓜多肽與多種藥物聯(lián)合應(yīng)用的情況。分析結(jié)果以連線呈現(xiàn)，連線越粗表示聯(lián)合使用的頻次越多、關(guān)系越密切，即同時出現(xiàn)的幾率最大。從圖1可知，常見合并用藥是倍他米松、玻璃酸鈉、骨肽注射劑、利多卡因注射液、七葉皂苷鈉注射液、紅花注射液、金天格膠囊。圖2則顯示與鹿瓜多肽聯(lián)合應(yīng)用的藥物主要類別中化學(xué)藥為補鈣藥、解熱鎮(zhèn)痛藥、抗骨質(zhì)疏松藥、抗生素、營養(yǎng)藥，中藥為補腎壯骨劑、活血化瘀劑。

2 靶點網(wǎng)絡(luò)預(yù)測

2.1 資料與方法

2.1.1 鹿瓜多肽有效成分獲取與靶點收集 利用CNKI（https://www.cnki.net/）[15-17]及TCMSP數(shù)據(jù)庫（http://tcmspw.com）獲取鹿瓜多肽主要組成藥材梅花鹿骨及甜瓜子的主要有效成分，然后使用SwissADME平臺（http://www.swissadme.ch/）限定有效成分的類藥性（drug likeness）評價不小于3個項目為“yes”，并將篩選出的有效成分導(dǎo)入SwissTargetPrediction平臺（http://www. swisstargetprediction.ch/）按照預(yù)測可能性大于0獲得藥物靶點。

圖1 與鹿瓜多肽聯(lián)合應(yīng)用藥物關(guān)聯(lián)分析

圖2 與鹿瓜多肽聯(lián)合應(yīng)用藥物類別關(guān)聯(lián)分析

2.1.2 骨關(guān)節(jié)炎疾病靶點收集 首先在GEO數(shù)據(jù)庫（https://www.ncbi.nlm.nih.gov/geo/）中檢索骨關(guān)節(jié)炎并限定數(shù)據(jù)類型（expression profiling by array）及物種（Homo sapiens）以獲得骨關(guān)節(jié)炎基因表達(dá)數(shù)據(jù)集，將之使用Perl及R語言進(jìn)行差異分析，從而獲得正常組及骨關(guān)節(jié)炎組存在差異的骨關(guān)節(jié)炎相關(guān)基因；然后在Genecards數(shù)據(jù)庫（https://www. genecards.org/）中檢索骨關(guān)節(jié)炎并限定相關(guān)性大于中位數(shù)獲取骨關(guān)節(jié)炎相關(guān)基因；分別在OMIM數(shù)據(jù)庫（https://omim.org/）及TTD數(shù)據(jù)庫（http://db. idrblab.net/ttd/）中檢索骨關(guān)節(jié)炎獲取骨關(guān)節(jié)炎相關(guān)基因；再將以上4個數(shù)據(jù)庫所得骨關(guān)節(jié)炎相關(guān)基因進(jìn)行合并、去重并與Unprot（http://eta.uniprot.org/）據(jù)庫中已驗證的人類基因相匹配，而得出骨關(guān)節(jié)炎疾病靶點。

2.1.3 鹿瓜多肽靶點與骨關(guān)節(jié)炎靶點交集獲得及“藥物-疾病-有效成分-靶點”網(wǎng)絡(luò)和蛋白相互作用（protein-protein interaction，PPI）網(wǎng)絡(luò)構(gòu)建 將“2.1.1”所得鹿瓜多肽有效成分靶點與“2.1.2”所得骨關(guān)節(jié)炎疾病靶點取交集獲得藥物直接作用于疾病的靶點；使用Cytoscape v.3.8.0以鹿瓜多肽、骨關(guān)節(jié)炎、有效成分、靶點為節(jié)點，其相互作用關(guān)系為邊構(gòu)建網(wǎng)絡(luò)并篩選度（degree）值不小于6倍中位數(shù)的有效成分作為核心成分，然后利用MCODE插件對交集靶點在String數(shù)據(jù)庫（https://string-db.org/）所得的PPI網(wǎng)絡(luò)進(jìn)行聚類分析并根據(jù)度值篩選每個聚類中大于中位數(shù)的靶點作為核心靶點。

2.1.4 基因本體（gene ontology，GO）富集分析和京都基因及基因組百科全書（Kyoto encyclopedia of genes and genomes，KEGG）通路富集分析 將“2.1.3”所得核心靶點導(dǎo)入Metascape平臺（https:// metascape.org/）分別進(jìn)行GO生物過程、分子功能、細(xì)胞成分分析及KEGG通路富集分析，再將結(jié)果使用易漢博平臺（http://www.ehbio.com/ImageGP/ index.php）繪制富集氣泡圖。

2.1.5 核心成分與核心靶點分子對接 將“2.1.3”中所得核心成分與核心靶點分別在Zinc數(shù)據(jù)庫（https://zinc.docking.org/）及Uniprot數(shù)據(jù)庫中查找其3D結(jié)構(gòu)，然后使用Autodock Vina進(jìn)行分子對接，其對接較好的結(jié)果則使用Pymol進(jìn)行可視化。

2.2 結(jié)果

2.2.1 鹿瓜多肽注射液有效成分與靶點 從CNKI及TCMSP數(shù)據(jù)庫收集到梅花鹿骨有效成分44種，相關(guān)靶點169個；甜瓜子有效成分33種，相關(guān)靶點182個，其具體信息見表7。

表7 鹿瓜多肽注射液有效成分信息

Table 7 Information of effective ingredients of Cervus and Cucumis Polypeptide Injection

有效成分類藥性SMILES來源 酪氨酸4C1=CC(=CC=C1CC(C(=O)O)N)O梅花鹿骨、甜瓜子 谷氨酸3C(CC(=O)O)C(C(=O)O)N梅花鹿骨、甜瓜子 蘇氨酸3CC(C(C(=O)O)N)O梅花鹿骨、甜瓜子 甘氨酸3C(C(=O)O)N梅花鹿骨、甜瓜子 絲氨酸3C(C(C(=O)O)N)O梅花鹿骨、甜瓜子 天冬氨酸3C(C(C(=O)O)N)C(=O)O梅花鹿骨、甜瓜子 異亮氨酸3CCC(C)C(C(=O)O)N梅花鹿骨、甜瓜子 賴氨酸3C(CCN)CC(C(=O)O)N梅花鹿骨、甜瓜子 苯丙氨酸4C1=CC=C(C=C1)CC(C(=O)O)N梅花鹿骨、甜瓜子 亮氨酸3CC(C)CC(C(=O)O)N梅花鹿骨、甜瓜子 精氨酸3C(CC(C(=O)O)N)CN=C(N)N梅花鹿骨、甜瓜子 脯氨酸3C1CC(NC1)C(=O)O梅花鹿骨、甜瓜子 丙氨酸3CC(C(=O)O)N梅花鹿骨、甜瓜子 纈氨酸3CC(C)C(C(=O)O)N梅花鹿骨、甜瓜子 胱氨酸1C(C(C(=O)O)N)SSCC(C(=O)O)N梅花鹿骨、甜瓜子 甲硫氨酸3CSCCC(C(=O)O)N梅花鹿骨 組氨酸3C1=C(NC=N1)CC(C(=O)O)N梅花鹿骨、甜瓜子 天冬氨酸3C(C(C(=O)O)N)C(=O)N梅花鹿骨 色氨酸5C1=CC=C2C(=C1)C(=CN2)CC(C(=O)O)N梅花鹿骨 羥脯氨酸3C1C(CNC1C(=O)O)O梅花鹿骨 蛋氨酸3CSCCC(C(=O)O)N甜瓜子 鈣N/AN/A梅花鹿骨、甜瓜子 鎂N/AN/A梅花鹿骨、甜瓜子 鉀N/AN/A梅花鹿骨、甜瓜子 鈉N/AN/A梅花鹿骨、甜瓜子 磷N/AN/A梅花鹿骨、甜瓜子 鐵N/AN/A梅花鹿骨、甜瓜子 銅N/AN/A梅花鹿骨、甜瓜子 鎳N/AN/A梅花鹿骨 鈷N/AN/A梅花鹿骨 鋅N/AN/A梅花鹿骨、甜瓜子 錳N/AN/A梅花鹿骨 硅N/AN/A梅花鹿骨 鈦N/AN/A梅花鹿骨 鎵N/AN/A梅花鹿骨 鍶N/AN/A梅花鹿骨 銫N/AN/A梅花鹿骨 軟骨素0CC(=O)NC1C(C(C(OC1O)CO)O)OC2C(C(C=C(O2)C(=O)O)O)O梅花鹿骨 維生素A4CC1=C(C(CCC1)(C)C)C=CC(=CC=CC(=CCO)C)C梅花鹿骨 維生素B21CC1=CC2=C(C=C1C)N(C3=NC(=O)NC(=O)C3=N2)CC(C(C(CO)O)O)O梅花鹿骨

續(xù)表7

SMILES-簡化分子線性輸入規(guī)范 N/A表示不存在類藥性及SMILES號的金屬成分的相關(guān)項目無法獲得

SMILES- simplified molecular input line entry system N/A indicates that related items of the metal components with no drug-likeness and SMILES numbers are not applicable

2.2.2 骨關(guān)節(jié)炎疾病靶點 在GEO數(shù)據(jù)庫檢索獲得骨關(guān)節(jié)炎相關(guān)數(shù)據(jù)集為GSE169077，其中包括正常樣品5例、骨關(guān)節(jié)炎樣品6例，差異分析則獲得骨關(guān)節(jié)炎相關(guān)基因345個，見圖3；Genecards、OMIM、TTD數(shù)據(jù)庫則分別獲得骨關(guān)節(jié)炎相關(guān)基因2366、6、23個，將以上數(shù)據(jù)整理后一共獲得1964個骨關(guān)節(jié)炎疾病靶點。

2.2.3 鹿瓜多肽注射液與骨關(guān)節(jié)炎交集靶點、“藥物-疾病-有效成分-靶點”網(wǎng)絡(luò)及PPI網(wǎng)絡(luò)構(gòu)建 將鹿瓜多肽注射液靶點與骨關(guān)節(jié)炎靶點取交集，如圖4所示，可得60個鹿瓜多肽直接作用于骨關(guān)節(jié)炎的靶點；通過構(gòu)建“鹿瓜多肽注射液-骨關(guān)節(jié)炎-有效成分-靶點”網(wǎng)絡(luò)并進(jìn)行網(wǎng)絡(luò)拓?fù)鋵W(xué)分析可得核心有效成分為棕櫚酸、肉豆蔻酸、色氨酸、維生素A、亮氨酸、谷氨酸、異亮氨酸、苯丙氨酸、賴氨酸、酪氨酸10種，見圖5；通過構(gòu)建PPI網(wǎng)絡(luò)并使用MCODE插件進(jìn)行聚類可分為3個聚類分型，再擇其度值較大者作為核心靶點，如圖6所示，包括過氧化物酶體增殖物激活受體γ（peroxisome proliferator-activated receptor gamma，PPARG）、髓過氧化物酶（myeloperoxidase，MPO）、受體型酪氨酸蛋白磷酸酶C（receptor-type tyrosine-protein phosphatase C，PTPRC）、絲裂原活化蛋白激酶1（mitogen-activated protein kinase 1，MAPK1）、絲裂原活化蛋白激酶14（mitogen-activated protein kinase 14，MAPK14）、雌激素受體β（ESR2）、基質(zhì)金屬蛋白酶-9（matrix metalloproteinase-9，MMP9）、前列腺素E合酶（prostaglandin E synthase，PTGES）、血管緊張素轉(zhuǎn)換酶2（angiotensin-converting enzyme，ACE）9個靶點。

圖4 鹿瓜多肽注射液與骨關(guān)節(jié)炎交集靶點韋恩圖

2.2.4 交集靶點GO富集分析和KEGG通路富集分析 通過將交集靶點導(dǎo)入Metascape平臺進(jìn)行GO分析，可得GO生物功能主要有細(xì)胞化學(xué)穩(wěn)態(tài)、血液循環(huán)的系統(tǒng)過程、激素應(yīng)答、免疫應(yīng)答調(diào)節(jié)、離子轉(zhuǎn)運調(diào)節(jié)、分泌調(diào)節(jié)、MAPK級聯(lián)途徑調(diào)節(jié)等；GO分子功能主要有氧化還原酶活性、肽酶活性、氨基結(jié)合、核受體活性等；GO細(xì)胞成分主要有膜側(cè)面、裂解空泡、樹突等；KEGG通路富集分析結(jié)果顯示，主要涉及癌癥中的通路、神經(jīng)活性配體-受體相互作用、炎癥介導(dǎo)的TRP通路調(diào)節(jié)、腫瘤壞死因子信號通路等，如圖7、8所示。

圖5 “鹿瓜多肽注射液-骨關(guān)節(jié)炎-有效成分-靶點”網(wǎng)絡(luò)

2.2.5 核心成分與核心靶點分子對接 將網(wǎng)絡(luò)分析中的核心成分與核心靶點進(jìn)行分子對接，可知多數(shù)核心成分與核心靶點間可形成穩(wěn)定結(jié)構(gòu)（結(jié)合能＜?20.93 kJ/mol）[18]。其中維生素A、色氨酸、酪氨酸及苯丙氨酸為對接較好的成分，MMP9、MPO、ESR2、MAPK1為對接較好的靶點，見圖9。其中對接最好的3個組合（維生素A與MAPK14、色氨酸與MMP9、維生素A與PTPRC）的分子對接構(gòu)象見圖10，由圖可知其多由氫鍵或疏水口袋結(jié)合。

3 討論

3.1 骨關(guān)節(jié)炎患者的人群特征差異

本研究結(jié)果顯示骨關(guān)節(jié)炎患者年齡多集中在46～65歲及66～75歲的中老年階段，與既往報道一致[19]。研究證實年齡是骨關(guān)節(jié)炎的獨立風(fēng)險因素[20]，其影響包括長時間的磨損、細(xì)胞外基質(zhì)隨年齡變化而導(dǎo)致的關(guān)節(jié)硬度增加、神經(jīng)肌肉的衰退、激素狀態(tài)變化、活動水平及飲食等方面[21]；而在細(xì)胞水平，細(xì)胞應(yīng)激與細(xì)胞損傷會導(dǎo)致慢性炎癥的發(fā)生，從而促進(jìn)與衰老相關(guān)的疾病，在骨關(guān)節(jié)炎中則通常表現(xiàn)為線粒體功能障礙及活性氧的產(chǎn)生、細(xì)胞對生長因子的反應(yīng)性降低、自噬與細(xì)胞死亡等[22]。同時本研究發(fā)現(xiàn)骨關(guān)節(jié)炎患者中女性占大多數(shù)（70.35%）且絕經(jīng)后女性骨關(guān)節(jié)炎發(fā)病率急劇增高，與流行病學(xué)報道相符[23]。研究認(rèn)為雌激素是骨關(guān)節(jié)炎的重要調(diào)節(jié)劑，可通過與細(xì)胞生長因子、黏附分子和細(xì)胞因子相互作用，在多個層面影響軟骨細(xì)胞功能[24]；但使用雌激素替代療法對骨關(guān)節(jié)炎的治療卻證據(jù)不足[25]。4050例骨關(guān)節(jié)炎患者中，高血壓病是最常見的合并病，多項研究表明高血壓與骨關(guān)節(jié)炎存在流行病學(xué)及分子生物學(xué)的相關(guān)性[26]；且二者與肥胖、高血脂、高血糖等代謝障礙常伴發(fā)[27]，該情況可能因類似的生活方式、相同的危險因素以及代謝紊亂患者全身的炎癥反應(yīng)相關(guān)[28]。除此之外的骨折、骨質(zhì)疏松癥、腰椎間盤突出癥等骨關(guān)節(jié)炎常見合并病則多與年齡及激素代謝所致的退行性變化相關(guān)。

圖6 交集靶點PPI網(wǎng)絡(luò)及核心靶點篩選

圖7 交集靶點GO富集分析

圖8 交集靶點KEGG通路富集分析

圖9 核心成分與核心靶點分子對接熱圖

圖10 對接較好的核心成分與核心靶點分子對接構(gòu)象圖

骨關(guān)節(jié)炎屬祖國醫(yī)學(xué)“骨痹”范疇，“骨痹”始載于《黃帝內(nèi)經(jīng)》的《素問·痹論篇》：“風(fēng)寒濕三氣雜至，合而為痹也?！远龃苏邽楣潜浴盵29]，可知骨痹多由風(fēng)、寒、濕等邪流駐肌腠經(jīng)絡(luò)，痹阻經(jīng)脈關(guān)節(jié)所致?！侗哉撈酚钟校骸肮潜圆灰?，復(fù)感于邪，內(nèi)舍于腎”[29]，認(rèn)為骨痹發(fā)病存在2個病理過程，一為“病久”，即久病正氣受戕，病邪乘虛而入；二為“復(fù)感”，即反復(fù)感受風(fēng)、寒、濕等外邪，邪盛而正虛?！端貑枴らL刺節(jié)論》則載：“病在骨，骨重不可舉，骨髓酸痛，寒氣至，名曰骨痹”[29]，可見骨痹主要癥狀為關(guān)節(jié)活動障礙及疼痛。通過對4050例骨關(guān)節(jié)炎患者分析發(fā)現(xiàn)最常見的中醫(yī)證候為肝腎虧虛證、濕熱證及氣滯血瘀證。因骨痹多發(fā)于老年人，其腎精多虛；女性多發(fā)，其先天之本常受累；且影響及筋骨，二者分屬肝腎所統(tǒng)，故肝腎虧虛多為骨痹之本?！毒霸廊珪诰雰?nèi)傷》曰：“脾主四肢，而勞倦過度則脾氣傷矣”[30]，因此骨痹亦與脾有較大關(guān)系，勞倦傷脾，釀生痰濕，痰濕生熱，變生慢性炎癥，膠結(jié)難解，則為骨痹。清代葉天士云：“初病脹痛無形，久則形堅似梗，是初為氣結(jié)在經(jīng)，久則血傷入絡(luò)”[31]，故骨痹中氣滯血瘀癥亦常見，且常與高血壓等心血管疾病并見。

以上分析可知，老年人尤其是女性在45歲以上的年齡應(yīng)當(dāng)對骨關(guān)節(jié)炎進(jìn)行預(yù)防，主要包括：注意飲食、控制體質(zhì)量；合理運動、避免損傷關(guān)節(jié)；控制基礎(chǔ)疾病、改變不良生活方式等；對于有相關(guān)癥狀的則要及時就醫(yī)，避免癥狀加重而造成功能障礙及日常生活能力下降；而對于肝腎不足、氣滯血瘀等體質(zhì)問題則可適當(dāng)采用中藥調(diào)節(jié)。

3.2 鹿瓜多肽注射液療效分析及聯(lián)合用藥分析

通過分析可見，鹿瓜多肽治療骨關(guān)節(jié)炎具有較高的總有效率，但仍存在少量患者療效不佳的情況，因此在使用鹿瓜多肽的同時應(yīng)當(dāng)注意辨證準(zhǔn)確并且適當(dāng)聯(lián)合藥物以進(jìn)一步提高療效；鑒于臨床上對鹿瓜多肽的使用多是聯(lián)合用藥，并且應(yīng)用的骨關(guān)節(jié)炎患者在年齡、性別、發(fā)病部位、合并病、住院時間、住院費用及治療情況等其他因素存在一定的差異，因此需排除其他因素，故本研究采用結(jié)局指標(biāo)結(jié)合傾向性分析排除混雜因素并構(gòu)建Logistic回歸模型以對鹿瓜多肽治療骨關(guān)節(jié)炎的療效進(jìn)行較準(zhǔn)確的判斷。從結(jié)果來看，各個組別均為鹿瓜多肽組療效較好，故臨床應(yīng)用鹿瓜多肽可提升骨關(guān)節(jié)炎的療效，應(yīng)予推廣。

鹿瓜多肽包含梅花鹿骨及甜瓜子兩藥，前者補虛羸、強筋骨[32]，后者化瘀、排膿、療傷止痛[33]，二者合用兼顧骨關(guān)節(jié)炎患者標(biāo)本兩端，并可對其疼痛癥狀及高血壓、骨折等常見兼證進(jìn)行治療，凸顯了中醫(yī)整體論治的特點。根據(jù)研究結(jié)果可知，骨關(guān)節(jié)炎患者使用鹿瓜多肽劑量約半數(shù)處于正常范圍，但仍存在較多超量用藥的情況；此外，患者的治療周期多為1～3 d，與說明書療程10～15 d差異較大，因此在使用該藥時應(yīng)當(dāng)注意合理用藥劑量以避免不良反應(yīng)并適當(dāng)延長用藥周期以更大程度發(fā)揮藥物療效并改善患者癥狀。根據(jù)鹿瓜多肽的聯(lián)合用藥分析可知，骨關(guān)節(jié)炎患者臨床常合并中藥補腎壯骨劑（金天格膠囊）、活血化瘀劑（脈血康膠囊）及化學(xué)藥抗炎鎮(zhèn)痛類（骨肽注射劑）藥物，前兩者分別針對肝腎虧虛之本以及氣滯血瘀之標(biāo)，而后者則針對骨關(guān)節(jié)炎表現(xiàn)出的疼痛及炎癥進(jìn)行對癥治療。另外部分骨關(guān)節(jié)炎患者由于關(guān)節(jié)腔滑液的缺失而出現(xiàn)相應(yīng)的臨床癥狀，玻璃酸鈉作為關(guān)節(jié)滑液的主要成分，可改善關(guān)節(jié)功能、緩解疼痛，但同時需指出玻璃酸鈉在軟骨保護(hù)和延緩骨關(guān)節(jié)炎疾病進(jìn)程中的作用尚存爭議[34]。除了對骨關(guān)節(jié)炎的患者以上的治療外，臨床常采用關(guān)節(jié)腔封閉治療，即利多卡因注射液、糖皮質(zhì)激素（倍他米松）。利多卡因注射液麻醉持續(xù)時間長、穿透性強，主要用于阻滯麻醉和硬膜外麻醉。倍他米松具有抗炎作用，能抑制炎癥，減輕充血，兩者配合使用會達(dá)到抗炎鎮(zhèn)痛的效果。從上可知化學(xué)藥治療主要針對疼痛、炎癥等臨床癥狀，而中藥則擅長標(biāo)本同治，兼顧肝腎之虛與氣滯血瘀之實，故臨床在使用鹿瓜多肽治療骨關(guān)節(jié)炎時應(yīng)當(dāng)根據(jù)患者情況進(jìn)行聯(lián)合用藥：若患者關(guān)節(jié)疼痛、炎癥癥狀較劇，可聯(lián)合化學(xué)藥的麻醉劑、抗生素、激素類藥物迅速緩解癥狀；反之，則可主要聯(lián)用中藥以調(diào)節(jié)骨關(guān)節(jié)炎患者證候并改善其癥狀。

3.3 鹿瓜多肽治療骨關(guān)節(jié)炎機制分析

根據(jù)文獻(xiàn)收集到的靶點可知鹿瓜多肽的主要成分包括金屬離子類、氨基酸類、維生素類及脂肪酸類等，而梅花鹿骨及甜瓜子存在較多的相同成分，因此可見二者配伍主要表現(xiàn)為協(xié)同作用。從靶點網(wǎng)絡(luò)分析所得的核心成分來看：棕櫚酸被體外實驗證實可抑制軟骨破壞從而具有改善骨關(guān)節(jié)炎的作用[35]；肉豆蔻酸在動物實驗中顯示出改善創(chuàng)傷骨關(guān)節(jié)炎大鼠疼痛癥狀的作用[36]；色氨酸、谷氨酸等功能性氨基酸被認(rèn)為因具有抗炎、抗氧化作用而在骨關(guān)節(jié)炎治療中存在多種益處[37]；維生素A血清水平則與髖關(guān)節(jié)骨關(guān)節(jié)炎發(fā)病呈負(fù)相關(guān)[38]。從PPI網(wǎng)絡(luò)所得的核心靶點來看：PPARG在實驗中顯示其表達(dá)增高可調(diào)節(jié)組織修復(fù)與炎癥消退[39]；MPO則可作為早期骨關(guān)節(jié)炎的診斷標(biāo)志物而存在于骨關(guān)節(jié)炎患者的滑膜液中[40]；細(xì)胞實驗表明PTPRC、MAPK14在骨關(guān)節(jié)炎軟骨細(xì)胞中表達(dá)顯著升高并可能在骨關(guān)節(jié)炎治療中發(fā)揮重要作用[41]；ESR2暫無研究表明其與骨關(guān)節(jié)炎之間存在直接關(guān)聯(lián)，但鑒于雌激素與骨關(guān)節(jié)炎之間在多方面存在相關(guān)性，因此其作用仍待進(jìn)一步研究[42]；MMP9在大鼠實驗中被證實其可抑制軟骨分化并促進(jìn)軟骨細(xì)胞凋亡而促進(jìn)骨關(guān)節(jié)炎加重[43]；MAPK1則在體外實驗中被證實通過促進(jìn)其高表達(dá)治療骨關(guān)節(jié)炎[44]；PTGES有報道稱其在骨關(guān)節(jié)炎發(fā)病機制中起重要作用[45]；研究顯示，ACE的表達(dá)與骨關(guān)節(jié)炎之間存在明顯相關(guān)性，其中腺嘌呤胞嘧啶單倍型是保護(hù)因素，鳥嘌呤胸腺嘧啶單倍型是危險因素[46]。GO及KEGG分析表明，除免疫炎癥、激素分泌、血液循環(huán)及MAPK通路等機制外，癌癥與骨關(guān)節(jié)炎同樣作為與衰老相關(guān)的疾病，在骨更新、血管生成和神經(jīng)生長等方面具有一定的相似性[47]。然而盡管有相關(guān)回顧性研究表明痛風(fēng)患者較之骨關(guān)節(jié)炎患者結(jié)直腸癌患病率低，其具體相關(guān)性仍有待進(jìn)一步研究[48]。此外，骨關(guān)節(jié)炎作為一種疼痛相關(guān)疾病，其涉及外周和中樞神經(jīng)機制[49]，因此KEGG通路分析中涉及相關(guān)機制。分子對接則進(jìn)一步限定鹿瓜多肽治療骨關(guān)節(jié)炎的可能機制應(yīng)主要涉及維生素A及氨基酸對MMP9、MPO、ESR2、MAPK1等靶點進(jìn)行調(diào)控從而發(fā)揮作用。

4 結(jié)論

本研究通過對4050例使用鹿瓜多肽治療骨關(guān)節(jié)炎的患者進(jìn)行人群特征、療效與聯(lián)合用藥分析，從而得出骨關(guān)節(jié)炎患者主要為女性中老年人，其常見合并病包括高血壓、糖尿病等代謝類疾病及骨折、腰椎間盤突出等骨關(guān)節(jié)疾病，其常見證型為肝腎虧虛、濕熱及氣滯血瘀證；其中藥治療多聯(lián)合給予補腎壯骨、活血化瘀的中成藥以求標(biāo)本同治，而化學(xué)藥治療常聯(lián)合抗生素、激素、麻醉藥、解熱鎮(zhèn)痛藥以快速緩解癥狀。而應(yīng)用鹿瓜多肽的骨關(guān)節(jié)炎患者療效較之未使用者更好，并且具有較高的總有效率。通過對鹿瓜多肽靶點網(wǎng)絡(luò)進(jìn)行分析得出其主要通過維生素A及氨基酸作用于MMP9、MPO、ESR2、MAPK1等靶點對免疫炎癥、激素分泌、血液循環(huán)及MAPK通路進(jìn)行調(diào)控從而發(fā)揮作用。該結(jié)論可為鹿瓜多肽進(jìn)一步研究及臨床應(yīng)用提供參考。

利益沖突 所有作者均聲明不存在利益沖突

[1] Das S K, Farooqi A. Osteoarthritis [J]., 2008, 22(4): 657-675.

[2] GBD Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990—2017: A systematic analysis for the Global Burden of Disease Study 2017 [J]., 2018, 392(10159): 1789-1858.

[3] Abramoff B, Caldera F E. Osteoarthritis: Pathology, diagnosis, and treatment options [J]., 2020, 104(2): 293-311.

[4] Glyn-Jones S, Palmer A J R, Agricola R,. Osteoarthritis [J]., 2015, 386(9991): 376-387.

[5] Sacitharan P K. Ageing and osteoarthritis [J]., 2019, 91: 123-159.

[6] 曾新雨, 王志飛, 黎元元, 等. 真實世界中鹿瓜多肽注射液對骨折人群腎功能的影響實效研究 [J]. 世界中醫(yī)藥, 2022, 17(2): 255-260.

[7] 彭昊, 汪喆, 李章華, 等. 鹿瓜多肽注射液促進(jìn)骨折愈合過程中血管內(nèi)皮生長因子表達(dá)的機制 [J]. 中國臨床康復(fù), 2006, 10(37): 81-83.

[8] 孫帥玲, 謝雁鳴, 黎元元, 等. 基于真實世界醫(yī)院使用鹿瓜多肽注射液患者臨床用藥特征分析 [J]. 中國中藥雜志, 2016, 41(21): 4060-4065.

[9] 莊嚴(yán), 謝邦鐵, 翁盛鑫, 等. 中藥上市后再評價HIS“真實世界”集成數(shù)據(jù)倉庫的構(gòu)建與實現(xiàn) [J]. 中國中藥雜志, 2011, 36(20): 2883-2887.

[10] 孫帥玲, 謝雁鳴, 楊靖, 等. 基于真實世界子宮內(nèi)膜異位癥患者臨床特征與用藥規(guī)律初探 [J]. 時珍國醫(yī)國藥, 2015, 26(12): 3064-3067.

[11] 北京協(xié)和醫(yī)院世界衛(wèi)生組織國際分類家族合作中心. 疾病和有關(guān)健康問題的國際統(tǒng)計分類(第3卷) [M]. 第2版第10次修訂本. 北京: 人民衛(wèi)生出版社, 2008: 507-545.

[12] 朱文鋒. 中醫(yī)診斷學(xué) [M]. 北京: 中國中醫(yī)藥出版社, 2002: 508-758.

[13] 傅春生. 新編藥物學(xué) [M]. 天津: 天津科學(xué)技術(shù)出版社, 2017: 1-424.

[14] 崔貫勛, 李梁, 王柯柯, 等. 關(guān)聯(lián)規(guī)則挖掘中Apriori算法的研究與改進(jìn) [J]. 計算機應(yīng)用, 2010, 30(11): 2952-2955.

[15] 張敏, 陳曉燕. 鹿骨粉制作工藝及營養(yǎng)成分分析 [J]. 食品科技, 2004, 29(1): 95-97.

[16] 李洪玲, 但建明, 王琪, 等. 甜瓜籽營養(yǎng)成分分析 [J]. 中國油脂, 2004, 29(10): 68-69.

[17] 張慧君, 李虎, 惠林沖, 等.不同甜瓜種子性狀及脂肪酸成分分析[J]. 中國油脂, 2022, doi: 10.19902/j.cnki. zgyz.1003-7969.210391.

[18] 劉毅, 王瑞敏, 彭雪峰, 等. 基于數(shù)據(jù)挖掘、網(wǎng)絡(luò)藥理學(xué)及臨床觀察探究賈躍進(jìn)治療女性失眠癥用藥規(guī)律及機制 [J]. 中草藥, 2021, 52(15): 4626-4640.

[19] Johnson V L, Hunter D J. The epidemiology of osteoarthritis [J]., 2014, 28(1): 5-15.

[20] Mandl L A. Osteoarthritis year in review 2018: Clinical [J]., 2019, 27(3): 359-364.

[21] Sacitharan P K, Vincent T L. Cellular ageing mechanisms in osteoarthritis [J]., 2016, 27(7/8): 421-429.

[22] Loeser R F. Aging and osteoarthritis [J]., 2011, 23(5): 492-496.

[23] Vina E R, Kwoh C K. Epidemiology of osteoarthritis: Literature update [J]., 2018, 30(2): 160-167.

[24] Gokhale J A, Frenkel S R, Dicesare P E. Estrogen and osteoarthritis [J]., 2004, 33(2): 71-80.

[25] Roman-Blas J A, Casta?eda S, Largo R,. Osteoarthritis associated with estrogen deficiency [J]., 2009, 11(5): 241.

[26] Ching K, Houard X, Berenbaum F,. Hypertension meets osteoarthritis-revisiting the vascular aetiology hypothesis [J]., 2021, 17(9): 533-549.

[27] Ran L, Chen Q, Zhang J Y,. The multimorbidity of hypertension and osteoarthritis and relation with sleep quality and hyperlipemia/hyperglycemia in China’s rural population [J]., 2021, 11(1): 17046.

[28] Courties A, Berenbaum F, Sellam J. The phenotypic approach to osteoarthritis: A look at metabolic syndrome-associated osteoarthritis [J]., 2019, 86(6): 725-730.

[29] 田代華整理. 黃帝內(nèi)經(jīng)素問 [M]. 北京: 人民衛(wèi)生出版社, 2005: 58-87.

[30] 明·張介賓著. 孫玉信, 朱平生主校. 景岳全書 [M]. 上海: 第二軍醫(yī)大學(xué)出版社, 2006: 331.

[31] 清·葉天士著. 蘇禮整理. 臨證指南醫(yī)案 [M]. 北京: 人民衛(wèi)生出版社, 2006: 340.

[32] 婁多峰, 婁玉鈐, 李滿意. 中華痹病大全 [M]. 北京: 中國醫(yī)藥科技出版社, 2019: 245.

[33] 劉建勛. 中華醫(yī)學(xué)百科全書[M]. 北京: 中國協(xié)和醫(yī)科大學(xué)出版社, 2020: 52.

[34] 曹彭凱, 王曉猛, 白偉俠, 等. 2018年版《骨關(guān)節(jié)炎診療指南》解讀 [J]. 河北醫(yī)科大學(xué)學(xué)報, 2018, 39(11): 1241-1243.

[35] Bastiaansen-Jenniskens Y M, Siawash M, van de Lest C H A,. Monounsaturated and saturated, but not n-6 polyunsaturated fatty acids decrease cartilage destruction under inflammatory conditions: A preliminary study [J]., 2013, 4(4): 321-328.

[36] Sekar S, Panchal S K, Ghattamaneni N K,. Dietary saturated fatty acids modulate pain behaviour in trauma-induced osteoarthritis in rats [J]., 2020, 12(2): E509.

[37] Li Y S, Xiao W F, Luo W,. Alterations of amino acid metabolism in osteoarthritis: Its implications for nutrition and health [J]., 2016, 48(4): 907-914.

[38] Qu Z H, Yang F K, Hong J Q,. Causal relationship of serum nutritional factors with osteoarthritis: A mendelian randomization study [J].(Oxford), 2021, 60(5): 2383-2390.

[39] Menarim B C, El-Sheikh Ali H, Loux S C,. Transcriptional and histochemical signatures of bone marrow mononuclear cell-mediated resolution of synovitis [J]., 2021, 12: 734322.

[40] Steinbeck M J, Nesti L J, Sharkey P F,. Myeloperoxidase and chlorinated peptides in osteoarthritis: Potential biomarkers of the disease [J]., 2007, 25(9): 1128-1135.

[41] Deng Z Q, Hu X T, Alahdal M,. High expression of MAPK-14 promoting the death of chondrocytes is an important signal of osteoarthritis process [J]., 2021, 9: e10656.

[42] Xu X, Li X F, Liang Y J,. Estrogen modulates cartilage and subchondral bone remodeling in an ovariectomized rat model of postmenopausal osteoarthritis [J]., 2019, 25: 3146-3153.

[43] Luo S M, Li W J, Wu W R,. Elevated expression of MMP8 and MMP9 contributes to diabetic osteoarthritis progression in a rat model [J]., 2021, 16(1): 64.

[44] Li Q, Wu M J, Fang G F,. microRNA?186?5p downregulation inhibits osteoarthritis development by targeting MAPK1 [J]., 2021, 23(4): 253.

[45] Kojima F, Naraba H, Miyamoto S,. Membrane-associated prostaglandin E synthase-1 is upregulated by proinflammatory cytokines in chondrocytes from patients with osteoarthritis [J]., 2004, 6(4): R355-R365.

[46] Qing Z, Ye J M. Association between ACE polymorphisms and osteoarthritis susceptibility [J]., 2015, 8(6): 7391-7396.

[47] Suri S, Walsh D A. Osteochondral alterations in osteoarthritis [J]., 2012, 51(2): 204-211.

[48] Slobodnick A, Krasnokutsky S, Lehmann R A,. Colorectal cancer among gout patients undergoing colonoscopy [J]., 2019, 25(8): 335-340.

[49] Perrot S. Osteoarthritis pain [J]., 2015, 29(1): 90-97.

Analysis on population characteristics, combined medication, curative effect and mechanism prediction of Cervus and Cucumis Polypeptide Injection in treatment of osteoarthritis based on real world and network of targets

LIU Yi1, XIE Yan-ming1, LI Yuan-yuan1, GAO Ming-hui2

1. Institute of Basic Research in Clinical Medicine, China Academy of Chinese Medical Sciences, Beijing 100700, China 2. School of Statistics, Renmin University of China, Beijing 100872, China

To explore the characteristics of the population and combination of drugs in the real world for the treatment of osteoarthritis with Cervus and Cucumis Polypeptide Injection (鹿瓜多肽注射液) and to predict its mechanism by network analysis.The information of 4050 patients with osteoarthritis treated with Cervus and Cucumis Polypeptide injection from the HIS database was retrieved, and SAS9.2 and R software were used to conduct frequency and rate statistics on their population characteristics and medication; By analyzing the outcomes of patients with osteoarthritis who used Cervus and Cucumis Polypeptide Injection to judge its effectiveness, and compared with the treatment outcomes of osteoarthritis patients who did not use the Cervus and Cucumis Polypeptide Injection; Clementine 12.0 was used to analyze the combined drugs were analyzed by association rules and sorted by confidence. Utilized CNKI, TCMSP database and Cytoscape software to construct “disease-drug-ingredient-target” network and protein-protein interaction network, and performed topological analysis to screen core ingredients and core targets; Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis was performed; Finally, performed molecular docking between the core ingredients and the core targets.Among the 4050 patients with osteoarthritis, the age group of 46—65 years was the most (2248 cases, 55.51%), followed by 66—75 years old (850 cases, 20.99%), among which female patients accounted for the majority(2657, 70.35%). The total effective rate of 4050 osteoarthritis patients treated with Cervus and Cucumis Polypeptide Injection was 97.88%, and the comparison of osteoarthritis patients with and without Cervus and Cucumis Polypeptide Injection showed that the efficacy of Cervus and Cucumis Polypeptide Injection was more effective. Association rules showed that the common concomitant drugs of Cervus and Cucumis Polypeptide Injection in the treatment of osteoarthritis are hormones (betamethasone), antibiotics (cefazolin), anti-inflammatory and analgesic (bone peptide injection), and traditional Chinese medicine to tonify kidney and strengthen bones (Jintiange Capsules), promote blood circulation and remove stasis (Maixuekang Capsule). Target-network analysis obtained 10 core ingredients including palmitic acid, myristic acid, tryptophan, vitamin A, etc; Nine core targets including: peroxisome proliferator-activated receptor gamma(PPARG), myeloperoxidase (MPO), receptor-type tyrosine-protein phosphatase C (PTPRC), etc, among which the combination with better molecular docking were: vitamin A and mitogen-activated protein kinase 14 (MAPK14), tryptophan and matrix metalloproteinase 9 (MMP9), vitamin A and PTPRC. The results of GO and KEGG analysis mainly involved: immune inflammation, hormone secretion, blood circulation and MAPK pathway.The patients with osteoarthritis treated by Cervus and Cucumis Polypeptide Injection are mainly female middle-aged and elderly people, and the common syndromes are liver and kidney deficiency, damp-heat, qi stagnation and blood stasis syndrome; Chinese medicine treatment is mostly combined with kidney-tonifying and bone-strengthening, blood circulation-promoting and stasis-removing drugs to treat both the incidental and fundamental aspects , and Western medicine treatment often concentrate in relieving symptoms. However, through target-network research, it was concluded that Cervus and Cucumis Polypeptide Injection mainly play a role in the regulation of,and other targets through vitamin A and amino acids. This conclusion can provide reference for further research and clinical application.

real world; target network; Cervus and Cucumis Polypeptide Injection; population characteristics; combined medication; osteoarthritis; hospital information system

R285.64

A

0253 - 2670(2022)17 - 5445 - 14

10.7501/j.issn.0253-2670.2022.17.021

2022-03-20

國家重點研發(fā)計劃項目（2018YFC1707400）；2018重點研發(fā)計劃：中醫(yī)“治未病”服務(wù)技術(shù)體系的構(gòu)建（2018YFC1704705）；2021年岐黃學(xué)者支持項目

劉 毅（1993—），男，在讀博士研究生，研究方向為中藥上市后再評價方法學(xué)研究。E-mail: 709473245@qq.com

謝雁鳴（1959—），女，教授，研究員，博士生導(dǎo)師，中國中醫(yī)科學(xué)院首席研究員，研究方向為中藥上市后再評價方法學(xué)研究、中西醫(yī)結(jié)合臨床、老年病學(xué)。E-mail: ktzu2018@163.com

黎元元，博士，副研究員，研究方向為中藥上市后再評價。Tel: (010)64014411-3351 E-mail: chibjyy@163.com

[責(zé)任編輯 潘明佳]