基因SDH2-3的生物信息学分析.docVIP

课程名称：生物信息学 主讲教师：张征锋 学 号： 2011211825 姓 名：肖静 成绩 ： 基因SDH2-3的生物信息学分析 前言 基因SDH2-3位于拟南芥5号染色体上，该基因的表达量很低，只占平均基因的9.5%。SDH2-3和SDH2-1、SDH2-2一起编码拟南芥线粒体呼吸链复合体Ⅱ琥珀酸脱氢酶的铁-硫中心。SDH2-3转录发生于种子成熟期，在干燥阶段持续转录，在干燥种子中富含SDH2-3转录产物，在萌芽期则明显减少[1-2]。对基因SDH2-3进行一系列生物信息学分析，了解该基因及其编码的蛋白质的结构并进行基因序列和蛋白质进化分析。 根据基因号SDH2-3，从NCBI中找到相应的DNA序列及蛋白质序列的fasta格式，之后进行一系列分析操作。 运用GENSCAN软件进行基因结构预测；运用CpGPlot来预测CpG岛；运用POLYAH来进行转录终止信号预测；运用PromoterScan来预测启动子区域。在Genebank中进行blastn同源性搜索，搜索同源性高的前10个基因序列，运用ClustalW2进行同源比对，运用TreeView构建进化树等完成基因序列进化分析。 运用ExPASy数据库中的ProtParam分析工具进行蛋白质的等电点、分子量、氨基酸组成等理化性质分析；运用ExPASy数据库中的ProtScale程序进行蛋白质的亲疏水性分析。运用PREDATOR工具来分析蛋白质的螺旋、折叠、转角、环等二级结构。在蛋白质结构数据库（PDB）中检索同源蛋白质的结构，运用Swiss-Model来预测蛋白质三级结构。运用TMHMM进行蛋白质跨膜区预测；运用SignalP 4.1 Server进行蛋白质信号肽分析。在Genebank中进行blastp同源性搜索，搜索同源性高的前10个蛋白质序列，运用ClustalW2进行同源比对，运用TreeView构建进化树等完成蛋白质序列进化分析。 一 过程及结果分析 1 DNA序列及蛋白质序列 根据基因号，从NCBI中找到相应的DNA序列及蛋白质序列的fasta格式。 网站如下： 2 基因序列分析 2.1 基因结构的预测分析 运用GENSCAN软件进行基因预测。 网站如下：/GENSCAN.html 如图3，在该基因中没有发现外显子。但在该基因中发现了SNG1基因，该基因始于997碱基，终于1194碱基，其后有polyA信号。SNG1是拟南芥54 个SCPL基因中的家族成员之一，属于拟南芥正常生长发育所必需的功能基因。SNG1编码芥子酰基葡萄糖: 苹果酸酰基转移酶。该酶的酰基受体底物是一些变化相对较小的氨基酸残基, 它与羟基肉桂酸的酯代谢有关。SNG1与芥子苹果酸的合成和苯丙次生代谢产物有关[3]。 2.2 CpG岛的预测分析 运用CpGPlot来预测CpG岛。 网站如下：http://www.ebi.ac.uk/Tools/seqstats/emboss_cpgplot/ 如图4，在该基因中没有发现CpG岛。 2.3 转录终止信号的预测分析 运用POLYAH来进行转录终止信号预测。 网站如下： /berry.phtml?topic=polyahgroup=programssubgroup=promoter 如图5，该基因有3个可能的polyA位点，分别位于189、200、929碱基处，权重分别为2.31、7.75、1.29。这与前述GENSCAN的预测结果不一致。可能是因为，真核生物基因组序列本身存在大量的重复序列，当以polyA位点预测基因终止信号位点时会出现较大比例的假阳性。 2.4 启动子区域的预测分析 运用PromoterScan来预测启动子区域。 网站如下：/molbio/proscan/ 如图6，没有发现可能的启动子区域。 2.5 基因序列进化分析 运用BLAST软件包搜索同源性高的前10个基因序列。 网站如下：/Blast.cgi 图7 同源序列 运用ClustalW2进行同源序列比对。 网站如下：http://www.ebi.ac.uk/Tools/msa/clustalw2/ 图8 同源序列比对 运用TreeView软件构建系统进化树。 图9 进化树 3 蛋白质序列分析 3.1 蛋白质一级结构分析 运用ExPASy数据库中的ProtParam分析工具进行蛋白质的等电点、分子量、氨基酸组成等理化性质分析。 网