999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

CuI催化溴代苯與L-纈氨酸的反應(yīng)機(jī)理的研究

2014-08-07 11:38:18潘曉曉鄒良明岳文喜李來(lái)才

潘曉曉, 鄒良明, 岳文喜, 李來(lái)才?

(1.四川師范大學(xué)化學(xué)與材料科學(xué)學(xué)院,四川成都610066; 2.四川工商職業(yè)技術(shù)學(xué)院,四川都江堰611837)

有機(jī)含氮化合物是一類非常重要的化合物,C—N鍵不僅在許多天然的富有生理活性的物質(zhì)中存在,在現(xiàn)代社會(huì)工業(yè)化生產(chǎn)出的合成物中也常常出現(xiàn).作為不可或缺的重要組分在一些結(jié)構(gòu)簡(jiǎn)單的含氮化合物基礎(chǔ)上構(gòu)建新的C—N鍵,對(duì)合成含有這種結(jié)構(gòu)單元的天然產(chǎn)物、藥物、農(nóng)藥、材料等非常重要,因此這一直是有機(jī)化學(xué)家重點(diǎn)研究的課題之一[1].C—N交叉偶聯(lián)起先是以Ni、Pd為催化劑[2],近年來(lái)人們發(fā)現(xiàn)Cu有很好的催化性能,并且具有更廉價(jià)更環(huán)保等優(yōu)勢(shì),銅催化的C—N交叉偶聯(lián)反應(yīng)已成為有機(jī)合成中構(gòu)建C—N鍵的主要途徑[3-4],如今經(jīng)過(guò)科學(xué)家們的不斷改進(jìn),銅催化的C—N交叉偶聯(lián)反應(yīng)所需的條件已變得更加溫和、簡(jiǎn)單,使其應(yīng)用范圍進(jìn)一步擴(kuò)大[5].然而,在這一領(lǐng)域中仍然存在著一些有待于更進(jìn)一步研究的課題.近年來(lái)對(duì)偶聯(lián)反應(yīng)機(jī)理的理論研究工作開(kāi)展得比較活躍[6-7],H.Z.Yu等[8]研究了亞銅催化的碳雜偶聯(lián)反應(yīng)機(jī)理,我們課題組研究過(guò)銅催化下N-(鄰氯苯基)苯甲酰胺分子內(nèi)O-芳基化偶聯(lián)反應(yīng)[9].它在對(duì)催化劑的優(yōu)化、反應(yīng)機(jī)理的描述和發(fā)現(xiàn)新的反應(yīng)通道方面有著很大的優(yōu)勢(shì),理論研究偶聯(lián)反應(yīng)微觀機(jī)理能夠更好地了解催化劑作用的本質(zhì),解釋實(shí)驗(yàn)現(xiàn)象.

2009年Masanari團(tuán)隊(duì)的F.Ullmann[10]研究發(fā)現(xiàn),用CuI催化劑能夠催化溴代苯和L-纈氨酸發(fā)生C—N偶聯(lián)反應(yīng).本文對(duì)這一典型的C—N偶聯(lián)反應(yīng)的微觀反應(yīng)特征進(jìn)行了詳細(xì)的機(jī)理研究,以了解這一合成反應(yīng)的微觀歷程,為相關(guān)的實(shí)驗(yàn)及新問(wèn)題的解決提供一些參考.用CuI催化劑能夠催化溴代苯和L-纈氨酸發(fā)生C—N偶聯(lián)反應(yīng),脫去HX,形成N-芳基-α-氨基酸,反應(yīng)式如圖1所示.

1 計(jì)算方法

采用密度泛函理論(DFT)[11-12],在B3LYP/6-311+G(d)基組水平上(Cu采用贗勢(shì)基組LanL2DZ,鹵素Cl、I采用贗勢(shì)機(jī)組MIDIX)對(duì)CuX催化L-纈氨酸(N-芳基-α-氨基酸)與溴苯反應(yīng)生成N-芳基-α-氨基酸過(guò)程中所有反應(yīng)物、中間體、過(guò)渡態(tài)和產(chǎn)物進(jìn)行了優(yōu)化,并在相同基組水平下進(jìn)行了頻率計(jì)算,各過(guò)渡態(tài)都有唯一虛頻,通過(guò)振動(dòng)分析進(jìn)一步確認(rèn)了過(guò)渡態(tài)的真實(shí)性.通過(guò)反應(yīng)中各物質(zhì)相對(duì)能量比較,分析了該反應(yīng)的速率決定步驟以及主要反應(yīng)通道.同時(shí)在此基組水平上運(yùn)用自然鍵軌道[13](NBO)分析方法分析了分子軌道間相互作用;用AIM 2000程序包[14]計(jì)算了相應(yīng)的成鍵臨界點(diǎn)(BCP)和成環(huán)臨界點(diǎn)(RCP)電荷密度,分析了成鍵特征,所有計(jì)算均采用Gaussian 09程序[15].

