基于圆周统计的蛋白质折叠的表征方法及系统技术方案

本发明专利技术涉及生物信息学技术领域，特别涉及一种基于圆周统计的蛋白质折叠的表征方法及系统。本发明专利技术的方法基于圆周统计对蛋白质折叠进行表征，高效、准确，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理。

【技术实现步骤摘要】
基于圆周统计的蛋白质折叠的表征方法及系统
本专利技术涉及生物信息学
，特别涉及一种基于圆周统计的蛋白质折叠的表征方法及系统。
技术介绍
神经退行性疾病给人类造成了巨大的社会和经济负担。比较常见的神经退行性疾病有老年痴呆、帕金森氏病、小脑萎缩症、多聚谷氨酰胺疾病等。美国500多万老年痴呆症患者每年的花费是1830亿美元；而中国的老年痴呆症患者也在500-600万之间。研究发现蛋白质折叠的研究对探索神经退行性疾病的致病机理意义重大；而且因有助于蛋白质分子药物的设计，蛋白质折叠研究对“折叠病”的预防和治疗也将起到重大的帮助作用。传统的蛋白质折叠的研究是基于笛卡尔坐标系；用该坐标系表示蛋白需要更多的变量；用角度坐标系来表示和研究蛋白质折叠也比较常用，但是由于角度周期性的特点，该方法的研究受到一些局限，效率和准确率均较低。因此，提出一种新的方法研究蛋白质折叠，提供效率及准确性，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理，成为亟待解决的问题。
技术实现思路
本专利技术的目的在于针对现有技术的上述缺陷，提供一种基于圆周统计的蛋白质折叠的表征方法，解决现有技术中蛋白质折叠研究效率和准确率较低的技术问题。本专利技术的目的可通过以下技术措施来实现：一种基于圆周统计的蛋白质折叠的表征方法，包括如下步骤：步骤1：建立第一圆周统计模型，用角度方差V(a)=1-R(a)表征热动力学涨落；建立第二圆周统计模型，用角度协方差表征热动力学涨落；步骤2：用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析；步骤3：将待研究蛋白的能量力场分割为侧链与主链的相互作用能量Es-b、不同侧链之间的相互作用能量Es-s、不同主链之间的相互作用能量Eb-b；分析主链上第m个氨基酸残基和第n个氨基酸残基之间的侧链相互作用能量的涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第n个氨基酸残基进行标记；步骤4：分别用第一圆周统计模型和第二圆周统计模型分析不同侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。优选地，在步骤1之前还包括如下步骤：导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值。优选地，所述方法还包括如下步骤：将待研究蛋白质的侧链分类亲水侧链和疏水侧链，并分别对亲水侧链和疏水侧链进行标记，分别用第一圆周统计模型和第二圆周统计模型分析不同亲水侧链之间的二面角相互作用热涨落、亲水侧链与疏水侧链之间的二面角相互作用热涨落、不同疏水侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。优选地，所述方法还包括：用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中不同亲水侧链之间的二面角相互作用、亲水侧链与疏水侧链之间的二面角相互作用、不同疏水侧链之间的二面角相互作用进行热动力学涨落分析。优选地，所述步骤3还包括：对待研究蛋白质序列上每个氨基酸残基进行编号，分析不同氨基酸残基之间的不同侧链、不同主链或侧链和主链之间的相互作用能量涨落。本专利技术还提供了一种基于圆周统计的蛋白质折叠的表征系统，包括：模型建立模块，用于建立第一圆周统计模型和第二圆周统计模型，所述第一圆周统计模型用角度方差V(a)=1-R(a)表征热动力学涨落；所述第二圆周统计模型用角度协方差表征热动力学涨落；化学键热涨落分析模块，用于用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析；能量热涨落分析模块，用于将待研究蛋白的能量力场分割为侧链与主链的相互作用能量Es-b、不同侧链之间的相互作用能量Es-s、不同主链之间的相互作用能量Eb-b；以及分析主链上第m个氨基酸残基和第n个氨基酸残基之间的侧链相互作用能量的涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第n个氨基酸残基进行标记；二面角相互作用热涨落分析模块，用于分别用第一圆周统计模型和第二圆周统计模型分析不同侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；以及当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。优选地，所述系统还包括：预处理模块，用于导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键或二硫键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值。优选地，所述二面角相互作用热涨落分析模块还用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中不同亲水侧链之间的二面角相互作用、亲水侧链与疏水侧链之间的二面角相互作用、不同疏水侧链之间的二面角相互作用进行热动力学涨落分析。本专利技术的有益效果在于，本专利技术的方法基于圆周统计对蛋白质折叠进行表征，高效、准确，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理。【附图说明】图1是本专利技术实施例的基于圆周统计的蛋白质折叠的表征方法的流程图；图2是本专利技术实施例的基于圆周统计的蛋白质折叠的表征系统的结构框图。【具体实施方式】为了使本专利技术的目的、技术方案及优点更加清楚明白，下面结合附图和具体实施例对本专利技术作进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。本专利技术利用圆周统计对蛋白质折叠进行表征，高效、准确，即可以更快研究蛋白质的折叠，也可以更准确描述蛋白质动态变化，从而研究相关神经退行性疾病的机理。本专利技术实施例还提供了一种基于圆周统计的蛋白质折叠的表征方法，如图1所示，包括如下步骤：步骤S101，预处理步骤：导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键热涨落、链相互作用能量热涨落及二面角相互作用热涨落的阈值；步骤S102，建立两种圆周统计模型：建立第一圆周统计模型，用角度方差V(a)=1-R(a)表征热动力学涨落；建立第二圆周统计模型，用角度协方差表征热动力学涨落；相比于线性统计学，圆周统计学更适用于解决蛋白质折叠表征的技术问题。但是圆周统计方差和协方差的数学模型在蛋白质折叠研究中还存在不一致性，具体分析如下：圆周统计中角度α可以用二维向量表示为角度方差可以用来研究蛋白质折叠过程中物理量的热动力学涨落，其公式为：v(a)=1-R(a)(E1.1)其中R2(a)=<cos(a)>2+<sin(a)>2＜1，尖括号表示平均值。圆周统计的协方差可以表示两个角度在蛋白质折叠过程中的相互影响，其数学模型表示为：由此导出的方差为：两个方差的公式E1.1和E1.3不一致，因此圆周统计方差本文档来自技高网...

