[论文解读] Tree Orbits under Permutation Group Action: Algorithm, Enumeration and Application to Viral Assembly

Miklós Bóna, Meera Sitharam+1 more|ArXiv.org|Jun 1, 2009

Bacteriophages and microbial interactions参考文献 23被引用 3

一句话总结

本文提出了一种组合与算法框架，用于在二十面体群作用下枚举并分析二十面体病毒衣壳的组装树。该研究提出了一种线性时间算法来计算树的稳定子，并推导出按大小计数树轨道的闭式公式，从而实现病毒自组装模型中组装路径概率的计算。

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_{60}$ 的子群格上的Möbius反演，计算每个子群 $G_i$ 所固定的树的数量，从而实现轨道枚举。
• 基于子群格结构与Möbius函数值，推导出被子群 $G_i$ 固定的树的生成函数 $f_{G_i}(x)$。
• 采用线性时间算法计算给定树 $\tau$ 在 $G$ 中的稳定子，基于置换群理论，且无需额外内存。
• 应用定理2与定理3，计算稳定子恰好为 $G_i$ 的树的数量，从而得出大小为 $i$ 的轨道数量。

实验结果

研究问题

• RQ1 对于具有60个单体的T=1衣壳，在二十面体群作用下，存在多少条不同的组装路径（树轨道）？
• RQ2 在群作用下，给定轨道大小，组装路径的概率分布如何？
• RQ3 能否设计一种高效算法，计算在二十面体群作用下给定树的稳定子？
• RQ4 几何稳定性和对称性因素如何影响特定组装路径的可能性？
• RQ5 在二十面体群的特定子群下仅被固定不变的树的组合结构是什么？

主要发现

• 对于 $T=1$ 情况且 $|X|=60$，稳定子恰好为 $G_1$（平凡群）的组装树数量为 $\overline{t}_{60}(G_1) = 1924465510132437394720184730922187571120346754532$。
• 稳定子为 $G_{60}$（完整二十面体群）的树的数量为 $\overline{t}_1(G_{60}) = 204$，对应最对称的组装路径。
• 稳定子为 $G_2$（指数30）的树的数量为 $\overline{t}_{30}(G_2) = 1670856367100496379411587456529324583988755126499875584$。
• 稳定子为 $G_3$（指数20）的树的数量为 $\overline{t}_{20}(G_3) = 10087157294451731428720995944759704$。
• 稳定子为 $G_5$（指数12）的树的数量为 $\overline{t}_{12}(G_5) = 20540071766413107840$。
• 尽管仅有204条对称路径（$G_{60}$-不变），但总不对称路径（$G_1$-不变）数量超过 $10^{99}$，表明可能的组装路径种类极为丰富。