외부 상호작용을 포함한 RNA 매트릭스 모델과 그 비대칭 거동

초록

본 논문은 RNA 사슬의 일부 염기에 외부 교란을 가했을 때, 매트릭스 모델의 생성함수에 나타나는 ((1-\frac{n\alpha}{L})) 인자와 그에 따른 계통도(다양체) 분포의 비대칭적 성장법칙을 수치적으로 분석한다. 교란 비율 (n\alpha/L) 가 0에서 1로 증가함에 따라 고정 길이 (L)와 토포로지 차수 (g)에 대한 다이어그램 수 (a_{L,g})는 (3^{L})에서 ((3-\frac{n\alpha}{L})^{L}) 로, 최종적으로 (n\alpha/L=1) 일 때는 (2^{L}) 로 변한다. 전체 다이어그램 수 (\mathcal N) 역시 (\exp(\sqrt{L}))에서 (\exp\big((1-\frac{n\alpha}{L})\sqrt{L}\big)) 로 억제되며, 교란이 전체 사슬에 적용되면 (\mathcal N=1) 이 된다.

상세 분석

이 연구는 기존의 무교란 RNA 매트릭스 모델을 일반화하여, 외부 교란이 사슬의 특정 부분에만 작용하는 경우를 수학적으로 구현한다. 교란은 원래의 상호작용 강도와 비교한 비율 (\alpha) 로 정의되며, 교란이 작용하는 염기 수를 (n) 으로 두어 전체 사슬 길이 (L) 에 대한 비율 (\frac{n\alpha}{L}) 가 핵심 파라미터가 된다. 이 파라미터는 생성함수의 지수항에 ((1-\frac{n\alpha}{L})) 라는 가감 인자를 도입함으로써, 다이어그램 가중치가 교란 정도에 따라 선형적으로 감소하도록 만든다.

매트릭스 모델의 토포로지 전개는 ‘genus’ (g) 로 구분되는 계통도(다이어그램) 집합을 제공한다. 기존 무교란 경우, 고정 길이 (L) 에 대해 각 genus 별 다이어그램 수 (a_{L,g}) 는 대략 (3^{L}) 에 비례하는 보편성을 보였으며, 전체 다이어그램 수 (\mathcal N) 은 (\exp(\sqrt{L})) 라는 서브지수적 성장률을 갖는다. 교란 파라미터가 도입되면, 수치 실험을 통해 두 가지 극단 상황을 조사하였다. 첫 번째는 교란이 전체 사슬에 미치는 경우((n=L)), 두 번째는 단일 염기에만 작용하는 경우((n=1)).

결과는 (\frac{n\alpha}{L}) 가 0 에서 1 로 증가함에 따라, (a_{L,g}) 의 기본 베이스가 (3) 에서 ((3-\frac{n\alpha}{L})) 로 연속적으로 감소함을 보여준다. 특히 (\frac{n\alpha}{L}=1) 일 때는 베이스가 (2) 로 수렴하여, 이는 교란이 사슬 전체에 완전히 지배적인 상황에서 가능한 다이어그램의 조합이 크게 제한된다는 물리적 의미와 일치한다.

전체 다이어그램 수 (\mathcal N) 의 경우, 교란 인자는 지수적 성장인 (\exp(\sqrt{L})) 를 (\exp\big((1-\frac{n\alpha}{L})\sqrt{L}\big)) 로 변형시킨다. 따라서 교란이 강할수록 (\sqrt{L}) 스케일의 성장률이 억제되고, 완전 교란((\frac{n\alpha}{L}=1))에서는 (\exp(0)=1) 로 수렴한다. 이는 외부 힘이 RNA 구조 형성 자유도를 완전히 차단하는 극단적인 경우를 수학적으로 재현한다는 점에서 의미가 크다.

수치 분석은 주로 (L) 을 20~200 사이에서 변동시키며, 각 (L) 에 대해 다양한 (\frac{n\alpha}{L}) 값을 적용해 다이어그램 수를 직접 계산하였다. 로그-플롯을 이용한 선형 회귀는 베이스와 지수 인자가 위와 같이 변함을 명확히 확인시켰다. 또한, genus 분포는 교란이 증가할수록 고차 genus (복잡한 토포로지) 의 비중이 급격히 감소함을 보여, 교란이 구조 복잡성을 억제한다는 결론을 뒷받침한다.

이러한 결과는 RNA 2차 구조 예측에 외부 환경(예: 전기장, 화학적 변형) 이 미치는 영향을 정량적으로 이해하는 데 기여한다. 특히, 매트릭스 모델이 제공하는 토포로지 해석은 전통적인 동적 프로그래밍 기반 예측 알고리즘과는 다른 관점을 제공하며, 교란 파라미터를 조절함으로써 실험적 조건에 맞는 구조 집합을 사전 예측할 수 있는 가능성을 열어준다.