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

    靜電排斥表面誘導(dǎo)溶菌酶分子站立

    2013-09-21 09:00:12
    物理化學(xué)學(xué)報 2013年4期
    關(guān)鍵詞:溶菌酶勢能電荷

    白 姝 李 浩 張 麟

    (天津大學(xué)化工學(xué)院生物工程系,天津300072; 天津大學(xué)系統(tǒng)生物工程教育部重點實驗室,天津300072)

    1 引言

    基因工程技術(shù)的飛速發(fā)展推動重組蛋白質(zhì)藥物在微生物體內(nèi)的大規(guī)模生產(chǎn).1-3但是,重組蛋白質(zhì)在原核細(xì)胞中的高水平表達(dá)往往導(dǎo)致其發(fā)生錯誤的折疊和聚集,4形成被稱為包含體(inclusion body)的無活性聚集體,需要通過蛋白質(zhì)復(fù)性過程使其恢復(fù)天然構(gòu)型和生物學(xué)活性.5,6因此,蛋白質(zhì)復(fù)性技術(shù)開發(fā)和理論完善對蛋白質(zhì)藥物生產(chǎn)具有重要意義.7

    蛋白質(zhì)復(fù)性的最大挑戰(zhàn)就是創(chuàng)造合適的復(fù)性條件,最大程度地抑制聚集體生成,從而提高復(fù)性收率.8而現(xiàn)階段重組蛋白質(zhì)復(fù)性仍存在收率低,操作條件較為嚴(yán)苛等問題.近期,Wang等9通過實驗研究發(fā)現(xiàn),帶有正電荷的離子交換介質(zhì)Q和DEAE Sepharose FF可以顯著促進(jìn)帶正電的溶菌酶復(fù)性,而帶有負(fù)電荷的離子交換介質(zhì)SP和CM Sepharose FF可以顯著促進(jìn)帶負(fù)電的牛血清白蛋白復(fù)性,并推斷與蛋白質(zhì)帶有同種電荷的離子交換介質(zhì)可以誘導(dǎo)蛋白質(zhì)在介質(zhì)表面形成定向排列.此定向排列促使蛋白質(zhì)之間的靜電排斥力增大,從而抑制蛋白質(zhì)的聚集而促進(jìn)復(fù)性.但是,通過實驗研究很難直接證實此機理.而分子動力學(xué)(MD)模擬10,11能夠直接提供系統(tǒng)的微觀結(jié)構(gòu)信息,并能通過數(shù)學(xué)統(tǒng)計等手段計算系統(tǒng)的宏觀性質(zhì),已經(jīng)成為與理論研究和實驗研究并行互補的科學(xué)研究手段,廣泛用于蛋白質(zhì)相關(guān)領(lǐng)域的研究.12迄今為止,已有大量研究者通過分子模擬研究界面吸引誘導(dǎo)的小分子13,14和蛋白質(zhì)15-23吸附行為,但很少研究界面排斥作用誘導(dǎo)的蛋白質(zhì)空間分布.

    因此,本文將構(gòu)建靜電排斥表面模型以模擬帶同種電荷的離子交換介質(zhì),通過分子動力學(xué)模擬考察蛋白質(zhì)與靜電排斥表面之間的相互作用過程,展示蛋白質(zhì)在靜電排斥表面上的空間取向和構(gòu)象轉(zhuǎn)換,考察表面所帶電荷數(shù)的影響規(guī)律,揭示同電荷離子交換介質(zhì)輔助蛋白質(zhì)復(fù)性的微觀機理,推動蛋白質(zhì)在荷電表面折疊和分子相互作用研究.

    2 模擬方法

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

    根據(jù)蛋白質(zhì)數(shù)據(jù)庫(http://www.rcsb.org/pdb/)中的晶體衍射結(jié)構(gòu)(PDB ID:3LYZ)構(gòu)建溶菌酶的全原子模型,24如圖1A所示.溶菌酶由129個殘基組成,pH=7時帶8個正電荷(包括17個帶正電荷的殘基和9個帶負(fù)電荷的殘基).

    本文構(gòu)建靜電排斥表面模型以模擬同電荷離子交換介質(zhì),包括基質(zhì)和配基兩部分(參比體系只含基質(zhì)).基質(zhì)由正六邊形網(wǎng)格結(jié)構(gòu)的平板表示,其中正六邊形的邊長為0.152 nm.配基的原子組成參照商業(yè)化介質(zhì)的化學(xué)結(jié)構(gòu)(https://www.gelifesciences.com/)確定,其鍵長、鍵角和二面角及其力常數(shù)等結(jié)構(gòu)參數(shù)由PRODRG 2.5(http://davapc1.bioch.dundee.ac.uk/cig-bin/prodrg_beta)25生成,而原子的部分電荷通過量子化學(xué)計算和類似結(jié)構(gòu)比對調(diào)整確定,標(biāo)示于圖1B中.其后,將配基連接到基質(zhì),根據(jù)實驗條件調(diào)整配基密度為0.22 mmol·mL-1,如圖1C所示.

    圖1 溶菌酶和靜電排斥表面的全原子模型Fig.1 All-atom models of the lysozyme and the surface with electrostatic repulsion

    構(gòu)建模擬體系時,首先將靜電排斥表面放入6.84 nm×6.91 nm×18.00 nm的長方體盒子底部(參比體系的盒子高度為10.00 nm).然后將1個溶菌酶分子放置于表面上方,其質(zhì)心距離表面約2 nm.進(jìn)而添加水分子,并添加適量反離子使模擬體系呈電中性.最后將模擬體系放入6.84 nm×6.91 nm×100 nm的長方體盒子中心進(jìn)行模擬.本文共構(gòu)建4個模擬體系,包括靜電排斥表面體系(標(biāo)記為Q),參比體系(標(biāo)記為control),介質(zhì)所帶電荷增大30%(標(biāo)記為1.3Q)和減小25%(標(biāo)記為0.75Q,體系高度為12 nm)的體系,其具體組成列于表1.

    2.2 分子動力學(xué)模擬

    分子動力學(xué)模擬利用GROMACS 4.0.5程序(http://www.gromacs.org/)26,27完成,使用GROMOS96 43A1力場.28模擬采用NVT系綜,通過v-rescale溫度偶聯(lián)法29控制溫度為298.15 K.采用半步蛙跳法積分,步長為2 fs.靜電相互作用計算采用PME(particle-mesh Ewald)算法,其截斷距離設(shè)為1.2 nm.范德華作用計算采用cut-off算法,其截斷半徑設(shè)為1.2 nm.近鄰原子列表截斷半徑設(shè)為1.2 nm.粒子運動的初始速度根據(jù)模擬溫度298.15 K下的Maxwell分布產(chǎn)生.對模擬盒子施加x,y,z的三個方向的周期性邊界條件.

    本文采用最速下降法進(jìn)行能量最小化,能量閥值設(shè)置為2000 kJ·mol-1.然后,進(jìn)行20 ns的分子動力學(xué)模擬,得到模擬軌跡用于后續(xù)分析.為了保證模擬結(jié)果的準(zhǔn)確性,本文中對體系Q執(zhí)行3次平行模擬.所有模擬均使用曙光TC2600刀片服務(wù)器完成.文中蛋白質(zhì)構(gòu)象圖利用RASMOL軟件30繪制.

    2.3 數(shù)據(jù)分析方法

    2.3.1 蛋白質(zhì)空間取向分析

    本文定義角度θ為蛋白質(zhì)的偶極和z坐標(biāo)軸之間的夾角,以定量描述蛋白質(zhì)在靜電排斥表面上的空間取向.當(dāng)θ=0°時,蛋白質(zhì)的偶極垂直于靜電排斥表面(平行于+z坐標(biāo)軸).當(dāng)θ為90°或-90°時,蛋白質(zhì)的偶極平行于靜電排斥表面(偶極指向+x方向定義為正).通過自編程序計算蛋白質(zhì)的θ值隨模擬時間的變化.

    表1 模擬體系的參數(shù)Table 1 Parameters of simulation systems

    2.3.2 分子間勢能分析

    利用GROMACS軟件包的g_energy程序進(jìn)行勢能分析,包括Lennard-Jones(LJ)勢能(標(biāo)記為ELJ)和庫侖勢能(標(biāo)記為EC).

    2.3.3 二級結(jié)構(gòu)分析

    為了描述MD模擬過程中蛋白質(zhì)的構(gòu)象轉(zhuǎn)化,利用Define Secondary Structure of Proteins(DSSP)方法31及GROMACS軟件包的do_dssp程序完成其二級結(jié)構(gòu)分析.

    2.3.4 均方根偏差計算

    利用GROMACS軟件包的g_rms程序計算溶菌酶分子相比于其天然結(jié)構(gòu)的均方根偏差(RMSD),如式(1)所示.RMSD值越小,表明溶菌酶天然結(jié)構(gòu)保持更好.

    2.3.5 回轉(zhuǎn)半徑計算

    利用GROMACS軟件包的g_gyrate程序計算溶菌酶的回轉(zhuǎn)半徑(Rg)以描述其結(jié)構(gòu)緊密程度,如式(2)所示.Rg值越小,表明蛋白質(zhì)的結(jié)構(gòu)越緊密.

    2.3.6 自由能表面圖

    參照文獻(xiàn)20,32計算蛋白質(zhì)處于特定θ和zcom的概率分布P(θ,zcom),其中zcom為蛋白質(zhì)質(zhì)心的z坐標(biāo)值.通過模擬軌跡采樣和自編程序計算完成.

    3 結(jié)果與討論

    3.1 溶菌酶空間取向分析

    首先通過系統(tǒng)快照(snapshots)直觀分析模擬體系的微觀構(gòu)象隨模擬時間的變化,包括溶菌酶的構(gòu)象,空間取向以及位置等,如圖2所示.此處,著重分析溶菌酶的偶極取向,偶極已放大以便清晰顯示,在圖中以棒狀標(biāo)出,其偶極方向由紅色指向藍(lán)色(印刷版中由黑色指向灰色).注意,此處僅顯示一條典型軌跡結(jié)果,其余平行模擬軌跡的分析結(jié)果一致,故略去.

    圖2 溶菌酶在靜電排斥表面上的構(gòu)象隨模擬時間變化圖Fig.2 Snapshots of lysozyme over the surface with electrostatic repulsion as a function of simulation time

    圖2A顯示,溶菌酶的初始位置靠近靜電排斥表面,其偶極方向接近平行方向,溶菌酶含有豐富的α-螺旋和β-折疊等二級結(jié)構(gòu).其后,蛋白質(zhì)受到表面的排斥作用而遠(yuǎn)離.2 ns時,蛋白質(zhì)已經(jīng)遠(yuǎn)離表面并呈現(xiàn)“站立”姿勢,其偶極方向接近垂直方向.此時蛋白質(zhì)整體結(jié)構(gòu)被拉伸,其二級結(jié)構(gòu)部分喪失而轉(zhuǎn)變?yōu)橄鄬θ嵝缘木砬Y(jié)構(gòu).6 ns時,蛋白質(zhì)依然遠(yuǎn)離表面,其偶極方向有輕微偏轉(zhuǎn),但依然接近垂直方向,即蛋白質(zhì)一直保持“站立”姿勢.其后,由于蛋白質(zhì)遠(yuǎn)離表面,且偶極旋轉(zhuǎn)會進(jìn)一步降低表面對其的排斥作用,蛋白質(zhì)可以再次靠近表面.例如,12 ns時,蛋白質(zhì)略微靠近表面,而其偶極旋轉(zhuǎn)到幾乎完全垂直的方向,表明蛋白質(zhì)靠近表面時,靜電相互作用增強利于誘導(dǎo)形成“站立”姿勢.但同時也需注意,20 ns時,蛋白質(zhì)進(jìn)一步接近表面,但其偶極反而偏離垂直方向,表明蛋白質(zhì)的取向還受到其他因素,如溶劑的影響,存在一定波動.

    在參比體系(圖2B)中,由于表面呈電中性,和蛋白質(zhì)之間不存在靜電相互作用.因此,蛋白質(zhì)一直保持在表面附近,在溶劑作用下存在一定的移動和轉(zhuǎn)動,但其偶極方向變化并無明顯規(guī)律.同時,在參比體系中,蛋白質(zhì)的二級結(jié)構(gòu)一直很好保持,表明蛋白質(zhì)自身結(jié)構(gòu)在溶液中比較穩(wěn)定.

    因此,微觀構(gòu)象分析表明,靜電排斥表面推動蛋白質(zhì)遠(yuǎn)離.在此過程中,表面誘導(dǎo)蛋白質(zhì)逐漸“站立”,使其偶極保持在垂直于表面的方向.盡管蛋白質(zhì)的偶極方向存在一定的擺動,但總體而言都接近垂直方向,即蛋白質(zhì)基本保持“站立”姿勢.

    進(jìn)而考察角度θ(定義和計算流程見2.3.1節(jié))的變化以定量描述蛋白質(zhì)在靜電排斥表面作用下的空間取向分布及其變化情況,同時通過蛋白質(zhì)質(zhì)心的z坐標(biāo)值定量描述蛋白質(zhì)與表面的相對位置變化,如圖3所示.此處,體系Q的數(shù)據(jù)為3次平行模擬結(jié)果的平均值.

    在靜電排斥表面上,蛋白質(zhì)的θ值從初始的-105°迅速變化到0°左右,即蛋白質(zhì)從“平躺”姿勢迅速轉(zhuǎn)變?yōu)椤罢玖ⅰ弊藙?同時其zcom迅速增大,表明蛋白質(zhì)迅速遠(yuǎn)離表面.相比而言,zcom變化更緩慢,即其變化速率小于θ的變化速率,表明蛋白質(zhì)在受到排斥而遠(yuǎn)離表面的過程中即已形成“站立”姿勢.在2 ns附近時,蛋白質(zhì)的zcom呈現(xiàn)峰值,此時蛋白質(zhì)的θ在0°附近波動,但都維持在-60°到60°的區(qū)間內(nèi),表明蛋白質(zhì)在遠(yuǎn)離表面時,由于靜電排斥減弱,蛋白質(zhì)的偶極方向存在一定的擺動.其后,蛋白質(zhì)的zcom值略微降低,在3.0到5.5 ns期間呈現(xiàn)平臺期,其zcom值維持在8 nm,而此時蛋白質(zhì)的θ值基本維持在0°,表明蛋白質(zhì)依然很好維持“站立”姿勢.此后,蛋白質(zhì)略微接近表面,而θ值的變化幅度增大,但依然維持在0°附近.到12 ns后,蛋白質(zhì)的θ值變化幅度進(jìn)一步增大,表明蛋白質(zhì)從遠(yuǎn)離靜電排斥表面的位置逐漸靠近表面時,由于靜電排斥作用增強,蛋白質(zhì)的偶極方向需要重新調(diào)整以達(dá)到新的平衡狀態(tài).其概率分布(圖3B)表明,蛋白質(zhì)的zcom主要介于5到10 nm之間,表明蛋白質(zhì)受到排斥而遠(yuǎn)離表面,其θ角主要分布于-60°到60°之間,確證蛋白質(zhì)形成“站立”姿勢.

    圖3 溶菌酶空間取向和位置隨模擬時間變化圖Fig.3 Orientation and the location of lysozyme as a function of simulation time

    在參比體系中,蛋白質(zhì)一直維持在表面附近,其zcom值變化較小,可以認(rèn)為蛋白質(zhì)只是在其初始位置附近波動而并無明顯移動(圖3C的統(tǒng)計結(jié)果給出明確證據(jù)).在模擬初始階段,蛋白質(zhì)的θ值也無明顯變化,只是在其初始值附近波動.2.5 ns后蛋白質(zhì)的θ值發(fā)生明顯變化,但分析表明都是無規(guī)的轉(zhuǎn)動導(dǎo)致,從圖2的系統(tǒng)快照中也可得到直觀證明.15 ns后,蛋白質(zhì)的偶極方向又恢復(fù)到初始狀態(tài)并存在輕微波動.

    因此,θ值和zcom值隨模擬時間的變化確證蛋白質(zhì)在靜電排斥表面作用下迅速遠(yuǎn)離,而在遠(yuǎn)離過程中其偶極旋轉(zhuǎn),形成特定的“站立”姿勢.此結(jié)果直接驗證實驗推測,表明蛋白質(zhì)確實在靜電排斥表面形成定向排布,有助于增強蛋白質(zhì)之間的靜電排斥作用而抑制溶菌酶聚集.但同時需注意,蛋白質(zhì)和表面間的距離,以及溶劑均影響蛋白質(zhì)的偶極方向使其存在一定的擺動.

    然后,通過蛋白質(zhì)和表面間相互作用勢能分析,考察蛋白質(zhì)形成特定“站立”姿勢的關(guān)鍵作用力,如圖4所示.此處,體系Q的數(shù)據(jù)為3次平行模擬結(jié)果的平均值.

    圖4 溶菌酶和介質(zhì)間的LJ勢能和庫侖勢能隨模擬時間變化圖Fig.4 Lennard-Jones(LJ)and Coulomb potential energies between lysozyme and resin as a function of simulation time

    僅在模擬初始階段觀察到靜電排斥表面和蛋白質(zhì)之間存在很小的LJ勢能,約-1.2 kJ·mol-1.其后,LJ勢能持續(xù)為0,這是因為蛋白質(zhì)被排斥遠(yuǎn)離表面(圖3)而LJ勢能是短程相互作用,隨距離迅速衰減.而蛋白質(zhì)和表面之間存在很大的靜電相互作用勢能,接近140 MJ·mol-1.靜電相互作用勢能隨模擬時間迅速減小,到2 ns附近時達(dá)到最小值83 MJ·mol-1,這是因為蛋白質(zhì)在靜電排斥作用下迅速遠(yuǎn)離表面(圖3),使靜電勢能迅速降低.其后,隨著蛋白質(zhì)再次靠近表面(圖3),蛋白質(zhì)和表面間的靜電相互作用勢能再次上升,最終達(dá)到105 MJ·mol-1.而在參比體系中,蛋白質(zhì)與表面之間的LJ勢能和靜電作用勢能都基本為0,因而表面對蛋白質(zhì)并無作用,蛋白質(zhì)呈現(xiàn)溶液主體行為,與圖2和圖3的結(jié)果完全吻合.

    因此,通過勢能分析可以確證,表面通過靜電排斥作用驅(qū)動蛋白質(zhì)遠(yuǎn)離并誘導(dǎo)其形成特定“站立”姿勢.為驗證此結(jié)論,本文還構(gòu)建帶負(fù)電的靜電表面,發(fā)現(xiàn)帶正電的溶菌酶分子在靜電吸引作用下直接吸附于負(fù)電表面(數(shù)據(jù)未顯示).如果通過調(diào)整溶液pH值使溶菌酶也帶負(fù)電(如pH=14),則溶菌酶同樣在靜電排斥作用下迅速遠(yuǎn)離,再次確證靜電排斥作用是關(guān)鍵驅(qū)動力.但同時也發(fā)現(xiàn),pH=14時,溶菌酶中只有帶負(fù)電的氨基酸殘基而無帶正電的氨基酸殘基,無法形成明顯偶極而導(dǎo)致無法形成特定“站立”姿勢(數(shù)據(jù)未顯示),表明正負(fù)電荷的不均勻分布是特定“站立”姿勢形成的關(guān)鍵因素.

    3.2 溶菌酶構(gòu)象分析

    利用DSSP方法確定蛋白質(zhì)二級結(jié)構(gòu)隨模擬時間的變化,如圖5所示.結(jié)果表明,溶菌酶的主要二級結(jié)構(gòu)是α-螺旋,間雜零星β-折疊結(jié)構(gòu).參比體系中,各顏色條帶隨模擬進(jìn)行都比較穩(wěn)定.例如,10號殘基、30號殘基和95號殘基處的α-螺旋在模擬過程中一直保持完好,而110號殘基處的α-螺旋結(jié)構(gòu)存在短暫消失的情況,但很快重新生成,且一直維持到模擬終點,表明蛋白質(zhì)的二級結(jié)構(gòu)得到很好保持.

    但在靜電排斥表面體系中,各顏色條帶均存在消失或增加的反復(fù)過程,總體呈現(xiàn)不連續(xù)分布狀態(tài),表明蛋白質(zhì)的二級結(jié)構(gòu)存在明顯變化.其變化細(xì)節(jié)表明蛋白質(zhì)中多處的α-螺旋結(jié)構(gòu)喪失.參照圖2的系統(tǒng)快照分析,在靜電排斥表面作用下,蛋白質(zhì)中帶正電荷的殘基受到排斥而試圖遠(yuǎn)離表面,但帶負(fù)電荷的殘基受到吸引而試圖靠近表面,因此蛋白質(zhì)受力不均而被輕微拉伸,其內(nèi)部原子相對位置的變化造成其二級結(jié)構(gòu)的破壞.從圖5中也可看出,帶電殘基處的二級結(jié)構(gòu)變化最為劇烈.例如,天然溶菌酶中ARG5-LYS13形成α-螺旋結(jié)構(gòu),在參比體系中此結(jié)構(gòu)很好保持(圖5B).而在靜電排斥表面體系中,由于ARG5和LYS13均帶電,受到表面的影響而導(dǎo)致此螺旋結(jié)構(gòu)不穩(wěn)定,甚至整體螺旋結(jié)構(gòu)被破壞(圖5A).同時需注意,隨著蛋白質(zhì)位置和偶極方向的調(diào)整,蛋白質(zhì)構(gòu)象也存在細(xì)微調(diào)整過程.

    因此,計算蛋白質(zhì)的RMSD和Rg值模擬時間的變化,以通過結(jié)構(gòu)偏差程度及其整體結(jié)構(gòu)緊密程度(詳見2.3節(jié)中數(shù)據(jù)分析方法描述)定量考察蛋白質(zhì)的構(gòu)象變化,如圖6所示.此處數(shù)據(jù)為3次平行模擬結(jié)果的平均值.

    參比體系中,蛋白質(zhì)的RMSD值在模擬開始時輕微增大,表明蛋白質(zhì)在溶液環(huán)境中存在細(xì)微結(jié)構(gòu)調(diào)整,其后則一直在0.10 nm上下波動.蛋白質(zhì)的Rg值隨模擬時間輕微降低,然后維持在1.35 nm上下波動.二者的數(shù)值都很小且波動幅度也很小,表明參比體系中蛋白質(zhì)并無明顯構(gòu)象變化,與圖2的系統(tǒng)快照和圖5的二級結(jié)構(gòu)分析結(jié)果一致.

    而在靜電排斥表面體系中,蛋白質(zhì)的RMSD值在模擬初始階段即迅速增大,并進(jìn)而維持在0.5 nm上下波動且持續(xù)到模擬終點,表明蛋白質(zhì)確實在模擬初始的排斥遠(yuǎn)離過程中受到很大的結(jié)構(gòu)擾動,使其結(jié)構(gòu)明顯偏離天然態(tài).蛋白質(zhì)的Rg值在模擬初始階段也迅速上升,高于參比體系,表明蛋白質(zhì)結(jié)構(gòu)確實被拉伸而使其整體結(jié)構(gòu)松散.其后,Rg值很快下降到與參比體系相當(dāng),進(jìn)而繼續(xù)下降并穩(wěn)定在1.30 nm上下波動,其穩(wěn)定值低于參比體系中的Rg值.綜合RMSD值和Rg值的變化表明蛋白質(zhì)被排斥遠(yuǎn)離表面后,由于距離增大(圖3)導(dǎo)致靜電排斥力減弱(圖4),蛋白質(zhì)受到外界干擾減小,開始構(gòu)象調(diào)整,呈現(xiàn)整體塌縮的變化趨勢.然而,蛋白質(zhì)的最終結(jié)構(gòu)仍然偏離天然結(jié)構(gòu),但其偏差數(shù)值并不大,表明其最終骨架結(jié)構(gòu)還能較好保持,與圖2的系統(tǒng)快照分析結(jié)果一致.因而,靜電排斥表面對溶菌酶的分子結(jié)構(gòu)有一定擾動,盡管整體骨架結(jié)構(gòu)變化不大,也可能造成其活性喪失.

    因此,綜合分析二級結(jié)構(gòu),RMSD和Rg值的結(jié)果表明,參比體系中,蛋白質(zhì)的構(gòu)象很好保持.而在靜電排斥表面體系中,蛋白質(zhì)由于受力不均,存在先拉伸后塌縮的變化過程,其二級結(jié)構(gòu)遭到破壞,但主體骨架結(jié)構(gòu)較好保持.

    3.3 表面所帶電荷數(shù)的影響分析

    圖5 分子動力學(xué)模擬中溶菌酶分子的二級結(jié)構(gòu)演化Fig.5 Secondary structure evolution of lysozyme during the molecular dynamics simulation

    上述結(jié)果表明,靜電排斥表面能夠誘導(dǎo)蛋白質(zhì)形成“站立”姿勢,而靜電排斥作用是主要驅(qū)動力.因此,表面所帶電荷數(shù)對蛋白質(zhì)的空間取向應(yīng)該具有重要影響.本文構(gòu)建電荷數(shù)增大的靜電排斥表面模型(1.3Q)以及電荷數(shù)減小(0.75Q)的模型,分析其與蛋白質(zhì)的相互作用及其對蛋白質(zhì)的空間取向和構(gòu)象變化的影響,如圖7、8、9所示.

    圖7顯示,在體系1.3Q中,蛋白質(zhì)的θ值一直維持在0°上下波動,且波動幅度很小,表明表面所帶電荷數(shù)增大使其靜電排斥作用增強,蛋白質(zhì)能夠更好形成并保持“站立”姿勢.蛋白質(zhì)的zcom值迅速增大,到6 nm處增大變緩,最終維持在10 nm上下波動.其概率分布(圖7B)表明,蛋白質(zhì)的θ角更集中于-60°到60°之間.因此,相比于體系Q(圖3),蛋白質(zhì)的排斥遠(yuǎn)離過程更快且最終穩(wěn)定距離更遠(yuǎn),但其“站立”姿勢也更穩(wěn)定,表明表面所帶電荷數(shù)增大確實有利于蛋白質(zhì)的“站立”姿勢形成和穩(wěn)定.

    圖6 溶菌酶中Cα原子的RMSD和Rg隨模擬時間變化圖Fig.6 RMSD and the Rgvalues for Cαatoms of lysozyme as a function of simulation time

    而在體系0.75Q中,蛋白質(zhì)的θ值主要維持在-50°上下波動且波動幅度很大.在17至20 ns期間,蛋白質(zhì)遠(yuǎn)離表面,其θ值發(fā)生明顯變化,在-120°至-60°區(qū)間內(nèi)波動,即偶極平行于表面,蛋白質(zhì)呈現(xiàn)“平躺”姿勢.蛋白質(zhì)的zcom值到1.6 ns時增大到5.7 nm,變化比體系1.3Q和體系Q中緩慢.其后,zcom值降低,到5 ns時再次上升,最終達(dá)到7.2 nm.因此,體系0.75Q中,模擬終點時zcom值小于體系1.3Q的情況,表明靜電排斥減弱,蛋白質(zhì)穩(wěn)定位置更靠近表面.但需要注意的是,雖然體系0.75Q中模擬終點的zcom值大于體系Q,但其概率分布(圖7C)表明,在zcom值較大區(qū)域蛋白質(zhì)的θ角偏離-60°到60°區(qū)域,表明蛋白質(zhì)此時并不存在“站立”姿勢.而其主要分布區(qū)域是在zcom=5 nm處,此時其θ=-60°,表明體系0.75Q中,蛋白質(zhì)也形成“站立”取向,但趨勢更弱.

    圖8顯示,在1.3Q和0.75Q兩體系中,蛋白質(zhì)和表面之間的LJ勢能都基本為0,與體系Q相同,表明在靜電排斥表面體系中,蛋白質(zhì)和表面間距離較遠(yuǎn),而LJ勢能為短程相互作用,隨距離增大迅速衰減.而靜電相互作用很強,在兩體系中均呈現(xiàn)迅速下降,然后趨于平穩(wěn)的變化趨勢,分別穩(wěn)定在83和19 MJ·mol-1.勢能變化表明蛋白質(zhì)的位置和空間取向變化受靜電排斥作用調(diào)控,且受到表面所帶電荷數(shù)影響.表面電荷數(shù)增大,蛋白質(zhì)的平衡位置更遠(yuǎn)離表面,且利于“站立”姿勢的形成和維持.

    圖9顯示,在1.3Q和0.75Q兩體系中,蛋白質(zhì)的RMSD值均隨模擬進(jìn)行而迅速上升,然后分別在0.70和0.15 nm上下波動.而兩體系中蛋白質(zhì)的Rg值略微不同,體系1.3Q中Rg值高于初始天然態(tài),且存在顯著增大再回落的過程(如2.6 ns時).而體系0.75Q中Rg值基本保持在天然態(tài)的Rg值.結(jié)合RMSD和Rg值分析表明,隨著所帶電荷數(shù)增大,表面對蛋白質(zhì)的結(jié)構(gòu)擾動加劇.反之,電荷數(shù)降低的表面對蛋白質(zhì)的結(jié)構(gòu)擾動減弱,使蛋白質(zhì)能夠更好保持其天然結(jié)構(gòu).

    圖7 介質(zhì)電荷數(shù)對溶菌酶的空間取向和位置的影響Fig.7 Effects of charge number of resin on the orientation and location of lysozyme

    圖8 介質(zhì)電荷數(shù)對溶菌酶和介質(zhì)間相互作用勢能的影響Fig.8 Effects of charge number of resin on the interaction energies between the lysozyme and the resin

    4 結(jié)論

    圖9 介質(zhì)電荷數(shù)對溶菌酶構(gòu)象的影響Fig.9 Effects of charge number of resin on the conformation of lysozyme

    以同電荷離子交換介質(zhì)促進(jìn)溶菌酶等蛋白質(zhì)復(fù)性的實驗結(jié)果為背景,通過靜電排斥表面模擬同電荷離子交換介質(zhì),構(gòu)建溶菌酶在靜電排斥表面體系中的全原子模型,利用分子動力學(xué)模擬考察蛋白質(zhì)在表面上的空間取向和構(gòu)象變化,并通過相互作用勢能分析深入剖析靜電排斥表面的作用機理和細(xì)節(jié).結(jié)果表明,蛋白質(zhì)在表面的靜電排斥作用下迅速遠(yuǎn)離,在此過程中逐漸形成其偶極垂直于表面的“站立”姿勢.此定向排布有助于增強蛋白質(zhì)之間的靜電排斥作用而抑制溶菌酶聚集.但是,由于蛋白質(zhì)中帶正電荷的殘基和帶負(fù)電荷的殘基受力方向相反,導(dǎo)致蛋白質(zhì)結(jié)構(gòu)拉伸而破壞其二級結(jié)構(gòu).研究發(fā)現(xiàn)表面所帶電荷數(shù)對蛋白質(zhì)的空間取向和構(gòu)象變化有顯著影響.電荷數(shù)增大使得表面對蛋白質(zhì)的靜電排斥作用增強,能夠誘導(dǎo)蛋白質(zhì)更快更好地形成“站立”姿勢.而表面所帶電荷數(shù)降低,則不利于蛋白質(zhì)維持“站立”姿勢.同時,本文的模擬中尚有四個影響因素需要在后續(xù)研究中考慮.其一,模擬體系中的大量溶劑分子會影響蛋白質(zhì)的構(gòu)象和空間取向調(diào)整.其二,本文研究的排斥性界面不同于吸附性界面.吸附性界面作用下,蛋白質(zhì)靠近界面使得相互作用增強,因而利于形成明顯的空間取向分布.而排斥性界面作用下,蛋白質(zhì)遠(yuǎn)離界面,使得相互作用減弱,不利于蛋白質(zhì)形成明顯空間取向,造成排斥性界面研究的困難.其三,限于計算量和模擬體系復(fù)雜度的考慮,本文僅模擬單個溶菌酶分子,能夠確證靜電排斥表面誘導(dǎo)蛋白質(zhì)形成定向排布,但無法直接闡釋表面對蛋白質(zhì)分子間相互作用的影響,將在后續(xù)研究中開展.其四,本文直接用晶體結(jié)構(gòu)建模計算,后續(xù)研究中可以先通過一定時間的模擬獲得溶菌酶在“自由”狀態(tài)下的構(gòu)象以消除溶液條件和力場的影響.總體而言,本文通過分子動力學(xué)模擬完整闡釋蛋白質(zhì)在靜電排斥表面上的空間取向和構(gòu)象變化過程,考察此過程的微觀細(xì)節(jié),并能確定表面電荷數(shù)等參數(shù)的影響規(guī)律,將有助于推動蛋白質(zhì)在荷電表面折疊和分子相互作用研究.

    (1) Durocher,Y.;Butler,M.Curr.Opin.Biotech.2009,20,700.doi:10.1016/j.copbio.2009.10.008

    (2)Manning,M.C.;Chou,D.K.;Murphy,B.M.;Payne,R.W.;Katayama,D.S.Pharm.Res.2010,27,544.doi:10.1007/s11095-009-0045-6

    (3) Ferrer-Miralles,N.;Domingo-Espin,J.;Corchero,J.L.;Vazquez,E.;Villaverde,A.Microb.Cell.Fact.2009,8,17.doi:10.1186/1475-2859-8-17

    (4) Fink,A.L.Fold.Des.1998,3,R9.

    (5) Middelberg,A.R.Trends Biotechnol.2002,20,437.doi:10.1016/S0167-7799(02)02047-4

    (6) Baldwin,R.L.J.Mol.Biol.2007,371,283.doi:10.1016/j.jmb.2007.05.078

    (7)Chen,Y.W.;Ding,F(xiàn).;Nie,H.F.;Serohijos,A.W.;Sharma,S.;Wilcox,K.C.;Yin,S.Y.;Dokholyan,N.V.Arch.Biochem.Biophys.2008,469,4.doi:10.1016/j.abb.2007.05.014

    (8) Lu,D.;Liu,Z.Annu.Rep.Prog.Chem.,Sect.C 2010,106,259.doi:10.1039/b903487k

    (9)Wang,G.;Dong,X.;Sun,Y.Biotechnol.Bioeng.2011,108,1068.doi:10.1002/bit.23038

    (10) Karplus,M.Biopolymers 2003,68,350.

    (11) Karplus,M.;McCammon,J.A.Nat.Struct.Biol.2002,9,646.doi:10.1038/nsb0902-646

    (12) Bai,H.J.;Lai,L.H.Acta Phys.-Chim.Sin.2010,26,1988.[白紅軍,來魯華.物理化學(xué)學(xué)報,2010,26,1988.]doi:10.3866/PKU.WHXB20100725

    (13)Zhang,X.R.;Wang,W.C.Acta Phys.-Chim.Sin.2002,18,680.[張現(xiàn)仁,汪文川.物理化學(xué)學(xué)報,2002,18,680.]doi:10.3866/PKU.WHXB20020803

    (14) Xiang,Z.H.;Wang,W.C.;Cao,D.P.Scientia Sinica Chimica 2012,42,235.[向中華,汪文川,曹達(dá)鵬.中國科學(xué):化學(xué),2012,42,235.]doi:10.1360/032011-297

    (15) Kubiak-Ossowska,K.;Mulheran,P.A.Langmuir 2010,26,15954.doi:10.1021/la102960m

    (16) Ravichandran,S.;Madura,J.D.;Talbot,J.J.Phys.Chem.B 2001,105,3610.doi:10.1021/jp010223r

    (17) Carlsson,F(xiàn).;Hyltner,E.;Arnebrant,T.;Malmsten,M.;Linse,P.J.Phys.Chem.B 2004,108,9871.doi:10.1021/jp0495186

    (18)Slater,G.W.;Holm,C.;Chubynsky,M.V.;de Haan,H.H.;Dube,A.;Grass,K.;Hickey,O.A.;Kingsburry,C.;Sean,D.;Shendruk,T.N.;Nhan,L.X.Electrophoresis 2009,30,792.doi:10.1002/elps.v30:5

    (19) Chandrasekhar,I.;Kastenholz,M.;Lins,R.D.;Oostenbrink,C.;Schuler,L.D.;Tieleman,D.P.;van Gunsteren,W.F.Eur.Biophys.J.Biophy.2003,32,67.

    (20) Zhang,L.;Zhao,G.F.;Sun,Y.J.Phys.Chem.B 2009,113,6873.doi:10.1021/jp809754k

    (21) Zhang,L.;Zhao,G.F.;Sun,Y.J.Phys.Chem.B 2010,114,2203.doi:10.1021/jp903852c

    (22)Kang,K.;Lu,D.N.;Liu,Z.Chin.J.Chem.Eng.2012,20,284.

    (23) Shen,L.M.;Lin,D.Q.;Mei,L.H.;Yao,S.J.Chinese Journal of Bioprocess Engineering 2006,4,29.[沈立民,林東強,梅樂和,姚善涇.生物加工過程,2006,4,29.]

    (24) Diamond,R.J.Mol.Biol.1974,82,371.doi:10.1016/0022-2836(74)90598-1

    (25)Schuttelkopf,A.W.;vanAaltern,D.M.F.Acta Crystallogr.D 2004,60,1355.doi:10.1107/S0907444904011679

    (26) Berendsen,H.J.;van der Spoel,D.;van Drunen,R.Comput.Phys.Commun.1995,91,43.doi:10.1016/0010-4655(95)00042-E

    (27) Lindahl,E.;Hess,B.;van der Spoel,D.J.Mol.Model.2001,7,306.

    (28)Mackerell,A.D.J.Comput.Chem.2004,25,1584.

    (29) Bussi,G.;Donadio,D.;Parrinello,M.J.Chem.Phys.2007,126,014101.doi:10.1063/1.2408420

    (30) Sayle,R.;Milnerwhite,E.Trends Biochem.Sci.1995,20,374.doi:10.1016/S0968-0004(00)89080-5

    (31) Kabsch,W.;Sander,C.Biopolymers 1983,22,2577.

    (32)Li,W.F.;Zhang,J.;Wang,J.;Wang,W.J.Am.Chem.Soc.2008,130,892.doi:10.1021/ja075302g

    猜你喜歡
    溶菌酶勢能電荷
    連續(xù)分布電荷體系電荷元的自能問題*
    物理通報(2024年4期)2024-04-09 12:41:28
    “動能和勢能”知識鞏固
    作 品:景觀設(shè)計
    ——《勢能》
    文化縱橫(2022年3期)2022-09-07 11:43:18
    “動能和勢能”知識鞏固
    “動能和勢能”隨堂練
    電荷知識知多少
    偶氮類食品著色劑誘惑紅與蛋溶菌酶的相互作用研究
    電荷守恒在化學(xué)解題中的應(yīng)用
    靜電現(xiàn)象有什么用?
    動物型溶菌酶研究新進(jìn)展
    天天添夜夜摸| 久久人妻福利社区极品人妻图片| 精品一品国产午夜福利视频| 两个人免费观看高清视频| 国产av一区二区精品久久| 精品人妻1区二区| 韩国精品一区二区三区| 国产在线观看jvid| 午夜福利,免费看| www.999成人在线观看| 国产精品久久久人人做人人爽| 黑人巨大精品欧美一区二区mp4| 一本大道久久a久久精品| 丰满人妻熟妇乱又伦精品不卡| 操出白浆在线播放| 日韩高清综合在线| 制服丝袜大香蕉在线| 国产高清激情床上av| 人人妻人人澡欧美一区二区 | 精品久久久精品久久久| 91老司机精品| 午夜成年电影在线免费观看| 男人操女人黄网站| av超薄肉色丝袜交足视频| 99riav亚洲国产免费| 欧美日韩一级在线毛片| 狂野欧美激情性xxxx| 亚洲欧洲精品一区二区精品久久久| 国产aⅴ精品一区二区三区波| 99国产精品99久久久久| 亚洲国产精品合色在线| 韩国精品一区二区三区| 欧美黑人精品巨大| 老汉色av国产亚洲站长工具| 欧美激情高清一区二区三区| 青草久久国产| 国产精品久久电影中文字幕| 中文字幕高清在线视频| 黄色成人免费大全| 久久久国产欧美日韩av| 亚洲精品美女久久久久99蜜臀| 亚洲中文av在线| 国产成人精品在线电影| 久久久久九九精品影院| 制服诱惑二区| 亚洲国产精品久久男人天堂| 免费一级毛片在线播放高清视频 | 欧美激情 高清一区二区三区| 国产精品免费一区二区三区在线| 国产成人精品在线电影| 视频在线观看一区二区三区| 桃色一区二区三区在线观看| 我的亚洲天堂| 欧美乱码精品一区二区三区| 免费在线观看影片大全网站| 啦啦啦免费观看视频1| 性色av乱码一区二区三区2| 一边摸一边抽搐一进一小说| 国产日韩一区二区三区精品不卡| 777久久人妻少妇嫩草av网站| 天天躁夜夜躁狠狠躁躁| 日韩精品免费视频一区二区三区| 国产一区二区在线av高清观看| 性色av乱码一区二区三区2| 男人操女人黄网站| 国产片内射在线| 欧美日本视频| 亚洲色图av天堂| 国产真人三级小视频在线观看| www.www免费av| 欧美黄色淫秽网站| 精品久久蜜臀av无| 精品国内亚洲2022精品成人| 精品电影一区二区在线| 涩涩av久久男人的天堂| 啪啪无遮挡十八禁网站| 国内精品久久久久久久电影| 国内久久婷婷六月综合欲色啪| 在线观看舔阴道视频| 成人永久免费在线观看视频| 国产xxxxx性猛交| 神马国产精品三级电影在线观看 | 他把我摸到了高潮在线观看| 国产精品九九99| 免费搜索国产男女视频| 大香蕉久久成人网| 91在线观看av| 好男人电影高清在线观看| 精品人妻在线不人妻| 国产精品av久久久久免费| 久久久久久久久免费视频了| 黑人巨大精品欧美一区二区mp4| 亚洲成人精品中文字幕电影| 久久久国产精品麻豆| 人成视频在线观看免费观看| 一进一出抽搐动态| 亚洲在线自拍视频| 动漫黄色视频在线观看| 一区福利在线观看| 久久久久久久午夜电影| 十八禁人妻一区二区| 久久这里只有精品19| 最新美女视频免费是黄的| 无限看片的www在线观看| 黄色 视频免费看| 男女做爰动态图高潮gif福利片 | 久久久久久久久久久久大奶| 精品午夜福利视频在线观看一区| 美国免费a级毛片| 成人av一区二区三区在线看| 欧美老熟妇乱子伦牲交| 日韩高清综合在线| 欧洲精品卡2卡3卡4卡5卡区| 在线十欧美十亚洲十日本专区| 麻豆久久精品国产亚洲av| 一区福利在线观看| 搡老妇女老女人老熟妇| 国产区一区二久久| 日日爽夜夜爽网站| 国产精品98久久久久久宅男小说| 亚洲片人在线观看| av超薄肉色丝袜交足视频| 狂野欧美激情性xxxx| 久久久久久久久中文| 亚洲 欧美 日韩 在线 免费| 精品国产美女av久久久久小说| 免费看美女性在线毛片视频| aaaaa片日本免费| 午夜免费鲁丝| 亚洲中文字幕日韩| 日韩精品免费视频一区二区三区| 欧美成人性av电影在线观看| 性少妇av在线| 99精品在免费线老司机午夜| 亚洲中文字幕一区二区三区有码在线看 | 一二三四在线观看免费中文在| 757午夜福利合集在线观看| a在线观看视频网站| av在线播放免费不卡| 国产男靠女视频免费网站| 女人高潮潮喷娇喘18禁视频| 纯流量卡能插随身wifi吗| 纯流量卡能插随身wifi吗| ponron亚洲| 97碰自拍视频| 精品久久久久久,| 亚洲精品久久成人aⅴ小说| 亚洲性夜色夜夜综合| 久久精品人人爽人人爽视色| 日韩欧美免费精品| 日韩欧美三级三区| 男女床上黄色一级片免费看| 国产一区二区三区在线臀色熟女| 国产高清有码在线观看视频 | 日韩中文字幕欧美一区二区| 黄片小视频在线播放| 在线天堂中文资源库| 精品国产美女av久久久久小说| 亚洲最大成人中文| 国产精品一区二区免费欧美| 亚洲av成人不卡在线观看播放网| 欧美老熟妇乱子伦牲交| 国产成人精品无人区| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久久久免费视频| 熟女少妇亚洲综合色aaa.| 国产成人欧美| 成人亚洲精品av一区二区| 久久久国产成人精品二区| 高清黄色对白视频在线免费看| 久久精品成人免费网站| 亚洲成av片中文字幕在线观看| 免费在线观看完整版高清| www日本在线高清视频| 一卡2卡三卡四卡精品乱码亚洲| 国产成人影院久久av| 亚洲中文字幕日韩| 一边摸一边做爽爽视频免费| 日韩大尺度精品在线看网址 | tocl精华| 国产成人精品无人区| 深夜精品福利| 91麻豆精品激情在线观看国产| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 国产成人一区二区三区免费视频网站| 99国产综合亚洲精品| 亚洲无线在线观看| 国产一区二区激情短视频| 波多野结衣高清无吗| 亚洲全国av大片| 十分钟在线观看高清视频www| 免费在线观看完整版高清| 久久精品国产清高在天天线| 亚洲av第一区精品v没综合| 一个人观看的视频www高清免费观看 | a级毛片在线看网站| 老司机午夜福利在线观看视频| 中文字幕另类日韩欧美亚洲嫩草| 丰满的人妻完整版| 日本精品一区二区三区蜜桃| 久久久久亚洲av毛片大全| 91字幕亚洲| 亚洲欧美精品综合一区二区三区| 老司机在亚洲福利影院| 亚洲 欧美一区二区三区| 91老司机精品| 久久午夜亚洲精品久久| 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美一区二区综合| 18美女黄网站色大片免费观看| 淫秽高清视频在线观看| 一区二区日韩欧美中文字幕| 久久国产乱子伦精品免费另类| 亚洲精品久久成人aⅴ小说| 亚洲av成人不卡在线观看播放网| 波多野结衣巨乳人妻| 亚洲色图综合在线观看| 级片在线观看| 国产精品电影一区二区三区| 亚洲最大成人中文| 日本免费一区二区三区高清不卡 | 成人国语在线视频| 国产精品久久久久久人妻精品电影| 国产一区二区在线av高清观看| 一a级毛片在线观看| 午夜福利,免费看| 欧美日韩亚洲综合一区二区三区_| 亚洲国产欧美一区二区综合| 视频在线观看一区二区三区| 欧美人与性动交α欧美精品济南到| 国产欧美日韩一区二区精品| 欧美人与性动交α欧美精品济南到| 婷婷丁香在线五月| 国内精品久久久久精免费| 欧美国产精品va在线观看不卡| 久久精品91蜜桃| 涩涩av久久男人的天堂| 欧美黄色片欧美黄色片| cao死你这个sao货| 久久香蕉激情| 精品高清国产在线一区| 黄色 视频免费看| 精品久久久久久久毛片微露脸| 国产97色在线日韩免费| 精品福利观看| 欧美日韩亚洲综合一区二区三区_| 精品欧美一区二区三区在线| 黑人巨大精品欧美一区二区mp4| 久久婷婷成人综合色麻豆| 亚洲国产精品成人综合色| 亚洲国产中文字幕在线视频| 精品一区二区三区视频在线观看免费| 亚洲av成人av| 亚洲国产中文字幕在线视频| 亚洲国产欧美日韩在线播放| 亚洲一码二码三码区别大吗| 国产麻豆69| 一本大道久久a久久精品| 丁香欧美五月| 国产成人精品在线电影| 啪啪无遮挡十八禁网站| 亚洲av五月六月丁香网| 久久精品成人免费网站| 欧美日韩乱码在线| 日本撒尿小便嘘嘘汇集6| 亚洲精品av麻豆狂野| 精品一区二区三区视频在线观看免费| 亚洲精品在线观看二区| 一进一出抽搐gif免费好疼| 老司机靠b影院| 日韩三级视频一区二区三区| 女性生殖器流出的白浆| 精品久久久久久久毛片微露脸| 久久精品aⅴ一区二区三区四区| 麻豆久久精品国产亚洲av| 久久人人97超碰香蕉20202| 少妇粗大呻吟视频| 亚洲全国av大片| 久久狼人影院| 国产精品日韩av在线免费观看 | 久久精品国产亚洲av高清一级| 天天添夜夜摸| 如日韩欧美国产精品一区二区三区| 一进一出抽搐动态| 免费一级毛片在线播放高清视频 | 国产又色又爽无遮挡免费看| 亚洲无线在线观看| 最近最新免费中文字幕在线| 波多野结衣巨乳人妻| 亚洲成人久久性| 99热只有精品国产| 制服丝袜大香蕉在线| 真人一进一出gif抽搐免费| 日韩大码丰满熟妇| 真人一进一出gif抽搐免费| 日本在线视频免费播放| 在线观看舔阴道视频| 一边摸一边做爽爽视频免费| 成人av一区二区三区在线看| 看黄色毛片网站| 久久婷婷成人综合色麻豆| 精品一品国产午夜福利视频| 国产精品精品国产色婷婷| 18禁裸乳无遮挡免费网站照片 | 99在线人妻在线中文字幕| 一级片免费观看大全| 日韩有码中文字幕| 最近最新免费中文字幕在线| 美女高潮喷水抽搐中文字幕| 久久精品成人免费网站| 99久久精品国产亚洲精品| 国产日韩一区二区三区精品不卡| 国产视频一区二区在线看| 国产免费av片在线观看野外av| 久久精品国产亚洲av香蕉五月| 在线天堂中文资源库| 亚洲精华国产精华精| 亚洲五月婷婷丁香| 天堂动漫精品| 国产亚洲精品久久久久5区| 精品欧美国产一区二区三| 香蕉久久夜色| 国产男靠女视频免费网站| 香蕉丝袜av| 别揉我奶头~嗯~啊~动态视频| 国产高清videossex| 一边摸一边抽搐一进一出视频| 亚洲国产欧美日韩在线播放| 亚洲精品在线观看二区| 亚洲久久久国产精品| 午夜福利在线观看吧| 国产成人av教育| 午夜日韩欧美国产| 国产高清视频在线播放一区| av欧美777| 99精品欧美一区二区三区四区| 国产精品av久久久久免费| bbb黄色大片| 国产亚洲精品久久久久5区| 成人亚洲精品一区在线观看| 久久久精品国产亚洲av高清涩受| 中文字幕久久专区| 女人高潮潮喷娇喘18禁视频| 视频区欧美日本亚洲| 非洲黑人性xxxx精品又粗又长| 午夜日韩欧美国产| 少妇粗大呻吟视频| 国产一卡二卡三卡精品| 亚洲人成电影观看| 亚洲免费av在线视频| 91大片在线观看| 亚洲午夜理论影院| 好男人电影高清在线观看| 国产真人三级小视频在线观看| 国产精品香港三级国产av潘金莲| 大香蕉久久成人网| av天堂久久9| 精品久久久精品久久久| 国产av在哪里看| 男女做爰动态图高潮gif福利片 | 精品久久久久久久久久免费视频| 麻豆成人av在线观看| 午夜福利一区二区在线看| 丝袜美腿诱惑在线| 午夜老司机福利片| 咕卡用的链子| 天天躁狠狠躁夜夜躁狠狠躁| 一个人免费在线观看的高清视频| 69精品国产乱码久久久| 窝窝影院91人妻| 99久久久亚洲精品蜜臀av| АⅤ资源中文在线天堂| 天天添夜夜摸| 人人妻人人澡人人看| 成人欧美大片| 国产午夜福利久久久久久| 免费高清视频大片| 男人的好看免费观看在线视频 | 欧美国产精品va在线观看不卡| 久久亚洲精品不卡| 18美女黄网站色大片免费观看| 亚洲欧洲精品一区二区精品久久久| 动漫黄色视频在线观看| 免费看美女性在线毛片视频| 亚洲男人的天堂狠狠| 夜夜看夜夜爽夜夜摸| 精品人妻1区二区| 大型av网站在线播放| www.熟女人妻精品国产| 妹子高潮喷水视频| 琪琪午夜伦伦电影理论片6080| 日本撒尿小便嘘嘘汇集6| 免费在线观看日本一区| 不卡av一区二区三区| 久久人人97超碰香蕉20202| 欧美午夜高清在线| 麻豆av在线久日| 国产日韩一区二区三区精品不卡| 一边摸一边抽搐一进一出视频| 黑人欧美特级aaaaaa片| 欧美日本视频| 男人舔女人下体高潮全视频| 日韩av在线大香蕉| 欧美黄色淫秽网站| 久久精品亚洲精品国产色婷小说| 久久午夜综合久久蜜桃| 日本黄色视频三级网站网址| 欧美黑人欧美精品刺激| 一进一出抽搐动态| 午夜视频精品福利| 国产精品精品国产色婷婷| 亚洲国产精品久久男人天堂| 午夜福利高清视频| 国产精品久久视频播放| 亚洲一区中文字幕在线| 国产精品精品国产色婷婷| 日日爽夜夜爽网站| 在线天堂中文资源库| 精品久久久久久成人av| 啦啦啦免费观看视频1| www国产在线视频色| 亚洲av成人不卡在线观看播放网| 在线观看日韩欧美| 亚洲精品国产区一区二| 亚洲成av人片免费观看| 亚洲五月色婷婷综合| 男女下面插进去视频免费观看| 国产亚洲精品av在线| 人人妻人人爽人人添夜夜欢视频| 在线av久久热| 亚洲av成人一区二区三| 精品电影一区二区在线| 国产精品免费一区二区三区在线| 在线观看日韩欧美| 久久精品亚洲精品国产色婷小说| aaaaa片日本免费| 一级毛片女人18水好多| 高潮久久久久久久久久久不卡| 欧美日本视频| 男女下面进入的视频免费午夜 | 精品一区二区三区视频在线观看免费| 亚洲av日韩精品久久久久久密| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人看的免费小视频| 久久久久久人人人人人| videosex国产| 性少妇av在线| 岛国在线观看网站| 黄色 视频免费看| 欧美日韩乱码在线| 亚洲男人的天堂狠狠| 少妇粗大呻吟视频| 满18在线观看网站| 欧美久久黑人一区二区| 无限看片的www在线观看| 日韩精品免费视频一区二区三区| 国产麻豆69| 国产免费男女视频| 无遮挡黄片免费观看| 亚洲精品美女久久av网站| www国产在线视频色| 国产亚洲av高清不卡| 国产一区二区三区在线臀色熟女| 国产av一区在线观看免费| 国产精品综合久久久久久久免费 | 午夜福利,免费看| 亚洲精品一区av在线观看| 国产伦人伦偷精品视频| 国产av一区在线观看免费| 多毛熟女@视频| 国产精品野战在线观看| 精品电影一区二区在线| 日日干狠狠操夜夜爽| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 极品教师在线免费播放| 国内久久婷婷六月综合欲色啪| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av第一区精品v没综合| 亚洲五月天丁香| 国内精品久久久久久久电影| 中文字幕精品免费在线观看视频| 一二三四社区在线视频社区8| www国产在线视频色| 久久久久国产精品人妻aⅴ院| 一区福利在线观看| 国产亚洲欧美在线一区二区| 老司机福利观看| 国产欧美日韩一区二区三区在线| 欧美日韩亚洲综合一区二区三区_| www.精华液| 午夜视频精品福利| 亚洲av熟女| 久久影院123| 婷婷精品国产亚洲av在线| 男女下面插进去视频免费观看| av片东京热男人的天堂| 日韩精品中文字幕看吧| 国产一卡二卡三卡精品| 日本vs欧美在线观看视频| 成人永久免费在线观看视频| 一进一出抽搐gif免费好疼| 国产亚洲精品久久久久久毛片| 欧美激情久久久久久爽电影 | 日本五十路高清| 国产精品久久久人人做人人爽| 国产精品九九99| 黄片大片在线免费观看| 亚洲av成人一区二区三| 国产私拍福利视频在线观看| 波多野结衣一区麻豆| 一本大道久久a久久精品| 天天躁狠狠躁夜夜躁狠狠躁| 两性夫妻黄色片| 手机成人av网站| 成在线人永久免费视频| 中文字幕人妻熟女乱码| 亚洲va日本ⅴa欧美va伊人久久| 精品乱码久久久久久99久播| 中文字幕人妻熟女乱码| 久热爱精品视频在线9| 午夜福利免费观看在线| 国产av一区在线观看免费| 亚洲精品一区av在线观看| 一本久久中文字幕| 欧美 亚洲 国产 日韩一| 长腿黑丝高跟| 法律面前人人平等表现在哪些方面| 国产精品影院久久| 日韩大码丰满熟妇| 色综合婷婷激情| 一区二区日韩欧美中文字幕| 午夜视频精品福利| 欧美精品啪啪一区二区三区| 国产主播在线观看一区二区| 久9热在线精品视频| 国产激情欧美一区二区| 欧美+亚洲+日韩+国产| 深夜精品福利| 免费观看人在逋| 美女扒开内裤让男人捅视频| 亚洲欧美激情在线| 国产精品av久久久久免费| 日本三级黄在线观看| 亚洲成人久久性| 亚洲av第一区精品v没综合| 久久久久精品国产欧美久久久| 又紧又爽又黄一区二区| 免费久久久久久久精品成人欧美视频| 美女国产高潮福利片在线看| 亚洲久久久国产精品| 成人亚洲精品av一区二区| 久久久久久人人人人人| 国产精品av久久久久免费| 精品久久蜜臀av无| 黑人操中国人逼视频| 国产精品免费一区二区三区在线| 成年女人毛片免费观看观看9| 久久欧美精品欧美久久欧美| 香蕉国产在线看| 不卡av一区二区三区| 国产蜜桃级精品一区二区三区| 一进一出抽搐动态| 妹子高潮喷水视频| 欧美成人性av电影在线观看| 操美女的视频在线观看| 一区二区三区激情视频| 亚洲av成人av| 国产亚洲欧美98| 国产极品粉嫩免费观看在线| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| 一a级毛片在线观看| 久久午夜亚洲精品久久| 欧美一级a爱片免费观看看 | 变态另类成人亚洲欧美熟女 | 99在线视频只有这里精品首页| 精品国产一区二区久久| 一边摸一边做爽爽视频免费| 免费av毛片视频| 国产麻豆69| 精品久久久久久久人妻蜜臀av | 欧美不卡视频在线免费观看 | 香蕉国产在线看| 国产色视频综合| 在线永久观看黄色视频| 人人澡人人妻人| www.999成人在线观看| 啦啦啦观看免费观看视频高清 | 国产亚洲欧美精品永久| 日本黄色视频三级网站网址| 亚洲狠狠婷婷综合久久图片| netflix在线观看网站| 欧美激情高清一区二区三区| 欧美精品啪啪一区二区三区| av超薄肉色丝袜交足视频| 日韩精品中文字幕看吧| 中文字幕另类日韩欧美亚洲嫩草| www日本在线高清视频| 夜夜爽天天搞| 国产一区二区三区视频了| 国产成人一区二区三区免费视频网站| 老司机靠b影院| 国产欧美日韩一区二区精品| 9热在线视频观看99| 69av精品久久久久久| 欧美精品亚洲一区二区| 日韩av在线大香蕉| 亚洲午夜精品一区,二区,三区| 国产成人精品无人区| 亚洲人成网站在线播放欧美日韩| 午夜福利视频1000在线观看 | 一边摸一边抽搐一进一小说| 老汉色∧v一级毛片|