受限空间中单链高分子动力学行为的计算机模拟

摘要

近年来，高分子链在受限环境中的动力学行为引起了人们极大的关注。对这一动力学行为的研究一方面有助于我们更好地理解细胞生物学中的很多现象；另一方面则因为其有着很广泛而又重要的科技应用前景。高分子链在受限环境中的动力学行为是一个极其复杂的过程，因而限制了很多实验手段的应用。另一方面，理论计算在处理高分子链单元的体积排斥效应时显得力不从心。由此，“计算机模拟”成为此领域中一种不可替代的重要研究手段。随着计算机硬件的发展和新的模拟技术和模型的提出，计算机模拟将会发挥更大的作用。
　　 相比连续空间的全原子模型，粗粒化格子模型虽然忽略了部分结构细节，但是从计算机模拟速度看更为高效，能够以不太长的计算机机时获得想要的结果。格子模型特别适合于进行Monte Carlo模拟。通过选择恰当的算法，MonteCarlo模拟还能对体系的动力学过程进行研究。
　　 本论文的主体部分采用格子空间的动态Monte Carlo模拟对单链高分子在受限空间中的动力学行为进行了研究。此外，为了便于与具有体积排斥作用(self-avoiding effect)的自规避行走(SAW)链所得到的结果进行比较，本文还对处在二维方形受限空间中的理想高分子链(高斯链)穿越纳米管道的动力学行为进行了理论研究。
　　 本博士论文的主要创新性贡献可以分为以下几个部分：
　　 1、理论研究了处在二维方形受限空间中的理想单链高分子穿越纳米管道进入另一尺寸较大的方形空间的动力学过程。计算了高分子链穿越过程中各状态下的自由能以及所需的穿越时间。链单元与管道之间的相互作用明显改变了穿越所需克服的势垒。管道长度增加时，穿越由受熵势垒控制转变为由链单元-管道相互作用控制。这种转变导致穿越时间随链长非单调变化。我们计算得到了穿越时间与高分子链长、管道长度以及链单元-管道相互作用的精确表达式。研究结果有助于人们理解各因素如何影响高分子链的穿越速度。
　　 2、使用Monte Carlo方法模拟研究了二维方形受限空间中的单链高分子穿越狭窄管道的动力学过程。模拟得到了高分子链在穿越过程所需克服的自由能势垒大小及各步所需时间随链长N以及管道长度M的变化关系。熵势垒的大小线性依赖于链长及管道长度的变化；链长增加时，第一和第二步穿越所需时间τ1和τ2线性增加，而第三步所需时间T3保持不变；如果保持链长不变，τ1～M2.62，τ2随M线性增加，τ3～M1.90。我们将扩散速度考虑进去，采用Fokker-Planck方程解析计算了穿越各步所需时间随N和M的依赖关系，所得结果与模拟结果很好地吻合。
　　 3、在前面工作的基础上，进一步研究了链单元-管道相互作用能ε对处在受限方形空间中的高分子链穿越纳米管道动力学行为的影响。模拟结果揭示，高分子链的整个穿越过程可分为两个阶段，即高分子链的受限阶段和成功穿越阶段，所需时间分别定义为τtrap以及τeac。log(τtrap)随链长N以及管道长度M的增加线性变化。对于前者，当链单元-管道相互作用能ε增加时，线性关系的斜率保持不变；而对于后者，线性关系的斜率随着ε的增加而变大。链单元-管道之间存在较强的排斥或吸附作用时，高分子链单元很难进入或离开管道，导致高分子链的穿越时间τesc迅速增大。
　　 4、采用理论计算和Monte Carlo模拟研究了单链高分子从一温度较低(T1)的半无限大二维空间穿过小孔进入另一温度较高(T2)的半无限大空间的动力学过程。两种方法所得结果吻合的很好。结果表明，当两侧空间存在温度差异时，正的空间化学势有利于高分子链从低温区域进入高温区域：高分子链穿越小孔所需时间τ与化学势△μ成反比，而与链长N成正比；温差△T(=T2-T1)越大，τ越小。
　　 5、采用“键涨落”模型(BFM)，对电场作用下单链DNA分子穿越深度不同交替排列的纳米管道的动力学行为进行了模拟研究。三维空间中，随着电场强度ε的增加，DNA链的迁移率逐渐增加并渐趋一饱和值；场强ε相同时，长链DNA的迁移率大于短链分子；DNA链进入窄管道前越过熵势垒所经历的等待时间即熵受限时间τtrap指数依赖于外加场强的倒数。模拟结果与实验结果很好地吻合。对高分子链在熵受限阶段的构象变化的细致研究揭示了装备对不同长度的高分子链的分离机制。DNA链在接近宽窄管道交界处时在电场作用下沿界面方向伸展变形，与短链分子相比，长链DNA有更多的链单元有机会接触界面，从而提高了越过熵势垒的几率，带动整条链进入窄管道，因此熵受限时间较短，平均迁移率较大。对DNA链在二维管道中的穿越过程的模拟研究所得结果证明了我们对上述分离机制的分析的正确性。
　　 6、将Monte Carlo方法和“键涨落”模型(BFM)相结合，模拟了受限于两平行平面间的单链DNA分子在外力作用下的伸展以及撤除外力后弛豫的动力学过程，研究了受限程度对DNA分子的伸展长度及弛豫过程的影响。模拟结果表明，随着受限程度的增加，DNA分子链的构象更加伸展，这主要是由于随着平面间距的减少，DNA分子不同链段之间流体力学相互作用将会被平面屏蔽所致，受限程度不同时DNA分子的弛豫过程进一步证实了这一点。DNA分子的伸展长度(即末端距)随着流速的变化关系与实验结果及其对此所做的理论分析是基本一致的。
　　 7、在电场作用下，DNA链的一部分链段进入受限程度较高的低熵区域而其余链段仍滞留在高熵的自由区域。电场撤除后，进入受限区域的链段将会自发回弹到自由区域。我们采用Monte Carlo模拟研究了这一动力学过程。模拟结果揭示，DNA链的回弹速度由慢到快逐渐增加，这与文献给出的实验现象及建立的理论模型是一致的，DNA链的回弹应归因于链单元在两个区域不同受限程度所产生的熵力的推动。模拟还研究了形成熵受限区域的障碍物的尺寸及间距对DNA链回弹速度的影响，进一步的分析表明链单元在受限区域的熵正比于链单元所能到达的自由空间占整个受限区域的比例。

