[論文レビュー] A robust approach to estimating rates from time-correlation functions

John D. Chodera, Phillip Elms+6 more|arXiv (Cornell University)|Aug 10, 2011

Force Microscopy Techniques and Applications参考文献 1被引用数 8

ひとこと要約

この論文は、分子動力学および単一分子実験における時間相関関数から力学的反応速度定数を推定するための頑健な手法を提示している。ノイズの多いデータの数値微分を回避することで、安定した速度定数推定が可能となる。この手法は反応フラックス理論に基づく仮想時間スケール法を用い、さまざまなサンプリング周波数や分離面の選択に対しても安定した推定結果を示し、RNAヘアピンの力分光法データにおいても成功が確認された。

ABSTRACT

While seemingly straightforward in principle, the reliable estimation of rate constants is seldom easy in practice. Numerous issues, such as the complication of poor reaction coordinates, cause obvious approaches to yield unreliable estimates. When a reliable order parameter is available, the reactive flux theory of Chandler allows the rate constant to be extracted from the plateau region of an appropriate reactive flux function. However, when applied to real data from single-molecule experiments or molecular dynamics simulations, the rate can sometimes be difficult to extract due to the numerical differentiation of a noisy empirical correlation function or difficulty in locating the plateau region at low sampling frequencies. We present a modified version of this theory which does not require numerical derivatives, allowing rate constants to be robustly estimated from the time-correlation function directly. We compare these approaches using single-molecule force spectroscopy measurements of an RNA hairpin.

研究の動機と目的

• 分子動力学および単一分子実験におけるノイズの多い有限時間の相関関数から、信頼性の高い速度定数を推定する課題に対処すること。
• 数値微分に依存する従来の手法がサンプリング周波数や測定ノイズに敏感であるという限界を克服すること。
• 明示的な微分やプラトー領域の特定を必要とせず、時間相関関数から直接速度定数を推定する手法を開発すること。
• 遷移経路理論および反応フラックス解析における分離面の選択に対する頑健性を向上させること。
• 理論的予測と実験結果の間の整合性を高めるために、多様な実験的・シミュレーション的条件下でも正確かつ一貫性のある速度定数推定を可能にすること。

提案手法

• 数値微分を回避するため、時間相関関数を直接用いる修正された反応フラックス理論を提案する。
• 時間相関関数の値が初期値の何分の一に減少するかに対応する時間スケールの逆数として定義される、仮想速度関数 $ k_{\mathrm{im}}(t) $ を導入する。
• 時間相関関数 $ C(t) = \langle h_A(0) h_A(t) \rangle $ を用いて、微分を必要とせずに $ k_{\mathrm{im}}(t) $ を計算する。
• RNAヘアピンの単一分子力分光法データにこの手法を適用し、さまざまなサンプリング周波数での性能を検証する。
• 状態AとBを定義するための二状態モデルを用い、分離面 $ x^{\ddagger} $ を導入し、$ h_A(t) $ および $ h_B(t) $ を指標関数とする。
• 平衡確率 $ \pi_A $, $ \pi_B $ を用いて、微視的反応速度定数を $ k_{A\rightarrow B} = \pi_B k $ および $ k_{B\rightarrow A} = \pi_A k $ として導出する。

実験結果

リサーチクエスチョン

• RQ1 ノイズの多いデータの数値微分が不安定な結果をもたらす場合、時間相関関数から速度定数を信頼性高く推定する方法は何か？
• RQ2 新しい手法は、実験的またはシミュレーション的データにおけるサンプリング周波数への感受性をどの程度低減するか？
• RQ3 反応フラックス関数のプラトー領域が特定しづらい場合でも、仮想速度関数 $ k_{\mathrm{im}}(t) $ が安定的かつ正確な速度定数推定を提供できるか？
• RQ4 分離面 $ x^{\ddagger} $ の選択に伴う速度定数推定の頑健性はどの程度か？
• RQ5 この手法は、バイオ分子反応速度論における理論的予測と実験測定の整合性を向上させられるか？

主な発見

• 仮想速度関数 $ k_{\mathrm{im}}(t) $ は1 kHzおよび50 kHzの両方のサンプリング周波数で安定した速度定数推定を示し、低周波数帯での反応フラックス法 $ k_{\mathrm{RF}}(t) $ よりも優れた性能を発揮した。
• 1 kHzのサンプリング周波数下では、$ k_{\mathrm{im}}(t) $ は $ \tau_{\mathrm{rxn}} $ の数倍にわたり安定したままであるのに対し、$ k_{\mathrm{RF}}(t) $ はノイズや分解能の低さにより早期に減少した。
• 仮想速度推定値 $ k_{\mathrm{im}}(\tau) $ は、単純な速度定数推定と比較して、分離面 $ x^{\ddagger} $ の選択に対して著しく感受性が低いことが示された。
• この手法は、さまざまなサンプリング周波数および分離面の配置に対しても頑健性を維持しており、実世界のデータに実用的であることが示された。
• 状態の占有率から推定された $ \pi_A $, $ \pi_B $ を用いて、$ k_{A\rightarrow B} = \pi_B k $ として微視的反応速度定数 $ k_{A\rightarrow B} $ および $ k_{B\rightarrow A} $ が一貫して導出された。
• この手法により、相関関数がノイズが多いか、間隔が広い場合でも信頼性の高い動力学的解析が可能となり、単一分子実験や長時間スケールのシミュレーションに適している。