分子内非互易性编程活性分子动力学

《Newton》：Programming active-molecule dynamics via intramolecular nonreciprocity

　　

专业术语翻译摘要 活性物质（active matter）由消耗能量以运动的自主驱动单元组成，从微生物到合成自推进颗粒均属此类。通常，单个活性颗粒的运动方式由其构造方式所固定。本研究表明，通过将少量活性颗粒连接成类分子链，可以产生更广泛的运动模式。链中含有A和B两类颗粒，其取向以非互易方式（nonreciprocally）耦合，即A对B的影响方式不同于B对A的影响方式。因此，仅改变链内A和B的排列顺序即可改变整个分子的运动方式，而无需改变颗粒本身或颗粒间的相互作用力。不同的序列可产生持续运动、圆周运动或蠕动运动（wriggling motion），非互易耦合强度（nonreciprocal coupling strength）提供了额外的调控手段来调节这些行为。单分子层面的差异也会影响多个分子的集体组织方式。本研究以运动诱导相分离（motility-induced phase separation, MIPS）——一种由持续运动而非吸引力引起的自发聚集现象——为例。即使仅由三个颗粒组成的链，改变内部序列也会使发生聚集所需的推进强度（propulsion strength）改变超过三个数量级。分析进一步表明，这种变化与每种分子在失去运动方向之前能够行进的距离有关。这些结果表明，分子的内部序列可以作为编程活性物质的简单“代码”，将单分子运动与大规模集体组织联系起来。

论文解读

一、研究背景与问题

活性物质是由持续消耗能量以驱动自身运动的单元组成的非平衡系统，涵盖从微生物到合成自推进颗粒的广泛范畴。传统上，单个活性颗粒的动力学行为由其微观构造所“硬编码”决定，包括其平移、旋转及对外界相互作用的响应方式。这意味着若要获得新的动力学行为，通常需要重新设计颗粒本身，这一限制严重制约了活性物质系统行为多样性的拓展。一个基本问题由此浮现：能否超越单颗粒设计规则，构建简单、可编程的结构，使其产生单个活性颗粒本质上无法实现的动力学响应？

近年来，非互易相互作用（nonreciprocal interactions）成为丰富活性物质动力学的新兴途径，即一个颗粒对另一个颗粒施加的影响并不被对称地返还。此类相互作用已通过化学信号、流体动力学耦合或外部施加的各向异性在实验上实现，并展现出手性运动、振荡和非互易相变等有趣的集体现象。然而，已有研究主要关注不同活性单元之间的非互易性如何塑造多体系统的集体行为，而非互易性能否被嵌入到具有内部有序结构的类分子构型中，并作为可编程的内部驱动力加以利用，仍基本未被探索。特别是，分子内非互易性与可调内部序列相结合，能否在单分子层面产生不同的动力学模式并系统性地影响集体组织，此前尚不清楚。

二、研究内容与核心结论

针对上述问题，研究人员开发了一个最小化框架，将内部序列作为主要设计变量，核心思想是在序列层面实现可编程性，而不改变任何颗粒层面的相互作用。研究人员构建了由两种自推进颗粒（A和B）按指定序列排列而成的短链活性分子，颗粒间通过非互易取向耦合（nonreciprocal orientational coupling）相互作用。同种相邻颗粒间为互易耦合，而异种相邻颗粒间为非互易耦合：一种颗粒倾向于与另一种颗粒对齐，而反向耦合则倾向于不同的取向响应，从而产生序列依赖的内部驱动力。

这一架构在两个层面产生了定性全新的动力学。在单分子层面，均质序列（如[A-A-A]）表现出标准的持续随机游走运动，而异质序列则产生独特的、序列特异的运动轨迹，这些轨迹是任一组成物种单独存在时无法实现的。例如，[A-B]二聚体表现出显著的圆周运动，[B-A-B]三聚体则表现出前向蠕动运动。在集体层面，以运动诱导相分离（MIPS）为代表，研究发现仅改变内部序列即可使MIPS发生的推进强度阈值移动多个数量级，而颗粒间相互作用保持不变。这些结果确立了分子内非互易性与序列可编程性相结合，作为跨尺度编码和控制活性物质动力学的最小化机制，为可编程活性物质提供了通用设计原则。

三、关键技术方法

