唇形科植物脚6基脚6基焦磷酸合酶编码基因及其氨基酸序列的生物信息学分析
开篇：润墨网以专业的文秘视角，为您筛选了一篇唇形科植物脚6基脚6基焦磷酸合酶编码基因及其氨基酸序列的生物信息学分析范文，如需获取更多写作素材，在线客服老师一对一协助。欢迎您的阅读与分享！
[摘要]脚6基脚6基焦磷酸合酶（geranylgeranyl pyrophosphate synthase，G
>> 唇形科植物脚6基脚6基焦磷酸合酶编码基因及其氨基酸序列的生物信息学分析 雷公藤脚6基焦磷酸合酶基因TwGPPS克隆与表达分析 FZ6基因及其蛋白的生物信息学分析 睡前泡脚的好处多 6个事项要注意 6招助你开发脚部健康潜能 冬日里，6种方法暖手脚 一天从脚“High”6次？这是病！ 杜仲法尼烯基焦磷酸合酶基因cDNA全长的克隆与序列分析 脚 牯牛降国家地质公园唇形科药用植物资源及其利用 腿与脚的外开在芭蕾基训中的作用 分级债基6天暴涨34% 债市“小牛”持续 32只债基回报超6% 债基赚钱效应发散 做好6件事冬季手脚不再冰凉 基TPICl6F690单片机的温湿度测量装置设计 脚的保健 我的“再生脚” 宝贝你的脚 妈妈的脚 月亮的脚在哪 常见问题解答 当前所在位置：l）进行亚细胞定位分析。用SOPMA（）观测其二级结构，功能域的预测用Pfam 27.0（http：///）和SMART （http：//smart.embl-heidelberg.de/）[20]进行；用SWISS-MODEL（http：//）完成GGPS 蛋白高级结构同源建模；MEGA7.0 软件构建系统进化树。
3 结果与分析
3.1 唇形科植物脚6基脚6基焦磷酸合酶核苷酸序列的结构及其氨基酸序列的理化性质
利用ORF Finder和ProtParam在线工具对唇形科9种植物GGPS氨基酸序列进行理化性质分析（表2）。可知其核苷酸序列的起始密码子均为ATG，终止密码子均为TGA。氨基酸残基（amino acids，aa）数在346～379 aa；各蛋白序列的相对分子质量为37 424.3～41 299.7 kDa，中位值为39 408.66 kDa；理论等电点均在6 PI左右，平均6.33 PI，提示GGPS蛋白为酸性蛋白。从GGPS氨基酸组成中可以看到，9种植物的GGPS蛋白除SmGGPS3外，所含酸性氨基酸残基比例均高于所含碱性氨基酸残基比例，进一步提示GGPS蛋白为酸性蛋白。各种植物GGPS蛋白中，含量最丰富的氨基酸残基主要集中在亮氨酸（Leu）、丙氨酸（Ala）、缬氨酸（Val），均不含硒半胱氨酸（Sec）、吡咯赖氨酸（Pyl）。总原子数，消光系数，基本一致。除SmGGPS1，PbGGPS，LcGGPS1的不稳定系数小于40，为稳定蛋白，其他均为不稳定蛋白。
3.2 脚6基脚6基焦磷酸合酶的信号肽、导肽，跨膜结构域，疏水性/亲水性和亚细胞定位特征
3.2.1 信号肽、导肽的预测和分析 信号肽（signal peptide）是分泌蛋白和膜蛋白以前体形式合成时在N端的15～30个氨基酸序列[21]。导肽（leader peptide）是一段引导新合成的肽链进入细胞器的识别序列[22]，导肽的预测与分析对蛋白质的功能分析、作用机制和作用途径等具有重要意义[23]。信号肽属于导肽的一部分，位于靠近N端的一段氨基酸序列，导肽功能的发挥需要信号肽的存在[24]。利用在线工具SignalP 4.1 Server的神经网络算法对9种唇形科植物的GGPS蛋白进行信号肽的预测，结果表明丹参GGPS氨基酸序列中不存在信号肽，毛喉鞘蕊花和米团花GGPS氨基酸序列进行信号肽预测也得到相类似的结果。通过在线预测工具TargetP 1.1Server，对唇形科植物GGPS氨基酸序列进行了预测。以SmGGPS1为例，预测可能性是4，即可能含有低相似度的N端叶绿体转运肽（chloroplast transit peptide）。转运肽序列长52个氨基酸，剪切位点位于Ser52～Ala53。无法确定SmGGPS1是否具有导肽，也未发现其导肽分裂位点。其他8种唇形科植物的GGPS预测分析结果显示，SmGGPS2的可靠性为5级，其余都在4级以上。LcGGPS4，LcGGPS5具有导肽分裂位点，具有导肽性，且它们的导肽很可能都是叶绿体转运肽，提示这些米团花中的GGPS蛋白合成后，可能转运到叶绿体中发挥作用。剩下与SmGGPS1相似，都不存在导肽分裂位点，不能确定具有何种导肽。
3.2.2 跨膜结构域特征 跨膜结构域一般由20个左右的疏水性氨基酸残基组成，主要形成α-螺旋，常由跨膜蛋白的效应区域所展F。利用在线工具TMHMM Server v.2.0对SmGGPS1蛋白进行跨膜结构进行预测，分析可知，其整条肽链都位于细胞膜之外，不存在跨膜结构。毛喉鞘蕊花和米团花GGPS蛋白跨膜结构域分析结果与丹参一致，提示本实验中的GGPS蛋白均不具跨膜结构域。信号肽是指导靶标蛋白质跨膜定位到膜上的N端氨基酸序列[25-26]，所以不含信号肽，理应无跨膜结构域，说明预测结果的合理性。
3.2.3 蛋白疏水性/亲水性的预测 蛋白质亲疏水性氨基酸组成是蛋白质折叠的主要驱动力[27]，用Protscale在线工具预测亲疏水性，结果表明，SmGGPS1的多肽链中第167位氨基酸有最低的亲水性分值-2.911。位于260位氨基酸疏水性最强，其分值为2.544。其中，亲水性氨基酸占65%，疏水性氨基酸占35%。两端多亲水性氨基酸，中间多疏水性氨基酸，推测折叠的蛋白为亲水性蛋白。其余8种GGPS合酶的疏水性/吸水性都与SmGGPS1类似，这也与跨膜结构预测的结果相吻合。
3.2.4 亚细胞定位特征 细胞中蛋白质在合成后被转运到特定的细胞器中，蛋白质的亚细胞定位分析及预测能极大的加速了蛋白质结构和功能的研究[28]。对9种唇形科植物的GGPS基因编码的氨基酸采用PSORT Prediction在线生物学工具进行亚细胞定位。结果表明（表3），SmGGPS1，PbGGPS，LcGGPS1，LcGGPS4，LcGGPS5位于膜结构上的可能性大于0.4；SmGGPS2，SmGGPS3位于线粒体基质上的可能性大于0.5；LcGGPS2，LcGGPS3位于细胞质的可能性大于0.4。
3.3 二级结构预测
蛋白质二级结构是指蛋白质多肽链氨基酸残基借助氢键折叠和盘绕形成的α-螺旋、β-折叠、无规则卷曲以及模体等组件，其中，α-螺旋和β-折叠是最常见的蛋白质二级结构。利用SOPMA对9种唇形科植物的GGPS合酶序列进行二级结构预测（表4），结果显示，唇形科GGPS合酶中均有α-螺旋、β-折叠、无规则卷曲和延伸链。以SmGGPS1为例，其主要结构元件是α-螺旋（45.33%）和无规则卷曲（30.22%），其次是延伸链（18.13%）和β-折叠（6.32%）。余下蛋白二级结构的主要结构元件与SmGGPS1完全一致。
3.4 蛋白质功能域的预测和分析
功能域（functional domain）又称结构域，是蛋白质分子中介于二级与三级结构之间的一种独立结构和功能单位，具有特定的生物学功能[29-30]。利用Smart在线软件对SmGGPS1蛋白的氨基酸序列进行功能域的预测和分析，结果表明（图1），SmGGPS1蛋白含有多聚异戊二烯基合成酶的活性结构域、2个天冬氨酸富集区结构域、活性位点残基盖结构域和镁离子结合位点结构域，其属于Isoprenoid_Biosyn_C1超家族，为类异戊二烯生物合成酶。对其他植物也进行了功能域的预测和分析，除SmGGPS3只有多聚异戊二烯基合成酶的活性结构域和底物结合位点外，其余唇形科植物的GGPS蛋白均存在上述结构域，这可能是由于SmGGPS3开放阅读框全长明显短于其他植物。
3.5 GGPS蛋白三级结构的预测和分析
蛋白质的三级结构是指蛋白质在其二级结构的基础上依靠氨基酸侧链之间的疏水相互作用、氢键、范德华力和静电作用等进一步盘绕、折叠所形成的天然构象。蛋白质的功能与其三级结构密切相关，对蛋白质高级结构的预测和分析，有助于理解蛋白质结构与功能之间的相关性[31-32]。利用同源建模工具SWISS-MODEL（http：//）完成蛋白质三级结构的预测和分析工作，找到了模板蛋白（ACCESSION：5E8L_A；Sequence Idenity：76.87%；GMQE：0.73），再用Swiss Pdb-Viewer工具显示丹参GGPS1结构域的3D结构[33-35]。结果显示：SmGGPS1蛋白由12个α-螺旋和一些不规则卷曲组成，与二级结构的预测结果一致-SmGGPS1主要结构单元为螺旋结构（图2）。
3.6 GGPS合酶的系统进化树分析
来源于同一祖先的不同植物在进化过程中的关系可以通过进化树来描述，通过构建植物进化树，可以了解一种植物在进化过程中的地位[36]。用MEGA 7.0软件对唇形科在内的25种有代表性的GGPS合酶蛋白序列进行系统进化树构建（图3）。结果显示，来源于植物，细菌，真菌，动物的GGPS合酶按照不同类群分为4群，在进化遗传学上亲缘越近的物种，在GGPS合酶的分子系统进化树上距离较近，依据氨基酸序列所得出的系统演化关系虽不能真实的反映植物在漫长历史长河中的自然进化关系，其结果对判断不同植物之间的亲缘关系仍具有一定的借鉴意义[37]。
4 讨论
GGPP是合成赤霉素类、二萜类、胡萝卜素类物质的起始前体物，而GGPS则是合成GGPP的关键酶，在植物的次生代谢中，调控处于代谢分支点前体的代谢方向。一般认为二萜类化合物以质体来源GGPP为前体，已证明，在拟南芥中，质体定位的GGPS蛋白可为赤霉素、类胡萝卜素、脱落酸和叶绿素等物质的合成提供GGPP前体[38]；在辣椒中，GGPS被证明分别在果实成熟期和花发育过程中提供类胡萝卜素的合成前体[39-40]，在烟草中，GGPS则为烟草抗虫（烟草天蛾）物质 HGL-DTGs（17-hydroxygeranyllinalool diterpenoid glycosides）的合成提供前体[41]。此外，茉莉酸甲酯（MJ），已经验证其在曼地亚红豆杉、加拿大红豆杉和番茄等植物中可上调GGPS合酶的基因表达，对二萜类物质的产生具有促进作用。
蛋白一级结构预测分析结构表明，GGPS蛋白为酸性，亲水性，多为不稳定蛋白，具有明显的疏水区和亲水区，不具有信号肽，可推测出GGPS可能不是分泌性蛋白，这与其都没有跨膜结构相对应。通过导肽分析发现，GGPS蛋白在细胞中的分布多样，说明其在细胞中广泛参与生物合成，参与的次生代谢是多样的。结合信号肽预测结果，可推知GGPS蛋白在游离核糖体上合成后，可能通过2种途径发挥作用，一是通过导肽进入叶绿体发挥作用；二是不进行蛋白转运，保留在细胞质基质中产生催化作用。GGPS蛋白的二级结构均以α-螺旋为主要结构，无规则卷曲、延伸链和β-折叠分布于整个肽链结构中。所有唇形科GGPS蛋白氨基酸序列中都含有多聚异戊二烯基合成酶的活性结构域和底物结合位点结构域，其属于Isoprenoid_Biosyn_C1超家族，为类异戊二烯生物合成酶。
利用生物信息学的方法对唇形科GGPS氨基酸序列的生理生化特性进行预测和分析，可以为GGPS蛋白及其编码基因的克隆提供可靠的依据；对其序列结构的预测和分析，可为其蛋白表达与修饰提供指导；并为更多物种GGPS蛋白及其编码基因的克隆提供可靠的依据。对其二级及高级结构的A测和分析有利于深入探讨该酶结构和功能之间的关系、作用机制和代谢过程。
[参考文献]
[1]张利平， 程克棣， 朱平. 紫杉烷类化合物的生物转化[J]. 药学学报， 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]张国华，张艳洁，丛日晨，等.赤霉素作用机制研究进展[J].西北植物学报，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]张蕾. 丹参脚6基脚6基焦磷酸含酶基因的克隆与功能研究[D]. 北京：中国人民军事医学科学院， 2009：3.
[8]何春年， 彭勇， 肖伟，等. 中国黄芩属植物传统药物学初步整理[J]. 中国现代中药， 2012， 14（1）：16.
[9]刘佳. 鞘蕊苏有效成分合成酶基因克隆及其含量与生态环境相关性研究[D]. 武汉：湖北中医药大学， 2015.
[10]龚海群. 中药金沸草和益母草化学成分的研究[D]. 兰州：兰州大学， 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]巴克沃尼什. 生物信息学基因和蛋白质分析的实用指南[M]. 北京：清华大学出版社， 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]董桑 周军， 辛培尧，等. 不同植物LDOX/ANS基因的生物信息学分析[J]. 基因组学与应用生物学， 2010， 29（5）： 815.
[23]龙芳， 李绍鹏， 李茂富. 7种植物ALAD基因的生物信息学分析[J]. 基因组学与应用生物学， 2013（6）：802.
[24]翟中和，王喜忠，丁明孝.细胞生物学[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基因的生物信息学分析[J]. 江西农业学报， 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]薛永常， 聂会忠， 刘长斌. 木质素合成酶C3H基因的生物信息学分析[J]. 生物信息学， 2009， 7（1）：13.
[30]王镜岩， 朱圣庚， 徐长法.生物化学[M]. 北京：高等教育出版社， 2002：533.
[31]蔡娜娜， 陈月辉， 李伟. 基于神经网络的蛋白质三级结构预测[J]. 计算机工程， 2010， 36（9）： 176.
[32]陈克克， 武雪. 植物查耳酮异构酶生物信息学分析[J]. 生物信息学， 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]李卿， 邸鹏， 陆文铨，等. 丹参柯巴基焦磷酸合酶的生物信息学分析[J]. 中草药， 2015， 46（6）：887.
[37]化文平. 丹参GGPP合酶基因的克隆及表_分析[D]. 西安：陕西师范大学， 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.