999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

星點設計-響應面法優化米拉貝隆緩釋片處方

2017-02-18 09:49:34浦雨偉丁逸梅盧定強南京工業大學藥學院南京20009江蘇省藥物研究所有限公司南京20009中國藥科大學藥學院南京20009
中國藥房 2017年1期
關鍵詞:優化

徐 斌,浦雨偉,周 衛,丁逸梅,2,盧定強,2(.南京工業大學藥學院,南京 20009;2.江蘇省藥物研究所有限公司,南京 20009;.中國藥科大學藥學院,南京 20009)

星點設計-響應面法優化米拉貝隆緩釋片處方

徐 斌1,2*,浦雨偉1,周 衛3,丁逸梅1,2,盧定強1,2(1.南京工業大學藥學院,南京 210009;2.江蘇省藥物研究所有限公司,南京 210009;3.中國藥科大學藥學院,南京 210009)

目的:優化米拉貝隆緩釋片處方。方法:采用聚氧乙烯(PEO)和羥丙甲基纖維素(HPMC K4M)為緩釋骨架材料,以粉末直壓法制備米拉貝隆緩釋片。以1、3、5、7 h的累積釋放度為指標,采用星點設計-響應面法優化處方中PEO、HPMC K4M和包衣液(OPADRY?)的用量,并進行驗證。比較所制緩釋片與原研片(MyrbetriqTM)在不同pH介質(水、pH 1.0的人工胃液、pH 4.5醋酸鹽緩沖液、pH 6.8磷酸鹽緩沖液)中不同轉速(100、50 r/min)下的累積釋放度,計算相似因子f2,擬合緩釋片的釋藥模型。結果:優化處方中每片米拉貝隆緩釋片含米拉貝隆25 mg、PEO 108.02 mg、HPMC K4M 21.69 mg、OPADRY?2.27%;其1、3、5、7 h的累積釋放度與預測值的相對誤差分別為4.78%、3.48%、0.69%、-1.41%。所制緩釋片與原研片在不同pH介質中和不同轉速下釋放曲線的f2均大于65;緩釋片釋藥符合零級釋藥特征。結論:優化處方制備的米拉貝隆緩釋片與原研片的釋放行為相似。

米拉貝??;緩釋片;星點設計-響應面法;處方優化;累積釋放度

米拉貝隆(Mirabegron)為白色粉末,易溶于甲醇、二甲基亞砜,幾乎不溶于水(0.082 mg/ml)。其屬于芳乙醇胺類新型選擇性β3腎上腺素受體激動藥,作用于膀胱逼尿平滑肌β3腎上腺素受體,使膀胱舒張,促進膀胱充盈并增加儲尿量,能有效減少排尿次數,改善膀胱活動過度引起的尿頻、尿急和尿失禁等,是此類中首個用于治療膀胱活動過度癥(OAB)的藥物。米拉貝隆口服50 mg/d后迅速吸收,tmax為3~5 h,cmax為24~26 ng/mL,終末半衰期約為50 h,口服生物利用度約為35%,效應呈劑量依賴性增加。該藥主要經水解酶代謝,部分經葡萄糖醛酸軛合反應及細胞色素P450(CYP)代謝,主要以原型經尿液和糞便排出體外,部分經呼吸排出[1]。米拉貝隆片劑由日本安斯泰來(Astellas)制藥公司開發,商品名MyrbetriqTM,上市規格為25 mg和50 mg,于2011年9月16日在日本上市,2012年6月28日經美國FDA批準用于治療成年人OAB。米拉貝隆為OAB患者提供了新的治療方案,且臨床藥效較好[2]。

米拉貝隆目前在國內尚無上市,對米拉貝隆的合成工藝研究已有文獻報道,但制劑研究報道較少。本研究采用星點設計-響應面法優化米拉貝隆緩釋片處方,以原研片劑MyrbetriqTM作為參比制劑,比較自制緩釋片與參比制劑體外釋放行為的相似性。

1 材料

1.1 儀器

C&C600B型單沖式壓片機(北京創博佳維科技有限公司);YD-20型片劑硬度儀、RC806型溶出試驗儀(天津天大天發科技有限公司);LC-2010AHT型高效液相色譜儀、2010A型紫外檢測器和島津工作站(日本Shimadzu公司)。

1.2 藥品與試劑

