核酸構造予測

核酸構造予測は、核酸の配列から二次構造および三次 構造を決定する計算手法です。二次構造は、1つまたは複数の核酸配列から予測できます。三次構造は、配列から、または相同配列の構造が既知の場合の比較モデリングによって予測できます。

核酸の二次構造予測は、主に塩基対形成塩基スタッキング相互作用に依存しています。多くの分子は複数の三次元構造を持つ可能性があるため、トランスファーRNA(tRNA)やマイクロRNA (miRNA)といった既知の核酸分子との明らかな配列および機能的類似性が認められない限り、これらの構造を予測することは困難です。多くの二次構造予測手法は動的計画法のバリエーションに依存しているため、擬似ノットを効率的に識別することができません。

RNAとDNAの構造予測手法は類似しているものの、そのアプローチには若干の違いがある。生体内では、DNA構造は2本の鎖が完全に相補的な二重鎖構造をとることが多いのに対し、RNA構造はリボソームスプライソソーム転移RNAなどの複雑な二次構造や三次構造に折り畳まれることが多い。これは、RNA中の余分な酸素が核酸骨格における水素結合の傾向を高めるためである。エネルギーパラメータも2つの核酸で異なる。構造予測手法は、完全に理論的なアプローチに従う場合もあれば、実験データを取り入れたハイブリッドなアプローチに従う場合もある。[ 1 ] [ 2 ]

単一配列構造予測

RNAを研究する研究者にとって共通の課題は、核酸配列のみに基づいて分子の三次元構造を決定することです。しかし、RNAの場合、最終的な構造の多くは分子の二次構造、つまり分子内塩基対形成によって決定されます。これは、多様な種において塩基対形成が高度に保存されていることからも明らかです。

熱力学モデル

小さなRNA分子の二次構造は、主に水素結合塩基スタッキングといった強い局所的相互作用によって決定されます。こうした相互作用の自由エネルギーを合計することで、与えられた構造の安定性を概算できます。与えられた二次構造の折り畳み自由エネルギーを予測するために、経験的な最近傍モデルが用いられます。最近傍モデルでは、各モチーフの自由エネルギー変化は、モチーフの配列とそれに最も近い塩基対の配列に依存します。(「近傍」という側面により、塩基スタッキングエネルギーを記述できます。)断片のエネルギーを合計することで、全体の自由エネルギーが得られます。最小の自由エネルギー(MFE)を持つ構造が最も安定です。[ 3 ]

RNA モデルの例には次のようなものがあります。

  • ターナー1999 [ 4 ]およびターナー2004 [ 5 ]。
  • Andronescu 2007 は、4つの標準塩基のいずれかを持つRNAのパラメータを示しています。これは、以前のウェットラボ実験と他の形式のRNA構造データから導き出されています。[ 6 ] [ 7 ]
  • Langdon 2018は、4つの標準塩基のいずれかを含むRNAのパラメータを示している。[ 8 ]

2020年現在、ほとんどのソフトウェアパッケージはTurner 2004をデフォルトとしています。2020年のベンチマークでは、RNA-RNA相互作用タスクにおいてTurner 2004は4つのアルゴリズムの中で最も性能が悪く、最も性能が高かったAndronescu 2007でも5ペア上回っています。構造予測においてはLangdon 2018が最も優れていることが知られています。[ 9 ]

最も安定した構造

MFE構造を見つける最も簡単な方法は、すべての可能な構造を生成し、その自由エネルギーを計算することですが、RNAの長さに応じて、配列の可能な構造の数は指数関数的に増加します。[ 10 ]分子が長くなると、可能な二次構造の数は膨大になります。100ヌクレオチドの配列には、10の25乗以上の二次構造が考えられますそのため、よりスマートな手法が必要になります。[ 3 ]number of secondary structures = (1,8)N, N- number of nucleotides

動的計画法アルゴリズム

RNAとDNAの二次構造を予測する最も一般的な方法は、動的計画法である。[ 11 ] [ 12 ] RNAの二次構造を予測する初期の試みの1つは、ルース・ヌシノフと同僚らによって行われた。彼らは、一連の「ブロック」(ポリヌクレオチド鎖)の長さと数を最大化する動的計画法ベースのアルゴリズムを開発した。[ 11 ]各「ブロック」には少なくとも2つのヌクレオチドが必要であったため、単一塩基マッチング手法に比べてアルゴリズムの記憶域要件が削減された。[ 11 ]ヌシノフらは後に、以前の折り畳みの結果を保存しながら徐々にサイズの大きいサブセクションを折り畳むことによってRNAのサイズ制限を約1,000塩基まで増やした、パフォーマンスが向上した適応型手法を発表した。これは現在ヌシノフアルゴリズムとして知られている。[ 12 ] 1981年にマイケル・ズーカーとパトリック・スティグラーは、ヌシノフらの解決策に匹敵する性能を持ちながら、「最適ではない」二次構造も見つけられるという追加の機能を備えた改良されたアプローチを提案した。[ 13 ]

