배제 부피와 국부 구조 협동성 그리고 단백질 접힘 속도의 고분자 물리학

초록

이 논문은 두 단계(two‑state) 접힘을 보이는 다양한 단백질을 대상으로, 배제 부피 효과와 미세 구조 협동성을 포함한 거시적 변분 모델을 이용해 접힘 장벽과 속도를 계산한다. 배제 부피 항을 넣으면 접힘 핵의 경계가 뚜렷해지고 자유에너지 장벽이 증가한다는 것이 확인되었으며, 장벽 높이는 접힘 네이티브 토폴로지를 나타내는 접촉 순서(contact order)와 강하게 상관된다. 또한, 전이 상태 집합의 강성 변화는 단백질마다 비균일하게 나타나지만, 프리팩터 k₀는 전체 집합에서 5배 이내의 변동만 보이며 평균 4 µs⁻¹ 정도이다. 모델은 안정도·점도 의존성, 다운힐 접힘 등 세부 현상까지 확장 가능함을 제시한다.

상세 분석

본 연구는 고분자 물리학적 접근을 통해 단백질 접힘 장벽을 정량화하려는 시도로, 기존의 단순한 Go‑model에 배제 부피(Excluded Volume)와 국부 구조 협동성(Cooperativity)을 명시적으로 포함한 변분(coarse‑grained variational) 모델을 구축하였다. 배제 부피 항은 비공유적 상호작용을 모사함으로써, 접힘 핵 내부와 주변 비접힌 영역 사이에 물리적 경계를 형성하게 만든다. 이 경계가 뚜렷해질수록 전이 상태(Transition State Ensemble, TSE)의 구조적 이질성이 증가하고, 자유에너지 장벽(ΔG‡)이 상승한다는 결과가 도출되었다. 특히, 장벽 높이와 네이티브 토폴로지를 나타내는 접촉 순서(Contact Order, CO) 사이의 상관계수는 0.85 이상으로, 토폴로지적 복잡성이 장벽을 결정짓는 주요 요인임을 재확인한다.

프리팩터 k₀는 고분자 사슬의 내부 마찰과 모노머 이완 시간(τₘ≈30 ns)을 기반으로 직접 계산되었으며, 전체 단백질 집합에 대해 1/k₀가 1–5 µs 범위, 평균 4 µs로 비교적 균일함을 보였다. 이는 “프리팩터는 모든 두 단계 단백질에 대해 거의 동일하다”는 기존 가정에 이론적 근거를 제공한다. 또한, 핵의 강성을 증가시키면 전이 상태의 국부 구조가 비균일하게 변형되는데, 이는 특정 단백질에서만 관찰되는 현상으로, 구조적 협동성의 정도가 단백질 고유의 토폴로지와 상호작용한다는 점을 시사한다.

모델은 또한 안정도(ΔGₙ)와 점도(η)의 변화를 통해 접힘 속도의 비선형 의존성을 예측할 수 있다. 예를 들어, ΔGₙ이 증가하면 장벽이 낮아져 속도가 급격히 상승하지만, 점도가 크게 증가하면 내부 마찰이 지배적으로 작용해 전체 속도가 감소한다. 이러한 예측은 실험적 점도 의존성 데이터와도 일치한다. 마지막으로, 장벽이 거의 사라지는 경우(ΔG‡≈k_BT 이하)에는 전통적인 장벽 통과 모델이 붕괴하고, 다운힐(folding downhill) 현상이 나타나는 영역을 모델이 자연스럽게 포착한다는 점에서, 고전적인 두 단계 모델을 넘어서는 확장성을 보여준다.

전반적으로, 배제 부피와 협동성을 포함한 변분 모델은 단백질 접힘 장벽과 속도를 토폴로지와 동역학적 파라미터에 기반해 정확히 예측하며, 프리팩터의 보편성, 전이 상태의 구조적 이질성, 그리고 환경 변수(안정도·점도)의 영향을 통합적으로 설명한다. 이는 향후 단백질 설계와 질병 관련 변이 분석에 있어, 물리적 기반의 정량적 예측 도구로 활용될 수 있는 중요한 기반을 제공한다.