[논문 리뷰] Tree Orbits under Permutation Group Action: Algorithm, Enumeration and Application to Viral Assembly

Miklós Bóna, Meera Sitharam외 1명|ArXiv.org|2009. 06. 01.

Bacteriophages and microbial interactions참고 문헌 23인용 수 3

한 줄 요약

이 논문은 이코사헤드랄 군의 작용 하에서 이코사헤드랄 바이러스 캡시드의 조립 트리들을 수세우고 분석하기 위한 조합적이고 알고리즘적인 프레임워크를 개발한다. 트리의 안정자(stabilizer)를 계산하는 선형 시간 알고리즘을 제시하고, 트리 궤도의 크기별로 닫힌 형식의 공식을 유도함으로써, 바이러스 자가조립 모델에서 조립 경로의 확률을 계산할 수 있다.

ABSTRACT

This paper uses combinatorics and group theory to answer questions about the assembly of icosahedral viral shells. Although the geometric structure of the capsid (shell) is fairly well understood in terms of its constituent subunits, the assembly process is not. For the purpose of this paper, the capsid is modeled by a polyhedron whose facets represent the monomers. The assembly process is modeled by a rooted tree, the leaves representing the facets of the polyhedron, the root representing the assembled polyhedron, and the internal vertices representing intermediate stages of assembly (subsets of facets). Besides its virological motivation, the enumeration of orbits of trees under the action of a finite group is of independent mathematical interest. If $G$ is a finite group acting on a finite set $X$, then there is a natural induced action of $G$ on the set $\mathcal{T}_X$ of trees whose leaves are bijectively labeled by the elements of $X$. If $G$ acts simply on $X$, then $|X| := |X_n| = n \cdot |G|$, where $n$ is the number of $G$-orbits in $X$. The basic combinatorial results in this paper are (1) a formula for the number of orbits of each size in the action of $G$ on $\mathcal{T}_{X_n}$, for every $n$, and (2) a simple algorithm to find the stabilizer of a tree $τ\in \mathcal{T}_X$ in $G$ that runs in linear time and does not need memory in addition to its input tree.

연구 동기 및 목표

• 조립 트리로 바이러스 캡시드 조립을 모델링하여, 단량체는 잎으로, 부분조립체는 내부 노드로 간주한다.
• 나무에 대한 군 작용을 통해 이코사헤드랄 대칭성 하에서 조립 경로의 확률 분포를 이해한다.
• 군 작용 하에서 주어진 트리의 안정자를 효율적으로 계산하는 알고리즘을 개발한다.
• 특히 60개의 면을 가진 T=1 캡시드에 대해 이코사헤드랄 군의 작용 하에서 트리 궤도를 크기별로 수세우는 것을 목표로 한다.
• 기하학적 및 대칭성 요소를 바탕으로 특정 바이러스 조립 경로의 가능성 확률을 추정하는 데 결과를 적용한다.

제안 방법

• 조립 트리로 바이러스 캡시드 조립을 모델링하며, 잎은 단량체(면), 내부 노드는 부분조립체, 루트는 완전한 캡시드로 간주한다.
• 이코사헤드랄 군 $G_{60}$을 레이블이 부여된 트리의 집합에 작용시켜 잎 레이블의 순열을 통해 트리 공간 $\mathcal{T}_X$에 작용을 유도한다.
• 군의 부분군 격자에 대한 모비우스 역전환을 사용하여 각 부분군 $G_i$에 의해 고정되는 트리의 수를 계산함으로써 궤도 수를 세는 데 사용한다.
• 부분군 $G_i$에 의해 고정되는 트리의 생성함수 $f_{G_i}(x)$를 유도하며, 이는 부분군 격자의 구조와 모비우스 함수 값에 기반한다.
• 추가 메모리가 필요 없이 순열군 이론에 기반한 선형 시간 알고리즘을 사용하여 주어진 트리 $\tau$의 안정자를 $G$에서 계산한다.
• 정리 2와 정리 3를 적용하여 정확히 $G_i$의 안정자를 가지는 트리의 수를 계산함으로써, 크기가 $i$인 궤도의 수를 산출한다.

실험 결과

연구 질문

• RQ1 T=1 캡시드에 대해 60개의 단량체를 가진 경우, 이코사헤드랄 군의 작용 하에서 서로 다른 조립 경로(트리 궤도)는 몇 개인가?
• RQ2 군 작용 하에서 궤도 크기에 따라 조립 경로의 확률 분포는 어떻게 되는가?
• RQ3 주어진 트리에 대해 이코사헤드랄 군 작용 하에서 안정자를 효율적으로 계산할 수 있는 알고리즘이 존재하는가?
• RQ4 기하학적 안정성과 대칭성 요소는 특정 조립 경로의 가능성에 어떻게 영향을 미치는가?
• RQ5 특정 이코사헤드랄 군의 부분군에 의해만 고정되는 트리의 조합적 구조는 어떠한가?

주요 결과

• T=1의 경우 $|X|=60$일 때, 정확히 $G_1$ (자명군)의 안정자를 가지는 조립 트리의 수는 $\overline{t}_{60}(G_1) = 1924465510132437394720184730922187571120346754532$ 이다.
• 완전한 이코사헤드랄 군 $G_{60}$의 안정자를 가지는 트리의 수는 $\overline{t}_1(G_{60}) = 204$ 이며, 이는 가장 대칭적인 경로에 해당한다.
• 인덱스 30인 $G_2$의 안정자를 가지는 트리의 수는 $\overline{t}_{30}(G_2) = 1670856367100496379411587456529324583988755126499875584$ 이다.
• 인덱스 20인 $G_3$의 안정자를 가지는 트리의 수는 $\overline{t}_{20}(G_3) = 10087157294451731428720995944759704$ 이다.
• 인덱스 12인 $G_5$의 안정자를 가지는 트리의 수는 $\overline{t}_{12}(G_5) = 20540071766413107840$ 이다.
• 오직 204개의 대칭적인 경로($G_{60}$-불변)가 존재하나, 비대칭적인 경로($G_1$-불변)의 총 수는 $10^{99}$를 초월하여 가능한 조립 경로의 광범위한 다양성을 보여준다.