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

    Coarse-grained simulations on interactions between spectrins and phase-separated lipid bilayers?

    2021-06-26 03:31:38XueguiLin林雪桂XiaojieChen陳曉潔andQingLiang梁清
    Chinese Physics B 2021年6期

    Xuegui Lin(林雪桂) Xiaojie Chen(陳曉潔) and Qing Liang(梁清)

    1Center for Statistical and Theoretical Condensed Matter Physics and Department of Physics,Zhejiang Normal University,Jinhua 321004,China

    2Department of Industrial Automation,Guangdong Polytechnic College,Zhaoqing 526100,China

    Keywords: protein,lipid bilayer,interaction,molecular dynamics simulation

    1. Introduction

    Erythrocyte is a specific cell without nucleus and most organelles.[1,2]The plasma membrane of erythrocyte is associated with a cytoskeleton composed of erythroid spectrins(referred to herein as spectrins) and some other membrane proteins.[2–4]The interaction between spectrins and erythrocyte membrane is considered as a key factor determining the erythrocyte structural stability and flexibility which are of great importance for the realization of the biological functions of erythrocyte.[5–11]Therefore,it is essential to understand the microscopic mechanism of interactions between spectrins and erythrocyte membrane from the molecular scale.

    The single chain of spectrin skeleton consists of two heterodimers, theαandβspectrins, which form antiparallel heterotetramers by laterally interacting with each other.[3,4,12]In the cytoskeleton of the heathy erythrocyte, theαandβspectrins interact with high affinity at the tetramerization site and form a stable network.[13,14]However,if some mutations of amino acids occurs at the tetramerization site, the affinity of theαandβspectrins may be considerably influenced and the stability of the erythrocyte cytoskeleton will be reduced, resulting in the deformation of erythrocyte and some blood diseases such as hereditary elliptocytosis (HE), hereditary pyropoikilocytosis (HPP), and hereditary spherocytosis(HS).[15–22]

    The predominant lipids of the plasma membrane of erythrocyte are cholesterol (CHOL) and phospholipids, and the molar ratio of CHOL: phospholipids is about 0.90. The main phospholipids include phosphatidylcholine (PC), phosphatidylethanolamine (PE), sphingomyelin (SM), and phosphatidylserine(PS).In addition,the lipids distribute asymmetrically in the two leaflets of the membrane. PS and PE almost exclusively distribute in the inner leaflet while the outer leaflet of the membranes is mainly composed of PC and SM.[23,24]

    The spectrin skeleton anchors to the erythrocyte membrane mainly though connecting with some integral/peripheral membrane proteins or protein complexes.[3,4,17,23,25–35]However,more and more evidences indicated that the direct interactions between spectrins and lipids could also be crucial for the structural stability and function of erythrocyte.[36–40]It has been shown that there are some specific binding sites on spectrins for PC,PE,and PS lipids in the inner leaflet of erythrocyte membrane,and theαandβspectrins interact differently with the same kind of lipids.[41–43]Additionally, CHOL and phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2, referred to as PIP2 in this work)also play important roles in the interactions of spectrins with erythrocyte membrane.[44–47]

    Furthermore, spectrins can also influence the nano-sized domain formation in the erythrocyte membrane and exhibit different affinities with the liquid-ordered(Lo)domains(e.g.,lipid rafts)and the liquid-disordered(Ld)domains in the erythrocyte membrane.[48–54]Therefore,the interactions between spectrins and the plasma membrane of erythrocyte are extremely complicated and still remain elusive. More efforts are necessary to reveal the microscopic mechanism of the interactions between spectrins and erythrocyte membrane, in particular,at the molecular scale.

    In this work, using coarse-grained molecular dynamics(MD) simulations, we investigate the interactions of spectrins with the phase-separated lipid bilayers composed of saturated lipids (e.g., dipalmitoyl-PC (DPPC), dipalmitoyl-PE (DPPE), or PIP2), unsaturated lipids (e.g., dilinoleoyl-PE (DIPE), dilinoleoyl-PC (DIPC), or palmitoyl-oleoyl-PS(POPS)),and CHOL.We try to reveal the microscopic mechanism of the affinities of spectrins with different lipids and Lo/Ld domains, and the effects of charged lipids and spectrin mutation on the spectrin–membrane interactions at the molecular scale. The simulations provide some new theoretical insights which are helpful for understanding the molecular structure of erythrocyte as well as the molecular mechanism of some relevant blood diseases.

    2. Model and methods

    2.1. System setup

    2.1.1. Lipid bilayers

    To capture the main structural and componential properties of erythrocyte membrane and simulate the lipid bilayers with relatively large size for time of microsecond(μs)length,we simplified the components of the plasma membrane of erythrocyte and built some ternary and quaternary lipid bilayers composed of the main components of erythrocyte membrane,PC, PE, PS, PIP2, and CHOL, with the coarse-grained Martini model.[55–57]The headgroups of PC and PE are zwitterionic while those of PS and PIP2 are anionic. PS carries one negative charge on the headgroup. PIP2 carries five negative charges in the unprotonated state since salt was added in our systems,[58]and its molecular structure was adopted from Ref. [59]. The coarse-grained molecular structures of the lipids used in our simulations are presented in Fig. 1(a)and the compositions of all the lipid bilayers are summarized in Table 1. Here, the tails of DPPC, DPPE,or PIP2 are fully saturated while those of DIPC, DIPE, or POPS are unsaturated. For all the lipid bilayers,the leaflet directly interacting with the spectrins is the inner leaflet and the other leaflet is the outer leaflet. The previous work has shown that the membrane composed of saturated lipids, unsaturated lipids, and cholesterol can spontaneously separate into Lo and Ld domains at 305 K.[60]In this work, the initial distribution of the lipids is homogeneous and the phase separation of the Lo and Ld domains spontaneously occurs in all the systems within several microseconds.

    Table 1. Overviews of the lipid compositions in the two leaflets of the lipid bilayers, the spectrin type in the corresponding system, and the simulation time. The spectrins are mutant with the mutation sites of W263R and W300R on the β17 spectrin in the last system marked with‘M’.

    2.1.2. Spectrin

    The structure of spectrin domains(PDB ID:3LBX)used in this work is shown in Fig. 1(b). The spectrins contain four homologous repeats: the repeats 0 and 1 of N-terminus of theαspectrins (α0 andα1, residues 19–158) and the repeats 16 and 17 of C-terminus of theβspectrins (β16 andβ17, residues 1900–1931 and 1934–2083 ).[12]The spectrin repeats interact head-to-head betweenα0 andβ17 through non-bonded interactions. The atomistic structure of the spectrins was coarse-grained with Martini model and the residues of the coarse-grained spectrins are renumbered from 1 to 322.Correspondingly,the residues 1–140 belong to theαspectrins while the residues 141–322 belong to theβspectrins in the coarse-grained model. Besides the wild-type spectrins, we also examined the interactions of the mutant spectrins with the lipid bilayer composed of DPPE/DIPC/CHOL.Here,we considered the mutant spectrins with the mutation sites of W263R and W300R (W2024R and W2061R in the original atomistic structure)on theβ17 spectrin as shown in Fig.1(c). Namely,the amino acids of tryptophan(W)at the positions of 263 and 300(2024 and 2061 in the original atomistic structure)in the wild-typeβ17 spectrin are replaced by arginine (R), respectively. This type of mutation is considered one of the causes of the blood disease of HE.[22]

    Fig. 1. (a) Coarse-grained (CG) molecular structures of the lipids used in this work. The lipids are shown in colors as follows: DPPC: red; DPPE:magenta;DIPC:pink;DIPE:blue;PIP2: yellow;POPS:green;CHOL:gray.This coloring scheme is used throughout the paper. (b)Atomistic and coarsegrained structures of the α (white)and β (cyan)spectrins. (c)Illustration of the mutation sites of W263R(yellow)and W300R(red)on the β17 spectrin.

    All the lipid bilayers with the spectrins were built using the script of INSANE provided by the Martini developers.[61]To save the simulation time for the adsorption of the spectrins onto the surface of the lipid bilayers,we put the spectrins near the surface of inner leaflets of the bilayers and the distances between them are about 1.4 nm–2.5 nm. The initial box size of the simulations is about 20 nm×20 nm×30 nm and periodic boundary conditions were used in all three directions. The size of the simulation box inzdirection was set relatively large to ensure the free adjustment of the spectrins.In all ternary lipid bilayers, the molar ratio of the saturated lipids, the unsaturated lipids, and CHOL molecules is 4:3:3.To examine the effect of charged lipids on the interactions between spectrins and the lipid bilayers,some saturated lipids in the inner leaflet of the ternary bilayers were substituted by the anionic lipids POPS or PIP2. For the convenience of comparison,we kept the charge number of the anionic lipids the same in the charged bilayers,i.e., the fractions of PIP2 and POPS are 2 mol% and 10 mol%, respectively, since the number of charged beads in the headgroup of PIP2 is five times larger than that of POPS.Each system contains about 8×104coarsegrained water beads to solvate the lipid bilayer. Furthermore,0.15-M salt was added into each system to mimic the physiological conditions and additional sodium ions were introduced to maintain electroneutrality of the systems containing anionic lipids(PIP2 or POPS).

    2.2. Coarse-grained MD simulations

    All the systems were simulated using Martini coarsegrained force field. Although ignoring the hydrogen-bonding interactions, some molecular details of lipids, and some secondary structural details of proteins,the Martini model,which used a four-to-one mapping scheme, is still one of the most successful models for predicting theμs-long dynamical properties of complex lipid bilayers and its interactions with proteins or nanoparticles.[56,62–65]In many cases, Martini model got good agreement with the atomistic models or experiments due to its good optimization on the parametrization of the lipid–lipid and lipid–protein interactions.[56,62,66]In the present simulations, the non-bonded interaction parameters were set as follows: the Coulomb interaction was shifted to zero between 0 nm and 1.2 nm and the Lennard–Jones interaction was shifted to zero between 0.9 nm and 1.2 nm. The NPT ensemble was used for the simulations. The temperature of all the simulations was maintained at 300 K using velocityrescaling coupling scheme with a relaxation time of 1 ps.[67]The bilayer(xy)plane and the normal(z)direction were separately coupled to a 1-bar pressure bath with a relaxation time of 2 ps by using the semi-isotropic Parrinello–Rahman pressure coupling scheme.[68,69]The compressibility of systems was set as 3×10?4bar?1(1 bar=105Pa)in both lateral and normal directions of the bilayer.After energy minimization and a 1-ns equilibrating simulation,except for the systems 3 and 4 which were simulated for 15μs,the other systems were simulated for 10μs. To ensure the reproducibility of the results, each simulation was run for 2–3 replicas with the same initial states of the lipid bilayers and proteins but different distances between them. The GROMACS package(version 2019.3)was used to perform the simulations with an integration time-step of 20 fs.

    2.3. Data analysis

    We used Visual Molecular Dynamics (VMD) to visualize the structures of the lipid bilayers and the spectrins.[70]To characterize the molecular details of the interactions between the spectrins and lipid bilayers and the structural stability of the spectrins,some tools from the GROMACS package were used to analyze the results. The tool ofgmx mindistwas adopted to calculate the numbers of contacts between the spectrins (or amino acids) and the headgroups of DPPC, DIPE,DPPE, DIPC, and CHOL molecules with a cutoff length of 0.6 nm as a function of simulation time. The tool ofgmx rmswas used to calculate the root-mean-square deviation(RMSD)of the spectrins. The RMSD is the most commonly-used function for characterizing the variation of the protein configuration at certain time(t)compared with the reference configuration of the protein at the initial time(t=0)and defined as

    Here,miis the mass of thei-th atom on the protein backbone,andri(t)is the position of thei-th atom at timet;Nis the total number of the atoms on the protein backbone;Mis the total mass of the atoms on the protein backbone andM=∑Ni=1mi.

    Additionally, to examine the effect of the spectrin binding on the structure of the lipid bilayer, we also analyzed the curvature of the lipid bilayers by projecting the coordinates of the lipid headgroups of the inner leaflet to a two-dimensional plane and fitting with a curve and the maximal curvature of the curve was used to estimate the curvature of the lipids.[71]

    3. Results and discussion

    3.1. Interactions between spectrins and lipid domains of the neutral lipid bilayers

    We first examine the interactions of the spectrins with the ternary bilayers only containing neutral lipids. We considered four lipid bilayers composed of DPPC/DIPC/CHOL, DPPC/DIPE/CHOL,DPPE/DIPC/CHOL, and DPPE/DIPE/CHOL, respectively,with the molar ratio of 4:4:3. Figure 2 shows the final structures of the 10-μs or 15-μs simulations on these lipid bilayers with the spectrins and the contact numbers of theαandβspectrins with the lipids during the simulations. The results show that the spectrins exhibit different interactions not only with PE and PC lipids but also with saturated and unsaturated lipids.

    For the ternary lipid bilayer composed of DPPC, DIPC,and CHOL (Fig. 2(a)), both theαandβspectrins contact weakly with the bilayer only near the end of the simulation time while they do not tend to contact with the bilayer during most of the simulation time, which indicates that the interaction between the spectrins and PC lipids is very weak. For the bilayers containing PE lipids (Figs. 2(b)–2(d)), the spectrins are steadily adsorbed on the surface of the bilayers. However,the spectrins interact with saturated and unsaturated PE lipids differently in the different lipid bilayers. In the lipid bilayer composed of DPPC/DIPE/CHOL, the spectrins are adsorbed on the surface of Ld domain enriched in unsaturated DIPE.In the bilayer composed of DPPE/DIPC/CHOL,the spectrins prefer to contact with the Lo domains enriched in the saturated DPPE. In the bilayer composed of DPPE/DIPE/CHOL,in most time of the simulation, the spectrins show no obvious preference with either DPPE or DIPE, and the contacts between the spectrins with DPPE and with DIPE are comparable. Only in the last 2μs,the number of contacts between the spectrins and DIPE becomes slightly larger than the number of contacts between the spectrins and DPPE, which implies that the spectrins tend to contact with the DIPE-enriched Ld domain.

    Fig. 2. Top views and lateral views of the final structures of the lipid bilayers composed of (a) DPPC/DIPC/CHOL, (b) DPPC/DIPE/CHOL,(c) DPPE/DIPC/CHOL, and (d) DPPE/DIPE/CHOL, with the spectrins, and the numbers (No.) of contacts between the α (β) spectrins and the lipid headgroups. For clarity,the structures of the lipid bilayers with the spectrins displayed and not displayed are both presented,and water and ions are not shown.

    The results of these four systems show that the spectrins exhibit a stronger interaction with PE lipids than with PC lipids,indicating that there are some specific binding sites for PE on the surface of the spectrins. This result agrees well with some experimental results.[41,42]Additionally, the affinity of the spectrins with the Lo and Ld domains is determined by the compositions of the lipid bilayers. For the bilayer of DPPE/DIPC/CHOL(Fig.2(c)),because of the stronger attraction between the spectrins and DPPE, the spectrins prefer to contact with the DPPE-enriched Lo domain. For the bilayers with DIPE-enriched Ld domains(Figs.2(b)and 2(d)),the spectrins tends to contact with the Ld domains. This is possibly because that the disordered packing of DIPE in the Ld domains facilitates the binding or insertion of the spectrins onto/into the lipid domains.[35]

    Furthermore, comparing the contact numbers of theαand theβspectrins with the lipids shown in Fig. 2, we find that the numbers of the contacts of theβspectrins with the lipids are obviously larger than those of theαspectrins with the lipids. Therefore, the binding ofβspectrins dominates the interactions between the spectrins and lipid bilayers,which also agrees well with the experimental result.[43]

    Fig. 3. Top views and lateral views of the final structures of the spectrin-bound lipid bilayers composed of (a) DPPC/DIPC/CHOL/PIP2,(b)DPPC/DIPC/CHOL/POPS,(c)DPPC/DIPE/CHOL/PIP2,and(d)DPPC/DIPE/CHOL/POPS,with the spectrins,and the numbers(No.) of contacts between the spectrins and the lipid headgroups.

    3.2. Effects of anionic lipids on the interactions between spectrins and lipid domains

    The anionic lipids such as PS and PI are important components of the inner leaflet of erythrocyte membrane. To reveal the effect of anionic lipids on the interactions between the spectrins and erythrocyte membrane, we substituted some of the saturated lipids with the anionic lipids PIP2 or POPS in the inner leaflet of the ternary lipid bilayers of DPPC/DIPC/CHOL and DPPC/DIPE/CHOL, but kept the compositions of the outer leaflets unchanged. We examined the interactions between the spectrins and the lipid bilayers with the inner leaflets composed of DPPC/DIPC/CHOL/PIP2,DPPC/DIPC/CHOL/POPS, DPPC/DIPE/CHOL/PIP2, and DPPC/DIPE/CHOL/POPS, respectively, as presented in Table 1 and the simulation results are shown in Fig.3.

    In the lipid bilayers with the inner leaflets of DPPC/DIPC/CHOL/PIP2 and of DPPC/DIPC/CHOL/POPS,after introducing the anionic lipids which mainly distribute in the Lo domains,the spectrins are attracted by the anionic lipids and adsorbed onto the surface of the PIP2-or POPS-containing Lo domains as shown in Figs. 3(a) and 3(b). We also calculated the numbers of contacts between the headgroups of the anionic lipids and the amino acids on the spectrin surface as shown in Fig.4. We find that the main binding sites of the anionic lipids on the spectrin surface are the positively charged(Lys, Arg) or polar (Thr, Ser) amino acids. Therefore, the binding of the spectrins on the surface of the Lo domains in Figs. 3(a) and 3(b) is dominated by the electrostatic interactions between the anionic lipids(PIP2 or POPS)and the positively charged or polar amino acids on the spectrin surface since the interactions between the spectrins and PC lipids is very weak as shown in Fig. 2(a). Additionally, as shown in Fig.4,PIP2 and POPS bind to the spectrins at different sites,indicating they show specific interactions with the spectrins,respectively.

    Fig. 4. The numbers of contacts between the headgroups of the anionic lipids (PIP2 and POPS) and the amino acids on the surface of the spectrins which are adsorbed on the surface of the lipid bilayers composed of DPPC/DIPC/CHOL/PIP2 and DPPC/DIPC/CHOL/POPS, respectively.Here,the last 1-μs trajectories of the MD simulations are used for analysis.

    However,in the bilayers with inner leaflets composed of DPPC, DIPE, CHOL, and PIP2 (or POPS) (Figs. 3(c) and 3(d)),the interaction between the spectrins and DIPE is strong and the molar ratio of DIPE is much larger than that of PIP2 or POPS, the interaction between the spectrins and the lipid bilayers is dominated by the attraction between the spectrins and DIPE.Therefore,the spectrins are still adsorbed on the surface of Ld domain just as in the neutral DPPC/DIPE/CHOL bilayer.Furthermore, because of the electrostatic interaction between the anionic lipids and the positively charged or polar amino acids on the spectrin surface, the distribution of the anionic lipids is disturbed by the spectrins. PIP2 lipids with a lower molar fraction mainly distribute in the Ld domain while POPS lipids with a higher molar fraction mostly distribute in the Lo domain and are partially partitioned into the Ld domain by the attraction of the spectrins.

    3.3. Effect of mutation on the interactions of spectrins with lipid domains

    It has been shown that the mutation of the spectrins may result in the erythrocyte deformation which associates with some blood diseases in previous work.[15–22]Here, we also considered the effect of mutation on the interactions between spectrins and erythrocyte membrane. We considered the mutations of W263R and W300R in theβ17 spectrin(Fig.1(c)),which can cause the blood disease of HE.[22]We examined the interactions of the mutant spectrins with the lipid bilayers composed of DPPE/DIPC/CHOL and the simulation results are shown in Fig.5. We find that, compared with the case of the wild-type spectrins(Fig.2(c)),the contacts between theαspectrins and the Lo domains significantly increases.Most importantly, the mutant spectrins become more strongly curved and the lipid bilayer is strongly bent by the spectrins correspondingly.

    Fig.5.Top views and lateral view of the final structure of the mutant spectrinbound lipid bilayer composed of DPPE/DIPC/CHOL,with the spectrins,and the numbers of contacts between the spectrins and the lipid headgroups.Here,the mutation sites are W263R and W300R on the β17 spectrin.

    We compared the RMSDs of the wild-type and the mutant spectrins(Fig.6). The results show that the RMSD of the mutant spectrins is larger than that of the wild-type spectrins,especially in the later stage of the simulation, indicating that the mutation may lead to the decrease of the structural stability of the spectrins. Because tryptophan is hydrophobic while arginine is positively-charged at physiological conditions,they have completely different properties and interact differently with the other amino acids in the spectrins. Thus, when the mutations of W263R and W300R occurs, the interaction between theαandβspectrins will be greatly influenced and the spectrin tetramer becomes more unstable compared to the wild-type spectrins. In the relevant experiment,[22]it has also found that these two mutations can lead to the decrease of the affinity of theαandβspectrins. Therefore, our simulation results agree very well with the experimental results.

    Furthermore,we also compared the curvature of the lipid bilayers interacting with the wild-type and the mutant spectrins in Fig.7. Due to its unstable structure,the mutant spectrins can fluctuate strongly on the surface of the lipid bilayer and result in the stronger curvature of the bilayer. Therefore,the maximal curvature of the bilayer bound with the mutant spectrins is obviously larger than that of the bilayer bound with the wild-type spectrins as shown in Fig. 7. This result demonstrates that, in the blood diseases such as HE or HPP,the deformation of erythrocyte is possibly induced by the mutant spectrins with an unstable structure.[22]

    Fig.6. The RMSDs of the wild-type spectrins and the mutant spectrins with the mutation sites of W263R and W300R adsorbed on the surface of the lipid bilayer composed of DPPE/DIPC/CHOL.

    Fig. 7. The maximal curvatures of the DPPE/DIPC/CHOL lipid bilayer bound by the wild-type and the mutant spectrins.

    4. Conclusion

    In this work,using coarse-grained MD simulation,we examined the interactions between the wild-type or mutant spectrins and lipid nanodomains in the multicomponent lipid bilayers. The results show that the spectrins exhibit a stronger interaction with PE lipids than with PC lipids, which agrees well with the experimental results. Additionally, the spectrins prefer to contact with the Ld domains in the ternary neutral lipid bilayers containing unsaturated DIPE. The anionic lipids can enhance the interactions between the spectrins and the lipid bilayers through the electrostatic interaction between the positively charged or polar amino acids on the surface of the protein and the anionic headgroups of the lipids, and the binding of the spectrins obviously influences the distribution of the anionic lipids in the lipid bilayers. Furthermore, we also considered the effect of mutation on the interactions between the spectrins and the lipid bilayer. It is found that the spectrin mutation reduces the structural stability of the spectrins and further leads to the strong deformation of the lipid bilayer,which could be attributed to the causes of some blood diseases. Our results may yield some theoretical insights into understanding the erythrocyte structure and the treatment of some spectrin-associated blood diseases.

    香蕉国产在线看| 国产老妇伦熟女老妇高清| 中文字幕制服av| 亚洲,欧美精品.| 大话2 男鬼变身卡| 99久久精品国产亚洲精品| 黄片播放在线免费| 日韩一区二区三区影片| 日韩视频在线欧美| 精品国产乱码久久久久久男人| 99热国产这里只有精品6| 成人18禁高潮啪啪吃奶动态图| 中文字幕精品免费在线观看视频| 国产在线一区二区三区精| 亚洲久久久国产精品| 这个男人来自地球电影免费观看| 男女国产视频网站| 麻豆乱淫一区二区| 欧美黄色片欧美黄色片| 在线观看免费日韩欧美大片| 丰满迷人的少妇在线观看| 美女国产高潮福利片在线看| 搡老岳熟女国产| 国产男女超爽视频在线观看| 女人精品久久久久毛片| av不卡在线播放| 十八禁高潮呻吟视频| 日韩熟女老妇一区二区性免费视频| 成人黄色视频免费在线看| 久久ye,这里只有精品| 日本欧美国产在线视频| av欧美777| 高清黄色对白视频在线免费看| 18禁裸乳无遮挡动漫免费视频| 亚洲av片天天在线观看| 精品人妻一区二区三区麻豆| 少妇裸体淫交视频免费看高清 | 亚洲国产av影院在线观看| 国产av国产精品国产| 老汉色∧v一级毛片| 亚洲国产精品999| 视频在线观看一区二区三区| 国产日韩一区二区三区精品不卡| 妹子高潮喷水视频| av有码第一页| 99久久精品国产亚洲精品| 真人做人爱边吃奶动态| 亚洲 国产 在线| 国产在线免费精品| 操美女的视频在线观看| 国产免费一区二区三区四区乱码| 国产一区有黄有色的免费视频| 99精国产麻豆久久婷婷| 久久久国产一区二区| 免费日韩欧美在线观看| 两个人免费观看高清视频| cao死你这个sao货| 熟女少妇亚洲综合色aaa.| 色婷婷久久久亚洲欧美| 久久亚洲精品不卡| 大香蕉久久成人网| 国产成人免费无遮挡视频| 久久久久久久大尺度免费视频| 91精品三级在线观看| 午夜av观看不卡| 又大又黄又爽视频免费| 在线观看免费高清a一片| 91字幕亚洲| 国产一区二区激情短视频 | 亚洲 欧美一区二区三区| 美女视频免费永久观看网站| 午夜福利在线免费观看网站| 大码成人一级视频| 亚洲欧美日韩另类电影网站| 日韩,欧美,国产一区二区三区| 又大又爽又粗| 丁香六月欧美| 免费黄频网站在线观看国产| 电影成人av| 男女边吃奶边做爰视频| a级毛片黄视频| 在线观看免费视频网站a站| 日韩人妻精品一区2区三区| 国产精品三级大全| 免费在线观看影片大全网站 | 男女国产视频网站| 国产高清视频在线播放一区 | 在线观看免费午夜福利视频| 丁香六月欧美| 免费少妇av软件| 99国产精品一区二区蜜桃av | 亚洲国产精品999| 五月开心婷婷网| 黄色毛片三级朝国网站| 亚洲成人手机| 秋霞在线观看毛片| 亚洲av成人不卡在线观看播放网 | 亚洲精品一二三| 曰老女人黄片| 一级a爱视频在线免费观看| 久久久久久久大尺度免费视频| 国产又色又爽无遮挡免| 国产三级黄色录像| 欧美精品一区二区免费开放| 亚洲av电影在线观看一区二区三区| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲欧美中文字幕日韩二区| 男人舔女人的私密视频| 国产xxxxx性猛交| 亚洲中文av在线| 日韩av不卡免费在线播放| 欧美少妇被猛烈插入视频| 欧美日韩国产mv在线观看视频| 97人妻天天添夜夜摸| 午夜精品国产一区二区电影| av视频免费观看在线观看| cao死你这个sao货| 欧美日韩亚洲国产一区二区在线观看 | 欧美人与善性xxx| 黄频高清免费视频| 婷婷色综合大香蕉| 亚洲人成电影免费在线| 久久ye,这里只有精品| 九色亚洲精品在线播放| 日本91视频免费播放| 国产欧美日韩精品亚洲av| 欧美日韩综合久久久久久| 亚洲精品自拍成人| 99久久精品国产亚洲精品| 日本av手机在线免费观看| 亚洲精品第二区| tube8黄色片| www.999成人在线观看| 男女下面插进去视频免费观看| 亚洲激情五月婷婷啪啪| 中文字幕人妻丝袜制服| 亚洲激情五月婷婷啪啪| 亚洲国产av新网站| 两个人看的免费小视频| 国产精品久久久久久精品电影小说| 考比视频在线观看| 国产成人av教育| 最近手机中文字幕大全| 国产精品久久久人人做人人爽| 亚洲五月婷婷丁香| 99久久综合免费| 精品亚洲成a人片在线观看| 欧美在线一区亚洲| 久久午夜综合久久蜜桃| 丰满饥渴人妻一区二区三| 九草在线视频观看| 成年女人毛片免费观看观看9 | 久久久久久久大尺度免费视频| 久久久国产精品麻豆| 欧美黄色淫秽网站| 国产精品一区二区精品视频观看| 国产精品国产av在线观看| 国产成人影院久久av| 高潮久久久久久久久久久不卡| av线在线观看网站| 天天操日日干夜夜撸| 国产亚洲欧美精品永久| 一级黄色大片毛片| 人人妻人人澡人人爽人人夜夜| 少妇人妻久久综合中文| 欧美精品亚洲一区二区| 免费高清在线观看日韩| 搡老乐熟女国产| 伊人久久大香线蕉亚洲五| 亚洲欧美一区二区三区黑人| 国产99久久九九免费精品| 亚洲少妇的诱惑av| 成年人黄色毛片网站| 精品国产国语对白av| 巨乳人妻的诱惑在线观看| 高清黄色对白视频在线免费看| 欧美少妇被猛烈插入视频| 在线亚洲精品国产二区图片欧美| 久久免费观看电影| 好男人视频免费观看在线| 成年动漫av网址| 亚洲熟女精品中文字幕| 精品国产国语对白av| 成人午夜精彩视频在线观看| av片东京热男人的天堂| 国产黄色视频一区二区在线观看| 国产亚洲欧美在线一区二区| 国产精品三级大全| videos熟女内射| 男女边摸边吃奶| 日日摸夜夜添夜夜爱| 免费少妇av软件| 青草久久国产| netflix在线观看网站| 一区福利在线观看| 亚洲免费av在线视频| 97人妻天天添夜夜摸| 午夜福利,免费看| bbb黄色大片| 色婷婷久久久亚洲欧美| 美女高潮到喷水免费观看| 少妇的丰满在线观看| 国产亚洲一区二区精品| 大话2 男鬼变身卡| netflix在线观看网站| 国产精品国产av在线观看| 日本wwww免费看| 丝袜脚勾引网站| 岛国毛片在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 性色av一级| 欧美日韩综合久久久久久| 久久久精品94久久精品| 少妇精品久久久久久久| 精品国产乱码久久久久久小说| 亚洲午夜精品一区,二区,三区| 看免费成人av毛片| 欧美精品av麻豆av| 欧美av亚洲av综合av国产av| 亚洲av成人精品一二三区| 午夜激情av网站| 90打野战视频偷拍视频| 亚洲精品久久午夜乱码| 国产又色又爽无遮挡免| 国产亚洲欧美精品永久| 亚洲成人免费av在线播放| 免费人妻精品一区二区三区视频| 男人爽女人下面视频在线观看| 在线精品无人区一区二区三| 久久久精品94久久精品| 久久国产精品男人的天堂亚洲| 国产精品成人在线| 丝袜脚勾引网站| 90打野战视频偷拍视频| 久久青草综合色| 精品视频人人做人人爽| 日韩大片免费观看网站| 韩国精品一区二区三区| 这个男人来自地球电影免费观看| 亚洲午夜精品一区,二区,三区| 女性生殖器流出的白浆| 我的亚洲天堂| 80岁老熟妇乱子伦牲交| 亚洲综合色网址| 91字幕亚洲| 日韩大码丰满熟妇| 9热在线视频观看99| 免费高清在线观看日韩| 中文字幕av电影在线播放| 赤兔流量卡办理| 色婷婷av一区二区三区视频| 国产成人欧美| 免费人妻精品一区二区三区视频| 亚洲欧美一区二区三区国产| 99国产精品一区二区蜜桃av | 高清黄色对白视频在线免费看| 日日爽夜夜爽网站| 亚洲久久久国产精品| 2018国产大陆天天弄谢| 国产在线一区二区三区精| 嫁个100分男人电影在线观看 | 久久性视频一级片| 一二三四在线观看免费中文在| 亚洲中文av在线| 久久影院123| 亚洲欧美色中文字幕在线| 国产欧美日韩综合在线一区二区| 老汉色∧v一级毛片| 美女大奶头黄色视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲久久久国产精品| 国产成人啪精品午夜网站| 精品亚洲成a人片在线观看| 欧美人与性动交α欧美精品济南到| 欧美日韩精品网址| 最近最新中文字幕大全免费视频 | 纵有疾风起免费观看全集完整版| 又黄又粗又硬又大视频| 免费高清在线观看日韩| 中文字幕高清在线视频| 亚洲国产日韩一区二区| 91老司机精品| 国产一卡二卡三卡精品| 国产视频首页在线观看| 国产真人三级小视频在线观看| 操出白浆在线播放| 免费在线观看影片大全网站 | 欧美精品啪啪一区二区三区 | 一区二区三区精品91| 中文乱码字字幕精品一区二区三区| 久久久国产精品麻豆| 色视频在线一区二区三区| 成年动漫av网址| 欧美黄色淫秽网站| 丁香六月天网| 最新在线观看一区二区三区 | 中国美女看黄片| 亚洲精品国产av成人精品| 国产亚洲欧美在线一区二区| 高清不卡的av网站| www.av在线官网国产| 欧美大码av| 日韩精品免费视频一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 十八禁高潮呻吟视频| 最黄视频免费看| 久久久精品区二区三区| 十八禁人妻一区二区| 欧美日韩av久久| 男女之事视频高清在线观看 | 老汉色av国产亚洲站长工具| 我的亚洲天堂| 精品少妇黑人巨大在线播放| 黄色怎么调成土黄色| 熟女av电影| 久久久精品94久久精品| 另类亚洲欧美激情| 亚洲欧洲日产国产| 欧美黑人精品巨大| 国产亚洲午夜精品一区二区久久| 国产主播在线观看一区二区 | 午夜福利,免费看| 激情五月婷婷亚洲| 另类亚洲欧美激情| av天堂久久9| 两个人免费观看高清视频| 又粗又硬又长又爽又黄的视频| 在线观看人妻少妇| 最新在线观看一区二区三区 | 日本色播在线视频| 日韩一区二区三区影片| 国产黄频视频在线观看| 激情视频va一区二区三区| 伊人亚洲综合成人网| 国产精品熟女久久久久浪| 欧美激情高清一区二区三区| av电影中文网址| 国产高清不卡午夜福利| av电影中文网址| 成人亚洲欧美一区二区av| 高清不卡的av网站| 国产免费现黄频在线看| 精品国产国语对白av| 90打野战视频偷拍视频| 午夜91福利影院| av有码第一页| 好男人视频免费观看在线| 国产亚洲欧美在线一区二区| 99久久综合免费| 女性生殖器流出的白浆| 久久精品国产综合久久久| 色94色欧美一区二区| 99国产精品99久久久久| 亚洲人成77777在线视频| 久久国产亚洲av麻豆专区| 捣出白浆h1v1| 国产1区2区3区精品| 亚洲精品美女久久av网站| 高潮久久久久久久久久久不卡| 深夜精品福利| 亚洲 欧美一区二区三区| 十八禁网站网址无遮挡| 亚洲色图综合在线观看| 人人妻人人澡人人爽人人夜夜| 91精品伊人久久大香线蕉| 日韩熟女老妇一区二区性免费视频| 中文字幕av电影在线播放| 国产在线免费精品| 王馨瑶露胸无遮挡在线观看| 久久中文字幕一级| 五月天丁香电影| 亚洲欧美一区二区三区黑人| a 毛片基地| 久久久久久久久久久久大奶| 又紧又爽又黄一区二区| 免费人妻精品一区二区三区视频| 美女视频免费永久观看网站| 人人妻人人澡人人看| 亚洲图色成人| 成人国语在线视频| 久久精品熟女亚洲av麻豆精品| 一级片'在线观看视频| 久久精品熟女亚洲av麻豆精品| 亚洲国产中文字幕在线视频| 国产欧美日韩一区二区三 | av在线播放精品| 国产三级黄色录像| 女人爽到高潮嗷嗷叫在线视频| 香蕉丝袜av| 热99国产精品久久久久久7| 国产免费福利视频在线观看| 老司机在亚洲福利影院| 天天躁夜夜躁狠狠久久av| 午夜激情av网站| 免费看十八禁软件| 欧美亚洲日本最大视频资源| 国产成人av教育| 十八禁高潮呻吟视频| 国产成人欧美在线观看 | 丰满迷人的少妇在线观看| 人人妻人人爽人人添夜夜欢视频| 日日摸夜夜添夜夜爱| 人人妻人人澡人人爽人人夜夜| 咕卡用的链子| cao死你这个sao货| 国产精品国产三级国产专区5o| 十八禁人妻一区二区| 免费看av在线观看网站| 18禁黄网站禁片午夜丰满| 午夜免费成人在线视频| 成人三级做爰电影| 精品人妻在线不人妻| 中国国产av一级| 99re6热这里在线精品视频| 国产99久久九九免费精品| videosex国产| 欧美性长视频在线观看| 中国美女看黄片| 在线观看www视频免费| 亚洲天堂av无毛| 在线亚洲精品国产二区图片欧美| 七月丁香在线播放| 久久精品熟女亚洲av麻豆精品| 国产片内射在线| 人妻一区二区av| 啦啦啦中文免费视频观看日本| 成人亚洲精品一区在线观看| 日韩免费高清中文字幕av| 永久免费av网站大全| 亚洲精品成人av观看孕妇| 制服诱惑二区| 亚洲图色成人| 又粗又硬又长又爽又黄的视频| 少妇被粗大的猛进出69影院| 国产成人av激情在线播放| 国产亚洲欧美在线一区二区| 国产精品久久久久成人av| 老鸭窝网址在线观看| 国产精品一二三区在线看| 亚洲一卡2卡3卡4卡5卡精品中文| av在线app专区| 中文精品一卡2卡3卡4更新| 午夜免费观看性视频| 中国国产av一级| 国产欧美亚洲国产| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品久久久久久| 国产黄色视频一区二区在线观看| 日本色播在线视频| 亚洲成人免费av在线播放| 亚洲国产中文字幕在线视频| 99国产精品免费福利视频| 中文字幕高清在线视频| 51午夜福利影视在线观看| 狂野欧美激情性bbbbbb| 久9热在线精品视频| 久久精品国产a三级三级三级| av又黄又爽大尺度在线免费看| 亚洲人成网站在线观看播放| 亚洲精品av麻豆狂野| 丝袜美腿诱惑在线| 侵犯人妻中文字幕一二三四区| 十分钟在线观看高清视频www| 成年av动漫网址| e午夜精品久久久久久久| avwww免费| 欧美少妇被猛烈插入视频| 精品人妻熟女毛片av久久网站| 精品国产国语对白av| 大型av网站在线播放| av天堂久久9| 黄色视频不卡| 久久精品亚洲熟妇少妇任你| 十八禁人妻一区二区| 午夜精品国产一区二区电影| 中文字幕色久视频| 国产伦理片在线播放av一区| 午夜福利视频精品| 欧美在线一区亚洲| 亚洲国产欧美日韩在线播放| 少妇的丰满在线观看| 国产欧美日韩一区二区三 | 热99国产精品久久久久久7| 婷婷色麻豆天堂久久| 久久久久久久大尺度免费视频| 欧美黄色淫秽网站| 欧美老熟妇乱子伦牲交| 精品卡一卡二卡四卡免费| 永久免费av网站大全| 免费看av在线观看网站| 一级毛片 在线播放| 国产精品久久久人人做人人爽| 日日夜夜操网爽| 性色av乱码一区二区三区2| 一二三四在线观看免费中文在| av福利片在线| 久久狼人影院| av国产精品久久久久影院| 性色av乱码一区二区三区2| 亚洲熟女毛片儿| 国产国语露脸激情在线看| 国产免费又黄又爽又色| 校园人妻丝袜中文字幕| 欧美xxⅹ黑人| 精品国产超薄肉色丝袜足j| 99热网站在线观看| 国产不卡av网站在线观看| 国产成人精品久久二区二区免费| 日韩人妻精品一区2区三区| 中文字幕制服av| 桃花免费在线播放| 亚洲第一青青草原| 50天的宝宝边吃奶边哭怎么回事| 国产在线视频一区二区| 丝袜美足系列| 视频在线观看一区二区三区| 久久久久久久大尺度免费视频| av福利片在线| 老熟女久久久| √禁漫天堂资源中文www| 国产精品九九99| a级毛片黄视频| 亚洲一卡2卡3卡4卡5卡精品中文| 美女午夜性视频免费| 色视频在线一区二区三区| 狠狠婷婷综合久久久久久88av| 亚洲中文字幕日韩| 99久久精品国产亚洲精品| 蜜桃在线观看..| av国产久精品久网站免费入址| 欧美在线一区亚洲| 亚洲专区中文字幕在线| 免费av中文字幕在线| 成人亚洲精品一区在线观看| 久久国产精品影院| 国产一区二区在线观看av| 午夜免费男女啪啪视频观看| 大片免费播放器 马上看| 日韩一本色道免费dvd| www.自偷自拍.com| 久久久久精品国产欧美久久久 | 欧美在线一区亚洲| 18禁观看日本| 男女高潮啪啪啪动态图| 黑丝袜美女国产一区| 亚洲三区欧美一区| 亚洲av男天堂| 亚洲第一青青草原| 老司机深夜福利视频在线观看 | 亚洲 欧美一区二区三区| 国产成人精品久久二区二区91| 欧美成人精品欧美一级黄| 熟女少妇亚洲综合色aaa.| 久久鲁丝午夜福利片| 下体分泌物呈黄色| 国产成人一区二区三区免费视频网站 | 成人亚洲欧美一区二区av| 国产视频一区二区在线看| 99久久精品国产亚洲精品| 亚洲第一青青草原| 2021少妇久久久久久久久久久| 日本色播在线视频| 午夜福利乱码中文字幕| 色播在线永久视频| 在线观看免费日韩欧美大片| 国产日韩欧美在线精品| 免费在线观看日本一区| 伊人亚洲综合成人网| av线在线观看网站| 国产一区二区 视频在线| 成年美女黄网站色视频大全免费| 欧美日韩亚洲高清精品| 欧美国产精品一级二级三级| 欧美日韩综合久久久久久| 啦啦啦啦在线视频资源| 中文精品一卡2卡3卡4更新| 国产一区二区激情短视频 | 国产91精品成人一区二区三区 | 国产成人a∨麻豆精品| 亚洲av国产av综合av卡| 一个人免费看片子| 国产av国产精品国产| 男女边摸边吃奶| videos熟女内射| 高潮久久久久久久久久久不卡| cao死你这个sao货| 亚洲伊人色综图| 欧美久久黑人一区二区| 高清黄色对白视频在线免费看| 自线自在国产av| 亚洲av电影在线观看一区二区三区| 精品免费久久久久久久清纯 | 国产男人的电影天堂91| 看十八女毛片水多多多| 秋霞在线观看毛片| 成年女人毛片免费观看观看9 | 亚洲国产精品一区三区| 久久精品国产亚洲av高清一级| 国产爽快片一区二区三区| 制服人妻中文乱码| 性色av乱码一区二区三区2| 狠狠精品人妻久久久久久综合| 男女高潮啪啪啪动态图| 国产伦人伦偷精品视频| 999久久久国产精品视频| 在线观看国产h片| 亚洲国产精品一区二区三区在线| 丝瓜视频免费看黄片| 制服人妻中文乱码| 久久99一区二区三区| 亚洲伊人色综图| 久久久久久久精品精品| 一级,二级,三级黄色视频| 欧美精品高潮呻吟av久久| 国产成人a∨麻豆精品|