CN106650305A - 一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法 - Google Patents
一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法 Download PDFInfo
- Publication number
- CN106650305A CN106650305A CN201610884824.9A CN201610884824A CN106650305A CN 106650305 A CN106650305 A CN 106650305A CN 201610884824 A CN201610884824 A CN 201610884824A CN 106650305 A CN106650305 A CN 106650305A
- Authority
- CN
- China
- Prior art keywords
- conformation
- individuals
- individual
- conformational
- new
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000000455 protein structure prediction Methods 0.000 title claims abstract description 17
- 230000001174 ascending effect Effects 0.000 claims abstract description 10
- 102000004169 proteins and genes Human genes 0.000 claims description 34
- 108090000623 proteins and genes Proteins 0.000 claims description 34
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000005457 optimization Methods 0.000 description 8
- 125000003275 alpha amino acid group Chemical group 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 230000002028 premature Effects 0.000 description 2
- 102000015636 Oligopeptides Human genes 0.000 description 1
- 108010038807 Oligopeptides Proteins 0.000 description 1
- 238000002441 X-ray diffraction Methods 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004883 computer application Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 102000038037 druggable proteins Human genes 0.000 description 1
- 108091007999 druggable proteins Proteins 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000324 molecular mechanic Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Biology (AREA)
- Chemical & Material Sciences (AREA)
- Medical Informatics (AREA)
- Crystallography & Structural Chemistry (AREA)
- Analytical Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法,首先,根据随机选择的构象个体、当前目标构象个体和能量值最低的构象构体生成三个不同的新构象个体;然后,分别计算当前种群中各构象个体到各新构象个体的距离,并根据距离进行升序排列;其次,计算离各新构象个体最近的部分构象个体的抽象凸支撑面,从而计算各新构象个体的能量估计值;最后,比较各新构象个体能量估计值,从而选取能量估计值最低的新构象个体进行能量函数评价。本发明提出一种预测精度高、搜索效率高的基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法。
Description
技术领域
本发明涉及一种生物学信息学、智能优化、计算机应用领域,尤其涉及的是,一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法。
背景技术
生物细胞中包含许多由20多种氨基酸所形成的长链折叠而成的蛋白质,蛋白质结构预测问题是当今计算生物学领域中的研究热点,不仅具有非常重要的理论机制,而且对新蛋白的设计、蛋白质之间相互作用建模、药物标靶蛋白的设计具有十分重要的指导意义。目前,最常用的蛋白质结构实验测定方法包括X-晶体衍射和核磁共振,但是上述两种蛋白质结构测定方法存在一定的缺陷,无法满足所有蛋白质结构预测的需求。例如,对于一些不易结晶的蛋白,无法使用X-晶体衍射方法进行测定;使用核磁共振房产测定一个蛋白的结构费用较大(大约15万美元),并且极其费时(大约需要半年)。因此,如何以计算机为工具,运用适当的算法,从氨基酸序列出发直接预测蛋白质的三维结构,成为当前生物信息学中一种重要的研究课题。
1965年,Anfinsen等提出了蛋白质的一级结构完全决定了其三维空间结构的著名论断,这一论断使得根据蛋白质得氨基酸序列从理论上预测其相应的空间结构成为了现实。从热力学角度来看,也就是说蛋白质的天然态结构通常对应着能量最低的构象。因此,从氨基酸序列出发,以能量函数来衡量构象的质量,采用适当的算法搜索能量模型中最小能量构象,从而预测蛋白质的天然态结构,已成为计算生物信息学中重要的研究课题之一。目前,最常用的预测方法可以划分为以下三类:针对高相似序列的同源建模方法;针对较低相似性序列的折叠识别方法;以及不依赖模板的从头预测方法(ab initio或de novo)。经过20多年的发展历程表明,对于序列相似度较高的情况(>50%),TBM预测精度能够达到左右;然而,序列相似度<30%的情况下并不理想;对于序列相似度<20%或寡肽(<10个残基的小蛋白)来说,从头预测方法是唯一的选择。
从头预测方法直接基于Anfinsen假说建立蛋白质物理或知识能量模型,然后设计适当优化算法求解最小能量构象。可以看出,从头预测方法必须考虑以下两个因素:(1)知识能量的构建;(2)构象空间搜索方法。第一个因素本质上属于分子力学问题,主要是为了能够计算得到每个蛋白质结构对应的能量值。第二个因素本质上属于全局优化问题,通过选择一种合适的优化方法,对构象空间进行快速搜索,得到能量最低构象。其中,蛋白质构象空间优化属于一类非常难解的NP-Hard问题。群体进化类算法是研究蛋白质分子构象优化的重要方法,主要包括差分进化算法(DE)、遗传算法(GA)、粒子群算法(PSO),这些算法不仅结构简单,易于实现,而且鲁棒性强,因此,经常被用于从头预测方法中的全局最小能量构象搜索。然而,对于上述群体算法,新构象生成策略的选择至关重要,不恰当的策略会导致算法搜索效率低、收敛速度慢,甚至陷入局部最优,出现早熟收敛现象,从而影响预测精度。
因此,现有的构象空间优化方法在搜索效率和预测精度方面存在着缺陷,需要改进。
发明内容
为了克服现有的群体蛋白质结构预测方法在预测精度和搜索效率方面的不足,本发明提出一种预测精度高、搜索效率高的基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法。
本发明解决其技术问题所采用的技术方案是:
一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法,所述方法包括以下步骤:
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)根据变异操作V1 i=Xa+F(Xb-Xc)生成变异构象个体V1 i;
3.3)根据变异操作生成变异构象个体
3.4)根据能量值对当前种群中的构象个体进行升序排列,找出能量最低的构象个体Xbest;
3.5)根据变异操作生成变异构象个体
3.6)分别对变异构象个体V1 i、和执行交叉操作生成新构象个体和其中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=1,2,…,N+1为所选构象个体的第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)。
进一步,所述步骤4)中,对种群中的每个构象个体都执行完步骤3)以后,迭代次数G=G+1,终止条件为迭代次数G达到预设最大迭代次数Gmax。
本发明的技术构思为:首先,根据随机选择的构象个体、当前目标构象个体和能量值最低的构象构体生成三个不同的新构象个体;然后,分别计算当前种群中各构象个体到各新构象个体的距离,并根据距离进行升序排列;其次,计算离各新构象个体最近的部分构象个体的抽象凸支撑面,从而计算各新构象个体的能量估计值;最后,比较各新构象个体能量估计值,从而选取能量估计值最低的新构象个体进行能量函数评价。
本发明的有益效果表现在:一方面,根据不同的策略生成多个不同的新构象个体,避免因策略选择不恰当而引起的早熟收敛,从而提高预测精度;另一方面,根据新构象个体的抽象凸能量估计值选择新构象构体,减少能量函数平均,从而降低计算代价,提高搜索效率。
附图说明
图1是基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法的流程图。
图2是基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法对蛋白质1AIL进行结构预测时的构象更新示意图。
图3是基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法对蛋白质1AIL进行结构预测时得到的构象分布图。
图4是基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法对蛋白质1AIL进行结构预测得到的三维结构图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图4,一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法,包括以下步骤:
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)根据变异操作V1 i=Xa+F(Xb-Xc)生成变异构象个体V1 i;
3.3)根据变异操作生成变异构象个体
3.4)根据能量值对当前种群中的构象个体进行升序排列,找出能量最低的构象个体Xbest;
3.5)根据变异操作生成变异构象个体
3.6)分别对变异构象个体V1 i、和执行交叉操作生成新构象个体和其中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=1,2,…,N+1为所选构象个体的第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)。
进一步,所述步骤4)中,对种群中的每个构象个体都执行完步骤3)以后,迭代次数G=G+1,终止条件为迭代次数G达到预设最大迭代次数Gmax。
本实施例序列长度为73的α折叠蛋白质1AIL为实施例,一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法,其中包含以下步骤:
1)选取蛋白质Rosetta Score3力场模型,即Rosetta Score3能量函数E(X),并输入待测蛋白质的序列信息;
2)初始化:设置种群规模NP=50,交叉概率CR=0.5,增益常数F=0.5,支撑面斜率控制因子M=1000,最大迭代次数Gmax=10000;根据输入的序列信息生成初始构象种群并计算每个构象个体的能量值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)根据变异操作V1 i=Xa+F(Xb-Xc)生成变异构象个体V1 i;
3.3)根据变异操作生成变异构象个体
3.4)根据能量值对当前种群中的构象个体进行升序排列,找出能量最低的构象个体Xbest;
3.5)根据变异操作生成变异构象个体
3.6)分别对变异构象个体V1 i、和执行交叉操作生成新构象个体和其中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=1,2,…,N+1为所选构象个体的第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)以后,迭代次数G=G+1,判断迭代次数G是否大于最大迭代次数Gmax,若G大于Gmax,则输出预测结果并退出,否则返回步骤3)。
以序列长度为73的α折叠蛋白质1AIL为实施例,运用以上方法得到了该蛋白质的近天然态构象,最小均方根偏差为平均均方根偏差为预测结构如图4所示。
以上说明是本发明以1AIL蛋白质为实例所得出的优化效果,并非限定本发明的实施范围,在不偏离本发明基本内容所涉及范围的的前提下对其做各种变形和改进,不应排除在本发明的保护范围之外。
Claims (2)
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。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610884824.9A CN106650305B (zh) | 2016-10-10 | 2016-10-10 | 一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610884824.9A CN106650305B (zh) | 2016-10-10 | 2016-10-10 | 一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106650305A true CN106650305A (zh) | 2017-05-10 |
CN106650305B CN106650305B (zh) | 2019-01-22 |
Family
ID=58853887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610884824.9A Active CN106650305B (zh) | 2016-10-10 | 2016-10-10 | 一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106650305B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108062457A (zh) * | 2018-01-15 | 2018-05-22 | 浙江工业大学 | 一种结构特征向量辅助选择的蛋白质结构预测方法 |
CN109285585A (zh) * | 2018-08-29 | 2019-01-29 | 浙江工业大学 | 一种基于动态抽象凸下界估计的群体蛋白质结构预测方法 |
CN109346126A (zh) * | 2018-08-29 | 2019-02-15 | 浙江工业大学 | 一种下界估计策略自适应蛋白质结构预测方法 |
CN109411013A (zh) * | 2018-08-29 | 2019-03-01 | 浙江工业大学 | 一种基于个体特定变异策略的群体蛋白质结构预测方法 |
CN110610742A (zh) * | 2019-09-20 | 2019-12-24 | 福建工程学院 | 一种基于蛋白质互作网络的功能模块检测方法 |
CN113035268A (zh) * | 2021-04-09 | 2021-06-25 | 上海交通大学 | 基于多目标分解优化策略的蛋白质结构优化方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0793286A (ja) * | 1993-09-20 | 1995-04-07 | Nec Corp | タンパク質立体構造予測方法 |
US5958784A (en) * | 1992-03-25 | 1999-09-28 | Benner; Steven Albert | Predicting folded structures of proteins |
CN103413067A (zh) * | 2013-07-30 | 2013-11-27 | 浙江工业大学 | 一种基于抽象凸下界估计的蛋白质结构预测方法 |
CN103970985A (zh) * | 2014-04-08 | 2014-08-06 | 浙江工业大学 | 一种基于局部抽象凸估计的群体全局优化方法 |
CN103984878A (zh) * | 2014-04-08 | 2014-08-13 | 浙江工业大学 | 一种基于树搜索和片段组装的蛋白质结构预测方法 |
CN104102955A (zh) * | 2014-07-23 | 2014-10-15 | 杭州中威智能科技有限公司 | 一种基于抽象凸估计选择差分进化算法的电力线路规划方法 |
CN104182653A (zh) * | 2014-08-25 | 2014-12-03 | 浙江工业大学 | 一种基于局部抽象凸低估支撑面的群体构象空间优化方法 |
CN105787292A (zh) * | 2014-12-18 | 2016-07-20 | 中国科学院大连化学物理研究所 | 蛋白质折叠的并行预测方法 |
US20160210399A1 (en) * | 2012-05-09 | 2016-07-21 | Memorial Sloan-Kettering Cancer Center | Methods and apparatus for predicting protein structure |
-
2016
- 2016-10-10 CN CN201610884824.9A patent/CN106650305B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5958784A (en) * | 1992-03-25 | 1999-09-28 | Benner; Steven Albert | Predicting folded structures of proteins |
JPH0793286A (ja) * | 1993-09-20 | 1995-04-07 | Nec Corp | タンパク質立体構造予測方法 |
US20160210399A1 (en) * | 2012-05-09 | 2016-07-21 | Memorial Sloan-Kettering Cancer Center | Methods and apparatus for predicting protein structure |
CN103413067A (zh) * | 2013-07-30 | 2013-11-27 | 浙江工业大学 | 一种基于抽象凸下界估计的蛋白质结构预测方法 |
CN103970985A (zh) * | 2014-04-08 | 2014-08-06 | 浙江工业大学 | 一种基于局部抽象凸估计的群体全局优化方法 |
CN103984878A (zh) * | 2014-04-08 | 2014-08-13 | 浙江工业大学 | 一种基于树搜索和片段组装的蛋白质结构预测方法 |
CN104102955A (zh) * | 2014-07-23 | 2014-10-15 | 杭州中威智能科技有限公司 | 一种基于抽象凸估计选择差分进化算法的电力线路规划方法 |
CN104182653A (zh) * | 2014-08-25 | 2014-12-03 | 浙江工业大学 | 一种基于局部抽象凸低估支撑面的群体构象空间优化方法 |
CN105787292A (zh) * | 2014-12-18 | 2016-07-20 | 中国科学院大连化学物理研究所 | 蛋白质折叠的并行预测方法 |
Non-Patent Citations (3)
Title |
---|
XIAO-HU HAO ET AL;: "《A Novel Method Using Abstract Convex Underestimation in Ab-Initio Protein Structure Prediction for Guiding Search in Conformational Feature Space》", 《IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS,》 * |
张贵军 等;: "《动态步长蛋白质构象空间搜索方法》", 《吉林大学学报(工学版)》 * |
张贵军 等;: "《基于广义凸下界估计的多模态差分进化算法》", 《软件学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108062457A (zh) * | 2018-01-15 | 2018-05-22 | 浙江工业大学 | 一种结构特征向量辅助选择的蛋白质结构预测方法 |
CN108062457B (zh) * | 2018-01-15 | 2021-06-18 | 浙江工业大学 | 一种结构特征向量辅助选择的蛋白质结构预测方法 |
CN109285585A (zh) * | 2018-08-29 | 2019-01-29 | 浙江工业大学 | 一种基于动态抽象凸下界估计的群体蛋白质结构预测方法 |
CN109346126A (zh) * | 2018-08-29 | 2019-02-15 | 浙江工业大学 | 一种下界估计策略自适应蛋白质结构预测方法 |
CN109411013A (zh) * | 2018-08-29 | 2019-03-01 | 浙江工业大学 | 一种基于个体特定变异策略的群体蛋白质结构预测方法 |
CN109346126B (zh) * | 2018-08-29 | 2020-10-30 | 浙江工业大学 | 一种下界估计策略自适应蛋白质结构预测方法 |
CN109411013B (zh) * | 2018-08-29 | 2020-10-30 | 浙江工业大学 | 一种基于个体特定变异策略的群体蛋白质结构预测方法 |
CN109285585B (zh) * | 2018-08-29 | 2021-05-18 | 浙江工业大学 | 一种基于动态抽象凸下界估计的群体蛋白质结构预测方法 |
CN110610742A (zh) * | 2019-09-20 | 2019-12-24 | 福建工程学院 | 一种基于蛋白质互作网络的功能模块检测方法 |
CN110610742B (zh) * | 2019-09-20 | 2023-12-19 | 福建工程学院 | 一种基于蛋白质互作网络的功能模块检测方法 |
CN113035268A (zh) * | 2021-04-09 | 2021-06-25 | 上海交通大学 | 基于多目标分解优化策略的蛋白质结构优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106650305B (zh) | 2019-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106650305B (zh) | 一种基于局部抽象凸支撑面的多策略群体蛋白质结构预测方法 | |
Senior et al. | Protein structure prediction using multiple deep neural networks in the 13th Critical Assessment of Protein Structure Prediction (CASP13) | |
AU2017338775B2 (en) | Phenotype/disease specific gene ranking using curated, gene library and network based data structures | |
CN106503484B (zh) | 一种基于抽象凸估计的多阶段差分进化蛋白质结构预测方法 | |
CN107609342B (zh) | 一种基于二级结构空间距离约束的蛋白质构象搜索方法 | |
Lin et al. | Efficient classification of hot spots and hub protein interfaces by recursive feature elimination and gradient boosting | |
CN109785901B (zh) | 一种蛋白质功能预测方法及装置 | |
Hoksza et al. | Efficient RNA pairwise structure comparison by SETTER method | |
Bartonek et al. | mRNA/protein sequence complementarity and its determinants: The impact of affinity scales | |
Wang et al. | GAEM: a hybrid algorithm incorporating GA with EM for planted edited motif finding problem | |
CN109360596B (zh) | 一种基于差分进化局部扰动的蛋白质构象空间优化方法 | |
Sun et al. | Protein function prediction using function associations in protein–protein interaction network | |
Dou et al. | Accurate identification of RNA D modification using multiple features | |
CN106096326A (zh) | 一种基于质心变异策略的差分进化蛋白质结构预测方法 | |
Kim et al. | Probabilistic models for semisupervised discriminative motif discovery in DNA sequences | |
Wang et al. | PLANNER: a multi-scale deep language model for the origins of replication site prediction | |
CN108920894B (zh) | 一种基于简约抽象凸估计的蛋白质构象空间优化方法 | |
Hao et al. | Conformational space sampling method using multi-subpopulation differential evolution for de novo protein structure prediction | |
CN109360597B (zh) | 一种基于全局和局部策略协作的群体蛋白质结构预测方法 | |
Yang et al. | PseKNC and Adaboost-based method for DNA-binding proteins recognition | |
Zhang et al. | A Fast Cluster Motif Finding Algorithm for ChIP‐Seq Data Sets | |
CN109390035B (zh) | 一种基于局部结构比对的蛋白质构象空间优化方法 | |
Basha Gutierrez et al. | A genetic algorithm for motif finding based on statistical significance | |
CN109411013B (zh) | 一种基于个体特定变异策略的群体蛋白质结构预测方法 | |
Li et al. | Improving domain-based protein interaction prediction using biologically-significant negative dataset |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |