舒蕊 劉健 章玉萍 張麗麗 陳明 代君君 范濤

（安徽省農(nóng)業(yè)科學院蠶桑研究所，合肥 230031）

蜜蜂是昆蟲膜翅目家族中的重要成員，不僅是重要的農(nóng)業(yè)經(jīng)濟昆蟲，也是科學研究中的模式生物。我國飼養(yǎng)的西方蜜蜂品種中，黃色蜂種占有較大比重，特別是以意大利蜜蜂為主的黃色蜂種。意蜂群中的個體表現(xiàn)出不同程度黃體色：雄蜂為黃色到暗棕色，絨毛淡黃色尤為明顯；工蜂的腹部具兩個黃色環(huán)帶的個體居多，絨毛淡黃色；蜂王的腹部多為黃色到暗棕色，尾部黑色，只有少數(shù)蜂王通體都是黃色。黃色蜂種及其雜交種遍布全國各地，在養(yǎng)蜂生產(chǎn)中發(fā)揮著重要作用[1]。然而，意蜂黃體色的形成機制一直未解。眾所周知，昆蟲黃色色素的形成源自蝶呤類色素，四氫生物蝶呤（BH4）合成途徑與蝶呤類色素的合成與代謝息息相關。BH4是芳香族氨基酸羥化酶的必需輔酶，也是一氧化氮合酶（NOS）的重要輔助因子，還能作為NOS的調控因子或直接作為自由基清除劑參與氧化應激作用[2,3]。BH4的合成不足會導致神經(jīng)遞質的缺乏以及細胞內皮功能障礙，引發(fā)多種神經(jīng)性生理代謝疾病，臨床醫(yī)學上把口服BH4作為治療其中某些疾病的有效方法[4-6]。現(xiàn)階段，BH4合成代謝網(wǎng)絡在雙翅目模式生物黑腹果蠅和鱗翅目模式生物家蠶中都有較多的研究，但是在以蜜蜂為代表的膜翅目昆蟲中幾乎沒有研究報道。研究蜜蜂體內BH4合成途徑對于意蜂乃至膜翅目的色素合成與代謝至關重要。

6-丙酮酰四氫蝶呤合酶（6-pyruvoyl tetrahydrobiopterin synthase，PTS）是BH4合成代謝途徑中的關鍵酶。PTS催化BH4從頭合成途徑中的第二步反應，即催化二氫新蝶呤三磷酸轉變?yōu)?-丙酮酰四氫蝶呤（PTP），該反應依賴于Mg2+和Zn2+[7-9]。人類和哺乳動物的PTS基因結構是高度保守的，它們的基因組序列有6～7kbp，包含6個外顯子。研究表明，PTS酶的功能異常導致BH4合成受阻，從而引起人類罹患BH4缺乏癥。此類患者體內分離發(fā)現(xiàn)了大量的PTS基因突變的位點，其中44個突變位點分布在PTS基因所有的6個外顯子和前面3個內含子的區(qū)域內。除了大量的突變位點，不同物種的PTS基因表達與調控方式也有不同。比如，果蠅的PTS基因轉錄了兩種類型的PTS mRNA：一類是受近端啟動子作用的，只在頭部特異出現(xiàn)，具有3個外顯子的序列；另一類是受遠端啟動子作用的，在5’端非翻譯區(qū)多了一段新的內含子序列。但是這兩類mRNA在果蠅體內表達的ORF序列是一樣的[10]。

生物信息學方法在功能基因的測序、識別、功能預測、結構分析等方面都發(fā)揮著重要作用。自意蜂基因組全序列測序完成后，越來越多的生物信息學方法應用到意蜂各類基因的研究工作中。生物信息學分析將促進人們對于意蜂6-丙酮酰四氫蝶呤合酶的了解和認識。

