단백질 접힘을 위한 통계역학 모델

초록

본 논문은 물 환경에서 단백질의 접힘·전개 전이를 기술하기 위해 두 단계 접힘 동역학을 기반으로 한 새로운 통계역학 공식화를 제시한다. 친수성·소수성 상호작용을 묘사한 수소화 탄화수소 용매 모델을 이용해 무한히 희석된 단백질 용액의 분배함수를 유도하고, 이를 통해 온도에 따른 단백질의 열용량 변화를 계산한다. 계산 결과는 황색포도상구균 뉴클레아제와 메트미오글로빈의 실험 데이터와 좋은 일치를 보인다.

상세 분석

이 연구는 단백질 접힘을 열역학적으로 정량화하기 위해 두 단계(two‑stage) 접힘 메커니즘을 가정한다. 첫 번째 단계는 고유 구조가 형성되는 ‘초기 접힘’이며, 두 번째 단계는 이미 형성된 구조가 물과의 상호작용을 통해 최종적인 안정 상태에 도달하는 ‘정착 단계’로 정의된다. 이러한 가정 하에 저자들은 각 단계별 자유에너지 기여를 분리하고, 이를 통합한 전체 분배함수 Z를 구성한다. 핵심은 물 속에서 소수성 효과를 정량화하기 위해 기존의 ‘hydrophobic hydrocarbon solvation model’를 차용한 점이다. 이 모델은 소수성 잔기가 물에 노출될 때 발생하는 엔트로피 감소와 엔탈피 변화를 파라미터화하여, 단백질 전체의 용매화 자유에너지 ΔG_solv를 계산한다. 무한히 희석된 용액 가정은 단백질 간 상호작용을 무시함으로써, 단백질 자체와 용매 사이의 상호작용만을 순수하게 분석할 수 있게 한다.

분배함수로부터 열용량 C_p(T) = (∂⟨E⟩/∂T)_V 를 도출하고, 온도 구간별로 두 단계의 전이 확률을 Boltzmann 인자 exp(−ΔG_i/kT) 로 표현한다. 결과적으로 C_p(T)는 저온에서의 ‘접힘 전이’와 고온에서의 ‘전개 전이’ 두 개의 피크를 보이며, 이는 실험적으로 관찰되는 열용량 곡선과 일치한다. 특히, 스태필로코커스 뉴클레아제와 메트미오글로빈에 대해 모델 파라미터를 조정하지 않고도 정량적인 일치를 얻은 점은 모델의 보편성을 시사한다.

이 논문의 주요 기여는 (1) 두 단계 접힘 메커니즘을 통계역학적으로 정식화한 점, (2) 소수성 용매 모델을 단백질 전체에 적용해 용매화 자유에너지를 정확히 계산한 점, (3) 실험 데이터와의 직접 비교를 통해 모델 검증을 수행한 점이다. 다만, 무한 희석 가정과 소수성 모델의 파라미터가 특정 단백질에 최적화되어 있지 않다는 점은 향후 연구에서 다중 단백질 시스템이나 고농도 용액에 대한 확장이 필요함을 암시한다.