唇形科植物牻牛兒基牻牛兒基焦磷酸合酶編碼基因及其氨基酸序列的生物信息學(xué)分析

陳宜均+榮齊仙+姜丹+申業(yè)+黃璐琦

[摘要]牻牛兒基牻牛兒基焦磷酸合酶（geranylgeranyl pyrophosphate synthase，GGPS）是二萜類物質(zhì)合成途徑中的一個(gè)關(guān)鍵酶。該文采用生物信息學(xué)方法對(duì)9種唇形科植物的牻牛兒基牻牛兒基焦磷酸合酶的核苷酸及其編碼的氨基酸序列的結(jié)構(gòu)、理化性質(zhì)、信號(hào)肽、導(dǎo)肽、跨膜結(jié)構(gòu)域、疏水性/親水性和亞細(xì)胞定位特征、二級(jí)結(jié)構(gòu)、蛋白質(zhì)功能域、三級(jí)結(jié)構(gòu)和進(jìn)化關(guān)系進(jìn)行了初步預(yù)測(cè)和分析，并構(gòu)建了牻牛兒基牻牛兒基焦磷酸合酶蛋白家族的系統(tǒng)進(jìn)化樹。結(jié)果表明，9種唇形科植物的GGPS氨基酸序列理化性質(zhì)基本一致，主要表現(xiàn)為親水性蛋白，有葉綠體轉(zhuǎn)運(yùn)肽，不存在信號(hào)肽，沒有跨膜結(jié)構(gòu)域；大多數(shù)定位于葉綠體，少數(shù)定位于線粒體。蛋白質(zhì)二級(jí)結(jié)構(gòu)中最主要的結(jié)構(gòu)元件是α-螺旋和無規(guī)則卷曲，含有多聚異戊二烯基合成酶的活性結(jié)構(gòu)域、2個(gè)天冬氨酸富集區(qū)結(jié)構(gòu)域、活性位點(diǎn)殘基蓋結(jié)構(gòu)域和鎂離子結(jié)合位點(diǎn)結(jié)構(gòu)域。研究結(jié)果將為GGPS的酶學(xué)特性及二萜類生物合成的分子機(jī)制研究提供理論參考。

[關(guān)鍵詞]唇形科； GGPS合酶； 二萜生物合成； 生物信息學(xué)

[Abstract]Geranylgeranyl pyrophosphate synthase enzyme is one of the key enzymes in the synthesis pathway of diterpenoid. Nine Lamiaceae genus GGPS synthase in Genebank was analyzed in this article. GGPS synthase the nucleic acid sequences and amino acid sequences， physicochemical properties， the signal peptide， leader peptides， transmembrane topological structure， hydrophobic， hydrophilic， subcellular localization， secondary structure， function domain， tertiary structure and evolutional relationship were predicted by using bioinformatics methods.Phylogenetic tree was constructed for the geranylgeranyl pyrophosphate synthase enzyme protein family. The results showed that GGPS amino acid sequence of the physical and chemical properties were basically identical， mainly hydrophilic protein， there existed chloroplast transit peptide， and no signal peptide and membrane structure domain， which mainly located in the chloroplast， the minor part located in mitochondria. The main secondary structures of the proteins are alpha helix and random coil. All these proteins have catalytic residues， aspartate-rich region， active site lid residues， substrate-Mg2+ binding site. The results provide theoretical reference for study on both the enzymatic characteristics of GGPS and the biosynthesis pathway of diterpenoid.

[Key words]Lamiaceae； geranylgeranyl pyrophosphate synthase （GGPS）； diterpene biosynthesis； bioinformatics

