프로토셀 안에서 촉매 반응망의 확률적 진화 모델

초록

이 논문은 반투과성 막을 가진 단순 지질 구체 안에서 촉매 반응망(CRS)의 동역학을 확률적 시뮬레이션으로 분석한다. Gillespie 알고리즘을 이용해 종의 생성·소멸과 반응 네트워크의 변화를 추적하고, 작은 프로토셀 내부에서 발생하는 확률적 효과와 자가촉매 집합(ACS)의 출현 조건을 탐구한다. 성장·분열은 제외했지만, 막 투과성, 분자 크기 제한, 초기 농도 분포 등이 네트워크 진화에 미치는 영향을 제시한다.

상세 분석

본 연구는 기존의 두 축, 즉 (i) 고정된 자가복제 분자 집합을 전제로 한 프로토셀 모델과 (ii) 무한히 큰 개방계(CSTR)에서 무작위 분자 집합의 촉매 반응망을 다루는 모델 사이의 간극을 메우고자 한다. 이를 위해 저자들은 가장 단순한 프로토셀 구조—반투과성 막이 존재하고, 모든 화학 반응이 내부 수용액에서 일어나는 구형 지질 구체—를 가정하고, 성장·분열을 무시한 고정 부피 시스템을 설정한다. 이때 막은 일정 길이 이하의 분자만을 자유롭게 통과시키며, 통과 가능한 종은 내부와 외부 농도가 즉시 평형을 이루는 ‘adiabatic’ 가정을 둔다. 이러한 가정은 실제 프로토셀에서 발생할 수 있는 삼투압, 막 표면 촉매 작용, 비대칭적인 물질 흐름 등을 단순화하지만, 핵심적인 질문—‘왜 살아있는 시스템은 구획화가 필요할까?’—에 집중할 수 있게 한다.

모델의 화학적 기반은 Kauffman의 무작위 폴리머 생성 규칙을 차용한다. 단일 문자(모노머)와 문자열(폴리머)로 이루어진 종 집합 X = {x₁,…,x_N}을 정의하고, 두 가지 기본 반응, 즉 ‘절단(cleavage)’과 ‘결합(condensation)’을 고려한다. 모든 반응은 최소 길이 이상의 종에 의해 촉매될 때만 발생하도록 제한함으로써, 촉매가 없는 반응은 에너지 장벽이 너무 높아 실질적으로 일어나지 않는다고 가정한다. 이때 촉매 종은 다중 반응을 촉진할 수 있으며, 하나의 반응도 다중 촉매에 의해 가속될 수 있다.

동역학은 Gillespie 알고리즘을 적용해 확률적 사건 발생 시간을 계산한다. 이는 전통적인 ODE 기반 모델이 놓치기 쉬운 ‘종의 소멸’과 ‘네트워크 축소’를 자연스럽게 포함한다. 특히, 작은 부피(1 µm³ 이하)에서는 분자 수가 수십 개 수준으로 떨어져, 통계적 변동이 크게 작용한다. 저자들은 이러한 저농도 상황에서 (1) 초기 조건의 무작위성, (2) 반응 순서에 따른 경로 의존성, (3) 서로 다른 화학 규칙(예: RAF 집합)의 존재가 네트워크 진화에 미치는 영향을 체계적으로 실험한다.

시뮬레이션 결과는 다음과 같은 핵심 통찰을 제공한다. 첫째, 반투과성 막이 존재함으로써 외부에서 유입 가능한 분자와 내부에 고정된 분자의 조합이 다양해져, 동일한 초기 화학 규칙이라도 서로 다른 프로토셀 집단에서 상이한 ACS가 형성될 확률이 높아진다. 둘째, 작은 프로토셀 내부에서는 ‘경로 의존성’이 두드러져, 특정 촉매가 초기 단계에서 우연히 생성되면 그 촉매가 주도하는 반응망이 급속히 확대되어 자가촉매 집합이 안정화된다. 반대로, 중요한 촉매가 늦게 나타나면 네트워크는 다른 방향으로 전개되거나 결국 소멸에 이른다. 셋째, 화학 규칙을 바꾸어 촉매 가능 최소 길이를 조정하거나, 반응률을 변형하면 ACS의 출현 빈도와 지속 시간이 크게 변한다. 이는 실제 원시 생명 시스템이 특정 화학 환경에 의존했을 가능성을 시사한다.

또한, 저자들은 모델의 한계도 명시한다. 고정 부피와 즉각적인 막 투과 평형 가정은 장기적인 성장·분열 과정에서 발생할 수 있는 물질 축적, 막의 물리적 변형, 삼투압 차이 등을 무시한다. 무한히 높은 농도 성장 가능성도 비현실적이므로, 시뮬레이션 시간은 실제 물리적 제한에 도달하기 전까지로 제한된다. 향후 연구에서는 이러한 요소들을 단계적으로 도입해, 프로토셀의 성장·분열과 화학 네트워크의 피드백을 통합한 모델을 구축할 계획이다.

전반적으로, 이 논문은 ‘프로토셀 안에서 촉매 반응망이 어떻게 확률적·동적 특성을 보이는가’에 대한 최초의 정량적 프레임워크를 제공한다. 작은 구획화된 환경이 화학적 다양성과 진화 가능성을 증폭시킨다는 가설을 수치적으로 뒷받침하며, 원시 생명체의 초기 단계에서 물리·화학적 제약이 어떻게 선택적 압력으로 작용했는지를 이해하는 데 중요한 단서를 제공한다.