基因组研究功能基因分析

11梅艳珍南京师范大学生命科学学院基因组学研究——功能基因分析现代生物学实验技术本文档共82页；当前第1页；编辑于星期六\2点32分2要求：掌握常用的序列比对工具能构建进化树能够预测蛋白质的二级结构、疏水区、跨膜区等能够进行简单的同源建模分析了解KEGG数据库的检索本文档共82页；当前第2页；编辑于星期六\2点32分3序列比对——BLAST应用本文档共82页；当前第3页；编辑于星期六\2点32分4

同源性(homology)：指从一些数据中推断出的两个基因或蛋白质序列具有共同祖先的结论，属于质的判断。A和B的关系上，是同源序列，或者非同源序列两种关系。而说A和B的同源性为80％都是不科学的。相似性(similarity)：是指一种直接的数量关系，如部分相同或相似的百分比或其它一些合适的度量。比如说，A序列和B序列的相似性是80％，或者4/5。生物序列的同源性序列间相似性越高，它们是同源序列的可能性就更高本文档共82页；当前第4页；编辑于星期六\2点32分5Blast程序评价序列相似性的两个数据Score：使用打分矩阵对匹配的片段进行打分，这是对各对氨基酸残基（或碱基）打分求和的结果，一般来说，匹配片段越长、相似性越高,则Score值越大。Evalue:在相同长度的情况下，两个氨基酸残基（或碱基）随机排列的序列进行打分，得到上述Score值的概率的大小。E值越小表示随机情况下得到该Score值的可能性越低。我们在获得一个Blast结果时需要看这两个指标。如果Blast获得的目标序列的Score值越高并且E-value越低表明结果越可信，反之越不可信.本文档共82页；当前第5页；编辑于星期六\2点32分6主要的BLAST程序（功能）本文档共82页；当前第6页；编辑于星期六\2点32分71.登陆blast主页组装的基因组序列库基本blast特定的BLAST所有的BLAST基因数据库本文档共82页；当前第7页；编辑于星期六\2点32分88核酸数据库中比对核酸序列蛋白质数据库中比对蛋白质序列蛋白质数据库中比对核酸序列蛋白质数据库中比对核酸序列核酸数据库中比对蛋白质序列本文档共82页；当前第8页；编辑于星期六\2点32分9Fasta格式文件本文档共82页；当前第9页；编辑于星期六\2点32分1010什么是fasta格式？怎么建立？新建一个txt文本文件，命名如:bph.txtFasta的格式：>序列名称序列本文档共82页；当前第10页；编辑于星期六\2点32分11本文档共82页；当前第11页；编辑于星期六\2点32分12121.序列信息部分填入查询（query）的序列序列范围（默认全部）选择搜索数据库如果接受其他参数默认设置，点击开始搜索本文档共82页；当前第12页；编辑于星期六\2点32分1313去冗余GenBank编码序列PDB+SwissProt+PIR+PRF本文档共82页；当前第13页；编辑于星期六\2点32分14

常用的检索数据库

14本文档共82页；当前第14页；编辑于星期六\2点32分15以下列蛋白序列为例，进行BLAST搜索：>P1MSDNGPQSNQRSAPRITFGGPTDSTDNNQNGGRNGARPKQRRPQGLPNNTASWFTALTQHGKEELRFPRGQGVPINTNSGPDDQIGYYRRATRRVRGGDGKMKELSPRWYFYYLGTGPEASLPYGANKEGIVWVATEGALNTPKDHIGTRNPNNNAATVLQLPQGTTLPKGFYAEGSRGGSQASSRSSSRSRGNSRNSTPGSSRGNSPARMASGGGETALALLLLDRLNQLESKVSGKGQQQQGQTVTKKSAAEASKKPRQKRTATKQYNVTQAFGRRGPEQTQGNFGDQDLIRQGTDYKHWPQIAQFAPSASAFFGMSRIGMEVTPSGTWLTYHGAIKLDDKDPQFKDNVILLNKHIDAYKTFPPTEPKKDKKKKTDEAQPLPQRQKKQPTVTLLPAADMDDFSRQLQNSMSGASADSTQA本文档共82页；当前第15页；编辑于星期六\2点32分16本文档共82页；当前第16页；编辑于星期六\2点32分17本文档共82页；当前第17页；编辑于星期六\2点32分1818基因名来源物种一致程度，登录号本文档共82页；当前第18页；编辑于星期六\2点32分19所选序列下载序列本文档共82页；当前第19页；编辑于星期六\2点32分20Cluster比对本文档共82页；当前第20页；编辑于星期六\2点32分21Clustalx的工作界面

