점 돌연변이가 리듬 수용체 RyR2 에너지 전도 경로를 파괴한다

초록

본 연구는 탄성 네트(Elastic Net) 모델을 이용해 심장형 Ryanodine 수용체(RyR2) N‑말단 도메인(173 aa)에서 에너지 전도 경로를 규명하고, 진화적으로 보존된 잔기들이 그 경로와 밀접히 연관됨을 확인하였다. 특히 E171·E173 두 글루탐산이 칼슘 결합 부위 후보로 제시되었으며, 질환을 일으키는 A77V 변이가 해당 경로를 붕괴시켜 에너지 응답성을 상실함을 보여준다.

상세 분석

이 논문은 먼저 RyR2의 N‑말단 도메인(173 aa)의 고해상도 X‑선 구조(PDB 3IM5)와 A77V 변이체(PDB 3IM7)를 기반으로, 단순화된 탄성 네트(Elastic Net) 모델을 적용하였다. 모델은 각 잔기의 열진동에 의해 인접 12개의 잔기와 거리‑에너지 상관관계를 계산하고, 전체 단백질에 대한 평균 에너지 응답 U_i를 도출한다. U_i 값이 높은 잔기들은 서로 강한 상관 U_ij를 보이며, 공간적으로 인접한 잔기들 사이에 연속적인 에너지 전도 경로를 형성한다.

계산 결과, 에너지 전도 경로에 포함된 잔기들은 대부분 기존 문헌에서 진화적으로 보존된 것으로 알려진 잔기들과 일치했으며, 이들은 두 개의 Mir 도메인(124‑178, 164‑217) 내부 혹은 그와 수소 결합으로 연결된 위치에 집중되어 있었다. 특히 E171과 E173은 경로의 말단에 위치하면서도 높은 에너지 응답을 보였는데, 이는 IP3R(1N4K)에서 동일한 위치에 있는 E283·E285가 칼슘 결합 부위임을 고려하면 RyR2에서도 칼슘 결합 부위 후보로 해석될 수 있다.

A77V 변이의 구조적 효과를 분석하면, 알라닌을 발린으로 치환함으로써 인접 부위에 두 개의 메틸기가 추가되어 국소적인 입체장애와 거리 변화를 초래한다. 차이 지도(ΔR_ij) 분석에서는 15,051개의 잔기 쌍 중 다수가 통계적으로 유의한 거리 변화를 보였으며, 이러한 변형이 경로상의 핵심 잔기들(예: 45, 78, 112 등)의 U_i와 U_ij 값을 급격히 감소시켜 에너지 전도 흐름을 끊었다. 변이체에서 붕괴된 경로는 그림 6의 빨간색 잔기로 표시되었으며, 이들 잔기는 모두 Mir 도메인과 직접적인 수소 결합 네트워크를 이루고 있었다.

IP3R에 대한 동일한 탄성 네트 분석을 수행한 결과, RyR2와 유사하게 에너지 전도 경로가 Mir 도메인 주변에 형성되고, 보존된 잔기들이 경로와 겹치는 현상이 재현되었다. 이는 두 채널이 구조적·기능적으로 공유하는 에너지 전도 메커니즘을 시사한다.

결론적으로, 탄성 네트 모델은 복잡한 대형 이온 채널에서 기능적 에너지 전도 경로를 효율적으로 추출할 수 있음을 입증했으며, A77V와 같은 질환 유발 점 돌연변이가 이 경로를 파괴함으로써 칼슘 누출 및 부정맥을 촉발하는 분자적 메커니즘을 설명한다.