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

        ?

        不對稱還原大位阻二芳香基甲酮的羰基還原酶的基因克隆及性質分析

        2013-11-12 02:20:48李哲劉衛(wèi)東陳曦賈士儒吳洽慶朱敦明馬延和
        生物工程學報 2013年1期

        李哲,劉衛(wèi)東,陳曦,賈士儒,吳洽慶,朱敦明,馬延和

        1 天津科技大學生物工程學院,天津 300457

        2 中國科學院天津工業(yè)生物技術研究所 工業(yè)酶國家工程實驗室,天津 300308

        光學純二芳香基甲醇是合成許多手性藥物的重要中間體,例如:(R)-nebenodine,(R)-orpheadrine和 (S)-carbinoxamine[1-2]。使用芳香基親核試劑對芳基醛進行加成是合成二芳香基甲醇常用的方法[1],但是芳香基親核試劑有許多缺點,如價格昂貴、不適合大規(guī)模應用或產物的 ee值和收率較低等[3]。對前手性二芳香基甲酮進行不對稱還原是獲得二芳香基甲醇最直接的方法。用化學合成的方法可以達到這種目的,但只對一些具有特殊結構的底物,如苯環(huán)的鄰位有取代基或者羰基的一邊是芳香雜環(huán)時,產物的ee值才較高[4-5],而對于羰基兩側的取代基差異性較小,如對位甲基取代的二苯基甲酮,則產物的光學純度非常低(ee<26%)。

        羰基還原酶具有高度的化學、區(qū)域和立體選擇性,常用于手性醇的合成,一些酶催化的反應已經(jīng)被用于工業(yè)生產[6-8]。雖然利用分離純化的羰基還原酶催化二芳香基酮不對稱還原已有報道[3,9],但因為該類底物的位阻大、羰基兩側的取代基差異性較小,常常使得反應的活性和立體選擇性不能滿足醫(yī)藥工業(yè)的需求。因此,對二芳香基甲酮的不對稱還原仍是生物催化中很有挑戰(zhàn)性的難題之一,尋找具有較高反應活性和立體選擇性的還原較大位阻二芳香基甲酮的羰基還原酶具有很好的理論和應用價值。

        豐富多樣的微生物為發(fā)現(xiàn)新的羰基還原酶提供了寶貴的資源,然而傳統(tǒng)的菌種篩選方法耗時、低效。近年來,基因組測序技術的發(fā)展使得可以通過基因數(shù)據(jù)分析快速發(fā)現(xiàn)新的目標酶[10-11]。畢赤酵母 GS115基因組序列已被測定[12],本文通過對畢赤酵母GS115的基因組數(shù)據(jù)進行分析,發(fā)現(xiàn)一個新的羰基還原酶,能夠不對稱還原二芳香基甲酮類化合物,并對其酶學性質進行了研究,為該酶的進一步應用奠定了基礎。

        1 材料與方法

        1.1 菌株與質粒

        畢赤酵母Pichia pastoris GS115為實驗室保存,大腸桿菌Escherichia coli Rosseta2 (DE3)、DH5α、表達質粒 pET32a (+)均購自 Novagen公司。

        1.2 培養(yǎng)基、主要試劑及儀器

        P. pastoris GS115采用YPD培養(yǎng)基在30 ℃進行培養(yǎng)。大腸桿菌DH5α和Rosseta2 (DE3)在37 ℃用Luria-Bertani培養(yǎng)基進行培養(yǎng)。

        限制性內切酶、T4 DNA連接酶、DNA和蛋白Marker購于Fermentas公司;高保真PCR擴增試劑盒購自大連寶生物工程有限公司;DNA凝膠回收試劑盒購于北京賽百盛生物技術公司。二芳香基甲酮類化合物購于Sigma公司;二芳香基醇消旋體通過NaBH4還原制得;蛋白純化儀:?KTA purifier 10,所用層析柱均為GE公司,高效液相色譜儀為 Agilent 1206型;手性分離柱OD-H及AD-H柱均購于DAICEL公司。其他所用試劑均為分析純。

        1.3 目的基因的克隆表達

        P. pastoris GS115中目的基因pascr序列從EMBL數(shù)據(jù)庫獲得 (GenBank Accession No.XM_002492630.1)。根據(jù)所選擇的蛋白對應的核酸序列設計引物,正向:5¢-CGCCATATG GTTTC TAAGGTTTTATTGACA-3¢,反向:5¢-CCGCTC GAGTTTATTAGCACGCAATATCT-3¢,下劃線部分分別為NdeⅠ和XhoⅠ的酶切位點。PCR擴增條件:94 ℃預變性 5 min;94 ℃變性 30 s,55 ℃退火30 s,72 ℃延伸50 s,30個循環(huán);72 ℃延伸15 min。PCR產物通過瓊脂糖凝膠電泳分離純化,核酸回收試劑盒回收目的片段,回收產物用NdeⅠ和 XhoⅠ進行雙酶切,酶切產物再經(jīng)瓊脂糖凝膠電泳回收后用T4 DNA連接酶連接到用相同內切酶切過的質粒pET32a (+),連接產物轉化到大腸桿菌DH5α感受態(tài)中,挑取陽性克隆進行測序驗證。將測序正確的重組質粒轉化Rosseta2(DE3)感受態(tài)細胞。

        將成功轉化的單克隆接種于 10 mL含有終濃度100 μg/mL氨芐青霉素以及34 μg/mL氯霉素的LB液體培養(yǎng)基中,37 ℃、200 r/min培養(yǎng)過夜,然后轉接到含有同樣抗生素的800 mL的LB培養(yǎng)基中,培養(yǎng)到OD600為1.0左右,向培養(yǎng)基中添加終濃度為0.1 mmol/L異丙基-β-D硫代半乳糖苷 (Isopropyl-β-D-thiogalactopyranoside),37 ℃誘導4 h,于4 ℃、7 000×g離心15 min收集菌體。

        1.4 目的蛋白的分離純化

        收集菌體用生理鹽水洗滌 2次后重懸于100 mL預冷的緩沖液A中 (20 mmol/L Tris-HCl,pH 8.0,5%甘油,500 mmol/L NaCl,20 mmol/L咪唑),高壓勻漿破碎,20 000×g離心30 min,棄去細胞碎片。上清用0.22 μm濾膜過濾后上樣,上樣流速為 3 mL/min,用 B液 (20 mmol/L Tris-HCl,pH 8.0,5%甘油,500 mmol/L NaCl,500 mmol/L咪唑)進行0~250 mmol/L咪唑線性洗脫,收集不同組分。用SDS-PAGE以及活力測定檢驗目的蛋白所在組分,合并純度較高的樣,使用超濾管 (Millipore截留分子量為10 kDa)進行超濾濃縮和脫鹽。蛋白質濃度用BCA試劑盒進行測定,使用牛血清白蛋白BSA作為標準蛋白。

        1.5 酶的輔酶偏好性及分子量測定

        羰基還原酶 PasCR的輔酶偏好性是以苯乙酮為底物使用96孔板在酶標儀上進行活力測定,反應總體積為200 μL,體系中含有酶液、底物以及輔酶NADPH或NADH,通過連續(xù)監(jiān)測340 nm吸光度下輔酶的消耗 (ε=6 220 L/(mol·cm))進行測定。每分鐘消耗1 μmol輔酶所需要的酶量定義為一個酶活單位。動力學參數(shù)在最適pH下按上述方法分別以苯乙酮和 4-氯-3-羰基丁酸乙酯為底物測定。

        羰基還原酶 PasCR在溶液中的聚集狀態(tài)通過結合凝膠層析以及SDS-PAGE來進行確定。所用層析柱為Superdex 200 10/300 GL。流動相為20 mmol/L磷酸緩沖液,150 mmol/L NaCl,pH 7.2,流速為0.4 mL/min。蛋白的分子量大小通過和標準蛋白的保留體積進行計算得到。標準蛋白為:Ovalbumin (43.0 kDa),Conalbumin (75.0 kDa),Conalbumin (158.0 kDa),F(xiàn)erritin (440.0 kDa),Thyroglobulin (669.0 kDa)。

        1.6 pH及溫度對酶活的影響

        PasCR的最適反應pH通過在不同pH的緩沖液中測定酶的活力確定。緩沖液分別為100 mmol/L的醋酸/醋酸鈉緩沖液 (pH 4.0~6.0),100 mmol/L的磷酸鈉緩沖液 (pH 6.0~8.0),100 mmol/L Tris-HCl緩沖液 (pH 8.0~9.0)。PasCR的最適反應溫度通過以苯乙酮為底物在不同溫度下測定酶的活力確定。

        pH穩(wěn)定性是將酶在不同pH的緩沖液中4 ℃孵育30 h,然后測定剩余活性,不同pH緩沖液如上所述。酶熱穩(wěn)定性是將純化后的酶放在含有100 mmol/L NaCl的Tris-HCl (20 mmol/L,pH 8.0)緩沖液中不同溫度下分別孵育不同時間,然后檢測殘余酶活力。

        1.7 不同化學物質對酶穩(wěn)定性的影響

        在等量的酶液中分別加入金屬螯合劑EDTA、變性劑等不同化學物質,室溫下孵育0.5 h后,測定其殘留活性,以不處理的酶液的活性作為對照。

        1.8 酶對二芳香基甲酮類化合物的不對稱還原

        葡萄糖脫氫酶偶聯(lián)的輔酶再生體系用于催化反應。1 mL磷酸鈉緩沖液 (100 mmol/L,pH 6.5)的反應體系中包含:底物100 μL (母液100 mmol/L溶于甲醇),葡萄糖18 mg,NADP+0.5 mg,葡萄糖脫氫酶2 mg,PasCR 2 mg,30 ℃、200 r/min搖床反應。過夜反應后,用相同體積的甲基叔丁基醚萃取,取出有機相層,除去溶劑后用HPLC流動相溶解,HPLC檢測ee值及轉化率,產物的構型通過和文獻比較確定[9]。

        2 結果與分析

        2.1 羰基還原酶PasCR序列比對及分析

        基因pascr位于畢赤酵母GS115的3號染色體上,編號為 chr3_0449,其所編碼蛋白的標注為“putative dihydrokaempferol 4-reductase”。Blast結果顯示 PasCR與短鏈脫氫酶家族的蛋白有較高的同源性。PasCR與小家鼠 Mus musculus delta-5-3-β-hydroxysteroid 脫氫酶 (3β-HSD Gene ID M58567)[13]、木蘭假絲酵母 Candida macedoniensis的羰基還原酶 (MER Gene ID AB183149)[14]、赭色擲孢酵母 Sporobolomyces salmonicolor AKU4429的羰基還原酶 (SSCR Gene ID AF160799)[15]、魯氏接合酵母Zygosaccharomyces rouxii的羰基還原酶(Zygosaccharomyces gene ID AF178079)[16]、釀酒酵母 Saccharomyces cerevisiae的羰基還原酶(YOL151w Gene ID Z74893)[17]分別有 25%、38%、31%、31%、43%的同源性,序列比對的結果如圖1所示。短鏈脫氫酶家族的蛋白主要有3個功能區(qū):輔酶結合區(qū)、反應催化區(qū)及底物結合區(qū)。從比對結果可以看出PasCR的N端非常保守,主要是結合輔酶,三聯(lián)催化基團 (Tyr174、Ser135、Lys178)這在短鏈脫氫酶家族中也是非常保守的;而C端同源性較低,主要決定了酶的底物譜及立體選擇性[18]。因此,與這些已知的短鏈脫氫酶家族中的酶相比,PasCR可能具有比較獨特的底物譜及立體選擇性。

        2.2 羰基還原酶PasCR的構建表達及純化

        以畢赤酵母GS115基因組為模板,通過PCR擴增的基因片段 (1 065 bp),經(jīng)測序驗證和目的基因一致。成功構建出工程菌pET32a(+)-pascr,并在大腸桿菌 Rosseta2 (DE3)中進行了重組表達,蛋白的C端帶有6×His tag標簽。經(jīng)Ni-NTA純化,得到了高純度的目的蛋白。蛋白表達及純化SDS-PAGE圖譜見圖2。PasCR單體的分子量為 39.2 kDa,凝膠層析結果顯示酶的分子量大約為80 kDa,說明其在溶液中以二聚體形式存在。

        2.3 酶動力學參數(shù)測定

        圖1 羰基還原酶PasCR序列比對結果Fig. 1 Amino acid sequence alignment of PasCR from P. pastoris GS115 with Mus musculus delta-5-3-β-hydroxysteroid dehydrogenase (3β-HSD), Candida macedoniensis carbonyl reductase (MER), Sporobolomyces salmonicolor AKU4429 carbonyl reductase (SSCR), Zygosaccharomyces rouxii carbonyl reductase(Zygosaccharomyces )and Saccharomyces cerevisiae carbonyl reductase (YOL151w).

        圖2 PasCR蛋白純化SDS-PAGE圖譜Fig. 2 SDS-PAGE analysis of the purified PasCR.1: molecular mass standard; 2: uninduced; 3: crude cell-free extract; 4: purified PasCR.

        PasCR專一性地利用NADPH作為輔酶,在以NADH作為輔酶時對苯乙酮沒有表現(xiàn)出活力。PasCR在最適pH條件下對苯乙酮、4-氯-3-羰基丁酸乙酯和輔酶NADPH的動力學參數(shù)列于表1。從表中可以看出羰基還原酶 PasCR對脂肪族的4-氯-3-羰基丁酸乙酯的親和力和催化效率都遠高于芳香族的苯乙酮,因此 PasCR對 4-氯-3-羰基丁酸乙酯具有更高的催化活性。

        2.4 酶的最適反應pH、溫度及pH和熱穩(wěn)定性

        以苯乙酮為底物,羰基還原酶 PasCR在pH 6.5的磷酸緩沖液中表現(xiàn)出了最高活性,結果如圖 3所示。在 pH 6.5的條件下羰基還原酶PasCR在35 ℃條件下表現(xiàn)出了最高活性,而到60 ℃時沒有任何活性 (圖 4),可能在該溫度下PasCR已失活。

        表1 PasCR對苯乙酮、4-氯-3-羰基丁酸乙酯和輔酶NADPH的動力學參數(shù)Table 1 Kinetic parameters of PasCR for acetophenone, ethyl 4-chloro-3-oxobutanoate and NADPH

        圖3 酶的最適反應pHFig. 3 Effect of pH on enzyme activity in different buffers, acetic acid/sodium acetate buffer (●), sodium phosphate buffer (■)and Tris-HCl buffer (▲). The maximum activity in pH 6.5 was taken as 100%.

        圖4 酶的最適反應溫度Fig. 4 Effect of temperature on enzyme activity. The optimal temperature of the carbonyl reductase PasCR was estimated at various temperatures (25 °C ? 60 °C)in phosphate buffer (100 mmol/L, pH 6.5).

        pH穩(wěn)定性研究表明,在pH 4.5~9.0范圍內PasCR有較好的穩(wěn)定性 (圖5)。熱穩(wěn)定研究表明,在30 ℃的條件下PasCR有很好的穩(wěn)定性,而在45 ℃和 50 ℃條件下 PasCR的穩(wěn)定性較差,在45 ℃時PasCR的半衰期約80 min,而在50 ℃的條件下PasCR的半衰期只有10 min (圖6)。

        圖5 酶的pH穩(wěn)定性Fig. 5 Effect of pH on enzyme stability in different pH buffers, acetic acid/sodium acetate buffer (pH 4.5?6.0,●), sodium phosphate buffer (pH 6.0?8.0, ■)and Tris-HCl buffer (pH 8.0?9.0, ▲).

        圖6 酶的溫度穩(wěn)定性Fig. 6 Effect of temperature on enzyme stability, the samples was withdrawn at different time intervals, the activity of enzyme without pre-incubation was taken as 100%.

        2.5 化學物質對酶的穩(wěn)定性影響

        不同化學物質對酶穩(wěn)定性影響分析結果見表2,從表中可以看出金屬螯合劑EDTA對酶的活力沒有明顯影響,說明該酶不需要金屬離子的輔助。β-巰基乙醇主要是破壞二硫鍵,對酶的活性也沒有明顯影響。一些非離子型表面活性劑如曲拉通X?100當添加1%時PasCR酶活喪失了近一半,而SDS濃度為0.5%時直接使酶失活。四氫呋喃可以使酶的活性提高,但當濃度大于10%時酶的活性明顯降低 (數(shù)據(jù)未顯示)。

        表2 不同化學物質對酶活的影響Table 2 Effect of different chemicals on PasCR activity

        2.6 對二芳香基甲酮的不對稱還原

        PasCR對二芳香基甲酮類化合物表現(xiàn)出了較好的活性,對它們進行催化反應的結果如表3所示。當?shù)孜餅?-氟苯甲酮時,僅有5%的轉化率和2.7%的ee值,隨著對位取代基團位阻的變大和供電子效應的增強,氯取代和甲基取代的二苯甲酮的轉化率分別提高至40%和80%,ee值分別提高至 67%和 85%,并且產物構型均由 R變成了 S。而鄰位取代基是甲基和氯時,產物的ee值分別為78%和80%,但是轉化率都較低。當?shù)孜餅榉謩e為2-苯甲酰吡啶和3-苯甲酰吡啶時,與之前的文獻[3,9,19-21]報道不同的是,產物的轉化率達90%,但是ee值僅有22%和14%。當?shù)孜餅?,2-二苯基乙酮時,產物的轉化率為88%,ee值達到了92%。說明PasCR對二苯基取代的一系列大位阻底物表現(xiàn)出了較好的轉化率和ee值。

        表3 羰基還原酶PasCR對二芳香甲酮化合物的還原Table 3 Reduction of diary ketones by carbonyl reductase PasCR

        3 結論

        本文成功從畢赤酵母GS115中克隆出1個新的羰基還原酶基因,將其在大腸桿菌中進行表達純化,并對其酶學性質進行了研究。結果表明,該酶以NADPH作為輔酶,能夠對一些體積較大的二芳香基甲酮類物質進行不對稱還原,當苯環(huán)上對位的取代基越大時產物的ee和轉化率越高。最適反應pH為6.5,最適反應溫度為35 ℃,在溶液中以二聚體形式存在。該酶在催化一些位阻較大的二芳香基甲酮類化合物,表現(xiàn)出與已經(jīng)報道的羰基還原酶不同的特性。文獻報道來源于赭色擲孢酵母Sporobolomyces salmonicolor的及一些商業(yè)來源的羰基還原酶對對位取代的二苯基甲酮表現(xiàn)出較低的立體選擇性,而還原苯甲?;拎r,產物的ee值可達99%[3,9]。固定化的釀酒酵母細胞不對稱還原 2-苯甲?;拎ど傻拇糴e值能達到有96%[22]。與此相反,PasCR不對稱還原4-甲基二苯甲酮和4-氯二苯甲酮時,得到產物的ee較高,而對2-苯甲酰吡啶和3-苯甲酰吡啶表現(xiàn)出較低的立體選擇性。由于對位取代的二苯基甲酮的不對稱還原具有更大的挑戰(zhàn)性,來源于畢赤酵母的羰基還原酶 PasCR將具有很好的應用價值和潛力。

        [1]Schmidt F, Stemmler RT, Rudolph J, et al.Catalytic asymmetric approaches towards enantiomerically enriched diarylmethanols and diarylmethylamines. Chem Soc Rev, 2006, 35(5):454?470.

        [2]Devalia J, De Vos C, Hanotte F, et al. A randomized, double-blind, crossover comparison among cetirizine, levocetirizine, and ucb 28557 on histamine-induced cutaneous responses in healthy adult volunteers. Allergy, 2001, 56(1): 50?57.

        [3]Truppo MD, Pollard D, Devine P. Enzymecatalyzed enantioselective diaryl ketone reductions.Org Lett, 2007, 9(2): 335?338.

        [4]Welch CJ, Grau B, Moore J, et al. Use of chiral HPLC-MS for rapid evaluation of the yeast-mediated enantioselective bioreduction of a diaryl ketone. J Org Chem, 2001, 66(20):6836?6837.

        [5]Corey EJ, Helal CJ. Reduction of carbonyl compounds with chiral oxazaborolidine catalysts: A new paradigm for enantioselective catalysis and a powerful new synthetic method. Angew Chem Int Ed, 1998, 37(15): 1986?2012.

        [6]Yang Y, Zhu DM, Piegat TJ, et al. Enzymatic ketone reduction: mapping the substrate profile of a short-chain alcohol dehydrogenase (YMR226c)from Saccharomyces cerevisiae. Tetrahedron:Asymmetry, 2007, 18(15): 1799?1803.

        [7]Moore JC, Pollard DJ, Kosjek B, et al. Advances in the enzymatic reduction of ketones. Acc Chem Res,2007, 40(12): 1412?1419.

        [8]Zhu DM, Wu QQ. Engineering the enantioselectivity of biocatalysts. Chin J Biotech,2009, 25(12): 1770?1778 (in Chinese).朱敦明, 吳洽慶. 生物催化劑立體選擇性的基因工程改造. 生物工程學報, 2009, 25(12):1770?1778.

        [9]Li H, Zhu DM, Hua L, et al. Enantioselective reduction of diaryl ketones catalyzed by a carbonyl reductase from sporobolomyces salmonicolor and its mutant enzymes. Adv Synth Catal, 2009, 351(4):583?588.

        [10]Zhu DM, Mukherjee C, Biehl ER, et al. Discovery of a mandelonitrile hydrolase from Bradyrhizobium japonicum USDA110 by rational genome mining. J Biotech, 2007, 129(4): 645?650.

        [11]Yu HM, Luo H, Shi Y, et al. Application of bioinformatics in researches of industrial biocatalysis. Chin J Biotech, 2004, 20(3): 325?330(in Chinese).于慧敏, 羅暉, 史悅, 等. 生物信息學在工業(yè)生物催化研究中的應用. 生物工程學報, 2004,20(3): 325?330.

        [12]De Schutter K, Lin YC, Tiels P, et al. Genome sequence of the recombinant protein production host Pichia pastoris. Nat Biotech, 2009, 27(6):561?566.

        [13]Bain P A, Yoo M, Clarke T, et al. Multiple forms of mouse 3 beta-hydroxysteroid dehydrogenase/delta 5-delta 4 isomerase and differential expression in gonads, adrenal glands, liver, and kidneys of both sexes. Proc Nat Acad Sci USA,1991, 88(20): 8870?8874.

        [14]Kataoka M, Hoshino-Hasegawa A, Thiwthong R, et al. Gene cloning of an NADPH-dependent menadione reductase from Candida macedoniensis,and its application to chiral alcohol production.Enzyme Microb Technol, 2006, 38(7): 944?951.

        [15]Kita K, Fukura T, Nakase KI, et al. Cloning,overexpression, and mutagenesis of the Sporobolomyces salmonicolor AKU4429 gene encoding a new aldehyde reductase, which catalyzes the stereoselective reduction of ethyl 4-chloro-3-oxobutanoate to ethyl(S)-4-chloro-3-hydroxybutanoate. Appl Environ Microbiol, 1999, 65(12): 5207?5211.

        [16]Costello CA, Payson RA, Menke MA, et al.Purification, characterization, cDNA cloning and expression of a novel ketoreductase from Zygosaccharomyces rouxii. Eur J Biochem, 2000,267(17): 5493?5501.

        [17]Choi YH, Choi HJ, Kim D, et al. Asymmetric synthesis of (S)-3-chloro-1-phenyl-1-propanol using Saccharomyces cerevisiae reductase with high enantioselectivity. Appl Microbiol Biotechnol,2010, 87(1): 185?193.

        [18]Joernvall H, Persson B, Krook M, et al. Short-chain dehydrogenases/reductases (SDR). Biochemistry,1995, 34(18): 6003?6013.

        [19]Takemoto M, Achiwa K. The synthesis of optically active pyridyl alcohols from the corresponding racemates by Catharanthus roseus cell cultures.Tetrahedron: Asymmet, 1995, 6(12): 2925?2928.

        [20]Salvi NA, Chattopadhyay S. Studies on Rhizopus arrhizus mediated enantioselective reduction of arylalkanones. Tetrahedron, 2001, 57(14):2833?2839.

        [21]Takemoto M, Tanaka K. Synthesis of optically active α-phenylpyridylmethanols by Camellia sinensis cell culture. J Mol Catal B: Enzym, 2001,15(4/6): 173?176.

        [22]Takemoto M, Achiwa K. Synthesis of optically active α-phenylpyridylmethanols with Baker's yeast. Chem Pharm Bull, 1994, 42(4): 802?805.

        a级特黄的片子| 国产黄色一区二区在线看| 无套无码孕妇啪啪| 色老板精品视频在线观看| 五月婷婷激情小说| 精品人妻久久av中文字幕| 大香焦av一区二区三区| 国模丽丽啪啪一区二区| 五月激情婷婷丁香| 丰满少妇一区二区三区专区 | 91色综合久久熟女系列| 国内精品久久久久久99| 亚洲美国产亚洲av| 亚洲AⅤ乱码一区二区三区| av在线高清观看亚洲| 色综合久久精品亚洲国产| 91av在线播放| av福利资源在线观看| 国产av一区二区三区性入口| 亚洲精品乱码久久久久久蜜桃不卡| 久久频道毛片免费不卡片| 亚洲va在线∨a天堂va欧美va| 久久频这里精品99香蕉| 青青草免费观看视频免费| 亚洲精品乱码久久久久久不卡 | 亚洲人成网站免费播放| 蜜桃成人永久免费av大| 亚洲一区第二区三区四区| 毛片无码国产| 亚洲熟妇一区无码| 少妇呻吟一区二区三区| 人人妻人人澡人人爽国产一区| √天堂中文官网8在线| 亚洲啊啊啊一区二区三区| 日本免费在线不卡一区二区| 亚洲色精品aⅴ一区区三区| 国产亚洲精品国产福利在线观看| 亚洲一区二区岛国高清| 成人一区二区免费中文字幕视频| 天天躁人人躁人人躁狂躁| 另类人妖在线观看一区二区 |