利用计算方法探究Alyteserin-1c在革兰氏阳性与革兰氏阴性细菌膜模型中作用机制的研究

《MicrobiologyOpen》：Exploring the Mechanism of Alyteserin-1c in Gram-Positive and Gram-Negative Bacterial Membrane Models Using a Computational Approach

　　

摘要：Alyteserin-1c是一种含23个氨基酸的阳离子抗菌肽（Antimicrobial Peptide, AMP），对革兰氏阴性菌的活性高于革兰氏阳性菌（Conlon et al. 2009），但其差异膜选择性的分子基础尚不完全清楚。研究人员采用全原子分子动力学（Molecular Dynamics, MD）模拟（每体系200 ns × 3次独立重复），考察了预组装的Alyteserin-1c六聚体与分别代表革兰氏阳性（75% PG/25% PE）和革兰氏阴性（75% PE/25% PG）细菌内膜的简化双层模型间的相互作用。六聚体组装体在所有模拟中均保持结构稳定性和大量α-螺旋含量，寡聚体凝聚力主要与残基L10、L13、V14、I17和F6的疏水堆积有关。阳离子残基K7和K15保持暴露于脂质并近乎持续接触脂质头部基团，提示其为主要膜锚定位点。比较膜分析表明革兰氏阴性膜模型具有更大的侧向脂质流动性、更低的双层密度及更强的局部酰基链扰动，与较低肽浓度下膜更易受影响相符。螺旋取向和孔几何分析不符合稳定桶板孔（barrel-stave）跨膜孔形成特征，而支持界面或瞬时环形孔（toroidal-like）膜扰动。革兰氏阳性体系中观察到瞬时单列水分子穿线，而革兰氏阴性膜中的扰动主要为表面局域化。这些发现为Alyteserin-1c的膜选择性提供了计算推导的机制框架，可为理性设计膜靶向抗菌肽提供参考。鉴于所用简化双层模型、预组装寡聚假设及可达模拟时间尺度，所提出机制应视为假说生成性质而非定论。

本研究发表于《MicrobiologyOpen》。

论文解读

Alyteserin-1c是从助产蟾蜍（Alytes obstetricans）皮肤分泌物中分离得到的一种23氨基酸阳离子α-螺旋抗菌肽（Antimicrobial Peptide, AMP），实验显示其对大肠杆菌（革兰氏阴性）的最小抑菌浓度（Minimum Inhibitory Concentration, MIC）为25 μM，对金黄色葡萄球菌（革兰氏阳性）MIC为100 μM，且对人红细胞溶血活性低（LC₅₀=145 μM）。传统抗生素因抗菌耐药（Antimicrobial Resistance, AMR）蔓延而效力下降，膜靶向AMP因其不依赖特定受体、通过膜溶解作用杀菌而受关注。两栖动物来源阳离子AMP因广谱活性和低真核细胞毒性成为研究热点，但多数AMP（包括alyteserin家族）对革兰氏阳/阴性菌存在2～4倍MIC差异，其选择性分子机制——源于差异结合亲和力、不同孔结构还是膜破坏动力学差异——因实验技术难以在原子分辨率下捕捉瞬态孔结构而仍不清楚。先前对Alyteserin-1c的单体和阳离子类似物模拟仅限500 ps且未涉及寡聚体协作成孔，而孔形成AMP通常需要四聚体至八聚体寡聚组装。因此研究人员通过全原子MD模拟考察预组装Alyteserin-1c六聚体在仿革兰氏阳性与阴性细菌内膜模型中的作用机制，以阐明其膜选择性的分子决定因素。

为开展研究，研究人员主要采用以下关键技术方法：从PDB（ID: 2L5R）获取Alyteserin-1c单体结构并经AMBER ff14SB力场能量最小化；用ClusPro 2.0平台通过PIPER刚体对接逐步构建二聚至七聚体寡聚体，并用ProFlex进行氢键稀释刚性分析筛选最适寡聚态；用CHARMM-GUI Membrane Builder构建革兰氏阳性内膜模型（75%磷脂酰甘油Phosphatidylglycerol, PG / 25%磷脂酰乙醇胺Phosphatidylethanolamine, PE）和革兰氏阴性内膜模型（75% PE / 25% PG），脂质采用AMBER LIPID17力场参数；六聚体经OPM服务器定向后插入双层中心，C端酰胺化、N端游离氨基，TIP3P水盒子、150 mM Na⁺/Cl^?中和，GROMACS 2024.4进行NPT系综（310 K, 1.0 bar）200 ns × 3次独立重复模拟；通过gmx energy分析体系能量，gmx rms/rmsd/gyrate及DSSP分析六聚体结构稳定性，自定义Python脚本（MDAnalysis）分析单体间氢键/盐桥/接触图、肽-脂质氢键及残基水平接触占有率，gmx density/MEMBPLUGIN计算膜质量密度及厚度，脂质序参数（Segmental Deuterium Order Parameter, S_CD）、倾斜角和Einstein关系求侧向扩散系数（Lateral Diffusion Coefficient, D_xy），HOLE2分析孔半径及水动力学（水占据数、偶极取向），自定义脚本计算螺旋倾斜角（Helix Tilt Angle, θ）。