米拉貝隆對照品(批號:141101,純度:99.71%)及其原料藥(批號:150105,純度:>98%)均由南京工業大學藥學院提供;MyrbetriqTM(日本Astellas Pharma Technologies公司,批號:F1400110,規格:每片25 mg);聚氧乙烯(PEO,型號:WSR N80)、羥丙甲纖維素(HPMC K4M)、歐巴代胃溶型薄膜包衣液(OPADRY?)均購自上海卡樂康包衣技術有限公司;硬脂酸鎂(安徽山河藥用輔料有限公司,藥用級);磷酸、磷酸二氫鉀為分析純;乙腈、庚烷磺酸鈉為色譜純;水為去離子水。

2 方法與結果

2.1 米拉貝隆緩釋片的制備

按處方量稱取已過80目篩的米拉貝隆原料藥與HPMC K4M、PEO,采用等量遞增法混合均勻,加入適量助流劑硬脂酸鎂。按每片片質量約250 mg、含米拉貝隆25 mg的規格,采用橢圓形8 mm×16 mm淺凹沖模,粉末直接壓片,硬度控制在8~10 kg/cm2,以OPADRY?包衣后,即得。

2.2 米拉貝隆含量測定方法的建立

2.2.1 對照品溶液的制備 精密稱取米拉貝隆對照品10 mg,置于100 mL量瓶中,用甲醇溶解并稀釋定容;精密量取1.0 mL,置于10 mL量瓶中,用釋放介質稀釋定容至刻度,即得。

2.2.2 色譜條件 色譜柱:Cosmosil C18(250 mm×4.6 mm,5 μm);流動相:10 mmol/L庚烷磺酸鈉+20 mmol/L磷酸二氫鉀溶液(用磷酸調節pH至3.0)-乙腈(75∶25,V/V),流速:1.0 mL/min;檢測波長:250 nm;柱溫:35℃;進樣量:20 μL。

2.2.3 方法學考察 按相關方法進行檢測,結果顯示,所用輔料均對米拉貝隆測定無干擾,其出峰時間為7.4 min。拉米貝隆的檢測質量濃度的線性范圍為2.818~36.63 μg/mL,峰面積(A)對質量濃度(c)的回歸方程為A=6.619 4c+1.682×103(r=0.999 9);精密度試驗中RSD=0.07%(n=6);穩定性試驗中樣品溶液室溫放置24 h的峰面積RSD=0.25%(n=6);低、中、高質量濃度(10.02、25.05、32.56 μg/mL)的質控溶液的平均回收率分別為99.78%、100.6%、98.97%,RSD分別為0.87%、1.14%、1.07%(n=3),均符合相關要求。

2.3 釋放度的測定方法

按2015年版《中國藥典》(四部)通則0931溶出度與釋放度測定第一法[3],以脫氣的pH 6.8磷酸鹽緩沖液(PBS)900 mL為釋放介質,轉速為100 r/min,溫度為(37±0.5)℃,分別于1、3、5、7、8.5、10、12 h取樣10 mL,同時補充等量等溫釋放介質。取樣液經0.45 μm微孔濾膜過濾,濾液按“2.2”項下方法測定米拉貝隆含量,計算累積釋放度(Q)。

2.4 星點設計-響應面法優化處方

2.4.1 釋放因素考察 本研究采用單因素試驗對處方各因素進行了預試驗。結果發現,米拉貝隆緩釋片體外釋放行為受PEO、HPMC和包衣液(OPADRY?)用量影響較大。根據原研片Q7h>90%,確定本研究以Q1h、Q3h、Q5h、Q7h為指標,將自制緩釋片與參比制劑(MyrbetriqTM)的體外釋藥行為采用相似因子法(f2)進行比較。f2=50·其中,Rt和Tt分別代表自制緩釋片和參比制劑第t點的平均累積釋放度,n為測試點數。當f2>50時認為二者釋放行為相似,f2的數值越大,釋放行為相似度越高[4]。

2.4.2 星點設計 采用星點設計,以X1(PEO,mg)、X2(HPMC K4M,mg)、X3(包衣液,%)的用量進行3因素5水平的優化試驗,在前期單因素基礎上確定每片用量優化范圍X1:70~140 mg,X2:10~40 mg,X3:1%~5%;以Q1h、Q3h、Q5h、Q7h為指標,采用Design Expert 8.0.6軟件進行分析。星點設計因素與水平見表1,星點設計與結果見表2。

