ZART: A novel multiresolution reconstruction algorithm with motion-blur correction for single particle analysis

One of the main purposes of CryoEM Single Particle Analysis is to reconstruct the three-dimensional structure of a macromolecule thanks to the acquisition of many particle images representing different poses of the sample. By estimating the orientation of each projected particle, it is possible to recover the underlying 3D volume by multiple 3D reconstruction methods, usually working either in Fourier or in real space. However, the reconstruction from the projected images works under the assumption that all particles in the dataset correspond to the same conformation of the macromolecule. Although this requisite holds for some macromolecules, it is not true for flexible specimens, leading to motion-induced artefacts in the reconstructed CryoEM maps. In this work, we introduce a new Algebraic Reconstruction Technique called ZART, which is able to include continuous flexibility information during the reconstruction process to improve local resolution and reduce motion blurring. The conformational changes are modelled through Zernike3D polynomials. Our implementation allows for a multiresolution description of the macromolecule adapting itself to the local resolution of the reconstructed map. In addition, ZART has also proven to be a useful algorithm in cases where flexibility is not so dominant, as it improves the overall aspect of the reconstructed maps by improving their local and global resolution.

Estimating conformational landscapes from Cryo-EM particles by 3D Zernike polynomials

The new developments in Cryo-EM Single Particle Analysis are helping us to understand how the macromolecular structure and function meet to drive biological processes. By capturing many states at the particle level, it is possible to address how macromolecules explore different conformations, information that is classically extracted through 3D classification. However, the limitations of classical approaches prevent us from fully understanding the complete conformational landscape due to the reduced number of discrete states accurately reconstructed. To characterize the whole structural spectrum of a macromolecule, we propose an extension of our Zernike3D approach, able to extract per-image continuous flexibility information directly from a particle dataset. Also, our method can be seamlessly applied to images, maps or atomic models, opening integrative possibilities. Furthermore, we introduce the ZART reconstruction algorithm, which considers the Zernike3D deformation fields to revert particle conformational changes during the reconstruction process, thus minimizing the blurring induced by molecular motions.

On bias, variance, overfitting, gold standard and consensus in single-particle analysis by cryo-electron microscopy

Single-particle analysis (SPA) by cryo-electron microscopy comprises the estimation of many parameters along its image-processing pipeline. Overfitting observed in SPA is normally due to misestimated parameters, and the only way to identify these is by comparing the estimates of multiple algorithms or, at least, multiple executions of the same algorithm.

Cryo-electron microscopy (cryoEM) has become a well established technique to elucidate the 3D structures of biological macromolecules. Projection images from thousands of macromolecules that are assumed to be structurally identical are combined into a single 3D map representing the Coulomb potential of the macromolecule under study. This article discusses possible caveats along the image-processing path and how to avoid them to obtain a reliable 3D structure. Some of these problems are very well known in the community. These may be referred to as sample-related (such as specimen denaturation at interfaces or non-uniform projection geometry leading to underrepresented projection directions). The rest are related to the algorithms used. While some have been discussed in depth in the literature, such as the use of an incorrect initial volume, others have received much less attention. However, they are fundamental in any data-analysis approach. Chiefly among them, instabilities in estimating many of the key parameters that are required for a correct 3D reconstruction that occur all along the processing workflow are referred to, which may significantly affect the reliability of the whole process. In the field, the term overfitting has been coined to refer to some particular kinds of artifacts. It is argued that overfitting is a statistical bias in key parameter-estimation steps in the 3D reconstruction process, including intrinsic algorithmic bias. It is also shown that common tools (Fourier shell correlation) and strategies (gold standard) that are normally used to detect or prevent overfitting do not fully protect against it. Alternatively, it is proposed that detecting the bias that leads to overfitting is much easier when addressed at the level of parameter estimation, rather than detecting it once the particle images have been combined into a 3D map. Comparing the results from multiple algorithms (or at least, independent executions of the same algorithm) can detect parameter bias. These multiple executions could then be averaged to give a lower variance estimate of the underlying parameters.