甲基丙烯酸酯聚合物如何作用于细菌外膜：粗粒化分子动力学给出的四阶段机制

本文信息

• 标题：通过粗粒化分子动力学模拟探索甲基丙烯酸酯聚合物与细菌外膜的相互作用
• 作者：Eduardo R. Almeida、Vinicius Firmino dos Santos、Madeleine Ramstedt、Thereza A. Soares
• 发表时间：2026年4月7日
• 单位：圣保罗大学（巴西）、于默奥大学（瑞典）、奥斯陆大学与 Hylleraas Centre（挪威）
• 引用格式：Almeida, E. R., Firmino dos Santos, V., Ramstedt, M., & Soares, T. A. (2026). Exploring the Interactions Between Methyl Methacrylate Polymers and the Bacterial Outer Membrane via Coarse-Grained Molecular Dynamics Simulations. Journal of Chemical Information and Modeling. https://doi.org/10.1021/acs.jcim.6c00729
• 源代码：https://github.com/BioMat-USP-RP/Input-files-for-CG-simulations-of-polymers-and-bacterial-outer-membrane

摘要

聚合物刷涂层为对抗医疗器械中的细菌黏附与生物膜形成提供了一种有前景的策略。然而，不同刷层化学组成如何与细菌膜相互作用的详细分子层面理解仍不完整。在本研究中，我们使用粗粒化分子动力学模拟（steered molecular dynamics 与 umbrella sampling），研究了四种甲基丙烯酸甲酯衍生聚合物——pDMAEMA（弱阳离子）、pMETAC（强阳离子）、pMEDSAH（两性离子）和 pSPMA（阴离子）——与大肠杆菌细菌外膜（OM）模型的相互作用与转运过程。模拟揭示了一个四步转运过程：接近、黏附、渗入和内化，并由不同的热力学与动力学特征所表征。阳离子聚合物与外膜表面表现出明显有利的黏附，尤其是与 LPS 分子的糖类内核结构域，这主要归因于有利的静电相互作用。在这些带正电聚合物的转运过程中，还可观察到 LPS 单元被拖拽至细菌外膜内叶的现象。相比之下，两性离子与阴离子聚合物表现出较不有利的黏附，这与其抗污行为一致。该方法提供了一个计算框架，可在分子细节层面解析与聚合物-膜相互作用相关的自由能图景与结构扰动，包括对动力学上不利过程的预测，例如聚合物向膜细胞内区域的渗入与内化。这些结果为理解水合、电荷与聚合物结构如何影响细菌膜相互作用提供了机制见解，并推动了抗污与抗菌表面涂层的分子设计。

核心结论

• 阳离子体系界面吸附更强：pMETAC 和 pDMAEMA 在黏附力与黏附自由能上均更占优。
• 水合作用决定抗黏附特征：pMEDSAH 与 pSPMA 在后期保持更高水合，降低深层膜耦合。
• 跨膜过程动力学受限明显：所有体系在膜中心附近都面临不同程度的渗入势垒。
• 强吸附不等于顺利穿膜：approach 与 adhesion 可自发发生，但 permeation 与 internalization 都是动力学不利步骤。
• 阳离子链会牵引膜组分重排：明确观察到 LPS 被正电聚合物拖向内侧，并在后期形成瞬态 nanopores。
• 材料设计应分目标优化：表面捕获能力和跨膜推进能力需要拆开设计，单一电荷指标无法覆盖全流程行为。

背景

分子刷涂层被广泛用于医疗器械表面的抗黏附与抗菌改性，因为它们可以通过化学组成调控在生理介质中的稳定性、生物相容性和界面功能。以亲水聚电解质刷为例，体系通常分为强聚电解质和弱聚电解质两类，前者在宽 pH 范围维持带电，后者则随环境 pH 改变电离状态并形成可切换界面。对产业端而言，这决定了导管、植入物和传感界面在血液、血清、唾液等复杂体系中的失效模式；对学术端而言，这意味着刷层化学、水合结构与生物相互作用之间需要可量化的分子机制映射。

