[论文解读] Small ionic radii limit time step in Martini 3 molecular dynamics simulations

Balázs Fábián, Sebastian Thallmair+1 more|arXiv (Cornell University)|Apr 11, 2022

Protein Structure and Dynamics参考文献 46被引用 16

一句话总结

本文指出，Martini 3 中使用 'tiny'（TQ5）粒子表示 Na+ 和 Cl− 离子，导致离子-水及离子-离子碰撞不稳定性，使最大稳定时间步长被限制在 25 fs。为实现更长的时间步长，作者提出增加离子质量或升级粒子尺寸至 'small'（SQ5）或 'regular'（RQ5），这两种方法均可在对结构和动力学性质影响极小的情况下实现 30–40 fs 的稳定模拟，显著提升采样效率。

ABSTRACT

Among other improvements, the Martini 3 coarse-grained force field provides a more accurate description of the solvation of protein pockets and channels through the consistent use of various bead types and sizes. Here, we show that the representation of Na$^+$ and Cl$^-$ ions as "tiny" (TQ5) beads limits the accessible time step to 25 fs. By contrast, with Martini 2, time steps of 30-40 fs were possible for lipid bilayer systems without proteins. This limitation is relevant for, e.g., phase separating lipid mixtures that require long equilibration times. We derive a quantitative kinetic model of time-integration instabilities in molecular dynamics (MD) as a function of time step, ion concentration and mass, system size, and simulation time. With this model, we demonstrate that ion-water interactions are the main source of instability at physiological conditions, followed closely by ion-ion interactions. We show that increasing the ionic masses makes it possible to use time steps up to 40 fs with minimal impact on static equilibrium properties and on dynamical quantities such as lipid and ion diffusion coefficients. Increasing the size of the bead representing the ions (and thus changing their hydration) also permits longer time steps. The use of larger time steps in Martini 3 simulations results in a more efficient exploration of configuration space. The kinetic model of MD simulation crashes can be used to determine the maximum allowed time step whenever sampling efficiency is critical.

研究动机与目标

• 识别在 Martini 3 模拟中使用 Na+ 和 Cl− 离子时时间步长不稳定的原因。
• 解决 Martini 3 中 25 fs 时间步长的限制，该限制阻碍了长时序模拟中的高效采样。
• 建立时间步长、离子浓度和离子质量作为函数的 MD 模拟崩溃率的定量动力学模型。
• 评估两种解决方案——增加离子质量与改变粒子类型——在不损害系统性质的前提下实现更长且稳定的模拟时间步长的可行性。
• 为高通量 MD 模拟中确定最优时间步长提供通用框架。

提出的方法

• 推导并全局拟合了 MD 模拟崩溃率 kcrash(cion, ∆t) 的动力学模型，以量化不稳定的来源。
• 将该模型应用于分析离子-水和离子-离子相互作用对时间积分不稳定性的影响。
• 通过均匀缩放（m′ = 2m 或 4m）增加离子质量，同时保持经典配分函数不变，并将动力学扰动降至最低。
• 将粒子类型从 'tiny'（TQ5）更改为 'small'（SQ5）或 'regular'（RQ5），以符合 Martini 3 的多尺度设计哲学，改变水合行为。
• 在不同时间步长和离子参数下，监测了系统性质——包括扩散系数、序参数、离子分布和膜流动性。
• 在多个体系中验证了该模型：NaCl 溶液、DPPC 膜以及含 PI-(4,5)P2 的膜，结果一致。

实验结果

研究问题

• RQ1 当在 Martini 3 模拟中使用 'tiny'（TQ5）粒子表示 Na+ 和 Cl− 离子时，是什么限制了时间步长？
• RQ2 离子-水和离子-离子相互作用如何导致 MD 时间积分中的数值不稳定性？
• RQ3 增加离子质量是否可实现超过 25 fs 的稳定模拟，且对平衡态和动力学性质影响极小？
• RQ4 将离子粒子尺寸从 TQ5 更改为 SQ5 或 RQ5 是否可实现更长的时间步长，同时保持系统行为不变？
• RQ5 所推导的崩溃率动力学模型是否可用于确定高通量 MD 模拟中的最优时间步长？

主要发现

• 使用 TQ5 粒子的 Martini 3 模拟中，最大稳定时间步长受限于离子-水和离子-离子碰撞引起的不稳定性，上限为 25 fs。
• 离子质量加倍可实现 30 fs 的稳定模拟，质量增至四倍则可实现 40 fs，对扩散系数和结构性质仅造成轻微影响。
• 将粒子类型从 TQ5 更改为 SQ5 或 RQ5 可实现 30–40 fs 的稳定模拟，与 Martini 3 的多尺度设计原则一致。
• 动力学崩溃模型能准确预测不稳定性趋势，可用于确定高通量模拟的最优时间步长。
• 脂质的扩散系数和序参数在不同时间步长及离子质量或粒子类型修改下基本保持不变，表明扰动极小。
• 40 fs 时间步长仅在使用 RQ5 粒子的体系中稳定，不推荐用于含 PI-(4,5)P2 的复杂体系。