[論文レビュー] Anisotropic Elastic Model for Short DNA Loops

G. Bijani, N. Hamedani Radja+2 more|ArXiv.org|May 3, 2006

DNA and Nucleic Acid Chemistry参考文献 1被引用数 4

ひとこと要約

本稿では、短いDNAループのための異方性弾性モデルを提案し、DNAが溝（ロール）に沿ってはるかに柔軟に曲がる（傾き方向より）という曲げの異方性が、鋭い端を持つ多角形のループ形状を誘導することを示している。主な結果は、長く制約のないDNAで用いられる調和平均とは異なり、有効な持久度は幾何平均 √(A₁A₂) に等しくなることである。これは、異方性そのものがミニサークルにおける高いループ形成確率を説明できないことを説明している。

ABSTRACT

Effect of bending anisotropy on a planar DNA loop, using energy minimization and neglecting entropic effects, is studied. We show that the anisotropy results in polygonal shape of the loop and increasing the anisotropy makes the edges sharper. Calculating the energy of such a loop lets us to find effective persistence length as the geometrical mean of hard and soft rigidities, which is quite different from harmonic mean for an unconstrained long DNA.

研究の動機と目的

• 曲げの異方性が平面的DNAループの形状およびエネルギーに与える影響を調査すること。
• 異方性のある曲げ剛性（A₁：傾き方向、A₂：ロール方向）がループの形状および安定性に与える影響を特定すること。
• エントロピー効果を無視してエネルギー最小化の下で短いDNAループの有効な持久度を導出すること。
• 異方性が実験的に観察されたミニサークルにおける高いサイクリゼーション確率を説明できるかどうかを評価すること。

提案手法

• DNAを異方性曲げ剛性を持つ平面的弾性ロッドとしてモデル化し、A₁とA₂を異なる曲げ方向（傾きとロール）に割り当てる。
• 平行で位相が一致するループ端の制約下でエネルギー最小化を用いて平衡状態のループ配置を特定する。
• 曲げとねじりの寄与を含む弾性エネルギー関数を導出する。ねじりエネルギーは「端が位相一致」する条件から生じる。
• 曲げエネルギーを円形ループのものと比較することで、有効な持久度を √(A₁A₂) として計算する。
• エネルギー式の妥当性を検証し、ねじり速度のゆらぎを推定するために等分配則を適用する。
• 数値推定を用いて、ねじり速度ωの変動が平均値から1%未満であることを確認し、モデルにおけるωの定数仮定が妥当であることを裏付ける。

実験結果

リサーチクエスチョン

• RQ1 エントロピー効果がない条件下で、曲げの異方性は短いDNAループの形状にどのように影響するか？
• RQ2 異方性曲げ剛性A₁およびA₂を有する短いDNAループの有効な持久度は何か？
• RQ3 短いループにおいて、幾何平均 √(A₁A₂) が調和平均よりもエネルギーの記述に適しているか？
• RQ4 異方性曲げだけがDNAミニサークルで観察される高いサイクリゼーション確率を説明できるか？
• RQ5 ループ内の多角形エッジの数は、ループ長および異方性比A₁/A₂とどのように関係するか？

主な発見

• 曲げの異方性により、DNAループは鋭い端を持つ多角形形状をとる。エッジは曲げ剛性が最小となる方向（ロール方向）に形成される。
• 多角形エッジの数は、ヘリカルターン数の2倍に等しく、ループ長に比例して増加する。
• 短いループの有効な持久度は幾何平均 √(A₁A₂) に等しく、長く制約のないDNAで用いられる調和平均（2/(1/A₁ + 1/A₂)）よりも大きい。
• ループエネルギーにおけるねじりエネルギー項は、「端が位相一致」する条件のため、L⁻¹の減衰因子を伴う振動的挙動を示す。
• 数値推定により、ねじり速度ωの変動は平均値から1%未満であることが判明し、モデルにおけるωの定数仮定が妥当であることが裏付けられた。
• 本モデルは、異方性そのものがミニサークルにおける高いループ形成確率を説明できないと予測している。なぜなら、√(A₁A₂) では実験的なJ係数を説明するには不十分であるからである。