脂质囊泡在全原子、MARTINI和耗散粒子动力学中的结构和动力学一致性

《RSC Advances》：Structural and dynamical consistency in lipid vesicles across all-atom, MARTINI, and dissipative particle dynamicsOpen Access

　　

研究人员开发并验证了一种介观耗散粒子动力学（dissipative particle dynamics, DPD）-4脂质模型，该模型基于已建立的MARTINI力场进行参数化。研究人员将该DPD模型和用于粗粒化（coarse-grained, CG）脂质囊泡的MARTINI模型与全原子（All-Atom）模拟进行基准比较，以评估它们在各自尺度上捕捉结构性质的能力。动力学性质则通过CG模型之间的相互比较以及与实验趋势的对比进行评估。研究发现，CG模型能够捕捉直径分别小于和大于50 nm的超小单层囊泡（ultrasmall unilamellar vesicles, USUVs）和小单层囊泡（small unilamellar vesicles, SUVs）的多尺度时空性质。具体而言，MARTINI和介观DPD脂质模型在捕捉USUVs的局部结构变化（如曲率依赖的脂质堆积）方面具有一致性。对于该特定系统，结构一致性已通过全原子模拟得到验证。在较短时间尺度上，动力学一致性已通过MARTINI与DPD模型之间的吻合得到验证。60 nm SUVs的CG动力学（其短时间尺度与全原子模拟匹配）表现出与相应实验中观察到的相同的拉伸指数弛豫行为。这种近多尺度验证凸显了在长度和时间上评估CG模型的更严格要求。通过这项工作，研究人员还展示了使用MARTINI和DPD模型来扩展囊泡模拟并在介观尺度区域提供动力学信息且具有显著计算节省的潜力。

脂质囊泡跨尺度模拟的结构与动力学一致性研究解读

在复杂化学系统中，非平衡行为往往需要在保留分子相互作用本质的同时进行表征，这促使多尺度——尤其是粗粒化（coarse-grained, CG）——表示方法的发展。CG方法通过将原子分组为珠子，扩展了构型空间的采样范围，但同时也丢失了被去除的精细相互作用效应，且粒子运动的加速会使动力学发生扭曲。因此，评估CG模型在不同尺度上是否保留结构一致性和动力学一致性至关重要。结构一致性指不同模型在局部和整体结构性质上的吻合，例如对关联函数g(r)的一致；动力学一致性指时间相关性质的一致，例如相关函数C(t)的一致。MARTINI力场由Marrink等人开发，基于四对一映射，广泛应用于脂质膜、蛋白质和烃类体系的模拟；耗散粒子动力学（dissipative particle dynamics, DPD）则采用软排斥势，可达到更长的时间和长度尺度。然而，这两种CG方法在脂质囊泡体系中的系统比较和验证仍不充分。为此，研究人员开展了本研究，以揭示MARTINI和DPD在复杂体系中的跨尺度一致性。

研究人员开发并验证了一种基于MARTINI力场参数化的介观DPD-4脂质模型，用于1,2-二肉豆蔻酰-sn-甘油-3-磷酸胆碱（DMPC）脂质囊泡，并与全原子（All-Atom）和MARTINI模型在12–60 nm尺寸范围内进行了系统比较。研究结果表明，CG模型能够捕捉超小单层囊泡（ultrasmall unilamellar vesicles, USUVs）和小单层囊泡（small unilamellar vesicles, SUVs）的多尺度时空性质；结构一致性通过全原子模拟验证，动力学一致性通过MARTINI与DPD模型间及实验趋势验证。60 nm SUVs的动力学遵循Zilman–Granek拉伸指数弛豫形式。该研究提供了跨尺度验证和参数化流程，展示了CG模型在介观尺度模拟中的潜力。论文发表在《RSC Advances》上。

关键技术方法方面，研究人员使用了三种模拟方法：全原子分子动力学（molecular dynamics, MD）模拟（CHARMM36力场，LAMMPS软件）采用预填充方法构建12 nm囊泡；MARTINI粗粒化模拟（MARTINI 2.2力场，P4水模型，GROMACS软件）通过CHARMM-GUI生成不同尺寸（12、17、25、60 nm）的DMPC囊泡；DPD介观模拟采用自行参数化的DPD-4模型（LAMMPS软件），该模型从MARTINI珠映射和键合项推导，并引入静电相互作用。样本为DMPC脂质囊泡，脂质数量从640到27,848不等。

研究结果包括以下几个方面：

3.1 Lipid counts in inner and outer layers（内外层脂质数量）：通过理论公式和模拟数据对比，研究人员发现囊泡外层与内层脂质数量之比随半径增大而非线性变化，模拟结果与理论预测高度一致，验证了模型对脂质分布的准确性。

3.2 Computational performance（计算性能）：对比了三种方法的计算成本。全原子模拟12 nm囊泡5 ns需127小时（144核），MARTINI模拟1000 ns需7.24小时（48核），DPD-4模拟137 ns需1.08小时（48核）。DPD-4模拟60 nm囊泡约需59 core-hours/ns，比全原子外推估计快约2700倍，凸显了CG方法的计算优势。

3.3 Vesicle structure（囊泡结构）：通过径向分布函数、密度分布和脂质尾部末端距（r_ee）等指标，发现三种模型在12 nm囊泡上的结构性质基本一致。Kolmogorov–Smirnov（KS）检验表明，DPD-4在头部基团（铵根和磷酸根）更接近全原子参考，MARTINI在尾部碳区更接近。随着囊泡尺寸增大，低r_ee的脂质群体减少，表明脂质堆积随曲率降低而更紧密。

3.4 Structural consistency（结构一致性）：利用面积/脂质值0.66 nm²构建60 nm囊泡，确定内外层脂质数量为12,673和15,175，并通过Packmol构建了27,848个脂质的囊泡，验证了理论表达式在更大尺寸下的适用性。

3.5 Dynamical consistency（动力学一致性）：计算了范霍夫函数G_s(r,t)和中间散射函数F_s(q,t)。DPD-4与全原子参考的吻合优于MARTINI（后者需时间重标定因子4）。两种模型均符合Zilman–Granek拉伸指数弛豫形式，但拟合前因子Γ与实验值（20–26 ?³ ns^?1）相差约4倍。将分析窗口延长至150 ns后，低q区域的拟合显著改善（q = 0.18 nm^?1处R²从?0.16升至0.64），Γ值更接近实验。

讨论部分指出，结构一致性在粗粒化映射中必然在长尺度上精确，而短时间尺度被平均化。动力学一致性因随机项的存在难以严格证明，但空间粗粒化等效于时间粗粒化，因此CG模型在长时间尺度上可以保持动力学一致性。研究结论可翻译为：本工作的关键发现是，从MARTINI珠映射和键合项参数化的DPD-4模型，与MARTINI具有结构和动力学一致性，并与全原子模型在脂质囊泡系统中具有结构一致性。因此，全原子和MARTINI基准模型可用于结构参数化DPD囊泡系统，以访问类似的长度和时间尺度。DPD模型还可应用于足够大、可被实验探测但其他模型难以模拟的USUVs。60 nm SUVs的膜动力学遵循实验观察到的Zilman–Granek拉伸指数弛豫，但定量前因子与实验差约4倍，需要进一步参数化。该方法可模拟12–60 nm的SUVs，未来扩展至巨型单层囊泡（GUVs）和最小全细胞模型需要验证多组分、含蛋白质的膜，并涉及更长的尺度和更粗的动力学，超出了本文验证的12–60 nm微秒量级。