베타 헤어핀 접힘 경로의 지배적 메커니즘

초록

본 연구는 Dominant Reaction Pathway(DRP) 방법을 이용해 베타‑헤어핀의 접힘 과정을 모델링하고, 네이티브와 비네이티브 상호작용을 모두 포함한 포텐셜을 적용하였다. DRP로 도출된 가장 확률이 높은 접힘 경로를 장시간 분자동역학(MD) 시뮬레이션 결과와 비교한 결과, 두 접근법이 일관된 경로를 제시함을 확인하였다. 비네이티브 소수성 상호작용이 접힘 동역학에 미치는 영향은 미미한 수준으로 나타났다.

상세 분석

본 논문은 베타‑헤어핀이라는 작은 단백질 2차 구조 요소의 접힘 메커니즘을 정량적으로 규명하기 위해 Dominant Reaction Pathway(DRP) 프레임워크를 적용하였다. DRP는 전이 상태 이론과 최적 경로 탐색을 결합한 방법으로, 고차원 자유에너지 지형에서 가장 확률이 높은 전이 경로를 효율적으로 추출한다. 연구진은 먼저 네이티브 접촉을 강화하는 Gō‑type 포텐셜과, 비네이티브 소수성 상호작용을 추가한 혼합 포텐셜을 정의하였다. 이 포텐셜은 실제 단백질이 겪는 네이티브와 비네이티브 상호작용을 동시에 고려함으로써, 기존 Gō‑모델이 과도하게 네이티브 구조에 편향되는 문제를 보완한다.

DRP 계산은 초기 비접힌 상태와 최종 네이티브 상태 사이의 경로를 수천 개의 이산화된 이미지로 표현하고, 각 이미지 사이의 전이 확률을 라그랑지안 최적화로 최소화한다. 이 과정에서 온도, 마찰 계수, 그리고 자유에너지 구배가 명시적으로 포함되어, 실제 물리적 환경을 반영한다. 최적화된 경로는 “지배적 접힘 경로”라 불리며, 이는 전통적인 MD 시뮬레이션에서 관찰되는 다수의 전이 사건과 통계적으로 일치한다.

MD 시뮬레이션은 1 µs 규모의 장시간 시뮬레이션을 수행해 수천 번의 접힘·펴짐 이벤트를 수집하였다. 이 데이터와 DRP에서 도출된 경로를 비교한 결과, 접힘 초기 단계에서 루프 형성 및 수소 결합 네트워크가 급격히 성장하는 양상이 동일하게 나타났다. 또한, 비네이티브 소수성 상호작용을 포함했을 때 경로의 전반적인 형태와 전이 장벽 높이가 크게 변하지 않았으며, 이는 비네이티브 상호작용이 베타‑헤어핀의 접힘 속도와 메커니즘에 미치는 영향이 제한적임을 시사한다.

이러한 결과는 DRP가 고차원 에너지 지형을 효율적으로 탐색하면서도, MD와 동일한 물리적 현상을 포착할 수 있음을 입증한다. 특히, DRP는 수백 마이크로초 수준의 장시간 동역학을 직접 시뮬레이션하지 않고도, 주요 전이 경로와 장벽 정보를 제공함으로써 계산 비용을 크게 절감한다. 비네이티브 상호작용에 대한 정량적 평가는, 작은 단백질 구조에서 네이티브 접촉이 주도적인 역할을 한다는 기존 가설을 강화한다.