林非凡，韓雅莉

(廣東工業(yè)大學(xué)輕工化工學(xué)院，廣東廣州510006)

細胞色素P450(Cytochrome P450，CYP450)是生物體內(nèi)一組結(jié)構(gòu)和功能密切相關(guān)的超家族基因編碼的同工酶，該酶系對內(nèi)、外源化合物的轉(zhuǎn)化、失活及活化起著至關(guān)重要的作用，包括參與類固醇激素、維生素、脂肪酸衍生物等內(nèi)源性物質(zhì)的合成和降解［1］，外來化合物如藥物、環(huán)境污染物、致癌物的代謝等［2］.CYP2E1是CYP450酶系的重要成員之一，其編碼的二甲基亞硝胺D-脫甲基酶是極其重要的I相代謝酶，其主要參與亞硝胺及其前體物和低分子量鹵代氰類化合物在體內(nèi)的代謝，也參與苯、乙醇、四氯化碳和煙草中的3’羥基吡啶等多種低分子量化合物的代謝［3-4］，特別是參與多種前致癌物和前毒物，如亞硝胺、苯胺等物質(zhì)的轉(zhuǎn)化為最終致癌活性物質(zhì)，因此，在多種癌癥易感性中起著重要作用.CYP2E1被認為是機體對環(huán)境、工業(yè)毒物或致癌物敏感程度的重要決定因素［5-6］.

本文采用生物信息學(xué)方法，對人CYP2E1基因序列特征、編碼的蛋白結(jié)構(gòu)、理化性質(zhì)、空間結(jié)構(gòu)、酶的活性中心和系統(tǒng)發(fā)育進行預(yù)測和分析，從分子水平揭示CYP2E1的作用機制，為其轉(zhuǎn)化環(huán)境化合物及致癌機理提供理論依據(jù)，也為污染物的防治提供新的思路與途徑.

1 材料與方法

1.1 序列來源

人CYP2E1基因及蛋白的GenBank的登錄號分別為NC000010.11和NM000773.3.

1.2 分析方法

利用NCBI(http:∥www.ncbi.nlm.nih.gov/blast)找到CYP2E1的開放閱讀框及編碼蛋白.利用DNAStar軟件的Editseq程序，對該開放閱讀框編碼氨基酸進行翻譯，并分析其氨基酸序列;利用Prot-Param(http:∥www.expasy.ch/tools/prot-param.html)在線工具預(yù)測人CYP2E1各理化性質(zhì);使用DNAStar軟件的Megalingn程序?qū)θ?、猴、犬、牛、豬、小鼠、大鼠的CYP2E1編碼蛋白序列進行系統(tǒng)進化分析;利用Tampered(www.ch.embnet.org/software/tepred-form.html)在線工具分析人CYP2E1的跨膜結(jié)構(gòu);采用SignalP(http:∥www.cbs.dtu.dk/services/SignalP)在線工具分析CYP2E1編碼蛋白的信號肽.利用DNAStar軟件中Protean程序和Chou-Fasman法(http:∥www.predictprotein.org/)預(yù)測人CYP2E1二級結(jié)構(gòu);利用WOLFPsort軟件預(yù)測人CYP2E1的亞細胞定位;利用SWISS-MODEL(http:∥swissmodel.expasy.org/)軟件預(yù)測CYP2E1三級結(jié)構(gòu)和酶的活性中心;采用ExPASY的ProtScale(http:∥expasy.org/tools/protscale.html)軟件分析人CYP2E1的親(疏)水性.

2 結(jié)果與分析

2.1 人CYP2E1基因基本信息

人CYP2E1基因定位于第10號染色體(10q24.3)，其上游為2 788 bp，下游為599 bp，由8個內(nèi)含子、9個外顯子和典型的TATA盒組成，其大小為11 413 bp，cDNA長度為1 497 bp［7］.該基因堿基組成:A=21.1%［397］，G=18.6%［350］，T=19.2%［361］，C=41.1%［362］，Ambigous=0.00%［0］，(A+T)=40.3%［758］，(C+G)=59.7%［712］.顯然，4種堿基比例不平衡.CG含量顯著大于AT含量.CYP2E1基因中存在2個多態(tài)位點，分別位于其5’端和內(nèi)含子6中［8］.

2.2 人CYP2E1的主要理化性質(zhì)和親(疏)水性預(yù)測

經(jīng)DNAStar分析人CYP2E1編碼的功能蛋白由493個氨基酸組成，相對分子質(zhì)量約為172 386.6 Da.氨基酸組成:堿性氨基酸15個(K，R)，占3.1%;酸性氨基酸14個(D，E)，占2.9%;疏水性氨基酸41個(I，L，F(xiàn)，W，V)，占8.3%;極性氨基酸210個(N，C，Q，S，T，Y)，占42.5%.等電點理論值7.04，表明人CYP2E1蛋白的極性氨基酸與疏水性氨基酸占有較大比例.經(jīng)ProtParam軟件分析，推測人CYP2E1的分子式為C6778H10986N2074O2429S374;在波長280 nm時吸光度(A值)為1.6;其不穩(wěn)定系數(shù)為44.09;其脂肪族氨基酸系數(shù)為43.79;其總平均親水性是0.583.該蛋白氨基酸組成中Gys、Gly、Thr、Leu含量較高，依次為6.2%、5.8%、4.8%和11.8%，而Trp含量僅為1.0%;其親(疏)水性如圖1所示，存在8個高疏水性區(qū)域，分別分布于8～23、64～81、168～185、209～214、286～300、308～317、381～396、437～459氨基酸位點處;主要的最小分值區(qū)域位于139、140、141、142、143、144、145和146這8個氨基酸位點，該區(qū)域為高親水性區(qū)域.因此，推測該蛋白是高疏水性蛋白.

