多糖分子对接

结合案例图分析，第一张案例图是一个蓝色椭圆表示的lectin domain，里面标注了几个氨基酸残基位点（aa）。一端是Glc-Man二糖链，G代表葡萄糖，M代表甘露糖，两个单糖通过糖苷键连接。绿色虚线表示二糖链与lectin domain之间形成了氢键。从图中可以看到，甘露糖（M）残基与lectin pocket中的关键残基靠近，而葡萄糖（G）残基则伸向溶剂区。这种”单锚点+溶剂暴露”的结合模式在糖-蛋白识别中很常见。

第二张案例图横轴是结合能（kcal/mol），纵轴是密度。橙色柱状图展示了一组对接构象的结合能分布，红色曲线是核密度估计。峰值在−8 kcal/mol左右，而最佳构象用蓝色虚线标在−9.3 kcal/mol。与主峰相比，最佳构象能量低约1.3 kcal/mol，说明它属于一个显著更优的结合模式。分布整体向左偏，说明大部分构象都是放能结合，lectin对这条二糖有较好的亲和力。

从结合模式示意图和结合能分布图的组合来看，这是一个多糖分子对接项目——体系是一条Glc-Man二糖链与一个凝集素结构域（lectin domain）的相互作用，客户想解释为什么这个二糖能特异性结合，以及结合能为−9.3 kcal/mol的可靠性。第一张案例图展示了结合模式的拓扑结构，第二张案例图给出了结合能的统计置信度。

多糖对接比普通小分子对接复杂。糖环上有多个羟基，旋转自由度大，而且糖苷键的构象空间很宽。我用AutoDock Vina对接时，把二糖的所有可旋转键都设为柔性，格点盒子大小设为80×80×80 Å以容纳长链糖。受体方面，lectin domain的钙离子被保留，因为很多凝集素需要Ca²⁺稳定结合口袋几何。

计算参数上，exhaustiveness=64，num_modes=20。对接后我按能量和RMSD做了聚类，最佳cluster的RMSD<1.5 Å，说明最佳构象周围有一个稳定的低能谷。氢键分析显示Man残基的3-OH和4-OH与lectin的两个残基形成氢键，这与已知的糖识别模式一致。

我遇到的两个坑。第一个是糖环构象：吡喃糖有4C1和1C4两种椅式构象，我先用RDKit生成两种构象分别对接，发现4C1构象的结合能低约1.5 kcal/mol，符合生物活性构象。第二个是溶剂效应：糖分子高度亲水，如果不考虑去溶剂化能，打分容易过强；我加了GBSA近似修正后，结合能从−10.5调到−9.3 kcal/mol，更接近实验IC50换算值。

这张图对客户的核心价值是：解释了Glc-Man二糖与lectin domain的特异性氢键网络，并给出−9.3 kcal/mol的定量结合能，支持后续的糖阵列实验设计。关于多糖分子对接的更多细节，我整理在多糖分子对接栏目里。