Advances in Xmipp for Cryo-Electron Microscopy: From Xmipp to Scipion.

Xmipp is an open-source software package consisting of multiple programs for processing data originating from electron microscopy and electron tomography, designed and managed by the Biocomputing Unit of the Spanish National Center for Biotechnology, although with contributions from many other developers over the world. During its 25 years of existence, Xmipp underwent multiple changes and updates. While there were many publications related to new programs and functionality added to Xmipp, there is no single publication on the Xmipp as a package since 2013. In this article, we give an overview of the changes and new work since 2013, describe technologies and techniques used during the development, and take a peek at the future of the package.

Experimental estimation of the longitudinal component of a highly focused electromagnetic field.

The detection of the longitudinal component of a highly focused electromagnetic beam is not a simple task. Although in recent years several methods have been reported in the literature, this measure is still not routinely performed. This paper describes a method that allows us to estimate and visualize the longitudinal component of the field in a relatively simple way. First, we measure the transverse components of the focused field in several planes normal to the optical axis. Then, we determine the complex amplitude of the two transverse field components: the phase is obtained using a phase recovery algorithm, while the phase difference between the two components is determined from the Stokes parameters. Finally, the longitudinal component is estimated using the Gauss's theorem. Experimental results show an excellent agreement with theoretical predictions.

Image Processing in Cryo-Electron Microscopy of Single Particles: The Power of Combining Methods.

Cryo-electron microscopy has established as a mature structural biology technique to elucidate the three-dimensional structure of biological macromolecules. The Coulomb potential of the sample is imaged by an electron beam, and fast semi-conductor detectors produce movies of the sample under study. These movies have to be further processed by a whole pipeline of image-processing algorithms that produce the final structure of the macromolecule. In this chapter, we illustrate this whole processing pipeline putting in value the strength of "meta algorithms," which are the combination of several algorithms, each one with different mathematical rationale, in order to distinguish correctly from incorrectly estimated parameters. We show how this strategy leads to superior performance of the whole pipeline as well as more confident assessments about the reconstructed structures. The "meta algorithms" strategy is common to many fields and, in particular, it has provided excellent results in bioinformatics. We illustrate this combination using the workflow engine, Scipion.

FSC-Q: a CryoEM map-to-atomic model quality validation based on the local Fourier shell correlation.

In recent years, advances in cryoEM have dramatically increased the resolution of reconstructions and, with it, the number of solved atomic models. It is widely accepted that the quality of cryoEM maps varies locally; therefore, the evaluation of the maps-derived structural models must be done locally as well. In this article, a method for the local analysis of the map-to-model fit is presented. The algorithm uses a comparison of two local resolution maps. The first is the local FSC (Fourier shell correlation) between the full map and the model, while the second is calculated between the half maps normally used in typical single particle analysis workflows. We call the quality measure "FSC-Q", and it is a quantitative estimation of how much of the model is supported by the signal content of the map. Furthermore, we show that FSC-Q may be helpful to detect overfitting. It can be used to complement other methods, such as the Q-score method that estimates the resolvability of atoms.

MicrographCleaner: A python package for cryo-EM micrograph cleaning using deep learning.

Cryo-EM Single Particle Analysis workflows require tens of thousands of high-quality particle projections to unveil the three-dimensional structure of macromolecules. Conventional methods for automatic particle picking tend to suffer from high false-positive rates, hampering the reconstruction process. One common cause of this problem is the presence of carbon and different types of high-contrast contaminations. In order to overcome this limitation, we have developed MicrographCleaner, a deep learning package designed to discriminate, in an automated fashion, between regions of micrographs which are suitable for particle picking, and those which are not. MicrographCleaner implements a U-net-like deep learning model trained on a manually curated dataset compiled from over five hundred micrographs. The benchmarking, carried out on approximately one hundred independent micrographs, shows that MicrographCleaner is a very efficient approach for micrograph preprocessing. MicrographCleaner (micrograph_cleaner_em) package is available at PyPI and Anaconda Cloud and also as a Scipion/Xmipp protocol. Source code is available at https://github.com/rsanchezgarc/micrograph_cleaner_em.

Automatic local resolution-based sharpening of cryo-EM maps.

