《Journal of Chemical Information and Modeling》：Synergistic Effects of pH and Temperature on Dengue Virus Envelope Dimers: Insights from Microsecond Molecular Dynamics Simulations
编辑推荐：
登革热发病机制依赖于低pH诱导的包膜糖蛋白（pE）构象重排。然而，pH与温度如何协同破坏pE二聚体形式——这一使融合活性三聚体形成得以启动的基础步骤——仍不清楚。研究人员针对医学相关性最高的登革血清型DENV-2与DENV-3开展研究。采用CHARMM36m力
登革热发病机制依赖于低pH诱导的包膜糖蛋白（pE）构象重排。然而，pH与温度如何协同破坏pE二聚体形式——这一使融合活性三聚体形成得以启动的基础步骤——仍不清楚。研究人员针对医学相关性最高的登革血清型DENV-2与DENV-3开展研究。采用CHARMM36m力场进行微秒尺度全原子分子动力学（MD）模拟，在28 °C、37 °C、40 °C及pH 5、6、7条件下开展，并辅以DENV-2在37 °C/pH 5、DENV-3在28 °C/pH 5及37 °C/pH 7条件下的扩展粗粒化分子动力学模拟。模拟表明，温度由28 °C升至37 °C倾向于在中性pH下扩展构象空间，可能通过调控界面相互作用网络削弱二聚体完整性。酸化以血清型特异方式放大该效应：模拟中DENV-2需pH 5，而DENV-3在pH 6即响应。距离与主成分分析揭示了非对称解离路径，称为“补偿性拥抱”（compensatory embrace），在全原子与粗粒化模型中均观察到。结构分析进一步突出血清型间组氨酸分布差异。此外，37 °C下的相互作用图谱识别出可能作为pH敏感闩锁的候选界面接触，在pH 7维持稳定而酸化时断裂。综上，结果提示可能存在物理化学触发因素的层级，促成二聚体向融合胜任单体转化。这些界面热点与“补偿性拥抱”构成结构知情的靶点，经天然病毒粒子条件验证后，可指导抗病毒药物与免疫原理性设计，凸显长时间尺度MD模拟作为药物发现中假设生成工具的价值。
研究背景与立项依据
登革病毒（DENV）是全球最普遍的虫媒病毒，2024年病例超1400万、死亡约9500人，重症与抗体依赖增强（ADE）密切相关。包膜糖蛋白（pE）介导受体结合与内体低pH膜融合，成熟病毒表面pE以反向平行同源二聚体平铺，低pH下二聚体失稳暴露域II融合环（FL）插入靶膜，进而形成融合三聚体。温度同样关键：28 °C（蚊媒）、37 °C（人体生理）、40 °C（发热）调节病毒粒子“呼吸”、表位暴露及二聚体稳定性，且具血清型依赖性。现有结构/计算研究多聚焦单血清型或孤立酸化，pH-温度联合效应未知。疫苗需应对四血清型并诱导平衡免疫，跨原聚体桥接的pE二聚体四级表位（EDE）是广谱中和热点，但二聚体稳态如何被物理化学变量协同破坏尚缺动态证据。该研究即在此背景下，用微秒全原子与粗粒化MD、主成分分析（PCA）与自由能景观（FEL）映射，考察DENV-2（PDB 1OAN）与DENV-3（PDB 7A3S）胞外域（ecpE）二聚体在pH 5/6/7 × 28/37/40 °C矩阵下的行为。
主要技术方法概述
研究人员从PDB获取DENV-2（1OAN）与DENV-3（7A3S）ecpE二聚体坐标，用PROPKA3分配pH 5/6/7组氨酸质子态，CHARMM-GUI重构缺失环，CHARMM36m力场建拓扑，TIP3P水盒150 mM离子强度，NAMD跑1 μs/体系全原子NPT生产；选取D2T37pH5、D3T28pH5、D3T37pH7做Martini 3粗粒化映射，GROMACS+PLUMED跑5副本×30 μs。轨迹用VMD/MDAnalysis处理，算Cα RMSD、界面RMSD（≤8 ?跨界面）、RMSF、质心距dM→M（60 ?解离阈）、PCA取PC1/PC2建FEL（ΔG=?kBT ln(P/Pmax)）、GetContacts分类氢键/π-阳离子/盐桥、粗粒化下算FL-SASA与FL-托槽距。统计上单轨迹非独立重复，仅作描述性趋势。
研究结果
3.1 节 DENV-2/3 pE胞外域结构与组氨酸分布比较：TM-align叠加TM-score 0.95、RMSD 1.93 ?，整体折叠保守但DI/DIII环区可变。DENV-2界面旁组氨酸His144/His149/His158/His317，DENV-3为His144/His149/His315/His345，His144与His317（D3为His315）文献定为DI-DIII界面主pH传感器；次级组氨酸可能微调血清型特异激活阈。静态PROPKA3质子态未含动态滴定（CpHMD），结论需谨慎。
3.2 节 pH与温度变异下二聚体稳定性：全局Cα RMSD显示DENV-2在pH 5/37 °C最高（5.6±1.4 ?），DENV-3在pH 6/37 °C最高（6.3±1.2 ?），印证DENV-2需更强酸化、DENV-3在 milder 酸度即响应。中性pH下升温28→37 °C两血清型均拓宽构象波动。
3.3 节 局部界面扰动与残基柔性：界面RMSD在DENV-2 pH5/37 °C达5.0±1.3 ?，DENV-3 pH6/37 °C达5.8±1.9 ?。RMSF显示一原体柔性增、另一原体降的链间不对称性，预示非对称解离。
3.4 节 pH与温度在1 μs内触发ecpE二聚体分离事件：dM→M直方中三轨迹超60 ?——D2pH5/37 °C、D3pH5/28 °C、D3pH7/37 °C；快照见“补偿性拥抱”：一侧DI/DIII撤、对侧DI推进补接触，总接触数不降反掩部分分离。粗粒化5副本中D3T28pH5与D3T37pH7各2副本完全解离，D2T37pH5未达阈但4/5副本FL-托槽距不对称；FL-SASA佐证。
3.5 节 pH与温度对ecpE二聚体本质动力学影响：中性pH升温28→37 °C，PCA云由单簇变双簇（D2）或三簇（D3），FEL由~2.3 kBT单盆变0.8 kBT双盆（D2）或<0.5 kBT浅盆（D3）。PC1为DIII转动、PC2为DI绕DI-DIII界面非对称摆，支撑拥抱机制。pH6呈“锁住”单盆，pH5/37 °C盆互联且DII向溶剂敞开，预示DIII释放。
3.6 节 温度对DENV pE二聚体界面分子相互作用的影响：跨温度共有接触中，28→37 °C蓝标接触频降、40 °C部分回弹。氢键（hbbb/hbsb/hbss）、π-阳离子（pc）、盐桥（sb）量化显示中性→酸化时界面盐桥/氢键网络被温度放大削弱；37 °C下pH7稳定、酸化断裂的候选闩锁接触被识别为血清型特异热点。
讨论与结论翻译
讨论指出，模拟提示酸pH与升温重塑ecpE构象景观，改变其介导入胞与免疫逃逸能力，并暗示小分子抑制融合活性的作用靶点。组氨酸分布差异提供血清型pH阈微调的结构基础；但静态质子态与可溶ectodomain脱离M蛋白/脂双层，限制向完整病毒粒子外推。粗粒化仅验证选中条件，非全矩阵统计确证。结论部分原文表述为：结果提示存在潜在物理化学触发因素层级，可能促成二聚体向融合胜任单体转化；所识界面热点与“补偿性拥抱”构成结构知情靶点，待天然病毒粒子条件验证后，可指导抗病毒与免疫原理性设计，并凸显微秒MD作为药物发现假设生成工具的价值。
（全文约1520汉字）