表1 星點設計因素與水平Tab 1 Factors and levels of central composite design

2.4.3 模型擬合 采用Design Expert 8.0.6軟件對試驗數據進行多元線性回歸,并以方程復相關系數(R2)和置信度(P)為標準進行判定。擬合結果顯示,二次多項式模型為最佳擬合模型,各方程如下:Q1h=46.474 89-0.465 73X1+0.017 019X2-3.084 70X3-0.000 188 552X1X2+0.013 427X1X3-0.001 648 64X2X3+0.001 555 3X12-0.001 266 22X22+0.128 59X32(R2=0.968,P=0.006 1);Q3h=84.043 76-0.799 69X1+1.288 41X2-3.565 54X3-0.002 841 75X1X2+0.081 469X1X3- 0.146 26X2X3+ 0.001 805 61X12-0.016 567X22-0.287 96X32(R2=0.970,P=0.005 8);Q5h=112.022 76-0.034 737X1-2.296 55X2+ 1.748 91X3-0.000 781 144X1X2+0.007 066 74X1X3-0.068 536X2X3-0.000 099 724 5X12+0.023 519X22-0.037 987X32(R2=0.992,P=0.003 3);Q7h=55.651 21+ 0.324 72X1+2.283 06X2+6.543 82X3-0.009 791 24X1X2-0.008 076 27X1X3-0.128 83X2X3-0.000 059 934 9X12-0.043 231X22-0.432 95X32(R2=0.967,P=0.004 8)。由優化數學模型的參數可知,3個因素中,釋放前、中期(Q1h、Q3h)X1和X3的系數均為負,說明二者對釋放度的影響是同向的,在釋放前期(Q1h)X1和X3對釋放度影響極顯著(P<0.01),X2對釋放度影響較小;釋放中期(Q3h)X1和X2對釋放度影響極顯著(P<0.01),X3對釋放度幾乎未見影響(P>0.05);釋放中后期(Q5h)X2對釋放度影響極顯著(P<0.01),X1影響較小;釋放后期(Q7h)除X2對釋放度有極顯著影響外(P<0.01),X1和X3對釋放度幾乎未見影響(P>0.05)。同時R2值(均大于0.9)和P值(均小于0.05)亦表明模型的擬合程度良好、代表性好,可用來分析和預測米拉貝隆緩釋片的處方。

表2 星點設計與結果Tab 2 Design and results of central composite

2.4.4 響應面分析 Q1h、Q3h、Q5h、Q7h與3個因素的響應面圖見圖1(無明顯影響的圖略)。

由1 h釋放情況可知,當X2不變時,響應面隨X1和X3的增加而下降,而X2變化對釋放速率影響不顯著。同理分析可知,3 h和5 h的釋放速率受X1和X2相互影響,要使Q達到最佳值,需要將X1和X2優化到合適的比例。7h響應面與5 h響應面相比,區分度明顯增強,說明7 h骨架材料的控釋能力有所減弱,可能是因為親水凝膠X1已基本全部溶蝕,整個骨架已變得松散。

圖1 Q1h、Q3h、Q5h、Q7h與3個因素的響應面圖Fig 1 Response surface plot of Q1h,Q3h,Q5h,Q7hto 3 factors

2.4.5 處方優化與驗證 為驗證所建立模型的預測能力,將MyrbetriqTM的Q1h、Q3h、Q5h、Q7h作為目標值,經Design Expert 8.0.6軟件處理,優選處方區域為每片米拉貝隆緩釋片中X1為91.53~108.02 mg,X2為21.50~23.06 mg,X3為2.25%~4.19%。對預測結果進行驗證,在Design Expert 8.0.6軟件設計的輔料配比中,X2、X3變化范圍較小,因此選擇對X1的高、中、低3個用量處方進行驗證。1號、2號、3號處方每片輔料(PEO、HPMC K4M、包衣液)配比依次分別為:108.02 mg、21.69 mg、2.27%,95.36 mg、22.29 mg、3.68%,91.53 mg、23.06 mg、4.19%。首先將PEO、HPMC K4M和主藥充分混合均勻,再加適量硬脂酸鎂混勻,粉末直接壓片,硬度8~10 kg,包衣,測定釋放度。采用f2法比較測定值與預測值的相似性,結果見表3。

表3 處方優化結果Tab 3 Results of validation test

