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

    部分含氯芳烴結(jié)構(gòu)與孔雀魚(yú)半數(shù)致死濃度(-logLC50)定量活性相關(guān)研究

    2017-10-13 04:01:44李建鳳黃茜雷光東
    生態(tài)毒理學(xué)報(bào) 2017年3期
    關(guān)鍵詞:含氯描述符芳烴

    李建鳳,黃茜,雷光東

    內(nèi)江師范學(xué)院 化學(xué)化工學(xué)院,內(nèi)江 641100

    部分含氯芳烴結(jié)構(gòu)與孔雀魚(yú)半數(shù)致死濃度(-logLC50)定量活性相關(guān)研究

    李建鳳,黃茜,雷光東*

    內(nèi)江師范學(xué)院 化學(xué)化工學(xué)院,內(nèi)江 641100

    構(gòu)建了不同類(lèi)型的分子頂點(diǎn)之間的電性關(guān)系作為結(jié)構(gòu)描述符,對(duì)部分含氯芳烴化合物結(jié)構(gòu)進(jìn)行了參數(shù)化表征,共得到7個(gè)與化合物結(jié)構(gòu)密切相關(guān)的結(jié)構(gòu)描述符。逐步回歸(SMR)篩選變量后,分別運(yùn)用多元線(xiàn)性回歸(MLR)和偏最小二乘回歸(PLS)建立了化合物結(jié)構(gòu)與孔雀魚(yú)半數(shù)致死濃度(-logLC50)之間的關(guān)系模型,兩模型建模相關(guān)系數(shù)(r2)分別為0.871、0.862;“留一法”交互檢驗(yàn)的相關(guān)系數(shù)(Q2)分別為0.808、0.589。結(jié)果表明分子結(jié)構(gòu)描述符能恰當(dāng)?shù)乇碚骰衔锝Y(jié)構(gòu)特征,所建模型具有良好的穩(wěn)定性和預(yù)測(cè)能力。

    含氯芳烴化合物;結(jié)構(gòu)描述符;半數(shù)致死濃度;結(jié)構(gòu)與性質(zhì)的關(guān)系

    Received14 November 2016accepted7 February 2017

    Abstract: The electrical relationship between the different types of molecular vertices was constructed as structural descriptors. Then, structures of some chlorinated aromatic compounds were characterized. Seven structural descriptors, which were closely related to the structure of the compounds, were obtained. Two models of the relationship between the structure and half lethal concentration (-logLC50) to guppy were established through multiple linear regression (MLR) and partial least squares regression (PLS) methods. The correlation coefficients (r2) of the two models were 0.871 and 0.862. The cross-validation correlation coefficients (Q2) of "leave one out method" were 0.808 and 0.589, respectively. The results showed that the molecular structural descriptors could properly characterize the structural characteristics of the compounds. The models had good stability and prediction ability.

    Keywords: chlorinated aromatic compounds; structural descriptors; half lethal concentration; structure-property relationship

    含氯芳烴化合物在化工、醫(yī)藥等行業(yè)被廣泛用作原料或溶劑,使用過(guò)程中易泄漏進(jìn)入環(huán)境而成為重要的污染物。含氯芳烴化合物進(jìn)入環(huán)境后危害動(dòng)植物的生長(zhǎng)、繁殖,通過(guò)食物鏈的傳遞最終對(duì)人體健康造成威脅。含氯芳烴化合物進(jìn)入人體后,損害人體肝臟和腎臟,甚至產(chǎn)生致癌的潛在風(fēng)險(xiǎn)。研究含氯芳烴化合物的環(huán)境毒性,對(duì)于規(guī)范其生產(chǎn)、應(yīng)用和排放具有重要意義。測(cè)定每一個(gè)化合物的環(huán)境參數(shù),工作量大,而且還有源源不斷產(chǎn)生的新的化合物。利用構(gòu)效關(guān)系估算化合物各種參數(shù)成為環(huán)境化學(xué)中一個(gè)重要的分支,例如在有機(jī)污染物的毒性[1-2]、分配性[3-5]、溶解性[6-7]、降解性[8-9]、富集特性[10-11]等方面都取得了較好的成果。本研究選取了部分含氯芳烴化合物為研究樣本,通過(guò)構(gòu)建的分子頂點(diǎn)電性關(guān)系描述符對(duì)化合物結(jié)構(gòu)進(jìn)行表征,運(yùn)用多元線(xiàn)性回歸(multiple linear regression, MLR)和偏最小二乘回歸(partial least squares regression)的方法建立化合物結(jié)構(gòu)與毒性的關(guān)系模型,為芳烴類(lèi)化合物的QSAR研究提供參考。

    1 材料與方法(Materials and methods)

    1.1 實(shí)驗(yàn)材料

    選取37個(gè)含氯芳烴化合物為研究樣本,化合物急性毒性以其對(duì)孔雀魚(yú)半數(shù)致死濃度負(fù)對(duì)數(shù)(-logLC50)表示,化合物及其對(duì)孔雀魚(yú)急性毒性值(-logLC50)取自文獻(xiàn)[12],列于表1。

    1.2 實(shí)驗(yàn)方法1.2.1 化合物分子結(jié)構(gòu)參數(shù)化表達(dá)

    在有機(jī)化合物的骨架圖中,每一個(gè)非氫原子可視為分子頂點(diǎn),認(rèn)為分子頂點(diǎn)之間的電性關(guān)系對(duì)化合物外在性質(zhì)產(chǎn)生影響。處于不同連接狀態(tài)的分子頂點(diǎn)對(duì)化合物外在性質(zhì)產(chǎn)生的影響可能不同,因而需要對(duì)化合物中的非氫原子(即分子頂點(diǎn))進(jìn)行分類(lèi)。分子中的氫原子對(duì)化合物的外在性質(zhì)的影響可以忽略,氫原子僅僅影響與其直接相連的分子頂點(diǎn)的分類(lèi)。參照文獻(xiàn)[13-14],化合物中的不同分子頂點(diǎn)按照式(1)分為4類(lèi)。

    k = 4 - hi(hi=1,2,3)

    (1)

    其中,k為非氫原子i所屬分子頂點(diǎn)類(lèi)型,hi為與非氫原子i直接相連的氫原子數(shù)。當(dāng)hi分別為3、2、1、0時(shí), 非氫原子i分別屬于第1、2、3、4類(lèi)分子頂點(diǎn)。

    分子頂點(diǎn)之間的關(guān)系以分子頂點(diǎn)之間的電性相關(guān)性表示,一般認(rèn)為2個(gè)分子頂點(diǎn)之間的相關(guān)程度隨著兩者的電負(fù)性增減呈正向變化、隨著兩者的距離的增減呈反向變化。只要滿(mǎn)足上述條件的函數(shù)關(guān)系式均可用于構(gòu)建分子頂點(diǎn)之間的電性關(guān)系,文獻(xiàn)[13-14]中是以倒數(shù)形距離關(guān)系函數(shù)計(jì)算分子頂點(diǎn)之間的電性關(guān)系,這里采用Gaussian形距離關(guān)系函數(shù)式(2)進(jìn)行計(jì)算。

    (2)

    qi、qj為分子頂點(diǎn)i、j的相對(duì)電負(fù)性(該原子的電負(fù)性與碳原子電負(fù)性之比值);dij為關(guān)系中的分子頂點(diǎn)i到分子頂點(diǎn)j的相對(duì)距離(兩者間最短路徑所經(jīng)化學(xué)鍵鍵長(zhǎng)之和/碳碳單鍵鍵長(zhǎng),如有多條路徑,則以最短的為準(zhǔn));n和l為分子頂點(diǎn)所屬類(lèi)型?;衔镏?類(lèi)分子頂點(diǎn)可以組合出10種關(guān)系項(xiàng):m11,m12,…,m44,簡(jiǎn)寫(xiě)為x1,x2,…,x10,如m13(即x3)表示第1類(lèi)分子頂點(diǎn)跟第3類(lèi)分子頂點(diǎn)之間的關(guān)系,依此類(lèi)推,這樣對(duì)于有機(jī)化合物最多將產(chǎn)生10個(gè)結(jié)構(gòu)描述符值來(lái)描述其結(jié)構(gòu)。

    1.2.2 QSAR建模與檢驗(yàn)

    首先采用逐步回歸(stepwise regression, SMR)依據(jù)變量顯著性篩選出變量組合,然后分別采用多元線(xiàn)性回歸(multiple linear regression, MLR)和偏最小二乘回歸(partial least squares regression)的方法建立模型,運(yùn)用“留一法”對(duì)模型穩(wěn)定性及預(yù)測(cè)能力進(jìn)行檢驗(yàn)。一個(gè)理想的模型一般需要滿(mǎn)足建模相關(guān)系數(shù)(r2)≥0.64、標(biāo)準(zhǔn)偏差(SD)與數(shù)值范圍的比例小于10%、交互檢驗(yàn)相關(guān)系數(shù)(Q2)≥0.50[15]。多元線(xiàn)性回歸模型要求樣本數(shù)(N)/變量數(shù)(n) ≥5;變量之間的共線(xiàn)性采用方差膨脹因子VIF進(jìn)行評(píng)價(jià),變量中最大VIF小于5,表明變量間沒(méi)有明顯共線(xiàn)性,方程可以接受。

    2 結(jié)果(Results)

    經(jīng)結(jié)構(gòu)表征得到化合物結(jié)構(gòu)描述符值,x1、x2、x5為全“0”項(xiàng),其余結(jié)構(gòu)描述符值列于表1。采用逐步回歸(SMR)篩選變量,以偏F檢驗(yàn)值對(duì)應(yīng)的顯著水平值P為依據(jù),當(dāng)候選變量中最大偏F檢驗(yàn)值的P≤0.05,則引入相應(yīng)變量,在已進(jìn)入方程的變量中,若其最小偏F檢驗(yàn)值的P≥0.1,則剔除相應(yīng)變量,變量有進(jìn)有出,最終逐步回歸共進(jìn)行了7步運(yùn)算,逐步回歸結(jié)果見(jiàn)表2。

    表1 化合物及其急性毒性值(-logLC50)Table 1 The compounds and their acute toxicity values (-logLC50)

    注:x3為第1類(lèi)分子頂點(diǎn)與第3類(lèi)分子頂點(diǎn)之間的關(guān)系;x4為第1類(lèi)分子頂點(diǎn)與第4類(lèi)分子頂點(diǎn)之間的關(guān)系;x6為第2類(lèi)分子頂點(diǎn)與第3類(lèi)分子頂點(diǎn)之間的關(guān)系;x7為第2類(lèi)分子頂點(diǎn)與第4類(lèi)分子頂點(diǎn)之間的關(guān)系;x8為第3類(lèi)分子頂點(diǎn)與第3類(lèi)分子頂點(diǎn)之間的關(guān)系;x9為第3類(lèi)分子頂點(diǎn)與第4類(lèi)分子頂點(diǎn)之間的關(guān)系;x10為第4類(lèi)分子頂點(diǎn)與第4類(lèi)分子頂點(diǎn)之間的關(guān)系;-logLC50為急性毒性值;EXP.、MLR.、PLS.分別表示急性毒性的實(shí)驗(yàn)值、MLR模型預(yù)測(cè)值、PLS模型預(yù)測(cè)值。

    Note: x3is the relationship between the 1st type of molecular vertices and the 3rd type of molecular vertices; x4is the relationship between the 1st type of molecular vertices and the 4th type of molecular vertices; x6is the relationship between the 2nd type of molecular vertices and the 3rd type of molecular vertices; x7is the relationship between the 2nd type of molecular vertices and the 4th type of molecular vertices; x8is the relationship between the 3rd type of molecular vertices; x9is the relationship between the 3rd type of molecular vertices and the 4th type of molecular vertices; x10is the relationship between the 4th type of molecular vertices; -logLC50is the acute toxicity value; EXP, MLR and PLS are the experimental value, MLR model predicted value and PLS model predicted value, respectively.

    由表2可以看出,當(dāng)逐步回歸到第6步時(shí)(此時(shí)共篩選出6個(gè)變量),建模相關(guān)系數(shù)(r2)達(dá)到最大值。當(dāng)逐步回歸到第7步時(shí)(此時(shí)共篩選出5個(gè)變量),建模相關(guān)系數(shù)(r2)與第6步相比有所降低,但變量數(shù)減少,模型復(fù)雜程度降低了,調(diào)整相關(guān)系數(shù)(r2)達(dá)到最大值,交互檢驗(yàn)的相關(guān)系數(shù)(Q2)達(dá)到最大值,標(biāo)準(zhǔn)偏差(SD)達(dá)到最小值,而此時(shí)變量中的最大VIF僅為2.408。綜合各方面考慮,應(yīng)該選擇第7步所篩選出來(lái)的變量組合進(jìn)行建模,多元線(xiàn)性回歸(MLR)模型如式(3)。

    -logLC50= -0.980+1.743×x3+1.577×x6+0.710×x8+0.262×x9+0.088×x10

    (3)

    N=37,r2=0.871,SD=0.185,F(xiàn)=41.798;Q2=0.808,F(xiàn)CV=26.035

    N為建模樣本數(shù),n為變量數(shù),r2為建模相關(guān)系數(shù),SD為標(biāo)準(zhǔn)偏差,F(xiàn)為顯著性檢驗(yàn)值;Q2為交互檢驗(yàn)的相關(guān)系數(shù),F(xiàn)CV為交互檢驗(yàn)的顯著性檢驗(yàn)值。上述模型符合樣本數(shù)(N)/變量數(shù)(n)≥5的經(jīng)驗(yàn)規(guī)則,建模相關(guān)系數(shù)(r2)達(dá)0.871,大于等于0.64,說(shuō)明模型擬合效果好;交互檢驗(yàn)的相關(guān)系數(shù)(Q2)達(dá)0.808,大于等于0.50,說(shuō)明模型具有良好的穩(wěn)健性和預(yù)測(cè)能力;標(biāo)準(zhǔn)偏差(SD)為0.185,而樣本急性毒性值(-logLC50)的數(shù)值范圍為2.12(最大值2.31-最小值0.19),0.185/2.12=8.73%,小于10%的標(biāo)準(zhǔn),說(shuō)明模型預(yù)測(cè)準(zhǔn)確性良好。由于模型中的變量x3和x6對(duì)于大部分樣本都為“0”,這樣可能導(dǎo)致模型結(jié)果存在一定的偶然性,下面用偏最小二乘回歸(PLS)建模,進(jìn)一步挖掘結(jié)構(gòu)描述符與化合物毒性之間的關(guān)系。

    以逐步回歸到第7步時(shí)篩選出5個(gè)結(jié)構(gòu)描述符為自變量X,化合物對(duì)孔雀魚(yú)急性毒性值(-logLC50)為因變量Y,建立偏最小二乘回歸(PLS)模型。建模相關(guān)系數(shù)(r2)及交互檢驗(yàn)相關(guān)系數(shù)(Q2)與建模主成分?jǐn)?shù)(A)的關(guān)系如圖1,由圖1可知應(yīng)該選擇3個(gè)主成分進(jìn)行建模,此時(shí)化合物及其對(duì)孔雀魚(yú)急性毒性值(-logLC50)與原始自變量回歸方程式為:-logLC50= -0.932+1.329×x3+1.257×x6+0.698×x8+0.265×x9+0.086×x10。此時(shí)建模相關(guān)系數(shù)(r2)為0.862,接近最大值,并且大于等于0.64的標(biāo)準(zhǔn);交互檢驗(yàn)的相關(guān)系數(shù)(Q2)達(dá)0.589,大于等于0.50的標(biāo)準(zhǔn);標(biāo)準(zhǔn)偏差(SD)為0.198,0.198/2.12=9.34%,小于10%的標(biāo)準(zhǔn),以上結(jié)果說(shuō)明PLS模型同樣具有良好的穩(wěn)定性和預(yù)測(cè)能力。

    圖1 r2/Q2隨主成分?jǐn)?shù)(A)變化情況Fig. 1 r2/Q2 changes with the principal components (A)

    表2 逐步回歸結(jié)果Table 2 Results of the stepwise regression

    圖2 樣本在前2個(gè)主成分得分分布Fig. 2 Compounds’ score distribution in the front two principal components

    圖3 偏最小二乘X空間規(guī)格化模型距離Fig. 3 Normalized distance to PLS model in X space

    圖2為37個(gè)樣本在PLS前2個(gè)主成分得分空間散點(diǎn)分布圖,所有的樣本點(diǎn)都落在95%置信度Hotelling T2橢圓置信圈內(nèi),沒(méi)有異常點(diǎn),說(shuō)明結(jié)構(gòu)描述符能夠恰當(dāng)表現(xiàn)含氯芳烴化合物分子結(jié)構(gòu)特征,并在統(tǒng)計(jì)模型中得到正確反映。為了考察樣本在X空間的擬合情況,對(duì)樣本在X空間的規(guī)格化模型距離作了分析(圖3),同樣可以發(fā)現(xiàn),所有樣本規(guī)格化模型距離都在95%的置信區(qū)間內(nèi),小于臨界值2.135,同樣沒(méi)有出現(xiàn)異常點(diǎn),說(shuō)明模型質(zhì)量良好。

    圖4為PLS中變量重要性(VIP)投影圖,通常認(rèn)為VIP大于1的變量與Y之間具有較大的相關(guān)性。圖4中可以發(fā)現(xiàn)x8、x10的VIP值大于1,說(shuō)明這2個(gè)變量對(duì)含氯芳烴化合物的急性毒性值(-logLC50)影響較大,而變量x8對(duì)應(yīng)于第3類(lèi)分子頂點(diǎn)之間的關(guān)系,變量x10對(duì)應(yīng)于第4類(lèi)分子頂點(diǎn)之間的關(guān)系,因而第3類(lèi)原子(與1個(gè)H相連)、第4類(lèi)原子(與0個(gè)H相連)對(duì)含氯芳烴化合物的急性毒性值(-logLC50)影響顯著,以上說(shuō)明取代基種類(lèi)、數(shù)量及分布情況對(duì)化合物的急性毒性值(-logLC50)有重要影響。取代基數(shù)量及分布情況與分子體積有關(guān),取代基數(shù)量越多,并且集中分布在分子的某一區(qū)域,可能會(huì)造成分子體積或局部體積大,分子難以通過(guò)脂質(zhì)膜而產(chǎn)生毒性,表現(xiàn)出毒性較小。例如1號(hào)化合物2,3,4,5-四氯苯胺與4號(hào)化合物2,3,4-三氯苯胺相比,取代基更多,分子體積更大,因而1號(hào)化合物更難通過(guò)脂質(zhì)膜而產(chǎn)生毒性,表現(xiàn)出1號(hào)化合物毒性比4號(hào)化合物毒性更小。

    圖4 變量重要性投影Fig. 4 Importance of variables

    兩模型對(duì)樣本的急性毒性值(-logLC50)進(jìn)行了預(yù)測(cè),預(yù)測(cè)值分別列于表1的MLR.列和PLS.列。圖5為模型預(yù)測(cè)值與實(shí)驗(yàn)值之間的相關(guān)圖,圖6為預(yù)測(cè)誤差分布圖(Err.1、Err.2分別對(duì)應(yīng)于MLR.、PLS.模型)。從圖5可以看出,絕大部分樣本點(diǎn)都分布在45°平分線(xiàn)附近,說(shuō)明總體預(yù)測(cè)效果好;圖5中還可以發(fā)現(xiàn)PLS.預(yù)測(cè)的樣本點(diǎn)較PLR.預(yù)測(cè)的樣本點(diǎn)略顯分散,同時(shí)圖6中也容易發(fā)現(xiàn)由Err.1組成的曲線(xiàn)略顯靠?jī)?nèi),Err.2組成的曲線(xiàn)略顯靠外,說(shuō)明MLR模型比PLS模型更為優(yōu)良。同樣對(duì)比兩模型的相關(guān)系數(shù)(r2/Q2),發(fā)現(xiàn)MLR模型相關(guān)系數(shù)(r2/Q2)比PLS模型大;對(duì)比兩模型的標(biāo)準(zhǔn)偏差(SD),發(fā)現(xiàn)MLR模型標(biāo)準(zhǔn)偏差(SD)比PLS模型小,這些也表明了MLR模型的預(yù)測(cè)準(zhǔn)確性?xún)?yōu)于PLS模型??紤]到化合物的急性毒性值(-logLC50)本身具有一定的誤差,因而兩模型預(yù)測(cè)結(jié)果都是可以接受的。

    圖5 計(jì)算值與實(shí)驗(yàn)值相關(guān)圖Fig. 5 Correlation between calculated and experimental values

    圖6 計(jì)算誤差分布Fig. 6 Distribution of calculated errors

    3 討論(Discussion)

    將化合物分子隱氫圖中的分子頂點(diǎn)按照其所連的氫原子數(shù)進(jìn)行了分類(lèi),進(jìn)而通過(guò)函數(shù)關(guān)系式構(gòu)建了不同類(lèi)型的分子頂點(diǎn)之間的電性關(guān)系作為結(jié)構(gòu)描述符。結(jié)構(gòu)描述符簡(jiǎn)單、易懂、計(jì)算方便,但是結(jié)構(gòu)描述符是基于分子二維平面結(jié)構(gòu)進(jìn)行計(jì)算得到的2D分子結(jié)構(gòu)描述符,不能反映有機(jī)化合物真實(shí)的三維立體結(jié)構(gòu)特征,無(wú)法區(qū)分順?lè)串悩?gòu)體、光學(xué)異構(gòu)體,因而有待改進(jìn)為基于分子三維立體結(jié)構(gòu)計(jì)算的3D分子結(jié)構(gòu)描述符。對(duì)部分含氯芳烴化合物結(jié)構(gòu)進(jìn)行了參數(shù)化表征,逐步回歸(SMR)篩選變量后分別運(yùn)用多元線(xiàn)性回歸(MLR)和偏最小二乘回歸(PLS)建立了化合物結(jié)構(gòu)-急性毒性關(guān)系模型,經(jīng)檢驗(yàn)?zāi)P途哂辛己玫姆€(wěn)定性和預(yù)測(cè)能力,相對(duì)而言多元線(xiàn)性回歸(MLR)模型較偏最小二乘回歸(PLS)模型更為優(yōu)良。模型在一定程度上揭示了影響化合物急性毒性的結(jié)構(gòu)因素,對(duì)于有機(jī)污染物的QSAR研究具有一定的參考價(jià)值。

    [1] 堵錫華, 王超. 神經(jīng)網(wǎng)絡(luò)法應(yīng)用于酚類(lèi)化合物對(duì)青海弧菌毒性的預(yù)測(cè)[J]. 生態(tài)毒理學(xué)報(bào), 2016, 11(4): 90-94

    Du X H, Wang C. Prediction of the toxicity of phenol derivatives to Vibrio qinghaiensis by neural network method [J]. Asian Journal of Ecotoxicology, 2016, 11(4): 90-94 (in Chinese)

    [2] Liao L M, Huang X, Li J. Structural characterization and acute toxicity simulation for nitroaromatic compounds [J]. Chinese Journal of Structural Chemistry, 2016, 35 (3): 449-456

    [3] 王甫洋, 張學(xué)勝, 劉輝. 鹵代苯化合物水溶解度及正辛醇/水分配系數(shù)的三維構(gòu)效關(guān)系研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2015, 35(4): 1032-1039

    Wang F Y, Zhang X S, Liu H. Three-dimensional quantitative structure property relationship (3D-QSPR) on aqueous solubility (-lgSw) and noctanol/water partition coefficient (lgKow) of halogeno-benzene [J]. Acta Scientiae Circumstantiae, 2015, 35(4): 1032-1039 (in Chinese)

    [4] 廖立敏, 李雪誼, 鄒寧, 等. 鹵代烷烴正辛醇/水分配系數(shù)的模擬[J]. 計(jì)算機(jī)與應(yīng)用化學(xué), 2014, 31(5): 619-622

    Liao L M, Li X Y, Zou N, et al. Octanol-water partition coefficients simulation for halogenated hydrocarbons [J]. Computers and Applied Chemistry, 2014, 31(5): 619-622 (in Chinese)

    [5] 廖立敏. 鹵代苯甲醚結(jié)構(gòu)表征及醇/水分配系數(shù)模擬[J]. 計(jì)算機(jī)與應(yīng)用化學(xué), 2015, 32(7): 889-892

    Liao L M. Structural characterization and octanol/water partition coefficient simulation for halogenated anisole compounds [J]. Computers and Applied Chemistry, 2015, 32(7): 889-892 (in Chinese)

    [6] 潘善飛, 胡桂香, 呂楊, 等. 離子液體中有機(jī)物溶解度的QSPR模型分析[J]. 物理化學(xué)學(xué)報(bào), 2010, 26(9): 2494-2502

    Pan S F, Hu G X, Lv Y, et al. QSPR model analysis on the solubility of organic compounds in ionic liquids [J]. Acta Physico-Chimica Sinica, 2010, 26(9): 2494-2502 (in Chinese)

    [7] 吳桂玲, 吳啟勛. 用分子連接性研究有機(jī)化合物溶解度的QSPR[J]. 武漢工程大學(xué)學(xué)報(bào), 2008, 30(4): 42-45

    Wu G L, Wu Q X. Study on QSPR of organiccompounds’solubility with molecular connectivity indexes [J]. Journal of Wuhan Institute of Technology, 2008, 30(4): 42-45 (in Chinese)

    [8] 程榮, 戚道鐸, 王建龍. 用納米鐵降解氯酚的分子結(jié)構(gòu)-性質(zhì)相關(guān)性[J]. 清華大學(xué)學(xué)報(bào): 自然科學(xué)版, 2010, 50(6): 873-876

    Cheng R, Qi D D, Wang J L. Relationship between molecular structure and the degradation of chlorophenols by iron nanoparticles [J]. Journal of Tsinghua University: Science and Technology, 2010, 50(6): 873-876 (in Chinese)

    [9] 秦良, 羅斯, 高樹(shù)梅, 等. 零價(jià)鐵降解氯代有機(jī)污染物的QSPR研究[J]. 環(huán)境化學(xué), 2009, 28(3): 400-403

    Qin L, Luo S, Gao S M, et al. The QSPR investigation of the degradation of chlorinated contamination by zero-valent iron [J]. Environmental Chemistry, 2009, 28(3): 400-403 (in Chinese)

    [10] 鄭玉婷, 喬顯亮, 楊先海, 等. 鹵代有機(jī)化合物生物富集因子的定量結(jié)構(gòu)-活性關(guān)系模型[J]. 生態(tài)毒理學(xué)報(bào), 2013, 8(5): 772-777

    Zheng Y T, Qiao X L, Yang X H, et al. Quantitative structure-activity relationship model for bioconcentration factors of halogenated organic compounds [J]. Asian Journal of Ecotoxicology, 2013, 8(5): 772-777 (in Chinese)

    [11] 馮長(zhǎng)君, 楊偉華, 沐來(lái)龍. 用分子形狀指數(shù)和苯環(huán)因子預(yù)測(cè)有機(jī)物在魚(yú)體中的生物富集因子[J]. 北京工業(yè)大學(xué)學(xué)報(bào), 2008, 34(6): 607-616

    Feng C J, Yang W H, Mu L L. Estimation and prediction of bio-concentration factors of nonionic organic chemicals in fish by molecular shap index and benzene ring factor [J]. Journal of Beijing University of Technology, 2008, 34(6): 607-616 (in Chinese)

    [12] Huang H, Wang X D, Dai X L, et al. Holographic quantitative structure-activity relationship for prediction acute toxicity of benzene derivative to the guppy (Peocilia reticulate) [J]. Journal of Environmental Sciences, 2004, 3(16): 423-427

    [13] 廖立敏, 梅虎, 李建鳳, 等. H-MEDV描述法對(duì)蒙椴樹(shù)葉揮發(fā)性組分的氣相色譜保留時(shí)間的預(yù)測(cè)[J]. 天然產(chǎn)物研究與開(kāi)發(fā), 2008, 20(1): 47-51

    Liao L M, Mei H, Li J F, et al. Prediction on retention times of components from essential oil of Tiliam ongolica leaves by H-MEDV [J]. Natural Product Research and Development, 2008, 20(1): 47-51 (in Chinese)

    [14] 周麗平, 李志良, 余瑜. HMEDV描述法對(duì)氨基喹啉類(lèi)抗瘧藥活性預(yù)測(cè)[J]. 重慶醫(yī)科大學(xué)學(xué)報(bào), 2005, 30(6): 794-798

    Zhou L P, Li Z L, Yu Y. Prediction on inhibitory activity of aminoquinolines on malaria by H-MEDV [J]. Journal of Chongqing Medical University, 2005, 30(6): 794-798 (in Chinese)

    [15] 顧云蘭, 陳鑫, 簡(jiǎn)美玲. 苯胺類(lèi)化合物結(jié)構(gòu)與毒性的密度泛函理論研究[J]. 化學(xué)研究與應(yīng)用, 2015, 27(8): 1139-1144

    Gu Y L, Cheng X, Jian M L. Study on the structure-toxicity relationship of aniline compounds by density functional theory [J]. Chemical Research and Application, 2015, 27(8): 1139-1144 (in Chinese)

    StudyontheCorrelationbetweenStructureandQuantitativeActivityofHalfLethalConcentration(-logLC50)toGuppyforSomeChlorinatedAromaticCompounds

    Li Jianfeng, Huang Xi, Lei Guangdong*

    College of Chemistry and Chemical Engineering, Neijiang Normal University, Neijiang 641100, China

    10.7524/AJE.1673-5897.20161114002

    2016-11-14錄用日期2017-02-07

    1673-5897(2017)3-732-07

    X171.5

    A

    雷光東(1959-),男,教授,主要從事有機(jī)物分子結(jié)構(gòu)與性質(zhì)的關(guān)系、分子設(shè)計(jì)及合成等方面研究,以第一作者或通訊作者在國(guó)內(nèi)外發(fā)表學(xué)術(shù)論文50余篇。

    四川省教育廳基金項(xiàng)目(14ZB0259);四川省科技廳應(yīng)用基礎(chǔ)項(xiàng)目(2015JY0077)

    李建鳳(1982-),女,碩士,實(shí)驗(yàn)師,研究方向?yàn)榉肿咏Y(jié)構(gòu)與性能,E-mail: afeng523@126.com

    *通訊作者(Corresponding author), E-mail: leigdnjtc@126.com

    李建鳳, 黃茜, 雷光東. 部分含氯芳烴結(jié)構(gòu)與孔雀魚(yú)半數(shù)致死濃度(-logLC50)定量活性相關(guān)研究[J]. 生態(tài)毒理學(xué)報(bào),2017, 12(3): 732-738

    Li J F, Huang X, Lei G D. Study on the correlation between structure and quantitative activity of half lethal concentration (-logLC50) to guppy for some chlorinated aromatic compounds [J]. Asian Journal of Ecotoxicology, 2017, 12(3): 732-738 (in Chinese)

    猜你喜歡
    含氯描述符芳烴
    2022年高考理綜化學(xué)模擬試題B卷參考答案
    基于結(jié)構(gòu)信息的異源遙感圖像局部特征描述符研究
    合理施用含氯化肥可提高獼猴桃產(chǎn)量和品質(zhì)
    關(guān)于重芳烴輕質(zhì)化與分離的若干思考
    含氯消毒劑,別僅憑名字辨別
    Linux單線(xiàn)程并發(fā)服務(wù)器探索
    利用CNN的無(wú)人機(jī)遙感影像特征描述符學(xué)習(xí)
    新型含氯(氟)的A2B型咔咯化合物的合成及其光學(xué)性質(zhì)
    輪胎填充油中8種多環(huán)芳烴的檢測(cè)
    高芳烴環(huán)保橡膠填充油量產(chǎn)
    日韩电影二区| 黄色配什么色好看| 丝瓜视频免费看黄片| 九色亚洲精品在线播放| av天堂久久9| 国产男女内射视频| 国产亚洲欧美精品永久| 亚洲欧洲国产日韩| 一区在线观看完整版| 亚洲天堂av无毛| 久久久久精品久久久久真实原创| 免费观看a级毛片全部| 亚洲精品,欧美精品| 久久99蜜桃精品久久| 久久精品亚洲av国产电影网| 这个男人来自地球电影免费观看 | 亚洲国产成人一精品久久久| 啦啦啦啦在线视频资源| 老司机影院成人| 亚洲精品国产av成人精品| 国精品久久久久久国模美| 美女福利国产在线| 久久久精品94久久精品| 色哟哟·www| 熟女少妇亚洲综合色aaa.| 成年女人毛片免费观看观看9 | 超色免费av| 一级毛片黄色毛片免费观看视频| 亚洲人成77777在线视频| 久久久国产精品麻豆| 黄色配什么色好看| 成人亚洲精品一区在线观看| 人人妻人人澡人人爽人人夜夜| 久久久久久久久久久免费av| 欧美日韩综合久久久久久| 女人久久www免费人成看片| www.av在线官网国产| 免费大片黄手机在线观看| 免费大片黄手机在线观看| 国产精品.久久久| 成人亚洲精品一区在线观看| 欧美日韩综合久久久久久| 亚洲成人av在线免费| 久久国产精品大桥未久av| 亚洲成人一二三区av| 老司机亚洲免费影院| 欧美国产精品va在线观看不卡| 99热全是精品| 韩国av在线不卡| 三级国产精品片| 超碰97精品在线观看| 久久久久久免费高清国产稀缺| 18禁动态无遮挡网站| 欧美亚洲日本最大视频资源| 亚洲精品国产av蜜桃| 国产精品久久久久久久久免| 妹子高潮喷水视频| 激情视频va一区二区三区| 女的被弄到高潮叫床怎么办| 美女视频免费永久观看网站| 日韩中文字幕欧美一区二区 | 国产成人精品久久二区二区91 | 天天躁狠狠躁夜夜躁狠狠躁| 精品视频人人做人人爽| 午夜免费男女啪啪视频观看| 日韩视频在线欧美| 男女午夜视频在线观看| 美女午夜性视频免费| 国产色婷婷99| 两个人免费观看高清视频| 午夜福利一区二区在线看| 国产又爽黄色视频| 成年动漫av网址| av在线老鸭窝| 国产午夜精品一二区理论片| 亚洲国产看品久久| 欧美变态另类bdsm刘玥| 欧美+日韩+精品| 在线精品无人区一区二区三| 亚洲男人天堂网一区| 伦理电影大哥的女人| 亚洲成人av在线免费| 中文字幕另类日韩欧美亚洲嫩草| 女性生殖器流出的白浆| 国产探花极品一区二区| 蜜桃在线观看..| 精品一品国产午夜福利视频| 久久精品aⅴ一区二区三区四区 | 香蕉丝袜av| 免费观看在线日韩| 一级,二级,三级黄色视频| 在线观看免费日韩欧美大片| 如日韩欧美国产精品一区二区三区| 国产精品成人在线| 91午夜精品亚洲一区二区三区| 黄片无遮挡物在线观看| a级毛片在线看网站| 成人黄色视频免费在线看| 波多野结衣一区麻豆| 精品国产一区二区久久| 十分钟在线观看高清视频www| 日韩大片免费观看网站| 亚洲精品国产一区二区精华液| 成人影院久久| 人人妻人人添人人爽欧美一区卜| 在线观看美女被高潮喷水网站| 久久99热这里只频精品6学生| 观看美女的网站| 亚洲欧洲精品一区二区精品久久久 | 少妇人妻 视频| 亚洲成国产人片在线观看| 久久久国产一区二区| 成人18禁高潮啪啪吃奶动态图| 日日摸夜夜添夜夜爱| 亚洲精品国产av成人精品| 一级a爱视频在线免费观看| 热re99久久精品国产66热6| 亚洲精品国产av蜜桃| 一个人免费看片子| 日韩精品免费视频一区二区三区| 午夜福利乱码中文字幕| 1024香蕉在线观看| www.av在线官网国产| 国产男女超爽视频在线观看| 欧美精品一区二区大全| 久久精品久久久久久噜噜老黄| 欧美另类一区| 91精品国产国语对白视频| 欧美精品亚洲一区二区| 久久久久久人人人人人| 亚洲精品久久成人aⅴ小说| 亚洲精品美女久久久久99蜜臀 | 麻豆乱淫一区二区| 99久久人妻综合| 欧美老熟妇乱子伦牲交| 90打野战视频偷拍视频| 亚洲一码二码三码区别大吗| 国产欧美日韩综合在线一区二区| 国产一区亚洲一区在线观看| 亚洲欧洲精品一区二区精品久久久 | 中文字幕人妻丝袜一区二区 | 国产白丝娇喘喷水9色精品| 亚洲精品视频女| 男女啪啪激烈高潮av片| av网站免费在线观看视频| 大香蕉久久成人网| 亚洲美女视频黄频| 国产精品三级大全| 午夜福利视频精品| 我要看黄色一级片免费的| 青春草国产在线视频| 老女人水多毛片| 国产午夜精品一二区理论片| 精品少妇黑人巨大在线播放| 久久久久久久久久人人人人人人| 久久久亚洲精品成人影院| 国产一区二区在线观看av| 久久久久人妻精品一区果冻| av在线播放精品| 一本大道久久a久久精品| 国产精品一二三区在线看| 亚洲少妇的诱惑av| 两个人免费观看高清视频| 免费高清在线观看日韩| 日产精品乱码卡一卡2卡三| 免费观看a级毛片全部| 中文精品一卡2卡3卡4更新| 亚洲精品久久午夜乱码| 曰老女人黄片| 乱人伦中国视频| 18+在线观看网站| 免费高清在线观看日韩| 国产精品偷伦视频观看了| 桃花免费在线播放| xxxhd国产人妻xxx| 男人爽女人下面视频在线观看| 成年av动漫网址| 午夜福利,免费看| 狠狠婷婷综合久久久久久88av| 观看av在线不卡| 欧美中文综合在线视频| 啦啦啦在线免费观看视频4| 亚洲国产精品成人久久小说| 国产深夜福利视频在线观看| 亚洲精华国产精华液的使用体验| 免费观看av网站的网址| 99热国产这里只有精品6| 免费av中文字幕在线| 亚洲精品国产av成人精品| 99热全是精品| 91成人精品电影| 91精品三级在线观看| 亚洲精品aⅴ在线观看| 亚洲av成人精品一二三区| 国产精品一区二区在线不卡| 女性被躁到高潮视频| 女的被弄到高潮叫床怎么办| 日本欧美视频一区| 久久女婷五月综合色啪小说| 久久99精品国语久久久| 成人国语在线视频| 啦啦啦视频在线资源免费观看| 亚洲一区中文字幕在线| 婷婷成人精品国产| h视频一区二区三区| 日韩熟女老妇一区二区性免费视频| 视频在线观看一区二区三区| 亚洲欧美成人综合另类久久久| 自拍欧美九色日韩亚洲蝌蚪91| 少妇猛男粗大的猛烈进出视频| 男人爽女人下面视频在线观看| 女的被弄到高潮叫床怎么办| 国产精品香港三级国产av潘金莲 | 国产一区亚洲一区在线观看| 欧美日韩视频高清一区二区三区二| 亚洲国产色片| 伊人久久国产一区二区| 日韩,欧美,国产一区二区三区| 老鸭窝网址在线观看| 极品人妻少妇av视频| 亚洲精品久久午夜乱码| 国产精品一二三区在线看| 成年动漫av网址| 欧美激情高清一区二区三区 | 亚洲av欧美aⅴ国产| 搡女人真爽免费视频火全软件| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产精品999| 欧美精品人与动牲交sv欧美| 国产97色在线日韩免费| 久久国内精品自在自线图片| 久久久久久久亚洲中文字幕| 国产精品av久久久久免费| 极品人妻少妇av视频| 一本大道久久a久久精品| 晚上一个人看的免费电影| 午夜免费男女啪啪视频观看| 国产成人精品久久久久久| 国产国语露脸激情在线看| 国产精品一区二区在线不卡| 国产精品熟女久久久久浪| 99久久人妻综合| 精品人妻熟女毛片av久久网站| 日韩人妻精品一区2区三区| 成人国产av品久久久| 日韩一区二区视频免费看| 成人毛片60女人毛片免费| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费视频播放在线视频| 欧美 日韩 精品 国产| 波多野结衣一区麻豆| 成人18禁高潮啪啪吃奶动态图| 色婷婷久久久亚洲欧美| videos熟女内射| 色网站视频免费| 亚洲国产av新网站| xxx大片免费视频| 91精品伊人久久大香线蕉| 最近手机中文字幕大全| 久久午夜综合久久蜜桃| 亚洲欧美精品自产自拍| 精品国产一区二区久久| 日本av手机在线免费观看| 一区二区三区激情视频| 赤兔流量卡办理| 亚洲精品视频女| 国产男女内射视频| a 毛片基地| av视频免费观看在线观看| 伊人亚洲综合成人网| 欧美精品av麻豆av| 国产精品国产三级国产专区5o| 18禁国产床啪视频网站| 久久久久久久久久久久大奶| xxx大片免费视频| 国产老妇伦熟女老妇高清| 男女无遮挡免费网站观看| 男女下面插进去视频免费观看| 国产亚洲午夜精品一区二区久久| 亚洲国产色片| 在线观看免费视频网站a站| 亚洲av中文av极速乱| 久久久久精品久久久久真实原创| av卡一久久| 精品人妻一区二区三区麻豆| 欧美日韩一区二区视频在线观看视频在线| av免费在线看不卡| 又黄又粗又硬又大视频| 午夜福利视频在线观看免费| 午夜福利在线观看免费完整高清在| 色吧在线观看| 精品久久久精品久久久| 免费黄频网站在线观看国产| 国产av码专区亚洲av| 国产97色在线日韩免费| 在线 av 中文字幕| 热99久久久久精品小说推荐| 交换朋友夫妻互换小说| av网站免费在线观看视频| 97在线视频观看| 最近的中文字幕免费完整| 18禁国产床啪视频网站| 汤姆久久久久久久影院中文字幕| 9色porny在线观看| 老汉色∧v一级毛片| 性高湖久久久久久久久免费观看| 精品一区二区三区四区五区乱码 | 如何舔出高潮| 人成视频在线观看免费观看| 亚洲av国产av综合av卡| 国产亚洲午夜精品一区二区久久| 亚洲三区欧美一区| 亚洲成av片中文字幕在线观看 | 天美传媒精品一区二区| 精品99又大又爽又粗少妇毛片| 欧美最新免费一区二区三区| 国产极品粉嫩免费观看在线| 亚洲精品av麻豆狂野| 亚洲国产成人一精品久久久| 成年人免费黄色播放视频| 一级毛片 在线播放| 在现免费观看毛片| 在线观看免费视频网站a站| 精品人妻偷拍中文字幕| 国产高清国产精品国产三级| 国产精品久久久久久久久免| 看十八女毛片水多多多| 日韩大片免费观看网站| 国产精品一二三区在线看| 大香蕉久久网| 亚洲综合色惰| 久久人人爽人人片av| 国产免费现黄频在线看| 夫妻午夜视频| 飞空精品影院首页| 涩涩av久久男人的天堂| 日本免费在线观看一区| 一边亲一边摸免费视频| 亚洲精品第二区| 亚洲国产日韩一区二区| 99久久中文字幕三级久久日本| 日本欧美国产在线视频| 亚洲人成网站在线观看播放| 亚洲熟女精品中文字幕| 毛片一级片免费看久久久久| 久久精品国产亚洲av天美| 国产福利在线免费观看视频| a级毛片黄视频| 国产亚洲av片在线观看秒播厂| 午夜福利在线观看免费完整高清在| av又黄又爽大尺度在线免费看| 十分钟在线观看高清视频www| 亚洲四区av| 国产av一区二区精品久久| 午夜日本视频在线| 高清不卡的av网站| 国产深夜福利视频在线观看| 天堂中文最新版在线下载| 欧美少妇被猛烈插入视频| 飞空精品影院首页| 最近中文字幕2019免费版| 夫妻午夜视频| 18禁动态无遮挡网站| 中文字幕人妻丝袜制服| 久久久国产一区二区| 少妇人妻精品综合一区二区| 精品酒店卫生间| 国产免费现黄频在线看| 日韩大片免费观看网站| 国产精品久久久久久精品电影小说| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩中文字幕视频在线看片| 国产精品免费视频内射| 亚洲国产日韩一区二区| 99热网站在线观看| 中文乱码字字幕精品一区二区三区| 尾随美女入室| 国产麻豆69| 夜夜骑夜夜射夜夜干| 少妇的逼水好多| 美女主播在线视频| 2021少妇久久久久久久久久久| 另类亚洲欧美激情| 男人操女人黄网站| 色播在线永久视频| 日本欧美视频一区| 只有这里有精品99| 久久国产亚洲av麻豆专区| 国产成人精品一,二区| 制服诱惑二区| 在线 av 中文字幕| 不卡视频在线观看欧美| 99久久中文字幕三级久久日本| 欧美成人午夜精品| videosex国产| 自拍欧美九色日韩亚洲蝌蚪91| 国产无遮挡羞羞视频在线观看| 免费女性裸体啪啪无遮挡网站| 丰满迷人的少妇在线观看| av在线app专区| 国产福利在线免费观看视频| 啦啦啦在线免费观看视频4| 一区福利在线观看| 少妇人妻 视频| 97精品久久久久久久久久精品| 波野结衣二区三区在线| 男女国产视频网站| 久久人人97超碰香蕉20202| 日韩av在线免费看完整版不卡| 卡戴珊不雅视频在线播放| 国产又爽黄色视频| 国产一区二区三区av在线| 中文字幕精品免费在线观看视频| 性高湖久久久久久久久免费观看| 黑人猛操日本美女一级片| 精品亚洲成国产av| 成人黄色视频免费在线看| 欧美日韩av久久| 伦理电影大哥的女人| 久久久久人妻精品一区果冻| 十八禁高潮呻吟视频| 大香蕉久久成人网| 国产精品嫩草影院av在线观看| av有码第一页| 国产精品三级大全| 极品少妇高潮喷水抽搐| 日本av免费视频播放| 人妻少妇偷人精品九色| 久久久国产精品麻豆| 国产高清国产精品国产三级| 国产又爽黄色视频| 久久久久精品久久久久真实原创| 国产伦理片在线播放av一区| 亚洲国产成人一精品久久久| 久久免费观看电影| 亚洲婷婷狠狠爱综合网| 宅男免费午夜| 亚洲内射少妇av| 亚洲第一av免费看| 成人国产麻豆网| 国产一级毛片在线| 汤姆久久久久久久影院中文字幕| 女性被躁到高潮视频| 久久女婷五月综合色啪小说| 最近中文字幕高清免费大全6| 亚洲视频免费观看视频| 大码成人一级视频| 久久精品国产a三级三级三级| 中文欧美无线码| 国产精品99久久99久久久不卡 | 一级毛片黄色毛片免费观看视频| 黄色怎么调成土黄色| 国产av一区二区精品久久| 久久精品国产自在天天线| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人aa在线观看| 中文字幕另类日韩欧美亚洲嫩草| 97在线视频观看| 狠狠精品人妻久久久久久综合| 91精品三级在线观看| 精品久久久久久电影网| 丝袜脚勾引网站| 欧美老熟妇乱子伦牲交| 777米奇影视久久| 一区二区三区激情视频| 亚洲成色77777| 久久久国产一区二区| av在线老鸭窝| 视频区图区小说| 黄色 视频免费看| 日韩av不卡免费在线播放| 欧美人与性动交α欧美软件| 宅男免费午夜| 大片电影免费在线观看免费| 国产精品一二三区在线看| videosex国产| 久久久国产一区二区| 97在线人人人人妻| 日韩av不卡免费在线播放| 国产成人a∨麻豆精品| 亚洲三级黄色毛片| 91aial.com中文字幕在线观看| 一级毛片我不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产无遮挡羞羞视频在线观看| 在线免费观看不下载黄p国产| 人人妻人人澡人人爽人人夜夜| xxxhd国产人妻xxx| 久久久a久久爽久久v久久| 国产成人精品无人区| 免费观看性生交大片5| 女人高潮潮喷娇喘18禁视频| 亚洲熟女精品中文字幕| 中文乱码字字幕精品一区二区三区| 国产av码专区亚洲av| 午夜免费男女啪啪视频观看| 亚洲综合色网址| 久久久久久久大尺度免费视频| 母亲3免费完整高清在线观看 | 综合色丁香网| av在线app专区| 制服诱惑二区| 国产精品久久久久久精品电影小说| 欧美+日韩+精品| av片东京热男人的天堂| 老汉色∧v一级毛片| 中文天堂在线官网| 我的亚洲天堂| 天天影视国产精品| 日产精品乱码卡一卡2卡三| 精品卡一卡二卡四卡免费| 久久久久国产网址| 亚洲内射少妇av| 欧美日韩亚洲高清精品| 91aial.com中文字幕在线观看| 麻豆乱淫一区二区| 欧美xxⅹ黑人| 欧美 亚洲 国产 日韩一| 伦精品一区二区三区| 国产成人精品婷婷| 欧美精品高潮呻吟av久久| 美女主播在线视频| 亚洲成人一二三区av| 免费观看在线日韩| 成人影院久久| 国产精品国产av在线观看| 免费在线观看黄色视频的| 18禁动态无遮挡网站| 伊人亚洲综合成人网| 欧美国产精品va在线观看不卡| 国产欧美亚洲国产| 国产成人av激情在线播放| 精品国产一区二区三区四区第35| 一级毛片 在线播放| 亚洲精品中文字幕在线视频| 国产黄频视频在线观看| 91精品国产国语对白视频| 18+在线观看网站| 一级毛片电影观看| 国产日韩欧美亚洲二区| 性少妇av在线| 久久久精品区二区三区| 超色免费av| 高清欧美精品videossex| 亚洲精品美女久久久久99蜜臀 | 国产精品熟女久久久久浪| 永久网站在线| 啦啦啦中文免费视频观看日本| 精品人妻在线不人妻| 欧美国产精品一级二级三级| 99久久人妻综合| 国产黄频视频在线观看| 麻豆精品久久久久久蜜桃| 天堂中文最新版在线下载| xxxhd国产人妻xxx| 日韩电影二区| 日韩av免费高清视频| 午夜免费观看性视频| 男女国产视频网站| 精品人妻在线不人妻| 精品一区在线观看国产| av视频免费观看在线观看| 97在线视频观看| kizo精华| 国产亚洲欧美精品永久| 天天影视国产精品| 成人国产麻豆网| 久久久精品区二区三区| 国产黄频视频在线观看| 欧美精品亚洲一区二区| 日本免费在线观看一区| 啦啦啦中文免费视频观看日本| 日韩电影二区| 亚洲av.av天堂| 婷婷色av中文字幕| 亚洲精品aⅴ在线观看| 久久久久国产精品人妻一区二区| 免费观看在线日韩| 91精品国产国语对白视频| 午夜福利,免费看| 日日爽夜夜爽网站| 亚洲成人av在线免费| 夜夜骑夜夜射夜夜干| 国产福利在线免费观看视频| 国产 精品1| 18在线观看网站| 老汉色∧v一级毛片| 亚洲国产av影院在线观看| 午夜福利视频在线观看免费| 国产精品免费大片| 两个人免费观看高清视频| 亚洲精华国产精华液的使用体验| 中文精品一卡2卡3卡4更新| 蜜桃在线观看..| 免费观看性生交大片5| 国产白丝娇喘喷水9色精品| 美国免费a级毛片| 男的添女的下面高潮视频| 97精品久久久久久久久久精品| 91aial.com中文字幕在线观看| 亚洲美女搞黄在线观看| 亚洲精品乱久久久久久| 午夜福利乱码中文字幕| 午夜免费鲁丝| 亚洲成人av在线免费| 男人爽女人下面视频在线观看| 精品人妻偷拍中文字幕| 久久久久久免费高清国产稀缺| 亚洲欧洲国产日韩| 亚洲国产毛片av蜜桃av| 男女高潮啪啪啪动态图| 日本黄色日本黄色录像| 国产免费一区二区三区四区乱码| 寂寞人妻少妇视频99o| 狂野欧美激情性bbbbbb|