[논문 리뷰] A robust approach to estimating rates from time-correlation functions

John D. Chodera, Phillip Elms외 6명|arXiv (Cornell University)|2011. 08. 10.

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} $ 가 일관되게 도출된다.
• 이 방법은 상관 함수가 노이즈가 있거나 희소하게 샘플링되는 경우에도 신뢰할 수 있는 동역학 분석을 가능하게 하여, 단일 분자 실험과 장시간 스케일 시뮬레이션에 적합하다.