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

大豆皂苷Ⅱ與人血清白蛋白相互作用機制研究

2016-09-14 12:24:52桑尚源光翠娥食品科學與技術國家重點實驗室江南大學江蘇無錫214122
食品工業科技 2016年8期
關鍵詞:大豆

桑尚源,光翠娥,江 波(食品科學與技術國家重點實驗室,江南大學,江蘇無錫214122)

大豆皂苷Ⅱ與人血清白蛋白相互作用機制研究

桑尚源,光翠娥*,江 波
(食品科學與技術國家重點實驗室,江南大學,江蘇無錫214122)

為研究吸收進入人體血循環的大豆皂苷Ⅱ與人血清白蛋白(HSA)相互作用機理。采用圓二色光譜、內源熒光光譜和分子模擬方法檢測表征皂苷Ⅱ與HSA之間的相互作用。圓二色光譜結果顯示,大豆皂苷Ⅱ與HSA的濃度比例為1∶2或1∶4時,HSA的二級結構α螺旋含量上升,β折疊含量下降。熒光光譜顯示,大豆皂苷Ⅱ對HSA的立體構象有影響。分子模擬顯示它們之間的結合作用力以氫鍵和疏水相互作用為主,HSA與皂苷糖基鏈中葡萄糖醛酸基(GlcA)和鼠李糖基(Rha)部分形成氫鍵作用的氨基酸殘基包括Glu153、Ser192、Gln196、Glu450、Asp451和Ser454等,與皂苷配基三萜烯部分形成疏水相互作用的氨基酸殘基主要有Ala、Leu、Val和Trp等。

大豆皂苷Ⅱ,人血清白蛋白,相互作用,圓二色光譜,分子模擬

大豆皂苷一般被認為是豆制品中的一種功能成分,其生物活性有很多。有益的作用包括:降低心血管疾病發生率、預防癌癥、保護肝臟和消炎等[1]。大豆皂苷Ⅱ是豆制品中皂苷的主要存在形式之一[2],化學結構見圖1。其具有多種生物活性,如抑制人體腎素活性[3]、抗皰疹病毒(HSV-1)[4]和作為免疫佐劑[5-6]。對其進行吸收、代謝與分布等生物利用度方面的研究,為今后研發出相關保健食品提供重要參考價值。

血清白蛋白(serum albumin,SA)是動物血漿中極易純化且含量又最豐富的蛋白質(0.6 mmol/L),因而為研究蛋白質在體內的生物學功能、理化性質、臨床應用和體內代謝等方面提供了較好的模式蛋白質。人血清白蛋白(HSA)由585個氨基酸殘基組成,分子量在67000 u左右。HSA的二級結構以α螺旋結構為主,在HSA的氨基酸全序列中含有35個半胱氨酸(Cys),其中只有一個自由的Cys未形成二硫鍵,其他17對Cys形成二硫鍵。HSA晶體的三維結構模型為“心型”,可分為三個結構域Ⅰ、Ⅱ和Ⅲ,每個結構域又分為亞結構域A和B,且每個亞結構域是由數個α螺旋鏈以圓筒狀組成,圓筒狀的亞結構域A和B又以槽口相對方向組合在一起。其中第1~195位氨基酸殘基組成結構域Ⅰ,第196~383位氨基酸組成結構域Ⅱ,余下的氨基酸組成結構域Ⅲ[7-8]。HSA具有結合和運輸內源性及外源性物質的功能,深度影響一些吸收進入人體的外源活性物質在體內的運送、分布、代謝及消除等過程。為此本文模擬生理條件,采用圓二色光譜、內源熒光光譜和分子模擬方法表征皂苷Ⅱ與HSA之間的相互作用機制,為將來開發含有高生物利用度大豆皂苷的保健食品提供理論依據。

圖1 大豆皂苷Ⅱ的化學結構Fig.1 Structure of soyasaponinⅡ

1 材料與方法

1.1 材料與儀器

