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

響應(yīng)曲面法優(yōu)化大蒜中總酚提取工藝及其抗氧化活性測定

2012-10-25 03:50:24杜俊娜陳書霞程智慧常燕霞
食品科學(xué) 2012年10期
關(guān)鍵詞:工藝

杜俊娜,陳書霞*,程智慧,常燕霞,周 靜

(西北農(nóng)林科技大學(xué)園藝學(xué)院,農(nóng)業(yè)部西北地區(qū)園藝作物生物學(xué)與種質(zhì)創(chuàng)制重點實驗室,陜西 楊凌 712100)

響應(yīng)曲面法優(yōu)化大蒜中總酚提取工藝及其抗氧化活性測定

杜俊娜,陳書霞*,程智慧,常燕霞,周 靜

(西北農(nóng)林科技大學(xué)園藝學(xué)院,農(nóng)業(yè)部西北地區(qū)園藝作物生物學(xué)與種質(zhì)創(chuàng)制重點實驗室,陜西 楊凌 712100)

采用響應(yīng)曲面法對料液比、提取時間、不同體積分數(shù)酸化甲醇提取劑3個因素進行優(yōu)化,以料液比、提取時間、不同體積分數(shù)酸化甲醇提取劑為自變量,總酚提取量為響應(yīng)值,利用Box-Behnken設(shè)計原理和響應(yīng)曲面法,研究各自變量及其交互作用對總酚提取的影響,模擬求得二次多項回歸方程的預(yù)測模型,并確定料液比1:14(g/mL)、提取時間31min、80%酸化甲醇為提取劑是最佳的提取條件。在此條件下,總酚提取量為0.67mg GAE/g,模型預(yù)測值偏差為6.9%,證明所選工藝條件為最佳工藝條件。同時,用5種方法對其抗氧化活性進行測定,結(jié)果表明大蒜具有良好的抗氧化活性。

大蒜;總酚提取量;抗氧化活性;響應(yīng)曲面法

大蒜(Allium sativumL.)為百合科蔥屬植物,在我國南北均有種植,是人們?nèi)粘I钪袀鹘y(tǒng)的調(diào)味料和傳統(tǒng)的藥食同源植物。因其含有維生素、硒和含硫化合物等多種生理活性成分,因此在抗動脈粥樣硬化[1]、抗癌[2]、保健和預(yù)防[3]、治療多種慢性疾病[4]等領(lǐng)域有著廣泛的應(yīng)用,而這些應(yīng)用與它們的抗氧化活性密不可分。據(jù)研究,成熟大蒜及其成分在內(nèi)皮細胞中皆表現(xiàn)出抗氧化性損傷的功能[5],并可在一定程度上抑制癌細胞的增長[6],因此,大蒜常常被當做是功能性的食品原料被加以利用[7]。

大蒜總酚不僅是一種天然無毒的抗氧化劑,也是理想的天然藥物,具有清除自由基、抗菌、抗病毒、抗氧化、抗腫瘤等作用。同時在醫(yī)藥、化妝品、農(nóng)用化學(xué)品等領(lǐng)域也具有重要的應(yīng)用。在目前大蒜多酚提取的工藝方法中,超聲波提取法具有提取效率高、不需高溫、能耗低、提取時間短等特點。但對于大蒜料液比、提取時間、不同體積分數(shù)提取劑等工藝參數(shù)尚有較少的研究對其進行優(yōu)化。

響應(yīng)曲面法(response surface methodology,RSM)是一種優(yōu)化反應(yīng)條件和加工工藝參數(shù)的有效方法,可用來進行建模、檢驗?zāi)P秃线m性、因子效應(yīng)的評估、考察以及尋求因子的最佳操作條件等[8]。近年來響應(yīng)曲面法已成功地應(yīng)用于天然化合物提取優(yōu)化實驗[9],與正交試驗設(shè)計相比,具有試驗周期短,回歸方程的精度高,能同時研究幾種因素間的交互作用等優(yōu)點[10]。

本研究在單因素試驗基礎(chǔ)上,利用響應(yīng)曲面法對大蒜中總酚的提取工藝進行優(yōu)化,分析料液比、提取時間、不同體積分數(shù)提取劑對總酚提取得率的影響,建立總酚提取的數(shù)學(xué)模型。并采用5種方法分析大蒜的抗氧化活性,為大蒜醫(yī)學(xué)功能的評價和相關(guān)產(chǎn)品的開發(fā)提供科學(xué)依據(jù)。

