同源建模和分子对接：从模板选择到对接验证的GPCR靶标完整流程

同源建模和分子对接的结合——先用序列相似性构建三维蛋白模型，再用分子对接预测配体结合模式——在结构生物学里是处理无晶体结构靶标的标配方案。但这套流程中，同源建模阶段的每个近似都会向下传导，最终在对接结果里被放大。理解这个传导链，比盲目跑通流程重要得多。

模板选择的序列一致性门槛

同源建模的核心假设是”序列相似=结构相似”。按Chothia和Lesk的经典结论，序列一致性>30%意味着骨架RMSD通常<1.5 Å，这是可靠的建模区。25-30%之间是”黄昏区”——模型可用但局部结构需要额外验证。低于25%的序列相似度，模型只能用来做粗粒度的域间排列分析，不能做分子对接。

团队处理过一个GPCR靶标——食欲素受体OX2R——无晶体结构可用。在PDB中搜索同源模板，找到人源β₂-肾上腺素受体（PDB 2RH1），序列一致性约31%，恰好跨过30%门槛。这个模板的选择对后续建模质量有决定性影响：β₂AR属于A类GPCR的典型结构，七根跨膜螺旋的排列高度保守，OX2R和β₂AR的跨膜区序列几乎无差距——关键差异出现在胞外loop 2（ECL2）区域，序列一致性只有约15%。

ECL2的差异至关重要，因为许多GPCR配体的结合口袋恰好位于ECL2和跨膜螺旋形成的裂隙中。模板在ECL2区域的序列差异意味着该段的三级结构不可靠——这是一个必须在对接阶段通过柔性对接或MD精修来补偿的信号。

Modeller的多模板策略与评分陷阱

同源建模工具Modeller提供了一种对抗模板偏差的手段：多模板建模。让Modeller同时基于β₂AR（2RH1）和腺苷A₂A受体（3EML，序列一致性约28%）两个模板生成约束，取两个模板在结构保守区的共识构象，降低单个模板在局部区域的偏差。

但多模板建模引入了另一个问题：DOPE评分和GA341评分的模板偏倚。两个评分都是基于统计势能，训练集本身偏向某些蛋白家族——如果两个模板都来自结构解析较多的A类GPCR，评分可能对某些局部构象存在系统性偏袒。团队的做法是不依赖单一评分，而是对Modeller生成的10个候选模型做Procheck拉氏图分析、Verify3D残基3D-1D兼容性检测，综合三项指标选出最优模型。最终选出的OX2R模型拉氏图完全允许区占93.5%（目标>90%），Verify3D评分80%以上残基通过（目标>80%）。

建模到对接的几何传导

同源建模和分子对接之间的连接点是对接口袋的几何定义。如果模型的口袋残基侧链取向与真实结构偏差较大，对接结果必然不可靠。团队在OX2R模型上，将对接口袋定义为β₂AR已知配体结合位点的同源残基对应区域——口袋中心坐标直接从模板映射，半径设定为15 Å以包含所有可能的结合构象。

对接前必须对口袋侧链做能量最小化。OX2R口袋中的Asp108、Tyr317和His350是关键残基，模型中Asp108的侧链χ₁角度偏离了β₂AR模板中同类残基的最优构象（~60°→~40°）。不做最小化直接对接，带酸性基团的配体就和Asp108形成了扭曲的盐桥——结合能被打分函数偏爱但构象物理上不合理。做了500步的侧链弛豫后，这个偏差被校正，后续对接pose的盐桥几何落在标准的2.8-3.2 Å距离范围内。

MD验证建模质量

同源模型的最终质量判断工具是MD。建模完成→口袋最优化→分子对接→对接pose MD验证，这条链条给出了最终答案。OX2R-配体复合物跑50 ns显式溶剂MD，蛋白骨架的Cα-RMSD在10 ns后稳定在~2.8 Å，跨膜螺旋保持完整，ECL2的波动较大（RMSF~3-4 Å）但在预期范围内。配体在口袋里的RMSD稳定在~2.0 Å——说明对接pose对应的构象在皮秒到纳秒尺度上是物理稳定的。

同源建模和分子对接之间的误差传导链，归根结底是”模型质量→口袋几何→对接构象→结合模式”的三级传导。模板选择是决定性的第一级，口袋弛豫补偿局部偏差，MD验证是整个链条守门员。更多蛋白建模和对接的细节讨论，站上有相关案例。

更多内容请访问 https://www.keyanxueshu.com/