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

四(鄰氟芐基)錫及一維鏈狀三芐基錫二茂鐵甲酸酯的合成、結構及性質

2015-12-05 07:28:02張復興鄺代治馮泳蘭王劍秋庾江喜蔣伍玖朱小明
無機化學學報 2015年6期
關鍵詞:結構

張復興 鄺代治 馮泳蘭 王劍秋 庾江喜 蔣伍玖 朱小明

(衡陽師范學院化學與材料科學系,功能金屬有機材料湖南省普通高等學校重點實驗室,衡陽 421008)

四(鄰氟芐基)錫及一維鏈狀三芐基錫二茂鐵甲酸酯的合成、結構及性質

張復興*鄺代治馮泳蘭王劍秋庾江喜蔣伍玖朱小明

(衡陽師范學院化學與材料科學系,功能金屬有機材料湖南省普通高等學校重點實驗室,衡陽421008)

合成了四(鄰氟芐基)錫(1)和一維鏈狀三芐基錫二茂鐵甲酸酯(2),經X射線衍射方法測定了化合物的晶體結構。化合物1屬四方晶系,空間群I41/a,晶體學參數a=1.967 94(16)nm,b=1.967 94(16)nm,c=0.593 16(5)nm,V=2.297 2(3)nm3,Z=4,Dc= 1.605 g·cm-3,μ(Mo Kα)=11.58 cm-1,F(000)=1 112,R1=0.020 8,wR2=0.057 6。化合物2屬單斜晶系,空間群為P21/n,晶體學參數:a=1.593 54(12)nm,b=1.007 23(8)nm,c=1.700 26(13)nm,β=91.001(10)°,V=2.728 6(4)nm3,Z=4,Dc=1.521 g·cm-3,μ(Mo Kα)=14.74 cm-1,F(000)=1 256,R1=0.038 4,wR2=0.095 2。晶體的中心錫原子分別呈四配位畸變四面體構型和五配位畸變三角雙錐構型。對其結構進行量子化學從頭計算,探討了配合物的穩定性、分子軌道能量以及部分前沿分子軌道的組成特征。測定了配合物的熱穩定性和體外抗癌活性。

四(鄰氟芐基)錫;三芐基錫二茂鐵甲酸酯;合成;結構;抗癌活性

烴基錫化合物由于具有結構的多樣性、豐富的反應性、較強的生物活性和催化活性,多年來一直引起人們的興趣[1-6]。烴基錫化合物中有四烴基錫、烴基鹵化錫、烴基氧化錫、烴基錫羧酸酯配合物等多種類型。近年來人們合成了一系列具有結構特點和特殊性能的有機錫化合物[7-14],芐基錫化合物特別是取代芐基錫化合物的合成,更豐富了烴基錫化合物的內容。為了更進一步探索取代芐基錫類化合物的結構與性能的關系,本文合成了四(鄰氟芐基)錫(1)和一維鏈狀三芐基錫二茂鐵甲酸酯(2),并通過元素分析、紅外光譜對其結構進行了表征。用X-射線單晶衍射測定了該化合物的晶體結構,對其結構進行量子化學從頭計算,探討了配合物分子的穩定性、分子軌道能量以及一些前沿分子軌道的組成特征。測定了配合物的熱穩定性和體外抗癌活性。

1 實驗部分

1.1試劑和儀器

日本島津IRPrestige-21紅外光譜儀(4 000~400 cm-1,KBr),PE-2400元素分析儀,Bruker SMART APEXⅡ單晶衍射儀,spectramax M5連續光譜酶標測試儀,TG209F3熱分析儀,X4數字顯微熔點,溫度計未經校正。所有試劑均為化學純。

1.2化合物的合成

化合物1:在三頸瓶中加入14.5 g(0.1 mol)鄰氟芐基氯,少量碘,60 mL正戊醇和40 mL鄰二甲苯,回流攪拌下分批加入5.93 g(0.05 mol)經活化處理的錫粉和1.0 g鎂條,快速攪拌回流反應6 h。趁熱過濾出未反應的錫粉,濾液倒入100 mL 5%HCl溶液中,充分攪拌。冷卻后分出水層,旋轉蒸發部分溶劑后,放置析出白色固體。抽濾分離出固體,用少量冷無水乙醇洗滌固體,再用苯重結晶,得到無色晶體10.84 g,收率78.15%。熔點:92~94℃。紅外光譜主要吸收峰:3 053(w),2 961(w),1 578(m),1 491 (s),14 541(m),519(m)cm-1。元素分析(C28H24F4Sn),計算值(%):C,60.57;H,4.33。實測值(%):C,60.22;H,4.24。

