[论文解读] 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

一句话总结

该论文提出了一种计算高效的基于粒子的粗粒化模型，用于在多尺度和多分辨率下模拟可调的脂质组装体。通过使用软短程成对势能和模块化粒子表示法（每脂质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.

研究动机与目标

• 开发一种计算高效、极简的模型，用于在多尺度和多分辨率下模拟多样化的脂质组装体。
• 与现有粗粒化脂质模型相比，提升数值稳定性和动力学性能，特别是针对使用硬核或Lennard-Jones势能的模型。
• 通过调节相互作用参数，实现膜性质的可调性，如液态与凝胶相、堆积比和弯曲行为。
• 提供一个模块化、可扩展的框架，用于模拟具有可定制化学特异性的膜相关生物过程。
• 展示该模型在重现大规模形态转变和膜物理特性方面的普适性。

提出的方法

• 该模型采用模块化粒子框架，每脂质包含2–5个粒子，分别表示头部、界面和尾部区域，并赋予不同的相互作用参数。
• 非键合相互作用由一种软短程成对势能控制：当 r ≤ r₀ 时，U_pair(r) = - (A/a) sin(r·a)；当 r₀ < r < r_c 时，U_pair(r) = - (B/b) sin(π/2 + (r - r_c)·c)，并在 r_c 处平滑衰减至零。
• 分子内作用力通过谐振键势（K_b(r - r₀)²）和键角势（K_θ(θ - θ₀)²）建模，以维持脂质链的刚性。
• 疏水效应通过界面与尾部粒子之间的内聚相互作用隐式归一化，而头部粒子仅具有排斥性相互作用。
• 该模型在LAMMPS中实现，并通过胶束、双分子层、囊泡以及体外拉伸实验的模拟进行了验证。
• 囊泡堆积采用恒定每脂质体积（VPL）策略，以避免不稳定性并减少翻转弛豫的需求。

实验结果

研究问题

• RQ1 一个采用软短程势能的极简粗粒化脂质模型能否自组装形成包括胶束、双分子层和囊泡在内的多样化形态？
• RQ2 与SDK或基于Lennard-Jones势能的成熟模型相比，该模型的动力学和数值稳定性如何？
• RQ3 该模型在多大程度上能再现关键膜物理特性，如相行为（液态与凝胶相）、弯曲性以及形态转变？
• RQ4 能否利用该模型的可调性，在计算成本降低的前提下模拟复杂的生物过程？
• RQ5 脂质堆积策略的选择（恒定每脂质面积APL vs. 恒定每脂质体积VPL）如何影响囊泡的稳定性和平衡化？

主要发现

• 该模型成功形成胶束、管状胶束和双分子层（包括周期性、双联和囊泡形式），其形态取决于脂质堆积比，且自组装过程无需外部约束。
• 液态双分子层相稳定，表现出真实的动力学行为，通过调节脂质柔性和相互作用参数可实现相行为的调控。
• 准单层模型（每脂质1.5个粒子）展现出极高的计算效率，支持对囊泡进行体外拉伸实验，揭示了膜的弯曲和机械特性。
• 与采用 $1/r^{12}$ 或Lennard-Jones势能的模型相比，该模型在高度粗粒化分辨率下展现出更优的数值稳定性和增强的动力学性能。
• 恒定每脂质体积（VPL）堆积策略相比恒定每脂质面积（APL）能产生更稳定的囊泡，显著减少翻转弛豫需求并降低结构伪影。
• 软短程势能确保力的平滑衰减，无需截断或位移处理，从而提升稳定性并支持长时间模拟。