已有研究已经说明，抗污能力与界面水合层强度高度相关。前期 MD 与 Monte Carlo 工作表明，两性离子刷层通常比 PEG 或非两性离子亲水聚合物形成更稳定的界面水网络，从而更有效抑制蛋白吸附；同时，碳间隔长度、偶极取向与局部溶剂化排斥会进一步放大这种差异。问题在于，这些结论主要建立在蛋白-刷层模型上，而细菌外膜（尤其是革兰阴性菌外膜）在成分、拓扑、电荷分布和疏水性上都远比单蛋白目标复杂，外层 LPS、离子桥联和膜不对称结构共同抬高了建模与解释难度，也抬高了机制外推门槛。

因此，这个方向的核心 gap 是缺少能同时解释黏附、渗入与内化全过程的分子级统一框架。尤其在带电单元比例变化时，实验已经观察到抗污与抗菌行为可被显著调制，但机制上仍不清楚：是静电吸附主导，还是水合屏障主导，或者两者在不同阶段交替主导。本文的意义就在于把问题从终点表征推进到过程分解，用粗粒化 SMD 与 US 自由能图景把“接近—黏附—渗入—内化”串成可比较路径，为后续刷层配方设计提供可执行判据。

关键科学问题

• 在 LPS 主导的外膜界面，聚合物最先被什么物理作用抓住，静电吸附和脱水代价谁先主导。
• 强黏附是否能自然转化成强渗入，还是会出现界面停留很强但向内推进困难的状态。
• 四类聚合物的差异能否被统一机制解释，并转换为可执行的设计参数。

创新点

• 同平台并行比较：四类聚合物在同一 OM 与同一模拟流程下比较，减少跨体系偏差。
• 路径与自由能联动：把 SMD 的时间分段与 US 的 PMF 结果联动解释，不只看单一指标。
• 指标体系更完整：力学、自由能、水合、离子分布和接触统计共同构成解释框架。

研究内容

方法详述：模型、参数与流程

本文采用 MARTINI 3 粗粒化框架，核心对象是四条长度一致的聚合物链（每条 96 个单体）与不对称 E. coli 外膜体系，形成同平台可比体系。

体系组成

SI 中还给出完整组分表（Table S2）：外膜外叶含 560 个 LPS 分子，内叶含 1260 个 DPPE 和 420 个 DPPG（比例 75:25），离子包括 3012 个 $\ce{Ca^{2+}}$ 和 4 个 $\ce{Cl^-}$，水珠 616312 个。$\ce{Ca^{2+}}$ 在这里主要起桥联作用，通过与带负电的 LPS 和 DPPG 相互作用稳定外膜结构——这在革兰阴性菌中是保守机制。

图1：四类甲基丙烯酸酯聚合物单体与粗粒化映射关系

• 图1A–图1D分别对应 pDMAEMA、pMETAC、pSPMA、pMEDSAH 的单体化学结构及其 CG bead 映射。
• 颜色说明：黑色线条为原子级化学结构，蓝色与粉色球表示映射后的粗粒化表示，灰色标记代表不同 bead 类型。

图1定义了后续相互作用分析的化学语义。后文看到的水合差异、黏附差异和离子相互作用，都由这些 bead 化学属性决定，属于参数驱动的结构结果。

四类聚合物参数如何构建（SI Table S1）

参数生成流程在 SI 中：四条链都设为 96-mer；总电荷取各单体电荷求和；体积由 GROMACS 2019.4 的 SASA 相关排除体积估算得到。这个参数化流程直接决定了后续黏附强度、去溶剂化代价和渗入势垒的排序，是全文的比较基线。

图2：细菌外膜模型与反应坐标定义

• 图2A–图2C给出 LPS、DPPE、DPPG 的结构与粗粒化表示。
• 图2D–图2E展示不对称外膜组装和沿膜法向推进的反应坐标。
• 颜色说明：lipid A 为粉色，LPS inner core 为紫色，outer core 为橙色，水珠为蓝色，黑色箭头表示反应坐标方向。

