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

    二苯甲酮類紫外防曬劑發(fā)光菌急性毒性及QSAR研究

    2017-10-13 04:02:16魏東斌趙慧敏杜宇國
    生態(tài)毒理學報 2017年3期
    關鍵詞:苯甲苯環(huán)甲氧基

    魏東斌,趙慧敏,杜宇國,2

    1. 中國科學院生態(tài)環(huán)境研究中心 環(huán)境化學與生態(tài)毒理學國家重點實驗室,北京1000852. 中國科學院大學,北京100049

    二苯甲酮類紫外防曬劑發(fā)光菌急性毒性及QSAR研究

    魏東斌1,2,*,趙慧敏1,杜宇國1,2

    1. 中國科學院生態(tài)環(huán)境研究中心 環(huán)境化學與生態(tài)毒理學國家重點實驗室,北京1000852. 中國科學院大學,北京100049

    二苯甲酮類化合物廣泛用于防曬劑、塑料添加劑、香味劑等。隨著紫外防曬產品的大量使用,其環(huán)境及健康風險越來越受到人們的關注。為揭示二苯甲酮類污染物的毒性特征,選擇了14種二苯甲酮類化合物作為目標化合物,測試了它們對發(fā)光細菌的急性毒性效應。分別運用二維、三維定量構效相關技術和分子對接技術探討了目標化合物的分子結構特征對毒性效應的影響。結果表明14種二苯甲酮類化合物對發(fā)光菌急性毒性的EC50值在17.67到243.82 mg·L-1范圍內,其中2-羥基-4-甲氧基-5-磺酸基二苯甲酮的急性毒性最低,2,2′,4,4′-四羥基二苯甲酮的毒性最高。羥基取代的二苯甲酮化合物的急性毒性隨著分子中羥基數(shù)量的增加而升高;具有相同羥基數(shù)量的二苯甲酮類化合物,羥基位于苯環(huán)4-位時毒性最高,3-位時次之,2-位時毒性最低。QSAR結果表明,運用靜電場、氫鍵受體場和氫鍵供體場能很好解釋這類化合物的毒性特征,若在苯環(huán)4-位引入帶正電荷的官能團、在苯環(huán)2-, 4-位引入氫鍵受體,都將導致毒性升高。上述研究結果將為科學評價該類化合物的潛在生態(tài)風險提供基礎數(shù)據(jù)。

    二苯甲酮;紫外防曬劑;發(fā)光細菌;急性毒性;定量構效相關

    Received23 January 2017accepted13 March 2017

    Abstract: Benzophenones (BPs) are widely used in UV filters, plastic additives, and fragrance enhancers. With the wide use of UV filters products, their potential risk to the ecological and human health has been concerned greatly. Fourteen BPs were selected as targets, and their acute toxicity was tested using the photobacterium assay. The related substituent contribution was disclosed by using 2- and 3-dimension quantitative structure-activity relationship (QSAR) and molecular docking methods. The results revealed that the EC50values of 14 BPs ranged from 17.67 to 243.82 mg·L-1, and the acute toxicity of 2-hydroxy-4-methoxy-5-sulfonic acid BP was the lowest, while that of 2,2’,4,4’-tetrahydroxy-BP was the highest. The acute toxicity increased with the increasing number of hydroxy group in BPs. For the BPs containing same number of hydroxy groups, BPs with hydroxy group at 4-position of benzene ring exhibited the highest toxicity, and lower toxicity was observed with hydroxy group at 3-position, while that at 2-position showed the lowest toxicity. Furthermore, 2- and 3- dimensional quantitative structure - activity relationships (2D- and 3D - QSAR) studies indicated that hydrogen-bond interactions and electrostatic effects were determinants for the acute toxicity of the BPs. The introduction of groups with positive charges at 4-position, or hydrogen-bond acceptor at 2-, 4-position of benzene ring, would lead to the increase of acute toxicity. The present results would provide useful information for evaluating the potential ecological and human health risk of BP-type UV-filters.

    Keywords: benzophenones; UV-filters; photobacterium assay; acute toxicity; quantitative structure - activity relationship

    紫外防曬劑作為添加劑已被廣泛用于個人護理用品、藥物、塑料、農用化學品等商業(yè)產品中,一項針對市售個人護理用品的調查中共檢出26種有機紫外防曬劑[1]。紫外防曬劑進入環(huán)境一般有2種途徑,其一是人們擦涂在身體上的防曬劑在游泳、洗浴等活動中直接進入環(huán)境水體,其二是生產和生活中使用的防曬劑經污水處理廠間接排入環(huán)境[2]。雖然每人每天防曬劑的使用量較小,但長期連續(xù)排放所帶來的環(huán)境風險不容忽視。二苯甲酮類化合物(BPs)是一類重要的紫外防曬劑,已有大量在環(huán)境介質、生物體中檢出的報道[3-7]。市政污水處理廠進水和二級處理出水中2-羥基-4-甲氧基-BP和2-羥基-4-甲氧基-5-磺酸基-BP的濃度范圍分別是1.5~19 μg·L-1和0.01~2.7 μg·L-1[8-12]。瑞士Glatt River河水中2-羥基-4-甲氧基-BP的累積濃度達到42 ng·L-1[6]。Wang等[13]研究了美國、韓國、日本、中國室內灰塵中2-羥基-4-甲氧基-BP、4-羥基-BP、2,4-二羥基-BP、2,2’-二羥基-4-甲氧基-BP和2,2’,4,4’-四羥基-BP的含量,5種防曬劑總含量的中位值范圍為78.3~612 ng·g-1。美國女性尿液中2-羥基-4-甲氧基-BP的含量為0.28~5 900 μg·L-1,2,4-二羥基-BP的含量為0.082~3 200 μg·L-1,4-羥基-BP的含量為0.082~22 μg·L-1[14]。

    盡管個人護理用品在上市前開展過相關的毒理學研究,但使用過程中也有部分報道顯示擦涂在皮膚上的二苯甲酮類化合物可能引起光致過敏[15]和細胞毒性[16],只是毒性效應并不十分顯著[17-19]。然而,隨著近年來毒理學研究的不斷深入,越來越多的體內、體外實驗表明BPs類化合物對內分泌系統(tǒng)有明顯的干擾效應[20-23]和發(fā)育毒性[24]。流行病學調查顯示男性暴露于2,2’,4,4’-四羥基-BP、4-羥基-BP等可導致夫妻生育能力分別下降31%和26%[25]。部分BPs類化合物已被列入“懷疑具有內分泌干擾效應的物質”[26]。除了內分泌干擾效應外,近年來也有研究者對該類化合物的遺傳毒性和急性毒性進行了評價,比如Zhao等[27]運用SOS/umu檢測了14種BPs類化合物的遺傳毒性,發(fā)現(xiàn)引起陽性效應的最低濃度范圍為81.97~>1 000 mg·L-1。SOS/umu和Ames實驗結果表明2,4-二羥基-BP在10 μg·mL-1暴露劑量下即顯示出活性[28]。2-羥基-4-甲氧基-BP和2-羥基-4-甲氧基-5-磺酸基-BP對大型蚤暴露的48 h-LC50分別為1.9和50 mg·L-1[12]。2-羥基-4-甲氧基-BP對鏈帶藻的72 h-IC10為0.56 mg·L-1[29]。Liu等[30]比較了系列BPs類化合物對大型蚤和其他水生生物的急性毒性,發(fā)現(xiàn)大型蚤作為受試生物更為敏感。

    細菌在生態(tài)系統(tǒng)中扮演非常重要的角色,因此通過研究污染物對細菌的毒性效應對評價污染物的生態(tài)風險具有重要意義?;诎l(fā)光細菌的急性毒性測試方法因具有簡便、快速、重現(xiàn)性好、與其他毒性測試方法所得結果相關性高等諸多優(yōu)點,被國內外研究者廣泛接受[31-32]。由于環(huán)境中化學污染物的種類繁多、毒性差異顯著,若對所有污染物一一進行各種毒性效應的檢測是不可能的。定量構效相關技術(QSAR)通過對部分化合物已有的毒性數(shù)據(jù)與其分子結構描述符進行統(tǒng)計建模的基礎上,實現(xiàn)對其他結構類似化合物毒性的科學預測,為環(huán)境中化學污染物的風險評價和管理提供了可能,也為揭示污染物的致毒機理提供參考[30,33-35]。

    因此,從BPs類防曬劑的分子結構出發(fā),選擇了14種具有代表性的化合物為研究對象,測定了它們對發(fā)光細菌的急性毒性效應。由于這些化合物都有共同的BP分子骨架,區(qū)別在于取代基的類型、位置和數(shù)量不同,因此,各化合物的毒性差異主要取決于母體骨架上取代基的不同。在定性分析各化合物毒性差異的基礎上,運用二維、三維定量構效相關方法對其毒性進行模擬預測,同時運用分子對接技術探討其毒性作用機制,為該類污染物的風險評價提供科學數(shù)據(jù)。

    1 材料與方法(Materials and methods)

    1.1 試劑與儀器

    14種BPs類紫外防曬劑(表1)、Na2HPO4·12H2O(純度99%)、NaH2PO4·2H2O(純度99%)、NaCl(純度99%)購自美國Sigma-Aldrich公司;胰蛋白胨、酵母浸膏購自BD公司;二甲亞砜(DMSO, 純度99.9%)購自Amresco公司。實驗所用儀器包括:Milli-Q型超純水儀(Millipore,美國)、Synergy2型多功能酶標儀(Bio-Tek,美國)、HPS-400型生化培養(yǎng)箱(哈爾濱東聯(lián))、LDZX-50FBS型高壓滅菌鍋(上海申安)、SW-CJ-1D型超凈工作臺(蘇州凈化)、MB100-2A型酶標板振蕩培養(yǎng)器(杭州奧盛)、移液器(Eppendorff,德國)。

    1.2 發(fā)光細菌急性毒性測試

    發(fā)光細菌毒性測試方法參照國家標準方法并有所改進[30]。實驗用菌株明亮發(fā)光桿菌T3小種(Photobacterium phosphorem T3)凍干粉由中國科學院南京土壤所提供,加入0.5 mL已滅菌的2%(質量分數(shù))NaCl溶液復蘇,取0.1 mL復蘇菌液加入20 mL濃度為3%(質量分數(shù))的NaCl溶液中作為實驗用菌液,21 ℃持續(xù)攪拌。將BPs母液用3%的NaCl溶液稀釋配制一系列不同濃度的工作液,進行毒性測試。取180 μL樣品溶液(或3% NaCl溶液)加入96孔板中,加入20 μL實驗菌液,在96孔板恒溫培養(yǎng)箱(21 ℃)中震蕩15 min,酶標儀測試各孔的發(fā)光強度。樣品組的發(fā)光強度記為Lsample,空白對照組的發(fā)光強度記為Lblank,每個樣品測定3次。發(fā)光抑制率IR (%) = (1- Lsample/Lblank) ×100%。分別以樣品濃度和相對發(fā)光抑制率為橫、縱坐標繪制半對數(shù)曲線,擬合求得相對抑制率為50%時所對應的濃度值,即為待測物的EC50值。SPSS軟件(Ver.12.0)用于毒性數(shù)據(jù)的統(tǒng)計分析。

    1.3 QSAR方法

    用ADRIANA.Code (Ver2.2.4)軟件計算受試化合物的分子結構參數(shù),共獲得29種分子描述符、8種形狀描述符和88種2D加權自相關拓撲描述符。以描述符為自變量、pEC50為因變量,用逐步多元回歸法建立定量構效相關方程。在此基礎上,運用Sybyl 7.0軟件的比較分子相似指數(shù)分析(CoMSIA)法進行三維構效關系研究。原子電荷用Gasteiger-Huckel方法計算,能量最小化使用鮑威爾共軛梯度算法。CoMSIA是由Sybyl中的QSAR可選操作實現(xiàn),在不作特殊聲明的情況下,其參數(shù)均用缺省值。在網格結點上以帶一個單位正電荷的sp3雜化碳原子為探針,每隔0.2 nm設一個測定點,用Tripos力場計算探針原子與化合物的相互作用能,空間場能和靜電場能閾值為30 kcal·mol-1。將14種BPs類化合物的急性毒性(pEC50)和CoMSIA值進行偏最小二乘(PLS)分析。以δ(2 kcal·mo1-1)作為列過濾閾值。交叉驗證得出最佳組分數(shù)和相關系數(shù),根據(jù)最佳組分數(shù)進行非交叉驗證,建立CoMSIA模型。

    1.4 分子對接Autodock方法

    用Autodock 4.0軟件將BPs分子與受體分子進行對接,按照幾何互補、能量互補、化學環(huán)境互補等原則實時評價配體與受體間的相互作用,識別2個分子間最佳結合模式。本研究采用半柔性對接方法,允許配體構象在一定范圍內變化。對接受體格點大小為40點×40點×40點,格點間距0.375 ?。計算中以半經驗勢函數(shù)為能量打分函數(shù),對小分子構象和位置進行全域搜索,其他參數(shù)為默認值。

    2 結果與討論(Results and discussion)

    14種BPs類防曬劑對發(fā)光菌的急性毒性測試結果列于表1,15 min EC50值范圍在17.67到243.82 mg·L-1之間。在這些化合物中,13# (2-羥基-4-甲氧基-5-磺酸基-BP)對發(fā)光菌的急性毒性最低,EC50值為243.82 mg·L-1,10# (2,2′,4,4′-四羥基-BP)的毒性最高,EC50值為17.67 mg·L-1。由于所有受試化合物都包含相同的骨架結構,其毒性差異主要取決于苯環(huán)上取代基位置和數(shù)量的變化。因此,可以通過研究取代基類型、位置和數(shù)量對毒性的影響以揭示其致毒機制。

    2.學習交互。目前有許多干部網絡學習平臺開設了學員論壇,但是交互模塊中學員之間的交流實效性較差,難以及時地解決問題,實現(xiàn)思維的碰撞,且個人的知識點多散落于論壇各版塊,難以形成知識的匯聚和共享。此外,對于學員發(fā)表內容的監(jiān)控也有一定難度。

    2.1 取代基對毒性的影響

    在所選擇的14種受試化合物中,大多數(shù)分子僅含有羥基官能團,因此選擇1#到10#化合物探討羥基數(shù)量對毒性的影響。從圖1a可見,BPs類化合物的急性毒性隨苯環(huán)上羥基數(shù)量的增多而升高。1#化合物無羥基,EC50值為67.09 mg·L-1,在10種BP類化合物中毒性最低;10#化合物有4個羥基,EC50值為17.67 mg·L-1,毒性最高。含有一個羥基的2#、3#和4#,2個羥基的5#、6#和7#,3個羥基的8#和10#都比不含羥基的1#毒性高。毒性從高到低遵循以下趨勢:4個羥基 > 3個羥基 > 2個羥基 > 1個羥基 > 沒有羥基。這一趨勢與其umu遺傳毒性規(guī)律一致[27,30]。

    表1 14種二苯甲酮類(BPs)紫外防曬劑結構及其對發(fā)光菌的急性毒性EC50值Table 1 The 14 tested benzophenone (BPs) compounds and their acute toxicity EC50 values for Photobacterium phosphorem T3

    除羥基數(shù)量之外,羥基在苯環(huán)上的位置對毒性也有影響。例如,2#、3#、4#化合物中僅有一個羥基,分別位于苯環(huán)的2-、3-、4-位,其毒性順序為:2- < 3- < 4-位羥基。與此類似,含有2個羥基的化合物2,2′-二羥基-BP (5#)和4,4′-二羥基-BP (6#),羥基在苯環(huán)4-位的6#化合物毒性大于羥基在苯環(huán)2-位的5#。這與SOS/umu遺傳毒性的規(guī)律相一致。此外,含有相同羥基數(shù)的化合物,羥基分布在2個苯環(huán)上的化合物比羥基分布在同一個苯環(huán)的毒性小。例如,2,3,4-三羥基-BP (9#)的毒性高于2,4,4′-三羥基-BP (8#);2,4-二羥基-BP (7#)比2,2′-二羥基-BP (5#)和4,4′-二羥基-BP (6#)有更高的毒性。這與該類化合物的SOS/umu遺傳毒性、內分泌干擾活性結果相似[21,27]。

    受試化合物分子中除羥基外,其他基團如甲氧基、辛氧基和磺酸基也對毒性有一定影響。比如5#和13#,苯環(huán)上4-位引入甲氧基降低了其毒性效應。再如7#、11#和12#,若羥基被甲氧基和辛氧基取代,毒性也明顯降低。若8#分子中的羥基被甲氧基取代即為14#,其毒性降低。此外,在11#分子結構的基礎上引入磺酸基即為13#,毒性亦顯著降低,可能是因為磺酸基的引入增加了化合物親水性,減少了其進入微生物細胞的量,導致毒性降低。若羥基被烷氧基取代,有可能減弱了原本BPs分子中的羥基與生物大分子之間的氫鍵作用,導致毒性降低。辛氧基的引入大大增加了分子體積,難以進入細胞產生毒性效應。這一結果與其對大型蚤、日本三角渦蟲急性毒性,以及SOS/umu遺傳毒性,雌激素受體活性等結果相似[21,27,30,36]。

    2.2 2D-QSAR研究

    用于多元線性回歸分析方法,建立了14種BP類化合物的發(fā)光細菌急性毒性pEC50與其分子結構描述符之間的定量構效相關模型:

    pEC50= 0.892CA - 7.281HLC - 53.328

    從方程式可以看出,BPs類紫外防曬劑的發(fā)光菌急性毒性與溶劑可接觸表面積呈正相關,與亨利常數(shù)呈負相關關系。這2個參數(shù)均與化合物分子的疏水性相關,溶劑可接觸分子表面積越大,化合物的疏水性越強;化合物的亨利常數(shù)越低,其疏水性越強,對發(fā)光細菌的毒性也越強。雖然2D-QSAR模型能夠較好地擬合14種化合物的急性毒性,但在模型驗證過程中發(fā)現(xiàn)其對其他結構類似化合物的毒性預測結果有較大出入。此外,通過統(tǒng)計分析方法自動入選的2個參數(shù)雖然有一定的理化意義,但仍難以揭示化合物的真正致毒機制。因此,有必要對BPs類化合物的毒性進行3D-QSAR研究。

    2.3 3D-QSAR研究

    運用CoMSIA分析方法,以8#化合物為模板進行分子疊合,經統(tǒng)計分析建立了BPs類紫外防曬劑對發(fā)光菌急性毒性(pEC50)與分子場信息的3D-QSAR模型。CoMSIA模型的統(tǒng)計參數(shù)(CoMSIA-1)如表2所示,立體場、靜電場、氫鍵受體場、氫鍵供體場和疏水場的貢獻比例分別為6.8%、23.9%、15.5%、43.8%和10.0%。由于立體場和疏水場在5個場中所占的比例最小,因此忽略立體場和疏水場后,基于靜電場、氫鍵供體場、氫鍵受體場重新構建了CoMSIA模型(CoMSIA-2)。靜電場、氫鍵受體場和氫鍵供體場的貢獻比例分別為31.0%、18.1%、50.9%。在3D-QSAR模型的3個場中,氫鍵供體場所占的比例是最大的,與CoMSIA-1模型相比,CoMSIA-2模型各參數(shù)都得到優(yōu)化,Q2由0.666增至0.737。

    圖1 取代基對BPs毒性的影響注:(a) 羥基數(shù)量對毒性的影響,(b) 烷氧基對毒性的影響。Fig. 1 The effects of substituents on acute toxicity of BPsNote: (a) effect of hydroxy numbers, (b) effects of alkoxy groups.

    圖2 CoMSIA模型色塊圖注:(a)靜電場,(b)氫鍵供體場,(c)氫鍵受體場。Fig. 2 Contour maps of CoMSIA modelsNote: (a) electrostatic field, (b) hydrogen-bond donor field, (c) hydrogen-bond acceptor field.

    圖3 BPs類化合物與熒光酶蛋白的分子對接模型Fig. 3 View of docked conformation including benzophenones and luciferase proteins

    CoMSIA-2中靜電場、氫鍵受體場和氫鍵供體場的色塊圖如圖2所示。通過分析色塊圖的顏色、大小及覆蓋區(qū)域,揭示取代基位置和性質對BPs類化合物發(fā)光菌急性毒性的影響。對于靜電場而言(圖2a),紅色塊和藍色塊的貢獻比例為80%和20%。紅色塊表示此處引入帶負電官能團毒性升高,藍色塊表示此處引入帶正電的官能團毒性升高。苯環(huán)4-位被藍色塊覆蓋,提示在苯環(huán)4-位上引入帶正荷的電官能團將增加發(fā)光菌急性毒性。此外,有紅色塊覆蓋在苯環(huán)附近,提示在苯環(huán)上引入供電子官能團可增加苯環(huán)上電子云密度,化合物的毒性也將增加。圖2b和2c顯示了氫鍵供體場和受體場的色塊圖。由于供體場和受體場所提示的信息是一致的,這里僅以氫鍵受體場(圖2c)為例進行解釋。紫紅色塊表示增加氫鍵受體毒性增強,紅色表示增加氫鍵受體毒性減弱。從圖6c看出,有紫紅色塊圍繞在苯環(huán)4-位,表示在該位置增加氫鍵受體將會增加發(fā)光菌毒性。比較化合物2#和10#、11#、12#,若在苯環(huán)的4-位引入含氧的羥基、甲氧基和辛氧基等氫鍵受體,其毒性顯著增強。由于所引入的3個氫鍵受體強弱順序為:羥基>甲氧基>辛氧基,該順序與化合物毒性降低順序一致。3塊紅色塊分布在苯環(huán)2-位、3-位和6-位,提示在這些位置引入氫鍵受體將降低發(fā)光菌急性毒性??赡苡捎谶@些位置上氫鍵供體的引入,增加了形成分子內氫鍵的可能性,降低了化合物與生物分子的相互作用,毒性降低。這與Liu等[30]針對BPs類化合物對大型蚤急性毒性的3D-QSAR結果相似。

    2.4 分子對接研究

    明亮發(fā)光桿菌T3小種(Photobacterium phosphorem T3)的發(fā)光機理是由還原態(tài)黃素單核苷酸(FMNH2)及長鏈脂肪醛在分子氧作用和胞內熒光酶催化下氧化為黃素單核苷酸(FMN)及長鏈脂肪酸,釋放出波長為450~490 nm藍綠光的過程。熒光素是導致生物體內催化或脂肪醛氧化發(fā)光的一類酶,是發(fā)光細菌發(fā)光的重要物質。熒光素酶是含有α、β兩個多肽亞基的單加氧酶,只有α、β兩個多肽亞基共存時才有活性。因此,本文采用計算機模擬熒光素酶的蛋白質結構,與BPs類化合物進行分子對接,以探討該類化合物的可能致毒機制。

    表3 BPs類化合物與熒光酶蛋白α-和β-亞基結合的吉普斯自由能(-dG)Table 3 The Gibbs free energy of benzophenones binding to α-subunit and β-subunit of luciferase proteins

    表2 CoMSIA模型參數(shù)Table 2 Summary of CoMSIA models’ parameters

    經Autodock 4.0軟件計算,分別得到BPs類化合物與熒光素蛋白酶α、β兩個多肽亞基的結合圖(如圖3)。由圖可知,無論是α還是β亞基,所有分子幾乎都結合在同一個位點上。對α亞基而言,與最多分子結合的蛋白氨基酸序列是TYR245,其次是ASP286;對β亞基而言,與最多分子結合的蛋白氨基酸序列是ASN9,其次是AVL15。事實上,在與BPs結合最多的2個氨基酸殘基色氨酸、天冬酰胺中含有羥基、酰胺基等,這些官能團容易與BPs化合物的官能團結合形成氫鍵,增強了化合物的毒性效應。在此基礎上,計算了BPs類化合物分別與熒光酶蛋白α和β亞基結合的吉普斯自由能(-dG)(列于表3),反映其與BPs類化合物的結合情況。-dG越小,表明蛋白質與化合物分子結合越緊密,化合物毒性越顯著。例如,毒性最低的13#化合物,其所對應的α和β亞基吉普斯自由能(-dG)分別為-6.06和-6.24,是所有化合物中-dG的最大值;毒性最強的10#化合物,其所對應的α和β亞基吉普斯自由能(-dG)分別是-7.25和-7.61,是所有化合物中-dG的最小值。

    綜上所述,通過對14種BPs類紫外防曬劑的發(fā)光菌急性毒性檢測和定量構效相關研究,發(fā)現(xiàn)該類化合物具有一定的急性毒性效應,15 min EC50值在幾十至幾百mg·L-1水平,長期排放可能導致潛在的生態(tài)風險。苯環(huán)上羥基取代基的數(shù)量越多,其毒性越高;羥基位于苯環(huán)4-位時毒性最高,位于苯環(huán)3-位時毒性次之,位于苯環(huán)2-位時毒性最低;若將羥基替換為烷氧基等取代基,毒性顯著降低。定量構效相關模型結果表明氫鍵供體場對毒性的貢獻最大,貢獻比例約占50%。在苯環(huán)4-位上引入帶正電荷的官能團將增加毒性,在苯環(huán)2-,4-位上引入氫鍵受體將顯著增加急性毒性。分子對接結果顯示,BPs類化合物容易分別與發(fā)光細菌熒光酶蛋白的α亞基和β亞基中的色氨酸、天冬酰胺殘基通過氫鍵結合,增強了化合物的毒性效應。本研究的結果為科學評價該類污染物的生態(tài)風險提供基礎數(shù)據(jù)。

    [1] Roelandts R, Vanhee J, Bonamie A, et al. A survey of ultraviolet absorbers in commercially available sun products [J]. International Journal of Dermatology, 1983, 22(4): 247-255

    [2] Tomson M B, Dauchy J, Hutchins S, et al. Groundwater contamination by trace level organics from a rapid infiltration site [J]. Water Research, 1981, 15(9): 1109-1116

    [3] Balmer M E, Buser H R, Muller M D, et al. Occurrence of some organic UV filters in wastewater, in surface waters, and in fish from Swiss lakes [J]. Environmental Science and Technology, 2005, 39(4): 953-962

    [4] Jeon H, Chung Y, Ryu J. Simultaneous determination of benzophenone type UV filters in water and soil by gas chromatography-mass spectrometry [J]. Journal of Chromatography A, 2006, 1131(1-2): 192-202

    [5] Loraine G A, Pettigrove M E. Seasonal variations in concentrations of pharmaceuticals and personal care products in drinking water and reclaimed wastewater in southern California [J]. Environmental Science and Technology, 2006, 40(3): 687-695

    [6] Negreira N, Rodriguez I, Rubi E, et al. Dispersive liquid-liquid microextraction followed by gas chromatography-mass spectrometry for the rapid and sensitive determination of UV filters in environmental water samples [J]. Analytical and Bioanalytical Chemistry, 2010, 398(2): 995-1004

    [7] Poiger T, Buser H R, Balmer M E. Occurrence of UV filter compounds from sunscreens in surface waters: Regional mass balance in two Swiss lakes [J]. Chemosphere, 2004, 55(7): 951-963

    [8] Vosough M, Mojdehi N R. Fast liquid chromatography-diode array detection assisted by chemometrics for quantification of seven ultraviolet filter in effluent wastewater [J]. Talanta, 2011, 85(4): 2175-2181

    [9] Li W H, Ma Y M, Guo C S. Occurrence and behavior of four of the most used sunscreen UV filters in a wastewater reclamation plant [J]. Water Research, 2007, 41(15): 3506-3512

    [10] Kupper T, Plagellat C, Braendli R C. Fate and removal of polycyclic musks, UV filters and biocides during wastewater treatment [J]. Water Research, 2006, 40(14): 2603-2612

    [11] Liu Y S, Ying G G, Shareef A, et al. Occurrence and removal of benzotriazoles and ultraviolet filters in a municipal wastewater treatment plant [J]. Environmental Pollution, 2012, 165: 225-232

    [12] Fent K, Kunz P Y, Zenker A. A tentative environmental risk assessment of the UV-filters 3-(4-methylbenzylidene-camphor), 2-ehyl-hexyl-4-trimethoxycinnamate, benzophenone-3, benzophenone-4 and 3-benzylidenne camphor [J]. Marine Environmental Research, 2010, 69(Suppl. 1): S4-S6

    [13] Wang L, Asimakopoulos A G, Moon H B, et al. Benzotriazole, benzothiazole and benzophenone compounds in indoor dust from the United States and east Asian countries [J]. Environmental Science and Technology, 2013, 47(9): 4752-4759

    [14] Kunisue T, Chen Z, Louis B G M, et al. Urinary concentrations of benzophenone-type UV filters in U.S. women and their association with endometriosis [J]. Environmental Science and Technology, 2012, 46(8): 4624-4632

    [15] Schauder S, Ippen H. Contact and photocontact sensitivity to sunscreens. Review of a 15 year experience and of the literature [J]. Contact Dermatiti, 1997, 37(5): 221-232

    [16] Matsumoto K, Akao Y, Kobayashi E,et al. Cytotoxic benzophenone derivatives from Garcinia species display a strong apoptosis-inducing effect against human leukemia cell lines [J]. Biological and Pharmaceutical Bulletin, 2003, 26(4): 569-571

    [17] Burdock G A, Pence D H, Ford R A. Safety evaluation of benzophenone [J]. Food and Chemical Toxicology, 1991, 29(11): 741-750

    [18] Dutta K, Das M, Rahman T. Toxicological impacts of benzophenone on the liver of guinea pig [J]. Bulletin of Environmental Contamination and Toxicology, 1993, 50(2): 282-285

    [19] Okereke C S, Barat S A, Abdel-Rahman M S. Safety evaluation of benzophenone-3 after dermal administration in rats [J]. Toxicology Letters, 1995, 80(1-3): 61-67

    [20] Schlumpf M, Schmid P, Durrer S, et al. Endocrine activity and developmental toxicity of cosmetic UV filters—An update [J]. Toxicology, 2004, 205(1-2): 113-122

    [21] Suzuki T, Kitamura S, Khota R, et al. Estrogenic and antiandrogenic activities of 17 benzophenone derivatives used as UV stabilizers and sunscreens [J]. Toxicology and Applied Pharmacology, 2005, 203(1): 9-17

    [22] Kunz P Y, Fent K. Estrogenic activity of UV filter mixtures [J]. Toxicology and Applied Pharmacology, 2006, 217(1): 86-99

    [23] Zucchi S, Bluthgen N, Ieronimo A.The UV-absorber benzophenone-4 alters transcripts of genes involved in hormonal pathways in zebrafish (Danio rerio) eleuthero-embryos and adult males [J]. Toxicology and Applied Pharmacology, 2011, 250(2): 137-146

    [24] Weisbrod C J, Kunz P Y, Zenker A K, et al. Effects of the UV filter benzophenone-2 on reproduction in fish [J]. Toxicology and Applied Pharmacology, 2007, 225(3): 255-266

    [25] Louis G M B, Kannan K, Sapra K J, et al. Urinary concentrations of benzophenone-type ultraviolet radiation filters and couples' fecundity [J]. American Journal of Epidemiology, 2014, 180(12): 1168-1175

    [26] Kerdivel G, Guevel R L, Habauzit D, et al. Estrogenic potency of benzophenone UV filters in breast cancer cells: Proliferative and transcriptional activity substantiated by docking analysis [J]. PLoS ONE, 2013, 8(4): e60567

    [27] Zhao H M, Wei D B, Li M, et al. Substituent contribution to the genotoxicity of benzophenone-type UV filters [J]. Ecotoxicology and Environmental Safety, 2013, 95: 241-246

    [28] Nakajima D, Asada S, Kageyama S, et al. Activity related to carcinogenicity of plastic additives in the benzophenone group [J]. Journal of University of Occupational and Environmental Health, 2006, 28(2): 143-156

    [29] Sieratowicz A, Kaiser D, Behr M, et al. Acute and chronic toxicity of four frequently used UV filter substances for Desmodesmus subspicatus and Daphnia magna [J]. Journal of Environmental Science and Health Part A, 2011, 46(12): 1311-1319

    [30] Liu H, Sun P, Liu H X, et al. Acute toxicity of benzophenone-type UV filters for Photobacterium phosphoreum and Daphnia magna: QSAR analysis, interspecies relationship and integrated assessment [J]. Chemosphere, 2015, 135: 182-188

    [31] Su L M, Zhang X J, Yuan X, et al. Evaluation of joint toxicity of nitroaromatic compounds and copper to Photobacterium phosphoreum and QSAR analysis [J]. Journal of Hazardous Materials, 2012, 241-242: 450-455

    [32] Zou X M, Lin Z F, Deng Z Q, et al. The joint effects of sulfonamides and their potentiator on Photobacterium phosphoreum: Differences between the acute and chronic mixture toxicity mechanisms [J]. Chemosphere, 2012, 86(1): 30-35

    [33] Li M, Wei D B, Zhao H M, et al. Genotoxicity of quinolones: Substituents contribution and transformation products QSAR evaluation using 2D and 3D models [J]. Chemosphere, 2014, 95: 220-226

    [34] Lyakurwa F S, Yang X H, Li X H, et al. Development of in silico models for predicting LSER molecular parameters and for acute toxicity prediction to fathead minnow (Pimephales promelas) [J]. Chemosphere, 2014, 108: 17-25

    [35] Yao Z F, Lin Z F, Wang T, et al. Using molecular docking-based binding energy to predict toxicity of binary mixture with different binding sites [J]. Chemosphere, 2013, 92(9): 1169-1176

    [36] Li M H. Acute toxicity of benzophenone-type UV filters and paraben preservatives to freshwater planarian, Dugesia japonica [J]. Toxicology and Environmental Chemistry, 2012, 94(3): 566-573

    AcuteToxicityandQSARStudiesonBenzophenone-typeUV-filterstoPhotobacterium

    Wei Dongbin1,2,*, Zhao Huimin1, Du Yuguo1,2

    1. State Key Laboratory of Environmental Chemistry and Ecotoxicology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China2. University of Chinese Academy of Sciences, Beijing 100049, China

    10.7524/AJE.1673-5897.20170123001

    2017-01-23錄用日期2017-03-13

    1673-5897(2017)3-234-09

    X171.5

    A

    魏東斌(1972-),男,環(huán)境科學博士,研究員,主要研究方向為有機污染化學,發(fā)表學術論文100余篇。

    國家自然科學基金(21377143, 21577154, 21590814)

    魏東斌(1972-),男,研究員,研究方向為有機污染化學,E-mail: weidb@rcees.ac.cn;

    魏東斌, 趙慧敏, 杜宇國. 二苯甲酮類紫外防曬劑發(fā)光菌急性毒性及QSAR研究[J]. 生態(tài)毒理學報,2017, 12(3): 234-242

    Wei D B, Zhao H M, Du Y G. Acute toxicity and QSAR studies on benzophenone-type UV-filters to photobacterium [J]. Asian Journal of Ecotoxicology, 2017, 12(3): 234-242 (in Chinese)

    猜你喜歡
    苯甲苯環(huán)甲氧基
    芳香族化合物同分異構體的書寫
    中學化學(2022年5期)2022-06-17 16:51:48
    HP-β-CD水相中4-甲氧基苯甲硫醚的選擇性氧化
    高效液相色譜法同時測定紡織品中11種二苯甲酮類紫外吸收劑
    世界上最苦的物質是什么?
    電腦報(2020年10期)2020-04-28 18:29:33
    對苯環(huán)上一元取代只有一種的同分異構體的系統(tǒng)分析
    2-(2-甲氧基苯氧基)-1-氯-乙烷的合成
    Z型三叉樹多肽聚苯環(huán)系統(tǒng)的Hosoya指標的計算公式
    DAD-HPLC法同時測定龍須藤總黃酮中5種多甲氧基黃酮
    中成藥(2017年4期)2017-05-17 06:09:50
    合成鄰甲氧基肉桂酸的兩步法新工藝
    具有極值hyper-Wiener指標的多聯(lián)苯鏈
    成人亚洲欧美一区二区av| 久久久国产成人免费| 69人妻影院| 成人漫画全彩无遮挡| 九九热线精品视视频播放| a级毛片免费高清观看在线播放| 2021少妇久久久久久久久久久| 内地一区二区视频在线| 日韩,欧美,国产一区二区三区 | 精品免费久久久久久久清纯| 午夜老司机福利剧场| 午夜精品国产一区二区电影 | 干丝袜人妻中文字幕| 精品久久久久久久人妻蜜臀av| 身体一侧抽搐| 日本猛色少妇xxxxx猛交久久| 综合色丁香网| 18禁动态无遮挡网站| 91狼人影院| 国产三级中文精品| 91aial.com中文字幕在线观看| 日韩强制内射视频| 一二三四中文在线观看免费高清| 美女内射精品一级片tv| 国产大屁股一区二区在线视频| www.av在线官网国产| 亚洲国产精品sss在线观看| 我要看日韩黄色一级片| 青春草国产在线视频| 国产精品野战在线观看| 又粗又爽又猛毛片免费看| 日本猛色少妇xxxxx猛交久久| 男人和女人高潮做爰伦理| 亚洲人与动物交配视频| 男女那种视频在线观看| 三级经典国产精品| 男女啪啪激烈高潮av片| 国产精品女同一区二区软件| 晚上一个人看的免费电影| 久久久久久久久久成人| 少妇被粗大猛烈的视频| av国产久精品久网站免费入址| 欧美高清性xxxxhd video| 岛国毛片在线播放| 婷婷六月久久综合丁香| 亚洲精品国产av成人精品| 高清毛片免费看| 国产白丝娇喘喷水9色精品| 国产一区二区在线观看日韩| 欧美xxxx黑人xx丫x性爽| 国产久久久一区二区三区| 久久久久久久久久久免费av| 美女国产视频在线观看| 在现免费观看毛片| 日韩一区二区视频免费看| 一区二区三区高清视频在线| 综合色av麻豆| 久久精品国产99精品国产亚洲性色| 尾随美女入室| 久久99蜜桃精品久久| 听说在线观看完整版免费高清| 国产中年淑女户外野战色| 午夜视频国产福利| 日本色播在线视频| 美女国产视频在线观看| 久久久国产成人免费| 一级毛片aaaaaa免费看小| 我的老师免费观看完整版| 亚洲av成人av| 韩国av在线不卡| 久久亚洲国产成人精品v| 亚洲精品,欧美精品| 中国国产av一级| 三级国产精品片| 国产又色又爽无遮挡免| 色噜噜av男人的天堂激情| 丝袜美腿在线中文| 国产精品麻豆人妻色哟哟久久 | 少妇的逼好多水| 久久久国产成人免费| 国产激情偷乱视频一区二区| 日韩欧美精品v在线| 天美传媒精品一区二区| 国产精品久久电影中文字幕| 白带黄色成豆腐渣| 久久99热这里只有精品18| videossex国产| 中文在线观看免费www的网站| 夜夜看夜夜爽夜夜摸| 国产精品人妻久久久影院| 1024手机看黄色片| 在线播放无遮挡| 深夜a级毛片| 69av精品久久久久久| 久久久久久久久久久丰满| 九九在线视频观看精品| 国产一区二区在线av高清观看| 中文字幕精品亚洲无线码一区| 免费av观看视频| 亚洲成人精品中文字幕电影| 最后的刺客免费高清国语| 一级爰片在线观看| 蜜臀久久99精品久久宅男| 国产精品女同一区二区软件| 精品人妻一区二区三区麻豆| 日韩av不卡免费在线播放| av在线天堂中文字幕| 欧美成人精品欧美一级黄| 成年版毛片免费区| 国产成人精品婷婷| 国产成人aa在线观看| 午夜免费男女啪啪视频观看| 欧美bdsm另类| 成人av在线播放网站| 国产又色又爽无遮挡免| 99久久精品热视频| 精品无人区乱码1区二区| 国产午夜精品论理片| av免费观看日本| 亚洲婷婷狠狠爱综合网| 秋霞伦理黄片| 色噜噜av男人的天堂激情| 日韩高清综合在线| 欧美高清成人免费视频www| 黄片无遮挡物在线观看| 欧美成人a在线观看| 国产亚洲av片在线观看秒播厂 | 少妇人妻一区二区三区视频| 亚洲精品国产av成人精品| 国产精品电影一区二区三区| 婷婷色综合大香蕉| 久久久久网色| 欧美激情久久久久久爽电影| 国产成人福利小说| 午夜精品国产一区二区电影 | 精品久久久噜噜| 亚洲精品自拍成人| 亚洲高清免费不卡视频| 国产伦精品一区二区三区四那| av播播在线观看一区| 99国产精品一区二区蜜桃av| 精品国内亚洲2022精品成人| 少妇被粗大猛烈的视频| 国产成人aa在线观看| 国产伦一二天堂av在线观看| 久久婷婷人人爽人人干人人爱| 国产精品熟女久久久久浪| 小蜜桃在线观看免费完整版高清| 成人国产麻豆网| 国产亚洲精品久久久com| 小蜜桃在线观看免费完整版高清| 日韩制服骚丝袜av| 精品人妻熟女av久视频| 成年版毛片免费区| 91久久精品国产一区二区三区| 亚洲精品456在线播放app| 小说图片视频综合网站| 国国产精品蜜臀av免费| av线在线观看网站| 亚洲精品日韩在线中文字幕| 欧美成人一区二区免费高清观看| 中文资源天堂在线| 搡女人真爽免费视频火全软件| 国产亚洲av嫩草精品影院| 天堂av国产一区二区熟女人妻| 97人妻精品一区二区三区麻豆| 午夜福利在线观看吧| 国产精品电影一区二区三区| 日本wwww免费看| 日本一二三区视频观看| 亚洲在久久综合| 色吧在线观看| 成年免费大片在线观看| 天堂中文最新版在线下载 | 日本与韩国留学比较| 久久精品久久精品一区二区三区| 秋霞伦理黄片| 精品一区二区免费观看| 麻豆精品久久久久久蜜桃| 九草在线视频观看| 有码 亚洲区| 18禁动态无遮挡网站| 国产午夜福利久久久久久| 久久久精品欧美日韩精品| 日本三级黄在线观看| 日韩欧美 国产精品| 青青草视频在线视频观看| 亚洲天堂国产精品一区在线| 中文欧美无线码| av国产久精品久网站免费入址| 亚洲av中文字字幕乱码综合| 久久99热这里只有精品18| 国产麻豆成人av免费视频| www.av在线官网国产| 亚洲av成人av| 亚洲精品一区蜜桃| 国产真实伦视频高清在线观看| 丰满少妇做爰视频| 国产免费男女视频| 亚洲欧洲日产国产| 老师上课跳d突然被开到最大视频| 欧美又色又爽又黄视频| 老司机影院成人| 久久久久久久久久成人| 国产一区二区亚洲精品在线观看| 国产亚洲av片在线观看秒播厂 | 99热网站在线观看| 晚上一个人看的免费电影| 我的老师免费观看完整版| 欧美性感艳星| 天堂网av新在线| 亚洲自拍偷在线| 国产在视频线精品| 国产精品一区二区三区四区免费观看| 噜噜噜噜噜久久久久久91| 欧美高清性xxxxhd video| 国产爱豆传媒在线观看| 美女内射精品一级片tv| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人精品婷婷| 精品久久久久久久久久久久久| 久久热精品热| 一本久久精品| 国产精品.久久久| 亚洲欧洲日产国产| 久久人妻av系列| 久久99精品国语久久久| 日韩亚洲欧美综合| 亚洲国产精品合色在线| 乱人视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 成人毛片60女人毛片免费| 亚洲欧美日韩卡通动漫| 亚洲,欧美,日韩| 又粗又硬又长又爽又黄的视频| 亚洲精品456在线播放app| 边亲边吃奶的免费视频| 日本一二三区视频观看| 99久久九九国产精品国产免费| 97在线视频观看| 久久精品国产鲁丝片午夜精品| 国产av在哪里看| 毛片一级片免费看久久久久| 欧美97在线视频| 亚洲第一区二区三区不卡| 99热网站在线观看| 亚洲国产精品sss在线观看| 特级一级黄色大片| 亚洲av二区三区四区| 国内精品美女久久久久久| 美女大奶头视频| 一边亲一边摸免费视频| 一区二区三区免费毛片| 免费看av在线观看网站| 秋霞伦理黄片| 国内精品美女久久久久久| 亚洲av.av天堂| 少妇猛男粗大的猛烈进出视频 | 欧美又色又爽又黄视频| 国产亚洲av嫩草精品影院| 赤兔流量卡办理| 菩萨蛮人人尽说江南好唐韦庄 | 寂寞人妻少妇视频99o| 日本爱情动作片www.在线观看| 免费观看精品视频网站| 我的女老师完整版在线观看| 99热这里只有精品一区| 91aial.com中文字幕在线观看| 男人和女人高潮做爰伦理| 亚洲欧美日韩东京热| 九色成人免费人妻av| 99热这里只有精品一区| 免费观看a级毛片全部| 国产人妻一区二区三区在| 亚洲成色77777| 丰满乱子伦码专区| 大又大粗又爽又黄少妇毛片口| 久久精品夜夜夜夜夜久久蜜豆| 大又大粗又爽又黄少妇毛片口| eeuss影院久久| 三级毛片av免费| 99视频精品全部免费 在线| 偷拍熟女少妇极品色| 成人美女网站在线观看视频| 精品久久久久久电影网 | 久久久成人免费电影| 建设人人有责人人尽责人人享有的 | 亚洲性久久影院| 亚洲国产欧美在线一区| 久久人人爽人人片av| 亚洲天堂国产精品一区在线| 韩国av在线不卡| 国产日韩欧美在线精品| 国产毛片a区久久久久| 日本免费一区二区三区高清不卡| 美女黄网站色视频| 久久久久网色| 看片在线看免费视频| 在线天堂最新版资源| 日韩中字成人| 久久综合国产亚洲精品| 国产精品一区www在线观看| 亚洲av不卡在线观看| 欧美人与善性xxx| 亚洲乱码一区二区免费版| 男女边吃奶边做爰视频| 中文乱码字字幕精品一区二区三区 | 国产视频首页在线观看| 狂野欧美白嫩少妇大欣赏| 精品久久久久久久人妻蜜臀av| 在线观看av片永久免费下载| 麻豆国产97在线/欧美| 成人av在线播放网站| 国产在线男女| 国产一区二区在线av高清观看| 国产在视频线在精品| 精品久久久久久久久亚洲| 人妻系列 视频| 国产女主播在线喷水免费视频网站 | 欧美色视频一区免费| 欧美最新免费一区二区三区| 99久久无色码亚洲精品果冻| 看十八女毛片水多多多| 欧美激情在线99| 国产伦在线观看视频一区| 床上黄色一级片| 最近最新中文字幕大全电影3| 欧美bdsm另类| 亚洲av电影不卡..在线观看| 国产一区二区在线观看日韩| 欧美极品一区二区三区四区| www日本黄色视频网| 国产亚洲一区二区精品| 蜜臀久久99精品久久宅男| 日韩 亚洲 欧美在线| 成年女人看的毛片在线观看| 男插女下体视频免费在线播放| h日本视频在线播放| 69人妻影院| 国产午夜精品论理片| 色综合站精品国产| 久久久久久久久中文| 亚洲国产欧美在线一区| 免费观看精品视频网站| 精品久久久久久久久av| 美女被艹到高潮喷水动态| 五月玫瑰六月丁香| 国产免费又黄又爽又色| 欧美精品国产亚洲| 成人欧美大片| 久久久久久久午夜电影| 波多野结衣巨乳人妻| 精品久久久久久久久av| 欧美变态另类bdsm刘玥| 激情 狠狠 欧美| 赤兔流量卡办理| 久久久精品94久久精品| 国产激情偷乱视频一区二区| 国产亚洲午夜精品一区二区久久 | 99久久成人亚洲精品观看| 国内揄拍国产精品人妻在线| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 观看免费一级毛片| 精品午夜福利在线看| 国产 一区 欧美 日韩| 九九久久精品国产亚洲av麻豆| 99热这里只有精品一区| 亚洲国产精品专区欧美| 老师上课跳d突然被开到最大视频| 草草在线视频免费看| 亚洲国产成人一精品久久久| 成人国产麻豆网| 亚洲美女搞黄在线观看| 99久久中文字幕三级久久日本| 国产乱人偷精品视频| 日本wwww免费看| 免费黄网站久久成人精品| 看非洲黑人一级黄片| 18禁动态无遮挡网站| 亚洲国产精品sss在线观看| 亚洲精品456在线播放app| 国产免费福利视频在线观看| 中文精品一卡2卡3卡4更新| 国产精品嫩草影院av在线观看| av黄色大香蕉| 日本免费在线观看一区| 美女国产视频在线观看| 性插视频无遮挡在线免费观看| 亚洲国产精品国产精品| 国产免费一级a男人的天堂| 我要看日韩黄色一级片| 国产一级毛片在线| 国产老妇女一区| 亚洲精品国产成人久久av| 亚洲av不卡在线观看| 国产午夜福利久久久久久| 舔av片在线| 亚洲精品一区蜜桃| 天天躁夜夜躁狠狠久久av| 免费观看人在逋| av在线亚洲专区| 乱人视频在线观看| 欧美激情国产日韩精品一区| 国产免费视频播放在线视频 | 网址你懂的国产日韩在线| 男女视频在线观看网站免费| 色5月婷婷丁香| 亚洲国产欧洲综合997久久,| 免费人成在线观看视频色| 国产黄a三级三级三级人| 亚洲丝袜综合中文字幕| 最近视频中文字幕2019在线8| 亚洲欧美成人精品一区二区| 日本爱情动作片www.在线观看| 岛国毛片在线播放| 一本久久精品| 久久久久久九九精品二区国产| 国产麻豆成人av免费视频| 国产亚洲5aaaaa淫片| 日本黄色片子视频| 久久精品久久久久久噜噜老黄 | 国产黄a三级三级三级人| 国产在线男女| 69人妻影院| 欧美最新免费一区二区三区| 国产老妇伦熟女老妇高清| 国产爱豆传媒在线观看| 成人二区视频| 精品一区二区免费观看| 国产成人精品一,二区| 麻豆国产97在线/欧美| 中文资源天堂在线| 欧美zozozo另类| 一级黄色大片毛片| 午夜爱爱视频在线播放| 亚洲乱码一区二区免费版| 一卡2卡三卡四卡精品乱码亚洲| 国产伦精品一区二区三区视频9| 级片在线观看| 国产成人精品婷婷| 亚洲欧美成人精品一区二区| 99久久无色码亚洲精品果冻| 寂寞人妻少妇视频99o| 三级毛片av免费| 深夜a级毛片| 精品99又大又爽又粗少妇毛片| 我的女老师完整版在线观看| 免费人成在线观看视频色| 亚洲18禁久久av| 天堂av国产一区二区熟女人妻| 亚洲最大成人av| 国产成人freesex在线| 国产精品久久久久久av不卡| 国产午夜福利久久久久久| 国产麻豆成人av免费视频| 少妇人妻一区二区三区视频| 中文字幕人妻熟人妻熟丝袜美| 中文精品一卡2卡3卡4更新| 18+在线观看网站| 搡老妇女老女人老熟妇| 国产av码专区亚洲av| 亚洲成人久久爱视频| 亚洲真实伦在线观看| 99热这里只有精品一区| 日日撸夜夜添| 亚洲av中文字字幕乱码综合| av视频在线观看入口| 99视频精品全部免费 在线| 国产精品一区www在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲av不卡在线观看| 免费不卡的大黄色大毛片视频在线观看 | 中国美白少妇内射xxxbb| 岛国在线免费视频观看| 亚洲中文字幕日韩| 麻豆精品久久久久久蜜桃| eeuss影院久久| 日本一本二区三区精品| 欧美性感艳星| 小说图片视频综合网站| 精品人妻熟女av久视频| 五月伊人婷婷丁香| 国产成人a区在线观看| 丰满乱子伦码专区| 三级经典国产精品| 简卡轻食公司| 久久精品国产自在天天线| 又黄又爽又刺激的免费视频.| 久久久久久久午夜电影| 伦理电影大哥的女人| 国内精品宾馆在线| 久久99蜜桃精品久久| 99久久人妻综合| 变态另类丝袜制服| 亚洲国产精品合色在线| 国产亚洲精品久久久com| 国产亚洲5aaaaa淫片| 99久久人妻综合| 男女边吃奶边做爰视频| 伦理电影大哥的女人| 99在线视频只有这里精品首页| 亚洲欧美日韩无卡精品| 国产又黄又爽又无遮挡在线| 青春草亚洲视频在线观看| 国产精品.久久久| 好男人在线观看高清免费视频| 日韩欧美精品v在线| 午夜久久久久精精品| 九九热线精品视视频播放| 女人久久www免费人成看片 | 亚洲成人中文字幕在线播放| 亚洲人成网站高清观看| 美女xxoo啪啪120秒动态图| 亚洲经典国产精华液单| 午夜爱爱视频在线播放| 久久久久精品久久久久真实原创| 日韩视频在线欧美| av视频在线观看入口| 三级国产精品欧美在线观看| 久久久国产成人精品二区| 色吧在线观看| 国产伦理片在线播放av一区| 18禁裸乳无遮挡免费网站照片| 欧美又色又爽又黄视频| 日韩欧美国产在线观看| 欧美性感艳星| 天堂中文最新版在线下载 | 一二三四中文在线观看免费高清| 国产高潮美女av| 亚洲精品456在线播放app| 精品午夜福利在线看| 亚洲在久久综合| 一本一本综合久久| www日本黄色视频网| 免费观看精品视频网站| 丝袜美腿在线中文| 日本免费一区二区三区高清不卡| 亚洲美女搞黄在线观看| 亚洲精品国产成人久久av| 老司机影院毛片| 久久久久性生活片| 国产成人免费观看mmmm| 亚洲无线观看免费| 禁无遮挡网站| 久久国产乱子免费精品| 国产免费又黄又爽又色| 岛国在线免费视频观看| 久久热精品热| 亚洲国产精品国产精品| 18禁在线播放成人免费| 国产一级毛片七仙女欲春2| 国产v大片淫在线免费观看| 国产av不卡久久| 岛国毛片在线播放| 51国产日韩欧美| 99久久精品热视频| 日韩在线高清观看一区二区三区| 亚洲欧美清纯卡通| 深夜a级毛片| 成人一区二区视频在线观看| 日日摸夜夜添夜夜爱| 国内少妇人妻偷人精品xxx网站| 久久精品国产亚洲av天美| 91久久精品国产一区二区三区| .国产精品久久| 国产精品一二三区在线看| 亚州av有码| 91精品一卡2卡3卡4卡| 午夜免费激情av| 日本色播在线视频| 免费一级毛片在线播放高清视频| 欧美成人一区二区免费高清观看| 亚洲乱码一区二区免费版| 夜夜爽夜夜爽视频| 色播亚洲综合网| 三级国产精品欧美在线观看| 国产中年淑女户外野战色| 哪个播放器可以免费观看大片| 久久精品国产亚洲av天美| 国产成人91sexporn| 中文在线观看免费www的网站| 国产精品人妻久久久影院| 国产精品野战在线观看| 99九九线精品视频在线观看视频| 亚洲,欧美,日韩| 少妇猛男粗大的猛烈进出视频 | 亚洲人成网站在线播| 我的女老师完整版在线观看| 精品国内亚洲2022精品成人| 中国国产av一级| 久久鲁丝午夜福利片| 欧美成人一区二区免费高清观看| 我的老师免费观看完整版| 久久久久性生活片| 午夜a级毛片| 国产色爽女视频免费观看| 男女国产视频网站| 黄色配什么色好看| 免费观看的影片在线观看| 91午夜精品亚洲一区二区三区| 嘟嘟电影网在线观看| 色综合亚洲欧美另类图片| 可以在线观看毛片的网站| 午夜亚洲福利在线播放| 成人亚洲欧美一区二区av| 亚洲精品456在线播放app| 一级毛片电影观看 | 美女内射精品一级片tv| 亚洲精品亚洲一区二区| 久久人妻av系列| 欧美zozozo另类| 国产三级中文精品| 精品久久久久久成人av| 久久人人爽人人爽人人片va|