3.1 Screening of Oligomeric States（寡聚态筛选）

通过ClusPro 2.0逐步对接获得二聚至七聚体，六聚体呈对称束状几何具潜在膜插入与破坏能力，七聚体外周无序被排除；结合刚性分析六聚体氢键指数最高（105）且破坏能最大（?4.326 kcal·mol^?1），研究人员选定六聚体为后续膜模拟用寡聚态。

3.2 Rigidity and Flexibility Analysis of Alyteserin-1c Oligomers（Alyteserin-1c寡聚体的刚性与柔性分析）

ProFlex氢键稀释分析显示从单体到六聚体分子内氢键数由7增至105，所需破坏能由?2.603增至?4.326 kcal·mol^?1，六聚体具最广泛氢键网络与最高机械稳定性，支持其为最适膜插入寡聚态。

3.3 Molecular Dynamics Simulations of Alyteserin-1c Hexamer（Alyteserin-1c六聚体的分子动力学模拟）

每体系3次200 ns独立重复NPT模拟，温度（GN: 309.999±0.001 K; GP: 309.998±0.001 K）与压力（GN: 0.942±0.08 bar; GP: 0.904±0.07 bar）全程稳定收敛，确认肽-膜体系热力学平衡。

3.4 System Energetics and Density（体系能量学与密度）

GN与GP体系势能分别稳定在?1.292±0.002×10⁶kJ/mol与?1.088±0.001×10⁶kJ/mol（差值源于体系原子数不同），Coulomb短程（Coulomb-SR）相互作用能GN（?1.864±0.003×10⁶kJ/mol）绝对值大于GP（?1.721±0.004×10⁶kJ/mol），反映GN含更多水分子及75% PE头区间氢键网络增强静电稳定。GP双层质量密度（1033.9±0.04 kg/m³）高于GN（1019.0±0.03 kg/m³），与PG富集膜紧密堆积一致。

3.5 Structural Stability of the Hexamer（六聚体结构稳定性）

GN与GP中六聚体骨架RMSD分别为0.387±0.013 nm与0.329±0.020 nm，GN略高反映较松散双层中构象灵活性更大，二者均达平稳平台无去折叠；回转半径（Radius of Gyration, R_g）GN: 1.491±0.018 nm，GP: 1.504±0.012 nm，维持紧凑组装；α-螺旋含量GN: 80.6%±0.8%，GP: 80.3%±0.5%，全程稳定，符合界面吸附carpet/环形孔机制的螺旋保留特征。

3.6 Inter-Monomer Cohesion and Peptide–Membrane Interactions（单体间内聚力与肽-膜相互作用）

单体间氢键（GN: 2.16±0.24; GP: 1.71±0.61/帧）与盐桥极少，六聚体稳定主要靠疏水堆积——L10?L13（95.2%）、V14?L13（95.2%）、V14?L10（90.5%）、I17?L10（85.7%）、I17?F6（81.0%）接触占有率>80%。阳离子残基K3、K7、K15不在界面，保持脂质暴露；残基水平接触占有率K7>98%、K15=100%（两体系），确定为首要静电锚定位点，K3约87%，His20最低（~43%–64%）位于C端溶剂侧。

肽-脂质氢键前80 ns由~16–20升至稳定平台（GN: 26.9±0.7/帧；GP: 28.5±1.5/帧），符合初始静电吸附→稳定嵌入两阶段过程；GP因75% PG具更多氢键供/受体故氢键略多。肽-水氢键由~160降至~145/帧，反映疏水面插入排开水，两体系终值相近，提示差异性源于脂质特异互作而非溶剂化差异。

3.7 Membrane Property Analysis（膜性质分析）

3.7.1 质量密度沿双层层法向呈典型双峰，GP头部区峰值更高（985±2 vs. 967±8 kg/m³），中心疏水区~710 kg/m³无连续水通道，与无稳定跨膜孔一致。

3.7.2 脂质RMSD持续上升属正常侧向扩散，GN终值更高（4.026 nm vs. GP 3.848 nm）表明GN侧向流动性更大；脂质R_g于75 ns内降至平台（GP略高），说明膜在插入六聚体周围压缩平衡。