二萜類化合物是由4個(gè)異戊二烯結(jié)構(gòu)單元組成的C20類萜類化合物，是以二萜類化合物的共同前體牻牛兒基牻牛兒基焦磷酸 （geranylgeranyl pyrophosphate， GGPP）為底物，在不同的二萜合酶的催化下形成的，廣泛地分布于自然界，在植物、動(dòng)物和海洋生物中均有二萜類化合物的存在。目前研究發(fā)現(xiàn)，二萜類化合物共有119種骨架，1萬多種化合物。主要的骨架類有克羅烷（cleordane）型、貝殼杉烷（kaurane）型、半日花烷（labdane）型和松香烷（abietane）型等。二萜類化合物廣泛存在于植物中，包括唇形科鼠尾草屬植物、巴豆屬植物、香茶菜屬植物、黃芩屬植物、番茄枝屬植物、Eurycoma 屬植物等，而且二萜類化合物具有重要的生理活性，如治療腫瘤的的紫杉醇^[1]；治療心血管系統(tǒng)疾病、抗菌消炎的丹參酮^[2]以及在農(nóng)業(yè)上用于高效調(diào)控植物生長(zhǎng)發(fā)育的赤霉素^[3]等。途徑（pyruvate/ glyceradehyde-3-P pathway），此途徑最初的前提物質(zhì)是丙酮酸和3-磷酸甘油醛（G3P），二者在1-脫氧-D木酮糖-5-磷酸合成酶（DXPS）的催化作用下聚合成1-脫氧-D-木酮糖-5-5磷酸（DXP），DXP經(jīng)異構(gòu)中間體、磷酸化、環(huán)化等步驟逐漸形成2-甲基赤蘚糖醇-4-磷酸（MEP）中間體，再經(jīng)過磷酸化、環(huán)化生成異戊烯基二磷酸（IPP）。IPP在異戊烯基二磷酸異構(gòu)酶（isopentenyl diphosphate isomerase，IDI）的作用下部分轉(zhuǎn)化為雙鍵異構(gòu)體DMAPP。高等植物體內(nèi)，IPP與DMAPP經(jīng)頭尾縮合形成具有C₁₀骨架的牻牛兒基焦磷酸（geranyl pyrophosphate，GPP），GPP加上第2個(gè)IPP單元形成具有C₁₅骨架的法呢基焦磷酸（farnesyl pyrophosphate，F(xiàn)PP），F(xiàn)PP由GGPP合成酶催化再加上第3個(gè)IPP單元縮合反應(yīng)產(chǎn)生具C₂₀骨架的GGPP。由于GGPP是多種初生和次生代謝物等共同前體，因此，GGPS能夠起到調(diào)節(jié)碳流的作用^[4]，從而成為初生和次生代謝途徑中的關(guān)鍵酶之一。

目前，有多種植物的GGPS已經(jīng)或正在被深入研究，包括冷杉Abies grandis、銀杏Ginkgo biloba和紅豆杉Taxus canadensis，T. wallichiana，T. wallichiana var. chinensis，T. x media）等裸子植物，以及以擬南芥Arabidopsis thaliana為代表的多種被子植物。已知以GGPP為前體的具有藥用價(jià)值的化合物有紫杉醇^[5]、銀杏酚和丹參酮等，且已證明GGPP是決定紅豆杉植物合成紫杉醇類物質(zhì)的關(guān)鍵酶之一。對(duì)模式植物擬南芥GGPS的研究表明^[6]，GGPS還參與植物的發(fā)育過程，并具有時(shí)空特異性。該基因在花期表達(dá)，而且其活性只有在花瓣和葉片中檢出，隨著花瓣的發(fā)育，酶活性逐漸降低。另外，在丹參過表達(dá)GGPS株系中隱丹參酮含量的顯著降低，RNAi株系中肉眼可見的顯著表型為根縮短、葉發(fā)黃，這提示赤霉素含量的升高及葉綠素和/或類胡蘿卜素含量的降低^[7]。

唇形科Lamiaceae植物通常含有豐富萜類化合物，廣泛應(yīng)用于醫(yī)藥和食品行業(yè)，如具有利濕、消腫止痛等功效的黃芩^[8-9]，具有強(qiáng)心、抗壓、平喘、抗血栓及降低眼內(nèi)壓等藥理作用的鞘蕊蘇^[10]和用于治療婦科、心血管及腎臟疾病的益母草屬，都具有很高的經(jīng)濟(jì)價(jià)值及實(shí)用價(jià)值。目前唇形科中已從丹參^[11]，毛喉鞘蕊花^[12]和米團(tuán)花^[13]中克隆了GGPS基因。其作為植物二萜類化合物生物合成過程中的關(guān)鍵酶之一，正在受到越來越多的關(guān)注，本文利用生物信息學(xué)方法，對(duì)唇形科植物中已克隆的9個(gè)GGPS合酶的cDNA序列及其編碼的氨基酸序列的組成、理化性質(zhì)、信號(hào)肽、導(dǎo)肽、跨膜結(jié)構(gòu)域、疏水性/親水性、亞細(xì)胞定位、二級(jí)結(jié)構(gòu)、功能結(jié)構(gòu)域、三級(jí)結(jié)構(gòu)和進(jìn)化關(guān)系進(jìn)行預(yù)測(cè)和分析，以期為今后深入研究更多植物GGPP合酶的功能和結(jié)構(gòu)特征以及其生理功能提供一些理論依據(jù)。

1 數(shù)據(jù)下載

