This research proposes to develop fast expectation-maximization (EM) algorithms for 3D reconstruction of high resolution reconstruction of small particles in Cryogenic Electron Microscopy (Cryo-EM). The proposed algorithm incorporates the contrast transfer function (CTF) of the microscope and can be used with any basis functions. Further, two strategies are proposed to speed up the algorithm by a factor of 1000. The two strategies are (1) domain reduction where the E step of the algorithm is sped up, and (2) implementation on commodity graphics cards, where the massive parallelism of the graphics processors speeds up the entire algorithm. Preliminary results using simulations based on ribosome and the IP3R Ca+ channel protein structure show that the E and M steps of the proposed algorithm are fast and accurate. The proposed research extends previous work by (1) further developing the domain reduction strategy into a multi-resolution strategy, (2) systematically parallelizing the M step, and (3) implementing and validating the complete algorithm. Further, systematic software design is planned for implementing the algorithm on graphics processor units. This software will be made available on the world-wide web for other researchers. The algorithm will be validated using simulations and real-world Cryo-EM images. The reconstruction accuracy will be measured by the average mean squared error in reconstruction, and the resolution will be measured by Fourier shell correlation. The accuracy will also be compared to the reconstruction accuracy of conventional algorithms. The speed of the algorithm will be measured as a function of number of images and signal to noise ratio.
Discerning the structure (detailed shape) of proteins is essential to the understanding of their function, and to the development of therapeutic interventions including drugs that affect their function. Cryo-EM is a powerful and flexible method for the determination of structures. We seek to increase the speed and usefulness of this method through the development of a fast computer algorithm.
|Huang, Chenxi; Tagare, Hemant D (2016) Robust w-Estimators for Cryo-EM Class Means. IEEE Trans Image Process 25:893-906|
|Dvornek, Nicha C; Sigworth, Fred J; Tagare, Hemant D (2015) SubspaceEM: A fast maximum-a-posteriori algorithm for cryo-EM single particle reconstruction. J Struct Biol 190:200-14|
|Kucukelbir, Alp; Sigworth, Fred J; Tagare, Hemant D (2014) Quantifying the local resolution of cryo-EM density maps. Nat Methods 11:63-5|
|Huang, Chenxi; Tagare, Hemant D (2014) Robust estimation for class averaging in cryo-EM Single Particle Reconstruction. Conf Proc IEEE Eng Med Biol Soc 2014:3329-32|
|Kucukelbir, Alp; Sigworth, Fred J; Tagare, Hemant D (2012) A Bayesian adaptive basis algorithm for single particle reconstruction. J Struct Biol 179:56-67|
|Barthel, Andrew C; Tagare, Hemant; Sigworth, Fred J (2011) Surface-Constrained 3D Reconstruction in Cryo-EM. Conf Rec Asilomar Conf Signals Syst Comput :1026-1030|