Coifman Ronald R, Shkolnisky Yoel, Sigworth Fred J, Singer Amit
Department of Mathematics, Program in Applied Mathematics, Yale University, 10 Hillhouse Ave. PO Box 208283, New Haven, CT 06520-8283 USA.
Appl Comput Harmon Anal. 2010 May;28(3):296-312. doi: 10.1016/j.acha.2009.11.003.
Recovering the three-dimensional structure of molecules is important for understanding their functionality. We describe a spectral graph algorithm for reconstructing the three-dimensional structure of molecules from their cryo-electron microscopy images taken at random unknown orientations.We first identify a one-to-one correspondence between radial lines in three-dimensional Fourier space of the molecule and points on the unit sphere. The problem is then reduced to determining the coordinates of points on the sphere given a subset of their pairwise geodesic distances. To recover those coordinates, we exploit the special geometry of the problem, as rendered by the Fourier projection-slice theorem, to construct a weighted graph whose vertices are the radial Fourier lines and whose edges are linked using the common line property. The graph organizes the radial lines on the sphere in a global manner that reveals the acquisition direction of each image. This organization is derived from a global computation of a few eigenvectors of the graph's sparse adjacency matrix. Once the directions are obtained, the molecule can be reconstructed using classical tomography methods.The presented algorithm is direct (as opposed to iterative refinement schemes), does not require any prior model for the reconstructed object, and is shown to have favorable computational and numerical properties. Moreover, the algorithm does not impose any assumption on the distribution of the projection orientations. Physically, this means that the algorithm is applicable to molecules that have unknown spatial preference.
恢复分子的三维结构对于理解其功能很重要。我们描述了一种光谱图算法,用于从以随机未知方向拍摄的冷冻电子显微镜图像中重建分子的三维结构。我们首先确定分子三维傅里叶空间中的径向线与单位球面上的点之间的一一对应关系。然后,问题就简化为在给定其成对测地距离子集的情况下确定球面上点的坐标。为了恢复这些坐标,我们利用傅里叶投影切片定理所呈现的问题的特殊几何结构,构建一个加权图,其顶点是径向傅里叶线,其边使用共线性质相连。该图以全局方式组织球面上的径向线,揭示每个图像的采集方向。这种组织源于对图的稀疏邻接矩阵的几个特征向量的全局计算。一旦获得方向,就可以使用经典断层扫描方法重建分子。所提出的算法是直接的(与迭代细化方案相对),不需要对重建对象有任何先验模型,并且显示出具有良好的计算和数值特性。此外,该算法对投影方向的分布不做任何假设。从物理意义上讲,这意味着该算法适用于具有未知空间偏好的分子。