[논문 리뷰] Identification and Analysis of Transition and Metastable Markov States

Linda Martini, Adam Kells외 3명|arXiv (Cornell University)|2016. 05. 13.

Protein Structure and DynamicsBiochemistry, Genetics and Molecular Biology참고 문헌 40인용 수 18

한 줄 요약

이 논문은 분자 동역학 궤적에서 전이 및 준안정 상태를 자동으로 식별하기 위한 새로운 알고리즘을 소개한다. 마코프 상태 모델(Markov state model) 접근법을 사용하여 집합 좌표의 집합으로부터 내재된 가장 느린 붕괴 속도를 캡처함으로써, 소형 주형(Ala5)과 더 큰 단백질(EGFR)에서 자유 에너지 장벽과 동적 경로를 정확하게 탐지한다. 이는 복잡한 생고분자 시스템에서 희귀한 구조적 전이를 연구하는 데 있어 강력하고 일반화 가능한 도구를 제공한다.

ABSTRACT

We present a new method that enables the identification and analysis of both transition and metastable conformational states from atomistic or coarse-grained molecular dynamics (MD) trajectories. Our algorithm is presented and studied by using both analytical and actual examples from MD simulations of the helix-forming peptide Ala5, and of a larger system, the epidermal growth factor receptor (EGFR) protein. In all cases, our method identifies automatically the corresponding transition states and metastable conformations in an optimal way, with the input of a set of relevant coordinates, by capturing accurately the intrinsic slowest relaxation rate. Our approach provides a general and easy to implement analysis method that provides unique insight into the molecular mechanism and the rare but crucial rate limiting conformational pathways occurring in complex dynamical systems such as molecular trajectories.

연구 동기 및 목표

• 분자 동역학 궤적에서 전이 및 준안정 상태를 식별하기 위한 일반적이고 자동화된 방법을 개발하는 것.
• 복잡한 생고분자 시스템에서 희귀하지만 동역학적으로 중요한 구조적 전이 경로를 탐지하는 데 도전하는 문제를 해결하는 것.
• 고차원 궤적 데이터로부터 자유 에너지 장벽과 가장 느린 붕괴 모드를 정확하게 식별하는 데 기여하는 것.
• 원자적 및 군집화된 시뮬레이션에 모두 적용 가능한 계산적으로 효율적이고 실행 가능한 프레임워크를 제공하는 것.

제안 방법

• 이 방법은 궤적 데이터의 구조적 동역학을 분석하기 위해 마코프 상태 모델(MSM) 프레임워크를 사용한다.
• 이 방법은 구조적 동역학의 운동성과 가장 느린 붕괴 모드에 기반하여 준안정 상태를 식별한다.
• 알고리즘은 상태 공간을 정의하고 전이 확률을 계산하기 위해 집합 좌표의 집합을 입력으로 사용한다.
• 전이 상태는 주로 반응 좌표에 위치하고 유속 축적량이 최소인 구성이 위치한 구조를 식별함으로써 탐지된다.
• 이 방법은 전이 행렬의 스펙트럼 성질을 활용하여 가장 느린 붕괴 속도를 추출하며, 이는 한계 속도 단계에 해당한다.
• 분석 모델과 실질적인 Ala5 및 EGFR의 MD 시뮬레이션을 통해 방법의 타당성을 검증하였으며, 강건성과 정확성을 입증하였다.

실험 결과

연구 질문

• RQ1 반응 좌표에 대한 사전 지식 없이 분자 동역학 궤적에서 전이 및 준안정 상태를 신뢰성 있게 식별할 수 있는가?
• RQ2 복잡한 생고분자 시스템에서 가장 느린 붕괴 모드가 한계 속도 결정 단계의 구조적 전이 경로를 어떻게 정의하는가?
• RQ3 기존 방법과 비교해 볼 때, 제안된 방법은 자유 에너지 장벽과 동적 중간체를 탐지하는 데 어떻게 성능을 발휘하는가?
• RQ4 고도로 복잡한 구조적 특성을 가진 시스템, 예를 들어 EGFR 단백질에서 전이 상태를 정확하게 식별할 수 있는가?
• RQ5 선택된 집합 좌표의 종류에 따라 이 방법의 성능이 얼마나 의존하는가?

주요 결과

• 이 방법은 알려진 헬릭스-코일 전이 행동과 일치하는 Ala5 펩타이드에서 전이 상태와 준안정 상태를 성공적으로 식별하였다.
• EGFR 단백질의 경우, 활성화 경로를 따라 뚜렷한 전이 상태를 탐지하였으며, 이는 실험 및 이전 시뮬레이션 결과와 일치하였다.
• 식별된 가장 느린 붕괴 속도는 내재된 동역학적 장벽과 매우 유사하여, 이 방법이 한계 속도 과정을 정확히 캡처하고 있음을 확인하였다.
• 다양한 집합 좌표 선택에 대해 높은 강건성을 확보하여 일반화 가능성을 입증하였다.
• 수동적인 편향 없이 자동으로 구조적 상태를 탐지하여 상태 할당의 주관성을 줄였다.
• 명확한 동적 경로와 잘 정의된 중간 상태를 드러내어 EGFR 활성화 메커니즘에 대한 기계적 통찰을 제공하였다.