图2明确了”聚合物在什么环境中前进”这个前提。LPS 的分层结构决定了电性匹配的空间选择性：从外到内依次是 outer core（glucose/galactose 糖残基，弱负电）、inner core（KDO 和 HEP 糖残基，强负电）、lipid A（疏水尾链）。$\ce{Ca^{2+}}$ 主要聚集在 inner core 区（z ≈ 2.0–2.5 nm），与强负电糖残基形成桥联。阳离子聚合物的 PMF 极小值也在这个位置，说明它们会被 inner core 的强负电”抓住”；而两性离子或阴离子聚合物的极小值更靠外（z ≈ 2.7 nm），在 outer core 区就更早受到排斥，难以继续深入。

还需要强调一点：外膜不对称性本身就是机制的一部分。如果把体系简化成对称磷脂双层，很多“先被外层糖基区捕获，再向疏水核心推进”的路径特征会被弱化，最终导致对抗菌刷层设计的判断偏乐观。

关键模拟设置

• 非键相互作用 cutoff 统一为 1.2 nm，静电使用 reaction field。
• 先做能量最小化，再做平衡，再进入 SMD 与 US。
• SMD 设置为沿膜法向的质心拉动，平均路径约 18 nm。
• 拉速为 $0.0001~\mathrm{nm/ps}$，弹簧常数为 $1000~\mathrm{kJ\cdot mol^{-1}\cdot nm^{-2}}$。
• 单条 SMD 轨迹约 320 ns，并进行 3 次重复。
• US 约 36 个窗口，间距约 0.3 nm，每个窗口采样 125–300 ns（不同体系时间不同），WHAM 重建 PMF，并做 bootstrap 误差估计。

如果把这套流程说得更直白一点，本文其实做了两件互补的事。SMD 更像是在给出一条可比较的推进路径，US/PMF 则负责把这条路径转换成自由能图景。

• 第一步，用 SMD 把一条自由聚合物链从体相水中缓慢推向外膜中心，记录这一路上受力、接触、水合和膜重排怎么变化；这一步解决的是“过程长什么样”。
• 第二步，从这条路径上挑出一系列代表性构型做 umbrella sampling，再用 WHAM 重建 PMF；这一步解决的是“哪一段热力学有利，哪一段动力学更难”。

方法流程图

graph TB
    subgraph S1["1.体系构建"]
    	direction TB
        A1["聚合物：96-mer"]
        A2["外层LPS + 内层DPPE/DPPG"]
        A3["310 K，1 bar，离子与水环境"]
    end

    subgraph S2["2.路径采样"]
        B1["平衡模拟"]
        B2["SMD三重复"]
        B3["时间分段：接近/黏附/渗入/内化"]
        B1 --> B2 --> B3
    end

    subgraph S3["3.热力学重建"]
    	direction LR
        C1["US分窗口采样"]
        C2["WHAM重建PMF"]
        C3["窗口重叠 + 分块收敛检查"]
        C1 --> C2 --> C3
    end

    S1 --> S2 --> S3 --> D["联动分析：受力-接触-水合-离子-自由能"]

结果一：四阶段路径的构象证据

SMD 轨迹把全过程稳定分为四段：approach（0–18 ns）、adhesion（18–34 ns）、permeation（34–112 ns）、internalization（112–320 ns），给出可重复的阶段边界。但构象快照只是第一条线索——要确认这个划分是否真实反映物理过程，还需要从受力、自由能、水合和离子分布等多个角度交叉验证。

图9：聚合物跨膜转运四阶段的构象快照。从 SMD 模拟中提取的聚合物跨细菌外膜（OM）转运的四个阶段快照：接近（A）、黏附（B）、渗入（C）和内化（D）。每个子图展示该阶段聚合物、外膜组分的典型构象及膜响应。

