KIMMDY：一种生物分子反应仿真器

观看分子如何形成与断裂键

生命依赖于细胞内无数微小的化学反应，但要详尽地观察这些反应的全过程几乎不可能。实验只能看到模糊的平均结果，而计算模型在同时追踪大型生物分子的运动以及化学键在长时间尺度上的形成与断裂方面也面临困难。本文介绍了 KIMMDY，这是一种新型的模拟器，它不试图追踪反应的每一个量子细节，而是模拟反应在复杂生物环境中会如何进行。这为我们观察组织在应力下如何老化、DNA 如何被光损伤以及相互竞争的反应通路如何塑造生物学现象打开了一扇窗。

模拟化学反应的新途径

传统分子模拟擅长展示蛋白质和 DNA 等大分子如何运动和折叠，但通常将化学键视为不可断裂。真正的反应性模拟确实存在，但它们的计算成本极高，难以覆盖生物学上相关的长时间尺度，尤其当涉及多种反应类型和位点时。KIMMDY（Kinetic Monte Carlo Molecular Dynamics 的缩写）采取了不同的路径。它让常规模拟处理分子的运动，然后使用一个附加模块搜索可能发生反应的部位，估算每种反应的速率，并随机选择下一步发生的反应。通过反复在运动、反应搜索和反应选择之间循环，KIMMDY 能在无需计算完整量子级反应路径的情况下，追踪从亚秒到更长时间尺度的复杂反应级联。

反应仿真器如何工作

KIMMDY 的核心是一个“反应仿真器”，它根据分子的瞬时构型预测反应速率。对于分子动力学模拟中的每一个快照，KIMMDY 识别可能的反应伙伴，并使用不同模型估算每个反应的势垒大小。在许多情况下，它采用图神经网络——这类机器学习模型将原子视为节点、化学键视为连边——根据先前的高水平量子计算来猜测能垒高度。在其他情况下，它则依赖更简单的基于物理的公式或来源于实验的经验规则。由这些能垒计算出速率，并将其输入动力学蒙特卡罗步骤，后者随机选择下一次反应、推进模拟时间，再在下一轮分子运动前更新分子结构和力场。

追踪蛋白质支架中的自由基损伤

为了验证 KIMMDY 的有效性，作者首先将其应用于实验数据可得的小有机自由基的一类简单反应。仿真器重现了哪些氢迁移反应最可能发生，尽管它低估了绝对速率。随后他们研究了更复杂的案例：胶原蛋白——赋予结缔组织强度的纤维状蛋白。在机械应力下，胶原中的键可能发生均裂，产生高度反应性的自由基，这些自由基通过沿链传递氢原子而从一处跳到另一处。在一个由 260 万个原子组成的胶原纤维中，KIMMDY 追踪了数百次这样的跳跃，显示自由基可以从起源处迁移几纳米远。模拟揭示了一种不寻常的交联（称为 PYD）和一种修饰氨基酸（DOPA）都充当了异常有效的自由基捕获点，这一结论得到了热力学计算的支持，并可通过重新解读早期的电子自旋共振谱得到印证。

当不同反应发生竞争时

很多生物体系可以以不止一种方式断裂同一键，导致要么形成自由基，要么产生更常见的“闭壳层”产物。例如在胶原中，肽键既可以对称断裂形成两个自由基（均裂），也可以与水反应并不对称断裂（水解）。在现实环境中对两种选项进行直接量子模拟几乎不可行。借助 KIMMDY，作者将一个用于力促均裂的物理模型与基于实验的启发式水解模型结合，让这两种反应在单个肽段和致密的胶原纤维内竞争。他们发现，在被拉伸的孤立链中，水解明显占优；但在拥挤、交联的纤维中，机械应力在某些区域高度集中，以致均裂变得具有竞争力甚至更快，这解释了为何在真实组织中会观测到自由基形成。

光驱动的 DNA 损伤与交联

KIMMDY 还阐明了紫外光如何改变 DNA。当相邻的两个胸腺嘧啶吸收紫外光时，它们可能形成环丁烷二聚体——一种与皮肤癌密切相关的损伤，同时也被用来在 DNA 纳米技术中交联链条。作者使用一个简单的几何规则，将反应键之间的距离和角度与反应发生的可能性关联起来，并根据小分子实验进行调优，借助 KIMMDY 估算了不同 DNA 构象的“量子产率”——每吸收一个光子形成二聚体的频率。仿真器预测，一些常用于 DNA 摆构装（如悬挂端和交叉点）的排列，与标准双链或有缺口的 DNA 相比，其产率出人意料地低。它还表明形成一个二聚体并不会显著改变附近再形成另一个二聚体的概率，暗示先前的损伤并不一定会使 DNA 更易形成更多二聚体。

这对生物学与设计为何重要

简而言之，KIMMDY 提供了一种快速且灵活的方式来演算活体物质内部化学的“如果……会怎样”场景。通过仿真而非完全从头模拟反应，它能够处理庞大的系统和长时间尺度，同时仍然考虑到分子不断变化的构型如何影响反应在何处以及何时发生。这一方法帮助发现了胶原中先前被忽视的一个自由基稳定位点，阐明了机械力如何改变断键通路之间的平衡，并提出了关于光照下 DNA 交联形成的新问题。随着该仿真器扩展到更多反应类型和更准确的速率模型，它有望成为理解生物分子如何老化、失效或为医学与纳米技术工程化的强大工具。

引用: Hartmann, E., Buhr, J., Riedmiller, K. et al. KIMMDY: a biomolecular reaction emulator. Nat Commun 17, 3500 (2026). https://doi.org/10.1038/s41467-026-71955-2

关键词: 生物分子模拟, 反应动力学, 分子建模, DNA 损伤, 机械化学