動的計画法アルゴリズムは、RNA二次構造の可能性のあるすべてのバリエーションを、構造を明示的に生成することなく暗黙的に検証する手段を提供します。まず、最短の断片から始めて、各可能な配列断片について、最も低いコンフォメーション自由エネルギーを決定します。次に、より長い断片について検討します。より長い断片については、短い配列に対して決定された最適な自由エネルギー変化を再帰的に計算することで、最も低い折り畳み自由エネルギーの決定を高速化します。配列全体の最低自由エネルギーが計算されると、RNA分子の正確な構造が決定されます。[ 3 ]

動的計画法アルゴリズムは、 「よくネストされた」塩基対形成パターン、すなわち、配列位置において互いに重なり合わない塩基とのみ水素結合を形成するパターンを検出するために一般的に用いられます。このカテゴリーに該当する二次構造には、二重らせんステムループ、そして転移RNA分子に見られる「クローバーリーフ」パターンの変異体が含まれます。これらの手法は、ワトソン・クリック型塩基対やフーグスティーン型塩基対など、特定の種類の塩基対形成相互作用に関連する自由エネルギーを推定する、事前に計算されたパラメータに依存しています。手法の複雑さに応じて、塩基スタッキングの影響を組み込むために、単一の塩基対や、2塩基または3塩基の短いセグメントが考慮される場合があります。この手法では、計算コストが非常に高くなるアルゴリズムの大幅な変更なしには、よくネストされていない擬似ノットを同定することはできません。 [ 14 ]

擬似結び目の予測

RNAの二次構造を予測する際の問題の1つは、標準的な自由エネルギー最小化法と統計的サンプリング法では擬似ノットを発見できないことである。[ 4 ]主な問題は、通常の動的計画法アルゴリズムでは、二次構造を予測する際に、最も近いヌクレオチド間の相互作用のみが考慮されるのに対し、擬似ノット構造は離れたヌクレオチド間の相互作用によって形成されることである。RivasとEddyは、擬似ノットを予測するための動的計画法アルゴリズムを発表した。[ 14 ]しかし、この動的計画法アルゴリズムは非常に遅い。自由エネルギー最小化のための標準的な動的計画法アルゴリズムは、時間がO(N 3 )(Nは配列中のヌクレオチドの数)でスケールするのに対し、RivasとEddyのアルゴリズムは時間がO(N 6 )でスケールする。このため、複数の研究者が擬似ノットのクラスを制限したアルゴリズムのバージョンを実装し、パフォーマンスの向上を実現した。たとえば、pknotsRGツールは、単純な再帰擬似ノットのクラスのみを含み、時間がO(N4)でスケールする。[ 15 ]

擬似結び目の予測という一般的な問題はNP完全であることが示されている。[ 16 ]

最適ではない構造

S. cerevisiae tRNA Pheの構造空間。真の生物学的構造はMFEに基づく構造とは異なり、これは主に多くの修飾塩基の存在による。エネルギーと構造はRNAsuboptを用いて計算され、構造距離はViennaRNA 1.xのRNAdistanceを用いて算出された。Bos taurus tRNA Pheも同様の状況にあるが、修飾塩基を適切に補正すればMFEとして正しい構造を復元できる(ViennaRNA 2.6.0 of 2023)。[ 17 ]

自由エネルギー最小化による 1 つの配列からの RNA 二次構造予測の精度は、いくつかの要因によって制限されます。

  1. 最近傍モデルの自由エネルギー値リストは不完全または不正確である可能性があります。
    • 修飾塩基は修飾されていない塩基と同じように折り畳まれると想定されることが多いが、実際には微妙な違い(自由エネルギーが0.5 kcal/mol程度変化する)から大きな違い(修飾により水素結合の形成が物理的に阻害される)までがある。[ 17 ]
  2. 既知の RNA のすべてが、熱力学的最小値に適合するように折り畳まれるわけではありません。
  3. 一部の RNA 配列には、複数の生物学的に活性なコンフォメーション(リボスイッチなど)があります。

