이차·삼차 구조 형성의 상호작용이 단백질 접힘 협동성에 미치는 영향

초록

본 연구는 고해상도 코스모-그레인 모델과 레플리카 교환 분자동역학을 이용해 짧은 α‑헬릭스, 긴 α‑헬릭스, 그리고 3‑헬릭스 번들을 조사하였다. 마이크로캐노니컬 분석을 통해 에너지 밀도 상태를 정확히 구하고, 엔트로피 곡선의 볼록성(ΔS)과 역온도 곡선, 구형 반경(Rg) 등을 에너지 기준으로 분석했다. 짧은 헬릭스는 두 상태 전이(두 단계) 특성을 보였고, 긴 헬릭스는 다운힐(연속) 전이, 3‑헬릭스 번들은 다시 두 상태 전이를 나타냈다. 핵심 결과는 이차 구조 형성(수소 결합)과 비원주성 삼차 접촉 손실이 전이의 성격을 결정한다는 점이다.

상세 분석

이 논문은 단백질 접힘 협동성을 미시적 열역학 관점에서 재조명한다. 전통적으로는 열용량 곡선의 칼로리메트릭 지표(δ)로 두 상태와 연속 전이를 구분했지만, 저자들은 마이크로캐노니컬 접근을 통해 엔트로피 S(E)=k_B ln Ω(E)와 그 1차 도함수인 역온도 β(E)=∂S/∂E를 직접 계산한다. ΔS(E)=H(E)−S(E) (H는 전이 구간에서의 접선) 를 이용해 엔트로피의 볼록성을 정량화함으로써, 중간 에너지 상태가 고갈되는 두 상태 전이에서는 ΔS가 뚜렷한 피크를, 연속 전이에서는 평탄한 곡선을 보인다.

시뮬레이션은 고해상도 코스모-그레인 모델(베라우·데스르노 2009)을 사용했으며, 레플리카 교환 MD와 WHAM(Weighted Histogram Analysis Method)으로 정확한 밀도 상태 Ω(E)를 얻었다. 구조적 지표로는 구형 반경 R_g(E)와 두 종류의 에너지 기여, 즉 수소 결합 에너지 E_hb(이차 구조)와 측쇄 상호작용 에너지 E_sc(삼차 구조)를 도입하고, 각각의 미분 dE_hb/dE, dE_sc/dE 를 통해 에너지 구간별 구조 형성 속도를 파악했다.

짧은 (AAQAA)_3 펩타이드는 ΔS에서 뚜렷한 피크와 비제로 잠열(ΔQ)을 보여 전형적인 두 상태 전이를 나타낸다. 이 구간에서 dE_hb/dE가 급격히 상승하고 dE_sc/dE는 거의 변하지 않아, 전이가 주로 이차 구조 형성에 의해 주도됨을 알 수 있다. 반면 (AAQAA)_15는 ΔS가 연속적이며 잠열이 0인 다운힐 전이를 보인다. 여기서는 R_g가 전이점 위에서 최소가 되며, 비원주성 비네이티브(compact non‑native) 상태에서 헬릭스로 재배열된다. dE_sc/dE가 음의 값을 갖는 구간이 존재하는데, 이는 삼차 접촉이 형성될 때 에너지 페널티가 발생함을 의미한다. 수소 결합 형성은 넓은 에너지 구간에 걸쳐 서서히 진행된다.

세 번째 시스템인 73‑아미노산 3‑헬릭스 번들(α3D, PDB 2A3D)은 ΔS에서 뚜렷한 피크와 비제로 잠열을 보여 두 상태 전이를 보이며, 전이 구간에서 R_g가 최소가 된다. 여기서는 dE_hb/dE와 dE_sc/dE 모두 전이 구간에 집중되어, 이차·삼차 구조가 동시에 급격히 변한다. 특히 측쇄‑측쇄 상호작용이 강해 비원주성 삼차 접촉이 크게 감소하면서 전이가 촉진된다.

이러한 결과를 종합하면, 전이의 협동성은 비원주성 삼차 접촉의 소실(dE_sc/dE < 0)과 이차 구조 형성(dE_hb/dE > 0)의 시점과 폭이 얼마나 겹치는가에 달려 있다. 두 상태 전이에서는 두 현상이 거의 동일한 에너지 구간에 집중되어 강한 협동성을 보이며, 다운힐 전이에서는 이차 구조 형성이 넓은 구간에 퍼져 있어 협동성이 약해진다. 이는 기존 격자 모델 및 이론적 헬릭스‑코일 전이 모델(예: Zimm‑Bragg)과 일치하면서도, 실제 연속적인 에너지-구조 관계를 정량적으로 제시한다는 점에서 의미가 크다.