以“Lamiaceae”和“geranylgeranyl diphosphate synthase”為搜索關(guān)鍵詞，從美國(guó)國(guó)立生物技術(shù)信息中心（National Center of Biotechnology Information，NCBI）下載來自丹參（Salvia miltiorrhiza，Sm）、毛喉鞘蕊花（Plectranthus barbatus，Pb）、米團(tuán)花（Leucosceptrum canum，Lc）的9種完整的唇形科GGPS核苷酸序列和氨基酸序列（表1）。

2 方法

利用美國(guó)國(guó)家生物技術(shù)信息中心網(wǎng)站（http：//www.ncbi.nlm.nih.gov） 對(duì)唇形科植物GGPS 核苷酸和氨基酸序列進(jìn)行在線分析^[14-15]。核苷酸及其編碼氨基酸序列的組成成分、相對(duì)分子質(zhì)量、等電點(diǎn)等理化性質(zhì)利用ProtParam 在線進(jìn)行分析，開放閱讀框的查找利用ORF Finder在線工具（http：//www.ncbi.nlm.nih.gov/gorf/orfig.cgi）。信號(hào)肽（signal peptide）預(yù)測(cè)利用SignalP 軟件3.0 版^[16]。轉(zhuǎn)運(yùn)肽的預(yù)測(cè)利用Target P1.1Server^[17]（http：//www.cbs.dtu.dk/services/TargetP/）；跨膜結(jié)構(gòu)域用TMHMMServer v. 2.0^[18]（http：//www.cbs.dtu.dk/services/ TMHMM-2.0/）進(jìn)行預(yù)測(cè)；用ProtScale（http：//web.expasy.org/protscale/）對(duì)疏水/親水性進(jìn)行預(yù)測(cè)；PSORT Prediction^[19] （http：//psort.hgc.jp/form.html）進(jìn)行亞細(xì)胞定位分析。用SOPMA（https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=/NPSA/npsa_sopma.html）觀測(cè)其二級(jí)結(jié)構(gòu)，功能域的預(yù)測(cè)用Pfam 27.0（http：//pfam.janelia.org/）和SMART （http：//smart.embl-heidelberg.de/）^[20]進(jìn)行；用SWISS-MODEL（http：//swissmodel.expasy.org）完成GGPS 蛋白高級(jí)結(jié)構(gòu)同源建模；MEGA7.0 軟件構(gòu)建系統(tǒng)進(jìn)化樹。

3 結(jié)果與分析

3.1 唇形科植物牻牛兒基牻牛兒基焦磷酸合酶核苷酸序列的結(jié)構(gòu)及其氨基酸序列的理化性質(zhì)

利用ORF Finder和ProtParam在線工具對(duì)唇形科9種植物GGPS氨基酸序列進(jìn)行理化性質(zhì)分析（表2）。可知其核苷酸序列的起始密碼子均為ATG，終止密碼子均為TGA。氨基酸殘基（amino acids，aa）數(shù)在346～379 aa；各蛋白序列的相對(duì)分子質(zhì)量為37 424.3～41 299.7 kDa，中位值為39 408.66 kDa；理論等電點(diǎn)均在6 PI左右，平均6.33 PI，提示GGPS蛋白為酸性蛋白。從GGPS氨基酸組成中可以看到，9種植物的GGPS蛋白除SmGGPS3外，所含酸性氨基酸殘基比例均高于所含堿性氨基酸殘基比例，進(jìn)一步提示GGPS蛋白為酸性蛋白。各種植物GGPS蛋白中，含量最豐富的氨基酸殘基主要集中在亮氨酸（Leu）、丙氨酸（Ala）、纈氨酸（Val），均不含硒半胱氨酸（Sec）、吡咯賴氨酸（Pyl）。總原子數(shù)，消光系數(shù)，基本一致。除SmGGPS1，PbGGPS，LcGGPS1的不穩(wěn)定系數(shù)小于40，為穩(wěn)定蛋白，其他均為不穩(wěn)定蛋白。

3.2 牻牛兒基牻牛兒基焦磷酸合酶的信號(hào)肽、導(dǎo)肽，跨膜結(jié)構(gòu)域，疏水性/親水性和亞細(xì)胞定位特征

