• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于網(wǎng)絡(luò)穩(wěn)健性的和血明目片潛在適應(yīng)證預(yù)測

    2022-06-02 13:11:58席玉婕高錦成郭非非楊洪軍
    世界中醫(yī)藥 2022年7期

    席玉婕 高錦成 郭非非 楊洪軍

    摘要 目的:基于藥物干擾疾病網(wǎng)絡(luò)穩(wěn)健性的方法來預(yù)測和血明目片潛在的適應(yīng)證。方法:使用中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫與分析平臺(TCMSP)、中醫(yī)藥綜合數(shù)據(jù)庫(TCMID)等數(shù)據(jù)庫搜集和血明目片中的活性成分,利用BATMAN-TCM預(yù)測成分的潛在靶點;通過DisGeNET和HPO數(shù)據(jù)庫獲取相應(yīng)的疾病蛋白靶標(biāo),基于靶點移除前后整體網(wǎng)絡(luò)拓?fù)涮卣鞯淖兓u價藥物靶點對于網(wǎng)絡(luò)穩(wěn)定性的影響,實現(xiàn)和血明目片干預(yù)不同疾病藥效的預(yù)測。通過建立“中藥成分-藥物靶點-疾病蛋白”關(guān)聯(lián)網(wǎng)絡(luò),從中篩選出潛在藥物作用靶點并對其進行富集分析,考察和血明目片干預(yù)疾病可能的生物學(xué)過程和信號通路,最終從體外細(xì)胞實驗的角度初步對藥效進行驗證。結(jié)果:搜集整理后得到和血明目片所含化學(xué)成分946個,預(yù)測高可信度靶點761個。通過考察藥物擾動疾病網(wǎng)絡(luò)穩(wěn)健性,發(fā)現(xiàn)和血明目片對翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化的疾病網(wǎng)絡(luò)有較強的擾動作用。通過構(gòu)建和血明目片干預(yù)不同疾病的“中藥成分-藥物靶點-疾病蛋白”異質(zhì)網(wǎng)絡(luò),獲得相應(yīng)的潛在藥物作用靶點,多項富集分析結(jié)果顯示和血明目片可能通過干預(yù)血管新生、白細(xì)胞遷移、脂多糖(LPS)應(yīng)答、核因子κB轉(zhuǎn)錄調(diào)控、凋亡通路、氧化應(yīng)激通路等過程來發(fā)揮作用。以上4種疾病潛在作用通路與和血明目片的潛在作用通路具有一致性,且體外細(xì)胞實驗很好地證明了和血明目片對于治療視網(wǎng)膜損傷具有一定的藥效。結(jié)論:和血明目片具有“多組成、多靶標(biāo)、多通路”的特點,對翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化等4種適應(yīng)證存在潛在的治療作用。

    關(guān)鍵詞 和血明目片;網(wǎng)絡(luò)穩(wěn)健性;富集分析;潛在適應(yīng)證

    Prediction of Potential Indications of Hexuemingmu Tablets Based on Network Robustness

    XI Yujie1,2,GAO Jincheng3,GUO Feifei1,YANG Hongjun2,4

    (1 Institute of Chinese Materia Medica,China Academy of Chinese Medical Sciences,Beijing 100700,China; 2 Tianjin University of Traditional Chinese Medicine,Tianjin 301617,China; 3 Shaanxi Provincial Hospital of Chinese Medicine,Xi′an 710003,China; 4 Experimental Research Center,China Academy of Chinese Medical Sciences,Beijing 100700,China)

    Abstract Objective:To predict the potential indications of Hexuemingmu tablets(HXMM tablets) based on the method of medicinal interfering on the network robustness of disease.Methods:The active components of HXMM tablets were collected by Traditional Chinese Medicine Integrated Database(TCMID) and Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform(TCMSP),and the potential targets of components were predicted by Bioinformatics Analysis Tool for Molecular Mechanism of Traditional Chinese Medicine(BATMAN-TCM).The corresponding protein targets of diseases were obtained from DisGeNET and Human Phenotype Ontology(HPO) databases.The effects of medicine targets on the network robustness were evaluated based on the changes of topological characteristics of the whole network after targets were removed,and the pharmacological effect of HXMM tablets interfering on different diseases were predicted.The correlation network of “Chinese medicine chemical component-medicine target-disease protein” was established,and the potential medicine action targets were screened out to conduct enrichment analysis.The possible biological processes and signal pathways of HXMM tablets interfering on diseases were investigated,and the pharmacological effect was eventually verified by in vitro cell experiments.Results:A total of 946 chemical components of HXMM tablets were obtained and 761 highly credible targets were predicted after collection.By investigating the medicinal interference on the network robustness of disease,this study found that HXMM tablets had a strong interfering effect on the disease networks of the four indications(pterygium,keratitis,ocular hypertension,and corneal neovascularization).The potential medicine active targets were obtained from the heterogeneous network of “Chinese medicine chemical component-medicine target-disease protein” of different diseases intervened by HXMM tablets.The results of multiple enrichment analyses showed that HXMM tablets may play a role by interfering on angiogenesis,leukocyte traffic,lipopolysaccharide(LPS) response,nuclear factor kappa B(NF-κB) transcriptional regulation,apoptosis pathway,oxidative stress pathway,and so on.The potential action pathways of the four indications and HXMM tablets were consistent.Conclusion:HXMM tablet has potential pharmacological effects on treating the four indications(pterygium,keratitis,ocular hypertension,and corneal neovascularization) with “multi-composition,multi-target,and multi-pathway” characteristics.

    Keywords Hexuemingmu tablets; Network robustness; Enrichment analysis; Potential indications

    中圖分類號:R285;R285.5文獻標(biāo)識碼:Adoi:10.3969/j.issn.1673-7202.2022.07.004

    和血明目片由蒲黃、丹參、地黃、墨旱蓮、菊花、黃芩(炭)、決明子、車前子、茺蔚子、女貞子、夏枯草、龍膽、郁金、木賊、赤芍、牡丹皮、山楂、當(dāng)歸、川芎共計19味藥組成。全方根據(jù)四物湯和二至丸加減[1-2],其中君藥生地黃、當(dāng)歸、赤芍、川芎補散結(jié)合,發(fā)揮活血化瘀、補血之功,補而不滯,行而不破,藥專力宏,墨旱蓮、女貞子以滋補肝腎,生精止血為主,為平補虛之良劑;臣藥中蒲黃、牡丹皮、丹參、茺蔚子既涼血止血又化瘀利水,和血效佳,山楂、郁金偏于行氣導(dǎo)滯,可消積散瘀;黃芩、菊花、龍膽、夏枯草、木賊為佐藥,上能疏散風(fēng)熱清利頭目,下可清熱燥濕瀉肝膽火;使藥決明子、車前子協(xié)同發(fā)揮清肝明目、利水通便之效。因此和血明目片具有涼血止血、滋陰化瘀、養(yǎng)肝明目的功效,用于陰虛肝旺,熱傷絡(luò)脈所引起的眼底出血[3]。

    和血明目片目前收錄在2020版《中華人民共和國藥典》中,在年齡相關(guān)性黃斑變性、糖尿病視網(wǎng)膜病變、視網(wǎng)膜靜脈阻塞等疾病中得到廣泛應(yīng)用[4-6]?,F(xiàn)代藥理學(xué)研究表明該藥在早期止血,中后期促進瘀血吸收,改善微循環(huán),擴張血管、抗血栓形成等方面療效較好[7-8]。基于本課題組開發(fā)的藥效預(yù)測平臺調(diào)研翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化、角膜不透明、錐形桿狀肌營養(yǎng)不良癥、結(jié)膜炎、瞼板腺炎、干燥性角膜結(jié)膜炎、色素性視網(wǎng)膜炎等相關(guān)數(shù)據(jù)。結(jié)果發(fā)現(xiàn)和血明目片干預(yù)翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化的擾動評分大于零,提示和血明目片可能對于以上4種適應(yīng)證具有潛在治療作用。

    1 資料與方法

    1.1 和血明目潛在藥物靶點預(yù)測 利用中藥系統(tǒng)藥理學(xué)數(shù)據(jù)庫與分析平臺(Traditional Chinese Medicine Systems Pharmacology,TCMSP)和中醫(yī)藥綜合數(shù)據(jù)庫(Traditional Chinese Medicine Integrated Database,TCMID)等數(shù)據(jù)庫獲取和血明目片中蒲黃、丹參、地黃、墨旱蓮、菊花、黃芩(炭)、決明子、車前子、茺蔚子、女貞子、夏枯草、龍膽、郁金、木賊、赤芍、牡丹皮、山楂、當(dāng)歸、川芎等中藥的化學(xué)成分信息[9-10],通過BATMAN-TCM數(shù)據(jù)庫[11](http://bionet.ncpsb.org/batman-tcm/)獲得藥物成分可能發(fā)生相互作用的靶標(biāo)蛋白,并選擇其中高可信度(打分值大于30)的蛋白作為藥物的潛在靶點。

    1.2 疾病數(shù)據(jù)的搜集 翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化、角膜不透明、錐形桿狀肌營養(yǎng)不良癥、結(jié)膜炎、瞼板腺炎、干燥性角膜結(jié)膜炎、色素性視網(wǎng)膜炎等10種疾病靶標(biāo)均來源于Human Phenotype Ontology數(shù)據(jù)庫(http://humanphenotype-ontology.github.io/)和DisGeNET(http://www.disgenet.org)數(shù)據(jù)庫[12-13]。

    1.3 靶標(biāo)的功能富集分析 藥物、疾病的靶標(biāo)功能分析中涉及的基因本體(Gene Ontology,GO)功能條目的富集分析是由Metascape軟件來實現(xiàn)的[14],具體的原理是基于超幾何分布設(shè)計的算法,而多重檢驗校正則使用Benjamini-Hochberg校正方法。

    1.4 基于網(wǎng)絡(luò)穩(wěn)健性的藥物干預(yù)疾病能力預(yù)測

    1.4.1 眼科疾病網(wǎng)絡(luò)的模 疾病網(wǎng)絡(luò)由疾病相關(guān)的基因之間的相互作用組成。采用來自于STRING[15](version 10)數(shù)據(jù)庫的置信度大于0.4的蛋白質(zhì)-蛋白質(zhì)相互作用(Protein-protein Interaction,PPI)網(wǎng)絡(luò)構(gòu)建(下載于2020年1月26日)。利用Cytoscape(version 3.4.0)對網(wǎng)絡(luò)進行可視化[16],計算藥物干預(yù)前后反映網(wǎng)絡(luò)穩(wěn)健性的拓?fù)涮卣?,包括整體網(wǎng)絡(luò)的平均最短路徑、平均連接度、連接度中心性和緊密度中心性[17]。

    1.4.2 藥物干預(yù)網(wǎng)絡(luò)穩(wěn)健性的評價方法 藥物干預(yù)前后疾病網(wǎng)絡(luò)拓?fù)涮卣鞯淖兓痆18],即藥物對疾病的擾動作用,可以通過在疾病網(wǎng)絡(luò)中刪除藥物靶點來表示。經(jīng)過R的iGraph包[19]進行計算后,平均最短路徑、網(wǎng)絡(luò)平均連接度、網(wǎng)絡(luò)連接度中心性和網(wǎng)絡(luò)緊密度中心性等4個指標(biāo)可用來作為評價干預(yù)前后疾病網(wǎng)絡(luò)的穩(wěn)健性。其中,網(wǎng)絡(luò)的平均連接度是網(wǎng)絡(luò)中所有節(jié)點連接度的平均值;網(wǎng)絡(luò)的連接度中心性為網(wǎng)絡(luò)中的所有節(jié)點的連接度計算而來;網(wǎng)絡(luò)的緊密度中心性也是由網(wǎng)絡(luò)中所有節(jié)點的緊密度計算得出,以上3個數(shù)值大小均與網(wǎng)絡(luò)穩(wěn)健性正相關(guān),數(shù)值越大網(wǎng)絡(luò)越穩(wěn)定,藥物對疾病網(wǎng)絡(luò)的干預(yù)作用相對越弱。而平均最短路徑代表網(wǎng)絡(luò)當(dāng)中所有最短路徑長度的平均值,數(shù)值越小網(wǎng)絡(luò)越穩(wěn)定,正擾動率說明藥物對網(wǎng)絡(luò)穩(wěn)定性破壞能力強。根據(jù)以上指標(biāo)相互作用關(guān)系,擾動率計算公式可表示為:校正前擾動率=(藥物干預(yù)后拓?fù)涮卣髦?藥物干預(yù)前拓?fù)涮卣髦担?藥物干預(yù)前拓?fù)涮卣髦?。而該擾動率通過零分布(Null distribution)校正后才可用于不同藥物的擾動率比較,校正方法見1.4.3部分。藥物擾動總分計算公式如下,公式中的擾動率均利用Null distribution校正:擾動總分=校正后的平均最短路徑擾動率-校正后平均連接度擾動率-校正后連接度中心性擾動率-校正后緊密度中心性擾動率。

    1.4.3 利用藥物干預(yù)隨機網(wǎng)絡(luò)的Null distribution作為總體分布校正藥物對于真實網(wǎng)絡(luò)的擾動率 由于藥物作用靶點大多屬于總體分布未知的小樣本數(shù)據(jù),因此在進行顯著性考察時常使用置換檢驗(Permutation test)。疾病網(wǎng)絡(luò)零分布(Null distribution)的實質(zhì)是在保留真實網(wǎng)絡(luò)節(jié)點數(shù)和邊數(shù)后生成隨機的網(wǎng)絡(luò)的拓?fù)涮卣髦怠kS機網(wǎng)絡(luò)構(gòu)成的Null distribution可作為基準(zhǔn),與真實疾病網(wǎng)絡(luò)的藥物作用進行比較,考察在相同的藥物擾動后,其擾動率的總體分布以及真實網(wǎng)絡(luò)和隨機網(wǎng)絡(luò)之間的網(wǎng)絡(luò)穩(wěn)健性是否存在顯著差異。通過與100個隨機網(wǎng)絡(luò)產(chǎn)生的Null distribution比較,真實網(wǎng)絡(luò)的擾動率能比較客觀地評價藥物在實際疾病網(wǎng)絡(luò)中引起顯著影響的干擾程度。此外,利用藥物在真實網(wǎng)絡(luò)擾動率在Null distribution中對應(yīng)的Z-score可用于評價校正后的網(wǎng)絡(luò)穩(wěn)健性的擾動率。Z-score越大,表明藥物干預(yù)后網(wǎng)絡(luò)穩(wěn)健性越差,藥物對疾病網(wǎng)絡(luò)的影響越大。

    1.5 “中藥成分-藥物靶點-疾病蛋白”異質(zhì)網(wǎng)絡(luò)的構(gòu)建 基于網(wǎng)絡(luò)穩(wěn)健性原理,將前期預(yù)測到的潛在藥物靶點及疾病數(shù)據(jù)導(dǎo)入中藥藥效預(yù)測平臺開展相關(guān)研究[10],得到相應(yīng)的擾動總分及“中藥成分-藥物靶點-疾病蛋白”的異質(zhì)網(wǎng)絡(luò),該網(wǎng)絡(luò)呈現(xiàn)了中藥成分與疾病之間以及疾病蛋白與藥物靶點之間的關(guān)聯(lián)關(guān)系,為中藥干預(yù)疾病的機制探索提供線索。該網(wǎng)絡(luò)包含了3種類型的節(jié)點,不同類型的結(jié)果分別使用不同的顏色和形狀來表示。其中粉色三角形表示中藥組成成分,藍色菱形表示藥物的潛在靶標(biāo),紫色圓形表示與適應(yīng)證相關(guān)的疾病蛋白。網(wǎng)絡(luò)還包含了2種聯(lián)系(邊),分別為化合物與靶標(biāo)的聯(lián)系(通過BATMAN-TCM來預(yù)測,打分大于30)、藥物靶點與疾病蛋白的PPI關(guān)系(來源于STRING數(shù)據(jù)庫中打分大于0.4的PPI)。為了考察和血明目片干預(yù)潛在適應(yīng)證可能參與的生物學(xué)過程,通過計算紅色矩形中藥物作用的潛在靶標(biāo)(疾病蛋白與藥物靶點的交集)的網(wǎng)絡(luò)拓?fù)涮卣髦担Y選連接度(Degree)大于中位數(shù)的值進行GO生物學(xué)富集分析。

    1.6 和血明目片干預(yù)潛在適應(yīng)證的作用特點探索? 根據(jù)上述和血明目片干預(yù)潛在適應(yīng)證的疾病蛋白富集分析結(jié)果,整理顯著性排名前5位的生物學(xué)過程。進行和血明目片作用于適應(yīng)證的環(huán)節(jié)探索?;谇捌陬A(yù)測的潛在適應(yīng)證藥物的潛在作用靶標(biāo),通過計算其網(wǎng)絡(luò)拓?fù)涮卣髦担Y選大于二倍Degree的靶標(biāo),認(rèn)為可能是疾病核心作用靶標(biāo)。將所有潛在適應(yīng)證的疾病核心作用靶標(biāo)取交集,并進行作用特點探索。

    1.7 分子生物學(xué)實驗驗證

    1.7.1 細(xì)胞 人視網(wǎng)膜上皮細(xì)胞ARPE-19(ATCC,美國,批號:20210603),培養(yǎng)條件為1640培養(yǎng)基添加10%胎牛血清(Fetal Bovine Serum,F(xiàn)BS)和1%青霉素-鏈霉素(Penicillin-Streptomycin,PS),二氧化碳培養(yǎng)箱內(nèi)培養(yǎng)(溫度:37 ℃,濕度:100%,5% CO2和95% O2)。

    1.7.2 藥物 和血明目片藥粉(西安碑林藥業(yè)股份有限公司,生產(chǎn)批號:20210428)。

    1.7.3 試劑與儀器 RPMI-1640培養(yǎng)基(貨號:11875093)、胎牛血清(貨號:1600044)、青-鏈霉素混合液(貨號:15070063)、胰蛋白酶-乙二胺四乙酸消化液(貨號:25200072)、磷酸鹽緩沖液(貨號:10010023)以上均購美國GIBCO公司;細(xì)胞凍存液(蘇州新賽美生物科技有限公司,貨號:C40100),CCK-8檢測試劑盒(Dojindo公司,日本,貨號:CK04),貝伐單抗(上海阿拉丁生化科技股份有限公司,貨號:B303060),氯化鈉(NaCl,北京化工廠,批號:20160223),氯化鉀(KCl,批號:20130612)、碳酸氫鈉(NaHCO3,批號:20150326)、磷酸二氫鈉(NaHPO3,批號:20090727)、氯化鎂(MgCl2,批號:F20080103)、氯化鈣(CaCl2,批號:F20090602)、葡萄糖(C6H12O6,批號:20150303)、以上均購自國藥集團化學(xué)試劑北京有限公司。CO2培養(yǎng)箱(SANYO公司,日本,型號:MCO-18AIC),超凈工作臺(蘇州凈化設(shè)備廠,型號:SW-CF-2FD),倒置顯微鏡(OLYMPUS公司,日本,型號:CKX41),酶聯(lián)免疫檢測儀(Thermo Fisher Scientific公司,美國,型號:MμltiskanTM FC),-86 ℃超低溫冰箱(青島海爾特種電器有限公司,型號:DW-86L728),離體腸灌流系統(tǒng)(上海奧爾科特生物科技有限公司,型號:ALC-M),電子天平(METTLER TOLEDO公司,瑞士,型號:ME54),臺式循環(huán)水式多用真空泵(鄭州長城科工貿(mào)有限公司,型號:SHB-Ⅲ),96孔板(型號:3599)、T25細(xì)胞培養(yǎng)瓶(型號:430639)以上購于美國Corning公司。

    1.7.4 含藥腸吸收液的制備

    1.7.4.1 臺式緩沖液的配制 NaCl 8.0 g,NaHCO3 1.0 g,KCl 0.28 g,MgCl20.1 g,NaH2PO4 0.05 g溶于500 mL蒸餾水中,密閉冷藏;CaCl2 0.2 g溶于500 mL蒸餾水中密閉冷藏;使用將前二者均勻攪拌,最后加入葡萄糖1.0 g混合即可。

    1.7.4.2 含藥腸吸收液的制備 將實驗前禁食12 h的大鼠稱重后,斷頸椎處死,沿腹中線和腹白線分別剪開皮膚與肌肉。迅速取出小腸,自胃幽門以下10 cm開始向下為空腸段,自胃幽門以下10 cm開始向下量取14 cm為腸段1;間隔10 cm量取14 cm為腸段2;再間隔10 cm量取14 cm為腸段3;從回盲瓣以上5 cm開始向上取14 cm為腸段4。將剪下的腸管放入0 ℃臺氏液中沖洗至無腸內(nèi)容物為止。將自制硅膠套管軟端插入腸管用絲線結(jié)扎,小心將腸道翻轉(zhuǎn)。用0 ℃臺氏液沖洗內(nèi)表面,將另一端用絲線結(jié)扎成囊狀。麥?zhǔn)显」苤蓄A(yù)先加入含藥供試藥液25 mL,通入混合氣體(95% O2,5% CO2),保持37 ℃恒溫。在腸管中注入臺氏液2 mL,放入麥?zhǔn)显」苤校?20 min后收取樣品。樣品放入干凈的離心管中,于-80 ℃密封直立保存,待用[20]。

    1.7.5 方法

    1.7.5.1 分組與模型制備 取人視網(wǎng)膜上皮細(xì)胞正常復(fù)蘇,使用含10%胎牛血清的1640培養(yǎng)基,于37 ℃,5% CO2,95%O2條件下進行正常培養(yǎng)。將細(xì)胞分為正常對照組、模型組、貝伐單抗-陽性藥組(1 mg/mL)、和血明目腸吸收液低濃度組(600 μg/mL)、中濃度組(300 μg/mL)、低濃度組(150 μg/mL)。待ARPE-19細(xì)胞密度達到70%~80%時,傾倒原培養(yǎng)液,改用無血清培養(yǎng)基饑餓處理12 h,不同濃度過氧化氫(80、100、150、200、300 μmol/L)處理2 h后終止培養(yǎng),每孔小心加入10 μL CCK-8溶液,盡量避免氣泡的產(chǎn)生,放入培養(yǎng)箱內(nèi)避光孵育1~4 h,酶標(biāo)儀于450 nm下測量孔板的吸光度(A)值,計算細(xì)胞存活率,公式為細(xì)胞成活率=(A觀察組-A空白組)/(A對照組-A空白組)×100%。

    1.7.5.2 給藥方法 常規(guī)培養(yǎng)ARPE-19細(xì)胞,待ARPE-19細(xì)胞密度達到對數(shù)成長期時配制成混懸液,接種在96孔板,接種密度為2×104個細(xì)胞/孔。過夜培養(yǎng)細(xì)胞后進行饑餓處理12 h,正常組與模型組加入無血清培養(yǎng)基,給藥組加入高中低不同濃度的和血明目腸吸收液,作用細(xì)胞24 h,隨后造模2 h,最后使用CCK-8法測定其細(xì)胞成活率。

    1.7.6 統(tǒng)計學(xué)方法 采用GraphPad Prism 8.3.0統(tǒng)計軟件進行數(shù)據(jù)分析,計量資料均以均數(shù)±標(biāo)準(zhǔn)差(±s)表示,多組間均數(shù)比較采用Kruskal-Wallis H檢驗,以P<0.05為差異有統(tǒng)計學(xué)意義。

    2 結(jié)果

    2.1 和血明目片潛在藥物靶點及潛在機制分析?? 基于上述方法得到和血明目片所含化學(xué)成分226個,高可信度靶標(biāo)854個。富集分析結(jié)果表明和血明目片的潛在靶標(biāo)富集的GO生物學(xué)功能條目包括血液循環(huán)(Blood Circulation)、脂質(zhì)生物合成(Lipid Biosynthetic Process)、氧化應(yīng)激(Response to Oxidative Stress,Reactive Oxygen Species Metabolic Process)、血管收縮(Positive Regulation of Vasoconstriction)等。見圖1。

    圖1 藥物潛在作用靶標(biāo)GO功能富集分析氣泡圖? 注:橫坐標(biāo)代表Enrichment,縱坐標(biāo)代表-Log(q-value),節(jié)點大小與之匹配的靶標(biāo)的數(shù)目有關(guān)

    2.2 和血明目片干預(yù)疾病網(wǎng)絡(luò)穩(wěn)健性分析 按照擾動總分對搜集獲得的9類適應(yīng)證由大至小進行排序。結(jié)果顯示,和血明目片干預(yù)9個潛在適應(yīng)證后,擾動總分大于0的共有4個,分別為翼狀胬肉、角膜炎、高眼壓癥和角膜新生血管化,擾動總分小于0的共有5個。見表2。

    2.3 和血明目片干預(yù)翼狀胬肉、角膜炎、高眼壓癥、角膜血管新生化的機制預(yù)測 和血明目片干預(yù)翼狀胬肉的“中藥成分-藥物靶點-疾病蛋白”異質(zhì)網(wǎng)絡(luò),包含了221個和血明目片的組成化合物、64個和血明目片的潛在靶標(biāo)以及122個與藥物關(guān)聯(lián)的疾病蛋白。其中,節(jié)點的Degree值排名前列的關(guān)鍵靶點有TP53、AKT1、EGFR、EGF、FN1等。見圖2A。和血明目片干預(yù)角膜炎的“中藥成分-藥物靶標(biāo)-疾病蛋白”異質(zhì)網(wǎng)絡(luò),包含了216個和血明目片的組成化合物、41個和血明目片的潛在靶標(biāo)以及95個與藥物關(guān)聯(lián)的疾病蛋白。其中,節(jié)點的Degree值排名前列的關(guān)鍵靶點有IL-10、IL-6、APOE及CAT等。見圖2B。和血明目片干預(yù)高眼壓癥的“中藥成分-藥物靶標(biāo)-疾病蛋白”異質(zhì)網(wǎng)絡(luò),包含了236個和血明目片的組成化合物、33個和血明目片的潛在靶標(biāo)以及47個與藥物關(guān)聯(lián)的疾病蛋白。其中,節(jié)點的Degree值排名前列的關(guān)鍵靶點有NR3C1、FAS、JUN、SHH及IL-1β等。見圖2C。和血明目片干預(yù)角膜新生血管化的“中藥成分-藥物靶標(biāo)-疾病蛋白”異質(zhì)網(wǎng)絡(luò),包含了189個和血明目片的組成化合物、32個和血明目片的潛在靶標(biāo)以及59個與藥物關(guān)聯(lián)的疾病蛋白。其中,節(jié)點的Degree值排名前列的關(guān)鍵靶點有CTNNB1、FGF2、IL-6、APOE及AKT1等。見圖2D。

    和血明目片干預(yù)翼狀胬肉、角膜炎、高眼壓癥、角膜血管新生化等適應(yīng)證中與藥物作用相關(guān)的疾病靶標(biāo)功能富集分析。其中,和血明目片干預(yù)翼狀胬肉顯著富集排名前5的生物學(xué)過程分別為血管新生(Blood Vessel Development,Angiogenesis)、細(xì)胞外基質(zhì)(Extracellular Matrix)、氧化損傷(Response to Wounding,Response to Oxygen Levels)、細(xì)胞增殖(Muscle Cell Proliferation)及轉(zhuǎn)移酶活性調(diào)節(jié)(Positive Regulation of Transferase activity)等。見圖3A。和血明目片干預(yù)角膜炎中位居前列顯著富集的生物學(xué)過程分別為炎癥反應(yīng)(Inflammation Response)、脂多糖反應(yīng)(Response to Lipopolysaccharide)、白細(xì)胞遷移(Leukocyte Migration)、核因子κB轉(zhuǎn)錄因子活性調(diào)節(jié)(Positive Regulation of NF-kappa B Transcription Factor Activity)及MAPK級聯(lián)的正向調(diào)節(jié)(Positive Regulation of MAPK Cascade)等。見圖3B。和血明目片干預(yù)高眼壓癥中位居前列顯著富集的生物學(xué)過程分別為凋亡信號通路(Apoptotic Signaling Pathway)、MAPK級聯(lián)的正向調(diào)節(jié)(Positive regulation of MAPK Cascade)、上皮細(xì)胞增殖(Epithlial Cell Proliferation)、血管新生(Blood Vessel Development)及細(xì)胞損傷(Response to Wounding)等。見圖3C。和血明目片干預(yù)角膜新生血管化中位居前列顯著富集的生物學(xué)過程分別為血管新生(Blood Vessel Development,Vasculature Development)、細(xì)胞遷移(Positive Regulation of Cell Migration)、炎癥反應(yīng)(Inflammation Response)、上皮細(xì)胞增殖(Epithlial Cell Proliferation)及MAPK級聯(lián)調(diào)節(jié)(Regulation of MAPK Cascade)等。見圖3D。

    2.4 和血明目干預(yù)4個潛在適應(yīng)證靶點作用特點? 血管新生和細(xì)胞增殖可能是和血明目片作用于翼狀胬肉,高眼壓癥和角膜血管新生化共有的2個環(huán)節(jié);MAPK級聯(lián)調(diào)節(jié)可能是角膜炎,高眼壓癥和角膜血管新生化的共有環(huán)節(jié);炎癥反應(yīng)可能是角膜炎和角膜血管新生化的共有環(huán)節(jié)。見表3。通過篩選疾病核心作用靶標(biāo)的網(wǎng)絡(luò)拓?fù)涮卣髦?,得到治療翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化疾病核心作用靶點分別為64、41、33、32個。見圖4。其中出現(xiàn)次數(shù)最多的共同核心靶點有MMP9、CXCR4、IL6和IL1B。翼狀胬肉與高眼壓癥的共有靶點有TP53、JUN,角膜炎和高眼壓癥共有靶點有核因子κB1,高眼壓癥與角膜血管新生化共有靶點有ALB,角膜炎與角膜血管新生化特異性共有靶標(biāo)有4個,分別是APOE、CCL2、IKBKB、TSPO。

    2.5 實驗驗證

    2.5.1 過氧化氫誘導(dǎo)ARPE-19細(xì)胞損傷 與正常組比較,過氧化氫濃度在150~300 μmol/L范圍內(nèi)的細(xì)胞成活率均有統(tǒng)計學(xué)意義(P<0.05),后續(xù)實驗均選擇接近半數(shù)抑制濃度的200 μmol/L的過氧化氫進行造模。見表4。

    2.5.2 和血明目方對ARPE-19細(xì)胞的保護作用?? 相較于模型組,預(yù)防給藥24 h后,濃度為600 μg/mL的和血明目方高劑量組的細(xì)胞活力較為顯著(P<0.01),300 μg/mL的和血明目中劑量組和1 mg/mL的貝伐單抗陽性藥的細(xì)胞活力差異有統(tǒng)計學(xué)意義(P<0.05)。見圖5。

    3 討論

    作為中華民族的文明瑰寶,中醫(yī)藥以其獨有的自身優(yōu)勢保持著強大的生命力,然而隨著現(xiàn)代科學(xué)飛速發(fā)展,國民健康意識的不斷提高,中藥的活性成分辨識、藥效物質(zhì)基礎(chǔ)研究、作用機制科學(xué)闡釋等復(fù)雜體系解析已成為當(dāng)代重大前沿科學(xué)問題。本研究通過構(gòu)建“中藥成分-藥物靶點-疾病蛋白”的關(guān)聯(lián)網(wǎng)絡(luò),考察和血明目片對于多種眼科疾病網(wǎng)絡(luò)的擾動作用,預(yù)測和血明目片潛在的適應(yīng)證,為中藥復(fù)雜機制闡釋提供了一定線索。

    根據(jù)和血明目片潛在藥物靶點富集分析結(jié)果顯示,其生物學(xué)功能主要與血液循環(huán)、氧化應(yīng)激應(yīng)答、血管收縮等調(diào)控密切相關(guān),這與和血明目片臨床應(yīng)用一致,說明該藥可能通過調(diào)控血液循環(huán)或氧化應(yīng)激等生物過程來干預(yù)眼科相關(guān)疾病?;谥兴幩幮ьA(yù)測平臺分析結(jié)果推測和血明目片對于翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化等4種疾病存在一定干預(yù)作用。表皮生長因子受體在抗血管生成、細(xì)胞增殖及細(xì)胞凋亡方面已取得公認(rèn)的療效[21],提示這些蛋白可能是HXMM治療翼狀胬肉的重要效應(yīng)蛋白。研究表明翼狀胬肉的形成受環(huán)境因素的影響,同時與細(xì)胞凋亡、血管新生、基質(zhì)金屬蛋白酶以及病毒感染有關(guān)[22]。因此和血明目可能通過調(diào)控氧化應(yīng)激等過程來影響細(xì)胞凋亡,同時也通過調(diào)控血管形成來干預(yù)疾病的發(fā)生。和血明目片還可能通過干預(yù)炎癥反應(yīng)白細(xì)胞遷移、脂多糖(LPS)應(yīng)答、核因子κB轉(zhuǎn)錄調(diào)控等功能來發(fā)揮作用。前人研究證實角膜炎的發(fā)生主要與免疫炎癥,細(xì)菌感染相關(guān)結(jié)論一致,由此預(yù)測和血明目可能通過調(diào)控免疫炎癥等過程來干預(yù)角膜炎的發(fā)生[23]。研究表明高眼壓癥的發(fā)病因素與全身系統(tǒng)相關(guān),主要涉及糖尿病、高血壓,以及阻塞性睡眠呼吸暫停綜合征等[24]。因此,和血明目片可能通過調(diào)控以上環(huán)節(jié)來干預(yù)高眼壓癥發(fā)生。在不同的疾病中,導(dǎo)致角膜新生血管形成的主要因素不同,例如角膜移植中的宿主免疫反應(yīng)、角膜緣干細(xì)胞缺乏中干細(xì)胞的維持和分化、接觸鏡相關(guān)眼病中的缺氧、單純皰疹角膜炎中的病毒感染[25]。因此,和血明目可能通過調(diào)控血管新生、炎癥反應(yīng)、氧化應(yīng)激以及凋亡等過程來干預(yù)疾病的發(fā)生。本研究得到4個共同核心靶點基質(zhì)金屬蛋白酶P9、趨化因子CXC亞家族受體4、白細(xì)胞介素6和白細(xì)胞介素-1β,這些靶點通常參與炎癥反應(yīng)、趨化因子表達等環(huán)節(jié),提示和血明目可能主要通過改善炎癥反應(yīng)治療其潛在適應(yīng)證;翼狀胬肉與高眼壓癥的共有靶點有腫瘤蛋白53、Jun,這些靶點的相關(guān)通路常與凋亡相關(guān)[26-27];角膜炎和高眼壓癥共有靶點有核因子κB1,提示這2種疾病的特異性環(huán)節(jié)可能與免疫調(diào)節(jié)介導(dǎo)的核因子κB信號通路密切相關(guān)[28];而白蛋白作為高眼壓癥與角膜血管新生化共有靶標(biāo)與凝血功能異常以及體液免疫密切相關(guān)[29-30];角膜炎與角膜血管新生化特異性共有靶標(biāo)有4個,分別是載脂蛋白E(Apolipoprotein E,ApoE)、C-C基序趨化因子配體2(C-C Motif Chemokine Ligand2,CCL2)、核因子κB激酶亞基B抑制劑(Inhibitor of Nuclear Factor Kappa B Kinase Subunit B,IKKbeta)、轉(zhuǎn)運體蛋白(Translocator Protein,TSPO),主要參與脂質(zhì)代謝、免疫應(yīng)答等[31-32]。而在得到和血明目干預(yù)4種適應(yīng)證的特異性靶標(biāo)后,通過文獻分析發(fā)現(xiàn)。翼狀胬肉代表性環(huán)節(jié)集中在血管新生、凋亡(表皮生長因子受體、FOS、JUN)[33];角膜炎則為炎癥反應(yīng)為主(IL10、胰島素樣生長因子、TAT)[34];干預(yù)高眼壓癥以凝血、氧化應(yīng)激為主要特點(載脂蛋白1、C3、超氧化物歧化酶1、β型轉(zhuǎn)化生長因子1)[35];角膜血管新生化以血管新生為藥效環(huán)節(jié)代表(血管緊張肽原、血纖蛋白溶解酶原、過氧化物酶體增殖物激活受體)[36-37]。以上結(jié)果進一步表明和血明目片可能針對血管新生、炎癥反應(yīng)、氧化應(yīng)激等環(huán)節(jié)進行潛在適應(yīng)證的干預(yù)。

    中醫(yī)學(xué)中認(rèn)為翼狀胬肉多因心肺二經(jīng)風(fēng)熱壅盛導(dǎo)致虛火上炎、氣滯血瘀[38-39],和血明目片可以清熱瀉火、理氣活血為主發(fā)揮藥效;角膜炎大多屬于外感風(fēng)熱或熱毒上攻[40-41],和血明目片可通過祛風(fēng)解表、清肝瀉火進行治療;高眼壓癥常與心火旺盛、肝腎陰虛導(dǎo)致的機體失調(diào)密切相關(guān)[24,42],臨床可利用和血明目片活血化瘀、清肝明目及滋補肝腎等功效進行治療;角膜血管新生化則屬于“血翳包睛”等范疇[43-44],系體內(nèi)臟腑熱毒壅盛,瘀血凝結(jié)所致,因此和血明目可以瀉火解毒、活血破瘀、滋補肝腎等為用藥依據(jù)。此外,和血明目片針對陰虛肝旺,熱傷絡(luò)脈等證型所引起的眼部病變的療效已經(jīng)得到臨床證實,與翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化的治則存在一定的關(guān)聯(lián),進一步提示和血明目片對于上述適應(yīng)證的潛在治療作用。本研究雖然不能直接確定其潛在適應(yīng)證,但為前期臨床探索指明了方向,提供了線索,在體外實驗的基礎(chǔ)上初步說明和血明目方對于視網(wǎng)膜損傷具有一定的保護作用,一定程度上快速確定實驗驗證的范圍,減少實驗驗證的工作量,節(jié)約時間和成本。后期計劃圍繞其潛在適應(yīng)證的作用機制對其中的具體作用環(huán)節(jié)進行精細(xì)化分析,尋找臨床中具有相似作用機制的西藥進行比較,針對疾病共有環(huán)節(jié)及特殊環(huán)節(jié)的作用靶標(biāo)進行相關(guān)的體內(nèi)外實驗驗證,其中和血明目片干預(yù)4種適應(yīng)證的共有環(huán)節(jié)為炎癥反應(yīng),可以重點關(guān)注其相關(guān)指標(biāo)白細(xì)胞介素6、白細(xì)胞介素1β、腫瘤壞死因子等指標(biāo)進行驗證,而翼狀胬肉、高眼壓癥、角膜血管新生化的特有環(huán)節(jié)為氧化應(yīng)激,則可以通過檢測超氧化物歧化酶、活性氧等金指標(biāo)進行藥效驗證。

    綜上所述,基于整合藥物和靶標(biāo)預(yù)測的相互作用關(guān)系,發(fā)現(xiàn)和血明目片對于翼狀胬肉、角膜炎、高眼壓癥、角膜新生血管化等疾病有較強的干預(yù)作用。隨后通過文獻考察相關(guān)適應(yīng)證的生物學(xué)過程及通路進一步驗證我們的結(jié)論。通過開發(fā)的藥物對于網(wǎng)絡(luò)穩(wěn)健性擾動的評價算法,進行藥物潛在適應(yīng)證的預(yù)測,為藥物的二次開發(fā)和精準(zhǔn)定位提供了思路,更為中藥復(fù)雜機制解析提供參考,希望未來可以應(yīng)用在更多中藥復(fù)方的研究當(dāng)中。

    參考文獻

    [1]孟超.和血明目片的方解及臨床應(yīng)用[A].第十屆全國中西醫(yī)結(jié)合眼科學(xué)術(shù)會議暨第五屆海峽眼科學(xué)術(shù)交流會論文匯編[C].福州:中國中西醫(yī)結(jié)合學(xué)會,2011:83-84.

    [2]王輝,高健生.和血明目片方義及臨床應(yīng)用[J].湖南中醫(yī)藥大學(xué)學(xué)報,2010,30(12):17-18.

    [3]國家藥典委員會.中華人民共和國藥典(一部)[S].北京:中國醫(yī)藥科技出版社,2015:133.

    [4]江一,苗彥.和血明目片應(yīng)用于濕性年齡相關(guān)性黃斑變性的臨床效果研究[J].中國處方藥,2020,18(12):98-99.

    [5]王中華,金暢.和血明目片聯(lián)合血塞通治療非缺血型視網(wǎng)膜靜脈阻塞的臨床觀察[J].臨床醫(yī)藥文獻(連續(xù)型電子期刊),2018,5(33):60-61.

    [6]武靜,楊林青,郝麗莎.和血明目片治療濕性年齡相關(guān)性黃斑變性的臨床觀察[J].中國中醫(yī)眼科雜志,2020,30(8):554-557.

    [7]魏凌君,鄭海濤.和血明目片治療視網(wǎng)膜靜脈阻塞性眼底出血療效分析[J].海峽藥學(xué),2017,29(1):148-150.

    [8]周黎紋.和血明目片在糖尿病視網(wǎng)膜病變性眼底出血治療中的有效運用[J].中國現(xiàn)代藥物應(yīng)用,2018,12(7):90-92.

    [9]Ru J,Li P,Wang J,et al.TCMSP:a database of systems pharmacology for drug discovery from herbal medicines[J].J Cheminform,2014,6:13.

    [10]Xu HY,Zhang YQ,Liu ZM,et al.ETCM:an encyclopaedia of traditional Chinese medicine[J].Nucleic Acids Res,2019,47(D1):D976-D982.

    [11]Liu Z,Guo F,Wang Y,et al.BATMAN-TCM:a Bioinformatics Analysis Tool for Molecular mechANism of Traditional Chinese Medicine[J].Sci Rep,2016,6:21146.

    [12]Groza T,Khler S,Moldenhauer D,et al.The Human Phenotype Ontology:Semantic Unification of Common and Rare Disease[J].Am J Hum Genet,2015,97(1):111-124.

    [13]Piero J,Bravo à,Queralt-Rosinach N,et al.DisGeNET:a comprehensive platform integrating information on human disease-associated genes and variants[J].Nucleic Acids Res,2017,45(D1):D833-D839.

    [14]Zhou Y,Zhou B,Pache L,et al.Metascape provides a biologist-oriented resource for the analysis of systems-level datasets[J].Nat Commun,2019,10(1):1523.

    [15]Szklarczyk D,Morris JH,Cook H,et al.The STRING database in 2017:quality-controlled protein-protein association networks,made broadly accessible[J].Nucleic Acids Res,2017,45(D1):D362-D368.

    [16]Shannon P,Markiel A,Ozier O,et al.Cytoscape:a software environment for integrated models of biomolecular interaction networks[J].Genome Res,2003,13(11):2498-504.

    [17]郭非非,張雨綺,唐仕歡,等.中藥抗COVID-19藥效預(yù)測分析平臺構(gòu)建及常用中藥潛在作用分析[J].中國中藥雜志,2020,45(10):2257-2264.

    [18]Guo F,Zhang W,Su J,et al.Prediction of Drug Positioning for Quan-Du-Zhong Capsules Against Hypertensive Nephropathy Based on the Robustness of Disease Network[J].Front Pharmacol,2019,10:49.

    [19]Null J U,Null L I,Null Y U,et al.iGraph:an incremental data processing system for dynamic graph[J].Frontiers of Computer Science in China,2016,6(3):15.

    [20]張旻昱,龔慕辛,楊洪軍.含藥腸吸收液:一種新的中藥體外藥理實驗方法[J].中草藥,2018,49(15):3457-3462.

    [21]Seto T,Kato T,Nishio M,et al.Erlotinib alone or with bevacizumab as first-line therapy in patients with advanced non-squamous non-small-cell lung cancer harbouring EGFR mutations(JO25567):an open-label,randomised,multicentre,phase 2 study[J].Lancet Oncol,2014,15(11):1236-1244.

    [22]靳懷運,王劍鋒.翼狀胬肉發(fā)病機制研究進展[J].國際眼科雜志,2016,16(6):1080-1083.

    [23]孫毅,鄧應(yīng)平.真菌性角膜炎研究進展[J].國際眼科雜志,2009,9(5):926-928.

    [24]高惠娟.高眼壓癥相關(guān)因素分析、檢查及處理對策[J].醫(yī)學(xué)理論與實踐,2019,32(13):2003-2005.

    [25]陳鶴,陸曉和.不同臨床背景下角膜新生血管發(fā)生機制的研究進展[J].眼科新進展,2019,39(4):394-397.

    [26]孫曉玲,張曉娟,張清泉,等.原發(fā)星形細(xì)胞瘤中TP53、PDGF、EGFR的表達及與臨床病理特征和預(yù)后的相關(guān)性[J].腫瘤防治研究,2021,48(2):154-158.

    [27]楊莉,黃鼎宇,馬蘭,等.鎘激活JNK/c-Jun信號通路誘導(dǎo)小鼠GC-2 spd精母細(xì)胞凋亡[J].公共衛(wèi)生與預(yù)防醫(yī)學(xué),2021,32(3):37-40.

    [28]石勤業(yè),郭劍,徐建紅.核轉(zhuǎn)錄因子κB及其抑制因子研究[J].醫(yī)學(xué)信息,2020,33(22):45-47,54.

    [29]Austin PC,Merlo J.Intermediate and advanced topics in multilevel logistic regression analysis[J].Stat Med,2017,36(20):3257-3277.

    [30]葉宗偉,楊毅寧.血清C反應(yīng)蛋白、白蛋白及其比值在炎癥相關(guān)性疾病中的進展[J].醫(yī)學(xué)綜述,2017,23(20):3979-3983,3988.

    [31]董政,婁煥堃,梅佳慧,等.載脂蛋白E基因多態(tài)性與炎癥的相關(guān)性研究進展[J].中國臨床實用醫(yī)學(xué),2020,11(5):74-77.

    [32]Chechneva OV,Deng W.Mitochondrial translocator protein(TSPO),astrocytes and neuroinflammation[J].Neural Regen Res,2016,11(7):1056-1057.

    [33]郭光明,王明榮.c-Jun和c-Fos及其與人類腫瘤[J].國外醫(yī)學(xué).遺傳學(xué)分冊,2005,28(4):207-211.

    [34]姚翠英,陳芳.炎癥和纖維標(biāo)志物在痤瘡增生性瘢痕形成中的表達[J].基因組學(xué)與應(yīng)用生物學(xué),2020,39(6):2890-2896.

    [35]孔德寧,杜玄一.益氣清熱降濁中藥治療慢性腎臟病3~4期療效及對轉(zhuǎn)化生長因子β1、氧化應(yīng)激指標(biāo)的影響[J].現(xiàn)代中西醫(yī)結(jié)合雜志,2017,26(22):2462-2464.

    [36]汪寧,溫昌明,張保朝.中西醫(yī)結(jié)合對急性心源性腦梗死患者血清PLG和神經(jīng)功能的影響[J].世界中醫(yī)藥,2019,14(4):895-898.

    [37]官勁帆,韓江全,田明巧,等.PPARγ與腦卒中關(guān)系的研究進展[J].齊齊哈爾醫(yī)學(xué)院學(xué)報,2021,42(6):494-498.

    [38]陳靜,李威.翼狀胬肉的發(fā)病機制及治療進展[J].現(xiàn)代中西醫(yī)結(jié)合雜志,2020,29(12):1364-1368.

    [39]曾志成,彭清華.退翳明目湯對翼狀胬肉切除術(shù)后基礎(chǔ)淚液分泌和淚膜穩(wěn)定性的影響[J].中醫(yī)雜志,2014,55(3):218-221.

    [40]田曼,劉昳.單純皰疹病毒性角膜炎中西醫(yī)研究進展[J].遼寧中醫(yī)藥大學(xué)學(xué)報,2019,21(4):215-221.

    [41]湯佳羽,蹇文淵.中醫(yī)藥治療單純皰疹病毒性角膜炎用藥規(guī)律分析[J].河南中醫(yī),2021,41(2):245-249.

    [42]李苑碧,彭清華.原發(fā)性開角型青光眼發(fā)病機制及其中西醫(yī)治療研究進展[J].遼寧中醫(yī)藥大學(xué)學(xué)報,2014,16(8):146-149.

    [43]莊家圓.角膜新生血管形成機制及治療的研究新進展[J].醫(yī)學(xué)食療與健康,2020,18(18):209-211.

    [44]王軍花,高桂平.角膜新生血管的發(fā)生與治療[J].中國組織工程研究與臨床康復(fù),2010,14(33):6223-6225.

    (2022-03-10收稿 本文編輯:王明)

    成人欧美大片| 欧美性猛交╳xxx乱大交人| 精品久久久久久久末码| 亚洲国产成人一精品久久久| 国产欧美亚洲国产| 亚洲精品乱码久久久久久按摩| 一本色道久久久久久精品综合| 日本av手机在线免费观看| 在现免费观看毛片| 国产综合懂色| 欧美激情国产日韩精品一区| av在线观看视频网站免费| 精品久久久久久久久av| 国产成人aa在线观看| 国产在线一区二区三区精| 日韩国内少妇激情av| 亚洲精品乱码久久久久久按摩| 国产亚洲91精品色在线| 欧美日韩国产mv在线观看视频 | 小蜜桃在线观看免费完整版高清| 人妻夜夜爽99麻豆av| 亚洲av.av天堂| 午夜福利视频精品| 欧美3d第一页| 国产探花极品一区二区| 哪个播放器可以免费观看大片| 舔av片在线| 熟女人妻精品中文字幕| 亚洲人成网站在线观看播放| 日日啪夜夜撸| 18禁在线播放成人免费| 三级国产精品片| a级一级毛片免费在线观看| 伊人久久精品亚洲午夜| 国产黄频视频在线观看| 青青草视频在线视频观看| 99热6这里只有精品| 国产日韩欧美在线精品| a级一级毛片免费在线观看| 成年女人看的毛片在线观看| 亚洲国产精品专区欧美| 性色avwww在线观看| 亚洲不卡免费看| 日本熟妇午夜| 黄片wwwwww| av一本久久久久| 在线观看美女被高潮喷水网站| 成人免费观看视频高清| 黄色日韩在线| 午夜福利在线观看免费完整高清在| 亚洲精品,欧美精品| 国产av码专区亚洲av| 下体分泌物呈黄色| 麻豆成人午夜福利视频| 午夜福利视频精品| 毛片女人毛片| 热99国产精品久久久久久7| 九九爱精品视频在线观看| 少妇人妻一区二区三区视频| videos熟女内射| 日韩伦理黄色片| 神马国产精品三级电影在线观看| 精品酒店卫生间| 日本三级黄在线观看| 在线 av 中文字幕| 综合色丁香网| 亚洲怡红院男人天堂| 国产精品.久久久| 97精品久久久久久久久久精品| 国产视频首页在线观看| 蜜臀久久99精品久久宅男| 亚洲经典国产精华液单| 成人美女网站在线观看视频| 久久久色成人| 国产成人a∨麻豆精品| av一本久久久久| 深夜a级毛片| 国产精品一区二区性色av| 精品少妇黑人巨大在线播放| 国产综合懂色| 欧美3d第一页| 99久久精品一区二区三区| 免费观看无遮挡的男女| 亚洲国产欧美人成| 国内揄拍国产精品人妻在线| 亚洲av欧美aⅴ国产| 99精国产麻豆久久婷婷| 亚洲精华国产精华液的使用体验| 中文字幕人妻熟人妻熟丝袜美| 国产精品女同一区二区软件| 一级毛片 在线播放| 久久99精品国语久久久| 男女边摸边吃奶| 人妻夜夜爽99麻豆av| 99视频精品全部免费 在线| av在线亚洲专区| 精品人妻一区二区三区麻豆| 国产精品一二三区在线看| 日韩一区二区三区影片| 国产美女午夜福利| 亚洲欧美清纯卡通| 欧美性感艳星| 色哟哟·www| 免费播放大片免费观看视频在线观看| 精品人妻偷拍中文字幕| 黄色日韩在线| 国产成人一区二区在线| 大又大粗又爽又黄少妇毛片口| 干丝袜人妻中文字幕| 人妻少妇偷人精品九色| 免费观看性生交大片5| 亚洲综合色惰| 日韩欧美 国产精品| 插逼视频在线观看| 国产精品99久久久久久久久| 麻豆成人av视频| 熟女av电影| 精品久久久噜噜| 性色av一级| 亚洲国产精品成人久久小说| 91久久精品国产一区二区成人| 丰满乱子伦码专区| 国产精品一及| 日本与韩国留学比较| 街头女战士在线观看网站| 欧美+日韩+精品| 国产日韩欧美亚洲二区| 97在线人人人人妻| 免费看av在线观看网站| 午夜福利视频精品| 一级毛片电影观看| 国产视频首页在线观看| 在线观看美女被高潮喷水网站| 久久久久久久久大av| 汤姆久久久久久久影院中文字幕| 精品酒店卫生间| 国产精品国产av在线观看| 国产乱来视频区| 国产人妻一区二区三区在| 丰满少妇做爰视频| 精品久久久久久久久av| 精品人妻偷拍中文字幕| 美女主播在线视频| 男的添女的下面高潮视频| 日本免费在线观看一区| 成年版毛片免费区| 十八禁网站网址无遮挡 | 亚洲内射少妇av| 中国三级夫妇交换| 久久久久久久久大av| 搡女人真爽免费视频火全软件| 成人黄色视频免费在线看| 男的添女的下面高潮视频| 国产高清不卡午夜福利| 国产精品.久久久| 天堂中文最新版在线下载 | 少妇的逼水好多| 国产成人a区在线观看| 又黄又爽又刺激的免费视频.| videos熟女内射| 老师上课跳d突然被开到最大视频| 伊人久久国产一区二区| 乱码一卡2卡4卡精品| 99热这里只有是精品50| 黄色欧美视频在线观看| 成人午夜精彩视频在线观看| 80岁老熟妇乱子伦牲交| 久久久欧美国产精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲人成网站高清观看| a级一级毛片免费在线观看| 哪个播放器可以免费观看大片| 国产乱人视频| 国产女主播在线喷水免费视频网站| 国产精品久久久久久久电影| 色综合色国产| 黄色配什么色好看| 99精国产麻豆久久婷婷| 国产有黄有色有爽视频| 色网站视频免费| 成人亚洲精品一区在线观看 | 精品少妇久久久久久888优播| 亚洲aⅴ乱码一区二区在线播放| 色婷婷久久久亚洲欧美| 国产亚洲一区二区精品| 国产 精品1| 简卡轻食公司| 久久精品国产亚洲网站| 秋霞在线观看毛片| 乱系列少妇在线播放| 久久久久九九精品影院| 亚洲精品国产色婷婷电影| 最近最新中文字幕免费大全7| 国产精品三级大全| 日本午夜av视频| 国内精品美女久久久久久| 一个人看的www免费观看视频| a级一级毛片免费在线观看| 高清视频免费观看一区二区| 网址你懂的国产日韩在线| 国产精品一区二区在线观看99| 波野结衣二区三区在线| 你懂的网址亚洲精品在线观看| 欧美区成人在线视频| 国产精品一区二区性色av| 日韩人妻高清精品专区| 神马国产精品三级电影在线观看| 精品国产乱码久久久久久小说| 日本色播在线视频| 爱豆传媒免费全集在线观看| 国产精品女同一区二区软件| 午夜福利视频精品| 午夜老司机福利剧场| 亚洲色图综合在线观看| 久久6这里有精品| 一区二区三区免费毛片| 国产在视频线精品| 简卡轻食公司| 久久精品久久久久久久性| 香蕉精品网在线| 极品少妇高潮喷水抽搐| 日本一本二区三区精品| 少妇的逼水好多| 十八禁网站网址无遮挡 | 亚洲内射少妇av| 18禁在线播放成人免费| 一本色道久久久久久精品综合| 大码成人一级视频| 国产一区二区亚洲精品在线观看| 国产成人a∨麻豆精品| 日本av手机在线免费观看| 日日撸夜夜添| 久久精品综合一区二区三区| 天美传媒精品一区二区| 大片免费播放器 马上看| 在现免费观看毛片| 婷婷色av中文字幕| av在线app专区| 欧美老熟妇乱子伦牲交| 女人久久www免费人成看片| 在线天堂最新版资源| 欧美区成人在线视频| 男人舔奶头视频| 99热网站在线观看| 校园人妻丝袜中文字幕| 天美传媒精品一区二区| 69av精品久久久久久| 久久久久网色| 久久久久网色| 国产精品久久久久久久久免| 国产91av在线免费观看| 夜夜爽夜夜爽视频| 另类亚洲欧美激情| 成人高潮视频无遮挡免费网站| 亚洲国产高清在线一区二区三| 免费观看av网站的网址| 欧美日韩一区二区视频在线观看视频在线 | 99久久精品热视频| 爱豆传媒免费全集在线观看| 国产精品精品国产色婷婷| 亚洲国产精品成人久久小说| 特大巨黑吊av在线直播| 97在线人人人人妻| 男女无遮挡免费网站观看| 欧美日本视频| 国产精品一区二区三区四区免费观看| 国产精品一区二区在线观看99| 777米奇影视久久| 26uuu在线亚洲综合色| 五月伊人婷婷丁香| 精品午夜福利在线看| 精品久久久久久久久av| 性色avwww在线观看| 欧美极品一区二区三区四区| 国产欧美日韩一区二区三区在线 | 国产69精品久久久久777片| 久久人人爽av亚洲精品天堂 | 国产亚洲午夜精品一区二区久久 | 最近最新中文字幕大全电影3| 国产免费又黄又爽又色| 色5月婷婷丁香| 日本免费在线观看一区| 日日啪夜夜撸| 国产伦精品一区二区三区视频9| 日本猛色少妇xxxxx猛交久久| 国模一区二区三区四区视频| 中文乱码字字幕精品一区二区三区| 成人鲁丝片一二三区免费| 男女那种视频在线观看| 18+在线观看网站| 国产高清国产精品国产三级 | 国产欧美另类精品又又久久亚洲欧美| 春色校园在线视频观看| 久久精品国产鲁丝片午夜精品| 国产精品一区二区三区四区免费观看| 日本黄大片高清| 欧美亚洲 丝袜 人妻 在线| 亚洲一级一片aⅴ在线观看| 一级毛片久久久久久久久女| 国产亚洲精品久久久com| 如何舔出高潮| 精品久久久久久久人妻蜜臀av| 人人妻人人澡人人爽人人夜夜| 天天躁夜夜躁狠狠久久av| 菩萨蛮人人尽说江南好唐韦庄| 噜噜噜噜噜久久久久久91| 一本久久精品| 日韩成人av中文字幕在线观看| 亚洲va在线va天堂va国产| 在线免费观看不下载黄p国产| 成年av动漫网址| 一级黄片播放器| 内地一区二区视频在线| a级毛色黄片| 日韩一区二区三区影片| 青春草视频在线免费观看| 黄色一级大片看看| 亚洲精品乱码久久久久久按摩| 高清av免费在线| 亚洲av电影在线观看一区二区三区 | 一级片'在线观看视频| 一二三四中文在线观看免费高清| 午夜爱爱视频在线播放| 久久久久久久久久久丰满| 97超碰精品成人国产| 国产精品一区二区性色av| 日韩中字成人| 男女边吃奶边做爰视频| 成人无遮挡网站| 国产中年淑女户外野战色| 新久久久久国产一级毛片| 精品视频人人做人人爽| 国产精品一区二区在线观看99| 黄片wwwwww| 亚洲精品国产成人久久av| 国产成人免费无遮挡视频| 亚洲精品第二区| 亚洲精品第二区| 高清在线视频一区二区三区| 日韩,欧美,国产一区二区三区| 两个人的视频大全免费| 3wmmmm亚洲av在线观看| 嫩草影院精品99| 2021少妇久久久久久久久久久| 女的被弄到高潮叫床怎么办| 91狼人影院| 人妻夜夜爽99麻豆av| 我的老师免费观看完整版| 国产精品久久久久久av不卡| 1000部很黄的大片| 内射极品少妇av片p| 亚洲丝袜综合中文字幕| 我的女老师完整版在线观看| 国产成人精品婷婷| 午夜激情久久久久久久| 91久久精品国产一区二区成人| 免费在线观看成人毛片| 国产色婷婷99| 免费黄网站久久成人精品| 欧美3d第一页| 国产v大片淫在线免费观看| 国产一级毛片在线| 韩国高清视频一区二区三区| 可以在线观看毛片的网站| 天堂中文最新版在线下载 | 日本色播在线视频| 黑人高潮一二区| 国产精品久久久久久av不卡| 中文字幕制服av| 美女脱内裤让男人舔精品视频| 国产欧美另类精品又又久久亚洲欧美| 男插女下体视频免费在线播放| 国产一区有黄有色的免费视频| 亚洲电影在线观看av| 国产精品国产三级国产av玫瑰| 欧美xxxx黑人xx丫x性爽| 国产爽快片一区二区三区| 少妇高潮的动态图| 最新中文字幕久久久久| 男男h啪啪无遮挡| 午夜精品国产一区二区电影 | 内射极品少妇av片p| 蜜臀久久99精品久久宅男| 日韩强制内射视频| 亚洲av不卡在线观看| 人体艺术视频欧美日本| 丰满乱子伦码专区| 欧美xxxx性猛交bbbb| 亚洲欧美成人精品一区二区| 一区二区三区精品91| 99九九线精品视频在线观看视频| 天堂俺去俺来也www色官网| 色5月婷婷丁香| 亚洲天堂av无毛| 久久97久久精品| 成人亚洲精品av一区二区| 成人亚洲精品一区在线观看 | 赤兔流量卡办理| 男女啪啪激烈高潮av片| 国产成人精品婷婷| 热99国产精品久久久久久7| av在线老鸭窝| 精品少妇黑人巨大在线播放| 国产黄频视频在线观看| 免费大片18禁| 久久精品国产亚洲网站| 亚洲人成网站在线观看播放| 22中文网久久字幕| 国产黄a三级三级三级人| 亚洲精品日本国产第一区| 国产探花在线观看一区二区| 成人毛片a级毛片在线播放| 有码 亚洲区| 少妇 在线观看| 日日摸夜夜添夜夜添av毛片| 亚洲欧美精品专区久久| 黄色怎么调成土黄色| 国产欧美日韩一区二区三区在线 | 国产淫语在线视频| 久久久午夜欧美精品| 我的女老师完整版在线观看| 卡戴珊不雅视频在线播放| 国产欧美日韩一区二区三区在线 | 丝袜喷水一区| 看十八女毛片水多多多| 3wmmmm亚洲av在线观看| 午夜老司机福利剧场| 精品久久久噜噜| 久久久午夜欧美精品| 国产av国产精品国产| 高清在线视频一区二区三区| 1000部很黄的大片| 美女xxoo啪啪120秒动态图| 久久久久久久国产电影| 在线a可以看的网站| 一级毛片 在线播放| 亚洲最大成人手机在线| 日韩电影二区| 街头女战士在线观看网站| 国产免费一区二区三区四区乱码| 成人综合一区亚洲| 精品人妻视频免费看| 91在线精品国自产拍蜜月| 欧美成人午夜免费资源| 国产午夜福利久久久久久| 中国国产av一级| 亚洲av在线观看美女高潮| 免费大片黄手机在线观看| 久久99热这里只有精品18| 中文乱码字字幕精品一区二区三区| 色5月婷婷丁香| 亚洲真实伦在线观看| 啦啦啦中文免费视频观看日本| 亚洲不卡免费看| av在线蜜桃| 日韩强制内射视频| 99热这里只有精品一区| 亚洲av在线观看美女高潮| 国产精品久久久久久久久免| 99久久九九国产精品国产免费| 久久99热这里只有精品18| 成年人午夜在线观看视频| 日韩在线高清观看一区二区三区| 亚洲不卡免费看| 成人亚洲欧美一区二区av| 久久久久久伊人网av| 人妻少妇偷人精品九色| 18+在线观看网站| 亚洲成人精品中文字幕电影| 久久午夜福利片| 美女被艹到高潮喷水动态| 国产精品成人在线| 欧美成人a在线观看| 大又大粗又爽又黄少妇毛片口| 男女下面进入的视频免费午夜| 亚洲国产最新在线播放| 在线播放无遮挡| 久久久成人免费电影| 欧美xxxx黑人xx丫x性爽| 久热久热在线精品观看| av在线亚洲专区| 国产黄色免费在线视频| 丰满乱子伦码专区| 男人爽女人下面视频在线观看| 亚洲av成人精品一区久久| 久久女婷五月综合色啪小说 | 岛国毛片在线播放| 亚洲欧美日韩无卡精品| 精品一区在线观看国产| 好男人在线观看高清免费视频| 中文字幕亚洲精品专区| 亚洲,欧美,日韩| 国产乱人视频| 久久人人爽人人爽人人片va| 久久国内精品自在自线图片| 联通29元200g的流量卡| 一区二区三区精品91| 精品久久久精品久久久| 午夜免费男女啪啪视频观看| 亚洲欧美日韩无卡精品| 嫩草影院新地址| 日本一本二区三区精品| 亚洲精品日韩在线中文字幕| 丰满乱子伦码专区| 天堂俺去俺来也www色官网| 亚洲精华国产精华液的使用体验| 一区二区三区精品91| 99热6这里只有精品| 国产永久视频网站| 乱系列少妇在线播放| 夫妻性生交免费视频一级片| 欧美性猛交╳xxx乱大交人| 精品一区二区三卡| 精品国产三级普通话版| 青春草国产在线视频| 制服丝袜香蕉在线| 91精品伊人久久大香线蕉| 亚洲精品一二三| 国产亚洲91精品色在线| av在线app专区| 久久国产乱子免费精品| 简卡轻食公司| 一个人看的www免费观看视频| 在线观看国产h片| 日韩视频在线欧美| 菩萨蛮人人尽说江南好唐韦庄| 三级国产精品片| 欧美xxⅹ黑人| 天天躁日日操中文字幕| 国产免费又黄又爽又色| 国产精品av视频在线免费观看| 久久99热6这里只有精品| 我的老师免费观看完整版| 国产免费视频播放在线视频| 国产午夜精品一二区理论片| 日韩三级伦理在线观看| 亚洲精品乱码久久久v下载方式| 国产黄色免费在线视频| 五月天丁香电影| 人体艺术视频欧美日本| 国产av国产精品国产| 99re6热这里在线精品视频| 国产精品精品国产色婷婷| 久久这里有精品视频免费| 亚洲精品乱久久久久久| 国产亚洲av片在线观看秒播厂| 男女下面进入的视频免费午夜| 美女高潮的动态| 亚洲精品久久久久久婷婷小说| 在线免费十八禁| 国产爱豆传媒在线观看| 2021少妇久久久久久久久久久| 99久久九九国产精品国产免费| 青青草视频在线视频观看| 人人妻人人看人人澡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 男人狂女人下面高潮的视频| 91午夜精品亚洲一区二区三区| 99久国产av精品国产电影| 99视频精品全部免费 在线| 六月丁香七月| 国产欧美亚洲国产| 精品少妇久久久久久888优播| 91狼人影院| 18禁在线播放成人免费| 人人妻人人看人人澡| 久久精品综合一区二区三区| 国产日韩欧美在线精品| 亚洲av福利一区| 日本一本二区三区精品| 我的女老师完整版在线观看| 精品亚洲乱码少妇综合久久| 日韩伦理黄色片| 国产精品一二三区在线看| 最近中文字幕2019免费版| 亚洲国产欧美人成| 18禁动态无遮挡网站| 国产成人精品久久久久久| 嫩草影院精品99| 精品人妻一区二区三区麻豆| 国产永久视频网站| 欧美zozozo另类| 成人综合一区亚洲| 男女下面进入的视频免费午夜| 麻豆成人午夜福利视频| 国产成人freesex在线| 97超视频在线观看视频| 国产视频内射| 九九在线视频观看精品| 免费大片18禁| 在线观看一区二区三区激情| freevideosex欧美| 一区二区三区免费毛片| 久久久久久久久久成人| 一边亲一边摸免费视频| 亚洲精品国产av蜜桃| 91久久精品电影网| 免费观看无遮挡的男女| 欧美潮喷喷水| 99久久精品热视频| 狂野欧美白嫩少妇大欣赏| 三级男女做爰猛烈吃奶摸视频| 亚洲内射少妇av| 欧美 日韩 精品 国产| 午夜福利网站1000一区二区三区| 69av精品久久久久久| 午夜免费观看性视频| 精品国产三级普通话版| 国产成人aa在线观看| 国产高清不卡午夜福利| 搡老乐熟女国产| 日日啪夜夜爽| 亚洲欧美日韩无卡精品| 亚洲精品日本国产第一区| 女人久久www免费人成看片|