扭转DNA孤子的能量正交性与构象恢复动力学

《Physica A: Statistical Mechanics and its Applications》：Energy orthogonality and conformational recovery dynamics of torsional DNA solitons

　　

摘要

DNA在局部机械应力下的结构完整性是基因组保真度的基础，然而支配其构象恢复的精确非线性机制仍知之甚少。在本研究中，我们利用受扰的正弦-Gordon方程，将双螺旋建模为受局部时变外力矩作用的非线性拓扑晶格。通过在扩展空间域上使用二阶St?rmer-Verlet有限差分格式进行积分，我们表征了外力停止后结构扰动场??(??,??)=??pert???unpert的弛豫动力学。扰动后的恢复阶段由复杂的、多频率的振铃包络所支配，我们通过全局??2范数‖??(??,??)‖2的双模多谐波阻尼优化拟合对其进行解析，追踪扭结内部摆动模态的激发。时空等高线分析证实了稳健的辐射泄漏机制，其中局域化势能转化为高频色散声子，沿由线性波速（??=1.0）定义的因果锥，系统地将结构应力排入远场。我们的结果表明，离散的局域化拓扑束缚态与连续辐射流形之间存在渐近动态解耦，防止了破坏性反向耦合进入核心基序。此外，触发这些结构转变所需的局域力矩阈值直接映射到能量超传输的非线性物理机制中，其中驱动频率在禁带光谱区域内过滤掉低振幅热噪声。最后，参数敏感性分析表明，结构恢复时间尺度根本上受声速传播限制约束，遵循相对于晶格扭转波速的逆幂律标度关系（??∝?????，??≈1.25）。因此，序列依赖性区域中升高的骨架刚度和碱基堆积刚度（指示高GC含量片段）表现出代数加速的稳定窗口。最终，该方法将基因组确立为一个高度优化的非平衡拓扑机器，能够干净地耗散机械冲击，同时稳健地保持结构信息。

引言

脱氧核糖核酸（DNA）是遗传信息的主要储存库，具有复杂的结构布局，促进遗传密码的保存、传递和复制[1]、[2]。由核苷酸构建模块组成，包含磷酸基团、糖环和含氮碱基，双螺旋通过精密的氢键网络和碱基堆积相互作用稳定[3]。DNA并非静态档案，而是一种高度动态的生物大分子，其结构组分不断发生构象变化，包括局部弯曲、拉伸和扭转，这些变化对复制、转录和蛋白质-DNA相互作用等关键生物功能至关重要[4]、[5]。

这种动态生物大分子运行于高度随机的非平衡细胞环境中。细胞内涨落、分子拥挤和环境热涨落引入了持续的机械扰动，威胁到局部碱基配对的稳定性[6]。因此，结构基序必须具备内在的物理驱动机制，以在持续的环境变化中维持构象稳定性。为了准确建模和解码这些非平衡生物趋势和稳定性边界，现代计算框架越来越多地利用预测性、物理约束的算法来绘制大分子韧性谱[7]。

底层的螺旋几何形状对这些参数高度敏感，结构变化可诱导局部缺陷或构象转变[8]。此外，DNA与外源性化合物（如小分子药物、嵌入剂和致癌物）的相互作用可扭曲螺旋规整性，显著改变细胞代谢和遗传稳定性[9]、[10]、[11]。在这些复杂的动力学中，扭转形变（扭转和反转）尤为重要，因为它们直接调节DNA开放态可及性和基因表达保真度[12]、[13]。

DNA中局部非线性激发的传播长期以来一直通过非线性偏微分方程（PDE）进行建模，典型示例包括Davydov孤子模型、Toda晶格的变体[14]和Peyrard-Bishop-Dauxois（PBD）模型[15]、[16]。这些模型有效地捕捉了大分子晶格中色散与非线性之间固有的关键平衡。当这些底层系统受到内源或外源扰动时，其数学处理需要对绝热演化和辐射阻尼进行严格理解，正如Kivshar和Malomed在近可积系统孤子动力学基础框架中所确立的那样[17]。

DNA通常作为活性分子马达的结构基质，其离散步进动力学产生局部非高斯机械力[18]。理解双螺旋如何从这些酶诱导的集中机械扰动中动态弛豫，对于量化全局基因组信息保真度至关重要。单分子实验和原子力显微镜已证实，DNA能够支持扭转波的传播，并在外部机械刺激下响应局部解旋[19]、[20]、[21]。

在数学上，Lie群方法已被用于验证DNA中的纵向位移表现为稳定的周期波[22]。最近，有理正弦-Gordon展开（RSGE）方法和改进辅助方程被用于推导精确的扭结孤子解，为耦合振子底层动力学提供了深刻见解[23]、[24]。在理想传播构型之外，近期研究已扩展到复杂大分子几何结构力学，特别关注支配DNA超螺旋在持续力矩下非线性动力学和拓扑转变的机制[25]。

正弦-Gordon（sG）方程已成为建模这些现象的稳健数学框架，从宏观粒子链[26]到工程机械超材料[27]、[28]都取得了有益的成果。虽然最近使用Boussinesq或Davey-Stewartson-Fokas（DSF）模型的研究探讨了时变系数如何改变孤子轮廓[24]、[29]，但类似策略也被应用于液晶中拓扑保护的缺陷，以理解对外部场的非互易响应[30]。

