WW 도메인 폴딩의 두 가지 지배 경로와 DRP 접근법

초록

본 연구는 원자 수준의 실험적 힘장을 이용해 WW 도메인 Fip35의 폴딩 메커니즘을 Dominant Reaction Pathways(DRP) 방법으로 분석하였다. 두 개의 주요 폴딩 경로가 존재함을 확인했으며, 각각은 두 개의 헤어핀 형성 순서가 서로 다르다. 온도와 초기 변성 상태에 따라 경로 비중이 변하고, 비원래 상호작용은 전이 구간에서 큰 영향을 미치지 않음이 밝혀졌다. DRP는 48 CPU 코어에서 1시간 내에 원자 수준의 폴딩 궤적을 생성할 수 있어 대규모 단백질 폴딩 연구에 실용적이다.

상세 분석

이 논문은 WW 도메인 Fip35의 폴딩을 원자 수준에서 정밀하게 기술하기 위해 AMBER ff99SB 힘장을 사용하고, 용매 효과는 Generalized Born 모델로 근사하였다. DRP 이론은 과잉 감쇠 Langevin 방정식에 기반해 온스아거-마촐프 액션을 최소화함으로써 가장 확률이 높은 전이 경로를 찾는다. 저자들은 전통적인 HJ 액션 직접 최소화가 10⁴~10⁵ 단계가 필요해 비현실적임을 지적하고, 대신 ratchet‑and‑pawl MD(rMD)로 다수의 시도 경로를 생성한 뒤, Eq. 3의 확률 가중치를 이용해 최적 경로를 선별하였다. rMD에서는 접촉 지도와 원자 간 거리 기반의 집합 변수 z를 정의하고, z가 감소할 때만 바이어스가 작용하도록 설계해 비물리적 힘을 최소화했다.

두 개의 지배 경로는 (1) 헤어핀 1이 먼저 완전하게 형성되고 이후 헤어핀 2가 접히는 경우와 (2) 헤어핀 2가 먼저 형성되는 경우이다. 온도 의존적인 자유에너지 계산을 두 개의 코스그레이드 모델(Go‑type 원래 상호작용 + 비원래 상호작용)에서 수행했으며, 고온에서는 경로 2의 비중이 증가함을 확인했다. 이는 비원래 상호작용이 전이 상태에서 크게 기여하지 않으며, 폴딩 경로의 형태는 주로 원래 상호작용에 의해 결정된다는 결론을 뒷받침한다.

또한 초기 변성 구조가 경로 선택에 영향을 미친다는 점을 강조한다. 여러 초기 구조를 무작위로 생성한 뒤 DRP 분석을 수행했을 때, 초기 구조가 헤어핀 1 쪽으로 이미 부분적으로 형성된 경우 경로 1이, 헤어핀 2 쪽으로 편향된 경우 경로 2가 우세하게 나타났다.

계산 효율성 측면에서, 48개의 CPU 코어를 이용해 단일 폴딩 궤적을 1시간 이내에 완성했으며, 이는 전통적인 전통 MD가 수개월~수년이 걸리는 상황과 비교해 획기적인 속도 향상이다. 따라서 DRP는 대규모 단백질 폴딩 네트워크를 탐색하는 데 실용적인 도구가 될 수 있다.