基于圆周统计的蛋白质折叠的表征方法及系统的制作方法
【专利摘要】本发明涉及生物信息学【技术领域】，特别涉及一种基于圆周统计的蛋白质折叠的表征方法及系统。本发明的方法基于圆周统计对蛋白质折叠进行表征，高效、准确，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理。
【专利说明】基于圆周统计的蛋白质折叠的表征方法及系统
【【技术领域】】
[0001]本发明涉及生物信息学【技术领域】，特别涉及一种基于圆周统计的蛋白质折叠的表征方法及系统。
【【背景技术】】
[0002]神经退行性疾病给人类造成了巨大的社会和经济负担。比较常见的神经退行性疾病有老年痴呆、帕金森氏病、小脑萎缩症、多聚谷氨酰胺疾病等。美国500多万老年痴呆症患者每年的花费是1830亿美元；而中国的老年痴呆症患者也在500-600万之间。研究发现蛋白质折叠的研究对探索神经退行性疾病的致病机理意义重大；而且因有助于蛋白质分子药物的设计，蛋白质折叠研究对“折叠病”的预防和治疗也将起到重大的帮助作用。
[0003]传统的蛋白质折叠的研究是基于笛卡尔坐标系；用该坐标系表示蛋白需要更多的变量；用角度坐标系来表示和研究蛋白质折叠也比较常用，但是由于角度周期性的特点，该方法的研究受到一些局限，效率和准确率均较低。因此，提出一种新的方法研究蛋白质折叠，提供效率及准确性，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理，成为亟待解决的问题。
【
【发明内容】
】
[0004]本发明的目的在于针对现有技术的上述缺陷，提供一种基于圆周统计的蛋白质折叠的表征方法，解决现有技术中蛋白质折叠研究效率和准确率较低的技术问题。
[0005]本发明的目的可通过以下技术措施来实现:
[0006]一种基于圆周统计的蛋白质折叠的表征方法，包括如下步骤:
[0007]步骤1:建立第一圆周统`计模型，用角度方差V(a)=l_R(a)表征热动力学涨落；建
立第二圆周统计模型，用角度协方差=《(S厂⑷).(+? 一〈〃,〉)>:表征热动力学涨落；
[0008]步骤2:用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析；
[0009]步骤3:将待研究蛋白的能量力场分割为侧链与主链的相互作用能量Es_b、不同侧链之间的相互作用能量Es_s、不同主链之间的相互作用能量Eb_b;分析主链上第m个氨基酸残基和第η个氨基酸残基之间的侧链相互作用能量的涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第η个氨基酸残基进行标记；
[0010]步骤4:分别用第一圆周统计模型和第二圆周统计模型分析不同侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。
[0011]优选地，在步骤I之前还包括如下步骤:
[0012]导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值。
[0013]优选地，所述方法还包括如下步骤:
[0014]将待研究蛋白质的侧链分类亲水侧链和疏水侧链，并分别对亲水侧链和疏水侧链进行标记，分别用第一圆周统计模型和第二圆周统计模型分析不同亲水侧链之间的二面角相互作用热涨落、亲水侧链与疏水侧链之间的二面角相互作用热涨落、不同疏水侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。
[0015]优选地，所述方法还包括:用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中不同亲水侧链之间的二面角相互作用、亲水侧链与疏水侧链之间的二面角相互作用、不同疏水侧链之间的二面角相互作用进行热动力学涨落分析。
[0016]优选地，所述步骤3还包括:
[0017]对待研究蛋白质序列上每个氨基酸残基进行编号，分析不同氨基酸残基之间的不同侧链、不同主链或侧链和主链之间的相互作用能量涨落。
[0018]本发明还提供了一种基于圆周统计的蛋白质折叠的表征系统，包括:
[0019]模型建立模块，用于建立第一圆周统计模型和第二圆周统计模型，所述第一圆周统计模型用角度方差V(a)=l_R(a)表征热动力学涨落；所述第二圆周统计模型用角度协方差= {{If - (Sd).(Sj -两> ft表征热动力学涨落；
[0020]化学键热涨落分析模块，用于用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析；
[0021]能量热涨落分析模块，用于将待研究蛋白的能量力场分割为侧链与主链的相互作用能量Es_b、不同侧链之间的相互作用能量Es_s、不同主链之间的相互作用能量Eb_b ；以及
[0022]分析主链上第m个氨基酸残基和第η个氨基酸残基之间的侧链相互作用能量 啲涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第
η个氨基酸残基进行标记；
[0023]二面角相互作用热涨落分析模块，用于分别用第一圆周统计模型和第二圆周统计模型分析不同侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；以及
[0024]当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。
[0025]优选地，所述系统还包括:
[0026]预处理模块，用于导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键或二硫键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值。
[0027]优选地，所述二面角相互作用热涨落分析模块还用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中不同亲水侧链之间的二面角相互作用、亲水侧链与疏水侧链之间的二面角相互作用、不同疏水侧链之间的二面角相互作用进行热动力学涨落分析。[0028]本发明的有益效果在于，本发明的方法基于圆周统计对蛋白质折叠进行表征，高效、准确，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理。
【【专利附图】
【附图说明】】
[0029]图1是本发明实施例的基于圆周统计的蛋白质折叠的表征方法的流程图；
[0030]图2是本发明实施例的基于圆周统计的蛋白质折叠的表征系统的结构框图。
【【具体实施方式】】
[0031]为了使本发明的目的、技术方案及优点更加清楚明白，下面结合附图和具体实施例对本发明作进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
[0032]本发明利用圆周统计对蛋白质折叠进行表征，高效、准确，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理。
[0033]本发明实施例还提供了一种基于圆周统计的蛋白质折叠的表征方法，如图1所示，包括如下步骤:
[0034]步骤SlOl，预处理步骤:
[0035]导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值；
[0036]步骤S102，建立两种圆周统计模型:
[0037]建立第一圆周统计模型，用角度方差V(a)=l_R(a)表征热动力学涨落；建立第二圆周统计模型，用角度协方差(6(--) = <(?,.——钩>))表征热动力学涨落；
[0038]相比于线性统计学，圆周统计学更适用于解决蛋白质折叠表征的技术问题。但是圆周统计方差和协方差的数学模型在蛋白质折叠研究中还存在不一致性，具体分析如下:
圆周统计中角度α可以用二维向量表示为S=角度方差可以用来研究蛋白质
折叠过程中物理量的热动力学涨落，其公式为:
[0039]
【权利要求】
1.一种基于圆周统计的蛋白质折叠的表征方法，其特征在于，包括如下步骤: 步骤1:建立第一圆周统计模型，用角度方差V(a)=l-R(a)表征热动力学涨落；建立第二圆周统计模型,用角度协方差(Ij) = φη - (Hi)Uffj - <々>)〉:表征热动力学涨落； 步骤2:用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析； 步骤3:将待研究蛋白的能量力场分割为侧链与主链的相互作用能量Es_b、不同侧链之间的相互作用能量Es_s、不同主链之间的相互作用能量Eb_b;分析主链上第m个氨基酸残基和第η个氨基酸残基之间的侧链相互作用能量啲涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第η个氨基酸残基进行标记； 步骤4:分别用第一圆周统计模型和第二圆周统计模型分析不同侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。
2.根据权利要求1所述的基于圆周统计的蛋白质折叠的表征方法，其特征在于，在步骤I之前还包括如下步骤: 导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号， 分别设定氢键或二硫键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值。
3.根据权利要求1所述的基于圆周统计的蛋白质折叠的表征方法，其特征在于，还包括如下步骤: 将待研究蛋白质的侧链分类亲水侧链和疏水侧链，并分别对亲水侧链和疏水侧链进行标记，分别用第一圆周统计模型和第二圆周统计模型分析不同亲水侧链之间的二面角相互作用热涨落、亲水侧链与疏水侧链之间的二面角相互作用热涨落、不同疏水侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。
4.根据权利要求3所述的基于圆周统计的蛋白质折叠的表征方法，其特征在于，还包括:用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中不同亲水侧链之间的二面角相互作用、亲水侧链与疏水侧链之间的二面角相互作用、不同疏水侧链之间的二面角相互作用进行热动力学涨落分析。
5.根据权利要求2所述的基于圆周统计的蛋白质折叠的表征方法，其特征在于，所述步骤3还包括: 对待研究蛋白质序列上每个氨基酸残基进行编号，分析不同氨基酸残基之间的不同侧链、不同主链或侧链和主链之间的相互作用能量涨落。
6.一种基于圆周统计的蛋白质折叠的表征系统，其特征在于，包括: 模型建立模块，用于建立第一圆周统计模型和第二圆周统计模型，所述第一圆周统计模型用角度方差V(a)=l_R(a)表征热动力学涨落；所述第二圆周统计模型用角度协方差 %.) = {{a； —— (?瓜表征热动力学涨落； 化学键热涨落分析模块，用于用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析； 能量热涨落分析模块，用于将待研究蛋白的能量力场分割为侧链与主链的相互作用能量Es_b、不同侧链之间的相互作用能量Es_s、不同主链之间的相互作用能量Eb_b ；以及 分析主链上第m个氨基酸残基和第η个氨基酸残基之间的侧链相互作用能量的涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第η个氨基酸残基进行标记； 二面角相互作用热涨落分析模块，用于分别用第一圆周统计模型和第二圆周统计模型分析不同侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；以及 当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。
7.根据权利要求6所述的基于圆周统计的蛋白质折叠的表征系统，其特征在于，还包括: 预处理模块，用于导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键或二硫键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值。
8.根据权利要求6所述的基于圆周统计的蛋白质折叠的表征系统，其特征在于，所述二面角相互作用热涨落分析模块还用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中不同亲水侧链之间的二面角相互作用、亲水侧链与疏水侧链之间的二面角相互作用、不同疏水侧链之间的二面角相互作用进行热动力学涨落分析。
【文档编号】G06F19/16GK103714264SQ201310683271
【公开日】2014年4月9日 申请日期:2013年12月12日 优先权日:2013年12月12日
【发明者】魏彦杰, 张慧玲, 彭丰斌, 孟金涛, 冯圣中 申请人:深圳先进技术研究院