【技术保护点】
一种基于圆周统计的蛋白质折叠的表征方法，其特征在于，包括如下步骤：步骤1：建立第一圆周统计模型，用角度方差V(a)=1?R(a)表征热动力学涨落；建立第二圆周统计模型，用角度协方差表征热动力学涨落；步骤2：用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析；步骤3：将待研究蛋白的能量力场分割为侧链与主链的相互作用能量Es?b、不同侧链之间的相互作用能量Es?s、不同主链之间的相互作用能量Eb?b；分析主链上第m个氨基酸残基和第n个氨基酸残基之间的侧链相互作用能量的涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第n个氨基酸残基进行标记；步骤4：分别用第一圆周统计模型和第二圆周统计模型分析不同侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。FDA0000436299540000011.jpg,FDA0000436299540000012.jpg

【技术特征摘要】
1.一种基于圆周统计的蛋白质折叠的表征方法，其特征在于，包括如下步骤：步骤1：建立第一圆周统计模型，用角度方差V(a)＝1-R(a)表征热动力学涨落，其中R2(a)＝<cos(a)>2+<sin(a)>2＜1；建立第二圆周统计模型，用角度协方差表征热动力学涨落；步骤2：用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中两侧链之间的氢键、两主链之间的氢键及两侧链之间的二硫键进行热动力学涨落分析；步骤3：将待研究蛋白质的能量力场分割为侧链与主链的相互作用能量Es-b、不同侧链之间的相互作用能量Es-s、不同主链之间的相互作用能量Eb-b；分析主链上第m个氨基酸残基和第n个氨基酸残基之间的侧链相互作用能量的涨落，当涨落值大于预设链相互作用能量热涨落阈值时，对第m个氨基酸残基和第n个氨基酸残基进行标记；步骤4：将待研究蛋白质的侧链分为亲水侧链和疏水侧链，并分别对亲水侧链和疏水侧链进行标记，分别用第一圆周统计模型和第二圆周统计模型分析不同亲水侧链之间的二面角相互作用热涨落、亲水侧链与疏水侧链之间的二面角相互作用热涨落、不同疏水侧链之间的二面角相互作用热涨落；当涨落值大于预设二面角相互作用热涨落阈值时，采用第一圆周统计模型表征；当涨落值小于预设二面角相互作用热涨落阈值时，采用第二圆周统计模型表征。2.根据权利要求1所述的基于圆周统计的蛋白质折叠的表征方法，其特征在于，在步骤1之前还包括如下步骤：导入一条待研究蛋白质序列的构型，给出氨基酸残基个数，识别主链及侧链，并分别对主链和侧链进行编号，分别设定氢键或二硫键热涨落的阈值、链相互作用能量热涨落的阈值及二面角相互作用热涨落的阈值。3.根据权利要求1所述的基于圆周统计的蛋白质折叠的表征方法，其特征在于，所述步骤4还包括：用第一圆周统计模型或第二圆周统计模型对待研究蛋白质序列构型中不同亲水侧链之间的二面角相互作用、亲水侧链与疏水侧链之间的二面角相互作用、不同疏水侧链之间的二面角相互作用进行热动力学涨落分析。4.根据权利要求2所述的基于圆周统计的蛋白质折叠的表征方法，其特征在于，所述步骤3还包括：对待研究蛋白质序列上每个氨基酸残基进行编号，分析不同氨基酸残基之间的不同侧链之间的相互作用能量...

【专利技术属性】
技术研发人员：魏彦杰，张慧玲，彭丰斌，孟金涛，冯圣中，
申请(专利权)人：深圳先进技术研究院，
类型：发明
国别省市：