圖1 CuI催化的L-纈氨酸與溴代苯C—N偶聯(lián)反應(yīng)Fig.1 Copper-catalyzed C—N coupling reaction of aryl halides with α-amino acids

圖2 CuI催化L-纈氨酸與溴苯的反應(yīng)機(jī)理Fig.2 Detailed route of CuI-catalyzed reaction between aryl halides and α-amino acids

2 計(jì)算結(jié)果與討論

圖2所示為CuI催化L-纈氨酸(α-氨基酸)與溴苯反應(yīng)生成N-芳基-α-氨基酸的反應(yīng)機(jī)理.由反應(yīng)物RE1與CuI結(jié)合開(kāi)始得到了2條不同的反應(yīng)通道,其中路徑IA是CuI連到了α-氨基酸的羧基氧上開(kāi)始的;而路徑IB則是CuI連到了α-氨基酸的氨基上開(kāi)始的.下面分別介紹這2條路徑的反應(yīng)機(jī)理.

2.1 路徑IA反應(yīng)機(jī)理分析表1為路徑IA中各化合物的能量及校正后各物質(zhì)的相關(guān)能量.該反應(yīng)機(jī)理中所有反應(yīng)物、中間體、過(guò)渡態(tài)和產(chǎn)物的優(yōu)化分子構(gòu)型以及部分鍵臨界點(diǎn)(BCP)電荷密度(圖中圓括號(hào)內(nèi))如圖3所示.

表1 反應(yīng)路徑IA中各駐點(diǎn)能量E、相對(duì)能量Erel和過(guò)渡態(tài)的振動(dòng)頻率vTable 1 Energies E,relative energies Erel,and frequencies v of the compounds in the reactions for IA

