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

    定點(diǎn)突變提高苯丙氨酸羥化酶的熱穩(wěn)定性

    2016-11-01 02:12:56葉雙雙周麗周哲敏
    生物工程學(xué)報(bào) 2016年9期
    關(guān)鍵詞:羥化酶苯丙氨酸熱穩(wěn)定性

    葉雙雙,周麗,周哲敏

    ?

    定點(diǎn)突變提高苯丙氨酸羥化酶的熱穩(wěn)定性

    葉雙雙,周麗,周哲敏

    江南大學(xué)生物工程學(xué)院,江蘇無(wú)錫 214122

    嗜熱菌中,蛋白質(zhì)存在Ala替換Gly以及Arg替換Lys的趨勢(shì)。為了提高紫色色桿菌來(lái)源的苯丙氨酸羥化酶的熱穩(wěn)定性,將該酶中所有Gly突變成Ala,Lys突變成Arg,篩選獲得熱穩(wěn)定性提高的突變體,并進(jìn)行組合突變,對(duì)突變酶的酶學(xué)性質(zhì)進(jìn)行研究。結(jié)果表明,突變酶K94R和G221A在50 ℃的半衰期分別為 26.2 min、16.8 min,比原始酶 (9.0 min) 分別提高了1.9倍、0.9倍,同時(shí)組合突變酶K94R/G221A在50 ℃處理1 h后仍保留65.6%的酶活,比原始酶 (8.6%) 高出6.6倍。圓二色譜結(jié)果顯示原始酶和突變酶K94R、G221A及K94R/G221A的m值分別為51.5 ℃、53.8 ℃、53.1 ℃和54.8 ℃。蛋白三維結(jié)構(gòu)模擬推測(cè)突變體熱穩(wěn)定性提高機(jī)理為:突變體K94R中Arg94與Ile95之間形成額外氫鍵,穩(wěn)定其所在的柔性區(qū)域;突變體G221A中Ala221與Leu281產(chǎn)生疏水作用,穩(wěn)定酶分子C-端柔性區(qū)。該研究結(jié)果為蛋白質(zhì)熱穩(wěn)定性改造提供了參考,也為苯丙氨酸羥化酶在功能性食品領(lǐng)域的應(yīng)用奠定了基礎(chǔ)。

    苯丙氨酸羥化酶,熱穩(wěn)定性,定點(diǎn)突變,丙氨酸突變,精氨酸突變

    苯丙氨酸在人體中主要代謝途徑是通過(guò)苯丙氨酸羥化酶 (Phenylalanine hydroxylase,PAH,EC 1.14.16.1) 的催化生成酪氨酸。若人體苯丙氨酸羥化酶存在缺陷,將導(dǎo)致血液中苯丙氨酸含量升高,最終形成苯丙酮尿癥[1]( Phenylketonuria,PKU)。目前,國(guó)際上苯丙酮尿癥的臨床治療方法都是攝入低苯丙氨酸食品或人工合成營(yíng)養(yǎng)食品[2]。這種治療方式雖然能在一定程度上控制苯丙酮尿癥,但是也會(huì)導(dǎo)致患者其他營(yíng)養(yǎng)素缺乏[3-6]。如果能將苯丙氨酸羥化酶添加到食物中,直接催化苯丙氨酸生成酪氨酸,可有效控制血液中苯丙氨酸含量,同時(shí)不影響其他營(yíng)養(yǎng)物質(zhì)的攝入。

    苯丙氨酸羥化酶是一種非血紅素鐵單加氧酶,廣泛存在于動(dòng)物、植物和微生物中,是苯丙氨酸代謝過(guò)程中的關(guān)鍵酶和限速酶??蒲泄ぷ髡邔?duì)多種來(lái)源的苯丙氨酸羥化酶進(jìn)行研究,發(fā)現(xiàn)與其他來(lái)源的苯丙氨酸羥化酶相比,紫色色桿菌來(lái)源的苯丙氨酸羥化酶 (CvPAH) 具有結(jié)構(gòu)簡(jiǎn)單[7]、與人類來(lái)源苯丙氨酸羥化酶催化區(qū)域幾乎完全一致[8]、能夠在人體溫度和pH條件下催化苯丙氨酸羥化反應(yīng)[7,9]、催化活性及穩(wěn)定性高[7]等特點(diǎn)。至2013年,Yew等[10]已通過(guò)紅細(xì)胞封裝的手段將CvPAH導(dǎo)入正常小鼠中,并降低小鼠血液中苯丙氨酸的水平,為其在功能食品領(lǐng)域的運(yùn)用提供了參考。

    本研究室在前期研究中,從紫色色桿菌中提取基因,實(shí)現(xiàn)了在大腸桿菌中高效表達(dá),并對(duì)其酶學(xué)性質(zhì)進(jìn)行了系統(tǒng)地研究[9]。結(jié)果表明該酶在37 ℃、pH 7.5時(shí)活性高于其他來(lái)源的苯丙氨酸羥化酶。然而該酶在高溫條件下迅速失活,增加了其醫(yī)藥應(yīng)用的成本。因此,提高CvPAH的熱穩(wěn)定性對(duì)其應(yīng)用至關(guān)重要。

    定向進(jìn)化能夠通過(guò)引入隨機(jī)突變提高酶的熱穩(wěn)定性,但是該方法需要構(gòu)建足夠大的突變體庫(kù),因此往往消耗大量的人力、物力及財(cái) 力[11–12];計(jì)算機(jī)模擬設(shè)計(jì)也能提高蛋白熱穩(wěn)定性,然而影響蛋白質(zhì)穩(wěn)定性的因素復(fù)雜,大大降低計(jì)算機(jī)預(yù)測(cè)的準(zhǔn)確性[13];基于B-factor的進(jìn)化策略能夠降低酶局部的柔性達(dá)到降低酶的熱穩(wěn)定性的目的,但是該方法需要建立在獲得晶體結(jié)構(gòu)的基礎(chǔ)上[12],限制了該方法的應(yīng)用。影響蛋白質(zhì)熱穩(wěn)定性的因素眾多[14-20],蛋白質(zhì)的氨基酸組成與其熱穩(wěn)定性關(guān)系密切。如果能夠以蛋白質(zhì)的一級(jí)結(jié)構(gòu)為基礎(chǔ),建立提高酶熱穩(wěn)定性的方法,不僅避免了高級(jí)結(jié)構(gòu)的限制,還能快速準(zhǔn)確地獲得突變體。

    1979年,Argos等[21]對(duì)嗜熱蛋白和嗜中溫蛋白氨基酸組成進(jìn)行了統(tǒng)計(jì)分析,發(fā)現(xiàn)相對(duì)于嗜中溫蛋白,嗜熱蛋白中氨基酸序列中有如下替換趨勢(shì):即甘氨酸 (Glycine,Gly) 替換為丙氨酸 (Alanine,Ala),賴氨酸 (Lysine,Lys) 替換為精氨酸 (Arginine,Arg)。Zhang 等[22]在T4溶菌酶中介入了單個(gè)或多個(gè)Ala突變,有效提高該酶的熱穩(wěn)定性。此外,Mrabet等[23]利用將Lys突變成Arg的策略有效提高木聚糖異構(gòu)酶的熱穩(wěn)定性。然而,尚未見將Ala和Arg突變策略組合的研究,以及應(yīng)用上述策略提高PAH熱穩(wěn)定性的研究。

    本實(shí)驗(yàn)選取CvPAH上所有的Gly及Lys,分別突變成Ala和Arg,構(gòu)建突變體,篩選熱穩(wěn)定性提高的突變酶,進(jìn)一步進(jìn)行組合突變,并研究突變體的酶學(xué)性質(zhì)。

    1 材料與方法

    1.1 材料

    1.1.1 菌種與培養(yǎng)基

    菌種:大腸桿菌JM109、BL21 (DE3) 購(gòu)自Novangen公司,BL21 (DE3)/pET24a-由本實(shí)驗(yàn)室構(gòu)建并保存[9]。

    LB培養(yǎng)基 (g/L):胰蛋白胨10,酵母膏5,NaCl 10,瓊脂粉20 (配制固體培養(yǎng)基)。

    2YT培養(yǎng)基 (g/L):胰蛋白胨16,酵母膏10,NaCl 10,瓊脂粉20 (配制固體培養(yǎng)基)。

    1.1.2 工具酶和生化試劑

    工具酶、DNA Marker和Protein Marker購(gòu)自寶生物工程 (大連) 有限公司;質(zhì)粒提取試劑盒購(gòu)自天根生化科技有限公司;膠回收試劑盒和PCR產(chǎn)物純化試劑盒購(gòu)自東洋紡 (上海) 生物科技有限公司;引物由上海生物工程公司合成;酵母提取物、蛋白胨購(gòu)自 Oxford 公司;IPTG、卡那霉素購(gòu)自上海生物工程公司;其余試劑均為國(guó)產(chǎn)分析純。

    1.2 方法

    1.2.1 突變位點(diǎn)選擇及引物設(shè)計(jì)

    根據(jù)嗜熱蛋白對(duì)Ala的偏好性明顯高于Gly、對(duì)Arg的偏好性高于Lys的特點(diǎn),選取CvPAH上所有Gly和Lys分別突變成Ala和Arg。根據(jù)GenBank提供的CvPAH基因序列 (Gene ID:AF146711),以重組質(zhì)粒pET24a-為模板,用軟件Primer 5.0設(shè)計(jì)引物 (表1)。

    表1 本文所用引物及序列

    續(xù)表1

    G290A upTGCCGGCGACCGCCAAGCCTGGGCGGATACC31 G290A downCTTCGGTATCCGCCCAGGCTTGGCGGTCGC C31 K16R upCCGACATCACCACCCGCCGCAATGTCGGACTG32 K16R downGGCTCAGTCCGACATTGCGGCGGGTGGTGATG32 K53R upGTACCAGCGCCAATGCCGCCTGCTGCCCGGC31 K53R downCGCGGCCGGGCAGCAGGCGGCATTGGCGCTG31 K82R upGGTGCCGGACTTCAATCGCCTCAACGAGAAG31 K82R downTCAGCTTCTCGTTGAGGCGATTGAAGTCCGG31 K86R upCAATAAGCTCAACGAGCGCCTGATGGCCGCC31 K86R downCGGTGGCGGCCATCAGGCGCTCGTTGAGCTT31 K94R upGGCCGCCACCGGCTGGAGGATCGTCGCGGTG31 K94R downCCGGCACCGCGACGATCCTCCAGCCGGTGGC31 K161R upCCTGGAGGCCTACGGCCGCGGCGGGGTGAAG31 K161R downTCGCCTTCACCCCGCCGCGGCCGTAGGCCTC31 K165R upCGGCAAGGGCGGGGTGCGCGCGAAGGCGCTG31 K165R downCGCCCAGCGCCTTCGCGCGCACCCCGCCCTTG32 K167R upGGGCGGGGTGAAGGCGCGCGCGCTGGGCGCG31 K205R upCGGCATCTTGTCCAGCCGCTCGGAATCCATC31 K205R downAGTAGATGGATTCCGAGCGGCTGGACAAGATG32 K239R upGATCGACACCTTCCAGCGCACCTACTTCGTC31 K239R downCGATGACGAAGTAGGTGCGCTGGAAGGTGTC31 K248R upCGTCATCGACAGCTTCCGGCAGCTGTTCGAC31 K248R downTGGCGTCGAACAGCTGCCGGAAGCTGTCGATG32

    The mutation sites are underlined.

    1.2.2 突變體的構(gòu)建

    PCR反應(yīng)體系:10 μL 5×PrimeSATAR Buffer (5 mmol/L Mg2+plus),4 μL濃度的dNTPs (2.5 mmol/L each),1 μL上游引物 (10 mmol/L),1 μL下游引物 (10 mmol/L),1 μL模板pET24a-,0.5 μL PrimeSTAR HS DNA Polymerase,超純水補(bǔ)至50 μL。

    PCR擴(kuò)增條件:95 ℃預(yù)變性5 min;98 ℃變性1 min,58 ℃退火30 s,72 ℃延伸6.5 min,25個(gè)循環(huán);72 ℃再延伸10 min。

    用Ⅰ消化PCR模板 (在37 ℃下反應(yīng)1 h),再進(jìn)行純化,獲得帶有突變基因的質(zhì)粒。將構(gòu)建的質(zhì)粒轉(zhuǎn)化感受態(tài)細(xì)胞JM109,涂布含有卡那霉素 (50 mg/L) 的抗性平板,篩選獲得陽(yáng)性克隆,提取質(zhì)粒轉(zhuǎn)化感受態(tài)細(xì)胞BL21 (DE3)。

    1.2.3 突變酶的表達(dá)與純化

    從平板上挑取突變菌種的單菌落,分別接種至5 mL含有卡那霉素 (50 mg/L) 的LB培養(yǎng)基,37 ℃、200 r/min過(guò)夜培養(yǎng)。取1 mL培養(yǎng)液接種至100 mL含有卡那霉素 (50 mg/L) 的2YT培養(yǎng)基中,37 ℃、200 r/min培養(yǎng)至菌濃度達(dá)到600= 0.6,加入終濃度0.6 mmol/L的IPTG,誘導(dǎo)表達(dá)24 h。

    離心收集菌體,用適當(dāng)體積20 mmol/L Tris-HCl (pH 7.5) 重懸后用超聲破碎儀破碎 (工作3 s,停7 s),離心(12 000 ×,4 ℃,30 min)取上清液, 0.45 μm微孔濾膜過(guò)濾,使用HiTrapTM1 mL DEAE FF于AKATA pure層析系統(tǒng)中純化蛋白,純化方法參見文獻(xiàn)[9]。采用Bradford法測(cè)定蛋白質(zhì)濃度[10]。

    1.2.4 CvPAH酶活測(cè)定

    CvPAH酶液稀釋至0.5 mg/mL,取50 μL加到250 μL 100 mmol/L HEPES (pH 7.5) 緩沖液中,然后在該體系中加入50 μL濃度為 10 mmol/L的L-苯丙氨酸 (L-phe)、50 μL濃度為50 mmol/L的DTT、50 μL濃度為50 μmol/L FeSO4、50 μL濃度為2 mmol/L的DMPH4,37 ℃反應(yīng)10 min,加入500 μL甲醇終止反應(yīng)。15 000 ×離心30 min,取上清,用0.22 μm微孔有機(jī)濾膜過(guò)濾,利用HPLC檢測(cè)產(chǎn)物生成量,計(jì)算酶活[9]。CvPAH活性測(cè)定反應(yīng)體系見表2。

    表2 CvPAH酶活測(cè)定反應(yīng)體系

    1.2.5 突變酶最適反應(yīng)溫度及熱穩(wěn)定性測(cè)定

    最適反應(yīng)溫度:根據(jù)表2的反應(yīng)體系,將突變酶分別置于20 ℃、30 ℃、35 ℃、40 ℃、45 ℃、50 ℃、60 ℃條件反應(yīng)10 min,測(cè)定殘余酶活。

    溫度穩(wěn)定性:突變酶在50 ℃水浴中熱處理15 min、30 min、45 min、60 min,立刻冰浴30 min。

    根據(jù)表2配制反應(yīng)體系,迅速振蕩混勻,置于37 ℃水浴鍋中反應(yīng)10 min,測(cè)定殘余酶活。

    1.2.6 突變酶m值測(cè)定

    用圓二色譜儀測(cè)定原始酶及突變酶m。溫度范圍為35–70 ℃,0.1 mg/L蛋白樣溶于 20 mmol/L Tris-HCl (pH 7.5),在222 nm的波長(zhǎng)條件下測(cè)定。用軟件Origin 8.0進(jìn)行平滑處理。

    1.2.7 突變酶動(dòng)力學(xué)參數(shù)測(cè)定

    以不同濃度 (0.2、0.4、0.6、0.8、1.0 mmol/L) L-phe為底物,在37 ℃、pH 7.5的條件下反應(yīng)3 min,測(cè)定突變酶的初始反應(yīng)速率。利用軟件GraphPad Prism 5進(jìn)行非線性擬合,計(jì)算動(dòng)力學(xué)參數(shù) (max、m、cat和catm)。

    1.2.8 突變酶的結(jié)構(gòu)分析

    取0.2 mg/mL蛋白溶解于20 mmol/L Tris-HCl緩沖液(pH 7.5) 中,使用圓二色譜分析儀,掃描波長(zhǎng)從190–250 nm,檢測(cè)蛋白質(zhì)的二級(jí)結(jié)構(gòu)。

    以CvPAH (PDB ID:1LTZ) 為模板,利用SWISS-MODEL對(duì)突變酶進(jìn)行同源建模。利用Pymol軟件對(duì)蛋白質(zhì)的三級(jí)結(jié)構(gòu)進(jìn)行分析。

    2 結(jié)果與分析

    2.1 突變位點(diǎn)保守性分析

    為了比較CvPAH中Gly和Lys的保守性,從動(dòng)物和植物中分別選取人苯丙氨酸羥化酶 (hPAH)和綠藻苯丙氨酸羥化酶 (CrPAH),與CvPAH進(jìn)行序列比對(duì) (圖1)。結(jié)果表明,CvPAH氨基酸序列中有19個(gè)Gly,其中有7個(gè)保守位點(diǎn),分別是92、100、142、160、186、198和200位;有11個(gè)Lys,并且都不保守。

    圖1 紫色色桿菌PAH與人和綠藻PAH氨基酸序列比對(duì)

    2.2 突變體篩選

    用粗酶液初步測(cè)定突變酶相對(duì)酶活和熱穩(wěn)定性。相對(duì)酶活 (圖2中A、B) 表明,Lys突變?yōu)锳rg對(duì)酶活影響較小,除了K86R酶活降至原始酶的45.5%,其余酶的相對(duì)酶活均保持在50%以上。Gly突變成Ala對(duì)CvPAH活性影響較大,其中G100A、G198A和G200A的酶活分別降至原始酶的10.4%、24.4%和28.1%。根據(jù)晶體結(jié)構(gòu) (數(shù)據(jù)未顯示) 分析可知G198和G200均處于活性中心周圍,對(duì)酶的活性部位的微環(huán)境有較大的影響,而G100直接參與輔酶 (6,7-dimethyltetrahydropterin,DMPH4)的結(jié)合,突變后直接影響輔酶結(jié)合,降低酶活。

    圖2 突變酶相對(duì)酶活

    選取Gly突變成Ala后相對(duì)酶活保持在50%以上的突變酶和所有由Lys突變成Arg的突變酶初步測(cè)定熱穩(wěn)定性 (在50 ℃保溫,測(cè)定15 min、30 min時(shí)的酶活保留率)。圖3表明,大部分突變酶的熱穩(wěn)定性下降或影響較小,只有G221A、K94R熱穩(wěn)定性明顯提高。因此,后續(xù)對(duì)G221A、K94R和組合突變酶K94R/G221A展開研究。

    圖3 突變酶熱穩(wěn)定性

    對(duì)各突變位點(diǎn)在晶體結(jié)構(gòu)中所處的位置進(jìn)行比較分析,Gly分布為:19、57、92、100、170、186、193、198、221、272、274、286、290均位于該酶的柔性區(qū),66、142、160、162、163、200均位于α-螺旋中,且19、92、170、193、221、272、274、286、290均位于酶分子表面。Lys分布為:所有點(diǎn)均分布在酶分子表面,且94、167、239位于柔性區(qū),其余位點(diǎn)位于α-螺旋中??梢?,蛋白表面氨基酸性質(zhì)對(duì)酶的熱穩(wěn)定性影響較大,對(duì)分子表面的柔性區(qū)進(jìn)行改造能有效提高突變位點(diǎn)選擇的效率,大大降低工作量。

    2.3 突變酶的酶學(xué)性質(zhì)

    2.3.1 突變酶純化

    突變酶經(jīng)過(guò)DEAE陰離子柱純化,進(jìn)行SDS-PAGE電泳檢測(cè) (圖4),突變酶的分子量大小與原始酶一致,均在35 kDa左右。

    圖4 突變酶純化結(jié)果

    2.3.2 突變酶的最適反應(yīng)溫度及熱穩(wěn)定性

    突變酶的活性研究結(jié)果表明,K94R、G221A和K94R/G221A的比酶活分別為473.8、413.6和294.3 U/mg,與原始酶 (562.7 U/mg) 相比,分別下降了15.8%、26.5%和47.7%。最適反應(yīng)溫度結(jié)果 (圖5) 表明,突變酶和原始酶均在40 ℃左右達(dá)到最高酶活,隨著溫度升高或降低,酶活均會(huì)降低,其中K94R/G221A在溫度降低過(guò)程中酶活降低速度最慢。

    圖5 突變酶最適反應(yīng)溫度

    熱穩(wěn)定性對(duì)比分析(圖6) 顯示,在50 ℃處理時(shí),原始酶的酶活下降速度明顯比突變酶快,尤其在前15 min,原始酶活性迅速喪失,最終只剩25.1%,而突變酶K94R、G221A和K94R/G221A酶活保留率分別達(dá)到65.5%、53.4%和80.2%。K94R和G221A在50 ℃的半衰期分別為26.2 min、16.8 min,比原始酶 (9.0 min)分別提高了1.9倍、0.9倍,同時(shí)K94R/G221A在50 ℃處理1 h后仍保留65.6%的酶活,比原始酶 (8.6%) 高出6.6倍。

    圖6 突變酶熱穩(wěn)定性比較

    2.3.3 突變酶的m值

    通過(guò)圓二色譜儀對(duì)突變酶K94R、G221A、K94R/G221A及原始酶進(jìn)行m測(cè)定,結(jié)果如圖7所示,原始酶的m值約51.5 ℃,而突變酶K94R、G221A及K94R/G221A分別達(dá)到53.8 ℃、53.1 ℃、54.8 ℃。因此,相比較原始酶,突變酶具有更高的耐熱性。

    圖7 圓二色譜測(cè)定突變酶Tm值

    2.3.4 突變酶的動(dòng)力學(xué)參數(shù)

    以不同濃度的L-phe為底物,測(cè)定突變酶的動(dòng)力學(xué)參數(shù)。如表3所示,與原始酶相比,突變后最大反應(yīng)速率max有所降低,其中K94R/G221A的max最低,降到5.8 μmol/(L·s)。突變酶對(duì)底物的米氏常數(shù)m均有不同程度下降,其中K94R對(duì)L-phe的m最低,為119.3 μmol/L,可見突變酶對(duì)底物的結(jié)合能力均有所提高,同時(shí)cat也有不同程度的下降。最終,催化效率 (catm) 僅略低于原始酶,即該突變對(duì)酶活的影響較小。結(jié)合穩(wěn)定性研究結(jié)果可知,相比原酶,突變酶更具有應(yīng)用價(jià)值。

    表3 突變酶動(dòng)力學(xué)參數(shù)

    2.4 突變酶的結(jié)構(gòu)分析

    二級(jí)結(jié)構(gòu)分析結(jié)果如圖8。突變酶與原始酶α-螺旋特征波長(zhǎng)208 nm和222 nm以及β-折疊特征波長(zhǎng)195 nm和216 nm處的橢圓率沒有顯著差別??梢奒94R、G221A兩個(gè)位點(diǎn)突變對(duì)酶的二級(jí)結(jié)構(gòu)幾乎沒有影響。

    圖8 突變酶二級(jí)結(jié)構(gòu)

    同源建模發(fā)現(xiàn),突變酶與原始酶的三級(jí)結(jié)構(gòu)沒有顯著差別,因此對(duì)氨基酸之間成鍵情況進(jìn)行分析。如圖9所示,94位氨基酸位于α-螺旋與β-折疊之間的柔性環(huán)上,由Lys突變成帶有更大側(cè)鏈基團(tuán)的Arg,與Ile95形成額外的氫鍵,有助于穩(wěn)定該環(huán)的結(jié)構(gòu),從而提高酶的穩(wěn)定性;此外,該位點(diǎn)位于蛋白質(zhì)表面區(qū)域,突變成Arg后胍基能夠與水分子之間形成偶極-電荷作用,使蛋白表面形成水化膜,從而提高蛋白質(zhì)熱穩(wěn)定性[24]。

    圖9 位點(diǎn)94突變前(A) 后(B) 三維結(jié)構(gòu)模擬

    如圖10所示,Gly221處于柔性區(qū)并且與酶分子的C-端柔性區(qū)接近,Gly沒有側(cè)鏈基團(tuán),構(gòu)象熵大。突變成Ala后,增加了一個(gè)疏水性基團(tuán) (甲基),拉近了與酶C-端Leu281的距離,與Leu281疏水側(cè)鏈出現(xiàn)重疊區(qū)域 (圖10),推測(cè)形成疏水作用[25–26],不僅增加了柔性環(huán)的穩(wěn)定性,同時(shí)起到錨定C-端柔性區(qū)的作用,提高了突變酶的熱穩(wěn)定性。

    圖10 位點(diǎn)22突變前(A) 后(B)三維結(jié)構(gòu)模擬

    3 結(jié)論

    通過(guò)Gly替換成Ala以及Lys替換成Arg的策略得到2個(gè)突變酶K94R和G221A,其熱穩(wěn)定性較原始酶分別提高了0.9倍和1.9倍。進(jìn)一步對(duì)兩個(gè)位點(diǎn)進(jìn)行組合突變,獲得K94R/ G221A,其熱穩(wěn)定性較原始酶提高了6.6倍。圓二色譜數(shù)據(jù)進(jìn)一步證明突變酶的具有更高的耐熱性。動(dòng)力學(xué)分析結(jié)果表明,突變酶對(duì)底物的親和力提高了,同時(shí)催化效率僅略低于原始酶。最后,通過(guò)蛋白三級(jí)結(jié)構(gòu)模擬發(fā)現(xiàn),K94R和G221A分別通過(guò)形成氫鍵和增強(qiáng)局部疏水作用提高了酶的熱穩(wěn)定性。

    REFERENCES

    [1] Flydal MI, Martinez A. Phenylalanine hydroxylase: function, structure, and regulation. IUBMB Life, 2013, 65(4): 341–349.

    [2] Williams RA, Mamotte CDS, Burnett JR. Phenylketonuria: an inborn error of phenylalanine metabolism. Clin Biochem Rev, 2008, 29(1): 31–41.

    [3] Przyrembel H, Bremer HJ. Nutrition, physical growth, and bone density in treated phenylketonuria. Eur J Pediatr, 2000, 159(2 Supl): S129–S135.

    [4] Chanoine JP. Selenium and thyroid function in infants, children and adolescents. Biofactors, 2003, 19(3/4): 137–143.

    [5] Beblo S, Reinhardt H, Demmelmair H, et al. Effect of fish oil supplementation on fatty acid status, coordination, and fine motor skills in children with phenylketonuria. J Pediatr, 2007, 150(5): 479–484.

    [6] Colomé C, Artuch R, Vilaseca MA, et al. Ubiquinone-10 content in lymphocytes of phenylketonuric patients. Clin Biochem, 2002, 35(1): 81–84.

    [7] Nakata H, Yamauchi T, Fujisawa H. Phenylalanine hydroxylase frompurification and characterization. J Biol Chem, 1979, 254(6): 1829–1833.

    [8] Onishi A, Liotta LJ, Benkovic SJ. Cloning and expression ofphenylalanine hydroxylase inand comparison of amino acid sequence with mammalian aromatic amino acid hydroxylases. J Biol Chem, 1991, 266(28): 18454–18459.

    [9] Cao SH, Zhou L, Cui WJ, et al. Heterologous expression and characterization of phenylalanine hydroxylase from. Biotechnol Bull, 2014, (8): 152–158 (in Chinese). 曹淑慧, 周麗, 崔文璟, 等. 紫色色桿菌苯丙氨酸羥化酶的異源表達(dá)及重組酶學(xué)性質(zhì)研究. 生物技術(shù)通報(bào), 2014, (8): 152–158.

    [10] Yew NS, Dufour E, Przybylska M, et al. Erythrocytes encapsulated with phenylalanine hydroxylase exhibit improved pharmacokinetics and lowered plasma phenylalanine levels in normal mice. Mol Genet Metab, 2013, 109(4): 339–344.

    [11] Turner NJ. Directed evolution drives the next generation of biocatalysts. Nat Chem Biol, 2009, 5(8): 56–573.

    [12] Reetz MT, Carballeira JD, Vogel A. Iterative saturation mutagenesis on the basis of B factors as a strategy for increasing protein thermostability. Angew Chem Int EdEngl, 2006, 45(46): 7745–7751.

    [13] Polizzi KM, Bommarius AS, Broering JM, et al. Stability of biocatalysts. Curr Opin Chem Biol, 2007, 11(2): 220–225.

    [14] Martin A, Kather I, Schmid FX. Origins of the high stability of an-selected cold-shock protein. J Mol Biol, 2002, 318(5): 1341–1349.

    [15] Bagautdinov B, Yutani K. Structure of indole-3-glycerol phosphate synthase fromHB8: implication for thermal stability. Acta Crystallogr D, 2011, 67(12): 1054–1064.

    [16] Chakravarty S, Varadarajan R. Elucidation of factors responsible for enhanced thermal stability of proteins: a structural genomics based study. Biochemistry, 2002, 41(25): 8152–8161.

    [17] Delboni LF, Mande SC, Rentier-Delrue F, et al. Crystal structure of recombinant triosephosphate isomerase from: an analysis of potential thermostability factors in 6 isomerases with known 3-dimensional structures points to the importance of hydrophobic interactions. Protein Sci, 1995, 4(12): 2594–2604.

    [18] Vogt G, Woell S, Argos P. Protein thermal stability, hydrogen bonds, and ion pairs. J Mol Biol, 1997, 269(4): 631–643.

    [19] Ma YF, Eglinton JK, Evans DE, et al. Removal of the four C-terminal glycine-rich repeats enhances the thermostability and substrate binding affinity of barley β-amylase. Biochemistry, 2000, 39(44): 13350–13355.

    [20] Mahanta P, Bhardwaj A, Kumar K, et al. Structural insights into N-terminal to C-terminal interactions and implications for thermostability of a (β/α)8-triosephosphate isomerase barrel enzyme. FEBS J, 2015, 282(18): 3543–3555.

    [21] Argos P, Rossman MG, Grau UM, et al. Thermal stability and protein structure. Biochemistry, 1979, 18(25): 5698–5703.

    [22] Zhang XJ, Baase WA, Matthews BW. Multiple alanine replacements within α-helix 126–134 of T4 lysozyme have independent, additive effects on both structure and stability. Protein Sci, 1992, 1(6): 761–776.

    [23] Mrabet NT, Van den Broeck A, Van den Brande I, et al. Arginine residues as stabilizing elements in proteins. Biochemistry, 1992, 31(8): 2239–2253.

    [24] Rahman RNZA, Fujiwara S, Nakamura H, et al. Ion pairs involved in maintaining a thermostable structure of glutamate dehydrogenase from a hyperthermophilic archaeon. Biochem Bioph Res Co, 1998, 248(3): 920–926.

    [25] Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem, 1976, 72(1/2): 248–254.

    [26] Martin A, Sieber V, Schmid FX.-selection of highly stabilized protein variants with optimized surface. J Mol Biol, 2001, 309(3): 717–726.

    (本文責(zé)編 郝麗芳)

    Thermal stability improvement for phenylalanine hydroxylase by site-directed mutagenesis

    Shuangshuang Ye, Li Zhou, and Zhemin Zhou

    School of Biotechnology, Jiangnan University, Wuxi 214122, Jiangsu, China

    In proteins of thermophilic bacteria, Gly is tend to be replaced by Ala and Lys is tend to be replaced by Arg to adapt the high temperature. In order to improve the thermal stability of phenylalanine hydroxylase (PAH) from, all the Gly on PAH were mutated to Ala and Lys to Arg. Positive mutant enzymes with improved thermal stability were selected, followed by combined mutation and characterization. The results revealed that half-lives of K94R and G221A mutants at 50 °C were 26.2 min and 16.8 min, which were increased by 1.9-times and 0.9-times than the parent enzyme (9.0 min). The residual activity of K94R/G221A mutant was improved to 65.6% after keeping at 50 °C for 1 h, which was 6.6 time higher than the parent enzyme (8.6%). Circular dichroism (CD) spectroscopy revealed thatmvalues of the parent enzyme, K94R, G221A and K94R/G221A were 51.5 ℃, 53.8 ℃, 53.1 ℃ and 54.8 ℃, respectively. According to the protein structure simulation, the two mutations were located on flexible loop. In the K94R mutant, the mutated Arg94 on the surface of the enzyme formed an extra hydrogen bond with Ile95, which stabilized the located loop. In the G221A mutant, the mutated Ala221 formed hydrophobic interaction with Leu281, which could stabilize both the loop and flexible area of the C-terminus of G221A. The results not only provided a reference for protein modification on thermal stability, but also laid the foundation for application of phenylalanine hydroxylase in the field of functional foods.

    phenylalanine hydroxylase, thermal stability, site-directed mutagenesis, alanine mutation, arginine mutation

    January 12, 2016; Accepted:May 30, 2016

    Zhemin Zhou. Tel: +86-510-85325210; E-mail: zhmzhou@jiangnan.edu.cn

    Supported by: National High Technology Research and Development Program of China (863 Program) (No. 2014AA021304),the Priority Academic Program Development of Jiangsu Higher-Education Institutions, the 111 Project, the Jiangsu Province “Collaborative Innovation Center for Advanced Industrial Fermentation” Industry Development Program (No. 111-2-06),The Fundamental Research Funds for the Central Universities (No. JUSRP51411B).

    國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃) (No. 2014AA021304),江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程自主項(xiàng)目,111引智計(jì)劃 (No. 111-2-06),江蘇省現(xiàn)代發(fā)酵工業(yè)協(xié)同創(chuàng)新中心,2014年基本科研-重點(diǎn)項(xiàng)目子項(xiàng)目 (No. JUSRP51411B) 資助。

    網(wǎng)絡(luò)出版時(shí)間:2016-06-13 網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/11.1998.Q.20160613.0956.002.html

    葉雙雙, 周麗, 周哲敏. 定點(diǎn)突變提高苯丙氨酸羥化酶的熱穩(wěn)定性. 生物工程學(xué)報(bào), 2016, 32(9): 1243–1254.

    Ye SS, Zhou L, Zhou ZM. Thermal stability improvement for phenylalanine hydroxylase by site-directed mutagenesis. Chin J Biotech, 2016, 32(9): 1243–1254.

    猜你喜歡
    羥化酶苯丙氨酸熱穩(wěn)定性
    CYP17A1基因His373Asn純合突變的17α-羥化酶/17,20碳鏈裂解酶缺陷癥合并糖尿病1例臨床分析
    2009~2019 年吉林省新生兒高苯丙氨酸血癥的發(fā)病率及治療效果分析
    橡膠樹F3’H基因克隆及其功能分析
    GEMOX方案聯(lián)合恩度治療肝癌的遠(yuǎn)期療效及對(duì)血清MBL、ASPH水平的影響
    淺談17-羥化酶缺乏癥及MDT治療
    PVC用酪氨酸鑭的合成、復(fù)配及熱穩(wěn)定性能研究
    酮苯丙氨酸鈣原料藥有關(guān)物質(zhì)的鑒定與含量測(cè)定
    提高有機(jī)過(guò)氧化物熱穩(wěn)定性的方法
    提高有機(jī)過(guò)氧化物熱穩(wěn)定性的方法
    可聚合松香衍生物的合成、表征和熱穩(wěn)定性?
    国产有黄有色有爽视频| 国产一区二区在线观看av| 国产极品粉嫩免费观看在线| 亚洲熟女精品中文字幕| 少妇人妻久久综合中文| 91精品国产国语对白视频| 少妇熟女欧美另类| 尾随美女入室| 国产免费一区二区三区四区乱码| 如何舔出高潮| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久久久久亚洲中文字幕| 99视频精品全部免费 在线| 五月伊人婷婷丁香| 多毛熟女@视频| 又黄又粗又硬又大视频| 青春草视频在线免费观看| 精品国产露脸久久av麻豆| 少妇被粗大猛烈的视频| 在线观看美女被高潮喷水网站| 最近的中文字幕免费完整| 一级爰片在线观看| 国产免费一级a男人的天堂| 日本wwww免费看| 91在线精品国自产拍蜜月| 人妻系列 视频| 巨乳人妻的诱惑在线观看| 国产成人午夜福利电影在线观看| 久久久久久久久久人人人人人人| 亚洲成国产人片在线观看| 青春草亚洲视频在线观看| 日韩免费高清中文字幕av| 国产精品嫩草影院av在线观看| av有码第一页| 国产一区二区激情短视频 | 久久国产亚洲av麻豆专区| 一边摸一边做爽爽视频免费| 极品少妇高潮喷水抽搐| 99久久人妻综合| 男人舔女人的私密视频| 内地一区二区视频在线| 国产xxxxx性猛交| 人体艺术视频欧美日本| av网站免费在线观看视频| 成人毛片a级毛片在线播放| 国产成人精品无人区| 中文天堂在线官网| 亚洲精品aⅴ在线观看| 大陆偷拍与自拍| 国产精品.久久久| 制服丝袜香蕉在线| 久久99蜜桃精品久久| 久久久久久伊人网av| 中文精品一卡2卡3卡4更新| 国产精品女同一区二区软件| 好男人视频免费观看在线| 丝袜脚勾引网站| 欧美 日韩 精品 国产| h视频一区二区三区| 国产一区二区三区综合在线观看 | 久久午夜综合久久蜜桃| 国产一区二区在线观看av| 少妇精品久久久久久久| 国产探花极品一区二区| 国产免费福利视频在线观看| 色婷婷av一区二区三区视频| 黄色怎么调成土黄色| 国产免费一区二区三区四区乱码| 韩国高清视频一区二区三区| 午夜影院在线不卡| 在线天堂中文资源库| 国产精品国产三级专区第一集| 777米奇影视久久| 久久97久久精品| 国产精品一区www在线观看| 国产爽快片一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产免费视频播放在线视频| 精品久久久精品久久久| 在线观看一区二区三区激情| 久久久亚洲精品成人影院| 97在线视频观看| 久久久久国产精品人妻一区二区| 亚洲成av片中文字幕在线观看 | 亚洲伊人久久精品综合| 人妻人人澡人人爽人人| 成年动漫av网址| 22中文网久久字幕| 久久热在线av| 老熟女久久久| 22中文网久久字幕| 国产综合精华液| 18禁动态无遮挡网站| 国产免费现黄频在线看| 一本—道久久a久久精品蜜桃钙片| 久久国内精品自在自线图片| 国产淫语在线视频| 国产一区二区三区av在线| 久久97久久精品| 欧美成人午夜免费资源| 九九在线视频观看精品| 夫妻午夜视频| 欧美成人午夜精品| 90打野战视频偷拍视频| 久久久久视频综合| 亚洲人与动物交配视频| 亚洲中文av在线| 精品一区二区三区四区五区乱码 | 国产成人精品福利久久| 日韩伦理黄色片| 精品卡一卡二卡四卡免费| 91成人精品电影| 亚洲欧美一区二区三区黑人 | 九色亚洲精品在线播放| 在线天堂中文资源库| 亚洲国产av影院在线观看| 国产在视频线精品| 18禁国产床啪视频网站| 搡女人真爽免费视频火全软件| 最后的刺客免费高清国语| 国产亚洲最大av| 97精品久久久久久久久久精品| 色5月婷婷丁香| 欧美日韩精品成人综合77777| 日日爽夜夜爽网站| 美女中出高潮动态图| 天天操日日干夜夜撸| 色哟哟·www| 麻豆精品久久久久久蜜桃| 久久鲁丝午夜福利片| 亚洲国产av影院在线观看| 国内精品宾馆在线| 亚洲,欧美精品.| 99香蕉大伊视频| 少妇猛男粗大的猛烈进出视频| 久久久国产欧美日韩av| 国产免费一级a男人的天堂| 搡老乐熟女国产| 亚洲精品456在线播放app| 人人妻人人澡人人爽人人夜夜| 成人国语在线视频| 三级国产精品片| 国产在线免费精品| 国产亚洲午夜精品一区二区久久| 九色亚洲精品在线播放| 国产成人精品在线电影| av一本久久久久| 亚洲av综合色区一区| 波野结衣二区三区在线| 久久久国产一区二区| 日韩一区二区三区影片| 免费高清在线观看日韩| 中国国产av一级| 欧美 日韩 精品 国产| 插逼视频在线观看| 免费av中文字幕在线| 在线观看www视频免费| 国产精品.久久久| 国产男人的电影天堂91| 晚上一个人看的免费电影| 人体艺术视频欧美日本| 在线观看美女被高潮喷水网站| 99久久人妻综合| 成人毛片a级毛片在线播放| 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| 两性夫妻黄色片 | 成年人免费黄色播放视频| 美女国产视频在线观看| a级毛色黄片| 日韩熟女老妇一区二区性免费视频| 国产精品人妻久久久影院| 在线免费观看不下载黄p国产| 亚洲精品美女久久久久99蜜臀 | av视频免费观看在线观看| 婷婷色综合大香蕉| 人妻少妇偷人精品九色| 波野结衣二区三区在线| 久久综合国产亚洲精品| 天天影视国产精品| 少妇被粗大猛烈的视频| 在线免费观看不下载黄p国产| 大香蕉久久网| av一本久久久久| 九九爱精品视频在线观看| 国产成人精品无人区| 人人妻人人添人人爽欧美一区卜| 美女中出高潮动态图| 波多野结衣一区麻豆| 国语对白做爰xxxⅹ性视频网站| 久久这里有精品视频免费| 99香蕉大伊视频| 日本欧美国产在线视频| 中国美白少妇内射xxxbb| 亚洲av成人精品一二三区| 欧美精品亚洲一区二区| 黄色一级大片看看| 亚洲内射少妇av| 美女视频免费永久观看网站| 免费女性裸体啪啪无遮挡网站| 成人影院久久| 国产成人精品一,二区| 日本与韩国留学比较| 日韩欧美精品免费久久| 国产午夜精品一二区理论片| 国产极品粉嫩免费观看在线| 国产 精品1| 免费大片黄手机在线观看| 又粗又硬又长又爽又黄的视频| 热re99久久精品国产66热6| 99国产综合亚洲精品| 国产精品一区二区在线不卡| 国产xxxxx性猛交| 69精品国产乱码久久久| 免费黄色在线免费观看| 国产福利在线免费观看视频| 国产1区2区3区精品| 久久精品熟女亚洲av麻豆精品| 波多野结衣一区麻豆| 女性被躁到高潮视频| 在线天堂中文资源库| 国产精品免费大片| 丝袜喷水一区| 久久久久精品久久久久真实原创| 欧美xxxx性猛交bbbb| videosex国产| 精品少妇久久久久久888优播| 欧美精品高潮呻吟av久久| 国产精品成人在线| 国产 精品1| 90打野战视频偷拍视频| 国产高清国产精品国产三级| 国产xxxxx性猛交| 精品国产乱码久久久久久小说| 黄色视频在线播放观看不卡| 日本-黄色视频高清免费观看| 天堂8中文在线网| 成人综合一区亚洲| 少妇人妻 视频| 亚洲五月色婷婷综合| 国产男女超爽视频在线观看| 亚洲精品av麻豆狂野| 两个人免费观看高清视频| 亚洲国产av影院在线观看| a级毛色黄片| 欧美xxⅹ黑人| 久久午夜福利片| 国产伦理片在线播放av一区| 精品久久国产蜜桃| 午夜影院在线不卡| 我的女老师完整版在线观看| 丝瓜视频免费看黄片| 丝袜人妻中文字幕| 婷婷色综合www| 久久鲁丝午夜福利片| 97在线视频观看| 亚洲少妇的诱惑av| 一区二区av电影网| 国产成人欧美| 黄片无遮挡物在线观看| 少妇 在线观看| 国产有黄有色有爽视频| 9191精品国产免费久久| 大香蕉久久网| 免费观看性生交大片5| 在线观看免费视频网站a站| 99香蕉大伊视频| 99热6这里只有精品| 免费女性裸体啪啪无遮挡网站| 一区二区三区四区激情视频| 伦精品一区二区三区| 啦啦啦在线观看免费高清www| 国产欧美日韩一区二区三区在线| 制服诱惑二区| 麻豆精品久久久久久蜜桃| 国产免费福利视频在线观看| 我要看黄色一级片免费的| 18禁国产床啪视频网站| 五月玫瑰六月丁香| 国产日韩欧美亚洲二区| 欧美人与性动交α欧美精品济南到 | 蜜臀久久99精品久久宅男| 国产精品熟女久久久久浪| 久久精品人人爽人人爽视色| 亚洲情色 制服丝袜| 精品一区二区三卡| 亚洲成国产人片在线观看| 国产爽快片一区二区三区| 一区二区三区精品91| av不卡在线播放| 亚洲成色77777| 精品第一国产精品| 久久这里有精品视频免费| 一本色道久久久久久精品综合| 香蕉精品网在线| 亚洲av中文av极速乱| 欧美人与性动交α欧美软件 | 亚洲精品日本国产第一区| 久久女婷五月综合色啪小说| 日本爱情动作片www.在线观看| 久久毛片免费看一区二区三区| 亚洲国产欧美在线一区| 极品少妇高潮喷水抽搐| 亚洲综合色惰| 精品人妻一区二区三区麻豆| 亚洲国产av新网站| 日韩制服丝袜自拍偷拍| 日韩精品有码人妻一区| 免费看av在线观看网站| 黄片无遮挡物在线观看| 另类精品久久| 国产黄频视频在线观看| 韩国高清视频一区二区三区| 国产精品成人在线| 亚洲成人一二三区av| 18在线观看网站| 午夜激情av网站| 中文字幕另类日韩欧美亚洲嫩草| 夜夜骑夜夜射夜夜干| 国产精品.久久久| 亚洲精品久久午夜乱码| kizo精华| 亚洲国产精品一区三区| 日韩成人伦理影院| 久久久久人妻精品一区果冻| 免费看不卡的av| 2018国产大陆天天弄谢| 狠狠精品人妻久久久久久综合| 亚洲美女黄色视频免费看| 欧美精品亚洲一区二区| 日韩一本色道免费dvd| 老司机影院成人| 美女脱内裤让男人舔精品视频| 精品少妇内射三级| 色网站视频免费| 热re99久久国产66热| 免费大片黄手机在线观看| 日韩成人av中文字幕在线观看| 蜜桃在线观看..| √禁漫天堂资源中文www| 亚洲第一av免费看| 777米奇影视久久| 国产高清三级在线| 久久久久视频综合| 午夜影院在线不卡| 中文精品一卡2卡3卡4更新| 国产色爽女视频免费观看| 亚洲av成人精品一二三区| 丝袜在线中文字幕| 国产极品天堂在线| 久久99热6这里只有精品| 国产精品 国内视频| 国产老妇伦熟女老妇高清| 色网站视频免费| 久久韩国三级中文字幕| 亚洲av电影在线进入| 一二三四中文在线观看免费高清| 欧美激情极品国产一区二区三区 | 久久久久久久久久人人人人人人| 大片免费播放器 马上看| 少妇高潮的动态图| 久久午夜福利片| 国产男女超爽视频在线观看| 婷婷色综合www| 香蕉丝袜av| 一本久久精品| 久久99一区二区三区| 大香蕉97超碰在线| 男的添女的下面高潮视频| 久久久久精品性色| 91午夜精品亚洲一区二区三区| 九色亚洲精品在线播放| 免费人妻精品一区二区三区视频| 国产成人欧美| 一区二区av电影网| av在线老鸭窝| 免费久久久久久久精品成人欧美视频 | 成年美女黄网站色视频大全免费| 美女国产高潮福利片在线看| 亚洲精品日韩在线中文字幕| 黄网站色视频无遮挡免费观看| 这个男人来自地球电影免费观看 | av女优亚洲男人天堂| 天天操日日干夜夜撸| 国产精品一区www在线观看| 久久久久久久久久久久大奶| 97精品久久久久久久久久精品| 1024视频免费在线观看| 建设人人有责人人尽责人人享有的| 日产精品乱码卡一卡2卡三| 日韩一区二区三区影片| 哪个播放器可以免费观看大片| 如何舔出高潮| 少妇人妻久久综合中文| 国产xxxxx性猛交| 国产激情久久老熟女| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 一区二区日韩欧美中文字幕 | 丁香六月天网| 少妇的逼好多水| 考比视频在线观看| 综合色丁香网| 国产精品一区二区在线观看99| 丝袜美足系列| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产精品成人久久小说| 欧美亚洲 丝袜 人妻 在线| 久久精品久久久久久久性| 欧美精品一区二区免费开放| 亚洲国产成人一精品久久久| av免费在线看不卡| 精品一品国产午夜福利视频| 欧美性感艳星| 国产白丝娇喘喷水9色精品| 亚洲成人手机| 久久人妻熟女aⅴ| 曰老女人黄片| 日产精品乱码卡一卡2卡三| 婷婷色av中文字幕| 亚洲欧洲国产日韩| 日本91视频免费播放| 中文精品一卡2卡3卡4更新| 午夜福利乱码中文字幕| 天堂8中文在线网| 国产在线免费精品| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| a级毛色黄片| 美女主播在线视频| 老女人水多毛片| 高清av免费在线| 国产av一区二区精品久久| 久久影院123| 99re6热这里在线精品视频| 高清不卡的av网站| 99精国产麻豆久久婷婷| 亚洲欧美中文字幕日韩二区| 超色免费av| 18禁动态无遮挡网站| 日本vs欧美在线观看视频| 国产色爽女视频免费观看| 亚洲内射少妇av| 日韩成人伦理影院| av卡一久久| 国产精品久久久久久久电影| 最近的中文字幕免费完整| 制服诱惑二区| 精品亚洲成a人片在线观看| 中文精品一卡2卡3卡4更新| 少妇的逼水好多| 久久综合国产亚洲精品| 国产男女超爽视频在线观看| 人妻少妇偷人精品九色| 亚洲天堂av无毛| 97超碰精品成人国产| 亚洲av电影在线观看一区二区三区| 久久鲁丝午夜福利片| 日韩欧美精品免费久久| 久久久亚洲精品成人影院| 热re99久久精品国产66热6| 人人妻人人澡人人看| 国产亚洲av片在线观看秒播厂| 午夜激情久久久久久久| 午夜福利,免费看| av卡一久久| 在线看a的网站| 韩国av在线不卡| av福利片在线| 久久国产精品男人的天堂亚洲 | 波野结衣二区三区在线| 丝袜脚勾引网站| 中文欧美无线码| 伊人久久国产一区二区| √禁漫天堂资源中文www| 精品少妇黑人巨大在线播放| 丰满少妇做爰视频| a级毛片在线看网站| 国产精品国产三级国产专区5o| freevideosex欧美| 侵犯人妻中文字幕一二三四区| 一边摸一边做爽爽视频免费| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久av不卡| 国产在线免费精品| 日本爱情动作片www.在线观看| 亚洲精华国产精华液的使用体验| 国国产精品蜜臀av免费| 国产亚洲精品第一综合不卡 | 亚洲人与动物交配视频| 美女大奶头黄色视频| 成年人午夜在线观看视频| 男女免费视频国产| 大香蕉97超碰在线| 亚洲成人av在线免费| 国产免费一区二区三区四区乱码| 女性生殖器流出的白浆| 国产精品久久久久久av不卡| 90打野战视频偷拍视频| 国产免费一区二区三区四区乱码| 国产成人精品在线电影| 精品国产国语对白av| 久久女婷五月综合色啪小说| 毛片一级片免费看久久久久| 看十八女毛片水多多多| 久久韩国三级中文字幕| 成人免费观看视频高清| 亚洲成人手机| 99精国产麻豆久久婷婷| 夜夜骑夜夜射夜夜干| 亚洲第一区二区三区不卡| 99久久精品国产国产毛片| 全区人妻精品视频| 久久99热6这里只有精品| 色婷婷av一区二区三区视频| av视频免费观看在线观看| 在线观看免费高清a一片| 全区人妻精品视频| 欧美3d第一页| 亚洲一级一片aⅴ在线观看| 国产精品不卡视频一区二区| 大片免费播放器 马上看| a级毛片在线看网站| 免费黄网站久久成人精品| 国产国拍精品亚洲av在线观看| 亚洲,欧美,日韩| 性色avwww在线观看| 蜜臀久久99精品久久宅男| 高清在线视频一区二区三区| 欧美精品av麻豆av| 九九在线视频观看精品| 久久精品久久久久久久性| 一本—道久久a久久精品蜜桃钙片| 国产伦理片在线播放av一区| 中文字幕精品免费在线观看视频 | 97在线人人人人妻| 久久午夜综合久久蜜桃| 香蕉丝袜av| 啦啦啦啦在线视频资源| 亚洲一级一片aⅴ在线观看| 久久av网站| 国产在视频线精品| 制服诱惑二区| 亚洲国产看品久久| 日日摸夜夜添夜夜爱| av网站免费在线观看视频| 久久久久久人妻| 久久久精品94久久精品| 黑丝袜美女国产一区| 99热这里只有是精品在线观看| 成人国语在线视频| 亚洲在久久综合| 亚洲av免费高清在线观看| 国产精品无大码| 国产又爽黄色视频| 91在线精品国自产拍蜜月| 永久免费av网站大全| 国产爽快片一区二区三区| 黑人猛操日本美女一级片| 老司机影院毛片| 天堂俺去俺来也www色官网| 日韩三级伦理在线观看| 捣出白浆h1v1| 丝袜脚勾引网站| 久久这里只有精品19| 亚洲欧美成人精品一区二区| 美国免费a级毛片| 91aial.com中文字幕在线观看| 在线观看一区二区三区激情| 亚洲,一卡二卡三卡| 成人毛片a级毛片在线播放| 成年av动漫网址| 亚洲欧洲精品一区二区精品久久久 | 99久国产av精品国产电影| 亚洲色图 男人天堂 中文字幕 | 国产精品秋霞免费鲁丝片| 在线观看免费视频网站a站| 亚洲成人手机| 久久99蜜桃精品久久| 精品一区二区三卡| 精品少妇内射三级| 国产黄频视频在线观看| 亚洲国产看品久久| 王馨瑶露胸无遮挡在线观看| 91国产中文字幕| 精品国产一区二区三区四区第35| 精品少妇内射三级| 成人手机av| 国产xxxxx性猛交| 91精品伊人久久大香线蕉| 婷婷色麻豆天堂久久| 国产精品麻豆人妻色哟哟久久| 日本av手机在线免费观看| 成人毛片60女人毛片免费| 亚洲图色成人| 激情五月婷婷亚洲| 99视频精品全部免费 在线| 国产成人欧美| 男女高潮啪啪啪动态图| 免费观看在线日韩| 久久久久久久大尺度免费视频| www.熟女人妻精品国产 | 麻豆乱淫一区二区| 在线观看三级黄色| 欧美人与善性xxx| 99九九在线精品视频| 久久久久久久久久久免费av| 久久精品国产自在天天线| 天天躁夜夜躁狠狠躁躁| av免费在线看不卡| 久久这里有精品视频免费| 观看美女的网站| 中文字幕另类日韩欧美亚洲嫩草| 日韩一区二区三区影片| 精品久久久久久电影网|