1 材料與方法

1.1 材料與試劑

大蒜(Allium sativumL.)鱗莖于2010年8月取自西北農(nóng)林科技大學(xué)園藝場大蒜資源圃。采新鮮鱗莖-20℃保存?zhèn)溆谩?/p>

福林肖卡試劑、1,1-二苯基-2-三硝基苯肼(1,1-diphenyl-2-picrylhydrazyl radical 2,2-diphenyl-1-(2,4,6-trinitrophenyl)hydrazyl,DPPH)、新亞銅、啡咯嗪、rolox標樣 美國Sigma公司。其他試劑均為分析純。

1.2 儀器與設(shè)備

SHIMADZU UV-1700型紫外分光光度計 日本島津公司;ZMQS 5001型超純水制備儀 法國Milipore公司;SORVAIL RC-SC-PLUS型高速冷凍臺式離心機 美國Kendro公司;KQ-300DE型數(shù)控超聲波機 昆山儀器有限公司。

1.3 方法

1.3.1 提取工藝

準確稱取一定量的大蒜鱗莖于研缽中,加入液氮研磨至粉末,轉(zhuǎn)移至50mL離心管中,再加入相應(yīng)量的不同提取液,在20℃條件下進行超聲波提取,10000r/min離心10min,轉(zhuǎn)移上清至絲口瓶中,-20℃貯存?zhèn)溆谩?/p>

1.3.2 提取條件優(yōu)化方案

根據(jù)單因素試驗,選取對試驗結(jié)果影響較大的因素,采用Box-Benhnken模型進行試驗設(shè)計,用Design Expert 7.0軟件進行響應(yīng)面優(yōu)化分析。

1.3.3 總酚提取量的測定

采用Folin-Ciocalteu(FC)法[11]。移取待測溶液300μL提取液,加入700μL蒸餾水,混合均勻后,加入200μL福林肖卡(FC)試劑,混合萃取1min后,加入2mL Na2CO3(10%)溶液,室溫避光靜置2h后,于波長765nm處測吸光度。以沒食子酸為標準品,繪制標準曲線,得回歸方程y=17.179x+0.0147,r2=0.9991,根據(jù)標準曲線方程求出提取液中總酚質(zhì)量濃度/(mg GAE/g,GAE(gallic acid equivalent)表示沒食子酸等效物)。

1.3.4 抗氧化活性的測定

1.3.4.1 DPPH清除力

通過檢測試驗樣品對有機自由基DPPH的清除能力可以評價其抗氧化性的強弱。試驗測定參照Mimica等[12]的方法。DPPH清除率按下式計算。

式中:Ai為試驗樣品吸光度;Ac為空白對照吸光度。

1.3.4.2 羥自由基清除力

對·OH的清除能力測定參照Sakanaka等[13]的方法。清除率按下式進行計算。

式中:Ai為試驗樣品吸光度;Ac為空白對照吸光度。

1.3.4.3 鐵氰化鉀還原力

樣品的還原能力與抗氧化能力呈正相關(guān),參照Oyaizu[14]的方法進行測定,還原能力的強弱以波長700nm處的吸光度表示。

1.3.4.4 銅離子還原力

取1mL提取液,分別加入1mL CuSO4·5H2O(5mmoL/L),1mL新亞銅(3.75mmoL/L),于37℃反應(yīng)30min,450nm波長處測吸光度。Trolox作標品,繪制標準曲線,得回歸方程y=0.0016x-0.0175,r2=0.998,根據(jù)標準曲線方程求出提取液中銅離子還原能力相當于trolox的含量的質(zhì)量濃度。結(jié)果用mg(trolox等同量)/g表示。

1.3.4.5 金屬螯合力

對亞鐵的螯合能力測定參照Ak等[15]的方法。螯合率按下式進行計算。

式中:Ai為試驗樣品吸光度;Ac為空白對照吸光度。

2 結(jié)果與分析

2.1 單因素試驗結(jié)果

2.1.1 料液比對大蒜中總酚提取量的影響

圖1 料液比對總酚提取量的影響Fig.1 Effect of material / liquid ratio on yield of total phenolic compounds