このため、類似した低自由エネルギー構造を予測する能力は、重要な情報を提供する可能性があります。このような構造は、準最適構造と呼ばれます。MFOLDは、準最適構造を生成するプログラムの一つです。[ 18 ]

構造の集合は、「ドットプロット」、つまり2つの塩基がペアになる頻度を含む正方行列を視覚化したものによって表すことができます。これは本質的にヒートマップの一種です。

ボルツマン集団

可能性のある構造の集合はアンサンブルを形成します。このアンサンブルからボルツマン分布(プログラムSFOLDで例示)に従ってサンプリングすることで、多くの可能性のある構造が得られます。[ 19 ] [ 20 ]

ボルツマンアンサンブルの状態は「ドットプロット」で表すこともできます。ここで、大きさはアンサンブル内で任意の2つの塩基がペアになる確率です。これらの大きさは暗黙的に熱力学データも表します。

比較二次構造予測

進化的に関連する複数のRNA配列は、部位間の共変異のパターンを共有しています。2つの部位が頻繁に同期して変化することが観察される場合、それらの部位間に構造上必要な水素結合が存在する可能性が高いと考えられます。これは、真の構造がどのようなものであるかについての新たな情報源となります。[ 16 ]

一般的に、アライメントの問題とコンセンサス構造の予測は密接に関連しています。コンセンサス構造の予測には、いくつかの異なるアプローチがあります。[ 21 ] [ 22 ]

  1. アライメントの折り畳み
  2. 特定の予測構造のア​​ライメント(場合によっては配列アライメントと併用)
  3. 予測された構造の配列と集合の同時アライメント

これらの各ステップは、反復的に組み合わせることもできます。[ 22 ]

揃えて折り畳む

実用的なヒューリスティックアプローチとしては、多重配列アライメントツールを使用して複数のRNA配列のアライメントを作成し、コンセンサス配列を見つけてフォールディングする方法があります。共変動は、相同性の高い複数のRNA配列と、関連はあるものの類似性のない配列の多重配列アライメントから抽出できます。アライメントの質がコンセンサス構造モデルの精度を決定します。コンセンサス配列は、個々の構造予測問題と同様に、様々なアプローチを用いてフォールディングされます。

  • PfoldはSCFGを実装しています。共変動は塩基対確率の「ドットプロット」に変換され、これは通常のPfoldの単一配列バージョンの上に重ねて考慮されます。[ 23 ]
  • ViennaスイートのRNAalifoldは、本質的には同スイートのRNAfoldのアライメントベースの亜種です。どちらも熱力学的フォールディングアプローチを採用しており、2つのモードで実行できます。1つは単一の構造を返すMFE(最小自由エネルギー)モード、もう1つは複数の構造候補を返すパーティショニングモードです。[ 24 ]コンセンサス構造のスコアは、平均熱力学的エネルギーと共変動から得られる追加ポイントの組み合わせです。[ 22 ]
    • 線形時間バリアントはLinearAlifoldとして発見されました。[ 25 ]
  • ILM(反復ループマッチング)は、熱力学相互情報量(共分散)のコンテンツスコアを組み合わせて使用​​します。その名前が示すように、これはヌシノフ型の予測構造間のループマッチングに基づいています。各反復において、予測された二次構造から「最適な」ヘリックス/ステムが発見され、残りの構造が再び折り畳まれます。これにより、擬似結び目構造の予測が可能になります。[ 26 ]

1つ折りにして位置合わせ/ジョイント位置合わせ

進化はしばしばRNA配列よりも機能的なRNA構造をより良く保存します。[ 24 ]そのため、2つ以上の大きく異なるが相同性のあるRNA配列から共通の構造を推定することは、生物学的によくある問題です。2つの配列の配列類似性が50%未満の場合、配列アライメントは適切ではなく、構造予測の精度向上には役立ちません。そのため、アライメントと二次構造予測を組み合わせたアプローチが必要になります。[ 22 ]

アライメントと二次構造予測を組み合わせる明白なアプローチは、単一配列構造予測法を用いて配列を折り畳み、得られた構造をアライメントすることです。2つの構造のアライメントは通常、2つの構造をツリー構造に解析し、2つのツリー間の最小編集距離を見つけることで行われます。[ 27 ]この方法は、 配列と構造のジョイントアライメントを実行する場合でもかなり高速に実行され、固定された、変化しない、潜在的に不正確なアライメントから開始するよりも優れた結果をもたらします。[ 28 ]

ドットプロットを折り曲げてジョイントを合わせる

