시토크롬 c 산화효소의 전자‑양성자 연동 펌핑 메커니즘

초록

본 연구는 시토크롬 c 산화효소에 네 개의 레독스 중심과 네 개의 양성자 결합 부위를 포함한 간단한 모델을 제시한다. 단일 전자 주입에 의해 시작되는 전자와 양성자 전달을 전기적 결합으로 기술하고, 마스터 방정식으로 확률 흐름을 기술한다. 합리적인 파라미터 범위에서 효율적인 양성자 펌핑이 가능함을 보이며, 실험에서 보고된 네 단계의 동역학적 과정을 모두 재현한다. 또한 막 전위와 온도에 따른 펌핑 효율 변화를 계산하고, 레독스 구동 양성자 이동 메커니즘을 논의한다.

상세 분석

이 논문은 시토크롬 c 산화효소(CcO)의 전자 전달과 양성자 펌핑을 동시에 기술할 수 있는 최소 모델을 구축한다. 네 개의 레독스 중심(heme a, heme a₃, Cu_A, Cu_B)과 네 개의 양성자 결합 부위(D, K, H, 그리고 외부 전위에 민감한 E 부위)를 도입함으로써, 전자와 양성자 사이의 정전기적 상호작용을 명시적으로 포함한다. 모델은 각 상태를 전자와 양성자 점유 여부로 정의하고, 전자 주입 후 시간에 따라 상태 전이 확률을 기술하는 마스터 방정식 집합을 도출한다. 전이율은 Marcus‑type 전자 전이 이론과 프로톤 전이의 전기화학적 구배를 결합해 계산되며, 전압(V_m)과 온도(T)의 함수로 표현된다.

주요 결과는 네 개의 실험적 동역학 단계(µs, ms, s, 수초 규모)가 파라미터 선택에 따라 자연스럽게 나타난다는 점이다. 초기 단계에서는 전자가 Cu_A에서 heme a₃로 빠르게 이동하면서 K 부위에 양성자를 결합시키고, 이어서 D 부위로 양성자 전달이 진행된다. 두 번째 단계에서는 전자와 양성자 사이의 정전기적 결합이 강화되어 양성자 펌핑 효율이 급격히 상승한다. 세 번째와 네 번째 단계는 막 전위에 의해 조절되는 느린 양성자 방출과 재결합 과정을 반영한다.

전압 의존성 분석에서는 양성자 펌핑 효율이 약 +100 mV에서 최댓값을 보이며, 과도한 전압(> +150 mV)에서는 역전 현상이 나타나 전자 흐름에 비해 양성자 역방향 이동이 증가한다. 온도 변화에 대해서는 활성화 에너지(E_a)가 약 12 kJ·mol⁻¹로 추정되어, 온도가 상승할수록 전이율이 가속화되지만 효율은 약간 감소한다는 점을 보여준다. 이러한 결과는 CcO가 전자 흐름을 전기화학적 에너지로 변환하는 과정에서 전자‑양성자 커플링이 핵심적인 역할을 함을 뒷받침한다.

모델의 한계는 실제 효소가 갖는 복잡한 단백질 구조와 물 분자 네트워크를 단순화했기 때문에, 정량적 파라미터(전이율 상수, 전기용량 등)의 정확성은 실험적 보정이 필요하다는 점이다. 그럼에도 불구하고, 최소 모델이 실험 데이터와 일치하는 네 단계 동역학을 재현한다는 점은 전자‑양성자 연동 메커니즘을 이해하는 데 강력한 이론적 토대를 제공한다.