하향형과 이중상태 단백질 접힘의 통계역학적 비교

초록

저자들은 간단한 통계역학 모델을 이용해 1BBL과 PIN1 WW 도메인의 접힘을 비교한다. 1BBL은 저온에서 에너지 장벽이 사라지는 ‘다운힐’ 접힘 특성을 보이지만, 일반적인 온도에서는 이중상태 거동과 유사한 거동을 나타낸다. 반면 PIN1은 전형적인 두 상태 폴더로, 온도 변화에 따라 명확한 자유에너지 장벽이 유지된다.

상세 분석

본 연구는 Wako‑Saitô‑Muñoz‑Eaton (WSME) 모델을 기반으로 한 통계역학적 프레임워크를 적용한다. WSME 모델은 각 아미노산 잔기의 네이티브 접촉을 이진 변수로 표현하고, 인접 잔기의 상호작용을 통해 전체 폴딩 자유에너지 지형을 계산한다. 저자들은 이 모델을 정확히 해석 가능한 형태로 단순화함으로써 평형 상태의 자유에너지 곡선과 열용량(Cp) 프로파일을 직접 구하였다. 1BBL과 PIN1의 구조적 파라미터(잔기 수, 네이티브 접촉 수, 접촉 강도 등)를 동일한 방식으로 입력하고, 온도 의존성을 조사하였다.

평형 분석 결과, 1BBL은 고온에서는 두 개의 뚜렷한 최소값을 갖는 이중상태 형태를 보이나, 온도가 낮아질수록 자유에너지 장벽이 급격히 감소한다. 특히 약 0.8 Tm 이하에서는 장벽이 사라져 연속적인 ‘다운힐’ 접힘이 가능함을 확인하였다. 반면 PIN1은 온도 구간 전반에 걸쳐 두 최소값 사이에 일정한 장벽을 유지하며, 전형적인 두 상태 전이를 나타낸다. 열용량 곡선에서도 차이가 뚜렷하게 드러난다. 1BBL은 낮은 온도에서 단일 피크가 약해지는 반면, PIN1은 온도 변화에 따라 뚜렷한 피크가 유지된다.

동역학적 측면에서는 마스터 방정식을 반정밀하게 풀어 접힘/펼침 속도 상수를 도출하였다. 저자들은 전이 경로를 두 단계(네이티브와 비네이티브)로 제한하고, 전이 행렬의 고유값을 분석함으로써 시간 상수의 온도 의존성을 평가했다. 1BBL은 고온에서 두 단계 지수함수적 감소를 보이지만, 저온에서는 단일 지수함수 형태로 전이 속도가 급격히 빨라진다. 이는 에너지 장벽이 사라지는 ‘다운힐’ 현상의 동역학적 증거이다. 반면 PIN1은 전 온도 구간에서 두 단계 지수함수 모델이 적합하며, 장벽에 의한 속도 제한이 지속된다.

이러한 결과는 ‘다운힐’ 폴딩이 반드시 온도에 독립적인 현상이 아니라, 특정 온도 이하에서만 실현될 수 있음을 시사한다. 또한, 모델 파라미터(특히 접촉 강도와 상호작용 범위)의 미세 조정이 폴딩 경로와 장벽 형성에 큰 영향을 미친다. 저자들은 실험적 열역학 데이터와 비교하여 모델의 정량적 정확성을 검증했으며, 1BBL의 경우 실험적으로 관찰된 비정상적인 열용량 곡선과 일치함을 확인하였다.

결론적으로, 통계역학적 모델은 다운힐 폴딩과 이중상태 폴딩 사이의 미세한 차이를 정량적으로 구분할 수 있는 강력한 도구이며, 온도 의존적인 자유에너지 지형과 동역학적 전이 특성을 동시에 설명한다는 점에서 의미가 크다.