생물세포의 기계·화학 결합 모델링 언어와 시뮬레이션 엔진

초록

본 논문은 세포 내·외부의 복합적인 기계·화학·전기 현상을 자연스럽게 기술할 수 있는 도메인 전용 모델링 언어(MML)와 이를 컴파일하여 메쉬‑프리 라그랑지안 시뮬레이터에서 실행할 수 있는 프레임워크를 제안한다. 객체와 프로세스라는 두 기본 개념을 도입해 입자, 링크, 농도, 영역 등 새로운 타입을 정의하고, 연속·이산 변환을 동시에 기술한다. 기존 SBML·Petri Net·MD 등과 비교해 동적 기하와 힘·스트레스 모델링을 자연스럽게 지원한다.

상세 분석

이 논문은 생물학적 세포를 활성 물질(active matter)로 보고, 세포가 경험하는 기계·화학·전기적 자극을 동시에 모델링할 필요성을 강조한다. 이를 위해 저자들은 ‘Mechanica Modeling Language (MML)’라는 새로운 도메인 특화 언어를 설계하였다. MML은 두 가지 기본 구성 요소인 **객체(objects)**와 프로세스(processes) 로 모델을 기술한다. 객체는 스칼라, 정수, 불리언 등 기본 타입뿐 아니라 입자(particle), 링크(link), 농도(concentration), 양(amount), 공간 영역(spatial region), 공간 필드(spatial field)와 같은 복합 타입을 포함한다. 각 객체는 고유 식별자를 가지며, 상태 공간(state space)을 통해 허용 가능한 값 범위를 정의한다. 특히 입자는 물리적 물질의 ‘패킷’으로서 위치와 질량, 부피, 방향 등을 가질 수 있고, 링크는 입자들 사이의 영구 혹은 가변적인 결합을 표현한다. 링크는 전통적인 MD의 결합 포텐셜을 확장해, 사용자가 임의의 힘 함수를 정의할 수 있게 함으로써 탄성 스프링부터 복잡한 다극자 상호작용까지 포괄한다.

프로세스는 전통적인 컴퓨팅 의미의 프로세스와 달리 연속 변환(continuous) 과 이산 변환(discrete) 두 형태를 제공한다. 연속 프로세스는 화학 반응식과 유사하게 반응물과 생성물 사이의 변환 속도(rate)를 정의하고, 입자에 부착된 농도·양 변수는 입자 이동에 따라 라그랑지안 방식으로 수송된다. 이산 프로세스는 SBML 이벤트와 유사하게 조건(predicate)과 확률(probability)을 통해 순간적인 상태 변화를 일으키며, 객체 생성·소멸도 지원한다. 이러한 설계는 기존 시스템 생물학 언어(SBML, Petri Net)가 공간적·기계적 변수를 다루지 못하는 한계를 극복한다.

컴파일러는 MML 스크립트를 파싱해 내부 중간 표현(IR)로 변환하고, 메쉬‑프리 라그랑지안 시뮬레이션 엔진에 적합한 실행 코드를 생성한다. 엔진은 입자 기반의 Lagrangian 방법을 사용해 물리적 힘과 화학 확산을 동시에 계산한다. 특히 MaterialRegion 개념을 도입해 복잡한 경계면(예: 세포막) 내부에 입자를 자동으로 채우고, 농도·양 변수를 영역 전체에 균일하게 할당하거나, 부피 변화에 따라 자동 스케일링한다. 이는 기존 입자 기반 시뮬레이터가 고정된 입자 수와 정적인 기하만 지원하던 문제를 해결한다.

논문은 기존 접근법을 체계적으로 비교한다. 전통적인 MD는 원자 수준에서만 가능하고, 반응성 MD는 입자 수 고정·힘 종류 제한이라는 제약이 있다. Stochastic particle 시뮬레이터(MCell, Smoldyn 등)는 고정된 기하와 화학 반응만 지원한다. tDPD는 입자에 화학 상태를 부여하지만 모델링 언어가 부재한다. 반면 MML은 tDPD의 장점을 계승하면서도 반응식·이산 이벤트·동적 링크·공간 영역을 하나의 언어로 통합한다.

핵심 기여는 다음과 같다. (1) 생물학적 개념에 기반한 도메인 전용 언어 설계, (2) 객체‑프로세스 기반 하이브리드 연속·이산 변환 모델링 프레임워크, (3) 메쉬‑프리 라그랑지안 엔진을 통한 효율적인 시뮬레이션, (4) 기존 시스템 생물학·MD·tDPD 등과의 기능적 차별성 입증. 이러한 설계는 세포 수준에서부터 조직·기관 수준까지 다양한 스케일의 기계·화학 결합 현상을 통합적으로 연구할 수 있는 기반을 제공한다.