Efficient Proton Relay Orchestrated by Covalent Bond Switching of Active Amino Acids in Protein Channels

기술 요약: 단백질 채널 내 활성 아미노산의 공유 결합 전환에 의한 효율적인 양성자 릴레이

문제 제기
단백질은 질량 수송 및 에너지 전환을 포함하여 다양한 기능을 수행하기 위해 20가지 표준 아미노산에 의존하며, 이는 종종 양성자 전달과 결합되어 있다. 단백질의 구조적 안정성은 펩타이드 결합 및 비공유 상호작용을 통해 잘 이해되어 왔으나, 아미노산 측쇄의 말단 X–H 공유 결합(여기서 X는 O, S, N 또는 C)의 역할은 전통적으로 정적인 것으로 간주되어 왔다. 공유 결합 전환(예: hHv1의 Glu119 및 M2 채널의 His에서 발생하는)에 의해 매개되는 양성자 전달의 특정 사례들은 확인되었지만, 20가지 모든 아미노산 잔기가 양자 역학적 결합 전환을 통해 양성자와 어떻게 상호작용하는지에 대한 체계적인 이해는 부족한 실정이다. 본 연구가 다루는 근본적인 질문은 아미노산을 양성자 전달을 매개하는 고유한 능력에 따라 분류할 수 있는 보편적인 메커니즘이 존재하는가 하는 점이다.

방법론
저자들은 인간 전압 개폐 양성자 채널(hHv1) 내에서 핵심 잔기인 Glu119를 나머지 19가지 표준 아미노산으로 각각 치환하는 체계적인 돌연변이 시뮬레이션을 수행하였다.

• 시뮬레이션 기법: 양성자 릴레이 역학을 모델링하기 위해 양자 역학/분자 역학(QM/MM) 분자 역학(MD) 시뮬레이션을 활용하였다.
• 자유 에너지 계산: 서로 다른 잔기에서의 X–H 결합 전환에 대한 자유 에너지 장벽($\Delta G^\ddagger$)을 계산하기 위해 확장 적응형 편향 힘(eABF) 방법을 적용하였다. 자발적인 전환이 관찰되지 않는 잔기의 경우, 이러한 장벽을 추정하기 위해 가상의 경로를 구축하였다.
• 기술자(Descriptor) 분석: 관찰된 활성의 물리적 기원을 식별하기 위해, 국부적 친핵성 지수($N_k$, $f^-_k$), 정전기 전위 극값(ESPE), 원자 전하, 결합 쌍극자 모멘트를 포함한 11가지 물리화학적 기술자와 $\Delta G^\ddagger$ 사이의 절대 스피어만 순위 상관관계를 계산하였다.
• 일반화: 연구 결과는 호흡 복합체 I에서 검증되었으며, 세 가지 서로 다른 단백질 채널(hHv1, Otop1 C 도메인, 인플루엔자 A 바이러스 M2)을 분석하여 활성 잔기와 비활성 잔기의 방사형 분포를 조사하였다.

주요 결과

1. 아미노산의 분류: 본 연구는 결합 전환 능력에 기초하여 20가지 표준 아미노산의 이진 분류를 확립하였다:
   • 활성 잔기 (13종): Glu, Asp, Ser, Tyr, Thr, Cys, Asn, His, Lys, Trp, Arg, Gln, Phe. 이 잔기들은 O–H, S–H, N–H 또는 $sp^2$ C–H로 끝나는 측쇄를 가진다. 이들은 낮은 에너지 장벽($\Delta G^\ddagger < 20$ kJ mol$^{-1}$)의 X–H 결합 전환을 통해 양성자 릴레이를 매개한다.
   • 비활성 잔기 (7종): Leu, Ala, Val, Gly, Ile, Met, Pro. 이들은 $sp^3$ C–H 결로 끝난다. 이들은 결합 전환을 일으키지 않으며, 대신 인접한 물 분자 와이어를 통해 양성자 전달을 지원한다.
2. 메커니즘적 변이:
   • O–H/S–H 전환: Glu 및 Asp와 같은 잔기는 거의 장벽 없는 전환을 보인다. Cys (S–H)는 하이드록실 함유 잔기의 메커니즘을 모방한다.
   • N–H 전환: Lys, Trp, Asn, Gln, Arg는 단일 질소 부위에서 전환된다. 이미다졸 고리를 포함하는 His는 두 개의 질소 부위에서 협동적 메커니즘을 사용하여 잠재적으로 효율적인 장거리 릴레이를 가능하게 한다.
   • C–H 전환: Phenylalanine (Phe, $sp^2$ C–H)은 C–H 결합 전환을 통해 전달을 매개하는 반면, Methionine (Met, $sp^3$ C–H)은 물 와이어를 형성하기 위한 회전 재배열에 의존한다.
3. 활성 결정 요인: 결합 전환을 위한 에너지 장벽은 주로 양성자 수용 원자의 국부적 친핵성($N_k$, $f^-_k$)과 **정전기 전위 극값(ESPE)**에 의해 결정된다.
   • 활성 잔기는 음의 ESPE ($< -0.02$ a.u.)와 높은 친핵성을 나타내어 하이드로늄 이온($H_3O^+$)의 인력을 촉진하고 에너지 장벽을 낮춘다.
   • 비활성 잔기는 양의 ESPE ($> 0.02$ a.u.)와 낮은 친핵성을 나타내어 결합 전환을 억제하는 에너지 장벽을 형성한다.
4. 공간적 분포: hHv1, Otop1, M2 채널에서 활성 잔기는 양성자 릴레이를 매개하기 위해 중앙 기공 내에 우선적으로 농축되는 반면, 비활성 잔기는 구조적 안정화를 제공하기 위해 주변부에 집중된다. 이러한 분포-함수 관계는 단백질 채널의 구체적인 구조에 따라 세부적으로 달라진다.
5. 일반성: X–H 결합 전환 메커니즘은 호흡 복합체 I(His292 및 Lys282에서)에서 확인되었으며, 이 부위들을 다른 활성 잔기(Gln, Asn)로 돌연변이시켜도 지속된다. 이는 생물학적 시스템에서 효율적인 양성자 전달을 위한 일반적인 메커니즘임을 시사한다.

의의 및 주장
본 논문은 전하, 극성 또는 소수성에 기반한 전통적인 분류를 넘어, 공유 결합 전환 활성에 기초한 아미노산의 새로운 분류 기준을 확립했다고 주장한다. 저자들은 이 연구가 생명이 20가지 유형의 아미노산을 어떻게 활용하는지에 대한 근본적인 물리적 통찰을 제공한다고 단언한다. 즉, 13가지 활성 잔기는 낮은 장벽의 양자 역학적 결합 전환을 통해 효율적인 양성자 전달을 매개하도록 구조적, 화학적으로 최적화되어 있는 반면, 7가지 비활성 잔기는 안정화 역할을 수행한다는 것이다. 핵심 기술자로서 정전기 전위 극값을 식별한 것은 생체 에너지 효소 및 채널의 양성자 전달 메커니즘을 이해하기 위한 예측 도구를 제공한다.