服務(wù)熱線:13718308763 13718987307       E-mail:[email protected]
        中國(guó)微生物菌種

        標(biāo)準(zhǔn)物質(zhì)查詢(xún)網(wǎng)
        China microbial strains and
        reference materials query network
        會(huì)員登錄 會(huì)員注冊(cè)
        我已閱讀并同意《服務(wù)協(xié)議》
        注冊(cè)
        微生物學(xué)通報(bào)】基于5′mRNA二級(jí)結(jié)構(gòu)優(yōu)化的鹵醇脫鹵酶高效表達(dá)及應(yīng)用
        發(fā)布日期:2025-12-19 15:26:29


        微生物學(xué)通報(bào)】基于5′mRNA二級(jí)結(jié)構(gòu)優(yōu)化的鹵醇脫鹵酶高效表達(dá)及應(yīng)用


        摘 要

        目的鹵醇脫鹵酶作為一種生物催化劑,兼具閉環(huán)與開(kāi)環(huán)反應(yīng)活性,在手性環(huán)氧化合物等合成領(lǐng)域應(yīng)用廣泛,其綠色高效制備具有重要意義。本研究通過(guò)對(duì)鹵醇脫鹵酶HheC編碼序列N端進(jìn)行基于mRNA二級(jí)結(jié)構(gòu)的序列優(yōu)化,實(shí)現(xiàn)HheC的高效表達(dá),并將其應(yīng)用于手性環(huán)氧氯丙烷的合成。

        方法采用mRNA預(yù)測(cè)工具預(yù)測(cè)5′mRNA二級(jí)結(jié)構(gòu)與熱力學(xué)性質(zhì),以降低mRNA二級(jí)結(jié)構(gòu)穩(wěn)定性、提高折疊自由能(ΔG)為目標(biāo),在不改變氨基酸序列的前提下設(shè)計(jì)高表達(dá)5′mRNA序列,并表征酶的表達(dá)效率和催化性能。

        結(jié)果基于5′mRNA序列設(shè)計(jì)獲得了HheC突變體,其蛋白表達(dá)水平從16.71%提升至33.39%,對(duì)1,3-二氯-2-丙醇的相對(duì)酶活提高了3倍。

        結(jié)論基于5′mRNA二級(jí)結(jié)構(gòu)優(yōu)化可顯著提高HheC表達(dá)水平,提升目標(biāo)產(chǎn)物合成效率,為構(gòu)建酶法催化合成手性環(huán)氧氯丙烷路線奠定了基礎(chǔ)。

        關(guān)鍵詞 鹵醇脫鹵酶;mRNA二級(jí)結(jié)構(gòu);翻譯表達(dá)調(diào)控

        鹵醇脫鹵酶(halohydrin dehalogenases, HHDHs, EC 4.5.1.X)是一類(lèi)重要的生物催化劑,它能夠催化鄰鹵醇羥基對(duì)鹵素進(jìn)行分子內(nèi)親核取代反應(yīng),生成對(duì)應(yīng)的環(huán)氧化物與鹵化氫;同時(shí),也能催化相應(yīng)的逆反應(yīng),即受N3-、CN-、NO2-等親核試劑進(jìn)攻,催化環(huán)氧化物開(kāi)環(huán),具有與轉(zhuǎn)氨酶、大部分脫氫酶類(lèi)似的雙向催化功能。根據(jù)序列同源性與結(jié)構(gòu)特征的差異,鹵醇脫鹵酶可分為A-G等多個(gè)亞型[1-2]。其中,HheC對(duì)短鏈鄰鹵醇具有高活性和高選擇性,在手性環(huán)氧氯丙烷、他汀側(cè)鏈、惡唑烷酮等合成中具有重要應(yīng)用,其高效制備也備受關(guān)注[3-4]。
        在現(xiàn)代生物技術(shù)領(lǐng)域,外源基因在宿主生物中的高效表達(dá)已成為眾多應(yīng)用的基石。影響外源基因在宿主中高效表達(dá)的關(guān)鍵因素是多方面的,包括基因序列的優(yōu)化、合適啟動(dòng)子及調(diào)控元件的選擇,以及對(duì)宿主細(xì)胞翻譯機(jī)制的調(diào)控。其中,mRNA結(jié)構(gòu),尤其是5′mRNA二級(jí)結(jié)構(gòu)被視為關(guān)鍵因素之一。mRNA二級(jí)結(jié)構(gòu)能夠顯著影響翻譯起始效率,進(jìn)而影響所編碼蛋白質(zhì)的整體表達(dá)水平[5-8]。適當(dāng)?shù)膍RNA二級(jí)結(jié)構(gòu)有助于蛋白質(zhì)的正確折疊[5],改善核糖體的分配,防止核糖體碰撞與堵塞,提高翻譯效率[6,9]。反之,高穩(wěn)定性的二級(jí)結(jié)構(gòu)可能阻礙核糖體的移動(dòng),使核糖體需要更長(zhǎng)時(shí)間解開(kāi)配對(duì)的堿基,導(dǎo)致“翻譯停滯”[8];核糖體下游高穩(wěn)定性的二級(jí)結(jié)構(gòu)甚至可能引發(fā)移碼突變[10]。此外,ΔG是衡量RNA二級(jí)結(jié)構(gòu)穩(wěn)定性的關(guān)鍵熱力學(xué)參數(shù),反映了RNA分子折疊成特定結(jié)構(gòu)時(shí)的能量變化。ΔG越低(即負(fù)值越大)表明形成該結(jié)構(gòu)需要釋放更多能量,二級(jí)結(jié)構(gòu)越穩(wěn)定。
        對(duì)于翻譯起始而言,翻譯起始區(qū)(translation initiation region, TIR)是啟動(dòng)翻譯的重要功能區(qū)域,包括核糖體結(jié)合位點(diǎn)(ribosome binding site, RBS)、SD序列(Shine-Dalgarno sequence)、起始密碼子及其下游的部分序列[11]。TIR序列對(duì)翻譯起始的影響可歸因于區(qū)域內(nèi)mRNA二級(jí)結(jié)構(gòu)引起的核糖體可及性差異。TIR內(nèi)一定程度的mRNA二級(jí)結(jié)構(gòu)能夠延長(zhǎng)其半衰期,防止其被過(guò)早降解[12-13],且由于30S亞基與RBS之間具有較強(qiáng)的親和力,并不影響二者的結(jié)合[14-15]。當(dāng)mRNA穩(wěn)定性達(dá)到一定閾值時(shí),其穩(wěn)定性與翻譯效率呈負(fù)相關(guān),穩(wěn)定的二級(jí)結(jié)構(gòu)會(huì)阻礙二者結(jié)合,從而抑制翻譯起始,顯著影響目的基因表達(dá)[16-17]。此外,mRNA二級(jí)結(jié)構(gòu)在調(diào)控目的蛋白表達(dá)水平中發(fā)揮著多維度作用,涉及核糖體結(jié)合效率、翻譯起始復(fù)合體組裝、mRNA降解動(dòng)力學(xué)及其與轉(zhuǎn)錄后調(diào)控網(wǎng)絡(luò)的耦合機(jī)制。結(jié)構(gòu)緊密的5′mRNA二級(jí)結(jié)構(gòu)可遮蔽核糖體結(jié)合位點(diǎn)(如起始密碼子AUG附近或SD序列),阻礙小亞基裝載,顯著抑制翻譯起始[18];還能通過(guò)限制eIF4E、eIF4G等起始因子對(duì)mRNA帽結(jié)構(gòu)的結(jié)合,抑制翻譯起始復(fù)合體的形成[19]。由此可見(jiàn),mRNA二級(jí)結(jié)構(gòu)的折疊方式與穩(wěn)定性能夠通過(guò)多種機(jī)制影響蛋白質(zhì)的合成速率,這為聚焦于5′mRNA二級(jí)結(jié)構(gòu)設(shè)計(jì),提升鹵醇脫鹵酶在大腸桿菌(Escherichia coli)中的異源表達(dá)水平提供了重要指導(dǎo)。
        本研究以HHDH亞型之一的鹵醇脫鹵酶HheC為研究對(duì)象。通過(guò)降低編碼序列5′mRNA二級(jí)結(jié)構(gòu)穩(wěn)定性在不改變氨基酸序列的前提下提高HheC在大腸桿菌中的異源表達(dá)水平,并表征酶的表達(dá)效率及其對(duì)1,3-二氯-2-丙醇(1,3-dichloro-2-propanol, 1,3-DCP)的催化性能,以期為手性環(huán)氧氯丙烷的酶法合成奠定基礎(chǔ)。

        1 材料與方法

        1.1 材料

        1.1.1 菌株
        本研究使用的表達(dá)宿主為E. coli BL21(DE3),表達(dá)載體為pET-28b(+),二者均由本實(shí)驗(yàn)室保藏;研究對(duì)象鹵醇脫鹵酶HheCL143N/A159T/P175S/W249P為實(shí)驗(yàn)室前期構(gòu)建,標(biāo)記為HheCM4。
        1.1.2 培養(yǎng)基
        LB培養(yǎng)基(g/L):蛋白胨10.0,酵母粉5.0,NaCl 10.0,pH 7.4。
        LB固體培養(yǎng)基:在LB液體培養(yǎng)基的基礎(chǔ)上加入20 g/L的瓊脂粉配制而成。

        1.2 mRNA二級(jí)結(jié)構(gòu)與熱力學(xué)性質(zhì)預(yù)測(cè)

        1.2.1 基于ViennaRNA的mRNA二級(jí)結(jié)構(gòu)與熱力學(xué)性質(zhì)分析
        為提高mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果的可信度,準(zhǔn)確反映熱力學(xué)性質(zhì)變化趨勢(shì),本研究采用3種不同算法的預(yù)測(cè)工具對(duì)mRNA二級(jí)結(jié)構(gòu)與熱力學(xué)性質(zhì)進(jìn)行分析。利用ViennaRNA在線工具(http://rna.tbi.univie.ac.at/)[20-21]對(duì)5′mRNA (基于起始密碼子的80 nt序列)的二級(jí)結(jié)構(gòu)、折疊自由能變化(ΔG)及核苷酸熱力學(xué)性質(zhì)進(jìn)行預(yù)測(cè)。其mRNA折疊算法包括最小自由能(minimum free energy, MFE)模型與配分函數(shù)(partition function)模型。MFE模型通過(guò)搜索自由能最低的構(gòu)象來(lái)預(yù)測(cè)序列最穩(wěn)定的折疊結(jié)構(gòu);配分函數(shù)模型基于統(tǒng)計(jì)力學(xué)框架計(jì)算所有可能構(gòu)象的熱力學(xué)加權(quán),從而獲得堿基配對(duì)概率與結(jié)構(gòu)多樣性信息,并設(shè)置避免孤對(duì)堿基對(duì)(avoid isolated base pairs)以提高預(yù)測(cè)結(jié)果的合理性。選擇Turner模型(2004)[22]的熱力學(xué)參數(shù),設(shè)定折疊溫度為37 ℃,在1 mol/L鹽濃度下進(jìn)行能量計(jì)算。
        1.2.2 基于Mfold的mRNA二級(jí)結(jié)構(gòu)與熱力學(xué)性質(zhì)分析
        源于UNAFold在線工具(https://www.unafold.org/mfold/applications/rna-folding-form.php)的Mfold模塊[23]同樣是基于最小自由能原理開(kāi)發(fā)的核酸二級(jí)結(jié)構(gòu)預(yù)測(cè)算法。其核心采用Zuker等提出的動(dòng)態(tài)規(guī)劃(dynamic programming)算法[24]與Turner模型(1999)[25]的熱力學(xué)參數(shù)系統(tǒng)性地搜索所有可能的堿基配對(duì)組合,從而精確識(shí)別最穩(wěn)定構(gòu)象及多個(gè)亞穩(wěn)態(tài)結(jié)構(gòu)。該模塊提供靈活的結(jié)構(gòu)約束功能,支持強(qiáng)制配對(duì)、禁配設(shè)定和多鏈雜交模擬,在寡核苷酸設(shè)計(jì)、引物篩選、RNA干擾等領(lǐng)域均有廣泛應(yīng)用。設(shè)定折疊溫度為37 ℃,在1 mol/L鹽濃度下進(jìn)行能量計(jì)算。
        1.2.3 基于UFold的mRNA二級(jí)結(jié)構(gòu)分析
        UFold[26]基于深度學(xué)習(xí)框架,通過(guò)端到端的圖像化神經(jīng)網(wǎng)絡(luò)模型實(shí)現(xiàn)高精度與高效率的mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)。其主體采用U-Net模型[27]全卷積神經(jīng)網(wǎng)絡(luò)架構(gòu),能夠在保持位置信息的同時(shí)逐層提取局部與全局特征。通過(guò)對(duì)L×L的配對(duì)矩陣進(jìn)行對(duì)稱(chēng)性處理,篩除不符合RNA二級(jí)結(jié)構(gòu)物理約束的配對(duì),并結(jié)合線性規(guī)劃進(jìn)行結(jié)構(gòu)優(yōu)化,最終獲得符合熱力學(xué)與幾何限制的穩(wěn)定二級(jí)結(jié)構(gòu)模型。該方法相比傳統(tǒng)自由能最小化方法具有更高的預(yù)測(cè)準(zhǔn)確性與處理效率,在處理假結(jié)預(yù)測(cè)與長(zhǎng)程配對(duì)識(shí)別任務(wù)中表現(xiàn)更優(yōu)秀。

        1.3 重組菌株構(gòu)建

        以HheCM4為模板,通過(guò)全質(zhì)粒PCR對(duì)目標(biāo)序列進(jìn)行突變。質(zhì)粒DNA由質(zhì)粒小提試劑盒(購(gòu)自北京擎科生物科技股份有限公司)提取獲得。引物序列見(jiàn)表1。PCR反應(yīng)體系(25 μL):2×Phanta Max Buffer (p505) 12.5 μL,dNTP Mix (10 mmol/L) 0.5 μL,Phanta Max Super-Fidelity DNA Polymerase 0.25 μL,均購(gòu)自南京諾唯贊生物科技股份有限公司,上、下游引物(10 μmol/L)各1 μL,DNA模板0.1 μL,ddH2O補(bǔ)足至25 μL。PCR反應(yīng)條件:95 ℃預(yù)變性5 min;95 ℃變性30 s,60 ℃退火30 s,72 ℃延伸3.5 min,30個(gè)循環(huán);72 ℃終延伸10 min。經(jīng)1.5%瓊脂糖凝膠電泳驗(yàn)證PCR產(chǎn)物后,用Dpn I于37 ℃消化模板,通過(guò)熱激法將產(chǎn)物轉(zhuǎn)化至E. coli BL21(DE3),涂布于含50 μg/mL卡那霉素的LB固體培養(yǎng)基,37 ℃培養(yǎng)過(guò)夜,挑取單菌落送北京擎科生物科技股份有限公司測(cè)序驗(yàn)證。

        表1   5′mRNA突變引物序列Table 1   Mutant primers of 5′mRNA sequences

        a:下游通用引物。

        a: Downstream universal primer.


        1.4 目的蛋白誘導(dǎo)表達(dá)與分析

        1.4.1 誘導(dǎo)表達(dá)
        挑取重組菌落至含50 μg/mL卡那霉素的LB液體培養(yǎng)基試管中,于37 ℃、200 r/min培養(yǎng)過(guò)夜。按2%體積分?jǐn)?shù)轉(zhuǎn)接至含50 μg/mL卡那霉素的LB液體培養(yǎng)基中,于37 ℃、180 r/min培養(yǎng)至OD600約為0.6-0.8,加入IPTG (終濃度0.1 mmol/L)進(jìn)行目標(biāo)蛋白的誘導(dǎo)表達(dá),28 ℃、180 r/min誘導(dǎo)12 h后,于8 000 r/min離心10 min,去除上清液后得到濕菌體,用于后續(xù)的分析與催化反應(yīng)。
        1.4.2 蛋白表達(dá)分析
        將50 g/L重組菌重懸于PB緩沖液(20 mmol/L, pH 7.4),在冰水浴下以180 W超聲破碎10 min (超聲時(shí)間2 s,間隔時(shí)間4 s)至澄清,4 ℃、12 000 r/min離心20 min后,取上清液用PB緩沖液稀釋3倍,并與5×loading buffer按4:1體積比混合,在沸水浴中放置5 min,取4 μL進(jìn)行SDS-PAGE分析。經(jīng)染色脫色儀10 min工作程序處理后,在凝膠成像儀下觀察電泳結(jié)果。采用圖像處理軟件ImageJ 1.54f對(duì)SDS-PAGE圖像進(jìn)行Rolling ball算法[28]背景校正后取反色,進(jìn)行目的蛋白表達(dá)量的半定量灰度分析。
        表達(dá)水平與酶活數(shù)據(jù)包含3個(gè)生物學(xué)平行重復(fù)組,以平均值±標(biāo)準(zhǔn)差(mean±SD)表示。組間顯著差異分析采用單因素方差分析(one-way ANOVA),使用Waller-Duncan多重比較法進(jìn)行差異檢測(cè)。顯著性水平設(shè)定為P<0.05,所有統(tǒng)計(jì)分析均使用SPSS 27.0軟件完成。

        1.5 鹵醇脫鹵酶催化1,3-DCP合成手性環(huán)氧氯丙烷

        1.5.1 反應(yīng)體系與反應(yīng)條件
        在20 mL Gly/NaOH緩沖液(600 mmol/L, pH 9.8)與乙酸乙酯體積比為1:1的反應(yīng)體系中加入含鹵醇脫鹵酶突變體的重組菌0.2 g (終濃度為10 g/L),400 mmol/L 1,3-DCP,于30 ℃孵育10 min后,180 r/min水浴振蕩反應(yīng)5 min后取樣,12 000 r/min離心1 min,取上層有機(jī)相經(jīng)無(wú)水硫酸鈉處理后進(jìn)行氣相色譜檢測(cè)。
        1.5.2 氣相檢測(cè)方法
        底物1,3-DCP與產(chǎn)物環(huán)氧氯丙烷(epichlorohydrin, ECH)使用GC (Agilent公司)進(jìn)行檢測(cè),具體檢測(cè)條件如下:色譜柱型號(hào)為BGB-175手性毛細(xì)管柱(0.25 mm×30 m×0.25 μm);柱箱條件為初始柱溫100 ℃,保持4 min后以2 ℃/min升溫至110 ℃,而后以30 ℃/min升高至180 ℃,保持3 min。進(jìn)樣口溫度設(shè)定為250 ℃,F(xiàn)ID檢測(cè)口溫度設(shè)定為250 ℃。載氣為高純氮,流量設(shè)定為1.5 mL/min,分流比為80:1。

        2 結(jié)果與分析

        2.1 5′mRNA序列設(shè)計(jì)

        基于翻譯起始區(qū)TIR中5′mRNA二級(jí)結(jié)構(gòu)穩(wěn)定性與目標(biāo)蛋白表達(dá)水平呈負(fù)相關(guān)的研究結(jié)論[5-8],本研究通過(guò)降低mRNA二級(jí)結(jié)構(gòu)穩(wěn)定性以提高翻譯起始效率,對(duì)HheC編碼序列的5′mRNA序列進(jìn)行優(yōu)化,從而增加HheC表達(dá)。在不改變氨基酸序列的前提下,分別以提高和降低ΔG為目標(biāo)對(duì)序列進(jìn)行設(shè)計(jì),并設(shè)計(jì)相應(yīng)引物對(duì)5′mRNA核苷酸序列進(jìn)行突變。其中,2號(hào)位丙氨酸插入突變ins2A為實(shí)驗(yàn)中意外所得,上標(biāo)為相應(yīng)密碼子編碼的丙氨酸,其余序列均編碼與HheCM4相同的氨基酸。
        以起始密碼子上游不同位點(diǎn)作為最小折疊自由能預(yù)測(cè)的起始位點(diǎn)均能得到相近的自由能預(yù)測(cè)結(jié)果(R2=0.938 5)[29]。因此,本研究以起始密碼子(AUG) 80 nt的序列(-20-60)作為研究對(duì)象,預(yù)測(cè)mRNA二級(jí)結(jié)構(gòu)并計(jì)算自由能。如表2所示,不同密碼子對(duì)應(yīng)的2號(hào)位丙氨酸插入突變體與出發(fā)酶HheCM4相比,mRNA二級(jí)結(jié)構(gòu)確實(shí)存在折疊自由能差異。由此,設(shè)計(jì)預(yù)期mRNA二級(jí)結(jié)構(gòu)穩(wěn)定性降低的高表達(dá)突變體HheCM4/H1 (ΔG=-13.20 kcal/mol)與HheCM4/H2 (ΔG=-12.00 kcal/mol),同時(shí)設(shè)計(jì)預(yù)期高穩(wěn)定性低表達(dá)的突變體HheCM4/L1 (ΔG=-21.70 kcal/mol)與HheCM4/L2 (ΔG=-24.10 kcal/mol)用于反向驗(yàn)證。在此基礎(chǔ)上進(jìn)一步表征了對(duì)應(yīng)突變體在宿主中的表達(dá)水平。

        表2   5′mRNA序列設(shè)計(jì)Table 2   Design of synonymous mutation sequences of 5′mRNA

        a:下劃線代表起始密碼;b:源于ViennaRNA;c:源于Mfold。

        a: The underline represents the start codon; b: Calculated by ViennaRNA; c: Calculated by Mfold.


        2.2 目標(biāo)蛋白表達(dá)水平與催化活性研究

        SDS-PAGE圖像與灰度掃描半定量分析結(jié)果如圖1所示,出發(fā)株HheCM4中HheC (28.7 kDa)占宿主細(xì)胞中可溶性蛋白的比例為16.71%。4種插入突變體中,ins2AGCA表達(dá)水平變化不顯著(16.61%),而ins2AGCU、ins2AGCC、ins2AGCG表達(dá)水平提升幅度相近,HheC占比分別為25.42%、27.36%與26.71%。序列設(shè)計(jì)結(jié)果中具有更高mRNA穩(wěn)定性的HheCM4/L2 (12.38%)表達(dá)水平低于HheCM4/L1 (14.98%),且較出發(fā)株均有所下降;而具有更低穩(wěn)定性的HheCM4/H2 (33.39%)表達(dá)水平高于HheCM4/H1 (31.78%),較出發(fā)株提升顯著,與設(shè)計(jì)預(yù)期一致。
        圖1   HheCM4及其突變體的表達(dá)水平Figure 1   Expression levels of HheCM4 and its mutants. A: SDS-PAGE analysis chart of HheCM4 and its mutants (Lane 1: Protein marker; Lane 2: HheCM4; Lane 3: ins2AGCA; Lane 4: ins2AGCU; Lane 5: ins2AGCC; Lane 6: ins2AGCG; Lane 7: HheCM4/L1; Lane 8: HheCM4/L2; Lane 9: HheCM4/H1; Lane 10: HheCM4/H2); B: Comparison of the soluble expression of HheCM4 and its mutants [Different lowercase letters indicate that there are significant differences in the expression level among different strains (P<0.05)].
        HheCM4及其突變體的全細(xì)胞催化活性測(cè)定結(jié)果如圖2所示,與表達(dá)水平測(cè)定結(jié)果對(duì)應(yīng),4種插入突變體中ins2AGCA具有最低的催化活性,較出發(fā)株HheCM4不存在顯著變化,而其余3種插入突變體催化活性均有顯著提升。序列設(shè)計(jì)結(jié)果中具有高mRNA穩(wěn)定性的HheCM4/L1與HheCM4/L2相對(duì)酶活均不及出發(fā)株,且具有更低ΔG的HheCM4/L2 (57.95%)催化活性低于HheCM4/L1 (72.18%);而mRNA穩(wěn)定性更低的HheCM4/H1 (278.85%)與HheCM4/H2 (300.78%)催化活性均存在顯著提升,且具有最高ΔG的HheCM4/H2表現(xiàn)出最高的催化活性,與表達(dá)水平以及mRNA結(jié)構(gòu)穩(wěn)定性預(yù)測(cè)結(jié)果一致。
        圖2   HheCM4及其突變體的相對(duì)酶活比較Figure 2   Comparison of the relative enzyme activities of HheCM4 and its mutants. Different lowercase letters indicate that there are significant differences in the relative activity among different strains (P<0.05).
        基于HheC表達(dá)水平與ΔG相關(guān)性繪制的散點(diǎn)圖如圖3所示。采用2種不同算法工具得到的ΔG結(jié)果均反映mRNA熱力學(xué)穩(wěn)定性與HheC表達(dá)量呈負(fù)相關(guān)性,即隨著自由能ΔG的升高(mRNA穩(wěn)定性降低) HheC表達(dá)量顯著增加,且2種算法的預(yù)測(cè)結(jié)果在分布趨勢(shì)與數(shù)據(jù)密度上基本一致。
        圖3   mRNA二級(jí)結(jié)構(gòu)折疊自由能(ΔG)與HheC表達(dá)水平相關(guān)性分析Figure 3   Correlation analysis of the folding free energy (ΔG) of mRNA secondary structure and the expression level of HheC.

        2.3 mRNA二級(jí)結(jié)構(gòu)分析

        MFE結(jié)構(gòu)是依據(jù)自由能最小化原則得到的、mRNA二級(jí)結(jié)構(gòu)最穩(wěn)定的構(gòu)象,代表了熱力學(xué)上最可能形成的結(jié)構(gòu)。質(zhì)心結(jié)構(gòu)(centroid structure)是通過(guò)概率模型計(jì)算所得的mRNA二級(jí)結(jié)構(gòu),表示熱力學(xué)構(gòu)象集合中最接近平均構(gòu)象的一種代表性結(jié)構(gòu)[30]。
        源于ViennaRNA的5′mRNA MFE二級(jí)結(jié)構(gòu)與質(zhì)心二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果如圖4所示。基于熵值對(duì)結(jié)構(gòu)進(jìn)行著色,反映了mRNA在該位置存在不同構(gòu)象的可能性。熵值趨于0 (暖色調(diào))表明該位點(diǎn)的結(jié)構(gòu)較不穩(wěn)定,可形成多種不同構(gòu)象,如環(huán)或單鏈結(jié)構(gòu);熵值趨于冷色調(diào)則表明該位點(diǎn)的結(jié)構(gòu)較穩(wěn)定,通常位于莖(stem)結(jié)構(gòu)中,更易于配對(duì)并形成穩(wěn)定的雙鏈結(jié)構(gòu)。HheCM4二級(jí)結(jié)構(gòu)的大部分區(qū)域形成了一個(gè)緊湊的發(fā)夾結(jié)構(gòu),且MFE結(jié)構(gòu)與質(zhì)心結(jié)構(gòu)具有較高的相似性(圖4A),表明該結(jié)構(gòu)具有一定程度的穩(wěn)定性。HheCM4/H1 (圖4H)與HheCM4/H2 (圖4I)的MFE結(jié)構(gòu)與質(zhì)心結(jié)構(gòu)相似度較低,存在更多不穩(wěn)定的非配對(duì)結(jié)構(gòu),因而易于解旋,有助于核糖體快速結(jié)合,通過(guò)提高翻譯效率增加HheC表達(dá)量。相比之下,HheCM4/L1 (圖4F)與HheCM4/L2 (圖4G)的MFE結(jié)構(gòu)與質(zhì)心結(jié)構(gòu)均呈高度相似的發(fā)夾結(jié)構(gòu),表明兩者具有高度穩(wěn)定的mRNA二級(jí)結(jié)構(gòu),從而通過(guò)阻礙翻譯起始降低HheC表達(dá)量。
        圖4   ViennaRNA mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果Figure 4   Prediction of mRNA structure by ViennaRNA. A: HheCM4; B: ins2AGCA; C: ins2AGCU; D: ins2AGCC; E: ins2AGCG; F: HheCM4/L1; G: HheCM4/L2; H: HheCM4/H1; I: HheCM4/H2.
        源于Mfold的mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果如圖5所示。由于Mfold輸出的mRNA結(jié)構(gòu)僅包括MFE結(jié)構(gòu)或熱力學(xué)能量接近MFE的次優(yōu)結(jié)構(gòu),因此大多呈現(xiàn)為相似度較高的連續(xù)莖環(huán)結(jié)構(gòu)。HheCM4 (圖5A)存在多個(gè)莖環(huán)結(jié)構(gòu),雙鏈配對(duì)區(qū)域較為集中;ins2AGCA (圖5B)結(jié)構(gòu)與其高度相似,配對(duì)密度略低;而ins2AGCU (圖5C)、ins2AGCC (圖5D)、ins2AGCG (圖5E)三者結(jié)構(gòu)高度相似,發(fā)夾短,環(huán)結(jié)構(gòu)分布更多,配對(duì)密度更低。HheCM4/L1 (圖5F)與HheCM4/L2 (圖5G)存在連續(xù)且緊密的長(zhǎng)莖結(jié)構(gòu),雙鏈配對(duì)密度顯著提升。HheCM4/H1 (圖5H)雖與HheCM4同樣具有較高的相似度,但配對(duì)密度進(jìn)一步下降,HheCM4/H2 (圖5I)則僅在少量堿基間形成較短的配對(duì)。

        圖5   Mfold mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果Figure 5   Prediction of mRNA structure by Mfold. A: HheCM4; B: ins2AGCA; C: ins2AGCU; D: ins2AGCC; E: ins2AGCG; F: HheCM4/L1; G: HheCM4/L2; H: HheCM4/H1; I: HheCM4/H2.
        源于UFold的mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果如圖6所示。HheCM4 (圖6A)主鏈存在連續(xù)的莖區(qū)雙鏈配對(duì),而插入突變體(圖6B-6E)莖區(qū)長(zhǎng)度降低,主鏈存在更多環(huán)結(jié)構(gòu)。HheCM4/L1 (圖6F)與HheCM4/L2 (圖6G)存在更長(zhǎng)的莖區(qū)和更致密的堿基配對(duì),且在HheCM4/L2中尤為顯著。相比之下,HheCM4/H2 (圖6I)配對(duì)結(jié)構(gòu)少,莖區(qū)結(jié)構(gòu)較短,發(fā)夾未完全閉合,同時(shí)更復(fù)雜的長(zhǎng)程相互作用表明其存在更大的三維構(gòu)象自由度。
        圖6   UFold mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果Figure 6   Prediction of mRNA structure by UFold. A: HheCM4; B: ins2AGCA; C: ins2AGCU; D: ins2AGCC; E: ins2AGCG; F: HheCM4/L1; G: HheCM4/L2; H: HheCM4/H1; I: HheCM4/H2.

        2.4 熱力學(xué)性質(zhì)分析

        山形圖(mountain plot)能夠基于mRNA二級(jí)結(jié)構(gòu)的MFE結(jié)構(gòu)自由能、質(zhì)心結(jié)構(gòu)自由能以及配對(duì)(partition function, PF)結(jié)構(gòu)展現(xiàn)序列特定區(qū)域的自由能變化。如圖7所示,HheCM4/L1 (圖7F)與HheCM4/L2 (圖7G)的三條曲線高度相似,且HheCM4/L2的MFE與Centroid曲線完全重合,這進(jìn)一步證明了其mRNA二級(jí)結(jié)構(gòu)具有高度穩(wěn)定性。HheCM4/H1 (圖7H)與HheCM4/H2 (圖7I)相較于HheCM4 (圖7A)自由能更低,并且HheCM4/H2中存在一個(gè)顯著的低能量區(qū)域,這或許更有利于mRNA二級(jí)結(jié)構(gòu)的部分解旋。

        圖7   mRNA二級(jí)結(jié)構(gòu)山形圖Figure 7   Mountain plot of mRNA secondary structure. A: HheCM4; B: ins2AGCA; C: ins2AGCU; D: ins2AGCC; E: ins2AGCG; F: HheCM4/L1; G: HheCM4/L2; H: HheCM4/H1; I: HheCM4/H2.
        熵(entropy)是量化二級(jí)結(jié)構(gòu)穩(wěn)定性變化的一個(gè)重要參數(shù),可用于衡量mRNA不同區(qū)域的結(jié)構(gòu)可變程度。低熵表示該位置的結(jié)構(gòu)穩(wěn)定,構(gòu)象選擇性低,在大多數(shù)情況下均能折疊成特定結(jié)構(gòu);而高熵值則相反,表明二級(jí)結(jié)構(gòu)高度動(dòng)態(tài),構(gòu)象變化較多。通過(guò)對(duì)各突變體的mRNA序列進(jìn)行熵值分析能夠進(jìn)一步揭示它們之間二級(jí)結(jié)構(gòu)的可變性和動(dòng)態(tài)特性。HheCM4、HheCM4/L1、HheCM4/L2中存在更多的低熵區(qū)域,表明其大部分區(qū)域構(gòu)象變化受限;而HheCM4/H1、HheCM4/H2中熵普遍較高,表明其二級(jí)結(jié)構(gòu)在多個(gè)區(qū)域存在較大的可變空間。

        2.5 高表達(dá)HheC催化1,3-DCP合成ECH反應(yīng)進(jìn)程

        本研究考察了HheCM4及其突變體(5 g/L濕細(xì)胞)催化100 mmol/L底物合成ECH的反應(yīng)進(jìn)程,結(jié)果如圖8所示。受限于鹵醇脫鹵酶的可逆反應(yīng),在該體系下的最大產(chǎn)率僅能達(dá)到約88%。HheCM4/H1與HheCM4/H2反應(yīng)5 min時(shí)產(chǎn)率即可分別達(dá)到14.31%與17.17%,且僅需約50 min即可達(dá)到最大產(chǎn)率;而出發(fā)株HheCM4反應(yīng)5 min時(shí)產(chǎn)率僅為7.54%,達(dá)到同一水平產(chǎn)率需將反應(yīng)時(shí)間延長(zhǎng)至約120 min。因此,在實(shí)際生產(chǎn)中本研究的高表達(dá)突變體能夠有效加速反應(yīng)進(jìn)程,在現(xiàn)有體系下完成底物的高效轉(zhuǎn)化,縮短生產(chǎn)周期,為鹵醇脫鹵酶高效催化1,3-DCP合成ECH奠定重要基礎(chǔ)。
        圖8   HheCM4及其突變體催化1,3-DCP合成ECH的反應(yīng)進(jìn)程Figure 8   Synthesis of ECH from 1,3-DCP catalyzed by HheCM4 and its mutants.

        3 討論與結(jié)論

        鹵醇脫鹵酶作為一種能夠催化鹵代醇與環(huán)氧化物相互轉(zhuǎn)化的生物催化劑,在環(huán)境修復(fù)與監(jiān)測(cè)[31-32]、合成藥物手性中間體[3-4]等領(lǐng)域均具有重要的應(yīng)用價(jià)值。另一方面,原核生物的基因表達(dá)受mRNA結(jié)構(gòu)的影響。在蛋白質(zhì)合成過(guò)程中,翻譯起始是其中限速且受調(diào)控程度最高的階段[11],易受到mRNA二級(jí)結(jié)構(gòu)穩(wěn)定性的影響。穩(wěn)定的mRNA二級(jí)結(jié)構(gòu)會(huì)阻礙核糖體30S亞基與RBS的結(jié)合,從而抑制翻譯起始,顯著影響目的基因的表達(dá)[9,16-17],而較低的mRNA二級(jí)結(jié)構(gòu)穩(wěn)定性有利于兩者結(jié)合,進(jìn)而高效啟動(dòng)翻譯[33-35]。因此,通過(guò)序列優(yōu)化降低翻譯起始區(qū)mRNA二級(jí)結(jié)構(gòu)的穩(wěn)定性能夠有效提高翻譯起始效率,實(shí)現(xiàn)外源基因在宿主細(xì)胞中的高效表達(dá)。
        本研究基于鹵醇脫鹵酶HheC編碼序列5′端mRNA二級(jí)結(jié)構(gòu)與折疊自由能,借助mRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)工具進(jìn)行計(jì)算,在維持氨基酸序列不變的前提下實(shí)現(xiàn)了核苷酸序列的優(yōu)化。最終,使HheC在E. coli BL21(DE3)中的表達(dá)水平由占宿主細(xì)胞可溶性蛋白的16.71%提升至33.39%,相對(duì)酶活提升至原來(lái)的3倍。此外,通過(guò)負(fù)向優(yōu)化得到的高穩(wěn)定性5′mRNA序列表現(xiàn)出低表達(dá)水平,進(jìn)一步驗(yàn)證了HheC 5′mRNA穩(wěn)定性與其表達(dá)水平的相關(guān)性。將所得成果應(yīng)用于催化1,3-DCP合成ECH,在本研究體積比為1:1的Gly/NaOH (600 mmol/L, pH 9.8)與乙酸異丁酯的雙相體系下,高表達(dá)HheC通過(guò)濕細(xì)胞催化100 mmol/L 1,3-DCP合成ECH達(dá)到產(chǎn)率上限的反應(yīng)時(shí)間由120 min縮短至50 min。HheC的高效表達(dá)不僅能夠提高目標(biāo)酶的產(chǎn)量,使微量活性差異更易被檢測(cè)與量化,從而提升蛋白質(zhì)工程中突變篩選的靈敏度和效率;而且對(duì)于生產(chǎn)應(yīng)用而言,更易于實(shí)現(xiàn)反應(yīng)體系的高效催化,顯著提高經(jīng)濟(jì)效益,降低生產(chǎn)成本,增強(qiáng)酶催化工藝路線的競(jìng)爭(zhēng)力。

        推薦閱讀

              久久久久亚洲精品乱码按摩 | 国产精品久久久久久久久夜色 | 97青草超碰久久国内精品91 | 男人和女人在做性视频 | 福利在线视频导航 | 色噜噜狠狠狠综合曰曰曰 | 国产爱福利 | 国产片婬乱一级毛片潮湿 | 亚洲欧美日本在线观看 | 操比视频免费 | 在线观看亚洲视频 | av黄色小说网 | 亚洲乱伦图片小说 | 欧洲成人在线视频 | AA片在线观看视频在线播放 | 男人天堂网站在线 | 免费爱爱网站 | 午夜黄色电影 | 天天想夜夜操 | 大香蕉尹人在线观看 | 91福利免费 | 国产精品 自拍偷拍 | 小早川怜子色情影片-ThePorn | 日逼短视频 | 久久丫精品 | 操嫩逼电影| 欧美一区二区三区二区 | 97超碰超碰 | 91超碰在线 | 婷婷国产天堂久久综合 | 国产精品精品国产婷婷这里Av | 色一情一乱一伦 | 91五个黑人玩一个女5p | 日本亚洲色大成网站www久久 | 大鸡吧插逼视频 | 久久婷婷五月综合色吧 | 国产美女精品人人做人人爽 | 97伊人超碰 | 欧美色综合网 | 久久久久久婷 | 中文字幕影片免费在线观看 | 无码插插 | 人人爽人人爽人人 | 亚洲国产成人精品无码电影 | 中文字幕久久青青 | 特黄特色大片免费播放叫疼 | 午夜久久久久久久久久影院 | 亚洲a色 我操大逼逼 | 丁香五月综合网 | 99青青草 | 少妇被又大又粗又爽毛片久久黑人 | 91超碰大香蕉 | www.污黄 | www.撸一撸 | 制服丝袜强奸乱伦 | 欧美丰满熟妇BBBBBB小短人 | 亚洲成人影视在线播放 | 99久热在线精品视频播放 | 蜜桃媒AV成人片免费看 | 孕交xxx | 在线看成人 | 国内毛片毛片毛片 | 人人摸人人插屄屄 | 99久久久久无码国产精品免费 | 奥门一级片 | 阿v视频在线观看 | 青青青视频在线 | 囯产精品久久久久久久久在饯观看 | 黄片无遮挡 | 97午夜影院 | 成人高清无码视频在线观看 | 留守少妇无力反抗呻吟 | 精品露脸国偷精品产拍 | 欧在线一二三四区 | 非洲人做爰全过程视频 | www.亚洲天堂 | 成人一二三视频 | 男人日女人下面视频 | 欧美黑人肏屄 | 国产美女被爽到高潮免费A片软件 | 狍与女人做爰在线视频 | 被戴sm刑具尽调教室虐男视频 | 国产一级a毛一级a爰片 | 寡妇遗像前被肉猛爆高潮 | 无码人妻精品一区二区蜜桃在 | 麻豆网视频免费观看 | 色哺乳xxxxhd牛奶a | 日韩aV在我线一区二区三反 | 久久女人| 亚洲情色在线 | 国产女主播AV在线 | 成人免费网站www污污污在线看 | 一级黄色录像带 | 自慰 免费观看 | 一级毛片AAAAAA免费看小说 | 日穴网站 | 最新无码中文字幕 | 日本裸体大乳mm77图片 | 91视频18| 四虎高清无码 | 爱爱一级视频 | 日本东京热无码视频免费 | 野外被三个男人躁一夜漫画 | 禁片天堂 | 影音先锋男人站 | 欧美性xxxxx极品娇小 | 色色色9999 | 床戏做爰无遮挡摸亲胸视频 | 全黄一级裸片视频 | 国产欧美7777 | 午夜影院亚洲 | 特黄片| 性Av播放| 无码性片 | 好屌网 | chinese性老太granny | 亚洲日韩在线免费观看 | 日韩精品一区二区成人 | 亚洲色图欧美 | 一道本视频在线 | 国产做爱xxxⅹ性视频播放量 | 开心激情网五月天 | 少妇bbb撒尿做饭爱bbw | 青青草视频www | 国产一级a毛一级a毛免费视频 | 国严人人操人人性生活 | 美女操B视频软件网战 | 五月丁香啪啪啪网 | 免费A级毛片无码免费120软件 | AV电影免费在线观看 | xxxx性开放按摩hd | 97在线无码精品秘 入口九色 | sarahtwain精品艳妇 一部一直喷奶水的av | 青青青国产精品一区二区 | 国产精品免费电影 | 亚洲综合在线一区 | 免费观看成人做爰A片免费看网站 | 欧美日韩小电影 | 国精产品一区一区三区mba· | 美女销魂免费一区二区 | 黄色电影网站在线观看 | 久久九九久精品国产免费直播 | 一级特黄AAA | 女上位av片在线观看免费 | 嗯啊用力好舒服 | 久热这里只有 | 秋霞7777鲁丝伊人久久影院 | 国产欧美日韩精品在线观看 | 91欧美日韩 | 久久久久久逼 | 少妇掰穴 | 亚洲免费小电影 | 欧美激情第八页 | 欧美午夜理伦三级在线观看 | 女fiee性zozo交体内谢 | 国产又粗又猛又大又爽的视频 | 91aaaaaa| 强烈淫药h多p | 国精产品一区一区三区mba· | 69色在线 | 操逼淫水视频 | 欧美乱伦图片 | 全黄裸体做爰视频 | brazzersxxxxhdhd1080 | 日本苍井空高潮无码a | 99re这里只有精品视频在线观看 | 成人免费黄色大片 | 亚洲无码中文字幕在线 | 扒开双腿疯狂进出爽爽爽电影 | 一级黄色高清视频 | 最新操逼网站 | 99视频免费观看 | 大查蕉在线视频色小妞 | 日本三级韩国三级欧美三级 | 91精品国产综合久久香蕉第1集 | 激情国产精品 | xxxx18hd日本hd教师 | 大鸡巴操小穴视频 | 国产成人免费1000部网站 | 国产一级特黄录像片 | 日日夜夜干 | 亚洲精品国产欧美在线观看 | 国产视频一区免费在线播放 | 久久艹天天艹 | 男女日皮免费视频 | 黑人玩弄人妻一区二区三区视频 | 9l视频自拍蝌蚪9l视频 | 麻豆3级片 | 日本一区二区三区四区五区六区 | 曰韩一级A片 | 日本一级婬片A片免费播放一 | 爆操黑丝美女 | 大鸡巴网免费视频在线 | 丰满护士巨好爽好大乳动漫 | 久久久久久久久免费观看完整版 | 日韩一区无码 | 特级西西444.WWW人体聚色 | 亚洲日、韩aⅴ | 人妖精品videosex性欧美 | 色涩视频 | 欧美精品成人网站 | 韩国三级久久电影 | 亚洲欧美校园春色 | 欧美性猛交xxxx乱大交喷浆 | 成人伊人综合 | 成人 免费 在线视频 黄 | 岳啊轻点我高潮了 | 国产综合色视频 | 特黄60分钟免费视频大全在线观看 | 国产对白农村老女人性视频对话 | 美女肏逼网 | 欧美人体大胆无码AA视频 | 极品尤物白丝拘束h | 天天日天天爽 | 三上悠亚被狂c躁到高潮失禁 | 国产女主播户外裸体勾搭 | 男生操女生下面的视频 | 黄片免费网址 | 黄色视频网站免费在线观看 | 久久久久久久久99精品情浪 | 一级卖婬片AAAA片爱奴 | 午夜探花 | 欧美13一16娇小xxxx | 网站黄色小视频 | 美女18禁秘 啪啪免费看 | 国产一级a毛一级a看高清视视频 | 骚逼尤物 | 日本三级bd高清 摸进她的内裤里疯狂揉她韩国电影 | 亚洲精品一区中文字幕 | 国产偷窥在线观看视频 | 女人18毛片aa**水真多 | 靠逼啪啪视频 | 免费超碰在线观看 | 精品无人国产偷自产在线 | 厨房里边做饭边啪啪爱爱小说 | 久久久久久免费毛片 | 在线播放青青草 | 久草视频福利 | 91插逼 | 青青草午夜 | 日韩激情网 | rihanhuangpian | 亚洲欧美日韩精品在线 | 国产一级特黄录像片 | 国产性网| 人人操久久| 口述岳下面好紧我受不 | 91资源在线观看 | 性处破与摘花过程videos | 亚洲AV三级片 | 超碰大香蕉在线 | 天天日小穴天天插小穴天天干小穴 | 亚洲中文字幕无码中文字无摸 | 国产精品久久久在线观看 | www.操操 | 一级a做一级a做片性视频爱小说 | 国产日本亚洲香蕉视频 | 欧美日逼| 在线欧美不卡 | 国产精品亚洲五月天丁香 | 欧美.—级 片内射欧美老妇 | 性爱在线视频 | 青娱乐99| 护士美妇疯狂迎合娇吟 | 日韩午夜无码中文字幕A片评价 | 国产精品久久久视频 | 青青青爽无码视频在线观看一 | 免费看男男gay啪啪 | 迷人女教师伦理k8 | 大鸡吧插入 | 欧美性受XXXX黑人XYX性爽冫 | 天天操比| 精品无码成人免费看视频 | 亚洲成人色情影视 | 人妻二区 | 轻点好疼你的太大了 | 黄色三级小说 | 高清无码成人免费在线 | 日本一级做a爰片毛片地址 | 天天综合网~永久入口 | 亚洲AⅤ无码一区二区波多野按摩 | www.三级片 | 大鸡巴伊人网 | 性色av蜜臀av高清久久苏然 | 蜜桃视频无码 | ass白嫩白嫩的少妇pics | 光棍影院艳母 | 99热这里只有精品1 | 超碰在线 人人 | 三级小毛片| 亚洲欧洲小说图片在线免费观看 | 超碰97综合 | 肏屄毛片| 亚洲成年AV在线观看 | 欧美大黑屄 | 黄色一级大片免费版 | 国产美女久久久 | 操比视频 | 黄色内射 | 91偷拍视频网 | 奇米影视狠狠 | 安去也俺来了在线观看免费 | 淫乱一级| www视频在线免费观看 | 水多视频| 色亚洲天堂| 亚洲大尺度吃奶做爰 | 扒开双腿抽打花蒂惩罚小说 | 精品91自产拍在线观看二区 | 1314全毛片 | 国产精品久久久. | 嗯∽啊~轻点禁视频漫画 | 超碰东京 | 天天日女人| 久久久久久久国产免费看 | 啊灬啊灬啊灬快灬深用力游戏 | 日本高清不卡视频 | 久久综合新金瓶梅一级黄大片 | 操操操操操操网 | 国产精品成人a免费观看 | 久久粉嫩 | 国产最猛性xxxx | 国产性猛交XXXX免费看久久 | 毛片一区二区 | 4438成人网五月天在线播放 | 夜夜撸天天干 | 欧美日韩一二三 | 亚洲国产精品自拍 | 激情五月天开心 | www免费视频在线观看 | 又粗又硬又爽免费看视频 | 人妻免费视频 | 粉嫩小泬BBBB免费看WWWW | 91人妻人人澡人人爽人人 | 他的硕大进出粉嫩h校园 | 欧美麻豆久久久久久中文 | 天天射天天干天天日 | 青青草大香蕉在线 | 亚洲A级毛片 | 九一操逼视频 | 国产日韩精品在线 | 国产人人爽 | aaa大香蕉 | 伸进老师小内裤疯狂揉摸激吻小说 | 老师的被到爽羞羞动漫 | 紧窄粉嫩被粗大撑开 | 爱福利视频 | 特级西西WWW444人体聚色 | 大菠萝导航系统福建省委 | 胡桃与旅行者深入交流视频 | 国产91视频免费看 | 亚洲电影一级片 | 好叼操这里只有精品 | 婷婷丁香亚洲 | 家庭伦乱视频 | 我挺进她的下面疯狂运动 | 欧美另类综合 | 国产精品午夜福利一区二区 | 少妇bbbb水多毛片人 | 成人电影91 | 人妻资源网 | 三级无码在线播放 | 少妇的呻吟 | 国产伦精品一区二区三区高清版禁 | 午夜羞羞视频在线观看 | 韩国面膜胸口一面膜下的作用 | 我被各种男人啪h网球王子 | 女生扒开让男生捅 | 免费看涩涩视频软件 | 被男狂揉吃奶胸60分钟 | 国产精品人人爽人人做Av片眉药 | 精品一二三 | 国产精品久久久久久久久夜色 | 操逼黄视频 | 色操网 | 欧美色逼 | 婷婷色一二三区波多野结衣 | 和岳每晚弄的高潮嗷嗷叫视频 | 国产精品久久精品 | 影音先锋日韩无码 | 日产mv免费观看 | 竹菊av一区二区三区四区五区 | 亚洲+日产+专区 | 成人无码小视频 | 99久久精品免费看国产 | 看美女操屁一级黄色片 | 豆花视频网站在线观看 | 免费国产羞羞网站视频 | 色噜噜在线播放 | 成人免费无遮挡黄H网站在线观看 | 欧美拍拍| 人人插人人摸 | 三级综合网 | 午夜免费观看一区二区三区视频 | 蜜臀久久99精品久久久老牛影视 | 三级片网站视频 | 中文字幕18页 | 91麻豆成人精品国产 | 影音先锋你懂得 | 嗯啊h客厅hh青梅h半夜偷爬窗 | 欧美巨大黑人极品精男 | 日本理伦手机在线观看免费 | 操逼爽视频 | 亚洲性爱在线播放 | se999se | 啪啪啪啪啪啪视频 | 肏屄视频网 | 99视频免费观看 | 成年视频网站在线观看 | 伊人春色网 | 精品aaa| 男人添女人下面高潮视频 | 佐佐木希av在线播放 | 日本一区二区三区在线视频 | 理论片午午伦夜理片 | 黄色网免费 | 69精品无码一区二区三区 | 不要钱的黄色app | 盗摄偷窥系列7777777 | 扒开腿挺进女警湿润的花苞视频 | 黄色AV免费 | 成人免费无码婬片在线观看免费 | 中国无码电影 | 亚洲五月天激情 | 毛产一级婬片A片AAA片A国 | 国产精品一区二区三区毛片 | 18禁免费看片 | 成人免费看AA片 | 天堂福利影院 | 伊人大香焦网 | 亚洲天堂AV在线播放 | 欧美色www| 欧美三级性爱 | 快点cao我要被cao烂了小说 | 黄污涩涩啪啪www免费下载 | 成人精品 | 公牛与女人又大又爽 | 人妻免费视频 | www.色在线 | 逼逼久久 | 成人欧美一区二区三区无码电影 | 国产人伦激情在线观看 | 免费国偷自产拍精品视频 | 玩弄少妇高潮A片水蜜桃网站 | 亚洲天天搞 | 北条麻妃办公室性爱 | 成人无码区免费A片视频QQ群 | 韩国三级黄色毛片 | 榴莲视频 | 99久久久 | 深爱激情丁香五月天 | 日日擼夜夜擼狠狠擼88 | 让人看了下面流水的视频 | 国产卡一卡二卡三无线乱码新区 | 久久久免费看片 | chinesexxxxhd老头 | 中文字幕乱码一区二区三区四区在线看 | 白丝爆羞羞视频 | 男女网站免费观看 | 国产成人精品国内自产色情按摩 | 婷婷亚洲五月***久久 | 国产成人无码区免费AⅤ片老人 | 美女露出让男生揉的电影 | 天天人人综合 | 国产激情的老师在线播放 | 亚洲黄频| chinese国产18粉嫩dh | 亚洲欧美日韩成人高清在线一区 | 大尺度吻戏原声版免费观看 | 吻床戏软件大尺度18软件 | 国产美女精品人人做人人爽 | 日本乱伦无码视频 | sesese999 | 免费 成人 美女姛免费看 | 小早川怜子无码入口处 | 久久无码爱爱 | 肉色丝袜好紧我要进去了 | 91叫床高潮娇喘嗯啊在线观看 | 97香蕉久久国产超碰青草软件 | 做逼视频 | 日韩无码AV一区二区三区 | 日韩后入 | 免费观看性生交大片性教育 | 操逼AV网站 | 偷偷操影院 | 黄色一级A| 亚洲AV无码成人精品区麻豆曹妞 | 久久国产不卡 | 男女脱了内裤互蹭的视频 | 国产日韩二区 | 亚州色色999 | 女人高潮一级毛片 | 黄网视频 | 卫生间大尺度做爰床戏 | 婷婷五月小说 | 国产十大尺度真做无删减电影 | 青青草操 | 大尺度做爰戏 | 午夜精品一区二区三区视频免费看 | 三级片_小视频 - 百度 | 亚洲综合在线播放 | 91亚洲国产成人久久精品麻豆 | 又粗又长又硬爽死我了视频 | 师尊双性精跪趴灌满h微博 | 国产亲子伦一区二区三区四区 | 综合色吧 | 中文无码日本一级A片久久影视 | 嗯灬啊灬快灬高潮了喷水视频 | 高清无码视频网站在线观看 | 国产精品老熟女视频一区二区 | 波多野吉衣中文字幕 | 少妇警花的高潮叫床在线小说 | 日本在线看 | A级大黄片 | 国产精品自拍在线观看 | 欧美老妇XX | a∨视频在线 | 99久久99久久免费精桃 | www.奇米色.com | 日韩国产精品一区二区 | 91AV在线电影 | 欧美一级生活片 | 99国产精品久久久久久久夜 | 亚洲无码在线精品 | 99热思思热 | 伊人婷婷影院 | 欧美日韩国产精 | 噜噜噜噜噜久久久久久91黄画 | 久久久精品电影 | 欧美一区二区三区艳史 | 拍真实国产伦偷精品 | 热久久av| 亚洲欧美成人毛片 | 国产AⅤ无码片毛片一级一区2 | 亚洲精品乱码久久久久久蜜桃动漫 | 色婷婷丁香 | 调教小骚逼 | 2025天天操夜夜操 | 在线淫色 | 教师的s肉续txt | 国产精品成人a免费观看 | 欧美日韩大黄片 | 日韩操逼图 | 图书室里娇喘的短裙白丝校花视频 | 极品人妻疯狂3p超刺激 | 國产一二三内射在线看片 | 三级在线网 | 国产在线一区不卡 | 国内精品毛片免费看 | 猛男大粗猛爽h男人味69XX | 乱lun合集小可的奶水 | 九九视频免费在线观看 | 把女人到嗷嗷嗷叫视频 | chinesefreehdxxxx高潮 | 色色开心网 | 色尼玛图| 丁香五月天在线 | 日韩小穴 | jizzjizz中国少妇中文 | 欧美操逼视频网 | 免费做爱视频在线观看视频网站 | 国内爱爱视频 | 97网站| 妈妈喊我小老公 | 亚洲AV黄片| 涩涩屋导航福利av导航 | 鸥美超逼视频 | 色噜噜狠狠狠狠色综合 | 中文字幕的1级毛片视频 | 国产精品精品国产婷婷这里Aⅴ | 中文字幕人妻一区二区三区高清 | 中国免费一级片 | 操操操啊啊啊 | 国內毛片 | 91精品日产乱码一二三区别 | 被男狂揉吃奶胸高潮视频观看 | 国产电影自拍 | 淫色人妻综合 | 五月天激情小说 | 城市美妇呻吟高潮泄身 | 娇妻张妍交换高潮 | 亚洲黄片免费 | 开心五月综合网 | 国产午夜高清 | 毛片毛片毛片毛片 | 国产无遮挡AAA片爽爽 | 黄色成人免费视频 | 午夜视频在线观看一区二区 | 国产网站操逼 | 成人无码做爰www免费 | 黄色免费在线观看视频 | 胡桃给旅行者特殊的奖励视频 | 深夜黄色小视频 | 蜜桃 免费毛片 | 成人开心激情 | 97在线无码精品秘 入口果冻 | 天天日天天摸 | 国产熟妇 码AV | 非洲人做爰全过程视频 | 国产精品国产三级国产专区55 | 免费看操逼的 | 美女黄网站在线观看 | 美女激情网站 | 午夜成人无码精品免费 | 奇米网久久| 欧美三级韩国三级日本三斤 | 亚洲精品粉嫩小泬18p图片 | 亚洲无码免费视频一区二区三区四虎 | 日本大奶子视频 | 猛烈顶弄h禁欲律师1v1视频 | 99久久久无码国产精品免费 | ThePorn人妻白浆 | 影音先锋国产在线 | 在线看逼 | 九九色网 | 一级片在线免费观看视频 | 四虎日韩 | 免费无码婬片AAAA片软件 | 2019年天天夜夜 | 欧美一级片在线 | 丰满秘书办公室猛烈进出在办公室 | 操逼黄网站| 日韩精品一区二区免费 | 中国老妇人性爱网站 | 亚洲国产成人精品女人 | 淫色插插 | 操女人逼小说 | 亚洲黄色视频在线免费观看 | 性感美女啪啪 | 中文字幕一区二区三区四虎在线 | 变态操逼视频 | 久久婷婷五月天 | 色色网页 | 成人三级无码久久 | 伊人黄色电影 | 国产片一区二区 | 三区视频在线播放 | 91午夜伦伦电影理论片老牛吃嫩草 | 国产在线第一页 | 8x8x成人免费视频入口 | 无套猛操| 思思热视频在线观看 | 18禁免费看片 | 国产激情啪啪 | 日本无码精品一区三d里番精品 | 中国明星裸体无遮挡 | 成人无码高清在线观看 | 久久综合热 | 干丝袜少妇 | 肏屄黄色视频 | 人妻精品无码中文字幕一区二区三区 | 欧美精品无码一区二区三区仓井松 | 精品久久久久久久久久久国产字幕 | 五月丁香婷婷综合网 | 国产美女粉嫩无套内谢 | 久久黄色免费网站 | 91亚洲国产成人精品一区二区三 | 欧美一级人妻免费视频 | www免费视频在线观看播放 | ~慢点~嗯哼抱着揉c小说健身房 | 一级A片在线免费观看 | 啊轻点灬吻浑圆嗯太深了用力 | 91性爱影院 | 国产丝袜人妻一区二区三区电影 | 91精品久久久久久久久中文字幕 | 18 爽爽精品 国产 | 欧美色视频免费 | 国产一级在线播放 | 夫妻性生活裸体视频 | 懂色av懂色av粉嫩av | 久久三级影院 | 国产伦精品一区二区三区照片91 | 日韩一级a视频 | 日本做床爱全过程激烈视频 | 天天日姘 | 亚洲精品秘 一区二区三区影 | 女人18片毛片90分钟免费 | 99福利视频 | 日本边添边摸边做边爱小说 | 色色综合色色 | 美女床上黄片 | 91熟女柏欣彤大黑牛自慰 | 囯产精品久久久久久久久久辛辛 | 希岛爱理avhd101观看 | sm在线看 | 色噜噜狠狠躁夜夜躁人人爽免费 | 五月天色区| 天堂8在线视频 | 337p大胆啪啪私拍人体 | 国产精品白丝av嫩草影院 | 99久久久国产精品无码网爆 | 黄片免费小视频 | 夜夜春亚到洲 | 久久精品三级电影 | 欧美日韩大黄片 | 黄色A毛片 | 中文字幕15页 | 国产馆av | www.亚洲 | 99热这里只有精品1 | 干逼AV| 色五月激情综合深爱 | 国产饥渴妇女在线 | 激情精品国产 | 天天干天天撸天天日 | www.操逼.con | 国产阿v视频 | 豆花视频最新网页在线观看 | 操粉逼 | 黑人一级毛片 | 亚洲v视频 | 韩国大尺度情事 | 在线无码区 | 国产一级一区 | 少妇被爽到 | xxnx国产麻豆 | 少妇做爱网站 | 亚洲无码综合在线 | 丁香五月久久 | 国产又粗又长又硬又爽的在线视频播放 | www.大鸡巴 | 骚视频在线观看 | 老牛精品无码AV一区 | 国产成人综合亚洲 | 操小逼片 | 婷婷开心色四房播播免费 | 女人18片毛片60分钟翻译 | 国产精品成人久久久久A级 | 91人妻人人澡人人爽人人精品 | 国产黑逼| 在线播放 神尾舞 | 国产成人精品久久 | 日日夜夜精品免费视频 | 天天拍天天操天天干 | 日韩中文字幕在线免费观看 | 差差漫画线上入口页面弹窗 | 蜜桃Av噜噜一区二区三区四区 | 干逼视频免费看 | 男人和女人插 | 久久福利一区 | 日本黄色视频在线看 | 好涨别揉小豆豆了要尿了 | 人人黄色| 麻烦操逼网 | 波多野结衣一品二品免费观看AV | www.肏屄 | 亚洲高清视频在线观看免费 | 日日騒av无码 | 麻豆一级A片久久久乱码 | 欧美日韩国产一级片 | 女人扒开屁股让我添 | a√农村老熟女又粗又大 | 国产AV天堂 | 现代黄色小说 | 少妇久久久久久 | 小男生被绑在机器上榨精文 | 麻豆传媒换妻操逼图 | 国产特一级A片一二三区 | 国产精品三级久久久人妻 | 色色色色综合 | 麻豆精品秘 国产传媒电影 | 啊轻点灬太粗嗯太深了用力太子妃 | 国产在线青青草 | 国产亚洲第一页 | 张慧敏陈丽的最裸体视频 | 妖精视频一区二区三区 | 最近2019中文字幕mv第三季歌词 | 免费的av网站 | 欧美影院一区二区 | 人人搞人人插 | 日韩电影一区二区三区 | 国产一级婬片A片A毛片18禁 | 欧美性猛交乱大交3 | 交换做爰h文 | 久久久国产精 | 无码精品一区二区三区同学聚会 | www.肏 | 老外把我添高潮了视频 | 精品一区二区三区在线成人 | 午夜理伦一区在线视频 | 00后rapper潮水仙踪林的歌词 | 国产成人无码精品久在线观看 | 白丝调教高潮失禁丨vk | 国产精品三级在线 | www.婷色情网一区 | 熟女性爱视频 | 欧美二三四区 | 一级黄色毛片电影 | 少妇一次50分钟高潮 | 黄片免费在线观看av | 91丨豆花丨国产熟女 | 欧美老妇性猛交XXXX免费 | 免费观看黄色视频 | 91国产偷拍 | 91人人澡 | 成人做爰黄 片免费观看 | 偷偷碰偷偷鲁免费视频一区 | 久久五月天婷婷综合 | 日韩精品人妻在线高清不卡一区二区 | 色99久久 | 免费人成在线观看视频高潮 | 亚洲乱伦欧美乱伦 | 国产精品久久久久久久久久久新郎 | 亚洲AV成人精品一区二区三区四区 | 国产夫妻在线视频 | 精品人妻一区二区三区奶水 | a在线观看视频 | 国产日产高清欧美一区二区三区 | 成人黄色一区二区三区真人做爽爽爽视频 | 少妇啪啪调情流白浆 | xxxx日本在线观看 | 美国少妇xxxx内射 | 久热精品视频免费在线观看 | 成人午夜精品无码区毛片四季视频 | 国产三级网址 | 私人成片免费观看 | 唐朝艳妇性史性话 | 国产麻豆操逼视频 | 大香蕉噜噜噜 | 大鸡巴操老逼 | 国产www视频 | 激情文学一区二区三区 | 日韩午| 国产福利一二三区 | 午夜福利成人视频 | 天堂网地址 | 嗯~用力啊~嗯~c我~啊哈校园 | 久久夜色撩人精品国产户外下载 | 美女av在线播放 888av | 麻豆tv在线入口 | 91肥老熟妇伦子伦 | 脱光x网 | 日本奴役捆绑受虐狂xxxx | 一道本无码视频 | 脔到她哭宫交h | 吾爱插逼网 | 中文字幕观看 | 中文在线а√天堂官网 | 欧美久久天堂 | 一边吻一边摸双乳的电影 | 日韩一级品电影 | 在线视频一区少妇露脸福利在线 | 亚洲清色 | 久热只有精品 | 黄色成人视频在线免费观看 | 懂色av蜜臀av粉嫩av分 国产人妖TS重口系列 | 欧美成人乱轮 | 退休后日大女儿的最新情况 | 做爱免费视频免费网址 | 中文字幕乱码亚洲无线码小说 | 抓住她的双乳大力揉弄 | 亚洲AV无码国产精品 | 娇妻淫伦3亲p视频 | 懂色.av | 国产一级做a爰片在线看免费 | 淫色人妻小说 | 叶子楣三级露全乳视频 | 色99在线| 男女做爱免费网站 | 操逼无码视频 | 欧美大香蕉免费看 | 成人无码www免费视频在线播放 | 免费一级婬片A片AAA毛私人玩物 | 干美女网站| 手机在线看片日韩 | 天堂国产在线视频 | 好大好硬好深 | 婷婷深爱激情 | 无码做爱视频在线观看 | 女女h文(h)小说 | 亚洲性爱小说 | 色情五月婷婷 | 制服诱惑一区二区三区 | 国产亚洲 久一区二区 | 午夜天天操性色 | 成人片无码观看在线视频 | 成人香蕉视频 | 免费少妇荡乳情欲视频 | 女人操逼毛片 | 一二三四专区人妇 | 青春草免费视频在线 | 黄色一级AV | 日韩中文字幕免费观看 | 另类TS人妖一区二区三区 | freexxxxhdvideo中国 | 欧美成人猛片 | 成人黄片网站 | 日韩 人妻 精品 无码 | 户外勾搭av片 | 成人18禁在线 | 女女同性av片免费网站 | 鸡巴视频在线观看 | 非洲毛片| 日韩AV中文字幕在线播放 | freejapanxxxxhdvideo | 无码毛片一区二区无码视频 | 北条麻妃激情视频 | 国产成人无码精品aa免费 | 又大又粗又硬又爽又黄毛片视频 | 国产日韩视频在线 | 大鸡吧色www | 九一综合色| 青青草视频在线网站 | 绿巨人破解版 | 日韩高清AV在线 | 999小视频 | 淫欲代价电影 | 日韩欧美日韩欧美 | 搡女人真爽免费视频大全 | 蜜桃视频高清无码 | 少妇高潮21p | www.younv | 澳门无码视频 | 国产精品免费一区久久电影 | 波多野结衣做爱 | 亚洲少妇五月天 | 久久91av | 丁香五月在线观看视频 | 男ji大巴进入女人的视频 | 护士美妇疯狂迎合娇吟 | 男插女电影| 夜间福利视频 | 九九精品在线观看 | 亚洲丁香婷婷 | 精品人妻少妇嫩草AV蜜桃 | 全黄一级裸体 | 日本三级午夜理伦三级三 | 成人无码A级毛片免费播放 | 中日韩在线视频 | 这里只有精品视频 | 娇h女嗯啊黑壮汉好猛hhh视频 | 北岛玲日韩一区二区三区 | 亚洲一区二区视频 | 无套内谢少妇高潮毛片 | 在线看91 | 免费观看成人网站 | 色999anyequ | 亚洲伦理视频在线观看 | 揉bbb玩bbb少妇bbb | 99精品久久久久久久免费看蜜月 | 人人天天操 | 五月天成人网址 | 午夜爱爱影院 | 国产性爱免费视频 | 黄黄的网站 | 超碰人人人人一夲道 | 69操逼视频 | 18男女拍拍拍久久精品 | 国产精品秘 久久久久久奇米影视 | 男人日女人嗯嗯嗯嗯免费网站 | 麻酥酥福利视频 | 国产CD系列TS人妖视频 | 精品视频在线播放 | 国产黄色录像视频 | 天天av天天翘天天综合网 | 91久久丝袜国产露脸动漫 | 韩国精品无码 | 搡老熟女大熟了88AV一区二区 | 精品一二三四区 | 被艹到爽| 欧美日韩国产精品一区二区 | 夜月直播 | 永久免费看A人片无码精 | 3级毛片 日韩三级片在线 | 女开嫩苞精品第一 | 日本少妇做爰免费视频软件 | 巨肉黄暴辣文高h文np短视频 | 国产色无码视频在线观看 | 亚洲欧美精品 | 九色丨老熟女丨91啦 | 9l农村站街老熟女 | 国产91高潮叫床ThePorn | 我与一位67老妇的性往事 | 日韩精品人妻中文字幕一二三区 | 三级黄带片| 公交车地铁高h荡肉呻吟 | 韩国黄色免费在线观看 | 美女嫩逼逼| 啊灬啊灬啊灬岳小说 | 午夜天堂 | 欲妇荡岳丰满少妇在线阅读 | 毛片肏逼| 日韩精品一区二区三区熟女人妻网站 | 97人人爽| 日本护士囗交吞精 | 60一本久熟女中文 | 操逼播放| 亚洲天堂免费在线观看视频 | 欧美国产午夜 | 999国产精品自拍 | 中文字幕一区二区久久人妻网站 | 太子妃被强行打开双腿h | 99热在线精品播放 | 俺去夜 | 色色色色色色色色网站 | 国产午夜成人在线 | 屄无码| 曰韩毛片 | 国产女人水真多18毛片18精品视频 | 国产亚洲视频在线观看 |