3.2.1 信號(hào)肽、導(dǎo)肽的預(yù)測(cè)和分析 信號(hào)肽（signal peptide）是分泌蛋白和膜蛋白以前體形式合成時(shí)在N端的15～30個(gè)氨基酸序列^[21]。導(dǎo)肽（leader peptide）是一段引導(dǎo)新合成的肽鏈進(jìn)入細(xì)胞器的識(shí)別序列^[22]，導(dǎo)肽的預(yù)測(cè)與分析對(duì)蛋白質(zhì)的功能分析、作用機(jī)制和作用途徑等具有重要意義^[23]。信號(hào)肽屬于導(dǎo)肽的一部分，位于靠近N端的一段氨基酸序列，導(dǎo)肽功能的發(fā)揮需要信號(hào)肽的存在^[24]。利用在線工具SignalP 4.1 Server的神經(jīng)網(wǎng)絡(luò)算法對(duì)9種唇形科植物的GGPS蛋白進(jìn)行信號(hào)肽的預(yù)測(cè)，結(jié)果表明丹參GGPS氨基酸序列中不存在信號(hào)肽，毛喉鞘蕊花和米團(tuán)花GGPS氨基酸序列進(jìn)行信號(hào)肽預(yù)測(cè)也得到相類似的結(jié)果。通過在線預(yù)測(cè)工具TargetP 1.1Server，對(duì)唇形科植物GGPS氨基酸序列進(jìn)行了預(yù)測(cè)。以SmGGPS1為例，預(yù)測(cè)可能性是4，即可能含有低相似度的N端葉綠體轉(zhuǎn)運(yùn)肽（chloroplast transit peptide）。轉(zhuǎn)運(yùn)肽序列長(zhǎng)52個(gè)氨基酸，剪切位點(diǎn)位于Ser52～Ala53。無法確定SmGGPS1是否具有導(dǎo)肽，也未發(fā)現(xiàn)其導(dǎo)肽分裂位點(diǎn)。其他8種唇形科植物的GGPS預(yù)測(cè)分析結(jié)果顯示，SmGGPS2的可靠性為5級(jí)，其余都在4級(jí)以上。LcGGPS4，LcGGPS5具有導(dǎo)肽分裂位點(diǎn)，具有導(dǎo)肽性，且它們的導(dǎo)肽很可能都是葉綠體轉(zhuǎn)運(yùn)肽，提示這些米團(tuán)花中的GGPS蛋白合成后，可能轉(zhuǎn)運(yùn)到葉綠體中發(fā)揮作用。剩下與SmGGPS1相似，都不存在導(dǎo)肽分裂位點(diǎn)，不能確定具有何種導(dǎo)肽。

3.2.2 跨膜結(jié)構(gòu)域特征 跨膜結(jié)構(gòu)域一般由20個(gè)左右的疏水性氨基酸殘基組成，主要形成α-螺旋，常由跨膜蛋白的效應(yīng)區(qū)域所展現(xiàn)。利用在線工具TMHMM Server v.2.0對(duì)SmGGPS1蛋白進(jìn)行跨膜結(jié)構(gòu)進(jìn)行預(yù)測(cè)，分析可知，其整條肽鏈都位于細(xì)胞膜之外，不存在跨膜結(jié)構(gòu)。毛喉鞘蕊花和米團(tuán)花GGPS蛋白跨膜結(jié)構(gòu)域分析結(jié)果與丹參一致，提示本實(shí)驗(yàn)中的GGPS蛋白均不具跨膜結(jié)構(gòu)域。信號(hào)肽是指導(dǎo)靶標(biāo)蛋白質(zhì)跨膜定位到膜上的N端氨基酸序列^[25-26]，所以不含信號(hào)肽，理應(yīng)無跨膜結(jié)構(gòu)域，說明預(yù)測(cè)結(jié)果的合理性。

3.2.3 蛋白疏水性/親水性的預(yù)測(cè) 蛋白質(zhì)親疏水性氨基酸組成是蛋白質(zhì)折疊的主要驅(qū)動(dòng)力^[27]，用Protscale在線工具預(yù)測(cè)親疏水性，結(jié)果表明，SmGGPS1的多肽鏈中第167位氨基酸有最低的親水性分值-2.911。位于260位氨基酸疏水性最強(qiáng)，其分值為2.544。其中，親水性氨基酸占65%，疏水性氨基酸占35%。兩端多親水性氨基酸，中間多疏水性氨基酸，推測(cè)折疊的蛋白為親水性蛋白。其余8種GGPS合酶的疏水性/吸水性都與SmGGPS1類似，這也與跨膜結(jié)構(gòu)預(yù)測(cè)的結(jié)果相吻合。