MOTIVATION: Recent technological advances and computational developments have allowed the reconstruction of Cryo-Electron Microscopy (cryo-EM) maps at near-atomic resolution. On a typical workflow and once the cryo-EM map has been calculated, a sharpening process is usually performed to enhance map visualization, a step that has proven very important in the key task of structural modeling. However, sharpening approaches, in general, neglects the local quality of the map, which is clearly suboptimal. RESULTS: Here, a new method for local sharpening of cryo-EM density maps is proposed. The algorithm, named LocalDeblur, is based on a local resolution-guided Wiener restoration approach of the original map. The method is fully automatic and, from the user point of view, virtually parameter-free, without requiring either a starting model or introducing any additional structure factor correction or boosting. Results clearly show a significant impact on map interpretability, greatly helping modeling. In particular, this local sharpening approach is especially suitable for maps that present a broad resolution range, as is often the case for membrane proteins or macromolecules with high flexibility, all of them otherwise very suitable and interesting specimens for cryo-EM. To our knowledge, and leaving out the use of local filters, it represents the first application of local resolution in cryo-EM sharpening. AVAILABILITY AND IMPLEMENTATION: The source code (LocalDeblur) can be found at https://github.com/I2PC/xmipp and can be run using Scipion (http://scipion.cnb.csic.es) (release numbers greater than or equal 1.2.1). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Validation of electron microscopy initial models via small angle X-ray scattering curves.

MOTIVATION: Cryo electron microscopy (EM) is currently one of the main tools to reveal the structural information of biological macromolecules. The re-construction of three-dimensional (3D) maps is typically carried out following an iterative process that requires an initial estimation of the 3D map to be refined in subsequent steps. Therefore, its determination is key in the quality of the final results, and there are cases in which it is still an open issue in single particle analysis (SPA). Small angle X-ray scattering (SAXS) is a well-known technique applied to structural biology. It is useful from small nanostructures up to macromolecular ensembles for its ability to obtain low resolution information of the biological sample measuring its X-ray scattering curve. These curves, together with further analysis, are able to yield information on the sizes, shapes and structures of the analyzed particles. RESULTS: In this paper, we show how the low resolution structural information revealed by SAXS is very useful for the validation of EM initial 3D models in SPA, helping the following refinement process to obtain more accurate 3D structures. For this purpose, we approximate the initial map by pseudo-atoms and predict the SAXS curve expected for this pseudo-atomic structure. The match between the predicted and experimental SAXS curves is considered as a good sign of the correctness of the EM initial map. AVAILABILITY AND IMPLEMENTATION: The algorithm is freely available as part of the Scipion 1.2 software at http://scipion.i2pc.es/.

Deep Consensus, a deep learning-based approach for particle pruning in cryo-electron microscopy.

Single-particle cryo-electron microscopy (cryo-EM) has recently become a mainstream technique for the structural determination of macromolecules. Typical cryo-EM workflows collect hundreds of thousands of single-particle projections from thousands of micrographs using particle-picking algorithms. However, the number of false positives selected by these algorithms is large, so that a number of different 'cleaning steps' are necessary to decrease the false-positive ratio. Most commonly employed techniques for the pruning of false-positive particles are time-consuming and require user intervention. In order to overcome these limitations, a deep learning-based algorithm named Deep Consensus is presented in this work. Deep Consensus works by computing a smart consensus over the output of different particle-picking algorithms, resulting in a set of particles with a lower false-positive ratio than the initial set obtained by the pickers. Deep Consensus is based on a deep convolutional neural network that is trained on a semi-automatically generated data set. The performance of Deep Consensus has been assessed on two well known experimental data sets, virtually eliminating user intervention for pruning, and enhances the reproducibility and objectivity of the whole process while achieving precision and recall figures above 90%.

Effect of linear polarizers on the behavior of partially coherent and partially polarized highly focused fields.

In this Letter, we describe the behavior of partially coherent, partially polarized focused vector beams after passing a linear polarizer placed at the focal plane of a high numerical aperture microscope lens. In particular, we develop a mathematical framework for such beams that helps the understanding of the performance of polarizers when interact with non-paraxial beams. The features of the focused field after the polarizer are numerically evaluated for some illustrative examples.