（多序列比对模式）本文档共82页；当前第21页；编辑于星期六\2点32分22Clustal的工作原理Clustal输入多个序列快速的序列两两比对，计算序列间的距离，获得一个距离矩阵。邻接法(NJ)构建一个树根据进化树，渐进比对多个序列。本文档共82页；当前第22页；编辑于星期六\2点32分23Clustalx的输出结果.aln格式文件这个文件是默认输出，可以转换成各种格式，而且很多软件都支持这种格式。.dnd格式文件引导树。就是根据两两序列相似值构建的一个指导后面多重联配的启发树不能做进化分析。进化分析要考虑的所有同源位点的一个综合效应，因此应该用.aln格式文件专门做进化分析。本文档共82页；当前第23页；编辑于星期六\2点32分24多序列比对实例输入文件的格式(fasta)：>KCC2_YEASTNYIFGRTLGAGSFGVVRQARKLSTN……>DMK_HUMANDFEILKVIGRGAFSEVAVVKMKQTGQVYAMKIMNK…….>KPRO_MAIZETRKFKVELGRGESGTVYKGVLEDDRHVAVKKLEN……>DAF1_CAEELQIRLTGRVGSGRFGNVSRGDYRGEAVAVKVFNALD……>1CSNHYKVGRRIGEGSFGVIFEGTNLLNN……不留空格本文档共82页；当前第24页；编辑于星期六\2点32分25第一步：输入序列文件。本文档共82页；当前第25页；编辑于星期六\2点32分26本文档共82页；当前第26页；编辑于星期六\2点32分27本文档共82页；当前第27页；编辑于星期六\2点32分28建议用treeview打开outtree，然后可以编辑本文档共82页；当前第28页；编辑于星期六\2点32分29本文档共82页；当前第29页；编辑于星期六\2点32分30建树软件-mega本文档共82页；当前第30页；编辑于星期六\2点32分31MEGA5可以识别fasta格式文件将重命名为17-RNASE1.fasta建树软件-mega本文档共82页；当前第31页；编辑于星期六\2点32分32本文档共82页；当前第32页；编辑于星期六\2点32分33ClustalW参数设置本文档共82页；当前第33页；编辑于星期六\2点32分34多序列联配后结果本文档共82页；当前第34页；编辑于星期六\2点32分35以.meg格式保存结果本文档共82页；当前第35页；编辑于星期六\2点32分36回到MEGA主窗口打开所保存的文件（.meg）本文档共82页；当前第36页；编辑于星期六\2点32分37点击按钮打开文件窗口本文档共82页；当前第37页；编辑于星期六\2点32分38显示保守位点显示变异位点本文档共82页；当前第38页；编辑于星期六\2点32分39回到MEGA主窗口构建进化树当前打开的文件选择邻接法建树本文档共82页；当前第39页；编辑于星期六\2点32分40选择Bootstrap检验本文档共82页；当前第40页；编辑于星期六\2点32分41本文档共82页；当前第41页；编辑于星期六\2点32分42蛋白质二级结构预测本文档共82页；当前第42页；编辑于星期六\2点32分43蛋白质结构为什么如此重要的？氨基酸序列只有折叠成特定的空间结构才具有相应的活性和相应的生物学功能DNA序列蛋白质序列蛋白质结构转录&翻译折叠本文档共82页；当前第43页；编辑于星期六\2点32分44为什么要研究蛋白质结构?生物体中许多重要的功能由蛋白质完成分析蛋白质结构、功能及其关系是蛋白质组计划中的一个重要组成部分分析蛋白质结构有助于药物设计研究有助于了解蛋白质相互作用，这对于生物学、医学和药学都是非常重要本文档共82页；当前第44页；编辑于星期六\2点32分45蛋白质二级结构α-helix(30-35%)

α-螺旋β-sheet/β-strand(20-25%)