由表3可知,f2較大的1號處方為米拉貝隆緩釋片最優處方:每片含米拉貝隆25 mg、PEO 108.02 mg、HPMC K4M 21.69 mg、OPADRY?2.27%。對3批自制緩釋片在不同介質(水、pH 1.0的人工胃液、pH 4.5醋酸鹽緩沖液、pH 6.8磷酸鹽緩沖液)中不同轉速(100、50 r/min)下的累積釋放度進行測定,取均值與MyrbetriqTM進行比較。結果顯示,在不同介質和不同轉速下自制緩釋片與市售MyrbetriqTM的f2均大于65,表明自制樣品的釋藥行為與MyrbetriqTM基本一致。2種樣品在不同pH介質中不同轉速下的釋放比較見表4。

表4 2種樣品在不同pH介質中不同轉速下的釋放比較(n=6)Tab 4 Comparison of the release of 2 samples in different pH mediums at different rotation rates(n=6)

2.5 釋藥機制

藥物釋放度曲線擬合的模型主要有4種[5]:零級釋放模型、一級釋放模型、Higuchi方程、Ritger-pappas方程。對米拉貝隆緩釋片的體外釋藥數據(Q1h、Q3h、Q5h、Q7h)進行擬合,結果見表5。

由表6顯示,米拉貝隆緩釋片按Ritger-Pappas方程擬合的r最大。Ritger-Pappas模型中的n描述了釋放機制:當n≤0.45時,藥物主要以Fick擴散機制釋放;當0.45<n<0.89時,藥物以擴散和溶蝕并存的方式釋放;當n>0.89時,藥物主要以溶蝕方式釋放。同時,n也可以反映藥物釋放動力學方面的情況:當n>0.66時,藥物以零級動力學釋放為主[6]。表5中n=1.140 8>0.89,表明米拉貝隆緩釋片釋藥機制主要以溶蝕方式釋放,呈現零級釋藥特征。

表5 米拉貝隆緩釋片的釋藥模型擬合結果Tab 5 Fitting result of drug release model of Mirabegron sustained-release tablet

3 討論

3.1 骨架材料的選擇

市售片的骨架材料主要以PEO為主,PEO作為親水凝膠骨架材料具有優良的流動性和可壓性,按分子質量大小有多種規格。通常對于水溶性較大的藥物,其釋藥方式主要以凝膠擴散機制為主;溶解度較小藥物,其釋藥方式主要依靠骨架的膨脹和溶蝕。米拉貝隆為難溶性藥物,在水中幾乎不溶,因此選擇低分子質量PEO為主要骨架材料,但單獨應用PEO釋藥行為難以達到接近零級,因此加入適量HPMC K4M聯合使用[7]。PEO和HPMC K4M在吸水性、膨脹性、溶蝕性和凝膠強度各項性質參數上差異較大,本研究通過星點設計效應面法優化調整處方中的PEO、HPMC K4M的用量比例,使藥物避免出現釋藥前期釋放較快的現象,在控制骨架片初始釋藥量的同時又能保證其最終釋藥完全,使骨架片的體外釋放更加平穩[8]。

3.2 處方工藝的優化

米拉貝隆親水凝膠骨架片中PEO用量較大,如果濕法制粒,制得軟材會比較黏連、易結塊,制粒時還會產生大量較大的硬顆粒,難以過篩,不便于操作。本研究選取的PEO為顆粒型,可壓性和流動性優良,可采用粉末直接壓片。這樣既有利于操作,又可以避免濕法制粒在干燥時可能產生的二次污染。PEO有一定的引濕性,市售片通過包薄膜衣保證了其穩定性,根據其薄膜衣內容,本研究選取與其相似的OPADRY?進行包衣,并通過星點設計考察了其用量對前期釋放的影響,獲得一個優化水平范圍。

3.3 釋藥行為與釋藥機制

