结合案例图分析，我把这个项目的完整过程、参数设置、结果解读和踩坑经验整理如下。
客户研究的是寡糖/多糖与蛋白的识别，他想知道糖链长度（聚合度 DP）从 1 到 8 时，结合亲和力如何变化，以及是否存在”结合饱和”现象。这对设计糖基药物很重要。
我用的软件是 AutoDock Vina 做对接，受体是凝集素类蛋白，配体是一系列从单糖到八糖的寡糖。由于糖分子柔性强、羟基多，我额外用 Glycam 力场做了预优化，再把低能构象送进对接。
参数设置：对接盒子 30×30×30 Å，覆盖已知糖结合位点；exhaustiveness=32，每个 DP 做 5 次重复对接取平均；打分后按 DP 做统计。
第一张案例图是 docking score 随聚合度 DP 的变化。DP 从 1 到 8，score 从约 -6.0 kcal·mol⁻¹ 逐步下降到 -8.0 kcal·mol⁻¹ 左右，但在 DP=6–8 之间基本平台化，标注为”binding saturation”（结合饱和）。这说明：增加糖单元能提升结合，但当糖链长度超过 6 个单元后，新增的糖单元与蛋白表面接触有限，自由能收益被构象熵损失抵消。
第二张案例图是氢键数量和配体效率（ligand efficiency）随 DP 的变化。红色曲线是分子间氢键数，蓝色虚线是配体效率（每重原子的结合能）。氢键数随 DP 单调增加，但配体效率在 DP=4–5 后明显下降，说明长链糖的”单位原子贡献”在降低。这个指标对药物设计更关键——我们不仅要结合强，还要单位质量效率高。
第一，糖的羟基质子化状态必须正确。我最初用默认 AutoDock 原子类型，部分羟基氢取向不对，导致氢键计数偏低；用 Glycam 重新分配后结果才合理。
第二，DP>6 后柔性问题严重。对接打分假设配体是刚性或半刚性，但八糖有几十个可旋转键，对接可能找不到全局最优。我建议客户对 DP≥6 的体系加 MD 抽样或 FEP 计算。
第三，配体效率比绝对 score 更有指导意义。客户一开始只关注 score 越来越负，我提醒他：药物分子通常希望单位重原子贡献高，长链糖虽然 score 好，但 MW 也高。
下次做多糖对接，我会把 DP 扫描和配体效率分析作为标准输出。对 DP 较小的体系用对接快速筛选，对 DP 较大的体系建议加 FEP 或 metadynamics。这样客户既能拿到趋势，也能知道哪个长度是”甜蜜点”。
计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验 需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。