在路徑IA中,反應(yīng)的第一步是催化劑CuI中的Cu原子與反應(yīng)物α-氨基酸中的O(7)活化絡(luò)合,得到中間體MA1.MA1中Cu—I鍵從反應(yīng)物的鍵長(zhǎng)0.230 1 nm增長(zhǎng)到0.237 2 nm,電荷密度從0.092 3 a.u.減少至0.071 4 a.u.,Cu—O鍵長(zhǎng)為0.198 9 nm,電荷密度為0.084 1 a.u.,H—O(1)鍵幾乎沒(méi)有變化.表明Cu與O有成鍵,RE1與CuI生成MA1.NBO計(jì)算分析得知,在中間體MA1中,BD(1)O(7)—H(17)→LP?(8)Cu、BD(1)C(4)—O(7)→LP?(8)Cu間的二階穩(wěn)定化能分別為35.55和20.31 kJ/mol,表明它們之間存在軌道間相互作用.從表1能量數(shù)據(jù)可以看出,中間體MA1的能量比α-氨基酸與催化劑CuI能量之和低82.89 kJ/mol,表明中間體MA1易形成并穩(wěn)定存在.MA1到TSA1過(guò)程中,H原子從O(1)上往I上遷移.H—I鍵長(zhǎng)為0.166 0 nm,電荷密度為0.134 3 a.u.;Cu—I鍵略有減弱,鍵長(zhǎng)增長(zhǎng)為0.247 9 nm.Cu—O(1)鍵長(zhǎng)為0.190 5 nm,電荷密度為0.098 6 a.u.,此鍵略有減弱.MA2中Cu—O(1)鍵長(zhǎng)為0.208 2 nm,電荷密度為0.065 1 a.u.;成鍵略有減弱.MA1→TSA1→MA2過(guò)程中,CuI連接到O(1)上,H從O上經(jīng)過(guò)過(guò)渡態(tài)到Cu上,進(jìn)而和I結(jié)合并脫去HI.RE2中Br—C(1)鍵長(zhǎng)為 0.191 9 nm,電荷密度為0.158 5 a.u..生成中間體MA3后,Cu—O(1)鍵長(zhǎng)為0.191 0 nm,電荷密度為0.098 1 a.u.,略有增強(qiáng),Br—C(1)鍵長(zhǎng)為0.196 0 nm,電荷密度為0.145 2 a.u.,此鍵略有減弱,Br—Cu鍵長(zhǎng)為0.233 5 nm,電荷密度為0.069 0 a.u.,說(shuō)明已經(jīng)形成了中間體MA3.MA3到TSA2過(guò)程中,C(1)—Cu鍵有一定的成鍵,鍵長(zhǎng)為0.196 9 nm,電荷密度為0.103 1 a.u.,C(1)—Br鍵長(zhǎng)為0.261 7 nm,此一步表現(xiàn)了苯環(huán)從 Br上往 Cu上過(guò)渡的狀態(tài).在TSA2到MA4過(guò)程中,C(1)—Cu鍵長(zhǎng)0.196 9 nm減小到0.193 6 nm,電荷密度增加到0.117 1 a.u.,此鍵進(jìn)一步增強(qiáng)成鍵,同時(shí)Cu—N也有一定的成鍵,鍵長(zhǎng)為0.203 8 nm,電荷密度為0.082 0 a.u.,形成MA4.在MA3→TSA2→MA4過(guò)程中,苯環(huán)從Br上過(guò)渡到 Cu上,同時(shí) Cu—N成鍵.在 MA4到TSA3過(guò)程中,Br已經(jīng)從Cu上往H(1)上遷移.H(1)—Br成鍵鍵長(zhǎng)為0.148 8 nm,N—H(1)鍵減弱到鍵長(zhǎng)為0.181 2 nm,電荷密度為0.044 2 a.u.,同時(shí)N—Cu鍵長(zhǎng)為0.190 9 nm,電荷密度0.114 4 a.u..在TSA3到MA5過(guò)程中,HBr的進(jìn)一步離去,O—Cu鍵長(zhǎng)從0.184 6 nm增長(zhǎng)到0.186 1 nm,電荷密度從0.115 8 a.u.減小到0.111 9 a.u.,此鍵有所減弱,Cu—C鍵長(zhǎng)從0.191 6 nm減小到0.188 0 nm,電荷密度從0.123 2 a.u.增加到0.134 4 a.u.,此鍵略有加強(qiáng),Cu—N鍵長(zhǎng)從0.190 9 nm減小到0.186 8 nm,電荷密度從0.114 9 a.u.增加到0.124 7 a.u.,此鍵有所加強(qiáng);進(jìn)一步驗(yàn)證了HBr的離去,利于Cu—C(1)和Cu—N的成鍵.在MA4→TSA3→MA5過(guò)程中,Cu—N鍵進(jìn)一步增強(qiáng),HBr離去.在 MA5到 TSA4過(guò)程中,C(1)—Cu鍵長(zhǎng)從0.188 0 nm增長(zhǎng)到0.190 4 nm,電荷密度從0.134 4 a.u.減小到0.116 0 a.u.,此鍵減弱.N—C(1)在中間體 MA5中未成鍵,到 TSA4有成鍵,鍵長(zhǎng)為0.205 4 nm,電荷密度為0.081 1 a.u..表明了苯環(huán)在向N上靠近,有一定的成鍵.在TSA4到MA6過(guò)程中,N—C(1)鍵長(zhǎng)0.205 4 nm減少到0.145 7 nm,電荷密度從0.081 1 a.u.增加到0.301 4 a.u.,說(shuō)明N與苯環(huán)已經(jīng)成鍵.O—Cu鍵長(zhǎng)從0.190 3 nm增加到0.212 1 nm,電荷密度從0.100 0 a.u.減小到0.066 4 a.u.,表面 Cu原子有脫離的趨勢(shì),在MA5→TSA4→MA6過(guò)程中,苯環(huán)從Cu原子轉(zhuǎn)移到N原子上,并且Cu—N鍵斷開(kāi).

2.2 路徑IB反應(yīng)機(jī)理分析表2為路徑IB上各化合物能量以及校正后各物質(zhì)的相關(guān)能量.該反應(yīng)機(jī)理中所有反應(yīng)物、中間體、過(guò)渡態(tài)和產(chǎn)物的優(yōu)化分子構(gòu)型以及部分鍵臨界點(diǎn)(BCP)電荷密度(圖中圓括號(hào)內(nèi))如圖4所示.