• A为接近阶段（0–18 ns）：聚合物位于膜外侧约 6.0 nm 范围内，开始去溶剂化但尚未与膜接触。
• B为黏附阶段（18–34 ns）：聚合物贴近 LPS 外层并形成稳定界面吸附，阳离子聚合物与 LPS 糖类内核区域相互作用更强。
• C为渗入阶段（34–112 ns）：聚合物向膜内推进，受力抬升，伴随局部膜结构重排、去溶剂化和瞬态缺陷。
• D为内化阶段（112–320 ns）：聚合物到达膜内叶，部分 LPS 分子被从外叶拖向内叶，膜表面形成纳米孔，随后膜结构自发重建。

四阶段划分来自构象快照、受力曲线、PMF 以及水合和离子分布的交叉一致性。approach 和 adhesion 可以自发发生，permeation 和 internalization 则对应动力学不利步骤。

结果二：膜结构的稳定性验证

密度分布结果显示，外膜宏观层状结构总体稳定，没有出现持续性大破裂。与此同时，进入内化阶段后，膜中心出现局部含水增强，中心水密度约为 $100~\mathrm{kg\cdot m^{-3}}$。这一观察与构象快照中的“瞬态缺陷”相互印证：聚合物推进确实扰动了膜，但扰动是局部的、动态的，而非整体破裂。

图3：四阶段中外膜与水的质量密度沿 z 轴分布

• 图3A–图3D分别对应接近、黏附、渗入、内化阶段。
• 线型说明：实线表示膜组分密度，虚线表示水密度，银色参考线为无聚合物扰动时的膜分布。

图3最核心的信息是整体结构稳定，但局部会被拉出动态缺陷。这类信号对应局部扰动窗口，支撑了本文对穿膜机制的保守判断：聚合物推进依赖局部重排与短时缺陷，体系不具备低阻力自由穿透通道。

结果三：离子重排与接触分析揭示静电匹配的空间选择性

$\ce{Ca^{2+}}$ 在外膜中本身承担桥联与稳定作用。在渗入和内化阶段，$\ce{Ca^{2+}}$ 分布发生明显重排，膜中心区域也出现增强信号（约 $4.0~\mathrm{kg\cdot m^{-3}}$），对应桥联环境重构。这是第三条验证线索：如果四阶段划分是真实的，那么离子分布和聚合物-膜接触应该在每个阶段呈现不同的特征模式。

图4：四阶段中钙离子沿膜法向的密度重排。

• 图4A–图4D对应接近、黏附、渗入、内化阶段下的 $\ce{Ca^{2+}}$ 分布变化。
• 颜色说明：不同颜色曲线对应不同聚合物体系，横轴为膜法向坐标

$\ce{Ca^{2+}}$ 密度分布呈现两个主峰：一个位于 LPS 区域（外膜外表面，距膜中心约 2.5 nm），另一个位于磷脂极性头基区域（外膜内表面，约 -2.5 nm）。当聚合物推进到深层时，$\ce{Ca^{2+}}$ 沿膜重新分布并渗透进入膜中心。本文明确指出，聚合物内化会导致 $\ce{Ca^{2+}}$ 离子沿膜重新分布，而这些离子通过与 LPS 和 DPPG 分子的有利相互作用得到维持。

接触分析揭示三类聚合物的不同相互作用模式

接触数统计（截断距离 0.6 nm）进一步揭示了聚合物与膜组分的特异性相互作用，呈现出三种截然不同的模式。

• 阴离子聚合物 pSPMA 与 $\ce{Ca^{2+}}$ 接触数高达 56.6 ± 1.4，与 LPS 接触数达 96.8 ± 9.3，表明它拖拽了这些离子向内推进。
• 阳离子聚合物 pDMAEMA 和 pMETAC 与 $\ce{Ca^{2+}}$ 几乎没有接触，但与 LPS 分子保持大量接触，渗入阶段分别达 205.6 ± 15.5 和 137.5 ± 5.4，解释了它们为什么更容易在界面站住脚。

本文明确观察到阳离子聚合物转运过程中 LPS 分子被从膜外叶拖向内叶，这一点与接触数统计高度一致：所有聚合物与磷脂（DPPE 和 DPPG）的接触数均为 0，说明转运过程中只有 LPS 分子被聚合物从外叶拖向内叶。pDMAEMA 甚至携带吸附的 LPS 到达膜内介质。