以80%酸化甲醇為提取劑,提取時間為30min的條件下,以料液比1:5、1:10、1:15、1:20、1:25、1:30(g/mL)進行試驗,不同料液比對大蒜總酚提取量的影響見圖1。由圖1可知,隨著液體比例的逐漸增加,總酚提取量呈增大趨勢。但當料液比達到一定量時,提取劑對有效成分的提取達到飽和。再增加提取劑的量,雜質(zhì)成分會競爭溶出而不利于有效成分的提取。而且會造成的浪費,成本增加。因此,在保證提取效果的同時,可減少提取劑的用量。綜合考慮,料液比以控制在1:15(g/mL)為宜。

2.1.2 提取時間對大蒜中總酚提取量的影響

圖2 提取時間對總酚提取量的影響Fig.2 Effect of extraction time on yield of total phenolic compounds

以80%酸化甲醇為提取劑,料液比1:10(g/mL)的條件下,以提取時間20、30、40、50、60min進行試驗,不同提取時間對總酚提取量的影響見圖2。

由圖2可知,隨著時間延長,總酚提取量逐漸增加,在30min前增加較明顯,說明時間越短提取越不充分,在30min后總酚提取量開始下降,之后回升,趨勢平緩,說明大蒜中抗氧化物質(zhì)已基本提取完,為縮短工時、節(jié)省能源,選取最佳提取時間為30min。

2.1.3 提取劑體積分數(shù)對大蒜中總酚提取量的影響

在料液比為1:10(g/mL),提取時間30min條件下,分別用50%、60%、70%、80%、90%、100%的酸化甲醇進行提取。不同體積分數(shù)提取劑對總酚提取量的影響見圖3。由圖3可知,80%酸化甲醇提取率最高,所以選擇80%酸化甲醇為最佳提取劑。

圖3 提取劑體積分數(shù)對總酚提取量的影響Fig.3 Effect of concentration of extracting agents on yield of total phenolic compounds

2.2 響應(yīng)面法提取條件的優(yōu)化

2.2.1 響應(yīng)面試驗設(shè)計及結(jié)果

在單因素試驗的基礎(chǔ)上,采用Box-Behnken中心組合試驗設(shè)計原理,以料液比、提取時間和不同體積分數(shù)提取劑3個外界因素為自變量,以總酚提取量為響應(yīng)值即因變量,因素編碼及水平見表1,試驗設(shè)計及結(jié)果見表2。

表1 響應(yīng)面試驗因素水平編碼表Table 1 Variables and levels in three-variable-three-level Box-Behnken experimental design

表2 Box-Behnken試驗設(shè)計及結(jié)果Table 2 Trials and results of the three-variable-three-level Box-Behnken experimental design

利用Design Expert 7.0軟件,對表2的試驗數(shù)據(jù)進行處理,得到回歸方程方差分析表,見表3。同時利用該軟件進行非線性回歸分析,得到預(yù)測模型如下:

表3 回歸方程方差分析表Table 3 Analysis of variances for the developed regression equation

由表3可以看出,F(xiàn)回歸=5.33,P=0.0191<0.05,表明模型效應(yīng)顯著,不同處理間的差異顯著;各因素中A2、B2、C2效應(yīng)均顯著;A、B、C、A B、A C、BC效應(yīng)不顯著,由此可見,各具體試驗因素對響應(yīng)值的影響不是簡單的線性關(guān)系。F失擬=0.64,說明失擬項與凈誤差沒有顯著的關(guān)聯(lián),失擬項P=0.6262>0.05,不顯著,說明回歸方程擬合程度良好。R2=0.8727,說明該回歸方程回歸效果較好;模型的調(diào)整決定系數(shù)R2Adj=0.7091,說明該模型能解釋70.91%響應(yīng)值的變化,因而該模型擬合程度良好,實驗誤差小,可信度高,該模型是合適的,可以用此模型對總酚提取量提取工藝進行分析和預(yù)測。

2.2.2 各因素及交互作用對總酚提取量的影響

圖4 各兩因素交互作用對總酚提取量影響的等高線及響應(yīng)面圖Fig.4 Response surface plots of the interactive effects of extraction time and concentration of extracting agents on yield of total phenolic compounds

