[发明专利]一种残基接触引导loop扰动的群体蛋白质结构预测方法

说明书

一种残基接触引导loop扰动的群体蛋白质结构预测方法，首先，采用片段重组和片段组装策略，广泛的搜索构象空间；然后，利用预测的残基接触图实时构建loop区域扰动模型，利用差分进化算法快速求解扰动量，生成结构多样化的候选构象，进一步改善二级结构间的位置关系；最终，提高整体结构的精度和效率。本发明提供一种预测精度较高的残基接触引导loop扰动的群体蛋白质结构预测方法。

技术领域

本发明涉及生物信息学、计算机应用领域，尤其涉及的是一种残基接触引导loop扰动的群体蛋白质结构预测方法。

背景技术

蛋白质是生命活动功能的执行者，从新陈代谢到疾病免疫，都离不开蛋白质的作用。蛋白质拥有复杂的结构，具有完整一级结构的多肽或蛋白质,只有当其折叠形成正确的三维结构才能具有正常的生物学功能。然而，蛋白质折叠成复杂的空间结构并不是那么容易，常常伴随着错误折叠。错误折叠的蛋白质不仅丧失了生物学功能，还会引起疾病，疯牛病、老年痴呆症、帕金森病等疾病就是因为蛋白质错误折叠引起的。因此高效获取蛋白质结构对于了解其生物功能、相关疾病的药物研发是非常关键的。

目前，主要是通过实验测定的方法获取蛋白质三维结构。实验测定蛋白质结构的方法主要包括X射线晶体衍射、核磁共振、冷冻电镜技术。这类方法十分复杂，所需的时间成本和资金投入极高，并且不能适用于所有类型的蛋白质。

已知蛋白质结构中蕴含着大量的有效信息，随着信息处理技术的快速发展这些信息不断被挖掘，这使得利用计算机模拟蛋白质折叠过程、预测蛋白质三维结构成为可能。全世界范围内有许多研究机构致力于研究利用生物大数据、人工智能、系统优化技术预测蛋白质三维结构，并逐渐应用于疾病诊断和药物设计，其中具有代表性的研究团队有美国华盛顿的David Baker实验室、美国密歇根大学的张阳实验室等。国内越来越多的高校和研究机构也加入到蛋白质结构预测的研究当中。

蛋白质结构中的loop区域连接着α螺旋和β折叠区域，结构十分灵活。Loop区域微小的变动有可能会造成整体结构大的变动。目前的蛋白质结构预测方法对loop区域的探索能力不够强，并且片段库在有效减小构象搜索空间的同时有可能没有覆盖到某些重要的区域，导致无法搜索到正确的结构。

因此，目前的蛋白质结构预测方法对loop区域的搜索能力不足，并且受到片段库的约束，需要改进。

发明内容

为了解决现有的蛋白质结构预测方法对loop区域探索能力不足导致整体结构精度不高的问题，本发明提出了一种残基接触引导loop扰动的群体蛋白质结构预测方法，增强对loop区域的探索，本发明提出的局部扰动策略能探索到片段库之外的构象空间，提高了loop区域结构的多样性，最终达到提高整体结构精度的目的。

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

一种残基接触引导loop扰动的群体蛋白质结构预测方法，包括以下步骤：

1)输入预测蛋白质的目标序列、9残基片段库和预测的残基接触图；

2)设置参数：种群规模NP、迭代次数G；

3)种群初始化：随机选择目标序列的滑动窗口进行片段组装，直至所有残基都至少被替换过一次，迭代此过程生成初始种群P＝{P₁,P₂,...,P_NP}，其中P_n表示种群P中的第n个蛋白质构象，n∈{1,2,...,NP}；

4)设g＝1，其中g∈{1,2,...,G}；

5)设n＝1，其中n∈{1,2,...,NP}；

6)片段重组，过程如下：