[論文レビュー] Statistical Mechanics of Semiflexible Chains: A Meanfield Variational Approach

D. Thirumalai, Bae‐Yeun Ha|ArXiv.org|May 20, 1997

Protein Structure and DynamicsBiochemistry, Genetics and Molecular Biology参考文献 7被引用数 21

ひとこと要約

本論文は、局所的な接ベクトルの制約を、計算が可能な統計力学的計算を可能にするために、全体の平均値制約に置き換える平均場変分法を導入し、半剛直ポリマー鎖のモデル化を行う。この手法は、DNA におけるシミュレーションや実験と定量的に一致する終点間距離分布および引張り応力下での弾性応答を正確に予測するが、対称性の破れを伴う強いネマチック場では、パーサシステンス長の過小評価を示す。これは、対称性の破れを捉えられていないためである。

ABSTRACT

We describe a simple meanfield variational approach to study a number of properties of intrinsically stiff chains which are appropriate models for a large class of biopolymers. We present the calculation of the distribution of end-to-end distance and the elastic response of stiff chains under tension using this approach. In the former example we find that the simple expression almost quantitatively fits the results of computer simulation. For the case of the stiff chain under tension we recover analytically all the known limits. We obtain quantitative agreement with recent experiments on the stretching of DNA. The limitations of our approach are also discussed.

研究の動機と目的

• DNA やバイオポリマーのような本質的に剛直なポリマーのための計算可能な統計力学フレームワークの開発。
• ワームラインチャインモデルにおける局所的接ベクトル制約を強制する困難さの解決。
• 終点間距離分布や引張り応力下での弾性応答といった、構造的性質の計算。
• DNA の引張り実験とシミュレーションデータに対する、この手法の精度の検証。
• 対称性の破れを伴う状況、例えば強いネマチック場における限界の特定。

提案手法

• ||u(s)|| = 1 という局所的制約を、接ベクトルの大きさに関する全体的平均値制約に置き換える。
• 真の統計力学的挙動を近似するために、変分自由エネルギー関数を用いる。
• 定常位相近似を適用し、有効なパーサシステンス長および相関関数を導出する。
• 作用に二次相互作用項を含めることで、外部場（例：ネマチック場や引張力）を組み込む。
• 弱場および強場のさまざまな極限において、得られた方程式を解析的に解き、物理的観測量を抽出する。
• 終点間距離および弾性に関する既知の解析的極限とシミュレーションデータとを比較して、結果の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 パーサシステンス長や終点間距離分布といった、主要な物理的性質を保持しつつ、半剛直鎖の統計力学をどのように簡略化できるか？
• RQ2 引張力下での剛直鎖の弾性応答は何か？ また、平均場アプローチは既知の解析的極限をどれほど正確に再現できるか？
• RQ3 ネマチック場は、半剛直鎖の構造的性質、特にパーサシステンス長にどのように影響を与えるか？
• RQ4 なぜ平均場変分アプローチは、強いネマチック場下でのパーサシステンス長の指数的増加を捉えられていないのか？
• RQ5 この手法は、僅かに曲げ剛性とパーサシステンス長をパラメータとして用いるだけで、実験的DNAの引張りデータを定量的に記述できるか？

主な発見

• 平均場変分法により、終点間距離分布を単純な式で得ることができ、シミュレーション結果とほぼ定量的に一致する。
• この手法は、引張力下での剛直鎖の弾性応答のすべての既知の極限を解析的に回復する。
• DNA の引張り実験では、僅かに曲げ剛性とパーサシステンス長を調整パラメータとして用いるだけで、実験データと定量的に一致する。
• 弱いネマチック場の極限では、場の方向に沿ってパーサシステンス長が (1 + 1/3 gl₀) の係数で線形に増加することが予測され、先行研究と一致する。
• 強いネマチック場の極限では、このアプローチは有限の増加（係数2）を予測するが、正確な結果では exp(√(glₚ)) のように指数関数的に増加する。これは、インスタントン効果を捉えられていないためである。
• この手法は、対称性の破れの領域で破綻する。これは、有効ポテンシャル内の degenerate な極小値間のトンネル効果を記述できないためであり、このような状況での限界を示している。