[发明专利]基于混合并行方式的蛋白质热力学分析高效随机模拟方法

说明书

【技术领域】

本发明涉及生物信息分析技术领域，特别是涉及一种基于混合并行方式的蛋白质热力学分析高效随机模拟方法。

【背景技术】

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

目前，研究蛋白质折叠的算法大多数都在分子动力学模拟和随机模拟中实现。一般而言，分子动力学模拟常用于研究蛋白质系统的动力学过程；而随机模拟则可以研究蛋白质系统的整个热动力学过程。针对使用高精确度的全原子蛋白质模型的模拟，需要计算成千上万个原子之间的多种相互作用力，对于分子动力学模拟只能模拟纳秒级的蛋白质折叠过程，故其在微秒到毫秒时间内的蛋白质折叠研究中具有很大的局限性；此外，分子动力学模拟也受一个初始实验构型的影响。而随机模拟不但能用于微秒到毫秒时间内的蛋白质折叠研究，而且不依赖于一个具体的初始构型，可以更广泛地搜索构型空间。

经典的WangLandau算法就是随机模拟领域最吸引人最有发展情景的新算法，它能解决生物信息学、统计物理学等多个领域的很多复杂问题。比如在蛋白质折叠研究中，该算法有两个最显著的优点：第一，蛋白质模拟不会局限在局部最小能量状态，因而能较好地在整个能量区间进行自由行走；第二，通过该算法可模拟和计算出蛋白质系统态密度，因而就能进一步求解得到宽广温度范围内的很多热动力学量如比热等，这样就能高效地分析和研究蛋白质折叠的整个热力学过程。但WangLandau算法在计算精度和速度上还有待于进一步提升。

鉴于此，克服该现有技术所存在的缺陷是本技术领域亟待解决的问题。

【发明内容】

本发明要解决的技术问题是提供一种基于混合并行方式的蛋白质热力学分析高效随机模拟方法。

本发明采用如下技术方案：

一种基于混合并行方式的蛋白质热力学分析高效随机模拟方法，包括：

步骤A：确定蛋白质能量模型和能量区间；

步骤B：确定所述蛋白质能量区间的分段方式；

步骤C：模拟及计算蛋白质系统态密度。

进一步地，所述步骤A进一步包括：

采用ECEPP蛋白质能量模型，ECEPP能量力场的表达形式为：

E_ECEPP=E_C+E_LJ+E_HB+E_Tor

其中，是两电荷之间的库伦作用力，r_ij表示原子i和j之间的距离；是两原子之间的兰纳-琼斯作用力；是氢键作用力；E_Tor=∑_lU_l(1±cos(n_lξ_l))是两面角旋转作用力，ξ_l是第l个两面角。

进一步地，所述步骤A进一步包括：

对所使用的蛋白质能量区间进行离散化处理，若取k个能量bin区间值，则对[E_min,E_max]平均划分k个bin区间，用每个bin区间中间的一个能量值代表能量区间值。

进一步地，所述步骤B进一步包括：

步骤B1：对能量区间平均分为M段，设相邻子能量区间之间的重合度等于Δ个bin区间，则每一段含有个bin区间；

步骤B2：依照当前计算得到的蛋白质系统态密度函数的对数S(E)分布特点，自适应地对能量区间分段，若某个子能量区间为[E_begin,E_end]，则

进一步地，所述步骤C进一步包括：

通过MPI的主从进程模式和OpenMP的多线程并行模式，模拟及计算蛋白质系统态密度。

进一步地，在所述主从进程模式的N个分进程中，分进程1为主进程，其余分进程均为子进程。

进一步地，所述主进程包括如下步骤：

步骤S11：初始化蛋白质系统态密度函数的对数S(E)=lng(E)=0，直方图H(E)=0（E_min≤E≤E_max），修正因子df=1(=lnf=lne)；