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

    小立碗蘚GPX家族的功能分化與催化特性研究

    2015-03-03 02:37:16楊志靈楊海靈
    西北植物學報 2015年1期

    王 煒,楊志靈,楊海靈

    (北京林業(yè)大學 生物科學與技術學院,北京 100083)

    ?

    小立碗蘚GPX家族的功能分化與催化特性研究

    王煒,楊志靈,楊海靈*

    (北京林業(yè)大學 生物科學與技術學院,北京 100083)

    摘要:谷胱甘肽過氧化物酶(GPX)在植物抵抗氧化脅迫中發(fā)揮重要作用。該研究從小立碗蘚(Physcomitrellapatens)基因組中挖掘到3個GPX基因,分別命名為PpGPX1、PpGPX2和PpGPX3。其中PpGPX1和PpGPX3只含有1個外顯子,而PpGPX2含有6個外顯子。表達模式分析發(fā)現(xiàn)PpGPX1和PpGPX2在檢測的所有條件下均表達,而PpGPX3在檢測的所有條件下均不表達。蛋白亞細胞定位分析發(fā)現(xiàn),PpGPX1蛋白定位在細胞質,而PpGPX2蛋白定位在葉綠體。在大腸桿菌中表達并純化了PpGPX1和PpGPX2蛋白,酶學性質分析發(fā)現(xiàn),PpGPX1和PpGPX2蛋白均只能利用Trx電子供體系統(tǒng),而不能利用GSH電子供體系統(tǒng);PpGPX2蛋白對過氧化物底物的催化活性和催化效率均高于PpGPX1。基因結構、表達模式、亞細胞定位和蛋白酶學性質的差異預示小立碗蘚GPX基因家族成員發(fā)生了功能分化,將PpGPX2蛋白的Pro158、Phe167和Phe172氨基酸殘基均突變?yōu)锳la,發(fā)現(xiàn)突變體蛋白對底物催化活性降低,說明這3個氨基酸位點對PpGPX2蛋白具有重要催化活性。

    關鍵詞:谷胱甘肽過氧化物酶;小立碗蘚;功能分化;亞細胞定位;酶學性質分析

    植物在受到干旱、低溫、高鹽、除草劑、病菌侵害等生物及非生物脅迫時,體內的活性氧(reactive oxygen species,ROS)含量會明顯上升,過多的活性氧積累會對植物的核酸、脂質、蛋白等生物大分子造成損傷[1]。植物在長期進化過程中,產生了一系列的酶促及非酶促機制來清除體內過剩的活性氧,例如抗壞血酸過氧化物酶(APX)、過氧化氫酶(CAT)、超氧化物歧化酶(SOD)和谷胱甘肽過氧化物酶(glutathione peroxidase,GPX)等酶促反應系統(tǒng),及谷胱甘肽、類胡蘿卜素、維生素E等非酶促反應系統(tǒng)[2]。谷胱甘肽過氧化物酶作為一種重要的活性氧清除劑,能夠直接與過氧化氫(H2O2)、叔丁基過氧化氫(tBOOH)和枯烯氫過氧化物(COOH)等發(fā)生反應,從而清除體內多余的活性氧,在植物抵御生物及非生物脅迫過程中發(fā)揮重要的作用[3]。例如,擬南芥AtGPX3基因的缺失突變體表現(xiàn)出活性氧清除能力和耐旱能力降低等表型[4];敲除AtGPX8基因后,在paraquat(一種除草劑)的處理下,擬南芥幼苗根部生長慢于野生型,而過表達植株根部生長快于野生型[5]。

    以往對植物GPX的研究主要集中在擬南芥[4-5]、水稻[6-7]、楊樹[8]、胡楊[9]、番茄[10]、甘藍型油菜[11]、百脈根[12]、鹽芥[13]、香蕉[14]等被子植物,而對植物中另一大類群-苔蘚植物的研究相對較少。小立碗蘚是目前完成全基因組測序的唯一苔蘚類植物,在分類上屬于葫蘆蘚科小立碗蘚屬,是典型的早期陸地苔蘚植物[15],作為植物分子生物學研究極具前景的模式植物已日益受到人們的重視。本研究以小立碗蘚為研究材料,從小立碗蘚基因組中挖掘到3個GPX基因,并對其序列特征、基因結構、表達模式、蛋白亞細胞定位、蛋白結構、蛋白生化性質及關鍵氨基酸位點進行了詳細的比較研究,揭示了小立碗蘚GPX基因家族的功能分化和分子特征。

    1材料和方法

    1.1實驗材料

    實驗材料小立碗蘚(Physcomitrellapatens)由北京大學郭紅衛(wèi)教授惠贈。

    1.2小立碗蘚GPX基因的鑒定

    以擬南芥(Arabidopsisthaliana)GPX蛋白序列(GenBank登錄號NP_198297.2)為模板,在Phytozome網(wǎng)站(http://www.phytozome.org)中,應用Blast工具檢索小立碗蘚全基因組數(shù)據(jù)庫,搜索得到小立碗蘚GPX基因,將其所編碼的蛋白在NCBI網(wǎng)站進行保守結構域分析(http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi),以確定其含有GPX結構域。

    1.3小立碗蘚GPX的系統(tǒng)發(fā)生分析

    將搜索得到的小立碗蘚GPX蛋白序列與其他陸地植物GPX蛋白序列通過Muscle(http://www.ebi.ac.uk/Tools/msa/muscle)進行序列比對,并進行手工校對,將比對結果導入MEGA5.05軟件,使用Neighbour-joining模型構建系統(tǒng)進化樹,Bootstrap設為1 000。

    1.4小立碗蘚GPX基因的克隆

    在BCDATG固體培養(yǎng)基培養(yǎng)小立碗蘚[16],培養(yǎng)條件為18 h光照、6 h暗培養(yǎng),培養(yǎng)溫度25 ℃。培養(yǎng)60 d后,用Aurum Total RNA Kits(Bio-Rad Laboratories)提取植株總RNA。將提取的RNA用無RNase的DNase I(Promega)處理后,用RNA PCR Kit(AMV) version 3.0 (TaKaRa)反轉錄成cDNA。根據(jù)小立碗蘚PpGPX1、PpGPX2和PpGPX3基因序列分別設計引物對PpGPX1-EX1/2、PpGPX2-EX1/2和PpGPX3-EX1/2(表1),引物均由Invitrogen公司合成,以小立碗蘚cDNA為模板分別進行PCR擴增。PCR反應程序為94 ℃預變性3 min;94 ℃變性30 s,55 ℃退火40 s,72 ℃延伸1 min,35個循環(huán);最后72 ℃延伸3 min。PCR產物經(jīng)1%瓊脂糖凝膠電泳檢測并用DNA純化回收試劑盒(北京鼎國昌盛生物技術有限公司)回收目的條帶,將純化產物連接入pEASY-T3載體(北京全式金生物),將連接產物轉化大腸桿菌Trans-T1感受態(tài)細胞,在含有X-gal/IPTG/Amp的LB平板上挑取多個陽性克隆進行雙向測序。最終得到測序正確的PpGPX/pEASY-T3重組質粒。

    1.5小立碗蘚GPX基因的表達模式分析

    選擇培養(yǎng)60 d后生長狀況良好的小立碗蘚,分別用0.5% H2O2、150 mmol/L NaCl、1 mmol/L水楊酸、0.1%莠去津和1 mmol/L CDNB噴灑。H2O2和莠去津處理12 h,其余方式每種處理24 h,每種處理設置3個重復。處理后的植株提取總RNA,用無RNase的DNaseI處理后反轉錄成cDNA。根據(jù)小立碗蘚基因序列設計3對引物(PpGPX1-SP1/2、PpGPX2-SP1/2和PpGPX3-SP1/2,表1),進行PCR擴增,以小立碗蘚Actin基因為內標(GenBank登錄號XM_001783849)。PCR反應程序為94 ℃預變性3 min;94 ℃變性30 s,60 ℃退火40 s,72 ℃延伸1 min,循環(huán)數(shù)分別設置24、26、28和30;最后72 ℃延伸3 min。PCR產物經(jīng)1%瓊脂糖凝膠電泳檢測并測序確定。

    1.6小立碗蘚GPX蛋白的亞細胞定位分析

    以小立碗蘚cDNA為模板,利用引物PpGPX1-SL1/2和PpGPX2-SL1/2分別進行PCR擴增。PCR產物經(jīng)1%瓊脂糖凝膠電泳檢測,并用DNA純化回收試劑盒(北京鼎國昌盛生物技術有限公司)回收目的條帶,將純化產物連接入pEASY-T3載體(北京全式金生物),將構建好的重組質粒轉化大腸桿菌Trans-T1感受態(tài)細胞,挑選多個陽性克隆進行測序。將測序正確的重組質粒用限制性內切酶BamHⅠ和MluⅠ(NEB)酶切,將酶切后的DNA片段分別連接到pPSN轉基因載體[17],并將轉基因重組質粒轉入事先制備好的擬南芥原生質體。原生質體的制備及轉化參考Sang-Dong Yoo等[18]的方法。轉化1 d后,用激光共聚焦顯微鏡(Olympus FV1000 MPE)進行觀察。GFP激發(fā)光波長為488 nm,葉綠體激發(fā)光波長為543 nm。

    1.7小立碗蘚PpGPX1和PpGPX2蛋白結構模擬

    目前,在植物中只有楊樹PtGPX5蛋白的X射線衍射晶體結構得到解析,因此我們以PtGPX5蛋白的晶體結構為模板,使用InsightII軟件中的Homology模塊對PpGPX1和PpGPX2的蛋白結構進行模擬。首先用Align 2D程序進行序列比對,然后使用Modeler程序進行蛋白質三維結構的模擬,利用profile-3D程序對所有模擬的結構進行校驗,選擇校驗值最高的蛋白結構。

    1.8小立碗蘚PpGPX1和PpGPX2蛋白的表達和純化

    將PpGPX1/pEASY-T3和PpGPX2/pEASY-T3重組質粒用限制性內切酶BamHⅠ和NotⅠ(NEB)酶切,目的片段經(jīng)DNA純化回收試劑盒回收后連接到改造過的pET30a蛋白表達載體中[19],使表達蛋白在N端含有6×His標簽。將構建好的重組質粒轉化大腸桿菌Tuner(DE3)感受態(tài)細胞,測序確定菌中質粒序列的正確性。在定點突變實驗中,PpGPX2的Pro159、Phe167和Phe172殘基分別用Ala替換,通過DNA測序確定所需要的位點進行了正確的突變,然后將突變后的PpGPX2基因克隆到改造后的pET30a蛋白表達載體中并轉化Tuner(DE3)感受態(tài)細胞。挑取陽性克隆菌在LB液體培養(yǎng)基中過夜培養(yǎng),將菌體1:100稀釋后,繼續(xù)培養(yǎng)至OD600≈0.5。加入終濃度為0.05 mmol/L IPTG誘導GPX蛋白表達。37 ℃誘導過夜后,離心收集菌體(1 000 r/min,3 min,4 ℃)。用緩沖液A(20 mmol/L磷酸鈉、0.5 mol/L NaCl和20 mmol/L咪唑,pH 7.4)重懸菌體,冰上超聲破碎裂解。將破碎液(10 000 r/min,10 min,4 ℃)離心,離心后上清和沉淀各取50 μL進行SDS-PAGE檢測(SDS-PAGE檢測采用12%分離膠和5%濃縮膠),剩余上清置于緩沖液A預平衡后的鎳離子親和層析柱(GE Healthcare Bio-Sciences)進行純化。雜蛋白用緩沖液A洗滌,目的蛋白用緩沖液B(20 mmol/L磷酸鈉、0.5 mol/L NaCl和500 mmol/L咪唑,pH 7.4)洗脫。

    表1 本研究中使用的引物

    1.9小立碗蘚PpGPX1和PpGPX2蛋白酶學性質測定及動力學分析

    小立碗蘚GPX蛋白的生化性質測定參照Zhao等[3]的方法,通過檢測340 nm波長處NADPH吸光度(A340)的降低來測定過氧化物酶活性,所有的反應都在25 ℃進行。反應體系包括:150 mmol/L Tris-HCl,pH 8.0,7.5 mmol/L EDTA,1 mmol/L H2O2,不同濃度的GPX蛋白及電子供體系統(tǒng)。電子供體系統(tǒng)包括Trx電子供體系統(tǒng)(0.2 mmol/L NADPH,9 μmol/LE.coliTrx,0.5 μmol/LE.coliTrx reductase)和GSH電子供體系統(tǒng)(0.2 mmol/L NADPH,2 mmol/L GSH和1.5 U谷胱甘肽過氧化物酶)。當測定GPX蛋白對tBOOH或COOH的活性時,將反應體系中的H2O2替換為叔丁基過氧化氫(tBOOH)或異丙基苯過氧化氫(COOH)。對照反應用等量的緩沖液B來替代GPX蛋白溶液,其他成分不變,每個反應設3次獨立重復實驗。

    動力學分析通過測量GPX蛋白在不同濃度的Trx和H2O2條件下的催化活性進行分析。在分析GPX蛋白對Trx的動力學時,固定H2O2的濃度為1.0 mmol/L,通過測量Trx濃度范圍為0.5~25 μmol/L時,GPX蛋白在不同Trx濃度下的催化活性獲得。在分析GPX蛋白對H2O2的動力學時,固定Trx濃度為5.0 μmol/L,通過測量Trx濃度范圍為0.005~1.0 mmol/L時,GPX蛋白在不同Trx濃度下的催化活性獲得。檢測的結果通過Hyper32程序(http://homepage.ntlworld.com/john.easterby/hyper32.html.)進行非線性回歸分析。

    2結果與分析

    2.1小立碗蘚GPX基因家族的挖掘和系統(tǒng)發(fā)育分析

    從小立碗蘚基因組中挖掘到3個可能的GPX基因,分別命名為PpGPX1、PpGPX2和PpGPX3,通過CDD保守結構域分析發(fā)現(xiàn),3個基因編碼的蛋白均具有GPX蛋白結構域,證明得到的基因為GPX基因。PpGPX1基因編碼170個氨基酸的蛋白,預測分子量為19.1 kD;PpGPX2基因編碼247個氨基酸的蛋白,預測分子量為26.8 kD;PpGPX3基因編碼184個氨基酸的蛋白,預測分子量為20.8 kD。將3個基因的cDNA序列和基因組序列進行比較以鑒定其基因結構,發(fā)現(xiàn)PpGPX1和PpGPX3基因都只含有1個外顯子,而PpGPX2基因有6個外顯子(圖1)。

    將小立碗蘚和蕨類植物江南卷柏(Selaginellamoellendorffii)、裸子植物油松(Pinustabuliformis)、單子葉植物水稻(Oryzasativa)、雙子葉植物毛果楊(Populustrichocarpa)和擬南芥(Arabidopsisthaliana)的GPX蛋白序列進行序列相似性分析(圖2),結果發(fā)現(xiàn),這些GPX間具有比較高的序列相似性,蛋白序列相似性范圍為39.6%~99.3%。將小立碗蘚與其他陸地植物的GPX進行系統(tǒng)發(fā)生關系分析,發(fā)現(xiàn)小立碗蘚GPX在最基部聚成一枝(圖4),說明小立碗蘚GPX基因家族在陸地植物分化的早期就產生了。

    2.2小立碗蘚GPX基因表達模式分析

    通過半定量RT-PCR檢測小立碗蘚GPX基因在正常生長以及H2O2、NaCl、水楊酸、莠去津和CDNB等脅迫下的表達模式,發(fā)現(xiàn)PpGPX1和PpGPX2在24~30個PCR擴增循環(huán)時,在不同生長條件下均表達,而PpGPX3在檢測的所有樣品中均不表達(圖3)。為了進一步確認PpGPX3基因的表達模式,用PpGPX3基因序列檢索NCBI中的小立碗蘚EST數(shù)據(jù)庫(含有382 587條小立碗蘚EST序列),未搜到任何匹配的EST序列,因此推測PpGPX3基因可能不表達,或可能在特定的發(fā)育時期或誘導條件下表達。

    圖1 小立碗蘚GPX基因的基因結構

    圖2 植物GPX蛋白序列比對

    圖3 小立碗蘚GPX基因在不同脅迫條件下的表達

    2.3小立碗蘚PpGPX1和PpGPX2蛋白的亞細胞定位分析

    為了研究小立碗蘚PpGPX1和PpGPX2蛋白的亞細胞定位,構建了PpGPX1和PpGPX2的C末端與綠色熒光蛋白eGFP相連的融合蛋白表達質粒,在擬南芥原生質體中瞬時表達后,用激光共聚焦顯微鏡進行觀察。結果發(fā)現(xiàn),PpGPX2-eGFP融合蛋白的綠色熒光與葉綠體自發(fā)熒光重疊,表明PpGPX2蛋白能夠定位到葉綠體。而PpGPX1-eGFP融合蛋白則與單獨表達eGFP時的熒光定位相似,說明PpGPX1定位到細胞質(圖5)。

    圖4 植物GPX基因的系統(tǒng)發(fā)育分析

    圖6 小立碗蘚PpGPX1和PpGPX2蛋白三維結構的模擬

    圖5 小立碗蘚GPX蛋白在擬南芥原生質體中的亞細胞定位

    2.4小立碗蘚PpGPX1和PpGPX2蛋白結構模擬

    目前,在植物中只有楊樹PtGPX5蛋白的X射線衍射晶體結構得到解析,因此我們以PtGPX5蛋白的晶體結構為模板,使用InsightⅡ軟件對小立碗蘚PpGPX1和PpGPX2的蛋白結構進行模擬(圖6,A),以profile-3D程序對模擬得到的結構進行校驗,校驗后得到的2個蛋白所有位點的校驗值均為正值,表明模擬所得到的結果具有較高的可信度(圖6,B)。PpGPX1和PpGPX2都含有4個α螺旋和6個β折疊,其中4個β折疊構成反向平行的β片層結構,將2個蛋白結構進行疊加,發(fā)現(xiàn)2個蛋白結構十分保守,除了loop區(qū)域有一些差異外,其余部分幾乎完全重疊。

    2.5小立碗蘚PpGPX1和PpGPX2蛋白的表達純化及酶學性質分析

    PpGPX2定位到葉綠體,其在N端含有一段信號肽,在植物細胞中加工成為成熟蛋白時,信號肽會被切除。為研究PpGPX1和PpGPX2蛋白的性質,將PpGPX1和去掉信號肽的PpGPX2在大腸桿菌中進行原核表達。在37 ℃誘導條件下,2個蛋白均能以可溶性蛋白的形式表達。經(jīng)過鎳離子親和層析純化后得到純度較高的蛋白,SDS-PAGE檢測結果顯示,純化蛋白的分子量大小與預測相一致(圖7)。

    先前的研究發(fā)現(xiàn),植物GPX能夠利用Trx和GSH還原系統(tǒng)還原過氧化物[8,20]。本研究分別以Trx和GSH作為還原劑檢測PpGPX1和PpGPX2蛋白對H2O2、tBOOH和COOH的催化活性。當用GSH電子供體系統(tǒng)時,PpGPX1和PpGPX2對3種底物均沒有催化活性,而應用Trx電子供體系統(tǒng)時,2個蛋白對3種底物均有催化活性(表2)。同一蛋白對3種底物的活性沒有明顯的差別,而2個蛋白之間的催化活性不同,PpGPX2蛋白對3種過氧化物底物的活性均高于PpGPX1蛋白,說明小立碗蘚GPX家族不同成員的活性產生了分化。

    分別以Trx和H2O2作為電子供體和過氧化物底物對PpGPX1和PpGPX2進行酶促動力學分析(表3)。結果表明,小立碗蘚GPX蛋白對Trx的Km值大于對H2O2的Km值,說明小立碗蘚GPX蛋白對H2O2的親和力高于對Trx的親和力(親和力越大Km值越小)。PpGPX2對Trx的親和力是PpGPX1的5.5倍,而對H2O2的親和力PpGPX1高于PpGPX2,PpGPX2對Trx和H2O2的催化效率均高于PpGPX1。

    2.6小立碗蘚GPX蛋白關鍵氨基酸位點功能分析

    根據(jù)先前的分子對接實驗[21],楊樹PtGPX5的Pro81、Phe90和Phe95氨基酸殘基可能在PtGPX5與Trx的互作中發(fā)揮作用。通過序列比對發(fā)現(xiàn),這3個氨基酸在小立碗蘚中分別為Pro159、Phe167和Phe172氨基酸殘基(圖2)。由于PpGPX2蛋白的催化活性和對Trx的親和力都高于PpGPX1蛋白,所以選擇PpGPX2蛋白研究這3個位點在小立碗蘚GPX中的作用。根據(jù)序列比對結果,在PpGPX2中這3個位點分別對應Pro158、Phe167和Phe172。將這3個氨基酸殘基分別突變?yōu)锳la后發(fā)現(xiàn),與野生型PpGPX2蛋白相比,3個突變體蛋白對3種過氧化物底物的活性均發(fā)生了不同程度的降低,其中PpGPX2(F172A)突變體蛋白對3種底物的催化活性下降最多,對H2O2、tBOOH和COOH的活性分別下降了84.8%、87.3%和88.3%。動力學分析得出PpGPX2(F172A)突變體蛋白對Trx和H2O2的kcat/Km值降低至野生型的15.80%和4.69%。這些結果說明Pro158、Phe167和Phe172氨基酸位點對于PpGPX2的催化具有重要作用。

    圖7 純化的小立碗蘚GPX蛋白SDS-PAGE電泳

    蛋白名稱ProteinnameKmTrx/(mmol·L-1)kcatTrx/(S-1)(kcat/Km)Trx/(L·mmol-1·S-1)KmH2O2/(mmol·L-1)kcatH2O2/(S-1)(kcat/Km)H2O2/(L·mmol-1·S-1)PpGPX112.140±0.5801.3310.1100.003±0.0010.13644.445PpGPX22.189±0.1750.7100.3240.005±0.0010.24654.348PpGPX2(F167A)0.411±0.1680.1160.3010.002±0.0000.10555.095PpGPX2(F172A)0.668±0.3070.0340.0510.007±0.0010.0152.548PpGPX2(P158A)2.911±0.7700.4030.1390.010±0.0010.10710.637

    表2 小立碗蘚GPX蛋白對不同底物的催化活性

    3討論

    GPX作為一種清除活性氧的關鍵酶,在植物抵御外界生物及非生物脅迫的過程中發(fā)揮著重要的作用。以往對GPX的研究主要集中于擬南芥、水稻、楊樹、胡楊、番茄、甘藍型油菜、百脈根、鹽芥、香蕉等被子植物,而對植物中另一大類群-苔蘚植物的研究相對較少,本實驗對苔蘚植物小立碗蘚GPX家族進行了深入的功能分化與分子特性研究,對于了解小立碗蘚抗氧化脅迫的機制具有重要的意義。

    雖然小立碗蘚GPX家族成員間有較高的序列相似性,但它們在基因結構、表達模式、蛋白亞細胞定位和蛋白酶學功能方面均有差異。例如,PpGPX1和PpGPX3都只有1個外顯子,而PpGPX2有6個外顯子;PpGPX1和PpGPX2在檢測的所有生長條件下均表達,而PpGPX3在所有檢測的生長條件下均不表達;PpGPX1蛋白定位到細胞質,而PpGPX2蛋白定位到葉綠體;PpGPX2蛋白對H2O2、tBOOH和COOH的催化活性及催化效率均高于PpGPX1蛋白。基因結構、表達模式、蛋白亞細胞定位和蛋白酶學功能方面的差異表明,小立碗蘚GPX家族成員間已經(jīng)發(fā)生了廣泛的功能分化。

    根據(jù)先前的研究,在動物和真菌中,GPX蛋白可以利用Trx電子供體系統(tǒng)和GSH電子供體系統(tǒng)進行過氧化物催化反應[22-23];水稻、油松、楊樹等高等植物則只能利用Trx電子供體系統(tǒng),不能使用GSH電子供體系統(tǒng)[20,24-25]。本研究利用Trx電子供體系統(tǒng)和GSH電子供體系統(tǒng)對小立碗蘚PpGPX1和PpGPX2蛋白進行酶學性質檢測,結果發(fā)現(xiàn),小立碗蘚GPX蛋白也只能利用Trx電子供體系統(tǒng)在催化反應中進行電子傳遞,而不能利用GSH電子供體系統(tǒng),說明陸地植物GPX是Trx依賴的過氧化物酶。

    先前的研究表明,植物GPX蛋白通過兩步反應完成一次催化循環(huán),首先GPX蛋白與過氧化物底物反應,將過氧化物還原,自身變?yōu)檠趸瘧B(tài),然后通過Trx電子供體系統(tǒng)將GPX蛋白還原,在這一反應過程中GPX蛋白和Trx蛋白會發(fā)生直接的相互作用[26]。根據(jù)先前的分子對接實驗[21],楊樹PtGPX5的Pro81、Phe90和Phe95氨基酸殘基可能在PtGPX5與Trx的互作中發(fā)揮作用。將楊樹PtGPX5蛋白第90位的Phe突變?yōu)镚lu后,PtGPX5突變體蛋白對Trx的親和力明顯下降[27]。根據(jù)序列比對分析,楊樹PtGPX5的Pro81、Phe90和Phe95氨基酸位點分別對應PpGPX2的Pro158、Phe167和Phe172氨基酸位點。通過蛋白的定點突變分析發(fā)現(xiàn),將PpGPX2的Pro158、Phe167和Phe172氨基酸殘基均突變?yōu)锳la后,3個突變體的催化活性均明顯降低,其中PpGPX2(F172A)突變體蛋白對H2O2、tBOOH和COOH的活性分別降低為原來的15.2%、12.7%和11.7%,說明這3個位點,尤其是第172位的Phe對小立碗蘚GPX蛋白的催化功能具有重要作用。

    本研究通過整合序列分析、基因表達模式分析、蛋白亞細胞定位分析、蛋白質結構模擬、蛋白酶學性質分析及關鍵氨基酸位點的研究,對苔蘚植物小立碗蘚GPX家族的功能分化及分子特性進行了深入的研究,為了解小立碗蘚抵抗氧化脅迫的機制提供了一定的信息。

    參考文獻:

    [1]MITTLER R,VANDERAUWERA S,GOLLERY M,etal.Reactive oxygen gene network of plants[J].TrendsinPlantScience,2004,9(10):490-498.

    [2]APEL K,HIRT H.Reactive oxygen species:Metabolism,oxidative stress,and signal transduction[J].AnnualReviewofPlantBiology,2004,55(1):373-399.

    [3]ZHAO L,HAN X M,WANG W,etal.Molecular and catalytic properties of glutathione peroxidase family proteins fromPinustabulaeformis[J].PlantMolecularBiologyReporter,2014,32(3):771-778.

    [4]MIAO Y,LV D,WANG P C,etal.AnArabidopsisglutathione peroxidase functions as both a redox transducer and a scavenger in abscisic acid and drought stress responses[J].PlantCell,2006,18(10):2 749-2 766.

    [5]GABER A,OGATA T,MARUTA T,etal.The involvement ofArabidopsisglutathione peroxidase 8 in the suppression of oxidative damage in the nucleus and cytosol[J].PlantandCellPhysiology,2012,53(9):1 596-1 606.

    [6]LI W J,FENG H,FAN J H,etal.Molecular cloning and expression of a phospholipid hydroperoxide glutathione peroxidase homolog inOryzasativa[J].BiochimicaEtBiophysicaActa-GeneStructureandExpression,2000,1 493(1-2):225-230.

    [7]WEI J CH(韋金池),YANG H L(楊海靈),Functional divergence of two glutathione peroxidase genes inOryzasativa[J].PlantScienceJournal(植物科學學報),2013,31(1):64-72(in Chinese).

    [8]NAVROT N,COLLIN V,GUALBERTO J,etal.Plant glutathione peroxidases are functional peroxiredoxins distributed in several subcellular compartments and regulated during biotic and abiotic stresses[J].PlantPhysiology,2006,142(4):1 364-1 379.

    [9]WANG F F(王菲菲),DING M Q(丁明全),etal.Cloning of glutathione peroxidase genePeGPXfromPopuluseuphraticaand the salt tolerance of the transformed plants[J].GenomicsandAppliedBiology(基因組學與應用生物學),2012,31(3):231-239(in Chinese).

    [10]DEPEGE N,DREVET J,BOYER N.Molecular cloning and characterization of tomato cDNAs encoding glutathione peroxidase-like proteins[J].EuropeanJournalofBiochemistry,1998,253(2):445-451.

    [11]HU M L(胡茂龍),LONG W H(龍衛(wèi)華),GAO J Q(高建芹),etal.Glutathione peroxidase geneBnGPX1 cloning from rapeseed and expression under abiotic stress[J].ChineseJournalofOilCropSciences(中國油料作物學報),2011,33(4):331-337(in Chinese).

    [12]RAMOS J,MATAMOROS M A,NAYA L,etal.The glutathione peroxidase gene family ofLotusjaponicus:characterization of genomic clones,expression analyses and immunolocalization in legumes[J].NewPhytologist,2009,181(1):103-114.

    [13]MA T T(馬亭亭),ZHOU Y J(周宜君),GAO F(高飛),etal.Cloning and expression analysis of glutathione peroxidase 6 (ThGPX6) inThellungiellahalophile[J].JournalofPlantGeneticResources(植物遺傳資源學報),2012,13(2):252-258(in Chinese).

    [14]ZHANG L L(張麗麗),XU B Y(徐碧玉),LIU J H(劉菊華),etal.Isolation and expression analysis of a cDNA encoding glutathione peroxidase from banana[J].ActaHorticulturaeSinica(園藝學報),2012,39(8):1 471-1 481(in Chinese).

    [15]ZHAO H(趙奐),ZHA X G(趙曉剛),HE Y K(何奕昆),etal.Physcomitrellapatens,a potential model system in plant molecular biology[J].ChineseBulletinofBotany(植物學通報),2004,21(2):129-138(in Chinese).

    [16]NISHIYAMA T,HIWATASHI Y,SAKAKIBARA K,etal.Tagged mutagenesis and gene-trap in the moss,Physcomitrellapatensby shuttle mutagenesis[J].DNAResearch,2000,7(1):9-17.

    [17]TANG Z X,YANG H L.Functional divergence and catalytic properties of dehydroascorbate reductase family proteins fromPopulustomentosa[J].MolecularBiologyReports,2013,40(8):5 105-5 114.

    [18]YOO S D,CHO Y H,SHEEN J.Arabidopsismesophyll protoplasts:a versatile cell system for transient gene expression analysis[J].NatureProtocols,2007,2(7):1 565-1 572.

    [19]YANG X,SUN W,LIU J P,etal.Biochemical and physiological characterization of a tau class glutathione transferase from rice (Oryzasativa)[J].PlantPhysiol.Biochem.,2009,47(11-12):1 061-1 068.

    [20]JUNG B G,LEE K O,LEE S S,etal.A Chinese cabbage cDNA with high sequence identity to phospholipid hydroperoxide glutathione peroxidases encodes a novel isoform of thioredoxin-dependent peroxidase[J].JournalofBiologicalChemistry,2002,277(15):12 572-12 578.

    [21]KOH C S,DIDIERJEAN C,NAVROT N,etal.Crystal structures of a poplar thioredoxin peroxidase that exhibits the structure of glutathione peroxidases:Insights into redox-driven conformational changes[J].JournalofMolecularBiology,2007,370(3):512-529.

    [22]CHANG C,SLESAK I,JORDA L,etal.Arabidopsischloroplastic glutathione peroxidases play a role in cross talk between photooxidative stress and immune responses[J].PlantPhysiology,2009,150(2):670-683.

    [23]FU L H,WANG X F,EYAL Y,etal.A selenoprotein in the plant kingdom.Mass spectrometry confirms that an opal codon (UGA) encodes selenocysteine inChlamydomonasreinhardtiiglutathione peroxidase[J].JournalofBiologicalChemistry,2002,277(29):25 983-25 991.

    [24]IQBAL A,YABUTA Y,TAKEDA T,etal.Hydroperoxide reduction by thioredoxin-specific glutathione peroxidase isoenzymes ofArabidopsisthaliana[J].FEBS,2006,273(24):5 589-5 597.

    [25]WANG Z B,WANG F,DUAN R,etal.Purification and physicochemical characterization of a recombinant phospholipid hydroperoxide glutathione peroxidase fromOryzasativa[J].JournalofBiochemistryandMolecularBiology,2007,40(3):412-418.

    [26]DELAUNAY A,PFLIEGER D,BARRAULT M B,etal.A thiol peroxidase is an H2O2receptor and redox-transducer in gene activation[J].Cell,2002,111(4):471-481.

    [27]SELLES B,HUGO M,TRUJILLO M,etal.Hydroperoxide and peroxynitrite reductase activity of poplar thioredoxin-dependent glutathione peroxidase 5:kinetics,catalytic mechanism and oxidative inactivation[J].BiochemicalJournal,2012,442(2):369-380.

    (編輯:宋亞珍)

    Functional Divergence and Catalytic Properties of Glutathione

    Peroxidase Family fromPhyscomitrellapatens

    WANG Wei,YANG Zhiling,YANG Hailing*

    (College of Life Sciences and Biotechnology,Beijing Forestry University,Beijing 100083,China)

    Abstract:In plants,glutathione peroxidase (GPX) plays important roles in protecting cells against oxidative stress.In this study,threeGPXs(PpGPX1/2/3) were identified from thePhyscomitrellapatensgenome.PpGPX1 andPpGPX3 only had one exon,whilePpGPX2 had six exons.Semiquantitative RT-PCR analysis revealed thatPpGPX1 andPpGPX2 were expressed in all tested conditions.However,PpGPX3 transcripts were not detected in any tested conditions.By transiently expressing C-terminal eGFP-GPX fusions inArabidopsisprotoplasts,we found that PpGPX1 and PpGPX2 were localized to cytoplasm and chloroplast,respectively.PpGPX1 and PpGPX2 were then over-expressed inE.coliand purified using a Ni Sepharose High Performance column.Enzymatic analysis found thatP.patensGPXs showed enzymatic activity towards peroxide substrates using Trx as electron donor,but not GSH.PpGPX2 showed higher catalytic activity and catalytic efficiency towards peroxide substrates than that of PpGPX1.The differences of gene structures,gene expression,protein subcellular localization and enzymatic characteristics amongP.patensGPXs indicating functional divergence.The function of Pro158,Phe167 and Phe172 of PpGPX2 were examined by site-directed mutagenesis.The mutant proteins showed decreased catalytic activity which indicated that these three residues were important for enzymatic activity.

    Key words:glutathione peroxidase;Physcomitrellapatens;functional divergence;subcellular localization;enzymatic characterization

    中圖分類號:Q786;Q789

    文獻標志碼:A

    作者簡介:王煒(1989-),男,在讀碩士研究生,主要從事蛋白質結構與生化性質研究。E-mail:wangw1989@163.com*通信作者:楊海靈,博士,教授,主要從事植物蛋白質結構與功能研究。E-mail:yhailing77@163.com

    基金項目:國家自然科學基金(31270641)

    收稿日期:2014-09-02;修改稿收到日期:2014-12-22

    文章編號:1000-4025(2015)01-0001-09

    doi:10.7606/j.issn.1000-4025.2015.01.0001

    亚洲精品国产精品久久久不卡| 黄色片一级片一级黄色片| 国产日韩欧美在线精品| 国产免费av片在线观看野外av| 久热这里只有精品99| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 飞空精品影院首页| 在线永久观看黄色视频| 久久免费观看电影| 人妻久久中文字幕网| 日韩免费高清中文字幕av| 91字幕亚洲| 亚洲精品美女久久久久99蜜臀| 老鸭窝网址在线观看| 国产精品亚洲av一区麻豆| 啦啦啦在线免费观看视频4| 韩国高清视频一区二区三区| 亚洲第一青青草原| 最近最新免费中文字幕在线| 中文字幕最新亚洲高清| 国产亚洲一区二区精品| 中文字幕人妻丝袜一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 青青草视频在线视频观看| 侵犯人妻中文字幕一二三四区| 亚洲精品一区蜜桃| av视频免费观看在线观看| 悠悠久久av| 亚洲人成电影观看| 高清黄色对白视频在线免费看| 久久天堂一区二区三区四区| 黄色视频,在线免费观看| 天天操日日干夜夜撸| 亚洲五月婷婷丁香| 一本色道久久久久久精品综合| 久久人人爽av亚洲精品天堂| 50天的宝宝边吃奶边哭怎么回事| 一区二区三区四区激情视频| 大片免费播放器 马上看| av电影中文网址| 狂野欧美激情性bbbbbb| 国产亚洲av高清不卡| 成人av一区二区三区在线看 | av视频免费观看在线观看| 国产一区二区三区综合在线观看| 日韩视频一区二区在线观看| 交换朋友夫妻互换小说| 美女主播在线视频| 色综合欧美亚洲国产小说| 国产av国产精品国产| 一二三四在线观看免费中文在| 国产淫语在线视频| 欧美亚洲 丝袜 人妻 在线| 丰满迷人的少妇在线观看| 麻豆乱淫一区二区| 日韩欧美一区二区三区在线观看 | 久久亚洲精品不卡| 嫁个100分男人电影在线观看| 正在播放国产对白刺激| 久久久久久人人人人人| 国产精品av久久久久免费| 亚洲国产av影院在线观看| 亚洲中文日韩欧美视频| 91字幕亚洲| 亚洲人成77777在线视频| 国产麻豆69| 亚洲成人手机| 国产精品久久久av美女十八| 一级,二级,三级黄色视频| av一本久久久久| 精品少妇一区二区三区视频日本电影| 国产伦理片在线播放av一区| 成人黄色视频免费在线看| 日本av手机在线免费观看| 不卡一级毛片| 精品一品国产午夜福利视频| 日韩三级视频一区二区三区| 韩国高清视频一区二区三区| 一个人免费看片子| 欧美日韩中文字幕国产精品一区二区三区 | 大香蕉久久成人网| 人人妻人人爽人人添夜夜欢视频| 亚洲九九香蕉| 99精品久久久久人妻精品| 一进一出抽搐动态| 一本—道久久a久久精品蜜桃钙片| 少妇人妻久久综合中文| 成人国产av品久久久| 丝袜美腿诱惑在线| 汤姆久久久久久久影院中文字幕| 桃花免费在线播放| 午夜福利视频精品| 一本大道久久a久久精品| 亚洲精品一卡2卡三卡4卡5卡 | 久久国产精品大桥未久av| 欧美亚洲 丝袜 人妻 在线| 国产区一区二久久| 9色porny在线观看| 黄色视频在线播放观看不卡| 欧美日韩亚洲高清精品| 丁香六月天网| 99国产精品免费福利视频| 亚洲欧美精品自产自拍| 成年人黄色毛片网站| 中文精品一卡2卡3卡4更新| 在线观看免费高清a一片| 久久人人爽人人片av| 999精品在线视频| 成年动漫av网址| 大陆偷拍与自拍| 欧美成狂野欧美在线观看| 精品欧美一区二区三区在线| 人人澡人人妻人| 如日韩欧美国产精品一区二区三区| 人人澡人人妻人| 精品少妇内射三级| 妹子高潮喷水视频| 国产一区有黄有色的免费视频| 日韩精品免费视频一区二区三区| 蜜桃在线观看..| 香蕉丝袜av| 中文字幕人妻熟女乱码| 欧美少妇被猛烈插入视频| 欧美性长视频在线观看| 成年av动漫网址| 宅男免费午夜| 国产成人精品无人区| 国产精品秋霞免费鲁丝片| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品第二区| h视频一区二区三区| 精品免费久久久久久久清纯 | av在线app专区| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区91| 亚洲免费av在线视频| 窝窝影院91人妻| 国产精品一二三区在线看| 日韩精品免费视频一区二区三区| 99精品久久久久人妻精品| 女人爽到高潮嗷嗷叫在线视频| 欧美在线黄色| 男人操女人黄网站| 热99国产精品久久久久久7| 两个人免费观看高清视频| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 精品国产一区二区三区四区第35| 欧美人与性动交α欧美精品济南到| 亚洲av片天天在线观看| 考比视频在线观看| 少妇人妻久久综合中文| 日韩 亚洲 欧美在线| www.999成人在线观看| av片东京热男人的天堂| 日韩一区二区三区影片| 最近中文字幕2019免费版| √禁漫天堂资源中文www| 亚洲精品在线美女| 免费观看av网站的网址| 国产成人av激情在线播放| 日韩制服丝袜自拍偷拍| 1024视频免费在线观看| 热re99久久国产66热| 在线观看一区二区三区激情| 成人黄色视频免费在线看| 美国免费a级毛片| 人妻 亚洲 视频| 久久香蕉激情| 777久久人妻少妇嫩草av网站| 人人妻人人爽人人添夜夜欢视频| 色综合欧美亚洲国产小说| 性色av一级| 国产精品一区二区在线观看99| 婷婷丁香在线五月| 伊人久久大香线蕉亚洲五| 91av网站免费观看| av在线播放精品| 国产在线观看jvid| 美女主播在线视频| 亚洲精品第二区| 色94色欧美一区二区| 高清av免费在线| 成人av一区二区三区在线看 | 在线观看www视频免费| 女性生殖器流出的白浆| 国产高清国产精品国产三级| 91字幕亚洲| 啪啪无遮挡十八禁网站| 王馨瑶露胸无遮挡在线观看| 黄片小视频在线播放| 操美女的视频在线观看| 亚洲精品国产区一区二| 精品福利永久在线观看| 丰满少妇做爰视频| 在线天堂中文资源库| 自拍欧美九色日韩亚洲蝌蚪91| 久久国产精品影院| 亚洲七黄色美女视频| 久久综合国产亚洲精品| 男女无遮挡免费网站观看| 国产欧美日韩一区二区精品| 少妇 在线观看| 午夜福利在线免费观看网站| 国产精品一区二区免费欧美 | 精品一区在线观看国产| 国产在线观看jvid| 亚洲五月婷婷丁香| 黄色a级毛片大全视频| 免费一级毛片在线播放高清视频 | 极品少妇高潮喷水抽搐| 久久国产精品男人的天堂亚洲| 午夜福利免费观看在线| 国产成人一区二区三区免费视频网站| 精品少妇内射三级| 日本五十路高清| 免费高清在线观看日韩| 午夜精品久久久久久毛片777| 啦啦啦在线免费观看视频4| 777久久人妻少妇嫩草av网站| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品成人久久小说| 少妇被粗大的猛进出69影院| 在线观看免费日韩欧美大片| 看免费av毛片| 欧美 日韩 精品 国产| 久久国产精品大桥未久av| 自拍欧美九色日韩亚洲蝌蚪91| 欧美精品亚洲一区二区| 日韩视频在线欧美| 最近最新中文字幕大全免费视频| 久久久久久亚洲精品国产蜜桃av| 18在线观看网站| 1024视频免费在线观看| 免费av中文字幕在线| 欧美精品啪啪一区二区三区 | 欧美老熟妇乱子伦牲交| 别揉我奶头~嗯~啊~动态视频 | 99精品欧美一区二区三区四区| 国产av又大| 50天的宝宝边吃奶边哭怎么回事| 亚洲专区中文字幕在线| 日韩大片免费观看网站| 在线 av 中文字幕| 精品亚洲成a人片在线观看| 最近中文字幕2019免费版| 午夜激情久久久久久久| 国产精品一二三区在线看| 美女高潮到喷水免费观看| av片东京热男人的天堂| 亚洲自偷自拍图片 自拍| 丝袜喷水一区| 日韩欧美免费精品| 亚洲成国产人片在线观看| 日韩免费高清中文字幕av| 女人被躁到高潮嗷嗷叫费观| 免费高清在线观看日韩| 日韩大片免费观看网站| 老司机影院成人| 午夜免费观看性视频| 在线看a的网站| 精品一区在线观看国产| 啦啦啦啦在线视频资源| 国产免费现黄频在线看| 欧美精品一区二区免费开放| 免费在线观看日本一区| 天天躁夜夜躁狠狠躁躁| 最近中文字幕2019免费版| 久久久国产一区二区| 久久久久国产精品人妻一区二区| 国产一卡二卡三卡精品| 成人黄色视频免费在线看| 国产免费一区二区三区四区乱码| 欧美av亚洲av综合av国产av| 交换朋友夫妻互换小说| 色视频在线一区二区三区| 自线自在国产av| 狠狠狠狠99中文字幕| 国产亚洲欧美精品永久| 色婷婷久久久亚洲欧美| 国产主播在线观看一区二区| 91国产中文字幕| 99国产精品99久久久久| 丁香六月天网| cao死你这个sao货| 亚洲avbb在线观看| av免费在线观看网站| 男女床上黄色一级片免费看| 热99re8久久精品国产| 亚洲国产欧美网| 午夜影院在线不卡| 国产成人影院久久av| 精品少妇内射三级| 欧美97在线视频| 黄片大片在线免费观看| 99国产精品99久久久久| 久久国产精品男人的天堂亚洲| 中国美女看黄片| 啦啦啦中文免费视频观看日本| 精品亚洲成a人片在线观看| 久久久久久久久久久久大奶| 国产亚洲午夜精品一区二区久久| 日本91视频免费播放| 欧美日韩一级在线毛片| 悠悠久久av| 91国产中文字幕| 少妇 在线观看| 永久免费av网站大全| 91精品国产国语对白视频| 亚洲久久久国产精品| 国产精品一区二区在线不卡| 少妇被粗大的猛进出69影院| 国产片内射在线| 国产淫语在线视频| 欧美中文综合在线视频| 美女午夜性视频免费| 国产97色在线日韩免费| av线在线观看网站| 国产精品秋霞免费鲁丝片| 12—13女人毛片做爰片一| 欧美日韩成人在线一区二区| 激情视频va一区二区三区| 麻豆av在线久日| 女性被躁到高潮视频| 91av网站免费观看| 不卡av一区二区三区| 久久人妻福利社区极品人妻图片| 丰满少妇做爰视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成人手机| 99精国产麻豆久久婷婷| 亚洲欧美一区二区三区久久| 久久国产精品大桥未久av| 男男h啪啪无遮挡| 黄色视频不卡| 50天的宝宝边吃奶边哭怎么回事| 十八禁网站网址无遮挡| 一区二区三区四区激情视频| 女人久久www免费人成看片| av在线播放精品| 日本黄色日本黄色录像| 免费一级毛片在线播放高清视频 | 亚洲五月婷婷丁香| 亚洲精品国产区一区二| 美女主播在线视频| 自线自在国产av| 操出白浆在线播放| 成年av动漫网址| 日本av免费视频播放| 成年人午夜在线观看视频| 亚洲欧美色中文字幕在线| 午夜免费成人在线视频| 俄罗斯特黄特色一大片| 成人三级做爰电影| 飞空精品影院首页| 18禁黄网站禁片午夜丰满| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩一级在线毛片| 日韩欧美国产一区二区入口| 日韩大片免费观看网站| 欧美97在线视频| 一级片'在线观看视频| 精品一区二区三卡| 多毛熟女@视频| 精品视频人人做人人爽| 亚洲国产中文字幕在线视频| 80岁老熟妇乱子伦牲交| 中文字幕最新亚洲高清| www.999成人在线观看| 久久久久久久精品精品| videosex国产| 亚洲国产欧美网| 大陆偷拍与自拍| 美国免费a级毛片| 国产亚洲精品久久久久5区| 国产精品久久久久成人av| 欧美激情高清一区二区三区| 男人爽女人下面视频在线观看| 满18在线观看网站| 国产精品免费视频内射| 青草久久国产| 在线亚洲精品国产二区图片欧美| 日韩一卡2卡3卡4卡2021年| 男女午夜视频在线观看| 色综合欧美亚洲国产小说| 精品视频人人做人人爽| 久久午夜综合久久蜜桃| 亚洲国产日韩一区二区| 中文字幕高清在线视频| 最近中文字幕2019免费版| 午夜福利影视在线免费观看| 人成视频在线观看免费观看| 久热爱精品视频在线9| 十八禁人妻一区二区| 啦啦啦啦在线视频资源| 人人妻人人澡人人看| 成人国产一区最新在线观看| 桃花免费在线播放| 国产淫语在线视频| 91大片在线观看| 男女床上黄色一级片免费看| 在线观看免费高清a一片| 亚洲人成77777在线视频| 一本综合久久免费| 天天操日日干夜夜撸| 9热在线视频观看99| 久久精品成人免费网站| 欧美日韩黄片免| 69精品国产乱码久久久| 新久久久久国产一级毛片| 在线av久久热| 两个人免费观看高清视频| 中国国产av一级| 女警被强在线播放| 亚洲少妇的诱惑av| 免费久久久久久久精品成人欧美视频| 久久亚洲国产成人精品v| 我的亚洲天堂| 女性生殖器流出的白浆| 老司机深夜福利视频在线观看 | 国产免费视频播放在线视频| 欧美激情高清一区二区三区| 狂野欧美激情性bbbbbb| 久9热在线精品视频| 午夜精品国产一区二区电影| 99久久精品国产亚洲精品| 日本91视频免费播放| 欧美另类亚洲清纯唯美| 人妻一区二区av| 国产av国产精品国产| 天天添夜夜摸| 无限看片的www在线观看| 交换朋友夫妻互换小说| 999久久久精品免费观看国产| 18禁黄网站禁片午夜丰满| 在线观看人妻少妇| 成年女人毛片免费观看观看9 | 一二三四在线观看免费中文在| 免费人妻精品一区二区三区视频| 在线亚洲精品国产二区图片欧美| 一本色道久久久久久精品综合| h视频一区二区三区| 久久国产精品人妻蜜桃| 成人亚洲精品一区在线观看| 免费日韩欧美在线观看| 一级,二级,三级黄色视频| 亚洲精品国产区一区二| 精品人妻在线不人妻| 夫妻午夜视频| 精品少妇黑人巨大在线播放| 咕卡用的链子| 天天操日日干夜夜撸| 建设人人有责人人尽责人人享有的| 欧美日韩视频精品一区| 精品亚洲成国产av| 国产片内射在线| 日韩三级视频一区二区三区| 日韩免费高清中文字幕av| 电影成人av| 999精品在线视频| 性少妇av在线| 国产又爽黄色视频| 人人澡人人妻人| 老司机影院毛片| 五月天丁香电影| 亚洲一区二区三区欧美精品| 飞空精品影院首页| 精品久久蜜臀av无| 亚洲精品一二三| 99国产精品99久久久久| 亚洲一码二码三码区别大吗| 国产精品久久久久久精品古装| 日本vs欧美在线观看视频| 亚洲第一av免费看| 中国国产av一级| 色94色欧美一区二区| 欧美激情久久久久久爽电影 | 免费在线观看日本一区| 肉色欧美久久久久久久蜜桃| 亚洲熟女毛片儿| 国产极品粉嫩免费观看在线| 美女午夜性视频免费| 国产精品一区二区免费欧美 | 久久国产精品影院| 99久久国产精品久久久| 欧美97在线视频| 淫妇啪啪啪对白视频 | 啦啦啦视频在线资源免费观看| 国产亚洲精品久久久久5区| 黑人巨大精品欧美一区二区mp4| 成人手机av| 女人高潮潮喷娇喘18禁视频| 国产高清videossex| 国产精品.久久久| 欧美性长视频在线观看| 香蕉国产在线看| 老熟女久久久| 黑人猛操日本美女一级片| 国产一区二区三区综合在线观看| 不卡av一区二区三区| 久久精品亚洲av国产电影网| 99国产极品粉嫩在线观看| 亚洲avbb在线观看| 亚洲精品国产精品久久久不卡| 欧美国产精品va在线观看不卡| 这个男人来自地球电影免费观看| 亚洲伊人色综图| 国产激情久久老熟女| 精品少妇黑人巨大在线播放| 热99国产精品久久久久久7| 久久性视频一级片| 少妇精品久久久久久久| 搡老乐熟女国产| 成在线人永久免费视频| 国产日韩欧美视频二区| 国产成人免费观看mmmm| 亚洲视频免费观看视频| 在线观看舔阴道视频| 久久久久国内视频| 女人高潮潮喷娇喘18禁视频| 搡老乐熟女国产| 国产精品自产拍在线观看55亚洲 | 一边摸一边抽搐一进一出视频| 国产熟女午夜一区二区三区| 欧美精品高潮呻吟av久久| 精品一区二区三区av网在线观看 | 亚洲精品日韩在线中文字幕| 国产精品影院久久| www.熟女人妻精品国产| 女人久久www免费人成看片| 丝袜脚勾引网站| 日韩有码中文字幕| avwww免费| 久久国产亚洲av麻豆专区| 一级片免费观看大全| www.自偷自拍.com| 亚洲欧洲精品一区二区精品久久久| videosex国产| 亚洲五月色婷婷综合| 日韩一卡2卡3卡4卡2021年| 91麻豆精品激情在线观看国产 | 亚洲精品日韩在线中文字幕| 新久久久久国产一级毛片| 丰满迷人的少妇在线观看| 午夜福利在线免费观看网站| 99国产精品一区二区三区| 亚洲av欧美aⅴ国产| 国产日韩欧美视频二区| 日本av免费视频播放| 999精品在线视频| 亚洲天堂av无毛| 老汉色∧v一级毛片| 精品国产乱码久久久久久小说| www.自偷自拍.com| 桃花免费在线播放| 精品国产一区二区三区四区第35| 丰满人妻熟妇乱又伦精品不卡| 18禁国产床啪视频网站| 久久综合国产亚洲精品| 性高湖久久久久久久久免费观看| 天天躁日日躁夜夜躁夜夜| 日本a在线网址| 欧美精品亚洲一区二区| 黄片大片在线免费观看| 女人被躁到高潮嗷嗷叫费观| 国产一区二区在线观看av| 亚洲第一av免费看| 欧美少妇被猛烈插入视频| 99精品欧美一区二区三区四区| 国产免费视频播放在线视频| 91av网站免费观看| 97人妻天天添夜夜摸| 两性夫妻黄色片| 成人国产一区最新在线观看| 中文字幕色久视频| 日韩欧美一区二区三区在线观看 | 国产精品.久久久| 黑人巨大精品欧美一区二区mp4| 日韩大片免费观看网站| 一二三四社区在线视频社区8| 免费日韩欧美在线观看| 黄色片一级片一级黄色片| 91麻豆精品激情在线观看国产 | 亚洲熟女精品中文字幕| 欧美黑人精品巨大| 王馨瑶露胸无遮挡在线观看| 中文欧美无线码| 精品人妻一区二区三区麻豆| 亚洲欧美日韩高清在线视频 | 高清在线国产一区| 少妇裸体淫交视频免费看高清 | 一区在线观看完整版| 亚洲国产看品久久| 男男h啪啪无遮挡| 国产精品成人在线| avwww免费| 天天躁日日躁夜夜躁夜夜| 国产熟女午夜一区二区三区| 欧美变态另类bdsm刘玥| 国产亚洲精品一区二区www | 久久性视频一级片| 久久狼人影院| 亚洲av成人一区二区三| 国产精品成人在线| 午夜免费观看性视频| 日韩免费高清中文字幕av| 欧美日韩福利视频一区二区| 肉色欧美久久久久久久蜜桃| 欧美人与性动交α欧美精品济南到| 水蜜桃什么品种好| 久久性视频一级片| videos熟女内射| 国产人伦9x9x在线观看| 国产在线免费精品|