低分子質量PEO水化和溶蝕速度快,膨脹性能優良。處方中其他條件不變時,PEO的用量對藥物的釋放有顯著影響,PEO用量越小,前期藥物釋放越快[9]。將PEO與高黏度的HPMC K4M合用既可以使阻滯釋放的凝膠層迅速形成,又可使擴散前沿和溶蝕前沿同時向片心移動,兩者的相互作用調節了藥物的釋放速度,避免了突釋。在釋放中期吸水飽和的PEO水凝膠與HPMC K4M一起形成相對穩定的水凝膠層,藥物的釋放隨著骨架的逐漸溶蝕而進行;隨著藥物濃度下降以及PEO水凝膠的逐漸溶蝕,保證了藥物的平穩釋放。在釋藥后期藥物濃度顯著下降,由于飽和的PEO水凝膠大部分已經溶蝕,加快了藥物的釋放速率,此時藥物的釋放主要受HPMC K4M控制,HPMC K4M用量對后期釋放的影響起主導作用。在整個藥物釋放過程中,PEO與HPMC K4M兩者的相互作用,避免了在釋放前期發生突釋,而溶蝕型骨架又使得藥物在后期可完全釋放[10]。Ritger-Pappas方程的擬合結果(r=0.997 9,n=1.140 8)亦表明,米拉貝隆緩釋片釋藥主要為溶蝕,釋藥行為呈現零級釋藥特征。

本文通過星點設計-響應面法設計安排試驗并對結果進行擬合分析,建立了PEO、HPMC K4M和包衣液用量與緩釋片在Q1h、Q3h、Q5h、Q7h之間的多元線性模型。經試驗驗證,優化處方為每片米拉貝隆25 mg、PEO 108.02 mg、HPMC K4M 21.69 mg、OPADRY?2.27%,其釋放度的測定值與預測值趨勢基本一致,所建立的模型有良好的預測性,在不同介質及轉速條件下與原研片釋放行為相似。為了便于操作,建議將OPADRY?包衣液增質量擬定在2.0%~3.0%范圍內。

[1] 李海元,周峰海,馬玉磊,等.膀胱過度活動癥治療新藥米拉貝隆[J].中國新藥雜志,2014,23(19):2215-2218.

[2] 唐啟東,宮平.米拉貝隆[J].中國藥物化學雜志,2012,22(6):544.

[3] 國家藥典委員會.中華人民共和國藥典:四部[S].2015年版.北京:中國醫藥科技出版社,2015:121-124.

[4] 牛劍釗,林蘭,張啟明.美國和日本溶出曲線相似性判定方法介紹[J].中國藥物評價,2013,30(2):68-69.

[5] 周衛,謝浩,嚴相平,等.阿昔洛韋緩釋片在不同介質中的釋放研究[J].中國藥房,2011,22(21):2002-2004.

[6] 陳秀毅,石海英,趙愛麗.甲磺酸二氫麥角堿緩釋片的制備和體外評價[J].中國新藥雜志,2013,22(12):1444-1447.

[7] 王福堅,張自強,周建平,等.Box-Behnken效應面法篩選聚氧乙烯骨架緩釋片處方和釋藥機制的研究[J].藥學與臨床研究,2013,21(4):339-343.

[8] 聶淑芳,劉輝,劉艷麗,等.新型緩釋骨架材料聚氧化乙烯與HPMC凝膠層結構的差異比較[J].中國醫藥工業雜志,2011,42(7):521-524.

[9] 梁靜閣,潘衛三,林威,等.星點設計-效應面法優選銀杏內酯滲透泵控釋片的制備工藝與處方[J].中國藥房,2014,25(27):2516-2520.

[10] 王博,張來華,李苑新,等.親水凝膠骨架緩釋片釋藥機制評價方法的研究進展[J].中國醫藥工業雜志,2009,40(10):782-785.

(編輯:鄒麗娟)

Optimization of Mirabegron Sustained-release Tablets by Central Composite Design-response Surface Method

XU Bin1,2,PU Yuwei1,ZHOU Wei3,DING Yimei1,2,LU Dingqiang1,2(1.School of Pharmacy,Nanjing Tech University,Nanjing 210009,China;2.Jiangsu Provincial Institute of Materia Medica,Nanjing 210009,China;3. School of Pharmacy,China Pharmaceutical University,Nanjing 210009,China)

