Viewing Angle Classification of Cryo-Electron Microscopy Images Using Eigenvectors

The cryo-electron microscopy (cryo-EM) reconstruction problem is to find the three-dimensional structure of a macromolecule given noisy versions of its two-dimensional projection images at unknown random directions. We introduce a new algorithm for identifying noisy cryo-EM images of nearby viewing angles. This identification is an important first step in three-dimensional structure determination of macromolecules from cryo-EM, because once identified, these images can be rotationally aligned and averaged to produce "class averages" of better quality. The main advantage of our algorithm is its extreme robustness to noise. The algorithm is also very efficient in terms of running time and memory requirements, because it is based on the computation of the top few eigenvectors of a specially designed sparse Hermitian matrix. These advantages are demonstrated in numerous numerical experiments.

[1]  D. Goldberg,et al.  Assessing experimentally derived interactions in a small world , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Stéphane Lafon,et al.  Diffusion maps , 2006 .

[3]  Manfredo P. do Carmo,et al.  Differential geometry of curves and surfaces , 1976 .

[4]  J. Frank.  Three-Dimensional Electron Microscopy of Macromolecular Assemblies: Visualization of Biological Molecules in Their Native State , 1996 .

[5]  Wen Jiang,et al.  Fast rotational matching of single-particle images. , 2005, Journal of structural biology.

[6]  Laurent Joyeux,et al.  Efficiency of 2D alignment methods. , 2002, Ultramicroscopy.

[7]  M. Heel,et al.  Angular reconstitution: a posteriori assignment of projection directions for 3D reconstruction. , 1987 .

[8]  John Milnor,et al.  Analytic Proofs of the “Hairy Ball Theorem” and the Brouwer Fixed Point Theorem , 1978 .

[9]  Amit Singer,et al.  Computing Steerable Principal Components of a Large Set of Images and Their Rotations , 2011, IEEE Transactions on Image Processing.

[10]  Herbert Solomon,et al.  Geometric Probability , 1978, CBMS-NSF regional conference series in applied mathematics.

[11]  A. Khorunzhy.  Sparse random matrices: spectral edge and statistics of rooted trees , 2001, Advances in Applied Probability.

[12]  M. van Heel.  Angular reconstitution: a posteriori assignment of projection directions for 3D reconstruction. , 1987, Ultramicroscopy.

[13]  J. Frank,et al.  A common-lines based method for determining orientations for N > 3 particle projections simultaneously. , 1996, Ultramicroscopy.

[14]  E. Wigner.  On the Distribution of the Roots of Certain Symmetric Matrices , 1958 .

[15]  Duncan J. Watts,et al.  Collective dynamics of ‘small-world’ networks , 1998, Nature.

[16]  M. Radermacher,et al.  Three-dimensional reconstruction from random projections: orientational alignment via Radon transforms. , 1994, Ultramicroscopy.

[17]  Yoram Bresler,et al.  Feasibility of tomography with unknown view angles , 1998, Proceedings 1998 International Conference on Image Processing. ICIP98 (Cat. No.98CB36269).

[18]  Yoel Shkolnisky,et al.  Three-Dimensional Structure Determination from Common Lines in Cryo-EM by Eigenvectors and Semidefinite Programming , 2011, SIAM J. Imaging Sci..

[19]  Charles L. Epstein,et al.  Introduction to the mathematics of medical imaging , 2003 .

[20]  R. MacKinnon.  Potassium Channels and the Atomic Basis of Selective Ion Conduction , 2004, Angewandte Chemie.

[21]  Oleksiy Khorunzhiy.  Rooted trees and moments of large sparse random matrices , 2003, DRW.

[22]  Yoel Shkolnisky,et al.  Reference Free Structure Determination through Eigenvectors of Center of Mass Operators. , 2010, Applied and computational harmonic analysis.

[23]  R. MacKinnon.  Potassium channels and the atomic basis of selective ion conduction (Nobel Lecture). , 2004 .

[24]  E. Wigner.  Characteristic Vectors of Bordered Matrices with Infinite Dimensions I , 1955 .

[25]  B. Chait,et al.  The structure of the potassium channel: molecular basis of K+ conduction and selectivity. , 1998, Science.

[26]  M. Heel,et al.  Single-particle electron cryo-microscopy: towards atomic resolution , 2000, Quarterly Reviews of Biophysics.

[27]  Ronald R. Coifman,et al.  Graph Laplacian Tomography From Unknown Random Projections , 2008, IEEE Transactions on Image Processing.

[28]  A. Singer,et al.  Representation theoretic patterns in three dimensional Cryo-Electron Microscopy I: The intrinsic reconstitution algorithm. , 2009, Annals of mathematics.

[29]  A. Goncharov,et al.  Determination of the spatial orientation of arbitrarily arranged identical particles of unknown structure from their projections , 1986 .

[30]  F. Natterer.  The Mathematics of Computerized Tomography , 1986 .

[31]  A. Singer.  Angular Synchronization by Eigenvectors and Semidefinite Programming. , 2009, Applied and computational harmonic analysis.

[32]  Yoram Bresler,et al.  Uniqueness of tomography with unknown view angles , 2000, IEEE Trans. Image Process..

[33]  Amit Singer,et al.  Representation Theoretic Patterns in Three-Dimensional Cryo-Electron Microscopy II—The Class Averaging Problem , 2011, Found. Comput. Math..