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

    溴氰菊酯農(nóng)藥殘留檢測的分子印跡預(yù)聚體系篩選及吸附性能

    2019-01-18 00:44:16周劉梅解新安
    農(nóng)業(yè)工程學(xué)報 2019年1期
    關(guān)鍵詞:溴氰菊酯印跡菊酯

    李 璐,周劉梅,解新安,李 雁

    ?

    溴氰菊酯農(nóng)藥殘留檢測的分子印跡預(yù)聚體系篩選及吸附性能

    李 璐,周劉梅,解新安,李 雁

    (華南農(nóng)業(yè)大學(xué)食品學(xué)院,廣州 510640)

    為快速篩選及確定合適的預(yù)聚體系,以高效制備對茶葉等食品中殘留的溴氰菊酯具有特異識別能力的分子印跡聚合物,以溴氰菊酯為模板,采用Materials Studio軟件和紫外光譜對預(yù)聚體系中功能單體和溶劑進行種類篩選及配比確定,輔助溴氰菊酯分子印跡聚合物的合成設(shè)計。對模板與不同功能單體預(yù)聚體系的結(jié)合能計算和光譜分析結(jié)果均表明溴氰菊酯與丙烯酰胺能形成更為穩(wěn)定的復(fù)合物,且以乙腈為溶劑更有利于溴氰菊酯分子印跡聚合物的制備。差示紫外光譜計算進一步得出當溴氰菊酯與丙烯酰胺的化學(xué)配位比=4時,可在乙腈溶液中形成穩(wěn)定的預(yù)聚體系,其結(jié)合常數(shù)為6.614×105。基于上述結(jié)果,采用沉淀聚合法正交優(yōu)化制備溴氰菊酯分子印跡聚合物,掃描電鏡和靜態(tài)平衡吸附試驗結(jié)果表明所制備聚合物具有良好的多孔結(jié)構(gòu)和特異吸附性,印跡因子為2.269,對3種結(jié)構(gòu)類似物的分離因子均大于3.2。研究結(jié)果表明分子模擬和光譜分析方法對分子印跡預(yù)聚體系篩選和聚合物合成設(shè)計具有重要的指導(dǎo)意義。

    農(nóng)藥;分子動力學(xué);紫外可見光譜

    0 引 言

    擬除蟲菊酯類農(nóng)藥憑借殺蟲譜廣、高效等優(yōu)點被廣泛應(yīng)用于農(nóng)業(yè)生產(chǎn)和公共衛(wèi)生領(lǐng)域[1]。但該類農(nóng)藥經(jīng)食物鏈富集進入人體內(nèi)后會引起頭暈、惡心和皮膚刺激等不良反應(yīng),而且其對機體的免疫力和心血管也有明顯的毒害作用[2-3]。因此,建立一種快速、靈敏、便捷的菊酯類農(nóng)藥殘留檢測方法對于保障食品安全很有必要。分子印跡聚合物(molecularly imprinted polymers,MIPs)作為一類性質(zhì)穩(wěn)定、成本低且對模板分子具有特異識別性能的新型材料現(xiàn)已被廣泛應(yīng)用于農(nóng)殘檢測的色譜分離、固相萃取和傳感器等領(lǐng)域[4-6]。分子印跡聚合中功能單體與模板分子形成的預(yù)聚體系及聚合溶劑環(huán)境對所制備MIPs的吸附和識別性能有舉足輕重的影響[7-8]。目前常用的功能單體和溶劑的篩選方法是試錯法,即結(jié)合參考文獻設(shè)計一系列的試驗,通過對比試驗結(jié)果,得到最優(yōu)的MIPs制備方案。該方法雖能夠優(yōu)化合成方案,但是盲目性強且試驗周期長、試錯成本高。因此,通過研究預(yù)聚體系中模板分子-功能單體間的相互作用情況以實現(xiàn)印跡條件的快速優(yōu)化成為分子印跡技術(shù)的研究熱點。

    本文以擬除蟲菊酯類農(nóng)藥中毒性強、應(yīng)用范圍廣的溴氰菊酯(decamethrin,DM)為模板分子,采用計算機分子動力學(xué)和紫外光譜法綜合研究溴氰菊酯與不同功能單體間的相互作用強弱,完成對功能單體種類及其與模板分子適宜比例的篩選,同時確定適宜的聚合溶劑。以前述預(yù)聚體的研究結(jié)果為基礎(chǔ),設(shè)計正交試驗優(yōu)化沉淀聚合法制備溴氰菊酯MIPs,并對其微觀結(jié)構(gòu)和吸附特性進行研究。研究結(jié)果將為分子印跡技術(shù)應(yīng)用于擬除蟲菊酯類農(nóng)藥的殘留檢測提供理論依據(jù)和數(shù)據(jù)支持。

    1 材料與方法

    1.1 基礎(chǔ)試劑與儀器

    溴氰菊酯、氯氰菊酯、高效氯氟氰菊酯和聯(lián)苯菊酯(質(zhì)量分數(shù)>98%),江蘇榮德試劑廠;丙烯酸(acrylic acid,AA)、丙烯酰胺(acrylamide,AM)、甲基丙烯酸(methacrylic acid,MAA)、乙二醇二甲基丙烯酸酯、偶氮二異丁腈,分析純,阿拉??;甲醇、溴化鉀,分析純,浙江晉巨化工有限公司;冰乙酸、乙腈,分析純,天津大茂化學(xué)試劑廠。

    UV-2500紫外-可見分光光度計,德國斯派克分析儀器公司;S-4800掃描電子顯微鏡,日本Hitachi公司;Gemini VII 2390全自動比表面及孔隙分析儀,美國麥克公司;Vertex 70傅里葉變換紅外光譜儀,德國布魯克公司。

    1.2 分子印跡預(yù)聚體系篩選

    1.2.1 分子動力學(xué)模擬

    計算模擬硬件采用Intel Pentium Dual-Core 2.8 GHz E6300處理器,軟件采用Materials Studio軟件的Amorphous Cell模塊和Discover模塊。

    模擬過程采用Compass力場,首先通過Amorphous Cell模塊建立包含1個模板分子和不同功能單體比例的預(yù)聚體系,然后用Discover模塊中的minimizer尋找到能量最低構(gòu)象后通過Discover進行動力學(xué)模擬。動力學(xué)模擬采用NVT系綜,模擬初始速度由Boltzmann分布隨機生成,初始溫度為298 K,模擬過程中采用Andersen恒溫浴以保持溫度不變。系統(tǒng)達到平衡后統(tǒng)計復(fù)合物、模板分子和功能單體的單點能,并按式(1)計算模板與功能單體相互作用結(jié)合能。

    式中為各組分的單點能,kJ/mol。

    分子動力學(xué)模擬中采用顯性溶劑模型研究溶劑化效應(yīng),模型采用Amorphous Cell模塊建立包含1個溴氰菊酯、經(jīng)確定的功能單體用量及10個溶劑分子的預(yù)聚體系,其余參數(shù)同上。對溶劑模型,結(jié)合能計算公式如式(2)。

    1.2.2 溴氰菊酯、功能單體及其混合物紫外吸收光譜掃描

    取一系列具塞試管,分別添加10L的1 mmol/L溴氰菊酯乙腈溶液、20L的10 mmol/L功能單體乙腈溶液和摩爾比為1∶2的溴氰菊酯-功能單體混合溶液,加入乙腈定容至10 mL,超聲(600 W,20 min),室溫靜置14 h。乙腈溶液為空白參比,在紫外光譜區(qū)(190~400 nm)對溴氰菊酯、功能單體及其混合物進行紫外吸收光譜掃描。

    1.2.3 不同摩爾比的模板分子-功能單體混合物預(yù)組裝體系紫外光譜掃描

    取一系列具塞試管,每根試管內(nèi)均加入1 mmol/L溴氰菊酯乙腈溶液10L,使得試管中溴氰菊酯的濃度固定為0.000 1 mmol/L;依次加入不同體積的10 mmol/L功能單體乙腈溶液,使得溴氰菊酯與功能單體的摩爾比為1∶2、1∶4、1∶6、1∶8、1∶10、1∶12、1∶14,加入乙腈定容至10 mL,超聲(600 W,20 min),室溫靜置14 h。以0.000 1 mmol/L溴氰菊酯乙腈溶液為空白參比,在紫外光譜區(qū)(190~400 nm)對模板分子-功能單體混合物預(yù)組裝體系進行紫外吸收光譜掃描。

    1.3 材料制備與表征

    1.3.1 溴氰菊酯分子印跡聚合物和非印跡聚合物的制備

    采用沉淀聚合法制備溴氰菊酯分子印跡聚合物。稱取1.010 g (2 mmol)溴氰菊酯、0.569 g (8 mmol)丙烯酰胺于錐形瓶中,再加入80 mL乙腈溶解,置于磁力攪拌器上常溫攪拌1 h,使模板分子與功能單體充分作用形成超分子配合物。然后加入7.920 g (40 mmol)交聯(lián)劑乙二醇二甲基丙烯酸酯和1.314 g (8 mmol)引發(fā)劑偶氮二異丁腈,繼續(xù)磁力攪拌10 min,使超分子配合物與交聯(lián)劑和引發(fā)劑充分接觸,混合均勻,形成透明體系。超聲脫氣10 min后將混合物轉(zhuǎn)移至三口燒瓶中,在60 ℃恒溫水浴鍋中機械攪拌反應(yīng)24 h,攪拌速度控制在100~200 r/min。

    對所得聚合物進行抽濾,并用無水乙醇和蒸餾水各清洗2次,于65 ℃烘箱中干燥12 h;以甲醇乙酸溶液(體積比1∶1)作為洗脫溶劑,85 ℃恒溫抽提10 h以除去模板分子,用甲醇洗去乙酸后再用蒸餾水反復(fù)清洗聚合物,最后將其置于65 ℃烘箱中烘干至質(zhì)量恒定,得到的白色粉末狀聚合物即為溴氰菊酯分子印跡聚合物。

    非印跡聚合物(non-imprinted polymers,NIPs)的制備及處理方法同MIPs,但不加入模板分子。

    1.3.2 溴氰菊酯分子印跡聚合物和非印跡聚合物的表征

    將溴氰菊酯MIPs和NIPs進行噴金處理后,用掃描電子顯微鏡觀測形貌。使用Gemini VII 2390分析儀測量MIPs和NIPs比表面積和孔隙參數(shù)。采用溴化鉀壓片法,用Vertex 70傅里葉變換紅外光譜儀對MIPs和NIPs進行紅外光譜分析(波數(shù)范圍:4 000~400 cm–1)。

    1.4 溴氰菊酯分子印跡聚合物吸附性能試驗設(shè)計

    1.4.1 靜態(tài)吸附試驗

    稱取25 mg聚合物,置于裝有25 mL濃度為100 mg/L 的溴氰菊酯乙腈溶液中,在室溫下水浴振蕩24 h (100 r/min)。靜置5 min后,過濾,取濾液稀釋后用紫外-可見光分光光度計在220 nm處測定吸光度,根據(jù)標準曲線(Abs=0.067 81+0.079 9,2=0.998 3,標準曲線適用濃度范圍為2~12 mg/L)計算平衡吸附后的物質(zhì)濃度,并由式(3)計算聚合物的吸附量。

    式中1為平衡吸附后溶液中目標物質(zhì)的濃度,mg/L;0為目標物質(zhì)的起始濃度,mg/L;為聚合物的質(zhì)量,g;e為靜態(tài)平衡吸附量,mg/g;為物質(zhì)溶液的體積,mL。

    1.4.2 選擇性吸附試驗

    稱取25 mg聚合物加入濃度為100 mg/L的溴氰菊酯溶液及其結(jié)構(gòu)類似物(即氯氰菊酯、高效氯氟氰菊酯和聯(lián)苯菊酯)溶液25 mL,在30 ℃下振蕩24 h后取濾液,于各菊酯類農(nóng)藥的特征波長下用紫外-可見光分光光度計測定吸光度,根據(jù)標準曲線計算平衡吸附后的物質(zhì)濃度,并根據(jù)式(3)、(4)、(5)和(6)計算溴氰菊酯MIPs對各菊酯類農(nóng)藥的吸附量(e)、靜態(tài)分配系數(shù)()、分離因子()和印跡因子(IF)。

    式中P為底物結(jié)合在聚合物上的濃度,mg/L;IF用于表征MIPs及NIPs的特異性結(jié)合能力,IF越大表明特異性越強;和分別為模板分子和競爭分子。為靜態(tài)分配系數(shù),越大表明聚合物的吸附能力越大;MIPs和NIPs分別為MIPs和NIPs的靜態(tài)分配系數(shù);>1說明溴氰菊酯MIPs能識別模板分子,且越大說明印跡效果越好。

    2 結(jié)果與討論

    2.1 計算機分子動力學(xué)模擬

    功能單體的選擇主要由模板分子的結(jié)構(gòu)特征確定,它必須能與模板分子作用且足夠強,以利于MIPs中選擇性識別位點的形成。溴氰菊酯分子中含有鹵素雙鍵、酯鍵、氰基等官能團,有3個不對稱的碳原子,在酸性介質(zhì)中較穩(wěn)定。因此初選甲基丙烯酸(methacrylic acid, MAA)和丙烯酸(acrylic acid, AA)2種酸性單體、丙烯酰胺(acrylamide, AM)中性單體為功能單體,設(shè)定模板分子與功能單體的比例為1∶2、1∶4、1∶6、1∶8,用1.2.1中所述模擬方法,在Compass力場和NVT系綜下對模板分子與功能單體構(gòu)型進行優(yōu)化,得到復(fù)合物平衡構(gòu)象,然后采用結(jié)合能Δ來表征溴氰菊酯和不同單體之間的相互作用,結(jié)果見圖1。

    圖1 溴氰菊酯與不同功能單體在不同比例下的結(jié)合能(DE)

    模板分子溴氰菊酯與功能單體形成復(fù)合物時結(jié)合能越低,說明模板與功能單體的結(jié)合作用力越大,復(fù)合物也相對越穩(wěn)定[18]。由圖1可知,3種復(fù)合物的結(jié)合能均隨著功能單體比例的增加而降低,當模板分子與功能單體的比例在1∶6~1∶8范圍時結(jié)合能的變化趨于平緩,其中DM-AM復(fù)合物在1∶2~1∶4的比例范圍內(nèi)結(jié)合能的下降趨勢最為明顯且結(jié)合能最低,這說明3種功能單體均能與溴氰菊酯相互作用,而DM與AM之間可形成更為穩(wěn)定的預(yù)聚體系。據(jù)此推測制備DM MIPs時AM為更合適的功能單體,DM和AM的摩爾比在1∶4~1∶6所得聚合物可能更穩(wěn)定。

    基于以上結(jié)果,采用Amorphous Cell模塊建立包含1個DM、4個AM及10個溶劑(溶劑種類為乙腈、丙酮、氯仿和正己烷)的預(yù)聚體系,對不同溶劑中溴氰菊酯和AM的結(jié)合能進行計算,結(jié)果如表1所示。

    表1 不同溶劑中DM-4AM復(fù)合物的結(jié)合能

    由表1可知,1個DM和4個AM復(fù)合物體系中,結(jié)合能的大小為正己烷>氯仿>丙酮>乙腈,說明復(fù)合物與乙腈的相互作用最小,有利于提供較好的聚合環(huán)境。用此配比在不同溶劑中按1.3.1所述方法進行沉淀聚合,發(fā)現(xiàn)以氯仿為聚合溶劑時,無沉淀析出,合成失??;以丙酮為聚合溶劑時,所制備的MIPs粒徑較大,且按1.4.1所述方式測定的靜態(tài)吸附量(23.81 mg/g)較低,而在乙腈溶劑中合成的溴氰菊酯MIPs具有較好的印跡效果,其靜態(tài)吸附量為51.04 mg/g,試驗結(jié)果表明可依據(jù)計算機模擬計算所得結(jié)合能快速篩選聚合所需溶劑。

    2.2 紫外光譜分析

    2.2.1 溴氰菊酯、功能單體以及混合物的紫外吸收光譜差異分析

    根據(jù)紫外光譜測定的原理可知,對于相互之間沒有任何相互作用力的2種化合物,其混合物在某一波長下的吸光度應(yīng)該等于2種化合物分別于該波長下的吸光度之和(即理論吸光度)[19];其中,在分子印跡中不同功能單體與模板分子相互作用力的大小可由模板分子-功能單體的實際吸光度與理論值的差示吸光度的大小作初步的判斷。模板分子-功能單體實際值與理論值的最大吸收峰處的差示吸光度越大,預(yù)示著模板分子與功能單體之間相互作用越強,所制備的分子印跡聚合物的選擇性識別能力越高[20]。本研究采用紫外光譜法研究乙腈溶液中模板分子溴氰菊酯與3種不同功能單體(AA、AM和MAA)的相對結(jié)合能力,結(jié)果如圖2所示。

    由圖2可知,溴氰菊酯-功能單體(AA、 AM和MAA)混合物在吸收峰處吸光度的實際值均大于理論值,說明溴氰菊酯與3種功能單體之間均存在可相互結(jié)合的作用力。在吸收峰(200 nm)處實際值與理論值的差示吸光度?具體為:?DM–AA=0.174、?DM–AM=0.187、?DM–MAA= 0.132;其中DM-AM混合物的紫外差示吸光度最大。由此可知,溴氰菊酯與AM之間的相互作用更強,即以溴氰菊酯為模板分子時,相比較于AA和MAA,AM是較為理想的功能單體,這一結(jié)果與分子動力學(xué)模擬結(jié)果相一致。

    圖2 溴氰菊酯、丙烯酸、丙烯酰胺、甲基丙烯酸及其混合物的紫外吸收光譜

    2.2.2 預(yù)聚體系差示紫外光譜分析

    差示紫外光譜法通常被用于研究分子印跡預(yù)聚體系中主客體的相互作用,該方法通過測定模板分子-功能單體復(fù)合物的差示吸光度?并對其進行擬合,繪制?對單體濃度b(=1、2、3……代表各濃度梯度)的化學(xué)配位數(shù)冪b之比??b–n與?的曲線,考察曲線線性相關(guān)性并確定模板分子和功能單體間的化學(xué)配位比和結(jié)合常數(shù),進而指導(dǎo)和優(yōu)化MIPs的合成。溴氰菊酯與3種功能單體在不同濃度比例下的差示紫外光譜如圖3所示。

    由圖3可知,整體而言,隨著功能單體濃度的增加,復(fù)合物的吸光度呈現(xiàn)上升的趨勢,當功能單體濃度達到一定值時,上升的趨勢趨于緩慢;同時,特征吸收峰向長波長移動發(fā)生了紅移現(xiàn)象。由此可知,在一定范圍內(nèi)增加功能單體的比例,有助于溴氰菊酯與功能單體形成更多的選擇性識別位點;但當功能單體濃度過量時,可能出現(xiàn)功能單體自身結(jié)締作用,造成非特異性識別位點增多。再對比圖3a、圖3b和圖3c的紅移程度,其中DM-AA復(fù)合物特征峰的吸光度紅移了2 nm,DM-AM復(fù)合物特征峰吸光度紅移了4 nm,DM-MAA復(fù)合物特征峰吸光度紅移2 nm;由此可知DM-AM復(fù)合物特征峰的位移量相對較大,即溴氰菊酯與AM的作用力更加穩(wěn)固。

    圖3 溴氰菊酯與不同比例的功能單體預(yù)聚體差示紫外光譜

    參考已報道的文獻[21],模板分子與功能單體間相互作用的公式可推導(dǎo)整理如式(7)。

    式中為化學(xué)配位比,在0≥0前提下,可取=1、2、3……等整數(shù);0為功能單體濃度;0為模板溴氰菊酯濃度;?為溴氰菊酯-功能單體復(fù)合物的差示吸光度;為結(jié)合常數(shù);?ε=(εε),其中εε分別為溴氰菊酯和溴氰菊酯-功能單體復(fù)合物的摩爾吸光系數(shù);為比色皿厚度。

    依據(jù)圖3數(shù)據(jù)以??0–n對?作圖時發(fā)現(xiàn),當溴氰菊酯與功能單體濃度比為1∶2、1∶4和1∶6時,擬合曲線相關(guān)性較差。分析原因是式(7)成立的假設(shè)條件是0≥0,即功能單體的添加量遠大于模板分子添加量,顯然這3種濃度并不滿足上述假設(shè)條件,因此只對比例為1∶8~1∶14時的最強吸收峰值進行擬合計算,結(jié)果如表2所示。

    從表2可以看出,對于DM-AA體系,當化學(xué)配位比為1、2、3時擬合曲線的2均大于0.96,體現(xiàn)出良好的線性相關(guān)性,但三者中=3時的結(jié)合常數(shù)最大,說明1個溴氰菊酯最有可能與3個AA作用,形成較為穩(wěn)定的復(fù)合物。對于DM-AM體系,5條擬合直線的線性相關(guān)性總體較高,其中=2、3和4的擬合曲線的2均大于0.99,表明1個溴氰菊酯與2個、3個或4個AM分子形成的復(fù)合物相對穩(wěn)定。同時可以看出,當溴氰菊酯與AM的配位比=4時,所得復(fù)合物的結(jié)合常數(shù)最大,為6.614′105,說明DM與AM形成的最主要復(fù)合物為DM-4AM型。而對于DM-MAA體系,5條擬合直線的線性都比較差,說明溴氰菊酯與MAA難以形成穩(wěn)定的復(fù)合物。

    表2 溴氰菊酯與不同功能單體預(yù)聚體系的紫外光譜擬合參數(shù)

    注: ?為溴氰菊酯-功能單體復(fù)合物的差示吸光度。

    Note: ?is differential absorbance of deltamethrin - functional monomer complex.

    上述紫外光譜研究結(jié)果表明,相比較于AA和MAA,溴氰菊酯與AM可形成更為穩(wěn)定的復(fù)合物,相互間的作用力更強;且當=4,即溴氰菊酯與AM的化學(xué)配位比為1∶4時,形成的預(yù)聚體復(fù)合物較為穩(wěn)定。

    綜合動力學(xué)模擬和紫外光譜法對預(yù)聚體系的研究結(jié)果可知,當以溴氰菊酯為模板分子時,選擇AM為功能單體且其與溴氰菊酯的化學(xué)配位比為1∶4,且以乙腈為聚合溶劑時,可制備具有較好印跡效果的分子印跡聚合物。

    2.3 多因素正交試驗

    根據(jù)2.1和2.2部分對溴氰菊酯與不同功能單體預(yù)聚體系的研究結(jié)果,選取功能單體AM用量、交聯(lián)劑乙二醇二甲基丙烯酸酯用量、溶劑乙腈用量作為3個因素,利用L9(34)設(shè)計正交試驗,其因素和水平如表3所示。

    表3 正交試驗因素水平

    固定溴氰菊酯與丙烯酰胺的摩爾比為1∶4,以偶氮二異丁腈為引發(fā)劑,采用沉淀聚合法在所設(shè)計的因素水平下制備溴氰菊酯分子印跡聚合物。正交試驗結(jié)果及方差分析見表4和表5。由表4的極差分析可以看出,各因素對印跡聚合物的靜態(tài)吸附性能影響順序為功能單體用量>溶劑用量>交聯(lián)劑用量,得出的最佳試驗組合是221,即功能單體AM用量4 mmol,交聯(lián)劑乙二醇二甲基丙烯酸酯用量20 mmol即(溴氰菊酯)、丙烯酰胺與乙二醇二甲基丙烯酸酯的摩爾比為1∶4∶20,乙腈溶劑用量40 mL,恰好是表4中的7號試驗組。參照此條件,對所得的優(yōu)化配方進行3次平行試驗驗證,所制備聚合物的平均靜態(tài)吸附量為68.61 mg/g。表5的方差分析結(jié)果也表明功能單體和溶劑用量對聚合物吸附特性影響顯著,再次印證對功能單體和溶劑進行預(yù)聚體系篩選的必要性。

    表4 正交試驗方案及結(jié)果

    注:1~3分別表示各因素各水平下值的總和。

    Note:1-3refer to sum of values of each factor at each level.

    表5 正交試驗方差分析

    注:*代表差異顯著,<0.05。

    Note: * represents test factor is significant different,<0.05.

    2.4 分子印跡聚合物的表征

    2.4.1 溴氰菊酯MIPs的微觀結(jié)構(gòu)

    聚合物的表面結(jié)構(gòu)、孔徑和孔容等對其吸附和識別性能有很大的影響。分子印跡聚合物應(yīng)具備一定的三維孔穴結(jié)構(gòu),使識別位點固定在聚合物的骨架中確保目標分子能夠快速的與識別位點結(jié)合,同時能夠保持洗去模板分子后孔穴的大小和形狀,使孔穴維持對模板分子的再結(jié)合能力[22]。

    以上述優(yōu)化比例制備溴氰菊酯分子印跡聚合物的表觀結(jié)構(gòu)及孔特性如圖4和表6所示。

    由圖4可以看出,所制備的溴氰菊酯MIPs形狀規(guī)整,粒徑較小且相對均勻,存在緊密交錯的網(wǎng)絡(luò)狀微孔;NIPs表面則呈現(xiàn)不規(guī)則的塊狀,且存在團聚現(xiàn)象。從表6的比表面積和孔容參數(shù)可知,溴氰菊酯MIPs的孔徑小于NIPs,這與掃描電鏡圖的結(jié)果一致,但具有更大的比表面積和孔容,這種三維孔穴結(jié)構(gòu)有利于底物與結(jié)合位點接觸,使得溴氰菊酯MIPs比NIPs具有更好的吸附能力[23]。

    表6 分子印跡聚合物和非印跡聚合物的比表面積與孔容參數(shù)

    2.4.2 傅里葉紅外光譜分析

    紅外光譜法是根據(jù)分子內(nèi)部原子間的相對振動和分子轉(zhuǎn)動等信息來確定物質(zhì)分子結(jié)構(gòu)和鑒別化合物的方法,可用于分析分子印跡聚合物中模板分子與功能單體結(jié)合過程中基團的變化以及結(jié)合位點的位置[24]。采用紅外光譜技術(shù)測定溴氰菊酯標準品、NIPs、溴氰菊酯MIPs和吸附溴氰菊酯后的MIPs的紅外譜圖,結(jié)果見圖5。

    注:a為DM;b為NIPs;c為DM MIPs;d為吸附DM后的MIPs。

    由圖5可知,所制備溴氰菊酯MIPs在 3 565和 3 454 cm–1處有吸收峰,且譜帶較強,為丙烯酰胺的N-H的伸縮振動。2 987和2 955 cm–1處為烷烴C-H的伸縮振動吸收峰。1 730 cm–1處為C=O的伸縮振動。1 635 cm–1處為C=C的伸縮振動,說明可能存在未反應(yīng)完全的C=C。1 259和1 157 cm–1分別為C-O-C的對稱和不對稱伸縮振動,說明交聯(lián)劑與功能單體等進行了交聯(lián)聚合[25]。溴氰菊酯MIPs和NIPs的譜圖相似,證明聚合物成功聚合和洗脫[26]。吸附溴氰菊酯后的MIPs其吸收峰比MIPs強;同時在3 200~3 600 cm–1范圍內(nèi),吸收峰向短波數(shù)方向移動,可能存在酰胺的N-H的伸縮振動和溴氰菊酯中C=O的伸縮振動,且N-H和C=O可形成氫鍵[27];烷烴C-H的伸縮振動峰強度增大且存在泛頻峰;1 048和880 cm–1處的C-N伸縮振動和乙烯基中C-H的面外彎曲振動吸收峰均有加強,表明溴氰菊酯成功吸附于MIPs;溴氰菊酯中的C=O質(zhì)子受體與3 454 cm–1處的N-H作用形成的氫鍵為模板與聚合物的主要結(jié)合點。

    2.5 溴氰菊酯MIPs的特異性吸附

    以溴氰菊酯及其結(jié)構(gòu)類似物氯氰菊酯、高效氯氟氰菊酯和聯(lián)苯菊酯為底物,底物質(zhì)量濃度為100 mg/L,采用靜態(tài)平衡結(jié)合方法,在30 ℃溫度下吸附24 h后測定所制備聚合物對各底物的靜態(tài)吸附量、分離因子和印跡因子(見表7)。從表7可以看出,溴氰菊酯MIPs對溴氰菊酯的吸附量及靜態(tài)分配系數(shù)明顯高于其結(jié)構(gòu)類似物,且相應(yīng)的分離因子均大于3.2;同時,溴氰菊酯MIPs和NIPs對各物質(zhì)的印跡因子均大于1,其中溴氰菊酯對應(yīng)的印跡因子為2.269,說明溴氰菊酯MIPs對模板溴氰菊酯表現(xiàn)出較高的特異吸附性。結(jié)合溴氰菊酯及其類似物的結(jié)構(gòu)可知,氯氰菊酯和高效氯氟氰菊酯與溴氰菊酯的結(jié)構(gòu)相似,但是2個溴原子的改變使得它們與聚合物的三維孔穴結(jié)構(gòu)的匹配度下降且親和性降低等;聯(lián)苯菊酯的結(jié)構(gòu)與溴氰菊酯差異較大,同時存在空間位阻而導(dǎo)致溴氰菊酯MIPs對其吸附性能較差。試驗結(jié)果表明計算機和紫外光譜分析輔助設(shè)計所制備的溴氰菊酯MIPs具有較高的特異吸附性。

    表7 溴氰菊酯MIPs和NIPs的選擇性參數(shù)

    3 結(jié) 論

    本文以溴氰菊酯為模板,采用計算機分子模擬和紫外光譜分析法篩選適宜的聚合溶劑和功能單體,同時確定功能單體與溴氰菊酯的比例;以此為基礎(chǔ)采用沉淀聚合法制備溴氰菊酯分子印跡聚合物并通過正交試驗優(yōu)化制備方案,同時對其吸附特性進行研究??傻贸鲆韵陆Y(jié)論:

    1)采用分子動力學(xué)和紫外光譜法對預(yù)聚體進行研究,結(jié)合能和差示紫外光譜法結(jié)果表明:丙烯酰胺為功能單體,且當溴氰菊酯與功能單體的摩爾比為1∶4時,可在乙腈溶液中得到較為穩(wěn)定的復(fù)合物預(yù)聚體。

    2)以預(yù)聚體的研究結(jié)果為基礎(chǔ)通過正交試驗對合成方案進行優(yōu)化,結(jié)果表明:模板分子(溴氰菊酯)、功能單體(丙烯酰胺)與交聯(lián)劑(乙二醇二甲基丙烯酸酯)的摩爾比為1∶4∶20,對應(yīng)40 mL乙腈聚合溶劑,制備的聚合物的選擇性和印跡效率最佳,平行驗證試驗確認其飽和吸附量為68.61 mg/g,印跡因子為2.269。

    3)微觀結(jié)構(gòu)和紅外光譜表征結(jié)果表明所制備MIPs具有較規(guī)整的交錯微孔結(jié)構(gòu),溴氰菊酯中的羰基與酰胺基團可作用形成氫鍵,是主要的特異性識別位點。

    論文研究結(jié)果表明計算機分子模擬和紫外光譜法對預(yù)聚體系的研究有助于提高優(yōu)化試驗條件效率,所制備溴氰菊酯MIPs具有較高的特異吸附性能。

    [1] Liu Xingyu, Zhang Qiuping, Li Shibao, et al. Developmental toxicity and neurotoxicity of synthetic organic insecticides in zebrafish (): A comparative study of deltamethrin, acephate, and thiamethoxam[J]. Chemosphere, 2018, 199(1): 16-25.

    [2] Liu Pengyan, Liu Yujie, Liu Qingxue, et al. Photodegradation mechanism of deltamethrin and fenvalerate[J]. Journal of Environmental Sciences, 2010, 22(7): 1123-1128.

    [3] Rehman H, Ali M, Atif F, et al. The modulatory effect of deltamethrin on antioxidants in mice[J]. Clinica Chimica Acta, 2006, 369(3): 61-65.

    [4] 趙風(fēng)年. 基于磁性納米分子印跡技術(shù)的三唑類農(nóng)藥多殘留檢測方法研究[D]. 北京:中國農(nóng)業(yè)科學(xué)院,2017. Zhao Fengnian. Study on the Multi-residues Detection of Triazole Pesticides Based on Magnetic Nano Molecularly Imprinted Technology[D]. Beijing: Chinese Academy of Agricultural Sciences Dissertation, 2017. (in Chinese with English abstract)

    [5] Ge Shengguang, Lu Juanjuan, Ge Lei, et al. Development of a novel deltamethrin sensor based on molecularly imprinted silica nanospheres embedded CdTe quantum dots[J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2011, 79(5): 1704-1709.

    [6] Shi Xizhi, Liu Jinghua, Sun Aili, et al. Group-selective enrichment and determination of pyrethroid insecticides in aquaculture seawater via molecularly imprinted solid phase extraction coupled with gas chromatography-electron capture detection[J]. Journal Chromatogr A, 2012, 1227: 60-66.

    [7] Mariusz M, Teobald K, Piotr P W, et al. Computational modeling of molecularly imprinted polymers as a green approach to the development of novel analytical sorbents[J]. Trends in Analytical Chemistry, 2018, 98: 64-78.

    [8] Han Yi, Gu Lili, Zhang Mengxiao, et al. Computer-aided design of molecularly imprinted polymers for recognition of atrazine[J]. Computational and Theoretical Chemistry, 2017, 1121: 29-34.

    [9] 黃雨杰,朱秋勁,趙曉聯(lián).紫外光譜與密度泛函理論在亞精胺分子印跡預(yù)組裝體系中的研究[J]. 北京化工大學(xué)學(xué)報:自然科學(xué)版,2014(4):48-54. Huang Yujie, Zhu Qiujin, Zhao Xiaolian. Ultraviolet spectroscopy and density functional theory study of spermidine molecularly imprinted pre-assembly systems[J]. Journal of Beijing University of Chemical Technology: Natural Science, 2014(4): 48-54. (in Chinese with English abstract)

    [10] 尹小英,衷友泉,江一帆,等. 分子印跡聚合反應(yīng)中功能單體與模板分子間作用力的光譜分析[J]. 光譜學(xué)與光譜分析,2010(8):2211-2214. Yin Xiaoying, Zhong Youquan, Jiang Yifan, et al. The spectroscopy analysis of intermolecular interaction between the template molecule and functional monomer before polymerization[J]. Spectroscopy and Spectral Analysis, 2010(8): 2211-2214. (in Chinese with English abstract)

    [11] 朱曉蘭. 久效磷分子印跡聚合物的研究與應(yīng)用[D]. 合肥:中國科學(xué)技術(shù)大學(xué),2006. Zhu Xiaolan. Monocrotophos-specific Molecular Imprinted Polymers and its Applications for Analysis of Organophosphorus Pesticides[D]. Hefei: University of Science & Technology of China, 2006. (in Chinese with English abstract)

    [12] Farrington K, Magner E, Regan F. Predicting the performance of molecularly imprinted polymers: Selective extraction of caffeine by molecularly imprinted solid phase extraction[J]. Analytica Chimica Acta, 2006, 566(1): 60-68.

    [13] Yang Wenming, Liu Lukuan, Zhou Zhiping, et al. Rational preparation of dibenzothiophene-imprinted polymers by surface imprinting technique combined with atom transfer radical polymerization[J]. Applied Surface Science, 2013, 282: 809-819.

    [14] Chen Lingxin, Wang Xiaoyan, Lu Wenhui, et al. Molecular imprinting: Perspectives and applications[J]. Chemical Society Reviews, 2016, 45(8): 2137-2211.

    [15] Martins N, Carreiro E P, Locati A, et al. Design and development of molecularly imprinted polymers for the selective extraction of deltamethrin in olive oil: An integrated computational-assisted approach[J]. Journal of Chromatography A, 2015, 1409: 1-10.

    [16] 郭逸蓉. 計算機輔助設(shè)計擬除蟲菊酯類農(nóng)藥分子印跡聚合物研制[D]. 杭州:浙江大學(xué),2010. Guo Yirong. Computer-assisted Design for the Preparation of Molecularly Imprinted Polymers for Pyrethroid Pesticides[D]. Hangzhou: Zhejiang University, 2010. (in Chinese with English abstract)

    [17] 戴芳芳. 溴氰菊酯電化學(xué)快速檢測及分子印跡預(yù)組裝體系研究[D]. 重慶:西南大學(xué),2016. Dai Fangfang. Research on Deltamethri`s Rapid Detection by Electrochemical Merhod and Molecular Imprinting Preassembly System[D]. Chongqing: Southwest University, 2016. (in Chinese with English abstract)

    [18] Bakas I, Oujji N B, Moczko E, et al. Computational and experimental investigation of molecular imprinted polymers for selective extraction of dimethoate and its metabolite omethoate from olive oil [J]. Journal of Choromatography A, 2013, 1274(1): 13-18.

    [19] 朱淮武. 有機分子結(jié)構(gòu)波譜解析[M]. 北京:化學(xué)工業(yè)出版社,2005.

    [20] 高文惠,高林,劉博,等. 制備戊唑醇分子印跡聚合物前功能單體的選擇[J]. 中國食品學(xué)報,2015,15(6):173-177. Gao Wenhui, Gao Lin, Liu Bo, et al. Selection of functional monomer before the preparation of tebuconazole molecularly imprinted polymer[J]. Journal of Chinese Institute of Food

    Science and Technology, 2015, 15(6): 173-177. (in Chinese with English abstract)

    [21] 張孝剛,朱秋勁,胡萍. 三聚氰胺分子印跡預(yù)組裝體系紫外光譜研究[J]. 食品科學(xué),2011,32(21):128-132.Zhang Xiaogang, Zhu Qiujin, Hu Ping. Ultraviolet spectroscopy investigation into molecular imprinting pre-assembly system[J]. Food Science, 2011, 32(21): 128-132. (in Chinese with English abstract)

    [22] Cormack P A G, Elorza A Z. Molecularly imprinted polymers: Synthesis and characterisation[J]. Journal of Chromatography B Analytical Technologies in the Biomedical and Life Sciences, 2004, 804(1): 173-182.

    [23] 張連科,劉心宇,王維大,等. 油料作物秸稈生物炭對水體中鉛離子的吸附特性與機制[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(7):218-226. Zhang Lianke, Liu Xinyu, Wang Weida, et al. Characteristics and mechanism of lead adsorption from aqueous solutions by oil crops straw-derived biochar[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(7): 218-226. (in Chinese with English abstract)

    [24] Chrzanowska A M, Poliwoda A, Wieczorek P P. Characterization of particle morphology of biochanin A molecularly imprinted polymers and their properties as a potential sorbent for solid-phase extraction[J]. Materials Science and Engineering: C, 2015, 49: 793-798.

    [25] Wang Hongwu, Liu Yanqing, Yao Su, et al. Selective recognization of dicyandiamide in bovine milk by mesoporous silica SBA-15 supported dicyandiamide imprinted polymer based on surface molecularly imprinting technique[J]. Food Chemistry, 2018, 240: 1262-1267.

    [26] 龐思思. 沉淀制備雙酚A π-π堆積自組裝印跡聚合物及其特異性識別機理研究[J]. 高分子學(xué)報,2014(1):49-55. Pang Sisi. Synthesis of bisphenol-A stacking self-assembly imprinted polymer by precipitation polymerization and study on specificity recognition mechanism[J]. Acta Polymerica Sinica, 2014(1): 49-55. (in Chinese with English abstract)

    [27] 佘永新,王淼,史曉梅,等. 類特異性分子印跡聚合物微球的識別機理及其光譜學(xué)研究[J]. 光譜學(xué)與光譜分析, 2010,30(11):3052-3055. She Yongxin, Wang Miao, Shi Xiaomei, et al. Spectroscopy study of binding mechanisms and molecular recognition of class-specific molecularly imprinted polymer beads[J]. Spectroscopy and Spectral Analysis, 2010, 30(11): 3052-3055. (in Chinese with English abstract)

    Screening of molecularly imprinted pre-assembly system for detection of deltamethrin pesticide residues and its specific adsorption properties

    Li Lu, Zhou Liumei, Xie Xinan, Li Yan

    (510640,)

    To efficiently prepare decamethrin molecular (DM) imprinted polymers, the molecular dynamics simulation and ultraviolet spectrum were involved to study the preassembled system of DM and functional monomers in this paper. The binding energy (?) between DM and three different functional monomers (acrylic acid, AA; acrylamide, AM; methacrylic acid, MAA) in different solvents (n-hexane, chloroform, acetone and acetonitrile) were calculated by molecular dynamics simulation. Results showed that the AM monomer and acetonitrile solvent had the highest stabilization energy for the pre-polymerization adducts. Additionally, the polymer obtained from the molar ratio of DM to AM at 1:4 to 1:6 may be more stable. Ultraviolet spectrum analysis results showed that there was strong interaction between DM and AM, the red shift degree of ultraviolet spectrum maximum absorption peak followed the order: AM > AA = MAA, and then indicated that AM was most suitable monomers for DM molecular imprinted polymers, which was consistent with the predictions based on the computer simulation. Differential ultraviolet spectrum showed that when acetonitrile used as solvent, one DM molecule could be combined with AA monomer in the ratio of 1:3 or AM monomer in the ratio of 1:4. However, it was difficult for DM to form stable complex with MAA because the fitting correlation coefficient of the equation was less than 0.86. Best performance was observed when chemical coordination ratio () between DM and AM was equal to 4, at this time, the fitting correlation coefficientof the equation was 0.997 3, and the binding constantwas 6.614×105. Based on the above computer simulation and ultraviolet spectrum analysis results, the molecularly imprinted polymers of DM were synthesized by the precipitation polymerization method with the AM as the functional monomer, and its synthetic formula was optimized through orthogonal experiment design. Then, the optimum imprinting conditions were established: 1 mmol DM as the template molecule, 4 mmol AM as the functional monomers, 40 mL acetonitrile as the solvent, 20 mmol ethylene dimethacrylate (EGDMA) as the cross-linker. Using these conditions, the DM molecular imprinted polymers had bigger saturated adsorption capacity of 68.61 mg/g than other conditions. The prepared molecular imprinted polymers were characterized by scanning electron microscopy (SEM), brunauer-emmett-teller (BET) and Fourier transform infrared spectroscopy (FTIR). SEM demonstrates that the DM molecular imprinted polymers formed staggered porous structure. The data of BET revealed that the DM molecular imprinted polymers had smaller pore size (15.83 nm) than non-imprinted polymers (18.84 nm), but with a larger specific surface area (95.58 m2/g) and pore volume (0.578 5 cm3/g). The FTIR results confirmed the successful preparation of the DM molecular imprinted polymers. And then the adsorption selectivity of prepared polymers on deltamethrin and its analogues cypermethrin, lambda-cyhalothrin and bifenthrin were studied. The experimental results showed that the static equilibrium adsorption capacity and static distribution coefficient of DM was much higher than that of other analogues. Its imprinting factor (IF) was 2.269, and the corresponding separation factor α was more than 3.2, which has exhibited obvious specific identification. All above results showed that computer simulation and ultraviolet spectrum analysis could help to design the better performance molecular imprinted polymers.

    pesticide; molecular dynamics; ultraviolet visible spectroscopy

    2018-06-09

    2018-10-19

    廣州市科技計劃項目(20160701140); 國家自然科學(xué)基金項目(31501471)

    李 璐,副教授,主要從事功能材料的制備與應(yīng)用研究。 Email:lulu_lee@scau.edu.cn

    10.11975/j.issn.1002-6819.2019.01.033

    TQ424.3

    A

    1002-6819(2019)-01-0269-09

    李 璐,周劉梅,解新安,李 雁.溴氰菊酯農(nóng)藥殘留檢測的分子印跡預(yù)聚體系篩選及吸附性能[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(1):269-277. doi:10.11975/j.issn.1002-6819.2019.01.033 http://www.tcsae.org

    Li Lu, Zhou Liumei, Xie Xinan, Li Yan.Screening of molecularly imprinted pre-assembly system for detection of deltamethrin pesticide residues and its specific adsorption properties[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(1): 269-277. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.01.033 http://www.tcsae.org

    猜你喜歡
    溴氰菊酯印跡菊酯
    馬 浩
    陶瓷研究(2022年3期)2022-08-19 07:15:18
    走進大美滇西·探尋紅色印跡
    云南畫報(2021年10期)2021-11-24 01:06:56
    溴氰菊酯的急性毒性效應(yīng)與毒性機理研究進展
    5%氯氟醚菊酯·高效氟氯氰菊酯懸浮劑正相高效液相色譜分析
    歐盟擬修訂高效氯氟氰菊酯在芹菜、茴香和大米中的最大殘留限量
    食品與機械(2019年1期)2019-03-30 01:14:36
    一株溴氰菊酯降解菌的分離鑒定及降解特性研究
    成長印跡
    亞臨界流體萃取技術(shù)在茶葉降低菊酯類農(nóng)殘方面的應(yīng)用
    印跡
    中國攝影(2014年12期)2015-01-27 13:57:04
    溴氰菊酯降解菌的分離與鑒定及其降解特性
    亚洲欧美成人综合另类久久久| 在线观看三级黄色| 欧美人与性动交α欧美精品济南到| 国产精品一区二区在线观看99| 1024香蕉在线观看| 亚洲人成网站在线观看播放| 波多野结衣av一区二区av| 丁香六月天网| tube8黄色片| 亚洲国产看品久久| 在线观看人妻少妇| 中文字幕亚洲精品专区| 中文字幕制服av| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠躁躁| 麻豆精品久久久久久蜜桃| 99热全是精品| 性色av一级| 婷婷色综合www| 天堂8中文在线网| xxxhd国产人妻xxx| 制服人妻中文乱码| 久久ye,这里只有精品| 搡老乐熟女国产| av国产精品久久久久影院| 久久韩国三级中文字幕| 啦啦啦在线免费观看视频4| av国产久精品久网站免费入址| 精品卡一卡二卡四卡免费| av卡一久久| av国产久精品久网站免费入址| 无限看片的www在线观看| 成年人午夜在线观看视频| 欧美人与性动交α欧美软件| 天堂8中文在线网| 亚洲综合精品二区| 高清欧美精品videossex| 色婷婷久久久亚洲欧美| 中文乱码字字幕精品一区二区三区| 老司机在亚洲福利影院| 亚洲精品国产区一区二| 精品亚洲乱码少妇综合久久| videos熟女内射| 少妇精品久久久久久久| 老汉色∧v一级毛片| 啦啦啦在线免费观看视频4| 久久精品国产综合久久久| 老司机深夜福利视频在线观看 | 精品久久蜜臀av无| 最新在线观看一区二区三区 | 中国三级夫妇交换| av又黄又爽大尺度在线免费看| 一边亲一边摸免费视频| 国产片内射在线| 一区二区三区四区激情视频| 午夜91福利影院| 2021少妇久久久久久久久久久| 侵犯人妻中文字幕一二三四区| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩一区二区三区在线| 久久人人97超碰香蕉20202| 日日啪夜夜爽| 久久国产精品男人的天堂亚洲| 青春草视频在线免费观看| 午夜免费男女啪啪视频观看| 久久精品国产亚洲av涩爱| 国产av码专区亚洲av| 一级,二级,三级黄色视频| 欧美人与性动交α欧美精品济南到| 日韩精品有码人妻一区| 亚洲免费av在线视频| 国产视频首页在线观看| 黑人欧美特级aaaaaa片| 丝袜美足系列| 丝袜脚勾引网站| 久久精品久久久久久久性| 伊人久久国产一区二区| 午夜免费观看性视频| 在线精品无人区一区二区三| 老司机影院毛片| 国产精品国产三级专区第一集| 精品国产一区二区久久| 亚洲欧洲日产国产| 久久青草综合色| 五月天丁香电影| 韩国av在线不卡| 女的被弄到高潮叫床怎么办| tube8黄色片| 国产精品二区激情视频| 久久女婷五月综合色啪小说| 黄频高清免费视频| 男女边摸边吃奶| 婷婷色综合大香蕉| 日韩精品有码人妻一区| 亚洲免费av在线视频| 搡老乐熟女国产| 中文字幕高清在线视频| 亚洲四区av| 日韩中文字幕欧美一区二区 | 免费黄频网站在线观看国产| 热99久久久久精品小说推荐| 成年人免费黄色播放视频| 国产免费福利视频在线观看| av在线观看视频网站免费| 久久精品久久久久久噜噜老黄| 欧美中文综合在线视频| 久久久久久久久免费视频了| 王馨瑶露胸无遮挡在线观看| 久久久久网色| 建设人人有责人人尽责人人享有的| 国产在线一区二区三区精| 男女无遮挡免费网站观看| 伦理电影大哥的女人| 午夜福利网站1000一区二区三区| 国产国语露脸激情在线看| 大片电影免费在线观看免费| 久久ye,这里只有精品| 免费在线观看视频国产中文字幕亚洲 | 在线观看www视频免费| 高清在线视频一区二区三区| 亚洲国产精品999| 国产女主播在线喷水免费视频网站| av免费观看日本| 啦啦啦在线观看免费高清www| 啦啦啦中文免费视频观看日本| 免费日韩欧美在线观看| 999精品在线视频| 午夜影院在线不卡| 亚洲国产日韩一区二区| 欧美精品人与动牲交sv欧美| 欧美日韩福利视频一区二区| 亚洲欧洲日产国产| 精品少妇内射三级| 99热全是精品| 人人妻,人人澡人人爽秒播 | 侵犯人妻中文字幕一二三四区| 国产欧美日韩一区二区三区在线| 久久99热这里只频精品6学生| 美女福利国产在线| 亚洲五月色婷婷综合| 丝袜喷水一区| 国产探花极品一区二区| 国产成人精品久久二区二区91 | 欧美精品高潮呻吟av久久| 新久久久久国产一级毛片| 国产不卡av网站在线观看| videosex国产| 久久天躁狠狠躁夜夜2o2o | 日本色播在线视频| 午夜久久久在线观看| 亚洲av国产av综合av卡| 日韩一本色道免费dvd| 久久久久网色| 国产国语露脸激情在线看| 黄色一级大片看看| 汤姆久久久久久久影院中文字幕| 亚洲成国产人片在线观看| 9191精品国产免费久久| 纯流量卡能插随身wifi吗| 日韩电影二区| 国产精品av久久久久免费| 香蕉丝袜av| 国产国语露脸激情在线看| 我的亚洲天堂| 97人妻天天添夜夜摸| 日韩精品免费视频一区二区三区| 黄色一级大片看看| 亚洲一区二区三区欧美精品| 欧美日韩成人在线一区二区| 午夜激情久久久久久久| 日本爱情动作片www.在线观看| 制服人妻中文乱码| 久久久久久久精品精品| 亚洲 欧美一区二区三区| 亚洲精品第二区| 亚洲成人国产一区在线观看 | 一级毛片黄色毛片免费观看视频| 在线观看免费日韩欧美大片| 成人国产av品久久久| 亚洲国产欧美在线一区| 久久精品久久精品一区二区三区| 精品少妇一区二区三区视频日本电影 | 国产一区二区在线观看av| 亚洲 欧美一区二区三区| 一区二区三区激情视频| 欧美黑人精品巨大| 国产福利在线免费观看视频| 成人国语在线视频| 午夜福利视频精品| 激情视频va一区二区三区| 亚洲精品美女久久av网站| 亚洲av中文av极速乱| 久久久久国产一级毛片高清牌| 亚洲欧美清纯卡通| 亚洲国产精品国产精品| 免费少妇av软件| 亚洲精品美女久久久久99蜜臀 | 美女脱内裤让男人舔精品视频| 亚洲国产av新网站| 久久久久久久精品精品| 啦啦啦啦在线视频资源| 国产成人一区二区在线| 视频在线观看一区二区三区| av网站在线播放免费| 欧美人与性动交α欧美软件| 极品少妇高潮喷水抽搐| 午夜福利免费观看在线| 婷婷色综合大香蕉| xxxhd国产人妻xxx| 麻豆精品久久久久久蜜桃| 久久久精品免费免费高清| 免费观看av网站的网址| 日韩中文字幕视频在线看片| 久久精品亚洲熟妇少妇任你| 亚洲欧美色中文字幕在线| 99国产综合亚洲精品| 精品国产一区二区久久| 国产成人免费无遮挡视频| 久久久久久人人人人人| 精品人妻熟女毛片av久久网站| 国产一区二区激情短视频 | 国语对白做爰xxxⅹ性视频网站| 欧美亚洲日本最大视频资源| 高清视频免费观看一区二区| 男的添女的下面高潮视频| 久久久久国产一级毛片高清牌| www日本在线高清视频| 久久国产精品男人的天堂亚洲| 国产男女内射视频| 菩萨蛮人人尽说江南好唐韦庄| 国产男女内射视频| 欧美精品亚洲一区二区| 国产av一区二区精品久久| 香蕉国产在线看| 搡老乐熟女国产| 在线免费观看不下载黄p国产| 热99久久久久精品小说推荐| 国产极品天堂在线| 丰满迷人的少妇在线观看| 啦啦啦啦在线视频资源| 亚洲精品美女久久久久99蜜臀 | 精品一区二区三区四区五区乱码 | 精品第一国产精品| 丰满饥渴人妻一区二区三| 欧美黑人欧美精品刺激| 中文欧美无线码| 天堂8中文在线网| 国产成人欧美| 老汉色av国产亚洲站长工具| 如日韩欧美国产精品一区二区三区| 两个人看的免费小视频| 一区二区三区激情视频| 亚洲,欧美精品.| 街头女战士在线观看网站| 亚洲精品久久午夜乱码| 欧美成人精品欧美一级黄| 国产免费现黄频在线看| 最新在线观看一区二区三区 | 久久久久久久久免费视频了| 亚洲精品在线美女| 两个人免费观看高清视频| 1024视频免费在线观看| 国产精品欧美亚洲77777| 国产精品免费视频内射| 国产成人一区二区在线| 99精品久久久久人妻精品| 热re99久久国产66热| 亚洲国产精品国产精品| 51午夜福利影视在线观看| 大码成人一级视频| 国产 精品1| 深夜精品福利| 久久这里只有精品19| 欧美乱码精品一区二区三区| 久久国产精品男人的天堂亚洲| 亚洲精品aⅴ在线观看| 国产精品欧美亚洲77777| 超碰97精品在线观看| videosex国产| av又黄又爽大尺度在线免费看| 欧美日韩精品网址| 人人妻人人爽人人添夜夜欢视频| 69精品国产乱码久久久| 91精品伊人久久大香线蕉| 国产在线一区二区三区精| 亚洲欧美一区二区三区黑人| 久久精品久久久久久噜噜老黄| 精品卡一卡二卡四卡免费| 欧美精品高潮呻吟av久久| 一级片免费观看大全| 欧美久久黑人一区二区| 国产一级毛片在线| 人人妻人人澡人人爽人人夜夜| 成人国语在线视频| 亚洲欧美中文字幕日韩二区| 国产毛片在线视频| 国产一区二区三区av在线| 国产av国产精品国产| 伊人久久国产一区二区| 肉色欧美久久久久久久蜜桃| 亚洲欧洲日产国产| 秋霞在线观看毛片| 亚洲,欧美,日韩| 99热网站在线观看| 女性被躁到高潮视频| 国产精品国产av在线观看| 欧美日韩福利视频一区二区| 夜夜骑夜夜射夜夜干| 久久婷婷青草| 老鸭窝网址在线观看| 免费观看性生交大片5| 日日啪夜夜爽| 亚洲第一区二区三区不卡| 久久午夜综合久久蜜桃| 亚洲精品中文字幕在线视频| 99国产精品免费福利视频| 亚洲伊人久久精品综合| 9热在线视频观看99| 日韩欧美一区视频在线观看| 午夜av观看不卡| 大片免费播放器 马上看| 秋霞伦理黄片| av免费观看日本| 老司机影院成人| 久久精品aⅴ一区二区三区四区| 国产黄色免费在线视频| 热re99久久国产66热| 免费少妇av软件| 欧美亚洲 丝袜 人妻 在线| 看免费av毛片| 国产成人免费无遮挡视频| 天天影视国产精品| 精品国产乱码久久久久久男人| 伦理电影免费视频| 中文字幕亚洲精品专区| 国产毛片在线视频| 亚洲av男天堂| 亚洲欧美一区二区三区国产| 韩国精品一区二区三区| 看免费成人av毛片| 亚洲综合色网址| 秋霞在线观看毛片| 久久久久精品人妻al黑| 国产一区亚洲一区在线观看| 妹子高潮喷水视频| 亚洲欧美一区二区三区黑人| 久久精品久久精品一区二区三区| 久久久久久人妻| 国产又色又爽无遮挡免| 99热国产这里只有精品6| 无限看片的www在线观看| 久久久欧美国产精品| 看非洲黑人一级黄片| 成年人午夜在线观看视频| 人人妻人人爽人人添夜夜欢视频| 99久久精品国产亚洲精品| a 毛片基地| a级片在线免费高清观看视频| 精品一区二区三区四区五区乱码 | 国产亚洲午夜精品一区二区久久| 国产精品国产三级国产专区5o| 久久av网站| 天天躁夜夜躁狠狠久久av| 精品一区在线观看国产| 99久久人妻综合| 国产成人精品无人区| 天美传媒精品一区二区| 亚洲国产毛片av蜜桃av| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美变态另类bdsm刘玥| 青草久久国产| 狠狠婷婷综合久久久久久88av| 久久ye,这里只有精品| 久久久久久人人人人人| 男的添女的下面高潮视频| 黄色怎么调成土黄色| 久久婷婷青草| 黄片无遮挡物在线观看| av国产久精品久网站免费入址| 无遮挡黄片免费观看| 成人手机av| 2021少妇久久久久久久久久久| 少妇精品久久久久久久| 一区二区三区精品91| 丰满少妇做爰视频| 18禁国产床啪视频网站| 亚洲欧美成人综合另类久久久| 亚洲自偷自拍图片 自拍| 国产一区二区三区综合在线观看| 久久久精品国产亚洲av高清涩受| 成年人午夜在线观看视频| 国产在线视频一区二区| 亚洲欧洲日产国产| 国产日韩欧美视频二区| 18禁国产床啪视频网站| 一级片'在线观看视频| 久久久久久久大尺度免费视频| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 欧美激情极品国产一区二区三区| 成年女人毛片免费观看观看9 | 亚洲av国产av综合av卡| 大片免费播放器 马上看| 国产又色又爽无遮挡免| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产区一区二| 国产在线视频一区二区| 国产一区二区在线观看av| 日韩 亚洲 欧美在线| 19禁男女啪啪无遮挡网站| 亚洲欧洲日产国产| 日日啪夜夜爽| 亚洲精品,欧美精品| 免费日韩欧美在线观看| 欧美黑人欧美精品刺激| 熟女少妇亚洲综合色aaa.| 伊人久久国产一区二区| 国产一区二区三区av在线| 老司机在亚洲福利影院| 国产精品.久久久| 国产视频首页在线观看| 丝袜喷水一区| 69精品国产乱码久久久| 9热在线视频观看99| 国产淫语在线视频| 91aial.com中文字幕在线观看| 99久国产av精品国产电影| 国产欧美亚洲国产| 久久久久精品国产欧美久久久 | 亚洲国产精品一区二区三区在线| 青春草亚洲视频在线观看| 嫩草影视91久久| 亚洲美女黄色视频免费看| 久久精品亚洲熟妇少妇任你| 国产在线免费精品| 一级片'在线观看视频| 桃花免费在线播放| 韩国精品一区二区三区| 男女下面插进去视频免费观看| 日本爱情动作片www.在线观看| 久久精品国产a三级三级三级| 男女下面插进去视频免费观看| 中文字幕高清在线视频| 亚洲,欧美精品.| 国产男女内射视频| videos熟女内射| 亚洲熟女精品中文字幕| 久久人人爽人人片av| 最近最新中文字幕大全免费视频 | 亚洲婷婷狠狠爱综合网| 欧美日韩亚洲综合一区二区三区_| 成人亚洲欧美一区二区av| 欧美久久黑人一区二区| 亚洲欧美色中文字幕在线| 日韩精品有码人妻一区| 黄网站色视频无遮挡免费观看| 国产女主播在线喷水免费视频网站| av卡一久久| 午夜福利,免费看| 久久久精品94久久精品| 青春草视频在线免费观看| 狂野欧美激情性xxxx| 亚洲av成人不卡在线观看播放网 | 搡老岳熟女国产| 日韩一卡2卡3卡4卡2021年| 纯流量卡能插随身wifi吗| 亚洲色图 男人天堂 中文字幕| 成年人午夜在线观看视频| 日韩熟女老妇一区二区性免费视频| 国产亚洲最大av| 成人午夜精彩视频在线观看| 色综合欧美亚洲国产小说| 99久久人妻综合| 国产极品粉嫩免费观看在线| 999精品在线视频| 又大又爽又粗| 777久久人妻少妇嫩草av网站| 男女免费视频国产| 人妻人人澡人人爽人人| 国产高清不卡午夜福利| 中文字幕av电影在线播放| 欧美黑人欧美精品刺激| 亚洲欧洲精品一区二区精品久久久 | 九九爱精品视频在线观看| 制服诱惑二区| 国产精品国产三级国产专区5o| 久久久亚洲精品成人影院| 不卡av一区二区三区| 99国产综合亚洲精品| 毛片一级片免费看久久久久| 久久久久久久久免费视频了| 午夜福利视频在线观看免费| 国产欧美日韩综合在线一区二区| 亚洲av综合色区一区| 精品福利永久在线观看| 又黄又粗又硬又大视频| 久久国产精品大桥未久av| 亚洲,一卡二卡三卡| 精品少妇黑人巨大在线播放| 国产一区二区 视频在线| 18禁国产床啪视频网站| 多毛熟女@视频| 美女中出高潮动态图| 最新在线观看一区二区三区 | 国产片特级美女逼逼视频| 天堂中文最新版在线下载| 午夜91福利影院| 最近中文字幕2019免费版| 国产成人啪精品午夜网站| 人人妻,人人澡人人爽秒播 | 人成视频在线观看免费观看| 国产在视频线精品| 精品国产露脸久久av麻豆| 高清不卡的av网站| 欧美黄色片欧美黄色片| 青青草视频在线视频观看| 男女床上黄色一级片免费看| 黄色怎么调成土黄色| 亚洲精品美女久久av网站| 一二三四中文在线观看免费高清| 日本黄色日本黄色录像| 久久久久久久精品精品| 亚洲精品美女久久久久99蜜臀 | bbb黄色大片| 国产爽快片一区二区三区| 少妇被粗大猛烈的视频| 亚洲三区欧美一区| 超色免费av| 久久久精品区二区三区| 国产黄频视频在线观看| 在线观看www视频免费| av网站免费在线观看视频| 国产精品免费视频内射| 这个男人来自地球电影免费观看 | 老鸭窝网址在线观看| 亚洲欧美精品综合一区二区三区| 色精品久久人妻99蜜桃| 九色亚洲精品在线播放| 视频在线观看一区二区三区| 亚洲精品第二区| 国产精品国产三级专区第一集| 一区二区三区激情视频| 日本色播在线视频| 999久久久国产精品视频| 亚洲国产av影院在线观看| 高清在线视频一区二区三区| 国产成人欧美在线观看 | 少妇 在线观看| av在线观看视频网站免费| 韩国高清视频一区二区三区| 欧美精品一区二区免费开放| 美女视频免费永久观看网站| 亚洲熟女精品中文字幕| 午夜精品国产一区二区电影| 午夜老司机福利片| 丝袜喷水一区| 精品久久蜜臀av无| 搡老岳熟女国产| 国产亚洲午夜精品一区二区久久| 午夜免费男女啪啪视频观看| 久久午夜综合久久蜜桃| 精品少妇黑人巨大在线播放| 日本猛色少妇xxxxx猛交久久| 999久久久国产精品视频| 精品久久久久久电影网| 亚洲精品国产av蜜桃| av在线观看视频网站免费| 国产 精品1| 亚洲av国产av综合av卡| 超碰97精品在线观看| 国产熟女午夜一区二区三区| 天天躁夜夜躁狠狠躁躁| 99热网站在线观看| 青春草视频在线免费观看| 久久ye,这里只有精品| 亚洲在久久综合| 91精品国产国语对白视频| 欧美激情高清一区二区三区 | 男的添女的下面高潮视频| 国产乱来视频区| 亚洲熟女精品中文字幕| 久久久久久人妻| 在线观看三级黄色| 久久毛片免费看一区二区三区| 色视频在线一区二区三区| 看十八女毛片水多多多| 亚洲欧美清纯卡通| 欧美精品人与动牲交sv欧美| 韩国av在线不卡| av福利片在线| 国产成人系列免费观看| 国产片内射在线| 国产欧美亚洲国产| 一本大道久久a久久精品| 国产福利在线免费观看视频| 国产乱来视频区| 美女主播在线视频| 精品人妻一区二区三区麻豆| 国产亚洲欧美精品永久| kizo精华| 午夜福利视频精品| 国产免费福利视频在线观看| 最近2019中文字幕mv第一页| 亚洲天堂av无毛| 涩涩av久久男人的天堂| av又黄又爽大尺度在线免费看| 久久毛片免费看一区二区三区| 最近中文字幕高清免费大全6| 美女脱内裤让男人舔精品视频| 水蜜桃什么品种好| 国产熟女午夜一区二区三区| 97在线人人人人妻| 日本猛色少妇xxxxx猛交久久| 男女床上黄色一级片免费看| 亚洲国产欧美网|