대사 최적화 원리로 보는 E. coli 필수 유전자 탐색

초록

이 논문은 질량 균형이나 플럭스 정착을 가정하지 않고, 전체 성장률을 최대로 하는 제약 기반 모델을 제시한다. E. coli 대사망에 적용한 결과, 기존의 정지 플럭스 해와 유사한 해를 얻으며 일부 대사물질은 비정상적인 생산을 보인다. 보존 대사 풀은 성장률과 플럭스 변동성에 핵심적인 역할을 하며, 변동 폭이 작은 ‘고정된’ 플럭스와 연관된 유전자는 실험적으로 필수 유전자로 확인된다.

상세 분석

본 연구는 기존의 Flux Balance Analysis(FBA)와 달리 질량 보존식과 플럭스 정지성을 강제하지 않는다. 대신, 전체 시스템의 성장률을 전역적인 목적 함수로 설정하고, 모든 반응 플럭스가 양의 실수값을 가질 수 있도록 제약을 부여한다. 이 접근법은 선형 계획법이 아닌 비선형 최적화 문제로 전환되며, 대사망의 토폴로지와 스토이키오메트리를 그대로 활용한다는 점에서 계산 효율성을 유지한다. E. coli의 대사망(약 1000개의 반응, 800여 개의 대사물질)을 대상으로 적용했을 때, 최적화된 해는 전통적인 FBA가 예측한 정지 플럭스와 거의 일치하지만, 특정 대사물질(예: 아미노산 전구체, 전자 수용체)의 경우 질량 균형이 깨져 순생산이 발생한다. 이는 실제 세포가 성장 과정에서 외부 환경에 따라 일정량의 대사산물을 비축하거나 배출한다는 생물학적 현상을 반영한다는 해석이 가능하다.

또한, 보존 대사 풀(conserved metabolic pools)이 모델 내에서 중요한 역할을 수행한다는 점을 강조한다. 보존 풀은 특정 대사물질들의 총량이 일정하게 유지되는 집합으로, 이들의 존재는 플럭스 공간을 크게 축소시켜 성장률에 대한 민감도를 높인다. 실험적으로 확인된 바와 같이, 보존 풀에 속한 반응은 플럭스 변동성이 낮아 ‘고정된(frozen)’ 플럭스로 분류된다. 이러한 고정된 플럭스와 연결된 유전자는 대사망에서 핵심적인 역할을 수행하며, 유전자 삭제 실험에서 필수성(essentiality)이 높은 것으로 나타난다.

논문은 또한 다양한 환경(탄소원, 질소원, 산소 가용성)에서 모델을 시뮬레이션하여, 플럭스 분포가 환경에 따라 어떻게 재배열되는지를 분석한다. 결과적으로, 환경 변화에 따라 일부 보존 풀은 유지되지만, 다른 풀은 해체되거나 새로운 풀로 재구성된다. 이는 대사망이 유연하게 재배열될 수 있음을 시사한다.

마지막으로, 모델이 예측한 플럭스 통계가 실험적 메타볼로믹스 데이터와 높은 상관관계를 보이며, 특히 고정된 플럭스에 해당하는 반응들의 흐름이 실험적으로 측정된 값과 거의 일치한다는 점을 강조한다. 이는 질량 균형을 강제하지 않으면서도 실제 세포의 대사 흐름을 정확히 포착할 수 있음을 입증한다.