表2 反應(yīng)路徑IB中各駐點(diǎn)能量E、相對(duì)能量Erel和過(guò)渡態(tài)的振動(dòng)頻率vTable 2 Energies E,relative energies Erel,and frequencies v of the compounds in the reactions for IB

圖3 路徑IA上化合物的構(gòu)型參數(shù)Fig.3 Geometric parameters of the compounds in pathway IA

圖4 路徑IB上化合物的構(gòu)型參數(shù)Fig.4 Geometric parameters of the compounds in pathway IB

在路徑IB中,反應(yīng)的第一步是催化劑CuI中的Cu原子與反應(yīng)物α-氨基酸的中的N(8)活化絡(luò)合,得到中間體 MB1.CuI中 Cu與 I鍵長(zhǎng)為0.230 1 nm,電荷密度為0.092 3 a.u.,形成中間體MB1時(shí)Cu與I鍵長(zhǎng)增長(zhǎng)到0.238 5 nm,電荷密度減少到0.089 8 a.u.,同時(shí)Cu與N鍵長(zhǎng)為0.198 3 nm,表面Cu與N之間有成鍵.NBO計(jì)算分析得知,在中間體MB1中,BD(1)N(8)—H(18)→LP?(6)Cu、BD(1)N(8)—H(19)→LP?(6)Cu間的二階穩(wěn)定化能分別為35.59和35.51 kJ/mol,表明它們之間存在軌道間相互作用.從表2能量數(shù)據(jù)可以看出,中間體MB1的能量比α-氨基酸與催化劑CuI能量之和低138.88 kJ/mol,表明中間體MB1易形成并穩(wěn)定存在,而且為后續(xù)反應(yīng)提供潛熱.中間體MB1上的H(1)原子經(jīng)過(guò)渡態(tài)TSB1遷移到I原子上并且脫去形成中間體MB2.在MB1到TSB1過(guò)程中,Cu—N鍵長(zhǎng)從0.198 3 nm減小到0.181 1 nm,電荷密度達(dá)到了0.140 8 a.u.,說(shuō)明了Cu與N進(jìn)一步成鍵,同時(shí)N上H(1)已經(jīng)轉(zhuǎn)移到Cu上,Cu—N鍵長(zhǎng)為0.145 1 nm,電荷密度為0.131 7 a.u.,表面了Cu—N配位鍵的形成.另外Cu—Br鍵長(zhǎng)從0.238 5 nm增長(zhǎng)到0.240 2 nm,電荷密度從0.089 8 a.u.減小到0.081 4 a.u.,表明Cu—Br鍵在進(jìn)一步減弱,隨著HI的離去,中間體MB2甲基上的H和Cu有輕微的成鍵,其電荷密度為0.997 7 a.u.,伴隨著Cu—N鍵的減弱,鍵長(zhǎng)增長(zhǎng)到0.184 1 nm.RE2的Br原子與MB2的Cu原子絡(luò)合成鍵形成MB3.形成MB3后Br—C(1)鍵長(zhǎng)伸長(zhǎng)到0.195 1 nm,電荷密度為0.148 7 a.u.,同時(shí) Br—Cu鍵長(zhǎng)為0.234 4 nm,電荷密度為0.066 4 a.u.,Cu—N鍵略有增強(qiáng),鍵長(zhǎng)為0.138 9 nm,電荷密度為0.129 9 a.u..在中間體MB3上的苯環(huán)經(jīng)過(guò)渡態(tài)TSB2從I原子上遷移到Cu原子上形成中間體MB4,在過(guò)渡態(tài)TSB2中Br—C(1)鍵增大到0.218 8 nm,電荷密度減少為0.090 7 a.u.,表明Br—C(1)鍵有所減弱,Br—Cu鍵長(zhǎng)從0.234 4 nm到0.233 5 nm,電荷密度從0.064 4 a.u.到0.069 8 a.u.,此鍵幾乎不變,Cu原子與C(1)有一定的成鍵,鍵長(zhǎng)為0.206 7 nm,電荷密度為0.083 1 a.u.,表明了苯環(huán)有從Br原子上往Cu原子上的遷移.在中間體MB4中Cu—C(1)鍵進(jìn)一步加強(qiáng),表明Cu—C(1)進(jìn)一步成鍵,苯環(huán)已經(jīng)連接到了Cu原子上.中間體MB4上的苯環(huán)經(jīng)過(guò)渡態(tài)TSB3從Cu原子上遷移到N原子上形成中間體MB5,隨后CuBr脫去形成產(chǎn)物.

