论文部分内容阅读
蛋白质折叠热力学假说中明确指出蛋白质一级结构可以完全决定其空间构型,天然结构下的蛋白质的自由能是全局最小值.因而如何从一级结构预测其三维构型以及如何找到蛋白质折叠能量函数的最小值就成为目前国内外的研究热点.蛋白质折叠过程是热力学过程和动力学过程的综合,其折叠的随机性、过程的复杂性可想而知.因而要对折叠过程进行研究,寻找能量函数的极小值必须借助简化模型和特殊方法进行.基于这种思想本文主要研究了以下内容:一、针对目前所出现的HP格点模型均为串行模型的特点,将元胞自动机并行演化的思想应用到HP格点模型中,提出了并行HP格点模型.这种模型不仅便于对生命现象进行更真实的模拟,而且也克服了串行模型运行速度慢、计算量大等的不足.二、理想的能量函数应该能够精确表达蛋白质的所有原子空间位置及其能量之间的关系.但是由于蛋白质分子本身的复杂性,直接综合其生物、物理、化学等各方面性质建立的模型自由度非常大,求解问题遇到很多困难.因此,求解时必须对能量函数进行适当的简化.基于这种目的,本文将已有能量函数与HP格子模型相结合,构建了一种简化的能量函数,以便于运用数学方法进行求解.三、鉴于区域分解算法在高维、不规则区域的优势,我们运用区域分解算法来寻求能量函数的极小值.即将能量表面分成小子块,并行地在每一子块上搜索极小值,然后相互交换信息.如果最小值恰恰在小子域的边界上,则还需要对边界进行预测.文中主要研究了预处理迭代的一些性质,并提出了局部正交化、近似LU分解等一些构造预处理的新方法.在构造预处理的过程中还发现,运用参数来讨论条件数的改善情况,其证明过程非常简便,文中还对参数的范围作了讨论.进行边界预测时本文采用虚拟边界预测法.实例分析结果证明了这种方法的可行性.