自组装分子动力学模拟：两亲性嵌段共聚物的相行为预测

嵌段共聚物的自组装比小分子表面活性剂复杂得多——它的每一个”构件”本身就是一个聚合物链段，链段的长度、比例、柔韧性都会影响最终的自组装形态。全原子模型能给出最精确的描述，但计算量限制了模拟的体系规模和时间尺度。粗粒化模型是解决这个矛盾的有效手段。

项目背景：Pluronic的温度响应行为

团队接手过一个PEO-PPO-PEO三嵌段共聚物（Pluronic F127）的温度响应性研究项目。F127是生物医药领域广泛使用的温敏性嵌段共聚物，在低温下以单分子链形式溶解于水，温度升高到临界胶束温度（CMT）以上时自发形成胶束。客户需要知道不同温度下的胶束尺寸和聚集数，用来优化药物递送体系的配方。

实验上可以用动态光散射（DLS）测量胶束的流体力学半径，但DLS只能给出平均值，无法区分球形胶束和棒状胶束的共存状态。全原子MD模拟在这种体系规模下（几百个嵌段共聚物分子，每个分子包含上百个原子）是不可行的——即使是粗粒化模型，也需要合理的参数映射策略。

Martini粗粒化建模

团队选择了Martini 3.0力场——这是目前应用最广泛的粗粒化力场，将4个重原子映射为一个粗粒化珠子。在Martini映射规则下，PEO链段用P5（极性）和P4（中等极性）珠子表示，PPO链段用SP1（非极性）珠子表示。

映射的具体方案是：每个EO单元映射为一个P5珠子，每个PO单元映射为一个SP1珠子。F127的分子式为PEO₁₀₀-PPO₆₅-PEO₁₀₀，对应粗粒化模型中约有230个珠子，全原子模型中约有1350个原子。粗粒化后体系规模缩小了约6倍，模拟速度提升了约100倍。

将200个F127分子随机分散在22×22×22 nm³的水盒子中，加水珠子和Na⁺离子（中和PO段的弱电离），在298 K和310 K两个温度下各运行1 μs的NPT模拟。

温度驱动的胶束化转变

298 K（低于CMT）的模拟结果：F127分子主要以单分子链的形式分布在水中，偶有小聚集体（3-5个分子）短暂形成又迅速解散。Rg（回旋半径）分布集中在3.5-4.2 nm范围，与文献中DLS测量的流体力学半径5.8 nm（无水合层时约为3.8 nm）基本一致。

310 K（高于CMT）的模拟结果发生了显著变化——200个分子在约200 ns内自组装形成了8-10个球形胶束，平均聚集数22±6。大胶束的Rg约6.5-7.5 nm，对应流体力学半径约9-10 nm。这个值与实验DLS在37°C下测得的25-28 nm有约3倍的差距，但这在Martini力场的预期范围内——粗粒化模型对链段柔性的描述偏高，导致胶束内部链段分布更松散。

更有价值的一个发现是：在310 K模拟的后期（500-1000 ns），出现了2个胶束融合成1个较大胶束的事件，以及1个胶束分裂成2个小胶束的事件。这种动态的融合-分裂过程解释了DLS数据中多分散指数（PDI）偏高（0.15-0.25）的原因——体系确实不是单分散的球形胶束，而是在动态平衡中存在大小分布。

粗粒化模型的适用边界

需要坦诚地说，Martini力场在这个项目中的表现是”定性正确、定量有偏差”。它能正确预测CMT附近胶束化的发生、胶束的大致尺寸和聚集数量级，但绝对值与实验的偏差在30-50%。

偏差的主要来源有两个：一是粗粒化珠子之间相互作用的平滑化处理，无法精确描述链段间的熵效应；二是4:1映射对PEO链段的水合行为描述不够精确。如果项目需要更高精度的尺寸预测，需要回退到全原子模型（但体系规模需要缩小到10-20个分子）或者使用迭代玻尔兹曼反演（IBI）方法重新拟合粗粒化参数。

客户对这个结果的评价是”定性趋势完全正确，定量精度够用于配方初筛”。他们用模拟筛选了3种不同分子量的Pluronic（F68、F127、F108），最终选定了F127用于后续实验——这个选择与模拟预测的胶束化温度和尺寸范围完全一致。粗粒化MD不是精度最高的工具，但在项目初筛阶段，它的效率优势是无可替代的。

如需了解更多的科研计算案例与项目经验，可以参考 https://www.keyanxueshu.com/ 上的技术文章。