세 나선 번들 단백질의 접힘 메커니즘: 보편성과 다양성

초록

본 연구는 전 원자 수준의 전이가능 전위 에너지를 이용해 단백질 A의 B 도메인과 빌린 헤드피스의 접힘 과정을 대규모 Monte Carlo 시뮬레이션으로 조사하였다. 클러스터링·그래프 이론 기반 분석과 p‑fold 검증을 통해 두 단백질 모두 접힘 초기에 2‑헬릭스 헤어핀 구조가 부분적으로 안정화되며, 전이 상태 집합(TSE) 전후의 구조적 차이를 원자 수준에서 규명하였다. 실험 Fi값과의 정량적 일치가 확인된 단백질 A와, 향후 실험 검증을 위한 Villin의 Fi값 예측을 제시한다.

상세 분석

본 논문은 전이가능(all‑atom transferable) 포텐셜을 기반으로 한 ‘ab initio’ Monte Carlo 접힘 시뮬레이션을 수천 회 수행함으로써, 3‑헬릭스 번들 단백질인 Protein A B‑도메인과 Villin headpiece의 접힘 경로를 원자 수준에서 정밀히 탐색하였다. 시뮬레이션 결과는 두 단백질 모두 초기 비정질 상태에서 2‑헬릭스 헤어핀(두 개의 연속된 α‑헬릭스와 그 사이의 루프) 형태가 빠르게 형성된 뒤, 이 구조가 부분적인 에너지 최소화 역할을 수행하며 전체 접힘 과정의 ‘프리‑핵심’ 단계가 된다. 이후 남은 헬릭스가 순차적으로 정렬되고, 최종적으로는 전통적인 3‑헬릭스 번들 형태가 완성된다.

전이 상태 집합(TSE)을 정의하기 위해 p‑fold(접힘 확률) 분석을 적용했으며, p‑fold≈0.5인 구조들을 대표적인 전이 상태로 선정하였다. 이들 TSE는 2‑헬릭스 헤어핀 구조가 유지된 채, 나머지 헬릭스가 부분적으로 형성된 ‘부분‑접힌’ 형태를 보인다. 특히, 두 단백질 모두 TSE에서 핵심적인 접촉 네트워크는 헤어핀의 말단 잔기와 나머지 헬릭스 사이의 장거리 상호작용으로 구성되어, 이는 전이 장벽을 결정짓는 주요 요인으로 작용한다는 점을 확인했다.

클러스터링은 그래프 이론 기반의 커뮤니티 탐지 알고리즘을 이용해 수행했으며, 각 클러스터는 에너지와 구조적 RMSD 기준으로 구분된다. 결과적으로, 초기 비정질 클러스터 → 2‑헬릭스 헤어핀 클러스터 → 전이 상태 클러스터 → 최종 네이티브 클러스터 순으로 전이 경로가 명확히 구분되었다. 이러한 단계적 구분은 실험적 φ‑값(Fi)과도 높은 상관관계를 보였으며, Protein A에 대해서는 시뮬레이션 Fi값이 실험값과 평균 절대 오차 <0.07로 일치하였다. Villin에 대해서는 아직 실험 데이터가 부족하지만, 시뮬레이션을 통해 도출한 Fi값을 향후 φ‑분석 실험에 대한 예측값으로 제시하였다.

또한, 본 연구는 이전에 수행된 Engrailed Homeodomain 연구와 비교했을 때, 3‑헬릭스 번들 단백질군 전반에 걸쳐 ‘2‑헬릭스 헤어핀 → 전이 상태 → 완전 번들’이라는 보편적인 접힘 메커니즘이 존재함을 제시한다. 이는 단백질 설계 및 변이체 분석에 있어, 초기 헤어핀 형성을 촉진하거나 방해하는 변이가 전체 접힘 속도와 안정성에 미치는 영향을 예측하는 데 유용한 통찰을 제공한다.