베타시트 단백질 접힘 전이상태와 Φ값 모델링

초록

본 연구는 작은 β‑시트 단백질인 PIN과 FBP WW 도메인의 전이상태를 두 개의 β‑헤어핀 중 하나가 형성된 두 개의 구조 혼합물로 가정하고, 변이로 인한 자유에너지 변화를 헤어핀 1, 헤어핀 2, 그리고 소수성 코어 세 부분으로 분해한다. 단일 파라미터(두 헤어핀의 전이상태 형성 정도)를 최적화함으로써 실험 Φ값을 정량적으로 재현하고, 특히 음의 Φ값을 자연스럽게 설명한다.

상세 분석

이 논문은 단일 장벽을 갖는 소형 단일 도메인 단백질의 접힘 메커니즘을 정량적으로 해석하려는 시도이다. 기존 Φ값 분석은 변이가 접힘 전이상태에 미치는 구조적 기여를 정성적으로만 제시했으며, 특히 음의 Φ값은 해석이 어려웠다. 저자들은 WW 도메인의 3‑스트랜드 β‑시트를 두 개의 β‑헤어핀(1,2)과 중앙의 소수성 코어로 분해하고, 전이상태를 “헤어핀 1만 형성된 상태”와 “헤어핀 2만 형성된 상태”라는 두 개의 미시적 구성으로 모델링한다. 이 가정은 최근 설계된 3‑스트랜드 β‑시트 단백질의 분자동역학 시뮬레이션에서 관찰된 전이상태와 일치한다.

변이에 의해 발생하는 자유에너지 변화 ΔΔG는 각 구조 요소별 기여(ΔΔG₁, ΔΔG₂, ΔΔG_core)로 분해된다. 전이상태의 전체 자유에너지 장벽은 두 구성의 가중 평균으로 표현되며, 가중치 p는 헤어핀 1이 전이상태에 차지하는 비율, (1‑p)는 헤어핀 2의 비율이다. Φ값은 실험적으로 측정된 접힘 속도 변화와 ΔΔG의 비율로 정의되는데, 모델식 Φ = p·(ΔΔG₁/ΔΔG) + (1‑p)·(ΔΔG₂/ΔΔG) + (ΔΔG_core/ΔΔG)·γ 로 전개된다. 여기서 γ는 코어가 전이상태에 어느 정도 형성되어 있는지를 나타내는 추가 파라미터이며, 본 연구에서는 γ를 0으로 두어 코어가 전이상태에 거의 기여하지 않음을 가정한다.

실험 데이터에 대해 p를 최적화한 결과, PIN 도메인에서는 p≈0.6, FBP WW 도메인에서는 p≈0.4가 최적값으로 도출되었다. 이는 두 단백질이 각각 헤어핀 1과 헤어핀 2가 전이상태에서 더 크게 기여한다는 것을 의미한다. 모델은 양의 Φ값뿐 아니라, 변이가 코어에 위치하면서 전이상태에서 구조가 약간 파괴되는 경우에 나타나는 음의 Φ값도 자연스럽게 재현한다. 이는 기존 모델이 전이상태를 단일 구조로 가정함으로써 놓쳤던 다중 경로성(multi‑pathway) 현상을 포착한 것이다.

또한, 변이 효과를 헤어핀별로 정량화함으로써 특정 잔기가 어느 헤어핀에 속하는지, 그리고 그 잔기의 친수성/소수성 특성이 전이상태 안정화에 어떻게 작용하는지를 명확히 알 수 있다. 예를 들어, 헤어핀 1에 위치한 친수성 잔기의 친수성 증가 변이는 전이상태에서 헤어핀 1이 부분적으로 형성된 상태를 방해해 Φ값을 감소시키는 반면, 같은 위치의 소수성 잔기 강화 변이는 Φ값을 상승시킨다.

이 모델의 강점은 파라미터 수를 최소화하면서도 실험 Φ값을 높은 정확도로 재현한다는 점이다. 단일 파라미터 p와 코어 기여 γ(=0)만으로 10여 개 이상의 변이에 대한 Φ값을 설명했으며, 이는 기존의 복잡한 구조 기반 시뮬레이션보다 효율적이다. 다만, 코어 기여를 0으로 가정한 것이 모든 β‑시트 단백질에 적용 가능한지는 추가 검증이 필요하다. 향후 다양한 β‑시트 구조와 더 큰 단백질에 대해 동일한 프레임워크를 적용하면 전이상태의 다중 경로성을 일반화할 수 있을 것으로 기대된다.