研究人员采用基于HOOMD-blue软件包实现的过阻尼布朗动力学（overdamped Brownian dynamics）模拟方法。模型采用过阻尼Langevin方程描述颗粒的平动和转动动力学，颗粒间通过谐振弹簧力维持分子内键合，所有颗粒对间通过Weeks-Chandler-Andersen（WCA）排斥势描述体积排除效应。非互易取向耦合编码在耦合矩阵中，其中同种颗粒间耦合强度为1，而异种颗粒间耦合强度由参数g控制。模拟中设置σ=γ=1，固定J_AA=1，旋转扩散系数D_r=10^?3，弹簧刚度k=60F₀，WCA能量尺度ε=0.5F₀，积分时间步长dt=3×10^?4/F₀。MIPS模拟中，三颗粒序列系统使用N=14,400个分子（共43,200个颗粒），单体系统使用N=43,200个分子，颗粒数密度固定为ρ=0.6。研究还通过净推进力自相关函数分析分子运动的有效参数，并利用最大团簇分数作为序参量定量表征相分离程度。

四、研究结果

模型与单分子动力学。 研究人员建立了活性分子模型，其中非互易相互作用仅作用于分子内部。在确定性极限（无噪声）下，[B-A-B]三聚体的取向动力学随非互易参数g呈现多个不同区域：g0时互易对齐占主导，分子弛豫到完全对齐的静态态；?2< />₀=arccos(1+g)的解析表达式，与模拟结果在?2< />

序列依赖的持久性运动。 在有限噪声条件下，研究人员通过均方位移（MSD）和净推进力自相关函数表征分子运动。对于[A-A-A]、[A-B-B]和[A-A-B]序列，净推进力自相关呈指数衰减形式C_Fnet(t)?F_eff²exp(?D_efft)，拟合得到有效参数分别为：[A-A-A]的F_eff=3F₀、D_eff=D_r/3；[A-B-B]的F_eff=3F₀、D_eff=3D_r；[A-A-B]的F_eff=F₀、D_eff=3D_r。利用这些参数代入解析MSD公式，与模拟结果吻合良好。相比之下，[B-A-B]序列的自相关函数围绕零值呈现显著振荡，反映其潜在的振荡取向动力学。这些结果表明，不同分子序列产生明显不同的力自相关和方向记忆水平，自相关衰减越慢意味着净推进方向寿命越长、持久长度越大。

序列依赖的MIPS阈值。 研究人员以最大团簇分数α?作为序参量，研究不同分子序列的MIPS行为。对于单体[A]分子，α?随推进强度F₀增加从接近零急剧上升至约0.9的最大值，随后在更大F₀下逐渐下降。代表性快照显示系统在小F₀时近乎均匀，中等F₀时强烈相分离，大F₀时大团簇碎裂导致α?降低。以α_th=0.6为MIPS起始的操作性阈值，四种三聚体序列的阈值推进强度F_0,th满足严格排序：F_0,th^[A-A-A]< />_0,th^[A-B-B]< />_0,th^[A-A-B]< />_0,th^[B-A-B]，跨越超过三个数量级。这一排序与单分子动力学一致：[A-A-A]具有最长持久性和最小有效旋转扩散，因此阈值最低；[B-A-B]因快速去相关强烈抑制持久性，需要最大的推进强度才能诱导MIPS。

持久距离与MIPS阈值的关联。 研究人员进一步定义了特征持久距离?_m=√MSD_F0=1(t_c)，其中t_c从净推进力自相关中提取，并定义阈值长度?_th=F_0,th?_m。尽管F_0,th在不同分子代码间变化超过三个数量级，对应的?_th值保持在同一数量级，包括对于无法定义常规持久长度的振荡性[B-A-B]序列也是如此。这一分析表明，MIPS阈值的大幅序列依赖变化可部分归因于分子持久距离的差异，[B-A-B]的大F_0,th与其更短的持久距离一致，而非意味着定性不同的聚集机制。

五、讨论与结论

研究人员研究了由两种自推进物种的短链组装而成的活性分子，其推进方向按指定内部序列非互易耦合。单分子层面，均质序列呈现标准持续随机游走动力学，而异质序列产生任一组成物种单独无法实现的序列特异轨迹。集体层面，序列重编程使MIPS起始阈值在推进强度上移动多个数量级，同时保持颗粒间相互作用不变。研究还扩展分析了二颗粒和四颗粒分子，提供了二颗粒编码的完整分类和所有对称不等价四颗粒编码的代表性轨迹。

实验上，预测的序列依赖MIPS起始偏移最可能在可编程活性机器人系统中直接检验，该系统允许直接实现指定的非互易对齐规则；也可在混合胶体微泳者组装或通过化学或泳动介导非互易耦合的活性颗粒系统中探索。更广泛地，相同的编码策略可扩展到活性物质中的其他集体现象。自然的下一个步骤是引入额外的分子内自由度（如内部偶极），以介导分子间的对齐和角动量交换，并探索由此产生的涌现相。总体而言，研究结果确立了序列编码的分子内非互易性作为跨尺度编程活性物质动力学的最小设计原则，从单分子运动到集体组织均适用。