pSPMA 体系还表现出 $\ce{Na^+}$ 在 LPS 叶中的优先积累，这是由阴离子聚合物骨架携带配位 $\ce{Na^+}$ 离子驱动的。

这里有一个容易忽略的细节：原文观察到，平均意义上 $\ce{Na^+}$ 和 $\ce{Cl^-}$ 不会穿膜扩散，后半程的主角是聚合物、LPS 和桥联离子的协同重排。

后半程阻力主要来自三部分：

1. 空间位阻效应。整个聚合物结构，包括其溶剂化壳，在膜内产生空间位阻。
2. LPS 拖拽。阳离子聚合物会携带 LPS 分子从外叶拖向内叶。
3. 局部瞬态膜缺陷。这会导致局部变薄和水渗透。

从设计角度看，这意味着增加电荷主要改善前半程（增强与LPS的静电吸附），但后半程推进仍面临上述三个挑战。

这里的核心物理图像是：界面吸附的热力学有利性不等于向内推进的动力学可行性。自由能计算允许我们将热力学有利过程（如吸附）与动力学受限过程（如渗入）分离开来，这在实验上很难区分。

结果四：水合层是两性离子与阴离子体系的重要缓冲器

为了理解为什么有些聚合物更容易“穿透”膜，需要回答一个基本问题：聚合物在推进过程中会失去多少水合壳层？这是第四条验证线索：如果两性离子聚合物真的“不愿意黏附”，那么它们应该在整个过程中保持更高的水合。

本文用两种互补的方法来回答这个问题：

• 图5给出径向分布函数 $g(r)$，定性分析“水还围在链周围有多紧密”。
• 图6给出配位数 CN，定量统计 4.0 nm 范围内有多少个水珠。

图5：四阶段中聚合物与水珠的径向分布函数。这里的径向分布函数记作 $g(r)$。

• 图5A–图5D对应接近、黏附、渗入、内化阶段的 polymer–TW 统计。
• 纵轴是 $g(r)$ 强度，反映水分子在聚合物周围的概率密度；横轴是与聚合物质心的距离。
• 关键信息：峰高降低对应去溶剂化，峰高保持对应水合壳保存。

图5显示了一个清晰的对比：

• pSPMA 在渗入和内化阶段，即图5C–5D，仍然保持较高的 $g(r)$ 峰，说明水合壳保存得更好。
• pDMAEMA 和 pMETAC 在后期阶段的峰强明显降低，说明去溶剂化更严重。

本文明确指出：pSPMA的强水合暗示其与OM的相互作用较弱，这与阴离子聚合物和LPS之间的静电排斥一致。

图6：四类聚合物在跨膜转运各阶段的水合配位数。CN 统计四种聚合物在 4.0 nm 径向范围内水珠（TW 珠，MARTINI 3 力场）的配位数，反映聚合物在转运各阶段的水合壳层稳定性。每个数值代表三次独立重复的平均值。

• 配位数的定义：在 4.0 nm 径向距离内的水珠数量，用来定量分析聚合物周围的水合结构。
• 四个阶段对比：依次展示接近、黏附、渗入、内化阶段的配位数变化。
• 颜色说明：不同颜色条形代表四种聚合物，高度反映配位数大小。

在 permeation 阶段，四类体系的配位数下降比例如下，显示了去溶剂化差异：

定量结论：四类聚合物在推进过程中都会失去一部分水合壳层，但 pSPMA 和 pMEDSAH 丢得更少。这说明它们在进入膜内时更不愿意完全脱水，也更不愿意和膜内部环境形成紧密耦合。这与它们较弱的深层耦合行为是一致的。

图5和图6的关系是：图5负责定性分析，看 $g(r)$ 峰高变化，直观感受“水合壳有多紧密”；图6负责定量统计，给出具体的配位数值。两者结合说明，pMEDSAH 和 pSPMA 在后期更容易保留一层含水外壳，这是它们抗污行为的关键机制。

