HIV-1脂质体全尺寸MD模拟

关键词：MARTINI粗粒化 | 脂质体 | 弯曲刚度 | 脂质翻转 | 有序结构域

一、研究背景

HIV-1作为一种包膜病毒，其脂质双分子层来源于宿主细胞质膜，是维持病毒结构完整性、介导膜融合和感染性的关键组分。脂质组学分析表明，HIV-1脂质囊泡富含胆固醇、鞘磷脂、磷脂酰丝氨酸和磷脂酰乙醇胺等24种脂质物种，且内外叶呈不对称分布。然而，全尺寸病毒脂质体在原子/粗粒化分辨率下的结构、力学和动力学特性仍缺乏系统表征。Bryer A.J.等人（PLOS Computational Biology, 2022）利用MARTINI粗粒化力场和GROMACS大规模并行分子动力学模拟，构建并模拟了直径为150nm的全尺寸不对称HIV-1脂质体模型（含约30万脂质分子、2460万粒子），进行了总计5.2微秒的生产模拟，系统揭示了其厚度分布、弯曲刚度、脂质侧向扩散与翻转以及纳米尺度有序结构域的动态行为。

图1 HIV-1脂质体全尺寸MARTINI粗粒化模型（150nm，24种脂质）（复现结果）

二、仿真方法

本研究采用MARTINI 2.1粗粒化力场，通过packmol初始布种结合"炼金术相变"生长程序构建150nm直径HIV-1脂质体。体系含24种脂质物种（POPS/DLPC/PRPS/PNSM/POPC/POPE/DOPC/PQPE/PGPS/PQPS/DPPC/PIPE/PRPE/CHOL/PAPE/PIPS/DPSM/PAPS/PUPE/DPGS/PUPC/PRPC/PGPE），按HIV-1病毒粒子脂质组学实验数据确定各物种丰度。内外叶不对称分配：负电荷、不饱和及乙醇胺头部基团富集于内叶，鞘磷脂富集于外叶。体系水合后含约2056万MARTINI水粒子、22465个Na+离子和5%抗冻粒子，总计约2459万粒子。模拟使用GROMACS 2018.1，10fs时间步长，Verlet截断方案，反应场静电，Berendsen恒压器（1bar），脂质与溶剂分别用Berendsen恒温器耦合323K。经过1.5微秒初始平衡后进行5.2微秒生产模拟。作为对照组，构建了46×46nm非对称平面膜体系并进行了7微秒模拟。疏水厚度通过泊松盘采样（1024点）计算内外叶C2B碳珠质心距离获得；侧向扩散系数由脂质质心均方位移（MSD）线性拟合得到；跨双层翻转速率基于自主研发的3D空间分类器统计；弯曲刚度通过球谐展开（74度149阶）结合Helfrich连续介质模型计算。

图2 脂质体疏水厚度分析（泊松盘采样1024点，平均26±1A）（复现结果）

图3 脂质跨双层翻转速率与侧向扩散系数（复现结果）

三、结果分析

模拟结果揭示了HIV-1脂质体的多项关键性质：①球形度在整个生产模拟中维持约0.9984，外径从145nm缓慢增至147nm。②疏水厚度分析显示平均值为26±1A，帧间标准差小于1A，厚度波动微小且均匀分布。③跨双层翻转事件显著——不饱和PE/PS头部基团脂质翻转速率>200（10-5事件/ns/分子），而外叶富集的鞘磷脂翻转速率<100，尽管翻转频繁但内外叶不对称组成在5.2微秒内维持稳定（±2脂质分子）。平面膜对照组则完全未观察到翻转事件，表明曲率与脂质翻转存在密切关联。④侧向扩散系数约为10-6 cm2/s量级，内叶扩散显著高于外叶，呈现异质性。⑤外叶观测到直径约10nm的瞬时低迁移率有序结构域，富集胆固醇、鞘磷脂和磷脂酰胆碱，这些"脂筏样"微区持续数纳秒后消散并重新形成。⑥球谐展开分析获得弯曲刚度kc=109kbT（298K），显著低于主相变温度以下的膜刚度值，表明脂质体处于流体相。⑦平面膜压缩模量分析显示外叶（155mN/m）显著高于内叶（103mN/m），与侧向迁移率差异一致。这些结果表明HIV-1脂质体的不对称脂质组成通过耦合曲率与翻转维持稳定性，同时有序结构域可能为包膜蛋白gp41的组织和融合提供动态平台。

图4 外叶脂质低迁移率有序结构域（~10nm，富集胆固醇和鞘磷脂）（复现结果）

图5 HIV-1脂质体弯曲刚度分析（球谐展开，kc=109kbT 298K）（复现结果）

四、我们提供的服务

我们团队在分子动力学模拟领域提供以下特色仿真服务：①基于GROMACS/NAMD的全原子与MARTINI粗粒化生物膜体系MD模拟，涵盖脂质体自组装、膜蛋白-脂质相互作用、跨膜转运与膜融合等过程。②MARTINI力场参数化与体系构建，支持复杂脂质组学数据驱动的非对称膜模型搭建与长时程稳定性评估。③脂质侧向扩散系数、跨双层翻转速率、膜厚度分布与弯曲刚度的定量分析（泊松盘采样、球谐展开、MSD拟合）。④脂质有序结构域（脂筏）识别与动态追踪，包括RMSF迁移率映射与富集组分分析。⑤平面膜压缩模量与弯曲模量的第一性原理计算，支持力场对比与方法学验证。

如有需要，欢迎通过公众号"320科技工作室"联系我们。