CM15抗菌肽与细菌模型膜GROMACS-MD模拟

Zaeifi D. et al. | Iran. J. Biotechnol. 2023, 21(2): e3344 | GROMACS + CHARMM36

一、研究背景与工程问题

抗菌肽（AMP）是应对细菌耐药性的重要方向，CM15（序列KWKLFKKIGAVLKVL）作为阳离埋线性肽，对革兰氏阳性菌和阴性菌均显示广谱抑制活性。但其穿透细菌膜的分子机制仍不明确——到底是桶晶模型、环形孔还是地毯模型？Zaeifi等人基于GROMACS 2016.3 + CHARMM36力场，对CM15与S. aureus（金黄色葡萄球菌）和E. coli（大肠杆菌）模型膜双层分别进行了120ns的全原子分子动力学模拟，您想知道结果怎么样吗？让我们一起来看。

图1 CM15在E. coli（A）和S. aureus（B）膜系统120ns结构演变——红色为脂质磷原子、棕色为肽链，肽段仅驻留在上层小叶

二、仿真方法与系统构建

两套膜系统均通过CHARMM-GUI构建：E. coli膜含PG:PE:CL=12:62:4、5330个TIP3水分子，表面积25.27nm²；S. aureus膜含PG:CL=46:32、8786个水分子，表面积35.94nm²。CM15（PDB: 2JMY）初始位置在膜上方30Å，在水中自由扩散后以无偏倚的取向接触膜双层。模拟参数：303.15K Langevin恒温器、Nose-Hoover热浴、1bar Parrinello-Rahman半同孔压力耦合（tau-p=2ps）、2fs时间步长、LINCS约束氢键、PME计算长程静电（网格间距1.6Å）、短程截断12Å（开关10Å切换）。系统加NaCl至0.1M中和。能量最小化后，膨文原子层平均面积约49.5Ų（E. coli）和58.7Ų（S. aureus），起始阶段均下降约5%。

图2 CM15倾斜角分布、质心距及与脂质接触原子数——N端先接触膜，LYS-6是两体系中接触最多的残基

三、核心发现与工程认识

CM15在两套系统中均以N端先接触膜，倾斜角接近平行于膜法线，但始终未穿透磷原子边界，也未观察到孔洞形成。这与之前POPG:POPC等体系中的插入结果形成对比，说明心碷脂（Cardiolipin/TOCL）的存在改变了CM15的穿透行为。关键接触残基：**LYS-1、LYS-6、VAL-14**，其中LYS-6在两体系中均为接触最多残基。S. aureus体系中膨文柔韧性更高（RMSD从0.29nm升至0.41nm）75ns处），自由能景观显示其稳定性高于E. coli。CM15在S. aureus中倾向于α-螺旋和π-螺旋、E. coli中为turn构像，提示膜脂组成导向了不同折叠模式。

图3 CM15骨架RMSD/RMSF/RASA分析——RMSD从0.29nm升至0.41nm（75ns），S. aureus体系中柔韧性更高

盐桥与氢键统计进一步揭示：S. aureus体系中盐桥分别占POPG 27%和TOCL1 26.2%，而E. coli中POPE单独贡献53.3%盐桥、但氢键仅POPE 6.2%和TOCL1 6.62%。这种差异表明：**TOCL心碷脂的存在大幅增强了S. aureus中的静电稳定化，而PE头基在E. coli中起主导作用**。静电相互作用能在z≥10Å处约为范德华力的3倍，证实阳离埋肽与阴离埋膜脂之间的长程静电力是驱动结合的主要因素。膜厚度和水分子数密度剖面进一步证实PG:TOCL比例越高、结构扰动越大——这对工程上设计针对特定细菌膜脂组成的抗菌肽有直接指导意义。

图4 自由能景观（FEL）与静电/范德华相互作用能——S. aureus稳定性更高，静电能在z≥10Å为范德华约3倍

图5 膜上层小叶脂质面积涨落、厚度涨落与数密度剖面——CM15未能穿透磷原子边界，仅引起局域扰动

四、我们提供的服务

① GROMACS分子动力学模拟服务——蛋白质-膜相互作用、抗菌肽穿膜机制、脂质双层自组装等全原子MD仿真 ② CHARMM-GUI膜系统构建与参数化——基于实验脂质组成的细菌膜模型构建、NaCl中和与平衡 ③ 抗菌肽-膜相互作用计算分析——RMSD/RMSF构像波动、RASA溶剂可及性、FEL自由能景观、静电/范德华力、盐桥/氢键统计 ④ 基于MD模拟的抗菌肽分子设计与优化——序列突变/残基替换、脂质选择性评估、穿膜机制预测 ⑤ 计算生物学全套方案咨询与项目合作 如有需要，欢迎通过公众号“320科技工作室”联系我们。