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

    有機磷類化合物大鼠急性毒性QSAR模型構(gòu)建與毒性機制研究

    2022-06-06 08:07:56鄭子廷閆賽紅查金苗
    生態(tài)毒理學報 2022年1期
    關鍵詞:實驗模型研究

    鄭子廷,閆賽紅,#,查金苗,*

    1. 中國科學院生態(tài)環(huán)境研究中心,環(huán)境水質(zhì)學國家重點實驗室,北京 100085 2. 中國科學院大學,北京 100049

    有機磷類化合物(organophosphorus compounds, OPs)包括有機磷阻燃劑(organophosphorus flame retardants, OPFRs)、有機磷殺蟲劑(organophosphorus pesticides, OPPs)和其他含磷的有機化合物。工業(yè)生產(chǎn)上,OFPRs作為溴化阻燃劑的低毒替代品被應用于家電、建筑和電子領域;農(nóng)業(yè)上,OPPs則被廣泛用于殺蟲和害蟲控制[1-2]。此外在醫(yī)學領域,一些抗病毒的藥物中也含有OPs[3]。因此,OPs容易通過磨損、揮發(fā)和溶解等方式進入各種環(huán)境介質(zhì)中。目前科研人員已在諸如河流、地下水、海水、空氣和室內(nèi)等環(huán)境介質(zhì)中檢測到了OPs[4-6]。甚至在東亞地區(qū)的人類母乳樣本中檢測到10余種OPs,最高濃度為70 ng·g-1[7]。越來越廣泛存在的OPs可能造成的危害和健康風險引起了研究者們的關注。

    OPs具有神經(jīng)毒性[8-9]、內(nèi)分泌毒性[10]和致癌性[11]。研究表明,接觸OPPs(如毒死蜱)會引起人類的神經(jīng)系統(tǒng)發(fā)育缺陷[12]。OPFRs如磷酸三甲苯(TMPP)、磷酸三苯酯(TPHP)和磷酸三(2,3-二溴丙基)酯(TDBPP)則顯示出對核受體(nuclear receptors, NRs)的親和力和內(nèi)分泌的干擾作用[13-14]。上述研究及更多OPs的毒性毒理研究往往都是通過傳統(tǒng)的體內(nèi)(invivo)或者體外(invitro)實驗進行的,這些方法需要耗費大量的時間與金錢。此外,實驗手段無法完整覆蓋一整類化合物的毒性數(shù)據(jù),并且測試在實驗室間存在誤差[15]。因此,迫切需要一種有效的替代方法,對化合物毒性數(shù)據(jù)進行評估。

    隨著計算科學的發(fā)展,計算毒理學已經(jīng)成為一種具有成本效益且快速的實驗替代方法[16]。計算毒理學整合各種來源的數(shù)據(jù),開發(fā)基于數(shù)學、統(tǒng)計科學和機器學習的模型,對化學品進行高效和高通量的定性或定量預測[17]。定量結(jié)構(gòu)活性關系(quantitative structure-activity relationship, QSAR)模型便是計算毒理學領域被廣泛利用的模型之一,其原理是建立化合物的化學結(jié)構(gòu)與各種毒理學終點效應和生物活性之間的函數(shù)關系[18]。QSAR模型被研究者們廣泛應用于藥物發(fā)現(xiàn)、毒理學和生物學等領域[19]。除了建立化學結(jié)構(gòu)與毒性數(shù)據(jù)之間的關系模型外,QSAR模型還能通過分子描述符的解釋從一定程度進行毒理機制的闡述。Alves等[20]使用QSAR方法為COVID-19病毒建立了抑制性藥物預測模型,并應用該模型從小分子數(shù)據(jù)庫中高通量地篩選了11個候選藥物,其中3個被美國轉(zhuǎn)化科學推動中心(National Center for Advancing Translational Sciences, NCATS)實驗證實為有效的病毒抑制劑。Jeon等[21]則收集大量實驗數(shù)據(jù),為魚胚胎急性毒性建立了一個通用的基線毒性QSAR模型,并使用該模型篩選毒性化合物。大量的研究表明QSAR模型用于化合物的毒性預測行之有效。

    分子對接是研究化合物與目標蛋白結(jié)合部位行為的工具,能夠通過計算手段探究化合物與特定靶點蛋白的結(jié)合特異性,從而從一定程度上反映化合物的毒性。Zou等[22]使用分子對接得到的配體-受體相互作用能量(binding energy)來反映化學品的毒性作用。Li等[23]通過對接羥基多溴聯(lián)苯醚(HO-PBDEs)與人類甲狀腺激素受體,從一定程度上揭示了HO-PBDEs的毒性機制。

    現(xiàn)有的一些針對大鼠急性口服毒性的QSAR模型存在沒有使用測試集[24]、選用的分子描述符數(shù)量不規(guī)范[25]等問題。在本研究中,完善地建立了基于逐步算法的線性回歸(SW-MLR)和基于遺傳算法的線性回歸(GA-MLR)QSAR模型用于預測大鼠的急性口服毒性。研究收集的數(shù)據(jù)集中含有9種沒有毒性實驗數(shù)據(jù)的OPs,最終開發(fā)完成的QSAR模型將用于預測這9種OPs的急性毒性。分子對接也將作為一個輔助工具來評估這9種OPs潛在的神經(jīng)毒性。

    1 材料與方法(Materials and methods)

    1.1 數(shù)據(jù)集的收集和劃分

    本研究首先收集了OPs對大鼠急性口服毒性的數(shù)據(jù)。檢索并從文獻中提取了62種環(huán)境中常見的OPs清單[26-27]。通過查詢EPA CompTox數(shù)據(jù)庫(https://comptox.epa.gov/dashboard),收集了同一實驗條件下獲得的53種OPs急性毒性數(shù)據(jù)LD50(導致一組實驗動物中剛好死亡一半的化合物濃度)。獲取的OPs的急性毒性數(shù)值使用log函數(shù)進行規(guī)范化處理。整個數(shù)據(jù)集按照4∶1的比例隨機劃分為訓練集(n=43)和測試集(n=10),沒有毒性實驗數(shù)據(jù)的9種OPs作為最后的預測集(n=9)。在之后的計算實驗中,訓練集僅用于建立模型,測試集則用于評估模型的預測能力,這種劃分方式有利于驗證QSAR模型的擬合和泛化能力。

    1.2 分子描述符的計算

    數(shù)據(jù)集中所有OPs的結(jié)構(gòu)文件(.sdf)來自PubChem(https://pubchem.ncbi.nlm.nih.gov/)數(shù)據(jù)庫。并使用PaDEL開源軟件[28]計算所有化合物的分子結(jié)構(gòu)描述符,共得到1 444個分子結(jié)構(gòu)描述符,包括(1)一維結(jié)構(gòu)描述符:如原子和化學鍵的數(shù)量、分子量等;(2)二維拓撲描述符:如拓撲指數(shù)、拓撲電荷和自相關系數(shù)等;(3)電荷描述符、原子電狀態(tài)描述符等。本研究中使用了PaDEL計算的一維和二維描述符,因為這些類型的描述符計算速度快,可重復性好,并且不需要分子結(jié)構(gòu)能量最小化過程,該過程中不同實驗室間可能會在計算上產(chǎn)生誤差。計算結(jié)束后,我們對得到的分子描述符數(shù)據(jù)進行了預處理,剔除了方差為0的單列分子描述符數(shù)據(jù),相關系數(shù)>0.99的多組描述符最終僅保留一組。經(jīng)過預處理后,最終保留了833組分子描述符用于后續(xù)分析。

    1.3 模型的構(gòu)建

    分別使用逐步算法(SW)和遺傳算法(GA)進行分子描述符的選擇。2種算法選擇的分子描述符變量通過基于普通最小二乘法(OLS)的多元線性回歸(MLR)方法關聯(lián)毒性終點并建立QSAR模型,之后進行模型的內(nèi)部和外部評估。研究中使用SPSS 21.0軟件建立了SW-MLR模型,使用DTC-QSAR開源軟件[29]建立了GA-MLR模型。所有參數(shù)的統(tǒng)計分析過程在EXCEL軟件中完成。本研究詳細的QSAR模型構(gòu)建及OPs的LD50預測工作流程如圖1所示。

    圖1 定量結(jié)構(gòu)活性關系(QSAR)模型構(gòu)建及有機磷類化合物(OPs)的急性毒性半致死量(LD50)預測的工作流程圖注:GA-MLR表示基于遺傳算法的多元線性回歸模型;SW-MLR表示基于逐步算法的多元線性回歸模型;R2表示決定系數(shù);MAE表示平均絕對誤差。Fig. 1 Workflow diagram for quantitative structure-activity relationship (QSAR) models construction and acute toxicity semi-lethal dose (LD50) prediction of organophosphorus compounds (OPs)Note: GA-MLR represents the multiple linear regression model based on genetic algorithm; SW-MLR represents the multiple linear regression model based on stepwise algorithm; R2 represents the coefficient of determination; MAE represents mean absolute error.

    1.4 模型評估

    1.5 模型應用域

    使用Williams圖(又稱杠桿方法)[30]對建立的QSAR模型進行應用域分析。基于帽子矩陣(一類投影矩陣)得出的值(h)與標準化殘差(δ)的比較。h大于警告杠桿值(h*)的化合物被認為是結(jié)構(gòu)異常值,標準化殘差>2.5的化合物則被認為是計算異常值。Williams圖的公式定義如下:

    (1)

    表1 本研究中用于模型評價的公式及參數(shù)Table 1 Equations and parameters used for model evaluation in this study

    (2)

    (3)

    1.6 分子對接

    我們先前的實驗研究證實了稀有鮈鯽的丁酰膽堿酯酶(BChE)能夠被一些OPs顯著抑制,從而引起神經(jīng)毒性[31]。這項研究中,借助分子對接工具探究預測集中9種OPs對哺乳動物類的BChE的結(jié)合與抑制能力。從蛋白質(zhì)數(shù)據(jù)庫中檢索到BChE的人類蛋白質(zhì)結(jié)構(gòu)(PDB ID: 5K5E)作為對接蛋白。通過AutoDockTools4程序進行預處理,去除水分、添加氫并分配電荷。對接網(wǎng)格的大小設置為52 ?×52 ?×52 ?,并使用AutoDock Vina程序[32]進行分子對接。對接前,從BChE蛋白中提取了共結(jié)晶化合物并重新對接至原蛋白進行對接程序的驗證。最終使用PLIP工具[33]分析了對接后化合物與蛋白質(zhì)的結(jié)合模式。

    2 結(jié)果(Results)

    2.1 預測模型

    逐步算法的參數(shù)設置為進入F值:<0.010,逐出F值:>0.050。建立的SW-MLR模型如下:

    logLD50=-0.154(±1.5150)+0.008(±0.0018)ATS3s-2.783(±0.5247)GATS3e+2.054(±0.4083)MATS8i-88.853(±13.1459)JGI7-2.871(±0.6230)MATS3i+1.257(±0.2907)maxssO-(±0.0553)0.180topoDiameter

    遺傳算法的參數(shù)設置為最大迭代次數(shù):100,方程參數(shù)數(shù)量:7,初始方程數(shù):100,每代方程數(shù):30,適應函數(shù):MAE。建立的GA-MLR模型如下;

    logLD50=5.608(±0.7354)+0.0067(±0.0049)AATSC7m-0.012(±0.0088)ATSC3i+0.0133(±0.0087)VR3_Dzp-3.1038(±0.9454)AVP_4-85.9765(±14.6601)JGI7+0.263(±0.1021)nsssCH-1.2293(±0.3834)GATS4p

    對于每個模型還執(zhí)行了2 000次Y-隨機測試。隨機測試結(jié)果顯示SW-MLR和GA-MLR的最大偽R2值分別為0.488和0.504,最大偽Q2值分別為0.299和0.234。所有的偽R2和偽Q2都很低,這表明原始模型的結(jié)果并不是偶然關聯(lián)的。SW-MLR和GA-MLR模型的擬合關系結(jié)果如圖2所示。

    圖2 SW-MLR(a)和GA-MLR(b)模型中OPs實驗值與預測值擬合關系Fig. 2 Fitting relationship between experimental and predicted values of OPs in SW-MLR (a) and GA-MLR (b) models

    2.2 分子描述符解釋

    在SW-MLR方程中,分子描述符ATS3s、MATS8i和maxssO與毒性終點logLD50正相關。ATS3s、MATS8i分別是基于原子的內(nèi)在狀態(tài)和電離電位來表征化合物拓撲結(jié)構(gòu)的分子描述符[34]。maxssO則代表分子結(jié)構(gòu)中—O—基團的最大電性拓撲狀態(tài),前人的研究中也選用maxssO來描述一些無毒的抗瘧化合物[35],這與我們研究中maxssO描述符在SW-MLR方程中的正系數(shù)一致。分子描述符GATS3e、JGI7、MATS3i和topoDiameter則與毒性終點logLD50負相關。GATS3e和JGI7描述符與化合物中高電負性的原子相關[36-37]。MATS3i和topoDiameter是特定拓撲結(jié)構(gòu)的分子描述符。

    在GA-MLR方程中,分子描述符AATSC7m、nsssCH和VR3_Dzp與毒性終點logLD50正相關。AATSC7m和nsssCH分別對應OPs分子結(jié)構(gòu)中R基上支鏈的長度[38]、支鏈或取代基烷基化的程度[39]。VR3_Dzp則是比較抽象的二維描述符,代表極化率加權的Barysz矩陣的特征向量[40]。分子描述符ATSC3i、AVP_4、JGI7和GATS4p則與毒性終點logLD50負相關。這4個描述符與分子的疏水性、極性和原子的電負性相關[36,41]。對于2個模型分子描述符的解釋從一定程度上揭示了OPs的化學結(jié)構(gòu)與毒性相關的內(nèi)在因素,對于識別OPs的毒性起到關鍵作用。

    2.3 模型應用域

    建立的SW-MLR模型和GA-MLR模型應用域表征Willimas圖如圖3所示。SW-MLR模型中,所有OPs類化合物都很好地落入了應用域范圍內(nèi)。GA-MLR模型中,僅有一種化合物(dimethoate, CAS: 60-51-5,h=0.570>h*=0.558)為離群點。結(jié)論表明本研究的2種模型能夠較好地預測OPs類化合物對大鼠的急性毒性數(shù)據(jù)。表2中給出了使用2種QSAR模型預測的沒有實驗數(shù)據(jù)的9種OPs的急性毒性結(jié)果。

    圖3 SW-MLR(a)和GA-MLR(b)的Williams應用域圖Fig. 3 Williams plotting for SW-MLR (a) and GA-MLR (b)

    2.4 模型比較

    2.5 分子對接

    AutoDock Vina重新對接BChE(5K5E)的原配體的結(jié)合能結(jié)果為-11.9 kcal·mol-1,均方根偏差(RMSD)為0.89 ?,說明對接程序適用于進一步的對接與研究。使用AutoDock Vina軟件對接預測集(9種OPs)與BChE(5K5E)的最低對接能量結(jié)果如表2所示。PLIP工具對OPs與BChE的結(jié)合模式可視化如圖4所示。

    AutoDock Vina對接的結(jié)果顯示其中8個OPs的對接結(jié)合能低于閾值[43]-6 kcal·mol-1,表明了潛在的結(jié)合能力以及進一步可能造成的神經(jīng)毒性。此外,取代基或支鏈烷基化程度較高的化合物(如磷酸三(3-氯丙基)酯)對接結(jié)合能也較高,這表明更低的結(jié)合能力與更低的神經(jīng)毒性,這一結(jié)果與GA-MLR模型中分子描述符nsssCH的含義一致。因此我們認為OPs的R基及取代基的烷基化能夠降低毒性。由圖4可知,BChE與OPs結(jié)合的關鍵氨基酸殘基有Gly116和Ser198(氫鍵);His428、Trp332和Phe329(鹽橋);Trp231、Phe329、Trp82、Tyr332和His438(π-堆疊)。這些關鍵氨基酸殘基與前人的研究較為一致[44-45]。

    3 討論(Discussion)

    研究采用了SW和GA算法,基于PaDEL分子描述符,對OPs化合物的大鼠口服急性毒性建立了多元線性回歸QSAR模型,SW-MLR模型具有更好的擬合能力,GA-MLR模型有著更好的預測能力。對模型分子描述符的分析顯示了OPs的R基的高烷基化和支鏈的延長能夠降低毒性。此外,通過模型預測了數(shù)據(jù)集中沒有實驗數(shù)據(jù)的9種OPs的毒性,并以分子對接方法檢測其與BChE的潛在結(jié)合能力。所構(gòu)建模型能夠方便有效地預測OPs類化合物的毒性,并為新化學品的合成與環(huán)境保護政策的制定提供幫助。

    表2 SW-MLR和GA-MLR對9種OPs的LD50預測值、沒有實驗值的9種OPs與丁酰膽堿酯酶(BChE)對接的能量與關鍵殘基Table 2 Predicted LD50 values by SW-MLR and GA-MLR for nine OPs, docking energy and key residue results by docking to butyrylcholinesterase (BChE)

    圖4 9種OPs與BChE(5K5E)分子對接結(jié)果的可視化注:A~I對應表2中化合物序號;疏水性相互作用為灰色虛線,平行π-堆疊為亮綠色虛線,垂直π-堆疊為深綠色虛線,氫鍵為藍色實線,鹽橋為黃色虛線。Fig. 4 Ligand interaction diagram of the main interactions of nine OPs in the active site of BChE (5K5E)Note: A~I represents the compounds in Table 2; hydrophobic interactions are represented by gray dashes; π-stackings are represented by bright green (parallel) and dark green (perpendicular) dashes; hydrogen bonds are represented by blue solid lines; salt bridges are represented by yellow dashes.

    通訊作者簡介:查金苗(1975—),男,博士,研究員,主要研究方向為水生模型生物體系的構(gòu)建與發(fā)展、水環(huán)境生物毒性測試方法、環(huán)境內(nèi)分泌干擾物的篩選技術研究、環(huán)境污染物對水生生物分子毒理機制和水生態(tài)系統(tǒng)完整性評估方法等。

    共同通訊作者簡介:閆賽紅(1988—),女,博士,助理研究員,主要研究方向為水生態(tài)毒理學。

    猜你喜歡
    實驗模型研究
    一半模型
    記一次有趣的實驗
    FMS與YBT相關性的實證研究
    遼代千人邑研究述論
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    視錯覺在平面設計中的應用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    做個怪怪長實驗
    EMA伺服控制系統(tǒng)研究
    3D打印中的模型分割與打包
    中国美白少妇内射xxxbb| 欧美一区二区亚洲| 国产淫片久久久久久久久| 精品99又大又爽又粗少妇毛片| 久久久欧美国产精品| 日本黄色视频三级网站网址| 日韩三级伦理在线观看| 插阴视频在线观看视频| 国产色爽女视频免费观看| 久久久久国产网址| 男女下面进入的视频免费午夜| 美女xxoo啪啪120秒动态图| 人妻少妇偷人精品九色| 欧美最黄视频在线播放免费| 在线观看av片永久免费下载| 淫秽高清视频在线观看| 人妻系列 视频| 国产伦精品一区二区三区四那| 欧美区成人在线视频| 一区福利在线观看| 久久6这里有精品| 亚洲一区高清亚洲精品| 国内精品一区二区在线观看| 久久精品国产自在天天线| 少妇裸体淫交视频免费看高清| 国产麻豆成人av免费视频| 观看免费一级毛片| 日本-黄色视频高清免费观看| 久久婷婷人人爽人人干人人爱| 少妇裸体淫交视频免费看高清| 日韩一区二区三区影片| 麻豆成人av视频| 神马国产精品三级电影在线观看| 青青草视频在线视频观看| 国产中年淑女户外野战色| 级片在线观看| 伦精品一区二区三区| 午夜福利在线观看免费完整高清在 | 97人妻精品一区二区三区麻豆| 中文欧美无线码| 精品久久国产蜜桃| 国产精品一区www在线观看| 黄色欧美视频在线观看| 我的老师免费观看完整版| 国产精品人妻久久久影院| 久久久久九九精品影院| 啦啦啦观看免费观看视频高清| 日本黄色视频三级网站网址| 亚洲欧洲日产国产| 成年免费大片在线观看| 99久久九九国产精品国产免费| 国产一区二区三区av在线 | 老熟妇乱子伦视频在线观看| 亚洲综合色惰| 国产人妻一区二区三区在| 欧美最黄视频在线播放免费| 国产乱人视频| 超碰av人人做人人爽久久| 欧美色欧美亚洲另类二区| 在线国产一区二区在线| 91av网一区二区| 校园人妻丝袜中文字幕| 亚洲欧美精品综合久久99| 天天一区二区日本电影三级| 久久亚洲精品不卡| 一本一本综合久久| 少妇的逼水好多| 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 51国产日韩欧美| kizo精华| 在线观看一区二区三区| 卡戴珊不雅视频在线播放| 一卡2卡三卡四卡精品乱码亚洲| 少妇熟女aⅴ在线视频| 成人三级黄色视频| 性插视频无遮挡在线免费观看| 国产精品久久电影中文字幕| 三级国产精品欧美在线观看| 成人漫画全彩无遮挡| 夜夜夜夜夜久久久久| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 97热精品久久久久久| 精品99又大又爽又粗少妇毛片| 波野结衣二区三区在线| 99在线人妻在线中文字幕| 日韩欧美精品免费久久| 亚洲最大成人手机在线| 免费观看a级毛片全部| 日韩一本色道免费dvd| 日本与韩国留学比较| 男女做爰动态图高潮gif福利片| 少妇的逼好多水| 两个人视频免费观看高清| 国产片特级美女逼逼视频| 亚洲av不卡在线观看| 国产探花在线观看一区二区| 亚洲国产精品合色在线| 我的老师免费观看完整版| 亚洲美女视频黄频| 激情 狠狠 欧美| 午夜福利视频1000在线观看| 麻豆av噜噜一区二区三区| 十八禁国产超污无遮挡网站| 亚洲最大成人中文| 国产av在哪里看| 麻豆国产97在线/欧美| 69人妻影院| 国产精品国产高清国产av| 亚洲在线观看片| 久久久欧美国产精品| 夫妻性生交免费视频一级片| 日日啪夜夜撸| 爱豆传媒免费全集在线观看| 亚洲av熟女| 美女 人体艺术 gogo| 国产亚洲av片在线观看秒播厂 | www日本黄色视频网| 成人特级av手机在线观看| 国内精品宾馆在线| 亚洲不卡免费看| 在线国产一区二区在线| 久久婷婷人人爽人人干人人爱| 日本一二三区视频观看| 99精品在免费线老司机午夜| 日日啪夜夜撸| 国产精品久久久久久av不卡| 一级毛片aaaaaa免费看小| 国产精品久久久久久亚洲av鲁大| 天美传媒精品一区二区| 亚洲在线自拍视频| 美女大奶头视频| 男女视频在线观看网站免费| 99精品在免费线老司机午夜| 亚洲精品日韩av片在线观看| 男人和女人高潮做爰伦理| 国产激情偷乱视频一区二区| 男人舔女人下体高潮全视频| 午夜福利视频1000在线观看| 欧美+日韩+精品| 九九久久精品国产亚洲av麻豆| 春色校园在线视频观看| 国产午夜精品一二区理论片| www日本黄色视频网| 国产成人影院久久av| 99精品在免费线老司机午夜| 成人国产麻豆网| 久久久久九九精品影院| 国产精品人妻久久久久久| 日韩中字成人| 久久久久久久久中文| 2022亚洲国产成人精品| 一区二区三区免费毛片| 欧美日本视频| 国产精品国产三级国产av玫瑰| 乱系列少妇在线播放| 伦精品一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲av免费高清在线观看| 欧美变态另类bdsm刘玥| 九九久久精品国产亚洲av麻豆| 久久久久网色| 国产精品.久久久| 精品一区二区免费观看| 99国产精品一区二区蜜桃av| 亚洲内射少妇av| 九色成人免费人妻av| 亚洲欧美日韩高清专用| 欧美一级a爱片免费观看看| 免费av不卡在线播放| 婷婷六月久久综合丁香| 天天躁日日操中文字幕| 国产综合懂色| 亚洲成人久久性| 中文字幕熟女人妻在线| 国产精品伦人一区二区| 亚洲在久久综合| 亚洲精品成人久久久久久| 亚洲七黄色美女视频| 蜜桃久久精品国产亚洲av| 欧美日韩综合久久久久久| 99热网站在线观看| 欧美日韩乱码在线| 亚洲av免费在线观看| 日韩欧美一区二区三区在线观看| 免费搜索国产男女视频| 久久精品国产99精品国产亚洲性色| 日韩在线高清观看一区二区三区| 精品人妻偷拍中文字幕| 国产亚洲91精品色在线| 亚洲欧美清纯卡通| 亚洲图色成人| 婷婷精品国产亚洲av| 最新中文字幕久久久久| 色播亚洲综合网| 一卡2卡三卡四卡精品乱码亚洲| 看免费成人av毛片| 麻豆国产97在线/欧美| 特大巨黑吊av在线直播| 日韩欧美 国产精品| 搡老妇女老女人老熟妇| 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 99热这里只有是精品在线观看| 美女脱内裤让男人舔精品视频 | 在线免费十八禁| 久久久久国产网址| 有码 亚洲区| 非洲黑人性xxxx精品又粗又长| 免费一级毛片在线播放高清视频| 12—13女人毛片做爰片一| 精品国产三级普通话版| 一区二区三区四区激情视频 | 三级国产精品欧美在线观看| 在线免费观看的www视频| 97热精品久久久久久| 尤物成人国产欧美一区二区三区| 在线a可以看的网站| 亚洲内射少妇av| 乱系列少妇在线播放| 免费一级毛片在线播放高清视频| 国内少妇人妻偷人精品xxx网站| 有码 亚洲区| 午夜激情欧美在线| 看黄色毛片网站| 亚洲av中文av极速乱| 黄色视频,在线免费观看| 成年免费大片在线观看| 看免费成人av毛片| a级毛片免费高清观看在线播放| 国产久久久一区二区三区| 男女啪啪激烈高潮av片| 18禁裸乳无遮挡免费网站照片| 欧美日本视频| 国产精品三级大全| 十八禁国产超污无遮挡网站| 亚洲欧美成人综合另类久久久 | 成人欧美大片| 亚洲一级一片aⅴ在线观看| 久久精品久久久久久久性| 亚洲中文字幕一区二区三区有码在线看| 日韩成人av中文字幕在线观看| 亚洲人成网站高清观看| 精品一区二区免费观看| 中文字幕免费在线视频6| 亚洲精品日韩在线中文字幕 | 国产精品久久久久久亚洲av鲁大| 亚洲国产色片| 人体艺术视频欧美日本| 国产精品女同一区二区软件| 久久这里有精品视频免费| 国产伦精品一区二区三区视频9| 久久久久性生活片| 国产高潮美女av| 99久久无色码亚洲精品果冻| 亚洲av中文字字幕乱码综合| 又粗又硬又长又爽又黄的视频 | 1000部很黄的大片| 国产视频首页在线观看| 久久精品久久久久久噜噜老黄 | 99久久中文字幕三级久久日本| 一本久久精品| 又黄又爽又刺激的免费视频.| 身体一侧抽搐| 99在线人妻在线中文字幕| 狂野欧美白嫩少妇大欣赏| 97超碰精品成人国产| 国产精品无大码| 男女视频在线观看网站免费| 乱码一卡2卡4卡精品| 两个人的视频大全免费| 久久久久久大精品| 男人狂女人下面高潮的视频| 亚洲色图av天堂| 亚洲综合色惰| 国产片特级美女逼逼视频| 日韩欧美在线乱码| 国产极品天堂在线| 丝袜美腿在线中文| 亚洲人成网站高清观看| 日韩一本色道免费dvd| 性插视频无遮挡在线免费观看| 国产成人一区二区在线| 欧美性猛交╳xxx乱大交人| 国产亚洲av片在线观看秒播厂 | 三级毛片av免费| 国产精品人妻久久久久久| 国产一区二区三区av在线 | 亚洲av二区三区四区| 色综合亚洲欧美另类图片| 欧美色视频一区免费| 成年版毛片免费区| 亚洲内射少妇av| 青春草国产在线视频 | 在线观看一区二区三区| 亚洲精品久久国产高清桃花| 一级毛片aaaaaa免费看小| 亚洲自拍偷在线| 国产精品久久视频播放| 亚洲三级黄色毛片| 少妇人妻一区二区三区视频| 国产精品久久视频播放| 中文字幕av在线有码专区| 亚洲欧洲国产日韩| 麻豆乱淫一区二区| 国产极品天堂在线| 亚洲精品久久国产高清桃花| kizo精华| 国产精品人妻久久久影院| 男插女下体视频免费在线播放| 欧美成人精品欧美一级黄| 全区人妻精品视频| 国产av不卡久久| 久久久久性生活片| 网址你懂的国产日韩在线| 99在线人妻在线中文字幕| a级毛色黄片| 国产白丝娇喘喷水9色精品| 在线免费观看的www视频| 少妇人妻精品综合一区二区 | 在线观看一区二区三区| 一本一本综合久久| av在线亚洲专区| 国产成人91sexporn| 麻豆成人午夜福利视频| 麻豆成人av视频| 色吧在线观看| 久久人人爽人人片av| 成年版毛片免费区| 久久久久久国产a免费观看| 成熟少妇高潮喷水视频| 国产午夜福利久久久久久| 六月丁香七月| 免费av不卡在线播放| 色播亚洲综合网| 91狼人影院| 日日撸夜夜添| 午夜激情福利司机影院| 欧美不卡视频在线免费观看| 国产一区亚洲一区在线观看| 韩国av在线不卡| 一级黄色大片毛片| 色哟哟·www| 91精品国产九色| 寂寞人妻少妇视频99o| 综合色丁香网| 日本一本二区三区精品| 亚洲欧美精品专区久久| 国产精华一区二区三区| 久久久久免费精品人妻一区二区| 久久精品国产清高在天天线| 久久精品国产亚洲网站| h日本视频在线播放| 丰满乱子伦码专区| 一本精品99久久精品77| 久久99蜜桃精品久久| 亚洲四区av| 伦理电影大哥的女人| 国产一区二区激情短视频| 久久久精品大字幕| 有码 亚洲区| 美女被艹到高潮喷水动态| 变态另类丝袜制服| 日韩视频在线欧美| 亚洲国产欧美人成| 午夜免费男女啪啪视频观看| 国产一区二区三区在线臀色熟女| 国产精品99久久久久久久久| 女人十人毛片免费观看3o分钟| 久久久久国产网址| 亚洲五月天丁香| 成年版毛片免费区| 一级毛片我不卡| 中文亚洲av片在线观看爽| 欧美xxxx黑人xx丫x性爽| h日本视频在线播放| 夜夜爽天天搞| 看十八女毛片水多多多| 99riav亚洲国产免费| 99视频精品全部免费 在线| 一区二区三区免费毛片| 国产午夜福利久久久久久| 亚洲av免费高清在线观看| 国产精品久久久久久久电影| 亚洲av免费高清在线观看| 身体一侧抽搐| 国产精华一区二区三区| 国产精品.久久久| 观看美女的网站| 久久久久免费精品人妻一区二区| 啦啦啦韩国在线观看视频| av女优亚洲男人天堂| 国产精品野战在线观看| 亚洲美女视频黄频| 成年免费大片在线观看| 插逼视频在线观看| 久久人人爽人人爽人人片va| 国产视频内射| 三级国产精品欧美在线观看| 18禁在线无遮挡免费观看视频| 国产极品天堂在线| 黑人高潮一二区| 又爽又黄无遮挡网站| 国产免费一级a男人的天堂| 18禁在线无遮挡免费观看视频| 老熟妇乱子伦视频在线观看| 久久久久九九精品影院| 一进一出抽搐gif免费好疼| 最近2019中文字幕mv第一页| 欧美日韩在线观看h| 亚洲精品456在线播放app| 欧美日韩精品成人综合77777| 女的被弄到高潮叫床怎么办| 日韩制服骚丝袜av| 日本一本二区三区精品| 国产精品日韩av在线免费观看| 欧美高清成人免费视频www| 成人漫画全彩无遮挡| 欧美激情国产日韩精品一区| 免费搜索国产男女视频| 精品熟女少妇av免费看| 最近的中文字幕免费完整| 最好的美女福利视频网| 精品久久久久久成人av| 黑人高潮一二区| 99久久精品一区二区三区| 国产视频首页在线观看| 人妻制服诱惑在线中文字幕| 嫩草影院入口| av在线蜜桃| 国产美女午夜福利| 精品人妻视频免费看| 校园春色视频在线观看| 99热只有精品国产| a级一级毛片免费在线观看| 亚洲最大成人av| 伊人久久精品亚洲午夜| 好男人视频免费观看在线| 国内精品美女久久久久久| 麻豆国产av国片精品| 91午夜精品亚洲一区二区三区| 国产精品一二三区在线看| 18禁黄网站禁片免费观看直播| 国内久久婷婷六月综合欲色啪| 日韩三级伦理在线观看| 中文字幕av在线有码专区| 观看美女的网站| 波野结衣二区三区在线| 国产片特级美女逼逼视频| 亚洲精品久久久久久婷婷小说 | 可以在线观看的亚洲视频| 只有这里有精品99| 伦精品一区二区三区| 最近视频中文字幕2019在线8| 国产精品乱码一区二三区的特点| 国产精品久久视频播放| 亚洲电影在线观看av| 亚洲精品久久国产高清桃花| 国产精品.久久久| av.在线天堂| 日韩强制内射视频| 男女那种视频在线观看| 日韩欧美一区二区三区在线观看| 亚洲人成网站在线播| 自拍偷自拍亚洲精品老妇| 99热6这里只有精品| 又粗又爽又猛毛片免费看| 精品欧美国产一区二区三| 小说图片视频综合网站| 性插视频无遮挡在线免费观看| 大香蕉久久网| 日日干狠狠操夜夜爽| 不卡视频在线观看欧美| 免费大片18禁| 天堂av国产一区二区熟女人妻| 国产精品一区二区三区四区久久| 国产av一区在线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产色片| 亚洲中文字幕日韩| 日韩欧美精品免费久久| 国产精品.久久久| av在线天堂中文字幕| 国产探花在线观看一区二区| 国产精品,欧美在线| 久久精品国产亚洲av天美| 深夜a级毛片| 夜夜夜夜夜久久久久| 亚洲av中文av极速乱| 久99久视频精品免费| 天堂中文最新版在线下载 | 国产精品一及| 久久人妻av系列| 午夜a级毛片| 99久国产av精品| 亚洲欧美日韩无卡精品| 中文字幕人妻熟人妻熟丝袜美| 又粗又硬又长又爽又黄的视频 | 亚洲国产欧洲综合997久久,| 国产成人一区二区在线| 久久久久九九精品影院| 亚洲真实伦在线观看| 日本在线视频免费播放| 欧美一区二区国产精品久久精品| 白带黄色成豆腐渣| 久99久视频精品免费| 亚洲精品成人久久久久久| av在线老鸭窝| 黄色视频,在线免费观看| av天堂中文字幕网| 非洲黑人性xxxx精品又粗又长| 亚洲欧洲国产日韩| 青春草视频在线免费观看| 日本av手机在线免费观看| 免费观看在线日韩| 婷婷色av中文字幕| 亚洲自偷自拍三级| 精品少妇黑人巨大在线播放 | 狂野欧美白嫩少妇大欣赏| 亚洲精品日韩在线中文字幕 | 国产精品人妻久久久影院| 久久久精品大字幕| 国产男人的电影天堂91| 日韩视频在线欧美| 亚洲熟妇中文字幕五十中出| 看片在线看免费视频| 亚洲精品粉嫩美女一区| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 国产av一区在线观看免费| 99国产极品粉嫩在线观看| 成人鲁丝片一二三区免费| 自拍偷自拍亚洲精品老妇| 少妇人妻精品综合一区二区 | 久久久久久大精品| 赤兔流量卡办理| 国产蜜桃级精品一区二区三区| 亚洲av不卡在线观看| 日韩精品有码人妻一区| 日本三级黄在线观看| 欧美bdsm另类| 中文字幕免费在线视频6| 26uuu在线亚洲综合色| 国产精品久久电影中文字幕| 亚洲欧美成人综合另类久久久 | 麻豆精品久久久久久蜜桃| 麻豆久久精品国产亚洲av| 亚洲国产精品成人久久小说 | 日本欧美国产在线视频| 亚洲国产高清在线一区二区三| 日韩亚洲欧美综合| 一进一出抽搐gif免费好疼| 热99re8久久精品国产| 国产av在哪里看| 蜜桃亚洲精品一区二区三区| 26uuu在线亚洲综合色| 国内久久婷婷六月综合欲色啪| 国产精品一及| 免费观看精品视频网站| 亚洲国产欧美人成| 三级经典国产精品| 99热全是精品| 久久久久久久久久黄片| 久久久久久久久久成人| 亚洲av免费高清在线观看| 黄片wwwwww| 国产激情偷乱视频一区二区| 特级一级黄色大片| 嫩草影院入口| 人人妻人人澡欧美一区二区| 亚洲av成人精品一区久久| 一区二区三区四区激情视频 | 亚洲一区高清亚洲精品| 中文精品一卡2卡3卡4更新| 国产精品一区二区三区四区免费观看| 国产黄色视频一区二区在线观看 | 午夜免费激情av| 亚洲国产精品sss在线观看| 99久久人妻综合| av又黄又爽大尺度在线免费看 | 内地一区二区视频在线| 99久久精品热视频| 国产色婷婷99| 欧美zozozo另类| 22中文网久久字幕| 精品久久久久久久久久久久久| 黄色视频,在线免费观看| 亚洲第一区二区三区不卡| 午夜福利成人在线免费观看| 麻豆国产97在线/欧美| 欧美日韩综合久久久久久| 啦啦啦啦在线视频资源| 九色成人免费人妻av| av天堂中文字幕网| 美女高潮的动态| 哪里可以看免费的av片| 国产精品一区二区在线观看99 | 国产亚洲91精品色在线| 人人妻人人看人人澡| 日日干狠狠操夜夜爽| 免费av不卡在线播放| 亚洲四区av| 国产极品天堂在线| 日本黄色片子视频| www.色视频.com| 亚洲欧美日韩高清在线视频| 欧美一区二区国产精品久久精品| 欧美+日韩+精品| 最后的刺客免费高清国语| 1000部很黄的大片| 女人被狂操c到高潮| 久久精品国产亚洲av香蕉五月| 大又大粗又爽又黄少妇毛片口| 听说在线观看完整版免费高清|