[論文レビュー] 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] よりも高速であり、特に並び替えを伴うケースでは著しく高い精度を発揮している。
• パラメータ選択にほとんど感受性がなく、側鎖を含めない限り、多様なタンパク質ファミリーにおいて安定した性能を発揮している。
• モンテカルロサンプリングを必要とせず、整合の確率的解釈を提供することで、解釈可能性と安定性が向上している。
• 構造的剛性制約により並び替えを扱えない従来の反復的動的計画法と比較して、本手法は並び替えを伴うケースで優れた性能を発揮している。