A General Computational Framework for Sequence-Level Thermodynamic Analysis of Modeled RNA–DNA and DNA–DNA Duplex States: A Proof-of-Concept Study in HIV-1

技術要約：モデル化されたRNA–DNAおよびDNA–DNA二重鎖状態の配列レベルの熱力学的解析のための一般的計算フレームワーク

問題と範囲
本研究は、ウイルスゲノムにおけるモデル化されたRNA–DNA二重鎖が、DNA–DNA二重鎖と比較して熱力学的にどの程度有利であるかを評価するための、再現可能な配列レベルの計算フレームワークの必要性に対処するものである。Rループ（RNAが一方のDNA鎖とハイブリダイズし、他方の鎖を排除する三本鎖構造）はゲノム生物学における確立された特徴であるが、直接的な配列レベルの熱力学計算は、細胞内でのRループ形成の証拠と混同されることが多い。著者は、このような計算は、クロマチンアクセシビリティ、転写状態、およびタンパク質による競合に依存する直接的な細胞内測定と同等に扱うべきではないと主張している。本研究では、HIV-1をプルーフ・オブ・コンセプト（概念実証）システムとして用い、生物学的な統合やRループ形成を証明すると主張することなく、実験的なフォローアップのための潜在的なウイルス配列コンテキストを優先順位付けするためのパイプラインを確立する。

方法論
本フレームワークはHXB2参照ゲノム（K03455.1）上で動作し、16個の定義済みゲノム領域（TAR、PBS、gag、pol、envなどの機能的要素を含む18〜132ヌクレオチドの範囲）を評価する。手法は以下の5段階のパイプラインで構成される：

1. 配列抽出： ロックされた参照ゲノムからDNA配列を抽出する。TをUに変換することでRNA表現を生成し、その相補的なDNA鎖をリバースコンプリメントとして生成する。
2. RNA–DNAモデリング： 特定のパラメータ（-d：オリゴに対するDNAパラメータ、および -m 3：ターゲット構造を無視する）を用いて、OligoWalk (RNAstructure 6.6) を使用してRNA–DNA二重鎖状態を計算する。出力される最小自由エネルギーを $\Delta G_{RD}$ と表記する。
3. DNA–DNA参照： DuplexFold (RNAstructure 6.6) を使用し、2本のDNA鎖に対するDNAパラメータ（-d）を用いて、二分子のDNA–DNA参照状態を計算する。出力される二分子MFEを $\Delta G_{DD}$ と表記する。
4. モデル定義による対比： 研究固有の熱力学的対比を $\Delta\Delta G_{model} = \Delta G_{RD} - \Delta G_{DD}$ と定義する。著者は、これら2つの観測値は異なるアルゴリズムから導出されたものであり、普遍的な熱力学定数ではないため、これが「モデル定義」の対比であることを明示している。
5. 統計解析： 本研究では、ゼロに対する両側一標本t検定、正確符号検定、および符号反転置換検定を採用している。感度分析では、重複領域（TARおよびR_5LTR）を除外している。また、パネル平均に対してブートストラップ信頼区間を算出している。

主な結果

• 熱力学的パターン： 16個の定義済み領域のうち、15個が負の $\Delta\Delta G_{model}$ 値を示した。これは、指定された手順の下では、モデル化されたRNA–DNA二重鎖が、対応するDNA–DNA二重鎖よりも熱力学的に有利（自由エネルギーがより負）であることを示している。U5領域のみが正の対比（+0.9 kcal/mol）を示した。
• 統計的有意性： パネルの平均対比は −8.00 kcal/mol (SD 5.45) であった。ゼロに対する両側一標本t検定の結果、$t_{15} = -5.87$、p値は $1.54 \times 10^{-5}$ であった。
• 堅牢性： 方向性の結果（負の対比）は、TAR、R_5LTR、またはその両方の重複領域を除外した感度分析においても維持された。パネル平均のブートストラップ95%信頼区間は [−10.55, −5.39] kcal/mol であった。
• 配列との関連性： 探索的解析により、正規化された対比（$\Delta\Delta G_{model}$/nt）とGC含有量の間に統計的に有意な負の相関（$r = -0.551, p = 0.027$）が認められた。これは、本モデルにおいてGC含有量が高いほど、より負の対比を示すことを示唆している。生の値に関する配列長との相関は、統計的に有意ではなかった。

意義および主張
本論文は、自らの貢献を生物学的な証明ではなく、検証可能な計算による予測として位置づけている。

• フレームワークの有用性： 主要な貢献は、RNA–DNAハイブリダイゼーションがDNA–DNAペアリングよりも熱力学的に有利となる配列コンテキストを特定するための、他のウイルスゲノムにも適用可能な再現可能なパイプラインであることである。
• 控えめな生物学的解釈： 著者は、これらの結果が細胞内のRループ形成、ゲノム統合部位、または因果的な生物学的メカニズムを確立するものではないことを明示している。本研究は、クロマチン・トポロジー、転写キネティクス、またはタンパク質による競合を考慮していない。
• 優先順位付けのシグナル： 本フレームワークは、制御されたRNA–DNAハイブリダイゼーション実験や細胞内Rループアッセイを検討すべき配列コンテキストを特定するための、優先順位付けツールとして機能する。
• 限界： 本研究は単一の参照ゲノム（HXB2）に限定されており、非独立な領域を使用し、比較される2つの状態に対して異なる計算手順に依存している。

結論として、本研究は、特定のRNAstructure 6.6の手順の下で、HIV-1 HXB2配列パネルが、DNA–DNA二重鎖よりもRNA–DNA二重鎖を支持する一貫した、統計的に有意なモデル化された熱力学的対比を示すことを実証している。これは、生体内（in vivo）のRループ生物学の確認ではなく、将来の経験的検証のための配列レベルの仮説として機能するものである。