[논문 리뷰] Enumeration of RNA structures by Matrix Models

Graziano Vernizzi, Henri Orland외 1명|HAL (Le Centre pour la Communication Scientifique Directe)|2004. 10. 31.

RNA and protein synthesis mechanisms인용 수 6

한 줄 요약

이 논문은 RNA 접촉 구조를 그들의 위상적 종수(genus)에 따라 정확하게 세는 행렬 모델을 제안하며, 긴 사슬에서는 비평면적 가짜접힘(pseudoknotted) 구조가 지배함을 드러낸다. 큰 N 행렬 적분 설정을 통해 정확한 생성함수를 유도하여 종수에 따라 구조를 분류하였으며, 종수가 시퀀스 길이에 비례하여 선형적으로 증가함(⟨g⟩ ≈ 0.23L)을 보였다.

ABSTRACT

We enumerate the number of RNA contact structures according to their genus, i.e. the topological character of their pseudoknots. By using a recently proposed matrix model formulation for the RNA folding problem, we obtain exact results for the simple case of an RNA molecule with an infinitely flexible backbone, in which any arbitrary pair of bases is allowed. We analyze the distribution of the genus of pseudoknots as a function of the total number of nucleotides along the phosphate-sugar backbone.

연구 동기 및 목표

• 평면적 이차 구조를 초월하는 RNA 접촉 구조를 세는 조합적 프레임워크를 개발하기 위해.
• 행렬 모델 설정을 통해 RNA 가짜접힘을 위상적 종수에 따라 분류하기 위해.
• RNA 길이와 종수에 따른 가짜접힘 복잡성의 분포를 정량화하기 위해.
• 무작위 RNA 서열에서 비평면적 구조가 통계적으로 지배적임을 드러내기 위해.
• 고정된 길이와 종수에서의 구조 수에 대한 정확한 해석적 결과를 제공하여, 엔트로피 및 위상공간 분석을 가능하게 하기 위해.

제안 방법

• 탄성적인 RNA 사슬의 분할 함수를 헤르미트 행렬 위의 행렬 적분으로 설정하고, 큰 N 전개를 사용한다.
• 가우시안 적분과 에르미트 다항식 항등식을 통해 생성함수 G(t,N)를 유도하며, 이는 종수에 따라 구조 수를 코딩한다.
• 생성함수를 t(염기 수)와 1/N²(종수)의 멱급수로 전개하고, 계수 a_{L,g}는 길이 L과 종수 g에서의 정확한 구조 수를 제공한다.
• 일반화된 라게르 다항식을 사용하여 정확한 해를 표현하고, 이는 구조 수의 재귀적 계산을 가능하게 한다.
• 모든 가능한 구조에 걸쳐 종수 분포를 계산하기 위해 총 다이어그램 수로 수를 정규화한다.
• 수치적으로 渐近적 행동을 추출하여, 종수 ⟨g⟩가 길이 L에 비례하여 선형적으로 증가함(⟨g⟩ ≈ 0.23L)을 보였다.

실험 결과

연구 질문

• RQ1 주어진 염기 수와 위상적 종수를 가진 RNA 접촉 구조의 정확한 수는 얼마인가?
• RQ2 RNA 길이가 증가함에 따라 가짜접힘 종수의 분포는 어떻게 변화하는가?
• RQ3 무작위 RNA 구조 중 평면적 구조와 비평면적 구조의 비율은 얼마이며, 이 비율은 길이에 따라 어떻게 변화하는가?
• RQ4 RNA 구조의 평균 종수는 시퀀스 길이에 따라 어떻게 渐近적으로 스케일링되는가?
• RQ5 행렬 모델 설정이 평면적 이차 구조를 초월하여 RNA 접힘의 위상적 복잡성(종수)을 포괄할 수 있는가?

주요 결과

• 길이 L과 종수 g를 가진 RNA 구조의 수는 생성함수 G(t,N)에서 추출된 계수 a_{L,g}에 의해 정확히 주어진다.
• L=4일 때 종수-2 구조의 수는 2이며, L=6일 때 v=1일 경우 종수-1 구조의 총 수는 25이다.
• RNA 구조의 평균 종수 ⟨g⟩는 길이에 비례하여 선형적으로 증가하며, 큰 L에 대해 ⟨g⟩ ≈ 0.23L이다.
• 구조의 최대 가능한 종수는 g ≤ L/4이며, 이는 종수-g의 구조가 최소 L = 4g 염기를 필요로 함을 의미한다.
• 모든 L ≫ 1에 대해 대부분의 구조는 비평면적이며, 종수는 특성값 ⟨g⟩ ≈ 0.23L 근처에 집중되어 있다.
• 모델은 천연 RNA 서열이 고종수 가짜접힘을 피하기 위해 진화적으로 선택되었음을 보여주며, 대부분의 서열이 종수 g ≤ 2를 가진다.