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

    魚腥藻脂肪氧合酶熱穩(wěn)定性提高的分子動(dòng)力學(xué)模擬

    2017-02-08 07:42:29陸兆新別小妹趙海珍呂鳳霞
    食品科學(xué) 2017年2期
    關(guān)鍵詞:熱穩(wěn)定性突變體活力

    錢 輝,陸兆新,張 充,別小妹,趙海珍,呂鳳霞*

    魚腥藻脂肪氧合酶熱穩(wěn)定性提高的分子動(dòng)力學(xué)模擬

    錢 輝,陸兆新,張 充,別小妹,趙海珍,呂鳳霞*

    (南京農(nóng)業(yè)大學(xué)食品科技學(xué)院,江蘇 南京 210095)

    為提高魚腥藻脂肪氧合酶(Anabaena sp. PCC 7160 lipoxygenase,Ana-LOX)的熱穩(wěn)定性,探索其熱穩(wěn)定性機(jī)制。通過(guò)同源建模預(yù)測(cè)了Ana-LOX的三維結(jié)構(gòu),并通過(guò)分子動(dòng)力學(xué)模擬技術(shù)手段分別模擬了Ana-LOX于298 K和320 K條件下在水溶液中的運(yùn)動(dòng)軌跡,計(jì)算均方根偏移與均方根漲落后,以均方根漲落為熱穩(wěn)定性評(píng)價(jià)指標(biāo),結(jié)合作用力分析,預(yù)測(cè)不穩(wěn)定區(qū)域。結(jié)果表明:Ana-LOX的T94、I96、A325~Q327區(qū)域柔性較大,是導(dǎo)致Ana-LOX熱不穩(wěn)定的關(guān)鍵區(qū)域。通過(guò)I96V定點(diǎn)突變后,半衰期提高到野生型的2.5 倍,驗(yàn)證了原預(yù)測(cè)結(jié)果的有效性,從而為脂肪氧合酶分子的熱穩(wěn)定性改造提供了理論依據(jù)。

    脂肪氧合酶;同源建模;分子動(dòng)力學(xué)模擬;熱穩(wěn)定性

    脂肪氧合酶(Lipoxygenase,LOX,EC1.13.1.13)屬于氧化還原酶,廣泛存在于動(dòng)物、植物、藻類、真菌以及少數(shù)細(xì)菌中[1]。LOX能夠?qū)R淮呋趸许?順-1,4-戊二烯結(jié)構(gòu)的多元不飽和脂肪酸,形成具有共軛雙鍵的脂肪酸氫過(guò)氧化物[2],這些氫過(guò)氧化物性質(zhì)活潑,在食品加工中具有重要的應(yīng)用價(jià)值[3]。此外,許多不飽和脂肪酸的氫過(guò)氧化物是食品風(fēng)味物質(zhì)的前體,能夠通過(guò)酶解進(jìn)一步轉(zhuǎn)化為酮、醛和醇類,較化學(xué)合成香料具有更高的商業(yè)價(jià)值[4]。

    目前市場(chǎng)上商品級(jí)LOX產(chǎn)品大多數(shù)以富含LOX的大豆粉作為酶源。大豆中存在多種LOX同工酶,導(dǎo)致LOX的產(chǎn)品質(zhì)量不穩(wěn)定[5]。因此,構(gòu)建高產(chǎn)重組菌成為發(fā)酵法生產(chǎn)LOX的關(guān)鍵。本實(shí)驗(yàn)室前期已克隆了魚腥藻(Anabaena sp. PCC 7120)LOX基因,并在Escherichia coli和Bacillus subtilis中獲得了表達(dá)[6-7];后續(xù)研究發(fā)現(xiàn)其在面粉加工中有著重要的應(yīng)用價(jià)值[8],并進(jìn)一步闡釋了其強(qiáng)筋漂白機(jī)理[9-10]。但魚腥藻脂肪氧合酶(Anabaena sp. PCC 7120 lipoxygenase,Ana-LOX)熱穩(wěn)定性較差,限制了其在食品工業(yè)中的應(yīng)用。因此,需要對(duì)其熱穩(wěn)定性機(jī)制進(jìn)行解析,從而制定出相應(yīng)的策略來(lái)提高其熱穩(wěn)定性,如loop環(huán)的改造[11]、定點(diǎn)突變[12]、融合自組裝雙親短肽[13]等。

    近年來(lái),大量蛋白質(zhì)晶體結(jié)構(gòu)被解析及同源建模算法的不斷完善,使得利用計(jì)算機(jī)模擬來(lái)完成蛋白質(zhì)結(jié)構(gòu)與功能的分析工作、更高效地研究熱穩(wěn)定性機(jī)制并對(duì)其進(jìn)行分子改造成為可能。Ana-LOX的晶體結(jié)構(gòu)尚無(wú)報(bào)道,不宜采用基于晶體結(jié)構(gòu)解析來(lái)分析其熱穩(wěn)定性。因此,應(yīng)用計(jì)算機(jī)輔助技術(shù),采取同源建模、分子動(dòng)力學(xué)(molecular dynamics,MD)模擬等方法,以均方根偏移(root mean square deviation,RMSD)與均方根漲落(root mean square fluctuation,RMSF)值來(lái)評(píng)估酶分子局部空間結(jié)構(gòu)的改變,探索其熱不穩(wěn)定區(qū)域[14]及酶熱穩(wěn)定性機(jī)制。本實(shí)驗(yàn)將MD用于LOX的研究,鎖定Ana-LOX的熱不穩(wěn)定區(qū)域,探索酶熱穩(wěn)定性機(jī)制,為L(zhǎng)OX的熱穩(wěn)定性改造提供理論依據(jù)。

    1 材料與方法

    1.1 材料與試劑

    pET32a-Ana-LOX重組質(zhì)粒由南京農(nóng)業(yè)大學(xué)酶工程實(shí)驗(yàn)室構(gòu)建。

    LB培養(yǎng)基:胰蛋白胨10 g/L、酵母提取物5 g/L、NaCl 10 g/L、pH 7.0,121 ℃條件下滅菌20 min。

    限制性內(nèi)切酶Dnp Ⅰ、primeSTAR聚合酶 寶生物工程(大連)有限公司;氨芐青霉素、異丙基-β-D-硫代半乳糖苷(isopropyl β-D-thiogalactoside,IPTG) 生工生物工程(上海)股份有限公司;亞油酸 美國(guó)Sigma公司;其他試劑均為國(guó)產(chǎn)分析純。

    1.2 儀器與設(shè)備

    JY91ⅡDN超聲波破碎儀 寧波新芝生物科技股份有限公司;UV-2600紫外-可見分光光度計(jì) 日本島津公司。

    1.3 方法

    1.3.1 Ana-LOX蛋白三維結(jié)構(gòu)構(gòu)建

    將Ana-LOX基因序列提交phyre2[15]數(shù)據(jù)庫(kù)(http:// www.sbg.bio.ic.ac.uk/phyre2),與晶體結(jié)構(gòu)獲得解析的同源蛋白質(zhì)進(jìn)行比對(duì),將一致性最高的序列作為同源建模的模板,構(gòu)建目的蛋白質(zhì)的三維結(jié)構(gòu)。將構(gòu)建好的模型提交PSVS服務(wù)器(http://psvs-1_5-dev.nesg.org/),用Procheck[16]程序進(jìn)行模型評(píng)估,借助Pymol軟件進(jìn)行結(jié)構(gòu)分析[17]。

    1.3.2 Amber 12軟件分子動(dòng)力學(xué)模擬步驟及參數(shù)設(shè)置

    1.3.2.1 模擬體系構(gòu)建

    利用Amber[18]程序包和ff99SB[19]力場(chǎng),加載小分子gaff力場(chǎng)與金屬離子參數(shù),使用Leap模塊添加氫原子,選用TIP3P八面體水分子盒模型,蛋白質(zhì)分子置于10 ?的水溶劑箱內(nèi),加入Na+平衡電荷。

    1.3.2.2 能量最小化

    用sander模塊進(jìn)行3 000步能量最小化,執(zhí)行最陡下降法500 步后轉(zhuǎn)為共軛梯度法,以消除不合理的能量勢(shì)壘。

    1.3.2.3 加熱

    對(duì)系統(tǒng)加熱使其從0 K升溫至298 K或320 K,并在保持體積不變的狀態(tài)下運(yùn)行50 ps帶有位置限制的動(dòng)力學(xué),將非鍵相互作用閾值設(shè)為8 ?,使用弱耦合的算法來(lái)控制溫度。

    1.3.2.4 平衡

    在常壓條件下系統(tǒng)分別處在298 K和320 K條件下,進(jìn)行500 ps的分子動(dòng)力學(xué)模擬,以平衡系統(tǒng)。

    1.3.2.5 動(dòng)力學(xué)模擬

    恒溫恒壓系統(tǒng)下進(jìn)行MD模擬,所用步長(zhǎng)為2 fs,以每1 000 步采集一次能量與坐標(biāo)信息。分析動(dòng)力學(xué)模擬過(guò)程中RMSD值變化,再選擇RMSD值平衡階段的數(shù)據(jù),計(jì)算RMSF值,形成軌跡文件,最后通過(guò)VMD軟件分析動(dòng)力學(xué)模擬軌跡。

    1.3.3 定點(diǎn)突變

    為了檢驗(yàn)MD預(yù)測(cè)的準(zhǔn)確性,針對(duì)模擬預(yù)測(cè)位點(diǎn),采取質(zhì)粒聚合酶鏈?zhǔn)椒磻?yīng)(polymerase chain reaction,PCR)突變法進(jìn)行定點(diǎn)突變。引物為I96V-f(5’-ACCGAT GTTCCTGACAGCTACTTTTTAG-3’)和I96V-r(5’-TGT CAGGAACATCGGTGGGATCAGTA-3’)。PCR程序?yàn)椋?6 ℃預(yù)變性5 min;98 ℃變性10 s,55 ℃退火5 s,72 ℃延伸90 s,循環(huán)25 次;72 ℃延伸5 min。全質(zhì)粒PCR產(chǎn)物經(jīng)限制性內(nèi)切酶DpnⅠ消化后轉(zhuǎn)化E. coli BL21(DE3)。

    1.3.4 LOX活力測(cè)定

    采用紫外-分光光度法測(cè)酶活力,反應(yīng)體系3 mL:取pH 9.0的Tris-HCl緩沖液2.79 mL、酶液10 μL、1.733 mmol/L亞油酸鈉底物200 μL混合均勻后放入30℃水浴中并開始計(jì)時(shí),以1 min內(nèi)3 mL反應(yīng)體系在234 nm波長(zhǎng)處的吸光度增加0.001為一個(gè)酶活力單位(U)。

    1.3.5 酶學(xué)性質(zhì)表征

    1.3.5.1 半衰期測(cè)定

    IPTG誘導(dǎo)表達(dá)目的蛋白,經(jīng)鎳柱親和層析純化后透析去除咪唑。將突變酶與野生型酶液調(diào)整至相同質(zhì)量濃度,在50 ℃條件下保溫,每隔5 min取樣,然后置于冰上冷卻10 min,按照1.3.4節(jié)測(cè)定LOX酶的殘余活力。再以ln(殘余活力)對(duì)時(shí)間作圖,并進(jìn)行線性擬合,所得斜率即為酶在該溫度條件下的一級(jí)失活常數(shù)(kd)。酶在一定溫度條件下的半衰期(t1/2)按以下公式計(jì)算。

    1.3.5.2 酶相對(duì)活力對(duì)溫度的變化曲線

    將緩沖液加熱至15、20、30、35、40、45、50、55、60 ℃,然后在各溫度條件下加入10 μL待測(cè)酶液,測(cè)定酶活力。以最適反應(yīng)溫度條件下的酶活力為100%,以不同溫度條件下酶的相對(duì)活力作圖,從而獲得該酶相對(duì)活力對(duì)溫度的變化曲線。

    1.3.5.3 酶動(dòng)力學(xué)參數(shù)的測(cè)定

    在3 mL反應(yīng)體系中以亞油酸為底物,以濃度梯度0.047、0.023、0.019、0.014、0.009 mmol/L添加底物,進(jìn)行酶活力測(cè)定。以摩爾消光系數(shù)ε234nm為2.5×104(mol/L·cm)-1,計(jì)算產(chǎn)物生成量。采用Lineweaver-Burk雙倒數(shù)法作圖,計(jì)算動(dòng)力學(xué)參數(shù)[20]。

    2 結(jié)果與分析

    2.1 蛋白模型建立與評(píng)估

    圖1 Ana-LOX結(jié)構(gòu)與模板結(jié)構(gòu)比對(duì)Fig. 1 Comparison between Ana-LOX and template structure

    采用同源建模技術(shù)預(yù)測(cè)Ana-LOX三維結(jié)構(gòu)。Ana-LOX序列提交phyre2,選擇藍(lán)桿藻(Cyanothece sp.)PCC 8801脂肪氧合酶晶體結(jié)構(gòu)(PDB ID:5EK8)為模板[21],構(gòu)建成如圖1所示的三維空間主體圖。將模板與Ana-LOX在pymol中比對(duì),定位出鐵離子活性中心。在建模過(guò)程中,Ana-LOX與模板序列一致性為20%,低于30%,因此選擇phyre2數(shù)據(jù)庫(kù)建模,phyre2數(shù)據(jù)庫(kù)采取隱馬可夫模型比對(duì),并結(jié)合了從頭折疊模擬,使得即使在序列一致性較低(15%~25%)的情況下,依然可以提供精準(zhǔn)的結(jié)構(gòu)信息[15]。

    圖2 Ana-LOX拉氏構(gòu)象圖Fig. 2 Ramachandran plot of Ana-LOX

    表1 Procheck的拉氏構(gòu)象圖概要Table1 Ramachandran plot summary from Procheck

    采用Procheck對(duì)建好的模型結(jié)構(gòu)進(jìn)行評(píng)估,檢查最終模型的二面角構(gòu)象和立體化學(xué)參數(shù),獲得其拉氏構(gòu)象圖(圖2)。結(jié)果表明,僅有1.9%的氨基酸位于不合理區(qū)域(表1),由此可以判斷模型質(zhì)量較高,該模型可以用于下一步的分子動(dòng)力學(xué)模擬實(shí)驗(yàn)。

    2.2 RMSD與RMSF值計(jì)算

    RMSD值是蛋白質(zhì)每個(gè)原子運(yùn)動(dòng)偏移的均值,反映蛋白質(zhì)整體的運(yùn)動(dòng)狀況。RMSD值平衡,反映蛋白質(zhì)此時(shí)運(yùn)動(dòng)處于平衡狀態(tài),此時(shí)即可抽取這一時(shí)間段來(lái)計(jì)算RMSF值。298 K條件下MD 1 ns后Ana-LOX就迅速進(jìn)入了平衡狀態(tài)(圖3a),取3~8 ns計(jì)算RMSF值。320 K條件下15 ns的MD過(guò)程中發(fā)現(xiàn)6~9 ns時(shí)RMSD值發(fā)生了躍遷,躍遷后又下降到達(dá)平衡狀態(tài)(圖3b)。RMSD值的躍遷說(shuō)明在6~9 ns時(shí),蛋白質(zhì)的構(gòu)象發(fā)生了較大的改變,然后又恢復(fù)平衡狀態(tài),在這一過(guò)程中發(fā)生構(gòu)象較大變化的位置可能是潛在的不穩(wěn)定的臨界位點(diǎn),分別取6~9 ns和10~15 ns計(jì)算RMSF值。

    圖4 Ana-LOX在298 K(3~8 ns )(a)與320 K(b)條件下RMSF值Fig. 4 RMSF values for Ana-LOX at 298 (a) and 320 K (b)

    如圖4所示,RMSF值高的區(qū)域振幅較大,可能是潛在柔性較大的不穩(wěn)定區(qū)域。N端的氨基酸RMSF值最高,但這是因N端位于蛋白質(zhì)末端的一側(cè)處于游離狀態(tài),導(dǎo)致天然振幅較大,不將其視為不穩(wěn)定區(qū)域。320 K條件下6~9 ns RMSD發(fā)生了躍遷,表明在MD過(guò)程中,存在部分位置氨基酸發(fā)生較大的位置偏移,使得RMSD發(fā)生躍遷。比較320 K條件下6~9 ns與平衡后的RMSF,選擇發(fā)生較大偏移的位點(diǎn)。選擇圖4a中出現(xiàn)RMSF峰值的氨基酸殘基位點(diǎn)與圖4b中6~9 ns條件下出現(xiàn)RMSF值躍遷的殘基位點(diǎn)作為候選,逐一分析各氨基酸與周圍氨基酸氫鍵、離子鍵、π鍵堆積等相互作用,確定T94、I96、A325~Q327區(qū)域(箭頭所示)為Ana-LOX不穩(wěn)定區(qū)域(圖5)。

    圖5 Ana-LOX不穩(wěn)定區(qū)域的位置Fig. 5 Location of the unstable regions of Ana-LOX

    2.3 突變體性質(zhì)鑒定

    圖6 Ana-LOX與I96V突變熱失活動(dòng)力學(xué)Fig. 6 Thermal inactivation kinetics of Ana-LOX and I96V mutant

    為了驗(yàn)證預(yù)測(cè)的準(zhǔn)確性,選擇MD預(yù)測(cè)的I96號(hào)位點(diǎn)突變成纈氨酸。突變體I96V經(jīng)純化,酶熱處理每隔5 min取樣測(cè)酶活力,以ln(殘余活力)(y)對(duì)時(shí)間(x)作圖,并進(jìn)行線性擬合(圖6)。野生型為y=-0.298 7x+4.515,R2=0.995 6;突變體I96V為y=-0.121 7x+4.686 6,R2=0.993。

    經(jīng)計(jì)算,野生型t1/2為2.3 min,突變體t1/2為5.7 min,是野生型的2.5 倍。由此判斷MD模擬對(duì)酶熱敏感區(qū)域的預(yù)測(cè)有效,可進(jìn)一步對(duì)預(yù)測(cè)出的T94、I96、A325~Q327區(qū)域進(jìn)行局部的剛性增強(qiáng),如定點(diǎn)突變引入更多的作用力或刪除柔性部分,以提高Ana-LOX的熱穩(wěn)定性。

    圖7 溫度對(duì)AAnnaa-LOX與I96V酶相對(duì)活力的影響Fig. 7 Effect of temperature on the activity of Ana-LOX and I96V mutant

    由圖7可知,I96位發(fā)生突變后,溫度對(duì)酶相對(duì)活力的影響發(fā)生了較大改變,在25~55 ℃范圍內(nèi)突變體可以保持80%以上的活力,而野生型對(duì)溫度的變化較為敏感,僅在33~37 ℃范圍內(nèi)酶相對(duì)活力可以維持80%以上。說(shuō)明定點(diǎn)突變拓寬了突變體I96V對(duì)溫度的適應(yīng)范圍。

    表2 動(dòng)力學(xué)參數(shù)測(cè)定Table2 Kinetic parameters of AAnnaa-LOX and I96V mutant

    與野生型相比,突變體I96V在熱穩(wěn)定性提高的同時(shí),酶的比活力也得到了提高,發(fā)生了活力與穩(wěn)定性的共進(jìn)化。如表2所示,突變體I96V的米氏常數(shù)(Km)相對(duì)野生型顯著降低,因此突變體I96V與底物更具有親和力,同時(shí)催化速率(kcat)均與kcat/Km高于野生型,比活力達(dá)到了野生型的2.48 倍。

    3 討 論

    影響蛋白質(zhì)熱穩(wěn)定性的因素有疏水相互作用、氫鍵、離子鍵、芳香環(huán)的相互作用、二硫鍵、金屬離子、氨基酸組成和外界因素等[22]。但酶的熱穩(wěn)定性往往受到諸多因素的協(xié)同影響,因此對(duì)于酶熱穩(wěn)定性的改造和機(jī)制研究一直是生物學(xué)和化學(xué)領(lǐng)域的難點(diǎn)。MD提供了在原子水平上探索蛋白質(zhì)動(dòng)力學(xué)特性的有力方法,作為一種計(jì)算機(jī)輔助手段,已被用于鑒定蛋白質(zhì)中可以作為增強(qiáng)穩(wěn)定性靶點(diǎn)的柔性區(qū)域?,F(xiàn)研究已表明,將分子動(dòng)力學(xué)模擬應(yīng)用到蛋白質(zhì)熱穩(wěn)定性的研究中,可以明顯提高蛋白質(zhì)的熱穩(wěn)定性[23]。Pikkemaat等[24]通過(guò)MD鑒定出脫鹵素酶的一段helix-loop-helix的高柔性區(qū)域,引入二硫鍵后使熱穩(wěn)定性及尿素變性都產(chǎn)生了巨大改變。高新星等[25]鑒定出枯草芽孢桿菌氨肽酶存在非催化的loop區(qū)域,并對(duì)其刪除,使半失活溫度(T50)從71 ℃提高到77 ℃。而有關(guān)魚腥藻脂肪氧合酶熱穩(wěn)定性提高的分子動(dòng)力學(xué)模擬研究國(guó)內(nèi)外尚未見報(bào)道。

    本研究為了比較298 K和370 K運(yùn)動(dòng)狀態(tài)下的構(gòu)象差異,分別在298 K和370 K條件下進(jìn)行了分子動(dòng)力學(xué)模擬;在MD過(guò)程中,當(dāng)活性中心缺少鐵離子時(shí),RMSD無(wú)法平衡,蛋白結(jié)構(gòu)發(fā)生崩塌(數(shù)據(jù)未給出),表明LOX活性中心的鐵離子對(duì)維持結(jié)構(gòu)穩(wěn)定、發(fā)揮催化活性不可缺少。而同源建模構(gòu)建的模型并沒有給出鐵離子信息,因此需要將模板與Ana-LOX模型比對(duì),定位出鐵離子位置。同時(shí),在MD過(guò)程中需要添加小分子力場(chǎng)與金屬離子參數(shù)。結(jié)果顯示,在298 K條件下1 ns后RMSD值就達(dá)到平衡,而320 K條件下6~9 ns之間RMSD值發(fā)生了躍遷后又恢復(fù)平衡(圖3)。整體上320 K條件下的RMSD高于298 K條件下的RMSD,RMSD變化較大證明在320 K 較高溫度下,主鏈Cα變化較為劇烈,分子柔性變大。比較兩個(gè)溫度下的RMSF,以及320K條件下6~9 ns和平衡時(shí)的RMSF,選取RMSF較大的氨基酸,并逐一結(jié)合蛋白質(zhì)結(jié)構(gòu)與作用力分析,預(yù)測(cè)出T94、I96、A3~Q327區(qū)域?yàn)锳na-LOX不穩(wěn)定區(qū)域(圖5),均位于Loop環(huán),且位于蛋白質(zhì)表面。研究已表明,影響蛋白質(zhì)熱穩(wěn)定性的氨基酸也主要位于蛋白質(zhì)表面[26],Loop環(huán)在蛋白質(zhì)中未形成有序二級(jí)結(jié)構(gòu)的部分,是柔性較大的區(qū)域,不利于蛋白質(zhì)的穩(wěn)定[27]。另外,從氨基酸組成角度分析,5 個(gè)氨基酸(T94、I96、A3、I326、Q327)中有4 個(gè)屬于疏水氨基酸,而疏水氨基酸一般分布于蛋白質(zhì)內(nèi)部,蛋白表面的疏水氨基酸不利于其穩(wěn)定性[28]。其中蘇氨酸(threonine,Thr)與谷氨酰胺(glutamine,Gln)也屬于熱不穩(wěn)定氨基酸,Thr可以和蛋白質(zhì)周圍的水分子發(fā)生相互作用,但在高溫條件下結(jié)合的水分子會(huì)釋放出來(lái),從而破壞了蛋白質(zhì)和水分子結(jié)合部位的局部結(jié)構(gòu),導(dǎo)致整個(gè)蛋白質(zhì)的失穩(wěn)[29]。而暴露狀態(tài)的Gln,在高溫條件下容易產(chǎn)生脫酰氨基反應(yīng),也不利于蛋白質(zhì)空間構(gòu)象保持穩(wěn)定[30]。

    綜上所述,本研究運(yùn)用Amber 12軟件,對(duì)Ana-LOX進(jìn)行了MD模擬,計(jì)算了蛋白質(zhì)分子在水溶液中的運(yùn)動(dòng)軌跡,發(fā)現(xiàn)Ana-LOX的T94、I96、A3~Q327區(qū)域柔性較大,為Ana-LOX不穩(wěn)定區(qū)域,并對(duì)I96號(hào)位進(jìn)行了突變,驗(yàn)證了預(yù)測(cè)的有效性,探索了增強(qiáng)酶熱穩(wěn)定性的機(jī)制,結(jié)果表明通過(guò)蛋白質(zhì)工程分子改造手段,可以增強(qiáng)上述區(qū)域的剛性,提高Ana-LOX熱穩(wěn)定性。

    [1] HANSEN J, GARRETA A, BENINCASA M, et al. Bacterial lipoxygenases, a new subfamily of enzymes? a phylogenetic approach[J]. Applied Microbiology and Biotechnology, 2013, 97(11): 4737-4747. DOI:10.1007/s003-013-4887-9.

    [2] ALBERTI J C, MARIANI M, de CARAFFA V B B, et al. A functional role identif ed for conserved charged residues at the active site entrance of lipoxygenase with double specificity[J]. Journal of Molecular Catalysis B: Enzymatic, 2016, 123: 167-173. DOI:10.1016/ j.molcatb.2015.10.011.

    [3] ZHANG Chong, ZHANG Shuang, LU Zhaoxin, et al. Effects of recombinant lipoxygenase on wheat flour, dough and bread properties[J]. Food Research International, 2013, 54(1): 26-32. DOI:10.1016/j.foodres.2013.05.0.

    [4] GOUNARIS Y. Biotechnology for the production of essential oils, flavours and volatile isolates. a review[J]. Flavour and Fragrance Journal, 2010, 25(5): 367-386. DOI:10.1002/ffj.1996.

    [5] HUGHES R K, WU Z C, ROBINSON D S, et al. Characterization of authentic recombinant pea-seed lipoxygenases with distinct properties and reaction mechanisms[J]. Biochemical Journal, 1998, 333: 33-43.

    [6] 張充, 周孝偉, 呂鳳霞, 等. 重組魚腥藻脂肪氧合酶基因的克隆表達(dá)、分離純化及活性分析[J]. 生物工程學(xué)報(bào), 2012, 28(4): 440-456. DOI:10.13345/j.cj b.2012.04.006.

    [7] ZHANG C, TAO T, YING Q, et al. Extracellular production of lipoxygenase from Anabaena sp. PCC 7120 in Bacillus subtilis and its effect on wheat protein[J]. Applied Microbiology and Biotechnology, 201 2, 94(4): 949-958. DOI:10.1007/s003-012-3895-5.

    [8] 任娣, 謝亞娟, 陸兆新, 等. 重組脂肪氧合酶對(duì)面團(tuán)流變性質(zhì)及面包品質(zhì)的影響[J]. 食品科學(xué), 2015, 36(13): 1-6. DOI:10.7506/spkx1002-6630-201513001.

    [9] WANG X M, LU F X, ZHANG C, et al. Peroxidation radical formation and regiospecif city of recombinated Anabaena sp. lipoxygenas e and its effect on modifying wheat proteins[J]. Journal of Agricultural and Food Chemistry, 2014, 62(7): 1713-1719. DOI:10.1021/jf503238h.

    [10] WANG X M, LU F X, ZHANG C, et al. Effects of recombinated Anabaena sp. lipoxygenase on the protein component and dough property of wheat f our[J]. Journal of Agricultural and Food Chemistry, 2014, 62(40): 9885-9892. DOI:10.1021/jf503238h.

    [11] LU X, LIU S, FENG Y, et al. Enhanced thermal stability of Pseudomonas aeruginosa lipoxygenase through modif cation of two highly f exible regions[J]. Applied Microbiology and Biotechnology, 2014, 98(4): 1663-1669. DOI:10.1007/s00253-013-5039-y.

    [12] GUO F F, ZHANG C F, BIE X M, et al. Improving the thermostability and activity of lipoxygenase from Anabaena sp. PCC 7120 by d irected evolution and site-directed mutagenesis[J]. Journal of Molecular Catalysis B: Enzymatic, 2014, 107: 23-30. DOI:10.1016/ j.molcatb.2014.05.016.

    [13] LU X, LIU S, ZHANG D, et al. Enhanced thermal stability and specif c activity of Pseudomonas aer uginosa lipoxygenase by fusing with self-assembling amphipathic peptides[J]. Applied Microbiology and Biotechnology, 2013, 97(21): 9419-9427. DOI:10.1007/s00253-013-4751-y.

    [14] CHRISTENSEN N J, KEPP K P. Stability mechanisms of a thermophilic laccase probed by molecular dynamics[J]. PLoS ONE, 2013, 8(4): e61985. DOI:10.1371/journal.pone.0061985.

    [15] KELLEY L A, MEZULIS S, YATES C M, et al. The Phyre2 web porta l for protein modeling, prediction and analysis[J]. Nature Protocols, 2015, 10(6): 845-858. DOI:10.1038/nprot.2015.053.

    [16] LASKOWSKI R A, MACARTHUR M W, MOSS D S, et al. Procheck-a program to check the stereochemical quality of protein struct ures[J]. Journal of Applied Crystallography, 1993, 26(2): 26283-26291. DOI:10.1107/S0021889892009944.

    [17] SEELIGER D, de GROOT B L. Ligand docking and binding site analysis with PyMOL and Autodock/Vina[J]. Journal of Computer-Aided Molecular Design, 2010, 24(5): 417-422. DOI:10.1007/s10822-010-9352-6.

    [18] SALOMON-FERRER R, CASE D A, WALKER R C. An overview of the Amber biomolecular simulation package[J]. Wiley Interdisciplinary Reviews-Computational Molecular Science, 2013, 3(2): 198-210. DOI:10.1002/wcms.1121.

    [19] L INDORFF-LARSEN K, PIANA S, PALMO K, et al. Improved sidechain torsion potentials for the Amber ff99SB protein force field[J]. Proteins-Structure Function and Bioinformatics, 2010, 78(8): 1950-1958. DOI:10.1002/prot.22711.

    [20] BRODHUN F, CRISTOBAL-SARRAMIAN A, ZABEL S, et al. An iron 13S-lipoxygenase with an α-linolenic acid specific hydroperoxidase activity from Fusarium oxysporum[J]. PLoS ONE, 2013, 8(5): e64919-e64919. DOI:10.1371/journal.pon e.0064919.

    [21] NEWIE J, ANDREOU A, NEUMANN P, et al. Crystal structure of a lipoxygenase from Cyanothece sp. may reveal novel features for substrate acquisition[J]. Journal of Lipid Res earch, 2016, 57(2): 276-287. DOI:10.1194/jlr.M064980.

    [22] 田健, 王平, 伍寧豐, 等. 理性設(shè)計(jì)提高蛋白質(zhì)熱穩(wěn)定性的研究進(jìn)展[J]. 生物技術(shù)進(jìn)展, 2012, 2(4): 233-239. DOI:10.3969/ j.issn.2095-2341.2012.04.01.

    [23] ZHANG C, ZHANG L, ZHANG Y, et al. Study on the stabil ity and enzymatic property improvement of porcine pancreas lipase modif ed by ionic liquids using molecular simulation[J]. Acta Chimica Sinica, 2016, 74(1): 74-80. DOI:10.6023/a15070496.

    [24] PIKKEMAAT M G, LINSSEN A B, BERENDSEN H J, et al. Molecular dynamics simulations as a tool for improving protein stability[J]. Protein Engineering, 2002, 15(3): 185-192. DOI:10.1093/ protein/15.3.185.

    [25] GAO X X, LIU Z M, CU I W J, et al. Enhanced thermal stability and hydrolytic ability of Bacillus subtilis aminopeptidase by removing the thermal sensitive domain in the non-catalytic region[J]. PLoS ONE, 2014, 9(3): e92357. DOI:10.1371/journal.pone.0092357.

    [26] MORLEY K L, KAZLAUSKAS R J. Improving enzyme properties: when are closer mutations better?[J]. Trends in Biotechnology, 2005, 23(5): 231-237. DOI:10.1016/j.tibtech.2005.03.005.

    [27] LI B C, YANG G Y, WU L, et al. Role of the N C-Loop in catalytic activity and stability in lipase from Fervidobacterium changbaicum[J]. PLoS ONE, 2012, 7(10): e46881. DOI:10.1371/journal.pone.0046881.

    [28] SCHWEHM J M, KRISTYANNE E S, BIGGERS C C, et al. Stability effects of increasing the hydrophobicity of solvent-exposed side chains in staphylococcal nuclease[J]. Biochemistry, 1998, 37(19): 6939-6948. DOI:10.1021/b i9725069.

    [29] NAGENDRA H G, SUKUMAR N, VIJAYAN M. Role of water in plasticity, stability, and action of proteins: the crystal structures of lysozyme at very low levels of hydration[J]. Proteins-Structure Function and Bioinformatics, 1998, 32(2): 229-240. DOI:10.1002/ (SICI)1097-0134(19980801)32:2<229::AID-PROT9>3.0.CO;2-F.

    [30] KANEKO H, MINAGAWA H, SHIMADA J. Rational design of thermostable lactate oxidase by analyzing quaternary structure and p revention of deamidation[J]. Biotechnology Letters, 2005, 27(22): 1777-1784. DOI:10.1007/s10529-005-3555-2.

    Molecular Dynamics Simulation of Thermal Stability Improvement of Lipoxygenase from Anabaena sp. PCC 7160

    QIAN Hui, LU Zhaoxin, ZHANG Chong, BIE Xiaomei, ZHAO Haizhen, L? Fengxia*
    (College of Food Science and Technology, Nanjing Agricultural University, Nanjing 210095, China)

    Homology modeling was used to predict the 3D structure of lipoxygenase from Anabaena sp. PCC 7160 (Ana-LOX) and to simulate its motion trajectory in aqueous solution by molecular dynamics simulation at 298 and 320 K, respectively. After root mean square deviation (RMSD) and root mean square f uctuation (RMSF) values were calculated, thermal stability was evaluated based on the RMSF value. Force analysis was also performed to predict the unstable region. Results indicated that the T94, I96, and A325 to Q327 areas of Ana-LOX were f exible and they played a crucial role in the thermal instability of Ana-LOX. The half-life of I96V mutant constructed by site-directed mutagenesis increased 2.5 folds compared with that of the wild-type enzyme, conf rming the validity of the prediction. Our research can provide a theoretical basis for the improvement of the thermal stability of lipoxygenase.

    lipoxygenase; homology modeling; molecular dynamics simulation; thermal stability

    10.7506/spkx1002-6630-201702001

    TS202.3

    A

    1002-6630(2017)02-0001-06

    錢輝, 陸兆新, 張充, 等. 魚腥藻脂肪氧合酶熱穩(wěn)定性提高的分子動(dòng)力學(xué)模擬[J]. 食品科學(xué), 2017, 38(2): 1-6. DOI:10.7506/spkx1002-6630-201702001. http://www.spkx.net.cn

    QIAN Hui, LU Zhaoxin, ZHANG Chong, et al. Molecular dynamics simulation of thermal stability improvement of lipoxygenase from Anabaena sp. PCC 7160[J]. Food Science, 2017, 38(2): 1-6. (in Chinese with English abstract)

    10.7506/spkx1002-6630-201702001. http://www.spkx.net.cn

    2016-03-30

    國(guó)家自然科學(xué)基金面上項(xiàng)目(31470095)

    錢輝(1990—),男,碩士研究生,主要從事食品生物技術(shù)研究。E-mail:2013108028@njau.edu.cn

    *通信作者:呂鳳霞(1963—),女,教授,博士,主要從事酶工程研究。E-mail:lufengxia@njau.edu.cn

    猜你喜歡
    熱穩(wěn)定性突變體活力
    活力
    CLIC1及其點(diǎn)突變體與Sedlin蛋白的共定位研究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    PVC用酪氨酸鑭的合成、復(fù)配及熱穩(wěn)定性能研究
    改制增添活力
    提高有機(jī)過(guò)氧化物熱穩(wěn)定性的方法
    收回編制 激發(fā)活力
    可聚合松香衍生物的合成、表征和熱穩(wěn)定性?
    Survivin D53A突變體對(duì)宮頸癌細(xì)胞增殖和凋亡的影響
    對(duì)羥基安息香醛苯甲酰腙的合成、表征及熱穩(wěn)定性
    九九久久精品国产亚洲av麻豆 | 精品久久久久久久久久久久久| 国产野战对白在线观看| 麻豆国产av国片精品| 少妇丰满av| 免费在线观看影片大全网站| 两性午夜刺激爽爽歪歪视频在线观看| 岛国在线免费视频观看| 婷婷六月久久综合丁香| 日韩国内少妇激情av| 中文字幕精品亚洲无线码一区| 少妇裸体淫交视频免费看高清| 五月伊人婷婷丁香| 悠悠久久av| 夜夜爽天天搞| 精品国产乱码久久久久久男人| 小说图片视频综合网站| 国产又黄又爽又无遮挡在线| 亚洲国产日韩欧美精品在线观看 | 两性午夜刺激爽爽歪歪视频在线观看| 国产精品亚洲av一区麻豆| 日本免费a在线| 丁香欧美五月| 婷婷亚洲欧美| svipshipincom国产片| 桃红色精品国产亚洲av| 久久中文字幕一级| 国产一级毛片七仙女欲春2| 特级一级黄色大片| 一级毛片高清免费大全| 一区二区三区国产精品乱码| 国产激情久久老熟女| 最近视频中文字幕2019在线8| 身体一侧抽搐| 日本黄大片高清| 亚洲人成网站在线播放欧美日韩| 757午夜福利合集在线观看| 国产成人精品久久二区二区免费| 黑人欧美特级aaaaaa片| 国产精华一区二区三区| 我要搜黄色片| 琪琪午夜伦伦电影理论片6080| 日韩精品青青久久久久久| 黑人欧美特级aaaaaa片| 日韩欧美国产在线观看| 老鸭窝网址在线观看| 午夜亚洲福利在线播放| 狠狠狠狠99中文字幕| 日韩欧美一区二区三区在线观看| 中出人妻视频一区二区| 亚洲午夜精品一区,二区,三区| 欧美日韩亚洲国产一区二区在线观看| 在线免费观看不下载黄p国产 | 黄色成人免费大全| 亚洲在线自拍视频| 亚洲在线自拍视频| 国产精品一区二区免费欧美| 日韩有码中文字幕| 亚洲中文av在线| 日本五十路高清| 高潮久久久久久久久久久不卡| 黄片大片在线免费观看| 亚洲精品美女久久久久99蜜臀| 国产成人av激情在线播放| 高清毛片免费观看视频网站| 18禁国产床啪视频网站| 国产91精品成人一区二区三区| 亚洲18禁久久av| 深夜精品福利| 男女午夜视频在线观看| 一区二区三区高清视频在线| 成人国产综合亚洲| 色老头精品视频在线观看| 免费无遮挡裸体视频| 99久久综合精品五月天人人| 成年人黄色毛片网站| 中文字幕av在线有码专区| 黑人巨大精品欧美一区二区mp4| 级片在线观看| 成人av在线播放网站| 成在线人永久免费视频| 国产又色又爽无遮挡免费看| 国产精品亚洲一级av第二区| 看免费av毛片| 久久久国产成人精品二区| 欧美成人一区二区免费高清观看 | 在线观看日韩欧美| 亚洲av第一区精品v没综合| 草草在线视频免费看| ponron亚洲| 欧美一区二区精品小视频在线| 午夜视频精品福利| 欧美日韩瑟瑟在线播放| 制服人妻中文乱码| 高清毛片免费观看视频网站| 亚洲精品国产精品久久久不卡| 丁香欧美五月| 欧美一区二区精品小视频在线| 午夜免费观看网址| www日本在线高清视频| 国产精品野战在线观看| www.熟女人妻精品国产| 精品午夜福利视频在线观看一区| 99久国产av精品| 精品午夜福利视频在线观看一区| 亚洲国产中文字幕在线视频| 欧美极品一区二区三区四区| 禁无遮挡网站| 久久亚洲真实| 久久久色成人| 免费搜索国产男女视频| 免费人成视频x8x8入口观看| 亚洲欧美日韩卡通动漫| 亚洲国产欧美人成| 老汉色av国产亚洲站长工具| 国产私拍福利视频在线观看| 999久久久国产精品视频| 麻豆久久精品国产亚洲av| 国产一区二区三区视频了| 老熟妇乱子伦视频在线观看| 免费看光身美女| 又黄又爽又免费观看的视频| 校园春色视频在线观看| 日本一二三区视频观看| 久久性视频一级片| 日本一二三区视频观看| 淫秽高清视频在线观看| 午夜两性在线视频| 久久久精品大字幕| 欧美中文日本在线观看视频| 午夜福利视频1000在线观看| or卡值多少钱| 久久中文看片网| 久久久国产欧美日韩av| 婷婷亚洲欧美| 国产成人精品久久二区二区免费| 色尼玛亚洲综合影院| 狂野欧美白嫩少妇大欣赏| 精品久久久久久久末码| 舔av片在线| 亚洲av免费在线观看| 男女下面进入的视频免费午夜| 亚洲中文字幕一区二区三区有码在线看 | 又紧又爽又黄一区二区| 国产精品野战在线观看| 亚洲av中文字字幕乱码综合| 亚洲成人久久爱视频| 1024手机看黄色片| 中国美女看黄片| 十八禁人妻一区二区| 熟女人妻精品中文字幕| aaaaa片日本免费| 桃色一区二区三区在线观看| 国产淫片久久久久久久久 | 精品无人区乱码1区二区| 欧美中文日本在线观看视频| 日本在线视频免费播放| 91麻豆av在线| 成人性生交大片免费视频hd| 桃色一区二区三区在线观看| e午夜精品久久久久久久| 久久久久性生活片| 久久精品国产99精品国产亚洲性色| 成人亚洲精品av一区二区| 老司机午夜十八禁免费视频| 成年免费大片在线观看| 法律面前人人平等表现在哪些方面| av片东京热男人的天堂| 精品午夜福利视频在线观看一区| 亚洲第一电影网av| 可以在线观看毛片的网站| 9191精品国产免费久久| 亚洲色图av天堂| 两性午夜刺激爽爽歪歪视频在线观看| 可以在线观看毛片的网站| 最近最新中文字幕大全电影3| 国产激情久久老熟女| 成年免费大片在线观看| 18禁黄网站禁片午夜丰满| 一本久久中文字幕| 女警被强在线播放| 日韩欧美免费精品| 午夜免费成人在线视频| 看免费av毛片| 午夜福利视频1000在线观看| 久久这里只有精品中国| 麻豆国产av国片精品| 久久久久亚洲av毛片大全| av中文乱码字幕在线| 亚洲成人精品中文字幕电影| 欧美在线黄色| 日韩欧美一区二区三区在线观看| 亚洲激情在线av| 久久中文字幕人妻熟女| 久久性视频一级片| www国产在线视频色| 久久久久久久久免费视频了| 亚洲天堂国产精品一区在线| 校园春色视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 99久久综合精品五月天人人| 精品久久久久久成人av| 亚洲国产精品久久男人天堂| 亚洲第一电影网av| 久久午夜综合久久蜜桃| 高潮久久久久久久久久久不卡| 亚洲av成人一区二区三| 国语自产精品视频在线第100页| 精品电影一区二区在线| 午夜免费观看网址| 狠狠狠狠99中文字幕| 久久久国产成人免费| 国产v大片淫在线免费观看| 亚洲第一电影网av| 欧美中文综合在线视频| 精品国内亚洲2022精品成人| 手机成人av网站| 亚洲专区国产一区二区| 我的老师免费观看完整版| 丁香欧美五月| 国产精品久久久久久精品电影| 男女做爰动态图高潮gif福利片| 国产午夜精品久久久久久| 中文字幕熟女人妻在线| 免费看美女性在线毛片视频| 久久久久国内视频| e午夜精品久久久久久久| 99热精品在线国产| 中文字幕久久专区| 国产精品99久久久久久久久| 一进一出抽搐gif免费好疼| 此物有八面人人有两片| 色在线成人网| 亚洲成人中文字幕在线播放| 午夜福利在线在线| 久久亚洲精品不卡| 嫩草影院精品99| 男插女下体视频免费在线播放| 又爽又黄无遮挡网站| 久久久久精品国产欧美久久久| 天堂影院成人在线观看| 久久热在线av| 日日夜夜操网爽| 国产亚洲精品综合一区在线观看| 国产伦一二天堂av在线观看| 中文资源天堂在线| 老司机深夜福利视频在线观看| 又黄又爽又免费观看的视频| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲av高清不卡| 天天一区二区日本电影三级| 免费无遮挡裸体视频| 午夜两性在线视频| 一进一出抽搐gif免费好疼| 午夜免费成人在线视频| 精华霜和精华液先用哪个| 久久精品国产亚洲av香蕉五月| 黄片小视频在线播放| 亚洲精品美女久久av网站| or卡值多少钱| 亚洲精品久久国产高清桃花| 一本久久中文字幕| 国产淫片久久久久久久久 | 两人在一起打扑克的视频| 精品国内亚洲2022精品成人| 小说图片视频综合网站| 国产主播在线观看一区二区| av视频在线观看入口| 亚洲人成电影免费在线| 亚洲av电影不卡..在线观看| 99热精品在线国产| 日本黄色视频三级网站网址| 在线免费观看的www视频| 欧美乱妇无乱码| 色尼玛亚洲综合影院| 天天一区二区日本电影三级| 美女大奶头视频| 亚洲狠狠婷婷综合久久图片| 99国产精品一区二区蜜桃av| 给我免费播放毛片高清在线观看| 久久久久精品国产欧美久久久| 又紧又爽又黄一区二区| 亚洲自拍偷在线| 婷婷六月久久综合丁香| 手机成人av网站| 在线a可以看的网站| 日韩欧美国产一区二区入口| 国产一区在线观看成人免费| cao死你这个sao货| 两人在一起打扑克的视频| 免费在线观看日本一区| 搡老熟女国产l中国老女人| 精品人妻1区二区| 国产精品一区二区三区四区久久| 国产成人影院久久av| 成年女人毛片免费观看观看9| 一个人看视频在线观看www免费 | 悠悠久久av| 嫩草影院精品99| 俄罗斯特黄特色一大片| 欧美xxxx黑人xx丫x性爽| 精品熟女少妇八av免费久了| 99久国产av精品| 99热精品在线国产| 岛国在线观看网站| 午夜福利高清视频| 婷婷六月久久综合丁香| 90打野战视频偷拍视频| 国模一区二区三区四区视频 | www日本在线高清视频| 免费在线观看影片大全网站| 一级a爱片免费观看的视频| 国产高清视频在线观看网站| 岛国在线免费视频观看| 欧美成狂野欧美在线观看| 男女那种视频在线观看| 国产视频一区二区在线看| 岛国视频午夜一区免费看| 男人和女人高潮做爰伦理| 丁香六月欧美| 久久久久久九九精品二区国产| 亚洲国产欧美人成| 在线观看美女被高潮喷水网站 | 丝袜人妻中文字幕| 性欧美人与动物交配| 啦啦啦韩国在线观看视频| 在线观看66精品国产| 国产精品久久久av美女十八| 国产成人福利小说| 午夜免费激情av| 精品久久久久久成人av| 无限看片的www在线观看| 男人舔女人的私密视频| 99精品欧美一区二区三区四区| 免费搜索国产男女视频| 亚洲欧美精品综合久久99| 激情在线观看视频在线高清| 久久久成人免费电影| 久久久色成人| 真人做人爱边吃奶动态| 美女大奶头视频| 最近最新中文字幕大全免费视频| 老司机福利观看| 制服人妻中文乱码| 成年版毛片免费区| 欧美av亚洲av综合av国产av| 亚洲熟妇中文字幕五十中出| 久久久国产成人免费| 成人精品一区二区免费| 动漫黄色视频在线观看| 国产主播在线观看一区二区| 国内精品美女久久久久久| 亚洲一区二区三区色噜噜| 88av欧美| 18禁黄网站禁片午夜丰满| 亚洲av成人精品一区久久| 九九在线视频观看精品| 19禁男女啪啪无遮挡网站| 18禁黄网站禁片午夜丰满| 中文字幕久久专区| 国产视频一区二区在线看| 国产免费av片在线观看野外av| 国产成人影院久久av| 日日夜夜操网爽| 国产又色又爽无遮挡免费看| 在线观看午夜福利视频| 少妇的丰满在线观看| 国内精品一区二区在线观看| 村上凉子中文字幕在线| 亚洲av熟女| 男女那种视频在线观看| 亚洲成人久久性| 琪琪午夜伦伦电影理论片6080| 男人的好看免费观看在线视频| 久久久久国产一级毛片高清牌| 丁香六月欧美| 国产高潮美女av| 国产欧美日韩精品亚洲av| 中国美女看黄片| 亚洲国产精品久久男人天堂| 又黄又爽又免费观看的视频| 亚洲欧美日韩无卡精品| 久久中文字幕一级| e午夜精品久久久久久久| 国产精品 欧美亚洲| 欧美绝顶高潮抽搐喷水| 99国产极品粉嫩在线观看| 亚洲中文av在线| 久久久久性生活片| 亚洲av中文字字幕乱码综合| 小说图片视频综合网站| 亚洲国产精品sss在线观看| 男女视频在线观看网站免费| 中文字幕av在线有码专区| 最好的美女福利视频网| 亚洲中文日韩欧美视频| 2021天堂中文幕一二区在线观| 三级国产精品欧美在线观看 | 亚洲精品美女久久av网站| 性色av乱码一区二区三区2| 在线播放国产精品三级| 色老头精品视频在线观看| 男人舔女人的私密视频| 欧美av亚洲av综合av国产av| 国产又色又爽无遮挡免费看| 国产精品女同一区二区软件 | 欧美三级亚洲精品| 少妇熟女aⅴ在线视频| 亚洲电影在线观看av| 亚洲熟妇熟女久久| 舔av片在线| 精品免费久久久久久久清纯| 色精品久久人妻99蜜桃| 亚洲狠狠婷婷综合久久图片| 日韩有码中文字幕| 国内久久婷婷六月综合欲色啪| 亚洲午夜精品一区,二区,三区| 在线观看免费视频日本深夜| 久久亚洲精品不卡| АⅤ资源中文在线天堂| 精品欧美国产一区二区三| 久99久视频精品免费| 最近最新中文字幕大全电影3| 一级黄色大片毛片| 久久久国产欧美日韩av| 亚洲性夜色夜夜综合| 男人和女人高潮做爰伦理| 99热这里只有精品一区 | 国内少妇人妻偷人精品xxx网站 | 脱女人内裤的视频| 日本免费一区二区三区高清不卡| 综合色av麻豆| 色播亚洲综合网| 99热这里只有精品一区 | 国产亚洲精品一区二区www| 欧美日韩黄片免| 精品国产乱码久久久久久男人| 国产精品 国内视频| 久久这里只有精品19| 怎么达到女性高潮| 我的老师免费观看完整版| 欧美日韩中文字幕国产精品一区二区三区| 国产免费av片在线观看野外av| 真人一进一出gif抽搐免费| 成熟少妇高潮喷水视频| 精品不卡国产一区二区三区| 在线视频色国产色| 麻豆成人午夜福利视频| 国产精品 欧美亚洲| 一二三四在线观看免费中文在| 999久久久国产精品视频| 欧美激情在线99| 天堂影院成人在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 日本 av在线| 亚洲色图 男人天堂 中文字幕| 黑人操中国人逼视频| 国产精品综合久久久久久久免费| 久久久久久大精品| 国模一区二区三区四区视频 | 国产精品一区二区免费欧美| 久久久久亚洲av毛片大全| 欧美激情在线99| 性欧美人与动物交配| 老熟妇乱子伦视频在线观看| 精品久久久久久,| 久久亚洲真实| 人妻久久中文字幕网| 欧美日韩综合久久久久久 | 中亚洲国语对白在线视频| 国产精品 欧美亚洲| 又粗又爽又猛毛片免费看| 舔av片在线| 啦啦啦韩国在线观看视频| 18禁黄网站禁片免费观看直播| 亚洲欧美一区二区三区黑人| 日本在线视频免费播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产欧美人成| 两性夫妻黄色片| 国产69精品久久久久777片 | 97人妻精品一区二区三区麻豆| 免费看美女性在线毛片视频| 久久久久性生活片| 色老头精品视频在线观看| www.自偷自拍.com| 高清毛片免费观看视频网站| 免费一级毛片在线播放高清视频| 免费在线观看影片大全网站| 性色av乱码一区二区三区2| 欧美极品一区二区三区四区| ponron亚洲| 国产一区二区在线av高清观看| 国内精品久久久久久久电影| 亚洲国产中文字幕在线视频| 在线观看66精品国产| 亚洲av成人不卡在线观看播放网| 亚洲七黄色美女视频| 亚洲在线观看片| 国产高清视频在线播放一区| 亚洲国产欧美一区二区综合| 国内久久婷婷六月综合欲色啪| 国产三级在线视频| 又大又爽又粗| 国产aⅴ精品一区二区三区波| 成人国产一区最新在线观看| 国产一区二区三区在线臀色熟女| 亚洲男人的天堂狠狠| 少妇熟女aⅴ在线视频| 国产一区二区激情短视频| 亚洲一区高清亚洲精品| 久久精品国产综合久久久| 宅男免费午夜| 久久这里只有精品19| 精品久久久久久久久久免费视频| 我的老师免费观看完整版| 视频区欧美日本亚洲| 亚洲欧美日韩高清专用| 亚洲无线观看免费| 午夜福利在线在线| a在线观看视频网站| 国产亚洲精品av在线| 久久久久亚洲av毛片大全| www.www免费av| 中文字幕久久专区| 一区福利在线观看| 久久久久久久午夜电影| 一本精品99久久精品77| 日韩欧美精品v在线| 日日干狠狠操夜夜爽| 淫妇啪啪啪对白视频| 嫩草影视91久久| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 欧美黑人欧美精品刺激| 国产亚洲欧美98| 一进一出抽搐gif免费好疼| 一个人免费在线观看的高清视频| 亚洲男人的天堂狠狠| 亚洲国产高清在线一区二区三| 色播亚洲综合网| 啦啦啦韩国在线观看视频| 国产精品免费一区二区三区在线| 久久久国产成人精品二区| 久久人妻av系列| 国产极品精品免费视频能看的| 舔av片在线| 看免费av毛片| ponron亚洲| 欧美xxxx黑人xx丫x性爽| 波多野结衣巨乳人妻| 在线永久观看黄色视频| 精品日产1卡2卡| 国产亚洲av高清不卡| 午夜亚洲福利在线播放| 国产 一区 欧美 日韩| 国产蜜桃级精品一区二区三区| 精品午夜福利视频在线观看一区| 我的老师免费观看完整版| 美女cb高潮喷水在线观看 | 在线观看免费视频日本深夜| 国产精品免费一区二区三区在线| 校园春色视频在线观看| 在线观看免费视频日本深夜| 天天躁日日操中文字幕| 韩国av一区二区三区四区| 精品国产三级普通话版| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 国产久久久一区二区三区| 最新中文字幕久久久久 | 欧美黄色淫秽网站| 又黄又粗又硬又大视频| 琪琪午夜伦伦电影理论片6080| 最新在线观看一区二区三区| 日韩免费av在线播放| 日本一二三区视频观看| 一级毛片精品| 97超视频在线观看视频| 国产成人av教育| 成人av一区二区三区在线看| 亚洲精品中文字幕一二三四区| 99久久精品国产亚洲精品| 老汉色∧v一级毛片| 我要搜黄色片| 看片在线看免费视频| 久久国产乱子伦精品免费另类| 老司机午夜十八禁免费视频| 国产高清视频在线播放一区| 51午夜福利影视在线观看| 成人永久免费在线观看视频| 噜噜噜噜噜久久久久久91| 精品久久蜜臀av无| 国产男靠女视频免费网站| 在线观看午夜福利视频| 黄色成人免费大全| 亚洲aⅴ乱码一区二区在线播放| 高潮久久久久久久久久久不卡| 香蕉av资源在线| 久99久视频精品免费| 国产精品一区二区三区四区免费观看 | 美女 人体艺术 gogo| 国内精品一区二区在线观看| 91av网站免费观看| 国产综合懂色| 宅男免费午夜| 制服人妻中文乱码| 亚洲 国产 在线| 国产麻豆成人av免费视频| 色老头精品视频在线观看| av视频在线观看入口| 亚洲午夜精品一区,二区,三区| 淫秽高清视频在线观看| 国产亚洲精品av在线|