[論文レビュー] 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におけるNa+およびCl−イオンに「小さな」(TQ5)ビーズが使用されることで、イオン-水およびイオン-イオンの衝突不安定性により、最大安定時間ステップが25 fsに制限されることを特定した。より長い時間ステップを可能にするために、著者らはイオン質量を増加するか、ビーズサイズを「小さな」(SQ5)または「通常の」(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−イオンを含むシミュレーションにおける時間ステップ不安定化の原因を特定すること。
• 25 fsの時間ステップ制限が、長時間スケールのシミュレーションにおける効率的サンプリングを妨げることを是正すること。
• 時間ステップ、イオン濃度、イオン質量を関数とするMDクラッシュレートの定量的運動論モデルを構築すること。
• より長い安定時間ステップを実現する2つの解決策—イオン質量の増加とビーズタイプの変更—の有効性を評価し、システムの性質に影響を与えないこと。
• ハイブリッドMDシミュレーションにおける最適時間ステップを決定する一般的なフレームワークを提供すること。

提案手法

• MDクラッシュレート kcrash(cion, ∆t) の運動論モデルを導出し、シミュレーションクラッシュ統計にグローバルにフィットさせ、不安定要因を定量化した。
• モデルを用いて、イオン-水およびイオン-イオン相互作用が時間積分不安定性に与える寄与を分析した。
• 古典的分配関数を保持しつつ、均一スケーリングによりイオン質量を増加（m′ = 2m または 4m）し、動的性質への摂動を最小限に抑えた。
• ビーズタイプを「小さな」(TQ5) から「小さな」(SQ5) もしくは「通常の」(RQ5) に変更し、Martini 3のマルチスケール哲学に従った水和状態の再定義を実施した。
• 時間ステップおよびイオンパラメータを変化させた際の、拡散係数、オーダーパラメータ、イオン分布、膜の流動性といったシステムの性質をモニタリングした。
• モデルは、NaCl溶液、DPPCバイレイヤー、PI-(4,5)P2含有膜を含む複数のシステムで検証され、一貫した結果が得られた。

実験結果

リサーチクエスチョン

• RQ1 Na+およびCl−イオンに「小さな」(TQ5) ビーズが使用されるMartini 3シミュレーションにおいて、時間ステップが何によって制限されるか。
• RQ2 イオン-水およびイオン-イオン相互作用は、MD時間積分における数値的不安定性にどのように寄与するか。
• RQ3 イオン質量を増加させることで、25 fsを超える時間ステップでも安定したシミュレーションが可能になり、平衡および動的性質への影響を最小限に抑えられるか。
• RQ4 イオンのビーズサイズをTQ5からSQ5またはRQ5に変更することで、システムの挙動を保持したままより長い時間ステップが可能になるか。
• RQ5 導出されたクラッシュレートの運動論モデルは、ハイブリッドMDキャンペーンにおける最適時間ステップを決定するために使用可能か。

主な発見

• TQ5イオンを用いたMartini 3シミュレーションにおける最大安定時間ステップは、イオン-水およびイオン-イオンの衝突に起因する不安定性により25 fsに制限される。
• イオン質量を2倍にすると30 fsでの安定シミュレーションが可能となり、4倍にすると40 fsでも安定するが、拡散係数や構造的性質への影響はわずかである。
• ビーズタイプをTQ5からSQ5またはRQ5に変更することで、30–40 fsでの安定シミュレーションが可能となり、Martini 3のマルチスケール設計と整合的である。
• 運動論的クラッシュモデルは不安定傾向を正確に予測でき、ハイブリッドシミュレーションにおける最適時間ステップの決定に利用可能である。
• 脂質の拡散係数およびオーダーパラメータは、時間ステップやイオン質量、ビーズタイプの変更に関わらずほとんど変化せず、摂動が小さいことが示された。
• 40 fsの時間ステップは、RQ5イオンを有する系でのみ安定であり、PI-(4,5)P2含有膜のような複雑な系には推奨されない。