V. Abrishami

Alignment of direct detection device micrographs using a robust Optical Flow approach

The introduction of direct detection devices in cryo-EM has shown that specimens present beam-induced motion (BIM). Consequently, in this work, we develop a BIM correction method at the image level, resulting in an integrated image in which the in-plane BIM blurring is compensated prior to particle picking. The methodology is based on a robust Optical Flow (OF) approach that can efficiently correct for local movements in a rapid manner. The OF works particularly well if the BIM pattern presents a substantial degree of local movements, which occurs in our data sets for Falcon II data. However, for those cases in which the BIM pattern corresponds to global movements, we have found it advantageous to first run a global motion correction approach and to subsequently apply OF. Additionally, spatial analysis of the Optical Flow allows for quantitative analysis of the BIM pattern. The software that incorporates the new approach is available in XMIPP (http://xmipp.cnb.csic.es).

Particle quality assessment and sorting for automatic and semiautomatic particle-picking techniques

Three-dimensional reconstruction of biological specimens using electron microscopy by single particle methodologies requires the identification and extraction of the imaged particles from the acquired micrographs. Automatic and semiautomatic particle selection approaches can localize these particles, minimizing the user interaction, but at the cost of selecting a non-negligible number of incorrect particles, which can corrupt the final three-dimensional reconstruction. In this work, we present a novel particle quality assessment and sorting method that can separate most erroneously picked particles from correct ones. The proposed method is based on multivariate statistical analysis of a particle set that has been picked previously using any automatic or manual approach. The new method uses different sets of particle descriptors, which are morphology-based, histogram-based and signal to noise analysis based. We have tested our proposed algorithm with experimental data obtaining very satisfactory results. The algorithm is freely available as a part of the Xmipp 3.0 package [http://xmipp.cnb.csic.es].

A fast iterative convolution weighting approach for gridding-based direct Fourier three-dimensional reconstruction with correction for the contrast transfer function.

We describe a fast and accurate method for the reconstruction of macromolecular complexes from a set of projections. Direct Fourier inversion (in which the Fourier Slice Theorem plays a central role) is a solution for dealing with this inverse problem. Unfortunately, the set of projections provides a non-equidistantly sampled version of the macromolecule Fourier transform in the single particle field (and, therefore, a direct Fourier inversion) may not be an optimal solution. In this paper, we introduce a gridding-based direct Fourier method for the three-dimensional reconstruction approach that uses a weighting technique to compute a uniform sampled Fourier transform. Moreover, the contrast transfer function of the microscope, which is a limiting factor in pursuing a high resolution reconstruction, is corrected by the algorithm. Parallelization of this algorithm, both on threads and on multiple CPUs, makes the process of three-dimensional reconstruction even faster. The experimental results show that our proposed gridding-based direct Fourier reconstruction is slightly more accurate than similar existing methods and presents a lower computational complexity both in terms of time and memory, thereby allowing its use on larger volumes. The algorithm is fully implemented in the open-source Xmipp package and is downloadable from http://xmipp.cnb.csic.es.

Cryo-EM and the elucidation of new macromolecular structures: Random Conical Tilt revisited

Cryo-Electron Microscopy (cryo-EM) of macromolecular complexes is a fundamental structural biology technique which is expanding at a very fast pace. Key to its success in elucidating the three-dimensional structure of a macromolecular complex, especially of small and non-symmetric ones, is the ability to start from a low resolution map, which is subsequently refined with the actual images collected at the microscope. There are several methods to produce this first structure. Among them, Random Conical Tilt (RCT) plays a prominent role due to its unbiased nature (it can create an initial model based on experimental measurements). In this article, we revise the fundamental mathematical expressions supporting RCT, providing new expressions handling all key geometrical parameters without the need of intermediate operations, leading to improved automation and overall reliability, essential for the success of cryo-EM when analyzing new complexes. We show that the here proposed RCT workflow based on the new formulation performs very well in practical cases, requiring very few image pairs (as low as 13 image pairs in one of our examples) to obtain relevant 3D maps.

Local analysis of strains and rotations for macromolecular electron microscopy maps

Macromolecular complexes perform their physiological functions by local rearrangements of their constituents and biochemically interacting with their reaction partners. These rearrangements may involve local rotations and the induction of local strains causing different mechanical efforts and stretches at the different areas of the protein. The analysis of these local deformations may reveal important insight into the way proteins perform their tasks. In this paper we introduce a method to perform this kind of local analysis using Electron Microscopy volumes in a fully objective and automatic manner. For doing so, we exploit the continuous nature of the result of an elastic image registration using B-splines as its basis functions. We show that the results obtained by the new automatic method are consistent with previous observations on these macromolecules.

An image processing approach to the simulation of electron microscopy volumes of atomic structures

The goal of this paper is to construct accurate three-dimensional models from atomic structures that can be compared to Electron Microscopy volumes. They, in their turn, can be used to solve fitting and docking problems as well as to generate ground-truth models for the development of image processing algorithms. We use Electron Atomic Scattering Factor (EASF) along with a careful filter design to greatly accelerate the current use of these functions and strongly reduce computer memory requirements. We show that our method is 3 orders of magnitude faster than existing implementations of EASFs in Fourier space as well as more accurate than their standard implementations. The algorithm is freely available in the software package Xmipp (http://xmipp.cnb.csic.es) [1].