论文部分内容阅读
近年来,随着分子生物学的飞速发展,人们将越来越多的目光集中于生物大分子这一研究领域。而DNA作为一种典型的生物大分子,其所具备的储存和传递生命遗传信息的能力使得它逐渐成为了分子生物领域的一个研究热点。由于DNA几何构型的平衡与稳定直接影响着遗传信息的表达,因此,对于DNA平衡几何构型及其稳定性的研究是了解并应用DNA分子链的基础。自1953年Watson和Crick借助X射线晶体衍射技术推测出DNA的双螺旋结构以来,关于 DNA分子链的理论基础研究不断取得突破。尽管对于 DN A分子链的内部分子结构的探索属于量子力学的研究范畴,但其在生物细胞中的宏观几何构型与外力势能作用下的力学性能的研究可以借助于经典力学的弹性细杆模型来实现。由此产生了一个经典弹性力学与分子生物学相互交叉的新的研究领域。其在本质上采用连续介质力学的概念与思路研究曲杆的几何形态与变形,而在方法上又可运用非线性力学和分析力学对DNA分子链的平衡构型与稳定性进行分析。因此,本文在连续介质力学的基础上,通过建立空间弹性细杆模型,来模拟计算DNA分子链在溶液中的平衡几何构型,并对其外力作用下的弹性响应与几何构型的稳定性做了研究分析。 第一部分简要介绍了 DNA弹性细杆的几何描述,概述了 Kirchhoff弹性杆理论的基本假定与其适用的前提条件和物理意义,给出了弹性细杆的Kirchhoff平衡方程;阐述了物质的界面层与界面张力、界面自由能的概念,通过 Young-Lap lac e方程推导出了固-液界面张力的表达式,将其引入Kirchhoff方程,并沿细杆截面周长对弧坐标进行积分得到了受界面牵引力作用的DNA分子链平衡控制方程。 第二部分则基于固-液界面的吸附原理,通过Gibbs吸附方程与Langmuir吸附方程,建立起溶液浓度与固-液界面张力的关系,并推导出其具体表达式;与此同时,借助Poisson-Boltzmann理论将溶液浓度引入长分子链的熵弹性,得到了溶液浓度与 DNA分子链弹性模量的直接关系式。将上述两式代入 DNA弹性细杆的平衡控制方程,应用龙格·库塔数值算法模拟计算出以圆柱形螺旋线与椭圆柱形螺旋线为初始构型的DNA链段在不同溶度浓度中的凝聚构型,并对其进行了初步的分析。 第三部分针对DNA在外力载荷作用下的边界条件,将端部力与界面张力共同作用下的DNA链段的几何构型与弹性响应的研究归结为求解Kirchhoff微分方程组的边值问题,并选用打靶法给出在其不同端部力作用下的数值解,从而确定DNA链段的几何构型。同时,利用DNA链段两端在主轴坐标系中的坐标计算得到不同端部力对应的末端距数值,拟合出力末端距曲线,并分析了不同溶液浓度对DNA链段的力末端距曲线的影响。 第四部分从分析力学的角度出发,通过 Kirchhoff动力学比拟,将动力学的时间变量 t置换为一维空间变量即弧坐标s,从而得到了基于弧坐标的哈密顿原理与Lagrange方程。并基于此推导出能量密度函数依赖于曲率、挠率及它们一阶导数的DNA螺旋线的Euler-Lagrange方程。 第五部分则通过引入截面的扭转角及其一阶导数,将曲线的Euler-Lagrange方程推广至曲杆的Euler-Lagrange方程组。在此基础上,通过与实验数据进行了对比,分析了采用圆截面弹性细杆模型与椭圆截面弹性细杆模型分别模拟A-,B-,Z-DNA几何构型的可行性,并拟合了两种模型的r0-h曲线。同时,分析了弹性细杆截面的几何特性对螺旋带模型几何构型的影响。 第六部分通过将外力势能项引入弹性细杆的能量密度函数,在弹性曲杆模型的Euler-Lagrange方程组基础上,推导出DNA弹性细杆在端部拉伸力作用下的Euler-Lagrange方程组,并基于此对不同初始曲率与初始扭率的DNA弹性细杆模型的拉伸稳定性进行了初步分析。