• <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í)踐與啟示
    亚洲最大成人中文| 欧美日韩乱码在线| 在线免费十八禁| 国产精品永久免费网站| 亚洲色图av天堂| 久久久午夜欧美精品| 天堂√8在线中文| 亚洲五月天丁香| 久久草成人影院| 日日摸夜夜添夜夜添小说| 亚洲专区国产一区二区| 国产高清激情床上av| 搡老妇女老女人老熟妇| 亚洲av不卡在线观看| 午夜视频国产福利| 99国产极品粉嫩在线观看| 99在线人妻在线中文字幕| a在线观看视频网站| xxxwww97欧美| 欧美日韩国产亚洲二区| av在线观看视频网站免费| 国产精品日韩av在线免费观看| 精品欧美国产一区二区三| 日韩欧美国产在线观看| 午夜福利在线观看免费完整高清在 | 成熟少妇高潮喷水视频| 国产伦在线观看视频一区| 欧美+亚洲+日韩+国产| 亚洲内射少妇av| 亚洲av免费高清在线观看| 日韩欧美国产在线观看| 国内精品美女久久久久久| 我的女老师完整版在线观看| 亚洲四区av| 最近最新中文字幕大全电影3| 少妇的逼水好多| 免费在线观看日本一区| 国产真实伦视频高清在线观看 | 日本精品一区二区三区蜜桃| a级毛片a级免费在线| 亚洲五月天丁香| 欧美色欧美亚洲另类二区| 久久国产精品人妻蜜桃| 亚洲国产日韩欧美精品在线观看| 人妻制服诱惑在线中文字幕| 色哟哟哟哟哟哟| 欧美日韩精品成人综合77777| 哪里可以看免费的av片| 男女视频在线观看网站免费| 成人午夜高清在线视频| 日本撒尿小便嘘嘘汇集6| 亚洲人成网站高清观看| 特级一级黄色大片| 嫩草影视91久久| 亚洲精品成人久久久久久| 美女高潮喷水抽搐中文字幕| 嫩草影院精品99| 一区二区三区免费毛片| 国产精品美女特级片免费视频播放器| 美女高潮的动态| 中文字幕人妻熟人妻熟丝袜美| 天堂影院成人在线观看| 午夜久久久久精精品| a级一级毛片免费在线观看| 国产亚洲精品综合一区在线观看| 国产成人一区二区在线| 亚洲精品一卡2卡三卡4卡5卡| 欧美三级亚洲精品| 大型黄色视频在线免费观看| 自拍偷自拍亚洲精品老妇| 日韩精品有码人妻一区| 国产精品免费一区二区三区在线| 91在线精品国自产拍蜜月| 真人一进一出gif抽搐免费| 日韩欧美在线二视频| 久久久久免费精品人妻一区二区| 亚洲色图av天堂| 18禁黄网站禁片免费观看直播| 亚洲内射少妇av| 91久久精品国产一区二区成人| 久久精品国产亚洲网站| 欧美极品一区二区三区四区| 女的被弄到高潮叫床怎么办 | 国产国拍精品亚洲av在线观看| 麻豆成人午夜福利视频| 欧美日本视频| 国产主播在线观看一区二区| 久久久色成人| 成人精品一区二区免费| 国产真实伦视频高清在线观看 | 亚洲电影在线观看av| 欧美黑人欧美精品刺激| 久久精品国产亚洲网站| 亚洲av二区三区四区| 亚洲综合色惰| 亚洲欧美清纯卡通| 亚洲中文日韩欧美视频| 午夜影院日韩av| 搡老妇女老女人老熟妇| 欧美高清性xxxxhd video| av天堂中文字幕网| 成年女人永久免费观看视频| 欧美色欧美亚洲另类二区| 欧美人与善性xxx| 午夜福利欧美成人| 亚洲四区av| 国产伦精品一区二区三区四那| 男人舔奶头视频| 欧美性猛交╳xxx乱大交人| 免费看av在线观看网站| 欧美日韩综合久久久久久 | 亚洲av电影不卡..在线观看| 欧美区成人在线视频| 欧美黑人欧美精品刺激| 日本a在线网址| 国产精品久久视频播放| 久久久久久久午夜电影| 黄色女人牲交| 国产黄色小视频在线观看| 韩国av在线不卡| 国产精品野战在线观看| 亚洲av五月六月丁香网| 中文亚洲av片在线观看爽| 精品人妻偷拍中文字幕| 久久精品国产亚洲av香蕉五月| 日本爱情动作片www.在线观看 | 国产高清视频在线观看网站| 简卡轻食公司| 黄色一级大片看看| 国产一区二区在线av高清观看| 午夜激情欧美在线| 成年版毛片免费区| 亚洲av免费在线观看| 99精品久久久久人妻精品| 韩国av一区二区三区四区| 日韩中字成人| www.www免费av| 久久久精品大字幕| 亚洲精品在线观看二区| 久久久久久久精品吃奶| 日日摸夜夜添夜夜添小说| 村上凉子中文字幕在线| 午夜免费激情av| 国产欧美日韩精品一区二区| 一a级毛片在线观看| 亚洲精品久久国产高清桃花| netflix在线观看网站| 欧美日韩精品成人综合77777| 欧美成人一区二区免费高清观看| 日韩欧美在线乱码| 免费av毛片视频| av女优亚洲男人天堂| 国产精品不卡视频一区二区| 美女黄网站色视频| 亚洲精品456在线播放app | 人妻少妇偷人精品九色| 日韩,欧美,国产一区二区三区 | 窝窝影院91人妻| 日本与韩国留学比较| 欧美zozozo另类| 国产高清三级在线| 精华霜和精华液先用哪个| 久久久久国内视频| 亚洲国产日韩欧美精品在线观看| avwww免费| 日本一二三区视频观看| 伊人久久精品亚洲午夜| 精品久久久久久久久av| 日本三级黄在线观看| 国产黄片美女视频| 国产高清激情床上av| 日韩 亚洲 欧美在线| 中文在线观看免费www的网站| 亚洲av中文字字幕乱码综合| 在线a可以看的网站| 男女下面进入的视频免费午夜| 男女视频在线观看网站免费| 国内久久婷婷六月综合欲色啪| av在线老鸭窝| 色综合婷婷激情| 午夜福利高清视频| 成人三级黄色视频| 色吧在线观看| 听说在线观看完整版免费高清| 日日撸夜夜添| 亚洲自偷自拍三级| 国产精品,欧美在线| 日韩强制内射视频| 免费观看人在逋| 很黄的视频免费| 亚洲在线观看片| 少妇的逼水好多| 久久精品国产亚洲av涩爱 | 制服丝袜大香蕉在线| 麻豆国产av国片精品| 色av中文字幕| 伊人久久精品亚洲午夜| 神马国产精品三级电影在线观看| 亚洲av成人av| 一区二区三区四区激情视频 | 少妇被粗大猛烈的视频| 国产一区二区三区在线臀色熟女| 99在线视频只有这里精品首页| 免费观看精品视频网站| 午夜福利在线观看免费完整高清在 | av在线老鸭窝| 国产真实伦视频高清在线观看 | 亚洲精品在线观看二区| 久久久色成人| av福利片在线观看| 日本免费a在线| 男女那种视频在线观看| 午夜福利欧美成人| 人妻少妇偷人精品九色| 亚洲中文字幕日韩| 丝袜美腿在线中文| av在线亚洲专区| 亚洲av免费在线观看| 成人二区视频| 日本免费一区二区三区高清不卡| 欧美成人一区二区免费高清观看| 国产精品乱码一区二三区的特点| or卡值多少钱| 免费在线观看日本一区| 又粗又爽又猛毛片免费看| 噜噜噜噜噜久久久久久91| 日本五十路高清| 久久久久久九九精品二区国产| 国产精品久久久久久精品电影| 亚洲av熟女| 很黄的视频免费| 亚洲av中文字字幕乱码综合| 欧美绝顶高潮抽搐喷水| 神马国产精品三级电影在线观看| 亚洲av美国av| 亚洲最大成人手机在线| 亚洲人成网站在线播| 午夜视频国产福利| av中文乱码字幕在线| 国产欧美日韩精品亚洲av| 给我免费播放毛片高清在线观看| 久久精品国产亚洲av天美| 日日摸夜夜添夜夜添av毛片 | 一个人观看的视频www高清免费观看| 中文字幕免费在线视频6| 精品福利观看| 日本精品一区二区三区蜜桃| 国产av一区在线观看免费| 国产av在哪里看| 毛片女人毛片| 波野结衣二区三区在线| 成人午夜高清在线视频| 久久久久久久久久成人| 国产精品爽爽va在线观看网站| 观看免费一级毛片| www.色视频.com| 久久99热6这里只有精品| 亚洲电影在线观看av| 亚洲欧美日韩卡通动漫| 亚洲性久久影院| 国产欧美日韩精品亚洲av| 九九爱精品视频在线观看| 人人妻,人人澡人人爽秒播| 国产亚洲av嫩草精品影院| 日韩国内少妇激情av| 久久久久久久午夜电影| 99精品在免费线老司机午夜| 国产成年人精品一区二区| 中文字幕免费在线视频6| 亚洲精品在线观看二区| 久9热在线精品视频| 波野结衣二区三区在线| 日韩av在线大香蕉| 国产探花极品一区二区| 欧美一区二区亚洲| 中文字幕高清在线视频| АⅤ资源中文在线天堂| 99热这里只有是精品50| 我的女老师完整版在线观看| 久久久久久久久久久丰满 | 亚洲国产色片| 日本-黄色视频高清免费观看| 欧美激情在线99| 亚洲午夜理论影院| 欧美一区二区国产精品久久精品| 成人一区二区视频在线观看| 国产男人的电影天堂91| 亚洲av成人精品一区久久| 午夜a级毛片| xxxwww97欧美| av中文乱码字幕在线| 精品一区二区三区视频在线| 午夜福利在线观看免费完整高清在 | 国产三级在线视频| 午夜激情福利司机影院| 亚洲 国产 在线| 久久国内精品自在自线图片| 欧美日韩综合久久久久久 | 亚洲国产日韩欧美精品在线观看| 床上黄色一级片| 国产综合懂色| 精品国产三级普通话版| 午夜福利成人在线免费观看| 精品久久久久久久久亚洲 | 99久久精品热视频| 别揉我奶头 嗯啊视频| 亚洲精品一区av在线观看| 国产精品一区www在线观看 | 欧美3d第一页| 村上凉子中文字幕在线| 麻豆av噜噜一区二区三区| 国产精品综合久久久久久久免费| 麻豆国产av国片精品| 熟女电影av网| 男插女下体视频免费在线播放| 国产乱人伦免费视频| 亚洲av第一区精品v没综合| 91在线观看av| 综合色av麻豆| 18+在线观看网站| 88av欧美| 精品久久久噜噜| 色吧在线观看| 又黄又爽又刺激的免费视频.| avwww免费| 极品教师在线视频| 黄片wwwwww| 女人被狂操c到高潮| 成人国产一区最新在线观看| 热99re8久久精品国产| 男插女下体视频免费在线播放| 国产成人一区二区在线| 观看美女的网站| 日日摸夜夜添夜夜添av毛片 | 成年版毛片免费区| 天天躁日日操中文字幕| 国产黄片美女视频| 久久午夜福利片| 白带黄色成豆腐渣| а√天堂www在线а√下载| 在线播放无遮挡| a级毛片免费高清观看在线播放| 联通29元200g的流量卡| 亚洲国产欧美人成| 波多野结衣巨乳人妻| 麻豆国产97在线/欧美| 亚洲欧美日韩高清在线视频| 免费高清视频大片| 国产精品野战在线观看| 九色成人免费人妻av| 黄色欧美视频在线观看| 久久久久久伊人网av| 内地一区二区视频在线| 国产不卡一卡二| 免费看日本二区| 日韩欧美精品v在线| 国产视频内射| 午夜福利18| 最近在线观看免费完整版| 少妇高潮的动态图| 乱人视频在线观看| 欧美潮喷喷水| 99久久中文字幕三级久久日本| 国产精品伦人一区二区| 国产熟女欧美一区二区| 国产精品日韩av在线免费观看| 亚洲第一电影网av| 99久久久亚洲精品蜜臀av| 久久久久久久精品吃奶| 午夜视频国产福利| 免费电影在线观看免费观看| 国产成人av教育| 国产黄片美女视频| www.www免费av| 国产老妇女一区| 最后的刺客免费高清国语| 欧美日韩综合久久久久久 | 两性午夜刺激爽爽歪歪视频在线观看| 成人午夜高清在线视频| 欧美最新免费一区二区三区| 春色校园在线视频观看| 麻豆久久精品国产亚洲av| 精品一区二区三区视频在线| 亚洲精品亚洲一区二区| 夜夜夜夜夜久久久久| 精品人妻1区二区| 在线看三级毛片| 午夜福利18| 国产精品久久久久久久久免| 成人国产麻豆网| 日本在线视频免费播放| 美女高潮喷水抽搐中文字幕| 窝窝影院91人妻| 欧美性猛交黑人性爽| 国产精品不卡视频一区二区| 最近中文字幕高清免费大全6 | 乱人视频在线观看| 日韩国内少妇激情av| 岛国在线免费视频观看| av福利片在线观看| 一本一本综合久久| 亚洲精品亚洲一区二区| 国产单亲对白刺激| 国产伦精品一区二区三区视频9| 在线播放国产精品三级| 免费不卡的大黄色大毛片视频在线观看 | 99热这里只有是精品50| 精品久久国产蜜桃| 色5月婷婷丁香| 亚洲av第一区精品v没综合| 99久久九九国产精品国产免费| 亚洲精华国产精华精| 中文亚洲av片在线观看爽| 久久午夜亚洲精品久久| 男人舔奶头视频| 干丝袜人妻中文字幕| 精品久久久久久久久亚洲 | 精品久久国产蜜桃| 日日摸夜夜添夜夜添av毛片 | 免费观看在线日韩| bbb黄色大片| 国产精品福利在线免费观看| 亚洲男人的天堂狠狠| 久久久成人免费电影| 亚洲天堂国产精品一区在线| 国产黄a三级三级三级人| 少妇的逼水好多| 日韩欧美免费精品| 国产黄片美女视频| 一本一本综合久久| 欧美又色又爽又黄视频| 欧美成人免费av一区二区三区| 最新中文字幕久久久久| 国产在线精品亚洲第一网站| 欧美性猛交黑人性爽| 欧美+日韩+精品| 免费搜索国产男女视频| 国产一区二区三区视频了| 亚洲人与动物交配视频| 成人三级黄色视频| 国产在视频线在精品| 日本精品一区二区三区蜜桃| 亚洲午夜理论影院| 免费人成视频x8x8入口观看| 全区人妻精品视频| 国产成人av教育| 国产一区二区三区av在线 | 美女高潮的动态| 欧美bdsm另类| 亚洲av成人精品一区久久| 欧美一区二区国产精品久久精品| 伊人久久精品亚洲午夜| 久久久久九九精品影院| 亚洲精品亚洲一区二区| 一卡2卡三卡四卡精品乱码亚洲| 国产伦精品一区二区三区视频9| 国产色爽女视频免费观看| 国产一区二区亚洲精品在线观看| 久久精品国产亚洲av涩爱 | 亚洲七黄色美女视频| 小说图片视频综合网站| 一进一出抽搐动态| 狠狠狠狠99中文字幕| 国产亚洲精品综合一区在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产老妇女一区| 国产人妻一区二区三区在| 最好的美女福利视频网| 色尼玛亚洲综合影院| 免费高清视频大片| 熟妇人妻久久中文字幕3abv| 色精品久久人妻99蜜桃| 午夜福利高清视频| 老司机福利观看| 精品乱码久久久久久99久播| 日韩在线高清观看一区二区三区 | 在线播放无遮挡| 免费人成视频x8x8入口观看| 悠悠久久av| 国产视频内射| 亚洲综合色惰| 赤兔流量卡办理| 91在线精品国自产拍蜜月| 欧美激情久久久久久爽电影| 亚洲av熟女| 又黄又爽又免费观看的视频| 毛片女人毛片| 久久精品人妻少妇| 美女 人体艺术 gogo| ponron亚洲| 如何舔出高潮| 国产精品一区二区三区四区久久| 中文资源天堂在线| а√天堂www在线а√下载| 全区人妻精品视频| 99热这里只有是精品在线观看| 搡女人真爽免费视频火全软件 | 亚洲成av人片在线播放无| 久久6这里有精品| 999久久久精品免费观看国产| 非洲黑人性xxxx精品又粗又长| av国产免费在线观看| 亚洲在线自拍视频| 观看美女的网站| 全区人妻精品视频| 日本黄色片子视频| 美女黄网站色视频| 国产午夜精品论理片| 在线观看午夜福利视频| 波多野结衣高清无吗| 久久精品国产亚洲av天美| 国产一区二区在线av高清观看| 欧美精品啪啪一区二区三区| 美女cb高潮喷水在线观看| 人人妻,人人澡人人爽秒播| 国产黄片美女视频| 熟女人妻精品中文字幕| 3wmmmm亚洲av在线观看| 精品国内亚洲2022精品成人| 成人午夜高清在线视频| 日韩欧美 国产精品| 一级av片app| 亚洲自偷自拍三级| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲四区av| 精品人妻视频免费看| 日本一本二区三区精品| 国产伦人伦偷精品视频| 国产一区二区三区视频了| 免费看av在线观看网站| 亚洲18禁久久av| 搞女人的毛片| 欧美色视频一区免费| 少妇人妻一区二区三区视频| 国内精品久久久久精免费| 最近中文字幕高清免费大全6 | 国模一区二区三区四区视频| 日日啪夜夜撸| 国产又黄又爽又无遮挡在线| 偷拍熟女少妇极品色| 亚洲国产高清在线一区二区三| 色视频www国产| 久久中文看片网| 亚洲无线观看免费| 亚洲va在线va天堂va国产| 国产精品福利在线免费观看| 久久久久久久亚洲中文字幕| 午夜亚洲福利在线播放| 国内精品美女久久久久久| 国产v大片淫在线免费观看| 久久人人爽人人爽人人片va| 国产精品一区二区三区四区免费观看 | 国产日本99.免费观看| 91麻豆精品激情在线观看国产| 中国美白少妇内射xxxbb| 日本一二三区视频观看| 91av网一区二区| 成人一区二区视频在线观看| 毛片一级片免费看久久久久 | 一本一本综合久久| 黄片wwwwww| 国产激情偷乱视频一区二区| 午夜福利在线观看吧| av在线观看视频网站免费| 国产色爽女视频免费观看| 校园人妻丝袜中文字幕| 人妻夜夜爽99麻豆av| 99热这里只有精品一区| 国产精品爽爽va在线观看网站| 99视频精品全部免费 在线| 欧美3d第一页| 国产av不卡久久| 亚洲经典国产精华液单| 麻豆精品久久久久久蜜桃| 国产精品免费一区二区三区在线| 熟妇人妻久久中文字幕3abv| 别揉我奶头~嗯~啊~动态视频| 草草在线视频免费看| 午夜日韩欧美国产| 俺也久久电影网| 日本色播在线视频| 日日夜夜操网爽| 久久天躁狠狠躁夜夜2o2o| 国产男靠女视频免费网站| 久久久精品大字幕| 少妇猛男粗大的猛烈进出视频 | 永久网站在线| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久午夜电影| 久久香蕉精品热| 如何舔出高潮| 丰满乱子伦码专区| 亚洲国产高清在线一区二区三| 婷婷丁香在线五月| 久久久久久国产a免费观看| 欧美区成人在线视频| 91精品国产九色| 欧美一区二区国产精品久久精品| 亚洲人成网站在线播放欧美日韩| 免费看av在线观看网站| 国产精品三级大全| 小说图片视频综合网站| 蜜桃亚洲精品一区二区三区| 国产一区二区激情短视频| 九九在线视频观看精品| 精品午夜福利视频在线观看一区| 九九热线精品视视频播放| av.在线天堂| 精品久久久久久成人av| 亚洲性夜色夜夜综合| 国产免费一级a男人的天堂| 在线播放国产精品三级| 一夜夜www| 国产高清视频在线播放一区| 日韩欧美国产在线观看|