진화적 복제·혁신·소실 모델의 평균장 이론

초록

이 논문은 단백질 도메인 클래스의 진화적 변화를 설명하기 위해 복제·혁신·소실 과정을 포함한 평균장 모델을 제시한다. 중국 레스토랑 과정(CRP)과 유사한 확률적 프레임워크를 구축하고, 시뮬레이션과 해석을 통해 도메인 클래스 규모와 수의 스케일링 법칙을 재현한다. 특히 요소 손실과 실제 도메인 클래스의 특수성을 모델에 반영함으로써 기존 모델과의 차이를 정량적으로 평가한다.

상세 분석

본 연구는 유전체 수준에서 관찰되는 단백질 도메인 클래스의 규모 분포가 멱법칙적 스케일링을 보인다는 실증적 사실에 착안한다. 이를 설명하기 위해 저자들은 ‘복제‑혁신‑소실’(duplication‑innovation‑loss) 과정을 기본 메커니즘으로 하는 확률적 모델을 설계한다. 모델의 핵심은 중국 레스토랑 과정(CRP)과 유사한 ‘중복’(복제) 규칙이다. 새로운 도메인이 등장할 확률은 현재 존재하는 클래스 수에 역비례하는 파라미터 α에 의해 조절되며, 기존 클래스에 새로운 원소가 추가되는 확률은 해당 클래스의 현재 크기에 비례한다. 이러한 기본 구조에 두 가지 중요한 변형을 도입한다. 첫째, ‘소실’ 메커니즘을 추가하여 각 클래스가 일정 확률로 원소를 잃거나 완전히 사라질 수 있게 한다. 이는 실제 유전체에서 도메인 손실이 빈번히 일어나는 현상을 반영한다. 둘째, 실제 도메인 클래스는 기능적·구조적 특성에 따라 서로 다른 유지·소실 확률을 가질 수 있다는 점을 모델에 포함한다. 이를 위해 클래스별 파라미터 β_i 를 도입해 클래스 특이성을 정량화한다.

평균장 접근법을 사용해 대규모 시스템의 동적 방정식을 근사한다. 저자들은 클래스 크기 분포 P(k)와 클래스 수 N_g(게놈 크기 g) 사이의 관계를 평균장 식으로 유도하고, 이를 통해 P(k)∝k^{-(1+α)} 형태의 멱법칙이 나타남을 보인다. 소실 파라미터 λ가 존재할 경우, 멱법칙의 지수는 α와 λ의 조합에 따라 미세하게 변형되며, 특히 큰 λ 값에서는 지수 감소가 관찰된다. 또한, 클래스 특이성을 반영한 경우에는 각 클래스별 지수가 다르게 나타나 전체 분포는 다중 멱법칙 또는 로그-정규 형태에 근접한다.

시뮬레이션 결과는 평균장 해석과 일치한다. 다양한 파라미터 조합에 대해 게놈 크기가 증가함에 따라 클래스 수가 g^{α} 스케일링을 따르고, 평균 클래스 크기가 g^{1‑α} 로 성장한다는 점이 확인된다. 특히, 소실을 포함한 모델은 실제 데이터에서 관찰되는 ‘희소 클래스’(크기가 1인 클래스)의 비율을 정확히 재현한다. 반면, 소실을 무시한 전통적 CRP 모델은 이러한 희소성을 과소평가한다.

결론적으로, 복제·혁신·소실을 동시에 고려한 평균장 모델은 단백질 도메인 레퍼토리의 스케일링 특성을 정량적으로 설명할 수 있다. 모델은 파라미터 α와 λ, 그리고 클래스 특이성 β_i 를 통해 다양한 진화적 시나리오를 탐색할 수 있는 유연성을 제공한다. 이는 유전체 진화 연구뿐 아니라, 복제와 혁신이 동시에 작용하는 다른 복합 시스템(예: 언어 어휘 성장, 기술 혁신 네트워크)에도 적용 가능한 일반적인 프레임워크를 제시한다.