到这里，正文已经把“膜怎么变”、“离子怎么变”、“水怎么变”讲清楚了。下一步要问的，就是这些结构变化最后会不会在受力曲线和自由能曲线上留下同样的排序。如果会，前面的结构解释才算真正闭环。

结果五：从受力到自由能的证据闭环

前面的 Figure 3–6 已经说明，推进过程会伴随膜重排、离子重排和去溶剂化。这是最后一条验证线索：这些结构变化，最后能不能在力学和自由能上闭合成一个一致的解释？如果四阶段划分是真实的，那么 SMD 受力曲线和 US 自由能曲线应该在同样的位置出现信号转折——这就是“证据闭环”。

图7：SMD 模拟中聚合物跨膜转运的受力时序曲线。展示四种聚合物在 320 ns SMD 模拟中跨膜转运时的受力 $F$ 随时间变化。每条曲线代表三次独立重复的平均值。聚合物从膜左侧接近并推向右侧。

• 受力曲线特征：不同颜色曲线对应四种聚合物，峰值位置反映各阶段的动力学阻力。
• 阶段分界：受力增长对应黏附建立，膜中心附近的峰值对应渗入和内化的势垒。
• 颜色说明：四种聚合物分别用不同颜色曲线表示，曲线高度反映该时刻施加的力大小。

图7的核心结论：最高力峰位于膜中心，说明 permeation 是阻力最大的步骤。阳离子聚合物（pDMAEMA 和 pMETAC）的力峰更高，说明它们与膜相互作用更强、推进时需要更大的力。相比之下，两性离子和阴离子聚合物（pMEDSAH 和 pSPMA）的力峰更低更窄，反映它们与膜的亲和力更弱——因为“不怎么粘”，所以通过更快，没有长时间持续的相互作用。

这里有一个容易混淆的细节：pMETAC 的黏附力（292.9）比 pDMAEMA（236.6）更高，但最大受力却相反（pMETAC 1855.7 < pDMAEMA 1940.0）。这并不矛盾——$F_{\mathrm{adh}}$ 反映的是“界面一旦接触，谁更容易被膜抓住”，而 $F_{\max}$ 反映的是“推进过程中哪里最难”。pMETAC 虽然更容易被界面捕获，但后续推进相对顺畅；pDMAEMA 虽然界面捕获稍弱，但一旦深入膜内，需要克服更大的阻力才能继续推进。这说明界面吸附和深层推进是两个独立的物理过程。

图8：聚合物向细菌外膜中心推进的势函数曲线。PMF 展示四种聚合物向细菌外膜（OM）中心转运过程的势函数。聚合物从左侧接近膜，沿反应坐标向右侧推进，黑色虚线标示细菌外膜中心位置 $z = 0$。

• 自由能极小值：反映黏附稳定性，阳离子聚合物在 LPS 内核区域约 $z = 2.0$ nm 处出现更深极小值。
• 中心势垒：膜中心区域的能量抬升反映渗入和内化阶段的自由能势垒。

图8的核心结论：阳离子聚合物在 LPS 内核区（z ≈ 2.0 nm）有更深的自由能极小值，说明它们与带负电的糖残基（KDO 和 HEP）的静电相互作用更强。相比之下，pMEDSAH 和 pSPMA 的极小值更靠外（z ≈ 2.7 nm），且数值更接近零，说明它们的黏附更不利。从黏附态到膜中心，所有体系都需要克服显著的自由能势垒（$\Delta G^{\ddagger}$），证实 permeation 是动力学不利步骤——尤其 pMEDSAH 的势垒高达 266.4 kJ/mol。

这里需要强调一个常被误解的点：图7的受力曲线和图8的 PMF 曲线是互补的，而不是重复。SMD 给出的是“沿着特定路径推进时遇到多大的阻力”，而 US/PMF 给出的是“去掉外力后，系统本身的热力学偏好”。两者结合才能完整回答“能不能继续往里走”：图7说动力学上多费力，图8说热力学上多不利——如果两个都说“很难”，那才是真的很难。

统一指标表：把排序和设计建议放在一张图景里

指标解释

