분자 프로그래밍의 강인성: 자연과의 게임 이론적 접근

초록

이 논문은 DNA 기반 화학 반응망(CRN)의 반응 속도 변동에 대한 강인성을 게임 이론으로 정량화한다. ‘자연’이 무작위로 반응 속도를 교란하는 상황을 모델링하고, 이를 통해 약 0.7 이상의 성공률을 유지하는 ‘근사 다수결’ CRN의 α‑강인성을 실험적으로 입증한다.

상세 분석

본 연구는 분자 프로그래밍에서 가장 실질적인 불확실성인 반응 속도 상수(k)의 변동을 ‘시계 속도(clock speed)’라는 개념으로 추상화하고, 이를 상대적 교란으로 보는 게임 이론적 프레임워크를 제시한다. 먼저, 무한 종 집합 S 위에 정의된 화학 반응망(CRN)을 정형화하고, Gillespie 알고리즘에 기반한 확률적 질량 작용(stochastic mass‑action) 의미론을 적용한다. CRN은 연속시간 마코프 체인으로 해석되며, 각 반응 ρ는 상태 x에서 rateₓ(ρ)=k·V^{1‑arity(ρ)}·∏_{species} x(species)^{arity} 로 표현된다.

게임 모델은 n명의 플레이어가 각각 CRN N_i를 전략으로 선택하는 n‑player 게임 G=(N,u)으로 정의된다. 여기서 N는 전략 공간의 직곱이며, u_i는 플레이어 i의 효용 함수로, 특정 트래젝터리 τ에 대해 u_i(τ)∈ℝ을 반환한다. ‘자연’ 플레이어는 효용이 0인 무관심 플레이어로 설정하고, 다른 플레이어들의 CRN이 촉매(C) 형태로만 상호작용하도록 제한한 ‘촉매 게임(catalytic game)’을 도입한다. 촉매는 반응 전후에 수량이 변하지 않으므로, 다른 플레이어는 촉매 농도를 통해 상대적 반응 속도만을 조절할 수 있다.

강인성 정의는 α‑robustness로, 플레이어 1이 자신의 CRN N₁을 사용했을 때, 다른 플레이어들의 존재가 기대 효용을 최소 α배 이하로 감소시키지 못한다는 조건이다. 수식적으로는
U₁( (N₁,N₂,…,N_n), ξ ) ≥ α·U₁( (N₁,∅,…,∅), ξ’ )
이며, ξ와 ξ’는 초기 상태 분포를 나타낸다.

구체적 사례로 ‘근사 다수결(approximate majority)’ 알고리즘을 구현한 CRN R을 사용한다. 원본 반응 R: 2X+Y→3X, X+2Y→3Y는 두 종 X와 Y가 서로 경쟁해 초기 과반수 종이 전체를 장악하도록 설계되었다. 이를 촉매 A, B를 도입해 R’으로 변형하고, A와 B의 초기 농도 a(0)=b(0)=100을 ‘자연’ CRN N_k: A↔B (전환율 k) 로 교란한다. k=10⁹인 경우, A와 B 사이의 전환이 매우 빠르게 일어나면서 a와 b가 무작위로 변동한다.

시뮬레이션은 MATLAB SimBiology를 이용해 10⁴개의 초기 조건(x(0),y(0))에 대해 10⁶번 반복 수행하였다. 결과는 ‘자연’이 존재하지 않을 때 성공 확률이 99%에 달하던 것이, a와 b가 교란될 경우에도 최소 70% 이상 유지됨을 보여준다. 즉, R’은 N_{10⁹}에 대해 α≈0.7‑robust함을 실증한다.

이 접근법의 장점은 (1) 강인성을 정량적 효용 함수로 명시함으로써 설계 목표와 직접 연결하고, (2) 악의적인 공격자뿐 아니라 자연적인 변동까지 포괄하는 일반화된 보안 모델을 제공한다는 점이다. 한계는 현재 효용이 이진 성공/실패에 국한되어 있어, 복합적인 성능 지표(예: 시간 복잡도, 에너지 소비)와의 통합이 필요하고, 다수 플레이어가 존재할 때의 내시 균형 분석이 미비하다는 점이다. 향후 연구는 전략적 균형 계산, 다중 목표 효용 설계, 그리고 실험실 수준에서의 DNA 구현 검증을 목표로 한다.