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

    菊粉酶降低N-羥乙酰神經(jīng)氨酸含量機制的分子模擬研究

    2020-10-09 08:35:12羅先錕梁美蓮朱秋勁
    核農(nóng)學報 2020年9期
    關(guān)鍵詞:外切菊粉殘基

    常 瑞 羅先錕 何 江 梁美蓮 朱秋勁

    (1貴州大學釀酒與食品工程學院/貴州省農(nóng)畜產(chǎn)品貯藏與加工重點實驗室,貴州 貴陽 550025;2 貴陽市農(nóng)產(chǎn)品質(zhì)量安全監(jiān)督檢驗測試中心,貴州 貴陽 550004)

    紅肉及其制品是供給人體蛋白質(zhì)、多不飽和脂肪酸、B 族維生素、硒等優(yōu)質(zhì)營養(yǎng)素的來源。但近來諸多流行性病學和免疫學研究表明,日常飲食紅肉及其加工制品而攝入的 N - 羥乙酰神經(jīng)氨酸 (Nglycolylneuraminic,Neu5Gc) 可被機體識別產(chǎn)物抗-Neu5Gc 抗體,這種抗體誘導機體處于低度炎癥狀態(tài),并與結(jié)直腸癌等炎癥的發(fā)生相關(guān)[1]。Neu5Gc 屬于吡喃酮糖衍生物[2],是N-乙酰神經(jīng)氨酸 (Nacetylneuraminic acid,Neu5Ac) 被胞苷單磷酸-N-乙酰神經(jīng)氨酸羥化酶 (CMP-Neu5Ac hydroxylase,CMAH) 羥化的產(chǎn)物,其在人體中不能正常合成,主要依靠外源飲食豬肉、牛肉等紅肉攝入[3]。2018年美國牛肉消費量為1 218 萬t,中國豬肉消費量為5 540 萬t,且依然持續(xù)增漲[4],因此,Neu5Gc 成為重要的食品風險因子[5],安全高效降低Neu5Gc 在紅肉中的含量關(guān)乎紅肉制品的安全。

    目前已有不同方式改變紅肉中Neu5Gc 存在狀態(tài)和含量的研究,蔣蕓等[6]通過對豬肉和牛肉進行水煮、微波加熱、有機酸腌制處理使得肉樣中Neu5Gc 含量降低70%;β-半乳糖苷酶處理肉樣后,雖然可使樣品中Neu5Gc 含量降低84%,但這種方式僅破壞了Neu5Gc 與生物分子相連的糖苷鍵,使結(jié)合態(tài)變成了游離態(tài)。鼠傷寒沙門氏菌LT2 唾液酸酶對4-甲基傘形-Neu5Ac 的效果好于4-甲基傘形-Neu5Gc[7],但該酶無法用于食品工業(yè)。梁美蓮等[8]以十多種食品加工常見酶制劑處理肉樣和Neu5Gc 標準品的效果研究表明,外切型菊粉酶能夠降低牛肉中Neu5Gc 的含量,對純標準品含量降低率可達50.52%[8],表明菊粉酶有望成為工業(yè)化降低游離態(tài)Neu5Gc 的有效工具。

    利用分子模擬技術(shù)獲得底物與酶作用的殘基信息,有助于酶的理性設(shè)計和作用機制解析。楊倩等[9]為提高米根霉α-淀粉酶(ROAmy) 的熱穩(wěn)定性,基于分子動力學模擬結(jié)果,對該酶中的3 個氨基酸殘基G128、K269 和G393 進行了突變,獲得了熱穩(wěn)定性更好的突變體。目前已有X-衍射晶體法和分子動力學模擬對外切菊粉酶作用果糖催化位點解析[10]、同源建模分析菊粉酶功能域和活性位點[11]、菊粉酶對果糖六磷酸和蔗果三糖相互作用模擬的報道[12]。但有關(guān)Neu5Gc 與菊粉酶相互作用機制的研究卻鮮有報道。本研究對外切型和內(nèi)切型菊粉酶與Neu5Gc 復合物進行分子動力學模擬,采用密度泛函理論分析相關(guān)殘基與Neu5Gc 分子間作用的類型、強弱和分布,以期獲得菊粉酶與Neu5Gc 相互作用的機制,為菊粉酶在紅肉加工中應(yīng)用提供參考。

    1 材料與方法

    1.1 試驗材料

    來自曲霉的外切型菊粉酶標準品(EC 3.2.1.80)、內(nèi)切型菊粉酶標準品(EC 3.2.1.7),購自愛爾蘭Megazyme 公司;Neu5Gc 標準品,購自美國Sigma 公司。菊粉酶結(jié)構(gòu)從蛋白質(zhì)PDB 數(shù)據(jù)庫獲取,均來源于曲霉,外切型菊粉酶(PDB 編號:1Y9G),內(nèi)切型菊粉酶(PDB 編號: 3SC7)。Neu5Gc 分子結(jié)構(gòu)取自Pubchem 數(shù)據(jù)庫(CID: 440001)。

    1.2 主要儀器與設(shè)備

    Waters Xevo TQ 沃特世三重四級桿液相色譜質(zhì)譜聯(lián)用儀,美國Waters 公司;Mill-Q 超純水儀,美國Milipore 公司。SYBYL-X(2.0)藥物設(shè)計平臺、開源分子動力學程序Gromacs (2 016.4)、Gaussian View6、Gaussian16 (Revision A.03)[13]程序由貴州大學云計算平臺支持。Discovery studio 3.5 客戶端(蛋白去水、去雜原子和對接空腔參數(shù)記錄)、對接程序AutoDock-Vina、波函數(shù)程序Mulitwfn3.6[14]、酶底復合物二維殘基相互作用可視化程序Ligpot+、氨基酸殘基虛擬突變工具Swiss PDB Viewer[15]、等值面圖繪制程序VMD(1.9)[16],均由個人電腦提供計算平臺。

    1.3 試驗方法

    1.3.1 菊粉酶作用效果檢測 取1 mL 10 mg·mL-1的標準Neu5Gc 溶液于液相瓶中,加入0.8%(w/w)菊粉酶,置于最適溫度50℃酶解30 min,酶解完成后用0.22 μm 膜過濾酶解液,備用。采用液相色譜質(zhì)譜聯(lián)用儀進行菊粉酶作用效果檢測,根據(jù)文獻[17]的方法略有修改。取1 μL 酶解液進樣分析,色譜條件:ACQUITY UPLC C18 色譜柱(1.7 μm,2.1 mm×100 mm),流動相A:100%乙腈,流動相B:0.1% 甲酸水溶液,梯度洗脫(5% A 0~3 min,10% A 3 ~5 min,100%A 5~6 min,5% A 6~7 min),流速0.2 mL·min-1,柱溫25℃,進樣量1 μL。質(zhì)譜條件: 離子源為電噴霧化離子源(ESI),源溫度100℃,探針溫度450℃,毛細管電壓3 000 V,錐孔電壓19 V,脫溶劑氣流速1 000 L·hr-1,負離子方式檢測,掃描方式為多重反應(yīng)監(jiān)測(multiple reaction monitoring,MRM)。定量分析離子對323.8/116.0 時,碰撞氣能量為17 V;定量分析離子對323.8/185.8 時,碰撞氣能量為15 V。

    1.3.2 Neu5Gc 與菊粉酶相互作用 參照文獻[18]的方法,采用分子對接、氨基酸殘基虛擬突變、分子動力學模擬研究Neu5Gc 與菊粉酶的分子間相互作用。為更加全面的考察Neu5Gc 與菊粉酶相互作用過程中涉及的蛋白殘基,同時采用AutoDock-Vina 和SYBYL(2.0)對接方法進行分子動力學模擬,研究實際水環(huán)境下的對接復合物穩(wěn)定性。

    AutoDock-Vina 對接: 對受體菊粉酶和配體Neu5Gc 均以AutodockTool-1.5.6 打開,進行定義原子類型和添加Casteiger 電荷處理,并保存為pdbqt 格式。從PDB 蛋白數(shù)據(jù)庫得到的菊粉酶結(jié)構(gòu)含小分子配體,這些小分子配體構(gòu)成了對接時的活性空腔,為獲得Neu5Gc 與菊粉酶不同配體空腔對接的最佳位置,根據(jù)Discovery studio 3.5 客戶端中顯示的配體空間參數(shù),采用AutoDock-Vina 進行不同配體位置的對接。設(shè)置對接時的最大打分構(gòu)象數(shù)目為9,對接盒子大小設(shè)置根據(jù)空間參數(shù)適當調(diào)整。對接完成后根據(jù)結(jié)合能大小和均方位移(root mean square deviation,RMSD)值作為最佳結(jié)合構(gòu)象選擇依據(jù)。對接后的最佳復合物結(jié)構(gòu)用Ligpot+進行氨基酸二維相互作用分析。

    Surflex-Dock (SFXC)-Geom 對接:受體菊粉酶加氫后刪除配體和水分子以及雜原子,加Casteiger 電荷,在Tripos 力場下采用最陡下降法進行1 000 步的能量最小化處理,其余參數(shù)均為默認。對接空腔由配體坐標產(chǎn)生,在其周圍氨基酸0.5 nm 范圍半徑內(nèi)產(chǎn)生活性空腔。最終對接配體可能的構(gòu)象數(shù)目設(shè)定為20種,配體的RMSD 最小控制在0.5。取對接后打分最高者進行對接分析。

    Gromacs 動力學模擬:依據(jù)結(jié)合能與打分排序?qū)σ陨蠈雍蟮膹秃衔镞M行篩選,對最佳Neu5Gc-菊粉酶復合物在水環(huán)境中進行動力學模擬以考察其結(jié)合穩(wěn)定性。受體分子菊粉酶的力場選擇AMBER99-ILDN,水盒子采用TIP3P 三點水模型,體系大小為盒子邊界距離蛋白文件2.5 nm。配體分子Neu5Gc 采用Amber Tools 17 中的GAFF 力場,產(chǎn)生itp 限制勢文件和gro拓撲文件,使用Gromacs 中的genrestr 工具生成體系距離限制勢文件,使用-DPOSRES 關(guān)鍵詞限制配體的位置。對體系加入抗衡離子平衡電荷。在模擬參數(shù)文件mdp 中先用最陡下降法進行2 000 步,步長為0.01 ps的能量最小化,鄰居列表生成選擇cutoff-scheme =Verlet,靜電作用為庫倫力,非鍵相互作用采用PME 截斷方式,LINCS 算法約束成鍵,每50 步輸出一次能量信息。然后進行500 ps 的控溫、控壓限制性動力學直到平衡溫度298.15 K,熱浴設(shè)定為velocity-rescale,控溫的時間常數(shù)為0.2 ps,壓浴設(shè)定為Berendsen,采用各項同性控壓,控壓時間常數(shù)為0.5 ps,只約束氫鍵,該過程每500 步輸入一次能量信息,每1 000 步輸出一次坐標信息。限制性動力學平衡后,對配體和受體設(shè)置索引和分組,放開限制進行100 ns 的成品模擬,設(shè)定每步0.002 ps,每1 000 步輸出一次能量信息。成品模擬平衡后,選擇平衡末端80~100 ns 處軌跡文件,每隔10幀從軌跡文件中提取平均構(gòu)象,采用g_mmpbsa 工具,基于分子力學泊松-波爾茲曼表面積(molecular mechanics poisson-boltzmann surface area,MM-PBSA)法進行蛋白和配體結(jié)合自由能△Gbind計算?!鱃bind值越低說明受體和配體間的親和力越高。溶劑中蛋白和配體結(jié)合自由能按照公式進行計算[19]:

    式中,Gcomplex是總的自由能;Gprotein和Gligand是溶劑中獨立的自由能。其中,溶劑中獨立的自由能GX按照公式進行計算:

    式中,X代表蛋白或配體或復合物;TS 是真空中構(gòu)象變化導致的熵變對自由能的貢獻,由于構(gòu)象熵對本研究變化影響不大,故不考慮;EMM是真空中的分子力學能;Gsolvation代表溶劑化能;按照公式進行計算:

    式中,Ebond代表成鍵作用,如鍵長、鍵角、旋轉(zhuǎn)等;Enobonded代表非鍵作用力,如范德華力和靜電作用力;Gpolar代表極性溶劑化能,由求解波茲曼方程得到;Gnonpolar由基于溶劑可及表面積SASA 模型、溶劑可及體積SAV 和WCA 模型計算。

    計算完成后,對動力學過程中重要殘基進行能量分解,以獲得殘基對結(jié)合自由能的貢獻。為考察模擬過程的穩(wěn)定性,進行蛋白和復合物的RMSD、殘基的RMSF 分析。

    1.3.3 殘基虛擬突變 當外切菊粉酶1Y9G 的活性殘基ASP189 用ASN 取代時,酶活性降低[9];以外切菊粉酶1Y9G 為模板同源建模構(gòu)建內(nèi)切型菊粉酶,當活性殘基ASP460 突變?yōu)锳LA 后,突變體催化活性喪失[20]。內(nèi)切型菊粉酶3SC7 定點突變表明ASP42 對催化活性有重要作用[21]。采用Swiss PDB Viewer 工具,選擇將外切型菊粉酶1Y9G 的活性殘基ASP41 突變?yōu)锳SN41;將內(nèi)切型菊粉酶3SC7 的活性殘基GLU233 突變?yōu)锳LA233。對突變后的菊粉酶與Neu5Gc 以AutoDock-Vina 對接并產(chǎn)生復合物結(jié)構(gòu),進行分子動力學模擬考察穩(wěn)定性和活性位點殘基變化。提取菊粉酶1Y9G 和3SC7 與Neu5Gc 動力學模擬平衡后最后一幀,在M062X/6-311G(d,p)水平下對活性殘基-Neu5Gc 復合物以分子中原子理論(atom in molecules,AIM) 進行鍵臨界關(guān)鍵點(bond critical point,BCP) 分析,用獨立梯度模型(independent gradient model,IGM)考察分子間弱相互作用。

    1.4 數(shù)據(jù)處理與分析

    試驗數(shù)據(jù)均用平均值表示,平行3 次。采用SPSS 23. 0 軟件進行統(tǒng)計分析,數(shù)據(jù)顯著性差異分析采用t檢驗,顯著性水平為P<0.05,采用OriginPro 2016 軟件制圖。

    2 結(jié)果與分析

    2.1 菊粉酶對Neu5Gc 的作用效果

    通過液相質(zhì)譜(liquid chromatography-mass spectrometry,LC-MS) 對2 種酶在最適溫度下與Neu5Gc 標準品作用后濾膜液進行Neu5Gc 含量檢測。由圖1可知,外切型菊粉酶對Neu5Gc 標準品中Neu5Gc 含量的平均降低率為55.33%,內(nèi)切型菊粉酶對Neu5Gc 標準品中Neu5Gc 含量的平均降低率為43.81%,兩組間差異不顯著(P>0.05),但均顯著高于對照組(1.73%)。表明外切型菊粉酶對吡喃糖衍生物Neu5Gc 的降低效果略優(yōu)于內(nèi)切菊粉酶。這可能與菊粉酶對吡喃糖衍生物,如蔗糖的水解能力差異有關(guān);外切型菊粉酶能將蔗糖轉(zhuǎn)化為單糖中的果糖;而大多數(shù)內(nèi)切酶轉(zhuǎn)化蔗糖的效率較低[22-23]。

    2.2 Neu5Gc 與菊粉酶不同位點分子對接

    由于2 種菊粉酶對Neu5Gc 作用效果相近,為探究作用過程中二者間活性殘基的差異,首先采用分子對接對二者進行初步探究。用Discovery studio 3.5 客戶端顯示菊粉酶結(jié)構(gòu),外切型菊粉酶1Y9G 共5 個配體,配體1 由NAG1001、NAG1002 組成,配體2 由NAG2001 組成,配體3 由NAG3001 組成,配體4 由NAG5001 組成,配體5 由FRU801 組成。3SC7 內(nèi)切型菊粉酶X 鏈共4 個配體,配體1 由MAN517 組成,配體2 由MAN518 和MAN519 組成,配體3 由NAG750、NAG751 和BMA752 組成,配體4 由EPE520 組成。對菊粉酶不同配體位置與Neu5Gc 進行AutoDock-Vina對接,結(jié)果如表1、2 所示,對接最佳復合物處Neu5Gc-菊粉酶二維活性殘基作用結(jié)果如圖2所示。Surflex-Dock 對接的結(jié)果如圖3所示。

    圖1 菊粉酶對Neu5Gc 標準品的作用效果Fig.1 Effect of inulinase on Neu5Gc standards

    由表1可知,1Y9G 在配體5 處對接結(jié)合能最低,為-30.96 kJ·mol-1,因此Neu5Gc 與1Y9G 在配體5 處對接結(jié)果最佳。由表2可知,3SC7 在配體2 處對接結(jié)合能最低,為-26.77 kJ·mol-1,因此Neu5Gc 與3SC7在其配體2 處對接結(jié)果最佳。由圖2可知,Neu5Gc 與1Y9G 配體5處的殘基ASN40、GLN57、SER103、GLN149、ASN186、ASP189、ASN265 和TRP335 相互作用產(chǎn)生10 個氫鍵,與殘基ASP41、PHE102、ARG188 和GLU241 產(chǎn)生疏水相互作用。由圖3中SYBYL 對接結(jié)果可知,Neu5Gc 與殘基 GLN14、ASN40、GLN57、SER103、GLU241、TRP335 產(chǎn)生8 個氫鍵。外切型菊粉酶1Y9G 催化果糖的活性位點殘基為ASN40、ASP41、TRP65、SER103、ARG188、ASP189、GLU241 和CYS242[10]。對比1Y9G 活性位點可知,Neu5Gc 均可與外切型菊粉酶殘基CYS242 外的其余活性殘基發(fā)生相互作用,且其主要與ASP 型、GLN 型殘基相互作用。

    Neu5Gc 與3SC7 配體2 周圍殘基ASN42、GLU43、GLU233、ASN320、GLY323 產(chǎn)生7 個氫鍵。與殘基VAL23、TRP40、PHE99、THR255、ASN265、ASP298 產(chǎn)生疏水相互作用。SYBYL 對接結(jié)果顯示,與3SC7 的殘基GLU258、ASN265、ARG295、ASN320、GLY323、ASP440 產(chǎn)生8 個氫鍵。內(nèi)切型菊粉酶催化果糖和戊酮糖時共有的殘基為TRP40、GLU43、GLN59、ASN61、TRP67、GLY68、ILE70、PHE99、ARG175、GLU233、ASN320 和GLY323[24]。對比3SC7 的活性位點可知,Neu5Gc 可與內(nèi)切型菊粉酶活性殘基TRP40、GLU43、PHE99、GLU233、ASN320 和GLY323 相互作用,表明其主要與GLU 型殘基作用。已知在外切和內(nèi)切型菊粉酶活性殘基中,ASP 和GLU 分別是重要的催化殘基[25]。綜合以上對接后相互作用殘基分析可知,Neu5Gc 主要與外切型菊粉酶活性殘基中的ASP 殘基和內(nèi)切型菊粉酶活性殘基中的GLU 殘基相互作用。這與Singh 等[12]采用分子對接得到外切型菊粉酶和果糖- 6 -磷酸主要結(jié)合殘基為ASP22、ASP128 和ASP179,內(nèi)切型菊粉酶與酮糖主要結(jié)合殘基為GLU34和GLU200 一致;也和動力學模擬果糖-6-磷酸和酮糖分別與外切型菊粉酶ASP285、GLN350 產(chǎn)生氫鍵作用的結(jié)論一致[26]。

    表1 1Y9G 與Neu5Gc AutoDock-Vina 對接結(jié)果Table1 1Y9G and Neu5Gc AutoDock-Vina docking results

    表2 3SC7 與Neu5Gc AutoDock-Vina 對接結(jié)果Table2 3SC7 and Neu5Gc AutoDock-Vina docking results

    2.3 菊粉酶與Neu5Gc 分子動力學模擬

    2.3.1 動力學模擬穩(wěn)定性評價 為獲得2 種菊粉酶與Neu5Gc 對接后最佳酶底復合物在水環(huán)境中的穩(wěn)定結(jié)合狀態(tài),提取最佳打分復合物進行100 ns 動力學模擬。為驗證動力學模擬中菊粉酶催化殘基對結(jié)合Neu5Gc 并穩(wěn)定其構(gòu)象的貢獻,對菊粉酶催化活性殘基進行虛擬突變。采用蛋白骨架波動、配體與受體蛋白的偏移程度、蛋白殘基波動進行模擬穩(wěn)定性評價,結(jié)果分別如圖4、5、6 所示。波動數(shù)值越大表明穩(wěn)定性越低。由圖4可知,1Y9G 突變后殘基蛋白穩(wěn)定性輕微提高,但突變后復合物穩(wěn)定性小于突變之前。由圖5可知,突變后蛋白殘基運動波動性相似,但均大于晶體結(jié)構(gòu)的殘基波動,這與蛋白穩(wěn)定性是一致的。對3SC7突變活性殘基后,蛋白穩(wěn)定性變小,殘基波動明顯高于晶體試驗值,復合物穩(wěn)定性明顯減弱。

    圖2 Neu5Gc 與菊粉酶1Y9G 和3SC7 AutoDock-Vina 對接結(jié)果二維模式圖Fig.2 Two-dimensional interaction model of Neu5Gc and inulinase 1Y9G and 3SC7 by AutoDock-Vina docking

    圖3 Neu5Gc 與菊粉酶1Y9G 和3SC7 Surflex-Dock 對接結(jié)果三維模式圖Fig.3 Three-dimensional interaction model of Neu5Gc and inulinase 1Y9G and 3SC7 by Surflex-Dock docking

    2.3.2 動力學模擬結(jié)合模式 為獲得穩(wěn)定對接狀態(tài)下Neu5Gc 與菊粉酶活性殘基相互作用過程中二者間的空間結(jié)合模式,對模擬平衡后20 ns 軌跡每隔10 ps進行提取,統(tǒng)計模擬過程中活性殘基與Neu5Gc 作用產(chǎn)生的氫鍵個數(shù)和貢獻權(quán)重,選取模擬平衡后的構(gòu)象進行結(jié)合模式分析。由圖7可知,Neu5Gc 與1Y9G 的殘基ASP41、GLN57、GLU241、ASN265 和TRP335 產(chǎn)生氫鍵。當突變1Y9G 的殘基ASP41 為ASN41 后,產(chǎn)生氫鍵的殘基為ASN40、GLN57、TRP65 和GLN149,其中活性殘基僅ASN40、TRP65 與Neu5Gc 穩(wěn)定相互作用,這是由于突變后ASN41 不能親核攻擊底物的異構(gòu)碳原子形成共價中間體,使得其活性喪失。根據(jù)外切型菊粉酶活性位點可知,在水溶液環(huán)境中Neu5Gc 可與外切型菊粉酶關(guān)鍵活性殘基ASP41、GLU241 以及GLN57 穩(wěn)定相互作用。由于殘基ASP41 和GLU241分別在菊粉酶催化底物過程中充當親核試劑、催化酸堿對[27],表明外切型菊粉酶可對Neu5Gc 時行催化作用。

    圖4 菊粉酶1Y9G(A)和1Y9G-Neu5Gc 復合物(B)動力學模擬蛋白骨架均方根位移Fig.4 Dynamics simulation RMSD of inulinase 1Y9G (A) and 1Y9G-Neu5Gc complex(B) protein skeleton

    圖5 菊粉酶1Y9G (A) 和3SC7 (B)動力學模擬殘基均方根漲落Fig.5 Dynamics simulation RMSF of inulinase 1Y9G(A) and 3SC7 (B) protein residues

    Neu5Gc 與3SC7 產(chǎn)生氫鍵的殘基為ASN42、GLU43、ASN265、ASP298 和ASN320,其中催化活性殘基ASP298 和GLU43 與晶體結(jié)構(gòu)分析結(jié)果一致[28]。由于Neu5Gc 是N-乙酰葡萄糖的衍生物,且有學者以3SC7 為模板和N -乙酰葡糖胺對接,表明殘基ASP285、PHE352、THR288、VAL284、LEU343、MET289、GLN350、LEU343、HIS351、PRO283、THR271可與其穩(wěn)定相互作用[29]。對比可知,對接獲得的殘基相互作用信息與水環(huán)境中動力學模擬得到的作用殘基結(jié)果差異較大,表明與Neu5Gc 產(chǎn)生較穩(wěn)定氫鍵作用的殘基對結(jié)合過程貢獻較大。當突變GLU233 為ALA233 后,Neu5Gc 僅與殘基ASP176、ARG175 和GLU173 產(chǎn)生氫鍵作用。根據(jù)內(nèi)切型菊粉酶活性位點可知,Neu5Gc 可與內(nèi)切型菊粉酶關(guān)鍵殘基GLU43、ASN320 穩(wěn)定相互作用。突變內(nèi)切型菊粉酶3SC7 的殘基GLU233 為ALA233 時,Neu5Gc 僅與3SC7 的活性殘基ARG175 產(chǎn)生氫鍵作用,這是由于GLU 較大的側(cè)鏈能夠造成底物結(jié)合過程的崩潰,由此影響了活性空腔的酸性環(huán)境。

    2.3.3 結(jié)合過程殘基能量分解 為獲得水環(huán)境下Neu5Gc 與菊粉酶結(jié)合過程中分子間相互作用類型、活性殘基穩(wěn)定Neu5Gc 構(gòu)象到活性空腔中的能量貢獻,對動力學模擬平衡后20 ns 軌跡每隔100 ps 進行提取,采用g_mmpbsa 工具對酶和底物結(jié)合能貢獻進行波爾茲曼能量分解,結(jié)果如表3所示。并對模擬平衡時菊粉酶與底物Neu5Gc 相互作用貢獻較大的殘基進行能量貢獻分解,結(jié)果如圖8所示。

    圖6 Neu5Gc 相對于菊粉酶3SC7(A) 和3SC7-Neu5Gc 復合物動力學模擬均方位移Fig.6 Dynamics simulation RMSD of Neu5Gc on inulinase 3SC7(A) and 3SC7-Neu5Gc(B)

    圖7 Neu5Gc 與菊粉酶1Y9G(A)、3SC7(B)、1Y9G 突變后(C)和3SC7 突變后(D)動力學模擬平衡結(jié)合模式圖Fig.7 Binding model of Neu5Gc and inulinase 1Y9G (A),3SC7 (B),mutated 1Y9G (C) and mutated 3SC7 (D)under dynamics simulation equilibrium

    表3 Neu5Gc 與菊粉酶體系動力學平衡后結(jié)合能貢獻分解Table3 Binding energy contribution decomposition of Neu5Gc and inulinase system under dynamics simulation equilibrium /(kJ·mol-1)

    由表3可知,對Neu5Gc 與菊粉酶結(jié)合過程起阻礙作用主要是極性溶劑化能,起促進作用的是范德華力、靜電相互作用和溶劑可及表面積。未突變時,外切菊粉酶1Y9G 與Neu5Gc 的結(jié)合能明顯小于內(nèi)切型菊粉酶3SC7,表明Neu5Gc 優(yōu)先與外切型菊粉酶相互作用。當突變1Y9G 的催化殘基后,范德華力增強,這可能是由于天冬酰胺使得結(jié)合位點的電荷由原來的負電荷變?yōu)檎姾?,靜電相互作用減少。這與Torabizadeh等[30]通過以1Y9M 建模的菊粉酶活性殘基經(jīng)磷酸吡哆醛修飾后,負電荷增加,靜電作用增強的結(jié)果相似。突變3SC7 的GLU233 為ALA233 后,總結(jié)合能增強,但靜電相互作用對結(jié)合過程的促進作用明顯減少,導致復合物穩(wěn)定性明顯降低。由圖8可知,結(jié)合過程中主要相關(guān)殘基的能量項主要由對結(jié)合過程起促進作用的分子力學能和起阻礙作用的極性溶劑化能組成。從有利于結(jié)合過程的分子力學能變化可知,1Y9G 中促進與Neu5Gc 結(jié)合的關(guān)鍵殘基是ASP41、GLU241 和TRP335,而當突變殘基ASP41 后主要為TRP65 殘基。3SC7 中促進與Neu5Gc 結(jié)合的關(guān)鍵殘基是GLU43、ASN265 和ASP298,當突變GLU233 后主要為ASP176殘基。表明突變關(guān)鍵活性殘基后,菊粉酶對于Neu5Gc穩(wěn)定到活性空腔的靜電能將下降,進一步證明了Neu5Gc 與菊粉酶主要以催化活性殘基相互穩(wěn)定作用。

    2.3.4 AIM 和IGM 分析 為考察動力學平衡后Neu5Gc 與活性殘基形成的氫鍵強度和分布,提取1Y9G 和3SC7 與Neu5Gc 動力學模擬平衡后關(guān)鍵結(jié)合殘基,選擇對弱相互作用計算較好的M062X 泛函,在6-311G(d,p)基組水平下產(chǎn)生波函數(shù)文件對臨界鍵關(guān)鍵點BCP 進行AIM 分析,結(jié)果如表4所示。采用獨立梯度約化泛函IGM 對氫鍵投影進行等值面分析,結(jié)果如圖9所示;分子間和分子內(nèi)相互作用散點圖如圖10所示。

    圖8 Neu5Gc 與菊粉酶1Y9G(A)、3SC7(B)、1Y9G 突變后(C)和3SC7 突變后(D)動力學平衡時殘基能量分解Fig.8 Residual energy decomposition of Neu5Gc and inulinase 1Y9G (A),3SC7 (B),mutated 1Y9G (C)and mutated 3SC7 (D) under dynamics simulation equilibrium

    表4 Neu5Gc 與菊粉酶1Y9G、3SC7 活性殘基分子間氫鍵AIM 分析(能量項ρ、▽2ρ、V、G、HBCP、EHB)Table4 AIM analysis of intermolecular hydrogen bonding between Neu5Gc and inulinase 1Y9G,3SC7 active residues (energy items ρ、▽2ρ、V、G、HBCP、EH B)

    在AIM 理論中,體系坐標函數(shù)的二階鞍點稱為鍵臨界點。鍵臨界點處的電子密度、拉普拉斯密度與鍵的強度和類型有關(guān)[31]。從鍵臨界點搜索結(jié)果可以看出,1Y9G 和3SC7 活性殘基與Neu5Gc 均產(chǎn)生4 個氫鍵臨界點。由Rozas 等[32]的氫鍵強度劃分標準:弱氫鍵(EHB<50.20 kJ·mol-1,▽2(rBCP) > 0,G(rBCP) +V(rBCP) >0),中等氫鍵(50.20 <EHB<100.4 kJ·mol-1,▽2(rBCP) > 0,G(rBCP) +V(rBCP) <0),強氫鍵(EHB>100.4 kJ·mol-1,▽2(rBCP) <0,G(rBCP) +V(rBCP) <0)。由表4可知,外切型菊粉酶1Y9G 活性殘基與Neu5Gc 作用的最大氫鍵強度為-9.65 kJ·mol-1,內(nèi)切型菊粉酶3SC7 活性殘基與Neu5Gc 作用的最大氫鍵強度為-11.53 kJ·mol-1,可知2 種菊粉酶與Neu5Gc均以弱氫鍵相互穩(wěn)定結(jié)合,其中,就氫鍵強度總和而言,外切型菊粉酶稍高于內(nèi)切型菊粉酶的作用效果,這與結(jié)合自由能分解得到Neu5Gc 更易于和內(nèi)切型菊粉酶在水環(huán)境下穩(wěn)定結(jié)合的結(jié)果是一致的。

    約化密度梯度(reduced density gradient,RDG)函數(shù)是基于電子密度拓撲的非共價相互作用分析方法,電子密度梯度降低的區(qū)域與非共價相互作用區(qū)域相關(guān),該區(qū)域用等值面圖可以將非共價作用的位置分布和類型進行可視化[33],以S 作為弱相互作用的等值面函數(shù)。

    圖9 Neu5Gc 與菊粉酶1Y9G(A)、3SC7(B)活性殘基分子間作用IGM 等值面圖Fig.9 Intramolecular IGM isosurface map of Neu5Gc with inulinase 1Y9G (A) and 3SC7 (B) active residues

    圖10 Neu5Gc 與菊粉酶1Y9G、3SC7 活性殘基分子間作用IGM 散點圖Fig.10 Intramolecular IGM scatter plot of Neu5Gc with inulinase 1Y9G and 3SC7 active residues

    式中,ρ 代表電子密度,▽ρ 代表電子密度梯度的模量。在IGM 理論中,當SIGM基于|▽ρ(r)IGM|時,SIGM和S 最大的差值代表鍵臨界點處的RDG 值,定義δg表示相互作用區(qū)域,劃分為片段內(nèi)相互作用的δgintra和片段間相互作用的δginter。其表達式如下:δg(r)=|▽ρ(r)IGM|-|▽ρ(r)|,δg(r)=δg(r)inter+δg(r)intra。將δg 描述符以符號函數(shù)sign(λ2)ρ 投影。當sign(λ2)ρ值為-0.04 時,δginter值顯示明顯的峰,代表了分子間相互作用的氫鍵;當sign(λ2)ρ 值為正值時,代表了空間排斥作用;當sign(λ2)ρ 值在-0.3 附近,δgintra有個較大峰,代表了分子內(nèi)相互作用。默認的sign(λ2)ρ 值在-0.05~0.05 范圍內(nèi)。等值面圖中藍色越強表示吸引作用越強,綠色的代表弱的范德華力作用[34]。

    由圖9可知,與Neu5Gc 穩(wěn)定結(jié)合的菊粉酶活性殘基和Neu5Gc 產(chǎn)生氫鍵的氮氧原子與受體之間有明顯的等值面分布。從等值面分布上看,1Y9G 活性殘基與Neu5Gc 產(chǎn)生的藍色區(qū)域明顯少于3SC7 的活性殘基,表明分子間靜電作用力弱于3SC7,這與結(jié)合能分解得到其靜電作用能較低相一致。由顏色分布可知范德華力是主要作用力,氮原子為供體產(chǎn)生的較強氫鍵僅為1Y9G 的GLN57 處,這與AIM 分析的結(jié)論一致。由3SC7 活性殘基和Neu5Gc 弱相互作用的等值面分布圖可知,產(chǎn)生氫鍵的位置主要為供體羥基氧原子和受體氧之間,與AIM 分析結(jié)果相一致。同樣,散點圖10中觀察到菊粉酶與Neu5Gc 分子間存在的弱氫鍵相互作用峰(sign(λ2)ρ -0.04 處)。以上結(jié)果表明Neu5Gc 羥基氧和酰胺氮與菊粉酶活性殘基提供的親核性質(zhì)是導致該結(jié)合過程穩(wěn)定存在的原因之一。

    3 討論

    菊粉酶廣泛用于功能性低聚糖的制備,但由于菊粉酶催化活性域具有一定保守性,不同底物與不同類型菊粉酶之間相互作用的過程和機制有一定共性。菊粉酶可以作用于單糖衍生物[12],對酮糖也有一定的降低效果[26],推測菊粉酶具有消減吡喃糖衍生物Neu5Gc 的潛在功能。因此,基于前期對多種食品酶制劑作用于Neu5Gc 的篩選結(jié)果進行深入研究,采用液相質(zhì)譜檢測以克服Neu5Gc 液相衍生化檢測干擾,其中,外切菊粉酶結(jié)果(55.33%)較液相法檢測結(jié)果(50.52%)提高了4.81%[7],這是由于衍生化試劑1,2-二氨基-4,5-亞甲二氧基苯(DMB) 會與菊粉酶活性中心競爭性結(jié)合,干擾檢測。外切型和內(nèi)切型菊粉酶活性域的差異性導致其對相同底物的作用效果不同,因此本研究觀察到內(nèi)切型菊粉酶對Neu5Gc 含量的降低率較外切型菊粉酶低6.71%左右。Neu5Gc 雖然不是菊粉酶的高效底物,但高達50%以上的降低率表明菊粉酶對Neu5Gc 具有一定的消減能力。分子動力學模擬廣泛用于酶底物結(jié)合過程細節(jié)的討論,故采用該方法研究2 種菊粉酶作用于Neu5Gc 的效果和差異機制。

    目前對菊粉酶有關(guān)的分子動力學模擬研究多集中在理性設(shè)計和底物識別方面[35],但對其于分子間弱相互作用的討論不夠深入。國內(nèi)武希箭[36]采用分子動力學模擬研究內(nèi)切型菊粉酶對菊粉的分子間相互作用,考察了結(jié)合能、殘基突變情況,發(fā)現(xiàn)主導菊粉酶發(fā)揮作用的關(guān)鍵分子間作用力是靜電力和范德華力,且活性位點殘基側(cè)鏈減少后,對底物穩(wěn)定到活性口袋的作用減弱,這與本研究結(jié)論一致,但其未定量研究活性殘基與底物結(jié)合的氫鍵強度和分布情況[36]。目前對于酶和底物分子對接的研究多采用單一平臺進行分析,形成復合物后采用動力學模擬進行驗證,但由于不同分子對接軟件的打分函數(shù)和評價規(guī)則不同,導致一些潛在的酶活性殘基不能被識別[37]。本研究綜合利用AutoDock-Vina 的高并行性和SYBYL 平臺的高精度分子力場優(yōu)勢進行對接,殘基與底物分子互作結(jié)果也表明了不同對接平臺在識別活性殘基方面的差異,采用多平臺對接可獲取更多酶底物結(jié)合過程信息。本研究在水環(huán)境下對于Neu5Gc 和菊粉酶對接復合物的動力學模擬表明,Neu5Gc 能夠與菊粉酶活性殘基以靜電作用力、范德華力和疏水作用力穩(wěn)定結(jié)合,氨基酸虛擬突變也證實了菊粉酶催化活性殘基ASP 和GLU 的重要貢獻,這與大多數(shù)研究結(jié)論一致[27]。但本研究通過分子中原子理論對活性殘基和Neu5Gc 復合物的分子間氫鍵定量分析發(fā)現(xiàn),外切型菊粉酶1Y9G 與Neu5Gc的分子間氫鍵強度較內(nèi)切型菊粉酶3SC7 高18.7 kJ·mol-1左右,這可能是二者對Neu5Gc 含量降低效果不同的原因之一。氫鍵強度差異對于酶催化過程中過渡態(tài)的形成和反應(yīng)勢壘具有一定影響,隨著量子力學/分子力學(quantum mechanics/molecular mechanics,QM/MM) 法在研究糖苷酶水解、糖基轉(zhuǎn)移酶的反應(yīng)機制方面的應(yīng)用[38],活性殘基對底物作用過程中熱力學和動力學變化逐漸被解析,但該過程涉及復雜力場參數(shù)的新建和優(yōu)化,較為困難。對于開發(fā)高效消減Neu5Gc 及其他食源性風險物質(zhì)的食品加工酶類,在獲知酶底物互作信息基礎(chǔ)上,利用分子生物學技術(shù)即可為實現(xiàn)酶的理性設(shè)計提供參考和支持。

    4 結(jié)論

    體外菊粉酶降低Neu5Gc 標準品含量的試驗結(jié)果表明,外切型菊粉酶和內(nèi)切型菊粉酶分別對Neu5Gc降低效果最佳為55.33%和43.81%。選擇兩類菊粉酶(PDB ID: 1Y9G 和3SC7)與Neu5Gc 穩(wěn)定結(jié)合的分子對接復合物進行分子動力學模擬和弱相互作用,結(jié)果表明,Neu5Gc 可以與外切型菊粉酶1Y9G 的活性殘基ASP41 和GLU241 通過分子間氫鍵產(chǎn)生穩(wěn)定的相互作用,內(nèi)切型菊粉酶3SC7 的活性殘基ASP298 和GLU43 對結(jié)合Neu5Gc 到其活性口袋過程有重要貢獻。突變1Y9G 中的關(guān)鍵活性殘基ASP41 為ASN41,3SC7 的活性殘基GLU233 為ALA233 后的結(jié)果顯示,活性殘基中ASP 的親核性和GLU 的酸堿對性質(zhì)對菊粉酶結(jié)合Neu5Gc 并導致其含量降低過程具有重要作用。殘基能量分解、AIM 分析和IGM 等值面分析表明,Neu5Gc 通過羥基氧原子和酰胺氮原子與菊粉酶活性殘基發(fā)生靜電作用主導的分子間弱相互作用,穩(wěn)定占據(jù)到菊粉酶的活性位點。紅肉是人體多種營養(yǎng)物質(zhì)的重要來源,本研究通過解析菊粉酶作用Neu5Gc 的活性殘基信息,為菊粉酶的理性設(shè)計和定向進化提供了新的理論參考。

    猜你喜歡
    外切菊粉殘基
    菊粉在凝固型發(fā)酵乳中的應(yīng)用
    基于各向異性網(wǎng)絡(luò)模型研究δ阿片受體的動力學與關(guān)鍵殘基*
    關(guān)于橢圓外切平行四邊形的一個幾何不變量
    “殘基片段和排列組合法”在書寫限制條件的同分異構(gòu)體中的應(yīng)用
    探究拋物線內(nèi)接、外切三角形的性質(zhì)
    橢圓內(nèi)接外切六邊形的幾何特性研討
    圓外切三角形與圓的關(guān)系
    降脂減肥用菊粉
    飲食科學(2016年9期)2016-11-18 15:59:36
    蛋白質(zhì)二級結(jié)構(gòu)序列與殘基種類間關(guān)聯(lián)的分析
    菊粉珍珠奶茶的研制
    精品久久蜜臀av无| 最近最新中文字幕大全免费视频 | 美国免费a级毛片| 国产一区亚洲一区在线观看| 国产精品 国内视频| 97精品久久久久久久久久精品| 可以免费在线观看a视频的电影网站 | 叶爱在线成人免费视频播放| 久久狼人影院| 国产乱来视频区| 国产成人免费无遮挡视频| 午夜老司机福利片| 狠狠婷婷综合久久久久久88av| 亚洲国产精品成人久久小说| 桃花免费在线播放| 国产成人一区二区在线| 欧美日韩亚洲高清精品| 我要看黄色一级片免费的| 男女高潮啪啪啪动态图| 国产精品一国产av| 免费高清在线观看日韩| 久久精品熟女亚洲av麻豆精品| 黄网站色视频无遮挡免费观看| 亚洲情色 制服丝袜| 亚洲美女视频黄频| 天天操日日干夜夜撸| 欧美人与善性xxx| 新久久久久国产一级毛片| 国产黄色视频一区二区在线观看| 午夜免费鲁丝| 少妇人妻久久综合中文| 亚洲欧美精品自产自拍| 午夜福利,免费看| 999久久久国产精品视频| 亚洲自偷自拍图片 自拍| 国产日韩一区二区三区精品不卡| 亚洲一卡2卡3卡4卡5卡精品中文| 51午夜福利影视在线观看| 午夜福利一区二区在线看| 日韩精品免费视频一区二区三区| 女的被弄到高潮叫床怎么办| 国产精品久久久久久精品电影小说| av国产精品久久久久影院| 又粗又硬又长又爽又黄的视频| 精品国产一区二区久久| 这个男人来自地球电影免费观看 | 免费av中文字幕在线| 亚洲国产av新网站| 亚洲成人av在线免费| 国产淫语在线视频| 久久精品亚洲av国产电影网| 十八禁网站网址无遮挡| 国产xxxxx性猛交| 在线亚洲精品国产二区图片欧美| 在线 av 中文字幕| 亚洲精品国产av成人精品| 国产视频首页在线观看| 亚洲国产中文字幕在线视频| 国产麻豆69| 成人手机av| 一级片免费观看大全| 在线天堂中文资源库| 久久毛片免费看一区二区三区| 久久久国产欧美日韩av| 丰满饥渴人妻一区二区三| av女优亚洲男人天堂| 99久久99久久久精品蜜桃| 天堂8中文在线网| 爱豆传媒免费全集在线观看| 国产熟女午夜一区二区三区| 如日韩欧美国产精品一区二区三区| 欧美精品人与动牲交sv欧美| 美女脱内裤让男人舔精品视频| 国产探花极品一区二区| 国产不卡av网站在线观看| 亚洲国产毛片av蜜桃av| 性色av一级| 久久久久久免费高清国产稀缺| 午夜福利视频精品| 久久99热这里只频精品6学生| 两个人看的免费小视频| 色婷婷av一区二区三区视频| 欧美日韩视频精品一区| 久久免费观看电影| 日韩熟女老妇一区二区性免费视频| 日韩,欧美,国产一区二区三区| 欧美人与善性xxx| 国产色婷婷99| 亚洲av在线观看美女高潮| videos熟女内射| 精品国产一区二区三区四区第35| 两个人免费观看高清视频| 91精品三级在线观看| 亚洲中文av在线| 免费高清在线观看日韩| 国产 一区精品| 亚洲熟女精品中文字幕| 亚洲国产欧美一区二区综合| 成年美女黄网站色视频大全免费| 女人久久www免费人成看片| 熟女av电影| 亚洲国产看品久久| videos熟女内射| 纯流量卡能插随身wifi吗| 在线免费观看不下载黄p国产| 免费在线观看视频国产中文字幕亚洲 | 在线观看三级黄色| 可以免费在线观看a视频的电影网站 | 免费在线观看黄色视频的| 色精品久久人妻99蜜桃| 亚洲精品aⅴ在线观看| 亚洲第一区二区三区不卡| 国产在线免费精品| 久久精品熟女亚洲av麻豆精品| 天天躁狠狠躁夜夜躁狠狠躁| 极品少妇高潮喷水抽搐| 99热全是精品| a 毛片基地| 男人爽女人下面视频在线观看| 九色亚洲精品在线播放| 侵犯人妻中文字幕一二三四区| 一二三四在线观看免费中文在| 一边摸一边抽搐一进一出视频| 热re99久久精品国产66热6| 亚洲av综合色区一区| 9热在线视频观看99| 欧美乱码精品一区二区三区| 91精品伊人久久大香线蕉| 精品少妇久久久久久888优播| 搡老乐熟女国产| 久久久精品94久久精品| 高清av免费在线| 亚洲国产欧美在线一区| 亚洲av男天堂| 精品国产一区二区久久| 看非洲黑人一级黄片| 精品国产乱码久久久久久小说| 免费人妻精品一区二区三区视频| 中文天堂在线官网| 欧美精品人与动牲交sv欧美| 亚洲欧美精品自产自拍| 狠狠婷婷综合久久久久久88av| 男的添女的下面高潮视频| 黑丝袜美女国产一区| 在线天堂最新版资源| 51午夜福利影视在线观看| 99国产综合亚洲精品| 高清黄色对白视频在线免费看| 久久免费观看电影| 制服诱惑二区| 国产一区二区在线观看av| 精品久久久久久电影网| xxx大片免费视频| 丰满乱子伦码专区| 精品第一国产精品| 国产一区二区 视频在线| 国产精品女同一区二区软件| 丰满迷人的少妇在线观看| 日本av手机在线免费观看| 在线观看免费高清a一片| 丰满少妇做爰视频| 热99国产精品久久久久久7| 老司机亚洲免费影院| 久久久久久久大尺度免费视频| 人人妻,人人澡人人爽秒播 | 最新在线观看一区二区三区 | 伦理电影大哥的女人| 熟女av电影| 人人妻人人爽人人添夜夜欢视频| 老汉色av国产亚洲站长工具| 亚洲国产日韩一区二区| 色94色欧美一区二区| 免费高清在线观看日韩| 男女之事视频高清在线观看 | 国产亚洲av片在线观看秒播厂| 亚洲成人免费av在线播放| 天天躁夜夜躁狠狠躁躁| 国产在线一区二区三区精| 91成人精品电影| avwww免费| 亚洲精华国产精华液的使用体验| 久久av网站| 欧美日本中文国产一区发布| 满18在线观看网站| 中文字幕人妻熟女乱码| 99热国产这里只有精品6| 99热网站在线观看| 国产精品人妻久久久影院| 国产精品欧美亚洲77777| 成人18禁高潮啪啪吃奶动态图| 人人澡人人妻人| 国产日韩一区二区三区精品不卡| 欧美日韩视频高清一区二区三区二| 亚洲国产精品一区三区| 人人澡人人妻人| 国产成人a∨麻豆精品| 欧美国产精品va在线观看不卡| 99久国产av精品国产电影| 最近2019中文字幕mv第一页| 亚洲成人一二三区av| av线在线观看网站| 国产色婷婷99| 国产一卡二卡三卡精品 | 欧美精品一区二区大全| 亚洲av中文av极速乱| 如日韩欧美国产精品一区二区三区| 久久精品久久久久久噜噜老黄| 免费不卡黄色视频| 国产亚洲欧美精品永久| 成人影院久久| 啦啦啦在线免费观看视频4| 咕卡用的链子| 丝袜脚勾引网站| 亚洲精品国产区一区二| 亚洲美女搞黄在线观看| 无遮挡黄片免费观看| 考比视频在线观看| 久久鲁丝午夜福利片| 狠狠婷婷综合久久久久久88av| 哪个播放器可以免费观看大片| 中文字幕最新亚洲高清| 曰老女人黄片| 纯流量卡能插随身wifi吗| 日韩视频在线欧美| 亚洲久久久国产精品| 51午夜福利影视在线观看| 街头女战士在线观看网站| 欧美日韩一级在线毛片| 在线观看三级黄色| 美女福利国产在线| 亚洲精品成人av观看孕妇| 欧美日韩亚洲高清精品| 国产伦理片在线播放av一区| 老司机影院毛片| 男的添女的下面高潮视频| 大香蕉久久成人网| 大码成人一级视频| 久久韩国三级中文字幕| 美女视频免费永久观看网站| 男女无遮挡免费网站观看| 亚洲中文av在线| 伊人亚洲综合成人网| 国产国语露脸激情在线看| svipshipincom国产片| 久久精品国产亚洲av涩爱| 国产女主播在线喷水免费视频网站| 久久99精品国语久久久| 少妇的丰满在线观看| 日本爱情动作片www.在线观看| www日本在线高清视频| 悠悠久久av| 尾随美女入室| 午夜福利网站1000一区二区三区| 亚洲国产av新网站| 观看美女的网站| 超碰成人久久| 在线亚洲精品国产二区图片欧美| 精品人妻在线不人妻| 亚洲成人av在线免费| 男人操女人黄网站| av天堂久久9| 高清欧美精品videossex| 蜜桃国产av成人99| 一区二区日韩欧美中文字幕| 午夜福利在线免费观看网站| 久久精品亚洲av国产电影网| 成人午夜精彩视频在线观看| 亚洲,一卡二卡三卡| 黑人巨大精品欧美一区二区蜜桃| 亚洲av在线观看美女高潮| 尾随美女入室| 天天操日日干夜夜撸| 亚洲av中文av极速乱| 日韩制服骚丝袜av| 七月丁香在线播放| 精品少妇久久久久久888优播| 高清在线视频一区二区三区| 国产精品 国内视频| 欧美日韩视频高清一区二区三区二| 日日啪夜夜爽| 久久精品亚洲av国产电影网| 曰老女人黄片| 午夜福利免费观看在线| 久久久久人妻精品一区果冻| 无遮挡黄片免费观看| 成人手机av| 免费日韩欧美在线观看| 捣出白浆h1v1| 国产精品三级大全| 亚洲伊人色综图| 免费人妻精品一区二区三区视频| 高清欧美精品videossex| 捣出白浆h1v1| 国产人伦9x9x在线观看| tube8黄色片| 美女主播在线视频| 秋霞在线观看毛片| 国产成人精品久久二区二区91 | 亚洲,欧美精品.| 久久久久国产一级毛片高清牌| 欧美成人精品欧美一级黄| av网站免费在线观看视频| 91精品伊人久久大香线蕉| 国产伦理片在线播放av一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产成人一精品久久久| 国产黄色免费在线视频| 超碰成人久久| videosex国产| 黑人欧美特级aaaaaa片| 午夜影院在线不卡| 69精品国产乱码久久久| 一二三四在线观看免费中文在| 丝袜美腿诱惑在线| 亚洲色图综合在线观看| 1024香蕉在线观看| 亚洲av福利一区| 老熟女久久久| 波多野结衣一区麻豆| 久久久久国产精品人妻一区二区| 国产伦理片在线播放av一区| √禁漫天堂资源中文www| 无遮挡黄片免费观看| 日韩不卡一区二区三区视频在线| 国产精品一二三区在线看| 人成视频在线观看免费观看| 亚洲国产最新在线播放| 久久久久精品人妻al黑| 成人影院久久| 久久狼人影院| 十分钟在线观看高清视频www| 在线亚洲精品国产二区图片欧美| av在线观看视频网站免费| 下体分泌物呈黄色| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦在线免费观看视频4| 亚洲,欧美精品.| 久久久久精品国产欧美久久久 | 五月天丁香电影| 黑丝袜美女国产一区| 一区二区三区四区激情视频| 熟女av电影| 国产亚洲午夜精品一区二区久久| 国产精品麻豆人妻色哟哟久久| 亚洲美女视频黄频| 另类精品久久| 国产成人精品无人区| av电影中文网址| 亚洲成色77777| 亚洲 欧美一区二区三区| 黄片小视频在线播放| 欧美另类一区| 日韩精品免费视频一区二区三区| 大码成人一级视频| 涩涩av久久男人的天堂| 国产麻豆69| 中文字幕亚洲精品专区| 国产精品无大码| 麻豆精品久久久久久蜜桃| 日韩人妻精品一区2区三区| 五月天丁香电影| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 天天添夜夜摸| 午夜福利视频精品| 国产精品.久久久| 国产99久久九九免费精品| 久久av网站| 国产视频首页在线观看| 午夜福利乱码中文字幕| 啦啦啦在线观看免费高清www| 久久久国产欧美日韩av| 国产成人91sexporn| 国产深夜福利视频在线观看| 欧美日韩福利视频一区二区| 老司机靠b影院| 国产一区二区三区av在线| 建设人人有责人人尽责人人享有的| 日韩av在线免费看完整版不卡| 亚洲欧美中文字幕日韩二区| 亚洲国产精品一区二区三区在线| 五月天丁香电影| 国产成人精品无人区| 亚洲一区中文字幕在线| 大话2 男鬼变身卡| 2021少妇久久久久久久久久久| 国产av码专区亚洲av| 久久精品国产综合久久久| 亚洲成色77777| 国产免费一区二区三区四区乱码| 国产高清不卡午夜福利| 亚洲国产欧美网| 人人妻人人澡人人爽人人夜夜| 国产成人免费无遮挡视频| 欧美激情极品国产一区二区三区| 亚洲色图综合在线观看| 女人久久www免费人成看片| 啦啦啦在线免费观看视频4| 伦理电影大哥的女人| 91精品国产国语对白视频| 精品酒店卫生间| 卡戴珊不雅视频在线播放| 美女高潮到喷水免费观看| 国产成人免费无遮挡视频| 久久女婷五月综合色啪小说| 王馨瑶露胸无遮挡在线观看| 久久女婷五月综合色啪小说| 精品久久久精品久久久| 国产有黄有色有爽视频| 一本大道久久a久久精品| 免费女性裸体啪啪无遮挡网站| 又黄又粗又硬又大视频| 欧美xxⅹ黑人| 国产精品嫩草影院av在线观看| 亚洲国产精品一区二区三区在线| 下体分泌物呈黄色| 丝袜人妻中文字幕| 国产精品一二三区在线看| 成人黄色视频免费在线看| 老司机靠b影院| 成人国产麻豆网| av线在线观看网站| 日日撸夜夜添| tube8黄色片| 一边亲一边摸免费视频| 秋霞伦理黄片| 国产视频首页在线观看| 亚洲在久久综合| 高清不卡的av网站| 久久女婷五月综合色啪小说| 又黄又粗又硬又大视频| av网站免费在线观看视频| 曰老女人黄片| 成人毛片60女人毛片免费| 一级a爱视频在线免费观看| 又大又爽又粗| 亚洲成人一二三区av| 青草久久国产| 一区二区三区激情视频| 男人操女人黄网站| 久久久久久人妻| 99久久精品国产亚洲精品| 久久国产精品男人的天堂亚洲| 2018国产大陆天天弄谢| 赤兔流量卡办理| 9色porny在线观看| 国产一区二区三区综合在线观看| 人人澡人人妻人| 精品少妇黑人巨大在线播放| 热re99久久国产66热| 在线天堂最新版资源| 久久久久网色| 视频在线观看一区二区三区| 丝袜在线中文字幕| 久久久精品国产亚洲av高清涩受| 国产亚洲精品第一综合不卡| 日本猛色少妇xxxxx猛交久久| 激情视频va一区二区三区| 电影成人av| 亚洲三区欧美一区| 免费在线观看黄色视频的| 亚洲国产精品国产精品| 又粗又硬又长又爽又黄的视频| 欧美变态另类bdsm刘玥| 看免费av毛片| 色网站视频免费| 欧美国产精品va在线观看不卡| 久久久久久久久久久久大奶| 日韩中文字幕欧美一区二区 | 久久青草综合色| 大陆偷拍与自拍| 亚洲欧美中文字幕日韩二区| 制服诱惑二区| 男女午夜视频在线观看| 精品午夜福利在线看| 另类亚洲欧美激情| 天美传媒精品一区二区| 大陆偷拍与自拍| 欧美日韩精品网址| 美女主播在线视频| 亚洲第一av免费看| 制服人妻中文乱码| 性高湖久久久久久久久免费观看| 亚洲自偷自拍图片 自拍| 天天操日日干夜夜撸| 免费看av在线观看网站| 菩萨蛮人人尽说江南好唐韦庄| 搡老岳熟女国产| 亚洲精品国产色婷婷电影| 捣出白浆h1v1| 亚洲四区av| 又大又爽又粗| 亚洲国产精品国产精品| 一级毛片电影观看| 亚洲av日韩精品久久久久久密 | 亚洲熟女毛片儿| 欧美日韩福利视频一区二区| 欧美在线黄色| 国产精品女同一区二区软件| 欧美日韩亚洲综合一区二区三区_| 欧美日韩av久久| 国产熟女午夜一区二区三区| 搡老岳熟女国产| 久久久久久久久免费视频了| 精品国产乱码久久久久久男人| 精品少妇一区二区三区视频日本电影 | 久久99一区二区三区| 欧美成人午夜精品| 午夜老司机福利片| 国产淫语在线视频| 色视频在线一区二区三区| 成人亚洲欧美一区二区av| 亚洲欧美中文字幕日韩二区| 少妇被粗大猛烈的视频| avwww免费| 黄色怎么调成土黄色| 人体艺术视频欧美日本| 亚洲精品第二区| 99热国产这里只有精品6| 一级爰片在线观看| 亚洲情色 制服丝袜| 久久国产精品大桥未久av| 九色亚洲精品在线播放| 中文字幕人妻熟女乱码| 亚洲第一av免费看| 最近中文字幕高清免费大全6| 国产精品av久久久久免费| 一二三四在线观看免费中文在| 成人国产麻豆网| 一区二区日韩欧美中文字幕| 丁香六月欧美| 欧美在线一区亚洲| 晚上一个人看的免费电影| 亚洲综合色网址| 国产精品二区激情视频| 久久精品亚洲熟妇少妇任你| 91aial.com中文字幕在线观看| tube8黄色片| 青青草视频在线视频观看| 免费女性裸体啪啪无遮挡网站| 妹子高潮喷水视频| 一本—道久久a久久精品蜜桃钙片| 国产毛片在线视频| 日韩免费高清中文字幕av| 亚洲人成电影观看| av在线app专区| 黑人欧美特级aaaaaa片| 看免费成人av毛片| 亚洲婷婷狠狠爱综合网| 香蕉丝袜av| 国产黄色免费在线视频| 最近的中文字幕免费完整| 亚洲精品自拍成人| 麻豆精品久久久久久蜜桃| 69精品国产乱码久久久| 韩国精品一区二区三区| 亚洲成国产人片在线观看| 国产免费又黄又爽又色| 一级毛片电影观看| 国产精品二区激情视频| 9191精品国产免费久久| 国产伦人伦偷精品视频| 99久国产av精品国产电影| 大话2 男鬼变身卡| 亚洲欧美一区二区三区国产| 视频在线观看一区二区三区| 亚洲成人一二三区av| 日本一区二区免费在线视频| 高清在线视频一区二区三区| 成年女人毛片免费观看观看9 | 一本—道久久a久久精品蜜桃钙片| 国产色婷婷99| 国产精品 欧美亚洲| 熟女av电影| 一二三四中文在线观看免费高清| 午夜福利免费观看在线| 欧美 亚洲 国产 日韩一| 美女大奶头黄色视频| 日韩一区二区三区影片| 99精品久久久久人妻精品| 精品亚洲乱码少妇综合久久| 午夜激情久久久久久久| 亚洲成人国产一区在线观看 | 亚洲熟女毛片儿| av又黄又爽大尺度在线免费看| 国产成人精品福利久久| 欧美日韩精品网址| 国产成人免费观看mmmm| 人妻一区二区av| 国产精品熟女久久久久浪| 日韩精品有码人妻一区| 久久毛片免费看一区二区三区| 亚洲欧美一区二区三区久久| 男人操女人黄网站| 国产av一区二区精品久久| 99九九在线精品视频| 亚洲精品,欧美精品| 青春草视频在线免费观看| 王馨瑶露胸无遮挡在线观看| 成人手机av| 日韩电影二区| 天天躁日日躁夜夜躁夜夜| 欧美人与性动交α欧美精品济南到| 久久婷婷青草| 一本久久精品| 日韩伦理黄色片| 国产又色又爽无遮挡免| 夜夜骑夜夜射夜夜干| 99精品久久久久人妻精品| kizo精华| 色婷婷久久久亚洲欧美| 日韩不卡一区二区三区视频在线| 国产成人一区二区在线| a 毛片基地| 国产成人欧美在线观看 |