颠覆性发现：聚合物中的“熵拔河”现象自发产生定向运动
在生物系统和人工材料中，含有活性片段的聚合物广泛存在，如染色质纤维、细胞骨架等。传统观点认为，要实现聚合物的定向运动，需要在模型中人为引入活性力的方向性关联。然而，这一假设面临着一个根本性挑战：生物系统中的活性力往往表现为各向同性的随机波动，如何产生定向运动一直是未解之谜。
维也纳大学Jan Smrek课题组通过构建一种活性-被动嵌段共聚物模型，发现当链段浓度足够高时，即使两个片段仅在波动幅度上存在差异（高温片段波动更强），没有任何人为引入的方向性关联，拓扑约束的存在会自发打破平移对称性，产生持久的定向运动。这一机制源于链末端熵力的不平衡，导致活性片段像“拔河”一样拉动整个聚合物链运动。该发现对理解致密活性聚合物体系的力学性质，特别是对生物功能至关重要的染色质构象和动力学具有重要意义。相关论文以“Entropic Tug of War: Topological Constraints Spontaneously Rectify the Dynamics of a Polymer with Heterogeneous Fluctuations”为题，发表在Phys. Rev. X上。并被Nature网站以“Polymers with purpose: molecules can squirm free of the pack”为题，作为研究亮点报道。
研究团队通过标度理论和分子动力学模拟揭示了这一现象的物理本质。在由不可穿透障碍物构成的网格中（图1a），聚合物链受到拓扑约束的限制。在平衡态下，链末端的熵力相互平衡；但当链的两端处于不同温度时，这种平衡被打破，产生与温度差ΔT成正比的净力Δf ~ kBΔT/a（其中a是网格尺寸）。这一力差驱动聚合物沿其原始路径定向运动。模拟结果证实，轮廓速度v与温度差ΔT呈线性关系（图1b），与网格尺寸a成反比（图1c），与理论预测完全吻合。
图1： (a) 模拟聚合物的二维示意图，其中包含处于温度Th的高温片段（橙色）和处于温度Tc的低温片段（蓝色），聚合物置于由不可穿透静态障碍物（黑色）构成的网格中，网格开口尺寸为a，配位数为z=4。(b) 和 (c) 轮廓速度（单位为每τ的网格常数）随温度差ΔT的变化关系（b，a=9时），以及随网格间距倒数1/a的变化关系（c，Th=3时）。虚线是对数据点（符号）的线性拟合。
进一步分析显示，这种定向运动对聚合物的扩散行为产生显著影响。研究发现，扩散系数D与轮廓速度v和每个网格单元的单体数Ne满足标度关系D ~ vNe b2（图2）。这一关系与从质心长时均方位移直接计算得到的扩散系数高度一致，验证了理论模型的可靠性。
图2： 扩散系数随网格间距倒数1/a的变化关系（Th=3.0时），蓝色点为由轮廓速度v b2 Ne计算得到，橙色点为直接从长时均方位移测量得到。插图：原始管道的每个网格单元包含的单体数Ne随间距a的变化关系（符号）。绿色虚线为幂律拟合结果Ne = 0.85a1.19。
在链段动力学方面，研究揭示了有趣的现象。活性片段和被动片段的中间单体的均方位移展现出复杂的标度行为（图3a、3b）。短时间尺度上表现为Rouse模型的亚扩散特征（指数0.5），而当活性片段探索多个网格单元后，被动片段被拉动并呈现超扩散行为（指数>1）。通过将时间和长度尺度按网格尺寸a进行归一化，不同系统的动力学行为呈现出普适的坍塌行为（图3c、3d），表明网格尺度是控制动力学行为的关键参数。
图3： (a) 高温（实线）和低温（虚线）片段中间单体的均方位移除以时间g1(t)/t随时间的变化，不同颜色代表不同网格间距a。(b) 均方位移的有效标度指数x1 = dln(g1)/dln(t)。点划线标记出指数1/2（Rouse模型）和1（扩散）。(c)、(d) 分别与(a)、(b)相同的数据，但按a2和t进行归一化（t定义见正文，虚线标记）。为清晰起见，未归一化的版本(a)、(b)中仅展示部分系统。
研究还发现，超扩散指数x1与链的构象指数2ν存在内在联系（图4）。当单体沿原始路径做弹道运动时，其均方位移的标度指数由路径的统计性质决定，即x1 = 2ν。这一发现解释了为何在不同参数下超扩散行为呈现多样性，并为实验观测提供了理论基础。
图4： 低温片段的指数2ν随轮廓距离s的变化关系。垂直虚线对应s（见正文）。插图：指数max(x1)与2νs的比较。a=3时的max(x1)值因采样频率较低而被低估。
质心运动的分析进一步揭示了系统的复杂动力学特性（图5a、5b）。活性片段、被动片段以及整个链的质心均方位移展现出从亚扩散到超扩散再到正常扩散的完整演化过程。值得注意的是，超扩散所能达到的最大指数x3随网格尺寸a变化（图5b插图），这与链的有限尺寸效应和横向Rouse弛豫密切相关。通过归一化处理，不同系统的质心动力学也呈现普适行为（图5c、5d）。
图5： (a) 高温片段、低温片段和整个链的质心均方位移除以时间g3(t)/t随时间的变化，不同颜色代表不同网格间距a。垂直虚线表示时间t̂，此时高温片段的g1(t̂)=2.5a2，标志着超扩散的开始。(b) 均方位移的有效标度指数x3 = dln(g3)/dln(t)。点划线标记出指数1/2（片段的Rouse动力学，总质心的蛇行运动）和1（扩散）。插图：最大指数x3随网格间距a的变化关系（符号）。拟合1使用公式(4)，拟合2直接使用v(a)的测量值（见附录B.4节）。(c)和(d)数据按a2和t*进行归一化（定义见轴标签）。为清晰起见，未归一化的版本(a)、(b)中仅展示部分系统。
研究团队进一步将理论推广到真实熔体体系。通过对柔性链和半柔性链熔体的模拟（图6），证实了上述方向性动力学行为的普遍存在。半柔性链体系的最大超扩散指数max(x3)=1.5与平均场理论预测值1.56高度吻合，而柔性链体系由于更强的约束释放效应，观测值1.22略低于平均场预测的1.47。
图6： (a) 高温（实线）和低温（虚线）片段中间单体的均方位移除以时间g1(t)/t随时间的变化，不同颜色代表不同链刚性。(b) 均方位移的有效标度指数x1 = dln(g1)/dln(t)。点划线标记出指数1/2（Rouse模型）和1（扩散）。(c) 高温片段、低温片段和整个链的质心均方位移除以时间g3(t)/t随时间的变化，不同颜色代表不同链刚性。(d) 均方位移的有效标度指数x3 = dln(g3)/dln(t)。点划线标记出指数1/2（片段的Rouse动力学，总质心的蛇行运动）和1（扩散）。
本研究不仅揭示了拓扑约束与异质波动协同产生定向运动的新机制，还为理解染色质动力学提供了重要理论框架。研究预测，活性片段的动力学特征可用于区分染色质折叠的被动机制（片段间特异性相互作用）和主动机制（不同活性水平）。未来的研究将进一步探索这一方向性运动对活性-被动共聚物熔体的稳态相行为的影响，以及其在功能性高分子材料开发中的应用前景。
热门跟贴