시토크롬 c의 동위원소 표지 아마이드 I 모드 진동 에너지 이완 이론적 연구

초록

본 연구는 물에 용해된 시토크롬 c에서 동위원소 표지된 아마이드 I 진동 모드의 에너지 이완(VER)을 시간 의존 섭동 이론으로 분석한다. 단백질과 물 각각이 기여하는 VER 경로를 분리하고, 방사형·각형 흥분 함수를 이용해 공명 정상 모드의 공간적 전이를 시각화한다. 결과는 아마이드 I 모드가 1 ps 이하의 시간 스케일로 이완되며, 이완 메커니즘이 해당 모드의 주변 환경에 민감함을 보여준다.

상세 분석

이 논문은 시토크롬 c(단백질)와 주변 수분이 결합된 시스템에서, 13C 또는 15N 등 동위원소로 표지된 아마이드 I 진동 모드의 에너지 이완 과정을 정량적으로 해석한다. 저자들은 시간 의존 섭동 이론(time‑dependent perturbation theory, TDPT)을 기반으로, 2차 및 3차 비선형 결합을 포함한 포텐셜 에너지 표면을 정상모드 기반으로 전개하였다. 이때, VER 속도는 Fermi’s golden rule 형태로 전이 행렬원소와 밀도 상태 함수의 곱으로 표현되며, 단백질 내부와 물 분자 각각이 제공하는 포톤(또는 포논) 모드와의 공명 조건을 명시적으로 계산한다.

핵심적인 기법은 “radial excitation function”과 “angular excitation function”을 도입해, 특정 아마이드 I 모드와 에너지 교환을 하는 정상모드들의 공간적 분포를 정량화한 점이다. 방사형 함수는 중심 아마이드와의 거리 r에 대한 에너지 전달 강도를, 각형 함수는 θ, φ 각도에 따른 방향성을 나타낸다. 이를 통해, 예를 들어 α‑헬릭스 내부에 위치한 아마이드 I는 주로 인접한 수소 결합 네트워크와 강하게 결합된 물 분자에 의해 이완되며, β‑시트 구조에서는 인접한 아미노산 잔기의 탄소‑탄소 결합 진동이 주요 에너지 수용체가 됨을 확인하였다.

또한, 단백질 기여와 물 기여를 각각 분리한 결과, 전체 VER 속도는 두 기여의 선형 합으로 근사될 수 있음을 보였다. 특히, 표지된 아마이드 I가 물에 노출된 경우 물의 저주파 회전·진동 모드와의 공명이 강화되어 이완 시간이 0.5 ps 수준으로 단축되는 반면, 내부에 매몰된 경우 단백질 내부의 고주파 정상모드와의 3차 비선형 결합이 지배적이며 이완 시간이 0.8–1.0 ps 범위에 머문다.

이러한 분석은 실험적 2D‑IR 스펙트로스코피에서 관찰되는 서브 피코초 이완 현상을 이론적으로 정당화한다. 특히, 환경 의존성(노출도, 주변 수소 결합 네트워크, 구조적 강성 등)이 VER 경로와 속도에 미치는 영향을 정량화함으로써, 단백질 내 국소적인 동역학을 파악하는 새로운 지표를 제공한다.