化合物2:在50 mL圓底燒瓶中,加入0.855 g(2 mmol)三芐基氯化錫、0.460 g(2 mmol)二茂鐵甲酸、40 mL苯和2.5 mmol三乙胺,在電磁攪拌下加熱回流反應5 h。趁熱過濾除去不溶性固體,濾液旋轉蒸發除去部分溶劑,放置析出棕色固體,用苯重結晶得棕紅色晶體0.816 3 g,產率65.73%。熔點:257~259℃。紅外光譜主要吸收峰:3 019(m),29 229(w),15 979(m),1 557(s),1 498(s),1 462(m),1 385(m),5 695(w),515(w)cm-1。元素分析(C32H30O2FeSn),計算值(%):C,61.89;H,4.87。實測值(%):C,61.78;H,4.82。

1.3晶體結構分析

分別選取大小為0.24mm×0.19 mm×0.17 mm(1)和0.29 mm×0.25 mm×0.23 mm(2)的晶體,在Bruker SMART APEXⅡCCD單晶衍射儀上,采用經石墨單色化的Mo Kα射線(λ=0.071 073 nm),于296(2) K,以φ~ω掃描方式收集數據。可觀察衍射點分別為1 256個和5 185個(I>2σ(I))用于結構分析和精修。全部數據經Lp因子和經驗吸收校正。晶體結構由直接法解出,非氫原子坐標通過數輪差值Fourier合成陸續確定,理論加氫法給出氫原子在晶胞中的位置坐標。對非氫原子坐標及其各向異性熱參數進行全矩陣最小二乘法修正。全部結構分析計算工作采用SHELX-97[15]程序系統完成。晶體學數據詳見表1。

CCDC:1056036,1;1056037,2。

表1 化合物的晶體學數據Table1 Crystallographic data of the title complexes

續表1

1.4熱穩定性和抗癌活性測定

利用TG209F3熱分析儀,在空氣氛中,升溫速率為20℃·min-1的實驗條件下測定配合物的熱穩定性。

采用MTT法[16]檢測化合物對人宮頸癌細胞(HeLa)的體外抗性腫瘤活性。將化合物用二甲基亞砜(DMSO)配成5.0 mg·mL-1的溶液,用RPMI-1640培養基分別稀釋成5,10,25,50,100 μg·mL-1。取處于指數生長期的HeLa細胞懸液加至96孔板中(細胞濃度為50 000個·mL-1,100 μL每孔),于37℃、5%CO2恒溫箱中培養16 h使細胞貼壁。去除上清,加100 μL不同濃度的上述化合物,每個濃度設 4個復孔,孵育24 h,棄去上清,每孔加入2.0 mg·mL-1的MTT溶液60 μL,繼續培養3h,除上清液后,每孔加入150 μL二甲基亞砜,低速振蕩10 min,使深藍色結晶溶解,用酶標儀在490 nm波長處測其吸光度值。按如下公式:抑制率r=(OD對照組-OD測試組)/ OD對照組×100%,計算各組對癌細胞的抑制率。每組實驗均重復3次,取其平均值。

2 結果與討論

2.1晶體結構描述

化合物的主要鍵長和鍵角分別列于表2、3,化合物的分子結構見圖1、圖2。

表2 化合物1的主要鍵長(nm)和鍵角(°)Table2 Selected distances(nm)and angles(°)of complex 1

表3 化合物2的主要鍵長(nm)和鍵角(°)Table3 Selected distances(nm)and angles(°)of complex 2

化合物1:由分子結構圖和結構參數可知,化合物中心錫原子與4個亞甲基碳原子相連形成四面體構型。4個Sn-C鍵鍵長相等,均為0.216 52(19) nm。錫原子和與之相連的碳原子之間所構成的鍵角中∠C(1i)-Sn(1)-C(1iii)與∠C(1ii)-Sn(1)-C(1)相等,為110.78(11)°,其余的4個鍵角相等,為108.82(6)°。分子中0.216 52 nm的Sn-C鍵鍵長,比四(對氟芐基)錫的Sn-C鍵平均鍵長(0.216 2 nm)[17]稍長,這是因為氟原子在鄰位,使得中心錫原子周圍的空間擁擠程度有所增加;而比三(鄰氟芐基)氯化錫的Sn-C鍵平均鍵長(0.214 3 nm)[17]要長得多,是因為氯原子被體積大很多的鄰氟芐基取代后,使得中心錫原子周圍的空間擁擠程度大大增加。