大豆皂苷Ⅱ(純度91.7%) 美國ChromaDex公司;人血清白蛋白(HSA) Sigma公司;三羥甲基氨甲烷(Tris)、氯化鈉 國藥集團試劑公司,分析純;Tris-HCl緩沖液 500 mL體積的0.1 mol/L Tris溶液中加入約420 mL體積的0.1 mol/L的鹽酸,調節pH至7.4,加入5.844 g氯化鈉,定容至1000 mL,4℃低溫保存待用。

MOS-450 AF-CD圓二色光譜儀 法國Bio-Logic公司;FLUORMAX-4熒光光譜儀 美國HORIBA公司;軟件Discover Studio 2.5(DS)、LigPlot+分析軟件美國Accelrys軟件公司。

1.2 實驗方法

1.2.1 人血清白蛋白的遠紫外圓二色光譜檢測 用Tris緩沖液配制HSA與大豆皂苷Ⅱ不同比例的混合液。HSA的最終濃度為1 μmol/L,大豆皂苷Ⅱ的終濃度依次為0、0.5、1、2、4 μmol/L,混合靜置30 min。分別吸取200 μL Tris-HCl緩沖溶液(pH7.4)和樣品注入2 mm比色皿中測定圓二色光譜,每次測量扣除對應皂苷濃度的Tris-HCl緩沖溶液的背景干擾。室溫下(293 K)光譜測定條件是200~250 nm波長范圍掃描,掃描速度為20 nm/min,響應時間為2 s。所得數據使用Selcon3方法計算得出血清白蛋白的二級結構含量。

1.2.2 人血清白蛋白的內源熒光光譜檢測 以295 nm為激發波長,狹縫為5 nm,室溫下(293 K),在300~450 nm波長范圍內掃描HSA樣品(與方法1.2.1相同的樣品)的發射光得到內源熒光光譜圖,重復掃描3次取平均值。空白為對應皂苷濃度的Tris-HCl緩沖溶液,熒光強度值為扣除空白后的相對值。