在这些受扰环境中，维持或恢复稳定的孤子构型需要精确干预。比例反馈和优化控制架构已成功用于稳定受扰sG扭结孤子，使其抵御外部变化[31]。此外，确保这些结构的运行寿命需要对轨道力学进行深入评估，特别是在持续振荡驱动场作用下受扰sG方程中周期轨道的稳定性[32]。

尽管已有大量关于孤立波精确解析解的文献[33]、[34]、[35]、[36]、[37]、[38]，但这些传统处理方法一般假设理想、可积或稳态情景。理解局域DNA扭结在瞬时、空间局域化机械扰动后非平衡结构恢复动力学的深刻物理鸿沟仍然存在。当前数学模型着重于静态结构轮廓或无摩擦传播，完全绕过了分子的时间弛豫历史。具体而言，活性外部驱动力（模拟酶力矩）与随后该势能不可逆地转化为辐射连续阻尼之间的非稳态相互作用尚未在多元分量能量空间中定量处理。这一遗漏限制了我们对其机械稳定性的理解。

本工作突破传统稳态分析，研究受扰sG缺陷空间在显式时变非齐次力矩作用下的结构弛豫和模态耦合。本研究的核心目标在于：

1. 通过实施双模多谐波阻尼优化，系统地映射瞬态弛豫阶段，提取特征结构恢复时间（??），评估局域DNA基序在局域应力下的谱鲁棒性。
2. 追踪系统的底层Hamilton空间，提供离散局域化拓扑束缚态（扭结核心）与连续辐射流形之间的渐近动态解耦。我们展示了基质如何沿严格的因果锥（??=1.0）将多余的局域化能量蜕入远场声学声子。
3. 建立序列依赖性标度定律，证明弛豫时间遵循相对于扭转波速的逆幂律标度关系（??∝?????）。这明确地将抽象孤子扰动参数与碱基对特异性扭转刚度（??）联系起来，绘制出不同基因组序列间构象恢复的精确机械加速图谱。

本文其余部分从理论构建到数值验证，全面分析DNA结构韧性。第2节建立研究的数学基础，提出受扰正弦-Gordon模型并定义支配双螺旋扭转动力学的物理参数。第3节详述计算方法，概述二阶St?rmer-Verlet积分格式、局域化力协议以及用于监测构象恢复的定量度量指标（如全局??2偏差范数）。第4节呈现详细的结果与讨论。本节利用时空热图演示辐射模态的动态解耦，并采用参数扫描建立DNA骨架刚度（??）与恢复时间尺度之间的逆幂律标度关系。最后，第5节提供结论，将结果综合到更广泛的生物学背景中，并讨论DNA作为遗传信息保真度拓扑机器的含义。

为了在非线连续力学与相应生物大分子结构性质之间建立清晰的数学桥梁，表1列出了本研究中主要使用的场变量、参数和谱度量。

模型构建

我们利用受扰sG方程对DNA双螺旋的非线性扭转动力学进行建模，该方程描述了核苷酸扭转角??沿连续生物大分子轴的时空演化。该数学表述扩展了配对DNA碱基的角振动与离散耦合摆链之间的经典类比，最初由Scott提出[39]。在未受扰的保守形式中，控制双曲偏微分方程的…

方法

为模拟拓扑DNA扭结的非平衡时空演化，方程(2)使用二阶St?rmer-Verlet（速度-Verlet）有限差分格式进行积分。选择这种几何辛积分器的原因在于其卓越的长期结构保持特性和在演化非线性Hamilton偏微分方程时极小的数值能量漂移。

连续空间域??∈[???/2,??/2]固定为渐近长度??=100，以…

讨论

通过将扭转形变建模为受扰sG方程中的拓扑孤子来研究DNA的机械韧性。我们的发现表明，DNA并非仅仅是静态的信息储存库，而是一个动态的非线性介质，能够进行精密的能量管理。通过对系统施加局域外部时空力??(??,??)，我们模拟了DNA结合蛋白（如RNA聚合酶或…）诱导的机械应力。

结论

在本工作中，我们利用受扰正弦-Gordon方程证明，DNA扭转动力学受复杂的非线性能量恢复机制支配。我们的数值模拟为双螺旋的结构韧性提供了关键见解，证实了即使在强烈的局域机械应力下，扭结孤子也能作为构象信息的稳健拓扑载体。

我们分析的核心——多谐波阻尼…

CRediT作者贡献声明

Mordecai Opoku Ohemeng：写作–审阅与编辑、写作–初稿、可视化、验证、监督、软件、资源、项目管理、方法论、调查研究、形式化分析、概念化。Joseph Ackora-Prah：写作–审阅与编辑、写作–初稿、可视化、验证、监督、项目管理、方法论、调查研究、形式化分析、概念化。Benedict Barnes：写作–审阅与编辑、写作–初稿、

利益冲突声明

作者声明他们没有已知的竞争性经济利益或个人关系，这些利益或关系可能会影响本文所述的工作。

Mordecai Opoku Ohemeng | Joseph Ackora-Prah | Benedict Barnes | Ishmael Takyi
爱达荷大学，工程学院，计算机系，莫斯科，爱达荷州，美国