[論文レビュー] Asymptotic Evolution of Protein-Protein Interaction Networks for General Duplication-Divergence Models

Kirill Evlampiev, Hervé Isambert|ArXiv.org|Nov 22, 2006

Bioinformatics and Genomic Networks参考文献 1被引用数 3

ひとこと要約

本稿では、6つのγパラメータを用いて突然変異由来相互作用の確率的選択を組み込んだ一般化された二重化・変異（GDD）モデルを提案し、タンパク質間相互作用（PPI）ネットワークの漸近的進化を研究する。解析的に長期間におけるネットワークトポロジーを導出し、生物学的に妥当な条件下では、保存的かつ非密なネットワーク—必然的にスケールフリーである—が出現することを示している。

ABSTRACT

Genomic duplication-divergence events, which are the primary source of new protein functions, occur stochastically at a wide range of genomic scales, from single gene to whole genome duplications. Clearly, this fundamental evolutionary process must have largely conditioned the emerging structure of protein-protein interaction (PPI) networks, that control many cellular activities. We propose and asymptotically solve a general duplication-divergence model of PPI network evolution based on the statistical selection of duplication-derived interactions. We also introduce a conservation index, that formally defines the statistical evolutionary conservation of PPI networks. Distinct conditions on microscopic parameters are then shown to control global conservation and topology of emerging PPI networks. In particular, conserved, non-dense networks, which are the only ones of potential biological relevance, are also shown to be necessary scale-free.

研究の動機と目的

• 遺伝子の二重化と変異を通じて、タンパク質間相互作用（PPI）ネットワークの長期的進化をモデル化する一般理論的枠組みの構築。
• PPIネットワークが進化的に保存的であり、トポロジカルに非密であること—生物学的意味を持つために必要な—条件の特定。
• 二重化由来相互作用の統計的選択に基づいて、PPIネットワークの保存性インデックスを形式的に定義・分析すること。
• 指数的二重化ダイナミクス下でのPPIネットワークの漸近的挙動を解明し、従来の時間線形モデルと対比すること。
• ドメインシャッフルや水平的遺伝子移動などの追加的進化過程を組み込むための理論的ベースラインを確立すること。

提案手法

• 各時刻ステップで一定割合 $ q $ の遺伝子が二重化され、その後非対称的変異によって「古（o）」と「新（n）」の二重体に分かれる一般化二重化・変異（GDD）モデルを提案。
• 二重化状態（s: 営利、o: 古、n: 新）に基づく遺伝子ペア間の相互作用の確率的保存をモデル化するための6つのγパラメータ（$ \gamma_{ij} $）を導入。
• ノードタイプと接続性を追跡する生成関数 $ H^{(n)}(x,y) $、$ T^{(n)}(x,y,z) $、$ F^{(n)}(x) $ を用いて、ネットワーク特徴（例：次数分布、クラスタリング、三角形）の再帰的関係を導出。
• すべての進化実現のアンサンブル平均 $ \langle Q^{(n)} \rangle $ を適用し、漸近的ネットワーク特性を研究。
• 全7パラメータモデルを、局所的（$ q \ll 1 $）および全ゲノム的（$ q = 1 $）二重化の3パラメータ極限に簡略化。
• 再帰的関係を漸近的に解き、PPIネットワークの長期的成長およびトポロジーを特定。特に、三角形数などの特徴の指数的増加を明らかにした。

実験結果

リサーチクエスチョン

• RQ1 γパラメータにどのような条件下で、PPIネットワークが生物学的妥当性を満たすように保存的かつ非密なままであるか？
• RQ2 PPIネットワークの漸近的トポロジー、特にスケールフリーデータ分布が、二重化・変異ダイナミクスからどのように生じるか？
• RQ3 指数的二重化ダイナミクスが、時間線形モデルと比較してPPIネットワーク進化にどのように寄与するか？
• RQ4 局所的二重化と全ゲノム二重化の異なる二重化レジームが、PPIネットワークの保存性および構造的特性に与える影響は何か？
• RQ5 二重化由来相互作用の統計的選択に基づいて、PPIネットワークの形式的保存性インデックスを定義できるか？

主な発見

• 保存的かつ非密なPPIネットワークは、γパラメータが特定の制約を満たす場合にのみ可能であり、これは二重化由来相互作用のうち一部のみが保持されることを保証する。
• 本モデルは、保存的かつ非密なネットワークが必然的にスケールフリーであると予測しており、指数的漸近的発展においてパワー則次数分布が出現することを示している。
• ネットワーク内の三角形の平均数は、時間とともに指数関数的に増加し、その増加率はγパラメータと二重化割合 $ q $ の組み合わせによって決定される。式（84）に示されている。
• Ispolatov *et al.* の時間線形二重化・変異モデルは、$ q \ll 1 $ かつ特定のγパラメータが設定されたGDDモデルの特殊極限として回復される。
• 全ゲノム二重化（$ q = 1 $）では、ネットワーク進化は3つのγパラメータ $ \gamma_{oo} $、$ \gamma_{on} $、$ \gamma_{nn} $ によって支配され、それぞれ古・古、古・新、新・新の二重体間の相互作用保存を制御する。
• ネットワーク特徴の再帰的関係は対称性を保ち、各ノードタイプ（s, o, n）が新しい接続およびモチーフに与える寄与を追跡することで導出され、正確な漸近的解析が可能になる。