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

    基于多特征融合預(yù)測(cè)蛋白質(zhì)相互作用界面

    2017-10-18 01:06:19陳心浩
    關(guān)鍵詞:殘基特征選擇分值

    陳心浩,胡 儉

    (中南民族大學(xué) 生物醫(yī)學(xué)工程學(xué)院,武漢 430074)

    基于多特征融合預(yù)測(cè)蛋白質(zhì)相互作用界面

    陳心浩,胡 儉

    (中南民族大學(xué) 生物醫(yī)學(xué)工程學(xué)院,武漢 430074)

    為高效準(zhǔn)確地預(yù)測(cè)蛋白質(zhì)相互作用界面,提取了傳統(tǒng)特征,并采用多種方法改進(jìn)進(jìn)化信息特征,利用特征選擇構(gòu)建了一個(gè)14維的預(yù)測(cè)模型.通過(guò)5折交叉驗(yàn)證和獨(dú)立測(cè)試,預(yù)測(cè)結(jié)果表明:該預(yù)測(cè)模型不僅顯著降低特征維度,而且選擇的特征組合具有較好的預(yù)測(cè)能力和較強(qiáng)的泛化能力.

    蛋白質(zhì)-蛋白質(zhì)界面;分類(lèi);進(jìn)化;特征選擇

    AbstractTo build a model of efficient and accurate classification of protein-protein interfaces, this study constructs two characteristics of traditional features and evolutionary information, a 14-dimensional feature model is constructed by feature selection.By cross-validation of the main data set and independent test set testing, results show that selects the features combination has better predictive ability and strong extension ability. Compared with the best models at the present stage, this study significantly reduce the dimensionality of the model case classification has improved.

    Keywordsprotein-protein interface; classification; evolutionary; feature selection

    區(qū)分蛋白質(zhì)晶體中的生物學(xué)相互作用界面(Biological interfaces)和無(wú)生物學(xué)意義的晶體學(xué)界面(Crystal interfaces),是結(jié)構(gòu)生物信息學(xué)中的一個(gè)重要研究方向.

    現(xiàn)有計(jì)算方法預(yù)測(cè)蛋白質(zhì)相互作用界面的特征主要分成兩大類(lèi):第一類(lèi)是以界面面積、疏水性和溫度因子等幾何特性和氨基酸理化特性為代表的傳統(tǒng)特征[1];第二類(lèi)則是以EPPIC方法為代表的進(jìn)化特征[2].為獲得良好的分類(lèi)效果,目前的主要策略是將上述特征進(jìn)行聯(lián)合.然而,這類(lèi)融合方法也存在弊端,如現(xiàn)階段分類(lèi)效果最好的Luo方法[3],該方法具有較高的特征維度(46維),且進(jìn)化信息計(jì)算復(fù)雜,不利于快速構(gòu)建本地分類(lèi)模型.因此,本文期望采用較為簡(jiǎn)便的方式計(jì)算進(jìn)化特征,融合傳統(tǒng)特征并使用特征選擇技術(shù),構(gòu)建一個(gè)低維高效的蛋白質(zhì)互作界面分類(lèi)模型.

    1 材料和方法

    1.1數(shù)據(jù)集,蛋白質(zhì)界面殘基、表面殘基的定義

    在構(gòu)建和測(cè)試模型過(guò)程中使用了三個(gè)數(shù)據(jù)集,Duarte數(shù)據(jù)集[2]作為主數(shù)據(jù)集用于構(gòu)建模型和優(yōu)化參數(shù),Bernauer[4]和Ponstingl[5]兩個(gè)經(jīng)典數(shù)據(jù)集作為獨(dú)立測(cè)試集.

    核心殘基(Core)位于互作界面中心,主要由疏水性氨基酸構(gòu)成.核心殘基周?chē)h(huán)繞著一圈殘基,此類(lèi)型殘基稱(chēng)之為環(huán)繞殘基(Rim).界面殘基、表面殘基、核心殘基與環(huán)繞殘基定義采用Proface方式定義[6].

    1.2傳統(tǒng)特征

    核心殘基(Core)與環(huán)繞殘基(Rim): 分別計(jì)算核心殘基與環(huán)繞殘基在界面殘基中的比例,構(gòu)成Core和Rim 這兩個(gè)特征.核心殘基數(shù)目(NoC) :每個(gè)蛋白質(zhì)復(fù)合體的核心殘基數(shù)構(gòu)成本特征.溫度因子(BF):將PDB中每個(gè)殘基溫度因子做Z-score歸一化,將歸一化后的界面殘基溫度因子平均值作為此蛋白質(zhì)復(fù)合體溫度因子.局部包裝密度( LD),熱點(diǎn)殘基數(shù)目(Nhs),氨基酸分布 (RP)定義方式來(lái)自Proface[6],界面疏水性 Hy)采用Jones定義[1].

    1.3進(jìn)化特征

    本文采用默認(rèn)參數(shù),使用PSI-BLAST程序?qū)δ繕?biāo)蛋白質(zhì)在NR數(shù)據(jù)庫(kù)中搜索其同源序列并構(gòu)建位置特異性矩陣.根據(jù)上述矩陣,采用Capra方法[7],對(duì)每一個(gè)殘基位置分別計(jì)算了SE(Shannon entropy of residues),SERP(Shannon entropy of residue properties),VNE(von Neumann entropy),RE(Relative Entropy)和JSD(Jensen-Shannon divergence score)5種保守性分值,并且對(duì)于計(jì)算出來(lái)的5種保守性分值采用3窗口平均,構(gòu)成另外5個(gè)保守性分值.計(jì)算公式如下:

    SEi=-∑α∈AAp(α) lg[p(α) ] ,

    (1)

    SERPi=-∑α∈Term)p(β) lg[p(β) ] ,

    (2)

    VNEi=-Tr(ρlg(ρ) ) ,

    ρ=diag[(p1,p2,…,p20)·BLUSUM62],

    (3)

    REi=-∑α∈AAp(α) lg[p(α)/q(α) ],

    (4)

    JSDi=λ∑α∈AAp(α) lg[p(α)/r(α) ]+

    (1-λ) ∑α∈AAq(α) lg[q(α)/r(α) ],

    (5)

    WindowScorei=0.5Entropyi+

    (6)

    公式(1)中,p(α)是20種常見(jiàn)氨基酸在位置i出現(xiàn)的概率,公式(2)中的p(β)則是根據(jù)Mirny研究[8]對(duì)氨基酸根據(jù)化學(xué)屬性分成6組,計(jì)算出的每一組在整體出現(xiàn)的概率,具體分組可見(jiàn)表1.VNE計(jì)算方法[9]如公式(3)所示,特點(diǎn)是將原始的概率得分使用BLUSUM62矩陣重新計(jì)算.RE的計(jì)算方式與SE接近,不同點(diǎn)是使用背景概率q(α)重新定義,其概率分布見(jiàn)表1.JSD是將RE做了背景頻率改進(jìn)[10],可以將保守性分?jǐn)?shù)歸一化0~1之間,在本文中λ=0.5.公式(6)即3窗口的算法,序列上第i個(gè)殘基與其鄰近的兩個(gè)殘基加權(quán)平均.將上述獲得的5個(gè)保守性分值和5個(gè)窗口保守性分值分別作Z-score變換,以消除不同蛋白質(zhì)復(fù)合體間差異.

    蛋白質(zhì)殘基保守性分值可以衡量殘基在進(jìn)化過(guò)程中變異程度,生物學(xué)界面殘基,特別是生物學(xué)界面上的核心殘基在進(jìn)化過(guò)程中相對(duì)保守.本文采用兩種方式計(jì)算核心殘基保守性分值[2],第一種是核心殘基-界面殘基保守性分值 (CI),計(jì)算核心殘基保守性分值平均值與界面殘基保守性分值平均值的差值,即將界面殘基保守性分值作為基準(zhǔn).第二種是核心殘基-表面殘基保守性分值( CS),計(jì)算核心殘基保守性分值平均值與表面殘基保守性分值平均值的差值.最終構(gòu)成20維進(jìn)化信息特征.

    表1 氨基酸屬性

    1.4特征選擇、分類(lèi)器與分類(lèi)評(píng)價(jià)

    增L去R選擇算法是一種改進(jìn)了的前向特征選擇方法[11].算法初始特征選擇從空集開(kāi)始,每輪先加入L維特征,然后從中除去R個(gè)特征,將每一輪AUC最高的特征組合挑選出來(lái)作為下一輪初始特征組合.

    分類(lèi)器采用R語(yǔ)言下隨機(jī)森林包,所涉及參數(shù)均采用默認(rèn)值.

    對(duì)單個(gè)特征和聯(lián)合特征測(cè)試均在Duarte數(shù)據(jù)集上完成,采用5折交叉驗(yàn)證.為排除隨機(jī)影響,5折交叉驗(yàn)證采用50次獨(dú)立分組取平均的結(jié)果,兩個(gè)獨(dú)立測(cè)試采用50次重復(fù)平均結(jié)果.分類(lèi)效果評(píng)價(jià)采用敏感度(SN)、特異度(SP)、準(zhǔn)確性、馬修相關(guān)性系數(shù)(MCC)、受試者工作曲線(ROC)及ROC曲線下面積(AUC)6個(gè)指標(biāo).MCC范圍是[-1, 1],當(dāng)MCC大于0代表正確的分類(lèi)效果,越接近1代表分類(lèi)效果越好.一般來(lái)說(shuō),當(dāng)MCC大于0.3表示有一定分類(lèi)效果,大約0.5時(shí)分類(lèi)效果較好.AUC也有類(lèi)似的評(píng)價(jià)標(biāo)準(zhǔn),當(dāng)AUC處于0.5到0.6之間表示只有微弱的分類(lèi)效果,當(dāng)AUC大于0.6表示此特征有一定的區(qū)分樣本能力,當(dāng)AUC大約0.8表示分類(lèi)效果很理想.

    2 結(jié)果與討論

    2.1特征分類(lèi)效果

    根據(jù)表2,在傳統(tǒng)特征中,Hy、Core、Rim、RP和Nhs5個(gè)特征的單獨(dú)使用分類(lèi)AUC均到達(dá)0.7以上,除Nhs每個(gè)特征的MCC都超過(guò)0.4,顯示出這些特征在生物學(xué)界面和晶體學(xué)界面上有較大的分布差異性.BF和NoC的AUC處于0.6到0.7之間,MCC大于0.3,有一定分類(lèi)效果.LD分類(lèi)效果較差,AUC不到0.6.

    表2 特征獨(dú)立使用分類(lèi)效果Tab.2 Independent feature classification results

    在進(jìn)化信息特征中,并非所有的保守性分值算法都適合本問(wèn)題,如CS-VNE的AUC小于0.5,產(chǎn)生相反的分類(lèi)效果.若以AUC為評(píng)價(jià)準(zhǔn)則,整體上來(lái)說(shuō),相同算法計(jì)算出的CS要略?xún)?yōu)于CI,這與Duarte得出的結(jié)論相同.在原始保守性分值與3維窗口計(jì)算出的保守性分?jǐn)?shù)比較中,不同種算法產(chǎn)生了不同的效果,如CI-SE-3WIN相比CI-SE分類(lèi)效果提升明顯,而CS-SE-3WIN相比CS-SE分類(lèi)效果卻變差.在20個(gè)進(jìn)化特征中,CS-SERPAUC達(dá)到0.758,MCC達(dá)到0.484,是28個(gè)特征中分類(lèi)效果最好的特征之一.

    圖1 特征選擇Fig.1 Feature selection

    2.2特征選擇

    以AUC為選擇標(biāo)準(zhǔn),本文采用增2去1選擇算法對(duì)28個(gè)特征做特征選擇.對(duì)于每一輪選擇出的特征組合,計(jì)算AUC和MCC,繪制的曲線如圖1所示,隨著特征數(shù)目的增加,AUC先快速上升,在第8輪特征選擇后達(dá)到頂點(diǎn),而后AUC緩慢下降;MCC上升速度相比于AUC較慢,而且在達(dá)到頂點(diǎn)后并沒(méi)有明顯的下降趨勢(shì).綜合AUC和MCC分值,最終選擇第14個(gè)特征組合,分別是Hy、Core、CS-SERP-3WIN、CI-SE-3WIN、RP、Nhs、CI-SE、BF、CI-RE、CI-SERP、CI-JSD-3WIN、CI-RE-3WIN、CI-SERP-3WIN、LD.選擇出的14個(gè)特征AUC為0.918,MCC為0.713,而全部28個(gè)特征AUC為0.901,MCC為0.706,可見(jiàn)本文在消減了一半特征維度情況下,AUC還是獲得了較大程度提升,說(shuō)明本文采用的特征選擇確實(shí)可以在保證預(yù)測(cè)準(zhǔn)確性條件下選擇出更有意義的特征組合.

    在特征選擇中沒(méi)有被選擇出來(lái)的特征,其中Rim是因?yàn)榕cCore成對(duì)偶關(guān)系,所包含的信息是完全一致的;NoC是因?yàn)樵诒疚闹卸鄠€(gè)特征涉及到核心殘基,信息上存在冗余因而沒(méi)有被選擇出來(lái).信息冗余同樣存在于20個(gè)進(jìn)化信息特征上,因此只有8個(gè)進(jìn)化信息特征被選擇出來(lái).雖然CS單個(gè)特征效果略好,但是在選擇出的8個(gè)進(jìn)化信息特征中只有一個(gè)CS,而獨(dú)立使用LD分類(lèi)效果較差卻可以被選擇出,說(shuō)明并非聯(lián)合較強(qiáng)特征一定會(huì)取得良好的分類(lèi)效果,還需要考慮各個(gè)特征之間的組合效應(yīng).

    2.3交叉驗(yàn)證與獨(dú)立測(cè)試效果

    表3所示的是Duarte數(shù)據(jù)集5折交叉驗(yàn)證結(jié)果和兩個(gè)獨(dú)立測(cè)試集的分類(lèi)效果,圖2所示的是相應(yīng)的ROC曲線.可以看到,本文在Duarte數(shù)據(jù)集上取得了AUC為0.918,MCC為0.713這樣良好的分類(lèi)效果.將本方法應(yīng)用于兩個(gè)獨(dú)立測(cè)試集上,Bernauer數(shù)據(jù)集AUC達(dá)到0.955,MCC達(dá)到0.745的MCC,Ponstingl數(shù)據(jù)集AUC為0.962,MCC為0.842,均獲得了良好的的分類(lèi)效果,可見(jiàn)本方法有較強(qiáng)的泛化能力.

    表3 Duarte數(shù)據(jù)集5折交叉驗(yàn)證和獨(dú)立測(cè)試集預(yù)測(cè)效果

    圖2 Duarte數(shù)據(jù)集5折交叉驗(yàn)證和獨(dú)立測(cè)試集ROC曲線Fig2 The ROC curves of 5-fold cross validation test and two independent datasets

    2.4與現(xiàn)有方法比較

    為更加全面地評(píng)價(jià)本方法,本文采用現(xiàn)階段分類(lèi)效果最好的兩個(gè)分類(lèi)器,即Luo方法和EPPIC方法對(duì)Duarte數(shù)據(jù)集做5折交叉驗(yàn)證,與本方法得到的結(jié)果進(jìn)行比較.EPPIC方法的預(yù)測(cè)效果直接取自文獻(xiàn)報(bào)道;對(duì)Luo使用的特征數(shù)據(jù),采用與本文相同的50次5折交叉驗(yàn)證進(jìn)行評(píng)價(jià).本方法與這兩種方法比較見(jiàn)圖3,從對(duì)比結(jié)果上來(lái)看,除SN本方法與現(xiàn)有方法相仿之外,SP、ACC和MCC本方法均有顯著提升,采用符號(hào)秩和檢驗(yàn)SP、ACC和MCC本方法差異達(dá)到5.24E-10、3.01E-09、5.19E-09,可以得出本方法在Duarte數(shù)據(jù)集上分類(lèi)效果優(yōu)于上述兩種方法的結(jié)論.

    圖3 本文方法與Luo方法、EPPIC比較Fig.3 Comparison of the performances of our method and Luo's Method and EPPIC

    3 結(jié)語(yǔ)

    本文提取了進(jìn)化特征和傳統(tǒng)特征,通過(guò)特征選擇構(gòu)建了一個(gè)高效的蛋白質(zhì)相互作用界面分類(lèi)模型.交叉驗(yàn)證和獨(dú)立測(cè)試的結(jié)果表明本方法可以達(dá)到較為理想的預(yù)測(cè)效果.與現(xiàn)有方法相比,本方法大幅度降低了特征維度,卻并沒(méi)有降低分類(lèi)效果.然而也有不完善的地方,如備選特征數(shù)目較少,對(duì)特征的生物學(xué)意義挖掘不深等,這些問(wèn)題將是作者下一步研究的重點(diǎn).

    [1] Jones S,JM Thornton. Analysis of protein-protein interaction sites using surface patches[J]. Journal of Molecular Biology, 1997, 272(1): 121-132.

    [2] Duarte J M, Srebniak A, Scharer, M A, et al. Protein interface classification by evolutionary analysis[J]. BMC Bioinformatics, 2012, 13(1): 334-334.

    [3] Luo J, Guo Y, Fu Y, et al. Effective discrimination between biologically relevant contacts and crystal packing contacts using new determinants[J]. Proteins, 2014, 82(11): 3090-3100.

    [4] Bernauer J, Bahadur R P, Rodier, et al. DiMoVo: a Voronoi tessellation-based method for discriminating crystallographic and biological protein-protein interactions[J]. Bioinformatics, 2008, 24(5): 652-658.

    [5] Ponstingl H, Kabir T, Thornton J M. Automatic inference of protein quaternary structure from crystals[J]. Journal of Applied Crystallography, 2003, 36(5): 1116-1122.

    [6] Saha R P, Bahadur R P, Pal A, et al. ProFace: a server for the analysis of the physicochemical features of protein-protein interfaces[J]. BMC Struct Biol, 2006, 6: 11.

    [7] Capra J A, Singh M. Predicting functionally important residues from sequence conservation[J]. Bioinformatics, 2007,23(15): 1875-1882.

    [8] Mirny L A, Shakhnovich E I. Universally conserved positions in protein folds: reading evolutionary signals about stability, folding kinetics and function[J]. Journal of Molecular Biology, 1999, 291(1): 177-196.

    [9] Caffrey D R, Somaroo S, Hughes J, et al. Are protein-protein interfaces more conserved in sequence than the rest of the protein surface[J]. Protein Science, 2004, 13(1): p. 190-202.

    [10] Lin J, Divergence measures based on the Shannon entropy[J]. IEEE Transactions on Information Theory, 1991,37(1): 145-151.

    [11] 姚 旭, 王曉丹, 張玉璽, 等. 特征選擇方法綜述[J]. 控制與決策,2012,27(2):161-166.

    StudyonProtein-ProteinInterfacialClassificationBasedonMulti-featureFusion

    ChenXinhao,HuJian

    (College of Biomedical Engineering, South-Central University for Nationalities, Wuhan 430074, China)

    Q811.4

    A

    1672-4321(2017)03-0080-04

    2017-03-30

    陳心浩(1968-),男,副教授,研究方向:醫(yī)學(xué)圖像處理與傳輸,E-mail: xinhaochen@mail.scuec.edu.cn

    國(guó)家自然科學(xué)基金資助項(xiàng)目(61002046);中央高校基本科研業(yè)務(wù)專(zhuān)項(xiàng)基金項(xiàng)目(CZP17025)

    猜你喜歡
    殘基特征選擇分值
    基于各向異性網(wǎng)絡(luò)模型研究δ阿片受體的動(dòng)力學(xué)與關(guān)鍵殘基*
    一起來(lái)看看交通違法記分分值有什么變化
    “殘基片段和排列組合法”在書(shū)寫(xiě)限制條件的同分異構(gòu)體中的應(yīng)用
    Kmeans 應(yīng)用與特征選擇
    電子制作(2017年23期)2017-02-02 07:17:06
    聯(lián)合互信息水下目標(biāo)特征選擇算法
    蛋白質(zhì)二級(jí)結(jié)構(gòu)序列與殘基種類(lèi)間關(guān)聯(lián)的分析
    基于支持向量機(jī)的蛋白質(zhì)相互作用界面熱點(diǎn)殘基預(yù)測(cè)
    基于特征選擇和RRVPMCD的滾動(dòng)軸承故障診斷方法
    宿遷城鎮(zhèn)居民醫(yī)保按病種分值結(jié)算初探
    病種分值結(jié)算模式下的醫(yī)療監(jiān)管之實(shí)踐與啟示
    亚洲情色 制服丝袜| 看黄色毛片网站| ponron亚洲| 韩国av一区二区三区四区| 精品卡一卡二卡四卡免费| 欧美最黄视频在线播放免费 | 超色免费av| 狠狠婷婷综合久久久久久88av| 性色av乱码一区二区三区2| 一级片免费观看大全| 国产精品偷伦视频观看了| 无人区码免费观看不卡| 国产男女内射视频| 亚洲国产看品久久| 99国产精品免费福利视频| 视频区欧美日本亚洲| 亚洲一码二码三码区别大吗| 国产精品秋霞免费鲁丝片| 一级片'在线观看视频| 电影成人av| 国产片内射在线| 看免费av毛片| 国产成人免费观看mmmm| 国产精品国产高清国产av | 亚洲精品在线观看二区| 两性夫妻黄色片| 麻豆成人av在线观看| av片东京热男人的天堂| 久久精品aⅴ一区二区三区四区| 9色porny在线观看| 男女下面插进去视频免费观看| 成人手机av| 婷婷丁香在线五月| 无限看片的www在线观看| 飞空精品影院首页| 国产高清激情床上av| 国产午夜精品久久久久久| 久久人妻福利社区极品人妻图片| 国产麻豆69| 久久午夜综合久久蜜桃| 国产视频一区二区在线看| tube8黄色片| 99re6热这里在线精品视频| 天天操日日干夜夜撸| 国产麻豆69| 脱女人内裤的视频| 国产99白浆流出| 国产又色又爽无遮挡免费看| 久久人妻熟女aⅴ| 别揉我奶头~嗯~啊~动态视频| 日韩欧美在线二视频 | 午夜亚洲福利在线播放| 操美女的视频在线观看| netflix在线观看网站| 国产99白浆流出| 两个人免费观看高清视频| 国内毛片毛片毛片毛片毛片| aaaaa片日本免费| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 黄片小视频在线播放| 国产精品久久久久久人妻精品电影| 18禁美女被吸乳视频| 99精品久久久久人妻精品| xxx96com| 一边摸一边做爽爽视频免费| 丰满的人妻完整版| 这个男人来自地球电影免费观看| 久久久久久久久免费视频了| 午夜精品国产一区二区电影| 人人妻人人澡人人爽人人夜夜| 一级作爱视频免费观看| 在线看a的网站| 欧美激情 高清一区二区三区| 欧美在线黄色| 黑人巨大精品欧美一区二区mp4| 美女福利国产在线| 成人影院久久| 一级片'在线观看视频| 国产亚洲欧美精品永久| 不卡av一区二区三区| 国产野战对白在线观看| 在线看a的网站| 欧美一级毛片孕妇| 黄色成人免费大全| 欧美在线黄色| 国产精品香港三级国产av潘金莲| 无限看片的www在线观看| 脱女人内裤的视频| 日日爽夜夜爽网站| 精品久久久精品久久久| 777久久人妻少妇嫩草av网站| 999精品在线视频| 亚洲色图av天堂| 国产亚洲av高清不卡| 国产精品98久久久久久宅男小说| 精品国产超薄肉色丝袜足j| 9色porny在线观看| 成年女人毛片免费观看观看9 | 国产激情欧美一区二区| 可以免费在线观看a视频的电影网站| 国产欧美亚洲国产| 丝瓜视频免费看黄片| 亚洲精品国产区一区二| 欧美日韩黄片免| 国产亚洲一区二区精品| 日韩中文字幕欧美一区二区| 亚洲精品粉嫩美女一区| 国产成人免费观看mmmm| 国产精品国产高清国产av | 9色porny在线观看| 日韩制服丝袜自拍偷拍| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成国产人片在线观看| 亚洲专区字幕在线| 日韩欧美一区二区三区在线观看 | 99国产精品99久久久久| 在线视频色国产色| 久久天躁狠狠躁夜夜2o2o| 成人18禁在线播放| 欧美 日韩 精品 国产| 国产午夜精品久久久久久| 一进一出抽搐动态| 99国产精品一区二区三区| 999久久久国产精品视频| 日韩熟女老妇一区二区性免费视频| 国产精品秋霞免费鲁丝片| 久久精品国产综合久久久| 久久久国产欧美日韩av| 亚洲精华国产精华精| av一本久久久久| 中文亚洲av片在线观看爽 | 精品久久久精品久久久| 久久久精品免费免费高清| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜添小说| 成人国产一区最新在线观看| 欧美在线一区亚洲| 国产精品一区二区免费欧美| av线在线观看网站| 夜夜爽天天搞| 精品国产乱子伦一区二区三区| 成人av一区二区三区在线看| 男女免费视频国产| 91老司机精品| 男女午夜视频在线观看| 中国美女看黄片| 男女高潮啪啪啪动态图| 精品电影一区二区在线| 亚洲av熟女| 国产精品欧美亚洲77777| 天堂俺去俺来也www色官网| 日本精品一区二区三区蜜桃| 亚洲成a人片在线一区二区| 一本大道久久a久久精品| 久久狼人影院| 欧美人与性动交α欧美精品济南到| 日日爽夜夜爽网站| 国产淫语在线视频| 制服诱惑二区| 高清欧美精品videossex| 一a级毛片在线观看| 这个男人来自地球电影免费观看| 亚洲av片天天在线观看| 国产免费av片在线观看野外av| 日韩人妻精品一区2区三区| 精品第一国产精品| 新久久久久国产一级毛片| 最近最新中文字幕大全免费视频| 亚洲精品久久午夜乱码| 国产在线观看jvid| 亚洲色图 男人天堂 中文字幕| 丁香六月欧美| 久热这里只有精品99| 又大又爽又粗| av电影中文网址| 99re在线观看精品视频| 麻豆国产av国片精品| 中文字幕高清在线视频| 一区二区日韩欧美中文字幕| 女警被强在线播放| 少妇 在线观看| 亚洲熟妇中文字幕五十中出 | 亚洲国产欧美网| 午夜福利欧美成人| 99久久国产精品久久久| 九色亚洲精品在线播放| 露出奶头的视频| 国产精品久久电影中文字幕 | 丝袜人妻中文字幕| 精品高清国产在线一区| 91麻豆精品激情在线观看国产 | 久久精品国产综合久久久| 国产高清激情床上av| 午夜精品久久久久久毛片777| 欧美乱色亚洲激情| 中文字幕人妻丝袜一区二区| 欧美另类亚洲清纯唯美| 大型黄色视频在线免费观看| 超色免费av| 搡老熟女国产l中国老女人| 国产精品久久久人人做人人爽| 别揉我奶头~嗯~啊~动态视频| 亚洲自偷自拍图片 自拍| 99热只有精品国产| 两人在一起打扑克的视频| 国产精品免费视频内射| 国产成人精品在线电影| 女警被强在线播放| 午夜福利影视在线免费观看| 成年人黄色毛片网站| 国产高清国产精品国产三级| 麻豆成人av在线观看| 男女床上黄色一级片免费看| 中国美女看黄片| 久久久久精品国产欧美久久久| 在线观看www视频免费| 久久久久久免费高清国产稀缺| 自拍欧美九色日韩亚洲蝌蚪91| 美女福利国产在线| 精品人妻1区二区| 久久中文看片网| 狠狠婷婷综合久久久久久88av| 国产欧美日韩综合在线一区二区| 一级毛片高清免费大全| av天堂在线播放| 国产有黄有色有爽视频| av网站免费在线观看视频| 大码成人一级视频| 丁香六月欧美| 在线观看免费视频日本深夜| 一区福利在线观看| 18禁裸乳无遮挡免费网站照片 | 一边摸一边抽搐一进一小说 | 美女午夜性视频免费| 久久中文字幕一级| 国产成人免费观看mmmm| 超碰成人久久| 日韩成人在线观看一区二区三区| 亚洲人成电影免费在线| 久久久精品区二区三区| 国产精品久久视频播放| 精品国产美女av久久久久小说| 99热网站在线观看| 久久久久久久国产电影| 麻豆av在线久日| 亚洲av成人一区二区三| 一二三四社区在线视频社区8| 中出人妻视频一区二区| 国产日韩欧美亚洲二区| 99久久国产精品久久久| 每晚都被弄得嗷嗷叫到高潮| av网站免费在线观看视频| 精品卡一卡二卡四卡免费| 日韩欧美三级三区| 国产精品久久电影中文字幕 | 中文亚洲av片在线观看爽 | 黄网站色视频无遮挡免费观看| 午夜免费观看网址| 久9热在线精品视频| 制服人妻中文乱码| 国产精品1区2区在线观看. | 欧美成狂野欧美在线观看| 人人妻人人添人人爽欧美一区卜| 精品国产超薄肉色丝袜足j| 国产不卡av网站在线观看| 亚洲国产欧美一区二区综合| 免费av中文字幕在线| 亚洲欧美精品综合一区二区三区| 亚洲国产看品久久| 亚洲欧美一区二区三区黑人| 大香蕉久久成人网| 人人澡人人妻人| 不卡一级毛片| 日本黄色日本黄色录像| 久久午夜亚洲精品久久| 黄色怎么调成土黄色| 午夜视频精品福利| 99re6热这里在线精品视频| 交换朋友夫妻互换小说| 久久草成人影院| 高清欧美精品videossex| 日本五十路高清| 一个人免费在线观看的高清视频| 久久午夜亚洲精品久久| 国产成人免费无遮挡视频| 亚洲美女黄片视频| 男女免费视频国产| 夫妻午夜视频| 欧美成人午夜精品| 久久久久精品人妻al黑| 夜夜夜夜夜久久久久| 久久久久国产精品人妻aⅴ院 | 一区在线观看完整版| bbb黄色大片| 9191精品国产免费久久| 欧美在线黄色| 亚洲精品美女久久久久99蜜臀| 国产麻豆69| 国产精品99久久99久久久不卡| 亚洲精品美女久久av网站| 一进一出抽搐动态| 王馨瑶露胸无遮挡在线观看| 91成人精品电影| 女性生殖器流出的白浆| x7x7x7水蜜桃| av国产精品久久久久影院| 99久久人妻综合| а√天堂www在线а√下载 | 久久国产精品男人的天堂亚洲| 国产成人精品久久二区二区91| 国产精品成人在线| 亚洲av片天天在线观看| 亚洲成人免费电影在线观看| 嫩草影视91久久| 在线观看免费视频网站a站| 久久99一区二区三区| 国产精品免费视频内射| 国产无遮挡羞羞视频在线观看| 热99久久久久精品小说推荐| 人成视频在线观看免费观看| 久久久久精品国产欧美久久久| 国产深夜福利视频在线观看| 免费人成视频x8x8入口观看| 啦啦啦视频在线资源免费观看| 国产精品免费视频内射| 亚洲熟妇熟女久久| 美女 人体艺术 gogo| 国精品久久久久久国模美| 亚洲色图av天堂| 19禁男女啪啪无遮挡网站| 少妇的丰满在线观看| 在线观看66精品国产| 村上凉子中文字幕在线| 久久久久久久精品吃奶| 91麻豆av在线| 一级毛片高清免费大全| 欧美av亚洲av综合av国产av| 少妇裸体淫交视频免费看高清 | 国产淫语在线视频| 欧美色视频一区免费| 又黄又粗又硬又大视频| 婷婷精品国产亚洲av在线 | 亚洲在线自拍视频| 下体分泌物呈黄色| 99re在线观看精品视频| 久久99一区二区三区| 91老司机精品| 久久久久久免费高清国产稀缺| 久久国产精品大桥未久av| 久久国产精品人妻蜜桃| 黄色成人免费大全| 久久精品国产亚洲av香蕉五月 | 中文亚洲av片在线观看爽 | 亚洲欧美一区二区三区久久| 精品电影一区二区在线| 亚洲aⅴ乱码一区二区在线播放 | 久久久国产成人免费| 国产99白浆流出| 国产色视频综合| 在线观看www视频免费| 男人的好看免费观看在线视频 | 国产91精品成人一区二区三区| 国产精品免费大片| 少妇被粗大的猛进出69影院| 在线播放国产精品三级| 国产一卡二卡三卡精品| 久久热在线av| 免费在线观看完整版高清| 成人18禁高潮啪啪吃奶动态图| 另类亚洲欧美激情| 黄片播放在线免费| 国产午夜精品久久久久久| 国产精品电影一区二区三区 | 99国产精品一区二区三区| 51午夜福利影视在线观看| 国产精品久久久久成人av| 久久人妻熟女aⅴ| av不卡在线播放| 国产亚洲精品久久久久5区| 精品久久久精品久久久| 最近最新中文字幕大全免费视频| 国产精品1区2区在线观看. | 水蜜桃什么品种好| 狠狠婷婷综合久久久久久88av| av不卡在线播放| 亚洲国产精品一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月天丁香| 变态另类成人亚洲欧美熟女 | 变态另类成人亚洲欧美熟女 | 热99re8久久精品国产| 国产精品永久免费网站| 久久精品亚洲熟妇少妇任你| 亚洲精品在线美女| 亚洲精品美女久久av网站| 国产熟女午夜一区二区三区| 国产在线精品亚洲第一网站| 一边摸一边做爽爽视频免费| 欧美日韩乱码在线| 国产精品偷伦视频观看了| 黄色毛片三级朝国网站| 精品福利观看| 久久香蕉国产精品| 亚洲综合色网址| 国产精品综合久久久久久久免费 | 91老司机精品| 国产高清激情床上av| 精品无人区乱码1区二区| 韩国精品一区二区三区| 热99re8久久精品国产| 中文亚洲av片在线观看爽 | 午夜精品国产一区二区电影| av片东京热男人的天堂| 久久精品熟女亚洲av麻豆精品| 国产激情久久老熟女| 他把我摸到了高潮在线观看| 淫妇啪啪啪对白视频| a级片在线免费高清观看视频| 黑人欧美特级aaaaaa片| 亚洲男人天堂网一区| 一级毛片女人18水好多| 国产精品自产拍在线观看55亚洲 | 丝袜人妻中文字幕| 每晚都被弄得嗷嗷叫到高潮| 两人在一起打扑克的视频| 亚洲第一av免费看| 亚洲精品国产色婷婷电影| 深夜精品福利| 亚洲av片天天在线观看| 亚洲av熟女| 亚洲欧美日韩另类电影网站| 香蕉国产在线看| 亚洲 国产 在线| 亚洲一区高清亚洲精品| 久久精品aⅴ一区二区三区四区| 国产精品 欧美亚洲| 欧美激情 高清一区二区三区| 精品少妇久久久久久888优播| videosex国产| 亚洲国产精品合色在线| 欧美日韩成人在线一区二区| 天天躁夜夜躁狠狠躁躁| 欧美一级毛片孕妇| 久久中文看片网| 亚洲欧美日韩高清在线视频| 免费在线观看完整版高清| 国产精品一区二区在线不卡| 亚洲国产看品久久| 色婷婷av一区二区三区视频| 激情视频va一区二区三区| 国产亚洲精品久久久久5区| www.精华液| 欧美日韩黄片免| 精品久久久久久久毛片微露脸| 最近最新免费中文字幕在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品在线观看二区| 久久中文字幕人妻熟女| 国产精品亚洲av一区麻豆| 亚洲五月色婷婷综合| 精品国产一区二区三区久久久樱花| 中文欧美无线码| 丰满迷人的少妇在线观看| 18禁美女被吸乳视频| 中亚洲国语对白在线视频| 欧美人与性动交α欧美软件| 午夜福利免费观看在线| 国产成人av激情在线播放| 精品一品国产午夜福利视频| 欧美成人免费av一区二区三区 | 亚洲欧美色中文字幕在线| 欧美成人免费av一区二区三区 | 亚洲五月婷婷丁香| 国产精品亚洲一级av第二区| 人成视频在线观看免费观看| 亚洲色图av天堂| 男女下面插进去视频免费观看| 女性被躁到高潮视频| а√天堂www在线а√下载 | 国产伦人伦偷精品视频| 脱女人内裤的视频| 99国产精品免费福利视频| 男女高潮啪啪啪动态图| 人妻 亚洲 视频| av片东京热男人的天堂| 韩国精品一区二区三区| 人人妻人人澡人人看| 精品电影一区二区在线| 国产欧美日韩综合在线一区二区| 国产无遮挡羞羞视频在线观看| 久9热在线精品视频| 欧美 日韩 精品 国产| 中文字幕人妻熟女乱码| 12—13女人毛片做爰片一| 欧美乱妇无乱码| 欧美精品啪啪一区二区三区| 国产有黄有色有爽视频| 精品国产国语对白av| 日韩熟女老妇一区二区性免费视频| 看片在线看免费视频| 亚洲av日韩精品久久久久久密| 一级毛片高清免费大全| 99精品欧美一区二区三区四区| 首页视频小说图片口味搜索| 婷婷精品国产亚洲av在线 | 久久人人97超碰香蕉20202| 少妇裸体淫交视频免费看高清 | 黄色 视频免费看| 男女免费视频国产| 日日摸夜夜添夜夜添小说| 视频区欧美日本亚洲| 一本大道久久a久久精品| 啪啪无遮挡十八禁网站| 国产精品国产高清国产av | 757午夜福利合集在线观看| 午夜亚洲福利在线播放| 精品福利观看| 69av精品久久久久久| 水蜜桃什么品种好| 亚洲国产欧美一区二区综合| 欧美日韩精品网址| 久久久久久久久久久久大奶| 国产精品98久久久久久宅男小说| 一区二区三区精品91| 亚洲精品国产区一区二| 欧美精品一区二区免费开放| 久久久久久久国产电影| 国产精品亚洲av一区麻豆| 欧美在线黄色| 国产精品免费大片| 亚洲精品一二三| 怎么达到女性高潮| av有码第一页| 老司机深夜福利视频在线观看| 国内毛片毛片毛片毛片毛片| 一边摸一边抽搐一进一出视频| 视频在线观看一区二区三区| 精品国产超薄肉色丝袜足j| 这个男人来自地球电影免费观看| 麻豆乱淫一区二区| 操出白浆在线播放| 精品亚洲成a人片在线观看| 欧美丝袜亚洲另类 | 国产精品久久久人人做人人爽| 最近最新中文字幕大全电影3 | 亚洲国产欧美一区二区综合| 国产亚洲一区二区精品| 一进一出抽搐动态| 婷婷精品国产亚洲av在线 | 国产日韩欧美亚洲二区| 国产男女内射视频| 国产精品免费一区二区三区在线 | 成人av一区二区三区在线看| 日本撒尿小便嘘嘘汇集6| 久久国产精品男人的天堂亚洲| 欧美大码av| 亚洲中文av在线| 99国产精品一区二区三区| 日韩人妻精品一区2区三区| 美女视频免费永久观看网站| 很黄的视频免费| 国产一区在线观看成人免费| 正在播放国产对白刺激| 精品人妻熟女毛片av久久网站| 免费观看精品视频网站| 丝袜美腿诱惑在线| 久久精品亚洲熟妇少妇任你| 亚洲av欧美aⅴ国产| 国产精品一区二区精品视频观看| 欧美色视频一区免费| 悠悠久久av| 成人国语在线视频| 波多野结衣av一区二区av| 成年女人毛片免费观看观看9 | 午夜精品在线福利| 国产成人欧美| 无人区码免费观看不卡| 亚洲国产精品sss在线观看 | 国产av又大| ponron亚洲| 国产精品秋霞免费鲁丝片| av天堂在线播放| 亚洲综合色网址| 在线观看一区二区三区激情| 女人精品久久久久毛片| tube8黄色片| 国产成人欧美在线观看 | 如日韩欧美国产精品一区二区三区| 一边摸一边做爽爽视频免费| 最近最新免费中文字幕在线| 国产乱人伦免费视频| 久久天躁狠狠躁夜夜2o2o| 又紧又爽又黄一区二区| 香蕉国产在线看| 成人特级黄色片久久久久久久| tocl精华| 久久久久久免费高清国产稀缺| 欧美日韩视频精品一区| 国产欧美日韩一区二区三区在线| 欧美在线一区亚洲| 免费在线观看影片大全网站| 国产欧美日韩一区二区精品| 无遮挡黄片免费观看| 搡老乐熟女国产| 午夜福利免费观看在线| 国产乱人伦免费视频| 别揉我奶头~嗯~啊~动态视频| 亚洲国产毛片av蜜桃av| 伦理电影免费视频| 久久精品人人爽人人爽视色| 色老头精品视频在线观看| 如日韩欧美国产精品一区二区三区| 国产日韩欧美亚洲二区| 欧美精品av麻豆av|