特異構造と樹形比較の利用には依然として弱点があります。真の(生理学的)構造は必ずしも最小自由エネルギー構造ではないことはよく知られています。また、樹形編集によるミスペアリングの修復能力には限界があることも分かっています。これらの問題に対処するには、可能性のある構造の集合全体を比較するという、さらに別のアプローチが必要になります。

可能な構造のアンサンブルの表現方法としてはドットプロットが既に存在し、ここでも有用である。アンサンブル構造を考慮しながら 2 つの配列をアラインメントするには、まず各配列の塩基対確率ドットプロットを予測する。次に問題は、配列とドットプロットの結合アラインメントとなる。この問題のプロトタイプアルゴリズムはサンコフアルゴリズム[ 29 ]であり、基本的には配列アラインメントとヌシノフ (最大対合) [ 11 ]フォールディング動的計画法を融合したものである。[ 30 ]オリジナルのサンコフアルゴリズムは理論的な演習であり、膨大な計算リソースを必要とする (時間でO (n3m) 、空間で O (n2m)、ここで n は配列長、m は配列数) ためである。

Sankoffのアルゴリズムの制限版を実装した注目すべき試みとしては、Foldalign、[ 31 ] [ 32 ] Dynalign、[ 33 ] [ 34 ] PMmulti/PMcomp、[ 30 ] Stemloc[ 35 ] Murlet、[ 36 ] LocARNA-P [ 37 ]などがある。これらの実装は実行には実用的であるが、アライメントの最大長や可能なコンセンサス構造のバリエーションなどの特性に追加の制限を課している。例えば、Foldalignは局所的なアライメントに焦点を当てており、配列アライメントの可能な長さを制限している。

CARNAでは、擬似ノットや複数の安定構造など、塩基対のより多様な配置を可能にするために、ジョイントアライメントに対する異なるアプローチが用いられています。これはMAX-SNP困難問題ですが、制約プログラミングとハードリミットにより、実行時間は実用的なレベルに保たれています。[ 38 ]

どのジョイントアライメント手法を用いる場合でも、一度に処理できるのは2つの配列(または2つの事前アライメント済み配列グループ)のみです。これは、ガイドツリーを用いた 通常の反復アライメントによって解決されます。

三次構造予測

RNAの二次構造が判明したら、次の課題は三次構造の予測です。最大の課題は、二本鎖らせん構造間の構造を決定することです。また、RNA分子はしばしば転写後修飾を受けたヌクレオシドを含んでおり、これが新たな非典型的な相互作用を引き起こす可能性があるため、三次構造予測には多くの困難をもたらします。[ 39 ] [ 40 ] [ 41 ] [ 42 ]

比較アプローチでは、関連する既知の構造をデータベース(例えば、3D構造のPDB)から抽出します。この構造をテンプレートとして用い、問題の配列をその上に移植し、調整します。[ 43 ]

de novoアプローチまたはab initioアプローチでは、既知の3D構造は使用されません。最初のステップは、上記のいずれかのアプローチを用いて二次構造を予測することです。[ 44 ]次に、二次構造は3D構造の断片に変換され、それらを結合してさらなる改良のための出発点を形成します。

  • より明白な方法は、3次元空間における全ての原子の座標を含む全原子表現を用いて一時的な3次元構造を表すことである。[ 45 ]このような構造の自由エネルギーは、分子動力学[ 46 ]やRNA特異的なコンフォメーションランドスケープのランダムサンプリング[ 47 ]などの既存の物理学に基づく方法を用いて改善することができ、その後、統計的スコアリングポテンシャルを用いてスクリーニングすることができる。[ 48 ]
  • 全原子構造は保存に大量のメモリを必要とし、改良には膨大な計算能力が必要となる。高解像度(全原子)精緻化段階に入る前に、まず粗粒度表現を用いて構造を真の形状に近づけることで[ 49 ]、高解像度段階での作業量を削減することができる。[ 50 ]

タンパク質フォールディングの問題と同様に、機械学習(ML)はRNA配列と2次元または3次元構造との間の高次関係を発見し、活用するためにも利用されてきました。MLはRNAフォールド全体の予測において従来の手法よりも優れていますが、2024年現在でも、分子内相互作用やリガンド結合部位のモデリングにおいては従来の手法(ML以外の手法)が依然として優位性を持っています。[ 51 ]

