[发明专利]基于随机优化的老年痴呆致病机理的分析方法

说明书

【技术领域】

本发明涉及生物信息学，特别涉及一种基于随机优化的老年痴呆致病机理的分析方法。

【背景技术】

蛋白质折叠研究蛋白质在短时间内从一级结构（亦即一维多肽链）折叠成天然三维结构，形成具有生命功能的大分子。生物体的遗传信息（DNA）通过RNA转录和翻译过程传递给蛋白质，因此蛋白质折叠也被称为第二遗传密码。通过蛋白质折叠研究可以揭示生命遗传信息的表达和功能传递的奥秘。蛋白质在从一级结构折叠到天然三维结构的过程中，可能发生误折叠或聚集，其结构和功能也受到破坏，从而引起‘折叠病’，例如老年痴呆症。

老年痴呆症又叫阿尔茨海默病，是一种中枢神经系统变性病。老年痴呆症起病隐袭，病程呈慢性进行性，是老年期痴呆常见的一种类型，主要表现为渐进性记忆障碍、认知功能障碍、人格改变及语言障碍等神经精神症状，严重影响社交、职业与生活功能。目前，美国500多万老年痴呆症患者每年的花费是1830亿美元；中国的老年痴呆症患者也有500-600万。蛋白质折叠研究有助于蛋白质分子药物的设计，因此，蛋白质折叠研究对探索老年痴呆症致病机理以及对老年痴呆症的预防和治疗将起重大的帮助作用。

对于一个有100个氨基酸的蛋白质，假设每个氨基酸有10个构象，其构象数目即有10¹⁰⁰；进一步假设搜索每个构象需要10-15秒，按照目前最快的千万亿次计算机处理速度，要搜索所有的蛋白质结构空间至少需3×10⁷⁷年。因此，搜索蛋白质全部结构空间不现实，这也对蛋白质结构高效搜索算法提出了非常高的要求。

蛋白质折叠算法依赖于一个基本假设，即蛋白质的天然结构是蛋白质自由能最低的结构。传统的分子动力学模拟和蒙特卡洛模拟方法在蛋白质折叠研究时，低温下的模拟会‘陷’在蛋白质能量的局部最优区间，很难‘跳’出来从而找到全局最优解。

经典的WangLandau算法属于随机优化领域，WangLandau算法能解决生物医学、统计物理等多个领域复杂的问题。例如，蛋白质折叠研究，该算法有两个最显著的优点：第一，该算法不会局限在局部最小能量状态，因而能搜索到全局最小能量状态；第二，该算法可模拟和计算出蛋白质系统的态密度，因而可进一步求解得到宽广温度范围内的热动力学量，如比热，实现对蛋白质系统进行定量分析。

【发明内容】

本发明要解决的技术问题在于提高随机优化的老年痴呆致病机理的高效方法精度和速度。

为此，本发明提供一种基于随机优化的老年痴呆致病机理的分析方法，包括以下步骤：

S1、确定能量模型：采用ECEPP能量力场模型以及角度坐标系；

S2、搜索最小自由能蛋白质构型：包括

主进程，

S211、局部最小化方法获得能量上界E_min和能量下界E_max，初始化蛋白质系统的态密度对数S(E)、修正因子df，df|_E＝(κΘ(E₀-E)+1)lnf，

其中，Θ(E₀-E)为Heaviside分段函数，κ、E₀、f是模型的参数；

S212、进行Metropolis式的随机变动，构建新构型，计算能量E_new、并修改E_min或E_max为E_new；

S213、重复步骤S212，完成进程间通信，计算全局S(E)；广播全局S(E)至所有从进程；

S214、执行步骤S212至步骤S213的迭代，E_min三次保持，迭代终止；

从进程，

S221、局部最小化方法获得能量上界E_min和下界E_max，初始化蛋白质系统的态密度对数S(E)、修正因子df，df|_E＝(κΘ(E₀-E)+1)lnf，

其中，Θ(E₀-E)为Heaviside分段函数，κ、E₀、f是模型的参数；

S222、进行Metropolis式的随机变动，构建新构型，计算能量E_new、并修改E_min或E_max为E_new；

S223、重复步骤S222，完成进程间通信，接收步骤S213计算的全局S(E)；