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

    蛋白質(zhì)分子機器結(jié)構(gòu)的特性分析及去折疊研究

    2022-10-13 08:45:56張莉莉蔣益鋒謝良旭
    數(shù)據(jù)采集與處理 2022年5期
    關(guān)鍵詞:均方殘基高斯

    張莉莉,蔣益鋒,謝良旭,孔 韌,常 珊

    (江蘇理工學(xué)院電氣信息工程學(xué)院生物信息與醫(yī)藥工程研究所,常州 213001)

    引 言

    在生物體內(nèi),由蛋白質(zhì)組成的各種類型的分子機器驅(qū)動著各種各樣的生命活動所必需的化學(xué)反應(yīng)[1]。深入理解和分析蛋白質(zhì)的結(jié)構(gòu)特征,可以闡明蛋白質(zhì)功能,解釋蛋白質(zhì)錯誤折疊引起的相關(guān)疾病起源,以及對于藥物設(shè)計工作具有重要意義[2]。蛋白質(zhì)GB1參與許多生理信號的檢測,包括激素、神經(jīng)遞質(zhì)和各種感覺刺激(光、氣味等物質(zhì))[3]。此外,它還與疾?。ɡ珉貌《竞桶柎暮D。┫嚓P(guān)的錯誤折疊狀態(tài)的存在以及β聚集(淀粉樣疾病)的研究相關(guān)[4],因此關(guān)于GB1蛋白的研究具有重要意義。理解蛋白質(zhì)GB1結(jié)構(gòu)的折疊機制和穩(wěn)定性是治療人類疾病的重要基礎(chǔ),也有助于蛋白質(zhì)的開發(fā)設(shè)計。近年來,關(guān)于蛋白質(zhì)結(jié)構(gòu)折疊的計算機模擬研究方法主要有兩種,即分子動力學(xué)模擬和彈性網(wǎng)絡(luò)模型。分子動力學(xué)模擬方法是一種細粒度方法,可以觀察到蛋白質(zhì)的折疊路徑、過渡態(tài)等,但是該方法計算復(fù)雜、耗時較長,就目前計算機的模擬水平,僅僅只對一些小蛋白質(zhì)分子的折疊結(jié)構(gòu)模擬效果較好[5]。而彈性網(wǎng)絡(luò)模型關(guān)鍵是給出適合簡化模型的勢函數(shù),計算簡單、耗時短,可模擬時間跨度大的去折疊過程,相對分子動力學(xué)模擬方法而言效率較高[6-7],能夠很好地再現(xiàn)蛋白質(zhì)的低頻運動(長時間動力學(xué)),提供關(guān)于它們的平衡動力學(xué)、天然結(jié)構(gòu)拓撲對它們穩(wěn)定性的影響、蛋白質(zhì)波動的定位特性或蛋白質(zhì)結(jié)構(gòu)域的定義的信息[8]。彈性網(wǎng)絡(luò)模型(Elastic network model,ENM)在蛋白質(zhì)結(jié)構(gòu)-功能關(guān)系研究中得到了廣泛應(yīng)用。先前有研究通過應(yīng)用ENM來評估生物分子整體編碼、蛋白質(zhì)功能性運動分析和關(guān)鍵位點識別等[9],此外,還有研究結(jié)果表明,應(yīng)用ENM有助于更好地理解轉(zhuǎn)運體系發(fā)揮生物學(xué)功能的分子機制。經(jīng)典的彈性網(wǎng)絡(luò)模型能夠提供蛋白質(zhì)在平衡態(tài)(通常為原生態(tài))附近的動態(tài)特性,因此它們被廣泛應(yīng)用于許多蛋白質(zhì)的系統(tǒng)比較。然而,蛋白質(zhì)折疊通常遠離平衡態(tài),所以一般的ENM不適合蛋白質(zhì)折疊的研究。蛋白質(zhì)結(jié)構(gòu)預(yù)測一直是生命科學(xué)里的一個重要問題,研究蛋白質(zhì)序列和結(jié)構(gòu)間關(guān)系的蛋白質(zhì)折疊問題是生物物理領(lǐng)域最重要的基礎(chǔ)問題之一。在2020年舉辦的第14屆蛋白質(zhì)結(jié)構(gòu)預(yù)測競賽CASP14(Critical assessment of protein structure prediction)中,Google DeepMind團隊使用AlphaFold2預(yù)測了多個物種中共30余萬個無實驗結(jié)構(gòu)的蛋白質(zhì)的結(jié)構(gòu)模型,并聯(lián)手EBI建立了結(jié)構(gòu)預(yù)測數(shù)據(jù)庫AFDB[10]。這一系列成果的出現(xiàn)吸引了科學(xué)界的大量關(guān)注。AlphaFold2等結(jié)構(gòu)預(yù)測方法目前僅能預(yù)測特定氨基酸序列的靜態(tài)構(gòu)象。蛋白質(zhì)在行使生物學(xué)功能時往往需要發(fā)生構(gòu)象變化。比如酶從失活狀態(tài)轉(zhuǎn)變?yōu)榛钚誀顟B(tài)、膜轉(zhuǎn)運蛋白需要通過構(gòu)象變化交替接觸膜兩側(cè)的溶液、蛋白和配體結(jié)合時發(fā)生構(gòu)象變化等等。高斯網(wǎng)絡(luò)模型(Gaussian network model,GNM)是經(jīng)典ENM方法的發(fā)展,是一種基于拓撲的、不依賴序列特異性的粗粒度模型。高斯網(wǎng)絡(luò)模型可以從晶體結(jié)構(gòu)提供蛋白質(zhì)構(gòu)象轉(zhuǎn)變的信息,不需要分子動力學(xué)模擬的高計算成本,是一種基于正態(tài)模式計算的迭代方法,被提出來用于研究蛋白質(zhì)折疊/去折疊過程。多年來,GB1在蛋白質(zhì)折疊的計算和實驗研究中被廣泛用作模型系統(tǒng)[11]。本文主要就是通過利用彈性網(wǎng)絡(luò)模型,模擬GB1蛋白結(jié)構(gòu)的展開過程,再現(xiàn)GB1的快運動與慢運動模式,同時研究它的拓撲結(jié)構(gòu)對自身穩(wěn)定性的影響。

    1 蛋白質(zhì)結(jié)構(gòu)

    本研究選擇分析的蛋白質(zhì)GB1(PDB代碼:6CHE)如圖1所示[12]。GB1是一種小球狀蛋白,由β折疊和α螺旋組成,共有56個殘基。8個殘基與W43形成天然接觸:其中4個殘基(F52、T53、V54和T55)位于相鄰的β折疊中,并與W43的骨架形成天然接觸,而其他4個殘基(L5、F30、K31和M34)與W43的側(cè)鏈相互 作 用。在GB1中,殘 基2-19形 成N端β折 疊,殘 基23-36形成α螺旋,殘基42-55形成C端β折疊[13]。

    2 本文方法

    2.1 高斯網(wǎng)絡(luò)模型

    在高斯網(wǎng)絡(luò)模型中,每個蛋白質(zhì)的三維結(jié)構(gòu)可以簡化為一個彈性網(wǎng)絡(luò),其中每個氨基酸(殘基)被看作為該網(wǎng)絡(luò)中的頂點,如果兩個頂點間距離小于截止距離,則用一根彈簧將其連接,所有彈簧的彈性系數(shù)都相同[14]?;谠摼W(wǎng)絡(luò)模型,網(wǎng)絡(luò)的總能量可以寫成

    式中:V為所有接觸殘基的總能量;γ為彈性系數(shù);{ΔR}為殘基漲落的N維列向量;Γ為N階對稱矩陣,在對稱矩陣中的元素可寫為

    式中:Rij為蛋白質(zhì)中第i個和第j個殘基之間的距離;Γc是截止距離(本研究中采用的截止距離為7.4?)。

    N階對稱矩陣Γ的逆矩陣可表示為

    式中:U為正交矩陣,其列向量Ui(1≤i≤N)是Γ的特征向量;Λ為對角矩陣,其對角線上的元素是Γ的特征值。

    蛋白質(zhì)中兩個殘基均方漲落的互相關(guān)性計算可表示為

    式中:i和j分別表示蛋白質(zhì)中第i個和第j個殘基;kB為玻爾茲曼常數(shù);T為絕對溫度。當(dāng)i=j時,第i個殘基的均方漲落計算式可表示為

    根據(jù)Debye-Waller理論,第i個殘基的B因子計算式可表示為

    在高斯網(wǎng)絡(luò)模型中,歸一化的互相關(guān)性系數(shù)可寫成[15]

    高斯網(wǎng)絡(luò)模型是建立在多聚體網(wǎng)絡(luò)的波動動力學(xué)基礎(chǔ)之上的,彈性網(wǎng)絡(luò)模型可以是原子層次上的粗?;P?,也可以是殘基層次上的粗?;P?。高斯網(wǎng)絡(luò)模型模擬方法可以把蛋白質(zhì)的功能性運動分解成為各個不同種運動模式的疊加,在不同種運動模式中,慢運動模式為對應(yīng)著與蛋白質(zhì)功能相關(guān)的大幅度集合運動[16]。通過與實驗數(shù)據(jù)的對比可以發(fā)現(xiàn)這種方法所得到的數(shù)據(jù)結(jié)果是可靠且有效的。

    2.2 去折疊原理

    為了研究蛋白質(zhì)的去折疊過程,本文提出了一種基于高斯網(wǎng)絡(luò)模型的迭代方法。所有殘基對之間距離的均方漲落都是基于高斯網(wǎng)絡(luò)模型計算的,第i個殘基和第j個殘基之間距離的均方漲落可表示為[17]

    式中:Rij和分別為殘基i和殘基j之間的瞬時和平衡分離向量。

    蛋白質(zhì)結(jié)構(gòu)去折疊過程的模擬方案如下:

    (1)基于式(8)和蛋白質(zhì)的天然拓撲結(jié)構(gòu)計算出結(jié)構(gòu)中所有殘基對之間距離的均方漲落值;

    (2)斷開距離均方漲落值最大的殘基對之間的天然接觸,得到對應(yīng)新Γ矩陣的結(jié)構(gòu)拓撲;

    (3)基于新的Γ矩陣,利用式(8)重新計算所有殘基對之間距離的均方漲落值;

    (4)重復(fù)上述兩個步驟,直到蛋白質(zhì)中所有的非共價接觸被斷開;

    (5)綜合由以上步驟得到的所有結(jié)構(gòu)拓撲信息,以獲取蛋白質(zhì)的去折疊過程。

    3 結(jié)果與討論

    3.1 實驗與模擬所得的B因子分析

    為了評價高斯網(wǎng)絡(luò)模型方法在本研究中應(yīng)用的可行性,計算了B因子,并與X射線(X-RAY)實驗數(shù)據(jù)對比。根據(jù)GNM模擬所得的數(shù)據(jù)與X-RAY實驗結(jié)果對比結(jié)果如圖2所示,其中,紅色曲線對應(yīng)基于GNM模擬所得的數(shù)據(jù),綠色曲線對應(yīng)X-RAY實驗數(shù)據(jù)??梢钥闯觯瑑蓷l曲線的峰值和谷值出現(xiàn)的位置幾乎相同,模擬所得的數(shù)據(jù)與實驗數(shù)據(jù)之間的相關(guān)系數(shù)為0.70。綜合以往的文獻研究得到,一般模擬數(shù)據(jù)與實驗數(shù)據(jù)之間的相關(guān)系數(shù)取值為0.53~0.89[6,14-15,18],本 次 實 驗結(jié) 果 所得 值 為0.70,在 此 范圍內(nèi),可見該方法是適用的,表明該模型適用于研究GB1蛋白的固有動力學(xué)。

    圖2 實驗與模擬所得的B因子對比Fig.2 Comparison of B factor between experiment and simulation

    3.2 蛋白質(zhì)的快運動與慢運動模式分析

    運動的快模式對應(yīng)于局部結(jié)構(gòu)中的幾何不規(guī)則性。以前的研究發(fā)現(xiàn),高頻波動殘基被認為是動力學(xué)關(guān)鍵殘基,對三級折疊的穩(wěn)定性至關(guān)重要[19]。圖3顯示了GB1蛋白的快運動模式。圖3中,橫坐標表示殘基序號,縱坐標表示殘基自身距離的均方漲落值,基于式(5)求得,單位為平方埃(?2)。從圖3可以看出,殘基Lys4、Ala26、Thr51和Val54(圖中已標注)是曲線中的峰值。本文結(jié)果與以前的研究[20]一致,表明這些殘基在蛋白質(zhì)的穩(wěn)定性中起著關(guān)鍵作用。

    圖3 GB1快運動模式結(jié)果Fig.3 The fastest mode shapes of GB1

    在蛋白質(zhì)的研究過程中,慢運動模式代表著蛋白質(zhì)結(jié)構(gòu)中編碼的長程集體運動,同時相關(guān)研究認為那些慢運動模式就相當(dāng)于大幅度的集體運動,而大幅度集體運動往往與蛋白質(zhì)運動相關(guān)[21]。圖4顯示了基于高斯網(wǎng)絡(luò)模型計算的GB1蛋白的最慢模式。從圖4可以看出,大多數(shù)殘基波動值較高,這意味著這些結(jié)構(gòu)相對而言不是很穩(wěn)定。同時,還可以從圖中看出,殘基Gln2、Tyr3和Thr18的波動值保持較低。

    圖4 GB1慢運動模式結(jié)果圖Fig.4 The slowest mode shapes of GB1

    3.3 蛋白質(zhì)去折疊過程分析

    為了詳細說明展開模擬過程中自然接觸的損失,構(gòu)建了不同快照中構(gòu)象的接觸圖,結(jié)果如圖5所示。圖5(a)顯示了GB1蛋白天然結(jié)構(gòu)的接觸圖,即當(dāng)兩個殘基之間的距離小于7.4?時,兩個殘基被定義為相互接觸。如果兩個殘基直接有接觸,則用*表示,圖5(b~f)分別展示了GB1蛋白的非共價接觸損失數(shù)(Loss number of noncovalent contact,LNNC)分別為20、50、100、130和170的結(jié)果。圖5(a)天然狀態(tài)下的接觸呈現(xiàn)結(jié)果與之前的相關(guān)研究一致[22]。結(jié)果表明,GB1蛋白的展開有一個優(yōu)先的過程,它顯示了一系列事件。

    圖5 GB1天然結(jié)構(gòu)以及非共價接觸損失數(shù)分別為20、50、100、130和170的接觸圖Fig.5 Contact maps of native conformation and conformations with LNNC of 20,50,100,130,170 for GB1

    由圖5(a)可以看出,在GB1的天然結(jié)構(gòu)中,碳末端折疊比氮末端折疊有更多更強的接觸(圖1),這可能導(dǎo)致碳末端區(qū)域更快的折疊。從圖5(b,c)的實驗結(jié)果可以看出,隨著殘基對之間非共價接觸損失個數(shù)的增加,GB1蛋白一開始主要是從β2折疊部分的殘基對之間的接觸先斷開,此外,從圖5(c)也可以看出β4在非共價接觸損失數(shù)為50左右的時候開始斷開了。繼而如圖5(d,e),α螺旋部分的殘基對之間的接觸再逐漸斷開,直至如圖5(f),最終幾乎所有接觸斷開,即GB1蛋白完全展開。該過程顯示了GB1蛋白的展開是從大量的α螺旋和β2折疊結(jié)構(gòu)元素的接觸損失開始,同時先保持了大部分其他β結(jié)構(gòu)的完整。本模擬結(jié)果與之前的實驗研究結(jié)果一致[13]。

    此外,折疊協(xié)同性被認為是蛋白質(zhì)折疊動力學(xué)的一個重要行為[23]。在本研究模型中,展開路徑是連續(xù)的,很難直接觀察展開過程中的協(xié)同性。事實上,這些高度合作的行為發(fā)生在這個迭代展開模型的近鄰步驟中。結(jié)果表明,解折疊路徑主要由其自身的拓撲結(jié)構(gòu)決定,迭代解折疊方法可以合理地描述GB1的去折疊過程。

    3.4 蛋白質(zhì)去折疊過程中的互相關(guān)分析

    此外,本文還研究了在GB1蛋白去折疊過程中殘基波動之間的相關(guān)性的變化。殘基波動之間的互相關(guān)用式(7)計算?;ハ嚓P(guān)值的取值范圍為-1到1。其中,正值表示殘基間運動方向相同,負值則表示它們之間運動方向相反。絕對互相關(guān)值越高,兩個殘基越相關(guān)(或反相關(guān))。另外,互相關(guān)值0意味著殘基的運動完全不相關(guān)[14]。圖6顯示了GB1蛋白的互相關(guān)圖。

    圖6 GB1天然結(jié)構(gòu)以及展開過程中非共價接觸損失數(shù)分別為20、50、100、130和170時的殘基互相關(guān)圖Fig.6 Cross-correlation maps calculated using all modes for native conformation and conformations with LNNC of 20,50,100,130,170 during the unfolding process of GB1

    如圖6(a)所示,沿著圖的對角線,有一些正相關(guān)的光塊,對應(yīng)α螺旋和β折疊的二級結(jié)構(gòu)。隨著殘基對之間非共價損失個數(shù)的增加,即隨著GB1蛋白的逐漸展開,如圖6(b,c),當(dāng)α螺旋和β折疊中的天然觸點開始丟失時,α螺旋和β折疊之間負相關(guān),β折疊之間的正相關(guān)性提高;隨著天然觸點丟失個數(shù)的增加,如圖6(d,e),α螺旋和β折疊僅部分保留,最后,如圖6(f),蛋白質(zhì)的結(jié)構(gòu)似乎被分成兩個方向相反的方向波動。該圖反映的是去折疊的最后狀態(tài),即蛋白質(zhì)結(jié)構(gòu)展開回到了最初未折疊的多肽鏈結(jié)構(gòu)。根據(jù)先前的研究發(fā)現(xiàn)[14],當(dāng)去折疊模擬到最后的階段時,蛋白質(zhì)的結(jié)構(gòu)也似乎被分為兩部分,上下波動方向相反,與本次實驗結(jié)果一致。

    4 結(jié)束語

    本研究基于GB1的拓撲結(jié)構(gòu),采用高斯網(wǎng)絡(luò)模型模擬了GB1的快運動與慢運動模式,并對其做了相應(yīng)的結(jié)果分析;同時,對其拓撲結(jié)構(gòu)做了展開過程的路徑研究;此外,還研究了GB1蛋白在去折疊過程中殘基波動之間相關(guān)性的變化。與相關(guān)實驗和分子動力學(xué)模擬數(shù)據(jù)吻合良好,表明彈性網(wǎng)絡(luò)模型的計算效率高,能夠準確模擬蛋白質(zhì)的動態(tài)和結(jié)構(gòu)特性,能夠很好地再現(xiàn)蛋白質(zhì)的運動特性,提供關(guān)于它們的平衡動力學(xué)、天然結(jié)構(gòu)拓撲對其穩(wěn)定性的影響、蛋白質(zhì)波動的定位特性或蛋白質(zhì)結(jié)構(gòu)域的信息,適用于對蛋白質(zhì)的研究。

    猜你喜歡
    均方殘基高斯
    小高斯的大發(fā)現(xiàn)
    一類隨機積分微分方程的均方漸近概周期解
    基于各向異性網(wǎng)絡(luò)模型研究δ阿片受體的動力學(xué)與關(guān)鍵殘基*
    “殘基片段和排列組合法”在書寫限制條件的同分異構(gòu)體中的應(yīng)用
    Beidou, le système de navigation par satellite compatible et interopérable
    天才數(shù)學(xué)家——高斯
    蛋白質(zhì)二級結(jié)構(gòu)序列與殘基種類間關(guān)聯(lián)的分析
    基于支持向量機的蛋白質(zhì)相互作用界面熱點殘基預(yù)測
    基于抗差最小均方估計的輸電線路參數(shù)辨識
    基于隨機牽制控制的復(fù)雜網(wǎng)絡(luò)均方簇同步
    精品免费久久久久久久清纯| 在线国产一区二区在线| 午夜视频精品福利| 国产精品,欧美在线| 搡老妇女老女人老熟妇| 亚洲第一欧美日韩一区二区三区| 特级一级黄色大片| 久久久国产欧美日韩av| 91av网站免费观看| 最新在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看 | 村上凉子中文字幕在线| 亚洲九九香蕉| xxxwww97欧美| 国产精品99久久久久久久久| 美女大奶头视频| 精品久久久久久久毛片微露脸| 亚洲国产色片| 免费一级毛片在线播放高清视频| 日韩精品青青久久久久久| 亚洲熟妇中文字幕五十中出| 亚洲五月婷婷丁香| 久久人妻av系列| 亚洲色图av天堂| 亚洲欧美日韩无卡精品| 我要搜黄色片| 久久精品国产99精品国产亚洲性色| 精品国内亚洲2022精品成人| 每晚都被弄得嗷嗷叫到高潮| 久久久水蜜桃国产精品网| 美女cb高潮喷水在线观看 | 老汉色av国产亚洲站长工具| 熟女人妻精品中文字幕| 免费一级毛片在线播放高清视频| 日本五十路高清| 欧美性猛交黑人性爽| 日韩有码中文字幕| 脱女人内裤的视频| www.自偷自拍.com| 精品国产乱码久久久久久男人| 在线观看日韩欧美| 精品久久久久久久久久免费视频| 国产精品99久久久久久久久| 国产精品99久久久久久久久| av在线天堂中文字幕| 久久精品综合一区二区三区| 久久精品影院6| 长腿黑丝高跟| 免费人成视频x8x8入口观看| 一边摸一边抽搐一进一小说| 麻豆国产av国片精品| 在线观看免费午夜福利视频| 成熟少妇高潮喷水视频| 国产成人av激情在线播放| 国产亚洲精品综合一区在线观看| 免费av不卡在线播放| 老汉色∧v一级毛片| 欧美+亚洲+日韩+国产| 午夜福利18| 国产91精品成人一区二区三区| АⅤ资源中文在线天堂| 亚洲av电影在线进入| 亚洲欧美日韩卡通动漫| 美女cb高潮喷水在线观看 | 国产精品99久久99久久久不卡| 国产一区二区三区视频了| 精品人妻1区二区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲熟妇熟女久久| www国产在线视频色| 淫秽高清视频在线观看| а√天堂www在线а√下载| 99riav亚洲国产免费| 两性夫妻黄色片| 久久久久性生活片| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久久久免费视频| 中文资源天堂在线| 亚洲午夜精品一区,二区,三区| 999精品在线视频| 日韩成人在线观看一区二区三区| АⅤ资源中文在线天堂| 色播亚洲综合网| 久9热在线精品视频| 精品午夜福利视频在线观看一区| 天堂网av新在线| 国产高清激情床上av| 亚洲精品粉嫩美女一区| 可以在线观看的亚洲视频| 少妇的逼水好多| 亚洲专区中文字幕在线| 免费看美女性在线毛片视频| 日韩欧美在线二视频| 在线a可以看的网站| 国产欧美日韩一区二区三| 又大又爽又粗| 亚洲av中文字字幕乱码综合| 欧美av亚洲av综合av国产av| 成人无遮挡网站| 久久午夜综合久久蜜桃| 97碰自拍视频| 91久久精品国产一区二区成人 | 欧美丝袜亚洲另类 | 69av精品久久久久久| 法律面前人人平等表现在哪些方面| 免费av不卡在线播放| 91麻豆精品激情在线观看国产| 男人的好看免费观看在线视频| 久久久精品大字幕| 亚洲 国产 在线| 久久久精品欧美日韩精品| xxxwww97欧美| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 久久久国产成人精品二区| 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 一本一本综合久久| 国产精品99久久久久久久久| 丰满人妻一区二区三区视频av | 色视频www国产| 不卡av一区二区三区| 久久精品aⅴ一区二区三区四区| 这个男人来自地球电影免费观看| 国产伦精品一区二区三区四那| 精品人妻1区二区| 最近最新中文字幕大全免费视频| 亚洲人成网站高清观看| 色在线成人网| 性色avwww在线观看| 中文在线观看免费www的网站| 国产熟女xx| av在线蜜桃| 91麻豆精品激情在线观看国产| 国产一区在线观看成人免费| 丁香欧美五月| 午夜福利免费观看在线| 国产亚洲av嫩草精品影院| 成年女人毛片免费观看观看9| 亚洲人成网站高清观看| 手机成人av网站| 亚洲av成人不卡在线观看播放网| 一进一出好大好爽视频| 一个人免费在线观看电影 | 国产成+人综合+亚洲专区| 两个人的视频大全免费| 久久久成人免费电影| 综合色av麻豆| 毛片女人毛片| 国产精品 欧美亚洲| 欧美精品啪啪一区二区三区| 国产男靠女视频免费网站| 国产午夜精品久久久久久| 欧美成人一区二区免费高清观看 | 热99在线观看视频| 亚洲欧美一区二区三区黑人| 国产精品久久视频播放| 成人特级av手机在线观看| 欧美一区二区国产精品久久精品| 国产成人aa在线观看| 欧美最黄视频在线播放免费| 久久国产精品影院| 在线免费观看的www视频| 亚洲国产色片| 中文字幕高清在线视频| 欧美+亚洲+日韩+国产| 亚洲18禁久久av| 免费观看的影片在线观看| 国产在线精品亚洲第一网站| 精品一区二区三区视频在线观看免费| 黄片小视频在线播放| 麻豆av在线久日| 制服人妻中文乱码| 国产午夜福利久久久久久| 国产精品一区二区三区四区久久| 久久久久九九精品影院| 欧美国产日韩亚洲一区| 亚洲av成人一区二区三| 精品电影一区二区在线| 免费高清视频大片| 精品一区二区三区视频在线 | 在线播放国产精品三级| 国产精品1区2区在线观看.| 少妇人妻一区二区三区视频| 久久中文字幕一级| 91麻豆av在线| 精品无人区乱码1区二区| 中出人妻视频一区二区| 特级一级黄色大片| 国产三级在线视频| 国产精品爽爽va在线观看网站| 成人av在线播放网站| 久久午夜综合久久蜜桃| 成人特级黄色片久久久久久久| 久久婷婷人人爽人人干人人爱| 亚洲黑人精品在线| 无遮挡黄片免费观看| 欧美激情久久久久久爽电影| av国产免费在线观看| 精品国产乱码久久久久久男人| 桃色一区二区三区在线观看| 搞女人的毛片| 久久精品国产亚洲av香蕉五月| 十八禁人妻一区二区| 最近最新中文字幕大全电影3| 亚洲欧美日韩卡通动漫| 亚洲五月天丁香| 亚洲精华国产精华精| 欧美日韩亚洲国产一区二区在线观看| 无限看片的www在线观看| 久久热在线av| 日本成人三级电影网站| 一a级毛片在线观看| 999久久久精品免费观看国产| 人妻夜夜爽99麻豆av| 天堂√8在线中文| 在线看三级毛片| 国产午夜精品论理片| 无遮挡黄片免费观看| 淫秽高清视频在线观看| 久久久久久九九精品二区国产| 亚洲18禁久久av| 日本与韩国留学比较| 色噜噜av男人的天堂激情| 久久中文字幕人妻熟女| 99久久国产精品久久久| 中文在线观看免费www的网站| 在线观看免费视频日本深夜| 国产又色又爽无遮挡免费看| 亚洲aⅴ乱码一区二区在线播放| 老司机福利观看| 又黄又爽又免费观看的视频| 欧美在线黄色| 成在线人永久免费视频| 国产精品亚洲av一区麻豆| 又黄又粗又硬又大视频| 少妇人妻一区二区三区视频| 岛国在线免费视频观看| 老司机福利观看| 国产黄片美女视频| 宅男免费午夜| 欧美乱妇无乱码| 精品国内亚洲2022精品成人| 久久精品人妻少妇| 叶爱在线成人免费视频播放| 亚洲欧美精品综合久久99| 老司机在亚洲福利影院| 国产精品免费一区二区三区在线| 中国美女看黄片| 国产久久久一区二区三区| 亚洲色图av天堂| 国产精品一及| 国产亚洲av嫩草精品影院| 国产精品野战在线观看| 国产野战对白在线观看| 亚洲av熟女| 亚洲第一欧美日韩一区二区三区| 色噜噜av男人的天堂激情| 最近最新中文字幕大全免费视频| 麻豆av在线久日| 波多野结衣巨乳人妻| 两个人的视频大全免费| 噜噜噜噜噜久久久久久91| 婷婷丁香在线五月| 老熟妇乱子伦视频在线观看| 在线免费观看的www视频| 非洲黑人性xxxx精品又粗又长| 成人高潮视频无遮挡免费网站| 亚洲,欧美精品.| 夜夜夜夜夜久久久久| 九九在线视频观看精品| 国产精品九九99| 性欧美人与动物交配| 国产精品国产高清国产av| 国产精品一区二区三区四区久久| 老司机午夜十八禁免费视频| 欧美日韩瑟瑟在线播放| 99久久精品国产亚洲精品| 真人做人爱边吃奶动态| 少妇人妻一区二区三区视频| 日韩欧美三级三区| 久久香蕉精品热| 国产黄a三级三级三级人| 老司机午夜十八禁免费视频| 91在线观看av| www.自偷自拍.com| 国产一区二区三区在线臀色熟女| 非洲黑人性xxxx精品又粗又长| netflix在线观看网站| 亚洲真实伦在线观看| 18禁美女被吸乳视频| 国产精品永久免费网站| 午夜福利高清视频| 香蕉丝袜av| 日韩欧美在线二视频| 久久久久久人人人人人| 偷拍熟女少妇极品色| 成人三级做爰电影| 身体一侧抽搐| 国内毛片毛片毛片毛片毛片| 蜜桃久久精品国产亚洲av| 99久久久亚洲精品蜜臀av| 一夜夜www| 亚洲av成人精品一区久久| 亚洲真实伦在线观看| 国产精品av久久久久免费| 黄片小视频在线播放| 噜噜噜噜噜久久久久久91| 精品99又大又爽又粗少妇毛片 | 精品久久蜜臀av无| 欧美日韩综合久久久久久 | 国产91精品成人一区二区三区| 女人被狂操c到高潮| 悠悠久久av| 岛国在线免费视频观看| 成人欧美大片| 在线免费观看的www视频| 欧美中文日本在线观看视频| 精品电影一区二区在线| 欧美黑人巨大hd| 午夜日韩欧美国产| 色吧在线观看| 国产成人欧美在线观看| 亚洲中文字幕一区二区三区有码在线看 | 我要搜黄色片| 欧美一级a爱片免费观看看| 国产又色又爽无遮挡免费看| x7x7x7水蜜桃| 身体一侧抽搐| 岛国在线免费视频观看| 成人亚洲精品av一区二区| 国产精品综合久久久久久久免费| 两个人视频免费观看高清| 国产精品电影一区二区三区| 久久精品aⅴ一区二区三区四区| 色视频www国产| 国产精品综合久久久久久久免费| 国产成人精品无人区| 精品久久久久久久人妻蜜臀av| 久久这里只有精品中国| 国产真实乱freesex| 欧美日本视频| 美女午夜性视频免费| 麻豆国产av国片精品| 亚洲国产欧美网| 国产高清视频在线观看网站| 97超级碰碰碰精品色视频在线观看| 全区人妻精品视频| 亚洲av第一区精品v没综合| 国产激情久久老熟女| 午夜影院日韩av| 国产成人aa在线观看| 三级男女做爰猛烈吃奶摸视频| 久久中文字幕人妻熟女| 99久久成人亚洲精品观看| 此物有八面人人有两片| 舔av片在线| 欧美大码av| 亚洲乱码一区二区免费版| 亚洲色图 男人天堂 中文字幕| 欧美不卡视频在线免费观看| АⅤ资源中文在线天堂| 国产激情久久老熟女| av天堂在线播放| 青草久久国产| 亚洲va日本ⅴa欧美va伊人久久| 一进一出好大好爽视频| 99在线视频只有这里精品首页| 婷婷亚洲欧美| 1024香蕉在线观看| 国产美女午夜福利| 香蕉久久夜色| 国内揄拍国产精品人妻在线| 午夜精品在线福利| 亚洲av片天天在线观看| 亚洲精品美女久久av网站| 别揉我奶头~嗯~啊~动态视频| 好看av亚洲va欧美ⅴa在| 国产精品一区二区三区四区久久| 久久伊人香网站| 女人高潮潮喷娇喘18禁视频| 久久性视频一级片| 一级毛片女人18水好多| 亚洲成人精品中文字幕电影| 亚洲黑人精品在线| 韩国av一区二区三区四区| 欧美黑人欧美精品刺激| 色尼玛亚洲综合影院| 亚洲成人精品中文字幕电影| 久久这里只有精品19| 精品国产超薄肉色丝袜足j| 中文字幕熟女人妻在线| 国产真实乱freesex| 国产黄片美女视频| 国产精品自产拍在线观看55亚洲| 最新美女视频免费是黄的| av国产免费在线观看| 在线看三级毛片| 久久久久久人人人人人| 18禁裸乳无遮挡免费网站照片| 不卡av一区二区三区| 国产精品日韩av在线免费观看| 最新美女视频免费是黄的| 美女午夜性视频免费| 97超级碰碰碰精品色视频在线观看| 午夜福利在线观看免费完整高清在 | 亚洲无线在线观看| 精品一区二区三区av网在线观看| 国产高清激情床上av| 国产单亲对白刺激| 成人国产一区最新在线观看| 久久久色成人| 变态另类成人亚洲欧美熟女| 男女之事视频高清在线观看| 精品国内亚洲2022精品成人| 亚洲激情在线av| 狂野欧美激情性xxxx| 免费在线观看影片大全网站| 日韩欧美免费精品| 日本在线视频免费播放| 嫩草影视91久久| 亚洲黑人精品在线| 窝窝影院91人妻| 精品国产三级普通话版| 精品一区二区三区四区五区乱码| 成人特级黄色片久久久久久久| 夜夜爽天天搞| 成人国产综合亚洲| 国产高清有码在线观看视频| 99精品在免费线老司机午夜| 男女做爰动态图高潮gif福利片| 精品一区二区三区视频在线 | 999久久久精品免费观看国产| 啪啪无遮挡十八禁网站| 免费看光身美女| 日韩欧美免费精品| 亚洲午夜理论影院| 国产精品爽爽va在线观看网站| 99精品在免费线老司机午夜| 宅男免费午夜| 99久国产av精品| 久久久久亚洲av毛片大全| 久9热在线精品视频| 99视频精品全部免费 在线 | 老司机福利观看| 岛国在线观看网站| 岛国视频午夜一区免费看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人久久性| 最近在线观看免费完整版| 亚洲人成网站在线播放欧美日韩| 最近最新中文字幕大全免费视频| 亚洲欧美日韩东京热| 亚洲av第一区精品v没综合| 搡老熟女国产l中国老女人| 亚洲av电影在线进入| 亚洲国产欧美网| 综合色av麻豆| 一级毛片女人18水好多| www.999成人在线观看| 国产在线精品亚洲第一网站| 亚洲18禁久久av| 国产又黄又爽又无遮挡在线| 久久亚洲精品不卡| 欧美成人免费av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 最新在线观看一区二区三区| 亚洲精品在线观看二区| 亚洲人与动物交配视频| 国产99白浆流出| 天天躁日日操中文字幕| 精品福利观看| 亚洲国产欧美一区二区综合| 老熟妇乱子伦视频在线观看| 夜夜爽天天搞| 国语自产精品视频在线第100页| 岛国在线免费视频观看| 成人永久免费在线观看视频| 香蕉丝袜av| 国产精品亚洲美女久久久| 在线观看美女被高潮喷水网站 | 19禁男女啪啪无遮挡网站| 国产探花在线观看一区二区| 亚洲av成人av| 脱女人内裤的视频| 变态另类丝袜制服| 小蜜桃在线观看免费完整版高清| 成人特级av手机在线观看| 免费在线观看日本一区| 国产欧美日韩精品亚洲av| 日本一本二区三区精品| 亚洲av美国av| 欧美一级a爱片免费观看看| 亚洲av成人av| 成人一区二区视频在线观看| 丰满人妻一区二区三区视频av | 99久国产av精品| 亚洲成人中文字幕在线播放| 婷婷丁香在线五月| 国产不卡一卡二| 亚洲第一电影网av| 久久天堂一区二区三区四区| 国产 一区 欧美 日韩| 亚洲美女视频黄频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久av美女十八| 一卡2卡三卡四卡精品乱码亚洲| 男人舔女人下体高潮全视频| 欧美极品一区二区三区四区| 热99在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 男女做爰动态图高潮gif福利片| 真实男女啪啪啪动态图| 99在线人妻在线中文字幕| 真人做人爱边吃奶动态| 十八禁网站免费在线| 亚洲欧美日韩高清在线视频| 床上黄色一级片| 欧美日本视频| 国产精品女同一区二区软件 | 91九色精品人成在线观看| 亚洲性夜色夜夜综合| 特级一级黄色大片| 国产精品久久久久久人妻精品电影| 曰老女人黄片| 最近最新免费中文字幕在线| av天堂中文字幕网| 久久天躁狠狠躁夜夜2o2o| 日韩精品中文字幕看吧| 色精品久久人妻99蜜桃| 国产一区二区三区视频了| 黄色片一级片一级黄色片| 午夜福利高清视频| 成人性生交大片免费视频hd| 好男人电影高清在线观看| 成人性生交大片免费视频hd| 黄色 视频免费看| 在线观看免费午夜福利视频| 波多野结衣高清无吗| 欧美3d第一页| 成人三级黄色视频| www国产在线视频色| 国产又色又爽无遮挡免费看| 日韩欧美国产在线观看| 床上黄色一级片| 搡老熟女国产l中国老女人| 亚洲欧美日韩卡通动漫| 久久国产精品影院| 香蕉国产在线看| 波多野结衣高清无吗| 欧美黄色淫秽网站| 国产一区在线观看成人免费| 国产精品永久免费网站| 男女床上黄色一级片免费看| 99国产精品99久久久久| 国产伦人伦偷精品视频| 欧美不卡视频在线免费观看| 亚洲av第一区精品v没综合| 日本熟妇午夜| 91久久精品国产一区二区成人 | 久久热在线av| 久久精品影院6| 99久久精品一区二区三区| 嫩草影院精品99| 老司机在亚洲福利影院| 日本 欧美在线| 精品不卡国产一区二区三区| 国产淫片久久久久久久久 | xxxwww97欧美| avwww免费| 国产 一区 欧美 日韩| 国产精品一区二区三区四区久久| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女| 不卡av一区二区三区| 麻豆av在线久日| 不卡一级毛片| 午夜激情欧美在线| 久久国产乱子伦精品免费另类| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区三区四区久久| 国产成人精品久久二区二区免费| 精品国产乱码久久久久久男人| 又黄又爽又免费观看的视频| www.www免费av| 日韩欧美 国产精品| aaaaa片日本免费| 最新中文字幕久久久久 | 久久国产精品影院| a在线观看视频网站| 亚洲成人久久性| 97超级碰碰碰精品色视频在线观看| 一级a爱片免费观看的视频| www.精华液| 男女床上黄色一级片免费看| 91在线观看av| 亚洲国产精品成人综合色| 午夜两性在线视频| 狂野欧美白嫩少妇大欣赏| av国产免费在线观看| 制服人妻中文乱码| 琪琪午夜伦伦电影理论片6080| 欧美绝顶高潮抽搐喷水| 色精品久久人妻99蜜桃| 欧美又色又爽又黄视频| 亚洲人成网站高清观看| 免费人成视频x8x8入口观看| 女人高潮潮喷娇喘18禁视频| 精品久久久久久久久久免费视频| 岛国在线免费视频观看| 神马国产精品三级电影在线观看| 精品一区二区三区av网在线观看| 99久久无色码亚洲精品果冻| 亚洲真实伦在线观看| 99热精品在线国产|