[发明专利]一种基于多样性指标的群体蛋白质构象空间优化方法

权利要求书

1.一种基于多样性指标的群体蛋白质构象空间优化方法，其特征在于，所述优化方法包括以下步骤：

1)给定输入序列信息；

2)设置参数：种群大小NP，重组因子RC，最大迭代代数G_max，多样性接收概率p；

3)利用PSIPRED预测查询序列的二级结构信息；

4)初始化：首先通过用Rosetta第一阶段进行片段组装生成一个规模为NP的初始种群P＝{x_i|i∈I}，其中i为种群个体编号，I为种群个体编号集合，I＝{1,2,...,NP}，x_i表示第i个个体；

5)采用Rosetta第二阶段的局部搜索策略对每个种群个体进行片段组装形成新的种群；

6)构象重组，过程如下：

6.1)从当前种群中选择两个不同的个体x_a和x_b，其中，a,b∈{1,2,...,NP},a≠b；根据二级结构信息随机选取一个loop区域，并在0和1之间随机生成一个均匀分布随机数rand1；

6.2)如果rand1RC，则针对步骤6.1)中选择的loop区域，交换构象x_a和x_b的二面角分别生成新的构象和

6.3)重复步骤6.1)-6.2)直到种群中的所有个体都进行了构象重组为止，生成新的种群

7)针对种群中的每一个构象i∈I进行变异，过程如下：

7.1)统计所有loop区域能够进行片段组装位置的数量W_L，并计算变异概率

7.2)生成一个0和1之间的均匀分布随机数rand2，如果rand2P_m，则对构象的所有loop区域进行一次片段组装生成变异个体

7.3)按照步骤7.1)-7.2)对当前种群中的所有构象进行变异后，生成新的种群

8)合并种群和P中的所有个体，即

9)针对中的每一个构象i∈I，计算多样性指标，过程如下：

9.1)根据每个构象的C_α原子坐标计算每个构象的质心

9.2)计算每个C_α原子到质心的欧氏距离，并记离质心最远的原子为离质心最近的原子为

9.3)计算每个C_α原子到的欧氏距离，并记离最远的原子为

9.4)分别计算步骤9.2)中所有距离的平均值偏差偏度

9.5)分别计算步骤9.3)中所有距离的平均值偏差偏度

9.6)计算每个C_α原子到的距离，并计算这些距离的平均值偏差偏度

9.7)计算每个C_α原子到的距离，并计算这些距离的平均值偏差偏度

9.8)将步骤9.4)-9.7)中的距离平均值、偏差和偏度合并得到构象的特征向量

9.9)计算多样性指标

其中，T＝{1,2,…,2NP}；表示向量U_i的第l维元素；

10)针对种群P中的每一个个体x_i,i∈I进行更新，过程如下：

10.1)分别根据Rosetta score2计算个体x_i和的能量score(x_i)和

10.2)如果则个体替换父代个体x_i，否则继续步骤10.3)；

10.3)如果D_i＞D′_i，则产生一个[0,1]的随机数rand3，如果rand3p，则个体替换父代个体x_i；

11)对步骤6)-10)进行迭代直到达到最大迭代代数G_max，然后继续步骤12)；

12)采用Rosetta第三阶段的局部搜索策略对每个种群个体进行片段组装形成新的种群，并根据步骤6)-10)进行迭代直到达到最大迭代代数G_max，其中能量函数采用Rosettasocre3；

13)采用Rosetta第四阶段的局部搜索策略对每个种群个体进行片段组装形成新的种群，并根据步骤6)-10)进行迭代直到达到最大迭代代数G_max，其中能量函数采用Rosettasocre4；

14)输出结果。