3.7.3 脂质序参数S_CD：sn-1棕榈酰（PA）链GN plateau（~0.231）>GP（~0.220）；sn-2油酰（OL）链C10处现特征极小值（GN: 0.031; GP: 0.025）验证LIPID17中Δ9顺式双键。GN中S_CD较高源于Lys7/K15强锚定限制邻近脂链构象——与GN整体更高D_xy并不矛盾，属局部有序/全局流动并存。

3.7.4 脂质倾斜角度GP各脂类均高于GN（头基PE: 41.0°±0.3° vs. 39.3°±0.4°；PGR: 42.0°±0.5° vs. 40.8°±0.4°；链PA: 28.3°±0.3° vs. 26.9°±0.2°；OL: 30.1°±0.3° vs. 28.6°±0.1°），与PG富集紧密堆积需倾斜填充空间相符。

3.7.5 双层厚度（P31峰间距）：GN 40.69±0.04 ?略大于GP 40.12±0.04 ?，两体系50 ns内达平衡且数值在PE/PG双层预期范围。

3.7.6 侧向扩散系数D_xy：GN 3.24±0.08×10^?8cm²/s > GP 3.05±0.06×10^?8cm²/s，均低于无肽POPE/POPG膜（~8.5×10^?7cm²/s），证实六聚体插入降低脂质侧向运动，GN整体流动性仍高于GP。

3.8 Helix Tilt Analysis（螺旋倾斜角分析）

六条单体相对于膜法向倾斜角GN: 80.3°±0.9°，GP: 81.1°±1.2°（即相对膜平面约9.7°与8.9°），近平行吸附于界面，不支持需垂直插入（~0°）的经典桶板孔模型，符合carpet/环形孔机制。GN中单体M6倾斜角最低（73.3°）提示稍深插入，GP各单体较均匀（77.0°–84.5°），反映GN较宽松双层允许六聚体构象灵活性更大。

3.9 Pore Radius and Water Dynamics（孔半径与水动力学）

3.9.1 孔半径：HOLE2给出平均最小孔半径GN 0.42±0.05 ?、GP 0.41±0.05 ?，远低于单水分子通透阈值（~1.4 ?），三重复一致，200 ns内未形成稳定跨膜水通透孔。

3.9.2 孔水占据数：GN 1.00±0.25分子/帧，GP 1.15±0.17分子/帧，无持续跨膜水流。

3.9.3 水偶极取向：GP中水偶极穿越全双层（Z=?8~+8 ?）上下叶反向排列（cosθ下叶≈?0.5，上叶≈+0.5），提示六聚体内可瞬时出现单列水穿线（非稳定开孔）；GN中偶极有序主要限于上叶界面（Z=0~+8 ?），为表面局域化膜扰动。研究人员认为GP致密双层约束六聚体为较直立对称排布促瞬时水穿线，GN较流动双层使六聚体位姿多变打断水线连续性——两者均为pre-pore界面扰动态，GN因更强静电锚定与局部脂链扰动可在更低肽浓度起效，与Conlon等报道4倍MIC差异相符。

讨论与结论总结

研究人员在结论中指出：本工作通过全原子MD模拟（每体系3×200 ns，总计1.2 μs）考察Alyteserin-1c六聚体于仿革兰氏阳性（75% PG/25% PE）与革兰氏阴性（75% PE/25% PG）细菌内膜模型中之行为。ProFlex刚性分析确定六聚体为最稳定紧凑寡聚态。插入膜后六聚体骨架RMSD稳定、α-螺旋含量~80%得以保持、回转半径无显著变化；寡聚体内聚力源于L10/L13/V14/I17/F6疏水堆积，K7（>98%）与K15（100%）暴露于脂质并作为主锚定位点与脂质头部基团持续接触。GN模型表现更大全局侧向脂质流动性（更高RMSD与D_xy、更低整体密度）及插入位点附近更强局部脂链有序化（更高S_CD），与GN膜在低肽浓度下更易被破坏相符。螺旋倾斜角~80°（近膜平行）及平均最小孔半径<0.5 ?不支持桶板孔模型，而支持carpet型或界面/瞬时环形孔（toroidal-like）膜破坏机制。水偶极分析显示GP中六聚体被致密 bilayer约束为较对称构型并出现瞬时单列水穿线信号，GN中扰动主要局限于上叶界面——两者均为pre-pore界面态，GN因强静电锚定与有效局部脂链扰动实现更低有效作用浓度，合理解释实验MIC差异。研究局限含：六聚体为真空刚体对接预测非实验验证寡聚态；膜模型为简化PE/PG双层未含革兰氏阴性外膜脂多糖（Lipopolysaccharide, LPS）或革兰氏阳性肽聚糖层。该计算机制框架可为理性设计具增强效价与选择性的膜靶向抗菌肽提供参考。