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.

Public Health Relevance

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.

Agency
National Institute of Health (NIH)
Institute
National Library of Medicine (NLM)
Type
Research Project (R01)
Project #
5R01LM010142-03
Application #
8318254
Study Section
Biomedical Library and Informatics Review Committee (BLR)
Program Officer
Ye, Jane
Project Start
2010-07-15
Project End
2014-07-14
Budget Start
2012-07-15
Budget End
2013-07-14
Support Year
3
Fiscal Year
2012
Total Cost
$389,256
Indirect Cost
$154,056
Name
Yale University
Department
Radiation-Diagnostic/Oncology
Type
Schools of Medicine
DUNS #
043207562
City
New Haven
State
CT
Country
United States
Zip Code
06520
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