[论文解读] Maximum likelihood for high-noise group orbit estimation and single-particle cryo-EM

Zhou Fan, Roy R. Lederman+3 more|arXiv (Cornell University)|Jul 2, 2021

Advanced Electron Microscopy Techniques and Applications参考文献 36被引用 4

一句话总结

本文为单颗粒冷冻电镜中高噪声群体轨道估计提出了一种最大似然框架，通过按代数不变量分层的费舍尔信息特征值分析对数似然率景观。证明了三阶矩足以在 SO(2) 和 SO(3) 模型中局部识别出旋转对称下的信号，通过低维蛋白质势能图中的经验噪声缩放验证，确认了理论预测。

ABSTRACT

Motivated by applications to single-particle cryo-electron microscopy (cryo-EM), we study several problems of function estimation in a high noise regime, where samples are observed after random rotation and possible linear projection of the function domain. We describe a stratification of the Fisher information eigenvalues according to transcendence degrees of graded pieces of the algebra of group invariants, and we relate critical points of the log-likelihood landscape to a sequence of moment optimization problems, extending previous results for a discrete rotation group without projections. We then compute the transcendence degrees and forms of these optimization problems for several examples of function estimation under $SO(2)$ and $SO(3)$ rotations, including a simplified model of cryo-EM as introduced by Bandeira, Blum-Smith, Kileel, Perry, Weed, and Wein. We affirmatively resolve conjectures that $3^ ext{rd}$-order moments are sufficient to locally identify a generic signal up to its rotational orbit in these examples. For low-dimensional approximations of the electric potential maps of two small protein molecules, we empirically verify that the noise-scalings of the Fisher information eigenvalues conform with our theoretical predictions over a range of SNR, in a model of $SO(3)$ rotations without projections.

研究动机与目标

• 理解在高噪声群体轨道恢复中，特别是单颗粒冷冻电镜中的最大似然估计的统计与几何结构。
• 根据群不变量的超越次数，表征费舍尔信息矩阵，将统计估计与代数几何联系起来。
• 解决关于三阶矩在 SO(2) 和 SO(3) 旋转下是否足以实现局部信号识别的猜想。
• 通过低维蛋白质势能图上的数值模拟，验证理论预测的费舍尔信息特征值的噪声缩放行为。

提出的方法

• 根据 SO(2) 和 SO(3) 群作用下不变量代数的分次部分的超越次数，对费舍尔信息特征值进行分层。
• 将对数似然率景观的临界点重新表述为一系列矩优化问题。
• 计算 SO(2) 和 SO(3) 下函数估计的超越次数与矩形式，包括简化的冷冻电镜模型。
• 使用正交 Procrustes 对齐和连续多参考对齐作为分析工具，研究旋转下的估计问题。
• 通过数值模拟评估不同信噪比下费舍尔信息特征值，与理论预测进行比较。
• 应用傅里叶切片定理与断层扫描投影模型，将冷冻电镜成像与从二维投影重建三维体积联系起来。

实验结果

研究问题

• RQ1 在 SO(2) 和 SO(3) 群轨道模型中，三阶矩是否足以局部识别出一般信号的旋转轨道？
• RQ2 费舍尔信息矩阵的结构如何与旋转群的代数不变量相关联？
• RQ3 在信噪比变化的冷冻电镜类模型中，费舍尔信息特征值的噪声缩放行为如何？
• RQ4 理论预测的费舍尔信息特征值与低维蛋白质势能图中的经验观测结果相比如何？
• RQ5 在高噪声、旋转对称的估计模型中，对数似然率景观能否通过矩优化问题来表征？

主要发现

• 本文证实，在 SO(2) 和 SO(3) 模型中，三阶矩足以局部识别出一般信号的旋转轨道，解决了先前工作中提出的猜想。
• 费舍尔信息特征值的噪声缩放行为与基于群不变量超越次数推导出的理论预测一致，在数值模拟的多种信噪比范围内得到验证。
• 根据分次不变量的超越次数对费舍尔信息特征值进行分层，为理解旋转对称模型中估计难度提供了一种新颖的代数几何框架。
• 对于两个小型蛋白质分子的低维近似模型，费舍尔信息特征值的噪声缩放行为与理论预测相符，证实了该模型在实际场景中的准确性。
• 证明了对数似然函数的临界点对应于一系列矩优化问题的解，将先前仅适用于离散群的结果扩展至连续群情形。
• 理论与数值结果支持在冷冻电镜的最大似然估计中使用高阶矩，即使在二阶方法失效的高噪声条件下亦成立。