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

    Box-Behnken設(shè)計(jì)-響應(yīng)面法結(jié)合基準(zhǔn)關(guān)聯(lián)度和AHP-熵權(quán)法優(yōu)化經(jīng)典名方小承氣湯的提取工藝

    2022-10-21 06:54:04趙玥瑛王昌海張澤康郭明雪楊雪穎張馨雨杜守穎
    中草藥 2022年20期
    關(guān)鍵詞:膏率名方蘆薈

    趙玥瑛,王昌海,張澤康,張 晴,郭明雪,孫 萌,席 鋮,楊雪穎,董 爽,張馨雨,杜守穎,陸 洋,白 潔

    Box-Behnken設(shè)計(jì)-響應(yīng)面法結(jié)合基準(zhǔn)關(guān)聯(lián)度和AHP-熵權(quán)法優(yōu)化經(jīng)典名方小承氣湯的提取工藝

    趙玥瑛,王昌海,張澤康,張 晴,郭明雪,孫 萌,席 鋮,楊雪穎,董 爽,張馨雨,杜守穎,陸 洋*,白 潔*

    北京中醫(yī)藥大學(xué)中藥學(xué)院,北京 102488

    應(yīng)用Box-Behnken設(shè)計(jì)-響應(yīng)面法,結(jié)合基準(zhǔn)關(guān)聯(lián)度和層次分析法(analytic hierarchy process,AHP)-熵權(quán)法優(yōu)化經(jīng)典名方小承氣湯的提取工藝。以蘆薈大黃素和大黃酚、厚樸酚、辛弗林含量、出膏率及指紋圖譜相似度為關(guān)鍵質(zhì)量屬性(critical quality attributes,CQAs),以加水倍量、提取時(shí)間、提取次數(shù)為關(guān)鍵工藝參數(shù)(critical process parameters,CPPs),采用單因素實(shí)驗(yàn)確定各因素水平,Box-Behnken設(shè)計(jì)優(yōu)化提取工藝參數(shù)并進(jìn)行驗(yàn)證。最后,計(jì)算不同提取參數(shù)的各評(píng)價(jià)指標(biāo)下各樣品與基準(zhǔn)樣品的基準(zhǔn)關(guān)聯(lián)度,采用AHP-熵權(quán)法確定各種評(píng)價(jià)指標(biāo)的權(quán)重系數(shù),進(jìn)行綜合評(píng)分,預(yù)測(cè)最佳提取工藝。方差分析結(jié)果顯示,2次項(xiàng)回歸模型的方差顯著(<0.01),且失擬值不顯著,表明所建模型具有統(tǒng)計(jì)學(xué)意義,各因素與響應(yīng)值之間的關(guān)系可以用所建模型進(jìn)行函數(shù)化。確定最佳提取工藝為全方飲片加6倍量水,提取30 min,提取1次。平行3次的驗(yàn)證試驗(yàn)綜合評(píng)分均值為95.37,RSD為2.31%,且符合基準(zhǔn)樣品質(zhì)控標(biāo)準(zhǔn)?;谫|(zhì)量源于設(shè)計(jì)(quality by design,QbD)理念、基準(zhǔn)關(guān)聯(lián)度和AHP-熵權(quán)法篩選的經(jīng)典名方小承氣湯提取工藝參數(shù)穩(wěn)健可靠,同時(shí)為其他經(jīng)典名方的制備工藝開發(fā)提供了參考。

    質(zhì)量源于設(shè)計(jì);經(jīng)典名方;小承氣湯;Box-Behnken設(shè)計(jì)-響應(yīng)面法;基準(zhǔn)關(guān)聯(lián)度;層次分析法-熵權(quán)法;蘆薈大黃素;大黃酚;厚樸酚;辛弗林;指紋圖譜;關(guān)鍵質(zhì)量屬性;關(guān)鍵工藝參數(shù)

    經(jīng)典名方小承氣湯(Xiaochengqi Decoction)出自《傷寒論》(東漢·張仲景),被收載于《古代經(jīng)典名方目錄(第一批)》[1]。全方由大黃、厚樸、枳實(shí)3味藥組成,具有瀉熱通便、消脹除滿的功效,主治痞、滿、實(shí)而不燥之陽(yáng)明腑實(shí)證,癥見潮熱、譫語(yǔ)、大便不通,效果顯著[2]?,F(xiàn)代應(yīng)用主要集中在便秘[3-6]、術(shù)后胃腸功能紊亂[7-10]、和術(shù)后早期炎性腸梗阻[11-12]。臨床上熱性感染性疾病中期,如大葉性肺炎、流行性腦膜炎等出現(xiàn)便秘,胸腹脹滿,高熱譫語(yǔ),口渴,舌苔老黃等實(shí)熱癥候時(shí),用小承氣湯酌加清熱解毒之品,收效頗著[13]。

    經(jīng)典名方中藥復(fù)方制劑的開發(fā)需基于基準(zhǔn)樣品的質(zhì)量標(biāo)準(zhǔn),應(yīng)保證關(guān)鍵質(zhì)量屬性與基準(zhǔn)樣品的一致性。但由于工業(yè)化生產(chǎn)時(shí)提取方式、提取容器等的變化以及受成本、效率的影響,提取時(shí)間可能發(fā)生較大改變,導(dǎo)致中試產(chǎn)品和基準(zhǔn)樣品的質(zhì)量不一致,因此,實(shí)現(xiàn)現(xiàn)代提取工藝與傳統(tǒng)制法的質(zhì)量一致性是經(jīng)典名方研究的關(guān)鍵所在,需統(tǒng)籌考慮提取工藝各因素對(duì)關(guān)鍵質(zhì)量屬性的影響。質(zhì)量源于設(shè)計(jì)(quality by design,QbD)是一種系統(tǒng)的研發(fā)方法,需要在可靠的科學(xué)和質(zhì)量風(fēng)險(xiǎn)管理基礎(chǔ)上,預(yù)先定義好目標(biāo),強(qiáng)調(diào)對(duì)產(chǎn)品和工藝的理解以及對(duì)工藝的控制[14-15]。本研究采用基準(zhǔn)關(guān)聯(lián)度(standard relation,SR)來(lái)評(píng)價(jià)現(xiàn)代提取工藝下各樣品與基準(zhǔn)樣品間質(zhì)量的一致性,SR值越接近100%,則該樣品與基準(zhǔn)樣品越接近。層次分析法(analytic hierarchy process,AHP)-熵權(quán)法將主觀賦權(quán)法和客觀賦權(quán)法相結(jié)合,用于確定不同測(cè)定指標(biāo)的權(quán)重系數(shù),所得評(píng)價(jià)結(jié)果更加真實(shí)科學(xué)[16]。本研究基于QbD理念篩選小承氣湯的最優(yōu)提取工藝參數(shù)。計(jì)算不同提取參數(shù)的各種評(píng)價(jià)指標(biāo)下各樣品與基準(zhǔn)樣品的基準(zhǔn)關(guān)聯(lián)度,采用AHP-熵權(quán)法確定各評(píng)價(jià)指標(biāo)的權(quán)重系數(shù),進(jìn)行綜合評(píng)分,結(jié)合Box-Behnken響應(yīng)面法進(jìn)行模型預(yù)測(cè),以保證在所篩選的提取工藝參數(shù)下制備的樣品與基準(zhǔn)樣品保持質(zhì)量一致,確保提取工藝的穩(wěn)定可靠,為小承氣湯的后續(xù)開發(fā)和工業(yè)化生產(chǎn)夯實(shí)基礎(chǔ)。

    1 儀器與試藥

    1.1 儀器

    WT3003型電子天平,杭州萬(wàn)特衡器有限公司;JM-B10002型電子天平,余姚市紀(jì)銘稱重校驗(yàn)設(shè)備有限公司;BSA 224S型電子分析天平、BT 125D型電子分析天平,賽多利斯科學(xué)儀器(北京)有限公司;HH-6型電熱恒溫水浴鍋,北京科偉永興儀器有限公司;DZF-6051型真空干燥器,北京利康達(dá)圣科技有限公司;BY-400C-1型醫(yī)用離心機(jī)、G20型醫(yī)用離心機(jī),北京白洋醫(yī)療器械有限公司;SHB-III型循環(huán)水式多用真空泵,鄭州長(zhǎng)城科工貿(mào)有限公司;DHG-9140 A型電熱恒溫鼓風(fēng)干燥箱,上海一恒科學(xué)儀器有限公司;KQ5200DA型數(shù)控超聲波清洗器,昆山市超聲儀器有限公司;DZTW型電子調(diào)溫電熱套500 mL,北京中儀泓瑞科技發(fā)展有限公司;MH-2000可調(diào)式電熱套2000 mL、MH-3000可調(diào)式電熱套3000 mL、MH-5000可調(diào)式電熱套5000 mL,北京科偉永興儀器有限公司;Thermo UItimate 3000高效液相色譜儀,DAD檢測(cè)器,CM7.2色譜工作站,賽默飛世爾科技(中國(guó))有限公司。

    1.2 試藥

    酒洗大黃飲片(批號(hào)1-21030301,產(chǎn)地四川阿壩縣)、姜厚樸飲片(批號(hào)1-20042001,產(chǎn)地陜西城固縣)、清炒枳實(shí)飲片(批號(hào)1-21030601,產(chǎn)地湖南漢壽縣),均購(gòu)自河北萬(wàn)歲藥業(yè)有限公司,分別經(jīng)北京中醫(yī)藥大學(xué)中藥鑒定系劉春生教授鑒定為蓼科大黃屬植物藥用大黃Baill.的干燥根和根莖,木蘭科木蘭屬植物厚樸Rehd. et Wils.的干燥干皮、根皮及枝皮和蕓香科柑橘屬植物酸橙L.的干燥幼果的炮制加工品。

    對(duì)照品蘆薈大黃素(批號(hào)110795-202011,質(zhì)量分?jǐn)?shù)97.5%)、大黃酚(批號(hào)110796-201922,質(zhì)量分?jǐn)?shù)99.4%)、厚樸酚(批號(hào)110729-202015,質(zhì)量分?jǐn)?shù)99.0%)、辛弗林(批號(hào)110727-202110,質(zhì)量分?jǐn)?shù)99.8%),中國(guó)食品藥品檢定研究院。

    乙腈、甲醇,色譜級(jí),美國(guó)Fisher公司;純凈水,杭州娃哈哈集團(tuán)有限公司;甲酸,色譜級(jí),天津市大茂化學(xué)試劑廠;去離子水。

    2 方法與結(jié)果

    2.1 小承氣湯提取液的制備[17]

    稱取酒洗大黃55.2 g、姜厚樸27.6 g、清炒枳實(shí)39.0 g,加入一定倍量的去離子水,加熱回流提取,用1層300目尼龍篩網(wǎng)濾過(guò),調(diào)整濾液體積至適宜體積即得小承氣湯。

    2.2 小承氣湯指標(biāo)成分的含量測(cè)定

    2.2.1 對(duì)照品溶液的制備 精密稱取蘆薈大黃素對(duì)照品3.63 mg、大黃酚對(duì)照品5.59 mg、厚樸酚對(duì)照品6.62 mg,用甲醇充分溶解后轉(zhuǎn)移至25 mL量瓶中;精密稱取辛弗林對(duì)照品9.47 mg,用甲醇充分溶解后轉(zhuǎn)移至25 mL量瓶中;分別加甲醇溶解定容至刻度,配成質(zhì)量濃度分別為141.57、222.26、262.15、378.04 μg/mL的對(duì)照品母液。再精密移取1 mL母液至25 mL量瓶中,加甲醇溶解定容至刻度,配成質(zhì)量濃度分別為5.66、8.89、10.49、15.12 μg/mL的對(duì)照品溶液。

    2.2.2 供試品溶液的制備[17]移取小承氣湯提取液3 mL于15 mL離心管內(nèi),加入7 mL甲醇,搖勻后,5000 r/min離心10 min,取上清液過(guò)0.45 μm微孔濾膜,即得供試品溶液1。精密移取上述離心后的上清液5 mL,通過(guò)聚酰胺柱(60~90目,2.5 g,內(nèi)徑1.5 cm,干法裝柱),用水20 mL洗脫,收集洗脫液于25 mL量瓶?jī)?nèi),加水定容搖勻,過(guò)0.45 μm微孔濾膜,得供試品溶液2。

    2.2.3 色譜條件[17]

    (1)供試品溶液1:色譜柱為Thermo Hypersil Gold C18(250 mm×4.6 mm,5 μm)柱;流動(dòng)相為乙腈-0.1%甲酸水溶液,洗脫梯度:0~10 min,15%乙腈;10~30 min,15%~20%乙腈;30~95 min,20%~30%乙腈;95~135 min,30%~60%乙腈;135~140 min,60%~95%乙腈;體積流量1 mL/min;檢測(cè)波長(zhǎng)為260 nm;柱溫30 ℃;進(jìn)樣體積20 μL。

    (2)供試品溶液2:色譜柱為Thermo Umisil C18(250 mm×4.6 mm,5 μm)柱;流動(dòng)相為甲醇-磷酸二氫鉀溶液(取磷酸二氫鉀0.6 g,十二烷基磺酸鈉1.0 g,冰醋酸1 mL,加水溶解并稀釋至1000 mL),等度洗脫:0~20 min,51%甲醇;體積流量為1 mL/min;檢測(cè)波長(zhǎng)為275 nm;柱溫30 ℃;進(jìn)樣體積30 μL。

    2.3 出膏率的測(cè)定

    將小承氣湯提取液濃縮到相對(duì)密度1.20~1.25,再轉(zhuǎn)移至已恒定質(zhì)量的250 mL蒸發(fā)皿中,水浴蒸干,60 ℃真空干燥72 h,即得小承氣湯干膏粉,計(jì)算出膏率(出膏率=/,為小承氣湯干膏粉的質(zhì)量,為全方飲片質(zhì)量)。

    2.4 指紋圖譜的測(cè)定[17]

    2.4.1 供試品溶液的制備 同“2.2.2”項(xiàng)下供試品溶液1的制備方法。

    2.4.2 色譜條件 同“2.2.3”項(xiàng)下供試品溶液1的色譜條件。

    2.5 提取工藝研究

    2.5.1 單因素考察

    (1)加水量:稱取全方飲片置于圓底燒瓶中,分別加水6、8、10、12、14、20倍,加熱回流1 h,依法測(cè)定蘆薈大黃素和大黃酚、厚樸酚、辛弗林含量及出膏率,平行操作3份,根據(jù)實(shí)驗(yàn)結(jié)果平均值確定加水倍量的低水平和高水平。由表1可知,加6~14倍量水時(shí),各指標(biāo)整體呈上升趨勢(shì),超過(guò)14倍量水后趨于平穩(wěn),故選加水量為6倍和14倍。

    (2)提取時(shí)間:稱取全方飲片置于圓底燒瓶中,加水8倍,分別回流提取30、60、90、120、150 min,依法測(cè)定蘆薈大黃素和大黃酚、厚樸酚、辛弗林含量及出膏率,平行操作3份,根據(jù)實(shí)驗(yàn)結(jié)果平均值確定提取時(shí)間的低水平和高水平。由表2可知,提取時(shí)間為30~120 min時(shí),除90 min的蘆薈大黃素和大黃酚含量低于60 min外,其他指標(biāo)均呈上升趨勢(shì),超過(guò)120 min后,蘆薈大黃素和大黃酚、厚樸酚的含量下降,辛弗林含量增加,出膏率無(wú)明顯變化,故選提取時(shí)間為30 min和120 min。

    表1 加水量對(duì)指標(biāo)成分含量和出膏率的影響

    Table 1 Effects of water addition on index component content and paste-forming rate

    加水量/倍質(zhì)量分?jǐn)?shù)/(mg?g?1)出膏率/% 蘆薈大黃素和大黃酚厚樸酚辛弗林 60.110.090.3414.55 80.160.120.3616.18 100.200.150.4417.66 120.200.190.4818.41 140.250.200.4818.37 200.250.250.5118.64

    表2 提取時(shí)間對(duì)指標(biāo)成分含量和出膏率的影響

    Table 2 Effects of extraction time on index component content and paste-forming rate

    提取時(shí)間/min質(zhì)量分?jǐn)?shù)/(mg?g?1)出膏率/% 蘆薈大黃素和大黃酚厚樸酚辛弗林 300.160.090.2814.82 600.160.120.3616.18 900.150.130.4417.19 1200.190.160.4920.17 1500.180.140.5920.18

    2.5.2 基準(zhǔn)樣品的制備 小承氣湯出自東漢張仲景所著的《傷寒論》,原文記載制法為“上三味,以水四升,煮取一升二合,去滓,分溫二服”[2]。綜合古籍和實(shí)際生產(chǎn)要求,確定全方的煎煮工藝為按處方量稱取3味飲片于陶瓷鍋中,加800 mL去離子水,武火(1600 W)煮沸轉(zhuǎn)文火,文火(600 W)煎煮約47 min至240 mL,用1層300目尼龍篩網(wǎng)過(guò)濾,調(diào)整濾液體積至240 mL,得水煎液。

    分別取5 mL水煎液于25 mL西林瓶中,放入?20 ℃冰箱中預(yù)冷凍12 h后,置于?80 ℃預(yù)冷2 h的冷凍干燥機(jī),凍干溫度為?80 ℃,真空度為(5±1)Pa,干燥72 h,即得小承氣湯基準(zhǔn)樣品[17]。

    2.5.3 Box-Behnken試驗(yàn)

    (1)關(guān)鍵質(zhì)量屬性(critical quality attributes,CQAs)和關(guān)鍵工藝參數(shù)(critical process parameters,CPPs)的確定:小承氣湯主治陽(yáng)明邪熱與腸中糟粕互結(jié)但燥結(jié)不甚之證,方中與此主治病癥相關(guān)的藥效成分為蘆薈大黃素、大黃酚、厚樸酚和辛弗林[18],同時(shí)這4個(gè)成分分別歸屬不同藥材,可以全面控制小承氣湯基準(zhǔn)樣品的質(zhì)量,故以4個(gè)指標(biāo)性成分含量(1為蘆薈大黃素和大黃酚含量,2為厚樸酚含量,3為辛弗林含量)、出膏率(4)和指紋圖譜相似度(5)作為小承氣湯制備工藝的CQAs[19-22]。基于經(jīng)驗(yàn)和已發(fā)表文獻(xiàn),選擇提取時(shí)間(1)、加水倍量(2)、提取次數(shù)(3)作為小承氣湯提取工藝的CPPs[23-25]。

    (2)Box-Behnken試驗(yàn)設(shè)計(jì):稱取與基準(zhǔn)樣品相同處方量的同一批飲片,通過(guò)單因素實(shí)驗(yàn)結(jié)果確定提取時(shí)間(1)和加水倍量(2)的范圍,依據(jù)經(jīng)驗(yàn)將提取次數(shù)(3)設(shè)計(jì)為1、2、3次,因素水平見表3。采用Design Expert 10.0軟件,進(jìn)行3因素3水平17個(gè)試驗(yàn)點(diǎn)的設(shè)計(jì),其中12個(gè)析因點(diǎn),5個(gè)中心點(diǎn)。

    (3)試驗(yàn)結(jié)果與響應(yīng)面分析:Box-Behnken試驗(yàn)設(shè)計(jì)及小承氣湯樣品檢測(cè)結(jié)果見表3。方差分析結(jié)果見表4。2次多項(xiàng)式模型優(yōu)于交互2項(xiàng)式和1項(xiàng)式模型,且總模型顯著(<0.001),失擬值不顯著,表明所建模型有統(tǒng)計(jì)學(xué)意義,各因素與響應(yīng)值之間的關(guān)系可以用此模型函數(shù)化。多元線性回歸方程為=?106.58-47.641-66.132-120.693-16.0012-22.9113-25.7323+24.9312+5.1422+5.1132,模型adj2為0.980 8,表明方程的擬合度較好,pred2為0.893 5,表明模型的預(yù)測(cè)性能良好,因此,可用此模型來(lái)預(yù)測(cè)各指標(biāo)的實(shí)際情況。

    2.6 基準(zhǔn)關(guān)聯(lián)度的原理方法

    2.6.1 確定評(píng)價(jià)對(duì)象 基準(zhǔn)關(guān)聯(lián)度是評(píng)價(jià)某樣品與基準(zhǔn)樣品相似度的關(guān)鍵參數(shù),本方法的基本路線如圖1所示。應(yīng)首先確定樣品評(píng)價(jià)的指標(biāo)數(shù)目和樣品數(shù)目。本實(shí)驗(yàn)將Box-Behnken設(shè)計(jì)的17個(gè)樣品和基準(zhǔn)樣品進(jìn)行對(duì)比,以蘆薈大黃素和大黃酚、厚樸酚、辛弗林的含量和出膏率以及指紋圖譜相似度5個(gè)指標(biāo)為評(píng)價(jià)指標(biāo),以X、S為評(píng)價(jià)對(duì)象,X表示Box-Behnken設(shè)計(jì)的第(=1,2,…,)個(gè)樣本的第(=1,2,…,)個(gè)指標(biāo)下的測(cè)量值,S表示基準(zhǔn)樣品的第(=1,2,…,)個(gè)指標(biāo)下的測(cè)量值。

    2.6.2 計(jì)算相對(duì)偏差(relative deviation,RD)值 RD表示第(=1,2,…,)個(gè)樣本的第(=1,2,…,)個(gè)指標(biāo)下的相對(duì)偏差,能夠更科學(xué)全面地分析不同樣品與基準(zhǔn)樣品之間的相似度。RD值越小,代表X相對(duì)于S的偏差越小,即該樣品與基準(zhǔn)樣品相似度越高。利用公式(1)計(jì)算X相對(duì)于S的RD值,計(jì)算結(jié)果見表5。結(jié)果表明,當(dāng)以蘆薈大黃素和大黃酚、厚樸酚、辛弗林含量以及指紋圖譜相似度為評(píng)價(jià)指標(biāo)時(shí),與基準(zhǔn)樣品質(zhì)量最為接近的試驗(yàn)號(hào)是5,以出膏率為評(píng)價(jià)指標(biāo)時(shí),與基準(zhǔn)樣品質(zhì)量接近的試驗(yàn)號(hào)為14。

    RD=|XS|/S(1)

    2.6.3 計(jì)算基準(zhǔn)關(guān)聯(lián)度(standard relation,SR) 以基準(zhǔn)樣品質(zhì)量為標(biāo)準(zhǔn)值(),以Box-Behnken試驗(yàn)設(shè)計(jì)下的各樣品質(zhì)量為測(cè)得值(),以不同的指標(biāo)為評(píng)價(jià)指標(biāo),按公式(2)計(jì)算所得結(jié)果,稱為基準(zhǔn)關(guān)聯(lián)度,結(jié)果見表6。SR表示第(=1,2,…,)個(gè)樣本的第(=1,2,…,)個(gè)指標(biāo)下的基準(zhǔn)關(guān)聯(lián)度。SR越接近100%,在該指標(biāo)下的工藝參數(shù)制得的樣品與基準(zhǔn)樣品相似度越高。

    表3 Box-Behnken實(shí)驗(yàn)設(shè)計(jì)與響應(yīng)值

    Table 3 Box-Behnken experimental design and response value

    試驗(yàn)編號(hào)X1/minX2/倍X3/次質(zhì)量分?jǐn)?shù)/(mg?g?1)出膏率/%指紋圖譜相似度綜合評(píng)分 Y1Y2Y3 175 (0)14 (+1)1 (?1)0.220.240.7518.650.976?17.61 230 (?1)10 (0)3 (+1)0.360.340.6436.510.969?123.42 375102 (0)0.290.330.8729.260.976?92.36 4301420.300.280.6222.250.969?66.21 5756 (?1)10.130.080.3414.240.98680.78 6120 (+1)1420.370.491.0535.150.964?201.46 775630.300.390.8231.920.961?123.58 81201030.440.550.9939.180.946?256.57 9751020.340.320.8228.750.969?106.30 10751430.440.710.8636.120.979?324.90 11751020.330.330.7828.790.967?106.20 12120620.260.280.6829.270.969?54.79 13751020.350.320.8228.220.970?109.97 14301010.170.130.5114.190.96557.67 15751020.290.390.8626.730.964?118.05 161201010.210.180.5320.880.96816.18 1730620.200.190.5418.560.98216.45 基準(zhǔn)樣品 0.110.100.3411.131.000

    表4 Box-Behnken實(shí)驗(yàn)設(shè)計(jì)方差分析

    Table 4 Variance analysis of Box-Behnken

    誤差來(lái)源模型F值 二次多項(xiàng)式模型交互二項(xiàng)式模型一項(xiàng)式模型 總模型91.95***65.02***71.58*** X184.22***40.38***22.98*** X2162.27***77.80***44.28*** X3540.45***259.10***147.49*** X1X24.752.28 X1X39.744.67 X2X312.285.89 X1212.14 X220.52 X320.51 失擬值4.498.0112.77 Radj20.980 80.960 00.929 7

    ***P<0.001

    圖1 方法流程圖

    SR=1-RD=1-|XS|/S(2)

    2.7 AHP-熵權(quán)法組合權(quán)重的計(jì)算

    2.7.1 AHP主觀賦權(quán)

    (1)構(gòu)建判斷矩陣:根據(jù)君臣佐使原則[26],小承氣湯方中大黃為君藥,味苦寒,主瀉熱祛實(shí),瀉火解毒,以蕩滌胃腸實(shí)熱;厚樸味苦辛溫,溫中行氣,除脹滿,為臣藥;枳實(shí)味苦微寒,破氣消積導(dǎo)滯,以消痞實(shí),為佐使,3藥配合,共奏瀉熱通便,寬中行氣之效[13]。本研究參考《中藥新藥質(zhì)量標(biāo)準(zhǔn)研究技術(shù)指導(dǎo)原則(試行)》要求[27],根據(jù)小承氣湯的處方組成,首選君藥中的有效成分為含量測(cè)定指標(biāo),兼顧與制備工藝、穩(wěn)定性的相關(guān)性,確定以大黃中的蘆薈大黃素和大黃酚(1)、厚樸的厚樸酚(2)、枳實(shí)的辛弗林(3)作為質(zhì)量控制指標(biāo),同時(shí)輔以出膏率(4)和指紋圖譜相似度(5)。指紋圖譜可反映小分子成分輪廓和定性特質(zhì),但與含量測(cè)定的指標(biāo)性成分相比,其所包含的質(zhì)量信息仍較少[28],故將其重要程度視為與出膏率同等重要。故按蘆薈大黃素和大黃酚>厚樸酚>辛弗林>出膏率≈指紋圖譜相似度的順序,再根據(jù)AHP理論判斷矩陣1~9標(biāo)度法,按同一層次內(nèi)個(gè)指標(biāo)的相對(duì)重要程度進(jìn)行打分(表7)。蘆薈大黃素和大黃酚為君藥的主要藥效成分,相對(duì)于厚樸中的厚樸酚較明顯重要,故標(biāo)為4,相對(duì)枳實(shí)中的辛弗林明顯重要,標(biāo)為5,相對(duì)于出膏率和指紋圖譜相似度明顯更為重要,標(biāo)為6,以此類推,根據(jù)公式(3)構(gòu)建判斷矩陣,數(shù)據(jù)見表8,用a表示因素相對(duì)于因素的比較結(jié)果。

    表5 Box-Behnken設(shè)計(jì)的RDi,j

    Table 5 RDij of Box-Behnken test

    試驗(yàn)號(hào)RDij試驗(yàn)號(hào)RDij試驗(yàn)號(hào)RDij Y1Y2Y3Y4Y5Y1Y2Y3Y4Y5Y1Y2Y3Y4Y5 11.001.401.210.680.0271.732.901.411.870.04132.182.201.411.540.03 22.272.400.882.280.0383.004.501.912.520.05140.550.300.500.270.04 31.642.301.561.630.0292.092.201.411.580.03151.642.901.531.400.04 41.731.800.821.000.03103.006.101.532.250.02160.910.800.560.880.03 50.180.200.000.280.01112.002.301.291.590.03170.820.900.590.670.02 62.363.902.092.160.04121.361.801.001.630.03

    表6 Box-Behnken設(shè)計(jì)的SRij

    Table6 SRij of Box-Behnken test

    試驗(yàn)號(hào)SRij試驗(yàn)號(hào)SRij Y1Y2Y3Y4Y5Y1Y2Y3Y4Y5 10.00?40.00?20.5932.4397.6010?200.00?510.00?52.94?124.5397.90 2?127.27?140.0011.76?128.0396.9011?100.00?130.00?29.41?58.6796.70 3?63.64?130.00?55.88?62.8997.6012?36.36?80.000.00?62.9896.90 4?72.73?80.0017.650.0996.9013?118.18?120.00?41.18?53.5597.00 581.8280.00100.0072.0698.601445.4570.0050.0072.5196.50 6?136.36?290.00?108.82?115.8196.4015?63.64?190.00?52.94?40.1696.40 7?72.73?190.00?41.18?86.7996.10169.0920.0044.1212.4096.80 8?200.00?350.00?91.18?152.0294.601718.1810.0041.1833.2498.20 9?109.09?120.00?41.18?58.3196.90

    (3)

    表7 AHP構(gòu)建判斷矩陣的標(biāo)準(zhǔn)[16]

    Table7 Criterion of constructing judgment matrix of AHP[16]

    標(biāo)度含義 1表示2個(gè)因素相比,具有相同重要性 3表示2個(gè)因素相比,前者比后者稍重要 5表示2個(gè)因素相比,前者比后者明顯重要 7表示2個(gè)因素相比,前者比后者強(qiáng)烈重要 9表示2個(gè)因素相比,前者比后者極端重要 2、4、6、8表示上述相鄰判斷的中間值 倒數(shù)若因素e與因素f的重要性之比為aef,那 么因素f與因素e重要性之比為afe=1/aef

    (2)計(jì)算權(quán)重:首先對(duì)判斷矩陣進(jìn)行幾何平均(方根法),按公式(4)計(jì)算得到初始權(quán)重系數(shù)W′,a表示第(=1,2,…,)個(gè)指標(biāo)因素對(duì)第1,2,…,個(gè)指標(biāo)因素的比較結(jié)果,然后按照公式(5)計(jì)算歸一化權(quán)重系數(shù)(主觀權(quán)重系數(shù),WS),得到5個(gè)指標(biāo)的WS,結(jié)果見表9。

    表8 指標(biāo)成對(duì)比較的判斷優(yōu)先矩陣

    Table 8 Decision matrix of paired comparison on indexes

    評(píng)價(jià)指標(biāo)Y1Y2Y3Y4Y5 Y114566 Y21/41344 Y31/51/3122 Y41/61/41/211 Y51/61/41/211

    W′=(a1a2…a)1/n(4)

    表9 組合權(quán)重系數(shù)

    Table 9 Combination weight coefficient

    評(píng)價(jià)指標(biāo)WjSWjOWj Y10.530.200.43 Y20.230.470.45 Y30.110.140.06 Y40.070.180.05 Y50.070.000.00

    WS=W′/W′ (5)

    (3)一致性檢驗(yàn):對(duì)矩陣進(jìn)行一致性檢驗(yàn)。首先根據(jù)公式(6)計(jì)算最大特征根max,再根據(jù)公式(7)計(jì)算一致性指標(biāo)(consistency index,CI),得到CI為0.04,一致性比例(consistency ratio,CR)(CR=CI/RI,RI為自由度指標(biāo))為0.03,均小于0.1,表明該矩陣具有一致性[16]。

    max=(aW′/W′)/(6)

    CI=(max-)/(-1) (7)

    2.7.2 信息熵客觀賦權(quán)

    (1)建立Box-Behnken試驗(yàn)原始數(shù)據(jù)矩陣:現(xiàn)有個(gè)樣本,個(gè)評(píng)價(jià)指標(biāo),原始數(shù)據(jù)矩陣=(X)×n,以蘆薈大黃素和大黃酚、厚樸酚、辛弗林含量,出膏率及指紋圖譜相似度5個(gè)指標(biāo)為評(píng)價(jià)指標(biāo),建立原始數(shù)據(jù)矩陣。

    (2)將原始數(shù)據(jù)轉(zhuǎn)化為概率矩陣:根據(jù)公式(8)將原始數(shù)據(jù)矩陣轉(zhuǎn)化為概率矩陣,P表示在第個(gè)指標(biāo)下第個(gè)樣本的概率。

    PX/X(8)

    (3)計(jì)算各指標(biāo)的信息熵和權(quán)重系數(shù):信息熵(H)越小,代表指標(biāo)下的數(shù)據(jù)離散程度較高,則其所提供的信息量就越大[16]。根據(jù)公式(9)計(jì)算出各指標(biāo)的H,依次為0.984 4、0.963 0、0.988 8、0.985 7、1.000 0,按照公式(10)計(jì)算各指標(biāo)客觀權(quán)重系數(shù)(WO)。

    H=?PlnP,=1/ln(9)

    WO=(1-H)/(1-H) (10)

    2.7.3 組合權(quán)重的確定 采用AHP得到WS,熵權(quán)法得到WO,根據(jù)公式(11)計(jì)算各指標(biāo)組合權(quán)重系數(shù)(W),結(jié)果見表9。

    WWSWO/WSWO(11)

    2.8 綜合評(píng)分

    通過(guò)AHP-熵權(quán)法得到的綜合權(quán)重系數(shù)可知,蘆薈大黃素和大黃酚指標(biāo)占比43%,厚樸酚占比45%,辛弗林占比6%,出膏率占比5%,因此根據(jù)公式(12)計(jì)算綜合評(píng)分,結(jié)果見表3。綜合評(píng)分為100時(shí),所對(duì)應(yīng)工藝參數(shù)下制得的樣品與基準(zhǔn)樣品的質(zhì)量最為一致。Box-Behnken模型預(yù)測(cè)的最佳工藝為加水6倍量,提取時(shí)間30 min,提取1次。

    綜合評(píng)分=SR,Y1×0.43+SR,Y2×0.45+SR,Y3×0.06+SR,Y4×0.05 (12)

    2.9 驗(yàn)證試驗(yàn)

    對(duì)優(yōu)化后的提取工藝參數(shù)進(jìn)行驗(yàn)證,重復(fù)試驗(yàn)3次,測(cè)定蘆薈大黃素和大黃酚、厚樸酚、辛弗林含量,出膏率及指紋圖譜相似度,計(jì)算綜合評(píng)分,結(jié)果見表10。驗(yàn)證試驗(yàn)的綜合評(píng)分均值為95.37(理論預(yù)測(cè)值為98.41),在課題組前期研究建立的小承氣湯基準(zhǔn)樣品質(zhì)控標(biāo)準(zhǔn)范圍內(nèi)(蘆薈大黃素和大黃酚總質(zhì)量分?jǐn)?shù)0.11~0.20 mg/g,厚樸酚質(zhì)量分?jǐn)?shù)0.08~0.15 mg/g,辛弗林質(zhì)量分?jǐn)?shù)0.23~0.40 mg/g;出膏率9.25%~11.25%;指紋圖譜與對(duì)照?qǐng)D譜相似度>0.9)[17],且綜合評(píng)分的RSD值為2.31%,表明本方法建立的數(shù)學(xué)模型具有良好的預(yù)測(cè)性,所確定的提取工藝條件重復(fù)性高。

    表10 驗(yàn)證試驗(yàn)結(jié)果

    Table 10 Validation of experimental results

    試驗(yàn)編號(hào)質(zhì)量分?jǐn)?shù)/(mg?g?1)Y4/%Y5綜合評(píng)分RSD/% Y1Y2Y3 10.120.100.3413.300.98594.122.31 20.110.100.3413.550.98997.91 30.120.100.3413.400.98694.07 基準(zhǔn)樣品0.110.100.3411.131.000

    3 討論

    在確定了處方、制法、用法等關(guān)鍵信息后,嚴(yán)格按照傳統(tǒng)工藝制得的基準(zhǔn)樣品質(zhì)量是衡量經(jīng)典名方小承氣湯中試產(chǎn)品及商業(yè)化規(guī)模產(chǎn)品質(zhì)量的標(biāo)準(zhǔn)。古籍記載小承氣湯以水為提取溶劑,最終根據(jù)藥液量來(lái)確定煎煮時(shí)間,與回流提取的實(shí)際工業(yè)生產(chǎn)情況不符。因此本研究對(duì)加水量、提取時(shí)間、提取次數(shù)進(jìn)行深入考察,以基準(zhǔn)樣品的質(zhì)控標(biāo)準(zhǔn)為指標(biāo)來(lái)對(duì)提取工藝進(jìn)行整體評(píng)價(jià),確定與傳統(tǒng)工藝基本一致的小承氣湯現(xiàn)代提取工藝參數(shù)。

    本研究基于QbD理念,以現(xiàn)行標(biāo)準(zhǔn)規(guī)范為參照,銜接古籍記載和現(xiàn)行規(guī)范,建立了以蘆薈大黃素和大黃酚、厚樸酚、辛弗林、出膏率和指紋圖譜相似度為主的質(zhì)量控制體系,確定以加水倍量、提取時(shí)間、提取次數(shù)作為關(guān)鍵工藝參數(shù),采用單因素實(shí)驗(yàn)確定各因素水平,進(jìn)一步應(yīng)用Box-Behnken設(shè)計(jì)優(yōu)化經(jīng)典名方小承氣湯提取工藝。AHP是人為確定各個(gè)評(píng)價(jià)指標(biāo)的權(quán)重系數(shù),受主觀影響較大;熵權(quán)法根據(jù)數(shù)據(jù)的變化規(guī)律判斷評(píng)價(jià)指標(biāo)的重要程度,但不能體現(xiàn)方劑各藥味君臣佐使的配伍關(guān)系。本研究結(jié)合基準(zhǔn)關(guān)聯(lián)度和AHP-熵權(quán)法對(duì)質(zhì)量評(píng)價(jià)體系進(jìn)行賦權(quán),明確5個(gè)指標(biāo)在質(zhì)量控制體系中的權(quán)重系數(shù),即可通過(guò)綜合評(píng)分及模型預(yù)測(cè)確定與基準(zhǔn)樣品質(zhì)量最為接近的提取工藝參數(shù),所得數(shù)據(jù)更加科學(xué)全面。

    基于QbD理念結(jié)合基準(zhǔn)關(guān)聯(lián)度和AHP-熵權(quán)法優(yōu)化得到的小承氣湯提取工藝參數(shù)穩(wěn)健可靠,符合經(jīng)典名方“傳承精華、古為今用、古今銜接”的基本研究原則,實(shí)現(xiàn)了其由古方向現(xiàn)代工業(yè)化生產(chǎn)轉(zhuǎn)化的關(guān)鍵一步,為其它經(jīng)典名方的現(xiàn)代提取工藝研究提供了理論參考,同時(shí)也為目前工業(yè)生產(chǎn)中如何保證中間體或制劑與基準(zhǔn)樣品間的一致性提供了評(píng)價(jià)方法。

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

    [1] 國(guó)家中醫(yī)藥管理局. 國(guó)家中醫(yī)藥管理局關(guān)于發(fā)布《古代經(jīng)典名方目錄(第一批)》的通知[EB/OL]. [2018-04-13]. http://kjs.satcm.gov.cn/zhengcewenjian/2018-04-16/7107. html.

    [2] 姜建國(guó). 傷寒論 [M]. 北京: 中國(guó)中醫(yī)藥出版社, 2004: 131, 189.

    [3] 胡時(shí)友. 小承氣湯配合中藥穴位貼敷治療中風(fēng)后便秘的臨床效果 [J]. 中外醫(yī)學(xué)研究, 2018, 16(2): 139-141.

    [4] 李梅, 歐陽(yáng)群, 吳小宇. 小承氣湯和厚樸三物湯藥理作用與臨床應(yīng)用探析 [J]. 亞太傳統(tǒng)醫(yī)藥, 2015, 11(12): 111-112.

    [5] 沈亞紅, 童樹洪. 小承氣湯與厚樸三物湯的藥理作用分析和臨床應(yīng)用比較 [J]. 中國(guó)藥業(yè), 2013, 22(6): 122-123.

    [6] 張紅, 陳孟溪. 小承氣湯治療癌痛鎮(zhèn)痛藥胃腸副作用40例 [J]. 新中醫(yī), 1996, 28(2): 52.

    [7] 張良清, 高海鴻, 馮時(shí)儐, 等. 胃腸手術(shù)前后應(yīng)用小承氣湯對(duì)康復(fù)的影響 [J]. 廣西醫(yī)學(xué), 2013, 35(2): 184-186.

    [8] 蔣旭華. 小承氣湯口服加穴位按摩促進(jìn)腹部術(shù)后腸蠕動(dòng)恢復(fù)的臨床護(hù)理 [J]. 中國(guó)醫(yī)學(xué)工程, 2011, 19(5): 87-88.

    [9] 姜玉芳. 小承氣湯預(yù)防術(shù)后病人腸脹氣的臨床觀察 [J]. 護(hù)理研究, 2003, 17(15): 908-909.

    [10] 龍期伯, 徐晉. 小承氣湯保留灌腸防治胃切除術(shù)后殘留排空延遲癥176例 [J]. 南京中醫(yī)藥大學(xué)學(xué)報(bào), 1997, 13(3): 15-16.

    [11] 王素利, 王存志. 中西醫(yī)結(jié)合治療術(shù)后早期炎性腸梗阻療效觀察 [J]. 中國(guó)中醫(yī)急癥, 2014, 23(4): 731-732.

    [12] 連愛(ài)霞. 小承氣湯保留灌腸輔助治療術(shù)后早期炎性腸梗阻的療效觀察 [J]. 中國(guó)中醫(yī)急癥, 2013, 22(12): 2128-2129.

    [13] 陳潮祖, 李大琦. 中醫(yī)方劑與治法 [M]. 成都: 四川科學(xué)技術(shù)出版社, 1984: 72.

    [14] 陳澤麒, 趙芳, 袁瑋, 等. 基于質(zhì)量源于設(shè)計(jì)(QbD)理念研究參麥注射液醇提水沉工藝 [J]. 中草藥, 2020, 51(17): 4433-4441.

    [15] 戴勝云, 徐冰, 史新元, 等. 中藥質(zhì)量源于設(shè)計(jì)方法和應(yīng)用: 產(chǎn)品和工藝設(shè)計(jì) [J]. 世界中醫(yī)藥, 2018, 13(3): 550-555.

    [16] 代珊, 李帥, 張愛(ài)軍, 等. 基于基準(zhǔn)關(guān)聯(lián)度和AHP-熵權(quán)法綜合評(píng)價(jià)經(jīng)典名方小續(xù)命湯古今提取工藝 [J]. 中草藥, 2022, 53(3): 726-734.

    [17] 齊琪, 趙玥瑛, 張晴, 等. 經(jīng)典名方小承氣湯的物質(zhì)基準(zhǔn)量值傳遞研究 [J]. 中草藥, 2021, 52(10): 2927-2937.

    [18] 謝臻. 大承氣湯方藥物質(zhì)基礎(chǔ)及其配伍規(guī)律研究 [D].廣州: 廣州中醫(yī)藥大學(xué), 2009.

    [19] 王琦, 盛增秀. 經(jīng)方應(yīng)用 [M]. 銀川: 寧夏人民出版社, 1981: 263.

    [20] 張植軒. 小承氣湯證探究 [J]. 北京中醫(yī), 1994(4): 46-47.

    [21] 羅思妮, 彭致鋮, 范倩, 等. 經(jīng)典名方小承氣湯中化學(xué)成分的UPLC-Q-Orbitrap-MS分析 [J]. 中國(guó)實(shí)驗(yàn)方劑學(xué)雜志, 2021, 27(23): 1-10.

    [22] 施理畫. Box-Behnken響應(yīng)面法用于小承氣湯的煎煮工藝優(yōu)化 [J]. 中國(guó)藥師, 2018, 21(5): 813-817.

    [23] 陳澤麒, 謝欣媛, 陳思銘, 等. 基于數(shù)據(jù)驅(qū)動(dòng)和機(jī)制模型的丹參提取動(dòng)力學(xué)研究 [J]. 中草藥, 2022, 53(1): 51-59.

    [24] 黃興國(guó), 張靜, 杜文慧, 等. QbD理念的經(jīng)典名方研究: 四物湯加減方的提取工藝優(yōu)化研究 [J]. 中國(guó)中藥雜志, 2019, 44(20): 4329-4335.

    [25] 楊艷玲, 李花花, 黃嘉怡, 等. 基于質(zhì)量源于設(shè)計(jì)(QbD)理念的經(jīng)典名方桃紅四物湯的提取工藝研究 [J]. 中草藥, 2022, 53(2): 403-412.

    [26] 李銳華, 孫增濤, 王永香, 等. 運(yùn)用AHP結(jié)合熵權(quán)法優(yōu)化疏風(fēng)定喘顆粒炙麻黃組的醇沉工藝 [J]. 中國(guó)實(shí)驗(yàn)方劑學(xué)雜志, 2016, 22(20): 1-5.

    [27] 國(guó)家藥品監(jiān)督管理局藥品審評(píng)中心. 中藥新藥質(zhì)量標(biāo)準(zhǔn)研究技術(shù)指導(dǎo)原則(試行) [EB/OL]. (2020-10-12) [2021-11-04]. https://www. cde. org. cn/zdyz/ domesticinfopage?zdyzIdCODE=b30261779e4848579b15cf34608a06e7.html.

    [28] 王智民, 劉菊妍, 王德勤, 等. 關(guān)于經(jīng)典名方研發(fā)的一些重要關(guān)鍵信息和科學(xué)問(wèn)題的幾點(diǎn)看法 [J]. 中國(guó)實(shí)驗(yàn)方劑學(xué)雜志, 2022, 28(1): 212-217.

    Optimization of extraction process for classical prescription Xiaochengqi Decoction based on Box-Behnken experimental design, standard relation and analytic hierarchy process combined with entropy method

    ZHAO Yue-ying, WANG Chang-hai, ZHANG Ze-kang, ZHANG Qing, GUO Ming-xue, SUN Meng, XI Cheng, YANG Xue-ying, DONG Shuang, ZHANG Xin-yu, DU Shou-ying, LU Yang, BAI Jie

    Beijing University of Chinese Medicine, School of Chinese Materia Medica, Beijing 102488, China

    The extraction process of the classical prescription Xiaochengqi Decoction (小承氣湯, XD) was optimized based on Box-Behnken experimental design, standard relation and analytic hierarchy process combined with entropy method.With the contents of aloe emodin and chrysophanol, magnolol, synephrine as well as the paste-forming rate and fingerprint similarity as critical quality attributes (CQAs). The amount of water, extraction time and frequency were used as the critical process parameters (CPPs). The single factor experiment design was used to determine the level of each factor, and the Box-Behnken experimental design was used to optimize the extraction process. Finally, the standard relation between each samples and standard decoction of various indexes of different extraction parameters was calculated, while AHP-entropy method was applied to determine the weight of each index for comprehensive score.The results of variance analysis showed that thevalue of the quadratic regression model was significant (< 0.01), and the lack-of-fit value was not significant, indicating that the model was statistically significant. The relationship between various factors and the response values could be functionalized by the established models. Through further optimization of design space, the optimal extraction process of XD was 6-fold volume of water, extracted once and extracted for 30 min. The average comprehensive evaluation of three parallel verification tests was 95.37, RSD value was 2.31%, and meet the quality requirements of standard decoction.The extraction process parameters of XD based on the quality by design (QbD), standard relation and analytic hierarchy process combined with entropy method was stable and reliable, which would provide guidance for the process development of other classical prescriptions.

    quality by design; classical prescription; Xiaochengqi Decoction; Box-Behnken design-response surface method; standard relation; analytic hierarchy process-entropy weight method; aloe emodin; chrysophanol; magnolol; sinflin;fingerprint; critical quality attributes; critical process parameters

    R283.6

    A

    0253 - 2670(2022)20 - 6472 - 09

    10.7501/j.issn.0253-2670.2022.20.019

    2022-04-12

    國(guó)家自然科學(xué)基金資助項(xiàng)目(82173989)

    趙玥瑛(1998—),女,碩士研究生,從事中藥新劑型開發(fā)與研究。Tel: 18197341874 E-mail: ll2820064026@163.com

    陸 洋,男,教授,博士生導(dǎo)師,從事中藥新劑型新技術(shù)研究。Tel: (010)53912123 E-mail: landocean28@163.com

    白 潔,女,副教授,碩士生導(dǎo)師,從事中藥新制劑與新技術(shù)研究。Tel: (010)84738657 E-mail: baijie22811@163.com

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

    猜你喜歡
    膏率名方蘆薈
    川貝止咳膏的制備工藝優(yōu)化
    我的植物朋友——蘆薈
    四物湯不同制備工藝出膏率比較研究※
    經(jīng)典名方溫膽湯的研究進(jìn)展
    蕪菁膏超聲提取工藝優(yōu)化及其抗氧化活性研究
    我愛(ài)蘆薈
    小讀者(2021年6期)2021-07-22 01:50:02
    逍遙散系列經(jīng)典名方防治阿爾茨海默病的研究進(jìn)展
    抗衰老名方“金匱腎氣丸”的妙用
    瘟疫苦難中走出的名醫(yī)名方
    女貞子提取工藝優(yōu)化研究
    欧美亚洲日本最大视频资源| 另类精品久久| 美女国产高潮福利片在线看| 中国国产av一级| 国产成人免费观看mmmm| 男女啪啪激烈高潮av片| 波野结衣二区三区在线| 伦精品一区二区三区| 蜜桃在线观看..| 亚洲欧美成人精品一区二区| 国产伦理片在线播放av一区| a级片在线免费高清观看视频| 精品人妻在线不人妻| 日韩一区二区视频免费看| 久热久热在线精品观看| 久久久久久久大尺度免费视频| 一边摸一边做爽爽视频免费| 亚洲欧美清纯卡通| 最近手机中文字幕大全| 亚洲四区av| 一区二区三区乱码不卡18| 欧美日韩国产mv在线观看视频| 久久99热这里只频精品6学生| 国产高清国产精品国产三级| 精品一区在线观看国产| 精品久久国产蜜桃| 在线观看三级黄色| 亚洲av福利一区| 午夜久久久在线观看| 高清av免费在线| 九色成人免费人妻av| 亚洲综合色网址| 亚洲精品国产色婷婷电影| 日本猛色少妇xxxxx猛交久久| 香蕉国产在线看| 国产极品天堂在线| 91国产中文字幕| 日本色播在线视频| 免费在线观看黄色视频的| 26uuu在线亚洲综合色| 精品一区二区免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 女人精品久久久久毛片| 日日撸夜夜添| 免费av中文字幕在线| 免费av中文字幕在线| 国产激情久久老熟女| 丰满乱子伦码专区| 国产亚洲最大av| 久久ye,这里只有精品| 亚洲精品久久午夜乱码| 国产成人免费无遮挡视频| 黄色怎么调成土黄色| 高清不卡的av网站| 精品卡一卡二卡四卡免费| 亚洲av免费高清在线观看| 亚洲国产最新在线播放| 国产精品熟女久久久久浪| 9热在线视频观看99| 国产xxxxx性猛交| tube8黄色片| 日韩精品免费视频一区二区三区 | 成年人午夜在线观看视频| 国产精品熟女久久久久浪| 五月天丁香电影| 久久人人爽人人片av| 亚洲精品456在线播放app| 激情视频va一区二区三区| 亚洲美女搞黄在线观看| 伊人久久国产一区二区| 久久久a久久爽久久v久久| 国产伦理片在线播放av一区| 99国产综合亚洲精品| 国产黄色免费在线视频| 看非洲黑人一级黄片| 久久这里只有精品19| xxx大片免费视频| 多毛熟女@视频| 国产高清国产精品国产三级| 国产成人免费无遮挡视频| 国产精品熟女久久久久浪| 18禁动态无遮挡网站| 亚洲精品一区蜜桃| 日韩欧美一区视频在线观看| 久久久久久久久久成人| 女性生殖器流出的白浆| 亚洲av日韩在线播放| 日韩人妻精品一区2区三区| 亚洲美女黄色视频免费看| 日韩精品有码人妻一区| 国产成人免费观看mmmm| 色哟哟·www| 欧美性感艳星| 王馨瑶露胸无遮挡在线观看| 9热在线视频观看99| 女的被弄到高潮叫床怎么办| 精品人妻一区二区三区麻豆| 国产淫语在线视频| 午夜视频国产福利| 女性生殖器流出的白浆| 交换朋友夫妻互换小说| av播播在线观看一区| 婷婷色综合大香蕉| 国产精品国产三级国产av玫瑰| 亚洲欧洲精品一区二区精品久久久 | 91aial.com中文字幕在线观看| 成人国产麻豆网| 精品久久蜜臀av无| 精品第一国产精品| 婷婷色综合www| 麻豆乱淫一区二区| 香蕉精品网在线| 亚洲国产毛片av蜜桃av| 少妇被粗大的猛进出69影院 | 亚洲第一区二区三区不卡| 午夜福利网站1000一区二区三区| 五月伊人婷婷丁香| 精品久久国产蜜桃| 精品国产国语对白av| 久久久精品区二区三区| 久久久久久久国产电影| 精品亚洲成国产av| 久久这里有精品视频免费| 欧美激情国产日韩精品一区| 天堂中文最新版在线下载| 日本av免费视频播放| 大片电影免费在线观看免费| 中文字幕av电影在线播放| 我要看黄色一级片免费的| 亚洲欧美一区二区三区国产| 久久久久久人人人人人| 精品国产一区二区三区四区第35| 欧美激情国产日韩精品一区| 国产精品久久久久久久电影| 国产黄色免费在线视频| 天天躁夜夜躁狠狠久久av| 国产精品无大码| 亚洲av福利一区| 亚洲,一卡二卡三卡| 男人爽女人下面视频在线观看| 欧美日韩亚洲高清精品| 午夜影院在线不卡| 成年人免费黄色播放视频| 永久网站在线| 中文欧美无线码| 777米奇影视久久| av卡一久久| 三级国产精品片| 涩涩av久久男人的天堂| 捣出白浆h1v1| 日韩av不卡免费在线播放| 国产麻豆69| 中文字幕最新亚洲高清| 91精品国产国语对白视频| 熟女人妻精品中文字幕| 国产在线一区二区三区精| 97超碰精品成人国产| 中文字幕人妻丝袜制服| 波多野结衣一区麻豆| 成人国语在线视频| 欧美日韩视频高清一区二区三区二| 欧美国产精品va在线观看不卡| 中文字幕制服av| 国产成人午夜福利电影在线观看| 国产av码专区亚洲av| 欧美日韩国产mv在线观看视频| 韩国高清视频一区二区三区| 欧美97在线视频| 国产有黄有色有爽视频| 男女边摸边吃奶| 18禁裸乳无遮挡动漫免费视频| 亚洲,一卡二卡三卡| 欧美人与性动交α欧美软件 | 国产精品一区二区在线观看99| 1024视频免费在线观看| 美女国产视频在线观看| 国产亚洲午夜精品一区二区久久| 18禁在线无遮挡免费观看视频| 三级国产精品片| 精品视频人人做人人爽| 90打野战视频偷拍视频| 亚洲av国产av综合av卡| 各种免费的搞黄视频| 日韩成人av中文字幕在线观看| 春色校园在线视频观看| 国产精品久久久久久精品古装| 精品国产一区二区三区四区第35| 大香蕉久久网| 91午夜精品亚洲一区二区三区| 精品久久蜜臀av无| 亚洲高清免费不卡视频| 欧美成人午夜精品| 国产精品久久久av美女十八| 亚洲成色77777| 在线观看免费日韩欧美大片| 九色成人免费人妻av| 亚洲综合精品二区| 亚洲国产成人一精品久久久| 两个人免费观看高清视频| 大话2 男鬼变身卡| 亚洲国产精品专区欧美| 看免费av毛片| 99热6这里只有精品| 一个人免费看片子| 黄色一级大片看看| 男人操女人黄网站| 午夜av观看不卡| 天天影视国产精品| 精品国产国语对白av| 日本黄大片高清| 久久青草综合色| 中文天堂在线官网| 亚洲久久久国产精品| 久久国产精品大桥未久av| 26uuu在线亚洲综合色| 日本av手机在线免费观看| 亚洲五月色婷婷综合| 精品人妻在线不人妻| 22中文网久久字幕| 免费av中文字幕在线| 五月天丁香电影| 亚洲欧美色中文字幕在线| 免费观看在线日韩| 亚洲成国产人片在线观看| 我要看黄色一级片免费的| 国产色爽女视频免费观看| 香蕉丝袜av| 少妇人妻久久综合中文| 久久97久久精品| 亚洲国产欧美日韩在线播放| 亚洲久久久国产精品| 免费观看a级毛片全部| 成人18禁高潮啪啪吃奶动态图| 三级国产精品片| 国产伦理片在线播放av一区| 久久99精品国语久久久| 久久热在线av| 超碰97精品在线观看| 人人妻人人爽人人添夜夜欢视频| 91午夜精品亚洲一区二区三区| 国产精品三级大全| 狂野欧美激情性bbbbbb| 赤兔流量卡办理| av播播在线观看一区| 观看美女的网站| 女人被躁到高潮嗷嗷叫费观| 国产一区二区激情短视频 | 欧美丝袜亚洲另类| 国产又色又爽无遮挡免| 国产成人精品无人区| 久久这里只有精品19| 99久久综合免费| 亚洲精品国产av蜜桃| 色94色欧美一区二区| 亚洲国产精品专区欧美| 最后的刺客免费高清国语| av在线老鸭窝| 久久人妻熟女aⅴ| 国产精品三级大全| 久久精品夜色国产| a级片在线免费高清观看视频| 91国产中文字幕| av国产精品久久久久影院| 欧美成人午夜免费资源| 国产免费现黄频在线看| 最近最新中文字幕免费大全7| 春色校园在线视频观看| av黄色大香蕉| 制服丝袜香蕉在线| 激情五月婷婷亚洲| 久久久久国产精品人妻一区二区| 一区二区日韩欧美中文字幕 | 人人妻人人澡人人看| 9191精品国产免费久久| 少妇精品久久久久久久| 亚洲一码二码三码区别大吗| 国产亚洲欧美精品永久| 亚洲欧美一区二区三区国产| 午夜免费男女啪啪视频观看| 蜜臀久久99精品久久宅男| 日韩成人av中文字幕在线观看| 亚洲精品456在线播放app| 伦理电影大哥的女人| 亚洲内射少妇av| 国产精品嫩草影院av在线观看| 亚洲丝袜综合中文字幕| 汤姆久久久久久久影院中文字幕| 夫妻午夜视频| 久久精品久久久久久久性| 免费大片黄手机在线观看| 日本wwww免费看| 久久精品国产综合久久久 | 最新中文字幕久久久久| 色视频在线一区二区三区| 中文字幕亚洲精品专区| 亚洲成av片中文字幕在线观看 | 伦理电影大哥的女人| av片东京热男人的天堂| 夜夜骑夜夜射夜夜干| 一区二区日韩欧美中文字幕 | 精品第一国产精品| 在线亚洲精品国产二区图片欧美| 波野结衣二区三区在线| av又黄又爽大尺度在线免费看| 国产成人aa在线观看| 亚洲久久久国产精品| 菩萨蛮人人尽说江南好唐韦庄| 国产片特级美女逼逼视频| 国产伦理片在线播放av一区| 免费在线观看黄色视频的| 春色校园在线视频观看| 亚洲精品久久成人aⅴ小说| 国产一区二区三区av在线| 嫩草影院入口| 美女中出高潮动态图| 免费高清在线观看视频在线观看| 国产成人av激情在线播放| 99久久人妻综合| 高清黄色对白视频在线免费看| 国产精品一区二区在线不卡| 欧美成人午夜免费资源| 在线看a的网站| 制服人妻中文乱码| 三上悠亚av全集在线观看| 亚洲第一区二区三区不卡| www.熟女人妻精品国产 | 中文天堂在线官网| 亚洲情色 制服丝袜| 精品99又大又爽又粗少妇毛片| 国产男女超爽视频在线观看| 精品午夜福利在线看| 人妻系列 视频| 女人久久www免费人成看片| 亚洲综合色惰| 婷婷色综合www| 热99久久久久精品小说推荐| 欧美人与性动交α欧美精品济南到 | 尾随美女入室| 精品人妻偷拍中文字幕| 亚洲情色 制服丝袜| 老熟女久久久| 亚洲av欧美aⅴ国产| 黄片播放在线免费| 涩涩av久久男人的天堂| 1024视频免费在线观看| av国产精品久久久久影院| 一级爰片在线观看| av在线观看视频网站免费| 嫩草影院入口| 国产av一区二区精品久久| 国产一区二区在线观看日韩| 美女国产高潮福利片在线看| 国产一区二区三区综合在线观看 | 国产深夜福利视频在线观看| 肉色欧美久久久久久久蜜桃| 一级a做视频免费观看| 女人被躁到高潮嗷嗷叫费观| 不卡视频在线观看欧美| 日日撸夜夜添| 韩国av在线不卡| 久久女婷五月综合色啪小说| 大片电影免费在线观看免费| 国产无遮挡羞羞视频在线观看| 大陆偷拍与自拍| 成人漫画全彩无遮挡| 在线观看人妻少妇| av播播在线观看一区| 高清欧美精品videossex| 久久亚洲国产成人精品v| 亚洲国产av影院在线观看| 一二三四在线观看免费中文在 | 午夜日本视频在线| 久久 成人 亚洲| 午夜福利乱码中文字幕| 亚洲天堂av无毛| 99久国产av精品国产电影| 中国美白少妇内射xxxbb| 亚洲成国产人片在线观看| 蜜桃国产av成人99| 纵有疾风起免费观看全集完整版| 在线观看美女被高潮喷水网站| 精品一区二区三卡| 一区二区三区精品91| 一本大道久久a久久精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产伦理片在线播放av一区| www.色视频.com| 少妇熟女欧美另类| 欧美日韩综合久久久久久| 亚洲国产精品一区三区| tube8黄色片| 国产精品一国产av| 色婷婷久久久亚洲欧美| 女人久久www免费人成看片| 欧美精品av麻豆av| 欧美亚洲日本最大视频资源| 一级毛片黄色毛片免费观看视频| 亚洲欧洲精品一区二区精品久久久 | av在线app专区| 午夜日本视频在线| 午夜91福利影院| 亚洲久久久国产精品| 一级a做视频免费观看| 亚洲伊人色综图| 成人手机av| 丝袜脚勾引网站| 国产白丝娇喘喷水9色精品| 搡老乐熟女国产| 亚洲欧美色中文字幕在线| 少妇被粗大猛烈的视频| 水蜜桃什么品种好| 色94色欧美一区二区| av天堂久久9| 久久久a久久爽久久v久久| 哪个播放器可以免费观看大片| 午夜激情av网站| 99视频精品全部免费 在线| 男女高潮啪啪啪动态图| 丁香六月天网| 亚洲熟女精品中文字幕| 久久久久精品性色| 国产成人精品在线电影| 亚洲av电影在线观看一区二区三区| 国产色爽女视频免费观看| 久久精品久久久久久噜噜老黄| 韩国av在线不卡| 满18在线观看网站| 男男h啪啪无遮挡| 18+在线观看网站| av免费在线看不卡| 欧美日韩亚洲高清精品| 男女啪啪激烈高潮av片| 久久国内精品自在自线图片| 岛国毛片在线播放| 制服人妻中文乱码| 日本欧美视频一区| 99热网站在线观看| 丰满饥渴人妻一区二区三| 永久网站在线| 王馨瑶露胸无遮挡在线观看| 亚洲成av片中文字幕在线观看 | 91aial.com中文字幕在线观看| 丝袜脚勾引网站| 插逼视频在线观看| 日日爽夜夜爽网站| 午夜福利,免费看| 多毛熟女@视频| 国产成人精品一,二区| 精品一区二区三区四区五区乱码 | 欧美 亚洲 国产 日韩一| 伦理电影大哥的女人| 免费久久久久久久精品成人欧美视频 | 狠狠精品人妻久久久久久综合| 丝袜人妻中文字幕| 又黄又粗又硬又大视频| 国产一区二区激情短视频 | 成人亚洲欧美一区二区av| 天堂8中文在线网| 美女中出高潮动态图| 99久久中文字幕三级久久日本| 国产日韩欧美视频二区| 日本黄色日本黄色录像| 久久影院123| 考比视频在线观看| 国产亚洲欧美精品永久| 老女人水多毛片| 亚洲伊人久久精品综合| 久久久久久人人人人人| 欧美精品一区二区免费开放| 99热全是精品| 天天躁夜夜躁狠狠久久av| www日本在线高清视频| 免费女性裸体啪啪无遮挡网站| 亚洲欧美清纯卡通| 国产成人精品在线电影| 少妇猛男粗大的猛烈进出视频| 精品少妇黑人巨大在线播放| 久久ye,这里只有精品| 成年动漫av网址| 我的女老师完整版在线观看| 国产爽快片一区二区三区| 男人操女人黄网站| 丰满饥渴人妻一区二区三| 精品人妻偷拍中文字幕| 在现免费观看毛片| 肉色欧美久久久久久久蜜桃| 亚洲综合精品二区| 国产高清不卡午夜福利| 91精品伊人久久大香线蕉| 五月开心婷婷网| 黄色毛片三级朝国网站| 欧美激情极品国产一区二区三区 | 国产精品人妻久久久久久| 精品久久国产蜜桃| 新久久久久国产一级毛片| 精品酒店卫生间| 欧美日韩视频高清一区二区三区二| 热99久久久久精品小说推荐| 一级片'在线观看视频| 一级爰片在线观看| 精品午夜福利在线看| 日韩视频在线欧美| 亚洲成国产人片在线观看| av线在线观看网站| 日韩中字成人| 亚洲av综合色区一区| 1024视频免费在线观看| 在线免费观看不下载黄p国产| 男的添女的下面高潮视频| 国产 精品1| 久久免费观看电影| 亚洲欧美一区二区三区国产| 美女大奶头黄色视频| 下体分泌物呈黄色| 欧美精品人与动牲交sv欧美| 久久久久网色| 天天操日日干夜夜撸| 久久久久网色| 国产永久视频网站| 亚洲欧美中文字幕日韩二区| av网站免费在线观看视频| 国产精品久久久久久精品电影小说| 永久免费av网站大全| 国产在视频线精品| 黑丝袜美女国产一区| 国产免费一区二区三区四区乱码| 国产精品.久久久| 丝袜美足系列| 亚洲成人av在线免费| 99九九在线精品视频| 亚洲精品美女久久久久99蜜臀 | 一区二区三区四区激情视频| 在线天堂中文资源库| 91精品伊人久久大香线蕉| 亚洲人成77777在线视频| videossex国产| 性高湖久久久久久久久免费观看| 国产av国产精品国产| 黄色配什么色好看| 国产免费视频播放在线视频| 日本色播在线视频| 在线亚洲精品国产二区图片欧美| 亚洲av成人精品一二三区| 久久精品国产亚洲av天美| 美女国产视频在线观看| 精品亚洲成国产av| 大码成人一级视频| 男人添女人高潮全过程视频| 日本猛色少妇xxxxx猛交久久| 欧美97在线视频| 一级a做视频免费观看| 秋霞伦理黄片| 亚洲欧洲精品一区二区精品久久久 | 日韩电影二区| 色哟哟·www| 成年动漫av网址| 久久久精品94久久精品| 免费在线观看完整版高清| 久久青草综合色| 亚洲经典国产精华液单| av国产久精品久网站免费入址| 国产男女超爽视频在线观看| 久久国内精品自在自线图片| 交换朋友夫妻互换小说| 国产乱来视频区| 一区二区三区精品91| 久热久热在线精品观看| 国产麻豆69| 18禁观看日本| 十分钟在线观看高清视频www| 亚洲第一av免费看| 有码 亚洲区| 色5月婷婷丁香| 亚洲色图综合在线观看| 波野结衣二区三区在线| 国产精品久久久久成人av| 人成视频在线观看免费观看| 一边摸一边做爽爽视频免费| 亚洲天堂av无毛| 99久久人妻综合| 国产一区亚洲一区在线观看| 一区二区日韩欧美中文字幕 | 下体分泌物呈黄色| 蜜臀久久99精品久久宅男| 香蕉精品网在线| 精品一区二区三区四区五区乱码 | 国产精品蜜桃在线观看| 亚洲欧洲国产日韩| av黄色大香蕉| 国产亚洲欧美精品永久| 日韩成人伦理影院| 观看美女的网站| 人人妻人人澡人人爽人人夜夜| av不卡在线播放| freevideosex欧美| 999精品在线视频| av片东京热男人的天堂| 少妇精品久久久久久久| 男女边吃奶边做爰视频| 久久精品国产亚洲av涩爱| 毛片一级片免费看久久久久| 国产精品久久久久久久电影| 欧美精品亚洲一区二区| 免费高清在线观看视频在线观看| 午夜精品国产一区二区电影| 亚洲一级一片aⅴ在线观看| 赤兔流量卡办理| 久久97久久精品| 美女脱内裤让男人舔精品视频| 欧美日韩av久久| 欧美精品一区二区免费开放| av线在线观看网站| 日本猛色少妇xxxxx猛交久久| 深夜精品福利| 亚洲国产最新在线播放|