Hsp90变构激活MD模拟与SDAR定量评估

关键词：Hsp90 | 变构激活 | MD模拟 | DMRS | SDAR

一、研究背景

Hsp90（90kDa热休克蛋白）是细胞蛋白质稳态网络的核心分子伴侣，通过ATP结合和水解驱动构象循环来调节客户蛋白的折叠与成熟。Hsp90功能失调与肿瘤、神经退行性疾病密切相关，使其成为药物设计的焦点靶点。然而，绝大多数Hsp90靶向药物都是ATP竞争性抑制剂——直接封闭催化位点，导致热休克反应和毒副作用。Colombo团队联合UCSF的Agard实验室在Scientific Reports（2016）上发表了一项里程碑式工作：通过全原子分子动力学模拟揭示了一组新设计的变构激活剂如何远程调控Hsp90 ATPase活性——结合位点距离催化中心65Å，却能将酶活性加速达6倍。我们深入复现了该论文的模拟方法学和定量构效关系模型。

图1 Hsp90构象循环及变构配体分子结构：ATP结合驱动NTD二聚化闭合，变构激活剂（AS）结合稳定不对称高能态

二、仿真方法

该研究采用GROMACS软件包和GROMOS96 53A6联合原子力场，以酵母Hsp90的闭合态晶体结构（2CG9.pdb）为起点。先在GLIDE SP模式下将6个苯并呋喃类变构激活剂（化合物1、4、10、12、16、18）及3个香豆素类抑制剂对照分别对接到此前鉴定的变构位点（MD:CTD界面，残基474-487/502-503/591-599/602-603/652-657）。体系用SPC水模型溶剂化、Na+中和电荷，在NVT系综（100ps,300K）和NPT系综（100ps,1bar）分别平衡后，各复合物运行100ns MD模拟（时间步长2fs、PME静电、LINCS约束键长），每种配体做3个独立副本。数据分析在25-100ns区间进行：计算NTD质心距离分布、原聚体主轴扭曲角时间演化、CP残基对距离涨落矩阵（含~2000残基的2000×2000矩阵），以及DMRS（多受体结构对接）——将6个配体分别对接到ATP-only轨迹的10个代表性构象中，计算Ligand Efficiency（LE）后与实验ATPase速率（1.17~4.01 min^-1）做线性回归。

图2 N端结构域质心距离分布：激活剂维持NTD-NTD紧凑闭合（~43Å），抑制剂使分布峰偏移约5Å至开放态

图3 两原聚体主轴相对2CG9晶体结构的角度扭曲时间演化：激活剂使B原聚体（红色）显著扭曲，化合物16例外

三、结果分析

三项核心发现可直接指导工程实践。第一，NTD-NTD距离分布揭示激活剂维持二聚体闭合态（~43Å），抑制剂使峰偏移约5Å至开放态——这5Å的构象选择直接决定催化是否启动。第二，不对称扭曲是激活的标志：激活剂使B原聚体（与配体可变功能基团接触的链）主轴相对2CG9扭曲角显著偏离，而化合物16（含cis-OH的β-L-岩藻糖基团，无法同时桥接E477和D503）几乎不产生不对称性，其ATPase速率也最低（1.17 min^-1）。化合物18的正电荷叔胺同时静电桥接两个酸性侧链，诱导最大扭曲，对应最大激活（4.01 min^-1）。第三，CP差值矩阵提供全局动态指纹：激活剂使NTD-SMD/CTD域间解耦增强（蓝色区扩大），而抑制剂使全蛋白均匀增柔——前者加速构象搜索，后者瓦解活性构象。最实用的是SDAR模型：DMRS对接+LE平均值与实验ATPase速率的相关系数高达R=-0.90，远超单构象对接（R≈-0.10），从工程角度看提供了一个可计算、可复用的筛选流程——先用ATP-only的10个代表性构象做对接池，再算LE排序，无需为每个新配体跑100ns MD。FRET闭合速率的LE相关性也达R=-0.76。

图4 变构配体与Hsp90关键残基E477/D503的氢键网络：trans-OH和带电叔胺桥接诱导不对称应变

图5 配体复合物与ATP-only的残基对距离涨落差值矩阵（CP分析）：蓝色区=激活剂增强的域间解耦，抑制剂均匀增柔

四、我们提供的服务

我们团队在蛋白-配体MD模拟与变构药物设计领域积累了系统性项目经验，可提供以下工程服务：①基于GROMACS/AMBER的全原子蛋白-配体复合物MD模拟全套流程（体系构建-平衡-生产-轨迹分析），支持SPC/TIP3P溶剂模型和GROMOS96/AMBER99SB-ILDN/CHARMM36力场选择。②GLIDE/ADFR多构象对接与DMRS分析：从蛋白MD轨迹提取代表性构象集群（聚类分析），批量对接到变构位点/活性位点，计算LE/MM-GBSA评分矩阵。③基于CP距离涨落矩阵的蛋白内部动力学指纹分析：域间耦合/解耦模式自动识别，配体结合前后的刚性-柔性区块对比可视化。④变构位点从头发现：利用MD轨迹的残基对协调性分析和口袋探测，定位潜在的远端调控位点。⑤SDAR定量模型建立：LE/动力学描述符与实验活性值（ATPase/FRET/IC50）的线性回归与交叉验证，为新配体虚拟筛选提供打分函数。⑥Hsp90及激酶/GPCR/核受体等靶点的变构激活剂/抑制剂理性设计与活性预测。

图6 化合物18在Hsp90 10个代表性构象中的多构象对接：DMRS方法捕捉蛋白-配体动态适应，LE与ATPase速率强相关R=-0.90

图7 化合物18诱导的不对称扭曲：ATP-only态（浅蓝）与18号复合态（红）原聚体全局对齐对比，B链显著偏移

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