本文選擇意蜂6-丙酮酰四氫蝶呤合酶，采用生物信息學方法，對其相應氨基酸序列的理化性質、疏水性、跨膜區(qū)、結構功能域以及空間二級結構和三級結構等方面進行分析和預測，并構建系統(tǒng)進化樹，旨在為日后深入研究意蜂6-丙酮酰四氫蝶呤合酶家族的結構和功能解析及應用等方面提供有力幫助和理論依據(jù)。

1 材料與方法

1.1 序列來源

來自NCBI數(shù)據(jù)庫（http://www.ncbi.nlm.nih.gov/）的意大利蜜蜂6-丙酮酰四氫蝶呤合酶（序列號XM-623992）。

1.2 方法

使用ProtScale分析蛋白質疏水性，使用ORFFinder和ProtParam對編碼意蜂6-丙酮酰四氫蝶呤合酶的氨基酸組分和理化性質進行分析，使用網(wǎng)站TMHMM預測蛋白質跨膜區(qū)，結構功能域分析使用SMART服務器，蛋白質二級結構預測使用PORTER完成，蛋白質三級建模使用SWISS-MODEL完成，使用軟件MEGA 6.0構建系統(tǒng)進化樹。

2 結果與分析

2.1 理化性質分析

使用ORF Finder和ProtParam分析意蜂6-丙酮酰四氫蝶呤合酶的蛋白組分和理化性質(表1)。意蜂PTS基因全長829bp，含有一個429bp的開放閱讀框，編碼142個氨基酸，氨基酸組成中，不含Pyl和Sec，Lys含量最高，達到11.3%。含有20個酸性氨基酸（Asp+Glu），20個堿性氨基酸（Arg+Lys），氨基酸序列的不穩(wěn)定系數(shù)為33.06，說明此蛋白是穩(wěn)定的。預測蛋白分子量為16.66027kDa，理論等電點（pI）為6.96，原子組成為C745H1179N195O219S9，脂肪族氨基酸指數(shù)為89.08，親水性指數(shù)為-0.499，表明其為親水性蛋白。

2.2 疏水性分析

圖1 意蜂6-丙酮酰四氫蝶呤合酶氨基酸序列的疏水性分析

ExPASy的ProtScale是常用的預測蛋白質疏水性的分析方法。使用ProtScale工具在線分析意蜂6-丙酮酰四氫蝶呤合酶的疏水性（圖1），其計算方法是基于K-D法的蛋白質疏水性。6-丙酮酰四氫蝶呤合酶氨基酸序列中負值的親水性氨基酸的比例高于正值的疏水性氨基酸，因而判斷其為親水性蛋白，與理化性質分析結果一致。

表1 意蜂6-丙酮酰四氫蝶呤合酶基因蛋白組分和理化性質

2.3 跨膜區(qū)分析

使用TMHMM工具對意蜂6-丙酮酰四氫蝶呤合酶進行跨膜區(qū)分析（圖2）?；赥MHMM計算的后驗概率，6-丙酮酰四氫蝶呤合酶的142個氨基酸都在膜外，該酶是一個典型的膜外蛋白，不是跨膜蛋白質。

圖2 意蜂6-丙酮酰四氫蝶呤合酶氨基酸序列的跨膜區(qū)分析

2.4 結構功能域分析

分析結構功能域對于蛋白質結構的分類和預測非常重要，使用SMART服務器進行意蜂6-丙酮酰四氫蝶呤合酶的結構功能域分析（圖3）。意蜂6-丙酮酰四氫蝶呤合酶在序列的7～140位點含有一個6-丙酮酰四氫蝶呤合成酶結構域。

圖3 意蜂6-丙酮酰四氫蝶呤合酶結構功能域分析

2.5 二級結構預測

