[논문 리뷰] Steps and bumps: precision extraction of discrete states of molecular machines using physically-based, high-throughput time series analysis

Max A. Little, Bradley C. Steel외 6명|arXiv (Cornell University)|2010. 04. 07.

ATP Synthase and ATPases Research참고 문헌 31인용 수 29

한 줄 요약

이 논문은 잡음이 있는 회전 또는 평행 이동 데이터로부터 분자의 기계에서 이산 상태를 정확하고 자동으로 추출할 수 있도록 하는 물리적 모델 기반의 고속 시계열 분석 프레임워크를 소개한다. 자기회귀 모델링과 최적 필터링을 푸리에 기반의 볼록부각 탐지와 융합함으로써, 기존 방법에 비해 정확도, 속도, 상관 잡음에 대한 강건성 면에서 뛰어나며, 박테리아의 편모 운동기와 F1-ATPase에서 숨겨진 대칭성과 비 포isson성 동역학을 성공적으로 규명한다.

ABSTRACT

We report new statistical time-series analysis tools providing significant improvements in the rapid, precision extraction of discrete state dynamics from large databases of experimental observations of molecular machines. By building physical knowledge and statistical innovations into analysis tools, we demonstrate new techniques for recovering discrete state transitions buried in highly correlated molecular noise. We demonstrate the effectiveness of our approach on simulated and real examples of step-like rotation of the bacterial flagellar motor and the F1-ATPase enzyme. We show that our method can clearly identify molecular steps, symmetries and cascaded processes that are too weak for existing algorithms to detect, and can do so much faster than existing algorithms. Our techniques represent a major advance in the drive towards automated, precision, highthroughput studies of molecular machine dynamics. Modular, open-source software that implements these techniques is provided at http://www.eng.ox.ac.uk/samp/members/max/software/

연구 동기 및 목표

• 상관된 열 잡음에 의해 오염된 분자의 기계에서 이산 상태 전이를 탐지하는 데 있어 기존 알고리즘의 한계를 극복하기 위해.
• 단일 분자 실험에서 얻은 대규모 시계열 데이터에 대한 고속, 자동화된 분석을 가능하게 하기 위해.
• 분자의 기계에서 숨겨진 대칭성, 체류 시간 분포, 다중 단계 반응 동역학을 정확히 복원하기 위해.
• 광범위한 생물물리학 및 시스템 생물학 분야에서의 활용을 위해 계산적으로 효율적이고 오픈소스 소프트웨어 프레임워크를 개발하기 위해.

제안 방법

• 이 방법은 물리적 분자의 운동 모델과 자기상관 구조를 통합한 L1-PWC 및 L1-PWC-AR1 필터를 사용하여 고정 윈도우 없이 시계열을 스무딩한다.
• 자기상관 잡음(AR1)을 포함한 일반화된 랑주반 역학을 사용하여 기반 역학을 모델링함으로써, 단계 전이를 유지하면서 잡음 억제를 향상시킨다.
• 스무딩된 신호에 대한 푸리에 계수 선택을 통해 볼록부각 탐지를 수행함으로써 히스토그램이나 커널 밀도 추정을 피하고 강건한 대칭성 탐지가 가능해진다.
• 스무딩된 시계열의 푸리에 스펙트럼에서 피크를 식별하여 상태 분포의 주요 주기성을 탐지한다.
• 추정된 상태 분포의 각 피크에 가장 가까운 시간 포인트에 각 시간 포인트를 할당하여 이산 상태를 분류한다.
• 이 접근법은 모듈식이며 오픈소스이며, 단일 분자 데이터의 실시간 및 고속 스트림 분석을 위한 소스 코드가 제공된다.

실험 결과

연구 질문

• RQ1 고도로 상관된 열 잡음에 묻혀 있는 약한 이산 상태 전이를 분자의 기계에서 탐지할 수 있는가?
• RQ2 변동하는 체류 시간과 잡음이 존재하는 조건에서 분자의 상태 전이의 대칭성과 주기성을 정확히 복원할 수 있는가?
• RQ3 시계열 분석을 통해 분자의 기계에서 비 포isson성 스텝 동작과 다중 단계 제한 속도 단계를 식별할 수 있는가?
• RQ4 잡음 또는 전이 지속 시간에 대한 사전 가정 없이도 대규모 데이터셋에서 고속, 자동화, 정밀한 상태 추출을 달성할 수 있는가?

주요 결과

• 모의 편모 운동기 데이터에서 L1-PWC-AR1 필터는 기존의 고전적 필터 대비 평균 절대 오차(MAE)를 최대 50% 감소시키고 상대 절대 거칠기(RAR)를 최대 70% 감소시켰다.
• ECF-Bump 방법은 20%의 체류 위치 비대칭성 또는 감마 분포를 가진 체류 시간이 존재하는 경우에도 95%의 시뮬레이션 테스트 케이스에서 진짜 대칭성(예: 26, 30, 40 상태)을 정확히 식별했다.
• 이 방법은 박테리아의 편모 운동기에서 비 포isson성 스텝 동작을 탐지하여 F1-ATPase에서 다중 단계 제한 속도 단계의 명확한 패턴을 높은 통계적 신뢰도로 드러냈다.
• 실제 데이터에서 알고리즘은 수천 개의 각도-시간 트레이스를 빠르게 처리했으며, 편모 운동기의 26단계 대칭성을 모순 없이 규명했다.
• 히스토그램-FFT 및 커널 밀도 추정과 피크 탐지 방법에 비해 볼록부각 탐지 방법이 모든 테스트 케이스에서 슈퍼어리어를 기록했으며, 특히 고잡음 또는 비 i.i.d. 체류 시간 조건에서 뛰어난 성능을 보였다.
• 오픈소스 소프트웨어 구현은 사용자 간섭 최소화로 재현 가능하고 고속 스트림 분석을 통해 단일 분자 동역학을 분석할 수 있도록 했다.