分子动力学模拟对接：把静态对接结果放进动态溶剂里重新审视

一个激酶抑制剂的对接项目，AutoDock Vina给出的Top 5构象里排名最高的那个，放进显式溶剂跑分子动力学模拟对接，10 ns就解离了——配体从结合口袋里漂出去，在溶剂里自由游荡。排名第二的构象反而稳定住了，50 ns的轨迹中配体RMSD始终在2 Å以内。静态分子对接的最大问题就在这里：它忽略了蛋白质和配体双方的柔性，打分函数在真空环境下做出的判断，在溶剂里不一定成立。这个项目差点提交了错误的优势构象——如果没做分子动力学模拟对接验证，结论就是错的。

静态对接的三个根问题

对接程序给一个固定的蛋白结构、配体旋转几个可旋转键、打分函数一算——排名最高的构象在实验上不一定是对的。问题的根在于三个被忽略的物理事实：蛋白侧链会重新取向以适配配体、活性位点的水分子会被挤出形成去溶剂化代价、配体结合前后的构象能差在静态打分中完全丢失。这三件事在静态框架里无法被正确处理——不是”处理得不好”，是”根本没有被处理”。分子动力学模拟对接把对接后的蛋白-配体复合物泡进显式溶剂里跑一段轨迹，让溶剂和蛋白柔性共同参与结合模式的自然筛选，是验证结合模式稳定性的标准做法。

操作流程与关键判据

操作上，对接给出Top 3-5个构象，分别建好复合物体系、加溶剂盒子、加抗衡离子、做能量最小化、NVT+NPT平衡，然后进入50-100 ns的生产模拟。看三样东西：配体RMSD是否稳定（2-3 Å以内说明没跑飞）、关键氢键的占据率（>60%才算有持续作用，低于30%的氢键在统计上几乎等同于噪声）、结合姿态是否在模拟中自发调整到更深的能量阱。这个项目中，排名第一的构象RMSD在8 ns后突破5 Å并持续攀升——配体已经离开了结合口袋，所谓”打分最高”只是真空环境下的局部最优，在真实溶剂里站不住脚。

MM/PBSA：从定性验证到半定量评估

MM/PBSA（分子力学-泊松-玻尔兹曼表面积）是做结合自由能计算的主流后处理方法。它从MD轨迹中反复提取快照（通常取轨迹尾部20-30 ns中的200-500帧），对每帧分别计算分子力学结合能、极性溶剂化自由能（PB求解）和非极性溶剂化自由能（SASA换算），取统计平均。MM/PBSA的绝对结合自由能和实验值之间的平均偏差在3-5 kcal/mol——不足以做精确排序，但区分强结合（<-8 kcal/mol）、中等结合（-5到-8）和弱结合（>-5）是可靠的。更精确的MM/GBSA方法使用广义Born模型替代PB，收敛更快但精度略低——作为高通量筛选的粗筛工具更合适。这个项目用MM/PBSA计算了五个构象的结合自由能，排名第一的构象ΔG=-3.2 kcal/mol（弱结合），排名第二的ΔG=-9.1 kcal/mol（强结合）——分子动力学模拟对接不仅纠正了定性判断，还给出了半定量的排名修正。

参数选择：内外介电常数的影响不可低估

参数方面，两件事值得注意。igb=2（GB模型的选择）和saltcon=0.15（离子强度150 mM的生理条件）是接近实验条件的默认组合。但PB计算中的内外介电常数设定对结果影响显著——蛋白内部介电常数选1（完全真空）还是2-4（考虑部分极化），结合能绝对值能差出一倍。认定ε_in=1是对蛋白内部极化效应的彻底否定——蛋白质内部的氢键网络和侧链极化对静电环境有真实的响应，ε_in=2是更合理的折中，既不完全忽略极化也不过度放大。这个项目比较了ε_in=1和ε_in=2两种设定，结合能差了6 kcal/mol——参数选择对结论的影响和构象选择一样大。

从验证到决策

回过头看这个对接项目，分子动力学模拟对接的价值不在于”算出更精确的数值”，在于”把虚假排名淘汰掉”。静态打分把真空环境下的局部最优推到了第一位，显式溶剂模拟用物理真实把它拉下来了。对接后不做MD验证就提交结论——这个做法在这个项目里被证明是不可接受的。每一条对接结果，在放进溶剂重新审视之前，都不应该被当作最终答案。