圖1 人CYP2E1疏水性分析Fig.1 Hydrophobicity analysis of CYP2E1

2.3 人CYP2E1跨膜結(jié)構(gòu)分析

預(yù)測與分析蛋白跨膜結(jié)構(gòu)域?qū)τ跍?zhǔn)確了解蛋白質(zhì)分子結(jié)構(gòu)與功能以及在細胞中的定位具有重要意義［9］.本文采用Tampered軟件分析了CYP2E1跨膜結(jié)構(gòu)，結(jié)果如圖2所示.根據(jù)預(yù)測模型，編碼蛋白產(chǎn)生了6次高峰，分別在44～62、99～118、205～223、248～264、327～347以及480～499氨基酸位點處，推測該蛋白存在6個跨膜區(qū)域，即一條多肽鏈跨膜折疊，形成由6個α螺旋段反平行裝配的球狀膜蛋白.推測該蛋白是跨膜蛋白.

圖2 人CYP2E1的跨膜區(qū)域分析Fig.2 Analysis of transmembrane domains of CYP2E1

2.4 人CYP2E1信號肽分析

采用SignalP3.0 Server軟件對CYP2E1進行信號肽預(yù)測，結(jié)果如圖3所示，C列為剪切位點得分，S列為信號肽得分，Y列為綜合剪切位點得分，典型信號肽的C值和Y值趨向于+1，S值在剪切位點之前高而在剪切位點之后變低，當(dāng)C、Y、S值計算結(jié)論是YES時表明預(yù)測結(jié)果較好［10］.分析結(jié)果表明，CYP2E1肽段分值曲線典型，C值為0.831，Y值為0.659，S值為0.653，C，Y，S均為最大值，其切割點位于第28位和第29氨基酸之間，該蛋白的信號肽位于第1-28位氨基酸.

2.5 人CYP2E1二級結(jié)構(gòu)與亞細胞定位預(yù)測

對于未知功能蛋白質(zhì)的分子結(jié)構(gòu)和結(jié)構(gòu)域或功能單位進行確認，有助于了解該蛋白的生物學(xué)功能［11］.本文采用Protean程序分析CYP2E1二級結(jié)構(gòu)，結(jié)果如圖4所示.從圖4可以看出該蛋白中有6段 α-螺旋，分布于9-26，34-79，92-96，110-202，251-279，375-389等氨基酸殘基位點處;有6處β-折疊，分布于210-213，233-237，295-326，349-379，386-401，402-405位點;還有32處 β-轉(zhuǎn)角，分布于30-33，63-159，193-221，230-261，276-303，318-367，382-396，408-411等位點.人CYP2E1有典型的細胞色素P450二級結(jié)構(gòu)，在二級結(jié)構(gòu)基礎(chǔ)上進一步發(fā)生緊密的卷曲與折疊，形成三維空間結(jié)構(gòu)，這樣有利于被催化底物可以沿著β1折疊進入催化中心部位.采用WOLFPsort軟件對預(yù)測該蛋白亞細胞定位，判斷該蛋白質(zhì)主要定位于膜結(jié)構(gòu)上，最可能分布于內(nèi)質(zhì)網(wǎng)膜、過氧化物酶體、溶酶體與高爾基體等膜上.

圖3 人CYP2E1信號肽預(yù)測Fig.3 Prediction of signal peptide of CYP2E1

圖4 人CYP2E1二級結(jié)構(gòu)圖Fig.4 Secondary structure of CYP2E1

2.6 人CYP2E1三級結(jié)構(gòu)和活性中心預(yù)測

經(jīng)SWISS-MODEL軟件模擬人CYP2E1的高級結(jié)構(gòu)和活性中心如圖5和圖6所示.預(yù)測出CYP2E1活性中心的體積約為0.19 nm3，外表為小球狀，整個分子有6個α-螺旋和1個β折疊，Asn202殘基、Ser298殘基以及Phe205殘基有可能與鐵卟啉環(huán)形成活性中心，CYP2E1蛋白血紅素鐵和底物催化位點之間的距離在0.3～0.6 nm之間，在空間上和底物結(jié)合為三明治式的結(jié)構(gòu)，參與催化反應(yīng).

圖5 人CYP2E1三級預(yù)測結(jié)構(gòu)Fig.5 Tertiary structure of CYP2E1

圖6 人CYP2E1的預(yù)測活性中心Fig.6 Prediction of the active center of CYP2E1

2.7 人CYP2E1的分子進化分析