圖2 化合物2的分子結構圖(橢球概率30%)Fig.2 Molecular structure of complex 2 with the ellipsoids drawn at the 30%probability level

化合物2:從分子結構圖和結構參數可知,三芐基錫通過羧基橋聯形成一維鏈狀結構。中心錫原子分別與3個芐基的亞甲基碳、不同的二茂鐵甲酸分子的2個羧基氧,形成五配位的三角雙錐結構。3個芐基碳原子C(1)、C(8)、C(15)占據了三角雙錐赤道平面的3個位置,2個氧原子O(1)、O(2i)則占據了赤道平面兩側的軸向位置。處于軸向位置的氧原子O(1)、O(2i)與處于赤道位置的碳原子的鍵角分別為:C(1)-Sn(1)-O(1)92.20(12)°、O(1)-Sn(1)-C(8)93.65(13)°、O(1) -Sn(1)-C(15)98.74(14)°、C(1)-Sn(1)-O(2i)79.74(12)°、C(8)-Sn(1)-O(2i)84.79(12)°、C(15)-Sn(1)-O(2i)90.42(13)°,都與90°有偏差,其偏差范圍在8.74°~10.26°。處于軸向位置的原子的鍵角O(1)-Sn(1)-O(2i)169.97(9)°,比180°線性角小了10.03°。中心錫原子與處于赤道平面的3個碳原子間的夾角分別為:C(1)-Sn(1)-C(15) 118.82(17)°、C(1)-Sn(1)-C(8)116.66(16)°、C(15)-Sn(1)-C(8)122.33(15)°,3個角之和為357.81°,與360°有較大的偏差。由此可知化合物為畸變程度較大的三角雙錐結構。

化合物2的結構中2個羧基氧原子以雙齒配位分別與2個不同的錫原子成鍵,其鍵長分別為:Sn(1)-O(1)0.214 6 nm和Sn(1i)-O(2)0.243 0 nm,由此而形成了一維鏈狀的分子結構,如圖3所示。

圖3 化合物2的一維鏈狀結構Fig.3 One-dimensional chain structure of complex 2

2.2量子化學研究

2.2.1分子的總能量和前沿分子軌道能量

根據晶體結構的原子坐標,運用Gaussian 03W程序和B3ylp/lanl2dz基組水平,計算得到分子的總能量和前沿分子軌道能量。

化合物1:ET=-373.141 486 9 a.u.,EHOMO=-0.384 88 a.u.,ELUMO=-0.314 06a.u.,ΔELUMO-HOMO=0.070 82 a.u.。從體系能量和前沿軌道的能量分析,化合物總能量較高,前沿占有軌道能量雖然低,但前沿未占軌道能量也較低,最高占據軌道與最低未占軌道的能量間隙ΔE很小,僅為0.070 82 a.u.,表明化合物分子結構穩定性有限,從氧化還原或電荷轉移的角度分析,化合物較容易失去電子而被氧化。

化合物2:ET=-1 513.896 132 3 a.u.,EHOMO=-0.198 50 a.u,ELUMO=-0.054 05 a.u.,ΔELUMO-HOMO=0.144 5 a.u.。從體系能量和前沿軌道的能量分析,化合物總能量低,前沿占有軌道能量較低,表明化合物分子結構有一定的穩定性。但前沿未占軌道能量也較低,最高占據軌道與最低未占軌道的能量間隙ΔE較小,只有0.144 5 a.u,因此,從氧化還原或電荷轉移的角度分析,化合物也較容易失去電子而被氧化。

2.2.2軌道成分分析

為探索化合物的電子結構與成鍵特征,對化合物分子軌道進行分析,用參與組合的各類原子軌道系數的平方和來表示該部分在分子軌道中的貢獻,并經歸一化。

化合物1:把化合物原子分為五部分:(a)錫原子Sn;(b)亞甲基碳原子C;(c)苯環碳原子C;(d)氟原子F;(e)氫原子H,前沿占有軌道和未占有軌道各取5個,計算結果見表4和圖4。

