• <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窗口
    固體超強酸催化合成丁基糖苷
    亚洲国产精品成人久久小说| 国产精品99久久久久久久久| 大话2 男鬼变身卡| 亚洲av.av天堂| 亚洲精品av麻豆狂野| 97精品久久久久久久久久精品| 欧美成人精品欧美一级黄| 亚洲人成网站在线观看播放| 久久99一区二区三区| 99热网站在线观看| 毛片一级片免费看久久久久| av线在线观看网站| 岛国毛片在线播放| 国产伦精品一区二区三区视频9| 午夜激情久久久久久久| 国产老妇伦熟女老妇高清| 伊人久久国产一区二区| 国产男女超爽视频在线观看| 久久久久久久久大av| 午夜免费鲁丝| 18禁动态无遮挡网站| 人妻一区二区av| 久久免费观看电影| 一级a做视频免费观看| 久久 成人 亚洲| 大片免费播放器 马上看| 国产欧美日韩一区二区三区在线 | 亚洲av男天堂| 国产精品久久久久成人av| 久久精品国产鲁丝片午夜精品| 久久毛片免费看一区二区三区| 大陆偷拍与自拍| 人人妻人人澡人人看| xxxhd国产人妻xxx| 国产精品久久久久久久久免| 女人久久www免费人成看片| 久久久国产成人免费| 91老司机精品| 国产高清视频在线播放一区| 亚洲精品美女久久久久99蜜臀| 麻豆av在线久日| 色尼玛亚洲综合影院| 老司机在亚洲福利影院| 欧美变态另类bdsm刘玥| 中亚洲国语对白在线视频| av网站在线播放免费| 国产老妇伦熟女老妇高清| 国产又色又爽无遮挡免费看| 老鸭窝网址在线观看| 成人三级做爰电影| 在线观看免费视频日本深夜| 99国产精品免费福利视频| 黄片大片在线免费观看| 久久99热这里只频精品6学生| 9热在线视频观看99| 久久人人97超碰香蕉20202| 国产亚洲欧美精品永久| 黑人猛操日本美女一级片| 99国产精品一区二区三区| 97在线人人人人妻| av视频免费观看在线观看| 午夜福利免费观看在线| 免费观看a级毛片全部| 国产精品免费大片| 亚洲精品国产色婷婷电影| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久精品久久久| 精品高清国产在线一区| 久久久久网色| 999精品在线视频| 免费少妇av软件| 99在线人妻在线中文字幕 | 日韩精品免费视频一区二区三区| 国产不卡av网站在线观看| 91成年电影在线观看| 免费女性裸体啪啪无遮挡网站| 母亲3免费完整高清在线观看| 亚洲精品自拍成人| av一本久久久久| 精品少妇一区二区三区视频日本电影| 欧美黄色片欧美黄色片| 免费在线观看黄色视频的| 在线观看免费高清a一片| 男人舔女人的私密视频| 免费女性裸体啪啪无遮挡网站| 国产97色在线日韩免费| 国产精品98久久久久久宅男小说| 国产不卡一卡二| av超薄肉色丝袜交足视频| 欧美中文综合在线视频| 悠悠久久av| av有码第一页| 日韩大片免费观看网站| 久久久久网色| 91精品三级在线观看| 久9热在线精品视频| 久久国产亚洲av麻豆专区| www.自偷自拍.com| 国产福利在线免费观看视频| 美女扒开内裤让男人捅视频| 变态另类成人亚洲欧美熟女 | 免费在线观看日本一区| 亚洲成人免费av在线播放| 国产免费福利视频在线观看| 久久ye,这里只有精品| 一区二区三区激情视频| 桃花免费在线播放| 国产av国产精品国产| 国产视频一区二区在线看| 欧美黑人欧美精品刺激| 色婷婷av一区二区三区视频| 色婷婷久久久亚洲欧美| 国产精品久久久人人做人人爽| 国产成人精品无人区| 免费av中文字幕在线| 日本黄色视频三级网站网址 | 青草久久国产| 美女福利国产在线| 日韩精品免费视频一区二区三区| 丝袜人妻中文字幕| 中文字幕高清在线视频| 成年动漫av网址| 18在线观看网站| 日本av手机在线免费观看| 大型黄色视频在线免费观看| 啪啪无遮挡十八禁网站| 日日爽夜夜爽网站| h视频一区二区三区| 亚洲精品自拍成人| 大片免费播放器 马上看| 在线观看人妻少妇| 最新的欧美精品一区二区| 看免费av毛片| 丁香欧美五月| 欧美亚洲 丝袜 人妻 在线| 怎么达到女性高潮| 亚洲少妇的诱惑av| 纵有疾风起免费观看全集完整版| 精品第一国产精品| 亚洲美女黄片视频| 中国美女看黄片| 黄色怎么调成土黄色| 在线亚洲精品国产二区图片欧美| 精品免费久久久久久久清纯 | 在线观看66精品国产| 中文字幕人妻丝袜一区二区| 亚洲国产av新网站| 在线观看66精品国产| 国产一区二区三区综合在线观看| 亚洲成国产人片在线观看| 757午夜福利合集在线观看| 两个人免费观看高清视频| 男人操女人黄网站| 一区二区三区乱码不卡18| av在线播放免费不卡| 国产黄频视频在线观看| 成人永久免费在线观看视频 | 久久人人爽av亚洲精品天堂| 色在线成人网| 波多野结衣一区麻豆| 亚洲国产欧美网| 久久狼人影院| 如日韩欧美国产精品一区二区三区| 日日爽夜夜爽网站| 另类亚洲欧美激情| 日韩有码中文字幕| 757午夜福利合集在线观看| 黄片小视频在线播放| 一进一出抽搐动态| 一边摸一边做爽爽视频免费| 叶爱在线成人免费视频播放| 久久人妻av系列| 日韩制服丝袜自拍偷拍| 夜夜爽天天搞| 男男h啪啪无遮挡| 国产又色又爽无遮挡免费看| 热99re8久久精品国产| 下体分泌物呈黄色| 亚洲人成77777在线视频| 午夜久久久在线观看| 久久婷婷成人综合色麻豆| 国内毛片毛片毛片毛片毛片| 国产人伦9x9x在线观看| www.自偷自拍.com| 久久性视频一级片| 欧美日韩国产mv在线观看视频| 久久国产精品男人的天堂亚洲| 一级片免费观看大全| 一级毛片女人18水好多| 97在线人人人人妻| 日韩熟女老妇一区二区性免费视频| 天天操日日干夜夜撸| 欧美人与性动交α欧美精品济南到| 热99re8久久精品国产| 欧美大码av| 色老头精品视频在线观看| 日韩成人在线观看一区二区三区| 成年动漫av网址| 亚洲av电影在线进入| 欧美成狂野欧美在线观看| www.熟女人妻精品国产| 欧美日韩精品网址| 男女无遮挡免费网站观看| 亚洲欧美一区二区三区黑人| 国产成人欧美在线观看 | 色综合婷婷激情| 每晚都被弄得嗷嗷叫到高潮| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 窝窝影院91人妻| 悠悠久久av| 国产精品.久久久| 国产成人av教育| 在线十欧美十亚洲十日本专区| 久久久欧美国产精品| 99久久国产精品久久久| 国产精品久久久久久精品古装| 黑丝袜美女国产一区| 久久久久精品国产欧美久久久| 国产高清视频在线播放一区| 母亲3免费完整高清在线观看| 黄色a级毛片大全视频| 欧美精品高潮呻吟av久久| 欧美日韩视频精品一区| 97人妻天天添夜夜摸| 人人妻人人澡人人爽人人夜夜| 亚洲精品国产一区二区精华液| 丰满饥渴人妻一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 黄色a级毛片大全视频| 欧美在线黄色| 国产麻豆69| 亚洲中文字幕日韩| 交换朋友夫妻互换小说| 在线观看免费高清a一片| 日韩欧美一区二区三区在线观看 | 岛国毛片在线播放| 亚洲五月婷婷丁香| 午夜成年电影在线免费观看| 色老头精品视频在线观看| 三上悠亚av全集在线观看| 日韩熟女老妇一区二区性免费视频| 中文欧美无线码| 免费日韩欧美在线观看| 国产男女超爽视频在线观看| 国产精品电影一区二区三区 | 欧美日韩福利视频一区二区| 久久ye,这里只有精品| 久热爱精品视频在线9| 在线永久观看黄色视频| 国产精品 国内视频| 男女高潮啪啪啪动态图| 麻豆乱淫一区二区| 蜜桃在线观看..| 制服诱惑二区| 91精品国产国语对白视频| 涩涩av久久男人的天堂| 免费观看人在逋| 午夜激情久久久久久久| 国产精品99久久99久久久不卡| 久久亚洲精品不卡| 欧美成狂野欧美在线观看| 黄色毛片三级朝国网站| 老汉色∧v一级毛片| 欧美日本中文国产一区发布| 国产高清国产精品国产三级| 一区福利在线观看| 亚洲精品乱久久久久久| 老熟女久久久| 后天国语完整版免费观看| a级毛片黄视频| 免费高清在线观看日韩| 美女福利国产在线| 亚洲精品粉嫩美女一区| 怎么达到女性高潮| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码| 亚洲精品在线美女| 超色免费av| 欧美激情久久久久久爽电影 | 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利影视在线免费观看| 19禁男女啪啪无遮挡网站| 黑人猛操日本美女一级片| 一个人免费看片子| 成人国产av品久久久| 国产精品偷伦视频观看了| 不卡av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 777久久人妻少妇嫩草av网站| av有码第一页| 中文字幕av电影在线播放| kizo精华| av一本久久久久| 国产又爽黄色视频| 在线av久久热| 嫩草影视91久久| 亚洲天堂av无毛| 亚洲欧美精品综合一区二区三区| 肉色欧美久久久久久久蜜桃| 精品少妇内射三级| 美女高潮到喷水免费观看| av天堂久久9| 国产主播在线观看一区二区| 国产精品久久电影中文字幕 | 午夜福利一区二区在线看| 一本色道久久久久久精品综合| 久久久精品区二区三区| 汤姆久久久久久久影院中文字幕| 9热在线视频观看99| 狠狠婷婷综合久久久久久88av| 久久国产精品人妻蜜桃| 欧美午夜高清在线| 在线播放国产精品三级| 国产精品久久久人人做人人爽| 在线av久久热| 午夜福利视频在线观看免费| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 欧美精品一区二区大全| 一边摸一边抽搐一进一出视频| 日韩熟女老妇一区二区性免费视频| 久久午夜综合久久蜜桃| 亚洲视频免费观看视频| 国产aⅴ精品一区二区三区波| 国产xxxxx性猛交| 99久久精品国产亚洲精品| 少妇粗大呻吟视频| 中文字幕人妻熟女乱码| 亚洲专区中文字幕在线| 精品国产一区二区久久| 国产在线一区二区三区精| 精品国产乱码久久久久久小说| 午夜两性在线视频| 青草久久国产| 免费在线观看影片大全网站| 在线观看免费午夜福利视频| 男人舔女人的私密视频| 国产精品美女特级片免费视频播放器 | 男女免费视频国产| 免费日韩欧美在线观看| 纵有疾风起免费观看全集完整版| 精品久久久精品久久久| 亚洲av电影在线进入| 在线观看免费高清a一片| 欧美日本中文国产一区发布| 免费一级毛片在线播放高清视频 | 精品久久久精品久久久| 日韩欧美免费精品| 高清欧美精品videossex| 欧美日韩亚洲综合一区二区三区_| 啦啦啦中文免费视频观看日本| av视频免费观看在线观看| 国产男女内射视频| 成年动漫av网址| 高清av免费在线| 午夜福利欧美成人| 成人18禁在线播放| 日韩欧美三级三区| 国产99久久九九免费精品| 99国产精品免费福利视频| 久久久久久久久免费视频了| 天堂动漫精品| 免费一级毛片在线播放高清视频 | 欧美性长视频在线观看| 69av精品久久久久久 | 高清毛片免费观看视频网站 | 黄色丝袜av网址大全| 80岁老熟妇乱子伦牲交| 最近最新中文字幕大全电影3 | 国产精品.久久久| 啦啦啦视频在线资源免费观看| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| videosex国产| 久久精品亚洲精品国产色婷小说| 欧美激情高清一区二区三区| 久久久久网色| 久久精品成人免费网站| 天堂俺去俺来也www色官网| 日韩大码丰满熟妇| 欧美亚洲日本最大视频资源| 欧美日韩一级在线毛片| 乱人伦中国视频| 精品久久久久久电影网| 免费黄频网站在线观看国产| 亚洲avbb在线观看| 少妇被粗大的猛进出69影院| 人人澡人人妻人| 国产在线视频一区二区| 精品少妇内射三级| 久久人妻熟女aⅴ| 在线 av 中文字幕| 亚洲av国产av综合av卡| 老汉色∧v一级毛片| 午夜精品久久久久久毛片777| 在线观看www视频免费| 久久性视频一级片| 999久久久国产精品视频| 99热国产这里只有精品6| 下体分泌物呈黄色| 最近最新中文字幕大全免费视频| 午夜福利视频精品| 午夜福利乱码中文字幕| 99热网站在线观看| 丝袜在线中文字幕| 最黄视频免费看| 国产深夜福利视频在线观看| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| 亚洲综合色网址| 久久精品国产亚洲av香蕉五月 | av免费在线观看网站| 亚洲精品自拍成人| 飞空精品影院首页| 精品国产一区二区久久| 日本欧美视频一区| 高清在线国产一区| 日本黄色日本黄色录像| 日日爽夜夜爽网站| 中文字幕另类日韩欧美亚洲嫩草| 老司机亚洲免费影院| 精品福利永久在线观看| 精品国产超薄肉色丝袜足j| 波多野结衣一区麻豆| 成人国产一区最新在线观看| 国产单亲对白刺激| 日韩欧美一区二区三区在线观看 | 成人精品一区二区免费| 精品免费久久久久久久清纯 | av欧美777| 在线播放国产精品三级| 亚洲熟妇熟女久久| 国产在视频线精品| a级毛片在线看网站| 久久久精品区二区三区| 国产成人一区二区三区免费视频网站| videos熟女内射| 久久久国产一区二区| 日韩欧美一区二区三区在线观看 | 99九九在线精品视频| 免费观看人在逋| 欧美日韩精品网址| 亚洲熟女精品中文字幕| 后天国语完整版免费观看| 久久久久久免费高清国产稀缺| 男人操女人黄网站| 国产在线一区二区三区精| 中文字幕色久视频| 宅男免费午夜| 免费一级毛片在线播放高清视频 | 色播在线永久视频| 18禁美女被吸乳视频| 亚洲专区国产一区二区| 中文字幕人妻熟女乱码| 亚洲欧美色中文字幕在线| 欧美中文综合在线视频| 国产在线视频一区二区| 日韩免费高清中文字幕av| 国产成人精品无人区| 建设人人有责人人尽责人人享有的| 亚洲熟女精品中文字幕| 精品亚洲成a人片在线观看| 成人三级做爰电影| 深夜精品福利| 一区二区三区乱码不卡18| 人成视频在线观看免费观看| 999久久久精品免费观看国产| 久久香蕉激情| 欧美精品高潮呻吟av久久| 黄色丝袜av网址大全| 国产欧美亚洲国产| 久久精品国产综合久久久| 岛国毛片在线播放| 91大片在线观看| 久久中文看片网| 成人精品一区二区免费| 交换朋友夫妻互换小说| 免费在线观看影片大全网站| 国产区一区二久久| 久久精品国产亚洲av香蕉五月 | 免费在线观看影片大全网站| 国产成人欧美在线观看 | 亚洲国产欧美一区二区综合| 亚洲国产成人一精品久久久| 777米奇影视久久| 1024香蕉在线观看| 99久久99久久久精品蜜桃| 久久久久久久精品吃奶| 亚洲精品国产区一区二| 狂野欧美激情性xxxx| 国产不卡av网站在线观看| 真人做人爱边吃奶动态| 欧美精品一区二区大全| 国产高清视频在线播放一区| 精品午夜福利视频在线观看一区 | av网站在线播放免费| 高清视频免费观看一区二区| av电影中文网址| 又大又爽又粗| 久久精品国产a三级三级三级| 久久精品国产99精品国产亚洲性色 | 满18在线观看网站| 国产精品久久久久久精品古装| 午夜老司机福利片| 久久精品国产亚洲av香蕉五月 | 国产高清视频在线播放一区| 亚洲综合色网址| 一二三四社区在线视频社区8| 亚洲av成人一区二区三| 视频在线观看一区二区三区| 天天躁夜夜躁狠狠躁躁| 大片电影免费在线观看免费| 国产在线视频一区二区| 久久久久网色| av福利片在线| 国产又色又爽无遮挡免费看| 窝窝影院91人妻| 男女下面插进去视频免费观看| 日本黄色视频三级网站网址 | 欧美乱码精品一区二区三区| 高清毛片免费观看视频网站 | 亚洲欧美一区二区三区久久| 这个男人来自地球电影免费观看| 国产欧美日韩一区二区三区在线| 中文字幕色久视频| 亚洲精品美女久久久久99蜜臀| 久久精品国产综合久久久| 亚洲一区二区三区欧美精品| 欧美+亚洲+日韩+国产| 亚洲欧洲日产国产| 色94色欧美一区二区| av免费在线观看网站| 热99re8久久精品国产| 久久精品国产亚洲av香蕉五月 | 嫩草影视91久久| 成人18禁在线播放| 天天添夜夜摸| 午夜两性在线视频| 一本色道久久久久久精品综合| 久久天堂一区二区三区四区| 亚洲av片天天在线观看| 国产成人影院久久av| 精品亚洲成a人片在线观看| 久久国产精品人妻蜜桃| 免费观看av网站的网址| 激情视频va一区二区三区| 亚洲精品乱久久久久久| 精品人妻熟女毛片av久久网站| 搡老熟女国产l中国老女人| 桃红色精品国产亚洲av| 精品亚洲成国产av| 天堂8中文在线网| 免费女性裸体啪啪无遮挡网站| 亚洲三区欧美一区| 日本wwww免费看| 少妇被粗大的猛进出69影院| 咕卡用的链子| 国产黄色免费在线视频| 免费黄频网站在线观看国产| 黄色怎么调成土黄色| 69av精品久久久久久 | 久久久久国产一级毛片高清牌| 亚洲,欧美精品.| 中文字幕精品免费在线观看视频| 免费人妻精品一区二区三区视频| 国产深夜福利视频在线观看| 国产高清国产精品国产三级| 中文字幕人妻丝袜制服| 丝袜在线中文字幕| 国产亚洲午夜精品一区二区久久| 日本黄色日本黄色录像| 精品福利观看| 99久久精品国产亚洲精品| 男女床上黄色一级片免费看| 色老头精品视频在线观看| 三级毛片av免费| 欧美老熟妇乱子伦牲交| 我要看黄色一级片免费的| 男人操女人黄网站| h视频一区二区三区| 亚洲午夜理论影院| 欧美成狂野欧美在线观看| 精品视频人人做人人爽| 午夜日韩欧美国产| 亚洲人成电影观看| 搡老岳熟女国产| 少妇猛男粗大的猛烈进出视频| 亚洲国产欧美一区二区综合| 新久久久久国产一级毛片| 国产成人av激情在线播放| 69av精品久久久久久 | 国产男女内射视频| 一个人免费在线观看的高清视频| 侵犯人妻中文字幕一二三四区| 国产一区二区在线观看av| 日本精品一区二区三区蜜桃| 久久久久网色| 老司机午夜十八禁免费视频| 免费在线观看完整版高清| 手机成人av网站| 亚洲国产av影院在线观看| 中文字幕人妻丝袜一区二区| 又紧又爽又黄一区二区| 久久久久久亚洲精品国产蜜桃av| 久久中文看片网| 色老头精品视频在线观看| 日本vs欧美在线观看视频| 老汉色av国产亚洲站长工具| 久久国产精品影院| 亚洲成国产人片在线观看| 国产精品.久久久| 精品国产超薄肉色丝袜足j|