Montesi Alberto, Pasquali Matteo, MacKintosh F C
Department of Chemical Engineering, Rice University, 6100 Main Street, Houston, Texas 77005, USA.
Phys Rev E Stat Nonlin Soft Matter Phys. 2004 Feb;69(2 Pt 1):021916. doi: 10.1103/PhysRevE.69.021916. Epub 2004 Feb 27.
We investigate the dynamics and pathways of the collapse of a single, semiflexible polymer in a poor solvent via three-dimensional Brownian Dynamics simulations. An example of this phenomenon is DNA condensation induced by multivalent cations. Earlier work indicates that the collapse of semiflexible polymers generically proceeds via a cascade through metastable racquet-shaped, long-lived intermediates towards the stable torus state. We investigate the rate of decay of uncollapsed states, analyze the preferential pathways of condensation, and describe the likelihood and lifespan of the different metastable states. The simulations are performed with a bead-stiff spring model with excluded volume interaction, bending stiffness, and exponentially decaying attractive potential. The semiflexible chain collapse is studied as a function of the three relevant length scales of the phenomenon, i.e., the total chain length L, the persistence length L(p), and the condensation length L(0)=square root of [k(B)TL(p)/u(0)], where u(0) is a measure of the attractive potential per unit length. Two dimensionless ratios, L/L(p) and L(0)/L(p), suffice to describe the dimensionless decay rate of uncollapsed states, which appears to scale as (L/L(p))(1/3)(L(0)/L(p)). The condensation sequence is described in terms of the time series of the well separated energy levels associated with each metastable collapsed state. The collapsed states are described quantitatively through the spatial correlation of tangent vectors along the chain. We also compare the results obtained with a locally inextensible bead-rod chain and with a phantom bead-spring chain. Finally, we show preliminary results on how steady shear flow influences the kinetics of collapse.
我们通过三维布朗动力学模拟研究了在不良溶剂中单个半柔性聚合物的塌缩动力学和途径。这种现象的一个例子是多价阳离子诱导的DNA凝聚。早期工作表明,半柔性聚合物的塌缩通常通过一系列过程进行,即通过亚稳态的球拍状、长寿命中间体向稳定的环形状态转变。我们研究了未塌缩状态的衰减速率,分析了凝聚的优先途径,并描述了不同亚稳态的可能性和寿命。模拟是使用具有排除体积相互作用、弯曲刚度和指数衰减吸引势的珠子-刚性弹簧模型进行的。作为该现象的三个相关长度尺度的函数,即总链长L、持久长度L(p)和凝聚长度L(0)=[k(B)TL(p)/u(0)]的平方根,其中u(0)是每单位长度吸引势的量度,研究了半柔性链的塌缩。两个无量纲比,L/L(p)和L(0)/L(p),足以描述未塌缩状态的无量纲衰减速率,其似乎按(L/L(p))(1/3)(L(0)/L(p))缩放。凝聚序列是根据与每个亚稳态塌缩状态相关的分离良好的能级的时间序列来描述的。塌缩状态通过沿链的切向量的空间相关性进行定量描述。我们还比较了使用局部不可伸长的珠子-杆链和虚珠子-弹簧链获得的结果。最后,我们展示了关于稳定剪切流如何影响塌缩动力学的初步结果。