Cryo-EM image alignment based on nonuniform fast Fourier transform.

In single particle analysis, two-dimensional (2-D) alignment is a fundamental step intended to put into register various particle projections of biological macromolecules collected at the electron microscope. The efficiency and quality of three-dimensional (3-D) structure reconstruction largely depends on the computational speed and alignment accuracy of this crucial step. In order to improve the performance of alignment, we introduce a new method that takes advantage of the highly accurate interpolation scheme based on the gridding method, a version of the nonuniform fast Fourier transform, and utilizes a multi-dimensional optimization algorithm for the refinement of the orientation parameters. Using simulated data, we demonstrate that by using less than half of the sample points and taking twice the runtime, our new 2-D alignment method achieves dramatically better alignment accuracy than that based on quadratic interpolation. We also apply our method to image to volume registration, the key step in the single particle EM structure refinement protocol. We find that in this case the accuracy of the method not only surpasses the accuracy of the commonly used real-space implementation, but results are achieved in much shorter time, making gridding-based alignment a perfect candidate for efficient structure determination in single particle analysis.

[1]  Pawel A Penczek,et al.  Gridding-based direct Fourier inversion of the three-dimensional ray transform. , 2004, Journal of the Optical Society of America. A, Optics, image science, and vision.

[2]  Jorge Nocedal,et al.  Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization , 1997, TOMS.

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

[4]  Pawel A Penczek,et al.  The Transform Class in SPARX and EMAN2. , 2007, Journal of structural biology.

[5]  M. J. D. Powell,et al.  On search directions for minimization algorithms , 1973, Math. Program..

[6]  J M Carazo,et al.  A multiresolution approach to orientation assignment in 3D electron microscopy of single particles. , 2004, Journal of structural biology.

[7]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[8]  G. Beylkin On the Fast Fourier Transform of Functions with Singularities , 1995 .

[9]  S Jonić,et al.  Spline-based image-to-volume registration for three-dimensional electron microscopy. , 2005, Ultramicroscopy.

[10]  Frank Natterer,et al.  Mathematical methods in image reconstruction , 2001, SIAM monographs on mathematical modeling and computation.

[11]  Hermann Schomberg Notes on direct and gridding-based Fourier inversion methods , 2002, Proceedings IEEE International Symposium on Biomedical Imaging.

[12]  H. Stark,et al.  Corrim-based alignment for improved speed in single-particle image processing. , 2003, Journal of structural biology.

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

[14]  Jan Timmer,et al.  The gridding method for image reconstruction by Fourier transformation , 1995, IEEE Trans. Medical Imaging.

[15]  A Leith,et al.  SPIDER and WEB: processing and visualization of images in 3D electron microscopy and related fields. , 1996, Journal of structural biology.

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

[17]  Jeffrey A. Fessler,et al.  Nonuniform fast Fourier transforms using min-max interpolation , 2003, IEEE Trans. Signal Process..

[18]  Vladimir Rokhlin,et al.  Fast Fourier Transforms for Nonequispaced Data , 1993, SIAM J. Sci. Comput..

[19]  E. Saff,et al.  Distributing many points on a sphere , 1997 .

[20]  J Frank,et al.  Classification of macromolecular assemblies studied as ‘single particles’ , 1990, Quarterly Reviews of Biophysics.

[21]  J M Carazo,et al.  XMIPP: a new generation of an open-source image processing package for electron microscopy. , 2004, Journal of structural biology.

[22]  J. Frank,et al.  The ribosome at improved resolution: new techniques for merging and orientation refinement in 3D cryo-electron microscopy of biological particles. , 1994, Ultramicroscopy.

[23]  Eric Westhof,et al.  Structure of the ribosome-bound cricket paralysis virus IRES RNA , 2006, Nature Structural &Molecular Biology.

[24]  Nikolaus Grigorieff,et al.  FREALIGN: high-resolution refinement of single particle structures. , 2007, Journal of structural biology.

[25]  J. D. O'Sullivan,et al.  A Fast Sinc Function Gridding Algorithm for Fourier Inversion in Computer Tomography , 1985, IEEE Transactions on Medical Imaging.

[26]  Willy Wriggers,et al.  2D fast rotational matching for image processing of biophysical data. , 2003, Journal of structural biology.

[27]  S. Provencher,et al.  Three-dimensional reconstruction from electron micrographs of disordered specimens. I. Method. , 1988, Ultramicroscopy.

[28]  Volker Rasche,et al.  Resampling of data between arbitrary grids using convolution interpolation , 1999, IEEE Transactions on Medical Imaging.

[29]  S. Provencher,et al.  Three-dimensional reconstruction from electron micrographs of disordered specimens. II. Implementation and results. , 1988, Ultramicroscopy.

[30]  Thierry Blu,et al.  Linear interpolation revitalized , 2004, IEEE Transactions on Image Processing.

[31]  Chao Yang,et al.  SPARX, a new environment for Cryo-EM image processing. , 2007, Journal of structural biology.

[32]  Chao Yang,et al.  Unified 3-D structure and projection orientation refinement using quasi-Newton algorithm. , 2005, Journal of structural biology.

[33]  J. Nocedal Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[34]  W. O. Saxton,et al.  The correlation averaging of a regularly arranged bacterial cell envelope protein , 1982, Journal of microscopy.

[35]  Jeffrey A. Fessler,et al.  Iterative tomographic image reconstruction using Fourier-based forward and back-projectors , 2004, IEEE Transactions on Medical Imaging.