双链DNA模型中孤子结构与动力学行为的神经与传统方法比较

为何DNA中的波动重要

在每个细胞内，DNA不断弯曲、扭转与振动。这些运动并非纯粹的随机抖动；它们可以帮助打开遗传编码的片段，从而使基因能够被读取或复制。本文探讨了某些类波扰动——称为孤子——如何沿着一个简化的双链DNA模型传播。通过将传统数学工具与受神经网络启发的方法结合，作者绘制出可存在的波型，分析它们的稳定性，以及何时其运动可能变得高度敏感甚至出现混沌。

双链DNA的简化图景

真实的DNA是扭曲的双螺旋，但在计算上常被表示为由弹簧连接的两条平行链。在本研究中，作者使用一个双链模型来追踪两类运动：沿链的伸展和链间间隙方向的运动。两个变量描述这些运动，模型包含表示材料刚度、张力、质量和链间距的参数。为简化问题，作者假定在小扰动下两条链同步运动。这将模型化简为一个有效的单方程，同时保留描述紧致、保持形状波所需的关键非线性特性。

搜索保持形状的波

工作的核心是系统地寻找化简方程的孤子解。孤子是随着传播保持形状而不是扩散的特殊波。使用基于里卡提方程的技术——一种简单但有力的数学关系——作者将原本复杂的方程转化为可用双曲、三角函数以及有理表达式等熟悉构件求解的形式。从这些构件中他们构造出对应明亮脉冲、暗凹与连接两种不同背景水平的楔形台阶族解。这些解析表达式大大扩展了该类DNA模型已知的可能波形目录。

基于神经网络的精确解的变体

接着，作者引入一种修改的方法，借用神经网络的架构但以纯符号方式使用。与其在数据上训练网络，他们直接从先前的里卡提解中选择第一隐藏层的激活函数。网络输出被视为对DNA位移的巧妙猜测。当将这一猜测形式代回控制方程时，使方程精确成立的要求转化为关于网络权重与偏置的一组代数条件。求解这些条件产生新的封闭形式波解。通过比较三维曲面、等高线图和截面图，作者展示了基于神经网络的方法能够重现并扩展传统里卡提方法发现的形状，包括楔-反楔对、明亮与暗色脉冲，甚至奇异的尖峰结构。

从稳态波到敏感与混沌运动

除了列举波形，研究还考察了这些结构的动态行为及其对初始条件的敏感性。将化简方程重写为平面动力学系统后，波在每一点的状态用二维相平面上的一个点表示。在保持物理参数固定的情况下改变初始条件，作者观察到所得轨迹发生巨大变化，表明即使是极小的初始差异也能强烈影响波的演化。随后他们加入一个模拟外界影响的周期驱动项，并计算Lyapunov指数——量化相邻轨迹发散速度的数值。正的Lyapunov指数表明，在某些参数区间内，DNA波的动力学可以变得混沌，长期行为实际上不可预测。

这对DNA与未来模型意味着什么

通俗地说，文章表明一个简单的双链DNA图景能够承载丰富多样的传播扰动，从平滑的台阶与局部凸起到极端的尖峰事件。混合解析与神经启发的方法提供了一套工具，可用于精确写出这些波并探查它们何时稳定或容易出现剧烈、混沌的变化。尽管理论性强，这项工作支持了能量可以以有组织的包携形式在DNA中移动并可能促进基因激活等生物过程的观点。作者建议相同策略可进一步扩展以包含更现实的效应——例如变化的环境、外力或过去运动的记忆——为更忠实地数学建模DNA的动态生活提供途径。

引用: Majid, S.Z., Sağlam, F.N.K. & Ullah, M.S. Comparison of neural and traditional techniques for soliton structures and dynamical behavior in a double-chain DNA model. Sci Rep 16, 10390 (2026). https://doi.org/10.1038/s41598-026-41229-4

关键词: DNA孤子, 非线性波, 神经网络方法, 混沌动力学, 数学生物物理学