一种基于树搜索和片段组装的蛋白质结构预测方法，包括以下步骤：A1、获取蛋白质的pdb格式的文件并清洗出所需数据；A2、生成片段库；A3、选取力场模型；A4、采用Rosetta的score3能量函数；A5、将整个能量图景离散化，并将每层进一步离散化，分成单独的分块区，每次搜索根据能量的权重随机选取一个能量层，并在该能量层依概率选取一个分块区，如果分块区中包含有构象，采用片段组装方法，随机选取序列上的一个片段，然后在片段库中随机选取一个片段，替换序列上的目标片段，用蒙特卡洛准则来判断是否接受这个构象，如果接受则将它放入一个集合当中。本发明专利技术可以大大减少计算量和计算时间，保证搜索到能量较低的构象。
【技术实现步骤摘要】
一种基于树搜索和片段组装的蛋白质结构预测方法
本专利技术涉及计算机应用领域，生物信息学，最优化理论，分子生物学等，特别涉及一种蛋白质三维结构预测方法，属于将现代智能优化方法应用到蛋白质三维结构预测。
技术介绍
生物信息学通过综合利用生物学，计算机科学和信息技术而揭示大量而复杂的生物数据所赋有的生物学奥秘。是当前研究的一个热点。生物信息学研究成果已经被广泛应用于序列比对，蛋白质比对，基因识别分析，分子进化，序列重叠群装配，遗传密码，药物设计，生物系统，蛋白质结构预测等。其中蛋白质结构预测是生物信息学领域的一个重要的分支。著名的诺贝尔化学奖获得者Anfinsen通过对核糖核酸酶A的经典研究表明去折叠的蛋白质在体外可以自发的进行再折叠，仅仅是序列本身已经包括了蛋白质正确折叠的所有信息，并提出蛋白质折叠的热力学假说，这一理论认为结构同源的蛋白质可以通过不同的折叠途径形成相似的天然构象，蛋白质的一级结构决定了其三维结构，即蛋白质的氨基酸序列决定了蛋白质的三维空间结构，从热力学分析，蛋白质的天然构型对应着其物理能量最小的状态。随着2000年基因组工作草图的完成，基因组的工作重心已经从测序转向了基因功能的识别，这一问题的解决显得尤为重要。人们通过研究蛋白质的内在结构来了解其内在机理，并找出蛋白质折叠的真正驱动力，这样可以按照人们的设想设计出符合特定需求的非天然蛋白质，对疾病进行有效的预测和控制。蛋白质结构从头预测也因此具有重大的意义。然后蛋白质的结构测定远远赶不上基因组测序的速度，单纯的依靠实验的方法测定所有的蛋白质的三维结构已经远远不能满足后基因时代对于蛋白质数据分析的需求，有必要发展一种可靠的理论预测方法，从而迅速、简便的获取蛋白质的结构信息。蛋白质的从头预测方法就是选择一种力场模型，将蛋白质结构抽象成一个数学模型，使其具有可计算性，然后进行算法的设计，从算法的时间和空间复杂度出发，利用全局以及局部的优化算法在能量曲面上搜索全局最优解，在数学上属于一种极其复杂的非凸多极值优化问题。利用蛋白质的一级序列，利用热力学方法，计算得出蛋白质的三维结构。针对力场模型的优化问题，2005年Bradley在结合蒙特卡洛优化方法开发了分子建模Rosetta程序，成功预测CASP6中T0281目标蛋白(包含70个残基)的Cα-RMSD为1.6对某些小目标蛋白(<85个残基)，甚至得到误差精度小于1.5的高分辨率三维解析结构；2006年，Zhan采用Basin-Paving(BP)方法得到了脑啡肽在ECEPP/2、ECEPP/3两种力场模型下最小能量构象；2009年Roosi提出几种BP改进策略，并给出了比较结果。结合随机全局优化算法，如遗传算法(GeneticAlgorithms,GA)、差分进化(DifferentialEvolution,DE)算法也应用于力场模型的优化，但是由于蛋白质是一种强柔性的大分子体系，其势函数表达式极其复杂，存在极多局部极小点，这些算法极容易陷入某个局优解；同时由于随机算法缺乏全局收敛理论依据及解的不确定性，使用了全局选择因子，算法只能收敛到全局最优解，而忽略了众多局部极值解，进一步限制了它们在实际问题中的应用。我们了解到蛋白质迅速折叠是因为随机热运动使得构象发生改变，而根据漏斗状的能量图景中所体现的原理，构象变化会使蛋白质沿着能量下降的方向搜索到天然结构。然而，如果采用穷举搜索的方式来遍历整个能量图景，那是非常不科学的方式，首先要想遍历整个能量图景，其计算量非常庞大，这个值与离散化能量图景后的最小单位值相关，离散得越精细，计算量相应越大，相应的时间成本也越高，需要的计算机资源也越多。
技术实现思路
为了克服现有的蛋白质结构预测方式的计算量较大、计算时间很长的不足，本专利技术提供一种基于树搜索和片段组装的蛋白质结构预测方法，采用将能量图景适度离散化并依概率随机选择能量层和分块区的方式，可以大大减少计算量，计算时间，同时保证搜索到能量较低的构象，达到预测蛋白质结构的目的。本专利技术解决其技术问题所采用的技术方案是：一种基于树搜索和片段组装的蛋白质结构预测方法，所述预测方法包括以下步骤：A1、获取蛋白质的pdb格式的文件并清洗出所需数据；A2、生成片段库；A3、选取力场模型，力场模型表示形式如下：Eprotein＝WinterrepEinterrep+WinteratrEinteratr+WsolvationEsolvation+Wbb/schbEbb/schb+Wbb/bbhbEbb/bbhb+Wsc/schbEsc/schb+WpairEpair+WdunbrackEdunbrack+WramaErama+WreferenceEreference（1）式中，Eprotein表示蛋白质的总能量，Einterrep表示范德华排斥力作用，Winterrep为Einterrep在整体的权重，Einteratr表示范德华吸引力作用，Winteratr为Einteratr在整体的权重，Esolvation为Lazarids和Karplus描述的隐含的溶解作用，Wsolvation为Esolvation在整体的权重，Ebb/schb、Ebb/bbhb、Esc/schb为依赖方向的氢键能量，Wbb/schb、Wbb/bbhb、Wsc/schb分别为其能量在整体的权重，Epair为残基对静电作用，Wpair为Epair在整体的权重，Edunbrack为氨基酸基于旋转异构体库的内部的能量，Wdunbrack为Edunbrack在整体的权重，Erama为参考特定位置的Ramachandrin骨链扭力，Wrama为Erama在整体的权重，Ereference为未折叠态的蛋白质的参考能量，Wreference为Ereference在整体的权重，Rosetta的能量函数就是将所有的能量项通过各自的权重线性相加；用于计算能量层的权重W(l)的公式：W(l)为第L层的权重值，Eprotein是蛋白质构象的能量值，该能量层被选中的概率P(l)为：某一个分块区被选中的的概率P(cell)为：其中nsel记录的是该分块区被选中的次数，nconfs记录的是绑定到该分块区上的构象的数目；A4、选取能量函数，采用Rosetta的score3能量函数；A5、蛋白质结构预测过程如下：Step1：初始化工作，导入各种程序运行所需的包和模块，氨基酸序列的初始化，创建相应的文件夹存放低能量的构象和相关数据和图表，创建MC对象和Mover对象等；Step2：进入迭代循环，给定迭代次数，由迭代次数控制循环次数；Step3：根据权重函数选择相应的能量层，权重函数列表由公式（2）动态生成并扩展，Eprotein由score3函数计算得出并放入一个能量列表中；Step4：依概率随机选中特定的分块区，概率由公式（4）计算得出；Step5：首先判断该分块区中是否存在构象，如果存在则取出该构象；Step6：本文档来自技高网...
【技术保护点】
一种基于树搜索和片段组装的蛋白质结构预测方法，其特征在于：所述预测方法包括以下步骤：A1、获取蛋白质的pdb格式的文件并清洗出所需数据；A2、生成片段库；A3、选取力场模型，力场模型表示形式如下：Eprotein＝Winter repEinter rep+Winter atrEinter atr+WsolvationEsolvation+Wbb/sc hbEbb/sc hb+Wbb/bb hbEbb/bb hb+Wsc/sc hbEsc/sc hb+WpairEpair+WdunbrackEdunbrack+WramaErama+WreferenceEreference （1）式中，Eprotein表示蛋白质的总能量，Einter rep表示范德华排斥力作用，Winter rep为Einter rep在整体的权重，Einter atr表示范德华吸引力作用，Winter atr为Einter atr在整体的权重，Esolvation为Lazarids和Karplus描述的隐含的溶解作用，Wsolvation为Esolvation在整体的权重，Ebb/sc hb、Ebb/bb hb、Esc/sc hb为依赖方向的氢键能量，Wbb/sc hb、Wbb/bb hb、Wsc/sc hb分别为其能量在整体的权重，Epair为残基对静电作用，Wpair为Epair在整体的权重，Edunbrack为氨基酸基于旋转异构体库的内部的能量，Wdunbrack为Edunbrack在整体的权重，Erama为参考特定位置的Ramachandrin骨链扭力，Wrama为Erama在整体的权重，Ereference为未折叠态的蛋白质的参考能量，Wreference为Ereference在整体的权重，Rosetta的能量函数就是将所有的能量项通过各自的权重线性相加；用于计算能量层的权重W(l)的公式：W(l)=Eprotein2+1.0222---(2)]]>W(l)为第L层的权重值，Eprotein是蛋白质构象的能量值，该能量层被选中的概率P(l)为：P(l)=W(l)Σl=1100W(l)---(3)]]>某一个分块区被选中的的概率P(cell)为：P(cell)=1.0[(1.0+nsel)*nconfs]---(4)]]>其中nsel记录的是该分块区被选中的次数，nconfs记录的是绑定到该分块区上的构象的数目；A4、选取能量函数，采用Rosetta的score3能量函数；A5、蛋白质结构预测过程如下：Step1：初始化工作，导入各种程序运行所需的包和模块，氨基酸序列的初始化，创建相应的文件夹存放低能量的构象和相关数据和图表，创建MC对象和Mover对象等；Step2：进入迭代循环，给定迭代次数，由迭代次数控制循环次数；Step3：根据权重函数选择相应的能量层，权重函数列表由公式（2）动态生成并扩展，Eprotein由score3函数计算得出并放入一个能量列表中；Step4：依概率随机选中特定的分块区，概率由公式（4）计算得出；Step5：首先判断该分块区中是否存在构象，如果存在则取出该构象；Step6：对取出的蛋白质构象进行扰动，也就是片段组装，随机选取该构象中的一段长度为3的片段，然后随机选取片段库中与其序列相同的片段，其本质就是替换该片段的3个二面角的值；Step7：用蒙特卡洛评判标准判断构象是否接受；如果接受则放入集合，不接受则跳转到Step2，进入下一次迭代，直到循环次数达到设定的循环值则停止。...
【技术特征摘要】
1.一种基于树搜索和片段组装的蛋白质结构预测方法，其特征在于：所述预测方法包括以下步骤：A1、获取蛋白质的pdb格式的文件并清洗出所需数据；A2、生成片段库；A3、选取力场模型，力场模型表示形式如下：Eprotein＝WinterrepEinterrep+WinteratrEinteratr+WsolvationEsolvation+Wbb/schbEbb/schb+Wbb/bbhbEbb/bbhb+Wsc/schbEsc/schb+WpairEpair+WdunbrackEdunbrack+WramaErama+WreferenceEreference(1)式中，Eprotein表示蛋白质的总能量，Einterrep表示范德华排斥力作用，Winterrep为Einterrep在整体的权重，Einteratr表示范德华吸引力作用，Winteratr为Einteratr在整体的权重，Esolvation为Lazarids和Karplus描述的隐含的溶解作用，Wsolvation为Esolvation在整体的权重，Ebb/schb、Ebb/bbhb、Esc/schb为依赖方向的氢键能量，Wbb/schb、Wbb/bbhb、Wsc/schb分别为其能量在整体的权重，Epair为残基对静电作用，Wpair为Epair在整体的权重，Edunbrack为氨基酸基于旋转异构体库的内部的能量，Wdunbrack为Edunbrack在整体的权重，Erama为参考特定位置的Ramachandrin骨链扭力，Wrama为Erama在整体的权重，Ereference为未折叠态的蛋白质的参考能量，Wreference为Ereference在整体的权重，Rosetta的能量函数就是将所有的能量项通过各自的权重线性相加；用于计算能量层的权重W(l)的公式：
【专利技术属性】
技术研发人员：张贵军，陈铭，秦传庆，郝小虎，周晓根，梅珊，李章维，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江;33
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。