3.2.4 亞細(xì)胞定位特征 細(xì)胞中蛋白質(zhì)在合成后被轉(zhuǎn)運(yùn)到特定的細(xì)胞器中，蛋白質(zhì)的亞細(xì)胞定位分析及預(yù)測(cè)能極大的加速了蛋白質(zhì)結(jié)構(gòu)和功能的研究^[28]。對(duì)9種唇形科植物的GGPS基因編碼的氨基酸采用PSORT Prediction在線生物學(xué)工具進(jìn)行亞細(xì)胞定位。結(jié)果表明（表3），SmGGPS1，PbGGPS，LcGGPS1，LcGGPS4，LcGGPS5位于膜結(jié)構(gòu)上的可能性大于0.4；SmGGPS2，SmGGPS3位于線粒體基質(zhì)上的可能性大于0.5；LcGGPS2，LcGGPS3位于細(xì)胞質(zhì)的可能性大于0.4。

3.3 二級(jí)結(jié)構(gòu)預(yù)測(cè)

蛋白質(zhì)二級(jí)結(jié)構(gòu)是指蛋白質(zhì)多肽鏈氨基酸殘基借助氫鍵折疊和盤繞形成的α-螺旋、β-折疊、無規(guī)則卷曲以及模體等組件，其中，α-螺旋和β-折疊是最常見的蛋白質(zhì)二級(jí)結(jié)構(gòu)。利用SOPMA對(duì)9種唇形科植物的GGPS合酶序列進(jìn)行二級(jí)結(jié)構(gòu)預(yù)測(cè)（表4），結(jié)果顯示，唇形科GGPS合酶中均有α-螺旋、β-折疊、無規(guī)則卷曲和延伸鏈。以SmGGPS1為例，其主要結(jié)構(gòu)元件是α-螺旋（45.33%）和無規(guī)則卷曲（30.22%），其次是延伸鏈（18.13%）和β-折疊（6.32%）。余下蛋白二級(jí)結(jié)構(gòu)的主要結(jié)構(gòu)元件與SmGGPS1完全一致。

3.4 蛋白質(zhì)功能域的預(yù)測(cè)和分析

功能域（functional domain）又稱結(jié)構(gòu)域，是蛋白質(zhì)分子中介于二級(jí)與三級(jí)結(jié)構(gòu)之間的一種獨(dú)立結(jié)構(gòu)和功能單位，具有特定的生物學(xué)功能^[29-30]。利用Smart在線軟件對(duì)SmGGPS1蛋白的氨基酸序列進(jìn)行功能域的預(yù)測(cè)和分析，結(jié)果表明（圖1），SmGGPS1蛋白含有多聚異戊二烯基合成酶的活性結(jié)構(gòu)域、2個(gè)天冬氨酸富集區(qū)結(jié)構(gòu)域、活性位點(diǎn)殘基蓋結(jié)構(gòu)域和鎂離子結(jié)合位點(diǎn)結(jié)構(gòu)域，其屬于Isoprenoid_Biosyn_C1超家族，為類異戊二烯生物合成酶。對(duì)其他植物也進(jìn)行了功能域的預(yù)測(cè)和分析，除SmGGPS3只有多聚異戊二烯基合成酶的活性結(jié)構(gòu)域和底物結(jié)合位點(diǎn)外，其余唇形科植物的GGPS蛋白均存在上述結(jié)構(gòu)域，這可能是由于SmGGPS3開放閱讀框全長(zhǎng)明顯短于其他植物。

3.5 GGPS蛋白三級(jí)結(jié)構(gòu)的預(yù)測(cè)和分析

蛋白質(zhì)的三級(jí)結(jié)構(gòu)是指蛋白質(zhì)在其二級(jí)結(jié)構(gòu)的基礎(chǔ)上依靠氨基酸側(cè)鏈之間的疏水相互作用、氫鍵、范德華力和靜電作用等進(jìn)一步盤繞、折疊所形成的天然構(gòu)象。蛋白質(zhì)的功能與其三級(jí)結(jié)構(gòu)密切相關(guān)，對(duì)蛋白質(zhì)高級(jí)結(jié)構(gòu)的預(yù)測(cè)和分析，有助于理解蛋白質(zhì)結(jié)構(gòu)與功能之間的相關(guān)性^[31-32]。利用同源建模工具SWISS-MODEL（http：//swissmodel.expasy.org）完成蛋白質(zhì)三級(jí)結(jié)構(gòu)的預(yù)測(cè)和分析工作，找到了模板蛋白（ACCESSION：5E8L_A；Sequence Idenity：76.87%；GMQE：0.73），再用Swiss Pdb-Viewer工具顯示丹參GGPS1結(jié)構(gòu)域的3D結(jié)構(gòu)^[33-35]。結(jié)果顯示：SmGGPS1蛋白由12個(gè)α-螺旋和一些不規(guī)則卷曲組成，與二級(jí)結(jié)構(gòu)的預(yù)測(cè)結(jié)果一致-SmGGPS1主要結(jié)構(gòu)單元為螺旋結(jié)構(gòu)（圖2）。