圖5 反應(yīng)過(guò)程各駐點(diǎn)能級(jí)變化圖Fig.5 The diagram of relative energies along the channels of the reactions

2.3 反應(yīng)過(guò)程能級(jí)變化分析反應(yīng)過(guò)程各駐點(diǎn)能級(jí)變化如圖5所示.從圖5中可以看出,路徑IA和路徑IB經(jīng)過(guò)L-纈氨酸(N-芳基-α-氨基酸)和催化劑CuI活化絡(luò)合成穩(wěn)定的中間體MA1和MB1,該過(guò)程中體系分別放出熱量82.9和138.9 kJ/mol,為后續(xù)反應(yīng)提供一定的潛熱.路徑IA中催化反應(yīng)速率控制步驟是MA4→TSA3→MA5過(guò)程,反應(yīng)活化能為282.2 kJ/mol;路徑IB中催化速率控制步驟是MB1→TSB1→MB2過(guò)程,反應(yīng)活化能為219.2 kJ/mol,由以上比較結(jié)果可以看出,IB通道具有較低的活化能,即IB通道為整個(gè)反應(yīng)的最優(yōu)反應(yīng)通道,通道路徑為RE1→MB1→TSB1→MB2→MB3→TSB2→MB4→TSB3→MB5→MB6.

3 結(jié)論

采用密度泛函理論在B3LYP/6-311+G?基組水平上對(duì)CuI催化L-纈氨酸(α-氨基酸)與溴苯反應(yīng)生成N-芳基-α-氨基酸的反應(yīng)機(jī)理進(jìn)行了理論研究,研究發(fā)現(xiàn)了2條可行的反應(yīng)通道,路徑IA:RE1→MA1→TSA1→MA2→MA3→TSA2→MA4→TSA3→MA5→TSA4→MA6;速率控制步驟為MA4→TSA3→MA5過(guò)程,反應(yīng)活化能為282.2 kJ/mol;路徑 IB:RE1→MB1→TSB1→MB2→MB3→TSB2→MB4→TSB3→MB5→MB6,速率控制步驟為MB1→TSB1→MB2過(guò)程,活化能為219.2 kJ/mol.通過(guò)比較各反應(yīng)通道中控制步驟的活化能,得出IB通道具有較低的活化能,即IB通道為整個(gè)反應(yīng)的最優(yōu)反應(yīng)通道,該通道的第一步活化絡(luò)合過(guò)程能量降低了138.9 kJ/mol,為其速率控制步驟MB1→TSB1→MB2提供了一定的潛熱,使反應(yīng)易于進(jìn)行,生成N-芳基-α-氨基酸.同時(shí)還計(jì)算了用CuCl催化的路徑IB的控制步驟活化能為229.4 kJ/mol,高于CuI催化的活化能(219.2 kJ/mol),反應(yīng)速率降低,這與文獻(xiàn)報(bào)道是一致的.

[1]王曄峰,曾京輝,崔曉瑞.銅催化C—N交叉偶聯(lián)反應(yīng)的研究進(jìn)展[J].有機(jī)化學(xué),2010,30(2):181-199.

[2]Murahashi S I,Davies S G.Transition metal catalyzed reactions[J].Blackwell Science,1999,5:99-131.

[3]Zhou Y,Zhao J,Liu L.Meta-selective transition-metal catalyzed arene C—H Bond functionalization[J].Angew Chem Int Ed Engl,2009,48:7126-7128.

[4]Armin D M,Diederich F.Metal-catalyzed Cross-coupling Reactions[M].2nd Ed.Weinheim:Wiley-VCH,2004.

[5]Beletskaya I P,Cheprakov A V.Copper in cross-coupling reactions[J].The Post-Ullmann Chemistry Coord Chem Rev,2004,248:2337-2364.

[6]Masanari K,Daisuke N,Masahiro F,et al.Rh-catalyzed reductive coupling reaction of aldehydes with conjugated dienes promoted by triethylborane[J].Org Lett,2009,17:3794.

[7]徐伯華.2-碘代硒苯與2-唑烷酮在CuI催化下的C—N偶聯(lián)反應(yīng)機(jī)理的理論研究[J].原子與分子物理學(xué)報(bào),2011,28(6):995.

