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的全原子分子动力学模拟，您想知道结果怎么样吗？让我们一起来看。
二、仿真方法与系统构建
两套膜系统均通过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%。
三、核心发现与工程认识
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构像，提示膜脂组成导向了不同折叠模式。
盐桥与氢键统计进一步揭示：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比例越高、结构扰动越大——这对工程上设计针对特定细菌膜脂组成的抗菌肽有直接指导意义。
