[发明专利]一种基于Loop区域高斯扰动的群体蛋白质结构预测方法

说明书

一种基于Loop区域高斯扰动的群体蛋白质结构预测方法，首先，通过种群初始化得到全局搜索的构象；然后通过交叉操作，交叉种群中个体Loop区域的部分；其次，对交叉后的种群和初始种群合并到一起挑选出能量较小的个体，使用高斯分布的二面角来对构象的Loop区域扰动，并使用玻尔兹曼概率判断是否接收；最后使用聚类算法对输出的最后一代种群聚类，得到最终的蛋白质三级结构。本发明提供一种预测精度较高的基于Loop区域高斯扰动的群体蛋白质结构预测方法。

技术领域

本发明涉及一种生物信息学、智能优化、计算机应用领域、尤其涉及的是，一种基于Loop区域高斯扰动的群体蛋白质结构预测方法。

背景技术

蛋白质是生物的物质基础，是有机大分子，是构成细胞的基本有机物，是生命活动的主要承担者。没有蛋白质就没有生命，氨基酸是蛋白质的基本组成单位，它是与生命及与各种形式的生命活动紧密联系在一起的物质。因为蛋白质特定的三维结构决定了蛋白质的功能，所以了解蛋白质的结构对于认识蛋白质的功能有十分重要的意义。了解蛋白质结构和功能之间的关系可以帮助我们设计具有特定功能的新型蛋白质、设计药物和疫苗解决人类的许多疾病，通过了解蛋白质折叠类疾病的病理，从而找出解决这类疾病的治疗方法。目前，实验室测定蛋白质结构的方法主要有X-射线衍射法和核磁共振法(NMR)，但是这两种方法所需时间长且耗资巨大，通常需要半年至一年时间，耗资几十万美元，而且有局限性，不能大规模地测定蛋白质三级结构，而使用计算机模拟的蛋白质结构预测是现在最有前景的蛋白质结构预测方法。

使用计算机模拟的蛋白质结构预测方法主要有同源建模法和从头预测法。同源建模法的出发点是认为序列相似的蛋白质结构也相似。基于这种认识，同源建模法计算目标蛋白质序列与模板之间的相似性，并以序列联配的形式呈现序列相似性计算结果，最后从联配出发构建目标蛋白质的空间结构。同源建模法适用于在模板库中存在同源序列的目标蛋白质，当目标蛋白质与结构模板之间的相似度30％时，同源建模法一般可以预测到较高精度的蛋白质三级结构。当序列相似度较低时，蛋白质结构预测精度较差。从头建模法也叫无模板建模方法，不依赖于已知结构数据库，因此即使在模板数据库中不存在与目标蛋白相似的结构时，从头预测法也能对目标蛋白质预测。

目前在从头蛋白质结构预测方面主要有两个挑战：一方面是巨大的构象搜索空间，另一方面是缺少一个精确的能量函数去确定接近天然态的蛋白质模型。伴随着能量函数的发展，已经有一些方法被提出解决构象采样的问题。片段组装方法的是一个非常有效地方法来减少巨大的搜索空间，片段组装仅使用蛋白质骨架原子的扭转角来简化的表示蛋白质的三维结构，首先将整条序列分割成多个9-mer或者3-mer片段，片段之间允许相互重叠，对每个片段，从PSI-BLAST计算出的相似序列对应的结构中，截取出相应位置的局部结构作为候选结构，然后使用MonteCarlo算法从每个片段的候选结构集中挑选出一个进行组合，使得组合成的全长结构能量最小。因此，在片段组装方法中，片段的质量是十分重要的。

在蛋白质结构从头预测方法中，由于能量函数的不精确以及采样能力的不足会导致构象空间的搜索效率低、收敛速度慢，局部搜索能力弱的问题，从而影响预测精度。比如蛋白质结构的Loop区域灵活多变，使用常规的片段组装方法不能对此部分进行充分的采样，因此，如何加强对蛋白质结构的Loop区域的局部搜索是我们需要解决的问题。

发明内容：

为了解决现有能量函数的不精确和构象空间搜索不充分导致的蛋白质结构预测精度较低的问题，本发明提出了一种预测精度较高的基于Loop区域高斯扰动的群体蛋白质结构预测方法。

本发明解决其技术问题所采用的技术方案是：

一种基于Loop区域高斯扰动的群体蛋白质结构预测方法，所述方法包括以下步骤：

1)参数设置：设置初始种群大小NP，交叉概率CR，种群交叉代数G，最大迭代次数G_max，温度参数KT,扰动参数a，设置能量函数；