3.6 GGPS合酶的系統(tǒng)進(jìn)化樹分析

來源于同一祖先的不同植物在進(jìn)化過程中的關(guān)系可以通過進(jìn)化樹來描述，通過構(gòu)建植物進(jìn)化樹，可以了解一種植物在進(jìn)化過程中的地位^[36]。用MEGA 7.0軟件對(duì)唇形科在內(nèi)的25種有代表性的GGPS合酶蛋白序列進(jìn)行系統(tǒng)進(jìn)化樹構(gòu)建（圖3）。結(jié)果顯示，來源于植物，細(xì)菌，真菌，動(dòng)物的GGPS合酶按照不同類群分為4群，在進(jìn)化遺傳學(xué)上親緣越近的物種，在GGPS合酶的分子系統(tǒng)進(jìn)化樹上距離較近，依據(jù)氨基酸序列所得出的系統(tǒng)演化關(guān)系雖不能真實(shí)的反映植物在漫長(zhǎng)歷史長(zhǎng)河中的自然進(jìn)化關(guān)系，其結(jié)果對(duì)判斷不同植物之間的親緣關(guān)系仍具有一定的借鑒意義^[37]。

4 討論

GGPP是合成赤霉素類、二萜類、胡蘿卜素類物質(zhì)的起始前體物，而GGPS則是合成GGPP的關(guān)鍵酶，在植物的次生代謝中，調(diào)控處于代謝分支點(diǎn)前體的代謝方向。一般認(rèn)為二萜類化合物以質(zhì)體來源GGPP為前體，已證明，在擬南芥中，質(zhì)體定位的GGPS蛋白可為赤霉素、類胡蘿卜素、脫落酸和葉綠素等物質(zhì)的合成提供GGPP前體^[38]；在辣椒中，GGPS被證明分別在果實(shí)成熟期和花發(fā)育過程中提供類胡蘿卜素的合成前體^[39-40]，在煙草中，GGPS則為煙草抗蟲（煙草天蛾）物質(zhì) HGL-DTGs（17-hydroxygeranyllinalool diterpenoid glycosides）的合成提供前體^[41]。此外，茉莉酸甲酯（MJ），已經(jīng)驗(yàn)證其在曼地亞紅豆杉、加拿大紅豆杉和番茄等植物中可上調(diào)GGPS合酶的基因表達(dá)，對(duì)二萜類物質(zhì)的產(chǎn)生具有促進(jìn)作用。

蛋白一級(jí)結(jié)構(gòu)預(yù)測(cè)分析結(jié)構(gòu)表明，GGPS蛋白為酸性，親水性，多為不穩(wěn)定蛋白，具有明顯的疏水區(qū)和親水區(qū)，不具有信號(hào)肽，可推測(cè)出GGPS可能不是分泌性蛋白，這與其都沒有跨膜結(jié)構(gòu)相對(duì)應(yīng)。通過導(dǎo)肽分析發(fā)現(xiàn)，GGPS蛋白在細(xì)胞中的分布多樣，說明其在細(xì)胞中廣泛參與生物合成，參與的次生代謝是多樣的。結(jié)合信號(hào)肽預(yù)測(cè)結(jié)果，可推知GGPS蛋白在游離核糖體上合成后，可能通過2種途徑發(fā)揮作用，一是通過導(dǎo)肽進(jìn)入葉綠體發(fā)揮作用；二是不進(jìn)行蛋白轉(zhuǎn)運(yùn)，保留在細(xì)胞質(zhì)基質(zhì)中產(chǎn)生催化作用。GGPS蛋白的二級(jí)結(jié)構(gòu)均以α-螺旋為主要結(jié)構(gòu)，無規(guī)則卷曲、延伸鏈和β-折疊分布于整個(gè)肽鏈結(jié)構(gòu)中。所有唇形科GGPS蛋白氨基酸序列中都含有多聚異戊二烯基合成酶的活性結(jié)構(gòu)域和底物結(jié)合位點(diǎn)結(jié)構(gòu)域，其屬于Isoprenoid_Biosyn_C1超家族，為類異戊二烯生物合成酶。

