효모와 식물 뿌리의 ZIP 수송체를 통한 아연 흡수 조절 모델링

초록

이 논문은 효모(Saccharomyces cerevisiae)와 식물 뿌리(Arabidopsis thaliana)에서 아연을 흡수하는 ZIP 수송체의 조절 메커니즘을 미분방정식 기반 모델로 구현한다. 효모에서는 ZRT1·ZRT2와 전사인자 ZAP1의 상호작용을, 식물에서는 여러 ZIP 수송체와 bZIP 전사인자를 포함한 세 가지 가설 모델(활성자 단독, 이합체형 활성자, 활성자/억제자)을 비교·검증한다. 파라미터 추정과 시뮬레이션 결과, 효모에서는 ZAP1의 양성 피드백이 고유입 상황에서 시스템을 안정화시키며, 식물에서는 활성자/억제자 모델이 가장 견고한 아연 항상성을 제공한다는 결론에 도달한다.

상세 분석

본 연구는 아연이 세포 내에서 필수 미량원임과 동시에 과다 시 독성을 나타내는 양면성을 고려해, 세포 외부 아연 농도 변화에 대한 내부 아연 농도와 수송체 발현 수준의 동적 조절을 정량적으로 설명하고자 한다. 이를 위해 일반적인 ODE(Ordinary Differential Equation) 프레임워크를 구축하고, 각각의 시스템에 맞게 변수와 파라미터를 정의하였다.

효모 모델에서는 세 가지 주요 변수(세포 내 아연 농도 Z, ZRT1·ZRT2 단백질 양 T₁, T₂, 그리고 ZAP1 전사인자 양 A)를 사용하였다. ZRT1은 고친화성, ZRT2는 저친화성 수송체로 알려져 있어, 각각의 Km 값이 다르게 설정되었다. ZAP1은 아연 결핍 시 활성화되어 ZRT1·ZRT2 유전자를 전사 촉진하고, 동시에 자체 발현에 대한 양성 피드백을 갖는다. 모델식은 아연 유입(수송체에 의한 Michaelis‑Menten 형태), 아연 배출·희석, ZAP1에 의한 전사 활성화(Hill 함수), 그리고 ZAP1 자체 합성·분해를 포함한다. 파라미터는 문헌값과 실험 데이터(아연 농도와 전사체 발현량)로부터 비선형 최소제곱법을 이용해 추정하였다. 추정 결과는 ZRT1이 저농도에서 주된 역할을, ZRT2가 고농도에서 보조적 역할을 수행한다는 기존 보고와 일치한다. 특히, ZAP1의 양성 피드백 계수 β가 0.3~0.5 범위에서 시스템이 과도한 아연 유입 시 발현 과잉을 방지하고, 안정적인 고정점을 유지함을 민감도 분석을 통해 확인하였다. 이는 피드백이 단순히 발현을 증폭하는 것이 아니라, 비선형적인 억제 메커니즘과 결합해 ‘스위치‑온/오프’ 형태의 이중 안정성을 제공한다는 점을 시사한다.

식물 뿌리 모델에서는 ZIP 수송체들의 발현 조절이 아직 명확히 규명되지 않아, 세 가지 가설 모델을 제시하였다. 첫 번째는 전사인자 활성자만 존재하는 단순 모델로, 아연 결핍 시 활성자가 직접 ZIP 유전자를 전사 촉진한다. 두 번째는 활성자가 이합체(dimer) 형태로 작용해 전사 활성도가 Hill 계수 n>1을 갖는 경우이며, 이는 협동적 결합을 통해 급격한 전이 현상을 재현한다. 세 번째는 활성자와 억제자 두 종류가 존재해, 활성자는 ZIP 전사를 촉진하고 억제자는 반대로 억제하는 복합 피드백 구조이다. 각 모델은 동일한 ODE 구조(아연 유입, 세포 내 아연 농도, 전사인자/억제자 농도)와 동일한 실험 데이터(다양한 외부 아연 농도에서 측정된 뿌리 내 아연 함량 및 ZIP mRNA 발현)로 피팅되었다.

모델 적합도는 AIC와 교차 검증을 통해 비교했으며, 활성자/억제자 모델이 가장 낮은 AIC 값을 보였다. 시뮬레이션 결과, 이 모델은 외부 아연 농도가 급격히 변하더라도 내부 아연 농도가 5~10 µM 범위 내에서 안정적으로 유지되는 ‘강인한 항상성’를 제공한다. 억제자는 고농도 아연에서 과도한 ZIP 발현을 억제함으로써 독성을 방지하고, 동시에 낮은 농도에서는 억제 효과가 약해져 활성자가 충분히 발현을 촉진한다. 반면, 단순 활성자 모델은 고농도에서 과잉 발현이 발생해 내부 아연이 급증하고, 이합체 모델은 파라미터 튜닝이 매우 민감해 실험적 변동성을 반영하기 어렵다.

전반적인 민감도 분석에서는 억제자 합성 속도 k_i와 억제자-활성자 결합 해리 상수 K_d가 시스템 안정성에 가장 큰 영향을 미치는 것으로 나타났다. 특히 k_i가 0.1–0.3 min⁻¹ 범위에서 내부 아연 농도가 5 µM 이하로 유지되었으며, K_d가 10⁻⁶ M 이하일 때 억제 효과가 충분히 강해 과잉 유입을 차단한다. 이러한 결과는 식물 뿌리에서 bZIP 전사인자와 함께 억제성 조절 인자가 존재한다는 가설을 뒷받침한다.

결론적으로, 효모에서는 ZAP1의 양성 피드백이 고유입 상황에서 시스템을 안정화시키는 역할을 하며, 식물에서는 활성자와 억제자가 공동으로 작용하는 복합 피드백이 아연 항상성을 가장 효과적으로 유지한다는 점을 모델링을 통해 입증하였다. 이러한 정량적 모델은 향후 유전자 편집이나 작물 개량 전략에서 아연 효율성을 최적화하는 데 활용될 수 있다.