亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于響應(yīng)面法優(yōu)化α-淀粉酶酶解鷹嘴豆淀粉工藝及建模的研究

        2015-08-02 03:58:37高義霞王淑芳王丹丹王廷璞
        食品工業(yè)科技 2015年15期
        關(guān)鍵詞:影響

        高義霞,王淑芳,王丹丹,王廷璞

        (天水師范學(xué)院生物工程與技術(shù)學(xué)院,甘肅天水 741001)

        基于響應(yīng)面法優(yōu)化α-淀粉酶酶解鷹嘴豆淀粉工藝及建模的研究

        高義霞,王淑芳,王丹丹,王廷璞*

        (天水師范學(xué)院生物工程與技術(shù)學(xué)院,甘肅天水 741001)

        以葡萄糖值(DE)為考察指標(biāo),利用單因素實(shí)驗(yàn)和響應(yīng)面法研究溫度、pH及酶與底物比(E/S)對(duì)α-淀粉酶酶解鷹嘴豆淀粉制糖的影響。應(yīng)用Lineweaver-Burk和Wilkinson統(tǒng)計(jì)法分別求解米氏常數(shù)(Km)和最大反應(yīng)速度(Vm),并建立相應(yīng)的動(dòng)力學(xué)模型。結(jié)果表明:α-淀粉酶酶解鷹嘴豆淀粉制糖的最佳參數(shù)為:pH 6.5,溫度55 ℃及E/S=0.12 U/mg,驗(yàn)證值為23.103%±0.454%。在30~50 ℃范圍內(nèi),Ea=18.977 kJ/moL,ΔH=19.624 kJ/mol。該研究為鷹嘴豆淀粉工業(yè)制糖提供理論依據(jù)。

        鷹嘴豆,酶解,動(dòng)力學(xué),響應(yīng)面

        鷹嘴豆(CicerarietinumL.),又稱桃豆、雞豌豆等,豆科鷹嘴豆屬植物[1],淀粉含量高達(dá)40%~60%[2],具有降血糖、降血脂[3]和抗腫瘤[4]等活性。鷹嘴豆淀粉消化率低、酶解過(guò)程與淀粉濃度、pH、溫度、酶類(lèi)型及離子強(qiáng)度等有關(guān),常規(guī)米氏方程難以真實(shí)描述反應(yīng)過(guò)程[5],而實(shí)際生產(chǎn)中,須明確酶促反應(yīng)條件及反應(yīng)動(dòng)力學(xué)特性,以最大化發(fā)揮酶催化作用。隨生淀粉酶等新型淀粉酶的出現(xiàn),直接選擇淀粉顆粒作為底物進(jìn)行水解成為一種趨勢(shì)。

        酶促反應(yīng)動(dòng)力學(xué)是研究酶促反應(yīng)速度及其影響因素的學(xué)科,在酶的結(jié)構(gòu)、功能、作用機(jī)制、反應(yīng)條件及代謝中作用等研究中,具有重要的理論和實(shí)踐意義。鐘昔陽(yáng)[6]利用α-淀粉酶對(duì)小麥面粉進(jìn)行了動(dòng)力學(xué)研究,結(jié)果表明,Wilkinson法預(yù)測(cè)值與實(shí)測(cè)值基本吻合。張國(guó)權(quán)[7]利用真菌淀粉酶對(duì)蕎麥淀粉進(jìn)行了動(dòng)力學(xué)研究,并計(jì)算了最大反應(yīng)速度(Vmax)和米氏常數(shù)(Km)。張麗華[8]利用響應(yīng)面法優(yōu)化蓮藕淀粉酶解工藝,得到最優(yōu)的酶量、溫度和時(shí)間。姜紹通[9]利用響應(yīng)面法對(duì)芋頭淀粉堿性蛋白酶酶解工藝參數(shù)進(jìn)行優(yōu)化,獲得最優(yōu)的時(shí)間、溫度、pH和酶量工藝參數(shù)。

        本實(shí)驗(yàn)以溫度、pH及酶與底物比(E/S)為因素,通過(guò)響應(yīng)面法優(yōu)化α-淀粉酶酶解鷹嘴豆淀粉制糖的最佳工藝,同時(shí)對(duì)酶解動(dòng)力學(xué)模型進(jìn)行研究,為鷹嘴豆淀粉工業(yè)制糖提供依據(jù)。

        1 材料與方法

        1.1 材料與儀器

        鷹嘴豆 淀粉含量為87%,新疆天山奇豆生物科技有限責(zé)任公司;α-淀粉酶(58.7 U/mg) Sigma產(chǎn)品(Lot:STBB0510,來(lái)自枯草芽孢桿菌);酒石酸鉀鈉、3,5-二硝基水楊酸、亞硫酸氫鈉、苯酚、葡萄糖等試劑為國(guó)產(chǎn)分析純 購(gòu)自生工生物工程股份有限公司。

        AL204型電子天平 梅特勒-托利多儀器有限公司;TGL-20M型高速臺(tái)式冷凍離心機(jī) 湖南湘儀離心機(jī)儀器有限公司;PHS-3D型 雷磁pH計(jì):上海精密科學(xué)有限公司;722型可見(jiàn)分光光度計(jì) 上海欣茂有限公司。

        1.2 實(shí)驗(yàn)方法

        1.2.1 鷹嘴豆淀粉的酶解 糖含量測(cè)定參考3,5-二硝基水楊酸(DNS)法[10]。稱取鷹嘴豆淀粉0.1 g,加入20 mL 0.05 mol/L pH 6.0磷酸緩沖液,攪拌,80 ℃糊化5 min,冷卻至50℃,按E/S=0.01 U/mg加入α-淀粉酶,振蕩混勻,50 min后迅速沸水滅活5 min,測(cè)定吸光度,計(jì)算葡萄糖值(DE)。

        DE(%)=還原糖量(mg)×100/[87%×淀粉量(mg)]

        式(1)

        式(2)

        式中:v為反應(yīng)速度(mg/(mL·min),m/V為還原糖含量(mg/mL);t為反應(yīng)時(shí)間(min)。

        1.2.2 單因素實(shí)驗(yàn) 在溫度50 ℃,pH6.0,E/S 0.04 U/mg條件下,考察時(shí)間分別為5、10、20、30、40、50、70,90、110、140、170、200 min時(shí)其對(duì)DE的影響;在pH6.0,E/S 0.04 U/mg,時(shí)間50 min條件下,考察溫度分別為30、40、50、60、70 ℃時(shí)其對(duì)DE的影響;在溫度50 ℃,E/S 0.04 U/mg,時(shí)間50 min條件下,考察pH分別為4、5、6、7、8時(shí)其對(duì)DE的影響;在溫度50 ℃、時(shí)間50 min、pH7.0條件下,考察E/S分別為0.01、0.02、0.04、0.08、0.16 U/mg時(shí)其對(duì)DE的影響。

        1.2.3 酶促反應(yīng)動(dòng)力學(xué)建模 Michaelis和Menten根據(jù)快速平衡學(xué)說(shuō)建立米氏方程,Briggs和Haldane用穩(wěn)態(tài)代替快速平衡態(tài),對(duì)米氏方程做了修正,得到更合理的酶促反應(yīng)動(dòng)力學(xué)過(guò)程。推導(dǎo)過(guò)程簡(jiǎn)述如下[11]:

        式(3)

        式(4)

        式(5)

        vm=k3E總

        式(6)

        式(7)

        式(8)

        式中,k1為ES生成速度;k2、k3分別為ES分解為E+S和P+E的速度;k4為E+P生成ES的速度(反映初期忽略不計(jì))。[S]為底物濃度(mg/mL);E總(U/mL)=[E]+[ES]。利用(8)式求Km和Vm[12]。

        1.2.4 酶量對(duì)反應(yīng)速度的影響 在反應(yīng)時(shí)間7 min、pH 7.0、底物濃度20 mg/mL條件下,分別加入0.2、0.4、0.6、0.8、1.0、1.2 mL U/mL酶,測(cè)定酶量對(duì)反應(yīng)速度的影響。

        1.2.5 響應(yīng)面法優(yōu)化酶解工藝 參考李忠海等[13],以pH、溫度及E/S為因子,各因子最優(yōu)點(diǎn)為基準(zhǔn)點(diǎn),確定適度步長(zhǎng),DE為響應(yīng)值,利用Design-Expert 8.06軟件設(shè)計(jì)Box-Benhnken實(shí)驗(yàn),實(shí)驗(yàn)見(jiàn)表1。在最優(yōu)工藝條件pH6.5,溫度55 ℃及E/S 0.12 U/mg進(jìn)行驗(yàn)證實(shí)驗(yàn)。

        表1 Box-Benhnken實(shí)驗(yàn)因素水平及編碼

        1.3 數(shù)據(jù)處理

        2 結(jié)果與分析

        2.1 標(biāo)準(zhǔn)曲線的制作

        DNS法測(cè)定糖含量的標(biāo)準(zhǔn)曲線方程為Y=1.629X-0.042,r=0.997(X為葡萄糖質(zhì)量,mg,Y為吸光度,A)。

        2.2 單因素實(shí)驗(yàn)

        時(shí)間對(duì)DE的影響見(jiàn)圖1,在0~50 min范圍內(nèi),隨時(shí)間的增加,DE逐漸增大,隨后出現(xiàn)下降,可能是復(fù)合糖生成速度逐漸接近、甚至超過(guò)淀粉水解速度,即存在較弱的復(fù)合反應(yīng),故DE的生成將趨于平穩(wěn)[14],據(jù)此,后續(xù)研究中時(shí)間選擇50 min。

        圖1 時(shí)間對(duì)DE的影響Fig.1 Effect of time on DE

        溫度對(duì)DE的影響見(jiàn)圖2,在30~50 ℃范圍內(nèi),隨溫度升高,酶與底物間碰撞次數(shù)增加,DE增加,50 ℃時(shí)達(dá)最大值,為最適溫度;另一方面,當(dāng)溫度增至50 ℃以上時(shí),酶逐漸失活,DE下降。

        圖2 溫度對(duì)DE的影響Fig.2 Effect of temperature on DE

        pH對(duì)DE的影響見(jiàn)圖3,pH7.0時(shí),DE達(dá)到最大值,為最適pH。pH影響酶活性中心必需基團(tuán)和底物的解離程度,從而影響酶與底物的結(jié)合,并最終影響DE的生成[15]。

        圖3 pH對(duì)DE的影響Fig.3 Effect of pH on DE

        E/S對(duì)DE的影響見(jiàn)圖4,隨E/S的增加,酶解速度幾乎呈線性增長(zhǎng),原因是反應(yīng)初期底物過(guò)量,酶促反應(yīng)速度與酶量成正比;隨E/S的進(jìn)一步增大,底物逐漸不足,反應(yīng)速度趨于平穩(wěn)[16]。因此E/S選擇0.1 U/mg。

        圖4 E/S對(duì)DE的影響Fig.4 Effect of E/S on DE

        2.3 響應(yīng)面分析

        Y=-140.651+30.884X1+2.233X2+115.114X3-0.367X1X2-50.268X1X3+3.439X2X3-0.821X12+2.436X22-610.919X32。

        Box-Benhnken結(jié)果及方差分析見(jiàn)表2和表3,方程具有顯著性(p=0.0263),表明模型顯著。失擬項(xiàng)p=0.213>0.05,失擬不顯著,回歸模型與實(shí)際情況擬合較好。pH極顯著,E/S影響顯著,溫度影響不顯著。各因素影響大小依次為pH、E/S及溫度。交互項(xiàng)影響均不顯著。最優(yōu)工藝為:pH 6.5,溫度55 ℃及E/S 0.12 U/mg。在此條件下進(jìn)行驗(yàn)證,結(jié)果見(jiàn)表4。

        表2 Box-Benhnken實(shí)驗(yàn)設(shè)計(jì)表及結(jié)果

        表3 回歸模型系數(shù)顯著性檢驗(yàn)及其方差分析

        表4 驗(yàn)證實(shí)驗(yàn)

        驗(yàn)證實(shí)驗(yàn)平均值為(23.103%±0.454%,n=6),RSD=1.964%,模型有效,可用于鷹嘴豆淀粉酶解參數(shù)的優(yōu)化。酶解效率較低,其原因一方面由于α-淀粉酶首先作用于較容易水解的無(wú)定形區(qū),隨后殘留的結(jié)晶區(qū)將部分地限制其水解[17],另一方面與淀粉分子中直鏈淀粉與支鏈淀粉的比例、顆粒大小、形狀、結(jié)晶結(jié)構(gòu)和結(jié)晶度等有關(guān)[18]。另外,鷹嘴豆中脂類(lèi)化合物易與直鏈淀粉形成復(fù)合物,該復(fù)合物對(duì)淀粉顆粒的糊化、膨脹和溶解具有較強(qiáng)的抑制作用,對(duì)酶水解耐受性強(qiáng)[19]。

        2.4 酶量對(duì)酶促反應(yīng)速度的影響

        酶量與反應(yīng)速度的關(guān)系見(jiàn)圖5,當(dāng)?shù)孜镞^(guò)量時(shí),酶促反應(yīng)速度與酶量成正比。時(shí)間與反應(yīng)速度的關(guān)系見(jiàn)圖6,在0~7 min時(shí)間與反應(yīng)速度線性較好,測(cè)定時(shí)間太短,反應(yīng)速度誤差較大,測(cè)定時(shí)間太長(zhǎng),則線性程度降低,綜合考慮,確定酶促反應(yīng)線性測(cè)定時(shí)間為7 min。

        圖5 酶量對(duì)酶促反應(yīng)速度的影響Fig.5 Effect of enzyme concentrationon reaction velocity

        圖6 反應(yīng)時(shí)間的確定Fig.6 Determination of reaction time

        2.5 Lineweaver-Buck法求解參數(shù)

        V-S曲線見(jiàn)圖7,以1/V為縱坐標(biāo),1/S為橫坐標(biāo),見(jiàn)圖8,Y=48.546X+5.054,r=0.992。計(jì)算得Vm=0.198 mg/(mL·min),Km=9.605 mg/mL。

        圖7 反應(yīng)速度與底物濃度的關(guān)系Fig.7 Relationship of reaction velocityand substrate concentration

        圖8 雙倒數(shù)作圖Table 8 Double-reciprocal plot

        2.6 Wilkinson統(tǒng)計(jì)法求解參數(shù)

        2.6.2 泰勒展開(kāi)式求精校解 Wilkinson法求得精確解為Vm=0.132 mg/mL·min,Km=5.736 mg/mL。

        2.6.3 Lineweaver-Buck法和Wilkinson統(tǒng)計(jì)法的比較 由表5可知,Lineweaver-Buck法和Wilkinson統(tǒng)計(jì)法計(jì)算結(jié)果有一定的差別。當(dāng)?shù)矸蹪舛容^低時(shí),Lineweaver-Buck法[20]取倒數(shù)后誤差較大,偏離直線較遠(yuǎn),且[S]適合范圍一般為在0.33 Km~2Km,因此影響Km和Vm的準(zhǔn)確性,但該法方便快速,計(jì)算量小[21]。Wilkinson統(tǒng)計(jì)法結(jié)果可靠[22]。綜合考慮各方面因素,最終采用Wilkinson統(tǒng)計(jì)法計(jì)算結(jié)果,最終得出Vm=0.132 mg/(mL·min),Km=5.736 mg/mL。

        表5 Lineweaver-Buck法和Wilkinson統(tǒng)計(jì)法數(shù)據(jù)比較

        2.7 Ea及ΔH計(jì)算

        在30~50 ℃范圍內(nèi),反應(yīng)速度隨溫度的增加而增大,以lnv對(duì)1/T(×103)作圖,見(jiàn)圖9,呈直線關(guān)系,用阿累尼烏斯(Arrhenius)方程擬合。

        式(9)

        式中:R為氣體常數(shù)8.31{J/(mol·K)};Ea為活化能(kJ/mol);B為積分常數(shù)。

        K3=A×exp(-Ea/RT)

        式(10)

        式中,A為指前因子{mg/(mL·min)}。

        以lnV為縱坐標(biāo),1/T(×103)為橫坐標(biāo),Y=4.654-2.284X,r=0.97。求得Ea=18.977 kJ/moL,A=104.962 mg/(mL·min)。

        lnKs=-ΔH/RT+lnC

        式(11)

        K3≥K2,Ks≈Km,以不同溫度下Km與1/T回歸,見(jiàn)圖10,以lnKs為縱坐標(biāo),1/T(×103)為橫坐標(biāo),Y=7.078-1.76X,r=0.99。求得ΔH=19.624kJ/mol。

        圖9 lnv與1/T的關(guān)系Fig.9 Relationship between lnv and 1/T

        圖10 lnKs與1/T的關(guān)系Table 10 Relationship between lnKs and 1/T

        3 結(jié)論

        本實(shí)驗(yàn)結(jié)合單因素實(shí)驗(yàn)和響應(yīng)面法優(yōu)化α-淀粉酶酶解鷹嘴豆淀粉的工藝參數(shù),并對(duì)酶解動(dòng)力學(xué)進(jìn)行研究。結(jié)果表明:酶解工藝為:pH 6.5,溫度55 ℃和E/S 0.12 U/mg。在此條件下,DE驗(yàn)證值為23.103%±0.454%,RSD=1.964%。利用Lineweaver-Burk及Wilkinson統(tǒng)計(jì)法進(jìn)行了數(shù)學(xué)回歸,并分析比較后選擇Wilkinson統(tǒng)計(jì)法,求解得動(dòng)力學(xué)參數(shù)Km=5.736 mg/mL,Vm=0.132 mg/(mL·min)。在30~50℃范圍內(nèi),Ea=18.977 kJ/moL,ΔH=19.624 kJ/mol,該模型推導(dǎo)過(guò)程中,各因子影響水解速度的數(shù)據(jù)擬合,方程相關(guān)性高,表明用米氏方程對(duì)α-淀粉酶酶解鷹嘴豆淀粉的動(dòng)力學(xué)過(guò)程進(jìn)行描述和模擬式有效的。

        [1]Torres-Fuentes C,Alaiz M,Vioque J. Affinity purification and characterization of chelating peptides from chickpea protein hydrolysates[J]. Food Chemistry,2011,129(2):485-490.

        [2]Chavan J K,Kadam S S,Salunkhe D K. Biochemistry and technology of chickpea(Cicer arietinum L)seeds[J]. CRC Critical Reviews in Food Science and Nutrition,1986,25:107-132.

        [3]張金波,苗昊翠,王威,等. 鷹嘴豆的應(yīng)用價(jià)值及其研究與利用[J]. 作物雜志,2011(1):10-12.

        [4]高捷,王華,寇曉虹,等. 鷹嘴豆肽抑制腫瘤作用和對(duì)免疫功能的影響[J]. 食品科學(xué),2012,33(3):215-219.

        [5]黎金,張國(guó)權(quán),羅勤貴. 蕎麥蛋白的堿性蛋白酶酶解動(dòng)力學(xué)研究[J]. 中國(guó)糧油學(xué)報(bào),2009,24(6):41-46.

        [6]鐘昔陽(yáng),楊積東,湯玉清,等. α-淀粉酶酶解小麥面粉動(dòng)力學(xué)模型研究[J]. 食品科學(xué),2012,33(7):96-100.

        [7]張國(guó)權(quán),史一一,魏益民,等. 蕎麥淀粉的真菌淀粉酶酶解動(dòng)力學(xué)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2007,23(5):42-46.

        [8]張麗華,韓永斌,孫晶,等. 響應(yīng)面法優(yōu)化蓮藕淀粉的酶解工藝參數(shù)研究[J]. 食品與發(fā)酵工業(yè),2005,31(10):81-84.

        [9]姜紹通,殷嘉憶,王華林,等. 響應(yīng)面法優(yōu)化酶法提取芋頭淀粉工藝參數(shù)[J]. 食品科學(xué),2014,35(6):24-29.

        [10]張惟杰. 糖復(fù)合物生化研究技術(shù)[M]. 杭州:浙江大學(xué)出版社(第二版),1999:10-11.

        [11]任省濤,程可可,宋安東,等. 蒸汽爆破玉米秸稈酶解動(dòng)力學(xué)[J]. 生物工程學(xué)報(bào),2011,27(4):592-597.

        [12]袁勤生,趙健. 酶與酶工程[M]. 上海:華東理工大學(xué)出版社,2005.

        [13]李忠海,劉建軍,鐘海雁,等. 響應(yīng)面法優(yōu)化甘薯淀粉酶解條件的研究[J]. 中國(guó)糧油學(xué)報(bào),2009,24(10):64-69.

        [14]齊明,房書(shū),沙同勝,等. 葡萄糖收率-糖化反應(yīng)的研究[J]. 牡丹江醫(yī)學(xué)院學(xué)報(bào),1993,(S1):29-31.

        [15]蔡麗華,馬美湖. 4種常用蛋白酶對(duì)牛骨蛋白的酶解動(dòng)力學(xué)研究[J]. 食品科學(xué),2010,31(5):150-154.

        [16]鄭寶東. 食品酶學(xué)[M]. 南京:東南大學(xué)出版社,2006,88.

        [17]二國(guó)二郎. 淀粉科學(xué)手冊(cè)[M]. 王薇青等譯. 北京:中國(guó)輕工業(yè)出版社,1990.

        [18]于軒,李兆豐,顧正彪,等. 淀粉結(jié)構(gòu)對(duì)其酶解性能影響的研究進(jìn)展[J]. 食品與發(fā)酵工業(yè),2012,38(4):137-140.

        [19]Holm J,Bjorck I,Ostrowska S,et al. Digestibility of amylase-lipid complexes in-vitro and in-vivo[J]. Starch,1983,35(9):294-297.

        [20]Komolprasert V,Ofoli R. Starch hydrolysis kinetics of Bacillus licheniformis α-amylase[J]. Journal Chem Tech Biotechnol,1991,51:209-223.

        [21]趙燕,涂勇剛,郭建軍,等. 玉米淀粉的β-淀粉酶酶解動(dòng)力學(xué)研究[J]. 食品科學(xué),2009,30(17):190-194.

        [22]鐘穗生,劉生玉,張林香. 馬鈴薯淀粉糖化動(dòng)力學(xué)的研究[J]. 中國(guó)糧油學(xué)報(bào),1998,13(1):21-25.

        Study on the processing optimization and modelling of enzymolysis for Chickpea starch based on response surface method

        GAO Yi-xia,WANG Shu-fang,WANG Dan-dan,WANG Ting-pu*

        (College of Bioengineering and Technology,Tianshui Normal University,Tianshui 741001,China)

        Using dextrose equivalent value(DE)as an index,the effect of temperature,pH and the ratio of enzyme to substrate concentration(E/S)on enzymolysis of chickpea starch with α-amylase were explored based on single experiment and response surface method. The Michaelis constant(Km),Maximum velocity(Vm)and corresponding kinetics modeling were also established by Lineweaver-Burk plotting and Wilkinson statistical method,respectively. The results indicated that the optimal parameters were pH6.5,temperature 55 ℃ and E/S value 0.12 U/mg,and the verification value was 23.103%±0.454%. Ea and ΔH were 18.977 kJ/mol and 19.624 kJ/mol in the temperature range of 30 ℃ to 50 ℃,respectively. The study might provide a theory basis for industrial production of sugar from chickpea starch.

        Chickpea;enzymolysis;kinetics;response surface method

        2014-10-21

        高義霞(1982-),女,碩士研究生,講師,從事食品化學(xué)研究,E-mail:egaoxy@yeah.net。

        *通訊作者:王廷璞(1965-),男,研究員,從事食品化學(xué)研究,E-mail:15293557513@163.com。

        國(guó)家自然科學(xué)基金(31360603)。

        TS201.1

        A

        1002-0306(2015)15-0087-05

        10.13386/j.issn1002-0306.2015.15.010

        猜你喜歡
        影響
        美食網(wǎng)紅如何影響我們吃什么
        是什么影響了滑動(dòng)摩擦力的大小
        哪些顧慮影響擔(dān)當(dāng)?
        影響大師
        沒(méi)錯(cuò),痛經(jīng)有時(shí)也會(huì)影響懷孕
        媽媽寶寶(2017年3期)2017-02-21 01:22:28
        擴(kuò)鏈劑聯(lián)用對(duì)PETG擴(kuò)鏈反應(yīng)與流變性能的影響
        基于Simulink的跟蹤干擾對(duì)跳頻通信的影響
        如何影響他人
        APRIL siRNA對(duì)SW480裸鼠移植瘤的影響
        国产成人精品午夜视频| 美女偷拍一区二区三区| 精品成人av人一区二区三区| 亚洲国产精品久久精品 | 欧洲亚洲视频免费| 国产一级av理论手机在线| 极品av一区二区三区| 精品人妻午夜一区二区三区四区| 国产精品久久久久久久成人午夜| 亚洲精品美女久久久久网站| 国产天堂av在线播放资源| 中国午夜伦理片| 夜爽8888视频在线观看| 青榴社区国产精品| 女人18毛片aa毛片免费| 男人扒开添女人下部免费视频 | 亚洲国产AⅤ精品一区二区不卡| 97成人精品在线视频| 国产果冻豆传媒麻婆精东| 久久精品视频在线看99| 亚洲欧美日韩精品香蕉| 毛片在线视频成人亚洲| 门卫又粗又大又长好爽| 日本韩无专砖码高清| 无码伊人久久大蕉中文无码| 亚洲成人一区二区av| 摸进她的内裤里疯狂揉她动图视频| 欧美精品一区二区性色a+v| 亚洲av噜噜狠狠蜜桃| 日韩亚洲无吗av一区二区| 一本一道波多野结衣av中文| 亚洲欧洲日产国产AV无码| 一区二区三区四区国产亚洲| 粉嫩小泬无遮挡久久久久久| 少妇白浆高潮无码免费区| 免费国产调教视频在线观看| 漂亮人妻出轨中文字幕| 亚洲va无码va在线va天堂| 一级片久久| 中文片内射在线视频播放| 国产欧美日韩va另类在线播放|