대사 흐름의 효소 비용은 로그 농도에서 볼록함

초록

대사 경로의 목표 플럭스를 달성하기 위해 필요한 효소 양을 로그 농도 변수로 표현하면, 효소 비용 함수가 해당 변수 공간의 다각형(폴리토프) 위에서 볼록함을 보인다. 따라서 효소 비용 최소화는 볼록 최적화 문제로 풀 수 있다.

상세 분석

본 논문은 효소 비용 최소화(Enzyme Cost Minimization, ECM)를 수학적으로 정형화하고, 그 핵심이 되는 비용 함수가 로그 농도 공간에서 볼록(convex)함을 증명한다. 먼저, 가역 Michaelis‑Menten 형태를 포함한 일반적인 가역 속도식을 ‘분리 가능한(rate separable)’ 형태 v = ε·k_cat·η_th·η_kin·η_reg 로 분해한다. 여기서 η_th는 열역학적 효율(역반응에 의한 손실), η_kin은 포화 및 알로스테릭 효과를 포함한 동역학적 효율을 나타낸다. 효소 요구량 ε는 목표 플럭스 v와 효소‑특이 비용 h를 곱한 형태 ε = v/(k_cat·η_th·η_kin·η_reg) 로 표현되며, 이를 비용 함수 y = h·ε 로 쓰면 y는 각 반응별로 동일한 구조를 가진 항들의 합이 된다.

다음으로, 메타볼라이트 농도 로그값 x_i = ln c_i 로 정의한 뒤, 모든 반응의 열역학적 제약(Θ_l·v_l > 0)과 물리적 농도 구간(x_min ≤ x ≤ x_max)을 결합하면, 가능한 메타볼라이트 프로파일이 다각형 P(폴리토프)로 형성된다. 이 폴리토프는 E‑face(반응이 평형에 도달하는 면)와 P‑face(농도 상한·하한에 닿는 면)로 구성된 볼록 집합이다.

핵심 증명은 η_th와 η_kin이 각각 로그 농도에 대해 로그-선형 혹은 로그-볼록 형태를 가지며, 그 역수(비용에 등장)는 볼록함을 유지한다는 점이다. 구체적으로 η_th = 1 − e^{−Θ}는 Θ가 로그 농도의 선형 결합이므로 η_th^{-1}는 볼록, η_kin은 분모가 양의 다항식이므로 역수도 볼록이다. 따라서 y(x) = ∑_l h_l·v_l/(k_cat,l·η_th,l·η_kin,l) 가 P 전체에서 볼록함을 만족한다.

볼록성은 최적점이 전역 최적임을 보장하고, 표준적인 볼록 최적화 알고리즘(예: interior‑point, gradient descent)으로 효소 수준과 메타볼라이트 농도를 동시에 계산할 수 있게 한다. 또한, 비용 함수가 엄격히 볼록하지 않을 경우(예: 일정 구간에서 상수) 정규화 항 y_reg = ∑ γ_i (x_i − x_i^0)^2 등을 추가해 엄격한 볼록성을 회복할 수 있다.

이론적 결과는 대사 공학에서 효소 과다 발현을 피하고, 진화적 관점에서 효소 비용이 최소화된 경로가 선택될 가능성을 설명한다. 또한, 기존 FBA와 달리 실제 효소 동역학과 열역학을 직접 반영함으로써 보다 현실적인 대사 설계가 가능해진다.