미세정준 분석으로 밝힌 헬리컬 단백질 접힘 협동성

초록

본 연구는 중간 해상도의 암시적 용매 모델을 이용한 코스그레인 시뮬레이션과 마이크로캐노니컬 분석을 결합해, 짧은 α‑헬릭스는 두 단계(two‑state) 협동성을 보이는 반면 사슬 길이가 늘어나면 다중 핵생성 부위가 형성돼 중간 상태가 풍부해지는 다운힐 폴더가 된다. 반면, 세 개의 헬릭스가 서로 상호작용하는 de novo 3‑헬릭스 번들은 여전히 두 단계 협동성을 유지한다.

상세 분석

이 논문은 전통적인 열역학적 지표(예: 열용량 곡선)만으로는 두 단계 협동성의 존재 여부를 확정하기 어려운 점을 지적하고, 밀도 상태(Ω(E))를 직접 추정해 엔트로피 S(E)=k_B ln Ω(E)를 마이크로캐노니컬 관점에서 분석하는 방법을 제시한다. 저자들은 replica‑exchange MD와 WHAM을 이용해 충분히 샘플링된 p(E)를 얻고, 이를 통해 ΔS(E)와 역마이크로캐노니컬 온도 T_μ⁻¹(E)=∂S/∂E를 계산하였다. ΔS(E)에서 볼 수 있는 ‘볼록 인트루더(convex intruder)’는 자유 에너지 장벽을 의미하며, 이는 비제로 마이크로캐노니컬 잠열 ΔQ와 역온도 곡선의 ‘백워딩(back‑bending)’ 현상으로 나타난다.

짧은 사슬 (AAQAA)_3은 ΔS(E)에서 뚜렷한 장벽과 ΔQ>0을 보이며, 헬리시티 θ(E)와 헬릭스 수 H(E)도 급격히 전이한다. 이는 전형적인 두 단계 폴더의 특징이다. 사슬 길이가 7로 늘어나면 장벽이 감소하지만 여전히 비제로 ΔQ가 존재한다. 그러나 n=10,15에서는 ΔS(E)의 볼록 인트루더가 사라지고 T_μ⁻¹(E)가 단조적으로 변하면서 ΔQ=0이 된다. 구조적으로는 전이 구간에서 H(E) > 1인 다중 헬릭스가 형성되며, 이는 여러 핵생성 부위가 동시에 존재함을 의미한다. 즉, 사슬이 길어질수록 ‘핵생성 경쟁’이 일어나 중간 상태가 풍부해져 다운힐 폴딩으로 전이한다.

세 개의 헬릭스가 서로 얽힌 de novo 3‑헬릭스 번들(α3D)은 전체 길이는 n=15와 비슷하지만, 헬리시티와 H(E)의 변화가 단조적이며 ΔS(E)에 여전히 볼록 인트루더와 비제로 ΔQ가 존재한다. 이는 헬릭스‑헬릭스 상호작용이 전이 장벽을 강화해 두 단계 협동성을 유지한다는 이론적 예측과 일치한다.

2차 구조 분석에서는 에너지와 잔기 위치에 따른 헬리시티 분포를 제시했으며, 짧은 사슬은 중심에서 핵생성이 시작돼 양쪽 말단으로 전파되는 반면, 긴 사슬은 중간에 두 개의 독립적인 헬릭스가 형성돼 각각이 독립적으로 녹는다. 3차 구조 측면에서는 반경(R_g)와 정규화 원통도(c) 변화를 조사했는데, 긴 사슬과 α3D 모두 전이 전후에 비단조적 변화를 보여 ‘체인 콜랩스’가 먼저 일어나고 이후 헬리컬 구조가 형성되는 순서를 확인했다. 또한 비국소 접촉 수와 네이티브 접촉 수를 비교했을 때, 전이 구간에서 비국소 접촉이 급증했다가 전이 후 감소하는 비대칭적인 패턴이 관찰돼, 네이티브 구조로 가기 위해 비국소 접촉을 깨뜨려야 함을 시사한다.

전반적으로 마이크로캐노니컬 분석은 엔트로피 곡선의 미세한 볼록성까지 포착함으로써, 전통적인 열역학적 지표가 놓칠 수 있는 ‘중간 상태의 풍부함’과 ‘핵생성 부위의 다중성’을 명확히 드러낸다. 이는 실험적 열용량 측정이나 CD 스펙트럼만으로는 구분하기 어려운 다운힐 폴더와 두 단계 폴더의 경계를 정량적으로 정의하는 데 큰 의미가 있다.