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微秒的生产模拟，系统揭示了其厚度分布、弯曲刚度、脂质侧向扩散与翻转以及纳米尺度有序结构域的动态行为。
二、仿真方法
本研究采用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连续介质模型计算。
三、结果分析
模拟结果揭示了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的组织和融合提供动态平台。