2025年12月、NvidiaはRNAの3D構造予測のためのMLベースのフレームワークであるRNAProをリリースしました。この手法は、スタンフォード大学が主催するKaggleコンペティションの上位入賞者の成果を統合したもので、入賞者はCASP16コンペティションの優勝者と統計誤差の範囲内のスコアを達成しました。CASP16コンペティションでは、主に物理学に基づいた従来の手法が人間の専門家によって指導されていました。これは MLモデルが従来のトップモデルの性能に匹敵する成果を達成した初めてのケースです。RNAProの手法は、テンプレートモデリング、多重配列アライメント(MSA)、事前学習済みのRNA基盤モデルRibonanzaNet2、およびAlphaFold3(Protenix)のオープンソース実装を組み合わせています。[ 53 ]

参照

参考文献

  1. ^ Ponce-Salvatierra, Almudena; Astha; Merdas, Katarzyna; Chandran, Nithin; Ghosh, Pritha; Mukherjee, Sunandan; Bujnicki, Janusz M (2019-01-22). 「実験データに基づくRNA 3D構造の計算モデリング」 . Bioscience Reports . 39 (2) BSR20180430. doi : 10.1042/bsr20180430 . ISSN  0144-8463 . PMC  6367127. PMID  30670629 .
  2. ^マグナス、マルシン;マテルスカ、ドロタ。ワハ、グジェゴシュ。チョジノフスキー、グジェゴシュ;ボニエキ、ミハル・J;プルタ、エルズビエタ。ウェイン・ドーソン。ドゥニン=ホルカヴィッツ、スタニスワフ。ブジニッキ、ヤヌシュ M (2014-04-23)。「実験的制約を利用した、RNA 3D 構造の計算モデリング」RNA 生物学11 (5): 522–536 .土井: 10.4161/rna.28826ISSN 1547-6286PMC 4152360PMID 24785264   
  3. ^ a b c Mathews DH (2006). 「RNA二次構造予測における革命」J. Mol. Biol . 359 (3): 526– 532. doi : 10.1016/j.jmb.2006.01.067 . PMID 16500677 . 
  4. ^ a b Mathews DH, Sabina J, Zuker M, Turner DH (1999). 「熱力学パラメータの拡張配列依存性がRNA二次構造の予測精度を向上させる」 . J Mol Biol . 288 (5): 911– 40. doi : 10.1006/jmbi.1999.2700 . PMID 10329189. S2CID 19989405 .  
  5. ^ Mathews DH, Disney MD, Childs JL, Schroeder SJ, Zuker M, Turner DH (2004). 「RNA二次構造予測のための動的計画法アルゴリズムへの化学修飾制約の組み込み」 . Proceedings of the National Academy of Sciences USA . 101 (19): 7287– 7292. Bibcode : 2004PNAS..101.7287M . doi : 10.1073/pnas.0401799101 . PMC 409911. PMID 15123812 .  
  6. ^ Andronescu, M; Condon, A; Hoos, HH; Mathews, DH; Murphy, KP (2007年7月1日). 「RNA二次構造予測のための効率的なパラメータ推定」.バイオインフォマティクス. 23 (13): i19-28. doi : 10.1093/bioinformatics/btm223 . PMID 17646296 . 
  7. ^ Chou, Fang-Chieh; Kladwang, Wipapat; Kappel, Kalli; Das, Rhiju (2016年7月26日). 「RNA最近接エネルギー予測のブラインドテスト」 . Proceedings of the National Academy of Sciences of the United States of America . 113 (30): 8430–5 . Bibcode : 2016PNAS..113.8430C . doi : 10.1073 / pnas.1523335113 . PMC 4968729. PMID 27402765 .  
  8. ^ Langdon, William B.; Petke, Justyna; Lorenz, Ronny (2018). 「RNAフォールド構造予測の進化と改善」.遺伝的プログラミング. コンピュータサイエンス講義ノート. 第10781巻. pp.  220– 236. doi : 10.1007/978-3-319-77553-1_14 . ISBN 978-3-319-77552-4
  9. ^ Raden, Martin; Müller, Teresa; Mautner, Stefan; Gelhausen, Rick; Backofen, Rolf (2020年12月). 「シード、アクセシビリティ、相互作用制約がsRNA標的予測に与える影響 - 体系的な評価」 . BMC Bioinformatics . 21 (1): 15. doi : 10.1186/s12859-019-3143-4 . PMC 6956497. PMID 31931703 .  
  10. ^ Zuker M.; Sankoff D. (1984). 「RNA二次構造とその予測」. Bull. Math. Biol . 46 (4): 591– 621. doi : 10.1016/s0092-8240(84)80062-2 (2025年7月1日現在非アクティブ). S2CID 189885784 . {{cite journal}}: CS1 maint: DOIは2025年7月時点で非アクティブです(リンク
  11. ^ a b c d Nussinov R, Piecznik G, Grigg JR and Kleitman DJ (1978)ループマッチングのアルゴリズム. SIAM Journal on Applied Mathematics.
  12. ^ a b Nussinov R, Jacobson AB (1980). 「一本鎖RNAの二次構造を予測するための高速アルゴリズム」 . Proc Natl Acad Sci USA . 77 (11): 6309–13 . Bibcode : 1980PNAS...77.6309N . doi : 10.1073/pnas.77.11.6309 . PMC 350273. PMID 6161375 .  
  13. ^ Zuker M, Stiegler P (1981). 「熱力学と補助情報を用いた大規模RNA配列の最適コンピュータフォールディング」 . Nucleic Acids Res . 9 (1): 133–48 . doi : 10.1093 / nar/9.1.133 . PMC 326673. PMID 6163133 .  
  14. ^ a b Rivas E, Eddy SR (1999). 「擬似ノットを含むRNA構造予測のための動的計画法アルゴリズム」J Mol Biol . 285 (5): 2053–68 . arXiv : physics/9807048 . doi : 10.1006/jmbi.1998.2436 . PMID 9925784 . S2CID 2228845 .  
  15. ^ Reeder J.; Giegerich R. (2004). 「熱力学に基づく実用的な擬似ノットフォールディングアルゴリズムの設計、実装、および評価」. BMC Bioinformatics . 5 104. doi : 10.1186/1471-2105-5-104 . PMC 514697. PMID 15294028 .  
  16. ^ a b Lyngsø RB, Pedersen CN (2000). 「エネルギーベースモデルによるRNAシュードノット予測」J Comput Biol . 7 ( 3–4 ): 409– 427. CiteSeerX 10.1.1.34.4044 . doi : 10.1089/106652700750050862 . PMID 11108471 .  
  17. ^ a b
  18. ^ Zuker M (2003). 「核酸のフォールディングとハイブリダイゼーション予測のためのMfoldウェブサーバー」 . Nucleic Acids Research . 31 (13): 3406– 3415. doi : 10.1093/nar/gkg595 . PMC 169194. PMID 12824337 .  
  19. ^ McCaskill JS (1990). 「RNA二次構造における平衡分配関数と塩基対結合確率」. Biopolymers . 29 ( 6–7 ): 1105–19 . doi : 10.1002/bip.360290621 . hdl : 11858 / 00-001M-0000-0013-0DE3-9 . PMID 1695107. S2CID 12629688 .  
  20. ^ Ding Y , Lawrence CE (2003). 「RNA二次構造予測のための統計的サンプリングアルゴリズム」 . Nucleic Acids Res . 31 (24): 7280–301 . doi : 10.1093/nar/gkg938 . PMC 297010. PMID 14654704 .  
  21. ^ Gardner PP; Giegerich, Robert (2004). 「RNA構造予測における比較アプローチの包括的比較」 . BMC Bioinformatics . 5 (1) 140. Bibcode : 2004BMCBi...5..140G . doi : 10.1186/1471-2105-5-140 . PMC 526219. PMID 15458580 .  
  22. ^ a b c d Bernhart SH, Hofacker IL (2009). 「コンセンサス構造予測からRNA遺伝子発見へ」 . Brief Funct Genomic Proteomic . 8 (6): 461– 71. doi : 10.1093/bfgp/elp043 . PMID 19833701 . 
  23. ^ Knudsen B, Hein J (2003). 「Pfold: 確率的文脈自由文法を用いたRNA二次構造予測」 . Nucleic Acids Res . 31 (13): 3423–8 . doi : 10.1093/nar/ gkg614 . PMC 169020. PMID 12824339 .  
  24. ^ a b Hofacker IL, Fekete M, Stadler PF (2002). 「アラインメントされたRNA配列の二次構造予測」J Mol Biol . 319 (5): 1059– 66. CiteSeerX 10.1.1.73.479 . doi : 10.1016/S0022-2836(02)00308-X . PMID 12079347 .  
  25. ^ Malik, A; Zhang, L; Gautam, M; Dai, N; Li, S; Zhang, H; Mathews, DH; Huang, L (2024年9月1日). 「LinearAlifold: RNAアライメントのための線形時間コンセンサス構造予測」 . Journal of Molecular Biology . 436 (17) 168694. arXiv : 2206.14794 . doi : 10.1016/j.jmb.2024.168694 . PMC 11377157. PMID 38971557 .  
  26. ^ Ruan, J., Stormo, GD & Zhang, W. (2004) ILM: 擬似ノットを用いたRNA二次構造予測のためのウェブサーバー. Nucleic Acids Research, 32(Web Server号), W146-149.
  27. ^ Shapiro BAとZhang K (1990)「ツリー比較を使用した複数のRNA二次構造の比較」Computer Applications in the Biosciences、vol. 6、no. 4、pp. 309–318。
  28. ^ Siebert, S; Backofen, R (2005年8月15日). 「MARNA:配列構造比較に基づくRNAの多重アライメントとコンセンサス構造予測」.バイオインフォマティクス. 21 (16): 3352–9 . doi : 10.1093/bioinformatics/bti550 . PMID 15972285 . 
  29. ^ Sankoff D (1985). 「RNAのフォールディング、アライメント、プロトシーケンス問題の同時解」. SIAM Journal on Applied Mathematics . 45 (5): 810– 825. CiteSeerX 10.1.1.665.4890 . doi : 10.1137/0145048 . 
  30. ^ a b Hofacker IL, Bernhart SH, Stadler PF (2004). 「RNA塩基対形成確率行列のアライメント」 .バイオインフォマティクス. 20 (14): 2222–7 . doi : 10.1093/bioinformatics/bth229 . PMID 15073017 . 
  31. ^ Havgaard JH, Lyngso RB, Stormo GD, Gorodkin J (2005). 「配列類似度が40%未満のRNA配列のペアワイズ局所構造アライメント」 .バイオインフォマティクス. 21 (9): 1815–24 . doi : 10.1093/bioinformatics/bti279 . PMID 15657094 . 
  32. ^ Torarinsson E, Havgaard JH, Gorodkin J. (2007) RNA配列の多重構造アライメントとクラスタリング.バイオインフォマティクス.
  33. ^ Mathews DH, Turner DH (2002). 「Dynalign: 2つのRNA配列に共通する二次構造を見つけるアルゴリズム」J Mol Biol . 317 (2): 191– 203. doi : 10.1006/jmbi.2001.5351 . PMID 11902836 . 
  34. ^ Harmanci AO、Sharma G、Mathews DH、(2007)、 Dynalignでの確率的アライメント制約を使用した効率的なペアワイズRNA構造予測、BMC Bioinformatics、8(130)。
  35. ^ Holmes I. (2005) RNA構造進化の加速確率推論BMCバイオインフォマティクス 2005年3月24日;6:73.
  36. ^ Kiryu H, Tabei Y, Kin T, Asai K (2007). 「Murlet: 構造RNA配列のための実用的なマルチアライメントツール」 .バイオインフォマティクス. 23 (13): 1588– 1598. doi : 10.1093/bioinformatics/btm146 . PMID 17459961 . 
  37. ^ Will, S; Joshi, T; Hofacker, IL; Stadler, PF; Backofen, R (2012年5月). 「LocARNA-P:構造RNAの正確​​な境界予測と検出精度の向上」 . RNA . 18 ( 5): 900–14 . doi : 10.1261/rna.029041.111 . PMC 3334699. PMID 22450757 .  
  38. ^ Sorescu, DA; Möhl, M; Mann, M; Backofen, R; Will, S (2012年7月). 「CARNA - RNA構造アンサンブルのアライメント」 . Nucleic Acids Research . 40 (Web Server版): W49-53. doi : 10.1093/nar/ gks491 . PMC 3394245. PMID 22689637 .  
  39. ^ Shapiro BA、Yingling YG、Kasprzak W、Bindewald E. (2007) RNA 構造予測のギャップを埋める。現在の構造バイオル。
  40. ^ Major F, Turcotte M, Gautheret D, Lapalme G, Fillion E, Cedergren R (1991年9月). 「RNAの3次元モデリングにおける記号計算と数値計算の融合」. Science . 253 (5025): 1255–60 . Bibcode : 1991Sci...253.1255F . doi : 10.1126/science.1716375 . PMID 1716375 . 
  41. ^ Major F, Gautheret D, Cedergren R (1993年10月). 「構造的制約からtRNA分子の3次元構造を再現する」 . Proc Natl Acad Sci USA . 90 (20): 9408–12 . Bibcode : 1993PNAS...90.9408M . doi : 10.1073 / pnas.90.20.9408 . PMC 47577. PMID 8415714 .  
  42. ^ Frellsen J, Moltke I, Thiim M, Mardia KV, Ferkinghoff-Borg J, Hamelryck T (2009). 「RNAコンフォメーション空間の確率モデル」 . PLOS Comput Biol . 5 (6) e1000406. Bibcode : 2009PLSCB...5E0406F . doi : 10.1371/journal.pcbi.1000406 . PMC 2691987. PMID 19543381 .  
  43. ^ Rother, Magdalena; Rother, Kristian; Puton, Tomasz; Bujnicki, Janusz M. (2011-02-07). 「ModeRNA:RNA 3D構造の比較モデリングツール」 . Nucleic Acids Research . 39 (10): 4007– 4022. doi : 10.1093 / nar/gkq1320 . ISSN 1362-4962 . PMC 3105415. PMID 21300639 .   
  44. ^ Neocles B Leontis、Eric Westhof編 (2012). RNA 3D構造解析と予測. Springer. ISBN 978-3-642-25740-7. OCLC  795570014 .
  45. ^ Zhao, Chenhan; Xu, Xiaojun; Chen, Shi-Jie (2017)、「VfoldによるRNA構造の予測」、Functional Genomics、Methods in Molecular Biology、vol. 1654、Springer New York、pp.  3– 15、doi : 10.1007/978-1-4939-7231-9_1ISBN 978-1-4939-7230-2PMC  5762135PMID  28986779
  46. ^ Vangaveti, Sweta; Ranganathan, Srivathsan V.; Chen, Alan A. (2016-10-04). 「RNA分子動力学の進歩:RNA力場シミュレータガイド」. Wiley Interdisciplinary Reviews: RNA . 8 (2) e1396. doi : 10.1002/wrna.1396 . ISSN 1757-7004 . PMID 27704698. S2CID 35501632 .   
  47. ^ Chen, Shi-Jie (2008年6月). 「RNAフォールディング:コンフォメーション統計、フォールディング速度論、イオン静電気」 . Annual Review of Biophysics . 37 (1): 197– 214. doi : 10.1146/annurev.biophys.37.032807.125957 . ISSN 1936-122X . PMC 2473866. PMID 18573079 .   
  48. ^ Laing, Christian; Schlick, Tamar (2011年6月). 「RNA構造予測、解析、設計への計算論的アプローチ」. Current Opinion in Structural Biology . 21 (3): 306– 318. doi : 10.1016/j.sbi.2011.03.015 . ISSN 0959-440X . PMC 3112238. PMID 21514143 .   
  49. ^ボニエツキ、ミハル・J.ラック、グジェゴシュ;ドーソン、ウェイン K.トマラ、コンラッド。ルカシュ、パヴェル。ソルティシンスキー、トマシュ。クリスチャン・M・ロザー; Bujnicki、Janusz M. (2015-12-19)。「SimRNA: RNA フォールディング シミュレーションと 3D 構造予測のための粗粒化手法」核酸研究44 (7): e63。土井10.1093/nar/gkv1479ISSN 0305-1048PMC 4838351PMID 26687716   
  50. ^ Stasiewicz, Juliusz; Mukherjee, Sunandan; Nithin, Chandran; Bujnicki, Janusz M. (2019-03-21). 「QRNAS:核酸構造の精密化のためのソフトウェアツール」 . BMC Structural Biology . 19 (1): 5. doi : 10.1186/s12900-019-0103-1 . ISSN 1472-6807 . PMC 6429776. PMID 30898165 .   
  51. ^ Nithin, Chandran; Kmiecik, Sebastian; Błaszczyk, Roman; Nowicka, Julita; Tuszyńska, Irina (2024-06-25). 「RNA 3D構造予測手法の比較分析:RNA–リガンド相互作用の強化モデリングに向けて」 . Nucleic Acids Research . 52 (13): 7465– 7486. doi : 10.1093/nar / gkae541 . ISSN 0305-1048 . PMC 11260495. PMID 38917327 .   
  52. ^ Zhang, Sicheng; Li, Jun; Zhou, Yuanzhe; Chen, Shi-Jie (2026年1月). 「CASP16におけるRNA 3D構造予測の強化:物理ベースモデリングと機械学習統合による予測精度の向上」 . Proteins: Structure, Function, and Bioinformatics . 94 (1): 239– 248. doi : 10.1002/prot.26856 . PMC 12354339. PMID 40488225 .  
  53. ^ Lee, Youhan; He, Shujun; Oda, Toshiyuki; Rao, G. John; Kim, Yehyun; Kim, Raehyun; et al. (2025-12-30). 「ブラインドコードコンペティションによるテンプレートベースのRNA構造予測の進歩」bioRxiv 10.64898/2025.12.30.696949 . 

さらに読む