[8]于海珠,傅堯,白小宇,等.亞銅催化的碳雜偶聯(lián)反應(yīng)機(jī)理[J].化學(xué)進(jìn)展,2010,22(4):557.

[9]李來(lái)才,蔡皖飛,汪曉慧,等.N-(鄰氯苯基)苯甲酰胺在CuX(X=I,Br)催化下分子內(nèi)O-芳基化反應(yīng)機(jī)理[J].物理化學(xué)學(xué)報(bào),2009,25:2101.

[10]Ullmann F.The Ullmann condensation and the synthesis of diarylamines[J].Ber Dtsch Chem Ges,1903,36:2382-2384.

[11]Kohn W,Sham L J.Self-Consistent equations including exchange and correlation effects[J].Phys Rev,1965,A140:1133.

[12]Hariharan P C,Pople J A.The influence of polarization functions on molecular orbital hydrogenation energies[J].Theory Chem Acta,1973,28:213.

[13]Reed A E,Weinhold F,Curtiss L A,et al.Natural bond orbital a nalysis of molecular interactions:theoretical studies of binary complexes of HF,H2O,NH3,N2,O2,F2,CO and CO2with HF,H2O and NH3[J].Chem Phys,1986,84:5687-5706.

[14]Biegler K F,Sch?nbohm J,Derdau R,et al.AIM 2000[S].Ver 2.0.Hamilton:McMaster University,2002.

[15]Frisch M J,Trucks G W,Schlegel H B,et al.Gaussian 09[S].Rev A.02.Pitsburgh:Gaussian,2009.

主站蜘蛛池模板: 亚洲精品国产成人7777| 亚洲人成在线免费观看| 99热这里只有精品国产99| 久久永久视频| 亚洲清纯自偷自拍另类专区| 国产一级片网址| 免费女人18毛片a级毛片视频| 国产一区二区丝袜高跟鞋| 2021最新国产精品网站| 午夜精品久久久久久久无码软件 | 久热99这里只有精品视频6| 亚洲va视频| 性色在线视频精品| 国产精品大尺度尺度视频| 亚洲婷婷六月| 三上悠亚精品二区在线观看| 精品一区二区三区视频免费观看| 91在线中文| 欧美视频在线不卡| 一本久道热中字伊人| 久草青青在线视频| AV无码无在线观看免费| 女人av社区男人的天堂| 精品超清无码视频在线观看| 中文字幕在线观看日本| 中文字幕中文字字幕码一二区| 青青草91视频| 亚洲国产亚综合在线区| 欧美成在线视频| 国产精品蜜臀| 奇米影视狠狠精品7777| 一区二区日韩国产精久久| 国产精品免费入口视频| 国产欧美网站| 亚洲黄色视频在线观看一区| 国产色偷丝袜婷婷无码麻豆制服| 伊人91在线| 国产成人久久综合777777麻豆 | 中文字幕在线观| 2020国产在线视精品在| 欧美成a人片在线观看| 伊人久久精品无码麻豆精品| 欧美亚洲国产精品第一页| 色哟哟精品无码网站在线播放视频| 国产极品嫩模在线观看91| 91成人在线免费观看| 四虎永久在线精品影院| 91久久偷偷做嫩草影院免费看| 毛片基地美国正在播放亚洲| 四虎永久在线| 国产人前露出系列视频| 久久久受www免费人成| 日本尹人综合香蕉在线观看| 国产午夜人做人免费视频中文| 国产另类乱子伦精品免费女| 国产中文一区a级毛片视频 | 婷婷综合亚洲| 国产成人久视频免费| 久久中文无码精品| 秋霞国产在线| 欧美三级日韩三级| 亚洲中文字幕无码mv| 亚洲成人77777| 国产成人综合网在线观看| 无码福利视频| 国产真实自在自线免费精品| 欧美亚洲中文精品三区| 中文成人在线| 亚洲国产欧美自拍| 国产精品欧美日本韩免费一区二区三区不卡 | 日韩福利视频导航| 欧美高清日韩| 亚洲激情99| 91亚洲国产视频| 久久久久久午夜精品| 国产屁屁影院| 欧美精品成人| 天堂成人av| 亚洲综合国产一区二区三区| 国产尤物视频在线| 伊伊人成亚洲综合人网7777| 色老二精品视频在线观看|