利用生物信息學(xué)的方法對(duì)唇形科GGPS氨基酸序列的生理生化特性進(jìn)行預(yù)測(cè)和分析，可以為GGPS蛋白及其編碼基因的克隆提供可靠的依據(jù)；對(duì)其序列結(jié)構(gòu)的預(yù)測(cè)和分析，可為其蛋白表達(dá)與修飾提供指導(dǎo)；并為更多物種GGPS蛋白及其編碼基因的克隆提供可靠的依據(jù)。對(duì)其二級(jí)及高級(jí)結(jié)構(gòu)的預(yù)測(cè)和分析有利于深入探討該酶結(jié)構(gòu)和功能之間的關(guān)系、作用機(jī)制和代謝過程。

[參考文獻(xiàn)]

[1]張利平， 程克棣， 朱平. 紫杉烷類化合物的生物轉(zhuǎn)化[J]. 藥學(xué)學(xué)報(bào)， 2004， 39（2）：153.

[2]Du G H， Zhang J T. The general situation and progress of the modern research of red sage root （Radix Salviae miltiorrhizae）[J]. Herald Med， 2004，23（7）：358.

[3]張國(guó)華，張艷潔，叢日晨，等.赤霉素作用機(jī)制研究進(jìn)展[J].西北植物學(xué)報(bào)，2009，29（2）：417.

[4]Laskaris G， Bounkhay M， Theodoridis G， et al. Induction of geranylgeranyl diphosphate synthase activity and taxane accumulation in Taxus baccata， cell cultures after elicitation by methyl jasmonate[J]. Plant Sci， 1999， 147（1）：1.

[5]Laskaris G， Jong C F D， Jaziri M， et al. Geranylgeranyl diphosphate synthase activity and taxane production in Taxus baccata cells[J]. Phytochemistry， 1999， 50（6）：939.

[6]Aharoni A， Giri A P， Deuerlein S， et al. Terpenoid metabolism in wild-type and transgenic Arabidopsis plants[J]. Plant Cell， 2004， 15（12）：2866.

[7]張蕾. 丹參牻牛兒基牻牛兒基焦磷酸含酶基因的克隆與功能研究[D]. 北京：中國(guó)人民解放軍軍事醫(yī)學(xué)科學(xué)院， 2009：3.

[8]何春年， 彭勇， 肖偉，等. 中國(guó)黃芩屬植物傳統(tǒng)藥物學(xué)初步整理[J]. 中國(guó)現(xiàn)代中藥， 2012， 14（1）：16.

[9]劉佳. 鞘蕊蘇有效成分合成酶基因克隆及其含量與生態(tài)環(huán)境相關(guān)性研究[D]. 武漢：湖北中醫(yī)藥大學(xué)， 2015.

[10]龔海群. 中藥金沸草和益母草化學(xué)成分的研究[D]. 蘭州：蘭州大學(xué)， 2011.

[11]Hua W， Song J， Li C， et al. Molecular cloning and characterization of the promoter of SmGGPPs and its expression pattern in Salvia miltiorrhiza[J]. Mol Biol Rep， 2012， 39（5）：5775.

[12]Andersen-Ranberg J， Kongstad K T， Nielsen M T， et al. Expanding the landscape of diterpene structural diversity through stereochemically controlled combinatorial biosynthesis[J]. Angew Chem， 2016， 55（6）：2142.

[13]Liu Y， Luo S H， Schmidt A， et al. A geranylfarnesyl diphosphate synthase provides the precursor for sesterterpenoid （C25） formation in the glandular trichomes of Leucosceptrum canum[J]. Plant Cell， 2016，28（3）： 804.

[14]巴克沃尼什. 生物信息學(xué)基因和蛋白質(zhì)分析的實(shí)用指南[M]. 北京：清華大學(xué)出版社， 2000：16.

[15]Mount D M. Bioinformatics： sequence and genome analysis[M]. Second Edition. New York：Cold Spring Harbor Laboratory Press， 2004：301.

[16]Bendtsen J D， Nielsen H， Von H G， et al. Improved prediction of signal peptides： SignalP 3.0[J]. J Mol Biol， 2004， 340（4）：783.

[17]Ehrbar K， Hapfelmeier S， Stecher B， et al. InvB is required for type Ⅲ-dependent secretion of SopA in Salmonella enterica serovar. typhimurium[J]. J Bacteriol， 2004， 186（4）：1215.

[18]Ikeda M， Arai M， Lao D M， et al. Transmembrane topology prediction methods： a re-assessment and improvement by a consensus method using a dataset of experimentally-characterized transmembrane topologies[J]. Silico Biol， 2002， 2（1）：19.

[19]Emanuelsson O， Nielsen H， Brunak S， et al. Predicting subcellular localization of proteins based on their N-terminal amino acid sequence.[J]. J Mol Biol， 2000， 300（4）：1005.