圖4 化合物1的前沿分子軌道示意圖Fig.4 Schematic diagram of frontier MO for complex 1

表4 化合物1的分子軌道組成Table4 Some calculated frontier molecular orbitals composition of complex 1

表4和圖4顯示化合物1分子的成鍵特征:①前沿占有分子軌道中,芐基苯環碳原子,對分子軌道的貢獻最大,達到了79.33%,并且在深層次軌道中均有較大的貢獻,說明苯環具有良好的共軛離域性和穩定性。②前沿占有分子軌道中,錫原子和芐基亞甲基碳原子對分子軌道貢獻都較小,分別為6.83%和6.26%,并且在深層次軌道中的貢獻也都較小,說明Sn-C鍵強度弱,分子的穩定有一定的限度。從鍵長來看,化合物1的Sn-C鍵比三(鄰氟芐基)氯化錫和四(對氟芐基)錫[17]的Sn-C鍵平均鍵長長,這與分子的穩定性相吻合。③比較HOMO與LUMO的各類原子軌道成份,可以看出,當電子從HOMO激發到LUMO軌道時,主要是鄰氟芐基上的電子向錫原子轉移。

化合物2:把化合物原子分為五部分:(a)二茂鐵基碳原子和鐵原子A;(b)羧基碳原子和氧原子B;(c)芐基碳原子C;(d)錫原子Sn;(e)氫原子H,前沿占有軌道和未占有軌道各取5個,計算結果見表5和圖5。

圖5 化合物2的前沿分子軌道示意圖Fig.5 Schematic diagram of frontier MO for complex 2

表5和圖5顯示化合物2分子的成鍵特征:①前沿占有分子軌道中,對分子軌道的貢獻幾乎來自二茂鐵基,達到了99.24%,并且在深層次軌道中仍有較大的貢獻,表明二茂鐵基具有很好的穩定性。②前沿占有分子軌道中,與錫原子直接相連的基團對分子軌道貢獻都很小,分別為羧基0.554 356%、芐基0.175 271%;而錫原子對分子軌道貢獻則更小,只有0.011 323%;表明它們之間鍵的強度較弱,分子穩定性有限,這與從能量的角度分析的結果是一致的。③比較HOMO與LUMO的各類原子軌道成份,可以看出,當電子從HOMO激發到LUMO軌道時,主要是二茂鐵基上的電子通過羧基向錫原子轉移,羧基既是電子轉移的橋梁,也是電子轉移的受體。

表5 化合物2的分子軌道組成Table 5 Some calculated frontier molecular orbitals composition of complex 2

2.3熱穩定性分析

化合物的熱重曲線如圖6、圖7所示。

圖6 化合物1的TG曲線Fig.6 TG curve of the complex 1

圖6顯示,隨溫度的升高,在初始階段40~187℃,化合物1基本沒有失重;在187~700℃區間,化合物出現明顯失重,總計72.75%,對應于配合物分子失去4個鄰氟芐基;最后穩定在約27.25%,殘余物可被視為SnO2,與27.14%的計算值吻合;上述熱分析結果表明化合物1在187℃之前是可以穩定存在的。

圖7 化合物2的TG曲線Fig.7 TG curve of the complex 2

圖7顯示,隨溫度的升高,在初始階段40~70℃,化合物2幾乎沒有任何失重;在70~530℃區間,化合物出現明顯失重,總計63.65%,對應于配合物分子發生熱分解失去有機基團;最后穩定在約36.34%,殘余物可被假定為SnO2+Fe3O4,與36.69%的計算值吻合;上述熱分析結果表明化合物2在70℃之前可以穩定存在,熱穩定性比化合物1低。

2.4化合物的體外抗癌活性

不同濃度梯度下化合物對肝癌細胞的抑制率數值見表6。

由表可以看出2個化合物對HeLa細胞株增殖的抑制情況,化合物1在高濃度時對HeLa人宮頸癌細胞有一定的抑制能力,但隨濃度降低其抑制率迅速變小,當濃度為5 μg·mL-1時,抑制率只有27.50%。而化合物2無論是在高濃度還是低濃度,都幾乎沒有抑制能力。

表6 不同濃度化合物對HeLa人宮頸癌細胞的抑制率Table 6 Growth inhibition rates for HeLa cell lines with title compounds at different concentrations

