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

    氨基酸團簇及磷脂在蒙脫石表面的電子傳遞計算模擬研究*

    2016-12-03 03:30:28董發(fā)勤孫仕勇郭玉婷霍婷婷
    功能材料 2016年4期
    關(guān)鍵詞:電子密度蒙脫石磷脂

    周 青,董發(fā)勤,邊 亮,孫仕勇,郭玉婷,4,霍婷婷

    (1. 西南科技大學 環(huán)境與資源學院,四川 綿陽 621010;2. 固體廢物處理與資源化教育部重點實驗室,四川 綿陽 621010;3. 中國科學院 新疆理化技術(shù)研究所,烏魯木齊 830011;4. 西南科技大學 生命科學與工程學院,四川 綿陽 621010)

    ?

    氨基酸團簇及磷脂在蒙脫石表面的電子傳遞計算模擬研究*

    周 青1,2,董發(fā)勤2,邊 亮3,孫仕勇1,2,郭玉婷2,4,霍婷婷2

    (1. 西南科技大學 環(huán)境與資源學院,四川 綿陽 621010;2. 固體廢物處理與資源化教育部重點實驗室,四川 綿陽 621010;3. 中國科學院 新疆理化技術(shù)研究所,烏魯木齊 830011;4. 西南科技大學 生命科學與工程學院,四川 綿陽 621010)

    分子模擬以其在分子尺度上的可視性已成為界面作用研究的一項重要手段。本文通過分子動力學研究氨基酸團簇及磷脂分子與可吸入顆粒物蒙脫石的吸附來模擬蒙脫石進入人體后與人體細胞膜膜蛋白及磷脂分子的相互作用過程。另外在計算體系中加入重金屬離子Cr(Ⅲ)、Cr(Ⅵ)、Cd(Ⅱ)、Pb(Ⅱ)以了解其對該過程的影響。模擬結(jié)果證明,氨基酸和磷脂均與蒙脫石具有較好的親和性,但氨基酸強于磷脂。重金屬離子則起到了促進作用。通過密度泛函理論研究作用前后分子的電子結(jié)構(gòu)發(fā)現(xiàn),相互作用過程中存在以蒙脫石為供體、氨基酸和磷脂為受體的電子傳遞。

    蒙脫石;氨基酸;磷脂;分子模擬;重金屬離子;電子傳遞

    0 引 言

    大氣顆粒物污染日趨嚴重,大氣顆粒物被吸入人體后會刺激巨噬細胞、呼吸道上皮細胞等產(chǎn)生自由基[1]。具有強氧化性的自由基由于含有未偶電子,不僅可以使細胞膜的不飽和脂肪酸發(fā)生脂質(zhì)過氧化,而且可以攻擊膜蛋白,核酸等生物分子使之發(fā)生鏈內(nèi)或鏈間交聯(lián)、斷裂,從而引起細胞的氧化損傷,甚至凋亡[2-6]。研究發(fā)現(xiàn),顆粒的成分及比表面積與其誘導產(chǎn)生細胞氧化損傷的能力有直接關(guān)系[3-4]。

    作為大氣顆粒物礦物成分的一種[4, 7],蒙脫石的片層具有極高的橫縱比(約500)和比表面積(約800 m2/g)[8]。盡管關(guān)于蒙脫石粉塵的致病案例早有報道,其致病機理的研究卻并不多見[9-11]。然而,由于在自然條件下蒙脫石層帶負電且有可交換的層間陽離子,其與生物分子如氨基酸、蛋白質(zhì)、脂多糖等的相互作用研究相對較多,尤其是蒙脫石吸附氨基酸、蛋白質(zhì)[12-15]。一定程度上,蒙脫石和氨基酸以及蛋白質(zhì)的相互作用是其二者親和力的一種反映,研究二者的吸附作用對于研究蒙脫石致病機理有一定的幫助。氨基酸、蛋白質(zhì)、磷脂等生物分子在蒙脫石表面的吸附普遍認為由離子交換、靜電作用、氫鍵、范德華力等完成的[14]。

    近年來,隨著計算模擬技術(shù)的發(fā)展,分子動力學模擬已廣泛應用于生物大分子結(jié)構(gòu)-功能關(guān)系的研究中,而應用于礦物與生物大分子的交互作用的較少[13, 16]。本文建立蒙脫石模型以及由多肽(谷氨酸、甘氨酸、絲氨酸、精氨酸)以及磷脂單分子層組成的細胞膜模型。通過計算機模擬二者在Cr(Ⅲ)、Cr(Ⅵ)、Cd(Ⅱ)、Pb(Ⅱ) 4種重金屬離子影響下的相互作用[17],以獲得二者的形變以及二者的電子傳輸機制等相關(guān)信息。通過分子模擬,從分子視角在一定程度上洞悉蒙脫石粉塵被吸入肺泡之后與人體細胞細胞膜的相互作用,以期為粉塵的致病機理提供一定的參考。

    1 計算模擬方法

    模擬計算軟件采用的是Materials Studio(Accelrys, USA)。采用分子力學建立谷氨酸(Glu)、甘氨酸(Gly)、絲氨酸(Ser)、精氨酸(Arg) 4種氨基酸團簇和磷脂(Phospholipid)等生物分子(BM)團簇模型[18](詳見表1)及蒙脫石(MMT)模型[19]。

    表1 生物分子模型的基本參數(shù)

    蒙脫石對4種氨基酸團簇的吸附位點及吸附量通過巨正則蒙脫卡洛計算獲得。計算步數(shù)為10 000 000步,計算方法采用的是Metropolis法,采用Universal力場,運動方程由Verlet算法整合[20]。

    隨后采用Universal力場進行分子動力學弛豫。先NPT系綜后NVT系綜,弛豫時間均為1 000 ps,步長設為1 fs,溫度設為310 K,壓力為標準大氣壓0.1 MPa。長程靜電作用計算采用Atom-based方法,短程范德華作用計算采用Ewald & group方法。弛豫過程中溫度采用Nose-Hoover熱浴調(diào)節(jié),壓力由Berendsen恒壓器調(diào)節(jié)[21]。

    圖1 模擬計算過程圖示

    最后利用密度泛函理論(DFT)進行量化計算。采用Perdew-Burke-Ernzerhof(PBE)廣義梯度近似(GGA)交換泛函計算生物分子與礦物結(jié)構(gòu)之間的電子傳遞[22]。平面波截斷能量收斂為500 eV,采用3×3×3的k點抽樣在布里淵區(qū)內(nèi)進行倒易空間積分。整個計算過程如圖1所示。

    2 結(jié)果討論

    2.1 相互作用強弱及相互作用位點

    蒙脫石與生物分子的親和性實際上表示為生物分子在蒙脫石表面的吸附。吸附作用的強弱在數(shù)值上表現(xiàn)為吸附能的不同,吸附能的計算公式為Eadsorption=E(adsobate+substrate)-Eadsorption-Esubstrate。吸附能為負值表示吸附行為是自發(fā)進行的,吸附后體系整體能量降低且體系穩(wěn)定,吸附過程是放熱的。吸附能的絕對值越大,表明吸附作用越強烈,蒙脫石與生物分子的相互作用也就越明顯。

    本文采用蒙脫卡洛計算得到了蒙脫石吸附生物分子的吸附能及吸附后生物分子的變形能,計算結(jié)果如表2所示。結(jié)果顯示,蒙脫石與5種生物分子之間的吸附能皆為負值,說明蒙脫石與5種生物分子均存在較強相互作用。其中,蒙脫石對4種氨基酸團簇的作用強于磷脂分子,4種氨基酸團簇則大體一致。較強的相互作用則導致生物分子較大的形變,數(shù)值上表達為變形能的絕對值較大,稍后將結(jié)合鍵能及扭轉(zhuǎn)能的變化綜合討論。

    表2 蒙脫石吸附生物分子的吸附能與變形能(kJ/mol)

    Table 2 Adsorption and deformation energy of biomolecules on montmorillonite surface (kJ/mol)

    能量PhospholipidGlySerGluArg吸附能-2052.56-2549.35-2549.97-2542.67-2552.49變形能-2185.72-2673.44-2670.71-2661.77-2668.83

    2.2 相互作用前后生物分子的能量變化

    生物分子的形變更加具體地反映為鍵能、扭轉(zhuǎn)能等的變化。本文不僅計算了包括鍵能、扭轉(zhuǎn)能在內(nèi)的價鍵能,還計算了包括靜電作用能以及范德華能在內(nèi)的非鍵合能,具體與圖3所示。由圖3可以看出,相似結(jié)構(gòu)的生物分子的價鍵能與分子的原子數(shù)呈一定的正相關(guān)。與蒙脫石作用后,5種生物分子的鍵能及扭轉(zhuǎn)能均升高。由于蒙脫石分子對生物分子界面原子的作用,使得生物分子不得不縮短鍵長以加強對表層原子的束縛,生物分子的鍵能相應變大。同理,由于蒙脫石的作用,緊縛的生物分子柔性減弱,鍵角旋轉(zhuǎn)、扭轉(zhuǎn)以及面外彎曲能力較弱,扭轉(zhuǎn)能升高。重金屬離子的存在則加強了蒙脫石與生物分子的作用,鍵能和扭轉(zhuǎn)能進一步升高。值得一提的是,在重金屬離子的協(xié)同作用下,個別生物分子的鍵能竟達873.80 kJ/mol。

    與價鍵能相比,非鍵合能的變化顯得似乎沒有規(guī)律可循。不同于價鍵能,非鍵合能的變化是生物分子內(nèi)、生物分子與生物分子之間以及生物分子與蒙脫石之間作用疊加的結(jié)果??偟膩砜?,長程靜電作用強于短程范德華作用,即在相互作用過程中,靜電作用力起主導作用。磷脂與Gly體系以范德華引力為主,其它生物分子體系以范德華斥力為主;磷脂與谷氨酸體系以靜電斥力為主,其它生物分子體系以引力為主。

    圖2 分子動力學弛豫之后蒙脫石與生物分子相對位置圖

    圖3 生物分子相互作用前后體系的能量變化

    2.3 蒙脫石與生物分子之間的電子傳遞

    確定了相互作用基團后,采用DFT理論進行量子化學計算,得到了相互作用前后體系分子的電子結(jié)構(gòu),從而得到二者之間的電子傳遞情況。

    分子的表面電勢是分子軌道及電子密度的反映,通過繪制分子的表面電勢等高線圖,可以定性地了解電子的傳遞方向,電勢低的區(qū)域的有相對自由的電子,而電勢高的區(qū)域則表示該區(qū)域?qū)﹄娮哟嬖谳^強的吸引力,電勢越高,得電子能力也就越強。從圖4可以看出,蒙脫石整體勢能較低,其中側(cè)面邊緣的硅氧四面體的Si原子區(qū)域存在明顯的負值區(qū),這主要是由于邊緣的Si的外層電子p軌道上的兩個電子未參與成鍵。未成鍵的孤對電子由于沒有鍵的束縛,能量又處于最高、次高,易于偏向電子受體,從而從蒙脫石分子中失去。而生物分子的外層密布著高電勢區(qū),這主要是由于外層H原子的電子由于受到內(nèi)層原子吸引,幾乎向外裸露的H原子核有較強的吸引電子的能力??傊瑥内厔萆?,蒙脫石易于失去電子,而生物分子易于得到電子。

    為進一步定量地了解二者之間的電子傳遞,根據(jù)量化計算結(jié)果,對蒙脫石與生物分子的作用基團進行Mulliken布局分析。Mulliken電荷是原子周圍電子密度高低的直接反映,電荷越高,電荷密度越低。Mulliken電荷的分子結(jié)果如表3所示,蒙脫石的Al、Si原子在作用后,電子密度變低,其中Al的變化尤為明顯;氧原子的電子密度增高。氨基酸—NH2和—COOH的電子密度增高。磷脂分子酯基—COO—電子密度增高,而親水端上—H2PO4的電子密度反而變低。

    圖4 蒙脫石與生物分子的表面電勢等值線圖

    作用基團作用前作用后PhospholipidGlySerGluArg蒙脫石Al1.5602.2292.2732.4032.2132.252Si2.0152.0682.1802.3362.3512.460O(四面體)-0.995-1.194-1.187-1.310-1.221-1.303O(八面體)-0.923-1.204-1.209-1.373-1.275-1.269氨基酸—NH2-0.102~-0.109—-0.173-0.158-0.139-0.119—COOH-0.481~-0.372—-0.527-0.659-0.572-0.648磷脂—COO—-0.279-0.366—————H2PO40.3460.418————

    與表面電勢顯示的趨勢較為一致,蒙脫石主要為電子供體,而生物分子為電勢受體,磷脂的親水基團—H2PO4除外。只是不同的是蒙脫石最大的電子供體為Al而非Si,這主要是因為孤對電子只在蒙脫石片層的側(cè)面Si原子中存在,而且Al的電子吸引能力比Si弱。Al、Si原子的電子經(jīng)由表層O原子傳遞給生物分子的作用基團,故蒙脫石表層氧原子的電子密度也有所增高。通過對磷脂作用基團—H2PO4的單個原子進行布局分析發(fā)現(xiàn)P原子為主要的電子供體(電荷由2.110 e升高為2.316 e)。由于重金屬離子是以蒙脫卡洛方法加入計算體系,以場的形式存在,其不存在電子得失,但重金屬離子加強了蒙脫石及生物分子的極化,并在傳遞過程中起到介導作用。

    3 結(jié) 論

    (1) 通過分子模擬計算發(fā)現(xiàn)蒙脫石在靜電力為主的作用力下,與蛋白質(zhì)和磷脂分子通過—NH2、—COOH、—H2PO4、—COO—相互作用,且蛋白質(zhì)的作用強于磷脂。在蒙脫石的作用下,生物分子發(fā)生較大的形變,鍵能增強,柔性變低。重金屬離子的存在則加強了該過程,發(fā)揮了協(xié)同作用。

    (2) 通過對相互作用后體系的電子結(jié)構(gòu)計算發(fā)現(xiàn)蒙脫石與生物分子的活性基團之間發(fā)生了以蒙脫石的Al、Si為供體、生物分子為受體的電子傳遞。

    [1] Nel A, Xia T, Madler L, et al. Toxic potential of materials at the nanolevel [J]. Science, 2006, 311(5761): 622-627.

    [2] Ciencewicki J, Trivedi S, Kleeberger S R. Oxidants and the pathogenesis of lung diseases [J]. J Allergy Clin Immun, 2008, 122(3): 456-468.

    [3] Dick C A J, Singh P, Daniels M, et al. Murine pulmonary inflammatory responses following instillation of size-fractionated ambient particulate matter [J]. Journal of Toxicology and Environmental Health-Part A, 2003, 66(23): 2193-2207.

    [4] Ghio A J, Carraway M S, Madden M C. Composition of air pollution particles and oxidative stress in cells, tissues, and living systems [J]. J Toxicol Env Heal B, 2012, 15(1): 1-21.

    [5] Gonzalez-Flecha B. Oxidant mechanisms in response to ambient air particles [J]. Molecular aspects of medicine, 2004, 25(1-2): 169-182.

    [6] Imrich A, Ning Y Y, Lawrence J, et al. Alveolar macrophage cytokine response to air pollution particles: oxidant mechanisms [J]. Toxicology and applied pharmacology, 2007, 218(3): 256-264.

    [7] Menendez I, Perez-Chacon E, Mangas J, et al. Dust deposits on La graciosa island (canary islands, spain): texture, mineralogy and a case study of recent dust plume transport [J]. Catena, 2014, 117: 133-144.

    [8] Gu C, Liu C, Johnston C T, et al. Pentachlorophenol radical cations generated on Fe(Ⅲ)-montmorillonite initiate octachlorodibenzo-p-dioxin formation in clays: density functional theory and fourier transform infrared studies [J]. Environmental Science & Technology, 2011, 45(4): 1399-1406.

    [9] Gibbs A R, Pooley F D. Fuller's earth (montmorillonite) pneumoconiosis [J]. Occupational and Environmental Medicine, 1994, 51(9): 644-646.

    [10] Elmore A R, Andersen F A, Pan C I R E. Final report on the safety assessment of aluminum silicate, calcium silicate, magnesium aluminum silicate, magnesium silicate, magnesium trisilicate, sodium magnesium silicate, zirconium silicate, attapulgite, bentonite, fuller's earth, hectorite, kaolin, lithium magnesium silicate, lithium magnesium sodium silicate, montmorillonite, pyrophyllite, and zeolite [J]. Int J Toxicol, 2003, 22: 37-102.

    [11] Fujimoto K, Muller N L, Kato S, et al. Pneumoconiosis in rush mat workers exposed to clay dye "sendo" dust: clinical, radiologic, and histopathologic features in seven patients [J]. Chest, 2004, 125(2): 737-743.

    [12] Kolman K, Makowski M M, GOLRIZ A A, et al. Adsorption, aggregation, and desorption of proteins on smectite particles [J]. Langmuir: the ACS Journal of Surfaces and Colloids, 2014, 30(39): 11650-11659.

    [13] Assifaoui A, Huault L, Maissiat C, et al. Structural studies of adsorbed protein (betalactoglobulin) on natural clay (montmorillonite) [J]. Rsc Adv, 2014, 4(105): 61096-61103.

    [14] Yu W H, Li N, Tong D S, et al. Adsorption of proteins and nucleic acids on clay minerals and their interactions: a review [J]. Applied Clay Science, 2013, 80-81: 443-452.

    [15] Tran A T T, James B J. A study the interaction forces between the bovine serum albumin protein and montmorillonite surface [J]. Colloid Surface A, 2012, 414: 104-114.

    [16] Liu X D, Lu X C, Meijer E J, et al. Atomic-scale structures of interfaces between phyllosilicate edges and water [J]. Geochim Cosmochim Ac, 2012, 81: 56-68.

    [17] Wang M, Deng Y Q, Dong F Q, et al. Speciation analysis of cadmium in dust fall about northern China towns [J]. Asian J Chem, 2013, 25(13): 7522-7526.

    [18] Bian L, Xu J B, Song M X, et al. Designing perovskite BFO(111) membrane as an electrochemical sensor for detection of amino acids: a simulation study [J]. Journal of Molecular Structure, 2015, 1099: 1-9.

    [19] Viani A, Gualtieri A F, Artioli G. The nature of disorder in montmorillonite by simulation of X-ray powder patterns [J]. American Mineralogist, 2002, 87(7): 966-975.

    [20] Joshi P, Shewale V, Pandey R, et al. Site specific interaction between ZnO nanoparticles and tryptophan: a first principles quantum mechanical study [J]. Physical Chemistry Chemical Physics, 2011, 13(2): 476-479.

    [21] Costa D, Tougerti A, Tielens F, et al. DFT study of the adsorption of microsolvated glycine on a hydrophilic amorphous silica surface [J]. Physical Chemistry Chemical Physics, 2008, 10(42): 6360-6368.

    [22] Almora-Barrios N, Austen K F, De Leeuw N H. Density functional theory study of the binding of glycine, proline, and hydroxyproline to the hydroxyapatite (0001) and (0110) surfaces [J]. Langmuir: the ACS Journal of Surfaces and Colloids, 2009, 25(9): 5018-5025.

    [23] Ramos M E, Huertas F J. Adsorption of glycine on montmorillonite in aqueous solutions [J]. Applied Clay Science, 2013, 80-81: 10-7.Electron transition of amino acids clusters and phospholipid on montmorillonite surface: a simulation research

    2. Key Laboratory of Solid Waste Treatment and Resource Recycle,

    Southwest University of Science and Technology, Mianyang 621010, China;

    3. Xinjiang Technical Institute of Physics and Chemistry, Chinese Academy of Sciences,

    Urumqi 830011, China;

    4. School of Life Science and Engineering, Southwest University of Science and Technology,

    Mianyang 621010, China)

    ZHOU Qing1,2,DONG Faqin2,BIAN Liang3,SUN Shiyong1,2,GUO Yuting2,4,HUO Tingting2

    (1. School of Environment and Resources, Southwest University of Science and Technology,Mianyang 621010, China;

    Molecular simulation has been an important tool in the study of interface interaction due to its visibility at the scale of molecular. To simulate the interaction between montmorillonite and cell membrane, the adsorption of the amino acids clusters and phospholipids on montmorillonite surface was studied by molecular dynamics simulation. Moreover, the heavy metal ions Cr(Ⅲ), Cr(Ⅵ), Cd(Ⅱ), Pb(Ⅱ) was added into the calculation system to acquire its effect on the adsorption process. The results showed that amino acids and phospholipids both had good affinity with montmorillonite. However, amino acids were stronger than phospholipids. Heavy metal ions then played a synergistic role. The electron transition between montmorillonite and the biomolecules was found by the calculation of electronic structure with Density Fuctional Theory. During the process, montmorillonite acted as donor and the biomolecules acceptors.

    montmorillonite; amino acids; phospholipid; molecular simulation; heavy metal ion; electron transition

    1001-9731(2016)04-04129-05

    國家重點基礎研究發(fā)展計劃資助項目(2014CB846003);國家自然科學基金重點資助項目(41130746,41472310)

    2015-12-06

    2016-03-05 通訊作者:董發(fā)勤,E-mail: fqdong@swust.edu.cn

    周 青 (1990-),男,河南信陽人,在讀碩士,師承董發(fā)勤教授,從事環(huán)境礦物學研究。

    X513;O561

    A

    10.3969/j.issn.1001-9731.2016.04.026

    猜你喜歡
    電子密度蒙脫石磷脂
    葛根固斂湯聯(lián)合蒙脫石散治療小兒腹瀉的臨床觀察
    顧及地磁影響的GNSS電離層層析不等像素間距算法*
    大黃酸磷脂復合物及其固體分散體的制備和體內(nèi)藥動學研究
    中成藥(2019年12期)2020-01-04 02:02:24
    不同GPS掩星電離層剖面產(chǎn)品相關(guān)性分析
    測繪通報(2019年11期)2019-12-03 01:47:34
    等離子體電子密度分布信息提取方法研究
    柚皮素磷脂復合物的制備和表征
    中成藥(2018年7期)2018-08-04 06:04:18
    一種適用于電離層電子密度重構(gòu)的AMART算法
    測繪學報(2018年1期)2018-02-27 02:23:07
    辣椒堿磷脂復合凝膠的制備及其藥動學行為
    中成藥(2017年12期)2018-01-19 02:06:31
    白楊素磷脂復合物的制備及其藥動學行為
    中成藥(2017年5期)2017-06-13 13:01:12
    草酸對蒙脫石的作用及機理研究
    性色av乱码一区二区三区2| 欧美亚洲 丝袜 人妻 在线| 久久久久久人人人人人| 婷婷丁香在线五月| 国产麻豆69| 中文字幕人妻丝袜制服| 精品视频人人做人人爽| 国产深夜福利视频在线观看| 美女国产高潮福利片在线看| 精品高清国产在线一区| svipshipincom国产片| 波多野结衣一区麻豆| 丝袜美腿诱惑在线| 国产91精品成人一区二区三区 | 久久天堂一区二区三区四区| 91精品国产国语对白视频| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产亚洲av高清一级| 精品一品国产午夜福利视频| 50天的宝宝边吃奶边哭怎么回事| 国产精品秋霞免费鲁丝片| 欧美亚洲 丝袜 人妻 在线| 一级毛片黄色毛片免费观看视频| 高清视频免费观看一区二区| 一个人免费看片子| 中文字幕制服av| 久久国产精品影院| 99久久99久久久精品蜜桃| 成人影院久久| 十分钟在线观看高清视频www| 日本av免费视频播放| 夫妻性生交免费视频一级片| 精品人妻在线不人妻| 亚洲少妇的诱惑av| 999久久久国产精品视频| 久久人人97超碰香蕉20202| 亚洲精品在线美女| 日本欧美国产在线视频| 午夜福利视频精品| 青春草亚洲视频在线观看| 中文精品一卡2卡3卡4更新| 两性夫妻黄色片| 美女大奶头黄色视频| 欧美日本中文国产一区发布| 日韩欧美一区视频在线观看| 女人高潮潮喷娇喘18禁视频| 午夜福利视频在线观看免费| www.999成人在线观看| 精品福利永久在线观看| 日韩熟女老妇一区二区性免费视频| 欧美日韩成人在线一区二区| 高清av免费在线| 在线观看免费高清a一片| 秋霞在线观看毛片| 在线观看人妻少妇| 19禁男女啪啪无遮挡网站| 男女之事视频高清在线观看 | 午夜久久久在线观看| av一本久久久久| 国产成人精品久久二区二区免费| 91老司机精品| 91九色精品人成在线观看| 老司机靠b影院| 999久久久国产精品视频| 国产又色又爽无遮挡免| 国产高清不卡午夜福利| 日本色播在线视频| 色综合欧美亚洲国产小说| 啦啦啦中文免费视频观看日本| av不卡在线播放| 久久久精品区二区三区| 国产精品麻豆人妻色哟哟久久| 人人妻人人爽人人添夜夜欢视频| 欧美xxⅹ黑人| 欧美在线一区亚洲| 欧美日韩亚洲综合一区二区三区_| 在线观看免费日韩欧美大片| 十分钟在线观看高清视频www| 激情五月婷婷亚洲| 真人做人爱边吃奶动态| 亚洲精品日本国产第一区| 精品少妇一区二区三区视频日本电影| 91字幕亚洲| 国产深夜福利视频在线观看| 好男人电影高清在线观看| 婷婷成人精品国产| 中文字幕亚洲精品专区| 亚洲熟女毛片儿| 五月天丁香电影| 热99国产精品久久久久久7| 大片免费播放器 马上看| 啦啦啦中文免费视频观看日本| 久久ye,这里只有精品| 又粗又硬又长又爽又黄的视频| 中文字幕色久视频| 亚洲成国产人片在线观看| 桃花免费在线播放| 十分钟在线观看高清视频www| 人人妻人人澡人人看| 波多野结衣一区麻豆| 国产黄频视频在线观看| 亚洲一区二区三区欧美精品| 老鸭窝网址在线观看| 在线看a的网站| 狠狠婷婷综合久久久久久88av| 国产熟女欧美一区二区| 亚洲av日韩精品久久久久久密 | 久久国产精品大桥未久av| 在线天堂中文资源库| 宅男免费午夜| 欧美老熟妇乱子伦牲交| 久久精品aⅴ一区二区三区四区| 亚洲伊人色综图| 黑人巨大精品欧美一区二区蜜桃| 久久久久久人人人人人| 热re99久久精品国产66热6| 亚洲欧美激情在线| 丝袜美足系列| 美女福利国产在线| 高潮久久久久久久久久久不卡| 91精品国产国语对白视频| 国产精品一国产av| 国产精品一区二区在线观看99| 久久av网站| 国产av一区二区精品久久| 巨乳人妻的诱惑在线观看| 韩国精品一区二区三区| 在线亚洲精品国产二区图片欧美| 亚洲精品国产区一区二| 另类精品久久| 亚洲熟女毛片儿| 99热全是精品| 少妇粗大呻吟视频| 一边摸一边抽搐一进一出视频| 真人做人爱边吃奶动态| av在线播放精品| 久久热在线av| 国产精品一区二区在线不卡| 国产激情久久老熟女| 91麻豆av在线| 成人三级做爰电影| 国产精品亚洲av一区麻豆| 18禁黄网站禁片午夜丰满| 欧美人与性动交α欧美精品济南到| 亚洲五月婷婷丁香| 国产成人影院久久av| 麻豆av在线久日| 久久久亚洲精品成人影院| 99国产精品一区二区蜜桃av | 亚洲男人天堂网一区| 亚洲精品一二三| 黄片小视频在线播放| 91精品三级在线观看| 捣出白浆h1v1| 国产福利在线免费观看视频| 在线看a的网站| 国产福利在线免费观看视频| 免费高清在线观看日韩| 女人久久www免费人成看片| 久久精品人人爽人人爽视色| 久久久国产精品麻豆| 性高湖久久久久久久久免费观看| www.av在线官网国产| 日本黄色日本黄色录像| 国产三级黄色录像| 女人久久www免费人成看片| 美女高潮到喷水免费观看| 色网站视频免费| 丝袜美腿诱惑在线| 最近最新中文字幕大全免费视频 | 日日摸夜夜添夜夜爱| 9191精品国产免费久久| 欧美日韩国产mv在线观看视频| 美国免费a级毛片| 亚洲自偷自拍图片 自拍| 欧美日韩精品网址| 水蜜桃什么品种好| 亚洲精品一卡2卡三卡4卡5卡 | www.精华液| 爱豆传媒免费全集在线观看| 777米奇影视久久| av一本久久久久| 日本av手机在线免费观看| av在线播放精品| 午夜日韩欧美国产| 午夜福利影视在线免费观看| 五月天丁香电影| 午夜激情av网站| 999精品在线视频| 亚洲少妇的诱惑av| 少妇精品久久久久久久| 亚洲精品久久成人aⅴ小说| 国产又色又爽无遮挡免| www.999成人在线观看| 免费高清在线观看日韩| 一级毛片黄色毛片免费观看视频| 日本a在线网址| 色精品久久人妻99蜜桃| 中国美女看黄片| 欧美日韩一级在线毛片| 免费不卡黄色视频| 大片电影免费在线观看免费| 免费av中文字幕在线| 黑人欧美特级aaaaaa片| 啦啦啦视频在线资源免费观看| 精品国产国语对白av| 中文乱码字字幕精品一区二区三区| 亚洲午夜精品一区,二区,三区| 久久久精品区二区三区| 好男人电影高清在线观看| 亚洲熟女精品中文字幕| 国产精品偷伦视频观看了| 婷婷色综合大香蕉| 国产97色在线日韩免费| 国产精品 欧美亚洲| 黑人欧美特级aaaaaa片| 少妇猛男粗大的猛烈进出视频| 9191精品国产免费久久| 欧美国产精品va在线观看不卡| 777久久人妻少妇嫩草av网站| 亚洲精品国产区一区二| 精品亚洲成a人片在线观看| 国产主播在线观看一区二区 | 婷婷丁香在线五月| 丰满迷人的少妇在线观看| 国产免费一区二区三区四区乱码| 免费日韩欧美在线观看| 丰满饥渴人妻一区二区三| 亚洲伊人色综图| 考比视频在线观看| 亚洲av日韩精品久久久久久密 | 日本色播在线视频| 老司机午夜十八禁免费视频| 熟女少妇亚洲综合色aaa.| 欧美久久黑人一区二区| 咕卡用的链子| 久久国产精品影院| 亚洲图色成人| h视频一区二区三区| 19禁男女啪啪无遮挡网站| 免费观看av网站的网址| 视频在线观看一区二区三区| 视频区图区小说| 亚洲欧洲日产国产| 亚洲色图综合在线观看| 美女国产高潮福利片在线看| 18禁观看日本| 国产一区有黄有色的免费视频| a级片在线免费高清观看视频| 免费高清在线观看日韩| 又大又爽又粗| 国产成人免费观看mmmm| 91老司机精品| 不卡av一区二区三区| 婷婷色av中文字幕| av在线老鸭窝| 亚洲精品美女久久久久99蜜臀 | 脱女人内裤的视频| 免费日韩欧美在线观看| 中文字幕亚洲精品专区| 日韩 欧美 亚洲 中文字幕| 久久鲁丝午夜福利片| 在线看a的网站| svipshipincom国产片| 成年人午夜在线观看视频| 精品亚洲乱码少妇综合久久| 蜜桃在线观看..| 久久国产精品大桥未久av| 欧美日韩一级在线毛片| √禁漫天堂资源中文www| 国产成人精品久久二区二区免费| 久久久久国产精品人妻一区二区| 国产爽快片一区二区三区| 精品亚洲成a人片在线观看| 久久鲁丝午夜福利片| 精品福利永久在线观看| 中文精品一卡2卡3卡4更新| 日日夜夜操网爽| 激情视频va一区二区三区| 精品人妻在线不人妻| 免费女性裸体啪啪无遮挡网站| bbb黄色大片| 人人澡人人妻人| 秋霞在线观看毛片| 精品亚洲成国产av| 看十八女毛片水多多多| 日本午夜av视频| svipshipincom国产片| 成年人午夜在线观看视频| 国产亚洲av片在线观看秒播厂| 国产高清不卡午夜福利| 国产成人av教育| 亚洲精品久久午夜乱码| 亚洲精品一卡2卡三卡4卡5卡 | 91精品伊人久久大香线蕉| 国产精品久久久久久精品古装| 午夜免费观看性视频| 日本a在线网址| 在线av久久热| 日韩免费高清中文字幕av| 老司机靠b影院| 国产91精品成人一区二区三区 | 精品少妇内射三级| 狂野欧美激情性bbbbbb| 国产成人系列免费观看| 日韩熟女老妇一区二区性免费视频| 嫁个100分男人电影在线观看 | 亚洲精品中文字幕在线视频| h视频一区二区三区| 叶爱在线成人免费视频播放| 亚洲一区二区三区欧美精品| 亚洲国产欧美在线一区| 国产av国产精品国产| 国产成人免费无遮挡视频| 美国免费a级毛片| 成年人午夜在线观看视频| 国产亚洲av片在线观看秒播厂| 一二三四社区在线视频社区8| 久久免费观看电影| 日韩av免费高清视频| 成人黄色视频免费在线看| 成人三级做爰电影| 免费日韩欧美在线观看| 男女午夜视频在线观看| 国产在线免费精品| 一区二区三区精品91| 黑人猛操日本美女一级片| 欧美国产精品va在线观看不卡| 97在线人人人人妻| 国产三级黄色录像| 一级片免费观看大全| 大码成人一级视频| 亚洲中文日韩欧美视频| 久久久久久免费高清国产稀缺| 蜜桃在线观看..| 久久精品国产亚洲av高清一级| 免费不卡黄色视频| av电影中文网址| 亚洲图色成人| 国产高清不卡午夜福利| 99热国产这里只有精品6| 亚洲,一卡二卡三卡| 纯流量卡能插随身wifi吗| 日本黄色日本黄色录像| 国产一区二区三区av在线| av在线播放精品| 亚洲av在线观看美女高潮| 中文字幕人妻丝袜制服| 亚洲天堂av无毛| 久久综合国产亚洲精品| 亚洲中文字幕日韩| 日韩,欧美,国产一区二区三区| 最黄视频免费看| 黄色a级毛片大全视频| 国产老妇伦熟女老妇高清| 涩涩av久久男人的天堂| videosex国产| 最新在线观看一区二区三区 | 欧美精品亚洲一区二区| 日韩电影二区| 亚洲人成电影观看| 丁香六月欧美| 老熟女久久久| 在线看a的网站| 激情五月婷婷亚洲| 韩国高清视频一区二区三区| 亚洲人成网站在线观看播放| 老鸭窝网址在线观看| 国产男女内射视频| 麻豆乱淫一区二区| 亚洲av美国av| 成人黄色视频免费在线看| 高清不卡的av网站| 亚洲国产精品国产精品| 叶爱在线成人免费视频播放| 在线av久久热| 一区二区三区乱码不卡18| 波多野结衣av一区二区av| 天堂中文最新版在线下载| 亚洲av男天堂| 九色亚洲精品在线播放| 午夜福利一区二区在线看| 男女边摸边吃奶| 青青草视频在线视频观看| 欧美日韩亚洲高清精品| 精品少妇黑人巨大在线播放| 69精品国产乱码久久久| 国产免费一区二区三区四区乱码| 女人高潮潮喷娇喘18禁视频| 狠狠婷婷综合久久久久久88av| 一级片'在线观看视频| 狂野欧美激情性bbbbbb| 久久天堂一区二区三区四区| 国产精品偷伦视频观看了| 黄色怎么调成土黄色| 午夜精品国产一区二区电影| 亚洲av电影在线进入| 各种免费的搞黄视频| 一边摸一边抽搐一进一出视频| 日韩精品免费视频一区二区三区| 免费看十八禁软件| 我的亚洲天堂| 亚洲图色成人| 成年美女黄网站色视频大全免费| www.熟女人妻精品国产| 久久久精品区二区三区| 欧美激情极品国产一区二区三区| 操美女的视频在线观看| 手机成人av网站| 熟女av电影| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 午夜影院在线不卡| 久久性视频一级片| 9191精品国产免费久久| 成年动漫av网址| 亚洲熟女精品中文字幕| 丰满少妇做爰视频| 一区福利在线观看| 国产黄频视频在线观看| 99国产精品一区二区三区| 亚洲国产精品999| 精品国产一区二区久久| 久久精品亚洲熟妇少妇任你| 亚洲第一青青草原| 极品人妻少妇av视频| 下体分泌物呈黄色| 国产一区二区三区综合在线观看| 狂野欧美激情性xxxx| 又紧又爽又黄一区二区| 国产精品九九99| av又黄又爽大尺度在线免费看| 国产黄色视频一区二区在线观看| 国产一区二区激情短视频 | 99精品久久久久人妻精品| 热99国产精品久久久久久7| 男女床上黄色一级片免费看| 日韩一卡2卡3卡4卡2021年| 欧美乱码精品一区二区三区| 国产午夜精品一二区理论片| 男人操女人黄网站| 亚洲欧美清纯卡通| 热99久久久久精品小说推荐| 视频区图区小说| 2021少妇久久久久久久久久久| 考比视频在线观看| 午夜老司机福利片| 777久久人妻少妇嫩草av网站| 日韩熟女老妇一区二区性免费视频| 一级片'在线观看视频| 久久久精品区二区三区| 少妇裸体淫交视频免费看高清 | 欧美另类一区| 日本av手机在线免费观看| 少妇精品久久久久久久| 日本av免费视频播放| 超碰成人久久| 成年人免费黄色播放视频| 一级片'在线观看视频| 亚洲精品一卡2卡三卡4卡5卡 | 精品亚洲乱码少妇综合久久| 亚洲人成电影观看| 丝袜美足系列| 久热这里只有精品99| bbb黄色大片| 一级黄色大片毛片| 黑丝袜美女国产一区| 中文字幕最新亚洲高清| 十八禁高潮呻吟视频| 最新的欧美精品一区二区| 亚洲国产欧美网| 亚洲人成网站在线观看播放| 午夜老司机福利片| 性色av一级| 中文字幕人妻丝袜制服| 久久久久久久久免费视频了| tube8黄色片| 欧美日韩视频精品一区| www日本在线高清视频| 黄色视频不卡| 夫妻午夜视频| 丁香六月天网| 十分钟在线观看高清视频www| 老熟女久久久| 婷婷色综合大香蕉| 精品国产超薄肉色丝袜足j| 亚洲一区二区三区欧美精品| 99久久99久久久精品蜜桃| 国产不卡av网站在线观看| e午夜精品久久久久久久| 亚洲成人免费电影在线观看 | 91老司机精品| 日韩 欧美 亚洲 中文字幕| 亚洲人成电影免费在线| 亚洲精品一卡2卡三卡4卡5卡 | 久久久久久久大尺度免费视频| 久久精品国产亚洲av高清一级| 久久青草综合色| 久久亚洲精品不卡| 99久久精品国产亚洲精品| 黄网站色视频无遮挡免费观看| 免费观看av网站的网址| 久久国产亚洲av麻豆专区| 尾随美女入室| 国产男人的电影天堂91| 999久久久国产精品视频| 久久免费观看电影| 精品久久久精品久久久| 精品少妇黑人巨大在线播放| 国产亚洲欧美精品永久| 中文字幕亚洲精品专区| 国产黄色免费在线视频| 老司机深夜福利视频在线观看 | a级片在线免费高清观看视频| e午夜精品久久久久久久| 妹子高潮喷水视频| 国产精品久久久久久人妻精品电影 | 久久久久视频综合| 999久久久国产精品视频| 十八禁人妻一区二区| 2021少妇久久久久久久久久久| 国产在线免费精品| 青春草亚洲视频在线观看| 国产在线观看jvid| 女警被强在线播放| 国产精品一二三区在线看| 国产女主播在线喷水免费视频网站| 精品少妇内射三级| 99re6热这里在线精品视频| 亚洲国产中文字幕在线视频| 日本五十路高清| 亚洲天堂av无毛| 国产亚洲av高清不卡| 国产福利在线免费观看视频| 国产成人精品久久二区二区免费| 高清视频免费观看一区二区| 国产男女内射视频| 亚洲国产毛片av蜜桃av| 亚洲欧美成人综合另类久久久| 韩国高清视频一区二区三区| 伊人亚洲综合成人网| 又黄又粗又硬又大视频| videosex国产| 国产99久久九九免费精品| 尾随美女入室| 久久久久精品国产欧美久久久 | 亚洲欧美一区二区三区久久| 亚洲成av片中文字幕在线观看| 男女高潮啪啪啪动态图| 亚洲国产看品久久| 首页视频小说图片口味搜索 | 久久性视频一级片| 亚洲精品成人av观看孕妇| 18在线观看网站| 99热全是精品| 亚洲天堂av无毛| 精品亚洲乱码少妇综合久久| 性少妇av在线| 免费黄频网站在线观看国产| 亚洲国产欧美在线一区| 无遮挡黄片免费观看| 欧美日韩一级在线毛片| 晚上一个人看的免费电影| 日本av免费视频播放| 一边摸一边做爽爽视频免费| 满18在线观看网站| 欧美 亚洲 国产 日韩一| 大陆偷拍与自拍| 午夜日韩欧美国产| 啦啦啦 在线观看视频| 丰满迷人的少妇在线观看| tube8黄色片| av片东京热男人的天堂| 乱人伦中国视频| 精品久久久久久电影网| 亚洲图色成人| 日韩制服骚丝袜av| 日韩制服丝袜自拍偷拍| 人人妻人人澡人人看| 国产欧美日韩一区二区三区在线| 欧美人与性动交α欧美精品济南到| 国产无遮挡羞羞视频在线观看| 大香蕉久久成人网| 国产日韩欧美视频二区| 免费高清在线观看视频在线观看| 亚洲欧美激情在线| 国产精品成人在线| 国产日韩欧美视频二区| 亚洲一区二区三区欧美精品| 99热全是精品| 视频区欧美日本亚洲| 久久ye,这里只有精品| 午夜福利,免费看| 欧美久久黑人一区二区| 久久av网站| 成人18禁高潮啪啪吃奶动态图| 人人妻人人爽人人添夜夜欢视频| 91老司机精品| 亚洲中文av在线| 好男人视频免费观看在线| 国产高清不卡午夜福利| 美国免费a级毛片| 狠狠婷婷综合久久久久久88av| 夜夜骑夜夜射夜夜干| 国产免费福利视频在线观看| 在线 av 中文字幕| 丝袜人妻中文字幕| 啦啦啦啦在线视频资源| 69精品国产乱码久久久| 亚洲精品国产av蜜桃| 日韩一卡2卡3卡4卡2021年| 亚洲成人免费av在线播放| 亚洲综合色网址| 宅男免费午夜| 日日夜夜操网爽|