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

    基于特征圖譜、化學(xué)計量學(xué)和分子對接的復(fù)方阿膠漿質(zhì)量標(biāo)志物研究

    2021-12-04 02:22:04任雪陽宋若蘭于啊香馬嘉慕折改梅
    中草藥 2021年23期
    關(guān)鍵詞:糖苷黨參批號

    許 嘯,張 淹,任雪陽,王 宇,董 英,宋若蘭,于啊香,魏 靜,馬嘉慕,折改梅*

    基于特征圖譜、化學(xué)計量學(xué)和分子對接的復(fù)方阿膠漿質(zhì)量標(biāo)志物研究

    許 嘯1,張 淹2,任雪陽1,王 宇1,董 英1,宋若蘭1,于啊香1,魏 靜1,馬嘉慕1,折改梅1*

    1. 北京中醫(yī)藥大學(xué)中藥學(xué)院,北京 102488 2. 山東東阿阿膠股份有限公司國家膠類中藥工程技術(shù)研究中心,山東 聊城 252201

    以質(zhì)量標(biāo)志物(quality markers,Q-Marker)“五原則”為指導(dǎo),結(jié)合特征圖譜、化學(xué)計量學(xué)、分子對接和含量測定等技術(shù),開展復(fù)方阿膠漿(Fufang E’jiao Jiang,F(xiàn)EJ)治療貧血癥的Q-Marker研究。色譜柱Agilent Zorbax SB AQ-C18(250 mm×4.6 mm,5 μm),流動相為0.1%甲酸乙腈溶液-0.1%甲酸水溶液,體積流量1.0 mL/min,柱溫30 ℃,檢測波長330 nm,焦地黃苯乙醇苷A1為參照峰,建立FEJ特征圖譜。采用層次聚類分析(hierarchical clustering analysis,HCA)、主成分分析(principal component analysis,PCA)、偏最小二乘法-判別分析(partial least square method-discriminant analysis,PLS-DA)對不同批次FEJ進(jìn)行分析。查閱《中國藥典》和文獻(xiàn),結(jié)合組方配伍環(huán)境(君、臣、佐使)、成分傳遞與溯源等多因素,收集篩選質(zhì)量控制6個關(guān)鍵成分。在GEO數(shù)據(jù)庫中檢索篩選貧血癥(anemia)的差異基因,采用Discovery Studio 2016 v16.1軟件,將成分與貧血癥相關(guān)靶點進(jìn)行對接驗證。以330、270 nm雙波長同時測定FEJ中異毛蕊花糖苷和黨參炔苷的含量。FEJ特征圖譜中指認(rèn)了咖啡酸、綠原酸、洋地黃葉苷C、阿魏酸、焦地黃苯乙醇苷A1、毛蕊花糖苷、焦地黃苯乙醇苷B1、異毛蕊花糖苷共8個特征成分。HCA、PCA和PLS-DA將16批FEJ分為3類,其中異毛蕊花糖苷、洋地黃葉苷C和焦地黃苯乙醇苷A1對FEJ質(zhì)量影響最為顯著。8個特征成分和6個關(guān)鍵成分的分子對接結(jié)果表明,異毛蕊花糖苷、黨參炔苷、綠原酸、毛蕊花糖苷4個成分能夠與激肽原-1(kininogen-1,KNG1)、泛素樣蛋白ISG15(ubiquitin-like protein ISG15,ISG15)、賴氨酸特異性脫甲基酶6A(lysine-specific demethylase 6A,KDM6A)、乙酰化酶動力蛋白-3(dynamin-3,DNM3)、蛋白酶激活受體2(proteinase-activated receptor 2,F(xiàn)2RL1)、半胱氨酰白三烯受體1(cysteinyl leukotriene receptor 1,CYSLTR1)6個靶蛋白通過氫鍵、疏水鍵、π-π鍵等作用力結(jié)合,具有較好的活性,可作為FEJ的Q-Marker。16批樣品中異毛蕊花糖苷、黨參炔苷質(zhì)量濃度分別為6.05~12.30、9.14~16.30 μg/mL。建立了FEJ特征圖譜,結(jié)合化學(xué)計量學(xué)和分子對接等技術(shù),篩選了4個Q-Marker異毛蕊花糖苷、黨參炔苷、毛蕊花糖苷、綠原酸,并建立了異毛蕊花糖苷、黨參炔苷的含量測定方法,為FEJ質(zhì)量控制提供參考。

    復(fù)方阿膠漿;質(zhì)量標(biāo)志物;特征圖譜;化學(xué)計量學(xué);分子對接;貧血癥;層次聚類分析;主成分分析;偏最小二乘法-判別分析;質(zhì)量控制;咖啡酸;綠原酸;洋地黃葉苷C;阿魏酸;焦地黃苯乙醇苷A1;毛蕊花糖苷;焦地黃苯乙醇苷B1;異毛蕊花糖苷

    復(fù)方阿膠漿(Fufang E’jiao Jiang,F(xiàn)EJ)包括阿膠、紅參、熟地黃、黨參、山楂共5味中藥。原方出自明代《景岳全書》中“兩儀膏”,按照“氣血互生,氣升血長”的中醫(yī)理論加味而成。臨床常用于氣血兩虛、頭暈?zāi)垦!⑿募率?、食欲不振及白?xì)胞減少癥和貧血等癥[1]?!吨袊幍洹?020年版FEJ質(zhì)量標(biāo)準(zhǔn)重點關(guān)注了阿膠的真?zhèn)魏蛢?yōu)劣。在中醫(yī)理論中,氣血相關(guān)理論與治療血管生成相關(guān)?,F(xiàn)代藥理研究表明,血管新生是涉及多種細(xì)胞如內(nèi)皮細(xì)胞和血管平滑肌細(xì)胞的多步驟復(fù)雜生物過程,是相關(guān)細(xì)胞因子、細(xì)胞外基質(zhì)等與細(xì)胞之間相互作用的結(jié)果。在缺血刺激下,新的血管和側(cè)支循環(huán)會逐漸生成,達(dá)到改善缺血的目的[2-3]。文獻(xiàn)報道,斑馬魚模型可用于研究血管生成[4-5]、微血管功能[6]等方面。本課題組前期研究發(fā)現(xiàn)FEJ的促血管生成活性優(yōu)于單味藥阿膠。因此紅參等4味植物來源中藥對FEJ功效有很重要的貢獻(xiàn),質(zhì)量標(biāo)準(zhǔn)應(yīng)該對其予以關(guān)注。

    2016年,劉昌孝院士[7-8]提出了中藥質(zhì)量標(biāo)志物(quality marker,Q-Marker)這一中藥質(zhì)量評價與質(zhì)量控制的核心概念,“五原則”從質(zhì)量傳遞與溯源、成分特有性、成分的有效性、成分可測性以及復(fù)方配伍環(huán)境5個方面,貫穿并指導(dǎo)了中藥Q-Marker的研究與發(fā)現(xiàn)[9]。中藥特征圖譜是一種綜合的、可量化的鑒別手段,能夠反映中藥多成分特點,已被廣泛應(yīng)用于中藥及中藥復(fù)方制劑的質(zhì)量控制[10-11]。將化學(xué)計量學(xué)方法如層次聚類分析(hierarchical cluster analysis,HCA)、主成分分析(principal components analysis,PCA)法、偏最小二乘法-判別分析(partial least squares-discrimination analysis,PLS-DA)等用于特征圖譜的數(shù)據(jù)分析可以更為有效的挖掘信息。有文獻(xiàn)報道[12-14],采用HCA分類結(jié)果,結(jié)合PCA和PLS-DA得到各成分的貢獻(xiàn)率,從而挖掘造成組間差異的主要成分。這些成分可以作為Q-Marker的候選成分。

    但是由于特征圖譜中實驗條件(如供試品制備方法、檢測波長等)的局限,存在不能全部體現(xiàn)化學(xué)成分信息的問題。這個也會影響到潛在Q-Marker的選擇,或者會造成潛在Q-Marker的丟失。因此,本研究又通過查閱《中國藥典》2020年版和參考文獻(xiàn),整理得到FEJ各藥味中活性成分,并從中藥組方的配伍環(huán)境(君、臣、佐、使)、成分的傳遞和溯源等方面對各單味藥材飲片進(jìn)行分析,從而有效地彌補特征圖譜成分信息不全,獲取質(zhì)量控制的關(guān)鍵成分。從而滿足Q-Marker五原則的質(zhì)量傳遞與溯源、成分特有性以及復(fù)方配伍環(huán)境等特征。

    分子對接技術(shù)是以結(jié)構(gòu)分子生物學(xué)和計算機輔助藥物設(shè)計相結(jié)合的方式,通過將配體與受體進(jìn)行對接,評價其對接體系的穩(wěn)定性,篩選與受體活性部位空間和電性特征相匹配的小分子化合物。該技術(shù)不僅周期短,操作性強,而且還避免了藥理實驗穩(wěn)定性差的缺點,已經(jīng)逐漸成為中藥活性成分篩選的重要手段[15-16]。本研究利用分子對接技術(shù),對FEJ中的特征成分和關(guān)鍵成分進(jìn)行治療貧血癥的活性篩選,從多種角度探索藥物與疾病的關(guān)聯(lián),同時確保了所選Q-Marker的有效性。最后,以Q-Marker為指標(biāo)成分,建立FEJ中熟地黃、黨參的多成分含量測定方法,為FEJ Q-Marker可測性提供了依據(jù)。該研究科學(xué)合理地提升了FEJ的質(zhì)量標(biāo)準(zhǔn),彌補了以總氮量為唯一含量測定方法。本研究的技術(shù)路線如圖1。

    圖1 FEJ Q-Marker研究的技術(shù)路線

    1 儀器與材料

    1.1 儀器

    華譜科儀Acchorm S6000高效液相色譜儀,配備在線脫氣機、四元泵、自動進(jìn)樣器、柱溫箱、DAD檢測器,華譜科儀(北京)科技有限公司;BT25S型十萬分之一電子天平,北京賽多利斯科學(xué)儀器有限公司;MX-S型渦旋振蕩器,北京大龍興創(chuàng)實驗儀器有限公司;3-18N型高速離心機,湖南恒諾儀器設(shè)備有限公司;8302型恒溫水浴鍋,上海利浦自動化儀表廠。

    1.2 材料與試劑

    FEJ(批號1905075、1907011、1907012、1907014、1907015、1907017、1907020、1907021、1907022、1909004、1909007、1909010、1909011、1909013、1909015、1909016,編號依次對應(yīng)為S1~S16),以及批號為1905075批次所對應(yīng)的紅參藥材(批號201812-001)、黨參藥材(批號201905-002)、熟地黃藥材(批號201905-005)、山楂藥材(批號201905-002)、紅參陰性樣品、黨參陰性樣品、熟地黃陰性樣品、山楂陰性樣品,均由東阿阿膠股份有限公司提供。對照品咖啡酸(批號110885-200102)、綠原酸(批號110753-200413)、阿魏酸(批號110773-201915),購于中國食品藥品檢定研究院;對照品焦地黃苯乙醇苷A1(批號DST190902-070)、毛蕊花糖苷(批號DST190422-061)、異毛蕊花糖苷(批號DST200513-060)、焦地黃苯乙醇苷B1(批號DST201203-071)、洋地黃葉苷C(批號DST201109-057),購于成都德思特生物技術(shù)有限公司;對照品黨參炔苷(批號PRF9121002),購自成都普瑞法科技開發(fā)有限公司;以上所有對照品的質(zhì)量分?jǐn)?shù)均≥95%。甲醇、乙腈為色譜級,美國Fishier Scientific公司,其余試劑均為分析純。超純水由實驗室Millipore Synergy超純水系統(tǒng)自制。

    2 方法與結(jié)果

    2.1 FEJ HPLC特征圖譜

    2.1.1 色譜條件 Agilent Zorbax SB-AQ C18色譜柱(250 mm×4.6 mm,5 μm),以0.1%甲酸乙腈溶液為流動相A,以0.1%甲酸水溶液為流動相B,梯度洗脫程序:0~20 min,2%~5% A;20~35 min,5%~12% A;35~45 min,12%~14% A;45~55 min,14%~20% A;55~65 min,20%~23% A;65~80 min,23%~35% A;體積流量1 mL/min;柱溫30 ℃;檢測波長330 nm;進(jìn)樣量10 μL。

    2.1.2 供試品溶液的制備 精密量取FEJ 2 mL,加入乙腈6 mL,渦旋至充分混勻,離心,取上清液,殘渣加入乙腈6 mL,渦旋至充分混勻,離心,合并上清液,蒸干,殘渣加30%甲醇溶解,并轉(zhuǎn)移至1 mL量瓶中,加30%甲醇至刻度,搖勻,即得供試品溶液。

    2.1.3 對照品溶液的制備 取咖啡酸、綠原酸、洋地黃葉苷C、阿魏酸、焦地黃苯乙醇苷A1、毛蕊花糖苷、焦地黃苯乙醇苷B1、異毛蕊花糖苷對照品適量,精密稱定,加甲醇制成分別含咖啡酸20 μg/mL、綠原酸20 μg/mL、洋地黃葉苷C 50 μg/mL、阿魏酸20 μg/mL、焦地黃苯乙醇苷A150 μg/mL、毛蕊花糖苷20 μg/mL、焦地黃苯乙醇苷B150 μg/mL、異毛蕊花糖苷50 μg/mL的混合溶液,即得混合對照品溶液。

    2.1.4 精密度試驗 取編號為S1的FEJ(批號1905075),按“2.1.2”項下方法制得供試品溶液,按“2.1.1”項下色譜條件連續(xù)進(jìn)樣6次,記錄色譜峰峰面積。以焦地黃苯乙醇苷A1為參照峰,各共有峰的相對峰面積RSD均小于1.70%,相對保留時間RSD均小于0.05%,表明儀器精密度良好。

    2.1.5 重復(fù)性試驗 取編號為S1的FEJ(批號1905075)6份,按“2.1.2”項下方法制得供試品溶液,按“2.1.1”項下色譜條件分別進(jìn)樣,記錄峰面積。以焦地黃苯乙醇苷A1為參照峰,各共有峰的相對峰面積RSD均小于2.99%,相對保留時間RSD均小于0.03%,表明方法重復(fù)性良好。

    2.1.6 穩(wěn)定性試驗 取編號為S1的FEJ(批號1905075),按“2.1.2”項下方法制得供試品溶液,按“2.1.1”項下色譜條件,分別在0、2、4、6、8、16、24 h進(jìn)樣,記錄峰面積。以焦地黃苯乙醇苷A1為參照峰,各共有峰的相對峰面積RSD均小于2.08%,相對保留時間RSD均小于0.09%,表明供試品溶液在24 h內(nèi)穩(wěn)定性良好。

    2.1.7 HPLC特征圖譜的生成與相似度評價 取16批FEJ(S1~S16),按“2.1.2”項下方法制得供試品溶液,按“2.1.1”項下色譜條件分別進(jìn)樣。采用《中藥指紋圖譜相似度評價評價系統(tǒng)(2012A版)》對所記錄圖譜的AIA文件進(jìn)行分析。以S1為參照圖譜,時間窗寬度設(shè)為0.2 min,對色譜峰進(jìn)行多點校正和全峰匹配,生成16批FEJ的特征圖譜匹配圖以及對照指紋圖譜,見圖2。

    其中,5號峰焦地黃苯乙醇苷A1分離度好,保留時間居中,選擇其為參照峰,計算S1~S16特征圖譜與對照特征圖譜的相似度,結(jié)果分別為0.998、0.994、0.996、0.998、0.998、0.997、0.998、0.994、0.993、0.998、0.997、0.994、0.989、0.991、0.996、0.997,均大于0.9,表明各批次樣品的化學(xué)成分穩(wěn)定。

    圖2 16批FEJHPLC特征圖譜 (S1~S16) 及其對照特征圖譜(R)

    2.1.8 共有峰的歸屬與指認(rèn) 根據(jù)對照品共指認(rèn)出8個化合物,1~8號共有峰分別為咖啡酸、綠原酸、洋地黃葉苷C、阿魏酸、焦地黃苯乙醇苷A1、毛蕊花糖苷、焦地黃苯乙醇苷B1、異毛蕊花糖苷,結(jié)果見圖3。

    取FEJ(批號1905075)及對應(yīng)批次的單味藥材飲片,按“2.1.2”項制得紅參、黨參、熟地黃、山楂的供試品溶液和FEJ供試品溶液,按“2.1.1”項下色譜條件分別進(jìn)樣。通過查閱各味中藥成分文獻(xiàn)報道及各成分對照品色譜和光譜圖,對FEJ特征圖譜色譜峰進(jìn)行歸屬,見圖4。結(jié)果顯示,1號峰咖啡酸為黨參[17]、熟地黃[18-19]中的共有成分;2號峰綠原酸為紅參[20-21]、山楂[22]中的共有成分;3號峰洋地黃葉苷C[23]為熟地黃中的成分;4號峰阿魏酸為黨參、熟地黃中的共有成分;5號峰焦地黃苯乙醇苷A1、6號峰毛蕊花糖苷、7號峰焦地黃苯乙醇苷B1、8號峰異毛蕊花糖苷為熟地黃中的成分[23-24]。

    1-咖啡酸 2-綠原酸 3-洋地黃葉苷C 4-阿魏酸 5-焦地黃苯乙醇苷A1 (S) 6-毛蕊花糖苷 7-焦地黃苯乙醇苷B1 8-異毛蕊花糖苷

    供試品溶液制備時使用乙腈沉淀了蛋白,因而在此色譜條件所得的HPLC特征圖譜中暫無共有峰來自阿膠。

    2.2 化學(xué)模式分析

    2.2.1 HCA 采用Origin 2019軟件進(jìn)行分析,運用組間連接法,Euclidean距離為測度,以FEJ特征圖譜中的8個特征峰峰面積為觀測量,繪制聚類分析圖,見圖5。16批樣品分成3大類:S1~S4為第1類,S5~S11為第2類、S12~S16為第3類,說明不同批次的合格藥材所生產(chǎn)的產(chǎn)品含量存在一定差異,此可能與飲片的原藥材來源、飲片的加工等過程有關(guān)。

    2.2.2 PCA和PLS-DA 為了進(jìn)一步分析造成不同批次樣品之間差異的化學(xué)成分,以樣品S1~S16為變量,8種成分的峰面積為變量,利用SIMCA 14.1軟件進(jìn)行PCA和PLS-DA處理,結(jié)果見圖6。將16批樣品分為3類,以PLS-DA模型中重要性投影(VIP)值>1.0篩選不同批次樣品間的差異性成分,VIP值的大小表明了對樣品分類的貢獻(xiàn)程度。結(jié)果顯示,PCA和PLS-DA共得到3個成分,分別為8號峰異毛蕊花糖苷、5號峰焦地黃苯乙醇苷A1和3號峰洋地黃葉苷C,表明這3種成分是影響不同批次樣品之間質(zhì)量差異貢獻(xiàn)較大的特征成分。提示在FEJ整個生產(chǎn)過程中可重點監(jiān)控這3個成分,減小批次間的差異性,為FEJ Q-Marker的選定提供參考。

    2.3 分子對接驗證

    2.3.1 基于GEO數(shù)據(jù)庫的貧血癥的靶點篩選 《中國藥典》2020年版中,F(xiàn)EJ功能主治為“補氣養(yǎng)血”,常用于治療氣血兩虛、白細(xì)胞減少癥和貧血等癥?!皻庋獌商摗痹趥鹘y(tǒng)中醫(yī)理論中涉及多個臟腑的氣血陰陽失調(diào),在西醫(yī)中并無對應(yīng)的病癥,但從臨床表現(xiàn)來看,以貧血癥和慢性疲勞綜合征的臨床癥狀與其類似[25]。因此對貧血癥進(jìn)行下一步研究。

    在GEO數(shù)據(jù)庫中進(jìn)行基因表達(dá)譜芯片篩選,以“anemia normal”為關(guān)鍵詞、實驗以人源作為樣本,得到基因表達(dá)譜芯片GSE32719。利用R軟件及aggy、limma、pheatmap等R軟件包,對芯片進(jìn)行表達(dá)值背景矯正和表達(dá)譜數(shù)據(jù)的歸一化預(yù)處理,包括原始數(shù)據(jù)格式的轉(zhuǎn)換,缺失值補充,背景矯正及用數(shù)據(jù)標(biāo)準(zhǔn)化;以|logFC|≥1且<0.05為篩選條件,選擇統(tǒng)計學(xué)檢驗后差異最為顯著的基因作為研究對象,獲取貧血癥正常組與實驗組差異表達(dá)基因,并將探針名轉(zhuǎn)化為標(biāo)準(zhǔn)基因名。

    圖4 紅參、黨參、熟地黃、山楂與FEJ供試品溶液的HPLC圖

    圖5 16批FEJ的HCA圖

    結(jié)果得到貧血癥相關(guān)658個差異基因,其中,以logFC呈現(xiàn)正值為上調(diào)基因463個,列出調(diào)控最顯著的前20位差異基因;以logFC呈現(xiàn)負(fù)值為下調(diào)基因有195個,列出調(diào)控最顯著的前20位差異基因并繪制差異基因火山圖和熱圖,見圖7。

    圖6 16批FEJ的PLS-DA模型及VIP值

    2.3.2 蛋白質(zhì)-蛋白質(zhì)相互作用(protein-protein interaction,PPI)網(wǎng)絡(luò)的構(gòu)建及拓?fù)浞治?將上述得到的差異基因?qū)隨tring平臺(https://string-db.org),選擇“multiple proteins”模式,蛋白種屬設(shè)為“Home sapiens”,篩選條件為最低相互作用閾值取“highest confidence(≥0.9)”構(gòu)建PPI網(wǎng)絡(luò),見圖8。網(wǎng)絡(luò)中節(jié)點(node)表示靶點蛋白,邊(edge)表示蛋白與蛋白之間的相互作用,共涉及536個節(jié)點和232條邊。平均節(jié)點度值為0.866 6,平均局部聚類系數(shù)為0.225。PPI網(wǎng)絡(luò)富集的值為0.007 47,表明該網(wǎng)絡(luò)具有顯著的相互作用。

    圖7 貧血癥差異表達(dá)基因火山圖(A)和熱圖(B)

    圖8 PPI網(wǎng)絡(luò)圖

    將PPI數(shù)據(jù)結(jié)果以*.tsv格式導(dǎo)出,導(dǎo)入Cytoscape 3.7.2軟件,利用CytoNCA 2.1.6插件進(jìn)行網(wǎng)絡(luò)拓?fù)浞治?,計算靶點拓?fù)鋮?shù)包括網(wǎng)絡(luò)節(jié)點度(degree centrality,DC)、介數(shù)中心度(betweenness centrality,BC)和緊密度(closeness centrality,CC),運用-score函數(shù)[-score=[-mean()]/std(),其中,代表拓?fù)鋮?shù)值,mean()為平均數(shù),std()為標(biāo)準(zhǔn)差]進(jìn)行標(biāo)準(zhǔn)化處理,將標(biāo)準(zhǔn)化處理后的拓?fù)鋮?shù)值相加。度、緊密度和介數(shù)可以用于評價節(jié)點在網(wǎng)絡(luò)中的中心性和重要性,節(jié)點的度、緊密度和介數(shù)越高,節(jié)點在交互網(wǎng)絡(luò)中的重要性越高[26]。

    選擇排名前20的靶點作為貧血癥的關(guān)鍵靶點,各靶點特征參數(shù)見表1。

    2.3.3 關(guān)鍵靶點(靶蛋白)信息的收集與前處理 選取貧血癥差異基因中排名前10的關(guān)鍵靶點,在PDB數(shù)據(jù)庫(http://www. rcsb.org)中搜索并下載靶蛋白的三維結(jié)構(gòu),選擇來源為人源(human)、解析度?。ǎ?.3 nm)、X-Ray解析的晶體結(jié)構(gòu)作為研究對象。將靶蛋白導(dǎo)入Discovery Studio 2016 v16.1軟件中,利用Clean Protein模塊和Prepare Protein模塊對靶蛋白進(jìn)行刪除多余蛋白質(zhì)構(gòu)象、刪去水分子、補全不完整殘基、加氫等蛋白預(yù)處理。選擇靶蛋白特有的小分子配體作為活性中心,刪除位于活性口袋位置的配體小分子(原配體),暴露活性口袋并將其定義為對接體系中受體分子。

    表1 貧血癥靶點特征參數(shù)

    2.3.4 原配體對接驗證 將靶蛋白中的原配體單獨提取,基于靶蛋白中的原配體的活性口袋位點,利用CDocker模塊,將原配體對接到對應(yīng)的活性口袋中,計算靶蛋白中原配體的分子構(gòu)象與對接后的分子構(gòu)象的均方根偏差(RMSD)。若RMSD<0.20 nm,則說明對接所得分子構(gòu)象可以在一定程度上還原配體與受體的結(jié)合模式,從而證明所選擇的對接方法與參數(shù)設(shè)置具有合理性,結(jié)果見表2。得到靶蛋白KNG1、ISG15、KDM6A、DNM3、F2RL1、CYSLTR1,這些靶蛋白的RMSD均小于0.2 nm,說明所選擇的對接方法及參數(shù)設(shè)置合理,可以用于下一步與成分對接。

    2.3.5 成分信息的收集與前處理 復(fù)方制劑藥味組成繁雜,化學(xué)成分眾多,某些成分在特征圖譜的單一檢測波長下吸收較弱,難以體現(xiàn)。因此,通過查閱《中國藥典》和相關(guān)文獻(xiàn),對紅參、熟地黃、黨參、山楂中可能成為Q-Marker的成分進(jìn)行收集篩選,擴充FEJ Q-Marker的選擇。篩選得到紅參、熟地黃、黨參、山楂中的關(guān)鍵成分人參皂苷Rg1、人參皂苷Re、人參皂苷Rb1、黨參炔苷、枸櫞酸、毛蕊花糖苷。同時結(jié)合特征圖譜研究中得到的8個特征成分,最終整理得到13個成分進(jìn)行分子對接驗證。

    表2 靶點活性口袋信息及計算結(jié)果

    在PubChem數(shù)據(jù)庫中搜索并下載成分的2D結(jié)構(gòu),保存為.sdf格式。將結(jié)構(gòu)導(dǎo)入Discovery Studio 2016 v16.1軟件中,利用Prepare Ligands模塊對小分子配體進(jìn)行加氫、產(chǎn)生三維結(jié)構(gòu)和異構(gòu)體等預(yù)處理,并利用Minimize Ligands模塊賦予CHARMm力場進(jìn)行能量最小化。

    2.3.6 成分與關(guān)鍵靶點(靶蛋白)對接 利用Discovery Studio 2016 v16.1軟件的CDocker模塊,將潛在活性成分對接到處理好的活性口袋中。SDS Site Sphere半徑根據(jù)原配體對接情況設(shè)定,Top Hit設(shè)為10,Pose Cluster Radius設(shè)為0.5,其他參數(shù)設(shè)置為默認(rèn)值,對接打分值選取-cdocker interaction energy值。對接打分值評價成分與靶點的相互作用,CDocker Energy值越低(即-cdocker interaction energy值越高),則說明化學(xué)成分與蛋白受體對接體系越穩(wěn)定,對接結(jié)果越可靠[27],其成為活性成分可能性越高。以原配體與靶蛋白對接后的對接打分值作為閾值,選取對接打分值大于閾值的成分作為FEJ治療貧血癥的活性成分,最終得到8個對接打分值大于閾值的成分,結(jié)果見表3。

    分子對接結(jié)果表明,異毛蕊花糖苷、黨參炔苷、綠原酸、毛蕊花糖苷等成分與各個靶點均有較好的結(jié)合,對接打分值較高。其中,異毛蕊花糖苷與靶蛋白DNM3和6RZ5的對接打分值最高;黨參炔苷與靶蛋白KDM6A的對接打分值最高;綠原酸與靶蛋白F2RL1對接打分值最高;毛蕊花糖苷與靶蛋白ISG15對接打分值最高。說明異毛蕊花糖苷、黨參炔苷、綠原酸、毛蕊花糖苷4個成分對貧血癥靶點能夠較好地作用于貧血癥靶點,可作為FEJ治療貧血癥的活性成分。

    表3 活性成分的對接打分值

    進(jìn)一步分析活性成分與靶蛋白的相互作用,以異毛蕊花糖苷與靶蛋白DNM3為例分析成分與靶蛋白之間的相互作用,見圖9。圖中左側(cè)為成分與靶蛋白對接圖,中間為配體與受體結(jié)合處相互作用3D圖,右側(cè)為配體與受體結(jié)合處相互作用2D圖。

    異毛蕊花糖苷結(jié)合于DNM3的蛋白RBD空腔中,與DNM3的關(guān)鍵氨基酸SER41、ASP208、ASN236、SER238、LYS206、GLY43、SER45、LYS44、ARG237發(fā)生氫鍵作用,與LEU209、VAL64發(fā)生疏水作用。結(jié)合的關(guān)鍵氨基酸殘基有酸性氨基酸殘基ASP208和堿性氨基酸殘基LYS44、ARG237,可形成氫鍵或鹽橋,構(gòu)成一個具有可成藥性的潛在活性位點,從而發(fā)生相互作用。

    圖9 異毛蕊花糖苷與靶蛋白DNM3的相互作用

    通過分子對接技術(shù),對質(zhì)量控制中的特征成分和關(guān)鍵成分進(jìn)行活性驗證,得到異毛蕊花糖苷、黨參炔苷、綠原酸、毛蕊花糖苷4個成分可作為FEJ的Q-Marker。其中,綠原酸在山楂、熟地等多個藥味中均有,專屬性差。毛蕊花糖苷在FEJ中含量很低,含量測定難度高。異毛蕊花糖苷是熟地黃的專屬性成分,黨參炔苷是黨參的專屬性成分。這兩個成分可溯性和有效性良好。因此將Q-Marker異毛蕊花糖苷和黨參炔苷作為指標(biāo)成分,進(jìn)行FEJ含量測定研究。

    2.4 FEJ多成分含量測定

    2.4.1 色譜條件 Thermo hypersil gold C18色譜柱(250 mm×4.6 mm,5 μm),以0.1%甲酸乙腈溶液為流動相A,以0.1%甲酸水溶液為流動相B梯度洗脫,梯度洗脫程序:0~15 min,2%~20% A;15~30 min,20% A;30~35 min,20%~22% A;35~40 min,22%~25% A;體積流量1 mL/min;柱溫30 ℃;檢測波長330、270 nm;進(jìn)樣量10 μL。

    2.4.2 供試品溶液的制備 精密量取FEJ 2 mL,加入乙腈6 mL,渦旋至充分混勻,離心,取上清液,殘渣加入乙腈6 mL,渦旋至充分混勻,離心,合并上清液,蒸干,殘渣加30%甲醇溶解,并轉(zhuǎn)移至1 mL量瓶中,加30%甲醇至刻度,搖勻,即得供試品溶液。同法制得熟地黃陰性樣品溶液和黨參陰性樣品溶液。

    2.4.3 對照品溶液的制備 取異毛蕊花糖苷、黨參炔苷對照品適量,精密稱定,加甲醇制成含異毛蕊花糖苷25 μg/mL、黨參炔苷50 μg/mL的混合溶液,即得對照品溶液。

    2.4.4 線性關(guān)系考察 取異毛蕊花糖苷對照品適量,精密稱定,置于5 mL量瓶中,加甲醇制成含有異毛蕊花糖苷0.484 mg/mL的對照品溶液。取黨參炔苷對照品適量,精密稱定,置于5 mL量瓶中,加甲醇制成含有黨參炔苷0.478 mg/mL的對照品溶液。分別精密量取異毛蕊花糖苷對照品溶液和黨參炔苷對照品溶液2.5 mL,配制成混合對照品溶液。取混合對照品溶液,用甲醇分別稀釋0、2、5、10、50倍,配制成系列質(zhì)量濃度的混合對照品溶液,分別進(jìn)樣。以質(zhì)量濃度()對色譜峰面積()繪制標(biāo)準(zhǔn)曲線,進(jìn)行線性回歸,得2種成分回歸方程分別為異毛蕊花糖苷=14 533-1 357.8,=0.999 6,線性范圍4.84~242.00 μg/mL;黨參炔苷=10 026+1 3224,=0.999 6,線性范圍4.78~239.00 μg/mL。

    2.4.5 專屬性試驗 對異毛蕊花糖苷和黨參炔苷的混合對照品溶液、FEJ供試品溶液(批號1905075),以及批號1905075 FEJ所對應(yīng)的熟地黃陰性樣品和黨參陰性樣品進(jìn)行測定。

    通過雙波長檢測,在330 nm下檢測異毛蕊花糖苷,在270 nm下檢測黨參炔苷。結(jié)果表明,檢測波長330 nm下熟地黃陰性樣品對FEJ供試品中異毛蕊花糖苷的測定無干擾;檢測波長270 nm下黨參陰性樣品對FEJ供試品中黨參炔苷的測定無干擾,見圖10,說明該方法專屬性強。

    2.4.6 精密度試驗 取FEJ樣品(批號1905075),按“2.4.2”項下方法制備供試品溶液,連續(xù)進(jìn)樣6次,測定異毛蕊花糖苷和黨參炔苷,結(jié)果各成分峰面積的RSD分別為0.86%、0.87%,滿足精密度限度(RSD)為6%[28],表明儀器精密度良好。

    圖10 FEJ供試品(a)、混合對照品(b)、陰性對照溶液(c) 在檢測波長330 nm (A) 和270 nm (B) 下的HPLC圖

    2.4.7 重復(fù)性試驗 取FEJ樣品(批號1905075),按“2.4.2”項下方法制備6份供試品溶液,測定異毛蕊花糖苷和黨參炔苷,結(jié)果各成分峰面積的RSD分別為2.01%、2.01%,滿足重復(fù)性限度(RSD)為6%[28],表明方法重復(fù)性良好。

    2.4.8 穩(wěn)定性試驗 取FEJ樣品(批號1905075),按“2.4.2”項下方法制備供試品溶液,分別于0、2、4、8、12、24 h進(jìn)樣,測定異毛蕊花糖苷和黨參炔苷,結(jié)果各成分峰面積的RSD分別為1.68%、1.67%,表明24 h內(nèi),供試品溶液穩(wěn)定性良好。

    2.4.9 加樣回收試驗 取FEJ樣品(批號1905075)6份,每份1 mL,按照接近1∶1的量計算所要添加的2種成分的對照品,按“2.4.2”項下方法制備供試品溶液,依法提取測定,計算各成分的回收率。結(jié)果表明,異毛蕊花糖苷和黨參炔苷的平均回收率分別為101.13%、103.19%,RSD分別為2.71%、2.90%,滿足回收率限度為80%~115%[28],表明方法的準(zhǔn)確度良好。

    2.4.10 樣品測定結(jié)果 取16個批次(S1~S16)的FEJ,每個批次平行制備3份供試品溶液,依次進(jìn)樣測定,結(jié)果見表5。16批FEJ中異毛蕊花糖苷和黨參炔苷的質(zhì)量濃度分別在6.05~12.30 μg/mL、9.14~16.30 μg/mL。

    表4 FEJ中異毛蕊花糖苷和黨參炔苷含量測定結(jié)果(, n = 3)

    3 討論

    3.1 FEJ Q-Marker的選擇及分析

    隨著中醫(yī)藥學(xué)的飛速發(fā)展,中成藥的普及性和應(yīng)用性越來越廣。2016年劉昌孝院士[7-8]提出了中藥Q-Marker概念,形成了中藥產(chǎn)品質(zhì)量控制的新理念。從生源途徑、化學(xué)物質(zhì)基礎(chǔ)、炮制加工過程、藥效研究和中醫(yī)理論等多個層面,密切關(guān)聯(lián)了中藥的物質(zhì)基礎(chǔ)、有效性和質(zhì)量控制,為全面提升中藥質(zhì)量標(biāo)準(zhǔn)提供理論和方法的支持[29-30]。

    本研究從特征圖譜中指認(rèn)了8個成分,保證Q-Marker的特有性和可溯性?;瘜W(xué)計量學(xué)分析提示了異毛蕊花糖苷、洋地黃葉苷C和焦地黃苯乙醇苷A1這3個成分應(yīng)被重點關(guān)注,為Q-Marker篩選提供參考。本研究通過查閱《中國藥典》和相關(guān)文獻(xiàn),對紅參、熟地黃、黨參、山楂中可能成為Q-Marker的成分進(jìn)行收集篩選,擴充FEJ Q-Marker的選擇。紅參和熟地黃是組方中的臣藥。其中《中國藥典》2020年版中使用人參皂苷Rb1、Re、Rg1作為紅參指標(biāo)成分,將這些指標(biāo)成分納入到后續(xù)研究中。黨參和山楂是組方中的佐使藥。黨參能夠提高造血機能、抗疲勞、增強機體免疫力等,其主要活性成分黨參炔苷[31-33]?!吨袊幍洹?020年版收載了黨參以黨參炔苷為對照的薄層鑒別方法,專屬性好,因此選定黨參炔苷進(jìn)一步研究。同時將山楂含量測定的指標(biāo)成分枸櫞酸納入后續(xù)研究中。

    分子對接技術(shù)作為一種快速的活性驗證方法,已被廣泛應(yīng)用于藥材飲片、中藥復(fù)方等活性成分的篩選和驗證中[34-35]。分子對接結(jié)果顯示,異毛蕊花糖苷、黨參炔苷、綠原酸、毛蕊花糖苷4個成分是FEJ治療貧血癥的活性成分,可作為Q-Marker。

    結(jié)合化合物來源的專屬性以及含量可測性等因素,選定異毛蕊花糖苷和黨參炔苷為FEJ的Q-Marker。這2個成分傳遞和溯源明確、具有活性且專屬性強,整體滿足中藥Q-Marker的五原則特征。

    3.2 提取條件考察

    FEJ的制備工藝采用紅參、黨參、熟地黃、山楂水煎煮的提取方法,化學(xué)成分組成復(fù)雜且多為大極性化合物。實驗前期,在特征圖譜和含量測定的供試品制備中,考察了不同提取溶劑,包括甲醇、乙醇、醋酸乙酯-正丁醇、乙腈對FEJ進(jìn)行提取,以及不同溶解定容溶劑,包括水及30%、70%、100%甲醇對FEJ上清液殘渣進(jìn)行溶解。研究表明使用甲醇作為提取溶劑時,影響部分成分的提取,且基線不穩(wěn)定;使用乙醇提取,色譜峰多,分離度差;使用醋酸乙酯-正丁醇(1∶1)和乙腈提取后,色譜峰分離度較好。由正丁醇沸點高,回收難度,因此采用乙腈提取的方法進(jìn)行制備。此外,研究發(fā)現(xiàn)不同溶解溶劑對色譜峰影響差異不明顯,其中30%甲醇復(fù)溶的色譜峰分離度較好,因此選用30%甲醇作為復(fù)溶溶劑。同時,優(yōu)化了提取溶劑的體積和提取次數(shù),保證含量測定中FEJ供試品能夠被提取完全。

    3.3 色譜條件考察

    本研究對特征圖譜和含量測定中的檢測波長、流動相組成等色譜條件進(jìn)行優(yōu)化。在特征圖譜研究中,利用二極管陣列檢測器在190~400 nm在線掃描,記錄并比較了不同波長下的特征圖譜。在203 nm下檢測到人參皂苷Rg1,然而由于響應(yīng)值低(含量少,推測可能與復(fù)方中紅參用量少有關(guān)),基線噪音影響較大,因此未選用203 nm作為檢測波長。通過查閱文獻(xiàn)以及DAD掃描顯示,黨參炔苷的最大吸收波長在270 nm左右[38],在檢測波長270 nm下,色譜圖顯示前15 min受干擾嚴(yán)重,苯乙醇苷類成分反映不完整,因此未選用270 nm作為檢測波長。在含量測定研究方法中,330 nm下色譜峰響應(yīng)值和各峰峰高比例較好,色譜峰數(shù)目較多且分離度較好。黨參炔苷在330 nm下紫外吸收極弱,因此,未在特征圖譜中對其進(jìn)行指認(rèn)。綜上,最終確定330 nm為檢測波長。

    本研究考察了5種流動相體系,包括乙腈-水、乙腈-0.05%甲酸水、乙腈-0.1%甲酸水、0.05%甲酸乙腈-0.05%甲酸水、0.1%甲酸乙腈-0.1%甲酸水。結(jié)果表明,當(dāng)0.1%甲酸乙腈-0.1%甲酸水作為流動相體系時,共有峰峰型良好,分離度好,且以甲酸水溶液為水相時可以進(jìn)一步改善各成分峰的峰形,效果最佳。此外,本實驗也考察了不同柱溫(25、30、35 ℃),不同體積流量(0.6、0.8、1.0 mL/min),結(jié)果發(fā)現(xiàn),當(dāng)柱溫為30 ℃、體積流量為1.0 mL/min,檢測波長為330 nm時,所得色譜峰的數(shù)目較多、峰形及分離度等方面均表現(xiàn)良好。

    在含量測定研究中,利用二極管陣列檢測器在190~400 nm在線掃描,得到異毛蕊花糖苷對照品在327 nm左右波長下有最大吸收;黨參炔苷對照品在268 nm左右波長下有最大吸收。最終選定在檢測波長330 nm下測定異毛蕊花糖苷,在檢測波長270 nm下測定黨參炔苷。對FEJ特征圖譜的洗脫條件進(jìn)行優(yōu)化,在保證成分達(dá)到分離效果的同時,縮短檢測時間。在該流動相洗脫條件下,異毛蕊花糖苷和黨參炔苷出峰情況良好,分離度達(dá)到要求。

    本研究提供了一種FEJ中Q-Marker研究的創(chuàng)新思路與方法,提升了FEJ質(zhì)量控制水平。

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

    [1] 中國藥典 [S]. 一部. 2020: 1322-1323.

    [2] 金永志, 李夢帆, 黃榮, 等. 當(dāng)歸補血湯促進(jìn)血管新生作用的研究進(jìn)展 [J]. 中醫(yī)藥導(dǎo)報, 2019, 25(13): 109-111,117.

    [3] 陳榮榮, 郭浩, 徐硯通, 等. 中藥復(fù)方和有效成分對血管新生促進(jìn)或抑制作用的研究進(jìn)展 [J]. 中草藥, 2013, 44(23): 3413-3421.

    [4] 史永平, 劉發(fā)生, 孔浩天, 等. 利用斑馬魚模型研究丹紅注射液及代表性成分的促進(jìn)血管生成作用 [J]. 時珍國醫(yī)國藥, 2019, 30(11): 2588-2590.

    [5] 王成, 龔莉虹, 郭朝成, 等. 模式生物斑馬魚在中藥藥效物質(zhì)篩選中的應(yīng)用進(jìn)展 [J]. 中草藥, 2019, 50(24): 6125-6134.

    [6] 陳維武, 白瑪卓瑪, 郭勝亞, 等. 二十五味珍珠片對斑馬魚微血管功能的影響 [J]. 中國實驗方劑學(xué)雜志, 2019, 25(13): 49-54.

    [7] 劉昌孝, 陳士林, 肖小河, 等. 中藥質(zhì)量標(biāo)志物(QMarker): 中藥產(chǎn)品質(zhì)量控制的新概念 [J]. 中草藥, 2016, 47(9): 1443-1457.

    [8] Liu C X, Cheng Y Y, Guo D A,. A new concept on quality marker for quality assessment and process control of Chinese medicines [J]., 2017, 9(1): 3-13.

    [9] 張鐵軍, 白鋼, 陳常青, 等. 基于“五原則”的復(fù)方中藥質(zhì)量標(biāo)志物(Q-marker)研究路徑 [J]. 中草藥, 2018, 49(1): 1-13.

    [10] 彭平, 熊少哲, 張蓓, 等. 基于全方特征圖譜質(zhì)量表征的芩連制劑質(zhì)量評價研究 [J]. 中草藥, 2021, 52(5): 1312-1322.

    [11] 許春芳, 崔莉莉, 鄢紅霞, 等. 昆仙膠囊HPLC特征指紋圖譜和多成分定量測定研究 [J]. 中草藥, 2020, 51(24): 6213-6220.

    [12] 張夢晨, 謝輝, 陸兔林, 等. 基于指紋圖譜、化學(xué)計量學(xué)、網(wǎng)絡(luò)藥理學(xué)的半夏湯洗前后質(zhì)量評價 [J]. 中草藥, 2021, 52(10): 2897-2908.

    [13] 高森, 王蘋, 唐鋮, 等. 基于HPLC指紋圖譜、多指標(biāo)成分含量測定及化學(xué)計量學(xué)的濕熱痹片質(zhì)量評價 [J]. 中草藥, 2020, 51(21): 5454-5461.

    [14] 陳曉璐, 郭振旺, 鄧家剛, 等. 基于肉桂質(zhì)量標(biāo)志物(Q-Marker)預(yù)測的質(zhì)量控制研究 [J]. 中草藥, 2021, 52(9): 2707-2718.

    [15] 李志立, 楊冰, 翟艷敏, 等. 基于分子對接技術(shù)虛擬篩選紅花中抗心肌缺血活性組分 [J]. 中國中藥雜志, 2020, 45(12): 2881-2890.

    [16] 王玉, 楊雪, 靳曉杰, 等. 基于中醫(yī)藥整合藥理學(xué)平臺、GEO數(shù)據(jù)庫芯片及分子對接探討大黃抗肝癌的作用機制 [J]. 中草藥, 2020, 51(20): 5207-5219.

    [17] 凈晶, 張吉發(fā), 李慧, 等. 長花黨參化學(xué)成分的研究 [J]. 時珍國醫(yī)國藥, 2013, 24(10): 2340-2342.

    [18] 馮蓬, 甘東浩, 王衛(wèi)國, 等. 基于網(wǎng)絡(luò)藥理學(xué)的熟地黃干預(yù)骨質(zhì)疏松癥分子機制研究 [J]. 山東中醫(yī)雜志, 2020, 39(2): 177-184.

    [19] Tao J H, Zhao M, Ling Y,. UPLC-Q-TOF/MS-based metabolic profiles of bioactive components inandherb pair by rat intestinal bacteria [J]., 2017, 9(2): 147-152.

    [20] Ali M B, Hahn E J, Paek K Y. CO2-induced total phenolics in suspension cultures ofC. A. Mayer roots: Role of antioxidants and enzymes [J]., 2005, 43(5): 449-457.

    [21] Chung I M, Lim J J, Ahn M S,. Comparative phenolic compound profiles and antioxidative activity of the fruit, leaves, and roots of Korean ginseng (Meyer) according to cultivation years [J]., 2016, 40(1): 68-75.

    [22] 趙博, 田恒康, 孫道涵, 等. 基于酚類特征圖譜的山楂質(zhì)量表征關(guān)聯(lián)分析研究 [J]. 北京中醫(yī)藥大學(xué)學(xué)報, 2017, 40(8): 684-692.

    [23] 鄒獻(xiàn)亮, 陳颋, 華臘, 等. 一測多評法同時測定熟地黃中4種苯乙醇苷 [J]. 中成藥, 2019, 41(5): 1085-1090.

    [24] Shen L J, Chen H L, Zhu Q F,. Identification of bioactive ingredients with immuno-enhancement and anti-oxidative effects from Fufang-Ejiao-Syrup by LC-MSncombined with bioassays [J]., 2016, 117: 363-371.

    [25] 張萬年. 益氣安神膏治療氣血兩虛型慢性疲勞綜合征的臨床研究 [D]. 廣州:廣州中醫(yī)藥大學(xué), 2018.

    [26] 馬婧, 陳茜, 邊亞倩, 等. 基于系統(tǒng)中藥學(xué)的丹參活血化瘀功效標(biāo)志物研究 [J]. 中國中藥雜志, 2020, 45(14): 3259-3265.

    [27] 劉美鳳, 蔣利榮, 劉華鼐, 等. 番石榴葉抗Ⅱ型糖尿病活性成分的虛擬篩選 [J]. 華南理工大學(xué)學(xué)報: 自然科學(xué)版, 2011, 39(3): 28-31.

    [28] 中國藥典 [S]. 四部. 2020: 480-483.

    [29] 彭任, 陸兔林, 胡立宏, 等. 中藥飲片質(zhì)量標(biāo)志物(Q-marker)研究進(jìn)展 [J]. 中草藥, 2020, 51(10): 2603-2610.

    [30] 陳麗華, 肖發(fā)林, 黃詩雨, 等. 中藥質(zhì)量評價研究思路及創(chuàng)新發(fā)展趨勢 [J]. 中草藥, 2021, 52(9): 2541-2547.

    [31] 廖江敏, 李林玉, 瞿云安, 等. 黨參藥材活性成分黨參炔苷的化學(xué)穩(wěn)定性研究 [J]. 天然產(chǎn)物研究與開發(fā), 2020, 32(3): 385-388.

    [32] Zhang P, Hu L H, Bai R B,. Structural characterization of a pectic polysaccharide fromand its immunomodulatory activitiesand[J]., 2017, 104(Pt A): 1359-1369.

    [33] Zhao X N, Hu Y L, Wang D Y,. The comparison of immune-enhancing activity of sulfated polysaccharidses fromand[J]., 2013, 98(1): 438-443.

    [34] 龔曼, 楊聯(lián)河, 朱利利, 等. 鄂西香茶菜甲素作為抗腫瘤先導(dǎo)物的發(fā)現(xiàn)與抗腫瘤活性驗證 [J]. 中國中藥雜志, 2021, 46(16): 4061-4068.

    [35] 郭錦晨, 王莖, 孫宇潔, 等. 基于網(wǎng)絡(luò)藥理學(xué)及分子對接探討桑葉-菊花治療高血壓的作用機制 [J]. 天然產(chǎn)物研究與開發(fā), 2021, 33(5): 847-858.

    [36] 鄭昆, 倫志彩, 羅兆亮. 超高效液相色譜-波長切換法同時測定健脾丸(濃縮丸)中4種成分的含量 [J]. 食品與藥品, 2021, 23(2): 142-146.

    Study of quality markers of Fufang E’jiao Jiang based on characteristic spectrum, chemometrics and molecular docking

    XU Xiao1, ZHANG Yan2, REN Xue-yang1, WANG Yu1, DONG Ying1, SONG Ruo-lan1, YU A-xiang1, WEI Jing1, MA Jia-mu1, SHE Gai-mei1

    1. School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing 100102, China 2. National Engineering Research Center for Gelatin-based Traditional Chinese Medicine, Dong-E-E-Jiao Co., Ltd., Liaocheng 252201, China

    Guided by the “Five Principles” of quality markers (Q-Marker), combined with characteristic spectrum, chemometrics and molecular docking and content determination, developed Q-Marker for the treatment of anemia with Fufang E'jiao Jiang (復(fù)方阿膠漿, FEJ).The characteristic spectrum of FEJ was determined by HPLC using Agilent Zorbax SB AQ-C18column. The mobile phase was 0.1% formic acid acetonitrile-0.1% formic acid water, and the flow rate was 1.0 mL/min, detection wavelength was 330 nm. Sixteen batches of FEJ characteristic spectrum were established with jionoside A1as the reference peak. Based on the results of common peak area determination of characteristic spectrums, the quality of different batches of FEJ was evaluated by hierarchica cluster analysis (HCA), principal components analysis (PCA) and partial least squares-discrimination analysis (PLS-DA). Six key components were collected from each herb inand literatures. The differential genes of anemia was searched in the GEO database, and Discovery Studio software was used to verify the components in FEJ. The content of two components in 16 batches was determined at the same time at 330 nm and 270 nm.There were eight components were identified, namely caffeic acid, chlorogenic acid, purpureaside C, ferulic acid, jionoside A1, acteoside, jionoside B1, isoacteoside. The results of HCA, PCA and PLS-DA showed that 15 batches of FEJ can be divided into three categories. Isoacteoside, purpureaside C, and jionoside A1were contributing to the quality difference between different batches of FEJ. The results of molecular docking showed that isoacteoside, lobetyolin, chlorogenic acid, and acteoside can bind to six key targets such as KNG1, ISG15, KDM6A, DNM3, F2RL1, and CYSLTR1 through hydrogen bonds, hydrophobic bonds, π-π bonds, etc. They can be used as Q-Markers of FEJ. The contents of isoacteoside and lobetyolin in 16 batches of FEJ were 6.05—12.30 μg/mL and 9.14—16.30 μg/mL, respectively.Q-Markers in FEJ was explored by characteristic spectrum, chemometrics and molecular docking, which can be used for quality control and evaluation of FEJ.

    Fufang E’jiao Jiang; quality markers (Q-Markers); characteristic spectrum; chemometrics; molecular docking; anemia; hierarchical cluster analysis; principal components analysis; partial least squares-discrimination analysis; quality control; caffeic acid; chlorogenic acid; purpureaside C; ferulic acid; jionoside A1; acteoside; jionoside B1; isoacteoside

    R283.6

    A

    0253 - 2670(2021)23 - 7148 - 14

    10.7501/j.issn.0253-2670.2021.23.009

    2021-06-16

    國家重點研發(fā)計劃(2018YFC1707300);國家重點研發(fā)計劃項目(2018YFC1706300)

    許 嘯(1995—),女,碩士研究生,研究方向為中(民族)藥藥效成分和新藥創(chuàng)制研究。Tel: 13521701679 E-mail: xux_26@163.com

    折改梅(1976—),女,研究員,博士生導(dǎo)師,研究方向為中(民族)藥藥效成分和新藥創(chuàng)制研究。Tel: (010)53912129 E-mail: shegaimei@126.com

    [責(zé)任編輯 鄭禮勝]

    猜你喜歡
    糖苷黨參批號
    一種JTIDS 信號批號的離線合批方法
    Anticancer Properties of Lobetyolin, an Essential Component of Radix Codonopsis (Dangshen)
    醫(yī)學(xué)科技期刊中藥品生產(chǎn)批號標(biāo)注探析
    中藥材批號劃分與質(zhì)量管理
    中成藥(2018年7期)2018-08-04 06:04:26
    復(fù)方黨參口服液免疫增強作用的配伍相關(guān)性
    中成藥(2018年4期)2018-04-26 07:12:58
    黨參中二氧化硫殘留量的測定
    中成藥(2018年1期)2018-02-02 07:20:31
    HPLC法同時測定黨參固本丸中7種成分
    中成藥(2017年10期)2017-11-16 00:50:06
    甜葉菊及其糖苷的研究與發(fā)展探索
    食品界(2016年4期)2016-02-27 07:36:47
    利用烷基糖苷遷移和擴張共軛亞油酸囊泡pH窗口
    固體超強酸催化合成丁基糖苷
    av欧美777| 国产亚洲欧美精品永久| 国产精品亚洲av一区麻豆| 97超级碰碰碰精品色视频在线观看| 国产精品综合久久久久久久免费 | 岛国视频午夜一区免费看| 久久久水蜜桃国产精品网| 国产成人精品久久二区二区91| 午夜精品在线福利| www国产在线视频色| 精品国产乱子伦一区二区三区| 丝袜美足系列| 成人av一区二区三区在线看| 婷婷六月久久综合丁香| 91av网站免费观看| 99久久精品国产亚洲精品| 制服人妻中文乱码| 一级a爱视频在线免费观看| 伊人久久大香线蕉亚洲五| 国产又色又爽无遮挡免费看| 黄频高清免费视频| 亚洲免费av在线视频| 国产精品一区二区免费欧美| 一二三四在线观看免费中文在| 脱女人内裤的视频| 国产伦一二天堂av在线观看| av中文乱码字幕在线| 亚洲精品粉嫩美女一区| 别揉我奶头~嗯~啊~动态视频| 欧美国产日韩亚洲一区| 老司机在亚洲福利影院| 国产主播在线观看一区二区| 老鸭窝网址在线观看| 亚洲一区高清亚洲精品| 国产成+人综合+亚洲专区| 女人被躁到高潮嗷嗷叫费观| 日韩三级视频一区二区三区| 日韩有码中文字幕| 99国产精品99久久久久| 日韩欧美国产一区二区入口| 国产一区二区三区在线臀色熟女| 日韩有码中文字幕| 99久久99久久久精品蜜桃| 悠悠久久av| 亚洲成人国产一区在线观看| 久久性视频一级片| 我的亚洲天堂| 国产精品综合久久久久久久免费 | 精品欧美一区二区三区在线| 国产精品国产高清国产av| 波多野结衣av一区二区av| 女生性感内裤真人,穿戴方法视频| 男人的好看免费观看在线视频 | 久久这里只有精品19| 免费在线观看日本一区| 久久久久久人人人人人| 亚洲无线在线观看| av片东京热男人的天堂| 搞女人的毛片| 国产又色又爽无遮挡免费看| 免费看十八禁软件| 9色porny在线观看| 亚洲成人久久性| 在线免费观看的www视频| 一级,二级,三级黄色视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av成人一区二区三| 欧美av亚洲av综合av国产av| 高清黄色对白视频在线免费看| 亚洲人成伊人成综合网2020| 亚洲专区字幕在线| 亚洲国产看品久久| 欧美成人一区二区免费高清观看 | 18禁国产床啪视频网站| 如日韩欧美国产精品一区二区三区| 高清在线国产一区| 激情视频va一区二区三区| 久久精品国产清高在天天线| 后天国语完整版免费观看| 欧美久久黑人一区二区| 每晚都被弄得嗷嗷叫到高潮| 国产成人精品在线电影| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 国产av一区二区精品久久| 日本 欧美在线| 校园春色视频在线观看| 女同久久另类99精品国产91| 亚洲精品美女久久久久99蜜臀| 桃红色精品国产亚洲av| 久久国产亚洲av麻豆专区| 亚洲 国产 在线| www.999成人在线观看| 亚洲成av人片免费观看| 男男h啪啪无遮挡| 自线自在国产av| 一个人免费在线观看的高清视频| 亚洲成国产人片在线观看| 美国免费a级毛片| 久久人人精品亚洲av| 欧美日韩亚洲综合一区二区三区_| 日韩一卡2卡3卡4卡2021年| 在线国产一区二区在线| 国产精品久久久av美女十八| 中亚洲国语对白在线视频| 亚洲欧洲精品一区二区精品久久久| 国产精品免费视频内射| 国产精品一区二区三区四区久久 | 桃色一区二区三区在线观看| 国产精品99久久99久久久不卡| 日韩精品免费视频一区二区三区| 男人的好看免费观看在线视频 | 美女扒开内裤让男人捅视频| 亚洲一区二区三区色噜噜| 亚洲熟妇熟女久久| 欧美丝袜亚洲另类 | 两人在一起打扑克的视频| 搡老熟女国产l中国老女人| 国产在线观看jvid| 欧美成人一区二区免费高清观看 | 色综合站精品国产| 中文字幕色久视频| 在线观看免费视频日本深夜| 色综合站精品国产| 亚洲精品国产一区二区精华液| 一边摸一边抽搐一进一出视频| 在线观看一区二区三区| 制服诱惑二区| 美国免费a级毛片| 巨乳人妻的诱惑在线观看| 两个人视频免费观看高清| 国产精品久久电影中文字幕| 99久久久亚洲精品蜜臀av| 色综合站精品国产| 国产成人啪精品午夜网站| 麻豆成人av在线观看| 精品乱码久久久久久99久播| 操美女的视频在线观看| 精品福利观看| 国产精品免费一区二区三区在线| 午夜老司机福利片| 9热在线视频观看99| 日本 av在线| 日韩大码丰满熟妇| 久久久久久久久中文| 久久人人精品亚洲av| 免费人成视频x8x8入口观看| 色尼玛亚洲综合影院| 黄网站色视频无遮挡免费观看| 淫妇啪啪啪对白视频| 一进一出抽搐gif免费好疼| 国产一区二区激情短视频| 欧美日韩福利视频一区二区| 国产亚洲精品一区二区www| 欧美日本亚洲视频在线播放| cao死你这个sao货| 亚洲在线自拍视频| 女人高潮潮喷娇喘18禁视频| 日本免费a在线| 国产亚洲欧美精品永久| 妹子高潮喷水视频| 嫩草影院精品99| 操美女的视频在线观看| 欧美精品亚洲一区二区| www.999成人在线观看| 一级作爱视频免费观看| 亚洲一区二区三区不卡视频| 久久精品国产亚洲av香蕉五月| 一级,二级,三级黄色视频| 欧美大码av| 女人被躁到高潮嗷嗷叫费观| 国产伦人伦偷精品视频| 级片在线观看| 国产精品久久久久久亚洲av鲁大| 国产成人精品久久二区二区91| 校园春色视频在线观看| 久久国产精品影院| videosex国产| 欧美中文综合在线视频| 自线自在国产av| 久久中文字幕一级| 美女扒开内裤让男人捅视频| 黄色 视频免费看| 国产午夜精品久久久久久| 三级毛片av免费| 国产亚洲精品av在线| 久久国产乱子伦精品免费另类| 久久中文字幕人妻熟女| 国产人伦9x9x在线观看| 亚洲av成人不卡在线观看播放网| 最近最新中文字幕大全免费视频| 精品久久久久久久久久免费视频| 亚洲精品在线观看二区| 国产精品久久视频播放| 真人做人爱边吃奶动态| 丰满的人妻完整版| 日本一区二区免费在线视频| 激情在线观看视频在线高清| 久久久久久久午夜电影| 亚洲第一电影网av| 成人三级黄色视频| 一区在线观看完整版| av电影中文网址| 一进一出好大好爽视频| 成年女人毛片免费观看观看9| 亚洲欧美日韩另类电影网站| 俄罗斯特黄特色一大片| 免费观看精品视频网站| 精品一品国产午夜福利视频| 国产精品秋霞免费鲁丝片| 手机成人av网站| 日本精品一区二区三区蜜桃| 婷婷丁香在线五月| 亚洲av成人不卡在线观看播放网| 十八禁人妻一区二区| 在线国产一区二区在线| 亚洲成av片中文字幕在线观看| 国产精品美女特级片免费视频播放器 | 欧美 亚洲 国产 日韩一| 精品人妻在线不人妻| 久久久国产成人免费| 久久精品亚洲熟妇少妇任你| 黄片小视频在线播放| aaaaa片日本免费| av免费在线观看网站| 视频区欧美日本亚洲| 人妻久久中文字幕网| 免费在线观看影片大全网站| 老鸭窝网址在线观看| 久久 成人 亚洲| 91成人精品电影| 国产精华一区二区三区| 色老头精品视频在线观看| 88av欧美| 成人18禁在线播放| 欧美日本亚洲视频在线播放| 99在线视频只有这里精品首页| 熟妇人妻久久中文字幕3abv| 国产一区二区在线av高清观看| 成人三级做爰电影| 天天躁夜夜躁狠狠躁躁| 一级,二级,三级黄色视频| 亚洲av片天天在线观看| 老汉色av国产亚洲站长工具| 国产av一区在线观看免费| av天堂在线播放| 亚洲五月婷婷丁香| 欧美日韩精品网址| а√天堂www在线а√下载| 欧美最黄视频在线播放免费| 18禁国产床啪视频网站| 亚洲欧美精品综合久久99| 老汉色∧v一级毛片| 国产成人精品久久二区二区91| 可以在线观看毛片的网站| 午夜a级毛片| 美女 人体艺术 gogo| 亚洲熟女毛片儿| 中文字幕另类日韩欧美亚洲嫩草| 亚洲色图综合在线观看| 国产成年人精品一区二区| 久久久久国产一级毛片高清牌| 亚洲精品国产区一区二| 亚洲黑人精品在线| 两性午夜刺激爽爽歪歪视频在线观看 | 搞女人的毛片| 久久久精品欧美日韩精品| 伦理电影免费视频| 搡老岳熟女国产| 可以在线观看的亚洲视频| 国产精品99久久99久久久不卡| 精品久久久久久久毛片微露脸| 中文字幕最新亚洲高清| av欧美777| 日本vs欧美在线观看视频| 亚洲中文字幕日韩| 久久精品国产亚洲av高清一级| 日日爽夜夜爽网站| 亚洲男人天堂网一区| 亚洲avbb在线观看| 国产在线观看jvid| 在线av久久热| 久久久久亚洲av毛片大全| 久久午夜亚洲精品久久| 午夜两性在线视频| 视频在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看 | 成在线人永久免费视频| 一区二区三区国产精品乱码| 国产一区二区在线av高清观看| 黄色丝袜av网址大全| 韩国精品一区二区三区| 69精品国产乱码久久久| 身体一侧抽搐| 日韩大码丰满熟妇| 久久国产精品影院| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲美女久久久| www.999成人在线观看| 国内久久婷婷六月综合欲色啪| 亚洲成人精品中文字幕电影| 操美女的视频在线观看| 日韩欧美国产一区二区入口| 麻豆久久精品国产亚洲av| 久久人妻熟女aⅴ| 精品久久久久久久人妻蜜臀av | 精品欧美一区二区三区在线| 久久久国产欧美日韩av| 日本撒尿小便嘘嘘汇集6| 亚洲精品中文字幕在线视频| 成人手机av| 岛国在线观看网站| 国内毛片毛片毛片毛片毛片| 国产亚洲精品综合一区在线观看 | 免费av毛片视频| 亚洲第一欧美日韩一区二区三区| 色在线成人网| 97超级碰碰碰精品色视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩一区二区精品| 精品一区二区三区视频在线观看免费| 老鸭窝网址在线观看| 美国免费a级毛片| 成人免费观看视频高清| 成人精品一区二区免费| 女性生殖器流出的白浆| 国产99久久九九免费精品| 他把我摸到了高潮在线观看| 正在播放国产对白刺激| 1024香蕉在线观看| av天堂久久9| 久久精品国产亚洲av香蕉五月| 久久草成人影院| 午夜影院日韩av| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区色噜噜| 亚洲人成电影观看| 久久香蕉精品热| 国产精品 国内视频| 久久久久国产精品人妻aⅴ院| 俄罗斯特黄特色一大片| 久久精品亚洲精品国产色婷小说| 大陆偷拍与自拍| 免费一级毛片在线播放高清视频 | 国产亚洲精品av在线| 国产亚洲欧美在线一区二区| 啦啦啦韩国在线观看视频| 亚洲男人天堂网一区| 高清在线国产一区| av欧美777| 老司机福利观看| 免费看a级黄色片| 精品一区二区三区四区五区乱码| www.自偷自拍.com| 美女 人体艺术 gogo| 日本免费一区二区三区高清不卡 | 久久人人精品亚洲av| 国产麻豆69| 精品电影一区二区在线| 日韩有码中文字幕| 每晚都被弄得嗷嗷叫到高潮| 操出白浆在线播放| 18禁观看日本| 很黄的视频免费| 欧美+亚洲+日韩+国产| 亚洲一卡2卡3卡4卡5卡精品中文| 丝袜在线中文字幕| 国产精品美女特级片免费视频播放器 | 国产av精品麻豆| 亚洲色图综合在线观看| 色综合站精品国产| 韩国av一区二区三区四区| 久久久久久久午夜电影| 丝袜美足系列| av在线天堂中文字幕| 欧美日本视频| 日本免费a在线| 又黄又粗又硬又大视频| 搡老妇女老女人老熟妇| 欧美 亚洲 国产 日韩一| 国产亚洲精品久久久久5区| 91九色精品人成在线观看| 国产一区二区三区在线臀色熟女| 咕卡用的链子| 中文字幕另类日韩欧美亚洲嫩草| 超碰成人久久| 国产成人精品久久二区二区91| 一区二区三区国产精品乱码| 国产av精品麻豆| 91九色精品人成在线观看| av在线天堂中文字幕| 女人爽到高潮嗷嗷叫在线视频| 黄色视频,在线免费观看| 亚洲中文日韩欧美视频| 90打野战视频偷拍视频| 欧美乱码精品一区二区三区| 18禁美女被吸乳视频| 欧美性长视频在线观看| cao死你这个sao货| 侵犯人妻中文字幕一二三四区| 男女下面插进去视频免费观看| 不卡av一区二区三区| 一夜夜www| 欧美精品啪啪一区二区三区| 久久精品影院6| 亚洲国产精品久久男人天堂| 日韩欧美国产一区二区入口| 日韩成人在线观看一区二区三区| 中文字幕人成人乱码亚洲影| 国产精华一区二区三区| 咕卡用的链子| 欧美黑人欧美精品刺激| 成人国语在线视频| 亚洲一区二区三区色噜噜| 丝袜美腿诱惑在线| 正在播放国产对白刺激| 国产野战对白在线观看| 国产精品 国内视频| 久久久久久久久免费视频了| 欧美日本亚洲视频在线播放| 一进一出抽搐动态| 好男人电影高清在线观看| 波多野结衣一区麻豆| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三区四区五区乱码| 免费看a级黄色片| 久久精品国产99精品国产亚洲性色 | 亚洲精品在线美女| 嫁个100分男人电影在线观看| 亚洲欧美日韩另类电影网站| 大香蕉久久成人网| 午夜精品国产一区二区电影| 精品无人区乱码1区二区| 亚洲少妇的诱惑av| 久久草成人影院| 国产一区在线观看成人免费| 欧美大码av| 一二三四社区在线视频社区8| 国产成人av激情在线播放| 老汉色av国产亚洲站长工具| 丝袜美足系列| 中文字幕另类日韩欧美亚洲嫩草| 国产乱人伦免费视频| 亚洲人成电影观看| 99久久久亚洲精品蜜臀av| 美女 人体艺术 gogo| 亚洲自拍偷在线| 国产亚洲欧美98| 无遮挡黄片免费观看| 亚洲男人的天堂狠狠| 色综合亚洲欧美另类图片| 亚洲情色 制服丝袜| 十八禁网站免费在线| 丁香六月欧美| av在线天堂中文字幕| 女性被躁到高潮视频| 日本 欧美在线| 黑人巨大精品欧美一区二区mp4| 桃红色精品国产亚洲av| 国产高清视频在线播放一区| 十分钟在线观看高清视频www| 97人妻精品一区二区三区麻豆 | 国产高清videossex| 中文字幕高清在线视频| 国产精品二区激情视频| 国产精品野战在线观看| 精品国产乱码久久久久久男人| 岛国在线观看网站| 黄片大片在线免费观看| 麻豆久久精品国产亚洲av| bbb黄色大片| 精品久久久久久,| 国产精品久久电影中文字幕| 日本精品一区二区三区蜜桃| 成人欧美大片| 国产亚洲精品久久久久5区| 国产精品98久久久久久宅男小说| 午夜福利视频1000在线观看 | 日韩有码中文字幕| 波多野结衣av一区二区av| 黄色丝袜av网址大全| 黄色视频不卡| 中文字幕精品免费在线观看视频| 人成视频在线观看免费观看| 悠悠久久av| 女人被狂操c到高潮| 一进一出好大好爽视频| 看黄色毛片网站| 狠狠狠狠99中文字幕| 一边摸一边做爽爽视频免费| bbb黄色大片| 美女大奶头视频| 国产高清videossex| 久久久久久久午夜电影| 成年版毛片免费区| 精品熟女少妇八av免费久了| 欧美在线黄色| 精品国产乱子伦一区二区三区| 亚洲性夜色夜夜综合| 老司机在亚洲福利影院| 又大又爽又粗| 精品不卡国产一区二区三区| 亚洲久久久国产精品| 18禁美女被吸乳视频| 久久久久久久久中文| 亚洲中文av在线| 国产在线精品亚洲第一网站| 国产av一区二区精品久久| 欧美另类亚洲清纯唯美| 他把我摸到了高潮在线观看| 黄色a级毛片大全视频| 叶爱在线成人免费视频播放| 亚洲人成网站在线播放欧美日韩| 欧美色视频一区免费| 久久国产乱子伦精品免费另类| 国产人伦9x9x在线观看| 成在线人永久免费视频| 久久国产精品影院| 日韩一卡2卡3卡4卡2021年| 久久亚洲精品不卡| 成人欧美大片| 久久亚洲真实| 日韩免费av在线播放| 一卡2卡三卡四卡精品乱码亚洲| 亚洲熟妇中文字幕五十中出| 一本综合久久免费| 成人手机av| 日韩成人在线观看一区二区三区| 久久热在线av| 国产伦人伦偷精品视频| 色哟哟哟哟哟哟| 久久草成人影院| 亚洲熟妇中文字幕五十中出| 国产精品乱码一区二三区的特点 | 美女国产高潮福利片在线看| 日韩有码中文字幕| 久久精品aⅴ一区二区三区四区| 中文字幕av电影在线播放| 国产成人精品无人区| 国产成人欧美| 精品欧美一区二区三区在线| 欧美日本中文国产一区发布| 手机成人av网站| bbb黄色大片| 国产av一区二区精品久久| 色播在线永久视频| 在线观看午夜福利视频| 男女午夜视频在线观看| 人人澡人人妻人| 亚洲欧美一区二区三区黑人| 亚洲国产欧美日韩在线播放| 美女扒开内裤让男人捅视频| 视频区欧美日本亚洲| 亚洲欧美精品综合一区二区三区| 18禁黄网站禁片午夜丰满| 婷婷六月久久综合丁香| netflix在线观看网站| 99riav亚洲国产免费| 国产av在哪里看| 中文字幕最新亚洲高清| 中文字幕另类日韩欧美亚洲嫩草| 亚洲 欧美 日韩 在线 免费| 精品欧美一区二区三区在线| 久久久水蜜桃国产精品网| 午夜激情av网站| 三级毛片av免费| 免费av毛片视频| 国产伦一二天堂av在线观看| 国产av又大| avwww免费| 热re99久久国产66热| 久久天躁狠狠躁夜夜2o2o| 最好的美女福利视频网| 亚洲在线自拍视频| 免费人成视频x8x8入口观看| 成人国产一区最新在线观看| 日韩欧美国产在线观看| 精品国产乱码久久久久久男人| 999久久久精品免费观看国产| 婷婷丁香在线五月| 国产亚洲欧美在线一区二区| 在线观看免费视频日本深夜| 在线av久久热| 亚洲精品国产精品久久久不卡| 欧美精品亚洲一区二区| 女人精品久久久久毛片| 亚洲熟妇中文字幕五十中出| 熟妇人妻久久中文字幕3abv| 在线观看日韩欧美| 久久精品国产亚洲av香蕉五月| 亚洲国产日韩欧美精品在线观看 | 亚洲国产中文字幕在线视频| 淫秽高清视频在线观看| 黄网站色视频无遮挡免费观看| 成年女人毛片免费观看观看9| 欧美日本中文国产一区发布| 久热爱精品视频在线9| 日本vs欧美在线观看视频| 日本精品一区二区三区蜜桃| 老司机在亚洲福利影院| 免费在线观看亚洲国产| 天堂动漫精品| 欧美日韩精品网址| 亚洲一区高清亚洲精品| 男女床上黄色一级片免费看| 97人妻天天添夜夜摸| 国产亚洲欧美98| 国产黄a三级三级三级人| 在线观看66精品国产| 亚洲第一欧美日韩一区二区三区| 国产一区二区在线av高清观看| 在线观看免费视频日本深夜| 亚洲第一欧美日韩一区二区三区| av有码第一页| 亚洲成人免费电影在线观看|