[20]Letunic I， Doerks T， Bork A P. SMART 7： recent updates to the protein domain annotation resource.[J]. Nucl Acid Res， 2012， 40（Databaseissue）：302.

[21]Von H G. Protein targeting signals[J]. Curr Opin Cell Biol， 1990， 2（4）：604.

[22]董嬌， 周軍， 辛培堯，等. 不同植物L(fēng)DOX/ANS基因的生物信息學(xué)分析[J]. 基因組學(xué)與應(yīng)用生物學(xué)， 2010， 29（5）： 815.

[23]龍芳， 李紹鵬， 李茂富. 7種植物ALAD基因的生物信息學(xué)分析[J]. 基因組學(xué)與應(yīng)用生物學(xué)， 2013（6）：802.

[24]翟中和，王喜忠，丁明孝.細(xì)胞生物學(xué)[M]. 北京：高等教育出版社，2000： 79.

[25]Petersen T N， Brunak S， Von H G， et al. SignalP 4.0： discriminating signal peptides from transmembrane regions[J]. Nat Methods， 2010， 8（10）：785.

[26]Seitz C， Eder C， Deiml B， et al. Cloning， functional identification and sequence analysis of flavonoid 3′-hydroxylase and flavonoid 3′，5′-hydroxylase cDNAs reveals independent evolution of flavonoid 3′，5′-hydroxylase in the Asteraceae family[J]. Plant Mol Biol， 2006， 61（3）：365.

[27]張弢. 不同植物查爾酮合成酶CHS基因的生物信息學(xué)分析[J]. 江西農(nóng)業(yè)學(xué)報(bào)， 2012， 24（6）：5.

[28]Herr D. Secretion of cellulase and beta-glucosidase by Trichoderma viride ITCC-1433 in submerged culture on different substrates[J]. Biotechnol Bioeng， 1979， 21（8）：1361.

[29]薛永常， 聶會(huì)忠， 劉長(zhǎng)斌. 木質(zhì)素合成酶C3H基因的生物信息學(xué)分析[J]. 生物信息學(xué)， 2009， 7（1）：13.

[30]王鏡巖， 朱圣庚， 徐長(zhǎng)法.生物化學(xué)[M]. 北京：高等教育出版社， 2002：533.

[31]蔡娜娜， 陳月輝， 李偉. 基于神經(jīng)網(wǎng)絡(luò)的蛋白質(zhì)三級(jí)結(jié)構(gòu)預(yù)測(cè)[J]. 計(jì)算機(jī)工程， 2010， 36（9）： 176.

[32]陳克克， 武雪. 植物查耳酮異構(gòu)酶生物信息學(xué)分析[J]. 生物信息學(xué)， 2009， 7（3）： 163.

[33]Arnold K， Bordoli L， Kopp J， et al. The SWISS-MODEL workspace： a web-based environment for protein structure homology modelling[J]. Bioinformatics， 2006， 22（2）：195.

[34]Schwede T， Kopp J， Guex N， et al. SWISS-MODEL： an automated protein homology-modeling server[J]. Nucl Acid Res， 2003， 31（13）：3381.

[35]Guex N， Peitsch M C. SWISS-MODEL and the Swiss-PdbViewer： an environment for comparative protein modeling[J]. Electrophoresis， 1997， 18（15）：2714.

[36]李卿， 邸鵬， 陸文銓，等. 丹參柯巴基焦磷酸合酶的生物信息學(xué)分析[J]. 中草藥， 2015， 46（6）：887.

[37]化文平. 丹參GGPP合酶基因的克隆及表達(dá)分析[D]. 西安：陜西師范大學(xué)， 2008：42.

[38]Okada K， Saito T， Nakagawa T， et al. Five geranylgeranyl diphosphate synthases expressed in different organs are localized into three subcellular compartments in Arabidopsis[J]. Plant Physiol， 2000， 122（4）：1045.

[39]Kuntz M， Rmer S， Suire C， et al. Identification of a cDNA for the plastid-located geranylgeranyl pyrophosphate synthase from Capsicum annuum： correlative increase in enzyme activity and transcript level during fruit ripening[J]. Plant J， 1992， 2（1）：25.

[40]Zhu C， Yamamura S， Koiwa H， et al. cDNA cloning and expression of carotenogenic genes during flower development in Gentiana lutea[J]. Plant Mol Biol， 2002， 48（3）：277.

[41]Jassbi A R， Gase K， Hettenhausen C， et al. Silencing geranylgeranyl diphosphate synthase in Nicotiana attenuata dramatically impairs resistance to Tobacco hornworm[J]. Plant Physiol， 2008， 146（3）：974.

[責(zé)任編輯 呂冬梅]