Orientation Estimation of Cryo-EM Images for Molecules with Arbitrary Symmetry by Semidefinite Programming
Organizations: Department of Mathematics, University of California, Davis · Department of Mathematics and Program in Applied and Computational Mathematics, Princeton University · Department of Statistics, University of Chicago · School of Mathematical Sciences, MOE-LSC, Institute of Natural Sciences, Shanghai Jiao Tong University
Abstract
This paper extends the non-unique games framework introduced in Bandeira et al. (2020) to address symmetric structures. Precisely, we consider the problem
where is a compact group, is a finite subgroup of symmetry elements, and is a loss function. A notable example is common lines based orientation estimation of cryo-EM images for molecules with point group symmetry , where is the rotation group . With a semidefinite programming relaxation, we turn the problem into a convex program that is solved with the alternating direction method of multipliers. Unlike existing approaches which are tailored to specific symmetry groups, ours provides a unified computational framework for addressing all point group symmetries. We demonstrate the applicability of the proposed algorithm on the cryo-EM orientation estimation problem on both simulated and experimental datasets. In particular, our algorithm is effective in generating high-quality ab initio models from only a few class averages computed from experimental images.
Figures & tables
| SNR | Clean | 1 | 1/2 | 1/4 | 1/8 |
| (without PR) | 0.030 (18.10) | 0.041 (30.53) | 0.034 (33.50) | 0.088 (44.10) | 0.570 (62.08) |
| Synchronization | 0.000 (7.58) | 0.005 (23.25) | 0.046 (33.27) | ||
| (without PR) | 0.032 (16.60) | 0.026 (25.61) | 0.047 (32.12) | 0.195 (49.83) | |
| (with PR) | 0.017 (12.73) | 0.017 (23.67) | 0.027 (29.58) | 0.043 (37.89) | |
| Synchronization | 0.521 (8.29) | 0.684 (40.00) | |||
| (without PR) | 0.026 (14.53) | 0.033 (24.38) | 0.096 (38.16) | 0.852 (60.91) |
| SNR | Clean | 1 | 1/2 | 1/4 | 1/8 |
| (without PR) | 0.022 (14.28) | 0.028 (21.76) | 0.039 (26.59) | 0.076 (34.91) | 0.398 (46.79) |
| Synchronization | 0.000 (7.54) | 0.003 (18.53) | 0.021 (24.55) | 0.150 (39.25) | |
| (without PR) | 0.019 (11.06) | 0.021 (20.58) | 0.043 (25.22) | 0.102 (35.73) | |
| (with PR) | 0.013 (8.92) | 0.017 (20.25) | 0.027 (22.75) | 0.035 (27.04) | |
| Synchronization | 0.003 (7.54) | 0.248 (24.35) | 0.661 (42.46) | ||
| (without PR) | 0.036 (14.64) | 0.032 (20.82) | 0.138 (33.92) | 0.568 (50.78) |
| SNR | Clean | 1 | 1/2 | 1/4 | 1/8 |
| (without PR) | 0.049 (35.52) | 0.035 (52.32) | 0.204 (61.37) | 0.684 (64.93) | |
| (with PR) | 0.040 (17.86) | 0.027 (51.54) | 0.123 (60.68) | 0.448 (60.07) | |
| (without PR) | 0.158 (60.60) | 0.220 (61.69) | 0.576 (43.23) | ||
| (with PR) | 0.082 (19.43) | 0.097 (50.67) | 0.511 (38.68) | ||
| (without PR) | 0.101 (12.52) | 0.180 (61.69) | |||
| (with PR) | 0.088 (12.15) | 0.167 (49.68) |
| SNR | Clean | 1 | 1/2 | 1/4 | 1/8 |
| (without PR) | 0.085 (12.72) | 0.064 (44.71) | 0.118 (52.61) | 0.903 (50.78) | |
| (with PR) | 0.082 (16.51) | 0.072 (28.61) | 0.096 (47.53) | 0.646 (58.90) | |
| (without PR) | 0.029 (14.66) | 0.160 (27.36) | 0.455 (62.33) | ||
| (with PR) | 0.024 (16.71) | 0.184 (26.25) | 0.456 (62.08) | ||
| (without PR) | 0.110 (12.65) | 0.185 (45.86) | 0.808 (53.78) | ||
| (with PR) | 0.109 (16.95) | 0.125 (44.83) | 0.771 (42.61) |
| 20 | 50 | 100 | 200 | 300 | |
|---|---|---|---|---|---|
| Time (minutes) | 2.0 | 6.8 | 21.3 | 78.7 | 176.2 |
Cryo-Electron Microscopy
Cryo-electron microscopy reconstructs biological structures from images of frozen specimens, from cellular volumes to atomic-resolution molecular models.
See all 27 papers