[논문 리뷰] Weighted Ensemble Path Sampling for Multiple Reaction Channels

Bin W. Zhang, David Jasnow외 1명|ArXiv.org|2009. 02. 16.

Advanced Chemical Sensor Technologies인용 수 5

한 줄 요약

이 논문은 전이 경로 샘플링(Transition Path Sampling, TPS)과 같은 전통적 방법과 달리 단일 경로에 갇히지 않고 분자의 시스템에서 다수의 반응 채널을 효과적으로 샘플링할 수 있는 가중 평균(Weighted Ensemble, WE) 경로 샘플링 방법을 입증한다. 이는 500K와 300K에서 알라닌 다펩타이드에 대해 고품질의 전이 경로 집합을 생성함으로써, 브루트 포스 시뮬레이션으로는 비용이 너무 많이 들기 때문에 불가능한 실온 조건에서도 채널의 비중을 정량적으로 분석할 수 있게 한다.

ABSTRACT

Finding and sampling multiple reaction channels for molecular transitions remains an important challenge in physical chemistry. Here we show that the weighted ensemble (WE) path sampling method can readily sample multiple channels. In a first test, both the WE and transition path sampling methods are applied to two-dimensional model potentials. The comparison explains why the weighted ensemble approach will not be trapped in one channel. The WE approach is then used to sample the full transition path ensemble in implicitly solvated alanine dipeptide at two different temperatures. The ensembles are of sufficient quality to permit quantification of the fractional importance of each channel, even at T=300K when brute-force simulation is prohibitively expensive.

연구 동기 및 목표

• 분자의 구조적 전이 과정에서 다수의 반응 채널을 샘플링하는 데 도전하는 것.
• 가중 평균(WE) 방법이 TPS 및 유사한 경로 샘플링 방법에서 흔히 발생하는 갇힘 문제를 피할 수 있음을 입증하는 것.
• 브루트 포스 시뮬레이션으로는 비용이 지나치게 많이 들기 때문에 실현 가능하지 않은 실온 조건(300K)에서 통계적으로 엄밀한 고품질의 전이 경로 집합을 생성하는 것.
• 전이 경로 집합 내에서 각기 다른 반응 채널의 비중을 정량화하는 것.

제안 방법

• WE 방법은 브루트 포스 경로의 비편향 복제 및 조합을 통해 특정 채널에 대한 편향 없이 경로 공간을 탐색한다.
• 이 방법은 가중치를 동적으로 조정하고 경로를 분할하거나 융합함으로써 통계적 정확성을 유지하는 가중치가 부여된 경로 집합을 유지한다.
• WE와 TPS의 다수의 채널을 샘플링하는 행동을 대비하기 위해 이차원 모델 퍼텐셜에 적용한다.
• 알라닌 다펩타이드에 대해는 암묵적 용매 모델을 사용하여 500K와 300K에서 C7eq와 C7ax 상태 간의 전이를 샘플링한다.
• 다이헤드럴 각도 Φ와 Ψ의 회전 방향을 기반으로 하여 확장된 다이헤드럴 평면에서 경로 집합을 분석하고, 서로 다른 채널을 식별 및 분류한다.
• 브루트 포스 시뮬레이션 결과와 비교하여 반응 속도와 경로 분포를 분석함으로써 WE 결과의 타당성을 검증한다.

실험 결과

연구 질문

• RQ1 TPS와 달리 WE 방법은 단일 경로에 갇히지 않고 다수의 반응 채널을 샘플링할 수 있는가?
• RQ2 실제 생체분자 시스템인 알라닌 다펩타이드의 고온 및 실온 조건에서 WE 방법은 전이 경로 샘플링에 얼마나 효과적인가?
• RQ3 300K에서 각기 다른 반응 채널이 전체 전이 경로 집합에 기여하는 비중은 얼마인가?
• RQ4 브루트 포스 시뮬레이션이 불가능한 조건에서 WE는 정확한 반응 속도와 경로 집합을 동시에 제공할 수 있는가?

주요 결과

• 500K에서 WE 방법은 브루트 포스 결과인 1.5×10⁻¹/ns와 뛰어난 일치를 보이는 반응 속도 1.4×10⁻¹/ns를 도출하였다.
• 300K에서 WE 방법은 반응 속도 1.6×10⁻³/ns를 계산하여 500K에 비해 장벽 높이가 100배 증가했음을 시사하였다.
• 300K에서의 경로 집합은 Φ와 Ψ의 시계방향 및 반시계방향 회전에 해당하는 네 가지 주요 경로 유형을 드러내었으며, '아래 오른쪽' C7ax 상태가 가장 우세한 채널임을 확인하였다.
• '아래 오른쪽' 채널은 Ψ=0과 Φ=0 근처에 높은 장벽으로 인해 두 개의 서브채널로 나뉘어 있었다.
• 경로 분포 분석 결과, 저온 조건에서 주요 채널이 더욱 두드러지게 나타나는 것으로 나타나, 이 채널이 동역학적으로 유리한 경향을 보임을 확인하였다.
• 본 연구는 알라닌 다펩타이드에 대해 통계적으로 엄밀한 전이 경로 집합을 처음으로 보고하였으며, 채널의 비중을 정량화할 수 있게 하였다.