[論文レビュー] Polynomial cases of the Discretizable Molecular Distance Geometry Problem

Leo Liberti, Carlile Lavor+2 more|arXiv (Cornell University)|Mar 7, 2011

Data Management and Algorithms被引用数 5

ひとこと要約

本稿は、タンパク質骨格の構造的制約を分析することで、離散化可能な分子距離幾何学問題（DMDGP）の多項式時間で解けるケースを同定する。枝刈りの際、特定の対角線または対角線以下のパターンに従って辺を刈り込む場合、探索木の幅が有界となり、最悪計算量が線形または二次時間に低下することを証明している。これは、実際のタンパク質においてBPアルゴリズムが実用的に効率的である理由を説明している。

ABSTRACT

An important application of distance geometry to biochemistry studies the embeddings of the vertices of a weighted graph in the three-dimensional Euclidean space such that the edge weights are equal to the Euclidean distances between corresponding point pairs. When the graph represents the backbone of a protein, one can exploit the natural vertex order to show that the search space for feasible embeddings is discrete. The corresponding decision problem can be solved using a binary tree based search procedure which is exponential in the worst case. We discuss assumptions that bound the search tree width to a polynomial size.

研究の動機と目的

• 一般にNP困難であるDMDGPが多項式時間で解ける条件を同定すること。
• 実際のタンパク質構造においてBranch-and-Prune（BP）アルゴリズムが実用的に効率的である理由を説明すること。
• BP探索木の幅を多項式サイズに有界にするグラフ構造（特に枝刈り辺のパターン）を特徴づけること。
• DMDGPにおける妥当な埋め込みの数が常に2の累乗であることを形式化し、証明すること。
• DMDGPおよびその変種（K-DMDGPとKを入力とするK-DMDGP）の計算量の上限を確立すること。

提案手法

• 各頂点をそのK個の直近の先行頂点からの距離に基づいて、最大2つの位置に再帰的に配置するBranch-and-Prune（BP）アルゴリズムを用いる。
• 非直近の先行頂点からの距離を検証することで、不適切な分岐を削除するための枝刈り辺を適用する。
• 探索プロセスを有向無閉路グラフ（DAG）としてモデル化し、ノードを探索木のレベル、エッジを枝刈り制約として表現する。
• 2つの枝刈り辺集合のクラスを定義・分析する：対角線パス（各頂点が一意の遠隔先行頂点を持つもの）と対角線以下のパス（複数の頂点が1つの先行頂点によって先行されるもの）。
• 頂点v > v₀が頂点u < v−Kに枝刈り辺を持つ場合、探索木の幅は2^{v₀−K}に有界であることを証明する。
• 16個のPDBタンパク質インスタンスに対して理論を実証的に検証し、12個が対角線条件を満たし、4個が対角線以下条件を満たすことが判明。両者ともv₀ = 4である。

実験結果

リサーチクエスチョン

• RQ1 グラフのどのような構造的条件下で、DMDGPのためのBranch-and-Pruneアルゴリズムが多項式時間で実行可能になるか？
• RQ2 一般にNP困難であるDMDGPに対して、実際のタンパク質骨格においてBPアルゴリズムが実用的に効率的であるのはなぜか？
• RQ3 枝刈り辺集合のどのようなグラフ理論的性質が、BPアルゴリズムにおける探索木幅の有界性を保証するか？
• RQ4 DMDGPにおける妥当な埋め込みの数が常に2の累乗であることを証明できるか、その仮定は何か？
• RQ5 K-DMDGPおよびKを入力とするK-DMDGPの計算量は何か？また、それはKにどのように依存するか？

主な発見

• K > 2の場合、DMDGPおよびその変種DMDGP_KとK-DMDGPはいずれもNP困難であることが示され、計算の困難性が確立された。
• 任意のDMDGPインスタンスにおける妥当な埋め込みの数は、常に2の累乗である。これはBPアルゴリズムの二分岐構造による証明である。
• すべての頂点v > v₀がu < v−Kを満たす頂点uに枝刈り辺を持つ場合、BP探索木の幅は2^{v₀−K}に有界となり、多項式時間計算量が保証される。
• 枝刈り辺が対角線または対角線以下のパターンに従うインスタンスでは、枝刈りの計算量が定数であれば、BPアルゴリズムは線形時間O(n)で実行される。
• 16個のPDBタンパク質インスタンスの実証的分析から、12個が対角線条件を満たし、4個が対角線以下条件を満たすことが判明。両者ともv₀ = 4であり、理論的上限を支持する。
• 枝刈り辺集合が対角線以下のパスを形成するか、対数的例外を含む場合、BPはO(n²)時間で実行され、依然として多項式時間である。