[論文レビュー] RNA-LEGO: Combinatorial Design of Pseudoknot RNA

Emma Yu Jin, Christian M. Reidys|ArXiv.org|Nov 9, 2007

RNA and protein synthesis mechanisms参考文献 16被引用数 13

ひとこと要約

本稿では、生成関数と転送定理を用いて、最小スタック長σを満たすk非交差RNA偽二次構造を列挙する組合せ的フレームワークRNA-LEGOを紹介する。精密な漸近公式を導出し、スタック長≥2の3および4非交差構造がそれぞれ~311.2470 × 4! / (n(n−1)⋯(n−4)) × 2.5881^nおよび~1.217×10^7 × n^(-21/2) × 3.0382^nの割合で増加することを示しており、これは、塩基対を分離する制約がかかると指数的増加が著しく低下することを示しており、RNA偽二次構造における中立的ネットワークの存在を支持する。

ABSTRACT

In this paper we enumerate $k$-noncrossing RNA pseudoknot structures with given minimum stack-length. We show that the numbers of $k$-noncrossing structures without isolated base pairs are significantly smaller than the number of all $k$-noncrossing structures. In particular we prove that the number of 3- and 4-noncrossing RNA structures with stack-length $\ge 2$ is for large $n$ given by $311.2470 \frac{4!}{n(n-1)...(n-4)}2.5881^n$ and $1.217\cdot 10^{7} n^{-{21/2}} 3.0382^n$, respectively. We furthermore show that for $k$-noncrossing RNA structures the drop in exponential growth rates between the number of all structures and the number of all structures with stack-size $\ge 2$ increases significantly. Our results are of importance for prediction algorithms for pseudoknot-RNA and provide evidence that there exist neutral networks of RNA pseudoknot structures.

研究の動機と目的

• k非交差RNA偽二次構造の最小スタック長制約を満たす組合せ的フレームワークの構築。
• 特にk ≥ 3の場合に、塩基対の分離が禁止された場合の構造的複雑さの低下を定量化すること。
• スタック長≥2の3および4非交差RNA構造の数に対する精密な漸近公式の導出。
• スタック長およびアーキン長制約に依存するk非交差構造の指数的増加率の調査。
• 構造的密度と増加率の分析を通じて、RNA偽二次構造における中立的ネットワークの存在に対する理論的証拠の提供。

提案手法

• 著者らは、最小スタック長σを満たすk非交差RNA構造の漸近的挙動を分析するため、生成関数と転送定理を用いる。
• タイプ(k,1)の構造の2変数生成関数を導出し、u₀ = (x²)^{σ−1} / ((x²)^σ − x² + 1) を含む変換により、より高いスタック長に一般化する。
• T_{k,σ}(n) の生成関数を、特定のバシレートテーブルの数を数えるf_k(2n,0) の生成関数で表現する。
• 重要なステップとして、Pringsheimの定理を適用し、T_{k,σ}(n) の生成関数の支配的特異点γ_{k,σ}を正の実数として特定し、指数的増加率を確認する。
• 転送定理を用いて、f_k(2n,0) の部分指数的要因をT_{k,σ}(n) に移すことで、精密な漸近展開を得る。
• このフレームワークは、k非交差構造の列挙を生成関数の特異展開と結びつけ、複素解析的手法による漸近的解析を可能にする。

実験結果

リサーチクエスチョン

• RQ1 スタック長≥2の3非交差RNA構造の漸近的数は何か？
• RQ2 塩基対の分離が禁止された場合、k非交差RNA構造の指数的増加率はどのように変化するか？
• RQ3 最小スタック長σ ≥ 2の4非交差RNA構造の漸近的挙動は何か？
• RQ4 RNA偽二次構造の組合せ的構造は、中立的ネットワークの存在を支持できるか？
• RQ5 最小アーキン長およびスタック長制約は、k非交差RNA構造の総数にどのように影響するか？

主な発見

• スタック長≥2の3非交差RNA構造の数は、十分に大きなnに対して漸近的に ~311.2470 × 4! / (n(n−1)(n−2)(n−3)(n−4)) × 2.5881^n に比例する。
• スタック長≥2の4非交差RNA構造の数は、十分に大きなnに対して漸近的に ~1.217×10^7 × n^(-21/2) × 3.0382^n に比例する。
• スタック長≥2が課された場合、k非交差RNA構造の指数的増加率は著しく低下し、構造的複雑さが顕著に減少することが示された。
• T_{k,σ}(n) の生成関数の支配的特異点γ_{k,σ}が正の実数であることが示され、指数的増加率が正当化された。
• 漸近的解析により、k=3の場合は増加率の部分指数的要因がn^{-1/2}、k=4の場合はn^{-21/2}であることが判明し、構造的制約が強化されていることが反映されている。
• これらの結果は、高い組合せ的複雑さにもかかわらず、利用可能な構造数が管理可能であるため、RNA偽二次構造における中立的ネットワークの存在に対する強い証拠を提供する。