分子动力学模拟自由能计算

结合案例图分析，本文以这个构象转变自由能项目为例，讲 umbrella sampling 的窗口与重叠。这是一个构象转变自由能测定项目。客户要研究某体系的构象转变或解离过程的动力学和热力学壁垒。第一张案例图是 umbrella sampling 得到的 PMF 曲线，翻越约 7.20 kcal/mol 能垒到 state B。第二张案例图是多个采样窗口分布，相邻窗口有重叠。

项目过程：构象转变自由能项目的执行链路

方法上自由能计算有热力学积分（TI）、自由能微扰（FEP）、umbrella sampling（US）几条路。客户要的是沿一条明确反应坐标的 PMF 剖面，US 最直观且能直接读出能垒，所以我选了 US。TI 更适合连续参数变化的自由能差，在本例的离散构象转变坐标上不如 US 直观；窗口数我定为 20 个，是收敛测试的结果。 我的方法：确定反应坐标，沿坐标均匀布 20 窗口，每窗口用谐振子偏置（k=1000 kJ/mol·nm²）限制坐标在某区间分别跑 10 ns MD，最后 WHAM 重构 PMF。客户要的是能垒数值和转变路径，这两张图正好给出。

案例分析：从 PMF 到活化能

图里读出的数要落地。第一步确认窗口充分重叠——相邻窗口分布重叠 20–40% 是 WHAM 收敛必要条件，重叠不足 PMF 会不连续。我对每个窗口算采样量，确保每个窗口 ≥ 5 ns 有效采样、重叠区统计误差 < 0.2 kcal/mol，否则 PMF 在接缝处抖动。 第二步用 WHAM 拼接，得到 ΔG‡=7.20±0.15 kcal/mol，并标注这是包含溶剂化的表观自由能。误差量级上，窗口数从 10 增到 20 时能垒从 6.1 升到 7.20 kcal/mol，少窗口会系统性低估约 1 kcal/mol。

总结教训：窗口重叠与偏置强度的坑

做的过程中我踩过实打实的坑。第一坑是窗口数：我画过 PMF 随窗口数的收敛曲线——10 个窗口时能垒只有 6.1 kcal/mol 且曲线锯齿状，15 个到 6.9，20 个稳定在 7.20±0.15 kcal/mol，证明窗口不够会系统性低估能垒。 第二坑是偏置强度 k，太弱窗口分布太宽（相邻重叠 > 60% 浪费采样），太强采样受限难越势垒，我试跑调 k 使相邻重叠维持在 20–40%。第三坑是反应坐标选错，初版用质心距离结果漏掉了一个扭转中间态，改成扭转角后才看到完整双势阱。

项目经验：自由能交付的可复用标准

作为可复用的项目经验，第一，反应坐标选择决定可靠性，选错会漏中间态。第二，WHAM 和 MBAR 各有优劣，WHAM 稳但假设窗口独立，MBAR 精但计算量大，按体系大小选；本例 MBAR 与 WHAM 偏差仅 0.1 kcal/mol。第三，误差带 ±0.15 kcal/mol 必须标注。客户据此估算转变速率常数。说到底，自由能项目的交付底线是 PMF 平滑可拼接；窗口不重叠，再漂亮的 7.20 kcal/mol 也只是碎在一地的曲线。

分子动力学模拟自由能计算

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验

需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。