Synthesis of light needles with tunable length and nearly constant irradiance.

We introduce a new method for producing optical needles with tunable length and almost constant irradiance based on the evaluation of the on-axis power content of the light distribution at the focal area. According to theoretical considerations, we propose an adaptive modulating continuous function that presents a large derivative and a zero value jump at the entrance pupil of the focusing system. This distribution is displayed on liquid crystal devices using holographic techniques. In this way, a polarized input beam is shaped and subsequently focused using a high numerical aperture (NA) objective lens. As a result, needles with variable length and nearly constant irradiance are produced using conventional optics components. This procedure is experimentally demonstrated obtaining a 53λ-long and 0.8λ-wide needle.

Polarisers in the focal domain: Theoretical model and experimental validation.

Polarisers are one of the most widely used devices in optical set-ups. They are commonly used with paraxial beams that propagate in the normal direction of the polariser plane. Nevertheless, the conventional projection character of these devices may change when the beam impinges a polariser with a certain angle of incidence. This effect is more noticeable if polarisers are used in optical systems with a high numerical aperture, because multiple angles of incidence have to be taken into account. Moreover, the non-transverse character of highly focused beams makes the problem more complex and strictly speaking, the Malus' law does not apply. In this paper we develop a theoretical framework to explain how ideal polarisers affect the behavior of highly focused fields. In this model, the polarisers are considered as birefringent plates, and the vector behaviour of focused fields is described using the plane-wave angular spectrum approach. Experiments involving focused fields were conducted to verify the theoretical model and a satisfactory agreement between theoretical and experimental results was found.

Optical encryption using photon-counting polarimetric imaging.

We present a polarimetric-based optical encoder for image encryption and verification. A system for generating random polarized vector keys based on a Mach-Zehnder configuration combined with translucent liquid crystal displays in each path of the interferometer is developed. Polarization information of the encrypted signal is retrieved by taking advantage of the information provided by the Stokes parameters. Moreover, photon-counting model is used in the encryption process which provides data sparseness and nonlinear transformation to enhance security. An authorized user with access to the polarization keys and the optical design variables can retrieve and validate the photon-counting plain-text. Optical experimental results demonstrate the feasibility of the encryption method.

Synthesis of highly focused fields with circular polarization at any transverse plane.

We develop a method for generating focused vector beams with circular polarization at any transverse plane. Based on the Richards-Wolf vector model, we derive analytical expressions to describe the propagation of these set of beams near the focal area. Since the polarization and the amplitude of the input beam are not uniform, an interferometric system capable of generating spatially-variant polarized beams has to be used. In particular, this wavefront is manipulated by means of spatial light modulators displaying computer generated holograms and subsequently focused using a high numerical aperture objective lens. Experimental results using a NA = 0.85 system are provided: irradiance and Stokes images of the focused field at different planes near the focal plane are presented and compared with those obtained by numerical simulation.

Experimental implementation of tightly focused beams with unpolarized transversal component at any plane.

The aim of this paper is to provide a formal framework for designing highly focused fields with specific transversal features when the incoming beam is partially polarized. More specifically, we develop a field with a transversal component that remains unpolarized in the focal area. Special attention is paid to the design of the input beam and the development of the experiment. The implementation of such fields is possible by using an interferometric setup combined with the use of digital holography techniques. Experimental results are compared with those obtained numerically.

Reconfigurable beams with arbitrary polarization and shape distributions at a given plane.

Methods for generating beams with arbitrary polarization based on the use of liquid crystal displays have recently attracted interest from a wide range of sources. In this paper we present a technique for generating beams with arbitrary polarization and shape distributions at a given plane using a Mach-Zehnder setup. The transverse components of the incident beam are processed independently by means of spatial light modulators placed in each path of the interferometer. The modulators display computer generated holograms designed to dynamically encode any amplitude value and polarization state for each point of the wavefront in a given plane. The steps required to design such beams are described in detail. Several beams performing different polarization and intensity landscapes have been experimentally implemented. The results obtained demonstrate the capability of the proposed technique to tailor the amplitude and polarization of the beam simultaneously.