1.一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法,其特征在于:所述蛋白质结构预测方法包括以下步骤:
1)选取蛋白质力场模型,即能量函数E(X),并输入待测蛋白质的序列信息;
2)初始化:设置种群规模NP,交叉概率CR,增益常数F,支撑面斜率控制因子M和最大迭代次数Gmax;根据输入的序列信息生成初始构象种群并计算每个构象个体的能量值f(Xi)=E(Xi),i=1,2,…,NP,其中,N表示维数, 表示第i个构象个体Xi的第N维元素,并初始化迭代次数G=0;
3)对种群中的每个构象Xi,i∈{1,2,…,NP}作如下处理:
3.1)从当前种群中随机选取四个互不相同的构象个体Xa,Xb,Xc,Xd,其中a≠b≠c≠d≠i;
3.2)根据变异操作 生成变异构象个体
3.3)根据变异操作 生成变异构象个体
3.4)根据能量值对当前种群中的构象个体进行升序排列,找出能量最低的构象个体Xbest;
3.5)根据变异操作 生成变异构象个体
3.6)分别对变异构象个体 和 执行交叉操作 生成新构象个体 和 其中m=1,2,3,j=1,2,…,N, 表示新构象个体 的第j维元素,表示变异构象个体 的第j维元素,表示目标构象个体Xi的第j维元素,R1表示0和1之间的随机小数,R2表示1和N之间的随机整数;
3.7)根据如下操作分别计算新构象个体 和 的能量估计值:
3.7.1)计算当前种群中每个构象个体到新构象个体 之间的距离表示第i个构象个体到第1个新构象个体 之间的距离;
3.7.2)根据距离对所有构象个体进行升序排列,选取距离最小的两个构象个体并记为并计算其抽象凸下界估计支撑面:
其中, 为所选构象个体 的能量函数值, 为所选构象个体 的第n维元素, 为辅助变量,M为支撑面斜率控制因子;
3.7.3)计算新构象个体 的能量估计值 其中max表示求最大值,min表示求最小值, 为支撑向量lt的第j维元素;
3.7.4)根据步骤3.7.1)-3.7.3)计算新构象个体 的能量估计值
3.7.5)根据步骤3.7.1)-3.7.3)计算新构象个体 的能量估计值
3.8)比较新构象个体 和 的能量估计值,并选出能量估计值最小的新构象个体,并记为Umin;
3.9)计算构象Umin的能量函数值f(Umin)=E(Umin),如果f(Umin)小于当前目标构象个体Xi的能量函数值f(Xi),则构象Umin替换构象Xi;
4)判断是否满足终止条件,若满足则输出结果并退出,否则返回步骤3)。
2.如权利要求1所述的一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法,其特征在于:所述步骤4)中,对种群中的每一个个体都执行完步骤3)以后,迭代次数G=G+1,终止条件为迭代次数G达到最大迭代次数Gmax。