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

    Diffusion of nucleotide excision repair protein XPA along DNA by coarse-grained molecular simulations?

    2021-10-28 07:18:44WeiweiZhang張偉偉andJianZhang張建
    Chinese Physics B 2021年10期
    關鍵詞:張建

    Weiwei Zhang(張偉偉) and Jian Zhang(張建)

    National Laboratory of Solid State Microstructures,School of Physics,Collaborative Innovation Center of Advanced Microstructures,Nanjing University,Nanjing 210093,China

    Keywords: nucleotide excision repair,XPA,one-dimensional diffusion along DNA,molecular simulation

    1. Introduction

    Nucleotide excision repair(NER)is the central DNA repair pathway widely observed in eubacteria, eukaryotes and archaea. It is responsible for the removal of DNA lesions caused by UV irradiation,chemical mutagens,and other stress damages.[1–5]The abnormalities of NER can lead to the xeroderma pigmentosum(XP)disease,which is featured by sensitivity to sunlight and high risk of developing skin cancers.[6,7]XPA (xeroderma pigmentosum complementation group A) is a key protein in NER,and mutation of XPA gene can lead to severe disease symptoms.[6,8]It has been showed that the XPA may play multiple roles in NER.[9–12]For example,in addition to preferentially binding to various types of DNA lesions and distorted substrates,[10,13]XPA has also been found to interact with several NER related proteins during the damage recognition and verification steps,acting as a scaffold.[11,14]

    The human XPA has the length of 273 amino acids and is composed of three domains, including the intrinsically disordered N-terminal region(N-region),the central DNA-binding domain(DBD),and the intrinsically disordered C-terminal region (C-region) (see Fig. 1). The flexibility of the two disordered regions were suggested to promote the binding of XPA to several partners,[15–19]while the central globular DNAbinding domain was found to bind with DNA and interact with numerous NER proteins as well.[20]

    Due to its biological importance in NER, XPA has received extensive attention.[7,21]Although previous investigations mostly focused on the interactions of the XPA at the DNA lesion site during the damage recognition and verification steps, recent experimental study showed that detecting the DNA lesions may also involve the one-dimensional diffusion of the XPA at the non-damage site of DNA.[22]Therefore,more detailed investigations for the molecular dynamics of the XPA diffusion along the non-damage site of DNA are important for deeper understanding the biological function of XPA in the NER.

    Although all-atom detailed molecular dynamics (MD)simulation has been widely used to directly simulate the details of the inter-molecular interactions involved in the protein-DNA recognitions, it is still extremely timeconsuming to study the whole diffusion process of proteins along DNA.[23–26]Recently, the combination of the atomic-interaction-based coarse-grained model(AICG2+)for proteins[27,28]and the three-site-per-nucleotide (3SPN.2C)model for DNA[29]has shown great success in studying the dynamical and structural features of the protein-DNA interactions.[30–33]Here,based on these models,we performed extensive coarse-grained molecular simulations to study the diffusion dynamics of XPA along the DNA.From our simulations,we found that XPA diffuses along DNA with the DNAbinding residues similar to that identified from crystal structure. It was also found that the diffusion of XPA along DNA is achieved through combination of one-dimensional and threedimensional mechanisms. At low salt concentrations,the onedimensional diffusion with strong rotational coupling is the dominant mechanism. At higher salt concentrations, the diffusion by three-dimensional mechanism becomes more probable. Meanwhile,the one-dimensional diffusion of XPA along DNA displays sub-diffusive behaviour at wide range of salt concentrations. In addition, the simulation showed that the XPA becomes more extended in the presence of DNA or at high salt concentration,suggesting that both DNA binding and increasing salt concentration may promote the binding of other target proteins to XPA by increasing the exposure extent of the binding sites.

    2. Model and method

    2.1. Coarse-grained model of protein-DNA system.

    Throughout this work, to model protein (XPA), the structure-based coarse-grained model was used, where each amino acid in protein was coarse-grained to a single bead located at Cαposition. The interactions between protein beads include AICG2+ potentials and Debye–H¨uckel-type electrostatic interactions. The AICG2+ energy function is given byVAICG2+=Vlocal+VGo+Vexv, whereVlocal,VGO, andVexvare local potentials, non-local structure-based potentials, and excluded volume potentials respectively. The details of these terms can be found in Refs. [28,34,35], and the similar energy function has also been successfully used in the studies of protein folding,enzyme catalysis,chromatin remodeler and nucleosome dynamics in recent works.[27,36–40]Note that the structure-based interactions were omitted for the disordered regions in XPA, and effects of salt concentration are implicitly modeled in the Debye–H¨uckel-type interactions with the details displayed as follows:

    whereλDis the Debye length,ε(T,C) is the dielectric constant,andIis ionic strength.

    To model dsDNA,the 3SPN.2C model,[29]developed by de Pablo’s group, was used. In this model each nucleotide is simplified to three CG particles, representing the phosphate group (P), sugar (S), and base (B). The interactions among these particles are composed of structure-based local potentials,base pairing,intrastrand base-stacking,interstrand crossstacking, excluded volume effects, and Debye–H¨uckel-type electrostatic interactions. The partial charge of the phosphate beads was set as?0.6e to consider the effect of local counterion condensation effect.[31,41]The 3SPN.2C model model is able to capture several properties of dsDNA and has been widely used to study the DNA-protein interactions.

    For the protein-DNA interactions,only the excluded volume effects and the electrostatic interactions were considered. As the electrostatic interaction plays a dominant role in protein-DNA interactions, we paid special attention to this term in the study of protein-DNA binding. For the structured motif in proteins, the partial charges of the CG beads were parametrized by the RESPAC method,[42]which can provide more accurate description of the electrostatic potential. In describing the electrostatic interactions between the protein and DNA,the partial charge of the phosphate bead was set as?1.0,following the work in Ref.[31].

    2.2. Simulation details

    The model for XPA was built based on the available crystal structure(PDB entry:6j44[8])with missing residues added by modeller9.23.[43]The sequence and structure of XPA are shown in Figs. 1(a), 1(b), and 1(c). For DNA, two random sequences with different lengths were prepared: R100, R150(the sequences and the structures are shown in Figs.1(d)and 1(e)). In our work, the reference structures for DNA were built with the 3DNA package.[44]For the protein-DNA system, the initial structures for DNA and protein in all the simulations were separated and placed randomly. During production runs,the center of mass(COM)of DNA was restrained to the original position and COM distance between XPA(DBD)and DNA was restrained to be smaller than 250 °A.The ionic strength was set to 50, 150, 250, 350, and 450 mM. At each ionic strength,24 simulations were carried out independently,each lasting for 5×108MD steps for XPA-R100 system and over 4×108MD steps for XPA-R150 system,respectively. To characterize the structure features of XPA in the absence of DNA,we also performed several simulations with cumulative length of approximately 1×109MD steps. For all the simulations,we recorded the coordinates of XPA and DNA at every 1×104MD steps. The first 2×107frames were omitted in the calculations of distributions during the posterior analysis.

    All the CG simulations were performed by the CafeMol3.2 package.[34]The simulations were conducted by Langevin dynamics with friction coefficientγ=0.02 at temperatureTof 300 K.The structures were visualized by using pymol.[45]

    2.3. Analysis

    To characterize the search modes of XPA on DNA, we defined a DNA recognition motif consisting of the residues 140, 141, 142, 151, 179, 207, 211,[46]which forms contacts with DNA in the crystal structure. Following the approach used in Ref. [41], the diffusion is considered as sliding if at least five of these residues maintain staying within 15 °A from the DNA.Correspondingly,the diffusion is considered as three-dimensional diffusion if all of the seven residues in the recognition motif leaves the DNA surface by at least 25 °A.In comparison,the protein is considered to perform the hopping dynamics if these residues are 15 °A–25 °A away from DNA.Based on above criteria, we classified both sliding and hopping as one-dimensional diffusion mechanism following the work in Ref[41].

    To describe the location of DBD searching site on DNA during one-dimensional diffusion,we used the mean position of the nucleotides that closest to the recognition motif(RM):

    In the above definition,tis simulation time,Lis 3.4 °A (one bp along the strand), RM is the recognition motif, composed of the residues 140,141,142,151,179,207,and 211,NRMis the number of residues in recognition motif,bp(j)is the base pair index of particlejin DNA;δkj=1,ifjis the DNA bead closest to the protein beadk;otherwiseδkj=0.[31]

    The protein bead and DNA bead form a contact if their distance is less than 7 °A.Additionally,we defined several local coordinate axes, following the procedure in Ref. [47], to characterize the rotational and translational motions of XPA around DNA.

    Fig.1. Sequence and structure of the XPA and DNA studied in this work. (a)Domain architecture of XPA.(b)The sequence of XPA protein with each domain shadowed by the same color as in panel(a). The secondary structure features of DBD domain is also shown. (c)The available crystal sturcture of DBD is shown in cartoon representation and the color scheme for secondary structures is the same as in panel(b).The disorderd regions are represented by dotted lines with the same color as in panel(a). The key DNA-binding residues(residue 140,141,142,151,179,207,211)are shown in ball-and-stick representation and colored differently. The sequences(d)and structures(e)of DNA in this study.

    3. Results

    3.1. XPA binding and diffusion along DNA

    Firstly, to quantify the structural properties of nonspecific interaction interface between XPA and DNA, based on two random DNA sequences with different lengths(R100,R150), we performed MD simulations for the XPA binding to the nonspecific dsDNA site at various salt concentrations and calculated the probabilities of contact formation between the DNA and each residue of the XPA at various salt concentrations. As shown in Fig. 2, there is a striking similarity in the landscape of contact probability between R100 and R150,indicating the conserved binding interface of XPA regardless of properties of DNA.Additionally,for both R100 and R150,one can see that, besides the DBD, residues in the N-region also have high contacting probabilities with DNA,which suggests that the N-region can be involved in the XPA–DNA interactions during the searching process of damage sites. Meanwhile, the results clearly show that increasing salt concentrations tends to decrease the contact probabilities, but without altering the peak locations of the contact probability distributions, which may suggest that the XPA uses the same set of residues for DNA binding at wide range of salt concentrations.Particularly, the residues with high contact probabilities have large overlap with the DNA-binding residues observed in crystal structure of the XPA–DNA complex,[46]demonstrating that the DNA binding surface in the DBD is conserved in the specific and non-specific binding. Such behaviour of maintaining binding interface was also observed in other DNA-binding proteins,[47]and was suggested to facilitate the search of the cognate site on DNA.[48,49]

    Fig. 2. Contact probabilities of the XPA residues with DNA at various salt concentrations. Residues in the DBD domain are shaded by pink backgound.The key DNA-binding residues in the crystal structure are highlighted by magenta circles. The result with R100(R150)is shown in the upper(lower)panel and a representative structure with DBD bound to R100 is also displayed.

    Fig.3. (a)and(b)Time series of the distance between the COM of N-region and DNA surface(green)and that between the COM of DBD and DNA surface (red) at 250 mM with sequence R100 (a) and R150 (b). When calculating the distance between the COM of N-terminal and DNA surface, only the residues ranging from 30 to 50, which have high probabilities of forming contacts with DNA, were considered. (c) and (d) Probabilities of onedimensional diffusion(green)and three-dimensional diffusion(pink)mechanism at different salt concentrations with sequence R100(c)and R150(d).

    Besides the properties of binding interface,we also characterized the binding and diffusion process of XPA along DNA. Considering the relatively higher probabilities of contact formation with DNA for the residues in the DBD and the N-region, we plotted the time series of the distances between the above two protein segments and the DNA surface to illustrate the binding dynamics. The centers of mass(COM)of the DBD and that of the residues 30–50 were used to represent the positions of the DBD and the N-region relative to DNA surface. From Figs.3(a)and 3(b), one can clearly see the attachment and detachment processes of the two XPA segments.Such attachment and detachment processes are consistent with experimental observation.[22]However, it is worthy to note that,because of the limited accessible time in our simulation,the repetition of attachment and detachment is only observed sporadically at relative high salt concentration (Fig. 3(a)).Additionally,inspection of trajectory of these two protein segments demonstrated that the attachment/detachment of the DBD and N-region are correlated, lacking of the anchoring and binding mechanism mediated by disordered region as typically observed in other DNA-binding proteins containing disordered regions.[50]After the simultaneous detachments of the N-region and the DBD,the whole XPA performs threedimensional diffusion and then rebinds to the DNA(Fig.3(a)).For many DNA binding proteins, rapidly finding their specific binding sites on DNA in response to external stimulus is vitally important. Thus, in addition to the random threedimensional diffusion driven by thermal diffusion,the protein also uses other searching mechanisms, such as sliding,[51–53]hopping or intersegmental transfer.[54]During sliding,the protein can tightly bind to DNA and diffuse along the DNA helix using the recognition region,[55]in which electrostatic interaction between protein and DNA plays dominant role. Whereas in the hopping diffusion process,the electrostatic interactions are much weaker, which allows the protein to diffuse more freely along the DNA surface and accelerates diffusion by skipping several DNA bases.[41]Unlike three-dimensional diffusion, in which the protein can diffuse far away from DNA,the protein remain staying close to the DNA during the sliding and hopping,which tends to significantly reduce the conformational space during target site search and therefore accelerate the searching of the target site in the DNA.[51]Previous studies showed that the combination of these different search modes is needed to achieve the maximal target-search efficiency.[54,56]Therefore, we also analyzed the probabilities of the one-dimensional diffusion mechanism and threedimensional diffusion mechanism (criteria about classifying different mechanisms can be found in Model and Method) at various salt concentrations for sequence R100 and R150, respectively (Figs. 3(c) and 3(d)). From these results, one can see that, increasing salt concentration can decrease the probability of one-dimensional searching mechanism and increase the three-dimensional searching mechanism for each DNA sequence. As shown in the Debye–Hc¨ukel model, the effect of salt concentration is mainly achieved by changing the length of Debye radius. With the increasing of salt concentrations,the screening effect of the implicit salt ions increases, leading to the decrease of electrostatic interactions. Such results again suggest the dominate role of electrostatic interaction in the XPA–DNA binding.

    Considering the important role of one-dimensional diffusion, we then analyzed the diffusion kinetics and calculated the mean square displacement(MSD)of the one-dimensional DBD translational motions along the DNA (MSD(?t) =〈(X(t+?t)?X(t))2〉, where〈〉represents the time average)as a function of time difference ?tat various salt concentrations, and results are shown in Fig.4. The fitting by a power law function MSD(?t) =K?tαwere also shown in the inset by log–log scale. From our results, one can clearly see that the MSD of XPA (DBD) increases upon increasing salt concentration during one-dimensional diffusion for both R100 and R150. The fitting of the diffusion kinetics by the above power law function for R100 and R150 both give theα<1.0,which suggests sub-diffusive motion of DBD along DNA(inset figure in Fig. 4) at wide range of salt concentrations and may facilitate the search of target DNA sequences as suggested by other works.[57–59]Furthermore,this kind of sub-diffusion may be related to the diverse strengths of DNA–XPA interactions,which may be ascribed to the variations of the DNA conformations and the fluctuations of the DNA-protein binding.This kind of connection is consistent with those in many other sub-diffusive dynamics,[59–69]and can be described with continuous time random walk(CTRW)model.[61,69]Note that the exponents of the sub-diffusive dynamics are almost invariant at different salt concentrations for each sequence. This observation suggests that the distributions of the diverse strengths are insensitive to the variations of salt concentrations, since the conformational fluctuations are intrinsic for the the DNAprotein system.This result bring us more insights on the DNAprotein interactions.

    Fig.4. Mean square displacement(MSD)of the one-dimensional XPA diffusion along the DNA(R100: upper panel,R150: lower panel)as a function of simulation time at different salt concentrations.The results of power function fitting α are also shown in the inset.

    Fig.5. Translational–rotational coupling in the one-dimensional diffusion of XPA along DNA.(a)Representative trajectory fragments of rotation(θ) and translation (y) of XPA along R100 (R150) at the salt concentrations of 50 mM, 150 mM, and 450 mM are shown in upper panel (lower panel). Representive structures with DBD bound to DNA at 50 mM and 450 mM are displayed. (b)The correlation plots between ?θ and ?y of XPA at 50 mM,150 mM,and 450 mM(R100: upper panel,R150: lower panel).

    For DNA-interacting proteins, several studies have suggested that diffusion along the DNA helical path is required to achieve efficient target search.[55,70]Meanwhile,the rotationcoupled sliding was proposed to be the general mechanism for the one-dimensional search of protein on DNA.[71]Therefore,we also investigated whether the one-dimensional diffusion of XPA is along helical pitch(rotation-coupled)or not(rotationuncoupled), and the results are shown in Fig. 5. To more appropriately describe this dynamical process,we defined several local axes following the Ref.[47]. The coordinates of the seven residues in the DNA-binding motif mentioned above were used to calculate the rotational angle and translational position along the DNA helical axis. In Fig. 5(a), we plotted the time series of rotational and translational motions of DBD along the R100 and R150 at various salt concentrations during the one-dimensional diffusion. From these trajectories, one can see that the protein can diffuse randomly along DNA at these different salt concentrations. At relatively low salt concentrations(50 mM,150 mM),where the electrostatic interactions are strong, the one-dimensional diffusion occurs dominantly by the rotation-coupled mechanism. While at high salt concentration(450 mM),the electrostatic interactions become weaker, which makes the protein move more freely along the DNA axis, leading to weak correlation between translational and rotational motions. To more clearly characterize the relationship between rotation and translation during one-dimensional diffusion, we also showed the rationaltranslational correlation plot in Fig.5(b)for all trajectories at three different salt concentrations. From these results,one can see an obvious correlation between the rotation and translation at low salt concentrations (50 mM), while it disappears at relatively high salt concentration (450 mM). Particularly,the rotational–translational correlation is still significant at the physiological relevant salt concentration(150 mM).

    3.2. Structural properties of XPA

    In addition to the diffusion features during this dynamical process, the structural properties of XPA were also studied,since its structure dynamics will surely affect the binding of XPA with other binding partners. The results for these two different DNA sequences are shown in Fig. 6. From our results, one can clearly see that both the DNA binding and varying salt concentration can change the structural properties of XPA and the difference in the effects of R100 and R150 on protein structure at each salt concentration is negligible.In detail, as shown in Fig. 6(a), at low salt concentration(50 mM), due to the strong intra-molecule electrostatic interactions in XPA, these two disordered regions pack against the DBD, leading to small radius of gyration (Rg) of XPA(less than 30 °A) in the absence of DNA, but it increases a lot (around 40 °A) upon binding DNA, which may be caused by the occupancy of intra-molecule interacting sites by DNA.While at high salt concentration (450 mM), the electrostatic interactions between intra and inter-molecule is decreased significantly and consequently the intra and inter molecule motions are more freely, leading to largeRgof XPA with or without DNA, and the effects of DNA binding is much deprecated as well (Fig. 6(b)). What is more, in Figs. 6(c)and 6(d), we also plotted the distribution of contact number between DBD and the disordered regions. In accordance with the above analysis, there are a large number of contacts between N(C)-terminal and DBD,and it decreases much unpon binding DNA at low salt concentration (Fig. 6(c)),while at high salt concentration,the contact number decreases significantly with or without DNA.Notably,the contact number in the absence of DNA is still less compared with binding DNA.From the analysis of the contact number,it becomes clear that both increasing salt concentrations and DNA binding can lead to the decrease of contact number,indicating more extended conformation of XPA. Based these observations, considering the dominate role of electrostatic interactions in describing the XPA–DNA interactioins and N(C)–DBD interactions in XPA, we can conclude that, both the DNA and salt concentration can act as modulating factors of the XPA structure.

    Fig.6. Effect of DNA binding and salt concentration on the conformation of XPA.(a)and(b)Rg distribution of XPA with(R100: blue,R150:purple)and without(black)DNA binding at the salt concentrations of 50 mM(a)and 450 mM(b). The representative structures of XPA are also shown,and the color scheme for domains is the same as in Fig.1(a). (c)and(d)Distribution of contact number formed between the DBD domain and the two disordered regions of XPA with(R100: blue,R150: purple)and without(black)DNA binding at the salt concentrations of 50 mM(c)and 450 mM(d).

    Nowdays, more and more studies have found that the conformation of protein is very sensitive to external stimulus, which play an important role in modulating the conformation of protein involed in various biological activities, like a cascade of reactions or signal transduction in cells,etc.[72–76]Interestingly, previously studies have also found that both salt concentration and binding DNA could affect the protein structure.[77–82]Hence,toghter with the findings that the XPA needs to cooperate with other binding partners during the NER pathway, and multiple binding sites for these partners are located at both DBD and disordered regions(N-terminal and Cterminal),[16,17,20,83–86]we can speculate that both the DNA itself and the salt concentration play a role in modulating the NER and at high salt concentration,the protein predominately adopt an extended conformation, ready for interacting with other partners,while at relative low salt concentration,the disordered regions tend to pack tightly with the DBD without DNA, however, upon DNA binding, the N(C)–DBD contacts are much reduced,which can lead to increased exposure of the binding sites in XPA for recruiting other partner proteins.

    4. Conclusion

    Here, we used a coarse-grained model to study the dynamical process of XPA diffusion along two different DNAs.Our results showed that the diffusion properties of XPA along these two DNAs are pretty similar,which indicates the robustness of dynamical behaviour of XPA along DNA.From more detailed perspective,our results showed that the XPA diffuses along the DNA with interface similar to the recognition motif found in the crystal structure, which may promote the search efficiency for the DNA target. Additionally, we also characterized the diffusion process of XPA along DNA. The XPA diffuses along DNA using combination of one-dimensional and three-dimensional search mechanisms,and increasing salt concentration can decrease the proportion of one-dimensional mechanism. During one-dimensional diffusion, the protein displays subdiffusive motions, which has also been observed in several biological system. Additionally, the correlation analysis of rotational and translational motions of DBD along DNA also showed that at relatively low salt concentration,the one-dimensional DBD translational motion is coupled to rotation,that is to say,sliding along helical pitch. Whereas at high salt concentration,the correlation disappears. All these results suggest that, like many other DNA binding proteins,[54,87,88]XPA maintains its position and orientation with respect to targeted DNA base pairs during one-dimensional diffusion at relatively low salt concentration, ensuring the efficient target search.

    In addition to binding to DNA, XPA has been found to act a key scaffold protein, interacting with many proteins involved in NER.[86]Therefore,we analyzed the structural features of XPA at various salt concentrations,in the presence and absence of DNA, respectively. The results showed that both DNA and salt concentration can modulate the protein structure. At low concentration,DNA binding can significantly increase theRgof protein and decrease the contact number between DBD and disordered regions,indicating more extended conformation of XPA.While at high concentration, the weak interaction between protein and DNA and intra-molecular interaction will let the protein to be in the extensible state as well. Based on these discussions, we can speculate that even at relatively low salt concentrations,when XPA bind to DNA,the contacts between DBD and other disordered regions are reduced, which tends to increase the exposure extent of the binding sites for other repair proteins,promoting the scaffold function of the XPA.

    Acknowledgment

    The authors thank the insightful discussions with Jun Wang,Wenfei Li,and Wei Wang.

    猜你喜歡
    張建
    鷸蚌相爭,漁人得利
    張建一:匠人匠心 征服世界的聲音
    綠色中國(2019年13期)2019-11-26 07:10:58
    本期焦點人物:張建一
    綠色中國(2019年13期)2019-11-26 07:10:50
    河鲀投毒案
    Les gar?ons sont-ils privilégiés?
    法語學習(2016年5期)2016-12-18 15:16:23
    那還用說
    那還用說
    那還用說
    愛你(2015年11期)2015-11-17 10:42:55
    那還用說
    肩關節(jié)生物力學
    亚洲欧美一区二区三区国产| 亚洲自偷自拍三级| 国产精品成人在线| 这个男人来自地球电影免费观看 | 日韩视频在线欧美| 在线看a的网站| 亚洲av男天堂| 99九九在线精品视频 | 3wmmmm亚洲av在线观看| 久久 成人 亚洲| 午夜激情福利司机影院| kizo精华| 91aial.com中文字幕在线观看| 精品一区二区免费观看| 色婷婷av一区二区三区视频| 成人黄色视频免费在线看| 菩萨蛮人人尽说江南好唐韦庄| 国产精品国产三级国产av玫瑰| 亚洲精品日韩在线中文字幕| 午夜av观看不卡| 最黄视频免费看| 亚洲天堂av无毛| 色婷婷av一区二区三区视频| 免费观看的影片在线观看| 亚洲精品一区蜜桃| 黑丝袜美女国产一区| 亚洲欧洲国产日韩| 最后的刺客免费高清国语| 亚洲人成网站在线播| 精品人妻熟女毛片av久久网站| 2022亚洲国产成人精品| 一级二级三级毛片免费看| 久久精品夜色国产| 亚洲性久久影院| 2022亚洲国产成人精品| 嫩草影院新地址| 日本黄色片子视频| 中文资源天堂在线| 中文欧美无线码| 国产精品一区二区在线不卡| 青春草国产在线视频| 免费看不卡的av| 天天操日日干夜夜撸| 国产欧美另类精品又又久久亚洲欧美| 麻豆成人午夜福利视频| 欧美变态另类bdsm刘玥| 自拍偷自拍亚洲精品老妇| av天堂久久9| 免费黄色在线免费观看| 亚洲精华国产精华液的使用体验| 国产真实伦视频高清在线观看| 青春草亚洲视频在线观看| 最近的中文字幕免费完整| 久久狼人影院| 在线看a的网站| 啦啦啦在线观看免费高清www| 高清欧美精品videossex| 国产亚洲av片在线观看秒播厂| av福利片在线观看| 欧美成人精品欧美一级黄| 中国美白少妇内射xxxbb| av在线观看视频网站免费| 免费久久久久久久精品成人欧美视频 | 国产黄频视频在线观看| 国产精品人妻久久久影院| 91久久精品国产一区二区三区| 欧美另类一区| 国产免费视频播放在线视频| 男人和女人高潮做爰伦理| 亚洲成人一二三区av| 少妇精品久久久久久久| 亚洲av欧美aⅴ国产| 国产成人精品婷婷| 性色avwww在线观看| 性色av一级| tube8黄色片| 国产午夜精品一二区理论片| 国产精品一区二区在线观看99| 中文天堂在线官网| 国产精品秋霞免费鲁丝片| 亚洲av.av天堂| 男女边吃奶边做爰视频| 校园人妻丝袜中文字幕| 久久午夜福利片| 日韩一区二区三区影片| 国产色爽女视频免费观看| 精品国产露脸久久av麻豆| 我的女老师完整版在线观看| 亚洲av免费高清在线观看| 欧美一级a爱片免费观看看| 国产免费福利视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 校园人妻丝袜中文字幕| 精品亚洲成a人片在线观看| 全区人妻精品视频| 亚洲精品国产成人久久av| 99热全是精品| 性高湖久久久久久久久免费观看| 26uuu在线亚洲综合色| 亚洲国产精品一区三区| 一本一本综合久久| 自线自在国产av| 午夜免费观看性视频| 免费看光身美女| 国产女主播在线喷水免费视频网站| 日韩伦理黄色片| 国产亚洲午夜精品一区二区久久| 日日爽夜夜爽网站| 在线观看人妻少妇| 这个男人来自地球电影免费观看 | 免费看av在线观看网站| 99精国产麻豆久久婷婷| 日韩中文字幕视频在线看片| 青春草国产在线视频| 免费看光身美女| 亚洲国产av新网站| 婷婷色综合www| 久久国产精品男人的天堂亚洲 | 欧美最新免费一区二区三区| 一本色道久久久久久精品综合| 亚洲国产精品成人久久小说| 人人妻人人爽人人添夜夜欢视频 | 国产中年淑女户外野战色| 色5月婷婷丁香| av免费观看日本| 精品少妇内射三级| 国产高清三级在线| 欧美日韩综合久久久久久| 亚洲精品日韩在线中文字幕| 久久久久久久久大av| 男女国产视频网站| 成人无遮挡网站| 啦啦啦啦在线视频资源| 国产欧美日韩精品一区二区| 国产有黄有色有爽视频| 丰满乱子伦码专区| 亚洲精品国产av蜜桃| 老司机影院毛片| 亚洲国产欧美在线一区| 一级黄片播放器| 内射极品少妇av片p| 成人国产麻豆网| 中文字幕精品免费在线观看视频 | 久久久久久久久久久丰满| av.在线天堂| √禁漫天堂资源中文www| 人妻系列 视频| 亚洲国产日韩一区二区| 亚洲自偷自拍三级| 嘟嘟电影网在线观看| 欧美最新免费一区二区三区| 国产视频首页在线观看| 三上悠亚av全集在线观看 | 99热全是精品| 你懂的网址亚洲精品在线观看| 免费久久久久久久精品成人欧美视频 | 国产一区二区三区av在线| 国产成人免费无遮挡视频| 亚洲av成人精品一区久久| 一级,二级,三级黄色视频| 久久久久久久精品精品| 国产淫语在线视频| 少妇人妻一区二区三区视频| 午夜老司机福利剧场| 中文字幕人妻熟人妻熟丝袜美| 熟女电影av网| 日韩中字成人| 一级毛片久久久久久久久女| av女优亚洲男人天堂| 中文字幕免费在线视频6| 久久av网站| 一级,二级,三级黄色视频| 国产亚洲最大av| 特大巨黑吊av在线直播| 九九在线视频观看精品| 丰满迷人的少妇在线观看| 狠狠精品人妻久久久久久综合| 在线天堂最新版资源| 女人精品久久久久毛片| 2021少妇久久久久久久久久久| 80岁老熟妇乱子伦牲交| 日本vs欧美在线观看视频 | 亚洲真实伦在线观看| 亚洲国产精品一区三区| 亚洲欧美一区二区三区国产| 丝袜喷水一区| 香蕉精品网在线| 少妇人妻精品综合一区二区| 亚洲av成人精品一区久久| 色5月婷婷丁香| 在线播放无遮挡| 大话2 男鬼变身卡| 成年美女黄网站色视频大全免费 | 免费av中文字幕在线| 免费看日本二区| 九九在线视频观看精品| 久久99热这里只频精品6学生| av在线app专区| 美女大奶头黄色视频| 亚洲精品久久午夜乱码| 99久久精品一区二区三区| 哪个播放器可以免费观看大片| 亚洲va在线va天堂va国产| 日韩欧美一区视频在线观看 | 久久久久久久亚洲中文字幕| 日韩中字成人| 国产日韩欧美视频二区| 亚洲经典国产精华液单| 一级毛片电影观看| 丁香六月天网| 少妇丰满av| 五月天丁香电影| 久久久久久久久久成人| 午夜久久久在线观看| 国产精品蜜桃在线观看| 亚洲精品乱码久久久久久按摩| 国产成人免费无遮挡视频| 国产淫语在线视频| 亚洲av.av天堂| 青春草国产在线视频| 九九在线视频观看精品| 最近中文字幕高清免费大全6| 国产一区有黄有色的免费视频| av在线观看视频网站免费| 日韩欧美精品免费久久| 亚洲国产日韩一区二区| av一本久久久久| 老熟女久久久| 免费观看a级毛片全部| 老司机影院成人| 国产成人freesex在线| 亚州av有码| 成人二区视频| 男女啪啪激烈高潮av片| 在现免费观看毛片| 亚洲经典国产精华液单| 青春草亚洲视频在线观看| 精品亚洲成国产av| 欧美日韩精品成人综合77777| 国产精品熟女久久久久浪| 看免费成人av毛片| 美女xxoo啪啪120秒动态图| 欧美精品高潮呻吟av久久| 一区在线观看完整版| 亚洲精品日本国产第一区| 国产在线免费精品| 午夜久久久在线观看| 久久精品国产亚洲av天美| 汤姆久久久久久久影院中文字幕| 亚洲综合精品二区| 女性生殖器流出的白浆| 日韩成人av中文字幕在线观看| 中文字幕久久专区| 免费播放大片免费观看视频在线观看| 国产69精品久久久久777片| 亚洲国产精品一区二区三区在线| 国产免费视频播放在线视频| 内射极品少妇av片p| 少妇的逼好多水| 久久精品国产a三级三级三级| 国产亚洲5aaaaa淫片| 一级毛片我不卡| 亚洲性久久影院| 久久久国产一区二区| 精品99又大又爽又粗少妇毛片| 午夜日本视频在线| 2021少妇久久久久久久久久久| 黄色配什么色好看| 在线观看美女被高潮喷水网站| 国产色婷婷99| 久久久久视频综合| 高清在线视频一区二区三区| 久久久久精品久久久久真实原创| 成人二区视频| 自拍偷自拍亚洲精品老妇| 五月天丁香电影| 亚洲成色77777| 美女xxoo啪啪120秒动态图| 亚洲三级黄色毛片| 亚洲综合色惰| 久久精品国产亚洲av天美| 街头女战士在线观看网站| 女性被躁到高潮视频| .国产精品久久| 9色porny在线观看| h视频一区二区三区| 极品人妻少妇av视频| 在线观看www视频免费| 成人毛片a级毛片在线播放| 好男人视频免费观看在线| 欧美激情国产日韩精品一区| 日韩强制内射视频| 日韩欧美一区视频在线观看 | kizo精华| 男女国产视频网站| 人妻 亚洲 视频| 女人精品久久久久毛片| 亚洲精品一区蜜桃| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频 | 中文字幕精品免费在线观看视频 | 久久亚洲国产成人精品v| 大片电影免费在线观看免费| 国产精品麻豆人妻色哟哟久久| 日韩伦理黄色片| 啦啦啦啦在线视频资源| 国产精品福利在线免费观看| 建设人人有责人人尽责人人享有的| 国产成人精品无人区| 亚洲va在线va天堂va国产| 少妇高潮的动态图| 日日爽夜夜爽网站| 黄色一级大片看看| 亚洲色图综合在线观看| 亚洲性久久影院| 伦理电影免费视频| 久久毛片免费看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 男女边吃奶边做爰视频| 日韩一区二区视频免费看| 18禁裸乳无遮挡动漫免费视频| 久久久久久久国产电影| 欧美日韩精品成人综合77777| 日韩一区二区三区影片| 婷婷色麻豆天堂久久| 亚洲精品国产av成人精品| 亚洲经典国产精华液单| av卡一久久| 日本黄色日本黄色录像| 啦啦啦中文免费视频观看日本| 尾随美女入室| 人人澡人人妻人| 日韩欧美一区视频在线观看 | 色5月婷婷丁香| 亚洲国产av新网站| 丰满饥渴人妻一区二区三| 一级,二级,三级黄色视频| 丝袜脚勾引网站| 大话2 男鬼变身卡| 国产无遮挡羞羞视频在线观看| 久久热精品热| 曰老女人黄片| 永久免费av网站大全| 高清欧美精品videossex| 国产av精品麻豆| 99久久综合免费| 女人久久www免费人成看片| 国产黄片视频在线免费观看| 免费看光身美女| 观看美女的网站| 尾随美女入室| 亚洲精品久久午夜乱码| 欧美激情极品国产一区二区三区 | 欧美国产精品一级二级三级 | 国产精品一区二区在线不卡| 日韩成人伦理影院| 在线观看美女被高潮喷水网站| 国产伦理片在线播放av一区| 欧美日韩视频高清一区二区三区二| 尾随美女入室| 久久青草综合色| 高清视频免费观看一区二区| 久久青草综合色| 免费看光身美女| 亚洲精品久久午夜乱码| 少妇高潮的动态图| 久久精品久久精品一区二区三区| www.av在线官网国产| 国产高清有码在线观看视频| 少妇高潮的动态图| 精品卡一卡二卡四卡免费| 亚洲欧美日韩东京热| 久久精品久久精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 日日爽夜夜爽网站| av天堂中文字幕网| 99久久精品国产国产毛片| av.在线天堂| 2021少妇久久久久久久久久久| 日本-黄色视频高清免费观看| 久久国内精品自在自线图片| 日韩欧美 国产精品| 一本一本综合久久| 日韩一本色道免费dvd| 精品人妻熟女av久视频| 日产精品乱码卡一卡2卡三| 麻豆乱淫一区二区| 国产亚洲5aaaaa淫片| 天堂俺去俺来也www色官网| 在线观看免费日韩欧美大片 | 在线精品无人区一区二区三| 免费观看无遮挡的男女| 成人亚洲欧美一区二区av| 亚洲国产最新在线播放| 制服丝袜香蕉在线| 国产黄色视频一区二区在线观看| 免费人妻精品一区二区三区视频| 国产精品久久久久久久久免| 日韩不卡一区二区三区视频在线| 精品亚洲成国产av| 免费久久久久久久精品成人欧美视频 | 在线观看国产h片| 国产伦理片在线播放av一区| av国产精品久久久久影院| 久久久久国产网址| 精品国产国语对白av| 免费av不卡在线播放| 成人无遮挡网站| 免费观看av网站的网址| 精品一区二区三区视频在线| 国产精品国产三级专区第一集| 丰满饥渴人妻一区二区三| 亚洲精品日韩av片在线观看| av线在线观看网站| 国产精品人妻久久久影院| 精品视频人人做人人爽| av天堂久久9| 在线精品无人区一区二区三| 精品少妇内射三级| 午夜福利在线观看免费完整高清在| 在现免费观看毛片| 汤姆久久久久久久影院中文字幕| 极品教师在线视频| 亚洲国产精品专区欧美| 日韩电影二区| 国国产精品蜜臀av免费| 内射极品少妇av片p| 少妇猛男粗大的猛烈进出视频| 大陆偷拍与自拍| 亚洲va在线va天堂va国产| 亚洲国产精品一区二区三区在线| 久久人人爽人人片av| 自拍偷自拍亚洲精品老妇| 如何舔出高潮| 女人精品久久久久毛片| 午夜免费观看性视频| 成人18禁高潮啪啪吃奶动态图 | 天堂中文最新版在线下载| 人妻少妇偷人精品九色| 亚洲真实伦在线观看| 特大巨黑吊av在线直播| 桃花免费在线播放| 精品一品国产午夜福利视频| 美女xxoo啪啪120秒动态图| www.av在线官网国产| 能在线免费看毛片的网站| 亚洲国产精品一区二区三区在线| 国产乱来视频区| 亚洲精品日韩av片在线观看| av在线app专区| 午夜激情久久久久久久| 日韩一区二区视频免费看| 久久人人爽av亚洲精品天堂| 免费av不卡在线播放| 男人和女人高潮做爰伦理| 国产精品一区二区在线观看99| 大片免费播放器 马上看| 国产黄片美女视频| 成年人午夜在线观看视频| 日日啪夜夜撸| 18禁动态无遮挡网站| 日韩制服骚丝袜av| 国产欧美日韩精品一区二区| 天美传媒精品一区二区| 人人妻人人澡人人爽人人夜夜| 美女cb高潮喷水在线观看| 日韩一本色道免费dvd| 男人狂女人下面高潮的视频| 我的女老师完整版在线观看| 赤兔流量卡办理| 啦啦啦中文免费视频观看日本| 人体艺术视频欧美日本| 51国产日韩欧美| 成人黄色视频免费在线看| 97在线人人人人妻| 欧美少妇被猛烈插入视频| 亚洲,欧美,日韩| 国模一区二区三区四区视频| 久久免费观看电影| 三上悠亚av全集在线观看 | a级片在线免费高清观看视频| 亚洲精品久久午夜乱码| 肉色欧美久久久久久久蜜桃| 亚洲精品成人av观看孕妇| 国产欧美日韩综合在线一区二区 | 国内揄拍国产精品人妻在线| 亚洲第一av免费看| 18禁在线播放成人免费| 亚洲成人av在线免费| 精品一区二区三区视频在线| 久久国产乱子免费精品| 又黄又爽又刺激的免费视频.| 黑人猛操日本美女一级片| 日韩,欧美,国产一区二区三区| 欧美日韩国产mv在线观看视频| 22中文网久久字幕| 三上悠亚av全集在线观看 | 亚洲,欧美,日韩| 亚洲国产最新在线播放| 亚洲av国产av综合av卡| 丝瓜视频免费看黄片| 日韩av免费高清视频| 人人妻人人爽人人添夜夜欢视频 | 人妻制服诱惑在线中文字幕| 久久国内精品自在自线图片| 街头女战士在线观看网站| 高清在线视频一区二区三区| 国产亚洲最大av| 色94色欧美一区二区| 丝袜脚勾引网站| 婷婷色综合大香蕉| 中文字幕久久专区| av国产精品久久久久影院| 制服丝袜香蕉在线| 亚洲av综合色区一区| 色吧在线观看| 亚洲国产精品999| 日韩欧美精品免费久久| 国产精品国产三级专区第一集| 99久久人妻综合| 一级毛片电影观看| 国产精品一区二区在线观看99| 水蜜桃什么品种好| 久久久久精品久久久久真实原创| 精品一品国产午夜福利视频| 一二三四中文在线观看免费高清| 成人毛片a级毛片在线播放| 色婷婷久久久亚洲欧美| 黑丝袜美女国产一区| 一区二区三区免费毛片| 欧美区成人在线视频| 春色校园在线视频观看| 婷婷色麻豆天堂久久| 久久精品国产a三级三级三级| 国国产精品蜜臀av免费| 色94色欧美一区二区| 狂野欧美激情性bbbbbb| 日韩制服骚丝袜av| 亚洲精华国产精华液的使用体验| 丰满少妇做爰视频| 国产精品偷伦视频观看了| 久久久国产欧美日韩av| 成人18禁高潮啪啪吃奶动态图 | 天天躁夜夜躁狠狠久久av| 国产爽快片一区二区三区| 色哟哟·www| 久久久国产精品麻豆| 国产 精品1| 精品国产一区二区久久| 蜜桃在线观看..| 插阴视频在线观看视频| 欧美日韩在线观看h| av在线老鸭窝| 少妇裸体淫交视频免费看高清| 日本免费在线观看一区| 免费看av在线观看网站| 啦啦啦在线观看免费高清www| 国产欧美亚洲国产| 国产美女午夜福利| 国产成人91sexporn| 欧美老熟妇乱子伦牲交| 免费观看的影片在线观看| 成人黄色视频免费在线看| 亚洲三级黄色毛片| 午夜福利视频精品| 91午夜精品亚洲一区二区三区| 伊人久久国产一区二区| 丰满迷人的少妇在线观看| 免费观看在线日韩| 欧美日韩亚洲高清精品| 免费久久久久久久精品成人欧美视频 | 国语对白做爰xxxⅹ性视频网站| 欧美亚洲 丝袜 人妻 在线| 国产亚洲av片在线观看秒播厂| 国产日韩一区二区三区精品不卡 | 久久久午夜欧美精品| 91成人精品电影| 精品少妇久久久久久888优播| 边亲边吃奶的免费视频| 99久久人妻综合| 美女xxoo啪啪120秒动态图| 国产免费福利视频在线观看| 国产91av在线免费观看| 狂野欧美激情性xxxx在线观看| 日本免费在线观看一区| 欧美日韩亚洲高清精品| a级毛片免费高清观看在线播放| av福利片在线观看| 另类精品久久| 一级片'在线观看视频| 免费看光身美女| 亚洲人与动物交配视频| 成人国产av品久久久| 久久99一区二区三区| 亚洲国产欧美日韩在线播放 | 成年人免费黄色播放视频 | 制服丝袜香蕉在线| 亚洲人成网站在线播| 国产高清国产精品国产三级| 亚洲电影在线观看av| 亚洲精品中文字幕在线视频 | 亚洲国产毛片av蜜桃av| 中文字幕av电影在线播放| 亚洲激情五月婷婷啪啪| 国产色爽女视频免费观看| 日本爱情动作片www.在线观看| 日韩人妻高清精品专区| 乱码一卡2卡4卡精品| 久久青草综合色| 男人舔奶头视频| 一级爰片在线观看| 美女福利国产在线| 欧美亚洲 丝袜 人妻 在线| 男人添女人高潮全过程视频| 精品国产国语对白av|