IF36.1｜全自动结构比对套件：从链映射到配体评分的完整方案

OpenStructure新套件实现大复合物全自动、可扩展的结构比对与评分，覆盖链映射、界面、配体。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

蛋白质结构预测领域长期依赖CASP、CAMEO、CAPRI等基准测试来评估模型准确性。随着AI方法如AlphaFold、RoseTTAFold的兴起，对大规模数据集和无需人工干预的稳健评分需求日益迫切。现有评分方法在处理大分子复合物时面临挑战：链映射问题随链数呈阶乘增长，传统枚举法不切实际；基于全局叠加的评分对结构域运动敏感，而局部叠加无关评分如LDDT、CAD-score则更具鲁棒性。此外，小分子配体姿态预测的评估需要自动化的配体识别、对称性校正和结合位点叠加。OpenStructure框架旨在整合这些需求，提供一致的自动化基准测试工具。

科学问题

卡在哪里： 现有链映射方法多依赖全局叠加，不适用于基于接触的评分；针对大复合物的局部叠加无关评分缺乏有效的链映射算法。配体评估方面，已有工具如Spyrmsd仅支持单一配体对，无法处理多配体复合物；基于指纹的评分依赖主观相互作用集和手动准备，难以复现。此外，缺乏一个统一框架同时覆盖蛋白质、核酸、小分子配体及复合物化学计量学的自动评分。

本文要回答： 本文旨在开发一套全自动、快速、可靠的评分工具，用于比较理论模型与实验参考结构，涵盖链映射、三级/四级结构、蛋白-蛋白界面及蛋白-配体复合物，并提供使用建议。

技术路线

作者主要用了：QSMap/QSMapR链映射、LDDT扩展、BiSyRMSD、LDDT-PLI、LDDT-LP、OpenStructure框架。

作者首先开发了QSMap和QSMapR两种链映射算法，分别优化基于接触的QS-score/LDDT和基于叠加的r.m.s.d.，以解决大复合物的自动映射问题。随后扩展LDDT实现以支持多链和核苷酸，并引入i-LDDT和bb-LDDT变体。针对小分子配体，开发了BiSyRMSD、LDDT-PLI和LDDT-LP三种评分，自动处理配体识别、对称性校正和结合位点叠加。所有方法集成于OpenStructure框架，提供命令行和Web接口，并整合了DockQ、ICS/IPS、GDT、TM-score等外部评分，形成一站式基准测试套件。

图 1　Fig. 1。图1展示了一个异源二聚体参考结构（含小分子配体）与模型比较的示意图。评估包括映射等效组分并评分差异，模型的不准确性可能涉及化学计量错误、三级结构差异、界面差异、结合位点大小形状错误以及配体构象或姿态错误。与“找不同”游戏不同，差异必须定量评估。

核心亮点

亮点 1｜QSMap/QSMapR高效处理大复合物链映射

QSMap优化QS-score或LDDT，QSMapR优化r.m.s.d.，均采用启发式贪心算法，从所有可能的初始链对开始扩展，避免阶乘枚举。在包含多达30条链的测试中，运行时间在100秒左右，而朴素枚举在链数超过10时已不可行。在CASP15装配模型数据集上，QSMap在接触类评分（QS-score、ICS、LDDT）上优于QSMapR，而QSMapR在全局拓扑评分（TM-score）上表现更好。与US-align、Foldseek-Multimer和AlphaFold-Multimer的链映射算法相比，QSMapR的映射质量与US-align相当，但速度快约一个数量级。

图 2　Fig. 2。图2a显示QSMap/QSMapR与朴素枚举的运行时间对比，朴素枚举在N>10时不可行，而QSMap/QSMapR可处理更大结构。图2b显示QSMap在QS-score上优于QSMapR，图2c显示QSMapR在全局叠加指标（如TM-score）上更优。阴影区域突出绝对差异≥0.1的数据点。

亮点 2｜LDDT扩展支持多链和核苷酸

新版LDDT实现与QSMap链映射紧密集成，能够处理大型装配体，并扩展至核苷酸。新增i-LDDT（仅考虑界面间距离）和bb-LDDT（仅考虑主链原子）变体。在CASP15 RNA建模挑战中进行了广泛测试，核苷酸的立体化学检查在CASP16中实现。该实现解决了CASP13中链映射算法处理大装配体不足的问题，为多链蛋白和RNA复合物提供了稳健的局部评分。

图 3　Fig. 3。图3a是LDDT-PLI与BiSyRMSD的散点图（n=31,905），两者强负相关，但存在离群值。图3b-d展示三个离群例子：侧链翻转导致LDDT-PLI略低；配体部分远离导致BiSyRMSD巨大而LDDT-PLI非零；镁离子微小位移对LDDT-PLI几乎无影响。

亮点 3｜三种配体评分互补评估姿态预测

BiSyRMSD是结合位点叠加后的对称校正r.m.s.d.，衡量配体绝对位置精度；LDDT-PLI是限制在聚合物-配体接触的LDDT，评估天然接触的再现并惩罚额外接触；LDDT-LP评估结合位点残基的准确性。在CASP15配体建模挑战的31,905个配体对上，LDDT-PLI与BiSyRMSD强负相关（Spearman ρ=-0.989），但存在离群值：如侧链翻转导致LDDT-PLI略低而BiSyRMSD极低；配体部分远离导致BiSyRMSD巨大而LDDT-PLI非零；镁离子微小位移对LDDT-PLI几乎无影响。额外接触惩罚能有效检测结合位点内的多余链或错误侧链。

图 4　Fig. 4。图4展示了LDDT-PLI中额外接触惩罚的效果。图4a中一条额外链穿过结合位点，导致LDDT-PLI降至0.53，而BiSyRMSD和LDDT-LP仍很高。图4b中末端环和Arg6侧链错位靠近配体，LDDT-PLI降至0.65，BiSyRMSD和LDDT-LP受影响较小。

亮点 4｜OpenStructure提供一站式自动评分套件

OpenStructure框架整合了上述所有评分，并纳入DockQ、ICS/IPS、GDT、TM-score、CAD-score等外部工具，提供命令行动作compare-structures和compare-ligand-structures，以及SWISS-MODEL Structure Assessment服务器和REST API。评分建议包括：局部叠加无关评分（LDDT、CAD-score）优于全局叠加评分以处理结构域运动；r.m.s.d.类评分对异常值敏感，应谨慎使用；接触类评分自然惩罚不完整模型，但r.m.s.d.类不惩罚；不完整参考结构下应使用QS-best或trimmed ICS/IPS变体；DockQ-ave比DockQ-wave更适合小界面重要的场景。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： QSMap/QSMapR链映射、LDDT扩展、BiSyRMSD、LDDT-PLI、LDDT-LP、OpenStructure框架
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事蛋白质结构预测、复合物建模或配体对接评估的生信研究人员和工具开发者。

带走一句

做复合物结构预测评估时，应优先采用局部叠加无关评分（LDDT、QS-score）并配合QSMap链映射；配体评估需同时使用BiSyRMSD和LDDT-PLI以捕捉不同错误类型。