[1]Chandrasekhar V,Thirumoorthi R,Metre R K,et al.J. Organomet.Chem.,2011,696:600-606

[2]Effendy,Marchetti F,Marinelli A,et al.Inorg.Chim.Acta, 2011,366(1):388-393

[3]Munoz-Flores B M,Santillan R,Farfan N,et al.J.Organomet. Chem.,2014,769:64-71

[4]Hanif M,Hussain M,Ali S,et al.Polyhedron,2010,29:613-619

[5]ZHANG Xiao-Yan(張曉燕),YANG Guang(楊光),ZHANG Jun(張俊),et al.Chem.J.Chinsese Universities(高等學校化學學報),2010,31(6):1162-1166

[6]Siddiqi Z A,Shahid M,Kumar S,et al.J.Organomet.Chem., 2009,694:3768-3774

[7]Mahmudov K T,Guedes da Silva M F C,Kopylovich M N, et al.J.Organomet.Chem.,2014,760:67-73

[8]Ruan B F,Tian Y U,Zhou H P,et al.Inorg.Chim.Acta, 2011,365:302-308

[9]YINHan-Dong(尹漢東),WANG Chuan-Hua(王傳華),WANG Yong(王勇),et al.Chinese J.Inorg.Chem.(無機化學學報), 2002,18(2):201-204

[10]YAN Wen-Hua(閆文華),KANG Wan-Li(康萬利),LI Jin-Huan(李金環).Chinese J.Appl.Chem.(應用化學),2007,24 (6):660-664

[11]Shujha S,Shah A,Rehman Z U,et al.Eur.J.Med.Chem., 2010,45:2902-2911

[12]ZHANG Fu-Xing(張復興),KUANG Dai-Zhi(鄺代治),WANG Jian-Qiu(王劍秋),et al.Chinese J.Org.Chem.(有機化學), 2008,28(8):1457-1461

[13]KUANG Dai-Zhi(鄺代治),JIANG Wu-Jiu(蔣伍玖),FENG Yong-Lan(馮泳蘭),et al.Chinese J.Inorg.Chem.(無機化學學報),2012,27(10):1981-1986

[14]WANGJian-Qiu(王劍秋),ZHANGFu-Xing(張復興),KUANG Dai-Zhi(鄺代治),et al.Chinese J.Inorg.Chem.(無機化學學報),2015,27(31):237-2426

[15]Sheldrick G M.SHELXL-97,Program for Crystal Structure Refinement,University of G?ttingen,Germany,1997.

[16]PENG Jun-Mei(彭俊梅),LI Wan(李婉),SHEN Kun(申坤), et al.Chem.J.Chinese Universities(高等學校化學學報), 2013,34(7):1646-1652

[17]ZHANGFu-Xing(張復興),WANGJian-Qiu(王劍秋),KUANG Dai-Zhi(鄺代治),et al.Chinese J.Org.Chem.(有機化學), 2003,23(4):368-371

Syntheses,Crystal Structures and Properties of the Tetra(o-fluorobenzyl)tin and the Tribenzyltin Ferrocenecarboxylate

ZHANG Fu-Xing*KUANG Dai-ZhiFENG Yong-LanWANG Jian-Qiu YU Jiang-XiJIANG Wu-JiuZHU Xiao-Ming
(Department of Chemistry and Material Science,Hengyang Normal University,Key Laboratory of Functional Organometallic Materials of Hengyang Normal University,College of Hunan Province,Hengyang,Hunan 421008,China)

The tetra(o-fluorobenzyl)tin(1)and the tribenzyltin ferrocenecarboxylate(2)have been synthesized.The crystal structures of the complexes were determined by X-ray diffraction.The crystal 1 belongs to tetragonal space group I41/a with a=1.967 94(16)nm,b=1.967 94(16)nm,c=0.593 16(5)nm,V=2.297 2(3)nm3,Z=4,Dc= 1.605 g·cm-3,μ(Mo Kα)=11.58 cm-1,F(000)=1 112,R1=0.020 8,wR2=0.057 6.The crystal 2 belongs to monoclinic space group P21/n with a=1.593 54(12)nm,b=1.007 23(8)nm,c=1.700 26(13)nm,β=91.001(10)°,V=2.728 6(4) nm3,Z=4,Dc=1.521 g·cm-3,μ(Mo Kα)=14.74 cm-1,F(000)=1 256,R1=0.038 4,wR2=0.095 2.The stabilities,some frontier molecular orbital energies and composition characteristics of some frontier molecular orbital of the complexes have been investigated by quantum chemistry calculation.Further more,thermal stability and anticancer activity of the complexes were tested.CCDC:1056036,1;1056037,2.

