[論文レビュー] Efficient Simulation of Tunable Lipid Assemblies Across Scales and Resolutions

John M. A. Grime, Jesper J. Madsen|arXiv (Cornell University)|Oct 11, 2019

Lipid Membrane Structure and BehaviorBiochemistry, Genetics and Molecular Biology参考文献 49被引用数 18

ひとこと要約

本稿では、複数スケールおよび解像度で調整可能な脂質アセンブリをシミュレートする計算的に効率的で粒子ベースの粗粒度モデルを提示する。ソフトで短距離作用の二体ポテンシャルとモジュラーなビード表現（1分子あたり2〜5ビード、および準モノレイヤーの1.5ビードモデル）を用いることで、モデルは調整可能な相挙動を示しながらミセル、バイレイヤー、バイスイクルを正確に形成する。既存のモデルと比較して数値的安定性と動的特性が向上し、曲げや相転移を含む膜物理学の主要な側面を捉えている。

ABSTRACT

We present a minimal model for simulating dynamics of assorted lipid assemblies in a computationally efficient manner. Our model is particle-based and consists of coarse-grained beads put together on a modular platform to give generic molecular lipids with tunable properties. The interaction between coarse-grained beads is governed by soft, short-ranged potentials that account for the renormalized hydrophobic effect in implicit solvent. The model faithfully forms micelles, tubular micelles, or bilayers (periodically infinite, bicelles, or vesicles) depending on the packing ratio of the lipid molecules. Importantly, for the self-assembled bilayer membranes it is straightforward to realize gel and fluid phases, the latter of which is most often of primary interest. We show that the emerging physics over a wide range of scales and resolutions is (with some restrictions) universal. The model, when compared to other popular lipid models, demonstrates improvements in the form of increased numerical stability and boosted dynamics. Possible strategies for customizing the models (e.g., adding chemical specificity) are briefly discussed. An implementation is available for the LAMMPS molecular dynamics simulator [Plimpton. J. Comp. Phys. 117, 1-19. (1995)] including illustrative input examples from the simulations we present.

研究の動機と目的

• 複数スケールおよび解像度で多様な脂質アセンブリをシミュレートする計算的に効率的で最小限のモデルの開発。
• ハードコアまたはレナードジョーンズポテンシャルを用いる既存の粗粒度脂質モデルと比較して、数値的安定性と動的特性の向上。
• 調整可能な相互作用パラメータを用いて、液体相とガラス相、パッケージング比、曲げ挙動などの調整可能な膜特性を実現すること。
• 化学的特異性をカスタマイズ可能な、モジュラーで拡張可能なフレームワークを提供し、膜関連の生物学的プロセスをシミュレートすること。
• 本モデルの普遍性を、大規模な形態的転移および膜物理学を再現する能力によって示すこと。

提案手法

• 1分子あたり2〜5ビードのモジュラー粒子ベースのフレームワークを用い、ヘッド、インターフェース、テール領域をそれぞれ異なる相互作用パラメータで表現する。
• 非共有相互作用は、ソフトで短距離作用の二体ポテンシャルで記述される：r ≤ r₀ の場合、U_pair(r) = - (A/a) sin(r·a)、r₀ < r < r_c の場合、- (B/b) sin(π/2 + (r - r_c)·c) で、r_c で滑らかにゼロに減少する。
• 分子内力は、調和的結合（K_b(r - r₀)²）および角度（K_θ(θ - θ₀)²）ポテンシャルによりモデル化され、脂質鎖の剛性を維持する。
• 疎水性効果は、インターフェースビードとテールビード間の凝集的相互作用に暗黙的に再正規化されるが、ヘッドビードは純粋に反発的相互作用を有する。
• モデルはLAMMPSに実装され、ミセル、バイレイヤー、バイスイクル、およびインシリコでの引き抜き実験のシミュレーションを通じて検証された。
• バイスイクルのパッケージングには、不安定性を回避しフリップフロップ緩和の必要を減らすために、脂質1分子あたりの体積（VPL）を一定に保つ戦略が用いられた。

実験結果

リサーチクエスチョン

• RQ1 ソフトで短距離作用のポテンシャルを有する最小限の粗粒度脂質モデルは、ミセル、バイレイヤー、バイスイクルといった多様な形態を自己組織化して形成できるか？
• RQ2 本モデルの動的特性および数値的安定性は、SDK やレナードジョーンズに基づくアプローチといった既存のモデルと比較してどのように異なるか？
• RQ3 本モデルは、相挙動（液体相対ガラス相）、曲げ、形態的転移といった膜物理学の主要な側面をどの程度正確に再現できるか？
• RQ4 本モデルの調整可能性を活用して、計算コストを低減しつつ、複雑な生物学的プロセスをシミュレートできるか？
• RQ5 脂質パッケージング戦略の選択（一定APL対VPL）が、バイスイクルの安定性および平衡化にどのように影響するか？

主な発見

• 本モデルは、脂質パッケージング比に応じてミセル、チューブ状ミセル、バイレイヤー（周期的、ビセルラ、バイスイクル形態を含む）を外部制約なしに自己組織化して形成する。
• 液体バイレイヤー相は安定しており、実際の動的挙動を示し、脂質の柔軟性と相互作用パラメータを調整することで相挙動を調整可能である。
• 準モノレイヤー・モデル（1分子あたり1.5ビード）は高い計算効率を示し、バイスイクルにおけるインシリコでの引き抜き実験を可能にし、膜の曲げおよび機械的性質を明らかにした。
• 特に極めて粗粒度化された解像度において、$1/r^{12}$ やレナードジョーンズポテンシャルを用いるモデルと比較して、数値的安定性が向上し、動的特性が向上している。
• 体積/分子あたり一定（VPL）のパッケージング戦略は、一定面積/分子あたり（APL）の戦略よりもより安定したバイスイクルを生成し、フリップフロップ緩和の必要を最小限に抑え、構造的アーチファクトを低減する。
• ソフトで短距離作用のポテンシャルにより、切断やシフトなしに滑らかに力が減少し、安定性が向上し、長時間シミュレーションが可能になる。