RESUMO
Image processing in cryogenic electron tomography (cryoET) is currently at a similar state as Single Particle Analysis (SPA) in cryogenic electron microscopy (cryoEM) was a few years ago. Its data processing workflows are far from being well defined and the user experience is still not smooth. Moreover, file formats of different software packages and their associated metadata are not standardized, mainly since different packages are developed by different groups, focusing on different steps of the data processing pipeline. The Scipion framework, originally developed for SPA (de la Rosa-Trevín et al., 2016), has a generic python workflow engine that gives it the versatility to be extended to other fields, as demonstrated for model building (Martínez et al., 2020). In this article, we provide an extension of Scipion based on a set of tomography plugins (referred to as ScipionTomo hereafter), with a similar purpose: to allow users to be focused on the data processing and analysis instead of having to deal with multiple software installation issues and the inconvenience of switching from one to another, converting metadata files, managing possible incompatibilities, scripting (writing a simple program in a language that the computer must convert to machine language each time the program is run), etcetera. Additionally, having all the software available in an integrated platform allows comparing the results of different algorithms trying to solve the same problem. In this way, the commonalities and differences between estimated parameters shed light on which results can be more trusted than others. ScipionTomo is developed by a collaborative multidisciplinary team composed of Scipion team engineers, structural biologists, and in some cases, the developers whose software packages have been integrated. It is open to anyone in the field willing to contribute to this project. The result is a framework extension that combines the acquired knowledge of Scipion developers in close collaboration with third-party developers, and the on-demand design of functionalities requested by beta testers applying this solution to actual biological problems.
Assuntos
Tomografia com Microscopia Eletrônica , Software , Algoritmos , Microscopia Crioeletrônica/métodos , Processamento de Imagem Assistida por Computador/métodos , Reprodutibilidade dos TestesRESUMO
Cryo-electron microscopy (cryo-EM) of frozen-hydrated preparations of isolated macromolecular complexes is the method of choice to obtain the structure of complexes that cannot be easily studied by other experimental methods due to their flexibility or large size. An increasing number of macromolecular structures are currently being obtained at subnanometer resolution but the interpretation of structural details in such EM-derived maps is often difficult because of noise at these high-frequency signal components that reduces their contrast. In this paper, we show that the method for EM density-map approximation using Gaussian functions can be used for denoising of single-particle EM maps of high (typically subnanometer) resolution. We show its denoising performance using simulated and experimental EM density maps of several complexes.
Assuntos
Microscopia Crioeletrônica/normas , Substâncias Macromoleculares/química , Razão Sinal-Ruído , Microscopia Crioeletrônica/métodos , Estrutura Molecular , Distribuição Normal , Conformação Proteica , Subunidades Ribossômicas Maiores/química , beta-Galactosidase/químicaRESUMO
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.
Assuntos
Substâncias Macromoleculares/química , Microscopia Eletrônica/métodos , Trifosfato de Adenosina/química , Algoritmos , Automação , Proteínas de Bactérias/química , Fenômenos Biomecânicos , Chaperonina 60/química , Proteínas de Choque Térmico/química , Humanos , Ribossomos Mitocondriais/química , Modelos Teóricos , Chaperonas Moleculares/química , Ligação Proteica , RotaçãoRESUMO
Electron microscopy (EM) provides access to structural information of macromolecular complexes in the 3-20 Å resolution range. Normal mode analysis has been extensively used with atomic resolution structures and successfully applied to EM structures. The major application of normal modes is the identification of possible conformational changes in proteins. The analysis can throw light on the mechanism following ligand binding, protein-protein interactions, channel opening and other functional macromolecular movements. In this article, we present a new web server, 3DEM Loupe, which allows normal mode analysis of any uploaded EM volume using a user-friendly interface and an intuitive workflow. Results can be fully explored in 3D through animations and movies generated by the server. The application is freely available at http://3demloupe.cnb.csic.es.
Assuntos
Substâncias Macromoleculares/ultraestrutura , Microscopia Eletrônica , Software , Internet , Substâncias Macromoleculares/química , Conformação ProteicaRESUMO
In this work we present a fast and automated algorithm for estimating the contrast transfer function (CTF) of a transmission electron microscope. The approach is very suitable for High Throughput work because: (a) it does not require any initial defocus estimation, (b) it is almost an order of magnitude faster than existing approaches, (c) it opens the way to well-defined extensions to the estimation of higher order aberrations, at the same time that provides defocus and astigmatism estimations comparable in accuracy to well established methods, such as Xmipp and CTFFIND3 approaches. The new algorithm is based on obtaining the wrapped modulating phase of the power spectra density pattern by the use of a quadrature filter. This phase is further unwrapped in order to obtain the continuous and smooth absolute phase map; then a Zernike polynomial fitting is performed and the defocus and astigmatism parameters are determined. While the method does not require an initial estimation of the defocus parameters or any non-linear optimization procedure, these approaches can be used if further refinement is desired. Results of the CTF estimation method are presented for standard negative stained images, cryo-electron microscopy images in the absence of carbon support, as well as micrographs with only ice. Additionally, we have also tested the proposed method with micrographs acquired from tilted and untilted samples, obtaining good results. The algorithm is freely available as a part of the Xmipp package [http://xmipp.cnb.csic.es].
Assuntos
Algoritmos , Ensaios de Triagem em Larga Escala/métodos , Processamento de Imagem Assistida por Computador/métodos , Microscopia Eletrônica de Transmissão/instrumentação , Microscopia Eletrônica de Transmissão/métodos , Modelos Teóricos , Calibragem , SoftwareRESUMO
Understanding how structure and function meet to drive biological processes is progressively shifting the cryoEM field towards a more advanced analysis of macromolecular flexibility. Thanks to techniques such as single-particle analysis and electron tomography, it is possible to image a macromolecule in different states, information that can subsequently be extracted through advanced image-processing methods to build a richer approximation of a conformational landscape. However, the interoperability of all of these algorithms remains a challenging task that is left to users, preventing them from defining a single flexible workflow in which conformational information can be addressed by different algorithms. Therefore, in this work, a new framework integrated into Scipion is proposed called the Flexibility Hub. This framework automatically handles intercommunication between different heterogeneity software, simplifying the task of combining the software into workflows in which the quality and the amount of information extracted from flexibility analysis is maximized.
Assuntos
Algoritmos , Software , Microscopia Crioeletrônica/métodos , Processamento de Imagem Assistida por Computador/métodos , Conformação Molecular , Substâncias Macromoleculares/químicaRESUMO
Single-particle analysis by electron microscopy is a well established technique for analyzing the three-dimensional structures of biological macromolecules. Besides its ability to produce high-resolution structures, it also provides insights into the dynamic behavior of the structures by elucidating their conformational variability. Here, the different image-processing methods currently available to study continuous conformational changes are reviewed.
Assuntos
Elétrons , Processamento de Imagem Assistida por Computador/estatística & dados numéricos , Imageamento Tridimensional/estatística & dados numéricos , Substâncias Macromoleculares/ultraestrutura , Microscopia Eletrônica/métodos , Proteínas/ultraestrutura , Algoritmos , Humanos , Substâncias Macromoleculares/química , Microscopia Eletrônica/instrumentação , Conformação Molecular , Simulação de Dinâmica Molecular , Análise de Componente Principal , Proteínas/química , TermodinâmicaRESUMO
Three-dimensional structure of a wide range of biological specimens can be computed from images collected by transmission electron microscopy. This information integrated with structural data obtained with other techniques (e.g., X-ray crystallography) helps structural biologists to understand the function of macromolecular complexes and organelles within cells. In this paper, we compare two three-dimensional transmission electron microscopy techniques that are becoming more and more related (at the image acquisition level as well as the image processing one): electron tomography and single-particle analysis. The first one is currently used to elucidate the three-dimensional structure of cellular components or smaller entire cells, whereas the second one has been traditionally applied to structural studies of macromolecules and macromolecular complexes. Also, we discuss possibilities for their integration with other structural biology techniques for an integrative study of living matter from proteins to whole cells.
Assuntos
Imageamento Tridimensional/métodos , Microscopia Eletrônica/métodos , Células/ultraestrutura , Substâncias MacromolecularesRESUMO
This paper presents an algorithm based on a continuous framework for a posteriori angular and translational assignment in three-dimensional electron microscopy (3DEM) of single particles. Our algorithm can be used advantageously to refine the assignment of standard quantized-parameter methods by registering the images to a reference 3D particle model. We achieve the registration by employing a gradient-based iterative minimization of a least-squares measure of dissimilarity between an image and a projection of the volume in the Fourier transform (FT) domain. We compute the FT of the projection using the central-slice theorem (CST). To compute the gradient accurately, we take advantage of a cubic B-spline model of the data in the frequency domain. To improve the robustness of the algorithm, we weight the cost function in the FT domain and apply a "mixed" strategy for the assignment based on the minimum value of the cost function at registration for several different initializations. We validate our algorithm in a fully controlled simulation environment. We show that the mixed strategy improves the assignment accuracy; on our data, the quality of the angular and translational assignment was better than 2 voxel (i.e., 6.54 angstroms). We also test the performance of our algorithm on real EM data. We conclude that our algorithm outperforms a standard projection-matching refinement in terms of both consistency of 3D reconstructions and speed.
Assuntos
Microscopia Eletrônica/métodos , Algoritmos , Bacteriorodopsinas/química , Bacteriorodopsinas/ultraestrutura , Processamento de Imagem Assistida por Computador/métodos , Processamento de Imagem Assistida por Computador/estatística & dados numéricos , Microscopia Eletrônica/estatística & dados numéricos , Modelos MolecularesRESUMO
Automatic prediction of gait events (e.g., heel contact, flat foot, initiation of the swing, etc.) and corresponding profiles of the activations of muscles is important for real-time control of locomotion. This paper presents three supervised machine learning (ML) techniques for prediction of the activation patterns of muscles and sensory data, based on the history of sensory data, for walking assisted by a functional electrical stimulation (FES). Those ML's are: 1) a multilayer perceptron with Levenberg-Marquardt modification of backpropagation learning algorithm; 2) an adaptive-network-based fuzzy inference system (ANFIS); and 3) a combination of an entropy minimization type of inductive learning (IL) technique and a radial basis function (RBF) type of artificial neural network with orthogonal least squares learning algorithm. Here we show the prediction of the activation of the knee flexor muscles and the knee joint angle for seven consecutive strides based on the history of the knee joint angle and the ground reaction forces. The data used for training and testing of ML's was obtained from a simulation of walking assisted with an FES system [39]. The ability of generating rules for an FES controller was selected as the most important criterion when comparing the ML's. Other criteria such as generalization of results, computational complexity, and learning rate were also considered. The minimal number of rules and the most explicit and comprehensible rules were obtained by ANFIS. The best generalization was obtained by the IL and RBF network.
Assuntos
Marcha/fisiologia , Redes Neurais de Computação , Algoritmos , Fenômenos Biomecânicos , Análise por Conglomerados , Entropia , Lógica Fuzzy , Articulação do Quadril/fisiologia , Humanos , Articulação do Joelho/fisiologia , Modelos Biológicos , Músculo Esquelético/fisiologia , Valores de ReferênciaRESUMO
In this paper we describe a novel method for determining synergies between joint motions in reaching movements by hierarchical clustering. A set of recorded elbow and shoulder trajectories is used in a learning algorithm to determine the relationships between angular velocities at elbow and shoulder joints. The learning algorithm is based on optimal criteria for obtaining the hierarchy of descriptions of movement trajectories. We show that this method finds complex synergism between optimal joint trajectories for a given set of data and angular velocities at the shoulder and elbow joints. Three other machine learning techniques (ML) are used for comparison with our method of hierarchical clustering of trajectories. These MLs are: (1) radial basis functions (RBF), (2) inductive learning (IL), and (3) adaptive-network-based fuzzy inference system (ANFIS). Better error characteristics were obtained using the method of hierarchical clustering in comparison with the other techniques. The advantage of the method of hierarchical clustering with respect to the other MLs is in integrating the spatial and temporal elements of reaching movements. Determination and analysis of spatio-temporal events of movement trajectories is a useful tool in designing control systems for functional electrical stimulation (FES) assisted manipulation.
Assuntos
Inteligência Artificial , Articulação do Cotovelo/fisiologia , Movimento , Articulação do Ombro/fisiologia , Algoritmos , Análise por Conglomerados , Lógica Fuzzy , Humanos , Redes Neurais de ComputaçãoRESUMO
Transmission electron microscopy, as most imaging devices, introduces optical aberrations that in the case of thin specimens are usually modeled in Fourier space by the so-called contrast transfer function (CTF). Accurate determination of the CTF is crucial for its posterior correction. Furthermore, the CTF estimation must be fast and robust if high-throughput three-dimensional electron microscopy (3DEM) studies are to be carried out. In this paper we present a robust algorithm that fits a theoretical CTF model to the power spectrum density (PSD) measured on a specific micrograph or micrograph area. Our algorithm is capable of estimating the envelope of the CTF which is absolutely needed for the correction of the CTF amplitude changes.
Assuntos
Astigmatismo/diagnóstico , Processamento de Imagem Assistida por Computador/métodos , Microscopia Eletrônica de Transmissão/métodos , Algoritmos , Desenho de Equipamento , Imageamento Tridimensional , Modelos Estatísticos , Distribuição Normal , Reprodutibilidade dos Testes , Software , Fatores de TempoRESUMO
We combine the self-assembly properties of amphiphilic molecules with the radiolysis method to produce specific sizes and shapes of metallic nano-objects. Radiolysis is used to synthesize core--shell structures consisting of nanometric linoleate spherical micelles as the core and silver as the shell. The validity of the technique is asserted by cryoelectron microscopy, which is an adequate technique for low density contrasts and core--shell structures. The shells are found to be homogeneous with a size of a few nanometers. Images are used to bring forward the hypothesis of the fabrication process.
RESUMO
In a context of automation of cryo-electron microscopy, we developed a novel method for improving visibility of diffraction rings in the power spectra of cryo-electron micrographs of vitreous ice (without carbon film or high concentration of diffracting material). We used these enhanced spectra to semi-automatically detect and remove micrographs and/or local areas introducing errors in the global 3D map (drifted and charged areas) or those unable to increase global signal-to-noise ratio (non-diffracting areas). Our strategy also allows a detection of micrographs/areas with a strong astigmatism. These images should be removed when using algorithms that do not correct astigmatism. Our sorting method is simple and fast since it uses the normalized cross-correlation between enhanced spectra and their copies rotated by 90 degrees. It owes its success mainly to the novel pre-processing of power spectra. The improved visibility also allows an easier visual check of accuracy of sorting. We show that our algorithm can even improve the visibility of diffraction rings of cryo-electron micrographs of pure water. Moreover, we show that this visibility depends strongly on ice thickness. This algorithm is implemented in the Xmipp (open-source image processing package) and is freely available for implementation in any other software package.
Assuntos
Algoritmos , Microscopia Crioeletrônica/métodos , Aumento da Imagem/métodos , Glutamato Sintase/química , Microscopia Eletrônica de Transmissão , Reconhecimento Automatizado de Padrão , Análise Espectral/métodosRESUMO
Transmission electron microscopy is a powerful technique for studying the three-dimensional (3D) structure of a wide range of biological specimens. Knowledge of this structure is crucial for fully understanding complex relationships among macromolecular complexes and organelles in living cells. In this paper, we present the principles and main application domains of 3D transmission electron microscopy in structural biology. Moreover, we survey current developments needed in this field, and discuss the close relationship of 3D transmission electron microscopy with other experimental techniques aimed at obtaining structural and dynamical information from the scale of whole living cells to atomic structure of macromolecular complexes.
Assuntos
Imageamento Tridimensional , Substâncias Macromoleculares , Microscopia Eletrônica de Transmissão/métodos , Algoritmos , Conformação MolecularRESUMO
We have developed an algorithm for the rigid-body registration of a CT volume to a set of C-arm images. The algorithm uses a gradient-based iterative minimization of a least-squares measure of dissimilarity between the C-arm images and projections of the CT volume. To compute projections, we use a novel method for fast integration of the volume along rays. To improve robustness and speed, we take advantage of a coarse-to-fine processing of the volume/image pyramids. To compute the projections of the volume, the gradient of the dissimilarity measure, and the multiresolution data pyramids, we use a continuous image/volume model based on cubic B-splines, which ensures a high interpolation accuracy and a gradient of the dissimilarity measure that is well defined everywhere. We show the performance of our algorithm on a human spine phantom, where the true alignment is determined using a set of fiducial markers.
RESUMO
Measuring the quality of three-dimensional (3D) reconstructed biological macromolecules by transmission electron microscopy is still an open problem. In this article, we extend the applicability of the spectral signal-to-noise ratio (SSNR) to the evaluation of 3D volumes reconstructed with any reconstruction algorithm. The basis of the method is to measure the consistency between the data and a corresponding set of reprojections computed for the reconstructed 3D map. The idiosyncrasies of the reconstruction algorithm are taken explicitly into account by performing a noise-only reconstruction. This results in the definition of a 3D SSNR which provides an objective indicator of the quality of the 3D reconstruction. Furthermore, the information to build the SSNR can be used to produce a volumetric SSNR (VSSNR). Our method overcomes the need to divide the data set in two. It also provides a direct measure of the performance of the reconstruction algorithm itself; this latter information is typically not available with the standard resolution methods which are primarily focused on reproducibility alone.
Assuntos
Processamento de Imagem Assistida por Computador/métodos , Imageamento Tridimensional , Microscopia Eletrônica , Algoritmos , Matemática , Modelos TeóricosRESUMO
Bipedal locomotion was simulated to generate a pattern of activating muscles for walking using electrical stimulation in persons with spinal cord injury (SCI) or stroke. The simulation presented in this study starts from a model of the body determined with user-specific parameters, individualized with respect to the lengths, masses, inertia, muscle and joint properties. The trajectory used for simulation was recorded from an able-bodied subject while walking with ankle-foot orthoses. A discrete mathematical model and dynamic programming were used to determine the optimal control. A cost function was selected as the sum of the squares of the tracking errors from the desired trajectories, and the weighted sum of the squares of agonist and antagonist activations of the muscle groups acting around the hip and knee joints. The aim of the simulation was to study plausible trajectories keeping in mind the limitations imposed by the spinal cord injury or stroke (e.g., spasticity, decreased range of movements in some joints, limited strength of paralyzed, externally activated muscles). If the muscles were capable of generating the movements required and the trajectory was achieved, then the simulation provided two kinds of information: 1) timing of the onset and offset of muscle activations with respect to the various gait events and 2) patterns of activation with respect to the maximum activation. These results are important for synthesizing a rule-based controller.
Assuntos
Transtornos Cerebrovasculares/reabilitação , Simulação por Computador , Modelos Biológicos , Contração Muscular/fisiologia , Paraplegia/reabilitação , Traumatismos da Medula Espinal/reabilitação , Caminhada/fisiologia , Algoritmos , Terapia por Estimulação Elétrica , Eletromiografia , Feminino , Marcha/fisiologia , Articulação do Quadril/fisiologia , Humanos , Articulação do Joelho/fisiologia , Aparelhos OrtopédicosRESUMO
Three-dimensional (3D) electron microscopy (3DEM) aims at the determination of the spatial distribution of the Coulomb potential of macromolecular complexes. The 3D reconstruction of a macromolecule using single-particle techniques involves thousands of 2D projections. One of the key parameters required to perform such a 3D reconstruction is the orientation of each projection image as well as its in-plane orientation. This information is unknown experimentally and must be determined using image-processing techniques. We propose the use of wavelets to match the experimental projections with those obtained from a reference 3D model. The wavelet decomposition of the projection images provides a framework for a multiscale matching algorithm in which speed and robustness to noise are gained. Furthermore, this multiresolution approach is combined with a novel orientation selection strategy. Results obtained from computer simulations as well as experimental data encourage the use of this approach.