tetra(o-fluorobenzyl)tin;tribenzyltin ferrocenecarboxylate;synthesis;crystal structure

O614.43+2

A

1001-4861(2015)06-1194-07

10.11862/CJIC.2015.157

2015-03-11。收修改稿日期:2015-03-27。

湖南省自然科學基金(No.13JJ3112)、湖南省科技計劃(No.2013TZ2025)、湖南省高校創新平臺開放基金(No.13K105)、湖南省重點學科基金、湖南省高校科技創新團隊支持計劃和湖南省高校重點實驗室開放基金資助。

*通訊聯系人。E-mail:zfx8056@163.com;會員登記號:S060018907M。

猜你喜歡
結構
DNA結構的發現
《形而上學》△卷的結構和位置
哲學評論(2021年2期)2021-08-22 01:53:34
論結構
中華詩詞(2019年7期)2019-11-25 01:43:04
新型平衡塊結構的應用
模具制造(2019年3期)2019-06-06 02:10:54
循環結構謹防“死循環”
論《日出》的結構
縱向結構
縱向結構
我國社會結構的重建
人間(2015年21期)2015-03-11 15:23:21
創新治理結構促進中小企業持續成長
現代企業(2015年9期)2015-02-28 18:56:50
主站蜘蛛池模板: 亚洲无码高清一区| 欧美一区二区精品久久久| 日韩黄色大片免费看| 中文无码精品A∨在线观看不卡| 亚洲天堂在线免费| 国产精品综合久久久| 欧美一区二区福利视频| 亚洲成年人网| 日韩二区三区无| 亚洲视频无码| 久久一色本道亚洲| 在线日本国产成人免费的| 九色国产在线| 中国美女**毛片录像在线| 一区二区三区国产| 手机永久AV在线播放| 色综合热无码热国产| 国产精品久久国产精麻豆99网站| 国产一级视频久久| 国产免费久久精品99re丫丫一| 日韩精品一区二区三区视频免费看| 久久一本精品久久久ー99| 欧美日韩国产综合视频在线观看| 国产亚洲高清视频| 中文字幕自拍偷拍| 国产在线高清一级毛片| 五月天婷婷网亚洲综合在线| 亚洲精品国产综合99久久夜夜嗨| 91精品综合| 精品一区二区三区中文字幕| 免费观看欧美性一级| 99视频在线免费看| 免费看黄片一区二区三区| 在线观看精品自拍视频| 久青草网站| 久久一级电影| 国产激爽爽爽大片在线观看| 亚洲激情区| 国产人成在线视频| 在线观看免费AV网| 69视频国产| 91精品aⅴ无码中文字字幕蜜桃| 污视频日本| 亚洲第一视频区| 91年精品国产福利线观看久久| 亚洲中文在线看视频一区| 久久午夜夜伦鲁鲁片不卡| 国产成人综合在线观看| 亚洲成人网在线播放| 欧美成人手机在线视频| 精品91自产拍在线| 国产精品99在线观看| 日韩成人免费网站| 亚洲精品国产综合99久久夜夜嗨| 国产网站在线看| 青青草91视频| 国产网友愉拍精品视频| 999福利激情视频| 97国产在线观看| 亚洲成a人在线播放www| 亚洲欧美日韩动漫| 国产一级毛片yw| 色天天综合久久久久综合片| 亚洲精品中文字幕午夜| 久久久久人妻一区精品| 女人爽到高潮免费视频大全| 色综合手机在线| 国产成人免费观看在线视频| 朝桐光一区二区| 亚洲日韩精品综合在线一区二区| 国产精品伦视频观看免费| 日韩免费成人| 岛国精品一区免费视频在线观看| 99热免费在线| 久久亚洲精少妇毛片午夜无码| 国产v精品成人免费视频71pao| 久久这里只有精品23| 久草国产在线观看| 成人av手机在线观看| 欧美激情视频一区二区三区免费| 精品久久综合1区2区3区激情| 国产色伊人|