由圖4a可知,以80%酸化甲醇為提取劑,料液比和提取時間對大蒜中總酚提取量的交互影響。在提取時間不變的條件下,隨著提取劑用量的逐漸增加,總酚提取量出現(xiàn)先上升后下降的趨勢,變化較不顯著。這可能是當料液比變大時,大蒜中的總酚類物質(zhì)與提取劑在單位時間內(nèi)存在較大的質(zhì)量濃度梯度,擴散系數(shù)大,擴散速率大。隨著時間的延長,有效組分在提取劑中的溶解基本趨于飽和,再加大料液比也不會促進提取。在料液比不變的條件下,隨著時間的逐漸延長,總酚提取量也出現(xiàn)了先上升后下降的趨勢,變化較為顯著,這說明總酚提取量隨著時間的延長而增加,當達到飽和時,再延長提取時間,反而會抑制總酚類物質(zhì)的提取。

由圖4b可知,在30min提取時間下,料液比和提取劑體積分數(shù)對大蒜中總酚提取量的交互影響。在料液比不變的條件下,總酚提取量隨提取劑體積分數(shù)的升高先增加后降低,提取劑體積分數(shù)大致在80%時,總酚提取量最高,70%和90%時,總酚提取量較低。在提取劑濃度不變的條件下,隨著提取劑用量的逐漸增加,總酚提取量出現(xiàn)先上升后下降的趨勢。這可能是因為料液比的增大會使大蒜中的總酚類物質(zhì)得到更充分地釋放,總酚提取量上升,當大蒜中的總酚提取接近完全時,再繼續(xù)增大提取劑用量,總酚提取量則出現(xiàn)下降。

由圖4c可知,在1:15(g/mL)料液比條件下,提取時間和提取劑體積分數(shù)對大蒜中總酚提取量的交互影響。在提取時間不變的條件下,同料液比和不同體積分數(shù)提取劑的交互作用,以80%酸化甲醇為提取劑的總酚提取量最高。在提取劑體積分數(shù)不變的條件下,隨著提取時間的增加,總酚提取量出現(xiàn)先上升后下降的趨勢。其原因可能是在剛開始提取的時候,大蒜中的總酚類物質(zhì)與提取劑在單位時間內(nèi)存在較大的質(zhì)量濃度差,擴散速度快,總酚提取量較高,隨著時間的延伸,部分組織的有效組分的結(jié)構(gòu)破壞,導(dǎo)致總酚提取量下降。

綜合圖4,在模型濃度范圍內(nèi)選擇出發(fā)點,按照模型(1)使用快速上升法進行提取工藝的優(yōu)化,提取的最佳條件為料液比1:13.76(g/mL)、提取時間30.91min、提取劑體積分數(shù)79.92%,在此條件下總酚提取量的理論值0.72mg GAE/g,可信度為0.7091。

2.2.3 最佳工藝參數(shù)的驗證

在最佳條工藝件下進行驗證實驗,修正工藝參數(shù)為料液比1:14(g/mL)、提取時間31min、用80%酸化甲醇為提取劑,實際測得的總酚提取量0.67mg GAE/g,模型預(yù)測值偏差6.9%,與理論值較為一致。因此,認為利用響應(yīng)面分析法得到的優(yōu)化提取工藝參數(shù)準確可靠,得到的提取條件具有一定的實際應(yīng)用價值。

2.3 最佳工藝提取液的抗氧化活性的測定

對采用最佳工藝所得提取液進行抗氧化活性測定。如圖5所示,DPPH清除率76.2%、羥自由基清除率52.96%、鐵氰化鉀還原力(實際吸光度)為0.289、銅離子還原力相當于1.888mg(trolox等同量)/g、而金屬螯合率為21.43%,說明大蒜鱗莖提取液具有很好的抗氧化活性。

圖 5 總酚提取量及抗氧化活性指標Fig.5 Total phenolic content extracted under optimized extraction conditions and antioxidant activity

3 結(jié) 論

通過單因素試驗和 Box-Behnken中心組合設(shè)計原理以及響應(yīng)面分析法對大蒜中總酚類物質(zhì)的提取工藝進行優(yōu)化,得到料液比、提取時間、不同體積分數(shù)提取劑3個因素與總酚提取量的回歸模型,經(jīng)檢驗證明該模型合理可靠,能較好地預(yù)測大蒜中總酚提取量。由該模型確定的最優(yōu)工藝條件為料液比1:14(g/mL)、提取時間31min、提取劑為80%酸化甲醇,在此條件下,得到總酚提取量為0.67mg GAE/g,與理論值0.72mg GAE/g較為一致。同時,抗氧化活性試驗結(jié)果表明大蒜具有很好的抗氧化活性。