采用Megalingn程序，將人與6物種(猴、犬、牛、豬、小鼠和大鼠)的CYP2E1氨基酸序列進行同源性分析，并構(gòu)建系統(tǒng)發(fā)育樹，如圖7所示.人類與猴相似度為96%;與大鼠、小鼠相似度均為80%，與牛相似度為71.2%，與犬相似度為72%，與豬相似度為70%.說明CYP2E1蛋白在進化上具有高度保守性.

圖7 人與6物種的系統(tǒng)發(fā)育樹Fig.7 Phylogenetic trees of human being and other 6 species

3 結(jié)論

本文通過對CYP2E1的一般性質(zhì)、親疏水性、二級結(jié)構(gòu)、三級結(jié)構(gòu)、活性中心、跨膜結(jié)構(gòu)、信號肽以及亞細胞定位等進行預(yù)測分析，依分析結(jié)果推測人CYP2E1處在細胞的膜結(jié)構(gòu)上，可能是膜結(jié)合蛋白，而且是一種跨膜蛋白，即一部分親水肽段在膜外，一部分疏水肽段埋在磷脂雙分子中與膜共價結(jié)合，催化功能位點則處于親水肽段部位［12］.本文還預(yù)測了信號肽剪切位點，表明該基因表達具有信號肽，推測該蛋白在信號肽的引導(dǎo)下通過內(nèi)質(zhì)網(wǎng)膜進入腔內(nèi)進行修飾成熟，然后整合在質(zhì)膜上.對CYP2E1酶活性中心分析表明酶活性中心體積較小，僅為0.19 nm3，這可能與代謝小分子化合物有關(guān);CYP2E1與底物結(jié)合形成類似三明治式結(jié)構(gòu)，如:其分子的Phe205和Phe293殘基可與氯唑沙宗的苯并唑環(huán)形成π－π共軛，Asp202殘基與氯唑沙宗羥基、Ser298殘基與氯唑沙宗氯化基團均以氫鍵形式結(jié)合.活性中心的這種結(jié)構(gòu)構(gòu)象也揭示了CYP2E1可與多種底物(如丁二烯、苯胺、甲苯、二甲基亞硝胺、乙醇、月桂酸等)結(jié)合的機理.

通過對人CYP2E1氨基酸序列與其他6種哺乳動物該蛋白序列的比較和聚類研究，表明CYP2E1具有高度保守性，人與猴的相似度達95%，與嚙齒動物小鼠和大鼠的相似度為80%，說明在生物進化過程中，該基因的保守序列變化很小，以至于其表達產(chǎn)物的生物活性類似，提示我們在開展CYP2E1的毒理學(xué)及藥理學(xué)方面的研究時，可以考慮以其他實驗動物(如小鼠、大鼠等)和其細胞系作為模型實驗材料，深入探討其作用特點與機理.

［1］Shark A H，Crawford M A，Reifen R.Update on alphalinelic acid［J］.Nutrition Reviews，2008，66(6):326-332.

［2］Brunner G.Mass transferfrom solid mater in gas extraction［J］.Ber Bunsenged Phys Chem，1984.887-891.

［3］Dwivedi C，Natarajan K，Matthees D P.Chemmopreventive effects of dietary flaxseed oil on colon tumor development［J］.Nutrition Cancer，2005，51(1):52-58.

［4］Yaniv Z.Schafferman D.Shamir I，et al.Cholesterol and triglyceride reduction in rats fed Matthiola incana seed oil rich in(n-3)fatty acids［J］.J Agric Food Chem，1999，47(2):637-642.

［5］Nguyen A D，Le Goffic F.Limits of wiped film short-path distiller［J］.Chemical Engineering Science，1997，52(16):2661-2666.

［6］Cvengros J，Pollak S，Micov M，et al.Film wiping in the molecular evapotator［J］.Chemical Engineering Journal，2001，81(1):9-14.

［7］David B G.A theoretical and experimental study of the centrifugal molecular still［J］.Aiche J，1972，18(2):269-276.

［8］Christon R，Haloue R，Durand G.Dietary polyunsaturated fatty acids and aging modylate lutathione-ralated antioxidants in rat liver［J］.J Nutr，1995，125:3062.

［9］Rui L M.Supercritical CO2extraction of γ-linolenic acid(GLA)from the cyanobacterium Arthrospira(Spirulina)maxma experinents and modeling［J］.Chemical Engineering Joumal，2005，10(5):147-152.

［10］Guil Guerrero J L，Campra-Madrid P，Belarbi E H.γ-Linolenic acid purification from seed oil sources by argentated silica gel chromatography colum［J］.Process Biochemistry 2000，36(1):341-354.

［11］Baba N H，Antoniades K，Habbal Z.Effects of dietary canola，olive，and lino lenic acid enriched olive oils on plasma lipids，lipid peroxidation and lipoprotein lipase activity in rats［J］.Nutrition Research，1999，19(4):601-612.

［12］Xu W K，Mong L S，Li Z Y.Isolation and identification of a bitter constituent from Luo-han-guo［J］.Sunripe Fruits Guihaia，1992(12):136-138.