단백질 응집을 위한 통계역학 모델링

초록

본 논문은 단백질·펩타이드 자체조립을 통계역학적으로 기술한다. 2‑state·3‑state 포츠 모델을 기반으로 Zimm‑Bragg 형태의 전이 행렬을 정확히 해석해 구역(핵) 형성, 상도표, 시트 함량 등을 계산한다. Aβ(1‑40)와 Curli 섬유에 적용해 실험과 좋은 일치를 보였다.

상세 분석

이 연구는 단백질 응집을 미시적인 자유에너지 모델로 전환함으로써 복잡한 다단계 집합 현상을 단순화한다. 저자들은 단백질을 코일, α‑헬릭스, β‑시트의 세 가지 구조 상태로 구분하고, 각 상태 사이의 상호작용을 포츠 모델의 q‑state 스핀 변수(t_i)로 표현한다. 여기서 q=2 혹은 q=3을 선택해 2‑state(시트‑코일) 혹은 3‑state(헬릭스‑시트‑코일) 시스템을 구축한다. 핵심은 전이 행렬을 이용해 그랜드 파티션 함수를 정확히 계산할 수 있다는 점이다. 전이 행렬은 Zimm‑Bragg 모델의 전파 파라미터(s)와 개시 파라미터(σ)를 일반화한 형태로, P_1, P_2(구조별 결합 에너지), K(구조와 무관한 결합), R_j(구조 경계 에너지) 등 여섯 개의 물리적 파라미터만으로 전체 자유에너지 지형을 기술한다.

특히, 용매와의 인터페이스를 포함한 1‑D 격자 모델을 확장해 핵(n_c) 형성에 필요한 자유에너지 A를 도입함으로써 핵형성 장벽과 성장 단계의 비대칭성을 정량화한다. 이때 n_i(0/1) 변수는 격자점이 용매인지 단백질인지 구분하고, χ 함수(δ의 보완)를 이용해 용매‑단백질, 단백질‑단백질 경계 에너지를 체계적으로 포함한다. 결과적으로, 핵이 형성되면 K와 P_j에 의해 빠른 선형 성장(elongation)이 가능하고, R_j는 구조 전이(코일‑시트, 헬릭스‑시트 등)의 억제 혹은 촉진 역할을 한다.

모델은 n×N 스트립 격자(예: 2×N, 3×N 등)로 확장돼 다중 필라멘트 간 횡방향 상호작용을 기술한다. 이러한 n×N 시스템은 실제 아밀로이드 섬유가 보이는 다중 체인 배열을 근사화한다. 전이 행렬을 n‑차원으로 일반화함으로써, 필라멘트 간 결합 에너지와 구조적 정렬을 포함한 전체 자유에너지 함수를 정확히 구할 수 있다.

실험 적용에서는 Aβ(1‑40)와 Curli 섬유의 시트 함량, 핵 크기, 농도 의존적 상전이를 모델 파라미터에 맞추어 시뮬레이션하였다. 파라미터 추정은 실험적 CD 스펙트럼, 전자현미경 이미지, 용액 농도 데이터와 비교해 수행되었으며, 모델이 예측한 시트 함량과 핵 성장 속도가 관측값과 높은 일치를 보였다. 이는 포츠‑전이 행렬 접근법이 복잡한 단백질 응집 현상을 정량적으로 설명할 수 있음을 증명한다.

이 논문의 주요 기여는 (1) 단백질 구조 상태를 포츠 스핀으로 일반화해 다중 상태 모델을 구축, (2) 전이 행렬을 이용해 그랜드 파티션 함수를 정확히 해석, (3) 용매‑핵 인터페이스와 다중 필라멘트 상호작용을 포함해 실제 섬유 구조를 재현, (4) 실험 데이터와 정량적 일치를 통해 모델의 실용성을 검증한 점이다. 향후에는 비평형 동역학, 변형된 격자(곡률·분기) 및 약물 결합 효과 등을 포함해 모델을 확장할 여지가 크다.