β-折叠Coil(40-50%)无规则卷曲Loop环β-turnβ-转角本文档共82页；当前第45页；编辑于星期六\2点32分4646蛋白质3D结构α螺旋Β折叠环或转角转角或卷曲本文档共82页；当前第46页；编辑于星期六\2点32分47本文档共82页；当前第47页；编辑于星期六\2点32分48JPred预测结果α螺旋β折叠本文档共82页；当前第48页；编辑于星期六\2点32分49二级结构预测网站PHDJPREDpbio.dundee.ac.uk/~www-jpred/PSIPREDNNPREDICT/~nomi/nnpredict.htmlChouandFassman本文档共82页；当前第49页；编辑于星期六\2点32分50预测蛋白质的理化性质本文档共82页；当前第50页；编辑于星期六\2点32分51部分预测工具ComputepI/Mw(ExPASy)计算蛋白序列的等电点和分子量TGREASE计算蛋白质序列疏水性工具TMHMM蛋白质跨膜区预测More…本文档共82页；当前第51页；编辑于星期六\2点32分52等电点，分子量预测工具本文档共82页；当前第52页；编辑于星期六\2点32分53本文档共82页；当前第53页；编辑于星期六\2点32分54本文档共82页；当前第54页；编辑于星期六\2点32分55TGREASE疏水性参数高正值的氨基酸具有更大的疏水性而低负值的氨基酸具有更强的亲水性本文档共82页；当前第55页；编辑于星期六\2点32分56本文档共82页；当前第56页；编辑于星期六\2点32分57蛋白质跨膜区预测(TMHMM)本文档共82页；当前第57页；编辑于星期六\2点32分58本文档共82页；当前第58页；编辑于星期六\2点32分59本文档共82页；当前第59页；编辑于星期六\2点32分60信号肽分析本文档共82页；当前第60页；编辑于星期六\2点32分61SignalP软件2.0版（）对信号肽分析。本文档共82页；当前第61页；编辑于星期六\2点32分62信号肽的作用一般是帮助蛋白质穿膜用的,跟蛋白质的细胞定位有关系。本文档共82页；当前第62页；编辑于星期六\2点32分63同源建模蛋白质三维结构预测本文档共82页；当前第63页；编辑于星期六\2点32分643D预测是可能的，因为：序列信息决定三级结构序列相似性(>30%)倾向于结构相似性3D预测是必须的，因为：DNA序列蛋白质序列空间结构本文档共82页；当前第64页；编辑于星期六\2点32分65本文档共82页；当前第65页；编辑于星期六\2点32分66本文档共82页；当前第66页；编辑于星期六\2点32分67蛋白质结构预测方法：同源建模法(Comparativehomologymodeling)依据蛋白序列与已经结构蛋白比对信息构建3D模型折叠识别法(Threadingfoldrecognition)寻找与未知蛋白最合适的模板，进行序列与结构比对，最终建立结构模型从头预测法(Abinitio/denovomethods)根据序列本身来从头预测蛋白质结构本文档共82页；当前第67页；编辑于星期六\2点32分68同源建模基本原理：1、一个蛋白质的结构由其氨基酸序列唯一的决定。由一级结构，在理论上，足以获取其二级、三级结构。2、三级结构的保守型远远大于一级结构的保守型。应用限制：模板蛋白和目标蛋白的序列一致性需要大于30%本文档共82页；当前第68页；编辑于星期六\2点32分69SWISS-MODELSWISS-MODEL:网址非专业人士应用最为广泛的一个在线建模服务器。特点：简单、自动化、对学术团队免费。Automatedmode:自动模式，可以称为是最傻瓜的方式提交自己的氨基酸序列+邮箱即可适用：一致性较高时本文档共82页；当前第69页；编辑于星期六\2点32分70本文档共82页；当前第70页；编辑于星期六\2点32分71邮箱模型命名氨基酸序列本文档共82页；当前第71页；编辑于星期六\2点32分72本文档共82页；当前第72页；编辑于星期六\2点32分73KEGG数据库本文档共82页；当前第73页；编辑于星期六\2点32分74本文档共82页；当前第74页；编辑于星期六\2点32分75特点KEGG是一个整合了基因组、化学和系统功能信息的数据库。把从已经完整测序的基因组中得到的基因目录与更高级别的细胞、物种和生态系统水平的系统功能关联起来是KEGG数据库的特色之一。人工创建了一个知识库，这个知识库是基于使用一种可计算的形式捕捉和组织实验得到的知识而形成的系统功能知识库。它是一个生物系统的计算机模拟。与