蛋白質的二級結構預測既是聯(lián)系一級結構和三級結構的紐帶，也為后續(xù)蛋白質三級結構的預測提供信息。使用PORTER服務器預測意蜂6-丙酮酰四氫蝶呤合酶蛋白質的二級結構（圖4）。意蜂6-丙酮酰四氫蝶呤合酶主要由4種二級結構組成，預測結果顯示58個氨基酸為α-螺旋，占比40.85%，其次是無規(guī)卷曲共41個氨基酸，占比28.87%，35個氨基酸是β-折疊，共占24.65%，β-轉角有8個氨基酸，占5.63%。

圖4 意蜂6-丙酮酰四氫蝶呤合酶二級結構預測

2.6 三級結構預測

圖5 意蜂6-丙酮酰四氫蝶呤合酶三級結構模擬

蛋白質的生物學功能在很大程度上是由其空間結構決定的，將蛋白質的高級結構分析清楚是理解其功能的關鍵。通過SWISS-MODEL同源建模意蜂6-丙酮酰四氫蝶呤合酶蛋白的三級結構，一般相似性超過50%即可獲得較精確的結構模型（ID：3i2b.1.A，序列一致性56.52%），同源建模如圖5。

2.7 系統(tǒng)進化樹分析

系統(tǒng)進化樹發(fā)育分析是生物信息學中的重要研究領域，通過它的計算分析可以判斷物種之間的親緣關系。本研究先在NCBI網(wǎng)站中搜索PTS序列信息，然后分別選擇膜翅目、雙翅目、鱗翅目、鞘翅目共15種昆蟲的PTS蛋白序列，構建系統(tǒng)進化樹（圖6）。分析結果表明，意大利蜜蜂與小蜜蜂、大蜜蜂以及中華蜜蜂的進化關系最密切。該結果為后續(xù)開展意蜂6-丙酮酰四氫蝶呤合酶家族功能基因的實驗研究提供了參考信息。

圖6 15種昆蟲6-丙酮酰四氫蝶呤合酶的系統(tǒng)進化樹

3 討論

蜜蜂在生物學上是全變態(tài)昆蟲，在農(nóng)業(yè)上是重要的授粉昆蟲，在科研中亦是重要的模式生物。目前對昆蟲6-丙酮酰四氫蝶呤合酶基因的研究多集中在果蠅和家蠶等其他種類的昆蟲上，而膜翅目意蜂6-丙酮酰四氫蝶呤合酶基因的研究卻鮮有報道。意蜂6-丙酮酰四氫蝶呤合酶基因在蛋白質結構和生物學功能方面相對于其他昆蟲是否有特別之處，需要科學工作者們做更多的研究與探索。本研究利用生物信息學方法對意大利蜜蜂6-丙酮酰四氫蝶呤合酶基因進行了系統(tǒng)分析，結果表明，意蜂6-丙酮酰四氫蝶呤合酶基因全長829bp，含有一個429bp的開放閱讀框，編碼142個氨基酸，屬于親水性蛋白；預測蛋白分子量為16.66027kD，理論等電點（pI）為6.96，原子組成為C745H1179N195O219S9；二級結構主要由α-螺旋構成；在序列的7～140位點含有一個6-丙酮酰四氫蝶呤合成酶結構域，含有相應的催化結構域和底物結合域；蛋白的二級結構主要由α-螺旋組成，其次是無規(guī)卷曲和β-折疊，同源建模了蛋白質的三級結構；使用15種6-丙酮酰四氫蝶呤合成酶氨基酸序列構建系統(tǒng)進化樹，分析結果與傳統(tǒng)的生物學分類結果一致。

本文的研究結果為深入開展意蜂6-丙酮酰四氫蝶呤合成酶的科學研究提供參考信息，為意蜂6-丙酮酰四氫蝶呤合成酶基因家族的分離、純化及克隆分析提供參考，也為進一步研究其高級結構與功能關系提供理論依據(jù)。本文是基于生物信息學工具對目標基因進行的系統(tǒng)研究，關于意蜂6-丙酮酰四氫蝶呤合成酶的具體生物學功能和表達機制等還需進一步研究。