[1] NENCINI C, CAVALLO F, CAPASSO A, et al. Evaluation of antioxidative properties ofAlliumspecies growing wild in Italy[J].Phytother Res, 2007, 21(9): 874-878.

[2] LEELARUNGRAYUB N, RATTANAPANONE V, CHANARAT N, et al. Quantitative evaluation of the antioxidant properties of garlic and shallot preparations[J]. Nutrition, 2006, 22(3): 266-274.

[3] TAPIERO H, TEW K D, BA G N, et al. Polyphenol:do they play a role in the prevention of human pathologies[J]. Biomed Pharmacother, 2002,56(4): 200-207.

[4] DIPLOCK A T, CHARLERX J L, CROZIER W G, et al. Functional food science and defence against reactive oxidative species[J]. British Journal of Nutrition, 1998, 80(Suppl 1): 77-112.

[5] IDE N, LAU B H. Garlic compounds protect vascular endothelial cells from oxidized low density lipoprotein-induced injury[J]. J Pharm Pharmacol, 1997, 49(9): 908-911.

[6] HU X, CAO BN , HU G, et al. Attenuation of cell migration and induction of cell death by aged garlic extract in rat sarcoma cells[J]. Int J Mol Med, 2002, 9(6): 641-643.

[7] JAE-HEE P, YOO K P, EUNJU P. Antioxidative and antigenotoxic effects of garlic(Allium sativumL.) prepared by different processing methods[J]. Plant Foods Hum Nutr, 2009, 64(4): 244-249.

[8] ANNADURAI G, SHEEJA R Y. Use of Box-Behken design of experiments for the adsorption of verofix-red using biopolymer[J]. Bioprocess Eng, 1998, 18(6): 463-466.

[9] 張囡, 康廷國, 尹海波. 星點設(shè)計-響應(yīng)面法優(yōu)化穿龍薯蕷多糖提取工藝[J]. 中藥材, 2011, 34(1): 123-126.

[10] 劉洋, 趙謀明, 楊寧. 響應(yīng)面分析法優(yōu)化仙人掌多糖提取工藝的研究[J]. 食品與機械, 2006, 22(6): 42-45.

[11] BILJANA B, NEDA M D, ISIDORA S, et al. Phenolics as antioxidants in garlic (Allium sativumL., Alliaceae)[J]. Food Chemistry, 2008, 111(4): 925-929.

[12] MIMICA D, BOZIN B, SOKOVICM, et al. Antimicrobial and antioxidant activities ofMelissa oficinalisL. (Lamiaceae) essential oil[J].Journal of Agricultural and Food Chemistry, 2004, 52(9): 2485-2489.

[13] SAKANAKA S, ISHIHARA Y. Comparison of antioxidant properties of persimmon vinegar and some other commercial vinegars in radicalscavenging assays and on lipid oxidation in tuna homogenates[J]. Food Chemistry, 2008, 107(2): 739-744.

[14] OYAIZU M. Studies on products of browning reaction: antioxidative activity of products of browning reaction[J]. Japanese Journal of Nutrition,1986, 44(6): 307-315.

[15] AK T, GLIN I. Antioxident and radical scavenging properties of curumin[J]. Chemico Biological Interactions, 2008, 174(1): 27-37.

Response Surface Methodology for Optimization of Extraction Conditions of Total Phenolic Compounds from Garlic and Determination of Antioxidant Activity

DU Jun-na,CHEN Shu-xia*,CHENG Zhi-hui,CHANG Yan-xia,ZHOU Jing
(Key Laboratory of Horticultural Plant Germplasm Resources Utilization in Northwest China, Ministry of Agriculture,College of Horticultural Science, Northwest A&F University, Yangling 712100, China)

Response Surface Methodology (RSM) was used to optimize extraction parameters (including material/liquid ratio,extraction time and concentration of extracting agents) so as to achieve the maximum extraction yield of total phenolic compounds.Studies were carried out on the respective independent variables (material/liquid ratio, extraction time and concentration of extracting agents) and the interactive effects of these three factors on the total phenolic compounds according to Box-Benhnken design. A mathematical quadratic polynomial regression equation reflecting the relationship between the total phenolic compounds and the above extraction parameters was established. The optimal extraction parameters were obtained as: material/liquid ratio 1:14 (g/mL) , extraction time 31 min, and 80% acidification methanol as extracting agent. Under the optimal condition,the total phenolic content was 0.67 mg GAE/g, and the model prediction deviation was 6.9%. Meanwhile, the antioxidant activity was determined with five different methods, and the results showed that garlic had good antioxidant activity.

