Subconductance states in NMDA receptor variants are Ca2+ impermeable and act at the M2 loop

기술 요약: NMDA 수용체 변이체의 하위 전도 상태는 Ca2+ 불투과성이며 M2 루프에서 작용한다

문제 제기
NMDA 수용체(NMDAR)의 데 노보(de novo) 미스센스 변이는 심각한 신경 발달 장애와 연관되어 있다. 이러한 변이체들은 종종 수용체의 게이팅이나 이온 투과성을 변화시키지만, 구체적으로 어떤 메커니즘을 통해 기능 저하(특히 $P_{Ca}$의 감소)를 일으키는지에 대해서는 여전히 불분명하다. 선행 연구에서는 베스티뷸(vestibule, 외부 M3 게이트와 M2 포어 루프 사이의 영역)에 위치한 GluN1-M641I 변이체가 $P_{Ca}$를 감소시킨다는 것을 확인했으나, 이 효과의 구조적 근거는 알려지지 않았다. 또한, 하위 전도 상태(중간 전류 수준)와 $Ca^{2+}$ 투과성 사이의 관계 역시 명확하게 확립되지 않았다.

연구 방법론
저자들은 전기생리학, 분자 역학(MD) 시뮬레이션, 통계 분석을 결합한 다각적인 접근 방식을 사용하였다:

• 전기생리학: 인간 GluN1/GluN2A, GluN2B, GluN2C, GluN2D를 발현하는 HEK293T 세포를 대상으로 단일 채널 온셀 패치 클램프(on-cell patch-clamp) 기록을 수행하였다. 테스트된 변이체는 GluN1-M641 (I, L, V, T), 상동적인 GluN2 변이체(GluN2A의 V639I, GluN2D의 V667I), 그리고 N-사이트 아스파라긴 변이체(GluN1-N616S/A, GluN2A-N614S/A)를 포함한다.
• $Ca^{2+}$ 투과성 분석: 0에서 10 mM의 세포 외 $Ca^{2+}$로 전환할 때의 역전 전위($\Delta E_{rev}$) 변화를 측정하여 Lewis 방정식을 통해 $P_{Ca}/P_{Na}$ 비율을 계산하였다.
• 클러스터 분석: 단일 채널 트레이스를 최대 간격 가능도(maximum interval likelihood, MIL) 알고리즘을 사용하여 분석함으로써, 클러스터 열림 확률($P_{cluster}$, M3 게이트의 지표)과 클러스터 내 열림 확률($Cluster P_{open}$, M2 게이트의 지표)을 결정하였다.
• 분자 역학(MD): 열린 상태(구조 9ARE 기반)의 쥐 GluN1/GluN2B 구조체에 대해 시뮬레이션을 수행하여, 야생형(Met641), Leucine(M641L), Isoleucine(M641I) 변이체의 M2 및 M3 세그먼트와의 접촉 빈도와 포어 내부의 물 분포를 분석하였다.

주요 결과

1. GluN1-M641I 효과의 특이성: GluN1-M641 및 GluN2 서브유닛의 상동 위치에 있는 변이체들 중, 오직 GluN1-M641I만이 $Ca^{2+}$ 투과성을 유의미하게 감소시켰다. 해당 위치의 다른 변이체들(L, V, T) 및 상동적인 GluN2 변이체들(V639I, V667I)은 $P_{Ca}$를 변화시키지 않았다.
2. 하위 전도 상태의 유도: GluN1-M641I는 하위 전도 상태의 평균 열림 시간을 증가시키는 특징을 가진 견고한 하위 전도 상태를 독특하게 유도하였다 (Sub Eq. $P_o \approx 0.14$ vs. WT 0.02). 이 효과는 이소류신 측쇄에 특이적이었으며, M641 위치의 다른 변이체들은 유의미한 하위 전도를 생성하지 않았다.
3. 서브유닛 의존적 게이팅 효과: GluN1-M641I는 모든 GluN2 서브유닛에서 하위 전도를 유도하였으나, 전체 열림 확률($Total P_o$)에 미치는 영향은 달랐다. GluN2A와 함께 있을 때는 $Total P_o$가 변하지 않았으나, GluN2B, GluN2C, GluN2D와 함께 있을 때는 주 전도(main conductance)의 열림 확률 감소와 평균 닫힘 시간 증가로 인해 $Total P_o$가 급격히 감소하였다.
4. 상동 GluN2 변이체: GluN1-M641I와 대조적으로, GluN2A(V639I) 및 GluN2D(V667I)의 상동 이소류신 변이체들은 수용체 게이팅을 강력하게 강화(GoF)하였으나, 견고한 하위 전도 상태를 유도하지는 않았다.
5. 분자적 메커니즘 (MD 시뮬레이션): MD 시뮬레이션 결과, Ile641 측쇄는 Met641 또는 Leu641에 비해 M2 포어 루프(특히 N-사이트 아스파라긴 및 인접한 세린)와 더 빈번하게 상호작용하며, M3 세그먼트와의 상호작용은 적었다. 이러한 상호작용은 M2 루프 주변의 물 분자 배출과 상관관계가 있었다.
6. N-사이트 변이체를 통한 메커니즘 확인: N-사이트 아스파라긴 변이체들(GluN1-N616S, GluN2A-N614S) 또한 강한 하위 전도 상태를 유도하였으며 $P_{Ca}$를 유의미하게 감소시켰다.
7. 하위 전도와 $Ca^{2+}$ 투과성 간의 상관관계: 하위 전도와 주 전도 열림 확률의 비율과 $P_{Ca}$ 사이에 강한 음의 상관관계가 관찰되었다. 하위 전도 점유율을 증가시키는 변이체들(N-사이트 돌연변이 및 화합물 EU1622-240 포함)은 일관되게 감소된 $Ca^{2+}$ 투과성을 나타냈다.
8. 게이트 분석: 클러스터 분석 결과, GluN1-M641I는 $P_{cluster}$ (M3 게이트 지표)를 변화시키지 않으면서 $Cluster P_{open}$ (M2 게이트 지표)을 유의미하게 감소시켰으며, 이는 하위 전도 상태가 M2 루프의 변화로부터 기인함을 시사한다.

의의 및 결론
본 논문은 NMDAR의 하위 전도 상태가 대부분 $Ca^{2+}$에 불투과하다는 결론을 내린다. 저자들은 GluN1-M641I 변이체가 포어를 직접 차단하는 것이 아니라, M2 포어 루프와의 상호작용을 증가시켜 하위 전도 상태를 유도함으로써 $Ca^{2+}$ 투과성을 감소시킨다고 제안한다. 이러한 상호작용은 물 분자를 밀어내고 M2 게이트의 역학을 변화시킨다. 본 연구는 M2 루프가 전도 수준과 $Ca^{2+}$ 선택성을 조절하는 핵심 부위임을 강조한다. 나아가, 본 연구 결과는 NMDAR 변이와 관련된 임상 표현형이 $Na^+$ 흥분과 $Ca^{2+}$ 유입의 디커플링(decoupling)에 의해 유발될 수 있음을 시사한다. 즉, 하위 전도 상태를 유도하는 변이체들은 전체적인 채널 활성을 유지하거나 변화시키면서도 실질적으로 $Ca^{2+}$ 신호를 감소시킨다.