단백질 G 접힘 궤적의 원자 수준 이정표

초록

본 연구는 래칫‑앤‑폴리 원리를 이용한 바이어스 알고리즘을 적용해, 고성능 컴퓨팅 없이도 명시적 용매 환경에서 단백질 G B1 도메인의 8가지 접힘 궤적을 원자 수준으로 재현하였다. 시뮬레이션 결과, 변성 상태에서 제한된 수의 핵심 네이티브 및 비네이티브 접촉이 다수의 다른 접촉을 유도하는 계층적 네트워크를 형성함을 보여준다. 전이 상태 집합은 실험 φ‑값과 0.69의 상관계수를 가지며, 구조적으로는 토폴로지는 원형이지만 측쇄와 수소결합이 완전히 형성되지 않은 상태이다.

상세 분석

이 논문은 전통적인 분자동역학(MD) 시뮬레이션이 요구하는 막대한 계산 비용을 극복하기 위해 ‘래칫‑앤‑폴리(Ratchet‑and‑Pawl)’ 원리를 기반으로 한 바이어스 알고리즘을 도입하였다. 이 방법은 시스템이 이미 형성된 접촉을 유지하도록 ‘래칫’ 역할을 하면서, 새로운 접촉이 형성될 때만 자유롭게 움직일 수 있게 함으로써, 실제 물리적 경로를 크게 왜곡하지 않으면서도 접힘 과정을 가속화한다. 결과적으로, 1.2 Å에서 2.5 Å 사이의 RMSD를 보이는 8개의 독립적인 접힘 궤적을 얻었으며, 이는 고성능 슈퍼컴퓨터 없이도 원자 수준의 정확도를 확보할 수 있음을 증명한다.

시뮬레이션에서 관찰된 가장 중요한 현상은 변성 상태에서 형성되는 접촉들의 인과관계 네트워크이다. 전체 접촉 중 소수의 로컬 및 비국소 네이티브 접촉이 ‘원인’ 역할을 하여 다수의 ‘결과’ 접촉을 유도한다는 점에서, 전통적인 ‘핵심 접촉’ 모델을 뒷받침한다. 특히, 실험적으로 약한 변성 조건에서 관찰된 NOE 신호와 일치하는 몇몇 핵심 접촉이 초기 단계에서 형성되는 것이 확인되었다. 흥미롭게도, 비네이티브 접촉도 일시적으로 형성되어 접힘 경로를 조정하는 데 기여한다는 점이 강조된다. 이는 비네이티브 상호작용이 단순히 ‘오류’가 아니라, 폴딩 라인드를 탐색하는 데 유용한 중간 상태를 제공한다는 새로운 관점을 제시한다.

전이 상태(Transition State, TS) 집합에 대한 φ‑값 분석에서는 실험값과 0.69의 높은 상관계수를 기록하였다. 구조적으로 TS는 전반적으로 토폴로지적으로 원형(native‑like)이며, 전반적인 2차 구조가 이미 형성돼 있다. 그러나 측쇄(side‑chain) 배치와 수소결합 네트워크는 아직 미완성 상태이며, 이는 실제 실험에서 관찰되는 ‘부분적으로 형성된’ 전이 상태와 일치한다. 이러한 결과는 TS가 ‘구조적 균일성’보다는 ‘에너지 장벽을 넘는 동역학적 특성’에 의해 정의된다는 기존 가설을 뒷받침한다.

전반적으로, 이 연구는 고전적인 MD 시뮬레이션이 직면한 계산적 한계를 뛰어넘어, 명시적 용매 환경에서 원자 수준의 폴딩 메커니즘을 정밀하게 재현할 수 있음을 보여준다. 또한, 네트워크 기반의 인과관계 분석과 비네이티브 접촉의 역할 규명은 단백질 폴딩 이론에 새로운 통찰을 제공한다.