garlic;total phenolic compounds;antioxidant activity;response surface methodology

S663.4

A

1002-6630(2012)10-0072-06

2011-05-21

西北農(nóng)林科技大學(xué)唐仲英育種基金項目(A212020911);西北農(nóng)林科技大學(xué)青年學(xué)術(shù)骨干支持計劃項目(01140303);國家公益性行業(yè)(農(nóng)業(yè))科研專項(200903018-7)

杜俊娜(1985—),女,碩士研究生,研究方向為蔬菜生理與分子育種。E-mail:dujunna2005@126.com

*通信作者:陳書霞(1971—),女,副教授,博士,研究方向為蔬菜生理與分子育種。E-mail:shuxiachen@nwsuaf.edu.cn

猜你喜歡
工藝
鋯-鈦焊接工藝在壓力容器制造中的應(yīng)用研究
金屬鈦的制備工藝
轉(zhuǎn)爐高效復(fù)合吹煉工藝的開發(fā)與應(yīng)用
山東冶金(2019年6期)2020-01-06 07:45:54
工藝的概述及鑒定要點
收藏界(2019年2期)2019-10-12 08:26:06
5-氯-1-茚酮合成工藝改進
螺甲螨酯的合成工藝研究
壓力缸的擺輾擠壓工藝及模具設(shè)計
模具制造(2019年3期)2019-06-06 02:11:00
石油化工工藝的探討
一段鋅氧壓浸出與焙燒浸出工藝的比較
FINEX工藝與高爐工藝的比較
新疆鋼鐵(2015年3期)2015-11-08 01:59:52
主站蜘蛛池模板: 久久国产精品电影| 久久精品国产999大香线焦| 激情视频综合网| 国产乱子伦精品视频| 找国产毛片看| 91精品国产自产在线老师啪l| 国产精品综合色区在线观看| 毛片免费视频| 国产美女无遮挡免费视频| 人妻丰满熟妇AV无码区| 亚洲天堂在线视频| 精品人妻AV区| 国产精品不卡永久免费| 国产自在线播放| 一级香蕉人体视频| 国产91无码福利在线| 国产极品粉嫩小泬免费看| 小13箩利洗澡无码视频免费网站| 在线观看视频一区二区| 伊人精品视频免费在线| 日韩一区二区在线电影| 日韩大乳视频中文字幕 | 伊人成人在线视频| 中国国产A一级毛片| 亚洲一区无码在线| 免费人欧美成又黄又爽的视频| 国模沟沟一区二区三区| 日本午夜精品一本在线观看| 波多野结衣一区二区三区88| 91麻豆国产精品91久久久| 亚洲欧美一级一级a| www.日韩三级| 99视频在线精品免费观看6| 久久午夜影院| 久久精品免费国产大片| 欧美自慰一级看片免费| 最新日本中文字幕| 亚洲香蕉在线| 一区二区三区在线不卡免费| av一区二区人妻无码| 国产精品深爱在线| 欧美日韩高清在线| 亚洲人成在线免费观看| 日韩毛片免费| 国产精品免费入口视频| 精品久久高清| 99热国产这里只有精品无卡顿"| 一级香蕉人体视频| 久久人搡人人玩人妻精品| 国产精品第| 成人福利免费在线观看| 在线播放真实国产乱子伦| 亚洲日韩精品无码专区97| 国产黄色免费看| 国产一区二区影院| 亚洲国产精品不卡在线| 亚洲人成网7777777国产| 91综合色区亚洲熟妇p| 一级片免费网站| 91在线国内在线播放老师| 精品欧美一区二区三区久久久| 国产亚洲欧美在线视频| 亚洲国产无码有码| 久久人与动人物A级毛片| 午夜日b视频| 成人免费黄色小视频| 极品尤物av美乳在线观看| 色135综合网| 亚洲成人一区二区三区| 日韩少妇激情一区二区| 免费看a毛片| 日本www在线视频| 精品亚洲国产成人AV| 国产91线观看| 97se亚洲综合在线天天| 91精品专区| AV无码无在线观看免费| 欧美日韩中文国产va另类| 91福利免费视频| 高潮毛片免费观看| 亚洲三级色| 日韩成人在线网站|