OBJECTIVE:To optimize the formulation of Mirabegron sustained-release tablets.METHODS:Using polyethylene oxide(PEO)and hydroxypropylmethyl cellulose(HPMC K4M)as the sustained release matrix,Mirabegron sustained-release tablets were prepared by powder direct compression technology.Using 1,3,5,7 h accumulative release rate as indexes,the amounts of PEO,HPMC K4M and OPADRY?were optimized by composite design-response surface method,and then validated.Accumulative release rates of sustained-release tablet and original tablet(MyrbetriqTM)were compared in different pH mediums(water,pH 1.0 simulated gastric fluid,pH 4.5 acetate buffer solution,pH 6.8 phosphate buffer solution)at different rotation rates(100,50 r/min),and similiar factor f2was calculated to fit drug release model of sustained-release tablet.RESULTS:In the optimized firmulation each Mirabegron sustained-release tablet contained mirabegron 25 mg,PEO 108.02 mg,HPMC K4M 21.69 mg,OPADRY?2.27%.Relative error of accumulative release rates at 1,3,5,7 h to predicted value were 4.78%,3.48%,0.69%and-1.41%,respectively.f2of release curves of sustained-release tablet and original tablet were higher than 65 in different pH medium at different rotation rates.The drug release of sustained-release tablet was fitted to zero-order release equation.CONCLUSIONS:Mirabegron sustained-release tablet by optimized technology is similar to original tablet in drug release behavior.

Mirabegron;Sustained-release tablet;Central composite design-response surface method;Formulation optimization;Accumulative release rate

R943

A

1001-0408(2017)01-0110-05

2016-07-06

2016-09-02)

*主管藥師,碩士。研究方向:藥物新制劑及新技術。電話:025-58139391。E-mail:nanjing5696@126.com

DOI10.6039/j.issn.1001-0408.2017.01.29

猜你喜歡
優化
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
PEMFC流道的多目標優化
能源工程(2022年1期)2022-03-29 01:06:28
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
由“形”啟“數”優化運算——以2021年解析幾何高考題為例
圍繞“地、業、人”優化產業扶貧
今日農業(2020年16期)2020-12-14 15:04:59
事業單位中固定資產會計處理的優化
消費導刊(2018年8期)2018-05-25 13:20:08
4K HDR性能大幅度優化 JVC DLA-X8 18 BC
幾種常見的負載均衡算法的優化
電子制作(2017年20期)2017-04-26 06:57:45
主站蜘蛛池模板: 老司机aⅴ在线精品导航| 伊人久综合| 日韩在线中文| 天天综合网色| 99久久性生片| 丁香婷婷激情网| 中美日韩在线网免费毛片视频| 999精品色在线观看| 亚洲国产精品久久久久秋霞影院 | 国产欧美日韩va另类在线播放| 国产网站一区二区三区| 高清无码一本到东京热| 制服丝袜国产精品| 在线网站18禁| 国产麻豆精品在线观看| 午夜福利免费视频| 国产成人永久免费视频| 国产福利2021最新在线观看| 日本日韩欧美| 亚洲天堂免费观看| 亚洲一区波多野结衣二区三区| 欧美成人综合视频| 一级全黄毛片| 青草精品视频| 无套av在线| 欧美日韩国产在线播放| 欧美国产在线看| 欧美va亚洲va香蕉在线| 91高清在线视频| 国产精品永久在线| 黄色网址免费在线| 国产精品密蕾丝视频| 国产乱子伦视频在线播放| 国产精品午夜福利麻豆| 欧美成人免费一区在线播放| 午夜福利亚洲精品| 91精品国产自产在线老师啪l| 91美女视频在线观看| 久久女人网| 亚洲爱婷婷色69堂| 久久这里只有精品国产99| 欧美a级在线| 日韩欧美中文在线| 一级香蕉人体视频| 97av视频在线观看| 波多野结衣视频一区二区 | 在线另类稀缺国产呦| 久久精品国产国语对白| 亚洲v日韩v欧美在线观看| 91久久偷偷做嫩草影院精品| 国产一区二区精品高清在线观看| 99在线观看国产| 538国产在线| 91啦中文字幕| 国产免费黄| 91成人在线观看视频| 日韩黄色在线| 午夜不卡视频| 国产尤物在线播放| 特黄日韩免费一区二区三区| 免费无码AV片在线观看中文| 亚洲无码91视频| 欧美午夜理伦三级在线观看| 亚洲免费人成影院| 再看日本中文字幕在线观看| 啪啪免费视频一区二区| 欧美三级不卡在线观看视频| 成人中文在线| 在线欧美日韩国产| 福利在线免费视频| 国产精品黄色片| 综合天天色| 国产精品99r8在线观看| 亚洲欧美在线综合图区| 色婷婷丁香| 精品久久高清| 国产乱子伦无码精品小说 | 午夜毛片免费看| 嫩草影院在线观看精品视频| 国产精品天干天干在线观看 | 91精品国产自产91精品资源| 亚洲一级毛片在线观播放|