• <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)定性?
    久久精品久久久久久久性| 欧美亚洲 丝袜 人妻 在线| 纯流量卡能插随身wifi吗| 久久久久精品国产欧美久久久 | 久久久久久久国产电影| av电影中文网址| 成人国产av品久久久| 99久国产av精品国产电影| 久久久久视频综合| 纵有疾风起免费观看全集完整版| 亚洲国产精品999| 男女无遮挡免费网站观看| 亚洲国产欧美日韩在线播放| 高清在线视频一区二区三区| 亚洲在久久综合| 国产男女超爽视频在线观看| 叶爱在线成人免费视频播放| 日韩精品有码人妻一区| 多毛熟女@视频| 久久 成人 亚洲| 亚洲欧美一区二区三区久久| 久久精品久久久久久久性| 亚洲精品一二三| 亚洲精品一二三| 日本91视频免费播放| 老汉色∧v一级毛片| 日韩欧美精品免费久久| 一边亲一边摸免费视频| 国产一区二区三区av在线| 99精国产麻豆久久婷婷| 高清在线视频一区二区三区| 天堂俺去俺来也www色官网| 日韩成人av中文字幕在线观看| 中国三级夫妇交换| 亚洲国产欧美一区二区综合| 欧美日本中文国产一区发布| 亚洲专区中文字幕在线 | 涩涩av久久男人的天堂| a 毛片基地| 亚洲精品中文字幕在线视频| 男人添女人高潮全过程视频| 人体艺术视频欧美日本| 日韩中文字幕视频在线看片| 午夜日韩欧美国产| 国产探花极品一区二区| 精品亚洲乱码少妇综合久久| 在现免费观看毛片| www.精华液| 桃花免费在线播放| 国产女主播在线喷水免费视频网站| 丝袜美腿诱惑在线| 亚洲成人av在线免费| 精品国产一区二区三区久久久樱花| 最近2019中文字幕mv第一页| 一区二区三区乱码不卡18| 一级毛片我不卡| 久久热在线av| 最黄视频免费看| 操出白浆在线播放| 国产97色在线日韩免费| 成人国语在线视频| 777米奇影视久久| 日韩伦理黄色片| 如日韩欧美国产精品一区二区三区| 在线看a的网站| 不卡视频在线观看欧美| 一本一本久久a久久精品综合妖精| 777久久人妻少妇嫩草av网站| 国产乱人偷精品视频| 欧美日韩一区二区视频在线观看视频在线| 侵犯人妻中文字幕一二三四区| 日韩不卡一区二区三区视频在线| 国产黄色免费在线视频| e午夜精品久久久久久久| 高清av免费在线| 精品一区在线观看国产| 看免费av毛片| 欧美日韩视频高清一区二区三区二| 菩萨蛮人人尽说江南好唐韦庄| 一级片'在线观看视频| 无遮挡黄片免费观看| 亚洲精品国产区一区二| 乱人伦中国视频| 色播在线永久视频| 女性被躁到高潮视频| 亚洲国产最新在线播放| 美女福利国产在线| 国产精品女同一区二区软件| 我的亚洲天堂| 久久ye,这里只有精品| 制服诱惑二区| 电影成人av| 国产精品av久久久久免费| 亚洲精品视频女| 天天躁夜夜躁狠狠久久av| 一级片免费观看大全| 大话2 男鬼变身卡| 欧美最新免费一区二区三区| 中文字幕人妻丝袜制服| 午夜福利网站1000一区二区三区| 高清黄色对白视频在线免费看| 色精品久久人妻99蜜桃| 久久人妻熟女aⅴ| 国产成人精品无人区| 天天操日日干夜夜撸| 国产熟女欧美一区二区| 你懂的网址亚洲精品在线观看| 久久免费观看电影| 美女大奶头黄色视频| 久久免费观看电影| 丝袜美腿诱惑在线| 在线亚洲精品国产二区图片欧美| av女优亚洲男人天堂| 99九九在线精品视频| 黄色毛片三级朝国网站| 中文字幕人妻丝袜制服| 亚洲精品自拍成人| 亚洲熟女精品中文字幕| 欧美亚洲 丝袜 人妻 在线| 久久精品久久久久久噜噜老黄| 日韩制服骚丝袜av| 一本色道久久久久久精品综合| 国产亚洲午夜精品一区二区久久| 国产人伦9x9x在线观看| 国产成人啪精品午夜网站| 日本爱情动作片www.在线观看| 最新的欧美精品一区二区| 国产欧美日韩综合在线一区二区| 国产精品成人在线| 国产极品粉嫩免费观看在线| 国产精品偷伦视频观看了| 桃花免费在线播放| 9色porny在线观看| 嫩草影院入口| 欧美成人午夜精品| 国产无遮挡羞羞视频在线观看| 老熟女久久久| 久久精品aⅴ一区二区三区四区| 国产一区有黄有色的免费视频| 丝瓜视频免费看黄片| 国产xxxxx性猛交| 人妻 亚洲 视频| 1024香蕉在线观看| 成人午夜精彩视频在线观看| 黄色视频在线播放观看不卡| 国产一级毛片在线| 丁香六月天网| 国产日韩欧美亚洲二区| 久久人人97超碰香蕉20202| 亚洲精品国产区一区二| 国产亚洲一区二区精品| 少妇人妻精品综合一区二区| 高清欧美精品videossex| 国产精品久久久久久精品电影小说| 国产精品99久久99久久久不卡 | 亚洲四区av| 国产成人91sexporn| 电影成人av| 成人手机av| 少妇被粗大猛烈的视频| 免费不卡黄色视频| 色综合欧美亚洲国产小说| 大码成人一级视频| 久久久久国产一级毛片高清牌| 男女边吃奶边做爰视频| 少妇人妻 视频| av片东京热男人的天堂| 精品一区二区免费观看| 日本猛色少妇xxxxx猛交久久| 少妇猛男粗大的猛烈进出视频| 女性被躁到高潮视频| 侵犯人妻中文字幕一二三四区| 精品亚洲成a人片在线观看| 日本欧美视频一区| 老汉色av国产亚洲站长工具| 老鸭窝网址在线观看| 极品人妻少妇av视频| 美女福利国产在线| 少妇人妻精品综合一区二区| 国产精品久久久av美女十八| 成人亚洲精品一区在线观看| 国产熟女午夜一区二区三区| 欧美 日韩 精品 国产| 五月开心婷婷网| 免费观看a级毛片全部| 午夜91福利影院| 在线亚洲精品国产二区图片欧美| 男女边摸边吃奶| 精品第一国产精品| 又黄又粗又硬又大视频| 国产亚洲精品第一综合不卡| 国产成人午夜福利电影在线观看| 视频在线观看一区二区三区| 一级毛片 在线播放| 丰满乱子伦码专区| 亚洲一区中文字幕在线| 一区二区日韩欧美中文字幕| 午夜福利网站1000一区二区三区| 日韩中文字幕视频在线看片| av卡一久久| 亚洲美女搞黄在线观看| 亚洲中文av在线| 一级,二级,三级黄色视频| 国产成人免费无遮挡视频| 久久狼人影院| 国产熟女午夜一区二区三区| 少妇 在线观看| 久久久久久久久久久免费av| 国产一区二区 视频在线| 精品亚洲成a人片在线观看| 丝袜人妻中文字幕| 99久久99久久久精品蜜桃| 两个人免费观看高清视频| 男人操女人黄网站| 9191精品国产免费久久| 人人妻,人人澡人人爽秒播 | 你懂的网址亚洲精品在线观看| 一区二区av电影网| 我的亚洲天堂| 亚洲成人一二三区av| av在线app专区| 女性生殖器流出的白浆| 国产黄频视频在线观看| 亚洲天堂av无毛| 国产片内射在线| 一边亲一边摸免费视频| 一级黄片播放器| 午夜免费男女啪啪视频观看| 欧美精品人与动牲交sv欧美| 嫩草影院入口| 国产在线一区二区三区精| 欧美激情高清一区二区三区 | 亚洲精品中文字幕在线视频| 国产免费现黄频在线看| 一本—道久久a久久精品蜜桃钙片| av卡一久久| 国产色婷婷99| 亚洲精品,欧美精品| 高清不卡的av网站| 另类亚洲欧美激情| av免费观看日本| av卡一久久| 啦啦啦中文免费视频观看日本| 女的被弄到高潮叫床怎么办| 久久免费观看电影| 中文字幕av电影在线播放| 国产亚洲精品第一综合不卡| 少妇被粗大猛烈的视频| 色婷婷久久久亚洲欧美| 老汉色av国产亚洲站长工具| 最近最新中文字幕免费大全7| 韩国高清视频一区二区三区| 久久毛片免费看一区二区三区| 欧美精品一区二区免费开放| 亚洲综合精品二区| 黄片小视频在线播放| 国产欧美日韩综合在线一区二区| 97人妻天天添夜夜摸| 欧美日韩综合久久久久久| 亚洲国产欧美在线一区| 欧美乱码精品一区二区三区| 午夜福利,免费看| bbb黄色大片| av不卡在线播放| 男女无遮挡免费网站观看| 精品一区二区三区四区五区乱码 | 在线观看三级黄色| 国产一区亚洲一区在线观看| 中文字幕高清在线视频| 欧美97在线视频| 一级片免费观看大全| 久久 成人 亚洲| 巨乳人妻的诱惑在线观看| 国产成人精品久久二区二区91 | 老司机影院成人| 久久鲁丝午夜福利片| 99热网站在线观看| 一区二区三区激情视频| 大话2 男鬼变身卡| 久久精品久久精品一区二区三区| 伊人久久国产一区二区| 婷婷色麻豆天堂久久| 美国免费a级毛片| 另类亚洲欧美激情| 国产成人欧美在线观看 | 精品少妇久久久久久888优播| 性色av一级| 乱人伦中国视频| 最近中文字幕2019免费版| a级毛片在线看网站| 视频在线观看一区二区三区| 在线观看免费午夜福利视频| 午夜久久久在线观看| netflix在线观看网站| 欧美亚洲日本最大视频资源| 国产精品麻豆人妻色哟哟久久| 久久av网站| 男的添女的下面高潮视频| 欧美日韩一级在线毛片| 国产精品欧美亚洲77777| 国产精品久久久久成人av| 亚洲伊人色综图| 久久女婷五月综合色啪小说| 日韩制服骚丝袜av| 亚洲av电影在线观看一区二区三区| 水蜜桃什么品种好| 丝袜喷水一区| 肉色欧美久久久久久久蜜桃| 女人精品久久久久毛片| 青春草亚洲视频在线观看| 午夜91福利影院| 欧美亚洲日本最大视频资源| 日韩欧美一区视频在线观看| 性高湖久久久久久久久免费观看| 18在线观看网站| 亚洲伊人久久精品综合| 国产成人免费无遮挡视频| 欧美在线黄色| 五月开心婷婷网| 精品第一国产精品| 一级爰片在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产老妇伦熟女老妇高清| 99精国产麻豆久久婷婷| 2018国产大陆天天弄谢| 亚洲伊人色综图| 永久免费av网站大全| 亚洲综合精品二区| 日韩熟女老妇一区二区性免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 另类精品久久| 亚洲av日韩精品久久久久久密 | 国产一区二区激情短视频 | 在线天堂中文资源库| 黄片无遮挡物在线观看| 丰满迷人的少妇在线观看| 青春草国产在线视频| 秋霞伦理黄片| 国产女主播在线喷水免费视频网站| 成人黄色视频免费在线看| 精品少妇黑人巨大在线播放| 大话2 男鬼变身卡| 精品国产一区二区久久| 久久久久人妻精品一区果冻| 少妇人妻精品综合一区二区| 搡老岳熟女国产| 黄频高清免费视频| 国产精品女同一区二区软件| 久久精品熟女亚洲av麻豆精品| 最近中文字幕2019免费版| 亚洲婷婷狠狠爱综合网| 亚洲人成77777在线视频| 日日爽夜夜爽网站| 色94色欧美一区二区| 国产免费视频播放在线视频| 老司机影院毛片| 99久久99久久久精品蜜桃| 国产成人精品在线电影| 国产免费一区二区三区四区乱码| 精品国产露脸久久av麻豆| 午夜日韩欧美国产| 国产熟女欧美一区二区| www日本在线高清视频| 亚洲第一区二区三区不卡| 一区二区三区乱码不卡18| 国产精品av久久久久免费| 人妻人人澡人人爽人人| 亚洲国产精品一区二区三区在线| 国产老妇伦熟女老妇高清| 满18在线观看网站| 欧美日韩av久久| 欧美精品人与动牲交sv欧美| 一边摸一边做爽爽视频免费| 日本av免费视频播放| 国产精品久久久久久精品古装| 日本av免费视频播放| 国产精品女同一区二区软件| 久久天堂一区二区三区四区| 2018国产大陆天天弄谢| 日日啪夜夜爽| 十分钟在线观看高清视频www| 女人被躁到高潮嗷嗷叫费观| 亚洲av中文av极速乱| 天美传媒精品一区二区| 亚洲男人天堂网一区| 免费久久久久久久精品成人欧美视频| 亚洲专区中文字幕在线 | 国产又爽黄色视频| 国产精品一区二区在线观看99| 精品一区二区三区av网在线观看 | 久久久久人妻精品一区果冻| 天天影视国产精品| 日本av免费视频播放| 免费黄色在线免费观看| 在线看a的网站| 亚洲成人手机| 国产男女超爽视频在线观看| 丁香六月欧美| 亚洲在久久综合| 久久久久久久久免费视频了| 欧美老熟妇乱子伦牲交| 久久久精品94久久精品| 国产精品久久久人人做人人爽| 大片电影免费在线观看免费| 国产日韩欧美亚洲二区| 777米奇影视久久| 国产亚洲精品第一综合不卡| 国产麻豆69| 丰满乱子伦码专区| 亚洲国产精品999| 国产亚洲午夜精品一区二区久久| 精品国产一区二区三区四区第35| 欧美少妇被猛烈插入视频| 国产xxxxx性猛交| 亚洲精品乱久久久久久| 超碰97精品在线观看| 永久免费av网站大全| 一级毛片我不卡| 黄频高清免费视频| 最近手机中文字幕大全| 久久女婷五月综合色啪小说| 亚洲伊人久久精品综合| 免费黄频网站在线观看国产| 亚洲美女视频黄频| 亚洲专区中文字幕在线 | 天堂俺去俺来也www色官网| 在线看a的网站| 9191精品国产免费久久| 久久人妻熟女aⅴ| 王馨瑶露胸无遮挡在线观看| videosex国产| 精品久久久精品久久久| 精品国产乱码久久久久久男人| 我的亚洲天堂| 亚洲av日韩精品久久久久久密 | 亚洲中文av在线| 好男人视频免费观看在线| 欧美激情极品国产一区二区三区| 久久99热这里只频精品6学生| 黄色一级大片看看| 日本av手机在线免费观看| 久久毛片免费看一区二区三区| 99热国产这里只有精品6| 超碰成人久久| 满18在线观看网站| 久久精品亚洲熟妇少妇任你| 成人毛片60女人毛片免费| 国产精品嫩草影院av在线观看| 亚洲人成77777在线视频| 中国三级夫妇交换| 在线观看免费视频网站a站| 国产精品亚洲av一区麻豆 | 男的添女的下面高潮视频| 亚洲欧美一区二区三区黑人| 五月天丁香电影| 中文字幕人妻丝袜一区二区 | 天天躁夜夜躁狠狠躁躁| 街头女战士在线观看网站| 狂野欧美激情性bbbbbb| 久久女婷五月综合色啪小说| 欧美日韩国产mv在线观看视频| 在线天堂中文资源库| 国产视频首页在线观看| 九草在线视频观看| 中文字幕高清在线视频| 久久精品国产亚洲av涩爱| 国产高清国产精品国产三级| 国产成人精品福利久久| 卡戴珊不雅视频在线播放| 老汉色av国产亚洲站长工具| 亚洲av福利一区| 国产午夜精品一二区理论片| 精品一区在线观看国产| 蜜桃在线观看..| 国产1区2区3区精品| 精品少妇黑人巨大在线播放| 高清不卡的av网站| 男女之事视频高清在线观看 | 亚洲人成网站在线观看播放| 久久精品久久久久久噜噜老黄| 国产精品久久久久成人av| 考比视频在线观看| 女人精品久久久久毛片| 女的被弄到高潮叫床怎么办| 各种免费的搞黄视频| 欧美av亚洲av综合av国产av | 久久久精品区二区三区| 美女主播在线视频| 999久久久国产精品视频| 十八禁人妻一区二区| 久久久久久久国产电影| 日本vs欧美在线观看视频| 国产乱人偷精品视频| 波多野结衣一区麻豆| 欧美日韩亚洲高清精品| 99精品久久久久人妻精品| 老汉色∧v一级毛片| 少妇人妻 视频| av在线app专区| 成人国产av品久久久| 久久人人爽人人片av| 久久久亚洲精品成人影院| 欧美日韩av久久| 日韩人妻精品一区2区三区| 婷婷色综合www| 亚洲精品自拍成人| 久久人人爽人人片av| 免费人妻精品一区二区三区视频| 秋霞在线观看毛片| 精品一区二区三区四区五区乱码 | 黄色视频在线播放观看不卡| 国产成人精品在线电影| 下体分泌物呈黄色| 国产在线免费精品| 午夜日韩欧美国产| 国产成人精品福利久久| 色播在线永久视频| 女人久久www免费人成看片| 热99国产精品久久久久久7| 国产女主播在线喷水免费视频网站| 下体分泌物呈黄色| 51午夜福利影视在线观看| www日本在线高清视频| 精品福利永久在线观看| 女人被躁到高潮嗷嗷叫费观| 性少妇av在线| av在线app专区| 免费观看a级毛片全部| 国产日韩欧美亚洲二区| 黄片小视频在线播放| 欧美乱码精品一区二区三区| av网站免费在线观看视频| 欧美成人午夜精品| 91aial.com中文字幕在线观看| 纯流量卡能插随身wifi吗| 精品国产一区二区久久| 欧美精品高潮呻吟av久久| 精品久久蜜臀av无| 亚洲欧美色中文字幕在线| 成人国语在线视频| 中文字幕人妻丝袜一区二区 | a级片在线免费高清观看视频| 国产男女内射视频| 无遮挡黄片免费观看| 七月丁香在线播放| 熟妇人妻不卡中文字幕| 国产熟女午夜一区二区三区| 久久人人97超碰香蕉20202| 最新的欧美精品一区二区| 卡戴珊不雅视频在线播放| 亚洲精品一二三| 亚洲情色 制服丝袜| 老汉色av国产亚洲站长工具| 中文字幕高清在线视频| 欧美激情 高清一区二区三区| 黄色毛片三级朝国网站| 免费女性裸体啪啪无遮挡网站| 丰满饥渴人妻一区二区三| 欧美日韩视频精品一区| 亚洲欧美激情在线| 视频区图区小说| 免费黄频网站在线观看国产| 青春草亚洲视频在线观看| 国产精品.久久久| 久久97久久精品| 青春草视频在线免费观看| 激情五月婷婷亚洲| 综合色丁香网| 精品一品国产午夜福利视频| 国产极品粉嫩免费观看在线| 国产精品久久久久久久久免| 色网站视频免费| 777米奇影视久久| 精品人妻在线不人妻| 中文精品一卡2卡3卡4更新| 热99久久久久精品小说推荐| 亚洲天堂av无毛| 日韩大片免费观看网站| 天天影视国产精品| 亚洲精品国产色婷婷电影| 国产精品久久久人人做人人爽| 亚洲情色 制服丝袜| 五月开心婷婷网| 最近手机中文字幕大全| 制服丝袜香蕉在线| 国产亚洲av高清不卡| 视频区图区小说| 国产精品香港三级国产av潘金莲 | 久久性视频一级片| 999久久久国产精品视频| 伦理电影免费视频| 国产伦人伦偷精品视频| av网站免费在线观看视频| 一级片'在线观看视频| 久久亚洲国产成人精品v| 在线观看免费高清a一片| 一级毛片黄色毛片免费观看视频| a 毛片基地| 街头女战士在线观看网站| 天天躁夜夜躁狠狠躁躁| 国产精品亚洲av一区麻豆 | 国产精品一国产av| 国产乱人偷精品视频| 一边亲一边摸免费视频| 高清视频免费观看一区二区| 丝袜美足系列| 纵有疾风起免费观看全集完整版| 精品国产超薄肉色丝袜足j| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美色中文字幕在线| 中文字幕最新亚洲高清| 女人久久www免费人成看片| 老司机亚洲免费影院| 大香蕉久久成人网| 亚洲图色成人|