Autocorrelation analysis for cryo-EM with sparsity constraints: Improved sample complexity and projection-based algorithms
arXiv:2209.10531 · doi:10.1073/pnas.2216507120
Abstract
The number of noisy images required for molecular reconstruction in single-particle cryo-electron microscopy (cryo-EM) is governed by the autocorrelations of the observed, randomly-oriented, noisy projection images. In this work, we consider the effect of imposing sparsity priors on the molecule. We use techniques from signal processing, optimization, and applied algebraic geometry to obtain new theoretical and computational contributions for this challenging non-linear inverse problem with sparsity constraints. We prove that molecular structures modeled as sums of Gaussians are uniquely determined by the second-order autocorrelation of their projection images, implying that the sample complexity is proportional to the square of the variance of the noise. This theory improves upon the non-sparse case, where the third-order autocorrelation is required for uniformly-oriented particle images and the sample complexity scales with the cube of the noise variance. Furthermore, we build a computational framework to reconstruct molecular structures which are sparse in the wavelet basis. This method combines the sparse representation for the molecule with projection-based techniques used for phase retrieval in X-ray crystallography.
arXiv comment31 pages, 5 figures, 1 movie
References in corpus (13)
- Relaxed Averaged Alternating Reflections for Diffraction Imaging
- Deep learning based mixed-dimensional GMM for characterizing variability in CryoEM
- Single-particle cryo-electron microscopy: Mathematical theory, computational challenges, and opportunities
- Guaranteed Blind Sparse Spikes Deconvolution via Lifting and Convex Optimization
- Reconstructing an Icosahedral Virus from Single-Particle Diffraction Experiments
- Method of moments for 3-D single particle ab initio modeling with non-uniform distribution of viewing angles
- Diagonal and Low-Rank Matrix Decompositions, Correlation Matrices, and Ellipsoid Fitting
- Common lines modeling for reference free ab-initio reconstruction in cryo-EM
- On random tomography with unobservable projection angles
- Solution of the phase problem for coherent scattering from a disordered system of identical particles
- Multi-target Detection with an Arbitrary Spacing Distribution
- Orthogonal Matrix Retrieval in Cryo-Electron Microscopy
- Sparse approximations of protein structure from noisy random projections