[論文レビュー] Modeling Bacterial DNA: Simulation of Self-avoiding Supercoiled Worm-Like Chains Including Structural Transitions of the Helix

Thibaut Lepage, Ivan Junier|arXiv (Cornell University)|Aug 7, 2017

DNA and Nucleic Acid Chemistry参考文献 43被引用数 5

ひとこと要約

本論文は、局所的なヘリカル構造遷移（例：変性バブルやZ-DNA）を明示的に含む自己回避的スーパокоイルド・ワーム・ライク・チェイン（ssWLC）としての細菌DNAをモデル化するためのモンテカルロシミュレーションフレームワークを提示する。この手法は、離散的ポリマー・モデリングと熱力学的モンテカルロサンプリングを組み合わせ、10–20 bp解像度で折りたたみと構造遷移をシミュレートし、負のスーパakoイルド状態下でのプレクトネーム形成とヘリックス遷移の研究を可能にする。

ABSTRACT

Under supercoiling constraints, naked DNA, such as a large part of bacterial DNA, folds into braided structures called plectonemes. The double-helix can also undergo local structural transitions, leading to the formation of denaturation bubbles and other alternative structures. Various polymer models have been developed to capture these properties, with Monte-Carlo (MC) approaches dedicated to the inference of thermodynamic properties. In this chapter, we explain how to perform such Monte-Carlo simulations, following two objectives. On one hand, we present the self-avoiding supercoiled Worm-Like Chain (ssWLC) model, which is known to capture the folding properties of supercoiled DNA, and provide a detailed explanation of a standard MC simulation method. On the other hand, we explain how to extend this ssWLC model to include structural transitions of the helix.

研究の動機と目的

• 細菌DNAにおけるスーパakoイルド誘発折りたたみ（プレクトネーム形成）と局所的構造遷移（例：変性バブル）を両方とも捉える計算モデルの開発を目的とする。
• 標準的な自己回避的スーパakoイルド・ワーム・ライク・チェイン（ssWLC）モデルを、負にスーパakoイルドされたDNAに関連するヘリカル構造遷移を含めるように拡張することを目的とする。
• 生理的条件下での、数kbペアに達する長さのDNA分子のマルチスケールシミュレーションを可能にする。これには、スーパakoイルドと構造的不均一性が含まれる。
• プレクトネーム形成と局所的DNA構造変化の競合を研究できる、実用的で実装可能なモンテカルロシミュレーションプロトコルを提供することを目的とする。
• 衝突検出、リンク数保存、けい放形成といった主要な課題に対処することで、数値的安定性を確保することを目的とする。

提案手法

• DNA鎖をN個のセグメントに離散化し、N個の接続点を介して結合させ、自己回避性を確保するため各セグメントを半径reの不透光なシリンダーとして表現する。
• 局所的な接線ベクトルと法線ベクトルを用いて、連続するセグメント間の曲げ角（θi）とねじれ角（φi）を定義し、B-DNAのPersistence lengthに一致するように離散的曲げモジュラスKdを調整する。
• 剛体的セグメント群（サイズM）を最大角度αだけ回転させるモンテカルロ移動を用い、全エネルギー（曲げ、ねじれ、自己回避ペナルティ）に基づいてメトロポリス・ハスティングス基準に従って受け入れる。
• 移動中のシリンダーと静的鎖セグメントとの間の交差を正確に検出するアルゴリズムを実装し、近接接近を交差とみなしてリンク数の誤りを防ぐ。
• ねじれと曲げの分解を用いてリンク数保存制約を適用し、回転中の位相的整合性を保証する。
• x = -1におけるアレクサンダー多項式を部分的不変量として用い、ねじれ構造を検出し、誤検出を最小限に抑える。

実験結果

リサーチクエスチョン

• RQ1 自己回避的スーパakoイルド・ワーム・ライク・チェインモデルは、細菌DNAにおける局所的構造遷移（例：変性バブルやZ-DNA）をどのように拡張できるか？
• RQ2 モンテカルロシミュレーションにおいて、受容率と構造的サンプリング効率の両立を図るための最適なブロックサイズと回転角度のバランスは何か？
• RQ3 生理的条件下におけるスーパakoイルドDNAの離散的ポリマー・モデルで、自己回避性とリンク数保存をどのように正確に実装できるか？
• RQ4 浮動小数点誤差に起因する誤ったリンク数変更を防ぐために、衝突検出における誤検出を最小限に抑える計算戦略は何か？
• RQ5 プレクトネーム形成と局所的構造遷移の競合を、10–20 bp解像度スケールで定量的にモデル化・シミュレートするにはどうすればよいか？

主な発見

• 離散的ssWLCモデルは、自己回避性とPersistence lengthに依存するプレクトネーム半径によって、負のスーパakoイルド下でのスーパakoイルドDNAの折りたたみを成功裏に再現した。
• 変性バブルやZ-DNAといった構造遷移の導入により、負にスーパakoイルドされたDNAにおける局所的ヘリカル不均一性のシミュレーションが可能になった。
• M ≈ 20ℓp および α ≈ 30° のブロック回転モンテカルロ法が、受容率と構造的サンプリング効率の両立に適していることが判明した。
• 衝突検査におけるニア・ミス検出（近接接近を交差とみなす）が、浮動小数点誤差に起因する誤ったリンク数変更やねじれ構造の生成を効果的に防止した。
• x = -1におけるアレクサンダー多項式の使用により、非自明なねじれ構造の効率的検出が可能となり、非ねじれ構造は誤分類される可能性が極めて低いことが確認された。
• シミュレーションフレームワークは、円形（例：プラスミド）および線形DNA分子の両方に対して安定しており、線形鎖の場合は適切な初期化と壁制約を適用することで成立する。