생체분자 시스템을 위한 모듈형 본드그래프 모델링과 분석

초록

본 논문은 본드그래프를 이용해 열역학적으로 일관된 계층적 생체분자 모델을 구축하고, 이를 통해 닫힌 시스템과 열린 시스템의 변환, 모듈성(계산적·행동적) 및 레트로액티비티를 정량적으로 분석한다. MAPK 캐스케이드를 사례로 제시하며, ATP‑ADP‑Pi와 같은 전원 공급이 레트로액티비티 감소에 필수임을 보인다.

상세 분석

본드그래프는 에너지 흐름을 ‘노력(effort)’과 ‘흐름(flow)’이라는 쌍으로 표현하고, 전력(p=ef)이라는 공통 통화를 통해 전기·기계·화학 등 이질적인 도메인을 통합한다. 화학 반응에서는 노력을 화학 퍼텐셜(µ), 흐름을 몰 흐름(v)으로 정의하고, C‑요소는 물질량 저장, Re‑요소는 두 포트 비선형 저항으로 반응 속도를 모델링한다. 스토이키오메트리 행렬(N_f, N_r)과 TF‑요소를 이용해 복잡한 네트워크를 하나의 연결 구조로 압축할 수 있다.

논문은 먼저 닫힌 시스템을 본드그래프로 기술한 뒤, ‘체스트랫(chemostat)’과 ‘플로우스트랫(flowstat)’이라는 외부 포트를 추가해 열린 시스템으로 전환한다. 이 과정에서 물질 보존과 에너지 보존이 각각 어떻게 유지되는지를 명확히 보여준다.

모듈성에 대해서는 두 가지 개념을 구분한다. 계산적 모듈성은 각 모듈이 물리·열역학적으로 정확히 연결될 수 있음을 의미하며, 이는 본드그래프의 포트 연결 규칙에 의해 자동 보장된다. 반면 행동적 모듈성은 모듈 고유의 동적 특성(예: 초감도, 스위칭) 이외부와 연결될 때도 유지되는지를 말한다. 생물학적 네트워크는 일반적으로 일방향 버퍼가 없기 때문에 레트로액티비티가 발생하고, 이는 모듈의 입력‑출력 특성을 왜곡한다. 본드그래프는 레트로액티비티를 전력 흐름의 불균형으로 해석하고, 이를 최소화하기 위해 ATP‑ADP‑Pi와 같은 전원 공급을 도입해야 함을 제시한다.

선형화 단계에서는 비선형 ODE를 테일러 전개해 작은 신호 근처에서 선형 근사 모델을 얻고, 이를 블록다이어그램 형태로 변환한다. 이렇게 하면 제어 이론의 피드백, 감도, 안정성 분석 도구를 직접 적용할 수 있다. 특히, 본드그래프 기반 선형 모델은 각 모듈의 전달함수와 상호 연결 구조를 명시적으로 드러내어 레트로액티비티가 피드백 루프 형태로 나타나는 것을 확인한다.

마지막으로 MAPK 캐스케이드를 사례로 삼아, 전통적인 미카엘리스‑멘텐 근사 대신 본드그래프 기반 완전 비선형 모델을 구축한다. 각 단계(Raf→MEK→ERK)를 독립적인 모듈로 정의하고, ATP‑ADP‑Pi 전원을 삽입해 레트로액티비티를 억제한다. 결과적으로 초감도 특성이 유지되면서도 전체 시스템이 열역학적으로 일관된 형태를 갖는다.

이 논문의 핵심 기여는 (1) 생화학 반응을 본드그래프 형식으로 체계화한 이론적 틀, (2) 닫힌·열린 시스템 전환을 위한 체스트랫·플로우스트랫 개념, (3) 계산적·행동적 모듈성을 구분하고 레트로액티비티를 전력 흐름 관점에서 해석한 점, (4) 선형화와 블록다이어그램을 통한 제어‑이론 적용 가능성을 제시한 점이다. 이러한 접근은 복잡한 신호전달 네트워크를 설계·분석하는 새로운 패러다임을 제공한다.