《生物信息学研究方法概述.ppt》由会员分享,可在线阅读,更多相关《生物信息学研究方法概述.ppt(35页珍藏版)》请在金锄头文库上搜索。
1、第二章第二章 生物信息学研究方法概述生物信息学研究方法概述马飞马飞南京师范大学南京师范大学比较基因组学与生物信息学实验室比较基因组学与生物信息学实验室一、生物信息学研究方法概述一、生物信息学研究方法概述1 生物信息学研究的三个层面生物信息学研究的三个层面初级层面初级层面中级层面中级层面高级层面高级层面初级层面初级层面 基于现有的生物信息数据库和资源,利用成熟的生物信息学工具基于现有的生物信息数据库和资源,利用成熟的生物信息学工具(专业网站、软件)解决生物信息学问题(专业网站、软件)解决生物信息学问题生物信息数据库(生物信息数据库(NCBI、EBI等)等)基因组序列分析、序列比对软件(基因组序列
2、分析、序列比对软件(GCG、BLAST、CLUSTAL等)等)系统发育树构造软件(系统发育树构造软件(PHYLIP、PALM、MEGA等)等)分子动力学模拟软件(分子动力学模拟软件(GROMACS、NAMD等)等)搜集、整理有特色的生物信息学数据集搜集、整理有特色的生物信息学数据集中级层面中级层面 利用数值计算方法、数理统计方法和相关的工具,研究生物利用数值计算方法、数理统计方法和相关的工具,研究生物信息学问题信息学问题概率、数理统计基础概率、数理统计基础科学计算基础科学计算基础现有的数理统计和科学计算工具(现有的数理统计和科学计算工具(EXCEL、SPSS、SAS、 MATLAB等)等)建立
3、有特色的生物信息学数据库建立有特色的生物信息学数据库高级层面高级层面 提出有重要意义的生物信息学问题;自主创新,发展新型方法,提出有重要意义的生物信息学问题;自主创新,发展新型方法,开发新型工具,引领生物信息学领域研究方向。开发新型工具,引领生物信息学领域研究方向。面向生物学领域,解决生物学问题面向生物学领域,解决生物学问题数学、物理、化学、计算科学等思想和方法数学、物理、化学、计算科学等思想和方法建立模型,发展算法建立模型,发展算法自行编程,开发软件,建立网页(自行编程,开发软件,建立网页(Linux系统、系统、C/C+、PERL、 数据库技术)数据库技术)http:/rich.yunda.
4、org/test/amphioxusest/ http:/rich.yunda.org/test/rg01/index.php 从事生物信息学研究应具备多方面的科学基础从事生物信息学研究应具备多方面的科学基础从事生物信息学研究应具备多方面的科学基础从事生物信息学研究应具备多方面的科学基础(1)、一定的计算能力,包括相应的软、硬设备。要有各种数据库或者能)、一定的计算能力,包括相应的软、硬设备。要有各种数据库或者能与国际、国内的数据库系统进行有效的交流。要有发达、稳定的互联网络与国际、国内的数据库系统进行有效的交流。要有发达、稳定的互联网络系统;系统;(2)、强有力的创新算法和软件。没有算法创新
5、,生物信息学就无法获得)、强有力的创新算法和软件。没有算法创新,生物信息学就无法获得持续的发展;持续的发展;(3)、与实验科学,特别是与自动化的大规模高通量的生物学研究方法与)、与实验科学,特别是与自动化的大规模高通量的生物学研究方法与平台技术建立广泛、紧密的联系。这些技术,既是产生生物信息数据的主平台技术建立广泛、紧密的联系。这些技术,既是产生生物信息数据的主要方法,又是验证生物信息学研究结果的关键手段。要方法,又是验证生物信息学研究结果的关键手段。 从事生物信息学研究的人员必须具备多学科交叉的知识。从事生物信息学研究的人员必须具备多学科交叉的知识。2 生物信息学的生物信息学的“降龙十八掌降
6、龙十八掌”第一式第一式第一式第一式 见龙在田见龙在田见龙在田见龙在田 (1 1)掌握生物信息数据库及其查掌握生物信息数据库及其查掌握生物信息数据库及其查掌握生物信息数据库及其查询搜索方法询搜索方法询搜索方法询搜索方法(Database & searching)(Database & searching)对分子生物信息数据库的种类以及某些具体数据库的掌握对分子生物信息数据库的种类以及某些具体数据库的掌握和了解和了解从现有数据库中熟练获得需要的数据信息(尤其是二级数从现有数据库中熟练获得需要的数据信息(尤其是二级数据库)据库)能熟练地进行数据库查询和数据库搜索(数据库查询系统能熟练地进行数据库查询
7、和数据库搜索(数据库查询系统Entrez、SRS;搜索工具;搜索工具BLAST等)等)数据库技术、互联网技术数据库技术、互联网技术第二式第二式第二式第二式 飞龙在天飞龙在天飞龙在天飞龙在天(2 2)学会生物信息学软件和工学会生物信息学软件和工学会生物信息学软件和工学会生物信息学软件和工具的应用具的应用具的应用具的应用(Software & application)(Software & application)利用成熟的生物信息学工具(专业网站、软件)解决生物信息学问题利用成熟的生物信息学工具(专业网站、软件)解决生物信息学问题基因组序列分析、序列比对软件(基因组序列分析、序列比对软件(GCG
8、、BLAST、CLUSTAL等)等)系统发育树构造软件(系统发育树构造软件(PHYLIP、PALM等等)基因芯片检测分析软件(商业软件基因芯片检测分析软件(商业软件ScanArray、Array-Pro等等 )分子动力学模拟软件(分子动力学模拟软件(GROMACS、NAMD等等)第三式第三式第三式第三式 鸿渐于陆鸿渐于陆鸿渐于陆鸿渐于陆(3 3)掌握概率论基础掌握概率论基础掌握概率论基础掌握概率论基础(Probability theory)(Probability theory)随机事件、概率随机事件、概率随机变量、概率分布随机变量、概率分布大数定律、中心极限定理大数定律、中心极限定理几乎用于
9、生物信息学的各个方面几乎用于生物信息学的各个方面“Most of the problems in computational sequence analysis are essentially statistical.”“Biological sequence analysis”第四式第四式第四式第四式 或跃在渊或跃在渊或跃在渊或跃在渊(4 4)掌握数理统计基础掌握数理统计基础掌握数理统计基础掌握数理统计基础(Statistical methods)(Statistical methods)样本和统计量(方差、均值样本和统计量(方差、均值)参数估计、假设检验参数估计、假设检验基本的统计分析(方
10、差分析、协方差分析、回归分析)基本的统计分析(方差分析、协方差分析、回归分析)常用统计软件的运用(常用统计软件的运用(SPSS、SAS)几乎用于生物信息学的各个方面几乎用于生物信息学的各个方面第五式第五式第五式第五式 羚羊触藩羚羊触藩羚羊触藩羚羊触藩 (5 5)熟悉基于频率的组分分析熟悉基于频率的组分分析熟悉基于频率的组分分析熟悉基于频率的组分分析方法和权重矩阵方法方法和权重矩阵方法方法和权重矩阵方法方法和权重矩阵方法(Composition analysis & (Composition analysis & weight matrix method)weight matrix method
11、)符号(如碱基)频率反映具有生物学意义的序列特征,符号(如碱基)频率反映具有生物学意义的序列特征,如内含子剪接位点的发现,如内含子剪接位点的发现,KOZAK规则的发现等规则的发现等核酸组分、氨基酸组分、密码子使用频率核酸组分、氨基酸组分、密码子使用频率k-tuples/k-mers频率分析频率分析权重矩阵权重矩阵(weight matrix)(weight matrix) 分析分析主要用于具有特定生物学意义的序列特征的分析主要用于具有特定生物学意义的序列特征的分析权重矩阵分析方法举例权重矩阵分析方法举例权重矩阵分析方法举例权重矩阵分析方法举例例:人类基因例:人类基因例:人类基因例:人类基因内含
12、子内含子内含子内含子/ /外显外显外显外显子剪接位点的子剪接位点的子剪接位点的子剪接位点的序列特征分析序列特征分析序列特征分析序列特征分析R = A or GR = A or GY = C or UY = C or UN = A, G, C or UN = A, G, C or U供体位点供体位点供体位点供体位点受体位点受体位点受体位点受体位点BayesianBayesian打分函数用于剪接位点预测的公式打分函数用于剪接位点预测的公式打分函数用于剪接位点预测的公式打分函数用于剪接位点预测的公式The likelihood that a property value v (of a new st
13、ructure) is drawn from the splicing site is:Score for the overall likelihood of the query sequence being a site is:Say we have a sequence S= S1S2Sn. Then one need to calculateSo to look for a donor site in the sequence, we might calculate 第六式第六式第六式第六式 潜龙勿用潜龙勿用潜龙勿用潜龙勿用 (6 6)信息论方法信息论方法信息论方法信息论方法(Infor
14、mation method)(Information method)信息熵信息熵信息的度量:是信息的度量:是信息符号信息符号出现何种出现何种状态状态的一种不确定性程度,信息的的一种不确定性程度,信息的获得要对不确定性进行否定。获得要对不确定性进行否定。生物信息的符号如生物信息的符号如ACGT四种符号,状态空间即其所有可能的排列四种符号,状态空间即其所有可能的排列用于结构预测用于结构预测信息熵信息熵H刻画了由刻画了由pi表示的随机试验结果的先验不确定性,或观察表示的随机试验结果的先验不确定性,或观察到输出时所获得的信息量。到输出时所获得的信息量。第七式第七式第七式第七式 利涉大川利涉大川利涉大川
15、利涉大川 (7 7)期望最大化(期望最大化(期望最大化(期望最大化(EMEM)方法)方法)方法)方法(Expectation Maximization)(Expectation Maximization)适用于具有隐变量的模型和问题,如神经网络模型中的隐节点和适用于具有隐变量的模型和问题,如神经网络模型中的隐节点和HMM模型中的隐状态等模型中的隐状态等用于结构的识别,用于结构的识别,如如Motif识别的识别的MEME方法、方法、HMM中的中的Baum-Welch算法算法第八式第八式第八式第八式 神龙摆尾神龙摆尾神龙摆尾神龙摆尾 (8 8)动态规划方法动态规划方法动态规划方法动态规划方法(Dyn
16、amic Programming)(Dynamic Programming)第九式第九式第九式第九式 密云不雨密云不雨密云不雨密云不雨(9 9)迭代方法迭代方法迭代方法迭代方法(Iteration)(Iteration)迭代的目的通常是在状态空间找到目标函数收敛的稳定解迭代的目的通常是在状态空间找到目标函数收敛的稳定解迭代的目的通常是在状态空间找到目标函数收敛的稳定解迭代的目的通常是在状态空间找到目标函数收敛的稳定解在运用模式识别方法时,对系统参数的学习通常要经过迭代来实现在运用模式识别方法时,对系统参数的学习通常要经过迭代来实现在运用模式识别方法时,对系统参数的学习通常要经过迭代来实现在运用
17、模式识别方法时,对系统参数的学习通常要经过迭代来实现迭代必须能够不断逼近稳定解迭代必须能够不断逼近稳定解迭代必须能够不断逼近稳定解迭代必须能够不断逼近稳定解第十式第十式第十式第十式 突如其来突如其来突如其来突如其来(1010)回归、拟合、相关性分析、回归、拟合、相关性分析、回归、拟合、相关性分析、回归、拟合、相关性分析、关联分析关联分析关联分析关联分析(Regression, fitting, (Regression, fitting, correlation & association)correlation & association)经典的统计分析方法经典的统计分析方法主要目的:描述和预
18、测自变量与因变量间的关系主要目的:描述和预测自变量与因变量间的关系第十一式第十一式第十一式第十一式 双龙取水双龙取水双龙取水双龙取水 (1111)判别分析方法判别分析方法判别分析方法判别分析方法( (DiscriminantDiscriminant analysis) analysis)第十二式第十二式第十二式第十二式 鱼跃于渊鱼跃于渊鱼跃于渊鱼跃于渊(1212)聚类分析方法聚类分析方法聚类分析方法聚类分析方法(Clustering method)(Clustering method)聚类分析(群分析)是实用多元统计分析的一个新分支,正处于发展聚类分析(群分析)是实用多元统计分析的一个新分支,
19、正处于发展阶段。理论上尚未完善,但应用十分广泛。实质上是一种分类问题,目的阶段。理论上尚未完善,但应用十分广泛。实质上是一种分类问题,目的是建立一种分类方法,将一批数据按照特征的亲疏、相似程度进行分类。是建立一种分类方法,将一批数据按照特征的亲疏、相似程度进行分类。条件:研究对象总体的类别数目未知,也不知总体样本的具体分类情条件:研究对象总体的类别数目未知,也不知总体样本的具体分类情况况目的:通过分析,选定描述个体相似程度的统计量、确定总体分类数目的:通过分析,选定描述个体相似程度的统计量、确定总体分类数目、建立分类方法;对研究对象给出合理的分类。目、建立分类方法;对研究对象给出合理的分类。(
20、“物以类聚物以类聚”是聚类是聚类分析的基本出发点分析的基本出发点 )定性、经验的分类的局限定性、经验的分类的局限定性、经验的分类的局限定性、经验的分类的局限 分类较粗、数据量小、凭借经验分类较粗、数据量小、凭借经验分类较粗、数据量小、凭借经验分类较粗、数据量小、凭借经验谱系聚类法(系统聚类法)、动态聚类法、模糊聚类法谱系聚类法(系统聚类法)、动态聚类法、模糊聚类法谱系聚类法(系统聚类法)、动态聚类法、模糊聚类法谱系聚类法(系统聚类法)、动态聚类法、模糊聚类法生物信息学中的聚类分析问题:生物信息学中的聚类分析问题:生物信息学中的聚类分析问题:生物信息学中的聚类分析问题: 根据根据DNA芯片获得的
21、基因表达数据进行基因聚类(数据量庞大)芯片获得的基因表达数据进行基因聚类(数据量庞大) 蛋白质相互作用网络的分类蛋白质相互作用网络的分类 根据不同物种的大分子序列进行相似性比较并构建系统发育树根据不同物种的大分子序列进行相似性比较并构建系统发育树X(1)X(2)X(3)X(4)X(5)GibbonSymphalangusHumanGorillaChimpanzee黑猩猩黑猩猩猩猩猩猩猿猴猿猴长臂猿长臂猿第十三式第十三式第十三式第十三式 震惊百里震惊百里震惊百里震惊百里(1313)MarkovMarkov模型的应用模型的应用模型的应用模型的应用(Markov model)(Markov mode
22、l)MarkovMarkov过程:从一种状态转移到另一种状态时,过程仅取决于前面过程:从一种状态转移到另一种状态时,过程仅取决于前面过程:从一种状态转移到另一种状态时,过程仅取决于前面过程:从一种状态转移到另一种状态时,过程仅取决于前面n n种状态,是一种有序种状态,是一种有序种状态,是一种有序种状态,是一种有序n n模型。模型。模型。模型。 n n是影响下一个状态选择的状态数。是影响下一个状态选择的状态数。是影响下一个状态选择的状态数。是影响下一个状态选择的状态数。最简单的最简单的最简单的最简单的MarkovMarkov过程是一阶过程,状态的选择完全取决于前一状态,过程是一阶过程,状态的选择
23、完全取决于前一状态,过程是一阶过程,状态的选择完全取决于前一状态,过程是一阶过程,状态的选择完全取决于前一状态,这种选择是依照概率来选择的。这种选择是依照概率来选择的。这种选择是依照概率来选择的。这种选择是依照概率来选择的。状态的选择是概率的,而非确定的。故状态的选择是概率的,而非确定的。故状态的选择是概率的,而非确定的。故状态的选择是概率的,而非确定的。故MarkovMarkov过程本质上是一种随机过程本质上是一种随机过程本质上是一种随机过程本质上是一种随机过程。过程。过程。过程。第十四式第十四式第十四式第十四式 损则有孚损则有孚损则有孚损则有孚(1414)隐隐隐隐MarkovMarkov模
24、型方法模型方法模型方法模型方法(HMM method)(HMM method)将核苷酸序列看成一个随将核苷酸序列看成一个随机序列,机序列,DNA序列的编码部分序列的编码部分与非编码部分在核苷酸的选用与非编码部分在核苷酸的选用频率上对应着不同的频率上对应着不同的Markov模模型。由于这些型。由于这些MarkovC+C+GG+ +C CGGC CGGC C GG C+C+C C GG+ +GG B BE E0.130.130.120.120.0340.0340.0100.0100.0120.0120.0030.0030.00320.00320.00020.0002模型的统计规律是未知的,模型的统
25、计规律是未知的,而而HMM能够自动寻找出它能够自动寻找出它们隐藏的统计规律。对于们隐藏的统计规律。对于高等生物这样复杂的高等生物这样复杂的DNA序列,序列,HMM必须学习不同必须学习不同的基因结构的信号。的基因结构的信号。隐隐Markov模型模型(HMM)语音识别语音识别语音识别语音识别(Speech recognition)(Speech recognition)光字符识别光字符识别光字符识别光字符识别(Optical character (Optical character recognition)recognition)生物序列分析生物序列分析生物序列分析生物序列分析(Biologica
26、l sequence (Biological sequence analysis)analysis)(1)序列比较与搜寻(尤其是多序列比对)序列比较与搜寻(尤其是多序列比对)(2)基因及信号的识别、预测(包括)基因及信号的识别、预测(包括DNA编码与非编码区的识别、真核编码与非编码区的识别、真核基因剪接位点信号识别、非编码区的转录调控信号识别、信号肽识别基因剪接位点信号识别、非编码区的转录调控信号识别、信号肽识别)(3)蛋白质二级结构、家族、超家族预测、分类等)蛋白质二级结构、家族、超家族预测、分类等生物特征识别生物特征识别生物特征识别生物特征识别(Biometrics)(Biometrics
27、)第十五式第十五式第十五式第十五式 时乘六龙时乘六龙时乘六龙时乘六龙(1515)感知器与人工神经网络方法感知器与人工神经网络方法感知器与人工神经网络方法感知器与人工神经网络方法( (PerceptronPerceptron & ANN & ANN method)method)第十六式第十六式第十六式第十六式 龙战于野龙战于野龙战于野龙战于野(1616)决策树、支持向量机及其决策树、支持向量机及其决策树、支持向量机及其决策树、支持向量机及其它模式识别方法它模式识别方法它模式识别方法它模式识别方法(Decision tree & SVM (Decision tree & SVM method)me
28、thod)模式识别是在输入样本中寻找特征并识别对象的一种方法。模式识别是在输入样本中寻找特征并识别对象的一种方法。模式识别主要有两种方法,一种是根据统计特征进行识别,另一种是模式识别主要有两种方法,一种是根据统计特征进行识别,另一种是根据对象的结构特征进行识别,而后者常用的方法为句法识别。根据对象的结构特征进行识别,而后者常用的方法为句法识别。在基因识别中,对于在基因识别中,对于DNA序列上的功能位点和特征信号的识别都需要序列上的功能位点和特征信号的识别都需要用到模式识别。用到模式识别。第十七式第十七式第十七式第十七式 履霜冰至履霜冰至履霜冰至履霜冰至(1717)微分方程的数值方法微分方程的数
29、值方法微分方程的数值方法微分方程的数值方法(Numerical methods)(Numerical methods)分子动力学模拟:研究生物大分子的构象,主要还是用基于半经验势分子动力学模拟:研究生物大分子的构象,主要还是用基于半经验势分子动力学模拟:研究生物大分子的构象,主要还是用基于半经验势分子动力学模拟:研究生物大分子的构象,主要还是用基于半经验势函数的分子动力学方法,而量子力学则在确定势函数的参数和研究局部性函数的分子动力学方法,而量子力学则在确定势函数的参数和研究局部性函数的分子动力学方法,而量子力学则在确定势函数的参数和研究局部性函数的分子动力学方法,而量子力学则在确定势函数的参
30、数和研究局部性质时起作用。对蛋白质进行动力学研究是利用计算机进行模拟实验的基础。质时起作用。对蛋白质进行动力学研究是利用计算机进行模拟实验的基础。质时起作用。对蛋白质进行动力学研究是利用计算机进行模拟实验的基础。质时起作用。对蛋白质进行动力学研究是利用计算机进行模拟实验的基础。分子动力学得到一组动力学微分方程,要求得到初值问题的解。分子动力学得到一组动力学微分方程,要求得到初值问题的解。分子动力学得到一组动力学微分方程,要求得到初值问题的解。分子动力学得到一组动力学微分方程,要求得到初值问题的解。微分方程的数值求解:有限差分法、有限元法微分方程的数值求解:有限差分法、有限元法微分方程的数值求解:有限差分法、有限元法微分方程的数值求解:有限差分法、有限元法十七式合一十七式合一十七式合一十七式合一 亢龙有悔亢龙有悔亢龙有悔亢龙有悔 (1818)最终要诀:各类方法综合最终要诀:各类方法综合最终要诀:各类方法综合最终要诀:各类方法综合运用运用运用运用All in one!All in one!综合运用不同的研究方法综合运用不同的研究方法始终面向生物学问题始终面向生物学问题知识和技能的学习方法知识和技能的学习方法文献的查阅和阅读方法文献的查阅和阅读方法中、英文论文的写作方法中、英文论文的写作方法