• $F_{\mathrm{adh}}$（黏附力）：聚合物到达膜表面时（反应坐标约4.4 nm处）的受力，反映界面捕获的难易程度。数值越大说明越容易被膜“抓住”。
• $F_{\max}$（最大力）：跨膜全程中的最大受力值，通常出现在膜中心附近，反映 permeation 阶段的动力学阻力。
• $\Delta G_{\mathrm{adh}}$（黏附自由能）：聚合物从体相到最稳定位置的自由能变化（PMF极小值），负值越大说明热力学上越有利于黏附。
• $\Delta G^{\ddagger}_{\mathrm{per}}$（渗入势垒，per = permeation）：PMF极小值与膜中心（z=0）自由能的差值——像从山脚爬到山顶的高度差，数值越大说明越难向内推进。

数据解读

这张表把全文最重要的排序浓缩在一起。pMETAC 和 pDMAEMA 的 $\Delta G_{\mathrm{adh}}$ 更负（-761.3 和 -585.6 kJ/mol），说明它们更容易在界面建立稳定吸附——这与阳离子-LPS 的静电匹配一致。但 $\Delta G^{\ddagger}_{\mathrm{per}}$ 呈现不同排序：pMEDSAH 的渗入势垒高达 266.4 kJ/mol，远高于其他体系，说明它虽然能到达界面，但很难继续深入。

这里有一个常被误解的细节：pMEDSAH 的”高势垒”不是缺陷，而是抗污优势。它的 $\Delta G_{\mathrm{adh}}$ 不够负（-440.4 kJ/mol），说明界面吸附不强；$\Delta G^{\ddagger}_{\mathrm{per}}$ 很高（266.4 kJ/mol），说明即便吸附了也不容易穿透。这种双重保守正是抗污材料需要的特性：既不愿意黏上去，就算黏上了也扒拉不进膜内。相比之下，阳离子聚合物虽然界面吸附更强，但渗入势垒相对较低（75.9 和 82.8 kJ/mol），说明它们一旦被膜”抓住”，后续推进反而相对容易——这对抗菌有利，但对膜扰动风险更高。

把这张表和前面的 Figure 7、Figure 8 合起来看，本文其实是在反复强调同一个结论：界面捕获能力和向内推进能力需要分开判断。

四阶段如何被交叉验证

本文最重要的方法论贡献，是用五类独立观测量交叉验证了四阶段的边界：

为什么这个交叉验证重要：

• 单一指标容易误判。如果只看受力曲线，可能认为“力大 = 难推进”；但结合 PMF 会发现，力大可能是因为相互作用强（热力学有利），也可能是因为势垒高（动力学不利）。
• 五类观测量互相印证。构象快照告诉你“聚合物在哪”，受力告诉你“推进多费力”，PMF 告诉你“热力学偏好什么”，水合告诉你“为什么两性离子不粘”，离子和接触告诉你“静电匹配在哪发生”。五条线索指向同一个四阶段边界，这才是可信的机制图景。
• 拆解了“界面吸附”和“向内推进”。黏附阶段的 PMF 极小值说明界面热力学有利，但渗入阶段的 PMF 势垒说明继续推进动力学不利——这两个过程的物理驱动力不同，不能用单一电荷指标覆盖。

设计建议的物理基础

结果逻辑图：从观测到结论

graph TB
    A("观测1：阳离子体系黏附更强") --> D("中间解释：LPS界面静电匹配更强")
    B("观测2：两性/阴离子体系水合保持更高") --> E("中间解释：水合壳层抑制深耦合")
    C("观测3：PMF显示中心区仍有势垒") --> F("中间解释：推进受结构重排与脱水代价限制")
    D --> G("结论：高黏附不等于高穿透")
    E --> G
    F --> G
    G --> H("设计建议：把界面捕获与深层推进分开优化")

SI 数据如何增强正文可信度

SI 提供主结论的统计支撑：

