[论文解读] A Novel Approach to Structure Alignment

Mattias Ohlsson, Carsten Peterson+2 more|ArXiv.org|Jun 19, 2000

Protein Structure and Dynamics参考文献 15被引用 5

一句话总结

该论文提出了一种新颖的结构对齐方法，通过迭代的平均场近似和精确的刚体变换，同时优化二值分配变量和实值原子坐标。该方法能够全面探索对齐空间，包括任意排列，实现了对距离度量和高效率的低敏感性，且结果具有概率可解释性，优于现有方法，在处理环状排列蛋白的挑战性案例中表现更优。

ABSTRACT

A novel approach for structure alignment is presented, where the key ingredients are: (1) An error function formulation of the problem simultaneously in terms of binary (Potts) assignment variables and real-valued atomic coordinates. (2) Minimization of the error function by an iterative method, where in each iteration a mean field method is employed for the assignment variables and exact rotation/translation of atomic coordinates is performed, weighted with the corresponding assignment variables. The approach allows for extensive search of all possible alignments, including those involving arbitrary permutations. The algorithm is implemented using a C_alpha representation of the backbone and explored on different protein structure categories using the Protein Data Bank (PDB) and is successfully compared with other algorithms. The approach performs very well with modest CPU consumption and is robust with respect to choice of parameters. It is extremely generic and flexible and can handle additional user-prescribed constraints easily. Furthermore, it allows for a probabilistic interpretation of the results.

研究动机与目标

• 开发一种全面探索所有可能对齐方式的结构对齐方法，包括二级结构任意排列的情况。
• 克服现有方法无法处理排列或依赖于严格假设的局限性。
• 在不使用昂贵的蒙特卡洛模拟的情况下，提供对齐结果的概率解释。
• 构建一个灵活的框架，可轻松扩展至侧链表示和用户自定义约束。
• 在多样化的蛋白质结构类别上评估性能，包括环状排列折叠等挑战性案例。

提出的方法

• 通过在二值Potts分配变量和实值原子坐标上的联合误差函数来表述结构对齐问题。
• 通过在分配变量的平均场近似和按分配加权的坐标刚体旋转/平移之间交替迭代最小化误差函数。
• 采用确定性平均场方法近似分配的后验分布，从而在无需随机采样的情况下实现概率解释。
• 通过在优化过程中隐式探索所有可能分配的完整空间，实现对二级结构的任意排列。
• 采用$C_{\alpha}$-表示法对蛋白质骨架进行简化，以降低复杂度同时保持结构保真度。
• 通过灵活的优化框架支持集成任意约束，包括自定义的间隙惩罚和函数形式。

实验结果

研究问题

• RQ1 一种结构对齐方法能否高效探索所有可能对齐方式的完整空间，包括二级结构的任意排列？
• RQ2 基于平均场的方法在准确性和计算成本方面与传统的迭代动态规划方法相比如何？
• RQ3 该方法在标准方法失效的环状排列蛋白中，能在多大程度上检测到结构相似性？
• RQ4 对齐得分的概率解释是否能提高对距离度量选择和异常值的鲁棒性？
• RQ5 该方法能否在不牺牲效率的前提下扩展至包含侧链取向和用户自定义约束？

主要发现

• 该方法成功对齐了1LEA和1XGNb4，即一对环状排列的翼状螺旋结构，识别出60个对齐残基，显著多于Vast的39个，证明其能够检测非传统结构相似性。
• 它对齐了Vast识别出的所有39个残基，同时还成功将1LEA的C端β片与1XGNb4的N端β片对齐，证明其具备检测排列折叠的能力。
• 该算法在平均速度上与Dali方法相当，且快于[5]，同时在排列案例中的准确性显著更高。
• 该方法在多样化的蛋白质家族中表现出稳健性能，对参数选择的敏感性极低，除非引入侧链。
• 该方法在无需蒙特卡洛采样的情况下提供对齐的概率解释，增强了可解释性和稳定性。
• 在涉及排列的案例中，该方法优于迭代动态规划方法，后者因结构刚性约束而无法处理此类情况。