目录

文摘

英文文摘

声明

第1章绪论

1.1计算机模拟

1.1.1 Monte Carlo方法的产生及发展

1.1.2 Monte Carlo方法在高分子科学领域的应用

1.1.3 Monte Carlo方法的基本思想及统计理论基础

1.1.4 高分子链的模型

1.1.5 高分子链的抽样方法

1.1.6键涨落模型(bond fluctuation model)

1.2外界势场下的理想链

1.3单链高分子在受限环境下的动力学行为

1.3.1 高分子链穿越纳米管道的动力学行为

1.3.2 高分子链在受限环境中的伸展

1.3.3 受限环境中的高分子链的自发回弹

1.4论文的主要内容

参考文献

第2章理想高分子链穿越纳米管道的解析计算

2.1引言

2.2 理论与算法

2.2.1 模型的建立

2.2.2 自由能的计算

2.2.3 平均穿越时间的计算

2.3结果与讨论

2.4结论

参考文献

第3章受限SAW链穿越纳米管道的计算机模拟

3.1引言

3.2模型和算法

3.3结果与讨论

3.3.1 穿越过程的自由能势垒

3.3.2 穿越时间与链长N的关系

3.3.3 穿越时间随管道长度M的变化

3.4结论

参考文献

第4章链-管道相互作用对受限SAW链穿越纳米管道的影响

4.1 引言

4.2模型和算法

4.3结果与讨论

4.3.1 链单元-管道相互作用能ε对Ttrap的影响

4.3.2 链单元-管道相互作用能ε对Tesc的影响

4.3.3 高分子链穿越的细致动力学过程的模拟

4.4结论

参考文献

第5章 温度差异对链穿越纳米管道动力学行为的影响

5.1 引言

5.2模型和算法

5.3 理论计算

5.4模拟结果与讨论

5.4.1 空间化学势μ对穿越时间T的影响

5.4.2 穿越时间T对链长N的依赖关系

5.4.3温度的影响

5.5结论

参考文献

第6章电场驱动下DNA链穿越熵受限管道的研究

6.1引言

6.2模型和算法

6.3理论准备

6.4结果与讨论

6.4.1 三维空间模拟

6.4.2 二维空间模拟

6.5结论

参考文献

第7章单链DNA在受限环境中伸展的计算机模拟

7.1引言

7.2模型和算法

7.3结果与讨论

7.3.1 无束缚时DNA链在外电场作用下的迁移

7.3.2 受限空间中一端固定的DNA链在外电场作用下的伸展

7.3.3 电场撤除后DNA链的弛豫

7.4结论

参考文献

第8章受限高分子链自发回弹的计算机模拟

8.1引言

8.2模型和算法

8.3结果与讨论

8.3.1 障碍物尺寸Dpillar对DNA链回弹的影响

8.3.2 障碍物间距Spillar对回弹的影响

8.4结论

参考文献

致谢

发表的学术论文