分子动力学粗粒化模拟

一、项目背景与科学问题

生物膜是细胞与细胞器功能的基础，其结构性质（如单脂面积、膜厚度、脂质有序度）与膜蛋白功能密切相关。全原子分子动力学（AA-MD）模拟虽然精度高，但计算量巨大，难以研究微秒-毫秒尺度的膜自组装、相分离与融合等集体行为。粗粒化（Coarse-Grained, CG）方法通过将若干原子合并为一个相互作用位点，将计算效率提升 2–3 个数量级，同时保留膜的主要物理化学特征。MARTINI 力场是生物膜粗粒化模拟的黄金标准，其将 DPPC 磷脂映射为约 12 个 CG 珠子，采用 Lennard-Jones 势与库仑势描述非键相互作用。

本项目采用 MARTINI 粗粒化模型构建二棕榈酰磷脂酰胆碱（DPPC）双层膜，开展多微秒级分子动力学模拟，分析单脂面积（APL）、膜厚度、密度剖面与序参数等结构性质。

二、计算方法与理论背景

2.1 MARTINI 力场

MARTINI 力场采用约 4:1 的映射比，将重原子划分为四类相互作用类型（极性 P、非极性 N、弱极性 C、类碳氢 Q），每类又细分为多个子类型，以不同 ε 与 σ 参数刻画相互作用强度。水由一个 CG 珠子代表 4 个水分子。非键相互作用采用 LJ 势 V(r) = 4ε[(σ/r)¹² − (σ/r)⁶]，键合相互作用采用谐波势。该参数化策略使 MARTINI 能很好地再现膜的实验性质（面积、厚度、有序度）。

2.2 模拟细节

• 系综：等温等压（NPT），温度 323 K（高于 DPPC 主相变温度 314 K，保证液晶相），压力 1 bar 采用 Berendsen/半各向异性恒压器控制。
• 时间步长：20 fs，总模拟时长 200 ns。
• 周期性边界条件，静电采用反应场（reaction-field）方法。

三、计算设置与模型构建

• 初始构型：预平衡的 DPPC 双层（每层 64 脂质分子，共 128 个），周围填充 MARTINI 水。
• 力场：MARTINI 2.2；软件：GROMACS。
• 能量最小化后分别进行 100 ps 位置限制模拟与 200 ns 自由模拟。
• 分析：每 1 ps 输出原子坐标，用于计算 APL、密度剖面、厚度等。

四、结果分析与案例图解读

图 7 的左面板给出单脂面积 APL 随时间演化曲线。蓝色细线为瞬时值，红色粗线为滑动平均（窗口 20 ns）。体系从初始构型快速弛豫后，APL 围绕约 0.64 nm² 平衡振荡。平衡值 0.64 nm² 与实验 X 射线散射得到的 DPPC 液晶相 APL（0.62–0.66 nm²）高度一致，也符合 MARTINI 文献值（约 0.64 nm²），表明力场与模拟协议可靠。曲线在 200 ns 内保持稳定平台，说明双层已充分平衡。

图 7 的右面板给出沿膜法线方向的密度剖面：红色为磷脂头部基团（胆碱+磷酸酯）密度，蓝色为尾链（碳氢）密度。两个头部分布峰位于 z ≈ ±1.7 nm，对应双层两个界面；尾部密度在膜中心达到最大值约 740 kg/m³，形成连续的疏水核心。头-头距离约 3.4 nm，即双层疏水厚度；结合 APL，可计算脂质堆积参数并评估膜的曲率弹性。密度剖面表明 DPPC 双层结构完整、无脂质外泄，模拟体系健康。

两面板共同回答了”粗粒化模型能否再现真实膜结构”这一核心问题：APL 与密度剖面均与实验/文献定量吻合，证明 MARTINI 粗粒化模拟可在大幅降低计算成本的前提下，忠实再现生物膜的关键结构参数。

五、结论与展望

本项目通过 MARTINI 粗粒化 MD 模拟获得了 DPPC 双层膜的平衡结构，APL ≈ 0.64 nm²、头-头距离 ≈ 3.4 nm，与实验一致。案例图以时间序列与空间剖面两种视角完整呈现膜的结构特征。

后续可拓展：(1) 加入胆固醇或蛋白质，研究膜组分对结构性质的影响；(2) 分析序参数 S_CD 与脂质横向扩散系数，表征膜流动性；(3) 拓展至多组分不对称膜与膜融合等集体过程。