CRN 공개 화학 반응 네트워크의 체계적 탐색

초록

**
Alloy 기반 형식 모델링을 이용해 작은 화학 반응 네트워크(CRN)를 자동 열거하고, 촉매·대사·시소우 구조와 같은 생물학적 제약을 적용한다. 특히 피드‑포워드·비경쟁 조건을 만족하는 비율‑독립 CRN을 탐색해 max, minmax, abs, ReLU와 같은 아날로그 연산을 수행하는 최소 규모의 회로를 발견한다.

**

상세 분석

**
이 논문은 형식 검증 도구인 Alloy를 화학 반응 네트워크(CRN)의 구조적·동적 제약을 기술하는 언어로 활용한다는 점에서 혁신적이다. 기존 CRN 설계는 수동적 직관이나 제한된 탐색에 의존했지만, 저자들은 Alloy의 선언적 제약과 SAT 기반 해석기를 이용해 지정된 범위 내에서 모든 가능한 네트워크를 완전 열거한다. 일반적인 CRN 모델을 ‘Species’와 ‘Reaction’이라는 두 개의 추상 시그니처로 정의하고, 반응물·생산물 시퀀스를 통해 다중 스토이키오메트리를 자연스럽게 표현한다. 기본 제약(반응물 최소 1개, 반응 중복 금지, 반응물·생산물 구분, 모든 종 사용) 외에도 ‘Elementary’, ‘Catalytic’, ‘Metabolic’ 등 계층적 서브클래스를 정의함으로써 생물학적 현실성을 단계적으로 강화한다. 특히 촉매 반응을 정의할 때는 반응물과 생산물에 공통 종이 존재하도록 하고, 대사 네트워크는 종이 촉매이면 항상 촉매로만 등장하도록 하는 이중 제약을 적용한다.

시소우(seesaw) 게이트와 같은 DNA 스트랜드 디스플레이스먼트 회로를 모델링할 때는 ‘Gate’와 ‘Strand’라는 추가 시그니처를 도입해 물리적 결합 구조를 추상화한다. 이렇게 계층화된 모델은 새로운 제약을 삽입하거나 기존 제약을 수정하는 것이 용이해, 연구자가 관심 있는 특정 설계 공간을 빠르게 정의하고 탐색할 수 있게 한다.

핵심적인 과학적 기여는 ‘비율‑독립’(rate‑independent) CRN의 구조적 조건을 명시하고, 이를 Alloy 제약으로 구현한 점이다. 저자들은 두 가지 구조적 속성, 즉 feed‑forward(반응 순서가 존재해 뒤의 반응이 앞의 반응이 만든 종을 소비하지 않음)와 non‑competitive(한 종이 소비되면 다른 반응의 반응물에 나타나지 않음)를 정의한다. 이 두 조건은 CRN이 반응 속도 상수에 무관하게 동일한 정적 평형에 수렴하도록 보장한다. 논문 부록에서 이론적 증명을 제공하며, 실제 탐색에서는 이러한 제약을 만족하는 최소 규모의 네트워크를 찾는다.

예시로 제시된 max 연산 CRN은 네 개의 반응으로 구성되며, 입력 A와 B의 초기 농도를 합산하고 최소값을 구한 뒤 차를 계산해 최대값을 출력한다. 이 구조는 기존 연구와 일치하지만, Alloy 기반 탐색을 통해 이것이 최소 규모임을 증명한다. 더 나아가, abs(절대값)와 minmax(최소·최대 동시 출력)와 같은 함수는 이전에 보고된 바 없으며, 저자들은 각각 5~6개의 반응으로 구현 가능한 최소 네트워크를 발견한다. ReLU 함수 역시 4개의 반응으로 구현 가능함을 보여, 비율‑독립 CRN이 인공 신경망의 핵심 활성화 함수까지 자연스럽게 구현할 수 있음을 시사한다.

성능 측면에서 Alloy는 탐색 공간이 급격히 커지는 경우에도 스코프(bound)를 조절해 실용적인 시간 안에 해를 찾는다. 예를 들어, 2종·2반응 스코프에서 수십 개의 후보를 빠르게 생성하고, 46종·46반응 스코프에서도 수백 개의 후보를 검증한다. 이는 전통적인 수치 시뮬레이션 기반 탐색보다 훨씬 효율적이며, 설계자가 직관적으로 정의한 제약을 그대로 유지한다는 장점이 있다.

마지막으로, 저자들은 구현한 Alloy 모델과 탐색 도구를 오픈소스로 공개하여, 다른 연구자들이 자신의 도메인에 맞는 제약을 추가하거나 새로운 기능을 탐색할 수 있도록 한다. 이는 형식 방법론을 화학·생물공학 설계에 적용하려는 커뮤니티에 큰 영향을 미칠 것으로 기대된다.

**