1.2.3 分子模擬 在蛋白質數據庫(PDB)中獲取HSA晶體的三維結構模型數據文件,登錄號為1h9z。使用DS2.5軟件刪除其中的配體小分子和水分子,并補入缺失的氫原子,得到HSA晶體的完整三維結構。大豆皂苷Ⅱ(soyasaponinⅡ)的立體結構數據文件從網站www.chemical-book.com中搜索下載,采用PRODRG網站(http://davapc1.bioch.dundee.ac.uk/cgibin/prodrg)對大豆皂苷Ⅱ的立體結構進行加電荷和構象優化等預處理[9]。使用DS軟件中的Libdock模塊進行剛性分子對接。以HSA的Trp214殘基為中心,半徑為13 ?的球體范圍定義為對接區域,對接程序參數設置如下:Docking Tolerance=0.25、Numbers of Hotspots= 100、Docking Preference=High Quality、Conformation Method=FAST、Parallel Processing=False。從對接結果大豆皂苷Ⅱ構象集中選取打分函數LibDockscore得分最高者,并將其對接復合物的文件導入Ligplot+軟件[10]中,得到蛋白與皂苷分子間相互作用網,并進一步分析得到氫鍵和疏水相互作用位點。

1.2.4 統計學分析 實驗所得數據都是平均值±標準差表示,并且用SPSS統計學軟件進行t檢驗,p<0.05表示差異有顯著性。

2 結果與討論

2.1 不同大豆皂苷Ⅱ濃度對HSA二級結構的影響

包含生色團的蛋白質大分子由于其構象的不對稱性,可引起左右圓偏振光具有不同光吸收的現象,即為圓二色性(circular dichroism,CD),這種光被吸收的程度與波長之間關系稱為圓二色光譜。遠紫外區190~250 nm波長范圍的圓二色光譜可反映蛋白質主鏈的構象,α螺旋結構在靠近192 nm有一正的譜帶,在遠紫外208 nm和222 nm處表現出兩個負的特征肩峰譜帶;β折疊的CD光譜在216 nm有一負譜帶[11]。

室溫下(293 K),HSA在不同大豆皂苷Ⅱ濃度下,進行遠紫外200~250 nm波長范圍掃描得到圓二色光譜,結果如圖2所示。從圖2中可以看出,HSA的CD光譜圖在200~220 nm波長范圍有一段負譜帶,說明其二級結構形式以α螺旋和β折疊為主。不同濃度的大豆皂苷Ⅱ未改變HSA的CD光譜峰型和峰位,但改變峰強度,表明大豆皂苷Ⅱ的存在引起HSA的二

級結構發生變化。

圖2 HSA-大豆皂苷Ⅱ混合溶液的遠紫外CD光譜圖(pH=7.4,T=293 K)Fig.2 Far-UV CD spectra of the HSA-soyasaponinⅡsystem (pH=7.4,T=293 K)

圓二色光譜數據通過Selcon3模型計算方法,計算得到HSA二級結構含量,如圖3所示。從圖3中可知,HSA與大豆皂苷Ⅱ的濃度比例為1∶1時,α螺旋和β折疊的含量未變化(p>0.05)。而在1∶0.5、1∶2和1∶4時,HSA的α螺旋含量由35.3%增加至46.4%、67.7%和56.5%,而β折疊含量由11.5%下降至8.8%、6.2%和6.7%(p<0.05)。大豆皂苷Ⅱ與HSA的結合可能改變了HSA的氫鍵骨架結構,導致HSA的二級結構從β折疊轉變為α螺旋。

圖3 不同蛋白-大豆皂苷Ⅱ比例下人血清白蛋白二級結構含量變化圖(pH=7.4,T=293 K)Fig.3 Effect of the HSA-soyasaponinⅡconcentration ratio on secondary structure constituents of HSA(pH=7.4,T=293 K)

2.2 不同大豆皂苷Ⅱ濃度對HSA立體構象的影響

當以295 nm為激發波長時,蛋白的內源熒光基團主要是由色氨酸(Trp)產生熒光發射光譜。加入一定量的小分子后,HSA與配體小分子相互作用后,受體蛋白的三維結構改變,或者蛋白中的色氨酸殘基附近微環境極性發生改變,進而表現出熒光光譜發生變化。HSA有唯一的Trp214,在蛋白中心疏水腔內部[7],即在亞結構域ⅡA與ⅢA之間,見圖5。

實驗得到不同蛋白-大豆皂苷Ⅱ濃度比例下的HSA內源熒光光譜,結果見表1。HSA蛋白在334 nm出現最大吸收峰,并隨大豆皂苷Ⅱ濃度增加,熒光光譜峰位發生藍移(從334~332 nm)。可能原因是HSA的三維結構改變,從而色氨酸殘基附近微環境極性減弱。從大豆皂苷Ⅱ對HSA內源熒光強度猝滅作用分析,大豆皂苷Ⅱ降低了HSA熒光特征峰的強度(比例為1∶4時,相對峰強度為98%),這是因為大豆皂苷Ⅱ結合在Trp214(位于HSA疏水腔)附近,影響了其微環境極性。

表1 不同蛋白-大豆皂苷Ⅱ比例下對熒光譜峰位及峰強度的影響Table1 Effect of different concentrations of soyasaponinⅡon fluorescence spectra profiles of HSA

2.3 大豆皂苷Ⅱ與HSA相互作用類型研究

為確定皂苷與HSA其他氨基酸殘基主要作用類型和相互作用位點,隨后進行計算機模擬分子對接實驗。DS工作站中Libdock模塊進行分子對接后的結果主要是從Conf Number、Absolute Energy、Relative Energy、LibDockscore、Poses、相對RMSD值等幾個方面對其進行評價。其中,對接吻合度打分函數LibDockscore是對幾何形狀、化學環境、能量的綜合評價參數,同時LibDockscore也是根據配體與受體對接時的對接模式(Conf Number)、親和力(Absolute Energy)和相對能量值(Relative Energy)等方面進行綜合打分。因此采用打分函數LibDockscore來評價對接結果較為合理,一般認為LibDockscore>100表示配體與受體能較好的結合,其值越大代表對接產生的復合物構象越穩定。Poses是與受體蛋白對接時所得配體的構象數。

大豆皂苷Ⅱ與HSA分子對接后,所得大豆皂苷Ⅱ立體構象總數(Poses)為8,其中LibDockscore所得最高分為187.092,表明兩者可穩定結合。選取打分函數(LibDockscore)評分最高的大豆皂苷Ⅱ構象,制作對接三維示意圖,見圖4。從圖4中可看出,函數LibDockscore評分值最高的大豆皂苷Ⅱ在HSA疏水腔內Trp214殘基附近,大豆皂苷Ⅱ與Trp214之間也許存在有利于結合的相互作用力。

圖4 大豆皂苷Ⅱ與HSA的對接結果三維示意圖Fig.4 The molecular docking schematic 3D-diagram of soyasaponinⅡand HSA

對圖4中復合物最佳構象結構運用Ligplot+軟件,得到分子間相互作用網二維圖,結果見圖5。Ligplot+軟件是一個用來分析配體和受體間相互作用的計算程序[10],它通過HBPLUS程序計算疏水和氫鍵作用,以二維圖形式表示出來,有效避免化學鍵交叉和重疊,比三維圖利于觀察和分析統計[12]。在圖5中,統計受體HSA中與配體皂苷小分子相互作用力類型主要為氫鍵和疏水相互作用,還可利用圖5對分子間作用位點進行詳細分析。

圖5 大豆皂苷Ⅱ與HSA分子間相互作用網二維圖分析Fig.5 Molecular interaction 2D-analysis between HSA and soyasaponinⅡ

使用Ligplot+軟件得到分子間相互作用網二維圖(圖5),現將體系中較為穩定的氫鍵作用統計在表2中。從表2中可看出,在HSA蛋白氨基酸殘基中,Glu153、Ser192、Gln196、Glu450、Asp451、Ser454等六個殘基與大豆皂苷Ⅱ有氫鍵相互作用。根據以上氫鍵作用分析的結果,使用DS工作站分析顯示出氫鍵作用的三維效果圖,得到圖6。在圖6中,大豆皂苷Ⅱ的第一個糖基(葡萄糖醛酸基,GlcA)、第三個糖基(鼠李糖基,Rha)與HSA存在氫鍵作用。

表2 大豆皂苷Ⅱ與人血清白蛋白分子間氫鍵分析Table2 Hydrogen bond analysis between HSA and soyasaponinⅡ

同樣地,統計分析二維圖(圖5)可以得到HSA與大豆皂苷Ⅱ的疏水相互作用詳細信息,參與疏水作用的氨基酸殘基主要有Lys195、Lys199、Ser202、Phe206、Ala210、Phe211、Trp214、His242、Leu347、Ser480、Leu481、Val482等。顯然HSA中具有內源熒光性的Trp214殘基參與了疏水相互作用。而且熒光檢測結果表2顯示,大豆皂苷Ⅱ可以減弱HSA中的Trp內源熒光強度。這也許是因為大豆皂苷Ⅱ與HSA中的Trp214殘基距離較近,發生了能量轉移。劉璐莎等[13]研究了土貝母皂苷Ⅱ與HSA的相互作用,發現土貝母皂苷Ⅱ與HSA中的Trp214殘基有疏水相互作用,且根據能量轉移原理計算出兩者的距離為4.95 nm,作者推斷受波長295 nm入射光激發的HSA能量向土貝母皂苷Ⅱ轉移,這也是造成土貝母皂苷Ⅱ對HSA熒光猝滅結果的原因之一。

圖6 HSA和大豆皂苷Ⅱ分子間氫鍵作用三維圖Fig.6 The schematic 3D-diagram of H-bonds between HSA and soyasaponinⅡ

在統計分析參與疏水相互作用的氨基酸殘基后,將這些氨基酸殘基通過DS 2.5軟件以三維圖形式顯示出來,見圖7。從圖7中可看出,參與疏水作用的氨基酸在HSA內部組成了一個疏水性腔體(即圖7中代表靜電中性的白色部分),大豆皂苷Ⅱ的強疏水性配基齊墩果烷伸進HSA的此疏水腔內,兩者產生了疏水相互作用。

圖7 人血清白蛋白和大豆皂苷Ⅱ分子間疏水作用三維圖Fig.7 The 3D-schematic diagrams of hydrophobic interaction between HSA and soyasaponinⅡ

3 結論

通過圓二色光譜、內源熒光光譜和分子建模方法研究了大豆皂苷Ⅱ與HSA之間的相互作用。圓二色光譜顯示大豆皂苷Ⅱ對HSA的二級結構含量產生重要影響,在HSA與大豆皂苷Ⅱ比例為1∶0.5、1∶2和1∶4時,HSA的α螺旋結構含量增加,而β折疊含量下降。HSA的內源熒光光譜表明蛋白與大豆皂苷Ⅱ結合后使其立體構象發生變化。計算機模擬分子對接表明

大豆皂苷Ⅱ與HSA的結合穩定性相對較好,主要作用力以氫鍵和疏水相互作用為主。HSA與皂苷糖基鏈GlcA和Rha部分形成氫鍵作用的氨基酸殘基包括Glu153、Ser192、Gln196、Glu450、Asp451和Ser454等,HSA與皂苷配基三萜烯部分形成疏水相互作用的氨基酸殘基主要有Ala、Leu、Val和Trp等。實驗結果為將來開發大豆皂苷保健食品提供了重要參考價值。

[1]Guang CE,Chen J,Sang SY,et al.Biological functionality of soyasaponins and soyasapogenols[J].Journal of Agricultural and Food Chemistry,2014,62(33):8247-8255.

[2]Hu J,Lee SO,Hendrich S,et al.Quantification of the group B soyasaponins by high-performance liquid chromatography[J].Journal of Agricultural and Food Chemistry,2002,50(9):2587-2594.

[3]Takahashi S,Hori K,Shinbo M,et al.Isolation of human renin inhibitor from soybean:soyasaponin I is the novel human renin inhibitor in soybean[J].Bioscience Biotechnology and Biochemistry,2008,72(12):3232-3236.

[4]Kinjo J,Yokomizo K,Hirakawa T,et al.Anti-herpes virus activity of fabaceous triterpenoidal saponins[J].Biological& Pharmaceutical Bulletin,2000,23(7):887-889.

[5]Oda K,Matsuda H,Murakami T,et al.Adjuvant and haemolytic activities of 47 saponins derived from medicinal and food plants [J].Biological Chemistry,2000,381(1):67-74.

[6]Oda K,Matsuda H,Murakami T,et al.Relationship between adjuvant activity and amphipathic structure of soyasaponins[J].Vaccine,2003,21(17-18):2145-2151.

[7]Majorek KA,Porebski PJ,Dayal A,et al.Structural and immunologic characterization of bovine,horse,and rabbit serum albumins[J].Molecular Immunology,2012,52(3-4):174-182.

[8]Sugio S,Kashima A,Mochizuki S,et al.Crystal structure of human serum albumin at 2.5 angstrom resolution[J].Protein Engineering,1999,12(6):439-446.

[9]張海玲,光翠娥,江波,等.皂苷類似物與腎素的分子對接和結合能分析[J].食品與生物技術學報,2014,33(10):1056-1062.

[10]Wallace A C,Laskowski R A,Thornton J M.Ligplot-a program to generate schematic diagrams of protein ligand interactions[J].Protein Engineering,1995,8(2):127-134.

[11]沈星燦,梁宏,何錫文,等.圓二色光譜分析蛋白質構象的方法及研究進展[J].分析化學,2004,32(3):388-394.

[12]Laskowski RA,Swindells MB.LigPlot+:Multiple ligandprotein interaction diagrams for drug discovery[J].Journal of Chemical Information and Modeling,2011,51(10):2778-2786.

[13]劉璐莎,樊君,胡春梅,等.土貝母皂苷II與人血清白蛋白相互作用機制的光譜研究[J].化學學報,2011,69(21):2589-2596.

[3]皇才進,韓魯佳,劉賢,等.基于近紅外光譜技術的秸稈工業分析[J].光譜學與光譜分析,2009,29(4):960-963.

[4]蔣錦鋒,李莉,趙明月.應用近紅外檢測技術快速測定煙葉主要化學成分[J].中國煙草學報,2006,12(2):8-12.

[5]Ozdemir D,Ozturk B.Near infrared spectroscopic determination of olive oil Adulteration with sunflower and corn oil[J].Journal of Food and Drug Analysis,2007,15(1):40-47.

[6]李楠,許韻華,宋雯雯,等.利用近紅外光譜技術快速檢測大豆氨基酸含量[J]植物遺傳資源學報,2012,13(6):1037-1044.

[7]齊龍,朱克衛,馬旭,等.近紅外光譜分析技術在大米檢測中的應用[J].農機化研究,2011(7):18-22.

[8]王傳梁,陳坤杰.基于近紅外漫反射技術的大米脂肪含量的研究[J].糧油加工,2007(2):62-64.

[9]GB/T 5497-1985,《糧食、油料檢驗水分測定法》[S].

[10]GB/T 5511-2008《谷物和豆類氮含量測定和粗蛋白質含量計算凱氏法》[S].

[11]GB/T 15684-1995《谷物制品脂肪酸值測定法》[S].

[12]GB/T 15683-2008《大米直鏈淀粉含量的測定》[S].

[13]嚴衍祿,陳斌,朱大洲,等.近紅外光譜分析的原理、技術與應用[M].北京:中國輕工業出版社,2013:162-168.

[14]范維燕,邢邯,林家永,等.近紅外光譜快速測定稻谷水分含量的研究[J].糧油食品科技,2008,16(5):19-52.

[15]Bagchi T B,Sharma S,Chattopadhyay K.Development of NIRS models to predict protein and amylose content of brown rice and proximate compositions of rice bran[J].Food Chemistry,2016,191(15):21-27.

[16]范維燕,邢邯,林家永,等.稻谷脂肪酸值近紅外光譜快速測定技術研究[J].食品科學,2009,30(24):347-350.

[17]Xie L H,Tang S Q,Chen N,et al.Optimisation of nearinfrared reflectance model in measuring protein and amylose content of rice flour[J].Food Chemistry,2014,142(1):92-100.

Interaction between soyasaponinⅡ and human serum albumin

SANG Shang-yuan,GUANG Cui-e*,JIANG Bo
(The State Key Laboratory of Food Science and Technology,Jiangnan University,Wuxi 214122,China)

Study on the interaction of the absorbed soyasaponinⅡ with human serum albumin(HSA)in the circulatory system was performed.The binding of soyasaponinⅡ to the protein was characterized by circular dichroism spectrometry,intrinsic fluorescence spectrometry and molecular.CD spectra displayed that with the concentration ratios of 1∶2 and 1∶4,alpha helical contents of HSA increased and the beta sheet contents decreased.Fluorescence spectra showed that the 3D structure of albumin after the binding of soyasaponinⅡ was changed.The molecular modeling revealed that the binding forces of soyasaponinⅡ to HSA were mainly hydrogen bond and hydrophobic interactions.GlcA and Rha parts of soyasaponin sugar chains and amino acid residues,including Glu153,Ser192,Gln196,Glu450,Asp451 and Ser454,formed hydrogen bonds.The aglycone of soyasaponinⅡ had the hydrophobic interaction with HSA,amino acid residues of which were mainly Leu,Val,Ala and Trp.

soyasaponinⅡ;HSA;interaction;circular dichroism;molecular modeling

TS201.1

A

1002-0306(2016)08-0130-05

10.13386/j.issn1002-0306.2016.08.018

2015-07-27

桑尚源(1990-),男,碩士,研究方向:食品營養與功能因子,E-mail:shangyuan.sang@qq.com。

*通訊作者:光翠娥(1976-),女,博士,副教授,研究方向:食品營養與功能因子,E-mail:guang1226@hotmail.com。

國家自然科學基金資助項目(31201289);食品科學與技術國家重點實驗室創新探索課題(SKLF-ZZA-201503)。

猜你喜歡
大豆
2022年國審大豆新品種
注意防治大豆點蜂緣蝽
今日農業(2022年16期)2022-11-09 23:18:44
大豆帶狀種植技術如何提升我國大豆產量
從大豆種植面積增長看我國糧食安全
巴西大豆播種順利
今日農業(2021年20期)2021-11-26 01:23:56
大豆的營養成分及其保健作用
今日農業(2021年14期)2021-10-14 08:35:34
用大豆“炸沉”軍艦
中國大豆供需平衡表
大豆農場
大豆搖鈴的夢(外三首)
讀寫算(中)(2015年10期)2015-11-07 07:24:12
主站蜘蛛池模板: 中文字幕啪啪| 日本妇乱子伦视频| 蜜臀AV在线播放| 日韩A级毛片一区二区三区| 亚洲精品少妇熟女| 国产69精品久久| 99视频在线观看免费| 国产精品视频白浆免费视频| 香蕉久久国产超碰青草| 中文字幕 91| 国产成人综合日韩精品无码首页| 亚洲欧美综合在线观看| 国产va免费精品观看| 18禁黄无遮挡免费动漫网站| 高清无码手机在线观看| 岛国精品一区免费视频在线观看| 国产精品网址在线观看你懂的| 国产成人在线小视频| 精品久久久久无码| 亚洲成肉网| 久久九九热视频| 国产v精品成人免费视频71pao| 亚洲色图综合在线| 国产精品网址你懂的| 亚洲国模精品一区| 欧日韩在线不卡视频| 老色鬼久久亚洲AV综合| 国产浮力第一页永久地址| aa级毛片毛片免费观看久| 欧美成人看片一区二区三区 | 五月激情综合网| 热99re99首页精品亚洲五月天| 中文国产成人久久精品小说| 日韩高清中文字幕| 欧美在线国产| 啪啪啪亚洲无码| 91精品国产丝袜| 精品国产欧美精品v| 怡春院欧美一区二区三区免费| 国产精彩视频在线观看| 五月婷婷丁香色| 久久一日本道色综合久久| 久久五月视频| 免费激情网址| 日韩欧美中文字幕在线韩免费| 亚洲欧美日韩中文字幕一区二区三区| 国内a级毛片| 99福利视频导航| 亚洲丝袜第一页| 日本五区在线不卡精品| 激情国产精品一区| 日本精品视频| 午夜精品影院| 九色在线视频导航91| 国产亚洲欧美在线人成aaaa| 99久久亚洲精品影院| 中文字幕无码中文字幕有码在线 | 亚洲综合在线最大成人| 国产精品爆乳99久久| 中文字幕调教一区二区视频| 波多野结衣一级毛片| 国产在线一区二区视频| 亚洲精品777| 亚洲Aⅴ无码专区在线观看q| 在线中文字幕日韩| 国内熟女少妇一线天| 午夜毛片免费观看视频 | 国产成人久久777777| 2019年国产精品自拍不卡| 成人日韩视频| 国产你懂得| 日韩欧美中文| 国产黄在线免费观看| 亚洲色图欧美视频| 呦女亚洲一区精品| 丝袜美女被出水视频一区| 亚洲美女久久| 国产人碰人摸人爱免费视频| 日韩免费毛片视频| 亚洲国产天堂久久综合| 国产剧情无码视频在线观看| 色欲不卡无码一区二区|