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

    Self-assembly behavior of disaccharide-containing supra-amphiphiles

    2023-03-14 06:52:32ZhenfeiGoTinnnWngZikunRoHuiYnRnZhngGuosongChen
    Chinese Chemical Letters 2023年1期

    Zhenfei Go,Tinnn Wng,Zikun Ro,Hui Yn,Rn Zhng,Guosong Chen,*

    a The State Key Laboratory of Molecular Engineering of Polymers and Department of Macromolecular Science,Fudan University,Shanghai 200433,China

    b Changchun Institute of Applied Chemistry,Chinese Academy of Sciences,Changchun 130022,China

    c School of Pharmaceutical Sciences,Liaocheng University,Liaocheng 252059,China

    Keywords:Supra-amphiphile Disaccharide Self-assembly Dynamic covalent bond Molecular dynamics simulation

    ABSTRACT In this paper,supra-amphiphilic compounds containing disaccharides and azobenzene ends have been constructed via dynamic covalent bond.It was found that the slight structural difference of the disaccharides made significant difference in the self-assembled morphologies.Namely,three kinds of azodisaccharide supra-amphiphiles were found to assemble into different morphologies,with the only difference in chemical structure from the disaccharides.More importantly,the structural difference between the disaccharides,including lactoside,maltoside and cellobioside was trivial.Molecular simulation revealed the packing of molecules was due to the different contribution from hydrogen bonds.The above results clearly indicated the contribution of saccharide packing,especially the related hydrogen bonding,to the final morphology of the assembled structures.

    Supra-amphiphiles,distinguished from amphiphiles,are constructed by dynamic covalent bonds or non-covalent interactions[1,2],such as hydrogen bond,π-πinteraction,host-guest conjugation and charge transfer interaction [3–7].Supra-amphiphiles can be easily constructed through controlled manner,but compared with ordinary amphiphiles,they are much more tending to self-assemble even under very low concentrations due to the participation of dynamic interactions [8].Therefore,developing the supra-amphiphiles can enrich the preparation of complex selfassemblies [9].Meanwhile,in order to obtain the microscopic properties of the assemblies,many modern experimental methods including TEM,SEM and AFM have been widely used [1–4,10].In recent years,due to the substantial increase in computing power,there have been many reports on the use of Monte Carlo (MC)and molecular dynamics (MD) simulation methods to study microscopic properties [11–14].These kinds of researches provide indepth information about dynamics and structural properties at the microscopic level.

    It is known that the chemical structures of saccharides are complex and the variation between different saccharides could be trivial.However,such trivial difference could lead to significant property difference when the saccharides were polymerized and even bulk materials were formed.A most typical example is between starch and cellulose.Both of them are glucose-based macromolecules.The basic repeating unit of starch is maltose (α1–4-linked) and that of cellulose is cellobiose (β1–4-linked).Therefore,the structural difference between starch and cellulose only lies in different glycosidic bonds between glucoside units (αandβglycosidic bonds),which lead to completely different physiochemical properties.Moreover,the structural isomerism of disaccharides will affect their specific binding with proteins [15–18].However,in most of the self-assembly studies,the contribution of saccharide structures has been more or less overlooked.Thus,it is demanding to construct saccharides-containing amphiphilic molecules and study their self-assembly behaviors in a relatively systematic manner.For this purpose,supra-amphiphile could be a choice since different saccharides could be installed to the same hydrophobic moietiesviadynamic covalent bonding.

    Herein,for the first time,the disaccharide-containing supraamphiphiles were constructed by arylhydrazone dynamic covalent bond to explore the possibility of controlling self-assembly behaviors by varying saccharide structures.Azobenzene,a common building block in supramolecular chemistry,was chosen as the model hydrophobic part.The three constructed supra-amphiphiles containing different disaccharides,including lactoside (Lac),cellobioside (Cel) and maltoside (Mal) were assembled using the same assembly method.We found that slight difference in chemical structure (Fig.1) could induce drastic change in assembly morphology.Vesicles were observed for Azo-Lac.For Azo-Cel,the morphology was like a sunflower.For Azo-Mal,only some irregular assemblies have been observed.To understand the above results better,molecular dynamics simulations on three azo-disaccharide supra-amphiphiles were performed.It was found that the aggregation degree of the three disaccharides were shown in the sequence Mal>Cel>Lac due to their slight difference in the configurations,which can be employed to explain the assembly morphologies.

    Fig.1.The difference of chemical structure between three disaccharides: lactose,cellobiose and maltose.Light blue ring shows Galactose.Deep blue ring shows glucose.The α glycosidic bond shows as red,the β glycosidic bond shows as pink.

    Three disaccharides with slight differences in configurations,which are lactoside (galactoside-(β1→4)-glucoside),maltoside(glucoside-(α1→4)-glucoside) and cellobioside (glucoside-(β1→4)-glucoside),were first employed to explore structural influence on assembled morphologies.As shown in Fig.1,these three disaccharides are stereo-isomers.Both of Mal and Cel are based on glucosides,the only difference lies in the configuration of glycosidic bond,Mal withα-bond while Cel withβ-bond.Although Lac is also connected byβ-glycosidic bond,one of the repeating glucosides was changed to galactoside which contains axial OH group at C-4 instead of the equatorial OH of glucoside.

    The synthetic procedures were shown in Scheme 1.Disaccharide-hydrazide was synthesized based on three disaccharides.Firstly,through the classical method of protecting hydroxyl groups on disaccharide,all hydroxyl groups were protected with acetyl groups (Ac-).Then the protected 1–hydroxyl group was deprotected by hydrazine acetate.Trichloroacetonitrile was used to activate the 1–hydroxyl group,followed by further reaction with ethyl 3-hydroxypropionate catalyzed by trimethylsilyl trifluoromethanesulfonate (TMSOTf) at -35 °C.The protective group of the hydroxyls was removed by hydrazine hydrate and disaccharide- hydrazide was finally obtained.After three kinds of disaccharide-hydrazide were synthesized,the Azo-CHO was synthesized from 4-hydroxybenzaldehyde and 4-phenylazophenol in 4 steps according to our previous works [18–21].At last,mixing Azo-CHO and disaccharide- hydrazide at the molar ratio of 1:1 in DMSO for 4 days until the complete formation of azo-disaccharide and then stored as original solution,during which,ultrasound or heating was used to accelerate the formation speed [22].The detailed synthetic procedures were reported in Supporting information.

    Fig.2.The TEM images (a–c) and DLS images (d–f) of Azo-Lac,Azo-Cel,Azo-Mal.

    The supra-amphiphile azo-disaccharide was assembled by a selective solvent method [22].The self-assembly study of azodisaccharide supra-amphiphile was performed in water.The selfassembled samples were prepared by adding 7.5 mL deionized water into 1 mL of azo-disaccharide (1 mg/mL in DMSO) under room temperature.Syringe pump was used,in order to maintain a constant dripping rate of 15 mm/h under vigorous stirring.Then the mixture was dialyzed (MWCO 1000) in deionized water to remove the extra DMSO.The final concentration of the samples was fixed at 0.1 mg/mL by adding deionized water.To visualize the morphologies of assemblies,TEM and DLS were used to characterize the morphologies.Tracking by TEM,vesicles with diameters of about 200–300 nm and wall thickness of about 50 nm was observed for Azo-Lac,indicating that a multi-layered vesicle structure has been obtained (Fig.2a).The DLS analysis showed an average particle size of about 260 nm,which was consistent with the TEM results (Fig.2d).For Azo-Cel,a surprising result was captured with TEM,i.e.,a sunflower-like morphology.It was a hollow spherical structure with a diameter ranging from 400 nm to 500 nm (Fig.2b).The DLS result showed the particle size was about 406 nm,which was consistent with the TEM observations (Fig.2e).As for Azo-Mal,the DLS result showed that there indeed existed assembly formation (Fig.2f),but only some irregular structures were observed from TEM (Fig.2c).It can be seen that the tiny differences between the three disaccharide-containing molecules have given three quite different morphologies of assemblies.From this point,the configurational isomerization effect should have a great influence on selfassembly behaviors.

    Scheme 1.Synthetic procedure of azo-disaccharide (Azo-Lac for example,Azo-Cel and Azo-Mal are shown in Supporting information).

    Fig.3.Molecular dynamics simulations of bilayer formed by azo-disaccharide.(a–c) Side-view of initial bilayer for Azo-Lac,Azo-Cel and Azo-Mal.(d–f) Side-view of simulation results by Azo-Lac,Azo-Cel and Azo-Mal after relaxation for 20 ns,water molecules are hidden for clarity.

    The Azo-Lac,Azo-Cel and Azo-Mal in aqueous solution have shown completely different self-assembly behaviors.As discussed above,the only difference between these three supra-amphiphiles is the slight difference in the configuration of the disaccharide,Mal and Cel have different glycosidic bonds,and Lac is formed by a different glycoside with slight difference.So,it is important to have a deep insight into the relations between disaccharide conformation and the assembled morphologies.The theoretical and computational methods have always been expected in predicting the required structural features for specific applications and unveiling the mechanism behind.In this work,molecular dynamics simulations have been performed to investigate the assembly morphologies of our supra-amphiphiles.

    For all three molecules,a pre-assembled lipid bilayer was constructed and placed in the center of a rectangular box with a dimension of 10×10×18 nm3.Each lipid bilayer consisting of 100 molecules was constructed using Packmol program [23].The disaccharide groups were initially set to point outward and the hydrophobic chains were set to point to the internal.Then,about 10,000 water molecules were filled into the vacuum at both sides of the bilayer.We constructed the initial structures of all lipid bilayers by mimicking the equilibrated and interleaved bilayer structure.The initial configurations of the three lipid bilayers were shown in Figs.3a–c and their top views are shown in Figs.S1d–f (Supporting information).The structures of azo-disaccharide were shown in Figs.S1a–c (Supporting information).The simulations were performed using GROMACS package [24–27] with the CHARMM36 force field [28].Water molecules were described using the simple point charge/extend (SPC/E) model [29].The energy minimization was firstly performed for each system using the steep descent method.Then,20 ns simulations under NPT ensemble were performed for the three systems at 298 K and 1 atm (the potential diagram see Fig.S2d in Supporting information).According to the classical method in simulating bilayer models,the semiisotropic pressure coupling was applied during the simulation,in which the simulation box in z and x,y directions was rescaled individually [30,31].

    Fig.4.Surface properties of bilayer MD simulation (a–c) Surface properties of Azo-Lac,Azo-Cel and Azo-Mal after relaxation for 20 ns.(d–f) Randomly selected configurations of the Azo-Lac,Azo-Cel and Azo-Mal molecules after relaxation for 20 ns.(g–i) Top-views of Figs.d–f.The disaccharide rings on the surface of Azo-Lac are shown in orange,the disaccharide rings on the surface of Azo-Cel are shown in red,and the disaccharide rings on the surface of Azo-Mal are shown in yellow.

    The variations of the box lengths during the simulations were shown in Figs.S1g–i (Supporting information).The fluctuation of the data in Figs.S5g–i,showed that all systems had achieved equilibrium.The simulated side-view results of the azo-disaccharides were shown in Figs.3d–f.At first glance,there were no evident differences among the three models.However,As described above,those three disaccharides show different conformations in water.And those differences in conformation may have a great influence on the spatial arrangement of disaccharide molecules.To further investigate the different assembly morphologies,the molecular surface model was shown for each model (Figs.4a–c).The disaccharide rings of Azo-Lac are mainly evenly distributed on the surface of the bilayer,while the aggregations of Azo-Cel and Azo-Mal disaccharide rings are obvious,especially for the Azo-Mal model.Figs.4d–i showed the detailed configuration of azo-disaccharides,which were selected randomly from the final configurations.It shows clearly that the disaccharide rings of Azo-Lac molecules distribute evenly and the distances among them are relatively large,indicating the weak interactions among the disaccharide rings.Whereas,the Azo-Cel and Azo-Mal molecules show obvious aggregation due to the strong interactions of disaccharide rings.

    The interactions between disaccharide rings could be related to the hydrogen bonds.Under the same simulation protocol,the hydrogen bond numbers could provide distinguishable information.So,the number variations of hydrogen bonds during simulation time were counted to confirm the difference in their interactions(Fig.5a).The number of hydrogen bonds in the three models remained stable during the last period of the simulation,which confirmed the equilibrium of the simulated models again.Obviously,the H–bond number of Azo-Mal was the largest,and the number of Azo-Lac was the lowest.The average number of hydrogen bonds on the surface of Azo-Lac,Azo-Cel and Azo-Mal are approximately 242,279 and 327,respectively (Fig.5b).In addition,the numbers of intramolecular and intermolecular hydrogen bonds in various azodisaccharide bilayers were also counted.The final results showed that the number of intramolecular hydrogen bonds in the three azo-disaccharide bilayer models were almost the same.But the number of intermolecular hydrogen bonds between the three azodisaccharide bilayers were quite different.The partially enlarged area also showed the same trend (Figs.S2a–c in Supporting information).The average number of intermolecular hydrogen bonds on the surface of Azo-Lac,Azo-Cel and Azo-Mal were approximately 81,119,and 164,respectively (Fig.6).Therefore,we could infer that,due to the subtle difference in the configuration of the three azo-disaccharide,there was a huge difference in the number of intermolecular hydrogen bonds on their surface.The integral lifetime of the hydrogen bonds in the three azo-disaccharide bilayers was also calculated.The integral lifetime of Azo-Lac,Azo-Cel and Azo-Mal are approximately 201.71 ps,228.47 ps and 265.22 ps,respectively (Table S1 in Supporting information).This indicated that,in Azo-Mal model,apart from displaying the largest number of Hbond,it also showed the strongest binding strength of H-bonds among the three amphiphiles.Therefore,even distribution of hydrogen bonds on the surface could lead to smoother surface in assemblies,and the greater the number of hydrogen bonds,the stronger the aggregation tendency.Based on those simulation results,Azo-Lac with fewer hydrogen bonds forms smooth vesicles,Azo-Cel with more hydrogen bonds has spikes on the surface of vesicles,and Azo-Mal,due to the strongest aggregation effect,is unable to form a regular morphology.It is consistent with the TEM data that Azo-Mal cannot give a regular morphology.

    Fig 5.(a) Numbers of hydrogen bond plotted versus simulation time.(b) Average number of hydrogen bonds on the surface of Azo-Lac,Azo-Cel and Azo-Mal.

    Fig 6.Numbers of intermolecular hydrogen bond and intramolecular hydrogen bond in Azo-Lac,Azo-Cel and Azo-Mal.

    The morphology of each disaccharide molecule was also studied after the simulation reached equilibrium (Fig.S2 in Supporting information).In a Lac molecule,the galactoside and glucoside rings prefer to form a bent structure (Scheme 2),resulting in a relatively large steric hindrance.This bent structure facilitates the formation of intramolecular H-bonds between the galactoside and glucoside rings,leading to the formation of intramolecular hydrogen bond.Therefore,due to the bent structure,it is hard for two Lac molecules to interact with each other.In addition,the conformation of the Lac exhibits a prominent hydrophobic interface that is inclined to interact with the hydrophobic blocksviahydrophobicity.The Lac molecules prefer to distribute evenly at the surface.Whereas,the two adjacent glucoside rings of Mal tend to be coplanar,which facilitates intermolecular interactions among Mal molecules.Thus,the number of intermolecular hydrogen bonds increases,which also promotes the saccharide packing at the surface.As to Azo-Mal model,the two glucoside rings in a Cel molecule arranged crosswise.The steric hindrance strength of Cel molecule lies between Lac and Mal.So,the aggregation degree of Azo-Cel model is slightly lower than that of Azo-Mal model but higher than that of Azo-Lac model.Based on the above analysis,we plotted Scheme 2 to illustrate the relationship between the microstructure and macro morphology.

    In the simulation part,we only paid attention to the influence of the disaccharide end on the morphology,but azobenzene can be transformed betweentransandcisconformations by UV or visible light.Previous works by us and some other groups have proved the assembly of azobenzene-containing molecules can be rapidly regulated by photo-irradiation [32,33].So,we expected possible morphology transformation upon irradiation of UV light.Under UV irradiation at 365 nm for 0.5 h,UV–vis spectroscopy was first used to investigate the photo responsiveness of the azo-disaccharide supraamphiphile.The UV data showed that azobenzene part had been transformed from the stabletransconformation to thecisconformation after 0.5 h of UV light irradiation (Figs.S3c and d in Supporting information).Meanwhile,the morphology after irradiation was characterized by TEM (Figs.S3a and b in Supporting information).The UV data showed that whencis-azobenzene was irradiated by visible light,it returned to thetransconformation (Figs.S3e and f in Supporting information).The TEM and UV–vis results showed that UV light irradiation had no significant effect on the morphology of the two assemblies.This indicated thecisortransconformations of the azobenzene group had no obvious influence on the morphologies,which also confirmed that the morphology of assemblies is mainly controlled by the disaccharide ends.

    In summary,the assemblies of disaccharide-based supraamphiphiles can be adjusted to diverse complex morphologies by varying the structures of disaccharide ends,which was demonstrated by both experiments and simulation.Three supraamphiphiles with different disaccharides have been prepared,which present significant differences in assembly morphologies due to minor structural variations of disaccharides.Although the three disaccharides have similar chemical structures,they have very different self-assembly behavior.Simulation results revealed that a subtle structural change could greatly affect the related Hbond formation,then the packing of disaccharide molecule,finally the assembly morphology.The morphologies were greatly controlled by the packing of the disaccharide moieties.

    Scheme 2.The cartoon illustration of the self-assembly process of the azo-disaccharide supra-amphiphile.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    G.Chen thanks the National Natural Science Foundation of China (NSFC,Nos.52125303,51721002,91956127 and 21975047)for financial support.R.Zhang thanks the National Natural Science Foundation of China (NSFC,Nos.21674114 and 91956127) for financial support.This work is also supported by the Shanghai Municipal Science and Technology Major Project (No.2018SHZDZX01) and ZJ Lab.

    Supplementary materials

    Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.cclet.2022.05.080.

    乱人视频在线观看| 高清毛片免费看| 亚洲图色成人| 亚洲精品一区蜜桃| 国产综合懂色| 美女被艹到高潮喷水动态| 国产又色又爽无遮挡免| 亚洲最大成人中文| 高清午夜精品一区二区三区| av网站免费在线观看视频 | 亚洲国产精品成人综合色| 少妇熟女欧美另类| 可以在线观看毛片的网站| 免费少妇av软件| 成人av在线播放网站| 亚洲精品日韩av片在线观看| 尤物成人国产欧美一区二区三区| 亚洲自偷自拍三级| 插逼视频在线观看| 亚洲av不卡在线观看| 青春草国产在线视频| 国产亚洲精品久久久com| av女优亚洲男人天堂| 中文字幕免费在线视频6| 中国美白少妇内射xxxbb| 三级国产精品欧美在线观看| 中文精品一卡2卡3卡4更新| 欧美xxxx黑人xx丫x性爽| 久热久热在线精品观看| 欧美+日韩+精品| 日韩一本色道免费dvd| 国产成人精品婷婷| 国产黄色视频一区二区在线观看| 国产欧美日韩精品一区二区| 国产伦精品一区二区三区四那| 久久精品国产亚洲av涩爱| 2022亚洲国产成人精品| 亚洲av电影在线观看一区二区三区 | 亚洲av男天堂| 五月玫瑰六月丁香| 国产精品久久久久久精品电影| 99re6热这里在线精品视频| 18禁在线无遮挡免费观看视频| 美女主播在线视频| 国产伦在线观看视频一区| 天堂影院成人在线观看| 99热全是精品| 成人午夜精彩视频在线观看| 人人妻人人澡欧美一区二区| 久久精品综合一区二区三区| 99热全是精品| 亚洲最大成人中文| 婷婷色综合大香蕉| 久久99热这里只频精品6学生| 干丝袜人妻中文字幕| 国产成人午夜福利电影在线观看| 深夜a级毛片| 欧美一区二区亚洲| 18禁动态无遮挡网站| 少妇熟女aⅴ在线视频| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 国产美女午夜福利| 欧美人与善性xxx| 日韩欧美精品v在线| 一级毛片电影观看| 日本爱情动作片www.在线观看| 美女黄网站色视频| 寂寞人妻少妇视频99o| 久久久久久久久中文| 日韩 亚洲 欧美在线| 国产精品三级大全| 热99在线观看视频| 自拍偷自拍亚洲精品老妇| 日韩制服骚丝袜av| 最近手机中文字幕大全| 啦啦啦韩国在线观看视频| 国产免费一级a男人的天堂| 午夜激情久久久久久久| 久久久久精品性色| 色5月婷婷丁香| 深夜a级毛片| 久久97久久精品| 极品教师在线视频| 精品不卡国产一区二区三区| 国产黄a三级三级三级人| 国产美女午夜福利| 日韩,欧美,国产一区二区三区| 亚洲熟女精品中文字幕| 国产精品久久视频播放| 国产淫语在线视频| 国产爱豆传媒在线观看| 日本免费a在线| 婷婷色av中文字幕| 女人久久www免费人成看片| 亚洲精品成人久久久久久| 26uuu在线亚洲综合色| 亚洲久久久久久中文字幕| 日韩欧美一区视频在线观看 | 精品国产一区二区三区久久久樱花 | 久久久久久久久久久丰满| 内地一区二区视频在线| 一级毛片黄色毛片免费观看视频| 国产白丝娇喘喷水9色精品| 秋霞伦理黄片| 一个人看视频在线观看www免费| 精品不卡国产一区二区三区| 亚洲成人中文字幕在线播放| 久久精品综合一区二区三区| 免费在线观看成人毛片| 91精品一卡2卡3卡4卡| 久久精品人妻少妇| av在线天堂中文字幕| 免费观看a级毛片全部| 成人特级av手机在线观看| 水蜜桃什么品种好| 国产高清三级在线| 中文字幕av成人在线电影| 亚洲自偷自拍三级| 久久久色成人| kizo精华| 婷婷色综合大香蕉| 最近手机中文字幕大全| 欧美xxxx性猛交bbbb| 亚洲国产最新在线播放| 伦精品一区二区三区| 美女被艹到高潮喷水动态| 蜜桃久久精品国产亚洲av| 国产色爽女视频免费观看| 国产成人午夜福利电影在线观看| 免费av观看视频| 波野结衣二区三区在线| 日本-黄色视频高清免费观看| 日日啪夜夜撸| 99久久精品热视频| 精品久久久久久久久av| 观看美女的网站| 亚洲va在线va天堂va国产| 亚洲国产精品成人久久小说| 青春草视频在线免费观看| 一区二区三区高清视频在线| 我的女老师完整版在线观看| 久久精品国产亚洲网站| 欧美人与善性xxx| 国产一区二区在线观看日韩| 99久国产av精品国产电影| 亚洲丝袜综合中文字幕| 亚洲熟女精品中文字幕| 国语对白做爰xxxⅹ性视频网站| 久久久久久久亚洲中文字幕| 免费看a级黄色片| 国产黄片视频在线免费观看| 丝袜喷水一区| 国产三级在线视频| 精品久久久精品久久久| 国产午夜精品一二区理论片| 久久久久精品久久久久真实原创| 老司机影院成人| 成人午夜精彩视频在线观看| 少妇被粗大猛烈的视频| 国模一区二区三区四区视频| 午夜免费男女啪啪视频观看| 国内揄拍国产精品人妻在线| 51国产日韩欧美| 国产色婷婷99| 欧美日韩精品成人综合77777| 国产精品久久视频播放| 在线播放无遮挡| 亚洲av成人精品一二三区| 天天躁日日操中文字幕| 十八禁国产超污无遮挡网站| 国产亚洲5aaaaa淫片| 午夜亚洲福利在线播放| 久久久久久久亚洲中文字幕| 日本午夜av视频| 丰满少妇做爰视频| 欧美日本视频| 免费不卡的大黄色大毛片视频在线观看 | 国产成人a∨麻豆精品| av在线观看视频网站免费| 国产精品无大码| 高清av免费在线| av在线播放精品| 乱系列少妇在线播放| 国产乱人偷精品视频| 国产伦一二天堂av在线观看| 天堂中文最新版在线下载 | 天天一区二区日本电影三级| 1000部很黄的大片| 26uuu在线亚洲综合色| 国产黄a三级三级三级人| 视频中文字幕在线观看| 天堂影院成人在线观看| 国产精品国产三级专区第一集| 日本爱情动作片www.在线观看| 高清在线视频一区二区三区| 只有这里有精品99| 大片免费播放器 马上看| 天堂网av新在线| 日本与韩国留学比较| 亚洲精品影视一区二区三区av| 黄色欧美视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品久久久com| 18+在线观看网站| 免费播放大片免费观看视频在线观看| 最近2019中文字幕mv第一页| 最近中文字幕2019免费版| 国产一区二区亚洲精品在线观看| 高清在线视频一区二区三区| 美女大奶头视频| 亚洲综合精品二区| 精品少妇黑人巨大在线播放| 国产国拍精品亚洲av在线观看| 婷婷色av中文字幕| 22中文网久久字幕| 久久6这里有精品| 联通29元200g的流量卡| 搡女人真爽免费视频火全软件| 久久久久网色| 三级男女做爰猛烈吃奶摸视频| 男女边吃奶边做爰视频| 国产亚洲av片在线观看秒播厂 | 国模一区二区三区四区视频| 天天一区二区日本电影三级| 日韩三级伦理在线观看| 噜噜噜噜噜久久久久久91| 亚洲va在线va天堂va国产| 特大巨黑吊av在线直播| 赤兔流量卡办理| 国产熟女欧美一区二区| 中文欧美无线码| 日日撸夜夜添| 免费不卡的大黄色大毛片视频在线观看 | 欧美另类一区| 国产成人精品婷婷| 精品熟女少妇av免费看| 亚洲欧美日韩无卡精品| 精华霜和精华液先用哪个| 高清日韩中文字幕在线| or卡值多少钱| 亚洲精品国产av成人精品| 久久久色成人| 中文字幕制服av| 人人妻人人看人人澡| 三级经典国产精品| 日本与韩国留学比较| 久久久久久久午夜电影| 老司机影院成人| 国产日韩欧美在线精品| 亚洲,欧美,日韩| 国产精品一区二区性色av| 久久99蜜桃精品久久| 少妇被粗大猛烈的视频| 日韩欧美三级三区| 99九九线精品视频在线观看视频| 卡戴珊不雅视频在线播放| 国产亚洲一区二区精品| a级一级毛片免费在线观看| 免费观看在线日韩| 亚洲怡红院男人天堂| 卡戴珊不雅视频在线播放| 女的被弄到高潮叫床怎么办| 国产成人午夜福利电影在线观看| 天天躁夜夜躁狠狠久久av| 亚洲av免费在线观看| 欧美激情在线99| 内地一区二区视频在线| 国产探花极品一区二区| 五月伊人婷婷丁香| 乱人视频在线观看| 一本一本综合久久| 亚洲精品国产av蜜桃| 亚洲人成网站在线观看播放| 熟女电影av网| 国产av在哪里看| xxx大片免费视频| 好男人视频免费观看在线| 国产探花在线观看一区二区| 免费av毛片视频| 人妻系列 视频| 欧美日韩一区二区视频在线观看视频在线 | 男女边摸边吃奶| 春色校园在线视频观看| 在线a可以看的网站| 欧美性感艳星| 午夜免费男女啪啪视频观看| 久久久久久国产a免费观看| 中文字幕人妻熟人妻熟丝袜美| 久久99精品国语久久久| 在线 av 中文字幕| 99久久精品一区二区三区| 一区二区三区高清视频在线| 18禁在线播放成人免费| 九九爱精品视频在线观看| 亚洲欧美精品自产自拍| 高清日韩中文字幕在线| 国产人妻一区二区三区在| 91精品一卡2卡3卡4卡| 欧美bdsm另类| 亚洲国产精品成人综合色| 久久久久久久亚洲中文字幕| 久久精品人妻少妇| 日日摸夜夜添夜夜爱| av线在线观看网站| 天天躁夜夜躁狠狠久久av| 日本av手机在线免费观看| 午夜福利成人在线免费观看| 免费av不卡在线播放| 日韩三级伦理在线观看| 水蜜桃什么品种好| 又爽又黄a免费视频| 色网站视频免费| 亚洲av免费在线观看| 久99久视频精品免费| 国产在线男女| 久久99蜜桃精品久久| 亚洲国产精品成人综合色| 老女人水多毛片| 国产免费一级a男人的天堂| 亚洲最大成人av| 嘟嘟电影网在线观看| .国产精品久久| 最近最新中文字幕免费大全7| 亚洲欧美中文字幕日韩二区| 中文字幕av在线有码专区| 久久99热这里只有精品18| 91精品国产九色| 亚洲乱码一区二区免费版| 国产人妻一区二区三区在| 亚洲欧美一区二区三区黑人 | 亚洲国产高清在线一区二区三| av卡一久久| 一个人看视频在线观看www免费| 日韩强制内射视频| 久久亚洲国产成人精品v| 丰满乱子伦码专区| av专区在线播放| 国产精品人妻久久久影院| 国产 一区 欧美 日韩| 美女主播在线视频| 午夜福利成人在线免费观看| 亚洲欧洲国产日韩| 色尼玛亚洲综合影院| 日韩中字成人| 亚洲不卡免费看| 亚洲欧美日韩无卡精品| 特大巨黑吊av在线直播| 久久久久精品久久久久真实原创| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区综合在线观看 | 日韩人妻高清精品专区| 三级毛片av免费| 午夜福利网站1000一区二区三区| 亚洲av免费高清在线观看| 高清午夜精品一区二区三区| 一本久久精品| 夜夜看夜夜爽夜夜摸| 欧美xxxx黑人xx丫x性爽| 免费无遮挡裸体视频| 91久久精品国产一区二区成人| 黄色配什么色好看| 天天躁日日操中文字幕| 精品国产露脸久久av麻豆 | 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 伦精品一区二区三区| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 在线观看美女被高潮喷水网站| av在线亚洲专区| 嫩草影院精品99| 午夜免费男女啪啪视频观看| 亚洲最大成人av| 在线播放无遮挡| 欧美高清成人免费视频www| 亚州av有码| 国产亚洲av片在线观看秒播厂 | 久久久久久久国产电影| 精品国产三级普通话版| 国产乱人偷精品视频| 日日撸夜夜添| 中文字幕久久专区| 91久久精品国产一区二区三区| 天堂av国产一区二区熟女人妻| 午夜激情久久久久久久| 80岁老熟妇乱子伦牲交| 国产精品精品国产色婷婷| 尾随美女入室| 视频中文字幕在线观看| 白带黄色成豆腐渣| 欧美成人一区二区免费高清观看| 欧美xxⅹ黑人| 97精品久久久久久久久久精品| 男人舔奶头视频| 久久精品国产亚洲网站| 日本午夜av视频| 精品99又大又爽又粗少妇毛片| 亚洲欧洲国产日韩| 又黄又爽又刺激的免费视频.| 亚洲一级一片aⅴ在线观看| 久久精品国产鲁丝片午夜精品| 日本色播在线视频| 内地一区二区视频在线| 爱豆传媒免费全集在线观看| 一级毛片aaaaaa免费看小| 精品久久久久久久久av| 中文资源天堂在线| 日本爱情动作片www.在线观看| 色综合色国产| 久久久久九九精品影院| 国产免费一级a男人的天堂| 乱系列少妇在线播放| 婷婷色麻豆天堂久久| 欧美丝袜亚洲另类| 午夜福利在线在线| 韩国高清视频一区二区三区| 日本午夜av视频| 九九爱精品视频在线观看| 久久久久久久久大av| 国产 一区 欧美 日韩| 天堂√8在线中文| 日韩欧美一区视频在线观看 | 看十八女毛片水多多多| 97热精品久久久久久| 精品国产露脸久久av麻豆 | 国产老妇女一区| 禁无遮挡网站| 亚洲av国产av综合av卡| 国产精品一二三区在线看| 亚洲国产精品sss在线观看| 18禁动态无遮挡网站| 日韩伦理黄色片| 极品少妇高潮喷水抽搐| 色哟哟·www| 久久久久免费精品人妻一区二区| 午夜福利视频精品| 国产有黄有色有爽视频| 简卡轻食公司| 亚洲熟女精品中文字幕| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲91精品色在线| av在线亚洲专区| 丝袜美腿在线中文| 国产熟女欧美一区二区| 最近的中文字幕免费完整| 只有这里有精品99| 丝瓜视频免费看黄片| 亚洲一级一片aⅴ在线观看| 欧美潮喷喷水| 一级二级三级毛片免费看| 国产成人精品福利久久| 亚洲最大成人手机在线| 91午夜精品亚洲一区二区三区| 女人被狂操c到高潮| 少妇丰满av| 极品少妇高潮喷水抽搐| 久久精品综合一区二区三区| 少妇熟女aⅴ在线视频| 日韩欧美精品免费久久| 永久免费av网站大全| 亚洲精品第二区| 深夜a级毛片| 亚洲不卡免费看| 中文精品一卡2卡3卡4更新| 亚洲av中文字字幕乱码综合| 欧美日韩视频高清一区二区三区二| 三级经典国产精品| 亚洲国产色片| 欧美日韩在线观看h| 亚洲精品乱码久久久v下载方式| 国产精品.久久久| 午夜亚洲福利在线播放| 插逼视频在线观看| 亚洲成人久久爱视频| 亚洲精品日韩av片在线观看| 在线观看人妻少妇| or卡值多少钱| 国产精品一区二区在线观看99 | 夫妻性生交免费视频一级片| 亚洲精品日本国产第一区| 2021天堂中文幕一二区在线观| 成人亚洲精品av一区二区| 九九久久精品国产亚洲av麻豆| 真实男女啪啪啪动态图| 国产精品精品国产色婷婷| 亚洲经典国产精华液单| 熟女人妻精品中文字幕| 美女大奶头视频| 国产 一区 欧美 日韩| 综合色av麻豆| 欧美不卡视频在线免费观看| 在现免费观看毛片| 亚洲精品,欧美精品| 国产在线一区二区三区精| 亚洲内射少妇av| 欧美人与善性xxx| 18禁裸乳无遮挡免费网站照片| 久久久亚洲精品成人影院| 久久草成人影院| 日韩欧美三级三区| 人妻一区二区av| 能在线免费看毛片的网站| 男女啪啪激烈高潮av片| 国产91av在线免费观看| 色吧在线观看| 亚洲精品一二三| 色播亚洲综合网| av.在线天堂| 建设人人有责人人尽责人人享有的 | 国产精品伦人一区二区| 久久99热这里只有精品18| 成人午夜精彩视频在线观看| 91精品伊人久久大香线蕉| 成人毛片a级毛片在线播放| 亚洲精品第二区| 日本wwww免费看| 免费人成在线观看视频色| 一级毛片电影观看| av福利片在线观看| 亚洲精品成人av观看孕妇| 亚洲国产av新网站| 午夜免费男女啪啪视频观看| videossex国产| 国产高潮美女av| 一本久久精品| 亚洲欧美精品专区久久| 99热这里只有是精品50| 搡老妇女老女人老熟妇| 夫妻性生交免费视频一级片| 免费在线观看成人毛片| av免费在线看不卡| 久久亚洲国产成人精品v| 最后的刺客免费高清国语| 国产av在哪里看| 老女人水多毛片| 免费大片18禁| 九九在线视频观看精品| 欧美丝袜亚洲另类| 爱豆传媒免费全集在线观看| 国产探花极品一区二区| 亚洲国产精品成人久久小说| 男人舔女人下体高潮全视频| a级一级毛片免费在线观看| 男女视频在线观看网站免费| 午夜日本视频在线| 老师上课跳d突然被开到最大视频| 欧美性猛交╳xxx乱大交人| 日本免费在线观看一区| 男人狂女人下面高潮的视频| 天堂√8在线中文| 99久久九九国产精品国产免费| 国产精品美女特级片免费视频播放器| 久久精品久久久久久噜噜老黄| 秋霞伦理黄片| 亚洲婷婷狠狠爱综合网| 午夜日本视频在线| 国产亚洲av片在线观看秒播厂 | 80岁老熟妇乱子伦牲交| 好男人视频免费观看在线| 亚洲精品久久久久久婷婷小说| 中文在线观看免费www的网站| 欧美一区二区亚洲| 肉色欧美久久久久久久蜜桃 | 国产精品一区二区三区四区免费观看| 亚洲国产色片| 日韩制服骚丝袜av| 国产午夜精品一二区理论片| 18禁裸乳无遮挡免费网站照片| 91aial.com中文字幕在线观看| 欧美一级a爱片免费观看看| 欧美日韩视频高清一区二区三区二| 国产黄片美女视频| 国产老妇伦熟女老妇高清| 午夜福利在线观看免费完整高清在| 午夜福利视频精品| 伊人久久国产一区二区| 国产欧美另类精品又又久久亚洲欧美| 成人综合一区亚洲| 亚洲内射少妇av| 国产探花极品一区二区| 最近手机中文字幕大全| 亚洲激情五月婷婷啪啪| 国产v大片淫在线免费观看| 亚洲av成人精品一二三区| 联通29元200g的流量卡| 深爱激情五月婷婷| 成人亚洲精品一区在线观看 | 国产精品不卡视频一区二区| 国精品久久久久久国模美| 肉色欧美久久久久久久蜜桃 | 三级毛片av免费| 只有这里有精品99| 蜜桃久久精品国产亚洲av| 美女黄网站色视频| 精品久久久噜噜| 国产视频首页在线观看| 免费大片18禁| 亚洲欧美日韩东京热| 六月丁香七月| 国产免费福利视频在线观看| 性色avwww在线观看| 1000部很黄的大片| 建设人人有责人人尽责人人享有的 | 亚洲精品视频女| 亚洲av不卡在线观看| 黑人高潮一二区| 成人午夜精彩视频在线观看| 69人妻影院| av天堂中文字幕网| 成年人午夜在线观看视频 | 一级毛片久久久久久久久女| 夫妻性生交免费视频一级片| 国产精品一区www在线观看| 国产黄频视频在线观看|