마이오트로핀 앙키린 반복 단백질의 균형과 동역학 정밀 분석

초록

본 연구는 Wako‑Saito‑Munoz‑Eaton(WSME) 모델을 이용해 작은 앙키린 반복 단백질인 마이오트로핀의 접힘 평형과 동역학을 정량적으로 조사한다. 모델은 이진 변수 기반의 네이티브‑센트릭 접근법으로, 기존 Ising 모델보다 미세한 구조 정보를 제공하면서도 정확히 해석 가능하다. 결과는 실험과 부분적으로 일치하여, 야생형 단백질이 약한 3상 평형을 보이면서도 비선형 자유에너지 지형에도 불구하고 2상‑유사 동역학을 나타냄을 보여준다. 돌연변이는 중간체를 안정화하거나 속도 제한 단계를 이동시켜 이러한 특성을 변화시킬 수 있다. 또한 N‑말단과 C‑말단을 통한 두 개의 대안 경로가 존재함을 확인했으며, 접힘과 풀림 경로가 서로 다르지만 밀접히 연관된다는 점을 제시한다.

상세 분석

WSME 모델은 각 아미노산 잔기의 접힌/펴진 상태를 0‑1 변수로 표현하고, 네이티브 접촉 에너지를 기반으로 전체 자유에너지 함수를 구성한다. 이 접근법은 Ising 모델이 단순히 반복 단위 간 상호작용만을 고려하는 반면, 개별 잔기의 구조적 기여와 비국소적인 상호작용을 동시에 반영한다는 점에서 우수하다. 마이오트로핀은 4개의 앙키린 반복으로 이루어진 작은 단백질이며, 실험적으로는 두 단계의 전이(중간체 존재)와 두 개의 병렬 경로(N‑말단, C‑말단)가 보고되었다. 모델 계산 결과는 자유에너지 프로파일이 다중 장벽을 가지고 있음에도 불구하고, 전체 접힘/펼침 반응은 실질적으로 두 단계(접힘과 풀림)로 지배된다는 것을 보여준다. 이는 자유에너지 지형이 ‘디자인’된 형태, 즉 중간 장벽이 충분히 낮아 반응 속도에 크게 기여하지 않는 경우에 해당한다. 돌연변이 시뮬레이션에서는 특정 반복을 안정화하거나 불안정화함으로써 중간체의 자유에너지 최소값을 상승·하강시켰으며, 결과적으로 속도 제한 단계가 N‑말단에서 C‑말단으로, 혹은 그 반대로 이동하는 현상이 관찰되었다. 이러한 변이는 실험에서 보고된 돌연변이 효과와 방향성은 일치하지만, 정량적 변화량은 모델이 과소평가하거나 과대평가하는 경향을 보였다. 또한, 접힘 과정에서는 C‑말단이 먼저 구조를 형성하고 N‑말단이 뒤따르는 경향이 있었으며, 풀림 과정에서는 그 순서가 역전되는 비대칭성이 드러났다. 이는 전통적인 ‘두‑상’ 해석이 실제 복잡한 경로를 단순화한 것이며, WSME 모델이 제공하는 미세한 동역학 정보를 통해 보다 정교한 메커니즘을 제시한다는 점에서 의미가 크다.