[发明专利]一种基于模拟退火算法的假结RNA结构预测方法

说明书

本发明涉及一种基于模拟退火算法的假结RNA结构预测方法，所述假结RNA结构预测方法涉及RNA二级结构及其自由能，所述RNA二级结构包含一种梯子结构，所述RNA二级结构的自由能能根据近邻热力学参数计算出来，随着RNA序列长度的增加，二级结构候选解的数量也呈指数量级增加。本发明能有效得到问题最优解，预测效率较高，收敛性较好，弥补了现有现有方法的缺陷。

技术领域

本发明涉及生物工程技术领域，特别涉及一种一种基于模拟退火算法的假结RNA结构预测方法。

背景技术

RNA(脱氧核糖核酸)是生物系统内最为重要的分子之一，它在生物体内具备多种功能。在最初的生物中心法则中，RNA被认为在表达遗传信息时作为蛋白质翻译，发挥了一部分短暂的作用。后来的研究发现，RNA除了翻译蛋白质，还具有多种其它功能，如调控基因表达、转运RNA、催化肽链形成和指导蛋白质合成等。随着科研工作者对RNA研究的逐步深入，RNA带给人们的形象也随之发生了变化。RNA被认为不仅仅是DNA到蛋白质之间的一种信息传递中介，它逐渐由功能单一的简单线性碱基序列，演变成种类多样、结构复杂、功能特异的生命核心物，同时RNA在中心法则中，取得了与DNA和蛋白质同样重要的地位。在RNA结构预测算法的研究上，近些年来虽然衍生了许多新的算法，但是从根本上而言，还是在原先经典算法的基础上加以改进的结果。因此从预测效率的角度来看其时间复杂度和空间复杂度并未取得大的改进。目前RNA二级结构预测的算法主要有比较序列分析和热力学最小自由能两种，由于比较序列分析法需要知道多个与之同源的序列及相应的二级结构信息因此相对于后者而言发展缓慢。当仅有一个RNA序列或者序列为差异很小的RNA小家族时，比较序列分析方法根本不适用。而后一种方法依据碱基配对和自由能分配规则，搜索有最小自由能量的二级结构。其基本假设是具有最小自由能量的RNA二级结构是最稳定的，并且最接近真实的RNA生物分子结构。以基于热力学的最小自由能算法是近年来RNA二级结构预测算法的代表，并在此算法基础之上得到了很大的改进和演进，出现很多其他的算法。

由于生物信息学领域新的数据呈现指数增长，海量的序列数据信息使得最小自由能算法的发展极其迅速。由于RNA分子具有降解速度快，难以结晶等特点，故通过x射线晶体衍射和NMR等实验方法去测定RNA分子的立体结构很不容易。计算机技术的快速发展使得各种复杂的预测算法在计算机上得以实现由此带动各种RNA预测算法的不断涌现。虽然现在各种各样的RNA算法层出不穷，但是大多数是在原先经典算法的基础之上加以改进以求计算精度和时间复杂度等，没有真正解决问题的难点。如，对于RNA二级结构预测问题有核磁共振、图形化暴力破解等传统物理方法可以较为准确的得到RNA结构，然而其结构并不稳定，使得这些物理方法都会消耗巨大的资源和时间且得到的结果也未必是准确的结构。但是通过RNA一级结构去预测其二级结构的过程是相当复杂的，其一，随着RNA分子长度的增加，其预测难度将呈指数级提升，较长的RNA分子十分难去预测；其二，在真实的RNA结构中存在假结，这种假结结构预测准确度更低。当下RNA二级结构主流的预测算法中，基于序列比较分析的算法需要大量的先验已知序列信息，如果信息不足其预测准确度不高；基于动态规划的算法计算成本较高；基于启发式算法，在有限的资源内，通过种群的寻优策略给出一个近似解决方案作为预测结果来预测RNA二级结构，但是这种方式得到最小自由能的近似值往往会导致预测结果录入的是局部最优的结果。RNA二级结构中茎区有时会形成一种特殊的位置关系，称之为假结(pseudoknot)。假结作为一种特殊的结构，是实现结构功能的重要因素。但是因为碱基相互重叠的特性，生物信息计算对假结的预测难以奏效。预测包含假结的RNA二级结构问题是NP难的，但要分析RNA的真实结构，对假结的准确预测是必须解决的问题。假结结构对RNA分子起着十分重要的作用，假结的相关研究实验也变成了RNA二级结构预测中的重中之重。