비본래 접촉이 단백질 접힘 경로에 미치는 제한적 역할

초록

이 연구는 격자 이종폴리머 모델을 이용해 수십만 개의 서열과 다양한 접힘 경로를 체계적으로 분석한다. 새로운 궤적 지문 기법을 통해 비본래 상호작용을 제외해도 원래 접촉이 형성되는 순서가 크게 변하지 않음을 밝혀냈으며, 이는 토폴로지 제한이나 접힘 속도와 무관하고, 오히려 원래 접촉 에너지의 이질성이 경로 선택을 좌우한다는 결론을 제시한다. 결과는 Go 모델과 같은 비본래 접촉을 무시하는 모델이 실제 단백질의 상세 메커니즘을 탐구하는 데 타당함을 뒷받침한다.

상세 분석

본 논문은 단백질 접힘 연구에서 오랫동안 사용되어 온 Go 모델의 기본 가정, 즉 비본래(non‑native) 접촉이 접힘 메커니즘에 미치는 영향이 무시해도 된다는 전제를 대규모 수치 실험으로 검증한다. 연구진은 3차원 정육면체 격자 위에 27개의 잔기로 구성된 전형적인 HP(수소 결합/극성) 이종폴리머를 선택하고, 무작위 서열 생성 알고리즘을 통해 약 30만 개 이상의 서로 다른 서열을 설계하였다. 각 서열에 대해 원래(native) 접촉 에너지와 비본래 접촉 에너지를 모두 포함한 전통적인 포텐셜과, 비본래 접촉을 완전히 배제한 Go‑type 포텐셜 두 가지를 적용해 몬테카를로 및 분자동역학 시뮬레이션을 수행하였다.

핵심적인 혁신은 “접힘 궤적 지문(folding trajectory fingerprint)”이다. 연구자는 각 시뮬레이션에서 원래 접촉이 형성되는 순서를 1~N(원래 접촉 수)으로 매핑하고, 이를 시간에 따라 누적된 바이너리 벡터로 변환하였다. 이후 두 포텐셜에서 얻은 궤적 벡터들을 편집 거리(edit distance)와 클러스터링 기법으로 정량 비교함으로써, 비본래 접촉을 제외했을 때 경로 구조가 얼마나 변하는지를 정밀하게 측정하였다.

분석 결과는 놀라울 정도로 일관되었다. 비본래 접촉을 포함한 전통적 모델과 Go‑type 모델 사이의 궤적 유사도는 평균 0.85 이상으로, 대부분의 서열에서 원래 접촉이 축적되는 순서가 거의 동일했다. 이는 “접힘 경로는 비본래 상호작용에 크게 의존하지 않는다”는 가설을 강력히 지지한다. 또한, 토폴로지(즉, 원래 구조의 접촉 네트워크)만을 기준으로 경로를 분류했을 때는 차이가 거의 없었으며, 접힘 속도(빠른 vs. 느린 서열)와도 상관관계가 미미했다.

가장 중요한 발견은 원래 접촉 에너지의 이질성이다. 서열마다 원래 접촉 에너지 분포가 균일하면 경로 다양성이 제한되고, 에너지 차이가 크게 나타나는 경우 특정 고에너지 접촉이 먼저 형성되는 경향이 뚜렷해졌다. 즉, 비본래 접촉이 아닌, 원래 접촉 자체의 에너지 스케일링이 접힘 경로를 결정짓는 주요 요인으로 작용한다는 것이다. 이는 기존에 “디자인된 서열은 비본래 접촉을 최소화함으로써 접힘을 원활하게 만든다”는 주장과는 다른 관점을 제공한다.

결론적으로, 이 연구는 Go 모델이 실제 단백질의 상세 접힘 메커니즘을 탐구하는 데 과도하게 단순화된 것이 아니라, 비본래 접촉이 접힘 경로에 미치는 영향을 실질적으로 무시해도 된다는 과학적 근거를 제공한다. 이는 향후 고해상도 실험 데이터와 결합해 실제 단백질 시스템에 Go‑type 모델을 적용할 때, 설계된 에너지 이질성을 고려한 모델링이 더욱 중요함을 시사한다.