• Table S1 给出四类聚合物净电荷与体积差异，为后续力学排序提供物理背景。
• Table S3 给出 1.0 μs 自由链的 RMSD 与 Rg，证明比较是在已平衡链构象上进行。
• Table S4 给出外膜面积与厚度收敛指标，证明膜基线结构可靠。
• Table S5 给出各阶段 CN 绝对值，避免只看百分比造成误判。
• FigS5 和 FigS6 给出 US 窗口重叠与分块分析，支撑 PMF 收敛。

图S6：PMF 分块分析与收敛性检查

• 图S6把每条 PMF 轨迹按采样时间分成多个 block，分别重建自由能曲线，用来判断不同时间块之间的轮廓是否一致。
• 读图重点是不同 block 的极小值位置、中心势垒高度和整体轮廓是否基本重合。
• 如果这些 block 曲线彼此接近，就说明 Figure 8 里的 PMF 不是某个短时间窗口偶然得到的结果；再结合 SI Figure S5 的窗口重叠情况，才能较有把握地说明 US 采样已经达到可接受收敛。

关键结论与批判性总结

核心结论

• 这篇工作最重要的结论，是把“界面吸附”和“向内推进”明确拆成两个物理问题。阳离子聚合物更容易在外膜表面建立有利黏附，但这并不等于更容易完成后半程渗入。
• 四阶段路径之所以有解释力，是因为构象、受力、自由能、水合和离子重排彼此能对上。文章把接近、黏附、渗入、内化这四步做成了可交叉验证的过程图谱。
• 两性离子和阴离子体系的高水合保持，是它们偏向抗污而非深层扰动的关键原因。这一点在 Figure 5、Figure 6 和 Table S5 中都有相互支撑的定量证据。
• 阳离子链推进过程中伴随 LPS 拖拽、$\ce{Ca^{2+}}$ 重排和局部 nanopore 形成，说明后半程的代价不只是聚合物自己进入膜内，还包括外膜组分被一起重构。

局限与边界

• 自由链不等于真实刷层。本文研究的是自由聚合物链与 OM 的相互作用，而不是显式接枝、显式高密度刷层，因此更接近“刷层末端链段如何与膜相遇”的机制上限，而不是完整表面体系的直接数值预言。
• 膜模型仍然是受控简化体系。虽然外膜做成了 LPS 外层 + DPPE/DPPG 内层的不对称模型，但真实菌膜中的蛋白、拥挤效应、剪切环境和多组分竞争吸附都还没进入模型。
• nanopore 与 LPS 拖拽主要来自模拟轨迹观察。这部分可以作为机制候选或结构后果来讨论，但还不能直接当作已经被实验独立证实的普适结论。
• 粗粒化力场的精度仍有限。MARTINI 3 能高效覆盖大体系和长时间尺度，但对氢键、局部构象和精细水合结构的描述不如全原子模型，某些定量值更适合被看作趋势判断。
• 动力学时间尺度还缺少直接实验标定。PMF 给出了后半程势垒，但真实渗入速率还受扩散系数和动力学前因子控制，这部分仍需和单细胞力谱、QCM、AFM 或荧光示踪实验对接。

后续方向

• 扩展到多链刷层体系：研究接枝密度、链长分布和刷层厚度对界面吸附与渗入行为的影响，建立从单链到刷层的跨尺度模型。
• 结合全原子模拟细化关键步骤：对黏附和渗入的关键构象用全原子MD进行精细化模拟，验证粗粒化模型的定量准确性，尤其是水合结构和氢键网络的细节。
• 拓展到其他细菌外膜组成：本文使用的是rough LPS（缺少O-抗原），而完整的smooth LPS具有更长的多糖链，可能显著改变聚合物的界面识别和渗入路径，需要系统比较不同LPS类型的影响。
• 与实验数据定量对接：结合单细胞力谱、石英晶体微天平（QCM）和原子力显微镜（AFM）等实验手段，验证计算预测的黏附力、黏附自由能和渗入势垒，建立计算-实验闭环验证体系。

本文结果表明，界面黏附强弱和后续渗入难度并不是同一个维度。approach 与 adhesion 可自发发生，而 permeation 与 internalization 仍受较大自由能势垒限制，因此材料设计不能只看单一电荷指标。