- 1、本文档共36页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
生物信息学bioinf4
* * 四、核苷酸和蛋白质序列为基础的数据库检索 Sequence-based database searching 序列对位排列(sequence alignment) 将两条或多条序列对位排列,突出相似的结构区域 序列1 序列2 两条DNA序列对位排列分析 两条蛋白质序列对位排列分析 分析功能 分析物种进化 检测突变、插入或缺失 序列延长 序列定位 基因表达谱分析 用途 序列对位排列排列分析的种类 序列对库对位排列分析 从数据库中寻找同源序列 主要涉及核苷酸数据库和蛋白质数据库 两序列对位排列分析 多序列对位排列分析 (一)序列对位排列分析的基本原理 1、记分矩阵(scoring matrix) 记分矩阵中含有两条序列对位排列时具体使用的分值 分数越高,两条序列匹配越好 DNA序列对位记分矩阵 序列1 A C G T T A G C 序列2 A C T T T G G C 记分 0.9 0.9 -0.1 0.9 0.9 -0.1 0.9 0.9 = 5.2 蛋白质序列对位排列分析记分复杂 一致氨基酸的记分不同 稀有氨基酸(C),分值高 普通氨基酸(S),分值低 相似氨基酸也记分,如D-E 序列1:TTYGAPPWCS 序列2: TGYAPPPWS * *** * 序列1:TTYGAPPWCS 序列2:TGYAPPPWS * * *** 多种记分矩阵 80年代建立的PAM矩阵(如PAM30、PAM70) 近年建立的BLOSUM矩阵(如BLOSUM62、BLOSUM80、BLOSUM45) 基于更敏感的对位排列分析 2、空位(间隔)罚分(gap penalty) 基因进化过程中产生突变 序列对位排列分析时允许插入空位 空位罚分涉及两个参数 插入 缺失 空位开放(gap opening) 空位延伸(gap extension) 蛋白质序列对位记分 序列1 V D S - C Y 序列2 V E S L C Y 记分 4 2 4 -11 9 7 BLAST FASTA Blitz (二)序列对库对位排列分析 主要检索体系 用待分析序列对数据库进行相似性分析 重复许多次的两两序列对位排列分析 从数据库中找出所有同源序列 1、基本概念 (1)Sequence identity 和 sequence similarity Identity: 两条序列在同一位点上的核苷酸或 氨基酸残基完全相同 Similarity (positive): 两条序列在同一位点上的 氨基酸残基的化学性质相似 Query: 1 IGQAQCSTFRGRIYNETNIDSAFATQRQANCP 32 IGQAQC TF+ RIYNET +AFAT +ANCP Sbjet: 2 IGQAQCGTFKDRIYNET---TAFATSLRANCP 32 (2)Global alignment 和 local alignment Query Subject Query Subject Query Subject Global alignment: 两条完整的序列相比较 Local alignment: 两条序列中相似程度最高的部分 相比较 (3)Gapped alignment 和 ungapped alignment Query Subject Query Subject Query Subject Query Subject Gapped alignment: 为达到最佳 alignment,序列中加入空位 Ungapped alignment:相比较序列的核苷酸或氨基酸 序列连续 (4)Alignment score 和 E (expect) value 衡量两条相比较序列相似程度的标准 raw score: 原始分,分值越大,两个比较序列相似程度越大 bit score: 采用统计学方法以原始分为基础计算的 E value: 期望得到的、完全由机会造成的、相当于或大于目 前分值的alignment 次数 E = 10, 5e - 46 = 5 ? 10-46,E
您可能关注的文档
- 牛津英语7A Unit 2-1-3 (Welcome to the unit).ppt
- 牛津英语7A Unit6 Welcome.ppt
- 牛皮癣(神经性皮炎).ppt
- 牛肉番茄面5029.ppt
- 牛津高中英语M2 Unit 3 Amazing people课件3 Word power(共36张PPT) .ppt.ppt
- 牛津高中英语Module5 Unit2 Task公开课.ppt
- 牛顿定律的初值敏感性(PPT 13) 【精编】.ppt
- 牛牛津小学英语3A Unit2 Nice to meet you第四课时课件.ppt
- 牡丹 [课件].ppt
- 物件导向技术—UML理论与实作(PPT-40).ppt
文档评论(0)